Towards a conceptual and operational union of bacterial systematics, ecology, and evolution

To completely understand the ecology of a bacterial community, we need to identify its ecologically distinct populations (ecotypes). The greatest promise for enumerating a community's constituent ecotypes is held by molecular approaches that identify bacterial ecotypes as DNA sequence clusters. These approaches succeed when ecotypes correspond with sequence clusters, but some models of bacterial speciation predict a one-to-many and others a many-to-one relationship between ecotypes and sequence clusters. A further challenge is that sequence-based phylogenies often contain a hierarchy of clusters and subclusters within clusters, and there is no widely accepted theory to guide systematists and ecologists to the size of cluster most likely to correspond to ecotypes. While present systematics attempts to use universal thresholds of sequence divergence to help demarcate species, the recently developed 'community phylogeny' approach assumes no universal thresholds, but demarcates ecotypes based on the analysis of a lineage's evolutionary dynamics. Theory-based approaches like this one can give a conceptual framework as well as operational criteria for hypothesizing the identity and membership of ecotypes from sequence data; ecology-based approaches can then confirm that the putative ecotypes are actually ecologically distinct. Bacterial ecotypes that are demonstrated to have a history of coexistence as ecologically distinct lineages (based on sequence analysis) and as a prognosis of future coexistence (based on ecological differences), are the fundamental units of bacterial ecology and evolution, and should be recognized by bacterial systematics.     

Bacterial species and speciation

Bacteria are profoundly different from eukaryotes in their patterns of genetic exchange. Nevertheless, ecological diversity is organized in the same way across all of life: individual organisms fall into more less discrete clusters on the basis of their phenotypic, ecological, and DNA sequence characteristics. Each sequence cluster in the bacterial world appears to correspond to an "ecotype," defined as a population of cells in the same ecological niche, which would all be out-competed by any adaptive mutant coming from the population. Ecotypes, so defined, share many of the dynamic properties attributed to eukaryotic species: genetic diversity within an ecotype is limited by a force of cohesion (in this case, periodic selection); different ecotypes are free to diverge without constraint from one another; and ecotypes are ecologically distinct. Also, ecotypes can be discovered and classified as DNA sequence clusters, even when we are ignorant of their ecology. Owing to the rarity and promiscuity of bacterial genetic exchange, speciation in the bacterial world is expected to be much less constrained than in the world of animals and plants.     

A species concept for bacteria based on adaptive divergence

Bacterial strains are currently grouped into species based on overall genomic similarity and sharing of phenotypes deemed ecologically important. Many believe this polyphasic taxonomy is in need of revision because it lacks grounding in evolutionary theory, and boundaries between species are arbitrary. Recent taxonomy efforts using multilocus sequence typing (MLST) data are based on the identification of distinct phylogenetic clusters. However, these approaches face the problem of deciding the phylogenetic level at which clusters are representative of evolutionary or taxonomically distinct units. In this review, I propose classifying two phylogenetic clusters as separate species only when they have statistically significantly diverged as a result of adaptive evolution. More than a method for classification, the concept of adaptive divergence can be used in a 'reverse ecology' approach to identify lineages that are in the process of speciation or genes involved in initial adaptive divergence.     

A phylogenetic perspective on habitat shifts and diversity in the North American Enallagma damselflies

Community ecologists are increasingly aware that the regional history of taxon diversification can have an important influence on community structure. Likewise, systematists recognize that ecological context can have an important influence on the processes of speciation and extinction that create patterns of descent. We present a phylogenetic analysis of 33 species of a North American radiation of damselflies (Zygoptera: Coenagrionidae: Enallagma Selys), which have been well studied ecologically, to elucidate the evolutionary mechanisms that have contributed to differences in diversity between larval habitats (lakes with and without fish predators). Analysis of molecular variation in 842 bp of the mitochondrial cytochrome oxidase I and II subunit and of the intervening Leu-tRNA and 37 morphological characters resulted in three well-defined clades that are only partially congruent with previous phylogenetic hypotheses. Molecular and morphological data partitions were significantly incongruent (p < .01). Lack of haplotype monophyly within species and small amounts of sequence divergence (< 1%) between related species in three of the four clades suggest that recent, and parallel, speciation has been an important source of community diversity. Reconstruction of habitat preference over the phylogeny suggests that the greater species diversity in fish-containing lake habitats reflects the recency of shifts into the fishless lake habit, although a difference in speciation or extinction rates between the two habitats is difficult to exclude as an additional mechanism.     

