Hierarchical and spatially explicit clustering of DNA sequences with BAPS software

SuperDCA for genome-wide epistasis analysis

The potential for genome-wide modelling of epistasis has recently surfaced given the possibility of sequencing densely sampled populations and the emerging families of statistical interaction models. Direct coupling analysis (DCA) has previously been shown to yield valuable predictions for single protein structures, and has recently been extended to genome-wide analysis of bacteria, identifying novel interactions in the co-evolution between resistance, virulence and core genome elements. However, earlier computational DCA methods have not been scalable to enable model fitting simultaneously to 10⁴-10⁵ polymorphisms, representing the amount of core genomic variation observed in analyses of many bacterial species. Here, we introduce a novel inference method (SuperDCA) that employs a new scoring principle, efficient parallelization, optimization and filtering on phylogenetic information to achieve scalability for up to 10⁵ polymorphisms. Using two large population samples of Streptococcus pneumoniae, we demonstrate the ability of SuperDCA to make additional significant biological findings about this major human pathogen. We also show that our method can uncover signals of selection that are not detectable by genome-wide association analysis, even though our analysis does not require phenotypic measurements. SuperDCA, thus, holds considerable potential in building understanding about numerous organisms at a systems biological level.

Defining and Evaluating a Core Genome Multilocus Sequence Typing Scheme for Genome-Wide Typing of Clostridium difficile

Clostridium difficile, recently renamed Clostridioides difficile, is the most common cause of antibiotic-associated nosocomial gastrointestinal infections worldwide. To differentiate endogenous infections and transmission events, highly discriminatory subtyping is necessary. Today, methods based on whole-genome sequencing data are increasingly used to subtype bacterial pathogens; however, frequently a standardized methodology and typing nomenclature are missing. Here we report a core genome multilocus sequence typing (cgMLST) approach developed for C. difficile Initially, we determined the breadth of the C. difficile population based on all available MLST sequence types with Bayesian inference (BAPS). The resulting BAPS partitions were used in combination with C. difficile clade information to select representative isolates that were subsequently used to define cgMLST target genes. Finally, we evaluated the novel cgMLST scheme with genomes from 3,025 isolates. BAPS grouping (n = 6 groups) together with the clade information led to a total of 11 representative isolates that were included for cgMLST definition and resulted in 2,270 cgMLST genes that were present in all isolates. Overall, 2,184 to 2,268 cgMLST targets were detected in the genome sequences of 70 outbreak-associated and reference strains, and on average 99.3% cgMLST targets (1,116 to 2,270 targets) were present in 2,954 genomes downloaded from the NCBI database, underlining the representativeness of the cgMLST scheme. Moreover, reanalyzing different cluster scenarios with cgMLST were concordant to published single nucleotide variant analyses. In conclusion, the novel cgMLST is representative for the whole C. difficile population, is highly discriminatory in outbreak situations, and provides a unique nomenclature facilitating interlaboratory exchange.

Phylogenetic Characterization of Salmonella enterica Serovar Typhimurium and Its Monophasic Variant Isolated from Food Animals in Japan Revealed Replacement of Major Epidemic Clones in the Last 4 Decades

Salmonella enterica serovar Typhimurium (Salmonella Typhimurium) and its monophasic variant (Salmonella 4,[5],12:i:-) are the major causes of gastroenteritis in both humans and animals. Pulsed-field gel electrophoresis and multilocus variable-number tandem-repeat analysis have been used widely as subtyping methods for these pathogens in molecular epidemiological analyses, but the results do not precisely reflect phylogenetic information. In this study, we performed a phylogenetic analysis of these serovars using whole-genome sequencing data and identified nine distinct genotypic clades. Then, we established an allele-specific PCR-based genotyping method detecting a clade-specific single nucleotide polymorphism to rapidly identify the clade of each isolate. Among a total of 815 isolates obtained from cattle in Japan between 1977 and 2017, clades 1, 7, and 9 contained 77% of isolates. Obvious replacement of the dominant clone was observed five times in this period, and clade 9, which mostly contains Salmonella 4,[5],12:i:-, is currently dominant. Among 140 isolates obtained from swine in Japan between 1976 and 2017, clades 3 and 9 contained 64% of isolates. Clade 9 is the latest clone as is the case in cattle isolates. Clade 9 is similar to an epidemic clone from Europe, which is characterized by sequence type 34 (ST34), chromosomal Salmonella genomic island 3, and a composite transposon containing antimicrobial resistance genes. The increased prevalence of clade 9 among food animals in Japan might be a part of the pandemic of the European Salmonella 4,[5],12:i:- clone.

Phylogenetic Characterization of Salmonella enterica Serovar Typhimurium and Its Monophasic Variant Isolated from Food Animals in Japan Revealed Replacement of Major Epidemic Clones in the Last 4 Decades

Salmonella enterica serovar Typhimurium (Salmonella Typhimurium) and its monophasic variant (Salmonella 4,[5],12:i:-) are the major causes of gastroenteritis in both humans and animals. Pulsed-field gel electrophoresis and multilocus variable-number tandem-repeat analysis have been used widely as subtyping methods for these pathogens in molecular epidemiological analyses, but the results do not precisely reflect phylogenetic information. In this study, we performed a phylogenetic analysis of these serovars using whole-genome sequencing data and identified nine distinct genotypic clades. Then, we established an allele-specific PCR-based genotyping method detecting a clade-specific single nucleotide polymorphism to rapidly identify the clade of each isolate. Among a total of 815 isolates obtained from cattle in Japan between 1977 and 2017, clades 1, 7, and 9 contained 77% of isolates. Obvious replacement of the dominant clone was observed five times in this period, and clade 9, which mostly contains Salmonella 4,[5],12:i:-, is currently dominant. Among 140 isolates obtained from swine in Japan between 1976 and 2017, clades 3 and 9 contained 64% of isolates. Clade 9 is the latest clone as is the case in cattle isolates. Clade 9 is similar to an epidemic clone from Europe, which is characterized by sequence type 34 (ST34), chromosomal Salmonella genomic island 3, and a composite transposon containing antimicrobial resistance genes. The increased prevalence of clade 9 among food animals in Japan might be a part of the pandemic of the European Salmonella 4,[5],12:i:- clone.

The impact of serotype-specific vaccination on phylodynamic parameters of Streptococcus pneumoniae and the pneumococcal pan-genome

In the United States, the introduction of the heptavalent pneumococcal conjugate vaccine (PCV) largely eliminated vaccine serotypes (VT); non-vaccine serotypes (NVT) subsequently increased in carriage and disease. Vaccination also disrupts the composition of the pneumococcal pangenome, which includes mobile genetic elements and polymorphic non-capsular antigens important for virulence, transmission, and pneumococcal ecology. Antigenic proteins are of interest for future vaccines; yet, little is known about how the they are affected by PCV use. To investigate the evolutionary impact of vaccination, we assessed recombination, evolution, and pathogen demographic history of 937 pneumococci collected from 1998-2012 among Navajo and White Mountain Apache Native American communities. We analyzed changes in the pneumococcal pangenome, focusing on metabolic loci and 19 polymorphic protein antigens. We found the impact of PCV on the pneumococcal population could be observed in reduced diversity, a smaller pangenome, and changing frequencies of accessory clusters of orthologous groups (COGs). Post-PCV7, diversity rebounded through clonal expansion of NVT lineages and inferred in-migration of two previously unobserved lineages. Accessory COGs frequencies trended toward pre-PCV7 values with increasing time since vaccine introduction. Contemporary frequencies of protein antigen variants are better predicted by pre-PCV7 values (1998-2000) than the preceding period (2006-2008), suggesting balancing selection may have acted in maintaining variant frequencies in this population. Overall, we present the largest genomic analysis of pneumococcal carriage in the United States to date, which includes a snapshot of a true vaccine-naïve community prior to the introduction of PCV7. These data improve our understanding of pneumococcal evolution and emphasize the need to consider pangenome composition when inferring the impact of vaccination and developing future protein-based pneumococcal vaccines.

Defining the Role of the Environment in the Emergence and Persistence of vanA Vancomycin-Resistant Enterococcus (VRE) in an Intensive Care Unit: A Molecular Epidemiological Study

OBJECTIVETo describe the transmission dynamics of the emergence and persistence of vanA vancomycin-resistant enterococcus (VRE) in an intensive care unit (ICU) using whole-genome sequencing of patient and environmental isolates.DESIGNRetrospective cohort study.SETTINGICU in a tertiary referral center.PARTICIPANTSPatients admitted to the ICU over an 11-month period.METHODS VanA VRE isolated from patients (n=31) were sequenced using the Illumina MiSeq platform. Environmental samples from bed spaces, equipment, and waste rooms were collected. All vanA VRE-positive environmental samples (n=14) were also sequenced. Data were collected regarding patient ward and bed movements.RESULTSThe 31 patient vanA VRE isolates were from screening (n=19), urine (n=4), bloodstream (n=3), skin/wound (n=3), and intra-abdominal (n=2) sources. The phylogeny from sequencing data confirmed several VRE clusters, with 1 group accounting for 38 of 45 isolates (84%). Within this cluster, cross-transmission was extensive and complex across the ICU. Directionality indicated that colonized patients contaminated environmental sites. Similarly, environmental sources not only led to patient colonization but also to infection. Notably, shared equipment acted as a conduit for transmission between different ICU areas. Infected patients, however, were not linked to further VRE transmission.CONCLUSIONSGenomic sequencing confirmed a predominantly clonal outbreak of VRE with complex transmission dynamics. The environmental reservoir, particularly from shared equipment, played a key role in ongoing VRE spread. This study provides evidence to support the use of multifaceted strategies, with an emphasis on measures to reduce bacterial burden in the environment, for successful VRE control.Infect Control Hosp Epidemiol 2018;39:668-675.

