Validation of two ribosomal RNA removal methods for microbial metatranscriptomics

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.

Diversity and composition of active and total bacteria in rhizospheric soil in response to continuous cropping years of Panax notoginseng

The synchronous research and analysis of total and active soil microbial communities can provide insight into how these communities are impacted by continuous cropping years and pathogen infection. The diversity of total and active bacteria in rhizospheric soil of 2-year-old and 3-year-old healthy and diseased Panax notoginseng can comprehensively reveal the bacterial response characteristics in continuous cropping practice. The results showed that 4916 operational taxonomic units (OTUs) were found in the rhizospheric soil bacterial community of P. notoginseng at the DNA level, but only 2773 OTUs were found at the RNA level. The rhizospheric environment had significant effects on the active and bacterial communities, as indicated by the number of OTUs, Shannon, Chao1, Faith's phylogenetic diversity (Faith's PD), and Simpson's diversity indexes. The DNA level can better show the difference in diversity level before and after infection with root rot. The bacterial Chao1 and Faith's PD diversity indexes of 2-year-old root rot-diseased P. notoginseng rhizospheric soil (D2) were higher than that of healthy plants, while the bacterial Shannon diversity index of 3-year-old root rot-diseased P. notoginseng rhizospheric soil (D3) was the lowest in the total bacteria. Principal coordinate analysis (PCoA) illustrated that the total bacterial species composition changed markedly after root rot disease. There were significant differences in the composition of active bacterial species between the 2-year and 3-year rhizospheres. In conclusion, the total and active edaphic rhizospheric bacterial communities could provide important opportunities to understand the responses of bacteria to continuous cropping of P. notoginseng.

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.

Dual RNA-seq in filarial nematodes and Wolbachia endosymbionts using RNase H based ribosomal RNA depletion

Lymphatic filariasis is caused by parasitic nematodes and is a leading cause of disability worldwide. Many filarial worms contain the bacterium Wolbachia as an obligate endosymbiont. RNA sequencing is a common technique used to study their molecular relationships and to identify potential drug targets against the nematode and bacteria. Ribosomal RNA (rRNA) is the most abundant RNA species, accounting for 80-90% of the RNA in a sample. To reduce sequencing costs, it is necessary to remove ribosomal reads through poly-A enrichment or ribosomal depletion. Bacterial RNA does not contain a poly-A tail, making it difficult to sequence both the nematode and Wolbachia from the same library preparation using standard poly-A selection. Ribosomal depletion can utilize species-specific oligonucleotide probes to remove rRNA through pull-down or degradation methods. While species-specific probes are commercially available for many commonly studied model organisms, there are currently limited depletion options for filarial parasites. Here, we performed total RNA sequencing from Brugia malayi containing the Wolbachia symbiont (wBm) and designed ssDNA depletion probes against their rRNA sequences. We compared the total RNA library to poly-A enriched, Terminator 5'-Phosphate-Dependent Exonuclease treated, NEBNext Human/Bacteria rRNA depleted and our custom nematode probe depleted libraries. The custom nematode depletion library had the lowest percentage of ribosomal reads across all methods, with a 300-fold decrease in rRNA when compared to the total RNA library. The nematode depletion libraries also contained the highest percentage of Wolbachia mRNA reads, resulting in a 16-1,000-fold increase in bacterial reads compared to the other enrichment and depletion methods. Finally, we found that the Brugia malayi depletion probes can remove rRNA from the filarial worm Dirofilaria immitis and the majority of rRNA from the more distantly related free living nematode Caenorhabditis elegans. These custom filarial probes will allow for future dual RNA-seq experiments between nematodes and their bacterial symbionts from a single sequencing library.

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.

A clinically validated human saliva metatranscriptomic test for global systems biology studies

The authors report here the development of a high-throughput, automated, inexpensive and clinically validated saliva metatranscriptome test that requires less than 100 μl of saliva. RNA is preserved at the time of sample collection, allowing for ambient-temperature transportation and storage for up to 28 days. Critically, the RNA preservative is also able to inactivate pathogenic microorganisms, rendering the samples noninfectious and allowing for safe and easy shipping. Given the unique set of convenience, low cost, safety and technical performance, this saliva metatranscriptomic test can be integrated into longitudinal, global-scale systems biology studies that will lead to an accelerated development of precision medicine, diagnostic and therapeutic tools.

Single-cell RNA sequencing in the context of neuropathic pain: progress, challenges, and prospects

Neuropathic pain, characterized by persistent or intermittent spontaneous pain as well as some unpleasant abnormal sensations, is one of the most prevalent health problems in the world. Ectopic nerve activity, central and peripheral nociceptive sensitization and many other potential mechanisms may participate in neuropathic pain. The complexity and ambiguity of neuropathic pain mechanisms result in difficulties in pain management, and existing treatment plans provide less-than-satisfactory relief. In recent years, single-cell RNA sequencing (scRNA-seq) has been increasingly applied and has become a powerful means for biological researchers to explore the complexity of neurobiology. This technique can be used to perform unbiased, high-throughput and high-resolution transcriptional analyses of neuropathic pain-associated cells, improving the understanding of neuropathic pain mechanisms and enabling individualized pain management. To date, scRNA-seq has been preliminarily used in neuropathic pain research for applications such as compiling a dorsal root ganglion atlas, identifying new cell types and discovering gene regulatory networks associated with neuropathic pain. Although scRNA-seq is a relatively new technique in the neuropathic pain field, there have been several studies based on animal models. However, because of the various differences between animals and humans, more attention should be given to translational medicine research. With the aid of scRNA-seq, researchers can further explore the mechanism of neuropathic pain to improve the clinical understanding of the diagnosis, treatment and management of neuropathic pain.

Pathogen detection and characterization from throat swabs using unbiased metatranscriptomic analyses

OBJECTIVES

Infectious diseases are common but are not easily or readily diagnosed with current methodologies. This problem is further exacerbated by the constant presence of mutated, emerging, and novel pathogens. One of the most common sites of infection by many pathogens is the human throat. However, there is no universal diagnostic test that can distinguish these pathogens. Metatranscriptomic (MT) analysis of the throat represents an important and novel development in infectious disease detection and characterization, because it is able to identify all pathogens using a fully unbiased approach.

METHODS

To test the utility of an MT approach to pathogen detection, throat samples were collected from participants before, during, and after an acute sickness.

RESULTS

Clear sickness-associated shifts in pathogenic microorganisms were detected in the patients. Important insights into microbial functions and antimicrobial resistance genes were obtained.

CONCLUSION

MT analysis of the throat represents an effective method for the unbiased identification and characterization of pathogens. Because MT data include all microorganisms in the sample, this approach should not only allow the identification of pathogens, but provide an understanding of the effects of the resident throat microbiome in the context of human health and disease.

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.

Oceanic Crustal Fluid Single Cell Genomics Complements Metagenomic and Metatranscriptomic Surveys With Orders of Magnitude Less Sample Volume

