New computational biology tools for the systematic analysis of the structure and expression of human genes

From the late 1980s, the automation of sequencing techniques and the computer spread gave rise to a flourishing number of new molecular structures and sequences and to proliferation of new databases in which to store them. Here are presented three computational approaches able to analyse the massive amount of publicly avalilable data in order to answer to important biological questions. The first strategy studies the incorrect assignment of the first AUG codon in a messenger RNA (mRNA), due to the incomplete determination of its 5' end sequence. An extension of the mRNA 5' coding region was identified in 477 in human loci, out of all human known mRNAs analysed, using an automated expressed sequence tag (EST)-based approach. Proof-of-concept confirmation was obtained by in vitro cloning and sequencing for GNB2L1, QARS and TDP2 and the consequences for the functional studies are discussed. The second approach analyses the codon bias, the phenomenon in which distinct synonymous codons are used with different frequencies, and, following integration with a gene expression profile, estimates the total number of codons present across all the expressed mRNAs (named here "codonome value") in a given biological condition. Systematic analyses across different pathological and normal human tissues and multiple species shows a surprisingly tight correlation between the codon bias and the codonome bias. The third approach is useful to studies the expression of human autism spectrum disorder (ASD) implicated genes. ASD implicated genes sharing microRNA response elements (MREs) for the same microRNA are co-expressed in brain samples from healthy and ASD affected individuals. The different expression of a recently identified long non coding RNA which have four MREs for the same microRNA could disrupt the equilibrium in this network, but further analyses and experiments are needed.

The handedness-associated PCSK6 locus spans an intronic promoter regulating novel transcripts

We recently reported the association of the PCSK6 gene with handedness through a quantitative genome-wide association study (GWAS; P < 0.5 × 10(-8)) for a relative hand skill measure in individuals with dyslexia. PCSK6 activates Nodal, a morphogen involved in regulating left-right body axis determination. Therefore, the GWAS data suggest that the biology underlying the patterning of structural asymmetries may also contribute to behavioural laterality, e.g. handedness. The association is further supported by an independent study reporting a variable number tandem repeat (VNTR) within the same PCSK6 locus to be associated with degree of handedness in a general population cohort. Here, we have conducted a functional analysis of the PCSK6 locus combining further genetic analysis, in silico predictions and molecular assays. We have shown that the previous GWAS signal was not tagging a VNTR effect, suggesting that the two markers have independent effects. We demonstrated experimentally that one of the top GWAS-associated markers, rs11855145, directly alters the binding site for a nuclear factor. Furthermore, we have shown that the predicted regulatory region adjacent to rs11855415 acts as a bidirectional promoter controlling the expression of novel RNA transcripts. These include both an antisense long non-coding RNA (lncRNA) and a short PCSK6 isoform predicted to be coding. This is the first molecular characterization of a handedness-associated locus that supports the role of common variants in non-coding sequences in influencing complex phenotypes through gene expression regulation.

La perturbation du locus Nr2f1-K12 entraine une différenciation gliale précoce dans un nouveau modèle murin de mégacôlon aganglionnaire

La maladie de Hirschsprung est une affection congénitale de la motilité intestinale caractérisée par un segment aganglionnaire dans le côlon terminal. Un criblage génétique par mutation insertionnelle aléatoire chez la souris nous a permis d’identifier la lignée transgénique Spot dont les homozygotes souffrent de mégacôlon aganglionnaire. L’analyse d’intestins d’embryons mutants a révélé une baisse de prolifération et un délai de migration des cellules de la crête neurale entériques (CCNe) progénitrices dus à leur différenciation gliale précoce, entrainant un défaut de colonisation de l’intestin et une aganglionose du côlon. Le séquençage du génome Spot indique que le transgène s’est inséré à l’intérieur du locus K12-Nr2f1 sur le chromosome 13, une région dépourvue de gènes préalablement associés à la maladie, perturbant également une séquence non-codante très conservée dans l’évolution. K12 est un gène d’ARN long non codant (ARNlnc) et antisens du gène Nr2f1, lui-même impliqué dans la gliogénèse du système nerveux central. Le séquençage du transcriptome des CCN a montré une surexpression de Nr2f1 et des formes courtes de K12 chez Spot et des essais luciférase ont révélé l’activité répressive de l’élément conservé. Nous avons observé l’expression de K12 dans les CCNe et sa localisation subcellulaire dans des zones transcriptionnellement actives du noyau. Avec l’émergence des ARNlnc régulateurs, ces données nous permettent de pointer deux nouveaux gènes candidats associés à une différenciation gliale prématurée du SNE menant au mégacôlon aganglionnaire, en supposant que la régulation de Nr2f1 se fait par son antisens, K12.

