Unveiling intraspecific diversity and evolutionary dynamics of the foodborne pathogen Bacillus paranthracis through high-quality pan-genome analysis

Understanding the evolutionary dynamics of foodborne pathogens throughout host-associated habitats is of utmost importance. Bacterial pan-genomes, as dynamic entities, are strongly influenced by ecological lifestyles. As a phenotypically diverse species in the Bacillus cereus group, Bacillus paranthracis is recognized as an emerging foodborne pathogen and a probiotic simultaneously. This poorly understood species is a suitable study model for adaptive pan-genome evolution. In this study, we determined the biogeographic distribution, abundance, genetic diversity, and genotypic profiles of key genetic elements of B. paranthracis. Metagenomic read recruitment analyses demonstrated that B. paranthracis members are globally distributed and abundant in host-associated habitats. A high-quality pan-genome of B. paranthracis was subsequently constructed to analyze the evolutionary dynamics involved in ecological adaptation comprehensively. The open pan-genome indicated a flexible gene repertoire with extensive genetic diversity. Significant divergences in the phylogenetic relationships, functional enrichment, and degree of selective pressure between the different components demonstrated different evolutionary dynamics between the core and accessory genomes driven by ecological forces. Purifying selection and gene loss are the main signatures of evolutionary dynamics in B. paranthracis pan-genome. The plasticity of the accessory genome is characterized by horizontal gene transfer (HGT), massive gene losses, and weak purifying or positive selection, which might contribute to niche-specific adaptation. In contrast, although the core genome dominantly undergoes purifying selection, its association with HGT and positively selected mutations indicates its potential role in ecological diversification. Furthermore, host fitness-related dynamics are characterized by the loss of secondary metabolite biosynthesis gene clusters (BGCs) and CAZyme-encoding genes and the acquisition of antimicrobial resistance (AMR) and virulence genes via HGT. This study offers a case study of pan-genome evolution to investigate the ecological adaptations reflected by biogeographical characteristics, thereby advancing the understanding of intraspecific diversity and evolutionary dynamics of foodborne pathogens.

Transporter Proteins as Ecological Assets and Features of Microbial Eukaryotic Pangenomes

Here we review two connected themes in evolutionary microbiology: (a) the nature of gene repertoire variation within species groups (pangenomes) and (b) the concept of metabolite transporters as accessory proteins capable of providing niche-defining "bolt-on" phenotypes. We discuss the need for improved sampling and understanding of pangenome variation in eukaryotic microbes. We then review the factors that shape the repertoire of accessory genes within pangenomes. As part of this discussion, we outline how gene duplication is a key factor in both eukaryotic pangenome variation and transporter gene family evolution. We go on to outline how, through functional characterization of transporter-encoding genes, in combination with analyses of how transporter genes are gained and lost from accessory genomes, we can reveal much about the niche range, the ecology, and the evolution of virulence of microbes. We advocate for the coordinated systematic study of eukaryotic pangenomes through genome sequencing and the functional analysis of genes found within the accessory gene repertoire.

Neutral evolution of cellular phenotypes

Eukaryotes exhibit a great diversity of cellular and subcellular morphologies, but their basic underlying architecture is fairly constant. All have a nucleus, Golgi, cytoskeleton, plasma membrane, vesicles, ribosomes, and all known lineages but one have mitochondrion-related organelles. Moreover, most eukaryotes undergo processes such as mitosis, meiosis, DNA recombination, and often perform feats such as phagocytosis, and amoeboid and flagellar movement. With all of these commonalities, it is obvious that eukaryotes evolved from a common ancestor, but it is not obvious how eukaryotes came to have their diverse structural phenotypes. Are these phenotypes adaptations to particular niches, their evolution dominated by positive natural selection? Or is eukaryotic cellular diversity substantially the product of neutral evolutionary processes, with adaptation either illusory or a secondary consequence? In this paper, we outline how a hierarchical view of phenotype can be used to articulate a neutral theory of phenotypic evolution, involving processes such as gene loss, gene replacement by homologues or analogues, gene duplication followed by subfunctionalization, and constructive neutral evolution. We suggest that neutral iterations of these processes followed by entrenchment of their products can explain much of the diversity of cellular, developmental, and biochemical phenotypes of unicellular eukaryotes and should be explored in addition to adaptive explanations.

An evolutionary ratchet leading to loss of elongation factors in eukaryotes

Background

The GTPase eEF1A is the eukaryotic factor responsible for the essential, universal function of aminoacyl-tRNA delivery to the ribosome. Surprisingly, eEF1A is not universally present in eukaryotes, being replaced by the paralog EFL independently in multiple lineages. The driving force behind this unusually frequent replacement is poorly understood.

Results

Through sequence searching of genomic and EST databases, we find a striking association of eEF1A replacement by EFL and loss of eEF1A’s guanine exchange factor, eEF1Bα, suggesting that EFL is able to spontaneously recharge with GTP. Sequence conservation and homology modeling analyses indicate several sequence regions that may be responsible for EFL’s lack of requirement for eEF1Bα.

Conclusions

We propose that the unusual pattern of eEF1A, eEF1Bα and EFL presence and absence can be explained by a ratchet-like process: if either eEF1A or eEF1Bα diverges beyond functionality in the presence of EFL, the system is unable to return to the ancestral, eEF1A:eEFBα-driven state.

