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

BACKGROUND: 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. RESULTS: Our statistical modeling is highly sensitive, accurately annotating the splicing defects of near-exon intronic variants, outperforming state-of-the-art predictive tools. We have incorporated the algorithm and branchpoint information into a web-based tool, SpliceAPP, to provide an interactive application. This user-friendly website allows users to upload any genetic variants with genome coordinates (e.g., chr15 74,687,208 A G), and the tool will output predictions for splicing error scores and evaluate the impact on nearby splice sites. Additionally, users can query branch site information within the region of interest. CONCLUSIONS: In summary, SpliceAPP represents a pioneering approach to screening pathogenic intronic variants, contributing to the development of precision medicine. It also facilitates the annotation of splicing motifs. SpliceAPP is freely accessible using the link https://bc.imb.sinica.edu.tw/SpliceAPP . Source code can be downloaded at https://github.com/hsinnan75/SpliceAPP .

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

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.

Identification of functional single nucleotide polymorphisms in the branchpoint site.

BACKGROUND: The human genome contains millions of single nucleotide polymorphisms (SNPs); many of these SNPs are intronic and have unknown functional significance. SNPs occurring within intron branchpoint sites, especially at the adenine (A), would presumably affect splicing; however, this has not been systematically studied. We employed a splicing prediction tool to identify human intron branchpoint sites and screened dbSNP for identifying SNPs located in the predicted sites to generate a genome-wide branchpoint site SNP database. RESULTS: We identified 600 SNPs located within branchpoint sites; among which, 216 showed a change in A. After scoring the SNPs by counting the As in the ± 10 nucleotide region, only four SNPs were identified without additional As (rs13296170, rs12769205, rs75434223, and rs67785924). Using minigene constructs, we examined the effects of these SNPs on splicing. The three SNPs (rs13296170, rs12769205, and rs75434223) with nucleotide substitution at the A position resulted in abnormal splicing (exon skipping and/or intron inclusion). However, rs67785924, a 5-bp deletion that abolished the branchpoint A nucleotide, exhibited normal RNA splicing pattern, presumably using two of the downstream As as alternative branchpoints. The influence of additional As on splicing was further confirmed by studying rs2733532, which contains three additional As in the ± 10 nucleotide region. CONCLUSIONS: We generated a high-confidence genome-wide branchpoint site SNP database, experimentally verified the importance of A in the branchpoint, and suggested that other nearby As can protect branchpoint A substitution from abnormal splicing.

Genetic epidemiological study doesn't support GLA IVS4+919G>A variant is a significant mutation in Fabry disease.

BACKGROUND: The GLA IVS4+919G>A which is linked to late-onset Fabry disease shows high frequency in Taiwan. METHODS: To determine whether IVS4+919G>A is a frequent cause of heart disease, we genotyped it in normal controls and other disease cohorts (type 2 diabetes, heart failure, ventricular tachycardia, atrial fibrillation and coronary artery disease). Normal controls and diabetes patients carrying the variant were evaluated for their cardiac condition. Minigene constructs were used to study GLA splicing patterns in different cell lines. RESULTS: GLA IVS4+919A was found in 4/1634 males (0.245%) and 2/1634 females (0.123%) in normal controls and in 4/2133 males (0.188%) and 7/1816 females (0.385%) in the type 2 diabetes cohort. Of all the 17 IVS4+919A carriers in these two groups, only two males reported heart-related disease (myocardial infarction and hypertensive heart disease). Furthermore, in the heart disease cohort (n=649), only one male carried the variant. Minigene constructs showed that the AGS (stomach) cell line showed a distinct GLA splicing pattern. CONCLUSION: Most subjects carrying GLA IVS4+919A did not show abnormal cardiac phenotypes. The near-absence of GLA IVS4+919A in heart disease cohort suggested that this variant is not a frequent cause of overt heart diseases in Taiwan and that the genotype-phenotype correlation and natural course of the disease need further investigation. We also showed that the GLA IVS4+919G>A nucleotide change did influence alternative splicing in a tissue-specific manner. SYNOPSIS: The GLA IVS4+919G>A variant is not a frequent cause of overt heart disease in Taiwan.

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

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.

Genome-wide association study of treatment refractory schizophrenia in Han Chinese.

We report the first genome-wide association study of a joint analysis using 795 Han Chinese individuals with treatment-refractory schizophrenia (TRS) and 806 controls. Three loci showed suggestive significant association with TRS were identified. These loci include: rs10218843 (P = 3.04 × 10(-7)) and rs11265461 (P = 1.94 × 10(-7)) are adjacent to signaling lymphocytic activation molecule family member 1 (SLAMF1); rs4699030 (P = 1.94 × 10(-6)) and rs230529 (P = 1.74 × 10(-7)) are located in the gene nuclear factor of kappa light polypeptide gene enhancer in B-cells 1 (NFKB1); and rs13049286 (P = 3.05 × 10(-5)) and rs3827219 (P = 1.66 × 10(-5)) fall in receptor-interacting serine/threonine-protein kinase 4 (RIPK4). One isolated single nucleotide polymorphism (SNP), rs739617 (P = 3.87 × 10(-5)) was also identified to be associated with TRS. The -94delATTG allele (rs28362691) located in the promoter region of NFKB1 was identified by resequencing and was found to associate with TRS (P = 4.85 × 10(-6)). The promoter assay demonstrated that the -94delATTG allele had a significant lower promoter activity than the -94insATTG allele in the SH-SY5Y cells. This study suggests that rs28362691 in NFKB1 might be involved in the development of TRS.

Novel DNA sequence variations of cytochrome P450 genes in the Han Chinese population.

AIM: Cytochrome P450 is an important monooxygenase responsible for the metabolism of a large variety of structurally diverse compounds. The aim of this study was to systematically investigate the DNA sequence variations in 14 cytochrome P450 genes relevant to drug metabolism in the Han Chinese population. MATERIALS & METHODS: We sequenced these 14 genes in 23 subjects and determined the allele frequencies. RESULTS: We identified a total of 312 genetic variants, which included 80 (25.6%) novel variants. These novel variants included 67 noncoding variants, four synonymous and nine nonsynonymous variants. Among these variants, an 11-nucleotide insertion in the 3-flanking site of intron 6 in CYP2E1 (c.968-77_-76insGATGGGTGGAT) had the highest allele frequency of 0.565. A total of 16 of these novel variants were predicted to have potential functional consequences; however, among them only c.-1299T>A in CYP2C18 and c.-498C>A in CYP2D6 reached a frequency of 4.9%. CONCLUSIONS: This study establishes a genetic database of cytochrome P450 genes in the Han Chinese population and suggests further genetic diversity throughout this important gene family.