RNA splicing analysis in genomic medicine

Spectrum and genotyping strategies of "dark" genetic matter in germline susceptibility genes of tumor syndromes

PURPOSE

Despite the widespread use of high-throughput genotyping strategies, certain mutation types remain understudied. We provide an overview of these often overlooked mutation types, with representative examples from common hereditary cancer syndromes.

METHODS

We conducted a comprehensive review of the literature and locus-specific variant databases to summarize the germline pathogenic variants discovered through non-routine genotyping methods. We evaluated appropriate detection and analysis methods tailored for these specific genetic aberrations. Additionally, we performed in silico splice predictions on deep intronic variants registered in the ClinVar database.

RESULTS

Our study suggests that, aside from founder mutations, most cases are sporadic. However, we anticipate a relatively high likelihood of splice effects for deep intronic variants. The findings underscore the significant clinical utility of genome sequencing techniques and the importance of applying relevant analysis methods.

Clinical and Genetic Analysis of 8 Children With Ornithine Transcarbamylase Deficiency: Two Novel Mutations

Background and Objectives

Cases and studies of neurologic symptoms in children caused by genetic metabolic diseases have been widely reported. Ornithine transcarbamylase deficiency (OTCD) is the most common urea cycle disorder, which is due to mutations in the OTC gene located on chromosome Xp21.1. In this study, we analyzed the clinical and genetic characteristics of 8 Chinese children diagnosed with OTCD.

Methods

A total of 8 patients (5 male and 3 female) were diagnosed with OTCD by biochemical and molecular analysis between 2015 and 2023. Clinical manifestations, biochemical features, and OTC gene sequencing analysis were reviewed retrospectively. The effect of c.664-1G>C on OTC mRNA synthesis was confirmed by a minigene splicing assay.

Results

All children were late-onset patients, with a median onset age of 3.6 years (range 1.8-17 years). Neurologic symptoms caused by hyperammonemia include vomiting, coma, dyssomnia, and seizures. The peak plasma ammonia levels ranged from 149 to 4,490 μmol/L, and alanine transaminase levels ranged from 20 to 1316 U/L. Four of them had received CRRT, and only 1 patient was admitted for liver transplantation. By December 2023, 4 patients had survived and 4 were deceased. Blood amino acids or urinary organic acids were detected in 7 cases. All patients underwent whole-exome sequencing, and 2 novel mutations were revealed (P1, c.617dupT and P2, c.664-1G>C). The alteration (c.664-1G>C) leads to the deletion of a 54-bp sequence in the exon 7 of the OTC gene (NM_000531.6: c.664_717del).

Discussion

Two novel pathogenic variants in the OTC gene were confirmed in 8 Chinese children by biochemical findings and genetic analysis. These findings will provide a better understanding of the diagnosis and treatment of pediatric patients with OTCD.

Disrupted Post-Transcriptional Regulation of Gene Expression as a Hallmark of Fatty Liver Progression

It is known that both transcriptional and post-transcriptional mechanisms control messenger RNA (mRNA) levels. Compared to transcriptional regulations, our understanding of how post-transcriptional regulations adapt during fatty liver progression at the whole-transcriptome level is unclear. While traditional RNA-seq analysis uses only reads mapped to exons to determine gene expression, recent studies support the idea that intron-mapped reads can be reliably used to estimate gene transcription. In this study, we analyzed differential gene expression at both the exon and intron levels using two liver RNA-seq datasets from mice that were fed a high-fat diet for seven weeks (mild fatty liver) or thirty weeks (severe fatty liver). We found that the correlation between gene transcription and mature mRNA levels was much lower in mice with mild fatty liver as compared with mice with severe fatty liver. This result indicates broad post-transcriptional regulations for early fatty liver and such regulations are compromised for severe fatty liver. Specifically, gene ontology analysis revealed that genes involved in synapse organization and cell adhesion were transcriptionally upregulated, while their mature mRNAs were unaffected in mild fatty liver. Further characterization of post-transcriptionally suppressed genes in early fatty liver revealed that their mRNAs harbor a significantly longer 3' UTR, one of the major features that may subject RNA transcripts to nonsense-mediated RNA decay (NMD). We further show that the expression of representative genes that were post-transcriptionally suppressed were upregulated in mice with a hepatocyte-specific defect of NMD. Finally, we provide data supporting a time-dependent decrease in NMD activity in the liver of a diet-induced metabolic-dysfunction-associated fatty liver disease mouse model. In summary, our study supports the conclusion that NMD is essential in preventing unwanted/harmful gene expression at the early stage of fatty liver and such a mechanism is lost due to decreased NMD activity in mice with severe fatty liver.

