16SリボソームRNA

16S リボソームRNA(16S rRNA)とは、シャイン・ダルガノ配列(Shine-Dalgarno sequence)に結合する原核生物リボソームの30S小サブユニットのコンポーネントである。このRNAをコードする遺伝子は16S rRNA遺伝子と呼ばれる。 16S rRNA遺伝子は、リボソームという生物の本質に関わる機能を持つRNAであるため、配列の保存性が高く、細菌や古細菌といった原核生物の間で高度に保存されている。そして、機能変化に伴う遺伝子の変異がこれからも起きる可能性が極めて低い。すなわち、遺伝子配列の進化速度が遅いことから、信頼できる分子時計として利用できる。また、遺伝子の長さが適当に長く（16S rRNAの場合、1600塩基対程度）、系統解析を行う上で十分な情報量を持つ。さらに、比較的変異しやすい部位も存在し、近縁な種でも比較が可能である。これらの特徴から、特に微生物系統学の分野において、この遺伝子配列は系統進化解析によく利用されている。カール・ウーズとジョージ・E・フォックス（英語版）が1977年に、系統学に16S rRNAを導入した。真核生物の場合に対応するものは18S rRNAなので、まとめてリボソーム小サブユニットRNA (Small Sub Unit rRNA、SSU rRNA) 系統解析と呼ばれることもある。 16S rRNAは23S rRNAと相互作用する。このRNA複合体は構造上、リボソームタンパク質の位置を決める足場として機能する役割を持ち、2つのリボソームサブユニット（50Sおよび30S）の結合を支援する。3 '末端には、mRNAのAUG開始コドンの上流に結合するShine-Dalgarno配列の相補鎖が含まれている。16S RNAの3 '末端は、タンパク質合成の開始に関与するS1およびS21タンパク質に結合する。 1492残基および1493残基でアデニン(A)が並んでいる箇所をAサイトと呼ぶ、このAサイトのアデニンが持つN1原子とmRNA骨格の2つのOH基の間に水素結合を形成し、正確なコドン-アンチコドンのペアリングを安定化している。 細菌の16S rRNA遺伝子には、リボソーム小サブユニットの二次構造に関与する、9つの超可変領域（V1-V9）が含まれており、これらの長さは約30-100塩基対である。保存の程度は超可変領域間で大きく異なり、より保存された領域は門や綱といったより高レベルの分類法に利用でき、一方で保存度の低い領域は属や種といったより低レベルの分類に利用される。16S rRNA配列全体をシーケンスすることで全超可変領域の比較が可能になるが、16S rRNAは約1,500塩基の長さを持つため、多様な細菌群集を満遍なくシーケンスするには費用がかかってしまう。そのため、細菌叢解析のような研究では通常、Illumina社製のゲノムシーケンス技術を利用しており、454パイロシーケンスやサンガーシーケンスよりもそれぞれ約50倍、12,000倍ほど安価にシーケンスすることができる。しかしながら、Illumina社製シーケンサーでは75〜250塩基（Illumina MiSeqでは最大300塩基）のリード長しか得られないため、細菌叢サンプルから16S rRNA遺伝子配列を完璧に組み立てることはできない。一方で、超可変領域はその短さのため、Illuminaシーケンサを1回実行するだけで配列解析を行えるため、この超可変領域は菌叢解析における理想的なターゲットになっている。 16S rRNA超可変領域は、細菌系統間で大きく配列が異なる場合があるが、全体としては16S rRNA遺伝子は真核生物(18SリボソームRNA)よりも良く均一性を維持しているため、アライメントが比較的容易である。さらに、16S rRNA遺伝子には超可変領域間の高度に保存された配列が含まれているため、異なる分類群にわたって同じ超可変領域を確実にPCR増幅できるユニバーサルプライマーの設計が可能である。すべての細菌系統をドメインから種に渡って正確に分類できる超可変領域は存在しないが、特定の分類レベルをほぼ確実に予測できるものもは知られている。多くの菌叢解析研究では、完全な16S rRNA遺伝子と同程度の正確性で門レベルの系統解析を行うことができるV4超可変領域を選択することが多い。保存度の低い領域は、高次の系統分類には不向きであるが、例えば特定の病原体を検出するような用途でよく利用される。2007年にChakravortyらが発表した研究では、どの超可変領域が疾患特異的かつ広範なアッセイに利用できるかを調べ、さまざまな病原体のV1-V8領域を報告している。また他の研究では、病原体の属の特定にはV3領域を利用することが最適であり、炭疽菌を含むテストされたすべてのCDC監視病原体においてはV6領域が種の区別に最も高い正確性を示した、と報告されている。 16S rRNA超可変領域をベースとした配列解析は、細菌系統の分類学的研究にとって有用であるが、ごく近縁の種同士を区別することは困難な場合がある。例えば腸内細菌科、クロストリジウム科、およびペプトストレプトコッカス科では、種間で16S rRNA遺伝子全体の最大99%の配列類似性をもつことが知られている。この場合、種間差異はV4配列中のほんの数塩基にしか出現しないため、特に低レベルの分類において、参照データベースに基づく手法では確実に分類することが困難である。また、利用する超可変領域の数を絞るほど、近縁な分類群の違いを観察できなくなり、サンプル全体の多様性の過小評価に繋がりうる。さらに、細菌のゲノムは、多様なV1、V2、V6領域の配列を持つ複数の16S rRNA遺伝子をマルチコピーで保持する場合がある。これらの理由から、16S rRNAの超可変領域に基づく解析は、細菌種を分類する完璧な方法とまでは言えない。しかしながらこのような欠点がありつつも、現実的には細菌群集研究に利用できる最も有用なツールの1つとして今日利用されている。 16S rRNA配列を解析する際は、ユニバーサルプライマーを用いてPCRによる増幅を行い、得られた増幅産物をシーケンスする方法が一般的である。シークエンシング反応を行わなくても群集構造の解析が可能なDGGE (Denaturing Gradient Gel Electrophoresis) 法や、顕微鏡で直接観察できるFISH法などの広い応用範囲も知られている。かつては制限酵素を用いたRFLPなどが使用されていた。 最も一般的なプライマーペアは、Weisburgらによって考案された、27F-1492Rと呼ばれているセットである。一部のアプリケーションでは、より短いアンプリコンが必要になる場合があり、たとえばチタンケミストリーを使用した454シーケンスでは、V1からV3をカバーするプライマーペア27F-534Rがよく選択される。また、27Fではなく8Fが使用される場合も多い。この2つのプライマーはほぼ同じであるが、27FはCではなくM（AとCの縮退塩基表記）を持っている。 16S rRNA遺伝子配列には、高度に保存されたプライマー結合部位に加えて、複数の超可変領域が含まれており、この領域の塩基配列を利用することで、細菌の系統的な同定を行うことができる。現在、16S rRNA遺伝子配列シーケンシングは、表現型をベースとした細菌同定法に代わる迅速で安価な代替手法として、医学の分野で広く普及している。また、細菌の識別のみならず、完全に新種な系統の発見や系統関係の再分類にも利用されている。未培養系統の新種記載においても利用される。次世代シーケンシング技術を活用することで、数千の16S rRNA配列を数時間程度で解析することが可能になっており、たとえば腸内細菌叢のメタゲノム研究などに利用されている。 細菌が持つ16S rRNA遺伝子配列は一つとは限らず、複数の16S rRNA遺伝子がゲノム中にマルチコピーで含まれることが多い。また例外として、一部の好熱性古細菌 (例えばThermoprotealesなど)には、16S rRNA遺伝子中にイントロンが含まれており、ユニバーサルプライマーのアニーリングに影響を与える可能性がある。また、サンプル中の真核生物に由来するミトコンドリアや葉緑体が持つ16S rRNAもPCRで増幅されることがある。また、ユニバーサルプライマーを用いた群集構造解析は、環境中に存在している16S rRNAをすべて増幅してしまうために、生存個体のみならず死亡して溶菌したようなRNAの残骸をも増幅しうる。 進化が垂直伝達によって駆動されるという仮定の下では、16S rRNA遺伝子は種特異的であるみなすことができ、原核生物間の系統関係を推測する確実な遺伝的マーカーであると長年考えられてきた。しかしながら、研究が進むに連れ、これらの遺伝子においても遺伝子の水平伝播が発生していることが分かってきた。このような遺伝子の転移性は、特別な大腸菌の遺伝子システムを用いた実験的によって確認されている。すなわち、大腸菌が本来持つ16S rRNA遺伝子を欠失させ、大腸菌とは綱あるいは門レベルで系統が異なる生物種由来の外来16S rRNA遺伝子を導入したところ、変異株として増殖することが示された。このような門レベルで異なる16S rRNA遺伝子の機能的互換性は、サーマスサーモフィルスでも確認されている。さらに、T. thermophilusでは、遺伝子全長の置換と部分的な置換の両方が観察された。部分的な置換は、宿主（T. thermophilus）と外来細菌の16S rRNA遺伝子間でさまざまなキメラが生成されることによる。このように、16S rRNA遺伝子は、垂直遺伝と水平遺伝子伝播を含む複数のメカニズムを通じて進化している可能性があり、特に後者については今まで考えられてきたよりもはるかに高い頻度で発生している可能性がある。 16S rRNA遺伝子はほぼ全ての微生物に存在し、適当に配列変化が起きるため、微生物の系統分類と同定に利用されてきた。ほとんどの細菌および古細菌のタイプ株が持つ16S rRNA遺伝子の配列情報は、NCBIなどの公共データベースから入手できる。ただし、これらのデータベースに格納された配列は、品質が検証されていないことがよくある。そのため、16S rRNA配列のみを収集する2次データベースが広く使用されている。使用されるケースが多い有名なデータベースは以下のとおりである。 EzBioCloudデータベースは、以前はEzTaxonと呼ばれていた。2018年9月の時点で15,290の有効な公開名を含む62,988の細菌と古細菌の系統を含んでおり、完全な階層分類システムで構成されている。最尤推定やOrthoANIなどに基づいた系統関係に基づいて、すべての種/亜種が少なくとも1つの16S rRNA遺伝子配列によって表されている。EzBioCloudデータベースは体系的に管理されており、定期的に更新されている。新しい候補種が登録されることもある。さらにWebサイト上では、ANIの計算やContEst16S、QIIMEおよびMothurパイプライン用の16S rRNA DBといったバイオインフォマティクスツールを提供している。 Ribosomal Database Project（RDP）は、関連するツール群やサービスと共にリボソームデータを提供するキュレーションデータベースである。系統的に纏められたリボソームRNA（rRNA）配列のアライメントや系統樹、rRNA二次構造図、およびアライメントの分析や表示をするさまざまなソフトウェアパッケージを提供している。 SILVAは、小サブユニット（16S/18S、SSU）と大サブユニット（23S/28S、LSU）のリボソームRNA配列を包括的に纏めたデータベースである。キュレーションを経たデータセットを定期的に更新している。 Greengenesは、品質管理された包括的な16Sリファレンスデータベースである。de novo系統に基づいて分類されており、標準的な操作上の分類単位を提供する。現在は積極的に維持されておらず、最後の更新は2013年である。 従来原核生物の分類は細胞の形態、分離の条件、染色法などで行っていたが、こうした表現型の形質では系統樹上の上下関係を説明するには至らなかった。しかし1970年代、シトクロム、フェレドキシン、5S rRNAなどの塩基配列を基にした系統分類が分子生物学の発展とともに徐々に活発化してきた。 遺伝子の一次構造に基づく系統分類は原核生物に対して特に有効であった。カール・ウーズらはリボソーム小サブユニットを構成するRNA、つまり16S rRNAの塩基配列を用いて原核生物の系統分類を行い、原核生物が真正細菌と古細菌という2つのドメインからなることを示した（1977年当時はオリゴヌクレオチドカタログ法を用いた）。 現在、16S rRNAを用いた系統解析は、系統樹の作成のみならず、任意の環境中における細菌・古細菌の群集構造の観測に役立っている。この方法を用いると、分離・培養が困難な難培養性の菌種を含めて網羅的に群集構造を明らかにできる他、新規の菌の存在を配列解析から明らかにする事ができる（1996年のBarnsによる；メタゲノミクス参照）。 ^ “Structure of functionally activated small ribosomal subunit at 3.3 angstroms resolution”. Cell 102 (5): 615–23. (September 2000). doi:10.1016/S0092-8674(00)00084-2. PMID 11007480.  ^ “16S ribosomal DNA amplification for phylogenetic study”. Journal of Bacteriology 173 (2): 697–703. (January 1991). doi:10.1128/jb.173.2.697-703.1991. PMC 207061. PMID 1987160. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC207061/.  ^ “Intragenomic heterogeneity between multiple 16S ribosomal RNA operons in sequenced bacterial genomes”. FEMS Microbiology Letters 228 (1): 45–9. (November 2003). doi:10.1016/S0378-1097(03)00717-1. PMID 14612235.  ^ “Comparative RNA function analysis reveals high functional similarity between distantly related bacterial 16 S rRNAs” (英語). Scientific Reports 7 (1): 9993. (August 2017). Bibcode: 2017NatSR...7.9993T. doi:10.1038/s41598-017-10214-3. PMC 5577257. PMID 28855596. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5577257/.  ^ “Phylogenetic structure of the prokaryotic domain: the primary kingdoms”. Proceedings of the National Academy of Sciences of the United States of America 74 (11): 5088–90. (November 1977). Bibcode: 1977PNAS...74.5088W. doi:10.1073/pnas.74.11.5088. PMC 432104. PMID 270744. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC432104/.  ^ “Towards a natural system of organisms: proposal for the domains Archaea, Bacteria, and Eucarya”. Proceedings of the National Academy of Sciences of the United States of America 87 (12): 4576–9. (June 1990). Bibcode: 1990PNAS...87.4576W. doi:10.1073/pnas.87.12.4576. PMC 54159. PMID 2112744. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC54159/.  ^ “Phylogenetic structure of the prokaryotic domain: the primary kingdoms”. Proceedings of the National Academy of Sciences of the United States of America 74 (11): 5088–90. (November 1977). Bibcode: 1977PNAS...74.5088W. doi:10.1073/pnas.74.11.5088. PMC 432104. PMID 270744. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC432104/.  ^ Czernilofsky, A. P.; Kurland, C. G.; Stöffler, G. (1975). “30S Ribosomal proteins associated with the 3′-terminus of 16S RNA”. FEBS Letters 58 (1): 281–284. doi:10.1016/0014-5793(75)80279-1. ISSN 0014-5793. PMID 1225593.  ^ a b “A quantitative map of nucleotide substitution rates in bacterial rRNA”. Nucleic Acids Research 24 (17): 3381–91. (September 1996). doi:10.1093/nar/24.17.3381. PMC 146102. PMID 8811093. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC146102/.  ^ “On the evolutionary descent of organisms and organelles: a global phylogeny based on a highly conserved structural core in small subunit ribosomal RNA”. Nucleic Acids Research 12 (14): 5837–52. (July 1984). doi:10.1093/nar/12.14.5837. PMC 320035. PMID 6462918. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC320035/.  ^ a b “Sensitivity and correlation of hypervariable regions in 16S rRNA genes in phylogenetic analysis”. BMC Bioinformatics 17 (1): 135. (March 2016). doi:10.1186/s12859-016-0992-y. PMC 4802574. PMID 27000765. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4802574/.  ^ “Generation of multimillion-sequence 16S rRNA gene libraries from complex microbial communities by assembling paired-end illumina reads”. Applied and Environmental Microbiology 77 (11): 3846–52. (June 2011). doi:10.1128/AEM.02772-10. PMC 3127616. PMID 21460107. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3127616/.  ^ a b “A method for high precision sequencing of near full-length 16S rRNA genes on an Illumina MiSeq”. PeerJ 4: e2492. (2016-09-20). doi:10.7717/peerj.2492. PMC 5036073. PMID 27688981. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5036073/.  ^ “The variability of the 16S rRNA gene in bacterial genomes and its consequences for bacterial community analyses”. PLOS ONE 8 (2): e57923. (2013-02-27). Bibcode: 2013PLoSO...857923V. doi:10.1371/journal.pone.0057923. PMC 3583900. PMID 23460914. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3583900/.  ^ a b “Sensitivity and correlation of hypervariable regions in 16S rRNA genes in phylogenetic analysis”. BMC Bioinformatics 17 (1): 135. (March 2016). doi:10.1186/s12859-016-0992-y. PMC 4802574. PMID 27000765. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4802574/.  ^ a b “A detailed analysis of 16S ribosomal RNA gene segments for the diagnosis of pathogenic bacteria”. Journal of Microbiological Methods 69 (2): 330–9. (May 2007). doi:10.1016/j.mimet.2007.02.005. PMC 2562909. PMID 17391789. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2562909/.  ^ a b “The variability of the 16S rRNA gene in bacterial genomes and its consequences for bacterial community analyses”. PLOS ONE 8 (2): e57923. (2013-02-27). Bibcode: 2013PLoSO...857923V. doi:10.1371/journal.pone.0057923. PMC 3583900. PMID 23460914. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3583900/.  ^ a b c “Characterization of the Gut Microbiome Using 16S or Shotgun Metagenomics”. Frontiers in Microbiology 7: 459. (2016-01-01). doi:10.3389/fmicb.2016.00459. PMC 4837688. PMID 27148170. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4837688/.  ^ “Intragenomic heterogeneity between multiple 16S ribosomal RNA operons in sequenced bacterial genomes”. FEMS Microbiology Letters 228 (1): 45–9. (November 2003). doi:10.1016/S0378-1097(03)00717-1. PMID 14612235.  ^ “16S ribosomal DNA amplification for phylogenetic study”. Journal of Bacteriology 173 (2): 697–703. (January 1991). doi:10.1128/jb.173.2.697-703.1991. PMC 207061. PMID 1987160. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC207061/.  ^ http://www.hmpdacc.org/tools_protocols.php#sequencing Archived 2010-10-30 at the Wayback Machine. ^ a b “Phylogenetic analysis of Aquaspirillum magnetotacticum using polymerase chain reaction-amplified 16S rRNA-specific DNA”. International Journal of Systematic Bacteriology 41 (2): 324–5. (April 1991). doi:10.1099/00207713-41-2-324. PMID 1854644.  ^ a b James, Greg (15 May 2018). “Universal Bacterial Identification by PCR and DNA Sequencing of 16S rRNA Gene”. PCR for Clinical Microbiology. Springer, Dordrecht. pp. 209–214. doi:10.1007/978-90-481-9039-3_28. ISBN 978-90-481-9038-6  ^ a b “Diversity of uncultured microorganisms associated with the seagrass Halophila stipulacea estimated by restriction fragment length polymorphism analysis of PCR-amplified 16S rRNA genes”. Applied and Environmental Microbiology 62 (3): 766–71. (March 1996). PMC 167844. PMID 8975607. http://aem.asm.org/cgi/reprint/62/3/766.pdf.  ^ “Microbial diversity in water and sediment of Lake Chaka, an athalassohaline lake in northwestern China”. Applied and Environmental Microbiology 72 (6): 3832–45. (June 2006). doi:10.1128/AEM.02869-05. PMC 1489620. PMID 16751487. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1489620/.  ^ “Identification of species by multiplex analysis of variable-length sequences”. Nucleic Acids Research 38 (22): e203. (December 2010). doi:10.1093/nar/gkq865. PMC 3001097. PMID 20923781. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3001097/.  ^ “Ribosomal DNA sequencing as a tool for identification of bacterial pathogens”. Current Opinion in Microbiology 2 (3): 299–305. (June 1999). doi:10.1016/S1369-5274(99)80052-6. PMID 10383862.  ^ “Impact of 16S rRNA gene sequence analysis for identification of bacteria on clinical microbiology and infectious diseases”. Clinical Microbiology Reviews 17 (4): 840–62, table of contents. (October 2004). doi:10.1128/CMR.17.4.840-862.2004. PMC 523561. PMID 15489351. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC523561/.  ^ “Reverse transcription of 16S rRNA to monitor ribosome-synthesizing bacterial populations in the environment”. Applied and Environmental Microbiology 75 (13): 4589–98. (July 2009). doi:10.1128/AEM.02970-08. PMC 2704851. PMID 19395563. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2704851/.  ^ “16S ribosomal DNA amplification for phylogenetic study”. Journal of Bacteriology 173 (2): 697–703. (January 1991). doi:10.1128/jb.173.2.697-703.1991. PMC 207061. PMID 1987160. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC207061/.  ^ “Burkholderia thailandensis sp. nov., a Burkholderia pseudomallei-like species”. International Journal of Systematic Bacteriology 48 Pt 1 (1): 317–20. (January 1998). doi:10.1099/00207713-48-1-317. PMID 9542103.  ^ Phylogenetic identification of uncultured pathogens using ribosomal RNA sequences. Methods in Enzymology. 235. (1994). pp. 205–222. doi:10.1016/0076-6879(94)35142-2. ISBN 978-0-12-182136-4. PMID 7520119. https://archive.org/details/bacterialpathoge0000unse_f8b6/page/205  ^ “Phylogenetic analysis of the bacterial communities in marine sediments”. Applied and Environmental Microbiology 62 (11): 4049–59. (November 1996). PMC 168226. PMID 8899989. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC168226/.  ^ “Next-generation sequencing of 16S ribosomal RNA gene amplicons”. Journal of Visualized Experiments (90). (August 2014). doi:10.3791/51709. PMC 4828026. PMID 25226019. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4828026/.  ^ “Use of 16S rRNA and rpoB genes as molecular markers for microbial ecology studies”. Applied and Environmental Microbiology 73 (1): 278–88. (January 2007). doi:10.1128/AEM.01177-06. PMC 1797146. PMID 17071787. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1797146/.  ^ “The distribution, diversity, and importance of 16S rRNA gene introns in the order Thermoproteales”. Biology Direct 10 (35): 35. (July 2015). doi:10.1186/s13062-015-0065-6. PMC 4496867. PMID 26156036. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4496867/.  ^ “Mutational robustness of 16S ribosomal RNA, shown by experimental horizontal gene transfer in Escherichia coli”. Proceedings of the National Academy of Sciences of the United States of America 109 (47): 19220–5. (November 2012). doi:10.1073/pnas.1213609109. PMC 3511107. PMID 23112186. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3511107/.  ^ “Comparative RNA function analysis reveals high functional similarity between distantly related bacterial 16 S rRNAs”. Scientific Reports 7 (1): 9993. (August 2017). doi:10.1038/s41598-017-10214-3. PMC 5577257. PMID 28855596. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5577257/.  ^ “Occurrence of randomly recombined functional 16S rRNA genes in Thermus thermophilus suggests genetic interoperability and promiscuity of bacterial 16S rRNAs”. Scientific Reports 9 (1): 11233. (August 2019). doi:10.1038/s41598-019-47807-z. PMC 6677816. PMID 31375780. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6677816/.  ^ “Uniting the classification of cultured and uncultured bacteria and archaea using 16S rRNA gene sequences”. Nature Reviews. Microbiology 12 (9): 635–45. (September 2014). doi:10.1038/nrmicro3330. PMID 25118885. https://www.researchgate.net/publication/264798721.  ^ 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 ^ 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. ^ 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. ^ “Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB”. Applied and Environmental Microbiology 72 (7): 5069–72. (July 2006). doi:10.1128/aem.03006-05. PMC 1489311. PMID 16820507. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1489311/.  ^ “An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea”. The ISME Journal 6 (3): 610–8. (March 2012). doi:10.1038/ismej.2011.139. PMC 3280142. PMID 22134646. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3280142/.