A High-Resolution View of Genome-Wide Pneumococcal Transformation

Transformation is an important mechanism of microbial evolution through which bacteria have been observed to rapidly adapt in response to clinical interventions; examples include facilitating vaccine evasion and the development of penicillin resistance in the major respiratory pathogen Streptococcus pneumoniae. To characterise the process in detail, the genomes of 124 S. pneumoniae isolates produced through in vitro transformation were sequenced and recombination events detected. Those recombinations importing the selected marker were independent of unselected events elsewhere in the genome, the positions of which were not significantly affected by local sequence similarity between donor and recipient or mismatch repair processes. However, both types of recombinations were sometimes mosaic, with multiple non-contiguous segments originating from the same molecule of donor DNA. The lengths of the unselected events were exponentially distributed with a mean of 2.3 kb, implying that recombinations are stochastically resolved with a fixed per base probability of 4.4×10⁻⁴ bp⁻¹. This distribution of recombination sizes, coupled with an observed under representation of large insertions within transferred sequence, suggests transformation has the potential to reduce the size of bacterial genomes, and is unlikely to act as an efficient mechanism for the uptake of accessory genomic loci.     

Pangenome Analysis of Burkholderia pseudomallei: Genome Evolution Preserves Gene Order despite High Recombination Rates

The pangenomic diversity in Burkholderia pseudomallei is high, with approximately 5.8% of the genome consisting of genomic islands. Genomic islands are known hotspots for recombination driven primarily by site-specific recombination associated with tRNAs. However, recombination rates in other portions of the genome are also high, a feature we expected to disrupt gene order. We analyzed the pangenome of 37 isolates of B. pseudomallei and demonstrate that the pangenome is ‘open’, with approximately 136 new genes identified with each new genome sequenced, and that the global core genome consists of 4568±16 homologs. Genes associated with metabolism were statistically overrepresented in the core genome, and genes associated with mobile elements, disease, and motility were primarily associated with accessory portions of the pangenome. The frequency distribution of genes present in between 1 and 37 of the genomes analyzed matches well with a model of genome evolution in which 96% of the genome has very low recombination rates but 4% of the genome recombines readily. Using homologous genes among pairs of genomes, we found that gene order was highly conserved among strains, despite the high recombination rates previously observed. High rates of gene transfer and recombination are incompatible with retaining gene order unless these processes are either highly localized to specific sites within the genome, or are characterized by symmetrical gene gain and loss. Our results demonstrate that both processes occur: localized recombination introduces many new genes at relatively few sites, and recombination throughout the genome generates the novel multi-locus sequence types previously observed while preserving gene order.     

Negative frequency-dependent selection is intensified at higher population densities in protist populations

Natural populations of free-living protists often exhibit high-levels of intraspecific diversity, yet this is puzzling as classic evolutionary theory predicts dominance by genotypes with high fitness, particularly in large populations where selection is efficient. Here, we test whether negative frequency-dependent selection (NFDS) plays a role in the maintenance of diversity in the marine flagellate Oxyrrhis marina using competition experiments between multiple pairs of strains. We observed strain-specific responses to frequency and density, but an overall signature of NFDS that was intensified at higher population densities. Because our strains were not selected a priori on the basis of particular traits expected to exhibit NFDS, these data represent a relatively unbiased estimate of the role for NFDS in maintaining diversity in protist populations. These findings could help to explain how bloom-forming plankton, which periodically achieve exceptionally high population densities, maintain substantial intraspecific diversity.     

Specificity and Strain-Typing Capabilities of Nanorod Array-Surface Enhanced Raman Spectroscopy for Mycoplasma pneumoniae Detection

Mycoplasma pneumoniae is a cell wall-less bacterial pathogen of the human respiratory tract that accounts for > 20% of all community-acquired pneumonia (CAP). At present the most effective means for detection and strain-typing is quantitative polymerase chain reaction (qPCR), which can exhibit excellent sensitivity and specificity but requires separate tests for detection and genotyping, lacks standardization between available tests and between labs, and has limited practicality for widespread, point-of-care use. We have developed and previously described a silver nanorod array-surface enhanced Raman Spectroscopy (NA-SERS) biosensing platform capable of detecting M. pneumoniae with statistically significant specificity and sensitivity in simulated and true clinical throat swab samples, and the ability to distinguish between reference strains of the two main genotypes of M. pneumoniae. Furthermore, we have established a qualitative lower endpoint of detection for NA-SERS of < 1 genome equivalent (cell/μl) and a quantitative multivariate detection limit of 5.3 ± 1 cells/μl. Here we demonstrate using partial least squares- discriminatory analysis (PLS-DA) of sample spectra that NA-SERS correctly identified M. pneumoniae clinical isolates from globally diverse origins and distinguished these from a panel of 12 other human commensal and pathogenic mycoplasma species with 100% cross-validated statistical accuracy. Furthermore, PLS-DA correctly classified by strain type all 30 clinical isolates with 96% cross-validated accuracy for type 1 strains, 98% cross-validated accuracy for type 2 strains, and 90% cross-validated accuracy for type 2V strains.     

