Search

Setting the pace of life and constraining the role of members in food webs, body size can affect the structure and dynamics of communities across multiple scales of biological organization (e.g., from the individual to the ecosystem). However, its effects on shaping microbial communities, as well as underlying assembly processes, remain poorly known. Here, we analyzed microbial diversity in the largest urban lake in China and disentangled the ecological processes governing microbial eukaryotes and prokaryotes using 16S and 18S amplicon sequencing. We found that pico/nano-eukaryotes (0.22−20 μm) and micro-eukaryotes (20−200 μm) showed significant differences in terms of both community composition and assembly processes even though they were characterized by similar phylotype diversity. We also found scale dependencies whereby micro-eukaryotes were strongly governed by environmental selection at the local scale and dispersal limitation at the regional scale. Interestingly, it was the micro-eukaryotes, rather than the pico/nano-eukaryotes, that shared similar distribution and community assembly patterns with the prokaryotes. This indicated that assembly processes of eukaryotes may be coupled or decoupled from prokaryotes’ assembly processes based on eukaryote cell size. While the results support the important influence of cell size, there may be other factors leading to different levels of assembly process coupling across size classes. Additional studies are needed to quantitatively parse the influence of cell size versus other factors as drivers of coordinated and divergent community assembly processes across microbial groups. Regardless of the governing mechanisms, our results show that there are clear patterns in how assembly processes are coupled across sub-communities defined by cell size. These size-structured patterns could be used to help predict shifts in microbial food webs in response to future disturbance.

ABSTRACT Arctic permafrost is thawing due to global warming, with unknown consequences on the microbial inhabitants or associated viruses. DNA viruses have previously been shown to be abundant and active in thawing permafrost, but little is known about RNA viruses in these systems. To address this knowledge gap, we assessed the composition of RNA viruses in thawed permafrost samples that were incubated for 97 days at 4°C to simulate thaw conditions. A diverse RNA viral community was assembled from metatranscriptome data including double-stranded RNA viruses, dominated by Reoviridae and Hypoviridae, and negative and positive single-stranded RNA viruses, with relatively high representations of Rhabdoviridae and Leviviridae, respectively. Sequences corresponding to potential plant and human pathogens were also detected. The detected RNA viruses primarily targeted dominant eukaryotic taxa in the samples (e.g., fungi, Metazoa and Viridiplantae) and the viral community structures were significantly associated with predicted host populations. These results indicate that RNA viruses are linked to eukaryotic host dynamics. Several of the RNA viral sequences contained auxiliary metabolic genes encoding proteins involved in carbon utilization (e.g., polygalacturosase), implying their potential roles in carbon cycling in thawed permafrost. IMPORTANCE Permafrost is thawing at a rapid pace in the Arctic with largely unknown consequences on ecological processes that are fundamental to Arctic ecosystems. This is the first study to determine the composition of RNA viruses in thawed permafrost. Other recent studies have characterized DNA viruses in thawing permafrost, but the majority of DNA viruses are bacteriophages that target bacterial hosts. By contrast RNA viruses primarily target eukaryotic hosts and thus represent potential pathogenic threats to humans, animals, and plants. Here, we find that RNA viruses in permafrost are novel and distinct from those in other habitats studied to date. The COVID-19 pandemic has heightened awareness of the importance of potential environmental reservoirs of emerging RNA viral pathogens. We demonstrate that some potential pathogens were detected after an experimental thawing regime. These results are important for understanding critical viral-host interactions and provide a better understanding of the ecological roles that RNA viruses play as permafrost thaws.

