Efficient and robust RNA-seq process for cultured bacteria and complex community transcriptomes

Advancing Gut Microbiome Research: The Shift from Metagenomics to Multi-Omics and Future Perspectives

The gut microbiome, a dynamic and integral component of human health, has co-evolved with its host, playing essential roles in metabolism, immunity, and disease prevention. Traditional microbiome studies, primarily focused on microbial composition, have provided limited insights into the functional and mechanistic interactions between microbiota and their host. The advent of multi-omics technologies has transformed microbiome research by integrating genomics, transcriptomics, proteomics, and metabolomics, offering a comprehensive, systems-level understanding of microbial ecology and host-microbiome interactions. These advances have propelled innovations in personalized medicine, enabling more precise diagnostics and targeted therapeutic strategies. This review highlights recent breakthroughs in microbiome research, demonstrating how these approaches have elucidated microbial functions and their implications for health and disease. Additionally, it underscores the necessity of standardizing multi-omics methodologies, conducting large-scale cohort studies, and developing novel platforms for mechanistic studies, which are critical steps toward translating microbiome research into clinical applications and advancing precision medicine.

Advances and applications in single-cell and spatial genomics

The applications of single-cell and spatial technologies in recent times have revolutionized the present understanding of cellular states and the cellular heterogeneity inherent in complex biological systems. These advancements offer unprecedented resolution in the examination of the functional genomics of individual cells and their spatial context within tissues. In this review, we have comprehensively discussed the historical development and recent progress in the field of single-cell and spatial genomics. We have reviewed the breakthroughs in single-cell multi-omics technologies, spatial genomics methods, and the computational strategies employed toward the analyses of single-cell atlas data. Furthermore, we have highlighted the advances made in constructing cellular atlases and their clinical applications, particularly in the context of disease. Finally, we have discussed the emerging trends, challenges, and opportunities in this rapidly evolving field.

RAG-seq: NSR-primed and Transposase Tagmentation-mediated Strand-specific Total RNA Sequencing in Single Cells

Single-cell RNA sequencing (scRNA-seq) has transformed our understanding of cellular diversity with unprecedented resolution. However, many current methods are limited in capturing full-length transcripts and discerning strand orientation. Here, we present RAG-seq, an innovative strand-specific total RNA sequencing technique that combines not-so-random (NSR) primers with Tn5 transposase-mediated tagmentation. RAG-seq overcomes previous limitations by delivering comprehensive transcript coverage and maintaining strand orientation, which are essential for accurate quantification of overlapping genes and detection of antisense transcripts. Through optimized reverse transcription with oligo-dT primers, rRNA depletion via Depletion of Abundant Sequences by Hybridization (DASH), and linear amplification, RAG-seq enhances sensitivity and reproducibility, especially for low-input samples and single cells. Application to mouse oocytes and early embryos highlights RAG-seq's superior performance in identifying stage-specific antisense transcripts, shedding light on their regulatory roles during early development. This advancement represents a significant leap in transcriptome analysis within complex biological contexts.

The application status of sequencing technology in global respiratory infectious disease diagnosis

Next-generation sequencing (NGS) has revolutionized clinical microbiology, particularly in diagnosing respiratory infectious diseases and conducting epidemiological investigations. This narrative review summarizes conventional methods for routine respiratory infection diagnosis, including culture, smear microscopy, immunological assays, image techniques as well as polymerase chain reaction(PCR). In contrast to conventional methods, there is a new detection technology, sequencing technology, and here we mainly focus on the next-generation sequencing NGS, especially metagenomic NGS(mNGS). NGS offers significant advantages over traditional methods. Firstly, mNGS eliminates assumptions about pathogens, leading to faster and more accurate results, thus reducing diagnostic time. Secondly, it allows unbiased identification of known and novel pathogens, offering broad-spectrum coverage. Thirdly, mNGS not only identifies pathogens but also characterizes microbiomes, analyzes human host responses, and detects resistance genes and virulence factors. It can complement targeted sequencing for bacterial and fungal classification. Unlike traditional methods affected by antibiotics, mNGS is less influenced due to the extended survival of pathogen DNA in plasma, broadening its applicability. However, barriers to full integration into clinical practice persist, primarily due to cost constraints and limitations in sensitivity and turnaround time. Despite these challenges, ongoing advancements aim to improve cost-effectiveness and efficiency, making NGS a cornerstone technology for global respiratory infection diagnosis.

A cross-systems primer for synthetic microbial communities

The design and use of synthetic communities, or SynComs, is one of the most promising strategies for disentangling the complex interactions within microbial communities, and between these communities and their hosts. Compared to natural communities, these simplified consortia provide the opportunity to study ecological interactions at tractable scales, as well as facilitating reproducibility and fostering interdisciplinary science. However, the effective implementation of the SynCom approach requires several important considerations regarding the development and application of these model systems. There are also emerging ethical considerations when both designing and deploying SynComs in clinical, agricultural or environmental settings. Here we outline current best practices in developing, implementing and evaluating SynComs across different systems, including a focus on important ethical considerations for SynCom research.

Unraveling the intricacies of host-pathogen interaction through single-cell genomics

Single-cell genomics provide researchers with tools to assess host-pathogen interactions at a resolution previously inaccessible. Transcriptome analysis, epigenome analysis, and immune profiling techniques allow for a better comprehension of the heterogeneity underlying both the host response and infectious agents. Here, we highlight technological advancements and data analysis workflows that increase our understanding of host-pathogen interactions at the single-cell level. We review various studies that have used these tools to better understand host-pathogen dynamics in a variety of infectious disease contexts, including viral, bacterial, and parasitic diseases. We conclude by discussing how single-cell genomics can advance our understanding of host-pathogen interactions.

Unraveling RNA contribution to the molecular origins of bacterial surface-enhanced Raman spectroscopy (SERS) signals

Surface-enhanced Raman spectroscopy (SERS) is widely utilized in bacterial analyses, with the dominant SERS peaks attributed to purine metabolites released during sample preparation. Although adenosine triphosphate (ATP) and nucleic acids are potential molecular origins of these metabolites, research on their exact contributions remains limited. This study explored purine metabolite release from E. coli and RNA integrity following various sample preparation methods. Standard water washing generated dominant SERS signals within 10 s, a duration shorter than the anticipated RNA half-lives under starvation. Evaluating RNA integrity indicated that the most abundant ribosomal RNA species remained intact for hours post-washing, whereas messenger RNA and transfer RNA species degraded gradually. This suggests that bacterial SERS signatures observed after the typical washing step could originate from only a small fraction of endogenous purine-containing molecules. In contrast, acid depurination led to degradation of most RNA species, releasing about 40 times more purine derivatives than water washing. Mild heating also instigated the RNA degradation and released more purine derivatives than water washing. Notably, differences were also evident in the dominant SERS signals following these treatments. This work provides insights into SERS-based studies of purine metabolites released by bacteria and future development of methodologies.

CleanUpRNAseq: An R/Bioconductor Package for Detecting and Correcting DNA Contamination in RNA-Seq Data

RNA sequencing (RNA-seq) has become a standard method for profiling gene expression, yet genomic DNA (gDNA) contamination carried over to the sequencing library poses a significant challenge to data integrity. Detecting and correcting this contamination is vital for accurate downstream analyses. Particularly, when RNA samples are scarce and invaluable, it becomes essential not only to identify but also to correct gDNA contamination to maximize the data's utility. However, existing tools capable of correcting gDNA contamination are limited and lack thorough evaluation. To fill the gap, we developed CleanUpRNAseq, which offers a comprehensive set of functionalities for identifying and correcting gDNA-contaminated RNA-seq data. Our package offers three correction methods for unstranded RNA-seq data and a dedicated approach for stranded data. Through rigorous validation on published RNA-seq datasets with known levels of gDNA contamination and real-world RNA-seq data, we demonstrate CleanUpRNAseq's efficacy in detecting and correcting detrimental levels of gDNA contamination across diverse library protocols. CleanUpRNAseq thus serves as a valuable tool for post-alignment quality assessment of RNA-seq data and should be integrated into routine workflows for RNA-seq data analysis. Its incorporation into OneStopRNAseq should significantly bolster the accuracy of gene expression quantification and differential expression analysis of RNA-seq data.

