Genes under positive selection in Escherichia coli

Contact-dependent growth inhibition (CDI) systems deploy a large family of polymorphic ionophoric toxins for inter-bacterial competition

Contact-dependent growth inhibition (CDI) is a widespread form of inter-bacterial competition mediated by CdiA effector proteins. CdiA is presented on the inhibitor cell surface and delivers its toxic C-terminal region (CdiA-CT) into neighboring bacteria upon contact. Inhibitor cells also produce CdiI immunity proteins, which neutralize CdiA-CT toxins to prevent auto-inhibition. Here, we describe a diverse group of CDI ionophore toxins that dissipate the transmembrane potential in target bacteria. These CdiA-CT toxins are composed of two distinct domains based on AlphaFold2 modeling. The C-terminal ionophore domains are all predicted to form five-helix bundles capable of spanning the cell membrane. The N-terminal "entry" domains are variable in structure and appear to hijack different integral membrane proteins to promote toxin assembly into the lipid bilayer. The CDI ionophores deployed by E. coli isolates partition into six major groups based on their entry domain structures. Comparative sequence analyses led to the identification of receptor proteins for ionophore toxins from groups 1 & 3 (AcrB), group 2 (SecY) and groups 4 (YciB). Using forward genetic approaches, we identify novel receptors for the group 5 and 6 ionophores. Group 5 exploits homologous putrescine import proteins encoded by puuP and plaP, and group 6 toxins recognize di/tripeptide transporters encoded by paralogous dtpA and dtpB genes. Finally, we find that the ionophore domains exhibit significant intra-group sequence variation, particularly at positions that are predicted to interact with CdiI. Accordingly, the corresponding immunity proteins are also highly polymorphic, typically sharing only ~30% sequence identity with members of the same group. Competition experiments confirm that the immunity proteins are specific for their cognate ionophores and provide no protection against other toxins from the same group. The specificity of this protein interaction network provides a mechanism for self/nonself discrimination between E. coli isolates.

Comparative genomics analysis and characterization of Shiga toxin-producing Escherichia coli O157:H7 strains reveal virulence genes, resistance genes, prophages and plasmids

Escherichia coli O157:H7 is a foodborne pathogen that has been linked to global disease outbreaks. These diseases include hemorrhagic colitis and hemolytic uremic syndrome. It is vital to know the features that make this strain pathogenic to understand the development of disease outbreaks. In the current study, a comparative genomic analysis was carried out to determine the presence of structural and functional features of O157:H7 strains obtained from 115 National Center for Biotechnology Information database. These strains of interest were analysed in the following programs: BLAST Ring Image Generator, PlasmidFinder, ResFinder, VirulenceFinder, IslandViewer 4 and PHASTER. Five strains (ECP19-198, ECP19-798, F7508, F8952, H2495) demonstrated a great homology with Sakai because of a few regions missing. Five resistant genes were identified, however, Macrolide-associated resistance gene mdf(A) was commonly found in all genomes. Majority of the strains (97%) were positive for 15 of the virulent genes (espA, espB, espF, espJ, gad, chuA, eae, iss, nleA, nleB, nleC, ompT, tccP, terC and tir). The plasmid analysis demonstrated that the IncF group was the most prevalent in the strains analysed. The prophage and genomic island analysis showed a distribution of bacteriophages and genomic islands respectively. The results indicated that structural and functional features of the many O157:H7 strains differ and may be a result of obtaining mobile genetic elements via horizontal gene transfer. Understanding the evolution of O157:H7 strains pathogenicity in terms of their structural and functional features will enable the development of detection and control of transmission strategies.

Lower viral evolutionary pressure under stable versus fluctuating conditions in subzero Arctic brines

BACKGROUND

Climate change threatens Earth's ice-based ecosystems which currently offer archives and eco-evolutionary experiments in the extreme. Arctic cryopeg brine (marine-derived, within permafrost) and sea ice brine, similar in subzero temperature and high salinity but different in temporal stability, are inhabited by microbes adapted to these extreme conditions. However, little is known about their viruses (community composition, diversity, interaction with hosts, or evolution) or how they might respond to geologically stable cryopeg versus fluctuating sea ice conditions.

RESULTS

We used long- and short-read viromics and metatranscriptomics to study viruses in Arctic cryopeg brine, sea ice brine, and underlying seawater, recovering 11,088 vOTUs (~species-level taxonomic unit), a 4.4-fold increase of known viruses in these brines. More specifically, the long-read-powered viromes doubled the number of longer (≥25 kb) vOTUs generated and recovered more hypervariable regions by >5-fold compared to short-read viromes. Distribution assessment, by comparing to known viruses in public databases, supported that cryopeg brine viruses were of marine origin yet distinct from either sea ice brine or seawater viruses, while 94% of sea ice brine viruses were also present in seawater. A virus-encoded, ecologically important exopolysaccharide biosynthesis gene was identified, and many viruses (~half of metatranscriptome-inferred "active" vOTUs) were predicted as actively infecting the dominant microbial genera Marinobacter and Polaribacter in cryopeg and sea ice brines, respectively. Evolutionarily, microdiversity (intra-species genetic variations) analyses suggested that viruses within the stable cryopeg brine were under significantly lower evolutionary pressures than those in the fluctuating sea ice environment, while many sea ice brine virus-tail genes were under positive selection, indicating virus-host co-evolutionary arms races.

CONCLUSIONS

Our results confirmed the benefits of long-read-powered viromics in understanding the environmental virosphere through significantly improved genomic recovery, expanding viral discovery and the potential for biological inference. Evidence of viruses actively infecting the dominant microbes in subzero brines and modulating host metabolism underscored the potential impact of viruses on these remote and underexplored extreme ecosystems. Microdiversity results shed light on different strategies viruses use to evolve and adapt when extreme conditions are stable versus fluctuating. Together, these findings verify the value of long-read-powered viromics and provide foundational data on viral evolution and virus-microbe interactions in Earth's destabilized and rapidly disappearing cryosphere. Video Abstract.

Panguiarchaeum symbiosum, a potential hyperthermophilic symbiont in the TACK superphylum

The biology of Korarchaeia remains elusive due to the lack of genome representatives. Here, we reconstruct 10 closely related metagenome-assembled genomes from hot spring habitats and place them into a single species, proposed herein as Panguiarchaeum symbiosum. Functional investigation suggests that Panguiarchaeum symbiosum is strictly anaerobic and grows exclusively in thermal habitats by fermenting peptides coupled with sulfide and hydrogen production to dispose of electrons. Due to its inability to biosynthesize archaeal membranes, amino acids, and purines, this species likely exists in a symbiotic lifestyle similar to DPANN archaea. Population metagenomics and metatranscriptomic analyses demonstrated that genes associated with amino acid/peptide uptake and cell attachment exhibited positive selection and were highly expressed, supporting the proposed proteolytic catabolism and symbiotic lifestyle. Our study sheds light on the metabolism, evolution, and potential symbiotic lifestyle of Panguiarchaeum symbiosum, which may be a unique host-dependent archaeon within the TACK superphylum.

Genomic profiling of Nitrospira species reveals ecological success of comammox Nitrospira

BACKGROUND

The discovery of microorganisms capable of complete ammonia oxidation to nitrate (comammox) has prompted a paradigm shift in our understanding of nitrification, an essential process in N cycling, hitherto considered to require both ammonia oxidizing and nitrite oxidizing microorganisms. This intriguing metabolism is unique to the genus Nitrospira, a diverse taxon previously known to only contain canonical nitrite oxidizers. Comammox Nitrospira have been detected in diverse environments; however, a global view of the distribution, abundance, and diversity of Nitrospira species is still incomplete.

RESULTS

In this study, we retrieved 55 metagenome-assembled Nitrospira genomes (MAGs) from newly obtained and publicly available metagenomes. Combined with publicly available MAGs, this constitutes the largest Nitrospira genome database to date with 205 MAGs, representing 132 putative species, most without cultivated representatives. Mapping of metagenomic sequencing reads from various environments against this database enabled an analysis of the distribution and habitat preferences of Nitrospira species. Comammox Nitrospira's ecological success is evident as they outnumber and present higher species-level richness than canonical Nitrospira in all environments examined, except for marine and wastewaters samples. The type of environment governs Nitrospira species distribution, without large-scale biogeographical signal. We found that closely related Nitrospira species tend to occupy the same habitats, and that this phylogenetic signal in habitat preference is stronger for canonical Nitrospira species. Comammox Nitrospira eco-evolutionary history is more complex, with subclades achieving rapid niche divergence via horizontal transfer of genes, including the gene encoding hydroxylamine oxidoreductase, a key enzyme in nitrification.

CONCLUSIONS

Our study expands the genomic inventory of the Nitrospira genus, exposes the ecological success of complete ammonia oxidizers within a wide range of habitats, identifies the habitat preferences of (sub)lineages of canonical and comammox Nitrospira species, and proposes that horizontal transfer of genes involved in nitrification is linked to niche separation within a sublineage of comammox Nitrospira. Video Abstract.

Purifying selection enduringly acts on the sequence evolution of highly expressed proteins in Escherichia coli

