Fatal tick-borne encephalitis virus infection in Dalmatian puppy-dogs after putative vector independent transmission

In a retrospective metatranscriptomics study, we identified tick-borne encephalitis virus (TBEV) to be the causative agent for a fatal non-suppurative meningoencephalitis in a three-week-old Dalmatian puppy in Switzerland. Further investigations showed that the two other littermates with similar signs and pathological lesions were also positive for TBEV. By using an unbiased approach of combining high-throughput sequencing (HTS) and bioinformatics we were able to solve the etiology and discover an unusual case of TBEV in three young puppies. Based on our findings, we suggest that a vector-independent transmission of TBEV occurred and that most likely an intrauterine infection led to the severe and fulminant disease of the entire litter. We were able to demonstrate the presence of TBEV RNA by in situ hybridization (ISH) in the brain of all three puppies. Furthermore, we were able to detect TBEV by RT-qPCR in total RNA extracted from formalin-fixed and paraffin embedded (FFPE) blocks containing multiple peripheral organs. Overall, our findings shed light on alternative vector-independent transmission routes of TBEV infections in dogs and encourage veterinary practitioners to consider TBEV as an important differential diagnosis in neurological cases in dogs.

Nematodes alter the taxonomic and functional profiles of benthic bacterial communities: A metatranscriptomic approach

Marine sediments cover 70% of the Earth's surface, and harbour diverse bacterial communities critical for marine biogeochemical processes, which affect climate change, biodiversity and ecosystem functioning. Nematodes, the most abundant and species-rich metazoan organisms in marine sediments, in turn, affect benthic bacterial communities and bacterial-mediated ecological processes, but the underlying mechanisms by which they affect biogeochemical cycles remain poorly understood. Here, we demonstrate using a metatranscriptomic approach that nematodes alter the taxonomic and functional profiles of benthic bacterial communities. We found particularly strong stimulation of nitrogen-fixing and methane-oxidizing bacteria in the presence of nematodes, as well as increased functional activity associated with methane metabolism and degradation of various carbon compounds. This study provides empirical evidence that the presence of nematodes results in taxonomic and functional shifts in active bacterial communities, indicating that nematodes may play an important role in benthic ecosystem processes.

Metatranscriptomics sheds light on the links between the functional traits of fungal guilds and ecological processes in forest soil ecosystems

Soil fungi belonging to different functional guilds, such as saprotrophs, pathogens, and mycorrhizal symbionts, play key roles in forest ecosystems. To date, no study has compared the actual gene expression of these guilds in different forest soils. We used metatranscriptomics to study the competition for organic resources by these fungal groups in boreal, temperate, and Mediterranean forest soils. Using a dedicated mRNA annotation pipeline combined with the JGI MycoCosm database, we compared the transcripts of these three fungal guilds, targeting enzymes involved in C- and N mobilization from plant and microbial cell walls. Genes encoding enzymes involved in the degradation of plant cell walls were expressed at a higher level in saprotrophic fungi than in ectomycorrhizal and pathogenic fungi. However, ectomycorrhizal and saprotrophic fungi showed similarly high expression levels of genes encoding enzymes involved in fungal cell wall degradation. Transcripts for N-related transporters were more highly expressed in ectomycorrhizal fungi than in other groups. We showed that ectomycorrhizal and saprotrophic fungi compete for N in soil organic matter, suggesting that their interactions could decelerate C cycling. Metatranscriptomics provides a unique tool to test controversial ecological hypotheses and to better understand the underlying ecological processes involved in soil functioning and carbon stabilization.

Unveiling the hidden viromes across the animal tree of life: insights from a taxonomic classification pipeline applied to invertebrates of 31 metazoan phyla

Invertebrates constitute the majority of animal species on Earth, including most disease-causing agents or vectors, with more diverse viromes when compared to vertebrates. Recent advancements in high-throughput sequencing have significantly expanded our understanding of invertebrate viruses, yet this knowledge remains biased toward a few well-studied animal lineages. In this study, we analyze invertebrate DNA and RNA viromes for 31 phyla using 417 publicly available RNA-Seq data sets from diverse environments in the marine-terrestrial and marine-freshwater gradients. This study aims to (i) estimate virome compositions at the family level for the first time across the animal tree of life, including the first exploration of the virome in several phyla, (ii) quantify the diversity of invertebrate viromes and characterize the structure of invertebrate-virus infection networks, and (iii) investigate host phylum and habitat influence on virome differences. Results showed that a set of few viral families of eukaryotes, comprising Retroviridae, Flaviviridae, and several families of giant DNA viruses, were ubiquitous and highly abundant. Nevertheless, some differences emerged between phyla, revealing for instance a less diverse virome in Ctenophora compared to the other animal phyla. Compositional analysis of the viromes showed that the host phylum explained over five times more variance in composition than its habitat. Moreover, significant similarities were observed between the viromes of some phylogenetically related phyla, which could highlight the influence of co-evolution in shaping invertebrate viromes.IMPORTANCEThis study significantly enhances our understanding of the global animal virome by characterizing the viromes of previously unexamined invertebrate lineages from a large number of animal phyla. It showcases the great diversity of viromes within each phylum and investigates the role of habitat shaping animal viral communities. Furthermore, our research identifies dominant virus families in invertebrates and distinguishes phyla with analogous viromes. This study sets the road toward a deeper understanding of the virome across the animal tree of life.

Composting reduces the risks of resistome in beef cattle manure at the transcriptional level

Transcriptomic evidence is needed to determine whether composting is more effective than conventional stockpiling in mitigating the risk of resistome in livestock manure. The objective of this study is to compare composting and stockpiling for their effectiveness in reducing the risk of antibiotic resistance in beef cattle manure. Samples collected from the center and the surface of full-size manure stockpiling and composting piles were subject to metagenomic and metatranscriptomic analyses. While the distinctions in resistome between stockpiled and composted manure were not evident at the DNA level, the advantages of composting over stockpiling were evident at the transcriptomic level in terms of the abundance of antibiotic resistance genes (ARGs), the number of ARG subtypes, and the prevalence of high-risk ARGs (i.e., mobile ARGs associated with zoonotic pathogens). DNA and transcript contigs show that the pathogen hosts of high-risk ARGs included Escherichia coli O157:H7 and O25b:H4, Klebsiella pneumoniae, and Salmonella enterica. Although the average daily temperatures for the entire composting pile exceeded 55°C throughout the field study, more ARG and ARG transcripts were removed at the center of the composting pile than at the surface. This work demonstrates the advantage of composting over stockpiling in reducing ARG risk in active populations in beef cattle manure.IMPORTANCEProper treatment of manure before land application is essential to mitigate the spread of antibiotic resistance in the environment. Stockpiling and composting are two commonly used methods for manure treatment. However, the effectiveness of composting in reducing antibiotic resistance in manure has been debated. This work compared the ability of these two methods to reduce the risk of antibiotic resistance in beef cattle manure. Our results demonstrate that composting reduced more high-risk resistance genes at the transcriptomic level in cattle manure than conventional stockpiling. This finding not only underscores the effectiveness of composting in reducing antibiotic resistance in manure but also highlights the importance of employing RNA analyses alongside DNA analyses.

Deep-sea in situ and laboratory multi-omics provide insights into the sulfur assimilation of a deep-sea Chloroflexota bacterium

Chloroflexota bacteria are abundant and globally distributed in various deep-sea ecosystems. It has been reported based on metagenomics data that two deep-sea Chloroflexota lineages (the SAR202 group and Dehalococcoidia class) have the potential to drive sulfur cycling. However, the absence of cultured Chloroflexota representatives is a significant bottleneck toward understanding their contribution to the deep-sea sulfur cycling. In this study, we find that Phototrophicus methaneseepsis ZRK33 isolated from deep-sea sediment has a heterotrophic lifestyle and can assimilate sulfate and thiosulfate. Using combined physiological, genomic, proteomic, and in situ transcriptomic methods, we find that strain ZRK33 can perform assimilatory sulfate reduction in both laboratory and deep-sea conditions. Metabolism of sulfate or thiosulfate by strain ZRK33 significantly promotes the transport and degradation of various macromolecules and thereby stimulates the energy production. In addition, metagenomic results show that genes associated with assimilatory and dissimilatory sulfate reduction are ubiquitously distributed in the metagenome-assembled genomes of Chloroflexota members derived from deep-sea sediments. Metatranscriptomic results also show that the expression levels of related genes are upregulated, strongly suggesting that Chloroflexota bacteria may play undocumented roles in deep-sea sulfur cycling.

IMPORTANCE

The cycling of sulfur is one of Earth's major biogeochemical processes and is closely related to the energy metabolism of microorganisms living in the deep-sea cold seep and hydrothermal vents. To date, some of the members of Chloroflexota are proposed to play a previously unrecognized role in sulfur cycling. However, the sulfur metabolic characteristics of deep-sea Chloroflexota bacteria have never been reported, and remain to be verified in cultured deep-sea representatives. Here, we show that the deep-sea Chloroflexota bacterium ZRK33 can perform sulfate assimilation in both laboratory and deep-sea conditions, which expands our knowledge of the sulfur metabolic potential of deep-sea Chloroflexota bacteria. We also show that the genes associated with assimilatory and dissimilatory sulfate reduction ubiquitously distribute in the deep-sea Chloroflexota members, providing hints to the roles of Chloroflexota bacteria in deep-sea sulfur biogeochemical cycling.

