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.

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.

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.

SPiP: Splicing Prediction Pipeline, a machine learning tool for massive detection of exonic and intronic variant effects on mRNA splicing.

Modeling splicing is essential for tackling the challenge of variant interpretation as each nucleotide variation can be pathogenic by affecting pre-mRNA splicing via disruption/creation of splicing motifs such as 5'/3' splice sites, branch sites, or splicing regulatory elements. Unfortunately, most in silico tools focus on a specific type of splicing motif, which is why we developed the Splicing Prediction Pipeline (SPiP) to perform, in one single bioinformatic analysis based on a machine learning approach, a comprehensive assessment of the variant effect on different splicing motifs. We gathered a curated set of 4616 variants scattered all along the sequence of 227 genes, with their corresponding splicing studies. The Bayesian analysis provided us with the number of control variants, that is, variants without impact on splicing, to mimic the deluge of variants from high-throughput sequencing data. Results show that SPiP can deal with the diversity of splicing alterations, with 83.13% sensitivity and 99% specificity to detect spliceogenic variants. Overall performance as measured by area under the receiving operator curve was 0.986, better than SpliceAI and SQUIRLS (0.965 and 0.766) for the same data set. SPiP lends itself to a unique suite for comprehensive prediction of spliceogenicity in the genomic medicine era. SPiP is available at: https://sourceforge.net/projects/splicing-prediction-pipeline/.

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.

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.

Missense RHD single nucleotide variants induce weakened D antigen expression by altering splicing and/or protein expression.

BACKGROUND: Although D variant phenotype is known to be due to genetic defects, including rare missense single nucleotide variants (SNVs), within the RHD gene, few studies have addressed the molecular and cellular mechanisms driving this altered expression. We and others showed previously that splicing is commonly disrupted by SNVs in constitutive splice sites and their vicinity. We thus sought to investigate whether rare missense SNVs located in "deep" exonic regions could also impair this mechanism. STUDY DESIGN AND METHODS: Forty-six missense SNVs reported within exons 6 and 7 were first selected from the Human RhesusBase. Their respective effect on splicing was assessed by using an in vitro assay. An RhD-negative cell model was further generated by using the CRISPR-Cas9 approach. RhD-mutated proteins were overexpressed in the newly created model, and cell membrane expression of the D antigen was measured by flow cytometry. RESULTS: Minigene splicing assay showed that 14 of 46 (30.4%) missense SNVs alter splicing. Very interestingly, further investigation of two missense SNVs, which both affect codon 338 and confer a weak D phenotype, showed various mechanisms: c.1012C>G (p.Leu338Val) disrupts splicing only, while c.1013T>C (p.Leu338Pro) alters only the protein structure, in agreement with in silico prediction tools and 3D protein structure visualization. CONCLUSION: Our functional data set suggests that missense SNVs damage quantitatively D antigen expression by, at least, two different mechanisms (splicing alteration and protein destabilization) that may act independently. These data thereby contribute to extend the current knowledge of the molecular mechanisms governing weakened D expression.

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.

Insights into the Role of the Discontinuous TM7 Helix of Human Ferroportin through the Prism of the Asp325 Residue.

The negatively charged Asp325 residue has proved to be essential for iron export by human (HsFPN1) and primate Philippine tarsier (TsFpn) ferroportin, but its exact role during the iron transport cycle is still to be elucidated. It has been posited as being functionally equivalent to the metal ion-coordinating residue His261 in the C-lobe of the bacterial homolog BbFpn, but the two residues arise in different sequence motifs of the discontinuous TM7 transmembrane helix. Furthermore, BbFpn is not subject to extracellular regulation, contrary to its mammalian orthologues which are downregulated by hepcidin. To get further insight into the molecular mechanisms related to iron export in mammals in which Asp325 is involved, we investigated the behavior of the Asp325Ala, Asp325His, and Asp325Asn mutants in transiently transfected HEK293T cells, and performed a comparative structural analysis. Our biochemical studies clearly distinguished between the Asp325Ala and Asp325His mutants, which result in a dramatic decrease in plasma membrane expression of FPN1, and the Asp325Asn mutant, which alters iron egress without affecting protein localization. Analysis of the 3D structures of HsFPN1 and TsFpn in the outward-facing (OF) state indicated that Asp325 does not interact directly with metal ions but is involved in the modulation of Cys326 metal-binding capacity. Moreover, models of the architecture of mammalian proteins in the inward-facing (IF) state suggested that Asp325 may form an inter-lobe salt-bridge with Arg40 (TM1) when not interacting with Cys326. These findings allow to suggest that Asp325 may be important for fine-tuning iron recognition in the C-lobe, as well as for local structural changes during the IF-to-OF transition at the extracellular gate level. Inability to form a salt-bridge between TM1 and TM7b during iron translocation could lead to protein instability, as shown by the Asp325Ala and Asp325His mutants.