Tailored antisense oligonucleotides designed to correct aberrant splicing reveal actionable groups of mutations for rare genetic disorders

Effective translation of rare disease diagnosis knowledge into therapeutic applications is achievable within a reasonable timeframe; where mutations are amenable to current antisense oligonucleotide technology. In our study, we identified five distinct types of abnormal splice-causing mutations in patients with rare genetic disorders and developed a tailored antisense oligonucleotide for each mutation type using phosphorodiamidate morpholino oligomers with or without octa-guanidine dendrimers and 2'-O-methoxyethyl phosphorothioate. We observed variations in treatment effects and efficiencies, influenced by both the chosen chemistry and the specific nature of the aberrant splicing patterns targeted for correction. Our study demonstrated the successful correction of all five different types of aberrant splicing. Our findings reveal that effective correction of aberrant splicing can depend on altering the chemical composition of oligonucleotides and suggest a fast, efficient, and feasible approach for developing personalized therapeutic interventions for genetic disorders within short time frames.

Characterization of pre-mRNA Splicing Defects Caused by CLCN5 and OCRL Mutations and Identification of Novel Variants Associated with Dent Disease

Dent disease (DD) is an X-linked renal tubulopathy characterized by low-molecular-weight proteinuria, hypercalciuria, nephrocalcinosis, nephrolithiasis and progressive renal failure. Two-thirds of cases are associated with inactivating variants in the CLCN5 gene (Dent disease 1, DD1) and a few present variants in the OCRL gene (Dent disease 2, DD2). The aim of the present study was to test the effect on the pre-mRNA splicing process of DD variants, described here or in the literature, and describe the clinical and genotypic features of thirteen unrelated patients with suspected DD. All patients presented tubular proteinuria, ten presented hypercalciuria and five had nephrolithiasis or nephrocalcinosis. CLCN5 and OCRL genes were analyzed by Sanger sequencing. Nine patients showed variants in CLCN5 and four in OCRL; eight of these were new. Bioinformatics tools were used to select fifteen variants with a potential effect on pre-mRNA splicing from our patients' group and from the literature, and were experimentally tested using minigene assays. Results showed that three exonic missense mutations and two intronic variants affect the mRNA splicing process. Our findings widen the genotypic spectrum of DD and provide insight into the impact of variants causing DD.

Addition of an affected family member to a previously ascertained autosomal recessive nonsyndromic hearing loss pedigree and systematic phenotype-genotype analysis of splice-site variants in MYO15A

Pathogenic variants in MYO15A are known to cause autosomal recessive nonsyndromic hearing loss (ARNSHL), DFNB3. We have previously reported on one ARNSHL family including two affected siblings and identified MYO15A c.5964+3G > A and c.8375 T > C (p.Val2792Ala) as the possible deafness-causing variants. Eight year follow up identified one new affected individual in this family, who also showed congenital, severe to profound sensorineural hearing loss. By whole exome sequencing, we identified a new splice-site variant c.5531+1G > C (maternal allele), in a compound heterozygote with previously identified missense variant c.8375 T > C (p.Val2792Ala) (paternal allele) in MYO15A as the disease-causing variants. The new affected individual underwent unilateral cochlear implantation at the age of 1 year, and 5 year follow-up showed satisfactory speech and language outcomes. Our results further indicate that MYO15A-associated hearing loss is good candidates for cochlear implantation, which is in accordance with previous report. In light of our findings and review of the literatures, 58 splice-site variants in MYO15A are correlated with a severe deafness phenotype, composed of 46 canonical splice-site variants and 12 non-canonical splice-site variants.

