「16SリボソームRNA」の版間の差分

削除された内容 追加された内容
空白文字の調整
(5人の利用者による、間の9版が非表示)
1行目:
[[ファイル:010_small_subunit-1FKA.gif|サムネイル|''[[サーマス・サーモフィルス|Thermus thermophilusの]]''30Sサブユニットの分子構造。タンパク質は青で、単一のRNA鎖はオレンジで示されている<ref name="Schluenzen">{{Cite journal|date=September 2000|title=Structure of functionally activated small ribosomal subunit at 3.3 angstroms resolution|journal=Cell|volume=102|issue=5|pages=615–23|DOI=10.1016/S0092-8674(00)00084-2|PMID=11007480}}</ref>。]]
{{出典の明記|date=2018年3月}}
'''16S rRNA系統解析リボソームRNA'''(16S('''16S rRNAけい''')うかいせき)は、[[リボソーム]]の小サブユニットの[[リボ核酸|RNA]][[塩基シャイン・ダルガノ配列]]を基にした[[微生物]]の進化系統を明らか(Shine-Dalgarno sequence)結合する方法の一つである。[[核生物]]の場合は18S [[リボソームRNA|rRNA]]で'''リボソーム30S小サブユニットのコンポーネントである。このRNA系統解析'''をコードする遺伝子は16S ('S'mall 'S'ub 'U'nit-rRNA、'''SSU-rRNA''') 遺伝子と呼ばれることもある。
 
16S rRNA遺伝子は、リボソームという生物の本質に関わる機能を持つRNAであるため、配列の保存性が高く、[[細菌]]や[[古細菌]]といった[[原核生物]]の間で高度に保存されている。そして、機能変化に伴う遺伝子の[[変異]]がこれからも起きる可能性が極めて低い。すなわち、遺伝子配列の[[進化]]速度が遅いことから、信頼できる[[分子時計]]として利用できる。また、遺伝子の長さが適当に長く(16S rRNAの場合、1600[[塩基対]]程度)、系統解析を行う上で十分な情報量を持つ。さらに、比較的変異しやすい部位も存在し、近縁な種でも比較が可能である。これらの特徴から、特に[[微生物]][[系統学]]の分野において、この遺伝子配列は系統進化解析によく利用されている<ref name="Weisburg">{{Cite journal|date=January 1991|title=16S ribosomal DNA amplification for phylogenetic study|journal=Journal of Bacteriology|volume=173|issue=2|pages=697–703|DOI=10.1128/jb.173.2.697-703.1991|PMID=1987160|PMC=207061}}</ref><ref name="pmid14612235">{{Cite journal|date=November 2003|title=Intragenomic heterogeneity between multiple 16S ribosomal RNA operons in sequenced bacterial genomes|journal=FEMS Microbiology Letters|volume=228|issue=1|pages=45–9|DOI=10.1016/S0378-1097(03)00717-1|PMID=14612235}}</ref><ref>{{Cite journal|date=August 2017|title=Comparative RNA function analysis reveals high functional similarity between distantly related bacterial 16 S rRNAs|journal=Scientific Reports|volume=7|issue=1|pages=9993|language=En|bibcode=2017NatSR...7.9993T|DOI=10.1038/s41598-017-10214-3|PMID=28855596|PMC=5577257}}</ref><ref name="woese1977">{{Cite journal|date=November 1977|title=Phylogenetic structure of the prokaryotic domain: the primary kingdoms|journal=Proceedings of the National Academy of Sciences of the United States of America|volume=74|issue=11|pages=5088–90|bibcode=1977PNAS...74.5088W|DOI=10.1073/pnas.74.11.5088|PMID=270744|PMC=432104}}{{オープンアクセス}}</ref>。[[カール・ウーズ]]と{{仮リンク|ジョージ・E・フォックス|en|George E. Fox|label=}}が1977年に、系統学に16S rRNAを導入した<ref name="Woese_1990">{{Cite journal|date=June 1990|title=Towards a natural system of organisms: proposal for the domains Archaea, Bacteria, and Eucarya|journal=Proceedings of the National Academy of Sciences of the United States of America|volume=87|issue=12|pages=4576–9|bibcode=1990PNAS...87.4576W|DOI=10.1073/pnas.87.12.4576|PMID=2112744|PMC=54159}}</ref><ref name="woese19772">{{Cite journal|date=November 1977|title=Phylogenetic structure of the prokaryotic domain: the primary kingdoms|journal=Proceedings of the National Academy of Sciences of the United States of America|volume=74|issue=11|pages=5088–90|bibcode=1977PNAS...74.5088W|DOI=10.1073/pnas.74.11.5088|PMID=270744|PMC=432104}}{{オープンアクセス}}</ref>。[[真核生物]]の場合は18S [[リボソームRNA|rRNA]]なので、まとめて'''リボソーム小サブユニットRNA'''(<u>S</u>mall <u>S</u>ub <u>U</u>nit rRNA、'''SSU rRNA''')系統解析と呼ばれることもある。
== 歴史 ==
従来[[原核生物]]の分類は[[細胞]]の形態、分離の条件、[[染色 (生物学)|染色法]]などで行っていたが、こうした[[表現型]]の形質では[[系統樹]]上の上下関係を説明するには至らなかった。しかし[[1970年代]]、[[チトクローム]]、[[フェレドキシン]]、5S rRNAなどの[[分子系統学|塩基配列を基にした系統分類]]が[[分子生物学]]の発展とともに徐々に活発化してきた。
 