Microbial signatures in human periodontal disease: a metatranscriptome meta-analysis

The characterization of oral microbial communities and their functional potential has been shaped by metagenomics and metatranscriptomics studies. Here, a meta-analysis of four geographically and technically diverse oral shotgun metatranscriptomics studies of human periodontitis was performed. In total, 54 subgingival plaque samples, 27 healthy and 27 periodontitis, were analyzed. The core microbiota of the healthy and periodontitis group encompassed 40 and 80 species, respectively, with 38 species being common to both microbiota. The differential abundance analysis identified 23 genera and 26 species, that were more abundant in periodontitis. Our results not only validated previously reported genera and species associated with periodontitis with heightened statistical significance, but also elucidated additional genera and species that were overlooked in the individual studies. Functional analysis revealed a significant up-regulation in the transcription of 50 gene families (UniRef-90) associated with transmembrane transport and secretion, amino acid metabolism, surface protein and flagella synthesis, energy metabolism, and DNA supercoiling in periodontitis samples. Notably, the overwhelming majority of the identified gene families did not exhibit differential abundance when examined across individual datasets. Additionally, 4 bacterial virulence factor genes, including TonB dependent receptor from P. gingivalis, surface antigen BspA from T. forsynthia, and adhesin A (PsaA) and Type I glyceraldehyde-3-phosphate dehydrogenase (GAPDH) from the Streptococcus genus, were also found to be significantly more transcribed in periodontitis group. Microbial co-occurrence analysis demonstrated that the periodontitis microbial network was less dense compared to the healthy network, but it contained more positive correlations between the species. Furthermore, there were discernible disparities in the patterns of interconnections between the species in the two networks, denoting the rewiring of the whole microbial network during the transition to the disease state. In summary, our meta-analysis has provided robust insights into the oral active microbiome and transcriptome in both health and disease.

Heterogeneous lineage-specific arginine deiminase expression within dental microbiome species

Arginine catabolism by the bacterial arginine deiminase system (ADS) has anticariogenic properties through the production of ammonia, which modulates the pH of the oral environment. Given the potential protective capacity of the ADS pathway, the exploitation of ADS-competent oral microbes through pre- or probiotic applications is a promising therapeutic target to prevent tooth decay. To date, most investigations of the ADS in the oral cavity and its relation to caries have focused on indirect measures of activity or on specific bacterial groups, yet the pervasiveness and rate of expression of the ADS operon in diverse mixed microbial communities in oral health and disease remain an open question. Here, we use a multivariate approach, combining ultra-deep metatranscriptomic sequencing with paired metataxonomic and in vitro citrulline quantification to characterize the microbial community and ADS operon expression in healthy and late-stage cavitated teeth. While ADS activity is higher in healthy teeth, we identify multiple bacterial lineages with upregulated ADS activity on cavitated teeth that are distinct from those found on healthy teeth using both reference-based mapping and de novo assembly methods. Our dual metataxonomic and metatranscriptomic approach demonstrates the importance of species abundance for gene expression data interpretation and that patterns of differential expression can be skewed by low-abundance groups. Finally, we identify several potential candidate probiotic bacterial lineages within species that may be useful therapeutic targets for the prevention of tooth decay and propose that the development of a strain-specific, mixed-microbial probiotic may be a beneficial approach given the heterogeneity of taxa identified here across health groups.

IMPORTANCE

Tooth decay is the most common preventable chronic disease, affecting more than two billion people globally. The development of caries on teeth is primarily a consequence of acid production by cariogenic bacteria that inhabit the plaque microbiome. Other bacterial strains in the oral cavity may suppress or prevent tooth decay by producing ammonia as a byproduct of the arginine deiminase metabolic pathway, increasing the pH of the plaque biofilm. While the benefits of arginine metabolism on oral health have been extensively documented in specific bacterial groups, the prevalence and consistency of arginine deiminase system (ADS) activity among oral bacteria in a community context remain an open question. In the current study, we use a multi-omics approach to document the pervasiveness of the expression of the ADS operon in both health and disease to better understand the conditions in which ADS activity may prevent tooth decay.

Potential PCR amplification bias in identifying complex ecological patterns: Higher species compositional homogeneity revealed in smaller-size coral reef zooplankton by metatranscriptomics

PCR-based high-throughput sequencing has permitted comprehensive resolution analyses of zooplankton diversity dynamics. However, significant methodological issues still surround analyses of complex bulk community samples, not least as in prevailing PCR-based approaches. Marine drifting animals-zooplankton-play essential ecological roles in the pelagic ecosystem, transferring energy and elements to higher trophic levels, such as fishes, cetaceans and others. In the present study, we collected 48 size-fractionated zooplankton samples in the vicinity of a coral reef island with environmental gradients. To investigate the spatiotemporal dynamics of zooplankton diversity patterns and the effect of PCR amplification biases across these complex communities, we first took metatranscriptomics approach. Comprehensive computational analyses revealed a clear pattern of higher/lower homogeneity in smaller/larger zooplankton compositions across samples respectively. Our study thus suggests changes in the role of dispersal across the sizes. Next, we applied in silico PCR to the metatranscriptomics datasets, in order to estimate the extent of PCR amplification bias. Irrespective of stringency criteria, we observed clear separations of size fraction sample clusters in both metatranscriptomics and in silico datasets. In contrast, the pattern-smaller-fractioned communities had higher compositional homogeneity than larger ones-was observed in the metatranscriptomics data but not in the in silico datasets. To investigate this discrepancy further, we analysed the mismatches of widely used mitochondrial CO1 primers and identified priming site mismatches likely driving PCR-based biases. Our results suggest the use of metatranscriptomics or, although less ideal, redesigning the CO1 primers is necessary to circumvent these issues.

Metagenomic next-generation sequencing of nasopharyngeal microbiota in COVID-19 patients with different disease severities

Throughout the COVID-19 pandemic, extensive research has been conducted on SARS-COV-2 to elucidate its genome, prognosis, and possible treatments. However, few looked at the microbial markers that could be explored in infected patients and that could predict possible disease severity. The aim of this study is to compare the nasopharyngeal microbiota of healthy subjects, moderate, under medication, and recovered SARS-COV-2 patients. In 2020, 38 nasopharyngeal swabs were collected from 6 healthy subjects, 14 moderates, 10 under medication and 8 recovered SARS-COV-2 patients at King Fahad Medical city. Metatranscriptomic sequencing was performed using Minion Oxford nanopore sequencing. No significant difference in alpha as well as beta diversity was observed among all four categories. Nevertheless, we have found that Streptococcus spp including Streptococcus pneumoniae and Streptococcus thermophilus were among the top 15 most abundant species detected in COVID-19 patients but not in healthy subjects. The genus Staphylococcus was found to be associated with COVID-19 patients compared to healthy subjects. Furthermore, the abundance of Leptotrichia was significantly higher in healthy subjects compared to recovered patients. Corynebacterium on the other hand, was associated with under-medication patients. Taken together, our study revealed no differences in the overall microbial composition between healthy subjects and COVID-19 patients. Significant differences were seen only at specific taxonomic level. Future studies should explore the nasopharyngeal microbiota between controls and COVID-19 patients while controlling for confounders including age, gender, and comorbidities; since these latter could affect the results and accordingly the interpretation.IMPORTANCEIn this work, no significant difference in the microbial diversity was seen between healthy subjects and COVID-19 patients. Changes in specific taxa including Leptotrichia, Staphylococcus, and Corynebacterium were only observed. Leptotrichia was significantly higher in healthy subjects, whereas Staphylococcus and Corynebacterium were mostly associated with COVID-19, and specifically with under-medication SARS-COV-2 patients, respectively. Although the COVID-19 pandemic has ended, the SARS-COV-2 virus is continuously evolving and the emergence of new variants causing more severe disease should be always kept in mind. Microbial markers in SARS-COV-2 infected patients can be useful in the early suspicion of the disease, predicting clinical outcomes, framing hospital and intensive care unit admission as well as, risk stratification. Data on which microbial marker to tackle is still controversial and more work is needed, hence the importance of this study.

Host adaptive radiation is associated with rapid virus diversification and cross-species transmission in African cichlid fishes

Adaptive radiations are generated through a complex interplay of biotic and abiotic factors. Although adaptive radiations have been widely studied in the context of animal and plant evolution, little is known about how they impact the evolution of the viruses that infect these hosts, which in turn may provide insights into the drivers of cross-species transmission and hence disease emergence. We examined how the rapid adaptive radiation of the cichlid fishes of African Lake Tanganyika over the last 10 million years has shaped the diversity and evolution of the viruses they carry. Through metatranscriptomic analysis of 2,242 RNA sequencing libraries, we identified 121 vertebrate-associated viruses among various tissue types that fell into 13 RNA and 4 DNA virus groups. Host-switching was commonplace, particularly within the Astroviridae, Metahepadnavirus, Nackednavirus, Picornaviridae, and Hepacivirus groups, occurring more frequently than in other fish communities. A time-calibrated phylogeny revealed that hepacivirus diversification was not constant throughout the cichlid radiation but accelerated 2-3 million years ago, coinciding with a period of rapid cichlid diversification and niche packing in Lake Tanganyika, thereby providing more closely related hosts for viral infection. These data depict a dynamic virus ecosystem within the cichlids of Lake Tanganyika, characterized by rapid virus diversification and frequent host jumping, and likely reflecting their close phylogenetic relationships that lower the barriers to cross-species virus transmission.

