article: Pain is necessary to alert us to actual or potential tissue damage. Specialized nerve cells in the body periphery, so called nociceptors, are fundamental to mediate pain perception and humans without pain perception are at permanent risk for injuries, burns and mutilations. Pain insensitivity can be caused by sensory neurodegeneration which is a hallmark of hereditary sensory and autonomic neuropathies (HSANs). Although mutations in several genes were previously associated with sensory neurodegeneration, the etiology of many cases remains unknown. Using next generation sequencing in patients with congenital loss of pain perception, we here identify bi-allelic mutations in the FLVCR1 (Feline Leukemia Virus subgroup C Receptor 1) gene, which encodes a broadly expressed heme exporter. Different FLVCR1 isoforms control the size of the cytosolic heme pool required to sustain metabolic activity of different cell types. Mutations in FLVCR1 have previously been linked to vision impairment and posterior column ataxia in humans, but not to HSAN. Using fibroblasts and lymphoblastoid cell lines from patients with sensory neurodegeneration, we here show that the FLVCR1-mutations reduce heme export activity, enhance oxidative stress and increase sensitivity to programmed cell death. Our data link heme metabolism to sensory neuron maintenance and suggest that intracellular heme overload causes early-onset degeneration of pain-sensing neurons in humans.
  Neurodegenerative disorders affecting peripheral sensory neurons lead to loss of pain perception as disease hallmark. The absence of protective behaviors towards noxious stimuli causes unintentional self-injuries and chronic ulcerations. Soft tissue infections and osteomyelitis, often requiring amputations, are common and complicate this disorder [1,2]. Autonomic dysfunction and motor deficits may be additional features of sensory and autonomic neuropathies (HSANs). Prominent loss of large and small myelinated fibers distinguishes sensory neuropathies from clinically similarly presenting channelopathy-associated pain insensitivity (CIP) [3]. Proteins which are involved in sensory neurodegeneration affect distinct molecular pathways: sphingolipid-metabolism, membrane-shaping of organelles, regulation of ion channels, endoplasmic reticulum turnover and axonal trafficking[1,4–8]. However, the molecular mechanisms underlying sensory neurodegeneration are still incompletely understood and disease-causing mutations remain to be identified in a substantial number of patients. Rapid progress in next-generation sequencing (NGS) technology has transformed the field of medical genomics leading to the identification of novel disease-genes[9,10]. In this study, next generation sequencing was performed in patients with HSAN but without mutations in the known genes associated with the disorder. Causative mutations were found in FLVCR1 (Feline Leukemia Virus subgroup C Receptor 1), a gene that has previously been associated to Posterior Column Ataxia and Retinitis Pigmentosa (PCARP) [11–16]. FLVCR1 is an ubiquitously expressed heme exporter[17,18], member of the Major Facilitator Superfamily (MFS) transporters[19]. Two different isoforms have been described. FLVCR1a resides in the plasma membrane and is responsible for heme detoxification in several cell types, such as erythroid progenitors, endothelial cells, hepatocytes, lymphocytes and intestinal cells[18,20–25]. FLVCR1b is located on mitochondria and is involved in the transport of newly synthesized heme from mitochondria to the cytosol[18]. The expression of FLVCR1a and FLVCR1b is needed to control the size of the cytoplasmic free-heme pool, which is essential for proper metabolic functions[21,23]. Heme is an essential co-factor involved in multiple biological processes: oxygen transport and storage, electron transfer, drug and steroid metabolism, signal transduction and microRNA processing[26]. However, excess free-heme is highly toxic due to its ability to promote oxidative stress, proteasome inhibition and mitochondrial dysfunction, that ultimately lead to cell death[26–28]. For this reason, the intracellular free-heme pool is finely regulated at multiple levels[26] and FLVCR1a-mediated heme export contributes to this process[21,22,26]. Our data show that primary fibroblasts and lymphoblastoid cell lines from patients with sensory neurodegeneration have reduced heme export activity due to FLVCR1-mutations. This results in enhanced oxidative stress and increased sensitivity to programmed cell death. These data add heme metabolism to the molecular pathways implicated in sensory neuron maintenance and pain processing.
  We studied trios with an affected child and healthy parents with the diagnosis of hereditary sensory and autonomic neuropathy (HSAN) (Fig 1 and Table 1). In an Italian boy with early-onset pain insensitivity (patient 1) and non-consanguineous parents whole-exome sequencing identified 559 variants with low prevalence (<0,01%) in dbSNP, the 1000-Genomes project, the Exome Variant Server or the ExAC browser. Variants were prioritized for de novo mutations and compound-heterozygous / homozygous variants (S1 Table). Amongst these variants compound-heterozygosity for mutations in FLVCR1 (NM_014053. 3) was detected (Fig 1B and 1C). The heterozygous mutation c. 574T>C; p. (Cys192Arg) has been previously reported as homozygous mutation in autosomal-recessive posterior column ataxia with retinitis pigmentosa (PCARP) [14]. The respective cytosine is located in the third transmembrane domain of the FLVCR1 protein (Fig 1B) and is evolutionarily conserved among species. The mother of patient 1 is heterozygous carrier of the c. 574T>C mutation. The second mutation was inherited from the father and is a heterozygous frameshift-mutation c. 610del; p. (Met204Cysfs*56). The deletion results in a premature stop codon, hence suggesting complete loss of protein function. Both mutations were considered as likely pathogenic accordingly to in-silico prediction programs (S2 Table). To forecast gene-specific pathogenicity of the variants in the context of the FLVCR1 gene, a Mutation Significance Cutoff (MSC) 23. 8 for CADD, 0. 697 for PolyPhen2 and 0. 04 for SIFT was applied based on calculations by the MSC server for HGMD and ClinVar variants[29]. If applicable, all tools predicted a deleterious impact of the specific mutations (S2 Table). Patient 1 was the only child of non-consanguineous parents and at 23 months of age presented with severe psychomotor delay, no language acquisition, absent reaction to painful stimuli, ulcero-mutilations of mouth and hands, joint hypermobility and scoliosis. Sensory nerve conduction was not evocable. The child had frequent crises resembling “dysautonomic crises”, characterized by agitation and tachycardia, spontaneously resolving after vomiting or stool emission. He had chronic macrocytic anemia, with normal ferritin, syderemia, transferrin, vitamin B12 and folic acid levels; bone marrow examination showed hyporegenerative anemia. At the age of 33 months he developed seizures. In the later disease course, optic atrophy and retinal degeneration were also noticed (Table 1). The clinical presentation thus showed a predominant pain insensitivity phenotype with clinical symptoms partially overlapping with PCARP, suggesting that distinct FLVCR1 mutations may result in different clinical outcomes. To substantiate the role of FLVCR1 in sensory neuropathies additional index cases and trios were subjected to next-generation sequencing using a gene panel consisting of 70 genes implicated in HSAN and similar syndromes, including FLVCR1. Amongst 165 index cases, we identified a second case with early-onset sensory neuropathy (patient 2) and compound heterozygous mutations c. 661C>T; p. (Pro221Ser) and c. 1324dup; p. (Tyr442Leufs*7) in FLVCR1 (Fig 1B and 1D). The non-consanguineous family of European origin (mother, father, index patient) was further subjected to trio exome-sequencing to evaluate mutations in other than the initially tested genes. No other candidate gene which was likely to account for the clinical phenotype was identified by this approach (S1 Table). Pro221 is an evolutionarily highly conserved amino acid residue of FLVCR1 and its mutation to serine is likely to interfere with FLVCR1 function. This missense exchange is extremely rare and accordingly to the Exome Aggregation Consortium (ExAC) has been observed in 1 of 120,990 alleles. Pro221Ser is also predicted as pathogenic (S2 Table). Again, the second FLVCR1 allele harbors a loss-of-function frameshift mutation, which has not been described until now. Patient 2 had mild developmental delay and started to walk with support at age of 24 months with a broad base gait. Speech was also delayed with first words at 22 months. At 12 months of age a non-healing wound at the tip of the right thumb with nail dystrophy was noticed. Despite local medication necrosis and acroosteolysis occurred. Similar lesions appeared on the 2nd and 3rd right fingers and left thumb with loss of nails and resorption of the underlying soft tissues. Wounds were slow healing, resulting in atrophic scarring. At examination, the patient presented hyperactive and clumsy, and showed marked reduction to painful stimuli. Affected fingers had stubby ends and nail dystrophy (Fig 1A). He had a positive histamine axonal flare response, and electroneurography showed a severe sensory neuropathy. Electroencephalogram and brain MRI gave normal results, while total spine MRI showed syrinx from T5 to T10, minor irregularities of the cervical central ependymal canal, and mild hyperintensity of the posterior columns (S1 Fig). An electroretinogram, performed at 3 years, showed reduced a-wave amplitude suggestive for initial retinal epithelium degeneration. Together, these data link FLVCR1 mutations to early-onset complicated sensory neuropathy. To investigate the functional consequences of FLVCR1 mutations, primary fibroblasts were derived from patient 1 whereas lymphoblastoid cell lines (LCLs) were generated from patient 2. qRT-PCR analyses showed that FLVCR1 mutations result in a specific decrease of FLVCR1a transcript whereas FLVCR1b mRNA levels were unaffected in both patients (Fig 2A and 2B and S2 Fig). Interestingly, immunoprecipitation and western blotting indicated that FLVCR1a is still expressed in patient fibroblasts and LCLs (Fig 2C and 2D). These data indicate that the mutations did not completely abolish FLVCR1a expression. Therefore, we hypothesized that the mutations could interfere with the heme export ability of FLVCR1, thus resulting in intracellular heme accumulation. To address this issue, we analyzed the expression of proteins involved in heme and iron metabolism. The expression of ALAS1 (δ-aminolevulinc acid synthase 1), the rate-limiting enzyme in the heme biosynthetic pathway, was similar in patient and control fibroblasts (Fig 2E and 2G). However, increased expression of the heme degrading enzyme HO1 (Heme oxygenase 1) was observed in patient compared to control fibroblasts (Fig 2E and 2G). HO1 catalyzes the degradation of heme into biliverdin, carbon monoxide and iron[30]. Consistent with the induction of HO1 in patient fibroblasts, increased mRNA levels of the iron exporter FPN (Ferroportin) and the iron-storage proteins FT-L (Ferritin L) and FT-H (Ferritin H) were observed in patient compared to control fibroblasts (Fig 2E). Contrary, patient LCLs were characterized by decreased ALAS1 and similar HO1 expression levels compared to healthy donors (Fig 2F, 2H and S3 Fig). Consistent with reduced ALAS1 expression, we observed a slight decrease of FPN, FT-L and FT-H mRNA levels in patient compared to control LCLs (Fig 2F). To definitively demonstrate that FLVCR1 mutations impair plasma membrane export, the amount of intracellular heme concentration was determined. Heme content was comparable between patient and control fibroblasts (Fig 2I) and LCLs (Fig 2J) under resting conditions, likely due to the induction of HO1 in patient fibroblasts and the reduction of ALAS1 in patient LCLs (Fig 2E–2H). However, following the stimulation of heme synthesis with ALA (δ-aminolevulinic acid), heme accumulation was observed in patient compared to control fibroblasts (Fig 2I). Moreover, heme content was higher in patient LCLs compared to the mean heme content of 4 different healthy donors LCLs (Fig 2J). Taken together these data indicate that FLVCR1 mutations impair plasma membrane heme export resulting in a transient accumulation of heme that activates cell-specific detoxifying mechanisms: excess-heme is rapidly catabolized by HO1 in patient-derived fibroblasts while heme downregulates its own synthesis in patient-derived LCLs. Excess free-heme is highly toxic due to its ability to promote oxidative stress and lipid peroxidation, leading to membrane injury and, ultimately, apoptosis[26]. To evaluate whether FLVCR1 mutations influence cellular oxidative status, the expression of genes involved in the antioxidant response was analyzed. Patient fibroblasts were characterized by decreased superoxide dismutase 1 (SOD1) and catalase mRNA levels, and increased superoxide dismutase 2 (SOD2) and thioredoxin transcripts (Fig 3A). Increased reactive oxygen species (ROS) were observed in patient compared to control fibroblasts (Fig 3B). Moreover, ROS levels were higher in patient LCLs compared to the mean ROS levels of 4 different healthy donors LCLs (Fig 3C). In order to investigate whether loss of FLVCR1a affects cell viability, the survival of cells was evaluated by Annexin V staining. The percentage of Annexin V-positive cells was comparable between patient and control fibroblasts under resting conditions (Fig 3D). However, following the stimulation with ALA, increased percentage of Annexin V-positive cells was detected in patient compared to control fibroblasts (Fig 3D). In addition, following ALA treatment, the percentage of Annexin V-positive cells was higher in patient LCLs compared to the mean of 4 different healthy donors LCLs (Fig 3E). The same result was observed following the stimulation of patient and controls LCLs with H2O2 (S4 Fig). In order to demonstrate that increased apoptosis in patient cells treated with ALA was due to heme overload, we performed rescue experiments by reducing intracellular heme loading. It has been previously reported that Hemopexin (HX) facilitates heme export through FLVCR1[31,32]. Thus, we added HX to the medium. HX treatment ameliorate cell survival in patient LCLs treated with ALA in a dose dependent manner (Fig 3F), thus suggesting that heme overload is responsible for the observed phenotype. Taken together these data indicate that the impairment of FLVCR1-mediated heme export triggers oxidative stress and increases the susceptibility to programmed cell death in patient-derived cells. The identification of FLVCR1 mutations in patients with loss of pain perception, as well as PCARP, indicate that FLVCR1 plays a critical role in the nervous system. Previous studies confirmed the expression of FLVCR1 in the mouse brain, retina and spinal cord[14]. We analyzed FLVCR1a mRNA levels in different mouse tissues. FLVCR1a mRNA levels are higher in the brain compared to other tissues characterized by elevated production of hemoproteins (S5 Fig). This data suggests that FLVCR1a mediated heme export is relevant for neurons. We previously reported that FLVCR1a is essential to prevent heme-induced oxidative stress in several cell types, including erythroid progenitors, hepatocytes and intestinal cells[21,22,33]. To get an insight into the role of FLVCR1a in neuronal cells, shRNA-mediated knockdown experiments were performed in SH-SY5Y neuroblastoma cells (Fig 4A and S2 Fig). Increased ROS levels were detected in FLVCR1a-downregulated SH-SY5Y cells compared to controls (Fig 4B), thus indicating a conserved role for FLVCR1a in neuroblastoma cells. Next, we investigated whether FLVCR1a is essential for the survival of these cells. Increased percentage of Annexin V-positive cells was detected in FLVCR1a-downregulated SH-SY5Y cells compared to controls, both under resting conditions and following the stimulation with ALA (Fig 4C). Together, these data indicate that heme export through FLVCR1a is important for the survival of neuronal cells.
  Different mechanisms leading to sensory neurodegeneration and pain loss have been described. For example mutations in voltage-gated sodium channels, membrane-shaping proteins or the autophagy receptor FAM134B, which triggers turnover of the endoplasmic reticulum, are causative for sensory neuron loss[1,5]. Altered sphingolipid metabolism, impaired neurotrophin signaling or epigenetic dysregulation likewise results in sensory neuron damage[34]. We here report patients with early-onset pain insensitivity because of mutations in the heme exporter FLVCR1. The finding suggests that heme metabolism is critically involved in sensory neuron maintenance. The FLVCR1 gene is under strong positive selection with Residual Variation Intolerance Score (RVIS) and Gene Damage Index (GDI) scores of 0. 02 (55. 45%) and 2854. 439, respectively. Thus, the genic intolerance is estimated to be rather average and when testing for selective pressure, the GDI is classified as “medium damage prediction”. As part of the Major-Facilitator-Superfamily (MFS) FLVCR1 and other family members share 12 transmembrane alpha helices connected by hydrophilic loops. The subRVIS-prediction tool indicates that the crucial MFS domain of FLVCR1, in which all four mutations reported here are located, has a low RVIS score of -1. 39 (3. 96%), indicating a strong negative selection for this particular domain. This is compatible with pathogenicity of mutations residing in the conserved MFS domain. In contrast, the less well conserved N- and C-terminal regions of FLVCR1 show higher subRVIS scores of 0. 07 (55. 86%) and 1. 12 (93. 99%), respectively, which explains the overall higher RVIS score of FLVCR1. FLVCR1 is a heme exporter involved in the control of intracellular heme levels in different cell types[20–25]. We demonstrated that sensory neuropathy-associated FLVCR1-mutations compromise the ability of FLVCR1a to export heme. Patient-derived primary cells showed altered expression of genes transcriptionally regulated by heme and accumulation of heme following the stimulation of endogenous heme synthesis. Thus, we propose that reduced heme export activity may cause heme accumulation and cellular damage by heme toxicity in sensory neurons. The analysis of patient-derived cells also indicates that FLVCR1 impairment triggers compensatory mechanisms to counteract the accumulation of heme. Indeed, increased heme degradation was observed in patient fibroblasts whereas decreased heme synthesis rates were found in patient LCLs. However, these mechanisms fail to completely protect the cells from heme toxicity, especially under stress conditions, resulting in oxidative stress and increased susceptibility to programmed cell death. These data are in agreement with previous studies indicating an essential role for FLVCR1a to prevent heme-induced oxidative stress in different cell types[21,22]. We here show a conserved role for FLVCR1a in neuronal cells. Silencing of FLVCR1a induces oxidative stress and affects the survival of neuroblastoma cells under resting conditions, thus suggesting that neuronal cells are particularly sensitive to the loss of FLVCR1a. We propose FLVCR1a as an important heme exporter in sensory neurons which regulates the oxidative stress response. Mutations in the human FLVCR1 gene were previously associated to Posterior Column Ataxia and Retinitis Pigmentosa (PCARP) [12]. PCARP is a childhood-onset, autosomal-recessive disorder with the clinical features of sensory ataxia and retinitis pigmentosa[11–14,35]. Interestingly, polyneuropathy has been reported in a patient affected by PCARP, however, not as leading clinical manifestation but as additional symptom at a later state of disease[11]. In the majority of patients affected by PCARP, homozygous missense mutations in the FLVCR1 gene have been reported[12]. These mutations affect highly conserved residues of FLVCR1 important for protein expression, localization, folding or regulation of heme efflux[12]. Both patients reported here harbour a missense-mutation on one allele and a severe frameshift-mutation on the other FLVCR1 alleles, the latter most likely resulting in a truncated FLVCR1a protein that may be rapidly degraded, acts as a dominant negative molecule or leads to nonsense-mediated mRNA decay (NMD). Such mutations are likely to cause severely impaired heme efflux rates, which may be particularly relevant for sensory neurons. This finding suggests that the type of mutation critically determines the phenotypic outcome of the disease. It is well known that heme is involved in multiple biological processes[26], yet the role of heme metabolism in the nervous system has been poorly addressed. The finding of FLVCR1 mutations in patients with peripheral sensory neuropathy suggests that heme metabolism via FLVCR1a plays an essential role for maintenance of these specialized neurons. Heme is essential for ATP production through oxidative phosphorylation. Indeed, heme is the fundamental cofactor of cytochromes of the electron transport chain[36–38] and reduced heme synthesis cause a reduction of complex IV activity and ATP synthesis[39]. Peripheral nerves need efficient energetic metabolism because of a considerable length of the axon. For this reason, sensory neurons may be particularly sensitive to disturbances in heme homeostasis and dependent on heme synthesis for oxidative phosphorylation, to sustain both anterograde and retrograde trafficking along the extraordinarily long axonal processes. Interestingly, the nerve growth factor signaling cascade (NGF/TrkA) is essential for sensory neuron survival and mutations in both NGF and its receptor, TrkA, encoded by NTRK1, lead to sensory and autonomic neuropathy. Since heme is known to regulate NGF signaling[40], it is tempting to speculate that altered NGF signaling provides a pathophysiological link to the FLVCR1 disorder. Heme levels are also directly involved in neuron excitability and nociception by modulating the activity of voltage-gated K+ channels. Voltage-gated K+ channels are a large family of K+ channels that open on membrane depolarization and contribute to maintain the resting potential and to determine the characteristics and frequency of action potentials[41]. Several studies highlight a pivotal role of K+ channels in pain processing[42,43]. Knocking-down specific K+ channels subunits results in alteration of nociception in mouse models[43]. Heme directly binds to the large-conductance calcium-dependent Slo1 BK channels and Kv1. 4 A-type K+ channels. The high affinity binding of heme to Slo1 decreases the frequency of channel opening leading to an inhibition of transmembrane K+ currents[44,45]. Its binding to the N-terminal domain of Kv1. 4 channels inhibits the fast inactivation of the channel, thus reducing cellular excitability[46]. Heme can also indirectly modulate the activity of K+ channels through the production of ROS and carbon monoxide. ROS modulation of K+ channel activity has already been reported[41] and CO plays multiple roles in nociception[47]. Moreover, it has been proposed that oxidative stress promotes the development of acquired forms of neuropathy, like diabetic neuropathy[48,49] and chemotherapy-induced neuropathy[50,51]. Taken together, our data emphasize a critical role of heme metabolism for homeostasis of sensory neurons and suggest FLVCR1 as a neuroprotective target.
  Patients were recruited by the specialized outpatient services for Rare Diseases at the Gaslini Hospital in Genova (Patient 1) and of Clinical/Medical Genetics at the San Camillo-Forlanini Hospital in Rome (Patient 2), Italy. Both patients are of European origin, and from non-consanguineous families. Patient 1 samples were obtained from the “Cell Line and DNA Biobank from patients affected by Genetic Diseases” (Istituto Giannina Gaslini), member of Telethon Network of Genetic Biobanks (project no. GTB12001). Patient 1 underwent peripheral blood sampling for molecular testing, and skin punch biopsy for fibroblast culture and functional studies. Patient 2 underwent peripheral blood sampling for molecular testing and lymphocyte culture. Because of slow healing wounds the family dissented with skin biopsies. Both patients were evaluated by full clinical and detailed instrumental investigations, comprising neurophysiological assessment, brain and total spine MRI. As controls, we used primary fibroblasts and LCLs derived from healthy donors. Control fibroblasts were obtained from the “Cell Line and DNA Biobank from patients affected by Genetic Diseases” (Istituto Giannina Gaslini), member of Telethon Network of Genetic Biobanks (project no. GTB12001). As control, we chose primary fibroblasts derived from an individual with comparable age to that of patient 1 (newborn child). Controls LCLs were derived from healthy donors at the Molecular Biotechnology Center, University of Torino, Italy. Control and patient samples derived from individuals belonging to the same genetic population (European origin). The study was approved by the local ethics board of the University Hospital Jena (research ethics approval 3244-09/11) and of the San Camillo Forlanini Hospital Rome (research ethics approval OSR83/2016). Written consent for the study was obtained from the participants or their legally authorized representatives. This research is in accordance with the Helsinki Declaration. Control and patient-derived cells, as well as SH-SY5Y cells, were collected, washed in PBS1X, resuspended in 10 mM Hepes, 150 mM NaCl, 5 mM CaCl2 buffer, and labelled with AnnexinV-FITC (BD Biosciences) for 15 minutes. Then, 2 μl propidium iodide (1 mg/ml) (PI; Sigma-Aldrich) was added. AnnexinV emission was detected in the green channel (525 nm) and propidium iodide in the red channel (575nm) on a FACS Calibur Cytometer (BD Biosciences). For each sample, 10’000 to 25’000 events were collected. Results were expressed as mean ± SEM. Statistical analyses were performed using one-way or two-way analysis of variance or Student’s t test. A P value of less than 0. 05 was considered significant.
