Whole-genome analysis of diverse Chlamydia trachomatis strains identifies phylogenetic relationships masked by current clinical typing

Unravelling Chlamydia trachomatis diversity in Amhara, Ethiopia: MLVA-ompA sequencing as a molecular typing tool for trachoma

Trachoma is the leading infectious cause of blindness worldwide and is now largely confined to around 40 low- and middle-income countries. It is caused by Chlamydia trachomatis (Ct), a contagious intracellular bacterium. The World Health Organization recommends mass drug administration (MDA) with azithromycin for treatment and control of ocular Ct infections, alongside improving facial cleanliness and environmental conditions to reduce transmission. To understand the molecular epidemiology of trachoma, especially in the context of MDA and transmission dynamics, the identification of Ct genotypes could be useful. While many studies have used the Ct major outer membrane protein gene (ompA) for genotyping, it has limitations. Our study applies a typing system novel to trachoma, Multiple Loci Variable Number Tandem Repeat Analysis combined with ompA (MLVA-ompA). Ocular swabs were collected post-MDA from four trachoma-endemic zones in Ethiopia between 2011-2017. DNA from 300 children with high Ct polymerase chain reaction (PCR) loads was typed using MLVA-ompA, utilizing 3 variable number tandem repeat (VNTR) loci within the Ct genome. Results show that MLVA-ompA exhibited high discriminatory power (0.981) surpassing the recommended threshold for epidemiological studies. We identified 87 MLVA-ompA variants across 26 districts. No significant associations were found between variants and clinical signs or chlamydial load. Notably, overall Ct diversity significantly decreased after additional MDA rounds, with a higher proportion of serovar A post-MDA. Despite challenges in sequencing one VNTR locus (CT1299), MLVA-ompA demonstrated cost-effectiveness and efficiency relative to whole genome sequencing, providing valuable information for trachoma control programs on local epidemiology. The findings suggest the potential of MLVA-ompA as a reliable tool for typing ocular Ct and understanding transmission dynamics, aiding in the development of targeted interventions for trachoma control.

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.

A retrospective evaluation of the Euroarray STI-11 multiplex system for the detection of eight STI causing agents

With an incidence of more than > 1,000,000/day, sexually transmitted diseases remain a major challenge for health care systems worldwide. To reduce disease burden, complications, and spread, rapid diagnosis permitting early therapy is pivotal. The range of pathogens is wide and co-infections are common. This complicates pre-analytics, which are based on different laboratory techniques with potentially long turnaround times, e.g., cultivation and multistep serologies. Multiplex PCR provides the opportunity to overcome these limitations. In this study, we evaluated a novel assay, the Euroarray STI-11 microarray (EA; Euroimmun Medizinische Labordiagnostika), for the detection of eight obligate or facultative pathogens. Three-hundred-thirteen clinical specimens, which had been tested and pre-characterized for STI causing agents as part of routine diagnostics, were used as cases and controls in this retrospective study. The EA detected 34/44 Chlamydia trachomatis, 48/50 HSV-1, 50/50 HSV-2, 48/48 Mycoplasma hominis, 45/47 Neisseria gonorrhoeae, 9/11 Treponema pallidum, 46/46 Ureaplasma parvum, and 49/49 Ureaplasma urealyticum infections, respectively. 293 samples were EA positive, with polymicrobial infections (positive for two to six microbial or viral agents) detected in 130/293 cases. Specificities were 100% in the respective control groups (n = 18-48 depending on targeted pathogen) except for N. gonorrhoeae (25/26) and U. urealyticum (44/45). The broad spectrum of obligate and facultative pathogens targeted by the EA makes it a valuable tool in the setting of STI diagnostics and surveillance. The test has the potential to diagnose diseases neglected or overlooked in routine clinical practice. Besides a low sensitivity for C. trachomatis, the EA demonstrated high performance for all analyzed parameters. Further studies are warranted in order to capture a larger variety of the tested pathogens.

Whole-genome sequencing of Chlamydia psittaci from Australasian avian hosts: A genomics approach to a pathogen that still ruffles feathers

Chlamydia psittaci is a globally distributed veterinary pathogen with zoonotic potential. Although C. psittaci infections have been reported in various hosts, isolation and culture of Chlamydia is challenging, hampering efforts to produce contemporary global C. psittaci genomes. This is particularly evident in the lack of avian C. psittaci genomes from Australia and New Zealand. In this study, we used culture-independent probe-based whole-genome sequencing to expand the global C. psittaci genome catalogue. Here, we provide new C. psittaci genomes from two pigeons, six psittacines, and novel hosts such as the Australian bustard (Ardeotis australis) and sooty shearwater (Ardenna grisea) from Australia and New Zealand. We also evaluated C. psittaci genetic diversity using multilocus sequence typing (MLST) and major outer membrane protein (ompA) genotyping on additional C. psittaci-positive samples from various captive avian hosts and field isolates from Australasia. We showed that the first C. psittaci genomes sequenced from New Zealand parrots and pigeons belong to the clonal sequence type (ST)24 and diverse 'pigeon-type' ST27 clade, respectively. Australian parrot-derived strains also clustered in the ST24 group, whereas the novel ST332 strain from the Australian bustard clustered in a genetically diverse clade of strains from a fulmar, parrot, and livestock. MLST and ompA genotyping revealed ST24/ompA genotype A in wild and captive parrots and a sooty shearwater, whilst 'pigeon-types' (ST27/35 and ompA genotypes B/E) were found in pigeons and other atypical hosts, such as captive parrots, a little blue penguin/Kororā (Eudyptula minor) and a zebra finch (Taeniopygia guttata castanotis) from Australia and New Zealand. This study provides new insights into the global phylogenomic diversity of C. psittaci and further demonstrates the multi-host generalist capacity of this pathogen.

Prevalence of sexually transmissible infections in adolescents treated in a family planning outpatient clinic for adolescents in the western Amazon

Sexually transmitted infections (STIs) are among the most common public health problems worldwide, especially among adolescents and young adults, who account for almost 50% of all STI patients. Studies on the subject in the western Amazon are limited. This study aimed to evaluate the prevalence of STIs (chlamydia, gonorrhea, trichomoniasis, herpes simplex virus, syphilis, human immunodeficiency virus [HIV], hepatitis B, and hepatitis C) in adolescents treated at a family planning outpatient clinic in the western Amazon: Porto Velho, Rondônia, Brazil. A total of 196 adolescents were enrolled. During the gynecological examination, endocervical samples were collected to test for four STIs (chlamydia, gonorrhea, trichomoniasis, and herpes simplex virus), and blood samples were collected for the detection of HIV, syphilis, and hepatitis B and C. The mean age was 17.3 ± 1.5 years, the age at sexarche was 14.4 ± 1.6 years, and 54.6% of participants had their first sexual intercourse at 14 years or younger. Only 1.0% of the adolescents used condoms in all sexual relations, and 19.9% had casual partner(s) in the last year. In the evaluation of prevalence, we found that 32% of the adolescents had at least one STI, with the most prevalent being chlamydia (23%), followed by trichomoniasis (5.6%), herpes simplex (4.6%), and gonorrhea (3.1%). No positive cases of hepatitis B, hepatitis C, or HIV were detected, but 1% of the adolescents tested positive for syphilis. These indicators will support more effective health care strategies aimed at improving the quality of life of populations in this region of the western Amazon. In conclusion, our findings demonstrated high rates of STIs in the studied patients, reinforcing the need to expand epidemiological studies to implement more appropriate public policies and intervention strategies to prevent STIs in adolescents and other vulnerable populations in the western Amazon.

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 Characterization of Non-Lymphogranuloma venereum Chlamydia trachomatis Indicates Distinct Infection Transmission Networks in Spain

