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

Chlamydia trachomatis: a model for intracellular bacterial parasitism

Chlamydia comprises a diverse group of obligate intracellular bacteria that cause infections in animals, including humans. These organisms share fascinating biology, including distinct developmental stages, non-canonical cell surface structures, and adaptations to intracellular parasitism. Chlamydia trachomatis is of particular interest due to its significant clinical importance, causing both ocular and sexually transmitted infections. The strain L2/434/Bu, responsible for lymphogranuloma venereum, is the most common strain used to study chlamydial molecular and cell biology because it grows readily in cell culture and is amenable to genetic manipulation. Indeed, this strain has enabled researchers to tackle fundamental questions about the molecular mechanisms underlying Chlamydia's developmental transitions and biphasic lifecycle and cellular adaptations to obligate intracellular parasitism, including characterizing numerous conserved virulence genes and defining immune responses. However, L2/434/Bu is not representative of C. trachomatis strains that cause urogenital infections in humans, limiting its utility in addressing questions of host tropism and immune evasion in reproductive organs. Recent research efforts are shifting toward understanding the unique attributes of more clinically relevant C. trachomatis genovars.

Whole-genome automated assembly pipeline for Chlamydia trachomatis strains from reference, in vitro and clinical samples using the integrated CtGAP pipeline

Whole genome sequencing (WGS) is pivotal for the molecular characterization of Chlamydia trachomatis (Ct)-the leading bacterial cause of sexually transmitted infections and infectious blindness worldwide. Ct WGS can inform epidemiologic, public health and outbreak investigations of these human-restricted pathogens. However, challenges persist in generating high-quality genomes for downstream analyses given its obligate intracellular nature and difficulty with in vitro propagation. No single tool exists for the entirety of Ct genome assembly, necessitating the adaptation of multiple programs with varying success. Compounding this issue is the absence of reliable Ct reference strain genomes. We, therefore, developed CtGAP-Chlamydia trachomatisGenome Assembly Pipeline-as an integrated 'one-stop-shop' pipeline for assembly and characterization of Ct genome sequencing data from various sources including isolates, in vitro samples, clinical swabs and urine. CtGAP, written in Snakemake, enables read quality statistics output, adapter and quality trimming, host read removal, de novo and reference-guided assembly, contig scaffolding, selective ompA, multi-locus-sequence and plasmid typing, phylogenetic tree construction, and recombinant genome identification. Twenty Ct reference genomes were also generated. Successfully validated on a diverse collection of 363 samples containing Ct, CtGAP represents a novel pipeline requiring minimal bioinformatics expertise with easy adaptation for use with other bacterial species.

Chlamydia trachomatis: From Urogenital Infections to the Pathway of Infertility

Chlamydia trachomatis (CT) is a major cause of sexually transmitted infections (STIs) worldwide, with significant implications for reproductive health. The bacterium's genome contains highly polymorphic regions, influencing both the type and severity of infections. These genetic variations, particularly those occurring in the major outer membrane protein (MOMP) gene, are critical for classifying the bacterium into distinct serovars and enable CT to adapt to diverse host environments, contributing to its immune evasion, persistence, and pathogenicity. Persistent or untreated urogenital infections can lead to chronic inflammation, tissue damage, and pelvic inflammatory disease, ultimately increasing the risk of ectopic pregnancy, spontaneous abortion, and infertility. This review consolidates current knowledge on the genetic diversity of CT, its potential role in modulating infection outcomes, and its immune evasion mechanisms. By integrating scientific evidence linking chlamydial infections to infertility, we underscore the urgent need for targeted research to address this critical public health challenge.

Chlamydia trachomatis genomes from rectal samples: description of a new clade comprising ompA-genotype L4 from Argentina