summary: Hereditary Sensory and Autonomic Neuropathy (HSAN) is a genetic disorder mainly characterized by the impairment of sensory neurons, which transmit information about sensations such as pain, temperature and touch. Therefore, unintentional self-injury, leading to ulcers and eventually amputations are common in affected individuals. Although mutations in several genes were previously associated with sensory neurodegeneration and pain insensitivity, the etiology of many cases remains unknown. We here identify mutations in the heme exporter protein FLVCR1 in patients with congenital inability to experience pain. We showed that FLVCR1 mutations results in reduced heme export activity, enhanced oxidative stress and increased sensitivity to programmed cell death. These data assign a surprising role for heme to sensory neuron maintenance.
keywords: heme
  cell death
  somatosensory system
  medicine and health sciences
  oxidative stress
  cell processes
  neuroscience
  fibroblasts
  mutation
  immunoprecipitation
  connective tissue cells
  frameshift mutation
  research and analysis methods
  sensory physiology
  animal cells
  proteins
  connective tissue
  biological tissue
  precipitation techniques
  biochemistry
  sensory neurons
  cellular neuroscience
  pain sensation
  cell biology
  post-translational modification
  anatomy
  physiology
  neurons
  genetics
  apoptosis
  biology and life sciences
  cellular types
  sensory systems
title: Mutations in the Heme Exporter FLVCR1 Cause Sensory Neurodegeneration with Loss of Pain Perception