Chlamydia trachomatis infection is an important public health problem. Our objective was to assess the dynamics of the transmission of this infection, analysing the distribution of circulating ompA genotypes and multilocus sequence types of C. trachomatis in Spain as a function of clinical and epidemiological variables. During 2018 and 2019, we genetically characterized C. trachomatis in tertiary hospitals in six areas in Spain (Asturias, Barcelona, Gipuzkoa, Mallorca, Seville and Zaragoza), with a catchment population of 3.050 million people. Genotypes and sequence types were obtained using polymerase chain reaction techniques that amplify a fragment of the ompA gene, and five highly variable genes (hctB, CT058, CT144, CT172 and pbpB), respectively. Amplicons were sequenced and phylogenetic analysis was conducted. We obtained genotypes in 636/698 cases (91.1%). Overall and by area, genotype E was the most common (35%). Stratifying by sex, genotypes D and G were more common among men, and genotypes F and I among women (p < 0.05). Genotypes D, G and J were more common in men who have sex with men (MSM) than in men who have sex with women (MSW), in whom the most common genotypes were E and F. The diversity index was higher in sequence typing (0.981) than in genotyping (0.791), and the most common sequence types were ST52 and ST108 in MSM, and ST30, ST148, ST276 and ST327 in MSW. Differences in genotype distribution between geographical areas were attributable to differences in population characteristics. The transmission dynamics varied with sexual behaviour: the predominant genotypes and most frequent sequence types found in MSM were different to those detected in MSW and women.

Specific high-resolution scheme to improve understanding of the spatio-temporal dispersion of lymphogranuloma venereum epidemic

Introduction

Lymphogranuloma venereum (LGV) is already endemic in vulnerable populations in several European countries; however, molecular epidemiology data with improved accuracy are necessary to better understand LGV epidemic in these countries. Current strategies to study the molecular epidemiology of LGV cases involve schemes based on a few genetic fragments of Chlamydia trachomatis, which have demonstrated limited discriminatory power for LGV. Therefore, this study aimed to propose a new combination of molecular markers based on the most variable genes of L-genotype genomes to improve the characterization of the current LGV epidemic in Madrid, Spain.

Methods

Four genes were selected according to their diversity index (CTLon_0054, CTLon_0087, CTLon_0243 and CTLon_0301) for use in combination with ompA. In silico and experimental studies were performed to compare the previously described multilocus sequence typing (MLST) schemes with our proposal. Moreover, the proposed scheme was applied (n = 68) to analyze the spatio-temporal spread of the LGV cases.

Results

Our proposal demonstrated higher diversity allowing the identification of three main groups compared to the previously published MLST based on hypervariable genes wherein only a single sequence type was identified. The temporal analysis showed that the major cluster was progressively diversifying, revealing a very active transmission chain. Furthermore, an L2b genome identical to that of the origin of the epidemic was detected, suggesting reintroductions or a low screening rate in vulnerable populations. The spatial distribution suggests that the selection and spread of new variants occurs from the central district to the peripheral regions.

Discussion

The scheme proposed in this study has proven to be useful for appropriate discrimination of LGV strains. This study, to our knowledge for the first time, demonstrates a spatio-temporal spread that increases our understanding and identifies areas with special susceptibility for maintenance of the endemic situation of LGV.

Genomic Analysis of MSM Rectal Chlamydia trachomatis Isolates Identifies Predicted Tissue-Tropic Lineages Generated by Intraspecies Lateral Gene Transfer-Mediated Evolution

Chlamydia trachomatis is an obligate intracellular bacterium that causes serious diseases in humans. Rectal infection and disease caused by this pathogen are important yet understudied aspects of C. trachomatis natural history. The University of Washington Chlamydia Repository has a large collection of male-rectal-sourced strains (MSM rectal strains) isolated in Seattle, USA and Lima, Peru. Initial characterization of strains collected over 30 years in both Seattle and Lima led to an association of serovars G and J with male rectal infections. Serovar D, E, and F strains were also collected from MSM patients. Genome sequence analysis of a subset of MSM rectal strains identified a clade of serovar G and J strains that had high overall genomic identity. A genome-wide association study was then used to identify genomic loci that were correlated with tissue tropism in a collection of serovar-matched male rectal and female cervical strains. The polymorphic membrane protein PmpE had the strongest correlation, and amino acid sequence alignments identified a set of PmpE variable regions (VRs) that were correlated with host or tissue tropism. Examination of the positions of VRs by the protein structure-predicting Alphafold2 algorithm demonstrated that the VRs were often present in predicted surface-exposed loops in both PmpE and PmpH protein structure. Collectively, these studies identify possible tropism-predictive loci for MSM rectal C. trachomatis infections and identify predicted surface-exposed variable regions of Pmp proteins that may function in MSM rectal versus cervical tropism differences.

Identification and characterization of mixed infections of Chlamydia trachomatis via high-throughput sequencing

Precise genotyping is necessary to understand epidemiology and clinical manifestations of Chlamydia trachomatis infection with different genotypes. Next-generation high-throughput sequencing (NGHTS) has opened new frontiers in microbial genotyping, but has been clinically characterized in only a few settings. This study aimed to determine C. trachomatis genotypes in particular mixed-genotype infections and their association with clinical manifestations and to characterize the sensitivity and accuracy of NGHTS. Cervical specimens were collected from 8,087 subjects from physical examination center (PEC), assisted reproductive technology center (ART) and gynecology clinics (GC) of Chenzhou Hospital of China. The overall prevalence of C. trachomatis was 3.8% (311/8087) whereas a prevalence of 2.8, 3.7 and 4.8% was found in PEC, ART and GC, respectively. The most frequent three C. trachomatis genotypes were E (27.4%, 83/303), F (21.5%, 65/303) and J (18.2%, 55/303). Moreover, NGHTS identified 20 (6.6%, 20/303) mixed-genotype infections of C. trachomatis. Genotype G was more often observed in the subjects with pelvic inflammatory disease than genotype E (adjusted OR = 3.61, 95%CI, 1.02–12.8, p = 0.046). Mixed-genotype infection was associated with severe vaginal cleanliness (degree IV) with an adjusted OR of 5.17 (95%CI 1.03–25.9, p = 0.046) whereas mixed-genotype infection with large proportion of minor genotypes was associated with cervical squamous intraepithelial lesion (SIL) with an adjusted OR of 5.51 (95%CI 1.17–26.01, p = 0.031). Our results indicated that NGHTS is a feasible tool to identity C. trachomatis mixed-genotype infections, which may be associated with worse vaginal cleanliness and cervical SIL.

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.

