Search

The Gardnerella genus, comprising at least 13 species, is associated with the polymicrobial disorder bacterial vaginosis (BV). However, the details of BV pathogenesis are poorly defined, and the contributions made by individual species, including Gardnerella spp., are largely unknown. We report here that colony phenotypes characterized by size (large and small) and opacity (opaque and translucent) are phase variable and are conserved among all tested Gardnerella strains, representing at least 10 different species. With the hypothesis that these different variants could be an important missing piece to the enigma of how BV develops in vivo , we characterized their phenotypic, proteomic, and genomic differences. Beyond increased colony size, large colony variants showed reduced vaginolysin secretion and faster growth rate relative to small colony variants. The ability to inhibit the growth of Neisseria gonorrhoeae and commensal Lactobacillus species varied by strain and, in some instances, differed between variants. Proteomics analyses indicated that 127-173 proteins were differentially expressed between variants. Proteins with increased expression in large variants of both strains were associated with amino acid and protein synthesis and protein folding, whereas those increased in small variants were related to nucleotide synthesis, phosphate transport, ABC transport, and glycogen breakdown. Furthermore, whole genome sequencing analyses revealed an abundance of genes associated with variable homopolymer tracts, implicating slipped strand mispairing in Gardnerella phase variation and illuminating the potential for previously unrecognized heterogeneity within clonal populations. Collectively, these results suggest that phase variants may be primed to serve different roles in BV pathogenesis.IMPORTANCEBacterial vaginosis is the most common gynecological disorder in women of childbearing age. Gardnerella species are crucial to the development of this dysbiosis, but the mechanisms involved in the infection are not understood. We discovered that Gardnerella species vary between two different forms, reflected in bacterial colony size. A slow-growing form makes large amounts of the toxin vaginolysin and is better able to survive in human cervix tissue. A fast-growing form is likely the one that proliferates to high numbers just prior to symptom onset and forms the biofilm that serves as a scaffold for multiple BV-associated anaerobic bacteria. Identification of the proteins that vary between different forms of the bacteria as well as those that vary randomly provides insight into the factors important for Gardnerella infection and immune avoidance., Competing Interests: The authors declare no conflict of interest.

Mycobacterium tuberculosis (M. tb), the causative agent of tuberculosis (TB), is a leading global cause of death from infectious disease. Biofilms are increasingly recognized as a relevant growth form during M. tb infection and may impede treatment by enabling bacterial drug and immune tolerance. M. tb has a complicated regulatory network that has been well-characterized for many relevant disease states, including dormancy and hypoxia. However, despite its importance, our knowledge of the genes and pathways involved in biofilm formation is limited. Here we characterize the biofilm transcriptomes of fully virulent clinical isolates and find that the regulatory systems underlying biofilm growth vary widely between strains and are also distinct from regulatory programs associated with other environmental cues. We used experimental evolution to investigate changes to the transcriptome during adaptation to biofilm growth and found that the application of a uniform selection pressure resulted in loss of strain-to-strain variation in gene expression, resulting in a more uniform biofilm transcriptome. The adaptive trajectories of transcriptomes were shaped by the genetic background of the M. tb population leading to convergence on a sub-lineage specific transcriptome. We identified widespread upregulation of non-coding RNA (ncRNA) as a common feature of the biofilm transcriptome and hypothesize that ncRNA function in genome-wide modulation of gene expression, thereby facilitating rapid regulatory responses to new environments. These results reveal a new facet of the M. tb regulatory system and provide valuable insight into how M. tb adapts to new environments., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Youngblom et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Severe outbreaks and deaths have been linked to the emergence and global spread of fluoroquinolone-resistant Clostridioides difficile over the past two decades. At the same time, metronidazole, a nitro-containing antibiotic, has shown decreasing clinical efficacy in treating C. difficile infection (CDI). Most metronidazole-resistant C. difficile exhibit an unusual resistance phenotype that can only be detected in susceptibility tests using molecularly intact heme. Here, we describe the mechanism underlying this trait. We find that most metronidazole-resistant C. difficile strains carry a T-to-G mutation (which we term PnimB G ) in the promoter of gene nimB, resulting in constitutive transcription. Silencing or deleting nimB eliminates metronidazole resistance. NimB is related to Nim proteins that are known to confer resistance to nitroimidazoles. We show that NimB is a heme-dependent flavin enzyme that degrades nitroimidazoles to amines lacking antimicrobial activity. Furthermore, occurrence of the PnimB G mutation is associated with a Thr82Ile substitution in DNA gyrase that confers fluoroquinolone resistance in epidemic strains. Our findings suggest that the pandemic of fluoroquinolone-resistant C. difficile occurring over the past few decades has also been characterized by widespread resistance to metronidazole., (© 2023. The Author(s).)