Prevalence, parameters, and pathogenic mechanisms for splice-altering acceptor variants that disrupt the AG exclusion zone

Predicting the pathogenicity of acceptor splice-site variants outside the essential AG is challenging, due to high sequence diversity of the extended splice-site region. Critical analysis of 24,445 intronic extended acceptor splice-site variants reported in ClinVar and the Leiden Open Variation Database (LOVD) demonstrates 41.9% of pathogenic variants create an AG dinucleotide between the predicted branchpoint and acceptor (AG-creating variants in the AG exclusion zone), 28.4% result in loss of a pyrimidine at the -3 position, and 15.1% result in loss of one or more pyrimidines in the polypyrimidine tract. Pathogenicity of AG-creating variants was highly influenced by their position. We define a high-risk zone for pathogenicity: > 6 nucleotides downstream of the predicted branchpoint and >5 nucleotides upstream from the acceptor, where 93.1% of pathogenic AG-creating variants arise and where naturally occurring AG dinucleotides are concordantly depleted (5.8% of natural AGs). SpliceAI effectively predicts pathogenicity of AG-creating variants, achieving 95% sensitivity and 69% specificity. We highlight clinical examples showing contrasting mechanisms for mis-splicing arising from AG variants: (1) cryptic acceptor created; (2) splicing silencer created: an introduced AG silences the acceptor, resulting in exon skipping, intron retention, and/or use of an alternative existing cryptic acceptor; and (3) splicing silencer disrupted: loss of a deep intronic AG activates inclusion of a pseudo-exon. In conclusion, we establish AG-creating variants as a common class of pathogenic extended acceptor variant and outline factors conferring critical risk for mis-splicing for AG-creating variants in the AG exclusion zone, between the branchpoint and acceptor.

A systematic analysis of splicing variants identifies new diagnoses in the 100,000 Genomes Project

BACKGROUND

Genomic variants which disrupt splicing are a major cause of rare genetic diseases. However, variants which lie outside of the canonical splice sites are difficult to interpret clinically. Improving the clinical interpretation of non-canonical splicing variants offers a major opportunity to uplift diagnostic yields from whole genome sequencing data.

METHODS

Here, we examine the landscape of splicing variants in whole-genome sequencing data from 38,688 individuals in the 100,000 Genomes Project and assess the contribution of non-canonical splicing variants to rare genetic diseases. We use a variant-level constraint metric (the mutability-adjusted proportion of singletons) to identify constrained functional variant classes near exon-intron junctions and at putative splicing branchpoints. To identify new diagnoses for individuals with unsolved rare diseases in the 100,000 Genomes Project, we identified individuals with de novo single-nucleotide variants near exon-intron boundaries and at putative splicing branchpoints in known disease genes. We identified candidate diagnostic variants through manual phenotype matching and confirmed new molecular diagnoses through clinical variant interpretation and functional RNA studies.

RESULTS

We show that near-splice positions and splicing branchpoints are highly constrained by purifying selection and harbour potentially damaging non-coding variants which are amenable to systematic analysis in sequencing data. From 258 de novo splicing variants in known rare disease genes, we identify 35 new likely diagnoses in probands with an unsolved rare disease. To date, we have confirmed a new diagnosis for six individuals, including four in whom RNA studies were performed.

CONCLUSIONS

Overall, we demonstrate the clinical value of examining non-canonical splicing variants in individuals with unsolved rare diseases.

Short amplicon reverse transcription-polymerase chain reaction detects aberrant splicing in genes with low expression in blood missed by ribonucleic acid sequencing analysis for clinical diagnosis