Chlamydia trachomatis Requires Functional Host-Cell Mitochondria and NADPH Oxidase 4/p38MAPK Signaling for Growth in Normoxia

Chlamydia trachomatis (Ct) is an intracellular energy-parasitic bacterium that requires ATP derived from infected cells for its growth. Meanwhile, depending on the O₂ concentration, the host cells change their mode of ATP production between oxidative phosphorylation in mitochondria (Mt) and glycolysis; this change depends on signaling via reactive oxygen species (ROS) produced by NADPH oxidases (NOXs) as well as Mt. It has been proposed that Ct correspondingly switches its source of acquisition of ATP between host-cell Mt and glycolysis, but this has not been verified experimentally. In the present study, we assessed the roles of host-cell NOXs and Mt in the intracellular growth of CtL2 (L2 434/Bu) under normoxia (21% O₂) and hypoxia (2% O₂) by using several inhibitors of NOXs (or the downstream molecule) and Mt-dysfunctional (Mt^d) HEp-2 cells. Under normoxia, diphenyleneiodonium, an inhibitor of ROS diffusion, abolished the growth of CtL2 and other Chlamydiae (CtD and C. pneumoniae). Both ML171 (a pan-NOX inhibitor) and GLX351322 (a NOX4-specific inhibitor) impaired the growth of CtL2 under normoxia, but not hypoxia. NOX4-knockdown cells diminished the bacterial growth. SB203580, an inhibitor of the NOX4-downstream molecule p38MAPK, also inhibited the growth of CtL2 under normoxia but not hypoxia. Furthermore, CtL2 failed to grow in Mt^d cells under normoxia, but no effect was observed under hypoxia. We conclude that under normoxia, Ct requires functional Mt in its host cells as an ATP source, and that this process requires NOX4/p38MAPK signaling in the host cells. In contrast to hypoxia, crosstalk between NOX4 and Mt via p38MAPK may be crucial for the growth of Ct under normoxia.

A Review of Next Generation Sequencing Methods and its Applications in Laboratory Diagnosis

Next-generation sequencing (NGS) is a new technology used to detect the sequence of DNA and RNA and to detect mutations or variations of significance. NGS generates large quantities of sequence data within a short time duration. The various types of sequencing includes Sanger Sequencing, Pyrosequencing, Sequencing by Synthesis (Illumina), Ligation (SoLID), Single molecule Fluorescent Sequencing (Helicos), Single molecule Real time Sequencing (Pacbio), Semiconductor sequencing (Ion torrent technology), Nanopore sequencing and fourth generation sequencing. These methods of sequencing have been modified and improved over the years such that it has become cost effective and accessible to diagnostic laboratories. Management of Outbreaks, rapid identification of bacteria, molecular case finding, taxonomy, detection of the zoonotic agents and guiding prevention strategies in HIV outbreaks are just a few of the many applications of Next Generation sequencing in clinical microbiology.

Trachoma

Trachoma is a neglected tropical disease caused by infection with conjunctival strains of Chlamydia trachomatis. It can result in blindness. Pathophysiologically, trachoma is a disease complex composed of two linked chronic processes: a recurrent, generally subclinical infectious-inflammatory disease that mostly affects children, and a non-communicable, cicatricial and, owing to trichiasis, eventually blinding disease that supervenes in some individuals later in life. At least 150 infection episodes over an individual's lifetime are needed to precipitate trichiasis; thus, opportunity exists for a just global health system to intervene to prevent trachomatous blindness. Trachoma is found at highest prevalence in the poorest communities of low-income countries, particularly in sub-Saharan Africa; in June 2021, 1.8 million people worldwide were going blind from the disease. Blindness attributable to trachoma can appear in communities many years after conjunctival C. trachomatis transmission has waned or ceased; therefore, the two linked disease processes require distinct clinical and public health responses. Surgery is offered to individuals with trichiasis and antibiotic mass drug administration and interventions to stimulate facial cleanliness and environmental improvement are designed to reduce infection prevalence and transmission. Together, these interventions comprise the SAFE strategy, which is achieving considerable success. Although much work remains, a continuing public health problem from trachoma in the year 2030 will be difficult for the world to excuse.

Whole-genome sequencing of Chlamydia trachomatis isolates from persistently infected patients

BACKGROUND

Current understanding of the causes of treatment failure in Chlamydia trachomatis is poor and antimicrobial susceptibility data are lacking. We used genome sequencing to seek evidence of antimicrobial resistance in isolates sourced from patients who were persistently infected.

METHODS

Genomic DNA was extracted from C. trachomatis isolates cultured in McCoy cell monolayers. Sequencing libraries were prepared using the SureSelect^XT Illumina paired-end protocol. Paired reads were mapped against a reference genome and single nucleotide variants (SNVs) were identified.

RESULTS

Seven isolates from persistently infected patients and five isolates from successfully treated patients were sequenced. No previously reported SNVs associated with antimicrobial resistance were found. A unique SNV was identified in the gyrA gene of one treatment failure isolate but was located outside of the quinolone resistance determining region; this SNV has been previously reported in other members of the Chlamydiaceae family.

CONCLUSION

No genomic evidence was found to explain the differences in clinical outcome for our two groups of patients. A mutation unrelated to antimicrobial susceptibility was found in an isolate from a persistently infected patient. The cause of these persistent infections with C. trachomatis remains unclear.

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.

Emergence of Novel Chlamydia trachomatis Sequence Types among Chlamydia Patients in the Republic of Belarus

Chlamydia trachomatis (CT) is a major cause of sexually transmitted diseases worldwide. The multilocus sequence typing (MLST) of clinical samples from random heterosexual chlamydia patients who were either asymptomatic or reported clinical manifestations of genital chlamydiosis (n = 63) in each of the seven major regions of the Republic of Belarus in 2017–2018 revealed 12 different CT sequence types (STs). We found seven known STs, ST4, ST6, ST9, ST13, ST38, ST95 and ST110, and five novel variants, namely ST271–ST275, which have not been detected elsewhere thus far. The ST4 variant was predominant (27/63, 42.9%) and detected in six out of seven regions. The two most common STs, ST9 and ST13, were regularly seen in four out of seven regions. In contrast, the remaining STs, ST6, ST38, ST95, ST110, and novel STs271-275, surfaced randomly in different parts of the country. The emergence of novel STs was registered in two regions, namely Minsk (ST271 and ST275) and Brest (ST271, ST272, ST273, and ST274). All the STs of detected CT strains were clustered into two Groups, I and III, which are characteristic of CT urogenital strains. No STs typical for Group II, specific to the LGV strains, were revealed. Our study contributes to better understanding the genetic diversity and molecular evolution of CT, one of the most important pathogens in public health worldwide.

Problems with Understanding Chlamydia trachomatis Immunology

The pay off for understanding Chlamydia. trachomatis (CT) immunology is the development of a vaccine. Two lines of research have contributed to our current understanding of CT immunology. The first is the Grayston model of type specific immunity and genus specific pathology which was elaborated by Caldwell and Morrison as the major outer membrane protein (MOMP) and heat shock protein 60 (HSP60) paradigm. The second is the murine model of C. muridarum (CM) infection which established the essential role of major histocompatibility complex class II and CD4 T cells in immunity. However neither approach has yielded a vaccine. I review these two lines of research and conclude with six problem areas in human CT immunology whose resolution may result in a vaccine.(122 words)

