Functional characterization of variants found in Japanese patients with hereditary hemorrhagic telangiectasia

Hereditary hemorrhagic telangiectasia (HHT) is an autosomal dominant form of vascular dysplasia. Genetic diagnosis is made by identifying loss-of-function variants in genes, such as ENG and ACVRL1. However, the causal mechanisms of various variants of unknown significance remains unclear. In this study, we analyzed 12 Japanese patients from 11 families who were clinically diagnosed with HHT. Sequencing analysis identified 11 distinct variants in ACVRL1 and ENG. Three of the 11 were truncating variants, leading to a definitive diagnosis, whereas the remaining eight were splice-site and missense variants that required functional analyses. In silico splicing analyses demonstrated that three variants, c.526-3C > G and c.598C > G in ACVRL1, and c.690-1G > A in ENG, caused aberrant splicing, as confirmed by a minigene assay. The five remaining missense variants were p.Arg67Gln, p.Ile256Asn, p.Leu285Pro, and p.Pro424Leu in ACVRL and p.Pro165His in ENG. Nanoluciferase-based bioluminescence analyses demonstrated that these ACVRL1 variants impaired cell membrane trafficking, resulting in the loss of bone morphogenetic protein 9 (BMP9) signal transduction. In contrast, the ENG mutation impaired BMP9 signaling despite normal cell membrane expression. The updated functional analysis methods performed in this study will facilitate effective genetic testing and appropriate medical care for patients with HHT.

Identification of a novel TSC1 gene variant in a patient with atypical vitiligo-like skin lesions: Unveiling the hidden tuberous sclerosis complex

BACKGROUND

Tuberous sclerosis complex (TSC), an autosomal-dominant disorder, is characterized by hamartomas affecting multiple organ systems. The underlying etiology of TSC is the pathogenic variations of the TSC1 or TSC2 genes. The phenotype variability of TSC could lead to missed diagnosis; therefore, the latest molecular diagnostic criteria for identifying a heterozygous pathogenic variant in either the TSC1 or TSC2 gene filled this gap. Furthermore, the pathogenicity of numerous variants remains unverified, potentially leading to misinterpretations of their functional consequences.

METHODS

In this study, a single patient presenting with atypical vitiligo-like skin lesions suspected to have TSC was enrolled. Targeted next-generation sequencing and Sanger sequencing were employed to identify a pathogenic variant. Additionally, a minigene splicing assay was conducted to assess the impact of TSC1 c.1030-2A>T, located in intron 10, on RNA splicing.

RESULTS

A novel TSC1: c.1030-2A>T heterozygosis variant was identified in intron 10. In vitro minigene assay revealed that the c.1030-2A>T variant caused exon 11 skipping, resulting in a frameshift in the absence of 112 base pairs of mature messenger RNA and premature termination after 174 base pairs (p.Ala344Glnfs*59).

CONCLUSION

The detection of this novel pathogenic TSC1 variant in the patient with atypical vitiligo-like skin lesions enrolled in our study ultimately resulted in the diagnosis of TSC. As a result, our study contributes to expanding the mutational spectrum of the TSC1 gene and refining the genotype-phenotype map of TSC.

Cohort analysis of novel SPAST variants in SPG4 patients and implementation of in vitro and in vivo studies to identify the pathogenic mechanism caused by splicing mutations

Introduction

Pure hereditary spastic paraplegia (SPG) type 4 (SPG4) is caused by mutations of SPAST gene. This study aimed to analyze SPAST variants in SPG4 patients to highlight the occurrence of splicing mutations and combine functional studies to assess the relevance of these variants in the molecular mechanisms of the disease.

Methods

We performed an NGS panel in 105 patients, in silico analysis for splicing mutations, and in vitro minigene assay.

Results and discussion

The NGS panel was applied to screen 105 patients carrying a clinical phenotype corresponding to upper motor neuron syndrome (UMNS), selectively affecting motor control of lower limbs. Pathogenic mutations in SPAST were identified in 12 patients (11.42%), 5 missense, 3 frameshift, and 4 splicing variants. Then, we focused on the patients carrying splicing variants using a combined approach of in silico and in vitro analysis through minigene assay and RNA, if available. For two splicing variants (i.e., c.1245+1G>A and c.1414-2A>T), functional assays confirm the types of molecular alterations suggested by the in silico analysis (loss of exon 9 and exon 12). In contrast, the splicing variant c.1005-1delG differed from what was predicted (skipping exon 7), and the functional study indicates the loss of frame and formation of a premature stop codon. The present study evidenced the high splice variants in SPG4 patients and indicated the relevance of functional assays added to in silico analysis to decipher the pathogenic mechanism.

