Flexible metagenome analysis using the MGX framework

Exploring Cereal Metagenomics: Unravelling Microbial Communities for Improved Food Security

Food security is an urgent global challenge, with cereals playing a crucial role in meeting the nutritional requirements of populations worldwide. In recent years, the field of metagenomics has emerged as a powerful tool for studying the microbial communities associated with cereal crops and their impact on plant health and growth. This chapter aims to provide a comprehensive overview of cereal metagenomics and its role in enhancing food security through the exploration of beneficial and pathogenic microbial interactions. Furthermore, we will examine how the integration of metagenomics with other tools can effectively address the adverse effects on food security. For this purpose, we discuss the integration of metagenomic data and machine learning in providing novel insights into the dynamic interactions shaping plant-microbe relationships. We also shed light on the potential applications of leveraging microbial diversity and epigenetic modifications in improving crop resilience and yield sustainability. Ultimately, cereal metagenomics has revolutionized the field of food security by harnessing the potential of beneficial interactions between cereals and their microbiota, paving the way for sustainable agricultural practices.

Long-term conservation tillage with reduced nitrogen fertilization intensity can improve winter wheat health via positive plant-microorganism feedback in the rhizosphere

Microbiome-based solutions are regarded key for sustainable agroecosystems. However, it is unclear how agricultural practices affect the rhizosphere microbiome, plant-microorganism interactions and crop performance under field conditions. Therefore, we installed root observation windows in a winter wheat field cultivated either under long-term mouldboard plough (MP) or cultivator tillage (CT). Each tillage practice was also compared at two nitrogen (N) fertilization intensities, intensive (recommended N-supply with pesticides/growth regulators) or extensive (reduced N-supply, no fungicides/growth regulators). Shoot biomass, root exudates and rhizosphere metabolites, physiological stress indicators, and gene expression were analyzed together with the rhizosphere microbiome (bacterial/archaeal 16S rRNA gene, fungal ITS amplicon, and shotgun metagenome sequencing) shortly before flowering. Compared to MP, the rhizosphere of CT winter wheat contained more primary and secondary metabolites, especially benzoxazinoid derivatives. Potential copiotrophic and plant-beneficial taxa (e.g. Bacillus, Devosia, and Trichoderma) as well as functional genes (e.g. siderophore production, trehalose synthase, and ACC deaminase) were enriched in the CT rhizosphere, suggesting that tillage affected belowground plant-microorganism interactions. In addition, physiological stress markers were suppressed in CT winter wheat compared to MP. In summary, tillage practice was a major driver of crop performance, root deposits, and rhizosphere microbiome interactions, while the N-fertilization intensity was also relevant, but less important.

Uncovering Microbiome Adaptations in a Full-Scale Biogas Plant: Insights from MAG-Centric Metagenomics and Metaproteomics

The current focus on renewable energy in global policy highlights the importance of methane production from biomass through anaerobic digestion (AD). To improve biomass digestion while ensuring overall process stability, microbiome-based management strategies become more important. In this study, metagenomes and metaproteomes were used for metagenomically assembled genome (MAG)-centric analyses to investigate a full-scale biogas plant consisting of three differentially operated digesters. Microbial communities were analyzed regarding their taxonomic composition, functional potential, as well as functions expressed on the proteome level. Different abundances of genes and enzymes related to the biogas process could be mostly attributed to different process parameters. Individual MAGs exhibiting different abundances in the digesters were studied in detail, and their roles in the hydrolysis, acidogenesis and acetogenesis steps of anaerobic digestion could be assigned. Methanoculleus thermohydrogenotrophicum was an active hydrogenotrophic methanogen in all three digesters, whereas Methanothermobacter wolfeii was more prevalent at higher process temperatures. Further analysis focused on MAGs, which were abundant in all digesters, indicating their potential to ensure biogas process stability. The most prevalent MAG belonged to the class Limnochordia; this MAG was ubiquitous in all three digesters and exhibited activity in numerous pathways related to different steps of AD.

RecruitPlotEasy: An Advanced Read Recruitment Plot Tool for Assessing Metagenomic Population Abundance and Genetic Diversity

Mapping of short metagenomic (or metatranscriptomic) read data to reference isolate or single-cell genomes or metagenome-assembled genomes (MAGs) to assess microbial population relative abundance and/or structure represents an essential task of many studies across environmental and clinical settings. The filtering for the quality of the read match and assessment of read mapping results are frequently performed without visual aids or with the assistance of visualizations produced through ad-hoc, in-house approaches. Here, we introduce RecruitPlotEasy, a fully automated, user-friendly pipeline for these purposes that integrates statistical approaches to quantify intra-population sequence and gene-content diversity and identify co-occurring relative populations in the sample. Hence, RecruitPlotEasy should also greatly facilitate population genetics studies.