Origins of bacterial diversity through horizontal genetic transfer and adaptation to new ecological niches

Horizontal genetic transfer (HGT) has played an important role in bacterial evolution at least since the origins of the bacterial divisions, and HGT still facilitates the origins of bacterial diversity, including diversity based on antibiotic resistance. Adaptive HGT is aided by unique features of genetic exchange in bacteria such as the promiscuity of genetic exchange and the shortness of segments transferred. Genetic exchange rates are limited by the genetic and ecological similarity of organisms. Adaptive transfer of genes is limited to those that can be transferred as a functional unit, provide a niche-transcending adaptation, and are compatible with the architecture and physiology of other organisms. Horizontally transferred adaptations may bring about fitness costs, and natural selection may ameliorate these costs. The origins of ecological diversity can be analyzed by comparing the genomes of recently divergent, ecologically distinct populations, which can be discovered as sequence clusters. Such genome comparisons demonstrate the importance of HGT in ecological diversification. Newly divergent populations cannot be discovered as sequence clusters when their ecological differences are coded by plasmids, as is often the case for antibiotic resistance; the discovery of such populations requires a screen for plasmid-coded functions. This paper reviews the features of bacterial genetics that allow HGT, the similarities between organisms that foster HGT between them, the limits to the kinds of adaptations that can be transferred, and amelioration of fitness costs associated with HGT; the paper also reviews approaches to discover the origins of new, ecologically distinct bacterial populations and the role that HGT plays in their founding.     

Evolving Concepts of Bacterial Species

The same evolutionary forces that cause diversification in sexual eukaryotes are expected to cause diversification in bacteria. However, in bacteria, the wider variety of mechanisms for gene exchange (or lack thereof) increases the range of expected diversity patterns compared to those of sexual organisms. Two parallel concepts for bacterial speciation have developed, based on ecological divergence or barriers to recombination in turn. Recent evidence from DNA sequence data shows that both processes can generate independently evolving groups that are equivalent to sexual species and that represent separate arenas within which recombination (when it occurs), selection and drift occur. It remains unclear, however, how often different processes act in concert to generate simple units of diversity, or whether a more complex model of diversity is required, specifying hierarchical levels at which different cohesive processes operate. We advocate an integrative approach that evaluates the effects of multiple evolutionary forces on diversity patterns. There is also great potential for laboratory studies of bacterial evolution that test evolutionary mechanisms inferred from population genetic analyses of multi-locus and genome sequence data.     

Ecological niche modelling as an exploratory tool for identifying species limits: an example based on Mexican muroid rodents

Niche conservatism theory suggests that recently diverged sister species share the same ecological niche. However, if the ecological niche evolves as part of the speciation process, the ecological pattern could be useful for recognizing cryptic species. In a broad sense systematists agree that the niche characters could be used for species differentiation. However, to date such characters have been ignored. We used the genetic algorithm for rule‐set production for modelling the ecological niche as a means of inferring ecological divergence in allopatric populations of muroid rodents for which taxonomic identity is uncertain. Our results show that niche differentiation is significant in most of the identified phylogroups. The differentiation is likely associated with natural evolutionary units, which can be identified by applying species concepts based on phylogenetic and ecological patterns (e.g. phylogenetic, cohesive, evolutionary). Even so, the role of the niche partition within phylogenetic reconstruction may be a limited one.     

Natural product biosynthetic gene diversity in geographically distinct soil microbiomes

The number of bacterial species estimated to exist on Earth has increased dramatically in recent years. This newly recognized species diversity has raised the possibility that bacterial natural product biosynthetic diversity has also been significantly underestimated by previous culture-based studies. Here, we compare 454-pyrosequenced nonribosomal peptide adenylation domain, type I polyketide ketosynthase domain, and type II polyketide ketosynthase alpha gene fragments amplified from cosmid libraries constructed using DNA isolated from three different arid soils. While 16S rRNA gene sequence analysis indicates these cloned metagenomes contain DNA from similar distributions of major bacterial phyla, we found that they contain almost completely distinct collections of secondary metabolite biosynthetic gene sequences. When grouped at 85% identity, only 1.5% of the adenylation domain, 1.2% of the ketosynthase, and 9.3% of the ketosynthase alpha sequence clusters contained sequences from all three metagenomes. Although there is unlikely to be a simple correlation between biosynthetic gene sequence diversity and the diversity of metabolites encoded by the gene clusters in which these genes reside, our analysis further suggests that sequences in one soil metagenome are so distantly related to sequences in another metagenome that they are, in many cases, likely to arise from functionally distinct gene clusters. The marked differences observed among collections of biosynthetic genes found in even ecologically similar environments suggest that prokaryotic natural product biosynthesis diversity is, like bacterial species diversity, potentially much larger than appreciated from culture-based studies.     

