Splicing analysis of 26 F8 nucleotide variations using a minigene assay

An intronic RNA element modulates Factor VIII exon-16 splicing

Pathogenic variants in the human Factor VIII (F8) gene cause Hemophilia A (HA). Here, we investigated the impact of 97 HA-causing single-nucleotide variants on the splicing of 11 exons from F8. For the majority of F8 exons, splicing was insensitive to the presence of HA-causing variants. However, splicing of several exons, including exon-16, was impacted by variants predicted to alter exonic splicing regulatory sequences. Using exon-16 as a model, we investigated the structure-function relationship of HA-causing variants on splicing. Intriguingly, RNA chemical probing analyses revealed a three-way junction structure at the 3'-end of intron-15 (TWJ-3-15) capable of sequestering the polypyrimidine tract. We discovered antisense oligonucleotides (ASOs) targeting TWJ-3-15 partially rescue splicing-deficient exon-16 variants by increasing accessibility of the polypyrimidine tract. The apical stem loop region of TWJ-3-15 also contains two hnRNPA1-dependent intronic splicing silencers (ISSs). ASOs blocking these ISSs also partially rescued splicing. When used in combination, ASOs targeting both the ISSs and the region sequestering the polypyrimidine tract, fully rescue pre-mRNA splicing of multiple HA-linked variants of exon-16. Together, our data reveal a putative RNA structure that sensitizes F8 exon-16 to aberrant splicing.

Whole F8 gene sequencing combined with splicing functional analyses led to a substantial increase of the molecular diagnosis yield for non-severe haemophilia A

INTRODUCTION

Conventional genetic investigation fails to identify the F8 causal variant in 2.5%-10% of haemophilia A (HA) patients with non-severe phenotypes. In these cases, F8 deep intronic variants could be causal.

AIM

To identify pathogenic F8 deep intronic variants in genetically unresolved families with non-severe HA analysed in the haematology laboratory of the Hospices Civils de Lyon.

METHODS

The whole F8 was analysed by next generation sequencing. The pathogenic impact of candidate variants identified was assessed using both in silico analysis (MaxEntScan and spliceAI) and functional analysis (RNA or minigene assay).

RESULTS

Sequencing was performed in 49/55 families included for which a DNA sample from a male propositus was available. In total, 33 candidate variants from 43 propositi were identified. These variants corresponded to 31 single nucleotide substitutions, one 173-bp deletion, and an 869-bp tandem triplication. No candidate variant was found in six propositi. The most frequent variants found were the association of [c.2113+1154G>C and c.5374-304C>T], identified in five propositi, and the c.2114-6529C>G identified in nine propositi. Four variants had been previously described as HA-causing. Splicing functional assay found a deleterious impact for 11 substitutions (c.671-94G>A, c.788-312A>G, c.2113+1154G>C, c.2114-6529C>G, c.5999-820A>T, c.5999-786C>A, c.5999-669G>T, c.5999-669G>A, c.5999-669G>C, c.6900+4104A>C, and c.6901-2992A>G). The HA-causing variant was identified in 33/49 (67%) cases. In total, F8 deep intronic variants caused 8.8% of the non-severe HA among the 1643 families analysed in our laboratory.

CONCLUSION

The results emphasise the value of whole F8 gene sequencing combined with splicing functional analyses to improve the diagnosis yield for non-severe HA.

Rescue of blood coagulation Factor VIII exon-16 mis-splicing by antisense oligonucleotides

The human Factor VIII ( F8 ) protein is essential for the blood coagulation cascade and specific F8 mutations cause the rare bleeding disorder Hemophilia A (HA). Here, we investigated the impact of HA-causing single-nucleotide mutations on F8 pre-mRNA splicing. We found that 14/97 (∼14.4%) coding sequence mutations tested in our study induced exon skipping. Splicing patterns of 4/11 (∼36.4%) F8 exons tested were especially sensitive to the presence of common disease-causing mutations. RNA-chemical probing analyses revealed a three-way junction structure at the 3' end of intron 15 (TWJ-3-15). TWJ-3-15 sequesters the polypyrimidine tract, a key determinant of 3' splice site strength. Using exon-16 of the F8 gene as a model, we designed specific antisense oligonucleotides (ASOs) that target TWJ-3-15 and identified three that promote the splicing of F8 exon-16. Interaction of TWJ-3-15 with ASOs increases accessibility of the polypyrimidine tract and inhibits the binding of hnRNPA1-dependent splicing silencing factors. Moreover, ASOs targeting TWJ-3-15 rescue diverse splicing-sensitive HA-causing mutations, most of which are distal to the 3' splice site being impacted. The TWJ-3-15 structure and its effect on mRNA splicing provide a model for HA etiology in patients harboring specific F8 mutations and provide a framework for precision RNA-based HA therapies.

