Culture-independent approaches to chlamydial genomics

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.

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.

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.

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.

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.

Ecology and evolution of chlamydial symbionts of arthropods

The phylum Chlamydiae consists of obligate intracellular bacteria including major human pathogens and diverse environmental representatives. Here we investigated the Rhabdochlamydiaceae, which is predicted to be the largest and most diverse chlamydial family, with the few described members known to infect arthropod hosts. Using published 16 S rRNA gene sequence data we identified at least 388 genus-level lineages containing about 14 051 putative species within this family. We show that rhabdochlamydiae are mainly found in freshwater and soil environments, suggesting the existence of diverse, yet unknown hosts. Next, we used a comprehensive genome dataset including metagenome assembled genomes classified as members of the family Rhabdochlamydiaceae, and we added novel complete genome sequences of Rhabdochlamydia porcellionis infecting the woodlouse Porcellio scaber, and of 'Candidatus R. oedothoracis' associated with the linyphiid dwarf spider Oedothorax gibbosus. Comparative analysis of basic genome features and gene content with reference genomes of well-studied chlamydial families with known host ranges, namely Parachlamydiaceae (protist hosts) and Chlamydiaceae (human and other vertebrate hosts) suggested distinct niches for members of the Rhabdochlamydiaceae. We propose that members of the family represent intermediate stages of adaptation of chlamydiae from protists to vertebrate hosts. Within the genus Rhabdochlamydia, pronounced genome size reduction could be observed (1.49-1.93 Mb). The abundance and genomic distribution of transposases suggests transposable element expansion and subsequent gene inactivation as a mechanism of genome streamlining during adaptation to new hosts. This type of genome reduction has never been described before for any member of the phylum Chlamydiae. This study provides new insights into the molecular ecology, genomic diversity, and evolution of representatives of one of the most divergent chlamydial families.

Evaluation of Urethrotropic-Clade Meningococcal Infection by Urine Metagenomic Shotgun Sequencing

Urethral infections caused by an emerging nongroupable (NG) urethrotropic clade of Neisseria meningitidis were first reported in the United States in 2015 (the "U.S. NmNG urethritis clade"). Here, we evaluate for the presence of other urethral pathogens in men with U.S. NmNG urethritis clade infection. We evaluated 129 urine specimens collected from men at a sexual health clinic, including 33 from patients with culture-confirmed or suspected urethral N. meningitidis infection and 96 specimens in which nucleic acid amplification test detected Neisseria gonorrhoeae, Chlamydia trachomatis, both pathogens, or neither pathogen. N. meningitidis was detected first by real-time PCR, followed by metagenomic shotgun sequencing of 91 specimens to identify coinfections. N. meningitidis genomes were sequenced following selective whole-genome amplification when possible. Metagenomic sequencing detected N. meningitidis in 16 of 17 specimens from culture-confirmed N. meningitidis cases, with no coinfection by other conventional urethral pathogens. Metagenomic sequencing also detected N. meningitidis in three C. trachomatis-positive specimens, one specimen positive for both N. gonorrhoeae and C. trachomatis, and nine specimens with negative N. gonorrhoeae and C. trachomatis results, eight of which had suspected Neisseria infections. N. meningitidis from culture-confirmed N. meningitidis cases belonged to the U.S. NmNG urethritis clade, while N. meningitidis identified in other specimens belonged to multiple clonal complexes. Additional urethral pathogens were predominant in non-N. meningitidis specimens, including N. gonorrhoeae, C. trachomatis, Mycoplasma genitalium, Ureaplasma urealyticum, and herpes simplex virus 2. Coinfection with other conventional urethral pathogens is rare in men with culture-confirmed U.S. NmNG urethritis clade infection and points to the strong association of this clade with disease.

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.

Pangenomics reveals alternative environmental lifestyles among chlamydiae

Chlamydiae are highly successful strictly intracellular bacteria associated with diverse eukaryotic hosts. Here we analyzed metagenome-assembled genomes of the "Genomes from Earth's Microbiomes" initiative from diverse environmental samples, which almost double the known phylogenetic diversity of the phylum and facilitate a highly resolved view at the chlamydial pangenome. Chlamydiae are defined by a relatively large core genome indicative of an intracellular lifestyle, and a highly dynamic accessory genome of environmental lineages. We observe chlamydial lineages that encode enzymes of the reductive tricarboxylic acid cycle and for light-driven ATP synthesis. We show a widespread potential for anaerobic energy generation through pyruvate fermentation or the arginine deiminase pathway, and we add lineages capable of molecular hydrogen production. Genome-informed analysis of environmental distribution revealed lineage-specific niches and a high abundance of chlamydiae in some habitats. Together, our data provide an extended perspective of the variability of chlamydial biology and the ecology of this phylum of intracellular microbes.

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.

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.

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.

Real-Time Selective Sequencing with RUBRIC: Read Until with Basecall and Reference-Informed Criteria

The Oxford MinION, the first commercial nanopore sequencer, is also the first to implement molecule-by-molecule real-time selective sequencing or “Read Until”. As DNA transits a MinION nanopore, real-time pore current data can be accessed and analyzed to provide active feedback to that pore. Fragments of interest are sequenced by default, while DNA deemed non-informative is rejected by reversing the pore bias to eject the strand, providing a novel means of background depletion and/or target enrichment. In contrast to the previously published pattern-matching Read Until approach, our RUBRIC method is the first example of real-time selective sequencing where on-line basecalling enables alignment against conventional nucleic acid references to provide the basis for sequence/reject decisions. We evaluate RUBRIC performance across a range of optimizable parameters, apply it to mixed human/bacteria and CRISPR/Cas9-cut samples, and present a generalized model for estimating real-time selection performance as a function of sample composition and computing configuration.

Metagenomic Analysis of Fish-Associated Ca. Parilichlamydiaceae Reveals Striking Metabolic Similarities to the Terrestrial Chlamydiaceae

Chlamydiae are an example of obligate intracellular bacteria that possess highly reduced, compact genomes (1.0-3.5 Mbp), reflective of their abilities to sequester many essential nutrients from the host that they no longer need to synthesize themselves. The Chlamydiae is a phylum with a very wide host range spanning mammals, birds, fish, invertebrates, and unicellular protists. This ecological and phylogenetic diversity offers ongoing opportunities to study intracellular survival and metabolic pathways and adaptations. Of particular evolutionary significance are Chlamydiae from the recently proposed Ca. Parilichlamydiaceae, the earliest diverging clade in this phylum, species of which are found only in aquatic vertebrates. Gill extracts from three Chlamydiales-positive Australian aquaculture species (Yellowtail kingfish, Striped trumpeter, and Barramundi) were subject to DNA preparation to deplete host DNA and enrich microbial DNA, prior to metagenome sequencing. We assembled chlamydial genomes corresponding to three Ca. Parilichlamydiaceae species from gill metagenomes, and conducted functional genomics comparisons with diverse members of the phylum. This revealed highly reduced genomes more similar in size to the terrestrial Chlamydiaceae, standing in contrast to members of the Chlamydiae with a demonstrated cosmopolitan host range. We describe a reduction in genes encoding synthesis of nucleotides and amino acids, among other nutrients, and an enrichment of predicted transport proteins. Ca. Parilichlamydiaceae share 342 orthologs with other chlamydial families. We hypothesize that the genome reduction exhibited by Ca. Parilichlamydiaceae and Chlamydiaceae is an example of within-phylum convergent evolution. The factors driving these events remain to be elucidated.