Single Nucleotide Polymorphism Phylogeny Research Articles

Phylogenomics of Acinetobacter species and analysis of antimicrobial resistance genes.

Non-baumannii Acinetobacter species are increasingly isolated in the clinical setting and the environment. The aim of the present study was to analyze a genome database of 837 Acinetobacter spp. isolates, which included 798 non-baumannii Acinetobacter genomes, in order to define the concordance of classification and discriminatory power of 7-gene MLST, 53-gene MLST, and single-nucleotide polymorphism (SNPs) phylogenies. Phylogenies were performed on Pasteur Multilocus Sequence Typing (MLST) or ribosomal Multilocus Sequence Typing (rMLST) concatenated alleles, or SNPs extracted from core genome alignment. The Pasteur MLST scheme was able to identify and genotype 72 species in the Acinetobacter genus, with classification results concordant with the ribosomal MLST scheme. The discriminatory power and genotyping reliability of the Pasteur MLST scheme were assessed in comparison to genome-wide SNP phylogeny on 535 non-baumannii Acinetobacter genomes assigned to Acinetobacter pittii, Acinetobacter nosocomialis, Acinetobacter seifertii, and Acinetobacter lactucae (heterotypic synonym of Acinetobacter dijkshoorniae), which were the most clinically relevant non-baumannii species of the A. baumannii group. The Pasteur MLST and SNP phylogenies were congruent at Robinson-Fould and Matching cluster tests and grouped genomes into four and three clusters in A. pittii, respectively, and one each in A. seifertii. Furthermore, A. lactucae genomes were grouped into one cluster within A. pittii genomes. The SNP phylogeny of A. nosocomialis genomes showed a heterogeneous population and did not correspond to the Pasteur MLST phylogeny, which identified two recombinant clusters. The antimicrobial resistance genes belonging to at least three different antimicrobial classes were identified in 91 isolates assigned to 17 distinct species in the Acinetobacter genus. Moreover, the presence of a class D oxacillinase, which is a naturally occurring enzyme in several Acinetobacter species, was found in 503 isolates assigned to 35 Acinetobacter species. In conclusion, Pasteur MLST phylogeny of non-baumannii Acinetobacter isolates coupled with in silico detection of antimicrobial resistance makes it important to study the population structure and epidemiology of Acinetobacter spp. isolates.

Precision metagenomics sequencing for food safety: hybrid assembly of Shiga toxin-producing Escherichia coli in enriched agricultural water.

Culture-independent metagenomic sequencing of enriched agricultural water could expedite the detection and virulotyping of Shiga toxin-producing Escherichia coli (STEC). We previously determined the limits of a complete, closed metagenome-assembled genome (MAG) assembly and of a complete, fragmented MAG assembly for O157:H7 in enriched agricultural water using long reads (Oxford Nanopore Technologies, Oxford), which were 107 and 105 CFU/ml, respectively. However, the nanopore assemblies did not have enough accuracy to be used in Single Nucleotide Polymorphism (SNP) phylogenies and cannot be used for the precise identification of an outbreak STEC strain. The present study aimed to determine the limits of detection and assembly for STECs in enriched agricultural water by Illumina MiSeq sequencing technology alone, followed by establishing the limit of hybrid assembly with nanopore long-read sequencing using three different hybrid assemblers (SPAdes, Unicycler, and OPERA-MS). We also aimed to generate a genome with enough accuracy to be used in a SNP phylogeny. The classification of MiSeq and nanopore sequencing identified the same highly abundant species. Using the totality of the MiSeq output and a precision metagenomics approach in which the E. coli reads are binned before assembly, the limit of detection and assembly of STECs by MiSeq were determined to be 105 and 107 CFU/ml, respectively. While a complete, closed MAG could not be generated at any concentration, a complete, fragmented MAG was produced using the SPAdes assembler with an STEC concentration of at least 107 CFU/ml. At this concentration, hybrid assembled contigs aligned to the nanopore-assembled genome could be accurately placed in a neighbor-joining tree. The MiSeq limit of detection and assembly was less sensitive than nanopore sequencing, which was likely due to factors including the small starting material (50 vs. 1 μg) and the dilution of the library loaded on the cartridge. This pilot study demonstrates that MiSeq sequencing requires higher coverage in precision metagenomic samples; however, with sufficient concentration, STECs can be characterized and phylogeny can be accurately determined.

High-Resolution Genomic Profiling of Carbapenem-Resistant Klebsiella pneumoniae Isolates: A Multicentric Retrospective Indian Study.

