Hierarchical and spatially explicit clustering of DNA sequences with BAPS software

Population genomics of Group B Streptococcus reveals the genetics of neonatal disease onset and meningeal invasion

Group B Streptococcus (GBS), or Streptococcus agalactiae, is a pathogen that causes preterm births, stillbirths, and acute invasive neonatal disease burden and mortality. Here, we investigate bacterial genetic signatures associated with disease onset time and meningeal tissue infection in acute invasive neonatal GBS disease. We carry out a genome-wide association study (GWAS) of 1,338 GBS isolates from newborns with acute invasive disease; the isolates had been collected annually, for 30 years, through a national bacterial surveillance program in the Netherlands. After controlling for the population structure, we identify genetic variation within noncoding and coding regions, particularly the capsule biosynthesis locus, statistically associated with neonatal GBS disease onset time and meningeal invasion. Our findings highlight the impact of integrating microbial population genomics and clinical pathogen surveillance, and demonstrate the effect of GBS genetics on disease pathogenesis in neonates and infants.

Traveler Mites: Population Genetic Structure of the Wing Mites Periglischrus paracaligus (Acari: Mesostigmata: Spinturnicidae)

Wing mites of the genus Periglischrus are ectoparasites exclusively associated with phyllostomid bats. These mites show high host specificity and have been studied to understand the evolutionary history of their bat hosts mainly by using a morphological variation. Through a phylogeographic approach, we analyzed the genetic diversity and population genetic structure of the ectoparasite Periglischrus paracaligus Herrin and Tipton which parasitizes Leptonycteris yerbabuenae Martínez and Villa (lesser long-nosed bat) in Mexico. By the implementation of a multilocus approach, we found that P. paracaligus populations were diverse for haplotype diversity, and had values ranging from 0.5 to 1. No genetic structuring in the P. paracaligus parasites was observed along with the distribution of the host, L. yerbabuenae, in Mexico, nor when populations or regions were compared, but our results revealed a process of historical demographic expansion in all the analyzed markers. We discuss possible scenarios that could explain the lack of population structure in the light of the data analyzed for the parasites and the biology of L. yerbabuenae, such as the interplay between parasite and host traits being responsible for the genetic make-up of parasite populations. We also inferred its phylogenetic position among wing mites parasitizing the two other species of Leptonycteris bats. Long-nosed bats' monophyly helps to explain the observed presence of distinctive clades in the wing mite's phylogeny in specific association with each long-nosed bat host species.

Use of Whole-Genome Sequencing to Unravel the Genetic Diversity of a Prevalent Mycobacterium bovis Spoligotype in a Multi-Host Scenario in Spain

Despite the efforts invested in the eradication of bovine tuberculosis in Spain, herd prevalence has remained constant in the country during the last 15 years (~1.5–1.9%) due to a combination of epidemiological factors impairing disease control, including between-species transmission. Here, our aim was to investigate the molecular diversity of Mycobacterium bovis isolates belonging to the highly prevalent SB0339 spoligotype in the cattle-wildlife interface in different regions of Spain using whole-genome sequencing (WGS). Genomic data of 136 M. bovis isolates recovered from different animal species (cattle, wild boar, fallow deer, and red deer) and locations between 2005 and 2018 were analyzed to investigate between- and within-species transmission, as well as within-herds. All sequenced isolates differed by 49–88 single nucleotide polymorphisms from their most recent common ancestor. Genetic heterogeneity was geographic rather than host species-specific, as isolates recovered from both cattle and wildlife from a given region were more closely related compared to isolates from the same species but geographically distant. In fact, a strong association between the geographic and the genetic distances separating pairs of M. bovis isolates was found, with a significantly stronger effect when cattle isolates were compared with wildlife or cattle-wildlife isolates in Spain. The same results were obtained in Madrid, the region with the largest number of sequenced isolates, but no differences depending on the host were observed. Within-herd genetic diversity was limited despite the considerable time elapsed between isolations. The detection of closely related strains in different hosts demonstrates the complex between-host transmission dynamics present in endemic areas in Spain. In conclusion, WGS results a valuable tool to track bTB infection at a high resolution and may contribute to achieve its eradication in Spain.

Convergence of MCR-8.2 and Chromosome-Mediated Resistance to Colistin and Tigecycline in an NDM-5-Producing ST656 Klebsiella pneumoniae Isolate From a Lung Transplant Patient in China

We characterized the first NDM-5 and MCR-8.2 co-harboring ST656 Klebsiella pneumoniae clinical isolate, combining with chromosomal gene-mediated resistance to colistin and tigecycline. The K. pneumoniae KP32558 was isolated from the bronchoalveolar lavage fluid from a lung transplant patient. Complete genome sequences were obtained through Illumina HiSeq sequencing and nanopore sequencing. The acquired resistance genes and mutations in chromosome-encoded genes associated with colistin and tigecycline resistance were analyzed. Comparative genomic analysis was conducted between mcr-8.2-carrying plasmids. The K. pneumoniae KP32558 was identified as a pan-drug resistant bacteria, belonging to ST656, and harbored plasmid-encoded blaNDM-5 and mcr-8.2 genes. The blaNDM-5 gene was located on an IncX3 type plasmid. The mcr-8.2 gene was located on a conjugative plasmid pKP32558-2-mcr8, which had a common ancestor with another two mcr-8.2-carrying plasmids pMCR8_020135 and pMCR8_095845. The MIC of KP32558 for colistin was 256 mg/L. The mcr-8.2 gene and mutations in the two-component system, pmrA and crrB, and the regulator mgrB, had a synergistic effect on the high-level colistin resistance. The truncation in the acrR gene, related to tigecycline resistance, was also identified. K. pneumoniae has evolved a variety of complex resistance mechanisms to the last-resort antimicrobials, close surveillance is urgently needed to monitor the prevalence of this clone.

Phylogeography of Nasutitermes ephratae (Termitidae: Nasutitermitinae) in neotropical region

The neotropical region ranks third in the number of termites and includes five different families. Of these, Termitidae is the most diverse and includes the species Nasutitermes ephratae, which is widespread in the neotropics. To date, only one study has been published about phylogeography in neotropical termites (N. corniger). Here, we explored the population genetic patterns of N. ephratae and also evaluated the phylogeographical processes involved in the evolutionary history of the species. We used the mitochondrial genes 16S rRNA and COII as molecular markers: these were sequenced for 128 samples of N. ephratae. We estimated the genetic diversity and divergence time as well as the demography and genetic structure. We also performed an ancestral area reconstruction and a haplotype network. The results showed high genetic variability, recent demographic expansion, and strong genetic structure. A dispersal route for the species, that occurred in both directions between South and Central America, was inferred. The results emphasize a temporary separation between the South and Central America populations that affected the origin of the current Central America populations. These populations were formed from different phylogeographic histories.

Genetic analyses reveal regional structure and demographic expansion of the predominant tea pest Empoasca onukii (Hemiptera: Cicadellidae) in China

BACKGROUND

The tea green leafhopper, Empoasca onukii Matsuda, is the most destructive insect pest of tea plantations in East Asia. Despite its economic importance and previous studies on this species, it remains unclear as to how this small-sized pest can have such wide range.

RESULTS

By sequencing three mitochondrial genes and 17 microsatellite loci, we revealed the regional structure and demographic expansion of 59 E. onukii populations in China. Bayesian analysis of population genetic structure (BAPS) on microsatellites identified four genetic groups with spatial discontinuities, while analysis on mitochondrial genes inferred five nested and differentiated clusters. Both the Mantel test and the generalized linear model indicated a significant pattern of isolation by geographic distance in E. onukii populations. Based on the approximate Bayesian computation approach, E. onukii was found to have originated from southwestern China and expanded northward and eastward. While MIGRATE-N and Bayesian stochastic search variable selection (BSSVS) procedure in BEAST confirmed the possible eastward and northward dispersal from Yunnan, they also detected more gene flow from the derived populations in central and southeastern China.

CONCLUSION

Our results suggest that the current distribution and structure of E. onukii is complicatedly influenced by human activities of cultivation, wide dissemination of tea in ancient China as well as recent transportation of tea seedlings for establishing new tea plantations. Insights into genetic differentiation and demographic expansion patterns from this study provide an important basis for the development of area-wide management of the E. onukii populations. © 2022 Society of Chemical Industry.

Epidemiological and Genetic Characteristics of Clinical Carbapenem-Resistant Acinetobacter baumannii Strains Collected Countrywide from Hospital Intensive Care Units (ICUs) in China

Acinetobacter baumannii is one of the key Gram-negative pathogens that can cause serious nosocomial infections. In China, a large proportion of clinical A. baumannii strains are multidrug resistant, among which strains resistant to carbapenems are particularly worrisome, as infections caused by such strains may limit the choice of existing antibiotics. We conducted a nationwide and genome-based surveillance on the prevalence and antibiotic susceptibility profile of carbapenem-resistant A. baumannii (CRAB) strains collected from intensive care units (ICUs) in hospitals in different provinces and investigated the routes of transmission and mechanism of resistance by whole-genome sequencing and phylogenetic analysis. We found that CRAB strains were prevalent in 71.4% (55/77) of the ICUs surveyed. Clonal spread of CRAB was found in 37.6% (29/77) of ICUs and a total of 22 different clones were identified. Most clones were transmissible within one ICU, but up to six clones could be detected in at least three hospitals. In addition, carbapenem-hydrolysing class D β-lactamases (CHDL) were found to be mainly responsible for carbapenem-resistance in A. baumannii and the ST2 global-clone is the predominant type of CRAB in China. Importantly, we found that CRAB isolates currently exhibited extremely low rate of resistance to colistin (0.4%) and tigecycline (2.5%), but a high rate of resistance to ceftazidime-avibactam (70.2%). Findings in this work shall facilitate development of appropriate antimicrobial regimens for treatment of CRAB infections. Further surveillance and research on the evolutionary and epidemiological features of clinical CRAB strains are necessary.

