Probing the ecological and evolutionary history of a thermophilic cyanobacterial population via statistical properties of its microdiversity

Genomic estimates of mutation and substitution rates contradict the evolutionary speed hypothesis of the latitudinal diversity gradient

The latitudinal diversity gradient (LDG) refers to a decrease in biodiversity from the equator to the poles. The evolutionary speed hypothesis, backed by the metabolic theory of ecology, asserts that nucleotide mutation and substitution rates per site per year are higher and thereby speciation rates are higher at higher temperatures, generating the LDG. However, prior empirical investigations of the relationship between the temperature and mutation or substitution rate were based on a few genes and the results were mixed. We here revisit this relationship using genomic data. No significant correlation between the temperature and mutation rate is found in 13 prokaryotes or in 107 eukaryotes. An analysis of 234 diverse trios of bacterial taxa indicates that the synonymous substitution rate is not significantly associated with the growth temperature. The same data, however, reveal a significant negative association between the nonsynonymous substitution rate and temperature, which is explainable by a larger fraction of detrimental nonsynonymous mutations at higher temperatures due to a stronger demand for protein stability. We conclude that the evolutionary speed hypothesis of the LDG is unsupported by genomic data and advise that future mechanistic studies of the LDG should focus on other hypotheses.

Hybridization breaks species barriers in long-term coevolution of a cyanobacterial population

Bacterial species often undergo rampant recombination yet maintain cohesive genomic identity. Ecological differences can generate recombination barriers between species and sustain genomic clusters in the short term. But can these forces prevent genomic mixing during long-term coevolution? Cyanobacteria in Yellowstone hot springs comprise several diverse species that have coevolved for hundreds of thousands of years, providing a rare natural experiment. By analyzing more than 300 single-cell genomes, we show that despite each species forming a distinct genomic cluster, much of the diversity within species is the result of hybridization driven by selection, which has mixed their ancestral genotypes. This widespread mixing is contrary to the prevailing view that ecological barriers can maintain cohesive bacterial species and highlights the importance of hybridization as a source of genomic diversity.

Multiple phase transitions shape biodiversity of a migrating population

In a wide variety of natural systems, closely related microbial strains coexist stably, resulting in high levels of fine-scale biodiversity. However, the mechanisms that stabilize this coexistence are not fully understood. Spatial heterogeneity is one common stabilizing mechanism, but the rate at which organisms disperse throughout the heterogeneous environment may strongly impact the stabilizing effect that heterogeneity can provide. An intriguing example is the gut microbiome, where active mechanisms affect the movement of microbes and potentially maintain diversity. We investigate how biodiversity is affected by migration rate using a simple evolutionary model with heterogeneous selection pressure. We find that the biodiversity-migration rate relationship is shaped by multiple phase transitions, including a reentrant phase transition to coexistence. At each transition, an ecotype goes extinct and dynamics exhibit critical slowing down (CSD). CSD is encoded in the statistics of fluctuations due to demographic noise-this may provide an experimental means for detecting and altering impending extinction.

Salt flat microbial diversity and dynamics across salinity gradient

Sabkhas are hypersaline, mineral-rich, supratidal mudflats that harbor microbes that are adapted to high salt concentration. Sabkha microbial diversity is generally studied for their community composition, but less is known about their genetic structure and heterogeneity. In this study, we analyzed a coastal sabkha for its microbial composition using 16S rDNA and whole metagenome, as well as for its population genetic structure. Our 16S rDNA analysis show high alpha diversity in both inner and edge sabkha than outer sabkha. Beta diversity result showed similar kind of microbial composition between inner and edge sabkha, while outer sabkha samples show different microbial composition. At phylum level, Bacteroidetes (~ 22 to 34%), Euryarchaeota (~ 18 to ~ 30%), unclassified bacteria (~ 24 to ~ 35%), Actinobacteria (~ 0.01 to ~ 11%) and Cyanobacteria (less than 1%) are predominantly found in both inside and edge sabkha regions, whereas Proteobacteria (~ 92 to ~ 97%) and Parcubacteria (~ 1 to ~ 2%) are predominately found in outer sabkha. Our 225 metagenomes assembly from this study showed similar bacterial community profile as observed in 16S rDNA-based analysis. From the assembled genomes, we found important genes that are involved in biogeochemical cycles and secondary metabolite biosynthesis. We observed a dynamic, thriving ecosystem that engages in metabolic activity that shapes biogeochemical structure via carbon fixation, nitrogen, and sulfur cycling. Our results show varying degrees of horizontal gene transfers (HGT) and homologous recombination, which correlates with the observed high diversity for these populations. Moreover, our pairwise population differentiation (Fst) for the abundance of species across the salinity gradient of sabkhas identified genes with strong allelic differentiation, lower diversity and elevated nonsynonymous to synonymous ratio of variants, which suggest selective sweeps for those gene variants. We conclude that the process of HGT, combined with recombination and gene specific selection, constitute the driver of genetic variation in bacterial population along a salinity gradient in the unique sabkha ecosystem.

