Quantitative genetic-by-soil microbiome interactions in a perennial grass affect functional traits.

Plants interact with diverse microbiomes that can impact plant growth and performance. Recent studies highlight the potential beneficial aspects of plant microbiomes, including the possibility that microbes facilitate the process of local adaptation in their host plants. Microbially mediated local adaptation in plants occurs when local host genotypes have higher fitness than foreign genotypes because of their affiliation with locally beneficial microbes. Here, plant adaptation results from genetic interactions of the host with locally beneficial microbes (e.g. host genotype-by-microbiome interactions). We used a recombinant inbred line (RIL) mapping population derived from upland and lowland ecotypes of the diploid C4 perennial bunch grass Panicum hallii to explore quantitative genetic responses to soil microbiomes focusing on functional root and shoot traits involved in ecotypic divergence. We show that the growth and development of ecotypes and their trait divergence depends on soil microbiomes. Moreover, we find that the genetic architecture is modified by soil microbiomes, revealing important plant genotype-by-microbiome interactions for quantitative traits. We detected a number of quantitative trait loci (QTL) that interact with the soil microbiome. Our results highlight the importance of microbial interactions in ecotypic divergence and trait genetic architecture in C4 perennial grasses.

EcoPLOT: dynamic analysis of biogeochemical data.

MOTIVATION: We have created EcoPLOT (parameterized linkage of omics-driven technologies), a web-app for the dynamic, interactive analysis of biogeochemical datasets that combines state-of-the-art analysis tools to statistically and graphically explore environmental, geochemical and microbiome datasets. Using the iterative random forest, a machine learning algorithm, EcoPLOT allows for the de novo discovery of drivers which exhibit significant impact on plant, microbial or soil dynamics. AVAILABILITY AND IMPLEMENTATION: EcoPLOT is built entirely within the R language. It can be accessed through any system where R is installed, including Windows, Mac and most Linux systems. EcoPLOT is free to use and can be accessed at https://github.com/cdsanchez18/EcoPLOT.

MeCorS: Metagenome-enabled error correction of single cell sequencing reads.

UNLABELLED: We present a new tool, MeCorS, to correct chimeric reads and sequencing errors in Illumina data generated from single amplified genomes (SAGs). It uses sequence information derived from accompanying metagenome sequencing to accurately correct errors in SAG reads, even from ultra-low coverage regions. In evaluations on real data, we show that MeCorS outperforms BayesHammer, the most widely used state-of-the-art approach. MeCorS performs particularly well in correcting chimeric reads, which greatly improves both accuracy and contiguity of de novo SAG assemblies. AVAILABILITY AND IMPLEMENTATION: https://github.com/metagenomics/MeCorS CONTACT: abremges@cebitec.uni-bielefeld.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Drought re-routes soil microbial carbon metabolism towards emission of volatile metabolites in an artificial tropical rainforest.

Drought impacts on microbial activity can alter soil carbon fate and lead to the loss of stored carbon to the atmosphere as CO2 and volatile organic compounds (VOCs). Here we examined drought impacts on carbon allocation by soil microbes in the Biosphere 2 artificial tropical rainforest by tracking 13C from position-specific 13C-pyruvate into CO2 and VOCs in parallel with multi-omics. During drought, efflux of 13C-enriched acetate, acetone and C4H6O2 (diacetyl) increased. These changes represent increased production and buildup of intermediate metabolites driven by decreased carbon cycling efficiency. Simultaneously,13C-CO2 efflux decreased, driven by a decrease in microbial activity. However, the microbial carbon allocation to energy gain relative to biosynthesis was unchanged, signifying maintained energy demand for biosynthesis of VOCs and other drought-stress-induced pathways. Overall, while carbon loss to the atmosphere via CO2 decreased during drought, carbon loss via efflux of VOCs increased, indicating microbially induced shifts in soil carbon fate.

Genomic potential of Marinobacter aquaeolei, a biogeochemical "opportunitroph".

