Rapid Gene Turnover as a Significant Source of Genetic Variation in a Recently Seeded Population of a Healthcare-Associated Pathogen

Evolutionary dynamics of the accessory genomes of Staphylococcus aureus

Staphylococcus aureus is a ubiquitous commensal and opportunistic bacterial pathogen that can cause a wide gamut of infections, which are exacerbated by the presence of multidrug-resistant and methicillin-resistant S. aureus. S. aureus is genetically heterogeneous and consists of numerous distinct lineages. Using 558 complete genomes of S. aureus, we aim to determine how the accessory genome content among phylogenetic lineages of S. aureus is structured and has evolved. Bayesian hierarchical clustering identified 10 sequence clusters, of which seven contained major sequence types (ST 1, 5, 8, 30, 59, 239, and 398). The seven sequence clusters differed in their accessory gene content, including genes associated with antimicrobial resistance and virulence. Focusing on the two largest clusters, BAPS8 and BAPS10, and each consisting mostly of ST5 and ST8, respectively, we found that the structure and connected components in the co-occurrence networks of accessory genomes varied between them. These differences are explained, in part, by the variation in the rates at which the two sequence clusters gained and lost accessory genes, with the highest rate of gene accumulation occurring recently in their evolutionary histories. We also identified a divergent group within BAPS10 that has experienced high gene gain and loss early in its history. Together, our results show highly variable and dynamic accessory genomes in S. aureus that are structured by the history of the specific lineages that carry them.IMPORTANCEStaphylococcus aureus is an opportunistic, multi-host pathogen that can cause a variety of benign and life-threatening infections. Our results revealed considerable differences in the structure and evolution of the accessory genomes of major lineages within S. aureus. Such genomic variation within a species can have important implications on disease epidemiology, pathogenesis of infection, and interactions with the vertebrate host. Our findings provide important insights into the underlying genetic basis for the success of S. aureus as a highly adaptable and resistant pathogen, which will inform current efforts to control and treat staphylococcal diseases.

Pandrug-resistant Acinetobacter baumannii from different clones and regions in Mexico have a similar plasmid carrying the blaOXA-72 gene

Background

Multidrug-resistant Acinetobacter baumannii is a common hospital-acquired pathogen. The increase in antibiotic resistance is commonly due to the acquisition of mobile genetic elements carrying antibiotic resistance genes. To comprehend this, we analyzed the resistome and virulome of Mexican A. baumannii multidrug-resistant isolates.

Methods

Six clinical strains of A. baumannii from three Mexican hospitals were sequenced using the Illumina platform, the genomes were assembled with SPAdes and annotated with Prokka. Plasmid SPAdes and MobRecon were used to identify the potential plasmid sequences. Sequence Type (ST) assignation under the MLST Oxford scheme was performed using the PubMLST database. Homologous gene search for known virulent factors was performed using the virulence factor database VFDB and an in silico prediction of the resistome was conducted via the ResFinder databases.

Results

The six strains studied belong to different STs and clonal complexes (CC): two strains were ST208 and one was ST369; these two STs belong to the same lineage CC92, which is part of the international clone (IC) 2. Another two strains were ST758 and one was ST1054, both STs belonging to the same lineage CC636, which is within IC5. The resistome analysis of the six strains identified between 7 to 14 antibiotic resistance genes to different families of drugs, including beta-lactams, aminoglycosides, fluoroquinolones and carbapenems. We detected between 1 to 4 plasmids per strain with sizes from 1,800 bp to 111,044 bp. Two strains from hospitals in Mexico City and Guadalajara had a plasmid each of 10,012 bp pAba78r and pAba79f, respectively, which contained the bla OXA-72 gene. The structure of this plasmid showed the same 13 genes in both strains, but 4 of them were inverted in one of the strains. Finally, the six strains contain 49 identical virulence genes related to immune response evasion, quorum-sensing, and secretion systems, among others.

Conclusion

Resistance to carbapenems due to pAba78r and pAba79f plasmids in Aba pandrug-resistant strains from different geographic areas of Mexico and different clones was detected. Our results provide further evidence that plasmids are highly relevant for the horizontal transfer of antibiotic resistance genes between different clones of A. baumannii.

The whole-genome molecular epidemiology of sequential isolates of Acinetobacter baumannii colonizing the rectum of patients in an adult intensive care unit of a tertiary hospital

IMPORTANCE

Acinetobacter baumannii is a multidrug-resistant nosocomial pathogen that colonizes and infects debilitated patients in the ICU. There is very little information on the genomic characteristics of colonizing strains. This information is important to understand the evolution of lineages of A. baumannii that develop resistance while patients receive antibiotic treatment in the ICU. Our study demonstrated different patterns of colonization of the rectum of ICU patients with different STs of A. baumannii while one ST colonized all patients. Some STs carried more antibiotic resistance genes compared to others. However, there was a correlation between ST and a particular resistance gene profile. Our results further elucidate the dynamics of enteric colonization of this opportunistic pathogen.

Genomic and phenotypic features of Acinetobacter baumannii isolated from oil reservoirs reveal a novel subspecies specialized in degrading hazardous hydrocarbons

The genus Acinetobacter encompasses biotechnologically relevant species and nosocomial pathogens. In this study, nine isolates recovered from different oil reservoir samples showed the ability to grow with petroleum as the only carbon source and possessed the ability to emulsify kerosene. The whole genomes of the nine strains were sequenced and analyzed. The average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values of all strains were compared to the reference strains, and the results were below the reference values (<97.88 and 82, respectively), suggesting that the isolates belong to a new subspecies of Acinetobacter baumannii. The name Acinetobacter baumannii oleum ficedula is proposed. A comparison of the whole genome repertoire of 290 Acinetobacter species indicated that the strains in this study resemble non-pathogenic Acinetobacter strains. However, the new isolates resemble A. baumannii when comparing virulence factors. The isolates in this study carry many genes involved in hydrocarbon degradation, indicating the potential to degrade most toxic compounds listed by environmental regulatory agencies such as ATSDR, EPA, and CONAMA. In addition, despite the absence of known biosurfactant or bioemulsifier genes, the strains showed emulsifying activity, suggesting the presence of new pathways or genes related to this process. This study investigated the genomic, phenotypic, and biochemical features of the novel environmental subspecies A. baumannii oleum ficedula, revealing their potential to degrade hydrocarbons and to produce biosurfactants or bioemulsifiers. Applying these environmental subspecies in bioaugmentation strategies sheds light on future approaches to bioremediation. The study shows the importance of genomic analysis of environmental strains and their inclusion in metabolic pathways databases, highlighting unique enzymes/alternative pathways for consuming hazardous hydrocarbons.

