1. Introduction
Mycotoxins are chemicals produced by fungi that have a highly toxic impact on human and animal health even at low concentrations [
1]. They contaminate at least 25% of the global food supply [
2], and that percentage is expected to increase with climate change [
3,
4]. The most common toxigenic fungal species are in the genera Aspergillus, Penicillium, and Fusarium [
5]. The population structure of these microorganisms in foods depends on factors such as temperature, gas composition and water availability [
5,
6]. Fumonisins are mycotoxins produced by Fusarium verticillioides and Fusarium proliferatum that contaminate some plants and their food products, mainly corn crops, and many other plants, such as sorghum, cowpea, asparagus, rice and wheat [
7]. Thus, regions where corn crops are used as dietary staples are highly exposed to fumonisin B1 (FB1), B2 (FB2) and B3 (FB3) [
8,
9].
It has been described that FB1 exposure is associated with esophageal cancer [
10,
11], porcine pulmonary edema [
12,
13], liver and kidney toxicity [
14,
15], and liver and kidney cancer [
16,
17] in different animal species. Additionally, pregnant animals fed FB1-contaminated feed have growth retardation, delayed or incomplete ossifications, cleft palate or hydrocephalus and fetal death [
18]. Neural tube defects (NTDs) due to fumonisin exposure, including spina bifida, exencephaly and craniorachischisis or meningomyelocele, have been reported [
6,
18].
To date, very few studies have focused on the effect of FB1 on female reproduction. These studies suggested the role played by FB1 in affecting female fertility, as dietary concentrations of FB1 reduced serum gonadotropin levels and lowered fertility without inducing histopathological changes in the ovaries [
19]. Furthermore, it has been reported that FB1 exposure leads to granulosa cell proliferation and affects progesterone production induced by FSH plus IGF-I [
20]. However, there is still a need to explore the interactive effects of FB1 on ovarian function and female fertility in detail, as the exact underlying mechanism is unfamiliar. In this study, we provide comprehensive knowledge of the intergenerational effect of FB1 on the structure and function of the ovary. In particular, histopathological changes, follicular development and folliculogenesis and steroidogenesis-related genes are investigated. We also illustrate the role of autophagy as a recovery mechanism in the ovaries of the second generation to protect against the prenatal effect of FB1.
3. Discussion
Mycotoxins are produced by more than one hundred fungi that belong to the genera Aspergillus, Penicillium, and Fusarium. They produce secondary metabolites with low molecular weights that contaminate various agricultural commodities [
23,
24], especially rice and corn, which are widely used ingredients in food in developing countries [
25,
26]. Their occurrence in the food chain is considered a serious problem worldwide [
27,
28]. In particular, fumonisin B1, which was first isolated from Fusarium [
29], mainly F. moniliforme, has been reported to be carcinogenic to humans [
30,
31] and may have a direct or indirect effect on the increased prevalence of some common diseases, such as female infertility, as notable numbers of women are undergoing fertility treatment [
32]. Consequently, it is vital to investigate the interactive effects of these mycotoxins on ovarian function.
The characteristics and quality of the ovary greatly affect its function in reproduction; therefore, any disruption to ovarian function can be associated with an ovarian disorder, such as polycystic ovary syndrome [
33]. Our results showed that the number of primordial follicles was significantly decreased in F1 and F2 females treated with 50 mg/kg fumonisin B1. It has been shown that the formation of primordial follicles can be influenced by several genes that control oocyte survival or apoptosis, thereby affecting the number of follicles formed [
34,
35]. Thus, FB1 affects the ovarian pool of healthy primordial follicles during early fetal life and is associated with a faster decline in ovarian function with aging. This finding is in agreement with previous studies reporting that prenatal exposure to environmentally relevant substances accelerates the age-related decline in reproductive capacity in the F1 generation [
22,
36]. Furthermore, the significant decrease in the number of primordial follicles in the 50 mg/kg treatment group of F2 females was associated with an increased number of tertiary follicles, suggesting that FB1 might cause early menarche by inducing early folliculogenesis. In fact, as the offspring of the second generation have a lower ovarian reserve of primordial follicles and therefore a shorter reproductive lifespan than those of the control, folliculogenesis was promoted early to maintain high fertility similar to the control, leading to a successful reproductive lifespan. This hypothesis seems consistent with previous studies [
21,
37,
38]. Thus, mothers of the first generation pass on a new phenotype that is better suited for the effect of FB1.
The high dose of FB1 significantly affected ovary weight among the two generations compared to the control. In fact, there was a significant increase in ovary weight at the higher doses for both the first and second generations compared to the control group. This finding is in agreement with previous reports showing that treatment with some substances, such as acrylamide, sodium fluoride, BPA and nonylphenol, led to an increase in ovary weight [
39,
40,
41,
42]. Although cyclic fluctuations are connected to ovarian function during the estrous cycle, it is known that ovarian weight does not significantly change in normal rats. Consequently, any change in ovarian weight, whether positive or negative, should be viewed as a sign of ovarian dysfunction that may be caused by a variety of conditions, including persistent polycystic ovaries, oocyte and follicle depletion, luteal cyst development, and reproductive aging. In particular, the dominant luteinized follicle shrinks and ruptures in a normal cycle, but in some circumstances, its cystic nature can persist during the luteal phase. This condition is known as "luteinized unruptured follicle syndrome (LUFS)", which is an abnormal condition in which the follicle does not rupture during the luteal phase, causing cysts and leading to an increase in ovarian weight and a decrease in fertility [
43,
44]. Luteal cyst development has been described in spontaneously cycling women with unexplained infertility [
45]. Thus, we suggest that high doses of FB1 may lead to the development of LUFS and probably to polycystic ovary syndrome (PCOS), the most common endocrine disorder in premenopausal women [
46]. This finding is concomitant with the results of the steroidogenesis- and folliculogenesis-related gene expression in both treatment groups of the first generation. In fact, PCOS is thought to be caused by many intraovarian disturbances in steroidogenesis, including deficiency in the activity of CYP19, the enzyme that catalyzes the rate-limiting step in the biosynthesis of estrogens from androgens [
47]. Thus, a decrease in the activity of this enzyme could be expected to result in a decrease in the expression of estrogen receptor beta (ESR2) and therefore increased ovarian androgen production and the development of PCOS [
48,
49]. In light of these facts, follicle growth and development of the first generation have been affected, leading to a decrease in the number of growing follicles, notably secondary and tertiary follicles. This may be explained by the degenerative follicles due to the vacuolized oocytes and disorganized granulosa cells with pyknotic nuclei, similar to the results reported by many previous studies [
22,
50,
51]. As a result, the fertility rate of the first generation significantly decreased in both the 20 and 50 mg/kg treatment groups.
However, the results in the ovaries from the treatment groups of the second generation were different and unexpected. Indeed, while ovary weight significantly increased with the high dose of treatment, the steroidogenesis- and folliculogenesis-related genes did not decrease as they had among the females of the first generation. While there was no significant variation in CYP19, ESR2 significantly increased. This result led us to think about the eventual role played by the autophagy process, as LC3 significantly increased in both groups of the second generation but significantly decreased in both groups of the first generation. The promoted autophagy among the second-generation ovaries may be considered a recovery mechanism from the effect of FB1. This finding is consistent with the stereological and histological results, as no significant difference was found in the fertility rate of the second generation compared to that of the control. Interestingly, the number of primary and secondary follicles of the treatment group was similar to that of the control group, and many vacuoles were observed within oocytes and between the layers of granulosa cells. These vacuoles could be interpreted as autophagosome structures related to autophagy. In fact, proper functional autophagy is required for the normal growth and development of follicles [
52], and its occurrence is an indicator of an adaptation to stress and can lead to reduced apoptotic cytotoxicity [
53,
54]. Its key role in maintaining normal cell homeostasis and preventing chronic cellular damage by removing toxins, damaged organelles, and unfolded proteins has been widely described [
55,
56]. Many previous studies have reported the association between autophagy and different ovarian disorders, including follicular cyst formation, metabolic abnormalities and PCOS [
57]. In particular, defective autophagy in ovarian cells leads to poor-quality oocytes, resulting in female infertility [
58].
4. Materials and Methods
4.1. Ethical Statement
This study was approved by the Scientific Research Ethics Committee at King Saud University, Riyadh, Saudi Arabia (Reference No: KSU-SE-22-41) and carried out in accordance with the approved guidelines. All experimental procedures are reported in accordance and compliance with the ARRIVE guidelines.
4.2. Study Design and Sampling
Thirty healthy pubertal virgin female Wistar-Albino rats (weighing 200–250 g) were housed separately in cages and kept in a facility with a standard laboratory chow diet and a 12–12 photoperiod at a temperature of 21–1°C. The rats were then kept with males, and when a white vaginal plug appeared on the cage flooring, mating was considered successful; that day was regarded as day 0 of gestation (GD 0). Thereafter, pregnant females were separated into three groups, and the following treatment regimen was followed from GD 6 to GD 21: 1) The first group of females (n = 10) received a gavage of distilled water and was designated the control group; FB1 was administered orally to the second group of females (n = 10) at a dose of 20 mg/kg; and FB1 was administered orally to the third group of females (n = 10) at a dose of 50 mg/kg. Following parturition, we obtained the first generation of offspring from mothers who had been treated with FB1, which are referred to as the animals of the first generation (F1); the control group is referred to as the control group of the first generation (CF1). When the female F1 and CF1 offspring were 4 weeks old (before puberty), a fraction of them was moved into a clear plastic box with a carbon dioxide tube attached and a flow rate of 10 L/h for 10 minutes in preparation for euthanasia. The ovaries were cleaned, rapidly measured, and assigned to groups based on where they came from. The remaining females from F1 and CF1 were allowed to mature sexually and were mated with males. The second generation of offspring was obtained from FB1-treated F1 mothers. These females are referred to as the FB1 offspring of the second generation (F2), while those obtained from the control group are referred to as the control group of the second generation (CF2). When the female F2 and CF2 offspring were 4 weeks old, they were euthanized, and their ovaries were sampled.
4.3. Histological Preparation
Neutral buffered formalin (NBF 10%) was used to fix ovary samples for a period of 24 h. The next day, they were embedded in paraffin and cut into sections with a thickness of 5–7 μm, collected on a hotplate and transferred to glass slides containing warm (30°C) water and albumin glycerol fixative for adhesion. Wrinkles were removed, and the sections were stained with hematoxylin and eosin for histological study [
21].
4.4. Immunofluorescence Staining and Confocal Microscopy
Immunostaining was performed as described in previous studies [
22]. We placed slides containing tissue sections on a hot plate (60 °C) and deparaffinized them with xylene. They were then rehydrated and washed twice with distilled water and three times with 1x PBS. The slides were removed from the wash solution and dried. After drying, the sections were placed in a suitable container, permeabilized with 0.1% Triton X-100 containing 0.1% sodium citrate, and treated with blocking buffer (1% BSA in PBS) at room temperature. Slides were placed in a humid box and incubated with primary antibody solution (anti-LC3) (dilution 1:500) from DGpeptidesCo. Ltd. overnight at 4°C on a flat surface in the dark. The next day, the slides were washed four times with 1x PBS and treated with the secondary antibody FITC (dilution 1:700, Abcam, USA) for 45 min at room temperature (RT) in the dark. Then, the slides were washed with PBS and TE buffer before adding Hoechst solution (diluted 1:15000, Hoechst 33342, Life Technologies, USA). Finally, the sections were placed in 50% glycerol/TE solution, and the edges were sealed with nail polish. Sections were observed and imaged for signal quantification with a spinning disk confocal microscope from Zeiss. The signal intensity for protein expression was analyzed by Zen 3.1 service (ZEN lite) and quantified using the GraphPad Prism 9 program (GraphPad Software).
4.5. Analysis of Gene Expression
RNA was extracted using the RNeasy Mini Kit (Qiagen, Westburg, The Netherlands) with DNase treatment on columns using an RNase-free DNase kit (Qiagen). Using a NanoDrop with a 260/280 nm ratio, we measured the quality and purity of the extract. Using RT–PCR and primer sets using an iScript ™ cDNA synthesis kit (Applied Biosystems, Carlsbad, CA) according to the manufacturer’s instructions, cDNA was reverse-transcribed from 0.1 to 0.5 µg of total RNA. Finally, real-time PCR (RT–PCR) was conducted by using SYBR green and an Applied Biosystems 7500 Fast RT–PCR system (Carlsbad, CA) with the gene-specific primers shown in Table 1. We calculated the relative amount for each gene transcript using the 2-ΔΔCT method and normalized by referencing to the gene GAPDH.
4.6. Statistical Analysis
The data were analyzed using GraphPad Prism version 9. One-way analysis of variance, followed by Tukey’s multiple comparison, was used for statistical comparisons. All values are presented as the mean ± standard deviation (SD). Significance was set at a P value < 0.05.