BackgroundCarbapenem-resistant Klebsiella pneumoniae (CRKP) is a threat to public health in India because of its high dissemination, mortality, and limited treatment options. Its genomic variability is reflected in the diversity of sequence types, virulence factors, and antimicrobial resistance (AMR) mechanisms. This study aims to characterize the clonal relationships and genetic mechanisms of resistance and virulence in CRKP isolates in India.Materials and MethodsWe characterized 344 retrospective K. pneumoniae clinical isolates collected from 8 centers across India collected in 2013–2019. Susceptibility to antibiotics was tested with VITEK 2. Capsular types, multilocus sequence type, virulence genes, AMR determinants, plasmid replicon types, and a single-nucleotide polymorphism phylogeny were inferred from their whole genome sequences.ResultsPhylogenetic analysis of the 325 Klebsiella isolates that passed quality control revealed 3 groups: K. pneumoniae sensu stricto (n = 307), K. quasipneumoniae (n = 17), and K. variicola (n = 1). Sequencing and capsular diversity analysis of the 307 K. pneumoniae sensu stricto isolates revealed 28 sequence types, 26 K-locus types, and 11 O-locus types, with ST231, KL51, and O1V2 being predominant. blaOXA-48-like and blaNDM-1/5 were present in 73.2% and 24.4% of isolates, respectively. The major plasmid replicon types associated with carbapenase genes were IncF (51.0%) and Col group (35.0%).ConclusionOur study documents for the first time the genetic diversity of K and O antigens circulating in India. The results demonstrate the practical applicability of genomic surveillance and its utility in tracking the population dynamics of CRKP. It alerts us to the urgency for longitudinal surveillance of these transmissible lineages.

Escherichia coli Causing Recurrent Urinary Tract Infections: Comparison to Non-Recurrent Isolates and Genomic Adaptation in Recurrent Infections.

Recurrent urinary tract infection (rUTI) remains a major problem for many women and therefore the pursuit for genomic and phenotypic traits which could define rUTI has been ongoing. The present study applied a genomic approach to investigate recurrent urinary tract infections by comparative analyses of recurrent and non-recurrent Escherichia coli isolates from general practice. From whole-genome sequencing data, phylogenetic clustering and genomic traits were studied on a collection of isolates which caused recurrent infection compared to non-recurrent isolates. In addition, genomic variation between the 1st and following infection was studied on a subset of the isolates. Evidence of limited adaptation between the recurrent infections based on single nucleotide polymorphism analyses with a range of 0–13 non-synonymous single nucleotide polymorphisms (SNPs) between the paired isolates. This included an overrepresentation of SNPs in metabolism genes. We identified several genes which were more common in rUTI isolates, including nine fimbrial genes, however, not significantly after false-discovery rate. Finally, the results show that recurrent isolates of the present dataset are not distinctive by variation in the core genome, and thus, did not cluster distinct from non-rUTI isolates in a SNP phylogeny.

Genomic Characterization of a New Biofilm-Forming and Adhesive ST398 Human-Adapted MSSA Lineage Causing Septic Knee Arthritis Following Surgical Reconstruction.

Methicillin-susceptible (MSSA) and methicillin-resistant Staphylococcus aureus (MRSA) is a pathogen commonly found in bone and joint infections, including septic arthritis. S. aureus virulence and the frailty of affected patients can cause several complications; a prompt and specific antibiotic treatment can positively affect the outcome of patients. We carried out an in-depth genomic characterization by Illumina whole genome sequencing and bioinformatics of two biofilm-producing M1 and M2 ST398 MSSA causing septic knee arthritis not-responding to antimicrobial therapy. The strains were characterized for antibiotic resistance, biofilm and adhesive properties as well as genomics, single nucleotide polymorphism phylogeny, resistomics and virulomics. Our results showed that M1 and M2 MSSA were ST398-t1451-agrI-Cap5, susceptible to cefoxitin and resistant to erythromycin and clindamycin, traits consistent with the lack of the SCCmec-locus and the presence of the sole blaZ and ermT. Furthermore, M1 and M2 were biofilm-producing and largely potentially adhesive strains, as indicated by the adhesion gene profile. Our data characterized a new human-adapted ST398 MSSA lineage, representing a “fusion” between the human-animal independent ST398 and the Livestock Associated (LA) ST398 lineages, forming biofilm and genomically predicted high adhesive, characterized by different genomic adaptation conferring a great ability to adhere to the host’s extracellular matrix causing septic knee arthritis.