== 分子生物学的機能 ==
[[遺伝子]]の一次構造に基づく系統分類は原核生物に対して特に有効であった。[[カール・ウーズ]]らはリボソーム小サブユニットを構成するRNA、つまり16S rRNAの塩基配列を用いて原核生物の系統分類を行い、原核生物が[[真正細菌]]と[[古細菌]]という2つの[[ドメイン (分類学)|ドメイン]]からなることを証明した([[1977年]]当時はオリゴヌクレオチドカタログ法を用いた)。
16S rRNAは23S rRNAと相互作用する。このRNA複合体は構造上、[[リボソーム]]タンパク質の位置を決める足場として機能する役割を持ち、2つのリボソームサブユニット(50Sおよび30S)の結合を支援する。3 '末端には、[[伝令RNA|mRNA]]のAUG[[開始コドン]]の上流に結合するShine-Dalgarno配列の相補鎖が含まれている。16S RNAの3 '末端は、[[タンパク質]]合成の開始に関与するS1およびS21タンパク質に結合する<ref>{{Cite journal|last=Czernilofsky|first=A. P.|last2=Kurland|first2=C. G.|last3=Stöffler|first3=G.|date=1975|title=30S Ribosomal proteins associated with the 3′-terminus of 16S RNA|journal=FEBS Letters|volume=58|issue=1|pages=281–284|DOI=10.1016/0014-5793(75)80279-1|ISSN=0014-5793|PMID=1225593}}</ref>。
 
