The transcriptional landscape of Chlamydia pneumoniae

MLDSPP: Bacterial Promoter Prediction Tool Using DNA Structural Properties with Machine Learning and Explainable AI

Bacterial promoters play a crucial role in gene expression by serving as docking sites for the transcription initiation machinery. However, accurately identifying promoter regions in bacterial genomes remains a challenge due to their diverse architecture and variations. In this study, we propose MLDSPP (Machine Learning and Duplex Stability based Promoter prediction in Prokaryotes), a machine learning-based promoter prediction tool, to comprehensively screen bacterial promoter regions in 12 diverse genomes. We leveraged biologically relevant and informative DNA structural properties, such as DNA duplex stability and base stacking, and state-of-the-art machine learning (ML) strategies to gain insights into promoter characteristics. We evaluated several machine learning models, including Support Vector Machines, Random Forests, and XGBoost, and assessed their performance using accuracy, precision, recall, specificity, F1 score, and MCC metrics. Our findings reveal that XGBoost outperformed other models and current state-of-the-art promoter prediction tools, namely Sigma70pred and iPromoter2L, achieving F1-scores >95% in most systems. Significantly, the use of one-hot encoding for representing nucleotide sequences complements these structural features, enhancing our XGBoost model's predictive capabilities. To address the challenge of model interpretability, we incorporated explainable AI techniques using Shapley values. This enhancement allows for a better understanding and interpretation of the predictions of our model. In conclusion, our study presents MLDSPP as a novel, generic tool for predicting promoter regions in bacteria, utilizing original downstream sequences as nonpromoter controls. This tool has the potential to significantly advance the field of bacterial genomics and contribute to our understanding of gene regulation in diverse bacterial systems.

A Functional Yeast-Based Screen Identifies the Host Microtubule Cytoskeleton as a Target of Numerous Chlamydia pneumoniae Proteins

Bacterial pathogens have evolved intricate ways to manipulate the host to support infection. Here, we systematically assessed the importance of the microtubule cytoskeleton for infection by Chlamydiae, which are obligate intracellular bacteria that are of great importance for human health. The elimination of microtubules in human HEp-2 cells prior to C. pneumoniae infection profoundly attenuated the infection efficiency, demonstrating the need for microtubules for the early infection processes. To identify microtubule-modulating C. pneumoniae proteins, a screen in the model yeast Schizosaccharomyces pombe was performed. Unexpectedly, among 116 selected chlamydial proteins, more than 10%, namely, 13 proteins, massively altered the yeast interphase microtubule cytoskeleton. With two exceptions, these proteins were predicted to be inclusion membrane proteins. As proof of principle, we selected the conserved CPn0443 protein, which caused massive microtubule instability in yeast, for further analysis. CPn0443 bound and bundled microtubules in vitro and co-localized partially with microtubules in vivo in yeast and human cells. Furthermore, CPn0443-transfected U2OS cells had a significantly reduced infection rate by C. pneumoniae EBs. Thus, our yeast screen identified numerous proteins encoded using the highly reduced C. pneumoniae genome that modulated microtubule dynamics. Hijacking of the host microtubule cytoskeleton must be a vital part of chlamydial infection.

Molecular dynamics simulation-based trinucleotide and tetranucleotide level structural and energy characterization of the functional units of genomic DNA

Genomes of most organisms on earth are written in a universal language of life, made up of four units - adenine (A), thymine (T), guanine (G), and cytosine (C), and understanding the way they are put together has been a great challenge to date. Multiple efforts have been made to annotate this wonderfully engineered string of DNA using different methods but they lack a universal character. In this article, we have investigated the structural and energetic profiles of both prokaryotes and eukaryotes by considering two essential genomic sites, viz., the transcription start sites (TSS) and exon-intron boundaries. We have characterized these sites by mapping the structural and energy features of DNA obtained from molecular dynamics simulations, which considers all possible trinucleotide and tetranucleotide steps. For DNA, these physicochemical properties show distinct signatures at the TSS and intron-exon boundaries. Our results firmly convey the idea that DNA uses the same dialect for prokaryotes and eukaryotes and that it is worth going beyond sequence-level analyses to physicochemical space to determine the functional destiny of DNA sequences.

López Font F, A. Zayas Rodríguez F. Persistence in Chlamydia

Chlamydia spp. are important causes of acute and persistent/chronic infections. All Chlamydia spp. display a unique biphasic developmental cycle alternating between an infectious elementary body (EB) and a replicative form, the reticulate body (RB), followed by the multiplication of RBs by binary fission and progressive differentiation back into EBs. During its intracellular life, Chlamydia employs multiple mechanisms to ensure its persistence inside the host. These include evasion of diverse innate immune responses, modulation of host cell structure and endocytosis, inhibition of apoptosis, activation of pro-signaling pathways, and conversion to enlarged, non-replicative but viable “aberrant bodies” (ABs). Early research described several systems for Chlamydial persistence with a significant number of variables that make a direct comparison of results difficult. Now, emerging tools for genetic manipulations in Chlamydia and advances in global microarray, transcriptomics, and proteomics have opened new and exciting opportunities to understand the persistent state of Chlamydia and link the immune and molecular events of persistence with the pathogenesis of recurrent and chronic Chlamydial infections. This chapter reviews our current understanding and advances in the molecular biology of Chlamydia persistence.

Expression and structure of the Chlamydia trachomatis DksA ortholog

Chlamydia trachomatis is a bacterial obligate intracellular parasite and a significant cause of human disease, including sexually transmitted infections and trachoma. The bacterial RNA polymerase-binding protein DksA is a transcription factor integral to the multicomponent bacterial stress response pathway known as the stringent response. The genome of C. trachomatis encodes a DksA ortholog (DksA_Ct) that is maximally expressed at 15–20 h post infection, a time frame correlating with the onset of transition between the replicative reticulate body (RB) and infectious elementary body (EB) forms of the pathogen. Ectopic overexpression of DksA_Ct in C. trachomatis prior to RB–EB transitions during infection of HeLa cells resulted in a 39.3% reduction in overall replication (yield) and a 49.6% reduction in recovered EBs. While the overall domain organization of DksA_Ct is similar to the DksA ortholog of Escherichia coli (DksA_Ec), DksA_Ct did not functionally complement DksA_Ec. Transcription of dksA_Ct is regulated by tandem promoters, one of which also controls expression of nrdR, encoding a negative regulator of deoxyribonucleotide biosynthesis. The phenotype resulting from ectopic expression of DksA_Ct and the correlation between dksA_Ct and nrdR expression is consistent with a role for DksA_Ct in the C. trachomatis developmental cycle.

Promotech: a general tool for bacterial promoter recognition

Promoters are genomic regions where the transcription machinery binds to initiate the transcription of specific genes. Computational tools for identifying bacterial promoters have been around for decades. However, most of these tools were designed to recognize promoters in one or few bacterial species. Here, we present Promotech, a machine-learning-based method for promoter recognition in a wide range of bacterial species. We compare Promotech's performance with the performance of five other promoter prediction methods. Promotech outperforms these other programs in terms of area under the precision-recall curve (AUPRC) or precision at the same level of recall. Promotech is available at https://github.com/BioinformaticsLabAtMUN/PromoTech .

Cross-species RNA-seq for deciphering host-microbe interactions

