# Characterization of XPR1/SLC53A1 variants located outside of the SPX domain in patients with primary familial brain calcification

Uriel López-Sánchez, Gaël Nicolas, Anne-Claire Richard, David Maltête,

Mahmoud Charif, Xavier Ayrignac, Cyril Goizet, Jawida Touhami, Gilles

Labesse, Jean-Luc Battini, et al.

► To cite this version:

Uriel López-Sánchez, Gaël Nicolas, Anne-Claire Richard, David Maltête, Mahmoud Charif, et al. Characterization of XPR1/SLC53A1 variants located outside of the SPX domain in patients with primary familial brain calcification. Scientific Reports, In press, 10.1038/s41598-019-43255-x. http://02109522

HAL Id: hal-02109522

https://hal.science/hal-02109522

Submitted on 3 Mar 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# Characterization of XPR1/SLC53A1 variants located outside of the SPX domain in patients with primary familial brain calcification

Received: 30 November 2018

Accepted: 17 April 2019

Published online: 01 May 2019

Uriel López-Sánchez<fnref n="1" />, Gaël Nicolas², Anne-Claire Richard², David Maltête³,⁴, Mahmoud Charif⁵, Xavier Ayrignac⁵, Cyril Goizet⁶, Jawida Touhami<fnref n="1" />, Gilles Labesse⁷, Jean-Luc Battini<fnref n="1" />,⁸ & Marc Sitbon<fnref n="1" />

Primary familial brain calcification (PFBC) is a rare neurological disease characterized by deposits of calcium phosphate in the basal ganglia and other regions of the brain. Pathogenic variants in the *XPR1/SLC53A1* gene, which encodes the only known inorganic phosphate exporter, cause an autosomal dominant form of PFBC. These variants are typically located in the SPX N-terminal domain of the protein. Here, we characterize three XPR1 variants outside of SPX in three PFBC patients with an apparently sporadic presentation: c.1375C > T p.(R459C), c.1855A > G p.(N619D) and c.1886T > G p.(I629S), with the latter identified as the first *XPR1/SLC53A1* de novo mutation to occur in a PFBC proband. When tested in an *in vitro* physiological complementation assay, the three XPR1 variants were impaired in phosphate export function, although they were normally expressed at the cell surface and could serve as functional receptors for retrovirus entry. Moreover, peripheral blood cells from the p.N619D patient could be assayed *ex vivo* and displayed significantly impaired phosphate export. Our results establish for the first time the clinical and molecular characteristics of XPR1 variants located outside the SPX domain and assert a direct link between these variants, deficient phosphate export, and PFBC. Moreover, we unveiled new structural features in XPR1 C-terminal domain that play a role in phosphate export and disease.

