Calibrating bacterial evolution

Multidimensional insights into the biodiversity of Streptomyces in soils of China: a pilot study

Streptomyces, a diverse group of filamentous bacteria found predominantly in soil, play a crucial role in nutrient cycling and produce many valuable secondary metabolites for the pharmaceutical industry. In this pilot study, we collected 19 soil samples from 14 provinces in China to preliminarily investigate the biodiversity and genetic structure of Streptomyces in soils of China from different dimensions, using recently developed cost-efficient amplicon and whole-genome library preparation methods. Amplicon analysis showed that Actinobacteria were among the most abundant bacteria, with 0.3% of amplicon sequence variants (ASVs) belonging to Streptomyces. Meanwhile, we successfully isolated 136 Streptomyces natural strains and assembled their genomes, including 26 previously unreported species, underscoring the need for further exploration of soil Streptomyces in China. Population genetics analysis revealed that homologous recombination may primarily drive the extensive genetic diversity observed in Streptomyces, as well as a complex population structure. Complementing this, pan-genome analysis shed light on gene diversity within Streptomyces and led to the discovery of rare genes, further emphasizing the vast genetic diversity of this genus. Additionally, multiple metabolic gene clusters were found in these Streptomyces strains, as well as some potentially unique or uncommon ones were found. These findings not only highlight the biological and metabolic diversity of Streptomyces but also provide a technical framework for future studies on the global biodiversity and evolution of this genus.

IMPORTANCE

Streptomyces, a prominent group of Actinobacteria, holds significant importance in ecosystems and biotechnology due to their diverse array of metabolic products. However, research on the biodiversity of soil Streptomyces across extensive geographical scales in China has been limited, and their genetic diversity has rarely been evaluated using modern population genetics principles. This pilot study successfully addresses these gaps by conducting a preliminary exploration on the biodiversity of Streptomyces in Chinese soils from multiple perspectives, providing valuable insights for a deeper understanding of their biodiversity and a novel technical framework for future large-scale explorations of its diversity.

A Short-Term View of Protein Sequence Evolution from Salmonella

Much of the study of protein sequence evolution is based on sequence changes inferred to have occurred in nature. The sequences compared for this purpose are usually sufficiently distant that purifying selection has had nearly its full effect and most of the changes inferred have been exposed to a variety of conditions. Here, I make use of large numbers of Salmonella genome sequences to study changes known to be of very recent origin because they are inferred from comparison of very closely related sequences. The effects of purifying selection are weak yet discernible on this short timescale: the ratio of nonsynonymous to synonymous changes is smaller than expected under selective neutrality, but only slightly so. Essential genes have lower rates of nonsynonymous change, as they do on a longer timescale, but much more of this association remains after controlling for expression level. Positive selection for nonsynonymous change is inferred for 151 genes. For nearly half of these, this is attributable to selection for loss of function. Other forms of positive selection inferred include selection for amino acid changes that make enzymes less sensitive to antibiotics and selection for activating changes to proteins involved in transcriptional regulation. Positively selected variants of many genes are likely favored only under unusual conditions and disfavored in the long term, making detection of the positive selection with more distant comparisons difficult or impossible. The short-term view provided by close comparisons complements the long-term view obtained from more distant comparisons such as those between species.

Geographically widespread and novel hemotropic mycoplasmas and bartonellae in Mexican free-tailed bats and sympatric North American bat species

Bacterial pathogens remain poorly characterized in bats, especially in North America. We describe novel (and in some cases panmictic) hemoplasmas (10.1% positivity) and bartonellae (25.6% positivity) across three colonies of Mexican free-tailed bats (Tadarida brasiliensis), a partially migratory species that can seasonally travel hundreds of kilometers. Molecular analyses identified three novel Candidatus hemoplasma species most similar to another novel Candidatus species in Neotropical molossid bats. We also detected novel hemoplasmas in sympatric cave myotis (Myotis velifer) and pallid bats (Antrozous pallidus), with sequences in the latter 96.5% related to Candidatus Mycoplasma haematohominis. We identified nine Bartonella genogroups, including those in cave myotis with 96.1% similarity to Candidatus Bartonella mayotimonensis. We also detected Bartonella rochalimae in migratory Mexican free-tailed bats, representing the first report of this human pathogen in the Chiroptera. Monthly sampling of migratory Mexican free-tailed bats during their North American occupancy period also revealed significant seasonality in infection for both bacterial pathogens, with prevalence increasing following spring migration, peaking in the maternity season, and declining into fall migration. The substantial diversity and seasonality of hemoplasmas and bartonellae observed here suggest that additional longitudinal, genomic, and immunological studies in bats are warranted to inform One Health approaches.

IMPORTANCE

Bats have been intensively sampled for viruses but remain mostly understudied for bacterial pathogens. However, bacterial pathogens can have significant impacts on both human health and bat morbidity and even mortality. Hemoplasmas and bartonellae are facultative intracellular bacteria of special interest in bats, given their high prevalence and substantial genetic diversity. Surveys have also supported plausible zoonotic transmission of these bacteria from bats to humans, including Candidatus Mycoplasma haematohominis and Candidatus Bartonella mayotimonensis. Greater characterization of these bacteria across global bat diversity (over 1,480 species) is therefore warranted to inform infection risks for both bats and humans, although little surveillance has thus far been conducted in North American bats. We here describe novel (and in some cases panmictic) hemoplasmas and bartonellae across three colonies of Mexican free-tailed bats and sympatric bat species. We find high genetic diversity and seasonality of these pathogens, including lineages closely related to human pathogens, such as Bartonella rochalimae.

Microbiome divergence of marine gastropod species separated by the Isthmus of Panama

The rise of the Isthmus of Panama separated the populations of many marine organisms, which then diverged into new geminate sister species currently living in the Eastern Pacific Ocean and the Caribbean Sea. However, we know very little about how such evolutionary divergences of host species have shaped the compositions of their microbiomes. Here, we compared the microbiomes of whole-body and shell-surface samples of geminate species of marine gastropods in the genera Cerithium and Cerithideopsis to those of congeneric outgroups. Our results suggest that the effects of ~3 million years of separation and isolation on microbiome composition varied among host genera and between sample types within the same hosts. In the whole-body samples, microbiome compositions of geminate species pairs tended to be similar, likely due to host filtering, although the strength of this relationship varied among the two groups and across similarity metrics. Shell-surface microbiomes show contrasting patterns, with co-divergence between the host taxa and a small number of microbial clades evident in Cerithideopsis but not Cerithium. These results suggest that (i) isolation of host populations after the rise of the Isthmus of Panama affected microbiomes of geminate hosts in a complex and host-specific manner, and (ii) host-associated microbial taxa respond differently to vicariance events than the hosts themselves.IMPORTANCEWhile considerable work has been done on evolutionary divergences of marine species in response to the rise of the Isthmus of Panama, which separated two previously connected oceans, how this event shaped the microbiomes of these marine hosts remains poorly known. Using whole-body and shell-surface microbiomes of closely related gastropod species from opposite sides of the Isthmus, we show that divergences of microbial taxa after the formation of the Isthmus are often not concordant with those of their gastropod hosts. Our results show that evolutionary responses of marine gastropod-associated microbiomes to major environmental perturbations are complex and are shaped more by local environments than host evolutionary history.

Geographically widespread and novel hemotropic mycoplasmas and bartonellae in Mexican free-tailed bats and sympatric North American bat species

Bacterial pathogens remain poorly characterized in bats, especially in North America. We describe novel (and in some cases panmictic) hemoplasmas (10.5% positivity) and bartonellae (25.5% positivity) across three colonies of Mexican free-tailed bats (Tadarida brasiliensis), a partially migratory species that can seasonally travel hundreds of kilometers. Molecular analyses identified three novel Candidatus hemoplasma species most similar to another novel Candidatus species in Neotropical molossid bats. We also detected novel hemoplasmas in sympatric cave myotis (Myotis velifer) and pallid bats (Antrozous pallidus), with sequences in the latter 96.5% related to C. Mycoplasma haemohominis. We identified nine Bartonella genogroups, including those in cave myotis with 96.7% similarity to C. Bartonella mayotimonensis. We also detected Bartonella rochalimae in migratory Mexican free-tailed bats, representing the first report of this human pathogen in the Chiroptera. The seasonality and diversity of these bacteria observed here suggest that additional longitudinal, genomic, and immunological studies in bats are warranted.

Host specialization and spatial divergence of bacteria associated with Peltigera lichens promote landscape gamma diversity

BACKGROUND