Genetic Traces of the Francisella tularensis Colonization of Spain, 1998-2020.

More than 1000 humans have acquired the febrile disease tularemia in Spain since the first notification of human cases in 1997. We here aimed to study the recent molecular evolution of the causative bacterium Francisella tularensis during disease establishment in Spain. Single-nucleotide polymorphisms (SNPs) and variable-number tandem repeats (VNTRs) were analyzed in whole-genome sequences (WGS) of F. tularensis. Short-read WGS data for 20 F. tularensis strains from humans infected in the periods 2014–2015 and 2018–2020 in Spain were generated. These data were combined with WGS data of 25 Spanish strains from 1998 to 2008 and two reference strains. Capillary electrophoresis data of VNTR genetic regions were generated and compared with the WGS data for the 11 strains from 2014 to 2015. Evolutionary relationships among strains were analyzed by phylogenetic methods. We identified 117 informative SNPs in a 1,577,289-nucleotide WGS alignment of 47 F. tularensis genomes. Forty-five strains from Spain formed a star-like SNP phylogeny with six branches emerging from a basal common node. The most recently evolved genomes formed four additional star-like structures that were derived from four branches of the basal common node. VNTR copy number variation was detected in two out of 10 VNTR regions examined. Genetic clustering of strains by VNTRs agreed with the clustering by SNPs. The SNP data provided higher resolution among strains than the VNTRs data in all but one cases. There was an excellent correlation between VNTR marker sizing by capillary electrophoresis and prediction from WGS data. The genetic data strongly support that tularemia, indeed, emerged recently in Spain. Distinct genetic patterns of local F. tularensis population expansions imply that the pathogen has colonized a previously disease-free geographical area. We also found that genome-wide SNPs provide higher genetic resolution among F. tularensis genomes than the use of VNTRs, and that VNTR copy numbers can be accurately predicted using short-read WGS data.

Whole-genome sequence-based comparison and profiling of virulence-associated genes of Bacillus cereus group isolates from diverse sources in Japan

BackgroundThe complete genome sequences of 44 Bacillus cereus group isolates collected from diverse sources in Japan were analyzed to determine their genetic backgrounds and diversity levels in Japan. Multilocus sequence typing (MLST) and core-genome single-nucleotide polymorphism (SNP) typing data from whole-genome sequences were analyzed to determine genetic diversity levels. Virulence-associated gene profiles were also used to evaluate the genetic backgrounds and relationships among the isolates.ResultsThe 44 B. cereus group isolates, including soil- and animal-derived isolates and isolates recovered from hospitalized patients and food poisoning cases, were genotyped by MLST and core-genome SNP typing. Genetic variation among the isolates was identified by the MLST and core-genome SNP phylogeny comparison against reference strains from countries outside of Japan. Exploratory principal component analysis and nonmetric multidimensional scaling (NMDS) analyses were used to assess the genetic similarities among the isolates using gene presence and absence information and isolate origins as the metadata. A significant correlation was seen between the principal components and the presence of genes encoding hemolysin BL and emetic genetic determinants in B. cereus, and the capsule proteins in B. anthracis. NMDS showed that the cluster of soil isolates overlapped with the cluster comprising animal-derived and clinical isolates.ConclusionsMolecular and epidemiological analyses of B. cereus group isolates in Japan suggest that the soil- and animal-derived bacteria from our study are not a significant risk to human health. However, because several of the clinical isolates share close genetic relationships with the environmental isolates, both molecular and epidemiological surveillance studies could be used effectively to estimate virulence in these important pathogens.

Genetic diversity and delineation of Salmonella Agona outbreak strains by next generation sequencing, Bavaria, Germany, 1993 to 2018.

BackgroundIn 2017, a food-borne Salmonella Agona outbreak caused by infant milk products from a French supplier occurred in Europe. Simultaneously, S. Agona was detected in animal feed samples in Bavaria.AimUsing next generation sequencing (NGS) and three data analysis methods, this study’s objectives were to verify clonality of the Bavarian feed strains, rule out their connection to the outbreak, explore the genetic diversity of Bavarian S. Agona isolates from 1993 to 2018 and compare the analysis approaches employed, for practicality and ability to delineate outbreaks caused by the genetically monomorphic Agona serovar.MethodsIn this observational retrospective study, three 2017 Bavarian feed isolates were compared to a French outbreak isolate and 48 S. Agona isolates from our strain collections. The later included human, food, feed, veterinary and environmental isolates, of which 28 were epidemiologically outbreak related. All isolates were subjected to NGS and analysed by: (i) a publicly available species-specific core genome multilocus sequence typing (cgMLST) scheme, (ii) single nucleotide polymorphism phylogeny and (iii) an in-house serovar-specific cgMLST scheme. Using additional international S. Agona outbreak NGS data, the cluster resolution capacity of the two cgMLST schemes was assessed.ResultsWe could prove clonality of the feed isolates and exclude their relation to the French outbreak. All approaches confirmed former Bavarian epidemiological clusters.ConclusionEven for S. Agona, species-level cgMLST can produce reasonable resolution, being standardisable by public health laboratories. For single samples or homogeneous sample sets, higher resolution by serovar-specific cgMLST or SNP genotyping can facilitate outbreak investigations.

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