Comparative genomic analyses of seventeen Streptococcus pneumoniae strains: insights into the pneumococcal supragenome

The distributed-genome hypothesis (DGH) states that pathogenic bacteria possess a supragenome that is much larger than the genome of any single bacterium and that these pathogens utilize genetic recombination and a large, noncore set of genes as a means of diversity generation. We sequenced the genomes of eight nasopharyngeal strains of Streptococcus pneumoniae isolated from pediatric patients with upper respiratory symptoms and performed quantitative genomic analyses among these and nine publicly available pneumococcal strains. Coding sequences from all strains were grouped into 3,170 orthologous gene clusters, of which 1,454 (46%) were conserved among all 17 strains. The majority of the gene clusters, 1,716 (54%), were not found in all strains. Genic differences per strain pair ranged from 35 to 629 orthologous clusters, with each strain's genome containing between 21 and 32% noncore genes. The distribution of the orthologous clusters per genome for the 17 strains was entered into the finite-supragenome model, which predicted that (i) the S. pneumoniae supragenome contains more than 5,000 orthologous clusters and (ii) 99% of the orthologous clusters (~3,000) that are represented in the S. pneumoniae population at frequencies of ≥0.1 can be identified if 33 representative genomes are sequenced. These extensive genic diversity data support the DGH and provide a basis for understanding the great differences in clinical phenotype associated with various pneumococcal strains. When these findings are taken together with previous studies that demonstrated the presence of a supragenome for Streptococcus agalactiae and Haemophilus influenzae, it appears that the possession of a distributed genome is a common host interaction strategy.     

High genetic diversity at the extreme range edge: nucleotide variation at nuclear loci in Scots pine (Pinus sylvestris L.) in Scotland

Nucleotide polymorphism at 12 nuclear loci was studied in Scots pine populations across an environmental gradient in Scotland, to evaluate the impacts of demographic history and selection on genetic diversity. At eight loci, diversity patterns were compared between Scottish and continental European populations. At these loci, a similar level of diversity (θ_sil=∼0.01) was found in Scottish vs mainland European populations, contrary to expectations for recent colonization, however, less rapid decay of linkage disequilibrium was observed in the former (ρ=0.0086±0.0009, ρ=0.0245±0.0022, respectively). Scottish populations also showed a deficit of rare nucleotide variants (multi-locus Tajima''s D=0.316 vs D=−0.379) and differed significantly from mainland populations in allelic frequency and/or haplotype structure at several loci. Within Scotland, western populations showed slightly reduced nucleotide diversity (π_tot=0.0068) compared with those from the south and east (0.0079 and 0.0083, respectively) and about three times higher recombination to diversity ratio (ρ/θ=0.71 vs 0.15 and 0.18, respectively). By comparison with results from coalescent simulations, the observed allelic frequency spectrum in the western populations was compatible with a relatively recent bottleneck (0.00175 × 4N_e generations) that reduced the population to about 2% of the present size. However, heterogeneity in the allelic frequency distribution among geographical regions in Scotland suggests that subsequent admixture of populations with different demographic histories may also have played a role.     

Classification of Bartonella Strains Associated with Straw-Colored Fruit Bats (Eidolon helvum) across Africa Using a Multi-locus Sequence Typing Platform

