Caitlin S. Pepperell

The Role of Selection in Shaping Diversity of Natural M. tuberculosis Populations

Mycobacterium tuberculosis (M.tb), the cause of tuberculosis (TB), is estimated to infect a new host every second. While analyses of genetic data from natural populations of M.tb have emphasized the role of genetic drift in shaping patterns of diversity, the influence of natural selection on this successful pathogen is less well understood. We investigated the effects of natural selection on patterns of diversity in 63 globally extant genomes of M.tb and related pathogenic mycobacteria. We found evidence of strong purifying selection, with an estimated genome-wide selection coefficient equal to −9.5×10−4 (95% CI −1.1×10−3 to −6.8×10−4); this is several orders of magnitude higher than recent estimates for eukaryotic and prokaryotic organisms. We also identified different patterns of variation across categories of gene function. Genes involved in transport and metabolism of inorganic ions exhibited very low levels of non-synonymous polymorphism, equivalent to categories under strong purifying selection (essential and translation-associated genes). The highest levels of non-synonymous variation were seen in a group of transporter genes, likely due to either diversifying selection or local selective sweeps. In addition to selection, we identified other important influences on M.tb genetic diversity, such as a 25-fold expansion of global M.tb populations coincident with explosive growth in human populations (estimated timing 1684 C.E., 95% CI 1620–1713 C.E.). These results emphasize the parallel demographic histories of this obligate pathogen and its human host, and suggest that the dominant effect of selection on M.tb is removal of novel variants, with exceptions in an interesting group of genes involved in transportation and defense. We speculate that the hostile environment within a host imposes strict demands on M.tb physiology, and thus a substantial fitness cost for most new mutations. In this respect, obligate bacterial pathogens may differ from other host-associated microbes such as symbionts.

Neurological manifestations of West Nile virus infection.

BACKGROUND Over the past four years, West Nile virus (WNV) has become a significant health issue in North America. In 2002, WNV infection made its first appearance in the human population in Canada. METHODS Patients who presented to the University Health Network and Mount Sinai Hospital in Toronto with neurological disease attributed to WNV infection were identified and followed by the neurology service. Clinical features and results of laboratory, electrodiagnostic, radiological and pathological studies are presented. RESULTS In August and September 2002, 26 patients were admitted with WNV infection; 14 presented with neurological illness. Encephalitis was the most common presentation (11 patients). Eleven patients developed neuromuscular disease; two at presentation and nine after encephalitis. While the majority had a motor process that localized to the anterior horn cell and/or motor neuron, two patients had evidence of a demyelinating neuropathy and one a sensorimotor axonal neuropathy. Less common manifestations included rhombencephalitis, ataxia, myelopathy and parkinsonism. Death occurred in four patients; two > 75 years of age, and two who were immunocompromised. CONCLUSIONS The most common neurological manifestation of WNV infection was encephalitis with subsequent neuromuscular involvement. The diversity of clinical and pathological findings, however, suggests widespread involvement of the central and peripheral nervous system. A poorer prognosis for neurological recovery and overall survival was seen in older and immunocompromised patients.

Low-Virulence Citrobacter Species Encode Resistance to Multiple Antimicrobials

ABSTRACT Citrobacter spp. are gram-negative commensal bacteria that infrequently cause serious nosocomial infections in compromised hosts. They are often resistant to cephalosporins due to overexpression of their chromosomal β-lactamase. During a recent study of multidrug-resistant Enterobacteriaceae (MDRE) in solid-organ transplant patients, we found that almost half of patients colonized with MDRE carried one or more cefpodoxime-resistant Citrobacter freundii, Citrobacter braakii, or Citrobacter amalonaticus strains. Pulsed-field gel electrophoresis showed that 36 unique strains of Citrobacter were present among 32 patients. Genetic and phenotypic analysis of the resistance mechanisms of these bacteria showed that the extended-spectrum β-lactamase (ESBL) SHV-5 or SHV-12 was encoded by 8 strains (26%) and expressed by 7 strains (19%). A number of strains were resistant to other drug classes, including aminoglycosides (28%), trimethoprim-sulfamethoxazole (31%), and fluoroquinolones (8%). PCR and DNA analysis of these multiresistant strains revealed the presence of class I integrons, including the first integrons reported for C. braakii and C. amalonaticus. The integrons encoded aminoglycoside resistance, trimethoprim resistance, or both. Despite the prevalence of MDR Citrobacter spp. in our solid-organ transplant patients, only a single infection with a colonizing strain was recorded over 18 months. Low-virulence Citrobacter spp., which can persist in the host for long periods, could influence pathogen evolution by accumulation of genes encoding resistance to multiple antimicrobial classes.

Dispersal of Mycobacterium tuberculosis via the Canadian fur trade