Lichens are micro-ecosystems relying on diverse microorganisms for nutrient cycling, environmental adaptation, and structural support. We investigated the spatial-scale dependency of factors shaping the ecological processes that govern lichen-associated bacteria. We hypothesize that lichens function as island-like habitats hosting divergent microbiomes and promoting landscape gamma-diversity. Three microenvironments -thalli, substrates, and neighboring soils- were sampled from four geographically overlapping species of Peltigera cyanolichens, spanning three bioclimatic zones in the Chilean Patagonia, to determine how bacterial diversity, assembly processes, ecological drivers, interaction patterns, and niche breadth vary among Peltigera microenvironments on a broad geographical scale.

RESULTS

The hosts' phylogeny, especially that of the cyanobiont, alongside climate as a secondary factor, impose a strong ecological filtering of bacterial communities within Peltigera thalli. This results in deterministically assembled, low diverse, and phylogenetically convergent yet structurally divergent bacterial communities. Host evolutionary and geographic distances accentuate the divergence in bacterial community composition of Peltigera thalli. Compared to soil and substrate, Peltigera thalli harbor specialized and locally adapted bacterial taxa, conforming sparse and weak ecological networks.

CONCLUSIONS

The findings suggest that Petigera thalli create fragmented habitats that foster landscape bacterial gamma-diversity. This underscores the importance of preserving lichens for maintaining a potential reservoir of specialized bacteria.

Discovery of disease-adapted bacterial lineages in inflammatory bowel diseases

Gut bacteria are implicated in inflammatory bowel disease (IBD), but the strains driving these associations are unknown. Large-scale studies of microbiome evolution could reveal the imprint of disease on gut bacteria, thus pinpointing the strains and genes that may underlie inflammation. Here, we use stool metagenomes of thousands of IBD patients and healthy controls to reconstruct 140,000 strain genotypes, revealing hundreds of lineages enriched in IBD. We demonstrate that these strains are ancient, taxonomically diverse, and ubiquitous in humans. Moreover, disease-associated strains outcompete their healthy counterparts during inflammation, implying long-term adaptation to disease. Strain genetic differences map onto known axes of inflammation, including oxidative stress, nutrient biosynthesis, and immune evasion. Lastly, the loss of health-associated strains of Eggerthella lenta was predictive of fecal calprotectin, a biomarker of disease severity. Our work identifies reservoirs of strain diversity that may impact inflammatory disease and can be extended to other microbiome-associated diseases.

Mutation Rate and Effective Population Size of the Model Cooperative Bacterium Myxococcus xanthus

Intrinsic rates of genetic mutation have diverged greatly across taxa and exhibit statistical associations with several other parameters and features. These include effective population size (Ne), genome size, and gametic multicellularity, with the latter being associated with both increased mutation rates and decreased effective population sizes. However, data sufficient to test for possible relationships between microbial multicellularity and mutation rate (µ) are lacking. Here, we report estimates of two key population-genetic parameters, Ne and µ, for Myxococcus xanthus, a bacterial model organism for the study of aggregative multicellular development, predation, and social swarming. To estimate µ, we conducted an ∼400-day mutation accumulation experiment with 46 lineages subjected to regular single colony bottlenecks prior to clonal regrowth. Upon conclusion, we sequenced one clonal-isolate genome per lineage. Given collective evolution for 85,323 generations across all lines, we calculate a per base-pair mutation rate of ∼5.5 × 10-10 per site per generation, one of the highest mutation rates among free-living eubacteria. Given our estimate of µ, we derived Ne at ∼107 from neutral diversity at four-fold degenerate sites across two dozen M. xanthus natural isolates. This estimate is below average for eubacteria and strengthens an already clear negative correlation between µ and Ne in prokaryotes. The higher and lower than average mutation rate and Ne for M. xanthus, respectively, amplify the question of whether any features of its multicellular life cycle-such as group-size reduction during fruiting-body development-or its highly structured spatial distribution have significantly influenced how these parameters have evolved.

Unveiling the co-phylogeny signal between plunderfish Harpagifer spp. and their gut microbiomes across the Southern Ocean

Understanding the factors that sculpt fish gut microbiome is challenging, especially in natural populations characterized by high environmental and host genomic complexity. However, closely related hosts are valuable models for deciphering the contribution of host evolutionary history to microbiome assembly, through the underscoring of phylosymbiosis and co-phylogeny patterns. Here, we propose that the recent diversification of several Harpagifer species across the Southern Ocean would allow the detection of robust phylogenetic congruence between the host and its microbiome. We characterized the gut mucosa microbiome of 77 individuals from four field-collected species of the plunderfish Harpagifer (Teleostei, Notothenioidei), distributed across three biogeographic regions of the Southern Ocean. We found that seawater physicochemical properties, host phylogeny, and geography collectively explained 35% of the variation in bacterial community composition in Harpagifer gut mucosa. The core microbiome of Harpagifer spp. gut mucosa was characterized by a low diversity, mostly driven by selective processes, and dominated by a single Aliivibrio Operational Taxonomic Unit (OTU) detected in more than 80% of the individuals. Nearly half of the core microbiome taxa, including Aliivibrio, harbored co-phylogeny signal at microdiversity resolution with host phylogeny, indicating an intimate symbiotic relationship and a shared evolutionary history with Harpagifer. The clear phylosymbiosis and co-phylogeny signals underscore the relevance of the Harpagifer model in understanding the role of fish evolutionary history in shaping the gut microbiome assembly. We propose that the recent diversification of Harpagifer may have led to the diversification of Aliivibrio, exhibiting patterns that mirror the host phylogeny.

IMPORTANCE

Although challenging to detect in wild populations, phylogenetic congruence between marine fish and its microbiome is critical, as it highlights intimate associations between hosts and ecologically relevant microbial symbionts. Our study leverages a natural system of closely related fish species in the Southern Ocean to unveil new insights into the contribution of host evolutionary trajectory on gut microbiome assembly, an underappreciated driver of the global marine fish holobiont. Notably, we unveiled striking evidence of co-diversification between Harpagifer and its microbiome, demonstrating both phylosymbiosis of gut bacterial communities and co-phylogeny of some specific bacterial symbionts, mirroring the host diversification patterns. Given Harpagifer's significance as a trophic resource in coastal areas and its vulnerability to climatic and anthropic pressures, understanding the potential evolutionary interdependence between the hosts and its microbiome provides valuable microbial candidates for future monitoring, as they may play a pivotal role in host species acclimatization to a rapidly changing environment.

Modeling the mosaic structure of bacterial genomes to infer their evolutionary history

The chronology and phylogeny of bacterial evolution are difficult to reconstruct due to a scarce fossil record. The analysis of bacterial genomes remains challenging because of large sequence divergence, the plasticity of bacterial genomes due to frequent gene loss, horizontal gene transfer, and differences in selective pressure from one locus to another. Therefore, taking advantage of the rich and rapidly accumulating genomic data requires accurate modeling of genome evolution. An important technical consideration is that loci with high effective mutation rates may diverge beyond the detection limit of the alignment algorithms used, biasing the genome-wide divergence estimates toward smaller divergences. In this article, we propose a novel method to gain insight into bacterial evolution based on statistical properties of genome comparisons. We find that the length distribution of sequence matches is shaped by the effective mutation rates of different loci, by the horizontal transfers, and by the aligner sensitivity. Based on these inputs, we build a model and show that it accounts for the empirically observed distributions, taking the Enterobacteriaceae family as an example. Our method allows to distinguish segments of vertical and horizontal origins and to estimate the time divergence and exchange rate between any pair of taxa from genome-wide alignments. Based on the estimated time divergences, we construct a time-calibrated phylogenetic tree to demonstrate the accuracy of the method.

Within-Host Evolution of Bacterial Pathogens in Acute and Chronic Infection

Bacterial pathogens undergo remarkable adaptive change in response to the selective forces they encounter during host colonization and infection. Studies performed over the past few decades have demonstrated that many general evolutionary processes can be discerned during the course of host adaptation, including genetic diversification of lineages, clonal succession events, convergent evolution, and balanced fitness trade-offs. In some cases, elevated mutation rates resulting from mismatch repair or proofreading deficiencies accelerate evolution, and active mobile genetic elements or phages may facilitate genome plasticity. The host immune response provides another critical component of the fitness landscapes guiding adaptation, and selection operating on pathogens at this level may lead to immune evasion and the establishment of chronic infection. This review summarizes recent advances in this field, with a special focus on different forms of bacterial genome plasticity in the context of infection, and considers clinical consequences of adaptive changes for the host.

Functional host-specific adaptation of the intestinal microbiome in hominids