Pathogenicity analysis and splicing rescue of a classical splice site variant (c.1343+1G>T) of CNOT1 gene associated with neurodevelopmental disorders

Mutations in the CNOT1 gene lead to an incurable rare neurological disorder mainly manifested as a clinical spectrum of intellectual disability, developmental delay, seizures, and behavioral problems. In this study, we investigated a classical splice site variant of CNOT1 (c.1343+1G>T) associated with neurodevelopmental disorders, which was a master regulator, orchestrating gene expression, RNA deadenylation, and protein ubiquitination. To link CNOT1 dysfunction with the neurodevelopmental phenotype observed in a patient, in vitro minigene assay was used to verify the effect of CNOT1 gene splice site variant c.1343+1G>T on mRNA splicing. We also explored the impact of transient transfection introducing modified U1 snRNA on correcting the splicing variant. Through minigene expression in mammalian cells, we demonstrated that the variant induced complete exon 12 skipping, which explained the patient's clinical condition and provided additional genetic diagnosis evidence for the clinical significance of the variant. Moreover, we confirmed that the aberrant splice pattern could be partially corrected by the modified U1 snRNA at the mRNA level, which provided strong evidence for the therapeutic potential of modified U1 snRNA in neutralizing the hazardous effect of incorrect splicing patterns.

Evaluation of Pathogenicity and Causativity of Variants in the MPZ and SH3TC2 Genes in a Family Case of Hereditary Peripheral Neuropathy

The implementation of NGS methods into clinical practice allowed researchers effectively to establish the molecular cause of a disorder in cases of a genetically heterogeneous pathology. In cases of several potentially causative variants, we need additional analysis that can help in choosing a proper causative variant. In the current study, we described a family case of hereditary motor and sensory neuropathy (HMSN) type 1 (Charcot-Marie-Tooth disease). DNA analysis revealed two variants in the SH3TC2 gene (c.279G>A and c.1177+5G>A), as well as a previously described variant c.449-9C>T in the MPZ gene, in a heterozygous state. This family segregation study was incomplete because of the proband's father's unavailability. To evaluate the variants' pathogenicity, minigene splicing assay was carried out. This study showed no effect of the MPZ variant on splicing, but the c.1177+5G>A variant in the SH3TC2 gene leads to the retention of 122 nucleotides from intron 10 in the RNA sequence, causing a frameshift and an occurrence of a premature stop codon (NP_078853.2:p.Ala393GlyfsTer2).

Biallelic Variants in TULP1 Are Associated with Heterogeneous Phenotypes of Retinal Dystrophy

Biallelic pathogenic variants in TULP1 are mostly associated with severe rod-driven inherited retinal degeneration. In this study, we analyzed clinical heterogeneity in 17 patients and characterized the underlying biallelic variants in TULP1. All patients underwent thorough ophthalmological examinations. Minigene assays and structural analyses were performed to assess the consequences of splice variants and missense variants. Three patients were diagnosed with Leber congenital amaurosis, nine with early onset retinitis pigmentosa, two with retinitis pigmentosa with an onset in adulthood, one with cone dystrophy, and two with cone-rod dystrophy. Seventeen different alleles were identified, namely eight missense variants, six nonsense variants, one in-frame deletion variant, and two splice site variants. For the latter two, minigene assays revealed aberrant transcripts containing frameshifts and premature termination codons. Structural analysis and molecular modeling suggested different degrees of structural destabilization for the missense variants. In conclusion, we report the largest cohort of patients with TULP1-associated IRD published to date. Most of the patients exhibited rod-driven disease, yet a fraction of the patients exhibited cone-driven disease. Our data support the hypothesis that TULP1 variants do not fold properly and thus trigger unfolded protein response, resulting in photoreceptor death.

Phenotype-Based Genetic Analysis Reveals Missing Heritability of KIF11-Related Retinopathy: Clinical and Genetic Findings

The purpose of this study was to detect the missing heritability of patients with KIF11-related retinopathy and to describe their clinical and genetic characteristics. We enrolled 10 individuals from 7 unrelated families harboring a pathogenic monoallelic variant in KIF11. All subjects underwent ophthalmic assessment and extraocular phenotype evaluations, as well as comprehensive molecular genetic analyses using next-generation sequencing. Minigene assays were performed to observe the effects of one novel deep intron variant (DIV) and one novel synonymous variant on pre-mRNA splicing. We detected 6 novel different disease-causing variants of KIF11 in the seven pedigrees. Co-segregation analysis and ultra-deep sequencing results indicated that 5 variants arose de novo in 5 families (71%). Functional validation revealed that the synonymous variant leads to an exon skip, while the DIV causes a pseudoexon (PE) inclusion. The patients presented with high variations in their phenotype, and two families exhibited incomplete penetrance. Ocular manifestations and characteristic facial features were observed in all patients, as well as microcephaly in seven patients, intellectual disability in five patients, and lymphedema in one patient. The key retinal features for KIF11-related retinopathy were retinal folds, tractional retinal detachment, and chorioretinal dysplasia. All seven probands had more severe visual detects than other affected family members. Our findings widen the genetic spectrum of KIF11 variants. DIV explained rare unresolved cases with KIF11-related retinopathy. The patients displayed a variable phenotype expressivity and incomplete penetrance, indicating the importance of genetic analysis for patients with KIF11-related retinopathy.

