A Proteomic Network Approach across the Kidney Stone Disease Reveals Endoplasmic Reticulum Stress and Crystal-Cell Interaction in the Kidney

Crystal-cell interactions are a vital step toward kidney stone formation. However, its mechanisms remained unclear. Here, a protein-protein interaction (PPI) network analysis of a kidney stone revealed that the proteins were enriched in a posttranslational protein modification process in the endoplasmic reticulum (ER). The in vitro study showed that the markers of ER stress, including Bip and CHOP, were upregulated, PERK and ATF6 were activated, and XBP-1 mRNA was spliced. An ER stress-specific protein, caspase-12, was activated in the apoptotic cells induced by calcium oxalate monohydrate (COM) crystals. The treatment with tunicamycin, an ER stress inducer, promoted the crystal-cell adhesion assayed by atomic absorption, reduced cell viability assayed by MTT, and downregulated the expression of proteins involved in the crystal formations. The treatment with salubrinal, an ER stress inhibitor, reversed the above effects for both tunicamycin and COM crystals. The aforementioned main observations were supported by in vivo study. These data demonstrated that ER stress was an essentially biological process of crystal-cell interactions. Our findings suggest that blocking ER stress may become a potential approach to preventing a kidney stone.

1.Introduction
The incidence of kidney stone disease is about 2–5% of the population in Asia and 8–15% in Europe and North Amer- ica [1]. Approximately 80% of kidney stone types are cal- cium oxalate stone. The recurrence rate of a kidney stone may be higher than 50% after five years [2]. The significant challenges for clinicians are to prevent recurrence of kidney stone patients. However, kidney stone formation is a com-
plex response of cells to the exposure to crystals; the precise mechanism causing crystal-cell interactions is still unclear.
Both a network analysis and a gene ontology (GO) anal- ysis are rapidly becoming powerful tools in complex disease studies [3, 4]. It provides a visual framework and protein enrichment for specific functional categories [5]. A large number of proteins are reported to be involved in the pro- cess of kidney stone formation. Some macromolecules, such as osteopontin (OPN) [6], matrix Gla protein (MGP) [7],
bikunin [8], and Tamm-Horsfall proteins [9], have been identified in both the urine and kidney stone matrix, and their gene variants have also been reported to affect the risk of kidney stone disease [10–12]. Recently, Wright et al. iden- tified more than 1000 proteins by a urinary proteome analy- sis [13]. These proteins may interact with each other to play a vital role in modulating crystal nucleation, growth, aggre- gation, and adhesion to renal epithelial cells. It also provides us with a significant amount of unprocessed information for understanding crystal-cell interactions in the process of kid- ney stone formation.

In the present study, we collected human protein candi- dates associated with a kidney stone to perform a bioinfor- matics analysis and found that the candidate proteins were significantly enriched in a posttranslational protein modifi- cation process in the ER. The in vitro studies showed that markers associated with ER stress, Bip/GRP78 and CHOP, were upregulated, suggesting that ER stress was directly involved in crystal-cell interactions. Caspase-12, an ER stress-specific caspase protein, was activated. Tunicamycin, an ER stress inducer, increased crystal-cell adhesion, reduced cell viability, and downregulated the expression of proteins associated with kidney stone formation. Salubrinal, an ER stress inhibitor, can reverse the above effects of both tunica- mycin and COM crystals. The aforementioned main obser- vations were supported by our in vivo study. These findings identified an essential mechanism of crystal-cell interaction in kidneys and provided a potential strategy for the preven- tion and treatment of a kidney stone, by targeting a cellular ER stress pathway.

2.Materials and Methods
Proteins that are associated with kidney stone formation were obtained from a proteome study and a search in PolySearch using the keywords “kidney stone” or “nephrolithiasis”. The proteome study performed a label- free nano-ultraperformance liquid chromatography between the 57 stone formers and 57 non-stone-forming controls [13]. The protein names were converted into one unified form using “DAVID Bioinformatics Resources”.The candi- date proteins were integrated into an “InWeb” PPI database [14]. A permutation test was performed to evaluate whether candidate proteins are significantly linked via PPI networks, rather than by chance. Highly connected proteins were obtained using a GeneNet Toolbox [15]. A backbone network was extracted from the top 10% hub proteins according to rank by the “node degree,” “betweenness centrality,” and “edge betweenness.” The PPI network was visualized using the Cytoscape software [16]. A GO analysis was performed to find protein enrichment using “DAVID Bioinformatics Resources” [17]. GO terms including biological processes (BP), cellular components (CC), and molecular functions (MF) were exam- ined.

