O-antigen diversity and lateral transfer of the wbe region among Vibrio splendidus isolates

Potential metabolic and genetic interaction among viruses, methanogen and methanotrophic archaea, and their syntrophic partners

The metabolism of methane in anoxic ecosystems is mainly mediated by methanogens and methane-oxidizing archaea (MMA), key players in global carbon cycling. Viruses are vital in regulating their host fate and ecological function. However, our knowledge about the distribution and diversity of MMA viruses and their interactions with hosts is rather limited. Here, by searching metagenomes containing mcrA (the gene coding for the α-subunit of methyl-coenzyme M reductase) from a wide variety of environments, 140 viral operational taxonomic units (vOTUs) that potentially infect methanogens or methane-oxidizing archaea were retrieved. Four MMA vOTUs (three infecting the order Methanobacteriales and one infecting the order Methanococcales) were predicted to cross-domain infect sulfate-reducing bacteria. By facilitating assimilatory sulfur reduction, MMA viruses may increase the fitness of their hosts in sulfate-depleted anoxic ecosystems and benefit from synthesis of the sulfur-containing amino acid cysteine. Moreover, cell-cell aggregation promoted by MMA viruses may be beneficial for both the viruses and their hosts by improving infectivity and environmental stress resistance, respectively. Our results suggest a potential role of viruses in the ecological and environmental adaptation of methanogens and methane-oxidizing archaea.

Phosphonate production by marine microbes: Exploring new sources and potential function

SignificancePhosphonates are a class of phosphorus metabolites characterized by a highly stable C-P bond. Phosphonates accumulate to high concentrations in seawater, fuel a large fraction of marine methane production, and serve as a source of phosphorus to microbes inhabiting nutrient-limited regions of the oligotrophic ocean. Here, we show that 15% of all bacterioplankton in the surface ocean have genes phosphonate synthesis and that most belong to the abundant groups Prochlorococcus and SAR11. Genomic and chemical evidence suggests that phosphonates are incorporated into cell-surface phosphonoglycoproteins that may act to mitigate cell mortality by grazing and viral lysis. These results underscore the large global biogeochemical impact of relatively rare but highly expressed traits in numerically abundant groups of marine bacteria.

Rapid evolutionary turnover of mobile genetic elements drives bacterial resistance to phages

[Figure: see text].

Vibrio-bivalve interactions in health and disease

In the marine environment, bivalve mollusks constitute habitats for bacteria of the Vibrionaceae family. Vibrios belong to the microbiota of healthy oysters and mussels, which have the ability to concentrate bacteria in their tissues and body fluids, including the hemolymph. Remarkably, these important aquaculture species respond differently to infectious diseases. While oysters are the subject of recurrent mass mortalities at different life stages, mussels appear rather resistant to infections. Thus, Vibrio species are associated with the main diseases affecting the worldwide oyster production. Here, we review the current knowledge on Vibrio-bivalve interaction in oysters (Crassostrea sp.) and mussels (Mytilus sp.). We discuss the transient versus stable associations of vibrios with their bivalve hosts as well as technical issues limiting the monitoring of these bacteria in bivalve health and disease. Based on the current knowledge of oyster/mussel immunity and their interactions with Vibrio species pathogenic for oyster, we discuss how differences in immune effectors could contribute to the higher resistance of mussels to infections. Finally, we review the multiple strategies evolved by pathogenic vibrios to circumvent the potent immune defences of bivalves and how key virulence mechanisms could have been positively or negatively selected in the marine environment through interactions with predators.

Vibrio splendidus O-antigen structure: a trade-off between virulence to oysters and resistance to grazers

