Predicting the impact of rare variants on RNA splicing in CAGI6.

Variants which disrupt splicing are a frequent cause of rare disease that have been under-ascertained clinically. Accurate and efficient methods to predict a variant's impact on splicing are needed to interpret the growing number of variants of unknown significance (VUS) identified by exome and genome sequencing. Here, we present the results of the CAGI6 Splicing VUS challenge, which invited predictions of the splicing impact of 56 variants ascertained clinically and functionally validated to determine splicing impact. The performance of 12 prediction methods, along with SpliceAI and CADD, was compared on the 56 functionally validated variants. The maximum accuracy achieved was 82% from two different approaches, one weighting SpliceAI scores by minor allele frequency, and one applying the recently published Splicing Prediction Pipeline (SPiP). SPiP performed optimally in terms of sensitivity, while an ensemble method combining multiple prediction tools and information from databases exceeded all others for specificity. Several challenge methods equalled or exceeded the performance of SpliceAI, with ultimate choice of prediction method likely to depend on experimental or clinical aims. One quarter of the variants were incorrectly predicted by at least 50% of the methods, highlighting the need for further improvements to splicing prediction methods for successful clinical application.

Noncoding variants are a rare cause of recessive developmental disorders in trans with coding variants.

PURPOSE: Identifying pathogenic noncoding variants is challenging. A single protein-altering variant is often identified in a recessive gene in individuals with developmental disorders (DD), but the prevalence of pathogenic noncoding "second hits" in trans with these is unknown. METHODS: In 4073 genetically undiagnosed rare-disease trio probands from the 100,000 Genomes project, we identified rare heterozygous protein-altering variants in recessive DD-associated genes. We identified rare noncoding variants on the other haplotype in introns, untranslated regions, promoters, and candidate enhancer regions. We clinically evaluated the top candidates for phenotypic fit and performed functional testing where possible. RESULTS: We identified 3761 rare heterozygous loss-of-function or ClinVar pathogenic variants in recessive DD-associated genes in 2430 probands. For 1366 (36.3%) of these, we identified at least 1 rare noncoding variant in trans. Bioinformatic filtering and clinical review, revealed 7 to be a good clinical fit. After detailed characterization, we identified likely diagnoses for 3 probands (in GAA, NPHP3, and PKHD1) and candidate diagnoses in a further 3 (PAH, LAMA2, and IGHMBP2). CONCLUSION: We developed a systematic approach to uncover new diagnoses involving compound heterozygous coding/noncoding variants and conclude that this mechanism is likely to be a rare cause of DDs.

Extending the phenotypes associated with TRIO gene variants in a cohort of 25 patients and review of the literature.

The TRIO gene encodes a rho guanine exchange factor, the function of which is to exchange GDP to GTP, and hence to activate Rho GTPases, and has been described to impact neurodevelopment. Specific genotype-to-phenotype correlations have been established previously describing striking differentiating features seen in variants located in specific domains of the TRIO gene that are associated with opposite effects on RAC1 activity. Currently, 32 cases with a TRIO gene alteration have been published in the medical literature. Here, we report an additional 25, previously unreported individuals who possess heterozygous TRIO variants and we review the literature. In addition, functional studies were performed on the c.4394A > G (N1465S) and c.6244-2A > G TRIO variants to provide evidence for their pathogenicity. Variants reported by the current study include missense variants, truncating nonsense variants, and an intragenic deletion. Clinical features were previously described and included developmental delay, learning difficulties, microcephaly, macrocephaly, seizures, behavioral issues (aggression, stereotypies), skeletal problems including short, tapering fingers and scoliosis, dental problems (overcrowding/delayed eruption), and variable facial features. Here, we report clinical features that have not been described previously, including specific structural brain malformations such as abnormalities of the corpus callosum and ventriculomegaly, additional psychological and dental issues along with a more recognizable facial gestalt linked to the specific domains of the TRIO gene and the effect of the variant upon the function of the encoded protein. This current study further strengthens the genotype-to-phenotype correlation that was previously established and extends the range of phenotypes to include structural brain abnormalities, additional skeletal, dental, and psychiatric issues.

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.

