Evolution in an oncogenic bacterial species with extreme genome plasticity: Helicobacter pylori East Asian genomes

<p>Abstract</p> <p>Background</p> <p>The genome of <it>Helicobacter pylori</it>, an oncogenic bacterium in the human stomach, rapidly evolves and shows wide geographical divergence. The high incidence of stomach cancer in East Asia might be related to bacterial genotype. We used newly developed comparative methods to follow the evolution of East Asian <it>H. pylori </it>genomes using 20 complete genome sequences from Japanese, Korean, Amerind, European, and West African strains.</p> <p>Results</p> <p>A phylogenetic tree of concatenated well-defined core genes supported divergence of the East Asian lineage (hspEAsia; Japanese and Korean) from the European lineage ancestor, and then from the Amerind lineage ancestor. Phylogenetic profiling revealed a large difference in the repertoire of outer membrane proteins (including <it>oipA</it>, <it>hopMN</it>, <it>babABC</it>, <it>sabAB </it>and <it>vacA-2</it>) through gene loss, gain, and mutation. All known functions associated with molybdenum, a rare element essential to nearly all organisms that catalyzes two-electron-transfer oxidation-reduction reactions, appeared to be inactivated. Two pathways linking acetyl~CoA and acetate appeared intact in some Japanese strains. Phylogenetic analysis revealed greater divergence between the East Asian (hspEAsia) and the European (hpEurope) genomes in proteins in host interaction, specifically virulence factors (<it>tipα</it>), outer membrane proteins, and lipopolysaccharide synthesis (human Lewis antigen mimicry) enzymes. Divergence was also seen in proteins in electron transfer and translation fidelity (<it>miaA, tilS</it>), a DNA recombinase/exonuclease that recognizes genome identity (<it>addA</it>), and DNA/RNA hybrid nucleases (<it>rnhAB</it>). Positively selected amino acid changes between hspEAsia and hpEurope were mapped to products of <it>cagA</it>, <it>vacA</it>, <it>homC </it>(outer membrane protein), <it>sotB </it>(sugar transport), and a translation fidelity factor (<it>miaA</it>). Large divergence was seen in genes related to antibiotics: <it>frxA </it>(metronidazole resistance), <it>def </it>(peptide deformylase, drug target), and <it>ftsA </it>(actin-like, drug target).</p> <p>Conclusions</p> <p>These results demonstrate dramatic genome evolution within a species, especially in likely host interaction genes. The East Asian strains appear to differ greatly from the European strains in electron transfer and redox reactions. These findings also suggest a model of adaptive evolution through proteome diversification and selection through modulation of translational fidelity. The results define <it>H. pylori </it>East Asian lineages and provide essential information for understanding their pathogenesis and designing drugs and therapies that target them.</p

Antibiotic treatment regimes as a driver of the global population dynamics of a major gonorrhea lineage

The Neisseria gonorrhoeae multilocus sequence type (ST) 1901 is among the lineages most commonly associated with treatment failure. Here, we analyze a global collection of ST-1901 genomes to shed light on the emergence and spread of alleles associated with reduced susceptibility to extended-spectrum cephalosporins (ESCs). The genetic diversity of ST-1901 falls into a minor and a major clade, both of which were inferred to have originated in East Asia. The dispersal of the major clade from Asia happened in two separate waves expanding from ∼1987 and 1996, respectively. Both waves first reached North America, and from there spread to Europe and Oceania, with multiple secondary reintroductions to Asia. The ancestor of the second wave acquired the penA 34.001 allele, which significantly reduces susceptibility to ESCs. Our results suggest that the acquisition of this allele granted the second wave a fitness advantage at a time when ESCs became the key drug class used to treat gonorrhea. Following its establishment globally, the lineage has served as a reservoir for the repeated emergence of clones fully resistant to the ESC ceftriaxone, an essential drug for effective treatment of gonorrhea. We infer that the effective population sizes of both clades went into decline as treatment schemes shifted from fluoroquinolones via ESC monotherapy to dual therapy with ceftriaxone and azithromycin in Europe and the United States. Despite the inferred recent population size decline, the short evolutionary path from the penA 34.001 allele to alleles providing full ceftriaxone resistance is a cause of concern

Genomic surveillance of Neisseria gonorrhoeae to investigate the distribution and evolution of antimicrobial-resistance determinants and lineages

