Hierarchical and spatially explicit clustering of DNA sequences with BAPS software

Phylogenomics of Enterococcus faecalis from wild birds: new insights into host-associated differences in core and accessory genomes of the species

Wild birds have been suggested to be reservoirs of antimicrobial resistant and/or pathogenic Enterococcus faecalis (Efs) strains, but the scarcity of studies and available sequences limit our understanding of the population structure of the species in these hosts. Here, we analysed the clonal and plasmid diversity of 97 Efs isolates from wild migratory birds. We found a high diversity, with most sequence types (STs) being firstly described here, while others were found in other hosts including some predominant in poultry. We found that pheromone-responsive plasmids predominate in wild bird Efs while 35% of the isolates entirely lack plasmids. Then, to better understand the ecology of the species, the whole genome of fivestrains with known STs (ST82, ST170, ST16 and ST55) were sequenced and compared with all the Efs genomes available in public databases. Using several methods to analyse core and accessory genomes (AccNET, PLACNET, hierBAPS and PANINI), we detected differences in the accessory genome of some lineages (e.g. ST82) demonstrating specific associations with birds. Conversely, the genomes of other Efs lineages exhibited divergence in core and accessory genomes, reflecting different adaptive trajectories in various hosts. This pangenome divergence, horizontal gene transfer events and occasional epidemic peaks could explain the population structure of the species.

The Population Structure of Pseudomonas aeruginosa Is Characterized by Genetic Isolation of exoU+ and exoS+ Lineages

The diversification of microbial populations may be driven by many factors including adaptation to distinct ecological niches and barriers to recombination. We examined the population structure of the bacterial pathogen Pseudomonas aeruginosa by analyzing whole-genome sequences of 739 isolates from diverse sources. We confirmed that the population structure of P. aeruginosa consists of two major groups (referred to as Groups A and B) and at least two minor groups (Groups C1 and C2). Evidence for frequent intragroup but limited intergroup recombination in the core genome was observed, consistent with sexual isolation of the groups. Likewise, accessory genome analysis demonstrated more gene flow within Groups A and B than between these groups, and a few accessory genomic elements were nearly specific to one or the other group. In particular, the exoS gene was highly overrepresented in Group A compared with Group B isolates (99.4% vs. 1.1%) and the exoU gene was highly overrepresented in Group B compared with Group A isolates (95.2% vs. 1.8%). The exoS and exoU genes encode effector proteins secreted by the P. aeruginosa type III secretion system. Together these results suggest that the major P. aeruginosa groups defined in part by the exoS and exoU genes are divergent from each other, and that these groups are genetically isolated and may be ecologically distinct. Although both groups were globally distributed and caused human infections, certain groups predominated in some clinical contexts.

Species delimitation and evolutionary reconstruction within an integrative taxonomic framework: A case study on Rhinolophus macrotis complex (Chiroptera: Rhinolophidae)

Species delimitation and evolutionary reconstruction remain challenging for non-model species that have experienced reticulate evolution and exhibit conflicting patterns of differentiation among multiple lines of evidence, such as mitochondrial and nuclear data and phenotypes. Here, we applied an integrative taxonomic approach to a case study of Rhinolophus macrotis complex, whose taxonomic status remains controversial, to provide insight into the systematics and evolutionary history of these species. By integrating traditional genetic markers with different modes of inheritance, genome-wide SNPs as well as phenotypic characteristics, we clarified the presence of three closely related species, R. episcopus, R. siamensis, and R. osgoodi, within this complex, and proposed a new taxonomic treatment for R. osgoodi. Our results suggested that hybridization and introgression are the main causes of low mtDNA divergence in these species. Combined with the demographic inference, we deduced that glacial-interglacial cycles drove geographic isolation and secondary contacts of these species, then promoted hybridization and lineage fusion among them, finally resulting in a reticulate evolutionary pattern. Overall, our study highlights the importance of combining multiple types of data to delimit species, especially those with conserved morphology, and to reveal the sophisticated processes of speciation.

Fast hierarchical Bayesian analysis of population structure

We present fastbaps, a fast solution to the genetic clustering problem. Fastbaps rapidly identifies an approximate fit to a Dirichlet process mixture model (DPM) for clustering multilocus genotype data. Our efficient model-based clustering approach is able to cluster datasets 10-100 times larger than the existing model-based methods, which we demonstrate by analyzing an alignment of over 110 000 sequences of HIV-1 pol genes. We also provide a method for rapidly partitioning an existing hierarchy in order to maximize the DPM model marginal likelihood, allowing us to split phylogenetic trees into clades and subclades using a population genomic model. Extensive tests on simulated data as well as a diverse set of real bacterial and viral datasets show that fastbaps provides comparable or improved solutions to previous model-based methods, while being significantly faster. The method is made freely available under an open source MIT licence as an easy to use R package at https://github.com/gtonkinhill/fastbaps.

Long-Term Intrahost Evolution of Methicillin Resistant Staphylococcus aureus Among Cystic Fibrosis Patients With Respiratory Carriage

Staphylococcus aureus is the most commonly identified airway colonizer of cystic fibrosis (CF) patients, and infections with methicillin-resistant S. aureus (MRSA) are associated with poor outcomes. Yet, little is known about the intrahost evolution of S. aureus among CF patients. We investigated convergent evolution and adaptation of MRSA among four CF patients with long-term respiratory carriage. For each patient, we performed whole-genome sequencing on an average of 21 isolates (range: 19–23) carried for a mean of 1,403 days (range: 903–1,679), including 25 pairs of isolates collected on the same day. We assessed intrahost diversity, population structure, evolutionary history, evidence of switched intergenic regions (IGRs), and signatures of adaptation in the context of patient age, antibiotic treatment, and co-colonizing microbes. Phylogenetic analysis delineated distinct multilocus sequence type ST5 (n = 3) and ST72 (n = 1) clonal populations in addition to sporadic, non-clonal isolates, and uncovered a putative transmission event. Variation in antibiotic resistance was observed within clonal populations, even among isolates collected on the same day. Rates of molecular evolution ranged from 2.21 to 8.64 nucleotide polymorphisms per year, and lineage ages were consistent with acquisition of colonization in early childhood followed by subsequent persistence of multiple sub-populations. Selection analysis of 1,622 core genes present in all four clonal populations (n = 79) found 11 genes variable in three subjects – most notably, ATP-dependent protease clpX, 2-oxoglutarate dehydrogenase odhA, fmtC, and transcription-repair coupling factor mfd. Only one gene, staphylococcal protein A (spa), was found to have evidence of gene-wide diversifying selection. We identified three instances of intrahost IGR switching events, two of which flanked genes related to quorum sensing. The complex microbial ecology of the CF airway poses challenges for management. We illustrate appreciable intrahost diversity as well as persistence of a dominant lineage. We also show that intrahost adaptation is a continual process, despite purifying selective pressure, and provide targets that should be investigated further for their function in CF adaptation.

Genomic and Functional Analysis of Emerging Virulent and Multidrug-Resistant Escherichia coli Lineage Sequence Type 648

The pathogenic extended-spectrum-beta-lactamase (ESBL)-producing Escherichia coli lineage ST648 is increasingly reported from multiple origins. Our study of a large and global ST648 collection from various hosts (87 whole-genome sequences) combining core and accessory genomics with functional analyses and in vivo experiments suggests that ST648 is a nascent and generalist lineage, lacking clear phylogeographic and host association signals. By including large numbers of ST131 (n = 107) and ST10 (n = 96) strains for comparative genomics and phenotypic analysis, we demonstrate that the combination of multidrug resistance and high-level virulence are the hallmarks of ST648, similar to international high-risk clonal lineage ST131. Specifically, our in silico, in vitro, and in vivo results demonstrate that ST648 is well equipped with biofilm-associated features, while ST131 shows sophisticated signatures indicative of adaption to urinary tract infection, potentially conveying individual ecological niche adaptation. In addition, we used a recently developed NFDS (negative frequency-dependent selection) population model suggesting that ST648 will increase significantly in frequency as a cause of bacteremia within the next few years. Also, ESBL plasmids impacting biofilm formation aided in shaping and maintaining ST648 strains to successfully emerge worldwide across different ecologies. Our study contributes to understanding what factors drive the evolution and spread of emerging international high-risk clonal lineages.

Joint sequencing of human and pathogen genomes reveals the genetics of pneumococcal meningitis

Streptococcus pneumoniae is a common nasopharyngeal colonizer, but can also cause life-threatening invasive diseases such as empyema, bacteremia and meningitis. Genetic variation of host and pathogen is known to play a role in invasive pneumococcal disease, though to what extent is unknown. In a genome-wide association study of human and pathogen we show that human variation explains almost half of variation in susceptibility to pneumococcal meningitis and one-third of variation in severity, identifying variants in CCDC33 associated with susceptibility. Pneumococcal genetic variation explains a large amount of invasive potential (70%), but has no effect on severity. Serotype alone is insufficient to explain invasiveness, suggesting other pneumococcal factors are involved in progression to invasive disease. We identify pneumococcal genes involved in invasiveness including pspC and zmpD, and perform a human-bacteria interaction analysis. These genes are potential candidates for the development of more broadly-acting pneumococcal vaccines.

International genomic definition of pneumococcal lineages, to contextualise disease, antibiotic resistance and vaccine impact

BACKGROUND