Blood RNA analysis can increase clinical diagnostic rate and resolve variants of uncertain significance.

PURPOSE: Diagnosis of genetic disorders is hampered by large numbers of variants of uncertain significance (VUSs) identified through next-generation sequencing. Many such variants may disrupt normal RNA splicing. We examined effects on splicing of a large cohort of clinically identified variants and compared performance of bioinformatic splicing prediction tools commonly used in diagnostic laboratories. METHODS: Two hundred fifty-seven variants (coding and noncoding) were referred for analysis across three laboratories. Blood RNA samples underwent targeted reverse transcription polymerase chain reaction (RT-PCR) analysis with Sanger sequencing of PCR products and agarose gel electrophoresis. Seventeen samples also underwent transcriptome-wide RNA sequencing with targeted splicing analysis based on Sashimi plot visualization. Bioinformatic splicing predictions were obtained using Alamut, HSF 3.1, and SpliceAI software. RESULTS: Eighty-five variants (33%) were associated with abnormal splicing. The most frequent abnormality was upstream exon skipping (39/85 variants), which was most often associated with splice donor region variants. SpliceAI had greatest accuracy in predicting splicing abnormalities (0.91) and outperformed other tools in sensitivity and specificity. CONCLUSION: Splicing analysis of blood RNA identifies diagnostically important splicing abnormalities and clarifies functional effects of a significant proportion of VUSs. Bioinformatic predictions are improving but still make significant errors. RNA analysis should therefore be routinely considered in genetic disease diagnostics.

Correction: Blood RNA analysis can increase clinical diagnostic rate and resolve variants of uncertain significance.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

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.

Identification of diagnostic candidates in Mendelian disorders using an RNA sequencing-centric approach.

BACKGROUND: RNA sequencing (RNA-seq) is increasingly being used as a complementary tool to DNA sequencing in diagnostics where DNA analysis has been uninformative. RNA-seq enables the identification of aberrant splicing and aberrant gene expression, improving the interpretation of variants of unknown significance (VUSs), and provides the opportunity to scan the transcriptome for aberrant splicing and expression in relevant genes that may be the cause of a patient's phenotype. This work aims to investigate the feasibility of generating new diagnostic candidates in patients without a previously reported VUS using an RNA-seq-centric approach. METHODS: We systematically assessed the transcriptomic profiles of 86 patients with suspected Mendelian disorders, 38 of whom had no candidate sequence variant, using RNA from blood samples. Each VUS was visually inspected to search for splicing abnormalities. Once aberrant splicing was identified in cases with VUS, multiple open-source alternative splicing tools were used to investigate if they would identify what was observed in IGV. Expression outliers were detected using OUTRIDER. Diagnoses in cases without a VUS were explored using two separate strategies. RESULTS: RNA-seq allowed us to assess 71% of VUSs, detecting aberrant splicing in 14/48 patients with a VUS. We identified four new diagnoses by detecting novel aberrant splicing events in patients with no candidate sequence variants from prior DNA testing (n = 32) or where the candidate VUS did not affect splicing (n = 23). An additional diagnosis was made through the detection of skewed X-inactivation. CONCLUSION: This work demonstrates the utility of an RNA-centric approach in identifying novel diagnoses in patients without candidate VUSs. It underscores the utility of blood-based RNA analysis in improving diagnostic yields and highlights optimal approaches for such analyses.

Temporal Whole-Transcriptomic Analysis of Characterized In Vitro and Ex Vivo Primary Nasal Epithelia.

