Microbiome characterization by high-throughput transfer RNA sequencing and modification analysis

Prokaryotic RNA N1-Methyladenosine Erasers Maintain tRNA m1A Modification Levels in Streptomyces venezuelae

tRNA modifications help maintain tRNA structure and facilitate translation and stress response. Found in all three kingdoms of life, m1A tRNA modification occurs in the T loop of many tRNAs, stabilizes tertiary tRNA structure, and impacts translation. M1A in the T loop is reversible by three mammalian demethylase enzymes, which bypasses the need of turning over the tRNA molecule to adjust its m1A levels in cells. However, no prokaryotic tRNA demethylase enzyme has been identified that acts on endogenous RNA modifications. Using Streptomyces venezuelae as a model organism, we confirmed the presence and quantitative m1A tRNA signatures using mass spectrometry and high-throughput tRNA sequencing. We identified two RNA demethylases that can remove m1A in tRNA and validated the activity of a previously annotated tRNA m1A writer. Using single-gene knockouts of these erasers and the m1A writer, we found dynamic changes of m1A levels in many tRNAs under stress conditions. Phenotypic characterization highlighted changes in their growth and altered antibiotic production. Our identification of the first prokaryotic tRNA demethylase enzyme paves the way for investigating new mechanisms of translational regulation in bacteria.

The modification landscape of Pseudomonas aeruginosa tRNAs

RNA modifications have a substantial impact on tRNA function, with modifications in the anticodon loop contributing to translational fidelity and modifications in the tRNA core impacting structural stability. In bacteria, tRNA modifications are crucial for responding to stress and regulating the expression of virulence factors. Although tRNA modifications are well-characterized in a few model organisms, our knowledge of tRNA modifications in human pathogens, such as Pseudomonas aeruginosa, remains limited. Here, we leveraged two orthogonal approaches to build a reference landscape of tRNA modifications in Escherichia coli, which enabled us to identify similar modifications in P. aeruginosa Our analysis supports a substantial degree of conservation between the two organisms, while also uncovering potential sites of tRNA modification in P. aeruginosa tRNAs that are not present in E. coli The mutational signature at one of these sites, position 46 of tRNAGln1(UUG) is dependent on the P. aeruginosa homolog of TapT, the enzyme responsible for the 3-(3-amino-3-carboxypropyl) uridine (acp3U) modification. Identifying which modifications are present on different tRNAs will uncover the pathways impacted by the different tRNA-modifying enzymes, some of which play roles in determining virulence and pathogenicity.

The modification landscape of P. aeruginosa tRNAs

RNA modifications have a substantial impact on tRNA function, with modifications in the anticodon loop contributing to translational fidelity and modifications in the tRNA core impacting structural stability. In bacteria, tRNA modifications are crucial for responding to stress and regulating the expression of virulence factors. Although tRNA modifications are well-characterized in a few model organisms, our knowledge of tRNA modifications in human pathogens, such as Pseudomonas aeruginosa, remains limited. Here we leveraged two orthogonal approaches to build a reference landscape of tRNA modifications in E. coli, which enabled us to identify similar modifications in P. aeruginosa. Our analysis revealed a substantial degree of conservation between the two organisms, while also uncovering potential sites of tRNA modification in P. aeruginosa tRNAs that are not present in E. coli. The mutational signature at one of these sites, position 46 of tRNAGln1(UUG) is dependent on the P. aeruginosa homolog of TapT, the enzyme responsible for the 3-(3-amino-3-carboxypropyl) uridine (acp3U) modification. Identifying which modifications are present on different tRNAs will uncover the pathways impacted by the different tRNA modifying enzymes, some of which play roles in determining virulence and pathogenicity.

Deciphering Microbe-Mediated Dissolved Organic Matter Reactome in Wastewater Treatment Plants Using Directed Paired Mass Distance

Understanding the reaction mechanism of dissolved organic matter (DOM) during wastewater biotreatment is crucial for optimal DOM control. Here, we develop a directed paired mass distance (dPMD) method that constructs a molecular network displaying the reaction pathways of DOM. It couples direction inference and PMD analysis to extract the substrate-product relationships and delta masses of potentially paired reactants directly from sequential mass spectrometry data without formula assignment. Using this method, we analyze the influent and effluent samples from the bioprocesses of 12 wastewater treatment plants (WWTPs) and build a dPMD network to characterize the core reactome of DOM. The network shows that the first step of the transformation triggers reaction cascades that diversify the DOM, but the highly overlapped subsequent reaction pathways result in similar effluent DOM compositions across WWTPs despite varied influents. Mass changes exhibit consistent gain/loss preferences (e.g., +3.995 and -16.031) but different occurrences across WWTPs. Combined with genome-centric metatranscriptomics, we reveal the associations among dPMDs, enzymes, and microbes. Most enzymes are involved in oxygenation, (de)hydrogenation, demethylation, and hydration-related reactions but with different target substrates and expressed by various taxa, as exemplified by Proteobacteria, Actinobacteria, and Nitrospirae. Therefore, a functionally diverse community is pivotal for advanced DOM degradation.

Precision run-on sequencing (PRO-seq) for microbiome transcriptomics

Bacteria respond to environmental stimuli through precise regulation of transcription initiation and elongation. Bulk RNA sequencing primarily characterizes mature transcripts, so to identify actively transcribed loci we need to capture RNA polymerase (RNAP) complexed with nascent RNA. However, such capture methods have only previously been applied to culturable, genetically tractable organisms such as E. coli and B. subtilis. Here we apply precision run-on sequencing (PRO-seq) to profile nascent transcription in cultured E. coli and diverse uncultured bacteria. We demonstrate that PRO-seq can characterize the transcription of small, structured, or post-transcriptionally modified RNAs, which are often absent from bulk RNA-seq libraries. Applying PRO-seq to the human microbiome highlights taxon-specific RNAP pause motifs and pause-site distributions across non-coding RNA loci that reflect structure-coincident pausing. We also uncover concurrent transcription and cleavage of CRISPR guide RNAs and transfer RNAs. We demonstrate the utility of PRO-seq for exploring transcriptional dynamics in diverse microbial communities.

A tRNA modification in Mycobacterium tuberculosis facilitates optimal intracellular growth