Pneumococcal conjugate vaccines have reduced the incidence of invasive pneumococcal disease, caused by vaccine serotypes, but non-vaccine-serotypes remain a concern. We used whole genome sequencing to study pneumococcal serotype, antibiotic resistance and invasiveness, in the context of genetic background.

METHODS

Our dataset of 13,454 genomes, combined with four published genomic datasets, represented Africa (40%), Asia (25%), Europe (19%), North America (12%), and South America (5%). These 20,027 pneumococcal genomes were clustered into lineages using PopPUNK, and named Global Pneumococcal Sequence Clusters (GPSCs). From our dataset, we additionally derived serotype and sequence type, and predicted antibiotic sensitivity. We then measured invasiveness using odds ratios that relating prevalence in invasive pneumococcal disease to carriage.

FINDINGS

The combined collections (n = 20,027) were clustered into 621 GPSCs. Thirty-five GPSCs observed in our dataset were represented by >100 isolates, and subsequently classed as dominant-GPSCs. In 22/35 (63%) of dominant-GPSCs both non-vaccine serotypes and vaccine serotypes were observed in the years up until, and including, the first year of pneumococcal conjugate vaccine introduction. Penicillin and multidrug resistance were higher (p < .05) in a subset dominant-GPSCs (14/35, 9/35 respectively), and resistance to an increasing number of antibiotic classes was associated with increased recombination (R2 = 0.27 p < .0001). In 28/35 dominant-GPSCs, the country of isolation was a significant predictor (p < .05) of its antibiogram (mean misclassification error 0.28, SD ± 0.13). We detected increased invasiveness of six genetic backgrounds, when compared to other genetic backgrounds expressing the same serotype. Up to 1.6-fold changes in invasiveness odds ratio were observed.

INTERPRETATION

We define GPSCs that can be assigned to any pneumococcal genomic dataset, to aid international comparisons. Existing non-vaccine-serotypes in most GPSCs preclude the removal of these lineages by pneumococcal conjugate vaccines; leaving potential for serotype replacement. A subset of GPSCs have increased resistance, and/or serotype-independent invasiveness.

Phylogeographic structure and ecological niche modelling reveal signals of isolation and postglacial colonisation in the European stag beetle

Lucanus cervus (L.), the stag beetle, is a saproxylic beetle species distributed widely across Europe. Throughout its distribution the species has exhibited pronounced declines and is widely considered threatened. Conservation efforts may be hindered by the lack of population genetic data and understanding of the spatial scale of population connectivity. To address this knowledge gap this research details the first broad scale phylogeographic study of L. cervus based on mitochondrial DNA (mtDNA) sequencing and microsatellite analysis of samples collected from 121 localities across Europe. Genetic data were complemented by palaeo-distribution models of spatial occupancy during the Last Glacial Maximum to strengthen inferences of refugial areas. A salient feature of the mtDNA was the identification of two lineages. Lineage I was widespread across Europe while lineage II was confined to Greece. Microsatellites supported the differentiation of the Greek samples and alongside palaeo-distribution models indicated this area was a glacial refuge. The genetic endemism of the Greek samples, and demographic results compatible with no signatures of spatial expansion likely reflects restricted dispersal into and out of the area. Lineage I exhibited a shallow star like phylogeny compatible with rapid population expansion across Europe. Demographic analysis indicated such expansions occurred after the Last Glacial Maximum. Nuclear diversity and hindcast species distribution models indicated a central Italian refuge for lineage I. Palaeo-distribution modelling results also suggested a western refuge in northern Iberia and south-west France. In conclusion the results provide evidence of glacial divergence in stag beetle while also suggesting high, at least on evolutionary timescales, gene flow across most of Europe. The data also provide a neutral genetic framework against which patterns of phenotypic variation may be assessed.

Comparative genomic analysis of Bacillus paralicheniformis MDJK30 with its closely related species reveals an evolutionary relationship between B. paralicheniformis and B. licheniformis

Background

Members of the genus Bacillus are important plant growth-promoting rhizobacteria that serve as biocontrol agents. Bacillus paralicheniformis MDJK30 is a PGPR isolated from the peony rhizosphere and can suppress plant-pathogenic bacteria and fungi. To further uncover the genetic mechanism of the plant growth-promoting traits of MDJK30 and its closely related strains, we used comparative genomics to provide insights into the genetic diversity and evolutionary relationship between B. paralicheniformis and B. licheniformis.

Results

A comparative genomics analysis based on B. paralicheniformis MDJK30 and 55 other previously reported Bacillus strains was performed. The evolutionary position of MDJK30 and the evolutionary relationship between B. paralicheniformis and B. licheniformis were evaluated by studying the phylogeny of the core genomes, a population structure analysis and ANI results. Comparative genomic analysis revealed various features of B. paralicheniformis that contribute to its commensal lifestyle in the rhizosphere, including an opening pan genome, a diversity of transport and the metabolism of the carbohydrates and amino acids. There are notable differences in the numbers and locations of the insertion sequences, prophages, genomic islands and secondary metabolic synthase operons between B. paralicheniformis and B. licheniformis. In particular, we found most gene clusters of Fengycin, Bacitracin and Lantipeptide were only present in B. paralicheniformis and were obtained by horizontal gene transfer (HGT), and these clusters may be used as genetic markers for distinguishing B. paralicheniformis and B. licheniformis.

Conclusions

This study reveals that MDJK30 and the other strains of lineage paralicheniformis present plant growth-promoting traits at the genetic level and can be developed and commercially formulated in agriculture as PGPR. Core genome phylogenies and population structure analysis has proven to be a powerful tool for differentiating B. paralicheniformis and B. licheniformis. Comparative genomic analyses illustrate the genetic differences between the paralicheniformis-licheniformis group with respect to rhizosphere adaptation.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5646-9) contains supplementary material, which is available to authorized users.

The history, genome and biology of NCTC 30: a non-pandemic Vibrio cholerae isolate from World War One

The sixth global cholera pandemic lasted from 1899 to 1923. However, despite widespread fear of the disease and of its negative effects on troop morale, very few soldiers in the British Expeditionary Forces contracted cholera between 1914 and 1918. Here, we have revived and sequenced the genome of NCTC 30, a 102-year-old Vibrio cholerae isolate, which we believe is the oldest publicly available live V. cholerae strain in existence. NCTC 30 was isolated in 1916 from a British soldier convalescent in Egypt. We found that this strain does not encode cholera toxin, thought to be necessary to cause cholera, and is not part of V. cholerae lineages responsible for the pandemic disease. We also show that NCTC 30, which predates the introduction of penicillin-based antibiotics, harbours a functional β-lactamase antibiotic resistance gene. Our data corroborate and provide molecular explanations for previous phenotypic studies of NCTC 30 and provide a new high-quality genome sequence for historical, non-pandemic V. cholerae.

Genomic Epidemiology of Global Carbapenemase-Producing Enterobacter spp., 2008-2014

We performed whole-genome sequencing on 170 clinical carbapenemase-producing Enterobacter spp. isolates collected globally during 2008-2014. The most common carbapenemase was VIM, followed by New Delhi metallo-β-lactamase (NDM), Klebsiella pneumoniae carbapenemase, oxacillin 48, and IMP. The isolates were of predominantly 2 species (E. xiangfangensis and E. hormaechei subsp. steigerwaltii) and 4 global clones (sequence type [ST] 114, ST93, ST90, and ST78) with different clades within ST114 and ST90. Particular genetic structures surrounding carbapenemase genes were circulating locally in various institutions within the same or between different STs in Greece, Guatemala, Italy, Spain, Serbia, and Vietnam. We found a common NDM genetic structure (NDM-GE-U.S.), previously described on pNDM-U.S. from Klebsiella pneumoniae ATCC BAA-214, in 14 different clones obtained from 6 countries spanning 4 continents. Our study highlights the importance of surveillance programs using whole-genome sequencing in providing insight into the molecular epidemiology of carbapenemase-producing Enterobacter spp.

The role of glacial-interglacial climate change in shaping the genetic structure of eastern subterranean termites in the southern Appalachian Mountains, USA

The eastern subterranean termite, Reticulitermes flavipes, currently inhabits previously glaciated regions of the northeastern U.S., as well as the unglaciated southern Appalachian Mountains and surrounding areas. We hypothesized that Pleistocene climatic fluctuations have influenced the distribution of R. flavipes, and thus the evolutionary history of the species. We estimated contemporary and historical geographic distributions of R. flavipes by constructing Species Distribution Models (SDM). We also inferred the evolutionary and demographic history of the species using mitochondrial (cytochrome oxidase I and II) and nuclear (endo-beta-1,4-glucanase) DNA sequence data. To do this, genetic populations were delineated using Bayesian spatial-genetic clustering, competing hypotheses about population divergence were assessed using approximate Bayesian computation (ABC), and changes in population size were estimated using Bayesian skyline plots. SDMs identified areas in the north with suitable habitat during the transition from the Last Interglacial to the Last Glacial Maximum, as well as an expanding distribution from the mid-Holocene to the present. Genetic analyses identified three geographically cohesive populations, corresponding with northern, central, and southern portions of the study region. Based on ABC analyses, divergence between the Northern and Southern populations was the oldest, estimated to have occurred 64.80 thousand years ago (kya), which corresponds with the timing of available habitat in the north. The Central and Northern populations diverged in the mid-Holocene, 8.63 kya, after which the Central population continued to expand. Accordingly, phylogeographic patterns of R. flavipes in the southern Appalachians appear to have been strongly influenced by glacial-interglacial climate change.

