Population genomics of Chlamydia trachomatis: insights on drift, selection, recombination, and population structure

Rosmarinic acid activates the Ras/Raf/MEK/ERK signaling pathway to regulate CD8+ T cells and autophagy to clear Chlamydia trachomatis in reproductive tract-infected mice

Chlamydia trachomatis (CT) is the leading cause of bacterial sexually transmitted diseases worldwide, which can cause diseases such as pelvic inflammatory disease, and cervical and fallopian tube inflammation, and poses a threat to human health. Rosmarinic acid (RosA) is an active ingredient of natural products with anti-inflammatory and immunomodulatory effects. This study aimed to investigate the role of RosA in inhibiting autophagy-regulated immune cells-CD8+ T cells via the Ras/Raf/MEK/ERK signaling pathway in a CT-infected mouse model. Mice were inoculated with CT infection solution vaginally, and the mechanistic basis of RosA treatment was established using H&E staining, flow cytometry, immunofluorescence, transmission electron microscopy, and western blot. The key factors involved in RosA treatment were further validated using the MEK inhibitor cobimetinib. Experimental results showed that both RosA and the reference drug azithromycin could attenuate the pathological damage to the endometrium caused by CT infection; flow cytometry showed that peripheral blood CD8+ T cells increased after CT infection and decreased after treatment with RosA and the positive drug azithromycin (positive control); immunofluorescence showed that endometrial CD8 and LC3 increased after CT infection and decreased after RosA and positive drug treatment; the results of transmission electron microscopy showed that RosA and the positive drug azithromycin inhibited the accumulation of autophagosomes; western bolt experiments confirmed the activation of autophagy proteins LC3Ⅱ/Ⅰ, ATG5, Beclin-1, and p62 after CT infection, as well as the inhibition of Ras/Raf/MEK/ERK signaling. RosA and azithromycin inhibition of autophagy proteins activates Ras/Raf/MEK/ERK signaling. In addition, the MEK inhibitor cobimetinib attenuated RosA's protective effect on endometrium by further activating CD8+ T cells on a CT-induced basis, while transmission electron microscopy, immunofluorescence, and western blots showed that cobimetinib blocked ERK signals activation and further induced phagocytosis on a CT-induced basis. These data indicated that RosA can activate the Ras/Raf/MEK/ERK signaling pathway to inhibit autophagy, and RosA could also regulate the activation of immune cells-CD8+T cells to protect the reproductive tract of CT-infected mice.

Genomic insights into local-scale evolution of ocular Chlamydia trachomatis strains within and between individuals in Gambian trachoma-endemic villages

Trachoma, a neglected tropical disease caused by Chlamydia trachomatis (Ct) serovars A-C, is the leading infectious cause of blindness worldwide. Africa bears the highest burden, accounting for over 86 % of global trachoma cases. We investigated Ct serovar A (SvA) and B (SvB) whole genome sequences prior to the induction of mass antibiotic drug administration in The Gambia. Here, we explore the factors contributing to Ct strain diversification and the implications for Ct evolution within the context of ocular infection. A cohort study in 2002-2003 collected ocular swabs across nine Gambian villages during a 6 month follow-up study. To explore the genetic diversity of Ct within and between individuals, we conducted whole-genome sequencing (WGS) on a limited number (n=43) of Ct-positive samples with an omcB load ≥10 from four villages. WGS was performed using target enrichment with SureSelect and Illumina paired-end sequencing. Out of 43 WGS samples, 41 provided sufficient quality for further analysis. ompA analysis revealed that 11 samples had highest identity to ompA from strain A/HAR13 (NC_007429) and 30 had highest identity to ompA from strain B/Jali20 (NC_012686). While SvB genome sequences formed two distinct village-driven subclades, the heterogeneity of SvA sequences led to the formation of many individual branches within the Gambian SvA subclade. Comparing the Gambian SvA and SvB sequences with their reference strains, Ct A/HAR13 and Ct B/Jali20, indicated an single nucleotide polymorphism accumulation rate of 2.4×10-5 per site per year for the Gambian SvA and 1.3×10-5 per site per year for SvB variants (P<0.0001). Variant calling resulted in a total of 1371 single nucleotide variants (SNVs) with a frequency >25 % in SvA sequences, and 438 SNVs in SvB sequences. Of note, in SvA variants, highest evolutionary pressure was recorded on genes responsible for host cell modulation and intracellular survival mechanisms, whereas in SvB variants this pressure was mainly on genes essential for DNA replication/repair mechanisms and protein synthesis. A comparison of the sequences between observed separate infection events (4-20 weeks between infections) suggested that the majority of the variations accumulated in genes responsible for host-pathogen interaction such as CTA_0166 (phospholipase D-like protein), CTA_0498 (TarP) and CTA_0948 (deubiquitinase). This comparison of Ct SvA and SvB variants within a trachoma endemic population focused on their local evolutionary adaptation. We found a different variation accumulation pattern in the Gambian SvA chromosomal genes compared with SvB, hinting at the potential of Ct serovar-specific variation in diversification and evolutionary fitness. These findings may have implications for optimizing trachoma control and prevention strategies.

Rapid expansion of lymphogranuloma venereum infections with fast diversification and spread of Chlamydia trachomatis L genovariants

IMPORTANCE

Numerous international organizations, including the World Health Organization, have been drawing attention to the global increase in sexually transmitted infections. Twenty years ago, lymphogranuloma venereum (LGV) was mainly considered a tropical disease; in recent decades, however, LGV has been increasingly present in high-income countries. This increase has been linked to men who have sex with men who participate in highly interconnected sexual networks, leading to a rapid spread of LGV. This study focuses on the spread of LGV, presenting the largest time series of LGV prevalence in Spain, which includes more than a thousand diagnosed cases in one large city. The number of LGV cases diagnosed was analyzed over time, and a selection of strains was subjected to molecular genotyping. The results indicate that the LGV epidemic is gradually evolving toward an increasingly complex diversification due to the selection of successful genovariants that have emerged by mutation and recombination events, suggesting that we are moving toward an unpredictable scenario.

Recombination in Bacterial Genomes: Evolutionary Trends

Bacterial organisms have undergone homologous recombination (HR) and horizontal gene transfer (HGT) multiple times during their history. These processes could increase fitness to new environments, cause specialization, the emergence of new species, and changes in virulence. Therefore, comprehensive knowledge of the impact and intensity of genetic exchanges and the location of recombination hotspots on the genome is necessary for understanding the dynamics of adaptation to various conditions. To this end, we aimed to characterize the functional impact and genomic context of computationally detected recombination events by analyzing genomic studies of any bacterial species, for which events have been detected in the last 30 years. Genomic loci where the transfer of DNA was detected pertained to mobile genetic elements (MGEs) housing genes that code for proteins engaged in distinct cellular processes, such as secretion systems, toxins, infection effectors, biosynthesis enzymes, etc. We found that all inferences fall into three main lifestyle categories, namely, ecological diversification, pathogenesis, and symbiosis. The latter primarily exhibits ancestral events, thus, possibly indicating that adaptation appears to be governed by similar recombination-dependent mechanisms.

Recent advances in genetic systems in obligate intracellular human-pathogenic bacteria

The ability to genetically manipulate a pathogen is fundamental to discovering factors governing host-pathogen interactions at the molecular level and is critical for devising treatment and prevention strategies. While the genetic "toolbox" for many important bacterial pathogens is extensive, approaches for modifying obligate intracellular bacterial pathogens were classically limited due in part to the uniqueness of their obligatory lifestyles. Many researchers have confronted these challenges over the past two and a half decades leading to the development of multiple approaches to construct plasmid-bearing recombinant strains and chromosomal gene inactivation and deletion mutants, along with gene-silencing methods enabling the study of essential genes. This review will highlight seminal genetic achievements and recent developments (past 5 years) for Anaplasma spp., Rickettsia spp., Chlamydia spp., and Coxiella burnetii including progress being made for the still intractable Orientia tsutsugamushi. Alongside commentary of the strengths and weaknesses of the various approaches, future research directions will be discussed to include methods for C. burnetii that should have utility in the other obligate intracellular bacteria. Collectively, the future appears bright for unraveling the molecular pathogenic mechanisms of these significant pathogens.

Genome organization and genomics in Chlamydia: whole genome sequencing increases understanding of chlamydial virulence, evolution, and phylogeny

The genus Chlamydia contains important obligate intracellular bacterial pathogens to humans and animals, including C. trachomatis and C. pneumoniae. Since 1998, when the first Chlamydia genome was published, our understanding of how these microbes interact, evolved and adapted to different intracellular host environments has been transformed due to the expansion of chlamydial genomes. This review explores the current state of knowledge in Chlamydia genomics and how whole genome sequencing has revolutionised our understanding of Chlamydia virulence, evolution, and phylogeny over the past two and a half decades. This review will also highlight developments in multi-omics and other approaches that have complemented whole genome sequencing to advance knowledge of Chlamydia pathogenesis and future directions for chlamydial genomics.

Patterns of within-host spread of Chlamydia trachomatis between vagina, endocervix and rectum revealed by comparative genomic analysis

