Modulation of Yorkie activity by alternative splicing is required for developmental stability.

The Hippo signaling pathway is a major regulator of organ growth, which controls the activity of the transcription coactivator Yorkie (Yki) in Drosophila and its homolog YAP in mammals. Both Yki and YAP proteins exist as alternatively spliced isoforms containing either one or two WW domains. The biological importance of this conserved alternative splicing event is unknown. Here, we identify the splicing factor B52 as a regulator of yki alternative splicing in Drosophila and show that B52 modulates growth in part through modulation of yki alternative splicing. Yki isoforms differ by their transcriptional activity as well as their ability to bind and bridge PPxY motifs-containing partners, and can compete in vivo. Strikingly, flies in which yki alternative splicing has been abrogated, thus expressing only Yki2 isoform, exhibit fluctuating wing asymmetry, a signal of developmental instability. Our results identify yki alternative splicing as a new level of modulation of the Hippo pathway, that is required for growth equilibration during development. This study provides the first demonstration that the process of alternative splicing contributes to developmental robustness.

The core spliceosomal factor U2AF1 controls cell-fate determination via the modulation of transcriptional networks.

Alternative splicing (AS) plays a central role during cell-fate determination. However, how the core spliceosomal factors (CSFs) are involved in this process is poorly understood. Here, we report the down-regulation of the U2AF1 CSF during stem cell differentiation. To investigate its function in stemness and differentiation, we downregulated U2AF1 in human induced pluripotent stem cells (hiPSCs), using an inducible-shRNA system, to the level found in differentiated ectodermal, mesodermal and endodermal cells. RNA sequencing and computational analysis reveal that U2AF1 down-regulation modulates the expression of development-regulating genes and regulates transcriptional networks involved in cell-fate determination. Furthermore, U2AF1 down-regulation induces a switch in the AS of transcription factors (TFs) required to establish specific cell lineages, and favours the splicing of a differentiated cell-specific isoform of DNMT3B. Our results showed that the differential expression of the core spliceosomal factor U2AF1, between stem cells and the precursors of the three germ layers regulates a cell-type-specific alternative splicing programme and a transcriptional network involved in cell-fate determination via the modulation of gene expression and alternative splicing of transcription regulators.

Both anti-inflammatory and antiviral properties of novel drug candidate ABX464 are mediated by modulation of RNA splicing.

ABX464 is a first-in-class, clinical-stage, small molecule for oral administration that has shown strong anti-inflammatory effects in the DSS-model for inflammatory bowel disease (IBD) and also prevents replication of the HIV virus. ABX464 which binds to cap binding complex (CBC) has demonstrated safety and efficacy in a phase 2a proof-of-concept clinical trial in patients with Ulcerative colitis. Previously, with limited technologies, it was not possible to quantify the effect of ABX464 on viral and cellular RNA biogenesis. Here, using RNA CaptureSeq and deep sequencing, we report that ABX464 enhances the splicing of HIV RNA in infected PBMCs from six healthy individuals and also the expression and splicing of a single long noncoding RNA to generate the anti-inflammatory miR-124 both ex vivo and in HIV patients. While ABX464 has no effect on pre-mRNA splicing of cellular genes, depletion of CBC complex by RNAi leads to accumulation of intron retention transcripts. These results imply that ABX464 did not inhibit the function of CBC in splicing but rather strengthens it under pathological condition like inflammation and HIV infection. The specific dual ability of ABX464 to generate both anti-inflammatory miR-124 and spliced viral RNA may have applicability for the treatment of both inflammatory diseases and HIV infection.

RNA sequencing analysis of activated macrophages treated with the anti-HIV ABX464 in intestinal inflammation.

RNA-Seq enables the generation of extensive transcriptome information providing the capability to characterize transcripts (including alternative isoforms and polymorphism), to quantify expression and to identify differential regulation in a single experiment. To reveal the capacity of new anti-HIV ABX464 candidate in modulating the expression of genes, datasets were generated and validated using RNA-seq approach. This comprehensive dataset will be useful to deepen the comprehensive understanding of the progression of human immunodeficiency virus (HIV) associated with mucosal damage in the gastrointestinal (GI) tract and subsequent inflammation, providing an opportunity to generate new therapies, diagnoses, and preventive strategies.

