AG-exclusion zone revisited: Lessons to learn from 91 intronic NF1 3' splice site mutations outside the canonical AG-dinucleotides

RNA variant assessment using transactivation and transdifferentiation

Understanding the impact of splicing and nonsense variants on RNA is crucial for the resolution of variant classification as well as their suitability for precision medicine interventions. This is primarily enabled through RNA studies involving transcriptomics followed by targeted assays using RNA isolated from clinically accessible tissues (CATs) such as blood or skin of affected individuals. Insufficient disease gene expression in CATs does however pose a major barrier to RNA based investigations, which we show is relevant to 1,436 Mendelian disease genes. We term these "silent" Mendelian genes (SMGs), the largest portion (36%) of which are associated with neurological disorders. We developed two approaches to induce SMG expression in human dermal fibroblasts (HDFs) to overcome this limitation, including CRISPR-activation-based gene transactivation and fibroblast-to-neuron transdifferentiation. Initial transactivation screens involving 40 SMGs stimulated our development of a highly multiplexed transactivation system culminating in the 6- to 90,000-fold induction of expression of 20/20 (100%) SMGs tested in HDFs. Transdifferentiation of HDFs directly to neurons led to expression of 193/516 (37.4%) of SMGs implicated in neurological disease. The magnitude and isoform diversity of SMG expression following either transactivation or transdifferentiation was comparable to clinically relevant tissues. We apply transdifferentiation and/or gene transactivation combined with short- and long-read RNA sequencing to investigate the impact that variants in USH2A, SCN1A, DMD, and PAK3 have on RNA using HDFs derived from affected individuals. Transactivation and transdifferentiation represent rapid, scalable functional genomic solutions to investigate variants impacting SMGs in the patient cell and genomic context.

A Comprehensive Overview of NF1 Mutations in Iranian Patients

Neurofibromatosis type 1 (NF1) is a genetic disorder caused by mutations in the NF1 gene. This disorder shows nearly complete penetrance and high phenotypic variability. We used the whole-exome sequencing technique to identify mutations in 32 NF1 cases from 22 Iranian families. A total of 31 variants, including 30 point mutations and one large deletion, were detected. In eight cases, variants were inherited, while they were sporadic in the remaining. Seven novel variants, including c.5576 T > G, c.6658_6659insC, c.2322dupT, c.92_93insAA, c.4360C > T, c.3814C > T, and c.4565_4566delinsC, were identified. The current study is the largest in terms of the sample size of Iranian NF1 cases with identified mutations. The results can broaden the spectrum of NF1 mutations and facilitate the process of genetic counseling in the affected families.

Splice-Switching Antisense Oligonucleotides Correct Phenylalanine Hydroxylase Exon 11 Skipping Defects and Rescue Enzyme Activity in Phenylketonuria

The PAH gene encodes the hepatic enzyme phenylalanine hydroxylase (PAH), and its deficiency, known as phenylketonuria (PKU), leads to neurotoxic high levels of phenylalanine. PAH exon 11 is weakly defined, and several missense and intronic variants identified in patients affect the splicing process. Recently, we identified a novel intron 11 splicing regulatory element where U1snRNP binds, participating in exon 11 definition. In this work, we describe the implementation of an antisense strategy targeting intron 11 sequences to correct the effect of PAH mis-splicing variants. We used an in vitro assay with minigenes and identified splice-switching antisense oligonucleotides (SSOs) that correct the exon skipping defect of PAH variants c.1199+17G>A, c.1199+20G>C, c.1144T>C, and c.1066-3C>T. To examine the functional rescue induced by the SSOs, we generated a hepatoma cell model with variant c.1199+17G>A using CRISPR/Cas9. The edited cell line reproduces the exon 11 skipping pattern observed from minigenes, leading to reduced PAH protein levels and activity. SSO transfection results in an increase in exon 11 inclusion and corrects PAH deficiency. Our results provide proof of concept of the potential therapeutic use of a single SSO for different exonic and intronic splicing variants causing PAH exon 11 skipping in PKU.