A fast likelihood solution to the genetic clustering problem

The investigation of genetic clusters in natural populations is an ubiquitous problem in a range of fields relying on the analysis of genetic data, such as molecular ecology, conservation biology and microbiology. Typically, genetic clusters are defined as distinct panmictic populations, or parental groups in the context of hybridisation. Two types of methods have been developed for identifying such clusters: model-based methods, which are usually computer-intensive but yield results which can be interpreted in the light of an explicit population genetic model, and geometric approaches, which are less interpretable but remarkably faster.Here, we introduce snapclust, a fast maximum-likelihood solution to the genetic clustering problem, which allies the advantages of both model-based and geometric approaches. Our method relies on maximising the likelihood of a fixed number of panmictic populations, using a combination of geometric approach and fast likelihood optimisation, using the Expectation-Maximisation (EM) algorithm. It can be used for assigning genotypes to populations and optionally identify various types of hybrids between two parental populations. Several goodness-of-fit statistics can also be used to guide the choice of the retained number of clusters.Using extensive simulations, we show that snapclust performs comparably to current gold standards for genetic clustering as well as hybrid detection, with some advantages for identifying hybrids after several backcrosses, while being orders of magnitude faster than other model-based methods. We also illustrate how snapclust can be used for identifying the optimal number of clusters, and subsequently assign individuals to various hybrid classes simulated from an empirical microsatellite dataset. snapclust is implemented in the package adegenet for the free software R, and is therefore easily integrated into existing pipelines for genetic data analysis. It can be applied to any kind of co-dominant markers, and can easily be extended to more complex models including, for instance, varying ploidy levels. Given its flexibility and computer-efficiency, it provides a useful complement to the existing toolbox for the study of genetic diversity in natural populations.

Phylogeny of the subfamily Stelliferinae suggests speciation in Ophioscion Gill, 1863 (Sciaenidae: Perciformes) in the western South Atlantic

Phylogenies based on morphological and molecular data confirm the monophyly of the subfamily Stelliferinae; however, there is no consensus on the intergeneric and interspecific relationships in the group. Previous studies suggested the non-monophyly of Ophioscion and Stellifer, and possible cryptic species in Ophioscion punctatissimus. Therefore, we used mitochondrial (16S rDNA and COI) and nuclear (Rhodopsin, EGR1, and RAG1) regions to examine phylogenetic relationships among species of this subfamily. Our results confirmed the monophyly of Stelliferinae and supports the close relationship among Bardiella, Corvula and Odontoscion, which form a sister group to Stellifer and Ophioscion. Notwithstanding, all the results support the non-monophyly of Stellifer and Ophioscion and we suggest that a taxonomic revision should consider Ophioscion as a junior synonym of Stellifer. Moreover, O. punctatissimus was grouped into two clades, with the O. punctatissimus lineage I (LI) being closer to O. scierus from the eastern Pacific than to the O. punctatissimus lineage II (LII). The most recent common ancestor (TMRCA) for the O. scierus and O. punctatissimus LI and O. punctatissimus LII clade dates from 7.2 (HPD: 4.3-10.5) Ma, whereas TMRCA for the O. scierus and O. punctatissimus LI clade dates from 5.3 (HPD: 2.4-8.6) Ma, indicating that speciation processes may be related to the rise of the Isthmus of Panama. Phylogeographic analyses corroborate the hypothesis of speciation in O. punctatissimus. These results suggest that lineages of O. punctatissimus originated from distinct ancestors and, by morphological similarity, were considered the same taxon. A taxonomic revision should be performed to validate the species status of such lineages.

Evaluation of phylogenetic reconstruction methods using bacterial whole genomes: a simulation based study

Background: Phylogenetic reconstruction is a necessary first step in many analyses which use whole genome sequence data from bacterial populations. There are many available methods to infer phylogenies, and these have various advantages and disadvantages, but few unbiased comparisons of the range of approaches have been made. Methods: We simulated data from a defined “true tree” using a realistic evolutionary model. We built phylogenies from this data using a range of methods, and compared reconstructed trees to the true tree using two measures, noting the computational time needed for different phylogenetic reconstructions. We also used real data from Streptococcus pneumoniae alignments to compare individual core gene trees to a core genome tree. Results: We found that, as expected, maximum likelihood trees from good quality alignments were the most accurate, but also the most computationally intensive. Using less accurate phylogenetic reconstruction methods, we were able to obtain results of comparable accuracy; we found that approximate results can rapidly be obtained using genetic distance based methods. In real data we found that highly conserved core genes, such as those involved in translation, gave an inaccurate tree topology, whereas genes involved in recombination events gave inaccurate branch lengths. We also show a tree-of-trees, relating the results of different phylogenetic reconstructions to each other. Conclusions: We recommend three approaches, depending on requirements for accuracy and computational time. Quicker approaches that do not perform full maximum likelihood optimisation may be useful for many analyses requiring a phylogeny, as generating a high quality input alignment is likely to be the major limiting factor of accurate tree topology. We have publicly released our simulated data and code to enable further comparisons.

Evaluation of phylogenetic reconstruction methods using bacterial whole genomes: a simulation based study

Background: Phylogenetic reconstruction is a necessary first step in many analyses which use whole genome sequence data from bacterial populations. There are many available methods to infer phylogenies, and these have various advantages and disadvantages, but few unbiased comparisons of the range of approaches have been made.

Methods: We simulated data from a defined 'true tree' using a realistic evolutionary model. We  built phylogenies from this data using a range of methods, and compared reconstructed trees to the true tree using two measures, noting the computational time needed for different phylogenetic reconstructions. We also used real data from Streptococcus pneumoniae alignments to compare individual core gene trees to a core genome tree.

Results: We found that, as expected, maximum likelihood trees from good quality alignments were the most accurate, but also the most computationally intensive. Using less accurate phylogenetic reconstruction methods, we were able to obtain results of comparable accuracy; we found that approximate results can rapidly be obtained using genetic distance based methods. In real data we found that highly conserved core genes, such as those involved in translation, gave an inaccurate tree topology, whereas genes involved in recombination events gave inaccurate branch lengths. We also show a tree-of-trees, relating the results of different phylogenetic reconstructions to each other.

Conclusions: We recommend three approaches, depending on requirements for accuracy and computational time. For the most accurate tree, use of either RAxML or IQ-TREE with an alignment of variable sites produced by mapping to a reference genome is best. Quicker approaches that do not perform full maximum likelihood optimisation may be useful for many analyses requiring a phylogeny, as generating a high quality input alignment is likely to be the major limiting factor of accurate tree topology.  We have publicly released our simulated data and code to enable further comparisons.

Infection patterns and molecular data reveal host and tissue specificity of Posthodiplostomum species in centrarchid hosts

Posthodiplostomum minimum utilizes a three-host life cycle with multiple developmental stages. The metacercarial stage, commonly known as 'white grub', infects the visceral organs of many freshwater fishes and was historically considered a host generalist due to its limited morphological variation among a wide range of hosts. In this study, infection data and molecular techniques were used to evaluate the host and tissue specificity of Posthodiplostomum metacercariae in centrarchid fishes. Eleven centrarchid species from three genera were collected from the Illinois portion of the Ohio River drainage and necropsied. Posthodiplostomum infection levels differed significantly by host age, host genera and infection locality. Three Posthodiplostomum spp. were identified by DNA sequencing, two of which were relatively common within centrarchid hosts. Both common species were host specialists at the genus level, with one species restricted to Micropterus hosts and the other preferentially infecting Lepomis. Host specificity is likely dictated by physiological compatibility and deviations from Lepomis host specificity may be related to host hybridization. Posthodiplostomum species also differed in their utilization of host tissues. Neither common species displayed strong genetic structure over the scale of this study, likely due to their utilization of bird definitive hosts.

Genomic epidemiology and population structure of Neisseria gonorrhoeae from remote highly endemic Western Australian populations

BACKGROUND

Neisseria gonorrhoeae causes gonorrhoea, the second most commonly notified sexually transmitted infection in Australia. One of the highest notification rates of gonorrhoea is found in the remote regions of Western Australia (WA). Unlike isolates from the major Australian population centres, the remote community isolates have low rates of antimicrobial resistance (AMR). Population structure and whole-genome comparison of 59 isolates from the Western Australian N. gonorrhoeae collection were used to investigate relatedness of isolates cultured in the metropolitan and remote areas. Core genome phylogeny, multilocus sequencing typing (MLST), N. gonorrhoeae multi-antigen sequence typing (NG-MAST) and N. gonorrhoeae sequence typing for antimicrobial resistance (NG-STAR) in addition to hierarchical clustering of sequences were used to characterize the isolates.