The intellectual disability of trisomy 21: differences in gene expression in a case series of patients with lower and higher IQ.

Trisomy 21 (T21), or Down syndrome (DS), is the most frequent and recognizable cause of intellectual disabilities. The level of disability, as evaluated by the intelligence quotient (IQ) test, varies considerably between patients independent of other factors. To determine the genetic or molecular basis of this difference, a high throughput transcriptomic analysis was performed on twenty T21 patients with high and low IQ, and 10 healthy controls using Digital Gene Expression. More than 90 millions of tags were sequenced in the three libraries. A total of 80 genes of potential interest were selected for the qPCR experiment validation, and three housekeeping genes were used for normalizing purposes. HLA DQA1 and HLA DRB1 were significantly downregulated among the patients with a low IQ, the values found in the healthy controls being intermediate between those noted in the IQ+ and IQ- T21 patients. Interestingly, the intergenic region between these genes contains a binding sequence for the CCCTC-binding factor, or CTCF, and cohesin (a multisubunit complex), both of which are essential for expression of HLA DQA1 and HLA DRB1 and numerous other genes. Our results might lead to the discovery of genes, or genetic markers, that are directly involved in several phenotypes of DS and, eventually, to the identification of potential targets for therapeutic interventions.

Transcriptome and proteome analysis of Pinctada margaritifera calcifying mantle and shell: focus on biomineralization.

BACKGROUND: The shell of the pearl-producing bivalve Pinctada margaritifera is composed of an organic cell-free matrix that plays a key role in the dynamic process of biologically-controlled biomineralization. In order to increase genomic resources and identify shell matrix proteins implicated in biomineralization in P. margaritifera, high-throughput Expressed Sequence Tag (EST) pyrosequencing was undertaken on the calcifying mantle, combined with a proteomic analysis of the shell. RESULTS: We report the functional analysis of 276 738 sequences, leading to the constitution of an unprecedented catalog of 82 P. margaritifera biomineralization-related mantle protein sequences. Components of the current "chitin-silk fibroin gel-acidic macromolecule" model of biomineralization processes were found, in particular a homolog of a biomineralization protein (Pif-177) recently discovered in P. fucata. Among these sequences, we could show the localization of two other biomineralization protein transcripts, pmarg-aspein and pmarg-pearlin, in two distinct areas of the outer mantle epithelium, suggesting their implication in calcite and aragonite formation. Finally, by combining the EST approach with a proteomic mass spectrometry analysis of proteins isolated from the P. margaritifera shell organic matrix, we demonstrated the presence of 30 sequences containing almost all of the shell proteins that have been previously described from shell matrix protein analyses of the Pinctada genus. The integration of these two methods allowed the global composition of biomineralizing tissue and calcified structures to be examined in tandem for the first time. CONCLUSIONS: This EST study made on the calcifying tissue of P. margaritifera is the first description of pyrosequencing on a pearl-producing bivalve species. Our results provide direct evidence that our EST data set covers most of the diversity of the matrix protein of P. margaritifera shell, but also that the mantle transcripts encode proteins present in P. margaritifera shell, hence demonstrating their implication in shell formation. Combining transcriptomic and proteomic approaches is therefore a powerful way to identify proteins involved in biomineralization. Data generated in this study supply the most comprehensive list of biomineralization-related sequences presently available among protostomian species, and represent a major breakthrough in the field of molluskan biomineralization.

Blood cells RNA biomarkers as a first long-term detection strategy for EPO abuse in horseracing.

Recombinant human erythropoietins (rHuEPOs) are glycoproteins drugs, produced by the pharmaceutical industry to restore production of red blood cells by stimulating human bone marrow for which this pathology has been diagnosed. It is suspected that these molecules are diverted as doping agents in horseracing to enhance oxygen transport and aerobic power in racehorses. Although indirect double-blotting or direct liquid chromatography-mass spectrometry (LC-MS) methods have been developed to confirm the presence of rHuEPO in a sample, the short detection time (48 h) is still a problem for doping control. In this context, gene profiling investigation through Serial Analysis of Gene Expression (SAGE) has been conducted on seven thoroughbreds treated with Eprex. This functional genomic method has been performed from total blood cells collected from each animal to assess the mRNA expression consecutive to rHuEPO injections. Sample pooling was chosen as a powerful, cost-effective, and rapid means of identifying the most common and specific changes in terms of gene expression profile and to eliminate individual variation. Consequently, three SAGE libraries were constructed, before, during, and after Eprex treatment. More than 71 440 mRNA signatures were observed and subjected to statistical analysis; 49 differentially expressed genes were identified and analyzed by qPCR. From the selected gene list, were defined as potential biomarkers in terms of their low inter-individual variation and capacity as strong markers of rHuEPO administration up to 60 days after the beginning of the doping period. In this paper, a new strategy is proposed to the horseracing industry to prevent rHuEPO abuse.