La perturbation du locus Nr2f1-K12 entraine une différenciation gliale précoce dans un nouveau modèle murin de mégacôlon aganglionnaire

La maladie de Hirschsprung est une affection congénitale de la motilité intestinale caractérisée par un segment aganglionnaire dans le côlon terminal. Un criblage génétique par mutation insertionnelle aléatoire chez la souris nous a permis d’identifier la lignée transgénique Spot dont les homozygotes souffrent de mégacôlon aganglionnaire. L’analyse d’intestins d’embryons mutants a révélé une baisse de prolifération et un délai de migration des cellules de la crête neurale entériques (CCNe) progénitrices dus à leur différenciation gliale précoce, entrainant un défaut de colonisation de l’intestin et une aganglionose du côlon. Le séquençage du génome Spot indique que le transgène s’est inséré à l’intérieur du locus K12-Nr2f1 sur le chromosome 13, une région dépourvue de gènes préalablement associés à la maladie, perturbant également une séquence non-codante très conservée dans l’évolution. K12 est un gène d’ARN long non codant (ARNlnc) et antisens du gène Nr2f1, lui-même impliqué dans la gliogénèse du système nerveux central. Le séquençage du transcriptome des CCN a montré une surexpression de Nr2f1 et des formes courtes de K12 chez Spot et des essais luciférase ont révélé l’activité répressive de l’élément conservé. Nous avons observé l’expression de K12 dans les CCNe et sa localisation subcellulaire dans des zones transcriptionnellement actives du noyau. Avec l’émergence des ARNlnc régulateurs, ces données nous permettent de pointer deux nouveaux gènes candidats associés à une différenciation gliale prématurée du SNE menant au mégacôlon aganglionnaire, en supposant que la régulation de Nr2f1 se fait par son antisens, K12.

TAZ role in lung cancer: resistance to BET inhibitors and functional interaction with lncRNA

Despite extensive research and introduction of innovative therapy, lung cancer prognosis remains poor, with a five years survival of only 17%. The success of pharmacological treatment is often impaired by drug resistance. Thus, the characterization of response mechanisms to anti-cancer compounds and of the molecular mechanisms supporting lung cancer aggressiveness are crucial for patient’s management. In the first part of this thesis, we characterized the molecular mechanism behind resistance of lung cancer cells to the Inhibitors of the Bromodomain and Extraterminal domain containing Proteins (BETi). Through a CRISPR/Cas9 screening we identified three Hippo Pathway members, LATS2, TAOK1 and NF2 as genes implicated in susceptibility to BETi. These genes confer sensitivity to BETi inhibiting TAZ activity. Conversely, TAZ overexpression increases resistance to BETi. We also displayed that BETi downregulate both YAP, TAZ and TEADs expression in several cancer cell lines, implying a novel BETi-dependent cytotoxic mechanism. In the second part of this work, we attempted to characterize the crosstalk between the TAZ gene and its cognate antisense long-non coding RNA (lncRNA) TAZ-AS202 in lung tumorigenesis. As for TAZ downregulation, TAZ-AS202 silencing impairs NSCLC cells proliferation, migration and invasion, suggesting a pro-tumorigenic function for this lncRNA during lung tumorigenesis. TAZ-AS202 regulates TAZ target genes without altering TAZ expression or localization. This finding implies an uncovered functional cooperation between TAZ and TAZ-AS202. Moreover, we found that the EPH-ephrin signaling receptor EPHB2 is a downstream effector affected by both TAZ and TAZ-AS202 silencing. EPHB2 downregulation significantly attenuates cells proliferation, migration and invasion, suggesting that, at least in part, TAZ-AS202 and TAZ pro-oncogenic activity depends on EPH-ephrin signaling final deregulation. Finally, we started to dissect the mechanism underlying the TAZ-AS202 regulatory activity on EPHB2 in lung cancer, which may involve the existence of an intermediate transcription factor and is the object of our ongoing research.