Transcriptome analysis: a powerful tool to understand individual microbial behaviors and interactions in ecosystems

Transcriptome analysis is a powerful tool for studying microbial ecology, especially individual microbial functions in an ecosystem and their interactions. With the development of high-throughput sequencing technology, great progress has been made in analytical methods for microbial communities in natural environments. 16S rRNA gene amplicon sequencing (ie microbial community structure analysis) and shotgun metagenome analysis have been widely used to determine the composition and potential metabolic capability of microorganisms in target environments without requiring culture. However, even if the types of microorganisms present and their genes are known, it is difficult to determine what they are doing in an ecosystem. Gene expression analysis (transcriptome analysis; RNA-seq) is a powerful tool to address these issues. The history and basic information of gene expression analysis, as well as examples of studies using this method to analyze microbial ecosystems, are presented.

Legume-rhizobia symbiosis: Translatome analysis

Leguminous plants can establish endosymbiotic relationships with nitrogen-fixing soil rhizobacteria. Bacterial infection and nodule organogenesis are two independent but highly coordinated genetic programs that are active during this interaction. These genetic programs can be regulated along all the stages of gene expression. Most of the studies, for both eukaryotes and prokaryotes, focused on the transcriptional regulation level determining the abundance of mRNAs. However, it has been demonstrated that mRNA levels only sometimes correlate with the abundance or activity of the coded proteins. For this reason, in the past two decades, interest in the role of translational control of gene expression has increased, since the subset of mRNA being actively translated outperforms the information gained only by the transcriptome. In the case of legume-rhizobia interactions, the study of the translatome still needs to be explored further. Therefore, this review aims to discuss the methodologies for analyzing polysome-associated mRNAs at the genome-scale and their contribution to studying translational control to understand the complexity of this symbiotic interaction. Moreover, the Dual RNA-seq approach is discussed for its relevance in the context of a symbiotic nodule, where intricate multi-species gene expression networks occur.

An improved bacterial mRNA enrichment strategy in dual RNA sequencing to unveil the dynamics of plant-bacterial interactions

BACKGROUND

Dual RNA sequencing is a powerful tool that enables a comprehensive understanding of the molecular dynamics underlying plant-microbe interactions. RNA sequencing (RNA-seq) poses technical hurdles in the transcriptional analysis of plant-bacterial interactions, especially in bacterial transcriptomics, owing to the presence of abundant ribosomal RNA (rRNA), which potentially limits the coverage of essential transcripts. Therefore, to achieve cost-effective and comprehensive sequencing of the bacterial transcriptome, it is imperative to devise efficient methods for eliminating rRNA and enhancing the proportion of bacterial mRNA. In this study, we modified a strand-specific dual RNA-seq method with the goal of enriching the proportion of bacterial mRNA in the bacteria-infected plant samples. The enriched method involved the sequential separation of plant mRNA by poly A selection and rRNA removal for bacterial mRNA enrichment followed by strand specific RNA-seq library preparation steps. We assessed the efficiency of the enriched method in comparison to the conventional method by employing various plant-bacterial interactions, including both host and non-host resistance interactions with pathogenic bacteria, as well as an interaction with a beneficial rhizosphere associated bacteria using pepper and tomato plants respectively.

RESULTS

In all cases of plant-bacterial interactions examined, an increase in mapping efficiency was observed with the enriched method although it produced a lower read count. Especially in the compatible interaction with Xanthmonas campestris pv. Vesicatoria race 3 (Xcv3), the enriched method enhanced the mapping ratio of Xcv3-infected pepper samples to its own genome (15.09%; 1.45-fold increase) and the CDS (8.92%; 1.49-fold increase). The enriched method consistently displayed a greater number of differentially expressed genes (DEGs) than the conventional RNA-seq method at all fold change threshold levels investigated, notably during the early stages of Xcv3 infection in peppers. The Gene Ontology (GO) enrichment analysis revealed that the DEGs were predominantly enriched in proteolysis, kinase, serine type endopeptidase and heme binding activities.

CONCLUSION

The enriched method demonstrated in this study will serve as a suitable alternative to the existing RNA-seq method to enrich bacterial mRNA and provide novel insights into the intricate transcriptomic alterations within the plant-bacterial interplay.

Characterization of non-invasive oropharyngeal samples and nucleic acid isolation for molecular diagnostics

Molecular diagnostics is an increasingly important clinical tool, especially in routine sampling. We evaluated two non-invasive methods (oral swabs and mouthwashes) for sampling nucleic acids from the oral/pharyngeal area. We created a workflow from sample collection (n = 59) to RT-qPCR based analysis. The samples were further characterized in terms of their cellular composition as well as the purity, degradation and microbial content of the derived DNA/RNA. We determined the optimal housekeeping genes applicable for these types of samples. The cellular composition indicated that mouthwashes contained more immune cells and bacteria. Even though the protocol was not specifically optimized to extract bacterial RNA it was possible to derive microbial RNA, from both sampling methods. Optimizing the protocol allowed us to generate stable quantities of DNA/RNA. DNA/RNA purity parameters were not significantly different between the two sampling methods. Even though integrity analysis demonstrated a high level of degradation of RNA, corresponding parameters confirmed their sequencing potential. RT-qPCR analysis determined TATA-Box Binding Protein as the most favorable housekeeping gene. In summary, we have developed a robust method suitable for multiple downstream diagnostic techniques. This protocol can be used as a foundation for further research endeavors focusing on developing molecular diagnostics for the oropharyngeal cavity.

Targeted rRNA depletion enables efficient mRNA sequencing in diverse bacterial species and complex co-cultures

IMPORTANCE

Microbes present one of the most diverse sources of biochemistry in nature, and mRNA sequencing provides a comprehensive view of this biological activity by quantitatively measuring microbial transcriptomes. However, efficient mRNA capture for sequencing presents significant challenges in prokaryotes as mRNAs are not poly-adenylated and typically make up less than 5% of total RNA compared with rRNAs that exceed 80%. Recently developed methods for sequencing bacterial mRNA typically rely on depleting rRNA by tiling large probe sets against rRNAs; however, such approaches are expensive, time-consuming, and challenging to scale to varied bacterial species and complex microbial communities. Therefore, we developed EMBR-seq+, a method that requires fewer than 10 short oligonucleotides per rRNA to achieve up to 99% rRNA depletion in diverse bacterial species. Finally, EMBR-seq+ resulted in a deeper view of the transcriptome, enabling systematic quantification of how microbial interactions result in altering the transcriptional state of bacteria within co-cultures.

The potential role of the microbiota in prostate cancer pathogenesis and treatment

The human body hosts a complex and dynamic population of trillions of microorganisms - the microbiota - which influences the body in homeostasis and disease, including cancer. Several epidemiological studies have associated specific urinary and gut microbial species with increased risk of prostate cancer; however, causal mechanistic data remain elusive. Studies have associated bacterial generation of genotoxins with the occurrence of TMPRSS2-ERG gene fusions, a common, early oncogenic event during prostate carcinogenesis. A subsequent study demonstrated the role of the gut microbiota in prostate cancer endocrine resistance, which occurs, at least partially, through the generation of androgenic steroids fuelling oncogenic signalling via the androgen receptor. These studies present mechanistic evidence of how the host microbiota might be implicated in prostate carcinogenesis and tumour progression. Importantly, these findings also reveal potential avenues for the detection and treatment of prostate cancer through the profiling and modulation of the host microbiota. The latter could involve approaches such as the use of faecal microbiota transplantation, prebiotics, probiotics, postbiotics or antibiotics, which can be used independently or combined with existing treatments to reverse therapeutic resistance and improve clinical outcomes in patients with prostate cancer.

Rational probe design for efficient rRNA depletion and improved metatranscriptomic analysis of human microbiomes

The microbiota that colonize the human gut and other tissues are dynamic, varying both in composition and functional state between individuals and over time. Gene expression measurements can provide insights into microbiome composition and function. However, efficient and unbiased removal of microbial ribosomal RNA (rRNA) presents a barrier to acquiring metatranscriptomic data. Here we describe a probe set that achieves efficient enzymatic rRNA removal of complex human-associated microbial communities. We demonstrate that the custom probe set can be further refined through an iterative design process to efficiently deplete rRNA from a range of human microbiome samples. Using synthetic nucleic acid spike-ins, we show that the rRNA depletion process does not introduce substantial quantitative error in gene expression profiles. Successful rRNA depletion allows for efficient characterization of taxonomic and functional profiles, including during the development of the human gut microbiome. The pan-human microbiome enzymatic rRNA depletion probes described here provide a powerful tool for studying the transcriptional dynamics and function of the human microbiome.