The incidence of gonorrhoea is increasing at an alarming pace, and therapeutic options continue to narrow as a result of worsening drug resistance. Neisseria gonorrhoeae is naturally competent, allowing the organism to adapt rapidly to selection pressures including antibiotics. A sub-population of N. gonorrhoeae carries the Gonococcal Genetic Island (GGI), which encodes a type IV secretion system (T4SS) that secretes chromosomal DNA. Previous research has shown that the GGI increases transformation efficiency in vitro , but the extent to which it contributes to horizontal gene transfer (HGT) during infection is unknown. Here we analysed genomic data from clinical isolates of N. gonorrhoeae to better characterize GGI+ and GGI- sub-populations and to delineate patterns of variation at the locus itself. We found the element segregating at an intermediate frequency (61%), and it appears to act as a mobile genetic element with examples of gain, loss, exchange and intra-locus recombination within our sample. We further found evidence suggesting that GGI+ and GGI- sub-populations preferentially inhabit distinct niches with different opportunities for HGT. Previously, GGI+ isolates were reported to be associated with more severe clinical infections, and our results suggest this could be related to metal-ion trafficking and biofilm formation. The co-segregation of GGI+ and GGI- isolates despite mobility of the element suggests that both niches inhabited by N. gonorrhoeae remain important to its overall persistence as has been demonstrated previously for cervical- and urethral-adapted sub-populations. These data emphasize the complex population structure of N. gonorrhoeae and its capacity to adapt to diverse niches.

The viridans group streptococci (VGS) are a large consortium of commensal streptococci that colonize the human body. Many species within this group are opportunistic pathogens causing bacteremia and infective endocarditis (IE), yet little is known about why some strains cause invasive disease. Identification of virulence determinants is complicated by the difficulty of distinguishing between the closely related species of this group. Here, we analyzed genomic data from VGS that were isolated from blood cultures in patients with invasive infections and oral swabs of healthy volunteers and then determined the best-performing methods for species identification. Using whole-genome sequence data, we characterized the population structure of a diverse sample of Streptococcus oralis isolates and found evidence of frequent recombination. We used multiple genome-wide association study tools to identify candidate determinants of invasiveness. These tools gave consistent results, leading to the discovery of a single synonymous single nucleotide polymorphism (SNP) that was significantly associated with invasiveness. This SNP was within a previously undescribed gene that was conserved across the majority of VGS species. Using the growth in the presence of human serum and a simulated infective endocarditis vegetation model, we were unable to identify a phenotype for the enriched allele in laboratory assays, suggesting a phenotype may be specific to natural infection. These data highlighted the power of analyzing natural populations for gaining insight into pathogenicity, particularly for organisms with complex population structures like the VGS. IMPORTANCE The viridians group streptococci (VGS) are a large collection of closely related commensal streptococci, with many being opportunistic pathogens causing invasive diseases, such as bacteremia and infective endocarditis. Little is known about virulence determinants in these species, and there is a distinct lack of genomic information available for the VGS. In this study, we collected VGS isolates from invasive infections and healthy volunteers and performed whole-genome sequencing for a suite of downstream analyses. We focused on a diverse sample of Streptococcus oralis genomes and identified high rates of recombination in the population as well as a single genome variant highly enriched in invasive isolates. The variant lies within a previously uncharacterized gene, nrdM , which shared homology with the anaerobic ribonucleoside triphosphate reductase, nrdD , and was highly conserved among VGS. This work increased our knowledge of VGS genomics and indicated that differences in virulence potential among S. oralis isolates were, at least in part, genetically determined.

In this issue of Cell, Saelens et al. describe a new function for mycobacterial Type VII secretion systems: manipulation of host cell migration. They find that a substantial proportion of global TB cases arise from bacteria lacking this function, raising questions about its role in pathoadaptation of Mycobacterium tuberculosis., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Mycobacterium tuberculosis is a globally distributed, lethal pathogen of humans. The virulence armamentarium of M. tuberculosis appears to have been developed on a scaffold of antiphagocytic defenses found among diverse, mostly free-living species of Mycobacterium . Pathoadaptation was further aided by the modularity, flexibility, and interactivity characterizing mycobacterial effectors and their regulators. During emergence of M. tuberculosis , novel genetic material was acquired, created, and integrated with existing tools. The major mutational mechanisms underlying these adaptations are discussed in this review, with examples. During its evolution, M. tuberculosis lost the ability and/or opportunity to engage in lateral gene transfer, but despite this it has retained the adaptability that characterizes mycobacteria. M. tuberculosis exemplifies the evolutionary genomic mechanisms underlying adoption of the pathogenic niche, and studies of its evolution have uncovered a rich array of discoveries about how new pathogens are made.