A novel method for sampling subgingival microbiome: a comparative metatranscriptomic study

The subgingival microbiome has been implicated in oral and systemic diseases such as periodontitis and Alzheimer's disease. However, subgingival sampling is challenging. We developed a novel method of sampling the subgingival microbiome by rotationally swabbing the supragingival area, named subgingival-P (for proxy) samples. We sampled and metatranscriptomically analyzed subgingival and subgingival-P samples of three different teeth in 20 individuals. The subgingival-P samples were comparable to the subgingival samples in the relative abundances of microorganisms and microbial gene expression levels. Our data demonstrate that the novel method of collecting and analyzing the subgingival-P samples can act as a proxy for the subgingiva, paving the way for large and diverse studies investigating the role of the subgingival microbiome in health and disease.

Syntrophic microbes involved in the oxidation of short-chain fatty acids in continuous-flow anaerobic digesters treating waste activated sludge with hydrochar

The rapid degradation of short-chain fatty acids (SCFAs) is an essential issue of anaerobic digestion (AD), in which SCFA oxidizers could generally metabolize in syntrophy with methanogens. The dynamic responses of active metagenome-assembled genomes to low concentrations of propionate and acetate were analyzed to identify specific syntrophic SCFA oxidizers and their metabolic characteristics in continuous-flow AD systems treating waste activated sludge with and without hydrochar. In this study, hydrochar increased methane production by 19%, possibly due to hydrochar enhancing acidification and methanogenesis processes. A putative syntrophic propionate oxidizer and two acetate oxidizers contributed substantially to the syntrophic degradation of SCFAs, and hydrochar positively regulated their functional gene expressions. A significant relationship was established between the replication rate of SCFA oxidizers and their stimulation-related transcriptional activity. Acetate was degraded in the hydrochar group, which might be mainly through the syntrophic acetate oxidizer from the genus Desulfallas and methanogens from the genus Methanosarcina.IMPORTANCEShort-chain fatty acid (SCFA) degradation is an important process in the methanogenic ecosystem. However, current knowledge of this microbial mechanism is mainly based on studies on a few model organisms incubated as mono- or co-cultures or in enrichments, which cannot provide appropriate evidence in complex environments. Here, this study revealed the microbial mechanism of a hydrochar-mediated anaerobic digestion (AD) system promoting SCFA degradation at the species level and identified key SCFA oxidizing bacteria. Our analysis provided new insights into the SCFA oxidizers involved in the AD of waste activated sludge facilitated by hydrochar.

Cover Crop Root Exudates Impact Soil Microbiome Functional Trajectories in Agricultural Soils

Background

Cover cropping is an agricultural practice that uses secondary crops to support the growth of primary crops through various mechanisms including erosion control, weed suppression, nutrient management, and enhanced biodiversity. Cover crops may elicit some of these ecosystem services through chemical interactions with the soil microbiome via root exudation, or the release of plant metabolites from roots. Phytohormones are one metabolite type exuded by plants that activate the rhizosphere microbiome, yet managing this chemical interaction remains an untapped mechanism for optimizing plant-soil microbiome interactions. Currently, there is limited understanding on the diversity of cover crop phytohormone root exudation patterns and how these chemical messages selectively enrich specific microbial taxa and functionalities in agricultural soils.

Results

Here, we link variability in cover crop root exudate composition to changes in soil microbiome functionality. Exudate chemical profiles from 4 cover crop species (Sorghum bicolor, Vicia villosa, Brassica napus, and Secale cereal) were used as the chemical inputs to decipher microbial responses. These distinct exudate profiles, along with a no exudate control, were amended to agricultural soil microcosms with microbial responses tracked over time using metabolomes and genome-resolved metatranscriptomes. Our findings illustrated microbial metabolic patterns were unique in response to cover crop exudate inputs over time, particularly by sorghum and cereal rye amended microcosms where we identify novel microbial members (at the genera and family level) who produced IAA and GA4 over time. We also identify broad changes in microbial nitrogen cycling in response chemical inputs.

Conclusions

We highlight that root exudate amendments alter microbial community function and phytohormone metabolisms, particularly in response to root exudates isolated from cereal rye and sorghum plants. Additionally, we constructed a soil microbial genomic catalog of microorganisms responding to commonly used cover crops, a public resource for agriculturally-relevant microbes. Many of our exudate-stimulated microorganisms are representatives from poorly characterized or novel taxa, highlighting the yet to be discovered metabolic reservoir harbored in agricultural soils. Our findings emphasize the tractability of high-resolution multiomics approaches to investigate processes relevant for agricultural soils, opening the possibility of targeting specific soil biogeochemical outcomes through biological precision agricultural practices that use cover crops and the microbiome as levers for enhanced crop production.

Quorum Sensing Enhances Direct Interspecies Electron Transfer in Anaerobic Methane Production

Direct interspecies electron transfer (DIET) provides an innovative way to achieve efficient methanogenesis, and this study proposes a new approach to upregulate the DIET pathway by enhancing quorum sensing (QS). Based on long-term reactor performance, QS enhancement achieved more vigorous methanogenesis with 98.7% COD removal efficiency. In the control system, methanogenesis failure occurred at the accumulated acetate of 7420 mg of COD/L and lowered pH of 6.04, and a much lower COD removal of 41.9% was observed. The more significant DIET in QS-enhancing system was supported by higher expression of conductive pili and the c-Cyts cytochrome secretion-related genes, resulting in 12.7- and 10.3-fold improvements. Moreover, QS enhancement also improved the energy production capability, with the increase of F-type and V/A-type ATPase expression by 6.3- and 4.2-fold, and this effect probably provided more energy for nanowires and c-Cyts cytochrome secretion. From the perspective of community structure, QS enhancement increased the abundance of Methanosaeta and Geobacter from 54.3 and 17.6% in the control to 63.0 and 33.8%, respectively. Furthermore, the expression of genes involved in carbon dioxide reduction and alcohol dehydrogenation increased by 0.6- and 7.1-fold, respectively. Taken together, this study indicates the positive effects of QS chemicals to stimulate DIET and advances the understanding of the DIET methanogenesis involved in environments such as anaerobic digesters and sediments.

Recent advances in applying omic technologies for studying acetic acid bacteria in industrial vinegar production: A comprehensive review

Vinegar and related bioproducts containing acetic acid as the main component are among the most appreciated fermented foodstuffs in numerous European and Asian countries because of their exceptional organoleptic and bio-healthy properties. Regarding the acetification process and obtaining of final products, there is still a lack of knowledge on fundamental aspects, especially those related to the study of biodiversity and metabolism of the present microbiota. In this context, omic technologies currently allow for the massive analysis of macromolecules and metabolites for the identification and characterization of these microorganisms working in their natural media without the need for isolation. This review approaches comprehensive research on the application of omic tools for the identification of vinegar microbiota, mainly acetic acid bacteria, with subsequent emphasis on the study of the microbial diversity, behavior, and key molecular strategies used by the predominant groups throughout acetification. The current omics tools are enabling both the finding of new vinegar microbiota members and exploring underlying strategies during the elaboration process. The species Komagataeibacter europaeus may be a model organism for present and future research in this industry; moreover, the development of integrated meta-omic analysis may facilitate the achievement of numerous of the proposed milestones. This work might provide useful guidance for the vinegar industry establishing the first steps towards the improvement of the acetification conditions and the development of new products with sensory and bio-healthy profiles adapted to the agri-food market.

Virome analysis of New Zealand's bats reveals cross-species viral transmission among the Coronaviridae

The lesser short-tailed bat (Mystacina tuberculata) and the long-tailed bat (Chalinolobus tuberculatus) are Aotearoa New Zealand's only native extant terrestrial mammals and are believed to have migrated from Australia. Long-tailed bats arrived in New Zealand an estimated two million years ago and are closely related to other Australian bat species. Lesser short-tailed bats, in contrast, are the only extant species within the Mystacinidae and are estimated to have been living in isolation in New Zealand for the past 16-18 million years. Throughout this period of isolation, lesser short-tailed bats have become one of the most terrestrial bats in the world. Through a metatranscriptomic analysis of guano samples from eight locations across New Zealand, we aimed to characterise the viromes of New Zealand's bats and determine whether viruses have jumped between these species over the past two million years. High viral richness was observed among long-tailed bats with viruses spanning seven different viral families. In contrast, no bat-specific viruses were identified in lesser short-tailed bats. Both bat species harboured an abundance of likely dietary- and environment-associated viruses. We also identified alphacoronaviruses in long-tailed bat guano that had previously been identified in lesser short-tailed bats, suggesting that these viruses had jumped the species barrier after long-tailed bats migrated to New Zealand. Of note, an alphacoronavirus species discovered here possessed a complete genome of only 22,416 nucleotides with entire deletions or truncations of several non-structural proteins, thereby representing what may be the shortest genome within the Coronaviridae identified to date. Overall, this study has revealed a diverse range of novel viruses harboured by New Zealand's only native terrestrial mammals, in turn expanding our understanding of bat viral dynamics and evolution globally.