RecruitPlotEasy is implemented in Python and R languages and is freely available open source software under the Artistic License 2.0 from https://github.com/KGerhardt/RecruitPlotEasy.

Isolation and Characterization of Barley (Hordeum vulgare) Extracellular Vesicles to Assess Their Role in RNA Spray-Based Crop Protection

The demonstration that spray-induced gene silencing (SIGS) can confer strong disease resistance, bypassing the laborious and time-consuming transgenic expression of double-stranded (ds)RNA to induce the gene silencing of pathogenic targets, was ground-breaking. However, future field applications will require fundamental mechanistic knowledge of dsRNA uptake, processing, and transfer. There is increasing evidence that extracellular vesicles (EVs) mediate the transfer of transgene-derived small interfering (si)RNAs in host-induced gene silencing (HIGS) applications. In this study, we establish a protocol for barley EV isolation and assess the possibilities for EVs regarding the translocation of sprayed dsRNA from barley (Hordeum vulgare) to its interacting fungal pathogens. We found barley EVs that were 156 nm in size, containing predominantly 21 and 19 nucleotide (nts) siRNAs, starting with a 5′-terminal Adenine. Although a direct comparison of the RNA cargo between HIGS and SIGS EV isolates is improper given their underlying mechanistic differences, we identified sequence-identical siRNAs in both systems. Overall, the number of siRNAs isolated from the EVs of dsRNA-sprayed barley plants with sequence complementarity to the sprayed dsRNA precursor was low. However, whether these few siRNAs are sufficient to induce the SIGS of pathogenic target genes requires further research. Taken together, our results raise the possibility that EVs may not be mandatory for the spray-delivered siRNA uptake and induction of SIGS.

Cross-Sectional Variations in Structure and Function of Coral Reef Microbiome With Local Anthropogenic Impacts on the Kenyan Coast of the Indian Ocean

Coral reefs face an increased number of environmental threats from anthropomorphic climate change and pollution from agriculture, industries and sewage. Because environmental changes lead to their compositional and functional shifts, coral reef microbial communities can serve as indicators of ecosystem impacts through development of rapid and inexpensive molecular monitoring tools. Little is known about coral reef microbial communities of the Western Indian Ocean (WIO). We compared taxonomic and functional diversity of microbial communities inhabiting near-coral seawater and sediments from Kenyan reefs exposed to varying impacts of human activities. Over 19,000 species (bacterial, viral and archaeal combined) and 4,500 clusters of orthologous groups of proteins (COGs) were annotated. The coral reefs showed variations in the relative abundances of ecologically significant taxa, especially copiotrophic bacteria and coliphages, corresponding to the magnitude of the neighboring human impacts in the respective sites. Furthermore, the near-coral seawater and sediment metagenomes had an overrepresentation of COGs for functions related to adaptation to diverse environments. Malindi and Mombasa marine parks, the coral reef sites closest to densely populated settlements were significantly enriched with genes for functions suggestive of mitigation of environment perturbations including the capacity to reduce intracellular levels of environmental contaminants and repair of DNA damage. Our study is the first metagenomic assessment of WIO coral reef microbial diversity which provides a much-needed baseline for the region, and points to a potential area for future research toward establishing indicators of environmental perturbations.

KAUST Metagenomic Analysis Platform (KMAP), enabling access to massive analytics of re-annotated metagenomic data

Exponential rise of metagenomics sequencing is delivering massive functional environmental genomics data. However, this also generates a procedural bottleneck for on-going re-analysis as reference databases grow and methods improve, and analyses need be updated for consistency, which require acceess to increasingly demanding bioinformatic and computational resources. Here, we present the KAUST Metagenomic Analysis Platform (KMAP), a new integrated open web-based tool for the comprehensive exploration of shotgun metagenomic data. We illustrate the capacities KMAP provides through the re-assembly of ~ 27,000 public metagenomic samples captured in ~ 450 studies sampled across ~ 77 diverse habitats. A small subset of these metagenomic assemblies is used in this pilot study grouped into 36 new habitat-specific gene catalogs, all based on full-length (complete) genes. Extensive taxonomic and gene annotations are stored in Gene Information Tables (GITs), a simple tractable data integration format useful for analysis through command line or for database management. KMAP pilot study provides the exploration and comparison of microbial GITs across different habitats with over 275 million genes. KMAP access to data and analyses is available at https://www.cbrc.kaust.edu.sa/aamg/kmap.start .

Genomic and proteomic profiles of biofilms on microplastics are decoupled from artificial surface properties