A broad-spectrum cloning vector that exists as both an integrated element and a free plasmid in Chlamydia trachomatis

Plasmid transformation of chlamydiae has created new opportunities to investigate host-microbe interactions during chlamydial infections; however, there are still limitations. Plasmid transformation requires a replicon derived from the native Chlamydia plasmid, and these transformations are species-specific. We explored the utility of a broad host-range plasmid, pBBR1MCS-4, to transform chlamydiae, with a goal of simplifying the transformation process. The plasmid was modified to contain chromosomal DNA from C. trachomatis to facilitate homologous recombination. Sequences flanking incA were cloned into the pBBR1MCS-4 vector along with the GFP:CAT cassette from the pSW2-GFP chlamydial shuttle vector. The final plasmid construct, pBVR2, was successfully transformed into C. trachomatis strain L2-434. Chlamydial transformants were analyzed by immunofluorescence microscopy and positive clones were sequentially purified using limiting dilution. PCR and PacBio-based whole genome sequencing were used to determine if the plasmid was maintained within the chromosome or as an episome. PacBio sequencing of the cloned transformants revealed allelic exchange events between the chromosome and plasmid pBVR2 that replaced chromosomal incA with the plasmid GFP:CAT cassette. The data also showed evidence of full integration of the plasmid into the bacterial chromosome. While some plasmids were fully integrated, some were maintained as episomes and could be purified and retransformed into E. coli. Thus, the plasmid can be successfully transformed into chlamydia without a chlamydial origin of replication and can exist in multiple states within a transformed population.

Completing the Genome Sequence of Chlamydia pecorum Strains MC/MarsBar and DBDeUG: New Insights into This Enigmatic Koala (Phascolarctos cinereus) Pathogen

Chlamydia pecorum, an obligate intracellular pathogen, causes significant morbidity and mortality in livestock and the koala (Phascolarctos cinereus). A variety of C. pecorum gene-centric molecular studies have revealed important observations about infection dynamics and genetic diversity in both koala and livestock hosts. In contrast to a variety of C. pecorum molecular studies, to date, only four complete and 16 draft genomes have been published. Of those, only five draft genomes are from koalas. Here, using whole-genome sequencing and a comparative genomics approach, we describe the first two complete C. pecorum genomes collected from diseased koalas. A de novo assembly of DBDeUG_2018 and MC/MarsBar_2018 resolved the chromosomes and chlamydial plasmids each as single, circular contigs. Robust phylogenomic analyses indicate biogeographical separation between strains from northern and southern koala populations, and between strains infecting koala and livestock hosts. Comparative genomics between koala strains identified new, unique, and shared loci that accumulate single-nucleotide polymorphisms and separate between northern and southern, and within northern koala strains. Furthermore, we predicted novel type III secretion system effectors. This investigation constitutes a comprehensive genome-wide comparison between C. pecorum from koalas and provides improvements to annotations of a C. pecorum reference genome. These findings lay the foundations for identifying and understanding host specificity and adaptation behind chlamydial infections affecting koalas.

Diagnosis of Chlamydia trachomatis genital infections in the era of genomic medicine

PURPOSE

Chlamydial genital infections constitute significant sexually transmitted infections worldwide. The often asymptomatic status of C. trachomatis (CT) infections leads to an increased burden on human reproductive health, especially in middle- and low-income settings. Early detection and management of these infections could play a decisive role in controlling this public health burden. The objective of this review is to provide an insight into the evolution of diagnostic methods for CT infections through the development of new molecular technologies, emphasizing on -omics' technologies and their significance as diagnostic tools both for effective patient management and control of disease transmission.

METHODS

Narrative review of the diagnostic methodologies of CT infections and the impact of the introduction of -omics' technologies on their diagnosis by review of the literature.

RESULTS

Various methodologies are discussed with respect to working principles, required specifications, advantages, and disadvantages. Implementing the most accurate methods in diagnosis is highlighted as the cornerstone in managing CT infections.

CONCLUSION

Diagnostics based on -omics' technologies are considered to be the most pertinent modalities in CT testing when compared to other available methods. There is a need to modify these effective and accurate diagnostic tools in order to render them more available and feasible in all settings, especially aiming on turning them to rapid point-of-care tests for effective patient management and disease control.

Comparison of genovars and Chlamydia trachomatis infection loads in ocular samples from children in two distinct cohorts in Sudan and Morocco

Trachoma is a blinding disease caused by repeated conjunctival infection with different Chlamydia trachomatis (Ct) genovars. Ct B genovars have been associated with more severe trachoma symptoms. Here, we investigated associations between Ct genovars and bacterial loads in ocular samples from two distinct geographical locations in Africa, which are currently unclear. We tested ocular swabs from 77 Moroccan children (28 with trachomatous inflammation-follicular (TF) and 49 healthy controls), and 96 Sudanese children (54 with TF and 42 healthy controls) with a Ct-specific real-time polymerase chain reaction (PCR) assay. To estimate bacterial loads, Ct-positive samples were further processed by multiplex real-time qPCR to amplify the chromosomal outer membrane complex B and plasmid open reading frame 2 of Ct. Genotyping was performed by PCR-based amplification of the outer membrane protein A gene (~1120 base pairs) of Ct and Sanger sequencing. Ct-positivities among the Moroccan and Sudanese patient groups were 60·7% and 31·5%, respectively. Significantly more Sudanese patients than Moroccan patients were genovar A-positive. In contrast, B genovars were significantly more prevalent in Moroccan patients than in Sudanese patients. Significantly higher Ct loads were found in samples positive for B genovars (598596) than A genovar (51005). Geographical differences contributed to the distributions of different ocular Ct genovars. B genovars may induce a higher bacterial load than A genovars in trachoma patients. Our findings emphasize the importance of conducting broader studies to elucidate if the noted difference in multiplication abilities are genovar and/or endemicity level dependent.

We investigated the association between different Ct genovars, the approximate load of infection, and the distribution of Chlamydia genovars by comparing samples from one trachoma-endemic area (i.e., the city of El-Gadaref in Al Qadarif, Sudan) and one previously endemic area (i.e., the Zagora Province in Morocco), currently considered as non-endemic. This study is the first to reveal a significant difference between the genome copy numbers of Ct genovar A and B/Ba in children with TF. Evidence that Ct is still circulating in rural foci of countries like Morocco that are no longer considered endemic implies that the continuation of the trachoma surveillance must be warranted in future to avoid further spreading of Ct. The clinical significance of different infectious loads in the development of sequelae has to be determined as well as whether these differences are genovar specific or related to the given endemicity level.

High-resolution genotyping of Lymphogranuloma Venereum (LGV) strains of Chlamydia trachomatis in London using multi-locus VNTR analysis-ompA genotyping (MLVA-ompA)

Background

Lymphogranuloma venereum (LGV) is caused by Chlamydia trachomatis strains with ompA genotypes L1 to L3. An LGV epidemic associated with the L2b genotype has emerged in the past few decades amongst men who have sex with men (MSM). C. trachomatis genotypes can be discriminated by outer membrane protein A gene (ompA) sequencing, however this method has limited resolution. This study employed a high-resolution genotyping method, namely, multi-locus tandem repeat (VNTR) analysis with ompA sequencing (MLVA-ompA), to assess the distribution of LGV MLVA-ompA genotypes amongst individuals attending genitourinary medicine (GUM) clinics in London.