RESULTS

Population structure analysis of the 59 isolates together with 72 isolates from an international collection, revealed six population groups suggesting that N. gonorrhoeae is a weakly clonal species. Two distinct population groups, Aus1 and Aus2, represented 63% of WA isolates and were mostly composed of the remote community isolates that carried no chromosomal AMR genotypes. In contrast, the Western Australian metropolitan isolates were frequently multi-drug resistant and belonged to population groups found in the international database, suggesting international transmission of the isolates.

CONCLUSIONS

Our study suggests that the population structure of N. gonorrhoeae is distinct between the communities in remote and metropolitan WA. Given the high rate of AMR in metropolitan regions, ongoing surveillance is essential to ensure the enduring efficacy of the empiric gonorrhoea treatment in remote WA.

Molecular Evidences of a Hidden Complex Scenario in Leporinus cf. friderici

The megadiversity of the neotropical ichthyofauna has been associated to recent diversification processes, reflecting in subtle or lacking morphological differentiation between species, challenging the classical taxonomic identification. Leporinus friderici occurs in several river basins of South America, and its nominal taxonomic validity has been questioned. Its wide distribution within the Brazilian Shield suggests that this species could be genetically structured among the hydrographic basins, despite a sharp morphological similarity. In this study, we performed phylogenetic analyses, based on three nuclear (recombination activating gene 1, RAG1, recombination activating gene 2, RAG2, and myosin heavy chain 6 cardiac muscle alpha gene, Myh6) and two mitochondrial (COI and Cytochrome b, Cytb) markers, in specimens morphologically similar to L. friderici and related species from different hydrographic basins in South America. Our phylogenetic tree identified four well-supported clades, which point out to the existence of taxonomic inconsistencies within this fish group. A clade named L. cf. friderici sensu stricto included eight Molecular Operational Taxonomic Units recently diversified in the Brazilian Shield basins. These results were also confirmed by a single-gene species delimitation analysis. It is suggested that this clade includes a species complex, characterizing taxonomic uncertainties. Another clade recovered only L. friderici from the Suriname rivers, validating this nominal species in its type locality. A third no-named clade, characterized by deeper species divergence, recovered five different nominal species interleaved with other undescribed forms previously also recognized as L. cf. friderici, indicating taxonomic errors. The fourth clade only included L. taeniatus. Our results showed a complex scenario involving the morphotype L. cf. friderici and allowed us to address aspects related to evolutionary diversification of this fish group and historical processes involved with, highlighting the importance of revealing hidden biodiversity for the taxonomy and conservationist action plans of these fish.

Origin, evolution, and distribution of the molecular machinery for biosynthesis of sialylated lipooligosaccharide structures in Campylobacter coli

Campylobacter jejuni and Campylobacter coli are the most common cause of bacterial gastroenteritis worldwide. Additionally, C. jejuni is the most common bacterial etiological agent in the autoimmune Guillain-Barré syndrome (GBS). Ganglioside mimicry by C. jejuni lipooligosaccharide (LOS) is the triggering factor of the disease. LOS-associated genes involved in the synthesis and transfer of sialic acid (glycosyltranferases belonging to family GT-42) are essential in C. jejuni to synthesize ganglioside-like LOS. Despite being isolated from GBS patients, scarce genetic evidence supports C. coli role in the disease. In this study, through data mining and bioinformatics analysis, C. coli is shown to possess a larger GT-42 glycosyltransferase repertoire than C. jejuni. Although GT-42 glycosyltransferases are widely distributed in C. coli population, only a fraction of C. coli strains (1%) are very likely able to express ganglioside mimics. Even though the activity of C. coli specific GT-42 enzymes and their role in shaping the bacterial population are yet to be explored, evidence presented herein suggest that loss of function of some LOS-associated genes occurred during agriculture niche adaptation.

Translating genomics into practice for real-time surveillance and response to carbapenemase-producing Enterobacteriaceae: evidence from a complex multi-institutional KPC outbreak

Background

Until recently, Klebsiella pneumoniae carbapenemase (KPC)-producing Enterobacteriaceae were rarely identified in Australia. Following an increase in the number of incident cases across the state of Victoria, we undertook a real-time combined genomic and epidemiological investigation. The scope of this study included identifying risk factors and routes of transmission, and investigating the utility of genomics to enhance traditional field epidemiology for informing management of established widespread outbreaks.

Methods

All KPC-producing Enterobacteriaceae isolates referred to the state reference laboratory from 2012 onwards were included. Whole-genome sequencing was performed in parallel with a detailed descriptive epidemiological investigation of each case, using Illumina sequencing on each isolate. This was complemented with PacBio long-read sequencing on selected isolates to establish high-quality reference sequences and interrogate characteristics of KPC-encoding plasmids.

Results

Initial investigations indicated that the outbreak was widespread, with 86 KPC-producing Enterobacteriaceae isolates (K. pneumoniae 92%) identified from 35 different locations across metropolitan and rural Victoria between 2012 and 2015. Initial combined analyses of the epidemiological and genomic data resolved the outbreak into distinct nosocomial transmission networks, and identified healthcare facilities at the epicentre of KPC transmission. New cases were assigned to transmission networks in real-time, allowing focussed infection control efforts. PacBio sequencing confirmed a secondary transmission network arising from inter-species plasmid transmission. Insights from Bayesian transmission inference and analyses of within-host diversity informed the development of state-wide public health and infection control guidelines, including interventions such as an intensive approach to screening contacts following new case detection to minimise unrecognised colonisation.

Conclusion

A real-time combined epidemiological and genomic investigation proved critical to identifying and defining multiple transmission networks of KPC Enterobacteriaceae, while data from either investigation alone were inconclusive. The investigation was fundamental to informing infection control measures in real-time and the development of state-wide public health guidelines on carbapenemase-producing Enterobacteriaceae surveillance and management.

Within-host evolution of Enterococcus faecium during longitudinal carriage and transition to bloodstream infection in immunocompromised patients

BACKGROUND

Enterococcus faecium is a leading cause of hospital-acquired infection, particularly in the immunocompromised. Here, we use whole genome sequencing of E. faecium to study within-host evolution and the transition from gut carriage to invasive disease.

METHODS

We isolated and sequenced 180 E. faecium from four immunocompromised patients who developed bloodstream infection during longitudinal surveillance of E. faecium in stool and their immediate environment.

RESULTS

A phylogenetic tree based on single nucleotide polymorphisms (SNPs) in the core genome of the 180 isolates demonstrated several distinct clones. This was highly concordant with the population structure inferred by Bayesian methods, which contained four main BAPS (Bayesian Analysis of Population Structure) groups. The majority of isolates from each patient resided in a single group, but all four patients also carried minority populations in stool from multiple phylogenetic groups. Bloodstream isolates from each case belonged to a single BAPS group, which differed in all four patients. Analysis of 87 isolates (56 from blood) belonging to a single BAPS group that were cultured from the same patient over 54 days identified 30 SNPs in the core genome (nine intergenic, 13 non-synonymous, eight synonymous), and 250 accessory genes that were variably present. Comparison of these genetic variants in blood isolates versus those from stool or environment did not identify any variants associated with bloodstream infection. The substitution rate for these isolates was estimated to be 128 (95% confidence interval 79.82 181.77) mutations per genome per year, more than ten times higher than previous estimates for E. faecium. Within-patient variation in vancomycin resistance associated with vanA was common and could be explained by plasmid loss, or less often by transposon loss.

CONCLUSIONS

These findings demonstrate the diversity of E. faecium carriage by individual patients and significant within-host diversity of E. faecium, but do not provide evidence for adaptive genetic variation associated with invasion.

Patterns of genetic variation among geographic and host-plant associated populations of the peach fruit moth Carposina sasakii (Lepidoptera: Carposinidae)

Background

Populations of herbivorous insects may become genetically differentiated because of local adaptation to different hosts and climates as well as historical processes, and further genetic divergence may occur following the development of reproductive isolation among populations. Here we investigate the population genetic structure of the orchard pest peach fruit moth (PFM) Carposina sasakii (Lepidoptera: Carposinidae) in China, which shows distinct biological differences when characterized from different host plants. Genetic diversity and genetic structure were assessed among populations from seven plant hosts and nine regions using 19 microsatellite loci and a mitochondrial sequence.

Results

Strong genetic differentiation was found among geographical populations representing distinct geographical regions, but not in host-associated populations collected from the same area. Mantel tests based on microsatellite loci indicated an association between genetic differentiation and geographical distance, and to a lesser extent environmental differentiation. Approximate Bayesian Computation analyses supported the scenario that PFM likely originated from a southern area and dispersed northwards before the last glacial maximum during the Quaternary.

Conclusions

Our analyses suggested a strong impact of geographical barriers and historical events rather than host plants on the genetic structure of the PFM; however, uncharacterized environmental factors and host plants may also play a role. Studies on adaptive shifts in this moth should take into account geographical and historical factors.

Electronic supplementary material

The online version of this article (10.1186/s12862-017-1116-7) contains supplementary material, which is available to authorized users.