When small non-coding RNAs meet the ribosome: Tuning the translational machinery

The functions of ribosomes in translation are complex and involve different types of activities critical for decoding the genetic code, linkage of amino acids via amide bonds to form polypeptide chains, as well as the release and proper targeting of the synthesized protein. Non-protein-coding RNAs (ncRNAs) have been recognized to be crucial in establishing regulatory networks.1 However all of the recently discovered ncRNAs involved in translation regulation target the mRNA rather than the ribosome. The main goal of this project is to identify potential novel ncRNAs that directly bind and possibly regulate the ribosome during protein biosynthesis. To address this question we applied various stress conditions to the archaeal model organism Haloferax volcanii and deep-sequenced the ribosome-associated small ncRNA interactome. In total we identified 6.250 ncRNA candidates. Significantly, we observed the emersed presence of tRNA-derived fragments (tRFs). These tRFs have been identified in all domains of life and represent a growing, yet functionally poorly understood, class of ncRNAs. Here we present evidence that tRFs from H. volcanii directly bind to ribosomes. In the presented genomic screen of the ribosome-associated RNome a 26 residue long fragment originating from the 5’ part of valine tRNA was by far the most abundant tRF. The Val-tRF is processed in a stress- dependent manner and was found to primarily target the small ribosomal subunit in vitro and in vivo. As a consequence of ribosome binding, Val-tRF reduces protein synthesis by interfering with peptidyl transferase activity. Therefore this tRF functions as ribosome-bound small ncRNA capable of regulating gene expression in H. volcanii under environmental stress conditions probably by fine-tuning the rate of protein production.2 Currently we are investigating the binding site of this tRF on the 30S subunit in more detail.

Post-transcriptional Regulation of Mammalian Gene Expression in Non-coding Region of Target RNA

Tumor Suppressor Candidate 2 (TUSC2) is a novel tumor suppressor gene located in the human chromosome 3p21.3 region. TUSC2 mRNA transcripts could be detected on Northern blots in both normal lung and some lung cancer cell lines, but no endogenous TUSC2 protein could be detected in a majority of lung cancer cell lines. Mechanisms regulating TUSC2 protein expression and its inactivation in primary lung cancer cells are largely unknown. We investigated the role of the 5’- and 3’-untranslated regions (UTRs) of the TUSC2 gene in the regulation of TUSC2 protein expression. We found that two small upstream open-reading frames (uORFs) in the 5’UTR of TUSC2 could markedly inhibit the translational initiation of TUSC2 protein by interfering with the “scanning” of the ribosome initiation complexes. Site-specific stem-loop array reverse transcription-polymerase chain reaction (SLA-RT-PCR) verified several micoRNAs (miRNAs) targeted at 3’UTR and directed TUSC2 cleavage and degradation. In addition, we used the established let-7-targeted high mobility group A2 (Hmga2) mRNA as a model system to study the mechanism of regulation of target mRNA by miRNAs in mammalian cells under physiological conditions. There have been no evidence of direct link between mRNA downregulation and mRNA cleavages mediated by miRNAs. Here we showed that the endonucleolytic cleavages on mRNAs were initiated by mammalian miRNA in seed pairing style. Let-7 directed cleavage activities among the eight predicted potential target sites have varied efficiency, which are influenced by the positional and the structural contexts in the UTR. The 5’ cleaved RNA fragments were mostly oligouridylated at their 3’-termini and accumulated for delayed 5’–3’ degradation. RNA fragment oligouridylation played important roles in marking RNA fragments for delayed bulk degradation and in converting RNA degradation mode from 3’–5’ to 5’–3’ with cooperative efforts from both endonucleolytic and non-catalytic miRNA-induced silencing complex (miRISC). Our findings point to a mammalian miRNA-mediated mechanism for the regulation of mRNA that miRNA can decrease target mRNA through target mRNA cleavage and uridine addition