Fine-scale knowledge of the changes in composition and function of the human gut microbiome compared that of our closest relatives is critical for understanding the evolutionary processes underlying its developmental trajectory. To infer taxonomic and functional changes in the gut microbiome across hominids at different timescales, we perform high-resolution metagenomic-based analyzes of the fecal microbiome from over two hundred samples including diverse human populations, as well as wild-living chimpanzees, bonobos, and gorillas. We find human-associated taxa depleted within non-human apes and patterns of host-specific gut microbiota, suggesting the widespread acquisition of novel microbial clades along the evolutionary divergence of hosts. In contrast, we reveal multiple lines of evidence for a pervasive loss of diversity in human populations in correlation with a high Human Development Index, including evolutionarily conserved clades. Similarly, patterns of co-phylogeny between microbes and hosts are found to be disrupted in humans. Together with identifying individual microbial taxa and functional adaptations that correlate to host phylogeny, these findings offer insights into specific candidates playing a role in the diverging trajectories of the gut microbiome of hominids. We find that repeated horizontal gene transfer and gene loss, as well as the adaptation to transient microaerobic conditions appear to have played a role in the evolution of the human gut microbiome.

Competition and evolutionary selection among core regulatory motifs in gene expression control

Gene products that are beneficial in one environment may become burdensome in another, prompting the emergence of diverse regulatory schemes that carry their own bioenergetic cost. By ensuring that regulators are only expressed when needed, we demonstrate that autoregulation generally offers an advantage in an environment combining mutation and time-varying selection. Whether positive or negative feedback emerges as dominant depends primarily on the demand for the target gene product, typically to ensure that the detrimental impact of inevitable mutations is minimized. While self-repression of the regulator curbs the spread of these loss-of-function mutations, self-activation instead facilitates their propagation. By analyzing the transcription network of multiple model organisms, we reveal that reduced bioenergetic cost may contribute to the preferential selection of autoregulation among transcription factors. Our results not only uncover how seemingly equivalent regulatory motifs have fundamentally different impact on population structure, growth dynamics, and evolutionary outcomes, but they can also be leveraged to promote the design of evolutionarily robust synthetic gene circuits.

Comparing different computational approaches for detecting long-term vertical transmission in host-associated microbiota

Long-term vertical transmissions of gut bacteria are thought to be frequent and functionally important in mammals. Several phylogenetic-based approaches have been proposed to detect, among species-rich microbiota, the bacteria that have been vertically transmitted during a host clade radiation. Applied to mammal microbiota, these methods have sometimes led to conflicting results; in addition, how they cope with the slow evolution of markers typically used to characterize bacterial microbiota remains unclear. Here, we use simulations to test the statistical performances of two widely-used global-fit approaches (ParaFit and PACo) and two event-based approaches (ALE and HOME). We find that these approaches have different strengths and weaknesses depending on the amount of variation in the bacterial DNA sequences and are therefore complementary. In particular, we show that ALE performs better when there is a lot of variation in the bacterial DNA sequences, whereas HOME performs better when there is not. Global-fit approaches (ParaFit and PACo) have higher type I error rates (false positives) but have the advantage to be very fast to run. We apply these methods to the gut microbiota of primates and our results suggest that only a small fraction of their gut bacteria is vertically transmitted.

Expanded diversity of novel hemoplasmas in rare and undersampled Neotropical bats

Hemotropic mycoplasmas are emerging as a model system for studying bacterial pathogens in bats, but taxonomic coverage of sampled host species remains biased. We leveraged a long-term field study in Belize to uncover novel hemoplasma diversity in bats by analyzing 80 samples from 19 species, most of which are infrequently encountered. PCR targeting the partial 16S rRNA gene found 41% of bats positive for hemoplasmas. Phylogenetic analyses found two novel host shifts of hemoplasmas, four entirely new hemoplasma genotypes, and the first hemoplasma detections in four bat species. One of these novel hemoplasmas (from Neoeptesicus furinalis) shared 97.6% identity in the partial 16S rRNA gene to a human hemoplasma (Candidatus Mycoplasma haemohominis). Additional analysis of the partial 23S rRNA gene allowed us to also designate two novel hemoplasma species, in Myotis elegans and Phyllostomus discolor, with the proposed names Candidatus Mycoplasma haematomyotis sp. nov. and Candidatus Mycoplasma haematophyllostomi sp. nov., respectively. Our analyses show that additional hemoplasma diversity in bats can be uncovered by targeting rare or undersampled host species.

Microbial symbionts are shared between ants and their associated beetles

The transmission of microbial symbionts across animal species could strongly affect their biology and evolution, but our understanding of transmission patterns and dynamics is limited. Army ants (Formicidae: Dorylinae) and their hundreds of closely associated insect guest species (myrmecophiles) can provide unique insights into interspecific microbial symbiont sharing. Here, we compared the microbiota of workers and larvae of the army ant Eciton burchellii with those of 13 myrmecophile beetle species using 16S rRNA amplicon sequencing. We found that the previously characterized specialized bacterial symbionts of army ant workers were largely absent from ant larvae and myrmecophiles, whose microbial communities were usually dominated by Rickettsia, Wolbachia, Rickettsiella and/or Weissella. Strikingly, different species of myrmecophiles and ant larvae often shared identical 16S rRNA genotypes of these common bacteria. Protein-coding gene sequences confirmed the close relationship of Weissella strains colonizing army ant larvae, some workers and several myrmecophile species. Unexpectedly, these strains were also similar to strains infecting dissimilar animals inhabiting very different habitats: trout and whales. Together, our data show that closely interacting species can share much of their microbiota, and some versatile microbial species can inhabit and possibly transmit across a diverse range of hosts and environments.

Exposing the small protein load of bacterial life

The ever-growing repertoire of genomic techniques continues to expand our understanding of the true diversity and richness of prokaryotic genomes. Riboproteogenomics laid the foundation for dynamic studies of previously overlooked genomic elements. Most strikingly, bacterial genomes were revealed to harbor robust repertoires of small open reading frames (sORFs) encoding a diverse and broadly expressed range of small proteins, or sORF-encoded polypeptides (SEPs). In recent years, continuous efforts led to great improvements in the annotation and characterization of such proteins, yet many challenges remain to fully comprehend the pervasive nature of small proteins and their impact on bacterial biology. In this work, we review the recent developments in the dynamic field of bacterial genome reannotation, catalog the important biological roles carried out by small proteins and identify challenges obstructing the way to full understanding of these elusive proteins.

Physiological and evolutionary contexts of a new symbiotic species from the nitrogen-recycling gut community of turtle ants

While genome sequencing has expanded our knowledge of symbiosis, role assignment within multi-species microbiomes remains challenging due to genomic redundancy and the uncertainties of in vivo impacts. We address such questions, here, for a specialized nitrogen (N) recycling microbiome of turtle ants, describing a new genus and species of gut symbiont-Ischyrobacter davidsoniae (Betaproteobacteria: Burkholderiales: Alcaligenaceae)-and its in vivo physiological context. A re-analysis of amplicon sequencing data, with precisely assigned Ischyrobacter reads, revealed a seemingly ubiquitous distribution across the turtle ant genus Cephalotes, suggesting ≥50 million years since domestication. Through new genome sequencing, we also show that divergent I. davidsoniae lineages are conserved in their uricolytic and urea-generating capacities. With phylogenetically refined definitions of Ischyrobacter and separately domesticated Burkholderiales symbionts, our FISH microscopy revealed a distinct niche for I. davidsoniae, with dense populations at the anterior ileum. Being positioned at the site of host N-waste delivery, in vivo metatranscriptomics and metabolomics further implicate I. davidsoniae within a symbiont-autonomous N-recycling pathway. While encoding much of this pathway, I. davidsoniae expressed only a subset of the requisite steps in mature adult workers, including the penultimate step deriving urea from allantoate. The remaining steps were expressed by other specialized gut symbionts. Collectively, this assemblage converts inosine, made from midgut symbionts, into urea and ammonia in the hindgut. With urea supporting host amino acid budgets and cuticle synthesis, and with the ancient nature of other active N-recyclers discovered here, I. davidsoniae emerges as a central player in a conserved and impactful, multipartite symbiosis.

Investigating the eco-evolutionary response of microbiomes to environmental change

Microorganisms are the primary engines of biogeochemical processes and foundational to the provisioning of ecosystem services to human society. Free-living microbial communities (microbiomes) and their functioning are now known to be highly sensitive to environmental change. Given microorganisms' capacity for rapid evolution, evolutionary processes could play a role in this response. Currently, however, few models of biogeochemical processes explicitly consider how microbial evolution will affect biogeochemical responses to environmental change. Here, we propose a conceptual framework for explicitly integrating evolution into microbiome-functioning relationships. We consider how microbiomes respond simultaneously to environmental change via four interrelated processes that affect overall microbiome functioning (physiological acclimation, demography, dispersal and evolution). Recent evidence in both the laboratory and the field suggests that ecological and evolutionary dynamics occur simultaneously within microbiomes; however, the implications for biogeochemistry under environmental change will depend on the timescales over which these processes contribute to a microbiome's response. Over the long term, evolution may play an increasingly important role for microbially driven biogeochemical responses to environmental change, particularly to conditions without recent historical precedent.

Vibrio vulnificus mutation rate: an in vitro approach

