Ribosomal RNA (rRNA): Structure, Function, and Taxonomic Profiling in Metagenomics

1. Introduction

Ribosomal RNA (rRNA) constitutes the most abundant class of RNA in living cells and serves as both the structural scaffold and the catalytic core of the ribosome, the macromolecular complex responsible for protein synthesis [20, 27]. Across all domains of life, rRNA molecules are organized into two ribosomal subunits. In prokaryotes, the small subunit (SSU) contains the 16S rRNA, while the large subunit (LSU) contains the 23S and 5S rRNAs [29]. In eukaryotes, the SSU houses the 18S rRNA, and the LSU contains the 28S, 5.8S, and 5S rRNAs [1, 22]. The genes encoding rRNA (rDNA) are organized in operons that are repeated in multiple copies within the genome, a feature that facilitates high levels of ribosome production and provides a rich target for molecular detection [1, 2].

The evolutionarily conserved nature of rRNA, combined with regions of sequence hypervariability, has made rRNA genes the marker of choice for phylogenetic reconstruction and taxonomic identification of microorganisms [3, 4, 5]. With the advent of metagenomics and environmental DNA (eDNA) analysis, rRNA amplicon sequencing has become a cornerstone of culture-independent microbial community profiling, enabling the characterization of complex microbiomes from environmental, veterinary, and agricultural samples [6, 7]. This article reviews the molecular structure and functional roles of rRNA, details the mechanisms of ribosome biogenesis and translational regulation, and provides a comprehensive overview of bioinformatic strategies for rRNA-based taxonomic profiling in metagenomic workflows.

2. Structural Organization of Ribosomal RNA

2.1. Prokaryotic rRNA Operons

In bacteria, the rRNA genes are typically organized in a single operon with the order 16S-23S-5S, transcribed by RNA polymerase as a single polycistronic precursor [8, 29]. The 16S rRNA is approximately 1,542 nucleotides in length in Escherichia coli and folds into a complex secondary structure comprising three major domains: the 5 prime domain, the central domain, and the 3 prime domain [29, 30]. The 23S rRNA is approximately 2,904 nucleotides long and forms six structural domains, many of which participate in the peptidyl transferase center (PTC) [30]. The 5S rRNA is a much shorter molecule of about 120 nucleotides that contributes to LSU stability [30].

The 16S rRNA contains conserved sequence motifs that are universal across bacteria, interspersed with hypervariable regions (V1 through V9) that provide phylogenetic signal at different taxonomic ranks [6, 4]. The 3 prime terminal region of the 16S rRNA contains the anti-Shine-Dalgarno (aSD) sequence (CCUCCU), which base-pairs with the Shine-Dalgarno motif upstream of the start codon on bacterial mRNAs to position the ribosome at the correct initiation site [9].

2.2. Eukaryotic rRNA Operons

Eukaryotic rRNA genes are organized in tandem repeats on chromosomes, with copy numbers ranging from dozens to thousands per genome [1, 2]. The primary transcript (pre-rRNA) is transcribed by RNA polymerase I and contains the 18S, 5.8S, and 28S rRNAs separated by internal transcribed spacers (ITS1 and ITS2) and flanked by external transcribed spacers (5 prime ETS and 3 prime ETS) [22, 28]. The 5S rRNA is transcribed separately by RNA polymerase III [1]. The 18S rRNA (approximately 1,800 nt in mammals) is the eukaryotic SSU counterpart, the 28S rRNA (approximately 4,700 nt) is the LSU component, and the 5.8S rRNA (approximately 160 nt) forms a base-paired complex with the 28S rRNA [1, 22]. The ITS regions, particularly ITS1, evolve rapidly and are highly variable, making them a preferred barcode marker for fungal identification [35].

2.3. Secondary and Tertiary Structure