Use of blood RNA sequencing (RNA-seq) as a splicing analysis tool for clinical interpretation of variants of uncertain significance (VUSs) found via whole-genome and exome sequencing can be difficult for genes that have low expression in the blood due to insufficient read count coverage aligned to specific genes of interest. Here, we present a short amplicon reverse transcription-polymerase chain reaction(RT-PCR) for the detection of genes with low blood expression. Short amplicon RT-PCR, is designed to span three exons where an exon harboring a variant is flanked by one upstream and one downstream exon. We tested short amplicon RT-PCRs for genes that have median transcripts per million (TPM) values less than one according to the genotype-tissue expression database. Median TPM values of genes analyzed in this study are SYN1 = 0.8549, COL1A1 = 0.6275, TCF4 = 0.4009, DSP = .2894, TTN = 0.2851, COL5A2 = 0.1036, TERT = 0.04452, NTRK2 = 0.0344, ABCA4 = 0.00744, PRPH = 0, and WT1 = 0. All these genes show insufficient exon-spanning read coverage in our RNA-seq data to allow splicing analysis. We successfully detected all genes tested except PRPH and WT1. Aberrant splicing was detected in SYN1, TCF4, NTRK2, TTN, and TERT VUSs. Therefore, our results show short amplicon RT-PCR is a useful alternative for the analysis of splicing events in genes with low TPM in blood RNA for clinical diagnostics.

Genetic variants of PKLR are associated with acute pain in sickle cell disease

Acute pain, the most prominent complication of sickle cell disease (SCD), results from vaso-occlusion triggered by sickling of deoxygenated red blood cells (RBCs). Concentration of 2,3-diphosphoglycerate (2,3-DPG) in RBCs promotes deoxygenation by preferentially binding to the low-affinity T conformation of HbS. 2,3-DPG is an intermediate substrate in the glycolytic pathway in which pyruvate kinase (gene PKLR, protein PKR) is a rate-limiting enzyme; variants in PKLR may affect PKR activity, 2,3-DPG levels in RBCs, RBC sickling, and acute pain episodes (APEs). We performed a candidate gene association study using 2 cohorts: 242 adult SCD-HbSS patients and 977 children with SCD-HbSS or SCD-HbSβ0 thalassemia. Seven of 47 PKLR variants evaluated in the adult cohort were associated with hospitalization: intron 4, rs2071053; intron 2, rs8177970, rs116244351, rs114455416, rs12741350, rs3020781, and rs8177964. All 7 variants showed consistent effect directions in both cohorts and remained significant in weighted Fisher's meta-analyses of the adult and pediatric cohorts using P < .0071 as threshold to correct for multiple testing. Allele-specific expression analyses in an independent cohort of 52 SCD adults showed that the intronic variants are likely to influence APE by affecting expression of PKLR, although the causal variant and mechanism are not defined.

RNA and protein expression analysis of HLA-DQB1*03:01:01:21Q allele: A null allele renamed as HLA-DQB1*03:01:01:21N

The characterization of the expression profile of HLA questionable alleles (Q) is clinically relevant in allogeneic hematopoietic stem cell transplantation (HSTC) because an aberrant expression of these alleles could lead to transplantation-related complications. HLA-DQB1*03:01:01:21Q shows a substitution at the donor splice site of intron 3 that potentially could affect the expression of this allele. In order to determine their expression profile at RNA and protein level, we analyzed the presence of the HLA-DQ7 molecule by complement-dependent cytotoxicity test (CDC) and flow cytometry, and their RNA processing by cDNA analyses and sequencing by Sanger methods. Our results reveal that HLA-DQ7 is not detectable by serological methods, this is confirmed by cDNA methods demonstrating the absence of specific HLA-DQB1*03:01:01:21Q mRNA, probably due to an intron 3 retention that creates a premature TGA stop codon, leading to mRNA degradation via nonsense-mediated decay (NMD). These findings demonstrate that the HLA-DQB1*03:01:01:21Q allele is nonexpressed, thus it has been renamed as DQB1*03:01:01:21N.

