massively-parallel-splicing-assay-to-examine-splicing-errors-caused

Massively Parallel Splicing Assay to Examine Splicing Errors Caused by Disease-Related Intronic Variants

Splicing errors represent 10-30% of the pathogenic mutations responsible for rare genetic disorders. RNA splicing ensures proper gene expression by ...

Academia Sinica

Academia Sinica and National Taiwan University

Translation in axons is required for growth cone chemotropic responses to many guidance cues. Although locally synthesized proteins are beginning to be identified, how specific mRNAs are selected for translation remains unclear. Control of poly(A) tail length by cytoplasmic polyadenylation element (CPE) binding protein 1 (CPEB1) is a conserved mechanism for mRNA-specific translational regulation that could be involved in regulating translation in axons.

Cytoplasmic polyadenylation and cytoplasmic polyadenylation element-dependent mRNA regulation are involved in Xenopus retinal axon development.

Full genome analysis of a young girl with deafness, dystonia, central hypomyelination, refractory seizure, and fluctuating liver function impairment revealed a heterozygous, de novo variant in the BCAP31 gene on chromosome Xq28 (NM_001256447.2:c.92G>A), mutations of which caused the X-linked recessive severe neurologic disorder deafness, dystonia, and cerebral hypomyelination. Reverse transcription-polymerase chain reaction of the patient's white blood cells showed the absence of wild-type BCAP31 messenger RNA (mRNA) but the presence of two novel BCAP31 mRNAs. The major alternatively spliced mRNA is due to Exon 2 skipping and the utilization of a new initiation site in Exon 3 that leads to a frameshift and truncated transcript while the minor novel mRNA has a 110 nucleotide insertion to Exon 2. Phasing studies showed that the de novo variant arose in the paternal X chromosome. X chromosome inactivation assay was done and confirmed that the patient's maternal X chromosome was preferentially inactivated, providing evidence that the mutated BCAP31 gene was the one predominantly expressed. According to the American College of Medical Genetics and Genomics guideline, this variant is deemed "pathogenic" (PS2, PS3, PM2, PP3, and PP4) and deleterious. This is the first reported female patient in BCAP31-related syndrome resulted from skewed X-inactivation and a de novo mutation in the active X chromosome.

De novo mutation and skewed X-inactivation in girl with BCAP31-related syndrome.

It is estimated that 10%-30% of disease-associated genetic variants affect splicing. Splicing variants may generate deleteriously altered gene product and are potential therapeutic targets. However, systematic diagnosis or prediction of splicing variants is yet to be established, especially for the near-exon intronic splice region. The major challenge lies in the redundant and ill-defined branch sites and other splicing motifs therein. Here, we carried out unbiased massively parallel splicing assays on 5,307 disease-associated variants that overlapped with branch sites and collected 5,884 variants across the 5' splice region. We found that strong splice sites and exonic features preserve splicing from intronic sequence variation. Whereas the splice-altering mechanism of the 3' intronic variants is complex, that of the 5' is mainly splice-site destruction. Statistical learning combined with these molecular features allows precise prediction of altered splicing from an intronic variant. This statistical model provides the identity and ranking of biological features that determine splicing, which serves as transferable knowledge and out-performs the benchmarking predictive tool. Moreover, we demonstrated that intronic splicing variants may associate with disease risks in the human population. Our study elucidates the mechanism of splicing response of intronic variants, which classify disease-associated splicing variants for the promise of precision medicine.

Mechanism and modeling of human disease-associated near-exon intronic variants that perturb RNA splicing.