Genome-wide detection of human intronic AG-gain variants located between splicing branchpoints and canonical splice acceptor sites

Human genetic variants that introduce an AG into the intronic region between the branchpoint (BP) and the canonical splice acceptor site (ACC) of protein-coding genes can disrupt pre-mRNA splicing. Using our genome-wide BP database, we delineated the BP-ACC segments of all human introns and found extreme depletion of AG/YAG in the [BP+8, ACC-4] high-risk region. We developed AGAIN as a genome-wide computational approach to systematically and precisely pinpoint intronic AG-gain variants within the BP-ACC regions. AGAIN identified 350 AG-gain variants from the Human Gene Mutation Database, all of which alter splicing and cause disease. Among them, 74% created new acceptor sites, whereas 31% resulted in complete exon skipping. AGAIN also predicts the protein-level products resulting from these two consequences. We performed AGAIN on our exome/genomes database of patients with severe infectious diseases but without known genetic etiology and identified a private homozygous intronic AG-gain variant in the antimycobacterial gene SPPL2A in a patient with mycobacterial disease. AGAIN also predicts a retention of six intronic nucleotides that encode an in-frame stop codon, turning AG-gain into stop-gain. This allele was then confirmed experimentally to lead to loss of function by disrupting splicing. We further showed that AG-gain variants inside the high-risk region led to misspliced products, while those outside the region did not, by two case studies in genes STAT1 and IRF7. We finally evaluated AGAIN on our 14 paired exome-RNAseq samples and found that 82% of AG-gain variants in high-risk regions showed evidence of missplicing. AGAIN is publicly available from https://hgidsoft.rockefeller.edu/AGAIN and https://github.com/casanova-lab/AGAIN.

Analysis of 200 unrelated individuals with a constitutional NF1 deep intronic pathogenic variant reveals that variants flanking the alternatively spliced NF1 exon 31 [23a] cause a classical neurofibromatosis type 1 phenotype while altering predominantly NF1 isoform type II

Neurofibromatosis type 1 results from loss-of-function NF1 pathogenic variants (PVs). Up to 30% of all NF1 PVs disrupt mRNA splicing, including deep intronic variants. Here, we retrospectively investigated the spectrum of NF1 deep intronic PVs in a cohort of 8,090 unrelated individuals from the University of Alabama at Birmingham (UAB) dataset with a molecularly confirmed neurofibromatosis type 1. All variants were identified through their effect on the NF1 transcript, followed by variant characterization at the DNA-level. A total of 68 distinct variants, which were ≥ 20 nucleotides away from the closest exon-intron junction, were identified in 2.5% unrelated individuals with NF1 (200/8,090). Nine different pathogenic splice variants, identified in 20 probands, led to exonization of different parts of intron 30 [23.2] or 31 [23a]. The two major NF1 transcript isoforms, distinguished by the absence (type I) or presence (type II) of the alternatively spliced cassette exon 31 [23a], are equally expressed in blood in control individuals without NF1 or NF1-affected individuals carrying their PV not in the introns flanking exon 31 [23a]. By fragment and cloning analysis we demonstrated that the exonization of intron 31 [23a] sequences due to deep intronic PV predominantly affects the NF1 isoform II. Seven additional (likely) pathogenic NF1 deep intronic variants not observed in the UAB dataset were found by classification of 36 variants identified by a literature search. Hence, the unique list of these 75 deep intronic (likely) PVs should be included in any comprehensive NF1 testing strategy.

High-yield identification of pathogenic NF1 variants by skin fibroblast transcriptome screening after apparently normal diagnostic DNA testing

