(57)【要約】【課題】口腔内細菌叢からHDPの発症予測を実現するための予測装置及び予測方法を提供することである。【解決手段】本発明の一態様は、妊婦の口腔内の細菌叢に関する細菌叢データを取得するデータ取得機能と、前記細菌叢データに基づき妊娠高血圧症候群（HDP）予測モデルを生成するステップと、前記HDP予測モデルに基づき、前記取得した細菌叢データから前記妊婦のHDP発症予測を実行するHDP予測機能と、前記HDP発症予測の結果を通知する結果通知機能と、を有し、前記データ取得機能は、前記妊婦の妊娠前後の何れかの時期から出産前後の何れかの時期までのデータ取得期間において、前記細菌叢データを複数の時点にて取得し、前記HDP予測機能は、前記複数の時点において取得された前記細菌叢データに基づき前記HDP発症予測を実行する予測装置に関する。【選択図】図３
2019041754.tif 000002
【特許請求の範囲】
【請求項１】
妊婦の口腔内の細菌叢に関する細菌叢データを取得するデータ取得機能と、
前記細菌叢データに基づき妊娠高血圧症候群（HDP）予測モデルを生成するHDP予測モデル生成機能と、
前記HDP予測モデルに基づき、前記取得した細菌叢データから前記妊婦のHDP発症予測を実行するHDP予測機能と、
前記HDP発症予測の結果を通知する結果通知機能と、
を有し、
前記データ取得機能は、前記妊婦の妊娠前後の何れかの時期から出産前後の何れかの時期までのデータ取得期間において、前記細菌叢データを複数の時点にて取得し、
前記HDP予測機能は、前記複数の時点において取得された前記細菌叢データに基づき前記HDP発症予測を実行する予測装置。
【請求項２】
前記HDP予測機能は、前記複数の時点で取得された細菌叢データの各細菌種を独立した変数として設定し、HDP群と非HDP群とをよく分ける上位の所定数の変数を抽出し、
前記HDP予測機能は、forward stepwise selection方式に従って、前記抽出した変数から前記HDP予測モデルを生成する、請求項１記載の予測装置。
【請求項３】
前記HDP予測モデル生成機能は、HDP群と非HDP群とをよく分ける細菌種として、Actinomyces sp., Staphylococcus sp., Capnocytophaga sp., Prevotella sp., Leptotrichia sp., GN02 sp., Streptococcus sp., Haemophilus sp., Corynebacterium sp., Aggregatibacter sp., Rothia sp., Leptotrichia sp., Lautropia sp.の１つ以上を抽出する、請求項１又は２記載の予測装置。
【請求項４】
前記HDP予測モデル生成機能は、HDP発症有無予測モデル、HDP発症日予測モデル及びHDPリスクモデルを生成し、
前記HDP予測機能は、前記生成されたHDP発症有無予測モデル、HDP発症日予測モデル及びHDPリスクモデルにそれぞれ基づき、HDP発症有無予測、HDP発症日予測及びHDPリスク算出を実行する、請求項１乃至３何れか一項記載の予測装置。
【請求項５】
妊婦の口腔内の細菌叢に関する細菌叢データを取得するステップと、
前記細菌叢データに基づき妊娠高血圧症候群（HDP）予測モデルを生成するステップと、
前記HDP予測モデルに基づき、前記取得した細菌叢データから前記妊婦のHDP発症予測を実行するステップと、
前記HDP発症予測の結果を通知するステップと、
を有する予測方法であって、
前記取得するステップは、前記妊婦の妊娠前後の何れかの時期から出産前後の何れかの時期までのデータ取得期間において、前記細菌叢データを複数の時点にて取得し、
前記実行するステップは、前記複数の時点において取得された前記細菌叢データに基づき前記HDP発症予測を実行する予測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、一般に医療診断技術に関し、より詳細には、非侵襲的に妊娠高血圧症候群の発症予測又は異常検出を行うための妊娠高血圧症候群発症有無予測のための予測装置及び予測方法に関する。
【背景技術】
【０００２】
妊娠合併症として、妊娠高血圧症候群（HDP: Hypertensive Disorders of Pregnancy）は、母体死亡、胎児発育不全、胎児死亡、常位胎盤早期剥離などを引き起こす母児の命に関わる周産期医療上の代表的な疾患であり、複雑な相互作用により発症する多因子疾患である。また、HDP罹患妊婦はその後の生活習慣病の発症リスクが高いことが知られている。妊婦の口腔内細菌叢と妊娠合併症との関連が報告されているが、その詳細なメカニズムは不明な点が多い。
【０００３】
非特許文献１によると、妊婦124名に対し、妊娠期間中あるいは産後3日以内にClinical attachment level (CAL)、Probing depth (PD)を測定し、PLBW(Preterm Low Birth Weight)だったケース群とそうでないコントロール群で評価したところ、ケース群（31名）のCAL平均値がコントロール群よりも有意に大きかった。また、CALが3mm以上の歯が60%以上を占めた妊婦においては、PLBW発症危険率が全体で5.9倍、初産で6.7倍だった。
【０００４】
また、非特許文献２によると、切迫早産から早産となった18名とそうでない正期産40名の第2トライメスターの歯垢細菌叢を比較したところ、歯周病関連細菌の一つであるTannerella forsythenesisの割合が有意に高かった。このTannerella forsythenesisは、歯周病に関連のある細菌と言われる３つのred complexのうちの一つである。
【０００５】
また、非特許文献３によると、腸内や膣内細菌叢は妊娠期間中大きく変動することが知られており、食事や歯磨きなどの影響を受ける口腔内細菌叢はさらに変動が激しいことが想定される。
【先行技術文献】
【非特許文献】
【０００６】
【非特許文献１】Offenbacher S, et al, Periodontal Infection as a Possible Risk Factor for Preterm Low Birth Weight, J Peridontol, 67, 1103-1113, 1996
【非特許文献２】Hasegawa K, et al, Associations between systemic status, oeruidibtak statysm seryn cytokine levels and delivery outcomes in pregnant women with a diagnosis of threatened premature labor (TPL), J Peridontol, 74, 1764-1770, 2003
【非特許文献３】Daniel B. DiGiulio, et al, Temporal and spatial variation of the human microbiota during pregnancy, PNAS, 112, 35, 11060-11065, 2015
【発明の概要】
【発明が解決しようとする課題】
【０００７】
このように、口腔内細菌叢と妊娠合併症との関連が報告されているが、口腔内細菌叢と妊娠合併症の1つであるHDPとの関連についての報告はほとんどない。また、口腔内細菌叢と妊娠合併症との関係には人種差があると考えられるが、日本人の妊婦を対象にした研究が少ない。また、口腔内細菌叢の経時的な変化による妊娠合併症への影響を考慮できていない。さらに、口腔内細菌叢から妊娠合併症を予測する際、個人内での細菌叢の揺らぎにより予測をうまく行えない可能性がある。
【０００８】
上述した問題点を鑑み、本発明の課題は、口腔内細菌叢からHDPの発症予測を実現するための予測装置及び予測方法を提供することである。
【課題を解決するための手段】
【０００９】
上記課題を解決するため、本発明の一態様は、妊婦の口腔内の細菌叢に関する細菌叢データを取得するデータ取得機能と、前記細菌叢データに基づき妊娠高血圧症候群（HDP）予測モデルを生成するHDP予測モデル生成機能と、前記HDP予測モデルに基づき、前記取得した細菌叢データから前記妊婦のHDP発症予測を実行するHDP予測機能と、前記HDP発症予測の結果を通知する結果通知機能と、を有し、前記データ取得機能は、前記妊婦の妊娠前後の何れかの時期から出産前後の何れかの時期までのデータ取得期間において、前記細菌叢データを複数の時点にて取得し、前記HDP予測機能は、前記複数の時点において取得された前記細菌叢データに基づき前記HDP発症予測を実行する予測装置に関する。
【発明の効果】
【００１０】
本発明によると、口腔内細菌叢からHDPの発症予測を実現するための予測装置及び予測方法を提供することができる。
【図面の簡単な説明】
【００１１】
【図１】図１は、本発明の一実施例による予測システムの構成を示す概略図である。
【図２】図２は、本発明の一実施例による予測装置の機能構成を示すブロック図である。
【図３】図３は、本発明の一実施例による予測処理を示すフローチャートである。
【図４】図４は、本発明の一実施例による特徴量抽出・モデル作成処理を示すフローチャートである。
【図５】図５は、本発明の他の実施例による特徴量抽出・モデル作成処理を示すフローチャートである。
【図６】図６は、本発明の他の実施例による特徴量抽出・モデル作成処理を示すフローチャートである。
【発明を実施するための形態】
【００１２】
以下、図面に基づいて本発明の実施の形態を説明する。
【００１３】
後述される実施例では、HDP発症予測のための予測装置が開示される。後述される実施例によると、複数の妊娠フェーズにおいて採取された妊婦の口腔内の細菌叢に関する細菌叢データを取得すると、予測装置は、HDPと相関を有する細菌種を抽出し、抽出した細菌種を用いて生成したHDP予測モデルに基づき、取得した細菌叢データから妊婦のHDP発症予測を実行する。これにより、人種間の差があると共に、経時的に変化すると考えられている細菌叢データからのHDP予測に対して、特定の人種に適したHDP予測モデルが生成可能になる。
【００１４】
まず、図１を参照して、本発明の一実施例による予測システムを説明する。図１は、本発明の一実施例による予測システムの構成を示す概略図である。
【００１５】
図１に示されるように、予測システム１０は、細菌叢計測装置５０及び予測装置１００を有する。
【００１６】
細菌叢計測装置５０は、採取された歯垢検体から口腔内の細菌叢を測定し、生成された細菌叢データを予測装置１００に提供する。具体的には、検査技師が妊婦の歯垢検体を採取し、細菌叢計測装置５０を用いて採取した歯垢検体に存在する細菌種を示す細菌叢データを取得し、取得した細菌叢データを予測装置１００に提供する。
【００１７】
人間の口腔内に存在する細菌種として、例えば、Porphyromonas gingivalis, Treponema denticola, Tabberella forsythensis, Bergeyella, Kingella denitrificans, Scardovia wiggsiae, Streptococus mutans, Selenomonas, Capnocytophaga, Treponema vincentiiなどを含む200種類以上の細菌種が知られており、特にPorphyromonas gingivalis, Treponema denticola, Tabberella forsythensisは、歯周病に関わるred complexとして分類され、妊娠合併症との相関が指摘されている。細菌叢計測装置５０は、口腔内の細菌種を計測するのに用いられる何れか既知のものであってもよい。
【００１８】
予測装置１００は、提供された細菌叢データからHDPを予測する。具体的には、予測装置１００は、以下で詳細に説明されるように、HDP群の細菌叢データ及び非HDP群の細菌叢データを含む既知のデータからHDP予測モデルを作成し、作成したHDP予測モデルに予測対象となる妊婦の細菌叢データを入力することによって、当該妊婦のHDPの発症の有無を予測する。なお、本発明によるHDP発症予測は、HDPの発症の有無に限定されず、後述されるように、HDP発症日の予測、HDPリスクの算出などを含むものであってもよい。
【００１９】
ここで、予測装置１００は、典型的には、コンピュータ、サーバ、スマートフォン、タブレットなどの通信機能を備えた情報処理装置により実現されてもよい。例えば、情報処理装置に搭載されたプロセッサが、メモリ装置に格納されたデータやプログラムを処理及び実行することによって、後述される各種機能及び処理を実行する。しかしながら、予測装置１００は、何れか特定のハードウェア構成に限定されるものでなく、適切なハードウェア構成により実現されてもよい。
【００２０】
次に、図２を参照して、本発明の一実施例による予測装置を説明する。図２は、本発明の一実施例による予測装置の機能構成を示すブロック図である。
【００２１】
図２に示されるように、予測装置１００は、データ取得機能１１０、データ保存機能１２０、HDP予測モデル生成機能１３０、HDP予測機能１４０及び結果通知機能１５０を有する。
【００２２】
データ取得機能１１０は、細菌叢計測装置５０から細菌叢データを取得する。具体的には、データ取得機能１１０は、妊婦の妊娠前後の何れかの時期から出産前後の何れかの時期までのデータ取得期間において、妊婦に関する細菌叢データを複数の時点にて取得する。例えば、データ取得機能１１０は、細菌叢計測装置５０から、異なる妊娠フェーズ（例えば、１２〜２３週と２４〜３４週）に採取された妊婦の口腔内の細菌叢データを取得してもよい。
【００２３】
データ保存機能１２０は、取得した妊婦の口腔内の細菌叢データを保存する。具体的には、データ保存機能１２０は、取得した細菌叢データを各妊娠フェーズで採取された細菌種（例えば、２００菌種など）のそれぞれを独立した変数として保存する。例えば、２つの妊娠フェーズでそれぞれ２００菌種の細菌種が採取された場合、経時的な複数時点の歯垢細菌組成の変化を加味したモデルを作成できるように、４００個（２００菌種×２妊娠フェーズ）の変数が設定されてもよい。
【００２４】
HDP予測モデル生成機能１３０は、取得された妊婦の口腔内の細菌叢データに基づき、HDPの発症の有無を予測するHDP予測モデルを生成する。具体的には、HDP予測モデル生成機能１３０は、HDPを発症しなかった妊婦（非HDP群）とHDPを発症した妊婦（HDP群）とを含む複数の妊婦から取得した口腔内の細菌叢データに基づき、HDP発症予測モデルを生成及び更新してもよい。例えば、以下で詳細に説明されるように、HDP予測モデル生成機能１３０は、取得した細菌叢データからHDP群と非HDP群とをよく分ける細菌種を特定し、変数を削減するため、二群をよく分ける細菌種として１０〜２０個の細菌種×妊娠フェーズを抽出してもよい。HDP群と非HDP群とをよく分ける細菌種として、red complexとして分類されるPorphyromonas gingivalis, Treponema denticola, Tabberella forsythensisが抽出されてもよい。また、HDP予測モデル生成機能１３０は、HDP群と非HDP群とをよく分ける細菌種として、Actinomyces sp., Staphylococcus sp., Capnocytophaga sp., Prevotella sp., Leptotrichia sp., GN02 sp., Streptococcus sp., Haemophilus sp., Corynebacterium sp., Aggregatibacter sp., Rothia sp., Leptotrichia sp., Lautropia sp.の１つ以上を抽出してもよい。
【００２５】
続いて、HDP予測モデル生成機能１３０は、抽出された細菌種を用いてHDP予測モデルを作成する。具体的な予測モデルの作成手法として、例えば、HDP予測モデル生成機能１３０は、Linear SVMを用いて、特徴量として用いられる細菌種をforward stepwise selectionによって選択してもよい。例えば、最初に２０個の変数（細菌種）について、HDP予測モデル生成機能１３０は、１つの変数だけを用いてHDP予測モデルM(1)〜M(20)を作成し、最も性能の良いM(A)（変数Aを使用）を特定する。次に、HDP予測モデル生成機能１３０は、当該変数を除いた残りの１９個の変数から１つの変数を選択してAと併せた２変数の予測モデルM(A,1)〜M(A,19)を作成し、M(A)を含めて最適なモデルを選択する。HDP予測モデル生成機能１３０は、上述した一連の処理を再帰的に実行し、最終的に最も性能の良い予測モデルM(A,B,C,...)を特定する。
【００２６】
なお、上述した実施例では、HDPの発症有無を示すHDP発症予測モデルについて説明したが、本発明によるHDP予測モデルは、これに限定されず、HDPの発症日を予測するHDP発症日予測及びHDPの発症リスクモデル又は発症確率を予測するHDP発症リスクモデルであってもよい。
【００２７】
HDP予測機能１４０は、HDP予測モデルに基づき、取得した妊婦の細菌叢データから当該妊婦のHDP発症予測を実行する。具体的には、HDP予測機能１４０は、複数の時点において取得した妊婦の細菌叢データをHDP予測モデルに投入し、複数の時点において取得した細菌叢データに基づくHDP予測モデルの出力結果に基づき当該妊婦のHDP発症有無予測、HDP発症日予測又はHDP発症リスク算出を実行してもよい。すなわち、HDP予測機能１４０は、当該妊婦にHDP発症の可能性があるか予測するか（HDP発症有無予測）、当該妊婦のHDP発症日を予測するか（HDP発症日予測）、あるいは、当該妊婦のHDP発症リスクを算出してもよい（HDP発症リスク算出）。
【００２８】
結果通知機能１５０は、HDP予測機能１４０によって実行されたHDP発症予測の結果を妊婦、医療関係者などに通知する。
【００２９】
次に、図３〜４を参照して、本発明の一実施例による予測処理を説明する。当該予測処理は、上述した予測装置１００によって実行される。図３は、本発明の一実施例による予測処理を示すフローチャートである。
【００３０】
図３に示されるように、ステップＳ１００において、予測装置１００は、細菌叢計測装置５０から細菌叢データを取得し、取得した細菌叢データを保存する。具体的には、細菌叢計測装置５０は、妊婦から採取された歯垢検体を計測し、当該歯垢検体に含まれる細菌種量を示す細菌叢データを予測装置１００に提供する。予測装置１００は、細菌叢計測装置５０から取得した細菌叢データを妊婦毎に保存する。本実施例によると、細菌叢データは、例えば、妊娠フェーズ毎など、妊娠前後の何れかの時期から出産前後の何れかの時期までのデータ取得期間において各妊婦に対して複数回採取される。
【００３１】
ステップｓ１０１において、予測装置１００は、取得した細菌叢データに対して異常値を検出及び除去し、欠測値を補完する。例えば、基準値から閾値範囲外の細菌種量が測定された場合、予測装置１００は、当該測定値を実際の値ではなく計測エラーとみなして、当該測定値を除去すると共に、他の細菌種量の値が類似する別の歯垢検体の当該測定値にランダムなノイズを加えて流用するか、他の細菌種量の値が類似する別の歯垢検体を複数抽出してそれらの当該測定値の分布からその分布に従うように１つの推定値を算出するなど何れかの方法により当該測定値を補完してもよい。例えば、検体Sの細菌種spA-2が閾値範囲外または欠測値である場合、まず検体S以外の検体群においてspA-2の細菌種量の分散Vを計算する。つぎに、細菌種群BALL={spA-1, spB-1,…, spA-2, spB-2}からspA-2を抜いた細菌種群BSUBの各細菌種量が検体Sと近い検体を1検体抽出し、そのspA-2の細菌種量に正規分布N(0,V)に従うランダムな１点を足し合わせた値を検体SのspA-2の推定量とする方法がある。このとき、検体Sと近い検体の抽出には、各細菌種量を各細菌種量の標準偏差で割って標準化した値を用いてユーグリッド距離を計算しても良いし、MSE（平均平方誤差）又はR2（決定係数）などを用いても良い。また別の補完方法として、細菌種群BSUBの各細菌種量が検体Sと近い検体を例えば10検体だけ抽出し、その平均値E2および分散V2を求め、正規分布N(E2,V2)に従うランダムな１点を検体SのspA-2の推定量とする方法がある。また、分布は必ずしも正規分布とせず、任意の分布で推定してもよい。
【００３２】
ステップｓ１０２において、予測装置１００は、既に出産して出産状況が確定している妊婦のデータを用いて、非HDP群とHDP群とをよく分ける変数を特徴量として抽出する。これには、ステューデントのt検定、ウェルチのt検定を用いた方法、Lasso回帰、Ridge回帰、Elastic Netなどにより有用な変数を限定する方法、主成分分析の寄与率などから変数を限定する方法などがあるが、それらに限定されるものでない。
【００３３】
ステップｓ１０３、ｓ１０５及びｓ１０７において、予測装置１００は、抽出した特徴量を用いてHDP発症予測モデルを生成する。具体的には、予測装置１００は、ステップｓ１０３において、HDPになるかどうかを予測する分類モデルとしてHDP発症予測モデルを生成する。また、予測装置１００は、ステップｓ１０５において、各妊婦のHDP発症日を予測する回帰モデルとしてHDP発症日予測モデルを生成する。また、予測装置１００は、ｓ１０７において、HDP発症のリスクを算出するHDP発症リスク算出モデルを生成する。回帰モデルには、線形回帰モデル、ロジスティック回帰モデル、SVMによる回帰モデル、多層ニューラルネットワークによる回帰モデルなどがあるが、それらに限定されるものでない。また、分類モデルにおいても、線形分類、ロジスティック回帰、ベイジアンネットワーク、SVM、k近傍法、ランダムフォレスト、多層ニューラルネットワークによる分類モデルなどがあるが、それらに限定されるものでない。HDP発症リスク算出モデルについては、オッズ比、寄与危険度、相対危険度などがあるが、それらに限定されるものでない。また、ステップｓ１０３、ｓ１０５、ｓ１０７は、ｓ１０２と切り離さずに一体で処理されてもよい。
【００３４】
ステップｓ１０４、ｓ１０６及びｓ１０８において、予測装置１００は、生成されたHDP発症予測モデル、HDP発症日予測モデル及びHDP発症リスク算出モデルに予測対象の妊婦の細菌叢データをそれぞれ投入し、HDP発症予測を実行する。具体的には、予測装置１００は、ステップｓ１０４において、HDP発症予測モデルを利用して予測対象の妊婦の細菌叢データから将来のHDP発症の有無を予測する。また、予測装置１００は、ステップｓ１０６において、HDP発症日予測モデルを利用して予測対象の妊婦の細菌叢データからHDP発症日を予測する。また、予測装置１００は、ステップｓ１０８において、HDP発症リスク算出モデルを利用して予測対象の妊婦のデータからHDPリスクを算出する。なお、HDP有無予測は、最終的なHDP発症有無の予測、この先n1週間（n1は任意の数値）でHDPを発症しない予測、この先n2週間（n2は任意の数値）でHDPを発症する予測、などを含むが、それらに限定されるものでない。
【００３５】
ステップｓ１０９において、予測装置１００は、HDP発症予測結果、HDP発症日予測結果及び／又はHDP発症リスク計算結果を妊婦又は医療機関に通知し、当該予測処理を終了する。
【００３６】
図４は、本発明の一実施例による特徴量抽出・モデル作成処理を示すフローチャートである。当該特徴量抽出・モデル作成処理は、上述された図３のステップｓ１０２，ｓ１０３において実行可能であり、後述されるステップｓ１０２−１〜５がステップｓ１０２の特徴量の抽出処理に対応し、ステップｓ１０３−１〜２がステップｓ１０３のHDP有無予測モデルの生成処理に対応する。当該特徴量抽出・モデル作成処理では、HDP予測モデル生成機能１３０は、取得したデータをHDP群と非HDP群とに分け、各変数のうち二群をよく分ける（すなわち、区別する、判別する、分離するなど）変数、例えば、二群をよく分ける上位の所定数の変数を抽出し、抽出された変数を用いてHDP予測モデルを作成し、最も性能の良い変数によるHDP予測モデルを特定する。次に、HDP予測モデル生成機能１３０は、当該変数を除いた残りの変数と特定された変数とを併せた２変数の予測モデルを作成し、これらのモデルを含めて最適なモデルを選択する。HDP予測モデル生成機能１３０は、上述した一連の処理を再帰的に実行し、最終的にHDP予測モデルを特定する。なお、上記はHDP発症予測モデルに利用する特徴量の抽出方法やHDP発症予測モデルの選定方法の一例であって、HDP発症予測モデルはこれに制限されない。例えば、多層ニューラルネットワークまたはLasso回帰などによって特徴量を抽出し、HDP発症予測モデルを生成する方法もある。
【００３７】
ステップｓ１０２−２において、予測装置１００は、モデル生成のために、既に出産済みの妊婦の細菌叢データを抽出する。
【００３８】
ステップｓ１０２−３において、予測装置１００は、妊娠フェーズ毎に採取された歯垢検体の細菌叢データについて、各妊娠フェーズにおいて採取された各細菌種のそれぞれを独立した変数として保存する。例えば、２つの妊娠フェーズでそれぞれ２００菌種の細菌種が採取された場合、４００個（２００菌種×２フェーズ）の変数が設定される。例えば、４００個のデータは、当該妊婦のＩＤ（userid）と関連付けて、"spA-1, spB-1, ...,spA-2, spB-2, ..."の変数列により表されてもよい。
【００３９】
ステップＳ１０２−４において、予測装置１００は、HDP群と非HDP群とに細菌叢データを分け、検定にてp値がa未満（例えば、a=0.01）となる二群間をよく分ける細菌種を抽出する。
【００４０】
ステップＳ１０３−１において、予測装置１００は、細菌叢データをHDP群D1と非HDP群D2とに分け、それぞれの比率e（例えば、e=0.9）でHDP群D1を学習用データD1a及び評価用データD1bに分け、非HDP群D2を学習用データD2a及び評価用データD2bに分ける。
【００４１】
ステップＳ１０３−２において、予測装置１００は、学習用データD1a, D2aを用いてforward stepwise selectionによって特徴量を選択し、HDP予測モデルMを作成する。例えば、最初に２０個の特徴量（細菌種）について、予測装置１００は、１つの変数だけを用いてHDP予測モデルM(1)〜M(20)を作成する。そして、予測装置１００は、評価データ用D1b, D2bを用いて、最も性能の良いM(A)（変数Aを使用）を特定する。性能スコアとして、例えば、AUC又はf-scoreが利用されてもよい。次に、予測装置１００は、学習用データD1a, D2aを用いて、当該変数を除いた残りの１９個の変数から１つの変数を選択してAと併せた２変数のHDP予測モデルM(A,1)〜M(A,19)を作成し、評価用データD1b, D2bを用いて、M(A)を含めて最適なモデルを選択する。予測装置１００は、上述した一連の処理を再帰的に実行し、最終的にHDP予測モデルM(A,B,C,...)を特定する。
【００４２】
なお、上述した実施例では、ＨＤＰ発症有無予測モデルを生成する処理を説明したが、同様の処理によってＨＤＰ発症日予測モデル及びＨＤＰ発症リスク算出モデルが生成可能であることは明らかであろう。なお、HDP予測モデルMは、SVM、線形ロジスティック回帰、ニューラルネットワーク、ランダムフォレストなどに基づき生成されてもよい。また、性能スコアは、AUC又は及びf-scoreに限定されず、sensitivityとspecificityとの調和平均などが用いられてもよい。また、forward step selectionに限定されず、backward stepwise selection又はbidirectional stepwise selectionが利用されてもよい。また、全細菌種を利用し、PCAやMDS、NMDSなどで次元削減した後、PC1, PC2, ..., PCNのN個の特徴量に変換した後対して、これらにforward/backward/bidirectional stepwise selectionを適用してもよい。
【００４３】
図５は、本発明の他の実施例による特徴量抽出・モデル作成処理を示すフローチャートである。当該特徴量抽出・モデル作成処理は、上述された図３のステップｓ１０２，ｓ１０５において実行可能であり、後述されるステップｓ１０２−１〜５がステップｓ１０２の特徴量の抽出処理に対応し、ステップｓ１０５−１〜３がステップｓ１０５の出産日予測モデルの生成処理に対応する。なお、ステップＳ１０２−１〜５は図４のステップＳ１０２−１〜５と同じであるため、その説明は省く。
【００４４】
ステップＳ１０５−１において、予測装置１００は、細菌叢データをHDP群D1と非HDP群D2とに分け、それぞれの比率e（例えば、e=0.9）でHDP群D1を学習用データD1a及び評価用データD1bに分け、非HDP群D2を学習用データD2a及び評価用データD2bに分ける。
【００４５】
ステップＳ１０５−２において、予測装置１００は、予測性能を向上させるため、学習用データD1a, D2aにおいてD1aをオーバーサンプリングしてもよく、及び／又はD2aをアンダーサンプリングしてもよい。
【００４６】
ステップＳ１０５−３において、予測装置１００は、学習用データD1a, D2aを用いてforward stepwise selectionによって特徴量を選択し、HDP発症日予測モデルMを作成する。例えば、最初に２０個の変数（細菌種）について、予測装置１００は、１つの変数だけを用いてHDP発症日予測モデルM(1)〜M(20)を作成する。そして、予測装置１００は、評価データ用D1b, D2bを用いて、最も性能の良いM(A)（変数Aを使用）を特定する。性能スコアとして、例えば、MSE（平均平方誤差）又はR2が利用されてもよい。次に、予測装置１００は、学習用データD1a, D2aを用いて、当該変数を除いた残りの１９個の変数から１つの変数を選択してAと併せた２変数のHDP発症日予測モデルM(A,1)〜M(A,19)を作成し、評価用データD1b, D2bを用いて、M(A)を含めて最適なモデルを選択する。予測装置１００は、上述した一連の処理を再帰的に実行し、最終的にHDP発症日予測モデルM(A,B,C,...)を特定する。
【００４７】
なお、HDP発症日予測モデルMは、SVM、線形ロジスティック回帰、ニューラルネットワークなどに基づき生成されてもよい。また、性能スコアは、MSE（平均平方誤差）又はR2に限定されなくてもよい。また、forward step selectionに限定されず、backward stepwise selection又はbidirectional stepwise selectionが利用されてもよい。また、全細菌種を利用し、PCAやMDS、NMDSなどで次元削減した後、PC1, PC2, ..., PCNのN個の特徴量に対してforward/backward/bidirectional stepwise selectionを適用してもよい。
【００４８】
図６は、本発明の他の実施例による特徴量抽出・モデル作成処理を示すフローチャートである。当該特徴量抽出・モデル作成処理は、上述された図３のステップｓ１０２，ｓ１０７において実行可能であり、後述されるステップｓ１０２−１〜５がステップｓ１０２の特徴量の抽出処理に対応し、ステップｓ１０７−１〜２がステップｓ１０７のHDPリスクモデルの生成処理に対応する。なお、ステップＳ１０２−１〜５は図４のステップＳ１０２−１〜５と同じであるため、その説明は省く。
【００４９】
ステップＳ１０７−１において、予測装置１００は、細菌種（spA-1, spB-1, ..., spA-2, spB-2, ...）毎に量の閾値範囲SspA-1=[AspA-1, Bspa-1]を設定し、HDPに対するリスク比(RRspA-1)を計算する。ここで、リスク比はRR=a(c+d)/c(a+b)であり、aはHDP群であって（細菌種の量が）閾値範囲内である歯垢検体数であり、bは非HDP群であって閾値範囲内である歯垢検体数であり、cはHDP群であって閾値範囲外である歯垢検体数であり、dは非HDP群であって閾値範囲外である歯垢検体数である。
【００５０】
ステップＳ１０７−２において、予測装置１００は、予測対象の妊婦の細菌種の量を（xspA-1, xspB-1,...）とした場合、
RR=ΠspRRspf(xsp)
f(xsp) = 1 if xsp is in Ssp
f(xsp) = -1 if xsp is not in Ssp
によって、HDPリスク予測モデルを生成してもよい。
【００５１】
なお、forward step selectionに限定されず、backward stepwise selection又はbidirectional stepwise selectionが利用されてもよい。また、全細菌種を利用し、PCAやMDS、NMDSなどで次元削減した後、PC1, PC2, ..., PCNのN個の特徴量に対してforward/backward/bidirectional stepwise selectionを適用してもよい。
上記の一例を実施した結果では、HDP群と非HDP群をよく分ける細菌種として、Actinomyces sp., Staphylococcus sp., Prevotella sp., GN02 sp. Staphylococcus sp.などが抽出され、これらを特徴量として用いたHDP発症有無予測モデルの性能は、感度0.60、特異度0.91となった。
【００５２】
従来は日本人の歯垢細菌叢からHDPを予測するモデルはなかったが、本発明によると、時系列に変化する細菌叢からHDPに関連する複数の細菌種を組み合わせることによって、日本人に適したHDP予測モデルを作成することが可能となり、HDPの事前の予防及び迅速な対応を可能にする。
【００５３】
以上、本発明の実施例について詳述したが、本発明は上述した特定の実施形態に限定されるものではなく、特許請求の範囲に記載された本発明の要旨の範囲内において、種々の変形・変更が可能である。
【符号の説明】
【００５４】
１０  予測システム
５０  細菌叢計測装置
１００  予測装置
１１０  データ取得機能
１２０  データ保存機能
１３０  HDP予測モデル生成機能
１４０  HDP予測機能
１５０  結果通知機能
【図１】
2019041754.tif 000003
【図２】
2019041754.tif 000004
【図３】
2019041754.tif 000005
【図４】
2019041754.tif 000006
【図５】
2019041754.tif 000007
【図６】
2019041754.tif 000008