Whole-genome analysis has provided insights into the evolution of Chlamydia trachomatis and, recently, into circulating strains that cause lymphogranuloma venereum (LGV). A large LGV outbreak of a new ompA-genotype, L2b, was first reported in Europe in the early 2000s, primarily affecting men who have sex with men (MSM), and then expanded globally. More recent work shows that this outbreak is diversifying into variants of described ompA-genotypes, with the same L2b genomic backbone. This study extends the investigation of LGV cases to Argentina and Finland. In 2017, an LGV outbreak was described in Argentina characterized by distinct genomic features shown by both ompA-genotyping and Multi-Locus Sequence Typing (MLST) analysis. We have obtained whole-genome sequences from cultured isolates and clinical samples via SureSelect (Agilent) target enrichment. Based on ompA and phylogenetic analyses, we describe further diversity within the ompA-genotype L2b clade, illustrating the transmission dynamics in Argentina and Finland. A key finding is that of a novel clade of Argentinian samples, characterized by a proposed new ompA-genotype L4. Additionally, we present the genome sequence of a non-LGV strain associated with anorectal proctitis. These findings contribute to the investigation of LGV evolution, particularly with the presence of the novel L4 lineage, and provide insights into genomic diversity and transmission dynamics of C. trachomatis.

Pangenome-Wide Association Study in the Chlamydiaceae Family Reveals Key Evolutionary Aspects of Their Relationship with Their Hosts

The Chlamydiaceae are a family of obligate intracellular bacteria known for their unique biphasic developmental cycle. Chlamydial are associated with various host organisms, including humans, and have been proposed as emerging pathogens. Genomic studies have significantly enhanced our understanding of chlamydial biology, host adaptation, and evolutionary processes. In this study, we conducted a complete pangenome association analysis (pan-GWAS) using 101 genomes from the Chlamydiaceae family to identify differentially represented genes in Chlamydia and Chlamydophila, revealing their distinct evolutionary strategies for interacting with eukaryotic hosts. Our analysis identified 289 genes with differential abundance between the two clades: 129 showed a strong association with Chlamydia and 160 with Chlamydophila. Most genes in Chlamydia were related to the type III secretion system, while Chlamydophila genes corresponded to various functional categories, including translation, replication, transport, and metabolism. These findings suggest that Chlamydia has developed a high dependence on mammalian cells for replication, facilitated by a complex T3SS for intracellular manipulation. In contrast, the metabolic and functional diversity in Chlamydophila allows it to colonize a broad range of hosts, such as birds, reptiles, amphibians, and mammals, making it a less specialized clade.

RIBAP: a comprehensive bacterial core genome annotation pipeline for pangenome calculation beyond the species level

Microbial pangenome analysis identifies present or absent genes in prokaryotic genomes. However, current tools are limited when analyzing species with higher sequence diversity or higher taxonomic orders such as genera or families. The Roary ILP Bacterial core Annotation Pipeline (RIBAP) uses an integer linear programming approach to refine gene clusters predicted by Roary for identifying core genes. RIBAP successfully handles the complexity and diversity of Chlamydia, Klebsiella, Brucella, and Enterococcus genomes, outperforming other established and recent pangenome tools for identifying all-encompassing core genes at the genus level. RIBAP is a freely available Nextflow pipeline at github.com/hoelzer-lab/ribap and zenodo.org/doi/10.5281/zenodo.10890871.

Design of a novel multi-epitope vaccine candidate against Chlamydia trachomatis using structural and nonstructural proteins: an immunoinformatics study

Chlamydia trachomatis (C. trachomatis) is an obligate intracellular bacterium which causes eye and sexually transmitted infections. During pregnancy, the bacterium is associated with preterm complications, low weight of neonates, fetal demise and endometritis leading to infertility. The aim of our study was design of a multi-epitope vaccine (MEV) candidate against C. trachomatis. After protein sequence adoption from the NCBI, potential epitopes toxicity, antigenicity, allergenicity, MHC-I and MHC-II binding, cytotoxic T lymphocytes (CTLs), Helper T lymphocytes (HTLs) and interferon-γ (IFN-γ)- induction were predicted. The adopted epitopes were fused together using appropriate linkers. In the next step, the MEV structural mapping and characterization, three-dimensional (3D) structure homology modeling and refinement were also performed. The MEV candidate interaction with the toll-like receptor 4 (TLR4) was also docked. The immune responses simulation was assessed using the C-IMMSIM server. Molecular dynamic (MD) simulation verified the structural stability of the TLR4-MEV complex. The Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approach demonstrated the MEV high affinity of binding to the TLR4, MHC-I and MHC-II. The MEV construct was also stable and water soluble and had enough antigenicity and lacked allergenicity with stimulation of T cells and B cells and INF-γ release. The immune simulation confirmed acceptable responses of both the humoral and cellular arms. It is proposed that in vitro and in vivo studies are needed to evaluate the findings of this study.Communicated by Ramaswamy H. Sarma.

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.