OPEN RESEARCH BADGES

This article has been awarded Open Materials, Open Data Badges. All materials and data are publicly accessible via the Open Science Framework at https://doi.org/10.5061/dryad.5hr7f31.

Phylogeographic diversification and postglacial range dynamics shed light on the conservation of the kelp Saccharina japonica

Studies of postglacial range shifts could enhance our understanding of seaweed species' responses to climate change and hence facilitate the conservation of natural resources. However, the distribution dynamics and phylogeographic diversification of the commercially and ecologically important kelp Saccharina japonica in the Northwest Pacific (NWP) are still poorly surveyed. In this study, we analyzed the evolutionary history of S. japonica using two mitochondrial markers and 24 nuclear microsatellites. A STRUCTURE analysis revealed two partially isolated lineages: lineage H, which is scattered along the coast of Japan; and lineage P, which occurs along the west coast of the Japan Sea. Ecological niche modeling projections to the Last Glacial Maximum (LGM) revealed that the southern coasts of the Japan Sea and the Pacific side of the Oshima and Honshu Peninsulas provided the most suitable habitats for S. japonica, implying that these regions served as ancient refugia during the LGM. Ancient isolation in different refugia may explain the observed divergence between lineages P and H. An approximate Bayesian computation analysis indicated that the two lineages experienced post-LGM range expansion and that postglacial secondary contact occurred in Sakhalin. Model projections into the year 2,100 predicted that S. japonica will shift northwards and lose its genetic diversity center on the Oshima Peninsula in Hokkaido and Shimokita Peninsula in Honshu. The range shifts and evolutionary history of S. japonica improve our understanding of how climate change impacted the distribution range and diversity of this species and provide useful information for the conservation of natural resources under ongoing environmental change in the NWP.

rPinecone: Define sub-lineages of a clonal expansion via a phylogenetic tree

The ability to distinguish different circulating pathogen clones from each other is a fundamental requirement to understand the epidemiology of infectious diseases. Phylogenetic analysis of genomic data can provide a powerful platform to identify lineages within bacterial populations, and thus inform outbreak investigation and transmission dynamics. However, resolving differences between pathogens associated with low-variant (LV) populations carrying low median pairwise single nucleotide variant (SNV) distances remains a major challenge. Here we present rPinecone, an R package designed to define sub-lineages within closely related LV populations. rPinecone uses a root-to-tip directional approach to define sub-lineages within a phylogenetic tree according to SNV distance from the ancestral node. The utility of this software was demonstrated using both simulated outbreaks and real genomic data of two LV populations: a hospital outbreak of methicillin-resistant Staphylococcus aureus and endemic Salmonella Typhi from rural Cambodia. rPinecone identified the transmission branches of the hospital outbreak and geographically confined lineages in Cambodia. Sub-lineages identified by rPinecone in both analyses were phylogenetically robust. It is anticipated that rPinecone can be used to discriminate between lineages of bacteria from LV populations where other methods fail, enabling a deeper understanding of infectious disease epidemiology for public health purposes.

Detection of vancomycin-resistant Enterococcus faecium hospital-adapted lineages in municipal wastewater treatment plants indicates widespread distribution and release into the environment

Vancomycin-resistant Enterococcus faecium (VREfm) is a leading cause of healthcare-associated infection. Reservoirs of VREfm are largely assumed to be nosocomial although there is a paucity of data on alternative sources. Here, we describe an integrated epidemiological and genomic analysis of E. faecium associated with bloodstream infection and isolated from wastewater. Treated and untreated wastewater from 20 municipal treatment plants in the East of England, United Kingdom was obtained and cultured to isolate E. faecium, ampicillin-resistant E. faecium (AREfm), and VREfm. VREfm was isolated from all 20 treatment plants and was released into the environment by 17/20 plants, the exceptions using terminal ultraviolet light disinfection. Median log₁₀ counts of AREfm and VREfm in untreated wastewater from 10 plants in direct receipt of hospital sewage were significantly higher than 10 plants that were not. We sequenced and compared the genomes of 423 isolates from wastewater with 187 isolates associated with bloodstream infection at five hospitals in the East of England. Among 481 E. faecium isolates belonging to the hospital-adapted clade, we observed genetic intermixing between wastewater and bloodstream infection, with highly related isolates shared between a major teaching hospital in the East of England and 9/20 plants. We detected 28 antibiotic resistance genes in the hospital-adapted clade, of which 23 were represented in bloodstream, hospital sewage, and municipal wastewater isolates. We conclude that our findings are consistent with widespread distribution of hospital-adapted VREfm beyond acute healthcare settings with extensive release of VREfm into the environment in the East of England.

Within-Species Genomic Variation and Variable Patterns of Recombination in the Tetracycline Producer Streptomyces rimosus

Streptomyces rimosus is best known as the primary source of the tetracycline class of antibiotics, most notably oxytetracycline, which have been widely used against many gram-positive and gram-negative pathogens and protozoan parasites. However, despite the medical and agricultural importance of S. rimosus, little is known of its evolutionary history and genome dynamics. In this study, we aim to elucidate the pan-genome characteristics and phylogenetic relationships of 32 S. rimosus genomes. The S. rimosus pan-genome contains more than 22,000 orthologous gene clusters, and approximately 8.8% of these genes constitutes the core genome. A large part of the accessory genome is composed of 9,646 strain-specific genes. S. rimosus exhibits an open pan-genome (decay parameter α = 0.83) and high gene diversity between strains (genomic fluidity φ = 0.12). We also observed strain-level variation in the distribution and abundance of biosynthetic gene clusters (BGCs) and that each individual S. rimosus genome has a unique repertoire of BGCs. Lastly, we observed variation in recombination, with some strains donating or receiving DNA more often than others, strains that tend to frequently recombine with specific partners, genes that often experience recombination more than others, and variable sizes of recombined DNA sequences. We conclude that the high levels of inter-strain genomic variation in S. rimosus is partly explained by differences in recombination among strains. These results have important implications on current efforts for natural drug discovery, the ecological role of strain-level variation in microbial populations, and addressing the fundamental question of why microbes have pan-genomes.

Genomic Analyses of >3,100 Nasopharyngeal Pneumococci Revealed Significant Differences Between Pneumococci Recovered in Four Different Geographical Regions

Understanding the structure of a bacterial population is essential in order to understand bacterial evolution. Estimating the core genome (those genes common to all, or nearly all, strains of a species) is a key component of such analyses. The size and composition of the core genome varies by dataset, but we hypothesized that the variation between different collections of the same bacterial species would be minimal. To investigate this, we analyzed the genome sequences of 3,118 pneumococci recovered from healthy individuals in Reykjavik (Iceland), Southampton (United Kingdom), Boston (United States), and Maela (Thailand). The analyses revealed a “supercore” genome (genes shared by all 3,118 pneumococci) of 558 genes, although an additional 354 core genes were shared by pneumococci from Reykjavik, Southampton, and Boston. Overall, the size and composition of the core and pan-genomes among pneumococci recovered in Reykjavik, Southampton, and Boston were similar. Maela pneumococci were distinctly different in that they had a smaller core genome and larger pan-genome. The pan-genome of Maela pneumococci contained several >25 Kb sequence regions (flanked by pneumococcal genes) that were homologous to genomic regions found in other bacterial species. Overall, our work revealed that some subsets of the global pneumococcal population are highly heterogeneous, and our hypothesis was rejected. This is an important finding in terms of understanding genetic variation among pneumococci and is also an essential point of consideration before generalizing the findings from a single dataset to the wider pneumococcal population.

Phylogeography of the Assassin Bug Sphedanolestes impressicollis in East Asia Inferred From Mitochondrial and Nuclear Gene Sequences

The assassin bug, Sphedanolestes impressicollis (Hemiptera: Reduviidae), is widely distributed in East Asia. It is an ideal model for evaluating the effects of climatic fluctuation and geographical events on the distribution patterns of East Asian reduviids. Here, we used two mitochondrial genes and one nuclear gene to investigate the phylogeographic pattern of the assassin bug based on comprehensive sampling in China, Japan, South Korea, Vietnam, and Laos. High levels of genetic differentiation were detected among the geographic populations classified into the northern and southern groups. A significant correlation was detected between genetic and geographical distances. The East China Sea land bridge served as a “dispersal corridor” during Pleistocene glaciation. The estimated divergence time indicated that the northern group may have separated from the eastern Chinese populations when the sea level rapidly rose during the “Ryukyu Coral Sea Stage” and the East China Sea land bridge was completely submerged. Demographic history and ecological niche modeling suggested that appropriate climatic conditions may have accounted for the rapid spread across the Korean Peninsula and Japan during the late Pleistocene. Our study underscores the pivotal roles of the Pleistocene sea level changes and climatic fluctuations in determining the distribution patterns of East Asian reduviids.

Genomic Diversity, Virulence, and Antimicrobial Resistance of Klebsiella pneumoniae Strains from Cows and Humans