Vibrio vulnificus is a multi-host pathogenic species currently subdivided into five phylogenetic lineages (L) plus one pathovar with the ability to infect fish due to a transmissible virulence plasmid. This plasmid (or a fragment of it) has been transmitted between lineages within the species, contributing to the evolution of V. vulnificus. This study aimed to provide an experimental approximation to the V. vulnificus mutation rate by determining spontaneous mutation rates from bacterial cultures of representants of the different lineages by whole-genome sequencing. To this purpose, synonymous SNP differences, i.e., spontaneous mutation not subjected to the evolutive forces, between initial and final culture after serial growth were evaluated and used for mutation rate calculation.

Evidence of cospeciation between termites and their gut bacteria on a geological time scale

Termites host diverse communities of gut microbes, including many bacterial lineages only found in this habitat. The bacteria endemic to termite guts are transmitted via two routes: a vertical route from parent colonies to daughter colonies and a horizontal route between colonies sometimes belonging to different termite species. The relative importance of both transmission routes in shaping the gut microbiota of termites remains unknown. Using bacterial marker genes derived from the gut metagenomes of 197 termites and one Cryptocercus cockroach, we show that bacteria endemic to termite guts are mostly transferred vertically. We identified 18 lineages of gut bacteria showing cophylogenetic patterns with termites over tens of millions of years. Horizontal transfer rates estimated for 16 bacterial lineages were within the range of those estimated for 15 mitochondrial genes, suggesting that horizontal transfers are uncommon and vertical transfers are the dominant transmission route in these lineages. Some of these associations probably date back more than 150 million years and are an order of magnitude older than the cophylogenetic patterns between mammalian hosts and their gut bacteria. Our results suggest that termites have cospeciated with their gut bacteria since first appearing in the geological record.

Neisseria montereyensis sp. nov., Isolated from Oropharynx of California Sea Lion (Zalophus californianus): Genomic, Phylogenetic, and Phenotypic Study

A novel Neisseria strain, designated CSL10203-ORH2^T, was isolated from the oropharynx of a wild California sea lion (Zalophus californianus) that was admitted to The Marine Mammal Center in California, USA. The strain was originally cultured from an oropharyngeal swab on BD Phenylethyl Alcohol (PEA) agar with 5% sheep blood under aerobic conditions. Phylogenetic analyses based on 16S rRNA, rplF, and rpoB gene sequences and the core genome sequences indicated that the strain was most closely related to only N. zalophi CSL 7565^T. The average nucleotide identity and digital DNA-DNA hybridization values between strain CSL10203-ORH2^T and the closely related species N. zalophi CSL 7565^T were 89.84 and 39.70%, respectively, which were significantly lower than the accepted species-defined thresholds for describing novel prokaryotic species at the genomic level. Both type strains were phenotypically similar but can be easily and unambiguously distinguished between each other by the analysis of their housekeeping genes, e.g., rpoB, gyrB, or argF. The major fatty acids in both type strains were C_12:0, C_16:0, C_16:1-c9, and C_18:1-c11. Based on the genomic, phenotypic, and phylogenetic properties, the novel strain represents a novel species of the genus Neisseria, for which the name Neisseria montereyensis sp. nov. with the type strain CSL10203-ORH2^T (= DSM 114706^T = CCUG 76428^T = NCTC 14721^T) is proposed. The genome G + C content is 45.84% and the complete draft genome size is 2,310,535 bp.

Widespread extinctions of co-diversified primate gut bacterial symbionts from humans

Humans and other primates harbour complex gut bacterial communities that influence health and disease, but the evolutionary histories of these symbioses remain unclear. This is partly due to limited information about the microbiota of ancestral primates. Here, using phylogenetic analyses of metagenome-assembled genomes (MAGs), we show that hundreds of gut bacterial clades diversified in parallel (that is, co-diversified) with primate species over millions of years, but that humans have experienced widespread losses of these ancestral symbionts. Analyses of 9,460 human and non-human primate MAGs, including newly generated MAGs from chimpanzees and bonobos, revealed significant co-diversification within ten gut bacterial phyla, including Firmicutes, Actinobacteriota and Bacteroidota. Strikingly, ~44% of the co-diversifying clades detected in African apes were absent from available metagenomic data from humans and ~54% were absent from industrialized human populations. In contrast, only ~3% of non-co-diversifying clades detected in African apes were absent from humans. Co-diversifying clades present in both humans and chimpanzees displayed consistent genomic signatures of natural selection between the two host species but differed in functional content from co-diversifying clades lost from humans, consistent with selection against certain functions. This study discovers host-species-specific bacterial symbionts that predate hominid diversification, many of which have undergone accelerated extinctions from human populations.

Rapid adaptations of Legionella pneumophila to the human host

Legionella pneumophila are host-adapted bacteria that infect and reproduce primarily in amoeboid protists. Using similar infection mechanisms, they infect human macrophages, and cause Legionnaires' disease, an atypical pneumonia, and the milder Pontiac fever. We hypothesized that, despite the similarities in infection mechanisms, the hosts are different enough that there exist high-selective value mutations that would dramatically increase the fitness of Legionella inside the human host. By comparing a large number of isolates from independent infections, we identified two genes, mutated in three unrelated patients, despite the short duration of the incubation period (2-14 days). One is a gene coding for an outer membrane protein (OMP) belonging to the OmpP1/FadL family. The other is a gene coding for an EAL-domain-containing protein involved in cyclic-di-GMP regulation, which in turn modulates flagellar activity. The clinical strain, carrying the mutated EAL-domain-containing homologue, grows faster in macrophages than the wild-type strain, and thus appears to be better adapted to the human host. As human-to-human transmission is very rare, fixation of these mutations into the population and spread into the environment is unlikely. Therefore, parallel evolution - here mutations in the same genes observed in independent human infections - could point to adaptations to the accidental human host. These results suggest that despite the ability of L. pneumophila to infect, replicate in and exit from macrophages, its human-specific adaptations are unlikely to be fixed in the population.

Detection of DNA replication errors and 8-oxo-dGTP-mediated mutations in E. coli by Duplex DNA Sequencing

Mutation is a phenomenon inescapable for all life-forms, including bacteria. While bacterial mutation rates are generally low due to the operation of error-avoidance systems, sometimes they are elevated by many orders of magnitude. Such a state, known as a hypermutable state, can result from exposure to stress or to harmful environments. Studies of bacterial mutation frequencies and analysis of the precise types of mutations can provide insights into the mechanisms by which mutations occur and the possible involvement of error-avoidance pathways. Several approaches have been used for this, like reporter assays involving non-essential genes or mutation accumulation over multiple generations. However, these approaches give an indirect estimation, and a more direct approach for determining mutations is desirable. With the recent development of a DNA sequencing technique known as Duplex Sequencing, it is possible to detect rare variants in a population at a frequency of 1 in 107 base pairs or less. Here, we have applied Duplex Sequencing to study spontaneous mutations in E. coli. We also investigated the production of replication errors by using a mismatch-repair defective (mutL) strain as well as oxidative-stress associated mutations using a mutT-defective strain. For DNA from a wild-type strain we obtained mutant frequencies in the range of 10-7 to 10-8 depending on the specific base-pair substitution, but we argue that these mutants merely represent a background of the system, rather than mutations that occurred in vivo. In contrast, bona-fide in vivo mutations were identified for DNA from both the mutL and mutT strains, as indicated by specific increases in base substitutions that are fully consistent with their established in vivo roles. Notably, the data reproduce the specific context effects of in vivo mutations as well as the leading vs. lagging strand bias among DNA replication errors.

Demographic Expansions and the Emergence of Host Specialization in Genetically Distinct Ecotypes of the Tick-Transmitted Bacterium Anaplasma phagocytophilum

In Europe, genetically distinct ecotypes of the tick-vectored bacterium Anaplasma phagocytophilum circulate among mammals in three discrete enzootic cycles. To date, potential ecological factors that contributed to the emergence of these divergent ecotypes have been poorly studied. Here, we show that the ecotype that predominantly infects roe deer (Capreolus capreolus) is evolutionarily derived. Its divergence from a host generalist ancestor occurred after the last glacial maximum as mammal populations, including roe deer, recolonized the European mainland from southern refugia. We also provide evidence that this host specialist ecotype's effective population size (N_e) has tracked changes in the population of its roe deer host. Specifically, both host and bacterium have undergone substantial increases in N_e over the past 1,500 years. In contrast, we show that while it appears to have undergone a major population expansion starting ~3,500 years ago, in the past 500 years, the contemporary host generalist ecotype has experienced a substantial reduction in genetic diversity levels, possibly as a result of reduced opportunities for transmission between competent hosts. IMPORTANCE The findings of this study reveal specific events important for the evolution of host specialization in a naturally occurring, obligately intracellular bacterial pathogen. Specifically, they show that host range shifts and the emergence of host specialization may occur during periods of population growth in a generalist ancestor. Our results also demonstrate the close correlation between demographic patterns in host and pathogen for a specialist system. These findings have important relevance for understanding the evolution of host range diversity. They may inform future work on host range dynamics, and they provide insights for understanding the emergence of pathogens that have human and veterinary health implications.