A flowchart was developed (Figure 1(a)).The HK2 (human kidney 2) cells (American Type Culture Collection) were exposed for three, six, 12, 24,and 48 hours to 100 μg/ml calcium oxalate monohydrate (COM) crystals, 1 μg/ml tunicamycin, and 1 μM salubrinal in DMEM, respectively. Tunicamycin, an inhibitor of N- linked glycosylation, can cause cell cycle arrest in a G1 phase to induce an unfolded protein response [18]. Salu- brinal is an inhibitor of ER stress, selectively inhibiting eIF2α dephosphorylation.Male Sprague-Dawley (SD) rats were divided into two groups, a treatment group (n = 25) and a control group (n = 5). The treatment group was administered 0.75% ethylene glycol (EG) in drinking water for up to eight weeks, and the control group was given plain water. Rats were dissected at zero-, two-, four-, six-, and eight-week intervals. The experiment was reviewed and approved by the Shenyang Medical College Animal Committee.The cells were homogenized using a homogenizer (Vibra Cell, Sonics, USA), and total RNA was extracted with an RNAiso kit (Takara, Japan). Total RNA was treated with DNase I for 30 minutes at 37°C, to avoid genomic DNA contamination.A semiquantitative PCR assay was performed using previously described methods [7]. The first-strand cDNA was synthesized from 2 μg of total RNA using an EasyScript First-Strand cDNA Synthesis SuperMix kit (TransGen Biotech, Beijing).

The glyceraldehyde-3- phosphate dehydrogenase gene (GAPDH) was used as an endogenous control. The primers were as follows: human XBP-1 (X-box binding protein 1) sense 5′-TTACGAGAGAAAACTCATGGCC-3′ and human XBP-1 antisense 5′-GGGTCCAAGTTGTCCAGAATGC-3′, rat XBP-1 sense 5′-TTACGAGAGAAAACTCATGGGC-3′ and rat XBP-1antisense 5′-GGGTCCAACTTGTCCAGAATGC-3′, human GAPDH sense 5′-GTCTCCTCTGACTTCAACAGCG-3′ andhuman GAPDH antisense 5′-ACCACCCTGTTGCTGTAGC CAA-3′, and rat GAPDH sense 5′-TGCTGAGTATGTCGTG GAGTCTA-3′, rat GAPDH antisense 5′-AGTGGGAGTTG CTGTTGAAATC-3′. The assay was performed in triplicate.Proteins were isolated from the culture cells by homogenization in RIPA buffer (Solarbio) with protease inhibitors. Primary antibodies—anti-caspase-12 (AB_2069200), anti-caspase-3 (AB_331439), anti- phospho-PERK (PKR-like endoplasmic reticulum kinase, AB_2095853), anti-phospho-P38 (AB_2139682), and anti- phospho-JNK (c-Jun N-terminal kinase, AB_331338) anti- bodies—were purchased from Cell Signaling (Beverly, MA). A mouse monoclonal anti-Bip/GRP78 (immunoglobulin heavy-chain binding protein, Grp78, AB_398292) antibody was purchased from BD Biosciences (Bedford, MA). Rabbit anti-CHOP (C/EBP homologous protein, AB_1122157) and mouse anti-OPN (AB_2194995) were purchased from Santa Cruz Biotechnology (CA, USA), rabbit polyclonal anti-MGP (cat no. TA339184) was purchased from OriGene, and mouse monoclonal anti-ATF6 (activating transcription fac- tor, AB_2058774) was purchased from Novus. Beta-actinantibodies were used to confirm equal protein loading among samples. All assays were repeated in triplicate.Apoptotic cells were detected using a terminal deoxynucleotidyl transferase-mediated X-dUTP nick end labeling (TUNEL) method, using the in situ apopto- sis detection kit (Takara Biomedicals), following the manufac- turer’s instructions. The cells were observed using a confocal laser-scanning microscope (Leica, Heidelberg, Germany).

