Hierarchical and spatially explicit clustering of DNA sequences with BAPS software

Population genomics provides insights into the evolution and adaptation to humans of the waterborne pathogen Mycobacterium kansasii

Mycobacterium kansasii can cause serious pulmonary disease. It belongs to a group of closely-related species of non-tuberculous mycobacteria known as the M. kansasii complex (MKC). Here, we report a population genomics analysis of 358 MKC isolates from worldwide water and clinical sources. We find that recombination, likely mediated by distributive conjugative transfer, has contributed to speciation and on-going diversification of the MKC. Our analyses support municipal water as a main source of MKC infections. Furthermore, nearly 80% of the MKC infections are due to closely-related M. kansasii strains, forming a main cluster that apparently originated in the 1900s and subsequently expanded globally. Bioinformatic analyses indicate that several genes involved in metabolism (e.g., maintenance of the methylcitrate cycle), ESX-I secretion, metal ion homeostasis and cell surface remodelling may have contributed to M. kansasii's success and its ongoing adaptation to the human host.

Zymomonas diversity and potential for biofuel production

BACKGROUND

Zymomonas mobilis is an aerotolerant α-proteobacterium, which has been genetically engineered for industrial purposes for decades. However, a comprehensive comparison of existing strains on the genomic level in conjunction with phenotype analysis has yet to be carried out. We here performed whole-genome comparison of 17 strains including nine that were sequenced in this study. We then compared 15 available Zymomonas strains for their natural abilities to perform under conditions relevant to biofuel synthesis. We tested their growth in anaerobic rich media, as well as growth, ethanol production and xylose utilization in lignocellulosic hydrolysate. We additionally compared their tolerance to isobutanol, flocculation characteristics, and ability to uptake foreign DNA by electroporation and conjugation.

RESULTS

Using clustering based on 99% average nucleotide identity (ANI), we classified 12 strains into four clusters based on sequence similarity, while five strains did not cluster with any other strain. Strains belonging to the same 99% ANI cluster showed similar performance while significant variation was observed between the clusters. Overall, conjugation and electroporation efficiencies were poor across all strains, which was consistent with our finding of coding potential for several DNA defense mechanisms, such as CRISPR and restriction-modification systems, across all genomes. We found that strain ATCC31821 (ZM4) had a more diverse plasmid profile than other strains, possibly leading to the unique phenotypes observed for this strain. ZM4 also showed the highest growth of any strain in both laboratory media and lignocellulosic hydrolysate and was among the top 3 strains for isobutanol tolerance and electroporation and conjugation efficiency.

CONCLUSIONS

Our findings suggest that strain ZM4 has a unique combination of genetic and phenotypic traits that are beneficial for biofuel production and propose investing future efforts in further engineering of ZM4 for industrial purposes rather than exploring new Zymomonas isolates.

The Population Structure of Borrelia lusitaniae Is Reflected by a Population Division of Its Ixodes Vector

Populations of vector-borne pathogens are shaped by the distribution and movement of vector and reservoir hosts. To study what impact host and vector association have on tick-borne pathogens, we investigated the population structure of Borrelia lusitaniae using multilocus sequence typing (MLST). Novel sequences were acquired from questing ticks collected in multiple North African and European locations and were supplemented by publicly available sequences at the Borrelia Pubmlst database (accessed on 11 February 2020). Population structure of B. lusitaniae was inferred using clustering and network analyses. Maximum likelihood phylogenies for two molecular tick markers (the mitochondrial 16S rRNA locus and a nuclear locus, Tick-receptor of outer surface protein A, trospA) were used to confirm the morphological species identification of collected ticks. Our results confirmed that B. lusitaniae does indeed form two distinguishable populations: one containing mostly European samples and the other mostly Portuguese and North African samples. Of interest, Portuguese samples clustered largely based on being from north (European) or south (North African) of the river Targus. As two different Ixodes species (i.e., I. ricinus and I. inopinatus) may vector Borrelia in these regions, reference samples were included for I. inopinatus but did not form monophyletic clades in either tree, suggesting some misidentification. Even so, the trospA phylogeny showed a monophyletic clade containing tick samples from Northern Africa and Portugal south of the river Tagus suggesting a population division in Ixodes on this locus. The pattern mirrored the clustering of B. lusitaniae samples, suggesting a potential co-evolution between tick and Borrelia populations that deserve further investigation.

Occurrence and Transmission of bla NDM-Carrying Enterobacteriaceae from Geese and the Surrounding Environment on a Commercial Goose Farm

We investigated the prevalence and transmission of NDM-producing Enterobacteriaceae in fecal samples of geese and environmental samples from a goose farm in southern China. The samples were cultivated on MacConkey agar plates supplemented with meropenem. Individual colonies were examined for bla _NDM, and bla _NDM-positive bacteria were characterized based on whole-genome sequencing (WGS) data from the Illumina and Oxford Nanopore Technologies (ONT) platforms. Of 117 samples analyzed, the carriage rates for New Delhi metallo-β-lactamase (NDM)-positive Enterobacteriaceae were 47.1, 18, and 50% in geese, inanimate environments (sewage, soil, fodder, and dust), and mouse samples, respectively. Two variants (bla _NDM-1 and bla _NDM-5, in 4 and 40 isolates, respectively) were found among 44 bla _NDM-positive Enterobacteriaceae; these variants belonged to eight species, and Escherichia coli was the most prevalent (50%). WGS analysis revealed that bla _NDM coexisted with diverse antibiotic resistance genes (ARGs). Population structure analysis showed that most E. coli and Enterobacter sp. isolates were highly heterogeneous, while most Citrobacter sp. and P. stuartii isolates possessed extremely high genetic similarities. In addition, bla _NDM-5-positive ST4358/ST48 E. coli isolates were found to be clonally spread between geese and the environment and were highly genetically similar to those reported from ducks, farm environments, and humans in China. Plasmid analysis indicated that IncX3 pHNYX644-1-like (n = 40) and untypeable pM2-1-like plasmids (n = 4) mediated bla _NDM spread. pM2-1-like plasmids possessed diverse ARGs, including bla _NDM-1, the arsenical and mercury resistance operons, and the maltose operon. Our findings revealed that the goose farm is a reservoir for NDM-positive Enterobacteriaceae The bla _NDM contamination of wild mice and the novel pM2-1-like plasmid described here likely adds to the risk for dissemination of bla _NDM and associated resistance genes.IMPORTANCE Carbapenem-resistant bacteria, in particular NDM-producing Enterobacteriaceae, have become a great threat to global public. These bacteria have been found not only in hospital and community environments but also among food animal production chains, which are recognized as reservoirs for NDM-producing Enterobacteriaceae However, the dissemination of NDM-producing bacteria in waterfowl farms has been less well explored. Our study demonstrates that the horizontal spread of bla _NDM-carrying plasmids and the partial clonal spread of bla _NDM-positive Enterobacteriaceae contribute to the widespread contamination of bla _NDM in the goose farm ecosystem, including mice. Furthermore, we found a novel and transferable bla _NDM-1-carrying multidrug resistance (MDR) plasmid that possessed multiple environmental adaptation-related genes. The outcomes of this study contribute to a better understanding of the prevalence and transmission of bla _NDM-carrying Enterobacteriaceae among diverse niches in the farm ecosystem.

Global phylogeography and invasion history of the spotted lanternfly revealed by mitochondrial phylogenomics

Biological invasion has been a serious global threat due to increasing international trade and population movements. Tracking the source and route of invasive species and evaluating the genetic differences in their native regions have great significance for the effective monitoring and management, and further resolving the invasive mechanism. The spotted lanternfly Lycorma delicatula is native to China and invaded South Korea, Japan, and the United States during the last decade, causing severe damages to the fruits and timber industries. However, its global phylogeographic pattern and invasion history are not clearly understood. We applied high-throughput sequencing to obtain 392 whole mitochondrial genome sequences from four countries to ascertain the origin, dispersal, and invasion history of the spotted lanternfly. Phylogenomic analyses revealed that the spotted lanternfly originated from southwestern China, diverged into six phylogeographic lineages, and experienced northward expansion across the Yangtze River in the late Pleistocene. South Korea populations were derived from multiple invasions from eastern China and Japan with two different genetic sources of northwestern (Loess Plateau) and eastern (East Plain) lineages in China, whereas the each of Japan and the United States had only one. The United States populations originated through single invasive event from South Korea, which served as a bridgehead of invasion. The environmental conditions, especially the distribution of host Ailanthus trees, and adaptability possibly account for the rapid spread of the spotted lanternfly in the native and introduced regions.

The global dissemination of hospital clones of Enterococcus faecium

BACKGROUND

The hospital-adapted A1 group of Enterococcus faecium remains an organism of significant concern in the context of drug-resistant hospital-associated infections. How this pathogen evolves and disseminates remains poorly understood.

METHODS

A large, globally representative collection of short-read genomic data from the hospital-associated A1 group of Enterococcus faecium was assembled (n = 973). We analysed, using a novel analysis approach, global diversity in terms of both the dynamics of the accessory genome and homologous recombination among conserved genes.

RESULTS

Two main modes of genomic evolution continue to shape E. faecium: the acquisition and loss of genes, including antimicrobial resistance genes, through mobile genetic elements including plasmids, and homologous recombination of the core genome. These events lead to new clones emerging at the local level, followed by the erosion of signals of clonality through recombination, and in some identifiable cases producing new clonal clusters. These patterns lead to new, emerging lineages which are able to spread globally over relatively short timeframes.

CONCLUSIONS