The evolutionary speed of a protein sequence is constrained by its expression level, with highly expressed proteins evolving relatively slowly. This negative correlation between expression levels and evolutionary rates (known as the E–R anticorrelation) has already been widely observed in past macroevolution between species from bacteria to animals. However, it remains unclear whether this seemingly general law also governs recent evolution, including past and de novo, within a species. However, the advent of genomic sequencing and high-throughput phenotyping, particularly for bacteria, has revealed fundamental gaps between the 2 evolutionary processes and has provided empirical data opposing the possible underlying mechanisms which are widely believed. These conflicts raise questions about the generalization of the E–R anticorrelation and the relevance of plausible mechanisms. To explore the ubiquitous impact of expression levels on molecular evolution and test the relevance of the possible underlying mechanisms, we analyzed the genome sequences of 99 strains of Escherichia coli for evolution within species in nature. We also analyzed genomic mutations accumulated under laboratory conditions as a model of de novo evolution within species. Here, we show that E–R anticorrelation is significant in both past and de novo evolution within species in E. coli. Our data also confirmed ongoing purifying selection on highly expressed genes. Ongoing selection included codon-level purifying selection, supporting the relevance of the underlying mechanisms. However, the impact of codon-level purifying selection on the constraints in evolution within species might be smaller than previously expected from evolution between species.

Shigella viruses Sf22 and KRT47 require outer membrane protein C for infection

Viruses rely on hosts for their replication: thus, a critical step in the infection process is identifying a suitable host cell. Bacterial viruses, known as bacteriophages or phages, often use receptor binding proteins to discriminate between susceptible and non-susceptible hosts. By being able to evade predation, bacteria with modified or deleted receptor-encoding genes often undergo positive selection during growth in the presence of phage. Depending on the specific receptor(s) a phage uses, this may subsequently affect the bacteria's ability to form biofilms, its resistance to antibiotics, pathogenicity, or its phenotype in various environments. In this study, we characterize the interactions between two T4-like phages, Sf22 and KRT47, and their host receptor S. flexneri outer membrane protein C (OmpC). Results indicate that these phages use a variety of surface features on the protein, and that complete resistance most frequently occurs when hosts delete the ompC gene in full, encode premature stop codons to prevent OmpC synthesis, or eliminate specific regions encoding exterior loops.

Mutational Pattern Induced by 5-Fluorouracil and Oxaliplatin in the Gut Microbiome

Chemotherapeutic agents, such as 5-fluorouracil (5-FU) and oxaliplatin (Oxi), can not only kill the cancer cell but also influence the proliferation of gut microbiota; however, the interaction between these drugs and gut microbiota remains poorly understood. In this study, we developed a powerful framework for taxonomy composition and genomic variation analysis to investigate the mutagenesis effect and proliferation influence of chemotherapeutic agents, such as 5-FU and Oxi, on gut microbiota and the interaction between these drugs and gut microbiota during chemotherapy. Using the gut microbiome data, we detected 1.45 million variations among the chemotherapy groups and found the drugs significantly affected mutation signatures of gut microbiota. Oxi notably increased transversion rate, whereas 5-FU reduced the rate. Traits related to cell division and nutrient mobilization showed evidence of strong selection pressure from chemotherapeutic agents. In addition, drug-associated bacteriome shift patterns and functional alterations were found: the metabolism changes in the 5-FU group implied that gut microbiota could provide additional nicotinamide adenine dinucleotide (NAD⁺) to inhibit cancer cell autophagy; in the Oxi group, the ribosome and lysine biosynthesis genes were obviously enriched. Our study provides a blueprint for characterizing the role of microbes and drug–microbe interaction in the gut microbiota response to chemotherapy.

A genome-wide scan of wastewater E. coli for genes under positive selection: focusing on mechanisms of antibiotic resistance

Antibiotic resistance is a global health threat and consequently, there is a need to understand the mechanisms driving its emergence. Here, we hypothesize that genes and mutations under positive selection may contribute to antibiotic resistance. We explored wastewater E. coli, whose genomes are highly diverse. We subjected 92 genomes to a statistical analysis for positively selected genes. We obtained 75 genes under positive selection and explored their potential for antibiotic resistance. We found that eight genes have functions relating to antibiotic resistance, such as biofilm formation, membrane permeability, and bacterial persistence. Finally, we correlated the presence/absence of non-synonymous mutations in positively selected sites of the genes with a function in resistance against 20 most prescribed antibiotics. We identified mutations associated with antibiotic resistance in two genes: the porin ompC and the bacterial persistence gene hipA. These mutations are located at the surface of the proteins and may hence have a direct effect on structure and function. For hipA, we hypothesize that the mutations influence its interaction with hipB and that they enhance the capacity for dormancy as a strategy to evade antibiotics. Overall, genomic data and positive selection analyses uncover novel insights into mechanisms driving antibiotic resistance.

Comparative Genomics Applied to Systematically Assess Pathogenicity Potential in Shiga Toxin-Producing Escherichia coli O145:H28

Shiga toxin-producing Escherichia coli (STEC) O145:H28 can cause severe disease in humans and is a predominant serotype in STEC O145 environmental isolates. Here, comparative genomics was applied to a set of clinical and environmental strains to systematically evaluate the pathogenicity potential in environmental strains. While the core genes-based tree separated all O145:H28 strains from the non O145:H28 reference strains, it failed to segregate environmental strains from the clinical. In contrast, the accessory genes-based tree placed all clinical strains in the same clade regardless of their genotypes or serotypes, apart from the environmental strains. Loss-of-function mutations were common in the virulence genes examined, with a high frequency in genes related to adherence, autotransporters, and the type three secretion system. Distinct differences in pathogenicity islands LEE, OI-122, and OI-57, the acid fitness island, and the tellurite resistance island were detected between the O145:H28 and reference strains. A great amount of genetic variation was detected in O145:H28, which was mainly attributed to deletions, insertions, and gene acquisition at several chromosomal “hot spots”. Our study demonstrated a distinct virulence gene repertoire among the STEC O145:H28 strains originating from the same geographical region and revealed unforeseen contributions of loss-of-function mutations to virulence evolution and genetic diversification in STEC.

Mosaic Evolution of Beta-Barrel-Porin-Encoding Genes in Escherichia coli

Bacterial porin-encoding genes are often found under positive selection. Local recombination has also been identified in a few of them to facilitate bacterial rapid adaptation, although it remains unknown whether it is a common evolutionary mechanism for the porins or outer membrane proteins in Gram-negative bacteria. In this study, we investigated the beta-barrel (β-barrel) porin-encoding genes in Escherichia coli that were reported under positive Darwinian selection. Besides fhuA that was found with ingenic local recombination previously, we identified four other genes, i.e., lamB, ompA, ompC, and ompF, all showing the similar mosaic evolution patterns. Comparative analysis of the protein sequences disclosed a list of highly variable regions in each family, which are mostly located in the convex of extracellular loops and coinciding with the binding sites of bacteriophages. For each of the porin families, mosaic recombination leads to unique combinations of the variable regions with different sequence patterns, generating diverse protein groups. Structural modeling indicated a conserved global topology among the different porins, with the extracellular surface varying a lot due to individual or combinatorial variable regions. The conservation of global tertiary structure would ensure the channel activity, while the wide diversity of variable regions may represent selection to avoid the invasion of phages, antibiotics or immune surveillance factors. Our study identified multiple bacterial porin genes with mosaic evolution. We hypothesize that this could be generalized strategy for outer membrane proteins to both maintain normal life processes and evade the attack of unfavored factors rapidly. IMPORTANCE Microevolution studies can disclose more elaborate evolutionary mechanisms of genes, appearing especially important for genes with multifaceted function such as those encoding outer membrane proteins. However, in most cases, the gene is considered as a whole unit, and the evolutionary patterns are disclosed. Here, we report that multiple bacterial porin proteins follow mosaic evolution, with local ingenic recombination combined with spontaneous mutations based on positive Darwinian selection, and conservation for most structural regions. This could represent a common mechanism for bacterial outer membrane proteins. The variable regions within each porin family showed large coincidence with the binding sites of bacteriophages, antibiotics, and immune factors and therefore would represent effective targets for the development of new antibacterial agents or vaccines.

Comparative genomics and evolutionary analysis of Lactococcus garvieae isolated from human endocarditis

Lactococcus garvieae is a well-known pathogen of fish, but is rarely involved in infections in humans and other mammals. In humans, the main clinical manifestation of L. garvieae infections is endocarditis usually related to the ingestion of contaminated food, such as undercooked fish and shellfish. This study presents the first complete genomic sequence of a clinical L. garvieae strain isolated from a patient with endocarditis and its comparative analysis with other genomes. This human isolate contains a circular chromosome of 2 099 060 bp and one plasmid of 50 557 bp. In comparison with other fully sequenced L. garvieae strains, the chromosomal DNA of L. garvieae Lg-Granada carries a low proportion of insertion sequence elements and a higher number of putative prophages. Our results show that, in general, L. garvieae is a highly recombinogenic species with an open pangenome in which almost 30 % of its genome has undergone horizontal transfers. Within the genus Lactococcus, L. lactis is the main donor of genetic components to L. garvieae but, taking Lg-Granada as a representative, this bacterium tends to import more genes from Bacilli taxa than from other Lactococcus species.

Evolutionary Ecology of Natural Comammox Nitrospira Populations