Endoribonuclease L (RNase L) regulates the myogenic and adipogenic potential of myogenic cells.

Skeletal muscle maintenance and repair involve several finely coordinated steps in which pluripotent stem cells are activated, proliferate, exit the cell cycle and differentiate. This process is accompanied by activation of hundreds of muscle-specific genes and repression of genes associated with cell proliferation or pluripotency. Mechanisms controlling myogenesis are precisely coordinated and regulated in time to allow the sequence of activation/inactivation of genes expression. Muscular differentiation is the result of the interplay between several processes such as transcriptional induction, transcriptional repression and mRNA stability. mRNA stability is now recognized as an essential mechanism of control of gene expression. For instance, we previously showed that the endoribonuclease L (RNase L) and its inhibitor (RLI) regulates MyoD mRNA stability and consequently muscle differentiation.We now performed global gene expression analysis by SAGE to identify genes that were down-regulated upon activation of RNase L in C2C12 myogenic cells, a model of satellite cells. We found that RNase L regulates mRNA stability of factors implicated in the control of pluripotency and cell differentiation. Moreover, inappropriate RNase L expression in C2C12 cells led to inhibition of myogenesis and differentiation into adipocytes even when cells were grown in conditions permissive for muscle differentiation. Conversely, over-expression of RLI allowed muscle differentiation of myogenic C2C12 cells even in non permissive conditions.These findings reveal the central role of RNase L and RLI in controlling gene expression and cell fate during myogenesis. Our data should provide valuable insights into the mechanisms that control muscle stem cell differentiation and into the mechanism of metaplasia observed in aging or muscular dystrophy where adipose infiltration of muscle occurs.

Simultaneous gene expression profiling in human macrophages infected with Leishmania major parasites using SAGE.

BACKGROUND: Leishmania (L) are intracellular protozoan parasites that are able to survive and replicate within the harsh and potentially hostile phagolysosomal environment of mammalian mononuclear phagocytes. A complex interplay then takes place between the macrophage (MPhi) striving to eliminate the pathogen and the parasite struggling for its own survival. To investigate this host-parasite conflict at the transcriptional level, in the context of monocyte-derived human MPhis (MDM) infection by L. major metacyclic promastigotes, the quantitative technique of serial analysis of gene expression (SAGE) was used. RESULTS: After extracting mRNA from resting human MPhis, Leishmania-infected human MPhis and L. major parasites, three SAGE libraries were constructed and sequenced generating up to 28,173; 57,514 and 33,906 tags respectively (corresponding to 12,946; 23,442 and 9,530 unique tags). Using computational data analysis and direct comparison to 357,888 publicly available experimental human tags, the parasite and the host cell transcriptomes were then simultaneously characterized from the mixed cellular extract, confidently discriminating host from parasite transcripts. This procedure led us to reliably assign 3,814 tags to MPhis' and 3,666 tags to L. major parasites transcripts. We focused on these, showing significant changes in their expression that are likely to be relevant to the pathogenesis of parasite infection: (i) human MPhis genes, belonging to key immune response proteins (e.g., IFNgamma pathway, S100 and chemokine families) and (ii) a group of Leishmania genes showing a preferential expression at the parasite's intra-cellular developing stage. CONCLUSION: Dual SAGE transcriptome analysis provided a useful, powerful and accurate approach to discriminating genes of human or parasitic origin in Leishmania-infected human MPhis. The findings presented in this work suggest that the Leishmania parasite modulates key transcripts in human MPhis that may be beneficial for its establishment and survival. Furthermore, these results provide an overview of gene expression at two developmental stages of the parasite, namely metacyclic promastigotes and intracellular amastigotes and indicate a broad difference between their transcriptomic profiles. Finally, our reported set of expressed genes will be useful in future rounds of data mining and gene annotation.