Divergent gut microbiota in two closely related house mouse subspecies under common garden conditions

The gastrointestinal microbiota (GM) is considered an important component of the vertebrate holobiont. GM-host interactions influence the fitness of holobionts and are therefore an integral part of evolution. The house mouse is a prominent model for GM-host interactions, and evidence suggests a role for GM in mouse speciation. However, previous studies based on short 16S rRNA GM profiles of wild house mouse subspecies failed to detect GM divergence, which is a prerequisite for the inclusion of GM in Dobzhansky-Muller incompatibilities. Here, we used standard 16S rRNA GM profiling in two mouse subspecies, Mus musculus musculus and M. m. domesticus, including the intestinal mucosa and content of three gut sections (ileum, caecum, and colon). We reduced environmental variability by sampling GM in the offspring of wild mice bred under semi-natural conditions. Although the breeding conditions allowed a contact between the subspecies, we found a clear differentiation of GM between them, in all three gut sections. Differentiation was mainly driven by several Helicobacters and two H. ganmani variants showed a signal of co-divergence with their hosts. Helicobacters represent promising candidates for studying GM-host co-adaptations and the fitness effects of their interactions.

Elevational Gradients Impose Dispersal Limitation on Streptomyces

Dispersal governs microbial biogeography, but the rates and mechanisms of dispersal remain poorly characterized for most microbial taxa. Dispersal limitation is driven by limits on dissemination and establishment, respectively. Elevation gradients create striking patterns of biogeography because they produce steep environmental gradients at small spatial scales, and these gradients offer a powerful tool to examine mechanisms of dispersal limitation. We focus on Streptomyces, a bacterial genus common to soil, by using a taxon-specific phylogenetic marker, the RNA polymerase-encoding rpoB gene. By targeting Streptomyces, we assess dispersal limitation at finer phylogenetic resolution than is possible using whole community analyses. We characterized Streptomyces diversity at local spatial scales (100 to 3,000 m) in two temperate forest sites located in the Adirondacks region of New York State: Woods Lake (<100 m elevation change), and Whiteface Mountain (>1,000 m elevation change). Beta diversity varied considerably at both locations, indicative of dispersal limitation acting at local spatial scales, but beta diversity was significantly higher at Whiteface Mountain. Beta diversity varied across elevation at Whiteface Mountain, being lowest at the mountain’s base. We show that Streptomyces taxa exhibit elevational preferences, and these preferences are phylogenetically conserved. These results indicate that habitat preferences influence Streptomyces biogeography and suggest that barriers to establishment structure Streptomyces communities at higher elevations. These data illustrate that Streptomyces biogeography is governed by dispersal limitation resulting from a complex mixture of stochastic and deterministic processes.

Metaeffector interactions modulate the type III effector-triggered immunity load of Pseudomonas syringae

The bacterial plant pathogen Pseudomonas syringae requires type III secreted effectors (T3SEs) for pathogenesis. However, a major facet of plant immunity entails the recognition of a subset of P. syringae’s T3SEs by intracellular host receptors in a process called Effector-Triggered Immunity (ETI). Prior work has shown that ETI-eliciting T3SEs are pervasive in the P. syringae species complex raising the question of how P. syringae mitigates its ETI load to become a successful pathogen. While pathogens can evade ETI by T3SE mutation, recombination, or loss, there is increasing evidence that effector-effector (a.k.a., metaeffector) interactions can suppress ETI. To study the ETI-suppression potential of P. syringae T3SE repertoires, we compared the ETI-elicitation profiles of two genetically divergent strains: P. syringae pv. tomato DC3000 (PtoDC3000) and P. syringae pv. maculicola ES4326 (PmaES4326), which are both virulent on Arabidopsis thaliana but harbour largely distinct effector repertoires. Of the 529 T3SE alleles screened on A. thaliana Col-0 from the P. syringae T3SE compendium (PsyTEC), 69 alleles from 21 T3SE families elicited ETI in at least one of the two strain backgrounds, while 50 elicited ETI in both backgrounds, resulting in 19 differential ETI responses including two novel ETI-eliciting families: AvrPto1 and HopT1. Although most of these differences were quantitative, three ETI responses were completely absent in one of the pathogenic backgrounds. We performed ETI suppression screens to test if metaeffector interactions contributed to these ETI differences, and found that HopQ1a suppressed AvrPto1m-mediated ETI, while HopG1c and HopF1g suppressed HopT1b-mediated ETI. Overall, these results show that P. syringae strains leverage metaeffector interactions and ETI suppression to overcome the ETI load associated with their native T3SE repertoires.

Landscape of mobile genetic elements and their antibiotic resistance cargo in prokaryotic genomes

Prokaryotic Mobile Genetic Elements (MGEs) such as transposons, integrons, phages and plasmids, play important roles in prokaryotic evolution and in the dispersal of cargo functions like antibiotic resistance. However, each of these MGE types is usually annotated and analysed individually, hampering a global understanding of phylogenetic and environmental patterns of MGE dispersal. We thus developed a computational framework that captures diverse MGE types, their cargos and MGE-mediated horizontal transfer events, using recombinases as ubiquitous MGE marker genes and pangenome information for MGE boundary estimation. Applied to ∼84k genomes with habitat annotation, we mapped 2.8 million MGE-specific recombinases to six operational MGE types, which together contain on average 13% of all the genes in a genome. Transposable elements (TEs) dominated across all taxa (∼1.7 million occurrences), outnumbering phages and phage-like elements (<0.4 million). We recorded numerous MGE-mediated horizontal transfer events across diverse phyla and habitats involving all MGE types, disentangled and quantified the extent of hitchhiking of TEs (17%) and integrons (63%) with other MGE categories, and established TEs as dominant carriers of antibiotic resistance genes. We integrated all these findings into a resource (proMGE.embl.de), which should facilitate future studies on the large mobile part of genomes and its horizontal dispersal.

Reassembling a cannon in the DNA defense arsenal: Genetics of StySA, a BREX phage exclusion system in Salmonella lab strains

Understanding mechanisms that shape horizontal exchange in prokaryotes is a key problem in biology. A major limit on DNA entry is imposed by restriction-modification (RM) processes that depend on the pattern of DNA modification at host-specified sites. In classical RM, endonucleolytic DNA cleavage follows detection of unprotected sites on entering DNA. Recent investigation has uncovered BREX (BacteRiophage EXclusion) systems. These RM-like activities employ host protection by DNA modification, but immediate replication arrest occurs without evident of nuclease action on unmodified phage DNA. Here we show that the historical stySA RM locus of Salmonella enterica sv Typhimurium is a variant BREX system. A laboratory strain disabled for both the restriction and methylation activity of StySA nevertheless has wild type sequence in pglX, the modification gene homolog. Instead, flanking genes pglZ and brxC each carry multiple mutations (μ) in their C-terminal domains. We further investigate this system in situ, replacing the mutated pglZμ and brxCμ genes with the WT counterpart. PglZ-WT supports methylation in the presence of either BrxCμ or BrxC-WT but not in the presence of a deletion/insertion allele, ΔbrxC::cat. Restriction requires both BrxC-WT and PglZ-WT, implicating the BrxC C-terminus specifically in restriction activity. These results suggests that while BrxC, PglZ and PglX are principal components of the BREX modification activity, BrxL is required for restriction only. Furthermore, we show that a partial disruption of brxL disrupts transcription globally.

Horizontal gene transfer is a major driver of evolution and adaptation in bacteria. Genes from outside may be beneficial or dangerous to the receiving cell. Benefits include new food sources such as sugars, or new homes by adhesion, or new resistances, as to antibiotics. Dangers are posed by bacteriophages--viruses that take over the cell machinery, multiply, and release progeny to kill sister cells. Host-dependent restriction-modification systems enable defense that distinguishes relatives from strangers: using a modification pattern (M) carried by DNA bases added by the host cell to prevent restriction (R). Sisters and cousin cells will have the same protective pattern on DNA, while DNA of foreign origin will have the wrong M pattern and be restricted (R, rejected). Typically, restriction involves nuclease digestion. Here we address the enigmatic StySA RM system, one of the earliest to be genetically characterized. It is a variant of the newly recognized defense mechanism, BREX. BREX systems also track DNA history via modification pattern, but restrict by a novel, uncharacterized mechanism. Like other BREX family systems, StySA-BREX modification requires multiple components. When StySA-BREX transcription is unbalanced, we find global disruption of gene transcription. The disruption pattern does not suggest SOS-inducing damage to DNA.

