Abstract
The impact of depression on spermatogenic function and protein expression was investigated by reference to a bioinformatics database for the prediction of key targets and pathways. Experimental validation was performed using a rat model established using the chronic restraint stress method. Organ coefficients and semen parameters were measured. Hematoxylin and eosin (HE) staining was performed to compare testicular tissue between model and control rats. Western blotting and RT–qPCR were conducted to evaluate protein and mRNA expression. As a result, depression and male infertility (MI) were found to share 81 common targets, and the PI3K/Akt/FOXO1 signaling pathway was involved in both conditions. Animal experiments showed testis and epididymis coefficients to be lower in the model group. HE staining showed damage to the seminiferous tubules of model rats. PI3K and p-Akt mRNA and protein were present at lower levels and FOXO1 at higher levels in the model group. Depression led to reduced spermatogenic function in rats and might be associated with differential expression of the PI3K/Akt/FOXO1 signaling pathway.
1. Introduction
Depression may be defined as a continuous and prolonged period of low mood, accompanied by symptoms such as loss of appetite, insomnia, severe fatigue, and a lack of sense of self-worth [1–3]. Epidemiological surveys have found that the lifetime prevalence of depression is 15%-18% [4], and 350 million people worldwide may suffer from the condition according to the World Health Organization (WHO) [5]. Indeed, depression is predicted to be a leading contributor to global disease burden by 2030 [6] and has been associated with diabetes, hypertension, and male infertility [7–9].
Male infertility (MI) refers to the inability of a couple to conceive naturally due to factors associated with the male partner after having regular sex for more than one year without the use of any contraception [10]. The global prevalence of infertility has been estimated at 7% with male factors accounting for up to 40% of cases [11]. 42.9% of male patients seeking fertility treatment have been diagnosed with depression [12]. And male patients are significantly more likely to experience depression after the failure of assisted reproductive technology (ART) than beforehand [13]. Men diagnosed with depression also had significantly lower semen volume, sperm count, and total motility than healthy men [14].
Many studies have sought to identify mechanisms which would explain the association between MI and depression. Chronic depression impaired the reproductive system through interactions between the hypothalamic–pituitary–adrenal (HPA) and hypothalamic–pituitary–gonadal (HPG) axes [15]. Moreover, depression could reduce sexual function in men, which in turn damage male fertility and exacerbate negative moods in males [16]. Although the use of antidepressants may improve depressive symptoms, the harmful effects on male reproductive function had gradually been discovered. Antidepressants have been reported to reduce sperm quality [17], induce sperm DNA fragmentation [18], and lower testosterone levels [19]. In summary, patients with MI are at higher risk of depression, but underlying mechanisms remain unclear.
Vast quantities of information relating to human health and biodiversity are available from biological databases and may be utilized to aid the understanding of the underlying biological mechanisms [20]. Therefore, the current study used bioinformatics analysis to predict the protein targets associated with depression and MI. Depression may affect the expression of the PI3K/Akt/FOXO1 pathway and related proteins, and these may also be associated with MI. Then, a rat model of depression was established, and spermatogenic function in testicular tissues and changes in signaling pathway proteins (PI3K, Akt, and FOXO1) were evaluated. The flow chart of the study design is shown in Figure 1.

