Accessory Genome Research Articles

Number of Candidate Effector Genes in Accessory Genomes Differentiates Pathogenic From Endophytic Fusarium oxysporum Strains.

The fungus Fusarium oxysporum (Fo) is widely known for causing wilt disease in over 100 different plant species. Endophytic interactions of Fo with plants are much more common, and strains pathogenic on one plant species can even be beneficial endophytes on another species. However, endophytic and beneficial interactions have been much less investigated at the molecular level, and the genetic basis that underlies endophytic versus pathogenic behavior is unknown. To investigate this, 44 Fo strains from non-cultivated Australian soils, grass roots from Spain, and tomato stems from United States were characterized genotypically by whole genome sequencing, and phenotypically by examining their ability to symptomlessly colonize tomato plants and to confer resistance against Fusarium Wilt. Comparison of the genomes of the validated endophytic Fo strains with those of 102 pathogenic strains revealed that both groups have similar genomes sizes, with similar amount of accessory DNA. However, although endophytic strains can harbor homologs of known effector genes, they have typically fewer effector gene candidates and associated non-autonomous transposons (mimps) than pathogenic strains. A pathogenic ‘lifestyle’ is associated with extended effector gene catalogs and a set of “host specific” effectors. No candidate effector genes unique to endophytic strains isolated from the same plant species were found, implying little or no host-specific adaptation. As plant-beneficial interactions were observed to be common for the tested Fo isolates, the propensity for endophytism and the ability to confer biocontrol appears to be a predominant feature of this organism. These findings allow prediction of the lifestyle of a Fo strain based on its genome sequence as a potential pathogen or as a harmless or even beneficial endophyte by determining its effectorome and mimp number.

Phylogenetic Analysis Reveals Distinct Evolutionary Trajectories of the Fluoroquinolones-Resistant Escherichia coli ST1193 From Fuzhou, China.

Escherichia coli (E. coli) ST1193 is an emerging fluoroquinolones-resistant and virulent lineage. Large gaps remain in our understanding of the evolutionary processes and differences of this lineage. Therefore, we used 76 E. coli ST1193 genomes to detect strain-level genetic diversity and phylogeny of this lineage globally. All E. coli ST1193 possessed fimH64, filCH5, and fumC14. There was 94.7% of isolates classified as O-type O75. There was 9.33% of E. coli ST1193 that possessed K5 capsular, while 90.67% of isolates possessed K1 capsular. The core genome analysis revealed that all isolates were divided into two phylogenetic clades (clade A and B). Clade A included 25 non-Chinese E. coli ST1193, and clade B contained all isolates collected from Fuzhou, China, respectively. The results of comparative genomics indicated Indels were identified in 150 clade-specific genes, which were enriched into the biological process and molecular function. Accessory genome phylogenetic tree showed a high degree of correlation between accessory genome clusters and core genome clades. There was significant difference in antibiotic resistance genes (ARGs) [blaCTX–M–55, blaTEM–1, sul2, tet(B), tet(R), APH(6)-Id, and AAC(3)-IId], virulence factors (cia, neuC, gad, and traT), and plasmid replicon types (IncQ1, Col156, and IncB/O/K/Z) between clade A (non-Chinese isolates) and clade B (Chinese isolates) (p < 0.05). Further analysis of the genetic environments of blaCTX–M–55 demonstrated that the flanking contexts of blaCTX–M–55 were diverse. In conclusion, our results reveal the distinct evolutionary trajectories of the spread of E. coli ST1193 in Fuzhou, China and non-China regions. This supports both global transmission and localized lineage expansion of this lineage following specific introductions into a geographic locality.

Determination of the two-component systems regulatory network reveals core and accessory regulations across Pseudomonas aeruginosa lineages.