Neurofibromatosis type 1 (NF1) is caused by inactivating mutations in NF1. Due to the size, complexity, and high mutation rate at the NF1 locus, the identification of causative variants can be challenging. To obtain a molecular diagnosis in 15 individuals meeting diagnostic criteria for NF1, we performed transcriptome analysis (RNA-seq) on RNA obtained from cultured skin fibroblasts. In each case, routine molecular DNA diagnostics had failed to identify a disease-causing variant in NF1. A pathogenic variant or abnormal mRNA splicing was identified in 13 cases: 6 deep intronic variants and 2 transposon insertions causing noncanonical splicing, 3 postzygotic changes, 1 branch point mutation and, in 1 case, abnormal splicing for which the responsible DNA change remains to be identified. These findings helped resolve the molecular findings for an additional 17 individuals in multiple families with NF1, demonstrating the utility of skin-fibroblast-based transcriptome analysis for molecular diagnostics. RNA-seq improves mutation detection in NF1 and provides a powerful complementary approach to DNA-based methods. Importantly, our approach is applicable to other genetic disorders, particularly those caused by a wide variety of variants in a limited number of genes and specifically for individuals in whom routine molecular DNA diagnostics did not identify the causative variant.

Combining genetic constraint with predictions of alternative splicing to prioritize deleterious splicing in rare disease studies

Background

Despite numerous molecular and computational advances, roughly half of patients with a rare disease remain undiagnosed after exome or genome sequencing. A particularly challenging barrier to diagnosis is identifying variants that cause deleterious alternative splicing at intronic or exonic loci outside of canonical donor or acceptor splice sites.

Results

Several existing tools predict the likelihood that a genetic variant causes alternative splicing. We sought to extend such methods by developing a new metric that aids in discerning whether a genetic variant leads to deleterious alternative splicing. Our metric combines genetic variation in the Genome Aggregate Database with alternative splicing predictions from SpliceAI to compare observed and expected levels of splice-altering genetic variation. We infer genic regions with significantly less splice-altering variation than expected to be constrained. The resulting model of regional splicing constraint captures differential splicing constraint across gene and exon categories, and the most constrained genic regions are enriched for pathogenic splice-altering variants. Building from this model, we developed ConSpliceML. This ensemble machine learning approach combines regional splicing constraint with multiple per-nucleotide alternative splicing scores to guide the prediction of deleterious splicing variants in protein-coding genes. ConSpliceML more accurately distinguishes deleterious and benign splicing variants than state-of-the-art splicing prediction methods, especially in “cryptic” splicing regions beyond canonical donor or acceptor splice sites.

Conclusion

Integrating a model of genetic constraint with annotations from existing alternative splicing tools allows ConSpliceML to prioritize potentially deleterious splice-altering variants in studies of rare human diseases.

Genome-wide detection of human variants that disrupt intronic branchpoints

Pre-messenger RNA splicing is initiated with the recognition of a single-nucleotide intronic branchpoint (BP) within a BP motif by spliceosome elements. Forty-eight rare variants in 43 human genes have been reported to alter splicing and cause disease by disrupting BP. However, until now, no computational approach was available to efficiently detect such variants in massively parallel sequencing data. We established a comprehensive human genome-wide BP database by integrating existing BP data and generating new BP data from RNA sequencing of lariat debranching enzyme DBR1-mutated patients and from machine-learning predictions. We characterized multiple features of BP in major and minor introns and found that BP and BP-2 (two nucleotides upstream of BP) positions exhibit a lower rate of variation in human populations and higher evolutionary conservation than the intronic background, while being comparable to the exonic background. We developed BPHunter as a genome-wide computational approach to systematically and efficiently detect intronic variants that may disrupt BP recognition. BPHunter retrospectively identified 40 of the 48 known pathogenic BP variants, in which we summarized a strategy for prioritizing BP variant candidates. The remaining eight variants all create AG-dinucleotides between the BP and acceptor site, which is the likely reason for missplicing. We demonstrated the practical utility of BPHunter prospectively by using it to identify a novel germline heterozygous BP variant of STAT2 in a patient with critical COVID-19 pneumonia and a novel somatic intronic 59-nucleotide deletion of ITPKB in a lymphoma patient, both of which were validated experimentally. BPHunter is publicly available from https://hgidsoft.rockefeller.edu/BPHunter and https://github.com/casanova-lab/BPHunter .

Comparison of In Silico Tools for Splice-Altering Variant Prediction Using Established Spliceogenic Variants: An End-User’s Point of View