Epidemiology of blaCTX-M-Positive Salmonella Typhimurium From Diarrhoeal Outpatients in Guangdong, China, 2010–2017

Salmonella enterica can lead to intestinal diarrhea, and the emergence and spread of cephalosporin-resistant Salmonella have brought great challenges to clinical treatment. Therefore, this study investigated the prevalence and transmission of bla_CTX-M genes among S. Typhimurium from diarrhoeal outpatients in Guangdong, China, from 2010 to 2017. A total of 221 bla_CTX-M-positive isolates were recovered from 1,263 S. Typhimurium isolates from the facal samples of diarrhoea patients in 45 general hospitals from 11 cities. The most popular CTX-M gene was bla_CTX-M-55 (39.6%, 72/182) in the CTX-M-1 group, followed by bla_CTX-M-14 (22.5%, 41/182) and bla_CTX-M-65 (19.2%, 35/182) in the CTX-M-9 group. The isolates that carried bla_CTX-M-9G had significantly higher resistance rates to multiple antibacterials compared with bla_CTX-M-1G (p < 0.01). Meanwhile, PFGE analysis not only showed the clonal transmission of bla_{CTX-M-55/14/65}-positve isolates of diarrhoeal outpatients’ origins from different hospitals in Guangdong province, but also the characteristic of bla_{CTX-M-55/14/65}-positve isolates’ bacterial persistence. Multilocus sequence typing (MLST) analysis indicated that these S. Typhimurium isolates possessed ST34 and ST19. Furthermore, genomic Beast phylogenomic analysis provided the evidence of a close relationship of bla_CTX-M-positive S. Typhimurium isolates between the outpatients and pork. Most bla_{CTX-M-55/14/65} genes were transmitted by non-typeable or IncI1/IncFII/IncHI2 plasmids with the size of ranging from ~80 to ~280 kb. Moreover, whole-genome sequencing (WGS) analysis further revealed that bla_{CTX-M-55/14/65} coexisted with other 25 types of ARGs, of which 11 ARGs were highly prevalent with the detection rates >50%, and it first reported the emergence of bla_TEM-141 in S. Typhimurium. This study underscores the importance of surveillance for bla_CTX-M-positive microbes in diarrhea patients.

Evolution of an Apomixis-Specific Allele Class in Supernumerary Chromatin of Apomictic Boechera

Asexual reproduction through seeds in plants (i.e., apomixis) is a heritable trait, and apomixis- linked loci have been identified in multiple species. However, direct identification of genomic elements is typically hindered as apomixis-linked loci and are commonly found in recombination-suppressed and repetitive regions. Heterochromatinized elements, such as B chromosomes and other supernumerary chromosomal DNA fragments have long been known to be associated with asexuality in both plants and animals and are prime candidate regions for the evolution of multiple apomixis factors controlling the individual elements of apomixis. Here, we examined molecular evolution, gene regulation, and chromosomal location of a male apomeiosis factor (UPG2), a long noncoding RNA gene, in sexual and apomictic Boechera with and without male apomeiosis (i.e., balanced and unbalanced apomicts). We revealed the origin of the gene in the apomixis genome on an apomixis-specific, supernumerary heterochromatic Boechera chromosome (Boe1). The UPG2 is active in the tapetum at male meiosis. We found allele classes specific to apomictic and sexual Boechera accessions and a third class that shares the features of both and points to a convergent transition state. Sex alleles are found only in some of the sexual accessions and have higher nucleotide divergence and lower transcriptional activity compared to apo alleles. These data demonstrate selective pressure to maintain the function of UPG2 for unreduced pollen formation in apomicts as the occasional transmission of the allele from unbalanced apomicts into sexual organisms that lead to pseudogenization and functional decay of copies in sexual organisms.

Helicobacter cinaedi is a human-adapted lineage in the Helicobacter cinaedi/canicola/'magdeburgensis' complex

Helicobacter cinaedi is an enterohepatic Helicobacter that causes bacteremia and other diseases in humans. While H. cinaedi-like strains are isolated from animals, including dog isolates belonging to a recently proposed H. canicola, little is known about the genetic differences between H. cinaedi and these animal isolates. Here, we sequenced 43 H. cinaedi- or H. canicola-like strains isolated from humans, hamsters, rats and dogs and collected 81 genome sequences of H. cinaedi, H. canicola and other enterohepatic Helicobacter strains from public databases. Genomic comparison of these strains identified four distinct clades (clades I–IV) in H. cinaedi/canicola/‘magderbugensis’ (HCCM) complex. Among these, clade I corresponds to H. cinaedi sensu stricto and represents a human-adapted lineage in the complex. We identified several genomic features unique to clade I. They include the accumulation of antimicrobial resistance-related mutations that reflects the human association of clade I and the larger genome size and the presence of a CRISPR-Cas system and multiple toxin-antitoxin and restriction-modification systems, both of which indicate the contribution of horizontal gene transfer to the evolution of clade I. In addition, nearly all clade I strains but only a few strains belonging to one minor clade contained a highly variable genomic region encoding a type VI secretion system (T6SS), which could play important roles in gut colonization by killing competitors or inhibiting their growth. We also developed a method to systematically search for H. cinaedi sequences in large metagenome data sets based on the results of genome comparison. Using this method, we successfully identified multiple HCCM complex-containing human faecal metagenome samples and obtained the sequence information covering almost the entire genome of each strain. Importantly, all were clade I strains, supporting our conclusion that H. cinaedi sensu stricto is a human-adapted lineage in the HCCM complex.

Evolution and dispersal of mitochondrial DNA haplogroup U5 in Northern Europe: insights from an unsupervised learning approach to phylogeography

Background

We combined an unsupervised learning methodology for analyzing mitogenome sequences with maximum likelihood (ML) phylogenetics to make detailed inferences about the evolution and diversification of mitochondrial DNA (mtDNA) haplogroup U5, which appears at high frequencies in northern Europe.

Methods

Haplogroup U5 mitogenome sequences were gathered from GenBank. The hierarchal Bayesian Analysis of Population Structure (hierBAPS) method was used to generate groups of sequences that were then projected onto a rooted maximum likelihood (ML) phylogenetic tree to visualize the pattern of clustering. The haplogroup statuses of the individual sequences were assessed using Haplogrep2.

Results

A total of 23 hierBAPS groups were identified, all of which corresponded to subclades defined in Phylotree, v.17. The hierBAPS groups projected onto the ML phylogeny accurately clustered all haplotypes belonging to a specific haplogroup in accordance with Haplogrep2. By incorporating the geographic source of each sequence and subclade age estimates into this framework, inferences about the diversification of U5 mtDNAs were made. Haplogroup U5 has been present in northern Europe since the Mesolithic, and spread in both eastern and western directions, undergoing significant diversification within Scandinavia. A review of historical and archeological evidence attests to some of the population interactions contributing to this pattern.

Conclusions

The hierBAPS algorithm accurately grouped mitogenome sequences into subclades in a phylogenetically robust manner. This analysis provided new insights into the phylogeographic structure of haplogroup U5 diversity in northern Europe, revealing a detailed perspective on the diversity of subclades in this region and their distribution in Scandinavian populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08572-y.

Climate Warming Since the Holocene Accelerates West-East Communication for the Eurasian Temperate Water Strider Species Aquarius paludum

Holocene climate warming has dramatically altered biological diversity and distributions. Recent human-induced emissions of greenhouse gases will exacerbate global warming and thus induce threats to cold-adapted taxa. However, the impacts of this major climate change on transcontinental temperate species are still poorly understood. Here, we generated extensive genomic datasets for a water strider, Aquarius paludum, which was sampled across its entire distribution in Eurasia and used these datasets in combination with ecological niche modeling (ENM) to elucidate the influence of the Holocene and future climate warming on its population structure and demographic history. We found that A. paludum consisted of two phylogeographic lineages that diverged in the middle Pleistocene, which resulted in a “west–east component” genetic pattern that was probably triggered by Central Asia-Mongoxin aridification and Pleistocene glaciations. The diverged western and eastern lineages had a second contact in the Holocene, which shaped a temporary hybrid zone located at the boundary of the arid–semiarid regions of China. Future predictions detected a potentially novel northern corridor to connect the western and eastern populations, indicating west–east gene flow would possibly continue to intensify under future warming climate conditions. Further integrating phylogeographic and ENM analyses of multiple Eurasian temperate taxa based on published studies reinforced our findings on the “west–east component” genetic pattern and the predicted future northern corridor for A. paludum. Our study provided a detailed paradigm from a phylogeographic perspective of how transcontinental temperate species differ from cold-adapted taxa in their response to climate warming.

Effect of childhood vaccination and antibiotic use on pneumococcal populations and genome-wide associations with disease among children in Nepal: an observational study

Background

Methods

Findings

Interpretation

Funding