A major debate in evolutionary biology is whether virulence is maintained as an adaptive trait and/or evolves to non-virulence. In the environment, virulence traits of non-obligatory parasites are subjected to diverse selective pressures and trade-offs. Here, we focus on a population of Vibrio splendidus that displays moderate virulence for oysters. A MARTX (Multifunctional-autoprocessing repeats-in-toxin) and a type-six secretion system (T6SS) were found to be necessary for virulence toward oysters, while a region (wbe) involved in O-antigen synthesis is necessary for resistance to predation against amoebae. Gene inactivation within the wbe region had major consequences on the O-antigen structure, conferring lower immunogenicity, competitive advantage and increased virulence in oyster experimental infections. Therefore, O-antigen structures that favour resistance to environmental predators result in an increased activation of the oyster immune system and a reduced virulence in that host. These trade-offs likely contribute to maintaining O-antigen diversity in the marine environment by favouring genomic plasticity of the wbe region. The results of this study indicate an evolution of V. splendidus towards moderate virulence as a compromise between fitness in the oyster as a host, and resistance to its predators in the environment.

Genomic analysis unveils important aspects of population structure, virulence, and antimicrobial resistance in Klebsiella aerogenes

Klebsiella aerogenes is an important pathogen in healthcare-associated infections. Nevertheless, in comparison to other clinically important pathogens, K. aerogenes population structure, genetic diversity, and pathogenicity remain poorly understood. Here, we elucidate K. aerogenes clonal complexes (CCs) and genomic features associated with resistance and virulence. We present a detailed description of the population structure of K. aerogenes based on 97 publicly available genomes by using both multilocus sequence typing and single-nucleotide polymorphisms extracted from the core genome. We also assessed virulence and resistance profiles using Virulence Finder Database and Comprehensive Antibiotic Resistance Database, respectively. We show that K. aerogenes has an open pangenome and a large effective population size, which account for its high genomic diversity and support that negative selection prevents fixation of most deleterious alleles. The population is structured in at least 10 CCs, including two novel ones identified here, CC9 and CC10. The repertoires of resistance genes comprise a high number of antibiotic efflux proteins as well as narrow- and extended-spectrum β-lactamases. Regarding the population structure, we identified two clusters based on virulence profiles because of the presence of the toxin-encoding clb operon and the siderophore production genes, irp and ybt. Notably, CC3 comprises the majority of K. aerogenes isolates associated with hospital outbreaks, emphasizing the importance of constant monitoring of this pathogen. Collectively, our results may provide a foundation for the development of new therapeutic and surveillance strategies worldwide.

Ancestral gene acquisition as the key to virulence potential in environmental Vibrio populations

Diseases of marine animals caused by bacteria of the genus Vibrio are on the rise worldwide. Understanding the eco-evolutionary dynamics of these infectious agents is important for predicting and managing these diseases. Yet, compared to Vibrio infecting humans, knowledge of their role as animal pathogens is scarce. Here we ask how widespread is virulence among ecologically differentiated Vibrio populations, and what is the nature and frequency of virulence genes within these populations? We use a combination of population genomics and molecular genetics to assay hundreds of Vibrio strains for their virulence in the oyster Crassostrea gigas, a unique animal model that allows high-throughput infection assays. We show that within the diverse Splendidus clade, virulence represents an ancestral trait but has been lost from several populations. Two loci are necessary for virulence, the first being widely distributed across the Splendidus clade and consisting of an exported conserved protein (R5.7). The second is a MARTX toxin cluster, which only occurs within V. splendidus and is for the first time associated with virulence in marine invertebrates. Varying frequencies of both loci among populations indicate different selective pressures and alternative ecological roles, based on which we suggest strategies for epidemiological surveys.

Aquatic Pseudomonads Inhibit Oomycete Plant Pathogens of Glycine max