Interactions between host epithelial cells and Acinetobacter baumannii promote the emergence of highly antibiotic resistant and highly mucoid strains

Acinetobacter baumannii is an important nosocomial pathogen. Upon colonizing a host, A. baumannii are subjected to selective pressure by immune defenses as they adapt to the host environment. However, the mechanism of this pathoadaptation is unknown. Here, we established an in vitro system to evolve A. baumannii driven by the continuous selective pressure exerted by epithelial cells, and we used a combination of experimental evolution, phenotypic characterization and multi-omics analysis to address the underlying mechanism. When continuously exposed to selective pressure by pulmonary epithelial cells, A. baumannii showed ptk mutation-mediated mucoid conversion (reduced adhesion and increased anti-phagocytic ability) by enhancement of capsular exopolysaccharide chain length; rsmG mutation-mediated deficiency of 7-methylguanosine modification in the 524th nucleotide of 16S rRNA, which increased ribosome translation efficiency; and rnaseI mutation-mediated changes in outer membrane permeability and efflux pump expression. Together, these mutations altered susceptibility to a variety of antimicrobial agents, including the novel antibiotic cefiderocol, by regulating siderophore and siderophore-receptor biosynthesis. In conclusion, pulmonary epithelial cells modulate A. baumannii pathoadaptation, implicating the host–microbe interaction in the survival and persistence of A. baumannii.

Using Genomes and Evolutionary Analyses to Screen for Host-Specificity and Positive Selection in the Plant Pathogen Xylella fastidiosa

Xylella fastidiosa infects several economically important crops in the Americas, and it also recently emerged in Europe. Here, using a set of Xylella genomes reflective of the genus-wide diversity, we performed a pan-genome analysis based on both core and accessory genes for two purposes: (i) to test associations between genetic divergence and plant host species and (ii) to identify positively selected genes that are potentially involved in arms-race dynamics. For the former, tests yielded significant evidence for the specialization of X. fastidiosa to plant host species. This observation contributes to a growing literature suggesting that the phylogenetic history of X. fastidiosa lineages affects the host range. For the latter, our analyses uncovered evidence of positive selection across codons for 5.3% (67 of 1,257) of the core genes and 5.4% (201 of 3,691) of the accessory genes. These genes are candidates to encode interacting factors with plant and insect hosts. Most of these genes had unknown functions, but we did identify some tractable candidates, including nagZ_2, which encodes a beta-glucosidase that is important for Neisseria gonorrhoeae biofilm formation; cya, which modulates gene expression in pathogenic bacteria, and barA, a membrane associated histidine kinase that has roles in cell division, metabolism, and pili formation.

IMPORTANCEXylella fastidiosa causes devasting diseases to several critical crops. Because X. fastidiosa colonizes and infects many plant species, it is important to understand whether the genome of X. fastidiosa has genetic determinants that underlie specialization to specific host plants. We analyzed genome sequences of X. fastidiosa to investigate evolutionary relationships and to test for evidence of positive selection on specific genes. We found a significant signal between genome diversity and host plants, consistent with bacterial specialization to specific plant hosts. By screening for positive selection, we identified both core and accessory genes that may affect pathogenicity, including genes involved in biofilm formation.

Pangenome evolution in environmentally transmitted symbionts of deep-sea mussels is governed by vertical inheritance

Microbial pangenomes vary across species; their size and structure are determined by genetic diversity within the population and by gene loss and horizontal gene transfer (HGT). Many bacteria are associated with eukaryotic hosts where the host colonization dynamics may impact bacterial genome evolution. Host-associated lifestyle has been recognized as a barrier to HGT in parentally transmitted bacteria. However, pangenome evolution of environmentally acquired symbionts remains understudied, often due to limitations in symbiont cultivation. Using high-resolution metagenomics, here we study pangenome evolution of two co-occurring endosymbionts inhabiting Bathymodiolus brooksi mussels from a single cold seep. The symbionts, sulfur-oxidizing (SOX) and methane-oxidizing (MOX) gamma-proteobacteria, are environmentally acquired at an early developmental stage and individual mussels may harbor multiple strains of each symbiont species. We found differences in the accessory gene content of both symbionts across individual mussels, which are reflected by differences in symbiont strain composition. Compared to core genes, accessory genes are enriched in genome plasticity functions. We found no evidence for recent horizontal gene transfer between both symbionts. A comparison between the symbiont pangenomes revealed that the MOX population is less diverged and contains fewer accessory genes, supporting that the MOX association with B. brooksi is more recent in comparison to that of SOX. Our results show that the pangenomes of both symbionts evolved mainly by vertical inheritance. We conclude that genome evolution of environmentally transmitted symbionts that associate with individual hosts over their lifetime is affected by a narrow symbiosis where the frequency of HGT is constrained.​.

Genome diversity of domesticated Acinetobacter baumannii ATCC 19606T strains

Acinetobacter baumannii has emerged as an important opportunistic pathogen worldwide, being responsible for large outbreaks for nosocomial infections, primarily in intensive care units. A. baumannii ATCC 19606^T is the species type strain, and a reference organism in many laboratories due to its low virulence, amenability to genetic manipulation and extensive antibiotic susceptibility. We wondered if frequent propagation of A. baumannii ATCC 19606^T in different laboratories may have driven micro- and macro-evolutionary events that could determine inter-laboratory differences of genome-based data. By combining Illumina MiSeq, MinION and Sanger technologies, we generated a high-quality whole-genome sequence of A. baumannii ATCC 19606^T, then performed a comparative genome analysis between A. baumannii ATCC 19606^T strains from several research laboratories and a reference collection. Differences between publicly available ATCC 19606^T genome sequences were observed, including SNPs, macro- and micro-deletions, and the uneven presence of a 52 kb prophage belonging to genus Vieuvirus. Two plasmids, pMAC and p1ATCC19606, were invariably detected in all tested strains. The presence of a putative replicase, a replication origin containing four 22-mer direct repeats, and a toxin-antitoxin system implicated in plasmid stability were predicted by in silico analysis of p1ATCC19606, and experimentally confirmed. This work refines the sequence, structure and functional annotation of the A. baumannii ATCC 19606^T genome, and highlights some remarkable differences between domesticated strains, likely resulting from genetic drift.