Dissecting the dominant hot spring microbial populations based on community-wide sampling at single-cell genomic resolution

With advances in DNA sequencing and miniaturized molecular biology workflows, rapid and affordable sequencing of single-cell genomes has become a reality. Compared to 16S rRNA gene surveys and shotgun metagenomics, large-scale application of single-cell genomics to whole microbial communities provides an integrated snapshot of community composition and function, directly links mobile elements to their hosts, and enables analysis of population heterogeneity of the dominant community members. To that end, we sequenced nearly 500 single-cell genomes from a low diversity hot spring sediment sample from Dewar Creek, British Columbia, and compared this approach to 16S rRNA gene amplicon and shotgun metagenomics applied to the same sample. We found that the broad taxonomic profiles were similar across the three sequencing approaches, though several lineages were missing from the 16S rRNA gene amplicon dataset, likely the result of primer mismatches. At the functional level, we detected a large array of mobile genetic elements present in the single-cell genomes but absent from the corresponding same species metagenome-assembled genomes. Moreover, we performed a single-cell population genomic analysis of the three most abundant community members, revealing differences in population structure based on mutation and recombination profiles. While the average pairwise nucleotide identities were similar across the dominant species-level lineages, we observed differences in the extent of recombination between these dominant populations. Most intriguingly, the creek's Hydrogenobacter sp. population appeared to be so recombinogenic that it more closely resembled a sexual species than a clonally evolving microbe. Together, this work demonstrates that a randomized single-cell approach can be useful for the exploration of previously uncultivated microbes from community composition to population structure.

Genomic Epidemiology of Global Carbapenemase-Producing Escherichia coli, 2015-2017

We describe the global molecular epidemiology of 229 carbapenemase-producing Escherichia coli in 36 countries during 2015-2017. Common carbapenemases were oxacillinase (OXA) 181 (23%), New Delhi metallo-β-lactamase (NDM) 5 (20%), OXA-48 (17%), Klebsiella pneumoniae carbapenemase 2 (15%), and NDM-1 (10%). We identified 5 dominant sequence types (STs); 4 were global (ST410, ST131, ST167, and ST405), and 1 (ST1284) was limited to Turkey. OXA-181 was frequent in Jordan (because of the ST410-B4/H24RxC subclade) and Turkey (because of ST1284). We found nearly identical IncX3-bla_OXA-181 plasmids among 11 STs from 12 countries. NDM-5 was frequent in Egypt, Thailand (linked with ST410-B4/H24RxC and ST167-B subclades), and Vietnam (because of ST448). OXA-48 was common in Turkey (linked with ST11260). Global K. pneumoniae carbapenemases were linked with ST131 C1/H30 subclade and NDM-1 with various STs. The global carbapenemase E. coli population is dominated by diverse STs with different characteristics and varied geographic distributions, requiring ongoing genomic surveillance.

Sub-Lineage Specific Phenolic Glycolipid Patterns in the Mycobacterium tuberculosis Complex Lineage 1

“Ancestral” Mycobacterium tuberculosis complex (MTBC) strains of Lineage 1 (L1, East African Indian) are a prominent tuberculosis (TB) cause in countries around the Indian Ocean. However, the pathobiology of L1 strains is insufficiently characterized. Here, we used whole genome sequencing (WGS) of 312 L1 strains from 43 countries to perform a characterization of the global L1 population structure and correlate this to the analysis of the synthesis of phenolic glycolipids (PGL) – known MTBC polyketide-derived virulence factors. Our results reveal the presence of eight major L1 sub-lineages, whose members have specific mutation signatures in PGL biosynthesis genes, e.g., pks15/1 or glycosyltransferases Rv2962c and/or Rv2958c. Sub-lineage specific PGL production was studied by NMR-based lipid profiling and strains with a completely abolished phenolphthiocerol dimycoserosate biosynthesis showed in average a more prominent growth in human macrophages. In conclusion, our results show a diverse population structure of L1 strains that is associated with the presence of specific PGL types. This includes the occurrence of mycoside B in one sub-lineage, representing the first description of a PGL in an M. tuberculosis lineage other than L2. Such differences may be important for the evolution of L1 strains, e.g., allowing adaption to different human populations.

Distinct increase in antimicrobial resistance genes among Escherichia coli during 50 years of antimicrobial use in livestock production in China

Antimicrobial use in livestock production is linked to the emergence and spread of antimicrobial resistance (AMR), but large-scale studies on AMR changes in livestock isolates remain scarce. Here we applied whole-genome sequence analysis to 982 animal-derived Escherichia coli samples collected in China from the 1970s to 2019, finding that the number of AMR genes (ARGs) per isolate doubled-including those conferring resistance to critically important agents for both veterinary (florfenicol and norfloxacin) and human medicine (colistin, cephalosporins and meropenem). Plasmids of incompatibility groups IncC, IncHI2, IncK, IncI and IncX increased distinctly in the past 50 years, acting as highly effective vehicles for ARG spread. Using antimicrobials of the same class, or even unrelated classes, may co-select for mobile genetic elements carrying multiple co-existing ARGs. Prohibiting or strictly curtailing antimicrobial use in livestock is therefore urgently needed to reduce the growing threat from AMR.

Genome-wide association, prediction and heritability in bacteria with application to Streptococcus pneumoniae

Whole-genome sequencing has facilitated genome-wide analyses of association, prediction and heritability in many organisms. However, such analyses in bacteria are still in their infancy, being limited by difficulties including genome plasticity and strong population structure. Here we propose a suite of methods including linear mixed models, elastic net and LD-score regression, adapted to bacterial traits using innovations such as frequency-based allele coding, both insertion/deletion and nucleotide testing and heritability partitioning. We compare and validate our methods against the current state-of-art using simulations, and analyse three phenotypes of the major human pathogen Streptococcus pneumoniae, including the first analyses of minimum inhibitory concentrations (MIC) for penicillin and ceftriaxone. We show that the MIC traits are highly heritable with high prediction accuracy, explained by many genetic associations under good population structure control. In ceftriaxone MIC, this is surprising because none of the isolates are resistant as per the inhibition zone criteria. We estimate that half of the heritability of penicillin MIC is explained by a known drug-resistance region, which also contributes a quarter of the ceftriaxone MIC heritability. For the within-host carriage duration phenotype, no associations were observed, but the moderate heritability and prediction accuracy indicate a moderately polygenic trait.

Whole genome sequences reveal the Xanthomonas perforans population is shaped by the tomato production system

Modern agricultural practices increase the potential for plant pathogen spread, while the advent of affordable whole genome sequencing enables in-depth studies of pathogen movement. Population genomic studies may decipher pathogen movement and population structure as a result of complex agricultural production systems. We used whole genome sequences of 281 Xanthomonas perforans strains collected within one tomato production season across Florida and southern Georgia fields to test for population genetic structure associated with tomato production system variables. We identified six clusters of X. perforans from core gene SNPs that corresponded with phylogenetic lineages. Using whole genome SNPs, we found genetic structure among farms, transplant facilities, cultivars, seed producers, grower operations, regions, and counties. Overall, grower operations that produced their own transplants were associated with genetically distinct and less diverse populations of strains compared to grower operations that received transplants from multiple sources. The degree of genetic differentiation among components of Florida's tomato production system varied between clusters, suggesting differential dispersal of the strains, such as through seed or contaminated transplants versus local movement within farms. Overall, we showed that the genetic variation of a bacterial plant pathogen is shaped by the structure of the plant production system.

The interplay between community and hospital Enterococcus faecium clones within health-care settings: a genomic analysis

BACKGROUND

The genomic relationships among Enterococcus faecium isolates are the subject of ongoing research that seeks to clarify the origins of observed lineages and the extent of horizontal gene transfer between them, and to robustly identify links between genotypes and phenotypes. E faecium is considered to form distinct groups-A and B-corresponding to isolates derived from patients who were hospitalised (A) and isolates from humans in the community (B). The additional separation of A into the so-called clades A1 and A2 remains an area of uncertainty. We aimed to investigate the relationships between A1 and non-A1 groups and explore the potential role of non-A1 isolates in shaping the population structure of hospital E faecium.

METHODS

We collected short-read sequence data from invited groups that had previously published E faecium genome data. This hospital-based isolate collection could be separated into three groups (or clades, A1, A2, and B) by augmenting the study genomes with published sequences derived from human samples representing the previously defined genomic clusters. We performed phylogenetic analyses, by constructing maximum-likelihood phylogenetic trees, and identified historical recombination events. We assessed the pan-genome, did resistome analysis, and examined the genomic data to identify mobile genetic elements. Each genome underwent chromosome painting by use of ChromoPainter within FineSTRUCTURE software to assess ancestry and identify hybrid groups. We further assessed highly admixed regions to infer recombination directionality.

FINDINGS

We assembled a collection of 1095 hospital E faecium sequences from 34 countries, further augmented by 33 published sequences. 997 (88%) of 1128 genomes clustered as A1, 92 (8%) as A2, and 39 (4%) as B. We showed that A1 probably emerged as a clone from within A2 and that, because of ongoing gene flow, hospital isolates currently identified as A2 represent a genetic continuum between A1 and community E faecium. This interchange of genetic material between isolates from different groups results in the emergence of hybrid genomes between clusters. Of the 1128 genomes, 49 (4%) hybrid genomes were identified: 33 previously labelled as A2 and 16 previously labelled as A1. These interactions were fuelled by a directional pattern of recombination mediated by mobile genetic elements. By contrast, the contribution of B group genetic material to A1 was limited to a few small regions of the genome and appeared to be driven by genomic sweep events.