Diverse chemical modifications fine-tune the function and metabolism of tRNA. Although tRNA modification is universal in all kingdoms of life, profiles of modifications, their functions, and physiological roles have not been elucidated in most organisms including the human pathogen, Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis. To identify physiologically important modifications, we surveyed the tRNA of Mtb, using tRNA sequencing (tRNA-seq) and genome-mining. Homology searches identified 23 candidate tRNA modifying enzymes that are predicted to create 16 tRNA modifications across all tRNA species. Reverse transcription-derived error signatures in tRNA-seq predicted the sites and presence of nine modifications. Several chemical treatments prior to tRNA-seq expanded the number of predictable modifications. Deletion of Mtb genes encoding two modifying enzymes, TruB and MnmA, eliminated their respective tRNA modifications, validating the presence of modified sites in tRNA species. Furthermore, the absence of mnmA attenuated Mtb growth in macrophages, suggesting that MnmA-dependent tRNA uridine sulfation contributes to Mtb intracellular growth. Our results lay the foundation for unveiling the roles of tRNA modifications in Mtb pathogenesis and developing new therapeutics against tuberculosis.

Apilactobacillus kunkeei releases RNA-associated membrane vesicles and proteinaceous nanoparticles

Extracellularly released particles, including membrane vesicles, have increasingly been recognized as important for bacterial community functions and host-interaction processes, but their compositions and functional roles differ between species and also between strains of the same species. In this study, we have determined the composition of membrane vesicles and protein particles identified in the cell-free pellets of two strains of Apilactobacillus kunkeei, a defensive symbiont of honeybees. The membrane vesicles were separated from the extracellular particles using density gradient ultracentrifugation. The peaks of the RNA and protein distributions were separated from each other and the highest concentration of RNA was observed in the fractions that contained the membrane vesicles while the highest protein concentration coincided with the fractions that contained extracellular particles. A comparative proteomics analysis by LC-MS/MS showed that 37 proteins with type-I signal peptides were consistently identified across the fractionated samples obtained from the cell-free pellets, of which 29 were orthologs detected in both strains. Functional predictions of the extracellular proteins revealed the presence of glycoside hydrolases, glycosyltransferases, giant proteins and peptidases. The extracellular transcriptomes mapped to a broad set of genes with a similar functional profile as the whole cell transcriptome. This study provides insights into the composition of membrane vesicles and extracellular proteins of a bee-associated symbiont.

A tRNA modification in Mycobacterium tuberculosis facilitates optimal intracellular growth

Diverse chemical modifications fine-tune the function and metabolism of tRNA. Although tRNA modification is universal in all kingdoms of life, profiles of modifications, their functions, and physiological roles have not been elucidated in most organisms including the human pathogen, Mycobacterium tuberculosis ( Mtb ), the causative agent of tuberculosis. To identify physiologically important modifications, we surveyed the tRNA of Mtb , using tRNA sequencing (tRNA-seq) and genome-mining. Homology searches identified 23 candidate tRNA modifying enzymes that are predicted to create 16 tRNA modifications across all tRNA species. Reverse transcription-derived error signatures in tRNA-seq predicted the sites and presence of 9 modifications. Several chemical treatments prior to tRNA-seq expanded the number of predictable modifications. Deletion of Mtb genes encoding two modifying enzymes, TruB and MnmA, eliminated their respective tRNA modifications, validating the presence of modified sites in tRNA species. Furthermore, the absence of mnmA attenuated Mtb growth in macrophages, suggesting that MnmA-dependent tRNA uridine sulfation contributes to Mtb intracellular growth. Our results lay the foundation for unveiling the roles of tRNA modifications in Mtb pathogenesis and developing new therapeutics against tuberculosis.

Fluorescent labeling of tRNA for rapid kinetic interaction studies with tRNA-binding proteins

Transfer RNA (tRNA) plays a critical role during translation and interacts with numerous proteins during its biogenesis, functional cycle and degradation. In particular, tRNA is extensively post-transcriptionally modified by various tRNA modifying enzymes which each target a specific nucleotide at different positions within tRNAs to introduce different chemical modifications. Fluorescent assays can be used to study the interaction between a protein and tRNA. Moreover, rapid mixing fluorescence stopped-flow assays provide insights into the kinetics of the tRNA-protein interaction in order to elucidate the tRNA binding mechanism for the given protein. A prerequisite for these studies is a fluorescently labeled molecule, such as fluorescent tRNA, wherein a change in fluorescence occurs upon protein binding. In this chapter, we discuss the utilization of tRNA modifications in order to introduce fluorophores at particular positions within tRNAs. Particularly, we focus on in vitro thiolation of a uridine at position 8 within tRNAs using the tRNA modification enzyme ThiI, followed by labeling of the thiol group with fluorescein. As such, this fluorescently labeled tRNA is primarily unmodified, with the exception of the thiolation modification to which the fluorophore is attached, and can be used as a substrate to study the binding of different tRNA-interacting factors. Herein, we discuss the example of studying the tRNA binding mechanism of the tRNA modifying enzymes TrmB and DusA using internally fluorescein-labeled tRNA.

Analysis of queuosine and 2-thio tRNA modifications by high throughput sequencing

Queuosine (Q) is a conserved tRNA modification at the wobble anticodon position of tRNAs that read the codons of amino acids Tyr, His, Asn, and Asp. Q-modification in tRNA plays important roles in the regulation of translation efficiency and fidelity. Queuosine tRNA modification is synthesized de novo in bacteria, whereas in mammals the substrate for Q-modification in tRNA is queuine, the catabolic product of the Q-base of gut bacteria. This gut microbiome dependent tRNA modification may play pivotal roles in translational regulation in different cellular contexts, but extensive studies of Q-modification biology are hindered by the lack of high throughput sequencing methods for its detection and quantitation. Here, we describe a periodate-treatment method that enables single base resolution profiling of Q-modification in tRNAs by Nextgen sequencing from biological RNA samples. Periodate oxidizes the Q-base, which results in specific deletion signatures in the RNA-seq data. Unexpectedly, we found that periodate-treatment also enables the detection of several 2-thio-modifications including τm5s2U, mcm5s2U, cmnm5s2U, and s2C by sequencing in human and E. coli tRNA. We term this method periodate-dependent analysis of queuosine and sulfur modification sequencing (PAQS-seq). We assess Q- and 2-thio-modifications at the tRNA isodecoder level, and 2-thio modification changes in stress response. PAQS-seq should be widely applicable in the biological studies of Q- and 2-thio-modifications in mammalian and microbial tRNAs.

Dietary Carbohydrate as Glycemic Load, Not Fat, Coupled with Genetic Permissiveness Favoring Rapid Growth and Extra Calories, Dictate Metabolic Syndrome and Diabetes Induction in Nile Rats (Arvicanthis niloticus)