Re-evaluating the pathogenicity of the c.783+2T>C BAP1 germline variant

BAP1 germline pathogenic sequence variants (PSVs) underlie a unique tumor predisposition syndrome (BAP1-TPDS) associated with an increased lifetime risk for developing primarily pleural and peritoneal mesothelioma and uveal and cutaneous melanoma. Overwhelmingly, BAP1 PSVs are unique, family-specific inactivating variants. We identified seven families, six of Jewish Iraqi origin, harboring an identical BAP1 splice variant (c.783+2T>C), currently assigned a "likely pathogenic" status. Given a nonclassical BAP1-TPDS tumor type clustering and low penetrance in these families, the pathogenicity of this variant was re-evaluated by a combined approach including literature analysis, revised bioinformatics analysis, allelic loss, effect on the transcript, and tumor protein expression patterns. None of the three available tumors showed an allelic loss, there was no discernable effect on alternative splicing based on reverse-transcription polymerase chain reaction, and there was no decrease or loss of somatic protein expression in 2/3 analyzed tumors. This led to assigning a Benign Strong (BS) criteria, BS4, supporting BS3 criteria, and weakening the Pathogenic Supporting (PP) criteria PP5. Combined, these data suggest that this sequence variant should be reclassified as a variant of unknown significance by American College of Medical Genetics (ACMG) criteria.

Splicing analysis of SLC40A1 missense variations and contribution to hemochromatosis type 4 phenotypes

Hemochromatosis type 4, or ferroportin disease, is considered as the second leading cause of primary iron overload after HFE-related hemochromatosis. The disease, which is predominantly associated with missense variations in the SLC40A1 gene, is characterized by wide clinical heterogeneity. We tested the possibility that some of the reported missense mutations, despite their positions within exons, cause splicing defects. Fifty-eight genetic variants were selected from the literature based on two criteria: a precise description of the nucleotide change and individual evidence of iron overload. The selected variants were investigated by different in silico prediction tools and prioritized for midigene splicing assays. Of the 15 variations tested in vitro, only two were associated with splicing changes. We confirm that the c.1402G>A transition (p.Gly468Ser) disrupts the exon 7 donor site, leading to the use of an exonic cryptic splicing site and the generation of a truncated reading frame. We observed, for the first time, that the p.Gly468Ser substitution has no effect on the ferroportin iron export function. We demonstrate alternative splicing of exon 5 in different cell lines and show that the c.430A>G (p.Asn144Asp) variant promotes exon 5 inclusion. This could be part of a gain-of-function mechanism. We conclude that splicing mutations rarely contribute to hemochromatosis type 4 phenotypes. An in-depth investigation of exon 5 auxiliary splicing sequences may help to elucidate the mechanism by which splicing regulatory proteins regulate the production of the full length SLC40A1 transcript and to clarify its physiological importance.

Dual Fluorescence Splicing Reporter Minigene Identifies an Antisense Oligonucleotide to Skip Exon v8 of the CD44 Gene

Splicing reporter minigenes are used in cell-based in vitro splicing studies. Exon skippable antisense oligonucleotide (ASO) has been identified using minigene splicing assays, but these assays include a time- and cost-consuming step of reverse transcription PCR amplification. To make in vitro splicing assay easier, a ready-made minigene (FMv2) amenable to quantitative splicing analysis by fluorescence microscopy was constructed. FMv2 was designed to encode two fluorescence proteins namely, mCherry, a transfection marker and split eGFP, a marker of splicing reaction. The split eGFP was intervened by an artificial intron containing a multicloning site sequence. Expectedly, FMv2 transfected HeLa cells produced not only red mCherry but also green eGFP signals. Transfection of FMv2CD44v8, a modified clone of FMv2 carrying an insertion of CD44 exon v8 in the multicloning site, that was applied to screen exon v8 skippable ASO, produced only red signals. Among seven different ASOs tested against exon v8, ASO#14 produced the highest index of green signal positive cells. Hence, ASO#14 was the most efficient exon v8 skippable ASO. Notably, the well containing ASO#14 was clearly identified among the 96 wells containing randomly added ASOs, enabling high throughput screening. A ready-made FMv2 is expected to contribute to identify exon skippable ASOs.