The promiscuous and highly mobile resistome of Acinetobacter baumannii

Antimicrobial resistance (AR) is a major global threat to public health. Understanding the population dynamics of AR is critical to restrain and control this issue. However, no study has provided a global picture of the whole resistome of Acinetobacter baumannii, a very important nosocomial pathogen. Here we analyse 1450+ genomes (covering >40 countries and >4 decades) to infer the global population dynamics of the resistome of this species. We show that gene flow and horizontal transfer have driven the dissemination of AR genes in A. baumannii. We found considerable variation in AR gene content across lineages. Although the individual AR gene histories have been affected by recombination, the AR gene content has been shaped by the phylogeny. Furthermore, many AR genes have been transferred to other well-known pathogens, such as Pseudomonas aeruginosa or Klebsiella pneumoniae. Despite using this massive data set, we were not able to sample the whole diversity of AR genes, which suggests that this species has an open resistome. Our results highlight the high mobilization risk of AR genes between important pathogens. On a broader perspective, this study gives a framework for an emerging perspective (resistome-centric) on the genomic epidemiology (and surveillance) of bacterial pathogens.

Accessory Genomic Epidemiology of Cocirculating Acinetobacter baumannii Clones

Acinetobacter baumannii has become one of the most important multidrug-resistant nosocomial pathogens all over the world. Nonetheless, very little is known about the diversity of A. baumannii lineages coexisting in hospital settings. Here, using whole-genome sequencing, epidemiological data, and antimicrobial susceptibility tests, we uncover the transmission dynamics of extensive and multidrug-resistant A. baumannii in a tertiary hospital over a decade. Our core genome phylogeny of almost 300 genomes suggests that there were several introductions of lineages from international clone 2 into the hospital. The molecular dating analysis shows that these introductions happened in 2006, 2007, and 2013. Furthermore, using the accessory genome, we show that these lineages were extensively disseminated across many wards in the hospital. Our results demonstrate that accessory genome variation can be a very powerful tool for conducting genomic epidemiology. We anticipate future studies employing the accessory genome along with the core genome as a powerful phylogenomic strategy to track bacterial transmissions over very short microevolutionary scales. IMPORTANCE Whole-genome sequencing for epidemiological investigations (genomic epidemiology) has been of paramount importance to understand the transmission dynamics of many bacterial (and nonbacterial) pathogens. Commonly, variation in the core genome, single nucleotide polymorphisms (SNPs), is employed to carry out genomic epidemiology. However, at very short periods of time, the core genome might not have accumulated enough variation (sufficient SNPs) to tell apart isolates. In this scenario, gene content variation in the accessory genome can be an option to conduct genomic epidemiology. Here, we used the accessory genome, as well as the core genome, to uncover the transmission dynamics of extensive and multidrug-resistant A. baumannii in a tertiary hospital for a decade. Our study shows that accessory genome variation can be a very powerful tool for conducting genomic epidemiology.

A novel vieuvirus from multidrug-resistant Acinetobacter baumannii

Bacteriophages are considered the most abundant biological entities on earth, and they are able to modulate the populations of their bacterial hosts. Although the potential of bacteriophages has been accepted as an alternative strategy to combat multidrug-resistant pathogenic bacteria, there still exists a considerable knowledge gap regarding their genetic diversity, which hinders their use as antimicrobial agents. In this study, we undertook a genomic and phylogenetic characterization of the phage Ab11510-phi, which was isolated from a multidrug-resistant Acinetobacter baumannii strain (Ab11510). We found that Ab11510-phi has a narrow host range and belongs to a small group of transposable phages of the genus Vieuvirus that have only been reported to infect Acinetobacter bacteria. Finally, we showed that Ab11510-phi (as well as other vieuvirus phages) has a high level of mosaicism. On a broader level, we demonstrate that comparative genomics and phylogenetic analysis are necessary tools for the proper characterization of phage diversity.

Whole-Genome Sequencing for Investigating a Health Care-Associated Outbreak of Carbapenem-Resistant Acinetobacter baumannii

Carbapenem-resistant Acinetobacter baumannii (CRAB) outbreaks in hospital settings challenge the treatment of patients and infection control. Understanding the relatedness of clinical isolates is important in distinguishing outbreak isolates from sporadic cases. This study investigated 11 CRAB isolates from a hospital outbreak by whole-genome sequencing (WGS), utilizing various bioinformatics tools for outbreak analysis. The results of multilocus sequence typing (MLST), single nucleotide polymorphism (SNP) analysis, and phylogenetic tree analysis by WGS through web-based tools were compared, and repetitive element polymerase chain reaction (rep-PCR) typing was performed. Through the WGS of 11 A. baumannii isolates, three clonal lineages were identified from the outbreak. The coexistence of bla_OXA-23, bla_OXA-66, bla_ADC-25, and armA with additional aminoglycoside-inactivating enzymes, predicted to confer multidrug resistance, was identified in all isolates. The MLST Oxford scheme identified three types (ST191, ST369, and ST451), and, through whole-genome MLST and whole-genome SNP analyses, different clones were found to exist within the MLST types. wgSNP showed the highest discriminatory power with the lowest similarities among the isolates. Using the various bioinformatics tools for WGS, CRAB outbreak analysis was applicable and identified three discrete clusters differentiating the separate epidemiologic relationships among the isolates.

Time-calibrated genomic evolution of a monomorphic bacterium during its establishment as an endemic crop pathogen

Horizontal gene transfer is of major evolutionary importance as it allows for the redistribution of phenotypically important genes among lineages. Such genes with essential functions include those involved in resistance to antimicrobial compounds and virulence factors in pathogenic bacteria. Understanding gene turnover at microevolutionary scales is critical to assess the pace of this evolutionary process. Here, we characterized and quantified gene turnover for the epidemic lineage of a bacterial plant pathogen of major agricultural importance worldwide. Relying on a dense geographic sampling spanning 39 years of evolution, we estimated both the dynamics of single nucleotide polymorphism accumulation and gene content turnover. We identified extensive gene content variation among lineages even at the smallest phylogenetic and geographic scales. Gene turnover rate exceeded nucleotide substitution rate by three orders of magnitude. Accessory genes were found preferentially located on plasmids, but we identified a highly plastic chromosomal region hosting ecologically important genes such as transcription activator-like effectors. Whereas most changes in the gene content are probably transient, the rapid spread of a mobile element conferring resistance to copper compounds widely used for the management of plant bacterial pathogens illustrates how some accessory genes can become ubiquitous within a population over short timeframes.