Pseudomonas aeruginosa possesses one of the most complex bacterial regulatory networks, which largely contributes to its success as a pathogen. However, most of its transcription factors (TFs) are still uncharacterized and the potential intra-species variability in regulatory networks has been mostly ignored so far. Here, we used DAP-seq to map the genome-wide binding sites of all 55 DNA-binding two-component systems (TCSs) response regulators (RRs) across the three major P. aeruginosa lineages. The resulting networks encompass about 40% of all genes in each strain and contain numerous new regulatory interactions across most major physiological processes. Strikingly, about half of the detected targets are specific to only one or two strains, revealing a previously unknown large functional diversity of TFs within a single species. Three main mechanisms were found to drive this diversity, including differences in accessory genome content, as exemplified by the strain-specific plasmid in IHMA87 outlier strain which harbors numerous binding sites of conserved chromosomally-encoded RRs. Additionally, most RRs display potential auto-regulation or RR-RR cross-regulation, bringing to light the vast complexity of this network. Overall, we provide the first complete delineation of the TCSs regulatory network in P. aeruginosa that will represent an important resource for future studies on this pathogen.

Comparative Genomic Analysis of Streptococcus dysgalactiae subspecies dysgalactiae Isolated From Bovine Mastitis in China.

Streptococcus dysgalactiae subsp. dysgalactiae (SDSD) is one of the most prevalent pathogens causing bovine mastitis worldwide. However, there is a lack of comprehensive information regarding genetic diversity, complete profiles of virulence factors (VFs), and antimicrobial resistance (AMR) genes for SDSD associated with bovine mastitis in China. In this study, a total of 674 milk samples, including samples from 509 clinical and 165 subclinical mastitis cases, were collected from 17 herds in 7 provinces in China from November 2016 to June 2019. All SDSD isolates were included in phylogenetic analysis based on 16S rRNA and multi-locus sequence typing (MLST). In addition, whole genome sequencing was performed on 12 representative SDSD isolates to screen for VFs and AMR genes and to define pan-, core and accessory genomes. The prevalence of SDSD from mastitis milk samples was 7.57% (51/674). According to phylogenetic analysis based on 16S rRNA, 51 SDSD isolates were divided into 4 clusters, whereas based on MLST, 51 SDSD isolates were identified as 11 sequence types, including 6 registered STs and 5 novel STs (ST521, ST523, ST526, ST527, ST529) that belonged to 2 distinct clonal complexes (CCs) and 4 singletons. Based on WGS information, 108 VFs genes in 12 isolates were determined in 11 categories. In addition, 23 AMR genes were identified in 11 categories. Pan-, core and accessory genomes were composed of 2,663, 1,633 and 699 genes, respectively. These results provided a comprehensive profiles of SDSD virulence and resistance genes as well as phylogenetic relationships among mastitis associated SDSD in North China.

Comparative Genomics of Mycobacterium avium Complex Reveals Signatures of Environment-Specific Adaptation and Community Acquisition.

ABSTRACTNontuberculous mycobacteria, including those in the Mycobacterium avium complex (MAC), constitute an increasingly urgent threat to global public health. Ubiquitous in soil and water worldwide, MAC members cause a diverse array of infections in humans and animals that are often multidrug resistant, intractable, and deadly. MAC lung disease is of particular concern and is now more prevalent than tuberculosis in many countries, including the United States. Although the clinical importance of these microorganisms continues to expand, our understanding of their genomic diversity is limited, hampering basic and translational studies alike. Here, we leveraged a unique collection of genomes to characterize MAC population structure, gene content, and within-host strain dynamics in unprecedented detail. We found that different MAC species encode distinct suites of biomedically relevant genes, including antibiotic resistance genes and virulence factors, which may influence their distinct clinical manifestations. We observed that M. avium isolates from different sources—human pulmonary infections, human disseminated infections, animals, and natural environments—are readily distinguished by their core and accessory genomes, by their patterns of horizontal gene transfer, and by numerous specific genes, including virulence factors. We identified highly similar MAC strains from distinct patients within and across two geographically distinct clinical cohorts, providing important insights into the reservoirs which seed community acquisition. We also discovered a novel MAC genomospecies in one of these cohorts. Collectively, our results provide key genomic context for these emerging pathogens and will facilitate future exploration of MAC ecology, evolution, and pathogenesis.IMPORTANCE Members of the Mycobacterium avium complex (MAC), a group of mycobacteria encompassing M. avium and its closest relatives, are omnipresent in natural environments and emerging pathogens of humans and animals. MAC infections are difficult to treat, sometimes fatal, and increasingly common. Here, we used comparative genomics to illuminate key aspects of MAC biology. We found that different MAC species and M. avium isolates from different sources encode distinct suites of clinically relevant genes, including those for virulence and antibiotic resistance. We identified highly similar MAC strains in patients from different states and decades, suggesting community acquisition from dispersed and stable reservoirs, and we discovered a novel MAC species. Our work provides valuable insight into the genomic features underlying these versatile pathogens.