Microplastics in marine ecosystems are colonized by diverse prokaryotic and eukaryotic communities. How these communities and their functional profiles are shaped by the artificial surfaces remains broadly unknown. In order to close this knowledge gap, we set up an in situ experiment with pellets of the polyolefin polymer polyethylene (PE), the aromatic hydrocarbon polymer polystyrene (PS), and wooden beads along a coastal to estuarine gradient in the Baltic Sea, Germany. We used an integrated metagenomics/metaproteomics approach to evaluate the genomic potential as well as protein expression levels of aquatic plastic biofilms. Our results suggest that material properties had a minor influence on the plastic-associated assemblages, as genomic and proteomic profiles of communities associated with the structurally different polymers PE and PS were highly similar, hence polymer-unspecific. Instead, it seemed that these communities were shaped by biogeographic factors. Wood, on the other hand, induced the formation of substrate-specific biofilms and served as nutrient source itself. Our study indicates that, while PE and PS microplastics may be relevant in the photic zone as opportunistic colonization grounds for phototrophic microorganisms, they appear not to be subject to biodegradation or serve as vectors for pathogenic microorganisms in marine habitats.

Legume-microbiome interactions unlock mineral nutrients in regrowing tropical forests

Legume trees form an abundant and functionally important component of tropical forests worldwide with N2-fixing symbioses linked to enhanced growth and recruitment in early secondary succession. However, it remains unclear how N2-fixers meet the high demands for inorganic nutrients imposed by rapid biomass accumulation on nutrient-poor tropical soils. Here, we show that N2-fixing trees in secondary Neotropical forests triggered twofold higher in situ weathering of fresh primary silicates compared to non-N2-fixing trees and induced locally enhanced nutrient cycling by the soil microbiome community. Shotgun metagenomic data from weathered minerals support the role of enhanced nitrogen and carbon cycling in increasing acidity and weathering. Metagenomic and marker gene analyses further revealed increased microbial potential beneath N2-fixers for anaerobic iron reduction, a process regulating the pool of phosphorus bound to iron-bearing soil minerals. We find that the Fe(III)-reducing gene pool in soil is dominated by acidophilic Acidobacteria, including a highly abundant genus of previously undescribed bacteria, Candidatus Acidoferrum, genus novus. The resulting dependence of the Fe-cycling gene pool to pH determines the high iron-reducing potential encoded in the metagenome of the more acidic soils of N2-fixers and their nonfixing neighbors. We infer that by promoting the activities of a specialized local microbiome through changes in soil pH and C:N ratios, N2-fixing trees can influence the wider biogeochemical functioning of tropical forest ecosystems in a manner that enhances their ability to assimilate and store atmospheric carbon.

Compared to conventional, ecological intensive management promotes beneficial proteolytic soil microbial communities for agro-ecosystem functioning under climate change-induced rain regimes

Projected climate change and rainfall variability will affect soil microbial communities, biogeochemical cycling and agriculture. Nitrogen (N) is the most limiting nutrient in agroecosystems and its cycling and availability is highly dependent on microbial driven processes. In agroecosystems, hydrolysis of organic nitrogen (N) is an important step in controlling soil N availability. We analyzed the effect of management (ecological intensive vs. conventional intensive) on N-cycling processes and involved microbial communities under climate change-induced rain regimes. Terrestrial model ecosystems originating from agroecosystems across Europe were subjected to four different rain regimes for 263 days. Using structural equation modelling we identified direct impacts of rain regimes on N-cycling processes, whereas N-related microbial communities were more resistant. In addition to rain regimes, management indirectly affected N-cycling processes via modifications of N-related microbial community composition. Ecological intensive management promoted a beneficial N-related microbial community composition involved in N-cycling processes under climate change-induced rain regimes. Exploratory analyses identified phosphorus-associated litter properties as possible drivers for the observed management effects on N-related microbial community composition. This work provides novel insights into mechanisms controlling agro-ecosystem functioning under climate change.

Illustration of the microbial community selected by optimized process and nutritional parameters resulting in enhanced biomethanation of rice straw without thermo-chemical pretreatment

Effects of different process and nutritional parameters on microbial community structure and function were investigated to enhance the biomethanation of rice straw without any thermochemical pre-treatment. The study was performed in a mesophilic anaerobic digester with cattle dung slurry as inoculum. The highest methane yield of 274 ml g^-1 volatile solids was obtained from particulate rice straw (1 mm size, 7.5% solids loading rate) at 37 °C, pH-7, when supplemented with urea (carbon: nitrogen ratio, 25:1) and zinc as trace element (100 µM) at 21 days hydraulic retention time. The optimization of conditions selected Clostridium, Bacteroides, and Ruminococcus as dominant hydrolytic bacteria and Methanosarcina as the methanogen. Analysis of metagenome and metatranscriptome revealed wide array of bacterial lignocellulolytic enzymes that efficiently hydrolyzed the rice straw. The methane yield was >80% of the theoretical yield, making this green process a sustainable choice for efficient extraction of energy from rice straw.

