Combining full-length gene assay and SpliceAI to interpret the splicing impact of all possible SPINK1 coding variants

BACKGROUND

Single-nucleotide variants (SNVs) within gene coding sequences can significantly impact pre-mRNA splicing, bearing profound implications for pathogenic mechanisms and precision medicine. In this study, we aim to harness the well-established full-length gene splicing assay (FLGSA) in conjunction with SpliceAI to prospectively interpret the splicing effects of all potential coding SNVs within the four-exon SPINK1 gene, a gene associated with chronic pancreatitis.

RESULTS

Our study began with a retrospective analysis of 27 SPINK1 coding SNVs previously assessed using FLGSA, proceeded with a prospective analysis of 35 new FLGSA-tested SPINK1 coding SNVs, followed by data extrapolation, and ended with further validation. In total, we analyzed 67 SPINK1 coding SNVs, which account for 9.3% of the 720 possible coding SNVs. Among these 67 FLGSA-analyzed SNVs, 12 were found to impact splicing. Through detailed comparison of FLGSA results and SpliceAI predictions, we inferred that the remaining 653 untested coding SNVs in the SPINK1 gene are unlikely to significantly affect splicing. Of the 12 splice-altering events, nine produced both normally spliced and aberrantly spliced transcripts, while the remaining three only generated aberrantly spliced transcripts. These splice-impacting SNVs were found solely in exons 1 and 2, notably at the first and/or last coding nucleotides of these exons. Among the 12 splice-altering events, 11 were missense variants (2.17% of 506 potential missense variants), and one was synonymous (0.61% of 164 potential synonymous variants). Notably, adjusting the SpliceAI cut-off to 0.30 instead of the conventional 0.20 would improve specificity without reducing sensitivity.

CONCLUSIONS

By integrating FLGSA with SpliceAI, we have determined that less than 2% (1.67%) of all possible coding SNVs in SPINK1 significantly influence splicing outcomes. Our findings emphasize the critical importance of conducting splicing analysis within the broader genomic sequence context of the study gene and highlight the inherent uncertainties associated with intermediate SpliceAI scores (0.20 to 0.80). This study contributes to the field by being the first to prospectively interpret all potential coding SNVs in a disease-associated gene with a high degree of accuracy, representing a meaningful attempt at shifting from retrospective to prospective variant analysis in the era of exome and genome sequencing.

Identification of protease-sensitive but not misfolding PNLIP variants in familial and hereditary pancreatitis

Mutations in the PNLIP gene have recently been implicated in chronic pancreatitis. Several PNLIP missense variants have been reported to cause protein misfolding and endoplasmic reticulum stress although genetic evidence supporting their association with chronic pancreatitis is currently lacking. Protease-sensitive PNLIP missense variants have also been associated with early-onset chronic pancreatitis although the underlying pathological mechanism remains enigmatic. Herein, we provide new evidence to support the association of protease-sensitive PNLIP variants (but not misfolding PNLIP variants) with pancreatitis. Specifically, we identified protease-sensitive PNLIP variants in 5 of 373 probands (1.3%) with a positive family history of pancreatitis. The protease-sensitive variants, p.F300L and p.I265R, were found to segregate with the disease in three families, including one exhibiting a classical autosomal dominant inheritance pattern. Consistent with previous findings, protease-sensitive variant-positive patients were often characterized by early-onset disease and invariably experienced recurrent acute pancreatitis, although none has so far developed chronic pancreatitis.

A homozygous splice variant in ATP5PO, disrupts mitochondrial complex V function and causes Leigh syndrome in two unrelated families