Fluids circulating through oceanic crust play important roles in global biogeochemical cycling mediated by their microbial inhabitants, but studying these sites is challenged by sampling logistics and low biomass. Borehole observatories installed at the North Pond study site on the western flank of the Mid-Atlantic Ridge have enabled investigation of the microbial biosphere in cold, oxygenated basaltic oceanic crust. Here we test a methodology that applies redox-sensitive fluorescent molecules for flow cytometric sorting of cells for single cell genomic sequencing from small volumes of low biomass (approximately 103 cells ml-1) crustal fluid. We compare the resulting genomic data to a recently published paired metagenomic and metatranscriptomic analysis from the same site. Even with low coverage genome sequencing, sorting cells from less than one milliliter of crustal fluid results in similar interpretation of dominant taxa and functional profiles as compared to 'omics analysis that typically filter orders of magnitude more fluid volume. The diverse community dominated by Gammaproteobacteria, Bacteroidetes, Desulfobacterota, Alphaproteobacteria, and Zetaproteobacteria, had evidence of autotrophy and heterotrophy, a variety of nitrogen and sulfur cycling metabolisms, and motility. Together, results indicate fluorescence activated cell sorting methodology is a powerful addition to the toolbox for the study of low biomass systems or at sites where only small sample volumes are available for analysis.

Gut Microbiome Activity Contributes to Prediction of Individual Variation in Glycemic Response in Adults

Limiting postprandial glycemic response (PPGR) is an important intervention in reducing the risk of chronic metabolic diseases and has been shown to impart significant health benefits in people with elevated levels of blood sugar. In this study, we collected gut microbiome activity data by assessing the metatranscriptome, and we measured the glycemic responses of 550 adults who consumed more than 30,000 meals, collectively, from omnivore or vegetarian/gluten-free diets. We demonstrate that gut microbiome activity, anthropometric factors, and food macronutrients modulate individual variation in glycemic response. We employ two predictive models, including a mixed-effects linear regression model (R = 0.77) and a gradient boosting machine model (Rtrain = 0.80/R2train = 0.64; Rtest = 0.64/R2test = 0.40), which demonstrate variation in PPGR between individuals when ingesting the same foods. All features in the final mixed-effects linear regression model were significant (p < 0.05) except for two features which were retained as suggestive: glutamine production pathways (p = 0.08) and the interaction between tyrosine metabolizers and carbs (p = 0.06). We introduce molecular functions as features in these two models, aggregated from microbial activity data, and show their statistically significant contributions to glycemic control. In summary, we demonstrate for the first time that metatranscriptomic activity of the gut microbiome is correlated with PPGR among adults.

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.

Prospects for multi-omics in the microbial ecology of water engineering

Advances in high-throughput sequencing technologies and bioinformatics approaches over almost the last three decades have substantially increased our ability to explore microorganisms and their functions - including those that have yet to be cultivated in pure isolation. Genome-resolved metagenomic approaches have enabled linking powerful functional predictions to specific taxonomical groups with increasing fidelity. Additionally, related developments in both whole community gene expression surveys and metabolite profiling have permitted for direct surveys of community-scale functions in specific environmental settings. These advances have allowed for a shift in microbiome science away from descriptive studies and towards mechanistic and predictive frameworks for designing and harnessing microbial communities for desired beneficial outcomes. Water engineers, microbiologists, and microbial ecologists studying activated sludge, anaerobic digestion, and drinking water distribution systems have applied various (meta)omics techniques for connecting microbial community dynamics and physiologies to overall process parameters and system performance. However, the rapid pace at which new omics-based approaches are developed can appear daunting to those looking to apply these state-of-the-art practices for the first time. Here, we review how modern genome-resolved metagenomic approaches have been applied to a variety of water engineering applications from lab-scale bioreactors to full-scale systems. We describe integrated omics analysis across engineered water systems and the foundations for pairing these insights with modeling approaches. Lastly, we summarize emerging omics-based technologies that we believe will be powerful tools for water engineering applications. Overall, we provide a framework for microbial ecologists specializing in water engineering to apply cutting-edge omics approaches to their research questions to achieve novel functional insights. Successful adoption of predictive frameworks in engineered water systems could enable more economically and environmentally sustainable bioprocesses as demand for water and energy resources increases.

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.

Next generation sequencing approaches to evaluate water and wastewater quality

The emergence of next generation sequencing (NGS) is revolutionizing the potential to address complex microbiological challenges in the water industry. NGS technologies can provide holistic insight into microbial communities and their functional capacities in water and wastewater systems, thus eliminating the need to develop a new assay for each target organism or gene. However, several barriers have hampered wide-scale adoption of NGS by the water industry, including cost, need for specialized expertise and equipment, challenges with data analysis and interpretation, lack of standardized methods, and the rapid pace of development of new technologies. In this critical review, we provide an overview of the current state of the science of NGS technologies as they apply to water, wastewater, and recycled water. In addition, a systematic literature review was conducted in which we identified over 600 peer-reviewed journal articles on this topic and summarized their contributions to six key areas relevant to the water and wastewater fields: taxonomic classification and pathogen detection, functional and catabolic gene characterization, antimicrobial resistance (AMR) profiling, bacterial toxicity characterization, Cyanobacteria and harmful algal bloom identification, and virus characterization. For each application, we have presented key trends, noteworthy advancements, and proposed future directions. Finally, key needs to advance NGS technologies for broader application in water and wastewater fields are assessed.

Environmental RNA: A Revolution in Ecological Resolution?

Current advancements in environmental RNA (eRNA) exploit its relatively fast turnover rate relative to environmental DNA (eDNA) to assess 'metabolically active' or temporally/spatially recent community diversity. However, this focus significantly underutilizes the trove of potential ecological information encrypted in eRNA. Here, we argue for pushing beyond current species-level eDNA detection capabilities by using eRNA to detect any organisms with unique eRNA profiles, potentially including different life-history stages, sexes, or even specific phenotypes within a species. We also discuss the future of eRNA as a means of assessing the physiological status of organisms and the ecological health of populations and communities, reflecting ecosystem-level conditions. We posit that eRNA has the potential to significantly improve the resolution of organism detection, biological monitoring, and biomonitoring applications in ecology.

Microbial single-cell RNA sequencing by split-pool barcoding

Single-cell RNA sequencing (scRNA-seq) has become an essential tool for characterizing gene expression in eukaryotes, but current methods are incompatible with bacteria. Here, we introduce microSPLiT (microbial split-pool ligation transcriptomics), a high-throughput scRNA-seq method for Gram-negative and Gram-positive bacteria that can resolve heterogeneous transcriptional states. We applied microSPLiT to >25,000 Bacillus subtilis cells sampled at different growth stages, creating an atlas of changes in metabolism and lifestyle. We retrieved detailed gene expression profiles associated with known, but rare, states such as competence and prophage induction and also identified unexpected gene expression states, including the heterogeneous activation of a niche metabolic pathway in a subpopulation of cells. MicroSPLiT paves the way to high-throughput analysis of gene expression in bacterial communities that are otherwise not amenable to single-cell analysis, such as natural microbiota.