The genus of Marinobacter is one of the most ubiquitous in the global oceans and assumed to significantly impact various biogeochemical cycles. The genome structure and content of Marinobacter aquaeolei VT8 was analyzed and compared with those from other organisms with diverse adaptive strategies. Here, we report the many "opportunitrophic" genetic characteristics and strategies that M. aquaeolei has adopted to promote survival under various environmental conditions. Genome analysis revealed its metabolic potential to utilize oxygen and nitrate as terminal electron acceptors, iron as an electron donor, and urea, phosphonate, and various hydrocarbons as alternative N, P, and C sources, respectively. Miscellaneous sensory and defense mechanisms, apparently acquired via horizontal gene transfer, are involved in the perception of environmental fluctuations and antibiotic, phage, toxin, and heavy metal resistance, enabling survival under adverse conditions, such as oil-polluted water. Multiple putative integrases, transposases, and plasmids appear to have introduced additional metabolic potential, such as phosphonate degradation. The genomic potential of M. aquaeolei and its similarity to other opportunitrophs are consistent with its cosmopolitan occurrence in diverse environments and highly variable lifestyles.

Impact of Harvest on Switchgrass Leaf Microbial Communities.

Switchgrass is a promising feedstock for biofuel production, with potential for leveraging its native microbial community to increase productivity and resilience to environmental stress. Here, we characterized the bacterial, archaeal and fungal diversity of the leaf microbial community associated with four switchgrass (Panicum virgatum) genotypes, subjected to two harvest treatments (annual harvest and unharvested control), and two fertilization levels (fertilized and unfertilized control), based on 16S rRNA gene and internal transcribed spacer (ITS) region amplicon sequencing. Leaf surface and leaf endosphere bacterial communities were significantly different with Alphaproteobacteria enriched in the leaf surface and Gammaproteobacteria and Bacilli enriched in the leaf endosphere. Harvest treatment significantly shifted presence/absence and abundances of bacterial and fungal leaf surface community members: Gammaproteobacteria were significantly enriched in harvested and Alphaproteobacteria were significantly enriched in unharvested leaf surface communities. These shifts were most prominent in the upland genotype DAC where the leaf surface showed the highest enrichment of Gammaproteobacteria, including taxa with 100% identity to those previously shown to have phytopathogenic function. Fertilization did not have any significant impact on bacterial or fungal communities. We also identified bacterial and fungal taxa present in both the leaf surface and leaf endosphere across all genotypes and treatments. These core taxa were dominated by Methylobacterium, Enterobacteriaceae, and Curtobacterium, in addition to Aureobasidium, Cladosporium, Alternaria and Dothideales. Local core leaf bacterial and fungal taxa represent promising targets for plant microbe engineering and manipulation across various genotypes and harvest treatments. Our study showcases, for the first time, the significant impact that harvest treatment can have on bacterial and fungal taxa inhabiting switchgrass leaves and the need to include this factor in future plant microbial community studies.

Conservation of Endophyte Bacterial Community Structure Across Two Panicum Grass Species.

Panicum represents a large genus of many North American prairie grass species. These include switchgrass (Panicum virgatum), a biofuel crop candidate with wide geographic range, as well as Panicum hallii, a close relative to switchgrass, which serves as a model system for the study of Panicum genetics due to its diploid genome and short growth cycles. For the advancement of switchgrass as a biofuel crop, it is essential to understand host microbiome interactions, which can be impacted by plant genetics and environmental factors inducing ecotype-specific phenotypic traits. We here compared rhizosphere and root endosphere bacterial communities of upland and lowland P. virgatum and P. hallii genotypes planted at two sites in Texas. Our analysis shows that sampling site predominantly contributed to bacterial community variance in the rhizosphere, however, impacted root endosphere bacterial communities much less. Instead we observed a relatively large core endophytic microbiome dominated by ubiquitously root-colonizing bacterial genera Streptomyces, Pseudomonas, and Bradyrhizobium. Endosphere communities displayed comparable diversity and conserved community structures across genotypes of both Panicum species. Functional insights into interactions between P. hallii and its root endophyte microbiome could hence inform testable hypotheses that are relevant for the improvement of switchgrass as a biofuel crop.