1492残基および1493残基で[[アデニン]](A)が並んでいる箇所をAサイトと呼ぶ、このAサイトのアデニンが持つN1原子とmRNA骨格の2つのOH基の間に[[水素結合]]を形成し、正確なコドン-アンチコドンのペアリングを安定化している。
現在、16S rRNAを用いた系統解析は、系統樹の作成のみならず、任意の環境中における細菌・古細菌の群集構造の観測に役立っている。この方法を用いると、分離・培養が困難な難培養性の菌種を含めて網羅的に群集構造を明らかにできる他、新規の菌の存在を配列解析から明らかにする事ができる(1996年のBarnsによる;[[メタゲノミクス]]参照)。
[[ファイル:16S_domain_structure.svg|右|サムネイル|835x835ピクセル|16S rRNAの二次構造<ref>{{Cite pmid|8811093}}</ref>。文字はすべての原核生物の保存的ヌクレオチドを示し、アスタリスクは細菌・古細菌において保存的なヌクレオチドを示す。他のすべてのヌクレオチドはドットで示される。ドメインIは5 '末端ドメイン、ドメインIIは中央ドメイン、ドメインIIIは大きな3'末端ドメイン、ドメインIVは小さな3 '末端ドメインに対応する。]]
 
== 超可変領域 ==
== 16S rRNA塩基配列が系統解析に適している点 ==
細菌の16S rRNA遺伝子には、リボソーム小サブユニットの二次構造に関与する、9つの超可変領域(V1-V9)が含まれており、これらの長さは約30-100[[塩基対]]である<ref name="pmid6462918">{{Cite journal|date=July 1984|title=On the evolutionary descent of organisms and organelles: a global phylogeny based on a highly conserved structural core in small subunit ribosomal RNA|journal=Nucleic Acids Research|volume=12|issue=14|pages=5837–52|DOI=10.1093/nar/12.14.5837|PMID=6462918|PMC=320035}}</ref>。保存の程度は超可変領域間で大きく異なり、より保存された領域は門や綱といったより高レベルの分類法に利用でき、一方で保存度の低い領域は属や種といったより低レベルの分類に利用される<ref name=":0">{{Cite journal|date=March 2016|title=Sensitivity and correlation of hypervariable regions in 16S rRNA genes in phylogenetic analysis|journal=BMC Bioinformatics|volume=17|issue=1|pages=135|DOI=10.1186/s12859-016-0992-y|PMID=27000765|PMC=4802574}}</ref>。16S rRNA配列全体をシーケンスすることで全超可変領域の比較が可能になるが、16S rRNAは約1,500塩基の長さを持つため、多様な細菌群集を満遍なくシーケンスするには費用がかかってしまう<ref name=":0" />。そのため、細菌叢解析のような研究では通常、Illumina社製のゲノムシーケンス技術を利用しており、454パイロシーケンスや[[サンガー法|サンガーシーケンス]]よりもそれぞれ約50倍、12,000倍ほど安価にシーケンスすることができる<ref>{{Cite journal|date=June 2011|title=Generation of multimillion-sequence 16S rRNA gene libraries from complex microbial communities by assembling paired-end illumina reads|journal=Applied and Environmental Microbiology|volume=77|issue=11|pages=3846–52|DOI=10.1128/AEM.02772-10|PMID=21460107|PMC=3127616}}</ref>。しかしながら、Illumina社製シーケンサーでは75〜250塩基(Illumina MiSeqでは最大300塩基)のリード長しか得られないため、細菌叢サンプルから16S rRNA遺伝子配列を完璧に組み立てることはできない<ref name=":1">{{Cite journal|date=2016-09-20|title=A method for high precision sequencing of near full-length 16S rRNA genes on an Illumina MiSeq|journal=PeerJ|volume=4|pages=e2492|DOI=10.7717/peerj.2492|PMID=27688981|PMC=5036073}}</ref>。一方で、超可変領域はその短さのため、Illuminaシーケンサを1回実行するだけで配列解析を行えるため、この超可変領域は菌叢解析における理想的なターゲットになっている<ref name=":1" />。
# リボソームという生物の本質に関わる機能を持ったRNAなので配列の保存性が高く、極めて関係の遠い生物同士でも配列の比較が可能である。
# 真核生物、原核生物問わずすべての種に存在し、機能変化に伴う遺伝子の[[変異]]がこれからも起きる可能性が極めて少ない。
# ゲノム内にコピーが複数個存在しても、塩基配列にほとんど差が無い。
# 遺伝子の長さが適当に長く(16S rRNAの場合、1600[[塩基対]]程度)、系統解析に十分な情報量を持つ。
# 比較的変異しやすい部位も存在し、近縁な種でも比較が可能である。
# 細胞内に大量に存在し、[[ポリメラーゼ連鎖反応|PCR]]の開発がなされる以前から塩基配列の比較が可能であった。
# 全生物にわたって完全に保存された部位が三箇所ほど存在し、そうした[[プライマー]]('''ユニバーサルプライマー''')を設計することにより塩基配列の決定が容易である。
以上の理由から、16S rRNAは微生物のみならず、最近は真核生物の系統分類にも使用されている。
なお、系統樹作成の際は、16S rRNA塩基配列のみならず、ほかの遺伝子も比較して構築していくことが一般的である。
 
16S rRNA超可変領域は、細菌系統間で大きく配列が異なる場合があるが、全体としては16S rRNA遺伝子は真核生物([[ 18SリボソームRNA|18SリボソームRNA]])よりも良く均一性を維持しているため、[[シーケンスアラインメント|アライメント]]が比較的容易である<ref name="pmid8811093">{{Cite journal|date=September 1996|title=A quantitative map of nucleotide substitution rates in bacterial rRNA|journal=Nucleic Acids Research|volume=24|issue=17|pages=3381–91|DOI=10.1093/nar/24.17.3381|PMID=8811093|PMC=146102}}</ref>。さらに、16S rRNA遺伝子には超可変領域間の高度に保存された配列が含まれているため、異なる[[タクソン|分類群]]にわたって同じ超可変領域を確実に[[ポリメラーゼ連鎖反応|PCR]]増幅できるユニバーサルプライマーの設計が可能である<ref name=":2">{{Cite journal|date=2013-02-27|title=The variability of the 16S rRNA gene in bacterial genomes and its consequences for bacterial community analyses|journal=PLOS ONE|volume=8|issue=2|pages=e57923|bibcode=2013PLoSO...857923V|DOI=10.1371/journal.pone.0057923|PMID=23460914|PMC=3583900}}</ref>。すべての細菌系統を[[ドメイン (分類学)|ドメイン]]から種に渡って正確に分類できる超可変領域は存在しないが、特定の分類レベルをほぼ確実に予測できるものもは知られている<ref name=":02">{{Cite journal|date=March 2016|title=Sensitivity and correlation of hypervariable regions in 16S rRNA genes in phylogenetic analysis|journal=BMC Bioinformatics|volume=17|issue=1|pages=135|DOI=10.1186/s12859-016-0992-y|PMID=27000765|PMC=4802574}}</ref>。多くの菌叢解析研究では、完全な16S rRNA遺伝子と同程度の正確性で門レベルの系統解析を行うことができるV4超可変領域を選択することが多い<ref name=":02" />。保存度の低い領域は、高次の系統分類には不向きであるが、例えば特定の病原体を検出するような用途でよく利用される。2007年にChakravortyらが発表した研究では、どの超可変領域が疾患特異的かつ広範なアッセイに利用できるかを調べ、さまざまな病原体のV1-V8領域を報告している<ref name=":3">{{Cite journal|date=May 2007|title=A detailed analysis of 16S ribosomal RNA gene segments for the diagnosis of pathogenic bacteria|journal=Journal of Microbiological Methods|volume=69|issue=2|pages=330–9|DOI=10.1016/j.mimet.2007.02.005|PMID=17391789|PMC=2562909}}</ref>。また他の研究では、病原体の属の特定にはV3領域を利用することが最適であり、[[炭疽症|炭疽菌]]を含むテストされたすべての[[米国の届出可能な疾患|CDC監視病原体]]においてはV6領域が種の区別に最も高い正確性を示した、と報告されている<ref name=":3" />。
== 実際の用法 ==
16S rRNAを利用する際は、ユニバーサルプライマーを用いて[[ポリメラーゼ連鎖反応|PCR]]による増幅を行い、得られた増幅産物をクローニングした後に[[シークエンス]]反応を行う方法が一般的である。ただ、最近はシークエンシング反応を行わなくても群集構造の解析が可能な[[DGGE]](Denaturing Gradient Gel Electrophoresis)、顕微鏡で直接観察できる[[蛍光 in situ ハイブリダイゼーション|FISH]]などの広い応用範囲がある。かつては[[制限酵素]]を用いた[[RFLP]]などが使用されていたが、現在はDGGEに取って代わられつつある。
 
16S rRNA超可変領域をベースとした配列解析は、細菌系統の分類学的研究にとって有用であるが、ごく近縁の種同士を区別することは困難な場合がある<ref name=":22">{{Cite journal|date=2013-02-27|title=The variability of the 16S rRNA gene in bacterial genomes and its consequences for bacterial community analyses|journal=PLOS ONE|volume=8|issue=2|pages=e57923|bibcode=2013PLoSO...857923V|DOI=10.1371/journal.pone.0057923|PMID=23460914|PMC=3583900}}</ref>。例えば''[[腸内細菌科]]''、''[[クロストリジウム科|クロストリジウム]][[腸内細菌科|科]]''、および''[[ペプトストレプトコッカス科|ペプトストレプトコッカス]]''科では、種間で16S rRNA遺伝子全体の最大99%の配列類似性をもつことが知られている<ref name=":4">{{Cite journal|date=2016-01-01|title=Characterization of the Gut Microbiome Using 16S or Shotgun Metagenomics|journal=Frontiers in Microbiology|volume=7|pages=459|DOI=10.3389/fmicb.2016.00459|PMID=27148170|PMC=4837688}}</ref>。この場合、種間差異はV4配列中のほんの数塩基にしか出現しないため、特に低レベルの分類において、参照データベースに基づく手法では確実に分類することが困難である<ref name=":4" />。また、利用する超可変領域の数を絞るほど、近縁な分類群の違いを観察できなくなり、サンプル全体の多様性の過小評価に繋がりうる<ref name=":22" />。さらに、細菌のゲノムは、多様なV1、V2、V6領域の配列を持つ複数の16S rRNA遺伝子をマルチコピーで保持する場合がある<ref name="pmid146122352">{{Cite journal|date=November 2003|title=Intragenomic heterogeneity between multiple 16S ribosomal RNA operons in sequenced bacterial genomes|journal=FEMS Microbiology Letters|volume=228|issue=1|pages=45–9|DOI=10.1016/S0378-1097(03)00717-1|PMID=14612235}}</ref>。これらの理由から、16S rRNAの超可変領域に基づく解析は、細菌種を分類する完璧な方法とまでは言えない。しかしながらこのような欠点がありつつも、現実的には細菌群集研究に利用できる最も有用なツールの1つとして今日利用されている<ref name=":4" />。
== 問題点 ==
16S rRNAクローンを用いた系統解析は、上記の7点に述べたような理由から正確であるとしているが、計算機の機能上昇に伴い16S rRNA以上の情報量を持つ遺伝子を用いている事もある。その一つが'''23S rRNA'''であり、この遺伝子を用いて作成した系統樹は極めて正確であるといわれている。
 
== PCRと配列シーケンシング ==
ただし、ユニバーサルプライマーを用いた群集構造解析は、環境中に存在している16S rRNAをすべて増幅してしまうために、生存個体のみを抽出しているわけではない。つまり、死亡して溶菌したようなRNAの残骸をも増幅しているということである。最近では一般的に弱いとされるRNAは環境中に長時間存在しているのではないかという研究結果も出ている。
 
=== PCR増幅 ===
16S rRNA配列を解析する際は、ユニバーサルプライマーを用いて[[ポリメラーゼ連鎖反応|PCR]]による増幅を行い、得られた増幅産物を[[シークエンス|シーケンス]]する方法が一般的である。シークエンシング反応を行わなくても群集構造の解析が可能な[[DGGE]] (Denaturing Gradient Gel Electrophoresis) 法や、顕微鏡で直接観察できる[[蛍光 in situ ハイブリダイゼーション|FISH]]法などの広い応用範囲も知られている。かつては[[制限酵素]]を用いた[[RFLP]]などが使用されていた。
 
最も一般的なプライマーペアは、Weisburg''らによって''考案された、27F-1492Rと呼ばれているセットである<ref name="Weisburg2">{{Cite journal|date=January 1991|title=16S ribosomal DNA amplification for phylogenetic study|journal=Journal of Bacteriology|volume=173|issue=2|pages=697–703|DOI=10.1128/jb.173.2.697-703.1991|PMID=1987160|PMC=207061}}</ref>''。''一部のアプリケーションでは、より短いアンプリコンが必要になる場合があり、たとえばチタンケミストリーを使用した454シーケンスでは、V1からV3をカバーするプライマーペア27F-534Rがよく選択される<ref>http://www.hmpdacc.org/tools_protocols.php#sequencing {{Webarchive|url=https://web.archive.org/web/20101030101644/http://hmpdacc.org/tools_protocols.php|date=2010-10-30}}</ref>。また、27Fではなく8Fが使用される場合も多い。この2つのプライマーはほぼ同じであるが、27FはCではなくMを持っている。
{| class="wikitable"
|+主なプライマー配列
!プライマー名
!シーケンス(5′–3 ′)
!{{Abbr|Ref.|Reference}}
|-
|8F
|AGA GTT TGA TCC TGG CTC AG
|<ref name="Eden2">{{Cite journal|date=April 1991|title=Phylogenetic analysis of Aquaspirillum magnetotacticum using polymerase chain reaction-amplified 16S rRNA-specific DNA|journal=International Journal of Systematic Bacteriology|volume=41|issue=2|pages=324–5|DOI=10.1099/00207713-41-2-324|PMID=1854644}}</ref> <ref name="James_G2">{{Cite book|title=PCR for Clinical Microbiology|first=Greg|last=James|date=15 May 2018|publisher=Springer, Dordrecht|pages=209–214|doi=10.1007/978-90-481-9039-3_28|chapter=Universal Bacterial Identification by PCR and DNA Sequencing of 16S rRNA Gene|isbn=978-90-481-9038-6}}</ref>
|-
|27F
|AGA GTT TGA TC '''M''' TGG CTC AG
|
|-
|U1492R
|GGT TAC CTT GTT ACG ACT T
|<ref name="Eden2" /> <ref name="James_G2" />
|-
|928F
|TAA AAC TYA AAK GAA TTG ACG GG
|<ref name="Weirdner2">{{Cite journal|date=March 1996|title=Diversity of uncultured microorganisms associated with the seagrass Halophila stipulacea estimated by restriction fragment length polymorphism analysis of PCR-amplified 16S rRNA genes|url=http://aem.asm.org/cgi/reprint/62/3/766.pdf|journal=Applied and Environmental Microbiology|volume=62|issue=3|pages=766–71|PMID=8975607|PMC=167844}}</ref>
|-
|336R
|ACT GCT GCS YCC CGT AGG AGT CT
|<ref name="Weirdner2" />
|-
|1100F
|YAA CGA GCG CAA CCC
|
|-
|1100R
|GGG TTG CGC TCG TTG
|
|-
|337F
|GAC TCC TAC GGG AGG CWG CAG
|
|-
|907R
|CCG TCA ATT CCT TTR AGT TT
|
|-
|785F
|GGA TTA GAT ACC CTG GTA
|
|-
|805R
|GAC TAC CAG GGT ATC TAA TC
|
|-
|533F
|GTG CCA GCM GCC GCG GTA A
|
|-
|518R
|GTA TTA CCG CGG CTG CTG G
|
|-
|1492R
|CGG TTA CCT TGT TAC GAC TT
|<ref name="HC_Jiang2">{{Cite journal|date=June 2006|title=Microbial diversity in water and sediment of Lake Chaka, an athalassohaline lake in northwestern China|journal=Applied and Environmental Microbiology|volume=72|issue=6|pages=3832–45|DOI=10.1128/AEM.02869-05|PMID=16751487|PMC=1489620}}</ref>
|}
 
=== NGSへの応用 ===
16S rRNA遺伝子配列には、高度に保存されたプライマー結合部位に加えて、複数の超可変領域が含まれており、この領域の塩基配列を利用することで、細菌の系統的な同定を行うことができる<ref>{{Cite journal|date=December 2010|title=Identification of species by multiplex analysis of variable-length sequences|journal=Nucleic Acids Research|volume=38|issue=22|pages=e203|DOI=10.1093/nar/gkq865|PMID=20923781|PMC=3001097}}</ref><ref>{{Cite journal|date=June 1999|title=Ribosomal DNA sequencing as a tool for identification of bacterial pathogens|journal=Current Opinion in Microbiology|volume=2|issue=3|pages=299–305|DOI=10.1016/S1369-5274(99)80052-6|PMID=10383862}}</ref>。現在、16S rRNA遺伝子配列シーケンシングは、[[表現型]]をベースとした細菌同定法に代わる迅速で安価な代替手法として、医学の分野で広く普及している<ref>{{Cite journal|date=October 2004|title=Impact of 16S rRNA gene sequence analysis for identification of bacteria on clinical microbiology and infectious diseases|journal=Clinical Microbiology Reviews|volume=17|issue=4|pages=840–62, table of contents|DOI=10.1128/CMR.17.4.840-862.2004|PMID=15489351|PMC=523561}}</ref>。また、細菌の識別のみならず、完全に新種な系統の発見や系統関係の再分類にも利用されている<ref>{{Cite journal|date=July 2009|title=Reverse transcription of 16S rRNA to monitor ribosome-synthesizing bacterial populations in the environment|journal=Applied and Environmental Microbiology|volume=75|issue=13|pages=4589–98|DOI=10.1128/AEM.02970-08|PMID=19395563|PMC=2704851}}</ref><ref name="Weisburg3">{{Cite journal|date=January 1991|title=16S ribosomal DNA amplification for phylogenetic study|journal=Journal of Bacteriology|volume=173|issue=2|pages=697–703|DOI=10.1128/jb.173.2.697-703.1991|PMID=1987160|PMC=207061}}</ref><ref>{{Cite journal|date=January 1998|title=Burkholderia thailandensis sp. nov., a Burkholderia pseudomallei-like species|journal=International Journal of Systematic Bacteriology|volume=48 Pt 1|issue=1|pages=317–20|DOI=10.1099/00207713-48-1-317|PMID=9542103}}</ref>。未培養系統の新種記載においても利用される<ref>{{Cite book|title=Phylogenetic identification of uncultured pathogens using ribosomal RNA sequences|volume=235|pages=[https://archive.org/details/bacterialpathoge0000unse_f8b6/page/205 205–222]|year=1994|pmid=7520119|doi=10.1016/0076-6879(94)35142-2|isbn=978-0-12-182136-4|series=Methods in Enzymology|url=https://archive.org/details/bacterialpathoge0000unse_f8b6/page/205}}</ref><ref>{{Cite journal|date=November 1996|title=Phylogenetic analysis of the bacterial communities in marine sediments|journal=Applied and Environmental Microbiology|volume=62|issue=11|pages=4049–59|PMID=8899989|PMC=168226}}</ref>。次世代シーケンシング技術を活用することで、数千の16S rRNA配列を数時間程度で解析することが可能になっており、たとえば[[腸内細菌|腸内細菌叢の]][[メタゲノム]]研究などに利用されている<ref name="pmid25226019">{{Cite journal|date=August 2014|title=Next-generation sequencing of 16S ribosomal RNA gene amplicons|journal=Journal of Visualized Experiments|volume=|issue=90|DOI=10.3791/51709|PMID=25226019|PMC=4828026}}</ref>。
 
== 解析における注意点 ==
細菌が持つ16S rRNA遺伝子配列は一つとは限らず、複数の16S rRNA遺伝子がゲノム中にマルチコピーで含まれることが多い<ref name="pmid170717872">{{Cite journal|date=January 2007|title=Use of 16S rRNA and rpoB genes as molecular markers for microbial ecology studies|journal=Applied and Environmental Microbiology|volume=73|issue=1|pages=278–88|DOI=10.1128/AEM.01177-06|PMID=17071787|PMC=1797146}}</ref>。また例外として、一部の[[好熱菌|好熱性]][[古細菌]] (例えば[[テルモプロテウス目|Thermoproteales]]など)には、16S rRNA遺伝子中に[[イントロン]]が含まれており、ユニバーサル[[プライマー (生物)|プライマー]]のアニーリングに影響を与える可能性がある<ref name="Jay">{{Cite journal|date=July 2015|title=The distribution, diversity, and importance of 16S rRNA gene introns in the order Thermoproteales|journal=Biology Direct|volume=10|issue=35|pages=35|DOI=10.1186/s13062-015-0065-6|PMID=26156036|PMC=4496867}}</ref>。また、サンプル中の真核生物に由来する[[ミトコンドリア]]や[[葉緑体]]が持つ16S rRNAもPCRで増幅されることがある。また、ユニバーサルプライマーを用いた群集構造解析は、環境中に存在している16S rRNAをすべて増幅してしまうために、生存個体のみならず死亡して溶菌したようなRNAの残骸をも増幅しうる。
 
=== 16S rRNA遺伝子の交雑 ===
進化が[[垂直感染|垂直伝達]]によって駆動されるという仮定の下では、16S rRNA遺伝子は種特異的であるみなすことができ、[[原核生物]]間の系統関係を推測する確実な遺伝的マーカーであると長年考えられてきた。しかしながら、研究が進むに連れ、これらの遺伝子においても遺伝子[[遺伝子の水平伝播|水平伝達]]が発生していることが分かってきた。このような遺伝子の転移性は、特別な[[大腸菌]]''の''遺伝子システムを用いた実験的によって確認もされている。すなわち、''大腸菌が本来持つ16S rRNA遺伝子を欠失させ''、''大腸菌''とは門レベルで系統が異なる生物種由来の外来16S rRNA遺伝子を導入したところ、変異株として増殖することが示された<ref>{{Cite journal|date=November 2012|title=Mutational robustness of 16S ribosomal RNA, shown by experimental horizontal gene transfer in Escherichia coli|journal=Proceedings of the National Academy of Sciences of the United States of America|volume=109|issue=47|pages=19220–5|DOI=10.1073/pnas.1213609109|PMID=23112186|PMC=3511107}}</ref><ref>{{Cite journal|date=August 2017|title=Comparative RNA function analysis reveals high functional similarity between distantly related bacterial 16 S rRNAs|journal=Scientific Reports|volume=7|issue=1|pages=9993|DOI=10.1038/s41598-017-10214-3|PMID=28855596|PMC=5577257}}</ref>。このような機能的互換性は、''[[サーマス・サーモフィルス|サーマスサーモフィルス]]で''も確認されている<ref>{{Cite journal|date=August 2019|title=Occurrence of randomly recombined functional 16S rRNA genes in Thermus thermophilus suggests genetic interoperability and promiscuity of bacterial 16S rRNAs|journal=Scientific Reports|volume=9|issue=1|pages=11233|DOI=10.1038/s41598-019-47807-z|PMID=31375780|PMC=6677816}}</ref>。さらに、''T. thermophilus''では、完全な遺伝子導入と部分的な遺伝子導入の両方が観察された。部分的な導入としては、宿主と外来細菌の遺伝子間で明らかにランダムな[[キメラ]]が自発的に生成されていた。したがって、16S rRNA遺伝子は、垂直遺伝と水平遺伝子伝達を含む複数のメカニズムを通じて進化している可能性があり、特に後者については今まで考えられてきたよりもはるかに高い頻度で発生している可能性がある。
 
== 16S rRNA配列データベース ==
16S rRNA遺伝子はほぼ全ての微生物に存在し、適当に配列変化が起きるため、微生物の系統分類と同定に利用されてきた<ref name="pmid25118885">{{Cite journal|date=September 2014|title=Uniting the classification of cultured and uncultured bacteria and archaea using 16S rRNA gene sequences|url=https://www.researchgate.net/publication/264798721|journal=Nature Reviews. Microbiology|volume=12|issue=9|pages=635–45|DOI=10.1038/nrmicro3330|PMID=25118885}}</ref>。ほとんどの細菌および古細菌のタイプ株が持つ16S rRNA遺伝子の配列情報は、[[アメリカ国立生物工学情報センター|NCBI]]などの公共データベースから入手できる。ただし、これらのデータベースに格納された配列は、品質が検証されていないことがよくある。そのため、16S rRNA配列のみを収集する2次データベースが広く使用されている。使用されるケースが多い有名なデータベースは以下のとおりである。
 
=== EzBioCloud ===
EzBioCloudデータベースは、以前は[[ EzTaxonデータベース|EzTaxon]]と呼ばれていた。2018年9月の時点で15,290の有効な公開名を含む62,988の細菌と古細菌の系統を含んでおり、完全な階層分類システムで構成されている。最尤推定やOrthoANIなどに基づいた系統関係に基づいて、すべての種/亜種が少なくとも1つの16S rRNA遺伝子配列によって表されている。EzBioCloudデータベースは体系的に管理されており、定期的に更新されている。新しい候補種が登録されることもある。さらにWebサイト上では、ANIの計算やContEst16S、QIIMEおよびMothurパイプライン用の16S rRNA DBといったバイオインフォマティクスツールを提供している<ref>Yoon, S. H., Ha, S. M., Kwon, S., Lim, J., Kim, Y., Seo, H. and Chun, J. (2017). Introducing EzBioCloud: A taxonomically united database of 16S rRNA and whole genome assemblies. Int J Syst Evol Microbiol. 67:1613–1617</ref>。
 
=== Ribosomal Database Project ===
<u>R</u>ibosomal <u>D</u>atabase <u>P</u>roject(RDP)は、関連するツール群やサービスと共にリボソームデータを提供するキュレーションデータベースである。系統的に纏められたリボソームRNA(rRNA)配列のアライメントや系統樹、rRNA二次構造図、およびアライメントの分析や表示をするさまざまなソフトウェアパッケージを提供している<ref>Larsen N, Olsen GJ, Maidak BL, McCaughey MJ, Overbeek R, Macke TJ, Marsh TL, Woese CR. (1993) The ribosomal database project. ''Nucleic Acids Res.'' Jul 1;21(13):3021-3.</ref>。
 
=== SILVA ===
[[ SILVAリボソームRNAデータベース|SILVA]]は、小サブユニット(16S/18S、SSU)と大サブユニット(23S/28S、LSU)のリボソームRNA配列を包括的に纏めたデータベースである。キュレーションを経たデータセットを定期的に更新している<ref>Elmar Pruesse, Christian Quast, Katrin Knittel, Bernhard M. Fuchs, Wolfgang Ludwig, Jörg Peplies, Frank Oliver Glöckner (2007) Nucleic Acids Res. SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. December; 35(21): 7188–7196.</ref>。
 
=== GreenGenes ===
Greengenesは、品質管理された包括的な16Sリファレンスデータベースである。''de novo''系統に基づいて分類されており、標準的な操作上の分類単位を提供する。現在は積極的に維持されておらず、最後の更新は2013年である<ref>{{Cite journal|date=July 2006|title=Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB|journal=Applied and Environmental Microbiology|volume=72|issue=7|pages=5069–72|DOI=10.1128/aem.03006-05|PMID=16820507|PMC=1489311}}</ref><ref>{{Cite journal|date=March 2012|title=An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea|journal=The ISME Journal|volume=6|issue=3|pages=610–8|DOI=10.1038/ismej.2011.139|PMID=22134646|PMC=3280142}}</ref>。
 
== 歴史 ==
従来[[原核生物]]の分類は[[細胞]]の形態、分離の条件、[[染色 (生物学)|染色法]]などで行っていたが、こうした[[表現型]]の形質では[[系統樹]]上の上下関係を説明するには至らなかった。しかし[[1970年代]]、[[シトクロム]]、[[フェレドキシン]]、5S rRNAなどの塩基配列を基にした系統分類が[[分子生物学]]の発展とともに徐々に活発化してきた。
 
[[遺伝子]]の一次構造に基づく系統分類は原核生物に対して特に有効であった。[[カール・ウーズ]]らはリボソーム小サブユニットを構成するRNA、つまり16S rRNAの塩基配列を用いて原核生物の系統分類を行い、原核生物が[[真正細菌]]と[[古細菌]]という2つの[[ドメイン (分類学)|ドメイン]]からなることを示した([[1977年]]当時はオリゴヌクレオチドカタログ法を用いた)。
 
現在、16S rRNAを用いた系統解析は、系統樹の作成のみならず、任意の環境中における細菌・古細菌の群集構造の観測に役立っている。この方法を用いると、分離・培養が困難な難培養性の菌種を含めて網羅的に群集構造を明らかにできる他、新規の菌の存在を配列解析から明らかにする事ができる(1996年のBarnsによる;[[メタゲノミクス]]参照)。
 
== 参考文献 ==
誤差要因としては基本的な進化速度の違いのほか、塩基配列をそのまま読み込むため、GC含率の違いに比較的強い影響を受ける。これは真核生物の初期進化や細菌の系統樹において大きな問題となった。
{{Reflist|30em}}
 
{{DEFAULTSORT:16S rRNAけいとうかいせき}}