Detection of large numbers of novel sequences in the metatranscriptomes of complex marine microbial communities

Decoding the molecular concerto: Toxicotranscriptomic evaluation of microplastic and nanoplastic impacts on aquatic organisms

The pervasive and steadily increasing presence of microplastics/nanoplastics (MPs/NPs) in aquatic environments has raised significant concerns regarding their potential adverse effects on aquatic organisms and their integration into trophic dynamics. This emerging issue has garnered the attention of (eco)toxicologists, promoting the utilization of toxicotranscriptomics to unravel the responses of aquatic organisms not only to MPs/NPs but also to a wide spectrum of environmental pollutants. This review aims to systematically explore the broad repertoire of predicted molecular responses by aquatic organisms, providing valuable intuitions into complex interactions between plastic pollutants and aquatic biota. By synthesizing the latest literature, present analysis sheds light on transcriptomic signatures like gene expression, interconnected pathways and overall molecular mechanisms influenced by various plasticizers. Harmful effects of these contaminants on key genes/protein transcripts associated with crucial pathways lead to abnormal immune response, metabolic response, neural response, apoptosis and DNA damage, growth, development, reproductive abnormalities, detoxification, and oxidative stress in aquatic organisms. However, unique challenge lies in enhancing the fingerprint of MPs/NPs, presenting complicated enigma that requires decoding their specific impact at molecular levels. The exploration endeavors, not only to consolidate existing knowledge, but also to identify critical gaps in understanding, push forward the frontiers of knowledge about transcriptomic signatures of plastic contaminants. Moreover, this appraisal emphasizes the imperative to monitor and mitigate the contamination of commercially important aquatic species by MPs/NPs, highlighting the pivotal role that regulatory frameworks must play in protecting all aquatic ecosystems. This commitment aligns with the broader goal of ensuring the sustainability of aquatic resources and the resilience of ecosystems facing the growing threat of plastic pollutants.

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.

Freshwater plastisphere: a review on biodiversity, risks, and biodegradation potential with implications for the aquatic ecosystem health

The plastisphere, a unique microbial biofilm community colonizing plastic debris and microplastics (MPs) in aquatic environments, has attracted increasing attention owing to its ecological and public health implications. This review consolidates current state of knowledge on freshwater plastisphere, focussing on its biodiversity, community assembly, and interactions with environmental factors. Current biomolecular approaches revealed a variety of prokaryotic and eukaryotic taxa associated with plastic surfaces. Despite their ecological importance, the presence of potentially pathogenic bacteria and mobile genetic elements (i.e., antibiotic resistance genes) raises concerns for ecosystem and human health. However, the extent of these risks and their implications remain unclear. Advanced sequencing technologies are promising for elucidating the functions of plastisphere, particularly in plastic biodegradation processes. Overall, this review emphasizes the need for comprehensive studies to understand plastisphere dynamics in freshwater and to support effective management strategies to mitigate the impact of plastic pollution on freshwater resources.

Environmental transcriptomics under heat stress: Can environmental RNA reveal changes in gene expression of aquatic organisms?

To safeguard biodiversity in a changing climate, taxonomic information about species turnover and insights into the health of organisms are required. Environmental DNA approaches are increasingly used for species identification, but cannot provide functional insights. Transcriptomic methods reveal the physiological states of macroorganisms, but are currently species-specific and require tissue sampling or animal sacrifice, making community-wide assessments challenging. Here, we test whether broad functional information (expression level of the transcribed genes) can be harnessed from environmental RNA (eRNA), which includes extra-organismal RNA from macroorganisms along with whole microorganisms. We exposed Daphnia pulex as well as phytoplankton prey and microorganism colonizers to control (20°C) and heat stress (28°C) conditions for 7 days. We sequenced eRNA from tank water (after complete removal of Daphnia) as well as RNA from Daphnia tissue, enabling comparisons of extra-organismal and organismal RNA-based gene expression profiles. Both RNA types detected similar heat stress responses of Daphnia. Using eRNA, we identified 32 Daphnia genes to be differentially expressed following heat stress. Of these, 17 were also differentially expressed and exhibited similar levels of relative expression in organismal RNA. In addition to the extra-organismal Daphnia response, eRNA detected community-wide heat stress responses consisting of distinct functional profiles and 121 differentially expressed genes across eight taxa. Our study demonstrates that environmental transcriptomics based on extra-organismal eRNA can noninvasively reveal gene expression responses of macroorganisms following environmental changes, with broad potential implications for the biomonitoring of health across the trophic chain.

Using metatranscriptomics to better understand the role of microbial nitrogen cycling in coastal sediment benthic flux denitrification efficiency

Spatial and temporal variability in benthic flux denitrification efficiency occurs across Port Phillip Bay, Australia. Here, we assess the capacity for untargeted metatranscriptomics to resolve spatiotemporal differences in the microbial contribution to benthic nitrogen cycling. The most abundant sediment transcripts assembled were associated with the archaeal nitrifier Nitrosopumilus. In sediments close to external inputs of organic nitrogen, the dominant transcripts were associated with Nitrosopumilus nitric oxide nitrite reduction (nirK). The environmental conditions close to organic nitrogen inputs that select for increased transcription in Nitrosopumilus (amoCAB, nirK, nirS, nmo, hcp) additionally selected for increased transcription of bacterial nitrite reduction (nxrB) and transcripts associated with anammox (hzo) but not denitrification (bacterial nirS/nirk). In sediments that are more isolated from external inputs of organic nitrogen dominant transcripts were associated with nitrous oxide reduction (nosZ) and changes in nosZ transcript abundance were uncoupled from transcriptional profiles associated with archaeal nitrification. Coordinated transcription of coupled community-level nitrification-denitrification was not well supported by metatranscriptomics. In comparison, the abundance of archaeal nirK transcripts were site- and season-specific. This study indicates that the transcription of archaeal nirK in response to changing environmental conditions may be an important and overlooked feature of coastal sediment nitrogen cycling.

Reverse engineering environmental metatranscriptomes clarifies best practices for eukaryotic assembly

BACKGROUND

Diverse communities of microbial eukaryotes in the global ocean provide a variety of essential ecosystem services, from primary production and carbon flow through trophic transfer to cooperation via symbioses. Increasingly, these communities are being understood through the lens of omics tools, which enable high-throughput processing of diverse communities. Metatranscriptomics offers an understanding of near real-time gene expression in microbial eukaryotic communities, providing a window into community metabolic activity.

RESULTS

Here we present a workflow for eukaryotic metatranscriptome assembly, and validate the ability of the pipeline to recapitulate real and manufactured eukaryotic community-level expression data. We also include an open-source tool for simulating environmental metatranscriptomes for testing and validation purposes. We reanalyze previously published metatranscriptomic datasets using our metatranscriptome analysis approach.

CONCLUSION

We determined that a multi-assembler approach improves eukaryotic metatranscriptome assembly based on recapitulated taxonomic and functional annotations from an in-silico mock community. The systematic validation of metatranscriptome assembly and annotation methods provided here is a necessary step to assess the fidelity of our community composition measurements and functional content assignments from eukaryotic metatranscriptomes.

The Epidermal Microbiome Within an Aggregation of Leopard Sharks (Triakis semifasciata) Has Taxonomic Flexibility with Gene Functional Stability Across Three Time-points

