Skip to main content


Springer Nature is making SARS-CoV-2 and COVID-19 research free. View research | View latest news | Sign up for updates

Biologically active substances-enriched diet regulates gonadotrope cell activation pathway in liver of adult and old rats


According to the Hippocrates’ theorem “Let food be your medicine and medicine be your food”, dietary interventions may induce changes in the metabolic and inflammatory state by modulating the expression of important genes involved in the chronic disorders. The aim of the present study was to evaluate the influence of long-term (14 months) use of biologically active substances-enriched diet (BASE-diet) on transcriptomic profile of rats’ liver. The experiment was conducted on 36 Sprague–Dawley rats divided into two experimental groups (fed with control or BASE-diet, both n = 18). Control diet was a semi-synthetic diet formulated according to the nutritional requirements for laboratory animals. The BASE-diet was enriched with a mixture of polyphenolic compounds, β-carotene, probiotics, and n-3 and n-6 polyunsaturated fatty acids. In total, n = 3,017 differentially expressed (DE) genes were identified, including n = 218 DE genes between control and BASE groups after 3 months of feeding and n = 1,262 after 14 months. BASE-diet influenced the expression of genes involved particularly in the gonadotrope cell activation pathway and guanylate cyclase pathway, as well as in mast cell activation, gap junction regulation, melanogenesis and apoptosis. Especially genes involved in regulation of GnRH were strongly affected by BASE-diet. This effect was stronger with the age of animals and the length of diet use. It may suggest a link between the diet, reproductive system function and aging.


The influence of diet on human health has been studied for a long time. In the fourth century BC, Hippocrates already stated: “Let food be your medicine and medicine be your food”. This theory was again recalled in the nineteenth century thanks to Ludwig Feuerbach who wrote “Man Is What He Eats”. Today’s scientific knowledge confirms the words of these great scientists. Recent studies have shown that our diet is an important factor which could contribute not only to the development but also to the inhibition of chronic diseases, including osteoporosis, diabetes, cancer, atherosclerosis, cardiovascular disease, neurodegenerative diseases and obesity (Virmani et al. 2006, 2013; Joseph et al. 2009; Lillycrop and Burdge 2012). It is possible due to nutritional factors, which may induce epigenetic changes also via direct influence on gene expression (Alam et al. 2012). In the late twentieth century, Nancy Fogg-Johnson and Alex Meroli created a new term which combines two fields of science: nutrition and genetics, named nutrigenomics. The aim of nutrigenomics is to study how various food ingredients affect the expression of specific genes and therefore provide tools to understand and control the worldwide epidemic of specific chronic diseases. It has been proven that these diseases more often arise from dysfunctional biological networks, instead of single common gene mutation (Liu et al. 2010; Ferguson et al. 2007a, b; Astley 2007).

Dietary interventions may induce changes in the metabolic and inflammatory state by modulating the expression of important genes involved in the chronic disorders (Lottenberg et al. 2012). This particularly applies to biologically active substances which have a protective influence on the body. For example, recent studies have shown that polyunsaturated fatty acids (PUFA) downregulate ATP-binding cassette transporter A-1 (ABCA-1), which results in reductions in HDL-C concentrations. PUFA suppress the liver X receptor/retinoid X receptor (LXR/RXR) gene responsible for ABCA-1 synthesis, which is a transporter involved in the HDL formation (Uehara et al. 2007). Furthermore, PUFA modulate expression of several genes involved in oxidative processes (such as PPAR-α), while impairing the sterol-regulatory element binding proteins (SREBPs) involved in lipogenesis (Hannah et al. 2001). A different study demonstrated that diet enriched with anti-inflammatory mixture containing resveratrol, green tea extract, α-tocopherol, vitamin C, omega-3 polyunsaturated fatty acids and tomato extract affected genes involved in inflammatory processes, oxidative stress and metabolism (Bakker et al. 2010). Also various dietary components including omega-3 fatty acids, plant flavonoids and carotenoids have been demonstrated to modulate gene expression, by decreasing iNOS and COX-2 gene expression induced by gliadin in RAW 264.7 macrophages stimulated with IFN-γ. Therefore, these compounds could preserve intestinal barrier integrity, play a protective role against toxicity of gliadin peptides and have a role in nutritional therapy of, for example, celiac disease (Ferretti et al. 2012; De Stefano et al. 2007).

Taking the latest data from the fields of nutrigenomics into consideration, the aim of the present study was to evaluate the influence of long-term (14 months) use of biologically active substances-enriched diet (BASE-diet) containing the mixture of polyphenolic compounds, β-carotene, probiotics and n-3 and n-6 polyunsaturated fatty acids on transcriptomic profile of rats liver. However, in the present study we did not want to investigate the mechanisms of action of the individual components, as these are already known. The main emphasis has been placed on the possible cumulative action of these compounds in the situation when they were added simultaneously as ingredients of semi-synthetic diet. Since 14-month period in case of rats is more than half of their life’s span, the obtained results could be of great value.

Materials and methods

Animal and diets

Experiment was carried out on thirty-six 8-week-old male Sprague–Dawley rats (Charles River Laboratories, Germany). The experiment was performed with the approval of the local ethical committee. Rats were kept in individual cages. Water and food were available ad libitum. The environment was regulated at 22 ± 0.5 °C, air humidity of 50 % on a 12 h/12 h L/D photoperiod throughout the entire experiment. The body weight and the feed consumption were recorded weekly. Animals were divided into two groups: a control group (n = 18) and an experimental group (n = 18).

Rats from both groups were receiving for 14 months semi-synthetic diets (control and BASE-diet) formulated according to the nutritional requirements for laboratory animals (NRC 1995) and are described in Table 1. Application of the semi-synthetic diet allowed the elimination of the additional impact of biologically active compounds contained in commercially available diets. BASE-diet was additionally enriched with the following biologically active compounds: 6 % of salmon fat replacing lard (to increase the level of unsaturated fatty acids), 8 % of hydrolysed water extract from small-leaved linden (Tilia cordata) inflorescence (as a source of antioxidant compounds), 8 % of puree from giant pumpkin (Curcubita macima) (the source of beta-carotene) and 1 % of two strains of bacteria with documented probiotic activity: Lactobacillus acidophilus LA-5 and Bifidobacterium animals ssp. lactis. The dry matter content of both diets was the same. Exact composition of macronutrients, micronutrients and vitamins mixtures included in the experimental diets is presented in Table 2.

