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

Resistome and virulome determination in Helicobacter pylori using next-generation sequencing with target-enrichment technology

The identification of Helicobacter pylori infection from gastric biopsy samples requires PCR or bacterial cultures. However, it is difficult to culture H. pylori because it is a fragile bacterium. Next-generation sequencing (NGS) allows direct assessment of the resistome and virulome. Here we describe a new NGS method for studying the resistome and virulome of H. pylori directly from gastric biopsies, based on enrichment analyses and targeted sequencing of H. pylori DNA. In all, 19 DNA samples from human gastric biopsies that tested positive for H. pylori were analyzed. The Agilent SureSelectXT target-enrichment protocol was used with a custom bait library prior to sequencing using the Agilent MagnisDx NGS Library Prep System. NGS sequencing was performed on the Illumina iSeq 100 sequencer using RNA probes for virulence, resistance, and molecular typing genes. The method yielded significant results with a limit of detection of around 1.8e5 CFU per mL H. pylori. Mutations in the 23S rDNA sequence associated with macrolide resistance and in the quinolone resistance-determining region of gyrase A associated with levofloxacin resistance were correctly identified. The results of MLST phylogeny analyses performed after target-enrichment were consistent with those obtained via conventional Sanger sequencing. Among the cagA-positive isolates, the gene was detected correctly, and the vacA genotype was determined. In conclusion, our enrichment method enables rapid assessment of the resistome and virulome of H. pylori directly from fresh gastric biopsies.IMPORTANCEHelicobacter pylori, a bacterium that infects at least 50% of the world population, is often treated by probabilistic antimicrobial therapies due to the lack of antimicrobial resistance data provided by clinical laboratories to clinicians. However, targeted antimicrobial therapies are increasingly recommended to achieve efficient eradication with a limited impact on the gut microbiota and with fewer adverse events for the patient. Recent advancements in next-generation sequencing strategies have opened new opportunities in the diagnosis of H. pylori infection. The significance of our research is the development of a novel next-generation sequencing strategy based on target-enrichment. This approach enables the identification of the resistome and the virulome of H. pylori directly from gastric biopsies, providing clinicians with a broad overview of therapeutic options.

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.

Evaluating methods for genome sequencing of Chlamydia trachomatis and other sexually transmitted bacteria directly from clinical swabs

Rates of bacterial sexually transmitted infections (STIs) are rising, and accessing their genomes provides information on strain evolution, circulating strains and encoded antimicrobial resistance (AMR). Notable pathogens include Chlamydia trachomatis (CT), Neisseria gonorrhoeae (NG) and Treponema pallidum (TP), globally the most common bacterial STIs. Mycoplasmoides (formerly Mycoplasma) genitalium (MG) is also a bacterial STI that is of concern due to AMR development. These bacteria are also fastidious or hard to culture, and standard sampling methods lyse bacteria, completely preventing pathogen culture. Clinical samples contain large amounts of human and other microbiota DNA. These factors hinder the sequencing of bacterial STI genomes. We aimed to overcome these challenges in obtaining whole-genome sequences and evaluated four approaches using clinical samples from Argentina (39), and Switzerland (14), and cultured samples from Finland (2) and Argentina (1). First, direct genome sequencing from swab samples was attempted through Illumina deep metagenomic sequencing, showing extremely low levels of target DNA, with under 0.01% of the sequenced reads being from the target pathogens. Second, host DNA depletion followed by Illumina sequencing was not found to produce enrichment in these very low-load samples. Third, we tried a selective long-read approach with the new adaptive sequencing from Oxford Nanopore Technologies, which also did not improve enrichment sufficiently to provide genomic information. Finally, target enrichment using a novel pan-genome set of custom SureSelect probes targeting CT, NG, TP and MG followed by Illumina sequencing was successful. We produced whole genomes from 64% of CT-positive samples, from 36% of NG-positive samples and 60% of TP-positive samples. Additionally, we enriched MG DNA to gain partial genomes from 60% of samples. This is the first publication to date to utilize a pan-genome STI panel in target enrichment. Target enrichment, though costly, proved essential for obtaining genomic data from clinical samples. These data can be utilized to examine circulating strains and genotypic resistance and guide public health strategies.

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.