Frequency-dependent selection in vaccine-associated pneumococcal population dynamics

Many bacterial species are composed of multiple lineages distinguished by extensive variation in gene content. These often cocirculate in the same habitat, but the evolutionary and ecological processes that shape these complex populations are poorly understood. Addressing these questions is particularly important for Streptococcus pneumoniae, a nasopharyngeal commensal and respiratory pathogen, because the changes in population structure associated with the recent introduction of partial-coverage vaccines have substantially reduced pneumococcal disease. Here we show that pneumococcal lineages from multiple populations each have a distinct combination of intermediate-frequency genes. Functional analysis suggested that these loci may be subject to negative frequency-dependent selection (NFDS) through interactions with other bacteria, hosts or mobile elements. Correspondingly, these genes had similar frequencies in four populations with dissimilar lineage compositions. These frequencies were maintained following substantial alterations in lineage prevalences once vaccination programmes began. Fitting a multilocus NFDS model of post-vaccine population dynamics to three genomic datasets using Approximate Bayesian Computation generated reproducible estimates of the influence of NFDS on pneumococcal evolution, the strength of which varied between loci. Simulations replicated the stable frequency of lineages unperturbed by vaccination, patterns of serotype switching and clonal replacement. This framework highlights how bacterial ecology affects the impact of clinical interventions.

Global phylogenetic analysis of Escherichia coli and plasmids carrying the mcr-1 gene indicates bacterial diversity but plasmid restriction

To understand the dynamics behind the worldwide spread of the mcr-1 gene, we determined the population structure of Escherichia coli and of mobile genetic elements (MGEs) carrying the mcr-1 gene. After a systematic review of the literature we included 65 E. coli whole genome sequences (WGS), adding 6 recently sequenced travel related isolates, and 312 MLST profiles. We included 219 MGEs described in 7 Enterobacteriaceae species isolated from human, animal and environmental samples. Despite a high overall diversity, 2 lineages were observed in the E. coli population that may function as reservoirs of the mcr-1 gene, the largest of which was linked to ST10, a sequence type known for its ubiquity in human faecal samples and in food samples. No genotypic clustering by geographical origin or isolation source was observed. Amongst a total of 13 plasmid incompatibility types, the IncI2, IncX4 and IncHI2 plasmids accounted for more than 90% of MGEs carrying the mcr-1 gene. We observed significant geographical clustering with regional spread of IncHI2 plasmids in Europe and IncI2 in Asia. These findings point towards promiscuous spread of the mcr-1 gene by efficient horizontal gene transfer dominated by a limited number of plasmid incompatibility types.

Integrated view of Vibrio cholerae in the Americas

Latin America has experienced two of the largest cholera epidemics in modern history; one in 1991 and the other in 2010. However, confusion still surrounds the relationships between globally circulating pandemic Vibrio cholerae clones and local bacterial populations. We used whole-genome sequencing to characterize cholera across the Americas over a 40-year time span. We found that both epidemics were the result of intercontinental introductions of seventh pandemic El Tor V. cholerae and that at least seven lineages local to the Americas are associated with disease that differs epidemiologically from epidemic cholera. Our results consolidate historical accounts of pandemic cholera with data to show the importance of local lineages, presenting an integrated view of cholera that is important to the design of future disease control strategies.

Genomic history of the seventh pandemic of cholera in Africa

The seventh cholera pandemic has heavily affected Africa, although the origin and continental spread of the disease remain undefined. We used genomic data from 1070 Vibrio cholerae O1 isolates, across 45 African countries and over a 49-year period, to show that past epidemics were attributable to a single expanded lineage. This lineage was introduced at least 11 times since 1970, into two main regions, West Africa and East/Southern Africa, causing epidemics that lasted up to 28 years. The last five introductions into Africa, all from Asia, involved multidrug-resistant sublineages that replaced antibiotic-susceptible sublineages after 2000. This phylogenetic framework describes the periodicity of lineage introduction and the stable routes of cholera spread, which should inform the rational design of control measures for cholera in Africa.

Distinct Campylobacter fetus lineages adapted as livestock pathogens and human pathobionts in the intestinal microbiota

Campylobacter fetus is a venereal pathogen of cattle and sheep, and an opportunistic human pathogen. It is often assumed that C. fetus infection occurs in humans as a zoonosis through food chain transmission. Here we show that mammalian C. fetus consists of distinct evolutionary lineages, primarily associated with either human or bovine hosts. We use whole-genome phylogenetics on 182 strains from 17 countries to provide evidence that C. fetus may have originated in humans around 10,500 years ago and may have "jumped" into cattle during the livestock domestication period. We detect C. fetus genomes in 8% of healthy human fecal metagenomes, where the human-associated lineages are the dominant type (78%). Thus, our work suggests that C. fetus is an unappreciated human intestinal pathobiont likely spread by human to human transmission. This genome-based evolutionary framework will facilitate C. fetus epidemiology research and the development of improved molecular diagnostics and prevention schemes for this neglected pathogen.

Population structure of Escherichia coli O26 : H11 with recent and repeated stx2 acquisition in multiple lineages

A key virulence factor of enterohaemorrhagic Escherichia coli (EHEC) is the bacteriophage-encoded Shiga toxin (Stx). Stxs are classified into two types, Stx1 and Stx2, and Stx2-producing strains are thought to cause more severe infections than strains producing only Stx1. Although O26 : H11 is the second most prevalent EHEC following O157 : H7, the majority of O26 : H11 strains produce Stx1 alone. However, Stx2-producing O26 strains have increasingly been detected worldwide. Through a large-scale genome analysis, we present a global phylogenetic overview and evolutionary timescale for E. coli O26 : H11. The origin of O26 has been estimated to be 415 years ago. Sequence type 21C1 (ST21C1), one of the two sublineages of ST21, the most predominant O26 : H11 lineage worldwide, emerged 213 years ago from one of the three ST29 sublineages (ST29C2). The other ST21 lineage (ST21C2) emerged 95 years ago from ST21C1. Increases in population size occurred in the late 20th century for all of the O26 lineages, but most remarkably for ST21C2. Analysis of the distribution of stx2-positive strains revealed the recent and repeated acquisition of the stx2 gene in multiple lineages of O26, both in ST21 and ST29. Other major EHEC virulence genes, such as type III secretion system effector genes and plasmid-encoded virulence genes, were well conserved in ST21 compared to ST29. In addition, more antimicrobial-resistance genes have accumulated in the ST21C1 lineage. Although current attention is focused on several highly virulent ST29 clones that have acquired the stx2 gene, there is also a considerable risk that the ST21 lineage could yield highly virulent clones.

A Supervised Statistical Learning Approach for Accurate Legionella pneumophila Source Attribution during Outbreaks

Public health agencies are increasingly relying on genomics during Legionnaires' disease investigations. However, the causative bacterium (Legionella pneumophila) has an unusual population structure, with extreme temporal and spatial genome sequence conservation. Furthermore, Legionnaires' disease outbreaks can be caused by multiple L. pneumophila genotypes in a single source. These factors can confound cluster identification using standard phylogenomic methods. Here, we show that a statistical learning approach based on L. pneumophila core genome single nucleotide polymorphism (SNP) comparisons eliminates ambiguity for defining outbreak clusters and accurately predicts exposure sources for clinical cases. We illustrate the performance of our method by genome comparisons of 234 L. pneumophila isolates obtained from patients and cooling towers in Melbourne, Australia, between 1994 and 2014. This collection included one of the largest reported Legionnaires' disease outbreaks, which involved 125 cases at an aquarium. Using only sequence data from L. pneumophila cooling tower isolates and including all core genome variation, we built a multivariate model using discriminant analysis of principal components (DAPC) to find cooling tower-specific genomic signatures and then used it to predict the origin of clinical isolates. Model assignments were 93% congruent with epidemiological data, including the aquarium Legionnaires' disease outbreak and three other unrelated outbreak investigations. We applied the same approach to a recently described investigation of Legionnaires' disease within a UK hospital and observed a model predictive ability of 86%. We have developed a promising means to breach L. pneumophila genetic diversity extremes and provide objective source attribution data for outbreak investigations.IMPORTANCE Microbial outbreak investigations are moving to a paradigm where whole-genome sequencing and phylogenetic trees are used to support epidemiological investigations. It is critical that outbreak source predictions are accurate, particularly for pathogens, like Legionella pneumophila, which can spread widely and rapidly via cooling system aerosols, causing Legionnaires' disease. Here, by studying hundreds of Legionella pneumophila genomes collected over 21 years around a major Australian city, we uncovered limitations with the phylogenetic approach that could lead to a misidentification of outbreak sources. We implement instead a statistical learning technique that eliminates the ambiguity of inferring disease transmission from phylogenies. Our approach takes geolocation information and core genome variation from environmental L. pneumophila isolates to build statistical models that predict with high confidence the environmental source of clinical L. pneumophila during disease outbreaks. We show the versatility of the technique by applying it to unrelated Legionnaires' disease outbreaks in Australia and the UK.

mvmapper: Interactive spatial mapping of genetic structures