Table 1 Composition of diets used for 14-month-long feeding of rats
Table 2 Composition of the mixtures of macronutrients, micronutrients and vitamins used in control and BASE-diet fed to rats for 14 months

The content of biologically active substances in feed was measured as follows. The fatty acid content was determined both in the raw ingredients and in the feed samples using gas chromatography (Table 3). Quantitative analysis of the total polyphenolic compounds in the hydrolysed water extract from small-leaved linden inflorescence was made using the method with Folin reagent (1 N) in the presence of Na2CO3 (20 %). The analysis of beta-carotene and other carotenoids content was performed using high-performance liquid chromatography coupled with electrochemical detection (HPLC–ECD).

Table 3 The fatty acid profiles in control and BASE-diets

After 3 and 14 months of experiment, animals from each experimental group (n = 9) were euthanized by exsanguination under general anaesthesia with the isoflurane. Liver samples from six randomly selected rats from each group were taken.


Liver samples were frozen with liquid nitrogen immediately after collection and stored at −80 °C until extraction. Total RNA was extracted with RNeasy Lipid Tissue Mini Kit (Qiagen, Germany) following the manufacturer’s recommended protocol. Subsequently, potential genomic DNA contamination was eliminated using Deoxyribonuclease I Amplification Grade (Sigma, USA) and RNeasyMinElute Cleanup Kit (Qiagen, Germany). RNA quantity and quality was measured using NanoDrop 2000 (NanoDrop Technologies, USA) and Bioanalyzer (Agilent Technologies, USA). To ensure optimal data quality, only RNA samples with RIN number ≥8.8 were included in the analysis.

The analysis of gene expression profile was performed using SurePrint G3 Rat Gene Expression Microarray, 8 × 60 K (Agilent Technologies, USA). The Low Input Quick Amp Labeling Kits (Agilent, USA) was used to amplify and label target RNA to generate complementary RNA (cRNA) for oligo microarrays used in gene expression profiling. Experiment was performed using a common reference design, where the common reference was RNA from 10-week-old healthy rats not participating in the experiment, housed for 2 weeks in the same room with experimental rats. On each two-colour microarray, 300 ng of cRNA from the control rats (labelled by Cy3) and 300 ng of cRNA from the rats fed with BASE-diet (labelled by Cy5) were hybridized. Microarray hybridization was performed with the Gene Expression Hybridization Kit (Agilent Technologies, USA), according to the manufacturer’s protocols. RNA Spike In Kit (Agilent Technologies, USA) was used as an internal control. In total, 24 microarrays, one for each animal, were done. Acquisition and analysis of hybridization intensities were performed using an Agilent DNA microarray scanner. The final analysis was carried out on 22 microarrays, which passed through the Agilent Feature Extraction’s and the Gene Spring’s control (one microarrays from the control group after 3 months of experiment and one from the experimental group after 3 months of experiment were rejected due to insufficient quality).

Signal detection and statistical analysis

Data were extracted and background was subtracted using the standard procedures contained in the Agilent Feature Extraction (FE) Software version FE performs also Lowess normalization. The statistical analysis was performed using Gene Spring 12 software (Agilent, USA). The samples underwent quality control and the results showed that each sample had similar QC metric profile. The next step was filtering probe sets by flags to remove poor-quality probes (absent flags). The statistical significance of the differences was evaluated using one-way ANOVA and Tukey’s HSD post hoc test (p < 0.05). A multiple testing correction was performed using Benjamini and Hochberg False Discovery Rate (FDR) <5 %. Microarray data were deposited at the Gene Expression Omnibus data repository under the number GSE51657 according to the MIAME requirements.

To identify the list of signalling pathways, the microarray data were analysed using Pathway Studio 6.0 (Ariadne Genomics).

Real-time PCR

To verify microarray results, the expressions of three randomly selected genes (GPX1, IRF7 and PRODH) were measured using real-time PCR method. The sequences of these genes were obtained from Ensembl database. Primers were designed using Primer-Blast software (NCBI database) and then checked for secondary structures using Oligo Calculator (free on-line access). The secondary structures of the amplicon were examined using m-fold Web Server (free on-line access). As a house keeping gene ACTB were used (Huang et al. 2013; Wang et al. 2013). The sequences of the primers are listed in Table 4.

Table 4 Primer sequences for real-time PCR verification of microarray results

cDNA was synthesized using Enhanced Avian HS RT-PCR Kit (Sigma-Aldrich, St. Louis, Missouri). All analyses were performed on individual samples of total RNA using a Brilliant III Ultra-Fast SYBR Green QPCR Master Mix kit (Agilent Technologies, USA) following the manufacturer’s protocol. Each sample was tested 3 times in a Stratagene Mx3005P Quantitative PCR instrument for RT-PCR. The relative expression of the target gene was calculated according to the following formula:

$$\Delta \varDelta {\text{CT}} = \Delta {\text{CT}}\left( {\text{sample}} \right) - \Delta {\text{CT}}\left( {\text{normal}} \right)$$

where ΔCT is the difference in CT between the targeted gene and housekeeping controls by minimizing the average CT of the controls. The fold-change calculated as: 2−ΔΔCT (Livak and Schmittgen 2001).


There was no difference in initial body weight and liver to body weight ratio among all experimental groups. Over the course of the 14-month study, all animals gained body weight in a time-dependent manner, regardless the treatment (Figs. 1, 2). There was also no difference in the amount of diet eaten between control and experimental rats (Fig. 3). Earlier published data showed that there were no statistically significant differences between control and BASE-diet groups in blood morphology parameters as well as biochemical liver function parameters (Oszkiel et al. 2014).

Fig. 1

Weight of rats participating in the experiment. CD 3—rats from control group after 3 months of experiment; BD 3—rats from BASE group after 3 months of experiment; CD 14—rats from control group after 14 months of experiment; BD 14—rats from BASE group after 14 months of experiment

Fig. 2

Liver to body weight ratio (%). CD 3—rats from control group after 3 months of experiment; BD 3—rats from BASE group after 3 months of experiment; CD 14—rats from control group after 14 months of experiment; BD 14—rats from BASE group after 14 months of experiment

Fig. 3

