Chlamydia pan-genomic analysis reveals balance between host adaptation and selective pressure to genome reduction

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

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

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.

New insight on chlamydiae

This article provides an overview of the current knowledge on chlamydiae, which are intracellular bacteria belonging to the Chlamydiaceae family. Whole-genome sequencing leads to great increases in the available data about Chlamydia spp. Recently, novel chlamydial taxons in various hosts living in different environments have been recognised. New species and taxons with Candidatus status have been recorded mainly in birds and reptiles. Chlamydia gallinacea is an emerging infectious agent in poultry with indirectly confirmed zoonotic potential. Recently, a new group of avian C. abortus strains with worldwide distribution in various wild bird families has been described. The definition of C. abortus species became outdated with the discovery of these strains and has been amended. It now includes two subgroups, mammalian and avian, the latter including all isolates hitherto referred to as atypical C. psittaci or C. psittaci/C. abortus intermediates.

Chromosome structure and DNA replication dynamics during the life cycle of the predatory bacterium Bdellovibrio bacteriovorus

Bdellovibrio bacteriovorus, an obligate predatory Gram-negative bacterium that proliferates inside and kills other Gram-negative bacteria, was discovered more than 60 years ago. However, we have only recently begun to understand the detailed cell biology of this proficient bacterial killer. Bdellovibrio bacteriovorus exhibits a peculiar life cycle and bimodal proliferation, and thus represents an attractive model for studying novel aspects of bacterial cell biology. The life cycle of B. bacteriovorus consists of two phases: a free-living nonreplicative attack phase and an intracellular reproductive phase. During the reproductive phase, B. bacteriovorus grows as an elongated cell and undergoes binary or nonbinary fission, depending on the prey size. In this review, we discuss: (1) how the chromosome structure of B. bacteriovorus is remodeled during its life cycle; (2) how its chromosome replication dynamics depends on the proliferation mode; (3) how the initiation of chromosome replication is controlled during the life cycle, and (4) how chromosome replication is spatiotemporally coordinated with the proliferation program.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

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

Recent Advances in Genomics-Based Approaches for the Development of Intracellular Bacterial Pathogen Vaccines

Infectious diseases continue to be a leading cause of morbidity and mortality worldwide. The majority of infectious diseases are caused by intracellular pathogenic bacteria (IPB). Historically, conventional vaccination drives have helped control the pathogenesis of intracellular bacteria and the emergence of antimicrobial resistance, saving millions of lives. However, in light of various limitations, many diseases that involve IPB still do not have adequate vaccines. In response to increasing demand for novel vaccine development strategies, a new area of vaccine research emerged following the advent of genomics technology, which changed the paradigm of vaccine development by utilizing the complete genomic data of microorganisms against them. It became possible to identify genes related to disease virulence, genetic patterns linked to disease virulence, as well as the genetic components that supported immunity and favorable vaccine responses. Complete genomic databases, and advancements in transcriptomics, metabolomics, structural genomics, proteomics, immunomics, pan-genomics, synthetic genomics, and population biology have allowed researchers to identify potential vaccine candidates and predict their effects in patients. New vaccines have been created against diseases for which previously there were no vaccines available, and existing vaccines have been improved. This review highlights the key issues and explores the evolution of vaccines. The increasing volume of IPB genomic data, and their application in novel genome-based techniques for vaccine development, were also examined, along with their characteristics, and the opportunities and obstacles involved. Critically, the application of genomics technology has helped researchers rapidly select and evaluate candidate antigens. Novel vaccines capable of addressing the limitations associated with conventional vaccines have been developed and pressing healthcare issues are being addressed.

Comparative Genomics and Pan-Genome Driven Prediction of a Reduced Genome of Akkermansia muciniphila