Identification of a Splicing Variant c.3813-3A>G in NPHP3 by Reanalysis of Whole Exome Sequencing in a Chinese Boy with Nephronophthisis

Nephronophthisis is an autosomal recessive cystic kidney disease characterized by tubular injury and commonly results in kidney failure. We reported a case of 4-year-old Chinese boy presented with severe anemia, kidney, and liver dysfunction. Whole exome sequencing (WES) was performed to identify the candidate variant with a negative result initially. After complete collection of clinical information, reanalysis of WES identified a homozygous NPHP3 variant c.3813-3A>G (NM_153240.4). The effect on mRNA splicing of the intronic variant was predicted through software (three in silico splice tools). Furthermore, in vitro minigene assay was conducted to validate the predicted deleterious effects of the intronic variant. All of the splice prediction programs and minigene assay indicated that the variant had an impact on the normal splicing pattern of NPHP3. Our study confirmed the effect of the c.3813-3A>G variant on NPHP3 splicing in vitro, which gives additional evidence for the clinical significance of the variant and provides a basis for genetic diagnosis of nephronophthisis 3. In addition, we think that it is essential to reanalyze WES data after the complete clinical information collection to avoid missing some important candidate variants.

Whole F9 gene sequencing identified deep intronic variations in genetically unresolved hemophilia B patients

BACKGROUND

The disease-causative variant remains unidentified in approximately 0.5% to 2% of hemophilia B patients using conventional genetic investigations, and F9 deep intronic variations could be responsible for these phenotypes.

OBJECTIVES

This study aimed to characterize deep intronic variants in hemophilia B patients for whom genetic investigations failed.

METHODS

We performed whole F9 sequencing in 17 genetically unsolved hemophilia B patients. The pathogenic impact of the candidate variants identified was studied using both in silico analysis (MaxEntScan and spliceAI) and minigene assay.

RESULTS

In total, 9 candidate variants were identified in 15 patients; 7 were deep intronic substitutions and 2 corresponded to insertions of mobile elements. The most frequent variants found were c.278-1806A>C and the association of c.278-1244A>G and c.392-864T>G, identified in 4 and 6 unrelated individuals, respectively. In silico analysis predicted splicing impact for 4 substitutions (c.278-1806A>C, c.392-864T>G, c.724-2385G>T, c.723+4297T>A). Minigene assay showed a deleterious splicing impact for these 4 substitutions and also for the c.278-1786_278-1785insLINE. In the end, 5 variants were classified as likely pathogenic using the American College of Medical Genetics and Genomics guidelines, and 4 as of unknown significance. Thus, the hemophilia B-causing variant was identified in 13/17 (76%) families.

CONCLUSION

We elucidated the causing defect in three-quarters of the families included in this study, and we reported new F9 deep intronic variants that can cause hemophilia B.

Presumed COL4A3/COL4A4 Missense/Synonymous Variants Induce Aberrant Splicing

Background

The incorrect interpretation of missense and synonymous variants can lead to improper molecular diagnosis and subsequent faulty genetic counselling. The aim of this study was to evaluate the pathogenicity of presumed COL4A3/COL4A4 missense and synonymous variants detected by next-generation sequencing to provide evidence for diagnosis and genetic counselling.

Methods

Patients' clinical findings and genetic data were analysed retrospectively. An in vitro minigene assay was conducted to assess the effect of presumed COL4A3/COL4A4 missense and synonymous variants on RNA splicing.

Results

Five unclassified COL4A3/COL4A4 variants, which were detected in five of 343 patients with hereditary kidney diseases, were analysed. All of them were predicted to affect splicing by Human Splicing Finder. The presumed COL4A3 missense variant c.4793T > G [p. (Leu1598Arg)] resulted in a loss of alternative full-length transcript during the splicing process. The COL4A3 transcript carried synonymous variant c.765G > A [p. (Thr255Thr)], led to an in-frame deletion of exon 13. Nevertheless, variants c.3566G > A [p. (Gly1189Glu)] in COL4A3 and c.3990G > A [p. (Pro1330Pro)], c.4766C > T [p. (Pro1589Leu)] in COL4A4 exhibited no deleterious effect on splicing. Among the five patients harbouring the abovementioned COL4A3/COL4A4 variants, three patients were genetically diagnosed with autosomal recessive Alport syndrome, one patient was highly suspected of having thin basement membrane nephropathy, and the other patient was clinically diagnosed with Alport syndrome.

Conclusions

COL4A3 presumed missense variant p. (Leu1598Arg) and synonymous variant p. (Thr255Thr) affect RNA splicing, which highlights the prime importance of transcript analysis of unclassified exonic sequence variants for better molecular diagnosis and genetic counselling. Meanwhile, the reliability of splicing predictions by predictive tools for exonic substitutions needs to be improved.