Objective: Whether dietary carbohydrate (CHO) or fat is more involved in type 2 diabetes (T2DM) induction uncomplicated by dietary fiber was addressed in a spontaneous diabetic model, the diurnal Nile rat that mimics the human condition. Methods: A total of 138 male Nile rats were fed plant-based and animal-based saturated fat where 10% energy as CHO and fat were exchanged across 5 diets keeping protein constant, from 70:10:20 to 20:60:20 as CHO:fat:protein %energy. Diabetes induction was analyzed by: 1. diet composition, i.e., CHO:fat ratio, to study the impact of diet; 2. quintiles of average caloric intake per day to study the impact of calories; 3. quintiles of diabetes severity to study the epigenetic impact on diabetes resistance. Results: High glycemic load (GLoad) was most problematic if coupled with high caloric consumption. Diabetes severity highlighted rapid growth and caloric intake as likely epigenetic factors distorting glucose metabolism. The largest weanling rats ate more, grew faster, and developed more diabetes when the dietary GLoad exceeded their gene-based metabolic capacity for glucose disposal. Diabetes risk increased for susceptible rats when energy intake exceeded 26 kcal/day and the GLoad was >175/2000 kcal of diet and when the diet provided >57% energy as CHO. Most resistant rats ate <25 kcal/day independent of the CHO:fat diet ratio or the GLoad adjusted to body size. Conclusion: Beyond the CHO:fat ratio and GLoad, neither the type of fat nor the dietary polyunsaturated/saturated fatty acid (P/S) ratio had a significant impact, suggesting genetic permissiveness affecting caloric and glucose intake and glucose disposition were key to modulating Nile rat diabetes. Fat became protective by limiting GLoad when it contributed >40% energy and displaced CHO to <50% energy, thereby decreasing the number of diabetic rats and diabetes severity.

Bacterial tRNA landscape revisited

The updated Wobble Hypothesis reasonably explains why some 40 tRNA species are sufficient to decode the 61 amino acid codons of the Universal Genetic Code. However, we still have no clue why eubacteria lack tRNA isoacceptors with ANN anticodons, whereas eukaryotes universally lack eight GNN anticodons, only one of which is also absent in bacteria. Direct tRNA sequencing could resolve the patterns of nucleoside modification that had been driving the divergent evolution in prokaryotes and eukaryotes, but this task will require the development of AI-supported base-callers that can recognize modified nucleosides without any subsequent analytical verification. Our knowledge of the bacterial tRNA landscape is moreover broadened by the recent discovery of antisense tRNAs and tRNA-derived fragments that should be examined in their roles for gene expression, translation, bacterial physiology or metabolism.

TRMT6/61A-dependent base methylation of tRNA-derived fragments regulates gene-silencing activity and the unfolded protein response in bladder cancer

RNA modifications are important regulatory elements of RNA functions. However, most genome-wide mapping of RNA modifications has focused on messenger RNAs and transfer RNAs, but such datasets have been lacking for small RNAs. Here we mapped N1-methyladenosine (m1A) in the cellular small RNA space. Benchmarked with synthetic m1A RNAs, our workflow identified specific groups of m1A-containing small RNAs, which are otherwise disproportionally under-represented. In particular, 22-nucleotides long 3' tRNA-fragments are highly enriched for TRMT6/61A-dependent m1A located within the seed region. TRMT6/61A-dependent m1A negatively affects gene silencing by tRF-3s. In urothelial carcinoma of the bladder, where TRMT6/61A is over-expressed, higher m1A modification on tRFs is detected, correlated with a dysregulation of tRF targetome. Lastly, TRMT6/61A regulates tRF-3 targets involved in unfolded protein response. Together, our results reveal a mechanism of regulating gene expression via base modification of small RNA.

tRNA modification dynamics from individual organisms to metaepitranscriptomics of microbiomes

tRNA is the most extensively modified RNA in cells. On average, a bacterial tRNA contains 8 modifications per molecule and a eukaryotic tRNA contains 13 modifications per molecule. Recent studies reveal that tRNA modifications are highly dynamic and respond extensively to environmental conditions. Functions of tRNA modification dynamics include enhanced, on-demand decoding of specific codons in response genes and regulation of tRNA fragment biogenesis. This review summarizes recent advances in the studies of tRNA modification dynamics in biological processes, tRNA modification erasers, and human-associated bacteria. Furthermore, we use the term "metaepitranscriptomics" to describe the potential and approach of tRNA modification studies in natural biological communities such as microbiomes. tRNA is highly modified in cells, and tRNA modifications respond extensively to environmental conditions to enhance translation of specific genes and produce tRNA fragments on demand. We review recent advances in tRNA sequencing methods, tRNA modification dynamics in biological processes, and tRNA modification studies in natural communities such as the microbiomes.

Site-Specific Profiling of 4-Thiouridine Across Transfer RNA Genes in Escherichia coli

The transfer RNA (tRNA) modification 4-thiouridine (s4U) acts as a near-ultraviolet (UVA) radiation sensor in Escherichia coli (E. coli), where it induces a growth delay upon exposure to the UVA radiation (∼310-400 nm). Herein, we report sequencing methodology for site-specific profiling of s4U modification in E. coli tRNAs. Upon the addition of iodoacetamide (IA) or iodoacetyl-PEG2-biotin (BIA), the nucleophilic sulfur of s4U forms a reaction product that is extensively characterized by liquid chromatography-mass spectrometry (LC-MS/MS) analysis. This method is readily applied to the alkylation of natively occurring s4U on E. coli tRNA. Next-generation sequencing of BIA-treated tRNA from E. coli revealed misincorporations at position 8 in 19 of the 20 amino acid tRNA species. Alternatively, tRNA from the ΔthiI strain, which cannot introduce the s4U modification, does not exhibit any misincorporation at the corresponding positions, directly linking the base transitions and the tRNA modification. Independently, the s4U modification on E. coli tRNA was further validated by LC-MS/MS sequencing. Nuclease digestion of wild-type and deletion strains E. coli tRNA with RNase T1 generated smaller s4U/U containing fragments that could be analyzed by MS/MS analysis for modification assignment. Furthermore, RNase T1 digestion of tRNAs treated either with IA or BIA showed the specificity of iodoacetamide reagents toward s4U in the context of complex tRNA modifications. Overall, these results demonstrate the utility of the alkylation of s4U in the site-specific profiling of the modified base in native cellular tRNA.