The majority of genic transcription is intronic. Introns are removed by splicing as branched lariat RNAs which require rapid recycling. The branch site is recognized during splicing catalysis and later debranched by Dbr1 in the rate-limiting step of lariat turnover. Through generation of a viable DBR1 knockout cell line, we find the predominantly nuclear Dbr1 enzyme to encode the sole debranching activity in human cells. Dbr1 preferentially debranches substrates that contain canonical U2 binding motifs, suggesting that branchsites discovered through sequencing do not necessarily represent those favored by the spliceosome. We find that Dbr1 also exhibits specificity for particular 5' splice site sequences. We identify Dbr1 interactors through co-immunoprecipitation mass spectrometry. We present a mechanistic model for Dbr1 recruitment to the branchpoint through the intron-binding protein AQR. In addition to a 20-fold increase in lariats, Dbr1 depletion increases exon skipping. Using ADAR fusions to timestamp lariats, we demonstrate a defect in spliceosome recycling. In the absence of Dbr1, spliceosomal components remain associated with the lariat for a longer period of time. As splicing is co-transcriptional, slower recycling increases the likelihood that downstream exons will be available for exon skipping.

The debranching enzyme Dbr1 regulates lariat turnover and intron splicing.

Splicing variants are a major class of pathogenic mutations, with their severity equivalent to nonsense mutations. However, redundant and degenerate splicing signals hinder functional assessments of sequence variations within introns, particularly at branch sites. We have established a massively parallel splicing assay to assess the impact on splicing of 11,191 disease-relevant variants. Based on the experimental results, we then applied regression-based methods to identify factors determining splicing decisions and their respective weights.

SpliceAPP: an interactive web server to predict splicing errors arising from human mutations.

Untranslated regions (UTRs) contain crucial regulatory elements for RNA stability, translation and localization, so their integrity is indispensable for gene expression. Approximately 3.7% of genetic variants associated with diseases occur in UTRs, yet a comprehensive understanding of UTR variant functions remains limited due to inefficient experimental and computational assessment methods. To systematically evaluate the effects of UTR variants on RNA stability, we established a massively parallel reporter assay on 6555 UTR variants reported in human disease databases. We examined the RNA degradation patterns mediated by the UTR library in two cell lines, and then applied LASSO regression to model the influential regulators of RNA stability. We found that UA dinucleotides and UA-rich motifs are the most prominent destabilizing element. Gain of UA dinucleotide outlined mutant UTRs with reduced stability. Studies on endogenous transcripts indicate that high UA-dinucleotide ratios in UTRs promote RNA degradation. Conversely, elevated GC content and protein binding on UA dinucleotides protect high-UA RNA from degradation. Further analysis reveals polarized roles of UA-dinucleotide-binding proteins in RNA protection and degradation. Furthermore, the UA-dinucleotide ratio of both UTRs is a common characteristic of genes in innate immune response pathways, implying a coordinated stability regulation through UTRs at the transcriptomic level. We also demonstrate that stability-altering UTRs are associated with changes in biobank-based health indices, underscoring the importance of precise UTR regulation for wellness. Our study highlights the importance of RNA stability regulation through UTR primary sequences, paving the way for further exploration of their implications in gene networks and precision medicine.

Article	Year
Cytoplasmic polyadenylation and cytoplasmic polyadenylation element-dependent mRNA regulation are involved in Xenopus retinal axon development. Neural development\| PubMed ID: 19254368	2009
De novo mutation and skewed X-inactivation in girl with BCAP31-related syndrome. Human mutation\| PubMed ID: 32652807	2020
Mechanism and modeling of human disease-associated near-exon intronic variants that perturb RNA splicing. Nature structural & molecular biology\| PubMed ID: 36303034	2022
The debranching enzyme Dbr1 regulates lariat turnover and intron splicing. Nature communications\| PubMed ID: 38816363	2024
SpliceAPP: an interactive web server to predict splicing errors arising from human mutations. BMC genomics\| PubMed ID: 38877417	2024
Multiplexed assays of human disease-relevant mutations reveal UTR dinucleotide composition as a major determinant of RNA stability. eLife\| PubMed ID: 39964837	2025

Chien-Ling Lin

Biography

JoVE Journal Publications

Other Publications