Assessing the impact of variants of unknown significance on splicing has become a critical issue and a bottleneck, especially with the widespread implementation of whole-genome or exome sequencing. Although multiple in silico tools are available, the interpretation and application of these tools are difficult and practical guidelines are still lacking. A streamlined decision-making process can facilitate the downstream RNA analysis in a more efficient manner. Therefore, we evaluated the performance of 8 in silico tools (Splice Site Finder, MaxEntScan, Splice-site prediction by neural network, GeneSplicer, Human Splicing Finder, SpliceAI, Splicing Predictions in Consensus Elements, and SpliceRover) using 114 NF1 spliceogenic variants, experimentally validated at the mRNA level. The change in the predicted score incurred by the variant of the nearest wild-type splice site was analyzed, and for type II, III, and IV splice variants, the change in the prediction score of de novo or cryptic splice site was also analyzed. SpliceAI and SpliceRover, tools based on deep learning, outperformed all other tools, with AUCs of 0.972 and 0.924, respectively. For de novo and cryptic splice sites, SpliceAI outperformed all other tools and showed a sensitivity of 95.7% at an optimal cut-off of 0.02 score change. Our results show that deep learning algorithms, especially those of SpliceAI, are validated at a significantly higher rate than other in silico tools for clinically relevant NF1 variants. This suggests that deep learning algorithms outperform traditional probabilistic approaches and classical machine learning tools in predicting the de novo and cryptic splice sites.

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.

Genotype-Phenotype Correlations in Neurofibromatosis Type 1: Identification of Novel and Recurrent NF1 Gene Variants and Correlations with Neurocognitive Phenotype

Neurofibromatosis type 1 (NF1) is one of the most common genetic tumor predisposition syndrome, caused by mutations in the NF1. To date, few genotype-phenotype correlations have been discerned in NF1, due to a highly variable clinical presentation. We aimed to study the molecular spectrum of NF1 and genotype-phenotype correlations in a monocentric study cohort of 85 NF1 patients (20 relatives, 65 sporadic cases). Clinical data were collected at the time of the mutation analysis and reviewed for accuracy in this investigation. An internal phenotypic categorization was applied. The 94% of the patients enrolled showed a severe phenotype with at least one systemic complication and a wide range of associated malignancies. Spine deformities were the most common complications in this cohort. We also reported 66 different NF1 mutations, of which 7 are novel mutations. Correlation analysis identified a slight significant inverse correlation between age at diagnosis and delayed acquisition of psychomotor skills with residual multi-domain cognitive impairment. Odds ratio with 95% confidence interval showed a higher prevalence of learning disabilities in patients carrying frameshift mutations. Overall, our results aim to offer an interesting contribution to studies on the genotype–phenotype of NF1 and in genetic management and counselling.

Mutation analysis reveals novel and known mutations in SAG gene in first two Egyptian families with Oguchi disease

Background

Oguchi disease is a rare type of congenital stationary night blindness associated with an abnormal fundus appearance. It is inherited in an autosomal recessive manner where two types exist according to the gene affected; type 1 associated with S-antigen (SAG) gene mutations and type 2 associated with rhodopsin kinase (GRK1) gene mutations.

Purpose

The aim of this work was to describe the clinical and genetic findings of the first two reported families of Oguchi disease in Egypt and African region.

Methods

Four members of two consanguineous Egyptian families with history of night blindness since childhood underwent complete ophthalmological examination, standard automated static perimetry, fundus color photography, fundus autofluorescence (FAF), fundus fluorescein angiography (FFA) in light-adapted state and spectral-domain optical coherence tomography (SD-OCT) of both the macula and the optic nerve head as well as central corneal thickness with repeated fundus photography following prolonged dark adaptation. Mutation screening of 7 coding exons of GRK1 gene and 15 coding exons of SAG gene as well as some flanking regions were performed using Sanger sequencing technique. The variants were tested for pathogenicity using different in silico functional analysis tools.

Results

The clinical examination and investigations confirmed Oguchi disease phenotype. One patient showed p.R193* (c.577C > T) which is a previously reported SAG gene mutation in a homozygous form. The other three patients from a different family showed (c.649–1 G > C), a novel canonical splice site SAG gene mutation in a homozygous form.

