Investigates how microbial communities in floodplain sediments and riparian soils drive biogeochemical cycling across the East River watershed, using metagenomics and 16S rRNA sequencing to link microbial diversity to hydrological connectivity and landscape position.
Floodplains are the low-lying lands along rivers that flood seasonally, and in mountain watersheds like the East River and Slate River near Crested Butte, Colorado, they act as biogeochemical hotspots — places where water, sediment, plant roots, and microorganisms interact to transform elements such as carbon, nitrogen, sulfur, and iron. The microbes living in these saturated soils control how nutrients are released to streams, how greenhouse gases like methane are produced or consumed, and how metals such as lead or arsenic move through the landscape. Because Gunnison Basin floodplains sit at the headwaters of the Colorado River, what happens in their sediments has consequences far downstream for water quality and for the carbon and nitrogen budgets of the entire region.
Understanding which microbes are present and what they do has been transformed by metagenomics — a technique in which all the DNA in a soil sample is sequenced together and then computationally sorted into draft genomes called metagenome-assembled genomes (MAGs). A MAG is essentially a reconstructed genetic blueprint for a single microbial type, assembled from environmental DNA without ever culturing the organism in a lab. With MAGs in hand, researchers can ask not just "who is there?" but "what metabolic jobs can they perform?" This is essential because most floodplain microbes have never been grown in culture. MAGs also reveal mobile genetic elements — pieces of DNA like plasmids or newly discovered giant linear elements that can shuttle metabolic capabilities between organisms through horizontal gene transfer, the swapping of genes between unrelated microbes.
A few other concepts run through this work. Biogeochemical cycling refers to the movement of elements between organisms and their environment, driven largely by microbial activity such as nitrification (oxidizing ammonia to nitrate) or methane oxidation. Redox transitions describe shifts between oxygen-rich (oxic) and oxygen-poor (anoxic) conditions as the water table rises and falls; these shifts determine which microbes are active and which minerals dissolve or precipitate. Finally, colloidal transport — the movement of nano-sized particles like iron oxides through groundwater — links microbial chemistry to the export of iron, organic matter, and nutrients from floodplain soils into rivers.
Early floodplain microbiome research in the East River watershed established that microbial community structure is tightly tied to landscape position and hydrology. A hillslope-to-riparian transect study found that proximity to the stream, depth to groundwater, and underlying weathered shale strongly shape both who lives where and what they do — with riparian sediments hosting distinct communities capable of carbon and nitrogen fixation, sulfate reduction, and selenium reduction that could affect water quality . Related geomicrobiology work outside the Gunnison Basin, including studies of rock varnish communities, helped develop the culture-independent molecular tools and conceptual frameworks now applied to floodplains .
Genomic bins reconstructed from metagenomic sequencing data representing individual microbial taxa
The cycling of chemical elements between living organisms and the physical environment, particularly carbon and nitrogen cycles
The interplay of belowground, microsite-scale biological, physical, and chemical processes that intersect to drive whole-ecosystem functioning
Transfer of genetic material between organisms through mechanisms other than vertical inheritance
DNA sequences that can move or be transferred between different positions within genomes or between genomes
Clustered viral sequences representing distinct viral taxa based on genomic similarity thresholds
Systematic collection of sediment samples from floodplain sites at regular depth intervals followed by DNA extraction and metagenomic sequencing to ch...
Standard metagenomics pipeline involving assembly with MEGAHIT, binning with multiple tools, dereplication, and taxonomic classification using GTDB-tk...
Standard Illumina library preparation using Nextera XT kit followed by paired-end sequencing on HiSeq 2500 platform to generate comprehensive metageno...
Sediment samples collected in 15-cm increments using 5 cm diameter soil core sampler with slide hammer, then 8.3 cm diameter bucket auger after reachi...
Comprehensive Hidden Markov Model analysis of metagenomic sequences to identify and quantify metabolic functions across biogeochemical pathways. Uses ...
Manual identification and curation of large genome fragments using GC content, coverage patterns, and taxonomic profiles to assemble complete extrachr...
Software for the correlation analysis featured in the above manuscript.
Climate change is driving vegetation shifts in mountain watersheds, with unknown impacts on biogeochemical cycles. We hypothesize that these shifts wi...
Drought is changing the American Mountain West at unprecedented rates with unknown consequences to soil microbiome composition and function. As a part...
This data package comprises analytical results and metadata from stream and groundwater samples collected from the Slate River, East River, Trail Cree...
Microorganisms play a key role in cycling nutrients and contaminants in the terrestrial environment depending on their genetic potential. Here, we pre...
Microorganisms play a key role in cycling nutrients and contaminants in the terrestrial environment depending on their genetic potential. Here, we pre...
Building on the transect approach, a watershed-scale study proposed that meander-bound segments of the floodplain serve as repeatable "scaling motifs" — units whose aggregate microbial capacities can be used to predict biogeochemical function across an entire river corridor (Carnevali et al., 2021). This idea — that complex microbial diversity can be summarized at a tractable landscape scale — has shaped much of the subsequent work in the Gunnison Basin.
A central result is that floodplain soils harbor a small but consistent core microbiome amid enormous overall diversity. Across the upper, middle, and lower East River, roughly one-third of reconstructed genomes were found at all three locations, and about 15% of species recurred across two consecutive years. This core was dominated by Betaproteobacteria with capacities for aerobic respiration, carbon monoxide oxidation, and thiosulfate oxidation, suggesting that despite high turnover, certain metabolic functions are reliably present (Carnevali et al., 2021). Importantly, gene and organism abundance did not predict transcription levels — rare microbes and rare genes can drive disproportionately large biogeochemical effects, with nitrification genes (amoCAB, nxrAB) and methanol/formate oxidation genes among the most highly expressed (Carnevali et al., 2021).
More recent depth-resolved sequencing at the Slate River floodplain near Crested Butte has refined the picture of nitrogen and carbon cycling. The dominant ammonia oxidizers are not the textbook ammonia-oxidizing bacteria — none were recovered — but rather Nitrosotalea-like archaea and comammox bacteria of the Palsa-1315 clade, both adapted to low-ammonia (oligotrophic) conditions and equipped to scavenge alternative nitrogen sources such as urea, cyanate, and biuret (Rasmussen et al., 2025). Methane cycling is similarly unconventional: methanogens appear only in deep anoxic sediments, conventional aerobic methanotrophs are absent, and instead the archaeon Candidatus Methanoperedens and uncultured bacterial lineages handle methane oxidation, sometimes reaching extreme local abundance (Rasmussen et al., 2024). Tied to these Methanoperedens archaea, researchers working in comparable wetland sediments discovered "Borgs" — giant linear extrachromosomal DNA elements up to about 1 million base pairs long that carry genes for redox reactions and energy conservation, potentially expanding the metabolic repertoire of their hosts through horizontal gene transfer (Al-Shayeb et al., 2022) (Dance, 2021).
Microbial activity also controls the fate of metals and iron-bearing particles. In contaminated floodplain soils, seasonal redox cycles repeatedly dissolve and reprecipitate iron oxides and sulfides, yet dissolved lead concentrations remain low because lead binds tightly to particulate organic matter that persists across these transitions (Dewey et al., 2021). In anoxic floodplain groundwater, up to 72% of dissolved iron actually travels as mobile colloids — silicon-coated ferrihydrite nanoparticles complexed with organic matter — that survive both oxic and anoxic conditions and can deliver iron, carbon, and nutrients to surface waters (Engel et al., 2023).
Research published since 2020 has shifted from cataloging community composition toward linking microbial function to hydrology in mechanistic, predictive ways. New work integrates depth-resolved metagenomics with measurements of water flow and chemistry to ask how hydrological transitions — snowmelt, drought, beaver-induced flooding — reshape microbial habitats in real time. A multi-year study of vertical exchange between biogeochemically active surface soils and underlying gravel beds showed that flooding events flush anoxic soil porewater downward into the gravel, while snowmelt and drought restore oxic conditions, with direct consequences for downstream water quality (Babey et al., 2024). Complementing this, machine-learning-based groundwater modeling is now being used to quantify how beaver ponds alter vertical fluxes between soil and gravel bed, with vertical-to-lateral flux ratios increasing roughly tenfold during ponded conditions (Wang et al., 2025).
In parallel, manual curation of large extrachromosomal elements and refined MAG recovery are revealing previously hidden players. The continuing discovery of Borgs, oligotrophic comammox bacteria, and uncultured nitrite oxidizers in the Nitrospirales suggests that even well-studied nitrogen and methane cycles are governed by lineages we are only beginning to recognize (Rasmussen et al., 2025) (Rasmussen et al., 2024).
Major questions remain about how to translate microbial gene content into ecosystem-scale predictions. We still do not know which rare microbes drive disproportionate fluxes of methane, nitrogen, and trace metals, nor how stable the "core" floodplain microbiome will be as snowmelt timing shifts, beaver populations change, and drought intensifies. The biological role of Borgs and other mobile genetic elements — whether they truly expand host metabolism in situ — is unresolved. Linking colloidal iron transport, redox-driven metal cycling, and microbial activity into integrated models of floodplain water quality is a promising frontier, as is understanding how root-regolith interactions and plant communities couple to the subsurface microbiome. Over the next decade, combining long-term metagenomic monitoring with hydrologic modeling and targeted cultivation of unconventional nitrifiers and methanotrophs will be key to forecasting how Gunnison Basin floodplains respond to a warming, more variable climate.
Al-Shayeb, B., et al. (2022). Borgs are giant genetic elements with potential to expand metabolic capacity. Nature. →
Babey, T., et al. (2024). Mountainous floodplain connectivity in response to hydrological transitions. Water Resources Research. →
Carnevali, P. B. M., et al. (2021). Meanders as a scaling motif for understanding of floodplain soil microbiome and biogeochemical potential at the watershed scale. Microbiome. →
Dance, A. (2021). Massive DNA "BORG" Structures Perplex Scientists. Nature. →
Dewey, C., et al. (2021). Porewater Lead Concentrations Limited by Particulate Organic Matter Coupled With Ephemeral Iron(III) and Sulfide Phases during Redox Cycles Within Contaminated Floodplain Soils. Environmental Science & Technology. →
Engel, M., et al. (2023). Structure and composition of natural ferrihydrite nano-colloids in anoxic groundwater. Water Research. →
Lavy, A., et al. (2019). Microbial communities across a hillslope-riparian transect shaped by proximity to the stream, groundwater table, and weathered bedrock. Ecology and Evolution. →
Northup, D. E., et al. (2010). Diversity of rock varnish bacterial communities from Black Canyon, New Mexico. Journal of Geophysical Research: Biogeosciences. →
Parchert, K., et al. (2012). Fungal Communities Associated with Rock Varnish in Black Canyon, New Mexico: Casual Inhabitants or Essential Partners? Geomicrobiology Journal. →
Rasmussen, A. N., et al. (2024). Diverse and unconventional methanogens, methanotrophs, and methylotrophs in metagenome-assembled genomes from subsurface sediments of the Slate River floodplain, Crested Butte, CO, USA. mSystems. →
Rasmussen, A. N., et al. (2025). Metagenome-assembled genomes for oligotrophic nitrifiers From a mountainous gravelbed floodplain. Environmental Microbiology. →
Wang, Y., et al. (2025). Quantifying groundwater response and uncertainty in beaver-influenced mountainous floodplains using machine learning-based model calibration. Water Resources Research. →
Microorganisms play a key role in cycling nutrients and contaminants in the terrestrial environment depending on their genetic potential. Here, we pre...
Microorganisms play a key role in cycling nutrients and contaminants in the terrestrial environment depending on their genetic potential. Here we pres...
This data package includes a time-series of field measurements from May to October 2020 in groundwater and surface water from the Slate River floodpla...
This data package includes a time series of soil sensor data (temperature, water content, bulk electrical conductivity, porewater dissolved oxygen and...
This data package includes processed and undiluted measurements for metal and anion concentrations from pore water (groundwater) samples from the Slat...
This dataset comprises paired environmental and genomic data for soil samples collected across meander-bound floodplains G (ERMG), L (ERML) and Z (ERM...
This dataset comprises paired environmental and genomic data for soil samples collected across meander-bound floodplains G (ERMG), L (ERML) and Z (ERM...
Microorganisms play a key role in cycling nutrients and contaminants in the terrestrial environment depending on their genetic potential. Here we pres...
A Maximum Likelihood phylogenetic tree of rpS3 gene dataset.<br>
87 HMMs for key metabolism enzymes used in the manuscript:<br>Lavy et. al., (2019) Microbial communities across a hillslope-riparian transect shaped b...
Metatranscriptomic data generated from soil collected at the East River watershed, Crested Butte, CO. Normalized by TPM. Annotations included, and lab...
Chemical concentration and soil property measurements taken from soil collected at the East River, Crested Butte, CO.
1660 Clusters of Ribosomal Protein S3 sequences.Files are archived within a ZIP file.<br>From publication:Lavy et al., (2019) Microbial communities ac...
Amino-acid sequences matching HMM search across samples.<br>Supplementary file for publication:Lavy et. al., (2019) Microbial communities across a hil...
Metatranscriptomic data generated from soil collected at the East River watershed, Crested Butte, CO. Normalized by TPM. Annotations included, and lab...
Chemical concentration and soil property measurements taken from soil collected at the East River, Crested Butte, CO.