Introduction

Chlamydia trachomatis, a gram-negative obligate intracellular bacterium, commonly causes sexually transmitted infections (STIs). Little is known about C. trachomatis transmission within the host, which is important for understanding disease epidemiology and progression.

Methods

We used RNA-bait enrichment and whole-genome sequencing to compare rectal, vaginal and endocervical samples collected at the same time from 26 study participants who attended Fijian Ministry of Health and Medical Services clinics and tested positive for C. trachomatis at each anatomic site.

Results

The 78 C. trachomatis genomes from participants resolved into two major clades of the C. trachomatis phylogeny (the "prevalent urogenital and anorectal" clade and "non-prevalent urogenital and anorectal" clade). For 21 participants, genome sequences were almost identical in each anatomic site. For the other five participants, two distinct C. trachomatis strains were present in different sites; in two cases, the vaginal sample was a mixture of strains.

Discussion

The absence of large numbers of fixed SNPs between C. trachomatis genomes within many of the participants could indicate recent acquisition of infection prior to the clinic visit without sufficient time to accumulate significant genetic variation in different body sites. This model suggests that many C. trachomatis infections may be resolved relatively quickly in the Fijian population, possibly reflecting common prescription or over-the-counter antibiotics usage.

One to host them all: genomics of the diverse bacterial endosymbionts of the spider Oedothorax gibbosus

Bacterial endosymbionts of the groups Wolbachia, Cardinium and Rickettsiaceae are well known for their diverse effects on their arthropod hosts, ranging from mutualistic relationships to reproductive phenotypes. Here, we analysed a unique system in which the dwarf spider Oedothorax gibbosus is co-infected with up to five different endosymbionts affiliated with Wolbachia, 'Candidatus Tisiphia' (formerly Torix group Rickettsia), Cardinium and Rhabdochlamydia. Using short-read genome sequencing data, we show that the endosymbionts are heterogeneously distributed among O. gibbosus populations and are frequently found co-infecting spider individuals. To study this intricate host-endosymbiont system on a genome-resolved level, we used long-read sequencing to reconstruct closed genomes of the Wolbachia, 'Ca. Tisiphia' and Cardinium endosymbionts. We provide insights into the ecology and evolution of the endosymbionts and shed light on the interactions with their spider host. We detected high quantities of transposable elements in all endosymbiont genomes and provide evidence that ancestors of the Cardinium, 'Ca. Tisiphia' and Wolbachia endosymbionts have co-infected the same hosts in the past. Our findings contribute to broadening our knowledge about endosymbionts infecting one of the largest animal phyla on Earth and show the usefulness of transposable elements as an evolutionary 'contact-tracing' tool.

Patterns of within-host spread of Chlamydia trachomatis between vagina, endocervix and rectum revealed by comparative genomic analysis

Chlamydia trachomatis , a gram-negative obligate intracellular bacterium, commonly causes sexually transmitted infections (STIs). Little is known about C. trachomatis transmission within the host, which is important for understanding disease epidemiology and progression. We used RNA-bait enrichment and whole-genome sequencing to compare rectal, vaginal and endocervical samples collected at the same time from 26 study participants who attended Fijian Ministry of Health and Medical Services clinics and tested positive for C. trachomatis at each anatomic site. The 78 C. trachomatis genomes from participants were from two major clades of the C. trachomatis phylogeny (the "prevalent urogenital and anorecta"l clade and "non-prevalent urogenital and anorectal" clade). For 21 participants, genome sequences were almost identical in each anatomic site. For the other five participants, two distinct C. trachomatis strains were present in different sites; in two cases, the vaginal sample was a mixture of strains. The absence of large numbers of fixed SNPs between C. trachomatis strains within many of the participants could indicate recent acquisition of infection prior to the clinic visit without sufficient time to accumulate significant variation in the different body sites. This model suggests that many C. trachomatis infections may be resolved relatively quickly in the Fijian population, possibly reflecting common prescription or over-the-counter antibiotics usage.

Importance

Chlamydia trachomatis is a bacterial pathogen that causes millions of sexually transmitted infections (STIs) annually across the globe. Because C. trachomatis lives inside human cells, it has historically been hard to study. We know little about how the bacterium spreads between body sites. Here, samples from 26 study participants who had simultaneous infections in their vagina, rectum and endocervix were genetically analyzed using an improved method to extract C. trachomatis DNA directly from clinical samples for genome sequencing. By analyzing patterns of mutations in the genomes, we found that 21 participants shared very similar C. trachomatis strains in all three anatomic sites, suggesting recent infection and spread. For five participants two C. trachomatis strains were evident, indicating multiple infections. This study is significant in that improved enrichment methods for genome sequencing provides robust data to genetically trace patterns of C. trachomatis infection and transmission within an individual for epidemiologic and pathogenesis interrogations.

The Impact of Lateral Gene Transfer in Chlamydia

Lateral gene transfer (LGT) facilitates many processes in bacterial ecology and pathogenesis, especially regarding pathogen evolution and the spread of antibiotic resistance across species. The obligate intracellular chlamydiae, which cause a range of diseases in humans and animals, were historically thought to be highly deficient in this process. However, research over the past few decades has demonstrated that this was not the case. The first reports of homologous recombination in the Chlamydiaceae family were published in the early 1990s. Later, the advent of whole-genome sequencing uncovered clear evidence for LGT in the evolution of the Chlamydiaceae, although the acquisition of tetracycline resistance in Chlamydia (C.) suis is the only recent instance of interphylum LGT. In contrast, genome and in vitro studies have shown that intraspecies DNA exchange occurs frequently and can even cross species barriers between closely related chlamydiae, such as between C. trachomatis, C. muridarum, and C. suis. Additionally, whole-genome analysis led to the identification of various DNA repair and recombination systems in C. trachomatis, but the exact machinery of DNA uptake and homologous recombination in the chlamydiae has yet to be fully elucidated. Here, we reviewed the current state of knowledge concerning LGT in Chlamydia by focusing on the effect of homologous recombination on the chlamydial genome, the recombination machinery, and its potential as a genetic tool for Chlamydia.

Chlamydia trachomatis Cross-Serovar Protection during Experimental Lung Reinfection in Mice

Chlamydia trachomatis causes most bacterial sexually transmitted diseases worldwide. Different major outer membrane proteins (MOMPs) define various serovars of this intracellular pathogen: In women, D to L3 can cause urethritis, cervicitis, salpingitis, and oophoritis, and, thus, infertility. Protective immunity might be serovar-specific since chlamydial infection does not appear to induce an effective acquired immunity and reinfections occur. A better understanding of induced cross-serovar protection is essential for the selection of suitable antigens in vaccine development. In our mouse lung infection screening model, we evaluated the urogenital serovars D, E, and L2 in this regard. Seven weeks after primary infection or mock-infection, respectively, mice were infected a second time with the identical or one of the other serovars. Body weight and clinical score were monitored for 7 days. Near the peak of the second lung infection, bacterial load, myeloperoxidase, IFN-γ, and TNF-α in lung homogenate, as well as chlamydia-specific IgG levels in blood were determined. Surprisingly, compared with mice that were infected then for the first time, almost independent of the serovar combination used, all acquired parameters of disease were similarly diminished. Our reinfection study suggests that efficient cross-serovar protection could be achieved by a vaccine combining chlamydial antigens that do not include nonconserved MOMP regions.

Ongoing evolution of Chlamydia trachomatis lymphogranuloma venereum: exploring the genomic diversity of circulating strains

Lymphogranuloma venereum (LGV), the invasive infection of the sexually transmissible infection (STI) Chlamydia trachomatis, is caused by strains from the LGV biovar, most commonly represented by ompA-genotypes L2b and L2. We investigated the diversity in LGV samples across an international collection over seven years using typing and genome sequencing. LGV-positive samples (n=321) from eight countries collected between 2011 and 2017 (Spain n=97, Netherlands n=67, Switzerland n=64, Australia n=53, Sweden n=37, Hungary n=31, Czechia n=30, Slovenia n=10) were genotyped for pmpH and ompA variants. All were found to contain the 9 bp insertion in the pmpH gene, previously associated with ompA-genotype L2b. However, analysis of the ompA gene shows ompA-genotype L2b (n=83), ompA-genotype L2 (n=180) and several variants of these (n=52; 12 variant types), as well as other/mixed ompA-genotypes (n=6). To elucidate the genomic diversity, whole genome sequencing (WGS) was performed from selected samples using SureSelect target enrichment, resulting in 42 genomes, covering a diversity of ompA-genotypes and representing most of the countries sampled. A phylogeny of these data clearly shows that these ompA-genotypes derive from an ompA-genotype L2b ancestor, carrying up to eight SNPs per isolate. SNPs within ompA are overrepresented among genomic changes in these samples, each of which results in an amino acid change in the variable domains of OmpA (major outer membrane protein, MOMP). A reversion to ompA-genotype L2 with the L2b genomic backbone is commonly seen. The wide diversity of ompA-genotypes found in these recent LGV samples indicates that this gene is under immunological selection. Our results suggest that the ompA-genotype L2b genomic backbone is the dominant strain circulating and evolving particularly in men who have sex with men (MSM) populations.