Comparing sediment DNA extraction methods for assessing organic enrichment associated with marine aquaculture

Marine sediments contain a high diversity of micro- and macro-organisms which are important in the functioning of biogeochemical cycles. Traditionally, anthropogenic perturbation has been investigated by identifying macro-organism responses along gradients. Environmental DNA (eDNA) analyses have recently been advocated as a rapid and cost-effective approach to measuring ecological impacts and efforts are underway to incorporate eDNA tools into monitoring. Before these methods can replace or complement existing methods, robustness and repeatability of each analytical step has to be demonstrated. One area that requires further investigation is the selection of sediment DNA extraction method. Environmental DNA sediment samples were obtained along a disturbance gradient adjacent to a Chinook (Oncorhynchus tshawytscha) salmon farm in Otanerau Bay, New Zealand. DNA was extracted using four extraction kits (Qiagen DNeasy PowerSoil, Qiagen DNeasy PowerSoil Pro, Qiagen RNeasy PowerSoil Total RNA/DNA extraction/elution and Favorgen FavorPrep Soil DNA Isolation Midi Kit) and three sediment volumes (0.25, 2, and 5 g). Prokaryotic and eukaryotic communities were amplified using primers targeting the 16S and 18S ribosomal RNA genes, respectively, and were sequenced on an Illumina MiSeq. Diversity and community composition estimates were obtained from each extraction kit, as well as their relative performance in established metabarcoding biotic indices. Differences were observed in the quality and quantity of the extracted DNA amongst kits with the two Qiagen DNeasy PowerSoil kits performing best. Significant differences were observed in both prokaryotes and eukaryotes (p < 0.001) richness among kits. A small proportion of amplicon sequence variants (ASVs) were shared amongst the kits (~3%) although these shared ASVs accounted for the majority of sequence reads (prokaryotes: 59.9%, eukaryotes: 67.2%). Differences were observed in the richness and relative abundance of taxonomic classes revealed with each kit. Multivariate analysis showed that there was a significant interaction between "distance" from the farm and "kit" in explaining the composition of the communities, with the distance from the farm being a stronger determinant of community composition. Comparison of the kits against the bacterial and eukaryotic metabarcoding biotic index suggested that all kits showed similar patterns along the environmental gradient. Overall, we advocate for the use of Qiagen DNeasy PowerSoil kits for use when characterizing prokaryotic and eukaryotic eDNA from marine farm sediments. We base this conclusion on the higher DNA quality values and richness achieved with these kits compared to the other kits/amounts investigated in this study. The additional advantage of the PowerSoil Kits is that DNA extractions can be performed using an extractor robot, offering additional standardization and reproducibility of results.

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.

Prokaryotic single-cell RNA sequencing by in situ combinatorial indexing

Despite longstanding appreciation of gene expression heterogeneity in isogenic bacterial populations, affordable and scalable technologies for studying single bacterial cells have been limited. Although single-cell RNA sequencing (scRNA-seq) has revolutionized studies of transcriptional heterogeneity in diverse eukaryotic systems1-13, the application of scRNA-seq to prokaryotes has been hindered by their extremely low mRNA abundance14-16, lack of mRNA polyadenylation and thick cell walls17. Here, we present prokaryotic expression profiling by tagging RNA in situ and sequencing (PETRI-seq)-a low-cost, high-throughput prokaryotic scRNA-seq pipeline that overcomes these technical obstacles. PETRI-seq uses in situ combinatorial indexing11,12,18 to barcode transcripts from tens of thousands of cells in a single experiment. PETRI-seq captures single-cell transcriptomes of Gram-negative and Gram-positive bacteria with high purity and low bias, with median capture rates of more than 200 mRNAs per cell for exponentially growing Escherichia coli. These characteristics enable robust discrimination of cell states corresponding to different phases of growth. When applied to wild-type Staphylococcus aureus, PETRI-seq revealed a rare subpopulation of cells undergoing prophage induction. We anticipate that PETRI-seq will have broad utility in defining single-cell states and their dynamics in complex microbial communities.

Draft genome sequence of scale drop disease virus (SDDV) retrieved from metagenomic investigation of infected barramundi, Lates calcarifer (Bloch, 1790)

Scale drop disease virus (SDDV) is a novel viral pathogen considered to be distributed in farmed barramundi (Lates calcarifer) in South-East Asia. Despite the severity of the disease, only limited genomic information related to SDDV is available. In this study, samples of SDDV-infected fish collected in 2019 were used. The microbiome of brain tissue was investigated using Illumina HiSeq DNA sequencing. Taxonomic analysis showed that SDDV was the main pathogen contained in the affected barramundi. De novo metagenome assembly recovered the SDDV genome, named isolate TH2019, 131 kb in length, and comprised of 135 ORFs. Comparison between this genome and the Singaporean SDDV reference genome revealed that the nucleotide identity within the aligned region was 99.97%. Missense, frameshift, insertion and deletion mutations were identified in 26 ORFs. Deletion of four deduced amino acid sequence in ORF_030L, identical to the SDDV isolate previously identified in Thailand, would be a potential biomarker for future strain classification. Interestingly, the genome of SDDV TH2019 harboured a unique 7,695-bp-long genomic region containing six hypothetical protein-encoded genes. Collectively, this study demonstrated that the SDDV genome can be sequenced directly, although with limited coverage depth, using metagenomic analysis of barramundi sample with severe infection.

Epigenetic mechanisms regulating COVID-19 infection

Coronavirus disease 2019 (COVID-2019) outbreak originating in December 2019 in Wuhan, China has emerged as a global threat to human health. The highly contagious SARS-CoV-2 infection and transmission presents a diversity of human host and increased disease risk with advancing age, highlighting the importance of in-depth understanding of its biological properties. Structural analyses have elucidated hot spots in viral binding domains, mutations, and specific proteins in the host such as the receptor angiotensin-converting enzyme 2 (ACE2) and the transmembrane protease serine 2 (TMPRSS2) to be implicated in cell entry and viral infectivity. Furthermore, epigenetic changes that regulate chromatin structure have shown a major impact in genome stabilization and maintenance of cellular homoeostasis and they have been implicated in the pathophysiology of the virus infection. Epigenetic research has revealed that global DNA methylation along with ACE2 gene methylation and post-translational histone modifications may drive differences in host tissue-, biological age- and sex-biased patterns of viral infection. Moreover, modulation of the host cells epigenetic landscape following infection represents a molecular tool used by viruses to antagonize cellular signalling as well as sensing components that regulate the induction of the host innate immune and antiviral defence programmes in order to enhance viral replication and infection efficiency. In this review, we provide an update of the main research findings at the interface of epigenetics and coronavirus infection. In particular, we highlight the epigenetic factors that interfere with viral replication and infection and may contribute to COVID-19 susceptibility, offering new ways of thinking in respect to host viral response.

A Simple, Cost-Effective, and Robust Method for rRNA Depletion in RNA-Sequencing Studies

