(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2024-09-26
(45)【発行日】2024-10-04
(54)【発明の名称】酵母または細菌を識別するための方法
(51)【国際特許分類】
C12Q 1/04 20060101AFI20240927BHJP
C12M 1/34 20060101ALI20240927BHJP
G01N 21/27 20060101ALI20240927BHJP
【FI】
C12Q1/04
C12M1/34 B
G01N21/27 A
(21)【出願番号】P 2020533199
(86)(22)【出願日】2018-12-20
(86)【国際出願番号】 FR2018053436
(87)【国際公開番号】W WO2019122734
(87)【国際公開日】2019-06-27
【審査請求日】2021-11-25
(32)【優先日】2017-12-21
(33)【優先権主張国・地域又は機関】FR
【前置審査】
(73)【特許権者】
【識別番号】304043936
【氏名又は名称】ビオメリュー
【氏名又は名称原語表記】BIOMERIEUX
(74)【代理人】
【識別番号】110002077
【氏名又は名称】園田・小林弁理士法人
(72)【発明者】
【氏名】ラレマンド, ジョーダン
(72)【発明者】
【氏名】ルルー, ドゥニ
(72)【発明者】
【氏名】プティ, マニュエル
【審査官】松田 芳子
(56)【参考文献】
【文献】特表2015-507182(JP,A)
【文献】特表2011-529187(JP,A)
【文献】特開2001-272399(JP,A)
【文献】Annual International Conference of the IEEE Engineering in Medicine and Biology Society,2015年,Vol.2015,p.759-762
【文献】Computers in Biology and Medicine,2017年09月01日,Vol.88,p.60-71
【文献】Proceedings of SPIE,2016年,vol.9887,98873L1-12
【文献】mBio,2013年,Vol.4, No.6,e00865-13
(58)【調査した分野】(Int.Cl.,DB名)
C12Q 1/00
JSTPlus/JMEDPlus/JST7580(JDreamIII)
CAplus/MEDLINE/EMBASE/BIOSIS(STN)
(57)【特許請求の範囲】
【請求項1】
微生物の特徴を明らかにするための方法であって、
- 波長範囲390nm~900nm
で、微生物
を照明することと、
- 前記照明した微生物によって反射した、または前記照明した微生物を通じて透過した光度を前記範囲内で獲得することと、
- 前記範囲内で獲得した光度に応じて、酵母または
細菌株であるものとして微生物を判定することと
を含み、
細菌株のグラムタイプおよび発酵特性の検出をさらに含み、
- 波長範囲390nm~900nm
で、前記
細菌株の少なくとも1つの細菌を照明することと、
- 前記照明した細菌によって反射した、または前記照明した細菌を通じて透過した光度を前記範囲内で獲得することと、
前記範囲内で獲得した光度に応じて、
細菌株のグラムタイプおよび発酵特性を判定することと
を含み、
グラムタイプおよび発酵特性の判定が、獲得した光度がグラム陰性タイプかつ発酵性の
細菌株の光度であるかどうかを予測する第1の分類の適用を含み、
グラムタイプおよび発酵特性の判定が、獲得した光度がグラム陽性タイプの
細菌株の光度であるかどうかを予測する第2の分類の適用、または、獲得した光度がグラム陰性タイプの
細菌株の光度であるかどうかを予測する第2の分類の適用を含み、
第1の分類が、波長範囲470nm~500nm
のスペクトルおよび波長範囲575nm~645nm
のスペクトルだけに応じて光度を区別し、
第2の分類が、波長範囲535nm~675nm
のスペクトルおよび波長範囲850nm~875nm
のスペクトルだけに応じて光度を区別する、
方法。
【請求項2】
微生物の特徴を明らかにするための方法であって、
- 波長範囲390nm~900nm
で、微生物
を照明することと、
- 前記照明した微生物によって反射した、または前記照明した微生物を通じて透過した光度を前記範囲内で獲得することと、
- 前記範囲内で獲得した光度に応じて、酵母または
細菌株であるものとして微生物を判定することと
を含み、
細菌株のグラムタイプおよび発酵特性の検出をさらに含み、
- 波長範囲390nm~900nm
で、前記
細菌株の少なくとも1つの細菌
を照明することと、
- 前記照明した細菌によって反射した、または前記照明した細菌を通じて透過した光度を前記範囲内で獲得することと、
前記範囲内で獲得した光度に応じて、
細菌株のグラムタイプおよび発酵特性を判定することと
を含み、
グラムタイプおよび発酵特性の判定が、獲得した光度がグラム陰性タイプかつ発酵性の
細菌株の光度であるかどうかを予測する第1の分類の適用を含み、
ここで、
- 照明および獲得が、
細菌株のコロニーと、血液寒天である、前記コロニーが成長した培地とを含む試料に対して直接的に行われ、
- 獲得した光度がグラム陰性タイプかつ発酵性の
細菌株の光度ではない場合、グラムタイプおよび発酵特性の判定が、獲得した光度がグラム陽性タイプの
細菌株の光度であるかどうかを予測する第2の分類の適用をさらに含み、
第1の分類が、波長範囲470nm~500nm
のスペクトルおよび波長範囲575nm~645nm
のスペクトルだけに応じて光度を区別し、
第2の分類が、波長範囲415nm~500nm
のスペクトルおよび波長範囲535nm~675nm
のスペクトルだけに応じて光度を区別する、
方法。
【請求項3】
獲得した光度がグラム陽性タイプの
細菌株の光度ではない場合、グラムタイプおよび発酵特性の判定が、獲得した光度が酵母の光度であるかどうかを予測する分類の適用をさらに含む、請求項2に記載の方法。
【請求項4】
- 獲得した光度がグラム陰性タイプかつ発酵性の
細菌株の光度であるかどうかを予測する分類が、グラム陰性タイプかつ発酵性の
細菌株の光度を、陰性タイプかつ非発酵性の
細菌株、グラム陽性タイプの
細菌株、および酵母から作り上げられたセットの光度と区別する分類であり、
- 獲得した光度がグラム陽性タイプの
細菌株の光度であるかどうかを予測する分類が、グラム陽性タイプの
細菌株の光度を、陰性タイプかつ非発酵性の
細菌株、および酵母から作り上げられたセットの光度と区別する分類であり、
- 獲得した光度が酵母の光度であるかどうかを予測する分類が、酵母の光度を、グラム陰性タイプかつ非発酵性の
細菌株から作り上げられたセットの光度と区別する分類である、
請求項3に記載の方法。
【請求項5】
獲得した光度が酵母の光度ではない場合、グラムタイプおよび発酵特性の判定が、獲得した光度がグラム陽性タイプの
細菌株の光度であるかどうかを予測する分類の適用をさらに含む、請求項1から4のいずれか一項に記載の方法。
【請求項6】
獲得した光度がグラム陽性タイプの
細菌株の光度ではない場合、グラムタイプおよび発酵特性の判定が、光度がグラム陰性タイプかつ発酵性の
細菌株の光度であるか、グラム陰性タイプかつ非発酵性の
細菌株の光度であるかどうかを予測する分類の適用をさらに含む、請求項5に記載の方法。
【請求項7】
- 獲得した光度が酵母の光度であるかどうかを予測する分類が、酵母の光度を、グラム陽性タイプの
細菌株、グラム陰性タイプかつ非発酵性の
細菌株、およびグラム陰性タイプかつ発酵性の
細菌株から作り上げられたセットの光度と区別する分類であり、
- 光度がグラム陽性タイプの
細菌株の光度であるかどうかを予測する分類が、グラム陽性タイプの
細菌株の光度を、グラム陰性タイプかつ非発酵性の
細菌株、およびグラム陰性タイプかつ発酵性の
細菌株から作り上げられたセットの光度と区別する分類であり、
- 光度がグラム陰性タイプかつ発酵性の
細菌株の光度であるか、グラム陰性タイプかつ非発酵性の
細菌株の光度であるかについての分類が、グラム陰性タイプかつ非発酵性の
細菌株の光度を、グラム陰性タイプかつ発酵性の
細菌株で作り上げられたグループの光度と区別する分類である、
請求項6に記載の方法。
【請求項8】
- 照明および獲得が、
細菌株のコロニーと、tryptone-soy寒天である、前記コロニーが成長した培地とを含む試料に対して直接的に行われ、
- 獲得した光度がグラム陽性タイプの
細菌株の光度ではない場合、グラムタイプおよび発酵特性の判定が、獲得した光度が酵母の光度であるかどうかを予測する分類の適用をさらに含む、
請求項1に記載の方法。
【請求項9】
獲得した光度が酵母の光度ではない場合、グラムタイプおよび発酵特性の判定が、光度がグラム陰性タイプかつ発酵性の
細菌株の光度であるかどうか、グラム陰性タイプかつ非発酵性の
細菌株の光度であるかどうかを予測する分類の適用をさらに含む、請求項8に記載の方法。
【請求項10】
- 光度がグラム陽性タイプの
細菌株の光度であるかどうかを予測する分類が、グラム陽性タイプの
細菌株の光度を、グラム陰性タイプかつ非発酵性の
細菌株、グラム陰性タイプかつ発酵性の
細菌株、および酵母から作り上げられたセットの光度と区別する分類であり、
- 獲得した光度が酵母の光度であるかどうかを予測する分類が、酵母の光度を、陰性タイプかつ非発酵性の
細菌株、および陰性タイプかつ発酵性の
細菌株から作り上げられたセットの強度と区別する分類であり、
- 光度がグラム陰性タイプかつ発酵性の
細菌株の光度であるか、グラム陰性タイプかつ非発酵性の
細菌株の光度であるかについての分類が、グラム陰性タイプかつ非発酵性の
細菌株の光度を、グラム陰性タイプかつ発酵性の
細菌株から作り上げられたグループの光度と区別する分類である、
請求項9に記載の方法。
【請求項11】
各分類が、波長範囲390nm~900nmのハイパースペクトル画像に対して、および、所定の精度の閾値、または所定の最大数のチャネルが取得されるまで、分類に使用されるスペクトルチャネルのセットを段階的に増加させることからなる手法に従って、訓練される、請求項1から10のいずれか一項に記載の方法。
【請求項12】
光度が、24個以下の複数のスペクトルチャネルに対して獲得される、請求項1から11のいずれか一項に記載の方法。
【請求項13】
光度が、第1および第2の分類のそれぞれに対して5以下、および好ましくは4以下の、複数のスペクトルチャネルに対して獲得される、請求項2から4のいずれか一項に記載の方法。
【請求項14】
光度の獲得が、株の細菌のコロニーのハイパースペクトルまたはマルチスペクトル画像の獲得を含み、光度が、コロニーに対応する前記画像の少なくとも1つのピクセルに応じて判定される、請求項1から13のいずれか一項に記載の方法。
【請求項15】
グラムタイプおよび発酵特性の判定が、コロニーのピクセルのセットの各ピクセルの光度に応じて行われ、グラムタイプおよび発酵特性が、前記第1の検出の結果の多数決によって判定される、請求項14に記載の方法。
【請求項16】
多数決が、ピクセルの70%以上の投票、または単純な多数決である、請求項15に記載の方法。
【請求項17】
微生物の特徴づけのためのシステムであって、
- 波長範囲390nm~900nmで微生物を照明するために構成された照明と、
- 前記照明した微生物によって反射した、または前記照明した微生物を通じて透過した光度を範囲390nm~900nmで獲得するために構成されたセンサと、
- 前記範囲内で獲得した光度に依存して、酵母または
細菌株であるものとして微生物を判定するために構成されたコンピュータユニットと
を備え、
請求項1から16のいずれか一項に記載の方法を実行するために構成された、
システム。
【請求項18】
微生物のコロニーを含む試料および前記コロニーが成長した培地を照明して、前記試料及び培地の画像を獲得するために構成された、請求項
17に記載のシステム。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、微生物分析の分野に関し、詳細には、微生物の特徴づけ、特に、酵母と細菌の識別、ならびに、後者の背景の中で、これらのグラムタイプおよびこれらの発酵または非発酵特性の識別に関する。
【0002】
都合のよいことに、本発明は、非発色、非蛍光発生、無着色の培地の中で成長した細菌のコロニーまたは酵母のハイパースペクトルまたはマルチスペクトル画像の分析に適用される。
【背景技術】
【0003】
病原微生物の分野では、微生物の特徴づけは、この微生物に感染した患者に対する処置を決定するために、微生物の種、および抗菌剤に対する微生物の感受性(または「アンチバイオグラム」)を識別することからなるのが好ましい。このために、複雑な微生物学的方法が、実験室において通常用いられ、前記方法は、ほとんどの場合、微生物の他の性質、特に、微生物の界(例えば酵母または細菌)、および、細菌の背景では、細菌のグラムタイプまたは細菌の発酵または非発酵特性、の予備知識を必要とする。実際、この情報は、特に、最終的に、微生物の種または微生物のアンチバイオグラムを判定するために、微生物に適した培地、または抗菌剤の種類を選べるようにする。例えば、出願人によって市販されている微生物API(登録商標)の識別のためのギャラリーを選ぶことは、微生物の界(例えば酵母対細菌)、または、識別されることになる菌株のグラムタイプの知識に基づく。同様に、出願人によって市販されているVitek(登録商標)2システムによって菌株のアンチバイオグラムを判定することは、前記株のグラムタイプおよび発酵または非発酵特性に応じてカードを選ぶことに基づく。また、識別されることになる微生物が酵母であるか細菌であるかに応じて異なるマトリックスを使用するMALDI-TOF質量分析による識別に言及してもよい。したがって、できるだけ早くこの情報を知ることが、特に識別を加速すること、または使用される消耗品の数を減らすことによって、微生物学的方法を最適化できるようにする。
【0004】
また、これらの性質の知識は、菌株の識別時に偽陽性を減らせるようにする。例として、出願人によって市販されているChromID(登録商標)Elite Mediumの背景では、試験される菌株の発酵特性の知識は、サルモネラ属の識別を強化する。具体的には、サルモネラ属である発酵細菌と、シュードモナス属である非発酵細菌は両方、発色性基質を変化させる。細菌が非発酵性であることが知られている場合、サルモネラ属は、追加の微生物学的アッセイがなくても、単純に排除されることが可能である。
【0005】
実験室における微生物学的方法を手引きすることを目的とした微生物の特徴づけに加えて、この情報は、臨床的にも有用である。特に、菌株のグラム分類は、例えば、菌のペプチドグリカンの割合といった、菌の細胞壁の特徴を明らかにできるようにし、細菌の分類において、または、第1の近似に対して、抗生物質に対する細菌の感受性を評価するために、使用される。したがって、細菌の2つの種類、すなわち、グラム「陽性」細菌とグラム「陰性」細菌が区分される。同様に、非発酵細菌、すなわちグルコースを異化できない細菌は、病原性細菌の中で特殊な位置を占めるということが観察される。実際、非発酵細菌は、抗生物質に対する高レベルの自然耐性があり、多くの院内感染に関与している。例として、緑膿菌(Pseudomonas aeruginosa)およびアシネトバクター属に言及してもよい。したがって、細菌の発酵または非発酵特性を素早く明らかにすることは、第1線の抗生物質治療をさらに効果的に指図できるようにし、多剤耐性株の拡散の速度を遅くする。
【0006】
歴史的に、上述の性質(界、グラム、および発酵性)のそれぞれは、専用の技法によって取得されてきた。例えば、菌株のグラムタイプは、「グラム染色」と呼ばれる手動的技法で判定されたが、この技法は、多くの手動工程(固定、染色、媒染、洗浄、過度染色、等)を含み、したがって、時間が長くかかる。したがって、細菌のグラムタイプの検出を自動化するため、特に、多くの試料を扱うために、様々な技法が開発されてきた。しかし、これらの技法は本質的に、細菌のグラムを簡単に観察できるようにするために、細菌の、または細菌の環境の、電磁応答を変化させ続ける。具体的には、第1の種類の技法は、顕微鏡のスライド上で細菌の細胞膜の染色を自動化することからなるが、グラムタイプについての最終判定のための工程は、常に、顕微鏡の下でスライドを観察する技術者によって行われる。したがって、この種類の技法は、完全に自動化されるわけではなく、その上、自動化しにくい。実際、グラム陽性細菌とグラム陰性細菌との間の色の差は微妙である可能性があり、このことが、実験室の技術者の介入が依然として必要である理由を説明する。第2の種類の技法は、細菌の細胞膜のペプチドグリカンによって開始された酵素反応によって分解される基質と、細菌を接触させることからなる。この反応は、濃縮がグラムの指標である発色団または蛍光物質を生み出す。これは、通常、細菌の発色性または蛍光発生性「ラベリング」と呼ばれる。先行技術からのこの種類の技法は、例えば、適切なデバイス(例えば分光計/蛍光計)を使用して発色団/蛍光物質の光度を測定し、次に、測定した強度を所定の閾値とコンピュータで比較することによって自動化されることが可能であるが、それでも、しばしば高価な、特殊な発色団または蛍光発生性の基質の設計を必要とする。さらに、どのような技法が使用されても、細菌は、その自然状態を変化させられ(例えば、細菌が染料を含む、発色性または蛍光性マーカでラベルを付けられる、等)、したがって、特徴づけのためのその後の試験(例えば、アンチバイオグラムの判定)のために、もはや使用されることが不可能である。
【0007】
細菌の発酵または非発酵特性の判定について、判定は、通常、試験される菌株の発酵または非発酵特性に応じて色を変化させる発色性媒質を用いる。例えば、「Kligler-Hajna」アッセイは、培地上で株を成長させることからなり、pH、ラクトース、グルコース、チオ硫酸塩、および第一鉄イオンに応じて色を変化させる比色指標を含む。この媒質は、グルコースの異化から細菌の発酵特性を検出し、pH指標の比色の色変化に反映される。また、菌株のトリブチリンエステラーゼの活動のためのアッセイ媒質に言及してもよく、非発酵性グラム陰性細菌の特徴を明らかにすることができる。
【発明の概要】
【0008】
本発明の目的は、微生物の特徴を明らかにするための方法を提案することであり、この方法は、自動的なものであり、微生物または微生物の培地に、これらの特性を判定するためにラベルを付ける、または染色する必要がない。
【0009】
このために、本発明は、微生物の特徴を明らかにするための方法に関し、
- 波長範囲390nm~900nmに自然の電磁応答を有する微生物を前記範囲で照明することと、
- 前記照明した微生物によって反射した、または前記照明した微生物を通じて透過した光度を前記範囲内で獲得することと、
- 前記範囲内で獲得した光度に応じて、酵母または菌株であるものとして微生物を判定することと、
を含む。
【0010】
「自然の電磁応答」は、本発明の意味では、関心のある少なくとも波長範囲での照明への微生物の電磁応答を変化させる要素(染料、色素原、フルオロゲン、等)によって、微生物である酵母または菌株が変化しないということを意味する。例えば、微生物のコロニーは、非発色かつ非蛍光培地の中で培養され、照明/獲得は、微生物の媒質の中に依然として存在するコロニーに対して直接的に行われる。
【0011】
言い換えれば、波長範囲390nm~900nmで、酵母および細菌を区別できる電磁特性を酵母および細菌が「自然に」有するということを、発明者は発見した。したがって、発色性または蛍光発生性の基質または染料を使用する必要はない。さらに、本発明による方法は、照明すること、スペクトルを測定すること、および、このスペクトルの処理、特にコンピュータ処理を実行すること、からなる限りにおいて速い。酵母または細菌レベルで微生物のこの特徴づけは、例えば、上述のような実験室の微生物学的方法を最適化できるようにする。
【0012】
酵母または細菌であるものとして微生物の特徴づけは、獲得した光度から直接的に行われ、微生物の種を事前に判定する必要がないということに留意されたい。特に、種のレベルで識別する必要がないことは、種のレベルが2つのクラスに限定されることが可能なので、予測モデルを大いに単純化するという長所をもたらす。
【0013】
都合のよいことに、方法は、微生物のコロニーが成長した栄養寒天培地を含むペトリ皿に応用される。例えば、栄養培地は、酵母または細菌を含む、または含むと疑われる生物学的試料、例えば尿、を使用して植え付けられ、次に、コロニーを成長させるために培養される。コロニーが培地で検出されるとすぐ、本発明の方法によって特徴が明らかにされる。したがって、方法は、培地の植え付け後に、いかなる物質移動または試薬追加も必要としない。例えば、コロニーは、一定間隔でペトリ皿の画像を撮影し、コロニー検出アルゴリズムを用いることによって、自動的に検出される。
【0014】
都合のよいことに、本発明による方法は、微生物の自己蛍光の分析ではなく、微生物の反射または微生物の吸光度の分析に基づく。具体的には、一般に照明は強すぎて、ハイパースペクトル画像またはマルチスペクトル画像で自己蛍光を観察することができない。
【0015】
1つの実施形態によれば、方法は、菌株のグラムタイプおよび発酵特性の検出をさらに含み、
- 波長範囲390nm~900nmに自然の電磁応答を有する前記株の少なくとも1つの細菌を前記範囲で照明することと、
- 前記照明した細菌によって反射した、または前記照明した細菌を通じて透過した光度を前記範囲内で獲得することと、
- 前記範囲内で獲得した光度に応じて、菌株のグラムタイプおよび発酵特性を判定することと
を含む。
【0016】
言い換えれば、細菌は、細菌のグラムタイプおよび細菌の発酵または非発酵特性に特徴的な電磁特性を「自然に」有する。したがって、本発明による方法は、この特性を測定し、次に、ここから、細菌のグラムタイプおよび発酵特性を抽出することからなる。具体的には、本発明により、菌株がグラム陽性であるか、グラム陰性かつ発酵性であるか、グラム陰性かつ非発酵性であるかを、範囲390nm~900nmにより、判定することができ、この情報を知っていることが、例えば、上述のような実験室の微生物学的方法を最適化できるようにする。
【0017】
また、本発明は、本発明による方法を実行するためのシステムを較正するための方法に関し、システムは、
- 波長範囲390nm~900nmで微生物を照明するために構成された照明と、
- 前記照明した微生物によって反射した、または前記照明した微生物を通じて透過した光度を範囲390nm~900nmで獲得するために構成されたセンサと、
- センサによって獲得した強度を分析するための命令を収めることができるコンピュータメモリ、および、コンピュータメモリに収められた分析命令を実行することができるマイクロプロセッサ、を備えるコンピュータユニットと
を含み、
較正の方法は、
- 範囲390nm~900nmで照明された細菌および酵母の前記範囲における光度を含む学習データベースを構築することと、
- 前記データベースに応じて、酵母または細菌の予測モデルの自動学習をコンピュータによって実行することと、
- 学習した予測モデルを実行するための分析命令をシステムのコンピュータメモリに格納することと
という工程を含む。
【0018】
また、本発明は、治療方法に関し、
- 酵母または細菌による感染症があると疑われる患者から試料を採取することと、
- 都合のよいことに、寒天培地に試料を植え付けること、前記植え付けた媒質を培養して微生物のコロニーを成長させること、および成長した1つまたは複数のコロニーを検出することによって、試料内にある1つまたは複数の微生物を検出することと、
- 前記範囲内で獲得した光度に応じて、酵母または菌株であるものとして微生物を判定することと、
- 前記判定の結果に応じて、1つまたは複数の抗菌物質を選択することと、
- 選択した1つまたは複数の抗菌物質を患者に投与することと
を含む。
【図面の簡単な説明】
【0019】
本発明は、純粋に例として示された以下の説明を読み、同一の参照番号が同一または類似の要素を示す添付の図面を参照すると、よりよく理解されるであろう。
【0020】
【
図1】本発明によるハイパースペクトルシステムの概略図である。
【
図2】本発明によるマルチスペクトルシステムの概略図である。
【
図3】
図2のシステムで使用される帯域通過フィルタの透過スペクトルの例を示す図である。
【
図4】
図1または
図2のシステムを使用して適用されるクラスY、GP、GNF、GNNを予測するための方法の流れ図である。
【
図5】「ステップフォワード」法によって特異なスペクトルチャネルを選択するための方法の流れ図である。
【
図6】クラスY、GP、GNF、GNNを予測するためのモデルを学習するための方法の流れ図である。
【
図7A】クラスY、GP、GNF、GNNについてのフラットな予測モデルを示す図である。
【
図7B】クラスY、GP、GNF、GNNの系統樹による階層式予測モデルを示す図である。
【
図7C】クラスY、GP、GNF、GNNの最適木による階層式予測モデルを示す図である。
【
図8】クラスY、GP、GNF、GNNの予測モデル学習のために使用される細菌種および酵母種を記述するテーブルである。
【
図9】予測モデルの較正および交差検証のために使用されるピクセル数、およびしたがってスペクトル数を記述する、COS媒質についてのテーブルである。
【
図10】予測モデルの較正および交差検証のために使用されるピクセル数、およびしたがってスペクトル数を記述する、TSA媒質についてのテーブルである。
【
図11】重みを付けられた予測レート、または「バランスのよい分類レート」の計算を示す図である。
【
図12】COS媒質の最適木のための主な特異なスペクトルチャネルの空間分散を示す図である。
【
図13】COS媒質の最適木のための各モデルの5つの主な特異なスペクトルチャネルの空間分散を示す図である。
【
図14】COS媒質の最適木の各モデルの5つの主な特異なチャネルを個別に示す図である。
【
図15】TSA媒質の最適木のための主な特異なスペクトルチャネルの空間分散を示す図である。
【
図16】COS媒質の最適木のための主な特異なスペクトルチャネルの空間分散を示す図である。
【
図17】TSA媒質の最適木の第1、第2、および第3の予測モデルとそれぞれ関連付けられた12個、8個、および11個の主な特異なチャネルをそれぞれ示す図である。
【
図18】系統樹の第1の予測モデルを示す図であり、図の一番上にある図がTSA媒質に対応し、図の一番下にある図がCOS媒質に対応する。
【
図19】系統樹の第2の予測モデルを示す図であり、図の一番上にある図がTSA媒質に対応し、図の一番下にある図がCOS媒質に対応する。
【
図20】系統樹の第3の予測モデルを示す図であり、図の一番上にある図がTSA媒質に対応し、図の一番下にある図がCOS媒質に対応する。
【発明を実施するための形態】
【0021】
以下に、表記法Ai,jは、マトリックスAのi行j列の要素に関する。
【0022】
図1を参照すると、ペトリ皿に注がれた寒天上で成長した酵母または細菌のコロニーの特徴づけのためのハイパースペクトルシステム10は、
- ハイパースペクトル画像獲得のためのデバイス12と、
- デバイス12を制御するため、および、デバイス12によって獲得された画像を受信して処理するために、(例えば、接続されたリンクまたはワイヤレスによって)デバイス12に接続されたコンピュータデータ処理ユニット14と
を備える。
【0023】
例えば、米国モンタナ州の会社Resononからの、リファレンス「Pika II」のハイパースペクトルイメージングシステムといった、デバイス12は、以下を含む。
- 波長範囲[λmin;λmax]=[390;900]ナノメートルで感度がよい、例えば、CCDまたはCMOSタイプのデジタルセンサといった、エレメンタリセンサのアレイを備えるデジタルセンサ、および、センサによって獲得されることになる波長を選択するための光散乱要素またはスペクトルグラフからなる、いわゆる「ハイパースペクトル」カメラ18。
- ハイパースペクトル画像が獲得されることになるペトリ皿22の光学画像の焦点を、カメラ18のデジタルセンサに合わせるための対物レンズ20。
- 範囲[λmin;λmax]で光を発することができ、ペトリ皿22の一様な前面照明を提供するための、例えば2つまたは4つのランプといった、例えば1つまたは複数のハロゲンランプからなる前面照明24。例えば、照明は、白色光のランプを含む。
- 範囲[λmin;λmax]でペトリ皿22の一様な背面照明を提供するための、例えば白色光LEDのマトリックスからなる背面照明26。
- スキャンすることによって皿22の完全な画像を取得するために、ペトリ皿22を保持し、対物レンズ20の前をペトリ皿22が通ることを可能にするキャリッジ28。
【0024】
このようにして、照明は、範囲[λmin;λmax]全体で提供される。
【0025】
例えば、デバイス12は、160マイクロメートルのサンプリングステップ(300マイクロメートルで推定される空間分解能)で、および、範囲[λmin;λmax]にわたる1.7ナノメートルのスペクトル分解能で、縦90ミリメートル横90ミリメートルの領域の画像を獲得するために構成される。
【0026】
したがって、デバイス12は、ペトリ皿によって反射された光の、N行M列のデジタル画像HSIを生み出し、ペトリ皿22は開いている(すなわち、ペトリ皿22のカバーがない)のが好ましい。
ピクセルの放射輝度は、一般に「光度」と呼ばれ、この場合、例えばデジタル写真撮影の分野から本質的に知られているように、露出時間中にカメラのセンサ18の対応する基本的な高感度部位の表面に入射する光の量に対応する。
【0027】
各ピクセルRad
i,j(λ)は、異なる波長[λ
min;λ
max]におけるピクセルに対応する皿22の放射輝度のデジタルスペクトルからなり、デジタルスペクトルは、以下の関係によって表現される。
ここで、Δλは、スペクトル分解能であり、pは、
に属する正の整数である。獲得波長λ
min+p×Δλは、通常、用語「チャネル」によって示される。
【0028】
データ処理ユニット14は、例えば、パーソナルコンピュータ、タブレット、スマートフォン、サーバ、スーパーコンピュータ、または、より一般には、獲得デバイス12によって生み出された画像HSIを処理するために構成された、特にDSP(「デジタルシグナルプロセッサ」)タイプのマイクロプロセッサに基づく、FPGAタイプの回路に基づく、これらのタイプの技術をミックスした回路に基づく、等のいずれかのシステムである。特に、ユニット14は、デバイス12によって生み出された画像、本発明による方法の実行のためのコンピュータ命令、前記実行に役立つパラメータを格納するため、ならびに、中間のおよび最終的な計算の結果を格納するための、メモリのセット(RAM、ROM、キャッシュ、大容量ストレージ、等)を提供される。ユニット14は、任意選択として、コロニーの特徴づけの最終的な結果、具体的には、グラムタイプおよび/もしくは発酵特性、ならびに/または、調査されるコロニーの細菌もしくは酵母の特徴の判定、を表示するためのディスプレイ画面を備える。単一の処理ユニットが説明されるが、本発明は、当然、いくつかの処理ユニット(例えば、画像HSIの事前処理を行うためのカメラ18の内部にマウントされたユニット、および、他の処理を実行するためのデバイス12の外部のユニット)によって行われる処理に適用される。さらに、システムは、例えば、オペレータが利用できるキーボード/マウスおよびドロップダウンメニュー、ペトリ皿の上にあり媒質についての情報を含むバーコード/QRコードを読み取るバーコード/QRコードリーダ、等によって、予測が媒質によって決まる場合に使用される試料に関するデータ、特に培地の種類をユニット14に入力できるようにするインターフェースを追加されてもよい。
【0029】
図1のハイパースペクトルシステムは、コロニーのクラスを予測するための種々のモデルに適合させることができ、多くのスペクトルチャネルを使用して予測精度を向上させることができるので、獲得波長に関して軽快であるという長所がある。高価格であることに加えて、このようなシステムは一般に、センサに入射する光の強度で画像を単に獲得することが目的の従来のCMOSまたはCCDカメラより空間分解能が小さい。
【0030】
図2を参照すると、マルチスペクトルシステム32は、例えば、対物レンズ20とカメラのセンサ32との間の、対物レンズ20の前に配列されたスペクトルフィルタ36のセットに連結された、高い空間分解能のCMOSまたはCCDカメラであることが好都合なカメラ32の分だけ、ハイパースペクトルシステム10とは異なる。フィルタ36のセットは、N
F個の別個の帯域通過フィルタからなり、それぞれは、範囲[λ
min;λ
max]のうちの範囲[λ
1;λ
2]の光だけを透過させるために構成され、スペクトルの半値全幅(FWHM)は50nm以下であり、好ましくは20nm以下である。例えば、420nmに中心を置く、会社Edmund Opticsからのフィルタといった、この種類のフィルタの透過スペクトルが、
図3に示される。セット36は、例えば、24個までの異なるフィルタを受けることができるフィルタホイールであり、ホイールはユニット14によって制御され、カメラの前に沿って前記フィルタを移動させ、フィルタのそれぞれに対する画像の撮影を制御するように作動する。
【0031】
マルチスペクトル画像HSI(λ)は、このようにして獲得され、画像の各ピクセルRad
i,j(λ)は、セット36によって取り除かれた種々のスペクトル帯域のピクセルに対応する皿22の放射輝度のデジタルスペクトルからなり、デジタルスペクトルは、以下の関係によって表現される。
ここで、λ
1,λ
2、…、
は、それぞれ、セット36のスペクトルフィルタの中心波長である。
【0032】
たった今説明されたシステムによる生物学的試料(例えば、尿、血液、気管支肺胞の試料、等)に含まれる微生物の特徴づけの方法40が、次に、
図4の流れ図に関して詳細に提示される。具体的には、この方法は、都合のよいことに、MALDI-TOF質量分析(例えば、出願人によって市販されているVitek(登録商標)MS)を使用して、前記試料に含まれる微生物を、および、(例えば、出願人によって市販されているVitek(登録商標)2プラットフォームを使用して)前記微生物のアンチバイオグラムを、識別するための、より一般的な方法の背景における応用を見つける。本質的に知られているように、これらの技法のそれぞれは、微生物の種類に関する媒質および/または試薬および/または特定の消耗品を選ぶ必要がある。例として、MALDI-TOF質量分析による酵母の識別は、都合のよいことに、この種類の技術において使用されるマトリックスにおけるギ酸の使用を伴う。さらに、Vitek(登録商標)2プラットフォームのカードの選択(カードは、アンチバイオグラムにおいて試験される成長培地および1つまたは複数の抗菌物質を含む)は、試験される微生物の細菌の特徴によって決まる。特に、発酵性グラム陰性細菌は、細菌のアンチバイオグラムを生成するための特定のカードを必要とする。
【0033】
都合のよいことに、以下で説明される特徴づけの方法40は、ペトリ皿上での最初の成長からスタートして、他の微生物学的方法のために微生物についての必要情報を取得すること、具体的には、成長したコロニーが酵母(「Y」)に相当するか細菌に相当するか、また、細菌の場合、この細菌がグラム陽性タイプ(「GP」)であるかグラム陰性タイプ(「GN」)であるか、また、グラム陰性タイプの細菌の背景では、この細菌が発酵性(「GNF」)であるか非発酵性(「GNN」)であるか、を知ることができるようにする。したがって、方法は、微生物のクラス、すなわち、クラスY、GP、GNF、またはGNNを予測できるようにする。
【0034】
方法の第1の工程42において、ペトリ皿は、ペトリ皿に堆積した栄養培地、または「培」地の表面で酵母または細菌のコロニーを成長させるように、例えば患者から採取された生物学的試料を植え付けられる。培地の主な目的は、前記コロニーを成長させること、および任意選択として、光の摂動を限定することによって特徴づけの精度を強化することである。反射した光度に応じたグラムタイプの検出について、培地は不透明であることが好ましく、検出の精度のレベルを向上させる。特に、不透明な媒質は、反射係数ρが10%以下であり、好ましくは5%以下であり、さらにいっそう好ましくは1%以下である。例えば、いわゆる「CPSO」寒天培地(「CPS」寒天は、媒質を不透明にするためにSiO2を含む)、「コロンビア」寒天(または「CNA」寒天)、ヒツジ血液の5%を含むコロンビア寒天(または「COS」寒天)、Man、Rogosa、Sharpe寒天(「MRSM」寒天)、チョコレート寒天(「PVX」寒天)、Tryptone-Soy寒天(「TSA」寒天)、等である。
【0035】
コロニーのこの種類の成長は従来のものであるので、より詳細に以下で説明されることはない。コロニーの成長は、都合のよいことに、オペレータによって手動で、または、本質的に知られた様式の自動的植え付け機を使用して自動的に、行われてもよい。都合のよいことに、調製は、コロニーに基づいて微生物の特徴づけが行われるコロニーが、間隔を空けられるように、および、デバイス12によって獲得した画像の複数のピクセルにコロニーのエリアが対応するように、行われる。これは、獲得した画像のピクセルのその後の識別、およびしたがって、画像処理アルゴリズムによるピクセルのセグメント化、または、ユーザによる、画像からのピクセルの抽出、を特に容易にすることができる。
【0036】
コロニーの成長が終わると、例えば、24時間、36時間、または48時間後、ペトリ皿は開けられ、キャリッジ28上に配列されるのが好ましく、照明24および26が点灯され、ペトリ皿の少なくとも1つのハイパースペクトル(または、マルチスペクトル)画像HSIが、44において、獲得デバイス12(または32)を使用して獲得され、コロニーを作り上げる微生物の種類を、獲得した画像から判定するためのコンピュータ処理を使用する処理ユニット14に格納される。
【0037】
ユニット14は、任意選択として、46において、以下の処理動作のうちの1つ、または、これらの処理動作の任意の組合せからなるノイズの事前処理を始める。
a. 本質的に知られた様式での、カメラのセンサのノイズ、特に、そのオフセット、その空間ノイズ、等の補正。
b. 画像HSIに「ハイライト」を形成する特に鏡のような、付随的に発生する好ましくない反射の処理。例えば、例えばピクセルがとることができる最大値の3分の2以上(すなわち、0から255の8ビットでエンコードされたピクセルの場合、170以上)といった、所定の閾値を超える値をとるピクセルを除去するために使用される閾値化。
c. 細菌の種類、および、使用される寒天の種類で変化しない波長で反射した光度ごとに画像HSIを分割することによって、照明の変化などの外部変動によって引き起こされる画像の変化を弱めることができるようにするレシオメトリック処理。
d. いくつかの画像HSIが獲得された場合、異常なピクセル値および/または獲得した画像の平均を見つけて削除すること。
【0038】
都合のよいことに、処理は、48において、ペトリ皿単独で生成された信号を抽出するために、ハイパースペクトルまたはマルチスペクトル反射画像における種々の波長の放射輝度の値を格納する事前処理された画像HSIの変換を続ける。これは、特に、照明24、26の光源の発光スペクトルの変動を取り除けるようにする。例えば、「フラットフィールド補正」(FFC)という種類の補正が、反射を取得するために使用され、これには、さらに、ピクセルごとのセンサの反応のばらつき(暗電流の分散、利得の分散、等)を補正するという長所がある。
ハイパースペクトル画像の背景において、この変換は、例えば、以下の関係による補正である。
ここで、Υ(λ)は、反射画像であり、Wは、例えば、90%を超える一様な反射率のシート(例えば、いわゆる「白い」シート、または、10%より小さい灰色のカードを伴う)といった、照明24、26によって照明された高反射率の中性オブジェクトの、ユニット14に格納されたハイパースペクトル画像であり、Bは、例えば対物レンズ20を覆う黒いキャップの画像といった低反射率の、また、m(λ
min+p×Δλ)=1であるか、マトリックスW(λ
min+p×Δλ)-B(λ
min+p×Δλ)の平均に等しい、中性オブジェクトの、ユニット14に格納されたハイパースペクトル画像である。
【0039】
同様に、マルチスペクトル画像の背景において、変換は、例えば、以下の関係による補正である。
ここで、Wは、例えば、90%を超える一様な反射率のシート(例えば、いわゆる「白い」シート、または、10%より小さい灰色のカードを伴う)といった、照明24、26によって照明された高反射率の中性オブジェクトの、ユニット14に格納されたマルチスペクトル画像であり、Bは、例えば対物レンズ20を覆う黒いキャップの画像といった低反射率の、また、m(λ
n)=1であるか、または、マトリックスW(λ
n)-B(λ
n)の平均に等しい、中性オブジェクトの、ユニット14に格納されたマルチスペクトル画像である。
【0040】
ユニット14は、50において、例えば画像HSI(λ)またはΥ(λ)から、細菌のコロニーを識別するためのアルゴリズムを、工程38に続いて、または、前述のステップと並行に、実行する。任意の従来の形状およびオブジェクト認識アルゴリズムが、コロニーに対応する、「Col(λ)」と呼ばれる画像のゾーンを抽出するために使用されてもよい。変形形態として、この選択はオペレータによって手動で行われ、オペレータは、ディスプレイ画面、および例えばマウスタイプのポインティングデバイスによってこのゾーンを選択する。例として、ゾーンCol(λ)は、コロニーに属するピクセルの座標のリストからなる。ピクセルの選択したゾーンは、ユニット14によって格納される。
【0041】
方法は、52において、所定の判定ルールを適用することによって、ゾーンCol(λ)の少なくとも1つのスペクトルに応じて、コロニーの微生物のクラスY、GP、GNF、またはGNNの予測を続けるが、この変形形態が以下で説明される。具体的には、この予測は、ゾーンCol(λ)の各ピクセル(i,j)のスペクトルΥi,j(λ)に応じて行われる。このために、クラスの最初の予測は、ゾーンCol(λ)の各ピクセル(i,j)に対して行われ、次に、クラスの最終的な予測のために多数決が使用される。第1の変形形態では、単純な多数決が使用され、すなわち、ゾーンCol(λ)の最大数のピクセルに基づいて予測されたクラスが、最後に保持されるクラスである。第2の変形形態では、クラスの予測の確実性を向上させるために、条件付きの投票が使用され、すなわち、最後に保持されるクラスは、ゾーンCol(λ)を作り上げるピクセル数のX%より大きいことに基づいて予測されるものであり、Xは、完全に50%より大きく、70%以上であることが好ましい。クラスがこの条件を満たさない場合、方法は、クラスの予測なしを返す。当然、コロニーのクラスは、例えば、ゾーンCol(λ)のスペクトルのセット{Υi,j(λ)}(i,j)∈Col(λ)の平均スペクトルΥcol(λ)といった、単一の値によって提供されてもよい。
【0042】
各コロニーのクラスY、GP、GNF、またはGNNは、所定の予測ルールを適用することによってユニット14によって見つけられるが、この変形形態が以下で説明される。予測されたクラスはユニット14に格納され、および/または、ユーザの注目のために画面に表示される。また、この予測は、都合のよいことに、識別のその後の工程、および/または、コロニーを形成した微生物のアンチバイオグラムのために、別の微生物分析器具に送られる。
【0043】
コロニーのピクセルのスペクトルΥ
i,j(λ)に応じてクラスY、GP、GNF、またはGNNを予測するための様々なモデル、特に、教師あり自動学習(「SML」:教師あり機械学習(supervised machine learning))に基づく予測モデルが次に説明される。使用されるSMLツールは、最初に、
図5を参照しながら説明され、次に、
図6および
図7に示された予測モデルの学習処理を通じて、予測モデルが下記で説明される。
【0044】
A.教師あり機械学習のためのツール
どのような学習が考えられるにしても、これは、学習データベースの構築で始まる。各クラスY、GP、GNF、およびGNNについて、細菌および酵母が選択され、これらのそれぞれは、ペトリ皿に注がれた寒天に植え付けられ、所定の時間、培養され、
図1で説明されたシステムによって、また必然的に、照明の同じ条件の下、および波長範囲390nm~900nmで、皿のハイパースペクトル画像が獲得される。寒天上で成長したコロニーのピクセルは、例えば方法40の工程46~50で説明された手法で抽出され、これらの関連付けられたスペクトルが学習データベースに格納される。したがって、学習データベースは、クラスY、GP、GNF、GNNそれぞれと関連付けられたスペクトル
、
、
、
の4つのセットを含む。これらのセットのそれぞれは、先行技術から本質的に知られているように、適切な学習のために使用される「較正」と呼ばれる第1の部分、および、計算した予測モデルの性能を評価するために使用される「交差検証」と呼ばれる第2の部分、という、2つに分けられる。
【0045】
第1の好ましい実施形態によれば、コンピュータによって行われる学習は、「ステップフォワード」と呼ばれ、予測モデルを訓練するために使用される。この種類の学習は、本質的に貧弱であり、
図2のシステムによって行われるような、マルチスペクトルアプリケーションを見つけるのに適しているように、最も特異なスペクトルチャネルの段階的選択に基づく。
【0046】
図5を参照すると、学習60は、初期化工程62で始まり、ここで、最大数R個の特異なチャネルが選択され、この数は、1と、スペクトルの獲得のために使用されるハイパースペクトルカメラのスペクトルチャネルの数Pとの間にある。選択された特異なチャネルのリストLは空にされ、候補チャネルのリストlは、ハイパースペクトルカメラのチャネルのセット{λ
1,λ
2,…,λ
P}に初期化される。
【0047】
後続の反復工程64において、リストlは、反復工程66を使用して、リストLからの最も特異なチャネルRで段階的に満たされる。より具体的には、工程64、工程66の所与の反復rの間、以下を行う。
- リストlから各チャネルλ
kを抽出すること。
- 68において、抽出したチャネルλ
kおよびリストLからのチャネル{λ
1,λ
2,…,λ
r-1}について、予測モデル
を判定すること。
- 70において、例えば、交差検証におけるスペクトルの妥当な分類のレベルといった、予測モデル
の性能尺度
を計算する。
【0048】
次に、工程64は、72において、最善の性能尺度を示す予測モデルの識別、および必然的に、リストLのチャネルと組み合わせたリストlの最も特異なチャネルλ
rの識別を続ける。74において、リストLは、次に、チャネルλ
rで完了し、後者は、工程64の次の反復r+1の間、リストlから削除される。R個の最も特異なチャネルが識別されると、学習方法は、74において、リストLの、および、後者、すなわち、工程72で識別された最後のモデルと関連付けられた予測超平面
の、格納で終わる。
【0049】
都合のよいことに、工程68で計算された予測モデルは、線形カーネルを伴い、ソフトマージンを伴う(「サポートベクトルマシン」のための)SVMタイプである「1対多(one against all)」のものである。この種類の学習は、較正スペクトルに応じて計算することからなり、超平面
は、下記で説明されることになるように、他のクラスの中からの1つ、2つ、または3つから形成されたセット(Cl)
のクラスCl(Cl=Y、GP、GNF、またはGNN)を区別する。例えば、モデルは、工程66の反復kおよび工程64の反復rの間、
という制約の下、
という関係に従って最適化問題を解決することによって訓練される。上記の式において、
- クラスClに、または、セット
に属する較正スペクトルΥ
m(λ)について、
は、リストL={λ
1,λ
2,…,λ
r-1}の、および、反復kの中でリストlから抽出されたチャネルλ
kの、スペクトルチャネルに対応するΥ
m(λ)の成分のベクトル、好ましくは、成分が、チャネルの値に応じて規則正しく配列されたベクトルに等しい。
-
およびβは、スペクトル
の次元に等しい次元の、および必然的に、rに等しい次元の、ベクトルである。
- Mは、クラスClに、または、1からMまで番号を付けられたセット
に、属する較正スペクトルΥ
m(λ)の数である。
-
.βは、ベクトル
とベクトルβのスカラ積である。
- ξ
mおよび
はスカラである。
- q
m∈{-1,1}であり、m番目の学習スペクトルがクラスClと関連付けられる場合、q
m=1であり、m番目のスペクトルが他のクラスのセット
と関連付けられる場合、q
m=-1である。
- Cは、所定のスカラである。
【0050】
クラスClにおけるピクセルΥ
i,j(λ)のスペクトルのメンバシップを予測するモデルは、したがって、以下の工程に従って生み出される。
- スペクトルΥ
i,j(λ)の、ベクトル
への変換。
- スペクトル
と超平面
との間の距離
の計算。
- 距離S
clに応じてクラスClを予測するためのルールの適用、例えば、スペクトルは、S
clの符号が正である場合、クラスClに属し、この符号が負である場合、セット
に属する。
【0051】
第2の実施形態によれば、学習は貧弱ではなく、チャネル全てを同時に使用することからなり、ここから生じる予測モデルは、
図1のシステムによるハイパースペクトルアプリケーションに特に適している。例えば、この学習は、線形カーネルおよびソフトマージンを伴うSVMタイプである「1対多」のものであり、
という制約の下、
という関係に従って最適化問題を解決することによって、他のクラスの中からの1つ、2つ、または3つから形成されたセット
のクラスCl(Cl=Y、GP、GNF、またはGNN)を区別する超平面
を、較正スペクトルに応じて計算することからなる。
これらの式において、
- Υ
m(λ)は、クラスClに、またはセット
に属する較正スペクトルΥ
m(λ)である
-
およびβは、較正スペクトル
の次元に等しい次元の、および必然的に、Pに等しい次元の、ベクトルである。
- Mは、クラスClに、または、1からMまで番号を付けられたセット
に、属する較正スペクトルΥ
m(λ)の数である。
- Υ
m(λ).βは、ベクトルβ倍したベクトルΥ
m(λ)のスカラ積である。
- ξ
mおよび
はスカラである。
- q
m∈{-1,1}であり、m番目の学習スペクトルがクラスClと関連付けられる場合、q
m=1であり、m番目のスペクトルが他のクラスのセット
と関連付けられる場合、q
m=-1である。
- Cは、所定のスカラである。
【0052】
クラスClへのピクセルΥ
i,j(λ)のスペクトルのメンバシップを予測するモデルは、したがって、以下の工程に従って生み出される。
- スペクトルΥ
i,j(λ)と超平面
との間の距離
の計算。
- 距離S
clに応じてクラスClを予測するためのルールの適用、例えば、スペクトルは、S
clの符号が正である場合、クラスClに属し、この符号が負である場合、セット
に属する。
【0053】
B.予測モデルを訓練する方法
図6を参照すると、予測モデルを訓練するための方法80は、82において、上述のような、クラスY、GP、GNF、およびGNNのための学習データベースの構築で始まる。方法80は、84において、予測モデルの構造の決定を続ける。具体的には、一方では
図7Aにおいて、また、その一方では
図6Bおよび
図6Cにおいて、それぞれ示されたものなど、2種類のモデルが可能である。
【0054】
第1の種類の予測モデルは、
図7Aに示され、「1対多」というタイプ90、72、74、76の4つの予測モデル、すなわち、クラスY対クラスGP、GNF、GNNの予測、クラスGP対クラスY、GNF、GNNの予測、クラスGNF対クラスY、GP、GNNの予測、およびクラスGNN対クラスY、GP、GNFの予測、を訓練することからなる。このために、
- 各予測モデルは、関係(6)または関係(7)について上述された学習ツールのうちの1つを使用して、学習データベースに応じて訓練され、クラスClは、Y、GP、GNF、GNNに等しく、セット
は、他の3つのクラスから作り上げられる。
- ピクセルのグラムタイプおよび発酵特性の予測(
図3の工程52)は、前記ピクセルから超平面
、
、
、
への、または、超平面
、
、
、
への、スペクトルΥ
i,j(λ)の距離S
Y、S
GP、S
GNF、S
GNNをそれぞれ計算し(工程90~工程96)、次に、88において、計算した距離に応じてピクセルのクラスを判定することによって取得される。特に、保持されるクラスは、最大距離に対応するものである。
【0055】
「フラット」と呼ばれる
図7Aに示された第1の種類の構造によれば、クラスY、GP、GNF、およびGNNは、等しく重要であるとみなされ、必然的に、識別は同様にエラーとなる。例えば、細菌GNNの代わりに酵母Yを識別することは、細菌GNNの代わりに細菌GNFを識別することと同様に深刻である。
図7Bに示された構造によれば、予測モデルは、系統発生学的分類木に従って編成され、その結果、異なるクラスが、等しく重要であるものとみなされる必要はもはやなく、情報を演繹的に、すなわち、スペクトルの形状に影響を及ぼす可能性のある進化情報を、導入できるようになる。より具体的には、この予測モデル木は、以下を含む。
- 酵母Yを細菌GP、GNF、GNFと区別することからなる第1のモデル100
- 細菌GPを細菌GNFおよびGNNと区別することからなる第2のモデル102
- 細菌GNFを細菌GNNと区別することからなる第3のモデル104
【0056】
モデル100~104のそれぞれは、関係(6)または関係(7)について上述された手法で取得され、クラスY、GP、GNF、およびGNNのうちの1つへのピクセルΥi,j(λ)のスペクトルのメンバシップの予測は、したがって、第1のモデル100の超平面への距離SYを計算することからなり、この距離の符号が正である場合、クラスYが予測される。そうでなければ、第2のモデル102の超平面への距離SGPが計算され、この距離の符号が正である場合、クラスGPが予測される。そうでなければ、第3のモデル104の超平面への距離SGNFが計算され、次に、クラスGNFが予測される。そうでなければ、クラスGNNが予測される。
【0057】
表現型モデルは、
図7Aに示されるような、フラットな予測構造に対する予測精度を改善できるようにするが、しかし、表現型木は、必ずしも、より妥当な結果を示す木ではないということに、発明者は気づいた。特に、特徴が明らかにされることになる微生物が成長してきた培地に依存する木が、好ましいことがあり、この媒質は、スペクトルの形状に影響を及ぼす。
図7Cに示されるような最適な予測構造が、較正スペクトルに応じて決定されるのが好ましい。より具体的には、第1の工程において、4つの予測モデルである、a)Y対GP、GNF、およびGNN、b)GP対Y、GNF、およびGNN、c)GNF対Y、GP、およびGNN、ならびに、d)GNN対Y、GP、およびGNFは上述のように計算され、最善の予測性能を有するモデルが、最適木の第1のモデル110として保持される。第2の工程において、第1のモデルのクラスは削除され、残りのクラスに対応する3つの予測モデルが計算される。例えば、第1のモデルが、クラスGNNに対応する場合、第2の工程の3つのモデルは、上述のように、a)Y対GPおよびGNF、b)GP対YおよびGNF、ならびに、c)GNF対YおよびGPである。3つのモデルの最善のものが、次に、最適木の第2のモデル112として保持される。第3の工程において、第1および第2のモデル110および112のクラスは削除され、2つの残りのクラスの間の予測モデルが上述のように計算され、最適木の第3のモデル114として保持される。クラスY、GP、GNF、およびGNNのうちの1つへの、ピクセルΥ
i,j(λ)のスペクトルのメンバシップの予測が、次に、
図7Bを参照しながら説明されたものと同様に、木を通過することによって取得される。
【0058】
図6に戻ると、予測モデルの構造が決定されると、学習処理80は、任意選択として、84において、予測のために使用されるスペクトルチャネルの数の削減を続け、この削減は、チャネルの選択および/または再グルーピングによって行われる。具体的には、
図7A、
図7B、および
図7Cを参照しながら上述された予測モデルが、
図5の「ステップフォワード」法に基づいて計算されるとき、使用されるチャネルの数は、パラメータRによって直接的に固定されてもよい。変形形態として、または追加として、モデルの性能の著しい向上を示さない追加のチャネルは、削除されてもよい。また、チャネルの再グルーピングは、変形形態または追加として、
図2および
図3について上述されたように、フィルタの幅に幅が対応する間隔に範囲390nm~900nmを分割することによって行われてもよい。単一のスペクトルチャネルが、次に、間隔ごとに保持される。λ
1,λ
2、…、
の場合のチャネルの最終的な数は、このようにして選択され、
図2のマルチスペクトルシステムのセット36のスペクトルフィルタの中心波長を定義する。下記で説明されることになるように、24個のチャネル、およびしたがって、24個のスペクトルフィルタだけを使用して高精度のクラスの予測を取得することができる。
【0059】
任意選択として、マルチスペクトルアプリケーションのための最終的なチャネルを選択し、マルチスペクトルシステムを同様に構築すると、
図2のシステムによるスペクトルの獲得に基づいて、新しい訓練動作が、予測モデルを洗練するために使用され、この訓練は、
図6および
図7について説明されたものに似ている。
【0060】
さらに、所定のR個の特異なチャネルの選択が説明されてきた。変形形態として、この数は、演繹的に固定されず、ギャップを見つけるための停止尺度は、チャネルの数に応じた性能の利得の停滞である。例えば、少なくとも1つのチャネルの追加が、X%を超えるまで、例えば下記で詳述されるBCRといった性能を向上させない場合、チャネルのサーチは停止され、例えばXは2%以下である。
【実施例】
【0061】
C.実施例
次に、たった今説明されたクラスY、GP、GNF、およびGNNの予測の応用を説明する。このために、21個の菌株および酵母株が使用され、これらの微生物種が、
図8に記述される。これらの種は、24時間、COS寒天およびTSA寒天上で培養され、これらの媒質のそれぞれに対する学習データベースを生じた。種および媒質のそれぞれについてのコロニー数およびピクセル数が、
図9(COS)および
図10(TSA)にそれぞれ記述され、ブロック1が較正データに対応し、ブロック2が交差検証データに対応する。
【0062】
クラスの予測の性能は、都合のよいことに、クラスの予測の感度の平均(よく分類されたスペクトルのレート)に等しくなるように計算される。この重みを付けられた尺度は、「バランス分類レート」すなわち「BCR:balance classification rate」とも呼ばれ、コロニーのサイズによるケースであり、種に応じて変化する、不安定なピクセルのセットを考慮できるようにする。BCRの計算は、
図11において思い出される。
【0063】
C.1.COS結果
C.1.1.フラットなモデル
下記のテーブル1は、
図7Aに示されたフラットな予測モデルに対する、および、関係(7)を使用して取得された予測モデルに対する、BCRを示す。
【0064】
種々の菌株が正確に予測されるということに、テーブル1を読むとすぐに気づくであろう。特に、特徴が明らかにされることになる微生物が細菌であると知っていると、第1の予測GP対GNNおよびGNF、ならびに第2の予測GNN対GPおよびGNFを使用して、細菌のグラムタイプ、および細菌の発酵または非発酵特性を予測することができる。この種類の予測は、出願人によって市販されているVitek(登録商標)2プラットフォームでアンチバイオグラムを生成するための消耗品を選択するのに特に有用である。
【0065】
C.1.2.最適木
図7Cに示され、関係(7)によって取得されたモデル110、112、114のBCRが、テーブル2に概説される。
【0066】
最適木は、系統樹とは著しく異なり、COS媒質の影響はおそらく、系統発生学上の差の影響より大きいということがわかる。
【0067】
図7Cに示され、
図4の「ステップフォワード」法、および、モデルのそれぞれに対してR=24である関係(6)によって取得されたモデル110、112、114のBCRが、テーブル3に概説される。
【0068】
性能の利得は、第1のモデルのための8番目のチャネルから、および、第3のモデルのための4番目のチャネルからスタートすることに限定されるということが、テーブル3を読むとわかる。市販のフィルタシステムに対応する24個のスペクトルフィルタを使用してマルチスペクトルアプリケーションを取得するために、都合のよいことに、第1、第2、および第3のモデル110、112、114のために8個のチャネル、14個のチャネル、および4個のチャネルがそれぞれ選択される。本実施形態のための性能値が、テーブル4に概説される。
【0069】
当然、
図2のシステムが利用できるスペクトルフィルタの数に応じて、他のチャネル数が選択されてもよい。
【0070】
「ステップフォワード」法は、クラスを予測するのに必要な情報を含むスペクトル範囲を決定できるようにするということにさらに留意されたい。各モデルの最初の5つのチャネルに限定すると、90%に近いまたは90%を超えるBCRが、それぞれ取得される。これらのチャネルのスペクトル分散が、
図12から
図14に示される。このようにして、50nmより大きく離れた異なるスペクトル帯域が、区別されることが可能である。具体的には、
A. 4つのクラスY、GP、GNF、およびGNNは、第1の範囲415nm~500nmおよび第2の範囲535nm~675nmにあるスペクトルデータだけを使用して効率的に予測されることが可能である。これらの範囲だけを使用すると、BCRは90%以上になる。575nm~675nmに第1の範囲を限定することによって、90%に近いまたは90%を超えるBCRについてモデルあたり4つのチャネルだけが使用される。任意選択として、第2のモデルのランク5のチャネルに対応する第3の範囲850nm~875nmが使用される。より具体的には、第1の範囲415nm~500nmでの予測は、単に、範囲415nm~440nmおよび470nm~495nmに基づいてもよい。したがって、本発明は、前記範囲のスペクトルを獲得し、前記範囲の前記スペクトルのみに応じてクラスを予測することからなるクラスY、GP、GNF、およびGNNを予測するためのいずれかの方法を含む。本発明が、クラスY、GP、GNF、およびGNNの間で区別できるようにすると、したがって、前述の範囲に含まれるスペクトルデータに基づいて酵母と細菌を区別できるようにするということにさらに留意されたい。したがって、本発明は、特徴が明らかにされることになる微生物の酵母または細菌の特徴を予測するための方法も含む。
B. クラスGNF対Y+GP+GNNの予測は、第1の範囲470nm~500nmおよび第2の範囲575nm~645nmに対してのみ効率的に行われてもよい。したがって、本発明は、前記範囲のスペクトルを獲得し、前記範囲の前記スペクトルのみに応じてクラスGNFを予測することからなるクラスGNFを予測するためのいずれかの方法を含む。特徴が明らかにされることになる微生物の細菌の特徴が既に知られているとき、予測は、クラスGNF対GPおよびGNNを予測することからなるということに留意されたい。したがって、本発明は、範囲470nm~500nmおよび575nm~645nmだけに基づくこの種類の予測も含む。
C. クラスGP対Y+GNNの予測は、第1の範囲535nm~675nm、および、より具体的には、範囲585nm~675nm、ならびに第2の範囲850nm~875nmだけに対して効率的に行われてもよい。したがって、本発明は、前記範囲のスペクトルを獲得し、前記範囲の前記スペクトルだけに応じてクラスGPを予測することからなるクラスGPを予測するためのいずれかの方法を含む。特徴が明らかにされることになる微生物の細菌の特徴が既に知られているとき、予測は、クラスGP対クラスGNFおよびGNNを、ならびに必然的に、クラスGP対グラム陰性細菌(GN)のクラスによって、予測することからなるということに留意されたい。したがって、本発明は、範囲535nm~675nm、およびより具体的には範囲585nm~675nm、ならびに範囲850nm~875nmだけに基づくこの種類の予測も含む。
D. 下記のポイントBおよびCに記述された予測を組み合わせると、したがって、3つの範囲と共に、また、特徴が明らかにされることになる微生物の細菌の特徴を知っていると、細菌のコロニーがGPであるか、GNFであるか、GNNであるかを判定することができるということがわかる。この種類の予測は、出願人によって市販されているVitek(登録商標)2プラットフォームでアンチバイオグラムを生成するための消耗品を選択するのに特に有用である。
【0071】
C.1.3.系統樹
図7Bに示され、関係(7)によって取得されたモデル100、102、104のBCRが、テーブル5に概説される。
【0072】
C.2.TSA結果
C.2.1.フラットなモデル
下記のテーブル1は、
図7Aに示されたフラットな予測モデルに対する、および、関係(7)によって取得された予測モデルに対する、BCRを示す。
【0073】
C.2.2.最適木
図7Cに示され、関係(7)によって取得されたモデル110、112、114に対するBCRが、テーブル7に概説される。
【0074】
最適木は、系統樹とは著しく異なり、COS媒質の影響はおそらく、系統発生学上の差の影響より大きいということがわかる。
【0075】
図7Cに示され、
図5の「ステップフォワード」法、および、モデルのそれぞれに対してR=24である関係(6)によって取得されたモデル110、112、114に対するBCRが、テーブル8に概説される。
【0076】
図15から
図17は、第1、第2、および第3のモデル110、112および114(
図16)に対する、ならびにモデル(
図17)ごとの、12個、8個、および11個のチャネルによる貧弱な手法での、全く同一の図(
図15)におけるブランチの主要なチャネルのスペクトル分散を示す。
【0077】
C.2.3.系統樹
図7Bに示され、関係(7)によって取得されたモデル100、102、104のBCRが、テーブル9に概説される。
【0078】
図7Bに示され、
図4の「ステップフォワード」法、および、モデルのそれぞれに対してR=24である関係(6)によって取得されたモデル100、102、104のBCRが、テーブル10に概説される。
【0079】
図18から
図20は、これらの図のそれぞれの一番上にTSA媒質があり、図のそれぞれの一番下にあるCOS媒質と比較して、モデルごとに主要なスペクトルチャネルのスペクトル分散を示す。