Methods

Clinical specimens were collected from individuals attending eight London-based GUM clinics. Specimens that tested positive for C. trachomatis by commercial nucleic acid amplification test (NAAT) were confirmed as LGV by pmpH real-time PCR. LGV-positive DNA extracts were subsequently genotyped using MLVA-ompA.

Results

Two hundred and thirty DNA extracts were confirmed as LGV, and 162 (70%) yielded complete MLVA-ompA genotypes. Six LGV MLVA-ompA genotypes were identified: 1.9.2b-L2, 1.9.3b-L2b, 1.9.2b-L2b, 1.9.2b-L2b/D, 1.4a.2b-L2b, and 5.9.2b-L1. The following LGV ompA genotypes were identified (in descending order of abundance): L2, L2b, L2b/D, and L1. Eight ompA sequences with the hybrid L2b/D profile were detected. The hybrid sequence was identical to the ompA of a recombinant L2b/D strain detected in Portugal in 2017.

Conclusions

The L2 ompA genotype was found to predominate in the London study population. The study detected an unusual hybrid L2b/D ompA profile that was previously reported in Portugal. We recommend further monitoring and surveillance of LGV strains within the UK population.

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.

Reproductive system infections in women: lower genital tract syndromes

Gynecological and obstetrical infectious diseases are an important component of women's health. A system approach to gynecological and obstetrical infection helps unify and classify microbial etiology and pathogenesis within a clinical anatomical framework of lower and upper genital tract syndromes. The reproductive system of women includes the vulva, vagina, cervix, uterus, fallopian tubes and ovaries. During pregnancy, additional tissues include the chorioamnion and placenta together with the fetus and amniotic fluid. We review in two parts reproductive system infection syndromes in women using selected research results to illustrate the clinical utility of the system approach in terms of diagnosis, treatment and prevention. We conclude that a reproductive system perspective will lead to improvements in understanding, management and prevention of these diseases.

The growing repertoire of genetic tools for dissecting chlamydial pathogenesis

Multiple species of obligate intracellular bacteria in the genus Chlamydia are important veterinary and/or human pathogens. These pathogens all share similar biphasic developmental cycles and transition between intracellular vegetative reticulate bodies and infectious elementary forms, but vary substantially in their host preferences and pathogenic potential. A lack of tools for genetic engineering of these organisms has long been an impediment to the study of their biology and pathogenesis. However, the refinement of approaches developed in C. trachomatis over the last 10 years, and adaptation of some of these approaches to other Chlamydia spp. in just the last few years, has opened exciting new possibilities for studying this ubiquitous group of important pathogens.

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.

HCK and ABAA: A Newly Designed Pipeline to Improve Fungi Metabarcoding Analysis

Introduction

The fungi ITS sequence length dissimilarity, non-specific amplicons, including chimaera formed during Polymerase Chain Reaction (PCR), added to sequencing errors, create bias during similarity clustering and abundance estimation in the downstream analysis. To overcome these challenges, we present a novel approach, Hierarchical Clustering with Kraken (HCK), to classify ITS1 amplicons and Abundance-Base Alternative Approach (ABAA) pipeline to detect and filter non-specific amplicons in fungi metabarcoding sequencing datasets.

Materials and Methods

We compared the performances of both pipelines against QIIME, KRAKEN, and DADA2 using publicly available fungi ITS mock community datasets and using BLASTn as a reference. We calculated the Precision, Recall, F-score using the True-Positive, False-positive, and False-negative estimation. Alpha diversity (Chao1 and Shannon metrics) was also used to evaluate the diversity estimation of our method.

Results

The analysis shows that ABAA reduced the number of false-positive with all metabarcoding methods tested, and HCK increases precision and recall. HCK, coupled with ABAA, improves the F-score and bring alpha diversity metric value close to that of the BLASTn alpha diversity values when compared to QIIME, KRAKEN, and DADA2.

Conclusion

The developed HCK-ABAA approach allows better identification of the fungi community structures while avoiding use of a reference database for non-specific amplicons filtration. It results in a more robust and stable methodology over time. The software can be downloaded on the following link: https://bitbucket.org/GottySG36/hck/src/master/.

Presence of Chlamydia trachomatis DNA in the amniotic fluid in women with preterm prelabor rupture of membranes

OBJECTIVE

The primary aim of this study was to assess the rate and load of amniotic fluid Chlamydia trachomatis DNA and their associations with intra-amniotic infection and intra-uterine inflammatory complications in women with preterm prelabor rupture of membranes (PPROM). The secondary aim was to assess the short-term morbidity of newborns from PPROM pregnancies complicated by amniotic fluid C. trachomatis DNA.

METHODS

A retrospective study of 788 women with singleton pregnancies complicated by PPROM between 24 + 0 and 36 + 6 weeks of gestation was performed. Transabdominal amniocenteses were performed at the time of admission. C. trachomatis DNA in the amniotic fluid was assessed by real-time polymerase chain reaction using a commercial AmpliSens® C. trachomatis/Ureaplasma/Mycoplasma hominis-FRT kit, and the level of Ct DNA was quantified.

RESULTS

Amniotic fluid C. trachomatis DNA complicated 2% (16/788) of the PPROM pregnancies and was present in very low loads (median 57 copies DNA/mL). In addition to amniotic fluid C. trachomatis DNA, other bacteria were detected in 62% (10/16) of the C. trachomatis DNA-complicated PPROM pregnancies. Amniotic fluid C. trachomatis DNA was associated with intra-amniotic infection, histologic chorioamnionitis (HCA), and funisitis in 31%, 47%, and 33%, respectively. The presence of C. trachomatis DNA accompanied by Ureaplasma species in the amniotic fluid was associated with a higher rate of HCA than the presence of amniotic fluid C. trachomatis DNA alone. The composite neonatal morbidity in newborns from PPROM pregnancies with amniotic fluid C. trachomatis DNA was 31%.

CONCLUSION

The presence of C. trachomatis DNA in the amniotic fluid is a relatively rare condition in PPROM. Amniotic fluid C. trachomatis DNA in PPROM is not related to intensive intra-amniotic and intr-auterine inflammatory responses or adverse short-term neonatal outcomes.

Whole-Genome Enrichment and Sequencing of Chlamydia trachomatis Directly from Patient Clinical Vaginal and Rectal Swabs