Molecular Epidemiology of Acinetobacter calcoaceticus-Acinetobacter baumannii Complex Isolated From Children at the Hospital Infantil de México Federico Gómez

The Acinetobacter calcoaceticus-baumannii (Acb) complex is regarded as a group of phenotypically indistinguishable opportunistic pathogens responsible for mainly causing hospital-acquired pneumonia and bacteremia. The aim of this study was to determine the frequency of isolation of the species that constitute the Acb complex, as well as their susceptibility to antibiotics, and their distribution at the Hospital Infantil de Mexico Federico Gomez (HIMFG). A total of 88 strains previously identified by Vitek 2®, 40 as Acinetobacter baumannii and 48 as Acb complex were isolated from 52 children from 07, January 2015 to 28, September 2017. A. baumannii accounted for 89.77% (79/88) of the strains; Acinetobacter pittii, 6.82% (6/88); and Acinetobacter nosocomialis, 3.40% (3/88). Most strains were recovered mainly from patients in the intensive care unit (ICU) and emergency wards. Blood cultures (BC) provided 44.32% (39/88) of strains. The 13.63% (12/88) of strains were associated with primary bacteremia, 3.4% (3/88) with secondary bacteremia, and 2.3% (2/88) with pneumonia. In addition, 44.32% (39/88) were multidrug-resistant (MDR) strains and, 11.36% (10/88) were extensively drug-resistant (XDR). All strains amplified the bla_OXA-51 gene; 51.13% (45/88), the bla_OXA-23 gene; 4.54% (4/88), the bla_OXA-24 gene; and 2.27% (2/88), the bla_OXA-58 gene. Plasmid profiles showed that the strains had 1–6 plasmids. The strains were distributed in 52 pulsotypes, and 24 showed identical restriction patterns, with a correlation coefficient of 1.0. Notably, some strains with the same pulsotype were isolated from different patients, wards, or years, suggesting the persistence of more than one clone. Twenty-seven sequence types (STs) were determined for the strains based on a Pasteur multilocus sequence typing (MLST) scheme using massive sequencing; the most prevalent was ST 156 (27.27%, 24/88). The Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-Cas I-Fb system provided amplification in A. baumannii and A. pittii strains (22.73%, 20/88). This study identified an increased number of MDR strains and the relationship among strains through molecular typing. The data suggest that more than one strain could be causing an infection in some patient. The implementation of molecular epidemiology allowed the characterization of a set of strains and identification of different attributes associated with its distribution in a specific environment.

High mortality in an outbreak of multidrug resistant Acinetobacter baumannii infection introduced to an oncological hospital by a patient transferred from a general hospital

OBJECTIVE

To describe the clinical features, outcomes, and molecular epidemiology of an outbreak of multidrug resistant (MDR) A. baumannii.

METHODS

We performed a retrospective analysis of all MDR A. baumannii isolates recovered during an outbreak from 2011 to 2015 in a tertiary care cancer hospital. Cases were classified as colonized or infected. We determined sequence types following the Bartual scheme and plasmid profiles.

RESULTS

There were 106 strains of A. baumannii isolated during the study period. Sixty-six (62.3%) were considered as infection and 40 (37.7%) as colonization. The index case, identified by molecular epidemiology, was a patient with a drain transferred from a hospital outside Mexico City. Ninety-eight additional cases had the same MultiLocus Sequence Typing (MLST) 758, of which 94 also had the same plasmid profile, two had an extra plasmid, and two had a different plasmid. The remaining seven isolates belonged to different MLSTs. Fifty-three patients (50%) died within 30 days of A. baumanniii isolation: 28 (20%) in colonized and 45 (68.2%) in those classified as infection (p<0.001). In multivariate regression analysis, clinical infection and patients with hematologic neoplasm, predicted 30-day mortality. The molecular epidemiology of this outbreak showed the threat posed by the introduction of MDR strains from other institutions in a hospital of immunosuppressed patients and highlights the importance of adhering to preventive measures, including contact isolation, when admitting patients with draining wounds who have been hospitalized in other institutions.

Genotyping and molecular characterization of clinical Acinetobacter baumannii isolates from a single hospital in Southwestern Iran

ACINETOBACTER BAUMANNII

(A. baumannii) is a pathogen responsible for nosocomial infections among the hospitalized patients. The aim of this study was to investigate genotyping and molecular characterization and to examine the biofilm formation ability of A. baumannii isolates. In total, 70 A. baumannii isolates were collected from patients admitted to Imam Khomeini Hospital in Ahvaz, Southwestern Iran. Minimum inhibitory concentrations (MIC) test was performed using Vitek 2 system. The presence of genes encoding metallo-β-lactamases, oxacillinases, and integrase and the biofilm formation ability were then evaluated. Multiple locus variable-number tandem repeat (VNTR) analysis (MLVA) typing and multiplex PCR were performed to determine the genetic relationships. The blaOXA-23-like gene had the highest prevalence. The frequency of genes encoding blaSPM, blaIMP, and blaVIM among MDR A. baumannii isolates were 12 (17.1%), 18 (25.7%), and 22 (31.4%), respectively. Moreover, 46 isolates (75.4%) harbored class I integron and 10 isolates (16.39%) carried class II integron. The number of weak, moderate and strong biofilm-producing isolates were 3 (4.3%), 7 (10%), and 55 (78.5%), respectively. The results showed that 70 A. baumannii isolates were grouped into 12 distinct MLVA types with five clusters and four singleton genotypes. In addition, 25 (35.7%) isolates were assigned to international clone (IC) variants, 37 (52.8%) isolates belonged to group 1 (IC II), and 8  (11.4%) isolates belonged to group 2 (IC I). Our findings revealed that the population structure of the A. baumannii isolates was genetically diverse. More focus on genetic variation and antibiotic resistance of A. baumannii isolates are recommended.