Plant compartment and genetic variation drive microbiome composition in switchgrass roots.

Switchgrass (Panicum virgatum) is a promising biofuel crop native to the United States with genotypes that are adapted to a wide range of distinct ecosystems. Various plants have been shown to undergo symbioses with plant growth-promoting bacteria and fungi, however, plant-associated microbial communities of switchgrass have not been extensively studied to date. We present 16S ribosomal RNA gene and internal transcribed spacer (ITS) data of rhizosphere and root endosphere compartments of four switchgrass genotypes to test the hypothesis that host selection of its root microbiota prevails after transfer to non-native soil. We show that differences in bacterial, archaeal and fungal community composition and diversity are strongly driven by plant compartment and switchgrass genotypes and ecotypes. Plant-associated microbiota show an enrichment in Alphaproteobacteria and Actinobacteria as well as Sordariales and Pleosporales compared with the surrounding soil. Root associated compartments display low-complexity communities dominated and enriched in Actinobacteria, in particular Streptomyces, in the lowland genotypes, and in Alphaproteobacteria, specifically Sphingobium, in the upland genotypes. Our comprehensive root analysis serves as a snapshot of host-specific bacterial and fungal associations of switchgrass in the field and confirms that host-selected microbiomes persist after transfer to non-native soil.

Soil exchange rates of COS and CO18O differ with the diversity of microbial communities and their carbonic anhydrase enzymes.

Differentiating the contributions of photosynthesis and respiration to the global carbon cycle is critical for improving predictive climate models. Carbonic anhydrase (CA) activity in leaves is responsible for the largest biosphere-atmosphere trace gas fluxes of carbonyl sulfide (COS) and the oxygen-18 isotopologue of carbon dioxide (CO18O) that both reflect gross photosynthetic rates. However, CA activity also occurs in soils and will be a source of uncertainty in the use of COS and CO18O as carbon cycle tracers until process-based constraints are improved. In this study, we measured COS and CO18O exchange rates and estimated the corresponding CA activity in soils from a range of biomes and land use types. Soil CA activity was not uniform for COS and CO2, and patterns of divergence were related to microbial community composition and CA gene expression patterns. In some cases, the same microbial taxa and CA classes catalyzed both COS and CO2 reactions in soil, but in other cases the specificity towards the two substrates differed markedly. CA activity for COS was related to fungal taxa and ß-D-CA expression, whereas CA activity for CO2 was related to algal and bacterial taxa and α-CA expression. This study integrates gas exchange measurements, enzyme activity models, and characterization of soil taxonomic and genetic diversity to build connections between CA activity and the soil microbiome. Importantly, our results identify kinetic parameters to represent soil CA activity during application of COS and CO18O as carbon cycle tracers.

Shotgun metagenome data of a defined mock community using Oxford Nanopore, PacBio and Illumina technologies.

Metagenomic sequence data from defined mock communities is crucial for the assessment of sequencing platform performance and downstream analyses, including assembly, binning and taxonomic assignment. We report a comparison of shotgun metagenome sequencing and assembly metrics of a defined microbial mock community using the Oxford Nanopore Technologies (ONT) MinION, PacBio and Illumina sequencing platforms. Our synthetic microbial community BMock12 consists of 12 bacterial strains with genome sizes spanning 3.2-7.2 Mbp, 40-73% GC content, and 1.5-7.3% repeats. Size selection of both PacBio and ONT sequencing libraries prior to sequencing was essential to yield comparable relative abundances of organisms among all sequencing technologies. While the Illumina-based metagenome assembly yielded good coverage with few misassemblies, contiguity was greatly improved by both, Illumina + ONT and Illumina + PacBio hybrid assemblies but increased misassemblies, most notably in genomes with high sequence similarity to each other. Our resulting datasets allow evaluation and benchmarking of bioinformatics software on Illumina, PacBio and ONT platforms in parallel.