The ability of A1 E. faecium to continually present new combinations of genes for potential selection suggests that controlling this pathogen will remain challenging but establishing a framework for understanding genomic evolution is likely to aid in tracking the threats posed by newly emerging lineages.

Evolutionary genomic and bacteria GWAS analysis of Mycobacterium avium subsp. paratuberculosis and dairy cattle Johne's disease phenotypes

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of Johne's disease in ruminants, which has important health consequences for dairy cattle. The Regional Dairy Quality Management Alliance (RDQMA) project is a multistate research program involving MAP isolates taken from three intensively studied commercial dairy farms in the northeastern United States, which emphasized longitudinal data collection of both MAP isolates and animal health in three regional dairy herds for a period of about 7 years. This paper reports the results of a pan-GWAS analysis involving 318 MAP isolates and dairy cow Johne's disease phenotypes, taken from these three farms. Based on our highly curated accessory gene count the pan-GWAS analysis identified several MAP genes associated with bovine Johne's disease phenotypes scored from these three farms, with some of the genes having functions suggestive of possible cause/effect relationships to these phenotypes. This paper reports a pan-genomic comparative analysis between MAP and Mycobacterium tuberculosis, assessing functional Gene Ontology category enrichments between these taxa. Finally, we also provide a population genomic perspective on the effectiveness of herd isolation, involving closed dairy farms, in preventing MAP inter-farm cross infection on a micro-geographic scale.IMPORTANCE Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of Johne's disease in ruminants, which has important health consequences for dairy cattle, and enormous economic consequences for the dairy industry. Understanding which genes in this bacterium are correlated with key disease phenotypes can lead to functional experiments targeting these genes and ultimately lead to improved control strategies. This study represents a rare example of a prolonged longitudinal study of dairy cattle where the disease was measured and the bacteria were isolated from the same cows. The genome sequences of over 300 MAP isolates were analyzed for genes that were correlated with a wide range of Johne's disease phenotypes. A number of genes were identified that were significantly associated with several aspects of the disease and suggestive of further experimental follow-up.

Evolutionary Dynamics Based on Comparative Genomics of Pathogenic Escherichia coli Lineages Harboring Polyketide Synthase (pks) Island

The genotoxin colibactin is a secondary metabolite produced by the polyketide synthase (pks) island harbored by extraintestinal pathogenic E. coli (ExPEC) and other members of the Enterobacteriaceae that has been increasingly reported to have critical implications in human health. The present study entails a high-throughput whole-genome comparison and phylogenetic analysis of such pathogenic E. coli isolates to gain insights into the patterns of distribution, horizontal transmission, and evolution of the island. For the current study, 23 pks-positive ExPEC genomes were newly sequenced, and their virulome and resistome profiles indicated a preponderance of virulence encoding genes and a reduced number of genes for antimicrobial resistance. In addition, 4,090 E. coli genomes from the public domain were also analyzed for large-scale screening for pks-positive genomes, out of which a total of 530 pks-positive genomes were studied to understand the subtype-based distribution pattern(s). The pks island showed a significant association with the B2 phylogroup (82.2%) and a high prevalence in sequence type 73 (ST73; n = 179) and ST95 (n = 110) and the O6:H1 (n = 110) serotype. Maximum-likelihood (ML) phylogeny of the core genome and intergenic regions (IGRs) of the ST95 model data set, which was selected because it had both pks-positive and pks-negative genomes, displayed clustering in relation to their carriage of the pks island. Prevalence patterns of genes encoding RM systems in the pks-positive and pks-negative genomes were also analyzed to determine their potential role in pks island acquisition and the maintenance capability of the genomes. Further, the maximum-likelihood phylogeny based on the core genome and pks island sequences from 247 genomes with an intact pks island demonstrated horizontal gene transfer of the island across sequence types and serotypes, with few exceptions. This study vitally contributes to understanding of the lineages and subtypes that have a higher propensity to harbor the pks island-encoded genotoxin with possible clinical implications.

Molecular Epidemiology and Clone Transmission of Carbapenem-Resistant Acinetobacter baumannii in ICU Rooms

Carbapenem-resistant Acinetobacter baumannii (CRAB) is a major cause of nosocomial infections and hospital outbreaks worldwide, remaining a critical clinical concern. Here we characterized and investigated the phylogenetic relationships of 105 CRAB isolates from an intensive care unit from one hospital in China collected over six years. All strains carried bla_OXA-23, bla_OXA-66 genes for carbapenem resistance, also had high resistance gene, virulence factor, and insertion sequence burdens. Whole-genome sequencing revealed all strains belonged to ST2, the global clone CC2. The phylogenetic analysis based on the core genome showed all isolates were dominated by a single lineage of three clusters and eight different clones. Two clones were popular during the collection time. Using chi-square test to identify the epidemiologically meaningful groupings, we found the significant difference in community structure only existed in strains from separation time. The haplotype and median-joining network analysis revealed genetic differences appeared among clusters and changes occurred overtime in the dominating cluster. Our results highlighted substantial multidrug-resistant CRAB burden in the hospital ICU environment demonstrating potential clone outbreak in the hospital.

Genome-Wide Association Studies for the Detection of Genetic Variants Associated With Daptomycin and Ceftaroline Resistance in Staphylococcus aureus

Background

As next generation sequencing (NGS) technologies have experienced a rapid development over the last decade, the investigation of the bacterial genetic architecture reveals a high potential to dissect causal loci of antibiotic resistance phenotypes. Although genome-wide association studies (GWAS) have been successfully applied for investigating the basis of resistance traits, complex resistance phenotypes have been omitted so far. For S. aureus this especially refers to antibiotics of last resort like daptomycin and ceftaroline. Therefore, we aimed to perform GWAS for the identification of genetic variants associated with DAP and CPT resistance in clinical S. aureus isolates.

Materials/methods

To conduct microbial GWAS, we selected cases and controls according to their clonal background, date of isolation, and geographical origin. Association testing was performed with PLINK and SEER analysis. By using in silico analysis, we also searched for rare genetic variants in candidate loci that have previously been described to be involved in the development of corresponding resistance phenotypes.

Results

GWAS revealed MprF P314L and L826F to be significantly associated with DAP resistance. These mutations were found to be homogenously distributed among clonal lineages suggesting convergent evolution. Additionally, rare and yet undescribed single nucleotide polymorphisms could be identified within mprF and putative candidate genes. Finally, we could show that each DAP resistant isolate exhibited at least one amino acid substitution within the open reading frame of mprF. Due to the presence of strong population stratification, no genetic variants could be associated with CPT resistance. However, the investigation of the staphylococcal cassette chromosome mec (SCCmec) revealed various mecA SNPs to be putatively linked with CPT resistance. Additionally, some CPT resistant isolates revealed no mecA mutations, supporting the hypothesis that further and still unknown resistance determinants are crucial for the development of CPT resistance in S. aureus.

Conclusion

We hereby confirmed the potential of GWAS to identify genetic variants that are associated with antibiotic resistance traits in S. aureus. However, precautions need to be taken to prevent the detection of spurious associations. In addition, the implementation of different approaches is still essential to detect multiple forms of variations and mutations that occur with a low frequency.

Wolbachia affects mitochondrial population structure in two systems of closely related Palaearctic blue butterflies

The bacterium Wolbachia infects many insect species and spreads by diverse vertical and horizontal means. As co-inherited organisms, these bacteria often cause problems in mitochondrial phylogeny inference. The phylogenetic relationships of many closely related Palaearctic blue butterflies (Lepidoptera: Lycaenidae: Polyommatinae) are ambiguous. We considered the patterns of Wolbachia infection and mitochondrial diversity in two systems: Aricia agestis/Aricia artaxerxes and the Pseudophilotes baton species complex. We sampled butterflies across their distribution ranges and sequenced one butterfly mitochondrial gene and two Wolbachia genes. Both butterfly systems had uninfected and infected populations, and harboured several Wolbachia strains. Wolbachia was highly prevalent in A. artaxerxes and the host's mitochondrial structure was shallow, in contrast to A. agestis. Similar bacterial alleles infected both Aricia species from nearby sites, pointing to a possible horizontal transfer. Mitochondrial history of the P. baton species complex mirrored its Wolbachia infection and not the taxonomical division. Pseudophilotes baton and P. vicrama formed a hybrid zone in Europe. Wolbachia could obscure mitochondrial history, but knowledge on the infection helps us to understand the observed patterns. Testing for Wolbachia should be routine in mitochondrial DNA studies.

Population genetics and evolutionary history of the endangered Eld's deer (Rucervus eldii) with implications for planning species recovery

Eld's deer (Rucervus eldii) with three recognised subspecies (R. e. eldii, R. e. thamin, and R. e. siamensis) represents one of the most threatened cervids found in Southeast Asia. The species has experienced considerable range contractions and local extinctions owing to habitat loss and fragmentation, hunting, and illegal trade across its distribution range over the last century. Understanding the patterns of genetic variation is crucial for planning effective conservation strategies. This study investigated the phylogeography, divergence events and systematics of Eld's deer subspecies using the largest mtDNA dataset compiled to date. We also analysed the genetic structure and demographic history of R. e. eldii using 19 microsatellite markers. Our results showed that R. e. siamensis exhibits two divergent mtDNA lineages (mainland and Hainan Island), which diverged around 0.2 Mya (95% HPD 0.1-0.2), possibly driven by the fluctuating sea levels of the Early Holocene period. The divergence between R. e. eldii and R. e. siamensis occurred around 0.4 Mya (95% HPD 0.3-0.5), potentially associated with the adaptations to warm and humid climate with open grassland vegetation that predominated the region. Furthermore, R. e. eldii exhibits low levels of genetic diversity and small contemporary effective population size (median = 7, 4.7-10.8 at 95% CI) with widespread historical genetic bottlenecks which accentuates its vulnerability to inbreeding and extinction. Based on the observed significant evolutionary and systematic distance between Eld's deer and other species of the genus Rucervus, we propose to classify Eld's deer (Cervus eldii) in the genus Cervus, which is in congruent with previous phylogenetic studies. This study provides important conservation implications required to direct the ongoing population recovery programs and planning future conservation strategies.