Panoramic Insights into Microevolution and Macroevolution of A Prevotella copri-containing Lineage in Primate Guts

Prevotella copri and its related taxa are widely detected in mammalian gut microbiomes and have been linked with an enterotype in humans. However, their microevolution and macroevolution among hosts are poorly characterized. In this study, extensively collected marker genes and genomes were analyzed to trace their evolutionary history, host specificity, and biogeographic distribution. Investigations based on marker genes and genomes suggest that a P. copri-containing lineage (PCL) harbors diverse species in higher primates. Firstly, P. copri in the human gut consisted of multiple groups exhibiting high genomic divergence and conspicuous but non-strict biogeographic patterns. Most African strains with high genomic divergence from other strains were phylogenetically located at the root of the species, indicating the co-evolutionary history of P. copri and Homo sapiens. Secondly, although long-term co-evolution between PCL and higher primates was revealed, sporadic signals of co-speciation and extensive host jumping of PCL members were suggested among higher primates. Metagenomic and phylogenetic analyses indicated that P. copri and other PCL species found in domesticated mammals had been recently transmitted from humans. Thirdly, strong evidence was found on the extensively horizontal transfer of genes (e.g., genes encoding carbohydrate-active enzymes) among sympatric P. copri groups and PCL species in the same primate host. Our study provides panoramic insights into the combined effects of vertical and horizontal transmission, as well as potential niche adaptation, on the microevolutionary and macroevolutionary history for an enterotype-representative lineage.

Interactions between strains govern the eco-evolutionary dynamics of microbial communities

Genomic data has revealed that genotypic variants of the same species, that is, strains, coexist and are abundant in natural microbial communities. However, it is not clear if strains are ecologically equivalent, and at what characteristic genetic distance they might exhibit distinct interactions and dynamics. Here, we address this problem by tracking 10 taxonomically diverse microbial communities from the pitcher plant Sarracenia purpurea in the laboratory for more than 300 generations. Using metagenomic sequencing, we reconstruct their dynamics over time and across scales, from distant phyla to closely related genotypes. We find that most strains are not ecologically equivalent and exhibit distinct dynamical patterns, often being significantly more correlated with strains from another species than their own. Although even a single mutation can affect laboratory strains, on average, natural strains typically decouple in their dynamics beyond a genetic distance of 100 base pairs. Using mathematical consumer-resource models, we show that these taxonomic patterns emerge naturally from ecological interactions between community members, but only if the interactions are coarse-grained at the level of strains, not species. Finally, by analyzing genomic differences between strains, we identify major functional hubs such as transporters, regulators, and carbohydrate-catabolizing enzymes, which might be the basis for strain-specific interactions. Our work suggests that fine-scale genetic differences in natural communities could be created and stabilized via the rapid diversification of ecological interactions between strains.

Limited Evidence for Microbial Transmission in the Phylosymbiosis between Hawaiian Spiders and Their Microbiota

The degree of similarity between the microbiotas of host species often mirrors the phylogenetic proximity of the hosts. This pattern, referred to as phylosymbiosis, is widespread in animals and plants. While phylosymbiosis was initially interpreted as the signal of symbiotic transmission and coevolution between microbes and their hosts, it is now recognized that similar patterns can emerge even if the microbes are environmentally acquired. Distinguishing between these two scenarios, however, remains challenging. We recently developed HOME (host-microbiota evolution), a cophylogenetic model designed to detect vertically transmitted microbes and host switches from amplicon sequencing data. Here, we applied HOME to the microbiotas of Hawaiian spiders of the genus Ariamnes, which experienced a recent radiation on the archipelago. We demonstrate that although Hawaiian Ariamnes spiders display a significant phylosymbiosis, there is little evidence of microbial vertical transmission. Next, we performed simulations to validate the absence of transmitted microbes in Ariamnes spiders. We show that this is not due to a lack of detection power because of the low number of segregating sites or an effect of phylogenetically driven or geographically driven host switches. Ariamnes spiders and their associated microbes therefore provide an example of a pattern of phylosymbiosis likely emerging from processes other than vertical transmission.

IMPORTANCE How host-associated microbiotas assemble and evolve is one of the outstanding questions of microbial ecology. Studies aiming at answering this question have repeatedly found a pattern of “phylosymbiosis,” that is, a phylogenetic signal in the composition of host-associated microbiotas. While phylosymbiosis was often interpreted as evidence for vertical transmission and host-microbiota coevolution, simulations have now shown that it can emerge from other processes, including host filtering of environmentally acquired microbes. However, distinguishing the processes driving phylosymbiosis in nature remains challenging. We recently developed a cophylogenetic method that can detect vertical transmission. Here, we applied this method to the microbiotas of recently diverged spiders from the Hawaiian archipelago, which display a clear phylosymbiosis pattern. We found that none of the bacterial operational taxonomic units is vertically transmitted. We show with simulations that this result is not due to methodological artifacts. Thus, we provide a striking empirical example of phylosymbiosis emerging from processes other than vertical transmission.

Vertical Inheritance Facilitates Interspecies Diversification in Biosynthetic Gene Clusters and Specialized Metabolites

While specialized metabolites are thought to mediate ecological interactions, the evolutionary processes driving chemical diversification, particularly among closely related lineages, remain poorly understood. Here, we examine the evolutionary dynamics governing the distribution of natural product biosynthetic gene clusters (BGCs) among 118 strains representing all nine currently named species of the marine actinobacterial genus Salinispora. While much attention has been given to the role of horizontal gene transfer (HGT) in structuring BGC distributions, we find that vertical descent facilitates interspecies BGC diversification over evolutionary timescales. Moreover, we identified a distinct phylogenetic signal among Salinispora species at both the BGC and metabolite level, indicating that specialized metabolism represents a conserved phylogenetic trait. Using a combination of genomic analyses and liquid chromatography–high-resolution tandem mass spectrometry (LC-MS/MS) targeting nine experimentally characterized BGCs and their small molecule products, we identified gene gain/loss events, constrained interspecies recombination, and other evolutionary processes associated with vertical inheritance as major contributors to BGC diversification. These evolutionary dynamics had direct consequences for the compounds produced, as exemplified by species-level differences in salinosporamide production. Together, our results support the concept that specialized metabolites, and their cognate BGCs, can represent phylogenetically conserved functional traits with chemical diversification proceeding in species-specific patterns over evolutionary time frames.

A holobiont view of island biogeography: Unravelling patterns driving the nascent diversification of a Hawaiian spider and its microbial associates

The diversification of a host lineage can be influenced by both the external environment and its assemblage of microbes. Here, we use a young lineage of spiders, distributed along a chronologically arranged series of volcanic mountains, to investigate how their associated microbial communities have changed as the spiders colonized new locations. Using the stick spider Ariamnes waikula (Araneae, Theridiidae) on the island of Hawai'i, and outgroup taxa on older islands, we tested whether each component of the "holobiont" (spider hosts, intracellular endosymbionts and gut microbial communities) showed correlated signatures of diversity due to sequential colonization from older to younger volcanoes. To investigate this, we generated ddRAD data for the host spiders and 16S rRNA gene amplicon data from their microbiota. We expected sequential colonizations to result in a (phylo)genetic structuring of the host spiders and in a diversity gradient in microbial communities. The results showed that the host A. waikula is indeed structured by geographical isolation, suggesting sequential colonization from older to younger volcanoes. Similarly, the endosymbiont communities were markedly different between Ariamnes species on different islands, but more homogeneous among A. waikula populations on the island of Hawai'i. Conversely, the gut microbiota, which we suspect is generally environmentally derived, was largely conserved across all populations and species. Our results show that different components of the holobiont respond in distinct ways to the dynamic environment of the volcanic archipelago. This highlights the necessity of understanding the interplay between different components of the holobiont, to properly characterize its evolution.

Additive Uncorrelated Relaxed Clock Models for the Dating of Genomic Epidemiology Phylogenies

Phylogenetic dating is one of the most powerful and commonly used methods of drawing epidemiological interpretations from pathogen genomic data. Building such trees requires considering a molecular clock model which represents the rate at which substitutions accumulate on genomes. When the molecular clock rate is constant throughout the tree then the clock is said to be strict, but this is often not an acceptable assumption. Alternatively, relaxed clock models consider variations in the clock rate, often based on a distribution of rates for each branch. However, we show here that the distributions of rates across branches in commonly used relaxed clock models are incompatible with the biological expectation that the sum of the numbers of substitutions on two neighboring branches should be distributed as the substitution number on a single branch of equivalent length. We call this expectation the additivity property. We further show how assumptions of commonly used relaxed clock models can lead to estimates of evolutionary rates and dates with low precision and biased confidence intervals. We therefore propose a new additive relaxed clock model where the additivity property is satisfied. We illustrate the use of our new additive relaxed clock model on a range of simulated and real data sets, and we show that using this new model leads to more accurate estimates of mean evolutionary rates and ancestral dates.