Riddle of the Sphinx: Emerging Role of Transfer RNAs in Human Cancer

The dysregulation of transfer RNA (tRNA) expression contributes to the diversity of proteomics, heterogeneity of cell populations, and instability of the genome, which may be related to human cancer susceptibility. However, the relationship between tRNA dysregulation and cancer susceptibility remains elusive because the landscape of cancer-associated tRNAs has not been portrayed yet. Furthermore, the molecular mechanisms of tRNAs involved in tumorigenesis and cancer progression have not been systematically understood. In this review, we detail current knowledge of cancer-related tRNAs and comprehensively summarize the basic characteristics and functions of these tRNAs, with a special focus on their role and involvement in human cancer. This review bridges the gap between tRNAs and cancer and broadens our understanding of their relationship, thus providing new insights and strategies to improve the potential clinical applications of tRNAs for cancer diagnosis and therapy.

Cell type-specific analysis by single-cell profiling identifies a stable mammalian tRNA-mRNA interface and increased translation efficiency in neurons

The correlation between codon and anticodon pools influences the efficiency of translation, but whether differences exist in these pools across individual cells is unknown. We determined that codon usage and amino acid demand are highly stable across different cell types using available mouse and human single-cell RNA-sequencing atlases. After showing the robustness of ATAC-sequencing measurements for the analysis of tRNA gene usage, we quantified anticodon usage and amino acid supply in both mouse and human single-cell ATAC-seq atlases. We found that tRNA gene usage is overall coordinated across cell types, except in neurons, which clustered separately from other cell types. Integration of these data sets revealed a strong and statistically significant correlation between amino acid supply and demand across almost all cell types. Neurons have an enhanced translation efficiency over other cell types, driven by an increased supply of tRNA^Ala (AGC) anticodons. This results in faster decoding of the Ala-GCC codon, as determined by cell type–specific ribosome profiling, suggesting that the reduction of tRNA^Ala (AGC) anticodon pools may be implicated in neurological pathologies. This study, the first such examination of codon usage, anticodon usage, and translation efficiency resolved at the cell-type level with single-cell information, identifies a conserved landscape of translation elongation across mammalian cellular diversity and evolution.

The occurrence order and cross-talk of different tRNA modifications

Chemical modifications expand the composition of RNA molecules from four standard nucleosides to over 160 modified nucleosides, which greatly increase the complexity and utility of RNAs. Transfer RNAs (tRNAs) are the most heavily modified cellular RNA molecules and contain the largest variety of modifications. Modification of tRNAs is pivotal for protein synthesis and also precisely regulates the noncanonical functions of tRNAs. Defects in tRNA modifications lead to numerous human diseases. Up to now, more than 100 types of modifications have been found in tRNAs. Intriguingly, some modifications occur widely on all tRNAs, while others only occur on a subgroup of tRNAs or even only a specific tRNA. The modification frequency of each tRNA is approximately 7% to 25%, with 5-20 modification sites present on each tRNA. The occurrence and modulation of tRNA modifications are specifically noticeable as plenty of interplays among different sites and modifications have been discovered. In particular, tRNA modifications are responsive to environmental changes, indicating their dynamic and highly organized nature. In this review, we summarized the known occurrence order, cross-talk, and cooperativity of tRNA modifications.

Coping with darkness: The adaptive response of marine picocyanobacteria to repeated light energy deprivation

The picocyanobacteria Prochlorococcus and Synechococcus are found throughout the ocean's euphotic zone, where the daily light:dark cycle drives their physiology. Periodic deep mixing events can, however, move cells below this region, depriving them of light for extended periods of time. Here, we demonstrate that members of these genera can adapt to tolerate repeated periods of light energy deprivation. Strains kept in the dark for 3 d and then returned to the light initially required 18-26 d to resume growth, but after multiple rounds of dark exposure they began to regrow after only 1-2 d. This dark-tolerant phenotype was stable and heritable; some cultures retained the trait for over 132 generations even when grown in a standard 13:11 light:dark cycle. We found no genetic differences between the dark-tolerant and parental strains of Prochlorococcus NATL2A, indicating that an epigenetic change is likely responsible for the adaptation. To begin to explore this possibility, we asked whether DNA methylation-one potential mechanism mediating epigenetic inheritance in bacteria-occurs in Prochlorococcus. LC-MS/MS analysis showed that while DNA methylations, including 6 mA and 5 mC, are found in some other Prochlorococcus strains, there were no methylations detected in either the parental or dark-tolerant NATL2A strains. These findings suggest that Prochlorococcus utilizes a yet-to-be-determined epigenetic mechanism to adapt to the stress of extended light energy deprivation, and highlights phenotypic heterogeneity as an additional dimension of Prochlorococcus diversity.

Adaptor Molecules Epitranscriptome Reprograms Bacterial Pathogenicity

The strong decoration of tRNAs with post-transcriptional modifications provides an unprecedented adaptability of this class of non-coding RNAs leading to the regulation of bacterial growth and pathogenicity. Accumulating data indicate that tRNA post-transcriptional modifications possess a central role in both the formation of bacterial cell wall and the modulation of transcription and translation fidelity, but also in the expression of virulence factors. Evolutionary conserved modifications in tRNA nucleosides ensure the proper folding and stability redounding to a totally functional molecule. However, environmental factors including stress conditions can cause various alterations in tRNA modifications, disturbing the pathogen homeostasis. Post-transcriptional modifications adjacent to the anticodon stem-loop, for instance, have been tightly linked to bacterial infectivity. Currently, advances in high throughput methodologies have facilitated the identification and functional investigation of such tRNA modifications offering a broader pool of putative alternative molecular targets and therapeutic avenues against bacterial infections. Herein, we focus on tRNA epitranscriptome shaping regarding modifications with a key role in bacterial infectivity including opportunistic pathogens of the human microbiome.

Lemon exosome-like nanoparticles enhance stress survival of gut bacteria by RNase P-mediated specific tRNA decay

Diet and bile play critical roles in shaping gut microbiota, but the molecular mechanism underlying interplay with intestinal microbiota is unclear. Here, we showed that lemon-derived exosome-like nanoparticles (LELNs) enhance lactobacilli toleration to bile. To decipher the mechanism, we used Lactobacillus rhamnosus GG (LGG) as proof of concept to show that LELNs enhance LGG bile resistance via limiting production of Msp1 and Msp3, resulting in decrease of bile accessibility to cell membrane. Furthermore, we found that decline of Msps protein levels was regulated through specific tRNAser UCC and tRNAser UCG decay. We identified RNase P, an essential housekeeping endonuclease, being responsible for LELNs-induced tRNAser UCC and tRNAser UCG decay. We further identified galacturonic acid-enriched pectin-type polysaccharide as the active factor in LELNs to increase bile resistance and downregulate tRNAser UCC and tRNAser UCG level in the LGG. Our study demonstrates a tRNA-based gene expression regulation mechanism among lactobacilli to increase bile resistance.