Metatranscriptomic Sequencing of Medically Important Mosquitoes Reveals Extensive Diversity of RNA Viruses and Other Microbial Communities in Western Australia

Mosquitoes harbor a wide diversity of microorganisms, including viruses that are human pathogens, or that are insect specific. We used metatranscriptomics, an unbiased high-throughput molecular approach, to describe the composition of viral and other microbial communities in six medically important mosquito species from across Western Australia: Aedes vigilax, Culex annulirostris, Cx. australicus, Cx. globocoxitus, Cx. pipiens biotype molestus, and Cx. quinquefasciatus. We identified 42 viral species, including 13 novel viruses, from 19 families. Culex mosquitoes exhibited a significantly higher diversity of viruses than Aedes mosquitoes, and no virus was shared between the two genera. Comparison of mosquito populations revealed a heterogenous distribution of viruses between geographical regions and between closely related species, suggesting that geography and host species may play a role in shaping virome composition. We also detected bacterial and parasitic microorganisms, among which Wolbachia bacteria were detected in three members of the Cx. pipiens complex, Cx. australicus, Cx. pipiens biotype molestus, and Cx. quinquefasciatus. In summary, our unbiased metatranscriptomics approach provides important insights into viral and other microbial diversity in Western Australian mosquitoes that vector medically important viruses.

Poly-omic risk scores predict inflammatory bowel disease diagnosis

Inflammatory bowel disease (IBD) is characterized by complex etiology and a disrupted colonic ecosystem. We provide a framework for the analysis of multi-omic data, which we apply to study the gut ecosystem in IBD. Specifically, we train and validate models using data on the metagenome, metatranscriptome, virome, and metabolome from the Human Microbiome Project 2 IBD multi-omic database, with 1,785 repeated samples from 130 individuals (103 cases and 27 controls). After splitting the participants into training and testing groups, we used mixed-effects least absolute shrinkage and selection operator regression to select features for each omic. These features, with demographic covariates, were used to generate separate single-omic prediction scores. All four single-omic scores were then combined into a final regression to assess the relative importance of the individual omics and the predictive benefits when considered together. We identified several species, pathways, and metabolites known to be associated with IBD risk, and we explored the connections between data sets. Individually, metabolomic and viromic scores were more predictive than metagenomics or metatranscriptomics, and when all four scores were combined, we predicted disease diagnosis with a Nagelkerke's R2 of 0.46 and an area under the curve of 0.80 (95% confidence interval: 0.63, 0.98). Our work supports that some single-omic models for complex traits are more predictive than others, that incorporating multiple omic data sets may improve prediction, and that each omic data type provides a combination of unique and redundant information. This modeling framework can be extended to other complex traits and multi-omic data sets.IMPORTANCEComplex traits are characterized by many biological and environmental factors, such that multi-omic data sets are well-positioned to help us understand their underlying etiologies. We applied a prediction framework across multiple omics (metagenomics, metatranscriptomics, metabolomics, and viromics) from the gut ecosystem to predict inflammatory bowel disease (IBD) diagnosis. The predicted scores from our models highlighted key features and allowed us to compare the relative utility of each omic data set in single-omic versus multi-omic models. Our results emphasized the importance of metabolomics and viromics over metagenomics and metatranscriptomics for predicting IBD status. The greater predictive capability of metabolomics and viromics is likely because these omics serve as markers of lifestyle factors such as diet. This study provides a modeling framework for multi-omic data, and our results show the utility of combining multiple omic data types to disentangle complex disease etiologies and biological signatures.

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

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

Decoding Microbial Responses to Ammonia Shock Loads in Biogas Reactors through Metagenomics and Metatranscriptomics

The presence of elevated ammonia levels is widely recognized as a significant contributor to process inhibition in biogas production, posing a common challenge for biogas plant operators. The present study employed a combination of biochemical, genome-centric metagenomic and metatranscriptomic data to investigate the response of the biogas microbiome to two shock loads induced by single pulses of elevated ammonia concentrations (i.e., 1.5 g NH4+/LR and 5 g NH4+/LR). The analysis revealed a microbial community of high complexity consisting of 364 Metagenome Assembled Genomes (MAGs). The hydrogenotrophic pathway was the primary route for methane production during the entire experiment, confirming its efficiency even at high ammonia concentrations. Additionally, metatranscriptomic analysis uncovered a metabolic shift in the methanogens Methanothrix sp. MA6 and Methanosarcina flavescens MX5, which switched their metabolism from the acetoclastic to the CO2 reduction route during the second shock. Furthermore, multiple genes associated with mechanisms for maintaining osmotic balance in the cell were upregulated, emphasizing the critical role of osmoprotection in the rapid response to the presence of ammonia. Finally, this study offers insights into the transcriptional response of an anaerobic digestion community, specifically focusing on the mechanisms involved in recovering from ammonia-induced stress.

Taxonomic and functional stability overrules seasonality in polar benthic microbiomes

Coastal shelf sediments are hot spots of organic matter mineralization. They receive up to 50% of primary production, which, in higher latitudes, is strongly seasonal. Polar and temperate benthic bacterial communities, however, show a stable composition based on comparative 16S rRNA gene sequencing despite different microbial activity levels. Here, we aimed to resolve this contradiction by identifying seasonal changes at the functional level, in particular with respect to algal polysaccharide degradation genes, by combining metagenomics, metatranscriptomics, and glycan analysis in sandy surface sediments from Isfjorden, Svalbard. Gene expressions of diverse carbohydrate-active enzymes changed between winter and spring. For example, β-1,3-glucosidases (e.g. GH30, GH17, GH16) degrading laminarin, an energy storage molecule of algae, were elevated in spring, while enzymes related to α-glucan degradation were expressed in both seasons with maxima in winter (e.g. GH63, GH13_18, and GH15). Also, the expression of GH23 involved in peptidoglycan degradation was prevalent, which is in line with recycling of bacterial biomass. Sugar extractions from bulk sediments were low in concentrations during winter but higher in spring samples, with glucose constituting the largest fraction of measured monosaccharides (84% ± 14%). In porewater, glycan concentrations were ~18-fold higher than in overlying seawater (1107 ± 484 vs. 62 ± 101 μg C l-1) and were depleted in glucose. Our data indicate that microbial communities in sandy sediments digest and transform labile parts of photosynthesis-derived particulate organic matter and likely release more stable, glucose-depleted residual glycans of unknown structures, quantities, and residence times into the ocean, thus modulating the glycan composition of marine coastal waters.

Megataxonomy and global ecology of the virosphere

Nearly all organisms are hosts to multiple viruses that collectively appear to be the most abundant biological entities in the biosphere. With recent advances in metagenomics and metatranscriptomics, the known diversity of viruses substantially expanded. Comparative analysis of these viruses using advanced computational methods culminated in the reconstruction of the evolution of major groups of viruses and enabled the construction of a virus megataxonomy, which has been formally adopted by the International Committee on Taxonomy of Viruses. This comprehensive taxonomy consists of six virus realms, which are aspired to be monophyletic and assembled based on the conservation of hallmark proteins involved in capsid structure formation or genome replication. The viruses in different major taxa substantially differ in host range and accordingly in ecological niches. In this review article, we outline the latest developments in virus megataxonomy and the recent discoveries that will likely lead to reassessment of some major taxa, in particular, split of three of the current six realms into two or more independent realms. We then discuss the correspondence between virus taxonomy and the distribution of viruses among hosts and ecological niches, as well as the abundance of viruses versus cells in different habitats. The distribution of viruses across environments appears to be primarily determined by the host ranges, i.e. the virome is shaped by the composition of the biome in a given habitat, which itself is affected by abiotic factors.

Biotic interactions explain seasonal dynamics of the alpine soil microbiome

While it is acknowledged that alpine soil bacterial communities are primarily driven by season and elevation, there is no consensus on the factors influencing fungi and protists. Here we used a holistic approach of the microbiome to investigate the seasonal dynamics in alpine grasslands, focusing on soil food web interactions. We collected 158 soil samples along elevation transects from three mountains in the Alps, in spring during snowmelt and in the following summer. Using metatranscriptomics, we simultaneously assessed prokaryotic and eukaryotic communities, further classified into trophic guilds. Our findings reveal that the consumers' pressure increases from spring to summer, leading to more diverse and evenly distributed prey communities. Consequently, consumers effectively maintain the diverse soil bacterial and fungal communities essential for ecosystem functioning. Our research highlights the significance of biotic interactions in understanding the distribution and dynamics of alpine microbial communities.

Host species and geography impact bee-associated RNA virus communities with evidence for isolation by distance in viral populations