Characterizing genetic structure across geographic space is a fundamental challenge in population genetics. Multivariate statistical analyses are powerful tools for summarizing genetic variability, but geographic information and accompanying metadata are not always easily integrated into these methods in a user-friendly fashion. Here, we present a deployable Python-based web-tool, mvmapper, for visualizing and exploring results of multivariate analyses in geographic space. This tool can be used to map results of virtually any multivariate analysis of georeferenced data, and routines for exporting results from a number of standard methods have been integrated in the R package adegenet, including principal components analysis (PCA), spatial PCA, discriminant analysis of principal components, principal coordinates analysis, nonmetric dimensional scaling and correspondence analysis. mvmapper's greatest strength is facilitating dynamic and interactive exploration of the statistical and geographic frameworks side by side, a task that is difficult and time-consuming with currently available tools. Source code and deployment instructions, as well as a link to a hosted instance of mvmapper, can be found at https://popphylotools.github.io/mvMapper/.

Methicillin Resistant Staphylococcus aureus Transmission in a Ghanaian Burn Unit: The Importance of Active Surveillance in Resource-Limited Settings

Objectives:Staphylococcus aureus infections in burn patients can lead to serious complications and death. The frequency of S. aureus infection is high in low- and middle-income countries presumably due to limited resources, misuse of antibiotics and poor infection control. The objective of the present study was to apply population genomics to precisely define, for the first time, the transmission of antibiotic resistant S. aureus in a resource-limited setting in sub-Saharan Africa.

Methods:Staphylococcus aureus surveillance was performed amongst burn patients and healthcare workers during a 7-months survey within the burn unit of the Korle Bu Teaching Hospital in Ghana.

Results: Sixty-six S. aureus isolates (59 colonizing and 7 clinical) were obtained from 31 patients and 10 healthcare workers. Twenty-one of these isolates were ST250-IV methicillin-resistant S. aureus (MRSA). Notably, 25 (81%) of the 31 patients carried or were infected with S. aureus within 24 h of admission. Genome comparisons revealed six distinct S. aureus clones circulating in the burn unit, and demonstrated multiple transmission events between patients and healthcare workers. Further, the collected S. aureus isolates exhibited a wide range of genotypic resistances to antibiotics, including trimethoprim (21%), aminoglycosides (33%), oxacillin (33%), chloramphenicol (50%), tetracycline (59%) and fluoroquinolones (100%).

Conclusion: Population genomics uncovered multiple transmission events of S. aureus, especially MRSA, within the investigated burn unit. Our findings highlight lapses in infection control and prevention, and underscore the great importance of active surveillance to protect burn victims against multi-drug resistant pathogens in resource-limited settings.

Pherotype Polymorphism in Streptococcus pneumoniae Has No Obvious Effects on Population Structure and Recombination

Natural transformation in the Gram-positive pathogen Streptococcus pneumoniae occurs when cells become "competent," a state that is induced in response to high extracellular concentrations of a secreted peptide signal called competence stimulating peptide (CSP) encoded by the comC locus. Two main CSP signal types (pherotypes) are known to dominate the pherotype diversity across strains. Using 4,089 fully sequenced pneumococcal genomes, we confirm that pneumococcal populations are highly genetically structured and that there is significant variation among diverged populations in pherotype frequencies; most carry only a single pherotype. Moreover, we find that the relative frequencies of the two dominant pherotypes significantly vary within a small range across geographical sites. It has been variously proposed that pherotypes either promote genetic exchange among cells expressing the same pherotype, or conversely that they promote recombination between strains bearing different pherotypes. We attempt to distinguish these hypotheses using a bioinformatics approach by estimating recombination frequencies within and between pherotypes across 4,089 full genomes. Despite underlying population structure, we observe extensive recombination between populations; additionally, we found significantly higher (although marginal) rates of genetic exchange between strains expressing different pherotypes than among isolates carrying the same pherotype. Our results indicate that pherotypes do not restrict, and may even slightly facilitate, recombination between strains; however, these marginal effects suggest the more likely possibility that the cause of CSP polymorphism lies outside of its effects on transformation. Our results suggest that the CSP balanced polymorphism does not causally underlie population differentiation. Therefore, when strains carrying different pherotypes encounter one another during cocolonization, genetic exchange can occur without restriction.

Population structuring of the ubiquitous stingless bee Tetragonisca angustula in southern Brazil as revealed by microsatellite and mitochondrial markers

Tetragonisca angustula is one of the most widespread stingless bees in the Neotropics. This species swarms frequently and is extremely successful in urban environments. In addition, it is one of the most popular stingless bee species for beekeeping in Latin America, so nest transportation and trading is common. Nest transportation can change the genetic structure of the host population, reducing inbreeding and increasing homogenization. Here, we evaluate the genetic structure of 17 geographic populations of T. angustula in southern Brazil to quantify the level of genetic differentiation between populations. Analyses were conducted on partially sequenced mitochondrial genes and 11 microsatellite loci of 1002 workers from 457 sites distributed on the mainland and on 3 islands. Our results show that T. angustula populations are highly differentiated as demonstrated by mitochondrial DNA (mtDNA) and microsatellite markers. Of 73 haplotypes, 67 were population-specific. MtDNA diversity was low in 9 populations but microsatellite diversity was moderate to high in all populations. Microsatellite data suggest 10 genetic clusters and low level of gene flow throughout the studied area. However, physical barriers, such as rivers and mountain ranges, or the presence or absence of forest appear to be unrelated to population clusters. Factors such as low dispersal, different ecological conditions, and isolation by distance are most likely shaping the population structure of this species. Thus far, nest transportation has not influenced the general population structure in the studied area. However, due to the genetic structure we found, we recommend that nest transportation should only occur within and between populations that are genetically similar.

Rapid Gene Turnover as a Significant Source of Genetic Variation in a Recently Seeded Population of a Healthcare-Associated Pathogen

Genome sequencing has been useful to gain an understanding of bacterial evolution. It has been used for studying the phylogeography and/or the impact of mutation and recombination on bacterial populations. However, it has rarely been used to study gene turnover at microevolutionary scales. Here, we sequenced Mexican strains of the human pathogen Acinetobacter baumannii sampled from the same locale over a 3 year period to obtain insights into the microevolutionary dynamics of gene content variability. We found that the Mexican A. baumannii population was recently founded and has been emerging due to a rapid clonal expansion. Furthermore, we noticed that on average the Mexican strains differed from each other by over 300 genes and, notably, this gene content variation has accrued more frequently and faster than the accumulation of mutations. Moreover, due to its rapid pace, gene content variation reflects the phylogeny only at very short periods of time. Additionally, we found that the external branches of the phylogeny had almost 100 more genes than the internal branches. All in all, these results show that rapid gene turnover has been of paramount importance in producing genetic variation within this population and demonstrate the utility of genome sequencing to study alternative forms of genetic variation.

Evidence for a primate origin of zoonotic Helicobacter suis colonizing domesticated pigs

Helicobacter suis is the second most prevalent Helicobacter species in the stomach of humans suffering from gastric disease. This bacterium mainly inhabits the stomach of domesticated pigs, in which it causes gastric disease, but it appears to be absent in wild boars. Interestingly, it also colonizes the stomach of asymptomatic rhesus and cynomolgus monkeys. The origin of modern human-, pig- or non-human primate-associated H. suis strains in these respective host populations was hitherto unknown. Here we show that H. suis in pigs possibly originates from non-human primates. Our data suggest that a host jump from macaques to pigs happened between 100 000 and 15 000 years ago and that pig domestication has had a significant impact on the spread of H. suis in the pig population, from where this pathogen occasionally infects humans. Thus, in contrast to our expectations, H. suis appears to have evolved in its main host in a completely different way than its close relative Helicobacter pylori in humans.

Comparative landscape genetics of pond-breeding amphibians in Mediterranean temporal wetlands: The positive role of structural heterogeneity in promoting gene flow

Comparative landscape genetics studies can provide key information to implement cost-effective conservation measures favouring a broad set of taxa. These studies are scarce, particularly in Mediterranean areas, which include diverse but threatened biological communities. Here, we focus on Mediterranean wetlands in central Iberia and perform a multi-level, comparative study of two endemic pond-breeding amphibians, a salamander (Pleurodeles waltl) and a toad (Pelobates cultripes). We genotyped 411 salamanders from 20 populations and 306 toads from 16 populations at 18 and 16 microsatellite loci, respectively, and identified major factors associated with population connectivity through the analysis of three sets of variables potentially affecting gene flow at increasingly finer levels of spatial resolution. Topographic, land use/cover, and remotely sensed vegetation/moisture indices were used to derive optimized resistance surfaces for the two species. We found contrasting patterns of genetic structure, with stronger, finer scale genetic differentiation in Pleurodeles waltl, and notable differences in the role of fine-scale patterns of heterogeneity in vegetation cover and water content in shaping patterns of regional genetic structure in the two species. Overall, our results suggest a positive role of structural heterogeneity in population connectivity in pond-breeding amphibians, with habitat patches of Mediterranean scrubland and open oak woodlands ("dehesas") facilitating gene flow. Our study highlights the usefulness of remotely sensed continuous variables of land cover, vegetation and water content (e.g., NDVI, NDMI) in conservation-oriented studies aimed at identifying major drivers of population connectivity.