Non-coding transcript in T cells (NTT): Antisense transcript activates PKR and NF-kappa B in human lymphocytes

T cell activation is a complex process involving many steps and the role played by the non-protein-coding RNAs (ncRNAs) in this phenomenon is still unclear. The non-coding T cells transcript (NTT) is differentially expressed during human T cells activation, but its function is unknown. Here, we detected a 426 m NTT transcript by RT-PCR using RNA of human lymphocytes activated with a synthetic peptide of HIV-1. After cloning, the sense and antisense 426 nt NTT transcripts were obtained by in vitro transcription and were sequenced. We found that both transcripts are highly structured and are able to activate PKR. A striking observation was that the antisense 426 nt NTT transcript is significantly more effective in activating PKR than the corresponding sense transcript. The transcription factor NF-kappa B is activated by PKR through phosphorylation and subsequent degradation of its inhibitor I-kappa B beta. We also found that the antisense 426 nt NTT transcript induces more efficiently the degradation Of I-kappa B beta than the sense transcript. Thus, this study suggests that the role played by NTT in the activation of lymphocytes can be mediated by PKR through NF-kappa B activation. However, the physiological significance of the activity of the antisense 426 nt NTT transcript remains unknown. (c) 2007 Elsevier Inc. All rights reserved.

Non-coding RNAs: the architects of eukaryotic complexity

Around 98% of all transcriptional output in humans is noncoding RNA. RNA-mediated gene regulation is widespread in higher eukaryotes and complex genetic phenomena like RNA interference, co-suppression, transgene silencing, imprinting, methylation, and possibly position-effect variegation and transvection, all involve intersecting pathways based on or connected to RNA signaling. I suggest that the central dogma is incomplete, and that intronic and other non-coding RNAs have evolved to comprise a second tier of gene expression in eukaryotes, which enables the integration and networking of complex suites of gene activity. Although proteins are the fundamental effectors of cellular function, the basis of eukaryotic complexity and phenotypic variation may lie primarily in a control architecture composed of a highly parallel system of trans-acting RNAs that relay state information required for the coordination and modulation of gene expression, via chromatin remodeling, RNA-DNA, RNA-RNA and RNA-protein interactions. This system has interesting and perhaps informative analogies with small world networks and dataflow computing.

Non-coding RNAs: multi-tasking molecules in the cell

In the last years it has become increasingly clear that the mammalian transcriptome is highly complex and includes a large number of small non-coding RNAs (sncRNAs) and long noncoding RNAs (lncRNAs). Here we review the biogenesis pathways of the three classes of sncRNAs, namely short interfering RNAs (siRNAs), microRNAs (miRNAs) and PIWI-interacting RNAs (piRNAs). These ncRNAs have been extensively studied and are involved in pathways leading to specific gene silencing and the protection of genomes against virus and transposons, for example. Also, lncRNAs have emerged as pivotal molecules for the transcriptional and post-transcriptional regulation of gene expression which is supported by their tissue-specific expression patterns, subcellular distribution, and developmental regulation. Therefore, we also focus our attention on their role in differentiation and development. SncRNAs and lncRNAs play critical roles in defining DNA methylation patterns, as well as chromatin remodeling thus having a substantial effect in epigenetics. The identification of some overlaps in their biogenesis pathways and functional roles raises the hypothesis that these molecules play concerted functions in vivo, creating complex regulatory networks where cooperation with regulatory proteins is necessary. We also highlighted the implications of biogenesis and gene expression deregulation of sncRNAs and lncRNAs in human diseases like cancer.