Microbes commonly exist in diverse and complex communities where species interact, and their genomic repertoires evolve over time. Our understanding of species interaction and evolution has increased during the last decades, but most studies of evolutionary dynamics are based on single species in isolation or in experimental systems composed of few interacting species. Here, we use the microbial ecosystem found in groundwater-fed sand filter as a model to avoid this limitation. In these open systems, diverse microbial communities experience relatively stable conditions, and the coupling between chemical and biological processes is generally well defined. Metagenomic analysis of 12 sand filters communities revealed systematic co-occurrence of at least five comammox Nitrospira species, likely promoted by low ammonium concentrations. These Nitrospira species showed intrapopulation sequence diversity, although possible clonal expansion was detected in a few abundant local comammox populations. Nitrospira species showed low homologous recombination and strong purifying selection, the latter process being especially strong in genes essential in energy metabolism. Positive selection was detected for genes related to resistance to foreign DNA and phages. We found that, compared to other habitats, groundwater-fed sand filters impose strong purifying selection and low recombination on comammox Nitrospira populations. These results suggest that evolutionary processes are more affected by habitat type than by species identity. Together, this study improves our understanding of species interaction and evolution in complex microbial communities and sheds light on the environmental dependency of evolutionary processes. IMPORTANCE Microbial species interact with each other and their environment (ecological processes) and undergo changes in their genomic repertoire over time (evolutionary processes). How these two classes of processes interact is largely unknown, especially for complex communities, as most studies of microbial evolutionary dynamics consider single species in isolation or a few interacting species in simplified experimental systems. In this study, these limitations are circumvented by examining the microbial communities found in stable and well-described groundwater-fed sand filters. Combining metagenomics and strain-level analyses, we identified the microbial interactions and evolutionary processes affecting comammox Nitrospira, a recently discovered bacterial type capable of performing the whole nitrification process. We found that abundant and co-occurrent Nitrospira populations in groundwater-fed sand filters are characterized by low recombination and strong purifying selection. In addition, by comparing these observations with those obtained from Nitrospira species inhabiting other environments, we revealed that evolutionary processes are more affected by habitat type than by species identity.

Evidence of episodic positive selection in Corynebacterium diphtheriae complex of species and its implementations in identification of drug and vaccine targets

BACKGROUND

Within the pathogenic bacterial species Corynebacterium genus, six species that can produce diphtheria toxin (C. belfantii, C. diphtheriae, C. pseudotuberculosis, C. rouxii, C. silvaticum and C. ulcerans) form a clade referred to as the C. diphtheria complex. These species have been found in humans and other animals, causing diphtheria or other diseases. Here we show the results of a genome scale analysis to identify positive selection in protein-coding genes that may have resulted in the adaptations of these species to their ecological niches and suggest drug and vaccine targets.

METHODS

Forty genomes were sampled to represent species, subspecies or biovars of Corynebacterium. Ten phylogenetic groups were tested for positive selection using the PosiGene pipeline, including species and biovars from the C. diphtheria complex. The detected genes were tested for recombination and had their sequences alignments and homology manually examined. The final genes were investigated for their function and a probable role as vaccine or drug targets.

RESULTS

Nineteen genes were detected in the species C. diphtheriae (two), C. pseudotuberculosis (10), C. rouxii (one), and C. ulcerans (six). Those were found to be involved in defense, translation, energy production, and transport and in the metabolism of carbohydrates, amino acids, nucleotides, and coenzymes. Fourteen were identified as essential genes, and six as virulence factors. Thirteen from the 19 genes were identified as potential drug targets and four as potential vaccine candidates. These genes could be important in the prevention and treatment of the diseases caused by these bacteria.

Assessment of selection pressure exerted on genes from complete pangenomes helps to improve the accuracy in the prediction of new genes

Bacterial genomes are massively sequenced, and they provide valuable data to better know the complete set of genes of a species. The analysis of thousands of bacterial strains can identify both shared genes and those appearing only in the pathogenic ones. Current computational gene finders facilitate this task but often miss some existing genes. However, the present availability of different genomes from the same species is useful to estimate the selective pressure applied on genes of complete pangenomes. It may assist in evaluating gene predictions either by checking the certainty of a new gene or annotating it as a gene under positive selection. Here, we estimated the selective pressure of 19 271 genes that are part of the pangenome of the human opportunistic pathogen Acinetobacter baumannii and found that most genes in this bacterium are subject to negative selection. However, 23% of them showed values compatible with positive selection. These latter were mainly uncharacterized proteins or genes required to evade the host defence system including genes related to resistance and virulence whose changes may be favoured to acquire new functions. Finally, we evaluated the utility of measuring selection pressure in the detection of sequencing errors and the validation of gene prediction.

Characterization of the first vaginal Lactobacillus crispatus genomes isolated in Brazil

BACKGROUND

Lactobacillus crispatus is the dominant species in the vaginal microbiota associated with health and considered a homeostasis biomarker. Interestingly, some strains are even used as probiotics. However, the genetic mechanisms of L. crispatus involved in the control of the vaginal microbiome and protection against bacterial vaginosis (BV) are not entirely known. To further investigate these mechanisms, we sequenced and characterized the first four L. crispatus genomes from vaginal samples from Brazilian women and used genome-wide association study (GWAS) and comparative analyses to identify genetic mechanisms involved in healthy or BV conditions and selective pressures acting in the vaginal microbiome.

METHODS

The four genomes were sequenced, assembled using ten different strategies and automatically annotated. The functional characterization was performed by bioinformatics tools comparing with known probiotic strains. Moreover, it was selected one representative strain (L. crispatus CRI4) for in vitro detection of phages by electron microscopy. Evolutionary analysis, including phylogeny, GWAS and positive selection were performed using 46 public genomes strains representing health and BV conditions.

RESULTS

Genes involved in probiotic effects such as lactic acid production, hydrogen peroxide, bacteriocins, and adhesin were identified. Three hemolysins and putrescine production were predicted, although these features are also present in other probiotic strains. The four genomes presented no plasmids, but 14 known families insertion sequences and several prophages were detected. However, none of the mobile genetic elements contained antimicrobial resistance genes. The genomes harbor a CRISPR-Cas subtype II-A system that is probably inactivated due to fragmentation of the genes csn2 and cas9. No genomic feature was associated with a health condition, perhaps due to its multifactorial characteristic. Five genes were identified as under positive selection, but the selective pressure remains to be discovered. In conclusion, the Brazilian strains investigated in this study present potential protective properties, although in vitro and in vivo studies are required to confirm their efficacy and safety to be considered for human use.

Diversification of OmpA and OmpF of Yersinia ruckeri is independent of the underlying species phylogeny and evidence of virulence-related selection

Yersinia ruckeri is the causative agent of enteric redmouth disease (ERM) which causes economically significant losses in farmed salmonids, especially Atlantic salmon (Salmo salar) and rainbow trout (Oncorhynchus mykiss, Walbaum). However, very little is known about the genetic relationships of disease-causing isolates in these two host species or about factors responsible for disease. Phylogenetic analyses of 16 representative isolates based on the nucleotide sequences of 19 housekeeping genes suggests that pathogenic Atlantic salmon and rainbow trout isolates represent distinct host-specific lineages. However, the apparent phylogenies of certain isolates has been influenced by horizontal gene transfer and recombinational exchange. Splits decomposition analysis demonstrated a net-like phylogeny based on the housekeeping genes, characteristic of recombination. Comparative analysis of the distribution of individual housekeeping gene alleles across the isolates demonstrated evidence of genomic mosaicism and recombinational exchange involving certain Atlantic salmon and rainbow trout isolates. Comparative nucleotide sequence analysis of the key outer membrane protein genes ompA and ompF revealed that the corresponding gene trees were both non-congruent with respect to the housekeeping gene phylogenies providing evidence that horizontal gene transfer has influenced the evolution of both these surface protein-encoding genes. Analysis of inferred amino acid sequence variation in OmpA identified a single variant, OmpA.1, that was present in serotype O1 and O8 isolates representing typical pathogenic strains in rainbow trout and Atlantic salmon, respectively. In particular, the sequence of surface-exposed loop 3 differed by seven amino acids to that of other Y. ruckeri isolates. These findings suggest that positive selection has likely influenced the presence of OmpA.1 in these isolates and that loop 3 may play an important role in virulence. Amino acid sequence variation of OmpF was greater than that of OmpA and was similarly restricted mainly to the surface-exposed loops. Two OmpF variants, OmpF.1 and OmpF.2, were associated with pathogenic rainbow trout and Atlantic salmon isolates, respectively. These OmpF proteins had very similar amino acid sequences suggesting that positive evolutionary pressure has also favoured the selection of these variants in pathogenic strains infecting both species.

Insights into the Metabolism and Evolution of the Genus Acidiphilium, a Typical Acidophile in Acid Mine Drainage