Tuberculosis (TB), caused by Mycobacterium tuberculosis ( M. tb ), is a leading cause of death due to infectious disease. TB is not traditionally associated with biofilms, but M. tb biofilms are linked with drug and immune tolerance and there is increasing recognition of their contribution to the recalcitrance of TB infections. Here, we used M. tb experimental evolution to investigate this complex phenotype and identify candidate loci controlling biofilm formation. We identified novel candidate loci, adding to our understanding of the genetic architecture underlying M. tb biofilm development. Under selective pressure to grow as a biofilm, regulatory mutations rapidly swept to fixation and were associated with changes in multiple traits, including extracellular matrix production, cell size, and growth rate. Genetic and phenotypic paths to enhanced biofilm growth varied according to the genetic background of the parent strain, suggesting that epistatic interactions are important in M. tb adaptation to changing environments., Competing Interests: TS, JK, MM, SF, LB, TM, MO The authors declare no conflict of interes, MY, CP No competing interests declared, (© 2022, Smith, Youngblom et al.)

Mycobacterium abscessus is a rapid growing, free-living species of bacterium that also causes lung infections in humans. Human infections are usually acquired from the environment; however, dominant circulating clones (DCCs) have emerged recently in both M. abscessus subsp. massiliense and subsp. abscessus that appear to be transmitted among humans and are now globally distributed. These recently emerged clones are potentially informative about the ecological and evolutionary mechanisms of pathogen emergence and host adaptation. The geographical distribution of DCCs has been reported, but the genomic processes underlying their transition from environmental bacterium to human pathogen are not well characterized. To address this knowledge gap, we delineated the structure of M. abscessus subspecies abscessus and massiliense using genomic data from 200 clinical isolates of M. abscessus from seven geographical regions. We identified differences in overall patterns of lateral gene transfer (LGT) and barriers to LGT between subspecies and between environmental and host-adapted bacteria. We further characterized genome reorganization that accompanied bacterial host adaptation, inferring selection pressures acting at both genic and intergenic loci. We found that both subspecies encode an expansive pangenome with many genes at rare frequencies. Recombination appears more frequent in M. abscessus subsp. massiliense than in subsp. abscessus , consistent with prior reports. We found evidence suggesting that phage are exchanged between subspecies, despite genetic barriers evident elsewhere throughout the genome. Patterns of LGT differed according to niche, with less LGT observed among host-adapted DCCs versus environmental bacteria. We also found evidence suggesting that DCCs are under distinct selection pressures at both genic and intergenic sites. Our results indicate that host adaptation of M. abscessus was accompanied by major changes in genome evolution, including shifts in the apparent frequency of LGT and impacts of selection. Differences were evident among the DCCs as well, which varied in the degree of gene content remodelling, suggesting they were placed differently along the evolutionary trajectory toward host adaptation. These results provide insight into the evolutionary forces that reshape bacterial genomes as they emerge into the pathogenic niche.

The type IV secretion system of Neisseria gonorrhoeae translocates single-stranded DNA into the extracellular space, facilitating horizontal gene transfer and initiating biofilm formation. Expression of this system has been observed to be low under laboratory conditions, and multiple levels of regulation have been identified. We used a translational fusion of lacZ to traD , the gene for the type IV secretion system coupling protein, to screen for increased type IV secretion system expression. We identified several physiologically relevant conditions, including surface adherence, decreased manganese or iron, and increased zinc or copper, which increase gonococcal type IV secretion system protein levels through transcriptional and/or translational mechanisms. These metal treatments are reminiscent of the conditions in the macrophage phagosome. The ferric uptake regulator, Fur, was found to repress traD transcript levels but to also have a second role, acting to allow TraD protein levels to increase only in the absence of iron. To better understand type IV secretion system regulation during infection, we examined transcriptomic data from active urethral infection samples from five men. The data demonstrated differential expression of 20 of 21 type IV secretion system genes during infection, indicating upregulation of genes necessary for DNA secretion during host infection.

Objective: A recent study reported a tuberculosis (TB) outbreak in which, among newly infected individuals, exposure to additional active infections was associated with a higher probability of developing active disease. Referred to as complex contagion , multiple reexposures to TB within a short period after initial infection is hypothesized to confer a greater likelihood of developing active infection in 1 y. The purpose of this article is to develop and validate an agent-based simulation model (ABM) to study the effect of complex contagion on population-level TB transmission dynamics., Methods: We built an ABM of a TB epidemic using data from a series of outbreaks recorded in the 20th century in Saskatchewan, Canada. We fit 3 dynamical schemes: base, with no complex contagion; additive, in which each reexposure confers an independent risk of activated infection; and threshold, in which a small number of reexposures confers a low risk and a high number of reexposures confers a high risk of activation., Results: We find that the base model fits the mortality and incidence output targets best, followed by the threshold and then the additive models. The threshold model fits the incidence better than the base model does but overestimates mortality. All 3 models produce qualitatively realistic epidemic curves., Conclusion: We find that complex contagion qualitatively changes the trajectory of a TB epidemic, although data from a high-incidence setting are reproduced better with the base model. Results from this model demonstrate the feasibility of using ABM to capture nuances in TB transmission.