Evident nuclear staining represented the apoptotic cells.Cell viability was mea- sured using the MTT assay. Briefly, 1× 106 cells were plated into a 96-well plate in 100 μl volume, and then, an appropri- ate concentration of COM crystals, tunicamycin, and salu- brinal was added to HK2 cells for two days of exposure, respectively. A labeling dye was added four hours before the endpoint. A stop solution was added to each well to stop the reaction and solubilize the cells. The absorbance was read at 490 nm.Atomic absorption was per- formed to assay the effect of COM crystal adhesion to cells. COM crystals were treated with ultrasound for 15 minutes and suspended in DMEM at a final concentration of 2 mg/ml. The cells were exposed to 100 μg/ml COM crystals for five minutes. The nonadherent crystals were removed by rapidly washing the plate three times with the PBS solution. The COM crystals that had adhered to cells were lysed with 5 ml of 6 M HCl. A quantitative anal- ysis of the COM crystals was conducted by measuring the calcium concentration of the supernatants using an atomic absorption analysis.Candidate proteins were collected from a PolySearch search with Z score ≥0.1. A permutation test was performed to evaluate network reliability using the GeneNet Toolbox software. A P value of <0.05 and a false discovery rate (FDR) of <10 were considered statistically sig-nificant for GO analysis. Comparisons were analyzed using a two-way analysis of variance (ANOVA) with Dunnett’s multiple comparison test. A P value of ≤0.05 was considered to be significant.

3.Results
Construction of the PPI Network of a Kidney Stone. Seven hundred seventy-five proteins were obtained from the prote- ome study, and 112 human proteins were found through a search conducted in PolySearch. A total of 839 candidate proteins were obtained after a name conversion and remov- ing duplicates, and they were then integrated into the InWeb database to construct a kidney stone PPI network. The net- work contained 340 nodes and 740 edges, which developed a permutation test at P < 0.001 (n = 1000 permutations) (Figure 1(b)). Thirty-one hub proteins were extracted to obtain a backbone PPI network of a kidney stone, which pre- sented an underlying biological network of kidney stone for- mation (Figure 1(c)).Identification of Biological Processes of Crystal-Cell Interaction. A GO analysis was performed for the thirty- one hub proteins. The GO_BP analysis demonstrated that the proteins had enriched eight biological processes (Table 1). Five of the processes were found to be associated with ER, such as “protein folding” (P = 1.33E − 05), “proteinstabilization” (P = 8.13E − 05), “ATF6-mediated unfoldedprotein response” (P = 1.03E − 04), and “protein folding in endoplasmic reticulum” (P = 2.22E − 04). The GO_CC analysis showed that the proteins were enriched in the “ER chaperone complex” (P = 1.34E − 04) and “ER lumen” (P = 2.32E − 04). The GO_MF analysis found that the pro- teins had enriched the “unfolded protein binding” (P = 1.32E − 06). An analysis of the KEGG pathway demon- strated that the proteins were enriched during “protein pro- cessing in ER” (P = 4.30E − 04).

These results suggest that a critical biological event in the ER might occur during the pro-cess of kidney stone formation. To confirm whether COM crystals activated ER stress, HK2 cells were exposed to crys- tals for different times and a western blot analysis was per- formed. The Bip/GRP78, an ER chaperone protein which interacts with initial molecules of unfolded protein responses [19], was meaningfully increased from three to 48 hours, and the peak value presented at six hours (Figures 2(a) and 2(c)). The CHOP, a downstream protein that targets ER stress [19], was significantly upregulated from 12 to 24 hours (Figures 2(a) and 2(d)). These results suggest that ER stress was initiated by exposure to COM crystals.To clarify which signaling pathway is involved in ER stress induced by COM crystals, we performed a western blot and semiquantitative RT-PCR. Three major ER stress response pathways were investigated, including PERK, inositol-requiring enzyme 1 (IRE1), and ATF6 pathways. PERK is a kinase that phosphorylates eIF2a, leading to a rapid reduction of translational initiation and repression of global protein synthesis [20]. XBP-1 is a transcription factor, and its mRNA is spliced by activated IRE1 resulting in the regulation of expression of some downstream genes [21].