Virus symbionts are important mediators of ecosystem function, yet we know little of their diversity and ecology in natural populations. The alarming decline of pollinating insects in many regions of the globe, especially the European honey bee, Apis mellifera, has been driven in part by worldwide transmission of virus pathogens. Previous work has examined the transmission of known honey bee virus pathogens to wild bee populations, but only a handful of studies have investigated the native viromes associated with wild bees, limiting epidemiological predictors associated with viral pathogenesis. Further, variation among different bee species might have important consequences in the acquisition and maintenance of bee-associated virome diversity. We utilized comparative metatranscriptomics to develop a baseline description of the RNA viromes associated with wild bee pollinators and to document viral diversity, community composition, and structure. Our sampling includes five wild-caught, native bee species that vary in social behavior as well as managed honey bees. We describe 26 putatively new RNA virus species based on RNA-dependent RNA polymerase phylogeny and show that each sampled bee species was associated with a specific virus community composition, even among sympatric populations of distinct host species. From 17 samples of a single host species, we recovered a single virus species despite over 600 km of distance between host populations and found strong evidence for isolation by distance in associated viral populations. Our work adds to the small number of studies examining viral prevalence and community composition in wild bees.

A metatranscriptomics strategy for efficient characterization of the microbiome in human tissues with low microbial biomass

The high background of host RNA poses a major challenge to metatranscriptome analysis of human samples. Hence, metatranscriptomics has been mainly applied to microbe-rich samples, while its application in human tissues with low ratio of microbial to host cells has yet to be explored. Since there is no computational workflow specifically designed for the taxonomic and functional analysis of this type of samples, we propose an effective metatranscriptomics strategy to accurately characterize the microbiome in human tissues with a low ratio of microbial to host content. We experimentally generated synthetic samples with well-characterized bacterial and host cell compositions, and mimicking human samples with high and low microbial loads. These synthetic samples were used for optimizing and establishing the workflow in a controlled setting. Our results show that the integration of the taxonomic analysis of optimized Kraken 2/Bracken with the functional analysis of HUMAnN 3 in samples with low microbial content, enables the accurate identification of a large number of microbial species with a low false-positive rate, while improving the detection of microbial functions. The effectiveness of our metatranscriptomics workflow was demonstrated in synthetic samples, simulated datasets, and most importantly, human gastric tissue specimens, thus providing a proof of concept for its applicability on mucosal tissues of the gastrointestinal tract. The use of an accurate and reliable metatranscriptomics approach for human tissues with low microbial content will expand our understanding of the functional activity of the mucosal microbiome, uncovering critical interactions between the microbiome and the host in health and disease.

Metatranscriptomics Reveals the RNA Virome of Ixodes Persulcatus in the China-North Korea Border, 2017

In recent years, numerous viruses have been identified from ticks, and some have been linked to clinical cases of emerging tick-borne diseases. Chinese northeast frontier is tick infested. However, there is a notable lack of systematic monitoring efforts to assess the viral composition in the area, leaving the ecological landscape of viruses carried by ticks not clear enough. Between April and June 2017, 7101 ticks were collected to perform virus surveillance on the China-North Korea border, specifically in Tonghua, Baishan, and Yanbian. A total of 2127 Ixodes persulcatus were identified. Further investigation revealed the diversity of tick-borne viruses by transcriptome sequencing of Ixodes persulcatus. All ticks tested negative for tick-borne encephalitis virus. Transcriptome sequencing expanded 121 genomic sequence data of 12 different virus species from Ixodes persulcatus. Notably, a new segmented flavivirus, named Baishan Forest Tick Virus, were identified, closely related to Alongshan virus and Harz mountain virus. Therefore, this new virus may pose a potential threat to humans. Furthermore, the study revealed the existence of seven emerging tick-borne viruses dating back to 2017. These previously identified viruses included Mudanjiang phlebovirus, Onega tick phlebovirus, Sara tick phlebovirus, Yichun mivirus, and three unnamed viruses (one belonging to the Peribunyaviridae family and the other two belonging to the Phenuiviridae family). The existence of these emerging tick-borne viruses in tick samples collected in 2017 suggests that their history may extend further than previously recognized. This study provides invaluable insights into the virome of Ixodes persulcatus in the China-North Korea border region, enhancing our ongoing efforts to manage the risks associated with tick-borne viruses.

Novel Tri-Segmented Rhabdoviruses: A Data Mining Expedition Unveils the Cryptic Diversity of Cytorhabdoviruses

Cytorhabdoviruses (genus Cytorhabdovirus, family Rhabdoviridae) are plant-infecting viruses with enveloped, bacilliform virions. Established members of the genus Cytorhabdovirus have unsegmented single-stranded negative-sense RNA genomes (ca. 10-16 kb) which encode four to ten proteins. Here, by exploring large publicly available metatranscriptomics datasets, we report the identification and genomic characterization of 93 novel viruses with genetic and evolutionary cues of cytorhabdoviruses. Strikingly, five unprecedented viruses with tri-segmented genomes were also identified. This finding represents the first tri-segmented viruses in the family Rhabdoviridae, and they should be classified in a novel genus within this family for which we suggest the name "Trirhavirus". Interestingly, the nucleocapsid and polymerase were the only typical rhabdoviral proteins encoded by those tri-segmented viruses, whereas in three of them, a protein similar to the emaravirus (family Fimoviridae) silencing suppressor was found, while the other predicted proteins had no matches in any sequence databases. Genetic distance and evolutionary insights suggest that all these novel viruses may represent members of novel species. Phylogenetic analyses, of both novel and previously classified plant rhabdoviruses, provide compelling support for the division of the genus Cytorhabdovirus into three distinct genera. This proposed reclassification not only enhances our understanding of the evolutionary dynamics within this group of plant rhabdoviruses but also illuminates the remarkable genomic diversity they encompass. This study not only represents a significant expansion of the genomics of cytorhabdoviruses that will enable future research on the evolutionary peculiarity of this genus but also shows the plasticity in the rhabdovirus genome organization with the discovery of tri-segmented members with a unique evolutionary trajectory.

Bridging the Gaps in Meta-Omic Analysis: Workflows and Reproducibility

The past few years have seen significant advances in the study of complex microbial communities associated with the evolution of sequencing technologies and increasing adoption of whole genome shotgun sequencing methods over the once more traditional Amplicon-based methods. Although these advances have broadened the horizon of meta-omic analyses in planetary health, human health, and ecology from simple sample composition studies to comprehensive taxonomic and metabolic profiles, there are still significant challenges in processing these data. First, there is a widespread lack of standardization in data processing, including software choices and the ease of installing and running attendant software. This can lead to several inconsistencies, making comparing results across studies and reproducing original results difficult. We argue that these drawbacks are especially evident in metatranscriptomic analysis, with most analyses relying on ad hoc scripts instead of pipelines implemented in workflow managers. Additional challenges rely on integrating meta-omic data, since methods have to consider the biases in the library preparation and sequencing methods and the technical noise that can arise from it. Here, we critically discuss the current limitations in metagenomics and metatranscriptomics methods with a view to catalyze future innovations in the field of Planetary Health, ecology, and allied fields of life sciences. We highlight possible solutions for these constraints to bring about more standardization, with ease of installation, high performance, and reproducibility as guiding principles.

Exploring the Role of Active Functional Microbiota in Flavor Generation by Integrated Metatranscriptomics and Metabolomics during Niulanshan Baijiu Fermentation

Active functional microbiota for producing volatile flavors is critical to Chinese baijiu fermentation. Microbial communities correlated with the volatile metabolites are generally explored using DNA-based sequencing and metabolic analysis. However, the active functional microbiota related to the volatile flavor compounds is poorly understood. In this study, an integrated metatranscriptomic and metabolomics analysis was employed to unravel the metabolite profiles comprehensively and the contributing active functional microbiota for flavor generation during Niulanshan baijiu fermentation. A total of 395, 83, and 181 compounds were annotated using untargeted metabolomics, including LC-MS, GC-MS, and HS-SPME-GC-MS, respectively. Significant variances were displayed in the composition of compounds among different time-point samples according to the heatmaps and orthogonal partial least-square discriminant analysis. The correlation between the active microbiota and the volatile flavors was analyzed based on the bidirectional orthogonal partial least squares discriminant analysis (O2PLS-DA) model. Six bacterial genera, including Streptococcus, Lactobacillus, Pediococcus, Campylobacter, Yersinia, and Weissella, and five fungal genera of Talaromyces, Aspergillus, Mixia, Rhizophagus, and Gloeophyllum were identified as the active functional microbiota for producing the volatile flavors. In summary, this study revealed the active functional microbial basis of unique flavor formation and provided novel insights into the optimization of Niulanshan baijiu fermentation.

Lower Airway Dysbiosis Augments Lung Inflammatory Injury in Mild-to-Moderate Chronic Obstructive Pulmonary Disease