Agricultural intensification and the evolution of host specialism in the enteric pathogen Campylobacter jejuni

Modern agriculture has dramatically changed the distribution of animal species on Earth. Changes to host ecology have a major impact on the microbiota, potentially increasing the risk of zoonotic pathogens being transmitted to humans, but the impact of intensive livestock production on host-associated bacteria has rarely been studied. Here, we use large isolate collections and comparative genomics techniques, linked to phenotype studies, to understand the timescale and genomic adaptations associated with the proliferation of the most common food-born bacterial pathogen (Campylobacter jejuni) in the most prolific agricultural mammal (cattle). Our findings reveal the emergence of cattle specialist C. jejuni lineages from a background of host generalist strains that coincided with the dramatic rise in cattle numbers in the 20th century. Cattle adaptation was associated with horizontal gene transfer and significant gene gain and loss. This may be related to differences in host diet, anatomy, and physiology, leading to the proliferation of globally disseminated cattle specialists of major public health importance. This work highlights how genomic plasticity can allow important zoonotic pathogens to exploit altered niches in the face of anthropogenic change and provides information for mitigating some of the risks posed by modern agricultural systems.

The Genomic Context for the Evolution and Transmission of Community-Associated Staphylococcus aureus ST59 Through the Food Chain

Sequence type 59 (ST59) is a predominant clonal lineage of community-acquired, methicillin-resistant Staphylococcus aureus (CA-MRSA) in Asia. Despite its increasing clinical relevance in China, the evolution and geographic expansion of ST59 has been relatively uncared for. Previous study has shown that ST59 was the predominant clone in food-related MRSA in China. This study compared the genomes of 87 clonal complex (CC) 59 S. aureus isolates sourced from food chain and infection cases to reconstruct the molecular evolution and geographical spread of ST59. Accordingly, three major sub-clades of ST59 were identified and these did not correlate with isolation source or location. Phylogenetic analysis estimated that ST59 in mainland China diverged from a most common recent ancestor around 1974, and most of the cases of cross-country transmission occurred between 1987 and 2000. Notably, two recent events of cross-country transmission through the food chain were observed, the isolates from these events diverged within relatively short time intervals. These isolates also showed high similarity in terms of their core genome, accessory genes, and antibiotic resistance patterns. These findings provide a valuable insight into the potential route of ST59 expansion in China and indicate a need for robust food chain surveillance to prevent the spread of this pathogen.

Phylogenomics Reveals Clear Cases of Misclassification and Genus-Wide Phylogenetic Markers for Acinetobacter

The Gram-negative Acinetobacter genus has several species of clear medical relevance. Many fully sequenced genomes belonging to the genus have been published in recent years; however, there has not been a recent attempt to infer the evolutionary history of Acinetobacter with that vast amount of information. Here, through a phylogenomic approach, we established the most up-to-date view of the evolutionary relationships within this genus and highlighted several cases of poor classification, especially for the very closely related species within the Acinetobacter calcoaceticus-Acinetobacter baumannii complex (Acb complex). Furthermore, we determined appropriate phylogenetic markers for this genus and showed that concatenation of the top 13 gives a very decent reflection of the evolutionary relationships for the genus Acinetobacter. The intersection between our top markers and previously defined universal markers is very small. In general, our study shows that, although there seems to be hardly any universal markers, bespoke phylogenomic approaches can be used to infer the phylogeny of different bacterial genera. We expect that ad hoc phylogenomic approaches will be the standard in the years to come and will provide enough information to resolve intricate evolutionary relationships like those observed in the Acb complex.

Origin of OXA-23 Variant OXA-239 from a Recently Emerged Lineage of Acinetobacter baumannii International Clone V

A. baumannii is a major cause of nosocomial infections all over the world. Although many isolates from developed countries have been studied in terms of their genome sequence, isolates from Latin America have been much less studied. In this study, using a population genomics approach considering the whole genomes of 148 isolates, we describe the recent emergence of the lineage ST758 endemic to Latin America and the inception of the OXA-239 carbapenemase. Our study highlights the urgent need to investigate recently emerged lineages of this species in Latin America and elsewhere, as these might harbor novel antibiotic resistance genes.

Over the last few decades, carbapenemase-producing Acinetobacter baumannii has become a major cause of nosocomial infections all over the world. However, the genome identity of lineages of this species in Latin America has not been studied as much as in developed countries. Here, through a population genomics approach considering the whole genomes of 148 isolates (almost 40 from Mexico and Honduras), we describe the recent emergence of the lineage sequence type 758 (ST758), which belongs to the international clone V and has spread out to Canada, Mexico, Honduras, and Colombia. Notably, this lineage was found to coexist with other A. baumannii lineages in hospitals in Mexico and Honduras. Isolates from this lineage show considerable variation in antibiotic resistance profiles, but most of them are resistant to carbapenems. Moreover, we found a variety of acquired oxacillinase (OXA) families within this lineage and tracked the very recent inception, and subsequent horizontal transmission, of the OXA-239 carbapenemase. This work highlights the urgent need to investigate recently emerged lineages of this species in Latin America and elsewhere, as these might harbor novel antibiotic resistance genes.

IMPORTANCEA. baumannii is a major cause of nosocomial infections all over the world. Although many isolates from developed countries have been studied in terms of their genome sequence, isolates from Latin America have been much less studied. In this study, using a population genomics approach considering the whole genomes of 148 isolates, we describe the recent emergence of the lineage ST758 endemic to Latin America and the inception of the OXA-239 carbapenemase. Our study highlights the urgent need to investigate recently emerged lineages of this species in Latin America and elsewhere, as these might harbor novel antibiotic resistance genes.

Risk Factors for Multidrug-Resistant Acinetobacter baumannii Infections in a Mass Burn Casualty Incident