Conclusion

The identification of the novel canonical splice site SAG gene variant in three members of the same family with clinically confirmed Oguchi disease reinforces its pathogenicity. A fourth patient from another family carried a previously reported mutation in the same gene. SAG gene variants may be the underlying genetic cause for Oguchi disease in Egypt. Our findings have expanded the spectrum of Oguchi disease-associated mutations in SAG gene and may serve as a basis for genetic diagnosis for Oguchi disease.

Identification of an elusive spliceogenic MYBPC3 variant in an otherwise genotype-negative hypertrophic cardiomyopathy pedigree

The finding of a genotype-negative hypertrophic cardiomyopathy (HCM) pedigree with several affected members indicating a familial origin of the disease has driven this study to discover causative gene variants. Genetic testing of the proband and subsequent family screening revealed the presence of a rare variant in the MYBPC3 gene, c.3331−26T>G in intron 30, with evidence supporting cosegregation with the disease in the family. An analysis of potential splice-altering activity using several splicing algorithms consistently yielded low scores. Minigene expression analysis at the mRNA and protein levels revealed that c.3331−26T>G is a spliceogenic variant with major splice-altering activity leading to undetectable levels of properly spliced transcripts or the corresponding protein. Minigene and patient mRNA analyses indicated that this variant induces complete and partial retention of intron 30, which was expected to lead to haploinsufficiency in carrier patients. As most spliceogenic MYBPC3 variants, c.3331−26T>G appears to be non-recurrent, since it was identified in only two additional unrelated probands in our large HCM cohort. In fact, the frequency analysis of 46 known splice-altering MYBPC3 intronic nucleotide substitutions in our HCM cohort revealed 9 recurrent and 16 non-recurrent variants present in a few probands (≤ 4), while 21 were not detected. The identification of non-recurrent elusive MYBPC3 spliceogenic variants that escape detection by in silico algorithms represents a challenge for genetic diagnosis of HCM and contributes to solving a fraction of genotype-negative HCM cases.

Intron mutations and early transcription termination in Duchenne and Becker muscular dystrophy

DMD pathogenic variants for Duchenne and Becker muscular dystrophy are detectable with high sensitivity by standard clinical exome analyses of genomic DNA. However, up to 7% of DMD mutations are deep intronic and analysis of muscle-derived RNA is an important diagnostic step for patients who have negative genomic testing but abnormal dystrophin expression in muscle. In this study, muscle biopsies were evaluated from 19 patients with clinical features of a dystrophinopathy, but negative clinical DMD mutation analysis. Reverse transcription-polymerase chain reaction or high-throughput RNA sequencing methods identified 19 mutations with one of three pathogenic pseudoexon types: deep intronic point mutations, deletions or insertions, and translocations. In association with point mutations creating intronic splice acceptor sites, we observed the first examples of DMD pseudo 3'-terminal exon mutations causing high efficiency transcription termination within introns. This connection between splicing and premature transcription termination is reminiscent of U1 snRNP-mediating telescripting in sustaining RNA polymerase II elongation across large genes, such as DMD. We propose a novel classification of three distinct types of mutations identifiable by muscle RNA analysis, each of which differ in potential treatment approaches. Recognition and appropriate characterization may lead to therapies directed toward full-length dystrophin expression for some patients.

A spotter's guide to SNPtic exons: The common splice variants underlying some SNP-phenotype correlations

BACKGROUND

Cryptic exons are typically characterised as deleterious splicing aberrations caused by deep intronic mutations. However, low-level splicing of cryptic exons is sometimes observed in the absence of any pathogenic mutation. Five recent reports have described how low-level splicing of cryptic exons can be modulated by common single-nucleotide polymorphisms (SNPs), resulting in phenotypic differences amongst different genotypes.

METHODS

We sought to investigate whether additional 'SNPtic' exons may exist, and whether these could provide an explanatory mechanism for some of the genotype-phenotype correlations revealed by genome-wide association studies. We thoroughly searched the literature for reported cryptic exons, cross-referenced their genomic coordinates against the dbSNP database of common SNPs, then screened out SNPs with no reported phenotype associations.