Application of Whole Genome Sequencing to Understand Diversity and Presence of Genes Associated with Sanitizer Tolerance in Listeria monocytogenes from Produce Handling Sources.

Recent listeriosis outbreaks linked to fresh produce suggest the need to better understand and mitigate L. monocytogenes contamination in packing and processing environments. Using whole genome sequencing (WGS) and phenotype screening assays for sanitizer tolerance, we characterized 48 L. monocytogenes isolates previously recovered from environmental samples in five produce handling facilities. Within the studied population there were 10 sequence types (STs) and 16 cgMLST types (CTs). Pairwise single nucleotide polymorphisms (SNPs) ranged from 0 to 3047 SNPs within a CT, revealing closely and distantly related isolates indicative of both sporadic and continuous contamination events within the facility. Within Facility 1, we identified a closely related cluster (0–2 SNPs) of isolates belonging to clonal complex 37 (CC37; CT9492), with isolates recovered during sampling events 1-year apart and in various locations inside and outside the facility. The accessory genome of these CC37 isolates varied from 94 to 210 genes. Notable genetic elements and mutations amongst the isolates included the bcrABC cassette (2/48), associated with QAC tolerance; mutations in the actA gene on the Listeria pathogenicity island (LIPI) 1 (20/48); presence of LIPI-3 (21/48) and LIPI-4 (23/48). This work highlights the potential use of WGS in tracing the pathogen within a facility and understanding properties of L. monocytogenes in produce settings.

PATO: Pangenome Analysis Toolkit.

We present the Pangenome Analysis Toolkit (PATO) designed to simultaneously analyze thousands of genomes using a desktop computer. The tool performs common tasks of pangenome analysis such as core-genome definition and accessory genome properties and includes new features that help characterize population structure, annotate pathogenic features and create gene sharedness networks. PATO has been developed in R to integrate with the large set of tools available for genetic, phylogenetic and statistical analysis in this environment. PATO can perform the most demanding bioinformatic analyses in minutes with an accuracy comparable to state-of-the-art software but 20-30× times faster. PATO also integrates all the necessary functions for the complete analysis of the most common objectives in microbiology studies. Finally, PATO includes the necessary tools for visualizing the results and can be integrated with other analytical packages available in R. The source code for PATO is freely available at https://github.com/irycisBioinfo/PATO under the GPLv3 license. Supplementary data are available at Bioinformatics online.

Expanding Characterized Diversity and the Pool of Complete Genome Sequences of Methylococcus Species, the Bacteria of High Environmental and Biotechnological Relevance.

The bacterial genus Methylococcus, which comprises aerobic thermotolerant methanotrophic cocci, was described half-a-century ago. Over the years, a member of this genus, Methylococcus capsulatus Bath, has become a major model organism to study genomic and metabolic basis of obligate methanotrophy. High biotechnological potential of fast-growing Methylococcus species, mainly as a promising source of feed protein, has also been recognized. Despite this big research attention, the currently cultured Methylococcus diversity is represented by members of the two species, M. capsulatus and M. geothermalis, while finished genome sequences are available only for two strains of these methanotrophs. This study extends the pool of phenotypically characterized Methylococcus strains with good-quality genome sequences by contributing four novel isolates of these bacteria from activated sludge, landfill cover soil, and freshwater sediments. The determined genome sizes of novel isolates varied between 3.2 and 4.0Mb. As revealed by the phylogenomic analysis, strains IO1, BH, and KN2 affiliate with M. capsulatus, while strain Mc7 may potentially represent a novel species. Highest temperature optima (45–50°C) and highest growth rates in bioreactor cultures (up to 0.3h−1) were recorded for strains obtained from activated sludge. The comparative analysis of all complete genomes of Methylococcus species revealed 4,485 gene clusters. Of these, pan-genome core comprised 2,331 genes (on average 51.9% of each genome), with the accessory genome containing 846 and 1,308 genes in the shell and the cloud, respectively. Independently of the isolation source, all strains of M. capsulatus displayed surprisingly high genome synteny and a striking similarity in gene content. Strain Mc7 from a landfill cover soil differed from other isolates by the high content of mobile genetic elements in the genome and a number of genome-encoded features missing in M. capsulatus, such as sucrose biosynthesis and the ability to scavenge phosphorus and sulfur from the environment.