INTERPRETATION

A2 and B isolates coming into the hospital form an important reservoir for ongoing A1 adaptation, suggesting that effective long-term control of the effect of E faecium could benefit from strategies to reduce these genomic interactions, such as a focus on reducing the acquisition of hospital A1 strains by patients entering the hospital.

FUNDING

Wellcome Trust.

A Phylogeographical Analysis of the Beetle Pest Species Callosobruchus chinensis (Linnaeus, 1758) in China

Callosobruchus chinensis (Coleoptera Bruchidae), is a pest of different varieties of legumes. In this paper, a phylogeographical analysis of C. chinensis was conducted to provide knowledge for the prevention and control of C. chinensis. A total of 224 concatenated mitochondrial sequences were obtained from 273 individuals. Suitable habitat shifts were predicted by the distribution modelling (SDM). Phylogeny, genetic structure and population demographic history were analyzed using multiple software. Finally, the least-cost path (LCP) method was used to identify possible dispersal corridors and genetic connectivity. The SDM results suggested that the distribution of C. chinensis experienced expansion and contraction with changing climate. Spatial distribution of mtDNA haplotypes showed there was partial continuity among different geographical populations of C. chinensis, except for the Hohhot (Inner Mongolia) population. Bayesian skyline plots showed that the population had a recent expansion during 0.0125 Ma and 0.025 Ma. The expansion and divergent events were traced back to Quaternary glaciations. The LCP method confirmed that there were no clear dispersal routes. Our findings indicated that climatic cycles of the Pleistocene glaciations, unsuitable climate and geographic isolation played important roles in the genetic differentiation of C. chinensis. Human activities weaken the genetic differentiation between populations. With the change in climate, the suitable areas of C. chinensis will disperse greatly in the future.

Frequencies and characteristics of genome-wide recombination in Streptococcus agalactiae, Streptococcus pyogenes, and Streptococcus suis

Streptococcus consists of ecologically diverse species, some of which are important pathogens of humans and animals. We sought to quantify and compare the frequencies and characteristics of within-species recombination in the pan-genomes of Streptococcus agalactiae, Streptococcus pyogenes and Streptococcus suis. We used 1081, 1813 and 1204 publicly available genome sequences of each species, respectively. Based on their core genomes, S. agalactiae had the highest relative rate of recombination to mutation (11.5743) compared to S. pyogenes (1.03) and S. suis (0.57). The proportion of the species pan-genome that have had a history of recombination was 12.85%, 24.18% and 20.50% of the pan-genomes of each species, respectively. The composition of recombining genes varied among the three species, and some of the most frequently recombining genes are implicated in adhesion, colonization, oxidative stress response and biofilm formation. For each species, a total of 22.75%, 29.28% and 18.75% of the recombining genes were associated with prophages. The cargo genes of integrative conjugative elements and integrative and mobilizable elements contained genes associated with antimicrobial resistance and virulence. Homologous recombination and mobilizable pan-genomes enable the creation of novel combinations of genes and sequence variants, and the potential for high-risk clones to emerge.

OUP accepted manuscript

The Zika virus (ZIKV) disease caused a public health emergency of international concern that started in February 2016. The overall number of ZIKV-related cases increased until November 2016, after which it declined sharply. While the evaluation of the potential risk and impact of future arbovirus epidemics remains challenging, intensified surveillance efforts along with a scale-up of ZIKV whole-genome sequencing provide an opportunity to understand the patterns of genetic diversity, evolution, and spread of ZIKV. However, a classification system that reflects the true extent of ZIKV genetic variation is lacking. Our objective was to characterize ZIKV genetic diversity and phylodynamics, identify genomic footprints of differentiation patterns, and propose a dynamic classification system that reflects its divergence levels. We analysed a curated dataset of 762 publicly available sequences spanning the full-length coding region of ZIKV from across its geographical span and collected between 1947 and 2021. The definition of genetic groups was based on comprehensive evolutionary dynamics analyses, which included recombination and phylogenetic analyses, within- and between-group pairwise genetic distances comparison, detection of selective pressure, and clustering analyses. Evidence for potential recombination events was detected in a few sequences. However, we argue that these events are likely due to sequencing errors as proposed in previous studies. There was evidence of strong purifying selection, widespread across the genome, as also detected for other arboviruses. A total of 50 sites showed evidence of positive selection, and for a few of these sites, there was amino acid (AA) differentiation between genetic clusters. Two main genetic clusters were defined, ZA and ZB, which correspond to the already characterized ‘African’ and ‘Asian’ genotypes, respectively. Within ZB, two subgroups, ZB.1 and ZB.2, represent the Asiatic and the American (and Oceania) lineages, respectively. ZB.1 is further subdivided into ZB.1.0 (a basal Malaysia sequence sampled in the 1960s and a recent Indian sequence), ZB.1.1 (South-Eastern Asia, Southern Asia, and Micronesia sequences), and ZB.1.2 (very similar sequences from the outbreak in Singapore). ZB.2 is subdivided into ZB.2.0 (basal American sequences and the sequences from French Polynesia, the putative origin of South America introduction), ZB.2.1 (Central America), and ZB.2.2 (Caribbean and North America). This classification system does not use geographical references and is flexible to accommodate potential future lineages. It will be a helpful tool for studies that involve analyses of ZIKV genomic variation and its association with pathogenicity and serve as a starting point for the public health surveillance and response to on-going and future epidemics and to outbreaks that lead to the emergence of new variants.

Phenotypic and genetic analyses of two Campylobacter fetus isolates from a patient with relapsed prosthetic valve endocarditis

Campylobacter fetus can cause intestinal and systemic disease in humans and are well established veterinary and economic pathogens. We report the complete genomic sequences of two C. fetus subsp. fetus (Cff) isolates recovered in 2017 (CITCf01) and 2018 (CITCf02) from a case of recurrent prosthetic valve endocarditis. Both were capable of growth aerobically. Their genomes were found to be highly conserved and syntenic with 99.97% average nucleotide identity (ANI) while differences in their respective sap loci defined the temporal separation of their genomes. Based on core genome phylogeny and ANI of 83 Cff genomes belonging to the previously described human-associated Cff lineage, CITCf01 and CITCf02 grouped in a clade of eleven sequence type (ST)3 Cff (including the Cff type strain NCTC 10842T). CITCf01 and CITCf02 were marked for their lack of unique genomic features when compared to isolates within the subspecies and the type strain in particular. We identified point mutations in oxidative stress response genes, among others, that may contribute to aerobiosis. We report a case of Cff causing relapsed prosthetic valve endocarditis and we highlight the sap island as a polymorphic site within the genetically stable ST3 lineage, central to pathogenicity.

A Severe Gastroenteritis Outbreak of Salmonella enterica Serovar Enteritidis Linked to Contaminated Egg Fried Rice, China, 2021

Salmonella contamination of eggs and egg shells has been identified as a public health problem worldwide. Here, we reported an outbreak of severe gastrointestinal symptoms caused by Salmonella enterica serovar Enteritidis (S. enteritidis) in China. We evaluated the outbreak by using epidemiological surveys, routine laboratory testing methods, and whole genome sequencing (WGS). This outbreak occurred in a canteen in Beijing, during March 9–11, 2021, 225 of the 324 diners who have eaten at the canteen showed gastrointestinal symptoms. The outbreak had characteristical epidemiological and clinical features. It caused a very high attack rate (69.4%) in a short incubation time. All patients developed diarrhea and high fever, accompanied by abdominal pain (62.3%), nausea (50.4%), and vomiting (62.7%). The average frequency of diarrhea was 12.4 times/day, and the highest frequency of diarrhea was as high as 50 times/day. The average fever temperature was 39.4°C, and the highest fever temperature was 42°C. Twenty strains of S. enteritidis were recovered, including 19 from the patients samples, and one from remained egg fried rice. Antibiotic susceptibility test showed that the 20 outbreak strains all had the same resistance pattern. PFGE results demonstrated that all 20 strains bore completely identical bands. Phylogenetic analysis based on WGS revealed that all 20 outbreak strains were tightly clustered together. So the pathogenic source of this food poisoning incident may was contaminated egg fried rice. Resistance gene analysis showed that the outbreak strains are all multi-drug resistant strains. Virulence gene analysis indicated that these outbreak strains carried a large number of virulence genes, including 2 types of Salmonella pathogenicity islands (SPI-1 and SPI-2). Other important virulence genes were also carried by the outbreak strains, such as pefABCD, rck and shdA. And the shdA gene was not in other strains located in the same evolutionary branch as the outbreak strain. We speculated that this is a significant reason for the serious symptoms of gastroenteritis in this outbreak. This outbreak caused by S. enteritidis suggested government should strengthen monitoring of the prevalence of outbreak clone strains, and take measures to mitigate the public health threat posed by contaminated eggs.

The global population structure and evolutionary history of the acquisition of major virulence factor-encoding genetic elements in Shiga toxin-producing Escherichia coli O121:H19