The human body is constantly exposed to microorganisms, which entails manifold interactions between human cells and diverse commensal or pathogenic bacteria. The cellular states of the interacting cells are decisive for the outcome of these encounters such as whether bacterial virulence programmes and host defence or tolerance mechanisms are induced. This Review summarizes how next-generation RNA sequencing (RNA-seq) has become a primary technology to study host-microbe interactions with high resolution, improving our understanding of the physiological consequences and the mechanisms at play. We illustrate how the discriminatory power and sensitivity of RNA-seq helps to dissect increasingly complex cellular interactions in time and space down to the single-cell level. We also outline how future transcriptomics may answer currently open questions in host-microbe interactions and inform treatment schemes for microbial disorders.

Acquisition of pcnB [poly(A) polymerase I] genes via horizontal transfer from the β, γ-Proteobacteria

Poly(A) polymerases (PAPs) and tRNA nucleotidyltransferases belong to a superfamily of nucleotidyltransferases and modify RNA 3'-ends. The product of the pcnB gene, PAP I, has been characterized in a few β-, γ- and δ-Proteobacteria. Using the PAP I signature sequence, putative PAPs were identified in bacterial species from the α- and ε-Proteobacteria and from four other bacterial phyla (Firmicutes, Actinobacteria, Bacteroidetes and Aquificae). Phylogenetic analysis, alien index and G+C content calculations strongly suggest that the PAPs in the species identified in this study arose by horizontal gene transfer from the β- and γ-Proteobacteria.

Transcriptional Landscape of Waddlia chondrophila Aberrant Bodies Induced by Iron Starvation

Chronic infections caused by obligate intracellular bacteria belonging to the Chlamydiales order are related to the formation of persistent developmental forms called aberrant bodies (ABs), which undergo DNA replication without cell division. These enlarged bacteria develop and persist upon exposure to different stressful conditions such as β-lactam antibiotics, iron deprivation and interferon-γ. However, the mechanisms behind ABs biogenesis remain uncharted. Using an RNA-sequencing approach, we compared the transcriptional profile of ABs induced by iron starvation to untreated bacteria in the Chlamydia-related species Waddliachondrophila, a potential agent of abortion in ruminants and miscarriage in humans. Consistent with the growth arrest observed following iron depletion, our results indicate a significant reduction in the expression of genes related to energy production, carbohydrate and amino acid metabolism and cell wall/envelope biogenesis, compared to untreated, actively replicating bacteria. Conversely, three putative toxin-antitoxin modules were among the most up-regulated genes upon iron starvation, suggesting that their activation might be involved in growth arrest in adverse conditions, an uncommon feature in obligate intracellular bacteria. Our work represents the first complete transcriptomic profile of a Chlamydia-related species in stressful conditions and sets the grounds for further investigations on the mechanisms underlying chlamydial persistence.

A novel method SEProm for prokaryotic promoter prediction based on DNA structure and energetics

MOTIVATION

Despite conservation in general architecture of promoters and protein-DNA interaction interface of RNA polymerases among various prokaryotes, identification of promoter regions in the whole genome sequences remains a daunting challenge. The available tools for promoter prediction do not seem to address the problem satisfactorily, apparently because the biochemical nature of promoter signals is yet to be understood fully. Using 28 structural and 3 energetic parameters, we found that prokaryotic promoter regions have a unique structural and energy state, quite distinct from that of coding regions and the information for this signature state is in-built in their sequences. We developed a novel promoter prediction tool from these 31 parameters using various statistical techniques.

RESULTS

Here, we introduce SEProm, a novel tool that is developed by studying and utilizing the in-built structural and energy information of DNA sequences, which is applicable to all prokaryotes including archaea. Compared to five most recent, diverged and current best available tools, SEProm performs much better, predicting promoters with an 'F-value' of 82.04 and 'Precision' of 81.08. The next best 'F-value' was obtained with PromPredict (72.14) followed by BProm (68.37). On the basis of 'Precision' value, the next best 'Precision' was observed for Pepper (75.39) followed by PromPredict (72.01). SEProm maintained the lead even when comparison was done on two test organisms (not involved in training for SEProm).

AVAILABILITY AND IMPLEMENTATION

The software is freely available with easy to follow instructions (www.scfbio-iitd.res.in/software/TSS_Predict.jsp).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Whole-genome Identification of Transcriptional Start Sites by Differential RNA-seq in Bacteria

Gene transcription in bacteria often starts some nucleotides upstream of the start codon. Identifying the specific Transcriptional Start Site (TSS) is essential for genetic manipulation, as in many cases upstream of the start codon there are sequence elements that are involved in gene expression regulation. Taken into account the classical gene structure, we are able to identify two kinds of transcriptional start site: primary and secondary. A primary transcriptional start site is located some nucleotides upstream of the translational start site, while a secondary transcriptional start site is located within the gene encoding sequence. Here, we present a step by step protocol for genome-wide transcriptional start sites determination by differential RNA-sequencing (dRNA-seq) using the enteric pathogen Shigella flexneri serotype 5a strain M90T as model. However, this method can be employed in any other bacterial species of choice. In the first steps, total RNA is purified from bacterial cultures using the hot phenol method. Ribosomal RNA (rRNA) is specifically depleted via hybridization probes using a commercial kit. A 5'-monophosphate-dependent exonuclease (TEX)-treated RNA library enriched in primary transcripts is then prepared for comparison with a library that has not undergone TEX-treatment, followed by ligation of an RNA linker adaptor of known sequence allowing the determination of TSS with single nucleotide precision. Finally, the RNA is processed for Illumina sequencing library preparation and sequenced as purchased service. TSS are identified by in-house bioinformatic analysis. Our protocol is cost-effective as it minimizes the use of commercial kits and employs freely available software.

A high-resolution transcriptome map identifies small RNA regulation of metabolism in the gut microbe Bacteroides thetaiotaomicron

Bacteria of the genus Bacteroides are common members of the human intestinal microbiota and important degraders of polysaccharides in the gut. Among them, the species Bacteroides thetaiotaomicron has emerged as the model organism for functional microbiota research. Here, we use differential RNA sequencing (dRNA-seq) to generate a single-nucleotide resolution transcriptome map of B. thetaiotaomicron grown under defined laboratory conditions. An online browser, called ‘Theta-Base’ (www.helmholtz-hiri.de/en/datasets/bacteroides), is launched to interrogate the obtained gene expression data and annotations of ~4500 transcription start sites, untranslated regions, operon structures, and 269 noncoding RNA elements. Among the latter is GibS, a conserved, 145 nt-long small RNA that is highly expressed in the presence of N-acetyl-D-glucosamine as sole carbon source. We use computational predictions and experimental data to determine the secondary structure of GibS and identify its target genes. Our results indicate that sensing of N-acetyl-D-glucosamine induces GibS expression, which in turn modifies the transcript levels of metabolic enzymes.

Bacteroides thetaiotaomicron is a human gut microbe and an emergent model organism. Here, Ryan et al. generate single-nucleotide resolution RNA-seq data for this bacterium and map transcription start sites and noncoding RNAs, one of which modulates expression of metabolic enzymes.

Dual RNA-seq of Orientia tsutsugamushi informs on host-pathogen interactions for this neglected intracellular human pathogen