Prophylactic Multi-Subunit Vaccine against Chlamydia trachomatis: In Vivo Evaluation in Mice

Chlamydia trachomatis is the most frequent sexually-transmitted disease-causing bacterium. Urogenital serovars of this intracellular pathogen lead to urethritis and cervicitis. Ascending infections result in pelvic inflammatory disease, salpingitis, and oophoritis. One of 200 urogenital infections leads to tubal infertility. Serovars A–C cause trachoma with visual impairment. There is an urgent need for a vaccine. We characterized a new five-component subunit vaccine in a mouse vaccination-lung challenge infection model. Four recombinant Pmp family-members and Ctad1 from C. trachomatis serovar E, all of which participate in adhesion and binding of chlamydial elementary bodies to host cells, were combined with the mucosal adjuvant cyclic-di-adenosine monophosphate. Intranasal application led to a high degree of cross-serovar protection against urogenital and ocular strains of C. trachomatis, which lasted at least five months. Critical evaluated parameters were body weight, clinical score, chlamydial load, a granulocyte marker and the cytokines IFN-γ/TNF-α in lung homogenate. Vaccine antigen-specific antibodies and a mixed Th1/Th2/Th17 T cell response with multi-functional CD4⁺ and CD8⁺ T cells correlate with protection. However, serum-transfer did not protect the recipients suggesting that circulating antibodies play only a minor role. In the long run, our new vaccine might help to prevent the feared consequences of human C. trachomatis infections.

Tryptophan Operon Diversity Reveals Evolutionary Trends among Geographically Disparate Chlamydia trachomatis Ocular and Urogenital Strains Affecting Tryptophan Repressor and Synthase Function

The obligate intracellular pathogen Chlamydia trachomatis (Ct) is the leading cause of bacterial sexually transmitted infections and blindness globally. To date, Ct urogenital strains are considered tryptophan prototrophs, utilizing indole for tryptophan synthesis within a closed-conformation tetramer comprised of two α (TrpA)- and two β (TrpB)-subunits. In contrast, ocular strains are auxotrophs due to mutations in TrpA, relying on host tryptophan pools for survival. It has been speculated that there is strong selective pressure for urogenital strains to maintain a functional operon. Here, we performed genetic, phylogenetic, and novel functional modeling analyses of 595 geographically diverse Ct ocular, urethral, vaginal, and rectal strains with complete operon sequences. We found that ocular and urogenital, but not lymphogranuloma venereum, TrpA-coding sequences were under positive selection. However, vaginal and urethral strains exhibited greater nucleotide diversity and a higher ratio of nonsynonymous to synonymous substitutions [Pi(a)/Pi(s)] than ocular strains, suggesting a more rapid evolution of beneficial mutations. We also identified nonsynonymous amino acid changes for an ocular isolate with a urogenital backbone in the intergenic region between TrpR and TrpB at the exact binding site for YtgR-the only known iron-dependent transcription factor in Chlamydia-indicating that selective pressure has disabled the response to fluctuating iron levels. In silico effects on protein stability, ligand-binding affinity, and tryptophan repressor (TrpR) affinity for single-stranded DNA (ssDNA) measured by calculating free energy changes (ΔΔG) between Ct reference and mutant tryptophan operon proteins were also analyzed. We found that tryptophan synthase function was likely suboptimal compared to other bacterial tryptophan prototrophs and that a diversity of urogenital strain mutations rendered the synthase nonfunctional or inefficient. The novel mutations identified here affected active sites in an orthosteric manner but also hindered α- and β-subunit allosteric interactions from distant sites, reducing efficiency of the tryptophan synthase. Importantly, strains with mutant proteins were inclined toward energy conservation by exhibiting an altered affinity for their respective ligands compared to reference strains, indicating greater fitness. This is not surprising as l-tryptophan is one of the most energetically costly amino acids to synthesize. Mutations in the tryptophan repressor gene (trpR) among urogenital strains were similarly detrimental to function. Our findings indicate that urogenital strains are evolving more rapidly than previously recognized with mutations that impact tryptophan operon function in a manner that is energetically beneficial, providing a novel host-pathogen evolutionary mechanism for intracellular survival.IMPORTANCE Chlamydia trachomatis (Ct) is a major global public health concern causing sexually transmitted and ocular infections affecting over 130 million and 260 million people, respectively. Sequelae include infertility, preterm birth, ectopic pregnancy, and blindness. Ct relies on available host tryptophan pools and/or substrates to synthesize tryptophan to survive. Urogenital strains synthesize tryptophan from indole using their intact tryptophan synthase (TS). Ocular strains contain a trpA frameshift mutation that encodes a truncated TrpA with loss of TS function. We found that TS function is likely suboptimal compared to other tryptophan prototrophs and that urogenital stains contain diverse mutations that render TS nonfunctional/inefficient, evolve more rapidly than previously recognized, and impact operon function in a manner that is energetically beneficial, providing an alternative host-pathogen evolutionary mechanism for intracellular survival. Our research has broad scientific appeal since our approach can be applied to other bacteria that may explain evolution/survival in host-pathogen interactions.

Homologous Recombination in Clostridioides difficile Mediates Diversification of Cell Surface Features and Transport Systems

Infections with C. difficile result in up to half a million illnesses and tens of thousands of deaths annually in the United States. The severity of C. difficile illness is dependent on both host and bacterial factors.

Illness caused by the pathogen Clostridioides difficile is widespread and can range in severity from mild diarrhea to sepsis and death. Strains of C. difficile isolated from human infections exhibit great genetic diversity, leading to the hypothesis that the genetic background of the infecting strain at least partially determines a patient’s clinical course. However, although certain strains of C. difficile have been suggested to be associated with increased severity, strain typing alone has proved insufficient to explain infection severity. The limited explanatory power of strain typing has been hypothesized to be due to genetic variation within strain types, as well as genetic elements shared between strain types. Homologous recombination is an evolutionary mechanism that can result in large genetic differences between two otherwise clonal isolates, and also lead to convergent genotypes in distantly related strains. More than 400 C. difficile genomes were analyzed here to assess the effect of homologous recombination within and between C. difficile clades. Almost three-quarters of single nucleotide variants in the C. difficile phylogeny are predicted to be due to homologous recombination events. Furthermore, recombination events were enriched in genes previously reported to be important to virulence and host-pathogen interactions, such as flagella, cell wall proteins, and sugar transport and metabolism. Thus, by exploring the landscape of homologous recombination in C. difficile, we identified genetic loci whose elevated rates of recombination mediated diversification, making them strong candidates for being mediators of host-pathogen interaction in diverse strains of C. difficile.

IMPORTANCE Infections with C. difficile result in up to half a million illnesses and tens of thousands of deaths annually in the United States. The severity of C. difficile illness is dependent on both host and bacterial factors. Studying the evolutionary history of C. difficile pathogens is important for understanding the variation in pathogenicity of these bacteria. This study examines the extent and targets of homologous recombination, a mechanism by which distant strains of bacteria can share genetic material, in hundreds of C. difficile strains and identifies hot spots of realized recombination events. The results of this analysis reveal the importance of homologous recombination in the diversification of genetic loci in C. difficile that are significant in its pathogenicity and host interactions, such as flagellar construction, cell wall proteins, and sugar transport and metabolism.

Reduction of non-typeable results using a plasmid oriented Lymfogranuloma venereum PCR for typing of Chlamydia trachomatis positive samples

Objectives

Typing of Chlamydia trachomatis (CT) is traditionally performed by characterising the ompA gene, resulting in more than a dozen different genovars, A to L. Type L is associated with Lymphogranuloma venereum (LGV) and commonly screened for using PCR, targeting the chromosomal pmpH gene. We aimed to develop and validate a new CT/LGV plasmid-based typing assay targeting the pgp3 gene, to increase sensitivity and thus reduce the number of non-typeable results.

Methods

The new pgp3 PCR assay using LNA probes to detect point mutations was analytically and prospectively validated in a routine diagnostic laboratory setting. For the analytical tests, quantified nucleotide constructs (gBlocks) were used to perform limit of detection analyses. Quality control panel samples from 2018 and 2019 for CT were also tested. For the clinical study patient samples which were collected in two months in 2018 were tested simultaneously using the pmpH PCR and the pgp3 PCR.

Results

Analytically, the assay proved to be 100% specific relative to the previously used LGV typing assay targeting the single copy pmpH gene but it was much more sensitive to detect non-LGV CT. In the quality control panel 2 nonLGV samples and 7 LGV samples were solely positive with the pgp3 PCR and not with the pmpH PCR. None of the samples from analytical specificity panels were positive, indicating 100% specificity. In a prospective panel of 152 clinical samples, 142 (93%) were successfully typed with the pgp3 PCR compared to 78% with the pmpH PCR. The pgp3 PCR was fully concordant with the pmpH PCR to identify all LGV subtypes and detected an increased number of clinical samples of non-LGV subtype.

Conclusion

We developed and validated a sensitive and specific plasmid-based typing assay to discriminate LGV from non-LGV CT subtypes. This is useful in a clinical setting to quickly determine the optimal treatment for Chlamydia trachomatis infections.