Here, we report three new Acidiphilium genomes, reclassified existing Acidiphilium species, and performed the first comparative genomic analysis on Acidiphilium in an attempt to address the metabolic potential, ecological functions, and evolutionary history of the genus Acidiphilium In the genomes of Acidiphilium, we found an abundant repertoire of horizontally transferred genes (HTGs) contributing to environmental adaption and metabolic expansion, including genes conferring photosynthesis (puf, puh), CO2 assimilation (rbc), capacity for methane metabolism (mmo, mdh, frm), nitrogen source utilization (nar, cyn, hmp), sulfur compound utilization (sox, psr, sqr), and multiple metal and osmotic stress resistance capacities (czc, cop, ect). Additionally, the predicted donors of horizontal gene transfer were present in a cooccurrence network of Acidiphilium Genome-scale positive selection analysis revealed that 15 genes contained adaptive mutations, most of which were multifunctional and played critical roles in the survival of extreme conditions. We proposed that Acidiphilium originated in mild conditions and adapted to extreme environments such as acidic mineral sites after the acquisition of many essential functions.IMPORTANCE Extremophiles, organisms that thrive in extreme environments, are key models for research on biological adaption. They can provide hints for the origin and evolution of life, as well as improve the understanding of biogeochemical cycling of elements. Extremely acidophilic bacteria such as Acidiphilium are widespread in acid mine drainage (AMD) systems, but the metabolic potential, ecological functions, and evolutionary history of this genus are still ambiguous. Here, we sequenced the genomes of three new Acidiphilium strains and performed comparative genomic analysis on this extremely acidophilic bacterial genus. We found in the genomes of Acidiphilium an abundant repertoire of horizontally transferred genes (HTGs) contributing to environmental adaption and metabolic ability expansion, as indicated by phylogenetic reconstruction and gene context comparison. This study has advanced our understanding of microbial evolution and biogeochemical cycling in extreme niches.

Characterization of Burkholderia cepacia Complex Core Genome and the Underlying Recombination and Positive Selection

Recombination and positive selection are two key factors that play a vital role in pathogenic microorganisms’ population adaptation and diversification. The Burkholderia cepacia complex (Bcc) represents bacterial species with high similarity, which can cause severe infections among cases suffering from the chronic granulomatous disorder and cystic fibrosis (CF). At present, no genome-wide study has been carried out focusing on investigating the core genome of Bcc associated with the two evolutionary forces. The general characteristics of the core genome of Bcc species remain scarce as well. In this study, we explored the core orthologous genes of 116 Bcc strains using comparative genomic analysis and studied the two adaptive evolutionary forces: recombination and positive selection. We estimated 1005 orthogroups consisting entirely of single copy genes. These single copy orthologous genes in some Cluster of Orthologous Groups (COG) categories showed significant differences in the comparison of several evolutionary properties, and the encoding proteins were relatively simple and compact. Our findings showed that 5.8% of the core orthologous genes strongly supported recombination; in the meantime, 1.1% supported positive selection. We found that genes involved in protein synthesis as well as material transport and metabolism are favored by selection pressure. More importantly, homologous recombination contributed more genetic variation to a large number of genes and largely maintained the genetic cohesion in Bcc. This high level of recombination between Bcc species blurs their taxonomic boundaries, which leads Bcc species to be difficult or impossible to distinguish phenotypically and genotypically.

Comparative genomics reveals different population structures associated with host and geographic origin in antimicrobial-resistant Salmonella enterica

Genetic variation in a pathogen, including the causative agent of salmonellosis, Salmonella enterica, can occur as a result of eco-evolutionary forces triggered by dissimilarities of ecological niches. Here, we applied comparative genomics to study 90 antimicrobial resistant (AMR) S. enterica isolates from bovine and human hosts in New York and Washington states to understand host- and geographic-associated population structure. Results revealed distinct presence/absence profiles of functional genes and pseudogenes (e.g., virulence genes) associated with bovine and human isolates. Notably, bovine isolates contained significantly more transposase genes but fewer transposase pseudogenes than human isolates, suggesting the occurrence of large-scale transposition in genomes of bovine and human isolates at different times. The high correlation between transposase genes and AMR genes, as well as plasmid replicons, highlights the potential role of horizontally transferred transposons in promoting adaptation to antibiotics. By contrast, a number of potentially geographic-associated single-nucleotide polymorphisms (SNPs), rather than geographic-associated genes, were identified. Interestingly, 38% of these SNPs were in genes annotated as cell surface protein-encoding genes, including some essential for antibiotic resistance and host colonization. Overall, different evolutionary forces and limited recent inter-population transmission appear to shape AMR S. enterica population structure in different hosts and geographic origins.

RHS-elements function as type II toxin-antitoxin modules that regulate intra-macrophage replication of Salmonella Typhimurium

RHS elements are components of conserved toxin-delivery systems, wide-spread within the bacterial kingdom and some of the most positively selected genes known. However, very little is known about how Rhs toxins affect bacterial biology. Salmonella Typhimurium contains a full-length rhs gene and an adjacent orphan rhs gene, which lacks the conserved delivery part of the Rhs protein. Here we show that, in addition to the conventional delivery, Rhs toxin-antitoxin pairs encode for functional type-II toxin-antitoxin (TA) loci that regulate S. Typhimurium proliferation within macrophages. Mutant S. Typhimurium cells lacking both Rhs toxins proliferate 2-times better within macrophages, mainly because of an increased growth rate. Thus, in addition to providing strong positive selection for the rhs loci under conditions when there is little or no toxin delivery, internal expression of the toxin-antitoxin system regulates growth in the stressful environment found inside macrophages.

Bacteria that reside and multiply inside of phagocytic cells are hard to treat with common antibiotics, partly because subpopulations of bacteria are non-growing. Very little is known about how bacteria regulate their growth in the phagocytic vesicle. We show that RHS elements, previously known to function as mobilizable toxins that inhibit growth of neighboring bacteria, also function as internally expressed toxin-antitoxin systems that regulate Salmonella Typhimurium growth in macrophages. RHS elements were discovered more than 30 years ago, but their role in biology has long remained unclear even though they are some of the most positively selected genes known. Our results suggest an explanation to why rhs genes are under such strong positive selection in addition to suggesting a novel function for these toxins in regulating bacterial growth.

Standardized phylogenetic and molecular evolutionary analysis applied to species across the microbial tree of life

There is growing interest in reconstructing phylogenies from the copious amounts of genome sequencing projects that target related viral, bacterial or eukaryotic organisms. To facilitate the construction of standardized and robust phylogenies for disparate types of projects, we have developed a complete bioinformatic workflow, with a web-based component to perform phylogenetic and molecular evolutionary (PhaME) analysis from sequencing reads, draft assemblies or completed genomes of closely related organisms. Furthermore, the ability to incorporate raw data, including some metagenomic samples containing a target organism (e.g. from clinical samples with suspected infectious agents), shows promise for the rapid phylogenetic characterization of organisms within complex samples without the need for prior assembly.

Residual Variation Intolerance Score Detects Loci Under Selection in Neuroinvasive Listeria monocytogenes

Listeria monocytogenes is a Gram-positive bacterium that can be found in a broad range of environments, including soil, food, animals, and humans. L. monocytogenes can cause a foodborne disease manifesting as sepsis and meningo-encephalitis. To evaluate signals of selection within the core genome of neuroinvasive L. monocytogenes strains, we sequenced 122 L. monocytogenes strains from cerebrospinal fluid (CSF) of Dutch meningitis patients and performed a genome-wide analysis using Tajima’s D and ω (dN/dS). We also evaluated the residual variation intolerance score (RVIS), a computationally less demanding methodology, to identify loci under selection. Results show that the large genetic distance between the listerial lineages influences the Tajima’s D and ω (dN/dS) outcome. Within genetic lineages we detected signals of selection in 6 of 2327 loci (<1%), which were replicated in an external cohort of 105 listerial CSF isolates from France. Functions of identified loci under selection were within metabolism pathways (lmo2476, encoding aldose 1-epimerase), putative antimicrobial resistance mechanisms (lmo1855, encoding PBPD3), and virulence factors (lmo0549, internalin-like protein; lmo1482, encoding comEC). RVIS over the two genetic lineages showed signals of selection in internalin-like proteins loci potentially involved in pathogen-host interaction (lmo0549, lmo0610, and lmo1290). Our results show that RVIS can be used to detect bacterial loci under selection.

Adaptive Evolution of Sphingobium hydrophobicum C1T in Electronic Waste Contaminated River Sediment

Electronic waste (e-waste) has caused a severe worldwide pollution problem. Despite increasing isolation of degradative microorganisms from e-waste contaminated environments, the mechanisms underlying their adaptive evolution in such habitats remain unclear. Sphingomonads generally have xenobiotic-degrading ability and may play important roles in bioremediation. Sphingobium hydrophobicum C1^T, characterized with superior cell surface hydrophobicity, was recently isolated from e-waste contaminated river sediment. To dissect the mechanisms driving its adaptive evolution, we evaluated its stress resistance, sequenced its genome and performed comparative genomic analysis with 19 other Sphingobium strains. Strain C1^T can feed on several kinds of e-waste-derived xenobiotics, exhibits a great resistance to heavy metals and possesses a high colonization ability. It harbors abundant genes involved in environmental adaptation, some of which are intrinsic prior to experiencing e-waste contamination. The extensive genomic variations between strain C1^T and other Sphingobium strains, numerous C1^T-unique genes, massive mobile elements and frequent genome rearrangements reflect a high genome plasticity. Positive selection, gene duplication, and especially horizontal gene transfer drive the adaptive evolution of strain C1^T. Moreover, presence of type IV secretion systems may allow strain C1^T to be a source of beneficial genes for surrounding microorganisms. This study provides new insights into the adaptive evolution of sphingomonads, and potentially guides bioremediation strategies.

Serotype-specific evolutionary patterns of antimicrobial-resistant Salmonella enterica

