Meanders as a scaling motif for understanding of floodplain soil microbiome and biogeochemical potential at the watershed scale

Background Biogeochemical exports from watersheds are modulated by the activity of microorganisms that function over micron scales. Here, we tested the hypothesis that meander-bound regions share a core microbiome and exhibit patterns of metabolic potential that broadly predict biogeochemical processes in floodplain soils along a river corridor. Results We intensively sampled the microbiomes of floodplain soils located in the upper, middle, and lower reaches of the East River, Colorado. Despite the very high microbial diversity and complexity of the soils, we reconstructed 248 quality draft genomes representative of subspecies. Approximately one third of these bacterial subspecies was detected across all three locations at similar abundance levels, and ~ 15% of species were detected in two consecutive years. Within the meander-bound floodplains, we did not detect systematic patterns of gene abundance based on sampling position relative to the river. However, across meanders, we identified a core floodplain microbiome that is enriched in capacities for aerobic respiration, aerobic CO oxidation, and thiosulfate oxidation with the formation of elemental sulfur. Given this, we conducted a transcriptomic analysis of the middle floodplain. In contrast to predictions made based on the prominence of gene inventories, the most highly transcribed genes were relatively rare amoCAB and nxrAB (for nitrification) genes, followed by genes involved in methanol and formate oxidation, and nitrogen and CO 2 fixation. Within all three meanders, low soil organic carbon correlated with high activity of genes involved in methanol, formate, sulfide, hydrogen, and ammonia oxidation, nitrite oxidoreduction, and nitrate and nitrite reduction. Overall, the results emphasize the importance of sulfur, one-carbon and nitrogen compound metabolism in soils of the riparian corridor. Conclusions The disparity between the scale of a microbial cell and the scale of a watershed currently limits the development of genomically informed predictive models describing watershed biogeochemical function. Meander-bound floodplains appear to serve as scaling motifs that predict aggregate capacities for biogeochemical transformations, providing a foundation for incorporating riparian soil microbiomes in watershed models. Widely represented genetic capacities did not predict in situ activity at one time point, but rather they define a reservoir of biogeochemical potential available as conditions change. Video abstract

[1]  Pedro M. Valero-Mora,et al.  ggplot2: Elegant Graphics for Data Analysis , 2010 .

[2]  G. King,et al.  Distribution, diversity and ecology of aerobic CO-oxidizing bacteria , 2007, Nature Reviews Microbiology.

[3]  T. Bolch,et al.  Importance and vulnerability of the world’s water towers , 2019, Nature.

[4]  J. Banfield,et al.  Consistent Metagenome-Derived Metrics Verify and Delineate Bacterial Species Boundaries , 2020, mSystems.

[5]  R. Macarthur The Problem of Pattern and Scale in Ecology: The Robert H. MacArthur Award Lecture , 2005 .

[6]  Pedro M. Coutinho,et al.  The carbohydrate-active enzymes database (CAZy) in 2013 , 2013, Nucleic Acids Res..

[7]  G. Likens,et al.  Complex response of the forest nitrogen cycle to climate change , 2012, Proceedings of the National Academy of Sciences.

[8]  James L. Melville The Uniform Manifold Approximation and Projection (UMAP) Method for Dimensionality Reduction [R package uwot version 0.1.10] , 2020 .

[9]  K. Williams,et al.  Taxonomically and metabolically distinct microbial communities with depth and across a hillslope to riparian zone transect , 2019, bioRxiv.

[10]  Anders F. Andersson,et al.  Binning metagenomic contigs by coverage and composition , 2014, Nature Methods.

[11]  Hadley Wickham,et al.  ggplot2 - Elegant Graphics for Data Analysis (2nd Edition) , 2017 .

[12]  V. Caron,et al.  United states. , 2018, Nursing standard (Royal College of Nursing (Great Britain) : 1987).

[13]  D. Dowrick,et al.  Water table elevation controls on soil nitrogen cycling in riparian wetlands along a European climatic gradient , 2004 .

[14]  Sarah L. Westcott,et al.  De novo clustering methods outperform reference-based methods for assigning 16S rRNA gene sequences to operational taxonomic units , 2015, PeerJ.

[15]  The UniProt Consortium,et al.  UniProt: a worldwide hub of protein knowledge , 2018, Nucleic Acids Res..

[16]  M. Wiley,et al.  Microbial Community Functional Potential and Composition Are Shaped by Hydrologic Connectivity in Riverine Floodplain Soils , 2016, Microbial Ecology.

