Genome Alert!: A standardized procedure for genomic variant reinterpretation and automated gene-phenotype reassessment in clinical routine.

PURPOSE: Retrospective interpretation of sequenced data in light of the current literature is a major concern of the field. Such reinterpretation is manual and both human resources and variable operating procedures are the main bottlenecks. METHODS: Genome Alert! method automatically reports changes with potential clinical significance in variant classification between releases of the ClinVar database. Using ClinVar submissions across time, this method assigns validity category to gene-disease associations. RESULTS: Between July 2017 and December 2019, the retrospective analysis of ClinVar submissions revealed a monthly median of 1247 changes in variant classification with potential clinical significance and 23 new gene-disease associations. Re-examination of 4929 targeted sequencing files highlighted 45 changes in variant classification, and of these classifications, 89% were expert validated, leading to 4 additional diagnoses. Genome Alert! gene-disease association catalog provided 75 high-confidence associations not available in the OMIM morbid list; of which, 20% became available in OMIM morbid list For more than 356 negative exome sequencing data that were reannotated for variants in these 75 genes, this elective approach led to a new diagnosis. CONCLUSION: Genome Alert! (https://genomealert.univ-grenoble-alpes.fr/) enables systematic and reproducible reinterpretation of acquired sequencing data in a clinical routine with limited human resource effect.

Kmerator Suite: design of specific k-mer signatures and automatic metadata discovery in large RNA-seq datasets.

The huge body of publicly available RNA-sequencing (RNA-seq) libraries is a treasure of functional information allowing to quantify the expression of known or novel transcripts in tissues. However, transcript quantification commonly relies on alignment methods requiring a lot of computational resources and processing time, which does not scale easily to large datasets. K-mer decomposition constitutes a new way to process RNA-seq data for the identification of transcriptional signatures, as k-mers can be used to quantify accurately gene expression in a less resource-consuming way. We present the Kmerator Suite, a set of three tools designed to extract specific k-mer signatures, quantify these k-mers into RNA-seq datasets and quickly visualize large dataset characteristics. The core tool, Kmerator, produces specific k-mers for 97% of human genes, enabling the measure of gene expression with high accuracy in simulated datasets. KmerExploR, a direct application of Kmerator, uses a set of predictor gene-specific k-mers to infer metadata including library protocol, sample features or contaminations from RNA-seq datasets. KmerExploR results are visualized through a user-friendly interface. Moreover, we demonstrate that the Kmerator Suite can be used for advanced queries targeting known or new biomarkers such as mutations, gene fusions or long non-coding RNAs for human health applications.

First Identification of RNA-Binding Proteins That Regulate Alternative Exons in the Dystrophin Gene.

The Duchenne muscular dystrophy (DMD) gene has a complex expression pattern regulated by multiple tissue-specific promoters and by alternative splicing (AS) of the resulting transcripts. Here, we used an RNAi-based approach coupled with DMD-targeted RNA-seq to identify RNA-binding proteins (RBPs) that regulate splicing of its skeletal muscle isoform (Dp427m) in a human muscular cell line. A total of 16 RBPs comprising the major regulators of muscle-specific splicing events were tested. We show that distinct combinations of RBPs maintain the correct inclusion in the Dp427m of exons that undergo spatio-temporal AS in other dystrophin isoforms. In particular, our findings revealed the complex networks of RBPs contributing to the splicing of the two short DMD exons 71 and 78, the inclusion of exon 78 in the adult Dp427m isoform being crucial for muscle function. Among the RBPs tested, QKI and DDX5/DDX17 proteins are important determinants of DMD exon inclusion. This is the first large-scale study to determine which RBP proteins act on the physiological splicing of the DMD gene. Our data shed light on molecular mechanisms contributing to the expression of the different dystrophin isoforms, which could be influenced by a change in the function or expression level of the identified RBPs.

RNA-Seq Analysis to Detect Abnormal Fusion Transcripts Linked to Chromothripsis.

RNA-Seq approach enables the detection and characterization of fusion or chimeric transcript associated to complex genome rearrangement. Until now, these events are classically identified at DNA level.Here we describe a complete procedure including a novel way of analyzing reads that combines genomic locations and local coverage to directly infer chimeric junctions with a high sensitivity and specificity, allowing identification of different classes of chimeric RNA events. We also recommend the best practices for the bioinformatics analysis and describe the experimental process for RNA validation using real-time PCR and sequencing.

Targeted RNA-Seq profiling of splicing pattern in the DMD gene: exons are mostly constitutively spliced in human skeletal muscle.

We have analysed the splicing pattern of the human Duchenne Muscular Dystrophy (DMD) transcript in normal skeletal muscle. To achieve depth of coverage required for the analysis of this lowly expressed gene in muscle, we designed a targeted RNA-Seq procedure that combines amplification of the full-length 11.3 kb DMD cDNA sequence and 454 sequencing technology. A high and uniform coverage of the cDNA sequence was obtained that allowed to draw up a reliable inventory of the physiological alternative splicing events in the muscular DMD transcript. In contrast to previous assumptions, we evidenced that most of the 79 DMD exons are constitutively spliced in skeletal muscle. Only a limited number of 12 alternative splicing events were identified, all present at a very low level. These include previously known exon skipping events but also newly described pseudoexon inclusions and alternative 3' splice sites, of which one is the first functional NAGNAG splice site reported in the DMD gene. This study provides the first RNA-Seq-based reference of DMD splicing pattern in skeletal muscle and reports on an experimental procedure well suited to detect condition-specific differences in this low abundance transcript that may prove useful for diagnostic, research or RNA-based therapeutic applications.

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.

Tau excess impairs mitosis and kinesin-5 function, leading to aneuploidy and cell death.

In neurodegenerative diseases such as Alzheimer's disease (AD), cell cycle defects and associated aneuploidy have been described. However, the importance of these defects in the physiopathology of AD and the underlying mechanistic processes are largely unknown, in particular with respect to the microtubule (MT)-binding protein Tau, which is found in excess in the brain and cerebrospinal fluid of affected individuals. Although it has long been known that Tau is phosphorylated during mitosis to generate a lower affinity for MTs, there is, to our knowledge, no indication that an excess of this protein could affect mitosis. Here, we studied the effect of an excess of human Tau (hTau) protein on cell mitosis in vivo. Using the Drosophila developing wing disc epithelium as a model, we show that an excess of hTau induces a mitotic arrest, with the presence of monopolar spindles. This mitotic defect leads to aneuploidy and apoptotic cell death. We studied the mechanism of action of hTau and found that the MT-binding domain of hTau is responsible for these defects. We also demonstrate that the effects of hTau occur via the inhibition of the function of the kinesin Klp61F, the Drosophila homologue of kinesin-5 (also called Eg5 or KIF11). We finally show that this deleterious effect of hTau is also found in other Drosophila cell types (neuroblasts) and tissues (the developing eye disc), as well as in human HeLa cells. By demonstrating that MT-bound Tau inhibits the Eg5 kinesin and cell mitosis, our work provides a new framework to consider the role of Tau in neurodegenerative diseases.

piRNAs and epigenetic conversion in Drosophila.

Transposable element (TE) activity is repressed in the Drosophila germline by Piwi-Interacting RNAs (piRNAs), a class of small non-coding RNAs. These piRNAs are produced by discrete genomic loci containing TE fragments. In a recent publication, we tested for the existence of a strict epigenetic induction of piRNA production capacity by a locus in the D. melanogaster genome. We used 2 lines carrying a transgenic 7-copy tandem cluster (P-lacZ-white) at the same genomic site. This cluster generates in both lines a local heterochromatic sector. One line (T-1) produces high levels of ovarian piRNAs homologous to the P-lacZ-white transgenes and shows a strong capacity to repress homologous sequences in trans, whereas the other line (BX2) is devoid of both of these capacities. The properties of these 2 lines are perfectly stable over generations. We have shown that the maternal transmission of a cytoplasm carrying piRNAs from the first line can confer to the inert transgenic locus of the second, a totally de novo capacity to produce high levels of piRNAs as well as the ability to induce homology-dependent silencing in trans. These new properties are stably inherited over generations (n>50). Furthermore, the converted locus has itself become able to convert an inert transgenic locus via cytoplasmic maternal inheritance. This results in a stable epigenetic conversion process, which can be performed recurrently--a phenomenon termed paramutation and discovered in Maize 60 y ago. Paramutation in Drosophila corresponds to the first stable paramutation in animals and provides a model system to investigate the epigenetically induced emergence of a piRNA-producing locus, a crucial step in epigenome shaping. In this Extra View, we discuss some additional functional aspects and the possible molecular mechanism of this piRNA-linked paramutation.

AutomiG, a biosensor to detect alterations in miRNA biogenesis and in small RNA silencing guided by perfect target complementarity.

Defects in miRNA biogenesis or activity are associated to development abnormalities and diseases. In Drosophila, miRNAs are predominantly loaded in Argonaute-1, which they guide for silencing of target RNAs. The miRNA pathway overlaps the RNAi pathway in this organism, as miRNAs may also associate with Argonaute-2, the mediator of RNAi. We set up a gene construct in which a single inducible promoter directs the expression of the GFP protein as well as two miRNAs perfectly matching the GFP sequences. We show that self-silencing of the resulting automiG gene requires Drosha, Pasha, Dicer-1, Dicer-2 and Argonaute-2 loaded with the anti-GFP miRNAs. In contrast, self-silencing of the automiG gene does not involve Argonaute-1. Thus, automiG reports in vivo for both miRNA biogenesis and Ago-2 mediated silencing, providing a powerful biosensor to identify situations where miRNA or siRNA pathways are impaired. As a proof of concept, we used automiG as a biosensor to screen a chemical library and identified 29 molecules that strongly inhibit miRNA silencing, out of which 5 also inhibit RNAi triggered by long double-stranded RNA. Finally, the automiG sensor is also self-silenced by the anti-GFP miRNAs in HeLa cells and might be easily used to identify factors involved in miRNA biogenesis and silencing guided by perfect target complementarity in mammals.

Profiles of piRNA abundances at emerging or established piRNA loci are determined by local DNA sequences.

Piwi-interacting RNAs (piRNAs) ensure transposable element silencing in Drosophila, thereby preserving genome integrity across generations. Primary piRNAs arise from the processing of long RNA transcripts produced in the germ line by a limited number of telomeric and pericentromeric loci. Primary piRNAs bound to the Argonaute protein Aubergine then drive the production of secondary piRNAs through the "ping-pong" amplification mechanism that involves an interplay with piRNAs bound to the Argonaute protein Argonaute-3. We recently discovered that clusters of P-element-derived transgenes produce piRNAs and mediate silencing of homologous target transgenes in the female germ line. We also demonstrated that some clusters are able to convert other homologous inactive transgene clusters into piRNA-producing loci, which then transmit their acquired silencing capacity over generations. This paramutation phenomenon is mediated by maternal inheritance of piRNAs homologous to the transgenes. Here we further mined our piRNA sequencing data sets generated from various strains carrying transgenes with partial sequence homology at distinct genomic sites. This analysis revealed that same sequences in different genomic contexts generate highly similar profiles of piRNA abundances. The strong tendency of piRNAs for bearing a U at their 5' end has long been recognized. Our observations support the notion that, in addition, the relative frequencies of Drosophila piRNAs are locally determined by the DNA sequence of piRNA loci.

Paramutation in Drosophila linked to emergence of a piRNA-producing locus.

A paramutation is an epigenetic interaction between two alleles of a locus, through which one allele induces a heritable modification in the other allele without modifying the DNA sequence. The paramutated allele itself becomes paramutagenic, that is, capable of epigenetically converting a new paramutable allele. Here we describe a case of paramutation in animals showing long-term transmission over generations. We previously characterized a homology-dependent silencing mechanism referred to as the trans-silencing effect (TSE), involved in P-transposable-element repression in the germ line. We now show that clusters of P-element-derived transgenes that induce strong TSE can convert other homologous transgene clusters incapable of TSE into strong silencers, which transmit the acquired silencing capacity through 50 generations. The paramutation occurs without any need for chromosome pairing between the paramutagenic and the paramutated loci, and is mediated by maternal inheritance of cytoplasm carrying Piwi-interacting RNAs (piRNAs) homologous to the transgenes. The repression capacity of the paramutated locus is abolished by a loss-of-function mutation of the aubergine gene involved in piRNA biogenesis, but not by a loss-of-function mutation of the Dicer-2 gene involved in siRNA production. The paramutated cluster, previously producing barely detectable levels of piRNAs, is converted into a stable, strong piRNA-producing locus by the paramutation and becomes fully paramutagenic itself. Our work provides a genetic model for the emergence of piRNA loci, as well as for RNA-mediated trans-generational repression of transposable elements.

A luminal glycoprotein drives dose-dependent diameter expansion of the Drosophila melanogaster hindgut tube.

An important step in epithelial organ development is size maturation of the organ lumen to attain correct dimensions. Here we show that the regulated expression of Tenectin (Tnc) is critical to shape the Drosophila melanogaster hindgut tube. Tnc is a secreted protein that fills the embryonic hindgut lumen during tube diameter expansion. Inside the lumen, Tnc contributes to detectable O-Glycans and forms a dense striated matrix. Loss of tnc causes a narrow hindgut tube, while Tnc over-expression drives tube dilation in a dose-dependent manner. Cellular analyses show that luminal accumulation of Tnc causes an increase in inner and outer tube diameter, and cell flattening within the tube wall, similar to the effects of a hydrostatic pressure in other systems. When Tnc expression is induced only in cells at one side of the tube wall, Tnc fills the lumen and equally affects all cells at the lumen perimeter, arguing that Tnc acts non-cell-autonomously. Moreover, when Tnc expression is directed to a segment of a tube, its luminal accumulation is restricted to this segment and affects the surrounding cells to promote a corresponding local diameter expansion. These findings suggest that deposition of Tnc into the lumen might contribute to expansion of the lumen volume, and thereby to stretching of the tube wall. Consistent with such an idea, ectopic expression of Tnc in different developing epithelial tubes is sufficient to cause dilation, while epidermal Tnc expression has no effect on morphology. Together, the results show that epithelial tube diameter can be modelled by regulating the levels and pattern of expression of a single luminal glycoprotein.

Important neuronal toxicity of microtubule-bound Tau in vivo in Drosophila.

The microtubule-associated protein Tau is found in large amount in axons of neurons and is involved in human neurodegenerative diseases called tauopathies, which include Alzheimer's disease. In these diseases, the Tau protein is abnormally hyperphosphorylated and one therapeutic strategy currently under consideration consists in inhibiting Tau phosphorylation. However, the consequences of an excess of hypophosphorylated Tau onto neuronal physiology have not been investigated in vivo. Here we studied how important is Tau phosphorylation for axonal transport and neurohormone release in vivo, using the Drosophila model. Surprisingly, our results demonstrate a stronger toxicity of hypophosphorylated Tau for neuronal function, when compared with normal or pseudophosphorylated Tau. This reveals a potential limit of the current therapeutic strategy aimed at inhibiting Tau phosphorylation.

Tenectin is a novel alphaPS2betaPS integrin ligand required for wing morphogenesis and male genital looping in Drosophila.

Morphogenesis of the adult structures of holometabolous insects is regulated by ecdysteroids and juvenile hormones and involves cell-cell interactions mediated in part by the cell surface integrin receptors and their extracellular matrix (ECM) ligands. These adhesion molecules and their regulation by hormones are not well characterized. We describe the gene structure of a newly described ECM molecule, tenectin, and demonstrate that it is a hormonally regulated ECM protein required for proper morphogenesis of the adult wing and male genitalia. Tenectin's function as a new ligand of the PS2 integrins is demonstrated by both genetic interactions in the fly and by cell spreading and cell adhesion assays in cultured cells. Its interaction with the PS2 integrins is dependent on RGD and RGD-like motifs. Tenectin's function in looping morphogenesis in the development of the male genitalia led to experiments that demonstrate a role for PS integrins in the execution of left-right asymmetry.

The endogenous siRNA pathway is involved in heterochromatin formation in Drosophila.

A new class of small RNAs (endo-siRNAs) produced from endogenous double-stranded RNA (dsRNA) precursors was recently shown to mediate transposable element (TE) silencing in the Drosophila soma. These endo-siRNAs might play a role in heterochromatin formation, as has been shown in S. pombe for siRNAs derived from repetitive sequences in chromosome pericentromeres. To address this possibility, we used the viral suppressors of RNA silencing B2 and P19. These proteins normally counteract the RNAi host defense by blocking the biogenesis or activity of virus-derived siRNAs. We hypothesized that both proteins would similarly block endo-siRNA processing or function, thereby revealing the contribution of endo-siRNA to heterochromatin formation. Accordingly, P19 as well as a nuclear form of P19 expressed in Drosophila somatic cells were found to sequester TE-derived siRNAs whereas B2 predominantly bound their longer precursors. Strikingly, B2 or the nuclear form of P19, but not P19, suppressed silencing of heterochromatin gene markers in adult flies, and altered histone H3-K9 methylation as well as chromosomal distribution of histone methyl transferase Su(var)3-9 and Heterochromatin Protein 1 in larvae. Similar effects were observed in dcr2, r2d2, and ago2 mutants. Our findings provide evidence that a nuclear pool of TE-derived endo-siRNAs is involved in heterochromatin formation in somatic tissues in Drosophila.

Tenectin, a novel extracellular matrix protein expressed during Drosophila melanogaster embryonic development.

During Drosophila embryonic development, various morphogenetic processes require the remodeling of the extracellular matrix. In a previous study, we have identified and characterized a cDNA encoding a novel putative extracellular matrix protein named tenebrin, in the beetle Tenebrio molitor. Here, we examine the expression of the Drosophila ortholog, referred to as Tenectin (Tnc), during embryonic development. Tnc is expressed in the majority of tissues of neuroectodermic origin such as hindgut, foregut, tracheal system, anal plate, and CNS. In the CNS, the Tnc transcript is restricted to a few cells, whereas the protein is located in the dorsal part of the axonal tracts. In the hindgut and the trachea, Tnc protein is expressed on the apical pole of the cells. Tnc is an extracellular matrix protein secreted in a polarized way in different organs of Drosophila embryos.