All rRNA molecules adopt highly conserved secondary structures stabilized by Watson-Crick base pairing and numerous non-canonical interactions [30]. The secondary structures of 16S and 23S rRNA were first inferred through comparative sequence analysis, a method that remains the gold standard for RNA structure prediction [30]. The Comparative RNA Web (CRW) site provides curated secondary structure models for 5S, 16S, and 23S rRNAs based on alignment of thousands of sequences [30].

At the tertiary level, rRNA folds into a compact, globular architecture within the ribosomal subunits. The functional centers of the ribosome, including the decoding site on the SSU and the PTC on the LSU, are composed entirely of rRNA, establishing the ribosome as a ribozyme [29, 30]. Cryo-electron microscopy (cryo-EM) studies have provided atomic-resolution models of the rRNA tertiary structure in various conformational states, revealing the precise positioning of rRNA helices, loops, and modified nucleotides [10]. These structural insights are fundamental to understanding ribosomal mechanics and are leveraged by computational methods such as RNA structure prediction algorithms.

3. Functional Roles of rRNA in Translation

3.1. Peptidyl Transferase Center

The PTC is located in domain V of the 23S rRNA and catalyzes the formation of peptide bonds between aminoacyl-tRNA in the A site and peptidyl-tRNA in the P site [30]. No protein components participate directly in catalysis; instead, the 23S rRNA positions the substrates and stabilizes the transition state through hydrogen bonding and electrostatic interactions [30]. The universality of this catalytic mechanism across all domains of life underscores the ancient origin of rRNA as the primordial enzyme [25].

3.2. Decoding Center and mRNA Interaction

The decoding center of the SSU is formed by conserved residues of the 16S rRNA (or 18S rRNA in eukaryotes) that monitor the geometry of codon-anticodon base pairing [29, 30]. The 3 prime end of the 16S rRNA, containing the aSD sequence, directly contacts the mRNA, stabilizing the initiation complex through base pairing with the Shine-Dalgarno sequence [9]. During elongation, the decoding center ensures that only the correct aminoacyl-tRNA is accepted, a fidelity mechanism that depends on rRNA conformational changes [10].

3.3. Programmed Ribosomal Frameshifting

Certain RNA viruses, including coronaviruses, utilize programmed ribosomal frameshifting (PRF) to regulate expression of overlapping open reading frames [10]. In PRF, a pseudoknot or other RNA secondary structure in the mRNA stalls the ribosome and induces a minus-1 shift in the reading frame. Cryo-EM structures of the mammalian ribosome primed for frameshifting on viral RNA reveal that the pseudoknot lodges at the entry to the mRNA channel, generating tension that promotes back slippage [10]. This mechanism is crucial for the expression of the viral RNA-dependent RNA polymerase and is a target for antiviral intervention [10]. Computational approaches for in silico modeling of ribosome stalling and programmed ribosomal frameshifting in RNA viruses provide a powerful means to predict and characterize these regulatory events.

3.4. Ribosome Heterogeneity and Regulatory Roles

Emerging evidence indicates that ribosomes are not uniform molecular machines. Variant rRNA alleles, encoded by different rDNA copies, are conserved across individuals and exhibit tissue-specific expression in mammals [2]. These variant rRNAs are incorporated into actively translating ribosomes and can influence the translational efficiency of specific mRNAs [2]. Additionally, rRNA undergoes extensive post-transcriptional modification, including 2 prime-O-methylation, pseudouridylation, and base methylation, which collectively contribute to ribosome heterogeneity and translational regulation [21, 32].

4. Ribosome Biogenesis and rRNA Modification

4.1. Pre-rRNA Transcription and Processing

Ribosome biogenesis begins with the transcription of pre-rRNA by RNA polymerase I (for the 18S, 5.8S, and 28S rRNAs) and RNA polymerase III (for the 5S rRNA) [1, 22]. The nascent pre-rRNA associates co-transcriptionally with assembly factors, including the UtpA and UtpB complexes, which chaperone the RNA and promote correct folding [28]. Processing of the pre-rRNA involves a series of endonucleolytic and exonucleolytic cleavages that remove the external and internal transcribed spacers, releasing the mature rRNAs [22]. The RNA helicase DHX37 is required for release of the U3 small nucleolar ribonucleoprotein (snoRNP) from pre-ribosomal particles, an essential step that enables formation of the central pseudoknot in the SSU [26].