Shiga toxin (Stx)-producing Escherichia coli (STEC) are foodborne pathogens causing serious diseases, such as haemorrhagic colitis and haemolytic uraemic syndrome. Although O157:H7 STEC strains have been the most prevalent, incidences of STEC infections by several other serotypes have recently increased. O121:H19 STEC is one of these major non-O157 STECs, but systematic whole genome sequence (WGS) analyses have not yet been conducted on this STEC. Here, we performed a global WGS analysis of 638 O121:H19 strains, including 143 sequenced in this study, and a detailed comparison of 11 complete genomes, including four obtained in this study. By serotype-wide WGS analysis, we found that O121:H19 strains were divided into four lineages, including major and second major lineages (named L1 and L3, respectively), and that the locus of enterocyte effacement (LEE) encoding a type III secretion system (T3SS) was acquired by the common ancestor of O121:H19. Analyses of 11 complete genomes belonging to L1 or L3 revealed remarkable interlineage differences in the prophage pool and prophage-encoded T3SS effector repertoire, independent acquisition of virulence plasmids by the two lineages, and high conservation in the prophage repertoire, including that for Stx2a phages in lineage L1. Further sequence determination of complete Stx2a phage genomes of 49 strains confirmed that Stx2a phages in lineage L1 are highly conserved short-tailed phages, while those in lineage L3 are long-tailed lambda-like phages with notable genomic diversity, suggesting that an Stx2a phage was acquired by the common ancestor of L1 and has been stably maintained. Consistent with these genomic features of Stx2a phages, most lineage L1 strains produced much higher levels of Stx2a than lineage L3 strains. Altogether, this study provides a global phylogenetic overview of O121:H19 STEC and shows the interlineage genomic differences and the highly conserved genomic features of the major lineage within this serotype of STEC.

Evolutionary Divergence of the Novel Staphylococcal Species Staphylococcus argenteus

Currently, invasive infections caused by Staphylococcus argenteus, which is a recently named staphylococcal species, are increasingly reported worldwide. However, only a few genomic studies of S. argenteus have offered comprehensive information regarding its genetic diversity, epidemiological characteristics, antimicrobial resistance genes (ARGs), virulence genes and other profiles. Here, we describe a comparative genomic analysis by population structure, pangenome, panmobilome, region-specific accessory genes confer an adaptive advantage in 153 S. argenteus strains which comprised 24 strains sequenced in this study and 129 strains whose genome sequences were available from GenBank. As a result, the population of S. argenteus comprised seven genetically distinct clades, including two major clades (C1 and C2), with distinct isolation source patterns. Pangenome analysis revealed that S. argenteus has an open pangenome composed of 7,319 genes and a core genome composed of 1,508 genes. We further determined the distributions of 75 virulence factors (VFs) and 30 known ARGs and identified at least four types of plasmids and 93 complete or partial putative prophages. It indicate that S. argenteus may show a similar level of pathogenicity to that of S. aureus. This study also provides insights into the evolutionary divergence of this pathogen, indicating that the geographical distribution was a potential driving force behind the evolutionary divergence of S. argenteus. The preferential horizontal acquisition of particular elements, such as staphylococcal cassette chromosome mec elements and plasmids, was observed in specific regions, revealing potential gene exchange between S. argenteus strains and local S. aureus strains. Moreover, multiple specific genes related to environmental adaptation were identified in strains isolated from East Asia. However, these findings may help promote our understanding of the evolutionary divergence of this bacterium at a high genetic resolution by providing insights into the epidemiology of S. argenteus and may help combat its spread.

Molecular evidence and ecological niche modeling reveal an extensive hybrid zone among three Bursera species (section Bullockia)

The genus Bursera, includes ~100 shrub and trees species in tropical dry forests with its center of diversification and endemism in Mexico. Morphologically intermediate individuals have commonly been observed in Mexican Bursera in areas where closely related species coexist. These individuals are assumed to result from interspecific hybridization, but no molecular evidence has supported their hybrid origins. This study aimed to investigate the existence of interspecific hybridization among three Mexican Bursera species (Bullockia section: B. cuneata, B. palmeri and B. bipinnata) from nine populations based on DNA sequences (three nuclear and four chloroplast regions) and ecological niche modeling for three past and two future scenario projections. Results from the only two polymorphic nuclear regions (PEPC, ETS) supported the hybrid origin of morphologically intermediate individuals and revealed that B. cuneata and B. bipinnata are the parental species that are genetically closer to the putative hybrids. Ecological niche modeling accurately predicted the occurrence of putative hybrid populations and showed a potential hybrid zone extending in a larger area (74,000 km2) than previously thought. Paleo-reconstructions showed a potential hybrid zone existing from the Last Glacial Maximum (~ 21 kya) that has increased since the late Holocene to the present. Future ecological niche projections show an increment of suitability of the potential hybrid zone for 2050 and 2070 relative to the present. Hybrid zone changes responded mostly to an increase in elevational ranges. Our study provides the first insight of an extensive hybrid zone among three Mexican Bursera species based on molecular data and ecological niche modeling.

Genetic diversity and transmission patterns of Burkholderia pseudomallei on Hainan island, China, revealed by a population genomics analysis

Burkholderia pseudomallei is a Gram-negative soil-dwelling bacillus that causes melioidosis, a frequently fatal infectious disease, in tropical and subtropical regions. Previous studies have identified the overall genetic and evolutionary characteristics of B. pseudomallei on a global scale, including its origin and transmission routes. However, beyond its known hyperendemicity foci in northern Australia and Southeast Asia, the distribution and genetic characteristics of B. pseudomallei in most tropical regions remain poorly understood, including in southern China. Here, we sequenced the genomes of 122 B. pseudomallei strains collected from Hainan, an island in southern China, in 2002–2018, to investigate the population structure, relationships with global strains, local epidemiology, and virulence and antimicrobial-resistance factors. A phylogenetic analysis and hierarchical clustering divided the Hainan strains into nine phylogenic groups (PGs), 80 % of which were concentrated within five major groups (group 1: corresponding to minor sequence types [STs], 12.3 %; group 3: ST46 and ST50, 31.1 %; group 9: ST58, 13.1 %; group 11: ST55, 8.2 %; group 15: mainly ST658, 15.6%). A phylogenetic analysis that included global strains suggested that B. pseudomallei in Hainan originated from Southeast Asian countries, transmitted in multiple historical importation events. We also identified several mutual transmission events between Hainan and Southeast Asian countries in recent years, including three importation events from Thailand and Singapore to Hainan and three exportation events from Hainan to Singapore, Malaysia, and Taiwan island. A statistical analysis of the temporal distribution showed that the Hainan strains of groups 3, 9, and 15 have dominated the disease epidemic locally in the last 5 years. The spatial distribution of the Hainan strains demonstrated that some PGs are distributed in different cities on Hainan island, and by combining phylogenic and geographic distribution information, we detected 21 between-city transmission events, indicating its frequent local transmission. The detection of virulence factor genes showed that 56 % of the Hainan strains in group 1 encode a B. pseudomallei-specific adherence factor, boaB, confirming the specific pathogenic characteristics of the Hainan strains in group 1. An analysis of the antimicrobial-resistance potential of B. pseudomallei showed that various kinds of alterations were identified in clinically relevant antibiotic resistance factors, such as AmrR, PenA and PBP3, etc. Our results clarify the population structure, local epidemiology, and pathogenic characteristics of B. pseudomallei in Hainan, providing further insight into its regional and global transmission networks and improving our knowledge of its global phylogeography.

The antimicrobial resistance landscape of Neisseria gonorrhoeae in New Zealand from November 2018 to March 2019 and the role of sexual orientation in transmission

The increasing use of culture independent diagnostic testing for the diagnosis of Neisseria gonorrhoeae infection has led to gaps in surveillance of antimicrobial resistance (AMR) rates due to limited availability of cultures. Our study reports the findings of a second national survey of N. gonorrhoeae in New Zealand, utilizing whole-genome sequencing (WGS) to study the population structure, prevalence of AMR, epidemiology and transmission of gonorrhoea isolates. We analysed 314 isolates and found a strong correlation between carriage of acquired resistance genes or chromosomal point mutations and phenotypic susceptibility testing results. Overall, the New Zealand rates of azithromycin resistance and decreased susceptibility to ceftriaxone remain lower than in most countries, which are part of the World Health Organization (WHO) Global Gonococcal Antimicrobial Surveillance Programme (GASP). The phylogeny provides evidence of a diverse population significantly associated with sexual behaviour groups. Transmission clustering with a ten single nucleotide polymorphism (SNP) cut-off identified 49 clusters, of which ten were solely associated with men who have sex with men (MSM), whereas remaining clusters included heterosexual patients, as well as MSM, suggesting that bridging of sexual networks is occurring. Utilizing pairwise SNP differences between isolates of the same sequence types we determined genetic variation for the three typing schemes used in this study [Multi locus sequence typing (MLST), multi-antigen sequence typing (NG-MAST), and sequence typing for antimicrobial resistance (NG-STAR)]. A median of 0.0 to 52.5 pairwise SNP differences within a single NG-STAR sequence type underlines previous findings of the superiority of the NG-STAR typing scheme in terms of genomic inherency. With our analysis incorporating epidemiological and genomic data, we were able to show a comprehensive overview of the N. gonorrhoeae population circulating in New Zealand, focussing on AMR and transmission within sexual networks. Regular surveillance studies to understand the origin, evolution and spread of AMR for gonorrhoea remain necessary to make informed decisions about public health guidelines, as the internationally rising rates of ceftriaxone and azithromycin resistance have already led to adaptation of current treatment guidelines in the UK and the USA, highlighting the importance of regular surveillance in individual countries.

