Microbiology of Big Soda Lake, a multi-extreme meromictic volcanic crater lake in the Nevada desert

Big Soda Lake, Nevada, is a multi-extreme meromictic lake, whose hypersaline hyperalkaline bottom waters feature permanent anoxia and high concentrations of arsenic, sulphide and ammonia. These properties make Big Soda Lake-and the adjacent Little Soda Lake-a fascinating system for exploring life's boundaries, discovering novel microbial taxa and identifying biotechnologically useful strains. To date, the taxonomic diversity and metabolic capabilities of microorganisms in this system remain largely unknown. Here, we fill this gap using microbiome surveys across the Big and Little Soda Lake water columns, including 16S rRNA sequencing, fungal ITS2 sequencing and gene- and genome-resolved metagenomics. We accompany these surveys with measurements of salinity, pH, temperature, oxygen, ammonium and ammonia concentrations. Our analyses reveal rich bacterial communities, taxonomically and functionally differentiated along Big Soda Lake's oxycline and, to lesser extent, between lakes. Fungal communities were dominated by a small number of families, while nearly no archaea were detected. Pathways related to perchlorate reduction, anoxygenic phototrophy, fermentation, dissimilatory metabolism of arsenite/arsenate, sulphur compounds, nitrogen compounds and hydrogen, were particularly prevalent. A total of 129 high-quality bacterial and archaeal metagenome-assembled genomes (completeness ≥ 80%, contamination ≤ 5%) were recovered, yielding insight into the taxonomic distribution of microbial metabolic pathways.

Trait biases in microbial reference genomes

Common culturing techniques and priorities bias our discovery towards specific traits that may not be representative of microbial diversity in nature. So far, these biases have not been systematically examined. To address this gap, here we use 116,884 publicly available metagenome-assembled genomes (MAGs, completeness ≥80%) from 203 surveys worldwide as a culture-independent sample of bacterial and archaeal diversity, and compare these MAGs to the popular RefSeq genome database, which heavily relies on cultures. We compare the distribution of 12,454 KEGG gene orthologs (used as trait proxies) in the MAGs and RefSeq genomes, while controlling for environment type (ocean, soil, lake, bioreactor, human, and other animals). Using statistical modeling, we then determine the conditional probabilities that a species is represented in RefSeq depending on its genetic repertoire. We find that the majority of examined genes are significantly biased for or against in RefSeq. Our systematic estimates of gene prevalences across bacteria and archaea in nature and gene-specific biases in reference genomes constitutes a resource for addressing these issues in the future.

Transport-limited reactions in microbial systems

Predicting microbial metabolic rates and emergent biogeochemical fluxes remains challenging due to the many unknown population dynamical, physiological and reaction-kinetic parameters and uncertainties in species composition. Here, we show that the need for these parameters can be eliminated when population dynamics and reaction kinetics operate at much shorter time scales than physical mixing processes. Such scenarios are widespread in poorly mixed water columns and sediments. In this 'fast-reaction-transport' (FRT) limit, all that is required for predictions are chemical boundary conditions, the physical mixing processes and reaction stoichiometries, while no knowledge of species composition, physiology or population/reaction kinetic parameters is needed. Using time-series data spanning years 2001-2014 and depths 180-900 m across the permanently anoxic Cariaco Basin, we demonstrate that the FRT approach can accurately predict the dynamics of major electron donors and acceptors (Pearson r ≥ 0.9 in all cases). Hence, many microbial processes in this system are largely transport limited and thus predictable regardless of species composition, population dynamics and kinetics. Our approach enables predictions for many systems in which microbial community dynamics and kinetics are unknown. Our findings also reveal a mechanism for the frequently observed decoupling between function and taxonomy in microbial systems.

The scaling of diversification rates with age is likely explained by sampling bias