Comprehensive Analysis of Alternative Splicing Signature in Gastric Cancer Prognosis Based on The Cancer Genome Atlas (TCGA) and SpliceSeq Databases

BACKGROUND Increasing evidence suggests that the alternative splicing (AS) signature plays a role in the carcinogenesis and prognosis of various cancers. However, the prognostic role of AS in gastric cancer is not clear and needs to be clarified. MATERIAL AND METHODS To identify the differentially expressed AS (DEAS) events, we performed a differential expression analysis between normal and tumor tissue. The DEAS event was further applied to construct a prognostic signature by performing univariate Cox regression analysis and least absolute shrinkage and selection operator (LASSO) analysis. The Kaplan-Meier curve analysis and receiver operating characteristic curve (ROC) analysis were used to evaluate the prognostic value of the AS signature. In addition, the network of the splicing events with splicing factors was constructed using the Cytoscape software. RESULTS A total of 30 005 alternative splicing (AS) events with 372 patients were retrieved from the SpliceSeq database and TCGA database. By performing differential expression analysis, a total of 419 alternative splicing events were screened out, including 56 upregulated and 363 downregulated. We further constructed an AS-related prognostic signature by conducting a series bioinformatics analyses. Moreover, we identified that the AS signature could serve as an independent predictor for the prognosis of GC. We also found that AS signature had a more robust and precise efficacy for prognostic prediction in GC patients. Interestingly, the areas under 3- and 5-year survival curves are similar, both of which are greater than 1-year survival curve, suggesting that the long-term predictive accuracy of our prognostic model built upon AS signature is superior. CONCLUSIONS We performed a comprehensive analysis of overall prognostic-associated AS events concerning GC and constructed a prognostic model to predict the long-term prognostic survival outcomes in GC patients. We also developed a network of splicing events with splicing factors to reveal new potential molecular diagnostic biomarkers and therapeutic targets for GC patients.

A synonymous UPF3B variant causing a speech disorder implicates NMD as a regulator of neurodevelopmental disorder gene networks

Loss-of-function mutations of the X-chromosome gene UPF3B cause male neurodevelopmental disorders (NDDs) via largely unknown mechanisms. We investigated initially by interrogating a novel synonymous UPF3B variant in a male with absent speech. In silico and functional studies using cell lines derived from this individual show altered UPF3B RNA splicing. The resulting mRNA species encodes a frame-shifted protein with a premature termination codon (PTC) predicted to elicit degradation via nonsense-mediated mRNA decay (NMD). UPF3B mRNA was reduced in the cell line, and no UPF3B protein was produced, confirming a loss-of-function allele. UPF3B is itself involved in the NMD mechanism which degrades both PTC-bearing mutant transcripts and also many physiological transcripts. RNAseq analysis showed that ~1.6% of mRNAs exhibited altered expression. These mRNA changes overlapped and correlated with those we identified in additional cell lines obtained from individuals harbouring other UPF3B mutations, permitting us to interrogate pathogenic mechanisms of UPF3B-associated NDDs. We identified 102 genes consistently deregulated across all UPF3B mutant cell lines. Of the 51 upregulated genes, 75% contained an NMD-targeting feature, thus identifying high-confidence direct NMD targets. Intriguingly, 22 of the dysregulated genes encoded known NDD genes, suggesting UPF3B-dependent NMD regulates gene networks critical for cognition and behaviour. Indeed, we show that 78.5% of all NDD genes encode a transcript predicted to be targeted by NMD. These data describe the first synonymous UPF3B mutation in a patient with prominent speech and language disabilities and identify plausible mechanisms of pathology downstream of UPF3B mutations involving the deregulation of NDD-gene networks.