4.2. Nucleotide Modifications

Mature rRNAs contain a dense array of chemically modified nucleotides. In eukaryotes, over 200 sites are modified, predominantly 2 prime-O-methylation and pseudouridylation, guided by box C/D and box H/ACA snoRNPs, respectively [21]. Base methylations, such as N6-methyladenosine (m6A) and 5-methylcytosine (m5C), are also present [11, 12, 32]. In eukaryotic 28S rRNA, m6A at position 4220 is catalyzed by the methyltransferase ZCCHC4, and loss of this modification reduces global translation and cell proliferation [12]. Similarly, the m5C methyltransferase NSUN5 (the functional orthologue of yeast Rcm1) modifies a conserved cytosine in the 25S/28S rRNA, and its ablation impairs translation fidelity and affects organismal growth and lifespan in metazoans [11, 34].

Ribosomal RNA methylation also has clinical implications in antimicrobial resistance. Acquired 16S rRNA methyltransferases in pathogenic bacteria confer high-level resistance to aminoglycoside antibiotics by methylating specific residues in the decoding center (e.g., m7G1405), preventing drug binding [24]. This mechanism is an emerging concern in veterinary medicine.

4.3. Small RNAs Derived from rRNA

Ribosomal RNA fragments (rRFs) are a class of small non-coding RNAs generated by endonucleolytic cleavage of mature rRNA [13]. Initially dismissed as degradation products but now recognized as potentially functional, rRFs can associate with Argonaute proteins and may regulate gene expression, similar to microRNAs [13]. The extent to which rRFs contribute to translational control or other cellular processes in veterinary species remains an active area of investigation.

5. Ribosomal RNA in Metagenomic Taxonomic Profiling

5.1. Rationale for rRNA as a Phylogenetic Marker

The rRNA genes possess several properties that make them ideal for taxonomic profiling. They are universally distributed across all cellular life, contain both conserved and variable regions, and are present in multiple copies per genome, which facilitates their amplification from samples with low biomass [3, 4, 5]. The 16S rRNA gene is the standard marker for bacteria and archaea [6], while the 18S rRNA gene is used for eukaryotic microbial profiling [14]. For fungi, the internal transcribed spacer (ITS) region of the rRNA operon, which evolves more rapidly than the rRNA genes themselves, provides greater discrimination at the species level [35].

5.2. Primer Selection and Amplification Bias

The accuracy of rRNA amplicon surveys depends critically on the choice of PCR primers, which must amplify a broad range of taxa while minimizing bias. Comprehensive in silico evaluations of 16S rRNA primers have assessed the coverage and phylum spectrum of 175 individual primers and 512 primer pairs against the SILVA reference database [6]. The primer pair 341F-785R (targeting the V3-V4 hypervariable region) was identified as having excellent coverage for both Bacteria and Archaea and is widely adopted for short-read amplicon sequencing [6]. For longer amplicons, primers targeting near-full-length 16S rRNA are used with long-read sequencing platforms [6].

Despite careful primer design, PCR amplification inevitably introduces bias. Differences in template secondary structure, GC content, and copy number can distort the observed community composition relative to the true composition [6]. These biases are documented and should be considered when interpreting amplicon data.

5.3. Amplicon Sequencing versus Metagenomic Shotgun

Two primary sequencing strategies are used for rRNA-based community profiling. Amplicon sequencing involves PCR amplification of a specific rRNA gene region followed by high-throughput sequencing of the resulting amplicons [6]. This approach is cost-effective and achieves deep coverage of the target gene, but it captures only taxonomic information and does not provide functional gene content. Metagenomic shotgun sequencing, in contrast, sequences all DNA in a sample and allows extraction of rRNA reads in silico, either by alignment to rRNA databases or through tools such as RNAmmer [8]. Shotgun sequencing provides both taxonomic and functional information but is more expensive and computationally demanding. A hybrid approach, using amplicon data for high-resolution taxonomy and shotgun data for functional annotation, is often employed.