Background: The discovery of antibiotics in the mid-twentieth century marked a major transition in tuberculosis (TB) treatment and control. There are few studies describing the duration of TB disease and its treatment from the pre-chemotherapy era and little data on how these treatments changed in response to the development of effective antibiotics. The goal of this research is to understand how inpatient treatment for high incidence populations, the First Nations peoples of Saskatchewan, Canada, changed in response to increasing availability of antibiotics effective against TB. We expected that as treatment regimens transitioned from convalescence-only to triple antibiotic therapy, the length of inpatient treatment would shorten., Methods: Analyses were performed on records of sanatoria admissions and discharges occurring between 1933 and 1959 in Saskatchewan, Canada. Year of antibiotic discovery was taken as a proxy for treatment regimen: no chemotherapy (pre-1944), mono-therapy (Streptomycin, 1944-1946), dual-therapy (Streptomycin and PAS, 1946-1952), and triple-therapy (Streptomycin, PAS, and INH 1952-). A pooled linear regression of log-transformed length of first admission as predicted by year of admission was modeled to assess the relationship between admission length and year of admission, corrected for clinical and demographic variables., Results: First admission length increased 19% in the triple-therapy era as compared to the pre-chemotherapy era, from 316 days (10.4 months) to 377 days (12.4 months). After the discovery of INH (1952), we find statistically significant increases in the proportion of successfully completed therapies (0.55 versus 0.60, p = 0.035), but also in patients who left hospital against medical advice (0.19 versus 0.29, p < 0.0001), indicating that as hospitalizations lengthened, more patients chose to discharge without the sanction of their physician. The readmission rate increased from 10 to 50% of all admissions while the province-level TB-specific death rate fell from 63.1 per 10,000 in 1933 to 4.7 per 10,000 in 1958., Conclusion: Counterintuitively, we find that the length of first admissions increased with the discovery of TB-treating antibiotics. Increasing admission volume and readmission rate indicate an intensification of inpatient TB treatment during this era. These analyses provide a novel estimate of the effect of changing treatment policy on sanatorium admissions in this population.

Evidence-based public health approaches that minimize the introduction and spread of new SARS-CoV-2 transmission clusters are urgently needed in the United States and other countries struggling with expanding epidemics. Here we analyze 247 full-genome SARS-CoV-2 sequences from two nearby communities in Wisconsin, USA, and find surprisingly distinct patterns of viral spread. Dane County had the 12 th known introduction of SARS-CoV-2 in the United States, but this did not lead to descendant community spread. Instead, the Dane County outbreak was seeded by multiple later introductions, followed by limited community spread. In contrast, relatively few introductions in Milwaukee County led to extensive community spread. We present evidence for reduced viral spread in both counties following the statewide "Safer at Home" order, which went into effect 25 March 2020. Our results suggest patterns of SARS-CoV-2 transmission may vary substantially even in nearby communities. Understanding these local patterns will enable better targeting of public health interventions.

Gardnerella spp. are pathognomonic for bacterial vaginosis, which increases the risk of preterm birth and the transmission of sexually transmitted infections. Gardnerella spp. are genetically diverse, comprising what have recently been defined as distinct species with differing functional capacities. Disease associations with Gardnerella spp. are not straightforward: patients with BV are usually infected with multiple species, and Gardnerella spp. are also found in the vaginal microbiome of healthy women. Genome comparisons of Gardnerella spp. show evidence of lateral gene transfer (LGT), but patterns of LGT have not been characterized in detail. Here we sought to define the role of LGT in shaping the genetic structure of Gardnerella spp. We analyzed whole genome sequencing data for 106 Gardnerella strains and used these data for pan genome analysis and to characterize LGT in the core and accessory genomes, over recent and remote timescales. In our diverse sample of Gardnerella strains, we found that both the core and accessory genomes are clearly differentiated in accordance with newly defined species designations. We identified putative competence and pilus assembly genes across most species; we also found them to be differentiated between species. Competence machinery has diverged in parallel with the core genome, with selection against deleterious mutations as a predominant influence on their evolution. By contrast, the virulence factor vaginolysin, which encodes a toxin, appears to be readily exchanged among species. We identified five distinct prophage clusters in Gardnerella genomes, two of which appear to be exchanged between Gardnerella species. Differences among species are apparent in their patterns of LGT, including their exchange with diverse gene pools. Despite frequent LGT and co-localization in the same niche, our results show that Gardnerella spp. are clearly genetically differentiated and yet capable of exchanging specific genetic material. This likely reflects complex interactions within bacterial communities associated with the vaginal microbiome. Our results provide insight into how such interactions evolve and are maintained, allowing these multi-species communities to colonize and invade human tissues and adapt to antibiotics and other stressors., (Copyright © 2020 Bohr, Mortimer and Pepperell.)