Lack of genetic structure in pinworm populations from New World primates in forest fragments

Microevolutionary processes in parasites are driven by factors related to parasite biology, host abundance and dispersal, and environmental conditions. Here, we test the prediction that isolation of host populations results in reduced genetic diversity and high differentiation among parasite populations. We conducted a population genetic analysis of two pinworms, Trypanoxyuris minutus and Trypanoxyuris atelis, commonly found parasitizing howler and spider monkeys in tropical rainforests across south-eastern Mexico, whose populations are currently isolated due to anthropogenic habitat loss and fragmentation. Mitochondrial DNA was employed to assess parasite genetic patterns, as well as to analyse their demography and population history. Both pinworm species showed high haplotype diversity but, unexpectedly, lower nucleotide diversity than that reported for other parasites. No genetic differentiation or population structure was detected in either pinworm species despite habitat loss, fragmentation and host isolation. Several scenarios are discussed that could help to explain the genetic panmixia found in both pinworm species, including higher than expected primate inter-fragment dispersal movements, and passive dispersal facilitating gene flow between parasite populations. The results suggest that large population sizes of parasites could be helping them to cope with the isolation and fragmentation of populations, delaying the effects of genetic drift. The present study highlights the complexity of the drivers that intervene in the evolutionary processes of parasites. Detailed genetic studies are needed, both in host and parasite populations, to assess the effects that habitat perturbation and environmental changes could have on the evolutionary dynamics of pinworms and primates.

Mechanisms of peripheral phylogeographic divergence in the indo-Pacific: lessons from the spiny lobster Panulirus homarus

Background

There is increasing recognition of the concordance between marine biogeographic and phylogeographic boundaries. However, it is still unclear how population-level divergence translates into species-level divergence, and what are the principal factors that first initiate that divergence, and then maintain reproductive isolation. This study examines the likely forces driving population and lineage divergences in the broadly-distributed Indo-Pacific spiny lobster Panulirus homarus, which has peripheral divergent lineages in the west and east. The study focuses particularly on the West Indian Ocean, which is emerging as a region of unexpected diversity. Mitochondrial control region (mtCR) and COI sequences as well as genotypes of 9 microsatellite loci were examined in 410 individuals from 17 locations grouped into 7 regions from South Africa in the west, and eastward across to Taiwan and the Marquesas Islands. Phylogenetic and population-level analyses were used to test the significance and timing of divergences and describe the genetic relationships among populations.

Results

Analyses of the mtCR revealed high levels of divergence among the seven regions (Ф_ST = 0.594, P < 0.001). Microsatellite analyses also revealed significant divergence among regions, but at a much lower level (F_ST = 0.066, P < 0.001). The results reveal different patterns of mtCR v. nDNA divergence between the two distinct peripheral lineages: a subspecies in South Africa and Madagascar, and a phylogeographically diverged population in the Marquesas. The results also expose a number of other more fine-scale population divergences, particularly in the Indian Ocean.

Conclusions

The divergence of peripheral lineages in the west and east of the species’ range appear to have been initiated and maintained by very different processes. The pattern of mitochondrial and nuclear divergence of the western lineage, implicates processes of parapatric isolation, secondary contact and introgression, and suggests possible maintenance through adaptation and behavioural reproductive isolation. In contrast, the eastern lineage appears to have diverged through a rare colonisation event, maintained through long-term isolation, and matches expectations of the core-periphery hypothesis. The process of active peripheral speciation may be a common force in the Indo-Pacific that helps drive some of the regions’ recognized biogeographic boundaries.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1050-8) contains supplementary material, which is available to authorized users.

Population genetic structure, antibiotic resistance, capsule switching and evolution of invasive pneumococci before conjugate vaccination in Malawi

INTRODUCTION

Pneumococcal infections cause a high death toll in Sub Saharan Africa (SSA) but the recently rolled out pneumococcal conjugate vaccines (PCV) will reduce the disease burden. To better understand the population impact of these vaccines, comprehensive analysis of large collections of pneumococcal isolates sampled prior to vaccination is required. Here we present a population genomic study of the invasive pneumococcal isolates sampled before the implementation of PCV13 in Malawi.

MATERIALS AND METHODS

We retrospectively sampled and whole genome sequenced 585 invasive isolates from 2004 to 2010. We determine the pneumococcal population genetic structure and assessed serotype prevalence, antibiotic resistance rates, and the occurrence of serotype switching.

RESULTS

Population structure analysis revealed 22 genetically distinct sequence clusters (SCs), which consisted of closely related isolates. Serotype 1 (ST217), a vaccine-associated serotype in clade SC2, showed highest prevalence (19.3%), and was associated with the highest MDR rate (81.9%) followed by serotype 12F, a non-vaccine serotype in clade SC10 with an MDR rate of 57.9%. Prevalence of serotypes was stable prior to vaccination although there was an increase in the PMEN19 clone, serotype 5 ST289, in clade SC1 in 2010 suggesting a potential undetected local outbreak. Coalescent analysis revealed recent emergence of the SCs and there was evidence of natural capsule switching in the absence of vaccine induced selection pressure. Furthermore, majority of the highly prevalent capsule-switched isolates were associated with acquisition of vaccine-targeted capsules.

CONCLUSIONS

This study provides descriptions of capsule-switched serotypes and serotypes with potential to cause serotype replacement post-vaccination such as 12F. Continued surveillance is critical to monitor these serotypes and antibiotic resistance in order to design better infection prevention and control measures such as inclusion of emerging replacement serotypes in future conjugate vaccines.

Comparative Genomics of Bacillus thuringiensis Reveals a Path to Specialized Exploitation of Multiple Invertebrate Hosts

Understanding the genetic basis of host shifts is a key genomic question for pathogen and parasite biology. The Bacillus cereus group, which encompasses Bacillus thuringiensis and Bacillus anthracis, contains pathogens that can infect insects, nematodes, and vertebrates. Since the target range of the essential virulence factors (Cry toxins) and many isolates is well known, this group presents a powerful system for investigating how pathogens can diversify and adapt to phylogenetically distant hosts. Specialization to exploit insects occurs at the level of the major clade and is associated with substantial changes in the core genome, and host switching between insect orders has occurred repeatedly within subclades. The transfer of plasmids with linked cry genes may account for much of the adaptation to particular insect orders, and network analysis implies that host specialization has produced strong associations between key toxin genes with similar targets. Analysis of the distribution of plasmid minireplicons shows that plasmids with orf156 and orf157, which carry genes encoding toxins against Lepidoptera or Diptera, were contained only by B. thuringiensis in the specialized insect clade (clade 2), indicating that tight genome/plasmid associations have been important in adaptation to invertebrate hosts. Moreover, the accumulation of multiple virulence factors on transposable elements suggests that cotransfer of diverse virulence factors is advantageous in terms of expanding the insecticidal spectrum, overcoming insect resistance, or through gains in pathogenicity via synergistic interactions between toxins.IMPORTANCE Population genomics have provided many new insights into the formation, evolution, and dynamics of bacterial pathogens of humans and other higher animals, but these pathogens usually have very narrow host ranges. As a pathogen of insects and nematodes, Bacillus thuringiensis, which produces toxins showing toxicity to many orders of insects and other invertebrates, can be used as a model to study the evolution of pathogens with wide host ranges. Phylogenomic analysis revealed that host specialization and switching occur at the level of the major clade and subclade, respectively. A toxin gene co-occurrence network indicates that multiple toxins with similar targets were accumulated by the same cell in the whole species. This accumulation may be one of the strategies that B. thuringiensis has used to fight against host resistance. This kind of formation and evolution of pathogens represents a different path used against multiple invertebrate hosts from that used against higher animals.

Pan-genome and phylogeny of Bacillus cereus sensu lato

BACKGROUND

Bacillus cereus sensu lato (s. l.) is an ecologically diverse bacterial group of medical and agricultural significance. In this study, I use publicly available genomes and novel bioinformatic workflows to characterize the B. cereus s. l. pan-genome and perform the largest phylogenetic and population genetic analyses of this group to date in terms of the number of genes and taxa included. With these fundamental data in hand, I identify genes associated with particular phenotypic traits (i.e., "pan-GWAS" analysis), and quantify the degree to which taxa sharing common attributes are phylogenetically clustered.

METHODS

A rapid k-mer based approach (Mash) was used to create reduced representations of selected Bacillus genomes, and a fast distance-based phylogenetic analysis of this data (FastME) was performed to determine which species should be included in B. cereus s. l. The complete genomes of eight B. cereus s. l. species were annotated de novo with Prokka, and these annotations were used by Roary to produce the B. cereus s. l. pan-genome. Scoary was used to associate gene presence and absence patterns with various phenotypes. The orthologous protein sequence clusters produced by Roary were filtered and used to build HaMStR databases of gene models that were used in turn to construct phylogenetic data matrices. Phylogenetic analyses used RAxML, DendroPy, ClonalFrameML, PAUP*, and SplitsTree. Bayesian model-based population genetic analysis assigned taxa to clusters using hierBAPS. The genealogical sorting index was used to quantify the phylogenetic clustering of taxa sharing common attributes.

RESULTS