Background

The emergence of antimicrobial-resistant (AMR) strains of the important human and animal pathogen Salmonella enterica poses a growing threat to public health. Here, we studied the genome-wide evolution of 90 S. enterica AMR isolates, representing one host adapted serotype (S. Dublin) and two broad host range serotypes (S. Newport and S. Typhimurium).

Results

AMR S. Typhimurium had a large effective population size, a large and diverse genome, AMR profiles with high diversity, and frequent positive selection and homologous recombination. AMR S. Newport showed a relatively low level of diversity and a relatively clonal population structure. AMR S. Dublin showed evidence for a recent population bottleneck, and the genomes were characterized by a larger number of genes and gene ontology terms specifically absent from this serotype and a significantly higher number of pseudogenes as compared to other two serotypes. Approximately 50% of accessory genes, including specific AMR and putative prophage genes, were significantly over- or under-represented in a given serotype. Approximately 65% of the core genes showed phylogenetic clustering by serotype, including the AMR gene aac (6′)-Iaa. While cell surface proteins were shown to be the main target of positive selection, some proteins with possible functions in AMR and virulence also showed evidence for positive selection. Homologous recombination mainly acted on prophage-associated proteins.

Conclusions

Our data indicates a strong association between genome content of S. enterica and serotype. Evolutionary patterns observed in S. Typhimurium are consistent with multiple emergence events of AMR strains and/or ecological success of this serotype in different hosts or habitats. Evolutionary patterns of S. Newport suggested that antimicrobial resistance emerged in one single lineage, Lineage IIC. A recent population bottleneck and genome decay observed in AMR S. Dublin are congruent with its narrow host range. Finally, our results suggest the potentially important role of positive selection in the evolution of antimicrobial resistance, host adaptation and serotype diversification in S. enterica.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1457-5) contains supplementary material, which is available to authorized users.

Identification of positively selected genes in human pathogenic treponemes: Syphilis-, yaws-, and bejel-causing strains differ in sets of genes showing adaptive evolution

BACKGROUND

Pathogenic treponemes related to Treponema pallidum are both human (causing syphilis, yaws, bejel) and animal pathogens (infections of primates, venereal spirochetosis in rabbits). A set of 11 treponemal genome sequences including those of five Treponema pallidum ssp. pallidum (TPA) strains (Nichols, DAL-1, Mexico A, SS14, Chicago), four T. p. ssp. pertenue (TPE) strains (CDC-2, Gauthier, Samoa D, Fribourg-Blanc), one T. p. ssp. endemicum (TEN) strain (Bosnia A) and one strain (Cuniculi A) of Treponema paraluisleporidarum ecovar Cuniculus (TPeC) were tested for the presence of positively selected genes.

METHODOLOGY/PRINCIPAL FINDINGS

A total of 1068 orthologous genes annotated in all 11 genomes were tested for the presence of positively selected genes using both site and branch-site models with CODEML (PAML package). Subsequent analyses with sequences obtained from 62 treponemal draft genomes were used for the identification of positively selected amino acid positions. Synthetic biotinylated peptides were designed to cover positively selected protein regions and these peptides were tested for reactivity with the patient's syphilis sera. Altogether, 22 positively selected genes were identified in the TP genomes and TPA sets of positively selected genes differed from TPE genes. While genetic variability among TPA strains was predominantly present in a number of genetic loci, genetic variability within TPE and TEN strains was distributed more equally along the chromosome. Several syphilitic sera were shown to react with some peptides derived from the protein sequences evolving under positive selection.

CONCLUSIONS/SIGNIFICANCE

The syphilis-, yaws-, and bejel-causing strains differed relative to sets of positively selected genes. Most of the positively selected chromosomal loci were identified among the TPA treponemes. The local accumulation of genetic variability suggests that the diversification of TPA strains took place predominantly in a limited number of genomic regions compared to the more dispersed genetic diversity differentiating TPE and TEN strains. The identification of positively selected sites in tpr genes and genes encoding outer membrane proteins suggests their role during infection of human and animal hosts. The driving force for adaptive evolution at these loci thus appears to be the host immune response as supported by observed reactivity of syphilitic sera with some peptides derived from protein sequences showing adaptive evolution.

Genome-wide analyses reveal genes subject to positive selection in Toxoplasma gondii

Toxoplasma gondii is an important protozoan pathogen that infects many wild and domestic animals and causes infections in immunocompromised humans. However, there has been little investigation of the molecular evolutionary trajectories of this pathogenic protozoa using comparative genomics data. Here, we employed a comparative evolutionary genomics approach to identify genes that are under site- and lineage-specific positive selection in nine strains of T. gondii, including two closely related species, Neospora caninum and Hammondia hammondi. Based on the analyses of five coccidian core genomes, 4.5% of the 5788 core genome genes showed strong signals for positive selection in the site model. In addition, the branch-site model analyses in the nine T. gondii core genomes indicated that 2 to 20 genes underwent significant positive selection along each lineage leading to T. gondii strains. Many of the protein products encoded by the positively selected genes are secretory or surface proteins that have previously been implicated in host pathogenesis. The adaptive changes in these positively selected genes might be related to dynamic interactions between the host immune systems and might play a crucial role in the infection and pathogenic processes of T. gondii.

Genomic analysis of Leptospira interrogans serovar Paidjan and Dadas isolates from carrier dogs and comparative genomic analysis to detect genes under positive selection

BACKGROUND

Leptospirosis is an emerging infectious disease worldwide that can cause high morbidity and mortality rates in humans and animals. The causative spirochetes have reservoirs in mammalian hosts, but there has been limited analysis of the genomes of isolates recovered from animals. The aims of this study were to characterize genomic features of two Leptospira interrogans strains recently isolated from asymptomatic dogs in Thailand (strains CUDO5 and CDUO8), and to perform comparative genome analyses with other strains. Molecular adaptive evolution in L. interrogans as signaled by positive selection also was analyzed.

RESULTS

Whole genome sequence analysis revealed that strains CUDO5 and CUDO8 had genome sizes of approximately 4.9 Mbp with 35.1% GC contents. Using monoclonal antibodies, strains CUDO5 and CUDO8 were identified as serovars Paidjan and Dadas, respectively. These strains harbored genes known to be associated with acute and chronic infections. Using Single Nucleotide Polymorphisms phylogeny (SNPs) with 97 L. interrogans strains, CUDO5 and CUDO8 had closest genetic relatedness with each other. Nevertheless, the serovar determinant region (rfb locus) showed variations in the genes encoding sugar biosynthesis. Amongst 13 representative L. interrogans strains examined for molecular adaptive evolution through positive selection under the site-model of Phylogenetic Analysis of Maximum Likelihood, genes responsible for iron acquisition (tlyA and hbpA), motility (fliN2, flgK, and flhB) and thermal adaptation (lpxD1) were under increased selective pressure.

CONCLUSIONS

L. interrogans serovar Paidjan strain CUDO5 and serovar Dadas strain CUDO8 had close genetic relatedness as analyzed by SNPs phylogeny. They contained genes with established roles in acute and chronic leptospirosis. The rfb locus in both serovars showed gene variation associated with sugar biosynthesis. Positive selection analysis indicated that genes encoding factors involved in motility, temperature adaptation, and iron acquisition were under strong positive selection in L. interrogans. These may be associated with adaptation in the early stages of infection.

Molecular Characterization of Carbapenem Resistant Klebsiella pneumoniae and Klebsiella quasipneumoniae Isolated from Lebanon

Klebsiella pneumoniae is a Gram-negative organism and a major public health threat. In this study, we used whole-genome sequences to characterize 32 carbapenem-resistant K. pneumoniae (CRKP) and two carbapenem-resistant K. quasipneumoniae (CRKQ). Antimicrobial resistance was assessed using disk diffusion and E-test, while virulence was assessed in silico. The capsule type was determined by sequencing the wzi gene. The plasmid diversity was assessed by PCR-based replicon typing to detect the plasmid incompatibility (Inc) groups. The genetic relatedness was determined by multilocus sequence typing, pan-genome, and recombination analysis. All of the isolates were resistant to ertapenem together with imipenem and/or meropenem. Phenotypic resistance was due to bla_OXA-48,bla_NDM-1, bla_NDM-7, or the coupling of ESBLs and outer membrane porin modifications. This is the first comprehensive study reporting on the WGS of CRKP and the first detection of CRKQ in the region. The presence and dissemination of CRKP and CRKQ, with some additionally having characteristics of hypervirulent clones such as the hypermucoviscous phenotype and the capsular type K2, are particularly concerning. Additionally, mining the completely sequenced K. pneumoniae genomes revealed the key roles of mobile genetic elements in the spread of antibiotic resistance and in understanding the epidemiology of these clinically significant pathogens.

Among older adults, age-related changes in the stool microbiome differ by HIV-1 serostatus

Background

HIV-1 infection and physiological aging are independently linked to elevated systemic inflammation and changes in enteric microbial communities (dysbiosis). However, knowledge of the direct effect of HIV infection on the aging microbiome and potential links to systemic inflammation is lacking.

Methods

In a cross-sectional study of older people living with HIV (PLWH) (median age 61.5 years, N = 14) and uninfected controls (median 58 years, n = 22) we compared stool microbiota, levels of microbial metabolites (short-chain fatty acid levels, SCFA) and systemic inflammatory biomarkers by HIV serostatus and age.

