(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公表特許公報(A)
(11)【公表番号】
(43)【公表日】2024-08-21
(54)【発明の名称】多重免疫蛍光イメージングを使用する組織学的染色のデジタル合成
(51)【国際特許分類】
G06T 5/60 20240101AFI20240814BHJP
G06T 1/00 20060101ALI20240814BHJP
H04N 1/46 20060101ALI20240814BHJP
G01N 33/48 20060101ALN20240814BHJP
G01N 33/483 20060101ALN20240814BHJP
G01N 21/64 20060101ALN20240814BHJP
【FI】
G06T5/60
G06T1/00 510
H04N1/46
G01N33/48 M
G01N33/483 C
G01N21/64 F
G01N33/48 P
【審査請求】有
【予備審査請求】未請求
(21)【出願番号】P 2023577709
(86)(22)【出願日】2022-06-07
(85)【翻訳文提出日】2023-12-15
(86)【国際出願番号】 US2022032468
(87)【国際公開番号】W WO2022265883
(87)【国際公開日】2022-12-22
(32)【優先日】2021-06-17
(33)【優先権主張国・地域又は機関】US
(81)【指定国・地域】
(71)【出願人】
【識別番号】507179346
【氏名又は名称】ベンタナ メディカル システムズ, インコーポレイテッド
(74)【代理人】
【識別番号】100118902
【氏名又は名称】山本 修
(74)【代理人】
【識別番号】100106208
【氏名又は名称】宮前 徹
(74)【代理人】
【識別番号】100196508
【氏名又は名称】松尾 淳一
(74)【代理人】
【識別番号】100138759
【氏名又は名称】大房 直樹
(72)【発明者】
【氏名】ワーン,シンウェイ
(72)【発明者】
【氏名】ヂャオ,ズオ
(72)【発明者】
【氏名】ロルサクル,アウラヌチ
(72)【発明者】
【氏名】ニエ,ヤオ
【テーマコード(参考)】
2G043
2G045
5B057
【Fターム(参考)】
2G043AA01
2G043AA03
2G043BA16
2G043EA01
2G043FA02
2G043KA02
2G043LA03
2G045AA24
2G045CB01
2G045FB03
2G045FB07
2G045FB12
2G045GC15
2G045JA01
5B057AA10
5B057BA02
5B057CA01
5B057CA08
5B057CA12
5B057CA16
5B057CB01
5B057CB08
5B057CB12
5B057CB16
5B057CC01
5B057CE08
5B057CE16
(57)【要約】
多重免疫蛍光(MPX)画像から合成組織化学的染色画像を取得するための技術は、組織切片のMPX画像のMチャンネルのそれぞれからの情報に基づくNチャンネル入力画像を生成することであって、MおよびNが正の整数であり、NがM以下である、Nチャンネル入力画像を生成することと、複数の画像対を含む訓練データセットを使用して訓練された生成器ネットワークを使用してNチャンネル入力画像を処理することによって合成画像を生成することと、を含み得る。合成画像は、少なくとも1つの組織化学的染色剤によって染色された組織切片を示す。複数の画像対の各画像対は、組織の第1の切片のMPX画像から生成されたNチャンネル画像と、少なくとも1つの組織化学的染色剤によって染色された組織の第2の切片の画像とを含む。
【選択図】
図8A
【特許請求の範囲】
【請求項1】
画像変換の方法であって、
組織切片の多重免疫蛍光(MPX)画像のMチャンネルのそれぞれからの情報に基づくNチャンネル入力画像を生成することであって、Mは正の整数であり、NはM以下の正の整数である、Nチャンネル入力画像を生成することと、
複数の画像対を含む訓練データセットを使用して訓練された生成器ネットワークを使用して前記Nチャンネル入力画像を処理することによって合成画像を生成することと
を含み、
前記合成画像が、少なくとも1つの組織化学的染色剤によって染色された組織切片を示し、
前記複数の画像対の各画像対について、前記対が、
組織の第1の切片のMPX画像から生成されたNチャンネル画像と、
前記少なくとも1つの組織化学的染色剤によって染色された前記組織の第2の切片の画像と
を含む、方法。
【請求項2】
前記Nチャンネル入力画像を生成することが、前記組織切片の前記MPX画像の前記Mチャンネルのそれぞれについて、前記チャンネルからの情報を前記Nチャンネル画像の前記Nチャンネルのうちの少なくとも1つにマッピングすることを含む、請求項1に記載の方法。
【請求項3】
前記マッピングすることが、前記組織切片の前記MPX画像の複数のチャンネルのそれぞれからの情報に基づく自己蛍光画像を生成することを含み、
前記Nチャンネル入力画像が、前記自己蛍光画像からの情報に基づく、請求項2に記載の方法。
【請求項4】
前記自己蛍光画像が、前記組織切片の前記MPX画像の前記複数のチャンネルの非線形結合に基づく、請求項3に記載の方法。
【請求項5】
前記自己蛍光画像が、前記複数のチャンネル間の自己蛍光のスペクトル分布に基づく、請求項3または4に記載の方法。
【請求項6】
前記MPX画像の前記複数のチャンネルが、前記組織切片の前記MPX画像の前記Mチャンネルである、請求項3から5のいずれか一項に記載の方法。
【請求項7】
前記複数の画像対の各画像対について、前記対の前記Nチャンネル画像が、前記第1の切片の前記MPX画像のMチャンネルのそれぞれからの情報に基づく、請求項1から6のいずれか一項に記載の方法。
【請求項8】
前記合成画像がNチャンネル画像である、請求項1から7のいずれか一項に記載の方法。
【請求項9】
前記Nチャンネル入力画像および前記合成画像のそれぞれがRGB画像である、請求項8に記載の方法。
【請求項10】
前記組織切片の前記MPX画像が暗視野画像であり、前記合成画像が明視野画像である、請求項1から9のいずれか一項に記載の方法。
【請求項11】
前記少なくとも1つの組織化学的染色剤がヘマトキシリンおよびエオシンである、請求項1から10のいずれか一項に記載の方法。
【請求項12】
前記Nチャンネル画像が、前記自己蛍光画像と、核対比染色に対応する前記組織切片の前記MPX画像のチャンネルとに基づく、請求項1から11のいずれか一項に記載の方法。
【請求項13】
前記Nチャンネル画像が、前記自己蛍光画像と、核対比染色に対応する前記組織切片の前記MPX画像のチャンネルとの線形結合に基づく、請求項1から11のいずれか一項に記載の方法。
【請求項14】
前記Nチャンネル画像が、前記線形結合に基づく光学濃度値の配列に基づく、請求項1から13のいずれか一項に記載の方法。
【請求項15】
Nが3に等しい、請求項1から14のいずれか一項に記載の方法。
【請求項16】
Mが少なくとも4である、請求項1から15のいずれか一項に記載の方法。
【請求項17】
前記生成器ネットワークが、敵対的生成ネットワークの一部として訓練された、請求項1から16のいずれか一項に記載の方法。
【請求項18】
前記生成器ネットワークが、U-Netとして実装される、請求項1から17のいずれか一項に記載の方法。
【請求項19】
前記生成器ネットワークが、エンコーダ-デコーダネットワークとして実装される、請求項1から18のいずれか一項に記載の方法。
【請求項20】
前記生成器ネットワークが、前記生成器ネットワークによる画像と予想される出力画像との間で測定されたL1損失を介して更新される、請求項1から19のいずれか一項に記載の方法。
【請求項21】
ユーザによって、対象の診断を決定する方法であって、前記決定することが、請求項1から20のいずれか一項に記載の前記合成画像に基づく、方法。
【請求項22】
前記ユーザによって、(i)前記合成画像、および/または(ii)前記対象の前記診断に基づく化合物による処置を施すことをさらに含む、請求項21に記載の方法。
【請求項23】
1つまたは複数のデータプロセッサと、
命令を含む非一時的コンピュータ可読記憶媒体であって、前記命令が、前記1つまたは複数のデータプロセッサ上で実行されると、前記1つまたは複数のデータプロセッサに、請求項1~20のいずれか一項に記載の画像変換の方法を実行させる、非一時的コンピュータ可読記憶媒体と、を備える、システム。
【請求項24】
1つまたは複数のデータプロセッサに請求項1~20のいずれか一項に記載の画像変換の方法を実行させるように構成された命令を含む、非一時的機械可読記憶媒体に有形に具現化されたコンピュータプログラム製品。
【発明の詳細な説明】
【技術分野】
【0001】
分野
本開示は、デジタル病理に関し、特に、多重免疫蛍光画像から合成組織学的染色画像を取得することを含む技術に関する。
【背景技術】
【0002】
背景
組織病理は、疾患の診断、治療に対する応答の評価、および/または疾患と戦うための薬理学的薬剤の開発などの様々な理由で、組織切片から準備されたスライドの検査を含み得る。組織切片およびその中の細胞は実質的に透明であることから、スライドの準備は、典型的には、関連する構造をより視認可能にするために組織切片を染色することを含む。デジタル病理は、デジタル画像を取得するために染色されたスライドをスキャンすることを含み得、デジタル画像は、その後、デジタル病理画像分析によって検査され得、および/またはヒト病理学者によって解釈され得る。
【0003】
組織切片の多重免疫蛍光(MPX)染色は、単一細胞レベルでの複数のバイオマーカーおよびそれらの同時発現の同時検出を可能にする。MPXは、免疫療法に対する応答に有意な影響を及ぼし得る腫瘍微小環境における免疫状況の特徴付けを可能にする。例えば、MPXの使用は、単一のスライドにおいてより多くのバイオマーカーの検出および共局在化を提供することができる。MPXは、効果的な治療法を発見し、新たな薬物を開発するための貴重なツールであることを約束するが、MPXの結果を基礎となる組織構造と正確に関連付けることに関していくつかの課題が残っている。
【発明の概要】
【0004】
概要
様々な実施形態では、組織切片の多重免疫蛍光(MPX)画像のM個のチャンネルのそれぞれからの情報に基づくNチャンネル入力画像を生成することであって、Mが正の整数であり、NがM以下の正の整数である、Nチャンネル入力画像を生成することと、複数の画像対を含む訓練データセットを使用して訓練された生成器ネットワークを使用してNチャンネル入力画像を処理することによって合成画像を生成することと、を含む、画像変換のコンピュータ実装方法が提供される。この方法では、合成画像は、少なくとも1つの組織化学的染色剤によって染色された組織切片を示し、複数対の画像の各対について、対は、組織の第1の切片のMPX画像から生成されたNチャンネル画像と、少なくとも1つの組織化学的染色剤によって染色された組織の第2の切片の画像とを含む。
【0005】
いくつかの実施形態では、Nチャンネル入力画像を生成することは、組織切片のMPX画像のM個のチャンネルのそれぞれについて、チャンネルからNチャンネル画像のN個のチャンネルのうちの少なくとも1つへのマッピング情報を含む。マッピングは、組織切片のMPX画像の複数のチャンネル(例えば、M個のチャンネルのそれぞれ)のそれぞれからの情報に基づく自己蛍光画像を作成することを含んでもよく、Nチャンネル入力画像は、自己蛍光画像からの情報に基づいてもよい。自己蛍光画像は、組織切片のMPX画像の複数のチャンネルの非線形な組み合わせに基づいてもよい。追加的または代替的に、自己蛍光画像は、複数のチャンネル間の自己蛍光のスペクトル分布に基づいてもよい。
【0006】
いくつかの実施形態では、合成画像は、Nチャンネル画像である。Nチャンネル入力画像および合成画像のそれぞれは、RGB画像であってもよい。追加的または代替的に、組織切片のMPX画像は、暗視野画像であってもよく、および/または合成画像は、明視野画像であってもよい。
【0007】
いくつかの実施形態では、少なくとも1つの組織化学的染色剤は、ヘマトキシリンおよびエオシンである。
【0008】
いくつかの実施形態では、Nチャンネル画像は、自己蛍光画像および核対比染色画像に基づく。例えば、Nチャンネル画像は、自己蛍光画像と核対比染色に対応する組織切片のMPX画像のチャンネルとの線形結合に基づいてもよい。Nチャンネル画像は、線形結合に基づく光学濃度値の配列に基づく。
【0009】
いくつかの実施形態では、Nは3に等しく、および/またはMは少なくとも4である。
【0010】
いくつかの実施形態では、生成器ネットワークは、敵対的生成ネットワークの一部として訓練されたものである。いくつかの実施形態では、生成器ネットワークは、U-Netおよび/またはエンコーダ-デコーダネットワークとして実装される。生成器ネットワークは、生成器ネットワークによる画像と予想される出力画像との間で測定されたL1損失を介して更新され得る。
【0011】
いくつかの実施形態では、コンピュータ実装方法の実施形態によって生成された合成画像に基づいて、ユーザによって対象の診断を決定することを含む方法が提供される。
【0012】
いくつかの実施形態では、ユーザによって、コンピュータ実装方法の実施形態によって生成された合成画像に基づいて対象の診断を決定することと、ユーザによって、(i)合成画像、および/または(ii)対象の診断に基づいて化合物による処置を施すことと、を含む、方法が提供される。
【0013】
いくつかの実施形態では、1つまたは複数のデータプロセッサと、1つまたは複数のデータプロセッサ上で実行されると、1つまたは複数のデータプロセッサに、本明細書に開示された1つまたは複数の方法の一部または全部を実行させる命令を含む非一時的コンピュータ可読記憶媒体とを含むシステムが提供される。
【0014】
いくつかの実施形態では、非一時的機械可読記憶媒体に有形に具現化され、1つまたは複数のデータプロセッサに、本明細書に開示された1つまたは複数の方法の一部または全部を実行させるように構成された命令を含む、コンピュータプログラム製品が提供される。
【0015】
本開示のいくつかの実施形態は、1つまたは複数のデータプロセッサを含むシステムを含む。いくつかの実施形態では、システムは、命令を含む非一時的コンピュータ可読記憶媒体であって、該命令が、1つまたは複数のデータプロセッサ上で実行されると、1つまたは複数のデータプロセッサに、本明細書に開示された1つまたは複数の方法の一部または全部および/または1つまたは複数のプロセスの一部または全部を実行させる、非一時的コンピュータ可読記憶媒体を含む。本開示のいくつかの実施形態は、1つまたは複数のデータプロセッサに、本明細書に開示された1つまたは複数の方法の一部または全部および/または1つまたは複数のプロセスの一部または全部を実行させるように構成された命令を含む非一時的機械可読記憶媒体において有形に具現化された、コンピュータプログラム製品を含む。
【0016】
使用された用語および表現は、限定ではなく説明の用語として使用され、そのような用語および表現の使用において、示されて説明された特徴のいかなる均等物またはその一部も除外する意図はないが、特許請求の範囲に記載された発明の範囲内で様々な変更が可能であることが認識される。したがって、特許請求の範囲に記載された本発明は、実施形態および任意の特徴によって具体的に開示されているが、本明細書に開示された概念の変更および変形は、当業者によってあてにされてもよく、そのような変更および変形は、添付の特許請求の範囲によって定義される本発明の範囲内にあると見なされることを理解されたい。
【図面の簡単な説明】
【0017】
本特許または出願ファイルは、カラーで作成された少なくとも1つの図面を含んでいる。カラー図面を含む本特許または特許出願公開の写しは、請求に応じて、必要な手数料を支払うことにより、特許庁によって提供されることになる。
【0018】
様々な実施形態の態様および特徴は、添付の図面を参照して例を説明することによってより明らかになるであろう。
【0019】
【
図1】デジタル病理学的ソリューションのワークフローの例示的な図を示している。
【
図2】多重免疫蛍光(MPX)イメージングワークフローの例示的な図および例示的なMPX画像を示している。
【
図3】MPXパネルからのスライドの例示的なMPX画像を示している。
【
図4】MPXパネルからのスライドの別の例示的なMPX画像を示している。
【
図5】H&E(明視野)画像(
図5B)に位置合わせされたMPX(暗視野)画像(
図5A)を示している。
【
図6】H&E染色された試料の画像を取得するためのワークフローの例を示している。
【
図7A】H&E画像および実際のH&E画像に類似するマッピング画像の例を示している。
【
図7B】H&E画像および実際のH&E画像に類似するマッピング画像の例を示している。
【
図8A】いくつかの実施形態にかかる例示的なプロセスのフローチャートを示している。
【
図8B】いくつかの実施形態にかかる別の例示的なプロセスのフローチャートを示している。
【
図9】いくつかの実施形態にかかる自己蛍光画像を生成するプロセスを示している。
【
図10A】いくつかの実施形態にかかる、MPX画像のM個のチャンネルの異なる対応するマッピングからそれぞれが取得された、3チャンネル(RGB)マッピングされた入力画像の例を示している。
【
図10B】いくつかの実施形態にかかる、MPX画像のM個のチャンネルの異なる対応するマッピングからそれぞれが取得された、3チャンネル(RGB)マッピングされた入力画像の例を示している。
【
図10C】いくつかの実施形態にかかる、MPX画像のM個のチャンネルの異なる対応するマッピングからそれぞれが取得された、3チャンネル(RGB)マッピングされた入力画像の例を示している。
【
図10D】いくつかの実施形態にかかる、MPX画像のM個のチャンネルの異なる対応するマッピングからそれぞれが取得された、3チャンネル(RGB)マッピングされた入力画像の例を示している。
【
図11】いくつかの実施形態にかかる、(例えば、深層学習ネットワークを訓練するため、および/または推論のために)3チャンネルマッピングされた入力画像を生成するプロセスを示している。
【
図12】いくつかの実施形態にかかる、3チャンネルマッピングされた入力画像を生成する別のプロセスを示している。
【
図13】いくつかの実施形態にかかる、3チャンネルマッピングされた入力画像を生成するさらなるプロセスを示している。
【
図14】いくつかの実施形態にかかる例示的なコンピューティング環境を示している。
【
図15】いくつかの実施形態にかかる条件付きGANモデルを示している。
【
図16】いくつかの実施形態にかかるPix2Pix GANを使用する条件付きGANモデルの実装の例を示している。
【
図17】いくつかの実施形態にかかる訓練画像を取得するプロセスについてのフローチャートを示している。
【
図18】いくつかの実施形態にかかる、全スライド画像からのタイルおよび画像パッチの抽出の例を示している。
【
図19】いくつかの実施形態にかかる、サイクルGANにおけるネットワークおよびネットワーク接続の表現を示している。
【
図20】いくつかの実施形態にかかる、サイクルGANを使用して画像を生成および識別するフローを示している。
【
図21】様々な実施形態にしたがって生成された入力画像、標的画像、および出力画像の例を示している。
【
図22】様々な実施形態にしたがって生成された入力画像、標的画像、および出力画像の例を示している。
【
図23】様々な実施形態にしたがって生成された入力画像、標的画像、および出力画像の例を示している。
【
図24】異なる4種類の組織における自己蛍光の基準分光分布の比較例を示している。
【
図25】無染色の組織試料のMPX画像の2つの例を示している。
【
図26】補正前後の膵臓組織試料のMPX画像の例を示している。
【発明を実施するための形態】
【0020】
詳細な説明
本明細書に開示されるシステム、方法およびソフトウェアは、多重免疫蛍光画像から合成組織学的染色画像を取得することを容易にする。特定の実施形態が説明されているが、これらの実施形態は例としてのみ提示されており、保護の範囲を限定することを意図するものではない。本明細書に記載された装置、方法、およびシステムは、様々な他の形態で具現化され得る。さらにまた、保護の範囲から逸脱することなく、本明細書に記載の例示的な方法およびシステムの形態の様々な省略、置換、および変更が行われてもよい。
【0021】
I.概要
デジタル病理は、対象を正確に診断し、治療上の意思決定を導くために、デジタル化された画像の解釈を伴い得る。デジタル病理学的ソリューションでは、関心のある生物学的オブジェクト(例えば、陽性/陰性腫瘍細胞など)を自動的に検出または分類するために、画像分析ワークフローが確立され得る。
図1は、デジタル病理学的ソリューションのワークフロー100の例示的な図を示している。デジタル病理学的ソリューションワークフロー100は、ブロック110において組織スライドを取得することと、ブロック120においてデジタル画像を取得するためにデジタル画像スキャナ(例えば、全スライド画像(WSI)スキャナ)によって予め選択された領域または組織スライドの全体をスキャンすることと、ブロック130において1つまたは複数の画像分析アルゴリズムを使用してデジタル画像に対して画像分析を実行することと、ブロック140において画像分析に基づいて関心のあるオブジェクトをスコアリング(例えば、陽性、陰性、中程度、弱いなどの定量的スコアリングまたは半定量的スコアリング)することと、を含む。
【0022】
例えば疾患によって引き起こされる組織変化の評価は、薄い組織切片を検査することによって行われ得る。一連の切片を取得するために組織試料(例えば、腫瘍の試料)がスライスされ得、各切片は、例えば、4~5ミクロンの厚さを有する。組織切片およびその中の細胞は実質的に透明であることから、スライドの準備は、典型的には、関連する構造をより視認可能にするために組織切片を染色することを含む。例えば、組織の異なる部分を1つまたは複数の異なる染色剤によって染色して、組織の異なる特徴を表し得る。
【0023】
各切片は、スライドに取り付けられてもよく、スライドは、その後、デジタル画像を作成するためにスキャンされ、デジタル画像は、その後、デジタル病理画像分析によって検査されおよび/または(例えば、画像ビューアソフトウェアを使用する)ヒト病理学者によって解釈され得る。病理学者は、(例えば、関心のある生物学的物体を検出および分類するために)有意な定量的測定値を抽出するために画像分析アルゴリズムの使用を可能にするために、スライドのデジタル画像をレビューし、手動で注釈を付け得る(例えば、腫瘍領域、壊死など)。従来、病理学者は、各連続する組織切片上の同じ態様を識別するために、組織試料からの複数の組織切片の各連続する画像に手動で注釈を付け得る。
【0024】
組織染色の1つのタイプは、組織構造を染色するために1つまたは複数の化学染料(例えば、酸性染料、塩基性染料)を使用する組織化学的染色である。組織化学的染色は、組織形態および/または細胞の微小解剖学的構造(例えば、細胞核を細胞質と区別するため、脂質滴を示すためなど)の一般的な態様を示すために使用され得る。組織化学的染色の一例は、ヘマトキシリンおよびエオシン(H&E)である。組織化学的染色の他の例は、3色染色(例えば、マッソントリクローム)、周期的酸シッフ(PAS)、銀染色および鉄染色を含む。組織化学的に染色された試料の画像は、典型的には、明視野顕微鏡法を使用して取得され、組織化学的に染色された試料のスライドは、典型的には、カラーチャンネルあたり画素あたり8ビット(例えば、赤、緑、青(RGB)のそれぞれ)の解像度でスキャンされる。
【0025】
組織試料のH&E染色では、ヘマトキシリンは、細胞核を青色に染色し、エオシンは、細胞外マトリックスおよび細胞質を桃色に染色し、他の構造は、異なる色調、色相、および/または桃色と青色の組み合わせを有するように染色され得る。しかしながら、H&E染色は、一般的な組織および細胞の解剖学的構造を識別するために有用であるが、下記に記載されるような免疫組織学的化学によって提供され得る、様々なタイプの癌を区別するために使用され得る情報(例えば、HER2スコアリング)などの、特定の診断評価を裏付けるために必要な特定の情報を提供することができない。
【0026】
組織染色の別のタイプは、目的の標的抗原(バイオマーカーとも呼ばれる)に特異的に結合する一次抗体を使用する免疫組織学的化学(IHC、「免疫染色」とも呼ばれる)である。IHCは、直接的または間接的であり得る。直接IHCでは、一次抗体は、標識(例えば、発色団またはフルオロフォア)に直接コンジュゲートされる。間接IHCでは、一次抗体が最初に標的抗原に結合し、次いで標識(例えば、発色団またはフルオロフォア)と結合した二次抗体が一次抗体に結合される。組織染色のためのIHCの使用は、典型的には、非常に高価な試薬の使用ならびに組織化学的染色よりも複雑な検査室機器および手順を必要とする。IHCフルオロフォアによって染色された試料の画像は、通常、暗視野顕微鏡を使用して取得される。
【0027】
多重免疫蛍光(MPX)イメージングでは、単一のスライドが複数のIHCフルオロフォアによって染色され、様々なフルオロフォアの発光スペクトルは互いに区別可能である。MPXパネル内の異なるフルオロフォアの数は、典型的には5つまたは6つであるが、それより少なくても多くてもよい。所望のバイオマーカーを標識するフルオロフォアに加えて、蛍光発生核対比染色剤、例えば4’,6-ジアミジノ-2-フェニルインドール(DAPI)が典型的に含まれる。MPX画像は、各フルオロフォアから暗視野画像の別個の対応するチャンネルへの発光(典型的には、チャンネルあたり画素あたり16ビットの解像度)を取り込むことによって取得され、MPX画像は、各チャンネルに異なる擬似カラーを割り当てることによって観察される。
【0028】
図2、
図3、および
図4は、IHC染色の3つの異なるパネルのMPX画像を示している。
図2の画像では、バイオマーカーCD8A、プログラム細胞死タンパク質1(PDCD1)、分化クラスター3-イプシロン(CD3E)、増殖マーカーKi-67(MKI67)、およびケラチン(KRT)に対応するチャンネルには、それぞれ偽色黄色、淡青色、紫色、赤色、および緑色が割り当てられている。
図3の画像では、核対比染色DAPIならびにバイオマーカーPanCK、CD31、線維芽細胞活性化タンパク質アルファ(FAP)および主要組織適合遺伝子複合体1(MHC1)に対応するチャンネルには、それぞれ偽色、暗青、薄青、緑、赤および桃色が割り当てられている。
図4の画像では、DAPIおよびバイオマーカーCD3(明るい青色)、プログラム死リガンド1(PD-L1)、PD1、PanCK、およびグランザイムB(Gzmb)に対応するチャンネルには、それぞれ偽色、暗い青色、明るい青色、緑色、黄色、赤色、および桃色が割り当てられている。
【0029】
IHCと比較して、(例えば、
図2~
図4に示すように)単一のMPXスライドにおける複数のバイオマーカーの検出は、バイオマーカーの共配置および共発現に関する貴重な情報を提供する。例えば、5つのバイオマーカーを同時に検出することは、単一のMPXスライドにおける多数(例えば、17から90)の表現型(定義されたバイオマーカーのセットの同時発現)の検出を可能にする。
【0030】
図2はまた、MPXイメージングワークフローの例を示している。そのようなワークフローは、アルゴリズム開発(例えば、クロストーク補正、自己蛍光補正、候補選択、病理学者のグラウンドトゥルース、特徴に関する機械学習、および/またはPD1強度クラスを含む)から始まり、アルゴリズム検証(例えば、アルゴリズムからの結果を病理学者の結果と比較する)、臨床報告および分析、ならびにアッセイ開発、続いてサンドイッチアッセイ検証を通して進行し得る。
【0031】
組織化学的に染色された(例えば、H&E)画像とは異なり、MPX画像は、典型的には、下にある組織構造の明確な描写を欠いている。その結果、さらなる分析のために関心領域(ROI)(例えば、腫瘍またはエピトープ領域)および/またはさらなる分析から除外される領域(例えば、壊死領域)を描写するためにMPX画像に注釈を付けるタスクは、単独では困難または不可能である可能性がある。そのような理由から、MPX画像に注釈を付ける、またはMPX画像に基づいて診断を確認する現在の手法は、一般に、必要な構造的詳細を提供するためにH&E染色された試料の近くの部分の対応する画像に依存する。H&E画像は、病理学者が(例えば、注釈付けを実行している間に)ディスプレイ上で2つの画像を容易に切り替えることができるように、MPX画像に位置合わせされる。
図6は、H&E画像を生成するためのワークフローの例を示しており、これは、固定、加工、包埋、切断、染色、および顕微鏡検査などの操作を含み得る。
【0032】
図5Aは、注釈付きMPX画像の例を示しており、
図5Bは、H&E染色され、MPX画像に位置合わせされた試料の近傍部分の対応画像の例を示している。基準H&E画像の作成には追加のコストがかかり、2つの画像のレビューには病理学者が追加の組織注釈をレビューする必要があり得る。一方の画像を他方に位置合わせすることは、スライド準備の手順中に各切片の固有の全体的および/または局所的な変形が生じる可能性があるため、時間がかかり、計算コストがかかる場合がある。位置合わせが不正確である場合、位置合わせされた画像が重ね合わされ、切り替えられ、または他の方法で一緒に見られるときに示されるバイオマーカーと下にある組織構造(例えば、細胞膜、核)との間の空間的関係も不正確であり得る。さらに、画像が最適に位置合わせされたとしても、2つの画像は、異なる部分からのものであるため、同じ組織空間を描写しない。一方の切片の構造は、他方の切片に存在しなくてもよく(例えば、そのサイズが各部分の厚さに匹敵する構造)、または他方の切片に異なるサイズおよび/または形態を有してもよく、その結果、ベースライン画像との位置合わせ画像のオーバーレイは、組織構造に対するバイオマーカーの位置を正確に取り込まないことがあり、これは、MPX画像の注釈付けおよび/またはレビュー中に困難をもたらし得る。
【0033】
MPX画像の1つまたは複数のチャンネルをRGB色空間にマッピングすることによって、H&E染色試料に似た画像を作成するために、従来の画像処理技術が使用され得る。
図7Aは、そのようなマッピング画像の例を示しており、
図7Bは、比較のための真のH&E染色画像を示している。
図7Aに示されるようなマッピングされた画像は、H&E画像に表面的に類似し得るが、日常的な臨床診療における病理学者は、そのような画像が、例えば、診断における実際の決定を支援するには不十分な画像の現実および品質を示すことを確認している。例えば、
図7Aおよび
図7B中の画像の右端に近い領域を比較することによって、
図7Aのマッピングされた画像は、
図7Bの中間レベルの値で明らかな構造的詳細の多くを欠いていることが分かる。
【0034】
これらの制限および他の制限を克服するために、MPX画像から、組織化学的に染色された試料の実際の画像との高度の構造類似性を有する組織化学的に染色された(例えば、H&E染色)試料を描写する合成(仮想的に染色された)画像を生成するための技術が本明細書に開示される。そのような合成画像の生成は、複雑で面倒な組織の調製および染色プロセスならびに位置合わせ操作の回避を支援し得る。臨床的には、そのような解決策は、MPX画像をオンおよびオフに切り替えることによってH&E染色組織構造を視覚化することを可能にするという利点を有し得、これは、病理学者の標準的な実践を満たし、多重イメージングの臨床的採用を促進し得る。開発的には、そのような解決策は、アルゴリズム開発のためのH&E分析のためのグラウンドトゥルースを生成する可能性を有し、バイオマーカーの発見を支援することができる別のMPX画像へのMPX画像のより容易な位置合わせを可能にし得る。
【0035】
合成画像の生成は、敵対的生成ネットワーク(GAN)を訓練する間に学習されたパラメータを含み得る訓練された生成器ネットワークによって実行され得る。GANは、入力画像が偽(すなわち、生成器ネットワークによって生成されている)であるか、または現実(すなわち、対象から収集された実際の画像を示す)であるかを予測するように構成された識別器ネットワークをさらに含み得る。これらの予測の精度に基づくフィードバックは、訓練中に生成器ネットワークに提供され得る。
【0036】
本開示の1つの例示的な実施形態は、組織部分のMPX画像のM個のチャンネルのそれぞれからの情報に基づくNチャンネル入力画像を生成することであって、Mが正の整数であり、NがM以下の正の整数である、Nチャンネル入力画像を生成することと、複数の画像対を含む訓練データセットを使用して訓練された生成器ネットワークを使用してNチャンネル入力画像を処理することによって合成画像を生成することと、を含む、画像変換の方法に関する。この方法では、合成画像は、少なくとも1つの組織化学的染色剤によって染色された組織切片を示し、複数対の画像の各対について、対は、組織の第1の切片のMPX画像から生成されたNチャンネル画像と、少なくとも1つの組織化学的染色剤によって染色された組織の第2の切片の画像とを含む。
【0037】
II.定義
本明細書で使用される場合、動作が何かに「基づく」場合、これは、動作が何かの少なくとも一部に少なくとも部分的に基づくことを意味する。
【0038】
本明細書で使用される場合、「実質的に(substantially)」、「およそ(approximately)」、および「約(about)」という用語は、当業者によって理解されるように、大部分が指定されるものであるが、必ずしも完全には指定されないもの(および完全に指定されるものを含む)として定義される。任意の開示された実施形態では、「実質的に」、「およそ」、または「約」という用語は、指定されたものの「[パーセンテージ]以内」で置き換えられてもよく、パーセンテージは、0.1、1、5、および10%を含む。
【0039】
本明細書で使用される場合、「試料」、「生物学的試料」または「組織試料」という用語は、ウイルスを含む任意の生物から取得される生物学的分子(タンパク質、ペプチド、核酸、脂質、炭水化物、またはそれらの組み合わせなど)を含む任意の試料を指す。生物の他の例は、哺乳類(ヒト、ネコ、イヌ、ウマ、ウシ、およびブタなどの獣医動物、ならびにマウス、ラット、霊長類などの実験動物など)、昆虫、環形動物、クモ形類動物、有袋類、爬虫類、両生類、細菌、および菌類などを含む。生物学的試料は、組織試料(組織切片や組織の針生検など)、細胞試料(Pap塗抹標本もしくは血液塗抹標本などの細胞学的塗抹標本、またはマイクロダイセクションによって取得された細胞の試料など)、または細胞分画、断片または細胞小器官(細胞を溶解し、遠心分離などによってそれらの成分を分離することによって取得される)を含む。生物学的試料の他の例は、血液、血清、尿、精液、糞便、脳脊髄液、間質液、粘膜、涙、汗、膿、生検組織(例えば、外科的生検または針生検によって取得される)、乳頭吸引物、耳垢、乳、膣液、唾液、ぬぐい液(頬スワブなど)、または最初の生物学的試料に由来する生物学的分子を含む任意の材料を含む。特定の実施形態では、本明細書で使用される「生物学的試料」という用語は、対象から得られた腫瘍またはその一部から調製された試料(均質化または液化された試料など)を指す。
【0040】
III.MPXイメージングを使用する組織学的染色のデジタル合成のための技術
本明細書に記載のMPXイメージングを使用する組織学的染色のデジタル合成の手法は、MPXの暗視野イメージングからH&E染色の明視野イメージングに変換するためにマルチモダリティ転写の技術を使用し得る。そのような手法は、異なる画像チャンネルからの情報を使用して、深層学習ネットワーク(例えば、敵対的生成ネットワークまたはGAN)に入力するための中間および前処理された訓練データを取得する。そのような手法は、入力訓練データセットを生成するためにMPX画像のチャンネルの多くまたは全てからの情報を使用することを含み得る。
【0041】
図8Aは、ソース画像(例えば、処理されるソース画像のセットからのソース画像)を標的画像と同様の特性を有する新たな画像(例えば、生成される新たな画像のセットの新たな画像)に変換する例示的なプロセス800のフローチャートを示している。プロセス800は、(例えば、
図14~
図16、
図19、および
図20に関して本明細書で説明するように)1つまたは複数のコンピューティングシステム、モデル、およびネットワークを使用して実行され得る。
【0042】
図8Aを参照すると、ブロック804において、組織部分のMPX画像のM個のチャンネルのそれぞれからの情報に基づくNチャンネル入力画像が生成され、Mは正の整数であり、NはM以下の正の整数である。ブロック808において、生成器ネットワークを使用してNチャンネル入力画像を処理することによって合成画像が生成される。合成画像は、少なくとも1つの組織化学的染色剤によって染色されている組織切片を示す。生成器ネットワークは、複数の画像対を含む訓練データセットを使用して訓練されており、各対は、(例えば、第1の切片のMPX画像のM個のチャンネルのそれぞれからの情報に基づいて)組織の第1の切片のMPX画像から生成されたNチャンネル画像と、少なくとも1つの組織化学的染色剤によって染色された組織の第2の切片の画像とを含む。
【0043】
ブロック804において生成された入力画像のチャンネル数Nは、ブロック808において生成された合成画像のチャンネル数に等しいことが望ましい場合がある。例えば、入力画像と同じ数のチャンネルを有する出力画像を生成するように構成されたPix2Pix GANなどの生成器ネットワークを使用してブロック808を実装することが望ましい場合がある。合成画像は、赤色チャンネル(R)、緑色チャンネル(G)、および青色チャンネル(B)を有するRGB画像などの3チャンネル画像であってもよい。
【0044】
図8Bは、ブロック816を含むプロセスブロック804の実装812を示している。ブロック816において、組織切片のMPX画像のM個のチャンネルのそれぞれについて、チャンネルからの情報が、Nチャンネル画像のNチャンネルのうちの少なくとも1つにマッピングされる。そのようなマッピングは、線形および/または非線形マッピングを含み得る。この例では、ブロック816は、ブロック820を含む。ブロック820において、組織切片のMPX画像の複数のチャンネルのそれぞれから(例えば、M個のチャンネルのそれぞれから)の情報に基づく自己蛍光画像が生成される。プロセスブロック804のこの実装では、Nチャンネル入力画像は、自己蛍光画像からの情報に基づく。ブロック812のインスタンスはまた、MPX画像のソースセットのうちの対応するものから訓練データセットの入力画像のそれぞれを生成するために実行され得る。
【0045】
プロセス800のいくつかの実施形態では、組織化学的染色剤は、ヘマトキシリンおよびエオシンである。
【0046】
プロセス800のいくつかの実施形態では、生成器ネットワークは、敵対的生成ネットワーク(例えば、cGAN、Pix2Pix GAN、またはサイクルGAN)の一部として訓練されたものである。
【0047】
自己蛍光画像は、(例えば、ブロック820において)MPX画像の複数のチャンネルの非線形な組み合わせとして生成され得る。一例では、自己蛍光画像の画素ごとに、MPX画像の複数のチャンネルの対応する画素の値のうちの最小値を選択することにより、MPX画像の複数のチャンネルが非線形に合成される。この演算は、例えば、全ての画素位置(x,y)について、b(x,y)=min[a1(x,y)、a2(x,y)、...、aM(x,y)]と表され得る(ここで、aはMPX画像を示し、bは自己蛍光画像を示す)。この動作は線形ではない。自己蛍光画像を生成することはまた、非線形結合から生じる画像に対して1つまたは複数の演算を実行することを含み得る(例えば、メディアンフィルタを適用するなどのノイズ低減演算)。
【0048】
図9は、いくつかの実施形態にかかる自己蛍光画像を生成するプロセスの例を示している。このような処理は、(例えば、様々なチャンネルで記録された組織の自己蛍光のレベル間の変動を補償するために)MPX画像のチャンネルを等化することを含むことが望ましい場合がある。
図9に示す例では、(例えば、上述したように最小画素値を選択することによって)MPX画像のM個のチャンネルを非線形合成する前に等値化して自己蛍光画像を生成している。
図9に示すように、MPX画像の様々なチャンネルに適用される等化は、自己蛍光の基準スペクトル分布にしたがって決定され得る。例えば、等化は、MPX画像の各チャンネルを対応する重み(例えば、スケーリング係数)で除算することによって実行されてもよく、様々なチャンネルのそれぞれの重みの値は、自己蛍光の所与の自己蛍光基準スペクトル分布にしたがって決定され得る。基準スペクトル分布は、MPX画像内の試料の組織タイプに応じて選択されてもよく、自己蛍光基準ベクトルの要素としての様々なチャンネルの重みの計算は、以下により詳細に説明される。
【0049】
図9に示す特定の例では、Mは6に等しく、6つのチャンネルのそれぞれは、6つのフルオロフォアの異なる1つに対応する(それらの対応するフルオロフォアおよび1つの蛍光発生核対比染色剤(DAPI)を有する5つのバイオマーカー)。そのようなMPX画像を作製するために使用される染色プロトコルは、5プレックスアッセイまたはパネルとも呼ばれ得る。他の例では、Mは、3、4、または5に等しくてもよく(例えば、それぞれ2プレックス、3プレックス、および4プレックスのアッセイまたはパネルに対応する)、または6よりも大きくてもよい(例えば、7、8、9、または10以上)。
図9に示す特定の例では、4とラベル付けされたチャンネルは、DAPIチャンネルである。
【0050】
図24は、異なる4種類の組織(乳癌、結腸直腸癌(CRC)、膵臓癌、および胃癌)についての自己蛍光の基準分光分布と、他の6チャンネルシステム(この図では、6とラベル付けされたチャンネルは、DAPIチャンネルである)についてのこれらの分布の平均との比較例を示している。基準スペクトル分布は、MPX画像に示される染色試料と同じ種類の組織(例えば、同じ腫瘍から)である無染色試料の基準MPX画像から取得され得る。
図25は、染色されていない組織試料のMPX画像の2つの例を示している。基準MPX画像は、例えば、染色された試料のMPX画像を取り込むために使用されるのと同じイメージングシステムを使用して取り込まれ得る。
【0051】
基準スペクトル分布は、自己蛍光基準ベクトル(例えば、M要素基準ベクトル)として取得され得、ベクトルの各要素は、基準MPX画像の異なるチャンネルに対応し、そのチャンネルにおける組織自己蛍光の相対的な程度を表す。相対的な程度は、例えば、線形、多項式、指数(例えば、大きさ、エネルギー、または電力)であり得る。
図9に示す例では、(上述したようにMPX画像の対応するチャンネルに重みとして適用され得る)自己蛍光基準ベクトルの6つの要素の値は、[0.0919185036025541 0.147589073045573 0.467635814185659 0.822229790703801 0.195309012162166 0.120506922298875]である。自己蛍光の相対的な程度が決定される基準MPX画像の特定の領域は、画像全体または画像の選択された関心領域(例えば、腫瘍領域または他の注釈付き領域)であり得る。自己蛍光の相対的な程度のそれぞれは、所望の画像領域にわたる対応するチャンネルの値の選択された統計量(例えば、算術平均、中央値などの平均)の値として計算され得る。
【0052】
統計値を計算する前に、所望の画像領域にわたって基準MPX画像の各チャンネルを正規化することが望ましい場合がある。そのような一例では、各チャンネルは、最小-最大スケーリングによって(すなわち、所望の画像領域にわたってチャンネルの最大および最小画素値を識別し、最小画素値を減算し、その結果を最小画素値と最大画素値との差で除算することによって領域内のチャンネルの各値を正規化することによって)0から1の範囲に別々に正規化される。
【0053】
ソースMPX画像内のチャンネル数がMよりも大きいことも可能であるが、MPX画像のチャンネルのうちの1つずつからの情報を使用してNチャンネル入力画像を生成し、同様にNチャンネル訓練画像のそれぞれについても同様にすることが望ましい場合がある。換言すれば、MがソースMPX画像内のチャンネルの総数に等しいことが望ましい場合がある。訓練データを生成するためにソースMPX画像の全てのチャンネルからの情報を最大限に使用することは、対応する標的画像との構造的類似度が最大の合成画像を生成するように訓練された生成器ネットワークを取得するために重要であり得る。同様に、自己蛍光画像を生成するために組み合わされるチャンネルの数はM未満であることが可能であるが、自己蛍光画像を生成するためにM個のチャンネルの全てを組み合わせることが望ましい場合がある。
【0054】
次に、マッピングされた入力画像の例について説明する。
図10A~
図10Dは、3チャンネル(RGB)画像の4つの例を示し、それぞれは、いくつかの実施形態にかかる異なる対応するマッピングにしたがって、IHC染色組織切片のMPX画像の6個のチャンネルからの情報に基づいて取得される。これらの4つの例のそれぞれは、自己蛍光画像とMPX画像のDAPI(または他の核対比染色)チャンネルからの画像との異なる線形結合に基づいている。本明細書の説明では、
図10Aに示すようなマッピングは、「AF+AF+DAPI」またはAADとして識別され、
図10Bに示すようなマッピングは、「AF+DAPI+DAPI」またはADDとして識別され、
図10Cに示すようなマッピングは、「擬似H&E」またはpHEとして識別され、
図10Dに示すようなマッピングは、「擬似H&E 2」またはpHE2として識別される。
【0055】
(例えば、
図11~
図13を参照すると)本明細書に記載の自己蛍光画像とDAPI画像との線形結合を実行する前に、自己蛍光画像とDAPI画像とを指定された範囲に再スケーリングする(あるいは、自己蛍光画像とMPX画像とを指定された範囲に再スケーリングする)ことが望ましい場合がある。一例では、再スケーリングは、画像の最小および最大画素値を識別し、(例えば、本明細書に記載の最小-最大スケーリングによって)画像の画素値を0から1の範囲に正規化することによって実行され得る正規化であり得る。別の例では、再スケーリングは、画像の各画素値に所望のスケーリング係数を適用することによって実行されてもよい(例えば、画像の画素値を0から65535の範囲から0から2000の範囲、または0から3000の範囲などに再スケーリングする)。
【0056】
(例えば、
図11~
図13を参照すると)本明細書で説明した自己蛍光画像とDAPI画像との線形結合を実行する前に、(例えば、組織の自己蛍光を補償するために)DAPI画像を補正したり、MPX画像を補正したりすることが望ましい場合がある。このような補正は、上述したように、DAPIやMPX画像の再スケーリングの前に実行されてもよく、後に実行されてもよい。DAPI画像の補正は、自己蛍光画像に自己蛍光基準ベクトルの対応する要素を乗算して加重自己蛍光画像を取得し、DAPI画像の各画素値から加重自己蛍光画像の対応する画素値を減算して補正DAPI画像を取得することによって実行されてもよい。MPX画像の各チャンネルの補正も同様に、自己蛍光画像に自己蛍光基準ベクトルの対応する要素を乗算して加重自己蛍光画像を取得し、チャンネルの各画素値から加重自己蛍光画像の対応する画素値を減算して補正チャンネル画像を取得することによって実行されてもよい。
図26は、補正前(左上)および補正後(右下)の膵臓組織試料のMPX画像の例を示している。
【0057】
図11および
図12に示すように、マッピングされた入力画像は、自己蛍光画像とDAPI画像との線形結合として生成され得る。
図11は、いくつかの実施形態にかかる、3チャンネルマッピングされた入力画像を生成するプロセスの例を示している。この例(
図10Aに示すようなマッピングされた入力画像「AAD」または「AF+AF+DAPI」に対応する)では、入力画像の第1のチャンネルが自己蛍光画像であり、入力画像の第2のチャンネルも自己蛍光画像であり、入力画像の第3のチャンネルがDAPI画像である。AAD画像の3つのチャンネルは、GANフレームワークの訓練データとして実際のH&E画像のカラーチャンネル(例えば、RGB)と対にされ得る。RGBの3つのチャンネルは、AAD画像の3つのチャンネルにマッピングされる。
【0058】
図12は、いくつかの実施形態にかかる、3チャンネルマッピングされた入力画像を生成するプロセスの別の例を示している。この例(
図10Bに示すようなマッピングされた入力画像「ADD」または「AF+DAPI+DAPI」に対応する)では、入力画像の第1のチャンネルが自己蛍光画像であり、入力画像の第2のチャンネルがDAPI画像であり、入力画像の第3のチャンネルもDAPI画像である。ADD画像の3つのチャンネルのそれぞれは、GANフレームワークのための訓練データとしての実際のH&E画像の対応するカラーチャンネルと対にされ得る。例えば、ADD画像の3つのチャンネルは、実際のH&E画像のRチャンネル、Gチャンネル、およびBチャンネルにそれぞれマッピングされ得る。
【0059】
あるいは、マッピングされた入力画像は、自己蛍光画像とDAPI画像との線形結合である光学濃度配列から生成されてもよい。そのような計算は、例えば、光が通って移動する材料による光の減衰を記述するランベルト・ベールの法則に基づいてもよい。所与の波長cおよび所与の染色剤sについて、この関係は、以下のように表され得る:[波長cでの透過光強度]=[波長cでの入射光強度]*exp(-[波長cでの染色組織の光学密度])、ここで、ODc=(組織中の染色剤sの量)*(波長cでの染色剤sの吸収係数)である。
【0060】
プロセスブロック804のそのような実施形態は、自己蛍光画像とDAPI画像との線形結合として光学濃度値の配列(例えば、光学密度値のNチャンネル配列)を生成することと、マッピングされた入力画像を取得するために光学濃度値の配列に非線形関数を適用すること(「光学濃度配列」とも呼ばれる)とを含み得る。自己蛍光画像とDAPI画像との線形結合として光学濃度配列を生成することは、例えば以下のように実行され得る:
OD_R=DAPIwt[1]*DAPI+afwt[1]*AF;
OD_G=DAPIwt[2]*DAPI+afwt[2]*AF;
OD_B=DAPIwt[3]*DAPI+afwt[3]*AF;
ここで、OD_R、OD_G、およびOD_Bは、光学濃度配列の赤、緑、および青のチャンネルを示し、DAPIは、DAPI画像を示す。AFは、自己蛍光画像を示す。DAPIwtおよびafwtは、3要素重みベクトルを示す。場合によっては、(例えば、隣接スペクトル内のブリードスルー密度を減算するために)光学密度配列の1つまたは複数のチャンネルを調整することが望ましい場合がある。
【0061】
DAPIwtおよびafwtの値は、実験的に選択され得る(例えば、得られたマッピングされた画像の目視検査およびこれらの画像と実際のH&E画像との比較によって)。そのような一例では、
図10Cに示すような入力画像pHEの光学濃度配列は、重みベクトルDAPIwt=[29.2500 48.7500 14.6250]およびafwt=[0.3900 23.5535 16.6717]を使用して取得される。
【0062】
あるいは、DAPIwtおよびafwtの値は、ヘマトキシリンおよびエオシンの色依存性吸収特性にしたがって選択され得る。そのような一例では、
図10Dに示すような入力画像pHE2の光学濃度配列は、重みベクトルDAPIwt=w_D*Abs_Hemおよびafwt=w_A*Abs_Eosを使用して取得され、式中、Abs-Hemは、その要素が赤色、緑色、および青色光(例えば、[0.65 0.70 0.29])に対するヘマトキシリンの相対吸光度を示す色ベクトルであり、Abs_Eosは、その要素が赤色、緑色、および青色光(例えば、[0.07 0.99 0.11])に対するエオシンの相対吸光度を示す色ベクトルであり、スカラー重みw_Dとw_Aとの和は1である(デフォルト値w_D=w_A=0.5)。
図13は、いくつかの実施形態にかかる、3チャンネルマッピングされた入力画像を生成するそのようなプロセスの例を示している。場合によっては、w_Aに対してw_Dの値を増加させることが望ましい場合がある(例えば、マッピングされた入力画像における細胞核の影響を増大させるために)。
【0063】
光学濃度値の配列からマッピングされた入力画像を取得することは、非線形関数(例えば、指数関数)を適用して、光学濃度値の配列を透過値の配列(「透過配列」とも呼ばれる)に変換することを含み得る。このような変換は、例えば、以下のように実行され得る:
T_R=exp(-OD_R);
T_G=exp(-OD_G);
T_B=exp(-OD_B);
ここで、T_R、T_G、およびT_Bは、透過値の配列の赤、緑、および青のチャンネルを示す。非線形関数の別の例では、指数関数の底は、ネイピア数eではなく10である。
【0064】
透過配列の値を所望の画像空間に再スケーリングして、マッピングされた入力画像を明視野画像として取得し得る。一例では、透過配列の各値は、マッピングされた入力画像を24ビットRGB画像として取得するために255倍にスケーリング(例えば、乗算)される。別の例では、透過配列の値は、最小-最大スケーリングによって所望の範囲(例えば、0から255の範囲)に正規化される。マッピングされた入力画像の他の出力画像空間は、Lab色空間(例えば、CIELAB、ハンターLabなど)またはYCbCr色空間(例えば、YCbCr、Y’CbCrなど)を含み得、透過配列からマッピングされた入力画像を生成することは、ある色空間から別の色空間への変換(例えば、RGBまたはsRGBからYCrCbまたはY’CrCb)を含み得る。
【0065】
図14は、様々な実施形態にかかる、少なくとも1つの組織化学的染色剤によって染色された組織切片を示す合成画像に組織切片のMPX画像を変換するための例示的なコンピューティング環境1400(すなわち、データ処理システム)を示している。
図14に示すように、この例においてコンピューティング環境1400によって実行される変換は、いくつかの段階、すなわち画像記憶段階1405、前処理段階1490、モデル訓練段階1410、変換段階1415、および分析段階1420を含む。画像記憶段階1410は、予め選択された領域または生物学的試料スライド(例えば、組織スライド)の全体からデジタル画像1430のソースセットおよびデジタル画像1435の標的セットを提供するためにアクセスされる(例えば、前処理段階1490によって)1つまたは複数のデジタル画像スキャナまたはデータベース1425を含み得る。
【0066】
モデル訓練段階1410は、他の段階によって使用されるべき1つまたは複数のモデル1440a~1440n(「n」は任意の自然数を表す)(本明細書では個別にモデル1440と呼ばれてもよく、まとめてモデル1440と呼ばれてもよい)を構築して訓練する。モデル1440は、畳み込みニューラルネットワーク(「CNN」)、初期ニューラルネットワーク、残差ニューラルネットワーク(「Resnet」)、U-Net、V-Net、シングルショットマルチボックス検出器(「SSD」)ネットワーク、リカレントニューラルネットワーク(「RNN」)、深層ニューラルネットワーク、正規化線形ユニット(「ReLU」)、ロングショートタームメモリ(「LSTM」)モデル、ゲーティッドリカレントユニット(「GRU」)モデルなど、またはそれらの任意の組み合わせを含み得る、機械学習(「ML」)モデルとすることができる。様々な実施形態では、生成モデルは、敵対的生成ネットワーク(「GAN」)、深層畳み込み敵対的生成ネットワーク(「DCGAN」)、変分自己符号化器(VAE)、隠れマルコフモデル(「HMM」)、ガウス混合モデル、ボルツマンマシンなど、またはそのような技術のうちの1つもしくは複数、例えばVAE-GANの組み合わせなどの、教師なし学習を使用して任意の種類のデータ分布を学習することができるモデル1440を訓練することによって学習されたパラメータによって構成される。コンピューティング環境1400は、ソース画像を生成画像に変換するために、同じ種類のモデルまたは異なる種類のモデルを使用し得る。特定の事例では、生成モデルは、出力画像が現実であるか偽であるかを分類しようとする損失関数によって構築されたGANであるモデル1440を訓練することによって学習されたパラメータによって構成され、同時にこの損失を最小限に抑えるように生成モデルを訓練する。
【0067】
図15に示す例示的な実施形態では、学習されたパラメータを提供するように訓練されたモデル1440は、GANモデルの拡張である条件付きGAN(「cGAN」)1500の一部であり、特定の条件または属性を有する画像を生成する。cGANは、出力の結合構成にペナルティを課す構造化損失を学習する。
図15を参照すると、cGAN1500は、生成器1510および識別器1515を含む。生成器1510は、ランダムに生成されたノイズベクトル1520および潜在特徴ベクトル(または1次元ベクトル)1525(条件、例えば、本例では、マッピングされた入力画像)を入力データおよび識別器1515からのフィードバックとして受け取り、できるだけ現実の標的画像1535に近い新たな画像1530を生成するニューラルネットワーク(例えば、CNN)である。識別器1515は、生成器1510からの生成画像1530が現実画像であるか偽画像であるかを識別するための分類器として構成されたニューラルネットワーク(例えば、CNN)である。潜在特徴ベクトル1525または条件は、入力画像または入力画像のセット1540(例えば、マッピングコントローラ1593によって提供される画像)から導出され、入力画像1540から期待されるクラス(例えば、選択されたマッピングにしたがって対応するMPX画像のM個のチャンネルから取得されたNチャンネル画像)または特定の特性のセットを符号化する。ランダムに生成されたノイズベクトル1520は、ガウス分布から生成されてもよく、ベクトル空間は、領域にとって重要であるが直接観測可能ではない潜在変数または隠れ変数から構成されてもよい。潜在特徴ベクトル1525およびランダムノイズベクトル1520は、生成器1510への入力1545として組み合わせられ得る。代替的または追加的に、ノイズは、ドロップアウト(例えば、入力を層に確率的に落とす)の形態で生成器1510内に追加されてもよい。
【0068】
生成器1510は、結合入力1545を受信し、問題領域(すなわち、組織化学的に染色された標的画像1535に関連する特徴の領域)内の潜在特徴ベクトル1525およびランダムノイズベクトル1520に基づいて画像1530を生成する。識別器1515は、標的画像1535および生成画像1530の両方を入力として条件画像分類を実行し、生成画像1530が標的画像1535の実際の翻訳であるか偽の翻訳であるかの可能性を予測する(1550)。識別器1515の出力は、生成画像1530のサイズに依存するが、単一の値または値の二乗活性化マップであってもよい。各値は、生成画像1530内のパッチが実在する尤度の確率である。これらの値を平均化して、必要に応じて全体的な尤度または分類スコアを得ることができる。生成器1510および識別器1515の両方の損失関数は、損失が、生成された画像1530が実際のものであるか、または標的画像1535の偽の翻訳であるかの可能性を予測する(1550)そのジョブを識別器1515がどれだけ良好に実行するかに依存するように構成され得る。十分な訓練の後、生成器1510は、より標的画像1535に似た生成画像1530を生成し始める。GAN1500の訓練は、事前定義された数の訓練インスタンスに対して進行し得、訓練または検証セットを使用して決定された1つまたは複数の性能測定基準(例えば、精度、適合率、および/または再現率)が対応する閾値を超える限り、結果として得られる学習されたパラメータが受け入れられ得る。あるいは、GAN1500の訓練は、最近の訓練反復に関連する1つまたは複数の性能測定基準が対応する閾値を超えるまで進行し得る。この時点で、生成された画像1530は、識別器がもはや現実と偽とを識別することができないように、標的画像1535に十分に類似し得る。生成器ネットワーク1510が訓練されると、入力されたNチャンネル画像のセット(選択されたマッピングにしたがって対応するMPX画像のM個のチャンネルから取得される)をGAN1500に入力して、入力された画像のセットを、組織化学的に染色されたスライドから取得された画像の標的セット(例えば、H&E染色)と同様の特性を有する新たな生成された画像のセットに変換し得る。その後、新たに生成された画像のセットは、病理学者によって、(例えば、MPX画像または合成H&E画像をオンおよびオフに切り替えることによって)MPX画像の注釈付け中に組織構造を視覚化するため、MPX画像に基づく診断を確認するため、および/または(例えば、バイオマーカーの複合パネルに基づく分析を支援するために)あるMPX画像を別のMPX画像に位置合わせするためなどに使用され得る。
【0069】
図14に戻って参照すると、この例ではモデル1440を訓練するために、前処理段階1490は、デジタル画像(デジタル画像1430のソースセットおよびデジタル画像1435の標的セット)を取得し、(例えば、マッピングコントローラ1493によって)選択されたマッピングにしたがってソース画像をNチャンネル画像にマッピングし、画像を訓練のための画像1445aのペアワイズサブセット(マッピングされた画像と標的画像の少なくとも一対)(例えば、90%)および検証のための画像1445bのペアワイズサブセット(例えば、10%)に分割し、画像1445aのペアワイズサブセットおよび画像1445bのペアワイズサブセットを前処理し、画像1445aのペアワイズサブセットを拡張し、場合によってはラベル1450によって画像1445aのペアワイズサブセットに注釈を付けることによって試料1445を生成する。画像のペアワイズサブセット1445aは、例えば、データベースまたは画像サーバなどのデータ記憶構造から取得され得る。各画像は、組織などの生物学的試料を示す。
【0070】
分割は、ランダムまたは擬似ランダム(例えば、90%/10%、80%/20%、または70%/30%を使用する)に実行されてもよく、または分割は、サンプリングバイアスおよびオーバーフィッティングを最小限に抑えるために、K倍交差検証、リーブワンアウト交差検証、リーブワングループアウト交差検証、ネステッド交差検証などのより複雑な検証技術にしたがって実行されてもよい。前処理は、各画像が単一の関心のあるオブジェクトのみを含むように画像をトリミングすることを含み得る。場合によっては、前処理は、全ての特徴を同じスケール(例えば、同じサイズスケールまたは同じカラースケールまたは彩度スケール)にするための標準化または再スケーリング(例えば、正規化)をさらに含み得る。特定の例では、画像は、所定の画素(例えば、2500画素)の最小サイズ(幅または高さ)または所定の画素(例えば、3000画素)の最大サイズ(幅または高さ)でサイズ変更され、元のアスペクト比を維持する。
【0071】
例えば、前処理段階1490は、ソースセットおよび標的セットからの複数のパッチ画像を訓練データ用の画像の1つまたは複数のペアワイズサブセットとして準備し得る。対の画像の準備は、ソース画像と標的画像の一致した対にアクセスすることを含み得、各一致した対のソース画像および標的画像は、同じ生物学的試料の近くの(例えば、隣接している、またはほぼ隣接している)切片のスライドからのものであり(例えば、腫瘍試料)、ソース画像に示される切片は、複数の(例えば、2、3、4、5、または6以上)選択されたIHCフルオロフォアによって染色されており、標的画像に示される切片は、1つまたは複数の選択された組織化学的染色剤によって染色されている。1つの非限定的な例では、ソース画像のそれぞれの切片は、6つのフルオロフォア(核対比染色剤、例えばDAPIを含む)によって染色されており、標的画像のそれぞれの切片は、H&Eによって染色されている。
図5Aおよび
図5Bは、IHC染色された組織切片の暗視野MPX画像(
図5A)と、H&E染色された同じ組織試料の近傍(例えば、隣接する)切片の明視野画像(
図5B)との一致した対の一例を示している。
【0072】
次に、前処理段階1490は、(例えば、マッピングコントローラ1493によって)選択されたマッピングにしたがってソース画像のそれぞれをマッピングして、対応する標的画像と対になるNチャンネル入力画像を取得し、対になった入力画像および標的画像(例えば、スライド画像全体)のそれぞれを所定のサイズ(例えば、128×128、256×256、または別のサイズ)のいくつかのパッチに分割して、訓練用のパッチの一致した対を生成し得る。例えば、精査する病理学者によって追加された腫瘍注釈など、画像内の関心領域からのパッチのみを使用することが望ましい場合がある。前処理段階1490(例えば、位置合わせコントローラ1495)は、画像がパッチに分割される前および/または後に、対になった入力画像および標的画像の位置合わせおよび/または位置合わせを実行し得る。位置合わせは、一方の画像を固定画像とも呼ばれる基準画像として指定することと、他方の画像が基準画像と位置合わせされるように、他方の画像に幾何学的変換または局所変位を適用することとを含み得る。組織化学的に染色された画像(すなわち、標的画像)は、ネットワークを訓練するためのグラウンドトゥルースを提供することから、位置合わせおよび位置合わせの目的で標的画像を基準画像として指定することが望ましい場合がある。入力セットおよび標的セットからの位置合わせされたパッチの対が選択され、このプロセスは、訓練データ用の画像の1つまたは複数のペアワイズサブセットをもたらす。前処理段階1490は、深層学習ネットワークを訓練するために、パッチ対をGANまたはcGANに入力され得る。
【0073】
図14に戻って参照すると、前処理段階1490は、データセット内の画像の修正バージョンを作成することによって、画像1445aのペアワイズサブセットのサイズを人工的に拡張するために拡張を使用し得る。画像データ拡張は、元の画像と同じクラスに属するデータセット内の画像の変換バージョンを作成することによって実行され得る。変換は、シフト、フリップ、ズームなどの画像操作の分野からの動作の範囲を含む。いくつかの例では、動作は、モデル1440が画像のペアワイズサブセット1445aから利用可能な状況外の状況下で実行することができることを保証するために、ランダム消去、シフト、輝度、回転、ガウスぼかし、および/または弾性変換を含む。
【0074】
モデル1440の訓練プロセスは、画像のペアワイズサブセット1445aからモデル1440に画像を入力して、モデル1440の1つまたは複数の損失関数または誤差関数(例えば、画像訓練データの確率を最大化するように識別器を訓練するための第1の損失関数と、生成器からサンプリングされた生成画像の確率を最小化するように識別器を訓練し、識別器がそれ自体の生成画像に割り当てる確率を最大化するように生成器を訓練するための第2の損失関数)を最小化するモデルパラメータのセット(例えば、重みおよび/またはバイアス)を見つける反復動作を実行することと、を含む。訓練の前に、モデル1440は、モデルの外部にあり、モデル1440の挙動を制御するように調整または最適化され得る設定であるハイパーパラメータの定義されたセットによって構成される。ほとんどのモデルでは、ハイパーパラメータは、例えばメモリまたは実行コストなど、モデルの異なる態様を制御するために明示的に定義される。しかしながら、モデルを特定のシナリオに適合させるために、追加のハイパーパラメータが定義され得る。例えば、ハイパーパラメータは、モデルの隠れユニットの数、モデルの学習率、畳み込みカーネル幅、またはモデルのカーネルの数を含み得る。訓練の各反復は、モデルパラメータのセットを使用する損失関数または誤差関数の値が、前の反復におけるモデルパラメータの異なるセットを使用する損失関数または誤差関数の値よりも小さくなるように、(ハイパーパラメータの対応する定義されたセットによって構成されるような)モデル1440のモデルパラメータのセットを見つけることを含むことができる。ハイパーパラメータとは対照的に、モデルパラメータは、モデルの内部にある変数であり、それらの値は、訓練データから推定(例えば、学習)され得る。損失関数または誤差関数は、モデル1440を使用して推論された出力と、ラベル1450を使用したグラウンドトゥルース標的画像との間の差を測定するように構築され得る。
【0075】
モデルパラメータのセットが識別されると、モデル1440は、訓練されており、画像のペアワイズサブセット1445b(試験または検証データセット)を使用して検証され得る。検証プロセスは、ハイパーパラメータを調整し、最終的に最適なハイパーパラメータのセットを見つけるために、K-分割交差検証、一個抜き交差検証、一群抜き交差検証、入れ子交差検証などの検証技術を使用して、画像のペアワイズサブセット1445bからモデル1440に画像を入力する反復動作を含む。最適なハイパーパラメータのセットが取得されると、画像のサブセット1445bからの画像の予約された試験セットをモデル1440に入力して出力(この例では、標的画像と同様の特性を有する生成画像)を取得し、Bland-Altman法およびスピアマンのランク相関係数などの相関技術を使用し、誤差、精度、適合率、再現率、受信者動作特性曲線(ROC)などの性能メトリックを計算して、出力がグラウンドトゥルース標的画像に対して評価される。
【0076】
理解されるべきであるように、他の訓練/検証機構が想定され、コンピューティング環境1400内に実装されてもよい。例えば、モデル1440は、訓練されてもよく、ハイパーパラメータは、画像1445aのペアワイズサブセットからの画像上で調整されてもよく、画像1445bのペアワイズサブセットからの画像は、モデル1440の性能を試験および評価するためにのみ使用されてもよい。
【0077】
モデル訓練段階1410は、1つまたは複数の訓練された変換モデル1460および任意に1つまたは複数の画像分析モデル1465を含む訓練されたモデルを出力する。場合によっては、第1のモデル1460aは、生物学的標本(組織切片)のソース画像1430から取得された入力画像を処理するように訓練されている。入力画像は、選択されたマッピングにしたがって生物学的試料のMPX画像のM個のチャンネルから取得されるNチャンネル画像である。入力画像は、変換段階1415内の変換コントローラ1470によって取得される。変換コントローラ1470は、1つまたは複数の訓練された変換モデル1460を使用して、入力画像を標的画像の特性を有する新たな画像1475に変換するためのプログラム命令を含む。標的画像の特徴は、1つまたは複数の選択された組織化学的染色剤(例えば、H&E)によって染色された組織切片の画像に関連付けられる。変換は、入力画像を訓練された生成器モデル(変換モデル1460の一部)に入力することと、生成器モデルによって新たな画像1475を生成することとを含む。
【0078】
場合によっては、新たな画像1475は、分析段階1420内の分析コントローラ1480に送信される。分析コントローラ1480は、1つまたは複数の画像分析モデル1465を使用して、新たな画像1475内の生物学的試料を分析し、分析に基づく分析結果1485を出力するためのプログラム命令を含む。新たな画像1475内の生物学的試料の分析は、新たな画像1475内の領域、新たな画像1475内の1つまたは複数の細胞、および/または細胞以外の新たな画像1475内のオブジェクトに基づいて測定値を抽出することを含み得る。面積ベースの測定は、最も基本的な評価、例えば、特定の染色(例えば、組織化学的染色)の面積(2次元)、脂肪小胞の面積、またはスライド上に存在する他の事象を定量化することを含み得る。細胞ベースの測定は、オブジェクト、例えば細胞の識別および列挙を目的とする。個々の細胞のこの特定は、細胞内区画のその後の評価を可能にする。最後に、個々の細胞から構成され得ない組織切片上に存在する事象またはオブジェクトを評価するために、アルゴリズムが利用され得る。特定の例では、イメージング分析アルゴリズムは、細胞または細胞内構造の位置を特定し、最終的に診断および予測を支援するために使用され得る細胞染色、形態、および/または構造の定量的表現を提供するように構成される。場合によっては、イメージング分析アルゴリズムは、標的画像(例えば、組織化学的に染色された切片の画像)の特性を有する画像の分析用に特に構成される。例えば、新たな画像1475の分析は、対応するMPXソース画像に対して、ROISおよび/または除外されるべき領域の自動描写を提供することを含み得る。別の例では、新たな画像1475の分析は、新たな画像1475を別の新たな画像1475に位置合わせすることと、位置合わせ結果を対応するMPXソース画像に適用することとを含み得る。
【0079】
明示的に示されていないが、コンピューティング環境1400は、開発者に関連付けられた開発者デバイスをさらに含み得ることが理解されよう。開発者デバイスからコンピューティング環境1400の構成要素への通信は、モデルに使用される入力画像のタイプ、使用されるモデルの数およびタイプ、各モデルのハイパーパラメータ、例えば学習率および隠れ層の数、データ要求をフォーマットする方法、使用される訓練データ(例えば、および訓練データへのアクセス方法)および使用される検証技術、および/またはコントローラプロセスを構成する方法を示し得る。
【0080】
生成器ネットワーク1510を訓練するために使用され得るcGANモデル1500の特定の一例は、Pix2Pix GANである。
図16は、選択されたマッピングにしたがってIHC-染色された組織切片の対応するMPX画像のM個のチャンネルからそれぞれ取得されたN-チャンネル画像を、組織化学的に染色された(例えば、H&E染色された)組織切片の合成画像に変換するように生成器ネットワーク1610を訓練するためにPix2Pix GAN1600を使用するcGANモデル1500の実装の例を示している。
図16に示すように、生成器ネットワーク1610は、U-Netアーキテクチャを使用して実装され、これは、入力をボトルネック層に漸進的にダウンサンプリングする層を有するエンコーダと、出力を生成するためにボトルネック出力を漸進的にアップサンプリングする層を有するデコーダとを含む。
図16に示すように、U-Netはまた、等しいサイズの特徴マップを有するエンコーダ層とデコーダ層との間のスキップ接続を含む。これらの接続は、エンコーダ層の特徴マップのチャンネルを対応するデコーダ層の特徴マップのチャンネルと連結する。特定の例では、生成器ネットワーク1610は、生成された画像と予想される出力画像(例えば、
図16の「予測画像」および「グラウンドトゥルース」)との間で測定されたL1損失を介して更新される。
【0081】
一般に、Pix2Pix GANを使用することは、生成器ネットワークを訓練するために使用されるべき画像パッチの一致した対が(例えば、画素レベルで)位置合わせされている必要がある。
図17は、(例えば、訓練および/または検証のために)本明細書に記載の一致した画像対から一致し、位置合わせされた画像パッチの対を生成する例示的なプロセス1700のフローチャートを示している。プロセス1700は、前処理段階1490によって(例えば、位置合わせコントローラ1495によって)実行され得る。
図17を参照すると、ブロック1704において、一致した画像対の低解像度バージョンが粗に位置合わせされる。一例では、そのような粗い位置合わせは、変換行列M
C(例えば、並進および/または回転を含む)を、IHC染色された切片のMPX画像のM個のチャンネルのNチャンネルマッピングI
Nmapの低解像度バージョンに適用して、それをH&E染色された一致する(例えば、近い)切片の画像I
H&Eに位置合わせすることによって実行され得る。変換行列Mcは、例えば、位置合わせされる各画像における組織の輪郭に基づいて自動的に算出され得る。
【0082】
ブロック1708において、(例えば、注釈をカバーする各画像上にグリッドを投影し、各画像から対応するタイルを抽出することによって)粗に位置合わせされた画像対の関心領域(ROI)からタイルが抽出される。
図18は、画像I
NmapのROIからのタイルT
Nmap(例えば、サイズが2048×2048画素)および画像I
H&Eからの対応するタイルT
H&E(例えば、サイズが2048×2048画素)の例を示している。ブロック1712において、抽出されたタイル対のフル解像度バージョンが精密に位置合わせされる。そのような精密な位置合わせは、例えば、一方の画像(例えば、I
Nmap)からのタイルを他方の画像(例えば、I
H&E(上述した基準画像))からの対応するタイルに位置合わせするために、一方の画像からのタイルのスケーリング、デスキュー、および/またはワープを含み得る。ブロック1716において、精密に位置合わせされたタイルは、位置合わせされた画像対を取得するために互いにスティッチングされ、ブロック1720において、位置合わせされた画像のそれぞれは、パッチにスライスされる(例えば、サイズが128×128画素、256×256画素、または別のサイズ)。そして、ブロック1724において、位置合わせされた各画像(例えば、
図18に示すようなパッチP
NmapおよびP
H&E)からの一致したパッチが組み合わされ(例えば、左右にスティッチングされる)、訓練画像を取得する。例えば、サイズ128×128画素のそれぞれである一対のマッチングパッチを互いにスティッチングして、サイズ128×256画素の訓練画像を取得し得る。または、サイズ256×256画素のそれぞれが、サイズ256×512画素の訓練画像を取得するために互いにスティッチングされてもよい。このようにして、位置合わせ画像対を取得するために各画像の精密に位置合わせされたタイルを互いにスティッチングし、訓練画像のセットを取得するために位置合わせ画像のそれぞれからの一致したパッチを組み合わせるプロセスは、対応するMPX画像のM個のチャンネルから取得されたNチャンネルマッピングされた画像を組織化学的に染色された部分の合成画像に変換するために生成器ネットワーク1510を訓練するためにPix2Pix GAN実装を使用するための訓練データのセットを取得するために実装され得る。
【0083】
生成器ネットワーク1510を訓練するために使用され得るcGANモデル1500の別の特定の例は、複数の生成器ネットワークおよび複数の識別器ネットワークを含むサイクルGANである。
図19は、サイクルGANにおける生成器ネットワークG
XおよびG
Yならびに識別器ネットワークD
XおよびD
Yの図を示している。この例では、Y領域は、組織化学的に染色された試料を表す画像に対応し、X領域は、選択されたマッピングにしたがって、IHC染色試料の対応するMPX画像のM個のチャンネルから得られたNチャンネル画像に対応する。
【0084】
図20は、本明細書で説明されるようなサイクルGANの適用における生成器ネットワークと識別器ネットワークとの間のフローを示している。サイクルGANは、X-Y生成器ネットワークG
Y2024(生成器ネットワーク1510として訓練される)を含み、組織化学的に染色された試料の画像をIHC染色試料のMPX画像のM個のチャンネルからのNチャンネルマッピングを描写する画像に変換するように構成および訓練されたY-X生成器ネットワークG
X2020も含む。生成器ネットワークG
X2020は、1つまたは複数の畳み込み層を含むことができ、U-netまたはV-netを含み得る。場合によっては、生成器ネットワークG
X2020は、特徴抽出エンコーダ、変換器、およびデコーダを含み、それぞれが1つまたは複数の畳み込み層を有する。生成器ネットワークG
X2020およびG
Y2024のアーキテクチャは同じであってもよい。
【0085】
サイクルGANは、IHC染色試料のMPX画像のM個のチャンネルからのNチャンネルマッピングを描写する実際の画像と偽の画像(例えば、実際のNチャンネル画像2012および偽のNチャンネル画像2016)とを識別する識別器ネットワークDX2032と、組織化学的染色試料を描写する偽の画像と実際の画像(例えば、実際の組織化学的染色画像2004および偽の組織化学的染色画像2008)とを識別する別の識別器ネットワークDY2028とを含む。識別器ネットワークDXおよびDYのそれぞれは、1つまたは複数の畳み込み層および活性化層を含み得、識別器ネットワークDXおよびDYのアーキテクチャは同じであり得る。
【0086】
サイクルGANの使用は、一致した画像の対の精密な位置合わせ(例えば、
図17および
図18を参照して本明細書に記載されるような画像I
NmapおよびI
H&Eの精密な位置合わせ)が訓練データを生成するために必要とされないという利点を有し得る。しかしながら、Pix2Pix GAN実装を使用して、対になった位置合わせ画像の画像パッチ上で生成器ネットワーク1510を訓練すると、より良好な結果が得られた。
【0087】
本開示にかかる方法は、MPX画像の画像をH&E染色試料の対応する合成画像に変換するために実装され得る。そのような方法は、例えば、病理学者がMPXソース画像の注釈付けを支援するために使用され得る。そのような方法は、実際のH&E染色を行うことなく、下にある組織構造の基準画像を提供するための高速スクリーニングプロセスの重要な部分として実装され得る。さらに、そのような「仮想染色」技術は、対応するMPX画像に基づく診断を確認するために使用されることもできる。さらに、本明細書に記載の画像変換の方法を使用して、バイオマーカーの拡張パネルの分析(例えば、共配置、共発現)を可能にするためのMPX画像の自動位置合わせを支援し得る。
【0088】
図21~
図23は、MPX画像を対応する合成H&E染色画像に変換するために(本明細書に記載のPix2Pix GANを使用して)プロセス800の実装を使用する3つの異なる例をそれぞれ示している。
図21~
図23のそれぞれは、
図10A~
図10Dに示されて本明細書に記載された4つの例示的なソース画像マッピング(pHE、pHE2、AAD、ADD)のそれぞれを使用して得られた結果を示している。
図21~
図23のそれぞれにおいて、上段は、Nチャンネルマッピングされたソース画像(「アプローチ」)を示し、中段は、対応する生成された合成H&E画像(「出力」)を示し、下段は、対応する実際のH&E画像(「標的」)を示す。出力画像は、標的画像と非常に類似していることが分かる。
【0089】
V.さらなる考察
本開示のいくつかの実施形態は、1つまたは複数のデータプロセッサを含むシステムを含む。いくつかの実施形態では、システムは、命令を含む非一時的コンピュータ可読記憶媒体であって、該命令が、1つまたは複数のデータプロセッサ上で実行されると、1つまたは複数のデータプロセッサに、本明細書に開示された1つまたは複数の方法の一部または全部および/または1つまたは複数のプロセスの一部または全部を実行させる、非一時的コンピュータ可読記憶媒体を含む。本開示のいくつかの実施形態は、1つまたは複数のデータプロセッサに、本明細書に開示された1つまたは複数の方法の一部または全部および/または1つまたは複数のプロセスの一部または全部を実行させるように構成された命令を含む非一時的機械可読記憶媒体において、有形に具現化されたコンピュータプログラム製品を含む。
【0090】
使用された用語および表現は、限定ではなく説明の用語として使用され、そのような用語および表現の使用において、示されて説明された特徴の均等物またはその一部を除外する意図はないが、特許請求の範囲に記載された発明の範囲内で様々な変更が可能であることが認識される。したがって、特許請求の範囲に記載された本発明は、実施形態および任意の特徴によって具体的に開示されているが、本明細書に開示された概念の変更および変形は、当業者によってあてにされてもよく、そのような変更および変形は、添付の特許請求の範囲によって定義される本発明の範囲内にあると見なされることを理解されたい。
【0091】
説明は、好ましい例示的な実施形態のみを提供し、本開示の範囲、適用可能性または構成を限定することを意図しない。むしろ、好ましい例示的な実施形態の説明は、様々な実施形態を実装するための可能な説明を当業者に提供する。添付の特許請求の範囲に記載の趣旨および範囲から逸脱することなく、要素の機能および配置に様々な変更が加えられ得ることが理解される。
【0092】
実施形態の完全な理解を提供するために、以下の説明において具体的な詳細が与えられる。しかしながら、これらの具体的な詳細なしで実施形態が実施され得ることが理解されよう。例えば、回路、システム、ネットワーク、プロセス、および他の構成要素は、実施形態を不必要に詳細に不明瞭にしないために、ブロック図形式の構成要素として示されてもよい。他の例では、実施形態を不明瞭にすることを避けるために、周知の回路、プロセス、アルゴリズム、構造、および技術が不必要な詳細なしに示されてもよい。
【手続補正書】
【提出日】2023-12-15
【手続補正1】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項1】
画像変換の方法であって、
組織切片の多重免疫蛍光(MPX)画像のMチャンネルのそれぞれからの情報に基づくNチャンネル入力画像を生成することであって、Mは正の整数であり、NはM以下の正の整数である、Nチャンネル入力画像を生成することと、
複数の画像対を含む訓練データセットを使用して訓練された生成器ネットワークを使用して前記Nチャンネル入力画像を処理することによって合成画像を生成することと
を含み、
前記合成画像が、少なくとも1つの組織化学的染色剤によって染色された組織切片を示し、
前記複数の画像対の各画像対について、前記対が、
組織の第1の切片のMPX画像から生成されたNチャンネル画像と、
前記少なくとも1つの組織化学的染色剤によって染色された前記組織の第2の切片の画像と
を含む、方法。
【請求項2】
前記Nチャンネル入力画像を生成することが、前記組織切片の前記MPX画像の前記Mチャンネルのそれぞれについて、前記チャンネルからの情報を前記Nチャンネル画像の前記Nチャンネルのうちの少なくとも1つにマッピングすることを含む、請求項1に記載の方法。
【請求項3】
前記マッピングすることが、前記組織切片の前記MPX画像の複数のチャンネルのそれぞれからの情報に基づく自己蛍光画像を生成することを含み、
前記Nチャンネル入力画像が、前記自己蛍光画像からの情報に基づく、請求項2に記載の方法。
【請求項4】
前記自己蛍光画像が、前記組織切片の前記MPX画像の前記複数のチャンネルの非線形結合に基づく、請求項3に記載の方法。
【請求項5】
前記自己蛍光画像が、前記複数のチャンネル間の自己蛍光のスペクトル分布に基づく、請求項
3に記載の方法。
【請求項6】
前記MPX画像の前記複数のチャンネルが、前記組織切片の前記MPX画像の前記Mチャンネルである、請求項
3に記載の方法。
【請求項7】
前記複数の画像対の各画像対について、前記対の前記Nチャンネル画像が、前記第1の切片の前記MPX画像のMチャンネルのそれぞれからの情報に基づく、請求項
1に記載の方法。
【請求項8】
前記合成画像がNチャンネル画像である、請求項
1に記載の方法。
【請求項9】
前記Nチャンネル入力画像および前記合成画像のそれぞれがRGB画像である、請求項8に記載の方法。
【請求項10】
前記組織切片の前記MPX画像が暗視野画像であり、前記合成画像が明視野画像である、請求項
1に記載の方法。
【請求項11】
前記少なくとも1つの組織化学的染色剤がヘマトキシリンおよびエオシンである、請求項
1に記載の方法。
【請求項12】
前記Nチャンネル画像が、前記自己蛍光画像と、核対比染色に対応する前記組織切片の前記MPX画像のチャンネルとに基づく、請求項
1に記載の方法。
【請求項13】
前記Nチャンネル画像が、前記自己蛍光画像と、核対比染色に対応する前記組織切片の前記MPX画像のチャンネルとの線形結合に基づく、請求項
1に記載の方法。
【請求項14】
前記Nチャンネル画像が、前記線形結合に基づく光学濃度値の配列に基づく、請求項
1に記載の方法。
【請求項15】
Nが3に等しい、請求項
1に記載の方法。
【請求項16】
Mが少なくとも4である、請求項
1に記載の方法。
【請求項17】
前記生成器ネットワークが、敵対的生成ネットワークの一部として訓練された、請求項
1に記載の方法。
【請求項18】
前記生成器ネットワークが、U-Netとして実装される、請求項
1に記載の方法。
【請求項19】
前記生成器ネットワークが、エンコーダ-デコーダネットワークとして実装される、請求項
1に記載の方法。
【請求項20】
前記生成器ネットワークが、前記生成器ネットワークによる画像と予想される出力画像との間で測定されたL1損失を介して更新される、請求項
1に記載の方法。
【請求項21】
ユーザによって、対象の診断を決定する方法であって、前記決定することが、請求項1から20のいずれか一項に記載の前記合成画像に基づく、方法。
【請求項22】
前記ユーザによって、(i)前記合成画像、および/または(ii)前記対象の前記診断に基づく化合物による処置を施すことをさらに含む、請求項21に記載の方法。
【請求項23】
1つまたは複数のデータプロセッサと、
命令を含む非一時的コンピュータ可読記憶媒体であって、前記命令が、前記1つまたは複数のデータプロセッサ上で実行されると、前記1つまたは複数のデータプロセッサに、請求項1~20のいずれか一項に記載の画像変換の方法を実行させる、非一時的コンピュータ可読記憶媒体と、を備える、システム。
【請求項24】
1つまたは複数のデータプロセッサに請求項1~20のいずれか一項に記載の画像変換の方法を実行させるように構成された命令を含む、非一時的機械可読記憶媒体に有形に具現化されたコンピュータプログラム製品。
【国際調査報告】