Insights into the role of glycerophospholipids on the iron export function of SLC40A1 and the molecular mechanisms of ferroportin disease.

SLC40A1 is the sole iron export protein reported in mammals. In humans, its dysfunction is responsible for ferroportin disease, an inborn error of iron metabolism transmitted as an autosomal dominant trait and observed in different ethnic groups. As a member of the major facilitator superfamily, SLC40A1 requires a series of conformational changes to enable iron translocation across the plasma membrane. The influence of lipids on protein stability and its conformational changes has been little investigated to date. Here, we combine molecular dynamics simulations of SLC40A1 embedded in membrane bilayers with experimental alanine scanning mutagenesis to analyze the specific role of glycerophospholipids. We identify four basic residues (Lys90, Arg365, Lys366, and Arg371) that are located at the membrane-cytosol interface and consistently interact with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) molecules. These residues surround a network of salt bridges and hydrogens bonds that play a critical role in stabilizing SLC40A1 in its basal outward-facing conformation. More deeply embedded in the plasma membrane, we identify Arg179 as a charged amino acid residue also tightly interacting with lipid polar heads. This results in a local deformation of the lipid bilayer. Interestingly, Arg179 is adjacent to Arg178, which forms a functionally important salt-bridge with Asp473 and is a recurrently associated with ferroportin disease when mutated to glutamine. We demonstrate that the two p.Arg178Gln and p.Arg179Thr missense variants have similar functional behaviors. These observations provide insights into the role of phospholipids in the formation/disruption of the SLC40A1 inner gate, and give a better understanding of the diversity of molecular mechanisms of ferroportin disease.

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.

Molecular model of the ferroportin intracellular gate and implications for the human iron transport cycle and hemochromatosis type 4A.

Ferroportin 1 (FPN1) is a major facilitator superfamily transporter that is essential for proper maintenance of human iron homeostasis at the systemic and cellular level. FPN1 dysfunction leads to the progressive accumulation of iron in reticuloendothelial cells, causing hemochromatosis type 4A (or ferroportin disease), an autosomal dominant disorder that displays large phenotypic heterogeneity. Although crystal structures have unveiled the outward- and inward-facing conformations of the bacterial homolog Bdellovibrio bacteriovorus Fpn (or Bd2019) and calcium has recently been identified as an essential cofactor, our molecular understanding of the iron transport mechanism remains incomplete. Here, we used a combination of molecular modeling, molecular dynamics simulations, and Ala site-directed mutagenesis, followed by complementary in vitro functional analyses, to explore the structural architecture of the human FPN1 intracellular gate. We reveal an interdomain network that involves 5 key amino acids and is likely very important for stability of the iron exporter facing the extracellular milieu. We also identify inter- and intradomain interactions that rely on the 2 Asp84 and Asn174 critical residues and do not exist in the bacterial homolog. These interactions are thought to play an important role in the modulation of conformational changes during the transport cycle. We interpret these results in the context of hemochromatosis type 4A, reinforcing the idea that different categories of loss-of-function mutations exist. Our findings provide an unprecedented view of the human FPN1 outward-facing structure and the particular function of the so-called "gating residues" in the mechanism of iron export.-Guellec, J., Elbahnsi, A., Le Tertre, M., Uguen, K., Gourlaouen, I., Férec, C., Ka, C., Callebaut, I., Le Gac, G. Molecular model of the ferroportin intracellular gate and implications for the human iron transport cycle and hemochromatosis type 4A.

TIMP1 intron 3 retention is a marker of colon cancer progression controlled by hnRNPA1.

We previously reported a 40-transcripts signature marking the normal mucosa to colorectal adenocarcinoma transition. Eight of these mRNAs also showed splicing alterations, including a specific intron 3 retention in tissue metalloprotease inhibitor I (TIMP1), which decreased during the early steps of colorectal cancer progression. To decipher the mechanism of intron 3 retention/splicing, we first searched for putative RNA binding protein binding sites onto the TIMP1 sequence. We identified potential serine arginine rich splicing factor 1 (SRSF1) and heterogeneous nuclear RiboNucleoProtein A1 (hnRNPA1) binding sites at the end of intron 3 and the beginning of exon 4, respectively. RNA immunoprecipitation showed that hnRNPA1, but not SRSF1 could bind to the corresponding region in TIMP1 pre-mRNA in live cells. Furthermore, using a TIMP1-based ex vivo minigene approach, together with a plasmon resonance in vitro RNA binding assay, we confirmed that hnRNPA1 could indeed bind to wild type TIMP1 exon 4 pre-mRNA and control TMP1 intron 3 splicing, the interaction being abolished in presence of a mutant sequence that disrupted this site. These results indicated that hnRNPA1, upon binding to TIMP1 exon 4, was a positive regulator of intron 3 splicing. We propose that this TIMP1-intron 3 + transcript belongs to the class of nuclear transcripts with "detained" introns, an abundant molecular class, including in cancer.

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.

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.

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.