They investigated the risk factors of multidrug-resistant Acinetobacter baumannii infections in 33 burn patients in a mass burn casualty. The independent variables included sex, burn size, vancomycin, ampicillin/sulbactam, cefazolin, meropenem, third-/fourth-generation cephalosporin, the number of classes of antibiotic used, and the number of days of antibiotic use. Multidrug-resistant A. baumannii infection was the outcome variable. Logistic regressions and structural equation model were used for statistical analysis. The average age was 21.7 years (range, 17-32 years; M = 11 [33%], F = 22 [67%]; mean burn area, 42% of the total body surface area; study period, June 28, 2015 to July 31, 2015). The incidence of infection was 39% (n = 13/33). For every percent increase in burn size, the odds of developing A. baumannii infections increased by 1.1 (P < .05). Similarly, the odds increased by 2.5 in every number increase in the classes of antibiotic used, by 1.2 in everyday increase in the days of antibiotic used, and by 12 in patients treated with vancomycin (P < .05). The correlations of these risk factors were demonstrated in a hypothetical structural equation model (P-value of chi-squared test > .05 and root mean square error of approximation < 0.05) in which burn size was the fundamental risk factor of multidrug-resistant A. baumannii infections. The model did not predict the emergence of multidrug-resistant A. baumannii infections. Conclusively, the risks of multidrug-resistant A. baumannii infections in burns are correlated with burn size, the number of classes of antibiotic used, the number of days of antibiotic use, and the prior use of vancomycin.

Carbapenem-Resistant Acinetobacter baumannii in Three Tertiary Care Hospitals in Mexico: Virulence Profiles, Innate Immune Response and Clonal Dissemination

Acinetobacter baumannii is one of the most important nosocomial pathogens distributed worldwide. Due to its multidrug-resistance and the propensity for the epidemic spread, the World Health Organization includes this bacterium as a priority health issue for development of new antibiotics. The aims of this study were to investigate the antimicrobial resistance profile, the clonal relatedness, the virulence profiles, the innate host immune response and the clonal dissemination of A. baumannii in Hospital Civil de Guadalajara (HCG), Hospital Regional General Ignacio Zaragoza (HRGIZ) and Pediatric ward of the Hospital General de México Eduardo Liceaga (HGM-P). A total of 252 A. baumannii clinical isolates were collected from patients with nosocomial infections in these hospitals between 2015 and 2016. These isolates showed a multidrug-resistant profile and most of them only susceptible to colistin. Furthermore, 83.3 and 36.9% of the isolates carried the bla_OXA–24 and bla_TEM–1 genes for resistance to carbapenems and β-lactam antibiotics, respectively. The clonal relatedness assessed by pulsed-field gel electrophoresis (PFGE) and by multi-locus sequence typing (MLST) demonstrated a genetic diversity. Remarkably, the ST136, ST208 and ST369 that belonged to the clonal complex CC92 and ST758 and ST1054 to the CC636 clonal complex were identified. The ST136 was a high-risk persistent clone involved in an outbreak at HCG and ST369 were related to the first carbapenem-resistant A. baumannii outbreak in HRGIZ. Up to 58% isolates were able to attach to A549 epithelial cells and 14.5% of them induced >50% of cytotoxicity. A549 cells infected with A. baumannii produced TNFα, IL-6 and IL-1β and the oxygen and nitrogen reactive species that contributes to the development of an inflammatory immune response. Up to 91.3% of clinical isolates were resistant to normal human serum activity. Finally, 98.5% of the clinical isolates were able to form biofilm over polystyrene tubes. In summary, these results demonstrate the increasingly dissemination of multidrug-resistant A. baumannii clones in three hospitals in Mexico carrying diverse bacterial virulence factors that could contribute to establishment of the innate immune response associated to the fatality risks in seriously ill patients.

Remarkable Genome Stability among emm1 Group A Streptococcus in Belgium over 19 Years

During the last two decades, there has been a public health concern of severe invasive infections caused by Group A Streptococcus (GAS) of the emm1 genotype. This study investigated the dynamics of emm1 GAS during 1994–2013 in Belgium. emm1 GAS isolated from blood, tissue, and wounds of patients with invasive infections (n = 23, S1–S23), and from patients with uncomplicated pharyngitis (n = 15, NS1–NS15) were subjected to whole-genome mapping (WGM; kpn) (Opgen). Whole-genome sequencing was performed on 25 strains (WGS; S1–S23 and NS6–NS7) (Illumina Inc.). Belgian GAS belonged to the M1T1 clone typified by the 36-kb chromosomal region encoding extracellular toxins, NAD+-glycohydrolase and streptolysin O. Strains from 1994–1999 clustered together with published strains (MGAS5005 and M1476). From 2001 onward, invasive GAS showed higher genomic divergence in the accessory genome and harbored on average 7% prophage content. Low evolutionary rate (2.49E-008; P > 0.05) was observed in this study, indicating a highly stable genome. The studied invasive and pharyngitis isolates were no genetically distinct populations based on the WGM and core genome phylogeny analyses. Two copies of the speJ superantigen were present in the 1999 and 2010 study strains (n = 3), one being chromosomal and one being truncated and associated with phage remnants.

This study showed that emm1 GAS in Belgium, compared with Canada and UK M1 strains, were highly conserved by harboring a remarkable genome stability over a 19-year period with variations observed in the accessory genome.

Update on the epidemiological typing methods for Acinetobacter baumannii

The outstanding ability of Acinetobacter baumannii to cause outbreaks and acquire multidrug resistance motivated the development of a plethora of typing techniques, which can help infection preventionists and hospital epidemiologists to more efficiently implement intervention controls. Nowadays, the world is witnessing a gradual transition from traditional typing methodology to whole genome sequencing-based approaches. Such approaches are opening new prospects and applications never achieved by existing typing methods. Herein, we provide the reader with an updated review on A. baumannii typing methods recapping the added value of well-established techniques previously applied for A. baumannii and detailing new ones (as clustered regularly interspaced short palindromic repeats-based typing) with a special focus on whole genome sequencing.

Inaccurate Multilocus Sequence Typing of Acinetobacter baumannii

Multilocus sequence typing has been useful for genotyping pathogens in surveillance and epidemiologic studies. However, it cannot reflect the true relationships of isolates for species with very dynamic genomes. Using a robust genome phylogeny, we demonstrated the limitations of this method for typing Acinetobacter baumannii.

Genome Analysis of Shigella flexneri Serotype 3b Strain SFL1520 Reveals Significant Horizontal Gene Acquisitions Including a Multidrug Resistance Cassette