Minigene splicing assays reveal new insights into exonic variants of the SLC12A3 gene in Gitelman syndrome

BACKGROUND

Gitelman syndrome (GS) is a type of salt-losing tubular disease, most of which is caused by SLC12A3 gene variants, and missense variants account for the majority. Recently, the phenomenon of exon skipping, in which variants disrupt normal pre-mRNA splicing, has been related to a variety of diseases. Therefore, we hypothesize that a certain proportion of SLC12A3 variants can result in disease via interfering with the normal splicing process.

METHODS

We analyzed 342 previously presumed SLC12A3 missense variants using bioinformatics programs and identified candidate variants that may alter the splicing of pre-mRNA through minigene assays.

RESULTS

Our study revealed that, among ten candidate variants, six variants (c.602G>A, c.602G>T, c.1667C>T, c.1925G>A, c.2548G>C, and c.2549G>C) led to complete or incomplete exon skipping by affecting exonic splicing regulatory elements and/or disturbing canonical splice sites.

CONCLUSION

It is worth mentioning that this is the largest study on pre-mRNA splicing of SLC12A3 exonic variants. In addition, our study emphasizes the importance of detecting splicing function at the mRNA level in GS and indicates that minigene analysis is a valuable tool for splicing functional assays of variants in vitro.

Molecular dynamics and minigene assay of new splicing variant c.4298-20T>A of COL4A5 gene that cause Alport syndrome