5.4. Bioinformatic Processing Pipelines

The bioinformatic processing of rRNA amplicon data follows a established workflow. Raw reads are quality-filtered, denoised to correct sequencing errors, and grouped into operational taxonomic units (OTUs) or amplicon sequence variants (ASVs). The traditional 97% sequence identity threshold for OTU clustering corresponds approximately to species-level distinction in many bacterial lineages, although this threshold has been critically re-evaluated [15]. ASV-based methods, such as DADA2, provide single-nucleotide resolution and are increasingly preferred.

Taxonomic assignment is performed by comparing representative sequences to reference databases. The metagenomics taxonomic classification Kraken2 and functional annotation pipelines article describes k-mer-based approaches for rapid classification. For rRNA-specific classifiers, the RDP classifier uses a naive Bayesian approach trained on 16S rRNA sequences. The SILVA database provides curated, aligned rRNA sequences with a consistent taxonomy that includes clades of uncultivated environmental sequences [16, 14, 7].

5.5. Reference Databases

The SILVA ribosomal RNA gene database is the most comprehensive and widely used resource for rRNA-based profiling [16, 14, 7]. SILVA provides quality-controlled, aligned SSU and LSU rRNA sequences from Bacteria, Archaea, and Eukaryota. The project originated from the ARB software workbench and has been maintained and expanded continuously [7]. The SILVA Incremental Aligner (SINA) enables accurate, high-throughput alignment of new sequences against the reference alignment [17]. The EzTaxon server provides a curated database of 16S rRNA sequences specifically for the identification of prokaryotes [18].

The Comparative RNA Web (CRW) site offers comparative sequence and structure models for rRNA and other RNAs, supporting both phylogenetic analysis and RNA structure research [30]. For fungal identification, the ITS region has been adopted as the universal DNA barcode, with dedicated reference databases such as UNITE [35].

6. Computational Methods for rRNA Analysis

6.1. Multiple Sequence Alignment

Accurate multiple sequence alignment (MSA) is essential for phylogenetic inference and taxonomic classification of rRNAs. SINA, the aligner developed for SILVA, combines k-mer searching with partial order alignment (POA) to achieve high accuracy at throughputs sufficient for large datasets [17]. In benchmarks, SINA reproduced reference MSAs from the BRAliBase III set with accuracies exceeding 96% and outperformed other high-throughput aligners such as PyNAST and mothur [17].

6.2. Gene Annotation in Genomes

For complete genome sequences, rRNA genes are annotated using predictors such as RNAmmer, which employs hidden Markov models trained on data from the 5S rRNA database and the European ribosomal RNA database project [8]. RNAmmer detects the major rRNA species (5S, 16S, and 23S for bacteria; 5S, 5.8S, 18S, and 28S for eukaryotes) with high sensitivity and specificity, enabling rapid annotation of even incomplete genomes [8].

6.3. Co-extraction of DNA and RNA for Combined Analysis

To distinguish active from dormant microbial community members, co-extraction of DNA and RNA from environmental samples allows parallel analysis of rDNA (taxonomic potential) and rRNA (taxonomic activity) [31]. Griffiths and colleagues described a rapid method for co-extraction of DNA and RNA from soils and sediments, yielding nucleic acids suitable for PCR and reverse transcription PCR amplification of rRNA genes [31]. This dual approach provides insight into the metabolically active fraction of a microbiome.