Assessment of branch point prediction tools to predict physiological branch points and their alteration by variants.

BACKGROUND: Branch points (BPs) map within short motifs upstream of acceptor splice sites (3'ss) and are essential for splicing of pre-mature mRNA. Several BP-dedicated bioinformatics tools, including HSF, SVM-BPfinder, BPP, Branchpointer, LaBranchoR and RNABPS were developed during the last decade. Here, we evaluated their capability to detect the position of BPs, and also to predict the impact on splicing of variants occurring upstream of 3'ss. RESULTS: We used a large set of constitutive and alternative human 3'ss collected from Ensembl (n = 264,787 3'ss) and from in-house RNAseq experiments (n = 51,986 3'ss). We also gathered an unprecedented collection of functional splicing data for 120 variants (62 unpublished) occurring in BP areas of disease-causing genes. Branchpointer showed the best performance to detect the relevant BPs upstream of constitutive and alternative 3'ss (99.48 and 65.84% accuracies, respectively). For variants occurring in a BP area, BPP emerged as having the best performance to predict effects on mRNA splicing, with an accuracy of 89.17%. CONCLUSIONS: Our investigations revealed that Branchpointer was optimal to detect BPs upstream of 3'ss, and that BPP was most relevant to predict splicing alteration due to variants in the BP area.

Novel diagnostic tool for prediction of variant spliceogenicity derived from a set of 395 combined in silico/in vitro studies: an international collaborative effort.

Variant interpretation is the key issue in molecular diagnosis. Spliceogenic variants exemplify this issue as each nucleotide variant can be deleterious via disruption or creation of splice site consensus sequences. Consequently, reliable in silico prediction of variant spliceogenicity would be a major improvement. Thanks to an international effort, a set of 395 variants studied at the mRNA level and occurring in 5' and 3' consensus regions (defined as the 11 and 14 bases surrounding the exon/intron junction, respectively) was collected for 11 different genes, including BRCA1, BRCA2, CFTR and RHD, and used to train and validate a new prediction protocol named Splicing Prediction in Consensus Elements (SPiCE). SPiCE combines in silico predictions from SpliceSiteFinder-like and MaxEntScan and uses logistic regression to define optimal decision thresholds. It revealed an unprecedented sensitivity and specificity of 99.5 and 95.2%, respectively, and the impact on splicing was correctly predicted for 98.8% of variants. We therefore propose SPiCE as the new tool for predicting variant spliceogenicity. It could be easily implemented in any diagnostic laboratory as a routine decision making tool to help geneticists to face the deluge of variants in the next-generation sequencing era. SPiCE is accessible at (https://sourceforge.net/projects/spicev2-1/).

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.

Functional analysis of novel RHD variants: splicing disruption is likely to be a common mechanism of variant D phenotype.

BACKGROUND: We previously showed that several variations in the RHD gene, including synonymous changes, can be classified as splice site variants and may play a direct role in D variant phenotype expression. We sought to extend our study to additional candidates, notably in the first and last exons of the gene, by engineering a novel universal splice reporting vector, i.e., minigene. STUDY DESIGN AND METHODS: Our previous plasmid construct was modified to allow subcloning of any exon(s) of interest for assessing effect of variations on splicing. Seventeen novel and/or uncharacterized variations of the RHD gene were selected for the study and tested in our novel model. RESULTS: We engineered and validated a novel universal minigene for assessing virtually any variations of interest for splicing defect. Of the 17 variants tested in the novel model, 11 were shown to alter splicing either totally or partially, including the silent c.1065C>T variation, which induces major skipping of exon 7, and may therefore be responsible for reducing D antigen expression. We also showed that while all three missense variations c.1154G>C, c.1154G>T, and c.1154G>A in exon 9 are splice site variants, splicing is differentially altered and D-negative phenotype observed in the presence of the latter substitution is likely due to a defect in RhD protein folding. CONCLUSION: Overall, we hypothesize that splicing alteration is likely to be a common mechanism of D phenotype variation that has been underestimated so far. Further large-scale studies are necessary to demonstrate this statement definitely.

The SLC40A1 R178Q mutation is a recurrent cause of hemochromatosis and is associated with a novel pathogenic mechanism.

Hemochromatosis type 4 is one of the most common causes of primary iron overload, after HFE-related hemochromatosis. It is an autosomal dominant disorder, primarily due to missense mutations in SLC40A1 This gene encodes ferroportin 1 (FPN1), which is the sole iron export protein reported in mammals. Not all heterozygous missense mutations in SLC40A1 are disease-causing. Due to phenocopies and an increased demand for genetic testing, rare SLC40A1 variations are fortuitously observed in patients with a secondary cause of hyperferritinemia. Structure/function analysis is the most effective way of establishing causality when clinical and segregation data are lacking. It can also provide important insights into the mechanism of iron egress and FPN1 regulation by hepcidin. The present study aimed to determine the pathogenicity of the previously reported p.Arg178Gln variant. We present the biological, clinical, histological and radiological findings of 22 patients from six independent families of French, Belgian or Iraqi decent. Despite phenotypic variability, all patients with p.Arg178Gln had elevated serum ferritin concentrations and normal to low transferrin saturation levels. In vitro experiments demonstrated that the p.Arg178Gln mutant reduces the ability of FPN1 to export iron without causing protein mislocalization. Based on a comparative model of the 3D structure of human FPN1 in an outward facing conformation, we argue that p.Arg178 is part of an interaction network modulating the conformational changes required for iron transport. We conclude that p.Arg178Gln represents a new category of loss-of-function mutations and that the study of "gating residues" is necessary in order to fully understand the action mechanism of FPN1.

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.

Analysis of long-range interactions in primary human cells identifies cooperative CFTR regulatory elements.

A mechanism by which control DNA elements regulate transcription over large linear genomic distances is by achieving close physical proximity with genes, and looping of the intervening chromatin paths. Alterations of such regulatory 'chromatin looping' systems are likely to play a critical role in human genetic disease at large. Here, we studied the spatial organization of a ≈790 kb locus encompassing the cystic fibrosis transmembrane conductance regulator (CFTR) gene. Dysregulation of CFTR is responsible for cystic fibrosis, which is the most common lethal genetic disorder in Caucasian populations. CFTR is a relatively large gene of 189 kb with a rather complex tissue-specific and temporal expression profile. We used chromatin conformation at the CFTR locus to identify new DNA sequences that regulate its transcription. By comparing 5C chromatin interaction maps of the CFTR locus in expressing and non-expressing human primary cells, we identified several new contact points between the CFTR promoter and its surroundings, in addition to regions featuring previously described regulatory elements. We demonstrate that two of these novel interacting regions cooperatively increase CFTR expression, and suggest that the new enhancer elements located on either side of the gene are brought together through chromatin looping via CTCF.

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.

Comprehensive functional annotation of 18 missense mutations found in suspected hemochromatosis type 4 patients.

Hemochromatosis type 4 is a rare form of primary iron overload transmitted as an autosomal dominant trait caused by mutations in the gene encoding the iron transport protein ferroportin 1 (SLC40A1). SLC40A1 mutations fall into two functional categories (loss- versus gain-of-function) underlying two distinct clinical entities (hemochromatosis type 4A versus type 4B). However, the vast majority of SLC40A1 mutations are rare missense variations, with only a few showing strong evidence of causality. The present study reports the results of an integrated approach collecting genetic and phenotypic data from 44 suspected hemochromatosis type 4 patients, with comprehensive structural and functional annotations. Causality was demonstrated for 10 missense variants, showing a clear dichotomy between the two hemochromatosis type 4 subtypes. Two subgroups of loss-of-function mutations were distinguished: one impairing cell-surface expression and one altering only iron egress. Additionally, a new gain-of-function mutation was identified, and the degradation of ferroportin on hepcidin binding was shown to probably depend on the integrity of a large extracellular loop outside of the hepcidin-binding domain. Eight further missense variations, on the other hand, were shown to have no discernible effects at either protein or RNA level; these were found in apparently isolated patients and were associated with a less severe phenotype. The present findings illustrate the importance of combining in silico and biochemical approaches to fully distinguish pathogenic SLC40A1 mutations from benign variants. This has profound implications for patient management.