Clinical Persistence of Chlamydia trachomatis Sexually Transmitted Strains Involves Novel Mutations in the Functional αββα Tetramer of the Tryptophan Synthase Operon

Clinical persistence of Chlamydia trachomatis (Ct) sexually transmitted infections (STIs) is a major public health concern. In vitro persistence is known to develop through interferon gamma (IFN-γ) induction of indoleamine 2,3-dioxygenase (IDO), which catabolizes tryptophan, an essential amino acid for Ct replication. The organism can recover from persistence by synthesizing tryptophan from indole, a substrate for the enzyme tryptophan synthase. The majority of Ct strains, except for reference strain B/TW-5/OT, contain an operon comprised of α and β subunits that encode TrpA and TrpB, respectively, and form a functional αββα tetramer. However, trpA mutations in ocular Ct strains, which are responsible for the blinding eye disease known as trachoma, abrogate tryptophan synthesis from indole. We examined serial urogenital samples from a woman who had recurrent Ct infections over 4 years despite antibiotic treatment. The Ct isolates from each infection episode were genome sequenced and analyzed for phenotypic, structural, and functional characteristics. All isolates contained identical mutations in trpA and developed aberrant bodies within intracellular inclusions, visualized by transmission electron microscopy, even when supplemented with indole following IFN-γ treatment. Each isolate displayed an altered αββα structure, could not synthesize tryptophan from indole, and had significantly lower trpBA expression but higher intracellular tryptophan levels compared with those of reference Ct strain F/IC-Cal3. Our data indicate that emergent mutations in the tryptophan operon, which were previously thought to be restricted only to ocular Ct strains, likely resulted in in vivo persistence in the described patient and represents a novel host-pathogen adaptive strategy for survival.IMPORTANCEChlamydia trachomatis (Ct) is the most common sexually transmitted bacterium with more than 131 million cases occurring annually worldwide. Ct infections are often asymptomatic, persisting for many years despite treatment. In vitro recovery from persistence occurs when indole is utilized by the organism's tryptophan synthase to synthesize tryptophan, an essential amino acid for replication. Ocular but not urogenital Ct strains contain mutations in the synthase that abrogate tryptophan synthesis. Here, we discovered that the genomes of serial isolates from a woman with recurrent, treated Ct STIs over many years were identical with a novel synthase mutation. This likely allowed long-term in vivo persistence where active infection resumed only when tryptophan became available. Our findings indicate an emerging adaptive host-pathogen evolutionary strategy for survival in the urogenital tract that will prompt the field to further explore chlamydial persistence, evaluate the genetics of mutant Ct strains and fitness within the host, and their implications for disease pathogenesis.

The Hidden Genomics of Chlamydia trachomatis

The application of whole-genome sequencing has moved us on from sequencing single genomes to defining unravelling population structures in different niches, and at the -species, -serotype or even -genus level, and in local, national and global settings. This has been instrumental in cataloguing and revealing a huge a range of diversity in this bacterium, when at first we thought there was little. Genomics has challenged assumptions, added insight, as well as confusion and glimpses of truths. What is clear is that at a time when we start to realise the extent and nature of the diversity contained within a genus or a species like this, the huge depth of knowledge communities have developed, through cell biology, as well as the new found molecular approaches will be more precious than ever to link genotype to phenotype. Here we detail the technological developments and insights we have seen during the relatively short time since we began to see the hidden genome of Chlamydia trachomatis.

Evolutionary Model of Cluster Divergence of the Emergent Marine Pathogen Vibrio vulnificus: From Genotype to Ecotype

Vibrio vulnificus is an emergent marine pathogen and is the cause of a deadly septicemia. However, the genetic factors that differentiate its clinical and environmental strains and its several biotypes remain mostly enigmatic. In this work, we investigated the underlying genomic properties and population dynamics of the V. vulnificus species to elucidate the traits that make these strains emerge as a human pathogen. The acquisition of different ecological determinants could have allowed the development of highly divergent clusters with different lifestyles within the same environment. However, we identified strains from both clusters in the mucosa of aquaculture species, indicating that manmade niches are bringing strains from the two clusters together, posing a potential risk of recombination and of emergence of novel variants. We propose a new evolutionary model that provides a perspective that could be broadly applicable to other pathogenic vibrios and facultative bacterial pathogens to pursue strategies to prevent their infections.

Vibrio vulnificus, an opportunistic pathogen, is the causative agent of a life-threatening septicemia and a rising problem for aquaculture worldwide. The genetic factors that differentiate its clinical and environmental strains remain enigmatic. Furthermore, clinical strains have emerged from every clade of V. vulnificus. In this work, we investigated the underlying genomic properties and population dynamics of the V. vulnificus species from an evolutionary and ecological point of view. Genome comparisons and bioinformatic analyses of 113 V. vulnificus isolates indicate that the population of V. vulnificus is made up of four different clusters. We found evidence that recombination and gene flow between the two largest clusters (cluster 1 [C1] and C2) have drastically decreased to the point where they are diverging independently. Pangenome and phenotypic analyses showed two markedly different lifestyles for these two clusters, indicating commensal (C2) and bloomer (C1) ecotypes, with differences in carbohydrate utilization, defense systems, and chemotaxis, among other characteristics. Nonetheless, we identified frequent intra- and interspecies exchange of mobile genetic elements (e.g., antibiotic resistance plasmids, novel “chromids,” or two different and concurrent type VI secretion systems) that provide high levels of genetic diversity in the population. Surprisingly, we identified strains from both clusters in the mucosa of aquaculture species, indicating that manmade niches are bringing strains from the two clusters together. We propose an evolutionary model of V. vulnificus that could be broadly applicable to other pathogenic vibrios and facultative bacterial pathogens to pursue strategies to prevent their infections and emergence.

CtGEM typing: Discrimination of Chlamydia trachomatis ocular and urogenital strains and major evolutionary lineages by high resolution melting analysis of two amplified DNA fragments

Chlamydia trachomatis infects the urogenital tract (UGT) and eyes. Anatomical tropism is correlated with variation in the major outer membrane protein encoded by ompA. Strains possessing the ocular ompA variants A, B, Ba and C are typically found within the phylogenetically coherent “classical ocular lineage”. However, variants B, Ba and C have also been found within three distinct strains in Australia, all associated with ocular disease in children and outside the classical ocular lineage. CtGEM genotyping is a method for detecting and discriminating ocular strains and also the major phylogenetic lineages. The rationale was facilitation of surveillance to inform responses to C. trachomatis detection in UGT specimens from young children. CtGEM typing is based on high resolution melting analysis (HRMA) of two PCR amplified fragments with high combinatorial resolving power, as defined by computerised comparison of 65 whole genomes. One fragment is from the hypothetical gene defined by Jali-1891 in the C. trachomatis B_Jali20 genome, while the other is from ompA. Twenty combinatorial CtGEM types have been shown to exist, and these encompass unique genotypes for all known ocular strains, and also delineate the TI and T2 major phylogenetic lineages, identify LGV strains and provide additional resolution beyond this. CtGEM typing and Sanger sequencing were compared with 42 C. trachomatis positive clinical specimens, and there were no disjunctions. CtGEM typing is a highly efficient method designed and tested using large scale comparative genomics. It divides C. trachomatis into clinically and biologically meaningful groups, and may have broad application in surveillance.

Horizontal gene transfer of Chlamydia: Novel insights from tree reconciliation

Recent comparative genomics studies have suggested that horizontal gene transfer (HGT) is one of the major processes in bacterial evolution. In this study, HGT events of 64 Chlamydia strains were investigated based on the pipeline employed in HGTree database constructed in our recent study. Tree reconciliation method was applied in order to calculate feasible HGT events. Following initial detection and an evaluation procedure, evidence of the HGT was identified in 548 gene families including 42 gene families transferred from outside of Chlamydiae phylum with high reliability. The donor species of inter-phylum HGT consists of 12 different bacterial and archaeal phyla, suggesting that Chlamydia might have even more various host range than in previous reports. In addition, each species of Chlamydia showed varying preference towards HGT, and genes engaged in HGT within Chlamydia and between other species showed different functional distribution. Also, examination of individual gene flows of niche-specific genes suggested that many of such genes are transferred mainly within Chlamydia genus. Our results uncovered novel features of HGT acting on Chlamydia genome evolution, and it would be also strong evidence that HGT is an ongoing process for intracellular pathogens. We expect that the results provide more insight into lineage- and niche-specific adaptations regarding their infectivity and pathogenicity.

Environmental Metagenomic Assemblies Reveal Seven New Highly Divergent Chlamydial Lineages and Hallmarks of a Conserved Intracellular Lifestyle