Shigella flexneri is a major etiological agent of shigellosis in developing countries, primarily occurring in children under 5 years of age. We have sequenced, for the first time, the complete genome of S. flexneri serotype 3b (strain SFL1520). We used a hybrid sequencing method––both long-read MinION Flow (Oxford Nanopore Technologies) and short-read MiSeq (Illumina) sequencing to generate a high-quality reference genome. The SFL1520 chromosome was found to be ∼4.58 Mb long, with 4,729 coding sequences. Despite sharing a substantial number of genes with other publicly available S. flexneri genomes (2,803), the SFL1520 strain contains 1,926 accessory genes. The phage-related genes accounted for 8% of the SFL1520 genome, including remnants of the Sf6 bacteriophage with an intact O-acetyltransferase gene specific to serotype 3b. The SFL1520 chromosome was also found to contain a multiple-antibiotic resistance cassette conferring resistance to ampicillin, chloramphenicol, streptomycin, and tetracycline, which was potentially acquired from a plasmid via transposases. The phylogenetic analysis based on core genes showed a high level of similarity of SFL1520 with other S. flexneri serotypes; however, there were marked differences in the accessory genes of SFL1520. In particular, a large number of unique genes were identified in SFL1520 suggesting significant horizontal gene acquisition in a relatively short time period. The major virulence traits of SFL1520 (such as serotype conversion and antimicrobial resistance) were associated with horizontal gene acquisitions highlighting the role of horizontal gene transfer in S. flexneri diversity and evolution.

Rapid Replacement of Acinetobacter baumannii Strains Accompanied by Changes in Lipooligosaccharide Loci and Resistance Gene Repertoire

The population structure of health care-associated pathogens reflects patterns of diversification, selection, and dispersal over time. Empirical data detailing the long-term population dynamics of nosocomial pathogens provide information about how pathogens adapt in the face of exposure to diverse antimicrobial agents and other host and environmental pressures and can inform infection control priorities. Extensive sequencing of clinical isolates from one hospital spanning a decade and a second hospital in the Cleveland, OH, metropolitan area over a 3-year time period provided high-resolution genomic analysis of the Acinetobacter baumannii metapopulation. Genomic analysis demonstrated an almost complete replacement of the predominant strain groups with a new, genetically distinct strain group during the study period. The new group, termed clade F, differs from other global clone 2 (GC2) strains of A. baumannii in several ways, including its antibiotic resistance and lipooligosaccharide biosynthesis genes. Clade F strains are part of a large phylogenetic group with broad geographic representation. Phylogenetic analysis of single-nucleotide variants in core genome regions showed that although the Cleveland strains are phylogenetically distinct from those isolated from other locations, extensive intermixing of strains from the two hospital systems was apparent, suggesting either substantial exchange of strains or a shared, but geographically restricted, external pool from which infectious isolates were drawn. These findings document the rapid evolution of A. baumannii strains in two hospitals, with replacement of the predominant clade by a new clade with altered lipooligosaccharide loci and resistance gene repertoires.IMPORTANCE Multidrug-resistant (MDR) A. baumannii is a difficult-to-treat health care-associated pathogen. Knowing the resistance genes present in isolates causing infection aids in empirical treatment selection. Furthermore, knowledge of the genetic background can assist in tracking patterns of transmission to limit the spread of infections in hospitals. The appearance of a new genetic background in A. baumannii strains with a different set of resistance genes and cell surface structures suggests that strong selective pressures exist, even in highly MDR pathogens. Because the new strains have levels of antimicrobial resistance similar to those of the strains that were displaced, we hypothesize that other features, including host colonization and infection, may confer additional selective advantages and contribute to their increased prevalence.

Whole-Genome Sequences of Five Acinetobacter baumannii Strains From a Child With Leukemia M2

Acinetobacter baumannii is an opportunistic pathogen and is one of the primary etiological agents of healthcare-associated infections (HAIs). A. baumannii infections are difficult to treat due to the intrinsic and acquired antibiotic resistance of strains of this bacterium, which frequently limits therapeutic options. In this study, five A. baumannii strains (810CP, 433H, 434H, 483H, and A-2), all of which were isolated from a child with leukemia M2, were characterized through antibiotic susceptibility profiling, the detection of genes encoding carbapenem hydrolyzing oxacillinases, pulsed-field gel electrophoresis (PFGE), multilocus sequence typing (MLST), adherence and invasion assays toward the A549 cell line, and the whole-genome sequence (WGS). The five strains showed Multidrug resistant (MDR) profiles and amplification of the bla_OXA-23 gene, belonging to ST758 and grouped into two PFGE clusters. WGS of 810CP revealed the presence of a circular chromosome and two small plasmids, pAba810CPa and pAba810CPb. Both plasmids carried genes encoding the Sp1TA system, although resistance genes were not identified. A gene-by-gene comparison analysis was performed among the A. baumannii strains isolated in this study and others A. baumannii ST758 strains (HIMFG and INCan), showing that 86% of genes were present in all analyzed strains. Interestingly, the 433H, 434H, and 483H strains varied by 8–10 single-nucleotide variants (SNVs), while the A2 and 810CP strains varied by 46 SNVs. Subsequently, an analysis using BacWGSTdb showed that all of our strains had the same resistance genes and were ST758. However, some variations were observed in relation to virulence genes, mainly in the 810CP strain. The genes involved in the synthesis of hepta-acylated lipooligosaccharides, the pgaABCD locus encoding poly-β-1-6-N-acetylglucosamine, the ompA gene, Csu pili, bap, the two-component system bfms/bfmR, a member of the phospholipase D family, and two iron-uptake systems were identified in our A. baumannii strains genome. The five A. baumannii strains isolated from the child were genetically different and showed important characteristics that promote survival in a hospital environment. The elucidation of their genomic sequences provides important information for understanding their epidemiology, antibiotic resistance, and putative virulence factors.

Genomic Analysis of Consecutive Acinetobacter baumannii Strains From a Single Patient

Acinetobacter baumannii is one of the most important nosocomial pathogens, and thus it is required to investigate how it disseminate in hospitals and infect patients. We performed whole genome sequencing for 24 A. baumannii strains isolated successively from the blood of a single patient to evaluate whether repeated infections were due to re-infection or relapse infection and to investigate within-host evolution. The whole genome of the first strain, BL1, was sequenced de novo using the PacBio RSII system. BL2-BL24, were sequenced with an Illumina Hiseq4000 and mapped to the genome sequences of BL1. We identified 42 single-nucleotide variations among the strains. The SNVs differentiated the strains into three groups, BL1, BL2-BL16, and BL17-BL24, indicating that the patient suffered from re-infections or co-infections by similar, but different strains. The results also showed that A. baumannii strains in each group were rather stable at the genomic level. Our study emphasizes the importance of intensive infection control.