Evolutionary dynamics of coding and non-coding transcriptomes.

Gene expression changes may underlie much of phenotypic evolution. The development of high-throughput RNA sequencing protocols has opened the door to unprecedented large-scale and cross-species transcriptome comparisons by allowing accurate and sensitive assessments of transcript sequences and expression levels. Here, we review the initial wave of the new generation of comparative transcriptomic studies in mammals and vertebrate outgroup species in the context of earlier work. Together with various large-scale genomic and epigenomic data, these studies have unveiled commonalities and differences in the dynamics of gene expression evolution for various types of coding and non-coding genes across mammalian lineages, organs, developmental stages, chromosomes and sexes. They have also provided intriguing new clues to the regulatory basis and phenotypic implications of evolutionary gene expression changes.

Chromatoid body mediated RNA regulation in mouse male germline

Male germ cell differentiation, spermatogenesis is an exceptional developmental process that produces a massive amount of genetically unique spermatozoa. The complexity of this process along with the technical limitations in the germline research has left many aspects of spermatogenesis poorly understood. Post-meiotic haploid round spermatids possess the most complex transcriptomes of the whole body. Correspondingly, efficient and accurate control mechanisms are necessary to deal with the huge diversity of transcribed RNAs in these cells. The high transcriptional activity in round spermatids is accompanied by the presence of an uncommonly large cytoplasmic ribonucleoprotein granule, called the chromatoid body (CB) that is conjectured to participate in the RNA post-transcriptional regulation. However, very little is known about the possible mechanisms of the CB function. The development of a procedure to isolate CBs from mouse testes was this study’s objective. Anti-MVH immunoprecipitation of cross-linked CBs from a fractionated testicular cell lysate was optimized to yield considerable quantities of pure and intact CBs from mice testes. This protocol produced reliable and reproducible data from the subsequent analysis of CB’s protein and RNA components. We found that the majority of the CB’s proteome consists of RNA-binding proteins that associate functionally with different pathways. We also demonstrated notable localization patterns of one of the CB transient components, SAM68 and showed that its ablation does not change the general composition or structure of the CB. CB-associated RNA analysis revealed a strong accumulation of PIWI-interacting RNAs (piRNAs), mRNAs and long non-coding RNAs (lncRNAs) in the CB. When the CB transcriptome and proteome analysis results were combined, the most pronounced molecular functions in the CB were related to piRNA pathway, RNA post-transcriptional processing and CB structural scaffolding. In addition, we demonstrated that the CB is a target for the main RNA flux from the nucleus throughout all steps of round spermatid development. Moreover, we provided preliminary evidence that those isolated CBs slice target RNAs in vitro in an ATPdependent manner. Altogether, these results make a strong suggestion that the CB functions involve RNA-related and RNA-mediated mechanisms. All the existing data supports the hypothesis that the CB coordinates the highly complex haploid transcriptome during the preparation of the male gametes for fertilization. Thereby, this study provides a fundamental basis for the future functional analyses of ribonucleoprotein granules and offers also important insights into the mechanisms governing male fertility.

Expression of Mitochondrial Non-coding RNAs (ncRNAs) Is Modulated by High Risk Human Papillomavirus (HPV) Oncogenes

The study of RNA and DNA oncogenic viruses has proved invaluable in the discovery of key cellular pathways that are rendered dysfunctional during cancer progression. An example is high risk human papillomavirus (HPV), the etiological agent of cervical cancer. The role of HPV oncogenes in cellular immortalization and transformation has been extensively investigated. We reported the differential expression of a family of human mitochondrial non-coding RNAs (ncRNAs) between normal and cancer cells. Normal cells express a sense mitochondrial ncRNA (SncmtRNA) that seems to be required for cell proliferation and two antisense transcripts (ASncmtRNAs). In contrast, the ASncmtRNAs are down-regulated in cancer cells. To shed some light on the mechanisms that trigger down-regulation of the ASncmtRNAs, we studied human keratinocytes (HFK) immortalized with HPV. Here we show that immortalization of HFK with HPV-16 or 18 causes down-regulation of the ASncmtRNAs and induces the expression of a new sense transcript named SncmtRNA-2. Transduction of HFK with both E6 and E7 is sufficient to induce expression of SncmtRNA-2. Moreover, E2 oncogene is involved in down-regulation of the ASncmtRNAs. Knockdown of E2 in immortalized cells reestablishes in a reversible manner the expression of the ASncmtRNAs, suggesting that endogenous cellular factors(s) could play functions analogous to E2 during non-HPV-induced oncogenesis.