Seedling root rot of soybeans caused by the host-specific pathogen Phytophthora sojae, and a large number of Pythium species, is an economically important disease across the Midwest United States that negatively impacts soybean yields. Research on biocontrol strategies for crop pathogens has focused on compounds produced by microbes from soil, however, recent studies suggest that aquatic bacteria express distinct compounds that efficiently inhibit a wide range of pathogens. Based on these observations, we hypothesized that freshwater strains of pseudomonads might be producing novel antagonistic compounds that inhibit the growth of oomycetes. To test this prediction, we utilized a collection of 330 Pseudomonas strains isolated from soil and freshwater habitats, and determined their activity against a panel of five oomycetes: Phytophthora sojae, Pythium heterothalicum, Pythium irregulare, Pythium sylvaticum, and Pythium ultimum, all of which are pathogenic on soybeans. Among the bacterial strains, 118 exhibited antagonistic activity against at least one oomycete species, and 16 strains were inhibitory to all pathogens. Antagonistic activity toward oomycetes was significantly more common for aquatic isolates than for soil isolates. One water-derived strain, 06C 126, was predicted to express a siderophore and exhibited diverse antagonistic profiles when tested on nutrient rich and iron depleted media suggesting that more than one compound was produced that effectively inhibited oomycetes. These results support the concept that aquatic strains are an efficient source of compounds that inhibit pathogens. We outline a strategy to identify other strains that express unique compounds that may be useful biocontrol agents.

Vibrio tapetis from wrasse used for ectoparasite bio-control in salmon farming: phylogenetic analysis and serotyping

So-called 'cleaner fish', including various wrasse (Labridae) species, have become increasingly popular in Norwegian salmon farming in recent years for biocontrol of the salmon louse Lepeophtheirus salmonis. Cleaner fish mortalities in salmon farms are, however, often high. Various bacterial agents are frequently associated with episodes of increased cleaner fish mortality, and Vibrio tapetis is regularly cultured from diseased wrasse. In the present study, we investigated the genetic relationships among 54 V. tapetis isolates (34 from wrasse species) by multilocus sequence analysis (MLSA; rpoD, ftsZ, pyrH, rpoA and atpA). In the resulting phylogenetic tree, all wrasse isolates belonged to sub-clusters within V. tapetis subsp. tapetis. Slide agglutination testing further confirmed the complete dominance amongst these isolates of 4 O-antigen serotypes, designated here as V. tapetis subsp. tapetis serotypes O1, O3, O4 and O5, respectively. A pilot challenge trial using serotypes O3, O4 and O5 did not indicate high pathogenicity towards ballan wrasse Labrus bergylta, thus questioning the role of V. tapetis as a primary pathogen of this fish species.

Environmental Pseudomonads Inhibit Cystic Fibrosis Patient-Derived Pseudomonas aeruginosa

Pseudomonas aeruginosa is an opportunistic pathogen which is evolving resistance to many currently used antibiotics. While much research has been devoted to the roles of pathogenic P. aeruginosa in cystic fibrosis (CF) patients, less is known of its ecological properties. P. aeruginosa dominates the lungs during chronic infection in CF patients, yet its abundance in some environments is less than that of other diverse groups of pseudomonads. Here, we sought to determine if clinical isolates of P. aeruginosa are vulnerable to environmental pseudomonads that dominate soil and water habitats in one-to-one competitions which may provide a source of inhibitory factors. We isolated a total of 330 pseudomonads from diverse habitats of soil and freshwater ecosystems and competed these strains against one another to determine their capacity for antagonistic activity. Over 900 individual inhibitory events were observed. Extending the analysis to P. aeruginosa isolates revealed that clinical isolates, including ones with increased alginate production, were susceptible to competition by multiple environmental strains. We performed transposon mutagenesis on one isolate and identified an ∼14.8-kb locus involved in antagonistic activity. Only two other environmental isolates were observed to carry the locus, suggesting the presence of additional unique compounds or interactions among other isolates involved in outcompeting P. aeruginosa. This collection of strains represents a source of compounds that are active against multiple pathogenic strains. With the evolution of resistance of P. aeruginosa to currently used antibiotics, these environmental strains provide opportunities for novel compound discovery against drug-resistant clinical strains.