The evolution of paralogous enzymes MAT and MATX within the Euglenida and beyond

Background

Methionine adenosyltransferase (MAT) is a ubiquitous essential enzyme that, in eukaryotes, occurs in two relatively divergent paralogues: MAT and MATX. MATX has a punctate distribution across the tree of eukaryotes and, except for a few cases, is mutually exclusive with MAT. This phylogenetic pattern could have arisen by either differential loss of old paralogues or the spread of one of these paralogues by horizontal gene transfer. Our aim was to map the distribution of MAT/MATX genes within the Euglenida in order to more comprehensively characterize the evolutionary history of MATX.

Results

We generated 26 new sequences from 23 different lineages of euglenids and one prasinophyte alga Pyramimonas parkeae. MATX was present only in photoautotrophic euglenids. The mixotroph Rapaza viridis and the prasinophyte alga Pyramimonas parkeae, which harbors chloroplasts that are most closely related to the chloroplasts in photoautotrophic euglenids, both possessed only the MAT paralogue. We found both the MAT and MATX paralogues in two photoautotrophic species (Phacus orbicularis and Monomorphina pyrum). The significant conflict between eukaryotic phylogenies inferred from MATX and SSU rDNA data represents strong evidence that MATX paralogues have undergone horizontal gene transfer across the tree of eukaryotes.

Conclusions

Our results suggest that MATX entered the euglenid lineage in a single horizontal gene transfer event that took place after the secondary endosymbiotic origin of the euglenid chloroplast. The origin of the MATX paralogue is unclear, and it cannot be excluded that it arose by a gene duplication event before the most recent common ancestor of eukaryotes.

Horizontal gene transfer in eukaryotes: the weak-link model

The significance of horizontal gene transfer (HGT) in eukaryotic evolution remains controversial. Although many eukaryotic genes are of bacterial origin, they are often interpreted as being derived from mitochondria or plastids. Because of their fixed gene pool and gene loss, however, mitochondria and plastids alone cannot adequately explain the presence of all, or even the majority, of bacterial genes in eukaryotes. Available data indicate that no insurmountable barrier to HGT exists, even in complex multicellular eukaryotes. In addition, the discovery of both recent and ancient HGT events in all major eukaryotic groups suggests that HGT has been a regular occurrence throughout the history of eukaryotic evolution. A model of HGT is proposed that suggests both unicellular and early developmental stages as likely entry points for foreign genes into multicellular eukaryotes.

Parallel re-modeling of EF-1α function: divergent EF-1α genes co-occur with EFL genes in diverse distantly related eukaryotes

BACKGROUND

Elongation factor-1α (EF-1α) and elongation factor-like (EFL) proteins are functionally homologous to one another, and are core components of the eukaryotic translation machinery. The patchy distribution of the two elongation factor types across global eukaryotic phylogeny is suggestive of a 'differential loss' hypothesis that assumes that EF-1α and EFL were present in the most recent common ancestor of eukaryotes followed by independent differential losses of one of the two factors in the descendant lineages. To date, however, just one diatom and one fungus have been found to have both EF-1α and EFL (dual-EF-containing species).

RESULTS

In this study, we characterized 35 new EF-1α/EFL sequences from phylogenetically diverse eukaryotes. In so doing we identified 11 previously unreported dual-EF-containing species from diverse eukaryote groups including the Stramenopiles, Apusomonadida, Goniomonadida, and Fungi. Phylogenetic analyses suggested vertical inheritance of both genes in each of the dual-EF lineages. In the dual-EF-containing species we identified, the EF-1α genes appeared to be highly divergent in sequence and suppressed at the transcriptional level compared to the co-occurring EFL genes.

CONCLUSIONS

According to the known EF-1α/EFL distribution, the differential loss process should have occurred independently in diverse eukaryotic lineages, and more dual-EF-containing species remain unidentified. We predict that dual-EF-containing species retain the divergent EF-1α homologues only for a sub-set of the original functions. As the dual-EF-containing species are distantly related to each other, we propose that independent re-modelling of EF-1α function took place in multiple branches in the tree of eukaryotes.

Climate change impacts on the tree of life: changes in phylogenetic diversity illustrated for acropora corals

The possible loss of whole branches from the tree of life is a dramatic, but under-studied, biological implication of climate change. The tree of life represents an evolutionary heritage providing both present and future benefits to humanity, often in unanticipated ways. Losses in this evolutionary (evo) life-support system represent losses in "evosystem" services, and are quantified using the phylogenetic diversity (PD) measure. High species-level biodiversity losses may or may not correspond to high PD losses. If climate change impacts are clumped on the phylogeny, then loss of deeper phylogenetic branches can mean disproportionately large PD loss for a given degree of species loss. Over time, successive species extinctions within a clade each may imply only a moderate loss of PD, until the last species within that clade goes extinct, and PD drops precipitously. Emerging methods of "phylogenetic risk analysis" address such phylogenetic tipping points by adjusting conservation priorities to better reflect risk of such worst-case losses. We have further developed and explored this approach for one of the most threatened taxonomic groups, corals. Based on a phylogenetic tree for the corals genus Acropora, we identify cases where worst-case PD losses may be avoided by designing risk-averse conservation priorities. We also propose spatial heterogeneity measures changes to assess possible changes in the geographic distribution of corals PD.