A Proteomic Network Approach across the Kidney Stone Disease Reveals Endoplasmic Reticulum Stress and Crystal-Cell Interaction in the Kidney.
The incidence of kidney stone disease is about 2-5% of the population in Asia and 8-15% in Europe and North America . Approximately 80% of kidney stone types are calcium oxalate stone. The recurrence rate of a kidney stone may be higher than 50% after five years . 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) analysis are rapidly becoming powerful tools in complex disease studies [3, 4]. It provides a visual framework and protein enrichment for specific functional categories . A large number of proteins are reported to be involved in the process of kidney stone formation. Some macromolecules, such as osteopontin (OPN) , matrix Gla protein (MGP) , bikunin , and Tamm-Horsfall proteins , 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. identified more than 1000 proteins by a urinary proteome analysis . These proteins may interact with each other to play a vital role in modulating crystal nucleation, growth, aggregation, 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 kidney stone formation.
In the present study, we collected human protein candidates associated with a kidney stone to perform a bioinformatics analysis and found that the candidate proteins were significantly enriched in a posttranslational protein modification 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 tunicamycin and COM crystals. The aforementioned main observations 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 prevention and treatment of a kidney stone, by targeting a cellular ER stress pathway.
2. Materials and Methods
2.1. Protein Dataset. 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 . The protein names were converted into one unified form using "DAVID Bioinformatics Resources".
2.2. PPI Network Construction and GO Analysis. The candidate proteins were integrated into an "InWeb" PPI database . 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 . 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 . A GO analysis was performed to find protein enrichment using "DAVID Bioinformatics Resources" . GO terms including biological processes (BP), cellular components (CC), and molecular functions (MF) were examined. A flowchart was developed (Figure 1(a)).
2.3. Cell Culture. The HK2 (human kidney 2) cells (American Type Culture Collection) were exposed for three, six, 12, 24, and 48 hours to 100 [micro]g/ml calcium oxalate monohydrate (COM) crystals, 1 [micro]g/ml tunicamycin, and 1 [micro]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 . Salubrinal is an inhibitor of ER stress, selectively inhibiting eIF2a dephosphorylation.
2.4. Animal Models. 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.
2.5. RNA Preparation. 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[degrees]C, to avoid genomic DNA contamination.
2.6. Semiquantitative PCR. A semiquantitative PCR assay was performed using previously described methods . The first-strand cDNA was synthesized from 2 [micro]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'-TTACGAGAG AAAACTCATGGCC-3' and human XBP-1 antisense 5'-GGGTCCAAGTTGTCCAGAATGC-3', rat XBP-1 sense 5'-TTACGAGAGAAAACTCATGGGC-3' and rat XBP-1 antisense 5'-GGGTCCAACTTGTCCAGAATGC-3', human GAPDH sense 5'-GTCTCCTCTGACTTCAACAGCG-3' and human 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.
2.7. Western Blot Analysis. 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), antiphospho-PERK (PKR-like endoplasmic reticulum kinase, AB_2095853), anti-phospho-P38 (AB_2139682), and anti-phospho-JNK (c-Jun N-terminal kinase, AB_331338) antibodies--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 factor, AB_2058774) was purchased from Novus. Beta-actin antibodies were used to confirm equal protein loading among samples. All assays were repeated in triplicate.
2.8. TUNEL Assay. Apoptotic cells were detected using a terminal deoxynucleotidyl transferase-mediated X-dUTP nick end labeling (TUNEL) method, using the in situ apoptosis detection kit (Takara Biomedicals), following the manufacturer's instructions. The cells were observed using a confocal laser-scanning microscope (Leica, Heidelberg, Germany). Evident nuclear staining represented the apoptotic cells.
2.9. Measurement of Cell Viability. Cell viability was measured using the MTT assay. Briefly, 1 x [10.sup.6] cells were plated into a 96-well plate in 100 [micro]l volume, and then, an appropriate concentration of COM crystals, tunicamycin, and salubrinal 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.
2.10. Atomic Absorption Assay. Atomic absorption was performed 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 [micro]/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 analysis of the COM crystals was conducted by measuring the calcium concentration of the supernatants using an atomic absorption analysis.
2.11. Statistical Analysis. Candidate proteins were collected from a PolySearch search with Z score [greater than or equal to] 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 significant for GO analysis. Comparisons were analyzed using a two-way analysis of variance (ANOVA) with Dunnett's multiple comparison test. A P value of [less than or equal to] 0.05 was considered to be significant.
3.1. Construction of the PPI Network of a Kidney Stone. Seven hundred seventy-five proteins were obtained from the proteome 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 removing duplicates, and they were then integrated into the InWeb database to construct a kidney stone PPI network. The network 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 presented an underlying biological network of kidney stone formation (Figure 1(c)).
3.2. 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), "protein stabilization" (P = 8.13E - 05), "ATF6-mediated unfolded protein 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 proteins had enriched the "unfolded protein binding" (P = 1.32E - 06). An analysis of the KEGG pathway demonstrated that the proteins were enriched during "protein processing in ER" (P = 4.30E - 04). These results suggest that a critical biological event in the ER might occur during the process of kidney stone formation.
3.3. COM Crystals Activated ER Stress through Three Response Signaling Pathways. To confirm whether COM crystals activated ER stress, HK2 cells were exposed to crystals for different times and a western blot analysis was performed. The Bip/GRP78, an ER chaperone protein which interacts with initial molecules of unfolded protein responses , 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 , 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 . XBP-1 is a transcription factor, and its mRNA is spliced by activated IRE1 resulting in the regulation of expression of some downstream genes . 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 protein 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 upregulated 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.
3.4. 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 crucial step 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 proteins 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 picked up 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 treatment 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 treatment with both COM crystals and tunicamycin reduced absorbance compared to the control cells (P < 0.05). Similarly, 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 crystals, 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.
3.6. 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 demonstrate that ER stress was also induced in rats with crystal formation, through all of the three central ER stress pathways.
In the present study, we performed a kidney stone network analysis and a GO analysis. The study found that the candidate 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 activated 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 formation. Furthermore, an inhibitor of ER stress can reverse cell viability and reduce crystal-cell adhesion. These findings suggest 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 compared with 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, network 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 evaluate 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 strategies have significantly improved the reliability of the subsequent 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 oxalate toxicity increased GRP78 and CHOP expression at an mRNA level in NRK49F cells . High oxalate concentrations are toxic to renal tubular cells, resulting in necrotic cell death and increasing the risk of stone formation . The current study clarified, for the first time, that ER stress is directly involved in crystal-cell interaction, and the apoptotic pathway in renal tubular cells induced by COM crystals. These findings suggest that ER stress occurs in the full process of kidney stone formation.
It is believed that the apoptosis of renal tubular epithelial cells might facilitate the formation of the nucleus of the crystal due to the attachment of calcium ions to the phosphatidylserine, which flips out of the cell membrane during the early stages of apoptosis . However, the precise mechanism regarding the apoptosis of renal tubular epithelial cells induced by COM remains unknown. Our results demonstrated that ER stress induced by COM crystals was associated with caspase-12 activation. Following its activation at the ER, caspase-12 directly cleaves procaspase-9, leading to caspase-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 . CHOP has been found to induce the expression of numerous proapoptotic factors, which promotes protein synthesis and oxidative stress . To confirm the effect of COM crystals on stress-related cell signaling, we investigated the activation of JNK. The activation of JNK significantly increased, compared 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 apoptosis, 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 calcification, 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 . MGP is a molecular determinant regulating calcification 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 regular MGP expression . 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 capability of the cell to manage COM.
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 viability, and increasing crystal-cell adhesion. Our findings provide a potential strategy for the prevention and treatment of kidney stones by targeting cellular ER stress pathways.
ER: Endoplasmic reticulum
PPI: Protein-protein interaction
MGP: Matrix Gla protein
HK2: Human kidney 2
COM: Calcium oxalate monohydrate
MTT: 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide
GO: Gene ontology
TUNEL: Terminal deoxynucleotidyl transferase-mediated X-dUTP nick end labeling
PCR: Polymerase chain reaction
Bip/GRP78: Immunoglobulin heavy-chain binding protein
GAPDH: Glyceraldehyde-3-phosphate dehydrogenase
PERK: PKR-like endoplasmic reticulum kinase
XBP-1: X-box binding protein 1
CHOP: C/EBP homologous protein
ATF6: Activating transcription factor
JNK: c-Jun N-terminal kinase
FDR: False discovery rate.
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
Baoyu Yang and Xiuli Lu contributed equally to this work.
This study was supported by the National Natural Science Foundation of China (no. 81570632), the Innovation Talent Project from the Education Department of Liaoning Province (nos. LR2016063 and LFW201701), and the Shenyang Science and Technology Bureau (nos. RC170272 and 18-007-0-02).
 C. Y. C. Pak, "Kidney stones," The Lancet, vol. 351, no. 9118, pp. 1797-1801, 1998.
 F. L. Coe, J. H. Parks, and J. R. Asplin, "The pathogenesis and treatment of kidney stones," New England Journal of Medicine, vol. 327, no. 16, pp. 1141-1152, 1992.
 V. Nair, C. V. Komorowsky, E. J. Weil et al., "A molecular morphometric approach to diabetic kidney disease can link structure to function and outcome," Kidney International, vol. 93, no. 2, pp. 439-449, 2018.
 N. Su, H. P. Choi, F. Wang et al., "Quantitative proteomic analysis of differentially expressed proteins and downstream signaling pathways in chronic bladder ischemia," Journal of Urology, vol. 195, no. 2, pp. 515-523, 2016.
 C. Y. Chen, A. Ho, H. Y. Huang, H. F. Juan, and H. C. Huang, "Dissecting the human protein-protein interaction network via phylogenetic decomposition," Scientific Reports, vol. 4, no. 1, 2014.
 K. Kohri, S. Nomura, Y. Kitamura et al., "Structure and expression of the mRNA encoding urinary stone protein (osteopontin)," Journal of Biological Chemistry, vol. 268, no. 20, pp. 15180-15184, 1993.
 X. Lu, B. Gao, T. Yasui et al., "Matrix Gla protein is involved in crystal formation in kidney of hyperoxaluric rats," Kidney and Blood Pressure Research, vol. 37, no. 1, pp. 15-23, 2013.
 F. Atmani and S. R. Khan, "Role of urinary bikunin in the inhibition of calcium oxalate crystallization," Journal of the American Society of Nephrology: JASN, vol. 10, Suppl 14, pp. S385-S388, 1999.
 L. Mo, H. Y. Huang, X. H. Zhu, E. Shapiro, D. L. Hasty, and X. R. Wu, "Tamm-Horsfall protein is a critical renal defense factor protecting against calcium oxalate crystal formation," Kidney International, vol. 66, no. 3, pp. 1159-1166, 2004.
 B. Gao, T. Yasui, Y. Itoh et al., "Association of osteopontin gene haplotypes with nephrolithiasis," Kidney International, vol. 72, no. 5, pp. 592-598, 2007.
 B. Gao, T. Yasui, Y. Itoh, K. Tozawa, Y. Hayashi, and K. Kohri, "A polymorphism of matrix Gla protein gene is associated with kidney stones," Journal of Urology, vol. 177, no. 6, pp. 2361-2365, 2007.
 D. F. Gudbjartsson, H. Holm, O. S. Indridason et al., "Association of Variants at UMOD with Chronic Kidney Disease and Kidney Stones--Role of Age and Comorbid Diseases," PLoS Genet, vol. 6, no. 7, p. e1001039, 2010.
 C. A. Wright, S. Howles, D. C. Trudgian et al., "Label-free quantitative proteomics reveals differentially regulated proteins influencing urolithiasis," Molecular & Cellular Proteomics, vol. 10, no. 8, p. M110.005686, 2011.
 K. Lage, N. T. Hansen, E. O. Karlberget al., "A large-scale analysis of tissue-specific pathology and gene expression of human disease genes and complexes," Proceedings of the National Academy of Sciences, vol. 105, no. 52, pp. 20870-20875, 2008.
 A. Taylor, J. Steinberg, T. S. Andrews, and C. Webber, "Gene-Net Toolbox for MATLAB: a flexible platform for the analysis of gene connectivity in biological networks," Bioinformatics, vol. 31, no. 3, pp. 442-444, 2015.
 P. Shannon, A. Markiel, O. Ozier et al., "Cytoscape: a software environment for integrated models of biomolecular interaction networks," Genome Research, vol. 13, no. 11, pp. 2498-2504, 2003.
 D. W. Huang, B. T. Sherman, and R. A. Lempicki, "Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources," Nature Protocols, vol. 4, no. 1, pp. 44-57, 2009.
 X. Lu, Y. Li, W. Wang et al., "3 [beta]-Hydroxysteroid-[DELTA] 24 reductase (DHCR24) protects neuronal cells from apoptotic cell death induced by endoplasmic reticulum (ER) stress," PLoS One, vol. 9, no. 1, p. e86753, 2014.
 J. Wu and R. J. Kaufman, "From acute ER stress to physiological roles of the unfolded protein response," Cell Death & Differentiation, vol. 13, no. 3, pp. 374-384, 2006.
 H. P. Harding, Y. Zhang, and D. Ron, "Protein translation and folding are coupled by an endoplasmic-reticulum- resident kinase," Nature, vol. 397, no. 6716, pp. 271-274, 1999.
 H. Yoshida, S. Nadanaka, R. Sato, and K. Mori, "XBP1 is critical to protect cells from endoplasmic reticulum stress: evidence from site-2 protease-deficient Chinese hamster ovary cells," Cell Structure and Function, vol. 31, no. 2, pp. 117-125, 2006.
 A. Abhishek, S. Benita, M. Kumari et al., "Molecular analysis of oxalate-induced endoplasmic reticulum stress mediated apoptosis in the pathogenesis of kidney stone disease," Journal of Physiology and Biochemistry, vol. 73, no. 4, article 587, pp. 561-573, 2017.
 M. S. J. Schepers, E. S. van Ballegooijen, C. H. Bangma, and C. F. Verkoelen, "Oxalate is toxic to renal tubular cells only at supraphysiologic concentrations," Kidney International, vol. 68, no. 4, pp. 1660-1669, 2005.
 V. N. Ratkalkar and J. G. Kleinman, "Mechanisms of stone formation," Clinical Reviews in Bone and Mineral Metabolism, vol. 9, no. 3-4, article 9104, pp. 187-197, 2011.
 Y. Ma, J. W. Brewer, J. Alan Diehl, and L. M. Hendershot, "Two distinct stress signaling pathways converge upon the CHOP promoter during the mammalian unfolded protein response," Journal of Molecular Biology, vol. 318, no. 5, pp. 1351-1365, 2002.
 H. Puthalakath, L. A. O'Reilly, P. Gunn et al., "ER stress triggers apoptosis by activating BH3-only protein Bim," Cell, vol. 129, no. 7, pp. 1337-1349, 2007.
 R. J. Davis, "Signal transduction by the JNK group of MAP kinases," Cell, vol. 103, no. 2, pp. 239-252, 2000.
 J. A. Wesson, R. J. Johnson, M. Mazzali et al., "Osteopontin is a critical inhibitor of calcium oxalate crystal formation and retention in renal tubules," Journal of the American Society of Nephrology, vol. 14, no. 1, pp. 139-147, 2003.