Numerous phylogenetic studies reported the existence of a pervasive scaling relationship between the ages of extant eukaryotic clades and their estimated diversification rates. The causes of this age-rate-scaling (ARS), whether biological and/or artifactual, remain unresolved. Here we fit diversification models to thousands of eukaryotic time-calibrated phylogenies to explore multiple potential causes of the ARS including parameter non-identifiability, model inadequacy, biases in taxonomic practice, and an important and ubiquitous form of sampling bias-preferentially analyzing larger extant clades. We distinguish between two mechanism by which such sampling biases can cause an ARS: First, by favoring clades that happen to be unusually large merely by chance (i.e., due to the stochastic nature of the cladogenic process), thus leading to rate overestimation, and second, by favoring clades that have truly higher diversification rates. We find that, of the proposed explanations, only sampling biases are likely to contribute to the observed ARS. We develop methods for fully correcting for sampling bias mechanism 1, and find that despite these corrections a substantial ARS remains. We then confirm using simulations that preferring trees with truly higher rates (mechanism 2) likely explains this residual ARS. Since we do not have a completely unbiased sample of clades, including extinct ones, for phylogenetic analyses, it is difficult to demonstrate unambiguously that sampling biases are the sole cause of the ARS. Sampling biases are, however, a parsimonious and plausible explanation for this widely observed macroevolutionary pattern, and this has implications for how we interpret the distribution of diversification rate estimates in extant clades.

Solutions in microbiome engineering: prioritizing barriers to organism establishment

Microbiome engineering is increasingly being employed as a solution to challenges in health, agriculture, and climate. Often manipulation involves inoculation of new microbes designed to improve function into a preexisting microbial community. Despite, increased efforts in microbiome engineering inoculants frequently fail to establish and/or confer long-lasting modifications on ecosystem function. We posit that one underlying cause of these shortfalls is the failure to consider barriers to organism establishment. This is a key challenge and focus of macroecology research, specifically invasion biology and restoration ecology. We adopt a framework from invasion biology that summarizes establishment barriers in three categories: (1) propagule pressure, (2) environmental filtering, and (3) biotic interactions factors. We suggest that biotic interactions is the most neglected factor in microbiome engineering research, and we recommend a number of actions to accelerate engineering solutions.

The rates of global bacterial and archaeal dispersal

The phylogenetic resolution at which microorganisms display geographic endemism, the rates at which they disperse at global scales, and the role of humans on global microbial dispersal are largely unknown. Answering these questions is necessary for interpreting microbial biogeography, ecology, and macroevolution and for predicting the spread of emerging pathogenic strains. To resolve these questions, I analyzed the geographic and evolutionary relationships between 36,795 bacterial and archaeal ("prokaryotic") genomes from ∼7000 locations around the world. I find clear signs of continental-scale endemism, including strong correlations between phylogenetic divergence and geographic distance. However, the phylogenetic scale at which endemism generally occurs is extremely small, and most "species" (defined by an average nucleotide identity ≥ 95%) and even closely related strains (average nucleotide identity ≥ 99.9%) are globally distributed. Human-associated lineages display faster dispersal rates than other terrestrial lineages; the average net distance between any two human-associated cell lineages diverging 50 years ago is roughly 580 km. These results suggest that many previously reported global-scale microbial biogeographical patterns are likely the result of recent or current environmental filtering rather than geographic endemism. For human-associated lineages, estimated transition rates between Europe and North America are particularly high, and much higher than for non-human associated terrestrial lineages, highlighting the role that human movement plays in global microbial dispersal. Dispersal was slowest for hot spring- and terrestrial subsurface-associated lineages, indicating that these environments may act as "isolated islands" of microbial evolution.

SARS-CoV-2 infections in 165 countries over time

BACKGROUND

Understanding the dynamics of the COVID-19 pandemic and evaluating the efficacy of control measures requires knowledge of the number of infections over time. This number, however, often differs from the number of confirmed cases because of a large fraction of asymptomatic infections and different testing strategies.

METHODS

This study uses death count statistics, age-dependent infection fatality risks, and stochastic modeling to estimate the prevalence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections among adults (aged 20 years or older) in 165 countries over time, from early 2020 until June 25, 2021. The accuracy of the approach is confirmed through comparison with previous nationwide seroprevalence surveys.

RESULTS

The estimates presented reveal that the fraction of infections that are detected vary widely over time and between countries, and hence confirmed cases alone often yield a false picture of the pandemic. As of June 25, 2021, the nationwide cumulative fraction of SARS-CoV-2 infections (cumulative infections relative to population size) was estimated as 98% (95% confidence interval [CI] 93-100%) for Peru, 83% (95% CI 61-94%) for Brazil, and 36% (95% CI 23-61%) for the United States.

CONCLUSIONS

The time-resolved estimates presented expand the possibilities to study the factors that influenced and still influence the pandemic's progression in 165 countries.

Fundamental Identifiability Limits in Molecular Epidemiology