Patterns of gene flow can have marked effects on the evolution of populations. To better understand the migration dynamics of Mycobacterium tuberculosis, we studied genetic data from European M. tuberculosis lineages currently circulating in Aboriginal and French Canadian communities. A single M. tuberculosis lineage, characterized by the DS6Quebec genomic deletion, is at highest frequency among Aboriginal populations in Ontario, Saskatchewan, and Alberta; this bacterial lineage is also dominant among tuberculosis (TB) cases in French Canadians resident in Quebec. Substantial contact between these human populations is limited to a specific historical era (1710–1870), during which individuals from these populations met to barter furs. Statistical analyses of extant M. tuberculosis minisatellite data are consistent with Quebec as a source population for M. tuberculosis gene flow into Aboriginal populations during the fur trade era. Historical and genetic analyses suggest that tiny M. tuberculosis populations persisted for ∼100 y among indigenous populations and subsequently expanded in the late 19th century after environmental changes favoring the pathogen. Our study suggests that spread of TB can occur by two asynchronous processes: (i) dispersal of M. tuberculosis by minimal numbers of human migrants, during which small pathogen populations are sustained by ongoing migration and slow disease dynamics, and (ii) expansion of the M. tuberculosis population facilitated by shifts in host ecology. If generalizable, these migration dynamics can help explain the low DNA sequence diversity observed among isolates of M. tuberculosis and the difficulties in global elimination of tuberculosis, as small, widely dispersed pathogen populations are difficult both to detect and to eradicate.

Diversity of Mycobacterium tuberculosis across Evolutionary Scales

Tuberculosis (TB) is a global public health emergency. Increasingly drug resistant strains of Mycobacterium tuberculosis (M.tb) continue to emerge and spread, highlighting adaptability of this pathogen. Most studies of M.tb evolution have relied on ‘between-host’ samples, in which each person with TB is represented by a single M.tb isolate. However, individuals with TB commonly harbor populations of M.tb numbering in the billions. Here, we use analyses of M.tb genomic data from within and between hosts to gain insight into influences shaping genetic diversity of this pathogen. We find that the amount of M.tb genetic diversity harbored by individuals with TB can vary dramatically, likely as a function of disease severity. Surprisingly, we did not find an appreciable impact of TB treatment on M.tb diversity. In examining genomic data from M.tb samples within and between hosts with TB, we find that genes involved in the regulation, synthesis, and transportation of immunomodulatory cell envelope lipids appear repeatedly in the extremes of various statistical measures of diversity. Many of these genes have been identified as possible targets of selection in other studies employing different methods and data sets. Taken together, these observations suggest that M.tb cell envelope lipids are targets of selection within hosts. Many of these lipids are specific to pathogenic mycobacteria and, in some cases, human-pathogenic mycobacteria. We speculate that rapid adaptation of cell envelope lipids is facilitated by functional redundancy, flexibility in their metabolism, and their roles mediating interactions with the host.

Genomic Signatures of Distributive Conjugal Transfer among Mycobacteria

Distributive conjugal transfer (DCT) is a newly described mechanism of lateral gene transfer (LGT) that results in a mosaic transconjugant structure, similar to the products of meiosis. We have tested popular LGT detection methods on whole-genome sequence data from experimental DCT transconjugants and used the best performing methods to compare genomic signatures of DCT with those of LGT through natural transformation, conjugative plasmids, and mobile genetic elements (MGE). We found that DCT results in transfer of larger chromosomal segments, that these segments are distributed more broadly around the chromosome, and that a greater proportion of the chromosome is affected by DCT than by other mechanisms of LGT. We used the best performing methods to characterize LGT in Mycobacterium canettii, the mycobacterial species most closely related to Mycobacterium tuberculosis. Patterns of LGT among M. canettii were highly distinctive. Gene flow appeared unidirectional, from lineages with minimal evidence of LGT to isolates with a substantial proportion (6–13%) of sites identified as recombinant. Among M. canettii isolates with evidence of LGT, recombinant fragments were larger and more evenly distributed relative to bacteria that undergo LGT through natural transformation, conjugative plasmids, and MGE. Spatial bias in M. canettii was also unusual in that patterns of recombinant fragment sharing mirrored overall phylogenetic structure. Based on the proportion of recombinant sites, the size of recombinant fragments, their spatial distribution and lack of association with MGE, as well as unidirectionality of DNA transfer, we conclude that DCT is the predominant mechanism of LGT among M. canettii.

Armed conflict and population displacement as drivers of the evolution and dispersal of Mycobacterium tuberculosis