2. Materials and Methods
2.1. Experimental Animals
Twenty male SPF Sprague–Dawley (SD) rats (aged 4-5 weeks; 200-220 g) were purchased from Beijing Vitonliol Laboratory Animal Technology Co., Ltd. (Animal License No. SCXK (Beijing) 2020-0033). Rats were housed in the SPF animal house of Beijing University of Chinese Medicine with 10-12 h ambient sunlight, 55-60% indoor humidity, and a temperature of 21-25°C. Animals were fed a solid diet and deionized water for 7 days of adaptation. Ethical approval was granted by the Experimental Animal Ethics Committee of Beijing University of Chinese Medicine (ethics number: BUCM-4-2020122303-4150), and guidelines for the ethical review of laboratory animal welfare under People’s Republic of China National Standard GB/T 35892-2018 were followed.
2.2. Drugs and Reagents
All reagents were purchased from GenePool (Beijing, China), including a sodium dodecyl sulfate (SDS) polyacrylamide gel electrophoresis (PAGE) Gel Kit (GPP1816), SDS–PAGE Loading Buffer (5x) (GPP1820), Protein Extraction Kit (GPP1815), Tris–Glycine Running Buffer (5x) (GPP1821), Total RNA Extraction Kit (DNase I) (GPQ1801), mRNA cDNA Synthesis Kit (GPQ1803), mRNA/lncRNA qPCR Kit (GPQ1808), and RNA Loading Buffer (5x) (GPQ181). Anti-PI3K antibody (1: 1000 dilution, Abcam, AB191606, Cambridge, UK) and GAPDH (1: 5000 dilution, AB181602) were purchased from Abcam; anti-p-Akt antibody (1: 1000 dilution, Bioss, bs-2720R, Beijing, China) and anti-FOXO1 antibody (1: 1000 dilution, Bioss, bs-2537R) were purchased from Beijing Bioss Biotechnology Co., Ltd. An electrophoresis instrument (CAVOY, PP-1150), spectrophotometer (Nano-Drop 2000, Thermo Scientific), real-time PCR instrument (LineGene 9600 Plus, Bioer Technology), and dual vertical electrophoresis tank (CAVOY, MP-8001) were also used.
2.3. Bioinformatics Analysis
2.3.1. Identification of Molecular Targets
Targets related to depression and MI were identified from the GeneCards database (https://www.genecards.org/), and a Venn diagram was constructed to indicate common targets. Network construction and analysis were performed based on the common targets.
2.3.2. Construction of the Protein–Protein Interaction (PPI) Network
Common targets were imported into the STRING database (https://string-db.org/, ver. 11.0), and a minimum interaction score value of 0.700 was selected as the screening criterion. Then, the PPI results were imported into Cytoscape software for visualization. The top ten proteins with the highest core degrees were analyzed by the CytoHubba plug-in and regarded as potential key targets.
2.3.3. GO and KEGG Pathway Enrichment Analyses
Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed on common target proteins using the DAVID database (https://david.ncifcrf.gov/) [21]. GO functional analysis was described in terms of biological processes (BPs), cellular components (CCs), and molecular functions (MFs). KEGG enrichment analysis revealed the signaling pathways enriched by the above targets.
2.4. In Vivo Experiments
2.4.1. Establishment of the Animal Model
Twenty SD rats were randomly divided by a random number table into control () and depression model () groups. Chronic restraint stress was used to establish the model of depression [22]. Model rats were placed inside a thin wire mesh to restrict movement on experimental day 1, and care was taken to adjust the mesh to ensure that rats could breathe freely. Restraint was imposed from 9:00 to 15:00 daily, during which time rats were not allowed food or water. At the end of the restraint period, the rats were returned to a single cage with food and water ad libitum. The modeling process lasted for 28 days. Control animals were reared normally with food and water ad libitum.
2.4.2. Behavioral Testing
The sucrose preference and open field tests were conducted for all animals, and significant differences between control and model groups () were considered to show the successful establishment of the model [23].
(1) Sucrose Preference Test (SPT). Rats were trained to drink sugary water prior to the experiment, and two bottles of 1% sucrose solution were given for the first 24 h. Thereafter, one bottle was replaced with sterile water. After adaptation, rats were deprived of water and food for 24 hours before being given access to 2 bottles of water: one of 1% sucrose and one of sterile water. Bottles were randomly placed on the left and right and weighed after 1 h. The sucrose preference ratio was determined using the following formula: .
(2) Open Field Test (OFT). The OFT device included a field reaction box () with the bottom divided into 25 equal compartments, and experiments were conducted in a dark, quiet room. Rats were placed in the reaction box, and the number of times they crossed the bottom (number of squares into which all four paws entered) was used to score horizontal movements. The number of times they stood upright (1 point for both feet leaving the bottom until they dropped) was used to score vertical movements. The experimental area was cleaned with 75% ethanol after each rat completed the test to remove odors. Two trained observers assessed rat locomotion over a period of 3 min before modeling and 1 day before sampling to avoid subjective errors.
2.4.3. Preparation of Tissue Samples
Rats were weighed and anesthetized by intraperitoneal injection of 50 mg/kg 1% sodium pentobarbital before sacrifice. Testes and epididymal tissues were removed and weighed using a precision electronic balance. Organ coefficients (testis or epididymis) were calculated as follows: .
The epididymis was transferred to prewarmed 0.9% saline, tissues minced with fine scissors, and incubated to assess semen quality (sperm concentration and viability). The right testes were used for Western blotting (WB) and reverse transcription polymerase chain reaction (RT–qPCR). The left testes were used for morphological observations by hematoxylin and eosin (HE) staining.
2.4.4. Analyses of Semen Quality
The epididymis was incubated in 0.9% saline at 37°C for 5 min, and sperm quality was analyzed by the Weili Sperm Quality Testing System (Beijing Xing Rong Technology Co., Ltd., China, No. WLJY-9000) at 37°C. Aliquots (5 μL) of semen were dropped onto the counting plate and covered with a coverslip for microscopic observations according to WHO standards [24]; 10 fields of view were randomly selected, and the number of spermatozoa in 10 large squares was counted. Mean values are presented as ×106 sperm per ml.
2.4.5. Histological Analysis
The left testis () was removed and cut along the largest surface. Half the sample was fixed with 4% paraformaldehyde for 24 h, trimmed and cut along the transverse axis, and fixed for 12 h. The sample was rinsed with running tap water for 24 h with deionized water for 2 h and dehumidified with absorbent paper before dehydration and degreasing with alcohol and xylene. Samples were embedded in paraffin, and 4-5 μm sections were cut. Three slices were taken from each testis, 10 fields of view were randomly selected for each slice, HE staining was performed, and histopathological changes were observed by electron microscopy (BK-FL4, OPTEC; Chongqing, China).
2.4.6. Western Blotting Assay
A portion of tissue from the right testis was lysed with RIPA lysis buffer to extract total proteins which were measured and separated by SDS–polyacrylamide gel electrophoresis (SDS–PAGE) and transferred to polyvinylidene difluoride (PVDF) membranes. Membranes were blocked in 5% skimmed milk for 1 h at room temperature and incubated with anti-PI3K, anti-p-Akt, anti-FOXO1, or anti-GAPDH antibodies overnight at 4°C. Membranes were washed three times with PBS, incubated with goat anti-rabbit secondary antibody for 1 h, and rewashed. Positive binding was visualized using an ECL developer, and the optical density of bands was analyzed using Quantity One software v.4.6.2 (Bio-Rad, Hercules, California, USA).
2.4.7. RT–qPCR Assay
100 mg testis tissues were homogenized and centrifuged, and 250 μL trichloromethane was added for 3 min before a second centrifugation. The supernatant was left at -20°C for 15 min and centrifuged at 4°C for 10 min. 20 μL RNase-free water was added to the precipitate to resuspend the RNA, and reverse transcription was performed by the cDNA reverse transcription kit according to the manufacturer’s instructions. PCR amplification was performed as follows: 35-45 cycles of predenaturation at 95°C for 2 min, denaturation at 95°C for 15 s, extension at 55-68°C for 30 s, and annealing at 45°C for 20 s. Primer sequences are shown in Table 1. Relative quantification (RQ) was performed by the method to determine differences in gene expression by comparing the ploidy of the target gene relative to the reference gene.
2.5. Statistical Analysis
All analyses were performed using SPSS software (version 26.0, Armonk, NY, USA). Continuous data are presented as . All data met the normal distribution and Chi-square requirement, and Student’s test was used to compare two groups. A value of was considered statistically significant.
3. Results
3.1. Target Screening
The top 500 depression targets and the top 500 MI targets were screened separately from the GeneCards database. The Venn diagram (Figure 2(a)) showed 81 common targets associated with depression and MI.

(a)

(b)

(c)

(d)
3.2. PPI Network Construction
Common targets for depression and MI were uploaded into the STRING database, and a network featuring the two diseases was constructed (Figure 2(b)) with a PPI network for visualization (Figure 2(c)) using Cytoscape software. The CytoHubba plug-in was used to analyze the PPI network and identify the top 10 targets, AKT1, VEGFA, IL6, TNF, IGF1, INS, TP53, ALB, PTPN11, and HRAS, according to their degree (Figure 2(d)).
3.3. GO and KEGG Enrichment Analyses
GO and KEGG enrichment analyses of the common targets were conducted via the DAVID database. The targets were involved in 410 BPs, 55 CCs, 54 MFs, and 138 signaling pathways. KEGG enrichment analysis revealed a possible association of the PI3K/Akt/FOXO1 signaling pathway with depression and MI. GO terms and KEGG pathways were visualized by the Omicshare website (https://www.omicshare.com/) according to their values (Figure 3).

(a)

(b)

(c)

(d)
3.4. Validation of a Rat Model of Depression
Model rats were sluggish, slow in movement, had an increased respiratory rate, reduced diet, and decreased body weight compared with controls. The SPT showed no significant differences in sugar-water consumption or sugar-water preference between the two groups on day 1 (Figure 4(a)). However, the sugar-water consumption and sucrose preference ratio were significantly lower in the model than in the control group by day 28 (). OFT scores for horizontal and vertical movements did not differ between the two groups on day 1, but scores were significantly lower for the model group than for controls by day 28 () (Figures 4(b) and 4(c)).

(a)

(b)

(c)
3.5. Organ Coefficients and Semen Quality
Testis and epididymis coefficients were reduced in model rats compared with controls (, ) (Figures 5(a) and 5(b)). Sperm concentration and motility were both lower in the model group () (Figures 5(c) and 5(d)), showing poorer semen quality in depressed rats than in controls.

(a)

(b)

(c)

(d)
3.6. HE Staining of Testis Tissues
Spermatozoa in the testes of the control group appeared normal in the structure under the light microscope with tightly arranged varicose seminiferous tubules and a large number of normal spermatozoa at all levels in the tubular lumen. Model spermatozoa showed an increased occurrence of abnormal structures of the seminiferous tubules. The arrangement of the seminiferous tubules was irregular, the lumen was reduced, and a large number of vacuolated structures were seen. In addition, an increase was seen in the proportion of abnormal spermatozoa and spermatocytes in tubules. Spermatogenic cells at all levels in the seminiferous tubules were reduced, the number of layers was disordered, and many spermatogenic cells were shed (Figure 6).

3.7. Protein Expression of PI3K, p-Akt, and FOXO1
By analyzing the gray value of the expression, the PI3K protein level in the control group was significantly higher than that in the model group (). Compared with the control group, the expression level of p-Akt in the model group was significantly reduced (). Compared with the control group, the expression levels of FOXO1 in the model group were significantly increased (). The results showed that the PI3K/Akt/FoxO1 pathway proteins were significantly changed in the samples from the rat depression model group (Figure 7).

(a)

(b)
3.8. Expression of PI3K, Akt, and FOXO1 mRNA
Compared to the control group, the expression level of PI3K mRNA in the model group was significantly reduced (). Compared to the control group, the expression level of Akt mRNA in the model group was significantly reduced (). The expression level of FOXO1 mRNA in the model group was significantly increased () (Figure 8).

4. Discussion
Infertile patients often feel low self-esteem, mistrust, lack of self-confidence, insecurity, and lack of faith [25, 26] with 25% to 60% reporting psychiatric symptoms and having higher depression levels than fertile couples [27]. Furthermore, depression tended to be more severe in patients with primary infertility than in those with secondary infertility [28]. Depression has a measurable impact on MI. Men with higher depression scores had significantly lower sperm concentrations and total sperm counts than nondepressed men [29]. Depression and MI thus form a vicious cycle, and disease progression is exacerbated.
The phosphatidylinositol 3 kinase (PI3K)/protein kinase B (Akt) signaling pathway is widespread and regulates such as proliferation, differentiation, apoptosis, and glucose transport [30]. PI3K/Akt has been closely associated with depression and is the target of antidepressants [31, 32]. PI3K/Akt signaling was considered to have the potential to improve depressive symptoms by regulating downstream targets, such as Bcl-2, mTOR, and FOXO1 [33, 34]. Akt knockout mice exhibited a depression-like phenotype [35], and levels of phosphorylated Akt in the ventral tegmental area were reduced in a mouse model of stress-induced depression, and fluoxetine reversed this effect [36]. The PI3k/Akt pathway inhibited glycogen synthase kinase-3β (GSK3β) [37] to enhance the synaptic transmission and antidepressant effects of ketamine [38]. FOXO1 is a member of the FOX family of transcription factors and is involved in antistress, apoptosis, cell proliferation, differentiation, and metabolism [39]. FOXO proteins were widely distributed in brain regions associated with emotional and stress regulation [40]. FOXO1 was mainly found in the striatum and hippocampus, and attenuation of its expression increased the differentiation of adult hippocampal neuronal progenitors [40]. Brain-derived neurotrophic factor (BDNF) is thought to be involved in the pathogenesis of mental disorders and regulates FOXO protein [41]. BDNF reduced FOXO transcriptional activity by stimulating tyrosine receptor kinase B (TrkB) and downstream kinases [42]. Activation of PI3K/AKT signaling promoted FOXO1 phosphorylation and inactivation with antidepressant effects [43]. Indeed, the FOXO pathway is the target for several antidepressant drugs, such as fluoxetine and venlafaxine, which regulate FOXO protein [44].
The PI3K/Akt pathway also regulates the spermatogenic function and promotes sperm proliferation and inhibits apoptosis. Activation of the downstream target proteins, mTOR, and p70S6K improved cyclin D3 expression in sperm cells, accelerating cell proliferation and differentiation [45]. FSH expression in supporting cells promoted PI3K/Akt activation and participated in mammalian spermatocyte proliferation, and PI3K inhibitors reduced supporting cell proliferation [46]. Moreover, cyclophosphamide inhibited the PI3K/Akt pathway to regulate FOXO1 expression and inhibit germ cell apoptosis [47]. Therefore, abnormal expression of the PI3K/Akt/FOXO1 pathway has been closely associated with depression and MI.
The current bioinformatic analysis and literature review have indicated the involvement of the PI3K/Akt/FOXO1 signaling pathway in the link between depression and MI. Organ coefficients and semen parameters were significantly decreased in depressed rats, and testicular spermatogenesis was impaired, indicating that depression had a negative effect on MI. PI3K and p-Akt protein and mRNA were decreased, while FOXO1 was increased, consistent with the findings of other scholars. Therefore, abnormal expression of the PI3K/Akt/FOXO1 signaling pathway has been implicated in the reduced spermatogenic function in the rat model of depression.
A potential mechanism to link depression and MI is presented as a basis for future experimental and clinical studies. However, we acknowledge some limitations of the present study. Many depression-MI-related pathways were predicted, and the PI3K/Akt/FOXO1 pathway was chosen for further study. Previous work has suggested this pathway to be involved with the two diseases. Future experiments concerning the effects of depression on reproductive hormones and the roles of apoptosis and autophagy on depression and MI are planned. The impact of depression on spermatogenesis in rats is in the experimental exploration stage, and refinements to experimental design and drug intervention investigations are required.
5. Conclusion
Bioinformatics analysis showed that depression and MI were correlated with the PI3K/Akt/FOXO1 signaling pathway. Spermatogenic function was decreased in a rat model of depression. Protein and mRNA levels of PI3K and p-Akt were reduced, and FOXO1 was increased. Depression may thus reduce spermatogenic function in rats, perhaps through a mechanism involving differential expression of PI3K/Akt/FOXO1 pathway proteins.
Data Availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.
Conflicts of Interest
The authors declare that there is no conflict of interest.
Authors’ Contributions
Qi Zhao, Fei Yan, and Chuan Ma contributed in conception and designed of the manuscript. Junlong Feng and Sheng Deng contributed in animal experiments and data collection. Cong Zhao and Huanan Zhang contributed in data analysis and manuscript modification. Haisong Li, Bin Wang, and Jisheng Wang oversaw the study. All authors read and approved the final manuscript. Qi Zhao, Fei Yan, and Chuan Ma contributed equally to this work and share first authorship.
Acknowledgments
This work was funded by the National Natural Science Foundation of China (no. 81973647) and the Capital Health Development Research Special Project (no. 2020-2-4191). The author sincerely thanked Beijing University of Chinese Medicine for the assistance and site support provided during this period.
Supplementary Materials
Figure S1: uncropped western blots. (Supplementary Materials)