Viral phylogenies provide crucial information on the spread of infectious diseases, and many studies fit mathematical models to phylogenetic data to estimate epidemiological parameters such as the effective reproduction ratio (Re) over time. Such phylodynamic inferences often complement or even substitute for conventional surveillance data, particularly when sampling is poor or delayed. It remains generally unknown, however, how robust phylodynamic epidemiological inferences are, especially when there is uncertainty regarding pathogen prevalence and sampling intensity. Here, we use recently developed mathematical techniques to fully characterize the information that can possibly be extracted from serially collected viral phylogenetic data, in the context of the commonly used birth-death-sampling model. We show that for any candidate epidemiological scenario, there exists a myriad of alternative, markedly different, and yet plausible "congruent" scenarios that cannot be distinguished using phylogenetic data alone, no matter how large the data set. In the absence of strong constraints or rate priors across the entire study period, neither maximum-likelihood fitting nor Bayesian inference can reliably reconstruct the true epidemiological dynamics from phylogenetic data alone; rather, estimators can only converge to the "congruence class" of the true dynamics. We propose concrete and feasible strategies for making more robust epidemiological inferences from viral phylogenetic data.

Why extinction estimates from extant phylogenies are so often zero

Time-calibrated phylogenies of extant species ("extant timetrees") are widely used to estimate historical speciation and extinction rates by fitting stochastic birth-death models.¹ These approaches have long been controversial, as many phylogenetic studies report zero extinction in many taxa, contradicting the high extinction rates seen in the fossil record and the fact that the majority of species ever to have existed are now extinct.^2-9 To date, the causes of this discrepancy remain unresolved. Here, we provide a novel and simple explanation for these "zero-inflated" extinction estimates, based on the recent discovery that there exist many alternative "congruent" diversification scenarios that cannot be distinguished based solely on extant timetrees.¹⁰ Due to such congruencies, estimation methods tend to converge to some scenario congruent to (i.e., statistically indistinguishable from) the true diversification scenario, but not necessarily to the true diversification scenario itself. This congruent scenario may exhibit negative extinction rates, a biologically meaningless but mathematically feasible situation, in which case estimators will tend to stick to the boundary of zero extinction. Based on this explanation, we make multiple testable predictions, which we confirm using analyses of simulated trees and 121 empirical trees. In contrast to other proposed mechanisms for erroneous extinction rate estimates,^5,11-14 our proposed mechanism specifically explains the zero inflation of previous extinction rate estimates in the absence of detectable model violations, even for large trees. Not only do our results likely resolve a long-standing mystery in phylogenetics, they demonstrate that model congruencies can have severe consequences in practice.

Unifying Phylogenetic Birth-Death Models in Epidemiology and Macroevolution

Birth–death stochastic processes are the foundations of many phylogenetic models and are widely used to make inferences about epidemiological and macroevolutionary dynamics. There are a large number of birth–death model variants that have been developed; these impose different assumptions about the temporal dynamics of the parameters and about the sampling process. As each of these variants was individually derived, it has been difficult to understand the relationships between them as well as their precise biological and mathematical assumptions. Without a common mathematical foundation, deriving new models is nontrivial. Here, we unify these models into a single framework, prove that many previously developed epidemiological and macroevolutionary models are all special cases of a more general model, and illustrate the connections between these variants. This unification includes both models where the process is the same for all lineages and those in which it varies across types. We also outline a straightforward procedure for deriving likelihood functions for arbitrarily complex birth–death(-sampling) models that will hopefully allow researchers to explore a wider array of scenarios than was previously possible. By rederiving existing single-type birth–death sampling models, we clarify and synthesize the range of explicit and implicit assumptions made by these models. [Birth–death processes; epidemiology; macroevolution; phylogenetics; statistical inference.]

Effects of forced taxonomic transitions on metabolic composition and function in microbial microcosms

Surveys of microbial systems indicate that in many situations taxonomy and function may constitute largely independent ('decoupled') axes of variation. However, this decoupling is rarely explicitly tested experimentally, partly because it is hard to directly induce taxonomic variation without affecting functional composition. Here we experimentally evaluate this paradigm using microcosms resembling lake sediments and subjected to two different levels of salinity (0 and 19) and otherwise similar environmental conditions. We used DNA sequencing for taxonomic and functional profiling of bacteria and archaea and physicochemical measurements to monitor metabolic function, over 13 months. We found that the taxonomic composition of the saline systems gradually but strongly diverged from the fresh systems. In contrast, the metabolic composition (in terms of proportions of various genes) remained nearly identical across treatments and over time. Oxygen consumption rates and methane concentrations were substantially lower in the saline treatment, however, their similarity either increased (for oxygen) or did not change significantly (for methane) between the first and last sampling time, indicating that the lower metabolic activity in the saline treatments was directly and immediately caused by salinity rather than the gradual taxonomic divergence. Our experiment demonstrates that strong taxonomic shifts need not directly affect metabolic rates.