RESULTS

This method discovered five probable SNPtic exons in the genes APC, FGB, GHRL, MYPBC3 and OTC. For four of these five exons, we observed that the phenotype associated with the SNP was compatible with the predicted splicing effect of the nucleotide change, whilst the fifth (in GHRL) likely had a more complex splice-switching effect.

CONCLUSION

Application of our search methods could augment the knowledge value of future cryptic exon reports and aid in generating better hypotheses for genome-wide association studies.

The intronic branch point sequence is under strong evolutionary constraint in the bovine and human genome

The branch point sequence is a cis-acting intronic motif required for mRNA splicing. Despite their functional importance, branch point sequences are not routinely annotated. Here we predict branch point sequences in 179,476 bovine introns and investigate their variability using a catalogue of 29.4 million variants detected in 266 cattle genomes. We localize the bovine branch point within a degenerate heptamer "nnyTrAy". An adenine residue at position 6, that acts as branch point, and a thymine residue at position 4 of the heptamer are more strongly depleted for mutations than coding sequences suggesting extreme purifying selection. We provide evidence that mutations affecting these evolutionarily constrained residues lead to alternative splicing. We confirm evolutionary constraints on branch point sequences using a catalogue of 115 million SNPs established from 3,942 human genomes of the gnomAD database.

Interpretable prioritization of splice variants in diagnostic next-generation sequencing

A critical challenge in genetic diagnostics is the computational assessment of candidate splice variants, specifically the interpretation of nucleotide changes located outside of the highly conserved dinucleotide sequences at the 5' and 3' ends of introns. To address this gap, we developed the Super Quick Information-content Random-forest Learning of Splice variants (SQUIRLS) algorithm. SQUIRLS generates a small set of interpretable features for machine learning by calculating the information-content of wild-type and variant sequences of canonical and cryptic splice sites, assessing changes in candidate splicing regulatory sequences, and incorporating characteristics of the sequence such as exon length, disruptions of the AG exclusion zone, and conservation. We curated a comprehensive collection of disease-associated splice-altering variants at positions outside of the highly conserved AG/GT dinucleotides at the termini of introns. SQUIRLS trains two random-forest classifiers for the donor and for the acceptor and combines their outputs by logistic regression to yield a final score. We show that SQUIRLS transcends previous state-of-the-art accuracy in classifying splice variants as assessed by rank analysis in simulated exomes, and is significantly faster than competing methods. SQUIRLS provides tabular output files for incorporation into diagnostic pipelines for exome and genome analysis, as well as visualizations that contextualize predicted effects of variants on splicing to make it easier to interpret splice variants in diagnostic settings.

Frequent LPA KIV-2 Variants Lower Lipoprotein(a) Concentrations and Protect Against Coronary Artery Disease

Background

Lipoprotein(a) (Lp(a)) concentrations are a major independent risk factor for coronary artery disease (CAD) and are mainly determined by variation in LPA. Up to 70% of the LPA coding sequence is located in the hyper-variable kringle IV type 2 (KIV-2) region. It is hardly accessible by conventional technologies, but may contain functional variants.

Objectives

This study sought to investigate the new, very frequent splicing variant KIV-2 4733G>A on Lp(a) and CAD.

Methods

We genotyped 4733G>A in the GCKD (German Chronic Kidney Disease) study (n = 4,673) by allele-specific polymerase chain reaction, performed minigene assays, identified proxy single nucleotide polymorphisms and used them to characterize its effect on CAD by survival analysis in UK Biobank (n = 440,234). Frequencies in ethnic groups were assessed in the 1000 Genomes Project.

Results

