Searching for species in haloarchaea

The Origin of Niches and Species in the Bacterial World

Niches are spaces for the biological units of selection, from cells to complex communities. In a broad sense, "species" are biological units of individuation. Niches do not exist without individual organisms, and every organism has a niche. We use "niche" in the Hutchinsonian sense as an abstraction of a multidimensional environmental space characterized by a variety of conditions, both biotic and abiotic, whose quantitative ranges determine the positive or negative growth rates of the microbial individual, typically a species, but also parts of the communities of species contained in this space. Microbial organisms ("species") constantly diversify, and such diversification (radiation) depends on the possibility of opening up unexploited or insufficiently exploited niches. Niche exploitation frequently implies "niche construction," as the colonized niche evolves with time, giving rise to new potential subniches, thereby influencing the selection of a series of new variants in the progeny. The evolution of niches and organisms is the result of reciprocal interacting processes that form a single unified process. Centrifugal microbial diversification expands the limits of the species' niches while a centripetal or cohesive process occurs simultaneously, mediated by horizontal gene transfers and recombinatorial events, condensing all of the information recovered during the diversifying specialization into "novel organisms" (possible future species), thereby creating a more complex niche, where the selfishness of the new organism(s) establishes a "homeostatic power" limiting the niche's variation. Once the niche's full carrying capacity has been reached, reproductive isolation occurs, as no foreign organisms can outcompete the established population/community, thereby facilitating speciation. In the case of individualization-speciation of the microbiota, its contribution to the animal' gut structure is a type of "niche construction," the result of crosstalk between the niche (host) and microorganism(s). Lastly, there is a parallelism between the hierarchy of niches and that of microbial individuals. The increasing anthropogenic effects on the biosphere (such as globalization) might reduce the diversity of niches and bacterial individuals, with the potential emergence of highly transmissible multispecialists (which are eventually deleterious) resulting from the homogenization of the microbiosphere, a possibility that should be explored and prevented.

Insights into gene expression changes under conditions that facilitate horizontal gene transfer (mating) of a model archaeon

Horizontal gene transfer is a means by which bacteria, archaea, and eukaryotes are able to trade DNA within and between species. While there are a variety of mechanisms through which this genetic exchange can take place, one means prevalent in the archaeon Haloferax volcanii involves the transient formation of cytoplasmic bridges between cells and is referred to as mating. This process can result in the exchange of very large fragments of DNA between the participating cells. Genes governing the process of mating, including triggers to initiate mating, mechanisms of cell fusion, and DNA exchange, have yet to be characterized. We used a transcriptomic approach to gain a more detailed knowledge of how mating might transpire. By examining the differential expression of genes expressed in cells harvested from mating conditions on a filter over time and comparing them to those expressed in a shaking culture, we were able to identify genes and pathways potentially associated with mating. These analyses provide new insights into both the mechanisms and barriers of mating in Hfx. volcanii.

Broad Environmental Tolerance for a Salicola Host-Phage Pair Isolated from the Cargill Solar Saltworks, Newark, CA, USA

Phages greatly influence the ecology and evolution of their bacterial hosts; however, compared to hosts, a relatively low number of phages, especially halophilic phages, have been studied. This study describes a comparative investigation of physicochemical tolerance between a strain of the halophilic bacterium, Salicola, isolated from the Cargill Saltworks (Newark, CA, USA) and its associated phage. The host grew in media between pH 6–8.5, had a salinity growth optimum of 20% total salts (ranging from 10%–30%) and an upper temperature growth limit of 48 °C. The host utilized 61 of 190 substrates tested using BIOLOG Phenotype MicroArrays. The CGφ29 phage, one of only four reported Salicola phages, is a DNA virus of the Siphoviridae family. Overall, the phage tolerated a broader range of environmental conditions than its host (salinity 0–30% total salts; pH 3–9; upper thermal limit 80 °C) and is the most thermotolerant halophilic phage ever reported. This study is the most comprehensive investigation to date of a Salicola host–phage pair and provides novel insights into extreme environmental tolerances among bacteriophages.

Within-Species Genomic Variation and Variable Patterns of Recombination in the Tetracycline Producer Streptomyces rimosus

Streptomyces rimosus is best known as the primary source of the tetracycline class of antibiotics, most notably oxytetracycline, which have been widely used against many gram-positive and gram-negative pathogens and protozoan parasites. However, despite the medical and agricultural importance of S. rimosus, little is known of its evolutionary history and genome dynamics. In this study, we aim to elucidate the pan-genome characteristics and phylogenetic relationships of 32 S. rimosus genomes. The S. rimosus pan-genome contains more than 22,000 orthologous gene clusters, and approximately 8.8% of these genes constitutes the core genome. A large part of the accessory genome is composed of 9,646 strain-specific genes. S. rimosus exhibits an open pan-genome (decay parameter α = 0.83) and high gene diversity between strains (genomic fluidity φ = 0.12). We also observed strain-level variation in the distribution and abundance of biosynthetic gene clusters (BGCs) and that each individual S. rimosus genome has a unique repertoire of BGCs. Lastly, we observed variation in recombination, with some strains donating or receiving DNA more often than others, strains that tend to frequently recombine with specific partners, genes that often experience recombination more than others, and variable sizes of recombined DNA sequences. We conclude that the high levels of inter-strain genomic variation in S. rimosus is partly explained by differences in recombination among strains. These results have important implications on current efforts for natural drug discovery, the ecological role of strain-level variation in microbial populations, and addressing the fundamental question of why microbes have pan-genomes.

Impact of Homologous Recombination on the Evolution of Prokaryotic Core Genomes

Microbial populations exchange genetic material through a process called homologous recombination. Although this process has been studied in particular organisms, we lack an understanding of its differential impact over the genome and across microbes with different life-styles. We used a common analytical framework to assess this process in a representative set of microorganisms. Our results uncovered important trends. First, microbes with different lifestyles are differentially impacted, with endosymbionts and obligate pathogens being those less prone to undergo this process. Second, certain genetic elements such as restriction-modification systems seem to be associated with higher rates of recombination. Most importantly, recombined genomes show the footprints of natural selection in which recombined regions preferentially contain genes that can be related to specific ecological adaptations. Taken together, our results clarify the relative contributions of factors modulating homologous recombination and show evidence for a clear a role of this process in shaping microbial genomes and driving ecological adaptations.