Current concepts, advances, and challenges in deciphering the human microbiota with metatranscriptomics

Metatranscriptomics refers to the analysis of the collective microbial transcriptome of a sample. Its increased utilization for the characterization of human-associated microbial communities has enabled the discovery of many disease-state related microbial activities. Here, we review the principles of metatranscriptomics-based analysis of human-associated microbial samples. We describe strengths and weaknesses of popular sample preparation, sequencing, and bioinformatics approaches and summarize strategies for their use. We then discuss how human-associated microbial communities have recently been examined and how their characterization may change. We conclude that metatranscriptomics insights into human microbiotas under health and disease have not only expanded our knowledge on human health, but also opened avenues for rational antimicrobial drug use and disease management.

Droplet-based high-throughput single microbe RNA sequencing by smRandom-seq

Bacteria colonize almost all parts of the human body and can differ significantly. However, the population level transcriptomics measurements can only describe the average bacteria population behaviors, ignoring the heterogeneity among bacteria. Here, we report a droplet-based high-throughput single-microbe RNA-seq assay (smRandom-seq), using random primers for in situ cDNA generation, droplets for single-microbe barcoding, and CRISPR-based rRNA depletion for mRNA enrichment. smRandom-seq showed a high species specificity (99%), a minor doublet rate (1.6%), a reduced rRNA percentage (32%), and a sensitive gene detection (a median of ~1000 genes per single E. coli). Furthermore, smRandom-seq successfully captured transcriptome changes of thousands of individual E. coli and discovered a few antibiotic resistant subpopulations displaying distinct gene expression patterns of SOS response and metabolic pathways in E. coli population upon antibiotic stress. smRandom-seq provides a high-throughput single-microbe transcriptome profiling tool that will facilitate future discoveries in microbial resistance, persistence, microbe-host interaction, and microbiome research.

ietanti-flamatoria or DAIN: A Crohn's disease management strategy tailored for Puerto Ricans

Diet has been increasingly shown to be of therapeutic benefit for patients with inflammatory bowel diseases (IBD), especially Crohn's disease (CD). Yet dietary guidelines are nonexistent. Moreover, diets tailored to Puerto Ricans with IBD living on the island, have not been developed and tested. The rising prevalence of IBD in Puerto Rico warrants exploring the use of diet as part of the treatment strategies for these patients [1]. Here, we describe the study design of "Dieta Anti-Inflamatoria" or DAIN, a parallel two-arm randomized pilot trial aiming at testing the efficacy of IBD-Anti-inflammatory diet (IBD-AID) adapted for adults with CD living in Puerto Rico (clinical trial registration number: NCT05627128). We tailored the IBD-AID to the local cuisine preferences and food availability by creating and adapting recipes consistent with the IBD-AID principles [2,3]. In focus groups with a Community Research Advisory Panel and one-on-one consultations with implementation experts, we identified several aspects of the intervention to adapt before the implementation. The objectives of the stakeholder/expert-informed adaptation were to improve feasibility and compliance while developing the culturally tailored dietary intervention. DAIN was designed for adults living in Puerto Rico with CD and geared to be affordable, appropriate, and acceptable for patients with mild-to-moderate CD. The significance of this work is the validation of culturally appropriate nutritional guidelines to help manage CD symptoms. DAIN provides a blueprint for a comprehensive nutritional program that can be adapted to regional preferences and local food availability allowing wider implementation of diet as an adjunct treatment in diverse clinical settings.

Optimization of RNA extraction for bacterial whole transcriptome studies of low-biomass samples

We developed a procedure for extracting maximal amounts of high-quality RNA from low-biomass producing (autotrophic) bacteria for experiments where sample volume is limited. Large amounts of high-quality RNA for downstream analyses cannot be obtained using larger quantities of culture volume. The performance of standard commercial silica-column based kit protocols and these procedures amended by ultrasonication or enzymatic lysis were assessed. The ammonium-oxidizing Nitrosomonas europaea and nitrite-oxidizing Nitrobacter winogradskyi were used as model organisms for optimization of the RNA isolation protocol. Enzymatic lysis through lysozyme digestion generated high-quality, high-yield RNA samples. Subsequent RNA-seq analysis resulted in qualitative data for both strains. The RNA extraction procedure is suitable for experiments with volume and/or biomass limitations, e.g., as encountered during space flight experiments. Furthermore, it will also result in higher RNA yields for whole transcriptome experiments where sample volume and/or biomass was increased to compensate the low-biomass characteristic of autotrophs.

•

Lysis through ultrasonication results in a high RNA yield but in a low RNA quality

•

Enzymatic lysis using lysozyme provides high-yield, high-quality RNA samples

•

RNA from low-biomass bacterial samples is suitable for downstream RNA-seq

Arginine Catabolism and Polyamine Biosynthesis Pathway Disparities Within Francisella tularensis Subpopulations

Francisella tularensis is a highly infectious zoonotic pathogen with as few as 10 organisms causing tularemia, a disease that is fatal if untreated. Although F. tularensis subspecies tularensis (type A) and subspecies holarctica (type B) share over 99.5% average nucleotide identity, notable differences exist in genomic organization and pathogenicity. The type A clade has been further divided into subtypes A.I and A.II, with A.I strains being recognized as some of the most virulent bacterial pathogens known. In this study, we report on major disparities that exist between the F. tularensis subpopulations in arginine catabolism and subsequent polyamine biosynthesis. The genes involved in these pathways include the speHEA and aguAB operons, along with metK. In the hypervirulent F. tularensis A.I clade, such as the A.I prototype strain SCHU S4, these genes were found to be intact and highly transcribed. In contrast, both subtype A.II and type B strains have a truncated speA gene, while the type B clade also has a disrupted aguA and truncated aguB. Ablation of the chromosomal speE gene that encodes a spermidine synthase reduced subtype A.I SCHU S4 growth rate, whereas the growth rate of type B LVS was enhanced. These results demonstrate that spermine synthase SpeE promotes faster replication in the F. tularensis A.I clade, whereas type B strains do not rely on this enzyme for in vitro fitness. Our ongoing studies on amino acid and polyamine flux within hypervirulent A.I strains should provide a better understanding of the factors that contribute to F. tularensis pathogenicity.

Storage media and RNA extraction approaches substantially influence the recovery and integrity of livestock fecal microbial RNA

Background

There is growing interest in understanding gut microbiome dynamics, to increase the sustainability of livestock production systems and to better understand the dynamics that regulate antibiotic resistance genes (i.e., the resistome). High-throughput sequencing of RNA transcripts (RNA-seq) from microbial communities (metatranscriptome) allows an unprecedented opportunity to analyze the functional and taxonomical dynamics of the expressed microbiome and emerges as a highly informative approach. However, the isolation and preservation of high-quality RNA from livestock fecal samples remains highly challenging. This study aimed to determine the impact of the various sample storage and RNA extraction strategies on the recovery and integrity of microbial RNA extracted from selected livestock (chicken and pig) fecal samples.

Methods

Fecal samples from pigs and chicken were collected from conventional slaughterhouses. Two different storage buffers were used at two different storage temperatures. The extraction of total RNA was done using four different commercially available kits and RNA integrity/quality and concentration were measured using a Bioanalyzer 2100 system with RNA 6000 Nano kit (Agilent, Santa Clara, CA, USA). In addition, RT-qPCR was used to assess bacterial RNA quality and the level of host RNA contamination.

Results

The quantity and quality of RNA differed by sample type (i.e., either pig or chicken) and most significantly by the extraction kit, with differences in the extraction method resulting in the least variability in pig feces samples and the most variability in chicken feces. Considering a tradeoff between the RNA yield and the RNA integrity and at the same time minimizing the amount of host RNA in the sample, a combination of storing the fecal samples in RNALater at either 4 °C (for 24 h) or -80 °C (up to 2 weeks) with extraction with PM kit (RNEasy Power Microbiome Kit) had the best performance for both chicken and pig samples.