Virulent Epidemic Pneumonia in Sheep Caused by the Human Pathogen Acinetobacter baumannii

The human pathogen Acinetobacter baumannii has emerged as a frequent cause of hospital-acquired infections, but infection of animals has rarely been observed. Here we analyzed an outbreak of epidemic pneumonia killing hundreds of sheep on a farm in Pakistan and identified A. baumannii as the infecting agent. A pure culture of strain AbPK1 isolated from lungs of sick animals was inoculated into healthy sheep, which subsequently developed similar disease symptoms. Bacteria re-isolated from the infected animals were shown to be identical to the inoculum, fulfilling Koch’s postulates. Comparison of the AbPK1 genome against 2283 A. baumannii genomes from the NCBI database revealed that AbPK1 carries genes for unusual surface structures, including a unique composition of iron acquisition genes, genes for O-antigen synthesis and sialic acid-specific acetylases of cell-surface carbohydrates that could enable immune evasion. Several of these unusual and otherwise rarely present genes were also identified in genomes of phylogenetically unrelated A. baumannii isolates from combat-wounded US military from Afghanistan indicating a common gene pool in this geographical region. Based on core genome MLST this virulent isolate represents a newly emerging lineage of Global Clone 2, suggesting a human source for this disease outbreak. The observed epidemic, direct transmission from sheep to sheep, which is highly unusual for A. baumannii, has important consequences for human and animal health. First, direct animal-to-animal transmission facilitates fast spread of pathogen and disease in the flock. Second, it may establish a stable ecological niche and subsequent spread in a new host. And third, it constitutes a serious risk of transmission of this hyper-virulent clone from sheep back to humans, which may result in emergence of contagious disease amongst humans.

Phylogenomics picks out the par excellence markers for species phylogeny in the genus Staphylococcus

Although genome sequencing has become a very promising approach to conduct microbial taxonomy, few labs have the resources to afford this especially when dealing with data sets of hundreds to thousands of isolates. The goal of this study was to identify the most adequate loci for inferring the phylogeny of the species within the genus Staphylococcus; with the idea that those who cannot afford whole genome sequencing can use these loci to carry out species assignation confidently. We retrieved 177 orthologous groups (OGs) by using a genome-based phylogeny and an average nucleotide identity analysis. The top 26 OGs showed topologies similar to the species tree and the concatenation of them yielded a topology almost identical to that of the species tree. Furthermore, a phylogeny of just the top seven OGs could be used for species assignment. We sequenced four staphylococcus isolates to test the 26 OGs and found that these OGs were far superior to commonly used markers for this genus. On the whole, our procedure allowed identification of the most adequate markers for inferring the phylogeny within the genus Staphylococcus. We anticipate that this approach will be employed for the identification of the most suitable markers for other bacterial genera and can be very helpful to sort out poorly classified genera.

In Silico Genome-Wide Analysis Reveals the Potential Links Between Core Genome of Acidithiobacillus thiooxidans and Its Autotrophic Lifestyle

The coinage “pan-genome” was first introduced dating back to 2005, and was used to elaborate the entire gene repertoire of any given species. Core genome consists of genes shared by all bacterial strains studied and is considered to encode essential functions associated with species’ basic biology and phenotypes, yet its relatedness with bacterial lifestyle of the species remains elusive. We performed the pan-genome analysis of sulfur-oxidizing acidophile Acidithiobacillus thiooxidans as a case study to highlight species’ core genome and its relevance with autotrophic lifestyle of bacterial species. The mathematical modeling based on bacterial genomes of A. thiooxidans species, including a novel strain ZBY isolated from Zambian copper mine plus eight other recognized strains, was attempted to extrapolate the expansion of its pan-genome, suggesting that A. thiooxidans pan-genome is closed. Further investigation revealed a common set of genes, many of which were assigned to metabolic profiles, notably with respect to energy metabolism, amino acid metabolism, and carbohydrate metabolism. The predicted metabolic profiles of A. thiooxidans were characterized by the fixation of inorganic carbon, assimilation of nitrogen compounds, and aerobic oxidation of various sulfur species. Notably, several hydrogenase (H₂ase)-like genes dispersed in core genome might represent the novel classes due to the potential functional disparities, despite being closely related homologous genes that code for H₂ase. Overall, the findings shed light on the distinguishing features of A. thiooxidans genomes on a global scale, and extend the understanding of its conserved core genome pertaining to autotrophic lifestyle.

Evolution of the Quorum network and the mobilome (plasmids and bacteriophages) in clinical strains of Acinetobacter baumannii during a decade

In this study, we compared eighteen clinical strains of A. baumannii belonging to the ST-2 clone and isolated from patients in the same intensive care unit (ICU) in 2000 (9 strains referred to collectively as Ab_GEIH-2000) and 2010 (9 strains referred to collectively as Ab_GEIH-2010), during the GEIH-REIPI project (Umbrella BioProject PRJNA422585). We observed two main molecular differences between the Ab_GEIH-2010 and the Ab_GEIH-2000 collections, acquired over the course of the decade long sampling interval and involving the mobilome: i) a plasmid harbouring genes for bla_{OXA 24/40} ß-lactamase and abKA/abkB proteins of a toxin-antitoxin system; and ii) two temperate bacteriophages, Ab105-1ϕ (63 proteins) and Ab105-2ϕ (93 proteins), containing important viral defence proteins. Moreover, all Ab_GEIH-2010 strains contained a Quorum functional network of Quorum Sensing (QS) and Quorum Quenching (QQ) mechanisms, including a new QQ enzyme, AidA, which acts as a bacterial defence mechanism against the exogenous 3-oxo-C12-HSL. Interestingly, the infective capacity of the bacteriophages isolated in this study (Ab105-1ϕ and Ab105-2ϕ) was higher in the Ab_GEIH-2010 strains (carrying a functional Quorum network) than in the Ab_GEIH-2000 strains (carrying a deficient Quorum network), in which the bacteriophages showed little or no infectivity. This is the first study about the evolution of the Quorum network and the mobilome in clinical strains of Acinetobacter baumannii during a decade.