Major genetic discontinuity and novel toxigenic species in Clostridioides difficile taxonomy

Clostridioides difficile infection (CDI) remains an urgent global One Health threat. The genetic heterogeneity seen across C. difficile underscores its wide ecological versatility and has driven the significant changes in CDI epidemiology seen in the last 20 years. We analysed an international collection of over 12,000 C. difficile genomes spanning the eight currently defined phylogenetic clades. Through whole-genome average nucleotide identity, and pangenomic and Bayesian analyses, we identified major taxonomic incoherence with clear species boundaries for each of the recently described cryptic clades CI–III. The emergence of these three novel genomospecies predates clades C1–5 by millions of years, rewriting the global population structure of C. difficile specifically and taxonomy of the Peptostreptococcaceae in general. These genomospecies all show unique and highly divergent toxin gene architecture, advancing our understanding of the evolution of C. difficile and close relatives. Beyond the taxonomic ramifications, this work may impact the diagnosis of CDI.

Streptococcus vicugnae sp. nov., isolated from faeces of alpacas (Vicugna pacos) and cattle (Bos taurus), Streptococcus zalophi sp. nov., and Streptococcus pacificus sp. nov., isolated from respiratory tract of California sea lions (Zalophus californianus)

Four novel independent strains of Streptococcus spp. were isolated from faeces of alpaca (SL1232^T), cattle (KCJ4950), and from respiratory tract of wild California sea lions (CSL7508^T, CSL7591^T). The strains were indole-, oxidase- and catalase-negative, non-spore-forming, non-motile Gram-positive cocci in short and long chains, facultative anaerobes. The 16S rRNA gene of SL1232^T and KCJ4950 shared 99.40-99.60% nucleotide similarity to strains of S. equinus, S. lutetiensis, S. infantarius, and the 16S rRNA gene of CSL7508^T and CSL7591^T demonstrated 98.72 and 98.92% similarity, respectively, to S. marimammalium. All other known Streptococcus species had the 16S rRNA gene sequence similarities of ≤95%. The genomes were sequenced for the novel strains. Average nucleotide identity (ANI) analysis for strains SL1232^T and KCJ4950, showed the highest similarity to S. equinus, S. lutetiensis, and S. infantarius with 85.21, 87.17, 88.47, 85.54, 87.47 and 88.89%, respectively, and strains CSL7508^T and CSL7591^T to S. marimammalium with 87.16 and 83.97%, respectively. Results of ANI were confirmed by pairwise digital DNA-DNA hybridization and phylogeny, which also revealed that the strains belong to three novel species of the genus Streptococcus. Phenotypical features of the novel species were in congruence with closely related members of the genus Streptococcus and gave negative reactions with the tested Lancefield serological groups (A-D, F and G). MALDI-TOF mass spectrometry supported identification of the species. Based on these data, we propose three novel species of the genus Streptococcus, for which the name Streptococcus vicugnae sp. nov. is proposed with the type strain SL1232^T (=NCTC 14341^T=DSM 110741^T=CCUG 74371^T), Streptococcus zalophi sp. nov. is proposed with the type strain CSL7508^T (=NCTC 14410^T=DSM 110742^T=CCUG 74374^T) and Streptococcus pacificus sp. nov. is proposed with the type strain CSL7591^T (=NCTC 14455^T=DSM 111148^T=CCUG 74655^T). The genome G+C content is 36.89, 34.85, and 35.34 % and draft genome sizes are 1906993, 1581094 and 1656080 bp for strains SL1232^T, CSL7508^T, and CSL7591^T, respectively.

Adaptive differentiation and rapid evolution of a soil bacterium along a climate gradient

Increasing evidence suggests that evolutionary processes frequently shape ecological patterns; however, most microbiome studies thus far have focused on only the ecological responses of these communities. By using parallel field experiments and focusing in on a model soil bacterium, we showed that bacterial “species” are differentially adapted to local climates, leading to changes in their composition. Furthermore, we detected strain-level evolution, providing direct evidence that both ecological and evolutionary processes operate on annual timescales. The consideration of eco-evolutionary dynamics may therefore be important to understand the response of soil microbiomes to future environmental change.

Microbial community responses to environmental change are largely associated with ecological processes; however, the potential for microbes to rapidly evolve and adapt remains relatively unexplored in natural environments. To assess how ecological and evolutionary processes simultaneously alter the genetic diversity of a microbiome, we conducted two concurrent experiments in the leaf litter layer of soil over 18 mo across a climate gradient in Southern California. In the first experiment, we reciprocally transplanted microbial communities from five sites to test whether ecological shifts in ecotypes of the abundant bacterium, Curtobacterium, corresponded to past adaptive differentiation. In the transplanted communities, ecotypes converged toward that of the native communities growing on a common litter substrate. Moreover, these shifts were correlated with community-weighted mean trait values of the Curtobacterium ecotypes, indicating that some of the trait variation among ecotypes could be explained by local adaptation to climate conditions. In the second experiment, we transplanted an isogenic Curtobacterium strain and tracked genomic mutations associated with the sites across the same climate gradient. Using a combination of genomic and metagenomic approaches, we identified a variety of nonrandom, parallel mutations associated with transplantation, including mutations in genes related to nutrient acquisition, stress response, and exopolysaccharide production. Together, the field experiments demonstrate how both demographic shifts of previously adapted ecotypes and contemporary evolution can alter the diversity of a soil microbiome on the same timescale.

Metabacillus schmidteae sp. nov., Cultivated from Planarian Schmidtea mediterranea Microbiota

Taxonogenomics combines phenotypic assays and genomic analysis as a means of characterizing novel strains. We used this strategy to study Marseille-P9898T strain, an aerobic, motile, Gram-negative, spore-forming, and rod-shaped bacterium isolated from planarian Schmidtea mediterranea. Marseille-P9898T is catalase-positive and oxidase-negative. The major fatty acids detected are 12-methyl-tetradecanoic acid, 13-methyl-tetradecanoic acid, and hexadecanoic acid. Marseille-P9898T strain shared more than 98% sequence similarity with the Metabacillus niabensis strain 4T19T (98.99%), Metabacillus halosaccharovorans strain E33T (98.75%), Metabacillus malikii strain NCCP-662T (98.19%), and Metabacillus litoralis strain SW-211T (97.15%). Marseille-P9898 strain belongs to Metabacillus genus. Genomic analysis revealed the highest similarities with Ortho-ANI and dDDH, 85.76% with Metabacillus halosaccharovorans, and 34.20% with Bacillus acidicola, respectively. These results show that the Marseille-P9898T strain is a novel bacterial species from Metabacillus genus, for which we propose the name of Metabacillus schmidteae sp. nov. (Type strain Marseille-P9898T = CSUR P9898T = DSM 111480T).

Pedobacter ghigonii sp. nov., Isolated from the Microbiota of the Planarian Schmidtea mediterranea

The planarian S. mediterranea is a platyhelminth with worldwide distribution that can regenerate any part of its body after amputation and has the capacity to eliminate a large spectrum of human bacterial pathogens. Surprisingly, the microbiota of S. mediterranea remains poorly investigated. Using the culturomics strategy to study the bacterial component of planarians, we isolated a new bacterial strain, Marseille-Q2390, which we characterized with the taxono-genomic approach that associates phenotypic assays and genome sequencing and analysis. Strain Marseille-Q2390 exhibited a 16S rRNA sequence similarity of 99.36% with Pedobacter kyungheensis strain THG-T17T, the closest phylogenetic neighbor. It is a white-pigmented, Gram-negative, and rod-shaped bacterium. It grows in aerobic conditions and belongs to the family Sphingobacteriaceae. The genome of strain Marseille-Q2390 is 5,919,359 bp-long, with a G + C content of 40.3%. By comparing its genome with other closely related strains, the highest Orthologous Average Nucleotide Identity (Ortho-ANI) and digital DNA-DNA hybridization (dDDH) values were 85.71% and 30.50%, respectively, which were found with Pedobacter soli strain 15-51T. We conclude that strain Marseille-Q2390T is sufficiently different from other nearby species to be classified within a new species for which we propose the name Pedobacter ghigonii sp. nov.

How sequence populations persist inside bacterial genomes

Compared to their eukaryotic counterparts, bacterial genomes are small and contain extremely tightly packed genes. Repetitive sequences are rare but not completely absent. One of the most common repeat families is REPINs. REPINs can replicate in the host genome and form populations that persist for millions of years. Here, we model the interactions of these intragenomic sequence populations with the bacterial host. We first confirm well-established results, in the presence and absence of horizontal gene transfer (hgt) sequence populations either expand until they drive the host to extinction or the sequence population gets purged from the genome. We then show that a sequence population can be stably maintained, when each individual sequence provides a benefit that decreases with increasing sequence population size. Maintaining a sequence population of stable size also requires the replication of the sequence population to be costly to the host, otherwise the sequence population size will increase indefinitely. Surprisingly, in regimes with high hgt rates, the benefit conferred by the sequence population does not have to exceed the damage it causes to its host. Our analyses provide a plausible scenario for the persistence of sequence populations in bacterial genomes. We also hypothesize a limited biologically relevant parameter range for the provided benefit, which can be tested in future experiments.