Homologous recombination (HR) enables the exchange of genetic material between and within species. Recent studies suggest that this process plays a major role in the microevolution of microbial genomes, contributing to core genome homogenization and to the maintenance of cohesive population structures. However, we still have a very poor understanding of the possible adaptive roles of intraspecific HR and of the factors that determine its differential impact across clades and lifestyles. Here we used a unified methodological framework to assess HR in 338 complete genomes from 54 phylogenetically diverse and representative prokaryotic species, encompassing different lifestyles and a broad phylogenetic distribution. Our results indicate that lifestyle and presence of restriction-modification (RM) machineries are among the main factors shaping HR patterns, with symbionts and intracellular pathogens having the lowest HR levels. Similarly, the size of exchanged genomic fragments correlated with the presence of RM and competence machineries. Finally, genes exchanged by HR showed functional enrichments which could be related to adaptations to different environments and ecological strategies. Taken together, our results clarify the factors underlying HR impact and suggest important adaptive roles of genes exchanged through this mechanism. Our results also revealed that the extent of genetic exchange correlated with lifestyle and some genomic features. Moreover, the genes in exchanged regions were enriched for functions that reflected specific adaptations, supporting identification of HR as one of the main evolutionary mechanisms shaping prokaryotic core genomes.

Introduction to Genome Biology and Diversity

Organisms display astonishing levels of cell and molecular diversity, including genome size, shape, and architecture. In this chapter, we review how the genome can be viewed as both a structural and an informational unit of biological diversity and explicitly define our intended meaning of genetic information. A brief overview of the characteristic features of bacterial, archaeal, and eukaryotic cell types and viruses sets the stage for a review of the differences in organization, size, and packaging strategies of their genomes. We include a detailed review of genetic elements found outside the primary chromosomal structures, as these provide insights into how genomes are sometimes viewed as incomplete informational entities. Lastly, we reassess the definition of the genome in light of recent advancements in our understanding of the diversity of genomic structures and the mechanisms by which genetic information is expressed within the cell. Collectively, these topics comprise a good introduction to genome biology for the newcomer to the field and provide a valuable reference for those developing new statistical or computation methods in genomics. This review also prepares the reader for anticipated transformations in thinking as the field of genome biology progresses.

The Biogeography of Great Salt Lake Halophilic Archaea: Testing the Hypothesis of Avian Mechanical Carriers

Halophilic archaea inhabit hypersaline ecosystems globally, and genetically similar strains have been found in locales that are geographically isolated from one another. We sought to test the hypothesis that small salt crystals harboring halophilic archaea could be carried on bird feathers and that bird migration is a driving force of these distributions. In this study, we discovered that the American White Pelicans (AWPE) at Great Salt Lake soak in the hypersaline brine and accumulate salt crystals (halite) on their feathers. We cultured halophilic archaea from AWPE feathers and halite crystals. The microorganisms isolated from the lakeshore crystals were restricted to two genera: Halorubrum and Haloarcula, however, archaea from the feathers were strictly Haloarcula. We compared partial DNA sequence of the 16S rRNA gene from our cultivars with that of similar strains in the GenBank database. To understand the biogeography of genetically similar halophilic archaea, we studied the geographical locations of the sampling sites of the closest-matched species. An analysis of the environmental factors of each site pointed to salinity as the most important factor for selection. The geography of the sites was consistent with the location of the sub-tropical jet stream where birds typically migrate, supporting the avian dispersal hypothesis.

Copy number variation is associated with gene expression change in archaea

Genomic instability, although frequently deleterious, is also an important mechanism for microbial adaptation to environmental change. Although widely studied in bacteria, in archaea the effect of genomic instability on organism phenotypes and fitness remains unclear. Here we use DNA segmentation methods to detect and quantify genome-wide copy number variation (CNV) in large compendia of high-throughput datasets in a model archaeal species, Halobacterium salinarum. CNV hotspots were identified throughout the genome. Some hotspots were strongly associated with changes in gene expression, suggesting a mechanism for phenotypic innovation. In contrast, CNV hotspots in other genomic loci left expression unchanged, suggesting buffering of certain phenotypes. The correspondence of CNVs with gene expression was validated with strain- and condition-matched transcriptomics and DNA quantification experiments at specific loci. Significant correlation of CNV hotspot locations with the positions of known insertion sequence (IS) elements suggested a mechanism for generating genomic instability. Given the efficient recombination capabilities in H. salinarum despite stability at the single nucleotide level, these results suggest that genomic plasticity mediated by IS element activity can provide a source of phenotypic innovation in extreme environments.

Consensus assessment of the contamination level of publicly available cyanobacterial genomes

Publicly available genomes are crucial for phylogenetic and metagenomic studies, in which contaminating sequences can be the cause of major problems. This issue is expected to be especially important for Cyanobacteria because axenic strains are notoriously difficult to obtain and keep in culture. Yet, despite their great scientific interest, no data are currently available concerning the quality of publicly available cyanobacterial genomes. As reliably detecting contaminants is a complex task, we designed a pipeline combining six methods in a consensus strategy to assess the contamination level of 440 genome assemblies of Cyanobacteria. Two methods are based on published reference databases of ribosomal genes (SSU rRNA 16S and ribosomal proteins), one is indirectly based on a reference database of marker genes (CheckM), and three are based on complete genome analysis. Among those genome-wide methods, Kraken and DIAMOND blastx share the same reference database that we derived from Ensembl Bacteria, whereas CONCOCT does not require any reference database, instead relying on differences in DNA tetramer frequencies. Given that all the six methods appear to have their own strengths and limitations, we used the consensus of their rankings to infer that >5% of cyanobacterial genome assemblies are highly contaminated by foreign DNA (i.e., contaminants were detected by 5 or 6 methods). Our results will help researchers to check the quality of publicly available genomic data before use in their own analyses. Moreover, we argue that journals should make mandatory the submission of raw read data along with genome assemblies in order to facilitate the detection of contaminants in sequence databases.