Bartonellae are facultative intracellular bacteria and are highly adapted to their mammalian host cell niches. Straw-colored fruit bats (Eidolon helvum) are commonly infected with several bartonella strains. To elucidate the genetic diversity of these bartonella strains, we analyzed 79 bartonella isolates from straw-colored fruit bats in seven countries across Africa (Cameroon, Annobon island of Equatorial Guinea, Ghana, Kenya, Nigeria, Tanzania, and Uganda) using a multi-locus sequencing typing (MLST) approach based on nucleotide sequences of eight loci (ftsZ, gltA, nuoG, ribC, rpoB, ssrA, ITS, and 16S rRNA). The analysis of each locus but ribC demonstrated clustering of the isolates into six genogroups (E1 – E5 and Ew), while ribC was absent in the isolates belonging to the genogroup Ew. In general, grouping of all isolates by each locus was mutually supportive; however, nuoG, gltA, and rpoB showed some incongruity with other loci in several strains, suggesting a possibility of recombination events, which were confirmed by network analyses and recombination/mutation rate ratio (r/m) estimations. The MLST scheme revealed 45 unique sequence types (ST1 – 45) among the analyzed bartonella isolates. Phylogenetic analysis of concatenated sequences supported the discrimination of six phylogenetic lineages (E1 – E5 and Ew) corresponding to separate and unique Bartonella species. One of the defined lineages, Ew, consisted of only two STs (ST1 and ST2), and comprised more than one-quarter of the analyzed isolates, while other lineages contained higher numbers of STs with a smaller number of isolates belonging to each lineage. The low number of allelic polymorphisms of isolates belonging to Ew suggests a more recent origin for this species. Our findings suggest that at least six Bartonella species are associated with straw-colored fruit bats, and that distinct STs can be found across the distribution of this bat species, including in populations of bats which are genetically distinct.     

Genomic Epidemiology of an Endoscope-Associated Outbreak of Klebsiella pneumoniae Carbapenemase (KPC)-Producing K. pneumoniae

Increased incidence of infections due to Klebsiella pneumoniae carbapenemase (KPC)-producing Klebsiella pneumoniae (KPC-Kp) was noted among patients undergoing endoscopic retrograde cholangiopancreatography (ERCP) at a single hospital. An epidemiologic investigation identified KPC-Kp and non-KPC-producing, extended-spectrum β-lactamase (ESBL)-producing Kp in cultures from 2 endoscopes. Genotyping was performed on patient and endoscope isolates to characterize the microbial genomics of the outbreak. Genetic similarity of 51 Kp isolates from 37 patients and 3 endoscopes was assessed by pulsed-field gel electrophoresis (PFGE) and multi-locus sequence typing (MLST). Five patient and 2 endoscope isolates underwent whole genome sequencing (WGS). Two KPC-encoding plasmids were characterized by single molecule, real-time sequencing. Plasmid diversity was assessed by endonuclease digestion. Genomic and epidemiologic data were used in conjunction to investigate the outbreak source. Two clusters of Kp patient isolates were genetically related to endoscope isolates by PFGE. A subset of patient isolates were collected post-ERCP, suggesting ERCP endoscopes as a possible source. A phylogeny of 7 Kp genomes from patient and endoscope isolates supported ERCP as a potential source of transmission. Differences in gene content defined 5 ST258 subclades and identified 2 of the subclades as outbreak-associated. A novel KPC-encoding plasmid, pKp28 helped define and track one endoscope-associated ST258 subclade. WGS demonstrated high genetic relatedness of patient and ERCP endoscope isolates suggesting ERCP-associated transmission of ST258 KPC-Kp. Gene and plasmid content discriminated the outbreak from endemic ST258 populations and assisted with the molecular epidemiologic investigation of an extended KPC-Kp outbreak.     

The acquisition of clinically relevant amoxicillin resistance in Streptococcus pneumoniae requires ordered horizontal gene transfer of four loci

Understanding how antimicrobial resistance spreads is critical for optimal application of new treatments. In the naturally competent human pathogen Streptococcus pneumoniae, resistance to β-lactam antibiotics is mediated by recombination events in genes encoding the target proteins, resulting in reduced drug binding affinity. However, for the front-line antibiotic amoxicillin, the exact mechanism of resistance still needs to be elucidated. Through successive rounds of transformation with genomic DNA from a clinically resistant isolate, we followed amoxicillin resistance development. Using whole genome sequencing, we showed that multiple recombination events occurred at different loci during one round of transformation. We found examples of non-contiguous recombination, and demonstrated that this could occur either through multiple D-loop formation from one donor DNA molecule, or by the integration of multiple DNA fragments. We also show that the final minimum inhibitory concentration (MIC) differs depending on recipient genome, explained by differences in the extent of recombination at key loci. Finally, through back transformations of mutant alleles and fluorescently labelled penicillin (bocillin-FL) binding assays, we confirm that pbp1a, pbp2b, pbp2x, and murM are the main resistance determinants for amoxicillin resistance, and that the order of allele uptake is important for successful resistance evolution. We conclude that recombination events are complex, and that this complexity contributes to the highly diverse genotypes of amoxicillin-resistant pneumococcal isolates.     

Identification of Zoonotic Genotypes of Giardia duodenalis