The 4733G>A variant (38.2% carrier frequency) was found in most isoform sizes. It reduces allelic expression without abolishing protein production, lowers Lp(a) by 13.6 mg/dL (95% CI: 12.5-14.7; P < 0.0001) and is the strongest variance-explaining factor after the smaller isoform. Splicing of minigenes was modified. Compound heterozygosity (4.6% of the population) for 4733G>A and 4925G>A, another KIV-2 splicing mutation, reduces Lp(a) by 31.8 mg/dL and most importantly narrows the interquartile range by 9-fold (from 42.1 to 4.6 mg/dL) when compared to the wild type. In UK Biobank 4733G>A alone and compound heterozygosity with 4925G>A reduced HR for CAD by 9% (95% CI: 7%-11%) and 12% (95% CI: 7%-16%) (both P < 0.001). Frequencies in ethnicities differ notably.

Conclusions

Functional variants in the previously inaccessible LPA KIV-2 region cooperate in determining Lp(a) variance and CAD risk. Even a moderate but lifelong genetic Lp(a) reduction translates to a noticeable CAD risk reduction. (J Am Coll Cardiol 2021;78:437–49)

Application of Combined Long Amplicon Sequencing (CoLAS) for Genetic Analysis of Neurofibromatosis Type 1: A Pilot Study

Elaborate analyses of the status of gene mutations in neurofibromatosis type 1 (NF1) are still difficult nowadays due to the large gene sizes, broad mutation spectrum, and the various effects of mutations on mRNA splicing. These problems cannot be solved simply by sequencing the entire coding region using next-generation sequencing (NGS). We recently developed a new strategy, named combined long amplicon sequencing (CoLAS), which is a method for simultaneously analysing the whole genomic DNA region and, also, the full-length cDNA of the disease-causative gene with long-range PCR-based NGS. In this study, CoLAS was specifically arranged for NF1 genetic analysis, then applied to 20 patients (five previously reported and 15 newly recruited patients, including suspicious cases) for optimising the method and to verify its efficacy and benefits. Among new cases, CoLAS detected not only 10 mutations, including three unreported mutations and one mosaic mutation, but also various splicing abnormalities and allelic expression ratios quantitatively. In addition, heterozygous mapping by polymorphisms, including introns, showed copy number monitoring of the entire NF1 gene region was possible in the majority of patients tested. Moreover, it was shown that, when a chromosomal level microdeletion was suspected from heterozygous mapping, it could be detected directly by breakpoint-specific long PCR. In conclusion, CoLAS not simply detect the causative mutation but accurately elucidated the entire structure of the NF1 gene, its mRNA expression, and also the splicing status, which reinforces its high usefulness in the gene analysis of NF1.

Pathogenic noncoding variants in the neurofibromatosis and schwannomatosis predisposition genes

Neurofibromatosis type 1 (NF1), type 2 (NF2), and schwannomatosis are a group of autosomal dominant disorders that predispose to the development of nerve sheath tumors. Pathogenic variants (PVs) that cause NF1 and NF2 are located in the NF1 and NF2 loci, respectively. To date, most variants associated with schwannomatosis have been identified in the SMARCB1 and LZTR1 genes, and a missense variant in the DGCR8 gene was recently reported to predispose to schwannomas. In spite of the high detection rate for PVs in NF1 and NF2 (over 90% of non-mosaic germline variants can be identified by routine genetic screening) underlying PVs for a proportion of clinical cases remain undetected. A higher proportion of non-NF2 schwannomatosis cases have no detected PV, with PVs currently only identified in around 70%-86% of familial cases and 30%-40% of non-NF2 sporadic schwannomatosis cases. A number of variants of uncertain significance have been observed for each disorder, many of them located in noncoding, regulatory, or intergenic regions. Here we summarize noncoding variants in this group of genes and discuss their established or potential role in the pathogenesis of NF1, NF2, and schwannomatosis.

Hybrid Minigene Assay: An Efficient Tool to Characterize mRNA Splicing Profiles of NF1 Variants