flowchart TD
    A[Sample Collection], > B[Nucleic Acid Extraction]
    B, > C{Sequencing Strategy}
    
    C, > D[Amplicon Sequencing\n(16S/18S/ITS PCR)]
    C, > E[Metagenomic Shotgun\nSequencing]
    
    D, > F[Quality Filtering\nand Denoising]
    E, > G[Quality Filtering\nand Host Read Depletion]
    
    F, > H[OTU/ASV Clustering\n(97% identity)]
    G, > I[Extraction of rRNA\nReads in Silico]
    
    H, > J[Taxonomic Assignment\nvia Reference Database\n(SILVA, RDP, EzTaxon)]
    I, > J
    
    J, > K[Community Composition\nProfiling]
    
    K, > L[Downstream Analysis\nAlpha/Beta Diversity\nDifferential Abundance]
    
    subgraph Databases
        M[SILVA Database]
        N[CW Site]
        O[EzTaxon]
        P[UNITE]
    end
    
    J -.-> M
    J -.-> N
    J -.-> O
    J -.-> P

7. Applications in Veterinary Microbiology

7.1. Microbiome Analysis in Production Animals

Ribosomal RNA profiling is extensively used to characterize the gut, respiratory, and urogenital microbiomes of livestock species, including cattle, swine, poultry, and horses. Changes in the composition of the gut microbiota have been associated with production traits, feed efficiency, and disease susceptibility. In poultry, 16S rRNA amplicon sequencing has been used to monitor the effects of feed additives on cecal microbial communities and to track the colonization of potential pathogens.

7.2. Pathogen Detection and Surveillance

Amplicon-based rRNA profiling can detect the presence of pathogenic taxa even when they are present at low abundance, as the universal primer approach does not require prior knowledge of the target. In veterinary diagnostic settings, this allows screening for bacterial pathogens such as Salmonella spp., Campylobacter spp., and Clostridium spp. without the need for culture or specific PCR assays. The computational approaches to understanding antimicrobial resistance (AMR) can be integrated with rRNA profiling to link taxonomic changes with resistance gene carriage.

7.3. Limitations and Considerations

A major limitation of rRNA-based profiling is the variable copy number of rRNA operons across bacterial species, ranging from 1 to 15 or more copies per genome. This copy number variation biases abundance estimates, as organisms with more rDNA copies are overrepresented in amplicon data. Computational normalization methods exist but cannot fully correct this bias. Additionally, horizontal gene transfer of rRNA genes, although rare, can confound phylogenetic inference [25]. Finally, the resolution of 16S rRNA is often insufficient to discriminate between closely related species or strains, necessitating the use of alternative markers such as the ITS region or whole-genome sequencing for fine-scale discrimination [35].

8. Conclusion

Ribosomal RNA is a molecule of extraordinary biological significance, serving as the catalytic and structural core of the ribosome while also providing the most widely used phylogenetic marker in microbial ecology. Its dual role as a functional RNA and as a molecular fossil of evolutionary history makes it an object of study across biochemistry, molecular biology, and bioinformatics. The chemical diversity of rRNA, conferred by site-specific modifications and variant alleles, contributes to ribosome heterogeneity and fine-tunes translational output. In metagenomics, rRNA amplicon and shotgun sequencing approaches enable comprehensive profiling of microbial communities from veterinary and environmental samples. Continued improvements in sequencing technology, reference databases such as SILVA [16, 14, 7], and computational tools for alignment [17] and classification will further extend the power of rRNA-based analyses.

References

[1] Hori Y, Engel C, Kobayashi T. Regulation of ribosomal RNA gene copy number, transcription and nucleolus organization in eukaryotes. Nat Rev Mol Cell Biol. 2023. https://www.semanticscholar.org/paper/e03f0724d16d7fd5b38a05bf0a462f2ef834d508

[2] Parks MM, Kurylo CM, Dass R, et al. Variant ribosomal RNA alleles are conserved and exhibit tissue-specific expression. Sci Adv. 2018. https://www.semanticscholar.org/paper/9b7f050267821a64cc357f5deb0f8631abf6a6e4

[3] White T. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. Journal. 1990. https://www.semanticscholar.org/paper/ea0d088835e68c5ef05556a03e8812e26a2d8380