Pathogenic deep intronic MTM1 variant activates a pseudo-exon encoding a nonsense codon resulting in severe X-linked myotubular myopathy

X-linked myotubular myopathy (XLMTM) is a severe congenital myopathy characterised by generalised weakness and respiratory insufficiency. XLMTM is associated with pathogenic variants in MTM1; a gene encoding the lipid phosphatase myotubularin. Whole genome sequencing (WGS) of an exome-negative male proband with severe hypotonia, respiratory insufficiency and centralised nuclei on muscle biopsy identified a deep intronic MTM1 variant NG_008199.1(NM_000252.2):c.1468-577A>G, which strengthened a cryptic 5' splice site (A>G substitution at the +5 position). Muscle RNA sequencing was non-diagnostic due to low read depth. Reverse transcription PCR (RT-PCR) of muscle RNA confirmed the c.1468-577A>G variant activates inclusion of a pseudo-exon encoding a premature stop codon into all detected MTM1 transcripts. Western blot analysis establishes deficiency of myotubularin protein, consistent with the severe XLMTM phenotype. We expand the genotypic spectrum of XLMTM and highlight benefits of screening non-coding regions of MTM1 in male probands with phenotypically concordant XLMTM who remain undiagnosed following exome sequencing.

Analysis of transcript-deleterious variants in Mendelian disorders: implications for RNA-based diagnostics

Background

At least 50% of patients with suspected Mendelian disorders remain undiagnosed after whole-exome sequencing (WES), and the extent to which non-coding variants that are not captured by WES contribute to this fraction is unclear. Whole transcriptome sequencing is a promising supplement to WES, although empirical data on the contribution of RNA analysis to the diagnosis of Mendelian diseases on a large scale are scarce.

Results

Here, we describe our experience with transcript-deleterious variants (TDVs) based on a cohort of 5647 families with suspected Mendelian diseases. We first interrogate all families for which the respective Mendelian phenotype could be mapped to a single locus to obtain an unbiased estimate of the contribution of TDVs at 18.9%. We examine the entire cohort and find that TDVs account for 15% of all “solved” cases. We compare the results of RT-PCR to in silico prediction. Definitive results from RT-PCR are obtained from blood-derived RNA for the overwhelming majority of variants (84.1%), and only a small minority (2.6%) fail analysis on all available RNA sources (blood-, skin fibroblast-, and urine renal epithelial cells-derived), which has important implications for the clinical application of RNA-seq. We also show that RNA analysis can establish the diagnosis in 13.5% of 155 patients who had received “negative” clinical WES reports. Finally, our data suggest a role for TDVs in modulating penetrance even in otherwise highly penetrant Mendelian disorders.

Conclusions

Our results provide much needed empirical data for the impending implementation of diagnostic RNA-seq in conjunction with genome sequencing.

Rational Design of an Activatable Reporter for Quantitative Imaging of RNA Aberrant Splicing In Vivo

Pre-mRNA splicing, the process of removing introns from pre-mRNA and the arrangement of exons to produce mature transcripts, is a crucial step in the expression of most eukaryote genes. However, the splicing kinetics remain poorly characterized in living cells, mainly because current methods cannot provide the dynamic information of splicing events. Here, we developed a genetically encoded bioluminescence reporter for real-time imaging of the pre-mRNA splicing process in living subjects. We showed that the bioluminescence reporter is able to visualize the pre-mRNA aberrant splicing process in living cells in a dose- and time-dependent manner. Moreover, this reporter could provide quantitative and longitudinal information of splicing activity in response to exogenous splicing inhibitors in living animals. Our data suggest that this activatable reporter could serve as a promising tool for the high-throughput screening of splicing modulators, which would facilitate the drug development for human diseases caused by the abnormal splicing of mRNA.