The first extensively drug resistant (XDR) Neisseria gonorrhoeae strain with high resistance to the extended-spectrum cephalosporin ceftriaxone was identified in 2009 in Japan, but no other strain with this antimicrobial-resistance profile has been reported since. However, surveillance to date has been based on phenotypic methods and sequence typing, not genome sequencing. Therefore, little is known about the local population structure at the genomic level, and how resistance determinants and lineages are distributed and evolve. We analysed the whole-genome sequence data and the antimicrobial-susceptibility testing results of 204 strains sampled in a region where the first XDR ceftriaxone-resistant N. gonorrhoeae was isolated, complemented with 67 additional genomes from other time frames and locations within Japan. Strains resistant to ceftriaxone were not found, but we discovered a sequence type (ST)7363 sub-lineage susceptible to ceftriaxone and cefixime in which the mosaic penA allele responsible for reduced susceptibility had reverted to a susceptible allele by recombination. Approximately 85 % of isolates showed resistance to fluoroquinolones (ciprofloxacin) explained by linked amino acid substitutions at positions 91 and 95 of GyrA with 99 % sensitivity and 100 % specificity. Approximately 10 % showed resistance to macrolides (azithromycin), for which genetic determinants are less clear. Furthermore, we revealed different evolutionary paths of the two major lineages: single acquisition of penA X in the ST7363-associated lineage, followed by multiple independent acquisitions of the penA X and XXXIV in the ST1901-associated lineage. Our study provides a detailed picture of the distribution of resistance determinants and disentangles the evolution of the two major lineages spreading worldwide

Rapid evolution of distinct Helicobacter pylori subpopulations in the Americas

For the last 500 years, the Americas have been a melting pot both for genetically diverse humans and for the pathogenic and commensal organisms associated with them. One such organism is the stomach-dwelling bacterium Helicobacter pylori, which is highly prevalent in Latin America where it is a major current public health challenge because of its strong association with gastric cancer. By analyzing the genome sequence of H. pylori isolated in North, Central and South America, we found evidence for admixture between H. pylori of European and African origin throughout the Americas, without substantial input from pre-Columbian (hspAmerind) bacteria. In the US, strains of African and European origin have remained genetically distinct, while in Colombia and Nicaragua, bottlenecks and rampant genetic exchange amongst isolates have led to the formation of national gene pools. We found three outer membrane proteins with atypical levels of Asian ancestry in American strains, as well as alleles that were nearly fixed specifically in South American isolates, suggesting a role for the ethnic makeup of hosts in the colonization of incoming strains. Our results show that new H. pylori subpopulations can rapidly arise, spread and adapt during times of demographic flux, and suggest that differences in transmission ecology between high and low prevalence areas may substantially affect the composition of bacterial populations

Recombination-Mediated Host Adaptation by Avian Staphylococcus aureus

Staphylococcus aureus are globally disseminated among farmed chickens causing skeletal muscle infections, dermatitis, and septicaemia. The emergence of poultry-associated lineages has involved zoonotic transmission from humans to chickens but questions remain about the specific adaptations that promote proliferation of chicken pathogens. We characterized genetic variation in a population of genome-sequenced S. aureus isolates of poultry and humanorigin. Genealogical analysis identified a dominant poultry-associated sequence cluster within the CC5 clonal complex. Poultry and human CC5 isolates were significantly distinct from each other and more recombination events were detected in the poultry isolates. We identified 44 recombination events in 33 genes along the branch extending to the poultry-specific CC5 cluster, and 47 genes were found more often in CC5 poultry isolates compared with those from humans. Many of these gene sequences were common in chicken isolates from other clonal complexes suggesting horizontal gene transfer among poultry associated lineages. Consistent with functional predictions for putative poultry-associated genes, poultry isolates showed enhanced growth at 42 degrees C and greater erythrocyte lysis on chicken blood agar in comparison with human isolates. By combining phenotype information with evolutionary analyses of staphylococcal genomes, we provide evidence of adaptation, following a human-to-poultry host transition. This has important implications for the emergence and dissemination of new pathogenic clones associated with modern agriculture.Peer reviewe

Repeated out-of-Africa expansions of Helicobacter pylori driven by replacement of deleterious mutations

