Patterns of gene flow define species of thermophilic Archaea - PubMed (original) (raw)

Patterns of gene flow define species of thermophilic Archaea

Hinsby Cadillo-Quiroz et al. PLoS Biol. 2012 Feb.

Abstract

Despite a growing appreciation of their vast diversity in nature, mechanisms of speciation are poorly understood in Bacteria and Archaea. Here we use high-throughput genome sequencing to identify ongoing speciation in the thermoacidophilic Archaeon Sulfolobus islandicus. Patterns of homologous gene flow among genomes of 12 strains from a single hot spring in Kamchatka, Russia, demonstrate higher levels of gene flow within than between two persistent, coexisting groups, demonstrating that these microorganisms fit the biological species concept. Furthermore, rates of gene flow between two species are decreasing over time in a manner consistent with incipient speciation. Unlike other microorganisms investigated, we do not observe a relationship between genetic divergence and frequency of recombination along a chromosome, or other physical mechanisms that would reduce gene flow between lineages. Each species has its own genetic island encoding unique physiological functions and a unique growth phenotype that may be indicative of ecological specialization. Genetic differentiation between these coexisting groups occurs in large genomic "continents," indicating the topology of genomic divergence during speciation is not uniform and is not associated with a single locus under strong diversifying selection. These data support a model where species do not require physical barriers to gene flow but are maintained by ecological differentiation.

PubMed Disclaimer

Conflict of interest statement

The authors have declared that no competing interests exist.

Figures

Figure 1

Figure 1. Phylogenetic relationships among S. islandicus from the Mutnovsky Volcano.

(A) ClonalFrame reconstruction based on seven loci from 97 S. islandicus strains from the Mutnovsky Volcano region of Kamchatka, Russia (details in Table S1). Strains in purple were isolated from spring M.16. The first number in each name indicates the spring from which strains were isolated, the second indicates the isolate number from that spring, and the third indicates year of isolation. * designates strains selected for genome sequencing and comparison. (B) ClonalFrame phylogeny based on the core genome alignment of 12 S. islandicus strains from hot spring M.16.

Figure 2

Figure 2. Heat map representation of homologous recombination frequency for every donor/recipient pair of branches of the core genome phylogeny of 12 S. islandicus strains.

Recombination frequency is measured relative to its expectation under the prior of the ClonalOrigin model and color coded according to the upper left color/magnitude legend (light blue and blue for the frequency of recombination events below a 1∶1 ratio and yellow to red for the frequency of recombination events above 1∶1). Light gray cells represent non-significant ratios with less than four observed and expected events. White shows number of events that match the prior expectations. Names of strains are color coded as Blue and Red groups.

Figure 3

Figure 3. Variation in recombination events between the Red and Blue groups through time.

For the Red (A) and Blue (B) recipient strains, the total proportion of events that could be assigned as originating from either donor Red (colored red) and Blue (colored blue) strains is shown as a function of coalescent time with 10 being the most recent divergence and 1 being the common ancestor of this set of strains. A coalescent unit of time is equal to the average length of a generation multiplied by the effective population size.

Figure 4

Figure 4. FST values between the Red and Blue groups along the chromosome of strain M.16.27.

(A) 10,000 bp windows on the M.16.27 genome where genome sequence is not present in all 10 strains from the Red and Blue groups. These positions highlight variable portions of the M.16.27 genome. *Indicates a recently integrated plasmid. (B) FST values were calculated for sliding windows of 10 kb moving in 5 kb steps. Empty windows where sequence from M.16.27 is not shared by all strains are not plotted. Shading highlights regions of the chromosome that are less differentiated beginning and ending with the first window with FST values lower than 0.5.

Figure 5

Figure 5. Growth of M.16 strains under standard heterotrophic conditions.

Lines are color coded for strains assigned to the Red and Blue groups. Negative control with no inoculum added is shown in grey. Error bars show the variation in growth among three independent replicate cultures.

Similar articles

Cited by

References

    1. Acinas S. G, Klepac-Ceraj V, Hunt D. E, Pharino C, Ceraj I, et al. Fine-scale phylogenetic architecture of a complex bacterial community. Nature. 2004;430:551–554. - PubMed
    1. Casamayor E. O, Pedros-Alio C, Muyzer G, Amann R. Microheterogeneity in 16S ribosomal DNA-defined bacterial populations from a stratified planktonic environment is related to temporal changes and to ecological adaptations. Appl Environ Microbiol. 2002;68:1706–1714. - PMC - PubMed
    1. Hunt D. E, David L. A, Gevers D, Preheim S. P, Alm E. J, et al. Resource partitioning and sympatric differentiation among closely related bacterioplankton. Science. 2008;320(5879):1081–1085. - PubMed
    1. Sikorski J, Nevo E. Adaptation and incipient sympatric speciation of Bacillus simplex under microclimatic contrast at “Evolution Canyons” I and II, Israel. Proc Natl Acad Sci U S A. 2005;102:15924–15929. - PMC - PubMed
    1. Oakley B. B, Carbonero F, van der Gast C. J, Hawkins R. J, Purdy K. J. Evolutionary divergence and biogeography of sympatric niche-differentiated bacterial populations. ISME Journal. 2010;4:488–497. - PubMed

Publication types

MeSH terms

LinkOut - more resources