Neurofibromatosis type 1 (NF1) is caused by heterozygous loss of function mutations in the NF1 gene. Although patients are diagnosed according to clinical criteria and few genotype-phenotype correlations are known, molecular analysis remains important. NF1 displays allelic heterogeneity, with a high proportion of variants affecting splicing, including deep intronic alleles and changes outside the canonical splice sites, making validation problematic. Next Generation Sequencing (NGS) technologies integrated with multiplex ligation-dependent probe amplification (MLPA) have largely overcome RNA-based techniques but do not detect splicing defects. A rapid minigene-based system was set up to test the effects of NF1 variants on splicing. We investigated 29 intronic and exonic NF1 variants identified in patients during the diagnostic process. The minigene assay showed the coexistence of multiple mechanisms of splicing alterations for seven variants. A leaky effect on splicing was documented in one de novo substitution detected in a sporadic patient with a specific phenotype without neurofibromas. Our splicing assay proved to be a reliable and fast method to validate novel NF1 variants potentially affecting splicing and to detect hypomorphic effects that might have phenotypic consequences, avoiding the requirement of patient's RNA.

WGS and RNA Studies Diagnose Noncoding DMD Variants in Males With High Creatine Kinase

Objective

To describe the diagnostic utility of whole-genome sequencing and RNA studies in boys with suspected dystrophinopathy, for whom multiplex ligation-dependent probe amplification and exomic parallel sequencing failed to yield a genetic diagnosis, and to use remnant normal DMD splicing in 3 families to define critical levels of wild-type dystrophin bridging clinical spectrums of Duchenne to myalgia.

Methods

Exome, genome, and/or muscle RNA sequencing was performed for 7 males with elevated creatine kinase. PCR of muscle-derived complementary DNA (cDNA) studied consequences for DMD premessenger RNA (pre-mRNA) splicing. Quantitative Western blot was used to determine levels of dystrophin, relative to control muscle.

Results

Splice-altering intronic single nucleotide variants or structural rearrangements in DMD were identified in all 7 families. Four individuals, with abnormal splicing causing a premature stop codon and nonsense-mediated decay, expressed remnant levels of normally spliced DMD mRNA. Quantitative Western blot enabled correlation of wild-type dystrophin and clinical severity, with 0%-5% dystrophin conferring a Duchenne phenotype, 10% ± 2% a Becker phenotype, and 15% ± 2% dystrophin associated with myalgia without manifesting weakness.

Conclusions

Whole-genome sequencing relied heavily on RNA studies to identify DMD splice-altering variants. Short-read RNA sequencing was regularly confounded by the effectiveness of nonsense-mediated mRNA decay and low read depth of the giant DMD mRNA. PCR of muscle cDNA provided a simple, yet informative approach. Highly relevant to genetic therapies for dystrophinopathies, our data align strongly with previous studies of mutant dystrophin in Becker muscular dystrophy, with the collective conclusion that a fractional increase in levels of normal dystrophin between 5% and 20% is clinically significant.

AG-exclusion zone revisited: Lessons to learn from 91 intronic NF1 3' splice site mutations outside the canonical AG-dinucleotides

Uncovering frequent motives of action by which variants impair 3′ splice site (3′ss) recognition and selection is essential to improve our understanding of this complex process. Through several mini‐gene experiments, we demonstrate that the pyrimidine (Y) to purine (R) transversion NM_000267.3(NF1):c.1722‐11T>G, although expected to weaken the polypyrimidine tract, causes exon skipping primarily by introducing a novel AG in the AG‐exclusion zone (AGEZ) between the authentic 3′ss AG and the branch point. Evaluation of 90 additional noncanonical intronic NF1 3′ss mutations confirmed that 63% of all mutations and 89% (49/55) of the single‐nucleotide variants upstream of positions ‐3 interrupt the AGEZ. Of these AGEZ‐interrupting mutations, 24/49 lead to exon skipping suggesting that absence of AG in this region is necessary for accurate 3′ss selection already in the initial steps of splicing. The analysis of 91 noncanonical NF1 3′ss mutations also shows that 90% either introduce a novel AG in the AGEZ, cause a Y>R transversion at position ‐3 or remove ≥2 Ys in the AGEZ. We confirm in a validation cohort that these three motives distinguish spliceogenic from splice‐neutral variants with 85% accuracy and, therefore, are generally applicable to select among variants of unknown significance those likely to affect splicing.