The B. cereus s. l. pan-genome currently consists of ≈60,000 genes, ≈600 of which are "core" (common to at least 99% of taxa sampled). Pan-GWAS analysis revealed genes associated with phenotypes such as isolation source, oxygen requirement, and ability to cause diseases such as anthrax or food poisoning. Extensive phylogenetic analyses using an unprecedented amount of data produced phylogenies that were largely concordant with each other and with previous studies. Phylogenetic support as measured by bootstrap probabilities increased markedly when all suitable pan-genome data was included in phylogenetic analyses, as opposed to when only core genes were used. Bayesian population genetic analysis recommended subdividing the three major clades of B. cereus s. l. into nine clusters. Taxa sharing common traits and species designations exhibited varying degrees of phylogenetic clustering.

CONCLUSIONS

All phylogenetic analyses recapitulated two previously used classification systems, and taxa were consistently assigned to the same major clade and group. By including accessory genes from the pan-genome in the phylogenetic analyses, I produced an exceptionally well-supported phylogeny of 114 complete B. cereus s. l. genomes. The best-performing methods were used to produce a phylogeny of all 498 publicly available B. cereus s. l. genomes, which was in turn used to compare three different classification systems and to test the monophyly status of various B. cereus s. l. species. The majority of the methodology used in this study is generic and could be leveraged to produce pan-genome estimates and similarly robust phylogenetic hypotheses for other bacterial groups.

Systematic longitudinal survey of invasive Escherichia coli in England demonstrates a stable population structure only transiently disturbed by the emergence of ST131

Escherichia coli associated with urinary tract infections and bacteremia has been intensively investigated, including recent work focusing on the virulent, globally disseminated, multidrug-resistant lineage ST131. To contextualize ST131 within the broader E. coli population associated with disease, we used genomics to analyze a systematic 11-yr hospital-based survey of E. coli associated with bacteremia using isolates collected from across England by the British Society for Antimicrobial Chemotherapy and from the Cambridge University Hospitals NHS Foundation Trust. Population dynamics analysis of the most successful lineages identified the emergence of ST131 and ST69 and their establishment as two of the five most common lineages along with ST73, ST95, and ST12. The most frequently identified lineage was ST73. Compared to ST131, ST73 was susceptible to most antibiotics, indicating that multidrug resistance was not the dominant reason for prevalence of E. coli lineages in this population. Temporal phylogenetic analysis of the emergence of ST69 and ST131 identified differences in the dynamics of emergence and showed that expansion of ST131 in this population was not driven by sequential emergence of increasingly resistant subclades. We showed that over time, the E. coli population was only transiently disturbed by the introduction of new lineages before a new equilibrium was rapidly achieved. Together, these findings suggest that the frequency of E. coli lineages in invasive disease is driven by negative frequency-dependent selection occurring outside of the hospital, most probably in the commensal niche, and that drug resistance is not a primary determinant of success in this niche.

Spatial and temporal genetic structure at the fourth trophic level in a fragmented landscape

A fragmented habitat becomes increasingly fragmented for species at higher trophic levels, such as parasitoids. To persist, these species are expected to possess life-history traits, such as high dispersal, that facilitate their ability to use resources that become scarce in fragmented landscapes. If a specialized parasitoid disperses widely to take advantage of a sparse host, then the parasitoid population should have lower genetic structure than the host. We investigated the temporal and spatial genetic structure of a hyperparasitoid (fourth trophic level) in a fragmented landscape over 50 × 70 km, using microsatellite markers, and compared it with the known structures of its host parasitoid, and the butterfly host which lives as a classic metapopulation. We found that population genetic structure decreases with increasing trophic level. The hyperparasitoid has fewer genetic clusters (K = 4), than its host parasitoid (K = 15), which in turn is less structured than the host butterfly (K = 27). The genetic structure of the hyperparasitoid also shows temporal variation, with genetic differentiation increasing due to reduction of the population size, which reduces the effective population size. Overall, our study confirms the idea that specialized species must be dispersive to use a fragmented host resource, but that this adaptation has limits.

Evolution of the Staphylococcus argenteus ST2250 Clone in Northeastern Thailand Is Linked with the Acquisition of Livestock-Associated Staphylococcal Genes

Staphylococcus argenteus is a newly named species previously described as a divergent lineage of Staphylococcus aureus that has recently been shown to have a global distribution. Despite growing evidence of the clinical importance of this species, knowledge about its population epidemiology and genomic architecture is limited. We used whole-genome sequencing to evaluate and compare S. aureus (n = 251) and S. argenteus (n = 68) isolates from adults with staphylococcal sepsis at several hospitals in northeastern Thailand between 2006 and 2013. The majority (82%) of the S. argenteus isolates were of multilocus sequence type 2250 (ST2250). S. aureus was more diverse, although 43% of the isolates belonged to ST121. Bayesian analysis suggested an S. argenteus ST2250 substitution rate of 4.66 (95% confidence interval [CI], 3.12 to 6.38) mutations per genome per year, which was comparable to the S. aureus ST121 substitution rate of 4.07 (95% CI, 2.61 to 5.55). S. argenteus ST2250 emerged in Thailand an estimated 15 years ago, which contrasts with the S. aureus ST1, ST88, and ST121 clades that emerged around 100 to 150 years ago. Comparison of S. argenteus ST2250 genomes from Thailand and a global collection indicated a single introduction into Thailand, followed by transmission to local and more distant countries in Southeast Asia and further afield. S. argenteus and S. aureus shared around half of their core gene repertoire, indicating a high level of divergence and providing strong support for their classification as separate species. Several gene clusters were present in ST2250 isolates but absent from the other S. argenteus and S. aureus study isolates. These included multiple exotoxins and antibiotic resistance genes that have been linked previously with livestock-associated S. aureus, consistent with a livestock reservoir for S. argenteus. These genes appeared to be associated with plasmids and mobile genetic elements and may have contributed to the biological success of ST2250.

In this study, we used whole-genome sequencing to understand the genome evolution and population structure of a systematic collection of ST2250 S. argenteus isolates. A newly identified ancestral species of S. aureus, S. argenteus has become increasingly known as a clinically important species that has been reported recently across various countries. Our results indicate that S. argenteus has spread at a relatively rapid pace over the past 2 decades across northeastern Thailand and acquired multiple exotoxin and antibiotic resistance genes that have been linked previously with livestock-associated S. aureus. Our findings highlight the clinical importance and potential pathogenicity of S. argenteus as a recently emerging pathogen.

Pre-vaccine serotype composition within a lineage signposts its serotype replacement - a carriage study over 7 years following pneumococcal conjugate vaccine use in the UK

Serotype replacement has been reported in carriage and disease after pneumococcal conjugate vaccine (PCV) introductions in the UK and globally. We previously described concurrent expansion and decline of sequence types associated with serotype replacement over 5 years following PCV introductions in the UK. Here we use whole-genome sequencing to fully characterise the population structure of pneumococcal isolates collected over seven winters encompassing PCV7 and PCV13 introductions in the UK, investigating the importance of lineages in serotype replacement. We analysed 672 pneumococcal genomes from colonised children of 4 years old or less. The temporal prevalence of 20 lineages, defined by hierarchical Bayesian analysis of population structure (BAPS), was assessed in the context of serotype replacement. Multiple serotypes were detected in the primary winter of sampling within three vaccine-type (VT) lineages BAPS4, BAPS10 and BAPS11, in which serotype replacement were observed. In contrast, serotype replacement was not seen in the remaining three VT lineages (BAPS1, BAPS13 and BAPS14), that expressed a single serotype (6B, 6A and 3, respectively) in the primary winter. One lineage, BAPS1 serotype 6B was undetectable in the population towards the end of the study period. The dynamics of serotype replacement, in this UK population, was preceded by the presence or absence of multiple serotypes within VT lineages, in the pre-PCV population. This observation could help predict which non-vaccine types (NVTs) may be involved in replacement in future PCV introductions here and elsewhere. It could further indicate whether any antibiotic resistance associated with the lineages is likely to be affected by replacement.

Dynamics and impact of homologous recombination on the evolution of Legionella pneumophila

Legionella pneumophila is an environmental bacterium and the causative agent of Legionnaires' disease. Previous genomic studies have shown that recombination accounts for a high proportion (>96%) of diversity within several major disease-associated sequence types (STs) of L. pneumophila. This suggests that recombination represents a potentially important force shaping adaptation and virulence. Despite this, little is known about the biological effects of recombination in L. pneumophila, particularly with regards to homologous recombination (whereby genes are replaced with alternative allelic variants). Using newly available population genomic data, we have disentangled events arising from homologous and non-homologous recombination in six major disease-associated STs of L. pneumophila (subsp. pneumophila), and subsequently performed a detailed characterisation of the dynamics and impact of homologous recombination. We identified genomic "hotspots" of homologous recombination that include regions containing outer membrane proteins, the lipopolysaccharide (LPS) region and Dot/Icm effectors, which provide interesting clues to the selection pressures faced by L. pneumophila. Inference of the origin of the recombined regions showed that isolates have most frequently imported DNA from isolates belonging to their own clade, but also occasionally from other major clades of the same subspecies. This supports the hypothesis that the possibility for horizontal exchange of new adaptations between major clades of the subspecies may have been a critical factor in the recent emergence of several clinically important STs from diverse genomic backgrounds. However, acquisition of recombined regions from another subspecies, L. pneumophila subsp. fraseri, was rarely observed, suggesting the existence of a recombination barrier and/or the possibility of ongoing speciation between the two subspecies. Finally, we suggest that multi-fragment recombination may occur in L. pneumophila, whereby multiple non-contiguous segments that originate from the same molecule of donor DNA are imported into a recipient genome during a single episode of recombination.

