Combining DGE and RNA-sequencing data to identify new polyA+ non-coding transcripts in the human genome.

Recent sequencing technologies that allow massive parallel production of short reads are the method of choice for transcriptome analysis. Particularly, digital gene expression (DGE) technologies produce a large dynamic range of expression data by generating short tag signatures for each cell transcript. These tags can be mapped back to a reference genome to identify new transcribed regions that can be further covered by RNA-sequencing (RNA-Seq) reads. Here, we applied an integrated bioinformatics approach that combines DGE tags, RNA-Seq, tiling array expression data and species-comparison to explore new transcriptional regions and their specific biological features, particularly tissue expression or conservation. We analysed tags from a large DGE data set (designated as 'TranscriRef'). We then annotated 750,000 tags that were uniquely mapped to the human genome according to Ensembl. We retained transcripts originating from both DNA strands and categorized tags corresponding to protein-coding genes, antisense, intronic- or intergenic-transcribed regions and computed their overlap with annotated non-coding transcripts. Using this bioinformatics approach, we identified ∼34,000 novel transcribed regions located outside the boundaries of known protein-coding genes. As demonstrated using sequencing data from human pluripotent stem cells for biological validation, the method could be easily applied for the selection of tissue-specific candidate transcripts. DigitagCT is available at http://cractools.gforge.inria.fr/softwares/digitagct.

New prognostic markers, determined using gene expression analyses, reveal two distinct subtypes of chronic myelomonocytic leukaemia patients.

Chronic myelomonocytic leukaemia (CMML) is a heterogeneous haematopoietic disorder characterized by myeloproliferative or myelodysplastic features. At present, the pathogenesis of this malignancy is not completely understood. In this study, we sought to analyse gene expression profiles of CMML in order to characterize new molecular outcome predictors. A learning set of 32 untreated CMML patients at diagnosis was available for TaqMan low-density array gene expression analysis. From 93 selected genes related to cancer and cell cycle, we built a five-gene prognostic index after multiplicity correction. Using this index, we characterized two categories of patients with distinct overall survival (94% vs. 19% for good and poor overall survival, respectively; P = 0·007) and we successfully validated its strength on an independent cohort of 21 CMML patients with Affymetrix gene expression data. We found no specific patterns of association with traditional prognostic stratification parameters in the learning cohort. However, the poor survival group strongly correlated with high-risk treated patients and transformation to acute myeloid leukaemia. We report here a new multigene prognostic index for CMML, independent of the gene expression measurement method, which could be used as a powerful tool to predict clinical outcome and help physicians to evaluate criteria for treatments.

A decrease in NAMPT activity impairs basal PARP-1 activity in cytidine deaminase deficient-cells, independently of NAD.

Cytidine deaminase (CDA) deficiency causes pyrimidine pool disequilibrium. We previously reported that the excess cellular dC and dCTP resulting from CDA deficiency jeopardizes genome stability, decreasing basal poly(ADP-ribose) polymerase 1 (PARP-1) activity and increasing ultrafine anaphase bridge (UFB) formation. Here, we investigated the mechanism underlying the decrease in PARP-1 activity in CDA-deficient cells. PARP-1 activity is dependent on intracellular NAD+ concentration. We therefore hypothesized that defects of the NAD+ salvage pathway might result in decreases in PARP-1 activity. We found that the inhibition or depletion of nicotinamide phosphoribosyltransferase (NAMPT), the rate-limiting enzyme in the NAD+ salvage biosynthesis pathway, mimicked CDA deficiency, resulting in a decrease in basal PARP-1 activity, regardless of NAD+ levels. Furthermore, the expression of exogenous wild-type NAMPT fully restored basal PARP-1 activity and prevented the increase in UFB frequency in CDA-deficient cells. No such effect was observed with the catalytic mutant. Our findings demonstrate that (1) the inhibition of NAMPT activity in CDA-proficient cells lowers basal PARP-1 activity, and (2) the expression of exogenous wild-type NAMPT, but not of the catalytic mutant, fully restores basal PARP-1 activity in CDA-deficient cells; these results strongly suggest that basal PARP-1 activity in CDA-deficient cells decreases due to a reduction of NAMPT activity.

Bloom syndrome protein restrains innate immune sensing of micronuclei by cGAS.

Cellular innate immune sensors of DNA are essential for host defense against invading pathogens. However, the presence of self-DNA inside cells poses a risk of triggering unchecked immune responses. The mechanisms limiting induction of inflammation by self-DNA are poorly understood. BLM RecQ-like helicase is essential for genome integrity and is deficient in Bloom syndrome (BS), a rare genetic disease characterized by genome instability, accumulation of micronuclei, susceptibility to cancer, and immunodeficiency. Here, we show that BLM-deficient fibroblasts show constitutive up-regulation of inflammatory interferon-stimulated gene (ISG) expression, which is mediated by the cGAS-STING-IRF3 cytosolic DNA-sensing pathway. Increased DNA damage or down-regulation of the cytoplasmic exonuclease TREX1 enhances ISG expression in BLM-deficient fibroblasts. cGAS-containing cytoplasmic micronuclei are increased in BS cells. Finally, BS patients demonstrate elevated ISG expression in peripheral blood. These results reveal that BLM limits ISG induction, thus connecting DNA damage to cellular innate immune response, which may contribute to human pathogenesis.