Phylogeographic Estimation and Simulation of Global Diffusive Dispersal

The analysis of time-resolved phylogenies (timetrees) and geographic location data allows estimation of dispersal rates, for example, for invasive species and infectious diseases. Many estimation methods are based on the Brownian Motion model for diffusive dispersal on a 2D plane; however, the accuracy of these methods deteriorates substantially when dispersal occurs at global scales because spherical Brownian motion (SBM) differs from planar Brownian motion. No statistical method exists for estimating SBM diffusion coefficients from a given timetree and tip coordinates, and no method exists for simulating SBM along a given timetree. Here, I present new methods for simulating SBM along a given timetree, and for estimating SBM diffusivity from a given timetree and tip coordinates using a modification of Felsenstein's independent contrasts and maximum likelihood. My simulation and fitting methods can accommodate arbitrary time-dependent diffusivities and scale efficiently to trees with millions of tips, thus enabling new analyses even in cases where planar BM would be a sufficient approximation. I demonstrate these methods using a timetree of marine and terrestrial Cyanobacterial genomes, as well as timetrees of two globally circulating Influenza B clades. My methods are implemented in the R package "castor." [Independent contrasts; phylogenetic; random walk; simulation; spherical Brownian motion.].

Simulating trees with millions of species

Motivation

The birth-death (BD) model constitutes the theoretical backbone of most phylogenetic tools for reconstructing speciation/extinction dynamics over time. Performing simulations of reconstructed trees (linking extant taxa) under the BD model in backward time, conditioned on the number of species sampled at present day and, in some cases, a specific time interval since the most recent common ancestor (MRCA), is needed for assessing the performance of reconstruction tools, for parametric bootstrapping and for detecting data outliers. The few simulation tools that exist scale poorly to large modern phylogenies, which can comprise thousands or even millions of tips (and rising).

Results

Here I present efficient software for simulating reconstructed phylogenies under time-dependent BD models in backward time, conditioned on the number of sampled species and (optionally) on the time since the MRCA. On large trees, my software is 1000–10 000 times faster than existing tools.

Availability and implementation

The presented software is incorporated into the R package ‘castor’, which is available on The Comprehensive R Archive Network (CRAN).

Supplementary information

Supplementary data are available at Bioinformatics online.

Microbial metabolite fluxes in a model marine anoxic ecosystem

Permanently anoxic regions in the ocean are widespread and exhibit unique microbial metabolic activity exerting substantial influence on global elemental cycles and climate. Reconstructing microbial metabolic activity rates in these regions has been challenging, due to the technical difficulty of direct rate measurements. In Cariaco Basin, which is the largest permanently anoxic marine basin and an important model system for geobiology, long-term monitoring has yielded time series for the concentrations of biologically important compounds; however, the underlying metabolite fluxes remain poorly quantified. Here, we present a computational approach for reconstructing vertical fluxes and in situ net production/consumption rates from chemical concentration data, based on a 1-dimensional time-dependent diffusive transport model that includes adaptive penalization of overfitting. We use this approach to estimate spatiotemporally resolved fluxes of oxygen, nitrate, hydrogen sulfide, ammonium, methane, and phosphate within the sub-euphotic Cariaco Basin water column (depths 150-900 m, years 2001-2014) and to identify hotspots of microbial chemolithotrophic activity. Predictions of the fitted models are in excellent agreement with the data and substantially expand our knowledge of the geobiology in Cariaco Basin. In particular, we find that the diffusivity, and consequently fluxes of major reductants such as hydrogen sulfide, and methane, is about two orders of magnitude greater than previously estimated, thus resolving a long-standing apparent conundrum between electron donor fluxes and measured dark carbon assimilation rates.

A General and Efficient Algorithm for the Likelihood of Diversification and Discrete-Trait Evolutionary Models