Shaping the Bacterial Epitranscriptome-5'-Terminal and Internal RNA Modifications

All domains of life utilize a diverse set of modified ribonucleotides that can impact the sequence, structure, function, stability, and the fate of RNAs, as well as their interactions with other molecules. Today, more than 160 different RNA modifications are known that decorate the RNA at the 5'-terminus or internal RNA positions. The boost of next-generation sequencing technologies sets the foundation to identify and study the functional role of RNA modifications. The recent advances in the field of RNA modifications reveal a novel regulatory layer between RNA modifications and proteins, which is central to developing a novel concept called "epitranscriptomics." The majority of RNA modifications studies focus on the eukaryotic epitranscriptome. In contrast, RNA modifications in prokaryotes are poorly characterized. This review outlines the current knowledge of the prokaryotic epitranscriptome focusing on mRNA modifications. Here, it is described that several internal and 5'-terminal RNA modifications either present or likely present in prokaryotic mRNA. Thereby, the individual techniques to identify these epitranscriptomic modifications, their writers, readers and erasers, and their proposed functions are explored. Besides that, still unanswered questions in the field of prokaryotic epitranscriptomics are pointed out, and its future perspectives in the dawn of next-generation sequencing technologies are outlined.

Effects of compound deer bone extract on osteoporosis model mice and intestinal microflora

The preventive and therapeutic mechanisms of CDBE on osteoporosis were studied by observing the serum bone-related biochemical indicators, bone trabecular micro-structure and intestinal flora in ovariectomized osteoporosis model mice, in order to provide a scientific theoretical basis for the further study on the effect of CDBE on osteoporosis, and the prevention and treatment of osteoporosis with clinical traditional Chinese medicines. The components in CDBE were detected by UHPLC-MS. A mouse osteoporosis model was established by the bilateral ovariectomy in female ICR mice. The biochemical indicators related to osteoporosis were detected, the right proximal tibia was scanned by Micro-CT, the intestinal microflora in the colon contents were examined, and the changes of microflora were taken as the main target to evaluate the effect of CDBE on the intestinal microflora in the model mice. A total of 16 compounds were obtained by the combined application of UHPLC-MS. CDBE could significantly increase the contents of E2, Ca²⁺ , CT, HyP, OCN, FOXP3, P1NP and CTX-II, in the model mice. CDBE could significantly improve the trabecular micro-structure, Tb.N, Tb.Sp, SMI and Conn.D. CDBE could make the intestinal flora of osteoporosis model mice tend to healthy mice in species and quantity. CDBE can improve the symptoms of postmenopausal osteoporosis in mice, with a positive effect on the intestinal flora of postmenopausal mice. Its mechanism of regulating the symptoms of osteoporosis may be related to the regulation of bone-related biochemical indicators in the serum of mice. PRACTICAL APPLICATIONS: This research has a positive impact on the development of functional food with anti-osteoporosis in the future.

Tissue-specific reprogramming of host tRNA transcriptome by the microbiome

Transfer RNAs (tRNAs) are essential for translation, and tRNA expression and modifications are regulated by many factors. However, the interplay between the microbiome and host tRNA profiles through host-microbiome interactions has not been explored. In this study, we investigated host-microbiome interactions via the tRNA profiling of four tissue types from germ-free and specific pathogen-free mice. Our analyses reveal that cytosolic and mitochondrial tRNA expression and tRNA modifications in the host are reprogrammed in a tissue-specific and microbiome-dependent manner. In terms of tRNA expression, the intestines and brains are more sensitive to the influence of the microbiome than the livers and kidneys. In terms of tRNA modifications, cytosolic tRNAs show more obvious changes in the livers and kidneys in the presence of the microbiome. Our findings reveal a previously unexplored relationship among the microbiome, tRNA abundance, and epitranscriptome in a mammalian host.

CaptureSeq: Hybridization-Based Enrichment of cpn60 Gene Fragments Reveals the Community Structures of Synthetic and Natural Microbial Ecosystems

Background. The molecular profiling of complex microbial communities has become the basis for examining the relationship between the microbiome composition, structure and metabolic functions of those communities. Microbial community structure can be partially assessed with “universal” PCR targeting taxonomic or functional gene markers. Increasingly, shotgun metagenomic DNA sequencing is providing more quantitative insight into microbiomes. However, both amplicon-based and shotgun sequencing approaches have shortcomings that limit the ability to study microbiome dynamics. Methods. We present a novel, amplicon-free, hybridization-based method (CaptureSeq) for profiling complex microbial communities using probes based on the chaperonin-60 gene. Molecular profiles of a commercially available synthetic microbial community standard were compared using CaptureSeq, whole metagenome sequencing, and 16S universal target amplification. Profiles were also generated for natural ecosystems including antibiotic-amended soils, manure storage tanks, and an agricultural reservoir. Results. The CaptureSeq method generated a microbial profile that encompassed all of the bacteria and eukaryotes in the panel with greater reproducibility and more accurate representation of high G/C content microorganisms compared to 16S amplification. In the natural ecosystems, CaptureSeq provided a much greater depth of coverage and sensitivity of detection compared to shotgun sequencing without prior selection. The resulting community profiles provided quantitatively reliable information about all three domains of life (Bacteria, Archaea, and Eukarya) in the different ecosystems. The applications of CaptureSeq will facilitate accurate studies of host-microbiome interactions for environmental, crop, animal and human health. Conclusions: cpn60-based hybridization enriched for taxonomically informative DNA sequences from complex mixtures. In synthetic and natural microbial ecosystems, CaptureSeq provided sequences from prokaryotes and eukaryotes simultaneously, with quantitatively reliable read abundances. CaptureSeq provides an alternative to PCR amplification of taxonomic markers with deep community coverage while minimizing amplification biases.

Analysis of RNA Modifications by Second- and Third-Generation Deep Sequencing: 2020 Update