Conclusion

Our findings provided a further emphasis on using a consistent methodology for sample storage, duration as well as a compatible RNA extraction approach. This is crucial as the impact of these technical steps can be potentially large compared with the real biological variability to be explained in microbiome and resistome studies.

Comparative Analysis of rRNA Removal Methods for RNA-Seq Differential Expression in Halophilic Archaea

Despite intense recent research interest in archaea, the scientific community has experienced a bottleneck in the study of genome-scale gene expression experiments by RNA-seq due to the lack of commercial and specifically designed rRNA depletion kits. The high rRNA:mRNA ratio (80–90%: ~10%) in prokaryotes hampers global transcriptomic analysis. Insufficient ribodepletion results in low sequence coverage of mRNA, and therefore, requires a substantially higher number of replicate samples and/or sequencing reads to achieve statistically reliable conclusions regarding the significance of differential gene expression between case and control samples. Here, we show that after the discontinuation of the previous version of RiboZero (Illumina, San Diego, CA, USA) that was useful in partially or completely depleting rRNA from archaea, archaeal transcriptomics studies have experienced a slowdown. To overcome this limitation, here, we analyze the efficiency for four different hybridization-based kits from three different commercial suppliers, each with two sets of sequence-specific probes to remove rRNA from four different species of halophilic archaea. We conclude that the key for transcriptomic success with the currently available tools is the probe-specificity for the rRNA sequence hybridization. With this paper, we provide insights into the archaeal community for selecting certain reagents and strategies over others depending on the archaeal species of interest. These methods yield improved RNA-seq sensitivity and enhanced detection of low abundance transcripts.

Comparison of rRNA depletion methods for efficient bacterial mRNA sequencing

Current methods of high-throughput RNA sequencing of prokaryotes, including transcriptome analysis or ribosomal profiling, need deep sequencing to achieve sufficient numbers of effective reads (e.g., mapping to mRNA) in order to also find weakly expressed genetic elements. The fraction of high-quality reads mapping to coding RNAs (i.e., mRNA) is mainly influenced by the large content of rRNA and, to a lesser extent, tRNA in total RNA. Thus, depletion of rRNA increases coverage and thus sequencing costs. RiboZero, a depletion kit based on probe hybridisation and rRNA-removal was found to be most efficient in the past, but it was discontinued in 2018. To facilitate comparability with previous experiments and to help choose adequate replacements, we compare three commercially available rRNA depletion kits also based on hybridization and magnetic beads, i.e., riboPOOLs, RiboMinus and MICROBExpress, with the former RiboZero. Additionally, we constructed biotinylated probes for magnetic bead capture and rRNA depletion in this study. Based on E. coli, we found similar efficiencies in rRNA depletion for riboPOOLs and the self-made depletion method; both comparable to the former RiboZero, followed by RiboMinus, succeeded by MICROBExpress. Further, our in-house protocol allows customized species-specific rRNA or even tRNA depletion or depletion of other RNA targets. Both, the self-made biotinylated probes and riboPOOLs, were most successful in reducing the rRNA content and thereby increasing sequencing depth concerning mRNA reads. Additionally, the number of reads matching to weakly expressed genes are increased. In conclusion, the self-made specific biotinylated probes and riboPOOLs are an adequate replacement for the former RiboZero. Both are very efficient in depleting rRNAs, increasing mRNA reads and thus sequencing efficiency.

Efficient and specific DNA oligonucleotide rRNA probe-based rRNA removal in Talaromyces marneffei

Emerging evidence showed that lncRNAs play important roles in a wide range of biological processes of fungi such as Saccharomyces cerevisiae. However, systemic identification of lncRNAs in non-model fungi is a challenging task as the efficiency of rRNA removal has been proved to be affected by mismatches of universal rRNA-targeting probes of commercial kits, which forces deeper sequencing depth and increases costs. Here, we developed a low-cost and simple rRNA depletion method (rProbe) that could efficiently remove more than 99% rRNA in both yeast and mycelium samples of Talaromyces marneffei. The efficiency and robustness of rProbe were demonstrated to outperform the Illumina Ribo-Zero kit. Using rProbe RNA-seq, we identified 115 differentially expressed lncRNAs and constructed lncRNA-mRNA co-expression network related to dimorphic switch of T. marneffei. Our rRNA removal method has the potential to be a useful tool to explore non-coding transcriptomes of non-model fungi by adjusting rRNA probe sequences species specifically.

Recent progress of circular RNAs in different types of human cancer: Technological landscape, clinical opportunities and challenges (Review)

Circular RNAs (circRNAs) are a novel class of endogenous non‑coding RNAs that have been recently regarded as functionally active. CircRNAs are remarkably stable and known to possess several biological functions such as microRNA sponging, regulating transcription and splicing and occasionally acting as polypeptide‑producing templates. CircRNAs show tissue‑specific expression and have been reported to be associated with the progression of several types of malignancies. Given the recent progress in genome sequencing and bioinformatics techniques, a rapid increment in the biological role of circRNAs has been observed. Concurrently, the patent search from different patent databases shows that the patent number of circRNA is increasing very quickly. These phenomena reveal a rapid development of the technological landscape. In the present review, the recent progress on circRNAs in various kinds of cancer has been investigated and their function as biomarkers or therapeutic targets and their technological landscape have been appreciated. A new insight into circRNAs structure and functional capabilities in cancer has been reviewed. Continually increasing knowledge on their critical role during cancer progression is projecting them as biomarkers or therapeutic targets for various kinds of cancer. Thus, recent updates on the functional role of circRNAs in terms of the technological landscape, clinical opportunities (biomarkers and therapeutic targets), and challenges in cancer have been illustrated.

The Small RNA AmiL Regulates Quorum Sensing-Mediated Virulence in Pseudomonas aeruginosa PAO1

Pseudomonas aeruginosa is an opportunistic and nosocomial pathogen of humans with hundreds of its virulence factors regulated by quorum sensing (QS) system. Small noncoding RNAs (sRNAs) are also key regulators of bacterial virulence. However, the QS regulatory sRNAs (Qrrs) that have been characterized in P. aeruginosa are still largely unknown. Here, sRNA AmiL (PA3366.1) in the amiEBCRS operon of PAO1 was identified as a novel Qrr by transcriptome sequencing (RNA-Seq). The expression of AmiL was negatively regulated by the las or rhl system, of which RhlR probably inhibited its transcription. AmiL deletion mutant and overexpressing strains were constructed in PAO1. Broad phenotypic changes were found, including reduced pyocyanin synthesis, elastase activity, biofilm formation, hemolytic activity, and cytotoxicity, as well as increased rhamnolipid production and swarming motility. AmiL appears to be a new regulator that influences diverse QS-mediated virulence. Furthermore, AmiL directly targeted PhzC, a key member of pyocyanin synthesis. AmiL also negatively regulated lasI expression in the early growth of PAO1, but predominantly increased rhlI expression and C4-HSL production in the middle and late stages. Therefore, a novel QS-sRNA signaling cascade of las/rhl (RhlR)-AmiL-PhzC/las/rhl was demonstrated, and it will help to shed new light on the virulence regulatory network of P. aeruginosa PAO1.

IMPORTANCEP. aeruginosa is a common nosocomial pathogen that causes diverse opportunistic infections in humans. The virulence crucial for infection is mainly regulated by QS. Small noncoding RNAs (sRNAs) involved in virulence regulation have also been identified in many bacteria. Recently, there is a growing interest in the new sRNA species, QS regulatory sRNAs (Qrrs). Understanding Qrrs-mediated regulation in P. aeruginosa virulence is therefore important to combat infection. In this study, a previously uncharacterized sRNA AmiL in PAO1 has been identified as a novel Qrr. It has been found to influence diverse QS-mediated virulence factors including pyocyanin, elastase, rhamnolipid, and hemolysin, as well as biofilm formation, swarming motility, and cytotoxicity. Furthermore, PhzC essential for pyocyanin synthesis was a direct target of AmiL. QS gene expression and C4-HSL production were also regulated by AmiL. This study provides insights into the roles of Qrr AmiL in modulating P. aeruginosa virulence.