Phylogenetic systematics of Butyrivibrio and Pseudobutyrivibrio genomes illustrate vast taxonomic diversity, open genomes and an abundance of carbohydrate-active enzyme family isoforms.

Butyrivibrio and Pseudobutyrivibrio dominate in anaerobic gastrointestinal microbiomes, particularly the rumen, where they play a key role in harvesting dietary energy. Within these genera, five rumen species have been classified ( Butyrivibrio fibrisolvens , Butyrivibrio hungatei , Butyrivibrio proteoclasticus , Pseudobutyrivibrio ruminis and Pseudobutyrivibrio xylanivorans ) and more recently an additional Butyrivibrio sp. group was added. Given the recent increase in available genomes, we re-investigated the phylogenetic systematics and evolution of Butyrivibrio and Pseudobutyrivibrio . Across 71 genomes, we show using 16S rDNA and 40 gene marker phylogenetic trees that the current six species designations ( P. ruminis , P. xylanivorans , B. fibrisolvens , Butyrivibrio sp., B. hungatei and B. proteclasticus) are found. However, pangenome analysis showed vast genomic variation and a high abundance of accessory genes (91.50–99.34 %), compared with core genes (0.66–8.50 %), within these six taxonomic groups, suggesting incorrectly assigned taxonomy. Subsequent pangenome accessory genomes under varying core gene cut-offs (%) and average nucleotide identity (ANI) analysis suggest the existence of 42 species within 32 genera. Pangenome analysis of those that still group within B. fibrisolvens , B. hungatei and P. ruminis , based on revised ANI phylogeny, also showed possession of very open genomes, illustrating the diversity that exists even within these groups. All strains of both Butyrivibrio and Pseudobutyrivibrio also shared a broad range of clusters of orthologous genes (COGs) (870), indicating recent evolution from a common ancestor. We also demonstrate that the carbohydrate-active enzymes (CAZymes) predominantly belong to glycosyl hydrolase (GH)2, 3, 5, 13 and 43, with numerous within family isoforms apparent, likely facilitating metabolic plasticity and resilience under dietary perturbations. This study provides a major advancement in our functional and evolutionary understanding of these important anaerobic bacteria.

Heuristic and Hierarchical-Based Population Mining of Salmonella enterica Lineage I Pan-Genomes as a Platform to Enhance Food Safety

The recent incorporation of bacterial whole-genome sequencing (WGS) into Public Health laboratories has enhanced foodborne outbreak detection and source attribution. As a result, large volumes of publicly available datasets can be used to study the biology of foodborne pathogen populations at an unprecedented scale. To demonstrate the application of a heuristic and agnostic hierarchical population structure guided pan-genome enrichment analysis (PANGEA), we used populations of S. enterica lineage I to achieve two main objectives: (i) show how hierarchical population inquiry at different scales of resolution can enhance ecological and epidemiological inquiries; and (ii) identify population-specific inferable traits that could provide selective advantages in food production environments. Publicly available WGS data were obtained from NCBI database for three serovars of Salmonella enterica subsp. enterica lineage I (S. Typhimurium, S. Newport, and S. Infantis). Using the hierarchical genotypic classifications (Serovar, BAPS1, ST, cgMLST), datasets from each of the three serovars showed varying degrees of clonal structuring. When the accessory genome (PANGEA) was mapped onto these hierarchical structures, accessory loci could be linked with specific genotypes. A large heavy-metal resistance mobile element was found in the Monophasic ST34 lineage of S. Typhimurium, and laboratory testing showed that Monophasic isolates have on average a higher degree of copper resistance than the Biphasic ones. In S. Newport, an extra sugE gene copy was found among most isolates of the ST45 lineage, and laboratory testing of multiple isolates confirmed that isolates of S. Newport ST45 were on average less sensitive to the disinfectant cetylpyridimium chloride than non-ST45 isolates. Lastly, data-mining of the accessory genomic content of S. Infantis revealed two cryptic Ecotypes with distinct accessory genomic content and distinct ecological patterns. Poultry appears to be the major reservoir for Ecotype 1, and temporal analysis further suggested a recent ecological succession, with Ecotype 2 apparently being displaced by Ecotype 1. Altogether, the use of a heuristic hierarchical-based population structure analysis that includes bacterial pan-genomes (core and accessory genomes) can (1) improve genomic resolution for mapping populations and accessing epidemiological patterns; and (2) define lineage-specific informative loci that may be associated with survival in the food chain.