Akkermanisia muciniphila imparts important health benefits and is considered a next-generation probiotic. It is imperative to understand the genomic diversity and metabolic potential of the species for safer applications as probiotics. As it resides with both health-promoting and pathogenic bacteria, understanding the evolutionary patterns are crucial, but this area remains largely unexplored. Moreover, pan-genome has previously been established based on only a limited number of strains and without careful strain selection. The pan-genomics have become very important for understanding species diversity and evolution. In the current study, a systematic approach was used to find a refined pan-genome profile of A. muciniphila by excluding too-diverse strains based on average nucleotide identity-based species demarcation. The strains were divided into four phylogroups using a variety of clustering techniques. Horizontal gene transfer and recombination patterns were also elucidated. Evolutionary patterns revealed that different phylogroups were expanding differently. Furthermore, a comparative evaluation of the metabolic potential of the pan-genome and its subsections was performed. Lastly, the study combines functional annotation, persistent genome, and essential genes to devise an approach to determine a minimal genome that can systematically remove unwanted genes, including virulent factors. The selection of one strain to be used as a chassis for the prediction of a reduced genome was very carefully performed by analyzing several genomic parameters, including the number of unique genes and the resistance and pathogenic potential of the strains. The strategy could be applied to other microbes, including human-associated microbiota, towards a common goal of predicting a minimal or a reduced genome.

A Reverse Genetic Approach for Studying sRNAs in Chlamydia trachomatis

sRNAs are noncoding transcripts that play critical roles in posttranscriptional regulation in prokaryotes. In the intracellular bacterium Chlamydia, sRNAs have been identified, but functional studies have been limited to an E. coli heterologous system. We have developed an inducible sRNA overexpression system in Chlamydia trachomatis and used it to screen putative sRNAs for effects on the Chlamydia developmental cycle, which involves conversion between replicating (RB) and infectious (EB) chlamydial forms. Overexpression of 4 of 13 C. trachomatis sRNAs decreased production of infectious EBs. We performed detailed characterization of CtrR3 and CtrR7, the two sRNAs that caused the largest progeny defects in our screen. By quantifying chlamydial number and infectious progeny, and by visualizing chlamydial forms using electron microscopy, we showed that overexpression of CtrR3 prevented RB-to-EB conversion, whereas CtrR7 overexpression blocked bacterial replication. We also describe a workflow that allowed us to identify the mRNA targets of CtrR3 in Chlamydia. We first used MS2 aptamer affinity purification coupled with RNA sequencing as an unbiased approach to isolate interacting mRNAs. We then prioritized candidates based on sequence complementarity to the CtrR3 target recognition sequence, which we had identified with bioinformatic and mutational analyses. Finally, we tested putative targets with translational fusion assays in E. coli and C. trachomatis. Using this integrated approach, we provide experimental evidence that YtgB and CTL0389 are mRNA targets of CtrR3 in Chlamydia. These findings demonstrate how our C. trachomatis sRNA overexpression system can be used to investigate the functions and mRNA targets of chlamydial sRNAs.

Genomic evolution and adaptation of arthropod-associated Rickettsia

Rickettsia species are endosymbionts hosted by arthropods and are known to cause mild to fatal diseases in humans. Here, we analyse the evolution and diversity of 34 Rickettsia species using a pangenomic meta-analysis (80 genomes/41 plasmids). Phylogenomic trees showed that Rickettsia spp. diverged into two Spotted Fever groups, a Typhus group, a Canadensis group and a Bellii group, and may have inherited their plasmids from an ancestral plasmid that persisted in some strains or may have been lost by others. The results suggested that the ancestors of Rickettsia spp. might have infected Acari and/or Insecta and probably diverged by persisting inside and/or switching hosts. Pangenomic analysis revealed that the Rickettsia genus evolved through a strong interplay between genome degradation/reduction and/or expansion leading to possible distinct adaptive trajectories. The genus mainly shared evolutionary relationships with α-proteobacteria, and also with γ/β/δ-proteobacteria, cytophagia, actinobacteria, cyanobacteria, chlamydiia and viruses, suggesting lateral exchanges of several critical genes. These evolutionary processes have probably been orchestrated by an abundance of mobile genetic elements, especially in the Spotted Fever and Bellii groups. In this study, we provided a global evolutionary genomic view of the intracellular Rickettsia that may help our understanding of their diversity, adaptation and fitness.

Chlamydia trachomatis induces lncRNA MIAT upregulation to regulate mitochondria-mediated host cell apoptosis and chlamydial development

