Dual RNA-seq analysis of in vitro infection multiplicity and RNA depletion methods in Chlamydia-infected epithelial cells.

Dual RNA-seq experiments examining viral and bacterial pathogens are increasing, but vary considerably in their experimental designs, such as infection rates and RNA depletion methods. Here, we have applied dual RNA-seq to Chlamydia trachomatis infected epithelial cells to examine transcriptomic responses from both organisms. We compared two time points post infection (1 and 24 h), three multiplicity of infection (MOI) ratios (0.1, 1 and 10) and two RNA depletion methods (rRNA and polyA). Capture of bacterial-specific RNA were greatest when combining rRNA and polyA depletion, and when using a higher MOI. However, under these conditions, host RNA capture was negatively impacted. Although it is tempting to use high infection rates, the implications on host cell survival, the potential reduced length of infection cycles and real world applicability should be considered. This data highlights the delicate nature of balancing host-pathogen RNA capture and will assist future transcriptomic-based studies to achieve more specific and relevant infection-related biological insights.

Chromatin accessibility dynamics of Chlamydia-infected epithelial cells.

Chlamydia are Gram-negative, obligate intracellular bacterial pathogens responsible for a broad spectrum of human and animal diseases. In humans, Chlamydia trachomatis is the most prevalent bacterial sexually transmitted infection worldwide and is the causative agent of trachoma (infectious blindness) in disadvantaged populations. Over the course of its developmental cycle, Chlamydia extensively remodels its intracellular niche and parasitises the host cell for nutrients, with substantial resulting changes to the host cell transcriptome and proteome. However, little information is available on the impact of chlamydial infection on the host cell epigenome and global gene regulation. Regions of open eukaryotic chromatin correspond to nucleosome-depleted regions, which in turn are associated with regulatory functions and transcription factor binding. We applied formaldehyde-assisted isolation of regulatory elements enrichment followed by sequencing (FAIRE-Seq) to generate temporal chromatin maps of C. trachomatis-infected human epithelial cells in vitro over the chlamydial developmental cycle. We detected both conserved and distinct temporal changes to genome-wide chromatin accessibility associated with C. trachomatis infection. The observed differentially accessible chromatin regions include temporally-enriched sets of transcription factors, which may help shape the host cell response to infection. These regions and motifs were linked to genomic features and genes associated with immune responses, re-direction of host cell nutrients, intracellular signalling, cell-cell adhesion, extracellular matrix, metabolism and apoptosis. This work provides another perspective to the complex response to chlamydial infection, and will inform further studies of transcriptional regulation and the epigenome in Chlamydia-infected human cells and tissues.

An Integrated Proteomic and Transcriptomic Analysis Reveals the Venom Complexity of the Bullet Ant Paraponera clavata.

A critical hurdle in ant venom proteomic investigations is the lack of databases to comprehensively and specifically identify the sequence and function of venom proteins and peptides. To resolve this, we used venom gland transcriptomics to generate a sequence database that was used to assign the tandem mass spectrometry (MS) fragmentation spectra of venom peptides and proteins to specific transcripts. This was performed alongside a shotgun liquid chromatography-mass spectrometry (LC-MS/MS) analysis of the venom to confirm that these assigned transcripts were expressed as proteins. Through the combined transcriptomic and proteomic investigation of Paraponera clavata venom, we identified four times the number of proteins previously identified using 2D-PAGE alone. In addition to this, by mining the transcriptomic data, we identified several novel peptide sequences for future pharmacological investigations, some of which conform with inhibitor cysteine knot motifs. These types of peptides have the potential to be developed into pharmaceutical or bioinsecticide peptides.

Early Transcriptional Landscapes of Chlamydia trachomatis-Infected Epithelial Cells at Single Cell Resolution.

Chlamydia are Gram-negative obligate intracellular bacterial pathogens responsible for a variety of disease in humans and animals worldwide. Chlamydia trachomatis causes trachoma in disadvantaged populations, and is the most common bacterial sexually transmitted infection in humans, causing reproductive tract disease. Antibiotic therapy successfully treats diagnosed chlamydial infections, however asymptomatic infections are common. High-throughput transcriptomic approaches have explored chlamydial gene expression and infected host cell gene expression. However, these were performed on large cell populations, averaging gene expression profiles across all cells sampled and potentially obscuring biologically relevant subsets of cells. We generated a pilot dataset, applying single cell RNA-Seq (scRNA-Seq) to C. trachomatis infected and mock-infected epithelial cells to assess the utility, pitfalls and challenges of single cell approaches applied to chlamydial biology, and to potentially identify early host cell biomarkers of chlamydial infection. Two hundred sixty-four time-matched C. trachomatis-infected and mock-infected HEp-2 cells were collected and subjected to scRNA-Seq. After quality control, 200 cells were retained for analysis. Two distinct clusters distinguished 3-h cells from 6- and 12-h. Pseudotime analysis identified a possible infection-specific cellular trajectory for Chlamydia-infected cells, while differential expression analyses found temporal expression of metallothioneins and genes involved with cell cycle regulation, innate immune responses, cytoskeletal components, lipid biosynthesis and cellular stress. We find that changes to the host cell transcriptome at early times of C. trachomatis infection are readily discernible by scRNA-Seq, supporting the utility of single cell approaches to identify host cell biomarkers of chlamydial infection, and to further deconvolute the complex host response to infection.