[4] Lane DJ, Pace B, Olsen GJ, et al. Rapid determination of 16S ribosomal RNA sequences for phylogenetic analyses. Proc Natl Acad Sci USA. 1985. https://www.semanticscholar.org/paper/f3e09ef23ec09015a1587123744c247057e38023

[5] Edwards U, Rogall T, Blöcker H, et al. Isolation and direct complete nucleotide determination of entire genes. Characterization of a gene coding for 16S ribosomal RNA. Nucleic Acids Res. 1989. https://www.semanticscholar.org/paper/a91c295ba98bd5cc6e09036f0d8ca9dbd1126242

[6] Klindworth A, Pruesse E, Schweer T, et al. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2012. https://www.semanticscholar.org/paper/2033072bcbe5655b243700a0a65134e9dc1b58d0

[7] Glöckner FO, Yilmaz P, Quast C, et al. 25 years of serving the community with ribosomal RNA gene reference databases and tools. J Biotechnol. 2017. https://www.semanticscholar.org/paper/3a55328adb3e0d3dec70edb90e7244eb29b82f41

[8] Lagesen K, Hallin P, Rødland E, et al. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007. https://www.semanticscholar.org/paper/51902625d9265e9aa897e1d2b0b6ef3f65dde23f

[9] Shine J, Dalgarno L. The 3'-terminal sequence of Escherichia coli 16S ribosomal RNA: complementarity to nonsense triplets and ribosome binding sites. Proc Natl Acad Sci USA. 1974. https://www.semanticscholar.org/paper/70bb49dd09b7498a50e37649b9d84c920ad1c7b7

[10] Bhatt PR, Scaiola A, Loughran G, et al. Structural basis of ribosomal frameshifting during translation of the SARS-CoV-2 RNA genome. Science. 2020. https://www.semanticscholar.org/paper/6510bcca802110a9ed6a24e80b0470b3ff3348ab

[11] Heissenberger C, Liendl L, Nagelreiter F, et al. Loss of the ribosomal RNA methyltransferase NSUN5 impairs global protein synthesis and normal growth. Nucleic Acids Res. 2019. https://www.semanticscholar.org/paper/3ead754092335e36bdd85b6ed89307a8a176397d

[12] Ma H, et al. N6-Methyladenosine methyltransferase ZCCHC4 mediates ribosomal RNA methylation. Nat Chem Biol. 2018. https://www.semanticscholar.org/paper/f7153d321a71d4babeeabedce51798302a4d638c

[13] Lambert M, Benmoussa A, Provost P. Small Non-Coding RNAs Derived from Eukaryotic Ribosomal RNA. Non-Coding RNA. 2019. https://www.semanticscholar.org/paper/b5a7b12f1b3c4b70f9f6c08f945fa03f9245352b

[14] Pruesse E, Quast C, Knittel K, et al. SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 2007. https://www.semanticscholar.org/paper/07c3f2cbbb351b25d7bea9ca82896c01e3484dba

[15] Edgar RC. Updating the 97% identity threshold for 16S ribosomal RNA OTUs. bioRxiv. 2017. https://www.semanticscholar.org/paper/54f4a6e0e29a4e0c4c0d3108494d14ce6b92c333

[16] Quast C, Pruesse E, Yilmaz P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2012. https://www.semanticscholar.org/paper/b204970b0503a923359bff532726666f5e0e971b

[17] Pruesse E, Peplies J, Glöckner FO. SINA: Accurate high-throughput multiple sequence alignment of ribosomal RNA genes. Bioinformatics. 2012. https://www.semanticscholar.org/paper/c59d5a7d079a88ac4fe1cf06600fba607d4720b6

[18] Chun J, Lee JH, Jung Y, et al. EzTaxon: a web-based tool for

[19] Holdt LM, Stahringer A, Sass K, et al. Circular non-coding RNA ANRIL modulates ribosomal RNA maturation and atherosclerosis in humans. Nat Commun. 2016. https://www.semanticscholar.org/paper/c0499eba8917a6493a77971a9be290c80ef00b6d