Chlamydia trachomatis persistent infection is the leading cause of male prostatitis and female genital tract diseases. Inhibition of host cell apoptosis is the key to maintaining Chlamydia survival in vivo, and long noncoding RNAs (lncRNAs) play important roles in its developmental cycle and pathogenesis. However, it is not clear how lncRNAs regulate persistent Chlamydia infection. Here, using a microarray method, we identified 1718 lncRNAs and 1741 mRNAs differentially expressed in IFN-γ-induced persistent C. trachomatis infection. Subsequently, 10 upregulated and 5 downregulated differentially expressed lncRNAs were verified by qRT-PCR to confirm the reliability of the chip data. The GO and KEGG analyses revealed that differentially regulated transcripts were predominantly involved in various signalling pathways related to host immunity and apoptosis response. Targeted silencing of three lncRNAs (MIAT, ZEB1-AS1 and IRF1) resulted in increased apoptosis rates. Furthermore, interference with lncRNA MIAT caused not only an obvious downregulation of the Bcl-2/Bax ratio but also a marked release of cytochrome c, resulting in a significantly elevated level of caspase-3 activation. Meanwhile, MIAT was involved in the regulation of chlamydial development during the persistent infection. Collectively, these observations shed light on the enormous complex lncRNA regulatory networks involved in mitochondria-mediated host cell apoptosis and the growth and development of C. trachomatis.

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.

CRISPR Interference To Inducibly Repress Gene Expression in Chlamydia trachomatis

The ability to inducibly repress gene expression is critical to the study of organisms, like Chlamydia, with reduced genomes in which the majority of genes are likely to be essential. We recently described the feasibility of a CRISPR interference (CRISPRi) system to inducibly repress gene expression in Chlamydia trachomatis. However, the initial system suffered from some drawbacks, primarily leaky expression of the anhydrotetracycline (aTc)-inducible dCas9 ortholog and plasmid instability, which prevented population-wide studies (e.g., transcript analyses) of the effects of knockdown. Here, we describe various modifications to the original system that have allowed us to measure gene expression changes within a transformed population of C. trachomatis serovar L2. These modifications include (i) a change in the vector backbone, (ii) the introduction of a weaker ribosome binding site driving dCas9 translation, and (iii) the addition of a degradation tag to dCas9 itself. With these changes, we demonstrate the ability to inducibly repress a target gene sequence, as measured by the absence of protein by immunofluorescence analysis and by decreased transcript levels. Importantly, the expression of dCas9 alone (i.e., without a guide RNA [gRNA]) had minimal impact on chlamydial growth or development. We also describe complementation of the knockdown effect by introducing a transcriptional fusion of the target gene 3' to dCas9. Finally, we demonstrate the functionality of a second CRISPRi system based on a dCas12 system that expands the number of potential chromosomal targets. These tools should provide the ability to study essential gene function in Chlamydia.

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.

Comparative Genome Analysis of 33 Chlamydia Strains Reveals Characteristic Features of Chlamydia Psittaci and Closely Related Species

To identify genome-based features characteristic of the avian and human pathogen Chlamydia(C.) psittaci and related chlamydiae, we analyzed whole-genome sequences of 33 strains belonging to 12 species. Using a novel genome analysis tool termed Roary ILP Bacterial Annotation Pipeline (RIBAP), this panel of strains was shown to share a large core genome comprising 784 genes and representing approximately 80% of individual genomes. Analyzing the most variable genomic sites, we identified a set of features of C. psittaci that in its entirety is characteristic of this species: (i) a relatively short plasticity zone of less than 30,000 nt without a tryptophan operon (also in C. abortus, C. avium, C. gallinacea, C. pneumoniae), (ii) a characteristic set of of Inc proteins comprising IncA, B, C, V, X, Y (with homologs in C. abortus, C. caviae and C. felis as closest relatives), (iii) a 502-aa SinC protein, the largest among Chlamydia spp., and (iv) an elevated number of Pmp proteins of subtype G (14 in C. psittaci, 14 in Cand. C. ibidis). In combination with future functional studies, the common and distinctive criteria revealed in this study provide important clues for understanding the complexity of host-specific behavior of individual Chlamydia spp.