Rationale: Chronic obstructive pulmonary disease (COPD) is associated with high morbidity, mortality, and healthcare costs. Cigarette smoke is a causative factor; however, not all heavy smokers develop COPD. Microbial colonization and infections are contributing factors to disease progression in advanced stages. Objectives: We investigated whether lower airway dysbiosis occurs in mild-to-moderate COPD and analyzed possible mechanistic contributions to COPD pathogenesis. Methods: We recruited 57 patients with a >10 pack-year smoking history: 26 had physiological evidence of COPD, and 31 had normal lung function (smoker control subjects). Bronchoscopy sampled the upper airways, lower airways, and environmental background. Samples were analyzed by 16S rRNA gene sequencing, whole genome, RNA metatranscriptome, and host RNA transcriptome. A preclinical mouse model was used to evaluate the contributions of cigarette smoke and dysbiosis on lower airway inflammatory injury. Measurements and Main Results: Compared with smoker control subjects, microbiome analyses showed that the lower airways of subjects with COPD were enriched with common oral commensals. The lower airway host transcriptomics demonstrated differences in markers of inflammation and tumorigenesis, such as upregulation of IL-17, IL-6, ERK/MAPK, PI3K, MUC1, and MUC4 in mild-to-moderate COPD. Finally, in a preclinical murine model exposed to cigarette smoke, lower airway dysbiosis with common oral commensals augments the inflammatory injury, revealing transcriptomic signatures similar to those observed in human subjects with COPD. Conclusions: Lower airway dysbiosis in the setting of smoke exposure contributes to inflammatory injury early in COPD. Targeting the lower airway microbiome in combination with smoking cessation may be of potential therapeutic relevance.

Reduced chemosymbiont genome in the methane seep thyasirid and the cooperated metabolisms in the holobiont under anaerobic sediment

Previous studies have deciphered the genomic basis of host-symbiont metabolic complementarity in vestimentiferans, bathymodioline mussels, vesicomyid clams and Alviniconcha snails, yet little is known about the chemosynthetic symbiosis in Thyasiridae-a family of Bivalvia regarded as an excellent model in chemosymbiosis research due to their wide distribution in both deep-sea and shallow-water habitats. We report the first circular thyasirid symbiont genome, named Candidatus Ruthturnera sp. Tsphm01, with a size of 1.53 Mb, 1521 coding genes and 100% completeness. Compared to its free-living relatives, Ca. Ruthturnera sp. Tsphm01 genome is reduced, lacking components for chemotaxis, citric acid cycle and de novo biosynthesis of small molecules (e.g. amino acids and cofactors), indicating it is likely an obligate intracellular symbiont. Nevertheless, the symbiont retains complete genomic components of sulphur oxidation and assimilation of inorganic carbon, and these systems were highly and actively expressed. Moreover, the symbiont appears well-adapted to anoxic environment, including capable of anaerobic respiration (i.e. reductions of DMSO and nitrate) and possession of a low oxygen-adapted type of cytochrome c oxidase. Analysis of the host transcriptome revealed its metabolic complementarity to the incomplete metabolic pathways of the symbiont and the acquisition of nutrients from the symbiont via phagocytosis and exosome. By providing the first complete genome of reduced size in a thyasirid symbiont, this study enhances our understanding of the diversity of symbiosis that has enabled bivalves to thrive in chemosynthetic habitats. The resources will be widely used in phylogenetic, geographic and evolutionary studies of chemosynthetic bacteria and bivalves.

Bioinformatic approaches for studying the microbiome of fermented food

High-throughput DNA sequencing-based approaches continue to revolutionise our understanding of microbial ecosystems, including those associated with fermented foods. Metagenomic and metatranscriptomic approaches are state-of-the-art biological profiling methods and are employed to investigate a wide variety of characteristics of microbial communities, such as taxonomic membership, gene content and the range and level at which these genes are expressed. Individual groups and consortia of researchers are utilising these approaches to produce increasingly large and complex datasets, representing vast populations of microorganisms. There is a corresponding requirement for the development and application of appropriate bioinformatic tools and pipelines to interpret this data. This review critically analyses the tools and pipelines that have been used or that could be applied to the analysis of metagenomic and metatranscriptomic data from fermented foods. In addition, we critically analyse a number of studies of fermented foods in which these tools have previously been applied, to highlight the insights that these approaches can provide.

Targeted curation of the gut microbial gene content modulating human cardiovascular disease

IMPORTANCE

One of the most-cited examples of the gut microbiome modulating human disease is the microbial metabolism of quaternary amines from protein-rich foods. By-products of this microbial processing promote atherosclerotic heart disease, a leading cause of human mortality globally. Our research addresses current knowledge gaps in our understanding of this microbial metabolism by holistically inventorying the microorganisms and expressed genes catalyzing critical atherosclerosis-promoting and -ameliorating reactions in the human gut. This led to the creation of an open-access resource, the Methylated Amine Gene Inventory of Catabolism database, the first systematic inventory of gut methylated amine metabolism. More importantly, using this resource we deliver here, we show for the first time that these gut microbial genes can predict human disease, paving the way for microbiota-inspired diagnostics and interventions.

Proton Pump Inhibitors Modulate Gene Expression Profile in Esophageal Mucosa and Microbiome

OBJECTIVE

Proton pump inhibitors (PPIs) are commonly used to manage children with upper gastrointestinal symptoms and without a formal diagnosis. We investigated the effect of PPIs on esophageal mucosal transcriptome and active microbiota in children with normal esophagi. Furthermore, we examined whether the differences in host esophageal mucosal gene expression were driven by an underlying esophageal epithelial cell type composition.

METHODS

Using metatranscriptomics, the host transcriptional and active microbial profiles were captured from 17 esophageal biopsy samples (PPI naïve [PPI-], n = 7; PPI exposed [PPI+], n = 10) collected from children without any endoscopic and histologic abnormalities in their esophagus (normal esophagus). Deconvolution computational analysis was performed with xCell to assess if the observed epithelial gene expression changes were related to the cell type composition in the esophageal samples.

RESULTS

The median (IQR) age of our cohort was 14 years (12-16) with female (63%) preponderance. Both groups were similar in terms of their demographics and clinical features. Compared with PPI-, the PPI+ had upregulation of 27 genes including the MUC genes. The cell type composition was similar between the PPI- and PPI+ groups. Prevotella sp and Streptococcus sp were abundant in PPI+ group.

CONCLUSIONS

In children with normal esophagus, PPI exposure can be associated with upregulation of esophageal mucosal homeostasis and epithelial cell function genes in a cell-type independent manner, and an altered esophageal microbiome. Additional studies are warranted to validate our findings and to investigate the causal effect of PPIs on the normal esophageal epithelium and microbial communities.

Rhizosphere and detritusphere habitats modulate expression of soil N-cycling genes during plant development

IMPORTANCE

Plant roots modulate microbial nitrogen (N) cycling by regulating the supply of root-derived carbon and nitrogen uptake. These differences in resource availability cause distinct micro-habitats to develop: soil near living roots, decaying roots, near both, or outside the direct influence of roots. While many environmental factors and genes control the microbial processes involved in the nitrogen cycle, most research has focused on single genes and pathways, neglecting the interactive effects these pathways have on each other. The processes controlled by these pathways determine consumption and production of N by soil microorganisms. We followed the expression of N-cycling genes in four soil microhabitats over a period of active root growth for an annual grass. We found that the presence of root litter and living roots significantly altered gene expression involved in multiple nitrogen pathways, as well as tradeoffs between pathways, which ultimately regulate N availability to plants.

Longitudinal Microbiome Changes in Supragingival Biofilm Transcriptomes Induced by Orthodontics

INTRODUCTION

Common oral diseases are known to be associated with dysbiotic shifts in the supragingival microbiome, yet most oral microbiome associations with clinical end points emanate from cross-sectional studies. Orthodontic treatment is an elective procedure that can be exploited to prospectively examine clinically relevant longitudinal changes in the composition and function of the supragingival microbiome.

METHODS

A longitudinal cohort study was conducted among 24 adolescent orthodontic patients who underwent saliva and plaque sampling and clinical examinations at time points: before fixed appliance bonding and at 1, 6, and 12 wk thereafter. Clinical indices included bleeding on probing (BOP), mean gingival index (GI), probing depths (PDs), and plaque index (PI). To study the biologically (i.e., transcriptionally) active microbial communities, RNA was extracted from plaque and saliva for RNA sequencing and microbiome bioinformatics analysis. Longitudinal changes in microbiome beta diversity were examined using PERMANOVA tests, and the relative abundance of microbial taxa was measured using Kruskal-Wallis tests, Wilcoxon rank-sum tests, and negative binomial and zero-inflated mixed models.

RESULTS

Clinical measures of oral health deteriorated over time-the proportion of sites with GI and PI ≥1 increased by over 70% between prebonding and 12 wk postbonding while the proportion of sites with PD ≥4 mm increased 2.5-fold. Streptococcus sanguinis, a health-associated species that antagonizes cariogenic pathogens, showed a lasting decrease in relative abundance during orthodontic treatment. Contrarily, caries- and periodontal disease-associated taxa, including Selenomonas sputigena, Leptotrichia wadei, and Lachnoanaerobaculum saburreum, increased in abundance after bonding. Relative abundances of Stomatobaculum longum and Mogibacterium diversum in prebonding saliva predicted elevated BOP 12 wk postbonding, whereas Neisseria subflava was associated with lower BOP.

CONCLUSIONS

This study offers insights into longitudinal community and species-specific changes in the supragingival microbiome transcriptome during fixed orthodontic treatment, advancing our understanding of microbial dysbioses and identifying targets of future health-promoting clinical investigations.