A novel method for the preparation of reproducible, stable, and non-infectious quality control materials for Chlamydia trachomatis nucleic acid detection

Chlamydia trachomatis (CT) is a significant sexually transmitted pathogen known to evoke severe complications, including infertility. Nucleic acid amplification tests (NAATs) are recommended by the World Health Organization to detect CT infection. Furthermore, the establishment of methods, performance validation, internal quality control, and external quality assessment for CT NAATs necessitate the utilization of quality control materials (QCs). QCs are specimens or solutions that are analyzed for quality control purposes in a test system. In this study, we established a novel cell line that stably integrates CT amplification target sequences for producing QCs for CT NAATs. Utilizing clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 technology, we integrated the CT plasmid-mediated sequence (comprising the full length of the cryptic plasmid and the major outer membrane protein gene, 9,136 bp) into the MUC4 gene of HEK293T cells. Positive clones were screened through flow cytometric sorting, single-cell culture, and PCR-based identification, followed by the establishment of stable cell lines. These cells were then processed using optimized cell preservation procedures to prepare QCs. The sequence insertion copy number was confirmed by real-time quantitative PCR. This novel CT QCs demonstrate excellent clinical applicability, non-infectiousness, quantifiability, and stability. With an integrated sequence exceeding 9 kb in length, it offers exceptional flexibility for adapting to new kit developments. Furthermore, maintaining a well-defined copy number and stable shelf life, the QCs closely aligns with the quality control requirements of CT NAATs. This study presents an innovative method for preparing QCs for CT nucleic acid detection, making a valuable contribution to improving the performance of CT NAATs.IMPORTANCEUntreated CT infections impose significant burdens on individuals and communities, underscoring the importance of early and accurate testing via CT NAATs for disease control. QCs are instrumental in identifying testing process issues. Hence, we developed a cell line integrating CT-amplified target sequences as readily accessible non-infectious QCs. These QCs boast several advantages: the integration of over 9 kb of CT sequence allows for broad applicability, allowing flexible adaptation to the development of new kits. Confirming the CT sequence copy number provides a reliable basis for QC concentration preparation and kit detection limit evaluation. Optimized preservation protocol enhances QC stability during storage, facilitating convenient shipment to clinical laboratories at ambient temperatures. In summary, our novel CT QCs offer a powerful tool for improving CT NAAT performance and present a fresh perspective on QC preparation for detecting nucleic acids from intracellular parasitic pathogens.

Strengthening Bordetella pertussis genomic surveillance by direct sequencing of residual positive specimens

Whole-genome sequencing (WGS) of microbial pathogens recovered from patients with infectious disease facilitates high-resolution strain characterization and molecular epidemiology. However, increasing reliance on culture-independent methods to diagnose infectious diseases has resulted in few isolates available for WGS. Here, we report a novel culture-independent approach to genome characterization of Bordetella pertussis, the causative agent of pertussis and a paradigm for insufficient genomic surveillance due to limited culture of clinical isolates. Sequencing libraries constructed directly from residual pertussis-positive diagnostic nasopharyngeal specimens were hybridized with biotinylated RNA "baits" targeting B. pertussis fragments within complex mixtures that contained high concentrations of host and microbial background DNA. Recovery of B. pertussis genome sequence data was evaluated with mock and pooled negative clinical specimens spiked with reducing concentrations of either purified DNA or inactivated cells. Targeted enrichment increased the yield of B. pertussis sequencing reads up to 90% while simultaneously decreasing host reads to less than 10%. Filtered sequencing reads provided sufficient genome coverage to perform characterization via whole-genome single nucleotide polymorphisms and whole-genome multilocus sequencing typing. Moreover, these data were concordant with sequenced isolates recovered from the same specimens such that phylogenetic reconstructions from either consistently clustered the same putatively linked cases. The optimized protocol is suitable for nasopharyngeal specimens with diagnostic IS481 Ct < 35 and >10 ng DNA. Routine implementation of these methods could strengthen surveillance and study of pertussis resurgence by capturing additional cases with genomic characterization.