As the size of phylogenetic trees and comparative data continue to grow and more complex models are developed to investigate the processes that gave rise to them, macroevolutionary analyses are becoming increasingly limited by computational requirements. Here, we introduce a novel algorithm, based on the “flow” of the differential equations that describe likelihoods along tree edges in backward time, to reduce redundancy in calculations and efficiently compute the likelihood of various macroevolutionary models. Our algorithm applies to several diversification models, including birth–death models and models that account for state- or time-dependent rates, as well as many commonly used models of discrete-trait evolution, and provides an alternative way to describe macroevolutionary model likelihoods. As a demonstration of our algorithm’s utility, we implemented it for a popular class of state-dependent diversification models—BiSSE, MuSSE, and their extensions to hidden-states. Our implementation is available through the R package $\texttt{castor}$. We show that, for these models, our algorithm is one or more orders of magnitude faster than existing implementations when applied to large phylogenies. Our algorithm thus enables the fitting of state-dependent diversification models to modern massive phylogenies with millions of tips and may lead to potentially similar computational improvements for many other macroevolutionary models.

Efficient comparative phylogenetics on large trees

Motivation

Biodiversity databases now comprise hundreds of thousands of sequences and trait records. For example, the Open Tree of Life includes over 1 491 000 metazoan and over 300 000 bacterial taxa. These data provide unique opportunities for analysis of phylogenetic trait distribution and reconstruction of ancestral biodiversity. However, existing tools for comparative phylogenetics scale poorly to such large trees, to the point of being almost unusable.

Results

Here we present a new R package, named 'castor', for comparative phylogenetics on large trees comprising millions of tips. On large trees castor is often 100-1000 times faster than existing tools.

Availability and implementation

The castor source code, compiled binaries, documentation and usage examples are freely available at the Comprehensive R Archive Network (CRAN).

Contact

louca.research@gmail.com.

Supplementary information

Supplementary data are available at Bioinformatics online.

Circumventing kinetics in biogeochemical modeling

Microbial metabolism drives biogeochemical fluxes in virtually every ecosystem. Modeling these fluxes is challenged by the incredible diversity of microorganisms, whose kinetic parameters are largely unknown. In poorly mixed systems, such as stagnant water columns or sediments, however, long-term bulk microbial metabolism may become limited by physical transport rates of substrates across space. Here we mathematically show that under these conditions, biogeochemical fluxes are largely predictable based on the system's transport properties, chemical boundary conditions, and the stoichiometry of metabolic pathways, regardless of the precise kinetics of the resident microorganisms. We formalize these considerations into a predictive modeling framework and demonstrate its use for the Cariaco Basin subeuphotic zone, one of the largest anoxic marine basins worldwide. Using chemical concentration data solely from the upper boundary (depth 180 m) and lower boundary (depth 900 m), but without a priori knowledge of metabolite fluxes, chemical depth profiles, kinetic parameters, or microbial species composition, we predict the concentrations and vertical fluxes of biologically important substances, including oxygen, nitrate, hydrogen sulfide, and ammonium, across the entire considered depth range (180-900 m). Our predictions largely agree with concentration measurements over a period of 14 years ([Formula: see text] = 0.78-0.92) and become particularly accurate during a period where the system was near biogeochemical steady state (years 2007-2009, [Formula: see text] = 0.86-0.95). Our work enables geobiological predictions for a large class of ecosystems without knowledge of kinetic parameters or geochemical depth profiles. Conceptually, our work provides a possible explanation for the decoupling between microbial species composition and bulk metabolic function, observed in various ecosystems.

A census-based estimate of Earth's bacterial and archaeal diversity

The global diversity of Bacteria and Archaea, the most ancient and most widespread forms of life on Earth, is a subject of intense controversy. This controversy stems largely from the fact that existing estimates are entirely based on theoretical models or extrapolations from small and biased data sets. Here, in an attempt to census the bulk of Earth's bacterial and archaeal ("prokaryotic") clades and to estimate their overall global richness, we analyzed over 1.7 billion 16S ribosomal RNA amplicon sequences in the V4 hypervariable region obtained from 492 studies worldwide, covering a multitude of environments and using multiple alternative primers. From this data set, we recovered 739,880 prokaryotic operational taxonomic units (OTUs, 16S-V4 gene clusters at 97% similarity), a commonly used measure of microbial richness. Using several statistical approaches, we estimate that there exist globally about 0.8–1.6 million prokaryotic OTUs, of which we recovered somewhere between 47%–96%, representing >99.98% of prokaryotic cells. Consistent with this conclusion, our data set independently "recaptured" 91%–93% of 16S sequences from multiple previous global surveys, including PCR-independent metagenomic surveys. The distribution of relative OTU abundances is consistent with a log-normal model commonly observed in larger organisms; the total number of OTUs predicted by this model is also consistent with our global richness estimates. By combining our estimates with the ratio of full-length versus partial-length (V4) sequence diversity in the SILVA sequence database, we further estimate that there exist about 2.2–4.3 million full-length OTUs worldwide. When restricting our analysis to the Americas, while controlling for the number of studies, we obtain similar richness estimates as for the global data set, suggesting that most OTUs are globally distributed. Qualitatively similar results are also obtained for other 16S similarity thresholds (90%, 95%, and 99%). Our estimates constrain the extent of a poorly quantified rare microbial biosphere and refute recent predictions that there exist trillions of prokaryotic OTUs.