Pharmacogenomic analysis of acute promyelocytic leukemia cells highlights CYP26 cytochrome metabolism in differential all-trans retinoic acid sensitivity.

Disease relapse sometimes occurs after acute promyelocytic leukemia (APL) therapy with all-trans retinoic acid (ATRA). Among the diagnostic parameters predicting relapse, heterogeneity in the in vitro differentiation rate of blasts is an independent factor. To identify biologic networks involved in resistance, we conducted pharmacogenomic studies in APL blasts displaying distinct ATRA sensitivities. Although the expression profiles of genes invested in differentiation were similarly modulated in low- and high-sensitive blasts, low-sensitive cells showed higher levels of transcription of ATRA-target genes, transcriptional regulators, chromatin remodelers, and transcription factors. In opposition, only high-sensitive blasts expressed the CYP26A1 gene, encoding the p450 cytochrome which is known to be involved in retinoic acid catabolism. In NB4 cells, ATRA treatment activates a novel signaling pathway, whereby interleukin-8 stimulates the expression of the homeobox transcription factor HOXA10v2, an effective enhancer of CYP26A1 transcription. These data were corroborated in primary APL cells, as maturation levels correlated with CYP26A1 expression. Treatment with a retinoic acid metabolism blocking agent (RAMBA) results in high-nucleoplasmic concentrations of retinoid and growth of NB4-resistant subclones. Hence, for APL blasts associated with poor prognosis, the low CYP26A1 expression may explain high risk of resistance installation, by increased retinoid pressure. Pharmacogenomic profiles of genes involved in retinoid acid metabolism may help to optimize anticancer therapies, including retinoids.

Use of the Serial Analysis of Gene Expression (SAGE) method in veterinary research: A concrete application in the study of the bovine trypanotolerance genetic control.

New postgenomic biotechnologies, such as transcriptome analyses, are now able to characterize the full complement of genes involved in the expression of specific biological functions. One of these is the Serial Analysis of Gene Expression (SAGE) technique, which consists of the construction of transcripts libraries for a quantitative analysis of the entire gene(s) expressed or inactivated at a particular step of cellular activation. Bioinformatic comparisons in the bovine genomic databases allow the identification of several up- and downregulated genes, expressed sequence tags, and unknown functional genes directly involved in the genetic control of the studied biological mechanism. We present and discuss the preliminary results in comparing the expressed genes in two total mRNA transcripts libraries obtained during an experimental Trypanosoma congolense infection in one trypanotolerant N'Dama animal cow. Knowing all the functional genes involved in the trypanotolerance control will permit validation of some results obtained with the quantitative trait locus approach, to set up specific microarrays sets for further metabolic and pharmacological studies, and to design field marker-assisted selection by introgression programs.

Mining SAGE data allows large-scale, sensitive screening of antisense transcript expression.

As a growing number of complementary transcripts, susceptible to exert various regulatory functions, are being found in eukaryotes, high throughput analytical methods are needed to investigate their expression in multiple biological samples. Serial Analysis of Gene Expression (SAGE), based on the enumeration of directionally reliable short cDNA sequences (tags), is capable of revealing antisense transcripts. We initially detected them by observing tags that mapped on to the reverse complement of known mRNAs. The presence of such tags in individual SAGE libraries suggested that SAGE datasets contain latent information on antisense transcripts. We raised a collection of virtual tags for mining these data. Tag pairs were assembled by searching for complementarities between 24-nt long sequences centered on the potential SAGE-anchoring sites of well-annotated human expressed sequences. An analysis of their presence in a large collection of published SAGE libraries revealed transcripts expressed at high levels from both strands of two adjacent, oppositely oriented, transcription units. In other cases, the respective transcripts of such cis-oriented genes displayed a mutually exclusive expression pattern or were co-expressed in a small number of libraries. Other tag pairs revealed overlapping transcripts of trans-encoded unique genes. Finally, we isolated a group of tags shared by multiple transcripts. Most of them mapped on to retroelements, essentially represented in humans by Alu sequences inserted in opposite orientations in the 3'UTR of otherwise different mRNAs. Registering these tags in separate files makes possible computational searches focused on unique sense-antisense pairs. The method developed in the present work shows that SAGE datasets constitute a major resource of rapidly investigating with high sensitivity the expression of antisense transcripts, so that a single tag may be detected in one library when screening a large number of biological samples.