Emergence and spread of a human-transmissible multidrug-resistant nontuberculous mycobacterium

Lung infections with Mycobacterium abscessus, a species of multidrug-resistant nontuberculous mycobacteria, are emerging as an important global threat to individuals with cystic fibrosis (CF), in whom M. abscessus accelerates inflammatory lung damage, leading to increased morbidity and mortality. Previously, M. abscessus was thought to be independently acquired by susceptible individuals from the environment. However, using whole-genome analysis of a global collection of clinical isolates, we show that the majority of M. abscessus infections are acquired through transmission, potentially via fomites and aerosols, of recently emerged dominant circulating clones that have spread globally. We demonstrate that these clones are associated with worse clinical outcomes, show increased virulence in cell-based and mouse infection models, and thus represent an urgent international infection challenge.

Indication for selfing in geographically separated populations and evidence for Pleistocene survival within the Alps: the case of Cylindrus obtusus (Pulmonata: Helicidae)

Background

Cylindrus obtusus is one of the most prominent endemic snail species of the Eastern Alps. It is restricted to alpine meadows and calcareous rocky habitats above 1500 m. Peculiar intraspecific differences have been observed in its genital tract in the eastern populations the two mucus glands associated with the love dart sac are highly variable, while almost no variation was observed in the western populations. This raises the question whether the mode and success of reproduction of the respective populations are different. To find out whether these anatomical differences reflect genetic differentiation, which might be an indication for distinct glacial refugia, we investigated a 650 bp fragment of the mitochondrial cytochrome oxidase subunit 1 gene (COI) (280 individuals) and 9 microsatellite loci from samples (487 individuals from 29 populations) covering the whole distribution range of the species.

Results

The COI sequences show a geographic differentiation between eastern, central and western populations. The westernmost localities, which were covered under ice sheets during glacial periods, are characterized by extreme low variability. Overall genetic distances among all individuals are small (max. 1.7% p-distance). The microsatellite analysis reveals a high differentiation between populations, implying restriction of gene flow. The highest genetic variability was found in the central populations. Remarkably, nearly all individuals from the eastern populations, which are more variable in their genital morphology, are homozygous in all microsatellite loci, although different alleles were found within populations.

Conclusions

The most peculiar outcome of the study is the strong evidence for selfing in C. obtusus as indicated by the microsatellite data in the easternmost populations. This finding is supported by the deformation of the mucus glands in the same populations. Since mucus glands play an important role in sexual reproduction, it seems plausible that in selfing organisms these structures are reduced. The phylogeographic structure revealed by COI sequences implies that the species has survived the ice ages within the Calcareous Alps. The small genetic distances among all individuals (max. 1.7%) suggest that C. obtusus has experienced severe bottlenecks in the past.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0977-0) contains supplementary material, which is available to authorized users.

Peripatric speciation in an endemic Macaronesian plant after recent divergence from a widespread relative

The Macaronesian Scrophularia lowei is hypothesized to have arisen from the widespread S. arguta on the basis of several phylogenetic studies of the genus, but sampling has been limited. Although these two annual species are morphologically distinct, the origin of S. lowei is unclear because genetic studies focused on this Macaronesian species are lacking. We studied 5 S. lowei and 25 S. arguta populations to determine the relationship of both species and to infer the geographical origin of S. lowei. The timing of S. lowei divergence and differentiation was inferred by dating analysis of the ITS region. A phylogenetic analysis of two nuclear (ITS and ETS) and two chloroplast (psbJ–petA and psbA–trnH) DNA regions was performed to study the relationship between the two species, and genetic differentiation was analysed by AMOVA. Haplotype network construction and Bayesian phylogeographic analysis were conducted using chloroplast DNA regions and a spatial clustering analysis was carried out on a combined dataset of all studied regions. Our results indicate that both species constitute a well-supported clade that diverged in the Miocene and differentiated in the Late Miocene-Pleistocene. Although S. lowei constitutes a well-supported clade according to nDNA, cpDNA revealed a close relationship between S. lowei and western Canarian S. arguta, a finding supported by the spatial clustering analysis. Both species have strong population structure, with most genetic variability explained by inter-population differences. Our study therefore supports a recent peripatric speciation of S. lowei—a taxon that differs morphologically and genetically at the nDNA level from its closest relative, S. arguta, but not according to cpDNA, from the closest Macaronesian populations of that species. In addition, a recent dispersal of S. arguta to Madeira from Canary Islands or Selvagens Islands and a rapid morphological differentiation after the colonization to generate S. lowei is the most likely hypothesis to explain the origin of the last taxon.

Whole genome sequencing reveals within-host genetic changes in paired meningococcal carriage isolates from Ethiopia

Background

Meningococcal colonization is a prerequisite for transmission and disease, but the bacterium only very infrequently causes disease while asymptomatic carriage is common. Carriage is highly dynamic, showing a great variety across time and space within and across populations, but also within individuals. The understanding of genetic changes in the meningococcus during carriage, when the bacteria resides in its natural niche, is important for understanding not only the carriage state, but the dynamics of the entire meningococcal population.

Results

Paired meningococcal isolates, obtained from 50 asymptomatic carriers about 2 months apart were analyzed with whole genome sequencing (WGS). Phylogenetic analysis revealed that most paired isolates from the same individual were closely related, and the average and median number of allelic differences between paired isolates defined as the same strain was 35. About twice as many differences were seen between isolates from different individuals within the same sequence type (ST). In 8%, different strains were detected at different time points. A difference in ST was observed in 6%, including an individual who was found to carry three different STs over the course of 9 weeks. One individual carried different strains from the same ST.

In total, 566 of 1605 cgMLST genes had undergone within-host genetic changes in one or more pairs. The most frequently changed cgMLST gene was relA that was changed in 47% of pairs. Across the whole genome, pilE, differed mostly, in 85% of the pairs. The most frequent mechanisms of genetic difference between paired isolates were phase variation and recombination, including gene conversion. Different STs showed variation with regard to which genes that were most frequently changed, mostly due to absence/presence of phase variation.

Conclusions

This study revealed within-host genetic differences in meningococcal isolates during short-term asymptomatic carriage. The most frequently changed genes were genes belonging to the pilin family, the restriction/modification system, opacity proteins and genes involved in glycosylation. Higher resolution genome-wide sequence typing is necessary to resolve the diversity of isolates and reveals genetic differences not discovered by traditional typing schemes, and would be the preferred choice of technology.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3806-3) contains supplementary material, which is available to authorized users.

Tracing the Enterococci from Paleozoic Origins to the Hospital

We examined the evolutionary history of leading multidrug resistant hospital pathogens, the enterococci, to their origin hundreds of millions of years ago. Our goal was to understand why, among the vast diversity of gut flora, enterococci are so well adapted to the modern hospital environment. Molecular clock estimation, together with analysis of their environmental distribution, phenotypic diversity, and concordance with host fossil records, place the origins of the enterococci around the time of animal terrestrialization, 425-500 mya. Speciation appears to parallel the diversification of hosts, including the rapid emergence of new enterococcal species following the End Permian Extinction. Major drivers of speciation include changing carbohydrate availability in the host gut. Life on land would have selected for the precise traits that now allow pathogenic enterococci to survive desiccation, starvation, and disinfection in the modern hospital, foreordaining their emergence as leading hospital pathogens.

Recombination-Mediated Host Adaptation by Avian Staphylococcus aureus

Staphylococcus aureus are globally disseminated among farmed chickens causing skeletal muscle infections, dermatitis, and septicaemia. The emergence of poultry-associated lineages has involved zoonotic transmission from humans to chickens but questions remain about the specific adaptations that promote proliferation of chicken pathogens. We characterized genetic variation in a population of genome-sequenced S. aureus isolates of poultry and human origin. Genealogical analysis identified a dominant poultry-associated sequence cluster within the CC5 clonal complex. Poultry and human CC5 isolates were significantly distinct from each other and more recombination events were detected in the poultry isolates. We identified 44 recombination events in 33 genes along the branch extending to the poultry-specific CC5 cluster, and 47 genes were found more often in CC5 poultry isolates compared with those from humans. Many of these gene sequences were common in chicken isolates from other clonal complexes suggesting horizontal gene transfer among poultry associated lineages. Consistent with functional predictions for putative poultry-associated genes, poultry isolates showed enhanced growth at 42 °C and greater erythrocyte lysis on chicken blood agar in comparison with human isolates. By combining phenotype information with evolutionary analyses of staphylococcal genomes, we provide evidence of adaptation, following a human-to-poultry host transition. This has important implications for the emergence and dissemination of new pathogenic clones associated with modern agriculture.