BackgroundLeptospirosis is an emerging infectious disease worldwide that can cause high morbidity and mortality rates in humans and animals. The causative spirochetes have reservoirs in mammalian hosts, but there has been limited analysis of the genomes of isolates recovered from animals. The aims of this study were to characterize genomic features of two Leptospira interrogans strains recently isolated from asymptomatic dogs in Thailand (strains CUDO5 and CDUO8), and to perform comparative genome analyses with other strains. Molecular adaptive evolution in L. interrogans as signaled by positive selection also was analyzed.ResultsWhole genome sequence analysis revealed that strains CUDO5 and CUDO8 had genome sizes of approximately 4.9 Mbp with 35.1% GC contents. Using monoclonal antibodies, strains CUDO5 and CUDO8 were identified as serovars Paidjan and Dadas, respectively. These strains harbored genes known to be associated with acute and chronic infections. Using Single Nucleotide Polymorphisms phylogeny (SNPs) with 97 L. interrogans strains, CUDO5 and CUDO8 had closest genetic relatedness with each other. Nevertheless, the serovar determinant region (rfb locus) showed variations in the genes encoding sugar biosynthesis. Amongst 13 representative L. interrogans strains examined for molecular adaptive evolution through positive selection under the site-model of Phylogenetic Analysis of Maximum Likelihood, genes responsible for iron acquisition (tlyA and hbpA), motility (fliN2, flgK, and flhB) and thermal adaptation (lpxD1) were under increased selective pressure.ConclusionsL. interrogans serovar Paidjan strain CUDO5 and serovar Dadas strain CUDO8 had close genetic relatedness as analyzed by SNPs phylogeny. They contained genes with established roles in acute and chronic leptospirosis. The rfb locus in both serovars showed gene variation associated with sugar biosynthesis. Positive selection analysis indicated that genes encoding factors involved in motility, temperature adaptation, and iron acquisition were under strong positive selection in L. interrogans. These may be associated with adaptation in the early stages of infection.

Redefining the differences in gene content between Yersinia pestis and Yersinia pseudotuberculosis using large-scale comparative genomics.

Yersinia pestis, the causative agent of plague, is best known for historical pandemics, but still actively causes disease in many parts of the world. Y. pestis is a recently derived clone of the pathogenic species Yersinia pseudotuberculosis, but is more associated with human infection. Numerous studies have documented genomic changes since the two species differentiated, although all of these studies used a relatively small sample set for defining these differences. In this study, we compared the complete genomic content between a diverse set of Y. pestis and Y. pseudotuberculosis genomes, and identified unique loci that could serve as diagnostic markers or for better understanding the evolution and pathogenesis of each group. Comparative genomics analyses also identified subtle variations in gene content between individual monophyletic clades within these species, based on a core genome single nucleotide polymorphism phylogeny that would have been undetected in a less comprehensive genome dataset. We also screened loci that were identified in other published studies as unique to either species and generally found a non-uniform distribution, suggesting that the assignment of these unique genes to either species should be re-evaluated in the context of current sequencing efforts. Overall, this study provides a high-resolution view into the genomic differences between Y. pestis and Y. pseudotuberculosis, demonstrating fine-scale differentiation and unique gene composition in both species.

Full genome SNP-based phylogenetic analysis reveals the origin and global spread of Brucella melitensis

