RNA-Seq-based analysis of the physiologic cold shock-induced changes in Moraxella catarrhalis gene expression

BACKGROUND Moraxella catarrhalis, a major nasopharyngeal pathogen of the human respiratory tract, is exposed to rapid downshifts of environmental temperature when humans breathe cold air. The prevalence of pharyngeal colonization and respiratory tract infections caused by M. catarrhalis is greatest in winter. We investigated how M. catarrhalis uses the physiologic exposure to cold air to regulate pivotal survival systems that may contribute to M. catarrhalis virulence. RESULTS In this study we used the RNA-seq techniques to quantitatively catalogue the transcriptome of M. catarrhalis exposed to a 26 °C cold shock or to continuous growth at 37 °C. Validation of RNA-seq data using quantitative RT-PCR analysis demonstrated the RNA-seq results to be highly reliable. We observed that a 26 °C cold shock induces the expression of genes that in other bacteria have been related to virulence a strong induction was observed for genes involved in high affinity phosphate transport and iron acquisition, indicating that M. catarrhalis makes a better use of both phosphate and iron resources after exposure to cold shock. We detected the induction of genes involved in nitrogen metabolism, as well as several outer membrane proteins, including ompA, m35-like porin and multidrug efflux pump (acrAB) indicating that M. catarrhalis remodels its membrane components in response to downshift of temperature. Furthermore, we demonstrate that a 26 °C cold shock enhances the induction of genes encoding the type IV pili that are essential for natural transformation, and increases the genetic competence of M. catarrhalis, which may facilitate the rapid spread and acquisition of novel virulence-associated genes. CONCLUSION Cold shock at a physiologically relevant temperature of 26 °C induces in M. catarrhalis a complex of adaptive mechanisms that could convey novel pathogenic functions and may contribute to enhanced colonization and virulence.

Evidence for a retroviral insertion in TRPM1 as the cause of congenital stationary night blindness and leopard complex spotting in the horse.

Leopard complex spotting is a group of white spotting patterns in horses caused by an incompletely dominant gene (LP) where homozygotes (LP/LP) are also affected with congenital stationary night blindness. Previous studies implicated Transient Receptor Potential Cation Channel, Subfamily M, Member 1 (TRPM1) as the best candidate gene for both CSNB and LP. RNA-Seq data pinpointed a 1378 bp insertion in intron 1 of TRPM1 as the potential cause. This insertion, a long terminal repeat (LTR) of an endogenous retrovirus, was completely associated with LP, testing 511 horses (χ(2)=1022.00, p<0.0005), and CSNB, testing 43 horses (χ(2)=43, p<0.0005). The LTR was shown to disrupt TRPM1 transcription by premature poly-adenylation. Furthermore, while deleterious transposable element insertions should be quickly selected against the identification of this insertion in three ancient DNA samples suggests it has been maintained in the horse gene pool for at least 17,000 years. This study represents the first description of an LTR insertion being associated with both a pigmentation phenotype and an eye disorder.

Identification and characterisation of ncRNAs in dormant bacteria

Recent years have led to increasing interest and appreciation of the possible importance of single cell heterogeneity in various biological processes. One of the examples of phenotypic heterogeneity in bacterial populations is antibiotic tolerant persister cells. Such an antibiotic tolerance phenotype is of considerable clinical relevance since dormant bacteria can re-establish infections rapidly after the antibiotic treatment has been terminated. Up to now mechanisms for establishing the persistence phenomenon in bacteria have remained largely enigmatic. Persisters are cells considered to be in a dormant state with down regulated gene expression. Only recently small regulatory RNAs (sRNAs) have been appreciated as important regulators of gene expression in response to environmental stimuli and several theoretical studies have suggested a possible involvement of sRNAs in the mechanisms of regulated heterogeneity in bacteria. We have experimentally addressed this potential link between sRNAs and persistence/dormancy in E. coli as an example of heterogeneity. Beside classical sRNAs we are focusing also on sRNAs directly associating with and possibly regulating the ribosome, the central enzyme of gene expression. The persister and dormant cell specific sRNA profile is studied by the comparative analysis of sRNA profile changes of the whole bacterial population after antibiotic killing. From RNA-Seq data ~ 25 000 potentially stable RNA fragments were identified and initial analysis predicted ~300 of them to be dormant/persister cell specific. After further evaluation the most prominent dormant/persister cell specific sRNAs are functionally characterized and their potential role in the persistence/dormancy will be evaluated by applying genetic, molecular and biochemical tools. The potential results of this project will provide a better understanding on the molecular mechanism of bacterial persistence/dormancy and on the role of ribosome-bound sRNA molecules in fine-tuning gene expression.