Introduction: Alport syndrome (AS; OMIM#308940) is a progressive hereditary kidney disease characterized by hearing loss and ocular abnormalities. According to the mode of inheritance, AS has three subtypes: X-linked (XL; OMIM#301050), autosomal recessive (AR; OMIM#203780), and autosomal dominant (AD; OMIM#104200). XLAS is caused by a pathogenic variant in COL4A5 (OMIM*303630) gene encoding type IV collagen (Col-IV) α5 chain, while ADAS and ARAS are consequences of a variant in COL4A3 (OMIM*120070) and COL4A4 (OMIM*120131) genes that encode Col-IV α3 and α4 chains, respectively. Usually, diagnosis of AS requires hereditary or pathological examinations. Splicing variants are hard to be determined as pathogenic or non-pathogenic based on the results of gene sequencing. Methods: This study focused on a splicing variant in COL4A5 gene, termed NM_000495.5: c.4298-20T>A, and to analyzed its authenticity and damaged α5 chain. In vitro minigene splicing assay was applied to investigate the effect of splicing variant, c.4298-20T>A, on COL4A5 mRNA synthesis. Molecular dynamics method was used to predict the capability of the responsible α5(IV) to form a triple helix. Results: The intron 46 of COL4A5 mRNA retained 18 bp, resulting in insertion of six amino acids behind the amino acid at position 1,433 of α5(IV). The predicted protein effect of this variant: p. (Pro1432_Gly1433insAspTyrPheValGluIle). As a consequence, the stability of α5(IV) secondary structure was impaired, probably leading to the unusual configuration of α345(IV). Discussion: Normally, splicing variant in COL4A5 gene can lead to phenotypes of XLAS, and the effect is associated with the extent of splicing. The patient reported here carried a c.4298-20T>A splicing variant in COL4A5 gene, and AS was highly suspected based on the pathology results. However, the patient did not manifest any ocular or ear abnormalities. We therefore present the c.4298-20T>A splicing variant in COL4A5 gene as likely-pathogenic splicing variant that leads to XLAS with mild phenotypes.

Genetic etiological analysis of auditory neuropathy spectrum disorder by next-generation sequencing

Objective

Auditory neuropathy spectrum disease (ANSD) is caused by both environmental and genetic causes and is defined by a failure in peripheral auditory neural transmission but normal outer hair cells function. To date, 13 genes identified as potentially causing ANSD have been documented. To study the etiology of ANSD, we collected 9 probands with ANSD diagnosed in the clinic and performed targeted next-generation sequencing.

Methods

Nine probands have been identified as ANSD based on the results of the ABR tests and DPOAE/CMs. Genomic DNA extracted from their peripheral blood was examined by next-generation sequencing (NGS) for a gene panel to identify any potential causal variations. For candidate pathogenic genes, we performed co-segregation among all family members of the pedigrees. Subsequently, using a mini-gene assay, we examined the function of a novel splice site mutant of OTOF.

Results

We analyzed nine cases of patients with ANSD with normal CMs/DPOAE and abnormal ABR, discovered three novel mutants of the OTOF gene that are known to cause ANSD, and six cases of other gene mutations including TBC1D24, LARS2, TIMM8A, MITF, and WFS1.

Conclusion

Our results extend the mutation spectrum of the OTOF gene and indicate that the genetic etiology of ANSD may be related to gene mutations of TBC1D24, LARS2, TIMM8A, MITF, and WFS1.

Characterization of a possible founder synonymous variant in TECTA in multiple individuals with autosomal recessive hearing loss

Synonymous variants have been shown to alter the correct splicing of pre-mRNAs and generate disease-causing transcripts. These variants are not an uncommon etiology of genetic disease; however, they are frequently overlooked during genetic testing in the absence of functional and clinical data. Here, we describe the occurrence of a synonymous variant [NM_005422.4 (TECTA):c.327C>T, p.(Gly109=)] in seven individuals with hearing loss from six unrelated families. The variant is not located near exonic/intronic boundaries but is predicted to impact splicing by activating a cryptic splicing donor site in exon 4 of TECTA. In vitro minigene assays show that the variant disrupts the reading frame of the canonical transcript, which is predicted to cause a premature termination codon 48 amino acids downstream of the variant, leading to nonsense-mediated decay. The variant is present in population databases, predominantly in Latinos of African ancestry, but is rare in other ethnic groups. Our findings suggest that this synonymous variant is likely pathogenic for TECTA-associated autosomal recessive hearing loss and seems to have arisen as a founder variant in this specific Latino subpopulation. This study demonstrates that synonymous variants need careful splicing assessment and support from additional testing methodologies to determine their clinical impact.

Detecting pathogenic deep intronic variants in Gitelman syndrome

Gitelman syndrome (GS) is a rare, autosomal recessive, salt-losing tubulopathy caused by loss of function in the SLC12A3 gene (NM_000339.2), which encodes the natrium chloride cotransporter. The detection of homozygous or compound heterozygous SLC12A3 variants is expected in GS, but 18%-40% of patients with clinical GS carry only one mutant allele. Previous reports identified some pathogenic deep intronic variants in SLC12A3. Here, we report the screening of SLC12A3 deep intronic variants in 13 patients with suspected GS carrying one mutated SLC12A3 allele. Variant screening used the HaloPlex Target Enrichment System Kit capturing whole introns and the promotor region of SLC12A3, followed by SureCall variant analysis. Rare intronic variants (<1% frequency) were identified, and pathogenicity evaluated by the minigene system. Deep intronic variant screening detected seven rare SLC12A3 variants from six patients. Only one variant showed pathogenicity in the minigene system (c.602-16G>A, intron 4) through activation of a cryptic acceptor site. No variants were detected in the promotor region. Deep intronic screening identified only one pathogenic variant in patients with suspected GS carrying monoallelic SLC12A3 variants. Our results suggest that deep intronic variants partially explain the cause of monoallelic variants in patients with GS.

Effects of 14 F9 synonymous codon variants on hemophilia B expression: Alteration of splicing along with protein expression

There is growing evidence that synonymous codon variants (SCVs) can cause disease through the disruption of different processes of protein production. The aim of the study is to investigate whether the 14 SCVs reported in the F9 variant database were the pathogenic causes of hemophilia B. The impacts of SCVs on splicing and protein expression were detected using a combination of in silico prediction, in vitro minigene splicing assay and cell expression detection. The splicing transcripts were identified and quantified by co-amplification fluorescent PCR. The mechanism of splicing was verified by a modified pU1snRNA and pU7snRNA approach. Aberrant splicing patterns were found in eight SCVs. Five of the 8 SCVs produced almost all aberrant splicing isoforms, which were expected to truncate protein, three of them presented a partial defect on both splicing and protein secretion, the overall effects were consistent with the residual Factor IX activity of the affected cases. Neither the pre-messenger RNA (mRNA) splicing process nor the protein function was impaired in the rest six SCVs. In conclusion, our study firstly revealed the pathogenic mechanism of the 14 F9 SCVs and highlighted the importance of performing mRNA splicing analysis and protein expression studies of SCVs in inherited disorders.

Novel OPN1LW/OPN1MW Exon 3 Haplotype-Associated Splicing Defect in Patients with X-Linked Cone Dysfunction

Certain combinations of common variants in exon 3 of OPN1LW and OPN1MW, the genes encoding the apo-protein of the long- and middle-wavelength sensitive cone photoreceptor visual pigments in humans, induce splicing defects and have been associated with dyschromatopsia and cone dysfunction syndromes. Here we report the identification of a novel exon 3 haplotype, G-C-G-A-T-T-G-G (referring to nucleotide variants at cDNA positions c.453, c.457, c.465, c.511, c.513, c.521, c.532, and c.538) deduced to encode a pigment with the amino acid residues L-I-V-V-A at positions p.153, p.171, p.174, p.178, and p.180, in OPN1LW or OPN1MW or both in a series of seven patients from four families with cone dysfunction. Applying minigene assays for all observed exon 3 haplotypes in the patients, we demonstrated that the novel exon 3 haplotype L-I-V-V-A induces a strong but incomplete splicing defect with 3-5% of residual correctly spliced transcripts. Minigene splicing outcomes were similar in HEK293 cells and the human retinoblastoma cell line WERI-Rb1, the latter retaining a cone photoreceptor expression profile including endogenous OPN1LW and OPN1MW gene expression. Patients carrying the novel L-I-V-V-A haplotype presented with a mild form of Blue Cone Monochromacy or Bornholm Eye Disease-like phenotype with reduced visual acuity, reduced cone electroretinography responses, red-green color vision defects, and frequently with severe myopia.

Case report: Altered pre-mRNA splicing caused by intronic variant c.1499 + 1G > A in the SLC4A4 gene

Proximal renal tubular acidosis (pRTA) with ocular abnormalities is an autosomal recessive disease caused by variants in the Solute Carrier Family 4 Member 4 (SLC4A4) gene. Patients present with metabolic acidosis and low plasma bicarbonate concentration (3∼17 mmol/L). In addition, they are often accompanied by ocular abnormalities, intellectual disability, and growth retardation. The patient underwent whole exome sequencing (WES) and bioinformatics analysis of variant pathogenicity in this study. Then, a minigene assay was conducted to analyze the splicing site variant further. Compound heterozygous variants in the SLC4A4 gene (NM_003759.3), c.145C > T (p.Arg49*) and c.1499 + 1G > A, were detected by WES. The minigene assay showed an mRNA splicing aberration caused by the c.1499 + 1G > A variant. Compared with the wild type, the mutant type caused 4-base insertion between exons 10 and 11 of SLC4A4 after expression in HEK293 cells. In conclusion, the c.1499 + 1G > A variant in the SLC4A4 gene may be one of the genetic causes in the patient. Moreover, our study provides the foundation for future gene therapy of such pathogenic variants.

Systematic Review of Genotype-Phenotype Correlations in Frasier Syndrome

Introduction

Frasier syndrome (FS) is a rare inherited kidney disease caused by intron 9 splicing variants of WT1. For wild-type WT1, 2 active splice donor sites in intron 9 cause a mixture of 2 essential transcripts (with or without lysine-threonine-serine [+/KTS or −KTS]), and imbalance of the +KTS/−KTS ratio results in the development of FS. To date, 6 causative intron 9 variants have been identified; however, detailed transcript analysis has not yet been conducted and the genotype-phenotype correlation also remains to be elucidated.

Methods

We conducted an in vitro minigene splicing assay for 6 reported causative variants and in vivo RNA sequencing to determine the +KTS/−KTS ratio using patients’ samples. We also performed a systematic review of reported FS cases with a description of the renal phenotype.

Results

The in vitro assay revealed that although all mutant alleles produced −KTS transcripts only, the wild-type allele produced both +KTS and −KTS transcripts at a 1:1 ratio. In vivo RNA sequencing showed that patients’ samples with all heterozygous variants produced similar ratios of +KTS to −KTS (1:3.2−1:3.5) and wild-type kidney showed almost a 1:1 ratio (1:0.85). A systematic review of 126 cases clarified that the median age of developing ESKD was 16 years in all FS patients, and there were no statistically significant differences between the genotypes or sex chromosome karyotypes in terms of the renal survival period.

Conclusion

Our study suggested no differences in splicing pattern or renal survival period among reported intron 9 variants causative of FS.

A Novel Splicing Variant of COL2A1 in a Fetus with Achondrogenesis Type II: Interpretation of Pathogenicity of In-Frame Deletions

Achondrogenesis type II (ACG2) is a lethal skeletal dysplasia caused by dominant pathogenic variants in COL2A1. Most of the variants found in patients with ACG2 affect the glycine residue included in the Gly-X-Y tripeptide repeat that characterizes the type II collagen helix. In this study, we reported a case of a novel splicing variant of COL2A1 in a fetus with ACG2. An NGS analysis of fetal DNA revealed a heterozygous variant c.1267-2_1269del located in intron 20/exon 21. The variant occurred de novo since it was not detected in DNA from the blood samples of parents. We generated an appropriate minigene construct to study the effect of the variant detected. The minigene expression resulted in the synthesis of a COL2A1 messenger RNA lacking exon 21, which generated a predicted in-frame deleted protein. Usually, in-frame deletion variants of COL2A1 cause a phenotype such as Kniest dysplasia, which is milder than ACG2. Therefore, we propose that the size and position of an in-frame deletion in COL2A1 may be relevant in determining the phenotype of skeletal dysplasia.

Usefulness of functional splicing analysis to confirm precise disease pathogenesis in Diamond-Blackfan anemia caused by intronic variants in RPS19

Diamond-Blackfan anemia (DBA) is mainly caused by pathogenic variants in ribosomal proteins and 22 responsible genes have been identified to date. The most common causative gene of DBA is RPS19 [NM_001022.4]. Nearly 180 RPS19 variants have been reported, including three deep intronic variants outside the splicing consensus sequence (c.72-92A > G, c.356 + 18G > C, and c.411 + 6G > C). We also identified one case with a c.412-3C > G intronic variant. Without conducting transcript analysis, the pathogenicity of these variants is unknown. However, it is difficult to assess transcripts because of their fragility. In such cases, in vitro functional splicing assays can be used to assess pathogenicity. Here, we report functional splicing analysis results of four RPS19 deep intronic variants identified in our case and in previously reported cases. One splicing consensus variant (c.411 + 1G > A) was also examined as a positive control. Aberrant splicing with a 2-bp insertion between exons 5 and 6 was identified in the patient samples and minigene assay results also identified exon 6 skipping in our case. The exon 6 skipping transcript was confirmed by further evaluation using quantitative RT-PCR. Additionally, minigene assay analysis of three reported deep intronic variants revealed that none of them showed aberrant splicing and that these variants were not considered to be pathogenic. In conclusion, the minigene assay is a useful method for functional splicing analysis of inherited disease.

Identification of compound heterozygous DNAH11 variants in a Han-Chinese family with primary ciliary dyskinesia

Primary ciliary dyskinesia (PCD) is a group of genetically and clinically heterogeneous disorders with motile cilia dysfunction. It is clinically characterized by oto‐sino‐pulmonary diseases and subfertility, and half of the patients have situs inversus (Kartagener syndrome). To identify the genetic cause in a Han‐Chinese pedigree, whole‐exome sequencing was conducted in the 37‐year‐old proband, and then, Sanger sequencing was performed on available family members. Minigene splicing assay was applied to verify the impact of the splice‐site variant. Compound heterozygous variants including a splice‐site variant (c.1974‐1G>C, rs1359107415) and a missense variant (c.7787G>A, p.(Arg2596Gln), rs780492669), in the dynein axonemal heavy chain 11 gene (DNAH11) were identified and confirmed as the disease‐associated variants of this lineage. The minigene expression in vitro revealed that the c.1974‐1G>C variant could cause skipping over exon 12, predicted to result in a truncated protein. This discovery may enlarge the DNAH11 variant spectrum of PCD, promote accurate genetic counselling and contribute to PCD diagnosis.

Identification of seven exonic variants in the SLC4A1, ATP6V1B1, and ATP6V0A4 genes that alter RNA splicing by minigene assay

Primary distal renal tubular acidosis (dRTA) is a rare tubular disease associated with variants in SLC4A1, ATP6V0A4, ATP6V1B1, FOXⅠ1, or WDR72 genes. Currently, there is growing evidence that all types of exonic variants can alter splicing regulatory elements, affecting the precursor messenger RNA (pre-mRNA) splicing process. This study was to determine the consequences of variants associated with dRTA on pre-mRNA splicing combined with predictive bioinformatics tools and minigene assay. As a result, among the 15 candidate variants, 7 variants distributed in SLC4A1 (c.1765C>T, p.Arg589Cys), ATP6V1B1 (c.368G>T, p.Gly123Val; c.370C>T, p.Arg124Trp; c.484G>T, p.Glu162* and c.1102G>A, p.Glu368Lys) and ATP6V0A4 genes (c.322C>T, p.Gln108* and c.1572G>A, p.Pro524Pro) were identified to result in complete or incomplete exon skipping by either disruption of exonic splicing enhancers (ESEs) and generation of exonic splicing silencers, or interference with the recognition of the classic splicing site, or both. To our knowledge, this is the first study on pre-mRNA splicing of exonic variants in the dRTA-related genes. These results highlight the importance of assessing the effects of exonic variants at the mRNA level and suggest that minigene analysis is an effective tool for evaluating the effects of splicing on variants in vitro.

Functional Evaluation of Splicing for Variants of Uncertain Significance in Patients with Inherited Retinal Diseases

Inherited retinal diseases (IRD) comprise a heterogeneous set of clinical and genetic disorders that lead to blindness. Given the emerging opportunities in precision medicine and gene therapy, it has become increasingly important to determine whether DNA variants with uncertain significance (VUS) are responsible for patients’ IRD. This research was performed to assess the functional consequence of six VUS identified in patients with IRD. Clinical assessments included an ophthalmic examination, best-corrected visual acuity, and kinetic perimetry. Imaging was acquired with the Optos ultra-widefield camera and spectral domain optical coherence tomography (SD-OCT). Genetic testing was performed by Molecular Vision Laboratories. VUS that were predicted to alter splicing were analyzed with a minigene assay, which revealed that VUS in the genes OPA1, CNGB1, and CLUAP1 altered spicing mechanisms. Due to emerging gene and cell therapies, these results expand the genotype-phenotype correlations for patients diagnosed with an IRD.

Congenital ataxia due to novel variant in ATP8A2

Congenital ataxias are a heterogeneous group of disorders characterized by congenital or early-onset ataxia. Here, we describe two siblings with congenital ataxia, who acquired independent gait by age 4 years. After 16 years of follow-up they presented near normal cognition, cerebellar ataxia, mild pyramidal signs, and dystonia. On exome sequencing, a novel homozygous variant (c.1580-18C > G - intron 17) in ATP8A2 was identified. A new acceptor splice site was predicted by bioinformatics tools, and functionally characterized through a minigene assay. Minigene constructs were generated by PCR-amplification of genomic sequences surrounding the variant of interest and cloning into the pCMVdi vector. Altered splicing was evaluated by expressing these constructs in HEK293T cells. The construct with the c.1580-18C > G homozygous variant produced an aberrant transcript, leading to retention of 17 bp of intron 17, by the use of an alternative acceptor splice site, resulting in a premature stop codon by insertion of four amino acids. These results allowed us to establish this as a disease-causing variant and expand ATP8A2-related disorders to include less severe forms of congenital ataxia.

Splicing Analysis of Exonic OCRL Mutations Causing Lowe Syndrome or Dent-2 Disease

Mutations in the OCRL gene are associated with both Lowe syndrome and Dent-2 disease. Patients with Lowe syndrome present congenital cataracts, mental disabilities and a renal proximal tubulopathy, whereas patients with Dent-2 disease exhibit similar proximal tubule dysfunction but only mild, or no additional clinical defects. It is not yet understood why some OCRL mutations cause the phenotype of Lowe syndrome, while others develop the milder phenotype of Dent-2 disease. Our goal was to gain new insights into the consequences of OCRL exonic mutations on pre-mRNA splicing. Using predictive bioinformatics tools, we selected thirteen missense mutations and one synonymous mutation based on their potential effects on splicing regulatory elements or splice sites. These mutations were analyzed in a minigene splicing assay. Results of the RNA analysis showed that three presumed missense mutations caused alterations in pre-mRNA splicing. Mutation c.741G>T; p.(Trp247Cys) generated splicing silencer sequences and disrupted splicing enhancer motifs that resulted in skipping of exon 9, while mutations c.2581G>A; p.(Ala861Thr) and c.2581G>C; p.(Ala861Pro) abolished a 5′ splice site leading to skipping of exon 23. Mutation c.741G>T represents the first OCRL exonic variant outside the conserved splice site dinucleotides that results in alteration of pre-mRNA splicing. Our results highlight the importance of evaluating the effects of OCRL exonic mutations at the mRNA level.

Analysis of the Impact of Known SPINK1 Missense Variants on Pre-mRNA Splicing and/or mRNA Stability in a Full-Length Gene Assay

It is increasingly appreciated that missense variants may not only alter protein structure and function but may also influence pre-mRNA splicing and/or mRNA stability. Here we explore this issue in the context of currently known SPINK1 missense variants using a full-length gene assay. We demonstrated that 4 (17%) out of 24 variants tested significantly reduced pre-mRNA splicing and/or stability as compared with the wild-type. However, since the strongest effect observed was a 23% reduction from normal, the contribution of SPINK1 missense variants to the clinical phenotype through an impact on mRNA processing alone may be relatively minor compared with their effects in relation to protein structure/function.

Functional analysis by minigene assay of putative splicing variants found in Bardet-Biedl syndrome patients

Bardet–Biedl syndrome (BBS) and Alström syndrome (ALMS) are rare diseases belonging to the group of ciliopathies. Although mutational screening studies of BBS/ALMS cohorts have been extensively reported, little is known about the functional effect of those changes. Thus, splicing variants are estimated to represent 15% of disease‐causing mutations, and there is growing evidence that many exonic changes are really splicing variants misclassified. In this study, we aimed to analyse for the first time several variants in BBS2,ARL6/BBS3,BBS4 and ALMS1 genes predicted to produce aberrant splicing by minigene assay. We found discordance between bioinformatics analysis and experimental data when comparing wild‐type and mutant constructs. Remarkably, we identified nonsense variants presumably resistant to nonsense‐mediated decay, even when a premature termination codon would be introduced in the second amino acid (p.(G2*) mutation in ARL6/BBS3 gene). As a whole, we report one of the first functional studies of BBS/ALMS1 variants using minigene assay, trying to elucidate their role in disease. Functional studies of variants identified in BBS and ALMS patients are essential for their proper classification and subsequent genetic counselling and could also be the start point for new therapeutic approaches, currently based only on symptomatic treatment.

Splicing defects caused by exonic mutations in PKD1 as a new mechanism of pathogenesis in autosomal dominant polycystic kidney disease

The correct splicing of precursor-mRNA depends on the actual splice sites plus exonic and intronic regulatory elements recognized by the splicing machinery. Surprisingly, an increasing number of examples reveal that exonic mutations disrupt the binding of splicing factors to these sequences or generate new splice sites or regulatory elements, causing disease. This contradicts the general assumption that missense mutations disrupt protein function and that synonymous mutations are merely polymorphisms. Autosomal dominant polycystic kidney disease (ADPKD) is a common inherited disorder caused mainly by mutations in the PKD1 gene. Recently, we analyzed a substantial number of PKD1 missense or synonymous mutations to further characterize their consequences on pre-mRNA splicing. Our results showed that one missense and 2 synonymous mutations induce significant defects in pre-mRNA splicing. Thus, it appears that aberrant splicing as a result of exonic mutations is a previously unrecognized cause of ADPKD.

Splicing analysis for exonic and intronic mismatch repair gene variants associated with Lynch syndrome confirms high concordance between minigene assays and patient RNA analyses

A subset of DNA variants causes genetic disease through aberrant splicing. Experimental splicing assays, either RT-PCR analyses of patient RNA or functional splicing reporter minigene assays, are required to evaluate the molecular nature of the splice defect. Here, we present minigene assays performed for 17 variants in the consensus splice site regions, 14 exonic variants outside these regions, and two deep intronic variants, all in the DNA mismatch-repair (MMR) genes MLH1, MSH2, MSH6, and PMS2, associated with Lynch syndrome. We also included two deep intronic variants in APC and PKD2. For one variant (MLH1 c.122A>G), our minigene assay and patient RNA analysis could not confirm the previously reported aberrant splicing. The aim of our study was to further investigate the concordance between minigene splicing assays and patient RNA analyses. For 30 variants results from patient RNA analyses were available, either performed by our laboratory or presented in literature. Some variants were deliberately included in this study because they resulted in multiple aberrant transcripts in patient RNA analysis, or caused a splice effect other than the prevalent exon skip. While both methods were completely concordant in the assessment of splice effects, four variants exhibited major differences in aberrant splice patterns. Based on the present and earlier studies, together showing an almost 100% concordance of minigene assays with patient RNA analyses, we discuss the weight given to minigene splicing assays in the current criteria proposed by InSiGHT for clinical classification of MMR variants.

A rare sequence variant in intron 1 of THAP1 is associated with primary dystonia

Although coding variants in THAP1 have been causally associated with primary dystonia, the contribution of noncoding variants remains uncertain. Herein, we examine a previously identified Intron 1 variant (c.71+9C>A, rs200209986). Among 1672 subjects with mainly adult-onset primary dystonia, 12 harbored the variant in contrast to 1/1574 controls (P < 0.01). Dystonia classification included cervical dystonia (N = 3), laryngeal dystonia (adductor subtype, N = 3), jaw-opening oromandibular dystonia (N = 1), blepharospasm (N = 2), and unclassified (N = 3). Age of dystonia onset ranged from 25 to 69 years (mean = 54 years). In comparison to controls with no identified THAP1 sequence variants, the c.71+9C>A variant was associated with an elevated ratio of Isoform 1 (NM_018105) to Isoform 2 (NM_199003) in leukocytes. In silico and minigene analyses indicated that c.71+9C>A alters THAP1 splicing. Lymphoblastoid cells harboring the c.71+9C>A variant showed extensive apoptosis with relatively fewer cells in the G2 phase of the cell cycle. Differentially expressed genes from lymphoblastoid cells revealed that the c.71+9C>A variant exerts effects on DNA synthesis, cell growth and proliferation, cell survival, and cytotoxicity. In aggregate, these data indicate that THAP1 c.71+9C>A is a risk factor for adult-onset primary dystonia.