Therapeutic target identification via differential genome analysis of antibiotic resistant Shigella sonnei and inhibitor evaluation against a selected drug target.

Shigella sonnei has been implicated in bloody diarrhea (accompanied by abdominal pain and fever) and is an emerging pathogen of concern, especially in developing countries. The major means of transmission is the fecal-oral route while sexual transmission has also been reported. In children, the impact might be stunted growth due to life-threatening illness. Resistance has been reported in this species for several types of antibiotics. In this study, we retrieved the antibiotic-resistant labeled whole genome sequences of the species from the PATRIC database and performed a pan-genome analysis to filter out core genes. Antibiotic resistance was studied in the core, accessory and unique genome. Core genes were utilized as seed substance for essentiality analysis and drug candidate assignment. Product of the gene aroG, i.e. chorismate biosynthetic process 3-deoxy-7-phosphoheptulonate synthase enzyme, responsible for aromatic amino acid family biosynthetic process, was taken for further downstream processing. Natural product libraries of flavonoids (n=178), ZINC database derived inhibitor compounds of the 3-deoxy-7-phosphoheptulonate synthase enzyme (n=112), and streptomycin compounds (n=737) were docked to find out potent inhibitors, followed by dynamics simulation of 50ns each for top compounds.. Physicochemical and ADMET profiling of the top compounds was done to analyze their safety for consumption. We propose that the top compounds: Phytoene from Streptomycin library and ZINC000036444158 (synonym:1,16-bis[(dihydroxyphosphinyl)oxy]hexadecane) from 3-deoxy-7-phosphoheptulonate synthase inhibitor library of ZINC database (and used as a control in this study) should be tested in vitro against Shigella sonnei, to fully determine their efficacy. This could add to the drying pipeline of potent drug molecules against emerging pathogens.

Genomic characterization of Enterobacter xiangfangensis STP-3: Application to real time petroleum oil sludge bioremediation

Sustainable treatment of petroleum oil sludge still remains as a major challenge to petroleum refineries. Bioremediation is the promising technology involving bacteria for simultaneous production of biosurfactant and followed by degradation of petroleum compounds. Complete genomic knowledge on such potential microbes could accentuate its successful exploitation. The present study discusses the genomic characteristics of novel biosurfactant producing petrophilic/ petroleum hydrocarbon degrading strain, Enterobacter xiangfangensis STP-3, isolated from petroleum refinery oil sludge contaminated soil. The genome has 4,584,462 bp and 4372 protein coding sequences. Functional analysis using the RAST and KEGG databases revealed the presence of biosynthetic gene clusters linked to glycolipid and lipopeptide production and multiple key candidate genes linked with the degradation pathway of petroleum hydrocarbons. Orthology study revealed diversity in gene clusters associated to membrane transport, carbohydrate, amino acid metabolism, virulence and defence mechanisms, and nucleoside and nucleotide synthesis. The comparative analysis with 27 other genomes predicted that the core genome contributes to its inherent bioremediation potential, whereas the accessory genome influences its environmental adaptability in unconventional environmental conditions. Further, experimental results showed that E. xiangfangensis STP-3 was able to degrade PHCs by 82 % in 14 days during the bioremediation of real time petroleum oil sludge with the concomitant production of biosurfactant and metabolic enzymes, To the best of our knowledge, no comprehensive genomic study has been previously reported on the biotechnological prospective of this species.

Genome Sequence Analysis of Clostridium chauvoei Strains of European Origin and Evaluation of Typing Options for Outbreak Investigations.