Chlamydia trachomatis, an obligately intracellular bacterium, is the most prevalent cause of bacterial sexually transmitted infections (STIs) worldwide. Numbers of U.S. infections of the urogenital tract and rectum have increased annually. Because C. trachomatis is not easily cultured, comparative genomic studies are limited, restricting our understanding of strain diversity and emergence among populations globally. While Agilent SureSelect^XT target enrichment RNA bait libraries have been developed for whole-genome enrichment and sequencing of C. trachomatis directly from clinical urine, vaginal, conjunctival, and rectal samples, public access to these libraries is not available. We therefore designed an RNA bait library (34,795 120-mer probes based on 85 genomes, versus 33,619 probes using 74 genomes in a previous one) to augment organism sequencing from clinical samples that can be shared with the scientific community, enabling comparison studies. We describe the library and limit of detection for genome copy input, and we present results of 100% efficiency and high-resolution determination of recombination and identical genomes within vaginal-rectal specimen pairs in women. This workflow provides a robust approach for discerning genomic diversity and advancing our understanding of the molecular epidemiology of contemporary C. trachomatis STIs across sample types, geographic populations, sexual networks, and outbreaks associated with proctitis/proctocolitis among women and men who have sex with men.IMPORTANCE Chlamydia trachomatis is an obligate intracellular bacterium that is not easily cultured, which limits our understanding of urogenital and rectal C. trachomatis transmission and impact on morbidity. To provide a publicly available workflow for whole-genome target enrichment and sequencing of C. trachomatis directly from clinical urine, vaginal, conjunctival, and rectal specimens, we developed and report on an RNA bait library to enrich the organism from clinical samples for sequencing. We demonstrate an increased efficiency in the percentage of reads mapping to C. trachomatis and identified recombinant and identical C. trachomatis genomes in paired vaginal-rectal samples from women. Our workflow provides a robust genomic epidemiologic approach to advance our understanding of C. trachomatis strains causing ocular, urogenital, and rectal infections and to explore geo-sexual networks, outbreaks of colorectal infections among women and men who have sex with men, and the role of these strains in morbidity.

One is not enough: On the effects of reference genome for the mapping and subsequent analyses of short-reads

Mapping of high-throughput sequencing (HTS) reads to a single arbitrary reference genome is a frequently used approach in microbial genomics. However, the choice of a reference may represent a source of errors that may affect subsequent analyses such as the detection of single nucleotide polymorphisms (SNPs) and phylogenetic inference. In this work, we evaluated the effect of reference choice on short-read sequence data from five clinically and epidemiologically relevant bacteria (Klebsiella pneumoniae, Legionella pneumophila, Neisseria gonorrhoeae, Pseudomonas aeruginosa and Serratia marcescens). Publicly available whole-genome assemblies encompassing the genomic diversity of these species were selected as reference sequences, and read alignment statistics, SNP calling, recombination rates, dN/dS ratios, and phylogenetic trees were evaluated depending on the mapping reference. The choice of different reference genomes proved to have an impact on almost all the parameters considered in the five species. In addition, these biases had potential epidemiological implications such as including/excluding isolates of particular clades and the estimation of genetic distances. These findings suggest that the single reference approach might introduce systematic errors during mapping that affect subsequent analyses, particularly for data sets with isolates from genetically diverse backgrounds. In any case, exploring the effects of different references on the final conclusions is highly recommended.

Mapping consists in the alignment of reads (i.e., DNA fragments) obtained through high-throughput genome sequencing to a previously assembled reference sequence. It is a common practice in genomic studies to use a single reference for mapping, usually the ‘reference genome’ of a species—a high-quality assembly. However, the selection of an optimal reference is hindered by intrinsic intra-species genetic variability, particularly in bacteria. It is known that genetic differences between the reference genome and the read sequences may produce incorrect alignments during mapping. Eventually, these errors could lead to misidentification of variants and biased reconstruction of phylogenetic trees (which reflect ancestry between different bacterial lineages). To our knowledge, this is the first work to systematically examine the effect of different references for mapping on the inference of tree topology as well as the impact on recombination and natural selection inferences. Furthermore, the novelty of this work relies on a procedure that guarantees that we are evaluating only the effect of the reference. This effect has proved to be pervasive in the five bacterial species that we have studied and, in some cases, alterations in phylogenetic trees could lead to incorrect epidemiological inferences. Hence, the use of different reference genomes may be prescriptive to assess the potential biases of mapping.

High-resolution multilocus sequence typing for Chlamydia trachomatis: improved results for clinical samples with low amounts of C. trachomatis DNA

Background

Several Multilocus Sequence Typing (MLST) schemes have been developed for Chlamydia trachomatis. Bom’s MLST scheme for MLST is based on nested PCR amplification and sequencing of five hypervariable genes and ompA. In contrast to other Chlamydia MLST schemes, Bom’s MLST scheme gives higher resolution and phylogenetic trees that are comparable to those from whole genome sequencing. However, poor results have been obtained with Bom’s MLST scheme in clinical samples with low concentrations of Chlamydia DNA.

Results

In this work, we present an improved version of the scheme that is based on the same genes and MLST database as Bom’s MLST scheme, but with newly designed primers for nested-1 and nested-2 steps under stringent conditions. Furthermore, we introduce a third primer set for the sequencing step, which considerably improves the performance of the assay. The improved primers were tested in-silico using a dataset of 141 Whole Genome Sequences (WGS) and in a comparative analysis of 32 clinical samples. Based on cycle threshold and melting curve analysis values obtained during Real-Time PCR of nested-1 & 2 steps, we developed a simple scoring scheme and flow chart that allow identification of reaction inhibitors as well as to predict with high accuracy amplification success. The improved MLST version was used to obtain a genovars distribution in patients attending an STI clinic in Tel Aviv.

Conclusions

The newly developed MLST version showed great improvement of assay results for samples with very low concentrations of Chlamydia DNA. A similar concept could be applicable to other MLST schemes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-020-02077-y.

Isolation and Identification of Chlamydia abortus from Aborted Ewes in Sulaimani Province, Northern Iraq

Abortion in small ruminants is a significant problem in Iraq and causes severe economic losses in sheep farms. Chlamydia abortus causes enzootic abortion in ewes and is associated with reproductive problems in sheep in Sulaimani province – Northern Iraq. During a lambing season in 2017, abortion was widespread among several sheep flocks in different regions of Sulaimani (Kalar, Said Sadiq, and Chamchamal), and C. abortus was one of the causes. Accordingly, we carried out this study to isolate and identify C. abortus in aborted ewes in these regions. We collected 30 samples of aborted fetuses from five herds in which abortions had been observed. The pathogen isolation was done by inoculation into embryonated chicken eggs and conventional PCR was used to identify C. abortus in clinical specimens. C. abortus was identified in one of the 30 aborted fetuses (3.33%) from the Kalar district, and all the remaining 29 samples (96.66%) were found positive to Brucella abortus. The gene ompA encoding the outer membrane protein of C. abortus was sequenced and got the accession number MK643153 in NCBI GenBank. The sequence was named C. abortus strain Sul/2017. Our isolate showed 99.79% homology with Sul/014 (accession No. KY399850) and differed from the latter by two amino acid substitutions at E115K and K259N. The topology of the phylogenetic tree based on the ompA gene showed that the isolate belongs to C. abortus and has a common ancestor with isolates of sheep in Iraq and Tunisia with accession numbers KY399850 and HQ62243, respectively.

Chlamydiae - What's New?