The epidermis of Chondrichthyan fishes consists of dermal denticles with production of minimal but protein-rich mucus that collectively, influence the attachment and biofilm development of microbes, facilitating a unique epidermal microbiome. Here, we use metagenomics to provide the taxonomic and functional characterization of the epidermal microbiome of the Triakis semifasciata (leopard shark) at three time-points collected across 4 years to identify links between microbial groups and host metabolism. Our aims include (1) describing the variation of microbiome taxa over time and identifying recurrent microbiome members (present across all time-points); (2) investigating the relationship between the recurrent and flexible taxa (those which are not found consistently across time-points); (3) describing the functional compositions of the microbiome which may suggest links with the host metabolism; and (4) identifying whether metabolic processes are shared across microbial genera or are unique to specific taxa. Microbial members of the microbiome showed high similarity between all individuals (Bray-Curtis similarity index = 82.7, where 0 = no overlap, 100 = total overlap) with the relative abundance of those members varying across sampling time-points, suggesting flexibility of taxa in the microbiome. One hundred and eighty-eight genera were identified as recurrent, including Pseudomonas, Erythrobacter, Alcanivorax, Marinobacter, and Sphingopxis being consistently abundant across time-points, while Limnobacter and Xyella exhibited switching patterns with high relative abundance in 2013, Sphingobium and Sphingomona in 2015, and Altermonas, Leeuwenhoekiella, Gramella, and Maribacter in 2017. Of the 188 genera identified as recurrent, the top 19 relatively abundant genera formed three recurrent groups. The microbiome also displayed high functional similarity between individuals (Bray-Curtis similarity index = 97.6) with gene function composition remaining consistent across all time-points. These results show that while the presence of microbial genera exhibits consistency across time-points, their abundances do fluctuate. Microbial functions however remain stable across time-points; thus, we suggest the leopard shark microbiomes exhibit functional redundancy. We show coexistence of microbes hosted in elasmobranch microbiomes that encode genes involved in utilizing nitrogen, but not fixing nitrogen, degrading urea, and resistant to heavy metal.

The oral microbiome in the pathophysiology of cardiovascular disease

Despite advances in our understanding of the pathophysiology of many cardiovascular diseases (CVDs) and expansion of available therapies, the global burden of CVD-associated morbidity and mortality remains unacceptably high. Important gaps remain in our understanding of the mechanisms of CVD and determinants of disease progression. In the past decade, much research has been conducted on the human microbiome and its potential role in modulating CVD. With the advent of high-throughput technologies and multiomics analyses, the complex and dynamic relationship between the microbiota, their 'theatre of activity' and the host is gradually being elucidated. The relationship between the gut microbiome and CVD is well established. Much less is known about the role of disruption (dysbiosis) of the oral microbiome; however, interest in the field is growing, as is the body of literature from basic science and animal and human investigations. In this Review, we examine the link between the oral microbiome and CVD, specifically coronary artery disease, stroke, peripheral artery disease, heart failure, infective endocarditis and rheumatic heart disease. We discuss the various mechanisms by which oral dysbiosis contributes to CVD pathogenesis and potential strategies for prevention and treatment.

Exploring Newer Biosynthetic Gene Clusters in Marine Microbial Prospecting

Marine microbes genetically evolved to survive varying salinity, temperature, pH, and other stress factors by producing different bioactive metabolites. These microbial secondary metabolites (SMs) are novel, have high potential, and could be used as lead molecule. Genome sequencing of microbes revealed that they have the capability to produce numerous novel bioactive metabolites than observed under standard in vitro culture conditions. Microbial genome has specific regions responsible for SM assembly, termed biosynthetic gene clusters (BGCs), possessing all the necessary genes to encode different enzymes required to generate SM. In order to augment the microbial chemo diversity and to activate these gene clusters, various tools and techniques are developed. Metagenomics with functional gene expression studies aids in classifying novel peptides and enzymes and also in understanding the biosynthetic pathways. Genome shuffling is a high-throughput screening approach to improve the development of SMs by incorporating genomic recombination. Transcriptionally silent or lower level BGCs can be triggered by artificially knocking promoter of target BGC. Additionally, bioinformatic tools like antiSMASH, ClustScan, NAPDOS, and ClusterFinder are effective in identifying BGCs of existing class for annotation in genomes. This review summarizes the significance of BGCs and the different approaches for detecting and elucidating BGCs from marine microbes.

Global Distribution and Richness of Armillaria and Related Species Inferred From Public Databases and Amplicon Sequencing Datasets

Armillaria is a globally distributed fungal genus most notably composed of economically important plant pathogens that are found predominantly in forest and agronomic systems. The genus sensu lato has more recently received attention for its role in woody plant decomposition and in mycorrhizal symbiosis with specific plants. Previous phylogenetic analyses suggest that around 50 species are recognized globally. Despite this previous work, no studies have analyzed the global species richness and distribution of the genus using data derived from fungal community sequencing datasets or barcoding initiatives. To assess the global diversity and species richness of Armillaria, we mined publicly available sequencing datasets derived from numerous primer regions for the ribosomal operon, as well as ITS sequences deposited on Genbank, and clustered them akin to metabarcoding studies. Our estimates reveal that species richness ranges from 50 to 60 species, depending on whether the ITS1 or ITS2 marker is used. Eastern Asia represents the biogeographic region with the highest species richness. We also assess the overlap of species across geographic regions and propose some hypotheses regarding the drivers of variability in species diversity and richness between different biogeographic regions.

Integrated metagenomic and metatranscriptomic analyses of ultraviolet disinfection effects on antibiotic resistance genes and bacterial communities during wastewater treatment

Ultraviolet (UV) disinfection is now widely implemented in wastewater treatment plants (WWTPs) worldwide, but its effect on antibiotic resistome of the surviving bacteria remains unclear. In this study, we employed high-throughput sequencing-based metagenomic and metatranscriptomic approaches to comprehensively elucidate the effects of UV disinfection on the shifts of bacterial community and antibiotic resistance genes (ARGs) on both DNA and mRNA levels in one WWTP. Metagenomic analyses revealed an insignificant change in the bacterial community after UV disinfection, while metatranscriptomic analyses showed that UV disinfection significantly changed the abundance of 13.79% of phyla and 10.32% of genera. In total, 38 ARG-like open reading frames (ORFs) and 327 ARG-like transcripts were identified in the DNA and RNA samples, respectively. The relative abundances of the total ARGs, each ARG type, and each ARG subtype also varied after UV disinfection. Additionally, UV disinfection significantly reduced the expression of total ARGs from 49.40 transcripts per kilobase of exon model per million mapped reads (TPM) to 47.62 TPM, and significantly changed the expression of 10.75% of ARG subtypes in wastewater (p < 0.05). Notably, the significant increase in the expression and obvious increase in the relative abundance of macrolide-lincosamide-streptogramin B (MLSB) resistance genes revealed that UV disinfection increases the potential health risk of MLSB resistance genes in wastewater. Moreover, potential host analyses of ARGs revealed the different preferences of antibiotic resistant bacteria (ARB) to ARGs. This study may shed new light on the underlying mechanism of the UV disinfection effect on antibiotic resistance.

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.

Lab-on-a-chip technology for in situ combined observations in oceanography

The oceans sustain the global environment and diverse ecosystems through a variety of biogeochemical processes and their complex interactions. In order to understand the dynamism of the local or global marine environments, multimodal combined observations must be carried out in situ. On the other hand, instrumentation of in situ measurement techniques enabling biological and/or biochemical combined observations is challenging in aquatic environments, including the ocean, because biochemical flow analyses require a more complex configuration than physicochemical electrode sensors. Despite this technical hurdle, in situ analyzers have been developed to measure the concentrations of seawater contents such as nutrients, trace metals, and biological components. These technologies have been used for cutting-edge ocean observations to elucidate the biogeochemical properties of water mass with a high spatiotemporal resolution. In this context, the contribution of lab-on-a-chip (LoC) technology toward the miniaturization and functional integration of in situ analyzers has been gaining momentum. Due to their mountability, in situ LoC technologies provide ideal instrumentation for underwater analyzers, especially for miniaturized underwater observation platforms. Consequently, the appropriate combination of reliable LoC and underwater technologies is essential to realize practical in situ LoC analyzers suitable for underwater environments, including the deep sea. Moreover, the development of fundamental LoC technologies for underwater analyzers, which operate stably in extreme environments, should also contribute to in situ measurements for public or industrial purposes in harsh environments as well as the exploration of the extraterrestrial frontier.