Black quarter caused by Clostridium (C.) chauvoei is an important bacterial disease that affects cattle and sheep with high mortality. A comparative genomics analysis of 64 C. chauvoei strains, most of European origin and a few of non-European and unknown origin, was performed. The pangenome analysis showed limited new gene acquisition for the species. The accessory genome involved prophages and genomic islands, with variations in gene composition observed in a few strains. This limited accessory genome may indicate that the species replicates only in the host or that an active CRISPR/Cas system provides immunity to foreign genetic elements. All strains contained a CRISPR type I-B system and it was confirmed that the unique spacer sequences therein can be used to differentiate strains. Homologous recombination events, which may have contributed to the evolution of this pathogen, were less frequent compared to other related species from the genus. Pangenome single nucleotide polymorphism (SNP) based phylogeny and clustering indicate diverse clusters related to geographical origin. Interestingly the identified SNPs were mostly non-synonymous. The study demonstrates the possibility of the existence of polymorphic populations in one host, based on strain variability observed for strains from the same animal and strains from different animals of one outbreak. The study also demonstrates that new outbreak strains are mostly related to earlier outbreak strains from the same farm/region. This indicates the last common ancestor strain from one farm can be crucial to understand the genetic changes and epidemiology occurring at farm level. Known virulence factors for the species were highly conserved among the strains. Genetic elements involved in Nicotinamide adenine dinucleotide (NAD) precursor synthesis (via nadA, nadB, and nadC metabolic pathway) which are known as potential anti-virulence loci are completely absent in C. chauvoei compared to the partial inactivation in C. septicum. A novel core-genome MLST based typing method was compared to sequence typing based on CRISPR spacers to evaluate the usefulness of the methods for outbreak investigations.

Niche Preference of Escherichia coli in a Peri-Urban Pond Ecosystem.

Escherichia coli comprises diverse strains with a large accessory genome, indicating functional diversity and the ability to adapt to a range of niches. Specific strains would display greatest fitness in niches matching their combination of phenotypic traits. Given this hypothesis, we sought to determine whether E. coli in a peri-urban pond and associated cattle pasture display niche preference. Samples were collected from water, sediment, aquatic plants, water snails associated with the pond, as well as bovine feces from cattle in an adjacent pasture. Isolates (120) were obtained after plating on Membrane Lactose Glucuronide Agar (MLGA). We used the uidA and mutS sequences for all isolates to determine phylogeny by maximum likelihood, and population structure through gene flow analysis. PCR was used to allocate isolates to phylogroups and to determine the presence of pathogenicity/virulence genes (stxI, stxII, eaeA, hlyA, ST, and LT). Antimicrobial resistance was determined using a disk diffusion assay for Tetracycline, Gentamicin, Ciprofloxacin, Meropenem, Ceftriaxone, and Azithromycin. Our results showed that isolates from water, sediment, and water plants were similar by phylogroup distribution, virulence gene distribution, and antibiotic resistance while both snail and feces populations were significantly different. Few of the feces isolates were significantly similar to aquatic ones, and most of the snail isolates were also different. Population structure analysis indicated three genetic backgrounds associated with bovine, snail, and aquatic environments. Collectively these data support niche preference of E. coli isolates occurring in this ecosystem.

Complete Genome Sequencing of Leptospira interrogans Isolates from Malaysia Reveals Massive Genome Rearrangement but High Conservation of Virulence-Associated Genes.

The ability of Leptospirae to persist in environments and animal hosts but to cause clinically highly variable disease in humans has made leptospirosis the most common zoonotic disease. Considering the paucity of data on variation in complete genomes of human pathogenic Leptospirae, we have used a combination of Single Molecule Real-Time (SMRT) and Illumina sequencing to obtain complete genome sequences of six human clinical L. interrogans isolates from Malaysia. All six contained the larger (4.28–4.56 Mb) and smaller (0.34–0.395 Mb) chromosome typical of human pathogenic Leptospirae and 0–7 plasmids. Only 24% of the plasmid sequences could be matched to databases. We identified a chromosomal core genome of 3318 coding sequences and strain-specific accessory genomes of 49–179 coding sequences. These sequences enabled detailed genomic strain typing (Genome BLAST Distance Phylogeny, DNA–DNA hybridization, and multi locus sequence typing) and phylogenetic classification (whole-genome SNP genotyping). Even though there was some shared synteny and collinearity across the six genomes, there was evidence of major genome rearrangement, likely driven by horizontal gene transfer and homologous recombination. Mobile genetic elements were identified in all strains in highly varying numbers, including in the rfb locus, which defines serogroups and contributes to immune escape and pathogenesis. On the other hand, there was high conservation of virulence-associated genes including those relating to sialic acid, alginate, and lipid A biosynthesis. These findings suggest (i) that the antigenic variation, adaption to various host environments, and broad spectrum of virulence of L. interrogans are in part due to a high degree of genomic plasticity and (ii) that human pathogenic strains maintain a core set of genes required for virulence.