Surprisingly extensive mixed phylogenetic and ecological signals among bacterial Operational Taxonomic Units

The lack of a consensus bacterial species concept greatly hampers our ability to understand and organize bacterial diversity. Operational taxonomic units (OTUs), which are clustered on the basis of DNA sequence identity alone, are the most commonly used microbial diversity unit. Although it is understood that OTUs can be phylogenetically incoherent, the degree and the extent of the phylogenetic inconsistency have not been explicitly studied. Here, we tested the phylogenetic signal of OTUs in a broad range of bacterial genera from various phyla. Strikingly, we found that very few OTUs were monophyletic, and many showed evidence of multiple independent origins. Using previously established bacterial habitats as benchmarks, we showed that OTUs frequently spanned multiple ecological habitats. We demonstrated that ecological heterogeneity within OTUs is caused by their phylogenetic inconsistency, and not merely due to ‘lumping’ of taxa resulting from using relaxed identity cut-offs. We argue that ecotypes, as described by the Stable Ecotype Model, are phylogenetically and ecologically more consistent than OTUs and therefore could serve as an alternative unit for bacterial diversity studies. In addition, we introduce QuickES, a new wrapper program for the Ecotype Simulation algorithm, which is capable of demarcating ecotypes in data sets with tens of thousands of sequences.     

Species matter: the role of competition in the assembly of congeneric bacteria

Interspecific competition is an important driver of community assembly in plants and animals, but phylogenetic evidence for interspecific competition in bacterial communities has been elusive. This could indicate that other processes such as habitat filtering or neutral processes are more important in bacterial community assembly. Alternatively, this could be a consequence of the lack of a consistent and meaningful species definition in bacteria. We hypothesize that competition in bacterial community assembly has gone undetected at least partly because overly broad measures of bacterial diversity units were used in previous studies. First, we tested our hypothesis in a simulation where we showed that how species are defined can dramatically affect whether phylogenetic overdispersion (a signal consistent with competitive exclusion) will be detected. Second, we demonstrated that using finer-scale Operational Taxonomic Units (OTUs) (with more stringent 16S rRNA sequence identity cutoffs or based on fast-evolving protein coding genes) in natural populations revealed previously undetected overdispersion. Finally, we argue that bacterial ecotypes, diversity units incorporating ecological and evolutionary theory, are superior to OTUs for the purpose of studying community assembly.     

Genomics, environmental genomics and the issue of microbial species

A microbial species concept is crucial for interpreting the variation detected by genomics and environmental genomics among cultivated microorganisms and within natural microbial populations. Comparative genomic analyses of prokaryotic species as they are presently described and named have led to the provocative idea that prokaryotes may not form species as we think about them for plants and animals. There are good reasons to doubt whether presently recognized prokaryotic species are truly species. To achieve a better understanding of microbial species, we believe it is necessary to (i) re-evaluate traditional approaches in light of evolutionary and ecological theory, (ii) consider that different microbial species may have evolved in different ways and (iii) integrate genomic, metagenomic and genome-wide expression approaches with ecological and evolutionary theory. Here, we outline how we are using genomic methods to (i) identify ecologically distinct populations (ecotypes) predicted by theory to be species-like fundamental units of microbial communities, and (ii) test their species-like character through in situ distribution and gene expression studies. By comparing metagenomic sequences obtained from well-studied hot spring cyanobacterial mats with genomic sequences of two cultivated cyanobacterial ecotypes, closely related to predominant native populations, we can conduct in situ population genetics studies that identify putative ecotypes and functional genes that determine the ecotypes' ecological distinctness. If individuals within microbial communities are found to be grouped into ecologically distinct, species-like populations, knowing about such populations should guide us to a better understanding of how genomic variation is linked to community function.     

Merging taxonomy with ecological population prediction in a case study of Vibrionaceae