Transmission in the Origins of Bacterial Diversity, From Ecotypes to Phyla

Any two lineages, no matter how distant they are now, began their divergence as one population splitting into two lineages that could coexist indefinitely. The rate of origin of higher-level taxa is therefore the product of the rate of speciation times the probability that two new species coexist long enough to reach a particular level of divergence. Here I have explored these two parameters of disparification in bacteria. Owing to low recombination rates, sexual isolation is not a necessary milestone of bacterial speciation. Rather, irreversible and indefinite divergence begins with ecological diversification, that is, transmission of a bacterial lineage to a new ecological niche, possibly to a new microhabitat but at least to new resources. Several algorithms use sequence data from a taxon of focus to identify phylogenetic groups likely to bear the dynamic properties of species. Identifying these newly divergent lineages allows us to characterize the genetic bases of speciation, as well as the ecological dimensions upon which new species diverge. Speciation appears to be least frequent when a given lineage has few new resources it can adopt, as exemplified by photoautotrophs, C1 heterotrophs, and obligately intracellular pathogens; speciation is likely most rapid for generalist heterotrophs. The genetic basis of ecological divergence may determine whether ecological divergence is irreversible and whether lineages will diverge indefinitely into the future. Long-term coexistence is most likely when newly divergent lineages utilize at least some resources not shared with the other and when the resources themselves will coexist into the remote future.

Genotypic and Lipid Analyses of Strains From the Archaeal Genus Halorubrum Reveal Insights Into Their Taxonomy, Divergence, and Population Structure

To gain a better understanding of how divergence occurs, and how taxonomy can benefit from studying natural populations, we isolated and examined 25 closely related Halorubrum strains obtained from different hypersaline communities and compared them to validly named species and other reference strains using five taxonomic study approaches: phylogenetic analysis using the 16S rRNA gene and multilocus sequencing analysis (MLSA), polar lipid profiles (PLP), average nucleotide identity (ANI) and DNA-DNA hybridization (DDH). 16S rRNA gene sequence could not differentiate the newly isolated strains from described species, while MLSA grouped strains into three major clusters. Two of those MLSA clusters distinguished candidates for new species. The third cluster with concatenated sequence identity equal to or greater than 97.5% was comprised of strains from Aran-Bidgol Lake (Iran) and solar salterns in Namibia and Spain, and two previously described species isolated from Mexico and Algeria. PLP and DDH analyses showed that Aran-Bidgol strains formed uniform populations, and that strains isolated from other geographic locations were heterogeneous and divergent, indicating that they may constitute different species. Therefore, applying only sequencing approaches and similarity cutoffs for circumscribing species may be too conservative, lumping concealed diversity into a single taxon. Further, our data support the interpretation that local populations experience unique evolutionary homogenization pressures, and once relieved of insular constraints (e.g., through migration) are free to diverge.

Comparative Analysis of Surface Layer Glycoproteins and Genes Involved in Protein Glycosylation in the Genus Haloferax

Within the Haloferax genus, both the surface (S)-layer protein, and the glycans that can decorate it, vary between species, which can potentially result in many different surface types, analogous to bacterial serotypes. This variation may mediate phenotypes, such as sensitivity to different viruses and mating preferences. Here, we describe S-layer glycoproteins found in multiple Haloferax strains and perform comparative genomics analyses of major and alternative glycosylation clusters of isolates from two coastal sites. We analyze the phylogeny of individual glycosylation genes and demonstrate that while the major glycosylation cluster tends to be conserved among closely related strains, the alternative cluster is highly variable. Thus, geographically- and genetically-related strains may exhibit diverse surface structures to such an extent that no two isolates present an identical surface profile.

Multilocus sequence analysis of homologous recombination and diversity in Arthrobacter sensu lato named species and glacier-inhabiting strains

Members of the bacterial genus Arthrobacter sensu lato are Gram-positive actinomycetes distributed worldwide and found in numerous environments including soil, water, glacier ice, and sewage. Homologous recombination is an important driving force in bacterial evolution, but its impact on Arthrobacter sensu lato evolution is poorly understood. We evaluated homologous recombination among 41 Arthrobacter sensu lato named species, using multilocus sequence analysis (MLSA). A high level of recombination was found, associated with strong diversification and a reticulate evolutionary pattern of Arthrobacter sensu lato. We also collected a total of 31 cold-adapted Arthrobacter sensu lato strains from two cold glaciers located in northwest China and two temperate glaciers in southwest China, and evaluated their diversity and population structure by MLSA. The glacier strains displayed high diversity, but rates of recombination among the four glacier groups were quite low, indicating that barriers to homologous recombination formed in the past among the populations on different glaciers. Our findings indicate that historical glaciation events shaped the contemporary distributions, taxonomic relationships, and phylogeographic patterns of Arthrobacter sensu lato species on glaciers.

A Rapid and Scalable Method for Multilocus Species Delimitation Using Bayesian Model Comparison and Rooted Triplets

Multilocus sequence data provide far greater power to resolve species limits than the single locus data typically used for broad surveys of clades. However, current statistical methods based on a multispecies coalescent framework are computationally demanding, because of the number of possible delimitations that must be compared and time-consuming likelihood calculations. New methods are therefore needed to open up the power of multilocus approaches to larger systematic surveys. Here, we present a rapid and scalable method that introduces 2 new innovations. First, the method reduces the complexity of likelihood calculations by decomposing the tree into rooted triplets. The distribution of topologies for a triplet across multiple loci has a uniform trinomial distribution when the 3 individuals belong to the same species, but a skewed distribution if they belong to separate species with a form that is specified by the multispecies coalescent. A Bayesian model comparison framework was developed and the best delimitation found by comparing the product of posterior probabilities of all triplets. The second innovation is a new dynamic programming algorithm for finding the optimum delimitation from all those compatible with a guide tree by successively analyzing subtrees defined by each node. This algorithm removes the need for heuristic searches used by current methods, and guarantees that the best solution is found and potentially could be used in other systematic applications. We assessed the performance of the method with simulated, published, and newly generated data. Analyses of simulated data demonstrate that the combined method has favorable statistical properties and scalability with increasing sample sizes. Analyses of empirical data from both eukaryotes and prokaryotes demonstrate its potential for delimiting species in real cases.