Klebsiella pneumoniae is a leading cause of severe infections in humans and dairy cows, and these infections are rapidly becoming untreatable due to the emergence of multidrug-resistant (MDR) strains. However, little is known about the relationship between bovine and human K. pneumoniae isolates at the genome population level. Here, we investigated the genomic structures, pangenomic profiles, virulence determinants, and resistomes of 308 K. pneumoniae isolates from humans and dairy cows, including 96 newly sequenced cow isolates. We identified 177 functional protein families that were significantly different across human and bovine isolates; genes expressing proteins related to metal ion (iron, zinc, and calcium) metabolism were significantly more prevalent among the bovine isolates. Siderophore systems were found to be prevalent in both the bovine and the human isolates. In addition, we found that the Klebsiella ferric uptake operon kfuABC was significantly more prevalent in clinical mastitis cases than in healthy cows. Furthermore, on two dairy farms, we identified a unique IncN-type plasmid, pC5, coharboring bla _CTX-M-1 and mph(A) genes, which confer resistance to cephalosporins and macrolides, respectively. We provide here the complete annotated sequence of this plasmid.IMPORTANCE We demonstrate here the genetic diversity of K. pneumoniae isolates from dairy cows and the mixed phylogenetic lineages between bovine and human isolates. The ferric uptake operon kfuABC genes were more prevalent in strains from clinical mastitis cows. Furthermore, we report the emergence of an IncN-type plasmid carrying the bla _CTX-M-1 and mph(A) genes among dairy farms in the United States. Our study evaluated the genomic diversity of the bovine and human isolates, and the findings uncovered different profiles of virulence determinants among bovine and human K. pneumoniae isolates at the genome population level.

Genetic homogeneity of the Pacific thread herring (Opisthonema libertate) (Günther, 1867) in the Eastern Pacific, inferred from mtDNA sequences

In the present study, the population genetic structure of the Pacific thread herring (Opisthonema libertate) was analyzed through mitochondrial DNA (mtDNA) control region sequences. Organisms were collected from June 2015 to July 2015 from four commercial landing sites (Bahia Magdalena, Guaymas and Mazatlan, Mexico, and Puntarenas, Costa Rica) and one artisanal sampling site (Puerto La Libertad, El Salvador). A total of 125 sequences were analyzed. High levels of haplotype (h = 0.990) and nucleotide (π = 0.030) diversity were found. Pairwise Φst comparisons indicated differences attributed mainly to the organisms from El Salvador. However, Bayesian inferences did not support the existence of different populations. The haplotype distribution between locations did not show a clear phylogeographic pattern. Mismatched distribution showed a unimodal pattern for the five sampled areas, indicative of sudden demographic expansion. These results were supported by Bayesian skyline plot. Our results do not support the hypothesis that the Pacific thread herring presents a population genetic structure. Future genetic comparisons should include a larger number of samples as well as more polymorphic molecular markers to further support our results.

Integrated analysis of population genomics, transcriptomics and virulence provides novel insights into Streptococcus pyogenes pathogenesis

Streptococcus pyogenes causes 700 million human infections annually worldwide, yet, despite a century of intensive effort, there is no licensed vaccine against this bacterium. Although a number of large-scale genomic studies of bacterial pathogens have been published, the relationships among the genome, transcriptome, and virulence in large bacterial populations remain poorly understood. We sequenced the genomes of 2,101 emm28 S. pyogenes invasive strains, from which we selected 492 phylogenetically diverse strains for transcriptome analysis and 50 strains for virulence assessment. Data integration provided a novel understanding of the virulence mechanisms of this model organism. Genome-wide association study, expression quantitative trait loci analysis, machine learning, and isogenic mutant strains identified and confirmed a one-nucleotide indel in an intergenic region that significantly alters global transcript profiles and ultimately virulence. The integrative strategy that we used is generally applicable to any microbe and may lead to new therapeutics for many human pathogens.

Limited contribution of non-intensive chicken farming to ESBL-producing Escherichia coli colonization in humans in Vietnam: an epidemiological and genomic analysis

OBJECTIVES

To investigate the risk of colonization with ESBL-producing Escherichia coli (ESBL-Ec) in humans in Vietnam associated with non-intensive chicken farming.

METHODS

Faecal samples from 204 randomly selected farmers and their chickens, and from 306 age- and sex-matched community-based individuals who did not raise poultry were collected. Antimicrobial usage in chickens and humans was assessed by medicine cabinet surveys. WGS was employed to obtain a high-resolution genomic comparison between ESBL-Ec isolated from humans and chickens.

RESULTS

The adjusted prevalence of ESBL-Ec colonization was 20.0% (95% CI 10.8%-29.1%) and 35.2% (95% CI 30.4%-40.1%) in chicken farms and humans in Vietnam, respectively. Colonization with ESBL-Ec in humans was associated with antimicrobial usage (OR = 2.52, 95% CI = 1.08-5.87) but not with involvement in chicken farming. blaCTX-M-55 was the most common ESBL-encoding gene in strains isolated from chickens (74.4%) compared with blaCTX-M-27 in human strains (47.0%). In 3 of 204 (1.5%) of the farms, identical ESBL genes were detected in ESBL-Ec isolated from farmers and their chickens. Genomic similarity indicating recent sharing of ESBL-Ec between chickens and farmers was found in only one of these farms.

CONCLUSIONS

The integration of epidemiological and genomic data in this study has demonstrated a limited contribution of non-intensive chicken farming to ESBL-Ec colonization in humans in Vietnam and further emphasizes the importance of reducing antimicrobial usage in both human and animal host reservoirs.

Genomic epidemiology and antimicrobial resistance of Neisseria gonorrhoeae in New Zealand

Background

Antimicrobial-resistant Neisseria gonorrhoeae is a major threat to public health. No studies to date have examined the genomic epidemiology of gonorrhoea in the Western Pacific Region, where the incidence of gonorrhoea is particularly high.

Methods

A population-level study of N. gonorrhoeae in New Zealand (October 2014 to May 2015). Comprehensive susceptibility testing and WGS data were obtained for 398 isolates. Relatedness was inferred using phylogenetic trees, and pairwise core SNPs. Mutations and genes known to be associated with resistance were identified, and correlated with phenotype.

Results

Eleven clusters were identified. In six of these clusters, >25% of isolates were from females, while in eight of them, >15% of isolates were from females. Drug resistance was common; 98%, 32% and 68% of isolates were non-susceptible to penicillin, ciprofloxacin and tetracycline, respectively. Elevated MICs to extended-spectrum cephalosporins (ESCs) were observed in 3.5% of isolates (cefixime MICs ≥ 0.12 mg/L, ceftriaxone MICs ≥ 0.06 mg/L). Only nine isolates had penA XXXIV genotypes, three of which had decreased susceptibility to ESCs (MIC = 0.12 mg/L). Azithromycin non-susceptibility was identified in 43 isolates (10.8%); two of these isolates had 23S mutations (C2611T, 4/4 alleles), while all had mutations in mtrR or its promoter.

Conclusions

The high proportion of females in clusters suggests transmission is not exclusively among MSM in New Zealand; re-assessment of risk factors for transmission may be warranted in this context. As elevated MICs of ESCs and/or azithromycin were found in closely related strains, targeted public health interventions to halt transmission are urgently needed.

Population genomics of pneumococcal carriage in Massachusetts children following introduction of PCV-13

The 13-valent pneumococcal conjugate vaccine (PCV-13) was introduced in the United States in 2010. Using a large paediatric carriage sample collected from shortly after the introduction of PCV-7 to several years after the introduction of PCV-13, we investigate alterations in the composition of the pneumococcal population following the introduction of PCV-13, evaluating the extent to which the post-vaccination non-vaccine type (NVT) population mirrors that from prior to vaccine introduction and the effect of PCV-13 on vaccine type lineages. Draft genome assemblies from 736 newly sequenced and 616 previously published pneumococcal carriage isolates from children in Massachusetts between 2001 and 2014 were analysed. Isolates were classified into one of 22 sequence clusters (SCs) on the basis of their core genome sequence. We calculated the SC diversity for each sampling period as the probability that any two randomly drawn isolates from that period belong to different SCs. The sampling period immediately after the introduction of PCV-13 (2011) was found to have higher diversity than preceding (2007) or subsequent (2014) sampling periods {Simpson’s D 2007: 0.915 [95 % confidence interval (CI) 0.901, 0.929]; 2011:  0.935 [0.927, 0.942]; 2014 :  0.912 [0.901, 0.923]}. Amongst NVT isolates, we found the distribution of SCs in 2011 to be significantly different from that in 2007 or 2014 (Fisher’s exact test P=0.018, 0.0078), but did not find a difference comparing 2007 to 2014 (Fisher’s exact test P=0.24), indicating greater similarity between samples separated by a longer time period than between samples from closer time periods. We also found changes in the accessory gene content of the NVT population between 2007 and 2011 to have been reduced by 2014. Amongst the new serotypes targeted by PCV-13, four were present in our sample. The proportion of our sample composed of PCV-13-only vaccine serotypes 19A, 6C and 7F decreased between 2007 and 2014, but no such reduction was seen for serotype 3. We did, however, observe differences in the genetic composition of the pre- and post-PCV-13 serotype 3 population. Our isolates were collected during discrete sampling periods from a small geographical area, which may limit the generalizability of our findings. Pneumococcal diversity increased immediately following the introduction of PCV-13, but subsequently returned to pre-vaccination levels. This is reflected in the distribution of NVT lineages, and, to a lesser extent, their accessory gene frequencies. As such, there may be a period during which the population is particularly disrupted by vaccination before returning to a more stable distribution. The persistence and shifting genetic composition of serotype 3 is a concern and warrants further investigation.

Comparative Genomic Analysis Reveals the Potential Risk of Vibrio parahaemolyticus Isolated From Ready-To-Eat Foods in China