A massive survey of Earth's Bacteria and Archaea reveals that their diversity is orders of magnitude lower than previously thought. The study also indicates that extinctions played an important role in prokaryotic evolution.

The global diversity of Bacteria and Archaea ("prokaryotes"), the most ancient and most widespread forms of life on Earth, is subject to high uncertainty. Here, to estimate the global diversity of prokaryotes, we analyzed a large number of 16S ribosomal RNA gene sequences, found in all prokaryotes and commonly used to catalogue prokaryotic diversity. Sequences were obtained from a multitude of environments across thousands of geographic locations worldwide. From this data set, we recovered 739,880 prokaryotic operational taxonomic units (OTUs), i.e., 16S gene clusters sharing 97% similarity, roughly corresponding to prokaryotic species. Using several statistical approaches and through comparison with existing databases and previous independent surveys, we estimate that there exist globally between 0.8 and 1.6 million prokaryotic OTUs. When restricting our analysis to the Americas, while controlling for the number of studies, we obtain similar estimates as for the global data set, suggesting that most OTUs are not restricted to a single continent but are instead globally distributed. Our estimates constrain the extent of a commonly hypothesized but poorly quantified rare prokaryotic biosphere and refute recent predictions that there exists trillions of prokaryotic OTUs. Our findings also indicate that, contrary to common speculation, extinctions may strongly influence global prokaryotic diversity.

Correcting for 16S rRNA gene copy numbers in microbiome surveys remains an unsolved problem

The 16S ribosomal RNA gene is the most widely used marker gene in microbial ecology. Counts of 16S sequence variants, often in PCR amplicons, are used to estimate proportions of bacterial and archaeal taxa in microbial communities. Because different organisms contain different 16S gene copy numbers (GCNs), sequence variant counts are biased towards clades with greater GCNs. Several tools have recently been developed for predicting GCNs using phylogenetic methods and based on sequenced genomes, in order to correct for these biases. However, the accuracy of those predictions has not been independently assessed. Here, we systematically evaluate the predictability of 16S GCNs across bacterial and archaeal clades, based on ∼ 6,800 public sequenced genomes and using several phylogenetic methods. Further, we assess the accuracy of GCNs predicted by three recently published tools (PICRUSt, CopyRighter, and PAPRICA) over a wide range of taxa and for 635 microbial communities from varied environments. We find that regardless of the phylogenetic method tested, 16S GCNs could only be accurately predicted for a limited fraction of taxa, namely taxa with closely to moderately related representatives (≲15% divergence in the 16S rRNA gene). Consistent with this observation, we find that all considered tools exhibit low predictive accuracy when evaluated against completely sequenced genomes, in some cases explaining less than 10% of the variance. Substantial disagreement was also observed between tools (R2<0.5) for the majority of tested microbial communities. The nearest sequenced taxon index (NSTI) of microbial communities, i.e., the average distance to a sequenced genome, was a strong predictor for the agreement between GCN prediction tools on non-animal-associated samples, but only a moderate predictor for animal-associated samples. We recommend against correcting for 16S GCNs in microbiome surveys by default, unless OTUs are sufficiently closely related to sequenced genomes or unless a need for true OTU proportions warrants the additional noise introduced, so that community profiles remain interpretable and comparable between studies.

Functional structure of the bromeliad tank microbiome is strongly shaped by local geochemical conditions