IMPORTANCE We demonstrate that clinical CF-derived isolates of P. aeruginosa are susceptible to competition in the presence of environmental pseudomonads. We observed that many diverse environmental strains exhibited varied antagonistic profiles against a panel of clinical P. aeruginosa isolates, suggesting the presence of distinct mechanisms of inhibition among these ecological strains. Understanding the properties of these antagonistic events offers the potential for discoveries of antimicrobial compounds or metabolic pathways important to the development of novel treatments for P. aeruginosa infections.

Lipopolysaccharide modification in Gram-negative bacteria during chronic infection

The Gram-negative bacterial lipopolysaccharide (LPS) is a major component of the outer membrane that plays a key role in host-pathogen interactions with the innate immune system. During infection, bacteria are exposed to a host environment that is typically dominated by inflammatory cells and soluble factors, including antibiotics, which provide cues about regulation of gene expression. Bacterial adaptive changes including modulation of LPS synthesis and structure are a conserved theme in infections, irrespective of the type or bacteria or the site of infection. In general, these changes result in immune system evasion, persisting inflammation and increased antimicrobial resistance. Here, we review the modifications of LPS structure and biosynthetic pathways that occur upon adaptation of model opportunistic pathogens (Pseudomonas aeruginosa, Burkholderia cepacia complex bacteria, Helicobacter pylori and Salmonella enterica) to chronic infection in respiratory and gastrointestinal sites. We also discuss the molecular mechanisms of these variations and their role in the host-pathogen interaction.

The O-antigen mediates differential survival of Salmonella against communities of natural predators

Antigenically distinct members of bacterial species can be differentially distributed in the environment. Predators known to consume antigenically distinct prey with different efficiencies are also differentially distributed. Here we show that antigenically distinct, but otherwise isogenic and physiologically indistinct, strains of Salmonella enterica show differential survival in natural soil, sediment and intestinal environments, where they would face a community of predators. Decline in overall cell numbers is attenuated by factors that inhibit the action of predators, including heat and antiprotozoal and antihelminthic drugs. Moreover, the fitness of strains facing these predators - calculated by comparing survival with and without treatments attenuating predator activity - varies between environments. These results suggest that relative survival in natural environments is arbitrated by communities of natural predators whose feeding preferences, if not species composition, vary between environments. These data support the hypothesis that survival against natural predators may drive the differential distribution of bacteria among microenvironments.

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.

Phylogenetic analysis and serotyping of Vibrio splendidus-related bacteria isolated from salmon farm cleaner fish

Cleaner fish, i.e. various wrasse (Labridae) species and lumpsucker Cyclopterus lumpus, are to an increasing extent used for biocontrol of the salmon louse Lepeophtheirus salmonis in European salmon farming. Although efficient de-licers, cleaner fish mortality levels in salmon farms are often high. Bacterial infections are common, and Vibrio splendidus-related strains are frequently identified during diagnostic investigations. The population structure of 112 V. splendidus-related isolates, derived primarily from wrasse species, was investigated by means of multilocus sequence analysis using 5 housekeeping genes (rpoD, ftsZ, pyrH, rpoA and atpA). Most isolates were found to be closely related to the V. splendidus type strain, yet displayed extensive genetic microdiversity. Slide agglutination testing using polyclonal rabbit antisera further indicated O-antigen variability. Intra-outbreak genetic and antigenic diversity suggests direct infection from seawater, rather than fish-to-fish transmission, as the main route of infection. The variable nature of isolates involved complicates qualified selection of representative candidate strains, e.g. for infection and vaccine trials.

Gene-specific selective sweeps in bacteria and archaea caused by negative frequency-dependent selection

BACKGROUND

Fixation of beneficial genes in bacteria and archaea (collectively, prokaryotes) is often believed to erase pre-existing genomic diversity through the hitchhiking effect, a phenomenon known as genome-wide selective sweep. Recent studies, however, indicate that beneficial genes spread through a prokaryotic population via recombination without causing genome-wide selective sweeps. These gene-specific selective sweeps seem to be at odds with the existing estimates of recombination rates in prokaryotes, which appear far too low to explain such phenomena.