Mitochondrial complex V plays an important role in oxidative phosphorylation by catalyzing the generation of ATP. Most complex V subunits are nuclear encoded and not yet associated with recognized Mendelian disorders. Using exome sequencing, we identified a rare homozygous splice variant (c.87+3A>G) in ATP5PO, the complex V subunit which encodes the oligomycin sensitivity conferring protein, in three individuals from two unrelated families, with clinical suspicion of a mitochondrial disorder. These individuals had a similar, severe infantile and often lethal multi-systemic disorder that included hypotonia, developmental delay, hypertrophic cardiomyopathy, progressive epileptic encephalopathy, progressive cerebral atrophy, and white matter abnormalities on brain MRI consistent with Leigh syndrome. cDNA studies showed a predominant shortened transcript with skipping of exon 2 and low levels of the normal full-length transcript. Fibroblasts from the affected individuals demonstrated decreased ATP5PO protein, defective assembly of complex V with markedly reduced amounts of peripheral stalk proteins, and complex V hydrolytic activity. Further, expression of human ATP5PO cDNA without exon 2 (hATP5PO-∆ex2) in yeast cells deleted for yATP5 (ATP5PO homolog) was unable to rescue growth on media which requires oxidative phosphorylation when compared to the wild type construct (hATP5PO-WT), indicating that exon 2 deletion leads to a non-functional protein. Collectively, our findings support the pathogenicity of the ATP5PO c.87+3A>G variant, which significantly reduces but does not eliminate complex V activity. These data along with the recent report of an affected individual with ATP5PO variants, add to the evidence that rare biallelic variants in ATP5PO result in defective complex V assembly, function and are associated with Leigh syndrome.

Expanding ACMG variant classification guidelines into a general framework

BACKGROUND

The American College of Medical Genetics and Genomics (ACMG)-recommended five variant classification categories (pathogenic, likely pathogenic, uncertain significance, likely benign, and benign) have been widely used in medical genetics. However, these guidelines are fundamentally constrained in practice owing to their focus upon Mendelian disease genes and their dichotomous classification of variants as being either causal or not. Herein, we attempt to expand the ACMG guidelines into a general variant classification framework that takes into account not only the continuum of clinical phenotypes, but also the continuum of the variants' genetic effects, and the different pathological roles of the implicated genes.

MAIN BODY

As a disease model, we employed chronic pancreatitis (CP), which manifests clinically as a spectrum from monogenic to multifactorial. Bearing in mind that any general conceptual proposal should be based upon sound data, we focused our analysis on the four most extensively studied CP genes, PRSS1, CFTR, SPINK1 and CTRC. Based upon several cross-gene and cross-variant comparisons, we first assigned the different genes to two distinct categories in terms of disease causation: CP-causing (PRSS1 and SPINK1) and CP-predisposing (CFTR and CTRC). We then employed two new classificatory categories, "predisposing" and "likely predisposing", to replace ACMG's "pathogenic" and "likely pathogenic" categories in the context of CP-predisposing genes, thereby classifying all pathologically relevant variants in these genes as "predisposing". In the case of CP-causing genes, the two new classificatory categories served to extend the five ACMG categories whilst two thresholds (allele frequency and functional) were introduced to discriminate "pathogenic" from "predisposing" variants.

CONCLUSION

Employing CP as a disease model, we expand ACMG guidelines into a five-category classification system (predisposing, likely predisposing, uncertain significance, likely benign, and benign) and a seven-category classification system (pathogenic, likely pathogenic, predisposing, likely predisposing, uncertain significance, likely benign, and benign) in the context of disease-predisposing and disease-causing genes, respectively. Taken together, the two systems constitute a general variant classification framework that, in principle, should span the entire spectrum of variants in any disease-related gene. The maximal compliance of our five-category and seven-category classification systems with the ACMG guidelines ought to facilitate their practical application.

Splicing Outcomes of 5' Splice Site GT>GC Variants That Generate Wild-Type Transcripts Differ Significantly Between Full-Length and Minigene Splicing Assays

Combining data derived from a meta-analysis of human disease-associated 5' splice site GT>GC (i.e., +2T>C) variants and a cell culture-based full-length gene splicing assay (FLGSA) of forward engineered +2T>C substitutions, we recently estimated that ∼15-18% of +2T>C variants can generate up to 84% wild-type transcripts relative to their wild-type counterparts. Herein, we analyzed the splicing outcomes of 20 +2T>C variants that generate some wild-type transcripts in two minigene assays. We found a high discordance rate in terms of the generation of wild-type transcripts, not only between FLGSA and the minigene assays but also between the different minigene assays. In the pET01 context, all 20 wild-type minigene constructs generated the expected wild-type transcripts; of the 20 corresponding variant minigene constructs, 14 (70%) generated wild-type transcripts. In the pSPL3 context, only 18 of the 20 wild-type minigene constructs generated the expected wild-type transcripts whereas 8 of the 18 (44%) corresponding variant minigene constructs generated wild-type transcripts. Thus, in the context of a particular type of variant, we raise awareness of the limitations of minigene splicing assays and emphasize the importance of sequence context in regulating splicing. Whether or not our findings apply to other types of splice-altering variant remains to be investigated.