KNOWLEDGE TRANSFER STATEMENT

Bonding braces was associated with subsequent changes in the oral microbiome characterized by increases in disease-associated species, decreases in health-associated species, and worsened clinical measures of oral health.

On the virome's hidden diversity: lessons from RHDV

Emerging infectious diseases are a major challenge to human and animal health. While predicting the emergence of pathogens is complex, the advent of high-throughput sequencing technologies has allowed the rapid identification of unknown microbiology diversity within organisms. Here, we discuss an example of a metatranscriptomics output to decipher viral evolution.

A Coculture of Photoautotrophs and Hydrolytic Heterotrophs Enables Efficient Upcycling of Starch from Wastewater toward Biomass-Derived Products: Synergistic Interactions Impacting Metabolism of the Consortium

Even with particular interest in sustainable development, due to the limited types of bioavailable carbon sources that could support heterotrophic/mixotrophic growth, microalgae-derived products still suffer from inconsistent yield and high costs. This study demonstrates a successful cocultivation of the photoautotroph Chlorella vulgaris with a hydrolytic-enzyme-abundant heterotroph, Saccharomycopsis fibuligera, enabling efficient starch upcycling from water/wastewater toward enhancing microalgae-dominant biomass and lipid production. The enzymatic activities of S. fibuligera contributed to the hydrolysis of starch into glucose, generating a 7-fold higher biomass through mixotrophic/heterotrophic growth of C. vulgaris. Further, scanning transmission electron microscopy (STEM) and quantitative analysis suggested a significantly induced accumulation of lipids in C. vulgaris. Results of meta-transcriptomics revealed the critical regulatory role of illumination in interaction shifting. Gene expression for glycolysis and lipid biosynthesis of C. vulgaris were highly activated during dark periods. Meanwhile, during illumination periods, genes coding for glucoamylase and the sulfur-related activities in S. fibuligera were significantly upregulated, leading to induced starch hydrolysis and potential increased competition for sulfur utilization, respectively. This study indicates that hydrolytic organisms could collaborate to make starch bioavailable for nonhydrolytic microalgae, thus broadening the substrate spectrum and making starch a novel biotechnological feedstock for microalgae-derived products, e.g., biofuels or single-cell protein.

Divergent hepaciviruses, delta-like viruses, and a chu-like virus in Australian marsupial carnivores (dasyurids)

Although Australian marsupials are characterised by unique biology and geographic isolation, little is known about the viruses present in these iconic wildlife species. The Dasyuromorphia are an order of marsupial carnivores found only in Australia that include both the extinct Tasmanian tiger (thylacine) and the highly threatened Tasmanian devil. Several other members of the order are similarly under threat of extinction due to habitat loss, hunting, disease, and competition and predation by introduced species such as feral cats. We utilised publicly available RNA-seq data from the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) database to document the viral diversity within four Dasyuromorph species. Accordingly, we identified fifteen novel virus sequences from five DNA virus families (Adenoviridae, Anelloviridae, Gammaherpesvirinae, Papillomaviridae, and Polyomaviridae) and three RNA virus taxa: the order Jingchuvirales, the genus Hepacivirus, and the delta-like virus group. Of particular note was the identification of a marsupial-specific clade of delta-like viruses that may indicate an association of deltaviruses with marsupial species. In addition, we identified a highly divergent hepacivirus in a numbat liver transcriptome that falls outside of the larger mammalian clade. We also detect what may be the first Jingchuvirales virus in a mammalian host-a chu-like virus in Tasmanian devils-thereby expanding the host range beyond invertebrates and ectothermic vertebrates. As many of these Dasyuromorphia species are currently being used in translocation efforts to reseed populations across Australia, understanding their virome is of key importance to prevent the spread of viruses to naive populations.

Inoculation with Lentilactobacillus buchneri alone or in combination with Lentilactobacillus hilgardii modifies gene expression, fermentation profile, and starch digestibility in high-moisture corn

Inoculants combining Lentilactobacillus buchneri and Lentilactobacillus hilgardii have been shown to improve the aerobic stability of high-moisture corn (HMC) and whole-plant corn silage, but the mode of action of this co-inoculation remains to be elucidated. This study used metatranscriptomics to evaluate the effects of inoculation with L. buchneri alone or combined with L. hilgardii on the bacterial community, gene expression, fermentation profile, and starch digestibility in HMC. High-moisture corn not inoculated (Control) or inoculated with L. buchneri NCIMB 40788 (LB) or L. buchneri NCIMB 40788 combined with L. hilgardii CNCM-I-4785 (Combo) was ensiled in mini silo bags for 30, 60, 120, and 180 days. The fermentation profile was evaluated at all time points. Metatranscriptomics was performed on samples collected on day 120. Combo had a greater alpha diversity richness index of contigs than LB and Control, and inoculation with Combo and LB modified the beta-diversity of contigs compared to Control. Out of 69 genes of interest, 20 were differentially expressed in LB compared to Control and 25 in Combo compared to Control. Of those differently expressed genes, 16 (10 of which were associated with carbohydrate metabolism and six with amino acid metabolism) were differently expressed in both LB and Combo compared to Control, and all those genes were upregulated in the inoculated silages. When we compared Combo and LB, we found seven genes expressed differently, four associated with carbohydrate metabolism and downregulated in Combo, and three associated with amino acid metabolism and upregulated in Combo. At day 120, the inoculated silages had more culturable lactic acid bacteria, higher Lactobacillus relative abundance, and lower Leuconostoc relative abundance than Control. The concentration of acetic acid remained low throughout ensiling in Control, but in LB and Combo, it increased up to day 60 and remained stable from day 60 to 180. The 1,2-propanediol was only detected in LB and Combo. Inoculation did not affect the concentration of starch, but starch digestibility was greater in Combo than in Control. Inoculation of HMC with Combo modified the gene expression and fermentation profile compared to Control and LB, improving starch digestibility compared to uninoculated HMC.

Insights into bacterial community metatranscriptome and metabolome in river water influenced by palm oil mill effluent final discharge

AIMS

This study aimed to investigate the effect of palm oil mill effluent (POME) final discharge on the active bacterial composition, gene expression, and metabolite profiles in the receiving rivers to establish a foundation for identifying potential biomarkers for monitoring POME pollution in rivers.

METHODS AND RESULTS

The POME final discharge, upstream (unpolluted by POME), and downstream (effluent receiving point) parts of the rivers from two sites were physicochemically characterized. The taxonomic and gene profiles were then evaluated using de novo metatranscriptomics, while the metabolites were detected using qualitative metabolomics. A similar bacterial community structure in the POME final discharge samples from both sites was recorded, but their composition varied. Redundancy analysis showed that several families, particularly Comamonadaceae and Burkholderiaceae [Pr(>F) = 0.028], were positively correlated with biochemical oxygen demand (BOD5) and chemical oxygen demand (COD). The results also showed significant enrichment of genes regulating various metabolisms in the POME-receiving rivers, with methane, carbon fixation pathway, and amino acids among the predominant metabolisms identified (FDR < 0.05, PostFC > 4, and PPDE > 0.95). This was further validated through qualitative metabolomics, whereby amino acids were detected as the predominant metabolites.

CONCLUSIONS

The results suggest that genes regulating amino acid metabolism have significant potential for developing effective biomonitoring and bioremediation strategies in river water influenced by POME final discharge, fostering a sustainable palm oil industry.

Virus discovery in dogs with non-suppurative encephalitis reveals a high incidence of tick-borne encephalitis virus infections in Switzerland

INTRODUCTION

Viral infections are a frequent cause of disseminated non-suppurative encephalitis in dogs. However, using routine diagnostic methods, the specific virus may remain unknown due to extensive or complete viral clearance or because the virus is unexpected or new. A metatranscriptomics-based approach of combining high-throughput sequencing (HTS) and bioinformatics analysis was used to investigate the viral etiology in archival cases of dogs with non-suppurative encephalitis. In formalin-fixed paraffin embedded (FFPE) brain material from the years 1976 to 2021 a high incidence of tick-borne encephalitis virus (TBEV) was detected. Moreover, canine distemper virus (CDV) was identified without typical demyelinating lesions and canine vesivirus (CaVV) was detected as an unexpected virus associated with non-suppurative encephalitis. We demonstrated the viral presence in brain tissues at the sites of inflammation by immunohistochemistry (IHC) and in situ hybridization (ISH). These results highlight the value of emerging sequencing technologies in veterinary diagnostics and expand our knowledge on the etiologies of encephalitis in dogs.

Into the microbial niche

The environmental niche concept describes the distribution of a taxon in the environment and can be used to understand community dynamics, biological invasions, and the impact of environmental changes. The uses and applications are still restricted in microbial ecology, largely due to the complexity of microbial systems and associated methodological limitations. The development of shotgun metagenomics and metatranscriptomics opens new ways to investigate the microbial niche by focusing on the metabolic niche within the environmental space. Here, we propose the metabolic niche framework, which, by defining the fundamental and realised metabolic niche of microorganisms, has the potential to not only provide novel insights into habitat preferences and the metabolism associated, but also to inform on metabolic plasticity, niche shifts, and microbial invasions.