Next Generation Sequencing in Newborn Screening in the United Kingdom National Health Service

Next generation DNA sequencing (NGS) has the potential to improve the diagnostic and prognostic utility of newborn screening programmes. This study assesses the feasibility of automating NGS on dried blood spot (DBS) DNA in a United Kingdom National Health Service (UK NHS) laboratory. An NGS panel targeting the entire coding sequence of five genes relevant to disorders currently screened for in newborns in the UK was validated on DBS DNA. An automated process for DNA extraction, NGS and bioinformatics analysis was developed. The process was tested on DBS to determine feasibility, turnaround time and cost. The analytical sensitivity of the assay was 100% and analytical specificity was 99.96%, with a mean 99.5% concordance of variant calls between DBS and venous blood samples in regions with ≥30× coverage (96.8% across all regions; all variant calls were single nucleotide variants (SNVs), with indel performance not assessed). The pipeline enabled processing of up to 1000 samples a week with a turnaround time of four days from receipt of sample to reporting. This study concluded that it is feasible to automate targeted NGS on routine DBS samples in a UK NHS laboratory setting, but it may not currently be cost effective as a first line test.

RNA-Seq Perspectives to Improve Clinical Diagnosis

In recent years, high-throughput next-generation sequencing technology has allowed a rapid increase in diagnostic capacity and precision through different bioinformatics processing algorithms, tools, and pipelines. The identification, annotation, and classification of sequence variants within different target regions are now considered a gold standard in clinical genetic diagnosis. However, this procedure lacks the ability to link regulatory events such as differential splicing to diseases. RNA-seq is necessary in clinical routine in order to interpret and detect among others splicing events and splicing variants, as it would increase the diagnostic rate by up to 10-35%. The transcriptome has a very dynamic nature, varying according to tissue type, cellular conditions, and environmental factors that may affect regulatory events such as splicing and the expression of genes or their isoforms. RNA-seq offers a robust technical analysis of this complexity, but it requires a profound knowledge of computational/statistical tools that may need to be adjusted depending on the disease under study. In this article we will cover RNA-seq analyses best practices applied to clinical routine, bioinformatics procedures, and present challenges of this approach.

Aberrant expression of alternative splicing variants in multiple sclerosis - A systematic review

OBJECTIVE

Alternative splicing is an important form of RNA processing that affects nearly all human genes. The differential expression of specific transcript and protein isoforms holds the potential of novel biomarkers for complex diseases. In this systematic review, we compiled the existing literature on aberrant alternative splicing events in multiple sclerosis (MS).

METHODS

A systematic literature search in the PubMed database was carried out and supplemented by screening the reference lists of the identified articles. We selected only MS-related original research studies which compared the levels of different isoforms of human protein-coding genes. A narrative synthesis of the research findings was conducted. Additionally, we performed a case-control analysis using high-density transcriptome microarray data to reevaluate the genes that were examined in the reviewed studies.

RESULTS

A total of 160 records were screened. Of those, 36 studies from the last two decades were included. Most commonly, peripheral blood samples were analyzed (32 studies), and PCR-based techniques were usually employed (27 studies) for measuring the expression of selected genes. Two studies used an exploratory genome-wide approach. Overall, 27 alternatively spliced genes were investigated. Nine of these genes appeared in at least two studies (CD40, CFLAR, FOXP3, IFNAR2, IL7R, MOG, PTPRC, SP140 and TNFRSF1A). The microarray data analysis confirmed differential alternative pre-mRNA splicing for 19 genes.

CONCLUSIONS

An altered RNA processing of genes mediating immune signaling pathways has been repeatedly implicated in MS. The analysis of individual exon-level expression patterns is stimulated by the advancement of transcriptome profiling technologies. In particular, the examination of genes encoded in MS-associated genetic regions may provide important insights into the pathogenesis of the disease and help to identify new biomarkers.