The profiling of gene expression by RNA sequencing (RNA-seq) has enabled powerful studies of global transcriptional patterns in all organisms, including bacteria. Because the vast majority of RNA in bacteria is rRNA, it is standard practice to deplete the rRNA from a total RNA sample such that the reads in an RNA-seq experiment derive predominantly from mRNA. One of the most commonly used commercial kits for rRNA depletion, the Ribo-Zero kit from Illumina, was recently discontinued abruptly and for an extended period of time. Here, we report the development of a simple, cost-effective, and robust method for depleting rRNA that can be easily implemented by any lab or facility. We first developed an algorithm for designing biotinylated oligonucleotides that will hybridize tightly and specifically to the 23S, 16S, and 5S rRNAs from any species of interest. Precipitation of these oligonucleotides bound to rRNA by magnetic streptavidin-coated beads then depletes rRNA from a complex, total RNA sample such that ∼75 to 80% of reads in a typical RNA-seq experiment derive from mRNA. Importantly, we demonstrate a high correlation of RNA abundance or fold change measurements in RNA-seq experiments between our method and the Ribo-Zero kit. Complete details on the methodology are provided, including open-source software for designing oligonucleotides optimized for any bacterial species or community of interest.IMPORTANCE The ability to examine global patterns of gene expression in microbes through RNA sequencing has fundamentally transformed microbiology. However, RNA-seq depends critically on the removal of rRNA from total RNA samples. Otherwise, rRNA would comprise upward of 90% of the reads in a typical RNA-seq experiment, limiting the reads coming from mRNA or requiring high total read depth. A commonly used kit for rRNA subtraction from Illumina was recently unavailable for an extended period of time, disrupting routine rRNA depletion. Here, we report the development of a "do-it-yourself" kit for rapid, cost-effective, and robust depletion of rRNA from total RNA. We present an algorithm for designing biotinylated oligonucleotides that will hybridize to the rRNAs from a target set of species. We then demonstrate that the designed oligonucleotides enable sufficient rRNA depletion to produce RNA-seq data with 75 to 80% of reads coming from mRNA. The methodology presented should enable RNA-seq studies on any species or metagenomic sample of interest.

Scalable and cost-effective ribonuclease-based rRNA depletion for transcriptomics

Bacterial RNA sequencing (RNA-seq) is a powerful approach for quantitatively delineating the global transcriptional profiles of microbes in order to gain deeper understanding of their physiology and function. Cost-effective bacterial RNA-seq requires efficient physical removal of ribosomal RNA (rRNA), which otherwise dominates transcriptomic reads. However, current methods to effectively deplete rRNA of diverse non-model bacterial species are lacking. Here, we describe a probe and ribonuclease based strategy for bacterial rRNA removal. We implemented the method using either chemically synthesized oligonucleotides or amplicon-based single-stranded DNA probes and validated the technique on three novel gut microbiota isolates from three distinct phyla. We further showed that different probe sets can be used on closely related species. We provide a detailed methods protocol, probe sets for >5000 common microbes from RefSeq, and an online tool to generate custom probe libraries. This approach lays the groundwork for large-scale and cost-effective bacterial transcriptomics studies.

Zymomonas mobilis metabolism: Novel tools and targets for its rational engineering

Zymomonas mobilis is an α-proteobacterium that interests the biofuel industry due to its perfect ethanol fermentation yields. From its first description as a bacterial isolate in fermented alcoholic beverages to date, Z. mobilis has been rigorously studied in directions basic and applied. The Z. mobilis powerful Entner-Doudoroff glycolytic pathway has been the center of rigorous biochemical studies and, aside from ethanol, it has attracted interest in terms of high-added-value chemical manufacturing. Energetic balances and the effects of respiration have been explored in fundamental directions as also in applications pursuing strain enhancement and the utilization of alternative carbon sources. Metabolic modeling has addressed the optimization of the biochemical circuitry at various conditions of growth and/or substrate utilization; it has been also critical in predicting desirable end-product yields via flux redirection. Lastly, stress tolerance has received particular attention, since it directly determines biocatalytical performance at challenging bioreactor conditions. At a genetic level, advances in the genetic engineering of the organism have brought forth beneficial manipulations in the Z. mobilis gene pool, e.g., knock-outs, knock-ins and gene stacking, aiming to broaden the metabolic repertoire and increase robustness. Recent omic and expressional studies shed light on the genomic content of the most applied strains and reveal landscapes of activity manifested at ambient or reactor-based conditions. Studies such as those reviewed in this work, contribute to the understanding of the biology of Z. mobilis, enable insightful strain development, and pave the way for the transformation of Z. mobilis into a consummate organism for biomass conversion.

Computational search for UV radiation resistance strategies in Deinococcus swuensis isolated from Paramo ecosystems

Ultraviolet radiation (UVR) is widely known as deleterious for many organisms since it can cause damage to biomolecules either directly or indirectly via the formation of reactive oxygen species. The goal of this study was to analyze the capacity of high-mountain Espeletia hartwegiana plant phyllosphere microorganisms to survive UVR and to identify genes related to resistance strategies. A strain of Deinococcus swuensis showed a high survival rate of up to 60% after UVR treatment at 800J/m2 and was used for differential expression analysis using RNA-seq after exposing cells to 400J/m2 of UVR (with >95% survival rate). Differentially expressed genes were identified using the R-Bioconductor package NOISeq and compared with other reported resistance strategies reported for this genus. Genes identified as being overexpressed included transcriptional regulators and genes involved in protection against damage by UVR. Non-coding (nc)RNAs were also differentially expressed, some of which have not been previously implicated. This study characterized the immediate radiation response of D. swuensis and indicates the involvement of ncRNAs in the adaptation to extreme environmental conditions.

A Robust Metatranscriptomic Technology for Population-Scale Studies of Diet, Gut Microbiome, and Human Health

A functional readout of the gut microbiome is necessary to enable precise control of the gut microbiome's functions, which support human health and prevent or minimize a wide range of chronic diseases. Stool metatranscriptomic analysis offers a comprehensive functional view of the gut microbiome, but despite its usefulness, it has rarely been used in clinical studies due to its complexity, cost, and bioinformatic challenges. This method has also received criticism due to potential intrasample variability, rapid changes, and RNA degradation. Here, we describe a robust and automated stool metatranscriptomic method, called Viomega, which was specifically developed for population-scale studies. Viomega includes sample collection, ambient temperature sample preservation, total RNA extraction, physical removal of ribosomal RNAs (rRNAs), preparation of directional Illumina libraries, Illumina sequencing, taxonomic classification based on a database of >110,000 microbial genomes, and quantitative microbial gene expression analysis using a database of ~100 million microbial genes. We applied this method to 10,000 human stool samples and performed several small-scale studies to demonstrate sample stability and consistency. In summary, Viomega is an inexpensive, high-throughput, automated, and accurate sample-to-result stool metatranscriptomic technology platform for large-scale studies and a wide range of applications.