The Chlamydiae phylum exclusively encompasses bacteria sharing a similar obligate intracellular life cycle. Existing 16S rDNA data support a high diversity within the phylum, however genomic data remain scarce owing to the difficulty in isolating strains using culture systems with eukaryotic cells. Yet, Chlamydiae genome data extracted from large scale metagenomic studies might help fill this gap. This work compares 33 cultured and 27 environmental, uncultured chlamydial genomes, in order to clarify the phylogenetic relatedness of the new chlamydial clades and to investigate the genetic diversity of the Chlamydiae phylum. The analysis of published chlamydial genomes from metagenomics bins and single cell sequencing allowed the identification of seven new deeply branching chlamydial clades sharing genetic hallmarks of parasitic Chlamydiae. Comparative genomics suggests important biological differences between those clades, including loss of many proteins involved in cell division in the genus Similichlamydia, and loss of respiratory chain and tricarboxylic acid cycle in several species. Comparative analyses of chlamydial genomes with two proteobacterial orders, the Rhizobiales and the Rickettsiales showed that genomes of different Rhizobiales families are much more similar than genomes of different Rickettsiales families. On the other hand, the chlamydial 16S rRNAs exhibit a higher sequence conservation than their Rickettsiales counterparts, while chlamydial proteins exhibit increased sequence divergence. Studying the diversity and genome plasticity of the entire Chlamydiae phylum is of major interest to better understand the emergence and evolution of this ubiquitous and ancient clade of obligate intracellular bacteria.

Genomic analyses of the Chlamydia trachomatis core genome show an association between chromosomal genome, plasmid type and disease

Background

Chlamydia trachomatis (Ct) plasmid has been shown to encode genes essential for infection. We evaluated the population structure of Ct using whole-genome sequence data (WGS). In particular, the relationship between the Ct genome, plasmid and disease was investigated.

Results

WGS data from 157 Ct isolates deposited in the Chlamydiales pubMLST database (http://pubMLST.org/chlamydiales/) were annotated with 902 genes including the core and accessory genome. Plasmid associated genes were annotated and a plasmid MLST scheme was defined allowing plasmid sequence types to be determined. Plasmid allelic variation was investigated. Phylogenetic relationships were examined using the Genome Comparator tool available in pubMLST. Phylogenetic analyses identified four distinct Ct core genome clusters and six plasmid clusters, with a strong association between the chromosomal genotype and plasmid. This in turn was linked to ompA genovars and disease phenotype. Horizontal genetic transfer of plasmids was observed for three urogenital-associated isolates, which possessed plasmids more commonly found in isolates resulting from ocular infections. The pgp3 gene was identified as the most polymorphic plasmid gene and pgp4 was the most conserved.

Conclusion

A strong association between chromosomal genome, plasmid type and disease was observed, consistent with previous studies. This suggests co-evolution of the Ct chromosome and their plasmids, but we confirmed that plasmid transfer can occur between isolates. These data provide a better understanding of the genetic diversity occurring across the Ct genome in association with the plasmid content.

Electronic supplementary material

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

Impact of Recombination on the Base Composition of Bacteria and Archaea

The mutational process in bacteria is biased toward A and T, and most species are GC-rich relative to the mutational input to their genome. It has been proposed that the shift in base composition is an adaptive process-that natural selection operates to increase GC-contents-and there is experimental evidence that bacterial strains with GC-rich versions of genes have higher growth rates than those strains with AT-rich versions expressing identical proteins. Alternatively, a nonadaptive process, GC-biased gene conversion (gBGC), could also increase the GC-content of DNA due to the mechanistic bias of gene conversion events during recombination. To determine what role recombination plays in the base composition of bacterial genomes, we compared the spectrum of nucleotide polymorphisms introduced by recombination in all microbial species represented by large numbers of sequenced strains. We found that recombinant alleles are consistently biased toward A and T, and that the magnitude of AT-bias introduced by recombination is similar to that of mutations. These results indicate that recombination alone, without the intervention of selection, is unlikely to counteract the AT-enrichment of bacterial genomes.

Chlamydia trachomatis Strain Types Have Diversified Regionally and Globally with Evidence for Recombination across Geographic Divides

Chlamydia trachomatis (Ct) is the leading cause of bacterial sexually transmitted diseases worldwide. The Ct Multi Locus Sequence Typing (MLST) scheme is effective in differentiating strain types (ST), deciphering transmission patterns and treatment failure, and identifying recombinant strains. Here, we analyzed 323 reference and clinical samples, including 58 samples from Russia, an area that has not previously been represented in Ct typing schemes, to expand our knowledge of the global diversification of Ct STs. The 323 samples resolved into 84 unique STs, a 3.23 higher typing resolution compared to the gold standard single locus ompA genotyping. Our MLST scheme showed a high discriminatory index, D, of 0.98 (95% CI 0.97-0.99) confirming the validity of this method for typing. Phylogenetic analyses revealed distinct branches for the phenotypic diseases of lymphogranuloma venereum, urethritis and cervicitis, and a sub-branch for ocular trachoma. Consistent with these findings, single nucleotide polymorphisms were identified that significantly correlated with each phenotype. While the overall number of unique STs per region was comparable across geographies, the number of STs was greater for Russia with a significantly higher ST/sample ratio of 0.45 (95% CI: 0.35-0.53) compared to Europe or the Americas (p < 0.009), which may reflect a higher level of sexual mixing with the introduction of STs from other regions and/or reassortment of alleles. Four STs were found to be significantly associated with a particular geographic region. ST23 [p = 0.032 (95% CI: 1-23)], ST34 [p = 0.019 (95% CI: 1.1-25)]; and ST19 [p = 0.001 (95% CI: 1.7-34.7)] were significantly associated with Netherlands compared to Russia or the Americas, while ST 30 [p = 0.031 (95% CI: 1.1-17.8)] was significantly associated with the Americas. ST19 was significantly associated with Netherlands and Russia compared with the Americans [p = 0.001 (95% CI: 1.7-34.7) and p = 0.006 (95% CI: 1.5-34.6), respectively]. Additionally, recombinant strains were ubiquitous in the data set [106 (32.8%)], although Europe had a significantly higher number than Russia or the Americas (p < 0.04), the majority of which were from Amsterdam [43 (87.8%) of 49)]. The higher number of recombinants in Europe indicates selective pressure and/or adaptive diversification that will require additional studies to elucidate.

Genome-wide profiling of humoral immunity and pathogen genes under selection identifies immune evasion tactics of Chlamydia trachomatis during ocular infection

The frequency and duration of Chlamydia trachomatis (Ct) ocular infections decrease with age, suggesting development of partial immunity. However, there is a lack of clear correlates of immunity to Ct infection in humans. We screened sera from a cohort of Gambian children followed for six-months against a Ct-proteome microarray. At genome sequence level, we detected signatures of selection from a population of ocular Ct isolates from Guinea-Bissau. Together these approaches allowed us to highlight the focus of humoral responses and hypothesise new modes of pathogen immune evasion. Children who were susceptible to frequent and/or prolonged Ct infection had a less focussed antibody response, including preferential recognition of forty-two antigens. There was evidence of positive and purifying selection across the genome, but little balancing selection. In contrast, most antigens that were associated with susceptibility were under neutral selection. These data suggest an evasion strategy in which Ct presents a large panel of irrelevant antigens to the immune system to block or misdirect protective responses. Development of a focused immune response, possibly induced through vaccination, may be an effective strategy to promote protection to Ct infection.

Molecular Genetic Analysis of Chlamydia Species

Species of Chlamydia are the etiologic agent of endemic blinding trachoma, the leading cause of bacterial sexually transmitted diseases, significant respiratory pathogens, and a zoonotic threat. Their dependence on an intracellular growth niche and their peculiar developmental cycle are major challenges to elucidating their biology and virulence traits. The last decade has seen tremendous advances in our ability to perform a molecular genetic analysis of Chlamydia species. Major achievements include the generation of large collections of mutant strains, now available for forward- and reverse-genetic applications, and the introduction of a system for plasmid-based transformation enabling complementation of mutations; expression of foreign, modified, or reporter genes; and even targeted gene disruptions. This review summarizes the current status of the molecular genetic toolbox for Chlamydia species and highlights new insights into their biology and new challenges in the nascent field of Chlamydia genetics.

Comprehensive global genome dynamics of Chlamydia trachomatis show ancient diversification followed by contemporary mixing and recent lineage expansion

Chlamydia trachomatis is the world's most prevalent bacterial sexually transmitted infection and leading infectious cause of blindness, yet it is one of the least understood human pathogens, in part due to the difficulties of in vitro culturing and the lack of available tools for genetic manipulation. Genome sequencing has reinvigorated this field, shedding light on the contemporary history of this pathogen. Here, we analyze 563 full genomes, 455 of which are novel, to show that the history of the species comprises two phases, and conclude that the currently circulating lineages are the result of evolution in different genomic ecotypes. Temporal analysis indicates these lineages have recently expanded in the space of thousands of years, rather than the millions of years as previously thought, a finding that dramatically changes our understanding of this pathogen's history. Finally, at a time when almost every pathogen is becoming increasingly resistant to antimicrobials, we show that there is no evidence of circulating genomic resistance in C. trachomatis.

Rapid detection and strain typing of Chlamydia trachomatis using a highly multiplexed microfluidic PCR assay

