(57)【要約】      （修正有）【課題】レンズレスイメージングシステムを用いて、微生物コロニーの散乱光パターンから微生物を判別する方法の提供。【解決手段】フォトセンサを用いたレンズレスイメージングシステムを用いて微生物を判別する方法であって、フォトセンサ上に載置した培養容器中で、微生物を含むサンプルを培養し、培養期間中に形成されるコロニーを経時的に撮像する工程、複数のコロニー形成画像から抽出される特徴から、少なくともコロニー最大増殖速度μｍａｘ及び相対平均輝度Ｉを含む３種以上、又は少なくともヒストグラムの偏差Ｇ及びドーナツ性Ｄ、又はドーナツ性Ｄ及び画像エントロピーＨを含む２種以上のパラメータを算出する工程、前記パラメータの値を、予め構築された既知微生物ごとに算出したパラメータのデータベースと照合・比較して微生物を特定する工程、を含む方法。【選択図】図１
2018033430.tif 000002
【特許請求の範囲】
【請求項１】
フォトセンサを用いたレンズレスイメージングシステムを用いて微生物を判別する方法であって、
フォトセンサ上に載置した培養容器中で、微生物を含むサンプルを培養し、培養期間中に形成されるコロニーを経時的に撮像する工程、
複数のコロニー形成画像から抽出される特徴から、コロニー最大増殖速度μｍａｘ、コロニー可視化時間ｔａ、相対平均輝度Ｉ、ヒストグラムの偏差Ｇ、ドーナツ性Ｄ、画像エントロピーＨ、画像エネルギー密度Ｅd、画像エネルギーＥ、重み付け中心の差異Ｗ、及び中心領域の平均輝度Ｃの１０種のパラメータのうち、少なくともコロニー最大増殖速度μｍａｘ及び相対平均輝度Ｉを含む３種以上、又は少なくともヒストグラムの偏差Ｇ及びドーナツ性Ｄ、又はドーナツ性Ｄ及び画像エントロピーＨを含む２種以上を算出する工程、
前記パラメータの値を、予め構築された既知微生物ごとに算出した当該パラメータのデータベースと照合・比較して微生物を特定する工程、を含む方法。
【請求項２】
フォトセンサがＣＭＯＳセンサである、請求項１記載の方法。
【請求項３】
撮像が、５分間隔毎に１８時間行われる、請求項１又は２記載の方法。
【請求項４】
微生物の特定が、前記パラメータを説明変数として多変量解析することにより行われる、請求項１〜３の何れか１項記載の方法。
【請求項５】
多変量解析が、判別分析、ロジスティック回帰分析、サポートベクトルマシン、決定木、ランダムフォレスト、ニューラルネットワーク、又はクラスター分析である請求項４記載の方法。
【請求項６】
微生物が、エシェリヒア・コリ、スタフィロコッカス・アウレウス、シュードモナス エルギノーザ、サルモネラ・エンテリカ、カンジダ・アルビカンス、スタフィロコッカス・エピデルミディス、スタフィロコッカス・ヘモリチカス、スタフィロコッカス・サプロフィティカス、及びスタフィロコッカス・シミュランスから選ばれる何れかであることを判別する請求項１〜５の何れか１項記載の方法。
【請求項７】
算出されるパラメータが、コロニー最大増殖速度μｍａｘ、コロニー可視化時間ｔａ及び相対平均輝度Ｉの３種、ヒストグラムの偏差Ｇ及びドーナツ性Ｄ、又はドーナツ性Ｄ及び画像エントロピーＨの２種、コロニー最大増殖速度μｍａｘ、可視化時間ｔａ、平均相対輝度Ｉ、及びヒストグラムの偏差Ｇの４種、又はコロニー最大増殖速度μｍａｘ、相対平均輝度Ｉ、ヒストグラムの偏差Ｇ、ドーナツ性Ｄ、画像エントロピーＨ、画像エネルギー密度Ｅd、重み付け中心の差異Ｗ及び中心領域の平均輝度Ｃの８種である請求項１〜６の何れか１項記載の方法。
【請求項８】
算出されるパラメータが、コロニー最大増殖速度μｍａｘ、コロニー可視化時間ｔａ及び相対平均輝度Ｉの３種であり、微生物がエシェリヒア・コリ、スタフィロコッカス・アウレウス、シュードモナス エルギノーザ、サルモネラ・エンテリカ、及びカンジダ・アルビカンスの５菌種の何れかであることを判別する請求項１〜７の何れか１項記載の方法。
【請求項９】
算出されるパラメータが、コロニー最大増殖速度μｍａｘ、コロニー可視化時間ｔａ、相対平均輝度Ｉ及びヒストグラムの偏差Ｇの４種であり、微生物がエシェリヒア・コリ、スタフィロコッカス・アウレウス、シュードモナス エルギノーザ、サルモネラ・エンテリカ、カンジダ・アルビカンス及びスタフィロコッカス・エピデルミディスの６菌種の何れかであることを判別する請求項１〜７の何れか１項記載の方法。
【請求項１０】
算出されるパラメータが、ヒストグラムの偏差Ｇ及びドーナツ性Ｄ、又はドーナツ性Ｄ及び画像エントロピーＨの２種であり、微生物がスタフィロコッカス・アウレウス及びスタフィロコッカス・エピデルミディスの何れかであることを判別する請求項１〜７の何れか１項記載の方法。
【請求項１１】
算出されるパラメータが、コロニー最大増殖速度μｍａｘ、相対平均輝度Ｉ、ヒストグラムの偏差Ｇ、ドーナツ性Ｄ、画像エントロピーＨ、画像エネルギー密度Ｅd、重み付け中心の差異Ｗ、及び中心領域の平均輝度Ｃの８種であり、微生物がスタフィロコッカス・アウレウス、スタフィロコッカス・エピデルミディス、スタフィロコッカス・ヘモリチカス、スタフィロコッカス・サプロフィティカス、スタフィロコッカス・シミュランスの５菌種の何れかであることを判別する請求項１〜７の何れか１項記載の方法。
【請求項１２】
微生物のコロニー形成画像をピクセルデータとして取得するフォトセンサと、微生物を含む試料を培養するための培養容器、微生物コロニーに対して光照明する光源、及び演算装置を含む、請求項１〜１１のいずれか１項記載の微生物の判別を行うための装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、光センシングによる微生物の判別方法に関する。
【背景技術】
【０００２】
微生物の検出・同定試験は、食品や医薬品、化粧品の製造工程の管理や、医療分野における臨床診断を目的として広く行われている。微生物による健康被害を最小化し、さらに微生物によってもたらされる疾病の蔓延を防ぐため、簡便かつ迅速な微生物検出・同定法の開発が求められ、例えば、医薬品または化粧品の製造現場において、特定微生物試験が行われている。これは、製品中に検出可能な特定の５菌種（Ｅ．ｃｏｌｉ、Ｓｔａｐｈｙｌｏｃｏｃｃｕｓ  ａｕｒｅｕｓ、Ｐｓｅｕｄｏｍｏｎａｓ  ａｅｒｕｇｉｎｏｓａ、Ｓａｌｍｏｎｅｌｌａ  ｅｎｔｅｒｉｃａ、Ｃａｎｄｉｄａ  ａｌｂｉｃａｎｓ）が存在しないことを判定する試験である。本試験では、製品の増菌培養（１８時間〜５日間）により得られた培養液を、鑑別培地上に塗布・培養（１８〜７２時間）し、得られたコロニーの色や形状から、特定５菌種であると疑わしいコロニーを検出する培養法が用いられている。本手法は、微生物の検出能が優れていることから、長きにわたり微生物検査の公定法として用いられているが、正確な種同定を可能とする一方で、多くの専門知識と熟練が要求される。また培養に基づいている手法であるため同定の結果を得るまでに２−３日を要するといった問題点がある。
【０００３】
公定法に代わる方法として、複数の生化学的試験を一枚のストリップ上に集積化させ、対象菌種を同定するキットも市販され、生化学的性状試験を大幅に簡便化できるが、培養に基づく手法であるため結果が出るまでに最大で４８時間を要し、迅速化は達成していない。また、近年ではＭａｔｒｉｘ  Ａｓｓｉｓｔｅｄ  Ｌａｓｅｒ  Ｄｅｓｏｒｐｔｉｏｎ／Ｉｏｎｉｚａｔｉｏｎ  Ｔｉｍｅ  ｏｆ  Ｆｌｉｇｈｔ  Ｍａｓｓ  Ｓｐｅｃｔｒｏｍｅｔｒｙ（ＭＡＬＤＩ−ＴＯＦ  ＭＳ）を用いて、菌体構成成分の網羅的な分析結果を種同定の手段として利用する手法が報告されている (非特許文献１)。しかし、この方法ではサンプル調製に煩雑な操作を要する、また特別な試薬を使用するため運営コストが高いといった問題点がある。
【０００４】
斯かる状況の下、コロニーの前処理が不要であり、運営コストの低い微生物同定として、微生物コロニーの画像パターンから菌種判別を行う手法が報告され、注目を浴びている  （非特許文献２）。その代表例であるＢＡｃｔｅｒｉａｌ  Ｒａｐｉｄ  Ｄｅｔｅｃｔｉｏｎ  ｕｓｉｎｇ  Ｏｐｔｉｃａｌ  ｓｃａｔｔｅｒｉｎｇ  Ｔｅｃｈｎｏｌｏｇｙ（ＢＡＲＤＯＴ）法は、ディッシュ上のコロニー（直径：１．１±０．２ｍｍ）に波長６３５  ｎｍのレーザー光を照射し、背面に投影された散乱光のパターンから菌種判別を行う手法である。本手法は、上述の培養法と同様、同一の菌種は同一の特徴を持つコロニーを形成するという前提に基づいた方法である。古典的な培養法では目視観察による菌種判別を試みているのに対し、本手法では散乱光パターンから得られる様々なコロニーの特徴（コロニーの大きさや厚さ、色や表面の粗さ等）の違いを捉え、独自のアルゴリズムを用いて定量的に評価しているため、目視よりも格段に高い精度で菌種判別が可能である。
【０００５】
これまでに、コロニー由来散乱光パターンの画像から抽出したＺｅｒｎｉｋｅモーメント（二値化した画像の円形範囲内の形状などの量を表すデータ）に基づく判別パラメータを使用することで、Ｌｉｓｔｅｒｉａ  ｍｏｎｏｃｙｔｏｇｅｎｅｓ、Ｓｔａｐｈｙｌｏｃｏｃｃｕｓ  ａｕｒｅｕｓ、Ｓａｌｍｏｎｅｌｌａ  ｔｙｐｈｉｍｕｒｉｕｍ、Ｖｉｂｒｉｏ  ｐａｒａｈａｅｍｏｌｙｔｉｃｕｓ及びＥｓｃｈｅｒｉｃｈｉａ  ｃｏｌｉの異なる５属の細菌を９０〜９９％で判別可能であることが報告されている（非特許文献３）。また同属異種の細菌の判別にも応用されており、Ｌｉｓｔｅｒｉａ属６種を正答率９１％以上で判別した報告（非特許文献２）やＶｉｂｒｉｏ属３種を正答率９３％以上で判別した報告（非特許文献４）がある。
【０００６】
以上のように、ＢＡＲＤＯＴ法は高い分類精度を有するが、分類対象が近縁種の場合は正答率が９０％程度にまで低下するといった問題点があり、さらなる高精度化が必要である。高精度化には判別パラメータの増加が必要であると考えられる。しかし、ＢＡＲＤＯＴ法はスナップショットに基づく判別法であるため、１枚の散乱光パターンから得られる情報量には限界がある。さらに、レーザー照射装置や投影スクリーン、レーザー照射装置とプレート上のコロニーの位置合わせを行うアクチュエーターなどが必要であるため、装置は大型・高価となり、大規模検査機関などでの利用が主である。したがって、迅速な対策が要求される、食品や医薬品、化粧品の製造工程や医療機関での微生物汚染においては、オンサイトで利用可能な小型・安価で簡便な菌種判定システムの開発が求められている。
【０００７】
一方、近年、安価なシステムでありながら、ハイスループットに散乱光パターンを取得する手法として、２次元フォトセンサを用いたレンズレスイメージングが報告されている。２次元フォトセンサは、数百〜数千万個の検出素子が平面上にアレイ化されており、広範囲の光情報を高い時間分解能をもって２次元画像として記録できるという特徴を有する。この２次元フォトセンサ直上に細胞を配置し、可干渉性の光照射下で撮像を行う手法が細胞のレンズレスイメージングである。本手法は、対物レンズによる像の拡大や結像処理を経ずに細胞を認識、画像化することで、広域に存在する細胞由来散乱光パターンを取得可能である。例えば、ＬＥＤ光照射下のセンサ撮像領域上に配置した個々の細胞に由来する影や散乱光パターンを取得することにより、一度に広視野をイメージングする手法がある（非特許文献５）。
【先行技術文献】
【非特許文献】
【０００８】
【非特許文献１】Welker M, Moore ER. 2011. Applications of whole-cell matrix-assisted laser-desorption/ionization time-of-flight mass spectrometry in systematic microbiology. Syst Appl Microbiol 34(1):2-11.
【非特許文献２】Banada PP, Guo SL, Bayraktar B, Bae E, Rajwa B, Robinson JP, Hirleman ED, Bhunia AK. 2007. Optical forward-scattering for detection of Listeria monocytogenes and other Listeria species. Biosens Bioelectron 22(8):1664-1671.
【非特許文献３】Banada PP, Huff K, Bae E, Rajwa B, Aroonnual A, Bayraktar B, Adil A, Robinson JP, Hirleman ED, Bhunia AK. 2009. Label-free detection of multiple bacterial pathogens using light-scattering sensor. Biosens Bioelectron 24(6):1685-1692.
【非特許文献４】Huff K, Aroonnual A, Littlejohn AEF, Rajwa B, Bae E, Banada PP, Patsekin V, Hirleman ED, Robinson JP, Richards GP and others. 2012. Light-scattering sensor for real-time identification of Vibrio parahaemolyticus, Vibrio vulnificus and Vibrio cholerae colonies on solid agar plate. Microbial Biotechnol 5(5):607-620.
【非特許文献５】Ozcan A, Demirci U. 2008. Ultra wide-field lens-free monitoring of cells on-chip. Lab Chip 8(1):98-106.
【発明の概要】
【発明が解決しようとする課題】
【０００９】
本発明は、レンズレスイメージングシステムを用いて、微生物コロニーの散乱光パターンから微生物を判別する方法を提供することに関する。
【課題を解決するための手段】
【００１０】
本発明者らは、微生物を培養することにより形成されるコロニーを、レンズレスイメージングシステムを用いて経時的に撮像し、得られたコロニー形成画像の特徴から算出される特定のパラメータを組み合わせて用いることにより、微生物の判別が高精度で可能であることを見出した。
【００１１】
すなわち、本発明は、以下の１）〜１２）に係るものである。
１）フォトセンサを用いたレンズレスイメージングシステムを用いて微生物を判別する方法であって、
フォトセンサ上に載置した培養容器中で、微生物を含むサンプルを培養し、培養期間中に形成されるコロニーを経時的に撮像する工程、
複数のコロニー形成画像から抽出される特徴から、コロニー最大増殖速度μｍａｘ、コロニー可視化時間ｔａ、相対平均輝度Ｉ、ヒストグラムの偏差Ｇ、ドーナツ性Ｄ、画像エントロピーＨ、画像エネルギー密度Ｅd、画像エネルギーＥ、重み付け中心の差異Ｗ、及び中心領域の平均輝度Ｃの１０種のパラメータのうち、少なくともコロニー最大増殖速度μｍａｘ及び相対平均輝度Ｉを含む３種以上、又は少なくともヒストグラムの偏差Ｇ及びドーナツ性Ｄ、又はドーナツ性Ｄ及び画像エントロピーＨを含む２種以上を算出する工程、
前記パラメータの値を、予め構築された既知微生物ごとに算出した当該パラメータのデータベースと照合・比較して微生物を特定する工程、を含む方法。
２）フォトセンサがＣＭＯＳセンサである、１）の方法。
３）撮像が、５分間隔毎に１８時間行われる、１）又は２）の方法。
４）微生物の特定が、前記パラメータを説明変数として多変量解析することにより行われる、１）〜３）の何れかの方法。
５）多変量解析が、判別分析、ロジスティック回帰分析、サポートベクトルマシン、決定木、ランダムフォレスト、ニューラルネットワーク、又はクラスター分析である４）の方法。
６）微生物が、エシェリヒア・コリ、スタフィロコッカス・アウレウス、シュードモナス エルギノーザ、サルモネラ・エンテリカ、カンジダ・アルビカンス、スタフィロコッカス・エピデルミディス、スタフィロコッカス・ヘモリチカス、スタフィロコッカス・サプロフィティカス、及びスタフィロコッカス・シミュランスから選ばれる何れかであることを判別する１）〜５）の何れかの方法。
７）算出されるパラメータが、コロニー最大増殖速度μｍａｘ、コロニー可視化時間ｔａ及び相対平均輝度Ｉの３種、ヒストグラムの偏差Ｇ及びドーナツ性Ｄ、又はドーナツ性Ｄ及び画像エントロピーＨの２種、コロニー最大増殖速度μｍａｘ、可視化時間ｔａ、平均相対輝度Ｉ、及びヒストグラムの偏差Ｇの４種、又はコロニー最大増殖速度μｍａｘ、相対平均輝度Ｉ、ヒストグラムの偏差Ｇ、ドーナツ性Ｄ、画像エントロピーＨ、画像エネルギー密度Ｅd、重み付け中心の差異Ｗ及び中心領域の平均輝度Ｃの８種である１）〜６）の何れかの方法。
８）算出されるパラメータが、コロニー最大増殖速度μｍａｘ、コロニー可視化時間ｔａ及び相対平均輝度Ｉの３種であり、微生物がエシェリヒア・コリ、スタフィロコッカス・アウレウス、シュードモナス エルギノーザ、サルモネラ・エンテリカ、及びカンジダ・アルビカンスの５菌種の何れかであることを判別する１）〜７）の何れかの方法。
９）算出されるパラメータが、コロニー最大増殖速度μｍａｘ、コロニー可視化時間ｔａ、相対平均輝度Ｉ及びヒストグラムの偏差Ｇの４種であり、微生物がエシェリヒア・コリ、スタフィロコッカス・アウレウス、シュードモナス エルギノーザ、サルモネラ・エンテリカ、カンジダ・アルビカンス及びスタフィロコッカス・エピデルミディスの６菌種の何れかであることを判別する１）〜７）の何れかの方法。
１０）算出されるパラメータが、ヒストグラムの偏差Ｇ及びドーナツ性Ｄ、又はドーナツ性Ｄ及び画像エントロピーＨの２種であり、微生物がスタフィロコッカス・アウレウス及びスタフィロコッカス・エピデルミディスの何れかであることを判別する１）〜７）の何れかの方法。
１１）算出されるパラメータが、コロニー最大増殖速度μｍａｘ、相対平均輝度Ｉ、ヒストグラムの偏差Ｇ、ドーナツ性Ｄ、画像エントロピーＨ、画像エネルギー密度Ｅd、重み付け中心の差異Ｗ、及び中心領域の平均輝度Ｃの８種であり、微生物がスタフィロコッカス・アウレウス、スタフィロコッカス・エピデルミディス、スタフィロコッカス・ヘモリチカス、スタフィロコッカス・サプロフィティカス、スタフィロコッカス・シミュランスの５菌種の何れかであることを判別する１）〜７）の何れかの方法。
１２）微生物のコロニー形成画像をピクセルデータとして取得するフォトセンサと、微生物を含む試料を培養するための培養容器、微生物コロニーに対して光照明する光源、及び演算装置を含む、１）〜１１）の何れかの微生物の判別を行うための装置。
【発明の効果】
【００１２】
本発明の方法によれば、安価で小型の装置を用いて、高精度に微生物を判別・同定可能なことから、化粧品、食品、医薬品等の製造分野や、臨床診断における微生物検査に導入することにより、微生物同定試験を大幅に簡便、迅速化できる。
【図面の簡単な説明】
【００１３】
【図１】本発明のレンズレスイメージングシステムの構成を示す断面図。
【図２】培養容器（チャンバー）の構築工程。
【図３】タイムラプス撮像された各微生物のコロニー形成画像。
【図４】３種のパラメータを用いた特定５菌種のクラスター分析結果。
【図５】２種のパラメータを用いた非階層的クラスター分析結果。
【図６】４−７種のパラメータを用いた非階層的クラスター分析（Ｗａｒｄ法）結果。
【図７】８種のパラメータを用いた線形判別分析の結果。
【発明を実施するための形態】
【００１４】
本発明の微生物の判別方法は、フォトセンサを用いたレンズレスイメージングシステムを用い、当該システムにより撮像して得られる画像（「レンズレス画像」）を解析するものである。
ここで、レンズレスイメージングシステムとは、対物レンズを介さない光センシングシステムを意味する。
本発明の微生物の判別方法において用いられるレンズレスイメージングシステムとしては、微生物のコロニー形成画像（散乱光パターン）をピクセルデータとして取得するフォトセンサと、微生物を含む試料を培養するための培養容器、微生物コロニーに対して光照明する光源、及び、コロニー形成画像のピクセルデータから定量的パラメータを抽出し、それらを用いて多変量解析を行うための演算装置から構成される。これに、必要に応じて培養温度を管理するインキュベーターを備えてなるデバイスは、本発明の微生物の判別方法を実行するための装置となり得る。図１に、本発明の方法を実行するための装置の一構成例を示す。
【００１５】
図１に示されるように、フォトセンサ２の上に、培養容器１が載置され、培養容器１の上方には、点光源３が配置されている。尚、培養容器はコロニー形成画像が撮像できるようにフォトセンサ上に載置されればよく、図１に示すようにフォトセンサ２に対して、微生物培養面（コロニー形成面）を向けるように載置してもよく、或いはコロニー形成面を上に向けて載置することでもよい。
【００１６】
ここで、培養容器としては、微生物の培養が可能で、培養面をフォトセンサに向けて積載可能なものであれば特に制限はされないが、例えば、顕微鏡観察用スライドガラス上にフレームシール（例えば、剥離フィルム及びカバーフィルムによって覆われているｉｎ  ｓｉｔｕ  ＰＣＲシール（例えば、厚さ：３００μｍ、Ｆｒａｍｅ−Ｓｅａｌ  Ｉｎｃｕｂａｔｉｏｎ  Ｃｈａｍｂｅｒｓ，ＢＩＯ  ＲＡＤ））を数枚（例えば２枚）重ねて接着形成されるチャンバーが好適に使用できる。図２に当該チャンバーの構築工程を示す。
【００１７】
フォトセンサとしては、形成された微生物コロニーに対して光照射した場合に、当該コロニー由来の光シグナル又は影シグナルを撮像できるものであればよく、例えば、ＣＣＤ（Ｃｈａｒｇｅ  Ｃｏｕｐｌｅｄ  Ｄｅｖｉｃｅ）、ＣＭＯＳ（Ｃｏｍｐｌｅｍｅｎｔａｒｙ  Ｍｅｔａｌ  Ｏｘｉｄｅ  Ｓｅｍｉｃｏｎｄｕｃｔｏｒ）、ダブルゲート構造のＴＦＴ（Ｔｈｉｎ  Ｆｉｌｍ  Ｔｒａｎｓｉｓｔｏｒ）、ＣＩＳ  Ｃｏｎｔａｃｔ Ｉｍａｇｅ Ｓｅｎｓｏｒ）等の撮像素子を用いることができるが安価であり、解像度が高い点からＣＭＯＳセンサを用いるのが好ましい。
ＣＭＯＳセンサとしては、検出素子サイズが２〜８μｍ四方であり、素子数が３０万〜１８１０万のものでカラーイメージングが可能なものが望ましく、例えば３．２μｍ四方の検出素子が２０４８×１５３６に配列し、撮像エリアのサイズが６．５５ｍｍ×４．９２ｍｍであるＤＦＫ６１ＢＵＣ０２（Ｉｍａｇｉｎｇ  Ｓｏｕｒｃｅ社，Ｇｅｒｍａｎｙ）などを挙げることができる。
【００１８】
照明用光源としては、ＬＥＤ、有機ＥＬ、蛍光灯、白熱球等が使用できるが、ＬＥＤを用いるのが好ましく、コロニー画像のコントラストを高くできる点から波長４００ｎｍ以上５００ｎｍ以下の範囲のＬＥＤを用いるのがより好ましい。
【００１９】
本発明の方法においては、フォトセンサアレイ上に載置した培養容器中で、微生物を含むサンプルが培養される。
ここで、「微生物を含む試料」は、菌種の判別・同定対象の微生物を含む試料であり、臨床試料、非臨床試料の何れでも良い。臨床試料としては、例えば血液、血清、血漿、血液分画、関節液、尿、精液、唾液、糞便、脳脊髄液、胃内容物、膣分泌物、組織ホモジェネート、骨髄穿刺液、骨ホモジェネート、痰、吸引液、ぬぐい液（swab）およびぬぐい液残滓（swab rinsate）、他の体液等が挙げられる。また、非臨床試料としては、食品、飲料、医薬品、化粧品、水、海水バラスト、空気、土壌、汚水、植物材料、血液製剤、ドナー臓器または組織試料等を含めた物質が挙げられる。
【００２０】
本発明において、「微生物」としては、コロニーを形成するものであれば特に限定されず、細菌、真菌等の何れでもよい。例えば、シュードモナス属、エシェリキア属、サルモネラ属、赤痢菌属、エンテロバクター属、クレブシエラ属、セラチア属、プロテウス属、カンピロバクター属、ヘモフィルス属、モルガネラ属、ビブリオ属、エルシニア属、アシネトバクター属、ステノトロフォモナス属、ブレブンディモナス属、ラルストニア属、アクロモバクター属、フゾバクテリウム属、プレボテラ属、ブランハメラ亜属、ナイセリア属、バークホルデリア属、シトロバクター属、ハフニア属、エドワードシエラ属、アエロモナス属、モラクセラ属、ブルセラ属、パスツレラ属、プロビデンシア属およびレジオネラ属等のグラム陰性細菌；腸球菌、連鎖球菌、ブドウ球菌、バチルス属、パエニバチルス属、乳酸桿菌属、リステリア属、ペプトストレプトコッカス属、プロピオン酸菌属、クロストリジウム属、バクテロイデス属、ガードネレラ属、コクリア属、ラクトコッカス属、ロイコノストック属、ミクロコッカス、マイコバクテリウム属およびコリネバクテリウム属等のグラム陽性細菌；カンジダ属、クリプトコックス属、ノカルジア属、アオカビ属、アルタナリア属、ロドトルラ属、アスペルギルス属、フザリウム属、サッカロミセス属およびトリコスポロン属等の真菌が挙げられる。
【００２１】
培養条件は適宜設定すればよく、例えば室温で、１８時間が挙げられる。撮像は、当該培養時間中、一定間隔、例えば１〜６０分間隔で行えばよいが、好ましくは５分間隔が挙げられる。
撮像条件は、コロニー形成画像の特徴が正確に抽出できれば特に限定されないが、例えば、露光時間は１／１００〜１／４秒、好ましくは１／１８〜１／１０秒、フレームレートは３．７５〜１０、好ましくは３.７５〜５とすることができ、ホワイトバランスは、例えばＲｅｄ：２５５、Ｇｒｅｅｎ：１５４、Ｂｌｕｅ：６４とすることが挙げられる。
また、ノイズの削減のため、画像を複数枚連続取得し、その平均画像を使用することが好ましい。
【００２２】
コロニー形成画像から抽出した特徴、例えば、面積、輝度、形状等から、コロニー最大増殖速度μｍａｘ、コロニー可視化時間ｔａ（Ｃｏｌｏｎｙ  ａｐｐｅａｒａｎｃｅ  ｔｉｍｅ）、相対平均輝度Ｉ、ヒストグラムの偏差Ｇ、ドーナツ性Ｄ、画像エントロピーＨ、及び画像エネルギー密度Ｅd、画像エネルギーＥ、重み付け中心の差異Ｗ、中心領域の平均輝度Ｃの計１０種のパラメータから選ばれる、少なくともコロニー最大増殖速度μｍａｘ及び相対平均輝度Ｉを含む３種以上、又は少なくともヒストグラムの偏差Ｇ及びドーナツ性Ｄ、又はドーナツ性Ｄ及び画像エントロピーＨを含む２種以上が算出される。
【００２３】
ここで、コロニー最大増殖速度μｍａｘは、時刻ｔにおけるコロニーの面積をＡｒｅａｔとすると、コロニー増殖速度μ（ｈ−１）は式（１）で表される。ｔ＝１０ｈまでにおける最大値をμｍａｘとした。
【００２４】
【数１】
2018033430.tif 000003
【００２５】
コロニー可視化時間ｔａは、レンズレス画像中のコロニー非存在領域の直線（１００ｐｘ）上の輝度値の最大値と最小値の差をＮとし、コロニー由来散乱光パターンの中心を通る直線（１００ｐｘ）上の最大値と最小値の差をＳとしたとき、Ｓ／Ｎ≧２となった最小の時刻ｔと定義される。
【００２６】
相対平均輝度Ｉは、光源やフォトセンサの個体差による輝度の標準化を行うために、コロニー由来散乱光パターン領域内の平均輝度値から各レンズレス画像のコロニー非存在領域の中から選んだ１００ｐｘ×１００ｐｘの領域の平均輝度値を引いた値と定義される。
【００２７】
ヒストグラムの偏差Ｇは、コロニーの散乱光パターン画像から輝度値を２０分割することで得たヒストグラムの各ピクセル数の標準偏差と定義される。
【００２８】
【数２】
2018033430.tif 000004
（式中、ｎは分割（ヒストグラムのビン）の数である２０、ｘiは各ビンに含まれるピクセル数、はｘ￣は各ビンに含まれるピクセル数の平均値を示す。）
【００２９】
ドーナツ性Ｄは、コロニー全体の平均相対輝度値Ｄｗと、コロニー中心領域（直径１／２）の平均相対輝度値Ｃの比率を示すパラメータであり、Ｄｗ／Ｃと定義される。
【００３０】
画像エントロピーＨは、一般的に画像の情報量を表す指標として用いられており、輝度勾配の共起行列と定義される。
【００３１】
【数３】
2018033430.tif 000005
（式中、Pijは２５６階調に補正した輝度値がiのピクセルに隣接するピクセルの輝度値がjである確率を示す。）
【００３２】
画像エネルギー密度Ｅdは、２次元離散フーリエ変換により求めた振幅スペクトルの平均として定義される。
【００３３】
画像エネルギーＥは、２次元離散フーリエ変換により求めた振幅スペクトルの総和として定義される。
【００３４】
重み付け中心の差異Ｗは、輝度値を考慮せずに求めたコロニーの中心と、輝度値を考慮して求めたコロニーの中心との距離と定義する。
【００３５】
中心領域の平均輝度Ｃは、コロニー中心領域（直径１／２）の平均相対輝度値である。
【００３６】
微生物の判別は、菌種不明の微生物のコロニー形成画像から算出された前記パラメータの値を、予め構築された既知微生物ごとの当該パラメータのデータベースと照合・比較することにより行われる。データベースは、多くの菌種について、上述した方法によりコロニー形成画像を撮像し、得られた画像ライブラリに基づいて前記パラメータを算出しデータベース化したものである。ここで、データベースを構築する際は、同一菌種においても複数の株からデータを収集するのが好ましい。
【００３７】
構築したデータベースから判別情報を得る方法としては、データベースの各パラメータを説明変数として多変量解析する方法が挙げられる。本発明で用いる多変量解析としては、特に制限はないが、判別分析、二次判別分析、ロジスティック回帰分析、サポートベクトルマシン、決定木、ランダムフォレスト、ニューラルネットワーク、又はクラスター分析などを採用できる。
具体的な例としては、構築したデータベース中のコロニー最大増殖速度μｍａｘ、コロニー可視化時間ｔａ、相対平均輝度Ｉ、ヒストグラムの偏差Ｇ、ドーナツ性Ｄ、画像エントロピーＨ、及び画像エネルギー密度Ｅd、画像エネルギーＥ、重み付け中心の差異Ｗ、中心領域の平均輝度Ｃからなる１０種のパラメータのうち、少なくともコロニー最大増殖速度μｍａｘ及び相対平均輝度Ｉを含む３種以上、又は少なくともヒストグラムの偏差Ｇ及びドーナツ性Ｄ、又はドーナツ性Ｄ及び画像エントロピーＨを含む２種以上のパラメータを説明変数として判別分析を行って判別関数を得ることにより菌種の判別を行うことができる。
【００３８】
判別に使用するために算出されるパラメータは、適宜増減可能であり、組み合わせにより判別精度が異なる。そのため、パラメータは分析方法や判別対象となる菌種を考慮して適宜決定すればよいが、少なくともコロニー最大増殖速度μｍａｘ、相対平均輝度Ｉ、及びコロニー可視化時間ｔａを含む３種以上を含むのがより好ましい。
例えば、エシェリヒア・コリ（Ｅ．ｃｏｌｉ）、スタフィロコッカス・アウレウス（Ｓ．ａｕｒｅｕｓ）、シュードモナス エルギノーザ（Ｐ．ａｅｒｕｇｉｎｏｓａ）、サルモネラ・エンテリカ（Ｓ．ｅｎｔｅｒｉｃａ）、カンジダ・アルビカンス（Ｃ．ａｌｂｉｃａｎｓ）の５菌種の判別においては、コロニー最大増殖速度μｍａｘ、コロニー可視化時間ｔａ及び相対平均輝度Ｉの３種を用いることによりその判別が可能である。また、当該コロニー最大増殖速度μｍａｘ、コロニー可視化時間ｔａ及び相対平均輝度Ｉに、ヒストグラムの偏差Ｇを加えた４種を使用することで、上記５菌種にスタフィロコッカス・エピデルミディス（Ｓ．ｅｐｉｄｅｒｍｉｄｉｓ）を加えた６菌種の判別が可能となる。
また、ヒストグラムの偏差Ｇ及びドーナツ性Ｄ、又はドーナツ性Ｄ及び画像エントロピーＨの２種を使用すると、スタフィロコッカス・アウレウス及びスタフィロコッカス・エピデルミディスの判別が可能である。
また、スタフィロコッカス属の近縁５種（スタフィロコッカス・アウレウス（Ｓ．ａｕｒｅｕｓ）、スタフィロコッカス・エピデルミディス（Ｓ．ｅｐｉｄｅｒｍｉｄｉｓ）、スタフィロコッカス・ヘモリチカス（Ｓ．ｈａｅｍｏｌｙｔｉｃｕｓ）、スタフィロコッカス・サプロフィティカス（Ｓ．ｓａｐｒｏｐｈｙｔｉｃｕｓ）、スタフィロコッカス・シミュランス （Ｓ．ｓｉｍｕｌａｎｓ））の判別が、コロニー最大増殖速度μｍａｘ、相対平均輝度Ｉ、ヒストグラムの偏差Ｇ、ドーナツ性Ｄ、画像エントロピーＨ、画像エネルギー密度Ｅd、重み付け中心の差異Ｗ及び中心領域の平均輝度Ｃの８種を用いることにより可能である。
【実施例】
【００３９】
以下、実施例を用いて本発明をより詳細に説明するが、本発明の技術的範囲はこれら実施例に限定されるものではない。
【００４０】
＜使用機器・試薬＞
撮像素子としてＣＭＯＳセンサ（ＤＦＫ６１ＢＵＣ０２，Ｉｍａｇｉｎｇ  Ｓｏｕｒｃｅ，Ｇｅｒｍａｎｙ）を使用した。本センサは、３．２μｍ四方の検出素子が２０４８×１５３６に配列した２次元フォトセンサであり、撮像エリアのサイズは６．５５ｍｍ×４．９２ｍｍである。ＣＭＯＳセンサによる撮像には同社の制御用ソフトウェア（ＩＣ  Ｃａｐｔｕｒｅ  ２．３）を利用した。
また、画像解析用ソフトウェアとしてＭＡＴＬＡＢ  （ＭａｔｈＷｏｒｋｓ  ｉｎｃ．，ＵＳＡ）を使用した。コロニーの面積、ラインプロファイル、輝度値の計測にはＩｍａｇｅ  Ｊを、クラスター分析には統計解析ソフトＲを使用した。
照明用光源として、Ｌｉｇｈｔ  Ｅｍｉｔｔｉｎｇ  Ｄｉｏｄｅ  （ＬＥＤ）（ＬＭ１−ＴＰＰ１−０１，ＣＯＴＣＯ，Ｃｈｉｎａ）を使用した。
また、寒天培地作製には、ＬＢ培地（Ｍｅｒｃｋ  Ｍｉｌｌｉｐｏｒｅ，Ｇｅｒｍａｎｙ）及び低融点アガロース（Ｉｎｖｉｔｒｏｇｅｎ，  ＵＳＡ）を使用した。
【００４１】
＜使用菌株および培養方法＞
細胞培養及びコロニー形成評価には、Ｅｓｃｈｅｒｉｃｈｉａ  ｃｏｌｉ  ＡＴＣＣ  ８７３９株、Ｓｔａｐｈｙｌｏｃｏｃｃｕｓ  ａｕｒｅｕｓ  ＡＴＣＣ６５３８株、Ｐｓｅｕｄｏｍｏｎａｓ  ａｅｒｕｇｉｎｏｓａ  ＡＴＣＣ９０２７株、Ｓａｌｍｏｎｅｌｌａ  ｅｎｔｅｒｉｃａ  ＮＢＲＣ１００７９７株、Ｃａｎｄｉｄａ  ａｌｂｉｃａｎｓ  ＡＴＣＣ１０２３１株、Ｓｔａｐｈｙｌｏｃｏｃｃｕｓ  ｅｐｉｄｅｒｍｉｄｉｓ  ＡＴＣＣ１４９９０株、Ｓ．ｈａｅｍｏｌｙｔｉｃｕｓ  ＡＴＣＣ２９９７０株、Ｓ．ｓａｐｒｏｐｈｙｔｉｃｕｓ  ＡＴＣＣ１５３０５株、Ｓ．ｓｉｍｕｌａｎｓ  ＡＴＣＣ２７８４８株の計９株を使用した。
グリセロール中に凍結保存されているＥ．ｃｏｌｉ、Ｓ．ａｕｒｅｕｓ、Ｐ．ａｅｒｕｇｉｎｏｓａ、Ｓ．ｅｎｔｅｒｉｃａ、Ｓ．ｅｐｉｄｅｒｍｉｄｉｓ、Ｓ．ｈａｅｍｏｌｙｔｉｃｕｓ、Ｓ．ｓａｐｒｏｐｈｙｔｉｃｕｓ、Ｓ．ｓｉｍｕｌａｎｓをＬＢ培地５ｍｌに融解後、３７℃にて１５−１８時間振とう培養することで菌体培養液を得た。Ｃ．ａｌｂｉｃａｎｓのグリセロールストックはＹＰＤ培地を用いて融解し、上記と同様に菌体培養液を得た。得られた菌体培養液中の菌体数は、バクテリアカウンター（Ｂａｃｔｅｒｉａ  Ｃｏｕｎｔｅｒ，Ｊａｐａｎ）及び位相差顕微鏡（ＥＣＬＩＰＳＥ  Ｅ２００，Ｎｉｋｏｎ，Ｊａｐａｎ）を用いてカウントした。その後、ＬＢ培地を溶媒とし、５×１０４  ｃｅｌｌｓ／ｍｌとなるように調製し、プレーティング用菌体懸濁液とした。
【００４２】
実施例１  レンズレスイメージングシステムの構築
（１）顕微鏡観察用スライドガラス（厚さ：１．２ｍｍ）に両面のフィルムを剥がしたｉｎ  ｓｉｔｕ  ＰＣＲシール（厚さ：３００μｍ、Ｆｒａｍｅ−Ｓｅａｌ  Ｉｎｃｕｂａｔｉｏｎ  Ｃｈａｍｂｅｒｓ，ＢＩＯ  ＲＡＤ）を貼り、その上にカバーフィルムのみを剥がした２枚目のｉｎ  ｓｉｔｕ  ＰＣＲシールを重ねて貼ることで、９ｍｍ角、高さ６００μｍのチャンバーを作製した（図２参照）。
作製したチャンバーに１．５％（ｗ／ｖ）低融点アガロース添加ＬＢ培地を流し込み、上からカバーガラス（厚さ：１７０μｍ）を載せ寒天培地表面を平らにした後、２０分間室温で放置することでＬＢ寒天培地を固化した。カバーガラス及び剥離フィルムを剥がした後、寒天培地上に菌体懸濁液１μｌ（５×１０４  ｃｅｌｌｓ／ｍｌ）を滴下し、再びカバーガラスを載せることで、プレーティングを行った。
【００４３】
（２）スペーサーシール（厚さ：３００μｍ）を貼ったＣＭＯＳセンサ保護ガラス上に、（１）で作製した培養チャンバーのカバーガラス面を下に向けて設置した（図１参照）。卓上インキュベーター（＃１−７６０５−０１，ＡＳＯＮＥ）を用いて培養チャンバーを３７°Ｃに保持し、培養チャンバーから高さ１０ｃｍの位置に置いたＬＥＤを点光源とした。
【００４４】
実施例２  特定５菌種コロニーの判別
（１）レンズレス画像ライブラリの構築
特定５菌種（Ｅ．ｃｏｌｉ、Ｓ．ａｕｒｅｕｓ、Ｐ．ａｅｒｕｇｉｎｏｓａ、Ｓ．ｅｎｔｅｒｉｃａ、Ｃ．ａｌｂｉｃａｎｓ）をそれぞれ播種した直後の培養チャンバーをＣＭＯＳセンサに載せ、３７°Ｃ条件下で５分間隔、約１８時間連続でタイムラプス撮像を行った。
撮像条件は、露光時間：１／１８秒、フレームレート：３．７５、ホワイトバランスＲｅｄ：２５５、Ｇｒｅｅｎ：１５４、Ｂｌｕｅ：６４とした。また、ノイズの削減のため、画像を４枚連続取得し、その平均画像を使用した。画像取得を開始した時間をｔ＝０ｈとした。
その結果、Ｅ．ｃｏｌｉ、Ｓ．ａｕｒｅｕｓ、Ｐ．ａｅｒｕｇｉｎｏｓａ、Ｓ．ｅｎｔｅｒｉｃａについては、プレーティング直後に像が観察されていなかった位置から、約１−５時間後に各菌種コロニー由来だと考えられる散乱光パターンが観察された（図３Ａ−Ｄ）。その後、培養時間の経過に伴い、レンズレス画像上のコロニー由来散乱光パターンが水平方向へ広がっていく変化が確認された。また、Ｃ．ａｌｂｉｃａｎｓは培養チャンバーへのプレーティング直後からパターンが観察され、その後水平方向への像の拡大が観察された（図３Ｅ）。これは、Ｅ．ｃｏｌｉ等の細菌のサイズが数μｍ程度であるのに対し、真核生物であるＣ．ａｌｂｉｃａｎｓの細胞サイズは直径約１０μｍと大きく、また出芽によって細胞のサイズが大きくなった状態で播種されていることに起因すると考えられた。２台のＣＭＯＳセンサを使用して得られたコロニーの画像の内、ランダムに選んだ計１５コロニー分のレンズレス画像をライブラリとした。
【００４５】
（２）判別パラメータに基づくクラスター分析
（１）で構築したライブラリ画像から、特定５菌種間においてコロニーのサイズの推移に差異が見られた。そこで、コロニーサイズの経時変化に基づく特定５菌種の判別パラメータとして、コロニー最大増殖速度μｍａｘ、コロニー可視化時間ｔａ（Ｃｏｌｏｎｙ  ａｐｐｅａｒａｎｃｅ  ｔｉｍｅ）、及び相対平均輝度Ｉ（定義式は前記のとおり）を用い、ｋ−ｍｅａｎｓ法によるクラスター分析を行い、特定５菌種の判別を試みた。
ｉ）判別パラメータの取得には、光源であるＬＥＤやイメージングを担うＣＭＯＳセンサの個体差による輝度の標準化を行うために、レンズレス画像上の輝度値の下位１％を０、上位１％を２５５とし、コントラストを調整した。さらに、レンズレス画像上の培地領域（コロニー非存在領域）と比較して、低輝度であるコロニーの散乱光パターン領域を白色として抽出するため、画像の白黒反転を行った。その後、画像の輝度値を１０段階で分類し、最大の値を閾値として設定することにより、２値化画像を作成した。次に、得られた２値化画像の膨張・収縮処理を行うことにより、オブジェクトの平滑化・穴埋め処理、及びノイズ除去を行い、コロニー領域を抽出した。その後、ＩｍａｇｅＪを用いてコロニー最大増殖速度μｍａｘ、コロニー可視化時間ｔａ、及び相対平均輝度Ｉを算出した。
ｉｉ）クラスター分析は、統計分析ソフトＲを用い、コロニー最大増殖速度μｍａｘ、コロニー可視化時間ｔａ、及び相対平均輝度Ｉの３つを含むデータセットに対し、値の標準化を行った。解析コロニー数をｎ、時刻ｔにおける各パラメータの値をｘｔとすると、標準化後の値Ｚは式（４）で定義される。
【００４６】
【数４】
2018033430.tif 000006
【００４７】
標準化した各値に対し、ｋ−ｍｅａｎｓ法によるクラスター分析を行った（クラスター数：５）。
ｉｉｉ）その結果、全７５個のデータセットを１５個からなる５つのクラスターに分類することができ、正答率は１００％（７５／７５）となった（図４）。
以上より、最大コロニー増殖速度μｍａｘ及び可視化時間ｔａ、平均相対輝度Ｉの３つの判別パラメータを組み合わせることによって、特定５菌種を判別可能であることが示唆された。
【００４８】
実施例３  同属近縁種コロニーの判別
同属近縁種であるＳ．ａｕｒｅｕｓ及びＳ．ｅｐｉｄｅｒｍｉｄｉｓコロニーの判別を試みた。
まず、Ｓ．ｅｐｉｄｅｒｍｉｄｉｓ由来コロニー形成過程のレンズレス画像を実施例２と同様にして取得し、レンズレス画像をライブラリに追加した。
実施例２で用いたμｍａｘ、ｔａ、Ｉの３つのパラメータを用いてＳ．ａｕｒｅｕｓ及びＳ．ｅｐｉｄｅｒｍｉｄｉｓの判別を試みた。その結果、Ｓ．ａｕｒｅｕｓコロニー５個及びＳ．ｅｐｉｄｅｒｍｉｄｉｓコロニー１個から構成されるクラスター１と、Ｓ．ａｕｒｅｕｓコロニー１０個及びＳ．ｅｐｉｄｅｒｍｉｄｉｓコロニー１４個から構成されるクラスター２に分類され、正答率は６３．３％となった。
そこで、Ｓ．ａｕｒｅｕｓ及びＳ．ｅｐｉｄｅｒｍｉｄｉｓの画像ライブラリに対して、ＭＡＴＬＡＢプログラムを用いて、ヒストグラムの偏差Ｇ、ドーナツ性Ｄ、画像エントロピーＨの３つを判別パラメータとして抽出した（定義式は前記のとおり）。ヒストグラムの偏差Ｇは、コロニーの散乱光パターンの輝度値の分布の違いを表すパラメータである。散乱光パターンのヒストグラムから得られたピクセル数の標準偏差として定義した。ヒストグラムの幅が狭い場合、Ｇは高い値を示すため、Ｓ．ｅｐｉｄｅｒｍｉｄｉｓが高い値を示すことが想定される。実際に、直径約２５０μｍのコロニー由来散乱光パターンに対し、Ｇを算出した結果、Ｓ．ａｕｒｅｕｓ、Ｓ．ｅｐｉｄｅｒｍｉｄｉｓはそれぞれ２８２．７±３２．６、４４９．５±９５．６となり、両菌種間で差異が確認された。ドーナツ性Ｄは、コロニー全体の平均相対輝度値と、コロニー中心領域（直径１／２）の輝度の比率を示すパラメータである。白黒反転後の散乱光パターンの中心領域の輝度値が高い場合、ドーナツ性Ｄは小さくなる。Ｓ．ａｕｒｅｕｓ、Ｓ．ｅｐｉｄｅｒｍｉｄｉｓはそれぞれ１．３０±０．１１、０．９４±０．０３となった。画像エントロピーＨは、パターン内部の輝度値の規則正しさ（ランダム性の低さ）を表すパラメータであり、輝度勾配の共起行列に基づき算出される。画像エントロピーＨを算出した結果、Ｓ．ａｕｒｅｕｓ、Ｓ．ｅｐｉｄｅｒｍｉｄｉｓはそれぞれ５．７３±０．１４、５．１２±０．２６となり、差異が確認された。
そこで、ヒストグラムの偏差Ｇとドーナツ性Ｄ、ドーナツ性Ｄと画像エントロピーＨの組み合わせでｋ−ｍｅａｎｓ法によるクラスター分析を行い、Ｓ．ａｕｒｅｕｓ及びＳ．ｅｐｉｄｅｒｍｉｄｉｓの判別を試みた（クラスター数：２）。その結果、Ｓ．ａｕｒｅｕｓコロニー１５個から構成されるクラスター１と、Ｓ．ｅｐｉｄｅｒｍｉｄｉｓコロニー１５個から構成されるクラスター２に分類され、正答率は１００％となった（図５）。
【００４９】
実施例４  クラスター分析法の検討
クラスター分析に使用するパラメータは、数の増加に伴い計算量が増加し、結果が出るまでの時間が長くなる。そのため、パラメータ数はより少ない方が望ましい。
そこで、一般的にクラスター分析手法の中で高精度であるとされているＷａｒｄ法を用いて階層的クラスター分析を行い、パラメータ数を減らすことが可能かどうか検証した。
７パラメータ（コロニー最大増殖速度μｍａｘ、可視化時間ｔａ、平均相対輝度Ｉ、ヒストグラムの偏差Ｇ、ドーナツ性Ｄ、画像エントロピーＨ、画像エネルギー密度Ｅd）を用いて、階層的クラスター分析を行った。このとき、Ｓ．ａｕｒｕｅｓコロニーのデータセットをナンバー１−１５、Ｓ．ｅｐｉｄｅｒｍｉｄｉｓコロニーのデータセットをナンバー１６−３０とした。その結果、４つの判別パラメータを使用した時点で１００％の精度で判別可能となった（図６）。
以上より、本発明の方法は、μｍａｘ、ｔａ、Ｉ、Ｇの４パラメータを用いて、Ｗａｒｄ法による階層的クラスター分析を行うことにより、Ｓ．ａｕｒｅｕｓ及びＳ．ｅｐｉｄｅｒｍｉｄｉｓを１００％の精度で判別可能である事が示唆された。
【００５０】
実施例５  判別分析を用いた微生物の判別
ＣＭＯＳイメージセンサ上にＳｔａｐｈｙｌｏｃｏｃｃｕｓ ａｕｒｅｕｓ、 Ｓｔａｐｈｙｌｏｃｏｃｃｕｓ ｅｐｉｄｅｒｍｉｄｉｓ、 Ｓｔａｐｈｙｌｏｃｏｃｃｕｓ ｈａｅｍｏｌｙｔｉｃｕｓ、 Ｓｔａｐｈｙｌｏｃｏｃｃｕｓ ｓａｐｒｏｐｈｙｔｉｃｕｓ、 Ｓｔａｐｈｙｌｏｃｏｃｃｕｓ ｓｉｍｕｌａｎｓを播種したＬＢ寒天培地を設置した。実施例２と同様にして、ＬＥＤ光照射下で画像取得を行った。取得した５菌種の画像を対象として、ＭＡＴＬＡＢを用いて、最大コロニー増殖速度μｍａｘ、ヒストグラムの偏差Ｇ、ドーナツ性Ｄ、画像エントロピーＨ、画像エネルギー密度Ｅd、平均輝度値Ｉ、重み付け中心の差異Ｗ、中心領域の平均輝度Ｃを算出した。各菌種をグループとし、得られたパラメータに基づいた線形判別分析を行った。
算出した８パラメータに基づき、Ｓｔａｐｈｙｌｏｃｏｃｃｕｓ属５種について１菌種１０コロニーを対象とし線形判別分析を行った。その結果、グループ間の境界となる判別関数ＬＤ１〜ＬＤ４を得た。この判別関数をもとにＳｔａｐｈｙｌｏｃｏｃｃｕｓ属５種の判別分析を行ったところ、１００％の精度で判別可能となった（図７）。
【産業上の利用可能性】
【００５１】
本発明の方法によれば、安価、小型の装置でありながら、微生物を属レベル（例えば、エシェリヒア・コリ、スタフィロコッカス・アウレウス、シュードモナス エルギノーザ、サルモネラ・エンテリカ、及びカンジダ・アルビカンスの５属）を、１０時間以内に正確に判別することができる。また、種レベル（例えば、スタフィロコッカス属の近縁５種（スタフィロコッカス・アウレウス、スタフィロコッカス・エピデルミディス、スタフィロコッカス・ヘモリチカス、スタフィロコッカス・サプロフィティカス、スタフィロコッカス・シミュランス）を正確に判別可能である。本発明のデバイスは安価、小型の装置でありながら高精度に微生物を同定可能なことから、化粧品や食品、医薬品等の製造分野や、臨床診断における微生物検査に導入することにより、微生物同定試験を大幅に簡便、迅速化できると考えられ、その工学的意義は大きい。
【符号の説明】
【００５２】
１  培養容器
２  フォトセンサ
３  光源
【図１】
2018033430.tif 000007
【図２】
2018033430.tif 000008
【図３】
2018033430.tif 000009
【図４】
2018033430.tif 000010
【図５】
2018033430.tif 000011
【図６】
2018033430.tif 000012
【図７】
2018033430.tif 000013