Vibrio parahaemolyticus is a major foodborne pathogen associated with the consumption of aquatic products. The presence of this bacterium in ready-to-eat (RTE) foods has recently been reported. However, the genomic features and potential risks of V. parahaemolyticus isolated from RTE foods are poorly understood. To help understand the genome-wide characteristics of RTE food isolates, the complete genomes of 27 RTE food isolates were sequenced and compared to those of 20 clinical and 19 other environmental (e.g., water and aquatic product source) isolates using a comparative genomics approach. Analysis revealed that V. parahaemolyticus RTE food isolates had higher numbers of genes on average and possessed more accessory genes than isolates from other sources. Most RTE food isolates were positive for some known virulence-associated genes and pathogenicity islands (PAIs), and some of these isolates were genetically homologous to clinical isolates. Genome-wide association analysis revealed 79 accessory genes and 78 missense single-nucleotide polymorphisms that affected 11 protein-coding genes were significantly associated with RTE food sources. These genes were mostly involved in defense mechanisms and energy production and conversion according to functional annotation in the COG database. KEGG Pathway analysis showed that these genes mainly affected the biofilm formation of V. parahaemolyticus, and subsequent experiments confirmed that nearly all RTE food isolates possessed the ability to form biofilm. The biofilm formation can facilitate the persistence of V. parahaemolyticus in RTE foods, and the presence of virulence-associated genes poses a pathogenic potential to humans. Our findings highlight the potential risk of V. parahaemolyticus in Chinese RTE foods and illustrate the genomic basis for the persistence of these isolates. This study will aid in re-evaluating the food safety threats conferred by this bacterium.

Evolutionary history and postglacial colonization of an Asian pit viper (Gloydius halys caucasicus) into Transcaucasia revealed by phylogenetic and phylogeographic analyses

It has been generally acknowledged that glacial climates at the time of the Pleistocene altered the patterns of species distributions, prompting latitudinal and altitudinal distribution shifts in several species, including poikilothermic species commonly known for their thermal sensitivity. However, the historical phylogeographic patterns of such species have remained largely unknown. Here, we present the historical biogeographic, phylogenetic, and phylogeographic relationships of the Caucasian pit viper, G. h. caucasicus, based on two mtDNA (cyt b and ND4) and one nDNA (c-mos) genes. This pit viper represents the westernmost member of the Crotalinae subfamily in the Palearctic and occurs in a variety of habitats, from 30 m to 3,000 m above sea level. In Iran, it is distributed on the northern and southern slopes of the Alborz Mountains, rendering it a target for phylogenetic and phylogeographic studies of a terrestrial poikilothermic animal. Our study identified four Iranian lineages of G. h. caucasicus along the northeastern to northwestern slopes of the Alborz Mountains and southern Azerbaijan (Talysh Mountains). Diversification of the Iranian lineages highlights population expansion and subsequent isolation into four plausible refugial areas during the Quaternary paleo-climatic oscillations, confirmed by our molecular dating and historical biogeographic analyses. The results of coalescence-based simulations support the incursion of the species from northeastern Iran to the western end of the Alborz, and then toward Transcaucasia via two directions: northern and southern slopes of the Alborz Mountains. Furthermore, our results clearly implied that G. h. caucasicus should be elevated to species rank and further referred to as G. caucasicus (Nikolsky, 1916).

Fast and flexible bacterial genomic epidemiology with PopPUNK

The routine use of genomics for disease surveillance provides the opportunity for high-resolution bacterial epidemiology. Current whole-genome clustering and multilocus typing approaches do not fully exploit core and accessory genomic variation, and they cannot both automatically identify, and subsequently expand, clusters of significantly similar isolates in large data sets spanning entire species. Here, we describe PopPUNK (Population Partitioning Using Nucleotide K -mers), a software implementing scalable and expandable annotation- and alignment-free methods for population analysis and clustering. Variable-length k-mer comparisons are used to distinguish isolates' divergence in shared sequence and gene content, which we demonstrate to be accurate over multiple orders of magnitude using data from both simulations and genomic collections representing 10 taxonomically widespread species. Connections between closely related isolates of the same strain are robustly identified, despite interspecies variation in the pairwise distance distributions that reflects species' diverse evolutionary patterns. PopPUNK can process 103-104 genomes in a single batch, with minimal memory use and runtimes up to 200-fold faster than existing model-based methods. Clusters of strains remain consistent as new batches of genomes are added, which is achieved without needing to reanalyze all genomes de novo. This facilitates real-time surveillance with consistent cluster naming between studies and allows for outbreak detection using hundreds of genomes in minutes. Interactive visualization and online publication is streamlined through the automatic output of results to multiple platforms. PopPUNK has been designed as a flexible platform that addresses important issues with currently used whole-genome clustering and typing methods, and has potential uses across bacterial genetics and public health research.

Independently evolved and gene flow-accelerated pesticide resistance in two-spotted spider mites

Pest species are often able to develop resistance to pesticides used to control them, depending on how rapidly resistance can emerge within a population or spread from another resistant population. We examined the evolution of bifenazate resistance in China in the two-spotted spider mite (TSSM) Tetranychus uticae Koch (Acari: Tetranychidae), one of the most resistant arthropods, by using bioassays, detection of mutations in the target cytb gene, and population genetic structure analysis using microsatellite markers. Bioassays showed variable levels of resistance to bifenazate. The cytb mutation G126S, which confers medium resistance in TSSM to bifenazate, had previously been detected prior to the application of bifenazate and was now widespread, suggesting likely resistance evolution from standing genetic variation. G126S was detected in geographically distant populations across different genetic clusters, pointing to the independent origin of this mutation in different TSSM populations. A novel A269V mutation linked to a low-level resistance was detected in two southern populations. Widespread resistance associated with a high frequency of the G126S allele was found in four populations from the Beijing area which were not genetically differentiated. In this case, a high level of gene flows likely accelerated the development of resistance within this local region, as well as into an outlying region distant from Beijing. These findings, therefore, suggest patterns consistent with both local evolution of pesticide resistance as well as an impact of migration, helping to inform resistance management strategies in TSSM.

Phylogeography and population genetics of the endemic Malagasy bat, Macronycteris commersoni s.s. (Chiroptera: Hipposideridae)

Macronycteris commersoni (Hipposideridae), a bat species endemic to Madagascar, is widespread across the island and utilizes a range of habitat types including open woodland, degraded habitats, and forested areas from sea level to 1,325 m. Despite being widely distributed, there is evidence that M. commersoni exhibits morphological and bioacoustic variation across its geographical range. We investigated the fine-scale phylogeographic structure of populations in the western half of the island using extensive spatial sampling and sequence data from two mitochondrial DNA regions. Our results indicated several lineages within M. commersoni. Individuals collected from northern Madagascar formed a single monophyletic clade (clade C). A second clade (clade B) included individuals collected from the south-western portion of the island. This second clade displayed more phylogeographical partitioning with differences in mtDNA haplotypes frequency detected between populations collected in different bioclimatic regions. Lineage dispersal, genetic divergence, and timing of expansion events of M. commersoni were probably associated with Pleistocene climate fluctuations. Our data suggest that the northern and the central western regions of Madagascar may have acted as refugia for this species during periods of cooler and drier climate conditions associated with the Pleistocene.

Population Structure, Antibiotic Resistance, and Uropathogenicity of Klebsiella variicola

Infections caused by antibiotic-resistant bacterial pathogens are a growing public health threat. Understanding of pathogen relatedness and biology is imperative for tracking outbreaks and developing therapeutics. Here, we detail the phylogenetic structure of 145 K. variicola genomes from different continents. Our results have important clinical ramifications as high-risk antibiotic resistance genes are present in K. variicola genomes from a variety of geographic locations and as we demonstrate that K. variicola clinical isolates can establish higher bladder titers than K. pneumoniae. Differential presence of these pilus genes inK. variicola isolates may indicate adaption for specific environmental niches. Therefore, due to the potential of multidrug resistance and pathogenic efficacy, identification of K. variicola and K. pneumoniae to a species level should be performed to optimally improve patient outcomes during infection. This work provides a foundation for our improved understanding of K. variicola biology and pathogenesis.

Klebsiella variicola is a member of the Klebsiella genus and often misidentified as Klebsiella pneumoniae or Klebsiella quasipneumoniae. The importance of K. pneumoniae human infections has been known; however, a dearth of relative knowledge exists for K. variicola. Despite its growing clinical importance, comprehensive analyses of K. variicola population structure and mechanistic investigations of virulence factors and antibiotic resistance genes have not yet been performed. To address this, we utilized in silico, in vitro, and in vivo methods to study a cohort of K. variicola isolates and genomes. We found that the K. variicola population structure has two distant lineages composed of two and 143 genomes, respectively. Ten of 145 K. variicola genomes harbored carbapenem resistance genes, and 6/145 contained complete virulence operons. While the β-lactam bla_LEN and quinolone oqxAB antibiotic resistance genes were generally conserved within our institutional cohort, unexpectedly 11 isolates were nonresistant to the β-lactam ampicillin and only one isolate was nonsusceptible to the quinolone ciprofloxacin. K. variicola isolates have variation in ability to cause urinary tract infections in a newly developed murine model, but importantly a strain had statistically significant higher bladder CFU than the model uropathogenic K. pneumoniae strain TOP52. Type 1 pilus and genomic identification of altered fim operon structure were associated with differences in bladder CFU for the tested strains. Nine newly reported types of pilus genes were discovered in the K. variicola pan-genome, including the first identified P-pilus in Klebsiella spp.