Target enrichment improves culture-independent detection of Neisseria gonorrhoeae and antimicrobial resistance determinants direct from clinical samples with Nanopore sequencing

Multi-drug-resistant Neisseria gonorrhoeae infection is a significant public health risk. Rapidly detecting N. gonorrhoeae and antimicrobial-resistant (AMR) determinants by metagenomic sequencing of urine is possible, although high levels of host DNA and overgrowth of contaminating species hamper sequencing and limit N. gonorrhoeae genome coverage. We performed Nanopore sequencing of nucleic acid amplification test-positive urine samples and culture-positive urethral swabs with and without probe-based target enrichment, using a custom SureSelect panel, to investigate whether selective enrichment of N. gonorrhoeae DNA improves detection of both species and AMR determinants. Probes were designed to cover the entire N. gonorrhoeae genome, with tenfold enrichment of probes covering selected AMR determinants. Multiplexing was tested in a subset of samples. The proportion of sequence bases classified as N. gonorrhoeae increased in all samples after enrichment, from a median (IQR) of 0.05 % (0.01-0.1 %) to 76 % (42-82 %), giving a corresponding median improvement in fold genome coverage of 365 times (112-720). Over 20-fold coverage, required for robust AMR determinant detection, was achieved in 13/15(87 %) samples, compared to 2/15(13 %) without enrichment. The four samples multiplexed together also achieved >20-fold genome coverage. Coverage of AMR determinants was sufficient to predict resistance conferred by changes in chromosomal genes, where present, and genome coverage also enabled phylogenetic relationships to be reconstructed. Probe-based target enrichment can improve N. gonorrhoeae genome coverage when sequencing DNA extracts directly from urine or urethral swabs, allowing for detection of AMR determinants. Additionally, multiplexing prior to enrichment provided enough genome coverage for AMR detection and reduces the costs associated with this method.

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

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

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

Introduction

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

Methods

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

Results

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

Discussion

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

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

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

Importance

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

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.

One clone to rule them all: Culture-independent genomics of Chlamydia psittaci from equine and avian hosts in Australia

Chlamydia psittaci is an avian pathogen with zoonotic potential. In Australia, C. psittaci has been well reported as a cause of reproductive loss in mares which subsequently have been the source of infection and illness in some in-contact humans. To date, molecular typing studies describe the predominant and clonal C. psittaci sequence type (ST)24 strains in horse, psittacine, and human infections. We sought to assess the clonality between ST24 strains and the emergence of equine ST24 with a comprehensive genomics approach. We used culture-independent probe-based and metagenomic whole-genome sequencing to investigate 13 C. psittaci genomes from horses, psittacines, and a pigeon from Australia. Published genomes of 36 C. psittaci strains were also used to contextualise our Australian dataset and investigate lineage diversity. We utilised a single-nucleotide polymorphism (SNP) based clustering and multi-locus sequence typing (MLST) approach. C. psittaci has four major phylogenetic groups (PG1-4) based on core-genome SNP-based phylogeny. PG1 contained clonal global and Australian equine, psittacine, and human ST24 genomes, with a median pairwise SNP distance of 68 SNPs. PG2, PG3, and PG4 had greater genomic diversity, including diverse STs collected from birds, livestock, human, and horse hosts from Europe and North America and a racing pigeon from Australia. We show that the clustering of C. psittaci by MLST was congruent with SNP-based phylogeny. The monophyletic ST24 clade has four major sub-lineages. The genomes of 17 Australian human, equine, and psittacine strains collected between 2008 and 2021 formed the predominant ST24 sub-lineage 1 (emerged circa 1979). Despite a temporal distribution of 13 years, the genomes within sub-lineage 1 had a median pairwise SNP distance of 32 SNPs, suggesting a recent population expansion or potential cross-host transmission. However, two C. psittaci genomes collected in 2015 from Victorian parrots clustered into distinct ST24 sub-lineage 4 (emerged circa 1965) with ovine strain C19/98 from Germany. This work describes a comprehensive phylogenomic characterisation of ST24 and identifies a timeline of potential bird-to-equine spillover events.

The role of tryptophan in Chlamydia trachomatis persistence