Identification and characterisation of ncRNAs in dormant bacteria

Recent years have led to increasing interest and appreciation of the possible importance of single cell heterogeneity in various biological processes. One of the examples of phenotypic heterogeneity in bacterial populations is antibiotic tolerant persister cells. Such an antibiotic tolerance phenotype is of considerable clinical relevance since dormant bacteria can re-establish infections rapidly after the antibiotic treatment has been terminated. Up to now mechanisms for establishing the persistence phenomenon in bacteria have remained largely enigmatic. Persisters are cells considered to be in a dormant state with down regulated gene expression. Only recently small regulatory RNAs (sRNAs) have been appreciated as important regulators of gene expression in response to environmental stimuli and several theoretical studies have suggested a possible involvement of sRNAs in the mechanisms of regulated heterogeneity in bacteria. We have experimentally addressed this potential link between sRNAs and persistence/dormancy in E. coli as an example of heterogeneity. Beside classical sRNAs we are focusing also on sRNAs directly associating with and possibly regulating the ribosome, the central enzyme of gene expression. The persister and dormant cell specific sRNA profile is studied by the comparative analysis of sRNA profile changes of the whole bacterial population after antibiotic killing. From RNA-Seq data ~ 25 000 potentially stable RNA fragments were identified and initial analysis predicted ~300 of them to be dormant/persister cell specific. After further evaluation the most prominent dormant/persister cell specific sRNAs are functionally characterized and their potential role in the persistence/dormancy will be evaluated by applying genetic, molecular and biochemical tools. The potential results of this project will provide a better understanding on the molecular mechanism of bacterial persistence/dormancy and on the role of ribosome-bound sRNA molecules in fine-tuning gene expression.

Identification and characterisation of ncRNAs in dormant bacteria

Recent years have led to increasing interest and appreciation of the possible importance of single cell heterogeneity in various biological processes. One of the examples of phenotypic heterogeneity in bacterial populations is antibiotic tolerant persister cells. Such an antibiotic tolerance phenotype is of considerable clinical relevance since dormant bacteria can re-establish infections rapidly after the antibiotic treatment has been terminated. Up to now mechanisms for establishing the persistence phenomenon in bacteria have remained largely enigmatic. Persisters are cells considered to be in a dormant state with down regulated gene expression. Only recently small regulatory RNAs (sRNAs) have been appreciated as important regulators of gene expression in response to environmental stimuli and several theoretical studies have suggested a possible involvement of sRNAs in the mechanisms of regulated heterogeneity in bacteria. We have experimentally addressed this potential link between sRNAs and persistence/dormancy in E. coli as an example of heterogeneity. Beside classical sRNAs we are focusing also on sRNAs directly associating with and possibly regulating the ribosome, the central enzyme of gene expression. The persister and dormant cell specific sRNA profile is studied by the comparative analysis of sRNA profile changes of the whole bacterial population after antibiotic killing. From RNA-Seq data ~ 25 000 potentially stable RNA fragments were identified and initial analysis predicted ~300 of them to be dormant/persister cell specific. After further evaluation the most prominent dormant/persister cell specific sRNAs are functionally characterized and their potential role in the persistence/dormancy will be evaluated by applying genetic, molecular and biochemical tools. The potential results of this project will provide a better understanding on the molecular mechanism of bacterial persistence/dormancy and on the role of ribosome-bound sRNA molecules in fine-tuning gene expression.

The transcriptome of equine peripheral blood mononuclear cells.