Quantitative models of nitrogen-fixing organisms

Nitrogen-fixing organisms are of importance to the environment, providing bioavailable nitrogen to the biosphere. Quantitative models have been used to complement the laboratory experiments and in situ measurements, where such evaluations are difficult or costly. Here, we review the current state of the quantitative modeling of nitrogen-fixing organisms and ways to enhance the bridge between theoretical and empirical studies.

FoodOmics as a new frontier to reveal microbial community and metabolic processes occurring on table olives fermentation

Table olives are considered the most widespread fermented food in the Mediterranean area and their consumption is expanding all over the world. This fermented vegetable can be considered as a natural functional food thanks to their high nutritional value and high content of bioactive compounds that contribute to the health and well-being of consumers. The presence of bioactive compounds is strongly influenced by a complex microbial consortium, traditionally exploited through culture-dependent approaches. Recently, the rapid spread of omics technologies has represented an important challenge to better understand the function, the adaptation and the exploitation of microbial diversity in different complex ecosystems, such as table olives. This review provides an overview of the potentiality of omics technologies to in depth investigate the microbial composition and the metabolic processes that drive the table olives fermentation, affecting both sensorial profile and safety properties of the final product. Finally, the review points out the role of omics approaches to raise at higher sophisticated level the investigations on microbial, gene, protein, and metabolite, with huge potential for the integration of table olives composition with functional assessments.

Phenome-Genome Profiling of Single Bacterial Cell by Raman-Activated Gravity-Driven Encapsulation and Sequencing

The small size and low DNA amount of bacterial cells have hindered establishing phenome-genome links in a precisely indexed, one-cell-per-reaction manner. Here, Raman-Activated Gravity-driven single-cell Encapsulation and Sequencing (RAGE-Seq) is presented, where individual cells are phenotypically screened via single-cell Raman spectra (SCRS) in an aquatic, vitality-preserving environment, then the cell with targeted SCRS is precisely packaged in a picoliter microdroplet and readily exported in a precisely indexed, "one-cell-one-tube" manner. Such integration of microdroplet encapsulation to Raman-activated sorting ensures high-coverage one-cell genome sequencing or cultivation that is directly linked to metabolic phenotype. For clinical Escherichia coli isolates, genome assemblies derived from precisely one cell via RAGE-Seq consistently reach >95% coverage. Moreover, directly from a urine sample of urogenital tract infection, metabolic-activity-based antimicrobial susceptibility phenotypes and genome sequence of 99.5% coverage are obtained simultaneously from precisely one cell. This single-cell global mutation map corroborates resistance phenotype and genotype, and unveils epidemiological features with high specificity and sensitivity. The ability to profile and correlate bacterial metabolic phenome and high-quality genome sequences at one-cell resolution suggests broad application of RAGE-Seq.

Copepod-Associated Gammaproteobacterial Alkaline Phosphatases in the North Atlantic Subtropical Gyre

Planktonic organisms may provide a niche to associated bacteria in the oligotrophic ocean. Bacterial fitness strategies in association with copepods - abundant planktonic crustaceans - were examined by sampling and incubation experiments in the North Atlantic Subtropical Gyre (NASG). The bacterial metatranscriptome was dominated by Gammaproteobacteria and showed expression of complete bacterial pathways including chemotaxis, cell signaling, and alkaline phosphatase activity. Quantitative PCR and reverse transcriptase qPCR revealed the consistent presence and expression of alkaline phosphatase genes primarily by Vibrio spp. in the copepod association. Copepod-associated bacteria appear to respond to prevailing phosphorus limitation by using alkaline phosphatases to break down organophosphoesters, presumably originating from the copepods. The results suggest that the basin-wide tendency for phosphorus limitation in the North Atlantic Ocean is occurring at microscales in these nitrogen-enriched copepod microenvironments. The bacterial communities and their fitness strategies supported by associations with these abundant mesozooplankton are unique from the surrounding seawater and could have large-scale implications for biogeochemical cycling, marine food web structuring, and copepod and ecosystem health.

Antimicrobial Resistance in Veterinary Medicine: An Overview

Antimicrobial resistance (AMR) represents one of the most important human- and animal health-threatening issues worldwide. Bacterial capability to face antimicrobial compounds is an ancient feature, enabling bacterial survival over time and the dynamic surrounding. Moreover, bacteria make use of their evolutionary machinery to adapt to the selective pressure exerted by antibiotic treatments, resulting in reduced efficacy of the therapeutic intervention against human and animal infections. The mechanisms responsible for both innate and acquired AMR are thoroughly investigated. Commonly, AMR traits are included in mobilizable genetic elements enabling the homogeneous diffusion of the AMR traits pool between the ecosystems of diverse sectors, such as human medicine, veterinary medicine, and the environment. Thus, a coordinated multisectoral approach, such as One-Health, provides a detailed comprehensive picture of the AMR onset and diffusion. Following a general revision of the molecular mechanisms responsible for both innate and acquired AMR, the present manuscript focuses on reviewing the contribution of veterinary medicine to the overall issue of AMR. The main sources of AMR amenable to veterinary medicine are described, driving the attention towards the indissoluble cross-talk existing between the diverse ecosystems and sectors and their cumulative cooperation to this warning phenomenon.

Bacterial metagenome analysis of Mytilus galloprovincialis collected from Istanbul and Izmir coastal stations of Turkey

