Refine
Has Fulltext
- yes (25)
Is part of the Bibliography
- yes (25)
Year of publication
Document Type
- Journal article (25)
Language
- English (25)
Keywords
- RNA-seq (4)
- escherichia coli (3)
- transcriptomics (3)
- Antisense RNA (2)
- Rhodobacter sphaeroides (2)
- bacteria (2)
- dRNA-Seq (2)
- gene expression (2)
- genome annotation (2)
- identification (2)
Institute
Sonstige beteiligte Institutionen
Background
Differential RNA-Seq (dRNA-Seq) is a recently developed method of performing primary transcriptome analyses that allows for the genome-wide mapping of transcriptional start sites (TSSs) and the identification of novel transcripts. Although the transcriptomes of diverse bacterial species have been characterized by dRNA-Seq, the transcriptome analysis of archaeal species is still rather limited. Therefore, we used dRNA-Seq to characterize the primary transcriptome of the model archaeon Haloferax volcanii.
Results
Three independent cultures of Hfx. volcanii grown under optimal conditions to the mid-exponential growth phase were used to determine the primary transcriptome and map the 5′-ends of the transcripts. In total, 4749 potential TSSs were detected. A position weight matrix (PWM) was derived for the promoter predictions, and the results showed that 64 % of the TSSs were preceded by stringent or relaxed basal promoters. Of the identified TSSs, 1851 belonged to protein-coding genes. Thus, fewer than half (46 %) of the 4040 protein-coding genes were expressed under optimal growth conditions. Seventy-two percent of all protein-coding transcripts were leaderless, which emphasized that this pathway is the major pathway for translation initiation in haloarchaea. A total of 2898 of the TSSs belonged to potential non-coding RNAs, which accounted for an unexpectedly high fraction (61 %) of all transcripts. Most of the non-coding TSSs had not been previously described (2792) and represented novel sequences (59 % of all TSSs). A large fraction of the potential novel non-coding transcripts were cis-antisense RNAs (1244 aTSSs). A strong negative correlation between the levels of antisense transcripts and cognate sense mRNAs was found, which suggested that the negative regulation of gene expression via antisense RNAs may play an important role in haloarchaea. The other types of novel non-coding transcripts corresponded to internal transcripts overlapping with mRNAs (1153 iTSSs) and intergenic small RNA (sRNA) candidates (395 TSSs).
Conclusion
This study provides a comprehensive map of the primary transcriptome of Hfx. volcanii grown under optimal conditions. Fewer than half of all protein-coding genes have been transcribed under these conditions. Unexpectedly, more than half of the detected TSSs belonged to several classes of non-coding RNAs. Thus, RNA-based regulation appears to play a more important role in haloarchaea than previously anticipated.
Background
Differential RNA-sequencing (dRNA-seq) is indispensable for determination of primary transcriptomes. However, using dRNA-seq data to map transcriptional start sites (TSSs) and promoters genome-wide is a bioinformatics challenge. We performed dRNA-seq of Bradyrhizobium japonicum USDA 110, the nitrogen-fixing symbiont of soybean, and developed algorithms to map TSSs and promoters.
Results
A specialized machine learning procedure for TSS recognition allowed us to map 15,923 TSSs: 14,360 in free-living bacteria, 4329 in symbiosis with soybean and 2766 in both conditions. Further, we provide proteomic evidence for 4090 proteins, among them 107 proteins corresponding to new genes and 178 proteins with N-termini different from the existing annotation (72 and 109 of them with TSS support, respectively). Guided by proteomics evidence, previously identified TSSs and TSSs experimentally validated here, we assign a score threshold to flag 14 % of the mapped TSSs as a class of lower confidence. However, this class of lower confidence contains valid TSSs of low-abundant transcripts. Moreover, we developed a de novo algorithm to identify promoter motifs upstream of mapped TSSs, which is publicly available, and found motifs mainly used in symbiosis (similar to RpoN-dependent promoters) or under both conditions (similar to RpoD-dependent promoters). Mapped TSSs and putative promoters, proteomic evidence and updated gene annotation were combined into an annotation file.
Conclusions
The genome-wide TSS and promoter maps along with the extended genome annotation of B. japonicum represent a valuable resource for future systems biology studies and for detailed analyses of individual non-coding transcripts and ORFs. Our data will also provide new insights into bacterial gene regulation during the agriculturally important symbiosis between rhizobia and legumes.
Genome-wide transcription start site profiling in biofilm-grown Burkholderia cenocepacia J2315
(2015)
Background: Burkholderia cenocepacia is a soil-dwelling Gram-negative Betaproteobacterium with an important role as opportunistic pathogen in humans. Infections with B. cenocepacia are very difficult to treat due to their high intrinsic resistance to most antibiotics. Biofilm formation further adds to their antibiotic resistance. B. cenocepacia harbours a large, multi-replicon genome with a high GC-content, the reference genome of strain J2315 includes 7374 annotated genes. This study aims to annotate transcription start sites and identify novel transcripts on a whole genome scale. Methods: RNA extracted from B. cenocepacia J2315 biofilms was analysed by differential RNA-sequencing and the resulting dataset compared to data derived from conventional, global RNA-sequencing. Transcription start sites were annotated and further analysed according to their position relative to annotated genes. Results: Four thousand ten transcription start sites were mapped over the whole B. cenocepacia genome and the primary transcription start site of 2089 genes expressed in B. cenocepacia biofilms were defined. For 64 genes a start codon alternative to the annotated one was proposed. Substantial antisense transcription for 105 genes and two novel protein coding sequences were identified. The distribution of internal transcription start sites can be used to identify genomic islands in B. cenocepacia. A potassium pump strongly induced only under biofilm conditions was found and 15 non-coding small RNAs highly expressed in biofilms were discovered. Conclusions: Mapping transcription start sites across the B. cenocepacia genome added relevant information to the J2315 annotation. Genes and novel regulatory RNAs putatively involved in B. cenocepacia biofilm formation were identified. These findings will help in understanding regulation of B. cenocepacia biofilm formation.
Campylobacter jejuni is currently the leading cause of bacterial gastroenteritis in humans. Comparison of multiple Campylobacter strains revealed a high genetic and phenotypic diversity. However, little is known about differences in transcriptome organization, gene expression, and small RNA (sRNA) repertoires. Here we present the first comparative primary transcriptome analysis based on the differential RNA–seq (dRNA–seq) of four C. jejuni isolates. Our approach includes a novel, generic method for the automated annotation of transcriptional start sites (TSS), which allowed us to provide genome-wide promoter maps in the analyzed strains. These global TSS maps are refined through the integration of a SuperGenome approach that allows for a comparative TSS annotation by mapping RNA–seq data of multiple strains into a common coordinate system derived from a whole-genome alignment. Considering the steadily increasing amount of RNA–seq studies, our automated TSS annotation will not only facilitate transcriptome annotation for a wider range of pro- and eukaryotes but can also be adapted for the analysis among different growth or stress conditions. Our comparative dRNA–seq analysis revealed conservation of most TSS, but also single-nucleotide-polymorphisms (SNP) in promoter regions, which lead to strain-specific transcriptional output. Furthermore, we identified strain-specific sRNA repertoires that could contribute to differential gene regulation among strains. In addition, we identified a novel minimal CRISPR-system in Campylobacter of the type-II CRISPR subtype, which relies on the host factor RNase III and a trans-encoded sRNA for maturation of crRNAs. This minimal system of Campylobacter, which seems active in only some strains, employs a unique maturation pathway, since the crRNAs are transcribed from individual promoters in the upstream repeats and thereby minimize the requirements for the maturation machinery. Overall, our study provides new insights into strain-specific transcriptome organization and sRNAs, and reveals genes that could modulate phenotypic variation among strains despite high conservation at the DNA level.
Converging evidence suggests a role of serotonin (5-hydroxytryptamine, 5-HT) and tryptophan hydroxylase 2 (TPH2), the rate-limiting enzyme of 5-HT synthesis in the brain, in modulating long-term, neurobiological effects of early-life adversity. Here, we aimed at further elucidating the molecular mechanisms underlying this interaction, and its consequences for socio-emotional behaviors, with a focus on anxiety and social interaction. In this study, adult, male Tph2 null mutant (Tph2\(^{-/-}\)) and heterozygous (Tph2\(^{+/-}\)) mice, and their wildtype littermates (Tph2\(^{+/+}\)) were exposed to neonatal, maternal separation (MS) and screened for behavioral changes, followed by genome-wide RNA expression and DNA methylation profiling. In Tph2\(^{-/-}\) mice, brain 5-HT deficiency profoundly affected socio-emotional behaviors, i.e., decreased avoidance of the aversive open arms in the elevated plus-maze (EPM) as well as decreased prosocial and increased rule breaking behavior in the resident-intruder test when compared to their wildtype littermates. Tph2\(^{+/-}\) mice showed an ambiguous profile with context-dependent, behavioral responses. In the EPM they showed similar avoidance of the open arm but decreased prosocial and increased rule breaking behavior in the resident-intruder test when compared to their wildtype littermates. Notably, MS effects on behavior were subtle and depended on the Tph2 genotype, in particular increasing the observed avoidance of EPM open arms in wildtype and Tph2\(^{+/-}\) mice when compared to their Tph2\(^{-/-}\) littermates. On the genomic level, the interaction of Tph2 genotype with MS differentially affected the expression of numerous genes, of which a subset showed an overlap with DNA methylation profiles at corresponding loci. Remarkably, changes in methylation nearby and expression of the gene encoding cholecystokinin, which were inversely correlated to each other, were associated with variations in anxiety-related phenotypes. In conclusion, next to various behavioral alterations, we identified gene expression and DNA methylation profiles to be associated with TPH2 inactivation and its interaction with MS, suggesting a gene-by-environment interaction-dependent, modulatory function of brain 5-HT availability.
Meningococcal meningitis is a severe central nervous system infection that occurs when Neisseria meningitidis (Nm) penetrates brain endothelial cells (BECs) of the meningeal blood-cerebrospinal fluid barrier. As a human-specific pathogen, in vivo models are greatly limited and pose a significant challenge. In vitro cell models have been developed, however, most lack critical BEC phenotypes limiting their usefulness. Human BECs generated from induced pluripotent stem cells (iPSCs) retain BEC properties and offer the prospect of modeling the human-specific Nm interaction with BECs. Here, we exploit iPSC-BECs as a novel cellular model to study Nm host-pathogen interactions, and provide an overview of host responses to Nm infection. Using iPSC-BECs, we first confirmed that multiple Nm strains and mutants follow similar phenotypes to previously described models. The recruitment of the recently published pilus adhesin receptor CD147 underneath meningococcal microcolonies could be verified in iPSC-BECs. Nm was also observed to significantly increase the expression of pro-inflammatory and neutrophil-specific chemokines IL6, CXCL1, CXCL2, CXCL8, and CCL20, and the secretion of IFN-γ and RANTES. For the first time, we directly observe that Nm disrupts the three tight junction proteins ZO-1, Occludin, and Claudin-5, which become frayed and/or discontinuous in BECs upon Nm challenge. In accordance with tight junction loss, a sharp loss in trans-endothelial electrical resistance, and an increase in sodium fluorescein permeability and in bacterial transmigration, was observed. Finally, we established RNA-Seq of sorted, infected iPSC-BECs, providing expression data of Nm-responsive host genes. Altogether, this model provides novel insights into Nm pathogenesis, including an impact of Nm on barrier properties and tight junction complexes, and suggests that the paracellular route may contribute to Nm traversal of BECs.
Microalga are of high relevance for the global carbon cycling and it is well-known that they are associated with a microbiota. However, it remains unclear, if the associated microbiota, often found in phycosphere biofilms, is specific for the microalga strains and which role individual bacterial taxa play. Here we provide experimental evidence that \(Chlorella\) \(saccharophila\), \(Scenedesmus\) \(quadricauda\), and \(Micrasterias\) \(crux-melitensis\), maintained in strain collections, are associated with unique and specific microbial populations. Deep metagenome sequencing, binning approaches, secretome analyses in combination with RNA-Seq data implied fundamental differences in the gene expression profiles of the microbiota associated with the different microalga. Our metatranscriptome analyses indicates that the transcriptionally most active bacteria with respect to key genes commonly involved in plant–microbe interactions in the Chlorella (Trebouxiophyceae) and Scenedesmus (Chlorophyceae) strains belong to the phylum of the α-Proteobacteria. In contrast, in the Micrasterias (Zygnematophyceae) phycosphere biofilm bacteria affiliated with the phylum of the Bacteroidetes showed the highest gene expression rates. We furthermore show that effector molecules known from plant-microbe interactions as inducers for the innate immunity are already of relevance at this evolutionary early plant-microbiome level.
Integrative "Omics"-Approach Discovers Dynamic and Regulatory Features of Bacterial Stress Responses
(2013)
Bacteria constantly face stress conditions and therefore mount specific responses to ensure adaptation and survival. Stress responses were believed to be predominantly regulated at the transcriptional level. In the phototrophic bacterium Rhodobacter sphaeroides the response to singlet oxygen is initiated by alternative sigma factors. Further adaptive mechanisms include post-transcriptional and post-translational events, which have to be considered to gain a deeper understanding of how sophisticated regulation networks operate. To address this issue, we integrated three layers of regulation: (1) total mRNA levels at different time-points revealed dynamics of the transcriptome, (2) mRNAs in polysome fractions reported on translational regulation (translatome), and (3) SILAC-based mass spectrometry was used to quantify protein abundances (proteome). The singlet oxygen stress response exhibited highly dynamic features regarding short-term effects and late adaptation, which could in part be assigned to the sigma factors RpoE and RpoH2 generating distinct expression kinetics of corresponding regulons. The occurrence of polar expression patterns of genes within stress-inducible operons pointed to an alternative of dynamic fine-tuning upon stress. In addition to transcriptional activation, we observed significant induction of genes at the post-transcriptional level (translatome), which identified new putative regulators and assigned genes of quorum sensing to the singlet oxygen stress response. Intriguingly, the SILAC approach explored the stress-dependent decline of photosynthetic proteins, but also identified 19 new open reading frames, which were partly validated by RNA-seq. We propose that comparative approaches as presented here will help to create multi-layered expression maps on the system level ("expressome"). Finally, intense mass spectrometry combined with RNA-seq might be the future tool of choice to re-annotate genomes in various organisms and will help to understand how they adapt to alternating conditions.
Phenotypic heterogeneity at the cellular level in response to various stresses, e.g., antibiotic treatment has been reported for a number of bacteria. In a clonal population, cell-to-cell variation may result in phenotypic heterogeneity that is a mechanism to survive changing environments including antibiotic therapy. Stenotrophomonas rnaltophilia has been frequently isolated from cystic fibrosis patients, can cause numerous infections in other organs and tissues, and is difficult to treat due to antibiotic resistances. S. maltophilia K279a produces the Li and L2 beta-lactamases in response to beta-lactam treatment. Here we report that the patient isolate S. rnaltophilia K279a diverges into cellular subpopulations with distinct but reversible morphotypes of small and big colonies when challenged with ampicillin. This observation is consistent with the formation of elongated chains of bacteria during exponential growth phase and the occurrence of mainly rod-shaped cells in liquid media. RNA-seq analysis of small versus big colonies revealed differential regulation of at least seven genes among the colony morphotypes. Among those, bleu and bla(L2) were transcriptionally the most strongly upregulated genes. Promoter fusions of b/a(L1) and b/a(L2) genes indicated that expression of both genes is also subject to high levels of phenotypic heterogeneous expression on a single cell level. Additionally, the comE homolog was found to be differentially expressed in homogenously versus heterogeneously bla(L2) expressing cells as identified by RNA(seq) analysis. Overexpression of cornE in S. maltophilia K279a reduced the level of cells that were in a bla(L2)-ON mode to 1% or lower. Taken together, our data provide strong evidence that S. maltophilia K279a populations develop phenotypic heterogeneity in an ampicillin challenged model. This cellular variability is triggered by regulation networks including b/a(L1), b/a(L2), and comE.
Background
Prokaryotes have relatively small genomes, densely-packed with protein-encoding sequences. RNA sequencing has, however, revealed surprisingly complex transcriptomes and here we report the transcripts present in the model hyperthermophilic Archaeon, Thermococcus kodakarensis, under different physiological conditions.
Results
Sequencing cDNA libraries, generated from RNA isolated from cells under different growth and metabolic conditions has identified >2,700 sites of transcription initiation, established a genome-wide map of transcripts, and consensus sequences for transcription initiation and post-transcription regulatory elements. The primary transcription start sites (TSS) upstream of 1,254 annotated genes, plus 644 primary TSS and their promoters within genes, are identified. Most mRNAs have a 5'-untranslated region (5'-UTR) 10 to 50 nt long (median = 16 nt), but ~20% have 5'-UTRs from 50 to 300 nt long and ~14% are leaderless. Approximately 50% of mRNAs contain a consensus ribosome binding sequence. The results identify TSS for 1,018 antisense transcripts, most with sequences complementary to either the 5'- or 3'-region of a sense mRNA, and confirm the presence of transcripts from all three CRISPR loci, the RNase P and 7S RNAs, all tRNAs and rRNAs and 69 predicted snoRNAs. Two putative riboswitch RNAs were present in growing but not in stationary phase cells. The procedure used is designed to identify TSS but, assuming that the number of cDNA reads correlates with transcript abundance, the results also provide a semi-quantitative documentation of the differences in T. kodakarensis genome expression under different growth conditions and confirm previous observations of substrate-dependent specific gene expression. Many previously unanticipated small RNAs have been identified, some with relative low GC contents (≤50%) and sequences that do not fold readily into base-paired secondary structures, contrary to the classical expectations for non-coding RNAs in a hyperthermophile.
Conclusion
The results identify >2,700 TSS, including almost all of the primary sites of transcription initiation upstream of annotated genes, plus many secondary sites, sites within genes and sites resulting in antisense transcripts. The T. kodakarensis genome is small (~2.1 Mbp) and tightly packed with protein-encoding genes, but the transcriptomes established also contain many non-coding RNAs and predict extensive RNA-based regulation in this model Archaeon.