Featuring BRCA1 and BRCA2 germline mutational landscape from Asturias (North Spain).

The singular BRCA1/2 mutational landscape of Asturias is updated 10 years after the first study. We analyzed BRCA1 and BRCA2 pathogenic variants in 1653 index cases. In total, 238 families were identified to carry a pathogenic variant, 163 families in BRCA1 and 75 families in BRCA2. This yielded a prevalence rate of 14.4%. Seven recurrent variants were found accounting for 55% of the cases. Among them, three are widely distributed (BRCA1 c.211A>G, c.470_471del and c.3331_3334del) and four had been reported as novel in Asturias: two in BRCA1 (c.1674del and c.2901_2902dup) and two in BRCA2 (c.2095C>T and c.4040_4035delinsC). A common haplotype was established for all recurrent variants indicating a shared ancestral origin. Three splicing analyses are shown: BRCA1:c.5152+3A>C and BRCA1:c.5333-3T>G that lead to skipping of exon 18, and 22 respectively, and BRCA1:c.5278-1G>T giving rise to two transcripts, one lacking exon 21 (p.Ille1760Glyfs*60) and one lacking the first 8 nucleotides of exon 21 (p.Phe1761Asnfs*14), supporting pathogenicity for these variants.

Target-capture full-length double-stranded cDNA long-read sequencing through Nanopore revealed novel intron retention in patient with tuberous sclerosis complex.

Tuberous sclerosis complex (TSC) is a relatively common autosomal dominant disorder characterized by multiple dysplastic organ lesions and neuropsychiatric symptoms caused by loss-of-function mutation of either TSC1 or TSC2. The genetic diagnosis of inherited diseases, including TSC, in the clinical field is widespread using next-generation sequencing. The mutations in protein-coding exon tend to be verified because mutations directly cause abnormal protein. However, it is relatively difficult to verify mutations in the intron region because it is required to investigate whether the intron mutations affect the abnormal splicing of transcripts. In this study, we developed a target-capture full-length double-stranded cDNA sequencing method using Nanopore long-read sequencer (Nanopore long-read target sequencing). This method revealed the occurrence of intron mutation in the TSC2 gene and found that the intron mutation produces novel intron retention splicing transcripts that generate truncated proteins. The protein-coding transcripts were decreased due to the expression of the novel intron retention transcripts, which caused TSC in patients with the intron mutation. Our results indicate that Nanopore long-read target sequencing is useful for the detection of mutations and confers information on the full-length alternative splicing of transcripts for genetic diagnosis.

One pipeline to predict them all? On the prediction of alternative splicing from RNA-Seq data.