Baoyu Yang, (1) Xiuli Lu, (1) Yang Li, (1,2) Yuanyuan Li, (3) Daojun Yu, (2,4) Weiwen Zhang, (3) Chaojie Duan, (3) Kazumi Taguchi, (4,5) Takahiro Yasui [ID], (4,5) Kenjiro Kohri, (4,5) and Bing Gao [ID] (2,3,4)
(1) Department of Biochemistry and Cell Biology, School of Life Sciences, Liaoning University, Shenyang 110036, China
(2) Key Laboratory of Renal Calcification Disease Prevention and Treatment, Shenyang Medical College, 146 Huanghe North Street, Shenyang 110034, China
(3) Department of Cell Biology and Genetics, Shenyang Medical College, 146 Huanghe North Street, Shenyang 110034, China
(4) China-Japan Kidney Stone Research Center, Shenyang Medical College, 146 Huanghe North Street, Shenyang 110034, China
(5) Department of Nephro-Urology, Nagoya City University Graduate School of Medical Sciences, 1 Kawasumi, Mizuho-ku, Nagoya 467-8601, Japan
Correspondence should be addressed to Bing Gao; email@example.com
Baoyu Yang and Xiuli Lu contributed equally to this work.
Received 31 May 2019; Revised 18 August 2019; Accepted 27 August 2019; Published 27 October 2019
Academic Editor: Luciana Mosca
Caption: Figure 1: Network analysis of a kidney stone. A flowchart of network construction and a GO analysis (a). The primary kidney stone PPI network was constructed from protein sets obtained from PolySearch and a proteomics data by Wright et al. , which consisted of 340 nodes connected via 740 edges (b). The backbone kidney stone PPI network was derived from the primary kidney stone PPI network, which consisted of 31 nodes connected via 51 edges (c). Bigger nodes represent genes with more links.
Caption: Figure 2: Markers and key sensors of ER stress were activated in HK2 cells under COM crystal exposure. The Bip/GRP78 significantly increased from three to 48 hours (a, c), and the CHOP was increased from 12 to 24 hours, after exposure to COM crystals (a, d). Phospho-PERK (P-PERK) was increased to three hours (a, e). Cleaved ATF6 was increased from three to 48 hours (a, f). XBP-1 was decreased to 12 hours and then increased to 24 hours (a, g). The spliced XBP-1 mRNA was upregulated from three hours to 12 hours (b, h). * P < 0.05, ** P < 0.01, and *** P < 0.001.
Caption: Figure 3: The ER stress-specific caspase was involved in apoptosis induced by COM crystals. Apoptotic cells were detected by the TUNEL and Hoechst assay (a). COM crystal exposure significantly increased apoptotic cell numbers (*** P <0.001) compared with the control group (b). The ER stress-specific caspase-12, caspase-3, phospho-JNK (P-JNK), and phospho-P38 (P-P38) expression was all detected in a time-dependent manner using a western blot (c-g). *** P < 0.001, * P < 0.05.
Caption: Figure 4: The effect of ER stress on crystal-cell interaction. (a) ER stress promoted crystal-cell adhesion. In the TM+COM group, the HK2 cells were treated with COM crystals for five minutes after exposing TM, an inducer of ER stress, for 24 hours. The COM crystals that adhered to the cells were lysed with 5 ml of 6M HCl. A quantitative analysis of the COM crystals was conducted, by measuring the concentration using an atomic absorption method. Crystal adhesion in the TM+COM group rose significantly compared with the COM group (*** P < 0.001), and SAL reversed it (** P < 0.01). SAL also reduced crystal-cell adhesion after exposure to COM crystals (* P < 0.05). (b) ER stress induced by TM affected the expression of proteins associated with the formation of a kidney stone. Protein expression showed a similar trend between TM exposure and COM exposure. *, (#) P < 0.05, **, (##) P < 0.01, and ***, (###) P < 0.001 vs. 0 h. (c) ER stress induced by TM reduced cell viability, and SAL reversed it in vitro. HK2 cells were exposed to 100 [micro]g/ml COM crystals and 1 [micro]g/ml TM for two days, respectively. COM crystals and TM significantly reduced cell viability (* P <0.05). HK2 cells were exposed to 100 [micro]g/ml COM crystals and 1 [micro]M SAL+100 [micro]g/ml COM crystals for two days, respectively. SAL significantly increased cell viability compared with the group of COM crystals (*** P < 0.001). (d) The number of cells was compared to that of the group of COM and the group of COM+SAL for 12, 24, and 48 hours. The number of cells in the COM+SAL group significantly increased at 24 hours (* P < 0.05). TM represented tunicamycin. SAL represented salubrinal.
Caption: Figure 5: Marker proteins of ER stress expressed changes in rats treated with 0.75% EG from day zero to eight weeks. Crystals were formed from the 2nd week. The red arrow points to the crystal (a). The Bip/GRP78, CHOP, phospho-PERK (P-PERK), ATF6, caspase-3, and caspase-12 were activated (b). XBP-1 was spliced (c).
Table 1: Gene Ontology analysis for the backbone network of a kidney stone. Term Count % P value GOTERM_BP_DIRECT GO:0022617~extracellular 5 16.129 8.45E--06 matrix disassembly GO:0006457~protein folding 6 19.3548 1.33E-05 GO:0050821~protein 5 16.129 8.13E-05 stabilization GO:0036500~ATF6-mediated 3 9.67742 1.03E-04 unfolded protein response GO:0006950~response to 4 12.9032 1.64E-04 stress GO:0006898~receptor- 5 16.129 1.85E-04 mediated endocytosis GO:0034975~protein folding 3 9.67742 2.22E-04 in endoplasmic reticulum GO:1900034~regulation of 4 12.9032 2.88E-04 cellular response to heat GOTERM_CC_DIRECT GO:0070062~extracellular 24 77.4194 3.39.E-14 exosome GO:0042470~melanosome 8 25.8065 1.85E-10 GO:0005925~focal adhesion 11 35.4839 2.87E-10 GO:0031012~extracellular 10 32.2581 5.97E-10 matrix GO:0071682~endocytic 4 12.9032 2.01E-06 vesicle lumen GO:0005829~cytosol 18 58.0645 2.02E-06 GO:0005576~extracellular 13 41.9355 2.69E-06 region GO:0005913~cell-cell 7 22.5807 1.00E-05 adherens junction GO:0016020~membrane 14 45.1613 1.18E-05 GO:0072562~blood 5 16.129 7.86E-05 microparticle GO:0034663~endoplasmic 3 9.67742 1.34E-04 reticulum chaperone complex GO:0009986~cell surface 7 22.5807 2.14E-04 GO:0005788~endoplasmic 5 16.129 2.32E-04 reticulum lumen GOTERM_MF_DIRECT GO:0005515~protein binding 27 87.0968 4.47E-08 GO:0001948~glycoprotein 6 19.3548 9.67E-08 binding GO:0051082~unfolded 6 19.3548 1.32E-06 protein binding GO:0098641~cadherin 7 22.5807 1.24E-05 binding involved in cell-cell adhesion GO:0019899~enzyme binding 7 22.5807 2.62E-05 GO:0023026~MHC class II 3 9.67742 3.86E-04 protein complex binding KEGG_PATHWAY hsa05200: pathways in 12 38.7097 1.00E-07 cancer hsa04151: PI3K-Akt 10 32.2581 3.66E-06 signaling pathway hsa05215: prostate cancer 6 19.3548 1.93E-05 hsa04114: oocyte meiosis 6 19.3548 5.45E-05 hsa04621: NOD-like 5 16.129 5.50E-05 receptor signaling pathway hsa04390: Hippo signaling 6 19.3548 2.55E-04 pathway hsa04141: protein 6 19.3548 4.30E-04 processing in endoplasmic reticulum hsa05219: bladder cancer 4 12.9032 5.14E-04 hsa04915: estrogen 5 16.129 5.40E-04 signaling pathway Term Benjamini FDR GOTERM_BP_DIRECT GO:0022617~extracellular 0.00493 0.012443 matrix disassembly GO:0006457~protein folding 0.00388 0.019533 GO:0050821~protein 0.01573 0.119559 stabilization GO:0036500~ATF6-mediated 0.01495 0.151433 unfolded protein response GO:0006950~response to 0.01897 0.240659 stress GO:0006898~receptor- 0.01786 0.271641 mediated endocytosis GO:0034975~protein folding 0.01839 0.326434 in endoplasmic reticulum GO:1900034~regulation of 0.02084 0.422888 cellular response to heat GOTERM_CC_DIRECT GO:0070062~extracellular 4.40E-12 3.94E-11 exosome GO:0042470~melanosome 1.20E-08 2.16E-07 GO:0005925~focal adhesion 1.24E-08 3.34E-07 GO:0031012~extracellular 1.94E-08 6.95E-07 matrix GO:0071682~endocytic 5.22E-05 0.002338 vesicle lumen GO:0005829~cytosol 4.38E-05 0.002354 GO:0005576~extracellular 4.99E-05 0.003128 region GO:0005913~cell-cell 1.63E-04 0.011684 adherens junction GO:0016020~membrane 1.70E-04 0.013736 GO:0072562~blood 0.00102 0.091484 microparticle GO:0034663~endoplasmic 0.00158 0.155484 reticulum chaperone complex GO:0009986~cell surface 0.00232 0.249403 GO:0005788~endoplasmic 0.00232 0.269618 reticulum lumen GOTERM_MF_DIRECT GO:0005515~protein binding 8.09.E-06 5.51E-05 GO:0001948~glycoprotein 8.75E-06 1.19E-04 binding GO:0051082~unfolded 7.96E-05 0.001626 protein binding GO:0098641~cadherin 5.62E-04 0.015318 binding involved in cell-cell adhesion GO:0019899~enzyme binding 9.50E-04 0.032354 GO:0023026~MHC class II 0.01158 0.47498 protein complex binding KEGG_PATHWAY hsa05200: pathways in 1.36E-05 1.17E-04 cancer hsa04151: PI3K-Akt 2.49.E-04 0.004295 signaling pathway hsa05215: prostate cancer 8.76E-04 0.022701 hsa04114: oocyte meiosis 0.00185 0.063994 hsa04621: NOD-like 0.0015 0.064592 receptor signaling pathway hsa04390: Hippo signaling 0.00577 0.299238 pathway hsa04141: protein 0.00833 0.503915 processing in endoplasmic reticulum hsa05219: bladder cancer 0.0087 0.601459 hsa04915: estrogen 0.00813 0.632287 signaling pathway Term Genes GOTERM_BP_DIRECT GO:0022617~extracellular BSG, A2M, MMP9, CDH1, matrix disassembly FN1 GO:0006457~protein folding HSP90AB1, P4HB, HSP90B1, HSP90AA1, APCS, CALR GO:0050821~protein HSP90AB1, HSP90AA1, stabilization APOA1, CALR, GAPDH GO:0036500~ATF6-mediated HSP90B1, HSPA5, CALR unfolded protein response GO:0006950~response to HSP90AB1, MAPK1, HSP90B1, stress HSP90AA1 GO:0006898~receptor- HSP90B1, HSP90AA1, APOA1, mediated endocytosis CALR, CD14 GO:0034975~protein folding HSP90B1, HSPA5, CALR in endoplasmic reticulum GO:1900034~regulation of HSP90AB1, MAPK1, cellular response to heat HSP90AA1, YWHAE GOTERM_CC_DIRECT GO:0070062~extracellular HSP90AB1, P4HB, A2M, BSG, exosome APCS, HSP90AA1, MMP9, YWHAB, CDH1, ITGB2, CALR, YWHAE, PCK1, JUP, MAPK1, HSP90B1, APOA1, EIF4A1, YWHAQ, HSPA5, EGF, GAPDH, CD 14, FN1 GO:0042470~melanosome HSP90AB1, P4HB, HSP90B1, BSG, HSP90AA1, YWHAB, HSPA5, YWHAE GO:0005925~focal adhesion JUP, MAPK1, P4HB, HSP90B1, BSG, YWHAB, YWHAQ, CDH1, HSPA5, CALR, YWHAE GO:0031012~extracellular JUP, P4HB, HSP90B1, matrix HSP90AA1, APCS, EIF4A1, HSPA5, CALR, GAPDH, FN1 GO:0071682~endocytic HSP90B1, HSP90AA1, APOA1, vesicle lumen CALR GO:0005829~cytosol HSP90AB1, AR, A2M, HSP90AA1, EXOC7, YWHAB, CALR, TAB1, YWHAE, PCK1, JUP, MAPK1, HSP90B1, APOA1, EIF4A1, YWHAQ, HBG2, GAPDH GO:0005576~extracellular P4HB, A2M, BMP2, APCS, region HSP90AA1, MMP9, CDH1, CALR, HSP90B1, APOA1, EGF, CD 14, FN1 GO:0005913~cell-cell HSP90AB1, JUP, BSG, adherens junction YWHAB, CDH1, HSPA5, YWHAE GO:0016020~membrane HSP90AB1, HSP90B1, BSG, HSP90AA1, EXOC7, EIF4A1, YWHAB, YWHAQ, CDH1, ITGB2, HSPA5, CALR, GAPDH, YWHAE GO:0072562~blood A2M, APOA1, APCS, HBG2, microparticle FN1 GO:0034663~endoplasmic P4HB, HSP90B1, HSPA5 reticulum chaperone complex GO:0009986~cell surface HSP90AB1, BMP2, HSP90AA1, APOA1, ITGB2, HSPA5, CALR GO:0005788~endoplasmic P4HB, HSP90B1, APOA1, reticulum lumen HSPA5, CALR GOTERM_MF_DIRECT GO:0005515~protein binding HSP90AB1, A2M, EXOC7, MMP9, ITGB2, CDH1, CALR, APOA1, HSPA5, EGF, GAPDH, FN1, P4HB, AR, BSG, BMP2, CTBP2, HSP90AA1, YWHAB, TAB1, YWHAE, JUP, MAPK1, HSP90B1, EIF4A1, YWHAQ, CD14 GO:0001948~glycoprotein HSP90AB1, HSP90AA1, CDH1, binding ITGB2, HSPA5, CALR GO:0051082~unfolded HSP90AB1, HSP90B1, protein binding HSP90AA1, APCS, HSPA5, CALR GO:0098641~cadherin HSP90AB1, JUP, BSG, binding involved in YWHAB, CDH1, HSPA5, YWHAE cell-cell adhesion GO:0019899~enzyme binding P4HB, AR, A2M, APOA1, YWHAB, HSPA5, YWHAE GO:0023026~MHC class II HSP90AB1, HSP90AA1, YWHAE protein complex binding KEGG_PATHWAY hsa05200: pathways in HSP90AB1, JUP, MAPK1, AR, cancer HSP90B1, BMP2, HSP90AA1, CTBP2, MMP9, CDH1, EGF, FN1 hsa04151: PI3K-Akt HSP90AB1, MAPK1, HSP90B1, signaling pathway HSP90AA1, YWHAB, YWHAQ, EGF, YWHAE, FN1, PCK1 hsa05215: prostate cancer HSP90AB1, MAPK1, AR, HSP90B1, HSP90AA1, EGF hsa04114: oocyte meiosis MAPK1, AR, YWHAB, YWHAQ, PPP3R2, YWHAE hsa04621: NOD-like HSP90AB1, MAPK1, HSP90B1, receptor signaling pathway HSP90AA1, TAB1 hsa04390: Hippo signaling BMP2, YWHAB, YWHAQ, CDH1, pathway ITGB2, YWHAE hsa04141: protein HSP90AB1, P4HB, HSP90B1, processing in endoplasmic HSP90AA1, HSPA5, CALR reticulum hsa05219: bladder cancer MAPK1, MMP9, CDH1, EGF hsa04915: estrogen HSP90AB1, MAPK1, HSP90B1, signaling pathway HSP90AA1, MMP9
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Yang, Baoyu; Lu, Xiuli; Li, Yang; Li, Yuanyuan; Yu, Daojun; Zhang, Weiwen; Duan, Chaojie; Taguchi, K|
|Publication:||Oxidative Medicine and Cellular Longevity|
|Date:||Nov 1, 2019|
|Previous Article:||Kalanchoe brasiliensis Cambess., a Promising Natural Source of Antioxidant and Antibiotic Agents against Multidrug-Resistant Pathogens for the...|
|Next Article:||Dexmedetomidine Attenuates Orthotopic Liver Transplantation-Induced Acute Gut Injury via [[alpha].sub.2]-Adrenergic Receptor-Dependent Suppression of...|