Plasmid-mediated virulence in Chlamydia

Chlamydia trachomatis infection of ocular conjunctiva can lead to blindness, while infection of the female genital tract can lead to chronic pelvic pain, ectopic pregnancy, and/or infertility. Conjunctival and fallopian tube inflammation and the resulting disease sequelae are attributed to immune responses induced by chlamydial infection at these mucosal sites. The conserved chlamydial plasmid has been implicated in enhancing infection, via improved host cell entry and exit, and accelerating innate inflammatory responses that lead to tissue damage. The chlamydial plasmid encodes eight open reading frames, three of which have been associated with virulence: a secreted protein, Pgp3, and putative transcriptional regulators, Pgp4 and Pgp5. Although Pgp3 is an important plasmid-encoded virulence factor, recent studies suggest that chlamydial plasmid-mediated virulence extends beyond the expression of Pgp3. In this review, we discuss studies of genital, ocular, and gastrointestinal infection with C. trachomatis or C. muridarum that shed light on the role of the plasmid in disease development, and the potential for tissue and species-specific differences in plasmid-mediated pathogenesis. We also review evidence that plasmid-associated inflammation can be independent of bacterial burden. The functions of each of the plasmid-encoded proteins and potential molecular mechanisms for their role(s) in chlamydial virulence are discussed. Although the understanding of plasmid-associated virulence has expanded within the last decade, many questions related to how and to what extent the plasmid influences chlamydial infectivity and inflammation remain unknown, particularly with respect to human infections. Elucidating the answers to these questions could improve our understanding of how chlamydia augment infection and inflammation to cause disease.

Genomic analysis of 61 Chlamydia psittaci strains reveals extensive divergence associated with host preference

BACKGROUND

Chlamydia (C.) psittaci, the causative agent of avian chlamydiosis and human psittacosis, is a genetically heterogeneous species. Its broad host range includes parrots and many other birds, but occasionally also humans (via zoonotic transmission), ruminants, horses, swine and rodents. To assess whether there are genetic markers associated with host tropism we comparatively analyzed whole-genome sequences of 61 C. psittaci strains, 47 of which carrying a 7.6-kbp plasmid.

RESULTS

Following clean-up, reassembly and polishing of poorly assembled genomes from public databases, phylogenetic analyses using C. psittaci whole-genome sequence alignment revealed four major clades within this species. Clade 1 represents the most recent lineage comprising 40/61 strains and contains 9/10 of the psittacine strains, including type strain 6BC, and 10/13 of human isolates. Strains from different non-psittacine hosts clustered in Clades 2- 4. We found that clade membership correlates with typing schemes based on SNP types, ompA genotypes, multilocus sequence types as well as plasticity zone (PZ) structure and host preference. Genome analysis also revealed that i) sequence variation in the major outer membrane porin MOMP can result in 3D structural changes of immunogenic domains, ii) past host change of Clade 3 and 4 strains could be associated with loss of MAC/perforin in the PZ, rather than the large cytotoxin, iii) the distinct phylogeny of atypical strains (Clades 3 and 4) is also reflected in their repertoire of inclusion proteins (Inc family) and polymorphic membrane proteins (Pmps).

CONCLUSIONS

Our study identified a number of genomic features that can be correlated with the phylogeny and host preference of C. psittaci strains. Our data show that intra-species genomic divergence is associated with past host change and includes deletions in the plasticity zone, structural variations in immunogenic domains and distinct repertoires of virulence factors.

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.

Genetic and clinical characteristics of genital Chlamydia trachomatis infection in Guangzhou, China

BACKGROUND

Genital Chlamydia trachomatis (CT) is one of the most common agents of sexually transmitted infections and can cause severe disorders. This study aimed to analyse the genetic and clinical characteristics of genital CT infection among women in Guangzhou, China.

METHODS