Host reproductive cycle influences the pouch microbiota of wild southern hairy-nosed wombats (Lasiorhinus latifrons)

Background

Marsupials are born much earlier than placental mammals, with most crawling from the birth canal to the protective marsupium (pouch) to further their development. However, little is known about the microbiology of the pouch and how it changes throughout a marsupial’s reproductive cycle. Here, using stringent controls, we characterized the microbial composition of multiple body sites from 26 wild Southern Hairy-nosed Wombats (SHNWs), including pouch samples from animals at different reproductive stages.

Results

Using qPCR of the 16S rRNA gene we detected a microbial community in the SHNW pouch. We observed significant differences in microbial composition and diversity between the body sites tested, as well as between pouch samples from different reproductive stages. The pouches of reproductively active females had drastically lower microbial diversity (mean ASV richness 19 ± 8) compared to reproductively inactive females (mean ASV richness 941 ± 393) and were dominated by gram positive bacteria from the Actinobacteriota phylum (81.7–90.6%), with the dominant families classified as Brevibacteriaceae, Corynebacteriaceae, Microbacteriaceae, and Dietziaceae. Three of the five most abundant sequences identified in reproductively active pouches had closest matches to microbes previously isolated from tammar wallaby pouches.

Conclusions

This study represents the first contamination-controlled investigation into the marsupial pouch microbiota, and sets a rigorous framework for future pouch microbiota studies. Our results indicate that SHNW pouches contain communities of microorganisms that are substantially altered by the host reproductive cycle. We recommend further investigation into the roles that pouch microorganisms may play in marsupial reproductive health and joey survival.

Supplementary Information

The online version contains supplementary material available at 10.1186/s42523-021-00074-8.

The biogeography of Streptomyces in New Zealand enabled by high-throughput sequencing of genus-specific rpoB amplicons

We evaluated Streptomyces biogeography in soils along a 1200 km latitudinal transect across New Zealand (NZ). Streptomyces diversity was examined using high-throughput sequencing of rpoB amplicons generated with a Streptomyces specific primer set. We detected 1287 Streptomyces rpoB operational taxonomic units (OTUs) with 159 ± 92 (average ± SD) rpoB OTUs per site. Only 12% (n = 149) of these OTUs matched rpoB sequences from cultured specimens (99% nucleotide identity cutoff). Streptomyces phylogenetic diversity (Faith's PD) was correlated with soil pH, mean annual temperature and plant community richness (Spearman's r: 0.77, 0.64 and -0.79, respectively; P < 0.05), but not with latitude. In addition, soil pH and plant community richness both explained significant variation in Streptomyces beta diversity. Streptomyces communities exhibited both high dissimilarity and strong dominance of one or a few species at each site. Taken together, these results suggest that dispersal limitation due to competitive interactions limits the colonization success of spores that relocate to new sites. Cultivated Streptomyces isolates represent a major source of clinically useful antibiotics, but only a small fraction of extant diversity within the genus have been identified and most species of Streptomyces have yet to be described.

Atribacteria Reproducing over Millions of Years in the Atlantic Abyssal Subseafloor

How microbial metabolism is translated into cellular reproduction under energy-limited settings below the seafloor over long timescales is poorly understood. Here, we show that microbial abundance increases an order of magnitude over a 5 million-year-long sequence in anoxic subseafloor clay of the abyssal North Atlantic Ocean. This increase in biomass correlated with an increased number of transcribed protein-encoding genes that included those involved in cytokinesis, demonstrating that active microbial reproduction outpaces cell death in these ancient sediments. Metagenomes, metatranscriptomes, and 16S rRNA gene sequencing all show that the actively reproducing community was dominated by the candidate phylum "Candidatus Atribacteria," which exhibited patterns of gene expression consistent with fermentative, and potentially acetogenic, metabolism. "Ca. Atribacteria" dominated throughout the 8 million-year-old cored sequence, despite the detection limit for gene expression being reached in 5 million-year-old sediments. The subseafloor reproducing "Ca. Atribacteria" also expressed genes encoding a bacterial microcompartment that has potential to assist in secondary fermentation by recycling aldehydes and, thereby, harness additional power to reduce ferredoxin and NAD⁺ Expression of genes encoding the Rnf complex for generation of chemiosmotic ATP synthesis were also detected from the subseafloor "Ca Atribacteria," as well as the Wood-Ljungdahl pathway that could potentially have an anabolic or catabolic function. The correlation of this metabolism with cytokinesis gene expression and a net increase in biomass over the million-year-old sampled interval indicates that the "Ca Atribacteria" can perform the necessary catabolic and anabolic functions necessary for cellular reproduction, even under energy limitation in millions-of-years-old anoxic sediments.IMPORTANCE The deep subseafloor sedimentary biosphere is one of the largest ecosystems on Earth, where microbes subsist under energy-limited conditions over long timescales. It remains poorly understood how mechanisms of microbial metabolism promote increased fitness in these settings. We discovered that the candidate bacterial phylum "Candidatus Atribacteria" dominated a deep-sea subseafloor ecosystem, where it exhibited increased transcription of genes associated with acetogenic fermentation and reproduction in million-year-old sediment. We attribute its improved fitness after burial in the seabed to its capabilities to derive energy from increasingly oxidized metabolites via a bacterial microcompartment and utilize a potentially reversible Wood-Ljungdahl pathway to help meet anabolic and catabolic requirements for growth. Our findings show that "Ca Atribacteria" can perform all the necessary catabolic and anabolic functions necessary for cellular reproduction, even under energy limitation in anoxic sediments that are millions of years old.

Selective Eradication of Staphylococcus aureus by the Designer Genetically Programmed Yeast Biocontrol Agent

Staphylococcus aureus is a common human pathogen that is particularly often associated with antibiotic resistance. The eradication of this ubiquitous infectious agent from its ecological niches and contaminated surfaces is especially complicated by excessive biofilm formation and persisting cells, which evade the antibacterial activity of conventional antibiotics. Here, we present an alternative view of the problem of specific S. aureus eradication. The constitutive heterologous production of highly specific bacteriolytic protease lysostaphin in yeast Pichia pastoris provides an efficient biocontrol agent, specifically killing S. aureus in coculture. A yeast-based anti-S. aureus probiotic was efficient in a high range of temperatures and target-to-effector ratios, indicating its robustness and versatility in eliminating S. aureus cells. The efficient eradication of S. aureus by live lysostaphin-producing P. pastoris was achieved at high scales, providing a simple, biocompatible and cost-effective strategy for S. aureus lysis in bioproduction and surface decontamination. Future biomedical applications based on designer yeast biocontrol agents require evaluation in in vivo models. However, we believe that this strategy is very promising since it provides highly safe, efficient and selective genetically programmed probiotics and targeted biocontrol agents.

A seventeenth-century Mycobacterium tuberculosis genome supports a Neolithic emergence of the Mycobacterium tuberculosis complex

BACKGROUND

Although tuberculosis accounts for the highest mortality from a bacterial infection on a global scale, questions persist regarding its origin. One hypothesis based on modern Mycobacterium tuberculosis complex (MTBC) genomes suggests their most recent common ancestor followed human migrations out of Africa approximately 70,000 years before present. However, studies using ancient genomes as calibration points have yielded much younger dates of less than 6000 years. Here, we aim to address this discrepancy through the analysis of the highest-coverage and highest-quality ancient MTBC genome available to date, reconstructed from a calcified lung nodule of Bishop Peder Winstrup of Lund (b. 1605-d. 1679).

RESULTS

A metagenomic approach for taxonomic classification of whole DNA content permitted the identification of abundant DNA belonging to the human host and the MTBC, with few non-TB bacterial taxa comprising the background. Genomic enrichment enabled the reconstruction of a 141-fold coverage M. tuberculosis genome. In utilizing this high-quality, high-coverage seventeenth-century genome as a calibration point for dating the MTBC, we employed multiple Bayesian tree models, including birth-death models, which allowed us to model pathogen population dynamics and data sampling strategies more realistically than those based on the coalescent.

CONCLUSIONS

The results of our metagenomic analysis demonstrate the unique preservation environment calcified nodules provide for DNA. Importantly, we estimate a most recent common ancestor date for the MTBC of between 2190 and 4501 before present and for Lineage 4 of between 929 and 2084 before present using multiple models, confirming a Neolithic emergence for the MTBC.