Trypanosoma cruzi mitochondrial maxicircles display species- and strain-specific variation and a conserved element in the non-coding region

Abstract Background The mitochondrial DNA of kinetoplastid flagellates is distinctive in the eukaryotic world due to its massive size, complex form and large sequence content. Comprised of catenated maxicircles that contain rRNA and protein-coding genes and thousands of heterogeneous minicircles encoding small guide RNAs, the kinetoplast network has evolved along with an extreme form of mRNA processing in the form of uridine insertion and deletion RNA editing. Many maxicircle-encoded mRNAs cannot be translated without this post-transcriptional sequence modification. Results We present the complete sequence and annotation of the Trypanosoma cruzi maxicircles for the CL Brener and Esmeraldo strains. Gene order is syntenic with Trypanosoma brucei and Leishmania tarentolae maxicircles. The non-coding components have strain-specific repetitive regions and a variable region that is unique for each strain with the exception of a conserved sequence element that may serve as an origin of replication, but shows no sequence identity with L. tarentolae or T. brucei. Alternative assemblies of the variable region demonstrate intra-strain heterogeneity of the maxicircle population. The extent of mRNA editing required for particular genes approximates that seen in T. brucei. Extensively edited genes were more divergent among the genera than non-edited and rRNA genes. Esmeraldo contains a unique 236-bp deletion that removes the 5'-ends of ND4 and CR4 and the intergenic region. Esmeraldo shows additional insertions and deletions outside of areas edited in other species in ND5, MURF1, and MURF2, while CL Brener has a distinct insertion in MURF2. Conclusion The CL Brener and Esmeraldo maxicircles represent two of three previously defined maxicircle clades and promise utility as taxonomic markers. Restoration of the disrupted reading frames might be accomplished by strain-specific RNA editing. Elements in the non-coding region may be important for replication, transcription, and anchoring of the maxicircle within the kinetoplast network.

Androgen responsive intronic non-coding RNAs

Abstract Background Transcription of large numbers of non-coding RNAs originating from intronic regions of human genes has been recently reported, but mechanisms governing their biosynthesis and biological functions are largely unknown. In this work, we evaluated the existence of a common mechanism of transcription regulation shared by protein-coding mRNAs and intronic RNAs by measuring the effect of androgen on the transcriptional profile of a prostate cancer cell line. Results Using a custom-built cDNA microarray enriched in intronic transcribed sequences, we found 39 intronic non-coding RNAs for which levels were significantly regulated by androgen exposure. Orientation-specific reverse transcription-PCR indicated that 10 of the 13 were transcribed in the antisense direction. These transcripts are long (0.5–5 kb), unspliced and apparently do not code for proteins. Interestingly, we found that the relative levels of androgen-regulated intronic transcripts could be correlated with the levels of the corresponding protein-coding gene (asGAS6 and asDNAJC3) or with the alternative usage of exons (asKDELR2 and asITGA6) in the corresponding protein-coding transcripts. Binding of the androgen receptor to a putative regulatory region upstream from asMYO5A, an androgen-regulated antisense intronic transcript, was confirmed by chromatin immunoprecipitation. Conclusion Altogether, these results indicate that at least a fraction of naturally transcribed intronic non-coding RNAs may be regulated by common physiological signals such as hormones, and further corroborate the notion that the intronic complement of the transcriptome play functional roles in the human gene-expression program.