ATF6 is an ER stress-regulated transmembrane transcription factor that is cleaved under ER stress, and the cleaved fragment subsequently moves to the nucleus and binds to DNA as a transcription factor, which causes the transcription of ER chaperones. The PERK protein was strongly activated as early as three hours after COM crystal exposure, and it was gradually decreased from 12 to 48 hours (Figures 2(a) and 2(e)). The full ATF6 protein presented a downward trend from three to 48 hours, whereas the cleaved ATF6 pro- tein increased from three hours and peaked at 12 hours (Figures 2(a) and 2(f)). The level of total XBP-1 protein expression was decreased, and it was then upregulated to 24 hours (Figures 2(a) and 2(g)). The semiquantitative PCR analysis showed that the spliced XBP-1 mRNA was upregu- lated from three hours to 12 hours (Figures 2(b) and 2(h)). These results demonstrate that COM crystals induce ER stress through at least three major pathways, including PERK, IRE1, and ATF6.ER Stress-Specific Caspase Was Activated in Apoptotic Cells Induced by COM Crystals. To understand whether the ER-related proteins were involved in crystal-induced apoptosis, HK2 cells were exposed to COM crystals for 24 hours, and apoptosis was detected by a TUNEL and Hoechst stain assay (Figure 3(a)). The number of apoptotic cells represented by the strong fluorescent signals in the COM-treated group was significantly more than what was found in the control group (Figure 3(b)).

Activated caspase- 12, an ER stress-specific caspase that mediates ER stress- induced apoptosis, increased markedly after three hours (Figures 3(c) and 3(d)). The cleaved caspase-3, a converging point of an apoptotic pathway, was increased at three and 12 hours, respectively (Figures 3(c) and 3(e)). The JNK and P38, stress-related proteins, significantly increased compared with the control group (Figures 3(c), 3(f), and 3(g)). These results suggest that ER stress was associated with COM crystal- induced apoptosis.3.5.The Effect of ER Stress on Crystal-Cell Interaction. To understand the effect of ER stress on crystal-cell adhesion, we analyzed the calcium concentration of crystals that adhered to HK2 cells treated with tunicamycin using an atomic absorption method. Crystal-cell adhesion is a crucialstep for kidney stone formation. Quantitative crystal analysis showed that the calcium content remaining on the cell surface of the tunicamycin group significantly increased, compared with that of the control group cells (Figure 4(a)). The treatment with salubrinal markedly reversed the effects of ER stress on crystal-cell adhesion (Figure 4(a)).

To investigate the effect of ER stress on the secretory pro- teins involved in crystal formation, we treated the HK2 cells with tunicamycin for three, six, 12, 24, and 48 hours and detected the expressions of OPN and MGP. OPN is the main stone matrix component, and MGP is a determined protein of extracellular matrix calcification. OPN expression pickedup slowly, and MGP expression increased temporarily during the first three hours, and then, it gradually decreased after exposure to tunicamycin (Figure 4(b)). The pattern of expression in either OPN or MGP was similar to the treat- ment with COM crystals and the treatment with tunicamycin (Figure 4(b)).To understand the effect of ER stress on cell viability, the MTT assay was performed. Figure 4(c) shows that the treat- ment with both COM crystals and tunicamycin reduced absorbance compared to the control cells (P < 0.05). Simi- larly, treatment with salubrinal meaningfully reversed the effects of COM crystals on cell viability (Figure 4(c)). Also, salubrinal significantly increased the number of viable cells, compared to the group that was only exposed to COM crys- tals, in the 24th hour (Figure 4(d)).

These results suggest that ER stress induced by both tunicamycin and calcium oxalate may promote crystal formation by increasing crystal-cell adhesion, affecting the expression of proteins involved in crystal formation such as OPN and MGP and reducing cell viability. The inhibition of ER stress can, therefore, reduce crystal formation.ER Stress Occurs in Rats with Crystal Formation. To confirm the above findings in vivo, SD rats were treated with 0.75% EG for two, four, six, and eight weeks. The formation of crystals was observed from the second week (Figure 5(a)). Similar to that found in vitro, the Bip/GRP78 and CHOP were increased from the second to the 8th week. The PERK and ATF6 were strongly activated from the second week,and caspase-12 and caspase-3 were cleaved from the second to the eighth week (Figure 5(b)). XBP-1 was spliced at the second to the 8th week (Figure 5(c)). These results demon- strate that ER stress was also induced in rats with crystal for- mation, through all of the three central ER stress pathways.

4.Discussion
In the present study, we performed a kidney stone network analysis and a GO analysis. The study found that the candi- date proteins were enriched in a posttranslational protein modification process in the ER. Our in vitro and in vivo study found that the markers of ER stress, Bip/GRP78 and CHOP, increased the expression in response to COM exposure, which is accompanied by the activation of all of the three major ER stress signaling pathways, including PERK, XBP- 1, and ATF6. An ER stress-specific caspase protein was acti- vated in apoptotic cells induced by COM crystals. We further investigated the effect of ER stress on crystal-cell interaction and found that ER stress led to reducing cell viability and affecting the expression of secretory proteins associated with crystal formation and increasing crystal-cell adhesion. These findings suggest that ER stress contributes to crystal forma- tion. Furthermore, an inhibitor of ER stress can reverse cell viability and reduce crystal-cell adhesion. These findings sug- gest that inhibitors to ER stress may become a promising strategy for the prevention and treatment of a kidney stone. Omics studies of a kidney stone have presented vast superiority in higher throughput and larger scale comparedwith traditional studies. Moreover, it may provide us with a panorama for interaction between the proteins and clues for finding a novel mechanism of the disease.