Nucleic acid amplification tests (NAATs) are recommended by the CDC for detection of Chlamydia trachomatis (Ct) urogenital infections. Current commercial NAATs require technical expertise and sophisticated laboratory infrastructure, are time-consuming and expensive, and do not differentiate the lymphogranuloma venereum (LGV) strains that require a longer duration of treatment than non-LGV strains. The multiplexed microfluidic PCR-based assay presented in this work simultaneously interrogates 13 loci to detect Ct and identify LGV and non-LGV strain-types. Based on amplified fragment length polymorphisms, the assay differentiates LGV, ocular, urogenital, and proctocolitis clades, and also serovars L1, L2, and L3 within the LGV group. The assay was evaluated in a blinded fashion using 95 clinical swabs, with 76 previously reported as urogenital Ct-positive samples and typed by ompA genotyping and/or Multi-Locus Sequence Typing. Results of the 13-plex assay showed that 51 samples fell within urogenital clade 2 or 4, 24 samples showed both clade 2 and 4 signatures, indicating possible mixed infection, gene rearrangement, or inter-clade recombination, and one sample was a noninvasive trachoma biovar (either a clade 3 or 4). The remaining 19 blinded samples were correctly identified as LGV clade 1 (3), ocular clade 3 (4), or as negatives (12). To date, no NAAT assay can provide a point-of-care applicable turnaround time for Ct detection while identifying clinically significant Ct strain types to inform appropriate treatment. Coupled with rapid DNA processing of clinical swabs (approximately 60 minutes from swab-in to result-out), the assay has significant potential as a rapid POC diagnostic for Ct infections.

European Chlamydia abortus livestock isolate genomes reveal unusual stability and limited diversity, reflected in geographical signatures

Background

Chlamydia abortus (formerly Chlamydophila abortus) is an economically important livestock pathogen, causing ovine enzootic abortion (OEA), and can also cause zoonotic infections in humans affecting pregnancy outcome. Large-scale genomic studies on other chlamydial species are giving insights into the biology of these organisms but have not yet been performed on C. abortus. Our aim was to investigate a broad collection of European isolates of C. abortus, using next generation sequencing methods, looking at diversity, geographic distribution and genome dynamics.

Results

Whole genome sequencing was performed on our collection of 57 C. abortus isolates originating primarily from the UK, Germany, France and Greece, but also from Tunisia, Namibia and the USA. Phylogenetic analysis of a total of 64 genomes shows a deep structural division within the C. abortus species with a major clade displaying limited diversity, in addition to a branch carrying two more distantly related Greek isolates, LLG and POS. Within the major clade, seven further phylogenetic groups can be identified, demonstrating geographical associations. The number of variable nucleotide positions across the sampled isolates is significantly lower than those published for C. trachomatis and C. psittaci. No recombination was identified within C. abortus, and no plasmid was found. Analysis of pseudogenes showed lineage specific loss of some functions, notably with several Pmp and TMH/Inc proteins predicted to be inactivated in many of the isolates studied.

Conclusions

The diversity within C. abortus appears to be much lower compared to other species within the genus. There are strong geographical signatures within the phylogeny, indicating clonal expansion within areas of limited livestock transport. No recombination has been identified within this species, showing that different species of Chlamydia may demonstrate different evolutionary dynamics, and that the genome of C. abortus is highly stable.

Electronic supplementary material

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

Molecular typing of Chlamydia trachomatis: An overview

Urogenital infection due to Chlamydia trachomatis (CT) is one of the most common bacterial sexually transmitted infections (STIs) and is a major public health problem worldwide. Molecular characterisation of CT is important for understanding the pathophysiological mechanisms of chlamydial disease and its transmission dynamics in sexual networks. Traditionally, strain typing of CT was based on serotyping methods characterising the major outer membrane protein (MOMP). With the advent of polymerase chain reaction and sequencing the era of molecular typing began. Molecular characterization of CT strains is based on sequence analysis of ompA gene encoding MOMP. However, in due course of time, improvements were made to enhance the discriminatory power of sequencing and quality of epidemiological information. New high-resolution genotyping methods using multiple loci such as multilocus sequence typing (MLST) and multiple loci variable number of tandem repeats (MLVA) were developed but were unable to differentiate mixed infections (MIs). The development of DNA-hybridisation methods emerged as a major breakthrough in detecting MIs. Although MLST and MLVA are more discriminative than other genotyping methods, they are laborious and expensive. DNA microarray technique is an affordable alternative for genotyping. Since recombination is widespread in the CT genome, ompA is not a reliable marker for phylogenetic studies; hence, whole genome sequencing may provide maximum phylogenetic resolution of CT strains. A descriptive review is provided of the various molecular CT typing methods. The vital information gained can be used for formulating screening programmes, targeted prevention and optimising therapeutic measures aiming to reduce disease transmission.

Spread of a new Chlamydia trachomatis variant from men who have sex with men to the heterosexual population after replacement and recombination in ompA and pmpH genes

OBJECTIVES

Sexually transmitted infections are frequently related to outbreaks in high-risk populations due to the dense sexual networks. We wanted to determine the dissemination of a Chlamydia trachomatis variant characterized by the pmpH-recombinant gene between L and G genotypes, which was previously described in a high-risk population.

METHODS

A total of 449 samples were analysed in two periods ranging from 2009 to 2015 for detection of the pmpH-recombinant gene. For those samples yielding positive amplification, a sampling was selected for phylogenetic reconstructions based on sequencing of five chromosomal genes.

RESULTS

Globally this variant was found in 113 of the 449 samples (25%). During the first years (2009-13), this variant was found almost exclusively in rectal samples (30/112 samples) of men who have sex with men and in only one non-rectal sample (1/63). In 2014, this variant was also found in urethral and pharyngeal samples (1/24 and 1/7, respectively). However, in 2015, an epidemiological change was observed as the proportion of this variant had increased in rectal samples (20/51; 39%) and non-rectal samples, including cervical samples (51/142; 36.4%). The molecular characterization revealed the replacement of the ompA gene belonging to subtype G in samples recovered from 2009 to 2013 by the ompA gene belonging to subtype J after 2013.

CONCLUSIONS

Our data would support the evidence that subtype J could be a 'subtype bridge' between different sexual networks, as subtype J has been found in men who have sex with men and heterosexual populations in similar proportions. This work reveals the necessity of implementing molecular surveillance in extra-rectal samples to help us understand the gaps in transmission.

The Chlamydia suis Genome Exhibits High Levels of Diversity, Plasticity, and Mobile Antibiotic Resistance: Comparative Genomics of a Recent Livestock Cohort Shows Influence of Treatment Regimes

Chlamydia suis is an endemic pig pathogen, belonging to a fascinating genus of obligate intracellular pathogens. Of particular interest, this is the only chlamydial species to have naturally acquired genes encoding for tetracycline resistance. To date, the distribution and mobility of the Tet-island are not well understood. Our study focused on whole genome sequencing of 29 C. suis isolates from a recent porcine cohort within Switzerland, combined with data from USA tetracycline-resistant isolates. Our findings show that the genome of C. suis is very plastic, with unprecedented diversity, highly affected by recombination and plasmid exchange. A large diversity of isolates circulates within Europe, even within individual Swiss farms, suggesting that C. suis originated around Europe. New World isolates have more restricted diversity and appear to derive from European isolates, indicating that historical strain transfers to the United States have occurred. The architecture of the Tet-island is variable, but the tetA(C) gene is always intact, and recombination has been a major factor in its transmission within C. suis. Selective pressure from tetracycline use within pigs leads to a higher number of Tet-island carrying isolates, which appear to be lost in the absence of such pressure, whereas the loss or gain of the Tet-island from individual strains is not observed. The Tet-island appears to be a recent import into the genome of C. suis, with a possible American origin.

Tet(C) Gene Transfer between Chlamydia suis Strains Occurs by Homologous Recombination after Co-infection: Implications for Spread of Tetracycline-Resistance among Chlamydiaceae

Chlamydia suis is a swine pathogen that has also recently been found to cause zoonotic infections of the human eye, pharynx, and gastrointestinal tract. Many strains contain a tetracycline class C gene [tet(C)] cassette that confers tetracycline resistance. The cassette was likely originally acquired by horizontal gene transfer from a Gram-negative donor after the introduction of tetracycline into animal feed in the 1950s. Various research groups have described the capacity for different Chlamydia species to exchange DNA by homologous recombination. Since over 90% of C. suis strains are tetracycline resistant, they represent a potential source for antibiotic-resistance spread within and between Chlamydiaceae species. Here, we examined the genetics of tet(C)-transfer among C. suis strains. Tetracycline-sensitive C. suis strain S45 was simultaneously or sequentially co-infected with tetracycline-resistant C. suis strains in McCoy cells. Potential recombinants were clonally purified by a harvest assay derived from the classic plaque assay. C. suis strain Rogers132, lacking transposases IS200 and IS605, was the most efficient donor, producing two unique recombinants detected in three of the 56 (5.4%) clones screened. Recombinants were found to have a minimal inhibitory concentration (MIC) of 8-16 μg/mL for tetracycline. Resistance remained stable over 10 passages as long as recombinants were initially grown in tetracycline at twice the MIC of S45 (0.032 μg/mL). Genomic analysis revealed that tet(C) had integrated into the S45 genome by homologous recombination at two unique sites depending on the recombinant: a 55 kb exchange between nrqF and pckG, and a 175 kb exchange between kdsA and cysQ. Neither site was associated with inverted repeats or motifs associated with recombination hotspots. Our findings show that cassette transfer into S45 has low frequency, does not require IS200/IS605 transposases, is stable if initially grown in tetracycline, and results in multiple genomic configurations. We provide a model for stable cassette transfer to better understand the capability for cassette acquisition by Chlamydiaceae species that infect humans, a matter of public health importance.