Capturing the genetic makeup of the active microbiome in situ.

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

BBMerge - Accurate paired shotgun read merging via overlap.

Merging paired-end shotgun reads generated on high-throughput sequencing platforms can substantially improve various subsequent bioinformatics processes, including genome assembly, binning, mapping, annotation, and clustering for taxonomic analysis. With the inexorable growth of sequence data volume and CPU core counts, the speed and scalability of read-processing tools becomes ever-more important. The accuracy of shotgun read merging is crucial as well, as errors introduced by incorrect merging percolate through to reduce the quality of downstream analysis. Thus, we designed a new tool to maximize accuracy and minimize processing time, allowing the use of read merging on larger datasets, and in analyses highly sensitive to errors. We present BBMerge, a new merging tool for paired-end shotgun sequence data. We benchmark BBMerge by comparison with eight other widely used merging tools, assessing speed, accuracy and scalability. Evaluations of both synthetic and real-world datasets demonstrate that BBMerge produces merged shotgun reads with greater accuracy and at higher speed than any existing merging tool examined. BBMerge also provides the ability to merge non-overlapping shotgun read pairs by using k-mer frequency information to assemble the unsequenced gap between reads, achieving a significantly higher merge rate while maintaining or increasing accuracy.

High-resolution phylogenetic microbial community profiling.

Over the past decade, high-throughput short-read 16S rRNA gene amplicon sequencing has eclipsed clone-dependent long-read Sanger sequencing for microbial community profiling. The transition to new technologies has provided more quantitative information at the expense of taxonomic resolution with implications for inferring metabolic traits in various ecosystems. We applied single-molecule real-time sequencing for microbial community profiling, generating full-length 16S rRNA gene sequences at high throughput, which we propose to name PhyloTags. We benchmarked and validated this approach using a defined microbial community. When further applied to samples from the water column of meromictic Sakinaw Lake, we show that while community structures at the phylum level are comparable between PhyloTags and Illumina V4 16S rRNA gene sequences (iTags), variance increases with community complexity at greater water depths. PhyloTags moreover allowed less ambiguous classification. Last, a platform-independent comparison of PhyloTags and in silico generated partial 16S rRNA gene sequences demonstrated significant differences in community structure and phylogenetic resolution across multiple taxonomic levels, including a severe underestimation in the abundance of specific microbial genera involved in nitrogen and methane cycling across the Lake's water column. Thus, PhyloTags provide a reliable adjunct or alternative to cost-effective iTags, enabling more accurate phylogenetic resolution of microbial communities and predictions on their metabolic potential.

Next generation sequencing data of a defined microbial mock community.

Generating sequence data of a defined community composed of organisms with complete reference genomes is indispensable for the benchmarking of new genome sequence analysis methods, including assembly and binning tools. Moreover the validation of new sequencing library protocols and platforms to assess critical components such as sequencing errors and biases relies on such datasets. We here report the next generation metagenomic sequence data of a defined mock community (Mock Bacteria ARchaea Community; MBARC-26), composed of 23 bacterial and 3 archaeal strains with finished genomes. These strains span 10 phyla and 14 classes, a range of GC contents, genome sizes, repeat content and encompass a diverse abundance profile. Short read Illumina and long-read PacBio SMRT sequences of this mock community are described. These data represent a valuable resource for the scientific community, enabling extensive benchmarking and comparative evaluation of bioinformatics tools without the need to simulate data. As such, these data can aid in improving our current sequence data analysis toolkit and spur interest in the development of new tools.

Similar Microbial Communities Found on Two Distant Seafloor Basalts.