Erratum in: Nat Commun. 2023 Mar 20;14(1):1539. doi: 10.1038/s41467-023-37302-5.Helicobacter pylori lives in the human stomach and has a population structure resembling that of its host. However, H. pylori fromEurope and the Middle East trace substantially more ancestry from modern African populations than the humans that carry them. Here, we use a collection of Afro-Eurasian H. pylori genomes to show that this African ancestry is due to at least three distinct admixture events. H. pylori from East Asia, which have undergone little admixture, have accumulated many more non-synonymous mutations than African strains. European and Middle Eastern bacteria have elevated African ancestry at the sites of these mutations, implying selection to remove them during admixture. Simulations show that population fitness can be restored after bottlenecks bymigration and subsequent admixture of small numbers of bacteria from non-bottlenecked populations. We conclude that recent spread of African DNA has been driven by deleterious mutations accumulated during the original out-of-Africa bottleneck.This work was supported by Sequencing Grants-in-aid for Scientific Research from the Ministry of Education, Culture, Sports, Science, and Technology (MEXT) of Japan (221S0002, 18KK0266, 19H03473, 21H00346 and 22H02871) to Y.Y. F.F.V. is financed by FCT through Assistant Researcher grant CEECIND/03023/2017 and a project grant PTDC/BTM-TEC/3238/ 2020. I.K. studentship was funded by the National Strategic Reference Framework Operational Program “Competitiveness, Entrepreneurship and Innovation” (NSRF 2014-2020, project No. MIS5002486) and sequencing of strains was supported by the InfeNeutra Project (NSRF 2007-2013, project no. MIS450598) of the Ministry of Culture and Edu- cation, Greece. K.T. and the sequencing of KI isolates was supported by Erik Philip-Sörensen Foundation grant G2016-08, and Swedish Society for Medical research (SSMF). All primary bioinformatics and parts of the comparative genomics were performed on resources provided by Swedish National Infrastructure for Computing (SNIC) through Uppsala Multidisciplinary Center for Advanced Computational Science (UPPMAX) under projects snic2018-8-24 and uppstore2017270. Work by S.S. was supported by the German Research Foundation (DFG, project number 158 989 968–SFB 900/A1) and by the Bavarian Ministry of Sci- ence and the Arts in the framework of the Bavarian Research Network “New Strategies Against Multi-Resistant Pathogens by Means of Digital Networking—bayresq.net”. D.F. was supported by Shanghai Municipal Science and Technology Major Project No. 2019SHZDZX02.info:eu-repo/semantics/publishedVersio

Epistasis, core-genome disharmony, and adaptation in recombining bacteria

Recombination of short DNA fragments via horizontal gene transfer (HGT) can introduce beneficial alleles, create genomic disharmony through negative epistasis, and create adaptive gene combinations through positive epistasis. For non-core (accessory) genes, the negative epistatic cost is likely to be minimal because the incoming genes have not co-evolved with the recipient genome and are frequently observed as tightly linked cassettes with major effects. By contrast, interspecific recombination in the core genome is expected to be rare because disruptive allelic replacement is likely to introduce negative epistasis. Why then is homologous recombination common in the core of bacterial genomes? To understand this enigma, we take advantage of an exceptional model system, the common enteric pathogens Campylobacter jejuni and C. coli that are known for very high magnitude interspecies gene flow in the core genome. As expected, HGT does indeed disrupt co-adapted allele pairings, indirect evidence of negative epistasis. However, multiple HGT events enable recovery of the genome’s co-adaption between introgressing alleles, even in core metabolism genes (e.g., formate dehydrogenase). These findings demonstrate that, even for complex traits, genetic coalitions can be decoupled, transferred, and independently reinstated in a new genetic background—facilitating transition between fitness peaks. In this example, the two-step recombinational process is associated with C. coli that are adapted to the agricultural niche

Disease-associated genotypes of the commensal skin bacterium Staphylococcus epidermidis

Some of the most common infectious diseases are caused by bacteria that naturally colonise humans asymptomatically. Combating these opportunistic pathogens requires an understanding of the traits that differentiate infecting strains from harmless relatives. Staphylococcus epidermidis is carried asymptomatically on the skin and mucous membranes of virtually all humans but is a major cause of nosocomial infection associated with invasive procedures. Here we address the underlying evolutionary mechanisms of opportunistic pathogenicity by combining pangenome-wide association studies and laboratory microbiology to compare S. epidermidis from bloodstream and wound infections and asymptomatic carriage. We identify 61 genes containing infection-associated genetic elements (k-mers) that correlate with in vitro variation in known pathogenicity traits (biofilm formation, cell toxicity, interleukin-8 production, methicillin resistance). Horizontal gene transfer spreads these elements, allowing divergent clones to cause infection. Finally, Random Forest model prediction of disease status (carriage vs. infection) identifies pathogenicity elements in 415 S. epidermidis isolates with 80% accuracy, demonstrating the potential for identifying risk genotypes pre-operatively.Peer reviewe

Genome-wide association of functional traits linked with<i> Campylobacter jejuni </i>survival from farm to fork

Campylobacter jejuni is a major cause of bacterial gastroenteritis worldwide, primarily associated with the consumption of contaminated poultry. C. jejuni lineages vary in host range and prevalence in human infection, suggesting differences in survival throughout the poultry processing chain. From 7,343 MLST-characterised isolates, we sequenced 600 C. jejuni and C. coli isolates from various stages of poultry processing and clinical cases. A genome-wide association study (GWAS) in C. jejuni ST-21 and ST-45 complexes identified genetic elements over-represented in clinical isolates that increased in frequency throughout the poultry processing chain. Disease-associated SNPs were distinct in these complexes, sometimes organised in haplotype blocks. The function of genes containing associated elements was investigated, demonstrating roles for cj1377c in formate metabolism, nuoK in aerobic survival and oxidative respiration, and cj1368-70 in nucleotide salvage. This work demonstrates the utility of GWAS for investigating transmission in natural zoonotic pathogen populations and provides evidence that major C. jejuni lineages have distinct genotypes associated with survival, within the host specific niche, from farm to fork. </p