Complete transcriptomic data at high resolution are available only for a few model organisms with medical importance. The gene structures of non-model organisms are mostly computationally predicted based on comparative genomics with other species. As a result, more than half of the horse gene models are known only by projection. Experimental data supporting these gene models are scarce. Moreover, most of the annotated equine genes are single-transcript genes. Utilizing RNA sequencing (RNA-seq) the experimental validation of predicted transcriptomes has become accessible at reasonable costs. To improve the horse genome annotation we performed RNA-seq on 561 samples of peripheral blood mononuclear cells (PBMCs) derived from 85 Warmblood horses. The mapped sequencing reads were used to build a new transcriptome assembly. The new assembly revealed many alternative isoforms associated to known genes or to those predicted by the Ensembl and/or Gnomon pipelines. We also identified 7,531 transcripts not associated with any horse gene annotated in public databases. Of these, 3,280 transcripts did not have a homologous match to any sequence deposited in the NCBI EST database suggesting horse specificity. The unknown transcripts were categorized as coding and noncoding based on predicted coding potential scores. Among them 230 transcripts had high coding potential score, at least 2 exons, and an open reading frame of at least 300 nt. We experimentally validated 9 new equine coding transcripts using RT-PCR and Sanger sequencing. Our results provide valuable detailed information on many transcripts yet to be annotated in the horse genome.

In search of epigenetic marks in testes and sperm cells of differentially fed boars

In search of transmittable epigenetic marks we investigated gene expression in testes and sperm cells of differentially fed F0 boars from a three generation pig feeding experiment that showed phenotypic differences in the F2 generation. RNA samples from 8 testes of boars that received either a diet enriched in methylating micronutrients or a control diet were analyzed by microarray analysis. We found moderate differential expression between testes of differentially fed boars with a high FDR of 0.82 indicating that most of the differentially expressed genes were false positives. Nevertheless, we performed a pathway analysis and found disparate pathway maps of development_A2B receptor: action via G-protein alpha s, cell adhesion_Tight junctions and cell adhesion_Endothelial cell contacts by junctional mechanisms which show inconclusive relation to epigenetic inheritance. Four RNA samples from sperm cells of these differentially fed boars were analyzed by RNA-Seq methodology. We found no differential gene expression in sperm cells of the two groups (adjusted P-value>0.05). Nevertheless, we also explored gene expression in sperm by a pathway analysis showing that genes were enriched for the pathway maps of bacterial infections in cystic fibrosis (CF) airways, glycolysis and gluconeogenesis p.3 and cell cycle_Initiation of mitosis. Again, these pathway maps are miscellaneous without an obvious relationship to epigenetic inheritance. It is concluded that the methylating micronutrients moderately if at all affects RNA expression in testes of differentially fed boars. Furthermore, gene expression in sperm cells is not significantly affected by extensive supplementation of methylating micronutrients and thus RNA molecules could not be established as the epigenetic mark in this feeding experiment.

NMD3 regulates both mRNA and rRNA nuclear export in African trypanosomes via an XPOI-linked pathway

Trypanosomes mostly regulate gene expression through post-transcriptional mechanisms, particularly mRNA stability. However, much mRNA degradation is cytoplasmic such that mRNA nuclear export must represent an important level of regulation. Ribosomal RNAs must also be exported from the nucleus and the trypanosome orthologue of NMD3 has been confirmed to be involved in rRNA processing and export, matching its function in other organisms. Surprisingly, we found that TbNMD3 depletion also generates mRNA accumulation of procyclin-associated genes (PAGs), these being co-transcribed by RNA polymerase I with the procyclin surface antigen genes expressed on trypanosome insect forms. By whole transcriptome RNA-seq analysis of TbNMD3-depleted cells we confirm the regulation of the PAG transcripts by TbNMD3 and using reporter constructs reveal that PAG1 regulation is mediated by its 5'UTR. Dissection of the mechanism of regulation demonstrates that it is not dependent upon translational inhibition mediated by TbNMD3 depletion nor enhanced transcription. However, depletion of the nuclear export factors XPO1 or MEX67 recapitulates the effects of TbNMD3 depletion on PAG mRNAs and mRNAs accumulated in the nucleus of TbNMD3-depleted cells. These results invoke a novel RNA regulatory mechanism involving the NMD3-dependent nuclear export of mRNA cargos, suggesting a shared platform for mRNA and rRNA export.