Findings

HIV and age were independently associated with distinct changes in the stool microbiome. For example, abundances of Enterobacter and Paraprevotella were higher and Eggerthella and Roseburia lower among PLWH compared to uninfected controls. Age-related microbiome changes also differed by HIV serostatus. Some bacteria with inflammatory potential (e.g. Escherichia) increased with age among PLWH, but not controls. Stool SCFA levels were similar between the two groups yet patterns of associations between individual microbial taxa and SCFA levels differed. Abundance of various genera including Escherichia and Bifidobacterium positively associated with inflammatory biomarkers (e.g. soluble Tumor Necrosis Factor Receptors) among PLWH, but not among controls.

Interpretation

The age effect on the gut microbiome and associations between microbiota and microbial metabolites or systemic inflammation differed based on HIV serostatus, raising important implications for the impact of therapeutic interventions, dependent on HIV serostatus or age.

Rapidly evolving changes and gene loss associated with host switching in Corynebacterium pseudotuberculosis

Phylogenomics and genome scale positive selection analyses were performed on 29 Corynebacterium pseudotuberculosis genomes that were isolated from different hosts, including representatives of the Ovis and Equi biovars. A total of 27 genes were identified as undergoing adaptive changes. An analysis of the clades within this species and these biovars, the genes specific to each branch, and the genes responding to selective pressure show clear differences, indicating that adaptation and specialization is occurring in different clades. These changes are often correlated with the isolation host but could indicate responses to some undetermined factor in the respective niches. The fact that some of these more-rapidly evolving genes have homology to known virulence factors, antimicrobial resistance genes and drug targets shows that this type of analysis could be used to identify novel targets, and that these could be used as a way to control this pathogen.

Evolution and Sequence Diversity of FhuA in Salmonella and Escherichia

The fhuACDB operon, present in a number of Enterobacteriaceae, encodes components essential for the uptake of ferric hydroxamate type siderophores. FhuA acts not only as a transporter for physiologically important chelated ferric iron but also as a receptor for various bacteriophages, toxins, and antibiotics, which are pathogenic to bacterial cells. In this research, fhuA gene distribution and sequence diversity were investigated in Enterobacteriaceae, especially Salmonella and Escherichia Comparative sequence analysis resulted in a fhuA phylogenetic tree that did not match the expected phylogeny of species or trees of the fhuCDB genes. The fhuA sequences showed a unique mosaic clustering pattern. On the other hand, the gene sequences showed high conservation for strains from the same serovar or serotype. In total, six clusters were identified from FhuA proteins in Salmonella and Escherichia, among which typical peptide fragment variations could be defined. Six fragmental insertions/deletions and two substitution fragments were discovered, for which the combination of polymorphism patterns could well classify the different clusters. Structural modeling demonstrated that all the six featured insertions/deletions and one substitution fragment are located at the apexes of the long loops present as part of the FhuA external pocket. These frequently mutated regions are likely under high selection pressure, with bacterial strains balancing escape from phage infection or toxin/antibiotics attack via fhuA gene mutations while maintaining the siderophore uptake activity essential for bacterial survival. The unusual fhuA clustering suggests that high-frequency exchange of fhuA genes has occurred between enterobacterial strains after distinctive species were established.

Multi-scale characterization of symbiont diversity in the pea aphid complex through metagenomic approaches

BACKGROUND

Most metazoans are involved in durable relationships with microbes which can take several forms, from mutualism to parasitism. The advances of NGS technologies and bioinformatics tools have opened opportunities to shed light on the diversity of microbial communities and to give some insights into the functions they perform in a broad array of hosts. The pea aphid is a model system for the study of insect-bacteria symbiosis. It is organized in a complex of biotypes, each adapted to specific host plants. It harbors both an obligatory symbiont supplying key nutrients and several facultative symbionts bringing additional functions to the host, such as protection against biotic and abiotic stresses. However, little is known on how the symbiont genomic diversity is structured at different scales: across host biotypes, among individuals of the same biotype, or within individual aphids, which limits our understanding on how these multi-partner symbioses evolve and interact.

RESULTS

We present a framework well adapted to the study of genomic diversity and evolutionary dynamics of the pea aphid holobiont from metagenomic read sets, based on mapping to reference genomes and whole genome variant calling. Our results revealed that the pea aphid microbiota is dominated by a few heritable bacterial symbionts reported in earlier works, with no discovery of new microbial associates. However, we detected a large and heterogeneous genotypic diversity associated with the different symbionts of the pea aphid. Partitioning analysis showed that this fine resolution diversity is distributed across the three considered scales. Phylogenetic analyses highlighted frequent horizontal transfers of facultative symbionts between host lineages, indicative of flexible associations between the pea aphid and its microbiota. However, the evolutionary dynamics of symbiotic associations strongly varied depending on the symbiont, reflecting different histories and possible constraints. In addition, at the intra-host scale, we showed that different symbiont strains may coexist inside the same aphid host.

CONCLUSIONS

We present a methodological framework for the detailed analysis of NGS data from microbial communities of moderate complexity and gave major insights into the extent of diversity in pea aphid-symbiont associations and the range of evolutionary trajectories they could take.

Genes under positive selection in the core genome of pathogenic Bacillus cereus group members

In this comparative genomics study our aim was to unravel genes under positive selection in the core genome of the Bacillus cereus group. Indeed, the members of this group share close genetic relationships but display a rather large phenotypic and ecological diversity, providing a unique opportunity for studying how genomic changes reflect ecological adaptation during the divergence of a bacterial group. For this purpose, we screened ten completely sequenced genomes of four pathogenic Bacillus species, finding that 254 out of 3093 genes have codon sites with d_N/d_S (ω) values above one. These results remained unchanged after having disentangled the confounding effects of recombination and selection signature in a Bayesian framework. The presumably adaptive nucleotide polymorphisms are distributed over a wide range of biological functions, such as antibiotic resistance, DNA repair, nutrient uptake, metabolism, cell wall assembly and spore structure. Our results indicate that adaptation to animal hosts, whether as pathogens, saprophytes or symbionts, is the major driving force in the evolution of the Bacillus cereus group. Future work should seek to understand the evolutionary dynamics of both core and accessory genes in an integrative framework to ultimately unravel the key networks involved in host adaptation.

PosiGene: automated and easy-to-use pipeline for genome-wide detection of positively selected genes

Many comparative genomics studies aim to find the genetic basis of species-specific phenotypic traits. A prevailing strategy is to search genome-wide for genes that evolved under positive selection based on the non-synonymous to synonymous substitution ratio. However, incongruent results largely due to high false positive rates indicate the need for standardization of quality criteria and software tools. Main challenges are the ortholog and isoform assignment, the high sensitivity of the statistical models to alignment errors and the imperative to parallelize large parts of the software. We developed the software tool PosiGene that (i) detects positively selected genes (PSGs) on genome-scale, (ii) allows analysis of specific evolutionary branches, (iii) can be used in arbitrary species contexts and (iv) offers visualization of the results for further manual validation and biological interpretation. We exemplify PosiGene's performance using simulated and real data. In the simulated data approach, we determined a false positive rate <1%. With real data, we found that 68.4% of the PSGs detected by PosiGene, were shared by at least one previous study that used the same set of species. PosiGene is a user-friendly, reliable tool for reproducible genome-wide identification of PSGs and freely available at https://github.com/gengit/PosiGene.

Genomic Comparisons Reveal Microevolutionary Differences in Mycobacterium abscessus Subspecies

Mycobacterium abscessus, a rapid-growing non-tuberculous mycobacterium, has been the cause of sporadic and outbreak infections world-wide. The subspecies in M. abscessus complex (M. abscessus, M. massiliense, and M. bolletii) are associated with different biologic and pathogenic characteristics and are known to be among the most frequently isolated opportunistic pathogens from clinical material. To date, the evolutionary forces that could have contributed to these biological and clinical differences are still unclear. We compared genome data from 243 M. abscessus strains downloaded from the NCBI ftp Refseq database to understand how the microevolutionary processes of homologous recombination and positive selection influenced the diversification of the M. abscessus complex at the subspecies level. The three subspecies are clearly separated in the Minimum Spanning Tree. Their MUMi-based genomic distances support the separation of M. massiliense and M. bolletii into two subspecies. Maximum Likelihood analysis through dN/dS (the ratio of number of non-synonymous substitutions per non-synonymous site, to the number of synonymous substitutions per synonymous site) identified distinct genes in each subspecies that could have been affected by positive selection during evolution. The results of genome-wide alignment based on concatenated locally-collinear blocks suggest that (a) recombination has affected the M. abscessus complex more than mutation and positive selection; (b) recombination occurred more frequently in M. massiliense than in the other two subspecies; and (c) the recombined segments in the three subspecies have come from different intra-species and inter-species origins. The results lead to the identification of possible gene sets that could have been responsible for the subspecies-specific features and suggest independent evolution among the three subspecies, with recombination playing a more significant role than positive selection in the diversification among members in this complex.

Genome-Wide Analyses Reveal Genes Subject to Positive Selection in Pasteurella multocida