Giardia duodenalis, originally regarded as a commensal organism, is the etiologic agent of giardiasis, a gastrointestinal disease of humans and animals. Giardiasis causes major public and veterinary health concerns worldwide. Transmission is either direct, through the faecal-oral route, or indirect, through ingestion of contaminated water or food. Genetic characterization of G. duodenalis isolates has revealed the existence of seven groups (assemblages A to G) which differ in their host distribution. Assemblages A and B are found in humans and in many other mammals, but the role of animals in the epidemiology of human infection is still unclear, despite the fact that the zoonotic potential of Giardia was recognised by the WHO some 30 years ago. Here, we performed an extensive genetic characterization of 978 human and 1440 animal isolates, which together comprise 3886 sequences from 4 genetic loci. The data were assembled into a molecular epidemiological database developed by a European network of public and veterinary health Institutions. Genotyping was performed at different levels of resolution (single and multiple loci on the same dataset). The zoonotic potential of both assemblages A and B is evident when studied at the level of assemblages, sub-assemblages, and even at each single locus. However, when genotypes are defined using a multi-locus sequence typing scheme, only 2 multi-locus genotypes (MLG) of assemblage A and none of assemblage B appear to have a zoonotic potential. Surprisingly, mixtures of genotypes in individual isolates were repeatedly observed. Possible explanations are the uptake of genetically different Giardia cysts by a host, or subsequent infection of an already infected host, likely without overt symptoms, with a different Giardia species, which may cause disease. Other explanations for mixed genotypes, particularly for assemblage B, are substantial allelic sequence heterogeneity and/or genetic recombination. Although the zoonotic potential of G. duodenalis is evident, evidence on the contribution and frequency is (still) lacking. This newly developed molecular database has the potential to tackle intricate epidemiological questions concerning protozoan diseases.     

Rapid Global Expansion of the Fungal Disease Chytridiomycosis into Declining and Healthy Amphibian Populations

The fungal disease chytridiomycosis, caused by Batrachochytrium dendrobatidis, is enigmatic because it occurs globally in both declining and apparently healthy (non-declining) amphibian populations. This distribution has fueled debate concerning whether, in sites where it has recently been found, the pathogen was introduced or is endemic. In this study, we addressed the molecular population genetics of a global collection of fungal strains from both declining and healthy amphibian populations using DNA sequence variation from 17 nuclear loci and a large fragment from the mitochondrial genome. We found a low rate of DNA polymorphism, with only two sequence alleles detected at each locus, but a high diversity of diploid genotypes. Half of the loci displayed an excess of heterozygous genotypes, consistent with a primarily clonal mode of reproduction. Despite the absence of obvious sex, genotypic diversity was high (44 unique genotypes out of 59 strains). We provide evidence that the observed genotypic variation can be generated by loss of heterozygosity through mitotic recombination. One strain isolated from a bullfrog possessed as much allelic diversity as the entire global sample, suggesting the current epidemic can be traced back to the outbreak of a single clonal lineage. These data are consistent with the current chytridiomycosis epidemic resulting from a novel pathogen undergoing a rapid and recent range expansion. The widespread occurrence of the same lineage in both healthy and declining populations suggests that the outcome of the disease is contingent on environmental factors and host resistance.     

Multi-locus inference of population structure: a comparison between single nucleotide polymorphisms and microsatellites

Although growing numbers of single nucleotide polymorphisms (SNPs) and microsatellites (short tandem repeat polymorphisms or STRPs) are used to infer population structure, their relative properties in this context remain poorly understood. SNPs and STRPs mutate differently, suggesting multi-locus genotypes at these loci might differ in ability to detect population structure. Here, we use coalescent simulations to measure the power of sets of SNPs and STRPs to identify population structure. To maximize the applicability of our results to empirical studies, we focus on the popular STRUCTURE analysis and evaluate the role of several biological and practical factors in the detection of population structure. We find that: (1) fewer unlinked STRPs than SNPs are needed to detect structure at recent divergence times <0.3 N_e generations; (2) accurate estimation of the number of populations requires many fewer STRPs than SNPs; (3) for both marker types, declines in power due to modest gene flow (N_em=1.0) are largely negated by increasing marker number; (4) variation in the STRP mutational model affects power modestly; (5) SNP haplotypes (θ=1, no recombination) provide power comparable with STRP loci (θ=10); (6) ascertainment schemes that select highly variable STRP or SNP loci increase power to detect structure, though ascertained data may not be suitable to other inference; and (7) when samples are drawn from an admixed population and one of its parent populations, the reduction in power to detect two populations is greater for STRPs than SNPs. These results should assist the design of multi-locus studies to detect population structure in nature.     

Genetic Barriers to Historical Gene Flow between Cryptic Species of Alpine Bumblebees Revealed by Comparative Population Genomics