Selective pressures imposed by pathogens have varied among human populations throughout their evolution, leading to marked inter-population differences at some genes mediating susceptibility to infectious and immune-related diseases. Here, we investigated the evolutionary history of a common polymorphism resulting in a Y529 versus C529 change in the cadherin related family member 3 (CDHR3) receptor which underlies variable susceptibility to rhinovirus-C infection and is associated with severe childhood asthma. The protective variant is the derived allele and is found at high frequency worldwide (69-95%). We detected genome-wide significant signatures of natural selection consistent with a rapid increase of the haplotypes carrying the allele, suggesting that non-neutral processes have acted on this locus across all human populations. However, the allele has not fixed in any population despite multiple lines of evidence suggesting that the mutation predates human migrations out of Africa. Using an approximate Bayesian computation method, we estimate the age of the mutation while explicitly accounting for past demography and positive or frequency-dependent balancing selection. Our analyses indicate a single emergence of the mutation in anatomically modern humans ~150 000 years ago and indicate that balancing selection has maintained the beneficial allele at high equilibrium frequencies worldwide. Apart from the well-known cases of the MHC and ABO genes, this study provides the first evidence that negative frequency-dependent selection plausibly acted on a human disease susceptibility locus, a form of balancing selection compatible with typical transmission dynamics of communicable respiratory viruses that might exploit CDHR3., (© The Author(s) 2019. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Mycobacterium tuberculosis ( Mtb ) is a globally distributed bacterial pathogen whose population structure has largely been shaped by the activities of its obligate human host. Oceania was the last major global region to be reached by Europeans and is the last region for which the dispersal and evolution of Mtb remains largely unexplored. Here, we investigated the evolutionary history of the Euro-American L4.4 sublineage and its dispersal to the South Pacific. Using a phylodynamics approach and a dataset of 236 global Mtb L4.4 genomes we have traced the origins and dispersal of L4.4 strains to New Zealand. These strains are predominantly found in indigenous Māori and Pacific people and we identify a clade of European, likely French, origin that is prevalent in indigenous populations in both New Zealand and Canada. Molecular dating suggests the expansion of European trade networks in the early 19th century drove the dispersal of this clade to the South Pacific. We also identify historical and social factors within the region that have contributed to the local spread and expansion of these strains, including recent Pacific migrations to New Zealand and the rapid urbanization of Māori in the 20th century. Our results offer new insight into the expansion and dispersal of Mtb in the South Pacific and provide a striking example of the role of historical European migrations in the global dispersal of Mtb ., (Copyright © 2019 Mulholland, Shockey, Aung, Cursons, O’Toole, Gautam, Brites, Gagneux, Roberts, Karalus, Cook, Pepperell and Arcus.)

Mycobacterium tuberculosis (M.tb) is a globally distributed, obligate pathogen of humans that can be divided into seven clearly defined lineages. An emerging consensus places the origin and global dispersal of M.tb within the past 6,000 years: identifying how the ancestral clone of M.tb spread and differentiated within this timeframe is important for identifying the ecological drivers of the current pandemic. We used Bayesian phylogeographic inference to reconstruct the migratory history of M.tb in Africa and Eurasia and to investigate lineage specific patterns of spread from a geographically diverse sample of 552 M.tb genomes. Applying evolutionary rates inferred with ancient M.tb genome calibration, we estimated the timing of major events in the migratory history of the pathogen. Inferred timings contextualize M.tb dispersal within historical phenomena that altered patterns of connectivity throughout Africa and Eurasia: trans-Indian Ocean trade in spices and other goods, the Silk Road and its predecessors, the expansion of the Roman Empire, and the European Age of Exploration. We found that Eastern Africa and Southeast Asia have been critical in the dispersal of M.tb. Our results further reveal that M.tb populations have grown through range expansion, as well as in situ, and delineate the independent evolutionary trajectories of bacterial subpopulations underlying the current pandemic., (© 2019 The Authors Molecular Ecology Published by John Wiley & Sons Ltd.)

Genetic differences are hypothesized to underlie ethnic disparities in incidence rates of the endemic systemic mycoses, including blastomycosis. Individuals of Hmong ancestry display elevated risk for this serious fungal infection. Here, we interrogated the genomes of Wisconsin (WI) Hmong blastomycosis patients using homozygosity mapping to uncover regions of the genome that are likely shared among the greater Hmong population and filtered for variants with high potential to affect disease susceptibility. This approach uncovered 113 candidate susceptibility variants, and among the most promising are those in genes involved in the interleukin-17 (IL-17) response. In particular, we identified 25 linked variants near the gene encoding IL-6 ( IL6 ). We validated differences in cytokine production between Hmong and European volunteers and formally demonstrated a critical role for IL-6 in the development of adaptive immunity to Blastomyces dermatitidis Our findings suggest that the dysregulation of IL-17 responses underlies a recently reported and poorly understood ethnic health disparity. IMPORTANCE Blastomycosis is a potentially life-threatening infection caused by the fungus Blastomyces dermatitidis As with related fungal diseases, blastomycosis is noted to affect some populations more than others. These patterns of illness are often not related to predisposing conditions or exposure risks; thus, genetic differences are thought to underlie these health disparities. People of Hmong ancestry in Wisconsin are at elevated risk of blastomycosis compared to the general population. We studied the genetic codes of Hmong blastomycosis patients and identified candidate sites in their genomes that may explain their susceptibility to this infection. We further studied one particular region of the genome that is involved with the immune processes that fight B. dermatitidis Our work revealed population differences in the response to fungi. A better understanding of the genetic underpinnings of susceptibility to infectious diseases has broader implications for community health, especially in the paradigm of personalized medicine., (Copyright © 2019 Merkhofer et al.)

Mycobacterium tuberculosis ( M. tb ), an obligate human pathogen and the etiological agent of tuberculosis (TB), remains a major threat to global public health. Comparative genomics has been invaluable for monitoring the emergence and spread of TB and for gaining insight into adaptation of M. tb . Most genomic studies of M. tb are based on single bacterial isolates that have been cultured for several weeks in vitro . However, in its natural human host, M. tb comprises complex, in some cases massive bacterial populations that diversify over the course of infection and cannot be wholly represented by a single genome. Recently, enrichment via hybridization capture has been used as a rapid diagnostic tool for TB, circumventing culturing protocols and enabling the recovery of M. tb genomes directly from sputum. This method has further applicability to the study of M. tb adaptation, as it enables a higher resolution and more direct analysis of M. tb genetic diversity within hosts with TB. Here we analyzed genomic material from M. tb and Mycobacterium bovis populations captured directly from sputum and from cultured samples using metagenomic and Pool-Seq approaches. We identified effects of sampling, patient, and sample type on bacterial genetic diversity. Bacterial genetic diversity was more variable and on average higher in sputum than in culture samples, suggesting that manipulation in the laboratory reshapes the bacterial population. Using outlier analyses, we identified candidate bacterial genetic loci mediating adaptation to these distinct environments. The study of M. tb in its natural human host is a powerful tool for illuminating host pathogen interactions and understanding the bacterial genetic underpinnings of virulence.

A recent study reported on a tuberculosis (TB) outbreak in a largely Inuit village. Among newly infected individuals, exposure to additional active cases was associated with an increasing probability of developing active disease within a year. Using binomial risk models, we evaluated two potential mechanisms by which multiple infections during the first year following initial infection could account for increasing disease risk with increasing exposures. In the reinfection model , each infectious contact confers an independent risk of an infection, and infections contribute independently to active disease. In the threshold model , disease risk follows a sigmoidal function with small numbers of infectious contacts conferring a low risk of active disease and large numbers of contacts conferring a high risk. To determine the dynamic impact of reinfection during the early phase of infection, we performed simulations from a modified Reed-Frost model of TB dynamics following spread from an initial number of cases. We parametrized this model with the maximum-likelihood estimates from the reinfection and threshold models in addition to the observed distribution of exposures among new infections. We find that both models can plausibly account for the observed increase in disease risk with increasing infectious contacts, but the threshold model confers a better fit than a nested model without a threshold ( p = 0.04). Our simulations indicate that multiple exposures to infectious individuals during this critical time period can lead to dramatic increases in outbreak size. In order to decrease TB burden in high-prevalence settings, it may be necessary to implement measures aimed at preventing repeated exposures, in addition to preventing primary infection., Competing Interests: The authors declare no competing interests.

A small number of high-burden countries account for the majority of tuberculosis cases worldwide. Detailed data are lacking from these regions. To explore the evolutionary history of Mycobacterium tuberculosis in China-the country with the third highest tuberculosis burden-we analysed a countrywide collection of 4,578 isolates. Little genetic diversity was detected, with 99.4% of the bacterial population belonging to lineage 2 and three sublineages of lineage 4. The deeply rooted phylogenetic positions and geographic restriction of these four genotypes indicate that their populations expanded in situ following a small number of introductions to China. Coalescent analyses suggest that these bacterial subpopulations emerged in China around 1,000 years ago, and expanded in parallel from the twelfth century onwards, and that the whole population peaked in the late eighteenth century. More recently, sublineage L2.3, which is indigenous to China and exhibited relatively high transmissibility and extensive global dissemination, came to dominate the population dynamics of M. tuberculosis in China. Our results indicate that historical expansion of four M. tuberculosis strains shaped the current tuberculosis epidemic in China, and highlight the long-term genetic continuity of the indigenous M. tuberculosis population.

On the basis of population genomic and phylogeographic analyses of 1669 Mycobacterium tuberculosis lineage 4 (L4) genomes, we find that dispersal of L4 has been completely dominated by historical migrations out of Europe. We demonstrate an intimate temporal relationship between European colonial expansion into Africa and the Americas and the spread of L4 tuberculosis (TB). Markedly, in the age of antibiotics, mutations conferring antimicrobial resistance overwhelmingly emerged locally (at the level of nations), with minimal cross-border transmission of resistance. The latter finding was found to reflect the relatively recent emergence of these mutations, as a similar degree of local restriction was observed for susceptible variants emerging on comparable time scales. The restricted international transmission of drug-resistant TB suggests that containment efforts at the level of individual countries could be successful.

Background: In Canada, tuberculosis disproportionately affects foreign-born and First Nations populations. Within First Nations' peoples, a high proportion of cases occur in association with outbreaks. Tuberculosis transmission in the context of outbreaks is thought to result from the convergence of several factors including characteristics of the cases, contacts, the environment, and the pathogen., Methods: We examined the epidemiological and genomic determinants of two well-characterized tuberculosis outbreaks attributed to two super-spreaders among First Nations in the province of Alberta. These outbreaks were associated with two distinct DNA fingerprints (restriction fragment-length polymorphisms or RFLPs 0.0142 and 0.0728). We compared outbreak isolates with endemic isolates not spatio-temporarily linked to outbreak cases. We extracted epidemiological variables pertaining to tuberculosis cases and contacts from individual public health records and the provincial tuberculosis registry. We conducted group analyses using parametric and non-parametric statistical tests. We carried out whole-genome sequencing and bioinformatic analysis using validated protocols., Results: We observed differences between outbreak and endemic groups in the mean number of total and child-aged contacts and the number of contacts with new positive and converted tuberculin skin tests in all group comparisons (p < 0.05). Differences were also detected in the proportion of cases with cavitation on a chest radiograph and the mean number of close contacts in selected group comparisons (p < 0.02). A phylogenetic network analysis of whole-genome sequencing data indicated that most outbreak and endemic strains were closely related to the source case for the 0.0142 fingerprint. For the 0.0728 fingerprint, the source case haplotype was circulating among endemic cases prior to the outbreak. Genetic and temporal distances were not correlated for either RFLP 0.0142 (r 2  = - 0.05) or RFLP 0.0728 (r 2  = 0.09) when all isolates were analyzed., Conclusions: We found no evidence that endemic strains acquired mutations resulting in their emergence in outbreak form. We conclude that the propagation of these outbreaks was likely driven by the combination of characteristics of the source cases, contacts, and the environment. The role of whole-genome sequencing in understanding mycobacterial evolution and in assisting public health authorities in conducting contact investigations and managing outbreaks is important and expected to grow in the future.

Tuberculosis (TB) is the leading cause of death by an infectious disease, and global TB control efforts are increasingly threatened by drug resistance in Mycobacterium tuberculosis . Unlike most bacteria, where lateral gene transfer is an important mechanism of resistance acquisition, resistant M. tuberculosis arises solely by de novo chromosomal mutation. Using whole-genome sequencing data from two natural populations of M. tuberculosis , we characterized the population genetics of known drug resistance loci using measures of diversity, population differentiation, and convergent evolution. We found resistant subpopulations to be less diverse than susceptible subpopulations, consistent with ongoing transmission of resistant M. tuberculosis . A subset of resistance genes ("sloppy targets") were characterized by high diversity and multiple rare variants; we posit that a large genetic target for resistance and relaxation of purifying selection contribute to high diversity at these loci. For "tight targets" of selection, the path to resistance appeared narrower, evidenced by single favored mutations that arose numerous times in the phylogeny and segregated at markedly different frequencies in resistant and susceptible subpopulations. These results suggest that diverse genetic architectures underlie drug resistance in M. tuberculosis and that combined approaches are needed to identify causal mutations. Extrapolating from patterns observed for well-characterized genes, we identified novel candidate variants involved in resistance. The approach outlined here can be extended to identify resistance variants for new drugs, to investigate the genetic architecture of resistance, and when phenotypic data are available, to find candidate genetic loci underlying other positively selected traits in clonal bacteria. IMPORTANCE Mycobacterium tuberculosis , the causative agent of tuberculosis (TB), is a significant burden on global health. Antibiotic treatment imposes strong selective pressure on M. tuberculosis populations. Identifying the mutations that cause drug resistance in M. tuberculosis is important for guiding TB treatment and halting the spread of drug resistance. Whole-genome sequencing (WGS) of M. tuberculosis isolates can be used to identify novel mutations mediating drug resistance and to predict resistance patterns faster than traditional methods of drug susceptibility testing. We have used WGS from natural populations of drug-resistant M. tuberculosis to characterize effects of selection for advantageous mutations on patterns of diversity at genes involved in drug resistance. The methods developed here can be used to identify novel advantageous mutations, including new resistance loci, in M. tuberculosis and other clonal pathogens.

Human-pathogenic bacteria are found in a variety of niches, including free-living, zoonotic, and microbiome environments. Identifying bacterial adaptations that enable invasive disease is an important means of gaining insight into the molecular basis of pathogenesis and understanding pathogen emergence. Staphylococcus saprophyticus , a leading cause of urinary tract infections, can be found in the environment, food, animals, and the human microbiome. We identified a selective sweep in the gene encoding the Aas adhesin, a key virulence factor that binds host fibronectin. We hypothesize that the mutation under selection ( aas &#95;2206A>C) facilitates colonization of the urinary tract, an environment where bacteria are subject to strong shearing forces. The mutation appears to have enabled emergence and expansion of a human-pathogenic lineage of S. saprophyticus . These results demonstrate the power of evolutionary genomic approaches in discovering the genetic basis of virulence and emphasize the pleiotropy and adaptability of bacteria occupying diverse niches. IMPORTANCE Staphylococcus saprophyticus is an important cause of urinary tract infections (UTI) in women; such UTI are common, can be severe, and are associated with significant impacts to public health. In addition to being a cause of human UTI, S. saprophyticus can be found in the environment, in food, and associated with animals. After discovering that UTI strains of S. saprophyticus are for the most part closely related to each other, we sought to determine whether these strains are specially adapted to cause disease in humans. We found evidence suggesting that a mutation in the gene aas is advantageous in the context of human infection. We hypothesize that the mutation allows S. saprophyticus to survive better in the human urinary tract. These results show how bacteria found in the environment can evolve to cause disease.

Mycobacteria have a distinct secretion system, termed type VII (T7SS), which is encoded by paralogous chromosomal loci (ESX) and associated with pathogenesis, conjugation, and metal homeostasis. Evolution of paralogous gene families is of interest because duplication is an important mechanism by which novel genes evolve, but there are potential conflicts between adaptive forces that stabilize duplications and those that enable evolution of new functions. Our objective was to delineate the adaptive forces underlying diversification of T7SS. Plasmid-borne ESX were described recently, and we found evidence that the initial duplication and divergence of ESX systems occurred on plasmids and was driven by selection for advantageous mutations. Plasmid conjugation has been linked to T7SS and type IV secretion systems (T4SS) in mycobacteria, and we discovered that T7SS and T4SS genes evolved in concert on the plasmids. We hypothesize that differentiation of plasmid ESX helps to prevent conjugation among cells harboring incompatible plasmids. Plasmid ESX appear to have been repurposed following migration to the chromosome, and there is evidence of positive selection driving further differentiation of chromosomal ESX. We hypothesize that ESX loci were initially stabilized on the chromosome by mediating their own transfer. These results emphasize the diverse adaptive paths underlying evolution of novelty, which in this case involved plasmid duplications, selection for advantageous mutations in the mobile and core genomes, migration of the loci between plasmids and chromosomes, and lateral transfer among chromosomes. We discuss further implications for the choice of model organism to study ESX functions in Mycobacterium tuberculosis., (© The Author(s) 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

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., Competing Interests: The authors declare that no competing interests exist.

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., Competing Interests: The authors declare no conflict of interest.

Lateral gene transfer is an important mechanism for evolution among bacteria. Here, genome-wide gene insertion and deletion rates are modeled in a maximum-likelihood framework with the additional flexibility of modeling potential missing data. The performance of the models is illustrated using simulations and a data set on gene family phyletic patterns from Gardnerella vaginalis that includes an ancient taxon. A novel application involving pseudogenization/genome reduction magnitudes is also illustrated, using gene family data from Mycobacterium spp. Finally, an R package called indelmiss is available from the Comprehensive R Archive Network at https://cran.r-project.org/package=indelmiss, with support documentation and examples., (Copyright © 2016 by the Genetics Society of America.)

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.

Late 19th century epidemics of tuberculosis (TB) in Western Canadian First Nations resulted in peak TB mortality rates more than six times the highest rates recorded in Europe. Using a mathematical modeling approach and historical TB mortality time series, we investigate potential causes of high TB mortality and rapid epidemic decline in First Nations from 1885 to 1940. We explore two potential causes of dramatic epidemic dynamics observed in this setting: first, we explore effects of famine prior to 1900 on both TB and population dynamics. Malnutrition is recognized as an individual-level risk factor for TB progression and mortality; its population-level effects on TB epidemics have not been explored previously. Second, we explore effects of heterogeneity in susceptibility to TB in two ways: modeling heterogeneity in susceptibility to infection, and heterogeneity in risk of developing disease once infected. Our results indicate that models lacking famine-related changes in TB parameters or heterogeneity result in an implausibly poor fit to both the TB mortality time series and census data; the inclusion of these features allows for the characteristic decline and rise in population observed in First Nations during this time period and confers improved fits to TB mortality data.

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., (© The Author(s) 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

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.

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 DS6(Quebec) 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.