Biallelic variants in MAATS1 encoding CFAP91, a calmodulin-associated and spoke-associated complex protein, cause severe astheno-teratozoospermia and male infertility

BACKGROUND

Multiple morphological abnormalities of the flagella (MMAF) consistently lead to male infertility due to a reduced or absent sperm motility defined as asthenozoospermia. Despite numerous genes recently described to be recurrently associated with MMAF, more than half of the cases analysed remain unresolved, suggesting that many yet uncharacterised gene defects account for this phenotype METHODS: Exome sequencing was performed on 167 infertile men with an MMAF phenotype. Immunostaining and transmission electron microscopy (TEM) in sperm cells from affected individuals were performed to characterise the ultrastructural sperm defects. Gene inactivation using RNA interference (RNAi) was subsequently performed in Trypanosoma.

RESULTS

We identified six unrelated affected patients carrying a homozygous deleterious variants in MAATS1, a gene encoding CFAP91, a calmodulin-associated and spoke-associated complex (CSC) protein. TEM and immunostaining experiments in sperm cells showed severe central pair complex (CPC) and radial spokes defects. Moreover, we confirmed that the WDR66 protein is a physical and functional partner of CFAP91 into the CSC. Study of Trypanosoma MAATS1's orthologue (TbCFAP91) highlighted high sequence and structural analogies with the human protein and confirmed the axonemal localisation of the protein. Knockdown of TbCFAP91 using RNAi impaired flagellar movement led to CPC defects in Trypanosoma as observed in humans.

CONCLUSIONS

We showed that CFAP91 is essential for normal sperm flagellum structure and function in human and Trypanosoma and that biallelic variants in this gene lead to severe flagellum malformations resulting in astheno-teratozoospermia and primary male infertility.

First estimate of the scale of canonical 5' splice site GT>GC variants capable of generating wild-type transcripts

It has long been known that canonical 5' splice site (5'SS) GT>GC variants may be compatible with normal splicing. However, to date, the actual scale of canonical 5'SSs capable of generating wild-type transcripts in the case of GT>GC substitutions remains unknown. Herein, combining data derived from a meta-analysis of 45 human disease-causing 5'SS GT>GC variants and a cell culture-based full-length gene splicing assay of 103 5'SS GT>GC substitutions, we estimate that ~15-18% of canonical GT 5'SSs retain their capacity to generate between 1% and 84% normal transcripts when GT is substituted by GC. We further demonstrate that the canonical 5'SSs in which substitution of GT by GC-generated normal transcripts exhibit stronger complementarity to the 5' end of U1 snRNA than those sites whose substitutions of GT by GC did not lead to the generation of normal transcripts. We also observed a correlation between the generation of wild-type transcripts and a milder than expected clinical phenotype but found that none of the available splicing prediction tools were capable of reliably distinguishing 5'SS GT>GC variants that generated wild-type transcripts from those that did not. Our findings imply that 5'SS GT>GC variants in human disease genes may not invariably be pathogenic.

Do pregnancies reduce iron overload in HFE hemochromatosis women? results from an observational prospective study

Background

HFE hemochromatosis is an inborn error of iron metabolism linked to a defect in the regulation of hepcidin synthesis. This autosomal recessive disease typically manifests later in women than men. Although it is commonly stated that pregnancy is, with menses, one of the factors that offsets iron accumulation in women, no epidemiological study has yet supported this hypothesis. The aim of our study was to evaluate the influence of pregnancy on expression of the predominant HFE p.[Cys282Tyr];[Cys282Tyr] genotype.

Methods

One hundred and forty p.Cys282Tyr homozygous women enrolled in a phlebotomy program between 2004 and 2011 at a blood centre in western Brittany (France) were included in the study. After checking whether the disease expression was delayed in women than in men in our study, the association between pregnancy and iron overload was assessed using multivariable regression analysis.

Results