Analysis of the bacteriorhodopsin-producing haloarchaea reveals a core community that is stable over time in the salt crystallizers of Eilat, Israel

Stability of microbial communities can impact the ability of dispersed cells to colonize a new habitat. Saturated brines and their halophile communities are presumed to be steady state systems due to limited environmental perturbations. In this study, the bacteriorhodopsin-containing fraction of the haloarchaeal community from Eilat salt crystallizer ponds was sampled five times over 3 years. Analyses revealed the existence of a constant core as several OTUs were found repeatedly over the length of the study: OTUs comprising 52 % of the total cloned and sequenced PCR amplicons were found in every sample, and OTUs comprising 89 % of the total sequences were found in more than one, and often more than two samples. LIBSHUFF and UNIFRAC analyses showed statistical similarity between samples and Spearman's coefficient denoted significant correlations between OTU pairs, indicating non-random patterns in abundance and co-occurrence of detected OTUs. Further, changes in the detected OTUs were statistically linked to deviations in salinity. We interpret these results as indicating the existence of an ever-present core bacteriorhodopsin-containing Eilat crystallizer community that fluctuates in population densities, which are controlled by salinity rather than the extinction of some OTUs and their replacement through immigration and colonization.

Bacterial Speciation: Genetic Sweeps in Bacterial Species

One theory of bacterial speciation states that bacterial and animal species share the property of cohesion, meaning that diversity within a species is constrained. A new study provides direct evidence that genome-wide sweeps can limit diversity within bacterial species.

Oysters and Vibrios as a Model for Disease Dynamics in Wild Animals

Disease dynamics in the wild are influenced by a number of ecological and evolutionary factors not addressed by traditional laboratory-based characterization of pathogens. Here we propose the oyster, Crassostrea gigas, as a model for studying the interaction of the environment, bacterial pathogens, and the host in disease dynamics. We show that an important first step is to ask whether the functional unit of pathogenesis is a bacterial clone, a population, or a consortium in order to assess triggers of disease outbreaks and devise appropriate monitoring tools. Moreover, the development of specific-pathogen-free (SPF) oysters has enabled assessment of the infection process under natural conditions. Finally, recent results show the importance of microbial interactions and host genetics in determining oyster health and disease.

Contributions of ancestral inter-species recombination to the genetic diversity of extant Streptomyces lineages

Streptomyces species produce many important antibiotics and have a crucial role in soil nutrient cycling. However, their evolutionary history remains poorly characterized. We have evaluated the impact of homologous recombination on the evolution of Streptomyces using multi-locus sequence analysis of 234 strains that represent at least 11 species clusters. Evidence of inter-species recombination is widespread but not uniform within the genus and levels of mosaicism vary between species clusters. Most phylogenetically incongruent loci are monophyletic at the scale of species clusters and their subclades, suggesting that these recombination events occurred in shared ancestral lineages. Further investigation of two mosaic species clusters suggests that genes acquired by inter-species recombination may have become fixed in these lineages during periods of demographic expansion; implicating a role for phylogeography in determining contemporary patterns of genetic diversity. Only by examining the phylogeny at the scale of the genus is apparent that widespread phylogenetically incongruent loci in Streptomyces are derived from a far smaller number of ancestral inter-species recombination events.

Recombination Does Not Hinder Formation or Detection of Ecological Species of Synechococcus Inhabiting a Hot Spring Cyanobacterial Mat