ABSTRACT Rivers have a significant role in global carbon and nitrogen cycles, serving as a nexus for nutrient transport between terrestrial and marine ecosystems. Although rivers have a small global surface area, they contribute substantially to worldwide greenhouse gas emissions through microbially mediated processes within the river hyporheic zone. Despite this importance, research linking microbial and viral communities to specific biogeochemical reactions is still nascent in these sediment environments. To survey the metabolic potential and gene expression underpinning carbon and nitrogen biogeochemical cycling in river sediments, we collected an integrated data set of 33 metagenomes, metaproteomes, and paired metabolomes. We reconstructed over 500 microbial metagenome-assembled genomes (MAGs), which we dereplicated into 55 unique, nearly complete medium- and high-quality MAGs spanning 12 bacterial and archaeal phyla. We also reconstructed 2,482 viral genomic contigs, which were dereplicated into 111 viral MAGs (vMAGs) of >10 kb in size. As a result of integrating gene expression data with geochemical and metabolite data, we created a conceptual model that uncovered new roles for microorganisms in organic matter decomposition, carbon sequestration, nitrogen mineralization, nitrification, and denitrification. We show how these metabolic pathways, integrated through shared resource pools of ammonium, carbon dioxide, and inorganic nitrogen, could ultimately contribute to carbon dioxide and nitrous oxide fluxes from hyporheic sediments. Further, by linking viral MAGs to these active microbial hosts, we provide some of the first insights into viral modulation of river sediment carbon and nitrogen cycling. IMPORTANCE Here we created HUM-V (hyporheic uncultured microbial and viral), an annotated microbial and viral MAG catalog that captures strain and functional diversity encoded in these Columbia River sediment samples. Demonstrating its utility, this genomic inventory encompasses multiple representatives of dominant microbial and archaeal phyla reported in other river sediments and provides novel viral MAGs that can putatively infect these. Furthermore, we used HUM-V to recruit gene expression data to decipher the functional activities of these MAGs and reconstruct their active roles in Columbia River sediment biogeochemical cycling. Ultimately, we show the power of MAG-resolved multi-omics to uncover interactions and chemical handoffs in river sediments that shape an intertwined carbon and nitrogen metabolic network. The accessible microbial and viral MAGs in HUM-V will serve as a community resource to further advance more untargeted, activity-based measurements in these, and related, freshwater terrestrial-aquatic ecosystems.

Understanding the mechanisms underlying the assembly of communities has long been the goal of many ecological studies. While several studies have evaluated community wide ecological assembly, fewer have focused on investigating the impacts of individual members within a community or assemblage on ecological assembly. Here, we adapted a previous null model β-nearest taxon index (βNTI) to measure the contribution of individual features within an ecological community to overall assembly. This new metric, called feature-level βNTI (βNTIfeat), enables researchers to determine whether ecological features (e.g., individual microbial taxa) contribute to divergence, convergence, or have insignificant impacts across spatiotemporally resolved metacommunities or meta-assemblages. Using βNTIfeat, we revealed that unclassified microbial lineages often contributed to community divergence while diverse groups (e.g., Crenarchaeota, Alphaproteobacteria, and Gammaproteobacteria) contributed to convergence. We also demonstrate that βNTIfeat can be extended to other ecological assemblages such as organic molecules comprising organic matter (OM) pools. OM had more inconsistent trends compared to the microbial community though CHO-containing molecular formulas often contributed to convergence, while nitrogen and phosphorus-containing formulas contributed to both convergence and divergence. A network analysis was used to relate βNTIfeat values from the putatively active microbial community and the OM assemblage and examine potentially common contributions to ecological assembly across different communities/assemblages. This analysis revealed that P-containing formulas often contributed to convergence/divergence separately from other ecological features and N-containing formulas often contributed to assembly in coordination with microorganisms. Additionally, members of Family Geobacteraceae were often observed to contribute to convergence/divergence in conjunction with both N- and P-containing formulas, suggesting a coordinated ecological role for family members and the nitrogen/phosphorus cycle. Overall, we show that βNTIfeat offers opportunities to investigate the community or assemblage members, which shape the phylogenetic or functional landscape, and demonstrate the potential to evaluate potential points of coordination across various community types.