Our study confirms that women with HFE hemochromatosis were diagnosed later than men cared for during the same period (52.6 vs. 47.4 y., P < 0.001). Compared to no pregnancy, having at least one pregnancy was not associated with lower iron markers. In contrast, the amount of iron removed by phlebotomies appeared significantly higher in women who had at least one pregnancy (e^β = 1.50, P = 0.047). This relationship disappeared after adjustment for confounding factors (e^β = 1.35, P = 0.088).

Conclusions

Our study shows that pregnancy status has no impact on iron markers level, and is not in favour of pregnancy being a protective factor in progressive iron accumulation. Our results are consistent with recent experimental data suggesting that the difference in disease expression observed between men and women may be explained by other factors such as hormones.

Heterozygous Mutations in BMP6 Pro-peptide Lead to Inappropriate Hepcidin Synthesis and Moderate Iron Overload in Humans

BACKGROUND & AIMS

Hereditary hemochromatosis is a heterogeneous group of genetic disorders characterized by parenchymal iron overload. It is caused by defective expression of liver hepcidin, the main regulator of iron homeostasis. Iron stimulates the gene encoding hepcidin (HAMP) via the bone morphogenetic protein (BMP)6 signaling to SMAD. Although several genetic factors have been found to cause late-onset hemochromatosis, many patients have unexplained signs of iron overload. We investigated BMP6 function in these individuals.

METHODS

We sequenced the BMP6 gene in 70 consecutive patients with a moderate increase in serum ferritin and liver iron levels who did not carry genetic variants associated with hemochromatosis. We searched for BMP6 mutations in relatives of 5 probands and in 200 healthy individuals (controls), as well as in 2 other independent cohorts of hyperferritinemia patients. We measured serum levels of hepcidin by liquid chromatography-tandem mass spectrometry and analyzed BMP6 in liver biopsy specimens from patients by immunohistochemistry. The functions of mutant and normal BMP6 were assessed in transfected cells using immunofluorescence, real-time quantitative polymerase chain reaction, and immunoblot analyses.

RESULTS

We identified 3 heterozygous missense mutations in BMP6 (p.Pro95Ser, p.Leu96Pro, and p.Gln113Glu) in 6 unrelated patients with unexplained iron overload (9% of our cohort). These mutations were detected in less than 1% of controls. p.Leu96Pro also was found in 2 patients from the additional cohorts. Family studies indicated dominant transmission. Serum levels of hepcidin were inappropriately low in patients. A low level of BMP6, compared with controls, was found in a biopsy specimen from 1 patient. In cell lines, the mutated residues in the BMP6 propeptide resulted in defective secretion of BMP6; reduced signaling via SMAD1, SMAD5, and SMAD8; and loss of hepcidin production.

CONCLUSIONS

We identified 3 heterozygous missense mutations in BMP6 in patients with unexplained iron overload. These mutations lead to loss of signaling to SMAD proteins and reduced hepcidin production. These mutations might increase susceptibility to mild-to-moderate late-onset iron overload.

Interlaboratory development and validation of a HRM method applied to the detection of JAK2 exon 12 mutations in polycythemia vera patients

Background

Myeloproliferative disorders are characterized by clonal expansion of normal mature blood cells. Acquired mutations giving rise to constitutive activation of the JAK2 tyrosine kinase has been shown to be present in the majority of patients. Since the demonstration that the V617F mutation in the exon 14 of the JAK2 gene is present in about 90% of patients with Polycythemia Vera (PV), the detection of this mutation has become a key tool for the diagnosis of these patients. More recently, additional mutations in the exon 12 of the JAK2 gene have been described in 5 to 10% of the patients with erythrocytosis. According to the updated WHO criteria the presence of these mutations should be looked for in PV patients with no JAK2 V617F mutation. Reliable and accurate methods dedicated to the detection of these highly variable mutations are therefore necessary.

Methods/Findings

For these reasons we have defined the conditions of a High Resolution DNA Melting curve analysis (HRM) method able to detect JAK2 exon 12 mutations. After having validated that the method was able to detect mutated patients, we have verified that it gave reproducible results in repeated experiments, on DNA extracted from either total blood or purified granulocytes. This HRM assay was further validated using 8 samples bearing different mutant sequences in 4 different laboratories, on 3 different instruments.