Recent studies of bacterial speciation have claimed to support the biological species concept—that reduced recombination is required for bacterial populations to diverge into species. This conclusion has been reached from the discovery that ecologically distinct clades show lower rates of recombination than that which occurs among closest relatives. However, these previous studies did not attempt to determine whether the more-rapidly recombining close relatives within the clades studied may also have diversified ecologically, without benefit of sexual isolation. Here we have measured the impact of recombination on ecological diversification within and between two ecologically distinct clades (A and B') of Synechococcus in a hot spring microbial mat in Yellowstone National Park, using a cultivation-free, multi-locus approach. Bacterial artificial chromosome (BAC) libraries were constructed from mat samples collected at 60°C and 65°C. Analysis of multiple linked loci near Synechococcus 16S rRNA genes showed little evidence of recombination between the A and B' lineages, but a record of recombination was apparent within each lineage. Recombination and mutation rates within each lineage were of similar magnitude, but recombination had a somewhat greater impact on sequence diversity than mutation, as also seen in many other bacteria and archaea. Despite recombination within the A and B' lineages, there was evidence of ecological diversification within each lineage. The algorithm Ecotype Simulation identified sequence clusters consistent with ecologically distinct populations (ecotypes), and several hypothesized ecotypes were distinct in their habitat associations and in their adaptations to different microenvironments. We conclude that sexual isolation is more likely to follow ecological divergence than to precede it. Thus, an ecology-based model of speciation appears more appropriate than the biological species concept for bacterial and archaeal diversification.

Halorubrum halodurans sp. nov., an extremely halophilic archaeon isolated from a hypersaline lake

Two extremely halophilic archaea, strains Cb34T and C170, belonging to the genus Halorubrum, were isolated from the brine of the hypersaline lake Aran-Bidgol in Iran. Cells of the two strains were motile, pleomorphic rods, stained Gram-variable and produced red-pigmented colonies. Strains Cb34T and C170 required 25 % (w/v) salts, pH 7.0 and 37 °C for optimal growth under aerobic conditions; 0.3 M Mg2+ was required. Cells of both isolates were lysed in distilled water and hypotonic treatment with < 10 % NaCl provoked cell lysis. Phylogenetic analysis based on 16S rRNA gene sequence similarities showed that these two strains were closely related to Halorubrum cibi B31T (98.8 %) and other members of the genus Halorubrum. In addition, studies based on the rpoB' gene revealed that strains Cb34T and C170 are placed among the species of Halorubrum and are closely related to Halorubrum cibi B31T, with rpoB' gene sequence similarity less than or equal to 95.7 %. The polar lipid patterns of both strains consisted of phosphatidylglycerol, phosphatidylglycerol phosphate methyl ester, phosphatidylglycerol sulfate and sulfated mannosyl glucosyl diether. The DNA G+C content was 62.1-62.4 mol%. DNA-DNA hybridization studies confirmed that strains Cb34T and C170 constitute a distinct species. Data obtained in this study show that the two strains represent a novel species, for which the name Halorubrum halodurans sp. nov. is proposed. The type strain is Cb34T ( = CECT 8745T = IBRC-M 10233T).

Microbial Speciation

What are species? How do they arise? These questions are not easy to answer and have been particularly controversial in microbiology. Yet, for those microbiologists studying environmental questions or dealing with clinical issues, the ability to name and recognize species, widely considered the fundamental units of ecology, can be practically useful. On a more fundamental level, the speciation problem, the focus here, is more mechanistic and conceptual. What is the origin of microbial species, and what evolutionary and ecological mechanisms keep them separate once they begin to diverge? To what extent are these mechanisms universal across diverse types of microbes, and more broadly across the entire the tree of life? Here, we propose that microbial speciation must be viewed in light of gene flow, which defines units of genetic similarity, and of natural selection, which defines units of phenotype and ecological function. We discuss to what extent ecological and genetic units overlap to form cohesive populations in the wild, based on recent evolutionary modeling and population genomics studies. These studies suggest a continuous "speciation spectrum," which microbial populations traverse in different ways depending on their balance of gene flow and natural selection.

Inferring Speciation Processes from Patterns of Natural Variation in Microbial Genomes

Microbial species concepts have long been the focus of contentious debate, fueled by technological limitations to the genetic resolution of species, by the daunting task of investigating phenotypic variation among individual microscopic organisms, and by a lack of understanding of gene flow in reproductively asexual organisms that are prone to promiscuous horizontal gene transfer. Population genomics, the emerging approach of analyzing the complete genomes of a multitude of closely related organisms, is poised to overcome these limitations by providing a window into patterns of genome variation revealing the evolutionary processes through which species diverge. This new approach is more than just an extension of previous multilocus sequencing technologies, in that it provides a comprehensive view of interacting evolutionary processes. Here we argue that the application of population genomic tools in a rigorous population genetic framework will help to identify the processes of microbial speciation and ultimately lead to a general species concept based on the unique biology and ecology of microorganisms.

Microbial diversity of hypersaline environments: a metagenomic approach

Recent studies based on metagenomics and other molecular techniques have permitted a detailed knowledge of the microbial diversity and metabolic activities of microorganisms in hypersaline environments. The current accepted model of community structure in hypersaline environments is that the square archaeon Haloquadratum waslbyi, the bacteroidete Salinibacter ruber and nanohaloarchaea are predominant members at higher salt concentrations, while more diverse archaeal and bacterial taxa are observed in habitats with intermediate salinities. Additionally, metagenomic studies may provide insight into the isolation and characterization of the principal microbes in these habitats, such as the recently described gammaproteobacterium Spiribacter salinus.

Horizontal gene transfer, dispersal and haloarchaeal speciation

The Halobacteria are a well-studied archaeal class and numerous investigations are showing how their diversity is distributed amongst genomes and geographic locations. Evidence indicates that recombination between species continuously facilitates the arrival of new genes, and within species, it is frequent enough to spread acquired genes amongst all individuals in the population. To create permanent independent diversity and generate new species, barriers to recombination are probably required. The data support an interpretation that rates of evolution (e.g., horizontal gene transfer and mutation) are faster at creating geographically localized variation than dispersal and invasion are at homogenizing genetic differences between locations. Therefore, we suggest that recurrent episodes of dispersal followed by variable periods of endemism break the homogenizing forces of intrapopulation recombination and that this process might be the principal stimulus leading to divergence and speciation in Halobacteria.

Biofilms formed by the archaeon Haloferax volcanii exhibit cellular differentiation and social motility, and facilitate horizontal gene transfer

BACKGROUND

Archaea share a similar microbial lifestyle with bacteria, and not surprisingly then, also exist within matrix-enclosed communities known as biofilms. Advances in biofilm biology have been made over decades for model bacterial species, and include characterizations of social behaviors and cellular differentiation during biofilm development. Like bacteria, archaea impact ecological and biogeochemical systems. However, the biology of archaeal biofilms is only now being explored. Here, we investigated the development, composition and dynamics of biofilms formed by the haloarchaeon Haloferax volcanii DS2.

RESULTS

Biofilms were cultured in static liquid and visualized with fluorescent cell membrane dyes and by engineering cells to express green fluorescent protein (GFP). Analysis by confocal scanning laser microscopy showed that H. volcanii cells formed microcolonies within 24 h, which developed into larger clusters by 48 h and matured into flake-like towers often greater than 100 μm in height after 7 days. To visualize the extracellular matrix, biofilms formed by GFP-expressing cells were stained with concanavalin A, DAPI, Congo red and thioflavin T. Stains colocalized with larger cellular structures and indicated that the extracellular matrix may contain a combination of polysaccharides, extracellular DNA and amyloid protein. Following a switch to biofilm growth conditions, a sub-population of cells differentiated into chains of long rods sometimes exceeding 25 μm in length, compared to their planktonic disk-shaped morphology. Time-lapse photography of static liquid biofilms also revealed wave-like social motility. Finally, we quantified gene exchange between biofilm cells, and found that it was equivalent to the mating frequency of a classic filter-based experimental method.

CONCLUSIONS

The developmental processes, functional properties and dynamics of H. volcanii biofilms provide insight on how haloarchaeal species might persist, interact and exchange DNA in natural communities. H. volcanii demonstrates some biofilm phenotypes similar to bacterial biofilms, but also has interesting phenotypes that may be unique to this organism or to this class of organisms, including changes in cellular morphology and an unusual form of social motility. Because H. volcanii has one of the most advanced genetic systems for any archaeon, the phenotypes reported here may promote the study of genetic and developmental processes in archaeal biofilms.

Genomic heterogeneity and ecological speciation within one subspecies of Bacillus subtilis

Closely related bacterial genomes usually differ in gene content, suggesting that nearly every strain in nature may be ecologically unique. We have tested this hypothesis by sequencing the genomes of extremely close relatives within a recognized taxon and analyzing the genomes for evidence of ecological distinctness. We compared the genomes of four Death Valley isolates plus the laboratory strain W23, all previously classified as Bacillus subtilis subsp. spizizenii and hypothesized through multilocus analysis to be members of the same ecotype (an ecologically homogeneous population), named putative ecotype 15 (PE15). These strains showed a history of positive selection on amino acid sequences in 38 genes. Each of the strains was under a different regimen of positive selection, suggesting that each strain is ecologically unique and represents a distinct ecological speciation event. The rate of speciation appears to be much faster than can be resolved with multilocus sequencing. Each PE15 strain contained unique genes known to confer a function for bacteria. Remarkably, no unique gene conferred a metabolic system or subsystem function that was not already present in all the PE15 strains sampled. Thus, the origin of ecotypes within this clade shows no evidence of qualitative divergence in the set of resources utilized. Ecotype formation within this clade is consistent with the nanoniche model of bacterial speciation, in which ecotypes use the same set of resources but in different proportions, and genetic cohesion extends beyond a single ecotype to the set of ecotypes utilizing the same resources.

Inteins as indicators of gene flow in the halobacteria

This research uses inteins, a type of mobile genetic element, to infer patterns of gene transfer within the Halobacteria. We surveyed 118 genomes representing 26 genera of Halobacteria for intein sequences. We then used the presence-absence profile, sequence similarity and phylogenies from the inteins recovered to explore how intein distribution can provide insight on the dynamics of gene flow between closely related and divergent organisms. We identified 24 proteins in the Halobacteria that have been invaded by inteins at some point in their evolutionary history, including two proteins not previously reported to contain an intein. Furthermore, the size of an intein is used as a heuristic for the phase of the intein's life cycle. Larger size inteins are assumed to be the canonical two domain inteins, consisting of self-splicing and homing endonuclease domains (HEN); smaller sizes are assumed to have lost the HEN domain. For many halobacterial groups the consensus phylogenetic signal derived from intein sequences is compatible with vertical inheritance or with a strong gene transfer bias creating these clusters. Regardless, the coexistence of intein-free and intein-containing alleles reveal ongoing transfer and loss of inteins within these groups. Inteins were frequently shared with other Euryarchaeota and among the Bacteria, with members of the Cyanobacteria (Cyanothece, Anabaena), Bacteriodetes (Salinibacter), Betaproteobacteria (Delftia, Acidovorax), Firmicutes (Halanaerobium), Actinobacteria (Longispora), and Deinococcus-Thermus-group.

Polyploidy in haloarchaea: advantages for growth and survival

The investigated haloarchaeal species, Halobacterium salinarum, Haloferax mediterranei, and H. volcanii, have all been shown to be polyploid. They contain several replicons that have independent copy number regulation, and most have a higher copy number during exponential growth phase than in stationary phase. The possible evolutionary advantages of polyploidy for haloarchaea, most of which have experimental support for at least one species, are discussed. These advantages include a low mutation rate and high resistance toward X-ray irradiation and desiccation, which depend on homologous recombination. For H. volcanii, it has been shown that gene conversion operates in the absence of selection, which leads to the equalization of genome copies. On the other hand, selective forces might lead to heterozygous cells, which have been verified in the laboratory. Additional advantages of polyploidy are survival over geological times in halite deposits as well as at extreme conditions on earth and at simulated Mars conditions. Recently, it was found that H. volcanii uses genomic DNA as genetic material and as a storage polymer for phosphate. In the absence of phosphate, H. volcanii dramatically decreases its genome copy number, thereby enabling cell multiplication, but diminishing the genetic advantages of polyploidy. Stable storage of phosphate is proposed as an alternative driving force for the emergence of DNA in early evolution. Several additional potential advantages of polyploidy are discussed that have not been addressed experimentally for haloarchaea. An outlook summarizes selected current trends and possible future developments.

Single-cell genomics reveals hundreds of coexisting subpopulations in wild Prochlorococcus

Extensive genomic diversity within coexisting members of a microbial species has been revealed through selected cultured isolates and metagenomic assemblies. Yet, the cell-by-cell genomic composition of wild uncultured populations of co-occurring cells is largely unknown. In this work, we applied large-scale single-cell genomics to study populations of the globally abundant marine cyanobacterium Prochlorococcus. We show that they are composed of hundreds of subpopulations with distinct "genomic backbones," each backbone consisting of a different set of core gene alleles linked to a small distinctive set of flexible genes. These subpopulations are estimated to have diverged at least a few million years ago, suggesting ancient, stable niche partitioning. Such a large set of coexisting subpopulations may be a general feature of free-living bacterial species with huge populations in highly mixed habitats.

DNA as a phosphate storage polymer and the alternative advantages of polyploidy for growth or survival

Haloferax volcanii uses extracellular DNA as a source for carbon, nitrogen, and phosphorous. However, it can also grow to a limited extend in the absence of added phosphorous, indicating that it contains an intracellular phosphate storage molecule. As Hfx. volcanii is polyploid, it was investigated whether DNA might be used as storage polymer, in addition to its role as genetic material. It could be verified that during phosphate starvation cells multiply by distributing as well as by degrading their chromosomes. In contrast, the number of ribosomes stayed constant, revealing that ribosomes are distributed to descendant cells, but not degraded. These results suggest that the phosphate of phosphate-containing biomolecules (other than DNA and RNA) originates from that stored in DNA, not in rRNA. Adding phosphate to chromosome depleted cells rapidly restores polyploidy. Quantification of desiccation survival of cells with different ploidy levels showed that under phosphate starvation Hfx. volcanii diminishes genetic advantages of polyploidy in favor of cell multiplication. The consequences of the usage of genomic DNA as phosphate storage polymer are discussed as well as the hypothesis that DNA might have initially evolved in evolution as a storage polymer, and the various genetic benefits evolved later.

Population and genomic analysis of the genus Halorubrum

The Halobacteria are known to engage in frequent gene transfer and homologous recombination. For stably diverged lineages to persist some checks on the rate of between lineage recombination must exist. We surveyed a group of isolates from the Aran-Bidgol endorheic lake in Iran and sequenced a selection of them. Multilocus Sequence Analysis (MLSA) and Average Nucleotide Identity (ANI) revealed multiple clusters (phylogroups) of organisms present in the lake. Patterns of intein and Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs) presence/absence and their sequence similarity, GC usage along with the ANI and the identities of the genes used in the MLSA revealed that two of these clusters share an exchange bias toward others in their phylogroup while showing reduced rates of exchange with other organisms in the environment. However, a third cluster, composed in part of named species from other areas of central Asia, displayed many indications of variability in exchange partners, from within the lake as well as outside the lake. We conclude that barriers to gene exchange exist between the two purely Aran-Bidgol phylogroups, and that the third cluster with members from other regions is not a single population and likely reflects an amalgamation of several populations.