From September 2020 to August 2021, a total of 8955 female patients were enrolled in this study. The presence of genital CT was detected by real-time PCR, and 273 positive samples were randomly selected for further genetic and clinical characteristics analysis.

RESULTS

The positive rate of genital CT infection was 7.5% (670/8955), with the highest rate in women aged 21-30 years. A total of 8 genotypes were identified: DH, J, K, and recombinant genotype Ba/D. The predominant genotype was J (n = 78, 28.6%), followed by E (n = 63, 23.1%), F (n = 48, 17.6%), and D (n = 38, 13.9%). Abnormal vaginal discharge (n = 165, 61.8%), cervical columnar epithelial ectopy (n = 124, 46.4%), vaginal itching (n = 77, 28.8%), and lower abdominal pain (n = 61, 22.8%) were the predominant symptoms. Additionally, genotype G infection exhibited a significantly higher rate of abnormal vaginal discharge (P = 0.03) and genotype D infection exhibited a higher white blood cell count (P = 0.01) than the other genotypes. Phylogenetic analysis revealed a total of 20 variants with 25 mutation positions and the H2 variant in four patients was first discovered in our study.

CONCLUSIONS

Genotypes J, E, F, and D were the major genotypes of genital CT in Guangzhou, and they manifested as abnormal vaginal discharge, cervical columnar epithelial ectopy, vaginal itching, and lower abdominal pain. The present study provides guidance for future integrated interventions to reduce the burden of genital CT infection and accelerate the development of vaccines.

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.

SARS-CoV-2 outbreak in a tri-national urban area is dominated by a B.1 lineage variant linked to a mass gathering event

The first case of SARS-CoV-2 in Basel, Switzerland was detected on February 26th 2020. We present a phylogenetic study to explore viral introduction and evolution during the exponential early phase of the local COVID-19 outbreak from February 26th until March 23rd. We sequenced SARS-CoV-2 naso-oropharyngeal swabs from 746 positive tests that were performed at the University Hospital Basel during the study period. We successfully generated 468 high quality genomes from unique patients and called variants with our COVID-19 Pipeline (COVGAP), and analysed viral genetic diversity using PANGOLIN taxonomic lineages. To identify introduction and dissemination events we incorporated global SARS-CoV-2 genomes and inferred a time-calibrated phylogeny. Epidemiological data from patient questionnaires was used to facilitate the interpretation of phylogenetic observations. The early outbreak in Basel was dominated by lineage B.1 (83·6%), detected first on March 2nd, although the first sample identified belonged to B.1.1. Within B.1, 68·2% of our samples fall within a clade defined by the SNP C15324T ('Basel cluster'), including 157 identical sequences at the root of the 'Basel cluster', some of which we can specifically trace to regional spreading events. We infer the origin of B.1-C15324T to mid-February in our tri-national region. The other genomes map broadly over the global phylogenetic tree, showing several introduction events from and/or dissemination to other regions of the world via travellers. Family transmissions can also be traced in our data. A single lineage variant dominated the outbreak in the Basel area while other lineages, such as the first (B.1.1), did not propagate. A mass gathering event was the predominant initial source of cases, with travel returners and family transmissions to a lesser extent. We highlight the importance of adding specific questions to epidemiological questionnaires, to obtain data on attendance of large gatherings and their locations, as well as travel history, to effectively identify routes of transmissions in up-coming outbreaks. This phylogenetic analysis in concert with epidemiological and contact tracing data, allows connection and interpretation of events, and can inform public health interventions. Trial Registration: ClinicalTrials.gov NCT04351503.

Coevolving Plasmids Drive Gene Flow and Genome Plasticity in Host-Associated Intracellular Bacteria

Plasmids are important in microbial evolution and adaptation to new environments. Yet, carrying a plasmid can be costly, and long-term association of plasmids with their hosts is poorly understood. Here, we provide evidence that the Chlamydiae, a phylum of strictly host-associated intracellular bacteria, have coevolved with their plasmids since their last common ancestor. Current chlamydial plasmids are amalgamations of at least one ancestral plasmid and a bacteriophage. We show that the majority of plasmid genes are also found on chromosomes of extant chlamydiae. The most conserved plasmid gene families are predominantly vertically inherited, while accessory plasmid gene families show significantly increased mobility. We reconstructed the evolutionary history of plasmid gene content of an entire bacterial phylum over a period of around one billion years. Frequent horizontal gene transfer and chromosomal integration events illustrate the pronounced impact of coevolution with these extrachromosomal elements on bacterial genome dynamics in host-dependent microbes.

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.