Comprehensive evaluation of methods for differential expression analysis of metatranscriptomics data

Understanding the function of the human microbiome is important but the development of statistical methods specifically for the microbial gene expression (i.e. metatranscriptomics) is in its infancy. Many currently employed differential expression analysis methods have been designed for different data types and have not been evaluated in metatranscriptomics settings. To address this gap, we undertook a comprehensive evaluation and benchmarking of 10 differential analysis methods for metatranscriptomics data. We used a combination of real and simulated data to evaluate performance (i.e. type I error, false discovery rate and sensitivity) of the following methods: log-normal (LN), logistic-beta (LB), MAST, DESeq2, metagenomeSeq, ANCOM-BC, LEfSe, ALDEx2, Kruskal-Wallis and two-part Kruskal-Wallis. The simulation was informed by supragingival biofilm microbiome data from 300 preschool-age children enrolled in a study of childhood dental disease (early childhood caries, ECC), whereas validations were sought in two additional datasets from the ECC study and an inflammatory bowel disease study. The LB test showed the highest sensitivity in both small and large samples and reasonably controlled type I error. Contrarily, MAST was hampered by inflated type I error. Upon application of the LN and LB tests in the ECC study, we found that genes C8PHV7 and C8PEV7, harbored by the lactate-producing Campylobacter gracilis, had the strongest association with childhood dental disease. This comprehensive model evaluation offers practical guidance for selection of appropriate methods for rigorous analyses of differential expression in metatranscriptomics. Selection of an optimal method increases the possibility of detecting true signals while minimizing the chance of claiming false ones.

Evolutionary conservation of Trichomonas-mycoplasma symbiosis across the host species barrier

Introduction

The protozoan parasite Trichomonas vaginalis is the most common cellular sexually transmitted disease in humans, and the closely related species Trichomonas gallinae is an avian parasite of ecological and economic importance. Phylogenetic evidence suggests T. vaginalis arose during bird to human transmission of a T. gallinae-like ancestor. Trichomonas vaginalis shares a strong clinical association with the independent sexually transmitted pathogen Metamycoplasma (formerly Mycoplasma) hominis, and the uncultured bacterium "Candidatus Malacoplasma (formerly Mycoplasma) girerdii," with the latter association being an order of magnitude stronger. Both bacterial species have been shown to profoundly influence T. vaginalis growth, energy production and virulence-associated mechanisms.

Methods

Evidence for a novel Malacoplasma sp. was discovered by in vivo Illumina metatranscriptomics sequencing of the T. gallinae-infected pigeon mouth. We leveraged published 16S rDNA profiling data from digestive tract of 12 healthy and 24 T. gallinae-infected pigeons to investigate association between the novel Malacoplasma sp. and T. gallinae. We utilised Illumina metagenomics sequencing targeted to pigeon oral and crop samples infected with the novel Malacoplasma sp. to generate its full-length genome sequence. Sequence similarity network analysis was used to compare annotated proteins from the novel Malacoplasma sp. with a range of other related species.

Results

Here we present evidence for a novel Malacoplasma species, related to "Ca. M. girerdii," that is strongly associated with T. gallinae in the upper digestive tract of domestic pigeons. Analysis of the genome sequence revealed gene features apparently specific to a Trichomonas-symbiotic Malacoplasma lineage.

Discussion

These data support a model of long-term association between Trichomonas and Malacoplasma spp. that has been conserved across diversification of the Trichomonas lineage and the host species barrier from birds to human.

The antibiotic resistance reservoir of the lung microbiome expands with age

Antimicrobial resistant lower respiratory tract infections (LRTI) are an increasing public health threat, and an important cause of global mortality. The lung microbiome influences LRTI susceptibility and represents an important reservoir for exchange of antimicrobial resistance genes (ARGs). Studies of the gut microbiome have found an association between age and increasing antimicrobial resistance gene (ARG) burden, however corollary studies in the lung microbiome remain absent, despite the respiratory tract representing one of the most clinically significant sites for drug resistant infections. We performed a prospective, multicenter observational study of 261 children and 88 adults with acute respiratory failure, ranging in age from 31 days to ≥ 89 years, admitted to intensive care units in the United States. We performed RNA sequencing on tracheal aspirates collected within 72 hours of intubation, and evaluated age-related differences in detectable ARG expression in the lung microbiome as a primary outcome. Secondary outcomes included number and classes of ARGs detected, proportion of patients with an ARG class, and composition of the lung microbiome. Multivariable logistic regression models (adults vs children) or continuous age (years) were adjusted for sex, race/ethnicity, LRTI status, and days from intubation to specimen collection. Detection of ARGs was significantly higher in adults compared with children after adjusting for sex, race/ethnicity, LRTI diagnosis, and days from intubation to specimen collection (adjusted odds ratio (aOR): 2.16, 95% confidence interval (CI): 1.10-4.22). A greater proportion of adults compared with children had beta-lactam ARGs (31% (CI: 21-41%) vs 13% (CI: 10-18%)), aminoglycoside ARGs (20% (CI: 13-30%) vs 2% (CI: 0.6-4%)), and tetracycline ARGs (14% (CI: 7-23%) vs 3% (CI: 1-5%)). Adults ≥70 years old had the highest proportion of these three ARG classes. The total bacterial abundance of the lung microbiome increased with age, and microbiome alpha diversity varied with age. Taxonomic composition of the lung microbiome, measured by Bray Curtis dissimilarity index, differed between adults and children (p = 0.003). The association between age and increased ARG detection remained significant after additionally including lung microbiome total bacterial abundance and alpha diversity in the multivariable logistic regression model (aOR: 2.38, (CI: 1.25-4.54)). Furthermore, this association remained robust when modeling age as a continuous variable (aOR: 1.02, (CI: 1.01-1.03) per year of age). Taken together, our results demonstrate that age is an independent risk factor for ARG detection in the lower respiratory tract microbiome. These data shape our understanding of the lung resistome in critically ill patients across the lifespan, which may have implications for clinical management and global public health.

Archaeome structure and function of the intestinal tract in healthy and H1N1 infected swine

Background

Methanogenic archaea represent a less investigated and likely underestimated part of the intestinal tract microbiome in swine.

Aims/Methods

This study aims to elucidate the archaeome structure and function in the porcine intestinal tract of healthy and H1N1 infected swine. We performed multi-omics analysis consisting of 16S rRNA gene profiling, metatranscriptomics and metaproteomics.

Results and discussion

We observed a significant increase from 0.48 to 4.50% of archaea in the intestinal tract microbiome along the ileum and colon, dominated by genera Methanobrevibacter and Methanosphaera. Furthermore, in feces of naïve and H1N1 infected swine, we observed significant but minor differences in the occurrence of archaeal phylotypes over the course of an infection experiment. Metatranscriptomic analysis of archaeal mRNAs revealed the major methanogenesis pathways of Methanobrevibacter and Methanosphaera to be hydrogenotrophic and methyl-reducing, respectively. Metaproteomics of archaeal peptides indicated some effects of the H1N1 infection on central metabolism of the gut archaea.

Conclusions/Take home message

Finally, this study provides the first multi-omics analysis and high-resolution insights into the structure and function of the porcine intestinal tract archaeome during a non-lethal Influenza A virus infection of the respiratory tract, demonstrating significant alterations in archaeal community composition and central metabolic functions.

Consistency across multi-omics layers in a drug-perturbed gut microbial community

Multi-omics analyses are used in microbiome studies to understand molecular changes in microbial communities exposed to different conditions. However, it is not always clear how much each omics data type contributes to our understanding and whether they are concordant with each other. Here, we map the molecular response of a synthetic community of 32 human gut bacteria to three non-antibiotic drugs by using five omics layers (16S rRNA gene profiling, metagenomics, metatranscriptomics, metaproteomics and metabolomics). We find that all the omics methods with species resolution are highly consistent in estimating relative species abundances. Furthermore, different omics methods complement each other for capturing functional changes. For example, while nearly all the omics data types captured that the antipsychotic drug chlorpromazine selectively inhibits Bacteroidota representatives in the community, the metatranscriptome and metaproteome suggested that the drug induces stress responses related to protein quality control. Metabolomics revealed a decrease in oligosaccharide uptake, likely caused by Bacteroidota depletion. Our study highlights how multi-omics datasets can be utilized to reveal complex molecular responses to external perturbations in microbial communities.

Microbiomes, Their Function, and Cancer: How Metatranscriptomics Can Close the Knowledge Gap

The interaction between the microbial communities in the human body and the onset and progression of cancer has not been investigated until recently. The vast majority of the metagenomics research in this area has concentrated on the composition of microbiomes, attempting to link the overabundance or depletion of certain microorganisms to cancer proliferation, metastatic behaviour, and its resistance to therapies. However, studies elucidating the functional implications of the microbiome activity in cancer patients are still scarce; in particular, there is an overwhelming lack of studies assessing such implications directly, through analysis of the transcriptome of the bacterial community. This review summarises the contributions of metagenomics and metatranscriptomics to the knowledge of the microbial environment associated with several cancers; most importantly, it highlights all the advantages that metatranscriptomics has over metagenomics and suggests how such an approach can be leveraged to advance the knowledge of the cancer bacterial environment.