Impaired Cell Cycle Regulation in a Natural Equine Model of Asthma.

Recurrent airway obstruction (RAO) is a common and potentially debilitating lower airway disease in horses, which shares many similarities with human asthma. In susceptible horses RAO exacerbation is caused by environmental allergens and irritants present in hay dust. The objective of this study was the identification of genes and pathways involved in the pathology of RAO by global transcriptome analyses in stimulated peripheral blood mononuclear cells (PBMCs). We performed RNA-seq on PBMCs derived from 40 RAO affected and 45 control horses belonging to three cohorts of Warmblood horses: two half-sib families and one group of unrelated horses. PBMCs were stimulated with hay dust extract, lipopolysaccharides, a recombinant parasite antigen, or left unstimulated. The total dataset consisted of 561 individual samples. We detected significant differences in the expression profiles between RAO and control horses. Differential expression (DE) was most marked upon stimulation with hay dust extract. An important novel finding was a strong upregulation of CXCL13 together with many genes involved in cell cycle regulation in stimulated samples from RAO affected horses, in addition to changes in the expression of several HIF-1 transcription factor target genes. The RAO condition alters systemic changes observed as differential expression profiles of PBMCs. Those changes also depended on the cohort and stimulation of the samples and were dominated by genes involved in immune cell trafficking, development, and cell cycle regulation. Our findings indicate an important role of CXCL13, likely macrophage or Th17 derived, and the cell cycle regulator CDC20 in the immune response in RAO.

Evidence for two protein coding transcripts at the Igf2as locus

The insulin-like growth factor 2 antisense (Igf2as) gene is part of the Ins-Igf2-H19 imprinted gene cluster. The function of the paternally expressed Igf2as is still elusive. In our previous work, we showed that Igf2as transcripts were located in the cytoplasm of C2C12 mouse myoblast cells, associated with polysomes and polyadenylated suggesting that Igf2as is protein coding. In the present work, the protein coding capacity of Igf2as was investigated. We demonstrate for the first time the existence of a polypeptide translated from an Igf2as construct. Furthermore, an RNA-Seq analysis was performed using RNA prepared from skeletal muscles of newborn wild-type and ∆ DMR1-U2 mice to further elucidate the function of Igf2as transcripts. We found no evidence for a regulatory role of Igf2as in the imprinted gene cluster. Interestingly, the RNA-Seq analysis indicated that Igf2as plays a role in the energy metabolism, the cell cycle, histone acetylation and muscle contraction pathways. Our Igf2as investigations further elucidated that there are two distinct Igf2as transcripts corresponding to two putative ORFs.

A sensitized RNA interference screen identifies a novel role for the PI3K p110γ isoform in medulloblastoma cell proliferation and chemoresistance

Medulloblastoma is the most common malignant brain tumor in children and is associated with a poor outcome. We were interested in gaining further insight into the potential of targeting the human kinome as a novel approach to sensitize medulloblastoma to chemotherapeutic agents. A library of small interfering RNA (siRNA) was used to downregulate the known human protein and lipid kinases in medulloblastoma cell lines. The analysis of cell proliferation, in the presence or absence of a low dose of cisplatin after siRNA transfection, identified new protein and lipid kinases involved in medulloblastoma chemoresistance. PLK1 (polo-like kinase 1) was identified as a kinase involved in proliferation in medulloblastoma cell lines. Moreover, a set of 6 genes comprising ATR, LYK5, MPP2, PIK3CG, PIK4CA, and WNK4 were identified as contributing to both cell proliferation and resistance to cisplatin treatment in medulloblastoma cells. An analysis of the expression of the 6 target genes in primary medulloblastoma tumor samples and cell lines revealed overexpression of LYK5 and PIK3CG. The results of the siRNA screen were validated by target inhibition with specific pharmacological inhibitors. A pharmacological inhibitor of p110γ (encoded by PIK3CG) impaired cell proliferation in medulloblastoma cell lines and sensitized the cells to cisplatin treatment. Together, our data show that the p110γ phosphoinositide 3-kinase isoform is a novel target for combinatorial therapies in medulloblastoma.

Investigation of a unique short open reading frame within the 3' untranslated region of the canine distemper virus matrix messenger RNA