The precise mapping and quantification of the numerous RNA modifications that are present in tRNAs, rRNAs, ncRNAs/miRNAs, and mRNAs remain a major challenge and a top priority of the epitranscriptomics field. After the keystone discoveries of massive m⁶A methylation in mRNAs, dozens of deep sequencing-based methods and protocols were proposed for the analysis of various RNA modifications, allowing us to considerably extend the list of detectable modified residues. Many of the currently used methods rely on the particular reverse transcription signatures left by RNA modifications in cDNA; these signatures may be naturally present or induced by an appropriate enzymatic or chemical treatment. The newest approaches also include labeling at RNA abasic sites that result from the selective removal of RNA modification or the enhanced cleavage of the RNA ribose-phosphate chain (perhaps also protection from cleavage), followed by specific adapter ligation. Classical affinity/immunoprecipitation-based protocols use either antibodies against modified RNA bases or proteins/enzymes, recognizing RNA modifications. In this survey, we review the most recent achievements in this highly dynamic field, including promising attempts to map RNA modifications by the direct single-molecule sequencing of RNA by nanopores.

The biological process of lysine-tRNA charging is therapeutically targetable in liver cancer

BACKGROUND & AIMS

Mature transfer RNAs (tRNA) charged with amino acids decode mRNA to synthesize proteins. Dysregulation of translational machineries has a fundamental impact on cancer biology. This study aims to map the tRNAome landscape in liver cancer patients and to explore potential therapeutic targets at the interface of charging amino acid with tRNA.

METHODS

Resected tumour and paired tumour-free (TFL) tissues from hepatocellular carcinoma (HCC) patients (n = 69), and healthy liver tissues from organ transplant donors (n = 21), HCC cell lines, and cholangiocarcinoma (CC) patient-derived tumour organoids were used.

RESULTS

The expression levels of different mature tRNAs were highly correlated and closely clustered within individual tissues, suggesting that different members of the tRNAome function cooperatively in protein translation. Interestingly, high expression of tRNA-Lys-CUU in HCC tumours was associated with more tumour recurrence (HR 1.1; P = .022) and worse patient survival (HR 1.1; P = .0037). The expression of Lysyl-tRNA Synthetase (KARS), the enzyme catalysing the charge of lysine to tRNA-Lys-CUU, was significantly upregulated in HCC tumour tissues compared to tumour-free liver tissues. In HCC cell lines, lysine deprivation, KARS knockdown or treatment with the KARS inhibitor cladosporin effectively inhibited overall cell growth, single cell-based colony formation and cell migration. This was mechanistically mediated by cell cycling arrest and induction of apoptosis. Finally, these inhibitory effects were confirmed in 3D cultured patient-derived CC organoids.

CONCLUSIONS

The biological process of charging tRNA-Lys-CUU with lysine sustains liver cancer cell growth and migration, and is clinically relevant in HCC patients. This process can be therapeutically targeted and represents an unexplored territory for developing novel treatment strategies against liver cancer.

Comprehensive analysis of the lysine acetylome and succinylome in the hippocampus of gut microbiota-dysbiosis mice

Introduction

Major depressive disorder is caused by gene–environment interactions, and the host microbiome has been recognized as an important environmental factor. However, the underlying mechanisms of the host–microbiota interactions that lead to depression are complex and remain poorly understood.

Objectives

The present study aimed to explore the possible mechanisms underlying gut microbiota dysbiosis-induced depressive-like behaviors.

Methods

We used high-performance liquid chromatography-tandem mass spectrometry to analyze alterations in the hippocampal lysine acetylome and succinylome in male mice that had received gut microbiota from fecal samples of either patients with major depressive disorder or healthy controls. This was followed by bioinformatic analyses.

Results

A total of 315 acetylation sites on 223 proteins and 624 succinylation sites on 494 proteins were differentially expressed in the gut microbiota-dysbiosis mice. The significantly acetylated proteins were primarily associated with carbon metabolism disruption and gene transcription suppression, while the synaptic vesicle cycle and protein translation were the most significantly altered functions for succinylated proteins. Additionally, our findings suggest that gut microbiota dysbiosis disturbs mitochondria-mediated biological processes and the MAPK signaling pathway through crosstalk between acetylation and succinylation on relevant proteins.

Conclusions

This is the first study to demonstrate modifications in acetylation and succinylation in gut microbiota-dysbiosis mice. Our findings provide new avenues for exploring the pathogenesis of gut microbiota dysbiosis-related depression, and highlight potential targets for depression treatment.

Noncanonical Roles of tRNAs: tRNA Fragments and Beyond

As one of the most abundant and conserved RNA species, transfer RNAs (tRNAs) are well known for their role in reading the codons on messenger RNAs and translating them into proteins. In this review, we discuss the noncanonical functions of tRNAs. These include tRNAs as precursors to novel small RNA molecules derived from tRNAs, also called tRNA-derived fragments, that are abundant across species and have diverse functions in different biological processes, including regulating protein translation, Argonaute-dependent gene silencing, and more. Furthermore, the role of tRNAs in biosynthesis and other regulatory pathways, including nutrient sensing, splicing, transcription, retroelement regulation, immune response, and apoptosis, is reviewed. Genome organization and sequence variation of tRNA genes are also discussed in light of their noncanonical functions. Lastly, we discuss the recent applications of tRNAs in genome editing and microbiome sequencing.

The Oral Microbiome of Healthy Japanese People at the Age of 90

For a healthy oral cavity, maintaining a healthy microbiome is essential. However, data on healthy microbiomes are not sufficient. To determine the nature of the core microbiome, the oral-microbiome structure was analyzed using pyrosequencing data. Saliva samples were obtained from healthy 90-year-old participants who attended the 20-year follow-up Niigata cohort study. A total of 85 people participated in the health checkups. The study population consisted of 40 male and 45 female participants. Stimulated saliva samples were obtained by chewing paraffin wax for 5 min. The V3–V4 hypervariable regions of the 16S ribosomal RNA (rRNA) gene were amplified by PCR. Pyrosequencing was performed using MiSeq. Operational taxonomic units (OTUs) were assigned on the basis of a 97% identity search in the EzTaxon-e database. Using the threshold of 100% detection on the species level, 13 species were detected: Streptococcus sinensis, Streptococcus pneumoniae, Streptococcus salivarius, KV831974_s, Streptococcus parasanguinis, Veillonella dispar, Granulicatella adiacens, Streptococcus_uc, Streptococcus peroris, KE952139_s, Veillonella parvula, Atopobium parvulum, and AFQU_vs. These species represent potential candidates for the core make-up of the human microbiome.