Reconstructing the Ancestral Relationships Between Bacterial Pathogen Genomes

Following recent developments in DNA sequencing technology, it is now possible to sequence hundreds of whole genomes from bacterial isolates at relatively low cost. Analyzing this growing wealth of genomic data in terms of ancestral relationships can reveal many interesting aspects of the evolution, ecology, and epidemiology of bacterial pathogens. However, reconstructing the ancestry of a sample of bacteria remains challenging, especially for the majority of species where recombination is frequent. Here, we review and describe the computational techniques currently available to infer ancestral relationships, including phylogenetic methods that either ignore or account for the effect of recombination, as well as model-based and model-free phylogeny-independent approaches.

Dynamics of genome change among Legionella species

Legionella species inhabit freshwater and soil ecosystems where they parasitize protozoa. L. pneumonphila (LP) serogroup-1 (Lp1) is the major cause of Legionnaires' Disease (LD), a life-threatening pulmonary infection that can spread systemically. The increased global frequency of LD caused by Lp and non-Lp species underscores the need to expand our knowledge of evolutionary forces underlying disease pathogenesis. Whole genome analyses of 43 strains, including all known Lp serogroups 1-17 and 17 emergent LD-causing Legionella species (of which 33 were sequenced in this study) in addition to 10 publicly available genomes, resolved the strains into four phylogenetic clades along host virulence demarcations. Clade-specific genes were distinct for genetic exchange and signal-transduction, indicating adaptation to specific cellular and/or environmental niches. CRISPR spacer comparisons hinted at larger pools of accessory DNA sequences in Lp than predicted by the pan-genome analyses. While recombination within Lp was frequent and has been reported previously, population structure analysis identified surprisingly few DNA admixture events between species. In summary, diverse Legionella LD-causing species share a conserved core-genome, are genetically isolated from each other, and selectively acquire genes with potential for enhanced virulence.

Chlamydia trachomatis Genital Infections

Etiology, transmission and protection: Chlamydia trachomatis is the leading cause of bacterial sexually transmitted infection (STI) globally. However, C. trachomatis also causes trachoma in endemic areas, mostly Africa and the Middle East, and is a leading cause of preventable blindness worldwide. Epidemiology, incidence and prevalence: The World Health Organization estimates 131 million new cases of C. trachomatis genital infection occur annually. Globally, infection is most prevalent in young women and men (14-25 years), likely driven by asymptomatic infection, inadequate partner treatment and delayed development of protective immunity. Pathology/Symptomatology: C. trachomatis infects susceptible squamocolumnar or transitional epithelial cells, leading to cervicitis in women and urethritis in men. Symptoms are often mild or absent but ascending infection in some women may lead to Pelvic Inflammatory Disease (PID), resulting in reproductive sequelae such as ectopic pregnancy, infertility and chronic pelvic pain. Complications of infection in men include epididymitis and reactive arthritis. Molecular mechanisms of infection: Chlamydiae manipulate an array of host processes to support their obligate intracellular developmental cycle. This leads to activation of signaling pathways resulting in disproportionate influx of innate cells and the release of tissue damaging proteins and pro-inflammatory cytokines. Treatment and curability: Uncomplicated urogenital infection is treated with azithromycin (1 g, single dose) or doxycycline (100 mg twice daily x 7 days). However, antimicrobial treatment does not ameliorate established disease. Drug resistance is rare but treatment failures have been described. Development of an effective vaccine that protects against upper tract disease or that limits transmission remains an important goal.

Tetracycline Selective Pressure and Homologous Recombination Shape the Evolution of Chlamydia suis: A Recently Identified Zoonotic Pathogen

Species closely related to the human pathogen Chlamydia trachomatis (Ct) have recently been found to cause zoonotic infections, posing a public health threat especially in the case of tetracycline resistant Chlamydia suis (Cs) strains. These strains acquired a tet(C)-containing cassette via horizontal gene transfer (HGT). Genomes of 11 Cs strains from various tissues were sequenced to reconstruct evolutionary pathway(s) for tet(C) HGT. Cs had the highest recombination rate of Chlamydia species studied to date. Admixture occurred among Cs strains and with Chlamydia muridarum but not with Ct Although in vitro tet(C) cassette exchange with Ct has been documented, in vivo evidence may require examining human samples from Ct and Cs co-infected sites. Molecular-clock dating indicated that ancestral clades of resistant Cs strains predated the 1947 discovery of tetracycline, which was subsequently used in animal feed. The cassette likely spread throughout Cs strains by homologous recombination after acquisition from an external source, and our analysis suggests Betaproteobacteria as the origin. Selective pressure from tetracycline may be responsible for recent bottlenecks in Cs populations. Since tetracycline is an important antibiotic for treating Ct, zoonotic infections at mutual sites of infection indicate the possibility for cassette transfer and major public health repercussions.

Molecular characteristics of the ompA gene of serotype B Chlamydia trachomatis in Qinghai Tibetan primary school students

To study the molecular characteristics of Chlamydia trachomatis, the major outer membrane protein gene (ompA) of C. trachomatis from primary school students with trachoma residing in the Qinghai Tibetan area was sequenced and compared with the same serotype in GenBank. In Jianshetang Primary School and Galeng Central Primary School in the Galeng Tibetan Township of Qinghai Haidong Sala Autonomous County, scraped samples were collected from the upper tarsal conjunctiva and lower conjunctival sac of both eyes of 45 students with trachoma, stored at 4°C, and transported to Beijing Tongren Hospital by air within 24 h. The samples were screened for C. trachomatis by real-time PCR. The ompA gene from the C. trachomatis-positive samples was amplified by nested PCR. The serotype was confirmed by National Center for Biotechnology Information (NCBI) BLAST search and homology analysis. The entire ompA gene sequence was compared with the corresponding gene sequences of serotype B strains available in GenBank. Of the 45 students aged 6-13 years with trachoma, 26 C. trachomatis-positive students were identified by the initial real-time PCR screening (average age, (9.09±1.63) years; sex ratio, 1.0), accounting for 57.78% (26/45). The cycle threshold values for real-time PCR were 16.79-37.77. Half (13/26) of C. trachomatis-positive students had a bacterial copy number of >10(5). The compliance rate of the ompA gene sequences with the C. trachomatis serotype B strains in GenBank was up to 99%. Two novel genetic mutations were found when the ompA gene was compared with those of the 11 serotype B strains in GenBank. The two non-synonymous mutations were located at (i) position 271 in the second constant domain, an adenine (A) to guanine (G) substitution (ACT→GCT), changing the amino acid at position 91 from threonine to alanine (Thr→Ala) in all 26 strains; and (ii) position 887 in the fourth variable domain, a cytosine (C) to thymine (T) substitution (GCA→GTA), changing the amino acid at residue 296 from alanine to valine (Ala→Val) in four of the 26 strains. Six mutations were identified relative to ATCC VR-573. The strains could be divided into two gene clusters according to the mutation at nucleotide position 887: CQZ-1 (China Qinghai Tibetan-1) and CQZ-2 (China Qinghai Tibetan-2). We thus detected two novel serotype B mutant strains of C. trachomatis among study subjects with trachoma.

Space: A Final Frontier for Vacuolar Pathogens

There is a fundamental gap in our understanding of how a eukaryotic cell apportions the limited space within its cell membrane. Upon infection, a cell competes with intracellular pathogens for control of this same precious resource. The struggle between pathogen and host provides us with an opportunity to uncover the mechanisms regulating subcellular space by understanding how pathogens modulate vesicular traffic and membrane fusion events to create a specialized compartment for replication. By comparing several important intracellular pathogens, we review the molecular mechanisms and trafficking pathways that drive two space allocation strategies, the formation of tight and spacious pathogen-containing vacuoles. Additionally, we discuss the potential advantages of each pathogenic lifestyle, the broader implications these lifestyles might have for cellular biology and outline exciting opportunities for future investigation.

Emancipating Chlamydia: Advances in the Genetic Manipulation of a Recalcitrant Intracellular Pathogen

Chlamydia species infect millions of individuals worldwide and are important etiological agents of sexually transmitted disease, infertility, and blinding trachoma. Historically, the genetic intractability of this intracellular pathogen has hindered the molecular dissection of virulence factors contributing to its pathogenesis. The obligate intracellular life cycle of Chlamydia and restrictions on the use of antibiotics as selectable markers have impeded the development of molecular tools to genetically manipulate these pathogens. However, recent developments in the field have resulted in significant gains in our ability to alter the genome of Chlamydia, which will expedite the elucidation of virulence mechanisms. In this review, we discuss the challenges affecting the development of molecular genetic tools for Chlamydia and the work that laid the foundation for recent advancements in the genetic analysis of this recalcitrant pathogen.