Linkage disequilibrium between rare mutations

The statistical associations between mutations, collectively known as linkage disequilibrium, encode important information about the evolutionary forces acting within a population. Yet in contrast to single-site analogues like the site frequency spectrum, our theoretical understanding of linkage disequilibrium remains limited. In particular, little is currently known about how mutations with different ages and fitness costs contribute to expected patterns of linkage disequilibrium, even in simple settings where recombination and genetic drift are the major evolutionary forces. Here, I introduce a forward-time framework for predicting linkage disequilibrium between pairs of neutral and deleterious mutations as a function of their present-day frequencies. I show that the dynamics of linkage disequilibrium become much simpler in the limit that mutations are rare, where they admit a simple heuristic picture based on the trajectories of the underlying lineages. I use this approach to derive analytical expressions for a family of frequency-weighted linkage disequilibrium statistics as a function of the recombination rate, the frequency scale, and the additive and epistatic fitness costs of the mutations. I find that the frequency scale can have a dramatic impact on the shapes of the resulting linkage disequilibrium curves, reflecting the broad range of time scales over which these correlations arise. I also show that the differences between neutral and deleterious linkage disequilibrium are not purely driven by differences in their mutation frequencies and can instead display qualitative features that are reminiscent of epistasis. I conclude by discussing the implications of these results for recent linkage disequilibrium measurements in bacteria. This forward-time approach may provide a useful framework for predicting linkage disequilibrium across a range of evolutionary scenarios.

Microcoleus (Cyanobacteria) form watershed-wide populations without strong gradients in population structure

The relative importance of separation by distance and by environment to population genetic diversity can be conveniently tested in river networks, where these two drivers are often independently distributed over space. To evaluate the importance of dispersal and environmental conditions in shaping microbial population structures, we performed genome‐resolved metagenomic analyses of benthic Microcoleus‐dominated cyanobacterial mats collected in the Eel and Russian River networks (California, USA). The 64 Microcoleus genomes were clustered into three species that shared >96.5% average nucleotide identity (ANI). Most mats were dominated by one strain, but minor alleles within mats were often shared, even over large spatial distances (>300 km). Within the most common Microcoleus species, the ANI between the dominant strains within mats decreased with increasing spatial separation. However, over shorter spatial distances (tens of kilometres), mats from different subwatersheds had lower ANI than mats from the same subwatershed, suggesting that at shorter spatial distances environmental differences between subwatersheds in factors like canopy cover, conductivity, and mean annual temperature decreases ANI. Since mats in smaller creeks had similar levels of nucleotide diversity (π) as mats in larger downstream subwatersheds, within‐mat genetic diversity does not appear to depend on the downstream accumulation of upstream‐derived strains. The four‐gamete test and sequence length bias suggest recombination occurs between almost all strains within each species, even between populations separated by large distances or living in different habitats. Overall, our results show that, despite some isolation by distance and environmental conditions, sufficient gene‐flow occurs among cyanobacterial strains to prevent either driver from producing distinctive population structures across the watershed.

Soil bacterial populations are shaped by recombination and gene-specific selection across a grassland meadow

