Search

Salmonella enterica serovar Enteritidis is one of the leading causes of salmonellosis in Australia. However, the majority of S. Enteritidis cases in Australia are travel-related with a small proportion being locally acquired. This study aimed to characterise the genomic features of Australian S. Enteritidis and compare them with international strains using multilevel genome typing (MGT). A total of 568 S. Enteritidis isolates from two Australian states across two consecutive years were analysed using the S. Enteritidis MGT scheme and database (MGTdb) - which contained 40,390 publicly available genomes from 99 countries. The Australian S. Enteritidis strains were divided into three phylogenetic clades (A, B and C). Clades A and C represented 16.4% and 3.5% of the total isolates, respectively, and were of local origin. Clade B accounted for 80.1% of the isolates which belonged to seven previously defined lineages but was dominated by the global epidemic lineage (MGT4-CC1). At MGT5 level, three out of five top sequence types (STs) in Australia were also top STs in Asia, suggesting that a fair proportion of Australian S. Enteritidis cases may be epidemiologically linked with Asian strains. In 2018, a large egg-associated local outbreak was caused by a recently defined clade B lineage prevalent in Europe and was closely related, but not directly linked, to three isolates from Europe. Additionally, antimicrobial-resistance genes were only found in Australian clade B isolates, with a predicted multidrug resistance (MDR) rate of 11.7%. Over half (54.8%) of the MDR isolates belonged to 10 MDR-associated MGT-STs, which were also frequent in Asian S. Enteritidis. IncX1 plasmids were frequently present in the Australian MDR isolates. Overall, this study investigated the genomic epidemiology of S. Enteritidis in Australia, including the first large local outbreak, using MGT. The open MGT platform enables a standardised and sharable nomenclature that can be effectively applied to public health for unified surveillance of S. Enteritidis nationally and globally.ImportanceSalmonella enterica serovar Enteritidis is a leading cause of foodborne infections. We previously developed a genomic typing database – MGTdb for S. Enteritidis to facilitate global surveillance of this pathogen. In this study we examined the genomic features of Australian S. Enteritidis using the MGTdb and found that Australian S. Enteritidis is mainly epidemiologically linked with Asian strains (especially strains carrying antimicrobial resistance genes) followed by European strains. The first large-scale egg-associated local outbreak in Australia was caused by a recently defined lineage prevalent in Europe, and three European isolates in the MGTdb were closely related but not directly linked to this outbreak. In summary, the S. Enteritidis MGTdb open platform is shown to be a potential powerful tool for national and global public health surveillance of this pathogen.

As a core microbe of the human gut ecosystem, Bacteroides vulgatus has been linked to multiple aspects of metabolic disorders in a collection of associative studies, which, while indicative, warrants more direct experimental evidence to verify. In this study, we experimentally demonstrated that oral administration of B. vulgatus Bv46 ameliorated the serum lipid profile and systemic inflammation of high-fat diet-induced hyperlipidemic rats in a microbiome-regulated manner, which appears to be associated with changes of bile acid metabolism, short-chain fatty acid biosynthesis, and serum metabolomic profile.

Listeria monocytogenes is a foodborne pathogen that can cause invasive disease with high mortality in immunocompromised individuals and can survive in a variety of food-associated environments for a long time. L. monocytogenes clonal complex (CC) 87 is composed of ST87 and three other STs and has been identified as the most common subgroup associated with both foods and human clinical infections in China. Therefore, the persistence of CC87 L. monocytogenes in food-associated environments poses a significant concern for food safety. In this study, 83 draft genomes of CC87 L. monocytogenes, including 60 newly sequenced genomes, were analyzed with all isolates from our previous surveillance in Zigong, Sichuang, China. Sixty-eight of the studied isolates were isolated from one retail market (M1 market), while the others were from seven other markets (M2-M8 markets) in the same city. Whole-genome multilocus sequence typing (wg-MLST) and the whole-genome single nucleotide polymorphism (wg-SNP) analysis were performed. Three persistent contamination routes were identified in the M1 market, caused by 2 clusters (A and B) and a wgST31 type. Cluster A isolates were associated with the persistent contamination in a raw meat stall (M1-S77), while Cluster B isolates caused a persistent contamination in aquatic foods stalls. Five wgST31 isolates caused persistent contamination in a single aquatic stall (M1-S65). A pLM1686-like plasmid was found in all Cluster A isolates. A novel plasmid, pLM1692, a truncated pLM1686 plasmid without the cadmium, and other heavy metal resistance genes were conserved in all wgST31 isolates. By comparing persistent and putative non-persistent isolates, four genes that were all located in the prophage comK might be associated with persistence. These findings enhanced our understanding of the underlying mechanisms of contamination and assist in formulating targeted strategies for the prevention and control of L. monocytogenes transmission from the food processing chain to humans.

Since acellular vaccines (ACV) were introduced in Australia, epidemic Bordetella pertussis strains changed from single nucleotide polymorphism (SNP) cluster II to SNP cluster I. Our previous proteomic analysis identified potential proteomic adaptations in the whole cell and secretome of SNP cluster I. Additionally, current ACVs were shown to be less efficacious against cluster I in mice models and there is a pressing need to discover new antigens to improve the ACV. One important source of novel antigens is the surfaceome. Therefore, in this study we established surface shaving in B. pertussis to compare the surfaceome of SNP cluster I (L1423) and II (L1191), and identify novel surface antigens for vaccine development. Surface shaving using 1 μg of trypsin for 5 min identified 126 proteins with the most abundant being virulence-associated and known outer membrane proteins. Cell viability counts showed minimal lysis from shaving. The proportion of immunogenic proteins was higher in the surfaceome than in the whole cell and secretome. Key differences in the surfaceome were identified between SNP cluster I and II, consistent with those identified in the whole cell proteome and secretome. These differences include unique transport proteins and decreased immunogenic proteins in L1423, and provides further evidence of proteomic adaptation in SNP cluster I. Finally, a comparison of proteins in each sub-proteome identified 22 common proteins. These included 11 virulence proteins (Prn, PtxA, FhaB, CyaA, TcfA, SphB1, Vag8, BrkA, BopD, Bsp22 and BipA) and 11 housekeeping proteins (TuF, CtpA, TsF, OmpH, GltA, SucC, SucD, FusA, GroEL, BP3330 and BP3561) which were immunogenic, essential and consistently expressed thus demonstrating their potential as future targets. This study established surface shaving in B. pertussis, confirmed key expression differences and identified unknown surface proteins which may be potential vaccine antigens.

Background Listeria monocytogenes consists of four lineages that occupy a wide variety of ecological niches. Sequence type (ST) 87 (serotype 1/2b), belonging to lineage I, is one of the most common STs isolated from food products, food associated environments and sporadic listeriosis in China. Here, we performed a comparative genomic analysis of the L. monocytogenes ST87 clone by sequencing 71 strains representing a diverse range of sources, different geographical locations and isolation years. Results The core genome and pan genome of ST87 contained 2667 genes and 3687 genes respectively. Phylogenetic analysis based on core genome SNPs divided the 71 strains into 10 clades. The clinical strains were distributed among multiple clades. Four clades contained strains from multiple geographic regions and showed high genetic diversity. The major gene content variation of ST87 genomes was due to putative prophages, with eleven hotspots of the genome that harbor prophages. All strains carry an intact CRISRP/Cas system. Two major CRISPR spacer profiles were found which were not clustered phylogenetically. A large plasmid of about 90 Kb, which carried heavy metal resistance genes, was found in 32.4% (23/71) of the strains. All ST87 strains harbored the Listeria pathogenicity island (LIPI)-4 and a unique 10-open read frame (ORF) genomic island containing a novel restriction-modification system. Conclusion Whole genome sequence analysis of L. monocytogenes ST87 enabled a clearer understanding of the population structure and the evolutionary history of ST87 L. monocytogenes in China. The novel genetic elements identified may contribute to its virulence and adaptation to different environmental niches. Our findings will be useful for the development of effective strategies for the prevention and treatment of listeriosis caused by this prevalent clone.

