# Fenugreek Stimulates the Expression of Genes Involved in Milk Synthesis and Milk Flow through Modulation of Insulin/GH/IGF-1 Axis and Oxytocin Secretion

Thomas Sevrin, Clair-Yves Boquien, Alexis Gandon, Isabelle Grit, Pierre de Coppet, Dominique Darmaun, Marie Cécile Alexandre-Gouabau

► To cite this version:

Thomas Sevrin, Clair-Yves Boquien, Alexis Gandon, Isabelle Grit, Pierre de Coppet, et al. Fenugreek Stimulates the Expression of Genes Involved in Milk Synthesis and Milk Flow through Modulation of Insulin/GH/IGF-1 Axis and Oxytocin Secretion. Genes, 2020, 11 (10), 10.3390/genes11101208. [https://hal-inrae.fr/hal-03209164](https://hal.inrae.fr/hal-03209164)

HAL Id: hal-03209164

https://hal.inrae.fr/hal-03209164

Submitted on 27 Apr 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Fenugreek Stimulates the Expression of Genes Involved in Milk Synthesis and Milk Flow through Modulation of Insulin/GH/IGF-1 Axis and Oxytocin Secretion

Thomas Sevrin<sup>1,2 ID</sup>, Clair-Yves Boquien<sup>2 ID</sup>, Alexis Gandon<sup>2</sup>, Isabelle Grit<sup>2</sup>, Pierre de Coppet<sup>2</sup>, Dominique Darmaun<sup>2,3</sup> and Marie-Cécile Alexandre-Gouabau<sup>2,* ID</sup>

Received: 18 September 2020; Accepted: 13 October 2020; Published: 16 October 2020

<figure><img src="image_2.png" /></figure>

check for
updates

**Abstract:** We previously demonstrated galactagogue effect of fenugreek in a rat model of lactation challenge, foreshadowing its use in women's breastfeeding management. To assess longitudinal molecular mechanisms involved in milk synthesis/secretion in dams submitted to fenugreek supplementation, inguinal mammary, pituitary glands and plasma were isolated in forty-three rats nursing large 12 pups-litters and assigned to either a control (CTL) or a fenugreek-supplemented (FEN) diet during lactation. RT-PCR were performed at days 12 and 18 of lactation (L12 and L18) and the first day of involution (Inv1) to measure the relative expression of genes related to both milk synthesis and its regulation in the mammary gland and lactogenic hormones in the pituitary gland. Plasma hormone concentrations were measured by ELISA. FEN diet induced 2- to 3-times higher fold change in relative expression of several genes related to macronutrient synthesis (*Fasn*, *Acaca*, *Fabp3*, *B4galt1*, *Lalba* and *Csn2*) and energy metabolism (*Cpt1a*, *Acads*) and in IGF-1 receptor in mammary gland, mainly at L12. Pituitary oxytocin expression and plasma insulin concentration (+77.1%) were also significantly increased. Altogether, these findings suggest fenugreek might extend duration of peak milk synthesis through modulation of the insulin/GH/IGF-1 axis and increase milk ejection by activation of oxytocin secretion.

**Keywords:** fenugreek; galactagogue; milk synthesis; gene expression; lactating mammary gland; pituitary gland; insulin; oxytocin

## 1. Introduction

The World Health Organization (WHO) recommends exclusive breastfeeding for infants up to 6 months of age, based on evidence of its clear health benefits on mother-infant dyad [1–3]. Yet, 6 months after delivery, less than 40% of mothers are still breastfeeding in several high-income countries of North America and Europe [4]. Although the early cessation of breastfeeding is multifactorial [5], the perception that their milk secretion is insufficient to support adequate infant growth is reported by about 35% of lactating women [6,7]. Though perceived milk insufficiency due to psychological

<footnote n="1">1 FRANCE Bébé Nutrition (FBN) Laboratory, 53000 Laval, France; thomas.sevrin@etu.univ-nantes.fr</footnote>

<footnote n="2">2 Mixed Research Unit 1280 Pathophysiology of Nutritional adaptations (UMR 1280 PhAN) Nantes University, Research Center in Human Nutrition-West (CRNH-O), Institute of Digestive Tract Diseases (IMAD), French National Research Institute for Agriculture, Food and Environment (INRAE), F-44000 Nantes, France; clair-yves.boquien@univ-nantes.fr (C.-Y.B.); alexis.gandon@univ-nantes.fr (A.G.); isabelle.grit@univ-nantes.fr (I.G.); Pierre.De-Coppet@univ-nantes.fr (P.d.C.); dominique.darmaun@univ-nantes.fr (D.D.)</footnote>

<footnote n="3">3 Nantes University Hospital (CHU) Nantes, F-44000 Nantes, France</footnote>

<footnote>* Correspondence: Marie-Cecile.Alexandre-Gouabau@univ-nantes.fr; Tel.: +33-(0)2-53-48-20-12</footnote>
issues likely is more frequent than true insufficient milk secretion [5,8], the latter can result from many
causes, ranging from true physiologic inability to lactate (5% of the cases) to suboptimal breastfeeding
management [8,9]. It is well established that mothers' milk production can often be increased through
psychological support or maternal breastfeeding counseling [5,9]. Nevertheless, more and more healthy
mothers try to enhance their milk supply by taking various nutritional supplements [10] presumed to
have a galactagogue effect.

In the absence of a reference galactagogue molecule, the use of plant extracts has grown rapidly [10], with nowadays 15% to 25% of mothers using herbal remedies to increase their milk supply in Western countries [11,12]. Fenugreek is the most consumed, accounting for 50% of the plant extracts taken during breastfeeding [11], although strong evidence of its effectiveness and safety is scarce [11,13]. In a recent study using the deuterium oxide enrichment method to accurately measure milk secretion in a model of rat nursing large litters [14], we found that dietary fenugreek at a dose of 1 g/kg body weight/day increased milk production by 16% and offspring growth by 11% without any evidence for adverse metabolic effect on either dams or offspring [15], demonstrating a promising potential of fenugreek in the treatment of breastfeeding deficiency. However, the mechanisms underlying its effect remain poorly understood [16] and deciphering the targets of fenugreek in the lactation process should allow initiating careful support to mothers with lactation difficulties.

Milk production is mainly achieved by complex interactions signals both from the endocrine
pathway and the local mammary gland environment. Indeed, in response to nipple areola stimulation
during suckling, a nervous message is sent to the hypothalamus that stimulates the secretion of both
prolactin and oxytocin, the main milk secretory hormones, at the pituitary level. Prolactin stimulates
the synthesis of main milk macronutrients by binding to its receptors located on the surface of lactocytes
(milk secretory cells). Oxytocin stimulates milk ejection by inducing the contraction of myoepithelial
cells surrounding mammary secretory alveoli [17,18]. Insulin also plays an essential role by stimulating
mainly protein and lactose synthesis in the mammary gland [19]. Finally, growth hormone (GH) and
its effector insulin-like growth factor 1 (IGF-1), and thyroid hormones, can stimulate milk production
directly in the mammary gland or by redirecting the nutrient flow toward breast tissues [20-22].

Among the therapeutic properties of its numerous bioactive constituents [23–25], such as trigonelline or diosgenin, fenugreek may increase prolactin secretion [26] due to an estrogenic action, which decreases the secretion of dopamine (a prolactin secretion inhibitor) in the hypothalamus [26–28]. Fenugreek may also increase food consumption through inhibition of leptin secretion [29]. Some fenugreek constituents had also been shown to increase GH secretion by rat pituitary cells [30]. In addition, fenugreek is a well-known hypoglycemic agent [25] due to its capacity to increase insulin secretion by β-pancreatic cells [24]. Finally, fenugreek is rich in antioxidant compounds [24,25] that could sustain an optimal mammary gland function. We hypothesized that fenugreek may increase milk production through an effect on the regulation of lactation.

The objectives of the present study were (i) to determine whether fenugreek impacted genes responsible for the synthesis of milk constituents and milk production, particularly candidate genes present in mammary and pituitary glands and responsible for the regulation of lactogenesis and milk ejection, and (ii) to determine the time course of the observed effects on the candidate genes at 3 time points: mid-lactation, late lactation and early mammary gland involution.

**2. Materials and Methods**

2.1. *Animal Experiment*

2.1.1. Housing and Diets

Pregnant Sprague-Dawley rats, purchased from Janvier Labs (Le Genest-Saint-Isle, France) at first gestational day (G1), were housed individually in cages with wood chips located on ventilated racks kept at a constant temperature of 22 ± 1 °C and at a relative humidity of 50% ± 3%. Cages were placed
in a room with a fixed 12 h light–dark cycle (light from 7:00 a.m. to 7:00 p.m.). Pregnant rats had access to water and food ad libitum.

During gestation, every dam received a control diet (CTL) based on AIN-93G diet with 20 g protein per 100 g of food. At parturition (G21 = L0), dams received ad libitum either the control diet (CTL) alone or the CTL diet supplemented with dry water extract of fenugreek seeds (Plantex, Sainte-Geneviève-des-Bois, France) (FEN) at 1 g/kg body weight/day as previously described [15]. The diets were manufactured by the "Experimental Food Preparation Unit" (INRAE-UPAE, Jouy-en-Josas, France).

2.1.2. Experimental Design

The experimental protocol was approved by the Animal Ethics committee of Pays de La Loire and registered under reference (protocol APAFIS 2018121018129789) (Angers, France, 10 December 2018). All animal procedures were carried out in accordance with current institutional guidelines on animal experimentation in the EU (directive 2010/63/EU), in France (French Veterinary Department—A44276) and the French Ministry of Research in compliance with the commonly-accepted "3Rs".

The study consisted of a succession of 3 separate experiments (Xp1, Xp2 and Xp3) with 22, 12 and 11 individual litters respectively, as shown in Figure 1. In the three experiments, all dams received CTL diet during gestation. Delivery occurred at G21 also considered as day 0 of lactation (L0) and dams were randomly assigned either to the CTL or FEN diet throughout lactation. Pups were pooled by sex and randomly adopted by CTL or FEN dams. Litters were homogenized to 12 pups with a sex ratio of 1/1 and balanced weight of 84 ± 2 g (Figure 1). The 3 experiments differed by the time at which animal follow up was stopped: post-natal day 21 for Xp1, L18 for Xp2 and L12 for Xp3. In Xp1 and Xp2, milk flow between dams and their pups was recorded between L11 and L18 using the deuterated water enrichment method and milk sample was collected at L18. In the Xp1, as previously described [14], pups were weaned at L20 to start the involution phase for dams’ mammary glands. At post-natal day 21, also assimilated to first involution day (Inv1), dams were fasted for 4 h before sacrifice and sampling of plasma and tissues. In Xp2, after milk sampling at L18, food and pups were removed from the cage for 2.5 h and 1 h respectively, before the beginning of dams’ sacrifice to avoid bias linked to hormone secretion during suckling. Finally, in Xp3, milk was sampled at L11 and dams were sacrificed at L12 after 1.5 h fasting and 1 h removal of pups (Figure 1).

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Experimental design: three series of animal experiments (Xp1, 2 and 3) were carried out, all with dams receiving either control (CTL) or fenugreek-supplemented diet (FEN), with different period of sampling for each experiment. G, L, PND and Inv correspond to the day of gestation, lactation, post-natal and involution days, respectively. Milk flow was measured by the deuterated water enrichment method [14,15] between L11 and L18. Tissues sampled were mammary and pituitary glands.</figcaption></figure>
### 2.1.3. Dam-Litter Dyad Follow-Up and Milk Flow Measurement

Throughout lactation period, the dam weight, food and water consumption and total litter weight were recorded every two days from L0 to L11, and every day from L11 to L18. Pup weight was calculated as litter weight divided by litter size (12) and weight gain by subtracting weight at L0 from daily weight.

Milk flow was measured by the water turnover method as previously described [14,15]. Briefly, at L11, after 4%, isoflurane anesthesia, mothers received an injection of 4.92 ± 0.10 g/kg body weight D₂O (99.9 mole% D₂-enrichment, Sigma-Aldrich, Lyon, France) into a tail vein. Plasma samples were collected from dams by tail snip from L11 to L15 and pooled urine samples from their own whole litters were collected from L12 to L18 by stimulating lower belly with an iced cotton bud. The D₂O enrichment of both plasma and urine samples was measured using the Fourier Transform infrared spectrophotometer α II® (Brucker, Rheinstetten, Germany). Milk flow calculation was performed using a bi-compartmental model [15] as no significant intra-litter sexual dimorphism was observed in previous experiments [14]. Flow constants of the model (h⁻¹) were then determined with the SAAM II® software (The Epsilon Group, Charlottesville, Virginia, US) and the dams' total body water (g) using the intercept of dams D₂O concentration curve. The water flow, from dams to their own litter (g/h), was calculated as the product of total body water by flow constant from dam to the litter and was expressed in g/day by multiplying by 24.

### 2.1.4. Milk, Plasma and Tissue Sampling

After 20 min dam/pups' separation and intraperitoneal injection of oxytocin (1 unit of Syntocinon®; Sigma-Tau, Ivry-sur-Seine, France) to stimulate milk ejection, dams were anaesthetized with 4% of isoflurane, and milk was manually collected and stored at -20 °C until analysis. Blood samples were collected before animal sacrifice via intra-cardiac puncture and placed into Ethylene Diamine Tetraacetic Acid (EDTA) tubes (Pfizer-Centravet, Plancoët, France), centrifuged at 1132 g for 15 min at 4 °C, and plasma were stored at -20 °C until analysis. Dams were sacrificed by intra-cardiac injection of 0.5 mL Exagon® (Richter pharma, Wels, Austria). Left inguinal mammary gland and pituitary gland were removed, sampled and immediately frozen in liquid nitrogen before storage at -80 °C.

## 2.2. Determination of Gene Expression in Mammary and Pituitary Glands Using RT-qPCR

### 2.2.1. RNA Extraction and Purification

The entire pituitary gland (weight 19 ± 5 mg) or 28 ± 8 mg (mean ± standard deviation (SD)) of the inguinal mammary gland were placed into 2 mL Eppendorf tubes with glass beads in carbonic ice, and were then crushed using Precellys Evolution Homogenizer® (Ozyme, Saint-Cyr l'Ecole, France) at 4000 rpm for 1 min. Total RNA was isolated from each tissue sample and purified following the NucleoSpin® RNA/Protein (Marchery-Nagel, Düren, Germany) extraction/clean up protocol. RNA was eluted in 60 or 100 µL of RNase-free water. RNA concentration was measured at 260 nm and sample quality determined by ratios of 260/280 and 260/230 nm (Supplementary Table S1) using a NanoVue Plus® (GE Healthcare, Little Chalfont, UK). As both ratios were in the range of 2.0-2.3, this guaranteed of sufficient pure RNA with low contamination by proteins and extraction solvents [31].

### 2.2.2. Reverse Transcription (RT)

Complementary DNA (cDNA) was prepared from 1.5 µg of purified total RNA using SuperScript® IV reverse transcriptase (Invitrogen, Villebon-Sur-Yvette, France) according to the manufacturer's instructions. Moreover, some samples were randomly chosen to perform negative RT (RT⁻) by replacing SuperScript® IV with water, to verify the absence of non-specific cDNA amplification in samples. Reverse transcription was performed with a SimplyAmp® thermal cycler (Applied Biosystem, Villebon-Sur-Yvette, France) with thermal protocol of 23 °C for 10 min, 50 °C for 10 min and 80 °C for 10 min. Finally, the cDNA obtained were diluted to 1:40 before quantitative analysis.
### 2.2.3. Primers Design and Validation

A total of 54 genes linked to milk constituents secretion and its regulation and lactogenic hormones secretion were studied [19,32,33] (Appendix A, Table A1). Primer pairs used for qPCR were designed based on the sequence obtained from the online gene section of National Center for Biotechnology Information (NCBI) (https://www.ncbi.nlm.nih.gov/gene/) and using open-source PerlPrimer v1.1.21 software (http://perlprimer.sourceforge.net.) [34] according to the following parameters: primer melting temperature ranging from 58 to 62 °C with a maximal difference of 2 °C, amplicon size ranging from 100 to 150 bp, primer length ranging from 20 to 24 bp, span intron/exon boundary and overlap intron/exon boundary by 7 bases to avoid the possibility of genomic DNA contamination. The selected primers had, if possible, no expandable dimer and a total dimerization with the lowest possible force (dG < 3). Then, to check for their selectivity and the absence of potential genomic DNA amplification, selected primers were blasted using Primer-BLAST of NCBI (https://www.ncbi.nlm.nih.gov/tools/primer-blast/index.cgi), against *Rattus norvegicus* "Refseq mRNA" and "Refseq representative genome". Finally, all the primers were validated by measuring their efficiency (E) to amplify mammary gland or pituitary gland cDNA and by ensuring they did not amplify RT-. Efficiency was measured by preparing a pool of all samples of 1:40 dilution of cDNA which was then serially diluted 5 times to 1:2. The curve of cycle quantification value (Cq) function of log2(dilution) was linearly fitted and PCR efficiency (E) was calculated according to the following formula:

$$E (\%) = (2^{(-1/\text{slope})} - 1) \times 100$$

Sample cDNA amplifications were compared to the five RT- amplifications by calculating ΔCq(RT-) = mean Cq(RT-)—mean Cq(sample), setting a Cq value of 50 assigned in case of no amplification (Supplementary Table S1). All genes with their official symbol (NCBI), main functions, primer sequences, amplicon length and efficiency are listed in Appendix A Table A1. Overall, the mean amplicon size was 124.7 ± 28.6 (mean ± SD). For mammary and pituitary glands, efficiencies were greater than 90 and ΔCq(RT-) greater than 10, indicating good amplification capacity and negligible genomic DNA amplification.

### 2.2.4. Quantitative Polymerase Chain Reaction (qPCR) Analysis

The qPCR was executed in Hard-Shell® 96-well PCR plates (Bio-Rad, Marne-la-Coquette, France) on a CFX connect™ Real-Time PCR detection instrument (Bio-Rad). Reactions were performed with a total volume of 15 µL, including 5 µL of 1:40 dilution of cDNA, 250 nM of forward and reverse gene-specific primers and iTAQ Universal SYBR® Green Supermix (Bio-Rad). The qPCR conditions were as follows: initial denaturation at 95 °C for 5 min followed by 40 cycles of 95 °C for 10 s, 60 °C for 30 s and 72 °C for 10 s. To check the specificity of PCR products, a melting curve was performed by increasing temperature from 65 to 95 °C with a 0.5 °C increment and 10 s per step.

Relative gene expression was calculated using the 2⁻ΔΔCq method as previously described [35]. For mammary gland, arithmetic mean Cq of three housekeeping genes (HKG): *Actb*, *Uxt* and *Rps9*, was used to calculate ΔCq (CqGene of interest—CqHousekeeping genes). *Actb* is a classically used HKG while *Uxt* and *Rps9* are known to be very stable genes in lactating mammary gland [36,37]. For pituitary gland, arithmetic mean Cq of *Actb*, *Vcp* and *Rps9* was used as HKG to calculate ΔCq. *Vcp* was shown to be very stable in the brain in an RNAseq study of the lab (unpublished data) and the Cq values of the three genes were strongly correlated (r from 0.68 to 0.81). ΔΔCq values (ΔCqTarget group—ΔCqReference group) were calculated using the CTL diet at L12 as a reference group to highlight the gene expression evolution over the lactation period. Values of −ΔΔCq were considered to be representative to the Log2(Fold change) assuming that mean primer efficiency was close to 100% and relatively similar between genes (CV of 4.7% and 5.9% in mammary and pituitary glands, respectively). Finally, fold change expression of the FEN diet was calculated for each lactation period with 2⁻ΔΔCq formula using, for each day, the CTL diet as a reference group.
## 2.3. Enzyme-Linked Immunosorbent Assay (ELISA) of Lactaogenic Hormones in Dam Plasma

Insulin (Rat Insulin ELISA kit®, ALPCO, Salem, NH, USA), IGF-1 (ELISA kit for Insulin Like Growth factor 1, Cloud-Clone Corp., Katy, TX, USA), prolactin (Prolactin rat ELISA, Demeditec Diagnosis GmbH, Kiel, Germany) and leptin (Rat Leptin ELISA Kit, Merck Milliport, Burlington, MA, USA) were measured in dams' plasma, following manufacturers' instructions. Plasma was diluted to 1:10 for the prolactin assay and to 1:100 for the IGF-1 assay. Optical density was read using a microplate reader Varioskan Lux® (ThermoFisher Scientific, Waltham, MA, USA).

## 2.4. Statistical Analysis

To study the global impact of fenugreek on genes responsible for milk synthesis and regulation in the mammary gland, multivariate analysis was performed, using SIMCA P® version 15.02 (Umetrics AB, Kinnelon, NJ, USA). The unsupervised Principal Component Analysis (PCA) and the supervised Partial Least Squares Discriminant Analysis (PLS-DA) were used [38] to cluster the mammary gland samples based on the selected gene expression profiles (−ΔΔCq values). PCA reduced the large set of variables (43 genes) into two principal components (PCA 1 and 2). PLS-DA enhanced the separation between the dietary groups (CTL and FEN) by summarizing the data into a few latent variables that maximized covariance between the response (dietary group) and the predictors (gene expression). The performance of the PLS-DA models was evaluated by the proportion of the total variance of the variables explained by the model (R2(X)), the proportion of the total variance of the response variable (dietary group) explained by the model (R2(Y)), the prediction ability of the model reflected by both a 7-fold cross-validation of the data (Q2) and the reliability of the PLS-DA model based on an analysis of variance (ANOVA) assessment of the cross-validatory (CV) predictive residuals (*p*-value of CV-ANOVA) of the model. The corresponding loading plot was used to determine the genes most responsible for separation in the PLS-DA score plot. Based on the PLS-DA results, genes were plotted according to their importance in separating the dietary groups in PLS-DA score plot and each gene received a value called the variable of importance in the projection (VIP). The variables with the highest VIP values (above 1) were the most powerful group of discriminators involved in the separation of CTL and FEN groups. Finally, in order to highlight diet-specific effect on mammary gland gene expression, lactation time effect was masked by building new PCA and PLS-DA models on a global matrix of normalized variables obtained using mean centering and standardization of variables by lactation point.

For univariate analysis, parametric tests were favored to maximize statistical power, such as the Student's test used for diet comparison. Their validity was admitted as the residual distribution of variables was not significantly different from normal distribution according to the Shapiro-Wilk test and variance was homogenous according to the Brown-Forsythe test. Relative genes' expressions in tissues and plasma hormone concentrations were analyzed using two way ANOVA with diet and period factors followed by Tukey's post-hoc test for period pairwise comparisons (3 levels) and Sidak's post-hoc test for diet pairwise comparisons (2 levels), as suggested by Kim [39]. Values of -ΔΔCq were chosen for statistical analysis of gene expression because distribution of -ΔΔCq values, unlike that of 2⁻⁻ΔΔCq values, was similar to normal distribution that enabled the use of parametric tests [40,41]. To determine the strength of the link between two variables, Pearson's correlation test was used. All univariate tests were performed with the GraphPad prism® software (San Diego, CA, USA), version 6. In all tests, *p* < 0.05 was considered as statistically significant.

# 3. Results

## 3.1. Fenugreek Supplementation Enhances Lactation Performances in a Model of Rat Nursing a Large Litter

Lactation performance was measured in experiments 1 and 2 (Xp1 and Xp2). As no significant interaction effect was found between experiments and diet effects using an ANOVA, data from Xp1 and Xp2 were pooled. Dam water intake was similar in control diet (CTL) and fenugreek-supplemented
diet (FEN) groups throughout the lactation period. Conversely, maternal food intake, pup final weight gain and milk flow from mother to pups all increased significantly, by 10.0%, 7.1% and 15.3%, respectively (Table 1).

<table><caption>**Table 1.** Effect of fenugreek supplementation on dams' lactation performances.</caption><thead><tr><th>Groups</th><th>N</th><th>Dam Water Intake (g/day)</th><th>Dam Food Intake (g/day)</th><th>Pup Total Weight Gain (g)</th><th>Milk Flow (g/day)</th></tr></thead><tbody><tr><td>CTL</td><td>17</td><td>53.93 ± 1.12</td><td>49.71 ± 0.79</td><td>36.19 ± 0.74</td><td>55.90 ± 2.36</td></tr><tr><td>FEN</td><td>16</td><td>55.22 ± 1.11</td><td>54.66 ± 1.00</td><td>38.75 ± 0.85</td><td>64.45 ± 2.35</td></tr><tr><td>p-value</td><td></td><td>0.422</td><td><strong>&lt;0.001</strong></td><td><strong>0.029</strong></td><td><strong>0.016</strong></td></tr></tbody></table>

Values are mean ± SEM (standard error of the mean) and were analyzed with Student's test. Data were obtained from pooled Xp1 and Xp2 experiments. CTL: control diet group and FEN: fenugreek-supplemented diet group throughout the lactation period. The values in bold correspond to significant *p*-values (<0.05).

## 3.2. Different Patterns of Gene Expression in the Mammary and Pituitary Glands and of Plasma Hormone Concentration throughout Lactation and Involution

Multivariate and univariate statistical analysis were performed to investigate the temporal pattern in gene expression in mammary gland from mid-lactation to the beginning of involution.

The genes with the highest expression in mammary gland were Csn2, Wap, Mtc01 and Lalba with ΔCq ≤ −6 and codes for major milk proteins and actors of mitochondrial oxidative phosphorylation (Tables 2 and 3). Plin2 and Fas, two genes involved in lipid synthesis, were also strongly expressed in the mammary gland with ΔCq < −1. Conversely, Cpt1a, Acads, Mtor, Igf1r Esr1, Insr, Lxra, Scd1, Pgm1 and Dgat2 were the least expressed genes in mammary gland over the entire lactation period, with ΔCq > 4 (Tables 2–5). In pituitary gland, Prl and Gh1 were both very intensively expressed, with ΔCq < −6 highlighting the high level of secretion of prolactin and GH during lactation (Table 6). In contrast, Oxt, Trhr and Vipr2 were the less expressed genes with ΔCq > 4.

<table><caption>**Table 2.** Effect of fenugreek supplementation on relative expression of milk constituent synthesis genes in mammary gland along lactation period.</caption><thead><tr><th rowspan="2">Genes</th><th rowspan="2">Period</th><th colspan="2">Log2 (Fold Change)</th><th rowspan="2">Fold Change CTL as Reference</th><th colspan="3">2-Way ANOVA</th></tr><tr><th>CTL at L12 as Reference CTL</th><th>FEN</th><th>Inter</th><th>Period</th><th>Diet</th></tr></thead><tbody><tr><td colspan="8" style="text-align: left; font-weight: bold;">Milk lipid synthesis</td></tr><tr><td rowspan="3">Lpl</td><td>L12 (n = 6)</td><td>0.00 ± 0.27 <sup>a,1</sup></td><td>0.97 ± 0.11 <sup>b,1</sup></td><td>1.95</td><td rowspan="3">0.034</td><td rowspan="3">&lt;0.001</td><td rowspan="3">0.045</td></tr><tr><td>L18 (n = 6)</td><td>0.52 ± 0.13 <sup>1</sup></td><td>0.90 ± 0.43 <sup>1</sup></td><td>1.31</td></tr><tr><td>Inv1 (n = 11)</td><td>-2.63 ± 0.16 <sup>2</sup></td><td>-2.83 ± 0.16 <sup>2</sup></td><td>0.87</td></tr><tr><td rowspan="3">Fasn</td><td>L12</td><td>0.00 ± 0.27 <sup>a,1</sup></td><td>1.10 ± 0.13 <sup>b,1</sup></td><td>2.14</td><td rowspan="3">0.192</td><td rowspan="3">&lt;0.001</td><td rowspan="3">0.049</td></tr><tr><td>L18</td><td>0.21 ± 0.17 <sup>1</sup></td><td>0.25 ± 0.37 <sup>1</sup></td><td>1.03</td></tr><tr><td>Inv1</td><td>-5.10 ± 0.23 <sup>2</sup></td><td>-4.84 ± 0.29 <sup>2</sup></td><td>1.19</td></tr><tr><td rowspan="3">Acaca</td><td>L12</td><td>0.00 ± 0.32 <sup>a,1</sup></td><td>1.08 ± 0.14 <sup>b,1</sup></td><td>2.12</td><td rowspan="3">0.078</td><td rowspan="3">&lt;0.001</td><td rowspan="3">0.012</td></tr><tr><td>L18</td><td>-0.18 ± 0.14 <sup>1</sup></td><td>0.16 ± 0.41 <sup>2</sup></td><td>1.26</td></tr><tr><td>Inv1</td><td>-4.24 ± 0.14 <sup>2</sup></td><td>-4.19 ± 0.18 <sup>3</sup></td><td>1.04</td></tr><tr><td rowspan="3">Acacb</td><td>L12</td><td>0.00 ± 0.16 <sup>1</sup></td><td>0.31 ± 0.08 <sup>1</sup></td><td>1.24</td><td rowspan="3">0.480</td><td rowspan="3">&lt;0.001</td><td rowspan="3">0.965</td></tr><tr><td>L18</td><td>0.52 ± 0.20 <sup>1</sup></td><td>0.43 ± 0.31 <sup>1</sup></td><td>0.94</td></tr><tr><td>Inv1</td><td>-2.12 ± 0.18 <sup>2</sup></td><td>-2.37 ± 0.26 <sup>2</sup></td><td>0.84</td></tr><tr><td rowspan="3">Fabp3</td><td>L12</td><td>0.00 ± 0.27 <sup>a,1</sup></td><td>1.84 ± 0.23 <sup>b,1</sup></td><td>3.59</td><td rowspan="3">0.009</td><td rowspan="3">&lt;0.001</td><td rowspan="3">&lt;0.001</td></tr><tr><td>L18</td><td>0.51 ± 0.16 <sup>1</sup></td><td>0.67 ± 0.38 <sup>2</sup></td><td>1.12</td></tr><tr><td>Inv1</td><td>-4.37 ± 0.16 <sup>2</sup></td><td>-3.95 ± 0.26 <sup>3</sup></td><td>1.34</td></tr><tr><td rowspan="3">Fabp4</td><td>L12</td><td>0.00 ± 0.24 <sup>1</sup></td><td>-0.03 ± 0.17 <sup>1</sup></td><td>0.98</td><td rowspan="3">0.959</td><td rowspan="3">&lt;0.001</td><td rowspan="3">0.759</td></tr><tr><td>L18</td><td>-1.73 ± 0.29 <sup>2</sup></td><td>-1.75 ± 0.48 <sup>2</sup></td><td>0.99</td></tr><tr><td>Inv1</td><td>-0.30 ± 0.16 <sup>1</sup></td><td>-0.44 ± 0.16 <sup>1</sup></td><td>0.91</td></tr><tr><td></td><td>L12</td><td>0.00 ± 0.71 <sup>1</sup></td><td>1.19 ± 0.16 <sup>1,2</sup></td><td>2.28</td><td></td><td></td><td></td></tr></tbody></table>
<table><caption>Table 2. *Cont.*</caption><thead><tr><th rowspan="2">Genes</th><th rowspan="2">Period</th><th colspan="2">Log2 (Fold Change)</th><th rowspan="2">Fold Change CTL as Reference</th><th colspan="3">2-Way ANOVA</th></tr><tr><th>CTL at L12 as Reference CTL</th><th>FEN</th><th>Inter</th><th>Period</th><th>Diet</th></tr></thead><tbody><tr><td rowspan="2">Scd1</td><td>L18</td><td>1.10 ± 0.18</td><td>2.10 ± 0.41 <sup>1</sup></td><td>2.00</td><td rowspan="2">0.399</td><td rowspan="2"><strong>0.013</strong></td><td rowspan="2">0.122</td></tr><tr><td>Inv1</td><td>0.04 ± 0.52</td><td>-0.03 ± 0.56 <sup>2</sup></td><td>0.95</td></tr><tr><td rowspan="3">Dgat1</td><td>L12</td><td>0.00 ± 0.15 <sup>1</sup></td><td>0.21 ± 0.03 <sup>1</sup></td><td>1.16</td><td rowspan="3">0.273</td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3">0.564</td></tr><tr><td>L18</td><td>1.04 ± 0.17 <sup>2</sup></td><td>0.85 ± 0.26 <sup>2</sup></td><td>0.88</td></tr><tr><td>Inv1</td><td>-0.50 ± 0.06 <sup>3</sup></td><td>-0.33 ± 0.08 <sup>3</sup></td><td>1.12</td></tr><tr><td rowspan="3">Dgat2</td><td>L12</td><td>0.00 ± 0.36 <sup>1,2</sup></td><td>-0.74 ± 0.20 <sup>1</sup></td><td>0.60</td><td rowspan="3">0.068</td><td rowspan="3"><strong>0.005</strong></td><td rowspan="3">0.322</td></tr><tr><td>L18</td><td>-0.33 ± 0.10 <sup>1</sup></td><td>-0.25 ± 0.42 <sup>1,2</sup></td><td>1.05</td></tr><tr><td>Inv1</td><td>0.65 ± 0.17 <sup>2</sup></td><td>0.12 ± 0.23 <sup>2</sup></td><td>0.69</td></tr><tr><td rowspan="3">Plin2</td><td>L12</td><td>0.00 ± 0.43 <sup>1</sup></td><td>0.50 ± 0.08 <sup>1</sup></td><td>1.41</td><td rowspan="3">0.267</td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3">0.069</td></tr><tr><td>L18</td><td>1.51 ± 0.16 <sup>2</sup></td><td>1.43 ± 0.26 <sup>2</sup></td><td>0.95</td></tr><tr><td>Inv1</td><td>-0.37 ± 0.13 <sup>1</sup></td><td>0.16 ± 0.16 <sup>1</sup></td><td>1.44</td></tr><tr><th colspan="8" style="text-align: left;">Milk lactose synthesis</th></tr><tr><td rowspan="3">Glut1</td><td>L12</td><td>0.00 ± 0.37 <sup>1</sup></td><td>0.49 ± 0.21 <sup>1</sup></td><td>1.40</td><td rowspan="3">0.310</td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3"><strong>0.031</strong></td></tr><tr><td>L18</td><td>-0.52 ± 0.14 <sup>1</sup></td><td>-0.06 ± 0.29 <sup>1</sup></td><td>1.38</td></tr><tr><td>Inv1</td><td>-2.08 ± 0.07 <sup>2</sup></td><td>-2.04 ± 0.08 <sup>2</sup></td><td>1.03</td></tr><tr><td rowspan="3">Pgm1</td><td>L12</td><td>0.00 ± 0.10 <sup>1</sup></td><td>0.33 ± 0.06 <sup>1</sup></td><td>1.26</td><td rowspan="3"><strong>0.032</strong></td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3">0.591</td></tr><tr><td>L18</td><td>-0.61 ± 0.11 <sup>2</sup></td><td>-0.60 ± 0.18 <sup>2</sup></td><td>1.01</td></tr><tr><td>Inv1</td><td>-0.05 ± 0.05 <sup>1</sup></td><td>-0.26 ± 0.08 <sup>3</sup></td><td>0.86</td></tr><tr><td rowspan="3">Ugp2</td><td>L12</td><td>0.00 ± 0.35 <sup>a,1</sup></td><td>1.03 ± 0.13 <sup>b,1</sup></td><td>2.05</td><td rowspan="3"><strong>0.042</strong></td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3"><strong>0.039</strong></td></tr><tr><td>L18</td><td>-0.67 ± 0.24 <sup>1</sup></td><td>-0.48 ± 0.44 <sup>2</sup></td><td>1.14</td></tr><tr><td>Inv1</td><td>-2.49 ± 0.11 <sup>2</sup></td><td>-2.56 ± 0.06 <sup>3</sup></td><td>0.95</td></tr><tr><td rowspan="3">B4galt1</td><td>L12</td><td>0.00 ± 0.26 <sup>a,1</sup></td><td>1.23 ± 0.13 <sup>b,1</sup></td><td>2.35</td><td rowspan="3"><strong>0.009</strong></td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3"><strong>0.001</strong></td></tr><tr><td>L18</td><td>0.74 ± 0.33 <sup>2</sup></td><td>0.88 ± 0.27 <sup>1</sup></td><td>1.11</td></tr><tr><td>Inv1</td><td>-3.22 ± 0.06 <sup>3</sup></td><td>-3.06 ± 0.09 <sup>2</sup></td><td>1.11</td></tr><tr><td rowspan="3">Lalba</td><td>L12</td><td>0.00 ± 0.58 <sup>a,1</sup></td><td>1.28 ± 0.12 <sup>b,1</sup></td><td>2.42</td><td rowspan="3">0.166</td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3"><strong>0.019</strong></td></tr><tr><td>L18</td><td>1.65 ± 0.15 <sup>2</sup></td><td>1.64 ± 0.30 <sup>1</sup></td><td>0.99</td></tr><tr><td>Inv1</td><td>-1.44 ± 0.28 <sup>3</sup></td><td>-0.86 ± 0.25 <sup>2</sup></td><td>1.49</td></tr><tr><th colspan="8" style="text-align: left;">Milk protein synthesis</th></tr><tr><td rowspan="3">Slc7a5</td><td>L12</td><td>0.00 ± 0.41 <sup>1</sup></td><td>0.83 ± 0.11 <sup>1</sup></td><td>1.78</td><td rowspan="3">0.086</td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3">0.164</td></tr><tr><td>L18</td><td>1.12 ± 0.25 <sup>2</sup></td><td>0.74 ± 0.24 <sup>1</sup></td><td>0.77</td></tr><tr><td>Inv1</td><td>-1.73 ± 0.16 <sup>3</sup></td><td>-1.31 ± 0.26 <sup>2</sup></td><td>1.34</td></tr><tr><td rowspan="3">Csn2</td><td>L12</td><td>0.00 ± 0.60 <sup>1</sup></td><td>1.18 ± 0.06 <sup>b,1</sup></td><td>2.27</td><td rowspan="3">0.104</td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3"><strong>0.013</strong></td></tr><tr><td>L18</td><td>1.48 ± 0.12 <sup>2</sup></td><td>1.36 ± 0.30 <sup>1</sup></td><td>0.92</td></tr><tr><td>Inv1</td><td>-2.13 ± 0.22 <sup>3</sup></td><td>-1.37 ± 0.24 <sup>2</sup></td><td>1.69</td></tr><tr><td rowspan="3">Wap</td><td>L12</td><td>0.00 ± 0.90 <sup>1</sup></td><td>1.47 ± 0.11 <sup>1</sup></td><td>2.77</td><td rowspan="3">0.329</td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3">0.156</td></tr><tr><td>L18</td><td>1.58 ± 0.05 <sup>1</sup></td><td>1.48 ± 0.28 <sup>1</sup></td><td>0.93</td></tr><tr><td>Inv1</td><td>-2.29 ± 0.54 <sup>2</sup></td><td>-1.91 ± 0.42 <sup>2</sup></td><td>1.30</td></tr><tr><th colspan="8" style="text-align: left;">Lactocyte water inflow</th></tr><tr><td rowspan="3">Aqp1</td><td>L12</td><td>0.00 ± 0.45 <sup>1</sup></td><td>0.58 ± 0.11</td><td>1.50</td><td rowspan="3"><strong>0.006</strong></td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3">0.458</td></tr><tr><td>L18</td><td>1.24 ± 0.10 <sup>2</sup></td><td>0.57 ± 0.26</td><td>0.63</td></tr><tr><td>Inv1</td><td>-0.39 ± 0.10 <sup>1</sup></td><td>0.05 ± 0.14</td><td>1.36</td></tr><tr><td rowspan="3">Aqp3</td><td>L12</td><td>0.00 ± 0.38 <sup>1</sup></td><td>0.63 ± 0.48 <sup>1</sup></td><td>1.55</td><td rowspan="3">0.751</td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3">0.240</td></tr><tr><td>L18</td><td>1.64 ± 0.35 <sup>2</sup></td><td>1.82 ± 0.28 <sup>1</sup></td><td>1.13</td></tr><tr><td>Inv1</td><td>-0.91 ± 0.30 <sup>1</sup></td><td>-0.75 ± 0.17 <sup>2</sup></td><td>1.12</td></tr></tbody></table>

Log2 (fold change) corresponded to -ΔΔCq values calculated with mean of 3 housekeeping genes (*Actb*, *Uxt*, *Rps9*) and with CTL at L12 as a reference group. Fold change correspond to 2⁻ΔΔCq values and were calculated each day for FEN diet with CTL of the day as a reference. Log2 (fold change) values were mean ± SEM, *n* = 5–6 at L12 and L18 and 10–11 at Inv1. Log2 (fold change) values were analyzed with 2-way ANOVA (*p*-values of diet, period and interaction effects reported in the last columns) followed by Sidak's *post-hoc* test for diet effect and Tukey's *post-hoc* test for period effect. Different letters and numbers represent significant difference (*p* < 0.05) between diets for each period and between periods for each diet, respectively. Green boxes represent significant overexpression in FEN group compared to CTL group at a given day. CTL: control diet group and FEN: fenugreek-supplemented diet group throughout the lactation period. L12 and L18: post-natal day of lactation 12 and 18, Inv1: first day of involution. The values in bold correspond to significant *p*-values (<0.05).
<table><caption>**Table 3.** Effect of fenugreek supplementation on relative expression of genes involved in energy metabolism and antioxidant activity in mammary gland along lactation period.</caption><thead><tr><th rowspan="2">Genes</th><th rowspan="2">Period</th><th colspan="2">Log2 (Fold Change)</th><th rowspan="2">Fold Change<br>CTL as Reference</th><th colspan="3">2-Way ANOVA</th></tr><tr><th>CTL at L12 as Reference<br>CTL</th><th>FEN</th><th>Inter</th><th>Period</th><th>Diet</th></tr></thead><tbody><tr><td colspan="8" style="text-align: left; font-weight: bold;">Lipid β-oxidation</td></tr><tr><td rowspan="3">Cpt1a</td><td>L12 (n = 6)</td><td>0.00 ± 0.75 <sup>a,1</sup></td><td>1.94 ± 0.64 <sup>b,1</sup></td><td>3.84</td><td rowspan="3">0.021</td><td rowspan="3">&lt;0.001</td><td rowspan="3">0.007</td></tr><tr><td>L18 (n = 6)</td><td>-0.26 ± 0.53 <sup>1</sup></td><td>0.29 ± 0.46 <sup>2</sup></td><td>1.46</td></tr><tr><td>Inv1 (n = 11)</td><td>3.28 ± 0.05 <sup>2</sup></td><td>3.22 ± 0.08 <sup>3</sup></td><td>0.96</td></tr><tr><td rowspan="3">Acads</td><td>L12</td><td>0.00 ± 0.52 <sup>a,1</sup></td><td>1.47 ± 0.40 <sup>b,1</sup></td><td>2.77</td><td rowspan="3">0.007</td><td rowspan="3">&lt;0.001</td><td rowspan="3">0.002</td></tr><tr><td>L18</td><td>0.19 ± 0.08 <sup>1</sup></td><td>0.39 ± 0.13 <sup>2</sup></td><td>1.15</td></tr><tr><td>Inv1</td><td>1.28 ± 0.06 <sup>2</sup></td><td>1.41 ± 0.07 <sup>1</sup></td><td>1.09</td></tr><tr><td colspan="8" style="text-align: left; font-weight: bold;">Krebs cycle initiation</td></tr><tr><td rowspan="3">Pdha1</td><td>L12</td><td>0.00 ± 0.53 <sup>1</sup></td><td>0.53 ± 0.23 <sup>1</sup></td><td>1.45</td><td rowspan="3">0.465</td><td rowspan="3">&lt;0.001</td><td rowspan="3">0.071</td></tr><tr><td>L18</td><td>0.32 ± 0.08 <sup>1</sup></td><td>0.59 ± 0.22 <sup>1</sup></td><td>1.21</td></tr><tr><td>Inv1</td><td>-1.46 ± 0.05 <sup>2</sup></td><td>-1.39 ± 0.06 <sup>2</sup></td><td>1.05</td></tr><tr><td rowspan="3">Cs</td><td>L12</td><td>0.00 ± 0.40 <sup>1</sup></td><td>0.70 ± 0.25 <sup>1</sup></td><td>1.62</td><td rowspan="3">0.246</td><td rowspan="3">&lt;0.001</td><td rowspan="3">0.013</td></tr><tr><td>L18</td><td>0.10 ± 0.07 <sup>1</sup></td><td>0.32 ± 0.20 <sup>1</sup></td><td>1.16</td></tr><tr><td>Inv1</td><td>-1.20 ± 0.04 <sup>2</sup></td><td>-1.03 ± 0.07 <sup>2</sup></td><td>1.12</td></tr><tr><td colspan="8" style="text-align: left; font-weight: bold;">Oxidative phosphorylation</td></tr><tr><td rowspan="3">Mtco1</td><td>L12</td><td>0.00 ± 0.53 <sup>1</sup></td><td>0.74 ± 0.25 <sup>1</sup></td><td>1.68</td><td rowspan="3">0.148</td><td rowspan="3">&lt;0.001</td><td rowspan="3">0.099</td></tr><tr><td>L18</td><td>0.52 ± 0.17 <sup>1</sup></td><td>0.35 ± 0.18 <sup>1</sup></td><td>0.89</td></tr><tr><td>Inv1</td><td>-0.90 ± 0.07 <sup>2</sup></td><td>-0.61 ± 0.12 <sup>2</sup></td><td>1.22</td></tr><tr><td rowspan="3">Atp5f1a</td><td>L12</td><td>0.00 ± 0.50</td><td>0.55 ± 0.23 <sup>1</sup></td><td>1.47</td><td rowspan="3">0.251</td><td rowspan="3">&lt;0.001</td><td rowspan="3">0.073</td></tr><tr><td>L18</td><td>-0.24 ± 0.15</td><td>0.07 ± 0.18 <sup>1</sup></td><td>1.24</td></tr><tr><td>Inv1</td><td>-0.61 ± 0.03</td><td>-0.64 ± 0.07 <sup>2</sup></td><td>0.98</td></tr><tr><td colspan="8" style="text-align: left; font-weight: bold;">Antioxidant enzymes</td></tr><tr><td rowspan="3">Sod1</td><td>L12</td><td>0.00 ± 0.12 <sup>1</sup></td><td>-0.24 ± 0.08</td><td>0.85</td><td rowspan="3">0.285</td><td rowspan="3">0.012</td><td rowspan="3">0.142</td></tr><tr><td>L18</td><td>-0.55 ± 0.14 <sup>2</sup></td><td>-0.48 ± 0.18</td><td>1.05</td></tr><tr><td>Inv1</td><td>-0.22 ± 0.08 <sup>1,2</sup></td><td>-0.48 ± 0.08</td><td>0.83</td></tr><tr><td rowspan="3">Cat</td><td>L12</td><td>0.00 ± 0.18</td><td>-0.18 ± 0.10</td><td>0.88</td><td rowspan="3">0.458</td><td rowspan="3">0.021</td><td rowspan="3">0.726</td></tr><tr><td>L18</td><td>-0.85 ± 0.13</td><td>-0.85 ± 0.21</td><td>1.00</td></tr><tr><td>Inv1</td><td>-0.51 ± 0.28</td><td>-0.11 ± 0.23</td><td>1.32</td></tr><tr><td rowspan="3">Gpx1</td><td>L12</td><td>0.00 ± 0.18</td><td>0.20 ± 0.13 <sup>1</sup></td><td>1.14</td><td rowspan="3">0.310</td><td rowspan="3">&lt;0.001</td><td rowspan="3">0.680</td></tr><tr><td>L18</td><td>0.14 ± 0.13</td><td>0.38 ± 0.23 <sup>1</sup></td><td>1.18</td></tr><tr><td>Inv1</td><td>-0.47 ± 0.17</td><td>-0.71 ± 0.17 <sup>2</sup></td><td>0.85</td></tr></tbody></table>

Log2 (fold change) corresponded to $-\Delta\Delta C_q$ values calculated with mean of 3 housekeeping genes (*Actb*, *Uxt*, *Rps9*) and with CTL at L12 as a reference group. Fold change correspond to $2^{-\Delta\Delta C_q}$ values and are calculated each day for FEN diet with CTL of the day as a reference. Log2 (fold change) values are mean ± SEM, *n* = 5–6 at L12 and L18 and 10–11 at Inv1. Log2 (fold change) values were analyzed with 2-way ANOVA (*p*-values of diet, period and interaction effects reported in the last columns) followed by Sidak's *post-hoc* test for diet effect and Tukey's *post-hoc* test for period effect. Different letters and numbers represent significant difference (*p* < 0.05) between diets for each period and between periods for each diet, respectively. Green boxes represent significant overexpression in FEN group compared to CTL group at a given day. CTL: control diet group and FEN: fenugreek-supplemented diet group throughout the lactation period. L12 and L18: post-natal day of lactation 12 and 18, Inv1: first day of involution. The values in bold correspond to significant *p*-values (<0.05).

<table><caption>**Table 4.** Effect of fenugreek on relative expression of genes related to milk synthesis regulatory factors in the mammary gland along lactation period.</caption><thead><tr><th rowspan="2">Genes</th><th rowspan="2">Period</th><th colspan="2">Log2 (Fold Change)</th><th rowspan="2">Fold Change<br>CTL as Reference</th><th colspan="3">2-Way ANOVA</th></tr><tr><th>CTL at L12 as Reference<br>CTL</th><th>FEN</th><th>Inter</th><th>Period</th><th>Diet</th></tr></thead><tbody><tr><td colspan="8" style="text-align: left; font-weight: bold;">Lipid metabolism regulatory factors</td></tr><tr><td rowspan="3">Srebf1</td><td>L12 (n = 6)</td><td>0.00 ± 0.37 <sup>1</sup></td><td>0.40 ± 0.52 <sup>1</sup></td><td>1.32</td><td rowspan="3">0.605</td><td rowspan="3">&lt;0.001</td><td rowspan="3">0.543</td></tr><tr><td>L18 (n = 6)</td><td>-2.11 ± 0.29 <sup>2</sup></td><td>-1.76 ± 0.62 <sup>2</sup></td><td>1.27</td></tr><tr><td>Inv1 (n = 11)</td><td>-3.03 ± 0.18 <sup>2</sup></td><td>-3.24 ± 0.21 <sup>3</sup></td><td>0.87</td></tr></tbody></table>
<table><caption>**Table 4.** *Cont.*</caption><thead><tr><th rowspan="2">Genes</th><th rowspan="2">Period</th><th colspan="2">Log2 (Fold Change)</th><th rowspan="2">Fold Change<br>CTL as Reference</th><th colspan="3">2-Way ANOVA</th></tr><tr><th>CTL at L12 as Reference<br>CTL</th><th>FEN</th><th>Inter</th><th>Period</th><th>Diet</th></tr></thead><tbody><tr><td rowspan="3">Pparg</td><td>L12</td><td>0.00 ± 0.46 <sup>1</sup></td><td>0.77 ± 0.24 <sup>1</sup></td><td>1.71</td><td rowspan="3">0.097</td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3">0.432</td></tr><tr><td>L18</td><td>-0.07 ± 0.28 <sup>1</sup></td><td>-0.12 ± 0.39 <sup>2</sup></td><td>0.97</td></tr><tr><td>Inv1</td><td>-0.93 ± 0.10 <sup>2</sup></td><td>-1.19 ± 0.10 <sup>3</sup></td><td>0.84</td></tr><tr><td rowspan="3">Spot14</td><td>L12</td><td>0.00 ± 0.49 <sup>1</sup></td><td>0.55 ± 0.38 <sup>1</sup></td><td>1.46</td><td rowspan="3">0.698</td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3">0.324</td></tr><tr><td>L18</td><td>-0.03 ± 0.40 <sup>1</sup></td><td>0.23 ± 0.39 <sup>1</sup></td><td>1.19</td></tr><tr><td>Inv1</td><td>-1.28 ± 24 <sup>2</sup></td><td>-1.28 ± 0.18 <sup>2</sup></td><td>1.00</td></tr><tr><td rowspan="3">Lxra</td><td>L12</td><td>0.00 ± 0.61 <sup>a,1</sup></td><td>1.24 ± 0.37 <sup>b</sup></td><td>2.38</td><td rowspan="3"><strong>0.024</strong></td><td rowspan="3">0.053</td><td rowspan="3">0.086</td></tr><tr><td>L18</td><td>0.52 ± 0.17 <sup>1-2</sup></td><td>0.44 ± 0.27</td><td>0.95</td></tr><tr><td>Inv1</td><td>1.08 ± 0.07 <sup>2</sup></td><td>1.00 ± 0.13</td><td>0.95</td></tr><tr><th colspan="8" style="text-align: left; background-color: #E0E0E0;">Protein synthesis regulatory factors</th></tr><tr><td rowspan="3">Akt1</td><td>L12</td><td>0.00 ± 0.48 <sup>1</sup></td><td>0.59 ± 0.24 <sup>1</sup></td><td>1.51</td><td rowspan="3">0.209</td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3">0.128</td></tr><tr><td>L18</td><td>-0.02 ± 0.08 <sup>1</sup></td><td>-0.10 ± 0.18 <sup>2</sup></td><td>0.94</td></tr><tr><td>Inv1</td><td>-0.81 ± 0.05 <sup>2</sup></td><td>-0.65 ± 0.05 <sup>2</sup></td><td>1.05</td></tr><tr><td rowspan="3">Mtor</td><td>L12</td><td>0.00 ± 0.56 <sup>a</sup></td><td>0.88 ± 0.23 <sup>b,1</sup></td><td>1.84</td><td rowspan="3"><strong>0.040</strong></td><td rowspan="3">0.189</td><td rowspan="3">0.223</td></tr><tr><td>L18</td><td>0.34 ± 0.13</td><td>0.09 ± 0.27 <sup>2</sup></td><td>0.84</td></tr><tr><td>Inv1</td><td>0.06 ± 0.05</td><td>0.06 ± 0.04 <sup>2</sup></td><td>1.00</td></tr><tr><td rowspan="3">Stat5</td><td>L12</td><td>0.00 ± 0.51</td><td>0.78 ± 0.34 <sup>1</sup></td><td>1.71</td><td rowspan="3">0.164</td><td rowspan="3"><strong>0.004</strong></td><td rowspan="3">0.212</td></tr><tr><td>L18</td><td>0.17 ± 0.13</td><td>0.09 ± 0.26 <sup>1,2</sup></td><td>0.94</td></tr><tr><td>Inv1</td><td>-0.36 ± 0.08</td><td>-0.35 ± 0.13 <sup>2</sup></td><td>1.01</td></tr></tbody></table>

Log2 (fold change) corresponded to -ΔΔCq values calculated with mean of 3 housekeeping genes (*Actb*, *Uxt*, *Rps9*) and with CTL at L12 as a reference group. Fold change correspond to 2⁻ΔΔCq values and are calculated each day for FEN diet with CTL of the day as a reference. Log2 (fold change) values are mean ± SEM, *n* = 5–6 at L12 and L18 and 10–11 at Inv1. Log2 (fold change) values were analyzed with 2-way ANOVA (*p*-values of diet, period and interaction effects reported in the last columns) followed by Sidak's post-hoc test for diet effect and Tukey's post-hoc test for period effect. Different letters and numbers represent significant difference (*p* < 0.05) between diets for each period and between periods for each diet, respectively. Green boxes represent significant overexpression in FEN group at a given day. CTL: control diet group and FEN: fenugreek-supplemented diet group throughout the lactation period. L12 and L18: post-natal day of lactation 12 and 18, Inv1: first day of involution. The values in bold correspond to significant *p*-values (<0.05).

Following the unsupervised Principal Component Analysis (PCA) on the un-normalized complete set of 43 genes, the PCA score plot, on components 1 and 2, revealed natural longitudinal classification with differences corresponding to the three lactation time points and a clear-cut separation between involution (Inv1) and lactation period (L12 and L18) (Figure 2a). The PCA loading scatter plot (Figure 2b) reflected relationships between variables in this model, deciphering three clusters that we labelled "kinetics 1, 2 and 3". In agreement with this longitudinal pattern, the univariate statistical analysis showed that expression of all genes evolves significantly between the middle and end of lactation, except for *Ghr* and *Mtor* (Tables 4 and 5).

Most of the genes studied exhibited a significant decrease in their expression between L18 and Inv1 ranging from 1.46-fold (*p* = 0.010) for *Atp5f1a* to 36.91-fold (*p* < 0.001) for *Fasn* and belonged to "kinetic 1" and "kinetic 2" patterns.

In the "kinetic 1" pattern, genes similarly expressed or only slightly modulated between L12 and L18, suggesting they had reached peak expression in this interval (Figure 2c). This pattern was represented mostly by genes coding for regulation of milk nutrient synthesis, energy metabolism, and galactose uptakes (*Stat5*, *Akt1*, *Spot14*, *Pparg*, *Prlr*, *Pdha1*, *Acacb*, *Mtco1*, *Cs*, *Ugp2* and *Atp5f1a*). All these genes were highly correlated with each other, as validated by Pearson's test (Figure 2f and Supplementary Figure S1) and mainly with *Prlr* (r ≥ 0.8), considered as the barycenter of the "kinetic 1" cluster in the PCA-loading scatter plot (Figure 2b). If not in the "kinetic 1" ellipse in the PCA (Figure 2b) genes related to de novo fatty acid synthesis and glucose uptake such as *Lpl*, *Fasn*, *Acaca*, *Fabp3* and *Glut1* were also strongly correlated with *Prlr* (r ≥ 0.7) (Supplementary Figure S2).
<table><caption>**Table 5.** Effect of fenugreek on relative expression of genes of lactogenic hormone receptors in mammary gland along lactation period.</caption><thead><tr><th rowspan="2">Genes</th><th rowspan="2">Period</th><th colspan="2">Log2 (Fold Change)</th><th rowspan="2">Fold Change<br>CTL as Reference</th><th colspan="3">2-Way ANOVA</th></tr><tr><th>CTL at L12 as Reference<br>CTL</th><th>FEN</th><th>Inter</th><th>Period</th><th>Diet</th></tr></thead><tbody><tr><td colspan="9" style="text-align: left; font-weight: bold;">Lactogenic hormone receptors</td></tr><tr><td rowspan="3">Prlr</td><td>L12 (n = 6)</td><td>0.00 ± 0.49 <sup>1</sup></td><td>0.40 ± 0.26 <sup>1</sup></td><td>1.32</td><td rowspan="3">0.499</td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3">0.642</td></tr><tr><td>L18 (n = 6)</td><td>-0.32 ± 0.30 <sup>1</sup></td><td>-0.11 ± 0.49 <sup>1</sup></td><td>1.15</td></tr><tr><td>Inv1 (n = 11)</td><td>-1.84 ± 0.16 <sup>2</sup></td><td>-2.09 ± 0.22 <sup>2</sup></td><td>0.84</td></tr><tr><td rowspan="3">Insr</td><td>L12</td><td>0.00 ± 0.64 <sup>1</sup></td><td>0.84 ± 0.32 <sup>1</sup></td><td>1.79</td><td rowspan="3">0.283</td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3">0.082</td></tr><tr><td>L18</td><td>2.29 ± 0.19 <sup>2</sup></td><td>2.35 ± 0.18 <sup>2</sup></td><td>1.04</td></tr><tr><td>Inv1</td><td>-0.85 ± 0.09 <sup>3</sup></td><td>-0.67 ± 0.11 <sup>3</sup></td><td>1.13</td></tr><tr><td rowspan="3">Oxtr</td><td>L12</td><td>0.00 ± 0.62 <sup>1</sup></td><td>0.79 ± 0.69 <sup>1</sup></td><td>1.72</td><td rowspan="3">0.649</td><td rowspan="3"><strong>0.025</strong></td><td rowspan="3">0.557</td></tr><tr><td>L18</td><td>-0.52 ± 0.37 <sup>1</sup></td><td>-0.58 ± 0.87 <sup>1</sup></td><td>0.96</td></tr><tr><td>Inv1</td><td>0.75 ± 0.22 <sup>1</sup></td><td>0.73 ± 0.26 <sup>1</sup></td><td>0.99</td></tr><tr><td rowspan="3">Ghr</td><td>L12</td><td>0.00 ± 0.56</td><td>0.89 ± 0.43</td><td>1.86</td><td rowspan="3">0.286</td><td rowspan="3">0.405</td><td rowspan="3">0.097</td></tr><tr><td>L18</td><td>-0.14 ± 0.23</td><td>0.34 ± 0.43</td><td>1.40</td></tr><tr><td>Inv1</td><td>0.07 ± 0.18</td><td>0.01 ± 0.18</td><td>0.96</td></tr><tr><td rowspan="3">Igf1r</td><td>L12</td><td>0.00 ± 0.59 <sup>a,1</sup></td><td>1.37 ± 0.52 <sup>b,1,2</sup></td><td>2.60</td><td rowspan="3">0.057</td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3"><strong>0.033</strong></td></tr><tr><td>L18</td><td>0.89 ± 0.11 <sup>1</sup></td><td>0.96 ± 0.38 <sup>1</sup></td><td>1.05</td></tr><tr><td>Inv1</td><td>1.87 ± 0.19 <sup>2</sup></td><td>1.97 ± 0.10 <sup>2</sup></td><td>1.07</td></tr><tr><td rowspan="3">Esr1</td><td>L12</td><td>0.00 ± 0.54 <sup>1</sup></td><td>0.34 ± 0.18 <sup>1</sup></td><td>1.27</td><td rowspan="3">0.204</td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3">0.547</td></tr><tr><td>L18</td><td>1.57 ± 0.14 <sup>2</sup></td><td>1.09 ± 0.32 <sup>1</sup></td><td>0.72</td></tr><tr><td>Inv1</td><td>-2.02 ± 0.24 <sup>3</sup></td><td>-1.44 ± 0.17 <sup>2</sup></td><td>1.49</td></tr></tbody></table>

Log2 (fold change) corresponded to -ΔΔCq values calculated with mean of 3 housekeeping genes (*Actb*, *Uxt*, *Rps9*) and with CTL at L12 as a reference group. Fold change correspond to 2⁻ΔΔCq values and are calculated each day for FEN diet with CTL of the day as a reference. Log2 (fold change) values are mean ± SEM, *n* = 5–6 at L12 and L18 and 10–11 at Inv1. Log2 (fold change) values were analyzed with 2-way ANOVA (*p*-values of diet, period and interaction effects reported in the last columns) followed by Sidak's *post-hoc* test for diet effect and Tukey's *post-hoc* test for period effect. Different letters and numbers represent significant difference (*p* < 0.05) between diets for each period and between periods for each diet, respectively. Green boxes represent significant overexpression in FEN group compared to CTL group at a given day. CTL: control diet group and FEN: fenugreek-supplemented diet group throughout the lactation period. L12 and L18: post-natal day of lactation 12 and 18, Inv1: first day of involution. The values in bold correspond to significant *p*-values (<0.05).

In "kinetic 2", overall gene expression rose from 1.70-fold (*p* = 0.042) for *Csn2* to 2.62-fold for *Aqp3* (*p* = 0.001) between L12 and L18, suggesting that peak expression of these genes occurred closer to L18 than L12 (Figure 2d). This cluster mainly included genes related to milk constituent secretion (*Slc7a5*, *Wap*, *Csn2*, *B4galt1*, *Lalba*, *Plin 2*, *Aqp3*, *InsR*, *Esr* and *Dgat1*) as illustrated in Figure 2b,d,g, with *Insr*, closest to the barycenter of the "kinetic 2" cluster (Figure 2b), and strongly correlated with all these genes (Figure 2g). Otherwise, *Esr1* was probably linked to protein expression as its expression was very strongly correlated with all main milk protein genes (*r* ≥ 0.85).

In "kinetic 3", contrary to other clusters, genes were significantly overexpressed at Inv1 compared to L18 (Figure 2e), with a rise ranging from 1.37-fold (*p* < 0.001) for *Pgm1* to 9.44-fold (*p* < 0.001) for *Cpt1a*. Several of these genes involved in fatty acid oxidation pathway (*Cpt1a*, *Acads*, *Lxra* and *Igf1r*) presented high correlation between each other (Figure 2h).

In the pituitary gland, the peak of expression was around L12 for Drd2, between L12 and L18 for Ghrhr and Prl, around L18 for Vipr2 and Trhr and around Inv1 for *Esr1* (Table 6).
<table><caption>**Table 6.** Effect of fenugreek on relative expression of lactogenic hormone-related genes in pituitary gland along lactation period.</caption><thead><tr><th rowspan="2">Genes</th><th rowspan="2">Period</th><th colspan="2">Log2 (Fold Change)</th><th rowspan="2">Fold Change<br>CTL as<br>Reference</th><th colspan="3">2-Way ANOVA</th></tr><tr><th>CTL at L12 as Reference<br>CTL</th><th>FEN</th><th>Inter</th><th>Period</th><th>Diet</th></tr></thead><tbody><tr><td colspan="9" style="text-align: left; font-weight: bold;">Pituitary lactogenic hormones</td></tr><tr><td rowspan="3">Prl</td><td>L12 (n = 5)</td><td>0.00 ± 0.06 <sup>1</sup></td><td>-0.15 ± 0.20 <sup>1</sup></td><td>0.90</td><td rowspan="3">0.291</td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3">0.652</td></tr><tr><td>L18 (n = 5)</td><td>-0.14 ± 0.15 <sup>1</sup></td><td>-0.20 ± 0.11 <sup>1</sup></td><td>0.95</td></tr><tr><td>Inv1 (n = 10)</td><td>-1.18 ± 0.06 <sup>2</sup></td><td>-0.93 ± 0.14 <sup>2</sup></td><td>1.19</td></tr><tr><td rowspan="3">Gh1</td><td>L12</td><td>0.00 ± 0.25</td><td>-0.34 ± 0.05 <sup>1</sup></td><td>0.79</td><td rowspan="3">0.143</td><td rowspan="3">0.110</td><td rowspan="3">0.925</td></tr><tr><td>L18</td><td>0.02 ± 0.19</td><td>0.21 ± 0.19 <sup>2</sup></td><td>1.14</td></tr><tr><td>Inv1</td><td>-0.24 ± 0.09</td><td>-0.06 ± 0.10 <sup>1-2</sup></td><td>1.14</td></tr><tr><td rowspan="3">Oxt</td><td>L12</td><td>0.00 ± 0.48</td><td>0.70 ± 0.12</td><td>1.62</td><td rowspan="3"><strong>0.050</strong></td><td rowspan="3">0.469</td><td rowspan="3"><strong>0.007</strong></td></tr><tr><td>L18</td><td>0.16 ± 0.44 <sup>a</sup></td><td>1.08 ± 0.10 <sup>b</sup></td><td>1.88</td></tr><tr><td>Inv1</td><td>0.42 ± 0.08</td><td>0.38 ± 0.12</td><td>0.97</td></tr><tr><td colspan="9" style="text-align: left; font-weight: bold;">Receptors of lactogenic hormone activators</td></tr><tr><td rowspan="3">Drd2</td><td>L12</td><td>0.00 ± 0.05 <sup>1</sup></td><td>-0.13 ± 0.06 <sup>1</sup></td><td>0.91</td><td rowspan="3">0.040</td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3">0.984</td></tr><tr><td>L18</td><td>-0.39 ± 0.04 <sup>2</sup></td><td>-0.45 ± 0.07 <sup>2</sup></td><td>0.96</td></tr><tr><td>Inv1</td><td>-0.46 ± 0.06 <sup>a,2</sup></td><td>-0.27 ± 0.06 <sup>b,1,2</sup></td><td>1.14</td></tr><tr><td rowspan="3">Trhr</td><td>L12</td><td>0.00 ± 0.16 <sup>1</sup></td><td>0.42 ± 0.21 <sup>1</sup></td><td>1.34</td><td rowspan="3">0.551</td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3"><strong>0.035</strong></td></tr><tr><td>L18</td><td>0.78 ± 0.10 <sup>2</sup></td><td>0.99 ± 0.10 <sup>2</sup></td><td>1.15</td></tr><tr><td>Inv1</td><td>0.55 ± 0.13 <sup>2</sup></td><td>0.67 ± 0.07 <sup>1-2</sup></td><td>1.09</td></tr><tr><td rowspan="3">Vipr2</td><td>L12</td><td>0.00 ± 0.25</td><td>0.35 ± 0.11 <sup>1</sup></td><td>1.27</td><td rowspan="3">0.375</td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3">0.061</td></tr><tr><td>L18</td><td>0.65 ± 0.16</td><td>1.07 ± 0.07 <sup>2</sup></td><td>1.34</td></tr><tr><td>Inv1</td><td>0.22 ± 0.15</td><td>0.25 ± 0.14 <sup>1</sup></td><td>1.02</td></tr><tr><td rowspan="3">Ghrhr</td><td>L12</td><td>0.00 ± 0.16 <sup>1</sup></td><td>-0.65 ± 0.44 <sup>1</sup></td><td>0.64</td><td rowspan="3">0.069</td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3">0.382</td></tr><tr><td>L18</td><td>0.46 ± 0.33 <sup>1</sup></td><td>1.20 ± 0.21 <sup>2</sup></td><td>1.68</td></tr><tr><td>Inv1</td><td>-2.31 ± 0.15 <sup>2</sup></td><td>-1.77 ± 0.27 <sup>3</sup></td><td>1.45</td></tr><tr><td rowspan="3">Ar</td><td>L12</td><td>0.00 ± 0.08</td><td>0.11 ± 0.09</td><td>1.08</td><td rowspan="3">0.420</td><td rowspan="3">0.152</td><td rowspan="3">0.565</td></tr><tr><td>L18</td><td>0.27 ± 0.08</td><td>0.18 ± 0.09</td><td>0.94</td></tr><tr><td>Inv1</td><td>0.11 ± 0.04</td><td>0.19 ± 0.07</td><td>1.06</td></tr><tr><td rowspan="3">Esr1</td><td>L12</td><td>0.00 ± 0.14 <sup>1</sup></td><td>0.26 ± 0.35 <sup>1</sup></td><td>1.20</td><td rowspan="3">0.508</td><td rowspan="3"><strong>&lt;0.001</strong></td><td rowspan="3">0.823</td></tr><tr><td>L18</td><td>0.56 ± 0.12 <sup>1-2</sup></td><td>0.69 ± 0.09 <sup>1-2</sup></td><td>1.09</td></tr><tr><td>Inv1</td><td>1.35 ± 0.16 <sup>2</sup></td><td>1.10 ± 0.27 <sup>2</sup></td><td>0.84</td></tr><tr><td colspan="9" style="text-align: left; font-weight: bold;">Estrogenic effect biomarker</td></tr><tr><td rowspan="3">C3</td><td>L12</td><td>0.00 ± 0.29</td><td>-0.09 ± 0.35 <sup>1</sup></td><td>0.94</td><td rowspan="3"><strong>0.017</strong></td><td rowspan="3"><strong>0.006</strong></td><td rowspan="3"><strong>0.006</strong></td></tr><tr><td>L18</td><td>0.24 ± 0.18 <sup>a</sup></td><td>1.31 ± 0.15 <sup>b,2</sup></td><td>2.10</td></tr><tr><td>Inv1</td><td>0.61 ± 0.26</td><td>1.00 ± 0.13 <sup>2</sup></td><td>1.31</td></tr></tbody></table>

Log2 (fold change) corresponded to -ΔΔCq values calculated with mean of 3 housekeeping genes (*Actb*, *Uxt*, *Rps9*) and with CTL at L12 as a reference group. Fold change correspond to 2⁻⁻ΔΔCq values and are calculated each day for FEN diet with CTL of the day as a reference. Log2 (fold change) values are mean ± SEM, *n* = 4-5 at L12 and L18 and 9-10 at Inv1 (for ESR1 *n* = 4 for CTL at L12). Log2 (fold change) values were analyzed with 2-way ANOVA (*p*-values of diet, period and interaction effects reported in the last columns) followed by Sidak's *post-hoc* test for diet effect and Tukey's *post-hoc* test for period effect. Different letters and numbers represent significant difference (*p* < 0.05) between diets for each period and between periods for each diet, respectively. Green boxes represent significant overexpression in FEN group compared to CTL group at a given day. CTL: control diet group and FEN: fenugreek-supplemented diet group throughout the lactation period. L12 and L18: post-natal day of lactation 12 and 18, Inv1: first day of involution. The values in bold correspond to significant *p*-values (<0.05).

Insulin, IGF-1 and leptin maternal plasma concentrations exhibited the same time course (Figure 3) with a significant decline between L12 and L18 (-65.9%, *p* = 0.030; -48.5%, *p* < 0.001; -79.7%, *p* < 0.001, respectively) followed by a sharp rise between L18 and Inv1 (+392.5%, *p* < 0.001; +156.9%, *p* < 0.001; +325.0%, *p* < 0.001, respectively), regardless of maternal diet. Plasma prolactin concentration decreased by 64.0% (*p* = 0.012) between L12 and L18. Finally, plasma estrogen concentration rose 37.1% (*p* = 0.029) between L18 and Inv1.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Temporal pattern of mammary gland gene expression in CTL and FEN groups through lactation and until involution. (a) Score plot on the first two dimensions of the PCA, accounting for 71% of the total variance for 43 un-normalized data (genes), (b) PCA loading plot presented clusters of variables highly associated to L12- and L18-groups (kinetic 1 and 2 clusters) or to Inv1 group (kinetic 3 cluster). Temporal pattern expression of representative genes of kinetic 1 (c), kinetic 2 (d) and kinetic 3 (e) clusters. Values were mean ± SEM and were analyzed with 2-way ANOVA followed by Tukey's post-hoc tests. #: *p* < 0.05, ##: *p* < 0.01, ###: *p* < 0.001 for comparison with previous period. Correlations between gene expression in cluster 1 (f), 2 (g) and 3 (h) were analyzed with Pearson tests. The shade of green corresponds to the strength of the correlation between the variables with a light green for a weak correlation and a dark green for a strong correlation. CTL: control diet group and FEN: fenugreek-supplemented diet group throughout the lactation period. L12 and L18: post-natal day of lactation 12 and 18, Inv1: first day of involution</figcaption></figure>

3.3. *Fenugreek Supplementation Impacts Mammary Gland Metabolic Pathways Differently during Lactation and at Involution*

A supervised Partial Least Squares-Discriminate Analysis (PLS-DA), built on the same original dataset as PCA, showed that the stage of lactation had more impact than diet on the first two components (Supplementary Figure S1). Thus, a mean-centering and standardization by lactation stage of all variables was performed to assess the effect of diet per se. The resulting PLS-DA score plot revealed a clear-cut separation of diet groups (Figure 4b) overall. Most of the genes studied contributed to diet separation with a VIP over 0.8. *Acacb*, *Fabp3*, *Plin2*, *Glut1*, *B4galt1*, *Acads*, *Cs*, *Scd1*, *Lalba*, *Csn2*, *Slc7a5*, *Mtco1* and *Igflr* were the genes the most associated to the FEN diet while *Sod1*, *Fabp4*, *Dgat2*, *Pparg*, *Prlr*, *C3* and *Pgm1* were the most associated to the CTL diet (Figure 4a,c).
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Effect of fenugreek on the concentration of insulin (a.1,a.2), prolactin (b), IGF-1 (c), leptin (d) and estrogens (e) in maternal plasma along the lactation period. Figure 3a.2 represents insulin concentration centered and reduced by period. Values are mean ± SEM with *n* = 4–6 at L12 and L18 and *n* = 11 at Inv1. Data were analyzed with 2-way ANOVA followed by Sidak's post-hoc test for diet effect and Tukey's post-hoc test for period effect: *p* < 0.1 and *: *p* < 0.05 for diet comparison; #: *p* < 0.05, and ###: *p* < 0.001 for comparison between a given period and the immediately preceding period. CTL: control diet group and FEN: fenugreek-supplemented diet group throughout the lactation period. L12 and L18: post-natal day of lactation 12 and 18, Inv1: first day of involution.</figcaption></figure>

Fenugreek modulated mammary gland gene expression mostly at L12 and Inv1. By focusing on each lactation time point with the PLS-DA, the predictive ability of the L12- (Figure 4e) and Inv1 (Figure 4h)-PLS-DA models was more important (R2(Y) cum and Q2 cum values above 0.9 and 0.7, respectively) than that of L18-PLS-DA model (Q2 cum = 0.39). At L12, the genes contributing the most to the separation of diets were *Acacb*, *Fabp3*, *Fasn*, *Glut1*, *Pgm1*, *B4galt1*, *Slc7a5*, associated with the FEN diet, but also *Fabp4*, *Dgat2* and C3 associated to the CTL diet (VIP > 0.8, Figure 4d,f). Moreover, most of the genes overexpressed in the FEN diet at L12 were related to the synthesis of milk constituents, energy metabolism and the regulation of milk secretion (Tables 2–5). At Inv1, some genes contribute strongly to the separation of both diets such as *Mtco1*, *Cs*, *Aqp1*, *Akt1*, *Cat*, *Plin2* and *Csn2* associated with the FEN diet, but also *Lpl*, *Dgat2*, *Pgm1*, *Ugp2*, *Pparg* and *Sod1* associated with the CTL diet (Figure 4g,i). Genes associated with the FEN diet were mostly related to milk secretion and energy metabolism.

At L18, diets were poorly separated from each other by PLS-DA and no difference was observed
in mammary gene expression between CTL and FEN diets.

3.3.1. Fenugreek Supplementation Increases the Expression of Mammary Genes Involved in Lipid, Lactose and Protein Synthesis

The effect of fenugreek was studied individually on genes involved in the synthesis of milk constituents (Table 2). The FEN diet significantly increased the overall expression of genes related to milk lipid synthesis (Lpl, Fas, Acaca, Fabp3), lactose synthesis (Glut1, Ugp2, B4galt1, Lalba) and protein synthesis (Lalba and Csn2) (Table 2). Moreover, Acacb, Fabp3, Plin2, Glut1, B4galt1, Scd1, Lalba, Csn2 and Slc7a5 were strongly associated with the FEN diet in the overall PLS-DA (VIP > 1, Figure 4a,c) and highly predictive of FEN group (regression coefficient ≥ 0.5).
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Fenugreek supplementation impacts mammary gland gene expression in CTL and FEN groups, regardless of lactation stage. PLS-DA score plot (on the first two components) built on all normalized variables (that have been mean-centered and scaled by lactation time point), revealed an effect of fenugreek supplementation on mammary gland gene expression during the entire period of lactation (b), at L12 (e) and Invl (h). PLS-DA loading plot presented clusters of variables highly associated to FEN and CTL groups' overall lactation (a), at L12 (d) and Invl (g) with high regression coefficient (i.e., high prediction vector) for FEN class overall lactation (c) or high VIP scores on the two first PLS-DA components, at L12 (f) and Invl (i). CTL: control diet group and FEN: fenugreek-supplemented diet group throughout the lactation period. L12 and L18: post-natal day of lactation 12 and 18, InV1: first day of involution.</figcaption></figure>
### 3.3.2. Fenugreek Supplementation Increases Expression of Genes Involved in Fuel Metabolism, Particularly Fatty Acid β-Oxidation

The overall expression of *Cpt1a*, *Acads* and *Cs* was significantly increased with the FEN diet compared to the CTL diet, with a particular overexpression of lipid β-oxidation-related genes, *Cpt1a* and *Acads*, at L12: 3.84-fold and 2.77-fold (*p* < 0.001), respectively (Table 3). Moreover, *Acads*, *Cs* and *Mtco1* were among the most powerful discriminators associated with the FEN diet in the overall PLS-DA (VIP > 0.8, Figure 4a,c).

### 3.3.3. Fenugreek Increases the Expression of Regulatory Factors Involved in Protein and Lipid Metabolism

The FEN diet had no significant overall effect on the gene expression of factors regulating milk protein and lipid synthesis (Table 4) but induced an early transitory significant overexpression of *Lxra* and *Mtor* at L12 (2.38-fold, *p* = 0.01 and 1.84-fold, *p* = 0.032, respectively) contributing, with *Pparg*, to discrimination of the FEN diet group in the PLS-DA model at L12 (VIP > 0.8, Figure 4d,f) Additionally, *Pparg* and *Akt1* were significantly upregulated at L12 compared to L18 (1.85-fold, *p* = 0.046 and 1.62-fold, *p* = 0.035, respectively) in the FEN diet, whereas similar expression was observed for the CTL diet (*p* ≥ 0.98).

### 3.3.4. Fenugreek Upregulates the Expression of Insulin, GH and IGF-1 Receptors

The effect of fenugreek was studied on receptors of the main lactogenic hormones (Table 5). *Igf1r* expression was the most enhanced by the FEN diet throughout lactation (*p* = 0.033; Figure 4a,c), specifically at L12 (2.60-fold; *p* = 0.012 and VIP > 0.8, Figure 4d), with a similar trend observed for *Insr* and *Ghr* expressions (VIP > 0.8, Figure 4d). Interestingly, normalized gene expression of *Igf1r* was strongly correlated with those of *Ghr* and *Insr* overall lactation (*r* = 0.76 and 0.71 respectively, Supplementary Figure S3).

### 3.4. Fenugreek Stimulates Oxytocin Expression at the Pituitary Level

The FEN diet increased the overall expression of *Oxt*, *C3* and *Trhr*, mostly at L18 for *Oxt* and *C3* (*p* = 0.029 and *p* = 0.004, respectively), and tended to increase the expression of *Vipr2*. In addition, *Drd2* was overexpressed in the FEN group at *Inv1* (*p* = 0.050), whereas it was on average under-expressed in this group during lactation (interaction: *p* = 0.040).

### 3.5. Fenugreek Supplementation Increases Plasma Insulin Concentration

The FEN diet had no significant overall impact on the concentration of these hormones (Figure 4). However, plasma estrogen concentration tended to be higher in the FEN compared to CTL group (*p* = 0.082), particularly at *Inv1* (+35.5%, *p* = 0.056). Moreover, variation of insulin concentration was not homogenous between lactation and involution (*p* = 0.035 for the Brown–Forsythe test), impairing our ability to detect any effect of supplementation during lactation. By centering and reducing insulin concentration by time period, a rise in insulin concentration was observed in the FEN group at L12 (+77.1%, *p* = 0.043), whereas no difference was observed at L18 and *Inv1* (interaction: *p* = 0.040). Finally, a significant interaction was observed for prolactin concentration explained by the precipitous drop between L12 and L18 in the CTL group (94.1%, *p* < 0.001), whereas it remained steady in the FEN group (*p* = 0.995). However, the large prolactin concentration variations observed during lactation (mean CV of 102.1% and interval of (2.87; 528.2) ng/mL) makes it difficult to draw any conclusion from these values.

## 4. Discussion

The findings of the current study confirm that fenugreek increases milk production in a model of rat nursing a large litter. They further demonstrate that the impact of dietary fenugreek supplementation
on milk production is associated with a stimulation of milk macronutrient synthesis and mammary fuel metabolism. These effects were mainly observed at L12 and Inv1 but not at L18, suggesting that fenugreek promotes an earlier peak of lactation and maintains milk production longer, whereas it does not seem to have further beneficial effect on milk macronutrient synthesis once peak lactation is reached. Our findings finally suggest that the effect of fenugreek on milk production at mid-lactation may be mediated by a stimulation of insulin secretion and a modulation of the insulin/GH/IGF-1 axis, while its action on maintaining lactation until the first day of mammary involution could be due to its estrogenic effect.

## 4.1. *The Time Course of the Gene Expression in Mammary Gland and Lactogenic Hormones Suggests Peak Lactation Is Closer to L18*

A noticeable result is the importance of lactation stage on mammary gland gene expression, as illustrated by the clear-cut separation of lactation time points in both PCA and PLS-DA models (Figure 2a and Supplementary Figure S1). The first day of involution (Inv1) presented clearly different gene expression pattern from those at mid- and late-lactation.

Most genes, particularly those following **kinetics 1** and 2 (Figure 2c,d), were downregulated during involution, implying they were overexpressed at mid- and late-lactation, highlighting their central role for milk production and regulation [32]. Genes following **kinetic 1** were related to the first steps of milk lipid and lactose synthesis [19], such as glucose and fatty acid uptake (*Glut1*, *Lpl*), fatty acid de novo synthesis (*Fasn*, *Acaca*) and galactose synthesis (*Ugp2*), to mammary energetic metabolism (*Cs*, *Mtco1*) and to milk synthesis regulatory factors (*Stat5*, *Akt1*, *Pparg*). These genes reached their peak expression between mid- and late-lactation and were probably under primary regulation of *Prlr* whose expression was highly correlated with all ($r > 0.8$). Genes following **kinetic 2** were associated with the final steps of milk constituent secretion [19], such as triglyceride synthesis (*Scd1*, *Dgat1*), lipid droplet formation (*Plin2*), lactose synthesis (*B4galt1*, *Lalba*), main milk protein synthesis (*Csn2*, *Wap*) and water inflow (*Aqp1*, *Aqp3*). They seemed to be mostly controlled by *Insr* ($r \ge 0.75$ with all) but also *Esr1*, mainly for milk protein genes ($r \ge 0.85$). To the contrary, genes following **kinetic 3** (Figure 2e) were upregulated at Inv1, suggesting that their expression was repressed during lactation. These genes are related to metabolic pathways repressed during lactation such as mitochondrial β-oxidation (*Acads*, *Cpt1a*) [32,33,42] or have a secondary role on regulation (*Igf1r*, *Lxra*) [33,43,44] or on lipid synthesis (*Fabp4*, *Dgat2*) [33,45] compared to other regulatory factors (*Prlr*, *Srebf1*) or other family members (*Fabp3*, *Dgat1*). Moreover, in the present study, the drop in the expression of de novo fatty acid synthesis genes at involution is correlated (Supplementary Figure S2) with the concomitant increase in the expression of β-oxidation genes likely related to energy requirements for mammary tissue remodeling. Finally, *Oxtr* was overexpressed in involuting mammary gland probably because myoepithelial cells were preserved longer than lactocytes and contracted to facilitate removal of degenerating epithelial cells [46].

Interestingly, genes related to milk constituent secretion followed the **kinetic 2** pattern, and as such, were overexpressed at L18 compared to L12. These genes are probably the more strongly correlated with actual milk secretion, suggesting that, among the three periods studied, milk production was likely at its maximum at L18. This is consistent with (i) the time course of plasma insulin, leptin and IGF-1, which all declined to their minimal concentration at L18. Indeed, the plasma concentration of these hormones is known to decline during lactation partly due to increased hormone extraction by the mammary gland [47–49] and they are negatively correlated with milk production [50]. (ii) The highest milk consumption by rat pups reported to occur at L15-L16 as they only start to consume solid food significantly at L17 and drinking water at L19 [51].

## 4.2. *Fenugreek Supplementation Stimulates Milk Macronutrient Synthesis*

Dietary fenugreek supplementation increased the overall expression of genes involved in fatty acid uptake (*Lpl*), de novo synthesis (*Fasn*, *Acaca*, *Acacb*), transport (*Fabp3*) and, to a lesser extent,
lipid secretion (*Plin2*). Fenugreek also stimulated the expression of genes involved in glucose uptake (*Glut1*), galactose (*Pgm1*, *Ugp2*) and lactose (*B4galt1*, *Lalba*) synthesis. Finally, fenugreek increased the expression of amino acid transporters (*Slc7a5*) and main milk proteins (*Csn2*, *Lalba*). These results suggest an increase in the synthesis of the 3 macronutrients of milk namely lipid, lactose and protein in the mammary gland [19] and are consistent with the stimulating effect of fenugreek on milk flow by 16% reported in our earlier study [15] and 15.3% in this study. More specifically, genes involved at every step of protein and lactose synthesis were significantly overexpressed, whereas only genes involved in the fatty acid synthesis were dramatically overexpressed for lipid synthesis. The expression of genes involved in triglyceride synthesis and lipid secretion was not or only slightly affected. This relative lack of conversion of fatty acids to triglycerides may account for the decrease in milk lipid concentration found in some fenugreek supplementation studies [52–54]. Otherwise, the stimulation by fenugreek of genes involved in mitochondrial energy metabolism (*Cs*, *Atp5f1a*, *Mtco1*) and particularly those involved in β-oxidation (*Cpt1a*, *Acads*) suggests that increased ATP synthesis took place in the mammary gland to meet the increased energy requirements associated with increased milk production [22]. Nevertheless, contrary to Liu et al. [55], we were unable to observe a significant increase in water inflow in response to galactagogue even though mean *Aqp1* and *Aqp3* expressions were increased at L12 and L21.

### 4.3. Fenugreek Supplementation Advances Peak Lactation and Maintains Milk Production

The lack of separation according to diet on the PLS-DA plot, and the lack of any difference in the expression of genes related to milk constituent synthesis and fuel metabolism in FEN and CTL groups at L18, suggests that the milk synthesis peak observed at late-lactation mainly in control dams was not impacted by fenugreek supplementation. However, the time course of some genes from L12 to L18 differed between diets: either the gene expression increased in the CTL diet and did not change in the FEN diet (*Lpl*, *B4galt1*, *Slc7a5*, *Csn2* and *Aqp1*), or the expression did not change in the CTL diet and decreased in the FEN diet (*Acaca*, *Fabp3*, *Ugp2*, *Cpt1a*, *Acads*). In both cases, gene expression was higher earlier in the FEN group. This is consistent with the clear separation between both diets observed in the PLS-DA model built on gene expression at L12 (Figure 4d) and the fact that greater overexpression was observed at L12 for most of the genes in the FEN group. Finally, the PLS-DA model built on gene expression at Inv1 also properly separated both diets and the expression of the genes *Lalba*, *Csn2*, *Plin2*, *Slc7a5*, *Aqp1* and *Esr1* was maintained to a slightly greater level in the FEN diet. Together these results suggest that fenugreek supplementation advanced the peak of milk production, observed near L18 for control, around L12 and maintain it to a level similar to control at L18. Thus, fenugreek supplementation may enable the dam to maintain a milk production slightly higher than control during involution.

### 4.4. Fenugreek Stimulates Milk Production Mainly by Modulating Insulin and Oxytocin Secretion

The main effect of fenugreek on the expression of genes related to milk synthesis was observed at L12. During the same period, fenugreek supplementation increased dams' plasma insulin by 77.1% and the expression of *Insr* by 1.79-fold, suggesting an effect of fenugreek on insulin secretion and its action in mammary gland. Receptors to GH and IGF-1 (GHR and IGF1R) were also remarkably overexpressed in the mammary gland at L12 (Table 5). These 3 hormones are closely linked by the insulin/GH/IGF-1 axis [56]. Briefly, GH is secreted by pituitary gland and stimulates IGF-1 secretion in the liver [56,57], which shares strong homology with insulin. Both hormones can bind to IGF1R and insulin receptors (INSR) and to hybrid formed by both receptors with more or less affinity. These receptors act through similar pathway mainly by activating the AKT/mTOR pathway [56]. In addition, insulin stimulates expression of GH receptor in peripheral tissues [58] but also modulates the expression of IGF1R and INSR [59,60].

We therefore hypothesize that fenugreek acts mainly by increasing insulin secretion by maternal pancreas β-cells. This is consistent with the hypoglycemic effect of fenugreek and its capacity to
stimulate insulin secretion through the potential action of trigonelline and 4-hydroxyisoleucin [23,24]. Together with a potential direct action of fenugreek compounds [24], the rise in plasma insulin under lactation-specific conditions, i.e., hypo-insulinemia and insulin-resistance in peripheral tissues except mammary gland, could explain the *Igf1r*, *Insr* and *Ghr* overexpression in the mammary gland [58,60,61]. The increase of IGF-1 and insulin binding to their receptor and to their hybrid receptor leads to an activation of the AKT/mTOR pathway, as suggested by overexpression of *Mtor* at L12 (by 1.84-fold) and of *Akt1* to a lesser extent (Table 4) and results in activation of lactose and protein synthesis [19]. GHR could also act secondarily on protein and lactose synthesis through the JAK2/STAT5 pathway [62]. Concerning lipid metabolism, *Lxra* and *Pparg* are the most overexpressed regulatory factors with fenugreek supplementation (Table 4). Both are stimulated by insulin [63–65] and likely IGF-1, and *Lxra* expression is also stimulated by GH [64]. *Lxra* has been shown to stimulate de novo fatty acid in the mammary gland [66] independently of SREBF1 but not fatty acid desaturation and triglyceride synthesis [44], which is consistent with fenugreek effect found only for fatty acid uptake and de novo synthesis. *Lxra* was also related to activation of lipolysis in adipocyte [63,64] and may be responsible for the maintenance of β-oxidation-related genes at L12, as suggested by the strong correlation with *Acads* and *Cpt1a* (*r* > 0.7, Figure 2e). *Pparg* had a secondary milk lipid regulatory role in rodent compared to ruminant [19] but promote de novo fatty acid synthesis, triglyceride synthesis and lipid secretion [19,45,65]. Finally, *Igf1* is known to be an important activator of fatty acid β-oxidation [56] that could counteract insulin inhibitory effect [20,42] on this pathway and notably at involution. It may act independently on the maintenance of *Acads* and *Cpt1a* expression, as suggested by the very strong correlation between *Igf1r*, *Acads* and *Cpt1a* expressions (*r* ≥ 0.84, Figure 2e). Suggested mechanisms associated to the galactologue effect of fenugreek are summarized in Figure 5.

Another mechanism by which fenugreek could increase milk flow may be its capacity to favor milk ejection by stimulating oxytocin secretion [22]. Indeed, fenugreek is a well-known oxytocic agent used to promote uterine contractions at delivery [23,25] and oxytocin expression was overexpressed in the pituitary gland during lactation in the fenugreek-supplemented diet, particularly at L18 (1.88-fold). Even if *Oxtr* expression was not significantly increased in mammary gland, fenugreek might maintain greater milk flow at L18 by promoting milk ejection (Figure 5).

Nevertheless, we could not observe increase in mother prolactin, as suggested in previous studies [26,53], because of strong variability in plasmatic prolactin concentration. These variations were probably due to the too-short dam/pup separation time before sacrifice (60 min) [67], which made prolactin concentration more likely representative to the time since last suckling.

### 4.5. Fenugreek Stimulates Estrogenic Activity at the End of Lactation and Increased Energy Availability for Milk Production

At L18 and Inv1, fenugreek-supplemented dams tended to have higher plasma estrogens, and in the pituitary, had higher expression of C3, which is known to be a strong estrogenic biomarker in the uterus [68]. Moreover, *Esr1* tended to be upregulated in the mammary gland at L21 (Figure 3g) and was strongly correlated with genes overexpressed at this period (*r* ≥ 0.74 with *Csn2*, *Lalba*, *Aqp1*, *Plin2* and *Akt1*, Figure 2g and Supplementary Figure S1). As fenugreek is known to contain several estrogenic compounds [23,27], those might contribute to maintaining the expression of milk secretory genes at Inv1 partly through the mediation of *Akt1* [69].

Finally, fenugreek supplementation increased global food consumption of dams by 10%, in line with earlier studies [52,70,71], principally at mid-lactation (Supplementary Figure S4) that could partly account for fenugreek lactogenic effect. Indeed, lactation is a period of high energy requirement characterized by an important hyperphagia, in particular in small animals where food intake can be increased by 450% [72]. Rats have little lipid storage during lactation and energy requirement for milk production almost exclusively comes from increased food consumption [22,73]. Thus, by enhancing food consumption, fenugreek might enhance energy availability for milk production. This increase in food intake observed after fenugreek supplementation might be mediated (i) by a family of its bioactive
compounds, steroid saponins, that stimulate eating motivation [74,75] and could act directly on the hypothalamus or through stimulation of prolactin, or (ii) by a greater peak of prolactin during suckling, as suggested by overexpression of receptors of prolactin secretion activators in the pituitary gland (thyrotropin releasing hormone and vasoactive intestinal peptide) [76]. Indeed, lactation hyperphagia is under the primary control of prolactin which decreases leptin secretion in adipocytes and, in the hypothalamus, stimulates orexigenic neuropeptide Y secretion and decreases its sensitivity to leptin [22]. Moreover, milk production is also known as an important determinant of food intake in rat [72,77]. Indeed, greater milk yield leads to higher dam's energy expenditure and likely, enhanced suckling stimulus, as mammary gland was reported to be emptied at each feeding for a rodent litter-size model of 12 pups [78], resulting in greater food intake. This is in accordance with the peak dietary intake observed when fenugreek dams started their peak milk production.

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Suggested mechanisms of action of fenugreek on milk synthesis and secretion. Green arrows denote positive activation, bold and thin arrows represent main and secondary fenugreek action pathways, respectively. Dotted lines represent hypothetical pathway of action. Dark and light green boxes represent significant or trend increase respectively of gene expression or plasma concentration overall for upper box and at each lactation period for lower boxes. White boxes depict the lack of modulation by fenugreek.</figcaption></figure>
Otherwise, we observed higher leptin and insulin level in plasma of fenugreek-supplemented dams likely caused by fenugreek stimulation of insulin secretion [79]. Yet, these hormones are known to exhibit important anorexigenic effect under normal physiological conditions [80]. However, during lactation, because of central resistance to their signals [22], insulin and leptin had only secondary roles in regulating lactation-associated hyperphagia compared to suckling stimulus and prolactin action [77,81]. Moreover, under normal conditions, orexigenic effect of fenugreek was observed despite an increase of plasma insulin [74]. Thus, effect of fenugreek on food intake was probably not related to the increase of both insulin and leptin plasma concentrations.

Whatever the cause, overall food increase in fenugreek-supplemented dams was of 10% and only of 5.3% (p = 0.041) related to dams' body weight, so it cannot alone account for the 15% increase in milk production observed in our study.

**5. Conclusions**

The current study confirms the ability of dietary fenugreek supplementation to promote milk production, and sheds light into its action at the molecular level. Indeed, milk flow was increased by 15.3% and most of the genes related to milk lipid, lactose and protein synthesis as well as energy metabolism were overexpressed. Moreover, our findings suggest that fenugreek may act by extending the duration of peak lactation until the beginning of involution rather than intensifying it. According to our data, the main mechanism of action of fenugreek on milk production is likely through activation of the insulin/GH/IGF-1 axis, which leads to a greater mammary sensitivity to these hormones. Consequently, milk protein and lactose synthesis could be mainly activated by the AKT/mTOR pathway, while milk lipid synthesis and energetic metabolism could be stimulated by LXRα, PPARγ and probably IGF1R itself. Moreover, fenugreek would stimulate milk ejection by an increase of pituitary oxytocin secretion. However, data from gene expression must be confirmed by proteomic and histologic analysis.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/2073-4425/11/10/1208/s1, Figure S1: PLS-DA score plot on non-normalized mammary gene expression separating groups according to the diet and the lactation period. Figure S2: Pearson correlation table between all mammary gland genes' expressions on row data. Figure S3: Pearson correlation table between all mammary gland genes' expressions on data normalized by the period of lactation. Figure S4: Food intake of dams supplemented or not with 1 g/kg/day fenugreek during lactation. Table S1: Validation of RNA extraction and primer design.

**Author Contributions:** Conceptualization, T.S., M.-C.A.-G., D.D. and C.-Y.B.; data curation, T.S. and M.-C.A.-G.; formal analysis, T.S.; funding acquisition, M.-C.A.-G., D.D. and C.-Y.B.; investigation, T.S., A.G. and I.G.; methodology, P.d.C.; project administration, C.-Y.B.; supervision, M.-C.A.-G. and C.-Y.B.; writing—original draft, T.S.; writing—review and editing, T.S., M.-C.A.-G., D.D., P.d.C. and C.-Y.B. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by Laboratoire FRANCE Bébé Nutrition (Laval, France), grant number Cifre 2017.0004.

**Acknowledgments:** The authors want to thank specifically Blandine Castellano (PhAN, Nantes, France), Sandrine Suzanne and Amandine Lefebvre (UTE, Nantes, France) for their precious help during the animal experiments. We also thank Anthony Pagniez (PhAN, Nantes, France) for his appropriate help with RT-qPCR experiments. Finally, we thank Patrick Nguyen (ONIRIS, Nantes, France) for giving us access to the infrared spectrometer that allowed us to carry out milk production measurements.

**Conflicts of Interest:** Thomas Sevrin received financial support through grants from FRANCE Bébé Nutrition. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and they agreed with the decision to publish the results. The authors declare no conflict of interest.

**Abbreviations**

GH: growth hormone; IGF-1: insulin-like growth factor 1; GHR: GH receptor; IGF1R: IGF-1 receptor; INSR: insulin receptor; PCA: principal component analysis; PLS-DA: partial least square discriminant analysis; AKT: serine/threonine kinase; mTOR: mechanistic target of rapamycin.
Appendix A List of the Primers' Sequences

<table><caption>Table A1. List of the primers' sequences used to measure each gene by qPCR with the amplicon size and the efficiency.</caption><thead><tr><th>Gene Name</th><th>Gene Symbol</th><th>Function</th><th>Accession Number</th><th>Primer Sequences (5' to 3')</th><th>Amplicon (pb)</th><th>Efficiency (%)</th></tr></thead><tbody><tr><td>Lipoprotein lipase</td><td>Lpl</td><td>Triglyceride hydrolysis and mediation of lipoprotein uptakes</td><td>NM_012598</td><td>F: AACTGCCACTTCAACCACAGC<br>R: CATACATTCCCTGTCACCGTCCA</td><td>70</td><td>95.8</td></tr><tr><td>Fatty acid synthase</td><td>Fasn</td><td>Fatty acid biosynthesis, mainly palmitate and long chain fatty acid synthesis</td><td>NM_017332</td><td>F: CGCCGTGGTGCTGGAGATTG<br>R: CTTGCCGAGGTTGGTGAGGAAG</td><td>142</td><td>98.3</td></tr><tr><td>Acetyl-CoA carboxylase α</td><td>Acaca,<br>Acc1</td><td>Catalyze the rate limiting step in fatty acid biosynthesis: malonyl-CoA synthesis</td><td>NM_022193</td><td>F: CGATTCCCATCCGCCTCTTCC<br>R: GGTCCTGCTTGTCTCCATACG</td><td>127</td><td>92.3</td></tr><tr><td>Acetyl-CoA carboxylase β</td><td>Acacb</td><td>Fatty acid biosynthesis, involved in inhibition of fatty acid oxidation</td><td>NM_053922</td><td>F: GACGCCAGCAACATCACTTCG<br>R: GCCTCTCGGTCCACCTTACG</td><td>107</td><td>85.5</td></tr><tr><td>Fatty acid binding protein 3</td><td>Fabp3</td><td>Uptake, intracellular metabolism and transport of fatty acids</td><td>NM_024162</td><td>F: ACAGGAAGGTCAAGTCGGTC<br>R: ATGGGTGAGAGTCAGGATGAG</td><td>127</td><td>93.3</td></tr><tr><td>Fatty acid binding protein 4</td><td>Fabp4</td><td>Lipid binding protein involved in fatty acid trafficking, main adipocyte FABP</td><td>NM_053365</td><td>F: TGTGGGGACCTGGAAACTCGT<br>R: CGAAGCCAACTCCCACTTCTTT</td><td>71</td><td>92.9</td></tr><tr><td>Stearoyl-CoA desaturase</td><td>Scd,<br>Scd1</td><td>Fatty acid desaturation, biosynthesis of oleic acid.</td><td>NM_139192</td><td>F: GTTGGGTGCCTTATCGCTTTCC<br>R: CTCCAGCCAGCCTCTTGTCTA</td><td>114</td><td>97.8</td></tr><tr><td>Diacylglycerol<br>O-acyltransferase 1</td><td>Dgat1</td><td>Catalyzes the synthesis of triglyceride from diglyceride and acyl-CoA</td><td>NM_053437</td><td>F: GGCATCATACTCCATCATCTTC<br>R: CCCACTGACCTTCTTCCC</td><td>115</td><td>97.1</td></tr><tr><td>Diacylglycerol<br>O-acyltransferase 2</td><td>Dgat2</td><td>Catalyzes the synthesis of triglyceride from diglyceride and acyl-CoA</td><td>NM_001012345</td><td>F: GGTCATCTCAGTCCTACAG<br>R: CTATCAGCCAGCAGTCAG</td><td>101</td><td>88.6</td></tr><tr><td>Perilipin 2</td><td>Plin2</td><td>Coating of intracellular lipid droplets and milk fat globules, involved in milk lipid secretion.</td><td>NM_001007144</td><td>F: CTTCTTCATTGACCTGCGAC<br>R: CACCCTTACCTTAAGTTTCTCCT</td><td>108</td><td>97.1</td></tr><tr><td>Solute carrier family 2,<br>facilitated glucose<br>transporter member 1</td><td>Slc2a1, Glut1</td><td>One of major glucose transporters, located in the cell membrane</td><td>NM_138827</td><td>F: CCAATATGTGGAGCAACTGTG<br>R: GGAAGCGATCTCATCGAAGG</td><td>129</td><td>108.2</td></tr><tr><td>Phosphoglucomutase 1</td><td>Pgm1</td><td>Interconvert glucose-1-phosphate and glucose-6-phosphate</td><td>NM_017033</td><td>F: CCCTTCACAGTGGAGATCGT<br>R: GTCGATGCGGATCTTCAGTC</td><td>122</td><td>96.1</td></tr><tr><td>UDP-glucose<br>pyrophosphorylase 2<br>β-1,4-</td><td>Ugp2</td><td>Convert glucose-1-phosphate to UDP-glucose</td><td>NM_001024743</td><td>F: CGGAAGATTCGATTCAACCC<br>R: GTTCCCAAACCACCATTGAG</td><td>111</td><td>93.0</td></tr><tr><td>galactosyltransferase 1</td><td>B4galt1</td><td>Catalytic subunit of lactose synthase complex, convert UDP-galactose and glucose into lactose</td><td>NM_053287</td><td>F: CGCTTTGTGTTCAGTGATGTGG<br>R: AACGTAAGGCAGGCTAAACC</td><td>122</td><td>100.7</td></tr><tr><td>Lacatalbumin α</td><td>Lalba</td><td>Regulatory subunit of lactose synthase complex, one of the main milk proteins</td><td>NM_012594</td><td>F: CTTGCTTGAATGGACCTGTG<br>R: GCTACTCTTGCACCAATTTCTG</td><td>126</td><td>92.7</td></tr><tr><td>Casein β</td><td>Csn2</td><td>Principal milk protein and primary source of essential amino acids for sucking young</td><td>NM_017120</td><td>F: AAACATCCAGCCTATTGCTC<br>R: CATCTGTTTGTGCTTGGGAA</td><td>118</td><td>96.6</td></tr><tr><td>Whey acidic protein</td><td>Wap</td><td>Whey protein associated with mammary gland, a major constituent of milk</td><td>NM_053751</td><td>F: CATGTCTTCAACTCAGTTCAGTCC<br>R: CATGTCATTCTGGGCACACTC</td><td>104</td><td>95.8</td></tr><tr><td>Solute carrier family 7<br>member 5</td><td>Slc7a5</td><td>Amino acid transporter, mediate the uptake of large neutral amino acids (branch-chain, aromatic)</td><td>NM_017353</td><td>F: GGGAAAGGACATAGGACAAGG<br>R: TCAGATAGTTCCATCCTCCGT</td><td>136</td><td>94.5</td></tr><tr><td>Aquaporin 1</td><td>Aqp1</td><td>Integral membrane protein, allow the passive transport of water along osmotic gradient.</td><td>NM_012778</td><td>F: GGCATTGAGATCATTGGCAC<br>R: GTGTAGTCAATGGCCAGCAG</td><td>142</td><td>97.8</td></tr><tr><td>Aquaporin 3</td><td>Aqp3</td><td>Passive water transporter located in the basolateral cell membrane</td><td>NM_031703</td><td>F: GGGCTCTACTATGATGCAATCTG<br>R: GCTGCTGTGCCTATGAACTG</td><td>148</td><td>95.5</td></tr></tbody></table>
<table><caption>Table A1. *Cont.*</caption><thead><tr><th>Gene Name</th><th>Gene Symbol</th><th>Function</th><th>Accession Number</th><th>Primer Sequences (5' to 3')</th><th>Amplicon (pb)</th><th>Efficiency (%)</th></tr></thead><tbody><tr><td>Carnitine palmitoyltransferase 1A</td><td>Cpt1a</td><td>Initiate β-oxidation by translocating long chain fatty acids across mitochondrial inner membrane</td><td>NM_031559</td><td>F: TGCCTGCCAGTTCCATTAGC<br>R: GTCTCACTCCTCTTGCCAACAG</td><td>143</td><td>90.1</td></tr><tr><td>Acyl-CoA dehydrogenase short chain</td><td>Acads</td><td>Catalyze the first step of mitochondrial fatty acid β-oxidation</td><td>NM_022512</td><td>F: CGGCAGAACAGGGTATCAG<br>R: AAGATGAGGTTAGCTGTAGATGAG</td><td>115</td><td>99.1</td></tr><tr><td>Pyruvate dehydrogenase E1 subunit α 1</td><td>Pdha1</td><td>Subunit of the pyruvate dehydrogenase complex that convert pyruvate to acetyl-CoA and CO2</td><td>NM_001004072</td><td>F: CACGGACCATCTCATCACTG<br>R: CTCCTCTTCGTCTGTTAATTCTG</td><td>100</td><td>91.2</td></tr><tr><td>Citrate synthase</td><td>Cs</td><td>Krebs cycle enzyme that catalyzes the synthesis of citrate from oxaloacetate and acetyl-CoA</td><td>NM_130755</td><td>F: GCTATAGTATCCCTGAGTGCCA<br>R: AGAGCCAAGAGACCTGTTCC</td><td>128</td><td>91.1</td></tr><tr><td>Cytochrome c oxidase I, mitochondrial</td><td>Mtco1, mtCox1</td><td>Subunit of the cytochrome c oxidase, catalyzes reduction of oxygen to water for oxidative phosphorylation</td><td>YP_665631</td><td>F: GCCTAGATGTAGACACCCGAG<br>R: TATTTCCTCCATGTAGTGTAGCGA</td><td>107</td><td>97.0</td></tr><tr><td>ATP synthase F1 subunit α</td><td>Atp5f1a</td><td>Subunit of the mitochondrial ATP synthase, catalyzes the ATP synthesis during oxidative phosphorylation</td><td>NM_023093</td><td>F: AGCGTTTCAATGATGGGACTG<br>R: GTACTTCATGGCATCTGCGTC</td><td>124</td><td>94.3</td></tr><tr><td>Sterol regulatory element binding transcription factor 1</td><td>Srebf1</td><td>Nuclear transcription factor that bind to the sterol regulatory element 1, regulates the transcription of genes important for sterol biosynthesis and lipogenesis</td><td>NM_001276707</td><td>F: GCAGCTGATGGAGACAGGGA<br>R: CGACAGCGTCAGAACAGCTA</td><td>76</td><td>99.3</td></tr><tr><td>Peroxisome proliferator activated receptor γ (PPARγ)</td><td>Pparg</td><td>Transcription factor that regulate expression of genes involved in lipid metabolism</td><td>NM_001145366</td><td>F: CCGTTCACAAGAGCTGACCC<br>R: TCGCACTTTGGTATTCTTGGAGC</td><td>73</td><td>94.7</td></tr><tr><td>Thyroid hormone responsive</td><td>Thrsp, Spot14</td><td>Play a role in regulation of lipogenesis and specially triglycerides with medium-chain fatty acids.</td><td>NM_012703</td><td>F: CACATCCTTACCCACCTGAC<br>R: CGTGTAAAGTGATCTTCCATAGAG</td><td>109</td><td>88.5</td></tr><tr><td>Nuclear receptor subfamily 1, group H, member 3 (LXRα)</td><td>Nr1h3, Lxra</td><td>Nuclear oxysterol receptor that regulate inflammation and homeostasis of lipids and cholesterol</td><td>NM_031627</td><td>F: CAGAGCCTACAGAACTTCGT<br>R: CAGCTCAGCACATTGTAATGG</td><td>123</td><td>101.0</td></tr><tr><td>AKT serine/threonine kinase 1 (AKT)</td><td>Akt1</td><td>Key signaling pathway protein notably linked to mTOR signaling, regulate metabolism, proliferation</td><td>NM_033230</td><td>F: CTACTATGCCATGAAGATCCTC<br>R: CTACTATGCCATGAAGATCCTC</td><td>124</td><td>93.9</td></tr><tr><td>Mechanistic target of rapamycin kinase (mTOR)</td><td>Mtor</td><td>Central regulator of cell metabolism, growth and survival, mediate nutrient and energy signal, key regulator of mRNA translation</td><td>NM_019906</td><td>F: ACCAATTATACTCGCTCCCTG<br>R: CATAGCAACCTCAAAGCAGTC</td><td>147</td><td>90.9</td></tr><tr><td>Signal transducer and activator of transcription 5A (STAT5)</td><td>Stat5a, Stat5</td><td>Involved in signal transduction and transcription activation in response to cytokine, growth factors and prolactin</td><td>NM_017064</td><td>F: AGAACACCCGCAATGATTACAG<br>R: GTGACATGTTTCTGAAGTGGG</td><td>129</td><td>97.8</td></tr><tr><td>Prolactin receptor</td><td>Prlr</td><td>Prolactin binding and signaling mainly through JAK2/STAT5 and Akt/mTOR pathways</td><td>NM_001034111</td><td>F: GAAATGCCAAATGACTTCACCT<br>R: TATAGCCCTTCAAAGCCACTG</td><td>111</td><td>91.5</td></tr><tr><td>Insulin receptor</td><td>Insr</td><td>Insulin binding and signaling mainly through AKT/mTOR and MAPK pathways</td><td>NM_017071</td><td>F: TTCATTCAGGAAGACCTTCGA<br>R: CAGGCCAGAGATGACAAGTGAC</td><td>259</td><td>95.4</td></tr><tr><td>Oxytocin receptor</td><td>Oxtr</td><td>G-protein coupled receptor binding oxytocin, activate PIP3-Ca2+ pathway</td><td>NM_012871</td><td>F: TTCATTCAGGAAGACCTTCGA<br>R: GAGTTCGTGGAAGAGGTGAC</td><td>146</td><td>88.3</td></tr><tr><td>Growth hormone receptor</td><td>Ghr</td><td>Growth hormone binding and signaling mainly through JAK2/STAT5 pathway</td><td>NM_017094</td><td>F: TTCTTCGTGCAGATGTGGAG<br>R: ATGGAAACCGGAAATCTTCTTCAC</td><td>115</td><td>102.1</td></tr><tr><td>Insulin-like growth factor 1 receptor</td><td>Igfr1</td><td>IGF-1 binding and signaling mainly through AKT/mTOR and MAPK pathways.</td><td>NM_052807</td><td>F: CCATAGAAAGAGGAATAACAGCAG<br>R: TACCTCCCATTCCATCAGGCA</td><td>108</td><td>88.9</td></tr></tbody></table>
<table><caption>Table A1. *Cont.*</caption><thead><tr><th>Gene Name</th><th>Gene Symbol</th><th>Function</th><th>Accession Number</th><th>Primer Sequences (5' to 3')</th><th>Amplicon (pb)</th><th>Efficiency (%)</th></tr></thead><tbody><tr><td>Estrogen receptor 1</td><td>Esr1, Esra</td><td>Nuclear receptor binding estrogens, act as a regulator of DNA transcription</td><td>NM_012689</td><td>F: CCACGATCAAGTTCACCTTCTG<br>R: CACATTTACCTTGATTCCTGTCCA</td><td>136</td><td>86.8 (MG)<br>92.1 (PG)</td></tr><tr><td>Complement C3</td><td>C3</td><td>Involved in innate immune response; central role in activating the complement system</td><td>NM_000064</td><td>F: CTTCTTCATTGACCTGCGAC<br>R: CACCCTTACCTTAAGTTTCTCCT</td><td>122</td><td>94.6 (MG)<br>98.4 (PG)</td></tr><tr><td>Superoxide dismutase 1</td><td>Sod1</td><td>Convert superoxide to hydrogen peroxide and oxygen, involved in response to oxidative stress</td><td>NM_017050</td><td>F: TACACAAGGCTGTACCACTG<br>R: CACACGATCTTCAATGGACAC</td><td>150</td><td>101.4</td></tr><tr><td>Catalase</td><td>Cat</td><td>Key antioxidant enzyme that convert hydrogen peroxide to water and oxygen, limit oxidative stress</td><td>NM_012520</td><td>F: CATAGCCAGAAGAGAAACCC<br>R:GAACAAGAAAGTAACCTGATGGAG</td><td>104</td><td>93.7</td></tr><tr><td>Glutathione peroxidase 1</td><td>Gpx1</td><td>Catalyze the reduction of hydroperoxides and H₂O₂ by glutathione to protect cells against oxidative damages</td><td>NM_030826</td><td>F: GTTCGGACATCAGGAGAATGG<br>R: GAAGGTAAAGAGCGGGTGAG</td><td>144</td><td>95.7</td></tr><tr><td>Prolactin</td><td>Prl</td><td>Pituitary hormone that primarily promote lactation. Play a role as growth regulator and in cell survival</td><td>NM_012629</td><td>F: AAACAGTATGTCCAAGATCGTGAG<br>R: ACTTCCGGAGGGACTTTCTG</td><td>112</td><td>99.8</td></tr><tr><td>Gh</td><td>Gh1</td><td>Pituitary hormone playing major role in growth control. Stimulate IGF1-1 secretion by the liver</td><td>NM_001034848</td><td>F: ACCTACAAAGAGTTCGAGCGT<br>R: GAAGCAATTCCATGTCAGTTCTC</td><td>147</td><td>94.2</td></tr><tr><td>Oxytocin/neurophysin I prepropeptide</td><td>Oxt</td><td>Hormone involved in contraction of smooth muscle during parturition and lactation</td><td>NM_012996</td><td>F: GGATATGCGCAAGTGTCTTCC<br>R: GAGGGCAGGTAGTTCTCCTC</td><td>140</td><td>106.8</td></tr><tr><td>Dopamine receptor D2</td><td>Drd2</td><td>G-protein coupled receptor of dopamine which inhibit adenylyl-cyclase and notably prolactin secretion</td><td>NM_012547</td><td>F: CCACTCAAGGGCAACTGTACC<br>R: ATGGGGCTATACCGGGTCC</td><td>179</td><td>91.5</td></tr><tr><td>Thyrotropin releasing hormone receptor</td><td>Trhr</td><td>G-protein coupled receptor of TRH that activate PIP3-Ca²⁺ pathway, promote release of prolactin</td><td>NM_013047</td><td>F: CAGACCGCTTTAGCACAGAG<br>R: ACTGGGTCCATTCTTCTCGG</td><td>112</td><td>90.1</td></tr><tr><td>Vasoactive intestinal peptide receptor 2</td><td>Vipr2</td><td>G-protein coupled receptor of VIP which activate adenylyl-cyclase and notably prolactin secretion</td><td>NM_017238</td><td>F: GCAGCCAAATGGAGAATCAC<br>R: TTTGCTTATGTTTCCTGGTCTG</td><td>145</td><td>104.0</td></tr><tr><td>Growth hormone releasing hormone receptor</td><td>Ghrhr</td><td>G-protein coupled receptor of GHRH that activate adenylyl-cyclase and stimulate GH transcription and secretion</td><td>NM_012850</td><td>F: GCTGTTTGCTACTTTCATCCT<br>R: CCTTGCACAGAATAGTGGAC</td><td>109</td><td>97.9</td></tr><tr><td>Adrogen receptor</td><td>Ar</td><td>Steroid hormone activated transcription factor, involved in cell proliferation and differentiation</td><td>NM_012502</td><td>F: CTCTGCCTCTGAAGTGTGGT<br>R: TACTGTCCAAACGCATGTCC</td><td>136</td><td>91.5</td></tr><tr><td>Actin β</td><td>Actb</td><td>Highly conserved proteins that are involved in cell motility, structure, integrity and signaling</td><td>NM_031144</td><td>F: CTATCGGCAATGAGCGGTTCC<br>R: GCACTGTGTTGGCATAGAGGTC</td><td>150</td><td>93.2 (PG)</td></tr><tr><td>Ubiquitously expressed, prefoldin-like chaperone</td><td>Uxt</td><td>Cofactor that modulates AR-dependent transcription, and also plays a critical role in TNF-induced apoptosis</td><td>NM_001006982</td><td>F: ATTGACCGAAAGAGTTCCCT<br>R: GTAGTTCTCTAAGTCCCTCTAGCA</td><td>108</td><td>91.3 (MG)<br>87.2 (PG)</td></tr><tr><td>Ribosomal protein S9</td><td>Rps9</td><td>One of 40S ribosomal subunit protein</td><td>NM_031108</td><td>F: GCTAAAGTTGATTGGAGAGTATGG<br>R: AGAGCGTTGCCTTCAACAG</td><td>146</td><td>99.0 (MG)<br>89.5 (PG)</td></tr><tr><td>Valosin-containing protein</td><td>Vcp</td><td>ATP binding protein involved in vesicle transport and fusion and endoplasmic reticulum function</td><td>NM_053864</td><td>F: TCAAGCGAGAGGATGAGGAG<br>R: TCACACCAATTGCCTTAAAGAG</td><td>140</td><td>92.3</td></tr></tbody></table>

Official name of each gene with official symbol given by NCBI gene (first symbol) as well as the most common symbol used if different. Functions were given by NCBI gene or GeneCards. Primers of *Actb*, *Lpl*, *Fasn*, *Acaca*, *Acacb*, *Scd1*, *Dgat1*, *Dgat2*, *Fabp4*, *Cpt1a*, *Srebf1*, *Pparg*, *Akt1*, *mTOR*, *Insr* and *Drd2* were already designed and validated in the lab, whereas other primers were designed using PerlPrimer. All primers were tested for efficiency and specificity (RT-) in mammary gland or/and pituitary gland. MG and PG: efficiency measured in mammary gland and pituitary gland respectively, if genes' expression was measured in both tissues.
## References

1. World Health Organization. *Indicators for Assessing Infant and Young Child Feeding Practices*; World Health Organization: Geneva, Switzerland, 2007.

2. Breastfeeding, S.O. Breastfeeding and the Use of Human Milk. *Pediatrics* **2012**, *129*, e827-e841. [CrossRef] [PubMed]

3. Mosca, F.; Giannì, M.L. Human milk: Composition and health benefits. *Pediatr. Med. Chir.* **2017**, *39*, 155. [CrossRef] [PubMed]

4. Victora, C.G.; Bahl, R.; Barros, A.J.D.; A França, G.V.; Horton, S.; Krasevec, J.; Murch, S.; Sankar, M.J.; Walker, N.; Rollins, N.C. Breastfeeding in the 21st century: Epidemiology, mechanisms, and lifelong effect. *Lancet* **2016**, *387*, 475-490. [CrossRef]

5. Thulier, D.; Mercer, J. Variables Associated With Breastfeeding Duration. *J. Obstet. Gynecol. Neonatal Nurs.* **2009**, *38*, 259-268. [CrossRef] [PubMed]

6. Gatti, L. Maternal Perceptions of Insufficient Milk Supply in Breastfeeding. *J. Nurs. Scholarsh.* **2008**, *40*, 355-363. [CrossRef] [PubMed]

7. Giannì, M.L.; Bettinelli, M.E.; Manfra, P.; Sorrentino, G.; Bezze, E.; Plevani, L.; Cavallaro, G.; Raffaeli, G.; Crippa, B.L.; Colombo, L.; et al. Breastfeeding Difficulties and Risk for Early Breastfeeding Cessation. *Nutrients* **2019**, *11*, 2266. [CrossRef]

8. Hillervik-Lindquist, C. Studies on Perceived Breast Milk Insufficiency A Prospective Study in a Group of Swedish Women. *Acta Paediatr.* **1991**, *80*, 1-29. [CrossRef]

9. Giugliani, E.R. Common problems during lactation and their management. *Jornal de Pediatria* **2004**, *80*, S147-S154. [CrossRef]

10. Bazzano, A.N.; Hofer, R.; Thibeau, S.; Gillispie, V.; Jacobs, M.; Theall, K.P. A Review of Herbal and Pharmaceutical Galactagogues for Breast-Feeding. *Ochsner J.* **2016**, *16*, 511-524.

11. Sim, T.F.; Sherriff, J.L.; Hattingh, H.L.; Parsons, R.; Tee, L.B. The use of herbal medicines during breastfeeding: A population-based survey in Western Australia. *BMC Complement. Altern. Med.* **2013**, *13*, 317. [CrossRef]

12. Abascal, K.; Yarnell, E. Botanical Galactagogues. *Altern. Complement. Ther.* **2008**, *14*, 288-294. [CrossRef]

13. Forinash, A.B.; Yancey, A.M.; Barnes, K.N.; Myles, T.D. The Use of Galactogogues in the Breastfeeding Mother. *Ann. Pharmacother.* **2012**, *46*, 1392-1404. [CrossRef]

14. Sevrin, T.; Alexandre-Gouabau, M.-C.; Darmaun, M.; Palvadeau, A.; André, A.; Nguyen, P.; Ouguerram, K.; Boquien, C.-Y. Use of water turnover method to measure mother's milk flow in a rat model: Application to dams receiving a low protein diet during gestation and lactation. *PLoS ONE* **2017**, *12*, e0180550. [CrossRef]

15. Sevrin, T.; Alexandre-Gouabau, M.-C.; Castellano, B.; Aguesse, A.; Ouguerram, K.; Nguyen, P.; Darmaun, D.; Boquien, C.-Y.; Gouabau, A. Impact of Fenugreek on Milk Production in Rodent Models of Lactation Challenge. *Nutrients* **2019**, *11*, 2571. [CrossRef] [PubMed]

16. Foong, S.C.; Tan, M.L.; Foong, W.C.; A Marasco, L.; Ho, J.J.; Ong, J.H. Oral galactagogues (natural therapies or drugs) for increasing breast milk production in mothers of non-hospitalised term infants. *Cochrane Database Syst. Rev.* **2020**. [CrossRef]

17. Truchet, S.; Honvo-Houéto, E. Physiology of milk secretion. *Best Pract. Res. Clin. Endocrinol. Metab.* **2017**, *31*, 367-384. [CrossRef]

18. E Grzeskowiak, L.; Wlodek, M.E.; Geddes, D.T. What Evidence Do We Have for Pharmaceutical Galactagogues in the Treatment of Lactation Insufficiency?-A Narrative Review. *Nutrients* **2019**, *11*, 974. [CrossRef]

19. Osorio, J.S.; Lohakare, J.; Bionaz, M. Biosynthesis of milk fat, protein, and lactose: Roles of transcriptional and posttranscriptional regulation. *Physiol. Genom.* **2016**, *48*, 231-256. [CrossRef]

20. Vernon, R.G. Lipid metabolism during lactation: A review of adipose tissue-liver interactions and the development of fatty liver. *J. Dairy Res.* **2005**, *72*, 460-469. [CrossRef]

21. Hull, K.; Harvey, S.; Hull, K.L. Growth hormone: Roles in female reproduction. *J. Endocrinol.* **2001**, *168*, 1-23. [CrossRef]

22. Anderson, S.M.; MacLean, P.S.; McManaman, J.L.; Neville, M.C. Lactation and its hormonal control. In *Knobil and Neill's Physiology of Reproduction*, 4th ed.; Elsevier: New York, NY, USA, 2015; pp. 2055-2105.

23. Bahmani, M.; Shirzad, H.; Mirhosseini, M.; Mesripour, A.; Rafieian-Kopaei, M. A Review on Ethnobotanical and Therapeutic Uses of Fenugreek (Trigonella foenum-graceum L). *J. Evid.-Based Integr. Med.* **2015**, *21*, 53-62. [CrossRef]
24. Fuller, S.; Stephens, J.M. Diosgenin, 4-Hydroxyisoleucine, and Fiber from Fenugreek: Mechanisms of Actions and Potential Effects on Metabolic Syndrome. *Adv. Nutr.* **2015**, *6*, 189–197. [CrossRef] [PubMed]

25. Wani, S.A.; Kumar, P. Fenugreek: A review on its nutraceutical properties and utilization in various food products. *J. Saudi Soc. Agric. Sci.* **2018**, *17*, 97–106. [CrossRef]

26. Abdou, R.M.; Fathey, M. Evaluation of early postpartum fenugreek supplementation on expressed breast milk volume and prolactin levels variation. *Egypt. Pediatr. Assoc. Gaz.* **2018**, *66*, 57–60. [CrossRef]

27. Anju, V.S.; Sreeharshan, S. In vitro estrogenic activities of fenugreek Trigonella foenum graecum seeds. *Indian J. Med. Res.* **2010**, *131*, 814–819.

28. Al-Chalabii, I.M.S. Diosgenin effect on rats and mice ovaries. *J. Fac. Med. Baghdad* **2005**, *47*, 296–301.

29. Kandhare, A.D.; Bandyopadhyay, D.; Thakurdesai, P.A. Low molecular weight galactomannans-based standardized fenugreek seed extract ameliorates high-fat diet-induced obesity in mice via modulation of FASn, IL-6, leptin, and TRIP-Br2. *RSC Adv.* **2018**, *8*, 32401–32416. [CrossRef]

30. Shim, S.H.; Lee, E.J.; Kim, J.S.; Kang, S.S.; Ha, H.; Lee, H.Y.; Kim, C.; Lee, J.-H.; Son, K.H. Rat Growth-Hormone Release Stimulators from Fenugreek Seeds. *Chem. Biodivers.* **2008**, *5*, 1753–1761. [CrossRef]

31. Koetsier, G.; Cantor, E. *A Practicle Guide to Analysing Nucleic Acid Concentration and Purity with Microvolume Spectrophotometer*; New England Biolabs Inc.: Ipswich, MA, USA, 2019.

32. Rudolph, M.C.; McManaman, J.L.; Hunter, L.; Phang, T.; Neville, M.C. Functional Development of the Mammary Gland: Use of Expression Profiling and Trajectory Clustering to Reveal Changes in Gene Expression During Pregnancy, Lactation, and Involution. *J. Mammary Gland Biol. Neoplasia* **2003**, *8*, 287–307. [CrossRef]

33. Rudolph, M.C.; McManaman, J.L.; Phang, T.; Russell, T.; Kominsky, D.J.; Serkova, N.J.; Stein, T.; Anderson, S.M.; Neville, M.C. Metabolic regulation in the lactating mammary gland: A lipid synthesizing machine. *Physiol. Genom.* **2007**, *28*, 323–336. [CrossRef]

34. Marshall, O.J. PerlPrimer: Cross-platform, graphical primer design for standard, bisulphite and real-time PCR. *Bioinformatics* **2004**, *20*, 2471–2472. [CrossRef] [PubMed]

35. Kasperczyk, A.; Machnik, G.; Dobrakowski, M.; Sypniewski, D.; Birkner, E.; Kasperczyk, S. Gene expression and activity of antioxidant enzymes in the blood cells of workers who were occupationally exposed to lead. *Toxicology* **2012**, *301*, 79–84. [CrossRef] [PubMed]

36. Bionaz, M.; Loor, J.J. Identification of reference genes for quantitative real-time PCR in the bovine mammary gland during the lactation cycle. *Physiol. Genom.* **2007**, *29*, 312–319. [CrossRef] [PubMed]

37. Kaur, R.; Sodhi, M.; Sharma, A.; Sharma, V.L.; Verma, P.; Swami, S.K.; Kumari, P.; Mukesh, M. Selection of suitable reference genes for normalization of quantitative RT-PCR (RT-qPCR) expression data across twelve tissues of riverine buffaloes (Bubalus bubalis). *PLoS ONE* **2018**, *13*, e0191558. [CrossRef] [PubMed]

38. Agnoux, A.M.; Antignac, J.-P.; Boquien, C.-Y.; David, A.; Desnots, E.; Ferchaud-Roucher, V.; Darmaun, D.; Parnet, P.; Alexandre-Gouabau, M.-C. Perinatal protein restriction affects milk free amino acid and fatty acid profile in lactating rats: Potential role on pup growth and metabolic status. *J. Nutr. Biochem.* **2015**, *26*, 784–795. [CrossRef]

39. Kim, H.-Y. Statistical notes for clinical researchers: post-hoc multiple comparisons. *Restor. Dent. Endod.* **2015**, *40*, 172–176. [CrossRef]

40. Hellemans, J.; Vandesompele, J. qPCR data analysis – unlocking the secret to successful results. In *PCR Troubleshooting and Optimization: The Essential Guide*; Caister Academic Press: Edinburgh, UK, 2011.

41. Yuan, J.S.; Reed, A.; Chen, F.; Stewart, J.C.N. Statistical analysis of real-time PCR data. *BMC Bioinform.* **2006**, *7*, 85. [CrossRef]

42. McManaman, J.L. Formation of milk lipids: A molecular perspective. *Clin. Lipidol.* **2009**, *4*, 391–401. [CrossRef]

43. Modha, G.; Blanchard, A.; Iwasiow, B.; Mao, X.J.; Troup, S.; Adeyinka, A.; Watson, P.; Shiu, R.; Myal, Y. Developmental Changes in Insulin-like Growth Factor I Receptor Gene Expression in the Mouse Mammary Gland. *Endocr. Res.* **2004**, *30*, 127–140. [CrossRef]

44. Oppi-Williams, C.; Suagee, J.; Corl, B. Regulation of lipid synthesis by liver X receptor α and sterol regulatory element-binding protein 1 in mammary epithelial cells. *J. Dairy Sci.* **2013**, *96*, 112–121. [CrossRef]

45. Bionaz, M.; Loor, J. Gene networks driving bovine milk fat synthesis during the lactation cycle. *BMC Genom.* **2008**, *9*, 366. [CrossRef] [PubMed]

46. Soloff, M.S. Oxytocin Receptors and Mammary Myoepithelial Cells. *J. Dairy Sci.* **1982**, *65*, 326–337. [CrossRef]
47. Bionaz, M.; Loor, J.J. Gene Networks Driving Bovine Mammary Protein Synthesis during the Lactation cycle. *Bioinform. Biol. Insights* **2011**, *5*, BBI-S7003. [CrossRef]

48. Jones, R.G.; Ilic, V.; Williamson, D.H. Physiological significance of altered insulin metabolism in the conscious rat during lactation. *Biochem. J.* **1984**, *220*, 455–460. [CrossRef] [PubMed]

49. Brogan, R.S.; Mitchell, S.E.; Trayhurn, P.; Smith, M.S. Suppression of Leptin During Lactation: Contribution of the Suckling Stimulus Versus Milk Production 1. *Endocrinology* **1999**, *140*, 2621–2627. [CrossRef]

50. Walsh, D.S.; Vesely, J.; Mahadevan, S. Relationship between Milk Production and Circulating Hormones in Dairy Cows. *J. Dairy Sci.* **1980**, *63*, 290–294. [CrossRef]

51. Ošt'ádalová, I.; Babický, A. Periodization of the early postnatal development in the rat with particular attention to the weaning period. *Physiol. Res.* **2012**, *61*, S1–S7.

52. Al-Sherwany, D.A.O. Feeding effects of fenugreek seeds on intake, milk yield, chemical composition of milk and some biochemical parameters in Hamdani ewes. *Al-Anbar J. Vet. Sci.* **2015**, *8*, 49–54.

53. Mahgoub, A.A.; Sallam, M.T. Effect of Extract Crushed Fenugreek Seeds as Feed Additive on some Blood Parameters, Milk Yield and Its Composition of Lactating Egyptian Buffaloes. *J. Anim. Poult. Prod.* **2016**, *7*, 269–273. [CrossRef]

54. Alamer, M.A.; Basiouni, G.F. Feeding effect of Fenugreek seed (Trigonella foenum-graecum L.) on lactation performance, some plasma constituent and growth hormone level in goats. *Pak. J. Biol. Sci.* **2005**, *8*, 1553–1556.

55. Liu, H.; Hua, Y.; Luo, H.; Shen, Z.; Tao, X.; Zhu, X. An Herbal Galactagogue Mixture Increases Milk Production and Aquaporin Protein Expression in the Mammary Glands of Lactating Rats. *Evid.-Based Complement. Altern. Med.* **2015**, *2015*, 1–6. [CrossRef] [PubMed]

56. Aguirre, G.A.; Ita, J.R.-D.; De La Garza, R.G.; Castilla-Cortázar, I. Insulin-like growth factor-1 deficiency and metabolic syndrome. *J. Transl. Med.* **2016**, *14*, 1–23. [CrossRef]

57. Qiu, H.; Yang, J.-K.; Chen, C. Influence of insulin on growth hormone secretion, level and growth hormone signalling. *Sheng Li Xue Bao Acta Physiol. Sin.* **2017**, *69*, 541–556.

58. Birzniece, V.; Sata, A.; Ho, K.K.Y. Growth hormone receptor modulators. *Rev. Endocr. Metab. Disord.* **2008**, *10*, 145–156. [CrossRef]

59. Chi, M.M.-Y.; Schlein, A.L.; Moley, K.H. High Insulin-Like Growth Factor 1 (IGF-1) and Insulin Concentrations Trigger Apoptosis in the Mouse Blastocyst via Down-Regulation of the IGF-1 Receptor\*This work was supported in part by NIH through Grants RO1-HD-38061-01A1 (to K.H.M.), P60-DK-30579 (to K.H.M.), and the Washington University Clinical Nutrition Research Unit Center Grant P30-DK-56341 (to K.H.M.); and by the Burroughs Wellcome Fund through a Career Award in the Biomedical Sciences (to K.H.M.). *Endocrinology* **2000**, *141*, 4784–4792. [CrossRef]

60. Payankaulam, S.; Raicu, A.-M.; Arnosti, D.N. Transcriptional Regulation of INSR, the Insulin Receptor Gene. *Genes* **2019**, *10*, 984. [CrossRef]

61. Rhodes, C.J.; White, M.F.; Leahy, J.L.; Kahn, S.E. Direct Autocrine Action of Insulin on β-Cells: Does It Make Physiological Sense? *Diabetes* **2013**, *62*, 2157–2163. [CrossRef]

62. Akers, R.M. Major Advances Associated with Hormone and Growth Factor Regulation of Mammary Growth and Lactation in Dairy Cows. *J. Dairy Sci.* **2006**, *89*, 1222–1234. [CrossRef]

63. Steffensen, K.R.; Gustafsson, J.-Å. Putative Metabolic Effects of the Liver X Receptor (LXR). *Diabetes* **2004**, *53*, S36–S42. [CrossRef]

64. Cruz-Garcia, L.; Sanchez-Gurmaches, J.; Gutiérrez, J.; Navarro, I.; Alvarez, I.N. Role of LXR in trout adipocytes: Target genes, hormonal regulation, adipocyte differentiation and relation to lipolysis. *Comp. Biochem. Physiol. Part A: Mol. Integr. Physiol.* **2012**, *163*, 120–126. [CrossRef] [PubMed]

65. Bionaz, M.; Chen, S.; Khan, M.J.; Loor, J. Functional Role of PPARs in Ruminants: Potential Targets for Fine-Tuning Metabolism during Growth and Lactation. *PPAR Res.* **2013**, *2013*, 1–28. [CrossRef] [PubMed]

66. Li, J.; Luo, J.; Zhu, J.; Sun, Y.; Yao, D.; Shi, H.; Wang, W. Regulation of the fatty acid synthase promoter by liver X receptor α through direct and indirect mechanisms in goat mammary epithelial cells. *Comp. Biochem. Physiol. Part B Biochem. Mol. Biol.* **2015**, *184*, 44–51. [CrossRef] [PubMed]

67. Wehrenberg, W.B.; Gaillard, R.C. Neuroendocrine Mechanisms Regulating Growth Hormone and Prolactin Secretion during Lactation. *Endocrinology* **1989**, *124*, 464–467. [CrossRef]

68. Christoffel, J.; Rimoldi, G.; Wuttke, W. Effects of 8-prenylnaringenin on the hypothalamo-pituitary-uterine axis in rats after 3-month treatment. *J. Endocrinol.* **2006**, *188*, 397–405. [CrossRef]
69. Lee, Y.-R.; Park, J.; Yu, H.-N.; Kim, J.-S.; Youn, H.J.; Jung, S.H. Up-regulation of PI3K/Akt signaling by 17β-estradiol through activation of estrogen receptor-α, but not estrogen receptor-β, and stimulates cell growth in breast cancer cells. *Biochem. Biophys. Res. Commun.* **2005**, *336*, 1221–1226. [CrossRef]

70. Degirmencioglu, T.; Unal, H.; Ozbilgin, S.; Kuraloglu, H. Effect of ground fenugreek seeds (Trigonella foenum-graecum) on feed consumption and milk performance in Anatolian water buffaloes. *Arch. Anim. Breed.* **2016**, *59*, 345–349. [CrossRef]

71. El-Nor, S.A.; Khattab, H.; Al-Alamy, H.; Salem, F.; Abdo, M. Effect of Some Medicinal Plants Seeds in the Rations on the Productive Performance of Lactating Buffaloes. *Int. J. Dairy Sci.* **2007**, *2*, 348–355. [CrossRef]

72. Hirscbberg, A.L. Hormonal regulation of appetite and food intake. *Ann. Med.* **1998**, *30*, 7–20. [CrossRef]

73. Hanwell, A.; Malcolm, P. Physiological effects of lactation on the mother. *Symp. Zool. Soc. Lond.* **1977**, *41*, 297–312.

74. Petit, P.; Sauvaire, Y.; Ponsin, G.; Manteghetti, M.; Fave, A.; Ribes, G. Effects of a fenugreek seed extract on feeding behaviour in the rat: Metabolic-endocrine correlates. *Pharmacol. Biochem. Behav.* **1993**, *45*, 369–374. [CrossRef]

75. Petit, P.R.; Sauvaire, Y.D.; Hillaire-Buys, D.M.; Leconte, O.M.; Baissac, Y.G.; Ponsin, G.R.; Ribes, G.R. Steroid saponins from fenugreek seeds: Extraction, purification, and pharmacological investigation on feeding behavior and plasma cholesterol. *Steroids* **1995**, *60*, 674–680. [CrossRef]

76. Matsushita, N.; Kato, Y.; Shimatsu, A.; Katakami, H.; Yanaihara, N.; Imura, H. Effects of VIP, TRH, GABA and dopamine on prolactin release from superfused rat anterior pituitary cells. *Life Sci.* **1983**, *32*, 1263–1269. [CrossRef]

77. Xu, J.; Kirigiti, M.A.; Grove, K.L.; Smith, M.S. Regulation of Food Intake and Gonadotropin-Releasing Hormone/Luteinizing Hormone during Lactation: Role of Insulin and Leptin. *Endocrinology* **2009**, *150*, 4231–4240. [CrossRef]

78. Morag, M.; Popliker, F.; Yagil, R. Effect of litter size on milk yield in the rat. *Lab. Anim.* **1975**, *9*, 43–47. [CrossRef] [PubMed]

79. Inui, A.; Tsai, M.; Amitani, H.; Asakawa, A. Stimulation of leptin secretion by insulin. *Indian J. Endocrinol. Metab.* **2012**, *16*, 543–548. [CrossRef]

80. Pliquett, R.; Fuhrer, D.; Falk, S.; Zysset, S.; Von Cramon, D.; Stumvoll, M. The Effects of Insulin on the Central Nervous System—Focus on Appetite Regulation. *Horm. Metab. Res.* **2006**, *38*, 442–446. [CrossRef] [PubMed]

81. Denis, R.G.P.; Williams, G.; Vernon, R.G. Regulation of serum leptin and its role in the hyperphagia of lactation in the rat. *J. Endocrinol.* **2003**, *176*, 193–203. [CrossRef]

<figure><img src="image_8.png" /></figure>

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

<footnote>**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.</footnote>