Phytotelmata in tank-forming Bromeliaceae plants are regarded as potential miniature models for aquatic ecology, but detailed investigations of their microbial communities are rare. Hence, the biogeochemistry in bromeliad tanks remains poorly understood. Here we investigate the structure of bacterial and archaeal communities inhabiting the detritus within the tanks of two bromeliad species, Aechmea nudicaulis and Neoregelia cruenta, from a Brazilian sand dune forest. We used metagenomic sequencing for functional community profiling and 16S sequencing for taxonomic profiling. We estimated the correlation between functional groups and various environmental variables, and compared communities between bromeliad species. In all bromeliads, microbial communities spanned a metabolic network adapted to oxygen-limited conditions, including all denitrification steps, ammonification, sulfate respiration, methanogenesis, reductive acetogenesis and anoxygenic phototrophy. Overall, CO2 reducers dominated in abundance over sulfate reducers, and anoxygenic phototrophs largely outnumbered oxygenic photoautotrophs. Functional community structure correlated strongly with environmental variables, between and within a single bromeliad species. Methanogens and reductive acetogens correlated with detrital volume and canopy coverage, and exhibited higher relative abundances in N. cruenta. A comparison of bromeliads to freshwater lake sediments and soil from around the world, revealed stark differences in terms of taxonomic as well as functional microbial community structure.

Taxonomic variability and functional stability in microbial communities infected by phages

Microbial communities can display large variation in taxonomic composition, yet this variation can coincide with stable metabolic functional structure and performance. The mechanisms driving the taxonomic variation within functional groups remain largely unknown. Biotic interactions, such as predation by phages, have been suggested as potential cause of taxonomic turnover, but the conditions for this scenario have not been rigorously examined. Further, it is unknown how predation by phages affects community function, and how these effects are modulated by functional redundancy in the communities. Here, we address these questions using a model for a methanogenic microbial community that includes several interacting metabolic functional groups. Each functional group comprises multiple competing clades, and each clade is attacked by a specialist lytic phage. Our model predicts that phages induce intense taxonomic turnover, resembling the variability observed in previous experiments. The functional structure and performance of the community are also disturbed by phage predation, but they become more stable as the functional redundancy in the community increases. The extent of this stabilization depends on the particular functions considered. Our model suggests mechanisms by which functional redundancy stabilizes community function and supports the interpretation that biotic interactions promote taxonomic turnover within microbial functional groups.

High taxonomic variability despite stable functional structure across microbial communities

Understanding the processes that are driving variation of natural microbial communities across space or time is a major challenge for ecologists. Environmental conditions strongly shape the metabolic function of microbial communities; however, other processes such as biotic interactions, random demographic drift or dispersal limitation may also influence community dynamics. The relative importance of these processes and their effects on community function remain largely unknown. To address this uncertainty, here we examined bacterial and archaeal communities in replicate 'miniature' aquatic ecosystems contained within the foliage of wild bromeliads. We used marker gene sequencing to infer the taxonomic composition within nine metabolic functional groups, and shotgun environmental DNA sequencing to estimate the relative abundances of these groups. We found that all of the bromeliads exhibited remarkably similar functional community structures, but that the taxonomic composition within individual functional groups was highly variable. Furthermore, using statistical analyses, we found that non-neutral processes, including environmental filtering and potentially biotic interactions, at least partly shaped the composition within functional groups and were more important than spatial dispersal limitation and demographic drift. Hence both the functional structure and taxonomic composition within functional groups of natural microbial communities may be shaped by non-neutral and roughly separate processes.

Integrating biogeochemistry with multiomic sequence information in a model oxygen minimum zone

Microorganisms are the most abundant lifeform on Earth, mediating global fluxes of matter and energy. Over the past decade, high-throughput molecular techniques generating multiomic sequence information (DNA, mRNA, and protein) have transformed our perception of this microcosmos, conceptually linking microorganisms at the individual, population, and community levels to a wide range of ecosystem functions and services. Here, we develop a biogeochemical model that describes metabolic coupling along the redox gradient in Saanich Inlet-a seasonally anoxic fjord with biogeochemistry analogous to oxygen minimum zones (OMZs). The model reproduces measured biogeochemical process rates as well as DNA, mRNA, and protein concentration profiles across the redox gradient. Simulations make predictions about the role of ubiquitous OMZ microorganisms in mediating carbon, nitrogen, and sulfur cycling. For example, nitrite "leakage" during incomplete sulfide-driven denitrification by SUP05 Gammaproteobacteria is predicted to support inorganic carbon fixation and intense nitrogen loss via anaerobic ammonium oxidation. This coupling creates a metabolic niche for nitrous oxide reduction that completes denitrification by currently unidentified community members. These results quantitatively improve previous conceptual models describing microbial metabolic networks in OMZs. Beyond OMZ-specific predictions, model results indicate that geochemical fluxes are robust indicators of microbial community structure and reciprocally, that gene abundances and geochemical conditions largely determine gene expression patterns. The integration of real observational data, including geochemical profiles and process rate measurements as well as metagenomic, metatranscriptomic and metaproteomic sequence data, into a biogeochemical model, as shown here, enables holistic insight into the microbial metabolic network driving nutrient and energy flow at ecosystem scales.