Non-coding RNAs as liquid biopsy biomarkers in cancer

Although non-coding RNAs have long been considered as non-functional "junk" RNAs, accumulating evidence in the past decade indicates that they play a critical role in pathogenesis of various cancers. In addition to their biological significance, the recognition that their expression levels are frequently dysregulated in multiple cancers have fueled the interest for exploiting their clinical potential as cancer biomarkers. In particular, microRNAs (miRNAs), a subclass of small non-coding RNAs that epigenetically modulate gene-transcription, have become one of the most well-studied substrates for the development of liquid biopsy biomarkers for cancer patients. The emergence of high-throughput sequencing technologies has enabled comprehensive molecular characterisation of various non-coding RNA expression profiles in multiple cancers. Furthermore, technological advances for quantifying lowly expressed RNAs in the circulation have facilitated robust identification of previously unrecognised and undetectable biomarkers in cancer patients. Here we summarise the latest progress on the utilisation of non-coding RNAs as non-invasive cancer biomarkers. We evaluated the suitability of multiple non-coding RNA types as blood-based cancer biomarkers and examined the impact of recent technological breakthroughs on the development of non-invasive molecular biomarkers in cancer.

Target-enriched sequencing enables accurate identification of bloodstream infections in whole blood

Bloodstream infections are within the top ten causes of death globally, with a mortality rate of up to 70%. Gold standard blood culture testing is time-consuming, resulting in delayed, but accurate, treatment. Molecular methods, such as RT-qPCR, have limited targets in one run. We present a new Ampliseq detection system (ADS) combining target amplification and next-generation sequencing for accurate identification of bacteria, fungi, and antimicrobial resistance determinants directly from blood samples. In this study, we included removal of human genomic DNA during nucleic acid extraction, optimized the target sequence set and drug resistance genes, performed antimicrobial resistance profiling of clinical isolates, and evaluated mock specimens and clinical samples by ADS. ADS successfully identified pathogens at the species-level in 36 h, from nucleic acid extraction to results. Besides pathogen identification, ADS can also present drug resistance profiles. ADS enabled detection of all bacteria and accurate identification of 47 pathogens. In 20 spiked samples and 8 clinical specimens, ADS detected at least 92.81% of reads mapped to pathogens. ADS also showed consistency with the three culture-negative samples, and correctly identified pathogens in four of five culture-positive clinical blood specimens. This Ampliseq-based technology promises broad coverage and accurate pathogen identification, helping clinicians to accurately diagnose and treat bloodstream infections.

Transcriptomic analysis of the Myxococcus xanthus FruA regulon, and comparative developmental transcriptomic analysis of two fruiting body forming species, Myxococcus xanthus and Myxococcus stipitatus

BACKGROUND

The Myxococcales are well known for their predatory and developmental social processes, and for the molecular complexity of regulation of these processes. Many species within this order have unusually large genomes compared to other bacteria, and their genomes have many genes that are unique to one specific sequenced species or strain. Here, we describe RNAseq based transcriptome analysis of the FruA regulon of Myxococcus xanthus and a comparative RNAseq analysis of two Myxococcus species, M. xanthus and Myxococcus stipitatus, as they respond to starvation and begin forming fruiting bodies.

RESULTS

We show that both species have large numbers of genes that are developmentally regulated, with over half the genome showing statistically significant changes in expression during development in each species. We also included a non-fruiting mutant of M. xanthus that is missing the transcriptional regulator FruA to identify the direct and indirect FruA regulon and to identify transcriptional changes that are specific to fruiting and not just the starvation response. We then identified Interpro gene ontologies and COG annotations that are significantly up- or down-regulated during development in each species. Our analyses support previous data for M. xanthus showing developmental upregulation of signal transduction genes, and downregulation of genes related to cell-cycle, translation, metabolism, and in some cases, DNA replication. Gene expression in M. stipitatus follows similar trends. Although not all specific genes show similar regulation patterns in both species, many critical developmental genes in M. xanthus have conserved expression patterns in M. stipitatus, and some groups of otherwise unstudied orthologous genes share expression patterns.

CONCLUSIONS

By identifying the FruA regulon and identifying genes that are similarly and uniquely regulated in two different species, this work provides a more complete picture of transcription during Myxococcus development. We also provide an R script to allow other scientists to mine our data for genes whose expression patterns match a user-selected gene of interest.

Comparative analysis of RNA enrichment methods for preparation of Cryptococcus neoformans RNA sequencing libraries

RNA sequencing (RNA-Seq) experiments focused on gene expression involve removal of ribosomal RNA (rRNA) because it is the major RNA constituent of cells. This process, called RNA enrichment, is done primarily to reduce cost: without rRNA removal, deeper sequencing must be performed to compensate for the sequencing reads wasted on rRNA. The ideal RNA enrichment method removes all rRNA without affecting other RNA in the sample. We tested the performance of three RNA enrichment methods on RNA isolated from Cryptococcus neoformans, a fungal pathogen of humans. We find that the RNase H depletion method is more efficient in depleting rRNA and more specific in recapitulating non-rRNA levels present in unenriched controls than the commonly-used Poly(A) isolation method. The RNase H depletion method is also more effective than the Ribo-Zero depletion method as measured by rRNA depletion efficiency and recapitulation of protein-coding RNA levels present in unenriched controls, while the Ribo-Zero depletion method more closely recapitulates annotated non-coding RNA (ncRNA) levels. Finally, we leverage these data to accurately map the C. neoformans mitochondrial rRNA genes, and also demonstrate that RNA-Seq data generated with the RNase H and Ribo-Zero depletion methods can be used to explore novel C. neoformans long non-coding RNA genes.

Mining the Microbiome and Microbiota-Derived Molecules in Inflammatory Bowel Disease

The intestinal microbiota is a complex community that consists of an ecosystem with a dynamic interplay between bacteria, fungi, archaea, and viruses. Recent advances in model systems have revealed that the gut microbiome is critical for maintaining homeostasis through metabolic digestive function, immune regulation, and intestinal barrier integrity. Taxonomic shifts in the intestinal microbiota are strongly correlated with a multitude of human diseases, including inflammatory bowel disease (IBD). However, many of these studies have been descriptive, and thus the understanding of the cause and effect relationship often remains unclear. Using non-human experimental model systems such as gnotobiotic mice, probiotic mono-colonization, or prebiotic supplementation, researchers have defined numerous species-level functions of the intestinal microbiota that have produced therapeutic candidates for IBD. Despite these advances, the molecular mechanisms responsible for the function of much of the microbiota and the interplay with host cellular processes remain areas of tremendous research potential. In particular, future research will need to unlock the functional molecular units of the microbiota in order to utilize this untapped resource of bioactive molecules for therapy. This review will highlight the advances and remaining challenges of microbiota-based functional studies and therapeutic discovery, specifically in IBD. One of the limiting factors for reviewing this topic is the nascent development of this area with information on some drug candidates still under early commercial development. We will also highlight the current and evolving strategies, including in the biotech industry, used for the discovery of microbiota-derived bioactive molecules in health and disease.

RiboRid: A low cost, advanced, and ultra-efficient method to remove ribosomal RNA for bacterial transcriptomics

RNA sequencing techniques have enabled the systematic elucidation of gene expression (RNA-Seq), transcription start sites (differential RNA-Seq), transcript 3′ ends (Term-Seq), and post-transcriptional processes (ribosome profiling). The main challenge of transcriptomic studies is to remove ribosomal RNAs (rRNAs), which comprise more than 90% of the total RNA in a cell. Here, we report a low-cost and robust bacterial rRNA depletion method, RiboRid, based on the enzymatic degradation of rRNA by thermostable RNase H. This method implemented experimental considerations to minimize nonspecific degradation of mRNA and is capable of depleting pre-rRNAs that often comprise a large portion of RNA, even after rRNA depletion. We demonstrated the highly efficient removal of rRNA up to a removal efficiency of 99.99% for various transcriptome studies, including RNA-Seq, Term-Seq, and ribosome profiling, with a cost of approximately $10 per sample. This method is expected to be a robust method for large-scale high-throughput bacterial transcriptomic studies.