Construction of Genetic Logic Gates Based on the T7 RNA Polymerase Expression System in Rhodococcus opacus PD630

Rhodococcus opacus PD630 (R. opacus) is a nonmodel, Gram-positive bacterium that holds promise as a biological catalyst for the conversion of lignocellulosic biomass to value-added products. In particular, it demonstrates both a high tolerance for and an ability to consume inhibitory lignin-derived aromatics, generates large quantities of lipids, exhibits a relatively rapid growth rate, and has a growing genetic toolbox for engineering. However, the availability of genetic parts for tunable, high-activity gene expression is still limited in R. opacus. Furthermore, genetic logic circuits for sophisticated gene regulation have never been demonstrated in Rhodococcus spp. To address these shortcomings, two inducible T7 RNA polymerase-based expression systems were implemented for the first time in R. opacus and applied to the construction of AND and NAND genetic logic gates. Additionally, three isopropyl β-d-1-thiogalactopyranoside (IPTG)-inducible promoters were created by inserting LacI binding sites into newly characterized constitutive promoters. Furthermore, four novel aromatic sensors for 4-hydroxybenzoic acid, vanillic acid, sodium benzoate, and guaiacol were developed, expanding the gene expression toolbox. Finally, the T7 RNA polymerase platform was combined with a synthetic IPTG-inducible promoter to create an IMPLY logic gate. Overall, this work represents the first demonstration of a heterologous RNA polymerase system and synthetic genetic logic in R. opacus, enabling complex and tunable gene regulation in this promising nonmodel host for bioproduction.

Effect of soil horizon stratigraphy on the microbial ecology of alpine paleosols

There is remarkable potential for research at the interface between the earth sciences and environmental microbiology that may lead to advances in our understanding of the role of bacterial communities in the surface or subsurface environment of our planet. One mainstay of sedimentary classification is the concept of differential soil and/or paleosol horizons being the result of primarily physical and chemical weathering, with relatively little understanding of how microbial communities between these stratified horizons differ, if at all. In this study we evaluate the differences in microbial community taxonomy and biogeochemical functional potential between stratified soil horizons in an alpine paleosol environment using next-generation sequencing (NGS) shotgun sequencing. Paleosols represent a unique environment to study the effect of differences soil horizon environments on the microbial community due to their relative isolation, and the fact that three distinct stratified soil horizons can be identified within the top 30 cm of the soil. This enables us to assess variation in microbial community composition that will be relatively distinct from variation due to distance alone. We test the hypothesis that variation in soil community composition is linked to variation in the physical and chemical parameters that define stratigraphy. Multivariate statistical analysis of sequencing reads from soil horizons across five sampling sites revealed that 1223 microbial genera vary significantly and consistently in abundance across stratified soil horizons at class level. Specifically Ktedonobacter, Bacilli and Betaproteobacteria responded most strongly to soil depth. Alpha diversity showed a positive correlation with soil depth. Beta diversity, however, did not differ significantly between horizons. Genes involved in carbohydrate and nitrogen metabolism were found to be more abundant in Ah horizon samples. Closer inspection of carbohydrate metabolism genes revealed that genes involved in CO₂ fixation, fermentation and saccharide metabolism decreased in abundance with depth while one‑carbon metabolism increased down profile.

Comparative metagenomic and metatranscriptomic analyses reveal the breed effect on the rumen microbiome and its associations with feed efficiency in beef cattle

BACKGROUND

Microorganisms are responsible for fermentation within the rumen and have been reported to contribute to the variation in feed efficiency of cattle. However, to what extent the breed affects the rumen microbiome and its association with host feed efficiency is unknown. Here, rumen microbiomes of beef cattle (n = 48) from three breeds (Angus, Charolais, Kinsella composite hybrid) with high and low feed efficiency were explored using metagenomics and metatranscriptomics, aiming to identify differences between functional potentials and activities of same rumen microbiomes and to evaluate the effects of host breed and feed efficiency on the rumen microbiome.

RESULTS

Rumen metagenomes were more closely clustered together and thus more conserved among individuals than metatranscriptomes, suggesting that inter-individual functional variations at the RNA level were higher than those at the DNA level. However, while mRNA enrichment significantly increased the sequencing depth of mRNA and generated similar functional profiles to total RNA-based metatranscriptomics, it led to biased abundance estimation of several transcripts. We observed divergent rumen microbial composition (metatranscriptomic level) and functional potentials (metagenomic level) among three breeds, but differences in functional activity (metatranscriptomic level) were less apparent. Differential rumen microbial features (e.g., taxa, diversity indices, functional categories, and genes) were detected between cattle with high and low feed efficiency, and most of them were breed-specific.

CONCLUSIONS

Metatranscriptomes represent real-time functional activities of microbiomes and have the potential to better associate rumen microorganisms with host performances compared to metagenomics. As total RNA-based metatranscriptomics seem to avoid potential biases caused by mRNA enrichment and allow simultaneous use of rRNA for generation of compositional profiles, we suggest their use for linking the rumen microbiome with host phenotypes in future studies. However, if exploration of specific lowly expressed genes is desired, mRNA enrichment is recommended as it will enhance the resolution of mRNA. Finally, the differential microbial features observed between efficient and inefficient steers tended to be specific to breeds, suggesting that interactions between host breed genotype and the rumen microbiome contribute to the variations in feed efficiency observed. These breed-associated differences represent an opportunity to engineer specific rumen microbiomes through selective breeding of the hosts.

The Ribosome as a Missing Link in Prebiotic Evolution III: Over-Representation of tRNA- and rRNA-Like Sequences and Plieofunctionality of Ribosome-Related Molecules Argues for the Evolution of Primitive Genomes from Ribosomal RNA Modules

We propose that ribosomal RNA (rRNA) formed the basis of the first cellular genomes, and provide evidence from a review of relevant literature and proteonomic tests. We have proposed previously that the ribosome may represent the vestige of the first self-replicating entity in which rRNAs also functioned as genes that were transcribed into functional messenger RNAs (mRNAs) encoding ribosomal proteins. rRNAs also encoded polymerases to replicate itself and a full complement of the transfer RNAs (tRNAs) required to translate its genes. We explore here a further prediction of our “ribosome-first” theory: the ribosomal genome provided the basis for the first cellular genomes. Modern genomes should therefore contain an unexpectedly large percentage of tRNA- and rRNA-like modules derived from both sense and antisense reading frames, and these should encode non-ribosomal proteins, as well as ribosomal ones with key cell functions. Ribosomal proteins should also have been co-opted by cellular evolution to play extra-ribosomal functions. We review existing literature supporting these predictions. We provide additional, new data demonstrating that rRNA-like sequences occur at significantly higher frequencies than predicted on the basis of mRNA duplications or randomized RNA sequences. These data support our “ribosome-first” theory of cellular evolution.

Emergence of Bias During the Synthesis and Amplification of cDNA for scRNA-seq