Mytilus galloprovincialis is a marine mollusk belonging to the Bivalvia class. It has been distributed largely in Turkish shores and worldwide aquatic environments. Besides being known as an environmental pollution indicator, it is highly consumed as a food and has a high economic value. Due to their nutritional mechanisms by filtering water, they are affected by pollution in seawater and mussels can host-microbial diversity of environmental origin as well as pathogenic bacteria. Therefore, in this study, bacterial species found in Mediterranean mussels collected from the coastal stations of Istanbul [Rumeli Kavagi (RK), Kucukcekmece (KC)], and Izmir [(Foca (MF), Urla (MU)] were investigated and compared with microbiological and metagenomic analyses. According to microbiological analysis results, 34 mussel-associated Enterobacteriaceae and Vibrionaceae family members were identified. As a result of the culture-independent metagenomic analysis, taxonomic groups for each station were identified and compared based on Operational Taxonomic Unit data. For all stations, the most abundant bacterial genera were the unclassified bacterial genera. The total number of mussel-related total richness identified in all groups was 4889 (RK = 1605; KC = 1930; MF = 1508; and MU = 1125). According to the metagenomic data obtained in this study, different relative amounts of Lachnospiraceae and Bacteroidetes taxa groups were reported for all stations. The pathogenic bacterial genera identified by metagenomic analyses which may be significant for the public health are Arcobacter, Clostridium, Aeromonas, Vibrio, Escherichia_Shigella, Klebsiella, Campylobacter, Helicobacter, Pseudomonas, Morganella, Serratia, Corynebacterium, Enterococcus, Staphylococcus, Yersinia, Mycoplasma, Brucellaceae_unclassified, Pantoea, and Proteus.

Microbial resolution of whole genome shotgun and 16S amplicon metagenomic sequencing using publicly available NEON data

Microorganisms are ubiquitous in the biosphere, playing a crucial role in both biogeochemistry of the planet and human health. However, identifying these microorganisms and defining their function are challenging. Widely used approaches in comparative metagenomics, 16S amplicon sequencing and whole genome shotgun sequencing (WGS), have provided access to DNA sequencing analysis to identify microorganisms and evaluate diversity and abundance in various environments. However, advances in parallel high-throughput DNA sequencing in the past decade have introduced major hurdles, namely standardization of methods, data storage, reproducible interoperability of results, and data sharing. The National Ecological Observatory Network (NEON), established by the National Science Foundation, enables all researchers to address queries on a regional to continental scale around a variety of environmental challenges and provide high-quality, integrated, and standardized data from field sites across the U.S. As the amount of metagenomic data continues to grow, standardized procedures that allow results across projects to be assessed and compared is becoming increasingly important in the field of metagenomics. We demonstrate the feasibility of using publicly available NEON soil metagenomic sequencing datasets in combination with open access Metagenomics Rapid Annotation using the Subsystem Technology (MG-RAST) server to illustrate advantages of WGS compared to 16S amplicon sequencing. Four WGS and four 16S amplicon sequence datasets, from surface soil samples prepared by NEON investigators, were selected for comparison, using standardized protocols collected at the same locations in Colorado between April-July 2014. The dominant bacterial phyla detected across samples agreed between sequencing methodologies. However, WGS yielded greater microbial resolution, increased accuracy, and allowed identification of more genera of bacteria, archaea, viruses, and eukaryota, and putative functional genes that would have gone undetected using 16S amplicon sequencing. NEON open data will be useful for future studies characterizing and quantifying complex ecological processes associated with changing aquatic and terrestrial ecosystems.

Comparative metatranscriptome analysis revealed broad response of microbial communities in two soil types, agriculture versus organic soil

BACKGROUND

Studying expression of genes by direct sequencing and analysis of metatranscriptomes at a particular time and space can disclose structural and functional insights about microbial communities. The present study reports comparative analysis of metatranscriptome from two distinct soil ecosystems referred as M1 (agriculture soil) and O1 (organic soil).

RESULTS

Analysis of sequencing reads revealed Proteobacteria as major dominant phyla in both soil types. The order of the top 3 abundant phyla in M1 sample was Proteobacteria > Ascomycota > Firmicutes, whereas in sample O1, the order was Proteobacteria > Cyanobacteria > Actinobacteria. Analysis of differentially expressed genes demonstrated high expression of transcripts related to copper-binding proteins, proteins involved in electron carrier activity, DNA integration, endonuclease activity, MFS transportation, and other uncharacterized proteins in M1 compared to O1. Of the particular interests, several transcripts related to nitrification, ammonification, stress response, and alternate carbon fixation pathways were highly expressed in M1. In-depth analysis of the sequencing data revealed that transcripts of archaeal origin had high expression in M1 compared to O1 indicating the active role of Archaea in metal- and pesticide-contaminated environment. In addition, transcripts encoding 4-hydroxyphenylpyruvate dioxygenase, glyoxalase/bleomycin resistance protein/dioxygenase, metapyrocatechase, and ring hydroxylating dioxygenases of aromatic hydrocarbon degradation pathways had high expression in M1. Altogether, this study provided important insights about the transcripts and pathways upregulating in the presence of pesticides and herbicides.

CONCLUSION

Altogether, this study claims a high expression of microbial transcripts in two ecosystems with a wide range of functions. It further provided clue about several molecular markers which could be a strong indicator of metal and pesticide contamination in soils. Interestingly, our study revealed that Archaea are playing a significant role in nitrification process as compared to bacteria in metal- and pesticide-contaminated soil. In particular, high expression of transcripts related to aromatic hydrocarbon degradation in M1 soil indicates their important role in biodegradation of pollutants, and therefore, further investigation is needed.

Advances and Challenges in Metatranscriptomic Analysis

Sequencing-based analyses of microbiomes have traditionally focused on addressing the question of community membership and profiling taxonomic abundance through amplicon sequencing of 16 rRNA genes. More recently, shotgun metagenomics, which involves the random sequencing of all genomic content of a microbiome, has dominated this arena due to advancements in sequencing technology throughput and capability to profile genes as well as microbiome membership. While these methods have revealed a great number of insights into a wide variety of microbiomes, both of these approaches only describe the presence of organisms or genes, and not whether they are active members of the microbiome. To obtain deeper insights into how a microbial community responds over time to their changing environmental conditions, microbiome scientists are beginning to employ large-scale metatranscriptomics approaches. Here, we present a comprehensive review on computational metatranscriptomics approaches to study microbial community transcriptomes. We review the major advancements in this burgeoning field, compare strengths and weaknesses to other microbiome analysis methods, list available tools and workflows, and describe use cases and limitations of this method. We envision that this field will continue to grow exponentially, as will the scope of projects (e.g. longitudinal studies of community transcriptional responses to perturbations over time) and the resulting data. This review will provide a list of options for computational analysis of these data and will highlight areas in need of development.

Phycosphere Microbial Succession Patterns and Assembly Mechanisms in a Marine Dinoflagellate Bloom

Given the ecological significance of microorganisms in algal blooming events, it is critical to understand the mechanisms regarding their distribution under different conditions. We tested the hypothesis that microbial community succession is strongly associated with algal bloom stages, and that the assembly mechanisms are cocontrolled by deterministic and stochastic processes. Community structures and underlying ecological processes of microbial populations (attached and free-living bacteria) at three algal bloom stages (pre-, during, and postbloom) over a complete dinoflagellate Scrippsiella trochoidea bloom were investigated. Both attached and free-living taxa had a strong response to the bloom event, and the latter was more sensitive than the former. The contribution of environmental parameters to microbial variability was 40.2%. Interaction analysis showed that complex positive or negative correlation networks exist in phycosphere microbes. These relationships were the potential drivers of mutualist and competitive interactions that impacted bacterial succession. Null model analysis showed that the attached bacterial community primarily exhibited deterministic processes at pre- and during-bloom stages, while dispersal-related processes contributed to a greater extent at the postbloom stage. In the free-living bacterial community, homogeneous selection and dispersal limitation dominated in the initial phase, which gave way to more deterministic processes at the two later stages. Relative contribution analyses further demonstrated that the community turnover of attached bacteria was mainly driven by environmental selection, while stochastic factors had partial effects on the assembly of free-living bacteria. Taken together, these data demonstrated that a robust link exists between bacterioplankton community structure and bloom progression, and phycosphere microbial succession trajectories are cogoverned by both deterministic and random processes.IMPORTANCE Disentangling the mechanisms shaping bacterioplankton communities during a marine ecological event is a core concern for ecologists. Harmful algal bloom (HAB) is a typical ecological disaster, and its formation is significantly influenced by alga-bacterium interactions. Microbial community shifts during the HAB process are relatively well known. However, the assembly processes of microbial communities in an HAB are not fully understood, especially the relative influences of deterministic and stochastic processes. We therefore analyzed the relative contributions of deterministic and stochastic processes during an HAB event. Both free-living and attached bacterial groups had a dramatic response to the HAB, and the relative importance of determinism versus stochasticity varied between the two bacterial groups at various bloom stages. Environmental factors and biotic interactions were the main drivers impacting the microbial shift process. Our results strengthen the understanding of the ecological mechanisms controlling microbial community patterns during the HAB process.

Epigenetic Regulation at the Interplay Between Gut Microbiota and Host Metabolism

Gut microbiota communities have coevolved for millions of years in a symbiotic relationship with their mammalian hosts. Elucidating and understanding the molecular mechanisms by which microbiota interacts with its host and how this contributes to the homeostasis of the host is crucial. One of these molecular relationships is the so-called chemical crosstalk between microbiota and host metabolisms, including the poorly explored epigenetic regulation of host tissues by the metabolic activity of gut microbiota in response to changes in diet. DNA methylation and histone modifications are epigenetic marks partly regulated by enzymes such as methylases and acetylases, whose activity depend on host and microbiota metabolites that act as substrates and cofactors for these reactions. However, providing a complete mechanistic description of the regulatory interactions between both metabolisms and the impact on the expression of host genes through an epigenetic modulation, remains elusive. This article presents our perspective on how metabolomic, metagenomic, transcriptomic, and epigenomic data can be used to investigate the "microbiota-nutrient metabolism-epigenetics axis." We also discuss the implications and opportunities this knowledge may have for basic and applied science, such as the impact on the way we structure future research, understand, and prevent diseases like type 2 diabetes or obesity.

Metagenomics analysis of cocoa bean fermentation microbiome identifying species diversity and putative functional capabilities

Fermentation of Theobroma cacao L. beans is the most critical stage in the production of cocoa products such as chocolates and its derivatives. There is a limited understanding of the complex response of microbial diversity during cocoa bean fermentation. The aim of the present study was to investigate microbial communities in the cocoa bean fermentation heap using a culture-independent approach to elucidate microbial diversity, structure, functional annotation and mapping unto metabolic pathways. Genomic DNA was extracted and purified from a sample of cocoa beans fermentation heap and was followed by library preparations. Sequence data was generated on Illumina Hiseq 2000 paired-end technology (Macrogen Inc). Taxonomic analysis based on genes predicted from the metagenome identified a high percentage of Bacteria (90.0%), Yeast (9%), and bacteriophages (1%) from the cocoa microbiome. Lactobacillus (20%), Gluconacetobacter (9%), Acetobacter (7%) and Gluconobacter (6%) dominated this study. The mean species diversity, measured by Shannon alpha-diversity index, was estimated at 142.81. Assignment of metagenomic sequences to SEED database categories at 97% sequence similarity identified a genetic profile characteristic of heterotrophic lactic acid fermentation of carbohydrates and aromatic amino acids. Metabolism of aromatic compounds, amino acids and their derivatives and carbohydrates occupied 0.6%, 8% and 13% respectively. Overall, these results provide insights into the cocoa microbiome, identifying fermentation processes carried out broadly by complex microbial communities and metabolic pathways encoding aromatic compounds such as phenylacetaldehyde, butanediol, acetoin, and theobromine that are required for flavour and aroma production. The results obtained will help develop targeted inoculations to produce desired chocolate flavour or targeted metabolic pathways for the selection of microbes for good aroma and flavour compounds formation.

A robust, cost-effective method for DNA, RNA and protein co-extraction from soil, other complex microbiomes and pure cultures

The soil microbiome is inherently complex with high biological diversity, and spatial heterogeneity typically occurring on the submillimetre scale. To study the microbial ecology of soils, and other microbiomes, biomolecules, that is, nucleic acids and proteins, must be efficiently and reliably co-recovered from the same biological samples. Commercial kits are currently available for the co-extraction of DNA, RNA and proteins but none has been developed for soil samples. We present a new protocol drawing on existing phenol-chloroform-based methods for nucleic acids co-extraction but incorporating targeted precipitation of proteins from the phenol phase. The protocol is cost-effective and robust, and easily implemented using reagents commonly available in laboratories. The method is estimated to be eight times cheaper than using disparate commercial kits for the isolation of DNA and/or RNA, and proteins, from soil. The method is effective, providing good quality biomolecules from a diverse range of soil types, with clay contents varying from 9.5% to 35.1%, which we successfully used for downstream, high-throughput gene sequencing and metaproteomics. Additionally, we demonstrate that the protocol can also be easily implemented for biomolecule co-extraction from other complex microbiome samples, including cattle slurry and microbial communities recovered from anaerobic bioreactors, as well as from Gram-positive and Gram-negative pure cultures.

Metatranscriptomic exploration of microbial functioning in clouds

Clouds constitute the uppermost layer of the biosphere. They host diverse communities whose functioning remains obscure, although biological activity potentially participates to atmospheric chemical and physical processes. In order to gain information on the metabolic functioning of microbial communities in clouds, we conducted coordinated metagenomics/metatranscriptomics profiling of cloud water microbial communities. Samples were collected from a high altitude atmospheric station in France and examined for biological content after untargeted amplification of nucleic acids. Living microorganisms, essentially bacteria, maintained transcriptional and translational activities and expressed many known complementary physiological responses intended to fight oxidants, osmotic variations and cold. These included activities of oxidant detoxification and regulation, synthesis of osmoprotectants/cryoprotectants, modifications of membranes, iron uptake. Consistently these energy-demanding processes were fueled by central metabolic routes involved in oxidative stress response and redox homeostasis management, such as pentose phosphate and glyoxylate pathways. Elevated binding and transmembrane ion transports demonstrated important interactions between cells and their cloud droplet chemical environments. In addition, polysaccharides, potentially beneficial for survival like exopolysaccharides, biosurfactants and adhesins, were synthesized. Our results support a biological influence on cloud physical and chemical processes, acting notably on the oxidant capacity, iron speciation and availability, amino-acids distribution and carbon and nitrogen fates.

Detection of ribonucleoside modifications by liquid chromatography coupled with mass spectrometry

A small set of ribonucleoside modifications have been found in different regions of mRNA including the open reading frame. Accurate detection of these specific modifications is critical to understanding their modulatory roles in facilitating mRNA maturation, translation and degradation. While transcriptome-wide next-generation sequencing (NGS) techniques could provide exhaustive information about the sites of one specific or class of modifications at a time, recent investigations strongly indicate cautionary interpretation due to the appearance of false positives. Therefore, it is suggested that NGS-based modification data can only be treated as predicted sites and their existence need to be validated by orthogonal methods. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) is an analytical technique that can yield accurate and reproducible information about the qualitative and quantitative characteristics of ribonucleoside modifications. Here, we review the recent advancements in LC-MS/MS technology that could help in securing accurate, gold-standard quality information about the resident post-transcriptional modifications of mRNA.