We synthesized population structure data from three studies that assessed the fine-scale distribution of Vibrionaceae among temporally and spatially distinct environmental categories in coastal seawater and animals. All studies used a dynamic model (AdaptML) to identify phylogenetically cohesive and ecologically distinct bacterial populations and their predicted habitats without relying on a predefined genetic cutoff or relationships to previously named species. Across the three studies, populations were highly overlapping, displaying similar phylogenetic characteristics (identity and diversity), and were predominantly congruent with taxonomic Vibrio species previously characterized as genotypic clusters by multilocus sequence analysis (MLSA). The environmental fidelity of these populations appears high, with 9 out of 12 reproducibly associating with the same predicted (micro)habitats when similar environmental categories were sampled. Overall, this meta-analysis provides information on the habitat predictability and structure of previously described species, demonstrating that MLSA-based taxonomy can, at least in some cases, serve to approximate ecologically cohesive populations.     

The bacterial species definition in the genomic era

The bacterial species definition, despite its eminent practical significance for identification, diagnosis, quarantine and diversity surveys, remains a very difficult issue to advance. Genomics now offers novel insights into intra-species diversity and the potential for emergence of a more soundly based system. Although we share the excitement, we argue that it is premature for a universal change to the definition because current knowledge is based on too few phylogenetic groups and too few samples of natural populations. Our analysis of five important bacterial groups suggests, however, that more stringent standards for species may be justifiable when a solid understanding of gene content and ecological distinctiveness becomes available. Our analysis also reveals what is actually encompassed in a species according to the current standards, in terms of whole-genome sequence and gene-content diversity, and shows that this does not correspond to coherent clusters for the environmental Burkholderia and Shewanella genera examined. In contrast, the obligatory pathogens, which have a very restricted ecological niche, do exhibit clusters. Therefore, the idea of biologically meaningful clusters of diversity that applies to most eukaryotes may not be universally applicable in the microbial world, or if such clusters exist, they may be found at different levels of distinction.     

Cyanobacterial ecotypes in the microbial mat community of Mushroom Spring (Yellowstone National Park, Wyoming) as species-like units linking microbial community composition, structure and function

We have investigated microbial mats of alkaline siliceous hot springs in Yellowstone National Park as natural model communities to learn how microbial populations group into species-like fundamental units. Here, we bring together empirical patterns of the distribution of molecular variation in predominant mat cyanobacterial populations, theory-based modelling of how to demarcate phylogenetic clusters that correspond to ecological species and the dynamic patterns of the physical and chemical microenvironments these populations inhabit and towards which they have evolved adaptations. We show that putative ecotypes predicted by the theory-based model correspond well with distribution patterns, suggesting populations with distinct ecologies, as expected of ecological species. Further, we show that increased molecular resolution enhances our ability to detect ecotypes in this way, though yet higher molecular resolution is probably needed to detect all ecotypes in this microbial community.     

Evolutionary criteria outperform operational approaches in producing ecologically relevant fungal species inventories

Analyses of the structure and function of microbial communities are highly constrained by the diversity of organisms present within most environmental samples. A common approach is to rely almost entirely on DNA sequence data for estimates of microbial diversity, but to date there is no objective method of clustering sequences into groups that is grounded in evolutionary theory of what constitutes a biological lineage. The general mixed Yule-coalescent (GMYC) model uses a likelihood-based approach to distinguish population-level processes within lineages from processes associated with speciation and extinction, thus identifying a distinct point where extant lineages became independent. Using two independent surveys of DNA sequences associated with a group of ubiquitous plant-symbiotic fungi, we compared estimates of species richness derived using the GMYC model to those based on operational taxonomic units (OTUs) defined by fixed levels of sequence similarity. The model predicted lower species richness in these surveys than did traditional methods of sequence similarity. Here, we show for the first time that groups delineated by the GMYC model better explained variation in the distribution of fungi in relation to putative niche-based variables associated with host species identity, edaphic factors, and aspects of how the sampled ecosystems were managed. Our results suggest the coalescent-based GMYC model successfully groups environmental sequences of fungi into clusters that are ecologically more meaningful than more arbitrary approaches for estimating species richness.     

Cohesive molecular genetic data delineate species diversity in the dinoflagellate genus Symbiodinium