The advent of single-cell omics technology has promoted our understanding of the genomic, epigenomic, and transcriptomic heterogeneity in individual cells. Compared to traditional sequencing studies using bulk cells, single-cell transcriptome technology is naturally more dynamic for in depth analysis of genomic variation resulting from cell division and is useful in unraveling the regulatory mechanisms of gene networks in many diseases. However, there are still some limitations of current single-cell RNA sequencing (scRNA-seq) protocols. Biases that arise during the RNA reverse transcription and cDNA pre-amplification steps are the most common problems and play pivotal roles in limiting the quantitative accuracy of scRNA-seq. In this review, we will describe how these biases emerge and impact scRNA-seq protocols. Moreover, we will introduce several current and convenient modified scRNA-seq methods that allow for bias to be decreased and estimated.

RNA profile diversity across arthropoda: guidelines, methodological artifacts, and expected outcomes

High-quality RNA is an important precursor for high-throughput RNA sequencing (RNAseq) and subsequent analyses. However, the primary metric used to assess RNA quality, the RNA Integrity Number (RIN), was developed based on model bacterial and vertebrate organisms. Though the phenomenon is not widely recognized, invertebrate 28S ribosomal RNA (rRNA) is highly prone to a form of denaturation known as gap deletion, in which the subunit collapses into two smaller fragments. In many nonmodel invertebrates, this collapse of the 28S subunit appears as a single band similar in size to the 18S rRNA subunit. This phenomenon is hypothesized to be commonplace among arthropods and is often misinterpreted as a "degraded" rRNA profile. The limited characterization of gap deletion in arthropods, a highly diverse group, as well as other nonmodel invertebrates, often biases RNA quality assessments. To test whether the collapse of 28S is a general pattern or a methodological artifact, we sampled more than half of the major lineages within Arthropoda. We found that the 28S collapse is present in ∼90% of the species sampled. Nevertheless, RNA profiles exhibit considerable diversity with a range of banding patterns. High-throughput RNAseq and subsequent assembly of high-quality transcriptomes from select arthropod species exhibiting collapsed 28S subunits further illustrates the limitations of current RIN proxies in accurately characterizing RNA quality in nonmodel organisms. Furthermore, we show that this form of 28S denaturation, which is often mistaken for true "degradation," can occur at relatively low temperatures.

Elucidating the impact of microbial community biodiversity on pharmaceutical biotransformation during wastewater treatment

In addition to removing organics and other nutrients, the microorganisms in wastewater treatment plants (WWTPs) biotransform many pharmaceuticals present in wastewater. The objective of this study was to examine the relationship between pharmaceutical biotransformation and biodiversity in WWTP bioreactor microbial communities and identify taxa and functional genes that were strongly associated with biotransformation. Dilution-to-extinction of an activated sludge microbial community was performed to establish cultures with a gradient of microbial biodiversity. Batch experiments were performed using the dilution cultures to determine biotransformation extents of several environmentally relevant pharmaceuticals. With this approach, because the communities were all established from the same original community, and using sequencing of the 16S rRNA and metatranscriptome, we identified candidate taxa and genes whose activity and transcript abundances associated with the extent of individual pharmaceutical biotransformation and were lost across the biodiversity gradient. Metabolic genes such as dehydrogenases, amidases and monooxygenases were significantly associated with pharmaceutical biotransformation, and five genera were identified whose activity significantly associated with pharmaceutical biotransformation. Understanding how biotransformation relates to biodiversity will inform the design of biological WWTPs for enhanced removal of chemicals that negatively impact environmental health.

Laboratory Methods in Molecular Epidemiology: Viral Infections

Viruses, which are the most abundant biological entities on the planet, have been regarded as the "dark matter" of biology in the sense that despite their ubiquity and frequent presence in large numbers, their detection and analysis are not always straightforward. The majority of them are very small (falling under the limit of 0.5 μm), and collectively, they are extraordinarily diverse. In fact, the majority of the genetic diversity on the planet is found in the so-called virosphere, or the world of viruses. Furthermore, the most frequent viral agents of disease in humans display an RNA genome, and frequently evolve very fast, due to the fact that most of their polymerases are devoid of proofreading activity. Therefore, their detection, genetic characterization, and epidemiological surveillance are rather challenging. This review (part of the Curated Collection on Advances in Molecular Epidemiology of Infectious Diseases) describes many of the methods that, throughout the last few decades, have been used for viral detection and analysis. Despite the challenge of having to deal with high genetic diversity, the majority of these methods still depend on the amplification of viral genomic sequences, using sequence-specific or sequence-independent approaches, exploring thermal profiles or a single nucleic acid amplification temperature. Furthermore, viral populations, and especially those with RNA genomes, are not usually genetically uniform but encompass swarms of genetically related, though distinct, viral genomes known as viral quasispecies. Therefore, sequence analysis of viral amplicons needs to take this fact into consideration, as it constitutes a potential analytic problem. Possible technical approaches to deal with it are also described here. *This article is part of a curated collection.

Clinical Metagenomic Next-Generation Sequencing for Pathogen Detection

Nearly all infectious agents contain DNA or RNA genomes, making sequencing an attractive approach for pathogen detection. The cost of high-throughput or next-generation sequencing has been reduced by several orders of magnitude since its advent in 2004, and it has emerged as an enabling technological platform for the detection and taxonomic characterization of microorganisms in clinical samples from patients. This review focuses on the application of untargeted metagenomic next-generation sequencing to the clinical diagnosis of infectious diseases, particularly in areas in which conventional diagnostic approaches have limitations. The review covers ( a) next-generation sequencing technologies and common platforms, ( b) next-generation sequencing assay workflows in the clinical microbiology laboratory, ( c) bioinformatics analysis of metagenomic next-generation sequencing data, ( d) validation and use of metagenomic next-generation sequencing for diagnosing infectious diseases, and ( e) significant case reports and studies in this area. Next-generation sequencing is a new technology that has the promise to enhance our ability to diagnose, interrogate, and track infectious diseases.

Using noninvasive metagenomics to characterize viral communities from wildlife

Microbial communities play an important role in organismal and ecosystem health. While high-throughput metabarcoding has revolutionized the study of bacterial communities, generating comparable viral communities has proven elusive, particularly in wildlife samples where the diversity of viruses and limited quantities of viral nucleic acid present distinctive challenges. Metagenomic sequencing is a promising solution for studying viral communities, but the lack of standardized methods currently precludes comparisons across host taxa or localities. Here, we developed an untargeted shotgun metagenomic sequencing protocol to generate comparable viral communities from noninvasively collected faecal and oropharyngeal swabs. Using samples from common vampire bats (Desmodus rotundus), a key species for virus transmission to humans and domestic animals, we tested how different storage media, nucleic acid extraction procedures and enrichment steps affect viral community detection. Based on finding viral contamination in foetal bovine serum, we recommend storing swabs in RNAlater or another nonbiological medium. We recommend extracting nucleic acid directly from swabs rather than from supernatant or pelleted material, which had undetectable levels of viral RNA. Results from a low-input RNA library preparation protocol suggest that ribosomal RNA depletion and light DNase treatment reduce host and bacterial nucleic acid, and improve virus detection. Finally, applying our approach to twelve pooled samples from seven localities in Peru, we showed that detected viral communities saturated at the attained sequencing depth, allowing unbiased comparisons of viral community composition. Future studies using the methods outlined here will elucidate the determinants of viral communities across host species, environments and time.