Pasteurella multocida, a Gram-negative opportunistic pathogen, has led to a broad range of diseases in mammals and birds, including fowl cholera in poultry, pneumonia and atrophic rhinitis in swine and rabbit, hemorrhagic septicemia in cattle, and bite infections in humans. In order to better interpret the genetic diversity and adaptation evolution of this pathogen, seven genomes of P. multocida strains isolated from fowls, rabbit and pigs were determined by using high-throughput sequencing approach. Together with publicly available P. multocida genomes, evolutionary features were systematically analyzed in this study. Clustering of 70,565 protein-coding genes showed that the pangenome of 33 P. multocida strains was composed of 1,602 core genes, 1,364 dispensable genes, and 1,070 strain-specific genes. Of these, we identified a full spectrum of genes related to virulence factors and revealed genetic diversity of these potential virulence markers across P. multocida strains, e.g., bcbAB, fcbC, lipA, bexDCA, ctrCD, lgtA, lgtC, lic2A involved in biogenesis of surface polysaccharides, hsf encoding autotransporter adhesin, and fhaB encoding filamentous haemagglutinin. Furthermore, based on genome-wide positive selection scanning, a total of 35 genes were subject to strong selection pressure. Extensive analyses of protein subcellular location indicated that membrane-associated genes were highly abundant among all positively selected genes. The detected amino acid sites undergoing adaptive selection were preferably located in extracellular space, perhaps associated with bacterial evasion of host immune responses. Our findings shed more light on conservation and distribution of virulence-associated genes across P. multocida strains. Meanwhile, this study provides a genetic context for future researches on the mechanism of adaptive evolution in P. multocida.

Can't you hear me knocking: contact-dependent competition and cooperation in bacteria

Microorganisms are in constant competition for growth niches and environmental resources. In Gram-negative bacteria, contact-dependent growth inhibition (CDI) systems link the fate of one cell with its immediate neighbor through touch-dependent, receptor-mediated toxin delivery. Though discovered for their ability to confer a competitive growth advantage, CDI systems also play significant roles in inter-sibling cooperation, promoting both auto-aggregation and biofilm formation. In this review, we detail the mechanisms of CDI toxin delivery and consider how toxin exchange between isogenic sibling cells could regulate gene expression.

CdiA Effectors Use Modular Receptor-Binding Domains To Recognize Target Bacteria

Contact-dependent growth inhibition (CDI) systems encode CdiA effectors, which bind to specific receptors on neighboring bacteria and deliver C-terminal toxin domains to suppress target cell growth. Two classes of CdiA effectors that bind distinct cell surface receptors have been identified, but the molecular basis of receptor specificity is not understood. Alignment of BamA-specific CdiA^EC93 from Escherichia coli EC93 and OmpC-specific CdiA^EC536 from E. coli 536 suggests that the receptor-binding domain resides within a central region that varies between the two effectors. In support of this hypothesis, we find that CdiA^EC93 fragments containing residues Arg1358 to Phe1646 bind specifically to purified BamA. Moreover, chimeric CdiA^EC93 that carries the corresponding sequence from CdiA^EC536 is endowed with OmpC-binding activity, demonstrating that this region dictates receptor specificity. A survey of E. coli CdiA proteins reveals two additional effector classes, which presumably recognize distinct receptors. Using a genetic approach, we identify the outer membrane nucleoside transporter Tsx as the receptor for a third class of CdiA effectors. Thus, CDI systems exploit multiple outer membrane proteins to identify and engage target cells. These results underscore the modularity of CdiA proteins and suggest that novel effectors can be constructed through genetic recombination to interchange different receptor-binding domains and toxic payloads.IMPORTANCE CdiB/CdiA two-partner secretion proteins mediate interbacterial competition through the delivery of polymorphic toxin domains. This process, known as contact-dependent growth inhibition (CDI), requires stable interactions between the CdiA effector protein and specific receptors on the surface of target bacteria. Here, we localize the receptor-binding domain to the central region of E. coli CdiA. Receptor-binding domains vary between CdiA proteins, and E. coli strains collectively encode at least four distinct effector classes. Further, we show that receptor specificity can be altered by exchanging receptor-binding regions, demonstrating the modularity of this domain. We propose that novel CdiA effectors are naturally generated through genetic recombination to interchange different receptor-binding domains and toxin payloads.

Genome-Wide Survey of Genes Under Positive Selection in Avian Pathogenic Escherichia coli Strains

The ability to obtain bacterial genomes from the same host has allowed for comparative studies that help in the understanding of the molecular evolution of specific pathotypes. Avian pathogenic Escherichia coli (APEC) is a group of extraintestinal strains responsible for causing colibacillosis in birds. APEC is also suggested to possess a role as a zoonotic agent. Despite its importance, APEC pathogenesis still has several cryptic pathogenic processes that need to be better understood. In this work, a genome-wide survey of eight APEC strains for genes with evidence of recombination revealed that ∼14% of the homologous groups evaluated present signs of recombination. Enrichment analyses revealed that nine Gene Ontology (GO) terms were significantly more represented in recombinant genes. Among these GO terms, several were noted to be ATP-related categories. The search for positive selection in these APEC genomes revealed 32 groups of homologous genes with evidence of positive selection. Among these groups, we found several related to cell metabolism, as well as several uncharacterized genes, beyond the well-known virulence factors ompC, lamB, waaW, waaL, and fliC. A GO term enrichment test showed a prevalence of terms related to bacterial cell contact with the external environment (e.g., viral entry into host cell, detection of virus, pore complex, bacterial-type flagellum filament C, and porin activity). Finally, the genes with evidence of positive selection were retrieved from genomes of non-APEC strains and tested as were done for APEC strains. The result revealed that none of the groups of genes presented evidence of positive selection, confirming that the analysis was effective in inferring positive selection for APEC and not for E. coli in general, which means that the study of the genes with evidence of positive selection identified in this study can contribute for the better understanding of APEC pathogenesis processes.

Evolutionary Genomics of an Ancient Prophage of the Order Sphingomonadales

The order Sphingomonadales, containing the families Erythrobacteraceae and Sphingomonadaceae, is a relatively less well-studied phylogenetic branch within the class Alphaproteobacteria. Prophage elements are present in most bacterial genomes and are important determinants of adaptive evolution. An “intact” prophage was predicted within the genome of Sphingomonas hengshuiensis strain WHSC-8 and was designated Prophage I^WHSC-8. Loci homologous to the region containing the first 22 open reading frames (ORFs) of Prophage I^WHSC-8 were discovered among the genomes of numerous Sphingomonadales. In 17 genomes, the homologous loci were co-located with an ORF encoding a putative superoxide dismutase. Several other lines of molecular evidence implied that these homologous loci represent an ancient temperate bacteriophage integration, and this horizontal transfer event pre-dated niche-based speciation within the order Sphingomonadales. The “stabilization” of prophages in the genomes of their hosts is an indicator of “fitness” conferred by these elements and natural selection. Among the various ORFs predicted within the conserved prophages, an ORF encoding a putative proline-rich outer membrane protein A was consistently present among the genomes of many Sphingomonadales. Furthermore, the conserved prophages in six Sphingomonas sp. contained an ORF encoding a putative spermidine synthase. It is possible that one or more of these ORFs bestow selective fitness, and thus the prophages continue to be vertically transferred within the host strains. Although conserved prophages have been identified previously among closely related genera and species, this is the first systematic and detailed description of orthologous prophages at the level of an order that contains two diverse families and many pigmented species.

Population Phylogenomics of Extraintestinal Pathogenic Escherichia coli

The emergence of genomics over the last 10 years has provided new insights into the evolution and virulence of extraintestinal Escherichia coli. By combining population genetics and phylogenetic approaches to analyze whole-genome sequences, it became possible to link genomic features to specific phenotypes, such as the ability to cause urinary tract infections. An E. coli chromosome can vary extensively in length, ranging from 4.3 to 6.2 Mb, encoding 4,084 to 6,453 proteins. This huge diversity is structured as a set of less than 2,000 genes (core genome) that are conserved between all the strains and a set of variable genes. Based on the core genome, the history of the species can be reliably reconstructed, revealing the recent emergence of phylogenetic groups A and B1 and the more ancient groups B2, F, and D. Urovirulence is most often observed in B2/F/D group strains and is a multigenic process involving numerous combinations of genes and specific alleles with epistatic interactions, all leading down multiple evolutionary paths. The genes involved mainly code for adhesins, toxins, iron capture systems, and protectins, as well as metabolic pathways and mutation-rate-control systems. However, the barrier between commensal and uropathogenic E. coli strains is difficult to draw as the factors that are responsible for virulence have probably also been selected to allow survival of E. coli as a commensal in the intestinal tract. Genomic studies have also demonstrated that infections are not the result of a unique and stable isolate, but rather often involve several isolates with variable levels of diversity that dynamically changes over time.

Evolution of bopA Gene in Burkholderia: A Case of Convergent Evolution as a Mechanism for Bacterial Autophagy Evasion