Genome-wide discovery of epistatic loci affecting antibiotic resistance in Neisseria gonorrhoeae using evolutionary couplings

Genome analysis should allow the discovery of interdependent loci that together cause antibiotic resistance. In practice, however, the vast number of possible epistatic interactions erodes statistical power. Here, we extend an approach that has been successfully used to identify epistatic residues in proteins to infer genomic loci that are strongly coupled. This approach reduces the number of tests required for an epistatic genome-wide association study of antibiotic resistance and increases the likelihood of identifying causal epistasis. We discovered 38 loci and 240 epistatic pairs that influence the minimum inhibitory concentrations of 5 different antibiotics in 1,102 isolates of Neisseria gonorrhoeae that were confirmed in a second dataset of 495 isolates. Many known resistance-affecting loci were recovered; however, the majority of associations occurred in unreported genes, such as murE. About half of the discovered epistasis involved at least one locus previously associated with antibiotic resistance, including interactions between gyrA and parC. Still, many combinations involved unreported loci and genes. While most variation in minimum inhibitory concentrations could be explained by identified loci, epistasis substantially increased explained phenotypic variance. Our work provides a systematic identification of epistasis affecting antibiotic resistance in N. gonorrhoeae and a generalizable approach for epistatic genome-wide association studies.

Landscape permeability and individual variation in a dispersal-linked gene jointly determine genetic structure in the Glanville fritillary butterfly

There is now clear evidence that species across a broad range of taxa harbor extensive heritable variation in dispersal. While studies suggest that this variation can facilitate demographic outcomes such as range expansion and invasions, few have considered the consequences of intraspecific variation in dispersal for the maintenance and distribution of genetic variation across fragmented landscapes. Here, we examine how landscape characteristics and individual variation in dispersal combine to predict genetic structure using genomic and spatial data from the Glanville fritillary butterfly. We used linear and latent factor mixed models to identify the landscape features that best predict spatial sorting of alleles in the dispersal-related gene phosphoglucose isomerase (Pgi). We next used structural equation modeling to test if variation in Pgi mediated gene flow as measured by Fst at putatively neutral loci. In a year when the population was recovering following a large decline, individuals with a genotype associated with greater dispersal ability were found at significantly higher frequencies in populations isolated by water and forest, and these populations showed lower levels of genetic differentiation at neutral loci. These relationships disappeared in the next year when metapopulation density was high, suggesting that the effects of individual variation are context dependent. Together our results highlight that (1) more complex aspects of landscape structure beyond just the configuration of habitat can be important for maintaining spatial variation in dispersal traits and (2) that individual variation in dispersal plays a key role in maintaining genetic variation across fragmented landscapes.

Genome-wide analysis of Streptococcus pneumoniae serogroup 19 in the decade after the introduction of pneumococcal conjugate vaccines in Australia

The decline in invasive pneumococcal disease (IPD), following the introduction of the 7-valent pneumococcal conjugate vaccination (PCV-7), was tempered by emergence of non-vaccine serotypes, particularly 19A. In Australia, three years after PCV-7 was replaced by PCV-13, containing 19A and 19F antigens, serogroup 19 was still a prominent cause of IPD in children under five. In this study we examined the evolution of serogroup 19 before and after introduction of paediatric vaccines in New South Wales (NSW), Australia. Genomes of 124 serogroup 19 IPD isolates collected before (2004) and after introduction of PCV-7 (2008) and PCV-13 (2014), from children under five in NSW, were analysed. Eleven core genome sequence clusters (cgSC) and 35 multilocus sequence types (ST) were identified. The majority (78/124) of the isolates belonged to four cgSCs: cgSC7 (ST199), cgSC11 (ST320), cgSC8 (ST63) and cgSC9 (ST2345). ST63 and ST2345 were exclusively serotype 19A and accounted for its predominantly intermediate penicillin resistance; these two clusters first appeared in 2008 and largely disappeared after introduction of PCV-13. Serogroup 19 was responsible for the highest proportion of vaccine failures in NSW. Relatively low immunogenicity of serogroup 19 antigens and Australia's three-dose vaccine schedule could affect the population dynamics of this serogroup.

Genomic Surveillance of Enterococcus faecium Reveals Limited Sharing of Strains and Resistance Genes between Livestock and Humans in the United Kingdom

The rise in rates of human infection caused by vancomycin-resistant Enterococcus faecium (VREfm) strains between 1988 to the 2000s in Europe was suggested to be associated with acquisition from livestock. As a result, the European Union banned the use of the glycopeptide drug avoparcin as a growth promoter in livestock feed. While some studies reported a decrease in VREfm in livestock, others reported no reduction. Here, we report the first livestock VREfm prevalence survey in the UK since 2003 and the first large-scale study using whole-genome sequencing to investigate the relationship between E. faecium strains in livestock and humans. We found a low prevalence of VREfm in retail meat and limited evidence for recent sharing of strains between livestock and humans with bloodstream infection. There was evidence for limited sharing of genes encoding antibiotic resistance between these reservoirs, a finding which requires further research.

Vancomycin-resistant Enterococcus faecium (VREfm) is a major cause of nosocomial infection and is categorized as high priority by the World Health Organization global priority list of antibiotic-resistant bacteria. In the past, livestock have been proposed as a putative reservoir for drug-resistant E. faecium strains that infect humans, and isolates of the same lineage have been found in both reservoirs. We undertook cross-sectional surveys to isolate E. faecium (including VREfm) from livestock farms, retail meat, and wastewater treatment plants in the United Kingdom. More than 600 isolates from these sources were sequenced, and their relatedness and antibiotic resistance genes were compared with genomes of almost 800 E. faecium isolates from patients with bloodstream infection in the United Kingdom and Ireland. E. faecium was isolated from 28/29 farms; none of these isolates were VREfm, suggesting a decrease in VREfm prevalence since the last UK livestock survey in 2003. However, VREfm was isolated from 1% to 2% of retail meat products and was ubiquitous in wastewater treatment plants. Phylogenetic comparison demonstrated that the majority of human and livestock-related isolates were genetically distinct, although pig isolates from three farms were more genetically related to human isolates from 2001 to 2004 (minimum of 50 single-nucleotide polymorphisms [SNPs]). Analysis of accessory (variable) genes added further evidence for distinct niche adaptation. An analysis of acquired antibiotic resistance genes and their variants revealed limited sharing between humans and livestock. Our findings indicate that the majority of E. faecium strains infecting patients are largely distinct from those from livestock in this setting, with limited sharing of strains and resistance genes.

Global emergence and population dynamics of divergent serotype 3 CC180 pneumococci

Streptococcus pneumoniae serotype 3 remains a significant cause of morbidity and mortality worldwide, despite inclusion in the 13-valent pneumococcal conjugate vaccine (PCV13). Serotype 3 increased in carriage since the implementation of PCV13 in the USA, while invasive disease rates remain unchanged. We investigated the persistence of serotype 3 in carriage and disease, through genomic analyses of a global sample of 301 serotype 3 isolates of the Netherlands3-31 (PMEN31) clone CC180, combined with associated patient data and PCV utilization among countries of isolate collection. We assessed phenotypic variation between dominant clades in capsule charge (zeta potential), capsular polysaccharide shedding, and susceptibility to opsonophagocytic killing, which have previously been associated with carriage duration, invasiveness, and vaccine escape. We identified a recent shift in the CC180 population attributed to a lineage termed Clade II, which was estimated by Bayesian coalescent analysis to have first appeared in 1968 [95% HPD: 1939-1989] and increased in prevalence and effective population size thereafter. Clade II isolates are divergent from the pre-PCV13 serotype 3 population in non-capsular antigenic composition, competence, and antibiotic susceptibility, the last of which resulting from the acquisition of a Tn916-like conjugative transposon. Differences in recombination rates among clades correlated with variations in the ATP-binding subunit of Clp protease, as well as amino acid substitutions in the comCDE operon. Opsonophagocytic killing assays elucidated the low observed efficacy of PCV13 against serotype 3. Variation in PCV13 use among sampled countries was not independently correlated with the CC180 population shift; therefore, genotypic and phenotypic differences in protein antigens and, in particular, antibiotic resistance may have contributed to the increase of Clade II. Our analysis emphasizes the need for routine, representative sampling of isolates from disperse geographic regions, including historically under-sampled areas. We also highlight the value of genomics in resolving antigenic and epidemiological variations within a serotype, which may have implications for future vaccine development.

Genomic Analysis of Multiresistant Staphylococcus capitis Associated with Neonatal Sepsis

Coagulase-negative staphylococci (CoNS), such as Staphylococcus capitis, are major causes of bloodstream infections in neonatal intensive care units (NICUs). Recently, a distinct clone of S. capitis (designated S. capitis NRCS-A) has emerged as an important pathogen in NICUs internationally.

