JP7224297B2 - 歯周ポケット炎症面積の推定方法 - Google Patents

歯周ポケット炎症面積の推定方法 Download PDF

Info

Publication number
JP7224297B2
JP7224297B2 JP2019550508A JP2019550508A JP7224297B2 JP 7224297 B2 JP7224297 B2 JP 7224297B2 JP 2019550508 A JP2019550508 A JP 2019550508A JP 2019550508 A JP2019550508 A JP 2019550508A JP 7224297 B2 JP7224297 B2 JP 7224297B2
Authority
JP
Japan
Prior art keywords
bacteria
streptococcus
pisa
amount
treponema
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2019550508A
Other languages
English (en)
Other versions
JPWO2019088271A1 (ja
Inventor
直之 外川
あい 原
伸也 村上
剛徳 野崎
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
GC Corp
Original Assignee
GC Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by GC Corp filed Critical GC Corp
Publication of JPWO2019088271A1 publication Critical patent/JPWO2019088271A1/ja
Application granted granted Critical
Publication of JP7224297B2 publication Critical patent/JP7224297B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6813Hybridisation assays
    • C12Q1/6834Enzymatic or biochemical coupling of nucleic acids to a solid phase
    • C12Q1/6837Enzymatic or biochemical coupling of nucleic acids to a solid phase using probe arrays or probe chips
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6876Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
    • C12Q1/6883Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B20/00ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
    • G16B20/20Allele or variant detection, e.g. single nucleotide polymorphism [SNP] detection
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B40/00ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
    • G16B40/20Supervised data analysis
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12NMICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N15/00Mutation or genetic engineering; DNA or RNA concerning genetic engineering, vectors, e.g. plasmids, or their isolation, preparation or purification; Use of hosts therefor
    • C12N15/09Recombinant DNA-technology
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6876Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
    • C12Q1/6888Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms
    • C12Q1/689Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms for bacteria
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B30/00ICT specially adapted for sequence analysis involving nucleotides or amino acids
    • G16B30/10Sequence alignment; Homology search
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H10/00ICT specially adapted for the handling or processing of patient-related medical or healthcare data
    • G16H10/40ICT specially adapted for the handling or processing of patient-related medical or healthcare data for data related to laboratory analysis, e.g. patient specimen analysis
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A50/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE in human health protection, e.g. against extreme weather
    • Y02A50/30Against vector-borne diseases, e.g. mosquito-borne, fly-borne, tick-borne or waterborne diseases whose impact is exacerbated by climate change

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Organic Chemistry (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Analytical Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Biotechnology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Genetics & Genomics (AREA)
  • Molecular Biology (AREA)
  • General Engineering & Computer Science (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Microbiology (AREA)
  • Medical Informatics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Bioethics (AREA)
  • Artificial Intelligence (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Evolutionary Computation (AREA)
  • Epidemiology (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
  • Investigating Or Analysing Biological Materials (AREA)

Description

本発明は、歯周ポケット炎症面積を推定する方法に関する。
歯周病の診断は歯周ポケットの測定、アタッチメントレベル、X線画像診断等によって行われている。しかしながら、これらの診断方法は被験体への負担が大きく、特に多人数について実施する場合は相当の時間を要する。さらに、これらの歯周病診断方法は操作手順が複雑であり、歯科医師の経験や技能に基づくため判定基準に個人差があるなどの問題がある。
そこで、これまでに、簡便な歯周病の診断方法が提案されている。例えば、特許文献1では、歯肉溝滲出液に含まれるタンパク質を歯周病疾患マーカーとして利用した歯周病診断方法が開示されている。また、特許文献2では、唾液中の複数のタンパク質を分析することにより、歯周ポケットの深さ(PPD)や歯肉出血指数(BOP)を予測する手法が開示されている。しかしながら、通常、PPDやBOPについては、測定箇所が168箇所(28本の歯×6点法)も存在するため、実際にどの箇所について予測がされているのか不明である。
他方、より一般的には、特許文献3に開示されているように、歯肉溝浸出液に存在する歯周病細菌の測定や、唾液中の歯周病細菌の測定が実施されている。唾液中に存在する歯周病細菌に対応する指標としては、歯周ポケットの深さ(PPD)と歯肉出血指数(BOP)から算出される、炎症面積(PISA値)が提唱されている(非特許文献1)。PISA値は、歯周ポケットの深さ(PPD)から計算されたPESA値(歯周ポケット内表面積)に対して歯肉出血指数(BOP)が乗算されることによって算出される。
また、これとは別に滝沢らからはCAPRS値が提唱されている(非特許文献2)。本指標は文献記載のとおり、アタッチメントレベルから臨床歯根表面積(歯肉縁より根尖側に位置する歯根の表面積の総和;Clinical Area of ToothRoot Surface; CARS)を求めた後、その値から有効歯根表面積の値を引いて算出される。実際には、歯肉縁の位置は全て解剖的歯頚線と一致しているケースに置き換えて算出され、これは非特許文献1のFig1.(b)に示されるPESAの値と同一であると考えられる。CAPRS値は、歯肉出血指数(BOP)が計算に加味されていないが、文献記載のとおり面積そのものをポケット内面の炎症領域と近似的に捉え、PISA値とともに「歯周ポケット炎症面積」と考えることができる。

これまで、これら「歯周ポケット炎症面積」を簡便に評価するために、唾液中の歯周病細菌数との関係について検討されたことがあるが、P.g菌(Porphyromonas gingivalis(ポルフィロモナス・ジンジバリス))の菌数とCAPRS値との間、及び、レッドコンプレックス菌(P.g菌、T.d菌(Treponema denticola(トレポネーマ・デンティコラ))、T.f菌(Tannerella forsythensis(タネレラ・フォーサイセンシス)))の菌数とCAPRS値との間には相関が無く(非特許文献2、図3a,b)、唾液サンプルから簡便に「歯周ポケット炎症面積」予測する方法は知られていなかった。
特開2004-51536号公報 特許第5869323号 特許第4917815号
Periodontal inflamed surface area: quantifying inflammatory burden.J Clin Periodontol. 2008 Aug;35(8):668-73. 「歯周病の全身疾患関連検査マーカーとしての歯周ポケット面積評価法の臨床的意義」日本口腔検査学会雑誌第1巻第1号:7-12、2009
このような状況下において、唾液中の細菌量の検出結果に基づいて、歯周ポケット炎症面積を推定する方法の提供が望まれていた。
また、このような状況下において、唾液中の細菌量の検出結果に基づいて、炎症面積(PISA値)等といった歯周組織の炎症度を簡便に予測する方法の提供が望まれていた。
本発明は、上記状況を考慮してなされたもので、以下に示す、歯周ポケット炎症面積を推定する方法、及び歯周組織の炎症度を包括的に推定する方法を提供するものである。
[1]唾液中の2種以上の細菌の細菌量を検出し、得られた検出結果を指標として歯周ポケット炎症面積を推定する方法であって、検出される細菌は、該細菌の細菌量と歯周ポケット炎症面積とが正の相関関係を示す細菌と、
該細菌の細菌量と歯周ポケット炎症面積とが負の相関関係を示す細菌とを含む、方法。
[2] 前記歯周ポケット炎症面積がPISA、又はCAPRSの値で表されるものである、[1]に記載の方法。
[3] 正の相関関係を示す細菌が、Treponema denticola、Tannerella forsythia、Fusobacterium nucleatum subsp. animalis、Porphyromonas gingivalis、Campylobacter rectus、Fusobacterium nucleatum subsp. nucleatum、Selenomonas noxia、Veillonella parvula、Streptococcus gordonii、Fusobacterium nucleatum subsp. vincentii、Streptococcus intermedius、Capnocytophaga ochracea、Capnocytophaga sputigena、Aggregatibacter actinomycetemcomitans、Fusobacterium nucleatum subsp. polymorphum、Fusobacterium periodonticum、SR1 sp. OT 345、Porphyromonas catoniae、Selenomonas sputigena、Neisseria flavescens、Streptococcus sobrinus、Parvimonas micra、Peptostreptococcus stomatis、Treponema socranskii、Eubacterium saphenum、Eubacterium nodatum、Treponema medium、Filifactor alocis、及びPorphyromonas endodontalisからなる群から選ばれる少なくとも1種である、[1]又は[2]に記載の方法。
[4] 負の相関関係を示す細菌が、Streptococcus mutans、Actinomyces odontolyticus、Streptococcus mitis bv 2、Streptococcus mitis、Campylobacter concisus、Capnocytophaga gingivalis、Prevotella pallens、Streptococcus salivarius、Eubacterium sulci、Rothia mucilaginosa、Prevotella denticola、Veillonella atypica、Prevotella histicola、Megasphaera micronuciformis、及びStreptococcus parasanguinisからなる群から選ばれる少なくとも1種である、[1]~[3]のいずれか1項に記載の方法。
[5] 以下の(1)~(4)工程を有する、[1]~[4]のいずれか1項に記載の方法。

(1)歯周ポケット炎症面積が既知の被験者の唾液サンプルから、唾液中の各種細菌の細菌量を検出する工程
(2)当該各種細菌の細菌量について、菌ごとに固有の歯周ポケット炎症面積との相関係数を求め、各種細菌の細菌量と歯周ポケット炎症面積との関係式を構築して予測モデルを作成する工程
(3)歯周ポケット炎症面積が未知の被験者の唾液サンプルから、唾液中の各種細菌の細菌量を検出する工程
(4)(3)で得られた各種細菌量を(2)で得られた関係式に入れて歯周ポケット炎症面積を推定する工程
[6] 前記予測モデルの作成方法が、線形回帰、回帰木、モデル木、ニューラルネットワーク、サポートベクターマシン、バギング、ブースティング、ランダムフォレストの機械学習アルゴリズムから選ばれる1種を使用する方法である、[5]に記載の方法。 [7] 唾液中の1種以上の細菌の細菌量を検出し、得られた検出結果を指標として、歯周組織の炎症度を包括的に推定する方法。
[8] 歯周組織の炎症度が、PISA又はCAPRSの値である、[7]に記載の方法。
[9] 検出される細菌の細菌量が、唾液中の該細菌のコピー数である、[7] 又は[8]に記載の方法。
[10] 検出される細菌の細菌量が、唾液中の該細菌の16S rRNA配列情報に基づく細菌量である、[7]~[8]のいずれか1項に記載の方法。
[11] 検出される細菌が、Porphyromonas属、Tannerella属、Treponema属、Prevotella属、Campylobacter属、Fusobacterium属、Streptococcus属、Aggregatibacter属、Capnocytophaga属、Eikenella属、Actinomyces属、Veillonella属及びSelenomonas属から選ばれる少なくとも1種の属に属する細菌である、[7]~[10]のいずれか1項に記載の方法。
[12] 検出される細菌が、Streptococcus mutans、Actinomyces odontolyticus、Streptococcus mitis bv 2、Streptococcus mitis、Campylobacter concisus、Prevotella intermedia、Campylobacter showae、Prevotella nigrescens、Eikenella corrodens、Capnocytophaga gingivalis、Actinomyces naeslundii II、Streptococcus constellatus、Campylobacter gracilis、Fusobacterium periodonticum、Fusobacterium nucleatum subsp. polymorphum、Aggregatibacter actinomycetemcomitans、Capnocytophaga sputigena、Capnocytophaga ochracea、Streptococcus intermedius、Fusobacterium nucleatum subsp. vincentii、Streptococcus gordonii、Veillonella parvula、Selenomonas noxia、Fusobacterium nucleatum subsp. nucleatum、Campylobacter rectus、Porphyromonas gingivalis、Fusobacterium nucleatum subsp. animalis、Tannerella forsythia及びTreponema denticolaから選ばれる少なくとも1種である、[7]~[11]のいずれか1項に記載の方法。
[13] 検出される細菌は、該細菌の細菌量と歯周組織の炎症度とが正の相関関係となり得る細菌を含む、[7]~[12]のいずれか1項に記載の方法。
[14] 前記正の相関関係となり得る細菌が、Fusobacterium periodonticum、Fusobacterium nucleatum subsp. polymorphum、Aggregatibacter actinomycetemcomitans、Capnocytophaga sputigena、Capnocytophaga ochracea、Streptococcus intermedius、Fusobacterium nucleatum subsp. vincentii、Streptococcus gordonii、Veillonella parvula、Selenomonas noxia、Fusobacterium nucleatum subsp. nucleatum、Campylobacter rectus、Porphyromonas gingivalis、Fusobacterium nucleatum subsp. animalis、Tannerella forsythia及びTreponema denticolaから選ばれる少なくとも1種である、[13]に記載の方法。
[15] 検出される細菌は、該細菌の細菌量と歯周組織の炎症度とが負の相関関係となり得る細菌を含む、[7]~[12]のいずれか1項に記載の方法。
[16] 前記負の相関関係となり得る細菌が、Streptococcus mutans、Actinomyces odontolyticus、Streptococcus mitis bv 2、Streptococcus mitis、Campylobacter concisus、Prevotella intermedia、Campylobacter showae、Prevotella nigrescens、Eikenella corrodens、Capnocytophaga gingivalis、Actinomyces naeslundii II、Streptococcus constellatus及びCampylobacter gracilisから選ばれる少なくとも1種である、[15]に記載の方法。
本発明によれば、唾液中の細菌量の検出結果に基づいて、歯周ポケット炎症面積を簡便に予測することができる。すなわち、精密な歯周病検査(ポケット実測や撮像)をしなくとも、採取した唾液を用いて、簡便に口腔全体の炎症度を推定することができる。
また、本発明によれば、唾液中の細菌量の検出結果に基づいて、炎症面積(PISA値、CAPRS値)といった歯周組織の炎症度を簡便に予測することができる。すなわち、精密な歯周病検査(ポケット実測や撮像)をしなくとも、採取した唾液を用いて、簡便に口腔全体の炎症度(炎症指標数値)を推定することができる。
各細菌量のSN比から算出したPISA予測のモデル木を示す図である。 各細菌量のSN比から算出したPISA予測のモデル木を示す図である(図1-1の続き。図1-1の四角で囲った番号の1番と図1-2の21番がつながっている。)。 図1のモデル木から予測されたPISA値(横軸)とPISA実測値(縦軸)との散布図を示す図である。 チップ間補正後の各細菌量のSN比からモデル木を作成し、そのモデル木から予測されたPISA値(横軸)とPISA実測値(縦軸)との散布図を示す図である。 全46データの内34データを用い、チップ間補正後の各細菌量のSN比から算出したPISA予測のモデル木を示す図である。 全46データの内34データを用い、チップ間補正後の各細菌量のSN比から算出したPISA予測のモデル木を示す図である(図4-1の続き。図4-1の四角で囲った番号の1番と図4-2の11番、及び図4-1の13番と図4-2の12番がつながっている。)。 図4のモデル木作成に使用した34データでの予測されたPISA値(横軸)とPISA実測値(縦軸)との散布図を示す図である。 全46データの内、モデル作成につかっていない残り12データのPISA実測値(縦軸)と図4のモデル木から予測されたPISA値(横軸)との散布図を示す図である。 PISA実測値と「総菌量に対する3菌の比率(公知)」との散布図(左)、PISA実測値とPISA予測値との散布図(右)の比較を示す図である。 PISA実測値と、正の相関を示す細菌と負の相関を示す細菌の比(「バランス指標」)との散布図を示す図である。 56サンプルで試験を行い、PISA値(横軸)と統計的に有意な相関が認められた細菌種の細菌量を説明変数とし、重回帰分析で予測モデル式を作成してPISA値を予測(縦軸)した結果の散布図である。
以下、本発明を詳細に説明する。本発明の範囲はこれらの説明に拘束されることはなく、以下の例示以外についても、本発明の趣旨を損なわない範囲で適宜変更し実施することができる。なお、本明細書において引用された全ての刊行物、例えば先行技術文献、及び公開公報、特許公報その他の特許文献は、参照として本明細書に組み込まれる。
本発明の第一の態様である、歯周ポケット炎症面積を推定する方法の発明(以下、「第一の発明群」ともいう。)は、以下の工程を包含するものである。

i)唾液中の2種以上の細菌の細菌量を検出する工程、及び、ii)得られた検出結果を指標として歯周ポケット炎症面積を推定する工程。
本発明の第二の態様である、歯周組織の炎症度を包括的に推定する方法(以下、「第二の発明群」ともいう。)の発明は、以下の工程を包含するものである。
i)唾液中(被験体、特にヒトの唾液中)の1種以上の細菌の細菌量を検出する工程、及び、ii)得られた検出結果を指標として歯周組織の炎症度を包括的に推定する工程。
1.唾液中の細菌の細菌量の検出に用いるオリゴヌクレオチドプローブ
本発明の方法においては、被験者から採取された唾液中から口腔内の細菌量を測定する際に、DNAチップを使用することができる。当該DNAチップには、例えば、以下のプローブ(a)(細菌特異的プローブ)を搭載することができ、さらに、プローブ(b)(総量指標プローブ)やプローブ(c)(絶対量指標プローブ)を搭載することもできる。
(a)細菌特異的プローブ:検出対象となる細菌の遺伝子(又はその一部)に特異的にハイブリダイズする核酸からなるプローブ
(b)総量指標プローブ:すべての細菌の遺伝子にハイブリダイズする核酸からなるプローブ
(c)絶対量指標プローブ:1種類又は複数種類の絶対量指標それぞれに特異的にハイブリダイズする核酸からなるプローブ
なお、一般的に、DNAチップとは、プローブが配置された基盤の総称である。また、DNAチップ及びDNAマイクロアレイ等の名称については、それぞれ区別はせず、同義語であるとする。
(1)検出対象となる唾液中の細菌
本発明の方法において、検出対象(細菌量の測定対象)となる唾液中の細菌としては、限定はされないが、例えば、以下に列挙する各属に属する細菌、すなわち、Porphyromonas属、Tannerella属、Treponema属、Prevotella属、Campylobacter属、Fusobacterium属、Streptococcus属、Aggregatibacter属、Capnocytophaga属、Eikenella属、Actinomyces属、Veillonella属、Selenomonas属、Eubacterium属、Parvimonas属、Filifactor属、Haemophilus属、Alloprevotella属、Solobacterium属、Rothia属、Peptostreptococcus属、Gemella属、Corynebacterium属、Neisseria属、Granulicatella属、Megasphaera属及びSR1門から選ばれる少なくとも1種の属に属する細菌であることが好ましい。
より具体的には、以下に列挙する各種細菌から選ばれる少なくも1種、又は2種以上を検出対象とすることがより好ましい。
Porphyromonas gingivalis
Tannerella forsythia
Treponema denticola
Campylobacter gracilis
Campylobacter rectus
Campylobacter showae
Fusobacterium nucleatum subsp. vincentii
Fusobacterium nucleatum subsp. polymorphum
Fusobacterium nucleatum subsp. animalis
Fusobacterium nucleatum subsp. nucleatum
Fusobacterium periodonticum
Prevotella intermedia
Prevotella nigrescens
Streptococcus constellatus
Aggregatibacter actinomycetemcomitans
Campylobacter concisus
Capnocytophaga gingivalis
Capnocytophaga ochracea
Capnocytophaga sputigena
Eikenella corrodens
Streptococcus gordonii
Streptococcus intermedius
Streptococcus mitis
Streptococcus mitis bv 2
Actinomyces odontolyticus
Veillonella parvula
Actinomyces naeslundii II
Selenomonas noxia
Streptococcus mutans
Eubacterium nodatum
Parvimonas micra
Filifactor alocis
Streptococcus sobrinus
Porphyromonas pasteri
Veillonella atypica
Haemophilus parainfluenzae
Alloprevotella spp. (A. rava,OT 308)
Streptococcus parasanguinis
Actinomyces israelii
Prevotella pallens
Prevotella loescheii
Prevotella histicola
Solobacterium moorei
Prevotella melaninogenica
Selenomonas sputigena
Rothia dentocariosa
Rothia mucilaginosa
Veillonella rogosae
Peptostreptococcus stomatis
Prevotella denticola
Porphyromonas endodontalis
Streptococcus salivarius
Actinomyces graevenitzii
Treponema medium
Treponema socranskii
Gemella sanguinis
Porphyromonas catoniae
Corynebacterium matruchotii
Eubacterium saphenum
Neisseria flavescens
Granulicatella adiacens
Eubacterium sulci
Megasphaera micronuciformis
Prevotella shahii
SR1 sp. OT 345
「第一の発明群」においては、細菌の細菌量と歯周ポケット炎症面積とが正の相関関係(細菌量が増えたときに、歯周ポケット炎症面積は増える関係)を示す細菌(以下、「正の相関関係を示す細菌」と略記する場合がある。)、及び細菌の細菌量と歯周ポケット炎症面積とが負の相関関係(細菌量が増えたときに、歯周ポケット炎症面積は減る関係)を示す細菌(以下、「負の相関関係を示す細菌」と略記する場合がある。)が用いられる。
歯周ポケット炎症面積とは、PISA(periodontal inflamed surface area:歯周炎症表面積)やCAPRS(conccaled area in periodontal pockt of tooth root surfase : 歯周ポケット内歯根表面積)を含む炎症を示す面積であり、同様の概念の指標があればその指標も含まれる。
PISA値は、口腔全体の歯周組織の炎症部位面積を平方ミリメートル(mm)で示したものであり、歯周ポケット表面積(PESA)とプロービング時の出血(bleeding on probing:BOP)の有無から算出することができる。歯周ポケット表面積(PESA)は、歯の種類ごとに事前に規定されている面積と歯周ポケットの深さ(PPD)から算出することができる。6点法の場合の自動計算スプレッドシート(エクセルファイル)が非特許文献1の追加情報として配布されており(https://www.parsprototo.info/pisa.html)、上記計算方法はエクセルファイルに記載の計算式をみれば誰でも確認できる。
CAPRS値は、文献記載のとおり、アタッチメントレベルから臨床歯根表面積(歯肉縁より根尖側に位置する歯根の表面積の総和;Clinical Area of ToothRoot Surface; CARS)を求めた後、その値から有効歯根表面積の値を引いて算出される。実際には、歯肉縁の位置は全て解剖的歯頚線と一致しているケースに置き換えて算出され、これは非特許文献1のFig1.(b)に示されるPESAの値と同一であると考えられる。CAPRS値は、歯肉出血指数(BOP)が計算に加味されていないが、文献記載のとおり面積そのものをポケット内面の炎症領域と近似的に捉え、PISA値とともに「歯周ポケット炎症面積」と考えることができる。 歯周ポケット炎症面積は、PISA、又はCAPRSの値で表されるものであることが好ましい。
正の相関関係を示す細菌と負の相関関係を示す細菌は、細菌量(またはSN比のような細菌量と比例する測定量)を測定できるツールにより確認することができる。ツールは特に限定はされず、例えばDNAチップを用いることができる。
DNAチップを用いて確認する場合、口腔内試料をDNAチップで測定し、その後歯周ポケット炎症面積と各細菌の細菌量またはSN比のような測定量との相関係数を算出し、相関係数が正の値となる細菌群と負の値となる細菌群として分類、特定することができる。これらの細菌は、測定回数が40回以上の場合に相関係数の絶対値が0.02以上であることが好ましく、0.1以上であることがより好ましく、0.2以上であることがさらに好ましく、0.4以上であることが特に好ましく、0.6以上であることが最も好ましい。
後述する、歯周ポケット炎症面積を推定するための予測モデルの作成に実験誤差補正後のデータを用いる際は、細菌群の分類にも実験誤差補正後のデータを用いる。
正の相関関係を示す細菌としては、例えば、以下に列挙する細菌が好ましく挙げられる。これらのうちの少なくとも1種、好ましくは2種以上を検出対象とすることがより好ましい。
Treponema denticola
Tannerella forsythia
Fusobacterium nucleatum subsp. animalis
Porphyromonas gingivalis
Campylobacter rectus
Fusobacterium nucleatum subsp. nucleatum
Selenomonas noxia
Veillonella parvula
Streptococcus gordonii
Fusobacterium nucleatum subsp. vincentii
Streptococcus intermedius
Capnocytophaga ochracea
Capnocytophaga sputigena
Aggregatibacter actinomycetemcomitans
Fusobacterium nucleatum subsp. polymorphum
Fusobacterium periodonticum

SR1 sp. OT 345
Porphyromonas catoniae
Selenomonas sputigena
Neisseria flavescens
Streptococcus sobrinus
Parvimonas micra
Peptostreptococcus stomatis
Treponema socranskii
Eubacterium saphenum
Eubacterium nodatum
Treponema medium
Filifactor alocis
Porphyromonas endodontalis
負の相関関係を示す細菌としては、例えば、以下に列挙する細菌が好ましく挙げられる。これらのうちの少なくとも1種、好ましくは2種以上を検出対象とすることがより好ましい。

Streptococcus mutans
Actinomyces odontolyticus
Streptococcus mitis bv 2
Streptococcus mitis
Campylobacter concisus
Capnocytophaga gingivalis
Prevotella pallens
Streptococcus salivarius
Eubacterium sulci
Rothia mucilaginosa
Prevotella denticola
Veillonella atypica
Prevotella histicola
Megasphaera micronuciformis
Streptococcus parasanguinis
「第二の発明群」においては、細菌の細菌量が、歯周組織の炎症度等(PISA値又はCAPRS値)に対して、正の相関関係を有する細菌、又は、負の相関関係を有する細菌が、それぞれ好ましく挙げられる。
上記正の相関関係を有する細菌としては、例えば、以下に列挙する細菌が好ましく挙げられ、これらのうちの少なくとも1種、好ましくは2種以上を検出対象とすることがより好ましい。
Treponema denticola
Tannerella forsythia
Fusobacterium nucleatum subsp. animalis
Porphyromonas gingivalis
Campylobacter rectus
Fusobacterium nucleatum subsp. nucleatum
Selenomonas noxia
Veillonella parvula
Streptococcus gordonii
Fusobacterium nucleatum subsp. vincentii
Streptococcus intermedius
Capnocytophaga ochracea
Capnocytophaga sputigena
Aggregatibacter actinomycetemcomitans
Fusobacterium nucleatum subsp. polymorphum
Fusobacterium periodonticum
また、上記負の相関関係を有する細菌としては、例えば、以下に列挙する細菌が好ましく挙げられ、これらのうちの少なくとも1種、好ましくは2種以上を検出対象とすることがより好ましい。
Streptococcus mutans
Actinomyces odontolyticus
Streptococcus mitis bv 2
Streptococcus mitis
Campylobacter concisus
Capnocytophaga gingivalis
(2)細菌特異的プローブ
本発明において、細菌特異的プローブとして使用され得るオリゴDNAは、唾液中の細菌に由来する核酸の塩基配列のうちの特定の領域中の塩基配列とハイブリダイズすることができるものである。ここで、当該核酸は、染色体DNAやプラスミドDNA等を含むDNA及びRNAのいずれでもよく限定はされないが、染色体DNAであることが好ましい。具体的には、本発明においてプローブとして使用されるオリゴヌクレオチドは、前記細菌の染色体DNA中の16S rRNA遺伝子の塩基配列とハイブリダイズすることができるものである。
本発明に用い得るプローブは、検出目的となる個々の細菌に特異的な塩基配列となるような領域を選択してその領域の塩基配列を設計することが好ましい。一般的に、プローブの設計の際には、特異的な領域を選択することに加え、融解温度(Tm)がそろっていて、二次構造を形成しにくいものである必要がある。
唾液中の細菌の各々の種に対応する特異的な塩基配列は、例えば、マルチプルアラインメントをとり、種間で異なる領域にプローブを設計するなどの手段により見出すことができる。アラインメントをとるためのアルゴリズムには、特に限定はないが、より具体的な解析プログラムとしては、例えば、ClustalX1.8等のプログラムを利用することができる。アラインメントをとる際のパラメータは、各プログラムのデフォルト状態で実行してもよいが、プログラムの種類などに応じて適宜調整することができる。
プローブの特異性は、属レベルの特異性を基に同じ属の細菌を一括に検出するものであってもよいし、個々の種レベルで検出可能な特異性であってもよく、検出目的に応じて適宜選択・設計が可能である。
本発明において用い得る細菌特異的なプローブの例を、下記の表Aに示す(配列番号1~29)。
(3)総量指標プローブ
総量指標プローブは、特定のプライマー対で増幅できた、検体中(唾液中)のすべての細菌を捕捉する目的のプローブである。細菌を検出する上では、検出対象細菌が、非検出対象細菌を含む全体の細菌の中でどの程度の割合であるのか、また、そもそも検体中にどれくらいの量の細菌が存在しているのかといった観点から、細菌の総量を検出することも重要となる。
非検出対象細菌は、存在や種類は分かっているが検出対象としなくてもよい細菌、及び存在や種類が不明である細菌の和(合計)として理解できる。

細菌の総量を検出するためには、例えば、DNAチップとは独立に細菌の総量を測定することも可能であるが、DNAチップ中に細菌の総量の指標となるプローブを搭載しておくことにより操作の簡便性が向上する。プローブについては、プライマー対によって増幅される塩基配列の中から、多種類の菌種に共通な塩基配列を使用してもよい。そのような配列が見つからない場合は、比較的共通な配列を複数設計し、それらを総合的に判断することで総量指標プローブとしてもよい。総量指標プローブは、好ましくは、検体中に含まれる細菌に由来する核酸にハイブリダイズするプローブ、詳しくは、前記特定のプライマー対により増幅される塩基配列のうちの、検出対象となる複数種類の細菌が共通に有する塩基配列を含むプローブである。総量指標プローブの例を、下記の表A(配列番号31)に示す。
総量指標は、個々の菌種特異的な増幅産物の合計量を表すため、一般的に量が多くなることから、目的のシグナル強度が、検出可能なシグナル強度の範囲を超えてしまうことがある。

そのような状況を防ぐためには、ハイブリダイゼーションに供する検体量を制限することが望ましい。又は、プローブを設計する際には、例えば当該プローブのTm値を低くする。具体的にはGC含量を少なくすることや、プローブの配列長自体を短くする方法が考えられる。

また、ハイブリダイゼーションに際して、増幅された核酸と総量指標プローブとのハイブリダイゼーションに対して競合的に作用するような核酸を添加することで、シグナル強度の低減化を図ることが可能である。このような核酸としては、例えば、総量指標プローブと全て又は部分的に同じ配列を有する核酸、又は総量指標プローブの相補配列を全て又は部分的に有する核酸などが挙げられる。
(4)絶対量指標プローブ
絶対量指標プローブは、絶対量指標の核酸にのみハイブリダイズするプローブである。
本明細書において、絶対量指標とは、増幅反応やハイブリダイゼーション反応の前に、検体中に一定量添加する核酸である。絶対量指標は、通常の増幅反応を行えば増幅反応が確実に行われる核酸であり、いわゆる陽性コントロールとしての役割を果たす。
従って、絶対量指標に特異的なプローブを、DNAチップに搭載しておけば、その検出結果から、増幅反応やハイブリダイゼーション等が適切に実施されたかを確認することができる。また、絶対量指標を1種類設定した場合、多少増幅効率やハイブリダイゼーション効率が増減した場合に、絶対量指標のシグナル強度を比較することにより補正係数を算出することができる。複数のDNAチップにおいて、補正したシグナル強度は比較することができる。
絶対量指標用プローブの例を、下記の表A(配列番号30)に示す。
に示す。また絶対量指標の例を、下記の配列番号74に示す。
絶対量指標用プローブ:
CTATTCGACCAGCGATATCACTACGTAGGC(配列番号30)
絶対量指標:
GTGAGAAGCCTACACAAACGTAACGTCAGGGCTAAGACAAACGCTAACGGTACACCCTAGATGGGAGCTTGTAGCTAGATCGCTAAGTCCTACCGACATGTAGGCATACTCACGAAGGCAATTCCCTGAAAGCCTCGTCTTATCCCGAACTTGGCATCTGCTGATACGTCAGGTTGAACGCGTACATTTACCTGTCATGCGTGGGCCTTCTCCGAATAGCCTACGTAGTGATATCGCTGGTCGAATAGGCGGATTGCTCATAAATGCACATTGGCTAAGGCCCACGGAACACGAATCACGTGAGATCACTTACTATTCGACGGAACTACTATACGCACCGGGACATGCAAGTAGCGTCCCACAAGCATAAGGAACTCTATACTCGCCATCTACGCAGCTACAGGGGATACACGTATGAGCGGTTACGAAGTAAAGCCGAGATAGAGCGGTCTTTAGAGAAAAAACAGGATTAGATACCCTGGTAGTCC(配列番号74)
絶対量指標を増幅反応前に添加するのであれば、特定のプライマー対にて増幅される核酸であること、すなわちプライマー対と相補な塩基配列を所持していること、かつ、ハイブリダイゼーションで検出するためには、検出対象細菌、非検出対象細菌いずれにおいても所持していない塩基配列を所持している必要がある。
特定のプライマーとは、増幅対象配列が限定されるという意味であり、プライマー対は必ずしも1対である必要はない。必要に応じて2対以上のプライマー対を用いるマルチプレックス手法も適用できる。プライマー対の例を、下記の表Bに示す。細菌増幅用プライマー対(配列番号32、33)や、絶対量指標用プライマー対(配列番号34、35)を利用することが可能である。
絶対量指標は、例えば、産業技術総合研究所が開発した定量解析用核酸標準物質を利用しても良いし、新たに設計しても良い。設計する場合には、例えば、ソフトウェア「EXCEL」(MICROSOFT社製)のRNDBETWEEN関数を使用し、1から4までの整数をランダムにX個(Xは任意の数)発生させ、それをつなげて1から4までの数値のみから構成されるX桁の数値とし、1をA、2をT、3をC、4をGと置き換えることにより、ATGCのX塩基によるランダム配列を多数得ることができる。

これらの配列につき、GとTの和がAとTの和と同数になる配列のみを抜粋し、抜粋された配列を、NCBIのGenBank等のデータベースに対してBlast検索し、生物由来の核酸に対し、類似配列の少ないものを選抜し、配列の両末端にプライマー配列を付加することで、設計可能である。また、設計された配列を適宜連結させて長くしたり、部分的に除去して短くしたりすることも可能である。
増幅反応時における反応効率をなるべく一定にするために、検出対象細菌にて増幅される塩基長と絶対量指標の増幅塩基長は大きな差のないようにすることが望ましい。例えば、検出対象細菌の増幅産物が500bp程度となるのであれば、絶対量指標の増幅産物は300bpから1000bp程度とすることが望ましい。

一方で、増幅後に電気泳動等で増幅鎖長を確認する場合においては、検出対象細菌とは異なる長さの増幅産物となるように設計した上で、絶対量指標由来の増幅産物を検出対象細菌のバンドとは異なる位置で検出し、ハイブリダイゼーションの前に増幅反応の成否を確認することも可能である。

最後に、検体中に含まれる絶対量指標はあまりにも濃度が高いと検出対象の細菌と増幅反応における競合が激しくなり、本来検出できるはずの検出対象細菌が検出できなくなる可能性もあるため、アプリケーションに応じて適宜濃度調整をする必要がある。
Figure 0007224297000001
Figure 0007224297000002
本発明に用いるプローブを設計する際は、ハイブリダイゼーションにおけるストリンジェンシーを考慮することが好ましい。ストリンジェンシーをある程度緊密にすることによって、各種細菌において各々の核酸中の特定の領域間で類似する塩基配列領域が存在しても、他の異なる領域を区別してハイブリダイズすることができる。また、当該特定の領域間の塩基配列がほとんど異なる場合は、ストリンジェンシーを緩やかに設定することができる。
このようなストリンジェンシーの条件としては、例えば、緊密条件の場合は50~60℃の条件下でのハイブリダイゼーションであり、ゆるやかな条件の場合は30~40℃の条件下でのハイブリダイゼーションである。ハイブリダイゼーションの条件において、ストリンジェントな条件としては、例えば、「0.24M Tris・HCl/0.24M NaCl/0.05% Tween-20、40℃」、「0.24M Tris・HCl/0.24M NaCl/0.05% Tween-20、37℃」、「0.24M Tris・HCl/0.24M NaCl/0.05% Tween-20、30℃」、よりストリンジェントな条件としては、例えば「0.24M Tris・HCl/0.24M NaCl/0.05% Tween-20、50℃」、「0.24M Tris・HCl/0.24M NaCl/0.05% Tween-20、55℃」、「0.06M Tris・HCl/0.06M NaCl/0.05% Tween-20、60℃」等の条件を挙げることができる。より詳細には、プローブを添加して1時間以上50℃に保ってハイブリッド形成させ、その後、0.24M Tris・HCl/0.24M NaCl/0.05% Tween-20中、50℃で20分の洗浄を4回、最後に、0.24M Tris・HCl/0.24M NaCl、50℃で10分の洗浄を1回行う方法もある。ハイブリダイゼーション、又は洗浄の際の温度を上げることにより、よりストリンジェントな条件を設定することができる。当業者であれば、このようなバッファーの塩濃度、温度等の条件に加えて、その他のプローブ濃度、プローブの長さ、反応時間等の諸条件を加味し、条件を設定することができる。ハイブリダイゼーション法の詳細な手順については、「Molecular Cloning, A Laboratory Manual 4th ed.」(Cold Spring Harbor Press (2012)、「Current Protocols in Molecular Biology」(John Wiley & Sons (1987-1997)) 等を参照することができる。
本発明に用いるプローブの長さは、限定はされないが、例えば、10塩基以上が好ましく、より好ましくは16~50塩基であり、さらに好ましくは18~35塩基である。プローブの長さが適切であれば(前記範囲内であれば)、非特異的なハイブリダイゼーション(ミスマッチ)を抑制し、特異的な検出に使用することができる。
本発明に用いるプローブの設計の際には、Tmを確認しておくことが好ましい。Tmとは、任意の核酸鎖の50%がその相補鎖とハイブリッド形成する温度を意味し、鋳型DNA又はRNAとプローブとが二本鎖を形成してハイブリダイズするためには、ハイブリダイゼーションの温度を最適化する必要がある。一方、この温度を下げすぎると非特異的な反応が起こりやすくなるため、温度は可能な限り高いことが望ましい。従って、設計しようとする核酸断片のTmはハイブリダイゼーションを行う上で重要な因子である。Tmの確認には、公知のプローブ設計用ソフトウェアを利用することができ、本発明で利用可能なソフトウェアとしては、例えばProbe Quest(登録商標;ダイナコム社)などが挙げられる。またTmの確認は、ソフトウェアを使わずに自ら計算することによっても行うことができる。その場合には、最近接塩基対法(Nearest Neighbor Method)、Wallance法、GC%法等に基づく計算式を利用することができる。本発明のプローブにおいては、限定はされないが、平均Tmが約35~70℃又は45~60℃であることが好ましい。なお、プローブとして特異的なハイブリダイズが可能な条件としては、その他にもGC含量等があり、その条件は当業者に周知である。
また、本発明に用いるプローブを構成するヌクレオチドは、DNA及びRNA、又はPNAのいずれであってもよく、DNA、RNA及びPNAの2種以上のハイブリッドであってもよい。
本発明に用いるプローブとしては、具体的には、例えば、以下の(d)又は(e)のDNAの塩基配列を含むものが好ましく挙げられる。例えば、表Bに示すプライマー(配列番号32~35)を用いて増幅する場合、前掲の表A(配列番号1~31)に示す配列をプローブとすることが可能であり、配列番号1~31に示される塩基配列から選ばれる少なくとも2つの配列を用いることが好ましい。また、配列番号1~31に示される塩基配列から選ばれる少なくとも2つの配列の相補配列であってもよく、配列番号1~31に示される塩基配列から選ばれる少なくとも2つの配列と実質的に同一の配列又は配列番号1~31に示される塩基配列から選ばれる少なくとも2つの配列の相補配列と実質的に同一の配列であってもよい。
ここで、実質的に同一とは、配列番号1~31に記載の配列又は相補配列に対してストリンジェントな条件下で特異的にハイブリダイズするものである。
(d) 配列番号1~31に示される塩基配列からなるDNA
(e) 上記(d)のDNAに対し相補的な塩基配列からなるDNAとストリンジェントな条件下でハイブリダイズし、かつ唾液中の細菌に由来する核酸の塩基配列のうちの少なくとも一部の塩基配列を検出し得る機能を有するDNA
上記(d)の各種DNAについて、それらの具体的な塩基配列、プローブ名、検出対象となる口腔内細菌については、前掲の表Aの記載が参照できる。

また、上記(e)のDNAは、上記(d)の各種DNA若しくはそれと相補的な塩基配列からなるDNA、又はこれらを断片化したものをプローブとして用い、コロニーハイブリダイゼーション、プラークハイブリダイゼーション、及びサザンブロット等の公知のハイブリダイゼーション法を実施し、cDNAライブラリーやゲノムライブラリーから得ることができる。ライブラリーは、公知の方法で作製されたものを利用してもよいし、市販のcDNAライブラリーやゲノムライブラリーを利用してもよく、限定はされない。ハイブリダイゼーション法の詳細な手順については、前記と同様のものを参照することができる。上記(e)のDNAに関し、「ストリンジェントな条件」とは、ハイブリダイゼーション時の条件であって、バッファーの塩濃度が24~390mM、温度が40~65℃、好ましくは塩濃度が48.8~195mM、温度が45~60℃の条件を意味する。具体的には、例えば97.5mMで50℃等の条件を挙げることができる。さらに、このような塩濃度や温度等の条件に加えて、プローブ濃度、プローブの長さ、反応時間などの諸条件も考慮し、前記(e)のDNAを得るための条件を適宜設定することができる。ハイブリダイズするDNAとしては、前記(d)のDNAの塩基配列に対して少なくとも60%以上の相同性を有する塩基配列であることが好ましく、より好ましくは80%以上、さらに好ましくは90%以上、さらに好ましくは95%以上、特に好ましくは98%以上、最も好ましくは99%以上である。
本発明に用いるプローブは、例えば、通常のオリゴヌクレオチド合成法を使用して化学合成する(精製はHPLC等により行う)ことにより作製することができる。そのようなプローブは、例えば、Probe Quest(登録商標:ダイナコム社製)により設計することができる。また、本発明のプローブには、例えば、タグ配列などの付加配列が含まれていてもよい。
本発明の方法において、検出目的となる前述した唾液中の細菌が有する核酸の塩基配列は、当該塩基配列そのものである必要はなく、塩基配列の一部が欠失、置換、挿入等により変異が生じたものであってもよい。したがって、検出目的の核酸の塩基配列は、当該塩基配列に相補的な配列と、ストリンジェントな条件下でハイブリダイズし、かつそれぞれの塩基配列に由来する機能や活性を有する変異型遺伝子も対象とすることができ、プローブは、このような変異型遺伝子の塩基配列を基礎として設計することもできる。ここで「ストリンジェントな条件」は、前記と同様の条件を適用することができる。
2.DNAチップ
前記したように、本発明の方法においては、唾液中の細菌量の検出・測定に、DNAチップを用いることができる。当該DNAチップは、、歯周組織の炎症度を包括的に推定するという用途に用いられるものであり、前記1.項で説明した各種オリゴヌクレオチドプローブが支持体となる基盤に複数配置されたものである。
支持体となる基盤の形態としては、平板(ガラス板、樹脂板、シリコン板等)、棒状、ビーズ等のいずれの形態のものも使用できる。支持体として、平板を使用する場合は、その平板上に、所定の間隔もって、所定のプローブを種類毎に固定することができる(スポッティング法等;Science 270, 467-470 (1995)等参照)。また、平板上の特定の位置で、所定のプローブを種類毎に逐次合成していくこともできる(フォトリソグラフィー法等;Science 251, 767-773 (1991)等参照)。他の好ましい支持体の形態としては、中空繊維を使用するものが挙げられる。支持体として中空繊維を使用する場合は、所定のプローブを種類毎に各中空繊維に固定し、すべての中空繊維を集束させ固定した後、繊維の長手方向で切断を繰り返すことにより得られるDNAチップ(以下「繊維型DNAチップ」と言う)が好ましく例示できる。このマイクロアレイは、貫通孔基板に核酸を固定化したタイプのものと説明することもでき、いわゆる「貫通孔型DNAチップ」とも言われる(特許第3510882号公報等参照)。
プローブの支持体への固定方法は限定されず、どのような結合様式でもよい。また、支持体に直接固定することに限定はされず、例えば、予め支持体をポリリジン等のポリマーでコーティング処理し、処理後の支持体にプローブを固定することもできる。さらに、支持体として中空繊維等の管状体を使用する場合は、管状体にゲル状物を保持させ、そのゲル状物にプローブを固定することもできる。
以下、DNAチップの一形態である繊維型DNAチップに関して詳細に説明する。このDNAチップは、例えば、下記(i)~(iv)の工程を経て作製することができる。
(i) 複数本の中空繊維を、中空繊維の長手方向が同一方向となるように3次元に配列して配列体を製造する工程
(ii) 前記配列体を包埋し、ブロック体を製造する工程
(iii) オリゴヌクレオチドプローブを含むゲル前駆体重合性溶液を前記ブロック体の各中空繊維の中空部に導入して重合反応を行い、プローブを含むゲル状物を中空部に保持させる工程
(iv) 中空繊維の長手方向と交差する方向で切断して、ブロック体を薄片化する工程
中空繊維に使用される材料としては、限定はされないが、例えば、特開2004-163211号公報等に記載の材料が好ましく挙げられる。
中空繊維は、その長手方向の長さが同一となるように3次元に配列される(工程(i))。配列方法としては、例えば、粘着シート等のシート状物に複数本の中空繊維を所定の間隔をもって平行に配置し、シート状とした後、このシートを螺旋状に巻き取る方法(特開平11-108928号公報参照)や、複数の孔が所定の間隔をもって設けられた多孔板2枚を孔部が一致するように重ね合わせ、それらの孔部に中空繊維を通過させ、その後2枚の多孔板の間隔を開いて仮固定し、2枚の多孔板間における中空繊維の周辺に硬化性樹脂原料を充満させて硬化させる方法(特開2001-133453号公報参照)などが挙げられる。
製造された配列体はその配列が乱れないように包埋される(工程(ii))。包埋の方法としては、ポリウレタン樹脂及びエポキシ樹脂等を繊維間の隙間に流し込む方法のほか、繊維どうしを熱融着により接着する方法等が好ましく挙げられる。
包埋された配列体には、各中空繊維の中空部に、オリゴヌクレオチドプローブを含むゲル前駆体重合性溶液(ゲル形成溶液)を充填し、中空部内で重合反応を行う(工程(iii))。これにより、各中空繊維の中空部に、プローブが固定されたゲル状物を保持させることができる。
ゲル前駆体重合性溶液とは、ゲル形成重合性モノマー等の反応性物質を含有する溶液であって、該モノマー等を重合、架橋させることにより該溶液がゲル状物となることが可能な溶液をいう。そのようなモノマーとしては、例えば、アクリルアミド、ジメチルアクリルアミド、ビニルピロリドン、メチレンビスアクリルアミド等が挙げられる。この場合、溶液には重合開始剤等が含まれていてもよい。
中空繊維内にプローブを固定した後、中空繊維の長手方向と交差する方向(好ましくは直交する方向)で、ブロック体を切断して薄片化する(工程(iv))。このようにして得られた薄片は、DNAチップとして使用できる。当該DNAチップの厚みは、0.01mm~1mm程度であることが好ましい。ブロック体の切断は、例えば、ミクロトーム及びレーザー等により行うことができる
前記した繊維型DNAチップとしては、例えば、三菱ケミカル社製DNAチップ(Genopal TM)等が好ましく挙げられる。
繊維型DNAチップでは、前記のように、プローブはゲル内で3次元的に配列され、3次元構造を維持することが可能となる。そのため、表面をコートしたスライドガラスにプローブを結合させた平面DNAチップに比べて、検出効率が上昇し、高感度で高再現性の検査をすることが可能となる。
また、DNAチップに配置されるプローブの種類の数は、1つのDNAチップに500種類以下、好ましくは250種類以下、さらに好ましくは100種類以下が好ましい。このように配置されたプローブ数(種類)をある程度制限することにより、目的の口腔内細菌をより高感度で検出することが可能となる。なお、プローブの種類は塩基配列によって区別される。従って、通常、同じ遺伝子に由来のプローブであっても塩基配列が1個でも異なれば別の種類として特定する。
3.唾液中の細菌の検出(細菌量の測定)
本発明の方法において、唾液中の細菌量を測定するために当該細菌を検出する方法は、例えば、下記の工程を含む方法である。
(i)被験者から採取した口腔内試料としての唾液を検体とし、検体中(唾液中)の核酸を抽出する工程
(ii)抽出した核酸を、前記した本発明のオリゴヌクレオチドプローブ又は本発明のDNAチップに接触させる工程
(iii)DNAチップから得られたシグナル強度から細菌量を算出する工程
以下に、当該検出方法の詳細を工程ごとに説明する。
(1)工程(i)について
本工程では、被験者又は被生物から採取した口腔内試料としての唾液を検体とし、検体中(唾液中)に含まれる細菌の核酸を抽出する。唾液を採取する方法は特には限定されず、例えば、市販の唾液採取キットを利用する方法、綿棒を口に含み唾液を採取する方法、唾液を容器に直接採取する方法等が挙げられる。
唾液を採取する被験者は、特に限定はされないが、例えば、歯周病などの口内炎症を患っている患者に加え、歯周病などの口内炎症の自覚が無い者、又は心疾患など歯周病と関連が示唆される全身疾患患者や妊婦などでもよいし、あるいは、歯周病の罹患の疑いのない健常者であってもよい。
次いで、得られた唾液中に存在する細菌の核酸抽出を行う。抽出の方法は限定されず、公知の方法を用いることができる。例えば、機器による自動抽出法、市販の核酸抽出キットを利用する方法、プロテイナーゼK処理後にフェノール抽出する方法、クロロホルムを利用する方法又は簡易抽出方法として、試料を加熱、溶解する方法等が挙げられる。また、特に検体中から核酸を抽出せず、次の工程に進んでも良い。
検体から得られた核酸は、そのままDNAチップ等に接触させてもよいし、PCR等により所望の塩基配列領域を増幅し、その増幅断片をDNAチップ等に接触させてもよく、限定はされない。得られた核酸をテンプレートとして増幅する領域は、本発明に用いるプローブ、又はDNAチップに配置したオリゴヌクレオチドの塩基配列を含む核酸領域をコードする部位である。増幅する所望の領域は、限定はされず、細菌の種を問わず保存性の高い領域の塩基配列を利用し、多種類の混合物を一度に増幅して得ることができる。このような増幅のための配列は、実験的に単離、精製し、単離されたポリヌクレオチドの塩基配列を解析し、その配列に基づいて決定してもよいし、また、塩基配列等の各種データベースで既知の塩基配列を検索し、アラインメントを取ることなどによって、In Silicoで決定してもよい。核酸又はアミノ酸などのデータベースは、特に限定されるものではないが、例えば、DDBJ(DNA Data Bank of Japan)、EMBL(European Molecular Biology Laboratry, EMBL nucleic acid sequence data library)、GenBank(Genetic sequence data bank)、NCBI(National Center for Biotechnology Information)のTaxonomyデータベース等を利用できる。
具体的に、増幅する所望の部位としては、細菌の染色体DNA中のリボソームRNA (16S rRNA)遺伝子であることが好ましい。当該領域の増幅に用い得るPCRプライマーとしては、例えば、前掲の表Bの配列番号32、33が好ましく挙げられる。なお、PCR法による核酸の増幅は、定法に従って行うことができる。
本工程において抽出した核酸及びその増幅断片は、適宜標識化し、ハイブリダイズさせた後の検出過程において利用することも可能である。具体的には、PCRプライマーの末端を各種レポーター色素で標識しておく方法、反応性のヌクレオチドアナログを逆転写反応時に取り込ませる方法、ビオチン標識したヌクレオチドを取り込ませる方法などが考えられる。さらに、調製後に蛍光標識試薬と反応させて標識することも可能である。蛍光試薬としては、例えば、各種レポーター色素(例えば、Cy5、Cy3、VIC、FAM、HEX、TET、フルオレセイン、FITC、TAMRA、Texas red、Yakima Yellow等)を用いることができる。
(2)工程(ii)について
本工程では、工程(i)で得た核酸又はその増幅断片を、本発明に用いるプローブ又はDNAチップに接触させるが、具体的には、当該核酸等を含むハイブリダイゼーション溶液を調製し、当該溶液中の核酸等を、DNAチップに搭載されたオリゴヌクレオチドプローブに結合(ハイブリダイズ)させる。ハイブリダイゼーション溶液は、SDSやSSC等の緩衝液を用いて、定法に従い、適宜調製することができる。
ハイブリダイゼーション反応は、ハイブリダイゼーション溶液中の核酸等が、DNAチップに搭載されたオリゴヌクレオチドプローブとストリンジェントな条件下でハイブリダイズし得るよう、反応条件(緩衝液の種類、pH、温度等)を適宜設定して行うことができる。なお、ここで言う「ストリンジェントな条件」とは、類似配列によるクロスハイブリダイゼーションを生じにくい、又は類似配列によってクロスハイブリダイゼーションした核酸を解離させる条件のことをいい、具体的には、ハイブリダイゼーション反応時又はハイブリダイゼーション後のDNAチップの洗浄条件を意味する。
例えば、ハイブリダイゼーション反応時の条件としては、反応温度は、35~70℃が好ましく、より好ましくは40~65℃であり、ハイブリダイズさせる際の時間は、約1分~16時間が好ましい。
また、ハイブリダイゼーション後のDNAチップの洗浄条件としては、洗浄液組成は、0.24M Tris・HCl/0.24M NaCl/0.05% Tween-20であることが好ましく、洗浄時の温度は、35~80℃又は40~65℃が好ましく、より好ましくは45~60℃である。より具体的には、塩(ナトリウム)濃度が48~780mMであり、温度が37~80℃である条件が好ましく、より好ましくは塩濃度が97.5~390mMであり、温度が45~60℃である条件である。
洗浄後は、プローブに結合した核酸等の標識を検出できる装置により、スポットごとに検出強度を測定する。例えば、前記核酸等を蛍光標識化していた場合は、各種蛍光検出装置、例えば、CRBIO(日立ソフトウェアエンジニアリング社製)、arrayWoRx(GE Healthcare社製)、Affymetrix 428 Array Scanner(Affymetrix,社製)、GenePix(Axon Instruments社製)、ScanArray(PerkinElmer社製)、ジェノパールリーダー(三菱ケミカル社製)などを用いて、蛍光強度を測定することができる。これらの装置については、蛍光スキャナーの場合は、例えば、レーザーの出力、検出部の感度を適宜調整してスキャンを行うことができ、CCDカメラ型のスキャナーの場合は、露光時間を適宜調節してスキャンを行うことができる。スキャンの結果に基づく定量方法は、定量ソフトウェアにより行う。定量ソフトウェアに特に限定はなく、スポットの蛍光強度の平均値、中央値等を用いて定量することができる。また、定量にあたっては、DNAフラグメントのスポット範囲の寸法精度などを考慮し、プローブを搭載していないスポットの蛍光強度をバックグラウンドとして用いるなど、調整を行うことが好ましい。
(3)工程(iii)について
本工程では、前記の手順で得られたシグナル強度より、検出対象菌種の細菌の細菌量を算出する。たとえば、検出対象細菌を検出するためのプローブのシグナル強度とバックグラウンドのシグナル強度からSN比として示す方法がある。または、あらかじめ細菌ごとに細菌の染色体DNAの濃度を変えて複数条件にて検出し、各濃度条件で得られるシグナル強度をもとに、細菌ごとに染色体DNA濃度を算出する換算係数(検量線)を取得しておき、それぞれの条件で得られたシグナル強度から染色体DNAの濃度を算出する方法などが好ましい。本発明においては、例えば、検出対象の細菌の16S rRNA配列情報に基づく、シグナル強度から、細菌量を算出することが好ましい。また、検出対象の細菌のゲノムコピー数を、細菌量とすることも好ましく採用できる。ゲノムコピー数は、DNAチップで検出されたシグナル強度に、事前に決定していた各細菌量算出係数を乗算し(且つ、必要に応じ、検出検体の希釈率を乗算し)、算出できる。各細菌量算出係数は、各細菌由来ゲノムDNAを検出したときのシグナル強度を測定し検量線を作成しておき、各細菌のシグナル強度から各細菌量を逆算する係数として求めておくことができる。
いずれの場合でも、各DNAチップの検出対象細菌のシグナル強度に補正係数を考慮することが好ましい。
4.歯周組織の炎症度の推定
本発明の方法は、唾液中の細菌の細菌量の検出結果を指標として、歯周ポケット炎症面積を推定する方法である。また、本発明の方法は、前記検出結果を指標として、歯周組織の炎症度を包括的に推定する方法である。
唾液中の細菌の細菌量の検出には、いかなるツールを用いてもよく、前記3.項で説明したようにDNAチップを用いる方法や、その他、酵素活性で菌の存在を確認する方法、電気的抵抗を測定し菌の総量を測定する方法、位相差顕微鏡と染色により細菌を数え上げる方法、培養して生菌数を測定する方法、リアルタイムPCRで個々の細菌数を定量する方法等が挙げられる。
当該推定においては、細菌量、または細菌量と比例する測定量として、DNAチップの蛍光強度やシグナル強度のSN比、リアルタイムPCRのCt値、酵素活性値、電気的に測定した場合の抵抗値、目視で数え上げた計数値等に基づいて、歯周ポケット炎症面積や歯周組織の炎症度を推定する。
歯周ポケット炎症面積を推定する具体的な方法としては、以下の方法が挙げられる。
(1)歯周ポケット炎症面積(PISA値、CAPRS値等)が既知の(実測されたPPD等から算出可能な)被験者の唾液サンプルから、唾液中の各種細菌の細菌量を検出する。
(2)当該各種細菌の細菌量について、菌ごとに固有の歯周ポケット炎症面積との相関係数を求め、各種細菌の細菌量と歯周ポケット炎症面積との関係式を構築して予測モデルを作成する。
(3)歯周ポケット炎症面積が未知の被験者の唾液サンプルから、唾液中の各種細菌の細菌量を検出する。
(4)(3)で得られた各種細菌量を(2)で得られた関係式に入れて歯周ポケット炎症面積を推定する。
予測モデルの作成方法は、特に限定はされないが、例えば、線形回帰、回帰木、モデル木、ニューラルネットワーク、サポートベクターマシン、バギング、ブースティング、ランダムフォレスト等の機械学習アルゴリズムなど、各種統計解析法を使用する方法があげられる。このうち、後述する実施例で示したモデル木では、あらかじめモデルを規定する必要が無く、具体的には、統計ソフトウェア「R」(R Development Core Team)の「caret」パッケージを用いた、「M5」メソッドでの最適化が、好ましく挙げられる。
歯周ポケット炎症面積が実測されている被験者の唾液サンプル(予測モデル作成に使用するデータ数)は、予測モデル作成に使用する細菌数(変数)以上の数を用いることが好ましい。データ数が蓄積されるごとに予測モデルを更新してもよい。
「第一の発明群」においては、検出される予測モデル作成に使用する細菌として、細菌の細菌量と歯周ポケット炎症面積と正の相関関係を示す細菌と、細菌の細菌量と歯周ポケット炎症面積と負の相関関係を示す細菌の、両方を用いる。歯周病の進行がない健康な人では、歯周ポケット炎症面積と正の相関関係を示す細菌が検出されないことも多く(細菌数0となってしまう)、正の相関関係を示す細菌のみでは、歯周ポケット炎症面積が小さい値の範囲について予測モデルが作成できないためである。
予測モデル作成に使用する細菌の例示は前記1.項に記載した通りであるが、最終的に予測精度が向上するように、細菌は、菌単独の相関係数の大きさだけでなく、歯周ポケット炎症面積の変動に応答するその菌の変化量の正確さや、菌同士の細菌量の相関係数を考慮(多重共線関係を避ける)して選択される。
正の相関を示す細菌として、好ましくは以下の細菌が挙げられる。
Porphyromonas gingivalis,Tannerella forsythia,Treponema denticola,Campylobacter rectus,Fusobacterium nucleatum subsp. vincentii,Fusobacterium nucleatum subsp. polymorphum,Fusobacterium nucleatum subsp. animalis,Fusobacterium nucleatum subsp. nucleatum,Fusobacterium periodonticum,Aggregatibacter actinomycetemcomitans,Capnocytophaga ochracea,Capnocytophaga sputigena,Streptococcus gordonii,Streptococcus intermedius,Veillonella parvula,Selenomonas noxia,Solobacterium moorei,Prevotella loescheii,Veillonella rogosae,Actinomyces israelii,Corynebacterium matruchotii,SR1 sp. OT 345,Porphyromonas catoniae,Selenomonas sputigena,Neisseria flavescens,Streptococcus sobrinus,Parvimonas micra,Peptostreptococcus stomatis,Treponema socranskii,Eubacterium saphenum,Eubacterium nodatum ,Treponema medium,Filifactor alocis,Porphyromonas endodontalis
正の相関を示す細菌として、さらに好ましくは以下の細菌が挙げられる。
,Porphyromonas gingivalis,Tannerella forsythia,Treponema denticola,Campylobacter rectus,Fusobacterium nucleatum subsp. animalis,Fusobacterium nucleatum subsp. nucleatum,Veillonella parvula,Selenomonas noxia,Eubacterium saphenum,Eubacterium nodatum ,Treponema medium,Filifactor alocis,Porphyromonas endodontalis
正の相関を示す細菌としては、1種以上を用いることが好ましく、4種以上を用いることがより好ましく、8種以上を用いることがさらに好ましく、12種以上を用いることが特に好ましい。また、100種以下を用いることが好ましく、75種以下を用いることがより好ましく、50種以下を用いることがさらに好ましく、25種以下を用いることが特に好ましい。

負の相関祖示す細菌として、好ましくは以下の細菌が挙げられる。Streptococcus mitis,Streptococcus mitis bv 2,Actinomyces odontolyticus,Streptococcus mutans,Campylobacter concisus,Capnocytophaga gingivalis,Prevotella pallens,Streptococcus salivarius,Eubacterium sulci,Rothia mucilaginosa,Prevotella denticola,Veillonella atypica,Prevotella histicola,Megasphaera micronuciformis,Streptococcus parasanguinis,Gemella sanguinis,Alloprevotella spp. (A. rava,OT 308),Prevotella melaninogenica,Actinomyces graevenitzii,Prevotella shahii,Rothia dentocariosa,Granulicatella adiacens,Porphyromonas pasteri,Haemophilus parainfluenzae
負の相関祖示す細菌として、さらに好ましくは以下の細菌が挙げられる。,Actinomyces odontolyticus,Streptococcus mutans,Prevotella pallens

負の相関を示す細菌としては、1種以上を用いることが好ましく、2種以上を用いることがより好ましく、4種以上を用いることがさらに好ましく、8種以上を用いることが特に好ましい。また、100種以下を用いることが好ましく、75種以下を用いることがより好ましく、50種以下を用いることがさらに好ましく、25種以下を用いることが特に好ましい。
ところで、本発明においては、予め規定された数の被験者(1次母標本)において統計解析処理を行い、データベースに記憶させておく。そして、歯周組織の炎症度と唾液中の細菌量との間の相関関係の解析結果から、被験者のそれぞれがどのような歯周組織の炎症度を有するのか、又は歯周ポケット内歯根表面積を有するのかを推定することができる。従って、被験者個人(一人)の検査を行う場合、上記複数の被験者由来のデータを母標本として、当該被験者個人のデータが、データベースに記憶させておいた母標本のデータのどこに位置するか又は当てはまるかを調べることによって、当該被験者個人に対する歯周組織の炎症度又は歯周ポケット内歯根表面積を推定することができる。なお、上記被験者個人のデータを母標本の値に組み込み、再度統計解析処理した後、当該被験者個人が母標本のどこに位置するかを調べるようにしてもよい。
本発明の方法によれば、唾液中の細菌種及びその細菌量に基づいて、歯周組織の炎症度(PISA値、CAPRS値)を推定・予測することができるため、従来、PISA値やCAPRS値を実測していた場合に比べて、格段に簡便に且つ大量の被験体に対しても、一定の算出基準の下で、歯周組織の炎症度を算出することができる。また、本発明においては、歯周組織の炎症度としては、陽性度だけではなく陰性度についても推定・予測の範囲に含まれる。よって、歯周組織の炎症度に対して、細菌量が逆相関する細菌種(常在菌)に着目した場合は、被験体の口腔内について健常状態が維持されているかどうかの推定・予測も行うことができる。さらには、歯周組織の炎症度に対して、細菌量が相関する細菌種、逆相関する細菌種の比を指標としてもよい。
以下に、実施例を挙げて本発明をより具体的に説明するが、本発明はこれらに限定されるものではない。
唾液サンプルの採取
大阪大学歯学部附属病院にて、歯周病治療中の20代~70代の男女被験者計46人から、唾液1ml以上、46サンプルを採取した。採取した唾液は使用前まで-20℃で凍結して保管した。
唾液サンプルと対応付けられたPISA値の算出
<臨床情報の取得>
すべてのサンプルについて下記の臨床情報を記録した。下記2項目は歯科において広く活用されている指標である。
(i)歯周ポケットの深さ(PPD):歯周プローブをポケットに挿入した際の,歯肉辺縁からプローブ先端までの距離を示す。6点法で頬側近心、頬側中央、頬側遠心、舌側近心、舌側中央、舌側遠心を測定し、1mm単位で数値化した。
(ii)プロービング時の出血(BOP):歯周プローブをポケットに挿入した際に出血の有無を示す。上記の6点法と対応した位置に対して出血がない場合を0、出血がある場合を1とした。
(i)歯周ポケットの深さ(PPD)と(ii)プロービング時の出血(BOP)を、表1にまとめた。
Figure 0007224297000003
Figure 0007224297000004
Figure 0007224297000005
Figure 0007224297000006
Figure 0007224297000007
Figure 0007224297000008
Figure 0007224297000009
Figure 0007224297000010
Figure 0007224297000011
Figure 0007224297000012
Figure 0007224297000013
Figure 0007224297000014
Figure 0007224297000015
Figure 0007224297000016
Figure 0007224297000017
Figure 0007224297000018
<PISA値の算出>
続いて、記録した6点法による歯周病検査結果から、公知文献(Nesse, W., Abbas, F., van der Ploeg, I., Spijkervet, F.K.,Dijkstra ,P.U., Vissink, A.: Periodontal inflamed surface area: quantifying inflammatory burden. Journal of Clinical Periodontology, 35: 668-673, 2008.)を参照し、表2に示すようにPISA(Periodontal inflamed surface area)を算出した。PISAは口腔全体の歯周組織の炎症部位面積を平方ミリメートル(mm2)で示したものである。
表3に、46名のPISA値をまとめて示した。

Figure 0007224297000019
Figure 0007224297000020
Figure 0007224297000021
唾液サンプル中の口腔内細菌の検出
<唾液からのDNAの抽出>
1-1 ビーズ破砕
2mlのチューブに前述の唾液のうち400μl、滅菌水400μlと0.1mmガラスビーズ0.4g(トミー精工社製#GB-01)を入れフタをしっかり締めた後、マイクロスマッシュ(トミー精工社製)にセットし、3200rpmで5分間ビーズ破砕した。その後、ビーズ等を落とすためにチューブを遠心機にセットし、1800gで5分間軽く遠心した。上清400μlを新しい1.5mlチューブに分注した。
1-2 酵素処理
1-1のサンプルに対してキアゲン社のBuffer ATL90μl、ProteinaseK 10μl、Buffer AL100μlを混合し、合計600μlとし、56℃で1時間インキュベートした。
1-3 カラム精製
インキュベート後、エタノール120μlを加えてボルテックスした。続いてQIAamp DNA Micro Kit(キアゲン社)を使用し、キットのカラムに720μl全量を入れて14000rpmで1分間遠心した後、溶出液を廃棄した。その後、キットのプロトコールに従い、Buffer AW1 500μl、Buffer AW2 500μlでカラムを洗浄した。続いて、空の状態で14000rpmで1分間遠心した後、滅菌水20μlでDNAを溶出した。46サンプルでそれぞれ約4~20ng/μlのDNA溶液を取得した。
<細菌DNAの増幅反応>
前記のDNA溶液を希釈して20pg/μlとし、これをPCRテンプレートとした。サンプル中の口腔内細菌の16SrRNAの検出対象領域の配列を増幅するために、以下の反応液組成及び反応条件でPCRを実施した。PCR用キットは、Premix Ex Taq(TM) Hot Start Version(Takara社製)を用い、ProFlex(TM)PCR System(Thermo Fisher Scientific社製)により行った。プライマーは下記の配列を有するプライマーを用いた。なお、フォワードプライマーは5’末端がCy5で標識化されているものを用いた。
フォワードプライマー(細菌増幅用):
5’-Cy5-TCCTACGGGAGGCAGCAGT-3’ (配列番号32)
リバースプライマー(細菌増幅用):
5’-CAGGGTATCTAATCCTGTTTGCTACC-3’ (配列番号33)
フォワードプライマー(絶対量指標増幅用):
5’-Cy5-GAGAAGCCTACACAAACGTAACGTC-3’ (配列番号34)
リバースプライマー(絶対量指標増幅用):
5’-CTCTAAAGACCGCTCTATCTCGG-3’(配列番号35)
<反応液組成>
2×Premix Ex Taq(登録商標) Hot Start Version 10μL
4μMフォワードプライマー(細菌増幅用) 1μL
4μMリバースプライマー(細菌増幅用) 1μL
4μMフォワードプライマー(絶対量指標増幅用) 1μL
4μMリバースプライマー(絶対量指標増幅用) 1μL
テンプレートDNA 5μL
絶対量指標 1μL
合計 20μL
<反応条件>
95℃で1分間加熱後、「解離:98℃(10sec)→アニーリング:55℃(30sec)→合成:72℃(20sec)」を1サイクルとして計40サイクル行い、4℃で冷却し、増幅産物を得た。
<DNAチップ:口腔内細菌検出用DNAチップの製造>
貫通孔型のDNAチップを、特開2007-74950号公報(メチル化DNA及び/又は非メチル化DNAの検出方法)の実施例1に記載の方法と同様の方法で製造を行った。
ただし、搭載させたオリゴヌクレオチドプローブは、表4に示す配列情報をもつプローブを用いた。
Figure 0007224297000022
<DNAチップへのハイブリダイゼーション>
以下のように各溶液を混合し、ハイブリダイゼーション溶液を調製した。
PCR後に得たDNA増幅産物 20μL
1M Tris-HCl 48μL
1M NaCl 48μL
0.5% Tween20 20μL
水 65μL
合計 200μL
DNAチップのハイブリダイゼーション、洗浄には自動ハイブリ洗浄装置(型式:AHF-200、三菱ケミカル社製)を用いた。
200μLのハイブリダイゼーション溶液を前記DNAチップに接触させ、50℃で16時間ハイブリダイゼーションした。
ハイブリダイゼーション後、下記の条件でDNAチップを洗浄した。0.24M Tris・HCl/0.24M NaCl/0.05% Tween-20溶液1000μLで220秒の洗浄を12回繰り返し、続いて、0.24M Tris・HCl/0.24M NaCl1000μLで220秒の洗浄を4回繰り返した。
洗浄終了後に、各チップを室温の0.24M Tris・HCl/0.24M NaCl混合溶液に移した。
<検出>
前記洗浄後、ジェノパールリーダー(型式:GR-S1、三菱ケミカル社製)を用い、下記条件でDNAチップの各スポットの蛍光強度を測定した。
<検出条件>
中心励起波長 :633nm
露光時間 :0.1、1、4、40秒
<結果>
検出対象細菌用プローブを搭載したスポットの蛍光強度を、バックグラウンド値(プローブを搭載していないスポットの蛍光強度の中央値)で除算し、ハイブリダイゼーションに由来する蛍光強度(以後、シグナル強度と呼ぶ)のSN比を算出した。
細菌量に基づいたPISA値の予測
<PISA値と細菌量の相関解析>
まず、サンプルごとに、PISA値および各細菌量を示すSN比のデータを対応させた。その結果を、表5に示した。
Figure 0007224297000023
Figure 0007224297000024
Figure 0007224297000025
Figure 0007224297000026
その後、エクセルの「データ分析」を用いてPISAと細菌のSN比との相関係数を求め、表6にまとめた。
Figure 0007224297000027
各細菌量を示すSN比と口腔内全体のPISA値との比較において、相関が高いものから順にTreponema denticola、Tannerella forsythia、Fusobacterium nucleatum subsp.animalis、Porphyromonas gingivalis、Campylobacter rectusであった。これらには、「Red Complex」が含まれており、口腔全体の炎症度の指標になると考えられた。一方、逆相関を示すものとしては、Streptococcus mutans、Actinomyces odontolyticus、Streptococcus mitis bv 2、Streptococcus mitis、Campylobacter concisusであった。これらは口腔全体の健康度の指標になると考えられた。
続いて、表5に示したデータを用いて、各細菌量のSN比をもとに機械学習の手法を用い、PISA値をモデル木で予測するモデルの作成を行った。解析には、統計ソフトウェア「R」(R Development Core Team)の「caret」パッケージを用い、「M5」メソッドでの最適化を実行した。
表5のデータを「bacteria」として読み込ませた後、次のコマンドを実行した。
m <- train(PISA~.,data=bacteria ,method="M5")
PISAを目的変数とし、説明変数を表5の細菌全種類としてモデル木を生成させるコマンドである。46サンプルデータすべてを使用した。
実行した最適モデルをm$finalModelと結果の出力コマンドを入力して出力させたところ、以下の予測モデルを得た。結果を、図1に示した。
M5 unpruned model tree:
(using smoothed linear models)

Actinomyces.odontolyticus <= 20.1 :
| Campylobacter.concisus <= 2.2 :
| | Capnocytophaga.sputigena <= 2.15 :
| | | Tannerella.forsythia <= 6.9 :
| | | | Streptococcus.mutans <= 5.3 : LM1 (2/12.282%)
| | | | Streptococcus.mutans > 5.3 : LM2 (2/32.179%)
| | | Tannerella.forsythia > 6.9 :
| | | | Streptococcus.mutans <= 2.95 : LM3 (3/16.447%)
| | | | Streptococcus.mutans > 2.95 : LM4 (3/2.075%)
| | Capnocytophaga.sputigena > 2.15 :
| | | Streptococcus.gordonii <= 35.45 : LM5 (2/1.965%)
| | | Streptococcus.gordonii > 35.45 : LM6 (3/26.323%)
| Campylobacter.concisus > 2.2 :
| | Treponema.denticola <= 15.35 :
| | | Fusobacterium.nucleatum.subsp..nucleatum <= 38.35 : LM7 (3/2.882%)
| | | Fusobacterium.nucleatum.subsp..nucleatum > 38.35 :
| | | | Actinomyces.odontolyticus <= 10.9 : LM8 (3/4.64%)
| | | | Actinomyces.odontolyticus > 10.9 : LM9 (2/2.948%)
| | Treponema.denticola > 15.35 : LM10 (3/39.648%)
Actinomyces.odontolyticus > 20.1 :
| Tannerella.forsythia <= 2.85 :
| | Streptococcus.gordonii <= 19.6 :
| | | Actinomyces.naeslundii.II <= 4.05 : LM11 (5/2.414%)
| | | Actinomyces.naeslundii.II > 4.05 : LM12 (3/3.304%)
| | Streptococcus.gordonii > 19.6 :
| | | Control.DNA <= 837.15 : LM13 (3/4.132%)
| | | Control.DNA > 837.15 : LM14 (2/0.098%)
| Tannerella.forsythia > 2.85 :
| | Streptococcus.constellatus <= 1.75 : LM15 (3/14.832%)
| | Streptococcus.constellatus > 1.75 :
| | | Porphyromonas.gingivalis <= 1.35 : LM16 (2/9.58%)
| | | Porphyromonas.gingivalis > 1.35 : LM17 (2/2.751%)

LM num: 1
.outcome =
5.9609 * Treponema.denticola
+ 2.719 * Streptococcus.gordonii
- 2.4637 * Actinomyces.odontolyticus
- 43.3698 * Streptococcus.mutans
+ 161.5714 * Capnocytophaga.sputigena
- 0.4833 * Control.DNA
+ 3.734 * Tannerella.forsythia
+ 1510.5973

LM num: 2
.outcome =
5.9609 * Treponema.denticola
+ 2.719 * Streptococcus.gordonii
- 2.4637 * Actinomyces.odontolyticus
- 43.3698 * Streptococcus.mutans
+ 161.5714 * Capnocytophaga.sputigena
- 0.4833 * Control.DNA
+ 3.734 * Tannerella.forsythia
+ 1507.5864

LM num: 3
.outcome =
5.9609 * Treponema.denticola
+ 2.719 * Streptococcus.gordonii
- 2.4637 * Actinomyces.odontolyticus
- 26.3717 * Streptococcus.mutans
+ 161.5714 * Capnocytophaga.sputigena
- 0.4833 * Control.DNA
+ 3.9077 * Tannerella.forsythia
+ 1392.5945

LM num: 4
.outcome =
5.9609 * Treponema.denticola
+ 2.719 * Streptococcus.gordonii
- 2.4637 * Actinomyces.odontolyticus
- 26.3717 * Streptococcus.mutans
+ 161.5714 * Capnocytophaga.sputigena
- 0.4833 * Control.DNA
+ 3.9077 * Tannerella.forsythia
+ 1394.1708

LM num: 5
.outcome =
5.9609 * Treponema.denticola
+ 1.2546 * Streptococcus.gordonii
- 2.4637 * Actinomyces.odontolyticus
- 26.3717 * Streptococcus.mutans
+ 201.9643 * Capnocytophaga.sputigena
- 0.4833 * Control.DNA
+ 5.5586 * Tannerella.forsythia
+ 1565.177

LM num: 6
.outcome =
5.9609 * Treponema.denticola
+ 1.336 * Streptococcus.gordonii
- 2.4637 * Actinomyces.odontolyticus
- 26.3717 * Streptococcus.mutans
+ 201.9643 * Capnocytophaga.sputigena
- 0.4833 * Control.DNA
+ 5.5586 * Tannerella.forsythia
+ 1555.611

LM num: 7
.outcome =
13.3672 * Treponema.denticola
+ 2.719 * Streptococcus.gordonii
- 2.4637 * Actinomyces.odontolyticus
- 30.4289 * Streptococcus.mutans
- 0.4833 * Control.DNA
+ 5.5586 * Tannerella.forsythia
- 0.3524 * Fusobacterium.nucleatum.subsp..nucleatum
+ 1314.4003

LM num: 8
.outcome =
13.3672 * Treponema.denticola
+ 2.719 * Streptococcus.gordonii
- 1.3817 * Actinomyces.odontolyticus
- 30.4289 * Streptococcus.mutans
- 0.4833 * Control.DNA
+ 5.5586 * Tannerella.forsythia
- 0.3171 * Fusobacterium.nucleatum.subsp..nucleatum
+ 1293.9573

LM num: 9
.outcome =
13.3672 * Treponema.denticola
+ 2.719 * Streptococcus.gordonii
- 1.3181 * Actinomyces.odontolyticus
- 30.4289 * Streptococcus.mutans
- 0.4833 * Control.DNA
+ 5.5586 * Tannerella.forsythia
- 0.3171 * Fusobacterium.nucleatum.subsp..nucleatum
+ 1294.0756

LM num: 10
.outcome =
15.293 * Treponema.denticola
+ 2.719 * Streptococcus.gordonii
- 2.4637 * Actinomyces.odontolyticus
- 30.4289 * Streptococcus.mutans
- 0.4833 * Control.DNA
+ 5.5586 * Tannerella.forsythia
+ 1320.5275

LM num: 11
.outcome =
3.3786 * Treponema.denticola
+ 5.1956 * Streptococcus.gordonii
- 2.8861 * Actinomyces.odontolyticus
- 0.5662 * Control.DNA
+ 29.9183 * Tannerella.forsythia
- 2.6289 * Actinomyces.naeslundii.II
+ 826.1665

LM num: 12
.outcome =
3.3786 * Treponema.denticola
+ 5.1956 * Streptococcus.gordonii
- 2.8861 * Actinomyces.odontolyticus
- 0.5662 * Control.DNA
+ 29.9183 * Tannerella.forsythia
- 2.9643 * Actinomyces.naeslundii.II
+ 825.1613

LM num: 13
.outcome =
3.3786 * Treponema.denticola
+ 5.4972 * Streptococcus.gordonii
- 2.8861 * Actinomyces.odontolyticus
- 0.6637 * Control.DNA
+ 29.9183 * Tannerella.forsythia
+ 0.4481 * Actinomyces.naeslundii.II
+ 906.7614

LM num: 14
.outcome =
3.3786 * Treponema.denticola
+ 5.4972 * Streptococcus.gordonii
- 2.8861 * Actinomyces.odontolyticus
- 0.6694 * Control.DNA
+ 29.9183 * Tannerella.forsythia
+ 0.4481 * Actinomyces.naeslundii.II
+ 910.6162
LM num: 15
.outcome =
-215.5803 * Porphyromonas.gingivalis
+ 3.3786 * Treponema.denticola
+ 3.1851 * Streptococcus.gordonii
- 2.8861 * Actinomyces.odontolyticus
- 56.1926 * Streptococcus.constellatus
- 0.5662 * Control.DNA
+ 36.302 * Tannerella.forsythia
+ 1327.7476

LM num: 16
.outcome =
-221.29 * Porphyromonas.gingivalis
+ 3.3786 * Treponema.denticola
+ 3.1851 * Streptococcus.gordonii
- 2.8861 * Actinomyces.odontolyticus
- 53.2351 * Streptococcus.constellatus
- 0.5662 * Control.DNA
+ 36.302 * Tannerella.forsythia
+ 1307.8776

LM num: 17
.outcome =
-221.29 * Porphyromonas.gingivalis
+ 3.3786 * Treponema.denticola
+ 3.1851 * Streptococcus.gordonii
- 2.8861 * Actinomyces.odontolyticus
- 53.2351 * Streptococcus.constellatus
- 0.5662 * Control.DNA
+ 36.302 * Tannerella.forsythia
+ 1307.1008
p <- predict(m,newdata=bacteria)とコマンドを入力し、作成した予測モデル「m」に対して、実測した各細菌量のSN比データの46サンプル分を代入し、46サンプルに対応する予測PISA値「p」得た。
引き続き、cor(p,bacteria$PISA)をコマンドを入力し、予測PISA値「p」と実測のPISA値の相関係数を計算させたところ、相関係数0.8759666となった。
予測PISA値「p」と実測のPISA値の散布図を図2に示す。
このモデルに示すようにDNAチップのSN比からPISA値を予測可能であることが示された。
特に表6に示した各細菌量を示すSN比と口腔内全体のPISA値との比較において、負の相関係数をもつ、Streptococcus mutans、Actinomyces odontolyticus、Campylobacter concisus、Actinomyces.naeslundii.II、Streptococcus.constellatusと、正の相関係数をもつCapnocytophaga sputigena、Tannerella forsythia、Streptococcus gordonii、Treponema denticola、Fusobacterium nucleatum subsp..nucleatum、Porphyromonas gingivalisの細菌群のSN比からPISA値を予測可能であることが示された。
実施例1と同一の検出結果を用いて細菌量に基づいたPISA値の予測を行った。
<PISA値と細菌量の相関解析>
まず、サンプルごとに、PISA値および各細菌量を示すSN比のデータを対応させた。その後、すべてのサンプルの絶対量指標プローブのSN比の中央値を算出し、個々のサンプルの絶対量指標プローブのSN比をこの中央値で除算した。これをチップ間補正値とする。続いて、個々のサンプルのすべての細菌のSN比データを、個々のサンプルのチップ間補正値で除算することにより、チップ間のSN比データを標準化した。その結果を、表7に示した。
Figure 0007224297000028
Figure 0007224297000029
Figure 0007224297000030
Figure 0007224297000031
その後、エクセルの「データ分析」を用いてPISAと細菌のSN比との相関係数を求め、表8に示した。
Figure 0007224297000032
各細菌量を示すSN比と口腔内全体のPISA値との比較において、相関が高いものから順にTannerella forsythia、Treponema denticola、Porphyromonas gingivalis、Fusobacterium nucleatum subsp.animalis、Selenomonas noxia、Streptococcus gordonii、Veillonella parvula、Campylobacter rectus、Fusobacterium nucleatum subsp. nucleatumであった。これらには、「Red Complex」が含まれており、口腔全体の炎症度の指標になると考えられた。一方、逆相関を示すものとしては、Streptococcus mutans、Actinomyces odontolyticus、Streptococcus mitis bv 2、Streptococcus mitis、Campylobacter concisusであった。これらは口腔全体の健康度の指標になると考えられた。
続いて、表7に示したデータを用いて、各細菌量のSN比をもとに機械学習の手法を用い、PISA値をモデル木で予測するモデルの作成を行った。解析には、統計ソフトウェア「R」(R Development Core Team)の「caret」パッケージを用い、「M5」メソッドでの最適化を実行した。
表7のデータを「bacteria」として読み込ませた後、次のコマンドを実行した。
m <- train(PISA~.,data=bacteria ,method="M5")
PISAを目的変数とし、説明変数を表7の細菌全種類としてモデル木を生成させるコマンドである。46サンプルデータすべてを使用した。
実行した最適モデルをm$finalModelと結果の出力コマンドを入力して出力させたところ、以下の予測モデルを得た。
M5 unpruned model tree:
(using smoothed linear models)


Treponema.denticola <= 6.1 :
| Campylobacter.concisus <= 1.45 :
| | Actinomyces.odontolyticus <= 53.95 :
| | | Streptococcus.mutans <= 15 :
| | | | Tannerella.forsythia <= 4.9 :
| | | | | Porphyromonas.gingivalis <= 1.45 : LM1 (2/25.547%)
| | | | | Porphyromonas.gingivalis > 1.45 : LM2 (2/24.859%)
| | | | Tannerella.forsythia > 4.9 : LM3 (3/3.438%)
| | | Streptococcus.mutans > 15 : LM4 (2/7.664%)
| | Actinomyces.odontolyticus > 53.95 :
| | | Fusobacterium.nucleatum.subsp..animalis <= 2.2 : LM5 (2/3.144%)
| | | Fusobacterium.nucleatum.subsp..animalis > 2.2 : LM6 (2/2.358%)
| Campylobacter.concisus > 1.45 :
| | Capnocytophaga.ochracea <= 1.15 :
| | | Fusobacterium.periodonticum <= 2.35 : LM7 (2/4.372%)
| | | Fusobacterium.periodonticum > 2.35 : LM8 (2/1.326%)
| | Capnocytophaga.ochracea > 1.15 :
| | | Total.bacteria <= 2252.6 : LM9 (4/4.488%)
| | | Total.bacteria > 2252.6 :
| | | | Capnocytophaga.ochracea <= 1.55 : LM10 (4/2.777%)
| | | | Capnocytophaga.ochracea > 1.55 :
| | | | | Tannerella.forsythia <= 2.9 :
| | | | | | Tannerella.forsythia <= 2.35 : LM11 (2/11.398%)
| | | | | | Tannerella.forsythia > 2.35 : LM12 (2/1.572%)
| | | | | Tannerella.forsythia > 2.9 : LM13 (3/17.057%)
Treponema.denticola > 6.1 :
| Campylobacter.rectus <= 6.4 :
| | Veillonella.parvula <= 1.25 :
| | | Streptococcus.gordonii <= 36.15 : LM14 (3/12.485%)
| | | Streptococcus.gordonii > 36.15 :
| | | | Tannerella.forsythia <= 11.5 : LM15 (2/8.205%)
| | | | Tannerella.forsythia > 11.5 : LM16 (2/1.572%)
| | Veillonella.parvula > 1.25 : LM17 (3/55.578%)
| Campylobacter.rectus > 6.4 :
| | Total.bacteria <= 1916.55 : LM18 (2/4.176%)
| | Total.bacteria > 1916.55 : LM19 (2/1.965%)

LM num: 1
.outcome =
0.1154 * Total.bacteria
- 87.6234 * Porphyromonas.gingivalis
+ 32.2177 * Tannerella.forsythia
+ 3.4246 * Treponema.denticola
+ 10.3816 * Fusobacterium.nucleatum.subsp..animalis
+ 4.5488 * Streptococcus.gordonii
- 2.8479 * Actinomyces.odontolyticus
- 4.1907 * Streptococcus.mutans
+ 353.0028

LM num: 2
.outcome =
0.1154 * Total.bacteria
- 87.6234 * Porphyromonas.gingivalis
+ 32.2177 * Tannerella.forsythia
+ 3.4246 * Treponema.denticola
+ 10.3816 * Fusobacterium.nucleatum.subsp..animalis
+ 4.5488 * Streptococcus.gordonii
- 2.8479 * Actinomyces.odontolyticus
- 4.1907 * Streptococcus.mutans
+ 354.0939

LM num: 3
.outcome =
0.1154 * Total.bacteria
- 87.6234 * Porphyromonas.gingivalis
+ 32.2177 * Tannerella.forsythia
+ 3.4246 * Treponema.denticola
+ 10.3816 * Fusobacterium.nucleatum.subsp..animalis
+ 4.5488 * Streptococcus.gordonii
- 2.8479 * Actinomyces.odontolyticus
- 4.1907 * Streptococcus.mutans
+ 353.0353

LM num: 4
.outcome =
0.1154 * Total.bacteria
- 87.6234 * Porphyromonas.gingivalis
+ 32.2177 * Tannerella.forsythia
+ 3.4246 * Treponema.denticola
+ 10.3816 * Fusobacterium.nucleatum.subsp..animalis
+ 4.5488 * Streptococcus.gordonii
- 2.8479 * Actinomyces.odontolyticus
- 5.4232 * Streptococcus.mutans
+ 347.5507

LM num: 5
.outcome =
0.1154 * Total.bacteria
- 87.6234 * Porphyromonas.gingivalis
+ 32.2177 * Tannerella.forsythia
+ 3.4246 * Treponema.denticola
+ 10.3816 * Fusobacterium.nucleatum.subsp..animalis
+ 4.5488 * Streptococcus.gordonii
- 3.3806 * Actinomyces.odontolyticus
+ 263.7849

LM num: 6
.outcome =
0.1154 * Total.bacteria
- 87.6234 * Porphyromonas.gingivalis
+ 32.2177 * Tannerella.forsythia
+ 3.4246 * Treponema.denticola
+ 10.3816 * Fusobacterium.nucleatum.subsp..animalis
+ 4.5488 * Streptococcus.gordonii
- 3.3806 * Actinomyces.odontolyticus
+ 262.54

LM num: 7
.outcome =
0.1154 * Total.bacteria
- 87.6234 * Porphyromonas.gingivalis
+ 27.5903 * Tannerella.forsythia
+ 3.4246 * Treponema.denticola
+ 10.3816 * Fusobacterium.nucleatum.subsp..animalis
+ 4.5488 * Streptococcus.gordonii
- 0.6784 * Actinomyces.odontolyticus
+ 33.5894

LM num: 8
.outcome =
0.1154 * Total.bacteria
- 87.6234 * Porphyromonas.gingivalis
+ 27.5903 * Tannerella.forsythia
+ 3.4246 * Treponema.denticola
+ 10.3816 * Fusobacterium.nucleatum.subsp..animalis
+ 4.5488 * Streptococcus.gordonii
- 0.6784 * Actinomyces.odontolyticus
+ 32.468

LM num: 9
.outcome =
0.0649 * Total.bacteria
- 87.6234 * Porphyromonas.gingivalis
+ 43.8294 * Tannerella.forsythia
+ 3.4246 * Treponema.denticola
+ 10.3816 * Fusobacterium.nucleatum.subsp..animalis
+ 4.5488 * Streptococcus.gordonii
- 0.6784 * Actinomyces.odontolyticus
+ 146.6484

LM num: 10
.outcome =
0.0785 * Total.bacteria
- 87.6234 * Porphyromonas.gingivalis
+ 44.0191 * Tannerella.forsythia
+ 3.4246 * Treponema.denticola
+ 10.3816 * Fusobacterium.nucleatum.subsp..animalis
+ 4.5488 * Streptococcus.gordonii
- 0.6784 * Actinomyces.odontolyticus
+ 93.2371

LM num: 11
.outcome =
0.0785 * Total.bacteria
- 87.6234 * Porphyromonas.gingivalis
+ 45.0694 * Tannerella.forsythia
+ 3.4246 * Treponema.denticola
+ 10.3816 * Fusobacterium.nucleatum.subsp..animalis
+ 4.5488 * Streptococcus.gordonii
- 0.6784 * Actinomyces.odontolyticus
+ 91.1073

LM num: 12
.outcome =
0.0785 * Total.bacteria
- 87.6234 * Porphyromonas.gingivalis
+ 45.0694 * Tannerella.forsythia
+ 3.4246 * Treponema.denticola
+ 10.3816 * Fusobacterium.nucleatum.subsp..animalis
+ 4.5488 * Streptococcus.gordonii
- 0.6784 * Actinomyces.odontolyticus
+ 91.0077

LM num: 13
.outcome =
0.0785 * Total.bacteria
- 87.6234 * Porphyromonas.gingivalis
+ 45.1623 * Tannerella.forsythia
+ 3.4246 * Treponema.denticola
+ 10.3816 * Fusobacterium.nucleatum.subsp..animalis
+ 4.5488 * Streptococcus.gordonii
- 0.6784 * Actinomyces.odontolyticus
+ 93.4098

LM num: 14
.outcome =
-0.0279 * Total.bacteria
- 142.0103 * Porphyromonas.gingivalis
+ 9.7177 * Tannerella.forsythia
+ 5.5503 * Treponema.denticola
+ 24.5588 * Campylobacter.rectus
+ 16.8253 * Fusobacterium.nucleatum.subsp..animalis
+ 10.5169 * Streptococcus.gordonii
+ 793.003 * Veillonella.parvula
- 497.9287

LM num: 15
.outcome =
-0.0279 * Total.bacteria
- 142.0103 * Porphyromonas.gingivalis
+ 9.7177 * Tannerella.forsythia
+ 5.5503 * Treponema.denticola
+ 24.5588 * Campylobacter.rectus
+ 16.8253 * Fusobacterium.nucleatum.subsp..animalis
+ 10.4819 * Streptococcus.gordonii
+ 793.003 * Veillonella.parvula
- 489.1573

LM num: 16
.outcome =
-0.0279 * Total.bacteria
- 142.0103 * Porphyromonas.gingivalis
+ 9.7177 * Tannerella.forsythia
+ 5.5503 * Treponema.denticola
+ 24.5588 * Campylobacter.rectus
+ 16.8253 * Fusobacterium.nucleatum.subsp..animalis
+ 10.4819 * Streptococcus.gordonii
+ 793.003 * Veillonella.parvula
- 489.4669

LM num: 17
.outcome =
-0.0279 * Total.bacteria
- 142.0103 * Porphyromonas.gingivalis
+ 9.7177 * Tannerella.forsythia
+ 5.5503 * Treponema.denticola
+ 24.5588 * Campylobacter.rectus
+ 16.8253 * Fusobacterium.nucleatum.subsp..animalis
+ 10.4024 * Streptococcus.gordonii
+ 853.2173 * Veillonella.parvula
- 506.1734

LM num: 18
.outcome =
-0.0958 * Total.bacteria
- 142.0103 * Porphyromonas.gingivalis
+ 9.7177 * Tannerella.forsythia
+ 5.5503 * Treponema.denticola
+ 32.3143 * Campylobacter.rectus
+ 16.8253 * Fusobacterium.nucleatum.subsp..animalis
+ 7.3721 * Streptococcus.gordonii
+ 686.8925 * Veillonella.parvula
+ 139.0553

LM num: 19
.outcome =
-0.0958 * Total.bacteria
- 142.0103 * Porphyromonas.gingivalis
+ 9.7177 * Tannerella.forsythia
+ 5.5503 * Treponema.denticola
+ 32.3143 * Campylobacter.rectus
+ 16.8253 * Fusobacterium.nucleatum.subsp..animalis
+ 7.3721 * Streptococcus.gordonii
+ 686.8925 * Veillonella.parvula
+ 143.569

Number of Rules : 19
引き続き、p <- predict(m,newdata=bacteria)とコマンドを入力し、作成した予測モデル「m」に対して実測した各細菌量のSN比データの46サンプル分を代入し46サンプルに対応する予測PISA値「p」を得た。
実際のPISA値と予測PISA値の相関係数を算出するためcor(p,bacteria$PISA)をコマンドを実行し、予測PISA値「p」と実測のPISA値の相関係数を計算させたところ、相関係数0.9291664を得た。予測PISA値「p」と実測のPISA値の散布図を、図3に示す。このモデルに示すようにDNAチップのSN比からPISA値を予測可能であることが示された。
特に、表8に示した各細菌量を示すSN比と口腔内全体のPISA値との比較において、負の相関係数をもつ、Streptococcus mutans、Actinomyces odontolyticus、Campylobacter concisusと、正の相関係数をもつTannerella forsythia、Treponema denticola、Porphyromonas gingivalis、Fusobacterium nucleatum subsp.animalis、Streptococcus gordonii、Veillonella parvula、Campylobacter rectus、Capnocytophaga ochracea、Fusobacterium periodonticumの細菌群のSN比からPISA値を予測可能であることが示された。
実施例2に記載の表7と同一データを用いて、各細菌量のSN比をもとに機械学習の手法を用い、PISA値をモデル木で予測するモデルの作成を行った。解析には、統計ソフトウェア「R」(R Development Core Team)の「caret」パッケージ「M5」メソッドを用いた。
予測モデルの構築では46サンプルのうち、ランダムに34サンプルを抽出してクロスバリデーション法に用いた。クロスバリデーション法におけるモデル構築訓練データと検証データの比率は75:25とし学習回数は10回とした。
モデル構築後、モデル構築には全く使用していない残りの12サンプルのデータを将来の未知データとして用い、検証に使用した。
より具体的には、表7の表をデータフレーム名「bacteria」とした後、以下のコマンドにより実行した。
fitControl <- trainControl(method="CV",p=0.75, number=10)
train.index <- sample(nrow(bacteria),nrow(bacteria)*0.75)
#モデル作成用データのセット
data.train <- bacteria[train.index,]
#テスト用データのセット
data.test <- bacteria[-train.index,]

data.m5 <- train(data.train[,-1],data.train$PISA,method="M5",trControl=fitControl)
(クロスバリデーション法におけるモデル構築訓練データと検証データの比率は75:25とし学習回数は10回とした。)

data.m5$results #構築モデルの表示

#未知データでの検証結果
PISA.pred <- predict(data.m5,newdata=data.test[,-1])
#相関係数の計算
cor(PISA.pred,data.test$PISA)
#散布図の出力
plot(PISA.pred,data.test$PISA)

#トレーニングデータでの結果
PISA.predtr <- predict(data.m5,newdata=data.train[,-1])
#相関係数の計算
cor(PISA.predtr,data.train$PISA)
#散布図の出力
plot(PISA.predtr,data.train$PISA)
その結果、構築モデルは以下の結果が得られた(図4)。
M5 unpruned model tree:
(using smoothed linear models)


Tannerella.forsythia <= 3.589 :
| Capnocytophaga.ochracea <= 1.085 : LM1 (3/3.579%)
| Capnocytophaga.ochracea > 1.085 :
| | Prevotella.intermedia <= 1.135 :
| | | Total.bacteria <= 2072.974 : LM2 (2/2.188%)
| | | Total.bacteria > 2072.974 : LM3 (3/2.889%)
| | Prevotella.intermedia > 1.135 :
| | | Treponema.denticola <= 2.289 : LM4 (3/0.856%)
| | | Treponema.denticola > 2.289 : LM5 (2/4.57%)
Tannerella.forsythia > 3.589 :
| Capnocytophaga.sputigena <= 1.786 :
| | Eikenella.corrodens <= 1.067 :
| | | Fusobacterium.periodonticum <= 1.502 : LM6 (3/2.054%)
| | | Fusobacterium.periodonticum > 1.502 :
| | | | Streptococcus.mitis <= 3.972 : LM7 (2/7.05%)
| | | | Streptococcus.mitis > 3.972 :
| | | | | Aggregatibacter.actinomycetemcomitans <= 1.334 : LM8 (3/5.601%)
| | | | | Aggregatibacter.actinomycetemcomitans > 1.334 : LM9 (2/21.879%)
| | Eikenella.corrodens > 1.067 : LM10 (3/7.804%)
| Capnocytophaga.sputigena > 1.786 :
| | Tannerella.forsythia <= 12.803 :
| | | Total.bacteria <= 2244.721 : LM11 (2/6.126%)
| | | Total.bacteria > 2244.721 : LM12 (2/103.221%)
| | Tannerella.forsythia > 12.803 :
| | | Total.bacteria <= 1907.395 : LM13 (2/25.331%)
| | | Total.bacteria > 1907.395 : LM14 (2/1.945%)

LM num: 1
.outcome =
0.1964 * Total.bacteria
+ 8.473 * Treponema.denticola
- 3.7524 * Streptococcus.mitis
+ 136.0844

LM num: 2
.outcome =
0.1872 * Total.bacteria
+ 8.473 * Treponema.denticola
- 102.9193 * Prevotella.intermedia
- 3.7524 * Streptococcus.mitis
+ 297.6459

LM num: 3
.outcome =
0.1877 * Total.bacteria
+ 8.473 * Treponema.denticola
- 102.9193 * Prevotella.intermedia
- 3.7524 * Streptococcus.mitis
+ 295.2116
LM num: 4
.outcome =
0.1964 * Total.bacteria
+ 10.9568 * Treponema.denticola
- 102.9193 * Prevotella.intermedia
- 3.7524 * Streptococcus.mitis
+ 261.0986

LM num: 5
.outcome =
0.1964 * Total.bacteria
+ 11.1029 * Treponema.denticola
- 102.9193 * Prevotella.intermedia
- 3.7524 * Streptococcus.mitis
+ 261.4416

LM num: 6
.outcome =
0.1528 * Total.bacteria
+ 6.7947 * Tannerella.forsythia
+ 6.5901 * Treponema.denticola
+ 29.1856 * Aggregatibacter.actinomycetemcomitans
+ 74.536 * Capnocytophaga.sputigena
- 197.1112 * Eikenella.corrodens
- 5.6233 * Streptococcus.mitis
+ 637.0725

LM num: 7
.outcome =
0.1528 * Total.bacteria
+ 6.7947 * Tannerella.forsythia
+ 6.5901 * Treponema.denticola
+ 32.1268 * Aggregatibacter.actinomycetemcomitans
+ 74.536 * Capnocytophaga.sputigena
- 197.1112 * Eikenella.corrodens
- 5.6233 * Streptococcus.mitis
+ 625.658

LM num: 8
.outcome =
0.1528 * Total.bacteria
+ 6.7947 * Tannerella.forsythia
+ 6.5901 * Treponema.denticola
+ 32.301 * Aggregatibacter.actinomycetemcomitans
+ 74.536 * Capnocytophaga.sputigena
- 197.1112 * Eikenella.corrodens
- 5.6233 * Streptococcus.mitis
+ 627.695

LM num: 9
.outcome =
0.1528 * Total.bacteria
+ 6.7947 * Tannerella.forsythia
+ 6.5901 * Treponema.denticola
+ 32.3574 * Aggregatibacter.actinomycetemcomitans
+ 74.536 * Capnocytophaga.sputigena
- 197.1112 * Eikenella.corrodens
- 5.6233 * Streptococcus.mitis
+ 627.9978

LM num: 10
.outcome =
0.1528 * Total.bacteria
+ 6.7947 * Tannerella.forsythia
+ 6.5901 * Treponema.denticola
+ 22.9812 * Aggregatibacter.actinomycetemcomitans
+ 74.536 * Capnocytophaga.sputigena
- 273.7656 * Eikenella.corrodens
- 5.6233 * Streptococcus.mitis
+ 673.9955

LM num: 11
.outcome =
0.1989 * Total.bacteria
+ 14.1331 * Tannerella.forsythia
+ 6.5901 * Treponema.denticola
+ 90.7394 * Capnocytophaga.sputigena
- 6.2112 * Streptococcus.mitis
+ 410.322

LM num: 12
.outcome =
0.1989 * Total.bacteria
+ 14.1331 * Tannerella.forsythia
+ 6.5901 * Treponema.denticola
+ 90.7394 * Capnocytophaga.sputigena
- 6.2112 * Streptococcus.mitis
+ 409.7214

LM num: 13
.outcome =
0.1989 * Total.bacteria
+ 14.1331 * Tannerella.forsythia
+ 6.5901 * Treponema.denticola
+ 90.7394 * Capnocytophaga.sputigena
- 6.2112 * Streptococcus.mitis
+ 497.6085

LM num: 14
.outcome =
0.1989 * Total.bacteria
+ 14.1331 * Tannerella.forsythia
+ 6.5901 * Treponema.denticola
+ 90.7394 * Capnocytophaga.sputigena
- 6.2112 * Streptococcus.mitis
+ 501.0283

Number of Rules : 14
構築したモデルに対しトレーニングデータを入力した結果において、実際のPISAの値と予測値の相関係数は0.9318094となった。その結果を、図5の散布図で示した。
一方、構築したモデルに対し未知データを入力した検証結果では、実際のPISAの値と予測値の相関係数は 0.5986115となった。その結果を、図6の散布図で示した。
この結果より、未知のデータに対してもPISAの値を予測できることが示された。
特に、表8に示した各細菌量を示すSN比と口腔内全体のPISA値との比較において、負の相関係数をもつ、Streptococcus mitisと、正の相関係数をもつTannerella forsythia、 Treponema denticola、Capnocytophaga ochracea、 Capnocytophaga sputigena、Eikenella corrodens 、Aggregatibacter actinomycetemcomitans、Prevotella intermedia、Fusobacterium periodonticumの細菌群のSN比からPISA値を予測可能であることが示された。
<既存の判定方法との比較>
唾液を検体とした歯周病の状態の判定については、歯周疾患者における抗菌療法の診療ガイドライン(日本歯周病学会編集)によると、Porphyromonas gingivalis、Tannerella forsythensis及びTreponema denticolaの3種類の細菌数の合計が、総菌数に対して0.05%未満の場合を「軽度」、総菌数に対して0.05%以上0.5%未満の場合を「中程度」、総菌数に対して0.5%以上の場合を「重度」であると判定することができる。
そこで、表7に示したSN比の内、レッドコンプレックスであるPorphyromonas gingivalis、Tannerella forsythia、Treponema denticolaの3菌種のSN比の合計値と、総菌量を示すTotal bacteriaのSN比の比率を算出して既存の指標とした。この指標により、前述の「軽度」、「中程度」、「重度」を判定した。
判定した結果とPISA実測値、実施例3の結果に基づくPISA予測値、3菌のSN比、Total bacteriaのSN比を、表9にまとめた。
Figure 0007224297000033
表9に示したPISA実測値をX軸とし、既存の「3菌の比率」をY軸にしたグラフを、図7の左グラフに示した。また同様に、PISA実測値をX軸とし、PISA予測値をY軸にしたグラフを、図7の右グラフに示した。PISA実測値との決定係数は既存の方法が0.41程度、本発明による方法が0.66程度となり、既存の方法と比べて精度よく歯周病の状態を予測できることが示された。また、既存の方法で「重度」と判定されたサンプルの「PISA予測」の平均値は、1424、「中程度」と判定されたサンプルのPISA予測の平均値は、736となり従来の指標との比較もできることが示された。
表7、8のデータを元に、細菌量がPISAと「正の相関」示す上位3菌種(Porphyromonas gingivalis、Treponema denticola、Tannerella forsythia)と細菌量がPISAと「負の相関」示す上位3菌種(Streptococcus mutans、Actinomyces odontolyticus、Streptococcus mitis bv 2)を選択し、それぞれの3種ずつの細菌のSN比の合計の比(正の相関を示す細菌と負の相関を示す細菌の比:「バランス指標」とする)を算出した。
算出した値(X軸)とPISA値(Y軸)との散布図を、図8に示した。その結果、XとY関係について、決定係数0.41程度の関係性があることが分かった。すなわち、バランス指標を用いた場合においてもPISAを推定できることが分かった。
近年新たに議論されている細菌種について検討するため実施例1と同様に新たに表10に記載の細菌プローブを搭載したDNAチップを作成した。
Figure 0007224297000034
実施例1と同一のサンプルに一部追加収集した56サンプルについて、PISAのデータも収集し、それらのサンプルについて、表10に記載の新たなDNAチップで蛍光強度データを取得した。取得時の実験条件は実施例1と同様に実施したが、以下の2点のみ変更した。
PCRに使用したプライマーを以下のように変更して実施した。
R、Yは混合塩基を示しており、RはAとG、YはCとTを示す。
フォワードプライマー(細菌増幅用):
5’-Cy5-TACGGGAGGCAGCAG-3’(配列番号72)
リバースプライマー(細菌増幅用):
5’-CRGGGTATCTAATCCYGTT-3’(配列番号73)
フォワードプライマー(絶対量指標増幅用):
5’-Cy5-GAGAAGCCTACACAAACGTAACGTC-3’(配列番号34)
リバースプライマー(絶対量指標増幅用):
5’-CTCTAAAGACCGCTCTATCTCGG-3’(配列番号35)
また、ハイブリダイゼーション温度、時間を50℃で16時間とした。
ひきつづき得られた蛍光強度を以下のように処理した。
検出対象細菌用プローブを搭載したスポットの蛍光強度を、バックグラウンド値(プローブを搭載していないスポットの蛍光強度の中央値)で減算し、ハイブリダイゼーションに由来するシグナル強度を算出した。このとき、シグナル強度がある閾値を下回るものについては、ノイズと判断し「0」とする。ここでは閾値として、プローブを搭載していないスポット30個の傾向強度のうち上下位5つを除いた20個の値の標準偏差の3倍の値を用いた。
さらに、全体の細菌に対する各細菌の相対比を、検出対象細菌のプローブのシグナル強度を総菌量指標用プローブのシグナル強度で除算して算出した。この後の解析には、これら総量に対する相対比をlog10で変換したものを用いた。ただし、値が「0」のものについては計算できないため、log10変換後の数値を―4に置き換えた。こうして、56検体すべてについてデータを得た。この結果とPISAについて検体ごとにまとめたものを表11に示した。
Figure 0007224297000035
Figure 0007224297000036
Figure 0007224297000037
Figure 0007224297000038
Figure 0007224297000039
Figure 0007224297000040
PISAと各細菌のlog10(総量に対する相対比)の値の相関係数を36種類すべてについて算出し、さらにそれら相関係数の絶対値が0.2より大きい細菌種を選抜した。相関係数が正の値を示すもの、負の値を示すものを以下に示す(表12)。
Figure 0007224297000041
相関係数が負の値を示す細菌群は、Prevotella pallens,Streptococcus salivarius,Eubacterium sulci,Rothia mucilaginosa,Prevotella denticola,Veillonella atypica,Prevotella histicola,Megasphaera micronuciformis,Streptococcus parasanguinisの9菌種とした。
相関係数が正の値を示す細菌群は、SR1 sp. OT 345,Porphyromonas catoniae,Selenomonas sputigena,Neisseria flavescens,Streptococcus sobrinus,Parvimonas micra,Peptostreptococcus stomatis,Treponema socranskii,Eubacterium saphenum,Eubacterium nodatum ,Treponema medium,Filifactor alocis,Porphyromonas endodontalis
の13菌種とした。
続いて、上記の22菌種のlog10変換後の数値を説明変数とし、PISA値を目的変数として重回帰分析で予測するモデルの作成を行った。解析には、統計ソフトウェア「R」(R Development Core Team)を用い、「lm」関数を用いて実行した。
表11のPISAと当該細菌のデータ(22種類、56サンプルデータ)を「data01」として読み込ませた後、次のコマンドを実行した。
res1 <- lm(PISA~.,data=data01)
(PISAを目的変数として、22種類の細菌を説明変数として、56サンプルのデータで重回帰分析を実施するコマンド)
実行した予測モデル式res1に各細菌量のデータを代入して得られる予測PISA値と実際のPISAとの相関係数は0.8618542となった。この散布図を図9に示した。
配列番号1~74:合成DNA

Claims (11)

  1. 唾液中の2種以上の細菌の細菌量を検出し、得られた検出結果を指標として歯周ポケット炎症面積を推定する方法であって、検出される細菌は、
    該細菌の細菌量と歯周ポケット炎症面積とが正の相関関係を示す細菌と、
    該細菌の細菌量と歯周ポケット炎症面積とが負の相関関係を示す細菌とを含む、方法。
  2. 前記歯周ポケット炎症面積がPISA、又はCAPRSの値で表されるものである、請求項1に記載の方法。
  3. 正の相関関係を示す細菌が、Treponema denticola、Tannerella forsythia、Fusobacterium nucleatum subsp. animalis、Porphyromonas gingivalis、Campylobacter rectus、Fusobacterium nucleatum subsp. nucleatum、Selenomonas noxia、Veillonella parvula、Streptococcus gordonii、Fusobacterium nucleatum subsp. vincentii、Streptococcus intermedius、Capnocytophaga ochracea、Capnocytophaga sputigena、Aggregatibacter actinomycetemcomitans、Fusobacterium nucleatum subsp. polymorphum、及びFusobacterium periodonticumからなる群から選ばれる少なくとも1種である、請求項1に記載の方法。
  4. 正の相関関係を示す細菌が、SR1 sp. OT 345、Porphyromonas catoniae、Selenomonas sputigena、Neisseria flavescens、Streptococcus sobrinus、Parvimonas micra、Peptostreptococcus stomatis、Treponema socranskii、Eubacterium saphenum、Eubacterium nodatum、Treponema medium、Filifactor alocis、及びPorphyromonas endodontalisからなる群から選ばれる少なくとも1種である、請求項1に記載の方法。
  5. 負の相関関係を示す細菌が、Streptococcus mutans、Actinomyces odontolyticus、Streptococcus mitis bv 2、Streptococcus mitis、Campylobacter concisus、及びCapnocytophaga gingivalisからなる群から選ばれる少なくとも1種である、請求項1に記載の方法。
  6. 負の相関関係を示す細菌が、Prevotella pallens、Streptococcus salivarius、Eubacterium sulci、Rothia mucilaginosa、Prevotella denticola、Veillonella atypica、Prevotella histicola、Megasphaera micronuciformis、及びStreptococcus parasanguinisからなる少なくとも1種である、請求項1に記載の方法。
  7. 負の相関係数を示す細菌が、Streptococcus mutans、Actinomyces odontolyticus、Campylobacter concisus、Actinomyces.naeslundii.II、及びStreptococcus.constellatusを含み、
    正の相関係数を示す細菌が、Capnocytophaga sputigena、Tannerella forsythia、Streptococcus gordonii、Treponema denticola、Fusobacterium nucleatum subsp. nucleatum、及びPorphyromonas gingivalisを含む、
    請求項1に記載の方法。
  8. 負の相関係数を示す細菌が、Streptococcus mutans、Actinomyces odontolyticus、Campylobacter concisusを含み、
    正の相関係数を示す細菌が、Tannerella forsythia、Treponema denticola、Porphyromonas gingivalis、Fusobacterium nucleatum subsp.animalis、Streptococcus gordonii、Veillonella parvula、Campylobacter rectus、Capnocytophaga ochracea、及びFusobacterium periodonticumを含む、
    請求項1に記載の方法。
  9. 負の相関係数を示す細菌が、Streptococcus mitisを含み、
    正の相関係数をもつTannerella forsythia、Treponema denticola、Capnocytophaga ochracea、 Capnocytophaga sputigena、Eikenella corrodens 、Aggregatibacter actinomycetemcomitans、Prevotella intermedia、及びFusobacterium periodonticumを含む、
    請求項1に記載の方法。
  10. 以下の(1)~(4)工程を有する、請求項1に記載の方法。
    (1)歯周ポケット炎症面積が既知の被験者の唾液サンプルから、唾液中の各種細菌の細菌量を検出する工程
    (2)当該各種細菌の細菌量について、菌ごとに固有の歯周ポケット炎症面積との相関係数を求め、各種細菌の細菌量と歯周ポケット炎症面積との関係式を構築して予測モデルを作成する工程
    (3)歯周ポケット炎症面積が未知の被験者の唾液サンプルから、唾液中の各種細菌の細菌量を検出する工程
    (4)(3)で得られた各種細菌量を(2)で得られた関係式に入れて歯周ポケット炎症面積を推定する工程
  11. 前記予測モデルの作成方法が、線形回帰、回帰木、モデル木、ニューラルネットワーク、サポートベクターマシン、バギング、ブースティング、ランダムフォレストの機械学習アルゴリズムから選ばれる1種を使用する方法である、請求項10に記載の方法。
JP2019550508A 2017-11-02 2018-11-02 歯周ポケット炎症面積の推定方法 Active JP7224297B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2017212403 2017-11-02
JP2017212403 2017-11-02
PCT/JP2018/040916 WO2019088271A1 (ja) 2017-11-02 2018-11-02 歯周ポケット炎症面積の推定方法

Publications (2)

Publication Number Publication Date
JPWO2019088271A1 JPWO2019088271A1 (ja) 2020-11-26
JP7224297B2 true JP7224297B2 (ja) 2023-02-17

Family

ID=66332003

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019550508A Active JP7224297B2 (ja) 2017-11-02 2018-11-02 歯周ポケット炎症面積の推定方法

Country Status (7)

Country Link
US (1) US20210164028A1 (ja)
EP (1) EP3705583A4 (ja)
JP (1) JP7224297B2 (ja)
AU (1) AU2018361590A1 (ja)
BR (1) BR112020008566A2 (ja)
CA (1) CA3084010A1 (ja)
WO (1) WO2019088271A1 (ja)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7271341B2 (ja) * 2019-06-28 2023-05-11 サンスター スイス エスエー 歯周病進行リスク測定方法及びキット
WO2023133083A1 (en) * 2022-01-07 2023-07-13 CellectGen, Inc. Biofluid-based diagnostic screening system and methods of use

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070269813A1 (en) 2005-11-03 2007-11-22 Dewhirst Floyd E Methods and arrays for identifying human microflora
JP2017023093A (ja) 2015-07-27 2017-02-02 三菱レイヨン株式会社 歯周病及び齲蝕のリスクの判定方法
WO2018012011A1 (ja) 2016-07-11 2018-01-18 三菱ケミカル株式会社 口腔内検査方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS4917815B1 (ja) 1970-11-17 1974-05-04
JPS51882B1 (ja) 1970-12-30 1976-01-12
JPH11108928A (ja) 1997-10-01 1999-04-23 Dainakomu:Kk 生体高分子配列シートの製造方法
JP4404328B2 (ja) 1999-08-26 2010-01-27 三菱レイヨン株式会社 生体高分子配列薄片の製造方法
JP2004051536A (ja) 2002-07-19 2004-02-19 Univ Nihon 歯周病疾患マーカー蛋白質及びそれを利用した歯周病診断方法
JP2004163211A (ja) 2002-11-12 2004-06-10 Mitsubishi Rayon Co Ltd 中空繊維へのゲル充填方法、ゲル充填中空繊維及びそれを用いた生体関連物質マイクロアレイ
JP5043320B2 (ja) 2005-09-13 2012-10-10 三菱レイヨン株式会社 メチル化dna及び/又は非メチル化dnaの検出方法
JP5869323B2 (ja) 2011-12-01 2016-02-24 花王株式会社 歯周病マーカー
US20170211131A1 (en) * 2015-11-02 2017-07-27 General Biologicals Corporation Method for evaluating oral health
JP2017085944A (ja) * 2015-11-06 2017-05-25 ライオン商事株式会社 イヌの歯周病の評価方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070269813A1 (en) 2005-11-03 2007-11-22 Dewhirst Floyd E Methods and arrays for identifying human microflora
JP2017023093A (ja) 2015-07-27 2017-02-02 三菱レイヨン株式会社 歯周病及び齲蝕のリスクの判定方法
WO2018012011A1 (ja) 2016-07-11 2018-01-18 三菱ケミカル株式会社 口腔内検査方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
NESSE, W. et al.,Periodontal inflamed surface area: quantifying inflammatory burden,J. Clin. Periodontol.,2008年,Vol. 35,pp. 668-673
PARK, Shin-Young et al.,Periodontal inflamed surface area as a novel numerical variable describing periodontal conditions,J. Periodontal Implant Sci.,2017年10月30日,Vol. 47,pp. 328-338

Also Published As

Publication number Publication date
CN111542617A (zh) 2020-08-14
WO2019088271A1 (ja) 2019-05-09
JPWO2019088271A1 (ja) 2020-11-26
AU2018361590A1 (en) 2020-06-18
BR112020008566A2 (pt) 2020-10-06
EP3705583A1 (en) 2020-09-09
CA3084010A1 (en) 2019-05-09
US20210164028A1 (en) 2021-06-03
EP3705583A4 (en) 2021-01-20

Similar Documents

Publication Publication Date Title
JP7122337B2 (ja) 口腔内検査方法
Shahi et al. Polymerase chain reaction (PCR)-based methods: promising molecular tools in dentistry
Santos et al. Reverse transcription and polymerase chain reaction: principles and applications in dentistry
KR101706070B1 (ko) 멀티플렉스 실시간 pcr을 이용한 다수의 구강질환 원인 세균의 동시 검출용 조성물 및 이를 이용한 검출 방법
Sakellari et al. Prevalence and distribution of Aggregatibacter actinomycetemcomitans serotypes and the JP2 clone in a Greek population
US20070269813A1 (en) Methods and arrays for identifying human microflora
JP2017023093A (ja) 歯周病及び齲蝕のリスクの判定方法
JP7224297B2 (ja) 歯周ポケット炎症面積の推定方法
JP7307683B2 (ja) 臨床的指標と関連性のある細菌群の情報を利用した口腔内検査方法
US11732310B2 (en) DNA chip for detecting dental caries bacteria
CN111757942B (zh) 利用与临床指标具有相关性的细菌群的信息的口腔内检查方法
CN111542617B (zh) 牙周袋炎症面积的推定方法
JP2007222136A (ja) 歯周病菌の検出方法
KR101336862B1 (ko) 실시간 역전사효소 중합반응을 이용한 결핵의 진단을 위한 정보제공방법 및 이를 위한 결핵 진단용 키트
RU2607046C2 (ru) Способ оценки обсемененности пародонта патогенными бактериями с применением полимеразной цепной реакции в реальном времени
JP2021010343A (ja) 口腔内細菌による健康状態の予測方法
WO2008137506A2 (en) Methods and kits for dog plaque disease
Indhumathi Microbial Analysis of Orange Complex Organisms of the Whole Saliva in Patients with Gingivitis and Gingival Recession using Next Generation Sequencing
Orrù et al. Clinical Practice & Epidemiology in Mental Health

Legal Events

Date Code Title Description
A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20201008

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210730

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20220712

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20220909

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20221107

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20230117

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230207

R150 Certificate of patent or registration of utility model

Ref document number: 7224297

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150