Effect of Long-Term Farming Practices on Agricultural Soil Microbiome Members Represented by Metagenomically Assembled Genomes (MAGs) and Their Predicted Plant-Beneficial Genes

To follow the hypothesis that agricultural management practices affect structure and function of the soil microbiome regarding soil health and plant-beneficial traits, high-throughput (HT) metagenome analyses were performed on Chernozem soil samples from a long-term field experiment designated LTE-1 carried out at Bernburg-Strenzfeld (Saxony-Anhalt, Germany). Metagenomic DNA was extracted from soil samples representing the following treatments: (i) plough tillage with standard nitrogen fertilization and use of fungicides and growth regulators, (ii) plough tillage with reduced nitrogen fertilization (50%), (iii) cultivator tillage with standard nitrogen fertilization and use of fungicides and growth regulators, and (iv) cultivator tillage with reduced nitrogen fertilization (50%). Bulk soil (BS), as well as root-affected soil (RS), were considered for all treatments in replicates. HT-sequencing of metagenomic DNA yielded approx. 100 Giga bases (Gb) of sequence information. Taxonomic profiling of soil communities revealed the presence of 70 phyla, whereby Proteobacteria, Actinobacteria, Bacteroidetes, Planctomycetes, Acidobacteria, Thaumarchaeota, Firmicutes, Verrucomicrobia and Chloroflexi feature abundances of more than 1%. Functional microbiome profiling uncovered, i.a., numerous potential plant-beneficial, plant-growth-promoting and biocontrol traits predicted to be involved in nutrient provision, phytohormone synthesis, antagonism against pathogens and signal molecule synthesis relevant in microbe-plant interaction. Neither taxonomic nor functional microbiome profiling based on single-read analyses revealed pronounced differences regarding the farming practices applied. Soil metagenome sequences were assembled and taxonomically binned. The ten most reliable and abundant Metagenomically Assembled Genomes (MAGs) were taxonomically classified and metabolically reconstructed. Importance of the phylum Thaumarchaeota for the analyzed microbiome is corroborated by the fact that the four corresponding MAGs were predicted to oxidize ammonia (nitrification), thus contributing to the cycling of nitrogen, and in addition are most probably able to fix carbon dioxide. Moreover, Thaumarchaeota and several bacterial MAGs also possess genes with predicted functions in plant-growth-promotion. Abundances of certain MAGs (species resolution level) responded to the tillage practice, whereas the factors compartment (BS vs. RS) and nitrogen fertilization only marginally shaped MAG abundance profiles. Hence, soil management regimes promoting plant-beneficial microbiome members are very likely advantageous for the respective agrosystem, its health and carbon sequestration and accordingly may enhance plant productivity. Since Chernozem soils are highly fertile, corresponding microbiome data represent a valuable reference resource for agronomy in general.

Metagenomic Composition Analysis of an Ancient Sequenced Polar Bear Jawbone from Svalbard

The sequencing of ancient DNA samples provides a novel way to find, characterize, and distinguish exogenous genomes of endogenous targets. After sequencing, computational composition analysis enables filtering of undesired sources in the focal organism, with the purpose of improving the quality of assemblies and subsequent data analysis. More importantly, such analysis allows extinct and extant species to be identified without requiring a specific or new sequencing run. However, the identification of exogenous organisms is a complex task, given the nature and degradation of the samples, and the evident necessity of using efficient computational tools, which rely on algorithms that are both fast and highly sensitive. In this work, we relied on a fast and highly sensitive tool, FALCON-meta, which measures similarity against whole-genome reference databases, to analyse the metagenomic composition of an ancient polar bear (Ursus maritimus) jawbone fossil. The fossil was collected in Svalbard, Norway, and has an estimated age of 110,000 to 130,000 years. The FASTQ samples contained 349 GB of nonamplified shotgun sequencing data. We identified and localized, relative to the FASTQ samples, the genomes with significant similarities to reference microbial genomes, including those of viruses, bacteria, and archaea, and to fungal, mitochondrial, and plastidial sequences. Among other striking features, we found significant similarities between modern-human, some bacterial and viral sequences (contamination) and the organelle sequences of wild carrot and tomato relative to the whole samples. For each exogenous candidate, we ran a damage pattern analysis, which in addition to revealing shallow levels of damage in the plant candidates, identified the source as contamination.