Deciphering the impact of coding and non-coding SCN1A gene variants on RNA splicing.

Variants that disrupt normal pre-mRNA splicing are increasingly being recognized as a major cause of monogenic disorders. The SCN1A gene, a key epilepsy gene that is linked to various epilepsy phenotypes, is no exception. Approximately 10% of all reported variants in the SCN1A gene are designated as splicing variants, with many located outside of the canonical donor and acceptor splice sites, and most have not been functionally investigated. However, given its restricted expression pattern, functional analysis of splicing variants in the SCN1A gene could not be routinely performed. In this study, we conducted a comprehensive analysis of all reported SCN1A variants and their potential to impact SCN1A splicing and conclude that splicing variants are substantially misannotated and under-represented. We created a splicing reporter system consisting of 18 splicing vectors covering all 26 protein-coding exons with different genomic contexts and several promoters of varying strengths in order to reproduce the wild-type splicing pattern of the SCN1A gene, revealing cis-regulatory elements essential for proper recognition of SCN1A exons. Functional analysis of 95 SCN1A variants was carried out, including all 68 intronic variants reported in the literature, located outside of the splice sites canonical dinucleotides; 21 exonic variants of different classes (synonymous, missense, nonsense and in-frame deletion) and six variants observed in patients with epilepsy. Interestingly, almost 20% of tested intronic variants had no influence on SCN1A splicing, despite being reported as causative in the literature. Moreover, we confirmed that the majority of predicted exonic variants affect splicing unravelling their true molecular mechanism. We used functional data to perform genotype-phenotype correlation, revealing distinct distribution patterns for missense and splice-affecting 'missense' variants and observed no difference in the phenotype severity of variants leading to in-frame and out-of-frame isoforms, indicating that the Nav1.1 protein is highly intolerant to structural variations. Our work demonstrates the importance of functional analysis in proper variant annotation and provides a tool for high-throughput delineation of splice-affecting variants in SCN1A in a whole-gene manner.

Clinical and Functional Characteristics of the E92K CFTR Gene Variant in the Russian and Turkish Population of People with Cystic Fibrosis.

The pathogenic variant E92K (c.274G > A) of the CFTR gene is rare in America and Europe, but it is common for people with cystic fibrosis from Russia and Turkey. We studied the effect of the E92K genetic variant on the CFTR function. The function of the CFTR channel was studied using the intestinal current measurements (ICM) method. The effects of CFTR modulators on the restoration of the CFTR function were studied in the model of intestinal organoids. To assess the effect of E92K on pre-mRNA splicing, the RT-PCR products obtained from patients' intestinal organoid cultures were analyzed. Patients with the genetic variant E92K are characterized by an older age of diagnosis compared to homozygotes F508del and a high frequency of pancreatic sufficiency. The results of the sweat test and the ICM method showed partial preservation of the function of the CFTR channel. Functional analysis of CFTR gene expression revealed a weak effect of the E92K variant on mRNA-CFTR splicing. Lumacaftor (VX-809) has been shown to restore CFTR function in an intestinal organoid model, which allows us to consider the E92K variant as a promising target for therapy with CFTR correctors.

Splice-site variant in the RPS7 5'-UTR leads to a decrease in the mRNA level and development of Diamond-Blackfan anemia.

Diamond-Blackfan anemia (DBA) is an inherited bone marrow failure syndrome characterized by erythroid aplasia. Pathogenic variants in ribosomal protein (RP) genes, GATA1, TSR2, and EPO, are considered to be the etiology of DBA. Variants in 5'-untranslated regions (UTRs) of these genes are poorly studied and can complicate the variant interpretation. We investigated the functional consequences NM_001011.4:c.-19 + 1G > T variant in the donor splice-site of the RPS7 5'-UTR. This variant was found in a family where two sons with DBA were carriers. Father, who also had this variant, developed myelodysplastic syndrome, which caused his death. Search for candidate causal variants and copy number variations in DBA-associated genes left RPS7 variant as the best candidate. Trio whole exome sequencing analysis revealed no pathogenic variants in other genes. Functional analysis using luciferase expression system revealed that this variant leads to disruption of splicing. Also, a decrease in the levels of mRNA and protein expression was detected. In conclusion, the established consequences of 5'-UTR splice-site variant c.-19 + 1G > T in the RPS7 gene provide evidence that it is likely pathogenic.

Health Characteristics of Patients with Cystic Fibrosis whose Genotype Includes a Variant of the Nucleotide Sequence c.3140-16T>A and Functional Analysis of this Variant.

Cystic fibrosis (CF) is the most common monogenic autosomal recessive disease, associated with pathogenic variants in the CFTR gene. The splicing variant c.3140-16T>A (3272-16T>A) has been described previously and, according to the Russian CF Patients Registry, occurs with a frequency of 0.34%. The phenotypic features of CF patients with the c.3140-16T>A variant were compared with those of patients with the genotype F508del/F508del. Patients with the allele c.3140-16T>A had higher average age and age at diagnosis, and the allele was present in a greater proportion of adults. Patients carrying the c.3140-16T>A allele were characterised by better physical development indicators, both in adults and in children, had preserved pancreatic function, as well as the absence of a number of complications, and required pancreatic enzyme replacement therapy less often than patients with the F508del/F508del genotype. Sweat test values also were lower in patients with the c.3140-16T>A genotype. According to the results of clinical and laboratory studies, the phenotype of patients with the genetic variant c.3140-16T>A can be considered "mild". Functional CFTR protein activity in the presence of c.3140-16T>A was evaluated using intestinal current measurements (ICM) and the forskolin-induced swelling assay on organoids obtained from patients' rectal biopsies. c.3140-16T>A had high residual CFTR channel activity and was amenable to effective pharmacological correction with thea VX-770 potentiator. To evaluate the effect of the variant on CFTR pre-mRNA splicing we performed a minigene assay, as well as RT-PCR analysis of RNA isolated from the nasal epithelium and rectal biopsy of patients. We showed that the c.3140-16T>A variant creates a novel acceptor AG dinucleotide within CFTR intron 19, resulting in a 14-nucleotide extension of exon 20. This frameshift produces a premature termination codon and triggers mRNA degradation by the nonsense-mediated decay (NMD) mechanism. Moreover, we observed that the c.3140-16T>A allele could produce a residual amount of normally spliced transcript, thus explaining the patient's mild phenotype.