Conclusion

The assay we have developed is thus a valid method, adapted to routine detection of JAK2 exon 12 mutations with highly reproducible results.

DNA damage-induced down-regulation of human Cdc25C and Cdc2 is mediated by cooperation between p53 and maintenance DNA (cytosine-5) methyltransferase 1

The Cdc25C phosphatase mediates cellular entry into mitosis in mammalian cells. Cdc25C activates Cdc2 for entry into mitosis by dephosphorylating Thr and Tyr at the site of inhibitory phosphorylation. The Cdc25C gene contains tumor suppressor p53 binding sites and is demonstrated to contribute to the p53-dependent cell cycle arrest upon DNA damage. Here we show that both Cdc25C and Cdc2 were down-regulated in wild-type HCT116 cells but not in p53-null, DNMT1-null or DNMT1and DNMT3b-null cells, upon p53 stabilization following doxorubicin-mediated DNA damage. Furthermore, zebularine, a drug that selectively traps and depletes nuclear DNMT1 and DNMT3b, relieved p53-mediated repression of endogenous Cdc25C and Cdc2. Methylation analysis of the Cdc25C and Cdc2 promoter displayed internal CG methylation proximal to the p53 binding site upon DNA damage in a p53-dependent manner. Chromatin immunoprecipitation of doxorubicin treated wild-type HCT116 cells showed the presence of DNMT1, p53, H3K9me2, and the transcriptional repressor HDAC1 on the Cdc25C and Cdc2 promoters, suggesting their involvement as repressive complexes in Cdc25C and Cdc2 gene silencing. Thus, the general mechanism of p53-mediated gene repression may involve recruitment of other repressive factors.

HAMP as a modifier gene that increases the phenotypic expression of the HFE pC282Y homozygous genotype

Hereditary hemochromatosis is a genetically heterogeneous disease of iron metabolism. The most common form of the disorder is an adult-onset form that has mainly been associated with the HFE pC282Y/pC282Y genotype. The phenotypic expression of this genotype is very heterogeneous and could be modulated by both environmental factors and modifier genes. The non-HFE hereditary hemochromatosis forms include a juvenile onset form associated with mutations in HAMP. From a cohort of 392 C282Y homozygous patients, we found 5 carriers of an additional HAMP mutation at the heterozygous state (pR59G, pG71D, or pR56X). We found that iron indices of these 5 patients were among the most elevated of the cohort. Moreover, we specified that the HAMP mutations were not detected in 300 control subjects. These results revealed that mutations in HAMP might increase the phenotypic expression of the pC282Y/pC282Y genotype. From a cohort of 31 patients with at least one chromosome lacking an HFE mutation, we further identified 4 males carrying a heterozygous HAMP mutation (pR59G or pG71D). Based on a digenic model of inheritance, these data suggest that the association of heterozygous mutations in the HFE and HAMP genes could lead, at least in some cases, to an adult-onset form of primary iron overload.

Complete Scanning of the Hereditary Hemochromatosis Gene (HFE) by Use of Denaturing HPLC

Background: Between 4% and 35% of hereditary hemochromatosis (HC) probands are C282Y or H63D heterozygotes or lack both of these two common HFE mutations, and 15 novel HFE mutations have been described recently. We evaluated denaturing HPLC (DHPLC) for screening of the whole HFE coding region and further defined whether HC probands with an incomplete HFE genotype carry uncommon mutations.

Methods: Analytical conditions for each coding exon were determined by a combination of computer melting profile predictions and experimental melting curves. To test accuracy for scanning the complete HFE coding region and optimize DHPLC running conditions, each melting domain was investigated with at least one mutation or one polymorphism as reference. We tested 100 DNA samples harboring the C282Y, H63D, or S65C mutations and 17 artificially created positive controls that carried either 1 of the 14 other known HFE mutations or 3 selected polymorphisms.

Results: Investigations on each of the coding exons 1, 2, 4, 5, and 6 could be performed at one analysis temperature. Coding exon 3 displayed a more complex melting profile and required two analysis temperatures. DHPLC detected all known HFE mutations as well as the three selected polymorphisms.

Conclusions: DHPLC can be used to scan the HFE gene in HC probands in whom at least one chromosome lacks an assigned mutation.