The diversity of symbiotic dinoflagellates ( Symbiodinium ) in pocilloporid corals originating from various reef habitats surrounding Heron Island, southern Great Barrier Reef, was examined by targeting ribosomal, mitochondrial, and chloroplast genes using six methods that analyse for sequence differences. The ability of each of 13 genetic analyses to characterize eight ecologically distinct Symbiodinium spp. was dependent on the level of conservation of the gene region targeted and the technique used. Other than differences in resolution, phylogenetic reconstructions using nuclear and organelle gene sequences were complementary and when combined produced a well-resolved phylogeny. Analysis of the ribosomal internal transcribed spacers using denaturing gradient gel electrophoresis fingerprinting in combination with sequencing of dominant bands provided a precise method for rapidly resolving and characterizing symbionts into ecologically and evolutionarily distinct units of diversity. Single-stranded conformation polymorphisms of the nuclear ribosomal large subunit (D1/D2 domain) identified the same number of ecologically distinct Symbiodinium spp., but profiles were less distinctive. The repetitive sequencing of bacterially cloned ITS2 polymerase chain reaction amplifications generated numerous sequence variants that clustered together according to the symbiont under analysis. The phylogenetic relationships between these clusters show how intragenomic variation in the ribosomal array diverges among closely related eukaryotic genomes. The strong correlation between phylogenetically independent lineages with different ecological and physiological attributes establishes a clear basis for assigning species designations to members of the genus Symbiodinium .     

Bighorn Sheep Genetic Structure in Wyoming Reflects Geography and Management

Aligning wildlife management boundaries with accurate biological units promotes effective conservation and management practices that reflect ecological and evolutionary processes. Neutral genetic markers allow for quantitative delineation of population structure without a priori assumptions or biases. In the United States, bighorn sheep (Ovis canadensis) are a charismatic component of Wyoming's biodiversity and a species that provides important viewing and hunting opportunities. Bighorn sheep abundances are relatively stable throughout Wyoming, and the species is managed by administrative units identified using expert knowledge, distribution and movement data, and geographic and administrative boundaries. We used a panel of 38 variable microsatellite loci and 512 base pairs of mitochondrial DNA sequence to identify the genetic structure throughout the state and in translocation source herds, quantify the extent of genetic diversity within each genetic cluster, and estimate the degree of gene flow among herds using blood and tissue samples collected 1989–2017. We identified genetic structure of Rocky Mountain bighorn sheep in the major mountain ranges of Wyoming, with strong support for ≥5 genetic clusters using microsatellite loci. These genetic clusters generally aligned with current management units, whereas mitochondrial data showed a more complex mosaic that was not geographically patterned. Genetic variation estimated from both markers was high within each herd and comparable among herds. The assignment of individuals reflected a combination of geographic isolation and translocation, which has been extensive. Our results provide a state-wide assessment of genetic diversity and structure that will enhance management by understanding the outcomes of translocation, identifying the source of unknown individuals, and parameterizing disease ecology models. © 2020 The Wildlife Society.     

A major shift in diversification rate helps explain macroevolutionary patterns in primate species diversity

Primates represent one of the most species rich, wide ranging, and ecologically diverse clades of mammals. What major macroevolutionary factors have driven their diversification and contributed to the modern distribution of primate species remains widely debated. We employed phylogenetic comparative methods to examine the role of clade age and evolutionary rate heterogeneity in the modern distribution of species diversity of Primates. Primate diversification has accelerated since its origin, with decreased extinction leading to a shift to even higher evolutionary rates in the most species rich family (Cercopithecidae). Older primate clades tended to be more diverse, however a shift in evolutionary rate was necessary to adequately explain the imbalance in species diversity. Species richness was also poorly explained by geographic distribution, especially once clade age and evolutionary rate shifts were accounted for, and may relate instead to other ecological factors. The global distribution of primate species diversity appears to have been strongly impacted by heterogeneity in evolutionary rates.     

The 3Cs provide a novel concept of bacterial species: messages from the genome as illustrated by Salmonella

A key issue troubling bacterial taxonomy and systematics is the lack of a biological species definition. Criteria to be used for defining bacterial species on genetic and biological bases should be able to reveal clear-cut boundaries among clusters of bacteria. To date, DNA–DNA re-association assays and ribosomal RNA sequence comparison have been useful in determining relative evolutionary distances among bacteria but the data are continuous and thus cannot define bacterial clusters as taxonomic units to be called species. Using Salmonella as models, we have looked for definite genetic and biologic uniqueness of clusters of bacteria. Based on our findings that each Salmonella lineage has a unique genome structure shared by strains of the same lineage but not overlapping with strains of other Salmonella lineages, we conclude that this is a result of genetic isolation following divergence of the bacteria. We propose that there should be genetic boundaries between different species of bacteria at the genomic level, which awaits further genomic information for validation.