Neural Plasticity
Volume 2019 (2019), Article ID 9653024, 12 pages
https://doi.org/10.1155/2019/9653024

Research Article

Functional Role of SIL1 in Neurodevelopment and Learning

Shilian Xu,¹ Jia Zhu,² Kai Mi,¹ Yan Shen,¹ and Xiaomin Zhang¹

¹Department of Physiology, School of Basic Medicine, Kunming Medical University, Kunming, Yunnan Province, China
²Department of Pathology, Institute of Diabetes and Nerve Disease, Jiaxing Hospital of Traditional Chinese Medicine, Jiaxing University, Jiaxing, Zhejiang Province, China

Correspondence should be addressed to Xiaomin Zhang; zhangxmtan@foxmail.com

Received 6 February 2019; Revised 2 June 2019; Accepted 1 July 2019; Published 18 August 2019

Academic Editor: Long-Jun Wu

Copyright © 2019 Shilian Xu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Background. Sil1 is the causative gene of Marinesco-Sjӧgren Syndrome (MSS). The mutated Sil1 generates shortened SIL1 protein which will form aggregation and be degraded rapidly. Mental retardation is a major symptom of MSS which suggests a role of SIL1 in the development of the central nervous system, but how SIL1 functions remains unclear. Objectives. The aim of this study is to explore the role of SIL1 in regulating cerebral development and its underlying molecular mechanism. Methods. The basic expression pattern of SIL1 in tissues and cultured cortical neurons is measured by immunostaining and Western blot. The expression of SIL1 is reduced in vitro and in vivo through RNA interference delivered by a lentivirus. The expression of NMDA receptor subunits and the function of the Reelin signaling pathway are then examined by surface biotinylation and Western blot subsequently. Finally, the spatial learning of young mice was assessed by the Barnes maze task. Results. SIL1 deficiency caused a diminished expression of both Reelin receptors and therefore impaired the Reelin signaling pathway. It then inhibited the developmental expression of GluN2A and impaired the spatial learning of 5-week-old mice. Conclusions. These results suggested that SIL1 is required for the development of the central nervous system which is associated with its role in Reelin signaling.

1. Introduction

SIL1 is an Endoplasmic Reticulum- (ER-) resident 54 kD protein that is composed of 461 amino acids [1]. SIL1 has an ER-targeting sequence in its amino terminus and an ER retention KDEL sequence in its carboxyl terminus [2]. SIL1 is the mammalian homolog of yeast Sls1p/Sil1p and functions as the nucleotide exchange factor of ER chaperone protein Bip [3]. It is well acknowledged that Bip is a member of the heat shock protein 70 family and it plays important roles in mediating folding and assembly of nascent proteins, as well as degradation of misfolded proteins [4–6]. Both binding and separating of Bip from its substrate protein require the assistance of cofactors [2]. As an adenine nucleotide exchange factor, SIL1 regulates the ATPase activity of Bip and promotes the release of the substrate protein [3, 7].

In 2005, two groups independently identified SIL1 (gene ID: 64374) as the causative gene of Marinesco-Sjӧgren Syndrome (MSS; OMIM 248800) [8, 9]. Single-gene mutation of SIL1 is enough to cause MSS. MSS is an autosomal recessive multisystem disorder, and its main symptoms include cerebellar ataxia, cataracts, mental retardation, myopathy, and short stature [9–11]. Different kinds of SIL1 mutation have been discovered in MSS patients, including missense mutation, in-frame deletion, and several single-nucleotide mutations that affect RNA-splicing sites [11–14]. These mutations cause codon shift or deletion inside exon 6 and exon 9 which result in abnormal expression of the protein [15]. As a consequence, the mutated SIL1 that could not bind with Bip or be stably retained in ER would be transported into the cytoplasm and degraded by the proteasome [12, 16].

As a cofactor of Bip, SIL1 expresses in all types of cells; however, only several organs are affected in MSS, especially the central nervous system. 90% of the MSS patients showed moderate to severe mental retardation [11], which indicates that SIL1 may play specific roles in the nervous system. It has been proved that the SIL1-mutated woozy mice developed ataxia which resulted from Purkinje cell loss in the anterior cerebellar lobules [17]. A recent research found that the migration and morphological maturation of cortical neurons during development are impaired after RNAi mediated gene silencing of SIL1. However, the cortical localization of neurons in adult mice was normal which suggested that the migration was only delayed but not irreversibly inhibited. These studies showed that SIL1 plays important roles in the development of the mammalian cortex and cerebellum, but the mechanism is still not clear yet. Furthermore, current evidence could not directly prove that these developmental defects caused by SIL1 deficiency are responsible for the intelligent disability of MSS patients. In our study, we aimed to provide new insights regarding these issues.

In our former study, we found that Bip was involved in the dynamic trafficking of the GluN2A-containing NMDA receptor induced by plasticity stimulation and the subunit-selective interaction between Bip and GluN2A was required for fear conditioning [18]. The NMDA receptor is a crucial excitatory receptor in the central nervous system which is a heterogeneous tetramer and has several subtypes. During early development, the GluN2B-containing NMDA receptor is the dominant subtype in the frontal cortex and then GluN2A increases extensively after birth and becomes the predominant subunit. This GluN2B to GluN2A dominancy shift is crucial for the intellectual development of mammals. Since SIL1 is the cofactor of Bip, we speculated that SIL1 may also be involved in the trafficking of GluN2A. Therefore, we inhibited SIL1 expression through RNA interference both in vitro and in vivo, and we found that the postnatal expression of the GluN2A-containing NMDA receptor was inhibited and the spatial learning of mice was impaired, both of which are caused by the impaired Reelin signaling pathway. Therefore, our data proved for the first time that SIL1 regulates the expression of GluN2A through the extracellular matrix-dependent signaling pathway which is essential for the neurodevelopment and spatial memory formation.

2. Materials and Methods

2.1. Animals

Male and female C57 mice were obtained and acclimated to the colony room for at least 2 weeks prior to mating. The colony room was maintained on a 14 hr light/10 h dark cycle with controlled temperature and humidity with food and water available freely. The subjects were the offspring from the pairing of a female with a male partner. Virus was injected at 3 weeks old and behavioral task was performed 2 weeks later.

2.2. Cortical Neuron Culture

