Comparative Genomics in Trachoma
Trachoma, a neglected tropical disease caused by Chlamydia trachomatis (Ct) serovar A-C, is the leading infectious cause of blindness worldwide. Its impact is substantial, with approximately 2.2 million individuals suffering from visual impairment, and about six million are irreversibly blind. Africa bears the highest burden, accounting for over 86% of global trachoma cases. We investigated Ct serovar A (SvA) and serovar 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 implications for Ct evolution within the context of ocular infection. A cohort study in 2002-2003 collected ocular swab samples across nine Gambian villages during a six-month follow-up study. To study 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. 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 and whole genome phylogenetic analysis revealed that 11 samples had a highest identity to Ct strain A/HAR13 (NC_007429) and 30 had a highest identity to Ct 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 a single nucleotide polymorphism accumulation rate of 2.4 ? 10^-5/site/year for the Gambian SvA and 1.3 ? 10^-5/site/year for SvB variants (p<0.0001). Variant calling resulted in a total of 1,371 single nucleotide variants (SNVs) with a frequency > 25% in SvA sequences, and 438 SNVs in SvB sequences, from which 740 (62.8%), and 241 (66.4%) were non-synonymous, respectively. 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_0498, CTA_0166, and CTA_0948. 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 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.