Metatranscriptome analysis of the intestinal microorganisms in Pardosa pseudoannulata in response to cadmium stress

Cadmium (Cd) generates a variety of physiological and ecological toxicity to spiders. However, little is known about the effects of Cd on symbiotic bacteria of spiders. Metatranscriptomics is increasing our knowledge of microorganisms in environment. To better understand the impact of Cd on the symbiotic bacteria of spiders, we generated and compared the metatranscriptomes of the intestinal microorganisms of Pardosa pseudoannulata with and without Cd stress. The community structure of intestinal microorganisms in P. pseudoannulata was composed of 4 kingdoms, namely bacteria, viruses, eukaryotes and archaea, including 46 phyla, 97 classes, 184 orders, 339 families, 470 genera, and 598 species. The abundance of eukaryotes, bacteria and viruses was decreased by 0.14%, 1.22% and 2.52% respectively while the archaea was increased by 99.16% when under Cd stress. We identified 1519 differentially expressed genes (DEGs), including 770 up-regulated and 749 down-regulated genes. The results of KEGG annotation revealed that the expression of genes that are involved in the carbon metabolism, protein and amino acid metabolism and synthesis, glucose metabolism, oxidative phosphorylation, and glutathione metabolism were influenced by Cd. Collectively, these findings showed that Cd significantly impacted the community structure and expression of related functional genes of intestinal microorganisms in P. pseudoannulata.

Selection of marker genes for genetic barcoding of microorganisms and binning of metagenomic reads by Barcoder software tools

Background

