Endoplasmic Reticulum Stress (ER Stress) and Unfolded Protein Response (UPR) Occur in a Rat Varicocele Testis Model
Using a surgically induced varicocele rat model, we show here strong evidence that the misfolded/unfolded protein response that is part of the stress response of the endoplasmic reticulum (ER) is activated in the varicocele testis (VCL), leading to the induction of apoptosis. To support this hypothesis, it is observed that the spliced variant of the X-box protein 1 (XBP1s), resulting from the activation of the inositol-requiring enzyme 1 (IRE1) membrane sensor, is significantly more represented in VCL testicular extracts. The activation of the IRE1/XBP1s pathway is also supported by the observation that the VCL testes show an increase phosphorylation of the c-Jun-kinase (JNK) known to be one intermediate of this pathway and an increased level of caspase-3, the terminal apoptotic effector, partly explaining the apoptotic status of the VCL testis.
Varicocele (VCL) is a large varicose dilatation of the spermatic cord veins located in the bursa above and around each testicle . With approximately 15% of the adult male population affected, varicocele is a relatively common disease that can occur at any age. The prevalence of varicocele in men with primary and secondary infertility is estimated to be 35% and 70-85%, respectively, making it the most common cause of male infertility [2, 3]. VCL sets in gradually and is accompanied by adverse effects on spermatogenesis, testicles, and spermatozoa. Although it is clear that higher scrotal temperature, inflammation, hypoxia, nutrient deprivation, oxidative stress, and increased apoptosis are characteristics of the varicocele testis, the full spectrum of the dysfunctional organ is not yet well defined . In this complex picture, the increased generation of reactive oxygen species (ROS) recorded in the damaged and dysfunctional VCL testis and their consequences on sperm structure and function are the most documented facets [5–8].
Hyperthermia, inflammation, oxidative stress, and nutrient deprivation are situations that are also known to promote stress of the cellular endoplasmic reticulum (ER stress) characterized by the unfolded protein response (the so-called UPR/ER stress response) [9–13]. In this particular stress response, an accumulation of misfolded/unfolded proteins in the ER of defective cells is encountered, resulting in the activation of an ER-associated protein degradation pathway (ERAD)  which will lead to ER-mediated apoptosis . In this context, the immunoglobulin binding protein (BiP), also known as the GRP78/HSPA5 protein linked to glucose (a member of the heat shock protein family), is a key player that orchestrates UPR and ER stress responses to delay cell death . Briefly, and as illustrated in Figure 1, the chaperone BiP/GRP78/HSPA5 blocks the ER stress response through its interaction with ER membrane sensor proteins including PERK, IRE1, and ATF6 [14, 16, 17]. When BiP is driven out of these membrane sensors by an excess of unfolded/misfolded ER proteins, ER stress signaling begins  and the 3 pathways are transiently activated concurringly to adjust cell survival and apoptotic signals (see Figure 1). When the UPR response fails, as may occur under established and/or prolonged ER stress, apoptosis triggered by specific inducers such as the homologous C/EBP protein known as CHOP transcription factor, the proapoptotic factors Bim, Bax, Bak, and the Jun-kinase/caspase-3 pathway wins (Figure 1) [14–17].
In the mammalian reproductive tract, ER chaperones were shown to be important for spermatogenesis as well as for posttesticular sperm maturation . In drosophila, ER stress in the male accessory organ was also reported to result in infertility . In addition, the UPR/ER-stress responses were shown to occur in mammalian testis following chemical or physical stresses [20–23]. These observations prompted us to suspect that the UPR/ER stress response could be part of VCL-induced testicular dysfunction. We therefore have chosen here to evaluate markers of UPR/ER pathways in a VCL-induced rat model.
2. Materials and Methods
This study was approved by the Ethics Committee of Royan Institute (IR. ACECR.ROYAN.REC.1397.224), and animal care was performed according to the Animal Ethics Committee of Royan Institute recommendations.
2.1. Chemicals and Antibodies
The CHOP mouse antibody (Cat. No. 2895) and phospho-JNK mouse antibody (Cat. No. 9255) were purchased from Cell Signaling Technology (Danvers, MA, USA). Antibodies against glutathione peroxidase 4 (GPx4: ab125066) and nuclear factor erythroid-related transcription factor 2 (NRF2: ab62352) were purchased from Abcam (Cambridge, UK). The mouse anti-actin antibody (Cat. No: A2228), goat and rabbit immunoglobulin G (IgG) coupled with horseradish peroxidase (HRP-Cat. No: Sc-2301), and goat and mouse IgG-HRP (Cat. No: P0447) were purchased from Sigma (Saint-Louis, MI, USA), Santa-Cruz (Dallas, TX, USA), and Dako-Agilent (Santa Clara, CA, USA), respectively, while the TUNEL kit was from Promega (Madison, WN, USA). The chromomycin A3 and acridine orange dyes were purchased from Sigma (Saint-Louis, MI, USA).
2.2. Experimental Design and Organ/Sperm Cell Collection
Twenty 8-week-old adult male Wistar rats (weight range: 180-220 grams) were obtained from the Institute of Biotechnology (Isfahan, Iran). For the induction of varicocele, rats were anaesthetised by intraperitoneal injection of ketamine (75 mg/kg) combined with xylazine (2.5 mg/kg) according to Molina et al. , and VCL was induced as originally described by Ko et al.  and Afiyani et al. . Briefly, a median incision of 2 cm was made and the diameter of the left renal vein was reduced to 1 mm; ligation of the left renal vein was performed midline at the junction of the adrenal and spermatic veins. Then, the anastomotic branch between the left testicular vein and the left common iliac vein was ligated with a 0-4 silk suture. Ten rats formed the varicocele-induced group, while ten other rats formed the control group. In this study, we did not form a sham group since on several occasions we have demonstrated that with the protocol used there were no differences between the sham group and the control group [26–28]. Rats were maintained in a controlled environment at 22-26°C, with 55-60% humidity and a 12-hour light/dark cycle. Two months after VCL induction, the body weight of the animals was measured and the animals were sacrificed as described earlier . The testes were evaluated based on their morphometric characteristics, including length, width, weight, thickness, and volume. The left testicle of each rat was divided into two parts; one part was fixed in a Bouin solution for histopathological evaluation while the other part was cryopreserved and stored prior to evaluation of the genes, proteins and molecules of interest.
The epididymides were also collected and divided into three parts. Sperm cell collection for evaluation of structural and functional spermatozoa parameters (sperm count, motility, morphology, lipid peroxidation, DNA damage, histone, and protamine content) was performed on epididymal tails (cauda) after chopping the tissue into small pieces which were then incubated at 37°C for 30 minutes in 1 ml VitaSperm medium (Inoclon, Iran).
2.3. Histological Analyses
Testicular sections fixed by Bouin were evaluated after staining with haematoxylin/eosin (HE) and parameters such as mean percentage of tubular diameter, epithelial height, meiotic index, percentage of spermatogenesis, and Johnsen’s score  were evaluated using an Olympus CX31 microscope (40x magnification). The TUNEL stain was used as previously indicated for the evaluation of apoptosis .
2.4. Spermatozoa Structural and Functional Evaluation
For the evaluation of spermatozoa mobility, 10 μl of sperm cell samples retrieved from the cauda epididymides (see above) were smeared on a slide and observed using an optical microscope (CX31 Olympus, Dubai, EAU). For each animal, at least 200 sperm cells were evaluated from different fields, and the mean percentage of mobile sperm was recorded. Sperm counts were estimated using a Makler counting chamber (sperm meter; sperm processor, Garkheda, Aurangabad, India) using a LABOMED CxL optical microscope (magnification: 20x) and data were expressed as million cells per ml. For morphology, 50 μl of cauda-retrieved spermatozoa samples were stained using eosine/nigrosine according to Afiyani et al. . The percentages of abnormal sperm morphology presented for each sample include abnormal head, neck, and flagellum. Sperm membrane lipid peroxidation level was monitored using BODIPY 581/591 C11 (Cayman Chemicals, Ann Arbor, MI, USA) as described earlier by Aitken et al. . Briefly, 2 million sperm cells were mixed with the BODIPY probe (5 mM final concentration) in the presence (positive control) or absence of H2O2 (10 mM). Controls and test tubes were then incubated for 30 min at 37°C. Samples were then washed twice using a phosphate-buffered saline (PBS) solution (1X) at 2000 rpm for 5 min at room temperature. A flow cytometer (FACSCalibur, Becton Dickinson, San Jose, CA, USA) was used to assess the percentage of sperm cells showing lipid peroxidation. For sperm DNA damage, the acridine orange (AO) staining method was used . Aliquots of 20 μl of cauda-retrieved sperm cells were smeared on slides and fixed with Carnoy’s solution at 4°C overnight. Then, the slides were stained using 150 μl of a freshly prepared AO solution (in citric acid 0.1 M, NaH2PO4 0.3 M, pH: 2.5) for 90 min in the dark at room temperature. The slides were then washed with PBS (1X) twice and observed using a fluorescent microscope (OPTIKA, Ponteranica, Italy). Two hundred sperm cells were counted on each slide, and the percentage of DNA damaged cells (cells having a red/orange nucleus) was calculated. For sperm chromatin integrity, both aniline blue (AB) and chromomycin A3 (CMA3) staining were used. AB indirectly evaluates the level of condensation of the sperm nucleus via its affinity to lysine-rich histones while CMA3 competes with protamines [33, 34]. To perform the AB staining, washed sperms were fixed with 3% glutaraldehyde in 0.2 M phosphate buffer () for 2 hours at 4°C. Then, two smears were prepared and air-dried for each condition. Aqueous AB (5%) in acetic acid (4%) was added to the slides for 5 minutes and the slides were washed once with phosphate-buffered saline (PBS 1X). The slides were observed using an optical microscope (OPTIKA, Ponteranica, Italy, magnification 100x). For each slide, 200 spermatozoa were counted, and the percentage of spermatozoa with immature chromatin (dark blue stained nucleus) was calculated (Figure 2). To perform the CMA3 staining, 30 μl washed sperm were fixed with 30 μl Carnoy’s solution for 5 minutes at 4°C. The slides were then rinsed with PBS (1X) and stained with CMA3 solution (0.25 mg/ml). The slides were then rinsed and covered. For each slide, 200 spermatozoa were counted using a fluorescent microscope (OPTIKA, Ponteranica, Italy, magnification 100x), and the percentage of CMA3-positive spermatozoa was calculated (Figure 2).
2.5. Total RNA Extraction and RT-PCR Analyses
In short, 50 mg of frozen testis tissues was used to extract RNA using the TRIzol reagent (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer’s protocol. RNA quantification was done by evaluating the absorbance ratio A260/A280. Complementary DNA was prepared using a commercial reverse transcription kit (Takara Bio Inc., Shiga, Japan). The qRT-PCR were performed in a total volume of 10 μl containing 1 μl cDNA, 2.8 μl H2O, 0.2 μl Rox, 1 μl specific primers, and 5 μl SYBR Green using a StepOnePlus thermal cycler (Applied Biosystems, Abi, Foster City, CA, USA). The 2-ΔΔCT method was used for data analysis. Primers were designed using the software Beacon Designer 7 and are presented in Table 1.
2.6. Western Blot Analyses
Testicular protein extracts (30 μg) were loaded on 12% SDS PAGE gel subsequently transferred on a PVDF membrane (Merck, Darmstadt, Germany). The membranes were then blocked and incubated with the different primary antibodies including: anti-phospho-SAPK/JNK (dilution: 1/1000), anti-GADD153/CHOP (dilution: 1/2000), anti-GPx4 (dilution: 1/200), anti-β-actin (dilution: 1/500), and anti-NRF2 (dilution: 1/500) overnight at 4°C. The membranes were then washed 3 times and appropriate secondary antibodies (anti-rabbit/anti-goat/anti-mouse IgG-HRP conjugated) depending on the primary antibody were used for 2 hours at room temperature. The ECL detection kit (Thermo Fisher Scientific, Waltham, MA, USA) was then used following the supplier recommendation. The β-actin signal was used to normalize sample loading differences. The ImageJ software (version 1.42q) was used to quantify the different signals.
2.7. Statistical Analysis
The statistical analysis was performed using IBM SPSS Statistics 25.0 and GraphPad Prism (version 7.00) for graphic design (GraphPad Software, San Diego, California, USA). The was represented for all data. An independent -test was used to compare the study tests between the control and varicocele groups.
3.1. Testis and Spermatozoon Characteristics in VCL-Induced Rats
The induction of VCL was first evaluated in the group of induced animals. The VCL induction protocol did not significantly alter the physiology of the animals, as evidenced by the fact that the mean body weight of the VCL and control animals was not significantly different (Table 2). As expected in a situation of VCL, the left testis volume reduction was found to be significant (Table 2). As expected also, and consistent with previous data, Table 3 shows that VCL induction resulted in a decrease in mean testicular tubule diameter and epithelium height associated with a decrease in the meiotic index, the percentage of spermatogenesis, and the Johnsen score.
Figures 3(a)–3(f) illustrate the classic structural defects observed in VCL testes, including the destructuring of the seminiferous epithelium, vacuolization, increased connective tissue, and Sertoli and Leydig cell abnormalities. In addition, Figure 3(g) shows that VCL testicular sections contain significantly more apoptotic cells than control testicular sections. Considering the effect of VCL induction on sperm parameters, we show in Figure 4(a) that the mean sperm concentration and motility were significantly decreased in semen samples taken from the tail of the epididymis of VCL animals. At the same time, we recorded a small but significant increase () in abnormal sperm morphologies in VCL animals. Upon closer examination of the integrity of the sperm nucleus in VCL animals, we recorded clear situations of protamine deficiency associated with increased levels of persistent histones and a very significant increase in sperm DNA damage, as shown by acridine orange staining (Figure 4(b)). These abnormalities reflecting spermatozoa nuclear immaturity and/or damage were also accompanied by an increased level of peroxidation of sperm membrane lipids, as shown by the higher BODIPY reactivity of VCL sperm samples (Figure 4(b)).
To assess the well-known situation of oxidative stress associated with VCL, we used real-time qPCR to compare control and VCL testes. We show in Figure 4(c) that the accumulation of NRF2 mRNA (a major transcription factor orchestrating the antioxidant response under oxidative stress) is significantly increased in VCL testes () compared to control testes () (). However, the increased accumulation of NRF2 mRNA in the VCL testes does not result in an increase in NRF2 protein content; this will be discussed later.
3.2. Assessment of the UPR/ER-Stress and Apoptotic Pathways in VCL Testis
To evaluate the induction of UPR/ER stress pathways in VCL testes, we first chose to monitor the level of the ER chaperone protein BiP/GRP78/HSAP5 using a qRT-PCR approach. We found no difference in the accumulation of BiP transcripts when we compared VCL testes extracts with control extracts (Figure 5(a)). We then examined the spliced form of the XBP1s protein, the transcription factor that performs UPR . As shown in Figure 1, the splicing of XBP1 mRNA is one of the benchmarks of UPR/ER stress response via the IRE1 membrane-associated protein . IRE1 oligomerizes and activates its ribonuclease domain by autophosphorylation. It then catalyzes the excision of a 26 bp intron on the mRNA of the ubiquitously expressed XBP1 protein (XBP1u). This causes a shift in the XBP1s coding sequence, resulting in the translation of a distinct XBP1s protein. Using this specific feature, we show in Figure 5(b) that the spliced form of the XBP1 transcript (XBP1s) is 6 times more represented in VCL testis extracts than in control testis extracts ( versus ; ), confirming the activation of the UPR pathway in the rat testis of VCL.
As indicated above (Figure 1), when the UPR response is prolonged, IRE1 stimulates the JNK/p-JNK pathway which induces apoptosis  via, in particular, downregulation of the antiapoptotic protein Bcl2 and upregulation of mitochondrial-associated proapoptotic actors including Bax and Bak, leading to induction of caspase-3. We show in Figure 5(c) that in the VCL testis, the JNK pathway is activated as revealed by the increased level of phospho-JNK (). We also show in Figure 5(d) that caspase-3 transcript accumulation is significantly greater in the VCL testis compared to controls ( versus , ). To further evaluate the apoptotic network in the VCL testis, we also monitored the accumulation of two classical proapoptotic factors associated with mitochondria, Bax and Bak (Figure 5(e)). As shown in Figure 5(d), the accumulation of Bax and Bak mRNA did not increase significantly in the VCL testis.
Apart from the proapoptotic IRE1/JNK pathway, the accumulation of unfolded/misfolded proteins in the ER has the potential to trigger the PERK/ATF4/CHOP pathway, which should contribute to proapoptotic signals (see Figure 1). We show in Figure 5(f) that VCL testicular extracts did not show a higher CHOP content when compared to control testis.
To test whether VCL testis could be concerned by UPR/ER stress mechanisms, we used our expertise in VCL-induced rat models [24–28]. After verifying that our VCL-induced animals had all the classical characteristics expected in terms of testicular damage and sperm distorted parameters [28, 38–40], we then focused our research on marker characteristics of UPR/ER stress pathways.
Although BiP/GRP78 expression was shown to be modulated in the rodent testis throughout spermatogenesis , in our rat model of surgically induced VCL, we show that there is no difference in testicular BiP content. There might be a rationale for that observation as it was shown that the expression level of BiP is mainly regulated at the posttranscriptional level through the control of protein degradation processes that may readily allow BiP to return to optimal levels . Nevertheless, despite this observation that could suggest that there is no-ER-stress response ongoing in the VCL testis, we further investigated the downstream effectors of the UPR/ER stress response starting with the XBP1s spliced variant mediated by IRE1 activation [43, 44]. We clearly show that the VCL testis extracts contain a large excess of XBP1s, which is indicative of the activation of the IRE1 membrane sensor. Our observation that the downstream proapoptotic JNK/phospho-JNK switch [45, 46] is also mobilized, as shown by the increased representation of phospho-JNK in testicular VCL extracts, confirms the activation of the IRE1 ER sensor. Indeed, it was reported that if processing of XBP1 mRNA by the RNase domain of IRE1α promotes survival of ER stress, activation of the mitogen-activated protein kinase JNK by IRE1 late in the ER stress response promotes apoptosis . Consistent with the activated IRE1 proapoptotic pathway in VCL testes, we logically found an increased level of caspase-3, which partly explains the apoptotic observations made on sections of VCL testes via the TUNEL assay (Figure 3). Positive regulation of caspase-3 in the rat VCL testis has been reported previously . Similarly, studies have reported IRE1-mediated apoptosis in germ cell lines following various stresses [22, 49–51]. Our data corroborates also earlier data reported elsewhere .
The second pathway by which UPR/ER stress can modulate apoptosis is the PERK/ATF4/NRF2/CHOP pathway. On the one hand, PERK activity induces the expression of the homologous protein of the C/EBP transcription factor (CHOP) (see Figure 1). We show here that the CHOP-dependent PERK pathway does not appear to be preponderant because CHOP levels do not change in the VCL testis. Consistent with the observation that CHOP level remains unchanged, neither Bim nor Bcl2 levels were changed in the VCL testis (Figure 5(e)). Normally, in a situation where CHOP is activated, one would expect downward regulation of antiapoptotic Bcl2 and upward regulation of proapoptotic Bim proteins . None of these events occur in the VCL testis. In addition, the mitochondrial-associated proapoptotic effectors, Bax and Bak, which are expected to partially respond to Bim, were not found to be altered in the VCL testis. On the other hand, PERK activity allows the nuclear translocation of NRF2 where it increases the expression of antioxidant genes . We cannot answer the question of whether there is more nuclear NRF2 in VCL testis cells because we have not specifically followed up on this question. We only show that in VCL testis extracts, we record a transcriptional induction of NRF2 that was not accompanied by an increase in NRF2 protein content or positive regulation of a NRF2 target gene (the glutathione peroxidase 4: GPx4 gene, not shown). This is rather difficult to explain, although a similar situation has been reported in chronic liver disease UPR/ER stress response where an increase in XBP1s has been shown to be associated with increased NRF2 degradation . Therefore, it is possible that the recorded increase in XBP1s in the VCL testes may mask the increase in NRF2 that would have been expected otherwise.
The ATF6 channel is the third branch through which it is commonly accepted that the UPR/ER stress response is mediated . Again, the sequestration of BiP/GRP78/HSPA5 by misfolded/unfolded ER proteins activates the ATF6 transcription factor. Early activation of ATF6 in the ER-stress response is expected in return to activate the transcription of BiP providing more chaperone aiming at alleviating the ER stress . It was recently reported that ATF6 exerts also its action through CHOP which, in turn, regulates the decision on the fate of the cells . Since neither BiP nor CHOP expression was shown upregulated in our model of VCL testis, at the time we evaluated the tissue response, it may suggest that we are probably late in the UPR/ER stress response kinetics as hypothesized elsewhere [9, 58, 59].
In conclusion, we show here that the UPR/ER response is most likely activated in the VCL testis as evidenced by the engagement of the IRE1/JNK/XBP1s pathway. In addition, caspase-3-mediated apoptosis, which is known to be the ultimate response to UPR/ER stress, is also activated in the VCL testis. The observation that the PERK/ATF4/ATF6/CHOP pathway does not appear to be activated in the VCL testis can be explained. Several authors have recently shown that the 3 pathways known to mediate the stress response of the UPR/ER work with a specific kinetics, with PERK being the first to be activated, then ATF6, and finally IRE1 [9, 58, 59]. Moreover, it is proposed that when IRE1 is activated, the XBP1s variant has been shown to exert negative control over the PERK channel. Therefore, our interpretation of the situation in the VCL testis rat model we studied is that we are likely in the late stages of the UPR/ER stress response as it is strongly suggested by the activation of the p-JNK-induced caspase-3 apoptotic signal. Consistent with our data and interpretation is the very recent publication that the ER stress markers p-JNK, p-IRE, cleaved-caspase-3, and Bax are induced in rat VCL testes, a situation that could be reversed by oral administration of anti-inflammatory and antioxidant herbal compounds . Our data also support Kim et al. , who indicated that IRE1 is the primary mediator of the stress response of the UPR/ER in testicular hyperthermia in mice. As already suggested by the recent work of Karna et al., the data collected so far in animal VCL models opens the way for additional therapeutic actions to reduce the impact of human VCL by means of antioxidant and anti-inflammatory molecules .
All reported data (with the exception of the GPx4 mRNA accumulation) are presented in the manuscript.
Conflicts of Interest
The authors state that there is no conflict of interest that could be perceived as prejudicial to the impartiality of this manuscript.
M.H.N.E. is responsible for the study design, data analysis and interpretation, manuscript writing, and final approval of the manuscript. M.T. is responsible for the study design, data collection, assembly and analysis, manuscript writing, and final approval of the manuscript. J.R.D. copy-edited and revised the manuscript before final approval.
This study was supported by the Royan Institute, Iran. This research did not receive any specific support from commercial/private sector funding bodies. We would like to express our gratitude to the staff of the Biotechnology Department of Royan Institute for their full support.
P. D. Kantartzi, Goulis ChD, G. D. Goulis, and I. Papadimas, “Male infertility and varicocele: myths and reality,” Hippokratia, vol. 11, no. 3, pp. 99–104, 2007.View at: Google Scholar
Z. C. Wang, J. F. Wang, Y. B. Li et al., “Involvement of endoplasmic reticulum stress in apoptosis of testicular cells induced by low-dose radiation,” Journal of Huazhong University of Science and Technology. Medical Sciences, vol. 33, no. 4, pp. 551–558, 2013.View at: Publisher Site | Google Scholar
K. W. Ko, J. Y. Chae, S. W. Kim et al., “The effect of the partial obstruction site of the renal vein on testis and kidney in rats: is the traditional animal model suitable for varicocele research,” Korean Journal of Urology, vol. 51, no. 8, pp. 565–571, 2010.View at: Publisher Site | Google Scholar
B. X. Du, J. Z. Zheng, B. Liu et al., “Tongjingling improves sperm DNA integrity and reduces oxidative stress in the testis of experimental varicocele rats,” Zhonghua Nan Ke Xue, vol. 23, no. 11, pp. 1025–1031, 2017.View at: Google Scholar
J. Zhang, N. Singh, K. S. Robinson-Taylor et al., “Hepatocyte autophagy is linked to C/EBP-homologous protein, Bcl2-interacting mediator of cell death, and BH3-interacting domain death agonist gene expression,” Journal of Surgical Research, vol. 195, no. 2, pp. 588–595, 2015.View at: Publisher Site | Google Scholar
K. K. Karna, B. R. Choi, J. H. You et al., “The ameliorative effect of monotropein, astragalin, and spiraeoside on oxidative stress, endoplasmic reticulum stress, and mitochondrial signaling pathway in varicocelized rats,” BMC Complementary and Alternative Medicine, vol. 19, no. 1, p. 333, 2019.View at: Publisher Site | Google Scholar
K. K. Karna, B. R. Choi, M. J. Kim, H. K. Kim, and J. K. Park, “The Effect of Schisandra chinensis Baillon on Cross-Talk between Oxidative Stress, Endoplasmic Reticulum Stress, and Mitochondrial Signaling Pathway in Testes of Varicocele-Induced SD Rat,” International Journal of Molecular Sciences, vol. 20, no. 22, p. 5785, 2019.View at: Publisher Site | Google Scholar