Soil microbial diversity is often studied from the perspective of community composition, but less is known about genetic heterogeneity within species. The relative impacts of clonal interference, gene-specific selection, and recombination in many abundant but rarely cultivated soil microbes remain unknown. Here we track genome-wide population genetic variation for 19 highly abundant bacterial species sampled from across a grassland meadow. Genomic inferences about population structure are made using the millions of sequencing reads that are assembled de novo into consensus genomes from metagenomes, as each read pair describes a short genomic sequence from a cell in each population. Genomic nucleotide identity of assembled genomes was significantly associated with local geography for over half of the populations studied, and for a majority of populations within-sample nucleotide diversity could often be as high as meadow-wide nucleotide diversity. Genes involved in metabolite biosynthesis and extracellular transport were characterized by elevated nucleotide diversity in multiple species. Microbial populations displayed varying degrees of homologous recombination and recombinant variants were often detected at 7-36% of loci genome-wide. Within multiple populations we identified genes with unusually high spatial differentiation of alleles, fewer recombinant events, elevated ratios of nonsynonymous to synonymous variants, and lower nucleotide diversity, suggesting recent selective sweeps for gene variants. Taken together, these results indicate that recombination and gene-specific selection commonly shape genetic variation in several understudied soil bacterial lineages.

Stabilization of extensive fine-scale diversity by ecologically driven spatiotemporal chaos

It has recently become apparent that the diversity of microbial life extends far below the species level to the finest scales of genetic differences. Remarkably, extensive fine-scale diversity can coexist spatially. How is this diversity stable on long timescales, despite selective or ecological differences and other evolutionary processes? Most work has focused on stable coexistence or assumed ecological neutrality. We present an alternative: extensive diversity maintained by ecologically driven spatiotemporal chaos, with no assumptions about niches or other specialist differences between strains. We study generalized Lotka-Volterra models with antisymmetric correlations in the interactions inspired by multiple pathogen strains infecting multiple host strains. Generally, these exhibit chaos with increasingly wild population fluctuations driving extinctions. But the simplest spatial structure, many identical islands with migration between them, stabilizes a diverse chaotic state. Some strains (subspecies) go globally extinct, but many persist for times exponentially long in the number of islands. All persistent strains have episodic local blooms to high abundance, crucial for their persistence as, for many, their average population growth rate is negative. Snapshots of the abundance distribution show a power law at intermediate abundances that is essentially indistinguishable from the neutral theory of ecology. But the dynamics of the large populations are much faster than birth-death fluctuations. We argue that this spatiotemporally chaotic "phase" should exist in a wide range of models, and that even in rapidly mixed systems, longer-lived spores could similarly stabilize a diverse chaotic phase.

Fine-Scale Haplotype Structure Reveals Strong Signatures of Positive Selection in a Recombining Bacterial Pathogen

Identifying genetic variation in bacteria that has been shaped by ecological differences remains an important challenge. For recombining bacteria, the sign and strength of linkage provide a unique lens into ongoing selection. We show that derived alleles <300 bp apart in Neisseria gonorrhoeae exhibit more coupling linkage than repulsion linkage, a pattern that cannot be explained by limited recombination or neutrality as these couplings are significantly stronger for nonsynonymous alleles than synonymous alleles. This general pattern is driven by a small fraction of highly diverse genes, many of which exhibit evidence of interspecies horizontal gene transfer and an excess of intermediate frequency alleles. Extensive simulations show that two distinct forms of positive selection can create these patterns of genetic variation: directional selection on horizontally transferred alleles or balancing selection that maintains distinct haplotypes in the presence of recombination. Our results establish a framework for identifying patterns of selection in fine-scale haplotype structure that indicate specific ecological processes in species that recombine with distantly related lineages or possess coexisting adaptive haplotypes.

Population Genetics in the Human Microbiome

While the human microbiome's structure and function have been extensively studied, its within-species genetic diversity is less well understood. However, genetic mutations in the microbiome can confer biomedically relevant traits, such as the ability to extract nutrients from food, metabolize drugs, evade antibiotics, and communicate with the host immune system. The population genetic processes by which these traits evolve are complex, in part due to interacting ecological and evolutionary forces in the microbiome. Advances in metagenomic sequencing, coupled with bioinformatics tools and population genetic models, facilitate quantification of microbiome genetic variation and inferences about how this diversity arises, evolves, and correlates with traits of both microbes and hosts. In this review, we explore the population genetic forces (mutation, recombination, drift, and selection) that shape microbiome genetic diversity within and between hosts, as well as efforts towards predictive models that leverage microbiome genetics.