Chlamydia trachomatis from Australian Aboriginal people with trachoma are polyphyletic composed of multiple distinctive lineages

Chlamydia trachomatis causes sexually transmitted infections and the blinding disease trachoma. Current data on C. trachomatis phylogeny show that there is only a single trachoma-causing clade, which is distinct from the lineages causing urogenital tract (UGT) and lymphogranuloma venerum diseases. Here we report the whole-genome sequences of ocular C. trachomatis isolates obtained from young children with clinical signs of trachoma in a trachoma endemic region of northern Australia. The isolates form two lineages that fall outside the classical trachoma lineage, instead being placed within UGT clades of the C. trachomatis phylogenetic tree. The Australian trachoma isolates appear to be recombinants with UGT C. trachomatis genome backbones, in which loci that encode immunodominant surface proteins (ompA and pmpEFGH) have been replaced by those characteristic of classical ocular isolates. This suggests that ocular tropism and association with trachoma are functionally associated with some sequence variants of ompA and pmpEFGH.

The Landscape of Realized Homologous Recombination in Pathogenic Bacteria

Recombination enhances the adaptive potential of organisms by allowing genetic variants to be tested on multiple genomic backgrounds. Its distribution in the genome can provide insight into the evolutionary forces that underlie traits, such as the emergence of pathogenicity. Here, we examined landscapes of realized homologous recombination of 500 genomes from ten bacterial species and found all species have "hot" regions with elevated rates relative to the genome average. We examined the size, gene content, and chromosomal features associated with these regions and the correlations between closely related species. The recombination landscape is variable and evolves rapidly. For example in Salmonella, only short regions of around 1 kb in length are hot whereas in the closely related species Escherichia coli, some hot regions exceed 100 kb, spanning many genes. Only Streptococcus pyogenes shows evidence for the positive correlation between GC content and recombination that has been reported for several eukaryotes. Genes with function related to the cell surface/membrane are often found in recombination hot regions but E. coli is the only species where genes annotated as "virulence associated" are consistently hotter. There is also evidence that some genes with "housekeeping" functions tend to be overrepresented in cold regions. For example, ribosomal proteins showed low recombination in all of the species. Among specific genes, transferrin-binding proteins are recombination hot in all three of the species in which they were found, and are subject to interspecies recombination.

Chlamydiaceae Genomics Reveals Interspecies Admixture and the Recent Evolution of Chlamydia abortus Infecting Lower Mammalian Species and Humans

Chlamydiaceae are obligate intracellular bacteria that cause a diversity of severe infections among humans and livestock on a global scale. Identification of new species since 1989 and emergence of zoonotic infections, including abortion in women, underscore the need for genome sequencing of multiple strains of each species to advance our knowledge of evolutionary dynamics across Chlamydiaceae. Here, we genome sequenced isolates from avian, lower mammalian and human hosts. Based on core gene phylogeny, five isolates previously classified as Chlamydia abortus were identified as members of Chlamydia psittaci and Chlamydia pecorum. Chlamydia abortus is the most recently emerged species and is a highly monomorphic group that lacks the conserved virulence-associated plasmid. Low-level recombination and evidence for adaptation to the placenta echo evolutionary processes seen in recently emerged, highly virulent niche-restricted pathogens, such as Bacillus anthracis. In contrast, gene flow occurred within C. psittaci and other Chlamydiaceae species. The C. psittaci strain RTH, isolated from a red-tailed hawk (Buteo jamaicensis), is an outlying strain with admixture of C. abortus, C. psittaci, and its own population markers. An average nucleotide identity of less than 94% compared with other Chlamydiaceae species suggests that RTH belongs to a new species intermediary between C. psittaci and C. abortus. Hawks, as scavengers and predators, have extensive opportunities to acquire multiple species in their intestinal tract. This could facilitate transformation and homologous recombination with the potential for new species emergence. Our findings indicate that incubator hosts such as birds-of-prey likely promote Chlamydiaceae evolution resulting in novel pathogenic lineages.

Genome expansion in bacteria: the curios case of Chlamydia trachomatis

Background

Recent findings indicated that a correlation between genomic % AT and genome size within strains of microbial species was predominantly associated with the uptake of foreign DNA. One species however, Chlamydia trachomatis, defied any explanation. In the present study 79 fully sequenced C. trachomatis genomes, representing ocular- (nine strains), urogenital- (36 strains) and lymphogranuloma venereum strains (LGV, 22 strains), in three pathogroups, in addition to 12 laboratory isolates, were scrutinized with the intent of elucidating the positive correlation between genomic AT content and genome size.

Results

The average size difference between the strains of each pathogroup was largely explained by the incorporation of genetic fragments. These fragments were slightly more AT rich than their corresponding host genomes, but not enough to justify the difference in AT content between the strains of the smaller genomes lacking the fragments. In addition, a genetic region predominantly found in the ocular strains, which had the largest genomes, was on average more GC rich than the host genomes of the urogenital strains (58.64 % AT vs. 58.69 % AT), which had the second largest genomes, implying that the foreign genetic regions cannot alone explain the association between genome size and AT content in C. trachomatis. 23,492 SNPs were identified for all 79 genomes, and although the SNPs were on average slightly GC rich (~47 % AT), a significant association was found between genome-wide SNP AT content, for each pathogroup, and genome size (p < 0.001, R² = 0.86) in the C. trachomatis strains.

Conclusions

The correlation between genome size and AT content, with respect to the C. trachomatis pathogroups, was explained by the incorporation of genetic fragments unique to the ocular and/or urogenital strains into the LGV- and urogential strains in addition to the genome-wide SNP AT content differences between the three pathogroups.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-015-1464-6) contains supplementary material, which is available to authorized users.

High-resolution typing of Chlamydia trachomatis: epidemiological and clinical uses

PURPOSE OF REVIEW

A state-of-the-art overview of molecular Chlamydia trachomatis typing methods that are used for routine diagnostics and scientific studies.

RECENT FINDINGS

Molecular epidemiology uses high-resolution typing techniques such as multilocus sequence typing, multilocus variable number of tandem repeats analysis, and whole-genome sequencing to identify strains based on their DNA sequence. These data can be used for cluster, network and phylogenetic analyses, and are used to unveil transmission networks, risk groups, and evolutionary pathways. High-resolution typing of C. trachomatis strains is applied to monitor treatment efficacy and re-infections, and to study the recent emergence of lymphogranuloma venereum (LGV) amongst men who have sex with men in high-income countries. Chlamydia strain typing has clinical relevance in disease management, as LGV needs longer treatment than non-LGV C. trachomatis. It has also led to the discovery of a new variant Chlamydia strain in Sweden, which was not detected by some commercial C. trachomatis diagnostic platforms.

SUMMARY

After a brief history and comparison of the various Chlamydia typing methods, the applications of the current techniques are described and future endeavors to extend scientific understanding are formulated. High-resolution typing will likely help to further unravel the pathophysiological mechanisms behind the wide clinical spectrum of chlamydial disease.

Chlamydia trachomatis In Vivo to In Vitro Transition Reveals Mechanisms of Phase Variation and Down-Regulation of Virulence Factors

Research on the obligate intracellular bacterium Chlamydia trachomatis demands culture in cell-lines, but the adaptive process behind the in vivo to in vitro transition is not understood. We assessed the genomic and transcriptomic dynamics underlying C. trachomatis in vitro adaptation of strains representing the three disease groups (ocular, epithelial-genital and lymphogranuloma venereum) propagated in epithelial cells over multiple passages. We found genetic features potentially underlying phase variation mechanisms mediating the regulation of a lipid A biosynthesis enzyme (CT533/LpxC), and the functionality of the cytotoxin (CT166) through an ON/OFF mechanism. We detected inactivating mutations in CT713/porB, a scenario suggesting metabolic adaptation to the available carbon source. CT135 was inactivated in a tropism-specific manner, with CT135-negative clones emerging for all epithelial-genital populations (but not for LGV and ocular populations) and rapidly increasing in frequency (~23% mutants per 10 passages). RNA-sequencing analyses revealed that a deletion event involving CT135 impacted the expression of multiple virulence factors, namely effectors known to play a role in the C. trachomatis host-cell invasion or subversion (e.g., CT456/Tarp, CT694, CT875/TepP and CT868/ChlaDub1). This reflects a scenario of attenuation of C. trachomatis virulence in vitro, which may take place independently or in a cumulative fashion with the also observed down-regulation of plasmid-related virulence factors. This issue may be relevant on behalf of the recent advances in Chlamydia mutagenesis and transformation where culture propagation for selecting mutants/transformants is mandatory. Finally, there was an increase in the growth rate for all strains, reflecting gradual fitness enhancement over time. In general, these data shed light on the adaptive process underlying the C. trachomatis in vivo to in vitro transition, and indicates that it would be prudent to restrict culture propagation to minimal passages and check the status of the CT135 genotype in order to avoid the selection of CT135-negative mutants, likely originating less virulent strains.