RNA splicing analysis contributes to reclassifying variants of uncertain significance and improves the diagnosis of monogenic disorders

BACKGROUND

Numerous variants of uncertain significance (VUSs) have been identified by whole exome sequencing in clinical practice. However, VUSs are not currently considered medically actionable.

OBJECTIVE

To assess the splicing patterns of 49 VUSs in 48 families identified clinically to improve genetic counselling and family planning.

METHODS

Forty-nine participants with 49 VUSs were recruited from the Reproductive and Genetic Hospital of CITIC-Xiangya. Bioinformatic analysis was performed to preliminarily predict the splicing effects of these VUSs. RT-PCR and minigene analysis were used to assess the splicing patterns of the VUSs. According to the results obtained, couples opted for different methods of reproductive interventions to conceive a child, including prenatal diagnosis and preimplantation genetic testing (PGT).

RESULTS

Eleven variants were found to alter pre-mRNA splicing and one variant caused nonsense-mediated mRNA decay, which resulted in the reclassification of these VUSs as likely pathogenic. One couple chose to undergo in vitro fertilisation with PGT treatment; a healthy embryo was transferred and the pregnancy is ongoing. Three couples opted for natural pregnancy with prenatal diagnosis. One couple terminated the pregnancy because the fetus was affected by short-rib thoracic dysplasia and harboured the related variant. The infants of the other two couples were born and were healthy at their last recorded follow-up.

CONCLUSION

RNA splicing analysis is an important method to assess the impact of sequence variants on splicing in clinical practice and can contribute to the reclassification of a significant proportion of VUSs. RNA splicing analysis should be considered for genetic disease diagnostics.

Dissection of pleiotropic effects of variants in and adjacent to F8 exon 19 and rescue of mRNA splicing and protein function

The pathogenic significance of nucleotide variants commonly relies on nucleotide position within the gene, with exonic changes generally attributed to quantitative or qualitative alteration of protein biosynthesis, secretion, activity, or clearance. However, these changes may exert pleiotropic effects on both protein biology and mRNA splicing due to the overlapping of the amino acid and splicing codes, thus shaping the disease phenotypes. Here, we focused on hemophilia A, in which the definition of F8 variants' causative role and association to bleeding phenotypes is crucial for proper classification, genetic counseling, and management of affected individuals. We extensively characterized a large panel of hemophilia A-causing variants (n = 30) within F8 exon 19 by combining and comparing in silico and recombinant expression analyses. We identified exonic variants with pleiotropic effects and dissected the altered protein features of all missense changes. Importantly, results from multiple prediction algorithms provided qualitative results, while recombinant assays allowed us to correctly infer the likely phenotype severity for 90% of variants. Molecular characterization of pathogenic variants was also instrumental for the development of tailored correction approaches to rescue splicing affecting variants or missense changes impairing protein folding. A single engineered U1snRNA rescued mRNA splicing of nine different variants and the use of a chaperone-like drug resulted in improved factor VIII protein secretion for four missense variants. Overall, dissection of the molecular mechanisms of a large panel of HA variants allowed precise classification of HA-affected individuals and favored the development of personalized therapeutic approaches.

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.

Splicing analysis of SLC40A1 missense variations and contribution to hemochromatosis type 4 phenotypes

Hemochromatosis type 4, or ferroportin disease, is considered as the second leading cause of primary iron overload after HFE-related hemochromatosis. The disease, which is predominantly associated with missense variations in the SLC40A1 gene, is characterized by wide clinical heterogeneity. We tested the possibility that some of the reported missense mutations, despite their positions within exons, cause splicing defects. Fifty-eight genetic variants were selected from the literature based on two criteria: a precise description of the nucleotide change and individual evidence of iron overload. The selected variants were investigated by different in silico prediction tools and prioritized for midigene splicing assays. Of the 15 variations tested in vitro, only two were associated with splicing changes. We confirm that the c.1402G>A transition (p.Gly468Ser) disrupts the exon 7 donor site, leading to the use of an exonic cryptic splicing site and the generation of a truncated reading frame. We observed, for the first time, that the p.Gly468Ser substitution has no effect on the ferroportin iron export function. We demonstrate alternative splicing of exon 5 in different cell lines and show that the c.430A>G (p.Asn144Asp) variant promotes exon 5 inclusion. This could be part of a gain-of-function mechanism. We conclude that splicing mutations rarely contribute to hemochromatosis type 4 phenotypes. An in-depth investigation of exon 5 auxiliary splicing sequences may help to elucidate the mechanism by which splicing regulatory proteins regulate the production of the full length SLC40A1 transcript and to clarify its physiological importance.