Calibration and analysis of genome-based models for microbial ecology

Microbial ecosystem modeling is complicated by the large number of unknown parameters and the lack of appropriate calibration tools. Here we present a novel computational framework for modeling microbial ecosystems, which combines genome-based model construction with statistical analysis and calibration to experimental data. Using this framework, we examined the dynamics of a community of Escherichia coli strains that emerged in laboratory evolution experiments, during which an ancestral strain diversified into two coexisting ecotypes. We constructed a microbial community model comprising the ancestral and the evolved strains, which we calibrated using separate monoculture experiments. Simulations reproduced the successional dynamics in the evolution experiments, and pathway activation patterns observed in microarray transcript profiles. Our approach yielded detailed insights into the metabolic processes that drove bacterial diversification, involving acetate cross-feeding and competition for organic carbon and oxygen. Our framework provides a missing link towards a data-driven mechanistic microbial ecology.

Assessing host extinction risk following exposure to Batrachochytrium dendrobatidis

Wildlife diseases are increasingly recognized as a major threat to biodiversity. Chytridiomycosis is an emerging infectious disease of amphibians caused by the fungus Batrachochytrium dendrobatidis (Bd). Using a mathematical model and simulations, we study its effects on a generic riparian host population with a tadpole and adult life stage. An analytical expression for the basic reproduction quotient, Qo, of the pathogen is derived. By sampling the entire relevant parameter space, we perform a statistical assessment of the importance of all considered parameters in determining the risk of host extinction, upon exposure to Bd. We find that Qo not only gives a condition for the initial invasion of the fungus, but is in fact the best predictor for host extinction. We also show that the role of tadpoles, which in some species tolerate infections, is ambivalent. While tolerant tadpoles may provide a reservoir for the fungus, thus facilitating its persistence or even amplifying its outbreaks, they can also act as a rescue buffer for a stressed host population. Our results have important implications for amphibian conservation efforts.

Salivary cortisol and psychological factors in women with chronic and acute oro-facial pain

The aim of this study was to compare the salivary cortisol level, pain intensity and psychological factors between patients with chronic and acute oro-facial pain (OP) and pain-free subjects. Twenty-seven females with chronic OP (a diagnosis of myofascial pain according to the Research Diagnostic Criteria for Temporomandibular Disorders with at least 6 months duration), 24 females with acute OP (<10 days duration) and 27 pain-free females participated. Morning saliva was collected from all participants for analyses of the cortisol level. The pain intensity was assessed on a 0-10 numeric rating scale. The participants were evaluated by the Symptom Checklist 90-revised for levels of depression and somatisation, and the Perceived Stress Scale. The cortisol levels among the three patient groups were similar with no significant group differences. The median (interquartile range) current pain level did not differ between chronic and acute OP and was, respectively, 5 (4) and 5 (3). Patients with chronic OP showed significantly higher scores for depression, somatisation and perceived stress compared with patients with acute OP (Ps < 0.001), but there were no significant differences between acute OP and controls. To conclude, there were no differences in cortisol level between groups, despite significant higher levels of depression, somatisation and perceived stress in patients with chronic OP. This shows that psychological distress has a more important role in chronic than in acute OP. However, the relation between pain, adreno-cortical activity and psychological distress is complex and warrants further investigation.

Specialized nursery pollination mutualisms as evolutionary traps stabilized by antagonistic traits

We examine the conditions for the transition from antagonism to mutualism between plants and their specialists nursery pollinators in a reference case which is the Trollius europaeus-Chiastocheta interaction. The mechanistic model we developed shows that a specialization of T. europaeus on Chiastocheta could be the result of an attempt to escape over-exploitation by closing its flower. The pressure for such an escape increases with the parasite's frequency and its pollination efficiency but decreases in the presence of alternative pollinators. The resulting specialization is a priori an unstable one, leading either to strong evolutionary oscillations, or to evolutionary suicide due to over-exploitation of the plants. It becomes stable if the plants develop a defense mechanism to regulate their parasite's population size and limit seed-exploitation. The development of a counter-measure by the latter can destabilize the mutualism depending on the costs linked to such a trait. On the other hand, we find that a specialization on a purely mutualistic basis would require a preexisting high diversity of flower-opening within the population.