Diversity in Chlamydial plasmids

Background

Evolutionary studies have been conducted that have investigated the chromosomal variance in the genus of Chlamydia. However, no all-encompassing genus-wide comparison has been performed on the plasmid. Therefore, there is a gap in the current knowledge on Chlamydia plasmid diversity.

Aims

This project is aimed to investigate and establish the nature and extent of diversity across the entire genus of Chlamydia, by comparing the sequences of all currently available plasmid carrying strains.

Methods

The PUBMED database was used to identify plasmid sequences from all available strains that met the set quality criteria for their inclusion in the study. Alignments were performed on the 51 strains that fulfilled the criteria using MEGA X software. Following that Maximum Likelihood estimation was used to construct 11 phylogenetic trees of the whole plasmid sequence, the individual 8 coding sequences, the iteron and a chromosomal gene ompA as a comparator.

Results

The genus-wide plasmid phylogeny produced three distinct lineages labelled as alpha, beta and gamma. Nineteen genotypes were found in the initial whole plasmid analysis. Their distribution was allocated as six C. pecorum, two C. pneumoniae, one C. gallinacea, one C. avium, one C. caviae, one C. felis, two C. psittaci, one C. trachomatis, one C. muridarum, and two C. suis. The chromosomal comparative gene ompA supported this distribution, with the same number of primary clades with the same species distribution. However, ompA sequence comparison resulted in fewer genotypes due to a reduced amount of available sequences (33 out of 51). All results were statistically significant.

Conclusion

The results of this study indicate that the common bacterial ancestor of all the species had a plasmid, which has diverged over time. Moreover, it suggests that there is a strong evolutionary selection towards these species retaining their plasmids due to its high level of conservation across the genus, with the notable exception of C. pneumoniae. Furthermore, the evolutionary analysis showed that the plasmid and the chromosome have co-evolved.

The Nature and Extent of Plasmid Variation in Chlamydia trachomatis

Chlamydia trachomatis is an obligate intracellular pathogen of humans, causing both the sexually transmitted infection, chlamydia, and the most common cause of infectious blindness, trachoma. The majority of sequenced C. trachomatis clinical isolates carry a 7.5-Kb plasmid, and it is becoming increasingly evident that this is a key determinant of pathogenicity. The discovery of the Swedish New Variant and the more recent Finnish variant highlight the importance of understanding the natural extent of variation in the plasmid. In this study we analysed 524 plasmid sequences from publicly available whole-genome sequence data. Single nucleotide polymorphisms (SNP) in each of the eight coding sequences (CDS) were identified and analysed. There were 224 base positions out of a total 7550 bp that carried a SNP, which equates to a SNP rate of 2.97%, nearly three times what was previously calculated. After normalising for CDS size, CDS8 had the highest SNP rate at 3.97% (i.e., number of SNPs per total number of nucleotides), whilst CDS6 had the lowest at 1.94%. CDS5 had the highest total number of SNPs across the 524 sequences analysed (2267 SNPs), whereas CDS6 had the least SNPs with only 85 SNPs. Calculation of the genetic distances identified CDS6 as the least variable gene at the nucleotide level (d = 0.001), and CDS5 as the most variable (d = 0.007); however, at the amino acid level CDS2 was the least variable (d = 0.001), whilst CDS5 remained the most variable (d = 0.013). This study describes the largest in-depth analysis of the C. trachomatis plasmid to date, through the analysis of plasmid sequence data mined from whole genome sequences spanning 50 years and from a worldwide distribution, providing insights into the nature and extent of existing variation within the plasmid as well as guidance for the design of future diagnostic assays. This is crucial at a time when single-target diagnostic assays are failing to detect natural mutants, putting those infected at risk of a serious long-term and life-changing illness.

Chlamydiae from Down Under: The Curious Cases of Chlamydial Infections in Australia