RNA-Seq employing a novel rRNA depletion strategy reveals a rich repertoire of snoRNAs in Euglena gracilis including box C/D and Ψ-guide RNAs targeting the modification of rRNA extremities

Previous mRNA transcriptome studies of Euglena gracilis have shown that this organism possesses a large and diverse complement of protein coding genes; however, the study of non-coding RNA classes has been limited. The natural extensive fragmentation of the E. gracilis large subunit ribosomal RNA presents additional barriers to the identification of non-coding RNAs as size-selected small RNA libraries will be dominated by rRNA sequences. In this study we have developed a strategy to significantly reduce rRNA amplification prior to RNA-Seq analysis thereby producing a ncRNA library allowing for the identification of many new E. gracilis small RNAs. Library analysis reveals 113 unique new small nucleolar (sno) RNAs and a large collection of snoRNA isoforms, as well as the first significant collection of nuclear tRNAs in this organism. A 3' end AGAUGN consensus motif and conserved structural features can now be defined for E. gracilis pseudouridine guide RNAs. snoRNAs of both classes were identified that target modification of the 3' extremities of rRNAs utilizing predicted base-pairing interactions with internally transcribed spacers (ITS), providing insight into the timing of steps in rRNA maturation. Cumulatively, this represents the most comprehensive analysis of small ncRNAs in Euglena gracilis to date.

Sizing Up the Uncultured Microbial Majority