Microbial communities in incipient soil systems serve as the only biotic force shaping landscape evolution. However, the underlying ecological forces shaping microbial community structure and function are inadequately understood. We used amplicon sequencing to determine microbial taxonomic assembly and metagenome sequencing to evaluate microbial functional assembly in incipient basaltic soil subjected to precipitation. Community composition was stratified with soil depth in the pre-precipitation samples, with surficial communities maintaining their distinct structure and diversity after precipitation, while the deeper soil samples appeared to become more uniform. The structural community assembly remained deterministic in pre- and post-precipitation periods, with homogenous selection being dominant. Metagenome analysis revealed that carbon and nitrogen functional potential was assembled stochastically. Sub-populations putatively involved in the nitrogen cycle and carbon fixation experienced counteracting assembly pressures at the deepest depths, suggesting the communities may functionally assemble to respond to short-term environmental fluctuations and impact the landscape-scale response to perturbations. We propose that contrasting assembly forces impact microbial structure and potential function in an incipient landscape; in situ landscape characteristics (here homogenous parent material) drive community structure assembly, while short-term environmental fluctuations (here precipitation) shape environmental variations that are random in the soil depth profile and drive stochastic sub-population functional dynamics.

The initial assembly of the plant microbiome may establish the trajectory of forthcoming microbiome states, which could determine the overall future health of the plant. However, while much is known about the initial microbiome assembly of grasses and agricultural crops, less is known about the initial microbiome of long-lived trees, such as poplar ( Populus

Predictive biogeochemical modeling requires data-model integration that enables explicit representation of the sophisticated roles of microbial processes that transform substrates. Data from high-resolution organic matter (OM) characterization are increasingly available and can serve as a critical resource for this purpose, but their incorporation into biogeochemical models is often prohibited due to an over-simplified description of reaction networks. To fill this gap, we proposed a new concept of biogeochemical modeling—termed substrate-explicit modeling—that enables parameterizing OM-specific oxidative degradation pathways and reaction rates based on the thermodynamic properties of OM pools. Based on previous developments in the literature, we characterized the resulting kinetic models by only two parameters regardless of the complexity of OM profiles, which can greatly facilitate the integration with reactive transport models for ecosystem simulations by alleviating the difficulty in parameter identification. The two parameters include maximal growth rate (μmax) and harvest volume (Vh) (i.e., the volume that a microbe can access for harvesting energy). For every detected organic molecule in a given sample, our approach provides a systematic way to formulate reaction kinetics from chemical formula, which enables the evaluation of the impact of OM character on biogeochemical processes across conditions. In a case study of two sites with distinct OM thermodynamics using ultra high-resolution metabolomics datasets derived from Fourier transform ion cyclotron resonance mass spectrometry analyses, our method predicted how oxidative degradation is primarily driven by thermodynamic efficiency of OM consistent with experimental rate measurements (as shown by correlation coefficients of up to 0.61), and how biogeochemical reactions can vary in response to carbon and/or oxygen limitations. Lastly, we showed that incorporation of enzymatic regulations into substrate-explicit models is critical for more reasonable predictions. This result led us to present integrative biogeochemical modeling as a unifying framework that can ideally describe the dynamic interplay among microbes, enzymes, and substrates to address advanced questions and hypotheses in future studies. Altogether, the new modeling concept we propose in this work provides a foundational platform for unprecedented predictions of biogeochemical and ecosystem dynamics through enhanced integration with diverse experimental data and extant modeling approaches.

ABSTRACT The ecological drivers that concurrently act upon both a virus and its host and that drive community assembly are poorly understood despite known interactions between viral populations and their microbial hosts. Hydraulically fractured shale environments provide access to a closed ecosystem in the deep subsurface where constrained microbial and viral community assembly processes can be examined. Here, we used metagenomic analyses of time-resolved-produced fluid samples from two wells in the Appalachian Basin to track viral and host dynamics and to investigate community assembly processes. Hypersaline conditions within these ecosystems should drive microbial community structure to a similar configuration through time in response to common osmotic stress. However, viral predation appears to counterbalance this potentially strong homogeneous selection and pushes the microbial community toward undominated assembly. In comparison, while the viral community was also influenced by substantial undominated processes, it assembled, in part, due to homogeneous selection. When the overall assembly processes acting upon both these communities were directly compared with each other, a significant relationship was revealed, suggesting an association between microbial and viral community development despite differing selective pressures. These results reveal a potentially important balance of ecological dynamics that must be in maintained within this deep subsurface ecosystem in order for the microbial community to persist over extended time periods. More broadly, this relationship begins to provide knowledge underlying metacommunity development across trophic levels. IMPORTANCE Interactions between viral communities and their microbial hosts have been the subject of many recent studies in a wide range of ecosystems. The degree of coordination between ecological assembly processes influencing viral and microbial communities, however, has been explored to a much lesser degree. By using a combined null modeling approach, this study investigated the ecological assembly processes influencing both viral and microbial community structure within hydraulically fractured shale environments. Among other results, significant relationships between the structuring processes affecting both the viral and microbial community were observed, indicating that ecological assembly might be coordinated between these communities despite differing selective pressures. Within this deep subsurface ecosystem, these results reveal a potentially important balance of ecological dynamics that must be maintained to enable long-term microbial community persistence. More broadly, this relationship begins to provide insight into the development of communities across trophic levels.

Primate simplex viruses, including Herpes simplex viruses 1 and 2, form a group of closely related herpesviruses, which establish latent infections in neurons of their respective host species. While neuropathogenic infections in their natural hosts are rare, zoonotic transmission of Macacine alphaherpesvirus 1 (McHV1) from macaques to humans is associated with severe disease. Human infections with baboon-derived Papiine alphaherpesvirus 2 (PaHV2) have not been reported, although PaHV2 and McHV1 share several biological properties, including neuropathogenicity in mice. The reasons for potential differences in PaHV2 and McHV1 pathogenicity are presently not understood, and answering these questions will require mutagenic analysis. Here, we report the development of a recombinant system, which allows rescue of recombinant PaHV2. In addition, we used recombineering to generate viruses carrying reporter genes (Gaussia luciferase or enhanced green fluorescent protein), which replicate with similar efficiency as wild-type PaHV2. We demonstrate that these viruses can be used to analyze susceptibility of cells to infection and inhibition of infection by neutralizing antibodies and antiviral compounds. In summary, we created a recombinant system for PaHV2, which in the future will be invaluable for molecular analyses of neuropathogenicity of PaHV2.

River corridor metabolomes reflect organic matter (OM) processing that drives aquatic biogeochemical cycles. Recent work highlights the power of ultrahigh-resolution mass spectrometry for understanding metabolome composition and river corridor metabolism. However, there have been no studies on the global chemogeography of surface water and sediment metabolomes using ultrahigh-resolution techniques. Here, we describe a community science effort from the Worldwide Hydrobiogeochemistry Observation Network for Dynamic River Systems (WHONDRS) consortium to characterize global metabolomes in surface water and sediment that span multiple stream orders and biomes. We describe the distribution of key aspects of metabolomes including elemental groups, chemical classes, indices, and inferred biochemical transformations. We show that metabolomes significantly differ across surface water and sediment and that surface water metabolomes are more rich and variable. We also use inferred biochemical transformations to identify core metabolic processes shared among surface water and sediment. Finally, we observe significant spatial variation in sediment metabolites between rivers in the eastern and western portions of the contiguous United States. Our work not only provides a basis for understanding global patterns in river corridor biogeochemical cycles but also demonstrates that community science endeavors can enable global research projects that are unfeasible with traditional research models.

ABSTRACT The Worldwide Hydrobiogeochemistry Observation Network for Dynamic River Systems (WHONDRS) aims to galvanize a global community to provide the scientific basis for improved management of dynamic river corridors. WHONDRS is a global research consortium working to understand connections among dynamic hydrology, biogeochemistry, and microbiology in river corridors from local to global scales. WHONDRS ascribes to the perspective that resources, knowledge, and data belong to the community as a whole and that science advances more rapidly and more robustly through community ownership. As such, WHONDRS provides free access to novel instrumentation, molecular analysis, and well-curated data associated with river corridor hydrology, biogeochemistry, and microbiology. There are a number of ways to be involved in WHONDRS, ranging from one-time surface water sampling to installation of WHONDRS-developed multiparameter sensors for continuous monitoring. WHONDRS hinges on broad involvement, and we encourage all interested parties to contact us and become part of the consortium.

ABSTRACT The Worldwide Hydrobiogeochemistry Observation Network for Dynamic River Systems (WHONDRS) aims to galvanize a global community to provide the scientific basis for improved management of dynamic river corridors. WHONDRS is a global research consortium working to understand connections among dynamic hydrology, biogeochemistry, and microbiology in river corridors from local to global scales. WHONDRS ascribes to the perspective that resources, knowledge, and data belong to the community as a whole and that science advances more rapidly and more robustly through community ownership. As such, WHONDRS provides free access to novel instrumentation, molecular analysis, and well-curated data associated with river corridor hydrology, biogeochemistry, and microbiology. There are a number of ways to be involved in WHONDRS, ranging from one-time surface water sampling to installation of WHONDRS-developed multiparameter sensors for continuous monitoring. WHONDRS hinges on broad involvement, and we encourage all interested parties to contact us and become part of the consortium.

ABSTRACT To improve predictions of ecosystem function in future environments, we need to integrate the ecological and environmental histories experienced by microbial communities with hydrobiogeochemistry across scales. A key issue is whether we can derive generalizable scaling relationships that describe this multiscale integration. There is a strong foundation for addressing these challenges. We have the ability to infer ecological history with null models and reveal impacts of environmental history through laboratory and field experimentation. Recent developments also provide opportunities to inform ecosystem models with targeted omics data. A major next step is coupling knowledge derived from such studies with multiscale modeling frameworks that are predictive under non-steady-state conditions. This is particularly true for systems spanning dynamic interfaces, which are often hot spots of hydrobiogeochemical function. We can advance predictive capabilities through a holistic perspective focused on the nexus of history, ecology, and hydrobiogeochemistry.

In a recent study of denitrification dynamics in hyporheic zone sediments, we observed a significant time lag (up to several days) in enzymatic response to the changes in substrate concentration. To explore an underlying mechanism and understand the interactive dynamics between enzymes and nutrients, we developed a trait-based model that associates a community's traits with functional enzymes, instead of typically used species guilds (or functional guilds). This enzyme-based formulation allows to collectively describe biogeochemical functions of microbial communities without directly parameterizing the dynamics of species guilds, therefore being scalable to complex communities. As a key component of modeling, we accounted for microbial regulation occurring through transcriptional and translational processes, the dynamics of which was parameterized based on the temporal profiles of enzyme concentrations measured using a new signature peptide-based method. The simulation results using the resulting model showed several days of a time lag in enzymatic responses as observed in experiments. Further, the model showed that the delayed enzymatic reactions could be primarily controlled by transcriptional responses and that the dynamics of transcripts and enzymes are closely correlated. The developed model can serve as a useful tool for predicting biogeochemical processes in natural environments, either independently or through integration with hydrologic flow simulators.

Community assembly processes generate shifts in species abundances that influence ecosystem cycling of carbon and nutrients, yet our understanding of assembly remains largely separate from ecosystem-level functioning. Here, we investigate relationships between assembly and changes in microbial metabolism across space and time in hyporheic microbial communities. We pair sampling of two habitat types (i.e., attached and planktonic) through seasonal and sub-hourly hydrologic fluctuation with null modeling and temporally-explicit multivariate statistics. We demonstrate that multiple selective pressures—imposed by sediment and porewater physicochemistry—integrate to generate changes in microbial community composition at distinct timescales among habitat types. These changes in composition are reflective of contrasting associations of Betaproteobacteria and Thaumarchaeota with ecological selection and with seasonal changes in microbial metabolism. We present a conceptual model based on our results in which metabolism increases when oscillating selective pressures oppose temporally-stable selective pressures. Our conceptual model is pertinent to both macrobial and microbial systems experiencing multiple selective pressures and presents an avenue for assimilating community assembly processes into predictions of ecosystem-level functioning.

The number of microbial operational taxonomic units (OTUs) within a community is akin to species richness within plant/animal (‘macrobial’) systems. A large literature documents OTU richness patterns, drawing comparisons to macrobial theory. There is, however, an unrecognized fundamental disconnect between OTU richness and macrobial theory: OTU richness is commonly estimated on a per-individual basis, while macrobial richness is estimated per-area. Furthermore, the range or extent of sampled environmental conditions can strongly influence a study’s outcomes and conclusions, but this is not commonly addressed when studying OTU richness. Here we (i) propose a new sampling approach that estimates OTU richness per-mass of soil, which results in strong support for species energy theory, (ii) use data reduction to show how support for niche conservatism emerges when sampling across a restricted range of environmental conditions, and (iii) show how additional insights into drivers of OTU richness can be generated by combining different sampling methods while simultaneously considering patterns that emerge by restricting the range of environmental conditions. We propose that a more rigorous connection between microbial ecology and macrobial theory can be facilitated by exploring how changes in OTU richness units and environmental extent influence outcomes of data analysis. While fundamental differences between microbial and macrobial systems persist (e.g., species concepts), we suggest that closer attention to units and scale provide tangible and immediate improvements to our understanding of the processes governing OTU richness and how those processes relate to drivers of macrobial species richness.

Candidate phylum OD1 bacteria (also referred to as Parcubacteria) have been identified in broad range of anoxic environments through community survey analysis. Although none of these species have been isolated in the laboratory, several genome sequences have been reconstructed from metagenomic sequence data and single-cell sequencing. The organisms have small (generally

Ecological community assembly is governed by a combination of (i) selection resulting from among-taxa differences in performance; (ii) dispersal resulting from organismal movement; and (iii) ecological drift resulting from stochastic changes in population sizes. The relative importance and nature of these processes can vary across environments. Selection can be homogeneous or variable, and while dispersal is a rate, we conceptualize extreme dispersal rates as two categories; dispersal limitation results from limited exchange of organisms among communities, and homogenizing dispersal results from high levels of organism exchange. To estimate the influence and spatial variation of each process we extend a recently developed statistical framework, use a simulation model to evaluate the accuracy of the extended framework, and use the framework to examine subsurface microbial communities over two geologic formations. For each subsurface community we estimate the degree to which it is influenced by homogeneous selection, variable selection, dispersal limitation, and homogenizing dispersal. Our analyses revealed that the relative influences of these ecological processes vary substantially across communities even within a geologic formation. We further identify environmental and spatial features associated with each ecological process, which allowed mapping of spatial variation in ecological-process-influences. The resulting maps provide a new lens through which ecological systems can be understood; in the subsurface system investigated here they revealed that the influence of variable selection was associated with the rate at which redox conditions change with subsurface depth.

Phototrophic microbial mats are compact ecosystems composed of highly interactive organisms in which energy and element cycling take place over millimeter-to-centimeter-scale distances. Although microbial mats are common in hypersaline environments, they have not been extensively characterized in systems dominated by divalent ions. Hot Lake is a meromictic, epsomitic lake that occupies a small, endorheic basin in north-central Washington. The lake harbors a benthic, phototrophic mat that assembles each spring and disassembles each fall and is subject to greater than tenfold variation in salinity (primarily Mg2+ and SO42-) and irradiation over the annual cycle. We examined spatiotemporal variation in the mat community at five time points throughout the annual cycle with respect to prevailing physicochemical parameters by sequencing the V4 region of the 16S rRNA gene coupled to near-full-length 16S RNA clone sequences. The composition of these microbial communities was relatively stable over the seasonal cycle and included dominant populations of Cyanobacteria, primarily a group IV cyanobacterium (Leptolyngbya), and Alphaproteobacteria (specifically, members of Rhodobacteraceae and Geminicoccus). Members of Gammaproteobacteria (e.g., Thioalkalivibrio and Halochromatium) and Deltaproteobacteria (e.g., Desulfofustis) that are likely to be involved in sulfur cycling peaked in summer and declined significantly by mid-fall, mirroring larger trends in mat community richness and evenness. Phylogenetic turnover analysis of abundant phylotypes employing environmental metadata suggests that seasonal shifts in light variability exert a dominant influence on the composition of Hot Lake microbial mat communities. The seasonal development and organization of these structured microbial mats provide opportunities for analysis of the temporal and physical dynamics that feed back to community function.