[17]  J. Banfield,et al.  The East River, Colorado, Watershed: A Mountainous Community Testbed for Improving Predictive Understanding of Multiscale Hydrological–Biogeochemical Dynamics , 2018 .

[18]  K. Williams,et al.  Predicting sedimentary bedrock subsurface weathering fronts and weathering rates , 2019, Scientific Reports.

[19]  K. Williams,et al.  LiDAR collection in August 2015 over the East River Watershed, Colorado, USA , 2017 .

[20]  Itai Sharon,et al.  Hydrogen-based metabolism as an ancestral trait in lineages sibling to the Cyanobacteria , 2019, Nature Communications.

[21]  Hilde van der Togt,et al.  Publisher's Note , 2003, J. Netw. Comput. Appl..

[22]  Marcel Brun,et al.  HMMER Cut-off Threshold Tool (HMMERCTTER): Supervised classification of superfamily protein sequences with a reliable cut-off threshold , 2017, bioRxiv.

[23]  J. Banfield,et al.  dRep: a tool for fast and accurate genomic comparisons that enables improved genome recovery from metagenomes through de-replication , 2017, The ISME Journal.

[24]  T. Rattei,et al.  Phylogenetic and environmental diversity of DsrAB-type dissimilatory (bi)sulfite reductases , 2014, The ISME Journal.

[25]  L. Hug,et al.  Aquatic Elusimicrobia are metabolically diverse compared to gut microbiome Elusimicrobia and some have novel nitrogenase-like gene clusters , 2019, bioRxiv.

[26]  K. Eriksson,et al.  Methanol formation during lignin degradation by Phanerochaete chrysosporium , 2004, Applied Microbiology and Biotechnology.

[27]  Brian C. Thomas,et al.  Unusual biology across a group comprising more than 15% of domain Bacteria , 2015, Nature.

[28]  R. Paine,et al.  Ecological Determinism in the Competition for Space: The Robert H. MacArthur Award Lecture , 1984 .

[29]  Dongwan D. Kang,et al.  MetaBAT, an efficient tool for accurately reconstructing single genomes from complex microbial communities , 2015, PeerJ.

[30]  G. King,et al.  Atmospheric carbon monoxide oxidation is a widespread mechanism supporting microbial survival , 2019, The ISME Journal.

[31]  Thomas K. F. Wong,et al.  ModelFinder: Fast Model Selection for Accurate Phylogenetic Estimates , 2017, Nature Methods.

[32]  Hiroyuki Ogata,et al.  KEGG: Kyoto Encyclopedia of Genes and Genomes , 1999, Nucleic Acids Res..

[33]  Microbial communities across a hillslope‐riparian transect shaped by proximity to the stream, groundwater table, and weathered bedrock , 2019, Ecology and evolution.

[34]  R. Maxwell,et al.  Snowmelt controls on concentration‐discharge relationships and the balance of oxidative and acid‐base weathering fluxes in an alpine catchment, East River, Colorado , 2017 .

[35]  K. Williams,et al.  Geochemical Exports to River From the Intrameander Hyporheic Zone Under Transient Hydrologic Conditions: East River Mountainous Watershed, Colorado , 2018, Water Resources Research.

[36]  J. Murrell,et al.  XoxF encoding an alternative methanol dehydrogenase is widespread in coastal marine environments. , 2015, Environmental microbiology.

[37]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[38]  K. Williams,et al.  Contrasting the hydrologic response due to land cover and climate change in a mountain headwaters system , 2016 .

[39]  A. von Haeseler,et al.  UFBoot2: Improving the Ultrafast Bootstrap Approximation , 2017, bioRxiv.

[40]  Kazutaka Katoh,et al.  A simple method to control over-alignment in the MAFFT multiple sequence alignment program , 2016, Bioinform..

[41]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[42]  M. Meybeck,et al.  Mountains of the world, water towers for humanity: Typology, mapping, and global significance , 2007 .

[43]  Brandi L. Cantarel,et al.  The Carbohydrate-Active EnZymes database (CAZy): an expert resource for Glycogenomics , 2008, Nucleic Acids Res..

[44]  Anne-Béatrice Dufour,et al.  The ade4 Package: Implementing the Duality Diagram for Ecologists , 2007 .

[45]  Wilfried Thuiller,et al.  Combining the fourth-corner and the RLQ methods for assessing trait responses to environmental variation. , 2014, Ecology.

[46]  W. Huber,et al.  Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 , 2014, Genome Biology.

[47]  C. Guertin,et al.  Land-use influences the distribution and activity of high affinity CO-oxidizing bacteria associated to type I-coxL genotype in soil , 2014, Front. Microbiol..