The oceanic crust forms two thirds of the Earth's surface and hosts a large phylogenetic and functional diversity of microorganisms. While advances have been made in the sedimentary realm, our understanding of the igneous rock portion as a microbial habitat has remained limited. We present the first comparative metagenomic microbial community analysis from ocean floor basalt environments at the Lo'ihi Seamount, Hawai'i, and the East Pacific Rise (EPR; 9°N). Phylogenetic analysis indicates the presence of a total of 43 bacterial and archaeal mono-phyletic groups, dominated by Alpha- and Gammaproteobacteria, as well as Thaumarchaeota. Functional gene analysis suggests that these Thaumarchaeota play an important role in ammonium oxidation on seafloor basalts. In addition to ammonium oxidation, the seafloor basalt habitat reveals a wide spectrum of other metabolic potentials, including CO2 fixation, denitrification, dissimilatory sulfate reduction, and sulfur oxidation. Basalt communities from Lo'ihi and the EPR show considerable metabolic and phylogenetic overlap down to the genus level despite geographic distance and slightly different seafloor basalt mineralogy.

Metagenomic insights into the dominant Fe(II) oxidizing Zetaproteobacteria from an iron mat at Lo´ihi, Hawai´l.

Zetaproteobacteria are among the most prevalent Fe(II)-oxidizing bacteria (FeOB) at deep-sea hydrothermal vents; however, knowledge about their environmental significance is limited. We provide metagenomic insights into an iron mat at the Lo´ihi Seamount, Hawai´l, revealing novel genomic information of locally dominant Zetaproteobacteria lineages. These lineages were previously estimated to account for ~13% of all local Zetaproteobacteria based on 16S clone library data. Biogeochemically relevant genes include nitrite reductases, which were previously not identified in Zetaproteobacteria, sulfide:quinone oxidases, and ribulose-1,5-bisphosphate carboxylase (RuBisCo). Genes assumed to be involved in Fe(II) oxidation correlate in synteny and share 87% amino acid similarity with those previously identified in the related Zetaproteobacterium Mariprofundus ferrooxydans PV-1. Overall, Zetaproteobacteria genes appear to originate primarily from within the Proteobacteria and the Fe(II)-oxidizing Leptospirillum spp. and are predicted to facilitate adaptation to a deep-sea hydrothermal vent environment in addition to microaerophilic Fe(II) and H2S oxidation. This dataset represents the first metagenomic study of FeOB from an iron oxide mat at a deep-sea hydrothermal habitat.

Mariprofundus ferrooxydans PV-1 the first genome of a marine Fe(II) oxidizing Zetaproteobacterium.

Mariprofundus ferrooxydans PV-1 has provided the first genome of the recently discovered Zetaproteobacteria subdivision. Genome analysis reveals a complete TCA cycle, the ability to fix CO(2), carbon-storage proteins and a sugar phosphotransferase system (PTS). The latter could facilitate the transport of carbohydrates across the cell membrane and possibly aid in stalk formation, a matrix composed of exopolymers and/or exopolysaccharides, which is used to store oxidized iron minerals outside the cell. Two-component signal transduction system genes, including histidine kinases, GGDEF domain genes, and response regulators containing CheY-like receivers, are abundant and widely distributed across the genome. Most of these are located in close proximity to genes required for cell division, phosphate uptake and transport, exopolymer and heavy metal secretion, flagellar biosynthesis and pilus assembly suggesting that these functions are highly regulated. Similar to many other motile, microaerophilic bacteria, genes encoding aerotaxis as well as antioxidant functionality (e.g., superoxide dismutases and peroxidases) are predicted to sense and respond to oxygen gradients, as would be required to maintain cellular redox balance in the specialized habitat where M. ferrooxydans resides. Comparative genomics with other Fe(II) oxidizing bacteria residing in freshwater and marine environments revealed similar content, synteny, and amino acid similarity of coding sequences potentially involved in Fe(II) oxidation, signal transduction and response regulation, oxygen sensation and detoxification, and heavy metal resistance. This study has provided novel insights into the molecular nature of Zetaproteobacteria.