Evidence from phylogenetic and genome fingerprinting analyses suggests rapidly changing variation in Halorubrum and Haloarcula populations

Halobacteria require high NaCl concentrations for growth and are the dominant inhabitants of hypersaline environments above 15% NaCl. They are well-documented to be highly recombinogenic, both in frequency and in the range of exchange partners. In this study, we examine the genetic and genomic variation of cultured, naturally co-occurring environmental populations of Halobacteria. Sequence data from multiple loci (~2500 bp) identified many closely and more distantly related strains belonging to the genera Halorubrum and Haloarcula. Genome fingerprinting using a random priming PCR amplification method to analyze these isolates revealed diverse banding patterns across each of the genera and surprisingly even for isolates that are identical at the nucleotide level for five protein coding sequenced loci. This variance in genome structure even between identical multilocus sequence analysis (MLSA) haplotypes indicates that accumulation of genomic variation is rapid: faster than the rate of third codon substitutions.

Ordering microbial diversity into ecologically and genetically cohesive units

We propose that microbial diversity must be viewed in light of gene flow and selection, which define units of genetic similarity, and of phenotype and ecological function, respectively. We discuss to what extent ecological and genetic units overlap to form cohesive populations in the wild, based on recent evolutionary modeling and on evidence from some of the first microbial populations studied with genomics. These show that if recombination is frequent and selection moderate, ecologically adaptive mutations or genes can spread within populations independently of their original genomic background (gene-specific sweeps). Alternatively, if the effect of recombination is smaller than selection, genome-wide selective sweeps should occur. In both cases, however, distinct units of overlapping ecological and genotypic similarity will form if microgeographic separation, likely involving ecological tradeoffs, induces barriers to gene flow. These predictions are supported by (meta)genomic data, which suggest that a 'reverse ecology' approach, in which genomic and gene flow information is used to make predictions about the nature of ecological units, is a powerful approach to ordering microbial diversity.