Air-liquid interface (ALI) cell culture of primary airway progenitors enables the differentiation and recapitulation of a pseudostratified epithelium in vitro, providing a highly useful tool for researching respiratory health and disease. Previous studies into gene expression in ALI-cultures compared to ex vivo nasal brushings have been limited in the number of time-points and/or the number of genes studied. In this study physiological and global transcriptomic changes were assessed in an extended in vitro 63-day human healthy nasal epithelium ALI-culture period and compared to ex vivo nasal brushing samples. Ex vivo nasal brushing samples formed distinct transcriptome clusters to in vitro ALI-cultured nasal epithelia, with from day 14 onwards ALI samples best matching the ex vivo samples. Immune response regulation genes were not expressed in the in vitro ALI-culture compared to the ex vivo nasal brushing samples, likely because the in vitro cultures lack an airway microbiome, lack airborne particles stimulation, or did not host an immune cell component. This highlights the need for more advanced co-cultures with immune cell representation to better reflect the physiological state. During the first week of ALI-culture genes related to metabolism and proliferation were increased. By the end of week 1 epithelial cell barrier function plateaued and multiciliated cell differentiation started, although widespread ciliation was not complete until day 28. These results highlight that time-points at which ALI-cultures are harvested for research studies needs to be carefully considered to suit the purpose of investigation (transcriptomic and/or functional analysis).

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.

Biallelic variants in COPB1 cause a novel, severe intellectual disability syndrome with cataracts and variable microcephaly.

BACKGROUND: Coat protein complex 1 (COPI) is integral in the sorting and retrograde trafficking of proteins and lipids from the Golgi apparatus to the endoplasmic reticulum (ER). In recent years, coat proteins have been implicated in human diseases known collectively as "coatopathies". METHODS: Whole exome or genome sequencing of two families with a neuro-developmental syndrome, variable microcephaly and cataracts revealed biallelic variants in COPB1, which encodes the beta-subunit of COPI (ß-COP). To investigate Family 1's splice donor site variant, we undertook patient blood RNA studies and CRISPR/Cas9 modelling of this variant in a homologous region of the Xenopus tropicalis genome. To investigate Family 2's missense variant, we studied cellular phenotypes of human retinal epithelium and embryonic kidney cell lines transfected with a COPB1 expression vector into which we had introduced Family 2's mutation. RESULTS: We present a new recessive coatopathy typified by severe developmental delay and cataracts and variable microcephaly. A homozygous splice donor site variant in Family 1 results in two aberrant transcripts, one of which causes skipping of exon 8 in COPB1 pre-mRNA, and a 36 amino acid in-frame deletion, resulting in the loss of a motif at a small interaction interface between ß-COP and ß'-COP. Xenopus tropicalis animals with a homologous mutation, introduced by CRISPR/Cas9 genome editing, recapitulate features of the human syndrome including microcephaly and cataracts. In vitro modelling of the COPB1 c.1651T>G p.Phe551Val variant in Family 2 identifies defective Golgi to ER recycling of this mutant ß-COP, with the mutant protein being retarded in the Golgi. CONCLUSIONS: This adds to the growing body of evidence that COPI subunits are essential in brain development and human health and underlines the utility of exome and genome sequencing coupled with Xenopus tropicalis CRISPR/Cas modelling for the identification and characterisation of novel rare disease genes.

Comparison of in silico strategies to prioritize rare genomic variants impacting RNA splicing for the diagnosis of genomic disorders.

The development of computational methods to assess pathogenicity of pre-messenger RNA splicing variants is critical for diagnosis of human disease. We assessed the capability of eight algorithms, and a consensus approach, to prioritize 249 variants of uncertain significance (VUSs) that underwent splicing functional analyses. The capability of algorithms to differentiate VUSs away from the immediate splice site as being 'pathogenic' or 'benign' is likely to have substantial impact on diagnostic testing. We show that SpliceAI is the best single strategy in this regard, but that combined usage of tools using a weighted approach can increase accuracy further. We incorporated prioritization strategies alongside diagnostic testing for rare disorders. We show that 15% of 2783 referred individuals carry rare variants expected to impact splicing that were not initially identified as 'pathogenic' or 'likely pathogenic'; one in five of these cases could lead to new or refined diagnoses.