Genomic Characterization of Endemic and Ecdemic Non-typhoidal Salmonella enterica Lineages Circulating Among Animals and Animal Products in South Africa

In Africa, the burden of illness caused by non-typhoidal Salmonella enterica is disproportionally high; however, whole-genome sequencing (WGS) efforts are overwhelmingly concentrated in world regions with lower burdens. While WGS is being increasingly employed in South Africa to characterize Salmonella enterica, the bulk of these efforts have centered on characterizing human clinical strains. Thus, very little is known about lineages circulating among animals in the country on a genomic scale. Here, we used WGS to characterize 63 Salmonella enterica strains isolated from livestock, companion animals, wildlife, and animal products in South Africa over a 60-year period. Genomes were assigned to serotypes Dublin, Hadar, Enteritidis, and Typhimurium (n = 18, 8, 13, and 24 strains, respectively) and sequence types (STs) ST10 (all S. Dublin), ST33 (all S. Hadar), ST11/ST366 (n = 12 and 1 S. Enteritidis, respectively), and ST19/ST34 (n = 23 and 1 S. Typhimurium, respectively; via seven-gene multi-locus sequence typing). Within-ST phylogenies were constructed using genomes sequenced in this study, plus publicly available genomes representative of each ST's (i) global (n = 2,802 and 1,569 S. Dublin and Hadar genomes, respectively) and (ii) African (n = 716 and 343 S. Enteritidis and Typhimurium genomes, respectively) population. For S. Dublin ST10, a largely antimicrobial-susceptible, endemic lineage circulating among humans, animals, and food in South Africa was identified, as well as a lineage that was likely recently introduced from the United States. For S. Hadar ST33, multiple South African lineages harboring streptomycin and tetracycline resistance-conferring genes were identified. African S. Enteritidis ST11 could be primarily partitioned into one largely antimicrobial-susceptible and one largely multidrug-resistant (MDR) clade, with South African isolates confined to the largely antimicrobial-susceptible clade. S. Typhimurium ST19/ST34 strains sequenced here were distributed across the African S. Typhimurium ST19/ST34 phylogeny, representing a diverse range of lineages, including numerous MDR lineages. Overall, this study provides critical insights into endemic and ecdemic non-typhoidal Salmonella enterica lineages circulating among animals, foods, and humans in South Africa and showcases the utility of WGS in characterizing animal-associated strains from a world region with a high salmonellosis burden.

High-resolution sweep metagenomics using fast probabilistic inference

Determining the composition of bacterial communities beyond the level of a genus or species is challenging because of the considerable overlap between genomes representing close relatives. Here, we present the mSWEEP pipeline for identifying and estimating the relative sequence abundances of bacterial lineages from plate sweeps of enrichment cultures. mSWEEP leverages biologically grouped sequence assembly databases, applying probabilistic modelling, and provides controls for false positive results. Using sequencing data from major pathogens, we demonstrate significant improvements in lineage quantification and detection accuracy. Our pipeline facilitates investigating cultures comprising mixtures of bacteria, and opens up a new field of plate sweep metagenomics.

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

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

Population Genetic Structure of the Invasive Spotted Alfalfa Aphid Therioaphis trifolii (Hemiptera: Aphididae) in China Inferred From Complete Mitochondrial Genomes

Biological invasions represent a natural rapid evolutionary process in which invasive species may present a major threat to biodiversity and ecosystem integrity. Analyzing the genetic structure and demographic history of invaded populations is critical for the effective management of invasive species. The spotted alfalfa aphid (SAA) Therioaphis trifolii is indigenous in the Mediterranean region of Europe and Africa and has invaded China, causing severe damages to the alfalfa industry. However, little is known about its genetic structure and invasion history. In this study, we obtained 167 complete mitochondrial genome sequences from 23 SAA populations across China based on high-throughput sequencing and performed population genetic and phylogenomic analyses. High haplotype diversity and low nucleotide diversity were found in SAA populations in China with distinct genetic structures, i.e., all populations diverged into three phylogenetic lineages. Demographic history analyses showed a recent expansion of the SAA population, consistent with the recent invasion history. Our study indicated that SAA may have invaded through multiple introduction events during commercial trades of alfalfa, although this needs further validation by nuclear markers.

A role of asynchrony of seasons in explaining genetic differentiation in a Neotropical toad

The process of diversification can be studied at the phylogeographic level by attempting to identify the environmental features that promote and maintain population divergence. Here we investigate diversification in Rhinella granulosa, a Neotropical toad from northeastern Brazil, by testing a range of hypotheses that encompass different putative mechanisms reducing gene flow among populations. We sequenced single nucleotide polymorphisms and examined individual predictions related to the role of geographic barriers (rivers), ecological gradients, historical habitat stability, and spatial variation in climate seasonality, also known as the asynchrony of seasons hypothesis. This hypothesis postulates that temporal asynchrony of wet and dry seasons over short distances causes parapatric populations to become isolated by time. After determining genetic structure, inferring past distributions, ranking demographic models, and estimating the power of monthly climatic variables, our results identified two populations that are not associated with geographic barriers, biome gradients, or historical refugia. Instead, they are predicted by spatial variation in monthly rainfall and minimum temperature, consistent with the asynchrony of seasons hypothesis, supported also by our comparative framework using multiple matrix regression and linear mixed effects modeling. Due to the toad's life history, climate likely mediates gene flow directly, with genetic differentiation being provoked by neutral mechanisms related to climate driven population isolation, and/or by natural selection against migrants from populations with different breeding times. The asynchrony of seasons hypothesis is seldom considered in phylogeographic studies, but our results indicate that it should be tested in systems where breeding is tightly coupled with climate.

Global Expansion of Linezolid-Resistant Coagulase-Negative Staphylococci

Coagulase-negative staphylococci (CoNS) for a long time were considered avirulent constituents of the human and warm-blooded animal microbiota. However, at present, S. epidermidis, S. haemolyticus, and S. hominis are recognized as opportunistic pathogens. Although linezolid is not registered for the treatment of CoNS infections, it is widely used off-label, promoting emergence of resistance. Bioinformatic analysis based on maximum-likelihood phylogeny and Bayesian clustering of the CoNS genomes obtained in the current study and downloaded from public databases revealed the existence of international linezolid-resistant lineages, each of which probably had a common predecessor. Linezolid-resistant S. epidermidis sequence-type (ST) 2 from Russia, France, and Germany formed a compact group of closely related genomes with a median pairwise single nucleotide polymorphism (SNP) difference of fewer than 53 SNPs, and a common ancestor of this lineage appeared in 1998 (1986-2006) before introduction of linezolid in practice. Another compact group of linezolid-resistant S. epidermidis was represented by ST22 isolates from France and Russia with a median pairwise SNP difference of 40; a common ancestor of this lineage appeared in 2011 (2008-2013). Linezolid-resistant S. hominis ST2 from Russia, Germany, and Brazil also formed a group with a high-level genome identity with median 25.5 core-SNP differences; the appearance of the common progenitor dates to 2003 (1996-2012). Linezolid-resistant S. hominis isolates from Russia demonstrated associated resistance to teicoplanin. Analysis of a midpoint-rooted phylogenetic tree of the group confirmed the genetic proximity of Russian and German isolates; Brazilian isolates were phylogenetically distant. repUS5-like plasmids harboring cfr were detected in S. hominis and S. haemolyticus.

Weak spatial-genetic structure in a native invasive, the southern pine beetle (Dendroctonus frontalis), across the eastern United States

The southern pine beetle, Dendroctonus frontalis, is a native pest of pine trees that has recently expanded its range into the northeastern United States. Understanding its colonization, dispersal, and connectivity will be critical for mitigating negative economic and ecological impacts in the newly invaded areas. Characterization of spatial-genetic structure can contribute to this; however, previous studies have reached different conclusions about regional population genetic structure, with one study reporting a weak east-west pattern, and the most recent reporting an absence of structure. Here we systematically assessed several explanations for the absence of spatial-genetic structure. To do this, we developed nine new microsatellite markers and combined them with an existing 24-locus data matrix for the same individuals. We then reanalyzed this full dataset alongside datasets in which certain loci were omitted with the goal of creating more favorable signal to noise ratios. We also partitioned the data based on the sex of D. frontalis individuals, and then employed a broad suite of genotypic clustering and isolation-by-distance (IBD) analyses. We found that neither inadequate information content in the molecular marker set, nor unfavorable signal-to-noise ratio, nor insensitivity of the analytical approaches could explain the absence of structure. Regardless of dataset composition, there was little evidence for clusters (i.e., distinct geo-genetic groups) or clines (i.e., gradients of increasing allele frequency differences over larger geographic distances), with one exception: significant IBD was repeatedly detected using an individual-based measure of relatedness whenever datasets included males (but not for female-only datasets). This is strongly indicative of broad-scale female-biased dispersal, which has not previously been reported for D. frontalis, in part owing to logistical limitations of direct approaches (e.g., capture-mark-recapture). Weak spatial-genetic structure suggests long-distance connectivity and that gene flow is high, but additional research is needed to understand range expansion and outbreak dynamics in this species using alternate approaches.

Streptococcus pneumoniae serotypes that frequently colonise the human nasopharynx are common recipients of penicillin-binding protein gene fragments from Streptococcus mitis