Pandora: nucleotide-resolution bacterial pan-genomics with reference graphs

We present pandora, a novel pan-genome graph structure and algorithms for identifying variants across the full bacterial pan-genome. As much bacterial adaptability hinges on the accessory genome, methods which analyze SNPs in just the core genome have unsatisfactory limitations. Pandora approximates a sequenced genome as a recombinant of references, detects novel variation and pan-genotypes multiple samples. Using a reference graph of 578 Escherichia coli genomes, we compare 20 diverse isolates. Pandora recovers more rare SNPs than single-reference-based tools, is significantly better than picking the closest RefSeq reference, and provides a stable framework for analyzing diverse samples without reference bias.

Whole-genome sequence-informed MALDI-TOF MS diagnostics reveal importance of Klebsiella oxytoca group in invasive infections: a retrospective clinical study

BackgroundKlebsiella spp. are opportunistic pathogens which can cause severe infections, are often multi-drug resistant and are a common cause of hospital-acquired infections. Multiple new Klebsiella species have recently been described, yet their clinical impact and antibiotic resistance profiles are largely unknown. We aimed to explore Klebsiella group- and species-specific clinical impact, antimicrobial resistance (AMR) and virulence.MethodsWe analysed whole-genome sequence data of a diverse selection of Klebsiella spp. isolates and identified resistance and virulence factors. Using the genomes of 3594 Klebsiella isolates, we predicted the masses of 56 ribosomal subunit proteins and identified species-specific marker masses. We then re-analysed over 22,000 Matrix-Assisted Laser Desorption Ionization - Time Of Flight (MALDI-TOF) mass spectra routinely acquired at eight healthcare institutions in four countries looking for these species-specific markers. Analyses of clinical and microbiological endpoints from a subset of 957 patients with infections from Klebsiella species were performed using generalized linear mixed-effects models.ResultsOur comparative genomic analysis shows group- and species-specific trends in accessory genome composition. With the identified species-specific marker masses, eight Klebsiella species can be distinguished using MALDI-TOF MS. We identified K. pneumoniae (71.2%; n = 12,523), K. quasipneumoniae (3.3%; n = 575), K. variicola (9.8%; n = 1717), “K. quasivariicola” (0.3%; n = 52), K. oxytoca (8.2%; n = 1445), K. michiganensis (4.8%; n = 836), K. grimontii (2.4%; n = 425) and K. huaxensis (0.1%; n = 12). Isolates belonging to the K. oxytoca group, which includes the species K. oxytoca, K. michiganensis and K. grimontii, were less often resistant to 4th-generation cephalosporins than isolates of the K. pneumoniae group, which includes the species K. pneumoniae, K. quasipneumoniae, K. variicola and “K. quasivariicola” (odds ratio = 0.17, p < 0.001, 95% confidence interval [0.09,0.28]). Within the K. pneumoniae group, isolates identified as K. pneumoniae were more often resistant to 4th-generation cephalosporins than K. variicola isolates (odds ratio = 2.61, p = 0.003, 95% confidence interval [1.38,5.06]). K. oxytoca group isolates were found to be more likely associated with invasive infection to primary sterile sites than K. pneumoniae group isolates (odds ratio = 2.39, p = 0.0044, 95% confidence interval [1.05,5.53]).ConclusionsCurrently misdiagnosed Klebsiella spp. can be distinguished using a ribosomal marker-based approach for MALDI-TOF MS. Klebsiella groups and species differed in AMR profiles, and in their association with invasive infection, highlighting the importance for species identification to enable effective treatment options.

Gene-gene relationships in an Escherichia coli accessory genome are linked to function and mobility.