Metagenomic approaches have revealed the complexity of environmental microbiomes with the advancement in whole genome sequencing displaying a significant level of genetic heterogeneity on the species level. It has become apparent that patterns of superior bioactivity of bacteria applicable in biotechnology as well as the enhanced virulence of pathogens often requires distinguishing between closely related species or sub-species. Current methods for binning of metagenomic reads usually do not allow for identification below the genus level and generally stops at the family level.

Results

In this work, an attempt was made to improve metagenomic binning resolution by creating genome specific barcodes based on the core and accessory genomes. This protocol was implemented in novel software tools available for use and download from http://bargene.bi.up.ac.za/. The most abundant barcode genes from the core genomes were found to encode for ribosomal proteins, certain central metabolic genes and ABC transporters. Performance of metabarcode sequences created by this package was evaluated using artificially generated and publically available metagenomic datasets. Furthermore, a program (Barcoding 2.0) was developed to align reads against barcode sequences and thereafter calculate various parameters to score the alignments and the individual barcodes. Taxonomic units were identified in metagenomic samples by comparison of the calculated barcode scores to set cut-off values. In this study, it was found that varying sample sizes, i.e. number of reads in a metagenome and metabarcode lengths, had no significant effect on the sensitivity and specificity of the algorithm. Receiver operating characteristics (ROC) curves were calculated for different taxonomic groups based on the results of identification of the corresponding genomes in artificial metagenomic datasets. The reliability of distinguishing between species of the same genus or family by the program was nearly perfect.

Conclusions