Removal of ribosomal RNAs, a major constituent (over 90%) of cellular RNA is a critical experimental step for transcriptomic studies that deal with messenger RNAs. In this manuscript, we describe a robust method to subtract ribosomal RNA from various RNA samples. The method is based on the enzymatic degradation of target RNA by short complementary DNA and RNA:DNA duplex specific nuclease. The method comprises carefully designed experimental procedures to minimize experimental bias and unwanted removal of messenger RNAs. We validate the method on various types of transcriptomic studies for seven diverse bacterial species. This method successfully removed ribosomal RNA with over 99% of efficiency and it was comparable to commercial systems even for degraded RNA samples at a fraction of a cost.

Comparison of transcriptional profiles of Treponema pallidum during experimental infection of rabbits and in vitro culture: Highly similar, yet different

Treponema pallidum ssp. pallidum, the causative agent of syphilis, can now be cultured continuously in vitro utilizing a tissue culture system, and the multiplication rates are similar to those obtained in experimental infection of rabbits. In this study, the RNA transcript profiles of the T. pallidum Nichols during in vitro culture and rabbit infection were compared to examine whether gene expression patterns differed in these two environments. To this end, RNA preparations were converted to cDNA and subjected to RNA-seq using high throughput Illumina sequencing; reverse transcriptase quantitative PCR was also performed on selected genes for validation of results. The transcript profiles in the in vivo and in vitro environments were remarkably similar, exhibiting a high degree of concordance overall. However, transcript levels of 94 genes (9%) out of the 1,063 predicted genes in the T. pallidum genome were significantly different during rabbit infection versus in vitro culture, varying by up to 8-fold in the two environments. Genes that exhibited significantly higher transcript levels during rabbit infection included those encoding multiple ribosomal proteins, several prominent membrane proteins, glycolysis-associated enzymes, replication initiator DnaA, rubredoxin, thioredoxin, two putative regulatory proteins, and proteins associated with solute transport. In vitro cultured T. pallidum had higher transcript levels of DNA repair proteins, cofactor synthesis enzymes, and several hypothetical proteins. The overall concordance of the transcript profiles may indicate that these environments are highly similar in terms of their effects on T. pallidum physiology and growth, and may also reflect a relatively low level of transcriptional regulation in this reduced genome organism.

Metatranscriptomics for the Human Microbiome and Microbial Community Functional Profiling

Shotgun metatranscriptomics (MTX) is an increasingly practical way to survey microbial community gene function and regulation at scale. This review begins by summarizing the motivations for community transcriptomics and the history of the field. We then explore the principles, best practices, and challenges of contemporary MTX workflows: beginning with laboratory methods for isolation and sequencing of community RNA, followed by informatics methods for quantifying RNA features, and finally statistical methods for detecting differential expression in a community context. In thesecond half of the review, we survey important biological findings from the MTX literature, drawing examples from the human microbiome, other (nonhuman) host-associated microbiomes, and the environment. Across these examples, MTX methods prove invaluable for probing microbe-microbe and host-microbe interactions, the dynamics of energy harvest and chemical cycling, and responses to environmental stresses. We conclude with a review of open challenges in the MTX field, including making assays and analyses more robust, accessible, and adaptable to new technologies; deciphering roles for millions of uncharacterized microbial transcripts; and solving applied problems such as biomarker discovery and development of microbial therapeutics.

Simple mapping-based quantification of a mock microbial community using total RNA-seq data

Most microbes in the natural environment are difficult to cultivate. Thus, culture-independent analysis for microbial community structure is important for the understanding of its ecological functions. An immense ribosomal RNA sequence collection is available from phylogenetic research on organisms in all domains. These sequences are available for use in genetic research. However, the amplicon-seq process using PCR requires the construction of a sequence library. Construction can introduce bias into quantitative analyses, and each domain of species needs its own primer set. Total RNA sequencing has the advantage of analyzing an entire microbial community, including bacteria, archea, and eukaryote, at once. Such analysis yields large amounts of ribosomal RNA sequences that can be used for analysis without PCR bias. Evaluation using total RNA-seq for quantitative analysis of microbial communities and comparison with amplicon-seq is still rare. In the present study, we developed a mapping-based total RNA-seq analysis to obtain quantitative information on microbial community structure and compared our results with ordinary amplicon-seq methods. We read total RNA sequences from a commercially available mock community (ATCC MSA-2003) and divided reads into small subunit ribosomal RNA (ssrRNA) origin reads and others, such as mRNA origin reads. We then mapped ssrRNA origin reads on annotated assembled contigs and obtained quantitative results under several analysis strategies. Removal of low complexity sequences, sorting ssrRNA with paired-in mode, and performing homology-based taxonomical assignments (BLAST+ or vsearch) showed superior outcomes to other strategies. Results with this approach showed a median relative abundance among ten mock community members of ~10%; ordinary amplicon-seq showed a much lower percentage. Thus, total RNA-seq can be a powerful tool for analyzing microbial community structure and is not limited to analyzing gene expression profiling of microbiomes.

Statistical approaches for differential expression analysis in metatranscriptomics

Motivation

Metatranscriptomics (MTX) has become an increasingly practical way to profile the functional activity of microbial communities in situ. However, MTX remains underutilized due to experimental and computational limitations. The latter are complicated by non-independent changes in both RNA transcript levels and their underlying genomic DNA copies (as microbes simultaneously change their overall abundance in the population and regulate individual transcripts), genetic plasticity (as whole loci are frequently gained and lost in microbial lineages) and measurement compositionality and zero-inflation. Here, we present a systematic evaluation of and recommendations for differential expression (DE) analysis in MTX.

Results

We designed and assessed six statistical models for DE discovery in MTX that incorporate different combinations of DNA and RNA normalization and assumptions about the underlying changes of gene copies or species abundance within communities. We evaluated these models on multiple simulated and real multi-omic datasets. Models adjusting transcripts relative to their encoding gene copies as a covariate were significantly more accurate in identifying DE from MTX in both simulated and real datasets. Moreover, we show that when paired DNA measurements (metagenomic data) are not available, models normalizing MTX measurements within-species while also adjusting for total-species RNA balance sensitivity, specificity and interpretability of DE detection, as does filtering likely technical zeros. The efficiency and accuracy of these models pave the way for more effective MTX-based DE discovery in microbial communities.

Availability and implementation

The analysis code and synthetic datasets used in this evaluation are available online at http://huttenhower.sph.harvard.edu/mtx2021.

Supplementary information

Supplementary data are available at Bioinformatics online.

Strain engraftment competition and functional augmentation in a multi-donor fecal microbiota transplantation trial for obesity

BACKGROUND

Donor selection is an important factor influencing the engraftment and efficacy of fecal microbiota transplantation (FMT) for complex conditions associated with microbial dysbiosis. However, the degree, variation, and stability of strain engraftment have not yet been assessed in the context of multiple donors.

METHODS

We conducted a double-blinded randomized control trial of FMT in 87 adolescents with obesity. Participants were randomized to receive multi-donor FMT (capsules containing the fecal microbiota of four sex-matched lean donors) or placebo (saline capsules). Following a bowel cleanse, participants ingested a total of 28 capsules over two consecutive days. Capsules from individual donors and participant stool samples collected at baseline, 6, 12, and 26 weeks post-treatment were analyzed by shotgun metagenomic sequencing allowing us to track bacterial strain engraftment and its functional implications on recipients' gut microbiomes.

RESULTS

Multi-donor FMT sustainably altered the structure and the function of the gut microbiome. In what was effectively a microbiome competition experiment, we discovered that two donor microbiomes (one female, one male) dominated strain engraftment and were characterized by high microbial diversity and a high Prevotella to Bacteroides (P/B) ratio. Engrafted strains led to enterotype-level shifts in community composition and provided genes that altered the metabolic potential of the community. Despite our attempts to standardize FMT dose and origin, FMT recipients varied widely in their engraftment of donor strains.

CONCLUSION

Our study provides evidence for the existence of FMT super-donors whose microbiomes are highly effective at engrafting in the recipient gut. Dominant engrafting male and female donor microbiomes harbored diverse microbial species and genes and were characterized by a high P/B ratio. Yet, the high variability of strain engraftment among FMT recipients suggests the host environment also plays a critical role in mediating FMT receptivity.

TRIAL REGISTRATION

The Gut Bugs trial was registered with the Australian New Zealand Clinical Trials Registry ( ACTRN12615001351505 ).