Average food consumption [g] of rats participating in the experiment. CD—rats from control group; BD—rats from BASE group

The consecutive steps of microarray data analysis have been presented in Fig. 4. Analysis started from identification of all differentially expressed (DE) genes in whole experiment and was conducted using ANOVA in Gene Spring software. In total, n = 3,017 DE genes were identified. Then, using Tukey’s HSD post hoc test the lists of differentially expressed genes between experimental groups were identified (Table 5). This analysis revealed n = 218 differentially expressed genes between control and BASE groups after 3 months of feeding, and n = 1,262 differentially expressed genes between control and BASE after 14 months of feeding. These two sets of differentially expressed genes were considered as diet-induced genes since the differences in expression resulted only from the diet fed to the animals.

Fig. 4

The analysis diagram of the data obtained in the experiment

Table 5 The number of differentially regulated genes between each experimental group (total number of differentially regulated genes n = 3,017)

Each of the two data sets was analysed using Pathway Studio 6.0 (Ariadne Genomics) software in order to identify statistically significant Ariadne Cell Signalling Pathways in which differentially expressed genes were involved (based on Gene Ontology). After 3 months of BASE-diet feeding, three Ariadne Cell Signalling Pathways were significantly regulated (p value <0.05) (Table 6). After 14 months of BASE-diet feeding, five Ariadne Cell Signalling Pathways were significantly regulated (p value <0.05) (Table 7). There were two signalling pathways which were regulated both after 3 and 14 months of BASE-diet administration, namely guanylate cyclase pathway and gonadotrope cell activation pathway.

Table 6 The list of Ariadne Cell Signalling Pathways significantly regulated in rats liver after 3 months of BASE-diet feeding
Table 7 The list of Ariadne Cell Signalling Pathways significantly regulated in rats liver after 14 months of BASE-diet feeding

Obtained microarray results were also validated using real-time PCR on randomly selected three genes. As shown in Fig. 5, gene expression for GPX1, IRF7, PRODH matched the expression obtained from microarray analysis.

Fig. 5

Expression of GPX1, IRF1 and PRODH genes in livers of control rats and rats fed with BASE-diet measured using DNA microarrays and real-time PCR (normalized vs. ACTB)


In the present study, we investigated the influence of biologically active substances-enriched diet (BASE-diet) fed for a long period of time (14 months) on the transcriptome of rat liver. It means that two factors could influence the results—the experimental diet and age of rats. Based on the numbers of genes regulated in individual comparisons (adult vs. older; control diet vs. BASE-diet), we could conclude that age factor was responsible for the regulation of higher number of genes (number of differentially regulated genes between adult and old BASE-diet fed rats was n = 2,255) than the diet factor (number of differentially regulated genes between control and BASE-diet fed rats was n = 218 after 3 months of experiment and n = 1,262 after 14 months of experiment) (see Table 5).

Our previous study has shown that BASE-diet, which was enriched with the mixture of polyphenols, beta-carotene, probiotics and polyunsaturated fatty acid, prevented from hepatic and systemic oxidative damage and was able to attenuate the development of some senile features in adult and old rats (Oszkiel et al. 2014). In rats fed with BASE-diet, the GSH/GSSG ratio has increased, GSH-Px and GSSG-R activities decreased, and SOD activity decreased when compared to control animals. These results suggest that BASE-diet positively influenced some parameters of antioxidant defence within the body. Because of that we investigated whether differentially regulated genes were involved in antioxidant activity. Analyses performed in Pathway Studio Software show relations between proteins encoded by genes involved in antioxidant activity and proteins encoded by genes differentially regulated in livers of rats fed with BASE-diet (Figs. 6, 7). It is clearly visible that in older rats the number of differentially regulated genes involved in antioxidant activity is significantly higher when compared to adult rats.

Fig. 6

Interactions between proteins involved in antioxidant activity (blue) and proteins encoded by differentially regulated genes (red) in livers of rats after 3 months of BASE-diet feeding

Fig. 7

Interactions between proteins involved in antioxidant activity (blue) and proteins encoded by differentially regulated genes (red) in livers of rats after 14 months of BASE-diet feeding

As expected, we were able to identify differentially expressed genes between control and BASE-diet both in adult (after 3 months of feeding) and older rats (after 14 months of feeding). Significant difference in the expression of these genes resulted in significant regulation of a few signalling pathways (three in the case of adult rats and five in the case of older rats). Two of these pathways were the same in adult and older rats, so we concluded that these two were purely dependent on the diet.

BASE-diet significantly influenced the expression of gonadotrope cell activation pathway and guanylate cyclase pathway. Gonadotrope cell activation pathway is involved in regulation of the reproduction system. It regulates the activity of the basophilic cells of the anterior pituitary gland specialized in secreting follicle-stimulating hormone (FSH) or luteinizing hormone (LH). FSH and LH are released by a trophic peptide hormone—gonadotropin-releasing hormone (GnRH). FSH stimulates development of ovarian follicles and regulates spermatogenesis in the testis. LH causes ovulation and formation of the corpus luteum in the ovary, stimulates production of oestrogen and progesterone by the ovary and stimulates testosterone production by the testis (Guyton and Hall 2006; Okamura et al. 2013). Previous studies reported that reproductive activity can be directly influenced by nutrients, including the biologically active compounds. Among such compounds are polyphenols which were included into BASE-diet. Chen et al. (2010) indicate that genistein and resveratrol can increase the ovarian follicular reserve and prolong the ovarian lifespan in rats. Gebrie et al. (2005) proved that the extract of Rumex steudelii (Ethiopian plant which roots contain phytosterols and polyphenols) prolonged significantly the oestrus cycle and the dioestrous phase. Studies of Colitti et al. (2007) revealed that polyphenols can protect from stress which activates the hypothalamo–pituitary–adrenal axis leading to enhanced glucocorticoid secretion and concurrently disrupts ovarian cycle. Thereby, polyphenols can improve reproduction activity (Matteri et al. 2001; Zhang et al. 2006; Colitti and Stefanon 2006; Sgorlon et al. 2006; Colitti et al. 2007).