The pangenome contains all genes encoded by a species, with the core genome present in all strains and the accessory genome in only a subset. Coincident gene relationships are expected within the accessory genome, where the presence or absence of one gene is influenced by the presence or absence of another. Here, we analysed the accessory genome of an Escherichia coli pangenome consisting of 400 genomes from 20 sequence types to identify genes that display significant co-occurrence or avoidance patterns with one another. We present a complex network of genes that are either found together or that avoid one another more often than would be expected by chance, and show that these relationships vary by lineage. We demonstrate that genes co-occur by function, and that several highly connected gene relationships are linked to mobile genetic elements. We find that genes are more likely to co-occur with, rather than avoid, another gene in the accessory genome. This work furthers our understanding of the dynamic nature of prokaryote pangenomes and implicates both function and mobility as drivers of gene relationships.

Comparative Genome Analysis and Spore Heat Resistance Assay Reveal a New Component to Population Structure and Genome Epidemiology Within Clostridium perfringens Enterotoxin-Carrying Isolates.

Clostridium perfringens causes a variety of human and animal enteric diseases including food poisoning, antibiotic-associated diarrhea, and necrotic enteritis. Yet, the reservoirs of enteropathogenic enterotoxin-producing strains remain unknown. We conducted a genomic comparison of 290 strains and a heat resistance phenotyping of 30 C. perfringens strains to elucidate the population structure and ecology of this pathogen. C. perfringens genomes shared a conserved genetic backbone with more than half of the genes of an average genome conserved in >95% of strains. The cpe-carrying isolates were found to share genetic context: the cpe-carrying plasmids had different distribution patterns within the genetic lineages and the estimated pan genome of cpe-carrying isolates had a larger core genome and a smaller accessory genome compared to that of 290 strains. We characterize cpe-negative strains related to chromosomal cpe-carrying strains elucidating the origin of these strains and disclose two distinct groups of chromosomal cpe-carrying strains with different virulence characteristics, spore heat resistance properties, and, presumably, ecological niche. Finally, an antibiotic-associated diarrhea isolate carrying two copies of the enterotoxin cpe gene and the associated genetic lineage with the potential for the emergence of similar strains are outlined. With C. perfringens as an example, implications of input genome quality for pan genome analysis are discussed. Our study furthers the understanding of genome epidemiology and population structure of enteropathogenic C. perfringens and brings new insight into this important pathogen and its reservoirs.

Multidrug-Resistant Lineage of Enterotoxigenic Escherichia coli ST182 With Serotype O169:H41 in Airline Waste.

Enterotoxigenic Escherichia coli (ETEC) is the primary aetiologic agent of traveller’s diarrhoea and a significant cause of diarrhoeal disease and death in developing countries. ETEC O169:H41 strains are known to cause both traveller’s diarrhoea and foodborne outbreaks in developed countries and are cause for concern. Here, whole-genome sequencing (WGS) was used to assemble 46 O169:H41 (ST182) E. coli draft genomes derived from two airplane waste samples sourced from a German international airport. The ST182 genomes were compared with all 84 publicly available, geographically diverse ST182 genomes to construct a core genome-based phylogenetic tree. ST182 isolates were all phylogroup E, the majority serotype O169:H41 (n = 121, 93%) and formed five major clades. The airplane waste isolates differed by an average of 15 core SNPs (range 0–45) but their accessory genome content was diverse. While uncommon in other ST182 genomes, all airplane-derived ST182 isolates carried: (i) extended-spectrum β-lactamase gene blaCTX–M–15 notably lacking the typical adjacent ISEcp1; (ii) qnrS1 and the S83L mutation in gyrA, both conferring resistance to fluoroquinolones; and (iii) a class 1 integron structure (IS26-intI1Δ648-dfrA17-aadA5-qacEΔ1-sul1-ORF-srpC-padR-IS6100-mphR-mrx-mphA-IS26) identified previously in major extraintestinal pathogenic E. coli STs but not in ETEC. ST182 isolates carried ETEC-specific virulence factors STp + CS6. Adhesin/invasin tia was identified in 89% of aircraft ST182 isolates (vs 23%) and was located on a putative genomic island within a hotspot region for various insertions including PAI I536 and plasmid-associated transposons. The most common plasmid replicons in this collection were IncFII (100%; F2:A-:B-) and IncB/O/K/Z (89%). Our data suggest that potentially through travel, E. coli ST182 are evolving a multidrug-resistant profile through the acquisition of class 1 integrons and different plasmids.