Predicting the total number of microbial cells on Earth and exploring the full diversity of life are fundamental research concepts that have undergone paradigm shifts in the genomic era. In this issue, Lloyd and colleagues (K.

Predicting the total number of microbial cells on Earth and exploring the full diversity of life are fundamental research concepts that have undergone paradigm shifts in the genomic era. In this issue, Lloyd and colleagues (K. G. Lloyd, A. D. Steen, J. L. Ladau, J. Yin, and L. Crosby, mSystems 3:e00055-18, https://doi.org/10.1128/mSystems.00055-18, 2018) present results that combine these two concepts by estimating the total diversity of all cells from Earth’s environments. Leveraging publicly available amplicon, metagenomic, and metatranscriptomic datasets, they determined that nearly all environments are dominated by uncultured lineages, with the exception of humans and human-associated habitats. They define a new concept: phylogenetically diverse noncultured cells (PDNC). Unlike viable but nonculturable cells (VBNC), PDNC are microorganisms for which traditional isolation techniques may never succeed. Lloyd et al. estimate that the majority of microorganisms in Earth’s ecosystems may be PDNC and conclude that culture-independent methods combined with innovative culturing techniques may be required to understand the ecology and physiology of these abundant and divergent microorganisms.

Quantitative RNA-seq Analysis Unveils Osmotic and Thermal Adaptation Mechanisms Relevant for Ectoine Production in Chromohalobacter salexigens

Quantitative RNA sequencing (RNA-seq) and the complementary phenotypic assays were implemented to investigate the transcriptional responses of Chromohalobacter salexigens to osmotic and heat stress. These conditions trigger the synthesis of ectoine and hydroxyectoine, two compatible solutes of biotechnological interest. Our findings revealed that both stresses make a significant impact on C. salexigens global physiology. Apart from compatible solute metabolism, the most relevant adaptation mechanisms were related to “oxidative- and protein-folding- stress responses,” “modulation of respiratory chain and related components,” and “ion homeostasis.” A general salt-dependent induction of genes related to the metabolism of ectoines, as well as repression of ectoine degradation genes by temperature, was observed. Different oxidative stress response mechanisms, secondary or primary, were induced at low and high salinity, respectively, and repressed by temperature. A higher sensitivity to H₂O₂ was observed at high salinity, regardless of temperature. Low salinity induced genes involved in “protein-folding-stress response,” suggesting disturbance of protein homeostasis. Transcriptional shift of genes encoding three types of respiratory NADH dehydrogenases, ATP synthase, quinone pool, Na⁺/H⁺ antiporters, and sodium-solute symporters, was observed depending on salinity and temperature, suggesting modulation of the components of the respiratory chain and additional systems involved in the generation of H⁺ and/or Na⁺ gradients. Remarkably, the Na⁺ intracellular content remained constant regardless of salinity and temperature. Disturbance of Na⁺- and H⁺-gradients with specific ionophores suggested that both gradients influence ectoine production, but with differences depending on the solute, salinity, and temperature conditions. Flagellum genes were strongly induced by salinity, and further induced by temperature. However, salt-induced cell motility was reduced at high temperature, possibly caused by an alteration of Na⁺ permeability by temperature, as dependence of motility on Na⁺-gradient was observed. The transcriptional induction of genes related to the synthesis and transport of siderophores correlated with a higher siderophore production and intracellular iron content only at low salinity. An excess of iron increased hydroxyectoine accumulation by 20% at high salinity. Conversely, it reduced the intracellular content of ectoines by 50% at high salinity plus high temperature. These findings support the relevance of iron homeostasis for osmoadaptation, thermoadaptation and accumulation of ectoines, in C. salexigens.

A Transcriptomic Approach to Identify Novel Drug Efflux Pumps in Bacteria

The core genomes of most bacterial species include a large number of genes encoding putative efflux pumps. The functional roles of most of these pumps are unknown, however, they are often under tight regulatory control and expressed in response to their substrates. Therefore, one way to identify pumps that function in antimicrobial resistance is to examine the transcriptional responses of efflux pump genes to antimicrobial shock. By conducting complete transcriptomic experiments following antimicrobial shock treatments, it may be possible to identify novel drug efflux pumps encoded in bacterial genomes. In this chapter we describe a complete workflow for conducting transcriptomic analyses by RNA sequencing, to determine transcriptional changes in bacteria responding to antimicrobials.

Burden-driven feedback control of gene expression

Cells use feedback regulation to ensure robust growth despite fluctuating demands for resources and differing environmental conditions. However, the expression of foreign proteins from engineered constructs is an unnatural burden that cells are not adapted for. Here we combined RNA-seq with an in vivo assay to identify the major transcriptional changes that occur in Escherichia coli when inducible synthetic constructs are expressed. We observed that native promoters related to the heat-shock response activated expression rapidly in response to synthetic expression, regardless of the construct. Using these promoters, we built a dCas9-based feedback-regulation system that automatically adjusts the expression of a synthetic construct in response to burden. Cells equipped with this general-use controller maintained their capacity for native gene expression to ensure robust growth and thus outperformed unregulated cells in terms of protein yield in batch production. This engineered feedback is to our knowledge the first example of a universal, burden-based biomolecular control system and is modular, tunable and portable.

Microbial effects of part-stream low-frequency ultrasonic pretreatment on sludge anaerobic digestion as revealed by high-throughput sequencing-based metagenomics and metatranscriptomics

BACKGROUND

Part-stream low-frequency ultrasound (LFUS) was one of the common practices for sludge disintegration in full-scale anaerobic digestion (AD) facilities. However, the effectiveness of part-stream LFUS treatment and its effect on AD microbiome have not been fully elucidated.

METHODS

Here we testified the effectiveness of part-stream LFUS pretreatment by treating only a fraction of feed sludge (23% and 33% total solid of the feed sludge) with 20 Hz LFUS for 70 s. State-of-the-art metagenomic and metatranscriptomic analysis was used to investigate the microbial process underpinning the enhanced AD performance by part-stream LFUS pretreatment.

RESULTS

By pretreating 33% total solid of the feed sludge, methane yield was increased by 36.5%, while the volatile solid reduction ratio remained unchanged. RNA-seq of the microbiome at stable stage showed that the continuous dosage of easy-degradable LFUS-pretreated feed sludge had gradually altered the microbial community by selecting Bacteroidales hydrolyzer with greater metabolic capability to hydrolyze cellulosic biomass without substrate attachment. Meanwhile, Thermotogales with excellent cell mobility for nutrient capturing was highly active within the community. Foremost proportion of the methanogenesis was contributed by the dominant Methanomicrobiales via carbon dioxide reduction. More interestingly, a perceivable proportion of the reverse electron flow of the community was input from Methanoculleus species other than syntrophic acetate-oxidizing bacteria. In addition, metagenomic binning retrieved several interesting novel metagenomic-assembled genomes (MAGs): MAG-bin6 of Alistipes shahii showed exceptional transcriptional activities towards protein degradation and MAG-bin11 of Candidatus Cloacimonetes with active cellulolytic GH74 gene detected.

CONCLUSIONS

In summary, despite the unchanged sludge digestibility, the applied part-stream LFUS pretreatment strategy was robust in adjusting the microbial pathways towards more effective substrate conversion enabled by free-living hydrolyser and beta-oxidation-capable methanogens.

Inhibition of anaerobic digestion processes: Applications of molecular tools

Inhibition of anaerobic digestion (AD) due to perturbation caused by substrate composition and/or operating conditions can significantly reduce performance. Such perturbations could be limited by elucidating microbial community response to inhibitors and devising strategies to increase community resilience. To this end, advanced molecular methods are increasingly being applied to study the AD microbiome, a diverse community of microbial populations with complex interactions. This literature review of AD inhibition studies indicates that inhibitory concentrations are highly variable, likely stemming from differences in community structure or activity profile and previous exposure to inhibitors. More recent molecular methods such as 'omics' tools, substrate mapping, and real-time sequencing are helping to unravel the complexity of AD inhibition by elucidating physiological and ecological significance of key microbial populations. The AD community must strive towards developing predictive abilities to avoid system failure (e.g., real-time tracking of an indicator species) to improve resilience of AD systems.

Modified RNA-seq method for microbial community and diversity analysis using rRNA in different types of environmental samples

RNA-seq-based SSU (small subunit) rRNA (ribosomal RNA) analysis has provided a better understanding of potentially active microbial community within environments. However, for RNA-seq library construction, high quantities of purified RNA are typically required. We propose a modified RNA-seq method for SSU rRNA-based microbial community analysis that depends on the direct ligation of a 5’ adaptor to RNA before reverse-transcription. The method requires only a low-input quantity of RNA (10–100 ng) and does not require a DNA removal step. The method was initially tested on three mock communities synthesized with enriched SSU rRNA of archaeal, bacterial and fungal isolates at different ratios, and was subsequently used for environmental samples of high or low biomass. For high-biomass salt-marsh sediments, enriched SSU rRNA and total nucleic acid-derived RNA-seq datasets revealed highly consistent community compositions for all of the SSU rRNA sequences, and as much as 46.4%-59.5% of 16S rRNA sequences were suitable for OTU (operational taxonomic unit)-based community and diversity analyses with complete coverage of V1-V2 regions. OTU-based community structures for the two datasets were also highly consistent with those determined by all of the 16S rRNA reads. For low-biomass samples, total nucleic acid-derived RNA-seq datasets were analyzed, and highly active bacterial taxa were also identified by the OTU-based method, notably including members of the previously underestimated genus Nitrospira and phylum Acidobacteria in tap water, members of the phylum Actinobacteria on a shower curtain, and members of the phylum Cyanobacteria on leaf surfaces. More than half of the bacterial 16S rRNA sequences covered the complete region of primer 8F, and non-coverage rates as high as 38.7% were obtained for phylum-unclassified sequences, providing many opportunities to identify novel bacterial taxa. This modified RNA-seq method will provide a better snapshot of diverse microbial communities, most notably by OTU-based analysis, even communities with low-biomass samples.

Capturing the genetic makeup of the active microbiome in situ

More than any other technology, nucleic acid sequencing has enabled microbial ecology studies to be complemented with the data volumes necessary to capture the extent of microbial diversity and dynamics in a wide range of environments. In order to truly understand and predict environmental processes, however, the distinction between active, inactive and dead microbial cells is critical. Also, experimental designs need to be sensitive toward varying population complexity and activity, and temporal as well as spatial scales of process rates. There are a number of approaches, including single-cell techniques, which were designed to study in situ microbial activity and that have been successively coupled to nucleic acid sequencing. The exciting new discoveries regarding in situ microbial activity provide evidence that future microbial ecology studies will indispensably rely on techniques that specifically capture members of the microbiome active in the environment. Herein, we review those currently used activity-based approaches that can be directly linked to shotgun nucleic acid sequencing, evaluate their relevance to ecology studies, and discuss future directions.

Feasibility of Metatranscriptome Analysis from Infant Gut Microbiota: Adaptation to Solid Foods Results in Increased Activity of Firmicutes at Six Months

Newborns are rapidly colonized by microbes and their intestinal tracts contain highly dynamic and rapidly developing microbial communities in the first months of life. In this study, we describe the feasibility of isolating mRNA from rapidly processed faecal samples and applying deep RNA-Seq analysis to provide insight into the active contributors of the microbial community in early life. Specific attention is given to the impact of removing rRNA from the mRNA on the phylogenetic and transcriptional profiling and its analysis depth. A breastfed baby was followed in the first six months of life during adaptation to solid food, dairy products, and formula. It was found that, in the weaning period, the total transcriptional activity of Actinobacteria, mainly represented by Bifidobacterium, decreased while that of Firmicutes increased over time. Moreover, Firmicutes and Actinobacteria, including the canonical Bifidobacteria as well as Collinsella, were found to be important contributors to carbohydrate fermentation and vitamin biosynthesis in the infant intestine. Finally, the expression of Lactobacillus rhamnosus-like genes was detected, likely following transfer from the mother who consumed L. rhamnosus GG. The study indicates that metatranscriptome analysis of the infant gut microbiota is feasible on infant stool samples and can be used to provide insight into the core activities of the developing community.