The results showed that the novel online tool BarcodeGenerator (http://bargene.bi.up.ac.za/) is an efficient approach for generating barcode sequences from a set of complete genomes provided by users. Another program, Barcoder 2.0 is available from the same resource to enable an efficient and practical use of metabarcodes for visualization of the distribution of organisms of interest in environmental and clinical samples.

Short-Term Transcriptional Response of Microbial Communities to Nitrogen Fertilization in a Pine Forest Soil

Numerous studies have examined the long-term effect of experimental nitrogen (N) deposition in terrestrial ecosystems; however, N-specific mechanistic markers are difficult to disentangle from responses to other environmental changes. The strongest picture of N-responsive mechanistic markers is likely to arise from measurements over a short (hours to days) time scale immediately after inorganic N deposition. Therefore, we assessed the short-term (3-day) transcriptional response of microbial communities in two soil strata from a pine forest to a high dose of N fertilization (ca. 1 mg/g of soil material) in laboratory microcosms. We hypothesized that N fertilization would repress the expression of fungal and bacterial genes linked to N mining from plant litter. However, despite N suppression of microbial respiration, the most pronounced differences in functional gene expression were between strata rather than in response to the N addition. Overall, ∼4% of metabolic genes changed in expression with N addition, while three times as many (∼12%) were significantly different across the different soil strata in the microcosms. In particular, we found little evidence of N changing expression levels of metabolic genes associated with complex carbohydrate degradation (CAZymes) or inorganic N utilization. This suggests that direct N repression of microbial functional gene expression is not the principle mechanism for reduced soil respiration immediately after N deposition. Instead, changes in expression with N addition occurred primarily in general cell maintenance areas, for example, in ribosome-related transcripts. Transcriptional changes in functional gene abundance in response to N addition observed in longer-term field studies likely result from changes in microbial composition.IMPORTANCE Ecosystems are receiving increased nitrogen (N) from anthropogenic sources, including fertilizers and emissions from factories and automobiles. High levels of N change ecosystem functioning. For example, high inorganic N decreases the microbial decomposition of plant litter, potentially reducing nutrient recycling for plant growth. Understanding how N regulates microbial decomposition can improve the prediction of ecosystem functioning over extended time scales. We found little support for the conventional view that high N supply represses the expression of genes involved in decomposition or alters the expression of bacterial genes for inorganic N cycling. Instead, our study of pine forest soil 3 days after N addition showed changes in microbial gene expression related to cell maintenance and stress response. This highlights the challenge of establishing predictive links between microbial gene expression levels and measures of ecosystem function.

Omics and multi-omics approaches to study the biosynthesis of secondary metabolites in microorganisms

Natural products produced by microorganisms represent the main source of bioactive molecules. The development of high-throughput (omics) techniques have importantly contributed to the renaissance of new antibiotic discovery increasing our understanding of complex mechanisms controlling the expression of biosynthetic gene clusters (BGCs) encoding secondary metabolites. In this context this review highlights recent progress in the use and integration of 'omics' approaches with focuses on genomics, transcriptomics, proteomics metabolomics meta-omics and combined omics as powerful strategy to discover new antibiotics.

Occurrence and seasonal dynamics of RNA viral genotypes in three contrasting temperate lakes

Decades of research have demonstrated the crucial importance of viruses in freshwater ecosystems. However, few studies have focused on the seasonal dynamics and potential hosts of RNA viruses. We surveyed microbial-sized (i.e. 5-0.2 μm) mixed community plankton transcriptomes for RNA viral genomes and investigated their distribution between microbial and macrobial plankton over a seasonal cycle across three temperate lakes by quantitative reverse transcriptase PCR (qRT-PCR). A total of 30 contigs bearing similarity to RNA viral genomes were recovered from a global assembly of 30 plankton RNA libraries. Of these, only 13 were found in >2 libraries and recruited >100 reads (of 9.13 x 107 total reads), representing several picornaviruses, two tobamoviruses and a reovirus. We quantified the abundance of four picornaviruses and the reovirus monthly from August 2014 to May 2015. Patterns of viral abundance in the >5 μm size fraction and representation in microbial-sized community RNA libraries over time suggest that one picornavirus genotype (TS24835) and the reovirus (TS148892) may infect small (<5 μm) eukaryotic microorganisms, while two other picornaviruses (TS24641 and TS4340) may infect larger (>5 μm) eukaryotic microorganisms or metazoa. Our data also suggest that picornavirus TS152062 may originate from an allochthonous host. All five viral genotypes were present in at least one size fraction across all 3 lakes during the year, suggesting that RNA viruses may easily disperse between adjacent aquatic habitats. Our data therefore demonstrate that RNA viruses are widespread in temperate lacustrine ecosystems, and may provide evidence of viral infection in larger eukaryotes (including metazoa) inhabiting the lakes.

Integrating omics to unravel the stress-response mechanisms in probiotic bacteria: Approaches, challenges, and prospects

Identifying the stress-response mechanism of probiotic bacteria has always captivated the interest of food producers. It is crucial to identify probiotic bacteria that have increased stress tolerance to survive during production, processing, and storage of food products. However, in order to achieve high resistance to environmental factors, there is a need to better understand stress-induced responses and adaptive mechanisms. With advances in bacterial genomics, there has been an upsurge in the application of other omic platforms such as transcriptomics, proteomics, metabolomics, and some more recent ones such as interactomics, fluxomics, and phenomics. These omic technologies have revolutionized the functional genomics and their application. There have been several studies implementing various omic technologies to investigate the stress responses of probiotic bacteria. Integrated omics has the potential to provide in-depth information about the mechanisms of stress-induced responses in bacteria. However, there remain challenges in integrating information from different omic platforms. This review discusses current omic techniques and challenges faced in integrating various omic platforms with focus on their use in stress-response studies.

RNA Sequencing and Analysis in Microorganisms for Metabolic Network Reconstruction

There is a strict interplay between metabolic networks and transcriptional regulation in bacteria; indeed, the transcriptome regulation, affecting the expression of large gene sets, can be used to predict the likely "on" or "off" state of metabolic genes as a function of environmental factors. Up to date, many bacterial transcriptomes have been studied by RNAseq, hundreds of experiments have been performed, and Giga bases of sequences have been produced. All this transcriptional information could potentially be integrated into metabolic networks in order to obtain a more comprehensive view of their regulation and to increase their prediction power.To get high-quality transcriptomic data, to be integrated into metabolic networks, it is paramount to clearly know how to produce highly informative RNA sequencing libraries and how to manage RNA sequencing data.In this chapter, we will get across the main steps of an RNAseq experiment: from removal of ribosomal RNAs, to strand-specific library preparation, till data analysis and integration. We will try to share our experience and know-how, to give you a precise protocol to follow, and some useful recommendations or tips and tricks to adopt in order to go straightforward toward a successful RNAseq experiment.

Metatranscriptomic and metagenomic description of the bacterial nitrogen metabolism in waste water wet oxidation effluents

Anaerobic digestion is a common method for reducing the amount of sludge solids in used waters and enabling biogas production. The wet oxidation process (WOX) improves anaerobic digestion by converting carbon into methane through oxidation of organic compounds. WOX produces effluents rich in ammonia, which must be removed to maintain the activity of methanogens. Ammonia removal from WOX could be biologically operated by aerobic granules. To this end, granulation experiments were conducted in 2 bioreactors containing an activated sludge (AS). For the first time, the dynamics of the microbial community structure and the expression levels of 7 enzymes of the nitrogen metabolism in such active microbial communities were followed in regard to time by metagenomics and metatranscriptomics. It was shown that bacterial communities adapt to the wet oxidation effluent by increasing the expression level of the nitrogen metabolism, suggesting that these biological activities could be a less costly alternative for the elimination of ammonia, resulting in a reduction of the use of chemicals and energy consumption in sewage plants. This study reached a strong sequencing depth (from 4.4 to 7.6 Gb) and enlightened a yet unknown diversity of the microorganisms involved in the nitrogen pathway. Moreover, this approach revealed the abundance and expression levels of specialised enzymes involved in nitrification, denitrification, ammonification, dissimilatory nitrate reduction to ammonium (DNRA) and nitrogen fixation processes in AS.

How to normalize metatranscriptomic count data for differential expression analysis

Background

Differential expression analysis on the basis of RNA-Seq count data has become a standard tool in transcriptomics. Several studies have shown that prior normalization of the data is crucial for a reliable detection of transcriptional differences. Until now it has not been clear whether and how the transcriptomic approach can be used for differential expression analysis in metatranscriptomics.

Methods

We propose a model for differential expression in metatranscriptomics that explicitly accounts for variations in the taxonomic composition of transcripts across different samples. As a main consequence the correct normalization of metatranscriptomic count data under this model requires the taxonomic separation of the data into organism-specific bins. Then the taxon-specific scaling of organism profiles yields a valid normalization and allows us to recombine the scaled profiles into a metatranscriptomic count matrix. This matrix can then be analyzed with statistical tools for transcriptomic count data. For taxon-specific scaling and recombination of scaled counts we provide a simple R script.

Results

When applying transcriptomic tools for differential expression analysis directly to metatranscriptomic data with an organism-independent (global) scaling of counts the resulting differences may be difficult to interpret. The differences may correspond to changing functional profiles of the contributing organisms but may also result from a variation of taxonomic abundances. Taxon-specific scaling eliminates this variation and therefore the resulting differences actually reflect a different behavior of organisms under changing conditions. In simulation studies we show that the divergence between results from global and taxon-specific scaling can be drastic. In particular, the variation of organism abundances can imply a considerable increase of significant differences with global scaling. Also, on real metatranscriptomic data, the predictions from taxon-specific and global scaling can differ widely. Our studies indicate that in real data applications performed with global scaling it might be impossible to distinguish between differential expression in terms of transcriptomic changes and differential composition in terms of changing taxonomic proportions.

Conclusions

As in transcriptomics, a proper normalization of count data is also essential for differential expression analysis in metatranscriptomics. Our model implies a taxon-specific scaling of counts for normalization of the data. The application of taxon-specific scaling consequently removes taxonomic composition variations from functional profiles and therefore provides a clear interpretation of the observed functional differences.

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.

Thriving or surviving? Evaluating active microbial guilds in Baltic Sea sediment

Microbial life in the deep subsurface biosphere is taxonomically and metabolically diverse, but it is vigorously debated whether the resident organisms are thriving (metabolizing, maintaining cellular integrity and expressing division genes) or just surviving. As part of Integrated Ocean Drilling Program Expedition 347: Baltic Sea Paleoenvironment, we extracted and sequenced RNA from organic carbon-rich, nutrient-replete and permanently anoxic sediment. In stark contrast to the oligotrophic subsurface biosphere, Baltic Sea Basin samples provided a unique opportunity to understand the balance between metabolism and other cellular processes. Targeted sequencing of 16S rRNA transcripts showed Atribacteria (an uncultured phylum) and Chloroflexi to be among the dominant and the active members of the community. Metatranscriptomic analysis identified methane cycling, sulfur cycling and halogenated compound utilization as active in situ respiratory metabolisms. Genes for cellular maintenance, cellular division, motility and antimicrobial production were also transcribed. This indicates that microbial life in deep subsurface Baltic Sea Basin sediments was not only alive, but thriving.

Strand-specific transcriptomes of Enterohemorrhagic Escherichia coli in response to interactions with ground beef microbiota: interactions between microorganisms in raw meat

BACKGROUND

Enterohemorrhagic Escherichia coli (EHEC) are zoonotic agents associated with outbreaks worldwide. Growth of EHEC strains in ground beef could be inhibited by background microbiota that is present initially at levels greater than that of the pathogen E. coli. However, how the microbiota outcompetes the pathogenic bacteria is unknown. Our objective was to identify metabolic pathways of EHEC that were altered by natural microbiota in order to improve our understanding of the mechanisms controlling the growth and survival of EHECs in ground beef.

RESULTS

Based on 16S metagenomics analysis, we identified the microbial community structure in our beef samples which was an essential preliminary for subtractively analyzing the gene expression of the EHEC strains. Then, we applied strand-specific RNA-seq to investigate the effects of this microbiota on the global gene expression of EHEC O2621765 and O157EDL933 strains by comparison with their behavior in beef meat without microbiota. In strain O2621765, the expression of genes connected with nitrate metabolism and nitrite detoxification, DNA repair, iron and nickel acquisition and carbohydrate metabolism, and numerous genes involved in amino acid metabolism were down-regulated. Further, the observed repression of ftsL and murF, involved respectively in building the cytokinetic ring apparatus and in synthesizing the cytoplasmic precursor of cell wall peptidoglycan, might help to explain the microbiota's inhibitory effect on EHECs. For strain O157EDL933, the induced expression of the genes implicated in detoxification and the general stress response and the repressed expression of the peR gene, a gene negatively associated with the virulence phenotype, might be linked to the survival and virulence of O157:H7 in ground beef with microbiota.

CONCLUSION

In the present study, we show how RNA-Seq coupled with a 16S metagenomics analysis can be used to identify the effects of a complex microbial community on relevant functions of an individual microbe within it. These findings add to our understanding of the behavior of EHECs in ground beef. By measuring transcriptional responses of EHEC, we could identify putative targets which may be useful to develop new strategies to limit their shedding in ground meat thus reducing the risk of human illnesses.

Viral to metazoan marine plankton nucleotide sequences from the Tara Oceans expedition

A unique collection of oceanic samples was gathered by the Tara Oceans expeditions (2009–2013), targeting plankton organisms ranging from viruses to metazoans, and providing rich environmental context measurements. Thanks to recent advances in the field of genomics, extensive sequencing has been performed for a deep genomic analysis of this huge collection of samples. A strategy based on different approaches, such as metabarcoding, metagenomics, single-cell genomics and metatranscriptomics, has been chosen for analysis of size-fractionated plankton communities. Here, we provide detailed procedures applied for genomic data generation, from nucleic acids extraction to sequence production, and we describe registries of genomics datasets available at the European Nucleotide Archive (ENA, www.ebi.ac.uk/ena). The association of these metadata to the experimental procedures applied for their generation will help the scientific community to access these data and facilitate their analysis. This paper complements other efforts to provide a full description of experiments and open science resources generated from the Tara Oceans project, further extending their value for the study of the world’s planktonic ecosystems.

Linking Genes to Microbial Biogeochemical Cycling: Lessons from Arsenic

The biotransformation of arsenic is highly relevant to the arsenic biogeochemical cycle. Identification of the molecular details of microbial pathways of arsenic biotransformation coupled with analyses of microbial communities by meta-omics can provide insights into detailed aspects of the complexities of this biocycle. Arsenic transformations couple to other biogeochemical cycles, and to the fate of both nutrients and other toxic environmental contaminants. Microbial redox metabolism of iron, carbon, sulfur, and nitrogen affects the redox and bioavailability of arsenic species. In this critical review we illustrate the biogeochemical processes and genes involved in arsenic biotransformations. We discuss how current and future metagenomic-, metatranscriptomic-, metaproteomic-, and metabolomic-based methods will help to decipher individual microbial arsenic transformation processes, and their connections to other biogeochemical cycle. These insights will allow future use of microbial metabolic capabilities for new biotechnological solutions to environmental problems. To understand the complex nature of inorganic and organic arsenic species and the fate of environmental arsenic will require integrating systematic approaches with biogeochemical modeling. Finally, from the lessons learned from these studies of arsenic biogeochemistry, we will be able to predict how the environment changes arsenic, and, in response, how arsenic biotransformations change the environment.

Towards a sustainable biobased industry - Highlighting the impact of extremophiles

The transition of the oil-based economy towards a sustainable economy completely relying on biomass as renewable feedstock requires the concerted action of academia, industry, politics and civil society. An interdisciplinary approach of various fields such as microbiology, molecular biology, chemistry, genetics, chemical engineering and agriculture in addition to cross-sectional technologies such as economy, logistics and digitalization is necessary to meet the future global challenges. The genomic era has contributed significantly to the exploitation of naturés biodiversity also from extreme habitats. By applying modern technologies it is now feasible to deliver robust enzymes (extremozymes) and robust microbial systems that are active at temperatures up to 120°C, at pH 0 and 12 and at 1000bar. In the post-genomic era, different sophisticated "omics" analyses will allow the identification of countless novel enzymes regardless of the lack of cultivability of most microorganisms. Furthermore, elaborate protein-engineering methods are clearing the way towards tailor-made robust biocatalysts. Applying environmentally friendly and efficient biological processes, terrestrial and marine biomass can be converted to high value products e.g. chemicals, building blocks, biomaterials, pharmaceuticals, food, feed and biofuels. Thus, further application of extremophiles has the potential to improve sustainability of existing biotechnological processes towards a greener biobased industry.

Effects of Gelling Agent and Extracellular Signaling Molecules on the Culturability of Marine Bacteria

Only 1% of marine bacteria are currently culturable using standard laboratory procedures, and this is a major obstacle for our understanding of the biology of marine microorganisms and for the discovery of novel microbial natural products. Therefore, the purpose of this study was to investigate if improved cultivation conditions, including the use of an alternative gelling agent and supplementation with signaling molecules, improve the culturability of bacteria from seawater. Replacing agar with gellan gum improved viable counts 3- to 40-fold, depending on medium composition and incubation conditions, with a maximum of 6.6% culturability relative to direct cell counts. Through V4 amplicon sequencing we found that culturable diversity was also affected by a change in gelling agent, facilitating the growth of orders not culturable on agar-based substrates. Community analyses showed that communities grown on gellan gum substrates were significantly different from communities grown on agar and that they covered a larger fraction of the seawater community. Other factors, such as incubation temperature and time, had less obvious effects on viable counts and culturable diversity. Supplementation with acylated homoserine lactones (AHLs) did not have a positive effect on total viable counts or a strong effect on culturable diversity. However, low concentrations of AHLs increased the relative abundance of sphingobacteria. Hence, with alternative growth substrates, it is possible to significantly increase the number and diversity of cultured marine bacteria.IMPORTANCE Serious challenges to human health, such as the occurrence and spread of antibiotic resistance and an aging human population in need of bioactive pharmaceuticals, have revitalized the search for natural microbial products. The marine environment, representing the largest ecosystem in the biosphere, harbors an immense and virtually untapped microbial diversity producing unique bioactive compounds. However, we are currently able to cultivate only a minute fraction of this diversity. The lack of cultivated microbes is hampering not only bioprospecting efforts but also our general understanding of marine microbes. In this study, we present a means to increase the number and diversity of cultured bacteria from seawater, showing that relatively simple changes to medium components may facilitate the isolation and growth of hitherto unknown bacterial orders.

Functional gene profiling through metaRNAseq approach reveals diet-dependent variation in rumen microbiota of buffalo (Bubalus bubalis)

Recent advances in next generation sequencing technology have enabled analysis of complex microbial community from genome to transcriptome level. In the present study, metatranscriptomic approach was applied to elucidate functionally active bacteria and their biological processes in rumen of buffalo (Bubalus bubalis) adapted to different dietary treatments. Buffaloes were adapted to a diet containing 50:50, 75:25 and 100:0 forage to concentrate ratio, each for 6 weeks, before ruminal content sample collection. Metatranscriptomes from rumen fiber adherent and fiber-free active bacteria were sequenced using Ion Torrent PGM platform followed by annotation using MG-RAST server and CAZYmes (Carbohydrate active enzymes) analysis toolkit. In all the samples Bacteroidetes was the most abundant phylum followed by Firmicutes. Functional analysis using KEGG Orthology database revealed Metabolism as the most abundant category at level 1 within which Carbohydrate metabolism was dominating. Diet treatments also exerted significant differences in proportion of enzymes involved in metabolic pathways for VFA production. Carbohydrate Active Enzyme(CAZy) analysis revealed the abundance of genes encoding glycoside hydrolases with the highest representation of GH13 CAZy family in all the samples. The findings provide an overview of the activities occurring in the rumen as well as active bacterial population and the changes occurring through different dietary treatments.

Exploring the Impacts of Anthropogenic Disturbance on Seawater and Sediment Microbial Communities in Korean Coastal Waters Using Metagenomics Analysis

The coastal ecosystems are considered as one of the most dynamic and vulnerable environments under various anthropogenic developments and the effects of climate change. Variations in the composition and diversity of microbial communities may be a good indicator for determining whether the marine ecosystems are affected by complex forcing stressors. DNA sequence-based metagenomics has recently emerged as a promising tool for analyzing the structure and diversity of microbial communities based on environmental DNA (eDNA). However, few studies have so far been performed using this approach to assess the impacts of human activities on the microbial communities in marine systems. In this study, using metagenomic DNA sequencing (16S ribosomal RNA gene), we analyzed and compared seawater and sediment communities between sand mining and control (natural) sites in southern coastal waters of Korea to assess whether anthropogenic activities have significantly affected the microbial communities. The sand mining sites harbored considerably lower levels of microbial diversities in the surface seawater community during spring compared with control sites. Moreover, the sand mining areas had distinct microbial taxonomic group compositions, particularly during spring season. The microbial groups detected solely in the sediment load/dredging areas (e.g., Marinobacter, Alcanivorax, Novosphingobium) are known to be involved in degradation of toxic chemicals such as hydrocarbon, oil, and aromatic compounds, and they also contain potential pathogens. This study highlights the versatility of metagenomics in monitoring and diagnosing the impacts of human disturbance on the environmental health of marine ecosystems from eDNA.