Evidence is accumulating that gene flow commonly occurs between recently diverged species, despite the existence of barriers to gene flow in their genomes. However, we still know little about what regions of the genome become barriers to gene flow and how such barriers form. Here, we compare genetic differentiation across the genomes of bumblebee species living in sympatry and allopatry to reveal the potential impact of gene flow during species divergence and uncover genetic barrier loci. We first compared the genomes of the alpine bumblebee Bombus sylvicola and a previously unidentified sister species living in sympatry in the Rocky Mountains, revealing prominent islands of elevated genetic divergence in the genome that colocalize with centromeres and regions of low recombination. This same pattern is observed between the genomes of another pair of closely related species living in allopatry (B. bifarius and B. vancouverensis). Strikingly however, the genomic islands exhibit significantly elevated absolute divergence (d_XY) in the sympatric, but not the allopatric, comparison indicating that they contain loci that have acted as barriers to historical gene flow in sympatry. Our results suggest that intrinsic barriers to gene flow between species may often accumulate in regions of low recombination and near centromeres through processes such as genetic hitchhiking, and that divergence in these regions is accentuated in the presence of gene flow.     

Recombination and linkage disequilibrium among mitochondrial genes in structured populations of the gynodioecious plant Silene vulgaris

The impact of intergenic recombination on the population genetics of plant mitochondrial genomes is unknown. In an effort to study this in the gynodioecious plant Silene vulgaris three-locus PCR/RFLP genotypes (based on the mitochondrial genes atpA, cox1, and cob) were determined for 239 individuals collected from 20 North American populations. Seventeen three-locus PCR/RFLP genotypes were found. Recombination was indicated by observation of each of the four two-locus genotypes possible when the two most common alleles are considered for each of two loci. Based on these common alleles the absolute values of standardized linkage disequilibrium |D'| between pairs of loci range from 0.17 to 0.78. This indicates modest disequilibrium, rather than the maximum value expected in the absence of recombination |D'=1|, or the linkage equilibrium expected if recombination is pervasive (D'=0). Values of D' did not depend on which pair of loci contributed alleles to the analysis. The direction of D' obtained for the common atpA and cox1 alleles was comparable in sign and magnitude to that obtained by examining similar information obtained in a prior study of European samples. All three loci indicated a high degree of population structure (average FST=0.63), which would limit the within-population genetic diversity required for intergenic recombination to create novel genotypes, if most mating is local. Thus, population structure acts as a constraint on the approach to linkage equilibrium.     

Haldane,Waddington and recombinant inbred lines: extension of their work to any number of genes

In the early 1930s, J. B. S. Haldane and C. H. Waddington collaborated on the consequences of genetic linkage and inbreeding. One elegant mathematical genetics problem solved by them concerns recombinant inbred lines (RILs) produced via repeated self or brother–sister mating. In this classic contribution, Haldane and Waddington derived an analytical formula for the probabilities of 2-locus and 3-locus RIL genotypes. Specifically, the Haldane–Waddington formula gives the recombination rate R in such lines as a simple function of the per generation recombination rate r. Interestingly, for more than 80 years, an extension of this result to four or more loci remained elusive. In 2015, we generalized the Haldane–Waddington self-mating result to any number of loci. Our solution used self-consistent equations of the multi-locus probabilities ‘for an infinite number of generations’ and solved these by simple algebraic operations. In practice, our approach provides a quantum leap in the systems that can be handled: the cases of up to six loci can be solved by hand while a computer program implementing our mathematical formalism tackles up to 20 loci on standard desktop computers.     

A recombination suppressor contributes to ecological speciation in OSTRINIA moths

Despite unparalleled access to species'' genomes in our post-genomic age, we often lack adequate biological explanations for a major hallmark of the speciation process—genetic divergence. In the presence of gene flow, chromosomal rearrangements such as inversions are thought to promote divergence and facilitate speciation by suppressing recombination. Using a combination of genetic crosses, phenotyping of a trait underlying ecological isolation, and population genetic analysis of wild populations, we set out to determine whether evidence supports a role for recombination suppressors during speciation between the Z and E strains of European corn borer moth (Ostrinia nubilalis). Our results are consistent with the presence of an inversion that has contributed to accumulation of ecologically adaptive alleles and genetic differentiation across roughly 20% of the Ostrinia sex chromosome (~4 Mb). Patterns in Ostrinia suggest that chromosomal divergence may involve two separate phases—one driving its transient origin through local adaptation and one determining its stable persistence through differential introgression. As the evolutionary rate of rearrangements in lepidopteran genomes appears to be one of the fastest among eukaryotes, structural mutations may have had a disproportionate role during adaptive divergence and speciation in Ostrinia and in other moths and butterflies.