Bioinformatic analysis of bacteria and host cell dual RNA-sequencing experiments.

Bacterial pathogens subvert host cells by manipulating cellular pathways for survival and replication; in turn, host cells respond to the invading pathogen through cascading changes in gene expression. Deciphering these complex temporal and spatial dynamics to identify novel bacterial virulence factors or host response pathways is crucial for improved diagnostics and therapeutics. Dual RNA sequencing (dRNA-Seq) has recently been developed to simultaneously capture host and bacterial transcriptomes from an infected cell. This approach builds on the high sensitivity and resolution of RNA sequencing technology and is applicable to any bacteria that interact with eukaryotic cells, encompassing parasitic, commensal or mutualistic lifestyles. Several laboratory protocols have been presented that outline the collection, extraction and sequencing of total RNA for dRNA-Seq experiments, but there is relatively little guidance available for the detailed bioinformatic analyses required. This protocol outlines a typical dRNA-Seq experiment, based on a Chlamydia trachomatis-infected host cell, with a detailed description of the necessary bioinformatic analyses with currently available software tools.

Genetic diversity in the plasticity zone and the presence of the chlamydial plasmid differentiates Chlamydia pecorum strains from pigs, sheep, cattle, and koalas.

BACKGROUND: Chlamydia pecorum is a globally recognised pathogen of livestock and koalas. To date, comparative genomics of C. pecorum strains from sheep, cattle and koalas has revealed that only single nucleotide polymorphisms (SNPs) and a limited number of pseudogenes appear to contribute to the genetic diversity of this pathogen. No chlamydial plasmid has been detected in these strains despite its ubiquitous presence in almost all other chlamydial species. Genomic analyses have not previously included C. pecorum from porcine hosts. We sequenced the genome of three C. pecorum isolates from pigs with differing pathologies in order to re-evaluate the genetic differences and to update the phylogenetic relationships between C. pecorum from each of the hosts. METHODS: Whole genome sequences for the three porcine C. pecorum isolates (L1, L17 and L71) were acquired using C. pecorum-specific sequence capture probes with culture-independent methods, and assembled in CLC Genomics Workbench. The pairwise comparative genomic analyses of 16 pig, sheep, cattle and koala C. pecorum genomes were performed using several bioinformatics platforms, while the phylogenetic analyses of the core C. pecorum genomes were performed with predicted recombination regions removed. Following the detection of a C. pecorum plasmid, a newly developed C. pecorum-specific plasmid PCR screening assay was used to evaluate the plasmid distribution in 227 C. pecorum samples from pig, sheep, cattle and koala hosts. RESULTS: Three porcine C. pecorum genomes were sequenced using C. pecorum-specific sequence capture probes with culture-independent methods. Comparative genomics of the newly sequenced porcine C. pecorum genomes revealed an increased average number of SNP differences (~11 500) between porcine and sheep, cattle, and koala C. pecorum strains, compared to previous C. pecorum genome analyses. We also identified a third copy of the chlamydial cytotoxin gene, found only in porcine C. pecorum isolates. Phylogenetic analyses clustered porcine isolates into a distinct clade, highlighting the polyphyletic origin of C. pecorum in livestock. Most surprising, we also discovered a plasmid in the porcine C. pecorum genome. Using this novel C. pecorum plasmid (pCpec) sequence, a) we developed a pCpec screening assay to evaluate the plasmid distribution in C. pecorum from different hosts; and b) to characterise the pCpec sequences from available previously sequenced C. pecorum genome data. pCpec screening showed that the pCpec is common in all hosts of C. pecorum, however not all C. pecorum strains carry pCpec. CONCLUSIONS: This study provides further insight into the complexity of C. pecorum epidemiology and novel genomic regions that may be linked to host specificity. C. pecorum plasmid characterisation may aid in improving our understanding of C. pecorum pathogenesis across the variety of host species this animal pathogen infects.

Early microRNA expression profile as a prognostic biomarker for the development of pelvic inflammatory disease in a mouse model of chlamydial genital infection.