Identification of Biomarkers to Construct a Competing Endogenous RNA Network and Establishment of a Genomic-Clinicopathologic Nomogram to Predict Survival for Children with Rhabdoid Tumors of the Kidney

Rhabdoid tumor of the kidney (RTK) is a rare and severely malignant tumor occurring in infancy and early childhood, with the overall outcomes remain poor. Neither gene regulatory networks nor biomarkers to predict the prognostic outcomes have been elucidated in RTK. In this study, RNA sequencing data were obtained to identify differentially expressed messenger RNAs (mRNAs), long noncoding RNAs (lncRNAs), and microRNAs (miRNAs) between RTK samples and normal samples. A total of 4217 mRNAs, 284 lncRNAs, and 286 miRNAs were screened out. Of those, 103 mRNAs, 80 lncRNAs, and 45 miRNAs were identified for a competing endogenous RNA (ceRNA) regulatory network, in which three significant modules were identified. A protein-protein interaction (PPI) network was constructed, and the hub-gene cluster consisted of four core genes (EXOSC2, PAK1IP1, WDR43, and POLR1D) was selected. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were also performed to analyze the functional characteristics of differentially expressed mRNAs. Subsequently, among 211 mRNAs, 8 lncRNAs, and 12 miRNAs associated with overall survival (OS) obtained by univariate Cox analysis, 5 mRNAs, 7 lncRNAs, and 7 miRNAs were identified and the risk score formulas were constructed correspondingly using the least absolute shrinkage and selection operator (LASSO) Cox regression model analysis. The log-rank tests and Kaplan-Meier analyses were performed to confirm the predictive value of the risk scores for OS in RTK patients. A genomic-clinicopathologic nomogram integrating the stage and risk scores based on RNAs was established and demonstrated high predictive accuracy and clinical value, which was validated through calibration curves, time-dependent receiver operating characteristic (ROC) curve analyses, and decision curve analysis (DCA). In conclusion, this study not only provided potential insights into the mechanisms underlying RTK, but also presented a practicable tool for predicting the prognosis in children with RTK.

Functions of Bacterial tRNA Modifications: From Ubiquity to Diversity

Modified nucleotides in tRNA are critical components of the translation apparatus, but their importance in the process of translational regulation had until recently been greatly overlooked. Two breakthroughs have recently allowed a fuller understanding of the importance of tRNA modifications in bacterial physiology. One is the identification of the full set of tRNA modification genes in model organisms such as Escherichia coli K12. The second is the improvement of available analytical tools to monitor tRNA modification patterns. The role of tRNA modifications varies greatly with the specific modification within a given tRNA and with the organism studied. The absence of these modifications or reductions can lead to cell death or pleiotropic phenotypes or may have no apparent visible effect. By linking translation through their decoding functions to metabolism through their biosynthetic pathways, tRNA modifications are emerging as important components of the bacterial regulatory toolbox.

Probing the diversity and regulation of tRNA modifications

Transfer RNAs (tRNAs) are non-coding RNAs essential for protein synthesis. tRNAs are heavily decorated with a variety of post-transcriptional modifications (tRNA modifications). Recent methodological advances provide new tools for rapid profiling of tRNA modifications and have led to discoveries of novel modifications and their regulation. Here, we provide an overview of the techniques for investigating tRNA modifications and of the expanding knowledge of their chemistry and regulation.

Comparative tRNA sequencing and RNA mass spectrometry for surveying tRNA modifications

Chemical modifications of the nucleosides that comprise transfer RNAs are diverse. However, the structure, location and extent of modifications have been systematically charted in very few organisms. Here, we describe an approach in which rapid prediction of modified sites through reverse transcription-derived signatures in high-throughput transfer RNA-sequencing (tRNA-seq) data is coupled with identification of tRNA modifications through RNA mass spectrometry. Comparative tRNA-seq enabled prediction of several Vibrio cholerae modifications that are absent from Escherichia coli and also revealed the effects of various environmental conditions on V. cholerae tRNA modification. Through RNA mass spectrometric analyses, we showed that two of the V. cholerae-specific reverse transcription signatures reflected the presence of a new modification (acetylated acp³U (acacp³U)), while the other results from C-to-Ψ RNA editing, a process not described before. These findings demonstrate the utility of this approach for rapid surveillance of tRNA modification profiles and environmental control of tRNA modification.

Adaptation of codon usage to tRNA I34 modification controls translation kinetics and proteome landscape

Codon usage bias is a universal feature of all genomes and plays an important role in regulating protein expression levels. Modification of adenosine to inosine at the tRNA anticodon wobble position (I34) by adenosine deaminases (ADATs) is observed in all eukaryotes and has been proposed to explain the correlation between codon usage and tRNA pool. However, how the tRNA pool is affected by I34 modification to influence codon usage-dependent gene expression is unclear. Using Neurospora crassa as a model system, by combining molecular, biochemical and bioinformatics analyses, we show that silencing of adat2 expression severely impaired the I34 modification levels for the ADAT-related tRNAs, resulting in major ADAT-related tRNA profile changes and reprogramming of translation elongation kinetics on ADAT-related codons. adat2 silencing also caused genome-wide codon usage-biased ribosome pausing on mRNAs and proteome landscape changes, leading to selective translational repression or induction of different mRNAs. The induced expression of CPC-1, the Neurospora ortholog of yeast GCN4p, mediates the transcriptional response after adat2 silencing and amino acid starvation. Together, our results demonstrate that the tRNA I34 modification by ADAT plays a major role in driving codon usage-biased translation to shape proteome landscape.

Modification of transfer RNA (tRNA) can have profound impacts on gene expression by shaping cellular tRNA pool. How codon usage bias and tRNA profiles synergistically regulate gene expression is unclear. By combining molecular, biochemical and bioinformatics analyses, we showed that the correlation between genome codon usage and tRNA I34 (inosine 34) modification modulates translation elongation kinetics and proteome landscape. Inhibition of tRNA I34 modification causes codon usage-dependent ribosome pausing on mRNAs during translation and changes cellular protein contents in a codon usage biased manner. Together, our results demonstrate that the tRNA I34 modification plays a major role in driving codon usage-dependent translation to determine proteome landscape in a eukaryotic organism.

Universally high transcript error rates in bacteria