This paper provides an overview of the current knowledge of chlamydiae. These intracellular microorganisms belonging to the Chlamydiaceae family are widely distributed throughout the world. Constant development of culture-independent approaches for characterisation of microbial genomes enables new discoveries in the field of Chlamydia. The number of new taxa is continuously increasing as well as the range of hosts. New species and genotypes are constantly being discovered, particularly new avian and reptilian agents, which are discussed in this article. Interestingly, wild animals are the main hosts for new Chlamydia species including different species of bird, turtle and snake. The availability of next-generation sequencing opens up a new prospect for research and leads to deeper knowledge of these interesting microorganisms about which much is still to discover.

Genomics of Ocular Chlamydia trachomatis after 5 years of SAFE interventions for trachoma in Amhara, Ethiopia

To eliminate trachoma as a public health problem, the WHO recommends the SAFE (Surgery, Antibiotics, Facial cleanliness, and Environmental improvement) strategy. As part of the SAFE strategy in the Amhara Region, Ethiopia, the Trachoma Control Program distributed over 124 million doses of antibiotic between 2007 and 2015. Despite this, trachoma remained hyperendemic in many districts and a considerable level of Chlamydia trachomatis (Ct) infection was evident. We utilised residual material from Abbott m2000 Ct diagnostic tests to sequence 99 ocular Ct samples from Amhara and investigated the role of Ct genomic variation in continued transmission of Ct. Sequences were typical of ocular Ct, at the whole-genome level and in tissue tropism-associated genes. There was no evidence of macrolide-resistance in this population. Polymorphism around ompA gene was associated with village-level trachomatous inflammation-follicular prevalence. Greater ompA diversity at the district-level was associated with increased Ct infection prevalence. We found no evidence for Ct genomic variation contributing to continued transmission of Ct after treatment, adding to evidence that azithromycin does not drive acquisition of macrolide resistance in Ct. Increased Ct infection in areas with more ompA variants requires longitudinal investigation to understand what impact this may have on treatment success and host immunity.

The simplified trachoma grading system, amended

A simplified grading system for trachoma was published by the World Health Organization (WHO) in 1987. Intended for use by non-specialist personnel working at community level, the system includes five signs, each of which can be present or absent in any eye: (i) trachomatous trichiasis; (ii) corneal opacity; (iii) trachomatous inflammation—follicular; (iv) trachomatous inflammation—intense; and (v) trachomatous scarring. Though neither perfectly sensitive nor perfectly specific for trachoma, these signs have been essential tools for identifying populations that need interventions to eliminate trachoma as a public health problem. In 2018, at WHO’s 4th global scientific meeting on trachoma, the definition of one of the signs, trachomatous trichiasis, was amended to exclude trichiasis that affects only the lower eyelid. This paper presents the amended system, updates its presentation, offers notes on its use and identifies areas of ongoing debate.

Effect of Short-Term Antimicrobial Therapy on the Tolerance and Antibiotic Resistance of Multidrug-Resistant Staphylococcus capitis

Background

Bacteria undergo adaptive mutation in the host. However, the specific effect of antimicrobial use on bacterial evolution and genome mutations related to bacterial survival within a patient is unclear.

Materials and Methods

Three S. capitis strains were sequentially isolated from cerebrospinal fluid of a clinical inpatient. Antimicrobial susceptibility, growth rate, biofilm formation and whole blood survival of these strains were measured. Relative fitness was calculated. The virulence was examined in the Galleria mellonella model. Whole-genome sequencing and in silico analysis were performed to explore the genetic mechanisms of the changes in antimicrobial resistance phenotype. Hypothetical proteins are cloned, expressed and characterized by detection the susceptibility to gentamycin.

Results

The first isolate was susceptible to rifampin (MIC=0.25 μg/mL), resistant to gentamicin (MIC=16 μg/mL), while the later two isolates were resistant to rifampin (MIC >64 μg/mL), susceptible to gentamicin (MIC=4 μg/mL). For the latter two strains, compared to the first, frameshift mutation in a hypothetical protein encoding gene and base substitutions (in genes saeR, moaA and rpoB) were discovered. The mutation of rpoB gene caused rifampicin resistance. Mutations in saeR, moaA and hypothetical gene are associated with changes in other biological traits. Amino acid sequence-based structure and function identification of the hypothetical protein indicated that a mutation in the encoding gene might be associated with altered aminoglycoside susceptibility. Growth curve showed that the later two isolates grew faster than the first isolate with a positive fitness advantage of 13.5%, and 14.8%, accordingly. Biofilm form ability and whole blood survival of the derivative mutants were also enhanced. No significant differences of virulence in the G. mellonella model were observed.

Conclusion

We report here for the first time that short-term clinical antibiotic use was associated with resistance mutations, collateral sensitivity, and positive in vivo fitness advantages to S. capitis during infection.

Chlamydiae in the Environment

Chlamydiae have been known for more than a century as major pathogens of humans. Yet they are also found ubiquitously in the environment where they thrive within protists and in an unmatched wide range of animals. This review summarizes recent advances in understanding chlamydial diversity and distribution in nature. Studying these environmental chlamydiae provides a novel perspective on basic chlamydial biology and evolution. A picture is beginning to emerge with chlamydiae representing one of the evolutionarily most ancient and successful groups of obligate intracellular bacteria.

Methicillin-Resistant Staphylococcus aureus: Risk for General Infection and Endocarditis Among Athletes

The first studies on Staphylococcus aureus (SA) infections in athletes were conducted in the 1980s, and examined athletes that perform in close physical contact, with particular attention to damaged or infected skin. Recent studies have used molecular epidemiology to shed light on the transmission of SA in professional athletes. These studies have shown that contact between athletes is prolonged and constant, and that these factors influence the appearance of infections caused by SA. These results support the need to use sanitary measures designed to prevent the appearance of SA infections. The factors triggering the establishment of SA within professional sports groups are the nasal colonization of SA, contact between athletes and sweating. Hence, there is a need to use the most modern molecular typing methods to evaluate the appearance of cutaneous SA disease. This review aims to summarize both the current SA infections known in athletes and the diagnostic methods employed for recognition, pointing to possible preventive strategies and the factors that can act as a springboard for the appearance of SA and subsequent transmission between athletes.

Sensitivity, specificity, inclusivity and exclusivity of the updated Aptima Combo 2 assay, which provides detection coverage of the new diagnostic-escape Chlamydia trachomatis variants

Background

Four new variants of Chlamydia trachomatis (nvCTs), detected in several countries, cause false-negative or equivocal results using the Aptima Combo 2 assay (AC2; Hologic). We evaluated the clinical sensitivity and specificity, as well as the analytical inclusivity and exclusivity of the updated AC2 for the detection of CT and Neisseria gonorrhoeae (NG) on the automated Panther system (Hologic).

Methods

We examined 1004 clinical AC2 samples and 225 analytical samples spiked with phenotypically and/or genetically diverse NG and CT strains, and other potentially cross-reacting microbial species. The clinical AC2 samples included CT wild type (WT)-positive (n = 488), all four described AC2 diagnostic-escape nvCTs (n = 170), NG-positive (n = 214), and CT/NG-negative (n = 202) specimens.

Results

All nvCT-positive samples (100%) and 486 (99.6%) of the CT WT-positive samples were positive in the updated AC2. All NG-positive, CT/NG-negative, Trichomonas vaginalis (TV)-positive, bacterial vaginosis-positive, and Candida-positive AC2 specimens gave correct results. The clinical sensitivity and specificity of the updated AC2 for CT detection was 99.7 and 100%, respectively, and for NG detection was 100% for both. Examining spiked samples, the analytical inclusivity and exclusivity were 100%, i.e., in clinically relevant concentrations of spiked microbe.