Transcriptome analysis of monocytic leukemia cell differentiation.

The human leukemia cell line U937 is a well-established model for studying monocytic cell differentiation. We used a modified protocol (SADE) of serial analysis of gene expression (SAGE) and developed a SADE linker-anchored PCR assay to investigate the pattern of expression of known genes and to identify new transcripts in proliferating cells and during cell growth arrest and differentiation. We implemented new informatic tools to compare expression profiles before and after exposure of cells to differentiation inducers. From the analysis of 47,388 tags, we identified 13,806 distinct transcripts, 265 of which showed significant variations (P<0.01). Among 1219 well-identified genes, major changes concerned transcription and translation components, cytoskeleton, and macrophage-specific genes. Nearly half of the tags, some of them expressed at high levels, matched partially characterized genes or ESTs, or revealed yet-unknown transcripts, providing a wealth of new candidate genes that may reveal novel aspects of terminal monocytic differentiation.

[Transcriptomes for serial analysis of gene expression]. / Etude des transcriptomes par analyse sérielle de l'expression des gènes.

The availability of the sequences for whole genomes is changing our understanding of cell biology. Functional genomics refers to the comprehensive analysis, at the protein level (proteome) and at the mRNA level (transcriptome) of all events associated with the expression of whole sets of genes. New methods have been developed for transcriptome analysis. Serial Analysis of Gene Expression (SAGE) is based on the massive sequential analysis of short cDNA sequence tags. Each tag is derived from a defined position within a transcript. Its size (14 bp) is sufficient to identify the corresponding gene and the number of times each tag is observed provides an accurate measurement of its expression level. Since tag populations can be widely amplified without altering their relative proportions, SAGE may be performed with minute amounts of biological extract. Dealing with the mass of data generated by SAGE necessitates computer analysis. A software is required to automatically detect and count tags from sequence files. Criterias allowing to assess the quality of experimental data can be included at this stage. To identify the corresponding genes, a database is created registering all virtual tags susceptible to be observed, based on the present status of the genome knowledge. By using currently available database functions, it is easy to match experimental and virtual tags, thus generating a new database registering identified tags, together with their expression levels. As an open system, SAGE is able to reveal new, yet unknown, transcripts. Their identification will become increasingly easier with the progress of genome annotation. However, their direct characterization can be attempted, since tag information may be sufficient to design primers allowing to extend unknown sequences. A major advantage of SAGE is that, by measuring expression levels without reference to an arbitrary standard, data are definitively acquired and cumulative. All publicly available data can thus be stored in a unique database, facilitating whole-genome analysis of differential expression between cell types, normal and diseased samples, or samples with and without drug treatment. SAGE data are readily amenable to statistical comparisons, allowing to determine the level of confidence of the observed variations. A major limitation of SAGE is that, because each analysis is obligatory performed on the whole set of expressed genes, it can hardly be performed on multiple samples, for example in kinetics studies or to compare the effects of large numbers of drugs. To overcome this limitation, high-throughput detection of a subset of mRNAs is more rapidly performed by parallel hybridization of mRNAs on arrays of nucleic acids immobilized on solid supports. From this point of view, a SAGE platform is a powerful instrument for selecting the most informative subset of genes, assembling them to design microarrays dedicated to a specific problem and calibrating measurement by comparison with a standard cell model for which SAGE data are available. This approach is an attractive alternative to strategies based exclusively on pangenomic arrays. A very large amount of SAGE data are already available and the problem is now to extract their biological meaning. Knowledge on metabolic pathways is already organized so that its successful integration in a SAGE platform can be undertaken. For other cell components and pathways, the problem lies on the lack of controlled vocabulary to describe gene activities, starting form a clear definition of the concept of biological function itself. Progress in gene and cell ontology is expected to facilitate computer-based extraction of biological knowledge from existing and forthcoming SAGE data.