Mode and dynamics of vanA-type vancomycin resistance dissemination in Dutch hospitals

BACKGROUND

Enterococcus faecium is a commensal of the gastrointestinal tract of animals and humans but also a causative agent of hospital-acquired infections. Resistance against glycopeptides and to vancomycin has motivated the inclusion of E. faecium in the WHO global priority list. Vancomycin resistance can be conferred by the vanA gene cluster on the transposon Tn1546, which is frequently present in plasmids. The vanA gene cluster can be disseminated clonally but also horizontally either by plasmid dissemination or by Tn1546 transposition between different genomic locations.

METHODS

We performed a retrospective study of the genomic epidemiology of 309 vancomycin-resistant E. faecium (VRE) isolates across 32 Dutch hospitals (2012-2015). Genomic information regarding clonality and Tn1546 characterization was extracted using hierBAPS sequence clusters (SC) and TETyper, respectively. Plasmids were predicted using gplas in combination with a network approach based on shared k-mer content. Next, we conducted a pairwise comparison between isolates sharing a potential epidemiological link to elucidate whether clonal, plasmid, or Tn1546 spread accounted for vanA-type resistance dissemination.

RESULTS

On average, we estimated that 59% of VRE cases with a potential epidemiological link were unrelated which was defined as VRE pairs with a distinct Tn1546 variant. Clonal dissemination accounted for 32% cases in which the same SC and Tn1546 variants were identified. Horizontal plasmid dissemination accounted for 7% of VRE cases, in which we observed VRE pairs belonging to a distinct SC but carrying an identical plasmid and Tn1546 variant. In 2% of cases, we observed the same Tn1546 variant in distinct SC and plasmid types which could be explained by mixed and consecutive events of clonal and plasmid dissemination.

CONCLUSIONS

In related VRE cases, the dissemination of the vanA gene cluster in Dutch hospitals between 2012 and 2015 was dominated by clonal spread. However, we also identified outbreak settings with high frequencies of plasmid dissemination in which the spread of resistance was mainly driven by horizontal gene transfer (HGT). This study demonstrates the feasibility of distinguishing between modes of dissemination with short-read data and provides a novel assessment to estimate the relative contribution of nested genomic elements in the dissemination of vanA-type resistance.

Metagenomics-Based Approach to Source-Attribution of Antimicrobial Resistance Determinants - Identification of Reservoir Resistome Signatures

Metagenomics can unveil the genetic content of the total microbiota in different environments, such as food products and the guts of humans and livestock. It is therefore considered of great potential to investigate the transmission of foodborne hazards as part of source-attribution studies. Source-attribution of antimicrobial resistance (AMR) has traditionally relied on pathogen isolation, while metagenomics allows investigating the full span of AMR determinants. In this study, we hypothesized that the relative abundance of fecal resistome components can be associated with specific reservoirs, and that resistomes can be used for AMR source-attribution. We used shotgun-sequences from fecal samples of pigs, broilers, turkeys- and veal calves collected across Europe, and fecal samples from humans occupationally exposed to livestock in one country (pig slaughterhouse workers, pig and broiler farmers). We applied both hierarchical and flat forms of the supervised classification ensemble algorithm Random Forests to classify resistomes into corresponding reservoir classes. We identified country-specific and -independent AMR determinants, and assessed the impact of country-specific determinants when attributing AMR resistance in humans. Additionally, we performed a similarity percentage analysis with the full spectrum of AMR determinants to identify resistome signatures for the different reservoirs. We showed that the number of AMR determinants necessary to attribute a resistome into the correct reservoir increases with a larger reservoir heterogeneity, and that the impact of country-specific resistome signatures on prediction varies between countries. We predicted a higher occupational exposure to AMR determinants among workers exposed to pigs than among those exposed to broilers. Additionally, results suggested that AMR exposure on pig farms was higher than in pig slaughterhouses. Human resistomes were more similar to pig and veal calves’ resistomes than to those of broilers and turkeys, and the majority of these resistome dissimilarities can be explained by a small set of AMR determinants. We identified resistome signatures for each individual reservoir, which include AMR determinants significantly associated with on-farm antimicrobial use. We attributed human resistomes to different livestock reservoirs using Random Forests, which allowed identifying pigs as a potential source of AMR in humans. This study thus demonstrates that it is possible to apply metagenomics in AMR source-attribution.

Early Signals of Vaccine-driven Perturbation Seen in Pneumococcal Carriage Population Genomic Data

Background

Pneumococcal conjugate vaccines (PCVs) have reduced pneumococcal diseases globally. Pneumococcal genomic surveys elucidate PCV effects on population structure but are rarely conducted in low-income settings despite the high disease burden.

Methods

We undertook whole-genome sequencing (WGS) of 660 pneumococcal isolates collected through surveys from healthy carriers 2 years from 13-valent PCV (PCV13) introduction and 1 year after rollout in northern Malawi. We investigated changes in population structure, within-lineage serotype dynamics, serotype diversity, and frequency of antibiotic resistance (ABR) and accessory genes.

Results

In children <5 years of age, frequency and diversity of vaccine serotypes (VTs) decreased significantly post-PCV, but no significant changes occurred in persons ≥5 years of age. Clearance of VT serotypes was consistent across different genetic backgrounds (lineages). There was an increase of nonvaccine serotypes (NVTs)—namely 7C, 15B/C, and 23A—in children <5 years of age, but 28F increased in both age groups. While carriage rates have been recently shown to remain stable post-PCV due to replacement serotypes, there was no change in diversity of NVTs. Additionally, frequency of intermediate-penicillin-resistant lineages decreased post-PCV. Although frequency of ABR genes remained stable, other accessory genes, especially those associated with mobile genetic element and bacteriocins, showed changes in frequency post-PCV.

Conclusions

We demonstrate evidence of significant population restructuring post-PCV driven by decreasing frequency of vaccine serotypes and increasing frequency of few NVTs mainly in children under 5. Continued surveillance with WGS remains crucial to fully understand dynamics of the residual VTs and replacement NVT serotypes post-PCV.

The most important fishery resource in the Amazon, the migratory catfish Brachyplatystoma vaillantii (Siluriformes: Pimelodidae), is composed by an unique and genetically diverse population in the Solimões-Amazonas River System

ABSTRACT The migratory catfish Brachyplatystoma vaillantii is one of the most important fishery resources in the Amazon. Intense capture occurs associated to its life cycle. In order to know the genetic status, we sequenced the mitochondrial DNA control region from 150 individuals of B. vaillantii, collected in five fishing landing locations, covering the length of the Solimões-Amazonas River in Brazil. Genetic diversity parameters suggest there is no genetic differentiation between the five localities. Population’s expansion indicated by R 2 and Fu’s Fs tests was also confirmed by the high number of unique haplotypes found. The Analyses of molecular variance indicated that nearly all variability was contained within locations (99.86%), and estimates of gene flow among B. vaillantii were high (F ST = 0.0014). These results suggest that Brachyplatystoma vaillantii forms a panmitic population along the Solimões-Amazonas River and, has greater genetic variability than other species of the Brachyplatystoma genus available so far. Although the influence of different tributaries on B. vaillantii migration patterns remains uncertain, a single population in the main channel should be consider in future policies for management of this resource. However, since the species’ life cycle uses habitats in several countries, its management and conservation depend greatly of internationally joined efforts.

Population Genomics of Bacterial Plant Pathogens

Population genomics is transforming our understanding of pathogen biology and evolution, and contributing to the prevention and management of disease in diverse crops. We provide an overview of key methods in bacterial population genomics and describe recent work focusing on three topics of critical importance to plant pathology: (i) resolving pathogen origins and transmission pathways during outbreak events, (ii) identifying the genetic basis of host specificity and virulence, and (iii) understanding how pathogens evolve in response to changing agricultural practices.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Evolution of Ciprofloxacin Resistance-Encoding Genetic Elements in Salmonella