For the primary neuron culture, we followed the methods of Zhang et al. [18]. Cortical neuronal cultures were prepared from one-day-old mice. Briefly, the cortex tissue was extracted and chopped to be size and then digested in 0.28% trypsin (Sigma, St. Louis, MO) for 15 min at 37°C with gentle shaking. Then, the sample was rinsed twice with ECS, and then the dissociated cells were plated at a density of in a 35 mm (for staining) dish on poly-L-lysine-coated coverslips or directly in a 60 mm dish. 2 ml or 4 ml of neurobasal plus medium (Invitrogen, Carlsbad, CA) containing 10% FBS, 2 mM glutamine, and 100 μg/ml gentamycin (all from Invitrogen Corporation, Carlsbad, CA) was added, and then cells were maintained at 37°C in 5% CO₂. After 24 hr, the medium was replaced with neurobasal medium containing 2% B27 supplement, 1% antibiotic, and 0.25% glutamine (Invitrogen). At DIV5, cytosine arabinofuranoside was added at a final concentration of 10 μM. Thereafter, half of the medium was replaced twice a week.

2.3. RNA Interference

The lentivirus of shRNA was designed to target 2 distinct coding sequences in mouse SIL1 according to a former literature [19] (shRNA#1: 5-GCTCCAACAAGAAGACAAA-3; shRNA#2: 5GGTTGCTGCGCTCTTTGAT-3) (constructed by GenePharma, Suzhou, China). For the cultured cortical neurons, shRNA#1, shRNA#2, or scramble shRNA was added into serum-free fresh culture medium at DIV6, and 48 hr after incubation with virus, the culture medium collected before infection was replaced. For in vivo experiments, shRNA#1 or shRNA#2 was delivered through bilateral intraventricular injection to 3-week-old mice and the behavioral tasks were performed 2 weeks later.

2.4. Cell Staining

For the neuron staining, we followed the methods of Zhang et al. [18]. For colocalization, cortical neurons on coverslips were fixed with 4% paraformaldehyde for 10 min after a brief rinse in prewarmed ECS at room temperature, then permeabilized, and blocked through incubating in PBS containing 0.5% Triton X-100 and 5% BSA for 0.5 hr. Then, the neurons were incubated with primary antibodies to SIL1 (rabbit anti-SIL1, Abcam) and synaptophysin (mouse anti-synaptophysin, Abcam) or SD95 (mouse anti-PSD95, Abcam) in PBS containing 0.5% Triton X-100 overnight at 4°C. After rinsing in PBS containing 0.5% Triton X-100 3 times, neurons were incubated with both Alexa 488-conjugated secondary antibody and Alexa 546-conjugated secondary antibody (donkey anti-rabbit or mouse secondary antibody, Abcam) for 1 h at room temperature. After rinsing with PBS three times, neurons were examined under a 60x, 1.4 numerical aperture oil-immersion objective on an Olympus confocal microscope.

2.5. PSD Fraction Separation

For the Western blot analysis of the PSD fraction, we followed the methods of Zhang et al. [18]. After rinsing twice with prewarmed ECS, the DIV14 neurons were homogenized in buffer A (320 mM sucrose, 10 mM HEPES, pH 7.4) through ultrasonication and centrifuged at 12 000×g for 20 min and the pellet was collected as a crude membrane fraction. The pellet was then dissolved in buffer B (4 mM HEPES, 1 mM EDTA, pH 7.4) and centrifuged at 12 000×g for another 20 min. After centrifugation, the pellet was further dissolved in buffer C (20 mM HEPES, 100 mM NaCl, 0.5% Triton X-100, pH 7.2) and incubated for 15 min at 4°C. The homogenate was then centrifuged at 12 000×g for 20 min, and the pellet was the Triton X-100-insoluble PSD fraction. This pellet was dissolved in buffer D (20 mM HEPES, 0.15 mM NaCl, 1% Triton X-100, 1% deoxycholic acid, 1% SDS, 1 mM DTT, pH 7.5) and incubated at 4°C for 1 hr. After centrifugation at 10 000×g for 15 min, a 4x sample buffer was added to the supernatant and the sample was boiled at 100°C for 10 min for Western blot.

2.6. Protein Extraction

For the protein extraction, we followed the methods of Zhang et al. [18]. DIV14 cortical cultures were rinsed three times in ECS, and then 500 μl RIPA buffer was added into each 60 mm dish, and the neurons were then collected and ultrasonicated twice for 8 sec each followed by incubation at 4°C for 2 hr. Homogenized neurons were then centrifuged at 12 000×g for 10 min at 4°C, and a 4x sample buffer was added to the supernatant and boiled at 100°C for 10 min, then used for SDS-PAGE. Adult mouse brain tissue was homogenized in RIPA buffer and then centrifuged at 700×g for 10 min at 4°C; the supernatant was incubated at 4°C for 2 hr and then centrifuged at 12 000×g for 10 min at 4°C. The supernatant was then mixed with a 4x sample buffer and boiled at 100°C for 10 min for Western blot.

2.7. Coimmunoprecipitation

For the coimmunoprecipitation, we followed the methods of Zhang et al. [18]. Briefly, after protein extraction from mouse cortical tissue, the protein concentration of the supernatant was adjusted to 1 μg/μl. 10 μg of the corresponding antibody was added to 200 μl of the sample and incubated overnight at 4°C. 45 μl of the sample was collected as total protein (input). Then, 30 μl of solubilized protein A-sepharose beads was added and incubated with the sample for another 2 hr at 4°C. Then, the mixture was rinsed three times with binding buffer and once with binding buffer containing 500 mM NaCl to remove nonspecific proteins. After final centrifugation, the supernatant was collected and the pellet was incubated with a 2x sample buffer and incubated at 90°C for 10 min. The input, supernatant, and pellet were then used for SDS-PAGE. For the input and supernatant, 15 μl of 4x sample buffer was added to 45 μl of the sample and then incubated at 90°C for 10 min before being used for Western blot.

2.8. cLTP Stimulation

For the cLTP stimulation, we followed the methods of Zhang et al. [18]. DIV14 cultured cortical neurons were rinsed twice in prewarmed ECS and then incubated in buffer A (0.5 μM TTX, 1 μM strychnine, 20 μM bicuculline in ECS, pH 7.4) for 5 min and then in buffer B (200 μM glycine in buffer A) for another 5 min. After rinsing in ECS for 5 min once, neurons were incubated in buffer A for another 20 min. For the control treatment, buffer B was replaced by ECS. The protein sample was then prepared as described in Protein Extraction and used for Western blot.

2.9. Biotinylation

For the biotinylation experiment, we followed the methods of Zhang et al. [18]. DIV14 cultured cortical neurons in a 60 mm culture dish were rinsed twice using ice-cold PBS with 1 mM MgCl₂ and 0.5 mM CaCl₂ and then incubated with PBS with 1 mg/ml sulfo-NHS-LC-biotin (Thermo Fisher Scientific) for 30 min at 4°C. 500 μl PBS with 100 μM glycine was then added to each 60 mm dish to quench the biotin binding for 15 min at 4°C. Neurons were then lysed in ice-cold RIPA buffer (0.1% SDS, 0.5% sodium deoxycholate, and 1% Triton X-100) after rinsing twice with ice-cold PBS and incubated at 4°C for another 30 min. Cells were concentrated at 12 000×g for 10 min, and the supernatant was incubated with NeutrAvidin beads (Thermo Fisher Scientific) for 2 hr at 4°C. After incubation, the beads were then washed three times with RIPA buffer and proteins were mixed with a 2x sample buffer and boiled at 100°C for 10 min before Western blot.

2.10. Barnes Maze Task

For the Barnes maze task, the mice have one habituation trial where it is guided to the shelter manually by the experimenter immediately after being positioned on the platform. This is followed by a 4-day acquisition period, where each day, the mice were allowed to navigate freely in the maze for three minutes, after which it is again manually guided to the shelter if it does not reach it during the exploration period. The last stage of the assay is one probe trial on the 5th day where the animal explores the maze with all holes closed. Learning and long-term memory are quantified by the time taken to find the shelter during the acquisition period and by the time spent near the closed shelter hole in the probe trial. Direct intraventricular injection of lentivirus of shRNA#1, shRNA#2, or scramble shRNA was performed when mice were 3 weeks old, and the Barnes maze task was performed 4 weeks later.

2.11. Statistical Analysis

All data were analyzed using GraphPad prism 5 software and are represented as . Data from multiple groups were quantified using one-way ANOVA, and comparisons of two groups were quantified by unpaired -test.

3. Results

3.1. Result 1: The Expression Profile of SIL1 In Vivo and In Vitro

In our previous study, we found that Bip is required for the dynamic synaptic insertion of the GluN2A-containing NMDA receptor which is crucial for fear conditioning [18]. SIL1 is a cofactor of Bip and has been found to play important roles in the central nervous system [9, 19, 20]. Therefore, in this study, we wanted to investigate the potential role of SIL1 in regulating GluN2A trafficking. We first examined the developmental expression pattern of SIL1 in the cortex during development and its expression in adult mouse tissues. In consistency with former studies [19], we found that the level of SIL1 was very low during early development in mouse cortical tissues (Figure 1(a)) and this is different from Bip which is highly expressed throughout development [21]. Furthermore, we found that the expression pattern of SIL1 was very similar with GluN2A that both proteins underwent a sudden increase between postnatal 7 and 10 days (Figure 1(a)). Thereafter, SIL1 maintained a relatively high level in the brain especially in the hippocampus and cerebellum compared to other tissues (Figure 1(b)). It is well acknowledged that the NMDA receptor has a subtype switch from GluN2B dominant to GluN2A dominant in early development during which the level of GluN2B maintains constant while GluN2A increases rapidly, and it is necessary for the proper function of the central nervous system. Therefore, this result suggested that SIL1 may participate in the developmental subtype switch of the NMDA receptor.

(a)

(b)

Figure 1: Expression of SIL1 in tissues. (a) Whole lysate of the cerebral cortex at various developmental stages (every 4 days since postnatal day 1 and the last sample was the maternal mice) was subjected to Western blot, and the expression of SIL1, GluN2A, GluN2B, and GAPDH was examined. (b) Whole lysate of adult mouse tissues was examined for the SIL1 expression, including the cerebellum, olfactory bulb, cortex, hippocampus, heart, kidney, lung, and liver.

Next, we examined the distribution of SIL1 in DIV14 cultured cortical neurons. The rabbit anti-SIL1 primary antibody and Alexa 488 conjunct anti-rabbit secondary antibody were used to show the distribution of SIL1, and mouse anti-PSD95 or synaptophysin and Alexa 594 conjunct anti-mouse secondary antibody were used to show PSD95 or synaptophysin. We found SIL1 to be extensively expressed in neurons; strong fluorescence could be detected both in soma and in extrusions. Besides, the partial colocalization of SIL1 with PSD95 suggested a localization of SIL1 in the distal dendrite (Figure 2(a)), and therefore, we further evaluated whether SIL1 is expressed within the synapse. We separated the extrasynaptic membrane fraction (S1), Triton X-100-soluble synaptic faction (S2), and Triton X-100-insoluble PSD fraction (S3) from the adult mouse cortex. In consistency with our former study [18], we found Bip to be mainly accumulated at the extrasynaptic membrane fraction and there was also a moderate portion of Bip in the PSD fraction (Figure 2(b)). However, the distribution of SIL1 was different from Bip but similar with GluN2A and GluN2B in which SIL1 mainly accumulated in the PSD fraction while a small amount could also be seen in the Triton X-100-soluble synaptic fraction (Figure 2(b)). Because SIL1 is an ER-resident protein, we assume that its attaching to the postsynaptic membrane depends on its binding with other synaptic proteins; thus, we evaluated the interaction of SIL1 with PSD95 and GluN2A through coimmunoprecipitation (CO-IP). We found mutual interaction between SIL1 and PSD95, as well as GluN2A (Figure 2(c)). Moreover, SIL1 had a stronger interaction with PSD95 compared to GluN2A (Figure 2(c)) while it had no interaction with GluN2B (Figure 2(d)). These results showed that SIL1 is distributed broadly in cortical neurons, and similar with GluN2A, a large portion of SIL1 is located in the postsynaptic compartment. Altogether, these data suggested a relationship between SIL1 and GluN2A, and a possible role of SIL1 is in the trafficking of GluN2A.

(a)

(b)

(c)

(d)

Figure 2: Expression of SIL1 in cultured cortical neurons. (a) Cellular localization of SIL1 was examined in cultured cortical neurons at DIV14. SIL1 was immunostained with anti-SIL1 antibody (green), and synaptophysin or PSD95 was immunostained with anti-synaptophysin or anti-PSD95 antibody (red). Scale bar equals to 10 μm. (b) Expression of SIL1 in different subneuronal compartments of the adult mouse cortex tissue was investigated. S1 indicated extrasynaptic membrane fraction, S2 indicated Triton X-100-soluble synaptic faction, and S3 indicated Triton X-100-insoluble PSD fraction. (c, d) Interactions between SIL1 and GluN2A, PSD95, Bip, and GluN2B were examined by coimmunoprecipitation. Anti-SIL1, anti-PSD95, anti-GluN2A, and anti-GluN2B antibodies were used to precipitate the other proteins. GAPDH was used as a negative control.

3.2. Result 2: Silencing of SIL1 Inhibited the Developmental Expression of GluN2A but Not Plasticity-Induced Transient Trafficking of GluN2A

Since we found that the developmental expression of SIL1 was very similar to that of GluN2A, and SIL1 interacted with GluN2A in cultured cortical neurons, we next investigated whether SIL1 was required for the normal expression of GluN2A. Therefore, we used shRNA to suppress the expression of SIL1. We designed two lentivirus-encoding shRNA targeting two distinct sequences of mouse SIL1 (shRNA#1 and shRNA#2). The cultured cortical neurons were incubated with either shRNA#1 or shRNA#2 for 48 hr from DIV6 and further cultured until DIV14, and then the protein level was evaluated through Western blot. Meanwhile, in the control group, the neurons were incubated with scramble shRNA. We found that both shRNA#1 and shRNA#2 effectively inhibited the expression of SIL1 (Figures 3(a) and 3(b)), as well as GluN2A, while GluN2B was not affected (Figures 3(a) and 3(b)). Next, we measured the membrane expression of GluN2A and GluN2B through surface biotinylation, and consistent with a former result, membrane expression of GluN2A, but not GluN2B, was reduced significantly after SIL1 silencing (Figures 3(c) and 3(d)). These results suggested that SIL1 is required for the normal developmental expression of GluN2A.

(a)

(b)

(c)

(d)

Figure 3: Expression of GluN2A was diminished after SIL1 silencing. (a) Two silencing vectors (shRNA#1 and shRNA#2) targeting SIL1 were designed and transfected into cultured cortical neurons at DIV6 through lentivirus transfection for two-day incubation, and whole cell lysates were harvested at DIV14. The scramble shRNA was used as the control. The expression of SIL1, GluN2A, GluN2B, and Bip was examined by Western blot. (b) Statistical analysis of (a): the expression of each protein was first normalized to GAPDH and then compared with the control. The levels of SIL1 and GluN2A were significantly reduced after shRNA#1 and shRNA#2 were reduced compared to the scramble shRNA-treated control group (SIL1: after shRNA#1 treatment, after shRNA#2 treatment; GluN2A: after shRNA#1 treatment, after shRNA#2 treatment). The levels of GluN2B and Bip were not altered (GluN2B: after shRNA#1 treatment, after shRNA#2 treatment; Bip: after shRNA#1 treatment, after shRNA#2 treatment). (c) The membrane expression of SIL1, GluN2A, GluN2B, and Bip was measured by biotinylation followed by Western blot (surface: membrane fraction; offer: total protein). (d) Statistical analysis of (c): the expression of each protein was first normalized to GAPDH and then compared with the control (GluN2A: after shRNA#1 treatment, after shRNA#2 treatment; GluN2B: after shRNA#1 treatment, after shRNA#2 treatment). All data were presented as . ##. for all treatments.

In our former study, we found that Bip regulates the dynamic synaptic insertion of the GluN2A-containing NMDA receptor after receiving plastic stimulation. Since SIL1 is a cofactor of Bip, we further investigated whether SIL1 also regulates the transient trafficking of GluN2A. Thus, we applied chemical long-term potentiation (cLTP) stimuli to the cortical neurons pretreated with silencing shRNA or scramble shRNA and then measured the dynamic membrane insertion of the GluN2A-containing NMDA receptor. In scramble shRNA treatment neurons, the membrane GluN2A was quickly elevated after cLTP as previously reported, while the membrane expression of GluN2A in shRNA#2-treated neurons was significantly lower than scramble shRNA-treated neurons (Figures 4(a) and 4(b)) but was comparable to the neurons without cLTP stimulation (Figures 4(a) and 4(b)). Although we could not directly compare the total and surface expression of GluN2A, it is obvious that the total level of GluN2A in the shRNA#2-treated group is lower than that in the shRNA#2-untreated neurons but the surface fraction of GluN2A in the shRNA#2-treated group is comparable to that in shRNA#2-untreated neurons (Figures 4(a) and 4(b)). Therefore, from this result, we can infer that the dynamic trafficking of GluN2A induced by cLTP stimulation was not fully abolished by SIL1 silencing and the seemingly reduced surface expression of GluN2A was caused by the decrease of overall expression of the protein. Altogether, these results suggested that SIL1 is required for the basic developmental expression of GluN2A but is not necessary for the transient dynamic trafficking of GluN2A in mature neurons.

(a)

(b)

Figure 4: Dynamic trafficking of GluN2A was not affected after SIL1 silencing. (a) Neurons transfected with silencing shRNA or scramble shRNA at DIV6 were subjected to chemical LTP (cLTP) at DIV14 (cLTP+shRNA#2 and cLTP+scramble). Immediately after cLTP treatment, the membrane expression (surface) and total protein level (offer) were assessed by surface biotinylation. Neurons transfected with scramble shRNA without cLTP stimulation were used as the control (con). (b) Statistical analysis of (a): the expression of each protein was first normalized to GAPDH and then compared with the control. The expression of GluN2A reduced after SIL1 silencing (, ) while the expression of GluN2B was not altered (, ). In the scramble shRNA-incubated neuron, GluN2A increased after cLTP stimulation (, ), while in shRNA#2 neurons, GluN2A was also increased after cLTP ( compared to , , ) which indicated that the dynamic membrane insertion of GluN2A was not impaired after SIL1 silencing. All data were shown as . ##. for all treatments.

3.3. Result 3: Silencing of SIL1 Impaired the Signaling of the Reelin Signaling Pathway

The above data suggested a novel role of SIL1 in regulating the developmental expression of GluN2A. But it seems that SIL1 was not directly involved in GluN2A maturation and assembly, because if this is the case, then the cLTP-induced GluN2A membrane insertion should also be inhibited after SIL1 silencing. Therefore, we assessed an alternative mechanism that may explain this phenomenon. It has been proved that the Reelin signaling pathway is important for the NMDA receptor composition change during neuron maturation [22, 23, 40]. And also, it was found that the deficiency of cortical neuron migration caused by SIL1 mutation is very similar with that of Reelin knockout mice [19, 24]. Thus, we speculate that silencing of SIL1 might have impaired the Reelin signaling pathway and in turn affected the expression of GluN2A. Reelin is an extracellular matrix protein and functions on neurons through both its receptors, VLDLR and ApoER2 [25]. Therefore, we first examined whether SIL1 was involved in the expression of VLDLR and ApoER2. Through coimmunoprecipitation, we identified a strong interaction between SIL1 and VLDLR as well as ApoER2 (Figure 5(a)). Furthermore, after SIL1 silencing, both the total amount and surface expression of VLDLR and ApoER2 were strongly reduced (Figures 5(b) and 5(c)). These results suggested that the expression of the two main receptors of Reelin signaling depended on SIL1. Because previous studies have shown that the activation of the Reelin signaling pathway is required for the developmental expression of GluN2A [23, 40], we further measured the phosphorylation of Dab1, which is a prominent hub molecule of Reelin signaling and has been known to regulate synaptic maturation [40]. As a result, the phosphorylation level but not the total amount of Dab1 significantly decreased after silencing of SIL1 (Figures 5(d) and 5(e)). The above data showed that SIL1 was required for the expression of both Reelin receptors and therefore played vital roles in the activation of the Reelin signaling pathway.

(a)

(b)

(c)

(d)

(e)

Figure 5: Reelin signaling was impaired after SIL1 silencing. (a) Interactions between SIL1 and PSD95, VLDLR, and ApoER2 were examined by coimmunoprecipitation. Anti-SIL1, anti-PSD95, anti-VLDLR, and anti-ApoER2 antibodies were used to coprecipitate SIL1. GAPDH was used as a negative control. (b) The membrane fraction and total expression of VLDLR and ApoER2 were examined after incubation of neurons with silencing shRNA (shRNA#1 and shRNA#2) or scramble shRNA (con). (c) Statistical analysis of (b): the expression of each protein was first normalized to GAPDH and then compared with the control; both total expression and surface fraction of ApoER2 and VLDLR were significantly reduced (ApoER2 total: after shRNA#1 treatment, after shRNA#2 treatment; ApoER2 surface: after shRNA#1 treatment, after shRNA#2 treatment; VLDLR total: after shRNA#1 treatment, after shRNA#2 treatment; and VLDLR surface: after shRNA#1 treatment, after shRNA#2 treatment). All data were shown as . ##. . (d) The phosphorylation of Dab1 after SIL1 silencing (shRNA#1 and shRNA#2) was assessed by antibody targeting phosphorylated Dab1. Neurons treated with scramble shRNA were used as the control (con). (e) Statistical analysis of (d): the phosphorylated Dab1 was first normalized to GAPDH and then to Dab1 and compared with the control (p-Dab1: after shRNA#1 treatment, after shRNA#2 treatment). All data were shown as . The phosphorylation of Dab1 was significantly inhibited by SIL1 silencing. ##. for all treatments.

3.4. Result 4: Silencing of SIL1 Impaired the Spatial Memory Which Was Rescued by Mimic Phosphorylation Dab1 Peptide

Sil1 is the only gene found to be responsible for the MSS which has a major symptom of mental retardation. Previous studies have shown that SIL1 silencing caused a delayed migration of neurons during neocortex development which is similar to Reelin knockout mice [19]. Our data also showed that SIL1 silencing resulted in an inactivation of Reelin signaling and obstructed synaptic NMDA receptor maturation. Other studies have shown that Reelin could modulate learning and memory by modulating NMDA receptor expression and function [40]. These evidence supported that the mental retardation caused by mutation of Sil1 may associate with malfunction of Reelin signaling and the subsequent abnormal expression of the NMDA receptor. Therefore, we examined the spatial memory formation of young mice after SIL1 silencing. The lentivirus-encoding shRNA#2 or scramble shRNA was administered to 3-week-old mice through bilateral ventricular injection. And the Barnes Maze task was performed 2 weeks later to examine the spatial learning of the mice. And we found that the escape latency of the SIL1 silencing mice was significantly longer than that of the control mice in the acquisition period since day 2 (Figure 6(a)) and the time that mice spent in the target quarter as well as the percentage of time in the target quarter in the probe trial is significantly reduced (Figures 6(b) and 6(c)). Since MSS patients show early onset cataract and muscle atrophy, we evaluated the performance of SIL1-silenced mice in two other behavior tasks, the beam-walking task and the rotarod task. The SIL1-deficient mice showed similar latency time compared to scramble shRNA-treated mice (Supplemental Fig. 1A and 1B), and this result ruled out the possibility that the impaired eyesight or reduced muscle coordination caused the low performance of mice in the Barnes maze task. In order to further prove the involvement of Reelin signaling in SIL1-mediated spatial learning, we used a mimic phosphorylation peptide of Dab1, pTyr-220, which is able to activate the Reelin signaling pathway [26]. We found that daily intraperitoneal injection of pTyr-220 after SIL1 silencing partly rescued the impaired performance of mice in the Barnes maze task (Figures 6(a) and 6(b)). In addition, the expression of GluN2A could also be rescued by pTyr-220 (Supplemental Fig. 2A and 2B). Altogether, these data provided direct evidence of SIL1 functioning in spatial learning which depended on Reelin signaling and was associated with its role in regulating the developmental maturation of the NMDA receptor.

(a)

(b)

(c)

Figure 6: Spatial memory formation of young mice was impaired after SIL1 silencing. (a) Mice were subjected to bilateral intraventrical injection of lentivirus-encoding shRNA#2 at 3 weeks old followed or not followed by daily intraperitoneal injection of mimic phosphorylation peptide of Dab1 (pTyr-220). Barnes maze task was conducted 2 weeks later to assess the spatial learning of mice. The escape latency was measured from day 1 to day 4 of the acquisition period (con: day 1: , day 2: , day 3: , and day 4: , ; shRNA#2: day 1: , day 2: , day 3: , and day 4: , ; and shRNA#2+pTyr-220: day 1: , day 2: , day 3: , and day 4: , ). All data were shown as . # and ## for the shRNA#2 treatment compared to the control; for shRNA#2 compared to shRNA#2+pTyr-220. (b) The time that mice spent in the target quarter was measured in the probe trial (con: , shRNA#2: , and shRNA#2+pTyr-220: , ). ## for shRNA#2 compared to the control and for shRNA#2 compared to shRNA#2+pTyr-220. and . (c) The percentage of time that mice spent in the target quarter was measured in the probe trial (con: , shRNA#2: , and shRNA#2+pTyr-220: ). All data were shown as . ## for shRNA#2 compared to the control and for shRNA#2 compared to shRNA#2+pTyr-220. and .

4. Discussion

Sil1 was identified as the causal gene of Marinesco-Sjӧgren Syndrome by two independent research groups in 2005 [8, 9]. Mental retardation is one of the major symptoms of MSS which indicated that SIL1 might play important roles in neurons. In this study, we intended to explore the molecular mechanism of SIL1 in regulating the development of the central nervous system.

We first detected the expression of SIL1 in brain tissues and neurons. We found that SIL1 expression was quite low during the early development which was consistent with other studies [18]. We also found that SIL1 significantly increased its expression among postnatal 7-10 days which was highly similar with the developmental expression pattern of GluN2A. In our previous study, we found that Bip mediated the subunit-selective increase of the GluN2A-containing NMDA receptor which made it a candidate for regulating the maturational change of the NMDA receptor subtype. However, Bip maintains constantly high expression throughout development which seems contradictory to a role in development and also suggested that other related proteins may be involved. Comparing to Bip, SIL1 possesses more characters that fit this function. Firstly, like GluN2A, SIL1 mainly accumulated in the PSD fraction of neurons and had direct interaction with PSD95. Secondly, SIL1 showed similar developmental expression with GluN2A. As a tetramer, each NMDA receptor is composed of two structural subunits of GluN1 and two functional subunits of GluN2 or GluN3. The subunit composition determines the function of different NMDA receptor subtypes; among them, GluN2A and GluN2B are two major types expressed in the neocortex and hippocampus [27–29]. A variety of studies have proved that the GluN2A-containing NMDA receptor is gradually becoming the dominant subtype during development which is mainly caused by the increase of GluN2A [30–32]. And as a result, the threshold for LTP is raised and this is the prerequisite for synaptic pruning and fine tuning of the neuronal circuit [33]. This shift of subtype dominancy is regarded as a symbol of synaptic maturation and is believed to be regulated by sensory stimulus [34, 35], but its molecular mechanism is not fully uncovered yet. Currently, it is known that the phosphorylation of the GluN2B PDZ domain by CK2 may participate in the synaptic clearance of GluN2B [36] and the repressor element 1-silencing transcription factor (REST) could also decrease the GluN2B expression though epigenetic modification of the gene of GluN2B [37]. And for GluN2A, it is believed that the extracellular matrix is very important for its developmental expression, especially the Reelin signaling pathway. Inhibition of the Reelin signaling pathway tempers the GluN2A/GluN2B switch [38, 39] and also impairs the learning and memory [40]. Consistently, we found in our study that SIL1 was not directly involved in the transient dynamic trafficking of GluN2A in mature neurons but it could regulate the developmental increase of GluN2A through the Reelin signaling pathway. After SIL1 silencing, the expression of both Reelin receptors decreased and the phosphorylation of their downstream hub molecule Dab1 was inhibited. According to the previous studies, this will in turn obstruct the expression of GluN2A and affect the learning ability of the mice [23, 38, 40, 41] as we have found in our current study. More importantly, we found that the mimic phosphorylation peptide of Dab1 was able to partially rescue the expression of GluN2A and the spatial learning ability of the mice. Therefore, our data suggested that the functional deficiency of SIL1 could impair spatial learning of mice because of the abnormal expression of GluN2A induced by malfunction of Reelin signaling.

Finally, our data combined with other studies indicated that compared to Bip, SIL1 potentially has more particular functions. SIL1 showed cell-type selectivity in the nervous system. The Amyotrophic lateral sclerosis (ALS)-vulnerable FF motor neurons have much less SIL1 compared to the disease-resistant S motor neurons, and this difference is the main reason for the different resistance of the two types of neurons against ALS [20]. And SIL1 also showed substrate selectivity. For example, the normal secretion of antibody in SIL1 knockout mice indicated that SIL1 is not involved in or it is not essential for the antibody assembly in B cells [42] which is a major function of Bip [43, 44]. But SIL1 is required for the insulin secretion from mouse pancreatic beta cells [45]. Actually, there is another nucleotide exchange factor of Bip which is GRP170 [46]. It has been proved that GRP170 interacts with unassembled antibody [47, 48]. Therefore, the two nucleotide exchange factors SIL1 and GRP170 may function separately. In our study, we discovered a novel substrate selection that SIL1 regulated the expression of GluN2A but not its homolog GluN2B. This functional selectivity is consistent with the finding that unlike Bip, the SIL1-mutated mice are viable [17]. It will be interesting to investigate the mechanism of the SIL1 selectivity and its physiological significance in the development of the central nervous system.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors have no conflict of interests to declare.

Authors’ Contributions

Xu Shilian performed the experiments, contributed reagents and materials, and wrote the paper. Zhu Jia performed the experiments and contributed reagents and materials. Mi Kai performed the experiments and analyzed the data. Shen Yan performed the experiments and analyzed the data. Zhang Xiaomin conceived and designed the experiments, performed the experiments, analyzed and interpreted the data, contributed reagents and materials, and wrote the paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant Nos. 31460258, 31360245, and 81501043), Applied Basic Research Foundation of Yunnan Province of China (Grant No. 2018FE001(-134)), One Hundred Talented Young Scholars Programme of Kunming Medical University (Grant No. 60117190404), and Foundation of Science and Technology Innovative Team Building of Kunming Medical University (Grant No. CXTD201807). The National Natural Science Foundation of China and Applied Basic Research Foundation of Yunnan Province of China provided the reagents, materials, and labor cost. The One Hundred Talented Young Scholars Programme of Kunming Medical University provided the traveling cost for the communication among authors.

Supplementary Materials

Mice were subjected to bilateral intraventricular injection of lentivirus-encoding shRNA#2 at 3 weeks old followed or not followed by daily intraperitoneal injection of pTyr-220. Mice that received scramble shRNA injection were used as the control (con). Then, the whole cortex lysate was obtained at 5 weeks and GluN2A was measured by Western blot. (B) Statistical analysis of (A). The expression of GluN2A was normalized to GAPDH first and then compared with the control (con) in which group neurons were transfected with scramble shRNA. After SIL1 silencing, GluN2A expression was reduced significantly (, , and ) but rescued after applying pTyr-220 (, compared to shRNA#2 treatment neurons and compared to the control). All data were shown as . ##. and for all treatments. (Supplementary Materials)

References

K. T. Chung, Y. Shen, and L. M. Hendershot, “BAP, a mammalian BiP-associated protein, is a nucleotide exchange factor that regulates the ATPase activity of BiP,” Journal of Biological Chemistry, vol. 277, no. 49, pp. 47557–47563, 2002.

J. H. Otero, B. Lizák, and L. M. Hendershot, “Life and death of a BiP substrate,” Seminars in Cell & Developmental Biology, vol. 21, no. 5, pp. 472–478, 2010.

J. R. Tyson and C. J. Stirling, “LHS1 and SIL1 provide a lumenal function that is essential for protein translocation into the endoplasmic reticulum,” The EMBO Journal, vol. 19, no. 23, pp. 6440–6452, 2000.

A. Bertolotti, Y. Zhang, L. M. Hendershot, H. P. Harding, and D. Ron, “Dynamic interaction of BiP and ER stress transducers in the unfolded-protein response,” Nature Cell Biology, vol. 2, no. 6, pp. 326–332, 2000.

B. Kleizen and I. Braakman, “Protein folding and quality control in the endoplasmic reticulum,” Current Opinion in Cell Biology, vol. 16, no. 4, pp. 343–349, 2004.

I. G. Haas, “BiP (GRP78), an essential hsp70 resident protein in the endoplasmic reticulum,” Experientia, vol. 50, no. 11-12, pp. 1012–1020, 1994.

J. Dudek, J. Benedix, S. Cappel et al., “Functions and pathologies of BiP and its interaction partners,” Cellular and Molecular Life Sciences, vol. 66, no. 9, pp. 1556–1569, 2009.

A. K. Anttonen, I. Mahjneh, R. H. Hämäläinen et al., “The gene disrupted in Marinesco-Sjögren syndrome encodes SIL1, an HSPA5 cochaperone,” Nature Genetics, vol. 37, no. 12, pp. 1309–1311, 2005.

J. Senderek, M. Krieger, C. Stendel et al., “Mutations in SIL1 cause Marinesco-Sjögren syndrome, a cerebellar ataxia with cataract and myopathy,” Nature Genetics, vol. 37, no. 12, pp. 1312–1314, 2005.

R. Herva, L. Wendt, G. Wendt, A. L. Saukkonen, J. Leisti, and V. Dubowitz, “A syndrome with juvenile cataract, cerebellar atrophy, mental retardation and myopathy,” Neuropediatrics, vol. 18, no. 3, pp. 164–169, 1987.

M. Krieger, A. Roos, C. Stendel et al., “SIL1 mutations and clinical spectrum in patients with Marinesco-Sjögren syndrome,” Brain, vol. 136, no. 12, pp. 3634–3644, 2013.

F. Ezgu, P. Krejci, S. Li et al., “Phenotype–genotype correlations in patients with Marinesco–Sjögren syndrome,” Clinical Genetics, vol. 86, no. 1, pp. 74–84, 2014.

M. Goto, M. Okada, H. Komaki et al., “A nationwide survey on Marinesco-Sjögren syndrome in Japan,” Orphanet Journal of Rare Diseases, vol. 9, no. 1, p. 58, 2014.

M. Horvers, A. K. Anttonen, A. E. Lehesjoki et al., “Marinesco-Sjögren syndrome due to SIL1 mutations with a comment on the clinical phenotype,” European Journal of Paediatric Neurology, vol. 17, no. 2, pp. 199–203, 2013.

A. K. Anttonen, E. Siintola, L. Tranebjaerg et al., “Novel SIL1 mutations and exclusion of functional candidate genes in Marinesco–Sjögren syndrome,” European Journal of Human Genetics, vol. 16, no. 8, pp. 961–969, 2008.

J. Howes, Y. Shimizu, M. J. Feige, and L. M. Hendershot, “C-terminal mutations destabilize SIL1/BAP and can cause Marinesco-Sjögren syndrome,” Journal of Biological Chemistry, vol. 287, no. 11, pp. 8552–8560, 2012.

L. Zhao, C. Longo-Guess, B. S. Harris, J. W. Lee, and S. L. Ackerman, “Protein accumulation and neurodegeneration in the woozy mutant mouse is caused by disruption of SIL1, a cochaperone of BiP,” Nature Genetics, vol. 37, no. 9, pp. 974–979, 2005.

X. M. Zhang, X. Y. Yan, B. Zhang et al., “Activity-induced synaptic delivery of the GluN2A-containing NMDA receptor is dependent on endoplasmic reticulum chaperone Bip and involved in fear memory,” Cell Research, vol. 25, no. 7, pp. 818–836, 2015.

Y. Inaguma, N. Hamada, H. Tabata et al., “SIL1, a causative cochaperone gene of Marinesco-Sjogren syndrome, plays an essential role in establishing the architecture of the developing cerebral cortex,” EMBO Molecular Medicine, vol. 6, no. 3, pp. 414–429, 2014.

A. F. de L'Etang, N. Maharjan, M. C. Braña et al., “Marinesco-Sjögren syndrome protein SIL1 regulates motor neuron subtype-selective ER stress in ALS,” Nature Neuroscience, vol. 18, no. 2, pp. 227–238, 2015.

W. Rupik, K. Jasik, J. Bembenek, and W. Widłak, “The expression patterns of heat shock genes and proteins and their role during vertebrate’s development,” Comparative Biochemistry and Physiology Part A: Molecular & Integrative Physiology, vol. 159, no. 4, pp. 349–366, 2011.

Y. Chen, M. S. Durakoglugil, X. Xian, and J. Herz, “ApoE4 reduces glutamate receptor function and synaptic plasticity by selectively impairing ApoE receptor recycling,” Proceedings of the National Academy of Sciences of the United States of America, vol. 107, no. 26, pp. 12011–12016, 2010.

L. Groc, D. Choquet, F. A. Stephenson, D. Verrier, O. J. Manzoni, and P. Chavis, “NMDA receptor surface trafficking and synaptic subunit composition are developmentally regulated by the extracellular matrix protein Reelin,” The Journal of Neuroscience, vol. 27, no. 38, pp. 10165–10175, 2007.

E. Förster, “Reelin, neuronal polarity and process orientation of cortical neurons,” Neuroscience, vol. 269, pp. 102–111, 2014.

H. H. Bock and P. May, “Canonical and non-canonical Reelin signaling,” Frontiers in Cellular Neuroscience, vol. 10, p. 166, 2016.

Y. Huang, S. Magdaleno, R. Hopkins, C. Slaughter, T. Curran, and L. Keshvara, “Tyrosine phosphorylated disabled 1 recruits Crk family adapter proteins,” Biochemical and Biophysical Research Communications, vol. 318, no. 1, pp. 204–212, 2004.

S. Cull-Candy, S. Brickley, and M. Farrant, “NMDA receptor subunits: diversity, development and disease,” Current Opinion in Neurobiology, vol. 11, no. 3, pp. 327–335, 2001.

P. Paoletti and J. Neyton, “NMDA receptor subunits: function and pharmacology,” Current Opinion in Pharmacology, vol. 7, no. 1, pp. 39–47, 2007.

P. Paoletti, C. Bellone, and Q. Zhou, “NMDA receptor subunit diversity: impact on receptor properties, synaptic plasticity and disease,” Nature Reviews Neuroscience, vol. 14, no. 6, pp. 383–400, 2013.

E. B. Roberts and A. S. Ramoa, “Enhanced NR2A subunit expression and decreased NMDA receptor decay time at the onset of ocular dominance plasticity in the ferret,” Journal of Neurophysiology, vol. 81, no. 5, pp. 2587–2591, 1999.

A. C. Flint, U. S. Maisch, J. H. Weishaupt, A. R. Kriegstein, and H. Monyer, “NR2A subunit expression shortens NMDA receptor synaptic currents in developing neocortex,” The Journal of Neuroscience, vol. 17, no. 7, pp. 2469–2476, 1997.

M. Sheng, J. Cummings, L. A. Roldan, Y. N. Jan, and L. Y. Jan, “Changing subunit composition of heteromeric NMDA receptors during development of rat cortex,” Nature, vol. 368, no. 6467, pp. 144–147, 1994.

M. Yasuda, E. M. Johnson-Venkatesh, H. Zhang, J. M. Parent, M. A. Sutton, and H. Umemori, “Multiple forms of activity-dependent competition refine hippocampal circuits in vivo,” Neuron, vol. 70, no. 6, pp. 1128–1142, 2011.

C. Bellone and R. A. Nicoll, “Rapid bidirectional switching of synaptic NMDA receptors,” Neuron, vol. 55, no. 5, pp. 779–785, 2007.

E. M. Quinlan, B. D. Philpot, R. L. Huganir, and M. F. Bear, “Rapid, experience-dependent expression of synaptic NMDA receptors in visual cortex in vivo,” Nature Neuroscience, vol. 2, no. 4, pp. 352–357, 1999.

A. Sanz-Clemente, J. A. Matta, J. T. R. Isaac, and K. W. Roche, “Casein kinase 2 regulates the NR2 subunit composition of synaptic NMDA receptors,” Neuron, vol. 67, no. 6, pp. 984–996, 2010.

A. Rodenas-Ruano, A. E. Chávez, M. J. Cossio, P. E. Castillo, and R. S. Zukin, “REST-dependent epigenetic remodeling promotes the developmental switch in synaptic NMDA receptors,” Nature Neuroscience, vol. 15, no. 10, pp. 1382–1390, 2012.

M. Sinagra, “Reelin, very-low-density lipoprotein receptor, and apolipoprotein E receptor 2 control somatic NMDA receptor composition during hippocampal maturation in vitro,” The Journal of Neuroscience, vol. 25, no. 26, pp. 6127–6136, 2005.

Y. Chen, U. Beffert, M. Ertunc et al., “Reelin modulates NMDA receptor activity in cortical neurons,” The Journal of Neuroscience, vol. 25, no. 36, pp. 8209–8216, 2005.

U. Beffert, E. J. Weeber, A. Durudas et al., “Modulation of synaptic plasticity and memory by Reelin involves differential splicing of the lipoprotein receptor Apoer2,” Neuron, vol. 47, no. 4, pp. 567–579, 2005.

S. Qiu and E. J. Weeber, “Reelin signaling facilitates maturation of CA1 glutamatergic synapses,” Journal of Neurophysiology, vol. 97, no. 3, pp. 2312–2321, 2007.

V. P. Ichhaporia, T. Sanford, J. Howes, T. N. Marion, and L. M. Hendershot, “Sil1, a nucleotide exchange factor for BiP, is not required for antibody assembly or secretion,” Molecular Biology of the Cell, vol. 26, no. 3, pp. 420–429, 2015.

M. R. Knittler and I. G. Haas, “Interaction of BiP with newly synthesized immunoglobulin light chain molecules: cycles of sequential binding and release,” The EMBO Journal, vol. 11, no. 4, pp. 1573–1581, 1992.

M. J. Feige, S. Groscurth, M. Marcinowski et al., “An unfolded CH1 domain controls the assembly and secretion of IgG antibodies,” Molecular Cell, vol. 34, no. 5, pp. 569–579, 2009.

A. A. Ittner, J. Bertz, T. Y. B. Chan, J. Van Eersel, P. Polly, and L. M. Ittner, “The nucleotide exchange factor SIL1 is required for glucose-stimulated insulin secretion from mouse pancreatic beta cells in vivo,” Diabetologia, vol. 57, no. 7, pp. 1410–1419, 2014.

A. Weitzmann, J. Volkmer, and R. Zimmermann, “The nucleotide exchange factor activity of Grp170 may explain the non-lethal phenotype of loss of Sil1 function in man and mouse,” FEBS Letters, vol. 580, no. 22, pp. 5237–5240, 2006.

H. Wang, A. M. Pezeshki, X. Yu, C. Guo, J. R. Subjeck, and X.-Y. Wang, “The endoplasmic reticulum chaperone GRP170: from immunobiology to cancer therapeutics,” Frontiers in Oncology, vol. 4, p. 377, 2015.

H. Y. Lin, P. Masso-Welch, Y. P. di, J. W. Cai, J. W. Shen, and J. R. Subjeck, “The 170-kDa glucose-regulated stress protein is an endoplasmic reticulum protein that binds immunoglobulin,” Molecular Biology of the Cell, vol. 4, no. 11, pp. 1109–1119, 1993.