Gonadotrope cell activation pathway could also be regulated by β-carotene. β-carotene is not only a retinol precursor, but also fulfils similarly function to vitamin E. It is a free radical scavenger, acting especially on singlet oxygen, and thus it is a potent antioxidant, similarly to other carotenoids. Moreover, many gene products linked to reproduction can be modulated by the product of retinol oxidation—retinoic acid. Therefore, it is considered that an optimal intake of β-carotene has a positive effect on fertility (Schweigert et al. 2003; Arechiga et al. 1998). Recent study has shown a direct relationship between β-carotene concentration at ovarian level and size and progesterone secretion by corpora lutea. It is suggested that β-carotene may have a positive effect on luteogenesis and luteal activity (Haliloglu et al. 2002; Arellano-Rodriguez et al. 2009).

Also the fat content in the diet affects the reproductive system. Study carried out on cows shows that fat supplementation is associated with increased dominant follicle diameter, greater progesterone concentrations, modulation of prostaglandin synthesis and improved oocyte and embryo quality. It collectively results in increasing likelihood of conception (Garcia-Bojalil et al. 1998; Lucy et al. 1993; Mattos et al. 2002).

However, there is some evidence that diet affects the reproductive system also on the transcriptome level. Majority of the studies focusing on transcriptional regulation of GnRH genes concerned GnRH-I; however, recently also a study concerning GnHR-II has been published (Lee et al. 2008). In our study, GnRH-I and GnRH-II genes, which are expressed mostly in the brain, were not directly regulated. However, we identified 38 genes involved in indirect regulation of GnRH-I and GnRH-II expression (Figs. 8, 9). The expression of some genes involved in regulation of GnRH can be modulated by food interventions. Plant polyphenols can alter the expression levels of oestrous cycle genes encoding PGHS-2 (upregulation), SOD2 and FOXO3 (downregulation) (Colitti et al. 2007). Dietary β-carotene can influence reproductive system on the transcriptome level as well. β-carotene is a precursor (inactive form) of vitamin A. It is currently believed that vitamin A regulates gene transcription also by retinol-binding proteins (RBP). Expression of the RBP genes is dependent on progesterone, which is dominant ovarian steroid hormone, known to be involved in the regulation of gonadotropin secretion. Progesterone regulates the GnRH-I gene through a feedback mechanism (Van Arnum 1998; Berry et al. 2011; Costello et al. 2010; Lee et al. 2008). Moreover, it has been demonstrated that unsaturated fatty acids may influence reproduction directly interfering with basal and GnRH-dependent gonadotrope activity (Garrel et al. 2011). Diet enriched with fish oil rich in 20:5 and 22:6 n-3 fatty acids modulate the hepatic expression of genes influencing reproductive performance: SREBF1, ASCL1 and FABP1 (Hutchinson et al. 2012). Dietary n-3 fatty acids also increase the progesterone receptor (PR) mRNA and oestrogen receptor-alpha (ER-1) expression (Bilby et al. 2006). Additionally, free fatty acids might directly modulate pituitary gonadotropin production by upregulating Lhb mRNA expression and suppressing Fshb mRNA expression (Sharma et al. 2013). None of the above-mentioned genes have been directly regulated in our experiment; however, a significant part of differentially expressed genes in the livers of BASE-diet fed rats was involved in regulation of all above-mentioned food-interventions-related proteins (Figs. 10, 11).

Fig. 8

Interactions between GnRH-I and GnRH-II (blue) and proteins encoded by differentially regulated genes (red) in livers of rats after 3 months of BASE-diet feeding

Fig. 9

Interactions between GnRH-I and GnRH-II (blue) and proteins encoded by differentially regulated genes (red) in livers of rats after 14 months of BASE-diet feeding

Fig. 10

Interactions between “food interventions-related proteins” involved in reproductive system regulation (green) and proteins encoded by differentially regulated genes (red) in livers of rats after 3 months of BASE-diet feeding

Fig. 11

Interactions between “food interventions-related proteins” involved in reproductive system regulation (green) and proteins encoded by differentially regulated genes (red) in livers of rats after 14 months of BASE-diet feeding

However, there are some latest data indicating that gonadotrope cell activation pathway gene could be expressed not only in reproductive tissues (Millar 2005). So far, the only well-known pathway connected with GnRH action was associated with fertility regulation in the hypothalamic–pituitary–gonadal axis. However, current research showed also that expression of genes involved in regulation of gonadotrope cell activation pathway was observed in other cells such as embryonic stem cell, adult brain, skeletal muscle, myocardium tissues, breast cancer stem cells, prostate cancer cell lines, pancreatic ductal adenocarcinoma, liver, kidney, thyroid cancer cells, mammary gland, prostate cancer, bone matrix and neuroblastoma cells (Ames et al. 2013; Chaudhry et al. 2013; Datta et al. 2012; Gandhi et al. 2013; Jones et al. 2007; Kim et al. 2014; Li et al. 2013; Ma et al. 2014; Polovlkova et al. 2013; Rosen et al. 2013; Savas et al. 2009; Singer et al. 2012; Turco et al. 2012; Zhang et al. 2009). Moreover, the latest studies demonstrate that GnRH are associated with regulation of proliferation, angiogenesis and inflammatory response, which are important in the pathogenesis of diseases associated with aging (Pincas et al. 2014). There are also some data about the negative relation between GnRH receptors and growth factors, which play a part in aging processes related to activity of stem cells, tissue regeneration, protein homeostasis and regulation of inflammatory processes leading to the age-related diseases (Cheung and Wong 2008). So far, there are only few reports about the influence of GnRH on the aging process and aging-associated diseases. De Haes et al. (2013) showed that silencing of the gene encoding the GnRH receptor resulted in prolongation of C. elegans life. In addition, Zhang et al. (2013) pointed out the relationship between the pituitary GnRH and immunoaging.