Coagulase-negative staphylococci (CoNS), such as Staphylococcus capitis, are major causes of bloodstream infections in neonatal intensive care units (NICUs). Recently, a distinct clone of S. capitis (designated S. capitis NRCS-A) has emerged as an important pathogen in NICUs internationally. Here, 122 S. capitis isolates from New Zealand (NZ) underwent whole-genome sequencing (WGS), and these data were supplemented with publicly available S. capitis sequence reads. Phylogenetic and comparative genomic analyses were performed, as were phenotypic assessments of antimicrobial resistance, biofilm formation, and plasmid segregational stability on representative isolates. A distinct lineage of S. capitis was identified in NZ associated with neonates and the NICU environment. Isolates from this lineage produced increased levels of biofilm, displayed higher levels of tolerance to chlorhexidine, and were multidrug resistant. Although similar to globally circulating NICU-associated S. capitis strains at a core-genome level, NZ NICU S. capitis isolates carried a novel stably maintained multidrug-resistant plasmid that was not present in non-NICU isolates. Neonatal blood culture isolates were indistinguishable from environmental S. capitis isolates found on fomites, such as stethoscopes and neonatal incubators, but were generally distinct from those isolates carried by NICU staff. This work implicates the NICU environment as a potential reservoir for neonatal sepsis caused by S. capitis and highlights the capacity of genomics-based tracking and surveillance to inform future hospital infection control practices aimed at containing the spread of this important neonatal pathogen.

Genome-wide identification of geographical segregated genetic markers in Salmonella enterica serovar Typhimurium variant 4,[5],12:i:

Salmonella enterica ser. Typhimurium monophasic variant 4,[5],12:i:- has been associated with food-borne epidemics worldwide and swine appeared to be the main reservoir in most of the countries of isolation. However, the monomorphic nature of this serovar has, so far, hindered identification of the source due to expansion of clonal lineages in multiple hosts and food producing systems. Since geographically structured genetic signals can shape bacterial populations, identification of biogeographical markers in S. 1,4,[5],12:i:- genomes can contribute to improving source attribution. In this study, the phylogeographical structure of 148 geographically and temporally related Italian S. 1,4,[5],12:i:- has been investigated. The Italian isolates belong to a large population of clonal S. Typhimurium/1,4,[5],12:i:- isolates collected worldwide in two decades showing up to 2.5% of allele differences. Phylogenetic reconstruction revealed that isolates from the same geographical origin form highly supported monophyletic groups, suggesting discrete geographical segregation. These monophyletic groups are characterized by the gene content of a large sopE-containing prophage. Within this prophage, genome-wide comparison identified several genes overrepresented in strains of Italian origin. This suggests that certain lineages may be characterized by the acquisition of specific accessory genetic markers useful for improving identification of the source in ongoing epidemics.

Multiple refugia from penultimate glaciations in East Asia demonstrated by phylogeography and ecological modelling of an insect pest

BACKGROUND

Refugial populations in Quaternary glaciations are critical to understanding the evolutionary history and climatic interactions of many extant species. Compared with the well-studied areas of Europe and Northern America, refugia of species in eastern Asia remain largely unknown. Here, we investigated the phylogeographic history of a globally important insect pest, the oriental fruit moth Grapholita molesta, in its native range of China.

RESULTS

Genetic structure analyses unveiled three distinct groups and a set of populations with admixture. Approximate Bayesian Computation (ABC) analyses support range expansion of this moth from southwest groups of Yunnan and Sichuan to northern and eastern China. A set of admixed populations was found around these two ancestral groups. This pattern of genetic structure points to two refugia located in the Yunnan region and Sichuan Basin. The split of the two refugia was dated to 329.2 thousand years ago in the penultimate glacial period. One of the lineages was exclusively found around the Sichuan Basin, indicating the formation of endemic populations in this refugium. Ecological niche model analysis suggested a shrinking distribution from the LIG period to the MID period in the Sichuan lineage but a wide and stable distribution in the other lineage.

CONCLUSIONS

Our results for the first time suggest that Yunnan and Sichuan jointly served as two large-scale refugia in eastern Asia in Quaternary glaciations, helping to maintain genetic diversity overall.

Population genetic analysis of Chadian Guinea worms reveals that human and non-human hosts share common parasite populations

Following almost 10 years of no reported cases, Guinea worm disease (GWD or dracunculiasis) reemerged in Chad in 2010 with peculiar epidemiological patterns and unprecedented prevalence of infection among non-human hosts, particularly domestic dogs. Since 2014, animal infections with Guinea worms have also been observed in the other three countries with endemic transmission (Ethiopia, Mali, and South Sudan), causing concern and generating interest in the parasites' true taxonomic identity and population genetics. We present the first extensive population genetic data for Guinea worm, investigating mitochondrial and microsatellite variation in adult female worms from both human and non-human hosts in the four endemic countries to elucidate the origins of Chad's current outbreak and possible host-specific differences between parasites. Genetic diversity of Chadian Guinea worms was considerably higher than that of the other three countries, even after controlling for sample size through rarefaction, and demographic analyses are consistent with a large, stable parasite population. Genealogical analyses eliminate the other three countries as possible sources of parasite reintroduction into Chad, and sequence divergence and distribution of genetic variation provide no evidence that parasites in human and non-human hosts are separate species or maintain isolated transmission cycles. Both among and within countries, geographic origin appears to have more influence on parasite population structure than host species. Guinea worm infection in non-human hosts has been occasionally reported throughout the history of the disease, particularly when elimination programs appear to be reaching their end goals. However, no previous reports have evaluated molecular support of the parasite species identity. Our data confirm that Guinea worms collected from non-human hosts in the remaining endemic countries of Africa are Dracunculus medinensis and that the same population of worms infects both humans and dogs in Chad. Our genetic data and the epidemiological evidence suggest that transmission in the Chadian context is currently being maintained by canine hosts.

Global spread of three multidrug-resistant lineages of Staphylococcus epidermidis

Staphylococcus epidermidis is a conspicuous member of the human microbiome, widely present on healthy skin. Here we show that S. epidermidis has also evolved to become a formidable nosocomial pathogen. Using genomics, we reveal that three multidrug-resistant, hospital-adapted lineages of S. epidermidis (two ST2 and one ST23) have emerged in recent decades and spread globally. These lineages are resistant to rifampicin through acquisition of specific rpoB mutations that have become fixed in the populations. Analysis of isolates from 96 institutions in 24 countries identified dual D471E and I527M RpoB substitutions to be the most common cause of rifampicin resistance in S. epidermidis, accounting for 86.6% of mutations. Furthermore, we reveal that the D471E and I527M combination occurs almost exclusively in isolates from the ST2 and ST23 lineages. By breaching lineage-specific DNA methylation restriction modification barriers and then performing site-specific mutagenesis, we show that these rpoB mutations not only confer rifampicin resistance, but also reduce susceptibility to the last-line glycopeptide antibiotics, vancomycin and teicoplanin. Our study has uncovered the previously unrecognized international spread of a near pan-drug-resistant opportunistic pathogen, identifiable by a rifampicin-resistant phenotype. It is possible that hospital practices, such as antibiotic monotherapy utilizing rifampicin-impregnated medical devices, have driven the evolution of this organism, once trivialized as a contaminant, towards potentially incurable infections.

Macroevolution of gastric Helicobacter species unveils interspecies admixture and time of divergence

Since the discovery of the human pathogen Helicobacter pylori, various other Helicobacter species have been identified in the stomach of domesticated and wild mammals. To better understand the evolutionary history of these ecologically similar but genetically distinct species, we analyzed 108 gastric Helicobacter genomes and included 54 enterohepatic Helicobacter genomes for comparison purposes. An admixture analysis supported the presence of an ecological barrier, preventing the genetic exchange between the gastric and enterohepatic Helicobacter species, and unraveled many gene flow events within and across species residing in the stomach. As pets can be colonized by multiple gastric Helicobacter species, the genetic exchange between the canine and feline strains was evident, with H. heilmannii and H. bizzozeronii showing the highest interspecies recombination. An admixture between H. pylori (in particular, the ancestral African strains), H. acinonychis from wild felines and H. cetorum from marine mammals was also identified. Because these latter species do not share the same host, this phenomenon is most likely a remaining signal of shared ancestry. A reconstruction of the time of divergence of the gastric Helicobacter spp. revealed that the domestic animal-related Helicobacter species evolved in parallel with H. pylori and its two closest relatives (H. acinonychis and H. cetorum), rather than together.

Some strains that have converged to infect Prunus spp. trees are members of distinct Pseudomonas syringae genomospecies and ecotypes as revealed by in silico genomic comparison

A complementary taxonomic and population genetic study was performed to delineate genetically and ecologically distinct species within the Pseudomonas syringae complex by assessing 16 strains including pathovar strains that have converged to infect Prunus spp. trees, and two outgroups. Both average nucleotide identity and genome-to-genome distance comparison methods revealed the occurrence of distinct genomospecies, namely 1, 2, 3 and 8 (sensu Gardan et al.), with the latter two being closely related. Strains classified as P. s. pv. morsprunorum clustered into two distinct genomospecies, namely 2 and 8. Both the AdaptML and hierarchical Bayesian analysis of population structure methods highlighted the presence of three ecotypes, and the taxonomically related genomospecies 3 and 8 strains were members of the same ecotype. The distribution of pathogenic and virulence-associated genetic traits among Pseudomonas strains did not reveal any distinct type III secretion system effector or phytotoxin distribution pattern that characterized single genomospecies and strains that infect Prunus spp. The complete WHOP (Woody HOst and Pseudomonas spp.) genomic region and the entire β-ketoadipate gene cluster, including the catBCA operon, were found only in the members of genomospecies 2 and in the two P. s. pv. morsprunorum strains of genomospecies 8. A reduced gene flow between the three ecotypes suggested that point mutations played a larger role during the evolution of the strains than recombination. Our data support the idea that Prunus trees can be infected by different strains of distinct Pseudomonas genomospecies/ecotypes through diverse mechanisms of host colonization and infection. Such strains may represent particular lineages that emerged from environments other than that of the infected plant upon acquiring genetic traits that gave them the ability to cause plant diseases. The complementary assessment of bacterial strains using both taxonomic approaches and methods that reveal ecologically homogeneous populations has proven useful in confirming the cohesion of bacterial clusters.