Increasing evidence suggest that the long "untranslated" region (UTR) between the matrix (M) and the fusion (F) proteins of morbilliviruses has a functional role. In canine distemper virus (CDV), the F 5' UTR was recently shown to code for a long F signal peptide (Fsp). Subsequently, it was reported that the M/F UTRs combined with the long Fsp were synergistically regulating the F mRNA and protein expression, thereby modulating virulence. Unique to CDV, a short putative open reading frame (ORF) has been identified within the wild-type CDV-M 3' UTR (termed M2). Here, we investigated whether M2 was expressed from the genome of the virulent and demyelinating A75/17-CDV strain. An expression plasmid encoding the M2 ORF tagged both at its N-terminal (HA) and C-terminal domains (RFP), was first constructed. Then, a recombinant virus with its putative M2 ORF replaced by HA-M2-RFP was successfully recovered from cDNA (termed recA75/17(green)-HA-M2-RFP). M2 expression in cells transfected or infected with these mutants was studied by immunoprecipitation, immunofluorescence, immunoblot and flow cytometry analyses. Although fluorescence was readily detected in HA-M2-RFP-transfected cells, absence of red fluorescence emission in several recA75/17(green)-HA-M2-RFP-infected cell types suggested lack of M2 biosynthesis, which was confirmed by the other techniques. Consistent with these data, no functional role of the short polypeptide was revealed by infecting various cell types with HA-M2-RFP over-expressing or M2-knockout recombinant viruses. Thus, in sharp contrast to the CDV-F 5' UTR reported to translate a long Fsp, our data provided evidence that the CDV-M 3' UTR does not express any polypeptides.

Generation and analysis of a mouse intestinal metatranscriptome through Illumina based RNA-sequencing

With the advent of high through-put sequencing (HTS), the emerging science of metagenomics is transforming our understanding of the relationships of microbial communities with their environments. While metagenomics aims to catalogue the genes present in a sample through assessing which genes are actively expressed, metatranscriptomics can provide a mechanistic understanding of community inter-relationships. To achieve these goals, several challenges need to be addressed from sample preparation to sequence processing, statistical analysis and functional annotation. Here we use an inbred non-obese diabetic (NOD) mouse model in which germ-free animals were colonized with a defined mixture of eight commensal bacteria, to explore methods of RNA extraction and to develop a pipeline for the generation and analysis of metatranscriptomic data. Applying the Illumina HTS platform, we sequenced 12 NOD cecal samples prepared using multiple RNA-extraction protocols. The absence of a complete set of reference genomes necessitated a peptide-based search strategy. Up to 16% of sequence reads could be matched to a known bacterial gene. Phylogenetic analysis of the mapped ORFs revealed a distribution consistent with ribosomal RNA, the majority from Bacteroides or Clostridium species. To place these HTS data within a systems context, we mapped the relative abundance of corresponding Escherichia coli homologs onto metabolic and protein-protein interaction networks. These maps identified bacterial processes with components that were well-represented in the datasets. In summary this study highlights the potential of exploiting the economy of HTS platforms for metatranscriptomics.

Base pairing and miscoding properties of 1,N(6)-ethenoadenine- and 3,N(4)-ethenocytosine-containing RNA oligonucleotides

Two RNA phosphoramidites containing the bases 1,N(6)-ethenoadenine (εA) and 3,N(4)-ethenocytosine (εC) were synthesized. These building blocks were incorporated into two 12-mer oligoribonucleotides for evaluation of the base pairing properties of these base lesions by UV melting curve (Tm) and circular dichroism measurements. The Tm data of the resulting duplexes with the etheno modifications opposing all natural bases showed a substantial destabilization compared to the corresponding natural duplexes, confirming their inability to form base pairs. The coding properties of these lesions were further investigated by introducing them into 31-mer oligonucleotides and assessing their ability to serve as templates in primer extension reactions with HIV, AMV, and MMLV reverse transcriptases (RT). Primer extension reactions showed complete arrest of the incorporation process using MMLV RT and AMV RT, while HIV RT preferentially incorporates dAMP opposite εA and dAMP as well as dTMP opposite εC. The properties of these RNA lesions are discussed in the context of its putative biological role.