The second pathway which was influenced by BASE-diet is guanylate cyclase pathway. Guanylyl cyclase, in response to calcium levels, synthesizes 3′,5′-cyclic guanosine monophosphate (cGMP) from guanosine triphosphate (GTP). cGMP is associated with processes, such as phototransduction, circadian entrainment, olfactory transduction, vascular smooth muscle contraction, gap junction, long-term depression, salivary secretion, borderline and spontaneous hypertension, platelet activation and learning ability (Cornilescu et al. 2014; Arshavsky and Wensel 2013; Chae et al. 2007; Chen et al. 2007; Golombek et al. 2003; Wang et al. 2013; Dismuke et al. 2013; Dam et al. 2014; Mao et al. 2013; Kameritsch et al. 2012; Robinson et al. 2012; Kawaguchi and Hirano 2013; Simon et al. 2013; Henkin et al. 2007, 2013; Henkin and Velicu 2009; Pechánová et al. 2009; Rukoyatkina et al. 2011; Rodrigo et al. 2006). However, none of these processes were earlier described as related to the liver functions and the metabolism of biologically active compounds. Some data suggest a direct link between guanylate cyclase pathway and gonadotrope cell activation, since there is a link between cGMP and genes involved in reproduction regulations, especially PRKG1 and NOS1. It has been demonstrated that cGMP analogues stimulate GnRH release, while a cGMP-dependent protein kinase (PRKG1) inhibitor conversely blocks nitric oxide-induced GnRH release. Therefore, nitric oxide may regulate the heme-containing signalling enzyme, guanylate cyclase, and thereby elevate the second messenger cGMP and facilitate GnRH secretion. Additionally, an effect of nitric oxide on the heme-containing enzyme should also be taken into account. Cyclooxygenase controls production of prostaglandins, known stimulators of GnRH secretion (Brann and Mahesh 1997).