RNA-Seq has become the standard approach to quantify and compare gene expression and alternative splicing in different conditions. In many cases the limiting factor is not the sequencing itself but the bioinformatic analysis. A variety of software tools exist that predict alternative splicing patterns from RNA-Seq data, but surprisingly, a systematic comparison of the predictions obtained from different pipelines has not been performed. Here we compare results from frequently used bioinformatic tools using a high-quality RNA-Seq dataset. We show that there is little overlap in the splicing changes predicted by different tools and that GO-term analysis of the splicing changes predicted by the individual targets yields very different results. Validation of bioinformatic predictions by RT-PCR suggest a high number of false positives in the splicing changes predicated by each pipeline, which probably dominates GO-term analysis. The validation rate is strongly increased for targets predicted by several tools, offering a strategy to reduce false positives. Based on these results we offer some guidelines that may contribute to make alternative splicing predictions more reliable and may thus increase the impact of conclusions drawn from RNA-Seq studies. Furthermore, we created rmappet, a nextflow pipeline that performs alternative splicing analysis using rMATS and Whippet with subsequent overlapping of the results, enabling robust splicing analysis with only one command (https://github.com/didrikolofsson/rmappet/).

A Patient with Corticobasal Syndrome and Progressive Non-Fluent Aphasia (CBS-PNFA), with Variants in ATP7B, SETX, SORL1, and FOXP1 Genes.

Our aim was to analyze the phenotypic-genetic correlations in a patient diagnosed with early onset corticobasal syndrome with progressive non-fluent aphasia (CBS-PNFA), characterized by predominant apraxia of speech, accompanied by prominent right-sided upper-limb limb-kinetic apraxia, alien limb phenomenon, synkinesis, myoclonus, mild cortical sensory loss, and right-sided hemispatial neglect. Whole-exome sequencing (WES) identified rare single heterozygous variants in ATP7B (c.3207C>A), SORL1 (c.352G>A), SETX (c.2385_2387delAAA), and FOXP1 (c.1762G>A) genes. The functional analysis revealed that the deletion in the SETX gene changed the splicing pattern, which was accompanied by lower SETX mRNA levels in the patient's fibroblasts, suggesting loss-of-function as the underlying mechanism. In addition, the patient's fibroblasts demonstrated altered mitochondrial architecture with decreased connectivity, compared to the control individuals. This is the first association of the CBS-PNFA phenotype with the most common ATP7B pathogenic variant p.H1069Q, previously linked to Wilson's disease, and early onset Parkinson's disease. This study expands the complex clinical spectrum related to variants in well-known disease genes, such as ATP7B, SORL1, SETX, and FOXP1, corroborating the hypothesis of oligogenic inheritance. To date, the FOXP1 gene has been linked exclusively to neurodevelopmental speech disorders, while our study highlights its possible relevance for adult-onset progressive apraxia of speech, which guarantees further study.

DJExpress: An Integrated Application for Differential Splicing Analysis and Visualization.

RNA-seq analysis of alternative pre-mRNA splicing has facilitated an unprecedented understanding of transcriptome complexity in health and disease. However, despite the availability of countless bioinformatic pipelines for transcriptome-wide splicing analysis, the use of these tools is often limited to expert bioinformaticians. The need for high computational power, combined with computational outputs that are complicated to visualize and interpret present obstacles to the broader research community. Here we introduce DJExpress, an R package for differential expression analysis of transcriptomic features and expression-trait associations. To determine gene-level differential junction usage as well as associations between junction expression and molecular/clinical features, DJExpress uses raw splice junction counts as input data. Importantly, DJExpress runs on an average laptop computer and provides a set of interactive and intuitive visualization formats. In contrast to most existing pipelines, DJExpress can handle both annotated and de novo identified splice junctions, thereby allowing the quantification of novel splice events. Moreover, DJExpress offers a web-compatible graphical interface allowing the analysis of user-provided data as well as the visualization of splice events within our custom database of differential junction expression in cancer (DJEC DB). DJEC DB includes not only healthy and tumor tissue junction expression data from TCGA and GTEx repositories but also cancer cell line data from the DepMap project. The integration of DepMap functional genomics data sets allows association of junction expression with molecular features such as gene dependencies and drug response profiles. This facilitates identification of cancer cell models for specific splicing alterations that can then be used for functional characterization in the lab. Thus, DJExpress represents a powerful and user-friendly tool for exploration of alternative splicing alterations in RNA-seq data, including multi-level data integration of alternative splicing signatures in healthy tissue, tumors and cancer cell lines.

Exploring New Functional Aspects of HTLV-1 RNA-Binding Protein Rex: How Does Rex Control Viral Replication?

After integration to the human genome as a provirus, human T-cell leukemia virus type 1 (HTLV-1) utilizes host T cell gene expression machinery for viral replication. The viral RNA-binding protein, Rex, is known to transport unspliced/incompletely spliced viral mRNAs encoding viral structural proteins out of the nucleus to enhance virus particle formation. However, the detailed mechanism of how Rex avoids extra splicing of unspliced/incompletely spliced viral mRNAs and stabilizes them for effective translation is still unclear. To elucidate the underlying molecular mechanism of Rex function, we comprehensively analyzed the changes in gene expression and splicing patterns in Rex-overexpressing T cells. In addition, we identified 81 human proteins interacting with Rex, involved in transcription, splicing, translation, and mRNA quality control. In particular, Rex interacts with NONO and SFPQ, which play important roles in the regulation of transcription and splicing. Accordingly, expression profiles and splicing patterns of a wide variety of genes are significantly changed in Rex-expressing T cells. Especially, the level of vPD-L1 mRNA that lacks the part of exon 4, thus encodes soluble PD-L1 was significantly increased in Rex-expressing cells. Overall, by integrated analysis of these three datasets, we showed for the first time that Rex intervenes the host gene expression machinery throughout the pathway, probably to escort viral unstable mRNAs from transcription (start) to translation (end). Upon exerting its function, Rex may alter the expression level and splicing patterns of various genes, thus influencing the phenotype of the host cell.

The First Korean Case of SLC12A3 Aberrant Skipping of Two Exons Detected by RNA Splicing Analysis.

Gitelman syndrome is a salt-losing tubular disorder that is transmitted as an autosomal recessive trait. Variants in the SLC12A3 gene are found in the majority of Gitelman syndrome patients. A 26-year-old woman visited the genetic counseling clinic. Her fiancé was a known Gitelman syndrome patient who was previously diagnosed with 2 pathogenic variants in SLC12A3. In advance of marriage and future family planning, she wanted to perform genetic testing of SLC12A3. A silent exonic variant c.1050G>A was found, and multiple splice site in silico algorithms predicted this variant to have potential alteration of splicing. This variant was classified as "variant of uncertain significance," and RNA splicing analysis was additionally performed. RNA splicing analysis showed aberrant splicing of exon 7-8 skipping. The result points out the potential pathogenicity of this variant, which should be considered a candidate of variant reclassification in the future. We highly recommend the performance of additional RNA splicing analysis, especially for silent variants predicted to have potential alteration of splicing.

Usefulness of functional splicing analysis to confirm precise disease pathogenesis in Diamond-Blackfan anemia caused by intronic variants in RPS19.

Diamond-Blackfan anemia (DBA) is mainly caused by pathogenic variants in ribosomal proteins and 22 responsible genes have been identified to date. The most common causative gene of DBA is RPS19 [NM_001022.4]. Nearly 180 RPS19 variants have been reported, including three deep intronic variants outside the splicing consensus sequence (c.72-92A > G, c.356 + 18G > C, and c.411 + 6G > C). We also identified one case with a c.412-3C > G intronic variant. Without conducting transcript analysis, the pathogenicity of these variants is unknown. However, it is difficult to assess transcripts because of their fragility. In such cases, in vitro functional splicing assays can be used to assess pathogenicity. Here, we report functional splicing analysis results of four RPS19 deep intronic variants identified in our case and in previously reported cases. One splicing consensus variant (c.411 + 1G > A) was also examined as a positive control. Aberrant splicing with a 2-bp insertion between exons 5 and 6 was identified in the patient samples and minigene assay results also identified exon 6 skipping in our case. The exon 6 skipping transcript was confirmed by further evaluation using quantitative RT-PCR. Additionally, minigene assay analysis of three reported deep intronic variants revealed that none of them showed aberrant splicing and that these variants were not considered to be pathogenic. In conclusion, the minigene assay is a useful method for functional splicing analysis of inherited disease.

OneStopRNAseq: A Web Application for Comprehensive and Efficient Analyses of RNA-Seq Data.

Over the past decade, a large amount of RNA sequencing (RNA-seq) data were deposited in public repositories, and more are being produced at an unprecedented rate. However, there are few open source tools with point-and-click interfaces that are versatile and offer streamlined comprehensive analysis of RNA-seq datasets. To maximize the capitalization of these vast public resources and facilitate the analysis of RNA-seq data by biologists, we developed a web application called OneStopRNAseq for the one-stop analysis of RNA-seq data. OneStopRNAseq has user-friendly interfaces and offers workflows for common types of RNA-seq data analyses, such as comprehensive data-quality control, differential analysis of gene expression, exon usage, alternative splicing, transposable element expression, allele-specific gene expression quantification, and gene set enrichment analysis. Users only need to select the desired analyses and genome build, and provide a Gene Expression Omnibus (GEO) accession number or Dropbox links to sequence files, alignment files, gene-expression-count tables, or rank files with the corresponding metadata. Our pipeline facilitates the comprehensive and efficient analysis of private and public RNA-seq data.

Rapid exome sequencing and adjunct RNA studies confirm the pathogenicity of a novel homozygous ASNS splicing variant in a critically ill neonate.

Rapid genomic diagnosis programs are transforming rare disease diagnosis in acute pediatrics. A ventilated newborn with cerebellar hypoplasia underwent rapid exome sequencing (75 h), identifying a novel homozygous ASNS splice-site variant (NM_133436.3:c.1476+1G>A) of uncertain significance. Rapid ASNS splicing studies using blood-derived messenger RNA from the family trio confirmed a consistent pattern of abnormal splicing induced by the variant (cryptic 5' splice-site or exon 12 skipping) with absence of normal ASNS splicing in the proband. Splicing studies reported within 10 days led to reclassification of c.1476+1G>A as pathogenic at age 27 days. Intensive care was redirected toward palliation. Cost analyses for the neonate and his undiagnosed, similarly affected deceased sibling, demonstrate that early diagnosis reduced hospitalization costs by AU$100,828. We highlight the diagnostic benefits of adjunct RNA testing to confirm the pathogenicity of splicing variants identified via rapid genomic testing pipelines for precision and preventative medicine.

A 3D Cell Culture Model Identifies Wnt/ß-Catenin Mediated Inhibition of p53 as a Critical Step during Human Hepatocyte Regeneration.

The liver is a highly regenerative organ. While mature hepatocytes under homeostatic conditions are largely quiescent, upon injury, they rapidly enter the cell cycle to recover the damaged tissue. In rodents, a variety of injury models have provided important insights into the molecular underpinnings that govern the proliferative activation of quiescent hepatocytes. However, little is known about the molecular mechanisms of human hepatocyte regeneration and experimental methods to expand primary human hepatocytes (PHH). Here, a 3D spheroid model of PHH is established to study hepatocyte regeneration and integrative time-lapse multi-omics analyses show that upon isolation from the native liver PHH acquire a regenerative phenotype, as seen in vivo upon partial hepatectomy. However, proliferation is limited. By analyzing global promoter motif activities, it is predicted that activation of Wnt/ß-catenin and inhibition of p53 signaling are critical factors required for human hepatocyte proliferation. Functional validations reveal that activation of Wnt signaling through external cues alone is sufficient to inhibit p53 and its proliferative senescence-inducing target PAI1 (SERPINE1) and drive proliferation of >50% of all PHH. A scalable 3D culture model is established to study the molecular and cellular biology of human hepatocyte regeneration. By using this model, an essential role of Wnt/ß-catenin signaling during human hepatocyte regeneration is identified.

Interactive Alternative Splicing Analysis of Human Stem Cells Using psichomics.

Alternative splicing (AS) generates functionally distinct transcripts and is involved in multiple cellular processes, including stem cell differentiation. Several epithelial-to-mesenchymal transition-related splicing factors have also been associated with pluripotency. Concomitantly with the interest in studying AS in stem cell biology, the advent of next-generation sequencing of RNA (RNA-seq) has increased the public availability of transcriptomic data and enabled genome-wide AS studies. To facilitate performing such analyses in large publicly available or user-provided transcriptomics datasets, the psichomics R package provides an easy-to-use interface and efficient data visualization tools for AS quantification and integrative analyses of AS and gene expression data.psichomics is employed herein to study AS changes between human stem cells and fibroblasts, based on dimensionality reduction, and median- and variance-based differential AS and gene expression analyses. Putative RNA-binding protein regulators involved in those alterations are then identified based on correlation analyses in large cohorts of human tissue transcriptomes. We identified several alterations, both novel and previously reported, in alternative splicing events and in the expression of their candidate regulators that are associated with stem cell differentiation into fibroblasts.

Comparative transcriptome analysis revealed genes involved in the fruiting body development of Ophiocordyceps sinensis.

Ophiocordyceps sinensis is a highly valued fungus that has been used as traditional Asian medicine. This fungus is one of the most important sources of income for the nomadic populations of the Tibetan Plateau. With global warming and excessive collection, the wild O. sinensis resources declined dramatically. The cultivation of O. sinensis hasn't been fully operational due to the unclear genetic basis of the fruiting body development. Here, our study conducted pairwise comparisons between transcriptomes acquired from different growth stages of O. sinensis including asexual mycelium (CM), developing fruiting body (DF) and mature fruiting body (FB). All RNA-Seq reads were aligned to the genome of O. sinensis CO18 prior to comparative analyses. Cluster analysis showed that the expression profiles of FB and DF were highly similar compared to CM. Alternative splicing analysis (AS) revealed that the stage-specific splicing genes may have important functions in the development of fruiting body. Functional enrichment analyses showed that differentially expressed genes (DEGs) were enriched in protein synthesis and baseline metabolism during fruiting body development, indicating that more protein and energy might be required for fruiting body development. In addition, some fruiting body development-associated genes impacted by ecological factors were up-regulated in FB samples, such as the nucleoside diphosphate kinase gene (ndk), ß subunit of the fatty acid synthase gene (cel-2) and the superoxide dismutase gene (sod). Moreover, the expression levels of several cytoskeletons genes were significantly altered during all these growth stages, suggesting that these genes play crucial roles in both vegetative growth and the fruiting body development. Quantitative PCR (qPCR) was used to validate the gene expression profile and the results supported the accuracy of the RNA-Seq and DEGs analysis. Our study offers a novel perspective to understand the underlying growth stage-specific molecular differences and the biology of O. sinensis fruiting body development.

Splicing dysregulation contributes to the pathogenicity of several F9 exonic point variants.

BACKGROUND: Pre-mRNA splicing is a complex process requiring the identification of donor site, acceptor site, and branch point site with an adjacent polypyrimidine tract sequence. Splicing is regulated by splicing regulatory elements (SREs) with both enhancer and suppressor functions. Variants located in exonic regions can impact splicing through dysregulation of native splice sites, SREs, and cryptic splice site activation. While splicing dysregulation is considered primary disease-inducing mechanism of synonymous variants, its contribution toward disease phenotype of non-synonymous variants is underappreciated. METHODS: In this study, we analyzed 415 disease-causing and 120 neutral F9 exonic point variants including both synonymous and non-synonymous for their effect on splicing using a series of in silico splice site prediction tools, SRE prediction tools, and in vitro minigene assays. RESULTS: The use of splice site and SRE prediction tools in tandem provided better prediction but were not always in agreement with the minigene assays. The net effect of splicing dysregulation caused by variants was context dependent. Minigene assays revealed that perturbed splicing can be found. CONCLUSION: Synonymous variants primarily cause disease phenotype via splicing dysregulation while additional mechanisms such as translation rate also play an important role. Splicing dysregulation is likely to contribute to the disease phenotype of several non-synonymous variants.

Functional Studies and In Silico Analyses to Evaluate Non-Coding Variants in Inherited Cardiomyopathies.

Point mutations are the most common cause of inherited diseases. Bioinformatics tools can help to predict the pathogenicity of mutations found during genetic screening, but they may work less well in determining the effect of point mutations in non-coding regions. In silico analysis of intronic variants can reveal their impact on the splicing process, but the consequence of a given substitution is generally not predictable. The aim of this study was to functionally test five intronic variants (MYBPC3-c.506-2A>C, MYBPC3-c.906-7G>T, MYBPC3-c.2308+3G>C, SCN5A-c.393-5C>A, and ACTC1-c.617-7T>C) found in five patients affected by inherited cardiomyopathies in the attempt to verify their pathogenic role. Analysis of the MYBPC3-c.506-2A>C mutation in mRNA from the peripheral blood of one of the patients affected by hypertrophic cardiac myopathy revealed the loss of the canonical splice site and the use of an alternative splicing site, which caused the loss of the first seven nucleotides of exon 5 (MYBPC3-G169AfsX14). In the other four patients, we generated minigene constructs and transfected them in HEK-293 cells. This minigene approach showed that MYBPC3-c.2308+3G>C and SCN5A-c.393-5C>A altered pre-mRNA processing, thus resulting in the skipping of one exon. No alterations were found in either MYBPC3-c.906-7G>T or ACTC1-c.617-7T>C. In conclusion, functional in vitro analysis of the effects of potential splicing mutations can confirm or otherwise the putative pathogenicity of non-coding mutations, and thus help to guide the patient's clinical management and improve genetic counseling in affected families.

Novel CDH1 germline mutations identified in Chinese gastric cancer patients.

AIM: To give a comprehensive report of E-cadherin gene (CDH1) variations in a population at a high risk for gastric cancer (GC). METHODS: The samples consisted of 178 men and 58 women with a mean age of 62.3 ± 9.4 years and an age range of 30-84 years. A total of 240 cancer-free controls were recruited (mean age of 61.8 ± 10.1 years, age range of 26-82 years). Samples were screened for CDH1 germline mutations by high-resolution melting analysis or directly sequencing. Luciferase reporter assay, RNA splicing assay and bioinformatic analysis were used to evaluate the effect of mutations. RESULTS: Four novel CDH1 sequence alterations were identified in GC patients including a G>T transition 49 bp before the start codon; a three-nucleotide deletion, c.44_46del TGC; one missense mutation, c.604G>A (V202I); and one variation in the intron, c.1320+7A>G. In addition, polymorphism frequencies were observed for CDH1-164delT, -161C>A, -73A>C, c.48+6C>T, c.48+62_48+63delinsCGTGCCCCAGCCC, c.894C>T (A298A), c.1224G>A (A408A), c.1888C>G (L630V), c.2076T>C (A692A), and c.2253C>T (N751N) which is similar to the data reported in http://www.ncbi.nlm.nih.gov/projects/SNP/. RNA splicing analysis suggested that the c.1320+7A>G and c.1224G>A variations did not affect exon splicing ability. Luciferase reporter assay demonstrated that the c.-49T variation might be helpful for E-cadherin transcription, though the increase in transcription activity is limited (only 33%). SIFT score and PolyPhen analysis both demonstrated that the L630V missense mutation probably damages protein function, while the V202I variant does not. CONCLUSION: This study reveals novel mutations in sporadic GC patients which had been poorly investigated for susceptibility genes.