Significance We used population genomic analyses to reconstruct the recent history and dispersal of a major clade of Mycobacterium tuberculosis in central Asia and beyond. Our results indicate that the fall of the Soviet Union and the ensuing collapse of public health systems led to a rise in M. tuberculosis drug resistance. We also show that armed conflict and population displacement is likely to have aided the export of this clade from central Asia to war-torn Afghanistan and beyond. The “Beijing” Mycobacterium tuberculosis (Mtb) lineage 2 (L2) is spreading globally and has been associated with accelerated disease progression and increased antibiotic resistance. Here we performed a phylodynamic reconstruction of one of the L2 sublineages, the central Asian clade (CAC), which has recently spread to western Europe. We find that recent historical events have contributed to the evolution and dispersal of the CAC. Our timing estimates indicate that the clade was likely introduced to Afghanistan during the 1979–1989 Soviet–Afghan war and spread further after population displacement in the wake of the American invasion in 2001. We also find that drug resistance mutations accumulated on a massive scale in Mtb isolates from former Soviet republics after the fall of the Soviet Union, a pattern that was not observed in CAC isolates from Afghanistan. Our results underscore the detrimental effects of political instability and population displacement on tuberculosis control and demonstrate the power of phylodynamic methods in exploring bacterial evolution in space and time.

Bacterial genetic signatures of human social phenomena among M.tuberculosis from an Aboriginal Canadian population

Despite a widespread global distribution and highly variable disease phenotype, there is little DNA sequence diversity among isolates of Mycobacterium tuberculosis. In addition, many regional population genetic surveys have revealed a stereotypical structure in which a single clone, lineage, or clade makes up the majority of the population. It is often assumed that dominant clones are highly adapted, that is, the overall structure of M. tuberculosis populations is the result of positive selection. In order to test this assumption, we analyzed genetic data from extant populations of bacteria circulating in Aboriginal communities in Saskatchewan, Canada. Demographic parameters of the bacterial population were estimated from archival epidemiological data collected over approximately 130 years since the onset of epidemic tuberculosis in the host communities. Bacterial genetic data were tested against neutral theory expectations and the local evolutionary history of M. tuberculosis investigated by phylogenetic analysis. Our findings are not consistent with positive selection on the bacterial population. Instead, we uncovered founder effects persisting over decades and barriers to gene flow within the bacterial population. Simulation experiments suggested that a combination of these neutral influences could result in the stereotypical structure of M. tuberculosis populations. Some aspects of population structure were suggestive of background selection, and data were on the whole consistent with combined effects of population bottlenecks, subdivision, and background selection. Neutral phenomena, namely, bottlenecks and partitions within populations, are prominent influences on the evolution of M. tuberculosis and likely contribute to restricted genetic diversity observed within this species. Given these influences, a complex evolutionary model will be required to define the relative fitness of different M. tuberculosis lineages and, ultimately, to uncover the genetic basis for its success as a pathogen.

A molecular portrait of maternal sepsis from Byzantine Troy

Pregnancy complications are poorly represented in the archeological record, despite their importance in contemporary and ancient societies. While excavating a Byzantine cemetery in Troy, we discovered calcified abscesses among a woman’s remains. Scanning electron microscopy of the tissue revealed ‘ghost cells’, resulting from dystrophic calcification, which preserved ancient maternal, fetal and bacterial DNA of a severe infection, likely chorioamnionitis. Gardnerella vaginalis and Staphylococcus saprophyticus dominated the abscesses. Phylogenomic analyses of ancient, historical, and contemporary data showed that G. vaginalis Troy fell within contemporary genetic diversity, whereas S. saprophyticus Troy belongs to a lineage that does not appear to be commonly associated with human disease today. We speculate that the ecology of S. saprophyticus infection may have differed in the ancient world as a result of close contacts between humans and domesticated animals. These results highlight the complex and dynamic interactions with our microbial milieu that underlie severe maternal infections. DOI: http://dx.doi.org/10.7554/eLife.20983.001

Local epidemic history as a predictor of tuberculosis incidence in Saskatchewan Aboriginal communities

BACKGROUND Average tuberculosis (TB) incidence rates are high in Canadian Aboriginal communities, but there is significant variability within this group. OBJECTIVE To determine whether local history of post-contact TB epidemics is predictive of contemporary epidemiology among Aboriginal communities in Saskatchewan, Canada. METHODS TB incidence, age-specific morbidity patterns and rates of clustering of TB genotypes from 1986 to 2004 were compared between two groups of communities: Group 1, in which post-contact epidemics of TB were established around 1870, and Group 2, in which they were delayed until after 1920. Concomitant effects of socio-economic and geographic variables were explored with multivariate models. RESULTS Group 2 communities were characterized by higher annual incidence of TB (median 431 per 100,000 population vs. 38/100,000). In multivariate models that included socio-economic and geographic variables, historical grouping remained a significant independent predictor of community incidence of TB. Clustering of TB genotypes was associated with Group 2 (OR 8.7, 95%CI 3.3-22.7) and age 10-34 years (OR 2.5, 95%CI 1.1-5.7). CONCLUSIONS TB transmission dynamics can vary significantly as a function of a populations historical experience with TB. Populations at different stages along the epidemic trajectory may be amenable to different types of interventions.