UNLABELLED: It is not currently possible to predict the probability of whether a woman with a chlamydial genital infection will develop pelvic inflammatory disease (PID). To determine if specific biomarkers may be associated with distinct chlamydial pathotypes, we utilized two Chlamydia muridarum variants (C. muridarum Var001 [CmVar001] and CmVar004) that differ in their abilities to elicit upper genital tract pathology in a mouse model. CmVar004 has a lower growth rate in vitro and induces pathology in only 20% of C57BL/6 mouse oviducts versus 83.3% of oviducts in CmVar001-infected mice. To determine if chemokine and cytokine production within 24 h of infection is associated with the outcome of pathology, levels of 15 chemokines and cytokines were measured. CmVar004 infection induced significantly lower levels of CXCL1, CXCL2, tumor necrosis factor alpha (TNF-α), and CCL2 in comparison to CmVar001 infection with similar rRNA (rs16) levels for Chlamydiae. A combination of microRNA (miRNA) sequencing and quantitative real-time PCR (qRT-PCR) analysis of 134 inflammation-related miRNAs was performed 24 h postinfection to determine if the chemokine/cytokine responses would also be reflected in miRNA expression profiles. Interestingly, 12 miRNAs (miR-135a-5p, miR298-5p, miR142-3p, miR223-3p, miR299a-3p, miR147-3p, miR105, miR325-3p, miR132-3p, miR142-5p, miR155-5p, and miR-410-3p) were overexpressed during CmVar004 infection compared to CmVar001 infection, inversely correlating with the respective chemokine/cytokine responses. To our knowledge, this is the first report demonstrating that early biomarkers elicited in the host can differentiate between two pathological variants of chlamydiae and be predictive of upper tract disease. IMPORTANCE: It is apparent that an infecting chlamydial population consists of multiple genetic variants with differing capabilities of eliciting a pathological response; thus, it may be possible to identify biomarkers specific for a given virulence pathotype. miRNAs are known to regulate genes that in turn regulate signaling pathways involved in disease pathogenesis. Importantly, miRNAs are stable and can reflect a tissue response and therefore have the potential to be biomarkers of disease severity. Currently, with respect to chlamydial infections, there is no way to predict whether an infected patient is more or less likely to develop PID. However, data presented in this study indicate that the expression of a specific miRNA profile associated with a virulent variant early in the infection course may be predictive of an increased risk of pelvic inflammatory disease, allowing more aggressive treatment before significant pathology develops.

Simultaneous transcriptional profiling of bacteria and their host cells.

We developed an RNA-Seq-based method to simultaneously capture prokaryotic and eukaryotic expression profiles of cells infected with intracellular bacteria. As proof of principle, this method was applied to Chlamydia trachomatis-infected epithelial cell monolayers in vitro, successfully obtaining transcriptomes of both C. trachomatis and the host cells at 1 and 24 hours post-infection. Chlamydiae are obligate intracellular bacterial pathogens that cause a range of mammalian diseases. In humans chlamydiae are responsible for the most common sexually transmitted bacterial infections and trachoma (infectious blindness). Disease arises by adverse host inflammatory reactions that induce tissue damage & scarring. However, little is known about the mechanisms underlying these outcomes. Chlamydia are genetically intractable as replication outside of the host cell is not yet possible and there are no practical tools for routine genetic manipulation, making genome-scale approaches critical. The early timeframe of infection is poorly understood and the host transcriptional response to chlamydial infection is not well defined. Our simultaneous RNA-Seq method was applied to a simplified in vitro model of chlamydial infection. We discovered a possible chlamydial strategy for early iron acquisition, putative immune dampening effects of chlamydial infection on the host cell, and present a hypothesis for Chlamydia-induced fibrotic scarring through runaway positive feedback loops. In general, simultaneous RNA-Seq helps to reveal the complex interplay between invading bacterial pathogens and their host mammalian cells and is immediately applicable to any bacteria/host cell interaction.

Visualization of comparative genomic analyses by BLAST score ratio.

BACKGROUND: The first microbial genome sequence, Haemophilus influenzae, was published in 1995. Since then, more than 400 microbial genome sequences have been completed or commenced. This massive influx of data provides the opportunity to obtain biological insights through comparative genomics. However few tools are available for this scale of comparative analysis. RESULTS: The BLAST Score Ratio (BSR) approach, implemented in a Perl script, classifies all putative peptides within three genomes using a measure of similarity based on the ratio of BLAST scores. The output of the BSR analysis enables global visualization of the degree of proteome similarity between all three genomes. Additional output enables the genomic synteny (conserved gene order) between each genome pair to be assessed. Furthermore, we extend this synteny analysis by overlaying BSR data as a color dimension, enabling visualization of the degree of similarity of the peptides being compared. CONCLUSIONS: Combining the degree of similarity, synteny and annotation will allow rapid identification of conserved genomic regions as well as a number of common genomic rearrangements such as insertions, deletions and inversions. The script and example visualizations are available at: http://www.microbialgenomics.org/BSR/.