Patterns of microbial diversity along a salinity gradient in the Guerrero Negro solar saltern, Baja CA Sur, Mexico

The goal of this study was to use environmental sequencing of 16S rRNA and bop genes to compare the diversity of planktonic bacteria and archaea across ponds with increasing salinity in the Exportadora de Sal (ESSA) evaporative saltern in Guerrero Negro, Baja CA S., Mexico. We hypothesized that diverse communities of heterotrophic bacteria and archaea would be found in the ESSA ponds, but that bacterial diversity would decrease relative to archaea at the highest salinities. Archaeal 16S rRNA diversity was higher in Ponds 11 and 12 (370 and 380 g l(-1) total salts, respectively) compared to Pond 9 (180 g l(-1) total salts). Both Pond 11 and 12 communities had high representation (47 and 45% of clones, respectively) by Haloquadratum walsbyi-like (99% similarity) lineages. The archaeal community in Pond 9 was dominated (79%) by a single uncultured phylotype with 99% similarity to sequences recovered from the Sfax saltern in Tunisia. This pattern was mirrored in bop gene diversity with greater numbers of highly supported phylotypes including many Haloquadratum-like sequences from the two highest salinity ponds. In Pond 9, most bop sequences, were not closely related to sequences in databases. Bacterial 16S rRNA diversity was higher than archaeal in both Pond 9 and Pond 12 samples, but not Pond 11, where a non-Salinibacter lineage within the Bacteroidetes >98% similar to environmental clones recovered from Lake Tuz in Turkey and a saltern in Chula Vista, CA was most abundant (69% of community). This OTU was also the most abundant in Pond 12, but only represented 14% of clones in the more diverse pond. The most abundant OTU in Pond 9 (33% of community) was 99% similar to an uncultured gammaproteobacterial clone from the Salton Sea. Results suggest that the communities of saltern bacteria and archaea vary even in ponds with similar salinity and further investigation into the ecology of diverse, uncultured halophile communities is warranted.

High level of intergenera gene exchange shapes the evolution of haloarchaea in an isolated Antarctic lake