Taken together, the present study showed that long-term (3 and 14 months) rats feeding with BASE-diet (diet enriches with polyphenolic compounds, β-carotene, probiotics and n-3 and n-6 polyunsaturated fatty acids) can affect liver genes expression. On transcriptomic level, the diet exerted its activity by modulating the expression of genes involved particularly in the gonadotrope cell activation pathway and guanylate cyclase pathway, as well as in mast cell activation, gap junction regulation, melanogenesis and apoptosis. The results indicating the strong influence of BASE-diet on genes involved in the gonadotrope cell activation pathway may suggest the impact of BASE-diet on reproduction system at the transcriptome level. BASE-diet especially strongly affects genes involved in regulation of GnRH which is responsible for the release of FSH and LH. This effect is stronger with the age of animals and the length of diet use. It can suggest a link between the diet, reproductive system function and aging. The reproductive-cell cycle theory assumes that the hormones that regulate reproduction, mainly LH and FSH, act in an antagonistic pleiotrophic manner to control aging through cell cycle signalling (Atwood and Bowen 2011). It means that the improvement of functioning of reproductive system may slow down the rate of senescence, thereby decelerating the rate of aging, and thus the lifespan. It allows us to draw the conclusion that the long-term use of biologically active substances-enriched diet can positively affect the reproductive system, and thus, as a consequence may delay the aging process. However, according to the latest scientific reports, the gonadotrope cell activation pathway plays a part not only in reproduction system, but also in regulation of the aging process. Therefore, the results indicating the strong influence of BASE-diet on genes involved in the gonadotrope cell activation pathway may also suggest the strong impact of BASE-diet on the aging process by the regulation of GnRH-dependent cells. These results indicate that it is highly probable that BASE-diet can modify the signalling pathways which control the aging process by changing the expression of genes involved in gonadotrope cell activation and as a result delay the flow of the aging process.


  1. Alam SE, Singh RB, Gupta S, Dherange P, De Meester F, Wilczynska A, Dharwadkar S, Wilson D, Hungin P (2012) Nutritional aspects of epigenetic inheritance. Can J Physiol Pharmacol 90:989–994

  2. Ames HM, Wang AA, Coughran A, Evaul K, Huang S, Graves CW, Soyombo AA, Ross TS (2013) Huntingtin-interacting protein 1 phosphorylation by receptor tyrosine kinases. Mol Cell Biol 33:3580–3593

  3. Arechiga CF, Staples CR, McDowell LR, Hansen PJ (1998) Effects of timed insemination and supplemental beta-carotene on reproduction and milk field of dairy cows under heat stress. J Dairy Sci 81:390–402

  4. Arellano-Rodriguez G, Meza-Herrera CA, Rodriguez-Martinez R, Dionisio-Tapia R, Hallford DM, Mellado M, Gonzalez-Bulnes A (2009) Short-term intake of beta-carotene-supplemented diets enhances ovarian function and progesterone synthesis in goats. J Anim Physiol Anim Nutr 93(6):710–715. doi:10.1111/j.1439-0396.2008.00859.x

  5. Arshavsky VY, Wensel TG (2013) Timing is everything: GTPase regulation in phototransduction. Invest Ophthalmol Vis Sci 54(12):7725–7733. doi:10.1167/iovs.13-13281

  6. Astley SB (2007) An introduction to nutrigenomics developments and trends. Genes Nutr 2:11–13

  7. Atwood CS, Bowen RL (2011) The reproductive-cell cycle theory of aging: an update. Exp Gerontol 46(2–3):100–107. doi:10.1016/j.exger.2010.09.007

  8. Bakker GC, van Erk MJ, Pellis L, Wopereis S, Rubingh CM, Cnubben NH, Kooistra T, van Ommen B, Hendriks HF (2010) An antiinflammatory dietary mix modulates inflammation and oxidative and metabolic stress in overweight men: a nutrigenomics approach. Am J Clin Nutr 91:1044–1059

  9. Berry DC, Jin H, Majumdar A, Noy N (2011) Signaling by vitamin A and retinol-binding protein regulates gene expression to inhibit insulin responses. Proc Natl Acad Sci USA 108(11):4340–4345. doi:10.1073/pnas.1011115108

  10. Bilby TR, Guzeloglu A, MacLaren LA, Staples CR, Thatcher WW (2006) Pregnancy, bovine somatotropin, and dietary n-3 fatty acids in lactating dairy cows: II. Endometrial gene expression related to maintenance of pregnancy. J Dairy Sci 89(9):3375–3385

  11. Brann DW, Mahesh VB (1997) Excitatory amino acids: evidence for a role in the control of reproduction and anterior pituitary hormone secretion. Endocr Rev 18(5):678–700

  12. Chae KS, Ko GY, Dryer SE (2007) Tyrosine phosphorylation of cGMP-gated ion channels is under circadian control in chick retina photoreceptors. Invest Ophthalmol Vis Sci 48(2):901–906

  13. Chaudhry AS, Thirumaran RK, Yasuda K, Yang X, Fan Y, Strom SC, Schuetz EG (2013) Genetic variation in aldo-keto reductase 1D1 (AKR1D1) affects the expression and activity of multiple cytochrome P450s. Drug metabolism and disposition: the biological fate of chemicals 41:1538–1547

  14. Chen SK, Ko GY, Dryer SE (2007) Somatostatin peptides produce multiple effects on gating properties of native cone photoreceptor cGMP-gated channels that depend on circadian phase and previous illumination. J Neurosci 27(45):12168–12175

  15. Chen ZG, Luo LL, Xu JJ, Zhuang XL, Kong XX, Fu YC (2010) Effects of plant polyphenols on ovarian follicular reserve in aging rats. Biochem Cell Biol 88(4):737–745. doi:10.1139/O10-012

  16. Cheung LW, Wong AS (2008) Gonadotropin-releasing hormone: GnRH receptor signaling in extrapituitary tissues. The FEBSJ 275:5479–5495

  17. Colitti M, Stefanon B (2006) Effect of natural antioxidants on superoxide dismutase and glutathione peroxidase mRNA expression in leukocytes from periparturient dairy cows. Vet Res Commun 30:19–27

  18. Colitti M, Sgorlon S, Stradaioli G, Farinacci M, Gabai G, Stefanon B (2007) Grape polyphenols affect mRNA expression of PGHS-2, TIS11b and FOXO3 in endometrium of heifers under ACTH-induced stress. Theriogenology 68(7):1022–1030

  19. Cornilescu CC, Cornilescu G, Burgie ES, Markley JL, Ulijasz AT, Vierstra RD (2014) Dynamic structural changes underpin photoconversion of a blue/green cyanobacteriochrome between its dark and photoactivated states. J Biol Chem 289(5):3055–3065

  20. Costello LM, O’Boyle P, Godkin JD, Diskin MG, Hynes AC, Morris DG (2010) Retinol-binding protein (RBP), retinol and β-carotene in the bovine uterus and plasma during the oestrouscycle and the relationship between systemic progesterone and RBP on day 7. Reprod Fertil Dev 22(8):1198–1205. doi:10.1071/RD10034

  21. Dam VS, Boedtkjer DM, Nyvad J, Aalkjaer C, Matchkov V (2014) TMEM16A knockdown abrogates two different Ca2+-activated Cl− currents and contractility of smoothmuscle in rat mesenteric small arteries. Pflugers Arch Eur J Physiol 466(7):1391–1409

  22. Datta K, Hyduke DR, Suman S, Moon BH, Johnson MD, Fornace AJ Jr (2012) Exposure to ionizing radiation induced persistent gene expression changes in mouse mammary gland. Radiation Oncology 7:205

  23. De Haes W, Frooninckx L, Van Rompay L, Temmerman L, Schoofs L (2013) Disrupting GNRH/AKH signalling in C. elegans promotes longevity. In: 17th International congress of comparative endocrinology

  24. De Stefano D, Maiuri MC, Simeon V, Grassia G, Soscia A, Cinelli MP, Carnuccio R (2007) Lycopene, quercetin and tyrosol prevent macrophage activation induced by gliadin and IFN-gamma. Eur J Pharmacol 566:192–199

  25. Dismuke WM, Liang J, Overby D, Stamer WD (2013) Concentration-related effects of nitric oxide and endothelin-1 on human trabecular meshwork cell contractility. Exp Eye Res. doi:10.1016/j.exer.2013.12.012

  26. Ferguson LR, Peterman I, Hubner C, Philpott M, Shellin AN (2007a) Uncoupling gene-diet interactions in inflammatory bowel disease (IBD). Genes Nutr 2:71–73

  27. Ferguson LR, Shelling AN, Lauren D (2007b) For Nutrigenomics New Zealand. Nutrigenomics and gut health: meeting report from an international conference in Auckland, New Zealand. Genes Nutr 2:157–160

  28. Ferretti G, Bacchetti T, Masciangelo S, Saturni L (2012) Celiac disease, inflammation and oxidative damage. A nutrigenetic approach. Nutrients 4(4):243–257. doi:10.3390/nu404024

  29. Gandhi AV, Saxena S, Relles D, Sarosiek K, Kang CY, Chipitsyna G, Sendecki JA, Yeo CJ, Arafat HA (2013) Differential expression of cytochrome P450 omega-hydroxylase isoforms and their association with clinicopathological features in pancreatic ductal adenocarcinoma. Ann Surg Oncol 20(Suppl 3):S636–S643

  30. Garcia-Bojalil CM, Staples CR, Risco CA, Savio JD, Thatcher WW (1998) Protein degradability and calcium salts of Long-chain fatty acids in the diets of lactating Dairy cows: reproductive responses. J Dairy Sci 81:1385–1395

  31. Garrel G, Simon V, Denoyelle C, Cruciani-Guglielmacci C, Migrenne S, Counis R, Magnan C, Cohen-Tannoudji J (2011) Unsaturated fatty acids stimulate LH secretion via novel PKCepsilon and -theta in gonadotrope cells and inhibit GnRH-induced LH release. Endocrinology 152(10):3905–3916. doi:10.1210/en.2011-1167

  32. Gebrie E, Makonnen E, Debella A, Zerihun L (2005) Phytochemical screening and pharmacological evaluations for the antifertility effect of the methanolic root extract of Rumex steudelii. J Ethnopharmacol 96(1–2):139–143

  33. Golombek DA, Ferreyra GA, Agostino PV, Murad AD, Rubio MF, Pizzio GA, Katz ME, Marpegan L, Bekinschtein TA (2003) From light to genes: moving the hands of the circadian clock. Front Biosci 8:s285–s293

  34. Guyton AC, Hall JE (2006) Textbook of medical physiology. Elsevier Saunders, Philadelphia

  35. Haliloglu S, Baspinar N, Serpek B, Erdem H, Bulut Z (2002) Vitamin A and beta-carotene levels in plasma, corpus luteum and follicular fluid of cyclic and pregnant cattle. Domest Anim Reprod 37:96–99

  36. Hannah VC, Ou J, Luong A, Goldstein JL, Brown MS (2001) Unsaturated fatty acids downregulate SREBP isoforms 1a and 1c by two mechanisms in HEK-293 cells. J Biol Chem 276:4365–4372

  37. Henkin RI, Velicu I (2009) Decreased parotid salivary cyclic nucleotides related to smell loss severity in patients with taste and smell dysfunction. Metabolism 58(12):1717–1723. doi:10.1016/j.metabol.2009.05.027

  38. Henkin RI, Velicu I, Papathanassiu A (2007) cAMP and cGMP in human parotid saliva: relationships to taste and smell dysfunction, gender, and age. Am J Med Sci 334(6):431–440

  39. Henkin RI, Levy LM, Fordyce A (2013) Taste and smell function in chronic disease: a review of clinical and biochemical evaluations of taste and smell dysfunction in over 5000 patients at The Taste and Smell Clinic in Washington, DC. Am J Otolaryngol 34(5):477–489. doi:10.1016/j.amjoto.2013.04.006

  40. Huang CZ, Huang WZ, Zhang G, Tang DL (2013) In vivo study on the effects of curcumin on the expression profiles of anti tumour genes (VEGF, CyclinD1 andCDK4) in liver of rats injected with DEN. Mol Biol Rep 40:5825–5831

  41. Hutchinson IA, Hennessy AA, Waters SM, Dewhurst RJ, Evans AC, Lonergan P, Butler ST (2012) Effect of supplementation with different fat sources on the mechanisms involved in reproductive performance in lactating dairy cattle. Theriogenology 78(1):12–27. doi:10.1016/j.theriogenology.2011.12.031

  42. Jones N, Hardy WR, Friese MB, Jorgensen C, Smith MJ, Woody NM, Burden SJ, Pawson T (2007) Analysis of a Shc family adaptor protein, ShcD/Shc4, that associates with muscle-specific kinase. Mol Cell Biol 27:4759–4773

  43. Joseph J, Cole G, Head E, Ingram D (2009) Nutrition, brain aging, and neurodegeneration. J Neurosci 29:12795–12801

  44. Kameritsch P, Pogoda K, Ritter A, Münzing S, Pohl U (2012) Gap junctional communication controls the overall endothelial calcium response to vasoactive agonists. Cardiovasc Res 93(3):508–515. doi:10.1093/cvr/cvr345

  45. Kawaguchi SY, Hirano T (2013) Gating of long-term depression by Ca2+/calmodulin-dependent protein kinase II through enhanced cGMPsignalling in cerebellar Purkinje cells. J Physiol 591(Pt 7):1707–1730. doi:10.1113/jphysiol.2012.245787

  46. Kim HK, Kim MG, Leem KH (2014) Collagen hydrolysates increased osteogenic gene expressions via a MAPK signaling pathway in MG-63 human osteoblasts. Food Funct 5:573–578

  47. Lee VH, Lee LT, Chow BK (2008) Gonadotropin-releasing hormone: regulation of the GnRH gene. FEBS J 275(22):5458–5478. doi:10.1111/j.1742-4658.2008.06676.x

  48. Li G, Zhao F, Cui Y (2013) Proteomics using mammospheres as a model system to identify proteins deregulated in breast cancer stem cells. Curr Mol Med 13:459–463

  49. Lillycrop KA, Burdge GC (2012) Epigenetic mechanisms lin king early nutrition to long term health. Best Pract Res Clin Endocrinol Metab 26:667–676

  50. Liu Q, Bengmark S, Qu S (2010) Nutrigenomics therapy of hepatisis C virus induced-hepatosteatosis. BMC Gastroenterol 10:49. doi:10.1186/1471-230X-10-49

  51. Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods 25:402–408

  52. Lottenberg AM, Afonso Mda S, Lavrador MS, Machado RM, Nakandakare ER (2012) The role of dietary fatty acids in the pathology of metabolic syndrome. J Nutr Biochem 23(9):1027–1040. doi:10.1016/j.jnutbio.2012.03.004

  53. Lucy MC, De La Sota RL, Staples CR, Thatcher WW (1993) Ovarian follicular populations in lactating Dairy cows treated with recombinant bovine somatotropin (sometribove) or saline and Fed Diets Differing in fat content and energy. J Dairy Sci 76:1014–1027

  54. Ma X, Pietsch J, Wehland M, Schulz H, Saar K, Hubner N, Bauer J, Braun M, Schwarzwalder A et al (2014) Differential gene expression profile and altered cytokine secretion of thyroid cancer cells in space. FASEB J 28:813–835

  55. Mao GK, Li JX, Bian FH, Han YY, Guo M, Xu BS, Zhang MJ, Xia GL (2013) Gap junction-mediated cAMP movement between oocytes and somatic cells. Front Biosci 5:755–767

  56. Matteri RL, Carroll JA, Dyer CJ (2001) Neuroendocrine response to stress. In: Mench JA, Moberg GP (eds) The biology of Animals stress. Basic principles and implications for animal welfare. CABI Publishing, New York, pp 43–76

  57. Mattos R, Staples CR, Williams J, Amorocho A, McGuire MA, Thatcher WW (2002) Uterine, ovarian, and production responses of lactating Dairy cows to increasing dietary concentrations of menhaden fish meal. J Dairy Sci 85:755–764

  58. Millar RP (2005) GnRHs and GnRH receptors. Anim Reprod Sci 88(1–2):5–28

  59. Okamura H, Tsukamura H, Ohkura S, Uenoyama Y, Wakabayashi Y, Maeda K (2013) Kisspeptin and GnRH pulse generation. Adv Exp Med Biol 784:297–323

  60. Oszkiel H, Wilczak J, Prostek A, Kamola D, Hulanicka M, Jank M (2014) Does long-term administration of biologically active substances-enriched diet influence general health status parameters in rats? J Aging Res Clin Pract 3(2):63–71

  61. Pechánová O, Jendeková L, Vranková S (2009) Effect of chronic apocynin treatment on nitric oxide and reactive oxygen species production In borderline and spontaneous hypertension. Pharmacol Rep 61(1):116–122

  62. Pincas H, Choi SG, Wang Q, Jia J, Turgeon JL, Sealfon SC (2014) Outside the box signaling: Secreted factors modulate GnRH receptor-mediated gonadotropin regulation. Mol Cell Endocrinol 385(1–2):56–61

  63. Polovlkova OG, Makeeva OA, Lezhnev AA, Goncharova IA, Kulish EV, Shipulin VM, Puzyrev VP (2013) Expression profile of calcineurin pathway genes in myocardium tissues in relation to ischemic heart remodeling in humans. Mol Biol 47:433–440

  64. Robinson JW, Zhang M, Shuhaibar LC, Norris RP, Geerts A, Wunder F, Eppig JJ, Potter LR, Jaffe LA (2012) Luteinizing hormone reduces the activity of the NPR2 guanylyl cyclase in mouse ovarian follicles, contributing to the cyclic GMP decrease that promotes resumption of meiosis in oocytes. Dev Biol 366(2):308–316. doi:10.1016/j.ydbio.2012.04.019

  65. Rodrigo R, Monfort P, Cauli O, Erceg S, Felipo V (2006) Pharmacological manipulation of cyclic GMP levels in brain re stores learning ability in animal models of hepatic encephalopathy: therapeutic implications. Neuropsychiatr Dis Treat 2(1):53–63

  66. Rosen MB, Das KP, Wood CR, Wolf CJ, Abbott BD, Lau C (2013) Evaluation of perfluoroalkyl acid activity using primary mouse and human hepatocytes. Toxicology 308:129–137

  67. Rukoyatkina N, Walter U, Friebe A, Gambaryan S (2011) Differentiation of cGMP-dependent and -independent nitric oxide effects on platelet apoptosis and reactive oxygen species production using platelets lacking soluble guanylyl cyclase. Thromb Haemost 106(5):922–933. doi:10.1160/TH11-05-0319

  68. Savas U, Machemer DE, Hsu MH, Gaynor P, Lasker JM, Tukey RH, Johnson EF (2009) Opposing roles of peroxisome proliferator-activated receptor alpha and growth hormone in the regulation of CYP4A11 expression in a transgenic mouse model. J Biol Chem 284:16541–16552

  69. Schweigert FJ, Steihagen B, Raila J, Seimann A, Peet D, Buscher U (2003) Concentrations of carotenoids, retinol and a-tocopherol in plasma and follicular fluid of woman undergoing IVF. Hum Reprod 18:1259–1264

  70. Sgorlon S, Stradaioli G, Zanin D, Stefanon B (2006) Biochemical and molecular responses to antioxidant supplementation in sheep. Small Rumin Res 64:143–151

  71. Sharma S, Morinaga H, Hwang V, Fan W, Fernandez MO, Varki N, Olefsky JM, Webster NJ (2013) Free fatty acids induce Lhb mRNA but suppress Fshb mRNA in pituitary LβT2 gonadotropes and diet-induced obesity reduces FSH levels in male mice and disrupts the proestrous LH/FSH surge in female mice. Endocrinology 154(6):2188–2199. doi:10.1210/en.2012-2218

  72. Simon A, Karbach S, Habermeier A, Closs EI (2013) Decoding the substrate supply to human neuronal nitric oxide synthase. PLoS One 8(7):e67707. doi:10.1371/journal.pone.0067707

  73. Singer A, Deuse Y, Koch U, Holscher T, Pfitzmann D, Jakob C, Hehlgans S, Baretton GB, Rentsch A et al (2012) Impact of the adaptor protein GIPC1/Synectin on radioresistance and survival after irradiation of prostate cancer. Strahlenther Onkol 188:1125–1132

  74. Subcommittee on Laboratory Animal Nutrition, Committee on Animal Nutrition, Board on Agriculture, National Research Council (1995) Nutrient requirements of laboratory animals, 4th edn. National Academy Press, Washington, DC

  75. Turco MY, Furia L, Dietze A, Fernandez Diaz L, Ronzoni S, Sciullo A, Simeone A, Constam D, Faretta M et al (2012) Cellular heterogeneity during embryonic stem cell differentiation to epiblast stem cells is revealed by the ShcD/RaLP adaptor protein. Stem Cells 30:2423–2436

  76. Uehara Y, Miura S, von Eckardstein A, Abe S, Fujii A, Matsuo Y et al (2007) Unsaturated fatty acids suppress the expression of the ATP-binding cassette transporter G1 (ABCG1) and ABCA1 genes via an LXR/RXR responsive element. Atherosclerosis 191:11–21

  77. Van Arnum SD (1998) Vitamin A. Vitamin A in Kirk-Othmer Encyclopedia of Chemical Technology, vol 45. John Wiley, New York, pp 99–107

  78. Virmani A, Binienda Z, Ali S, Gaetani F (2006) Links between nutrition, drug abuse, and the metabolic syndrome. Ann N Y Acad Sci 1074:303–314

  79. Virmani A, Pinto L, Binienda Z, Ali S (2013) Food, nutrigenomics, and neurodegeneration-neuroprotection by what you eat! Mol Neurobiol 48(2):353–362. doi:10.1007/s12035-013-8498-3

  80. Wang L, Espinoza HM, Gallagher EP (2013) Brief exposure to copper induces apoptosis and alters mediators of olfactory signal transduction in coho salmon. Chemosphere 93(10):2639–2643. doi:10.1016/j.chemosphere.2013.08.044

  81. Zhang WY, Liu HQ, Xie LL, Li Y, Kwik-Uribe CL, Zhu XZ (2006) Procyanidin dimer B2 [epicatechin-(4beta-8)-epicatechin] suppresses the expression of cyclooxygenase-2 in endotoxin-treated monocytic cells. Biochem Biophys Res Commun 345:508–515

  82. Zhang Q, Li Y, Liu G, Xu X, Song X, Liang B, Li R, Xie J, Du M et al (2009) Alternative transcription initiation and splicing variants of the DHRS4 gene cluster. Biosci Rep 29:47–56

  83. Zhang G, Li J, Purkayastha S, Tang Y, Zhang H, Yin Y, Li B, Liu G, Cai D (2013) Hypothalamic programming of systemic ageing involving IKK-β, NF-κB and GnRH. Nature 497(7448):211–216

Download references


This work was supported by Grant No. N N312 337939 from the Ministry of Sciences and Higher Education.

Conflict of interest

Hanna Oszkiel, Jacek Wilczak and Michał Jank declare that they have no conflict of interest.

Ethical standard

All institutional and national guidelines for the care and use of laboratory animals were followed.

Author information

Correspondence to Michał Jank.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Oszkiel, H., Wilczak, J. & Jank, M. Biologically active substances-enriched diet regulates gonadotrope cell activation pathway in liver of adult and old rats. Genes Nutr 9, 427 (2014).

Download citation


  • Rat
  • Liver
  • Biologically active substances
  • Genes expression
  • GnRH