In Australia, the most researched and perhaps the most successful chlamydial species are the human pathogen Chlamydia trachomatis, animal pathogens Chlamydia pecorum and Chlamydia psittaci. C. trachomatis remains the leading cause of sexually transmitted infections in Australians and trachoma in Australian Indigenous populations. C. pecorum is globally recognised as the infamous koala and widespread livestock pathogen, whilst the avian C. psittaci is emerging as a horse pathogen posing zoonotic risks to humans. Certainly not innocuous, the human infections with Chlamydia pneumoniae seem to be less prevalent that other human chlamydial pathogens (namely C. trachomatis). Interestingly, the complete host range for C. pecorum and C. psittaci remains unknown, and infections by other chlamydial organisms in Australian domesticated and wildlife animals are understudied. Considering that chlamydial organisms can be encountered by either host at the human/animal interface, I review the most recent findings of chlamydial organisms infecting Australians, domesticated animals and native wildlife. Furthermore, I also provide commentary from leading Australian Chlamydia experts on challenges and future directions in the Chlamydia research field.

Whole-genome sequencing of ocular Chlamydia trachomatis isolates from Gadarif State, Sudan

BACKGROUND

Trachoma, caused by ocular Chlamydia trachomatis, is the leading infectious cause of blindness worldwide. Sudan first reported trachoma in the 1930s and has since been consistently endemic. Ocular C. trachomatis previously isolated from trachoma patients in Sudan in 1963 was antigenically identical to an isolate from Saudi Arabia (A/SA1). No contemporary ocular C. trachomatis whole genome sequences have been reported from Sudan.

METHODS

This study sequenced twenty ocular C. trachomatis isolates to improve understanding of pathogen diversity in North-East Africa and examine for genomic variation specific to Sudan, possibly related to the persistence of trachoma in surveyed communities. High quality, whole genome sequences were obtained from 12/20 isolates.

RESULTS

All isolates were serovar A and had tarP and trpA sequences typical of classical, ocular C. trachomatis isolates. The Sudanese isolates formed a closely related subclade within the T2-trachoma clade of C. trachomatis phylogeny distinct from geographically disparate ocular isolates, with little intra-population diversity. We found 333 SNPs that were conserved in Sudanese ocular isolates but rare compared to other ocular C. trachomatis populations, which were focused in two genomic loci (CTA0172-CTA0173 and CTA0482).

CONCLUSIONS

Limited intra-population diversity and geographical clustering of ocular C. trachomatis suggests minimal transmission between and slow diversification within trachoma-endemic communities. However, diversity may have been higher pre-treatment in these communities. Over-representation of Sudan-specific SNPs in three genes suggests they may have an impact on C. trachomatis growth and transmission in this population.

Multilocus Sequence Typing (MLST) of Chlamydiales

Developed two decades ago as a molecular method to provide definite characterization of a bacterial isolate, Multilocus Sequence Typing (MLST) is today globally adopted as a universal fine-detailed molecular typing tool and has been applied to numerous pathogenic and nonpathogenic bacterial as well eukaryotic organisms. MLST utilizes DNA sequence of several conserved housekeeping (HK) genes which are assigned an allelic number, which then collectively constitute an allelic profile or sequence type (ST), a "molecular barcode" of the interrogated bacterial strain or a eukaryotic organism. Here, we describe the principles and molecular approaches for generating MLST data for an analysis of a bacteria in the order Chlamydiales, using a Chlamydia pecorum-specific MLST scheme as an example.

Genetic Transformation of a C. trachomatis Ocular Isolate With the Functional Tryptophan Synthase Operon Confers an Indole-Rescuable Phenotype