The challenge of genetically unresolved haemophilia A patients: Interest of the combination of whole F8 gene sequencing and functional assays

BACKGROUND

The causative variant remains unidentified in 2%-5% of haemophilia A (HA) patients despite an exhaustive sequencing of the full F8 coding sequence, splice consensus sequences, 5'/3' untranslated regions and copy number variant (CNV) analysis. Next-generation sequencing (NGS) has provided significant improvements for a complete F8 analysis.

AIM

The aim of this study was to identify and characterize pathogenic non-coding variants in F8 of 15 French and Canadian HA patients genetically unresolved, through the use of NGS, mRNA sequencing and functional confirmation of aberrant splicing.

METHODS

We sequenced the entire F8 gene using an NGS capture method. We analysed F8 mRNA in order to detect aberrant transcripts. The pathogenic effect of candidate intronic variants was further confirmed using a minigene assay.

RESULTS

After bioinformatic analysis, 11 deep intronic variants were identified in 13 patients (8 new variants and 3 previously reported). Three variants were confirmed to be likely pathogenic with the presence of an aberrant transcript during mRNA analysis and minigene assay. We also found a small intronic deletion in 6 patients, recently described as causing mild HA.

CONCLUSION

With this comprehensive work combining NGS and functional assays, we report new deep intronic variants that cause HA through splicing alteration mechanism. Functional analyses are critical to confirm the pathogenic effect of these variants and will be invaluable in the future to study the large number of variants of uncertain significance that may affect splicing that will be found in the human genome.

Deep intronic F8 c.5999-27A>G variant causes exon 19 skipping and leads to moderate hemophilia A

: Hemophilia A, an X-linked recessive bleeding disorder, is caused by mutations of F8 gene. In about 2% hemophilia A patients, no exonic mutation of F8 gene was found. We aimed to identify deep intronic mutations of F8 gene. We reanalyzed the next-generation sequencing data of six hemophilia A patients with negative F8 variant in either coding region or splice site. Deep intronic F8 c.5999-27A>G variant (NM_000132.3) was found in two unrelated moderate hemophilia A patients from different region, and one patient's mother was mild hemophilia A patient. Splice site prediction algorithms showed no impact of this variant on F8 mRNA splicing of exon 19, including Human Splicing Finder 3.1, NNSPLICE 0.9, NetGene2, and Transcript-inferred Pathogenicity score. Exonic splicing enhancer was predicted by ESEfinder, and no difference was found between the wild type and mutant sequence. The branch point predicted by SVM-BPfinder suggested that F8 c.5999-27A>G variant may disrupt the branch point in intron 18 and affect the acceptor site splicing of F8 exon 19. Sanger sequencing of F8 cDNA from peripheral blood mononuclear cells confirmed that F8 c.5999-27A>G variant caused F8 exon 19 skipping in proband and his mother. Skewed X chromosome inactivation was found in another X chromosome of this mother, combined with F8 c.5999-27A>G variant in trans. In conclusion, our study suggests that deep intronic F8 c.5999-27A>G variant may be responsible for F8 exon 19 skipping and lead to moderate hemophilia A. Systematic reanalysis of next-generation sequencing data could promote the diagnostic yields.

The highly prevalent deletions in F8 intron 13 found in French mild hemophilia A patients result from both founder effect and recurrent de novo events

BACKGROUND

Recently, our group has reported a 13-bp deletion in a poly(T)-track in the F8 intron 13 as the causative variant in approximately 6% of all cases of mild haemophilia A (HA) in France. The systematic screening of mild HA patients for this deletion identified individuals carrying deletions from 9 to 14-bp in the same region.

AIMS

To demonstrate that these highly prevalent deletions could result from a recurrent molecular mechanism and to determine the clinical significance of deletions other than 13-bp in size.

METHODS

Haplotype analysis using five polymorphic markers was performed in 71 unrelated French mild hemophilia A patients. Minigene analysis was performed to study the splicing impact of deletions from 1 to 14-bp.

RESULTS

A peculiar haplotype (H1) was identified in 22.5% of patients carrying the 13-bp deletion. Haplotypes differing from H1 only for the two most distal markers were found in more than the half of patients. These results confirmed the founder effect origin for the 13-bp deletion. However, the 9 patients carrying other sizes of deletion had a different haplotype suggesting that these deletions arose independently. Supporting the recurrent mechanism hypothesis, similar deletions were also found in 3/19 genetically unresolved mild Canadian patients. In vitro splicing analysis confirmed that deletions larger than 9-bp had a deleterious impact on splicing of F8 transcript.

CONCLUSION

We demonstrated that the poly(T)-track in F8 intron 13 is a deletion hotspot. We recommend that deletions in this region should be specifically investigated in all genetically unresolved mild HA patients.