Soil Bacterial Populations Are Shaped by Recombination and Gene-Specific
Total Page:16
File Type:pdf, Size:1020Kb
The ISME Journal (2020) 14:1834–1846 https://doi.org/10.1038/s41396-020-0655-x ARTICLE Soil bacterial populations are shaped by recombination and gene-specific selection across a grassland meadow 1 2 2 2 Alexander Crits-Christoph ● Matthew R. Olm ● Spencer Diamond ● Keith Bouma-Gregson ● Jillian F. Banfield 2,3,4,5 Received: 17 September 2019 / Revised: 22 March 2020 / Accepted: 2 April 2020 / Published online: 23 April 2020 © The Author(s) 2020. This article is published with open access Abstract 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 fi 1234567890();,: 1234567890();,: identity of assembled genomes was signi cantly 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. Introduction Soil microbial communities play key biogeochemical roles in terrestrial ecosystems [1]. Within a single hectare of temperate grassland soil, there can be over 1000 kg of microbial biomass and corresponding large microbial Supplementary information The online version of this article (https:// doi.org/10.1038/s41396-020-0655-x) contains supplementary population sizes [2]. Recent progress has been made in material, which is available to authorized users. cataloging the diversity of 16S rRNA genes in soils [3], which is useful for understanding microbial community * fi Jillian F. Ban eld composition, but this technique is incapable of discerning jbanfi[email protected] most genetic variation within populations [4]. In addition, 1 Department of Plant and Microbial Biology, University of many of the most common soil microorganisms, such as the California, Berkeley, CA, USA highly abundant Acidobacteria, Verrucomicrobia, and 2 Department of Earth and Planetary Science, University of Gemmatimonadetes phyla, are underrepresented in or nearly California, Berkeley, CA, USA absent from culture collections and genomic databases, 3 Department of Environmental Science, Policy, and Management, even at the level of class or phylum [1, 5, 6]. For these University of California, Berkeley, CA, USA reasons, the processes of recombination and selection in 4 Earth Sciences Division, Lawrence Berkeley National Laboratory, many of the most globally prolific soil microbial phyla Berkeley, CA, USA remain unstudied. However, genome-resolved metage- 5 Chan Zuckerberg Biohub, San Francisco, CA, USA nomics, in which shotgun sequencing of metagenomic Soil bacterial populations are shaped by recombination and gene-specific selection across a grassland. 1835 DNA is assembled and binned into draft genomes, has [27, 28]. In soil microbial populations, comparative fre- recently resulted in whole genome characterization of these quencies of gene-specific sweeps versus genome-wide clo- rarely cultivated but widespread soil bacteria [7–11]. nal strain competition and replacement are largely Sequence variation within genome-resolved metagenomic unknown, as are the spatial scales at which these evolu- datasets can therefore be used to track changes in allele tionary processes can occur. frequencies, and to infer the operation of evolutionary for- Previously, we conducted a large scale genome-resolved ces of genetic drift, natural selection, and homologous metagenomics study of soils from a grassland meadow in recombination in these natural populations [12, 13]. the Angelo Coast Range Reserve in northern California that Homologous recombination can vary dramatically in its established a dataset of 896 phylogenetically diverse importance relative to the other processes for different microbial genomes dereplicated by species, and reported on microbial species. Analyses of reference genomes have community composition [29]. The soil at the site is a sandy shown that homologous recombination frequently occurs in loam mixture of ~45% clay-, ~45% silt-, and 10% sand- bacteria populations, both globally and locally [14–16]. For sized particles [8], with pH values in the range of 4.6–4.9. example, certain populations of hotspring Cyanobacteria The mineralogy of the soils at the site was previously found approach panmixia, where recombination is so frequent that to be predominantly vermiculite, with plagioclase and alkali individual cells are unlinked random mixtures of alleles feldspars and minor apatite [30]. The soils were classified as [17]. In other species-like oceanic Vibrio, recombination is Ultic Haploxeralfs of the Holohan-Hollowtree-Casabonne high but large blocks of alleles important for ecological complex, and at a depth of 30 cm were found to have a bulk niche differentiation are co-inherited and may remain linked density of ~2.0 g/cm3, a cation exchange capacity between due to selection [18, 19]. In soils, Streptomyces flavogriseus 17 and 19 meq per 100 g soil, C:N ratios from 10–12, and isolates were also found to approach a freely recombining total C concentrations of 10 mg g−1 up to 18 mg g−1 [30]. panmixia [20]. In contrast, Myxoccocus xanthus isolates The grassland is dominated by annual Mediterranean recovered from a series of soil samples were distinct but grasses and forbs [8, 31]. The meadow has been part of a highly clonal, implying recombination between strains was rainfall amendment climate change study ongoing for over low [21]. However, despite their frequent cultivation, 17 years and has also been studied in the context of plant cultivation-independent studies show that those taxa are diversity and productivity [32], invertebrate herbivores and comparatively rare in soils compared with undercultivated predators [31], fungal communities [33], soil organic matter members of the Acidobacteria, Gemmatimonadetes, and [30], metabolomics, and metaproteomics [8]. By analyzing Verrucomicrobia [1]. The degree of recombination in these the population genomics of 19 highly abundant bacterial rarely cultivated but abundant soil bacterial lineages has not species across this meadow, we found high nucleotide been investigated, but recombination may be widespread, as diversity within samples, intrapopulation genetic structure high cell densities could promote the sharing of genetic often shifting over local spatial scales, varying degrees of material via transformation [22], conjugation [23], or the homologous recombination for different species, and gene- uptake of extracellular vesicles [24]. specific population differentiation partially driven by When recombination rates are low or selection is extre- selection. mely strong, several clonal strains compete until one or more beneficial alleles is highly selected for, resulting in a single clonal genotype increasing in abundance or even Materials and methods sweeping to fixation. However, when recombination unlinks gene variants within a population, beneficial alleles Sampling, genome sequencing, and metagenomic can sweep through a population independent of genomic assembly context in a selective sweep [25] with gene/locus-specific effects. One genome-resolved metagenomic study observed The sampling scheme, local soil characteristics, and study a single clonal sweep over a 9 year period for one Chlor- design that were utilized in this analysis have been pre- obium population in a freshwater lake [26], while most of viously described [29]. Previously, 60 soil samples were the other bacterial populations studied possessed genomic collected at depths of 10–20, 20–30, and 30–40 cm near the loci with unusually few SNPs (single nucleotide poly- respective centers of six 10 m diameter plots, spaced 5 m morphisms), an observation interpreted as evidence for apart. Sampling plots were spatially arranged in “blocks” of gene-specific selective sweeps. However, positive selection two plots across the meadow, and one of the two plots in acting on a genomic locus in a recombining population can each block received extended spring rainfall [31]. Samples also leave additional locus-specific signals, including higher were collected over a period of 2 months before and fol- linkage disequilibrium (a strong association between lowing autumn rainfall, resulting in 10 samples per plot alleles) and different allele frequencies between populations