Deep Lake in Antarctica is a globally isolated, hypersaline system that remains liquid at temperatures down to -20 °C. By analyzing metagenome data and genomes of four isolates we assessed genome variation and patterns of gene exchange to learn how the lake community evolved. The lake is completely and uniformly dominated by haloarchaea, comprising a hierarchically structured, low-complexity community that differs greatly to temperate and tropical hypersaline environments. The four Deep Lake isolates represent distinct genera (∼85% 16S rRNA gene similarity and ∼73% genome average nucleotide identity) with genomic characteristics indicative of niche adaptation, and collectively account for ∼72% of the cellular community. Network analysis revealed a remarkable level of intergenera gene exchange, including the sharing of long contiguous regions (up to 35 kb) of high identity (∼100%). Although the genomes of closely related Halobacterium, Haloquadratum, and Haloarcula (>90% average nucleotide identity) shared regions of high identity between species or strains, the four Deep Lake isolates were the only distantly related haloarchaea to share long high-identity regions. Moreover, the Deep Lake high-identity regions did not match to any other hypersaline environment metagenome data. The most abundant species, tADL, appears to play a central role in the exchange of insertion sequences, but not the exchange of high-identity regions. The genomic characteristics of the four haloarchaea are consistent with a lake ecosystem that sustains a high level of intergenera gene exchange while selecting for ecotypes that maintain sympatric speciation. The peculiarities of this polar system restrict which species can grow and provide a tempo and mode for accentuating gene exchange.

From green to red: horizontal gene transfer of the phycoerythrin gene cluster between Planktothrix strains

Horizontal gene transfer is common in cyanobacteria, and transfer of large gene clusters may lead to acquisition of new functions and conceivably niche adaption. In the present study, we demonstrate that horizontal gene transfer between closely related Planktothrix strains can explain the production of the same oligopeptide isoforms by strains of different colors. Comparison of the genomes of eight Planktothrix strains revealed that strains producing the same oligopeptide isoforms are closely related, regardless of color. We have investigated genes involved in the synthesis of the photosynthetic pigments phycocyanin and phycoerythrin, which are responsible for green and red appearance, respectively. Sequence comparisons suggest the transfer of a functional phycoerythrin gene cluster generating a red phenotype in a strain that is otherwise more closely related to green strains. Our data show that the insertion of a DNA fragment containing the 19.7-kb phycoerythrin gene cluster has been facilitated by homologous recombination, also replacing a region of the phycocyanin operon. These findings demonstrate that large DNA fragments spanning entire functional gene clusters can be effectively transferred between closely related cyanobacterial strains and result in a changed phenotype. Further, the results shed new light on the discussion of the role of horizontal gene transfer in the sporadic distribution of large gene clusters in cyanobacteria, as well as the appearance of red and green strains.

Multilocus sequence analysis of Thermoanaerobacter isolates reveals recombining, but differentiated, populations from geothermal springs of the Uzon Caldera, Kamchatka, Russia

Thermal environments have island-like characteristics and provide a unique opportunity to study population structure and diversity patterns of microbial taxa inhabiting these sites. Strains having ≥98% 16S rRNA gene sequence similarity to the obligately anaerobic Firmicutes Thermoanaerobacter uzonensis were isolated from seven geothermal springs, separated by up to 1600 m, within the Uzon Caldera (Kamchatka, Russian Far East). The intraspecies variation and spatial patterns of diversity for this taxon were assessed by multilocus sequence analysis (MLSA) of 106 strains. Analysis of eight protein-coding loci (gyrB, lepA, leuS, pyrG, recA, recG, rplB, and rpoB) revealed that all loci were polymorphic and that nucleotide substitutions were mostly synonymous. There were 148 variable nucleotide sites across 8003 bp concatenates of the protein-coding loci. While pairwise F_ST values indicated a small but significant level of genetic differentiation between most subpopulations, there was a negligible relationship between genetic divergence and spatial separation. Strains with the same allelic profile were only isolated from the same hot spring, occasionally from consecutive years, and single locus variant (SLV) sequence types were usually derived from the same spring. While recombination occurred, there was an “epidemic” population structure in which a particular T. uzonensis sequence type rose in frequency relative to the rest of the population. These results demonstrate spatial diversity patterns for an anaerobic bacterial species in a relative small geographic location and reinforce the view that terrestrial geothermal springs are excellent places to look for biogeographic diversity patterns regardless of the involved distances.

Quantifying homologous replacement of loci between haloarchaeal species

In vitro studies of the haloarchaeal genus Haloferax have demonstrated their ability to frequently exchange DNA between species, whereas rates of homologous recombination estimated from natural populations in the genus Halorubrum are high enough to maintain random association of alleles between five loci. To quantify the effects of gene transfer and recombination of commonly held (relaxed core) genes during the evolution of the class Halobacteria (haloarchaea), we reconstructed the history of 21 genomes representing all major groups. Using a novel algorithm and a concatenated ribosomal protein phylogeny as a reference, we created a directed horizontal genetic transfer (HGT) network of contemporary and ancestral genomes. Gene order analysis revealed that 90% of testable HGTs were by direct homologous replacement, rather than nonhomologous integration followed by a loss. Network analysis revealed an inverse log-linear relationship between HGT frequency and ribosomal protein evolutionary distance that is maintained across the deepest divergences in Halobacteria. We use this mathematical relationship to estimate the total transfers and amino acid substitutions delivered by HGTs in each genome, providing a measure of chimerism. For the relaxed core genes of each genome, we conservatively estimate that 11–20% of their evolution occurred in other haloarchaea. Our findings are unexpected, because the transfer and homologous recombination of relaxed core genes between members of the class Halobacteria disrupts the coevolution of genes; however, the generation of new combinations of divergent but functionally related genes may lead to adaptive phenotypes not available through cumulative mutations and recombination within a single population.

Horizontal gene transfer and the evolution of bacterial and archaeal population structure

Many bacterial and archaeal lineages have a history of extensive and ongoing horizontal gene transfer and loss, as evidenced by the large differences in genome content even among otherwise closely related isolates. How ecologically cohesive populations might evolve and be maintained under such conditions of rapid gene turnover has remained controversial. Here we synthesize recent literature demonstrating the importance of habitat and niche in structuring horizontal gene transfer. This leads to a model of ecological speciation via gradual genetic isolation triggered by differential habitat-association of nascent populations. Further, we hypothesize that subpopulations can evolve through local gene-exchange networks by tapping into a gene pool that is adaptive towards local, continuously changing organismic interactions and is, to a large degree, responsible for the observed rapid gene turnover. Overall, these insights help to explain how bacteria and archaea form populations that display both ecological cohesion and high genomic diversity.