Chlamydia trachomatis is the leading cause of preventable blindness and the most common bacterial sexually transmitted infection. Different strains are associated with ocular or urogenital infections, and a proposed mechanism that may explain this tissue tropism is the active tryptophan biosynthesis pathway encoded by the genomic trpRBA operon in urogenital strains. Here we describe genetic complementation studies that are essential to confirm the role of tryptophan synthase in the context of an ocular C. trachomatis genomic background. Ocular strain A2497 was transformed with the (urogenital) pSW2::GFP shuttle vector showing that there is no strain tropism barrier to this plasmid vector; moreover, transformation had no detrimental effect on the growth kinetics of A2497, which is important given the low transformation efficiency of C. trachomatis. A derivative of the pSW2::GFP vector was used to deliver the active tryptophan biosynthesis genes from a urogenital strain of C. trachomatis (Soton D1) to A2497 with the aim of complementing the truncated trpA gene common to most ocular strains. After confirmation of intact TrpA protein expression in the transformed A2497, the resulting transformants were cultivated in tryptophan-depleted medium with and without indole or tryptophan, showing that complementation of the truncated trpA gene by the intact and functional urogenital trpRBA operon was sufficient to bestow an indole rescuable phenotype upon A2497. This study proves that pSW2::GFP derived vectors do not conform to the cross-strain transformation barrier reported for other chlamydia shuttle vectors, suggesting these as a universal vector for transformation of all C. trachomatis strains. This vector promiscuity enabled us to test the indole rescue hypothesis by transforming ocular strain A2497 with the functional urogenital trpRBA operon, which complemented the non-functional tryptophan synthase. These data confirm that the trpRBA operon is necessary and sufficient for chlamydia to survive in tryptophan-limited environments such as the female urogenital tract.

Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications

The PubMLST.org website hosts a collection of open-access, curated databases that integrate population sequence data with provenance and phenotype information for over 100 different microbial species and genera.  Although the PubMLST website was conceived as part of the development of the first multi-locus sequence typing (MLST) scheme in 1998 the software it uses, the Bacterial Isolate Genome Sequence database (BIGSdb, published in 2010), enables PubMLST to include all levels of sequence data, from single gene sequences up to and including complete, finished genomes.  Here we describe developments in the BIGSdb software made from publication to June 2018 and show how the platform realises microbial population genomics for a wide range of applications.  The system is based on the gene-by-gene analysis of microbial genomes, with each deposited sequence annotated and curated to identify the genes present and systematically catalogue their variation.  Originally intended as a means of characterising isolates with typing schemes, the synthesis of sequences and records of genetic variation with provenance and phenotype data permits highly scalable (whole genome sequence data for tens of thousands of isolates) means of addressing a wide range of functional questions, including: the prediction of antimicrobial resistance; likely cross-reactivity with vaccine antigens; and the functional activities of different variants that lead to key phenotypes.  There are no limitations to the number of sequences, genetic loci, allelic variants or schemes (combinations of loci) that can be included, enabling each database to represent an expanding catalogue of the genetic variation of the population in question.  In addition to providing web-accessible analyses and links to third-party analysis and visualisation tools, the BIGSdb software includes a RESTful application programming interface (API) that enables access to all the underlying data for third-party applications and data analysis pipelines.

Unveiling the Multilocus Sequence Typing (MLST) Schemes and Core Genome Phylogenies for Genotyping Chlamydia trachomatis

Multilocus sequence typing (MLST) has become a useful tool for studying the genetic diversity of important public health pathogens, such as Chlamydia trachomatis (Ct). Four MLST schemes have been proposed for Ct (data available from Chlamydiales MLST databases). However, the lack of a sole standardized scheme represents the greatest limitation regarding typing this species. This study was thus aimed at evaluating the usefulness of the four MLST schemes available for Ct, describing each molecular marker's pattern and its contribution toward a description of intra-specific genetic diversity and population structure. The markers for each scheme, showed a variable power of dicrimination, exhibiting in some cases over estimation in the determination of Sequence Types (STs). However, individual analysis of each locus's typing efficiency and discrimination power led to identifying 8 markers as having a suitable pattern for intra-specific typing. analyzing the 8 candidate markers gave a combination of 3 of these loci as an optimal scheme for identifying a large amount of STs, maximizing discrimination power whilst maintaining suitable typing efficiency. One scheme was compared against core genome phylogenies, finding a higher typing resolution through the last approach. These results confirm once again that although complete genome data, in particular from core genome MLST (cgMLST) allow a high resolution clustering for Ct isolates. There are combinations of molecular markers that could generate equivalent results, with the advantage of representing an easy implementation strategy and lower costs leading to contribute to the monitoring and molecular epidemiology of Ct.