[48]  K. Williams,et al.  Investigating Microtopographic and Soil Controls on a Mountainous Meadow Plant Community Using High‐Resolution Remote Sensing and Surface Geophysical Data , 2019, Journal of Geophysical Research: Biogeosciences.

[49]  Andrew C. Tolonen,et al.  Functional Diversity of Carbohydrate-Active Enzymes Enabling a Bacterium to Ferment Plant Biomass , 2014, PLoS genetics.

[50]  B. Keel,et al.  Comparison of Burrows-Wheeler Transform-Based Mapping Algorithms Used in High-Throughput Whole-Genome Sequencing: Application to Illumina Data for Livestock Genomes1 , 2018, Front. Genet..

[51]  Donovan H. Parks,et al.  A complete domain-to-species taxonomy for Bacteria and Archaea , 2020, Nature Biotechnology.

[52]  Toni Gabaldón,et al.  trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses , 2009, Bioinform..

[53]  Zhenglu Yang,et al.  dbCAN2: a meta server for automated carbohydrate-active enzyme annotation , 2018, Nucleic Acids Res..

[54]  D. Dušková,et al.  Fermentation of pectin and glucose, and activity of pectin‐degrading enzymes in the rumen bacterium Lachnospira multiparus , 2001, Letters in applied microbiology.

[55]  P. Hugenholtz,et al.  Genome sequences of rare, uncultured bacteria obtained by differential coverage binning of multiple metagenomes , 2013, Nature Biotechnology.

[56]  Cindy J. Castelle,et al.  The rise of diversity in metabolic platforms across the Candidate Phyla Radiation , 2019, bioRxiv.

[57]  Huub J. M. Op den Camp,et al.  PQQ-dependent methanol dehydrogenases: rare-earth elements make a difference , 2014, Applied Microbiology and Biotechnology.

[58]  Siu-Ming Yiu,et al.  IDBA-UD: a de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth , 2012, Bioinform..

[59]  Alexander J Probst,et al.  Recovery of genomes from metagenomes via a dereplication, aggregation and scoring strategy , 2017, Nature Microbiology.

[60]  R. Dahlgren,et al.  Convergent evidence for widespread rock nitrogen sources in Earth’s surface environment , 2018, Science.

[61]  Stéphane Dray,et al.  Improved testing of species traits--environment relationships in the fourth-corner problem. , 2012, Ecology.

[62]  Stephen D. Sebestyen,et al.  Sources, transformations, and hydrological processes that control stream nitrate and dissolved organic matter concentrations during snowmelt in an upland forest , 2008 .

[63]  P. Legendre,et al.  vegan : Community Ecology Package. R package version 1.8-5 , 2007 .

[64]  R. Irizarry ggplot2 , 2019, Introduction to Data Science.

[65]  Blake A. Simmons,et al.  MaxBin 2.0: an automated binning algorithm to recover genomes from multiple metagenomic datasets , 2016, Bioinform..

[66]  K. Williams,et al.  Integrated imaging of above and below ground properties and their interactions: A case study in East River Watershed, Colorado , 2018 .

[67]  Brian C. Thomas,et al.  Mediterranean grassland soil C–N compound turnover is dependent on rainfall and depth, and is mediated by genomically divergent microorganisms , 2019, Nature Microbiology.

[68]  Mark A. Miller,et al.  Creating the CIPRES Science Gateway for inference of large phylogenetic trees , 2010, 2010 Gateway Computing Environments Workshop (GCE).

[69]  P. Nico,et al.  Shale as a Source of Organic Carbon in Floodplain Sediments of a Mountainous Watershed , 2020, Journal of Geophysical Research: Biogeosciences.

[70]  Brian C. Thomas,et al.  Thousands of microbial genomes shed light on interconnected biogeochemical processes in an aquifer system , 2016, Nature Communications.

[71]  P. Boeckx,et al.  Assessment of the importance of dissimilatory nitrate reduction to ammonium for the terrestrial nitrogen cycle , 2011 .

[72]  Miriam L. Land,et al.  Trace: Tennessee Research and Creative Exchange Prodigal: Prokaryotic Gene Recognition and Translation Initiation Site Identification Recommended Citation Prodigal: Prokaryotic Gene Recognition and Translation Initiation Site Identification , 2022 .

[73]  Hiroyuki Ogata,et al.  KofamKOALA: KEGG Ortholog assignment based on profile HMM and adaptive score threshold , 2019, bioRxiv.

[74]  P. Legendre,et al.  Associations between species and groups of sites: indices and statistical inference. , 2009, Ecology.

[75]  S. Levin The problem of pattern and scale in ecology , 1992 .