Studying emerging or neglected pathogens is often challenging due to insufficient information and absence of genetic tools. Dual RNA-seq provides insights into host-pathogen interactions, and is particularly informative for intracellular organisms. Here we apply dual RNA-seq to Orientia tsutsugamushi (Ot), an obligate intracellular bacterium that causes the vector-borne human disease scrub typhus. Half the Ot genome is composed of repetitive DNA, and there is minimal collinearity in gene order between strains. Integrating RNA-seq, comparative genomics, proteomics, and machine learning to study the transcriptional architecture of Ot, we find evidence for wide-spread post-transcriptional antisense regulation. Comparing the host response to two clinical isolates, we identify distinct immune response networks for each strain, leading to predictions of relative virulence that are validated in a mouse infection model. Thus, dual RNA-seq can provide insight into the biology and host-pathogen interactions of a poorly characterized and genetically intractable organism such as Ot.

Complete genome sequence and annotation of the laboratory reference strain Shigella flexneri serotype 5a M90T and genome-wide transcriptional start site determination

Background

Shigella is a Gram-negative facultative intracellular bacterium that causes bacillary dysentery in humans. Shigella invades cells of the colonic mucosa owing to its virulence plasmid-encoded Type 3 Secretion System (T3SS), and multiplies in the target cell cytosol. Although the laboratory reference strain S. flexneri serotype 5a M90T has been extensively used to understand the molecular mechanisms of pathogenesis, its complete genome sequence is not available, thereby greatly limiting studies employing high-throughput sequencing and systems biology approaches.

Results