TRIAL PROTOCOL

The trial protocol is available at https://bmjopen.bmj.com/content/9/4/e026174 . Video Abstract.

The Microbiome in Obstructive Sleep Apnea

Recent evidence has highlighted important associations between obstructive sleep apnea and the microbiome. Although the intricacies of the pathophysiologic mechanisms are not well understood, available evidence suggests a bidirectional relationship between OSA and microbiota composition. Sleep fragmentation, intermittent hypoxia, and intermittent hypercapnia all play significant roles in altering the microbiome, and initial evidence has shown that alterations of the microbiota affect sleep patterns. Animal model evidence strongly supports the idea that the microbiome mediates disease states associated with OSA including hypertension, atherosclerosis, and obesity. The majority of evidence focuses on changes in the gut microbiome, which may result from OSA as well as contribute to sleep pattern changes, OSA-related CVD, and obesity. Meanwhile, a developing body of work suggests changes in the upper airway microbiome may be associated with OSA and periodontitis-related oral cavity microbiome changes may have significance in OSA-related CVD. Lastly, while evidence is limited, several studies suggest there may be a role for treatment of OSA and OSA-related comorbidities through alteration of the microbiome with probiotics, prebiotics, and microbiota transplantation. These early animal and human studies begin to characterize the interrelationships of the microbiome and OSA and may lead to new avenues for treatment.

Factors Associated With the Microbiome in Moderate-Late Preterm Babies: A Cohort Study From the DIAMOND Randomized Controlled Trial

The gut microbiota of preterm infants is affected by perinatal factors and, in turn, may impact upon infant health. In this study, we collected fecal samples at Day-10 (D10) and 4-months corrected-age (4M) from 227 moderate–late preterm (MLPT) babies enrolled in a randomized controlled trial of nutritional management. A total of 320 samples underwent 16S amplicon sequencing, and shotgun metagenomic sequencing was performed on 94 samples from the 4M time point. The microbiome of babies whose families lived in lower socioeconomic status (SES) areas exhibited a significantly higher microbial alpha diversity at D10 (Wilcoxon test, p = 0.021), greater abundance of Bifidobacterium (linear model, q = 0.020) at D10 and Megasphaera (q = 0.031) at 4M. Hospital of birth explained 5.2% of the observed variance in 4M samples (PERMANOVA, p = 0.038), with Staphylococcus aureus more abundant in fecal samples from babies born in Middlemore hospital (linear model, q = 0.016). Maternal antibiotic (Wilcoxon test, p = 0.013) and probiotic (p = 0.04) usage within the four-week period before sample collection was associated with a reduction in the alpha diversity of D10 samples. Infant probiotic intake explained 2.1% (PERMANOVA, p = 0.021) of the variance in the D10 microbial profile with increased Lactobacillus (linear model, q = 1.1 × 10⁻¹⁰) levels. At 4M, the microbiome of infants who were breastmilk fed had reduced alpha diversity when compared to non-breastmilk fed infants (Wilcoxon test, p < 0.05). Although causality cannot be inferred within our study, we conclude that in MLPT babies, maternal socioeconomic factors, as well as the perinatal medical environment and nutrition impact on the development of the newborn microbiome.

The Multiomics Analyses of Fecal Matrix and Its Significance to Coeliac Disease Gut Profiling

Gastrointestinal (GIT) diseases have risen globally in recent years, and early detection of the host’s gut microbiota, typically through fecal material, has become a crucial component for rapid diagnosis of such diseases. Human fecal material is a complex substance composed of undigested macromolecules and particles, and the processing of such matter is a challenge due to the unstable nature of its products and the complexity of the matrix. The identification of these products can be used as an indication for present and future diseases; however, many researchers focus on one variable or marker looking for specific biomarkers of disease. Therefore, the combination of genomics, transcriptomics, proteomics and metabonomics can give a detailed and complete insight into the gut environment. The proper sample collection, sample preparation and accurate analytical methods play a crucial role in generating precise microbial data and hypotheses in gut microbiome research, as well as multivariate data analysis in determining the gut microbiome functionality in regard to diseases. This review summarizes fecal sample protocols involved in profiling coeliac disease.

Cross-species RNA-seq for deciphering host-microbe interactions

The human body is constantly exposed to microorganisms, which entails manifold interactions between human cells and diverse commensal or pathogenic bacteria. The cellular states of the interacting cells are decisive for the outcome of these encounters such as whether bacterial virulence programmes and host defence or tolerance mechanisms are induced. This Review summarizes how next-generation RNA sequencing (RNA-seq) has become a primary technology to study host-microbe interactions with high resolution, improving our understanding of the physiological consequences and the mechanisms at play. We illustrate how the discriminatory power and sensitivity of RNA-seq helps to dissect increasingly complex cellular interactions in time and space down to the single-cell level. We also outline how future transcriptomics may answer currently open questions in host-microbe interactions and inform treatment schemes for microbial disorders.

The diversity of Shine-Dalgarno sequences sheds light on the evolution of translation initiation

Shine-Dalgarno (SD) sequences, the core element of prokaryotic ribosome-binding sites, facilitate mRNA translation by base-pair interaction with the anti-SD (aSD) sequence of 16S rRNA. In contrast to this paradigm, an inspection of thousands of prokaryotic species unravels tremendous SD sequence diversity both within and between genomes, whereas aSD sequences remain largely static. The pattern has led many to suggest unidentified mechanisms for translation initiation. Here we review known translation-initiation pathways in prokaryotes. Moreover, we seek to understand the cause and consequence of SD diversity through surveying recent advances in biochemistry, genomics, and high-throughput genetics. These findings collectively show: (1) SD:aSD base pairing is beneficial but nonessential to translation initiation. (2) The 5' untranslated region of mRNA evolves dynamically and correlates with organismal phylogeny and ecological niches. (3) Ribosomes have evolved distinct usage of translation-initiation pathways in different species. We propose a model portraying the SD diversity shaped by optimization of gene expression, adaptation to environments and growth demands, and the species-specific prerequisite of ribosomes to initiate translation. The model highlights the coevolution of ribosomes and mRNA features, leading to functional customization of the translation apparatus in each organism.

Metatranscriptomic analyses of the oral microbiome

Although the composition of the oral human microbiome is now well studied, regulation of genes within oral microbial communities remains mostly uncharacterized. Current concepts of periodontal disease and caries highlight the importance of oral biofilms and their role as etiological agents of those diseases. Currently, there is increased interest in exploring and characterizing changes in the composition and gene-expression profiles of oral microbial communities. These efforts aim to identify changes in functional activities that could explain the transition from health to disease and the reason for the chronicity of those infections. It is now clear that the functions of distinct species within the subgingival microbiota are intimately intertwined with the rest of the microbial community. This point highlights the relevance of examining the expression profile of specific species within the subgingival microbiota in the case of periodontal disease or caries lesions, in the context of the other members of the biofilm in vivo. Metatranscriptomic analysis of the oral community is the starting point for identifying environmental signals that modulate the shift in metabolism of the community from commensal to dysbiotic. These studies give a snapshot of the expression patterns of microbial communities and also allow us to determine triggers to diseases. For example, in the case of caries, studies have unveiled a potential new pathway of sugar metabolism, namely the use of sorbitol as an additional source of carbon by Streptococcus mutans; and in the case of periodontal disease, high levels of extracellular potassium could be a signal of disease. Longitudinal studies are needed to identify the real markers of the initial stages of caries and periodontal disease. More information on the gene-expression profiles of the host, along with the patterns from the microbiome, will lead to a clearer understanding of the modulation of health and disease. This review presents a summary of these initial studies, which have opened the door to a new understanding of the dynamics of the oral community during the dysbiotic process in the oral cavity.

Efficient and cost-effective bacterial mRNA sequencing from low input samples through ribosomal RNA depletion

BACKGROUND