New Variant of Multidrug-Resistant Salmonella enterica Serovar Typhimurium Associated with Invasive Disease in Immunocompromised Patients in Vietnam

Salmonella Typhimurium is a major diarrheal pathogen and associated with invasive nontyphoid Salmonella (iNTS) disease in vulnerable populations. We present the first characterization of iNTS organisms in Southeast Asia and describe a different evolutionary trajectory from that of organisms causing iNTS in sub-Saharan Africa. In Vietnam, the globally distributed monophasic variant of Salmonella Typhimurium, the serovar I:4,[5],12:i:− ST34 clone, has reacquired a phase 2 flagellum and gained a multidrug-resistant plasmid to become associated with iNTS disease in HIV-infected patients. We document distinct communities of S. Typhimurium and I:4,[5],12:i:− in animals and humans in Vietnam, despite the greater mixing of these host populations here. These data highlight the importance of whole-genome sequencing surveillance in a One Health context in understanding the evolution and spread of resistant bacterial infections.

Nontyphoidal Salmonella (NTS), particularly Salmonella enterica serovar Typhimurium, is among the leading etiologic agents of bacterial enterocolitis globally and a well-characterized cause of invasive disease (iNTS) in sub-Saharan Africa. In contrast, S. Typhimurium is poorly defined in Southeast Asia, a known hot spot for zoonotic disease with a recently described burden of iNTS disease. Here, we aimed to add insight into the epidemiology and potential impact of zoonotic transfer and antimicrobial resistance (AMR) in S. Typhimurium associated with iNTS and enterocolitis in Vietnam. We performed whole-genome sequencing and phylogenetic reconstruction on 85 human (enterocolitis, carriage, and iNTS) and 113 animal S. Typhimurium isolates isolated in Vietnam. We found limited evidence for the zoonotic transmission of S. Typhimurium. However, we describe a chain of events where a pandemic monophasic variant of S. Typhimurium (serovar I:4,[5],12:i:− sequence type 34 [ST34]) has been introduced into Vietnam, reacquired a phase 2 flagellum, and acquired an IncHI2 multidrug-resistant plasmid. Notably, these novel biphasic ST34 S. Typhimurium variants were significantly associated with iNTS in Vietnamese HIV-infected patients. Our study represents the first characterization of novel iNTS organisms isolated outside sub-Saharan Africa and outlines a new pathway for the emergence of alternative Salmonella variants into susceptible human populations.

Pneumococcal vaccine impacts on the population genomics of non-typeable Haemophilus influenzae

The implementation of pneumococcal conjugate vaccines (PCVs) has led to a decline in vaccine-type disease. However, there is evidence that the epidemiology of non-typeable Haemophilus influenzae (NTHi) carriage and disease can be altered as a consequence of PCV introduction. We explored the epidemiological shifts in NTHi carriage using whole genome sequencing over a 5-year period that included PCV13 replacement of PCV7 in the UK’s National Immunization Programme in 2010. Between 2008/09 and 2012/13 (October to March), nasopharyngeal swabs were taken from children <5 years of age. Significantly increased carriage post-PCV13 was observed and lineage-specific associations with Streptococcus pneumoniae were seen before but not after PCV13 introduction. NTHi were characterized into 11 discrete, temporally stable lineages, congruent with current knowledge regarding the clonality of NTHi. The increased carriage could not be linked to the expansion of a particular clone and different co-carriage dynamics were seen before PCV13 implementation when NTHi co-carried with vaccine serotype pneumococci. In summary, PCV13 introduction has been shown to have an indirect effect on NTHi epidemiology and there exists both negative and positive, distinct associations between pneumococci and NTHi. This should be considered when evaluating the impacts of pneumococcal vaccine design and policy.

RhierBAPS: An R implementation of the population clustering algorithm hierBAPS

Identifying structure in collections of sequence data sets remains a common problem in genomics. hierBAPS, a popular algorithm for identifying population structure in haploid genomes, has previously only been available as a MATLAB binary. We provide an R implementation which is both easier to install and use, automating the entire pipeline. Additionally, we allow for the use of multiple processors, improve on the default settings of the algorithm, and provide an interface with the ggtree library to enable informative illustration of the clustering results. Our aim is that this package aids in the understanding and dissemination of the method, as well as enhancing the reproducibility of population structure analyses.

Late Pleistocene genetic diversification and demographic expansion in the widely distributed neotropical ant Neoponera villosa (Ponerinae)

Phylogeographic studies of continent-wide distributed species are key to understand population dynamics processes that occurred at large geographical scales. Here, we examined two mitochondrial (mt) DNA sequence (COI, Cyt b) and eight nuclear microsatellites markers to investigate the cohesiveness, genetic diversity and demographic history of Neoponera villosa (Fabricius), a ponerine ant species widely distributed along most part of the Neotropics and southern Nearctic. The reconstructed phylogeny and mt variation supported the cohesiveness of the examined populations of N. villosa. The species probably originated in South America during the late Pliocene/middle Pleistocene and subsequently dispersed to Central America and the Transitional Nearctic-Neotropical zone during the late Pleistocene, with an increase in its population size ca. 30 thousand years ago. The limited phylogeographic structure observed in N. villosa supports its late Pleistocene range expansion and gene flow among distant geographic areas in central and southern Mexico and Central America.

Anthropogenic and environmental factors associated with high incidence of mcr-1 carriage in humans across China

MCR-1-positve Escherichia coli (MCRPEC) have been reported in humans worldwide; however, thus far, their prevalence is low and potential sources for human mcr-1 carriage have not yet been identified. Here, we analyse a nationwide epidemiological dataset on MCRPEC in humans throughout China and assess the factors associated with MCRPEC carriage using natural and national anthropogenic data. We identified 774 non-duplicate MCRPEC isolates from 774 stool samples collected from 5,159 healthy individuals in 30 provinces and municipalities in 2016, with a prevalence of MCRPEC ranging from 3.7 to 32.7% (average: 15.0%)-substantially higher than previously reported. MCRPEC carriage was associated with provincial regions, the production of sheep and freshwater aquaculture, annual consumption of total meat, pork and mutton, and daily intake of aquaculture products. MCRPEC was significantly more prevalent in provinces with higher aquaculture industries. Whole-genome sequencing analysis revealed that the MCRPEC isolates were clustered into four distinct lineages, two of which were dominant and harboured most of the MCRPEC isolates. The high prevalence of MCRPEC in the community poses a substantial risk for colistin usage in clinical practice and suggests the need for intestinal screening of mcr-1 carriers in intensive care units in Chinese hospitals. Furthermore, our data suggest that aquaculture is a significant reservoir of mcr-1.

Invasion of the assassin bug Agriosphodrus dohrni (Hemiptera: Reduviidae) to Japan: Source estimation inferred from mitochondrial and nuclear gene sequences

A large-sized assassin bug Agriosphodrus dohrni (Signoret), has been recorded from India, Vietnam, China and Japan. It is one of the potential biological control agents against some important agricultural and forest pests. This species is speculated to have invaded Japan from its native range in China about 60 years ago. We used three mitochondrial gene fragments (COI, Cytb, and ND5) and one nuclear gene fragment (EF-1α) to clarify the invasion history of A. dohrni and assess the effects of geographic events and associated ecological adaptation on the distribution pattern. The native populations of A. dohrni in China are divided into three distinct groups, which might be molded by the Early Pleistocene glaciation event and diverged during the Calabrian Stage. However, consistent with the hypothesis of a recent invasion, extremely low level of genetic variation was detected in the Japanese populations, with only two haplotypes for the combined mitochondrial genes. Both the splits network and the ML/BI phylogenetic trees revealed that haplotypes of Japan were more closely-related to those from eastern China. Therefore, we postulate that there has been only one introduction event, probably from somewhere around the Nanjing (NJ) and Lin'an (LA) populations of eastern China.

Evolution of HIV-1 within untreated individuals and at the population scale in Uganda

HIV-1 undergoes multiple rounds of error-prone replication between transmission events, resulting in diverse viral populations within and among individuals. In addition, the virus experiences different selective pressures at multiple levels: during the course of infection, at transmission, and among individuals. Disentangling how these evolutionary forces shape the evolution of the virus at the population scale is important for understanding pathogenesis, how drug- and immune-escape variants are likely to spread in populations, and the development of preventive vaccines. To address this, we deep-sequenced two regions of the HIV-1 genome (p24 and gp41) from 34 longitudinally-sampled untreated individuals from Rakai District in Uganda, infected with subtypes A, D, and inter-subtype recombinants. This dataset substantially increases the availability of HIV-1 sequence data that spans multiple years of untreated infection, in particular for different geographical regions and viral subtypes. In line with previous studies, we estimated an approximately five-fold faster rate of evolution at the within-host compared to the population scale for both synonymous and nonsynonymous substitutions, and for all subtypes. We determined the extent to which this mismatch in evolutionary rates can be explained by the evolution of the virus towards population-level consensus, or the transmission of viruses similar to those that establish infection within individuals. Our findings indicate that both processes are likely to be important.