RESULTS

We use mathematical modeling to investigate potential solutions to this apparent paradox. Most microbes in nature evolve in heterogeneous, dynamic communities, in which ecological interactions can substantially impact evolution. Here, we focus on the effect of negative frequency-dependent selection (NFDS) such as caused by viral predation (kill-the-winner dynamics). The NFDS maintains multiple genotypes within a population, so that a gene beneficial to every individual would have to spread via recombination, hence a gene-specific selective sweep. However, gene loci affected by NFDS often are located in variable regions of microbial genomes that contain genes involved in the mobility of selfish genetic elements, such as integrases or transposases. Thus, the NFDS-affected loci are likely to experience elevated rates of recombination compared with the other loci. Consequently, these loci might be effectively unlinked from the rest of the genome, so that NFDS would be unable to prevent genome-wide selective sweeps. To address this problem, we analyzed population genetic models of selective sweeps in prokaryotes under NFDS. The results indicate that NFDS can cause gene-specific selective sweeps despite the effect of locally elevated recombination rates, provided NFDS affects more than one locus and the basal rate of recombination is sufficiently low. Although these conditions might seem to contradict the intuition that gene-specific selective sweeps require high recombination rates, they actually decrease the effective rate of recombination at loci affected by NFDS relative to the per-locus basal level, so that NFDS can cause gene-specific selective sweeps.

CONCLUSION

Because many free-living prokaryotes are likely to evolve under NFDS caused by ubiquitous viruses, gene-specific selective sweeps driven by NFDS are expected to be a major, general phenomenon in prokaryotic populations.

Horizontal transfers and gene losses in the phospholipid pathway of bartonella reveal clues about early ecological niches

Bartonellae are mammalian pathogens vectored by blood-feeding arthropods. Although of increasing medical importance, little is known about their ecological past, and host associations are underexplored. Previous studies suggest an influence of horizontal gene transfers in ecological niche colonization by acquisition of host pathogenicity genes. We here expand these analyses to metabolic pathways of 28 Bartonella genomes, and experimentally explore the distribution of bartonellae in 21 species of blood-feeding arthropods. Across genomes, repeated gene losses and horizontal gains in the phospholipid pathway were found. The evolutionary timing of these patterns suggests functional consequences likely leading to an early intracellular lifestyle for stem bartonellae. Comparative phylogenomic analyses discover three independent lineage-specific reacquisitions of a core metabolic gene—NAD(P)H-dependent glycerol-3-phosphate dehydrogenase (gpsA)—from Gammaproteobacteria and Epsilonproteobacteria. Transferred genes are significantly closely related to invertebrate Arsenophonus-, and Serratia-like endosymbionts, and mammalian Helicobacter-like pathogens, supporting a cellular association with arthropods and mammals at the base of extant Bartonella spp. Our studies suggest that the horizontal reacquisitions had a key impact on bartonellae lineage specific ecological and functional evolution.

Explaining microbial genomic diversity in light of evolutionary ecology

Comparisons of closely related microorganisms have shown that individual genomes can be highly diverse in terms of gene content. In this Review, we discuss several studies showing that much of this variation is associated with social and ecological interactions, which have an important role in the population biology of wild populations of bacteria and archaea. These interactions create frequency-dependent selective pressures that can either stabilize gene frequencies at intermediate levels in populations or promote fast gene turnover, which presents as low gene frequencies in genome surveys. Thus, interpretation of gene-content diversity requires the delineation of populations according to cohesive gene flow and ecology, as micro-evolutionary changes arise in response to local selection pressures and population dynamics.

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.

Comparative genomics of the classical Bordetella subspecies: the evolution and exchange of virulence-associated diversity amongst closely related pathogens

Background

The classical Bordetella subspecies are phylogenetically closely related, yet differ in some of the most interesting and important characteristics of pathogens, such as host range, virulence and persistence. The compelling picture from previous comparisons of the three sequenced genomes was of genome degradation, with substantial loss of genome content (up to 24%) associated with adaptation to humans.