BackgroundBrucellosis is an important zoonotic disease that affects both humans and animals. We sequenced the full genome and characterised the genetic diversity of two Brucella melitensis isolates from Malaysia and the Philippines. In addition, we performed a comparative whole-genome single nucleotide polymorphism (SNP) analysis of B. melitensis strains collected from around the world, to investigate the potential origin and the history of the global spread of B. melitensis.ResultsSingle sequencing runs of each genome resulted in draft genome sequences of MY1483/09 and Phil1136/12, which covered 99.85% and 99.92% of the complete genome sequences, respectively. The B. melitensis genome sequences, and two B. abortus strains used as the outgroup strains, yielded a total of 13,728 SNP sites. Phylogenetic analysis using whole-genome SNPs and geographical distribution of the isolates revealed spatial clustering of the B. melitensis isolates into five genotypes, I, II, III, IV and V. The Mediterranean strains, identified as genotype I, occupied the basal node of the phylogenetic tree, suggesting that B. melitensis may have originated from the Mediterranean regions. All of the Asian B. melitensis strains clustered into genotype II with the SEA strains, including the two isolates sequenced in this study, forming a distinct clade denoted here as genotype IId. Genotypes III, IV and V of B. melitensis demonstrated a restricted geographical distribution, with genotype III representing the African lineage, genotype IV representing the European lineage and genotype V representing the American lineage.ConclusionWe showed that SNPs retrieved from the B. melitensis draft full genomes were sufficient to resolve the interspecies relationships between B. melitensis strains and to discriminate between the vaccine and endemic strains. Phylogeographic reconstruction of the history of B. melitensis global spread at a finer scale by using whole-genome SNP analyses supported the origin of all B. melitensis strains from the Mediterranean region. The possible global distribution of B. melitensis following the ancient trade routes was also consistent with whole-genome SNP phylogeny. The whole genome SNP phylogenetics analysis, hence is a powerful tool for intraspecies discrimination of closely related species.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1294-x) contains supplementary material, which is available to authorized users.

Tracing melioidosis back to the source: using whole-genome sequencing to investigate an outbreak originating from a contaminated domestic water supply.

Melioidosis, a disease of public health importance in Southeast Asia and northern Australia, is caused by the Gram-negative soil bacillus Burkholderia pseudomallei. Melioidosis is typically acquired through environmental exposure, and case clusters are rare, even in regions where the disease is endemic. B. pseudomallei is classed as a tier 1 select agent by the Centers for Disease Control and Prevention; from a biodefense perspective, source attribution is vital in an outbreak scenario to rule out a deliberate release. Two cases of melioidosis within a 3-month period at a residence in rural northern Australia prompted an investigation to determine the source of exposure. B. pseudomallei isolates from the property's groundwater supply matched the multilocus sequence type of the clinical isolates. Whole-genome sequencing confirmed the water supply as the probable source of infection in both cases, with the clinical isolates differing from the likely infecting environmental strain by just one single nucleotide polymorphism (SNP) each. For the first time, we report a phylogenetic analysis of genomewide insertion/deletion (indel) data, an approach conventionally viewed as problematic due to high mutation rates and homoplasy. Our whole-genome indel analysis was concordant with the SNP phylogeny, and these two combined data sets provided greater resolution and a better fit with our epidemiological chronology of events. Collectively, this investigation represents a highly accurate account of source attribution in a melioidosis outbreak and gives further insight into a frequently overlooked reservoir of B. pseudomallei. Our methods and findings have important implications for outbreak source tracing of this bacterium and other highly recombinogenic pathogens.

Scalable SNP Analyses of 100+ Bacterial or Viral Genomes

With the flood of whole genome finished and draft microbial sequences, analysts need faster, more scalable bioinformatics tools for sequence comparison. An algorithm is described to find single nucleotide polymorphisms (SNPs) in whole genome data. It scales to hundreds of bacterial or viral genomes, and can be used for finished and/ or draft genomes available as unassembled contigs. The method is fast to compute, finding SNPs and building a SNP phylogeny in seconds to hours. It identified thousands of putative SNPs from all publicly available Filoviridae, Poxviridae, foot-and-mouth disease virus, Bacillus, and Escherichia coli genomes and plasmids. The SNP-based trees it generated were consistent with known taxonomy and trees determined in other studies. The approach described can handle as input hundreds of megabases of sequence in a single run. The algorithm kSNP is based on k-mer analysis using suffix arrays and requires no multiple sequence alignment. uncover novel regions that correlate with phenotype outside of wellcharacterized genes or non-coding sequence. It should also be useful in horizontal gene transfer studies, since one can examine SNPs across the entire genome. Although beyond the scope of this paper, microarrays with probes designed for all putative SNPs can be used to experimentally validate SNP alleles, identify sequencing errors, and characterize SNP alleles in unsequenced isolates to place them on a phylogeny (manuscript in preparation).