We have sequenced, assembled, annotated and manually curated the full genome of S. flexneri 5a M90T. This yielded two complete circular contigs, the chromosome and the virulence plasmid (pWR100). To obtain the genome sequence, we have employed long-read PacBio DNA sequencing followed by polishing with Illumina RNA-seq data. This provides a new hybrid strategy to prepare gapless, highly accurate genome sequences, which also cover AT-rich tracks or repetitive sequences that are transcribed. Furthermore, we have performed genome-wide analysis of transcriptional start sites (TSS) and determined the length of 5′ untranslated regions (5′-UTRs) at typical culture conditions for the inoculum of in vitro infection experiments. We identified 6723 primary TSS (pTSS) and 7328 secondary TSS (sTSS). The S. flexneri 5a M90T annotated genome sequence and the transcriptional start sites are integrated into RegulonDB (http://regulondb.ccg.unam.mx) and RSAT (http://embnet.ccg.unam.mx/rsat/) databases to use their analysis tools in the S. flexneri 5a M90T genome.

Conclusions

We provide the first complete genome for S. flexneri serotype 5a, specifically the laboratory reference strain M90T. Our work opens the possibility of employing S. flexneri M90T in high-quality systems biology studies such as transcriptomic and differential expression analyses or in genome evolution studies. Moreover, the catalogue of TSS that we report here can be used in molecular pathogenesis studies as a resource to know which genes are transcribed before infection of host cells. The genome sequence, together with the analysis of transcriptional start sites, is also a valuable tool for precise genetic manipulation of S. flexneri 5a M90T. Further, we present a new hybrid strategy to prepare gapless, highly accurate genome sequences. Unlike currently used hybrid strategies combining long- and short-read DNA sequencing technologies to maximize accuracy, our workflow using long-read DNA sequencing and short-read RNA sequencing provides the added value of using non-redundant technologies, which yield distinct, exploitable datasets.

Structure of a Core Promoter in Bifidobacterium longum NCC2705

Bacterial promoters consist of core sequence motifs termed -35 and -10 boxes. The consensus motifs are TTGACA and TATAAT, respectively, which were identified from leading investigations on Escherichia coli However, the consensus sequences are not likely to fit genetically divergent bacteria. The sigma factor of the genus Bifidobacterium has a characteristic polar domain in the N terminus, suggesting the possibility of specific promoter recognition. We reevaluated the structure of Bifidobacterium longum NCC2705 promoters and compared them to other bacteria. Transcriptional start sites (TSSs) of the B. longum NCC2705 strain were identified using transcriptome sequencing (RNA-Seq) analysis to extract promoter regions. Conserved motifs of a bifidobacterial promoter were determined using regions upstream of TSSs and a hidden Markov model. As a result, consensus motifs of the -35 and -10 boxes were TTGTGC and TACAAT, respectively. To assess each base of both motifs, we constructed 37 plasmids based on pKO403-TPCTcon, including the hup promoter connected with a chloramphenicol acetyltransferase as a reporter gene. This reporter assay showed two optimal motifs of the -35 and -10 boxes, namely, TTGNNN and TANNNT, respectively. We further analyzed spacer lengths between the -35 and -10 boxes via a bioinformatics approach. The spacer lengths predominant in bacteria have been generally reported to be approximately 17 bp. In contrast, the predominant spacer lengths in the genus Bifidobacterium and related species were 11 bp, in addition to 17 bp. A reporter assay to assess the spacer lengths indicated that the 11-bp spacer length produced unusually high activity.IMPORTANCE The structures of sigma factors vary among bacterial strains, indicating that recognition rules may also vary. Therefore, we investigated the promoter structure of Bifidobacterium longum NCC2705 using a bioinformatics approach and wet analyses. The most frequent and optimal motifs were similar to other bacterial consensus motifs. The optimal spacer length between the two boxes was reported to be 17 bp. It is widely applied to a bioinformatics approach for other bacteria. Unexpectedly, conserved spacer lengths were 11 bp as well as 17 bp in the genus Bifidobacterium Moreover, the sigma factor of the genus Bifidobacterium has a characteristic domain in the N terminus which may contribute to the additional functions. Hence, it would be valuable to reevaluate the promoter in other organisms.

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.

NK Cell-Mediated Processing Of Chlamydia psittaci Drives Potent Anti-Bacterial Th1 Immunity

Natural killer (NK) cells are innate immune cells critically involved in the early immune response against various pathogens including chlamydia. Here, we demonstrate that chlamydia-infected NK cells prevent the intracellular establishment and growth of the bacteria. Upon infection, they display functional maturation characterized by enhanced IFN-γ secretion, CD146 induction, PKCϴ activation, and granule secretion. Eventually, chlamydia are released in a non-infectious, highly immunogenic form driving a potent Th1 immune response. Further, anti-chlamydial antibodies generated during immunization neutralize the infection of epithelial cells. The release of chlamydia from NK cells requires PKCϴ function and active degranulation, while granule-associated granzyme B drives the loss of chlamydial infectivity. Cellular infection and bacterial release can be undergone repeatedly and do not affect NK cell function. Strikingly, NK cells passing through such an infection cycle significantly improve their cytotoxicity. Thus, NK cells not only protect themselves against productive chlamydial infections but also actively trigger potent anti-bacterial responses.

Regulatory RNAs in Virulence and Host-Microbe Interactions

Bacterial regulatory RNAs are key players in adaptation to changing environmental conditions and response to diverse cellular stresses. However, while regulatory RNAs of bacterial pathogens have been intensely studied under defined conditions in vitro, characterization of their role during the infection of eukaryotic host organisms is lagging behind. This review summarizes our current understanding of the contribution of the different classes of regulatory RNAs and RNA-binding proteins to bacterial virulence and illustrates their role in infection by reviewing the mechanisms of some prominent representatives of each class. Emerging technologies are described that bear great potential for global, unbiased studies of virulence-related RNAs in bacterial model and nonmodel pathogens in the future. The review concludes by deducing common principles of RNA-mediated gene expression control of virulence programs in different pathogens, and by defining important open questions for upcoming research in the field.

Floxed-Cassette Allelic Exchange Mutagenesis Enables Markerless Gene Deletion in Chlamydia trachomatis and Can Reverse Cassette-Induced Polar Effects

As obligate intracellular bacteria, Chlamydia spp. have evolved numerous, likely intricate, mechanisms to create and maintain a privileged intracellular niche. Recent progress in elucidating and characterizing these processes has been bolstered by the development of techniques enabling basic genetic tractability. Florescence-reported allelic exchange mutagenesis (FRAEM) couples chromosomal gene deletion with the insertion of a selection cassette encoding antibiotic resistance and green fluorescent protein (GFP). Similar to other bacteria, many chlamydial genes exist within polycistronic operons, raising the possibility of polar effects mediated by insertion cassettes. Indeed, FRAEM-mediated deletion of Chlamydia trachomatis tmeA negatively impacts the expression of tmeB We have adapted FRAEM technology by employing a gfp-bla cassette flanked by loxP sites. Conditional expression of Cre recombinase in Chlamydia tmeA containing a floxed cassette resulted in deletion of the marker and restoration of tmeB expression.IMPORTANCE C. trachomatis infections represent a significant burden to human health. The ability to genetically manipulate Chlamydia spp. is overcoming historic confounding barriers that have impeded rapid progress in understanding overall chlamydial pathogenesis. The current state of genetic manipulation in Chlamydia spp. requires further development, including mechanisms to generate markerless gene disruption. We leveraged a stepwise Cre-lox approach to excise selection marker genes from a deleted gene locus. We found this process to be efficient, and the removal of extraneous elements resulted in the reversal of a negative polar effect on a downstream gene. This technique facilitates a more direct assessment of gene function and adds to the Chlamydia molecular toolbox by facilitating the deletion of genes within operons.

Toward a Universal Structural and Energetic Model for Prokaryotic Promoters

With almost no consensus promoter sequence in prokaryotes, recruitment of RNA polymerase (RNAP) to precise transcriptional start sites (TSSs) has remained an unsolved puzzle. Uncovering the underlying mechanism is critical for understanding the principle of gene regulation. We attempted to search the hidden code in ∼16,500 promoters of 12 prokaryotes representing two kingdoms in their structure and energetics. Twenty-eight fundamental parameters of DNA structure including backbone angles, basepair axis, and interbasepair and intrabasepair parameters were used, and information was extracted from x-ray crystallography data. Three parameters (solvation energy, hydrogen-bond energy, and stacking energy) were selected for creating energetics profiles using in-house programs. DNA of promoter regions was found to be inherently designed to undergo a change in every parameter undertaken for the study, in all prokaryotes. The change starts from some distance upstream of TSSs and continues past some distance from TSS, hence giving a signature state to promoter regions. These signature states might be the universal hidden codes recognized by RNAP. This observation was reiterated when randomly selected promoter sequences (with little sequence conservation) were subjected to structure generation; all developed into very similar three-dimensional structures quite distinct from those of conventional B-DNA and coding sequences. Fine structural details at important motifs (viz. -11, -35, and -75 positions relative to TSS) of promoters reveal novel to our knowledge and pointed insights for RNAP interaction at these locations; it could be correlated with how some particular structural changes at the -11 region may allow insertion of RNAP amino acids in interbasepair space as well as facilitate the flipping out of bases from the DNA duplex.

Virulence-related comparative transcriptomics of infectious and non-infectious chlamydial particles

BACKGROUND

Members of the phylum Chlamydiae are obligate intracellular pathogens of humans and animals and have a serious impact on host health. They comprise several zoonotic species with varying disease outcomes and prevalence. To investigate differences in virulence, we focused on Chlamydia psittaci, C. abortus and Waddlia chondrophila. Most threatening is C. psittaci, which frequently infects humans and causes psittacosis associated with severe pneumonia. The closest relative of C. psittaci is C. abortus, which shares the vast majority of genes but less frequently infects humans, and causes stillbirth and sepsis. W. chondrophila is more distantly related, and occasional human infections are associated with respiratory diseases or miscarriage. One possible explanation for differences in virulence originate from species-specific genes as well as differentially expressed homologous virulence factors.

RESULTS

RNA-sequencing (RNA-Seq) was applied to purified infectious elementary bodies (EBs) and non-infectious reticulate bodies (RBs) in order to elucidate the transcriptome of the infectious and replicative chlamydial states. The results showed that approximately half of all genes were differentially expressed. For a descriptive comparison, genes were categorised according to their function in the RAST database. This list was extended by the inclusion of inclusion membrane proteins, outer membrane proteins, polymorphic membrane proteins and type III secretion system effectors. In addition, the expression of fifty-six known and a variety of predicted virulence and immunogenic factors with homologs in C. psittaci, C. abortus and W. chondrophila was analysed. To confirm the RNA-Seq results, the expression of nine factors was validated using real-time quantitative polymerase chain reaction (RT-qPCR). Comparison of RNA-Seq and RT-qPCR results showed a high mean Pearson correlation coefficient of 0.95.

CONCLUSIONS

It was shown that both the replicative and infectious chlamydial state contained distinctive transcriptomes and the cellular processes emphasised in EBs and RBs differed substantially based on the chlamydial species. In addition, the very first interspecies transcriptome comparison is presented here, and the considerable differences in expression of homologous virulence factors might contribute to the differing infection rates and disease outcomes of the pathogens. The RNA-Seq results were confirmed by RT-qPCR and demonstrate the feasibility of interspecies transcriptome comparisons in chlamydia.

WikiGenomes: an open web application for community consumption and curation of gene annotation data in Wikidata

With the advancement of genome-sequencing technologies, new genomes are being sequenced daily. Although these sequences are deposited in publicly available data warehouses, their functional and genomic annotations (beyond genes which are predicted automatically) mostly reside in the text of primary publications. Professional curators are hard at work extracting those annotations from the literature for the most studied organisms and depositing them in structured databases. However, the resources don’t exist to fund the comprehensive curation of the thousands of newly sequenced organisms in this manner. Here, we describe WikiGenomes (wikigenomes.org), a web application that facilitates the consumption and curation of genomic data by the entire scientific community. WikiGenomes is based on Wikidata, an openly editable knowledge graph with the goal of aggregating published knowledge into a free and open database. WikiGenomes empowers the individual genomic researcher to contribute their expertise to the curation effort and integrates the knowledge into Wikidata, enabling it to be accessed by anyone without restriction.

Database URL: www.wikigenomes.org

A genetic system for targeted mutations to disrupt and restore genes in the obligate bacterium, Ehrlichia chaffeensis

Obligate intracellular bacteria (obligates) belonging to Rickettsiales and Chlamydiales cause diseases in hundreds of millions of people worldwide and in many animal species. Lack of an efficient system for targeted mutagenesis in obligates remains a major impediment in understanding microbial pathogenesis. Challenges in creating targeted mutations may be attributed to essential nature of majority of the genes and intracellular replication dependence. Despite success in generating random mutations, a method that works well in creating mutations in specific genes of interest followed by complementation remains problematic for obligates and is a highly sought-after goal. We describe protocols to generate stable targeted mutations by allelic exchange in Ehrlichia chaffeensis, an obligate intracellular tick-borne bacterium responsible for human monocytic ehrlichiosis. Targeted mutations in E. chaffeensis were created to disrupt two genes, and also to restore one gene by another allelic exchange mutation leading to the restoration of transcription and protein expression from the inactivated gene and the recovered organisms also express mCherry, which distinguishes from the wild type. We expect that the methods developed are broadly applicable to other obligates, particularly to rickettsial pathogens, to routinely perform targeted mutations to enable studies focused on protein structure-function analyses, host-pathogen interactions and in developing vaccines.

Small RNAs in Bacterial Virulence and Communication

Bacterial pathogens must endure or adapt to different environments and stresses during transmission and infection. Posttranscriptional gene expression control by regulatory RNAs, such as small RNAs and riboswitches, is now considered central to adaptation in many bacteria, including pathogens. The study of RNA-based regulation (riboregulation) in pathogenic species has provided novel insight into how these bacteria regulate virulence gene expression. It has also uncovered diverse mechanisms by which bacterial small RNAs, in general, globally control gene expression. Riboregulators as well as their targets may also prove to be alternative targets or provide new strategies for antimicrobials. In this article, we present an overview of the general mechanisms that bacteria use to regulate with RNA, focusing on examples from pathogens. In addition, we also briefly review how deep sequencing approaches have aided in opening new perspectives in small RNA identification and the study of their functions. Finally, we discuss examples of riboregulators in two model pathogens that control virulence factor expression or survival-associated phenotypes, such as stress tolerance, biofilm formation, or cell-cell communication, to illustrate how riboregulation factors into regulatory networks in bacterial pathogens.

Biphasic Metabolism and Host Interaction of a Chlamydial Symbiont

Chlamydiae are obligate intracellular bacteria comprising well-known human pathogens and ubiquitous symbionts of protists, which are characterized by a unique developmental cycle. Here we comprehensively analyzed gene expression dynamics of Protochlamydia amoebophila during infection of its Acanthamoeba host by RNA sequencing. This revealed a highly dynamic transcriptional landscape, where major transcriptional shifts are conserved among chlamydial symbionts and pathogens. Our data served to propose a time-resolved model for type III protein secretion during the developmental cycle, and we provide evidence for a biphasic metabolism of P. amoebophila during infection, which involves energy parasitism and amino acids as the carbon source during initial stages and a postreplicative switch to endogenous glucose-based ATP production. This fits well with major transcriptional changes in the amoeba host, where upregulation of complex sugar breakdown precedes the P. amoebophila metabolic switch. The biphasic chlamydial metabolism represents a unique adaptation to exploit eukaryotic host cells, which likely contributed to the evolutionary success of this group of microbes. IMPORTANCE Chlamydiae are known as major bacterial pathogens of humans, causing the ancient disease trachoma, but they are also frequently found in the environment where they infect ubiquitous protists such as amoebae. All known chlamydiae require a eukaryotic host cell to thrive. Using the environmental chlamydia Protochlamydia amoebophila within its natural host, Acanthamoeba castellanii, we investigated gene expression dynamics in vivo and throughout the complete chlamydial developmental cycle for the first time. This allowed us to infer how a major virulence mechanism, the type III secretion system, is regulated and employed, and we show that the physiology of chlamydiae undergoes a complete shift regarding carbon metabolism and energy generation. This study provides comprehensive insights into the infection strategy of chlamydiae and reveals a unique adaptation to life within a eukaryotic host cell.

If it transcribes, we can sequence it: mining the complexities of host-pathogen-environment interactions using RNA-seq

Host-pathogen interactions are exceedingly complex because they involve multiple host tissues, often occur in the context of normal microflora, and can span diverse microenvironments. Although decades of gene expression studies have provided detailed insights into infection processes, technical challenges have restricted experiments to single pathogenic species or host tissues. RNA-sequencing (RNA-seq) has revolutionized the study of gene expression because in addition to quantifying transcriptional output, it allows detection and characterization of all transcripts in a genome. Here, we review how refined approaches to RNA-seq are used to map the transcriptional networks that control host-pathogen interactions. These enhanced techniques include dRNA-seq and term-seq for the fine-scale mapping of transcriptional start and termination sites, and dual RNA-seq for simultaneous sequencing of host and bacterial pathogen transcriptomes. Dual RNA-seq experiments are currently limited to in vitro infection systems that do not fully reflect the complexities of the in vivo environment, thus a challenge is to develop in vivo model systems and experimental approaches that address the biological heterogeneity of host environments, followed by the integration of RNA-seq with other genome-scale datasets to identify the transcriptional networks that mediate host-pathogen interactions.

Evolutionary Cell Biology of Division Mode in the Bacterial Planctomycetes-Verrucomicrobia- Chlamydiae Superphylum

Bacteria from the Planctomycetes, Verrucomicrobia, and Chlamydiae (PVC) superphylum are exceptions to the otherwise dominant mode of division by binary fission, which is based on the interaction between the FtsZ protein and the peptidoglycan (PG) biosynthesis machinery. Some PVC bacteria are deprived of the FtsZ protein and were also thought to lack PG. How these bacteria divide is still one of the major mysteries of microbiology. The presence of PG has recently been revealed in Planctomycetes and Chlamydiae, and proteins related to PG synthesis have been shown to be implicated in the division process in Chlamydiae, providing important insights into PVC mechanisms of division. Here, we review the historical lack of observation of PG in PVC bacteria, its recent detection in two phyla and its involvement in chlamydial cell division. Based on the detection of PG-related proteins in PVC proteomes, we consider the possible evolution of the diverse division mechanisms in these bacteria. We conclude by summarizing what is known and what remains to be understood about the evolutionary cell biology of PVC division modes.

Interrogating Genes That Mediate Chlamydia trachomatis Survival in Cell Culture Using Conditional Mutants and Recombination

Intracellular bacterial pathogens in the family Chlamydiaceae are causes of human blindness, sexually transmitted disease, and pneumonia. Genetic dissection of the mechanisms of chlamydial pathogenicity has been hindered by multiple limitations, including the inability to inactivate genes that would prevent the production of elementary bodies. Many genes are also Chlamydia-specific genes, and chlamydial genomes have undergone extensive reductive evolution, so functions often cannot be inferred from homologs in other organisms. Conditional mutants have been used to study essential genes of many microorganisms, so we screened a library of 4,184 ethyl methanesulfonate-mutagenized Chlamydia trachomatis isolates for temperature-sensitive (TS) mutants that developed normally at physiological temperature (37°C) but not at nonphysiological temperatures. Heat-sensitive TS mutants were identified at a high frequency, while cold-sensitive mutants were less common. Twelve TS mutants were mapped using a novel markerless recombination approach, PCR, and genome sequencing. TS alleles of genes that play essential roles in other bacteria and chlamydia-specific open reading frames (ORFs) of unknown function were identified. Temperature-shift assays determined that phenotypes of the mutants manifested at distinct points in the developmental cycle. Genome sequencing of a larger population of TS mutants also revealed that the screen had not reached saturation. In summary, we describe the first approach for studying essential chlamydial genes and broadly applicable strategies for genetic mapping in Chlamydia spp. and mutants that both define checkpoints and provide insights into the biology of the chlamydial developmental cycle.

IMPORTANCE

Study of the pathogenesis of Chlamydia spp. has historically been hampered by a lack of genetic tools. Although there has been recent progress in chlamydial genetics, the existing approaches have limitations for the study of the genes that mediate growth of these organisms in cell culture. We used a genetic screen to identify conditional Chlamydia mutants and then mapped these alleles using a broadly applicable recombination strategy. Phenotypes of the mutants provide fundamental insights into unexplored areas of chlamydial pathogenesis and intracellular biology. Finally, the reagents and approaches we describe are powerful resources for the investigation of these organisms.

Purification of infectious and non-infectious chlamydial particles using iodixanol for density gradient preparation

Chlamydiae are obligate intracellular bacteria with two distinct morphological stages, the infectious elementary bodies (EBs) and non-infectious reticulate bodies (RBs). Here we describe a rapid and straightforward protocol for the purification of EBs and RBs involving special density gradients. It has been successfully applied to three chlamydial species.

Sequestration of host metabolism by an intracellular pathogen

For intracellular pathogens, residence in a vacuole provides a shelter against cytosolic host defense to the cost of limited access to nutrients. The human pathogen Chlamydia trachomatis grows in a glycogen-rich vacuole. How this large polymer accumulates there is unknown. We reveal that host glycogen stores shift to the vacuole through two pathways: bulk uptake from the cytoplasmic pool, and de novo synthesis. We provide evidence that bacterial glycogen metabolism enzymes are secreted into the vacuole lumen through type 3 secretion. Our data bring strong support to the following scenario: bacteria co-opt the host transporter SLC35D2 to import UDP-glucose into the vacuole, where it serves as substrate for de novo glycogen synthesis, through a remarkable adaptation of the bacterial glycogen synthase. Based on these findings we propose that parasitophorous vacuoles not only offer protection but also provide a microorganism-controlled metabolically active compartment essential for redirecting host resources to the pathogens.

Regulatory (pan-)genome of an obligate intracellular pathogen in the PVC superphylum

Like other obligate intracellular bacteria, the Chlamydiae feature a compact regulatory genome that remains uncharted owing to poor genetic tractability. Exploiting the reduced number of transcription factors (TFs) encoded in the chlamydial (pan-)genome as a model for TF control supporting the intracellular lifestyle, we determined the conserved landscape of TF specificities by ChIP-Seq (chromatin immunoprecipitation-sequencing) in the chlamydial pathogen Waddlia chondrophila. Among 10 conserved TFs, Euo emerged as a master TF targeting >100 promoters through conserved residues in a DNA excisionase-like winged helix-turn-helix-like (wHTH) fold. Minimal target (Euo) boxes were found in conserved developmentally-regulated genes governing vertical genome transmission (cytokinesis and DNA replication) and genome plasticity (transposases). Our ChIP-Seq analysis with intracellular bacteria not only reveals that global TF regulation is maintained in the reduced regulatory genomes of Chlamydiae, but also predicts that master TFs interpret genomic information in the obligate intracellular α-proteobacteria, including the rickettsiae, from which modern day mitochondria evolved.

One Face of Chlamydia trachomatis: The Infectious Elementary Body

The lifestyle of Chlamydiae is unique: the bacteria alternate between two morphologically distinct forms, an infectious non-replicative elementary body (EB), and a replicative, non-infectious reticulate body (RB). This review focuses on recent advances in understanding the structure and function of the infectious form of the best-studied member of the phylum, the human pathogen Chlamydia trachomatis. Once considered as an inert particle of little functional capacity, the EB is now perceived as a sophisticated entity that encounters at least three different environments during each infectious cycle. We review current knowledge on its composition and morphology, and emerging metabolic activities. These features confer resistance to the extracellular environment, the ability to penetrate a host cell and ultimately enable the EB to establish a niche enabling bacterial survival and growth. The bacterial and host molecules involved in these processes are beginning to emerge.

A Coming of Age Story: Chlamydia in the Post-Genetic Era

Chlamydia spp. are ubiquitous, obligate, intracellular Gram-negative bacterial pathogens that undergo a unique biphasic developmental cycle transitioning between the infectious, extracellular elementary body and the replicative, intracellular reticulate body. The primary Chlamydia species associated with human disease are C. trachomatis, which is the leading cause of both reportable bacterial sexually transmitted infections and preventable blindness, and C. pneumoniae, which infects the respiratory tract and is associated with cardiovascular disease. Collectively, these pathogens are a significant source of morbidity and pose a substantial financial burden on the global economy. Past efforts to elucidate virulence mechanisms of these unique and important pathogens were largely hindered by an absence of genetic methods. Watershed studies in 2011 and 2012 demonstrated that forward and reverse genetic approaches were feasible with Chlamydia and that shuttle vectors could be selected and maintained within the bacterium. While these breakthroughs have led to a steady expansion of the chlamydial genetic tool kit, there are still roads left to be traveled. This minireview provides a synopsis of the currently available genetic methods for Chlamydia along with a comparison to the methods used in other obligate intracellular bacteria. Limitations and advantages of these techniques will be discussed with an eye toward the methods still needed, and how the current state of the art for genetics in obligate intracellular bacteria could direct future technological advances for Chlamydia.

dRNA-Seq Reveals Genomewide TSSs and Noncoding RNAs of Plant Beneficial Rhizobacterium Bacillus amyloliquefaciens FZB42

Bacillus amyloliquefaciens subsp. plantarum FZB42 is a representative of Gram-positive plant-growth-promoting rhizobacteria (PGPR) that inhabit plant root environments. In order to better understand the molecular mechanisms of bacteria-plant symbiosis, we have systematically analyzed the primary transcriptome of strain FZB42 grown under rhizosphere-mimicking conditions using differential RNA sequencing (dRNA-seq). Our analysis revealed 4,877 transcription start sites for protein-coding genes, identified genes differentially expressed under different growth conditions, and corrected many previously mis-annotated genes. We also identified a large number of riboswitches and cis-encoded antisense RNAs, as well as trans-encoded small noncoding RNAs that may play important roles in the gene regulation of Bacillus. Overall, our analyses provided a landscape of Bacillus primary transcriptome and improved the knowledge of rhizobacteria-host interactions.

Following the Footsteps of Chlamydial Gene Regulation

Regulation of gene expression ensures an organism responds to stimuli and undergoes proper development. Although the regulatory networks in bacteria have been investigated in model microorganisms, nearly nothing is known about the evolution and plasticity of these networks in obligate, intracellular bacteria. The phylum Chlamydiae contains a vast array of host-associated microbes, including several human pathogens. The Chlamydiae are unique among obligate, intracellular bacteria as they undergo a complex biphasic developmental cycle in which large swaths of genes are temporally regulated. Coupled with the low number of transcription factors, these organisms offer a model to study the evolution of regulatory networks in intracellular organisms. We provide the first comprehensive analysis exploring the diversity and evolution of regulatory networks across the phylum. We utilized a comparative genomics approach to construct predicted coregulatory networks, which unveiled genus- and family-specific regulatory motifs and architectures, most notably those of virulence-associated genes. Surprisingly, our analysis suggests that few regulatory components are conserved across the phylum, and those that are conserved are involved in the exploitation of the intracellular niche. Our study thus lends insight into a component of chlamydial evolution that has otherwise remained largely unexplored.

Discovery of putative small non-coding RNAs from the obligate intracellular bacterium Wolbachia pipientis

Wolbachia pipientis is an endosymbiotic bacterium that induces a wide range of effects in its insect hosts, including manipulation of reproduction and protection against pathogens. Little is known of the molecular mechanisms underlying the insect-Wolbachia interaction, though it is likely to be mediated via the secretion of proteins or other factors. There is an increasing amount of evidence that bacteria regulate many cellular processes, including secretion of virulence factors, using small non-coding RNAs (sRNAs), but sRNAs have not previously been described from Wolbachia. We have used two independent approaches, one based on comparative genomics and the other using RNA-Seq data generated for gene expression studies, to identify candidate sRNAs in Wolbachia. We experimentally characterized the expression of one of these candidates in four Wolbachia strains, and showed that it is differentially regulated in different host tissues and sexes. Given the roles played by sRNAs in other host-associated bacteria, the conservation of the candidate sRNAs between different Wolbachia strains, and the sex- and tissue-specific differential regulation we have identified, we hypothesise that sRNAs may play a significant role in the biology of Wolbachia, and in particular in its interactions with its host.

The role of peptidoglycan in chlamydial cell division: towards resolving the chlamydial anomaly

Chlamydiales are obligate intracellular bacteria including some important pathogens causing trachoma, genital tract infections and pneumonia, among others. They share an atypical division mechanism, which is independent of an FtsZ homologue. However, they divide by binary fission, in a process inhibited by penicillin derivatives, causing the formation of an aberrant form of the bacteria, which is able to survive in the presence of the antibiotic. The paradox of penicillin sensitivity of chlamydial cells in the absence of detectable peptidoglycan (PG) was dubbed the chlamydial anomaly, since no PG modified by enzymes (Pbps) that are the usual target of penicillin could be detected in Chlamydiales. We review here the recent advances in this field with the first direct and indirect evidences of PG-like material in both Chlamydiaceae and Chlamydia-related bacteria. Moreover, PG biosynthesis is required for proper localization of the newly described septal proteins RodZ and NlpD. Taken together, these new results set the stage for a better understanding of the role of PG and septal proteins in the division mechanism of Chlamydiales and illuminate the long-standing chlamydial anomaly. Moreover, understanding the chlamydial division mechanism is critical for the development of new antibiotics for the treatment of chlamydial chronic infections.

Evidence for horizontal gene transfer between Chlamydophila pneumoniae and Chlamydia phage

Chlamydia-infecting bacteriophages, members of the Microviridae family, specifically the Gokushovirinae subfamily, are small (4.5–5 kb) single-stranded circles with 8–10 open-reading frames similar to E. coli phage ϕX174. Using sequence information found in GenBank, we examined related genes in Chlamydophila pneumoniae and Chlamydia-infecting bacteriophages. The 5 completely sequenced C. pneumoniae strains contain a gene orthologous to a phage gene annotated as the putative replication initiation protein (PRIP, also called VP4), which is not found in any other members of the Chlamydiaceae family sequenced to date. The C. pneumoniae strain infecting koalas, LPCoLN, in addition contains another region orthologous to phage sequences derived from the minor capsid protein gene, VP3. Phylogenetically, the phage PRIP sequences are more diverse than the bacterial PRIP sequences; nevertheless, the bacterial sequences and the phage sequences each cluster together in their own clade. Finally, we found evidence for another Microviridae phage-related gene, the major capsid protein gene, VP1 in a number of other bacterial species and 2 eukaryotes, the woodland strawberry and a nematode. Thus, we find considerable evidence for DNA sequences related to genes found in bacteriophages of the Microviridae family not only in a variety of prokaryotic but also eukaryotic species.

Non-coding RNA regulation in pathogenic bacteria located inside eukaryotic cells

Intracellular bacterial pathogens have evolved distinct lifestyles inside eukaryotic cells. Some pathogens coexist with the infected cell in an obligate intracellular state, whereas others transit between the extracellular and intracellular environment. Adaptation to these intracellular lifestyles is regulated in both space and time. Non-coding small RNAs (sRNAs) are post-transcriptional regulatory molecules that fine-tune important processes in bacterial physiology including cell envelope architecture, intermediate metabolism, bacterial communication, biofilm formation, and virulence. Recent studies have shown production of defined sRNA species by intracellular bacteria located inside eukaryotic cells. The molecules targeted by these sRNAs and their expression dynamics along the intracellular infection cycle remain, however, poorly characterized. Technical difficulties linked to the isolation of “intact” intracellular bacteria from infected host cells might explain why sRNA regulation in these specialized pathogens is still a largely unexplored field. Transition from the extracellular to the intracellular lifestyle provides an ideal scenario in which regulatory sRNAs are intended to participate; so much work must be done in this direction. This review focuses on sRNAs expressed by intracellular bacterial pathogens during the infection of eukaryotic cells, strategies used with these pathogens to identify sRNAs required for virulence, and the experimental technical challenges associated to this type of studies. We also discuss varied techniques for their potential application to study RNA regulation in intracellular bacterial infections.

Relaxed selection drives a noisy noncoding transcriptome in members of the Mycobacterium tuberculosis complex

Related species are often used to understand the molecular underpinning of virulence through examination of a shared set of biological features attributable to a core genome of orthologous genes. An important but insufficiently studied issue, however, is the extent to which the regulatory architectures are similarly conserved. A small number of studies have compared the primary transcriptomes of different bacterial species, but few have compared closely related species with clearly divergent evolutionary histories. We addressed the impact of differing modes of evolution within the genus Mycobacterium through comparison of the primary transcriptome of M. marinum with that of a closely related lineage, M. bovis. Both are thought to have evolved from an ancestral generalist species, with M. bovis and other members of the M. tuberculosis complex having subsequently undergone downsizing of their genomes during the transition to obligate pathogenicity. M. marinum, in contrast, has retained a large genome, appropriate for an environmental organism, and is a broad-host-range pathogen. We also examined changes over a shorter evolutionary time period through comparison of the primary transcriptome of M. bovis with that of another member of the M. tuberculosis complex (M. tuberculosis) which possesses an almost identical genome but maintains a distinct host preference.

Our comparison of the transcriptional start site (TSS) maps of M. marinum and M. bovis uncovers a pillar of conserved promoters, noncoding RNA (NCRNA), and a genome-wide signal in the −35 promoter regions of both species. We identify evolutionarily conserved transcriptional attenuation and highlight its potential contribution to multidrug resistance mediated through the transcriptional regulator whiB7. We show that a species population history is reflected in its transcriptome and posit relaxed selection as the main driver of an abundance of canonical −10 promoter sites in M. bovis relative to M. marinum. It appears that transcriptome composition in mycobacteria is driven primarily by the availability of such sites and that their frequencies diverge significantly across the mycobacterial clade. Finally, through comparison of M. bovis and M. tuberculosis, we illustrate that single nucleotide polymorphism (SNP)-driven promoter differences likely underpin many of the transcriptional differences between M. tuberculosis complex lineages.

Single-cell RNA-seq: advances and future challenges

Phenotypically identical cells can dramatically vary with respect to behavior during their lifespan and this variation is reflected in their molecular composition such as the transcriptomic landscape. Single-cell transcriptomics using next-generation transcript sequencing (RNA-seq) is now emerging as a powerful tool to profile cell-to-cell variability on a genomic scale. Its application has already greatly impacted our conceptual understanding of diverse biological processes with broad implications for both basic and clinical research. Different single-cell RNA-seq protocols have been introduced and are reviewed here—each one with its own strengths and current limitations. We further provide an overview of the biological questions single-cell RNA-seq has been used to address, the major findings obtained from such studies, and current challenges and expected future developments in this booming field.

Differential RNA-seq: the approach behind and the biological insight gained

RNA-sequencing has revolutionized the quantitative and qualitative analysis of transcriptomes in both prokaryotes and eukaryotes. It provides a generic approach for gene expression profiling, annotation of transcript boundaries and operons, as well as identifying novel transcripts including small noncoding RNA molecules and antisense RNAs. We recently developed a differential RNA-seq (dRNA-seq) method which in addition to the above, yields information as to whether a given RNA is a primary or processed transcript. Originally applied to describe the primary transcriptome of the gastric pathogen Helicobacter pylori, dRNA-seq has since provided global maps of transcriptional start sites in diverse species, informed new biology in the CRISPR-Cas9 system, advanced to a tool for comparative transcriptomics, and inspired simultaneous RNA-seq of pathogen and host.

Riboregulators and the role of Hfq in photosynthetic bacteria

Anoxygenic and oxygenic bacteria directly convert solar energy into biomass using photosynthesis. The formation and composition of photosynthetic complexes has to be tightly controlled in response to environmental conditions, as exposure to sunlight can be harmful due to the generation of reactive oxygen species and the damaging effects of UV irradiation. Therefore, photosynthetic bacteria are exposed to a particular set of regulatory challenges in addition to those that also affect other bacteria, requiring sophisticated regulatory systems. Indeed, hundreds of potential regulatory RNAs have been identified in photosynthetic model bacteria as well as antisense RNAs (asRNAs) of up to several kb in length that protect certain mRNAs from degradation. The trans-acting small non-coding RNAs (sRNAs), PcrZ and PsrR1, control pigment and photosystem biogenesis in Rhodobacter sphaeroides and cyanobacteria, respectively. The asRNAs IsrR and As1_flv4 act as negative regulators and the asRNAs PsbA2R and PsbA3R as positive effectors of photosynthesis gene expression in Synechocystis 6803.

Structural properties of prokaryotic promoter regions correlate with functional features

The structural properties of the DNA molecule are known to play a critical role in transcription. In this paper, the structural profiles of promoter regions were studied within the context of their diversity and their function for eleven prokaryotic species; Escherichia coli, Klebsiella pneumoniae, Salmonella Typhimurium, Pseudomonas auroginosa, Geobacter sulfurreducens Helicobacter pylori, Chlamydophila pneumoniae, Synechocystis sp., Synechoccocus elongates, Bacillus anthracis, and the archaea Sulfolobus solfataricus. The main anchor point for these promoter regions were transcription start sites identified through high-throughput experiments or collected within large curated databases. Prokaryotic promoter regions were found to be less stable and less flexible than the genomic mean across all studied species. However, direct comparison between species revealed differences in their structural profiles that can not solely be explained by the difference in genomic GC content. In addition, comparison with functional data revealed that there are patterns in the promoter structural profiles that can be linked to specific functional loci, such as sigma factor regulation or transcription factor binding. Interestingly, a novel structural element clearly visible near the transcription start site was found in genes associated with essential cellular functions and growth in several species. Our analyses reveals the great diversity in promoter structural profiles both between and within prokaryotic species. We observed relationships between structural diversity and functional features that are interesting prospects for further research to yet uncharacterized functional loci defined by DNA structural properties.

Ultra deep sequencing of Listeria monocytogenes sRNA transcriptome revealed new antisense RNAs

Listeria monocytogenes, a gram-positive pathogen, and causative agent of listeriosis, has become a widely used model organism for intracellular infections. Recent studies have identified small non-coding RNAs (sRNAs) as important factors for regulating gene expression and pathogenicity of L. monocytogenes. Increased speed and reduced costs of high throughput sequencing (HTS) techniques have made RNA sequencing (RNA-Seq) the state-of-the-art method to study bacterial transcriptomes. We created a large transcriptome dataset of L. monocytogenes containing a total of 21 million reads, using the SOLiD sequencing technology. The dataset contained cDNA sequences generated from L. monocytogenes RNA collected under intracellular and extracellular condition and additionally was size fractioned into three different size ranges from <40 nt, 40–150 nt and >150 nt. We report here, the identification of nine new sRNAs candidates of L. monocytogenes and a reevaluation of known sRNAs of L. monocytogenes EGD-e. Automatic comparison to known sRNAs revealed a high recovery rate of 55%, which was increased to 90% by manual revision of the data. Moreover, thorough classification of known sRNAs shed further light on their possible biological functions. Interestingly among the newly identified sRNA candidates are antisense RNAs (asRNAs) associated to the housekeeping genes purA, fumC and pgi and potentially their regulation, emphasizing the significance of sRNAs for metabolic adaptation in L. monocytogenes.

TSSer: an automated method to identify transcription start sites in prokaryotic genomes from differential RNA sequencing data

MOTIVATION

Accurate identification of transcription start sites (TSSs) is an essential step in the analysis of transcription regulatory networks. In higher eukaryotes, the capped analysis of gene expression technology enabled comprehensive annotation of TSSs in genomes such as those of mice and humans. In bacteria, an equivalent approach, termed differential RNA sequencing (dRNA-seq), has recently been proposed, but the application of this approach to a large number of genomes is hindered by the paucity of computational analysis methods. With few exceptions, when the method has been used, annotation of TSSs has been largely done manually.

RESULTS

In this work, we present a computational method called 'TSSer' that enables the automatic inference of TSSs from dRNA-seq data. The method rests on a probabilistic framework for identifying both genomic positions that are preferentially enriched in the dRNA-seq data as well as preferentially captured relative to neighboring genomic regions. Evaluating our approach for TSS calling on several publicly available datasets, we find that TSSer achieves high consistency with the curated lists of annotated TSSs, but identifies many additional TSSs. Therefore, TSSer can accelerate genome-wide identification of TSSs in bacterial genomes and can aid in further characterization of bacterial transcription regulatory networks.

AVAILABILITY

TSSer is freely available under GPL license at http://www.clipz.unibas.ch/TSSer/index.php

A new metabolic cell-wall labelling method reveals peptidoglycan in Chlamydia trachomatis

Peptidoglycan (PG), an essential structure in the cell walls of the vast majority of bacteria, is critical for division and maintaining cell shape and hydrostatic pressure. Bacteria comprising the Chlamydiales were thought to be one of the few exceptions. Chlamydia harbour genes for PG biosynthesis and exhibit susceptibility to 'anti-PG' antibiotics, yet attempts to detect PG in any chlamydial species have proven unsuccessful (the 'chlamydial anomaly'). We used a novel approach to metabolically label chlamydial PG using d-amino acid dipeptide probes and click chemistry. Replicating Chlamydia trachomatis were labelled with these probes throughout their biphasic developmental life cycle, and the results of differential probe incorporation experiments conducted in the presence of ampicillin are consistent with the presence of chlamydial PG-modifying enzymes. These findings culminate 50 years of speculation and debate concerning the chlamydial anomaly and are the strongest evidence so far that chlamydial species possess functional PG.

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.