The incidence of ciprofloxacin resistance in Salmonella has increased dramatically in the past decade. To track the evolutionary trend of ciprofloxacin resistance-encoding genetic elements during this period, we surveyed the prevalence of Salmonella in food products in Shenzhen, China, during the period of 2012 to 2017 and performed whole-genome sequencing and genetic analysis of 566 ciprofloxacin-resistant clinical Salmonella strains collected during this survey. We observed that target gene mutations have become much less common, with single gyrA mutation currently detectable in Salmonella enterica serovar Typhimurium only. Multiple plasmid-mediated quinolone resistance (PMQR) genes located in the chromosome and plasmids are now frequently detectable in ciprofloxacin-resistant Salmonella strains of various serotypes. Among them, the qnrS1 gene was often harbored by multiple plasmids, with p10k-like plasmids being the most dominant. Importantly, p10k-like plasmids initially were not conjugative but became transmissible with the help of a helper plasmid. Ciprofloxacin resistance due to combined effect of carriage of the qnrS1 gene and other resistance mechanisms is common. In S Typhimurium, carriage of qnrS1 is often associated with a single gyrA mutation; in other serotypes, combination of qnrS1 and other PMQR genes located in the chromosomal fragment or plasmid is observed. Another major mechanism of ciprofloxacin resistance, mainly observable in S Derby, involves a chromosomal fragment harboring the qnrS2-aac(6')lb-cr-oqxAB elements. Intriguingly, this chromosomal fragment, flanked by IS26, could form a circular intermediate and became transferrable. To conclude, the increase in the incidence of various PMQR mobile genetic elements and their interactions with other resistance mechanism contribute to a sharp increase in the prevalence of ciprofloxacin-resistant clinical Salmonella strains in recent years.IMPORTANCE Resistance of nontyphoidal Salmonella to fluoroquinolones such as ciprofloxacin is known to be mediated by target mutations. This study surveyed the prevalence of Salmonella strains recovered from 2,989 food products in Shenzhen, China, during the period 2012 to 2017 and characterized the genetic features of several PMQR gene-bearing plasmids and ciprofloxacin resistance-encoding DNA fragments. The emergence of such genetic elements has caused a shift in the genetic location of ciprofloxacin resistance determinants from the chromosomal mutations to various mobile genetic elements. The distribution of these PMQR plasmids showed that they exhibited high serotype specificity, except for the p10k-like plasmids, which can be widely detected and efficiently transmitted among Salmonella strains of various serotypes by fusing to a new conjugative helper plasmid. The sharp increase in the prevalence of ciprofloxacin resistance in recent years may cause a predisposition to the emergence of multidrug-resistant Salmonella strains and pose huge challenges to public health and infection control efforts.

Meningococcal Disease-Associated Prophage-Like Elements Are Present in Neisseria gonorrhoeae and Some Commensal Neisseria Species

Neisseria spp. possess four genogroups of filamentous prophages, termed Nf1 to 4. A filamentous bacteriophage from the Nf1 genogroup termed meningococcal disease-associated phage (MDA φ) is associated with clonal complexes of Neisseria meningitidis that cause invasive meningococcal disease. Recently, we recovered an isolate of Neisseria gonorrhoeae (ExNg63) from a rare case of gonococcal meningitis, and found that it possessed a region with 90% similarity to Nf1 prophages, specifically, the meningococcal MDA φ. This led to the hypothesis that the Nf1 prophage may be more widely distributed amongst the genus Neisseria. An analysis of 92 reference genomes revealed the presence of intact Nf1 prophages in the commensal species, Neisseria lactamica and Neisseria cinerea in addition to the pathogen N. gonorrhoeae. In N. gonorrhoeae, Nf1 prophages had a restricted distribution but were present in all representatives of MLST ST1918. Of the 160 phage integration sites identified, only one common insertion site was found between one isolate of N. gonorrhoeae and N. meningitidis. There was an absence of any obvious conservation of the receptor for prophage entry, PilE, suggesting that the phage may have been obtained by natural transformation. An examination of the restriction modification systems and mutated mismatch repair systems with prophage presence suggested that there was no obvious preference for these hosts. A timed phylogeny inferred that N. meningitidis was the donor of the Nf1 prophages in N. lactamica and N. gonorrhoeae. Further work is required to determine whether Nf1 prophages are active and can act as accessory colonization factors in these species.

Mitochondrial DNAs provide insight into trypanosome phylogeny and molecular evolution

Background

Trypanosomes are single-celled eukaryotic parasites characterised by the unique biology of their mitochondrial DNA. African livestock trypanosomes impose a major burden on agriculture across sub-Saharan Africa, but are poorly understood compared to those that cause sleeping sickness and Chagas disease in humans. Here we explore the potential of the maxicircle, a component of trypanosome mitochondrial DNA to study the evolutionary history of trypanosomes.

Results

We used long-read sequencing to completely assemble maxicircle mitochondrial DNA from four previously uncharacterized African trypanosomes, and leveraged these assemblies to scaffold and assemble a further 103 trypanosome maxicircle gene coding regions from published short-read data. While synteny was largely conserved, there were repeated, independent losses of Complex I genes. Comparison of pre-edited and non-edited genes revealed the impact of RNA editing on nucleotide composition, with non-edited genes approaching the limits of GC loss. African tsetse-transmitted trypanosomes showed high levels of RNA editing compared to other trypanosomes. The gene coding regions of maxicircle mitochondrial DNAs were used to construct time-resolved phylogenetic trees, revealing deep divergence events among isolates of the pathogens Trypanosoma brucei and T. congolense.

Conclusions

Our data represents a new resource for experimental and evolutionary analyses of trypanosome phylogeny, molecular evolution and function. Molecular clock analyses yielded a timescale for trypanosome evolution congruent with major biogeographical events in Africa and revealed the recent emergence of Trypanosoma brucei gambiense and T. equiperdum, major human and animal pathogens.

Unexpected Gene-Flow in Urban Environments: The Example of the European Hedgehog

Simple Summary

An urban environment holds many barriers for mammals with limited mobility such as hedgehogs. These barriers appear often unsurmountable (e.g., rivers, highways, fences) and thus hinder contact between hedgehogs, leading to genetic isolation. In our study we tested whether these barriers affect the hedgehog population of urban Berlin, Germany. As Berlin has many of these barriers, we were expecting a strong genetic differentiation among hedgehog populations. However, when we looked at unrelated individuals, we did not see genetic differentiation among populations. The latter was only detected when we included related individuals too, a ‘family clan’ structure that is referred to as gamodemes. We conclude that the high percentage of greenery in Berlin provides sufficient habitat for hedgehogs to maintain connectivity across the city.

Abstract

We use the European hedgehog (Erinaceus europaeus), a mammal with limited mobility, as a model species to study whether the structural matrix of the urban environment has an influence on population genetic structure of such species in the city of Berlin (Germany). Using ten established microsatellite loci we genotyped 143 hedgehogs from numerous sites throughout Berlin. Inclusion of all individuals in the cluster analysis yielded three genetic clusters, likely reflecting spatial associations of kin (larger family groups, known as gamodemes). To examine the potential bias in the cluster analysis caused by closely related individuals, we determined all pairwise relationships and excluded close relatives before repeating the cluster analysis. For this data subset (N = 65) both clustering algorithms applied (Structure, Baps) indicated the presence of a single genetic cluster. These results suggest that the high proportion of green patches in the city of Berlin provides numerous steppingstone habitats potentially linking local subpopulations. Alternatively, translocation of individuals across the city by hedgehog rescue facilities may also explain the existence of only a single cluster. We therefore propose that information about management activities such as releases by animal rescue centres should include location data (as exactly as possible) regarding both the collection and the release site, which can then be used in population genetic studies.

Culture-free genome-wide locus sequence typing (GLST) provides new perspectives on Trypanosoma cruzi dispersal and infection complexity

Analysis of genetic polymorphism is a powerful tool for epidemiological surveillance and research. Powerful inference from pathogen genetic variation, however, is often restrained by limited access to representative target DNA, especially in the study of obligate parasitic species for which ex vivo culture is resource-intensive or bias-prone. Modern sequence capture methods enable pathogen genetic variation to be analyzed directly from host/vector material but are often too complex and expensive for resource-poor settings where infectious diseases prevail. This study proposes a simple, cost-effective 'genome-wide locus sequence typing' (GLST) tool based on massive parallel amplification of information hotspots throughout the target pathogen genome. The multiplexed polymerase chain reaction amplifies hundreds of different, user-defined genetic targets in a single reaction tube, and subsequent agarose gel-based clean-up and barcoding completes library preparation at under 4 USD per sample. Our study generates a flexible GLST primer panel design workflow for Trypanosoma cruzi, the parasitic agent of Chagas disease. We successfully apply our 203-target GLST panel to direct, culture-free metagenomic extracts from triatomine vectors containing a minimum of 3.69 pg/μl T. cruzi DNA and further elaborate on method performance by sequencing GLST libraries from T. cruzi reference clones representing discrete typing units (DTUs) TcI, TcIII, TcIV, TcV and TcVI. The 780 SNP sites we identify in the sample set repeatably distinguish parasites infecting sympatric vectors and detect correlations between genetic and geographic distances at regional (< 150 km) as well as continental scales. The markers also clearly separate TcI, TcIII, TcIV and TcV + TcVI and appear to distinguish multiclonal infections within TcI. We discuss the advantages, limitations and prospects of our method across a spectrum of epidemiological research.

Phylogenomic Analysis of Campylobacter fetus Reveals a Clonal Structure of Insertion Element ISCfe1 Positive Genomes

Subspecies of the species Campylobacter fetus are associated with specific host niches including mammals and reptiles. Campylobacter fetus subsp. fetus is a zoonotic pathogen infecting humans. Infections can vary from an acute intestinal illness to severe systemic infections, with sheep and cattle as major reservoirs. In contrast, Campylobacter fetus subsp. venerealis causes bovine genital campylobacteriosis, which leads to abortion in cattle and a high economic burden for the farmers. Therefore, high-quality molecular subtyping is indispensable for interventional epidemiology. We used whole-genome sequencing (WGS) data of 283 Campylobacter fetus strains from 18 countries and compared several methods for Campylobacter fetus subtyping, including WGS, multilocus sequence typing, PCR assays, and the presence of the insertion element ISCfe1. We identified a highly clonal clade (designated as clade 1) that harbors the insertion sequence ISCfe1. The presence of this insertion sequence is an essential diagnostic tool for the identification of the subspecies Campylobacter fetus subsp. venerealis, serving as a target for several PCR assays. However, we have found a high sequence variability for the ISCfe1 besides the presence of ISCfe1-paralogues in certain other genomes (n = 7) which may cause incorrect diagnostic results. Clade 1 seems to be the cattle-specific clade of this species. We propose that only this clade might be designated as Campylobacter fetus subsp. venerealis as it harbors the ISCfe1 marker sequence, which is a major target for molecular methods currently used for Campylobacter fetus subspecies identification. Fostering this proposal, we defined eleven stable nucleotide markers specific for this clade. Additionally, we developed a bioinformatics toolbox for the fast identification of this clade based on WGS data. In conclusion, our results demonstrate that WGS can be used for Campylobacter fetus subtyping overcoming limitations of current PCR and MLST protocols.

Host adaptation and microbial competition drive Ralstonia solanacearum phylotype I evolution in the Republic of Korea

Bacterial wilt caused by the Ralstonia solanacearum species complex (RSSC) threatens the cultivation of important crops worldwide. We sequenced 30 RSSC phylotype I (R. pseudosolanacearum) strains isolated from pepper (Capsicum annuum) and tomato (Solanum lycopersicum) across the Republic of Korea. These isolates span the diversity of phylotype I, have extensive effector repertoires and are subject to frequent recombination. Recombination hotspots among South Korean phylotype I isolates include multiple predicted contact-dependent inhibition loci, suggesting that microbial competition plays a significant role in Ralstonia evolution. Rapid diversification of secreted effectors presents challenges for the development of disease-resistant plant varieties. We identified potential targets for disease resistance breeding by testing for allele-specific host recognition of T3Es present among South Korean phyloype I isolates. The integration of pathogen population genomics and molecular plant pathology contributes to the development of location-specific disease control and development of plant cultivars with durable resistance to relevant threats.

Genetic structure of recently fragmented suburban populations of European stag beetle

Habitat loss and fragmentation due to urbanization can negatively affect metapopulation persistence when gene flow among populations is reduced and population sizes decrease. Inference of patterns and processes of population connectivity derived from spatial genetic analysis has proven invaluable for conservation and management. However, a more complete account of population dynamics may be obtained by combining spatial and temporal sampling. We, therefore, performed a genetic study on European stag beetle (Lucanus cervus L.) populations in a suburban context using samples collected in three locations and during the period 2002-2016. The sampling area has seen recent landscape changes which resulted in population declines. Through the use of a suite of F ST, clustering analysis, individual assignment, and relatedness analysis, we assessed fine scale spatiotemporal genetic variation within and among habitat patches using 283 individuals successfully genotyped at 17 microsatellites. Our findings suggested the three locations to hold demographically independent populations, at least over time scales of relevance to conservation, though with higher levels of gene flow in the past. Contrary to expectation from tagging studies, dispersal appeared to be mainly female-biased. Although the life cycle of stag beetle suggests its generations to be discrete, no clear temporal structure was identified, which could be attributed to the varying duration of larval development. Since population bottlenecks were detected and estimates of effective number of breeders were low, conservation actions are eminent which should include the establishment of suitable dead wood for oviposition on both local and regional scales to increase (re)colonization success and connectivity among current populations.

Comparative Genomic Analysis Reveals Genetic Mechanisms of the Variety of Pathogenicity, Antibiotic Resistance, and Environmental Adaptation of Providencia Genus

The bacterial genus Providencia is Gram-negative opportunistic pathogens, which have been isolated from a variety of environments and organisms, ranging from humans to animals. Providencia alcalifaciens, Providencia rettgeri, and Providencia stuartii are the most common clinical isolates, however, these three species differ in their pathogenicity, antibiotic resistance and environmental adaptation. Genomes of 91 isolates of the genus Providencia were investigated to clarify their genetic diversity, focusing on virulence factors, antibiotic resistance genes, and environmental adaptation genes. Our study revealed an open pan-genome for the genus Providencia containing 14,720 gene families. Species of the genus Providencia exhibited different functional constraints, with the core genes, accessory genes, and unique genes. A maximum-likelihood phylogeny reconstructed with concatenated single-copy core genes classified all Providencia isolates into 11 distant groups. Comprehensive and systematic comparative genomic analyses revealed that specific distributions of virulence genes, which were highly homologous to virulence genes of the genus Proteus, contributed to diversity in pathogenicity of Providencia alcalifaciens, Providencia rettgeri, and Providencia stuartii. Furthermore, multidrug resistance (MDR) phenotypes of isolates of Providencia rettgeri and Providencia stuartii were predominantly due to resistance genes from class 1 and 2 integrons. In addition, Providencia rettgeri and Providencia stuartii harbored more genes related to material transport and energy metabolism, which conferred a stronger ability to adapt to diverse environments. Overall, our study provided valuable insights into the genetic diversity and functional features of the genus Providencia, and revealed genetic mechanisms underlying diversity in pathogenicity, antibiotic resistance and environmental adaptation of members of this genus.

Are population isolations and declines a threat to island endemic water striders? A lesson from demographic and niche modelling of Metrocoris esakii (Hemiptera: Gerridae)

Genetic stochasticity and bottlenecking in the course of Pleistocene glaciations have been identified as threatening the survival of local endemics. However, the mechanisms by which local endemic species balance the influences of these two events remain poorly understood. Here, we generated a double-digest restriction site-associated DNA sequencing (ddRAD-seq) data set, mined mitochondrial sequences and constructed ecological niche models for the island endemic water strider Metrocoris esakii (Hemiptera: Gerridae). We found that M. esakii comprised three divergent lineages (i.e., north, central and south) isolated by geographical barriers and generally experienced population declines with the constriction of suitable areas during the Last Glacial Maximum (LGM). Further demographic model testing and stairway plots revealed a history of recent gene flow among the neighbouring lineages and rapid recovery at the end of the LGM, indicating that M. esakii at least had the potential for an adaptive response to population fragmentation and bottlenecking. The northern lineage did not show genetic bottlenecking during the LGM, which was probably due to its large effective population size (N_e ) from migration, which improved its adaptive potential. Relative to the ddRAD-seq data set, the demographic results based on mitochondrial sequences were less conclusive, showing weak differentiation and oversimplified demographic trajectories for the three genetic lineages. Overall, this study provides some degree of optimism for the survival of island endemic water striders from a demographic perspective, but further evaluation of their extinction risk under the impacts of human activities is required.

Evolutionary history and current distribution of the West Mediterranean lineage of Brucella melitensis in Italy

Ovine and caprine brucellosis, caused by Brucella melitensis, is one of the world’s most widespread zoonoses and is a major cause of economic losses in domestic ruminant production. In Italy, the disease remains endemic in several southern provinces, despite an ongoing brucellosis eradication programme. In this study, we used whole-genome sequencing to detail the genetic diversity of circulating strains, and to examine the origins of the predominant sub-lineages of B. melitensis in Italy. We reconstructed a global phylogeny of B. melitensis, strengthened by 339 new whole-genome sequences, from Italian isolates collected from 2011 to 2018 as part of a national livestock surveillance programme. All Italian strains belonged to the West Mediterranean lineage, which further divided into two major clades that diverged roughly between the 5th and 7th centuries. We observed that Sicily serves as a brucellosis burden hotspot, giving rise to several distinct sub-lineages. More than 20 putative outbreak clusters of ovine and caprine brucellosis were identified, several of which persisted over the 8 year survey period despite an aggressive brucellosis eradication campaign. While the outbreaks in Central and Northern Italy were generally associated with introductions of single clones of B. melitensis and their subsequent dissemination within neighbouring territories, we observed weak geographical segregation of genotypes in the southern regions. Biovar determination, recommended in routine analysis of all Brucella strains by the World Organisation for Animal Health (OIE), could not discriminate among the four main global clades. This demonstrates a need for updating the guidelines used for monitoring B. melitensis transmission and spread, both at the national and international level, and to include whole-genome-based typing as the principal method for identification and tracing of brucellosis outbreaks.

A relic of the past: current genetic patterns of the palaeoendemic tree Nothofagus macrocarpa were shaped by climatic oscillations in central Chile

BACKGROUND AND AIMS

The Mediterranean-type forest of central Chile is considered a 'biodiversity hotspot' and a relic of a wider ancient distribution produced by past climatic oscillations. Nothofagus macrocarpa, commonly known as 'roble de Santiago', is a threatened palaeoendemic of this forest, poorly represented in the protected area system. This tree has been repeatedly misidentified as the sister species N. obliqua, which has affected its recognition and protection. Only a few populations of N. macrocarpa remain within a matrix of intensive land use that has been affected by recent forest fires. We tested the hypothesis that current populations of N. macrocarpa are a relic state of a previously widespread range, with the aim of contributing to its identification, its biogeographical history and the design of conservation measures using genetic information.

METHODS

We analysed remnant N. macrocarpa forests using nuclear (nDNA) and chloroplast DNA (cpDNA) sequences, conducted phylogenetic and phylogeographical analyses to reconstruct its biogeographical history, and assessed microsatellites [simple sequence repeats (SSRs)] to determine contemporary patters of diversity within and among all remnant populations. We also examined the degree of past, current and potential future isolation of N. macrocarpa populations using ecological niche models (ENMs).

KEY RESULTS

The species N. macrocarpa was confirmed by nDNA sequences, as previously suggested by chromosomal analysis. Small isolated populations of N. macrocarpa exhibited moderate to high genetic diversity according to SSRs. cpDNA analysis revealed a marked past latitudinal geographical structure, whereas analysis of SSRs did not find such current structure. ENM analyses revealed local expansion-contraction of the N. macrocarpa range during warmer periods, particularly in the northern and central ranges where basal-most cpDNA haplotypes were detected, and recent expansion to the south of the distribution.

CONCLUSIONS

Genetic patterns confirm that N. macrocarpa is a distinct species and suggest a marked latitudinal relic structure in at least two evolutionarily significant units, despite contemporary among-population gene flow. This information must be considered when choosing individuals (seeds and/or propagules) for restoration purposes, to avoid the admixture of divergent genetic stocks.

Occurrence of Brucella ceti in striped dolphins from Italian Seas

Brucella ceti infections have been increasingly reported in cetaceans, although a very limited characterization of Mediterranean Brucella spp. isolates has been previously reported and relatively few data exist about brucellosis among cetaceans in Italy. To address this gap, we studied 8 cases of B. ceti infection in striped dolphins (Stenella coeruleoalba) stranded along the Italian coastline from 2012 to 2018, investigated thanks to the Italian surveillance activity on stranded cetaceans. We focused on cases of stranding in eastern and western Italian seas, occurred along the Apulia (N = 6), Liguria (N = 1) and Calabria (N = 1) coastlines, through the analysis of gross and microscopic findings, the results of microbiological, biomolecular and serological investigations, as well as the detection of other relevant pathogens. The comparative genomic analysis used whole genome sequences of B. ceti from Italy paired with the publicly available complete genomes. Pathological changes consistent with B. ceti infection were detected in the central nervous system of 7 animals, showing non-suppurative meningoencephalitis. In 4 cases severe coinfections were detected, mostly involving Dolphin Morbillivirus (DMV). The severity of B. ceti-associated lesions supports the role of this microbial agent as a primary neurotropic pathogen for striped dolphins. We classified the 8 isolates into the common sequence type 26 (ST-26). Whole genome SNP analysis showed that the strains from Italy clustered into two genetically distinct clades. The first clade comprised exclusively the isolates from Ionian and Adriatic Seas, while the second one included the strain from the Ligurian Sea and those from the Catalonian coast. Plotting these clades onto the geographic map suggests a link between their phylogeny and topographical distribution. These results represent the first extensive characterization of B. ceti isolated from Italian waters reported to date and show the usefulness of WGS for understanding of the evolution of this emerging pathogen.

"Colonization routes and demographic history of Chirostoma humboldtianum in the central Mexican plateau"

In the present study we evaluate the population structure and potential colonization routes of the silverside Chirostoma humboldtianum through approximate Bayesian computations. Six microsatellite loci were amplified in a total of 288 individuals from six different locations covering the complete geographic distribution of the species. Additionally, two mitochondrial DNA markers, a D loop control region and cytochrome b were amplified in a subset of 107 individuals. The results found with microsatellites allow recovering well-structured populations that have experienced a drastic reduction in the effective population size. On the other hand, mtDNA sequences showed a moderate phylogeographic structure with shared haplotypes between geographic localities and signalsof a slight increase in the effective population size. Finally, the approximate Bayesian computation analysis performed with both datasets suggested a west-to-east colonization route for the species in Central Mexico.

Phylogeography and the population genetic structure of flowering cherry Cerasus serrulata (Rosaceae) in subtropical and temperate China

Cerasus serrulata (Rosaceae) is an important flowering cherry resource which is valuable for developing new cultivars of flowering cherries. It is broadly distributed and possesses abundant variations. In this study, phylogeographic analysis was conducted to reveal the evolutionary history to better understand the genetic diversity and genetic structure of C. serrulata so as to provide more accurate molecular insights into better conservation and utilization of the germplasm resources. A total of 327 individuals from 18 wild populations were collected. Three chloroplast DNA (cpDNA) fragments (matK, trnD-E, and trnS-G) and the nuclear internal transcribed spacer (ITS) were utilized. The results showed a high genetic diversity at both species level and population level of C. serrulata. High genetic differentiation and the existence of the phylogeographic structure were detected. No significant expansion events were discovered. Two geographic lineages were inferred. One was confined to the Qinling Mountains and the Taihang Mountains. The other was from the Wuling Mountains to the Jiangnan Hilly Regions and then went northeast to the coast of Asia. In addition, some taxonomic treatments of the C. serrulata complex are discussed and reconsidered. Conservation and utilization strategies of wild C. serrulata germplasm resources were recommended.

Frequency-dependent selection can forecast evolution in Streptococcus pneumoniae

Predicting how pathogen populations will change over time is challenging. Such has been the case with Streptococcus pneumoniae, an important human pathogen, and the pneumococcal conjugate vaccines (PCVs), which target only a fraction of the strains in the population. Here, we use the frequencies of accessory genes to predict changes in the pneumococcal population after vaccination, hypothesizing that these frequencies reflect negative frequency-dependent selection (NFDS) on the gene products. We find that the standardized predicted fitness of a strain, estimated by an NFDS-based model at the time the vaccine is introduced, enables us to predict whether the strain increases or decreases in prevalence following vaccination. Further, we are able to forecast the equilibrium post-vaccine population composition and assess the invasion capacity of emerging lineages. Overall, we provide a method for predicting the impact of an intervention on pneumococcal populations with potential application to other bacterial pathogens in which NFDS is a driving force.

Using affinity propagation clustering for identifying bacterial clades and subclades with whole-genome sequences of Francisella tularensis

By combining a reference-independent SNP analysis and average nucleotide identity (ANI) with affinity propagation clustering (APC), we developed a significantly improved methodology allowing resolving phylogenetic relationships, based on objective criteria. These bioinformatics tools can be used as a general ruler to determine phylogenetic relationships and clustering of bacteria, exemplary done with Francisella (F.) tularensis. Molecular epidemiology of F. tularensis is currently assessed mostly based on laboratory methods and molecular analysis. The high evolutionary stability and the clonal nature makes Francisella ideal for subtyping with single nucleotide polymorphisms (SNPs). Sequencing and real-time PCR can be used to validate the SNP analysis. We investigate whole-genome sequences of 155 F. tularensis subsp. holarctica isolates. Phylogenetic testing was based on SNPs and average nucleotide identity (ANI) as reference independent, alignment-free methods taking small-scale and large-scale differences within the genomes into account. Especially the whole genome SNP analysis with kSNP3.0 allowed deciphering quite subtle signals of systematic differences in molecular variation. Affinity propagation clustering (APC) resulted in three clusters showing the known clades B.4, B.6, and B.12. These data correlated with the results of real-time PCR assays targeting canSNPs loci. Additionally, we detected two subtle sub-clusters. SplitsTree was used with standard-setting using the aligned SNPs from Parsnps. Together APC, HierBAPS, and SplitsTree enabled us to generate hypotheses about epidemiologic relationships between bacterial clusters and describing the distribution of isolates. Our data indicate that the choice of the typing technique can increase our understanding of the pathogenesis and transmission of diseases with the eventual for prevention. This is opening perspectives to be applied to other bacterial species. The data provide evidence that Germany might be the collision zone where the clade B.12, also known as the East European clade, overlaps with the clade B.6, also known as the Iberian clade. Described methods allow generating a new, more detailed perspective for F. tularensis subsp. holarctica phylogeny. These results may encourage to determine phylogenetic relationships and clustering of other bacteria the same way.

First Report on the Finding of Listeria mnocytogenes ST121 Strain in a Dolphin Brain

Listeria monocytogenes (Lm) is a ubiquitous bacterium that causes the foodborne illness, listeriosis. Clonal complexes (CC), such as CC121, are overrepresented in the food production industry, and are rarely reported in animals and the environment. Working within a European-wide project, we investigated the routes by which strains are transmitted from environments and animals to food and the food production environment (FPE). In this context, we report, for the first time, the occurrence of a ST121 (CC121) strain isolated from a dolphin brain. The genome was compared with the genomes of 376 CC121 strains. Genomic comparisons showed that 16 strains isolated from food were the closest to the dolphin strain. Like most of the food strains analyzed here, the dolphin strain included genomic features (transposon Tn6188, plasmid pLM6179), both described as being associated with the strain’s adaptation to the FPE. Like all 376 strains, the dolphin strain contained a truncated actA gene and inlA gene, both described as being associated with attenuated virulence. Despite this fact, the strain was able to cross blood-brain barrier in immunosuppressed dolphin exposed polychlorinated biphenyl and invaded by parasites. Our data suggest that the dolphin was infected by a food-related strain released into the Mediterranean Sea.

Co-circulation of Multidrug-resistant Shigella Among Men Who Have Sex With Men in Australia

BACKGROUND

In urban Australia, the burden of shigellosis is either in returning travelers from shigellosis-endemic regions or in men who have sex with men (MSM). Here, we combine genomic data with comprehensive epidemiological data on sexual exposure and travel to describe the spread of multidrug-resistant Shigella lineages.

METHODS

A population-level study of all cultured Shigella isolates in the state of Victoria, Australia, was undertaken from 1 January 2016 through 31 March 2018. Antimicrobial susceptibility testing, whole-genome sequencing, and bioinformatic analyses of 545 Shigella isolates were performed at the Microbiological Diagnostic Unit Public Health Laboratory. Risk factor data on travel and sexual exposure were collected through enhanced surveillance forms or by interviews.

RESULTS

Rates of antimicrobial resistance were high, with 17.6% (95/541) and 50.6% (274/541) resistance to ciprofloxacin and azithromycin, respectively. There were strong associations between antimicrobial resistance, phylogeny, and epidemiology. Specifically, 2 major MSM-associated lineages were identified: a Shigellasonnei lineage (n = 159) and a Shigella flexneri 2a lineage (n = 105). Of concern, 147/159 (92.4%) of isolates within the S. sonnei MSM-associated lineage harbored mutations associated with reduced susceptibility to recommended oral antimicrobials: namely, azithromycin, trimethoprim-sulfamethoxazole, and ciprofloxacin. Long-read sequencing demonstrated global dissemination of multidrug-resistant plasmids across Shigella species and lineages, but predominantly associated with MSM isolates.

CONCLUSIONS

Our contemporary data highlight the ongoing public health threat posed by resistant Shigella, both in Australia and globally. Urgent multidisciplinary public health measures are required to interrupt transmission and prevent infection.

Twentieth-century emergence of antimicrobial resistant human- and bovine-associated Salmonella enterica serotype Typhimurium lineages in New York State

Salmonella enterica serotype Typhimurium (S. Typhimurium) boasts a broad host range and can be transmitted between livestock and humans. While members of this serotype can acquire resistance to antimicrobials, the temporal dynamics of this acquisition is not well understood. Using New York State (NYS) and its dairy cattle farms as a model system, 87 S. Typhimurium strains isolated from 1999 to 2016 from either human clinical or bovine-associated sources in NYS were characterized using whole-genome sequencing. More than 91% of isolates were classified into one of four major lineages, two of which were largely susceptible to antimicrobials but showed sporadic antimicrobial resistance (AMR) gene acquisition, and two that were largely multidrug-resistant (MDR). All four lineages clustered by presence and absence of elements in the pan-genome. The two MDR lineages, one of which resembled S. Typhimurium DT104, were predicted to have emerged circa 1960 and 1972. The two largely susceptible lineages emerged earlier, but showcased sporadic AMR determinant acquisition largely after 1960, including acquisition of cephalosporin resistance-conferring genes after 1985. These results confine the majority of AMR acquisition events in NYS S. Typhimurium to the twentieth century, largely within the era of antibiotic usage.

Diversification in immunogenicity genes caused by selective pressures in invasive meningococci

We studied population genomics of 486 Neisseria meningitidis isolates causing meningitis in the Netherlands during the period 1979-2003 and 2006-2013 using whole-genome sequencing to evaluate the impact of a hyperendemic period of serogroup B invasive disease. The majority of serogroup B isolates belonged to ST-41/44 (41 %) and ST-32 complex (16 %). Comparing the time periods, before and after the decline of serogroup B invasive disease, there was a decrease of ST-41/44 complex sequences (P=0.002). We observed the expansion of a sub-lineage within ST-41/44 complex sequences being associated with isolation from the 1979-2003 time period (P=0.014). Isolates belonging to this sub-lineage expansion within ST-41/44 complex were marked by four antigen allele variants. Presence of these allele variants was associated with isolation from the 1979-2003 time period after correction for multiple testing (Wald test, P=0.0043 for FetA 1-5; P=0.0035 for FHbp 14; P=0.012 for PorA 7-2.4 and P=0.0031 for NHBA two peptide allele). These sequences were associated with 4CMenB vaccine coverage (Fisher's exact test, P<0.001). Outside of the sub-lineage expansion, isolates with markedly lower levels of predicted vaccine coverage clustered in phylogenetic groups showing a trend towards isolation in the 2006-2013 time period (P=0.08). In conclusion, we show the emergence and decline of a sub-lineage expansion within ST-41/44 complex isolates concurrent with a hyperendemic period in meningococcal meningitis. The expansion was marked by specific antigen peptide allele combinations. We observed preliminary evidence for decreasing 4CMenB vaccine coverage in the post-hyperendemic period.

Strong Gene Flow Undermines Local Adaptations in a Host Parasite System

Simple Summary

The co-evolution of hosts and parasites depends on their ability to adapt to each other’s defense and counter-defense mechanisms. The strength of selection on those mechanisms may vary among populations, resulting in a geographical mosaic of co-evolution. The boreo-montane paper wasp Polistes biglumis and its parasite Polistes atrimandibularis exemplify this type of co-evolutionary system. Here, we used genetic markers to examine the genetic population structures of these wasps in the western Alps. We found that both host and parasite populations displayed similar levels of genetic variation. In the host species, populations located near to each other were genetically similar; in both the host and the parasite species populations farther apart were significantly different. Thus, apparent dispersal barriers (i.e., high mountains) did not seem to restrict gene flow across populations as expected. Furthermore, there were no major differences in gene flow between the two species, perhaps because P. atrimandibularis parasitizes both alpine and lowland host species and annually migrates between alpine and lowland populations. The presence of strong gene flow in a system where local populations experience variable levels of selection pressure challenges the classical hypothesis that restricted gene flow is required for local adaptations to evolve.

Abstract

The co-evolutionary pathways followed by hosts and parasites strongly depend on the adaptive potential of antagonists and its underlying genetic architecture. Geographically structured populations of interacting species often experience local differences in the strength of reciprocal selection pressures, which can result in a geographic mosaic of co-evolution. One example of such a system is the boreo-montane social wasp Polistes biglumis and its social parasite Polistes atrimandibularis, which have evolved local defense and counter-defense mechanisms to match their antagonist. In this work, we study spatial genetic structure of P. biglumis and P. atrimandibularis populations at local and regional scales in the Alps, by using nuclear markers (DNA microsatellites, AFLP) and mitochondrial sequences. Both the host and the parasite populations harbored similar amounts of genetic variation. Host populations were not genetically structured at the local scale, but geographic regions were significantly differentiated from each other in both the host and the parasite in all markers. The net dispersal inferred from genetic differentiation was similar in the host and the parasite, which may be due to the annual migration pattern of the parasites between alpine and lowland populations. Thus, the apparent dispersal barriers (i.e., high mountains) do not restrict gene flow as expected and there are no important gene flow differences between the species, which contradict the hypothesis that restricted gene flow is required for local adaptations to evolve.

Mass drug administration with azithromycin for trachoma elimination and the population structure of Streptococcus pneumoniae in the nasopharynx

Objective

Mass drug administration (MDA) with azithromycin for trachoma elimination reduces nasopharyngeal carriage of Streptococcus pneumoniae in the short term. We evaluated S. pneumoniae carried in the nasopharynx before and after a round of azithromycin MDA to determine whether MDA was associated with changes in pneumococcal population structure and resistance.

Methods

We analysed 514 pneumococcal whole genomes randomly selected from nasopharyngeal samples collected in two Gambian villages that received three annual rounds of MDA for trachoma elimination. The 514 samples represented 293 participants, of which 75% were children aged 0–9 years, isolated during three cross-sectional surveys (CSSs) conducted before the third round of MDA (CSS-1) and at 1 (CSS-2) and 6 (CSS-3) months after MDA. Bayesian Analysis of Population Structure (BAPS) was used to cluster related isolates by capturing variation in the core genome. Serotype and multilocus sequence type were inferred from the genotype. Antimicrobial resistance determinants were identified from assemblies, including known macrolide resistance genes.

Results

Twenty-seven BAPS clusters were assigned. These consisted of 81 sequence types (STs). Two BAPS clusters not observed in CSS-1 (n = 109) or CSS-2 (n = 69), increased in frequency in CSS-3 (n = 126); BAPS20 (8.73%, p 0.016) and BAPS22 (7.14%, p 0.032) but were not associated with antimicrobial resistance. Macrolide resistance within BAPS17 increased after treatment (CSS-1 n = 0/6, CSS-2/3 n = 5/5, p 0.002) and was carried on a mobile transposable element that also conferred resistance to tetracycline.

Discussion

Limited changes in pneumococcal population structure were observed after the third round of MDA, suggesting treatment had little effect on the circulating lineages. An increase in macrolide resistance within one BAPS highlights the need for antimicrobial resistance surveillance in treated villages.

Whole-Genome-Based Survey for Polyphyletic Serovars of Salmonella enterica subsp. enterica Provides New Insights into Public Health Surveillance

Serotyping has traditionally been considered the basis for surveillance of Salmonella, but it cannot distinguish distinct lineages sharing the same serovar that vary in host range, pathogenicity and epidemiology. However, polyphyletic serovars have not been extensively investigated. Public health microbiology is currently being transformed by whole-genome sequencing (WGS) data, which promote the lineage determination using a more powerful and accurate technique than serotyping. The focus in this study is to survey and analyze putative polyphyletic serovars. The multi-locus sequence typing (MLST) phylogenetic analysis identified four putative polyphyletic serovars, namely, Montevideo, Bareilly, Saintpaul, and Muenchen. Whole-genome-based phylogeny and population structure highlighted the polyphyletic nature of Bareilly and Saintpaul and the multi-lineage nature of Montevideo and Muenchen. The population of these serovars was defined by extensive genetic diversity, the open pan genome and the small core genome. Source niche metadata revealed putative existence of lineage-specific niche adaptation (host-preference and environmental-preference), exhibited by lineage-specific genomic contents associated with metabolism and transport. Meanwhile, differences in genetic profiles relating to virulence and antimicrobial resistance within each lineage may contribute to pathogenicity and epidemiology. The results also showed that recombination events occurring at the H1-antigen loci may be an important reason for polyphyly. The results presented here provide the genomic basis of simple, rapid, and accurate identification of phylogenetic lineages of these serovars, which could have important implications for public health.

Chromosome Architecture and Gene Content of the Emergent Pathogen Acinetobacter haemolyticus

Acinetobacter haemolyticus is a Gammaproteobacterium that has been involved in serious diseases frequently linked to the nosocomial environment. Most of the strains causing such infections are sensitive to a wide variety of antibiotics, but recent reports indicate that this pathogen is acquiring very efficiently carbapenem-resistance determinants like the blaNDM-1 gene, all over the world. With this work we contribute with a collection set of 31 newly sequenced nosocomial A. haemolyticus isolates. Genome analysis of these sequences and others collected from RefSeq indicates that their chromosomes are organized in 12 syntenic blocks that contain most of the core genome genes. These blocks are separated by hypervariable regions that are rich in unique gene families, but also have signals of horizontal gene transfer. Genes involved in virulence or encoding different secretion systems are located inside syntenic regions and have recombination signals. The relative order of the synthetic blocks along the A. haemolyticus chromosome can change, indicating that they have been subject to several kinds of inversions. Genomes of this microorganism show large differences in gene content even if they are in the same clade. Here we also show that A. haemolyticus has an open pan-genome.

High-resolution characterisation of ESBL/pAmpC-producing Escherichia coli isolated from the broiler production pyramid

The presence of extended-spectrum β-lactamase (ESBL) or plasmid-mediated AmpC β-lactamase (pAmpC)-producing Escherichia coli (ESBL/pAmpC-EC) in livestock is a public health risk given the likelihood of their transmission to humans via the food chain. We conducted whole genome sequencing on 100 ESBL/pAmpC-EC isolated from the broiler production to explore their resistance and virulence gene repertoire, characterise their plasmids and identify transmission events derived from their phylogeny. Sequenced isolates carried resistance genes to four antimicrobial classes in addition to cephalosporins. Virulence gene analysis assigned the majority of ESBL/pAmpC-EC to defined pathotypes. In the complex genetic background of ESBL/pAmpC-EC, clusters of closely related isolates from various production stages were identified and indicated clonal transmission. Phylogenetic comparison with publicly available genomes suggested that previously uncommon ESBL/pAmpC-EC lineages could emerge in poultry, while others might contribute to the maintenance and dissemination of ESBL/pAmpC genes in broilers. The majority of isolates from diverse E. coli lineages shared four dominant plasmids (IncK2, IncI1, IncX3 and IncFIB/FII) with identical ESBL/pAmpC gene insertion sites. These plasmids have been previously reported in diverse hosts, including humans. Our findings underline the importance of specific plasmid groups in the dissemination of cephalosporin resistance genes within the broiler industry and across different reservoirs.

Phylogeographical Analyses and Antibiotic Resistance Genes of Acinetobacter johnsonii Highlight Its Clinical Relevance

Acinetobacter johnsonii has been severely understudied and its population structure and the presence of antibiotic resistance genes (ARGs) are very much uncertain. Our phylogeographical analysis shows that intercontinental transmission has occurred frequently and that different lineages are circulating within single countries; notably, clinical and nonclinical strains are not well differentiated from one another. Importantly, in this species recombination is a significant source of single nucleotide polymorphisms.

Acinetobacter johnsonii has been severely understudied and its population structure and the presence of antibiotic resistance genes (ARGs) are very much uncertain. Our phylogeographical analysis shows that intercontinental transmission has occurred frequently and that different lineages are circulating within single countries; notably, clinical and nonclinical strains are not well differentiated from one another. Importantly, in this species recombination is a significant source of single nucleotide polymorphisms. Furthermore, our results show this species could be an important reservoir of ARGs since it has a significant amount of ARGs, and many of them show signals of horizontal gene transfer. Thus, this study clearly points out the clinical importance of A. johnsonii and the urgent need to better appreciate its genomic diversity.

Unexpected conservation and global transmission of agrobacterial virulence plasmids

The accelerated evolution and spread of pathogens are threats to host species. Agrobacteria require an oncogenic Ti or Ri plasmid to transfer genes into plants and cause disease. We developed a strategy to characterize virulence plasmids and applied it to analyze hundreds of strains collected between 1927 and 2017, on six continents and from more than 50 host species. In consideration of prior evidence for prolific recombination, it was surprising that oncogenic plasmids are descended from a few conserved lineages. Characterization of a hierarchy of features that promote or constrain plasticity allowed inference of the evolutionary history across the plasmid lineages. We uncovered epidemiological patterns that highlight the importance of plasmid transmission in pathogen diversification as well as in long-term persistence and the global spread of disease.

Increasing incidence of group B streptococcus neonatal infections in the Netherlands is associated with clonal expansion of CC17 and CC23

Group B streptococcus (GBS) is the leading cause of neonatal invasive disease worldwide. In the Netherlands incidence of the disease increased despite implementation of preventive guidelines. We describe a genomic analysis of 1345 GBS isolates from neonatal (age 0–89 days) invasive infections in the Netherlands reported between 1987 and 2016. Most isolates clustered into one of five major lineages: CC17 (39%), CC19 (25%), CC23 (18%), CC10 (9%) and CC1 (7%). There was a significant rise in the number of infections due to isolates from CC17 and CC23. Phylogenetic clustering analysis revealed that this was caused by expansion of specific sub-lineages, designated CC17-A1, CC17-A2 and CC23-A1. Dating of phylogenetic trees estimated that these clones diverged in the 1960s/1970s, representing historical rather than recently emerged clones. For CC17-A1 the expansion correlated with acquisition of a new phage, carrying gene encoding a putative cell-surface protein. Representatives of CC17-A1, CC17-A2 and CC23-A1 clones were identified in datasets from other countries demonstrating their global distribution.

Evolution of Salmonella enterica serotype Typhimurium driven by anthropogenic selection and niche adaptation

Salmonella enterica serotype Typhimurium (S. Typhimurium) is a leading cause of gastroenteritis and bacteraemia worldwide, and a model organism for the study of host-pathogen interactions. Two S. Typhimurium strains (SL1344 and ATCC14028) are widely used to study host-pathogen interactions, yet genotypic variation results in strains with diverse host range, pathogenicity and risk to food safety. The population structure of diverse strains of S. Typhimurium revealed a major phylogroup of predominantly sequence type 19 (ST19) and a minor phylogroup of ST36. The major phylogroup had a population structure with two high order clades (α and β) and multiple subclades on extended internal branches, that exhibited distinct signatures of host adaptation and anthropogenic selection. Clade α contained a number of subclades composed of strains from well characterized epidemics in domesticated animals, while clade β contained multiple subclades associated with wild avian species. The contrasting epidemiology of strains in clade α and β was reflected by the distinct distribution of antimicrobial resistance (AMR) genes, accumulation of hypothetically disrupted coding sequences (HDCS), and signatures of functional diversification. These observations were consistent with elevated anthropogenic selection of clade α lineages from adaptation to circulation in populations of domesticated livestock, and the predisposition of clade β lineages to undergo adaptation to an invasive lifestyle by a process of convergent evolution with of host adapted Salmonella serotypes. Gene flux was predominantly driven by acquisition and recombination of prophage and associated cargo genes, with only occasional loss of these elements. The acquisition of large chromosomally-encoded genetic islands was limited, but notably, a feature of two recent pandemic clones (DT104 and monophasic S. Typhimurium ST34) of clade α (SGI-1 and SGI-4).

Contaminated in-house environment contributes to the persistence and transmission of NDM-producing bacteria in a Chinese poultry farm

While carbapenem use is prohibited in the poultry production chain and carbapenem-resistant Enterobacteriaceae (CRE) are absent from hatchery farms, New Delhi metallo-β-lactamase-producing CRE contamination of commercial broiler chicken farms (grow-out farms) can occur via living hosts such as flies. However, it is not known whether the inanimate factors from in-house environment play a role in the persistence of CRE on commercial farms. Herein, we monitored one typical broiler house in Hebei Province, China, from January 2017 to April 2018. We collected 350 cloacal samples from four broiler batches along with 582 environmental samples (194 in the raising period and 388 in the vacancy period) from sites including the surfaces of drooping boards, feeding troughs, nipple drinkers, corridor floors, sewage trenches, and air. All samples were screened for bla_NDM and cultured for NDM-producing isolates. The resistance profiles, genotypes, and genetic context of bla_NDM in CRE isolates were further characterized. Results showed that 1-day-old broilers, which were transferred from a hatchery farm and negative for CRE, acquired bla_NDM within 24 h of transfer (2 days of age), with a detection rate of up to 18.6%. High bla_NDM detection rates (26.8%-31.4%) were obtained among all environmental samples except air after standard cleaning and disinfection during the vacancy period. bla_NDM carriage rates (52.9%-72.9%) within the flocks remain stable and high across the next three broiler batches. Overall, 279 NDM-producing bacteria, including 259 Enterobacteriaceae (8 species), 14 Morganellaceae (3 species), three Alcaligenes faecalis and three Pseudomonas putida isolates, were recovered from 85 (24.3%) cloacal and 101 (17.4%) environmental samples. Three NDM variants, including NDM-5 (n = 181), NDM-1 (n = 92), and NDM-9 (n = 3), and a novel NDM-like-metallo-β-lactamase (NLM, n = 3) were identified among the samples. The predominant NDM-producing CRE species among the samples were Klebsiella pneumoniae (CRKP; 32.6%, n = 91) and Escherichia coli (CREC; 27.2%, n = 76). Both clonal and horizontal transmission of bla_NDM and an overlap of sequence types (STs) were observed in both CREC and CRKP from chicken and environmental samples. Notably, ST6751 CREC and ST37 CRKP persisted throughout the 16-month surveillance period. IncX3 (n = 197, 7 species), IncA/C₂ (n = 41, 5 species), and IncFII (n = 8, E. coli) were the three major bla_NDM-carrying plasmid types among the isolates. Although routine cleaning and disinfection procedures and "all-in/all-out" management were performed, once introduced to the farm environment, a diverse range of NDM-positive isolates may survive and persist, becoming an important reservoir of NDM-positive CRE for broiler chickens. Therefore, cleaning and disinfection procedures should be improved on poultry farms to avoid cross-contamination of NDM-producing bacteria between different batches of chickens, as well as further downstream in the poultry production chain.