Streptococcus pneumoniae is an important global pathogen that causes bacterial pneumonia, sepsis and meningitis. Beta-lactam antibiotics are the first-line treatment for pneumococcal disease, however, their effectiveness is hampered by beta-lactam resistance facilitated by horizontal genetic transfer (HGT) with closely related species. Although interspecies HGT is known to occur among the species of the genus Streptococcus, the rates and effects of HGT between Streptococcus pneumoniae and its close relatives involving the penicillin binding protein (pbp) genes remain poorly understood. Here we applied the fastGEAR tool to investigate interspecies HGT in pbp genes using a global collection of whole-genome sequences of Streptococcus mitis, Streptococcus oralis and S. pneumoniae. With these data, we established that pneumococcal serotypes 6A, 13, 14, 16F, 19A, 19F, 23F and 35B were the highest-ranking serotypes with acquired pbp fragments. S. mitis was a more frequent pneumococcal donor of pbp fragments and a source of higher pbp nucleotide diversity when compared with S. oralis. Pneumococci that acquired pbp fragments were associated with a higher minimum inhibitory concentration (MIC) for penicillin compared with pneumococci without acquired fragments. Together these data indicate that S. mitis contributes to reduced β-lactam susceptibility among commonly carried pneumococcal serotypes that are associated with long carriage duration and high recombination frequencies. As pneumococcal vaccine programmes mature, placing increasing pressure on the pneumococcal population structure, it will be important to monitor the influence of antimicrobial resistance HGT from commensal streptococci such as S. mitis.

Occurrence of Colibacillosis in Broilers and Its Relationship With Avian Pathogenic Escherichia coli (APEC) Population Structure and Molecular Characteristics

Avian pathogenic Escherichia coli (APEC) causes colibacillosis, the disease with the highest economic loss for the broiler industry. However, studies focusing on the prevalence and population structure of APEC in the broiler production pyramid are scarce. Here, we used genotyping and serotyping data to elucidate the APEC population structure and its changes in different broiler production stages along with whole-genome sequencing (WGS) in a subset of APEC isolates to determine transmission patterns amongst dominant APEC sequence types (STs) and characterize them in detail. Comparison of genotypes encountered in both APEC and avian fecal E. coli (AFEC) provided further insights. Overall, APEC-related mortality, as the proportion of the total sampled mortality in the broiler production, was high (35%), while phylogroup C and serogroup O78 were predominant amongst APEC isolates. We found a low (34.0%) and high (53.3%) incidence of colibacillosis in chicks and end-cycle broilers, respectively, which may be related to a shift in APEC genotypes, suggesting a trend from commensalism to pathogenicity across different broiler production stages. Despite considerable APEC genotypic diversity, there was substantial genotype overlap (40.9%, overall) over the production stages and convergence of STs to the four clusters. Within these clusters, WGS data provided evidence of clonal transmission events and revealed an enriched virulence and resistance APEC repertoire. More specifically, sequenced APEC were assigned to defined pathotypes based on their virulence gene content while the majority (86%) was genotypically multi-drug resistant. Interestingly, WGS-based phylogeny showed that a subset of APEC, which are cephalosporin-resistant, may originate directly from cephalosporin-resistant AFEC. Finally, exploration of the APEC plasmidome indicated that the small fraction of the APEC virulome carried by IncF plasmids is pivotal for the manifestation of the APEC pathotype; thus, plasmid exchange can promote pathogenicity in strains that are at the edge of the commensal and pathogenic states.

Large-scale genomic analysis of antimicrobial resistance in the zoonotic pathogen Streptococcus suis

BACKGROUND

Antimicrobial resistance (AMR) is among the gravest threats to human health and food security worldwide. The use of antimicrobials in livestock production can lead to emergence of AMR, which can have direct effects on humans through spread of zoonotic disease. Pigs pose a particular risk as they are a source of zoonotic diseases and receive more antimicrobials than most other livestock. Here we use a large-scale genomic approach to characterise AMR in Streptococcus suis, a commensal found in most pigs, but which can also cause serious disease in both pigs and humans.

RESULTS

We obtained replicated measures of Minimum Inhibitory Concentration (MIC) for 16 antibiotics, across a panel of 678 isolates, from the major pig-producing regions of the world. For several drugs, there was no natural separation into 'resistant' and 'susceptible', highlighting the need to treat MIC as a quantitative trait. We found differences in MICs between countries, consistent with their patterns of antimicrobial usage. AMR levels were high even for drugs not used to treat S. suis, with many multidrug-resistant isolates. Similar levels of resistance were found in pigs and humans from regions associated with zoonotic transmission. We next used whole genome sequences for each isolate to identify 43 candidate resistance determinants, 22 of which were novel in S. suis. The presence of these determinants explained most of the variation in MIC. But there were also interesting complications, including epistatic interactions, where known resistance alleles had no effect in some genetic backgrounds. Beta-lactam resistance involved many core genome variants of small effect, appearing in a characteristic order.

CONCLUSIONS

We present a large dataset allowing the analysis of the multiple contributing factors to AMR in S. suis. The high levels of AMR in S. suis that we observe are reflected by antibiotic usage patterns but our results confirm the potential for genomic data to aid in the fight against AMR.

Genomic Approaches to Plant-Pathogen Epidemiology and Diagnostics

Diseases have a significant cost to agriculture. Findings from analyses of whole-genome sequences show great promise for informing strategies to mitigate risks from diseases caused by phytopathogens. Genomic approaches can be used to dramatically shorten response times to outbreaks and inform disease management in novel ways. However, the use of these approaches requires expertise in working with big, complex data sets and an understanding of their pitfalls and limitations to infer well-supported conclusions. We suggest using an evolutionary framework to guide the use of genomic approaches in epidemiology and diagnostics of plant pathogens. We also describe steps that are necessary for realizing these as standard approaches in disease surveillance.

Kinship networks of seed exchange shape spatial patterns of plant virus diversity

By structuring farmers’ informal networks of seed exchange, kinship systems play a key role in the dynamics of crop genetic diversity in smallholder farming systems. However, because many crop diseases are propagated through infected germplasm, local seed systems can also facilitate the dissemination of seedborne pathogens. Here, we investigate how the interplay of kinship systems and local networks of germplasm exchange influences the metapopulation dynamics of viruses responsible for the cassava mosaic disease (CMD), a major threat to food security in Africa. Combining anthropological, genetic and plant epidemiological data, we analyzed the genetic structure of local populations of the African cassava mosaic virus (ACMV), one of the main causal agents of CMD. Results reveal contrasted patterns of viral diversity in patrilineal and matrilineal communities, consistent with local modes of seed exchange. Our results demonstrate that plant virus ecosystems have also a cultural component and that social factors that shape regional seed exchange networks influence the genetic structure of plant virus populations.

This study combines ethnobotanical and epidemiological data to understand how social networks of seed exchange influence the genetic structure of the African cassava mosaic virus in Gabon. Results reveal contrasted patterns of viral diversity in patrilineal and matrilineal communities, consistent with cultural differences in modes of seed exchange.

The global emergence of a novel Streptococcus suis clade associated with human infections

Streptococcus suis, a ubiquitous bacterial colonizer in pigs, has recently extended host range to humans, leading to a global surge of deadly human infections and three large outbreaks since 1998. To better understand the mechanisms for the emergence of cross-species transmission and virulence in human, we have sequenced 366 S. suis human and pig isolates from 2005 to 2016 and performed a large-scale phylogenomic analysis on 1,634 isolates from 14 countries over 36 years. We show the formation of a novel human-associated clade (HAC) diversified from swine S. suis isolates. Phylogeographic analysis identified Europe as the origin of HAC, coinciding with the exportation of European swine breeds between 1960s and 1970s. HAC is composed of three sub-lineages and contains several healthy-pig isolates that display high virulence in experimental infections, suggesting healthy-pig carriers as a potential source for human infection. New HAC-specific genes are identified as promising markers for pathogen detection and surveillance. Our discovery of a human-associated S. suis clade provides insights into the evolution of this emerging human pathogen and extend our understanding of S. suis epidemics worldwide.

Identification of a Recently Dominant Sublineage in Salmonella 4,[5],12:i:- Sequence Type 34 Isolated From Food Animals in Japan

Salmonella enterica subsp. enterica serovar Typhimurium sequence type 34 (ST34) and its monophasic variant (Salmonella 4,[5],12:i:-) are among the most frequently isolated clones from both humans and animals worldwide. Our previous study demonstrated that Salmonella Typhimurium/4,[5],12:i:- strains isolated in Japan could be classified into nine clades and that clade 9 consisted of ST34 strains. In Japan, ST34/clade 9 was first found in the 1990s and has become predominant among food animals in recent years. In the present study, we analyzed the whole genome-based phylogenetic relationships and temporal information of 214 Salmonella Typhimurium/4,[5],12:i:- ST34/clade 9 strains isolated from 1998 to 2017 in Japan. The 214 strains were classified into two sublineages: the newly identified clade 9–2 diverged from clade 9 in the early 2000s and has predominated in recent years. Clonally expanding subclades in clades 9–1 or 9–2 lacked Gifsy-1 or HP1 prophages, respectively, and some strains in these subclades acquired plasmids encoding antimicrobial resistance genes. Additional genome reduction around the fljB gene encoding the phase 2-H antigen was generated by an IS26-mediated deletion adjacent to the transposon in clade 9–2. Although most of the clade 9 strains were isolated from cattle in Japan, the clonally expanding subclades in clade 9–2 (i.e., all and 24% strains of subclades 9–2a and 9–2b, respectively) were isolated from swine. The spread of clade 9 in recent years among food animals in Japan was responsible for the emergence of multiple host-adapted sublineages involving the clonally expanding subclades generated by mobile genetic element-mediated microevolution.