Chlamydia trachomatis (C. trachomatis) is the most common etiological agent of bacterial sexually transmitted infections (STIs) and a worldwide public health issue. The natural course with C. trachomatis infection varies widely between individuals. Some infections clear spontaneously, others can last for several months or some individuals can become reinfected, leading to severe pathological damage. Importantly, the underlying mechanisms of C. trachomatis infection are not fully understood. C. trachomatis has the ability to adapt to immune response and persist within host epithelial cells. Indoleamine-2,3-dioxygenase (IDO) induced by interferon-gamma (IFN-γ) degrades the intracellular tryptophan pool, to which C. trachomatis can respond by converting to a non-replicating but viable state. C. trachomatis expresses and encodes for the tryptophan synthase (TS) genes (trpA and trpB) and tryptophan repressor gene (trpR). Multiple genes interact to regulate tryptophan synthesis from exogenous indole, and persistent C. trachomatis can recover its infectivity by converting indole into tryptophan. In this review, we discuss the characteristics of chlamydial infections, biosynthesis and regulation of tryptophan, the relationship between tryptophan and C. trachomatis, and finally, the links between the tryptophan/IFN-γ axis and C. trachomatis persistence.

Metagenomic Shotgun Sequencing of Endocervical, Vaginal, and Rectal Samples among Fijian Women with and without Chlamydia trachomatis Reveals Disparate Microbial Populations and Function across Anatomic Sites: a Pilot Study

Chlamydia trachomatis is a sexually transmitted pathogen and a global public health concern. Little is known about the microbial composition and function across endocervical, vaginal, and rectal microbiomes in the context of C. trachomatis infection. We evaluated the microbiomes of 10 age-matched high-risk Fijian women with and without C. trachomatis using metagenomic shotgun sequencing (MSS). Lactobacillus iners and Lactobacillus crispatus dominated the vagina and endocervix of uninfected women. Species often found in higher relative abundance in bacterial vaginosis (BV)-Mageeibacillus indolicus, Prevotella spp., Sneathia spp., Gardnerella vaginalis, and Veillonellaceae spp.-were dominant in C. trachomatis-infected women. This combination of BV pathogens was unique to Pacific Islanders compared to previously studied groups. The C. trachomatis-infected endocervix had a higher diversity of microbiota and microbial profiles that were somewhat different from those of the vagina. However, community state type III (CST-III) and CST-IV predominated, reflecting pathogenic microbiota regardless of C. trachomatis infection status. Rectal microbiomes were dominated by Prevotella and Bacteroides, although four women had unique microbiomes with Gardnerella, Akkermansia, Bifidobacterium, and Brachyspira. A high level of microbial similarity across microbiomes in two C. trachomatis-infected women suggested intragenitorectal transmission. A number of metabolic pathways in the endocervix, driven by BV pathogens and C. trachomatis to meet nutritional requirements for survival/growth, 5-fold higher than that in the vagina indicated that endocervical microbial functions are likely more diverse and complex than those in the vagina. Our novel findings provide the impetus for larger prospective studies to interrogate microbial/microbiome interactions that promote C. trachomatis infection and better define the unique genitorectal microbiomes of Pacific Islanders. IMPORTANCE Chlamydia trachomatis is the primary cause of bacterial sexually transmitted infections worldwide, with a disturbing increase in annual rates. While there is a plethora of data on healthy and pathogenic vaginal microbiomes-defining microbial profiles and associations with sexually transmitted infections (STIs)-far fewer studies have similarly examined the endocervix or rectum. Further, vulnerable populations, such as Pacific Islanders, remain underrepresented in research. We investigated the microbial composition, structure, and function of these anatomic microbiomes using metagenomic shotgun sequencing among a Fijian cohort. We found, primarily among C. trachomatis-infected women, unique microbial profiles in endocervical, vaginal, and rectal microbiomes with an increased diversity and more complex microbial pathways in endocervical than vaginal microbiomes. Similarities in microbiome composition across sites for some women suggested intragenitorectal transmission. These novel insights into genitorectal microbiomes and their purported function require prospective studies to better define Pacific Islander microbiomes and microbial/microbiome interactions that promote C. trachomatis infection.