Autophagy is an important defense mechanism targeting intracellular bacteria to restrict their survival and growth. On the other hand, several intracellular pathogens have developed an antiautophagy mechanism to facilitate their own replication or intracellular survival. Up to now, no information about the origin or evolution of the antiautophagic genes in bacteria is available. BopA is an effector protein secreted by Burkholderia pseudomallei via the type three secretion system, and it has been shown to play a pivotal role in their escape from autophagy.  The evolutionary origin of bopA was examined in this work. Sequence similarity searches for BopA showed that no homolog of BopA was detected in eukaryotes. However, eukaryotic linear motifs were detected in BopA. The phylogenetic tree of the BopA proteins in our analysis is congruent with the species phylogeny derived from housekeeping genes. Moreover, there was no obvious difference in GC content values of bopA gene and their respective genomes. Integrated information on the taxonomic distribution, phylogenetic relationships, and GC content of the bopA gene of Burkholderia revealed that this gene was acquired via convergent evolution, not from eukaryotic host through horizontal gene transfer (HGT) event. This work has, for the first time, characterized the evolutionary mechanism of bacterial evasion of autophagy. The results of this study clearly demonstrated the role of convergent evolution in the evolution of how bacteria evade autophagy.

Adaptation of Escherichia coli traversing from the faecal environment to the urinary tract

The majority of extraintestinal pathogenic Escherichia coli (ExPEC) causing urinary tract infections (UTI) are found in the patient's own gut flora, but only limited knowledge is available on the potential adaptation that may occur in the bacteria in order to traverse the perineum and successfully infect the urinary tract. Here, matching pairs of faecal and UTI isolates from 42 patients were compared pairwise using in-depth whole-genome sequencing to investigate whether genetic changes were evident for successful colonization in these two different environments. The identified non-synonymous mutations (0-12 substitutions in each pair) were primarily associated to genes encoding virulence factors and nutrient metabolism; and indications of parallel evolution were observed in genes encoding the major phase-variable protein antigen 43, a toxin/antitoxin locus and haemolysin B. No differences in virulence potential were observed in a mouse UTI model for five matching faecal and UTI isolates with or without mutations in antigen 43 and haemolysin B. Variations in plasmid content were observed in only four of the 42 pairs. Although, we observed mutations in known UTI virulence genes for a few pairs, the majority showed no detectable differences with respect to mutations or mobilome when compared to their faecal counterpart. The results show that UPECs are successful in colonizing both the bladder and gut without adaptation.

CdiA Effectors from Uropathogenic Escherichia coli Use Heterotrimeric Osmoporins as Receptors to Recognize Target Bacteria

Many Gram-negative bacterial pathogens express contact-dependent growth inhibition (CDI) systems that promote cell-cell interaction. CDI⁺ bacteria express surface CdiA effector proteins, which transfer their C-terminal toxin domains into susceptible target cells upon binding to specific receptors. CDI⁺ cells also produce immunity proteins that neutralize the toxin domains delivered from neighboring siblings. Here, we show that CdiA^EC536 from uropathogenic Escherichia coli 536 (EC536) uses OmpC and OmpF as receptors to recognize target bacteria. E. coli mutants lacking either ompF or ompC are resistant to CDI^EC536-mediated growth inhibition, and both porins are required for target-cell adhesion to inhibitors that express CdiA^EC536. Experiments with single-chain OmpF fusions indicate that the CdiA^EC536 receptor is heterotrimeric OmpC-OmpF. Because the OmpC and OmpF porins are under selective pressure from bacteriophages and host immune systems, their surface-exposed loops vary between E. coli isolates. OmpC polymorphism has a significant impact on CDI^EC536 mediated competition, with many E. coli isolates expressing alleles that are not recognized by CdiA^EC536. Analyses of recombinant OmpC chimeras suggest that extracellular loops L4 and L5 are important recognition epitopes for CdiA^EC536. Loops L4 and L5 also account for much of the sequence variability between E. coli OmpC proteins, raising the possibility that CDI contributes to the selective pressure driving OmpC diversification. We find that the most efficient CdiA^EC536 receptors are encoded by isolates that carry the same cdi gene cluster as E. coli 536. Thus, it appears that CdiA effectors often bind preferentially to "self" receptors, thereby promoting interactions between sibling cells. As a consequence, these effector proteins cannot recognize nor suppress the growth of many potential competitors. These findings suggest that self-recognition and kin selection are important functions of CDI.

Bacterial pathogens often live in crowded communities where cells reside in close contact with one another. Many of these bacteria possess contact-dependent growth inhibition (CDI) systems, which allow cells to touch and inhibit each other using toxic CdiA proteins. CDI⁺ bacteria also produce immunity proteins that specifically protect the cell from the CdiA toxins of neighboring sibling cells. The CDI system from Escherichia coli EC93 was the first to be characterized and its CdiA toxin recognizes a receptor (BamA) that is identical in virtually all E. coli isolates. Here, we describe a different CDI system from uropathogenic E. coli 536, which causes urinary tract infections. In contrast to E. coli EC93, CdiA from E. coli 536 binds to receptor proteins (OmpC/OmpF) that vary widely between different E. coli isolates. Thus, uropathogenic E. coli preferentially bind and deliver toxins into sibling cells and other closely related E. coli strains. These results suggest that CDI systems distinguish between "self" and "non-self" cells. Moreover, because sibling cells are immune to CdiA-mediated growth inhibition, these findings raise the possibility that toxin exchange may be used for communication and cooperative behavior between genetically identical bacteria.

Comparative Genomic Analysis among Four Representative Isolates of Phytophthora sojae Reveals Genes under Evolutionary Selection

Comparative genomic analysis is useful for identifying genes affected by evolutionary selection and for studying adaptive variation in gene functions. In Phytophthora sojae, a model oomycete plant pathogen, the related study is lacking. We compared sequence data among four isolates of P. sojae, which represent its four major genotypes. These isolates exhibited >99.688%, >99.864%, and >98.981% sequence identities at genome, gene, and non-gene regions, respectively. One hundred and fifty-three positive selection and 139 negative selection candidate genes were identified. Between the two categories of genes, the positive selection genes were flanked by larger intergenic regions, poorly annotated in function, and less conserved; they had relatively lower transcription levels but many genes had increased transcripts during infection. Genes coding for predicted secreted proteins, particularly effectors, were overrepresented in positive selection. Several RxLR effector genes were identified as positive selection genes, exhibiting much stronger positive selection levels. In addition, candidate genes with presence/absence polymorphism were analyzed. This study provides a landscape of genomic variation among four representative P. sojae isolates and characterized several evolutionary selection-affected gene candidates. The results suggest a relatively covert two-speed genome evolution pattern in P. sojae and will provide clues for identification of new virulence factors in the oomycete plant pathogens.

Molecular Evolution of the Yersinia Major Outer Membrane Protein C (OmpC)

The genus Yersinia includes species with a wide range of eukaryotic hosts (from fish, insects, and plants to mammals and humans). One of the major outer membrane proteins, the porin OmpC, is preferentially expressed in the host gut, where osmotic pressure, temperature, and the concentrations of nutrients and toxic products are relatively high. We consider here the molecular evolution and phylogeny of Yersinia ompC. The maximum likelihood gene tree reflects the macroevolution processes occurring within the genus Yersinia. Positive selection and horizontal gene transfer are the key factors of ompC diversification, and intraspecies recombination was revealed in two Yersinia species. The impact of recombination on ompC evolution was different from that of another major porin gene, ompF, possibly due to the emergence of additional functions and conservation of the basic transport function. The predicted antigenic determinants of OmpC were located in rapidly evolving regions, which may indicate the evolutionary mechanisms of Yersinia adaptation to the host immune system.

Gene flow, recombination, and positive selection in Stenotrophomonas maltophilia: mechanisms underlying the diversity of the widespread opportunistic pathogen

Stenotrophomonas maltophilia is a global multidrug-resistant human opportunistic pathogen in clinical environments. Stenotrophomonas maltophilia is also ubiquitous in aqueous environments, soil, and plants. Various molecular typing methods have revealed that S. maltophilia exhibits high levels of phenotypic and genotypic diversity. However, information regarding the genomic diversity within S. maltophilia and the corresponding genetic mechanisms resulting in said diversity remain scarce. The genome sequences of 17 S. maltophilia strains were selected to investigate the mechanisms contributing to genetic diversity at the genome level. The core and large pan-genomes of the species were first estimated, resulting in a large, open pan-genome. A species phylogeny was also reconstructed based on 344 orthologous genes with one copy per genome, and the contribution of four evolutionary mechanisms to the species genome diversity was quantified: 15%-35% of the genes showed evidence for recombination, 0%-25% of the genes in one genome were likely gained, 0%-44% of the genes in some genomes were likely lost, and less than 0.3% of the genes in a genome were under positive selection pressures. We observed that, among the four main mechanisms, homologous recombination plays a key role in maintaining diversity in S. maltophilia. In this study, we provide an overview of evolution in S. maltophilia to provide a better understanding of its evolutionary dynamics and its relationship with genome diversity.

Kin Recognition in Bacteria

The ability of bacteria to recognize kin provides a means to form social groups. In turn these groups can lead to cooperative behaviors that surpass the ability of the individual. Kin recognition involves specific biochemical interactions between a receptor(s) and an identification molecule(s). Recognition specificity, ensuring that nonkin are excluded and kin are included, is critical and depends on the number of loci and polymorphisms involved. After recognition and biochemical perception, the common ensuing cooperative behaviors include biofilm formation, quorum responses, development, and swarming motility. Although kin recognition is a fundamental mechanism through which cells might interact, microbiologists are only beginning to explore the topic. This review considers both molecular and theoretical aspects of bacterial kin recognition. Consideration is also given to bacterial diversity, genetic relatedness, kin selection theory, and mechanisms of recognition.