Molecular epidemiology of carbapenemase-producing Escherichia coli from duck farms in south-east coastal China

OBJECTIVES

To determine the dissemination and molecular characteristics of NDM-producing Escherichia coli strains from duck farms in south-east coastal China and their threats to human health.

METHODS

A total of 232 NDM-producing E. coli were recovered from 1505 samples collected from 25 duck farms and their surrounding environments in five provinces in China. Resistance genes were confirmed using PCR. Genomic characteristics of the carbapenemase-producing isolates were determined by WGS and bioinformatic analysis.

RESULTS

The rate of NDM-positive E. coli detected in samples from the five provinces ranged from 3.7% to 28.5%. There was substantial variation in the prevalence of NDM-positive E. coli from different duck farms in each province studied. Three variants (blaNDM-1, blaNDM-4 and blaNDM-5) were found in 232 NDM-positive E. coli; blaNDM-5 (94.8%, 220/232) was the most prevalent. WGS analysis indicated that ST746, ST48, ST1011 and ST167 E. coli isolates were prevalent in the current study and poultry was likely the primary reservoir for NDM-positive ST746 and ST48 E. coli in China. Phylogenomic analysis showed that NDM-positive E. coli isolates from ducks were closely related to those of human origin. In addition, WGS analysis further revealed that blaNDM co-existed with other antibiotic resistance genes, conferring resistance to nine classes of antimicrobials.

CONCLUSIONS

This study revealed that ducks farm in China are an important reservoir for NDM-positive E. coli and STs of the isolates showed obvious distinctive diversities in geographical distribution. The distribution and spread of NDM-positive E. coli in duck farms poses a threat to public health.

Molecular Epidemiology of Fosfomycin Resistant E. coli from a Pigeon Farm in China

We determined the prevalence and molecular characteristics of fosfomycin-resistant Escherichia coli from a domestic pigeon farm. A total of 79 samples collected from pigeons and their surrounding environments were screened for the presence of fosfomycin resistant isolates and these included 49 E. coli isolates that displayed high-level resistance (MIC ≥ 256 mg L⁻¹) and carried the fosA3 gene on plasmids with sizes ranging from 80 to 370 kb. MLST analysis of these fosA3-positive E. coli isolates indicated the presence of nine sequence types (ST6856, ST8804, ST457, ST746, ST533, ST165, ST2614, ST362 and ST8805) of which ST6856 was the most prevalent (24.5%, 12/49). PFGE combined with genomic context comparative analyses indicated that the fosA3 gene was spread by horizontal transfer as well as via clonal transmission between E. coli in the pigeon farm, and IS26 played an important role in fosA3 transmission. The high prevalence of fosA3 in the pigeon farm and the high similarity of the fosA3 genomic environment between E. coli isolates from humans and pigeons indicated that the pigeon farm served as a potential reservoir for human infections. The pigeon farm was found to be an important reservoir for the fosA3 gene and this should be further monitored.

A Whole-Genome-Based Gene-by-Gene Typing System for Standardized High-Resolution Strain Typing of Bacillus anthracis

Whole-genome sequencing (WGS) has been established for bacterial subtyping and is regularly used to study pathogen transmission, to investigate outbreaks, and to perform routine surveillance. Core-genome multilocus sequence typing (cgMLST) is a bacterial subtyping method that uses WGS data to provide a high-resolution strain characterization. This study aimed at developing a novel cgMLST scheme for Bacillus anthracis, a notorious pathogen that causes anthrax in livestock and humans worldwide. The scheme comprises 3,803 genes that were conserved in 57 B. anthracis genomes spanning the whole phylogeny. The scheme has been evaluated and applied to 584 genomes from 50 countries. On average, 99.5% of the cgMLST targets were detected. The cgMLST results confirmed the classical canonical single-nucleotide-polymorphism (SNP) grouping of B. anthracis into major clades and subclades. Genetic distances calculated based on cgMLST were comparable to distances from whole-genome-based SNP analysis with similar phylogenetic topology and comparable discriminatory power. Additionally, the application of the cgMLST scheme to anthrax outbreaks from Germany and Italy led to a definition of a cutoff threshold of five allele differences to trace epidemiologically linked strains for cluster typing and transmission analysis. Finally, the association of two clusters of B. anthracis with human cases of injectional anthrax in four European countries was confirmed using cgMLST. In summary, this study presents a novel cgMLST scheme that provides high-resolution strain genotyping for B. anthracis. This scheme can be used in parallel with SNP typing methods to facilitate rapid and harmonized interlaboratory comparisons, essential for global surveillance and outbreak analysis. The scheme is publicly available for application by users, including those with little bioinformatics knowledge.

Whole-genome analysis uncovers loss of blaZ associated with carriage isolates belonging to methicillin-resistant Staphylococcus aureus (MRSA) clone ST5-VI in Cape Verde

•

One of the first whole genome analyses of Staphylococcus aureus carriage isolates in an African country.

•

Genome data allowed to place S. aureus isolates from Cape Verde in a phylogenetic context.

•

Loss of blaZ-carrying plasmids and transposons is not rare, which is also evident in other international MRSA isolates.

•

Steady increase in antimicrobial drug resistance in Cape Verde.

•

Data provide genomic information for the design of intervention measures to decrease antimicrobial resistance.

Objectives

Surveillance studies for Staphylococcus aureus carriage are a primary tool to survey the prevalence of methicillin-resistant S. aureus (MRSA) in the general population, patients and healthcare workers. We have previously reported S. aureus carriage in various African countries, including Cape Verde.

Methods

Whole-genome sequences of 106 S. aureus isolates from Cape Verde were determined.

Results

Staphylococcus aureus carriage isolates in Cape Verde show high genetic variability, with the detection of 27 sequence types (STs) and three primary genetic clusters associated with ST152, ST15 and ST5. One transmission event with less than eight core-genome single nucleotide polymorphisms (cgSNP) differences was detected among the ST5-VI MRSA lineage. Genetic analysis confirmed the phenotypic resistance and allowed the identification of six independent events of plasmid or transposon loss associated with the deletion of blaZ in nine isolates. In the four ST5 MRSA isolates, loss of the blaZ plasmid coincided with the acquisition of SCCmec type VI and an unusual penicillin phenotype with a minimum inhibitory concentration (MIC) at the breakpoint, indicating an adaptation trend in this endemic lineage. Similar events of blaZ plasmid loss, with concomitant acquisition SCCmec elements, were detected among ST5 isolates from different geographical origins.

Conclusion

Overall, the genome data allowed to place isolates in a phylogenetic context and to identify different blaZ gene deletions associated with plasmid or transposon loss. Genomic analysis unveiled adaptation and evolution trends, namely among emerging MRSA lineages in the country, which deserve additional consideration in the design of future infection control protocols.

ProkEvo: an automated, reproducible, and scalable framework for high-throughput bacterial population genomics analyses

Whole Genome Sequence (WGS) data from bacterial species is used for a variety of applications ranging from basic microbiological research, diagnostics, and epidemiological surveillance. The availability of WGS data from hundreds of thousands of individual isolates of individual microbial species poses a tremendous opportunity for discovery and hypothesis-generating research into ecology and evolution of these microorganisms. Flexibility, scalability, and user-friendliness of existing pipelines for population-scale inquiry, however, limit applications of systematic, population-scale approaches. Here, we present ProkEvo, an automated, scalable, reproducible, and open-source framework for bacterial population genomics analyses using WGS data. ProkEvo was specifically developed to achieve the following goals: (1) Automation and scaling of complex combinations of computational analyses for many thousands of bacterial genomes from inputs of raw Illumina paired-end sequence reads; (2) Use of workflow management systems (WMS) such as Pegasus WMS to ensure reproducibility, scalability, modularity, fault-tolerance, and robust file management throughout the process; (3) Use of high-performance and high-throughput computational platforms; (4) Generation of hierarchical-based population structure analysis based on combinations of multi-locus and Bayesian statistical approaches for classification for ecological and epidemiological inquiries; (5) Association of antimicrobial resistance (AMR) genes, putative virulence factors, and plasmids from curated databases with the hierarchically-related genotypic classifications; and (6) Production of pan-genome annotations and data compilation that can be utilized for downstream analysis such as identification of population-specific genomic signatures. The scalability of ProkEvo was measured with two datasets comprising significantly different numbers of input genomes (one with ~2,400 genomes, and the second with ~23,000 genomes). Depending on the dataset and the computational platform used, the running time of ProkEvo varied from ~3-26 days. ProkEvo can be used with virtually any bacterial species, and the Pegasus WMS uniquely facilitates addition or removal of programs from the workflow or modification of options within them. To demonstrate versatility of the ProkEvo platform, we performed a hierarchical-based population structure analyses from available genomes of three distinct pathogenic bacterial species as individual case studies. The specific case studies illustrate how hierarchical analyses of population structures, genotype frequencies, and distribution of specific gene functions can be integrated into an analysis. Collectively, our study shows that ProkEvo presents a practical viable option for scalable, automated analyses of bacterial populations with direct applications for basic microbiology research, clinical microbiological diagnostics, and epidemiological surveillance.