Conclusions

The updated AC2, including two CT targets and one NG target, showed a high sensitivity, specificity, inclusivity and exclusivity for the detection of CT WT, nvCTs, and NG. The updated AC2 on the fully automated Panther system offers a simple, rapid, high-throughput, sensitive, and specific diagnosis of CT and NG, which can easily be combined with detection of Mycoplasma genitalium and TV.

Protein-protein interactions of HPV-Chlamydia trachomatis-human and their potential in cervical cancer

Aim: HPV is an important cause of cervical cancer, but Chlamydia trachomatis (CT) is suspiciously involved in this disease ranging from direct to its involvement as a cofactor with HPV. We performed this study to understand the interaction of HPV and C. trachomatis with humans and its contribution to cervical cancer. Materials & methods: Host-pathogen and pathogen-pathogen protein-protein interaction maps of HPV/CT/human were prepared and compared to analyze interactions during single/coinfection of C. trachomatis and HPV. The interacting human proteins were detected by their involvement in cervical cancer. Results: C. trachomatis may interact with several cancer associated proteins while HPV and C. trachomatis largely interact with different human proteins, suggesting different pathogenesis. Conclusion: C. trachomatis coinfection with HPV may modulate cervical cancer development.

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.

Chlamydia trachomatis: when the virulence-associated genome backbone imports a prevalence-associated major antigen signature

Chlamydia trachomatis is the most prevalent sexually transmitted bacterium worldwide and the causative agent of trachoma. Its strains are classified according to their ompA genotypes, which are strongly linked to differential tissue tropism and disease outcomes [ocular disease, urogenital disease and lymphogranuloma venereum (LGV)]. While the genome-based species phylogenetic tree presents four main clades correlating with tropism/prevalence, namely ocular, LGV, urogenital T1 (more prevalent genotypes) and urogenital T2 (less prevalent genotypes), inter-clade exchange of ompA is considered a rare phenomenon probably mediating marked tropism alterations. An LGV epidemic, associated with the clonal expansion of the L2b genotype, has emerged in the last few decades, raising concerns particularly due to its atypical clinical presentation (ulcerative proctitis) and circulation among men who have sex with men (MSM). Here, we report an LGV outbreak, mostly affecting human immunodeficiency virus-positive MSM engaging in high-risk sexual practices, caused by an L2b strain with a rather unique non-LGV ompA signature that precluded the laboratory notification of this outbreak as LGV. C. trachomatis whole-genome capture and sequencing directly from clinical samples was applied to deeply characterize the genomic backbone of this novel LGV outbreak-causing clone. It revealed a chimeric genome structure due to the genetic transfer of ompA and four neighbouring genes from a serovar D/Da strain, likely possessing the genomic backbone associated with the more prevalent urogenital genotypes (T1 clade), to an LGV (L2b) strain. The hybrid L2b/D-Da strain presents the adhesin and immunodominant antigen MOMP (major outer membrane protein) (encoded by ompA) with an epitope repertoire typical of non-invasive genital strains, while keeping the genome-dispersed virulence fingerprint of a classical LGV strain. As previously reported for inter-clade ompA exchange among non-LGV clades, this novel C. trachomatis genomic mosaic involving a contemporary epidemiologically and clinically relevant LGV strain may have implications on its transmission, tissue tropism and pathogenic capabilities. The emergence of variants with epidemic and pathogenic potential highlights the need for more focused surveillance strategies to capture C. trachomatis evolution in action.

Bioinformatics approach to understand the mode of microbial pathogenesis of Chlamydia trachomatis and their implications in gynecologic malignancy

Chlamydia trachomatis has a say on the target gene i.e., modulating the expression of target gene in the host so that it is given protection from the immune cells and so its survival and replication are not arrested by the host. The current study reports a wide range of C. trachomatis proteins that target the cellular as well as sub-cellular components of the host in gynecologic malignancy. Various bioinformatics tools was used to conduct an in-depth analysis on nuclear and eukaryotic sub cellular localization signal to find the sequences of the predicted proteins of C. trachomatis strain G. A total of 411 proteins was identified with 79.54% maximum expected accuracy and 51.02% least expected accuracy. There were uneven prediction of proteins along with redundancies between BaCeILo and HSLpred in the determination of sub-cellular localization of the CT proteins. The highest molecular weight proteins (>80 kDa) were observed to be the targeted proteins to nucleus of host cell. There was no constant patterns observed in the values of isoelectric point (pI) in case of mitochondrial targeting. The expression of eight proteins were significant with different fold changes. The in-silico study provided much detailed insights for further research in gynecological cancer. However, further experiments should be conducted to validate the specificity and confirmatory roles played by these predicted proteins in carcinogenesis.

Whole-Genome Sequencing of Chlamydia trachomatis Directly from Human Samples

Whole-genome sequencing is a powerful, high-resolution tool that can be used to generate accurate data on bacterial population structure, phylogeography and mutations associated with antimicrobial resistance. The ability to sequence pathogen genomes directly from clinical specimens, without the requirement for in vitro culturing, is attractive in terms of time- and labor-saving, especially in the case of slow growing, or obligate intracellular pathogens, such as Chlamydia trachomatis. However clinical samples typically contain too low levels of pathogen nucleic acid, plus relatively high levels of human and natural microbiota DNA/RNA, to make this a viable option. Using a combination of whole-genome enrichment and deep sequencing, which has been proven to be a nonmutagenic approach, we can capture all known variation found within C. trachomatis genomes. The method is a consistent and sensitive tool that enables rapid whole-genome sequencing of C. trachomatis directly from clinical samples and has the potential to be adapted to other pathogens with a similar clonal nature.

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.

Genome-based classification of Burkholderia cepacia complex provides new insight into its taxonomic status

Background

Accurate classification of different Burkholderia cepacia complex (BCC) species is essential for therapy, prognosis assessment and research. The taxonomic status of BCC remains problematic and an improved knowledge about the classification of BCC is in particular needed.

Methods

We compared phylogenetic trees of BCC based on 16S rRNA, recA, hisA and MLSA (multilocus sequence analysis). Using the available whole genome sequences of BCC, we inferred a species tree based on estimated single-copy orthologous genes and demarcated species of BCC using dDDH/ANI clustering.

Results

We showed that 16S rRNA, recA, hisA and MLSA have limited resolutions in the taxonomic study of closely related bacteria such as BCC. Our estimated species tree and dDDH/ANI clustering clearly separated 116 BCC strains into 36 clusters. With the appropriate reclassification of misidentified strains, these clusters corresponded to 22 known species as well as 14 putative novel species.

Conclusions

This is the first large-scale and systematic study of the taxonomic status of the BCC and could contribute to further insights into BCC taxonomy. Our study suggested that conjunctive use of core phylogeny based on single-copy orthologous genes, as well as pangenome-based dDDH/ANI clustering would provide a preferable framework for demarcating closely related species.

Reviewer

This article was reviewed by Dr. Xianwen Ren.