RNA sequencing is a powerful approach to quantify the genome-wide distribution of mRNA molecules in a population to gain deeper understanding of cellular functions and phenotypes. However, unlike eukaryotic cells, mRNA sequencing of bacterial samples is more challenging due to the absence of a poly-A tail that typically enables efficient capture and enrichment of mRNA from the abundant rRNA molecules in a cell. Moreover, bacterial cells frequently contain 100-fold lower quantities of RNA compared to mammalian cells, which further complicates mRNA sequencing from non-cultivable and non-model bacterial species. To overcome these limitations, we report EMBR-seq (Enrichment of mRNA by Blocked rRNA), a method that efficiently depletes 5S, 16S and 23S rRNA using blocking primers to prevent their amplification.

RESULTS

EMBR-seq results in 90% of the sequenced RNA molecules from an E. coli culture deriving from mRNA. We demonstrate that this increased efficiency provides a deeper view of the transcriptome without introducing technical amplification-induced biases. Moreover, compared to recent methods that employ a large array of oligonucleotides to deplete rRNA, EMBR-seq uses a single or a few oligonucleotides per rRNA, thereby making this new technology significantly more cost-effective, especially when applied to varied bacterial species. Finally, compared to existing commercial kits for bacterial rRNA depletion, we show that EMBR-seq can be used to successfully quantify the transcriptome from more than 500-fold lower starting total RNA.

CONCLUSIONS

EMBR-seq provides an efficient and cost-effective approach to quantify global gene expression profiles from low input bacterial samples.

The Potential of Circular RNAs as Cancer Biomarkers

Circular RNA (circRNA) is a covalently closed RNA structure that has several proposed functions related to cancer development. Recently, cancer-specific and tissue-specific circRNAs have been identified by high-throughput sequencing and are curated in publicly available databases. CircRNAs have features that are ideal properties of biomarkers, including conservation, abundance, and stability in plasma, saliva, and urine. Many circRNAs with predictive and prognostic significance in cancer have been described, and functional mechanisms for some circRNAs have been suggested. CircRNA also has great potential as a noninvasive biomarker for early cancer detection, although further investigation is necessary before clinical application is feasible.See all articles in this CEBP Focus section, "NCI Early Detection Research Network: Making Cancer Detection Possible."

Improved bacterial RNA-seq by Cas9-based depletion of ribosomal RNA reads

A major challenge for RNA-seq analysis of gene expression is to achieve sufficient coverage of informative nonribosomal transcripts. In eukaryotic samples, this is typically achieved by selective oligo(dT)-priming of messenger RNAs to exclude ribosomal RNA (rRNA) during cDNA synthesis. However, this strategy is not compatible with prokaryotes in which functional transcripts are generally not polyadenylated. To overcome this, we adopted DASH (depletion of abundant sequences by hybridization), initially developed for eukaryotic cells, to improve both the sensitivity and depth of bacterial RNA-seq. DASH uses the Cas9 nuclease to remove unwanted cDNA sequences prior to library amplification. We report the design, evaluation, and optimization of DASH experiments for standard bacterial short-read sequencing approaches, including software for automated guide RNA (gRNA) design for Cas9-mediated cleavage in bacterial rDNA sequences. Using these gRNA pools, we effectively removed rRNA reads (56%-86%) in RNA-seq libraries from two different model bacteria, the Gram-negative pathogen Salmonella enterica and the anaerobic gut commensal Bacteroides thetaiotaomicron DASH works robustly, even with subnanogram amounts of input RNA. Its efficiency, high sensitivity, ease of implementation, and low cost (∼$5 per sample) render DASH an attractive alternative to rRNA removal protocols, in particular for material-constrained studies where conventional ribodepletion techniques fail.

A period of 10 weeks of increased protein consumption does not alter faecal microbiota or volatile metabolites in healthy older men: a randomised controlled trial

Diet has a major influence on the composition and metabolic output of the gut microbiome. Higher-protein diets are often recommended for older consumers; however, the effect of high-protein diets on the gut microbiota and faecal volatile organic compounds (VOC) of elderly participants is unknown. The purpose of the study was to establish if the faecal microbiota composition and VOC in older men are different after a diet containing the recommended dietary intake (RDA) of protein compared with a diet containing twice the RDA (2RDA). Healthy males (74⋅2 (sd 3⋅6) years; n 28) were randomised to consume the RDA of protein (0⋅8 g protein/kg body weight per d) or 2RDA, for 10 weeks. Dietary protein was provided via whole foods rather than supplementation or fortification. The diets were matched for dietary fibre from fruit and vegetables. Faecal samples were collected pre- and post-intervention for microbiota profiling by 16S ribosomal RNA amplicon sequencing and VOC analysis by head space/solid-phase microextraction/GC-MS. After correcting for multiple comparisons, no significant differences in the abundance of faecal microbiota or VOC associated with protein fermentation were evident between the RDA and 2RDA diets. Therefore, in the present study, a twofold difference in dietary protein intake did not alter gut microbiota or VOC indicative of altered protein fermentation.

Compositional and functional characterisation of biomass-degrading microbial communities in guts of plant fibre- and soil-feeding higher termites

BACKGROUND

Termites are among the most successful insect lineages on the globe and are responsible for providing numerous ecosystem services. They mainly feed on wood and other plant material at different stages of humification. Lignocellulose is often a principal component of such plant diet, and termites largely rely on their symbiotic microbiota and associated enzymes to decompose their food efficiently. While lower termites and their gut flagellates were given larger scientific attention in the past, the gut lignocellulolytic bacteria of higher termites remain less explored. Therefore, in this study, we investigated the structure and function of gut prokaryotic microbiomes from 11 higher termite genera representative of Syntermitinae, Apicotermitinae, Termitidae and Nasutitermitinae subfamilies, broadly grouped into plant fibre- and soil-feeding termite categories.

RESULTS

Despite the different compositional structures of the studied termite gut microbiomes, reflecting well the diet and host lineage, we observed a surprisingly high functional congruency between gut metatranscriptomes from both feeding groups. The abundance of transcripts encoding for carbohydrate active enzymes as well as expression and diversity profiles of assigned glycoside hydrolase families were also similar between plant fibre- and soil-feeding termites. Yet, dietary imprints highlighted subtle metabolic differences specific to each feeding category. Roughly, 0.18% of de novo re-constructed gene transcripts were shared between the different termite gut microbiomes, making each termite gut a unique reservoir of genes encoding for potentially industrially applicable enzymes, e.g. relevant to biomass degradation. Taken together, we demonstrated the functional equivalence in microbial populations across different termite hosts.

CONCLUSIONS

Our results provide valuable insight into the bacterial component of the termite gut system and significantly expand the inventory of termite prokaryotic genes participating in the deconstruction of plant biomass. Video Abstract.

Differences in Compositions of Gut Bacterial Populations and Bacteriophages in 5-11 Year-Olds Born Preterm Compared to Full Term

Preterm infants are exposed to major perinatal, post-natal, and early infancy events that could impact on the gut microbiome. These events include infection, steroid and antibiotic exposure, parenteral nutrition, necrotizing enterocolitis, and stress. Studies have shown that there are differences in the gut microbiome during the early months of life in preterm infants. We hypothesized that differences in the gut microbial composition and metabolites in children born very preterm persist into mid-childhood. Participants were healthy prepubertal children aged 5-11 years who were born very preterm (≤32 weeks of gestation; n = 51) or at term (37-41 weeks; n = 50). We recorded the gestational age, birth weight, mode of feeding, mode of birth, age, sex, and the current height and weight of our cohort. We performed a multi'omics [i.e., 16S rRNA amplicon and shotgun metagenomic sequencing, SPME-GCMS (solid-phase microextraction followed by gas chromatography-mass spectrometry)] analysis to investigate the structure and function of the fecal microbiome (as a proxy of the gut microbiota) in our cross-sectional cohort. Children born very preterm were younger (7.8 vs. 8.3 years; p = 0.034), shorter [height-standard deviation score (SDS) 0.31 vs. 0.92; p = 0.0006) and leaner [BMI (body mass index) SDS -0.20 vs. 0.29; p < 0.0001] than the term group. Children born very preterm had higher fecal calprotectin levels, decreased fecal phage richness, lower plasma arginine, lower fecal branched-chain amino acids and higher fecal volatile (i.e., 3-methyl-butanoic acid, butyrolactone, butanoic acid and pentanoic acid) profiles. The bacterial microbiomes did not differ between preterm and term groups. We speculate that the observed very preterm-specific changes were established in early infancy and may impact on the capacity of the very preterm children to respond to environmental changes.