Proteomics analysis of a human brain sample from a mucolipidosis type IV patient reveals pathophysiological pathways
Orphanet Journal of Rare Diseases volume 16, Article number: 39 (2021)
Mucolipidosis type IV (MLIV), an ultra-rare neurodevelopmental and neurodegenerative disorder, is caused by mutations in the MCOLN1 gene, which encodes the late endosomal/lysosomal transient receptor potential channel TRPML1 (mucolipin 1). The precise pathophysiogical pathways that cause neurological disease in MLIV are poorly understood. Recently, the first post-mortem brain sample became available from a single MLIV patient, and in the current study we performed mass spectrometry (MS)-based proteomics on this tissue with a view to delineating pathological pathways, and to compare with previously-published data on MLIV, including studies using the Mcoln1−/− mouse.
A number of pathways were altered in two brain regions from the MLIV patient, including those related to the lysosome, lipid metabolism, myelination, cellular trafficking and autophagy, mTOR and calmodulin, the complement system and interferon signaling. Of these, levels of some proteins not known previously to be associated with MLIV were altered, including APOD, PLIN4, ATG and proteins related to interferon signaling. Moreover, when proteins detected by proteomics in the human brain were compared with their orthologs detected in the Mcoln1−/− mouse by RNAseq, the results were remarkably similar. Finally, analysis of proteins in human and mouse CSF suggest that calbindin 1 and calbindin 2 might be useful as biomarkers to help chart the course of disease development.
Despite the sample size limitations, our findings are consistent with the relatively general changes in lysosomal function previously reported in MLIV, and shed light on new pathways of disease pathophysiology, which is required in order to understand the course of disease development and to determine the efficacy of therapies when they become available for this devastating disease.
Mucolipidosis type IV (MLIV), an ultra-rare lysosomal storage disease (LSD) primarily found in the Ashkenazi Jewish population , is caused by mutations in the MCOLN1 gene, which encodes the late endosomal/lysosomal transient receptor potential channel TRPML1 (mucolipin 1) . Over two dozen mutations in MCOLN1 have been identified which lead to MLIV . TRPML1 is a non-selective cation channel that regulates lysosomal ion content [4, 5]. TRPML1 also plays a role in autophagy regulation, intracellular trafficking and in the mechanistic target of rapamycin (mTOR)/transcription factor EB (TFEB) signaling axis [6, 7]. MLIV is a neurodevelopmental and neurodegenerative disorder with severe psychomotor developmental delay, impaired gastric function, mental retardation, motor deficits and progressive visual impairment [1, 8]. At present there is no therapy for MLIV and our understanding of brain pathophysiology is limited due to the lack of human tissue samples available for research.
We recently obtained brain tissue and cerebrospinal fluid (CSF) from a 36-year-old male who died from MLIV-related complications. This is the first time that post-mortem brain tissue has been made available from a MLIV patient. We performed mass spectrometry (MS)-based proteomics on this tissue with a view to delineating pathological pathways. Moreover, we compared our data with that obtained using a mouse model of MLIV, namely the Mcoln1−/− mouse [9, 10], which further validated that the Mcoln1−/− mouse is a genuine animal model of MLIV, and of more importance, allowed us to identify some novel pathological components and pathways in the human MLIV brain tissue, such as activation of complement and interferon (IFN) pathways and lipid droplet and lipoprotein metabolism. Along with the identification of a potential biomarker (calbindin 2) for early disease development in the CSF, our study paves the way for further understanding of the molecular pathways that lead to disease pathophysiology in MLIV, a pre-requisite for development of putative therapies.
Human brain and CSF samples
Human brain and CSF samples were obtained from the NIH NeuroBioBank Brain and Tissue Repository at the University of Maryland, Baltimore. The post mortem MLIV patient was a 36 year-old male and the control was a 34-year-old healthy male who died from unrelated causes.
Brain tissue MS-based proteomics
Brain tissue was lysed using a GentleMACS dissociator (Miltenyi Biotec, Bergisch Gladbach, Germany) in 50 mM Tris HCl containing 5% sodium dodecyl sulfate (SDS) and supplemented with a protease inhibitor cocktail (1:200, Sigma-Aldrich, Missouri, USA). Homogenates were incubated with 5 mM dithiothreitol for 1 h at 56 °C followed by 10 mM iodoacetamide in the dark at room temperature for 45 min. Samples were loaded onto S-Trap microcolumns (Protifi, New-York, USA) and subjected to tryptic digestion [11, 12]. Samples were vacuum dried and stored at -80 °C. The resulting peptides were analyzed by nanoflow liquid chromatography (NanoAcquity, Waters, Milford, MA, USA) coupled to high resolution, high mass accuracy mass spectrometry (Fusion Lumos, Thermo Scientific, Waltham, MA, USA) and sequentially analyzed in discovery mode. Raw data was processed with MaxQuant v18.104.22.168 . Data was analyzed with the Andromeda search engine against the human proteome database (www.uniprot.com) and appended with common laboratory protein contaminants. Quantitative comparisons were calculated using Perseus v22.214.171.124. Pathway analysis was performed using Gene Analytics . Proteins were categorized based on the Human Protein Atlas  (http://www.proteinatlas.org).
CSF MS-based proteomics
Samples were lysed with 8 M urea and incubated at room temperature for 30 min prior to incubation with 5 mM dithiothreitol (Sigma-Aldrich) at room temperature for 1 h, followed by alkylation with 10 mM iodoacetamide (Sigma-Aldrich) in the dark at room temperature for 45 min. Samples were diluted in 2 M urea with 50 mM ammonium bicarbonate. Proteins were digested with trypsin (50:1; w/vol) (Promega, Madison, WI, USA) overnight at 37 °C followed by a second trypsin digestion for 4 h. Digestion was terminated using trifluoracetic acid (1% final concentration). Peptides were desalted using an Oasis hydrophilic lipophilic balance column (Waters) and then vacuum dried and stored at -80 °C. The resulting peptides were analyzed using nanoflow liquid chromatography (nanoACQUITY, Waters) coupled to high resolution, high mass accuracy mass spectrometry (Q Exactive HFX for mouse CSF or Fusion Lumos for human CSF, Thermo Scientific). Samples were analyzed in random order in discovery mode. The raw data was processed with MaxQuant v126.96.36.199  and analyzed with the Andromeda search engine against the mouse or human proteome database (www.uniprot.com), and appended with common laboratory protein contaminants. Quantitative comparisons were calculated using Perseus v188.8.131.52. Pathway analysis was done using Gene Analytics .
Mcoln1−/− mice on a C57Bl/6 J background  were used. Littermates (Mcoln1+/+ or Mcoln1±) were used as controls and both females and males were used in all studies. Genotyping was performed by PCR using genomic DNA extracted from mouse tails. Mice were maintained in the Experimental Animal Center of the Weizmann Institute of Science. Animal experiments were approved by the Weizmann Institute Institutional Animal Care and Use Committee.
1-, 2-, 3- and 7-month-old Mcoln± or Mcoln1−/− mice were anesthetized using 200 mg/ml sodium pentobarbital and CSF aspirated from the cisterna magna using a glass needle. CSF samples were kept at -80 °C.
Brain tissue was homogenized using a GentleMACS dissociator and mRNA was isolated using the RNeasy mini kit (Qiagen GmbH, Hilden, Germany). RNA concentration (ratio 260/230 and 260/280 nm) was measured using a NanoDrop ND-1000 (Thermo Scientific, Waltham, MA, USA) and RNA integrity evaluated using an RNA screen tape on a Tapestation 2200 (Agilent, Santa Clara, California, USA). A bulk variation of MARSseq  was used to construct RNAseq libraries. Sequencing was performed using an Illumina Nextseq-500 75 cycle high output kit (Illumina, San Diego, California, USA; paired end sequencing). Raw reads were mapped to the Mus musculus genome (mm10) using STAR. Only reads which mapped uniquely to genes were considered for further analysis. Differentially-expressed genes (DEGs) were selected using a 1.5-fold change cutoff between two populations and adjusted p value for multiple gene testing of < 0.05 . Pathway analysis was performed using Gene Analytics . Ortholog conversion was done using HUGO Gene Nomenclature Committee.
Enzyme-linked immunosorbent assay
Calbindin 2 was quantified using a Calertinin ELISA kit (Abcam, Cambridge, MA) according to manufacturer’s instructions. Samples were analyzed in triplicate.
MS-based proteomics identified 3968 proteins in the cerebellum (CB) (with more than one unique peptide) and 3936 in the cerebral cortex (CRB). Of these, 959 were increased and 1512 decreased in MLIV compared to the control in the CB, and 824 proteins increased and 1168 decreased in the CRB (ratio > 1.5). 2113 proteins were detected in the CSF, with 1018 elevated (ratio > 1.5) and 526 proteins decreased in MLIV compared to the control. Due to the minimal sample size, label-free, relative quantification was calculated by dividing protein intensity in the MLIV sample by protein intensity in the control [18, 19]. When a ratio could not be calculated, i.e., protein levels were below the level of detection in one or other of the samples, changes are indicated by ‘ + ’ for a protein that was detected only in the MLIV sample, and by ‘-’ when the protein was detected only in the control sample (see Tables 1, 2, 3, 4, 5, 6 and 7 and 9).
To obtain an overview of changes in the MLIV brain, we first interrogated the proteomics data by pathway analysis. Identical pathways (i.e. sphingolipid metabolism, interferon signaling, the lysosome and others, see Fig. 1a, b) were upþregulated in both the CB and the CRB and likewise similar pathways were down-regulated in the CB and CRB (Fig. 1c, d). The latter includes the GABAergic pathway (including GABRA1, GABRA3, GABRB1, GABRB2, GABBR2) and other neuronal genes (GRIN2B, HPCA, NEGR1, GPM6A, and NPTX1; not shown), along with proteins associated with L1 cell adhesion (L1CAM, Fig. 1.c, d), which plays a role in nervous system development . The identification of similar pathways in two distinct brain regions adds validity to the accuracy of the data, which is taken from one MLIV brain and one control brain, but analyzed experimentally in two different brain areas from each sample. Moreover, the up- and down-regulated pathways are largely consistent with what is known to date on MLIV pathophysiology in Mcoln1−/− mice [3, 7].
We next examined individual proteins associated with some of these pathways. Levels of most of the lysosomal proteins that were detected by MS increased in MLIV brain including lysosomal hydrolases, cathepsins and a number of other lysosomal proteins (Table 1). Less lysosomal proteins were detected in the CSF than in the CB and CRB, but levels of all detected proteins were similarly increased. The up-regulation of lysosomal proteins is consistent with the broad lysosomal defects resulting from the loss of TRPML1 function that underlies MLIV pathology .
Levels of a number of proteins in the sphingolipid metabolic pathway increased in MLIV CB and CRB (although most were not detected in the CSF), as were some proteins associated with cholesterol metabolism (MVD, MVK, MSMO1, SEC14L2), which were elevated in the CB but not CRB (Table 2), also consistent with some of the changes in lipid metabolism previously reported in MLIV [22, 23]. Levels of PLIN4, a protein associated with lipoprotein metabolism and lipid droplets increased, as did APOD (Table 2).
Changes in lipid metabolism are often related to disrupted myelination  which is also detected in MLIV [25,26,27]. Interestingly, levels of proteins associated with myelination were reduced in the CB, but not in the CRB in the MLIV brain, where their number was increased (Table 3). The reason for the latter is not entirely clear. In the Mcoln1−/− mouse, myelination defects are found primarily in the CRB [9, 28]; however, our data suggest a regional variability in hypomyelination in different human brain regions, and reduced levels of myelin proteins may reflect degeneration in the cerebellum in this patient as has been reported for other older MLIV patients .
TRPML1 has been suggested to indirectly regulate membrane traffic [23, 29]. A relatively large number of proteins associated with trafficking between the ER, Golgi and the lysosome decreased, mainly in the CB, with many of the proteins detected only in the control, but not in the MLIV sample (Table 4). Likewise, changes were detected in a number of autophagy-related proteins (Table 4), consistent with the constitutive block of autophagy reported in MLIV . Several proteins related to autophagy were elevated in MLIV including the autophagosome markers MAP1LC3A and MAP1LC3B , indicative of the inhibited autophagosome degradation of lysosomes, whereas a number of ATG proteins, which regulate autophagosome formation  were down-regulated (Table 4), supporting the notion that lysosomes are unable to fuse with autophagosomes in MLIV . Among the key regulators of autophagy is TFEB . The activity and localization of TFEB is regulated by mechanistic target of rapamycin (mTOR) phosphorylation. TRPML1 has been suggested to be part of a feedback loop between mTOR and TFEB [34, 35] which could explain the elevation of LAMTOR proteins in MLIV CB and in MLIV CRB (Table 5). Moreover, lysosomal calcium release through TRPML1 is required for mTORC1 activation  and TRPML1 is essential for MTORC1 recruitment and activation in a Ca2 + -calmodulin-dependent manner , consistent with the decrease in levels of several calmodulin-associated proteins in MLIV (Table 5).
Neuroinflammation in MLIV has been reported, with activation of microglia and astrocytes in Mcoln1−/− mouse brain [25, 38], and we observed elevation in a number of markers such as MAC2, LGALS3BP, GFAP and TGFBI. Unexpectedly, levels of a number of proteins associated with the complement system were also elevated in the MLIV CB and CRB, and also in the CSF (Table 6), suggesting involvement of the classical complement pathway . Levels of proteins related to the type 1 interferon (IFN) response (Fig. 1a, b) were also increased, with 17 of the 19 detected IFN-related proteins increased in MLIV in the CB and 14 elevated in the CRB (Table 7). These proteins include IFN-stimulated proteins such as ISG15, GBP1, IFIT1 and IFIT3. Only 5 of the IFN-related proteins were detected in the CSF, and all were similarly increased, indicating activation of the IFN response in MLIV patients.
In addition to proteomics analysis of human brain, we also performed RNAseq on the brain of Mcoln1−/− mice  at the end-stage of disease, i.e. at ~ 7 months of age (Fig. 1e). 184 genes were differentially expressed (DE) (fold-change ≥ 1.5, p adjusted ≤ 0.05) with 174 up-regulated and 10 down-regulated. When genes detected in the mouse by RNAseq overlapped with orthologous proteins detected by proteomics in the human brain, the results were remarkably similar. Thus, levels of RNA encoding complement genes were elevated in 7-month old Mcoln1−/− mouse (Table 6), as were levels of IFN-related genes (Table 7), confirming the validity of the human proteomics data. The overlap between the patient and mouse datasets also highlights the strength of the Mcoln1−/− mouse as an accurate pre-clinical model of MLIV. A number of additional pathways were up-regulated in the mouse RNAseq data, which were not detected in the proteomics analysis of human brain, including inflammatory genes such as cytokines and chemokines. Moreover, many of the genes and pathways that were detected in 7-month-old Mcoln1−/− mice by RNAseq (or detected by proteomics analysis of human brain) were also detected at earlier pre-symptomatic stages, i.e. 1 month and early symptomatic (2/3-month old mice). Thus, even at early pre-symptomatic stages, levels of genes associated with the complement and IFN pathways were up-regulated. Likewise, microglia activation also begins early in disease development  as RNAseq demonstrated that Mpeg1, a macrophage marker expressed also in microglia, was DE in 1-month Mcoln1−/− mice, as was Cxcl10, whose levels are also altered in models of experimental autoimmune encephalomyelitis .
Levels of a number of other genes, including cathepsins and genes related to sphingolipid metabolism, were also up-regulated at early stages of disease in Mcoln1−/− mice, as were genes related to microglia activation  which included Trem2 and Tyrobp, key disease-associated microglia (DAM) genes . Genes associated with defective myelination were down-regulated  (Table 8). Defective myelination is a hallmark of MLIV pathology [25, 27] and this data indicates that changes in the expression of myelination genes, or in levels of proteins associated with myelination, are most significant at the active stage of brain myelination in early development.
Finally, we performed MS-based proteomics on CSF from the human brain patient and from CSF obtained from the Mcoln1−/− mice. In the latter, 2092 proteins were identified and quantified. 67 proteins were DE (Mcoln1−/− versus Mcoln1±) at 1 month, 99 were DE at 2 months, 49 at 3 months and 162 at 7-months of age (fold-change ≥ 2, p < 0.05, number of peptides ≥ 2; not shown). While the same proteins were not detected at each of these ages, the deregulated proteins belonged to similar pathways. Thus, proteins related to a number of pathways increased, including lysosomal proteins (i.e. LAMP1, GNS, PPT1, TPP1, CTSD, CTSB, CTSL, CTSZ, PRCP), Ca+-binding proteins and proteins that regulate Ca+ levels (i.e. RCN1, SLC3A2, PPP3R1), proteins related to the immune system (i.e. LAG3, MAC2, LGALSL, LYZ2, RELT, TREM2 and complement C4B, C1QB, CFP), and several proteins associated with cellular trafficking (LMAN1, COPE, TOM1, VPS26A, STX7). Levels of proteins related to apolipoprotein metabolism were decreased (APOM, APOBEC2, APOA4). As no common proteins were detected, we then compared all Mcoln1−/−mice to all Mcoln1± mice, independent of their age. 63 proteins were DE (48 up-regulated and 15 down-regulated). Levels of the Purkinje cell protein PCP4 were increased, as were levels of most lysosomal proteins (Table 9). Levels of A1BG (alpha-1B-glycoprotein, Table 9) decreased by the largest amount in mouse CSF. A1BG levels were also decreased in the MLIV patient CSF (Table 9). A1BG has been suggested as a biomarker [42,43,44].
Three of the top 5 up-regulated proteins in the CSF from Mcoln1−/− mice are related, namely calbindin 1 (CALB1), calbindin 2 (CALB2) and PCP4 (Table 9). Strikingly, levels of CALB2 and PCP4 were elevated in the CSF from the human MLIV patient; CALB2 and PCP4 were not detected in the CSF of the control human brain. Calbindins have been suggested to act as biomarkers in diseases with cerebellar involvement  including LSDs [46, 47]. To determine whether either of these proteins might be a suitable biomarker in human MLIV, we analyzed levels of CALB2 in CSF; note that due to the limited amount of CSF available to us, we were only able to evaluate levels of one of these two proteins. By ELISA, CALB2 levels were 2247.5 pg/ml CSF in the control patient and 2838.5 pg/ml in the CSF from the MLIV patient. Likewise, levels of CALB2 were elevated in mouse CSF as early as 1 month of age, and detected in 69% of Mcoln1−/− mice but in only 12% of Mcoln1± mice (Table 10), suggesting that CALB2 might indeed act as a biomarker for MLIV which could also be used to test the efficacy of therapies as and when they become available.
Until the current study, the vast majority of information on MLIV neuropathology was obtained from clinical observation , studies in cell culture  or from study of the Mcoln1−/− mouse . The current study extends this data by analyzing the first available human brain sample from an MLIV patient. Having said that, our study is somewhat unconventional inasmuch as we make conclusions about pathophysiological pathways in MLIV based on a proteomics study comparing only one human brain sample with one control, although proteomics was performed on two different brain areas, namely the cortex and cerebellum, and on the CSF. Since this is the first time that such tissue has been made available for research purposes, and despite the obvious limitations of this approach, we are of the opinion that such a study is warranted and of value for both the research, clinical and patient communities associated with MLIV. Nevertheless, conclusions must be drawn with suitable caution taking into account the minimal sample size.
A broad overview of the pathophysiology in MLIV suggests that impairment of lysosome-related cellular pathways leads to brain pathology and changes in neurological function (Fig. 2). Since TRMPL1 encodes for a late endosomal/lysosomal channel involved in regulation of a broad array of lysosomal processes, the effect of reducing its activity might be more general than that seen in LSDs which are caused, for instance, by the loss of a specific lysosomal hydrolase. Thus, loss of TRPML1 leads to more wide-ranging changes in lysosomal function (supported by changes in levels of a large number of lysosomal proteins in the current study, along with defective autophagy (Fig. 2)). Since the brain tissue was taken post-mortem, it is obviously not possible to define which pathways are primary or secondary, but since TRMPL1 directly impinges upon the lysosome, it is reasonable to assume that the lysosome is first affected. Inhibition of autophagy in MLIV is likely a result of inhibited lysosomal function, and in turn may lead to dysregulation of upstream mTOR and TFEB signaling (Fig. 2), although recent studies  on the role of TRPML1 in regulation of TFEB suggest that changes in TFEB signaling could also be a direct result of TRPML1 loss. The temporal relationship between these changes and mitochondrial abnormalities, along with changes in lipoprotein metabolism, cannot be directly established in the current study, but the relationships between autophagy and mitochondrial defects are well documented [50, 51]. These altered cellular pathways lead to a complex course of neuropathophysiological changes, including aberrant neuronal function and neuroinflammation, with the latter including the interferon response, complement activation, microgliosis and astrocytosis (Fig. 2). The time-course of these events cannot be determined in autopsy MLIV samples, although a number of these pathways are altered at the early stages in disease development, at least in the mouse model, as ascertained by the RNAseq data, and are therefore likely to be the direct consequence of the loss of TRPML1.
While a number of these pathways (Fig. 2) were previously implicated in MLIV pathophysiology, we now provide a more comprehensive picture of some of them and identify a number of new proteins that were not previously known to be identified with MLIV. For instance, levels of a number of lysosomal hydrolases were increased. Earlier studies had shown similar findings, including studies in fibroblasts derived from MLIV patients [23, 52]. The elevation we report in cathepsin levels is similar to that observed in Mcoln1−/− mice microglia , in TRPML1-deficient cells  and in astrocytic cultures , as is the elevation in LAMP1 [9, 23, 38]. Thus, the increase in levels of lysosomal proteins is consistent with the enlarged lysosomes reported in MLIV retina .
The changes in apolipoproteins are much more extensive than previously reported, since only changes in ApoE were reported in MLIV microglia  and in astrocytes from Mcoln1−/− mice . We now report a significant increase in levels of APOD in the CB and CRB. APOD is mainly expressed in glia and neurons  and is elevated in some other neurological disorders such as Alzheimer’s disease, specifically in the CSF , in multiple sclerosis , in the aging brain  and in neuronal tissues from Niemann Pick C disease . APOD and APOE have also been implicated in Parkinson’s disease . APOD and APOE, are involved in the removal of lipids during nerve cell degeneration , which may be of importance in demyelination. Moreover, levels of two PLIN proteins (PLIN 3 and 4) were elevated in the MLIV brain and in neuronal forms of Gaucher disease [63, 64]. Lipid droplets have been implicated in neurodegeneration [65, 66] including in Parkinson’s disease, which suggest that they play a broad role in neurodegenerative diseases.
Likewise, we observed down-regulation of proteins associated with trafficking between the ER, Golgi and the lysosome, in line with the defects in retrograde trafficking of late endosome and lysosome to Golgi retrograde trafficking that were reported in TRPML1-deficient cells [67, 68].
Defective intracellular trafficking can be reflected in changes in autophagic pathways. Previous studies in neuronal cultures from Mcoln1−/− mice demonstrated increased levels of the classic autophagic marker, LC3-II , but our current data allows us to further dissect autophagy pathway inasmuch as a number of the autophagosome markers, i.e. ATG proteins, which control autophagosome formation , are decreased. While interpretation of the LC3 readout in the post-mortem tissue is limited due to inability to track the kinetics of the autophagy process, increased levels of LC3 in the context of inhibited lysosomal function and decreased expression of the proteins involved in the upstream stages of autophagosome formation, clearly indicate block of autophagy . Increased levels of LC3 in this case result from diminished degradation of autophagosomes on lysosomes due to inhibited lysosomal function. TRPML1 is known to induce autophagosome biogenesis through a signaling pathway involving activation of calcium/calmodulin-dependent protein kinase β , consistent with our observation of decreased levels of CAMKK proteins in the MLIV patient brain. Likewise, TRPML1 overexpression leads to mTOR inhibition and activation of its targets, including TFEB and ULK1 ; however, we cannot validate changes in TFEB activation in our study since it was not detected by the proteomics analysis.
Neuroinflammation, activation of the complement system, and defective myelination were also observed in the MLIV brain. For the former, noticeable pathways of neuroinflammation include activation of the interferon pathway, as previously demonstrated in Mcoln1−/− microglia  and astrocytes , and similar to that previously observed in another LSD, Gaucher disease [63, 71]. However, in Gaucher disease, the interferon pathway was categorized as a secondary disease pathway, since eliminating this pathway had no effect on the course of the disease in mice ; whether this is also the case in MLIV is currently unknown.
Since we were also able to obtain small amounts of CSF from the human MLIV patient, we performed a preliminary analysis to determine whether changes in the brain might be reflected in the CSF, and thus indicate that analysis of CSF might be useful to identify biomarkers. The CSF proteomics data was largely consistent with the brain data, but less so with CSF aspirated from Mcoln1−/− mice. However, one specific protein detected in both human and mouse CSF, namely CALB2, was significantly elevated in both the proteomics analysis and was confirmed by ELISA, similar to that observed in late infantile neuronal ceroid lipofuscinosis type 2 (CLN2) and CLN3 patients . Due to the low amount of CSF available to us, we were unable to verify whether CALB1 might also act as a biomarker, but CALB1 has also been suggested as a biomarker in Niemann-Pick type C1 . Thus, while not yet being fully validated, our data suggest that further analysis of CALB1 and CALB2 might be worthwhile in light of their potential as a biomarker for MLIV.
We analyzed the proteome in two brain regions from an MLIV patient and a healthy aged- and gender-matched control, and annotated the observed changes in protein levels with changes previously observed in in vitro and in vivo models of MLIV. Our data confirm the usefulness of the Mcoln1−/− mouse as a genuine model of MLIV, and more importantly, shed light on some new putative players in MLIV pathology, such as apolipoproteins, lipid droplets, the IFN pathway, disease-associated microglia activation and autophagy. Finally, analysis of CSF from both the human patient and from Mcoln1−/− mice suggest that CALB1 and CALB2 might be useful as biomarkers to help chart the course of disease development and the outcomes of pre-clinical and clinical trials.
Availability of data and materials
The RNAseq dataset generated during the current study was deposited in the Gene Expression Omnibus (GEO) database, www.ncbi.nlm.nih.gov/geo (accession no. GSE157825). Mass spectrometry data was deposited to the ProteomeXchange Consortium via the PRIDE  partner repository with the dataset identifiers PXD020490, PXD020491 and PXD020494.
Cyclic adenosine monophosphate
Central nervous system
Lysosomal storage disease
Mechanistic target of rapamycin
Mucolipidosis type IV
Transcription factor EB
Two pore channels
Transient receptor potential channel
Altarescu G, Sun M, Moore DF, et al. The neurogenetics of mucolipidosis type IV. Neurology. 2002;59:306–13.
Sun M, Goldin E, Stahl S, et al. Mucolipidosis type IV is caused by mutations in a gene encoding a novel transient receptor potential channel. Hum Mol Genet. 2000;9:2471–8.
Wakabayashi K, Gustafson AM, Sidransky E, et al. Mucolipidosis type IV: an update. Mol Genet Metab. 2011;104:206–13.
Colletti GA, Kiselyov K. TRPML1. In: Transient Receptor Potential Channels.Vol 704. Advances in Experimental Medicine and Biology. Dordrecht: Springer, Dordrecht; 2011:209–219.
Wang W, Zhang X, Gao Q, et al. TRPML1: an ion channel in the lysosome. Handb Exp Pharmacol. 2014;222:631–45.
Abe K, Puertollano R. Role of TRP channels in the regulation of the endosomal pathway. Physiology. 2011;26:14–22.
Boudewyn LC, Walkley SU. Current concepts in the neuropathogenesis of mucolipidosis type IV. Kielian T, ed. J Neurochem 2019;148:669–689.
Chitayat D, Meunier CM, Hodgkinson KA, et al. Mucolipidosis type IV: Clinical manifestations and natural history. Am J Med Genet. 1991;41:313–8.
Micsenyi MC, Dobrenis K, Stephney G, et al. Neuropathology of the Mcoln1−/− knockout mouse model of mucolipidosis type IV. J Neuropathol Exp Neurol. 2009;68:125–35.
Venugopal B, Browning MF, Curcio Morelli C, et al. Neurologic, gastric, and opthalmologic pathologies in a murine model of mucolipidosis type IV. Am J Hum Genet. 2007;81:1070–83.
Zougman A, Selby PJ, Banks RE. Suspension trapping (STrap) sample preparation method for bottom-up proteomics analysis. Proteomics. 2014;14:1006–1000.
Elinger D, Gabashvili A, Levin Y. Suspension trapping (S-Trap) is compatible with typical protein extraction buffers and detergents for bottom-up proteomics. J Proteome Res. 2019;18:1441–5.
Cox J, Mann M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat Biotechnol 2008;26:1367–1372.
Ben-Ari Fuchs S, Lieder I, Stelzer G, et al. GeneAnalytics: an integrative gene set analysis tool for next generation sequencing. RNAseq Microarray Data OMICS. 2016;20:139–51.
Uhlén M, Fagerberg L, Hallström BM, et al. Tissue-based map of the human proteome. Science. 2015;347:1260419–1260419.
Jaitin DA, Kenigsberg E, Keren-Shaul H, et al. Massively parallel single-cell RNA-seq for marker-free decomposition of tissues into cell types. Science. 2014;343:776–9.
Kohen R, Barlev J, Hornung G, et al. UTAP: user-friendly transcriptome analysis pipeline. BMC Bioinformatics. 2019;20:1–7.
Cox J, Hein MY, Luber CA, et al. Accurate proteome-wide label-free quantification by delayed normalization and maximal peptide ratio extraction. Termed MaxLFQ Mol Cell Proteomics. 2014;13:2513–26.
Tyanova S, Temu T, Sinitcyn P, et al. The Perseus computational platform for comprehensive analysis of (prote)omics data. Nat Methods. 2016;13:731–40.
Linneberg C, Toft CLF, Kjaer-Sorensen K, et al. L1cam-mediated developmental processes of the nervous system are differentially regulated by proteolytic processing. Sci Rep. 2019;9:1–14.
Kogot-Levin A, Zeigler M, Ornoy A, et al. Mucolipidosis Type IV: The effect of increased lysosomal pH on the abnormal lysosomal storage. Pediatr Res. 2009;65:686–90.
Chen C-S, Bach G, Pagano RE. Abnormal transport along the lysosomal pathway in Mucolipidosis, type IV disease. Proc Natl Acad Sci U S A. 1998;95:6373–8.
Miedel MT, Rbaibi Y, Guerriero CJ, et al. Membrane traffic and turnover in TRP-ML1–deficient cells: a revised model for mucolipidosis type IV pathogenesis. J Exp Med. 2008;205:1477–90.
Chrast R, Saher G, Nave KA, et al. Lipid metabolism in myelinating glial cells: lessons from human inherited disorders and mouse models. J Lipid Res. 2011;52:419–34.
Grishchuk Y, Sri S, Rudinskiy N, et al. Behavioral deficits, early gliosis, dysmyelination and synaptic dysfunction in a mouse model of mucolipidosis IV. Acta Neuropathol Commun. 2014;2:133.
Frei KP, Patronas NJ, Crutchfield KE, et al. Mucolipidosis type IV: Characteristic MRI findings. Neurology. 1998;51:565–9.
Grishchuk Y, Peña KA, Coblentz J, et al. Impaired myelination and reduced brain ferric iron in the mouse model of mucolipidosis IV. Dis Model Mech 2015;8:dmm.021154–1601.
Mepyans M, Andrzejczuk L, Sosa J, et al. Early evidence of delayed oligodendrocyte maturation in the mouse model of mucolipidosis type IV. Dis Model Mech 2020;13:dmm044230.
LaPlante JM, Ye CP, Quinn SJ, et al. Functional links between mucolipin-1 and Ca2+-dependent membrane trafficking in mucolipidosis IV. Biochem Biophys Res Commun. 2004;322:1384–91.
Vergarajauregui S, Connelly PS, Daniels MP, et al. Autophagic dysfunction in mucolipidosis type IV patients. Hum Mol Genet. 2008;17:2723–37.
Mizushima N, Yoshimori T, Levine B. Methods in mammalian autophagy research. Cell. 2010;140:313–26.
Mizushima N, Yoshimori T, Ohsumi Y. The role of Atg proteins in autophagosome formation. Annu Rev Cell Dev Biol. 2011;27:107–32.
Napolitano G, Ballabio A. TFEB at a glance. J Cell Sci. 2016;129:2475–81.
Wong C-O, Li R, Montell C, et al. Drosophila TRPML Is required for TORC1 activation. Curr Biol. 2012;22:1616–21.
Onyenwoke RU, Sexton JZ, Yan F, et al. The mucolipidosis IV Ca2+ channel TRPML1 (MCOLN1) is regulated by the TOR kinase. Biochem J. 2015;470:331–42.
Li R-J, Xu J, Fu C, et al. Regulation of mTORC1 by lysosomal calcium and calmodulin. eLife 2016;5:487.
Sun X, Yang Y, Zhong XZ, et al. A negative feedback regulation of MTORC1 activity by the lysosomal Ca2+ channel MCOLN1 (mucolipin 1) using a CALM (calmodulin)-dependent mechanism. Autophagy. 2018;14:38–52.
Weinstock LD, Furness AM, Herron SS, et al. Fingolimod phosphate inhibits astrocyte inflammatory activity in mucolipidosis IV. Hum Mol Genet. 2018;27:2725–38.
Holers VM. Complement and its receptors: new insights into human disease. Annu Rev Immunol. 2014;32:433–59.
Ko EM, Ma JH, Guo F, et al. Deletion of astroglial CXCL10 delays clinical onset but does not affect progressive axon loss in a murine autoimmune multiple sclerosis model. J Neuroinflammation. 2014;11:1–11.
Keren-Shaul H, Spinrad A, Weiner A, et al. A unique microglia type associated with restricting development of Alzheimer’s disease. Cell. 2017;169(1276–1290):e17.
Canales NAG, Marina VM, Castro JS, et al. A1BG and C3 are overexpressed in patients with cervical intraepithelial neoplasia III. Oncol Lett. 2014;8:939–47.
Piyaphanee N, Ma Q, Kremen O, et al. Discovery and initial validation of α 1‐B glycoprotein fragmentation as a differential urinary biomarker in pediatric steroid‐resistant nephrotic syndrome. Mischak H, Thongboonkerd V, Schanstra JP, et al., eds. Proteomics Clin Appl 2011;5:334–342.
Sugihara T, Tanaka S, Braga-Tanaka I III, et al. Screening of biomarkers for liver adenoma in low-dose-rate γ-ray-irradiated mice. Int J Radiat Biol. 2018;94:315–26.
Kiyosawa K, Mokuno K, Murakami N, et al. Cerebrospinal fluid 28-kDa calbindin-D as a possible marker for Purkinje cell damage. J Neurol Sci. 1993;118:29–33.
Bradbury A, Bagel J, Sampson M, et al. Cerebrospinal fluid calbindin D concentration as a biomarker of cerebellar disease progression in Niemann-pick type C1 disease. J Pharmacol Exp Ther. 2016;358:254–61.
Sleat DE, Tannous A, Sohar I, et al. Proteomic analysis of brain and cerebrospinal fluid from the three major forms of neuronal ceroid lipofuscinosis reveals potential biomarkers. J Proteome Res. 2017;16:3787–804.
Schiffmann R, Mayfield J, Swift C, et al. Quantitative neuroimaging in mucolipidosis type IV. Mol Genet Metab. 2014;111:147–51.
Rosato AS, Montefusco S, Soldati C, et al. TRPML1 links lysosomal calcium to autophagosome biogenesis through the activation of the CaMKKβ/VPS34 pathway. Nat Commun. 2019;10:1–16.
Jennings JJ, Zhu J-H, Rbaibi Y, et al. Mitochondrial aberrations in mucolipidosis Type IV. J Biol Chem. 2006;281:39041–50.
Kiselyov K, Jennings JJ Jr, Rbaibi Y, et al. Autophagy, mitochondria and cell death in lysosomal storage diseases. Autophagy. 2007;3:259–62.
Bozzato A, Barlati S, Borsani G. Gene expression profiling of mucolipidosis type IV fibroblasts reveals deregulation of genes with relevant functions in lysosome physiology. Biochim Biophys Acta. 2008;1782:250–8.
Cougnoux A, Drummond RA, Fellmeth M, et al. Unique molecular signature in mucolipidosis type IV microglia. J Neuroinflammation. 2019;16:1–12.
Colletti GA, Miedel MT, Quinn J, et al. Loss of lysosomal ion channel transient receptor potential channel mucolipin-1 (TRPML1) leads to cathepsin B-dependent apoptosis. J Biol Chem. 2012;287:8082–91.
Chandra M, Zhou H, Li Q, et al. A role for the Ca2+ channel TRPML1 in gastric acid secretion, based on analysis of knockout mice. Gastroenterology. 2011;140(857–867):e1.
Dassati S, Waldner A, Schweigreiter R. Apolipoprotein D takes center stage in the stress response of the aging and degenerative brain. Neurobiol Aging. 2014;35:1632–42.
Terrisse L, Poirier J, Bertrand P, et al. Increased Levels of Apolipoprotein D in Cerebrospinal Fluid and Hippocampus of Alzheimer’s Patients. J Neurochem. 1998;71:1643–50.
Reindl M, Knipping G, Wicher I, et al. Increased intrathecal production of apolipoprotein D in multiple sclerosis. J Neuroimmunol. 2001;119:327–32.
Kalman K, McConathy W, Araoz C, et al. Apolipoprotein D in the aging brain and in Alzheimer’s dementia. Neurol Res. 2016;22:330–6.
Suresh S, Yan Z, Patel RC, et al. Cellular cholesterol storage in the niemann-pick disease type C mouse is associated with increased expression and defective processing of apolipoprotein D. J Neurochem. 1998;70:242–51.
Li L, Liu M-S, Li G-Q, et al. Relationship between apolipoprotein superfamily and Parkinson’s disease. Chin Med J. 2017;130:2616–23.
Van Giau V, Bagyinszky E, An SSA, et al. Role of apolipoprotein E in neurodegenerative diseases. Neuropsychiatr Dis Treat. 2015;11:1723–37.
Vardi A, Ben-Dor S, Cho SM, et al. Mice defective in interferon signaling help distinguish between primary and secondary pathological pathways in a mouse model of neuronal forms of Gaucher disease. J Neuroinflammation. 2020;17:1–13.
Blumenreich S, Yaacobi C, Vardi A, et al. Substrate reduction therapy using Genz‐667161 reduces levels of pathogenic components in a mouse model of neuronopathic forms of Gaucher disease. J Neurochem 2020:jnc.15136.
Han X, Zhu J, Zhang X, et al. Plin4-dependent lipid droplets hamper neuronal mitophagy in the MPTP/p-induced mouse model of Parkinson’s disease. Front Neurosci. 2018;12:1257.
Wang C-W. Lipid droplets, lipophagy, and beyond. Biochim Biophys Acta. 2016;1861:793–805.
Pryor PR, Reimann F, Gribble FM, et al. Mucolipin-1 Is a Lysosomal Membrane Protein Required for Intracellular Lactosylceramide Traffic. Traffic. 2006;7:1388–98.
Shen D, Wang X, Li X, et al. Lipid storage disorders block lysosomal trafficking by inhibiting a TRP channel and lysosomal calcium release. Nat Commun. 2012;3:1–11.
Curcio Morelli C, Charles FA, Micsenyi MC, et al. Macroautophagy is defective in mucolipin-1-deficient mouse neurons. Neurobiol Dis. 2010;40:370–7.
García-Arencibia M, Hochfeld WE, Toh PPC, et al. Autophagy, a guardian against neurodegeneration. Semin Cell Dev Biol. 2010;21:691–8.
Vitner EB, Farfel-Becker T, Ferreira NS, et al. Induction of the type I interferon response in neurological forms of Gaucher disease. J Neuroinflammation. 2016;13:1–15.
Perez-Riverol Y, Csordas A, Bai J, et al. The PRIDE database and related tools and resources in 2019: improving support for quantification data. Nucleic Acids Res. 2019;47:D442–50.
Wyant GA, Abu-Remaileh M, Frenkel EM, et al. NUFIP1 is a ribosome receptor for starvation-induced ribophagy. Science. 2018;360:751–8.
Chacon J, Rosas L, Cuajungco MP. ZnT3 expression levels are down-regulated in the brain of Mcoln1 knockout mice. Mol Brain. 2019;12:24.
We thank Dr. Yishai Levin and Dr. Tali Shalit from The De Botton Protein Profiling Institute and The Mantoux Bioinformatics Institute of the Nancy and Stephen Grand Israel National Center for Personalized Medicine, Weizmann Institute of Science, for help with proteomics, and Dr. Soo Min Cho for help with the initials stages of this work. We thank Dr. Nick Platt for fruitful discussions. We thank Rebecca Oberman for her contribution to the MLIV community. We acknowledge the NIH NeuroBioBank for providing brain tissue and CSF for this study. A.H. Futerman is the Joseph Meyerhoff Professor of Biochemistry at the Weizmann Institute of Science.
This work was supported by the MLIV Foundation.
The Maryland Repository's consent process, documents, and collection procedures are approved annually by the University of Maryland Institutional Review Board. Mice were maintained under specific pathogen-free conditions and handled according to protocols approved by the Weizmann Institute Animal Care Committee according to international guidelines.
Consent for publication
All authors have read and agreed to the published version of the manuscript.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Vardi, A., Pri-Or, A., Wigoda, N. et al. Proteomics analysis of a human brain sample from a mucolipidosis type IV patient reveals pathophysiological pathways. Orphanet J Rare Dis 16, 39 (2021). https://doi.org/10.1186/s13023-021-01679-7