However, net- work analysis may also provide us with false information due to candidate proteins linked via PPIs by chance. In the present study, we performed a permutation test to eval- uate whether candidate proteins are significantly linked via PPIs, rather than discovering an outcome by chance. This approach provided us with a network that was as realistic as possible. We also extracted the backbone network from the primary kidney stone network to obtain more essential proteins associated with a kidney stone. Hence, these strate- gies have significantly improved the reliability of the subse- quent GO analysis.ER stress has been found in various diseases such as diabetes, inflammation, and neurodegenerative disorders. However, concerning a kidney stone, there is little evidence. Hyperoxaluria is a major metabolic promoter of kidney stone formation. A previous study found that the oxalatetoxicity increased GRP78 and CHOP expression at an mRNA level in NRK49F cells [22]. High oxalate concentra- tions are toxic to renal tubular cells, resulting in necrotic cell death and increasing the risk of stone formation [23]. The current study clarified, for the first time, that ER stress is directly involved in crystal-cell interaction, and the apopto- tic pathway in renal tubular cells induced by COM crystals. These findings suggest that ER stress occurs in the full pro- cess of kidney stone formation.It is believed that the apoptosis of renal tubular epithelial cells might facilitate the formation of the nucleus of the crys- tal due to the attachment of calcium ions to the phosphati- dylserine, which flips out of the cell membrane during the early stages of apoptosis [24].

However, the precise mecha- nism regarding the apoptosis of renal tubular epithelial cells induced by COM remains unknown. Our results demon- strated that ER stress induced by COM crystals was associ- ated with caspase-12 activation. Following its activation at the ER, caspase-12 directly cleaves procaspase-9, leading tocaspase-9-dependent activation of caspase-3. Also, ER stress-induced apoptosis is regulated by both the CHOP and the ATF6 pathway. CHOP, a downstream of the PERK pathway, is a transcription factor that inhibits the expression of Bcl-2 and promotes apoptosis [25]. CHOP has been found to induce the expression of numerous proapoptotic factors, which promotes protein synthesis and oxidative stress [26]. To confirm the effect of COM crystals on stress-related cell signaling, we investigated the activation of JNK. The activation of JNK significantly increased, com- pared with the control cells. Activated JNK induces the expression of inflammatory genes by phosphorylating the transcription factor activator protein to promote apoptosis [19, 27]. Constant ER stress can lead to further cell apopto- sis, inflammation, and oxidant production.Increasing evidence shows that some proteins may play critical modulating roles in the development of a kidney stone [7, 28].

OPN is a natural inhibitor of abnormal calcifi- cation, and it plays a vital inhibitory role in the phases of stone formation, including crystallization, crystal retention, crystal congregation, and stone formation in vitro or in vivo [28]. MGP is a molecular determinant regulating calcifica- tion of the extracellular matrix. Crystals were observed to form in the renal tubules with a lack of MGP expression, whereas crystals barely formed in the renal tubules with reg- ular MGP expression [7]. In the current study, we analyzed the expressions of these two proteins in ER stress cells induced by tunicamycin. Our results showed that ER stress induced by tunicamycin finally affected these proteins’ expression and promoted crystal formation. We believe that both TM and COM could initially upregulate the expression of MGP and OPN to avoid crystal formation, and continuous exposure to TM or COM may induce feedback that increases the ER stress thus inducing apoptosis and reducing the capa- bility of the cell to manage COM.

5.Conclusions
The crystal-cell interactions are a crucial step toward kidney stone formation. This study seeks to explore what biological processes contribute to crystal-cell interaction. We find that ER stress may increase the risk of crystal formation by being involved in cell apoptosis, attenuating the protein expression associated with kidney stone formation, reducing cell viabil- ity, and increasing Salubrinal crystal-cell adhesion. Our findings pro- vide a potential strategy for the prevention and treatment of kidney stones by targeting cellular ER stress pathways.