Primary familial brain calcification (PFBC) is a rare neurological disease characterized by the presence of calcium phosphate deposits in the microvessels of the basal ganglia and other brain regions. Clinical signs may start at any age<fnref n="1" /> and comprise mostly movement disorders, cognitive impairment and psychiatric symptoms, and other neurological symptoms². However, up to 42% of the patients can be asymptomatic<fnref n="1" />. PFBC is typically inherited as an autosomal dominant trait with four causal genes identified: *PiT2/SLC20A2*³, *PDGFRB*⁴, *PDGFB*⁵ and *XPR1*⁶, which has recently been assigned to the solute carrier (SLC) family of transporters as *SLC53A1* (https://www.genenames.org/cgi-bin/genefamilies/set/752/). Most recently, *MYORG* bi-allelic mutations have also been identified as causing autosomal recessive PFBC⁷.

Interestingly, both *PiT2/SLC20A2*⁸,⁹ and *XPR1/SLC53A1* encode phosphate transporters, with the latter encoding for XPR1, the sole known inorganic phosphate exporter in humans and other metazoans¹⁰. Five *XPR1/SLC53A1* missense variants have been identified to be causal for PFBC in seven families⁶,¹¹ and functional analyses have shown that the five XPR1 corresponding proteins have reduced cell surface expression and/or were

<footnote n="1">¹Institut de Génétique Moléculaire de Montpellier, University of Montpellier, CNRS, Montpellier, France. ²Normandie Univ, UNIROUEN, Inserm U1245, Rouen University Hospital, Department of Genetics and CNR-MAJ, Normandy Center for Genomic and Personalized Medicine, Rouen, France. ³Department of Neurology, Rouen University Hospital and University of Rouen, Rouen, France. ⁴INSERM U1239, Laboratory of Neuronal and Neuroendocrine Differentiation and Communication, Mont-Saint-Aignan, France. ⁵Department of Neurology, Montpellier University Hospital, Montpellier, France. ⁶INSERM U1211, Univ Bordeaux, Laboratoire Maladies Rares, Génétique et Métabolisme; CHU Bordeaux, Service de Génétique Médicale, Bordeaux, France. ⁷Centre de Biochimie Structurale, University of Montpellier, CNRS, Montpellier, France. ⁸Institut de Recherche en Infectiologie de Montpellier, University of Montpellier, CNRS, Montpellier, France. Correspondence and requests for materials should be addressed to J.-L.B. (email: jean-luc.battini@irim.cnrs.fr) or M.S. (email: marc.sitbon@igmm.cnrs.fr)</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Pedigrees and CT-scans of patients harboring *XPR1/SLC53A1* variants located outside of the SPX-encoding domain. (a) Reduced pedigrees of the probands harboring *XPR1/SLC53A1* variants. Open symbol: unaffected; filled symbol: affected (presence of brain calcifications); symbol with a diagonal line: deceased; arrow: proband; numbers in parentheses: age at death; Star: DNA available; WT: wild type. (b) CT scans of the probands in axial sections, arrows indicate the different areas of calcifications; Le: Lenticular; T: Thalamic; SCWM: Subcortical white matter; CbH: Cerebellar hemisphere; Ve: Vermis; Midb: Midbrain. (c) Schematic representation of the *XPR1/SLC53A1* gene with the SPX-encoding exons indicated in orange, with the changes corresponding to the three PFBC variants indicated on exons 11 and 14 (arrows). (d) Schematic representation of the XPR1 protein structure with the SPX domain highlighted in orange, and the three residues corresponding to the two XPR1/SLC53A1 variants p.(R459C and p.N619D) and the de novo mutant (p.I629S) shown in green. The topological structure was generated using Protter program³⁹. (e) Alignment of XPR1 orthologs, exclusively from metazoan species, showing the evolutionary conservation of the regions wherein PFBC variants were located.</figcaption></figure>

significantly impaired in phosphate export⁶,¹¹. Interestingly, all five damaging variants identified and functionally characterized so far are located in the N-terminal cytoplasmic region of XPR1, with four located within the SPX domain⁶,¹¹. SPX, which is highly conserved among eukaryote XPR1, has been shown to regulate phosphate intracellular levels in yeast and plants¹². Here, we characterized three novel *XPR1/SLC53A1* missense variants that are located outside of the N-terminal SPX-encoding region, of which two were located in the XPR1 cytoplasmic C-terminal encoding region. We show that the c.1855A > G p.(N619D) C-terminal variant, which has been previously reported as of unknown significance¹³, can now be considered as causal for PFBC and used in a clinical setting. Moreover, DNA sequencing of both parents of the proband allowed us to confirm parenthood and identify the c.1886T > G p.(I629S) C-terminal variant, reported here for the first time, as the first PFBC mutation documented to occur *de novo* in *XPR1/SLC53A1*. Furthermore, we assessed the functional impact of the corresponding mutations on cell surface expression, retroviral receptor function and phosphate export, unveiling for the first time a key and specific role of domains outside the N-terminal regions and the SPX domain in XPR1 phosphate export functions, and PFBC.

## Results

**XPR1/SLC53A1 variants outside of the region encoding the N-terminal cytoplasmic domain.**

We report here the phenotypes associated with three *XPR1/SLC53A1* variants. Two were recorded by the French PFBC study group: c.1375C > T p.(R459C) and c.1855A > G p.(N619D); and the third one is a novel c.1886T > G p.(I629S) variant (Fig. 1a,b). The two former variants were mentioned in a recent international study aiming at reporting diagnostic yields¹³, while the latter, identified after the data freeze of the aforementioned study, is entirely novel.

All three variants were heterozygous (Fig. S1). The p.R459C variant maps to exon 11 (NG_050964, NM_004736.3), corresponding to the predicted largest extracellular loop of XPR1/SLC53A1, while variants p.N619D and p.I629S, encoded by exon 14, are both located in the XPR1/SLC53A1 cytoplasmic C-terminal domain (Fig. 1c,d). No particular function for XPR1/SLC53A1 C-term has been described yet. Remarkably, all
three variants are located in highly evolutionary conserved regions of XPR1 (Fig. 1e), and are predicted damaging by the MutationTaster¹⁴, PolyPhen-2¹⁵, and SIFT¹⁶ bioinformatics tools. All three variants are absent from control databases, including the genome aggregation database (gnomAD) that gathers over 130,000 individuals (http://gnomad.broadinstitute.org/)¹⁷.

**Phenotype of the patients harboring the new XPR1/SLC53A1 variants.** *Patient ROU 5059, c.1375C > T p.(R459C) variant.* The patient presented at age 58 years with parkinsonism that started at age 55 by a left upper limb tremor. He was born preterm and presented learning difficulties requiring education in a school for special needs. A diagnosis of mild intellectual disability was performed. In addition, he was treated for high blood pressure, presented a history of arthrosis from age 15 and experienced idiopathic peripheral facial palsy with partial spontaneous recovery at age 58. Upon examination at age 59, he presented bilateral akinetic hypertonic syndrome with left predominance associated to a resting tremor. This parkinsonism was considered as L-Dopa-responsive.

His computed tomography (CT) scan revealed severe bilateral lenticular calcifications as well as moderate bilateral cerebellar hemisphere, bilateral faint supratentorial white matter and faint left midbrain calcifications, with a total calcification score (TCS) of 20/80 (Fig. 1). Brain scintigraphy with ¹²³I-Iofupane showed marked bilateral dopaminergic neuron loss with right predominance. The patient had a positive family history of Parkinson disease for the maternal grandfather (no CT scan), while his mother was asymptomatic and presented a normal CT scan. His father died at the age of 67, with no history of neuropsychiatric disease.

We identified a c.1375C > T p.(R459C) variant in the proband, which was absent in his mother's DNA. This variant was classified as likely pathogenic (class 4) following the ACMG-AMP recommendations¹³.

*Patient EXT 1219, c.1855A > G p.(N619D) variant.* The patient presented with sudden deafness of unknown cause and unstable gait at the age of 69. He had a negative personal medical history. Upon examination, he presented a mild static cerebellar syndrome. The patient's CT scan revealed bilateral severe lenticular calcifications as well as bilateral moderate supratentorial white matter, bilateral faint calcification of the thalamus, bilateral severe cerebellar calcification, and faint vermian calcification (TCS = 29/80).

The patient's 71-year-old brother was reported to present a gait disorder since he was 50. However, neither CT scan nor DNA samples were available for the brother. The father died at the age of 50 after committing suicide.

We identified a c.1855A > G p.(N619D) *XPR1/SLC53A1* variant in this patient. This variant has been initially recorded as of unknown significance (class 3), following the ACMG-AMP recommendations¹³.

*Patient ROU 5075, c.1886T > G p.(I629S) variant.* The patient was 42 years of age when she presented with fluc tuant paresthesia of the left upper limb. Her medical history was marked by a diagnosis of celiac disease. She also had a history of migraine without aura. She presented several episodes of paresthesia, sometimes associated with light paresis of the limb or ipsilateral hemiparesis that were not related to her headache episodes. She also com plained of overall fatigability and anxiety. Upon examination at age 45, she still presented subjective hypoesthesia of the left hemibody together with left motor fatigability. However, no objective sign was observed, including the absence of pyramidal, extrapyramidal, or cerebellar sign. Brain CT scan of the patient revealed the presence of bilateral moderate lenticular and cerebellar hemisphere calcifications (TCS = 12/80) (Fig. 1).

The patient's family history was negative. Both parents presented a normal neurological examination with normal CT scans. She presented though a heterozygous c.1886T > G p.(I629S) variant that we could demonstrate to have occurred de novo, as neither parent harbored this variant (Fig. S1). We confirmed parenthood using informative microsatellites. Hence, the p.I629S variant could be classified as likely pathogenic (class 4), following ACMG-AMP recommendations.

**All XPR1/SLC53A1 PFBC variants are functionally expressed at the cell surface.** To evaluate the impact of *XPR1/SLC53A1* variants on XPR1 expression, we cloned and assayed *WT XPR1*, as well as mutants encoding the p.R459C, p.N619D or p.I629S variants, or the artificial XPR1 deletion construct, p.(L612_T696del), in which the C-terminal cytoplasmic domain was entirely deleted. All constructs were cloned in the pCHIX expression vector. Expression, presence at the cell surface, and retroviral receptor and transport functions were assessed for the PFBC variants, in comparison with those of the WT XPR1 and the p.(L612_T696del) construct.

WT and all three PFBC variants of XPR1/SLC53A1 were detected at the cell surface by flow cytometry using the XRBD ligand, an XPR1-specific ligand previously described¹⁰. In contrast, the p.(L612_T696del) construct was not detectable at the cell surface (Fig. 2a) despite efficient expression. While all naturally occurring variants of XPR1/SLC53A1 were equivalently and efficiently expressed at the cell surface, they appeared to accumulate at different intracellular levels (Fig. 2b). Also, expression of the PFBC variants had no detectable effect on the cell surface expression of the phosphate importers PiT1/SLC20A1 and PiT2/SLC20A2 (Fig. 2a).

XPR1/SLC53A1 is the entry receptor of mouse retroviruses, in particular the xenotropic murine leukemia viruses (X-MLV)¹⁸⁻²⁰. We therefore monitored PFBC XPR1/SLC53A1 variants for their capacity to serve as cell surface receptor for X-MLV entry and infection. For this purpose, we generated CHO hamster cells, which are refractory to X-MLV entry, as stably expressing either WT human XPR1 or the different PFBC variants. Similar cell surface expression of the WT XPR1 and the three PFBC variants was also observed on hamster CHO cells, while the p.(L612_T696del) construct was not detected at the cell surface (Fig. 3a, upper panel). Of note, the p.(L612_696Tdel) construct reached a significantly lower overall cellular expression in CHO cells when compared to WT and the three PFBC variants (Fig. 3b). All cell lines showed equivalent levels of expression of the glucose transporter GLUT1/SLC2A1 at the cell surface, demonstrating no generic defects or differences in SLC cell surface expression among the different CHO cell lines (Fig. 3a, lower panel). When challenged for retrovirus infection mediated by the X-MLV envelope glycoprotein (Env), we observed that PFBC variants served as
<figure><img src="image_4.png" /><figcaption>**Figure 2.** PFBC *XPR1/SLC53A1* variants and XPR1, PiT1 and PiT2 cell surface expression. (a) Cell surface detection of XPR1/SLC53A1, PiT1/SLC20A1 and PiT2/SLC20A2 with RBD specific ligands derived from mouse xenotropic MLV, koala endogenous retrovirus, and mouse amphotropic-MLV, respectively (open histograms); HEK293T positive and negative control cells for XPR1/SLC53A1 expression were transfected with either an anti-*luciferase* siRNA (siLUC) or an anti-*XPR1/SLC53A1* siRNA (siXPR1), respectively. WT XPR1/SLC53A1 and variants were assayed for cell surface expression upon co-transfection of siXPR1 alone (−) or in combination with an expression vector coding either for the HA-tagged WT XPR1 (WT), or the PFBC XPR1 C-terminal variants (R459C, N619D, or I629S), or the p.(L612_696Tdel) artificial XPR1 mutant (Δ612–696). Numbers indicate the delta mean fluorescence intensity of a representative experiment (n = 3) compared to the non-specific staining obtained with the secondary IgG antibody (grey histograms). (b) Representative immunoblot of HA-tagged XPR1/SLC53A1-containing cell lysates with an anti-HA antibody (upper panel), or with an anti β-actin antibody used as loading control (lower panel).</figcaption></figure>

cell surface retroviral receptor as efficiently as WT XPR1 (Fig. 3c), rendering CHO cells as efficiently infectable as H293T susceptible cells (Fig. 3e). As expected, cells expressing the p.(L612_696Tdel) XPR1 construct that was not expressed at the cell surface, remained uninfected (Fig. 3c). All CHO cell lines, including the latter, remained equally permissive to retroviral infection mediated by the VSV-G Env which enter cells through a different receptor, indicating that all cell lines had similar abilities to support general retroviral infection past viral entry (Fig. 3d), and that defects in infection by the p.(L612_696Tdel) construct was solely due to the absence of this artificial mutant at the cell surface.

**Deficiency in phosphate export function of the three novel XPR1/SLC53A1 PFBC variants.**

XPR1 is a conserved and ubiquitously expressed SLC, which we previously identified as the only known metazoan inorganic phosphate exporter¹⁰. Phosphate export mediated by XPR1 can be efficiently assessed by a complementation assay, in which the knockdown of endogenous XPR1 by an siRNA directed against *XPR1 3'UTR* and the ensuing decreased phosphate efflux are compensated by the *de novo* transfection of an XPR1 expression vector that is not recognized by the siRNA¹⁰. Upon introduction of WT XPR1, siRNA-mediated decreased phosphate efflux was fully complemented, while complementation by the p.(L612_696Tdel) construct, which was not expressed at the cell surface, was similar to that of an empty control vector (Fig. 4a, left panel). Noticeably, the three PFBC variants had different levels of complementation for the phosphate exporter function, with the p.N619D variant showing the most marked defect (Fig. 4a, left panel). Importantly, this defect was observed only for phosphate efflux, as phosphate uptake remained unchanged for all constructions, thus confirming the specificity of XPR1 for phosphate export (Fig. 4a right panel).

We had access to peripheral blood mononuclear cells (PBMC) from the patient harboring the p.N619D variant. When assayed *ex vivo* for phosphate efflux, PBMC from the p.N619D patient were significantly impaired in phosphate efflux (Fig. 4b left panel), while phosphate uptake remained unaffected in the patient PBMC (Fig. 4b right panel).
<figure><img src="image_5.png" /><figcaption>**Figure 3.** PFBC XPR1/SLC53A1 variants do not alter virus entry. (a) Cell surface detection of XPR1/SLC53A1 (upper panel) and GLUT1/SLC2A1 (lower panel) in CHO hamster cell populations stably expressing the control parental expression vector lacking XPR1/SLC53A1 (LNCX), or encoding either WT XPR1/SLC53A1, or variants p.R459C, p.N619D or p.I629S, or the p.(L612_\_T696del) deletion mutant. Cells were assayed for XPR1/SLC53A1 and GLUT1/SLC2A1 surface expression with their corresponding XRBD and H2RBD ligands, respectively. Numbers indicate the specific delta mean fluorescence intensity of a representative experiment (n = 3), as compared to non-specific staining with the secondary IgG (grey histograms). (b) Representative immunoblot of HA-tagged XPR1/SLC53A1-containing cell lysates from (a) with an anti-HA antibody (upper panel), or with an anti β-actin antibody used as loading control (lower panel). (c) Analysis of viral infection in CHO cells stably expressing either of WT or XPR1/SLC53A1 variants. Cells were assayed for infection with retroviral vectors pseudotyped with X-MLV Env or (d) VSV-G. Two days after infection, cells were analyzed by flow cytometry. Graphics are average values from 3 different wells in one representative experiment (n = 3). (e) Infectivity of viral supernatants on control HEK293T susceptible cells.</figcaption></figure>

Discussion

All XPR1 variants demonstrated so far to be causative of PFBC have been mapped within the N-terminal cyto-
plasmic domain, mostly in the highly conserved SPX domain. Here, we described and characterized for the
first time three PFBC-associated variants outside of SPX, one in extracellular loop 3 and two in the C-terminal
cytoplasmic region. Moreover, with the *XPR1/SLC53A1* p.(I629S) variant, we report the first mutation that
occurred *de novo*. The only non equivocal *de novo* mutations previously identified in PFBC patients were found
in PiT2/SLC20A2<sup>21</sup> and *PDGFB*<sup>22</sup>. Moreover, we showed that all three XPR1 variants studied here have impaired
phosphate exporter function, while maintaining efficient cell surface expression and retroviral receptor function.
Interestingly, a more profound defect in Pi export as tested *in vitro* was consistently observed with the two vari-
ants located in the cytoplasmic region.

The three patients reported here showed lenticular and cerebellar hemisphere calcifications that appeared to
be less severe than those observed in patients with XPR1 changes located in the SPX-harboring N-terminal cyto-
plasmic region. However, the number of patients remains too low to perform subgroup analyses that would take
into account age as a covariable. Nevertheless, our present characterization and assessment of these three XPR1
variants, in exons 11 and 14, allowed the reclassification of c.1375C > T p.(R459C) and c.1886T > G p.(I629S) as
pathogenic variants (class 5), with an impaired phosphate exporter function likely to play a direct role in patho-
genicity. More importantly, our evaluation of variant c.1855A > G p.(N619D), which has been considered so far
as of unknown significance, allows now its classification as likely pathogenic (class 4), and therefore its use for
genetic counseling.

Depletion of *Xpr1* in mice is lethal, emphasizing the key role of this gene in phosphate homeostasis<sup>23</sup> and we
have previously shown that mutations in the cytoplasmic N-terminal domain of XPR1 associated to PFBC lead
to a defective phosphate export *in vitro* and *ex vivo*<sup>6,11</sup>. Also, a *Pit2/Slc20a2* conditional knockout mouse model
led to dysregulation of phosphate levels in cerebrospinal fluid and to brain vascular calcifications<sup>24-26</sup>. Although
calcium phosphate depositions in PFBC appeared to be due to an alteration of phosphate metabolism, the
mechanisms underlying the respective roles of the five genes identified as causative for PFBC remain unsolved.
Interestingly, heterozygous variants in *PiT2/SLC20A2*, *PDGFRB*, *PDGFB*, and *XPR1/SLC53A1* are associated with
disease development, while only bi-allelic variants in *MYORG* have been described in PFBC<sup>7</sup>.

The SPX domain has been shown to bind inositol pyrophosphate molecules with high affinity¹². Whether XPR1 C-terminal region also interacts with cytosolic partners or with SPX remains to be studied. Interestingly, XPR1 C-terminal region includes a conserved segment from residues 600 to 635 (numbering is using human XPR1 as a reference), with the cytoplasmic domain predicted to start at L612 (Figs 1d and S2). Conserved region
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Effect of PFBC mutations on XPR1/SLC53A1-mediated phosphate efflux. (a) Efflux and uptake of [³³P]-phosphate in HEK293T cells transfected with siRNA directed against *luciferase* (siLUC), or *XPR1/SLC53A1* (siXPR1), or with the combination of siXPR1 and an expression vector coding either for the HA-tagged WT XPR1/SLC53A1 (WT), or the R459C, N619D, or I629S PFBC variants, or the artificial XPR1/SLC53A1 mutant Δ612–696. Results are shown as means ± s.e.m. in a representative experiment (n = 3); \*P ≤ 0.05; \*\*P ≤ 0.001; \*\*\*P ≤ 0.0001, in comparison to the control (siLUC). (b) Efflux and uptake of [³³P] phosphate in PBMC isolated either from a healthy donor (circles) or from the PFBC patient harboring the p.N619D variant (squares). Efflux is measured as % of released to total cellular [³³P] phosphate; and uptake is expressed as nmol/mg of protein in 30 min. Bars represent means ± s.e.m \*\*P ≤ 0.01.</figcaption></figure>

600–635 is predicted to be mainly in helical conformation (PSIPRED²⁷, DISOPRED²⁸, MEMSAT²⁹), while the very C-terminus sequence, downstream of residue M646, is highly variable and predicted to be natively unfolded (not shown). Isoleucine residue at position 629 is strictly conserved in metazoan (Fig. 1e) and often substituted in other species to the similar and hydrophobic valine (Fig. S2). Conservation of a hydrophobic residue at this position supports that a substitution toward a polar and flexible residue, such as seen with the PFBC p.1629S variant, is damaging. With regard to the PFBC p.N619D variant, it is remarkable that no substitution of N619 toward D619 was observed across species, although Asp residue is the Asn most similar and isosteric residue, mainly differing by a negative charge harbored by Asp. It is therefore likely that the additional charge brought by Asp in PFBC p.N619D prevents proper local folding or interaction with other domains of the protein core or (a) cytoplasmic partner(s). As tested *ex vivo* on the proband PBMC, the p.N619D variant led to a defect of phosphate export in the proband's cells (Fig. 4b).

Interestingly, the p.(L612_T696del) construct in which the C-terminal cytoplasmic domain has been deleted can be expressed, although it did not appear to reach the cell surface. A similar role of the C-terminal cytoplasmic domain in the plasma membrane trafficking of GLUT1/SLC2A1 has been observed³⁰; unpublished observations. On the other hand, this C-terminal cytoplasmic domain also harbors multiple predicted phosphorylation and ubiquitin consensus sites between positions 660 and 690 that may explain trafficking defect. Our data however suggest that partner interaction of the XPR1 C-terminal domain plays a role beyond vesicular trafficking. Therefore, the search for XPR1/SLC53A1 C-terminal interacting partners is likely to further our understanding of the role of this region in phosphate homeostasis and calcification diseases.

The three XPR1 C-terminal variants analyzed here were efficiently expressed, despite that XPR1 p.N619D showed lower level expression. Moreover, when introduced into an XPR1 expression vector, the three mutations did not alter the capability of XPR1 to serve as retroviral receptor for X-MLV. Therefore, the XPR1 variants reported here seem to exclusively alter phosphate export, although at different levels, with impact neither on expression of the phosphate importers PiT1/SLC20A1 and PiT2/SLC20A2, nor on phosphate uptake. In this sense, the XPR1 C-terminal variants are similar to the PFBC variants of the SPX domain⁶,¹¹. XPR1/SLC53A1 heterozygous variants are not the most common cause of PFBC²; however, all the functional analyses of damaging variants resulted in phosphate efflux impairment. Importantly, this study is the first one to report a *de novo* XPR1 mutation causing PFBC p.(I629S). Moreover, we characterized clinically and functionally for the first time XPR1/SLC53A1 variants located outside of the region that encodes the SPX-harboring N-terminal cytoplasmic domain, unveiling a specific role of the XPR1 C-terminal domain in phosphate transport and PFBC. The
nature and evolutionary residue conservation of the two variants of the C-terminal domain, p.N619D and p.I629S (Figs 1e and S2), as well as their specific impact on phosphate export, suggested that the former mutation likely impacted on a direct protein interaction, either with cellular partner(s) or intramolecularly, whereas the latter more likely leads to changes in a local structural conformation.

## Patients, Materials and Methods

**Guidelines and regulations.** All patients or legal guardians, and unaffected relatives when involved, provided informed written consent for genetic analyses and experimental protocols that were approved by the Comité de Protection des Personnes (CPP) Ile de France II ethics committee. All experiments with human and animal samples and cell lines were performed in accordance with guidelines and regulations set by the French Centre National de la Recherche Scientifique (CNRS), and the Institut National de la Santé et la Recherche Médicale (INSERM).

**Patients.** Patients were recruited by the French PFBC study group from multiple French centers as previously described³¹. Medical charts and DNA blood samples were sent to Inserm U1245, Rouen, France, for genetic analyses, with the aforementioned ethics committee approval. XPR1 screening was performed by Sanger sequencing in patients formerly negatively screened for SLC20A2 and PDGFB variants. Variants were interpreted following the ACMG-AMP recommendations³².

After the identification of the I629S variant in patient ROU 5075, which was absent from the parental sample, we confirmed parenthood using a set of four informative microsatellites.

CT scans were analyzed using a previously described method that allowed the rating and computing of patient calcifications as a total calcification score (TCS) on a scale of 0–80³¹.

**Cells.** Human HEK293T and Chinese hamster ovary (CHO) cells were cultured in DMEM supplemented with 10% fetal bovine serum (FBS) (PAN-Biotech) and non-essential amino acids. Cells were incubated at 37 °C in a 5% CO2 and humid atmosphere. For phosphate-free experiments, cells were cultured in phosphate-free DMEM supplemented with 10% dialyzed FBS, as previously described⁶,¹⁰,¹¹.

**Plasmids and siRNAs.** The p.R459C, p.N619D and p.I629S mutations were generated by site-direct mutagenesis using recombinant overlapping PCR and introduced in the pCHIX expression vector, which contains two copies of the HA tag³³. HA-tagged versions of the XPR1 variants were also inserted into a pLNCX retroviral vector³⁴. Small interfering RNAs (siRNA) (Integrated DNA Technologies) targeting the 3' UTR of human XPR1/SLC53A1 were as follows: 5'-ugauauAACuccugugcaattt-3' and 5'-ccuuaacagcuguaagcuautt-3'. siRNA directed against the firefly luciferase gene was used as control. HEK293T cells grown on poly-D-lysine-coated 6-well plates were transfected with 50 pmol siRNA per well, using the calcium phosphate method along with either empty or XPR1/SLC53A1-harboring expression vectors.

**Phosphate fluxes.** Phosphate uptake and efflux assays in HEK293T cells were carried out as previously described¹⁰. Amount of phosphate uptake was calculated from the concentration of cold phosphate in the medium multiplied by the ratio of cellular [³³P]phosphate to total [³³P]phosphate supplemented within a period of 30 minutes. Percentage of phosphate efflux was calculated as the ratio of released [³³P]phosphate to total cellular [³³P]phosphate. Phosphate uptake and efflux assays in peripheral blood mononuclear cells (PBMC) were performed as previously described⁶,¹¹.

**Immunoblotting.** Whole cell extracts were separated by 12% SDS-PAGE under reducing conditions and transferred to PVDF membranes. Proteins were detected with antibodies against HA (3F10, Roche; 1:5000) or β-actin (A5441, Sigma Aldrich; 1:5000). Proteins of interest were detected with horseradish peroxidase (HRP)-conjugated anti-mouse or anti-rat antibodies (Southern Biotech; 1:5000). Visualization was performed with Pierce ECL western blotting substrate (Thermo Scientific).

**Flow cytometry.** Cell surface expression of phosphate transporters was monitored on HEK293T cells with soluble ligands derived from the receptor-binding domain (RBD) of different retroviral Envs. RBD derived from mouse X-MLV, koala endogenous retrovirus, and mouse amphotropic-MLV Envs, were used to detect XPR1/SLC53A1, PiT1/SLC20A1 and PiT2/SLC20A2, respectively, and binding assays performed as previously described¹⁰,³⁵,³⁶. Control expression of the GLUT1 glucose transporter was monitored with the GLUT1-specific H2RBD ligand, derived from the human T-cell leukemia virus type 2 Env³⁰,³³. Briefly, 5 × 10⁵ cells were incubated in PBA (PBS complemented with 2% FBS) containing the adequate RBD, for 30 min at 37 °C under agitation, followed by two washes with PBA and incubation with an Alexa Fluor 488-conjugated anti-mouse IgG1 antibody (Life Technologies; 1:5000) for 20 min at 4 °C. Cells were promptly analyzed on NovoCyte flow cytometer (Becton Dickinson), and data were analyzed with the FlowJo package. All RBD ligands were produced as previously described¹⁰,³⁵, or obtained from METAFORA-biosystems.

**Virus production.** EGFP virus vectors were produced in 2 × 10⁶ HEK293T cells in 10 cm dishes, co-transfected using calcium phosphate with a combination of 10 µg of the MLV-based LNCG retroviral vector carrying the EGFP reporter gene³⁷, with 5 µg of the MLV Gag/Pol expression vector (pC57GPBEB)³⁸, and 5 µg of either the vesicular stomatitis virus (VSV) G protein, or the X-MLV Env expression vectors. Virion-containing media were harvested 2 days later, filtered through 0.45 µm (pore size) filter and stored at -80 °C before use. LNCX virus vectors carrying WT and mutated XPR1/SLC53A1 genes were produced in the same conditions, except that the LNCG retroviral vector was replaced by the various XPR1/SLC53A1 LNCX vectors.
**Selection of cells stably producing XPR1 variants, and retrovirus infection.** CHO cells stably expressing XPR1 constructs were generated by transducing CHO cells with the pLNC(XPR1), pLNC(XPR1R459C), pLNC(XPR1N619D) or pLNC(XPR1I629S) vectors, or the empty pLNCX vector, and selected the next day with medium containing 1.5 mg per ml of G418 (active fraction). G418-resistant clones were pooled after 2 weeks of selection before further experiments.

CHO cells that stably expressed WT, p.R459C, p.N619D or p.I629S XPR1/SLC53A1 were plated in 12-well plates with 2 × 10⁴ cells per well. Infections were performed the following day with serial dilutions of replication-defective LNCG retroviral vector pseudotyped with either X-MLV Env or VSV-G glycoproteins. Percentages of EGFP-positive cells were analyzed by flow cytometry.

**Statistical analysis.** The Student's t test was applied using GraphPad Prism 6 software; P values were as follows: *P* ≤ 0.05; **P** ≤ 0.01; ***P** ≤ 0.001; ****P** ≤ 0.0001.

**References**

1. Nicolas, G. et al. Brain calcification process and phenotypes according to age and sex: Lessons from SLC20A2, PDGFB, and PDGFRB mutation carriers. *Am J Med Genet B Neuropsychiatr Genet* **168**, 586–594, https://doi.org/10.1002/ajmg.b.32336 (2015).

2. Ramos, E. M., Oliveira, J., Sobrido, M. J. & Coppola, G. In *GeneReviews(R)* (eds M. P. Adam *et al*. ) (2004).

3. Wang, C. et al. Mutations in SLC20A2 link familial idiopathic basal ganglia calcification with phosphate homeostasis. *Nat Genet* **44**, 254–256, https://doi.org/10.1038/ng.1077 (2012).

4. Nicolas, G. et al. Mutation of the PDGFRB gene as a cause of idiopathic basal ganglia calcification. *Neurology* **80**, 181–187, https://doi.org/10.1212/WNL.0b013e31827ccf34 (2013).

5. Keller, A. et al. Mutations in the gene encoding PDGF-B cause brain calcifications in humans and mice. *Nat Genet* **45**, 1077–1082, https://doi.org/10.1038/ng.2723 (2013).

6. Legati, A. et al. Mutations in XPR1 cause primary familial brain calcification associated with altered phosphate export. *Nat Genet* **47**, 579–581, https://doi.org/10.1038/ng.3289 (2015).

7. Yao, X. P. et al. Biallelic Mutations in MYORG Cause Autosomal Recessive Primary Familial Brain Calcification. *Neuron* **98**, 1116–1123 e1115, https://doi.org/10.1016/j.neuron.2018.05.037 (2018).

8. Bottger, P. et al. Characterization of transport mechanisms and determinants critical for Na+-dependent Pi symport of the PiT family paralogs human PiT1 and PiT2. *Am J Physiol Cell Physiol* **291**, C1377–1387, https://doi.org/10.1152/ajprenal.zh2-4590-corr.2006 (2006).

9. Kavanaugh, M. P. et al. Cell-surface receptors for gibson ape leukemia virus and amphotropic murine retrovirus are inducible sodium-dependent phosphate symporters. *Proc Natl Acad Sci USA* **91**, 7071–7075 (1994).

10. Giovannini, D., Touhami, J., Charnet, P., Sitbon, M. & Battini, J. L. Inorganic phosphate export by the retrovirus receptor XPR1 in metazoans. *Cell Rep* **3**, 1866–1873, https://doi.org/10.1016/j.celrep.2013.05.035 (2013).

11. Anheim, M. et al. XPR1 mutations are a rare cause of primary familial brain calcification. *J Neurol* **263**, 1559–1564, https://doi.org/10.1007/s00415-016-8166-4 (2016).

12. Wild, R. et al. Control of eukaryotic phosphate homeostasis by inositol polyphosphate sensor domains. *Science* **352**, 986–990, https://doi.org/10.1126/science.aad9858 (2016).

13. Ramos, E. M. et al. Primary brain calcification: an international study reporting novel variants and associated phenotypes. *Eur J Hum Genet*, https://doi.org/10.1038/s41431-018-0185-4 (2018).

14. Schwarz, J. M., Cooper, D. N., Schuelke, M. & Seelow, D. MutationTaster2: mutation prediction for the deep-sequencing age. *Nat Methods* **11**, 361–362, https://doi.org/10.1038/nmeth.2890 (2014).

15. Adzhubei, I. A. et al. A method and server for predicting damaging missense mutations. *Nat Methods* **7**, 248–249, https://doi.org/10.1038/nmeth0410-248 (2010).

16. Ng, P. C. & Henikoff, S. Predicting deleterious amino acid substitutions. *Genome Res* **11**, 863–874, https://doi.org/10.1101/gr.176601 (2001).

17. Lek, M. et al. Analysis of protein-coding genetic variation in 60,706 humans. *Nature* **536**, 285–291, https://doi.org/10.1038/nature19057 (2016).

18. Tailor, C. S., Nouri, A., Zhao, Y., Takeuchi, Y. & Kabat, D. A sodium-dependent neutral-amino-acid transporter mediates infections of feline and baboon endogenous retroviruses and simian type D retroviruses. *Journal of virology* **73**, 4470–4474 (1999).

19. Battini, J. L., Rasko, J. E. & Miller, A. D. A human cell-surface receptor for xenotropic and polytropic murine leukemia viruses: possible role in G protein-coupled signal transduction. *Proc Natl Acad Sci USA* **96**, 1385–1390 (1999).

20. Yang, Y. L. et al. Receptors for polytropic and xenotropic mouse leukaemia viruses encoded by a single gene at Rmc1. *Nat Genet* **21**, 216–219, https://doi.org/10.1038/6005 (1999).

21. Ferreira, J. B. et al. First report of a de novo mutation at SLC20A2 in a patient with brain calcification. *J Mol Neurosci* **54**, 748–751, https://doi.org/10.1007/s12031-014-0357-9 (2014).

22. Nicolas, G. et al. A de novo nonsense PDGFB mutation causing idiopathic basal ganglia calcification with laryngeal dystonia. *Eur J Hum Genet* **22**, 1236–1238, https://doi.org/10.1038/ejhg.2014.9 (2014).

23. Ansermet, C. et al. Renal Fanconi Syndrome and Hypophosphatemic Rickets in the Absence of Xenotropic and Polytropic Retroviral Receptor in the Nephron. *J Am Soc Nephrol* **28**, 1073–1078, https://doi.org/10.1681/ASN.2016070726 (2017).

24. Jensen, N., Autzen, J. K. & Pedersen, L. Slc20a2 is critical for maintaining a physiologic inorganic phosphate level in cerebrospinal fluid. *Neurogenetics* **17**, 125–130, https://doi.org/10.1007/s10048-015-0469-6 (2016).

25. Jensen, N. et al. Loss of function of Slc20a2 associated with familial idiopathic Basal Ganglia calcification in humans causes brain calcifications in mice. *J Mol Neurosci* **51**, 994–999, https://doi.org/10.1007/s12031-013-0085-6 (2013).

26. Wallingford, M. C. et al. SLC20A2 Deficiency in Mice Leads to Elevated Phosphate Levels in Cerebrospinal Fluid and Glymphatic Pathway-Associated Arteriolar Calcification, and Recapitulates Human Idiopathic Basal Ganglia Calcification. *Brain Pathol* **27**, 64–76, https://doi.org/10.1111/bpa.12362 (2017).

27. McGuffin, L. J., Bryson, K. & Jones, D. T. The PSIPRED protein structure prediction server. *Bioinformatics* **16**, 404–405 (2000).

28. Ward, J. J., Sodhi, J. S., McGuffin, L. J., Buxton, B. F. & Jones, D. T. Prediction and functional analysis of native disorder in proteins from the three kingdoms of life. *J Mol Biol* **337**, 635–645, https://doi.org/10.1016/j.jmb.2004.02.002 (2004).

29. Jones, D. T. Improving the accuracy of transmembrane protein topology prediction using evolutionary information. *Bioinformatics* **23**, 538–544, https://doi.org/10.1093/bioinformatics/btl677 (2007).

30. Manel, N., Battini, J. L. & Sitbon, M. Human T cell leukemia virus envelope binding and virus entry are mediated by distinct domains of the glucose transporter GLUT1. *J Biol Chem* **280**, 29025–29029, https://doi.org/10.1074/jbc.M504549200 (2005).

31. Nicolas, G. et al. Phenotypic spectrum of probable and genetically-confirmed idiopathic basal ganglia calcification. *Brain* **136**, 3395–3407, https://doi.org/10.1093/brain/awt255 (2013).

32. Richards, S. et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. *Genet Med* **17**, 405–424, https://doi.org/10.1038/gim.2015.30 (2015).
33. Manel, N. et al. The ubiquitous glucose transporter GLUT-1 is a receptor for HTLV. *Cell* **115**, 449–459 (2003).

34. Miller, A. D. & Rosman, G. J. Improved retroviral vectors for gene transfer and expression. *Biotechniques* **7**, 980–982, 984–986, 989–990 (1989).

35. Petit, V. et al. Optimization of tumor xenograft dissociation for the profiling of cell surface markers and nutrient transporters. *Lab Invest* **93**, 611–621, https://doi.org/10.1038/labinvest.2013.44 (2013).

36. Laval, J. et al. Metabolic adaptation of neutrophils in cystic fibrosis airways involves distinct shifts in nutrient transporter expression. *J Immunol* **190**, 6043–6050, https://doi.org/10.4049/jimmunol.1201755 (2013).

37. Rasko, J. E., Battini, J. L., Gottschalk, R. J., Mazo, I. & Miller, A. D. The RD114/simian type D retrovirus receptor is a neutral amino acid transporter. *Proc Natl Acad Sci USA* **96**, 2129–2134 (1999).

38. Lassaux, A., Sitbon, M. & Battini, J. L. Residues in the murine leukemia virus capsid that differentially govern resistance to mouse Fv1 and human Ref1 restrictions. *Journal of virology* **79**, 6560–6564, https://doi.org/10.1128/JVI.79.10.6560-6564.2005 (2005).

39. Omasits, U., Ahrens, C. H., Muller, S. & Wollscheid, B. Protter: interactive protein feature visualization and integration with experimental proteomic data. *Bioinformatics* **30**, 884–886, https://doi.org/10.1093/bioinformatics/btt607 (2014).

## Acknowledgements

We are grateful to the patients and their families for their contribution to this study, to all collaborators who sent blood samples and medical charts, and to members of our labs for data and reagents sharing, helpful discussions and constant support. We acknowledge the imaging facility MRI, member of the national infrastructure France-BioImaging supported by the French National Research Agency (ANR-10-INBS-04, «Investments for the future»). ULS was supported by successive grants from the Labex EpiGenMed (ANR-10-LABX-12-01), Consejo Nacional de Ciencia y Tecnología (CONACYT - 243449/440597), and the Fondation pour la Recherche Médicale (FRM: FDT20170436884); this work was supported by grants from the Fondation pour la Recherche Médicale (DBI201312285579 to MS); from the French National Research Agency (CALIPHOS, ANR-17-CE14-0008-03, to JLB, MS, and GN); and Labex GR-Ex (ANR-11-LABX-0051); Labex is funded by the program ‘Investissements d'Avenir' of the French National Research Agency. This study was co-supported by European Union and Région Normandie. Europe gets involved in Normandie with the European Regional Development Fund (ERDF). JLB and MS are supported by INSERM.

## Author Contributions

Conceived and designed the experiments: U.L.S., J.L.B., M.S. Patient selection, follow-up and clinical sample collection: A.C.R., D.M., M.C., X.A., C.G., G.N. Performed and interpreted the molecular genetics experiments: A.C.R., G.N. Clinical coordination: G.N. Performed most *ex vivo* and *in vitro* experiments U.L.S.; RBD and transporter detection assays: U.L.S., J.T. Protein modeling: G.L. Analyzed the data, drafted and wrote the manuscript: U.L.S., G.N., J.L.B., M.S.

## Additional Information

**Supplementary information** accompanies this paper at https://doi.org/10.1038/s41598-019-43255-x.

**Competing Interests:** J.L.B. and M.S. are inventors on patents describing the use of RBD ligands; M.S. is the co-founder of METAFORA-biosystems, a start-up company that focuses on metabolite transporters under physiological and pathological conditions.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_7.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019