Results

For a more comprehensive picture of lineage evolution, we employed comparative genomic and phylogenomic analyses using seven additional diverse, newly sequenced Bordetella isolates. Genome-wide single nucleotide polymorphism (SNP) analysis supports a reevaluation of the phylogenetic relationships between the classical Bordetella subspecies, and suggests a closer link between ovine and human B. parapertussis lineages than has been previously proposed. Comparative analyses of genome content revealed that only 50% of the pan-genome is conserved in all strains, reflecting substantial diversity of genome content in these closely related pathogens that may relate to their different host ranges, virulence and persistence characteristics. Strikingly, these analyses suggest possible horizontal gene transfer (HGT) events in multiple loci encoding virulence factors, including O-antigen and pertussis toxin (Ptx). Segments of the pertussis toxin locus (ptx) and its secretion system locus (ptl) appear to have been acquired by the classical Bordetella subspecies and are divergent in different lineages, suggesting functional divergence in the classical Bordetellae.

Conclusions

Together, these observations, especially in key virulence factors, reveal that multiple mechanisms, such as point mutations, gain or loss of genes, as well as HGTs, contribute to the substantial phenotypic diversity of these versatile subspecies in various hosts.

Unique core genomes of the bacterial family vibrionaceae: insights into niche adaptation and speciation

Background

The criteria for defining bacterial species and even the concept of bacterial species itself are under debate, and the discussion is apparently intensifying as more genome sequence data is becoming available. However, it is still unclear how the new advances in genomics should be used most efficiently to address this question. In this study we identify genes that are common to any group of genomes in our dataset, to determine whether genes specific to a particular taxon exist and to investigate their potential role in adaptation of bacteria to their specific niche. These genes were named unique core genes. Additionally, we investigate the existence and importance of unique core genes that are found in isolates of phylogenetically non-coherent groups. These groups of isolates, that share a genetic feature without sharing a closest common ancestor, are termed genophyletic groups.

Results

The bacterial family Vibrionaceae was used as the model, and we compiled and compared genome sequences of 64 different isolates. Using the software orthoMCL we determined clusters of homologous genes among the investigated genome sequences. We used multilocus sequence analysis to build a host phylogeny and mapped the numbers of unique core genes of all distinct groups of isolates onto the tree. The results show that unique core genes are more likely to be found in monophyletic groups of isolates. Genophyletic groups of isolates, in contrast, are less common especially for large groups of isolate. The subsequent annotation of unique core genes that are present in genophyletic groups indicate a high degree of horizontally transferred genes. Finally, the annotation of the unique core genes of Vibrio cholerae revealed genes involved in aerotaxis and biosynthesis of the iron-chelator vibriobactin.

Conclusion

The presented work indicates that genes specific for any taxon inside the bacterial family Vibrionaceae exist. These unique core genes encode conserved metabolic functions that can shed light on the adaptation of a species to its ecological niche. Additionally, our study suggests that unique core genes can be used to aid classification of bacteria and contribute to a bacterial species definition on a genomic level. Furthermore, these genes may be of importance in clinical diagnostics and drug development.

Population genomics of early events in the ecological differentiation of bacteria

Genetic exchange is common among bacteria, but its effect on population diversity during ecological differentiation remains controversial. A fundamental question is whether advantageous mutations lead to selection of clonal genomes or, as in sexual eukaryotes, sweep through populations on their own. Here, we show that in two recently diverged populations of ocean bacteria, ecological differentiation has occurred akin to a sexual mechanism: A few genome regions have swept through subpopulations in a habitat-specific manner, accompanied by gradual separation of gene pools as evidenced by increased habitat specificity of the most recent recombinations. These findings reconcile previous, seemingly contradictory empirical observations of the genetic structure of bacterial populations and point to a more unified process of differentiation in bacteria and sexual eukaryotes than previously thought.

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.