Shigella and enteroinvasive Escherichia coli (EIEC) cause human bacillary dysentery with similar invasion mechanisms and share similar physiological, biochemical and genetic characteristics. Differentiation of Shigella from EIEC is important for clinical diagnostic and epidemiological investigations. However, phylogenetically, Shigella and EIEC strains are composed of multiple clusters and are different forms of E. coli , making it difficult to find genetic markers to discriminate between Shigella and EIEC. In this study, we identified 10 Shigella clusters, seven EIEC clusters and 53 sporadic types of EIEC by examining over 17000 publicly available Shigella and EIEC genomes. We compared Shigella and EIEC accessory genomes to identify cluster-specific gene markers for the 17 clusters and 53 sporadic types. The cluster-specific gene markers showed 99.64% accuracy and more than 97.02% specificity. In addition, we developed a freely available in silico serotyping pipeline named Shigella EIEC Cluster Enhanced Serotype Finder (ShigEiFinder) by incorporating the cluster-specific gene markers and established Shigella and EIEC serotype-specific O antigen genes and modification genes into typing. ShigEiFinder can process either paired-end Illumina sequencing reads or assembled genomes and almost perfectly differentiated Shigella from EIEC with 99.70 and 99.74% cluster assignment accuracy for the assembled genomes and read mapping respectively. ShigEiFinder was able to serotype over 59 Shigella serotypes and 22 EIEC serotypes and provided a high specificity of 99.40% for assembled genomes and 99.38% for read mapping for serotyping. The cluster-specific gene markers and our new serotyping tool, ShigEiFinder (installable package: https://github.com/LanLab/ShigEiFinder, online tool: https://mgtdb.unsw.edu.au/ShigEiFinder/), will be useful for epidemiological and diagnostic investigations.

Since 1817, cholera, caused by Vibrio cholerae, has been characterized by seven distinct pandemics. The ongoing seventh pandemic (7P) began in 1961. In this study, we developed a Multilevel Genome Typing (MGT) tool for classifying the V. cholerae species with a focus on the 7P. MGT is based on multilocus sequence typing (MLST), but the concept has been expanded to include a series of MLST schemes that compare population structure from broad to fine resolutions. The V. cholerae MGT consists of eight levels, with the lowest, MGT1, composed of 7 loci and the highest, MGT8, consisting of the 7P core genome (3,759 loci). We used MGT to analyze 5,771 V. cholerae genomes. The genetic relationships revealed by lower MGT levels recapitulated previous findings of large-scale 7P transmission across the globe. Furthermore, the higher MGT levels provided an increased discriminatory power to differentiate subgroups within a national outbreak. Additionally, we demonstrated the usefulness of MGT for non-7P classification. In a large non-7P MGT1 type, MGT2 and MGT3 described continental and regional distributions, respectively. Finally, MGT described trends of 7P in virulence, and MGT2 to MGT3 sequence types (STs) grouped isolates of the same ctxB, tcpA, and ctxB-tcpA genotypes and characterized their trends over the pandemic. MGT offers a range of resolutions for typing V. cholerae. The MGT nomenclature is stable, transferable, and directly comparable between investigations. The MGT database (https://mgtdb.unsw.edu.au/) can accept and process newly submitted samples. MGT allows tracking of existing and new isolates and will be useful for understanding future spread of cholera. IMPORTANCE In 2017, the World Health Organization launched the “Ending Cholera” initiative to reduce cholera-related deaths by 90% by 2030. This strategy emphasized the importance of the speed and accessibility of newer technologies to contain outbreaks. Here, we present a new tool named Multilevel Genome Typing (MGT), which classifies isolates of the cholera-causing agent, Vibrio cholerae. MGT is a freely available online database that groups genetically similar V. cholerae isolates to quickly indicate the origins of outbreaks. We validated the MGT database retrospectively in an outbreak setting, showcasing rapid confirmation of the Nepalese origins for the 2010 Haiti outbreak. In the past 5 years, thousands of V. cholerae genomes have been submitted to the NCBI database, which underscores the importance of and need for proper genome data classification for cholera epidemiology. The V. cholerae MGT database can assist in early decision making that directly impacts controlling both the local and global spread of cholera.

Salmonella enterica serovar Enteritidis is a major cause of foodborne Salmonella infections and outbreaks in humans. Effective surveillance and timely outbreak detection are essential for public health control. Multilevel genome typing (MGT) with multiple levels of resolution has been previously demonstrated as a promising tool for this purpose. In this study, we developed MGT with nine levels for S. Enteritidis and characterised the genomic epidemiology of S. Enteritidis in detail. We examined 26 670 publicly available S. Enteritidis genome sequences from isolates spanning 101 years from 86 countries to reveal their spatial and temporal distributions. Using the lower resolution MGT levels, globally prevalent and regionally restricted sequence types (STs) were identified; avian associated MGT4-STs were found that were common in human cases in the USA; temporal trends were observed in the UK with MGT5-STs from 2014 to 2018 revealing both long lived endemic STs and the rapid expansion of new STs. Using MGT3 to MGT6, we identified multidrug resistance (MDR) associated STs at various MGT levels, which improves precision of detection and global tracking of MDR clones. We also found that the majority of the global S. Enteritidis population fell within two predominant lineages, which had significantly different propensity of causing large scale outbreaks. An online open MGT database has been established for unified international surveillance of S. Enteritidis. We demonstrated that MGT provides a flexible and high-resolution genome typing tool for S. Enteritidis surveillance and outbreak detection.

Purpose. Human-adapted Bordetella parapertussis is one of the causative agents of whooping cough; however, there are currently no genotyping systems with high discriminatory power for this bacterial pathogen. We therefore aimed to develop a multilocus variable-number tandem repeat analysis (MLVA) for human-adapted B. parapertussis.Methodology. Four highly polymorphic variable number tandem repeat (VNTR) loci in the B. parapertussis genome were selected and amplified by multiplex PCR. MLVA was performed based on the number of tandem repeats at VNTR loci. The discriminatory power of MLVA was evaluated with three laboratory reference strains and 50 human isolates of B. parapertussis.Results. Multiplex PCR-based MLVA characterized 53 B. parapertussis reference strains and isolates into 25 MLVA types and the Simpson diversity index was 0.91 (95 % confidence interval, 0.86-0.97). The three reference strains exhibited different MLVA types. Thirty-one Japanese isolates, ten French isolates and three Taiwanese isolates belonged to fourteen, nine and three MLVA types, respectively. In contrast, all five Australian isolates belonged to the same type. Two Japanese isolates collected from patients with known epidemiological links had the same type.Conclusion. Our novel MLVA method has high discriminatory power for genotyping human B. parapertussis. Regarding this organism, this genotyping system is a promising tool for epidemiological surveillance and investigating outbreaks.

Macrolides such as erythromycin are the empirical treatment of Bordetella pertussis infections. China has experienced an increase in erythromycin-resistant B. pertussis isolates since they were first reported in 2013. Here, we undertook a genomic study on Chinese B. pertussis isolates from 2012 to 2015 to elucidate the origins and phylogenetic relationships of erythromycin-resistant B. pertussis isolates in China. A total of 167 Chinese B. pertussis isolates were used for antibiotic sensitivity testing and multiple locus variable-number tandem repeat (VNTR) analysis (MLVA). All except four isolates were erythromycin-resistant and of the four erythromycin-sensitive isolates, three were non-ptxP1. MLVA types (MT), MT55, MT104 and MT195 were the predominant types. Fifty of those isolates were used for whole genome sequencing and phylogenetic analysis. Genome sequencing and phylogenetic analysis revealed three independent erythromycin-resistant lineages and all resistant isolates carried a mutation in the 23S rRNA gene. A novel fhaB3 allele was found uniquely in Chinese ptxP1 isolates and these Chinese ptxP1-ptxA1-fhaB3 had a 5-fold higher mutation rate than the global ptxP1-ptxA1 B. pertussis population. Our results suggest that the evolution of Chinese B. pertussis is likely to be driven by selection pressure from both vaccination and antibiotics. The emergence of the new non-vaccine fhaB3 allele in Chinese B. pertussis population may be a result of selection from vaccination, whereas the expansion of ptxP1-fhaB3 lineages was most likely to be the result of selection pressure from antibiotics. Further monitoring of B. pertussis in China is required to better understand the evolution of the pathogen.

Escherichia albertii is a newly recognized species in the genus Escherichia that causes diarrhea. The population structure, genetic diversity and genomic features has not been fully examined. Here, 169 E. albertii isolates from different sources and regions in China were sequenced and combined with 312 publicly available genomes for phylogenetic and genomic analyses. The E. albertii population was divided into 2 clades and 8 lineages, with lineage 3 (L3), L5 and L8 more common in China. Clinical isolates were observed in all clades/lineages. Virulence genes were found to be distributed differently among lineages: subtypes of the intimin encoding gene eae and the cytolethal distending toxin (Cdt) gene cdtB were lineage associated, the second type three secretion system (ETT2) island was truncated in L3 and L6. Seven new eae subtypes and 1 new cdtB subtype (cdtB-VI) were found. Alarmingly, 85.9% of the Chinese E. albertii isolates were predicted to be multidrug resistant (MDR) with 35.9% harboured genes capable of conferring resistance to 10 to 14 different drug classes. By in silico multi-locus sequence typing, majority of the MDR isolates belonged to 4 STs (ST4638, ST4479, ST4633 and ST4488). Thirty-four intact plasmids carrying MDR and virulence genes, and 130 intact prophages were identified from 17 complete E. albertii genomes. Ten plasmid replicon types were found to be significantly associated with MDR. The 130 intact prophages were clustered into 5 groups, with group 5 prophages harbouring more virulence genes. Our findings provided fundamental insights into the population structure, virulence variation and MDR of E. albertii.Impact statementE. albertii is newly recognized foodborne pathogen causing diarrhea. Elucidation of its genomic features is important for surveillance and control of E. albertii infections. In this work, 169 E. albertii genomes from difference sources and regions in China were collected and sequenced, which contributed to the currently limited genomic data pool of E. albertii. In combination with 312 publicly available genomes from 14 additional countries, the population structure of E. albertii was defined. The presence and subtypes of virulence genes in different lineages were significantly different, indicating potential pathogenicity variation. Additionally, the presence of multidrug resistance (MDR) genes was alarmingly high in the Chinese dominated lineages. MDR associated STs and plasmid subtypes were identified, which could be used as sentinels for MDR surveillance. Moreover, the subtypes of plasmids and prophages were distributed differently across lineages, and were found to contribute to the acquisition of virulence and MDR genes of E. albertii. Altogether, this work reveals the diversity of E. albertii and characterized its genomic features in unprecedented detail.Data SummaryAll newly sequenced data in this work were deposited in National Center for Biotechnology Information (NCBI) under the BioProject of PRJNA693666, including 6 complete genomes and raw reads of 164 E. albertii isolates.

The Bordetella genus is divided into two groups: classical and non-classical. Bordetella pertussis, Bordetella bronchiseptica and Bordetella parapertussis are known as classical bordetellae, a group of important human pathogens causing whooping cough or whooping cough-like disease and hypothesized to have evolved from environmental non-classical bordetellae. Bordetella infections have increased globally driving the need to better understand these pathogens for the development of new treatments and vaccines. One unexplored component in Bordetella is the role of serine, threonine and tyrosine phosphorylation. Therefore, this study characterized the phosphoproteome of classical bordetellae and examined its potential role in Bordetella biology and virulence. Applying strict identification of localization criteria, this study identified 70 unique phosphorylated proteins in the classical bordetellae group with a high degree of conservation. Phosphorylation was a key regulator of Bordetella metabolism with proteins involved in gluconeogenesis, TCA cycle, amino acid and nucleotide synthesis significantly enriched. Three key virulence pathways were also phosphorylated including type III secretion system, alcaligin synthesis and the BvgAS master transcriptional regulatory system for virulence genes in Bordetella. Seven new phosphosites were identified in BvgA with 6 located in the DNA binding domain. Of the 7, 4 were not present in non-classical bordetellae. This suggests that serine/threonine phosphorylation may play an important role in stabilizing/destabilizing BvgA binding to DNA for fine-tuning of virulence gene expression and that BvgA phosphorylation may be an important factor separating classical from non-classical bordetellae. This study provides the first insight into the phosphoproteome of classical Bordetella species and the role that Ser/Thr/Tyr phosphorylation may play in Bordetella biology and virulence.

Most human infectious diseases originate from animals. Thus, how to reduce or prevent pandemic zoonoses before they emerge in people is becoming a critical issue., Our ability to predict evolutionary trajectories of pathogens is one of the promising leverages to fight against the pandemic disease, yet few studies have addressed this question in situ, due to the difficulty in monitoring the milestone evolutionary events for a given pathogen and in understanding the evolutionary strategies. In this study, we monitored the real-time evolution of Vibrio parahaemolyticus in response to successive antibiotic treatment in three shrimp farms in North China from 2011 to 2018 by whole-genome sequencing. Results showed that the stepwise emergence of resistance was associated with the antibiotic usage. Genomic analysis of resistant isolates showed that the acquisition of the resistant mobile genetic elements flanked by an insertion sequence (ISVal1) closely mirrored the antibiotics used in shrimp farms since 2014. Next, we also identified 50 insertion sites of ISVal1 in the chromosome, which facilitated the formation of pathogenicity islands (PAIs) and fitness islands in the following years. Further, horizontal transfers of a virulent trh-nik-ure genomic island (GI) and two GIs improving the fitness have been observed in two farms since 2016. In this case study, we proposed that the insertion sequence triggered four major evolutionary events during the outbreaks of shrimp disease in three farms, including horizontal transfer of transposon (HTT) (stage 1), the formation of resistance islands (stage 2) and the PAIs (stage 3), and horizontal transfer of the PAIs (stage 4). This study presented the first in vivo evolutionary trajectories for a given bacterial pathogen, which helps us to understand the emergence mechanisms of new genotypes. IMPORTANCE Most human infectious diseases originate from animals. Thus, how to reduce or prevent pandemic zoonoses before they emerge in people is becoming a critical issue. Continuous genomic surveillance of the evolutionary trajectories of potential human pathogens on farms is a promising strategy to realize early warning. Here, we conducted an 8-year surveillance of Vibrio parahaemolyticus in three shrimp farms. The results showed that the use of antibiotics and horizontal transfer of transposons (HTT) drove the evolution of V. parahaemolyticus, which could be divided into four stages: HTT, formation of resistance islands, formation of pathogenicity islands (PAIs), and horizontal transfer of PAIs. This study presented the first in vivo monitoring of evolutionary trajectories for a given bacterial pathogen, providing valuable information for the prevention of pandemic zoonoses.

Here we investigated nationwide clinical Bordetella pertussis isolated during 2005–2017 from different provinces of Iran, a country with more than 50 years whole cell vaccine immunisation history. Our results revealed the ongoing increase in the population of ptxP3/fim3–2 B. pertussis isolates in different provinces which were differentiated into nine clades. The largest clade (clade 8) which was previously found to be prevalent in Tehran was also prevalent across the country and clade 5 with ptxP3/prn9 genotype has also increased in frequency (14% of all ptxP3 isolates) in recent years. Furthermore, we detected the first ptxP3 B. pertussis isolates that does not express filamentous hemagglutinin (FhaB) as one of the major antigens of the pathogen and a key component of the acellular pertussis vaccine.

Campylobacter concisusis an emerging enteric pathogen that is associated with several gastrointestinal diseases, such as inflammatory bowel disease (IBD), which includes Crohn’s disease (CD) and ulcerative colitis (UC). Currently, only three completeC. concisusgenomes are available and more completeC. concisusgenomes are needed in order to better understand the genomic features and pathogenicity of this emerging pathogen. DNA extracted from 22C.concisusstrains were subjected to Oxford Nanopore genome sequencing. Complete genome assembly was performed using Nanopore genome data in combination with previously reported short-read Illumina data. Genome features of completeC. concisusgenomes were analysed using bioinformatic tools. The enteric disease associations ofC. concisusplasmids were examined using 239C.concisusstrains and confirmed using PCRs. Proteomic analysis was used to examine T6SS secreted proteins. We successfully obtained 13 completeC. concisusgenomes in this study. Analysis of 16 completeC. concisusgenomes (3 from public databases) identified multiple novel plasmids. pSma1 plasmid was found to be associated with severe UC. Sec-SRP, Tat and T6SS were found to be the main secretion systems inC. concisusand proteomic data showed a functional T6SS despite the lack of ClpV. T4SS was found in 25% of completeC. concisusgenomes. This study also found that GS2 strains had larger genomes and higher GC content than GS1 strains and more often had plasmids. In conclusion, this study provides fundamental genomic data for understandingC. concisusplasmids, genomospecies features, evolution, secretion systems and pathogenicity.

Summary Background Corynebacteritum straitum has been considered as an emerging multi-drug-resistant (MDR) pathogen. Isolation of MDR C. striatum as the only organism from respiratory samples from hospitalized patients is increasing in China. Aim To elucidate the genomic epidemiology and evolution of C. striatum in China. Methods A total of 260 isolates from 2016 to 2018 were collected from three hospitals in three regions of China. Antibiotic sensitivity testing was performed on all isolates. Whole-genome sequencing was applied to all isolates to assess their genomic diversity and relationships and detect the presence of antimicrobial resistance genes (ARG) and ARG cassettes. Findings Almost all isolates (96.2%, 250/260) showed multi-drug-resistance. Genome sequencing revealed four major lineages with lineage IV emerging as the epidemic lineage. Most of the diversity was developed in the last 6 years. Each hospital has its own predominant clones with potential spread between Hebei and Guangdong hospitals. Genomic analysis further revealed multiple antimicrobial resistance genes. Conclusions Our results suggested that four lineages of C. striatum have spread in parallel across China, causing persistent and extensive transmissions within hospitals. MDR C. striatum infection has become a national epidemic. Antibiotic-driven selection pressure may have played significant roles in forming persistent and predominant clones. Our data provide the basis for surveillance and prevention strategies to control the epidemic caused by MDR C. striatum.

Global outbreaks of shrimp acute hepatopancreatic necrosis disease (AHPND) caused by V. parahaemolyticus represent an urgent issue for the shrimp industry. This study revealed that the transmission mode of AHPND consists of two steps, the transregional dissemination of V. parahaemolyticus and the horizontal transfer of an AHPND-associated plasmid. Surprisingly, the introduction of the AHPND-associated plasmid also offers a novel mechanism of genetic exchange mediated by insertion sequences, and it improved the fitness of V. parahaemolyticus in a harsh environment. The results presented herein suggest that current shrimp farming practices promote genetic mixture between endemic and oceanic V. parahaemolyticus populations, which introduced the plasmid and accelerated bacterial adaptation by the acquisition of ecologically important functions. This entails a risk of the emergence of new virulent populations both for shrimp and humans. This study improves our understanding of the global dissemination of the AHPND-associated plasmid and highlights the urgent need to improve biosecurity for shrimp farming., Vibrio parahaemolyticus is an important foodborne pathogen and has recently gained particular notoriety because it causes acute hepatopancreatic necrosis disease (AHPND) in shrimp, which has caused significant economic loss in the shrimp industry. Here, we report a whole-genome analysis of 233 V. parahaemolyticus strains isolated from humans, diseased shrimp, and environmental samples collected between 2008 and 2017, providing unprecedented insight into the historical spread of AHPND. The results show that V. parahaemolyticus is genetically diverse and can be divided into 84 sequence types (STs). However, genomic analysis of three STs of V. parahaemolyticus identified seven transmission routes in Asia since 1996, which promoted the transfer of an AHPND-associated plasmid. Notably, the insertion sequence (ISVal1) from the plasmid subsequently mediated the genetic exchange among V. parahaemolyticus STs and resulted in the deletion of an 11-kb region regulating cell mobility and the production of capsular polysaccharides. Phenotype assays confirmed that this deletion enhanced biofilm formation, providing a novel mechanism for environmental adaptation. We conclude that the transmission mode of AHPND consists of two steps, the transmission of V. parahaemolyticus and the subsequent horizontal transfer of the AHPND-associated plasmid. This plasmid allows ISVal1 to mediate genetic exchange and improve pathogen fitness in shrimp ponds. Current shrimp farming practices promoted such genetic exchanges, which highlighted a risk of the emergence of new virulent populations, with potentially devastating consequences for both aquaculture and human health. This study addressed the basic questions regarding the transmission mechanism of AHPND and provided novel insights into shrimp and human disease management. IMPORTANCE Global outbreaks of shrimp acute hepatopancreatic necrosis disease (AHPND) caused by V. parahaemolyticus represent an urgent issue for the shrimp industry. This study revealed that the transmission mode of AHPND consists of two steps, the transregional dissemination of V. parahaemolyticus and the horizontal transfer of an AHPND-associated plasmid. Surprisingly, the introduction of the AHPND-associated plasmid also offers a novel mechanism of genetic exchange mediated by insertion sequences, and it improved the fitness of V. parahaemolyticus in a harsh environment. The results presented herein suggest that current shrimp farming practices promote genetic mixture between endemic and oceanic V. parahaemolyticus populations, which introduced the plasmid and accelerated bacterial adaptation by the acquisition of ecologically important functions. This entails a risk of the emergence of new virulent populations both for shrimp and humans. This study improves our understanding of the global dissemination of the AHPND-associated plasmid and highlights the urgent need to improve biosecurity for shrimp farming.

We sequenced the genomes of 14 sequential Salmonella enterica serovar Typhimurium isolates obtained over a five year period from a patient with persistent Salmonella bacteriuria. The isolates formed five distinct lineages; two of which co-existed over four years. We inferred that the observed within-patient variation resulted from mutation events.

Shiga toxin producing Escherichia coli O157:H7 (STEC O157) is naturally found in the gastrointestinal tract of cattle and can cause severe disease in humans. There is limited understanding of the population dynamics and microevolution of STEC O157 at herd level. In this study, isolates from a closed beef herd of 23 cows were used to examine the population turnover in the herd. Of the nine STEC O157 clades previously described, clade 7 was found in 162 of the 169 isolates typed. Multiple locus variable number tandem repeat analysis (MLVA) differentiated 169 isolates into 33 unique MLVA types. Five predominant MLVA types were evident with most of the remaining types containing only a single isolate. MLVA data suggest that over time clonal replacement occurred within the herd. Genome sequencing of 18 selected isolates found that the isolates were divided into four lineages, representing four different 'clones' in the herd. Genome data confirmed clonal replacement over time and provided evidence of cross transmission of strains between cows. The findings enhanced our understanding of the population dynamics of STEC O157 in its natural host that will help developing effective control measures to prevent the spread of the pathogen to the human population.

Streptococcus suis is an important pathogen of pigs and may cause serious disease in humans. Serotyping is an important tool for detection and epidemiological studies of S. suis . Thirty-three reference serotypes and nine novel cps loci (NCLs) are recognized in S. suis . To gain a better understanding of the prevalence and genetic characteristics of NCLs, we investigated the serotype identity of 486 isolates isolated between 2013 and 2015 in China by capsular gene typing methods. Two hundred seventy-six isolates carried NCLs belonging to 16 groups, 8 of which appear to have not been reported previously. These isolates showed autoagglutination, polyagglutination, or nonagglutination with reference antisera and thus were nonserotypeable. Almost all isolates carrying the unknown NCLs were encapsulated, with various capsular thicknesses, indicating that they are most likely novel serotypes. To simultaneously identify the currently recognized 17 NCLs, an 18-plex detection system using the Luminex xTAG universal array technology was developed. Our data also provide valuable genetic information for monitoring the variations within NCLs by investigating the genetic characteristics of different subtypes within NCLs. IMPORTANCE Nonserotypeable Streptococcus suis isolates have been reported in many studies, and 9 novel cps loci (NCLs) have already been identified in nonserotypeable isolates. Moreover, novel cps loci are continually being found. The main purpose of this study was to investigate the prevalence and characteristics of NCLs in S. suis isolates recovered between 2013 and 2015 in China. This study provides valuable genetic information for monitoring the variations within NCLs. Meanwhile, a fast and cost-effective 18-plex detection system that can simultaneously identify the currently recognized 17 NCLs was developed in this study. This system will serve as a valuable tool for detecting known and identifying additional novel cps loci among nonserotypeable S. suis isolates.

De novo assembly of bacterial genomes from next-generation sequencing (NGS) data allows a reference-free discovery of single nucleotide polymorphisms (SNP). However, substantial rates of errors in genomes assembled by this approach remain a major barrier for the reference-free analysis of genome variations in medically important bacteria. The aim of this report was to improve the quality of SNP identification in bacterial genomes without closely related references. We developed a bioinformatics pipeline (SnpFilt) that constructs an assembly using SPAdes and then removes unreliable regions based on the quality and coverage of re-aligned reads at neighbouring regions. The performance of the pipeline was compared against reference-based SNP calling for Illumina HiSeq, MiSeq and NextSeq reads from a range of bacterial pathogens including Salmonella, which is one of the most common causes of food-borne disease. The SnpFilt pipeline removed all false SNP in all test NGS datasets consisting of paired-end Illumina reads. We also showed that for reliable and complete SNP calls, at least 40-fold coverage is required. Analysis of bacterial isolates associated with epidemiologically confirmed outbreaks using the SnpFilt pipeline produced results consistent with previously published findings. The SnpFilt pipeline improves the quality of de-novo assembly and precision of SNP calling in bacterial genomes by removal of regions of the assembly that may potentially contain assembly errors. SnpFilt is available from https://github.com/LanLab/SnpFilt.

Streptococcus suis is an important zoonotic pathogen. Serotype 2 and sequence type (ST) 1 are the most frequently reported strains in both infected humans and pigs. ST7 is only endemic to China, and it was responsible for outbreaks in 1998 and 2005 in China. In the present study, 38 sporadic ST7 S. suis strains, which mostly caused sepsis, were collected from patients in the Guangxi Zhuang Autonomous Region (GX) between 2007 and 2018. Of 38 sporadic ST7 strains, serotype 14 was the most frequent (27 strains, 71.1%), followed by serotype 2 (11 strains, 28.9%). The phylogenetic structure of the ST7 population, including epidemic and sporadic ST7 strains, was constructed using mutational single-nucleotide polymorphisms (SNPs). High diversity within the ST7 population was revealed and divided into five lineages. Only one sporadic ST7 strain, GX14, from a Streptococcal toxic-shock-like syndrome (STSLS) patient was clustered into the same lineage as the epidemic strains. GX14 and the epidemic strains diverged in 1974. The sporadic ST7 strains of GX were mainly clustered into lineage 5, which emerged in 1980. Comparing to genome of epidemic strain, the major differences in genome of sporadic ST7 strains of GX was the absence of 89 kb pathogenicity island (PAI) specific to epidemic strain and insertion of 128 kb ICE_phage tandem MGE or ICE portion of the MGE. These mobile elements play a significant role in the horizontal transfer of antibiotic resistance genes in sporadic ST7 strains. Our results enhanced the understanding of the evolution of the ST7 strains and their ability to cause life-threatening infections in humans.

Salmonella enterica serovar Typhimurium is the leading cause of salmonellosis in Australia, and the ability to identify outbreaks and their sources is vital to public health. Here, we examined the utility of whole-genome sequencing (WGS), including complete genome sequencing with Oxford Nanopore technologies, in examining 105 isolates from an endemic multi-locus variable number tandem repeat analysis (MLVA) type over 5 years. The MLVA type was very homogeneous, with 90 % of the isolates falling into groups with a five SNP cut-off. We developed a new two-step approach for outbreak detection using WGS. The first clustering at a zero single nucleotide polymorphism (SNP) cut-off was used to detect outbreak clusters that each occurred within a 4 week window and then a second clustering with dynamically increased SNP cut-offs were used to generate outbreak investigation clusters capable of identifying all outbreak cases. This approach offered optimal specificity and sensitivity for outbreak detection and investigation, in particular of those caused by endemic MLVA types or clones with low genetic diversity. We further showed that inclusion of complete genome sequences detected no additional mutational events for genomic outbreak surveillance. Phylogenetic analysis found that the MLVA type was likely to have been derived recently from a single source that persisted over 5 years, and seeded numerous sporadic infections and outbreaks. Our findings suggest that SNP cut-offs for outbreak cluster detection and public-health surveillance should be based on the local diversity of the relevant strains over time. These findings have general applicability to outbreak detection of bacterial pathogens.

Salmonella enterica subspecies enterica is a highly diverse subspecies with more than 2600 serovars and the ability to distinguish serovars within this group is vital for surveillance. With the development of whole-genome sequencing technology, serovar prediction by traditional serotyping is being replaced by molecular serotyping. Existing in silico serovar prediction approaches utilise surface antigen encoding genes, core genome MLST and serovar-specific gene markers or DNA fragments for serotyping. However, these serovar-specific gene markers or DNA fragments only distinguished a smaller number of serovars. In this study, we compared 2258 Salmonella accessory genomes to identify 414 candidate serovar-specific or lineage specific gene markers for 106 serovars which includes 24 polyphyletic serovars and the paraphyletic serovar Enteritidis. A combination of several lineage-specific gene markers can be used for the clear identification of the polyphyletic serovars and paraphyletic serovars. We designed and evaluated an in silico serovar prediction approach by screening 1089 genomes representing 106 serovars against a set of 131 serovar-specific gene markers. The presence or absence of one or more serovar-specific gene markers was used to predict the serovar of an isolate from genomic data. We show that serovar-specific gene markers have comparable accuracy to other in silico serotyping methods with 84.8% of isolates assigned to the correct serovar with no false positives and false negatives and 10.5% of isolates assigned to a small subset of serovars containing the correct serovar with varied false positives. Combined, 95.3% of genomes were correctly assigned to a serovar. This approach would be useful as diagnosis moves to culture-independent and metagenomic methods as well as providing a third alternative to confirm other genome-based analyses. The identification of a set of gene markers may also be useful in the development of more cost-effective molecular assays designed to detect specific gene markers of the all major serovars in a region. These assays would be useful in serotyping isolates where cultures are no longer obtained and traditional serotyping is therefore impossible.

Pertussis remains one of the least controlled vaccine-preventable diseases despite high vaccine coverage in many countries. There are ongoing debates about the causes of its resurgence. Major changes have occurred in the Bordetella pertussis population since the introduction of vaccination. Currently circulating strains in Australia and many other developed countries are grouped in single nucleotide polymorphism (SNP) cluster I (also known as ptxP3 strains). The emergence and expansion of SNP cluster I has been associated with two major genetic changes in B. pertussis: a change in its pertussis toxin promoter (to ptxP3) which leads to increased pertussis toxin production and the change of the acellular vaccine pertactin gene allele to prn2. More recently, strains that lack pertactin have emerged independently in different lineages. The resurgence of pertussis in highly vaccinated populations can be, at least in part, explained by genetic changes that increase the fitness of circulating B. pertussis strains.

Five strains of a Gram-stain-positive, catalase-negative, α-haemolytic, coccus-shaped chain-forming organism were isolated separately from the lower respiratory tracts of five animals of Marmota himalayana in the endemic area of plague, the Qinghai-Tibet Plateau, China. Based on their morphological characteristics, biochemical features and molecular phylogenetic studies, the strains were placed as representing a new member of the genus Streptococcus . Comparative 16S rRNA gene sequence studies indicated that strain HTS5T shared 96.5, 96.2 and 96.0 % similarity with Streptococcus gallinaceus CCUG 42692T, Streptococcus parasanguinis ATCC 15912T and Streptococcus suis ATCC 43765T, respectively. Sequence analysis of its rpoB and sodA genes showed that strain HTS5T was most closely related to Streptococcus cuniculi CCUG 65085T with 9.2 and 10.9 % interspecies divergence, respectively. The whole genome phylogenetic tree based on 339 core genes of 65 Streptococcus genomes confirmed that HTS5T belongs to a distinct lineage that is well separated from recognized species of the genus Streptococcus . In silico DNA–DNA hybridization using 65 available genomes from GenBank showed that HTS5T displayed less than 70 % DNA–DNA relatedness with the other 65 species of the genus Streptococcus deposited in the GenBank database. The genome of strain HTS5T (2 322 791 bp) contained 2377 genes and had a G+C content of 41.6 mol%. Therefore, the five strains are considered to represent a novel species of the genus Streptococcus for which the name Streptococcus marmotae sp. nov. is proposed. The type strain is HTS5T (=DSM 101995T=CGMCC 1.15534T).

Purpose. Both Shigella and enteroinvasive Escherichia coli (EIEC) can cause enterocolitis, but they have a distinct epidemiology and public health relevance. Current culture-independent testing (CIT) methods to identify Shigella in faecal samples rely on the ipaH gene as the target, which is also found in EIEC genomes. The aim of this study was to design an assay that can identify EIEC in cultures from CIT ipaH-positive samples. Methodology. Shigella and EIEC genomes were screened to find unique regions present in EIEC genomes using a comparative genomics approach and differentiating genetic loci that are suitable PCR targets were identified. The primers for these loci were designed and tested in 6501 and 104 genomes of Shigella and EIEC, respectively. Results. An assay with two sets of multiplex PCR reactions that differentiates Shigella and EIEC based on the presence/absence of at least two out of six loci was developed and evaluated. The majority of Shigella genomes lacked all six loci, while at least two loci were present in most EIEC genomes. This assay successfully differentiated clinical EIEC from Shigella with a limit of detection of 105 cells ml−1. The sensitivity and specificity were over 95 and 99%, respectively. The assay can further subtype EIEC genomes into their genetic lineages. Conclusion. This new highly specific assay can assist in the identification of EIEC in ipaH PCR-positive samples and augment the public health laboratory surveillance of EIEC and shigellosis.

Salmonella enterica serovar Typhimurium is the most common Salmonella serovar causing foodborne infections in Australia and many other countries. Twenty-one S . Typhimurium strains from Salmonella reference collection A (SARA) were analyzed using Illumina high-throughput genome sequencing. Single nucleotide polymorphisms (SNPs) in 21 SARA strains ranged from 46 to 11,916 SNPs, with an average of 1,577 SNPs per strain. Together with 47 strains selected from publicly available S . Typhimurium genomes, the S . Typhimurium core genes (STCG) were determined. The STCG consist of 3,846 genes, a set that is much larger than that of the 2,882 Salmonella core genes (SCG) found previously. The STCG together with 1,576 core intergenic regions (IGRs) were defined as the S . Typhimurium core genome. Using 93 S . Typhimurium genomes from 13 epidemiologically confirmed community outbreaks, we demonstrated that typing based on the S . Typhimurium core genome (STCG plus core IGRs) provides superior resolution and higher discriminatory power than that based on SCG for outbreak investigation and molecular epidemiology of S . Typhimurium. STCG and STCG plus core IGR typing achieved 100% separation of all outbreaks compared to that of SCG typing, which failed to separate isolates from two outbreaks from background isolates. Defining the S . Typhimurium core genome allows standardization of genes/regions to be used for high-resolution epidemiological typing and genomic surveillance of S. Typhimurium.

Shiga toxin-producing Escherichia coli (STEC) O157 is an important foodborne pathogen. The aims of this study were to determine genetic relatedness of STEC O157 isolated from foods in China. STEC O157 isolates from food were characterized by virulence gene typing, antibiotyping, Multi-locus sequence typing (MLST), pulsed-field gel electrophoresis analysis (PFGE), Multi-locus variable number tandem repeat analysis (MLVA), and the single nucleotide polymorphism (SNP) clade typing. Of the 30 STEC O157 isolates analyzed, all isolates harbored eae, exhA, stx1 and/or stx2 genes with stx2c subtype predominating. By MLST, they were relatively homogenous with only 4 STs. PFGE and MLVA generated 22 pulsotyples and 23 patterns, respectively, which showed considerable diversity. Only one clade 8 isolate was detected. These results indicate that STEC O157 isolates from foods in China were heterogeneous. There was no correlation between genotypic characteristics and sources of isolates. Since different subtyping methods often give different discriminatory powers, the use of more than one subtyping approach is necessary in providing a more accurate picture of the genetic diversity of STEC O157. Four isolates were from ready-to-eat meat or salads, underscores the risk of infections. There is a need for surveillance of STEC O157 in foods and clinically and implementation of prevention strategies to prevent outbreaks of STEC O157 infections in China.

Shiga-toxin producing Escherichia coli (STEC) O157 is an important foodborne pathogen that can be transmitted to humans both directly and indirectly from the feces of beef cattle, its primary reservoir. Numerous studies have investigated the shedding dynamics of E. coli O157 by beef cattle, however the spatiotemporal trends of shedding are still not well understood. Molecular tools can increase the resolution through the use of strain typing to explore transmission dynamics within and between herds and identify strain-specific characteristics that may influence pathogenicity and spread. Previously, the shedding dynamics and molecular diversity, through the use of multilocus variable number of tandem repeat analysis (MLVA) of STEC O157, were separately investigated in an Australian beef herd over a 9-month study period. Variation in shedding was observed over time and 33 MLVA types were identified. The study presented here combines the two datasets previously published with an aim to clarify the relationship between epidemiological variables and strain types. Three major genetic clusters were identified that were significantly associated with the location of the cattle in different paddocks. No significant association between genetic clusters and individual cow was observed. Results from this molecular epidemiological study provide evidence for herd-level clonal replacement over time that may have been triggered by movement to a new paddock. In conclusion, this study has provided further insight into STEC O157 shedding dynamics and pathogen transmission. Knowledge gaps remain regarding the relationship of strain types and the shedding dynamics of STEC O157 by beef cattle that could be further clarified through the use of whole genome sequencing.

Salmonella enterica subsp enterica serovar Typhimurium (S. Typhimurium) is a serovar with broad host range. To determine the genomic diversity of S. Typhimurium, we sequenced 39 isolates (37 Australian and 2 UK isolates) representing 14 Repeats Groups (RGs) determined primarily by clustered regularly interspaced short palindromic repeats (CRISPR). Analysis of single nucleotide polymorphisms (SNPs) among the 39 isolates yielded an average of 1,232 SNPs per isolate, ranging from 128 SNPs to 11,339 SNPs relative to the reference strain LT2. Phylogenetic analysis of the 39 isolates together with 66 publicly available genomes divided the 105 isolates into five clades and 19 lineages, with the majority of the isolates belonging to clades I and II. The composition of CRISPR profiles correlated well with the lineages, showing progressive deletion and occasional duplication of spacers. Prophage genes contributed nearly a quarter of the S. Typhimurium accessory genome. Prophage profiles were found to be correlated with lineages and CRISPR profiles. Three new variants of HP2-like P2 prophage, several new variants of P22 prophage and a plasmid-like genomic island StmGI_0323 were found. This study presents evidence of horizontal transfer from other serovars or species and provides a broader understanding of the global genomic diversity of S. Typhimurium.

Backgroud: The emergence and spread of antimicrobial resistance has become a major global public health concern. A component element of this is the spread of the plasmid-encoded extended-spectrum β-lactamase (ESBL) genes, conferring resistance to third-generation cephalosporins. The purpose of this study was to investigate the molecular characteristics of ESBL-encoding genes identified in Escherichia coli cultured from diarrheic patients in China from 2013 to 2014. Methods: A total of 51 Escherichia coli were confirmed as ESBL producers by double-disk synergy testing of 912 Escherichia coli isolates studied. Polymerase chain reaction (PCR) and DNA sequencing were performed to identify the corresponding ESBL genes. Susceptibility testing was tested by the disk diffusion method. Plasmids were typed by PCR-based replicon typing and their sizes were determined by S1-nuclease pulsed-field gel electrophoresis. Multi-locus sequence typing (MLST) and phylogrouping were also performed. Broth mating assays were carried out for all isolates to determine whether the ESBL marker could be transferred by conjugation. Results: Of the 51 ESBL-positive isolates identified, blaCTX-M, blaTEM, blaOXA and blaSHV were detected in 51, 26, 3, 1 of these isolates respectively. Sequencing revealed that 7 blaCTX-M subtypes were detected, with blaCTX-M-14 being the most common, followed by blaCTX-M-79 and blaCTX-M-28. Of the 26 TEM-positive isolates identified, all of these were blaTEM-1 genotypes. All isolates contained one to three large plasmids and ten replicon types were detected. Of these, IncFrep (n = 50), IncK/B (n = 31), IncFIB (n = 26), IncB/O (n = 14), and IncI1-Ir (n = 8) replicon types were the predominating incompatibility groups. Twenty-six isolates demonstrated the ability to transfer their cefotaxime resistance marker at high transfer rates. MLST typing identified 31 STs and phylogenetic grouping showed that twelve of the 51 donor strains belonged to phylogroup B2. Conclusions: This study highlights the diversity of the ESBL producing E. coli and also the diversity of ESBL genes and plasmids carrying these genes in China, which poses a threat to public health.

Streptococcus suis , an important zoonotic pathogen, is a highly diverse species with only a subset of strains that cause disease in humans. Our previous study proposed a minimum core genome (MCG) sequence typing method and defined seven MCG groups, with MCG group 1 as the prevalent group causing human infections. In this study, we identified a set of 10 single nucleotide polymorphisms (SNPs) distributed in six genes that were used to identify the seven MCG groups. The 10 SNPs were typed for 179 S. suis isolates collected from slaughtered pigs. The most prevalent groups among the tested isolates were MCG groups 6 and 7. Most of the isolates (147/179) were genotyped as mrp negative, epf negative, sly negative, and CDS2157 positive. The 179 isolates were also typed by multilocus sequence typing (MLST) and divided into 115 sequence types (STs), 111 of which were new. The 6 serotypes (29, 11, 5, 12, 30, and 2) represented 72.3% of the serotyped isolates. Our data show that the typing assay facilitates the application of genome data to the surveillance of S. suis .

Salmonella enterica serovar Typhimurium is the most common Salmonella serovar causing human infections in Australia and many other countries. A total of 12,112 S. Typhimurium isolates from New South Wales were analyzed by multi-locus variable number of tandem repeat (VNTR) analysis (MLVA) using five VNTRs from 2007 to 2014. We found that mid ranges of repeat units of 8–14 in VNTR locus STTR5, 6–13 in STTR6, and 9–12 in STTR10 were always predominant in the population (>50%). In vitro passaging experiments using MLVA type carrying extreme length alleles found that the majority of long length alleles mutated to short ones and short length alleles mutated to longer ones. Both data suggest directional mutability of VNTRs toward mid-range repeats. Sequencing of 28 isolates from a newly emerged MLVA type and its five single locus variants revealed that single nucleotide variation between isolates with up to two MLVA differences ranged from 0 to 12 single nucleotide polymorphisms (SNPs). However, there was no relationship between SNP and VNTR differences. A population genetic model of the joint distribution of VNTRs and SNPs variations was used to estimate the mutation rates of the two markers, yielding a ratio of 1 VNTR change to 6.9 SNP changes. When only one VNTR repeat difference was considered, the majority of pairwise SNP difference between isolates were 4 SNPs or fewer. Based on this observation and our previous findings of SNP differences of outbreak isolates, we suggest that investigation of S. Typhimurium community outbreaks should include cases of 1 repeat difference to increase sensitivity. This study offers new insights into the short-term VNTR evolution of S. Typhimurium and its application for epidemiological typing.

Evolutionary studies using single nucleotide polymorphisms (SNPs) have separated Bordetella pertussis isolates into six major clusters, with recent isolates forming cluster I. The expansion of cluster I isolates was characterised by changes in genes encoding antigenic components in acellular vaccines, including pertactin (Prn). Here, we determined the initial emergence of the pertussis toxin promoter allele, ptxP3, from an evolutionary perspective. This allele was previously shown in a study from the Netherlands to be associated with increased pertussis toxin production as a result of a single base mutation in the ptxP. The ptxP region of 313 worldwide isolates was sequenced, including 208 isolates from Australia collected over a 40 year period. Eight alleles were identified, of which only two predominated: ptxP1 and ptxP3. One novel allele was also found. ptxP3 was only found in SNP cluster I of B. pertussis and its emergence is concurrent with the change to the non-vaccine prn2 allele. Our results suggest that the globally distributed cluster I of B. pertussis has the ability to evade vaccine induced selection pressure.

Despite high vaccine coverage, pertussis incidence has increased substantially in recent years in many countries. A significant factor that may be contributing to this increase is adaptation to the vaccine by Bordetella pertussis, the causative agent of pertussis. In this study, we first assessed the genetic diversity of B. pertussis by microarray-based comparative genome sequencing of 10 isolates representing diverse genotypes and different years of isolation. We discovered 171 single nucleotide polymorphisms (SNPs) in a total of 1.4 Mb genome analyzed. The frequency of base changes was estimated as one per 32 kb per isolate, confirming that B. pertussis is one of the least variable bacterial pathogens. We then analyzed an international collection of 316 B. pertussis isolates using a subset of 65 of the SNPs and identified 42 distinct SNP profiles (SPs). Phylogenetic analysis grouped the SPs into six clusters. The majority of recent isolates belonged to clusters I-IV and were descendants of a single prevaccine lineage. Cluster I appeared to be a major clone with a worldwide distribution. Typing of genes encoding acellular vaccine (ACV) antigens, ptxA, prn, fhaB, fim2, and fim3 revealed the emergence and increasing incidence of non-ACV alleles occurring in clusters I and IV, which may have been driven by ACV immune selection. Our findings suggest that B. pertussis, despite its high population homogeneity, is evolving in response to vaccination pressure with recent expansion of clones carrying variants of genes encoding ACV antigens.

Multilocus variable-number tandem-repeat analysis (MLVA) of 316 Bordetella pertussis isolates collected over 40 years from Australia and 3 other continents identified 66 MLVA types (MTs), including 6 predominant MTs. Typing of genes encoding acellular vaccine antigens showed changes that may be vaccine driven in 2 MTs prevalent in Australia.

The genus Salmonella consists of two species S. enterica and S. bongori. S. enterica has a well defined subspecies structure with seven subspecies consistently delineated by sequence variation. Frequency of recombination between subspecies and within a subspecies is markedly different. Subspecies I undergoes frequent recombination as demonstrated recently, demystifying the long-held belief that Salmonella is a highly clonal organism. The majority of disease causing serovars are from subspecies I with the most important serovars in human health being Typhimurium and Typhi. Typhimurium has developed considerable diversity and may be a very old serovar. The majority of the isolates belong to a single clonal complex by multilocus sequence typing. Typhimurium isolates are divided into phage types and some of the phage types do not have a single origin as determined using mutational changes. Phage type DT104 is heterogeneous and represented in multiple sequence types, with its multidrug-resistant variant most successful causing epidemics in many parts of the world. Typhi, a human restricted serovar, is relatively young compared to Typhimurium, and has a low level of sequence variation. Single nucleotide polymorphisms (SNPs) have been shown to be very useful for typing and resolving relationships within Typhi. Genome sequences of 19 isolates revealed more than 1700 SNPs. The fully resolved phylogenetic tree allows one to trace the mutational changes occurred during clonal diversification. Genome wide SNPs have greatly enhanced our understanding of the evolution of Salmonella clones.

Multilocus variable-number tandem repeats (VNTRs) are widely used as molecular markers to differentiate isolates of homogenous pathogenic clones. We explored the genomes of Salmonella enterica serovar Typhi strains CT18 and Ty2 for potential VNTRs. Among the 43 potential VNTRs screened, 2 were found to be polymorphic. Together with seven polymorphic VNTRs from previous studies, they were used to type 73 global serovar Typhi isolates. A total of 70 multilocus VNTR analysis (MLVA) profiles were found, distinguishing all except three pairs of isolates into individual profiles. The discriminatory power was 0.999. Phylogenetic analysis showed that the MLVA profiles can be divided into seven clusters. However, except for the closely related isolates, the relationships derived were in conflict with those inferred from single nucleotide polymorphism (SNP) typing using 38 SNPs done previously. We concluded that MLVA can resolve the relationships only among closely related isolates. A combination of SNP typing and MLVA typing offers the best approach for local and global epidemiology and the evolutionary analysis of serovar Typhi. We suggest that seven of the nine most polymorphic VNTRs be used as a standardized typing scheme for epidemiological typing.

The Vibrio pathogenicity island (VPI) encodes the toxin-coregulated pilus and other virulence factors for Vibrio cholerae to colonize the human intestine to cause cholera. We assessed the level of genetic variation of VPI in nine nonpandemic isolates, and compared them with the sixth and seventh pandemic strains by sequencing c. 5 kb each from the start, middle and end regions of the VPI. Variation is similar among the three regions at around 2%, except for the tcpA gene, which has a much higher level of variation (23%). Numerous recombination segments were identified with sizes up to 2177 bp. Nearly all VPI genes sequenced have a ratio of synonymous to nonsynonymous substitutions considerably lower than that for housekeeping genes, suggesting that VPI genes are under positive selection pressure for change. The tagA gene was deleted or damaged in six isolates, which is likely to affect the efficiency of colonization of the human intestine. Two genes, orf2 and acfD, previously found to be translated differently in the sixth and seventh pandemic strains, were determined to be mutant in the seventh and sixth pandemic strains, respectively. These findings enhance our understanding of variation in the VPI, and of the pathogenic potential of VPI-positive environmental isolates.

Shigella flexneri, which causes shigellosis in humans, evolved from Escherichia coli. The sequencing of Shigella genomes has revealed that a large number of insertion sequence (IS) elements (over 200 elements) reside in the genome. Although the presence of these elements has been noted previously and summarized, more detailed analyses are required to understand their evolutionary significance. Here, the genome of S. flexneri strain 2457T is used to investigate the spatial distribution of IS copies around the chromosome and the location of elements with respect to genes. It is found that most IS isoforms occur essentially randomly around the genome. Two exceptions are IS91 and IS911, which appear to cluster due to local hopping. The location of IS elements with respect to genes is biased, however, revealing the action of natural selection. The non-coding regions of the genome (no more than 21%) carry disproportionally more IS elements (at least 28%) than the coding regions, implying that selection acts against insertion into genes. Of the genes disrupted by ISs, those involved in signal transduction, intracellular trafficking, and cell motility are most commonly targeted, suggesting selection against genes in these categories.

Summary Objectives Despite high pertussis vaccination coverage, Australia experienced a prolonged epidemic in 2008–2012. The predominant Bordetella pertussis genotype harboured pertussis toxin promoter allele, ptxP3 , and pertactin gene allele, prn2 . The emergence and expansion of prn non-expressing isolates (Prn negative), were also observed. We aimed to investigate the microevolution and genomic diversity of epidemic B. pertussis isolates. Methods We sequenced 22 B . pertussis isolates collected in 2008–2012 from two states of Australia which are geographically widely separated. Ten of the 22 were Prn negative isolates with three different modes of silencing of prn ( prn ::IS 481 F, prn ::IS 481 R and prn ::IS 1002 ). Five pre-epidemic isolates were also sequenced for comparison. Results Five single nucleotide polymorphisms were common in the epidemic isolates and differentiated them from pre-epidemic isolates. The Australian epidemic isolates can be divided into five lineages (EL1–EL5) with EL1 containing only Prn negative isolates. Comparison with global isolates showed that three lineages remained geographically and temporally distinct whereas two lineages mixed with isolates from 2012 UK outbreak. Conclusion Our results suggest significant diversification and the microevolution of B. pertussis within the 2008–2012 Australian epidemic.

Campylobacter concisus is an oral bacterium that has been shown to be associated with inflammatory bowel disease (IBD). In this study we examined clusters of oral C. concisus strains isolated from patients with IBD and healthy controls by analysing six housekeeping genes. In addition, we investigated the population structure of C. concisus strains. Whether oral and enteric strains form distinct clusters based on the sequences of these housekeeping genes was also investigated. The oral C. concisus strains were found to contain two genomospecies, which belong to the two genomospecies previously found in enteric C. concisus strains. C. concisus clusters formed based on the sequences of a single aspA gene were the same as that formed by using previously reported MLST schemes. The analysis of combined oral and enteric C. concisus strains found that enteric C. concisus strains did not form distinct clusters. Genetic structure analysis identified five subpopulations of C. concisus and showed that genetic recombination between C. concisus strains was common. However, genetic recombination was significantly less in oral strains isolated from patients with IBD than from healthy individuals. Previously reported oral and enteric intestinal epithelial invasive C. concisus strains were in cluster II and subpopulation III. Furthermore, this study shows that there are no distinct enteric C. concisus strain clusters or subpopulations.

Whole-genome next-generation sequencing (NGS) was used to retrospectively examine 57 isolates from five epidemiologically confirmed community outbreaks (numbered 1 to 5) caused by Salmonella enterica serovar Typhimurium phage type DT170. Most of the human and environmental isolates confirmed epidemiologically to be involved in the outbreaks were either genomically identical or differed by one or two single nucleotide polymorphisms (SNPs), with the exception of those in outbreak 1. The isolates from outbreak 1 differed by up to 12 SNPs, which suggests that the food source of the outbreak was contaminated with more than one strain while each of the other four outbreaks was caused by a single strain. In addition, NGS analysis ruled in isolates that were initially not considered to be linked with the outbreak, which increased the total outbreak size by 107%. The mutation process was modeled by using known mutation rates to derive a cutoff value for the number of SNP difference to determine whether or not a case was part of an outbreak. For an outbreak with less than 1 month of ex vivo / in vivo evolution time, the maximum number of SNP differences between isolates is two or four using the lowest or highest mutation rate, respectively. NGS of S . Typhimurium significantly increases the resolution of investigations of community outbreaks. It can also inform a more targeted public health response by providing important supplementary evidence that cases of disease are or are not associated with food-borne outbreaks of S . Typhimurium.

Salmonella enterica serovar Typhimurium is a food-borne pathogen and a leading cause of gastroenteritis in humans. Recently, we sequenced a phage-type DT108 strain (L945) and found reads with high similarity to both Salmonella typhi strain CT18 plasmid pHCM2 and bacteriophage SSU5. In this study, we completely sequenced the novel phage-like plasmid which was designated as pSTM\_Φ. The presence of this phage-like plasmid was examined in a collection of 284 Salmonella Typhimurium isolates using PCR of the parB gene and only one other isolate (L946) was found to carry the phage-like plasmid suggesting that it is infrequently present amongst Salmonella Typhimurium isolates. pSTM\_Φ is a circular phage-like plasmid of 107.7 kb encoding 132 coding regions (ORFs) with the majority of the ORFs encoding hypothetical proteins. Comparative analysis with other closely related phage-like plasmids and the SSU5 phage revealed that there were four divergent lineages of phage-like plasmids found in the family of Enterobacteriaceae . In conclusion, pSTM_Φ is a new member of an emerging family of phage-like plasmids.

Streptococcus suis is an important pathogen of pigs and may cause serious disease in humans. Serotyping is one of the important diagnostic tools and is used for the epidemiological study of S. suis . Nontypeable S. suis strains have been reported in many studies; however, the capsular polysaccharide (CPS) synthesis cps loci of nontypeable strains have not been analyzed. In this study, we investigated the genetic characteristics of cps loci in 78 nontypeable strains isolated from healthy pigs. Eight novel cps loci (NCLs) were found, and all of them were located between the orfZ-orfX region and the glf gene. All NCLs possess the wzy and wzx genes, strongly suggesting that the CPSs of these NCLs were synthesized using the Wzx/Wzy-dependent pathway. The cps genes found in the 78 isolates were assigned to 96 homology groups (HGs), 55 of which were NCL specific. The encapsulation of the 78 isolates was also examined using transmission electron microscopy. Fifty-three isolates were found to have a capsule, and these were of varied thicknesses. Our data enhance our understanding of the cps gene cluster diversity of nontypeable S. suis strains and provide insight into the evolution of the S. suis capsular genes.

Summary Pathogenic Escherichia coli cause over 160 million cases of dysentery and one million deaths per year, whereas non-pathogenic E. coli constitute part of the normal intestinal flora of healthy mammals and birds. The evolutionary pathways underlying this dichotomy in bacterial lifestyle were investigated by multilocus sequence typing of a global collection of isolates. Specific pathogen types (enterohaemorrhagic E. coli , enteropathogenic E. coli , enteroinvasive E. coli , K1 and Shigella ) have arisen independently and repeat- edly in several lineages, whereas other lineages con- tain only few pathogens. Rates of evolution have accelerated in pathogenic lineages, culminating in highly virulent organisms whose genomic contents are altered frequently by increased rates of homolo- gous recombination; thus, the evolution of virulence is linked to bacterial sex. This long-term pattern of evolution was observed in genes distributed through- out the genome, and thereby is the likely result of episodic selection for strains that can escape the host immune response.

The genetic relationship and population structure of Salmonella enterica subspecies I strains were analysed using nucleotide sequences of four genes (mglA, proV, torC and speC). Fifteen strains from the Salmonella reference collection B (SARB), belonging to 13 serovars, were analysed. Sequence data of two housekeeping genes, mdh and mutS, of the same 15 strains reported by Brown et al. (2003) (Proc Natl Acad Sci U S A 100, 15676-15681) were also included in the analyses. Phylogenetic analysis revealed that there was a lack of congruence among the six gene trees. Split decomposition analysis resolved only five strains with a network structure, while others showed a star phylogeny. Compatibility values for the SARB strains were the lowest in comparison to those for strains representing different subspecies of S. enterica. These results showed that the genes studied have undergone frequent recombination, suggesting a low level of clonality within subspecies I of S. enterica.