A role for Tau protein in maintaining ribosomal DNA stability and cytidine deaminase-deficient cell survival.

Cells from Bloom's syndrome patients display genome instability due to a defective BLM and the downregulation of cytidine deaminase. Here, we use a genome-wide RNAi-synthetic lethal screen and transcriptomic profiling to identify genes enabling BLM-deficient and/or cytidine deaminase-deficient cells to tolerate constitutive DNA damage and replication stress. We found a synthetic lethal interaction between cytidine deaminase and microtubule-associated protein Tau deficiencies. Tau is overexpressed in cytidine deaminase-deficient cells, and its depletion worsens genome instability, compromising cell survival. Tau is recruited, along with upstream-binding factor, to ribosomal DNA loci. Tau downregulation decreases upstream binding factor recruitment, ribosomal RNA synthesis, ribonucleotide levels, and affects ribosomal DNA stability, leading to the formation of a new subclass of human ribosomal ultrafine anaphase bridges. We describe here Tau functions in maintaining survival of cytidine deaminase-deficient cells, and ribosomal DNA transcription and stability. Moreover, our findings for cancer tissues presenting concomitant cytidine deaminase underexpression and Tau upregulation open up new possibilities for anti-cancer treatment.Cytidine deaminase (CDA) deficiency leads to genome instability. Here the authors find a synthetic lethal interaction between CDA and the microtubule-associated protein Tau deficiencies, and report that Tau depletion affects rRNA synthesis, ribonucleotide pool balance, and rDNA stability.

New chimeric RNAs in acute myeloid leukemia.

Background: High-throughput next generation sequencing (NGS) technologies enable the detection of biomarkers used for tumor classification, disease monitoring and cancer therapy. Whole-transcriptome analysis using RNA-seq is important, not only as a means of understanding the mechanisms responsible for complex diseases but also to efficiently identify novel genes/exons, splice isoforms, RNA editing, allele-specific mutations, differential gene expression and fusion-transcripts or chimeric RNA (chRNA). Methods: We used Crac, a tool that uses genomic locations and local coverage to classify biological events and directly infer splice and chimeric junctions within a single read. Crac's algorithm extracts transcriptional chimeric events irrespective of annotation with a high sensitivity, and CracTools was used to aggregate, annotate and filter the chRNA reads. The selected chRNA candidates were validated by real time PCR and sequencing.  In order to check the tumor specific expression of chRNA, we analyzed a publicly available dataset using a new tag search approach. Results:  We present data related to acute myeloid leukemia (AML) RNA-seq analysis. We highlight novel biological cases of chRNA, in addition to previously well characterized leukemia chRNA. We have identified and validated 17 chRNAs among 3 AML patients: 10 from an AML patient with a translocation between chromosomes 15 and 17 (AML-t(15;17), 4  from patient with normal karyotype (AML-NK) 3 from a patient with chromosomal 16 inversion (AML-inv16). The new fusion transcripts can be classified into four groups according to the exon organization. Conclusions:  All groups suggest complex but distinct synthesis mechanisms involving either collinear exons of different genes, non-collinear exons, or exons of different chromosomes. Finally, we check tumor-specific expression in a larger RNA-seq AML cohort and identify new AML biomarkers that could improve diagnosis and prognosis of AML.

Identification of a 20-gene expression-based risk score as a predictor of clinical outcome in chronic lymphocytic leukemia patients.

Despite the improvement in treatment options, chronic lymphocytic leukemia (CLL) remains an incurable disease and patients show a heterogeneous clinical course requiring therapy for many of them. In the current work, we have built a 20-gene expression (GE)-based risk score predictive for patients overall survival and improving risk classification using microarray gene expression data. GE-based risk score allowed identifying a high-risk group associated with a significant shorter overall survival (OS) and time to treatment (TTT) (P ≤ .01), comprising 19.6% and 13.6% of the patients in two independent cohorts. GE-based risk score, and NRIP1 and TCF7 gene expression remained independent prognostic factors using multivariate Cox analyses and combination of GE-based risk score together with NRIP1 and TCF7 gene expression enabled the identification of three clinically distinct groups of CLL patients. Therefore, this GE-based risk score represents a powerful tool for risk stratification and outcome prediction of CLL patients and could thus be used to guide clinical and therapeutic decisions prospectively.

Development of gene expression-based risk score in cytogenetically normal acute myeloid leukemia patients.

Patients with normal karyotype represent the single largest cytogenetic group of acute myeloid leukemia (AML), with highly heterogeneous clinical and molecular characteristics. In this study, we sought to determine new prognostic biomarkers in cytogenetically normal (CN)-AML patients. A gene expression (GE)-based risk score was built, summing up the prognostic value of 22 genes whose expression is associated with a bad prognosis in a training cohort of 163 patients. GE-based risk score allowed identifying a high-risk group of patients (53.4%) in two independent cohorts of CN-AML patients. GE-based risk score and EVI1 gene expression remained independent prognostic factors using multivariate Cox analyses. Combining GE-based risk score with EVI1 gene expression allowed the identification of three clinically different groups of patients in two independent cohorts of CN-AML patients. Thus, GE-based risk score is powerful to predict clinical outcome for CN-AML patients and may provide potential therapeutic advances.