Errors can occur at any level during the replication and transcription of genetic information. Genetic mutations derived mainly from replication errors have been extensively studied. However, fundamental details of transcript errors, such as their rate, molecular spectrum, and functional effects, remain largely unknown. To globally identify transcript errors, we applied an adapted rolling-circle sequencing approach to Escherichia coli, Bacillus subtilis, Agrobacterium tumefaciens, and Mesoplasma florum, revealing transcript-error rates 3 to 4 orders of magnitude higher than the corresponding genetic mutation rates. The majority of detected errors would result in amino-acid changes, if translated. With errors identified from 9929 loci, the molecular spectrum and distribution of errors were uncovered in great detail. A G→A substitution bias was observed in M. florum, which apparently has an error-prone RNA polymerase. Surprisingly, an increased frequency of nonsense errors towards the 3' end of mRNAs was observed, suggesting a Nonsense-Mediated Decay-like quality-control mechanism in prokaryotes.

sMETASeq: Combined Profiling of Microbiota and Host Small RNAs

Understanding microbial communities' roles in human health and disease requires methods that accurately characterize the microbial composition and their activity and effects within human biological samples. We present sMETASeq (small RNA Metagenomics by Sequencing), a novel method that uses sequencing of small RNAs to jointly measure host small RNA expression and create metagenomic profiles and detect small bacterial RNAs. We evaluated the performance of sMETASeq on a mock bacterial community and demonstrated its use on different human samples, including colon cancer, oral leukoplakia, cervix cancer, and a panel of human biofluids. In all datasets, the detected microbes reflected the biology of the different sample types.

A simplified qPCR method revealing tRNAome remodeling upon infection by genotype 3 hepatitis E virus

The landscape of tRNA-viral codons regulates viral adaption at the translational level, presumably through adapting to host codon usage or modulating the host tRNAome. We found that the major zoonotic genotype of hepatitis E virus (HEV) has not adapted to host codon usage, prompting exploration of the effects of HEV infection on the host tRNAome. However, tRNAome quantification is largely impeded by the extremely short sequences of tRNAs and redundancy of tRNA genes. Here, we present a length-extension and stepwise simplified qPCR method that utilizes a universal DNA/RNA hybrid tRNA adaptor and degenerate primers. Using this novel methodology, we observe that HEV infection dramatically reprograms the hepatic tRNAome, which is likely to facilitate translation of viral RNAs. This tRNAome quantification method bears broad implications for future tRNA research and possibly tRNA-based diagnostics.

LOTTE-seq (Long hairpin oligonucleotide based tRNA high-throughput sequencing): specific selection of tRNAs with 3'-CCA end for high-throughput sequencing

Transfer RNAs belong to the most abundant type of ribonucleic acid in the cell, and detailed investigations revealed correlations between alterations in the tRNA pool composition and certain diseases like breast cancer. However, currently available methods do not sample the entire tRNA pool or lack specificity for tRNAs. A specific disadvantage of such methods is that only full-length tRNAs are analysed, while tRNA fragments or incomplete cDNAs due to RT stops at modified nucleosides are lost. Another drawback in certain approaches is that the tRNA fraction has to be isolated and separated from high molecular weight RNA, resulting in considerable labour costs and loss of material. Based on a hairpin-shaped adapter oligonucleotide selective for tRNA transcripts, we developed a highly specific protocol for efficient and comprehensive high-throughput analysis of tRNAs that combines the benefits of existing methods and eliminates their disadvantages. Due to a 3'-TGG overhang, the adapter is specifically ligated to the tRNA 3'-CCA end. Reverse transcription prior to the ligation of a second adapter allows to include prematurely terminated cDNA products, increasing the number of tRNA reads. This strategy renders this approach a powerful and universal tool to analyse the tRNA pool of cells and organisms under different conditions in health and disease.

Competitive endogenous RNA (ceRNA) regulation network of lncRNAs, miRNAs, and mRNAs in Wilms tumour

Background

Competitive endogenous RNAs (ceRNAs) have revealed a new mechanism of interaction between RNAs. However, an understanding of the ceRNA regulatory network in Wilms tumour (WT) remains limited.

Methods

The expression profiles of mRNAs, miRNAs and lncRNAs in Wilms tumour samples and normal samples were obtained from the Therapeutically Applicable Research to Generate Effective Treatment (TARGET) database. The EdgeR package was employed to identify differentially expressed lncRNAs, miRNAs and mRNAs. Functional enrichment analyses via the ClusterProfile R package were performed, and the lncRNA–miRNA–mRNA interaction ceRNA network was established in Cytoscape. Subsequently, the correlation between the ceRNA network and overall survival was analysed.

Results

A total of 2037 lncRNAs, 154 miRNAs and 3609 mRNAs were identified as differentially expressed RNAs in Wilms tumour. Of those, 205 lncRNAs, 26 miRNAs and 143 mRNAs were included in the ceRNA regulatory network. The results of Gene Ontology (GO) analysis revealed that the differentially expressed genes (DEGs) were mainly enriched in terms related to response to mechanical stimuli, transcription factor complexes, and transcription factor activity (related to RNA polymerase II proximal promoter sequence-specific DNA binding). The results of the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that the DEGs were mainly enriched in pathways related to the cell cycle. The survival analysis results showed that 16 out of the 205 lncRNAs, 1 out of 26 miRNAs and 5 out of 143 mRNAs were associated with overall survival in Wilms tumour patients (P < 0.05).

Conclusions

CeRNA networks play an important role in Wilms tumour. This finding might provide effective, novel insights for further understanding the mechanisms underlying Wilms tumour.

Reproductive Microbiomes in Wild Animal Species: A New Dimension in Conservation Biology

Communities of microbes have coevolved in animal organisms and are found in almost every part of the body. Compositions of those communities (microbiota) as well as their genomes and genes (microbiomes) are critical for functional regulations of the body organ systems-the digestive or 'gut' microbiome being the most described so far. Based on extensive research in humans, microbiomes in the reproductive tract may play a role in reproductive functions and pregnancy. However, in wild animal species, those microbiomes have been poorly studied, and as a result, little is known about their involvement in fertility or parental/offspring health. This emerging research area is highly relevant to conservation biology from captive breeding management to successful reintroduction or maintenance of wild populations. The objective of this chapter is to review current knowledge about reproductive microbiomes in healthy wild animal species. While recognizing the current technical limits of microbial identification in all animal species, we also explore the link between microbial communities (within female or male reproductive systems) and fertility, from conception to birth outcome. In addition, it is critical to understanding how reproductive microbiomes are affected by environmental factors (including captivity, contact with other individuals, or changes in the ecosystem) to optimize conservation efforts. Thus, reproductive microbiomes represent a novel dimension in conservation biology that will likely gain importance in the future.