Skip to main content

Serum metabolomics study of the association between dairy intake and the anti-müllerian hormone annual decline rate



Dairy intake has been implicated in later ovarian aging but mechanism underlying the association is unknown. This study aimed to investigate (1) associations between dairy intake and metabolites previously shown related to anti-müllerian hormone (AMH) decline rate; (2) mediating roles of these metabolites in the prospective association of total dairy consumption with odds of AMH fast decline rate.


The participants comprised 186 reproductive-aged women randomly selected from the Tehran Lipid and Glucose Study. AMH was measured at baseline (1999–2001) and the 5th follow-up (2014–2017), and dietary data was collected at the second follow-up (2005–2008) using a food frequency questionnaire. Untargeted metabolomics was performed by gas chromatography–mass spectrometry using fasting-serum samples of the second follow-up. We analyzed dairy intake in association with the eight metabolites linked to the higher odds of AMH fast decline rate using linear regression with the Benjamini–Hochberg false discovery correction. Mediatory roles of the metabolites were assessed by bootstrapping.


Mean age and BMI of the participants at metabolomics assessment were 44.7 ± 5.87 years and 28.8 ± 4.88 kg/m2, respectively. Phosphate, branched-chain amino acids (BCAAs), and proline decreased significantly from the first to the third tertile of total dairy intake. Total dairy as a continuous variable inversely associated with phosphate (beta = −0.166; p value = 0.018), valine (beta = −0.176; p value = 0.016), leucine (beta = −0.226; p value = 0.002), proline (beta = −0.219; p value = 0.003), and urea (beta = −0.156; p = 0.035) after accounting for all potential covariates and correction for multiplicity (q-value < 0.1). Fermented dairy showed similar results, but milk did not associate with any of the metabolites. Simple mediation showed significant indirect effects for phosphate, proline, and BCAAs but not urea. Entering the sum of phosphate, proline, and BCAAs as a mediator, the metabolites' total indirect effects were significant [β = −0.12 (95% CIs − 0.26, − 0.04)]. In contrast, the direct association of total dairy intake with the fast decline in AMH was non-significant [β = −0.28 (95% CIs − 0.67, 0.10)].


Total dairy was inversely associated with AMH decline rate-related metabolites. Inverse association of dairy intakes with the odds of AMH fast decline rate was indirectly mediated by lower phosphate, proline, and BCAAs.


During women's reproductive lifespan, fertility reduces with the gradual exhaustion of ovarian follicular pools till menopause occurs [1]. The reduction of this ovarian reserve is highly age-dependent, but the depletion rate varies among chronological age individuals [2,3,4]. Environmental factors, including nutrition, are thought to influence reproductive aging [5, 6]. A prospective study from the Nurses' Health Study (NHS) showed modestly later timing of menopause in women aged < 51 years with higher consumption of low-fat dairy and skim milk [7]. Higher intakes of low-fat dairy and skim milk in another prospective study from NHS ӀӀ, also related to the lower risk of early menopause [8]. Recently, we reported inverse prospective associations between total dairy, milk, and fermented dairy intakes and odds of anti-müllerian hormone (AMH) fast decline rate, as assessed by two-time point measurements of AMH during the study [9]. Since AMH is the best biomarker of ovarian reserve [10, 11], our findings reinforce the possible association of dairy intakes with ovarian aging.

Metabolomics is an emerging tool that can improve the understanding of the biological process, metabolic pathway related to nutrition, identifying biomarkers for diseases or specific exposures, including nutrition [12, 13]. We conducted a prospective study on associations of metabolites quantified by untargeted metabolomics approach using gas chromatography–mass spectrometry (GC–MS) with AMH decline rate to understand the metabolic pathway related to ovarian reserve reduction. The investigation suggests eight metabolites that are positively related to the odds of AMH fast decline rate [14]. In the present study, we aimed to assess the cross-sectional associations between dairy intakes and the eight metabolites related to the AMH fast decline rate. Furthermore, we aimed to explore whether the AMH decline rate-related metabolites mediate the association between dairy intake and AMH decline rate. Therefore, mediation analysis was used to evaluate the direct/indirect associations between total dairy intake and odds of AMH fast decline rate, considering the metabolites as mediators.



We conducted this research in the Tehran Lipid and Glucose Study (TLGS). The TLGS is a population-based cohort study in which 15,005 individuals aged ≥ 3 years recruited between 1999 and 2001 from residents of district 13 of Tehran, Iran's capital. Data on lifestyle factors and health of the participants have been reassessed every 3 years after the baseline survey [15, 16].

AMH concentrations of 1015 women aged 20–50 years who met eligibility criteria were measured at baseline. The eligibility criteria were having regular menstrual cycles and proven natural fertility, not having a history of hormonal disorders, i.e., polycystic ovarian syndrome, hysterectomy, or ovarian surgery. Among which 245 women were randomly selected to re-measure AMH concentrations at the 5th follow-up examination (2014–2017). As previously reported, the demographic and anthropometric characteristics of the 245 women, except for age that was higher, were similar to the other women who did not include in this study  [9]. Metabolomics analysis for 186 randomly selected women were done after excluding those with a history of hysterectomy or ovarian surgery (n = 7) during the follow-up and those without dietary data (n = 11). The random case selection function in SPSS was used to determine the participants in each step mentioned above. More details for the participants' selection were reported elsewhere [14]. The Ethics committee of the Research Institute for Endocrine Sciences at Shahid Beheshti University of Medical Sciences (IR.SBMU.ENDOCRINE.REC.A396.534) approved the study, and written informed consent was obtained from all the participants.

AMH annual decline rate

Serum AMH concentrations were measured at the baseline and 5th follow-up examination by enzyme immunoassay (EIA) with Gen II kit (Beckman Coulter, Inc., CA, USA) and the Sunrise ELISA reader (Tecan Co. Salzburg, Austria). AMH decline rate was determined by subtracting the AMH values at the 5th follow-up examination from the baseline values, divided by the baseline values. This total AMH decline rate was divided by the duration of follow-up (in years) to estimate the AMH annual decline rate and expressed in percentages. The AMH fast decline rate was defined as yearly decline rates ≥ of 5.93% based on this variable's third tertile.

Metabolomics analysis

Untargeted metabolomics analysis was conducted using the serum samples collected during the second follow-up examination (2005–2008). Sample preparation and derivatization, data acquisition, and data processing have been previously reported in detail. [14] In brief, iced methanol was added to 100 µl serum sample to extract metabolites and precipitate proteins. After drying the supernatant by blowing the pure nitrogen, oxiamtion with 50 µg methoxyamine hydrochloride in pyridine and then trimethylsilylation with 50 ml of N-methyl-N-trimethylsilyltrifluoroacetamide (MSTFA) with 1% trimethylchlorosilane were performed. The derivatized sample was injected in splitless mode into the DB-5MS column of gas chromatography to separate the metabolites. The single quadrupole mass spectrometer was used to detect mass using the electron impact (EI) ionization technique at 70 eV. Mass spectra were recorded at 0.4 scans/s with a range of m/z 45–450. Raw data were processed using the mass spectrometry-data independent analysis software (MS-DIAL, v.3.98). Metabolites were annotated using accurate mass and retention index (RI) with Mass Bank of North America (MoNA;, and manually checked by National Institute of Standards and Technology (NIST) library, the Golm Metabolome Database (GMD), and the Human Metabolome Database (HMDB). Peak height using the quantifier mass ion was used for quantification, which was normalized using succinic acid d-4 as an internal standard (peak height metabolite/peak height succinic acid d-4); the relative peak intensity was used in the statistical analyses. We previously found eight metabolites positively associated with odds of AMH fast decline rate, including N-Acetyl-d-glucosamine, phosphate, valine, leucine, isoleucine, proline, pyroglutamic acid, and urea [14]. In this study, we assessed the associations between the AMH decline rate-related metabolites and dairy intakes.

Dietary assessment

Dietary data that were collected during the second follow-up (2005–2008) examination by the valid food frequency questionnaire (FFQ) was used in the current study. Daily dietary intakes of participants during the past year were estimated based on the frequency and amount of the 168 food items of the questionnaire. “low-fat milk”, “high-fat milk”, and “chocolate milk” were food items included for milk consumption. “Yoghurt”, “cheese”, “dough”, “kashk”, and “ice-cream” were food items included for estimated intakes of fermented dairy. Total dairy intakes were the sum of milk and fermented dairy consumptions. Intakes of these dairy groups were calculated as grams per energy intake of 1000 kcal (g/1000 kcal). A validity study suggested correlation coefficients of 0.59 in women for dairy products intakes assessed by the FFQ and multiple dietary recalls [17].

Other variables

Demographic and reproductive information was gathered by interview using a predefined questionnaire. Anthropometric variables were measurements by trained staff at the TLGS unit, and body mass index (BMI) was calculated as weight (kg)/height (m2). Blood pressure was measured twice on the right arm with a mercury sphygmomanometer, and the mean of the two measures was used for the participants' blood pressure. Biochemical analyses were done using overnight fasting serum samples. Fasting blood glucose (FBG), triglycerides, total cholesterol, and high-density lipoprotein cholesterol (HDL-C) were measured by the enzymatic colorimetric method using Pars Azmon kits (Pars Azmon Inc., Tehran, Iran). In this study, data collected during the second follow-up examination (2005–2008) were used as covariates.

Statistical analysis

Participants' characteristics at the second follow-up examination, when metabolomics analyses were done, were compared by the tertiles of total dairy intakes using analysis of variance (ANOVA) for normally distributed variables, the Kruskal–Wallis test for non-normally distributed variables, and chi-square for categorical variables.

Relative peak intensities of metabolites were transformed into the natural log (Ln) or square root and then autoscaled before statistical analyses. Relative peak intensities of metabolites across tertiles of dairy groups were compared using the general linear model (GLM), and results were reported as differences in least-square means considering the first tertile as a reference group. Association between metabolites (as continuous variables) and dairy groups were also examined using linear regression. Statistical analyses were performed using IBM SPSS Statistics (version 20; IBM Corp, Armonk, NY). The Benjamini–Hochberg procedure was applied to reduce the false discovery rate (FDR) due to multiplicity [18]. p value ≤ 0.05 and FDR-adjusted q-value ≤ 0.2 were defined as statistically significant.

Those metabolites were identified to be associated with dairy intake were selected. The mediating effects of the metabolites on the association between total dairy intakes and AMH fast decline rate were examined by PROCESS for SPSS macro (Version 3.5). The dependent variable was AMH fast decline (as a dichotomous variable), the independent was total dairy intakes, and mediators were the metabolites. The indirect effect was statistically tested with the bootstrap method based on 5000 simulations and 95% confidence intervals (CIs). The indirect effect of each metabolite was investigated separately; since the metabolites were highly correlated, we could not find their independent indirect (mediation) effects using multiple mediation analysis. Therefore, the sum of phosphate, proline, and BCAAs was calculated and used as the mediator to estimate the total indirect effects of the 5 metabolites. Standardized β (95% CIs) were reported for direct and indirect effects. In all statistical analyses, age, BMI, dietary intakes of energy, grains (g/1000 kcal), and legumes (g/1000 kcal) at the time of metabolomics analysis were accounted for as covariates.


Participants’ characteristics

At the time of metabolomics analysis, the mean age and BMI of women were 44.7 ± 5.87 years and 28.8 ± 4.55 kg/m2, respectively, and 2.2% were smokers. Characteristics of the participants by tertiles of total dairy intake are shown in Table 1. No significant difference in age, BMI, waist circumference, and biochemical measurements were observed across the tertile of total dairy intake; the baseline value of AMH and the AMH annual decline rate were not significantly different across the tertiles. Furthermore, the frequency of individuals with diabetes did not significantly differ across tertiles of total dairy intake at metabolomics assessment. Mean intake of protein was higher, and grains were lower in women with higher intakes of total dairy intakes; the other dietary variables did not differ significantly.

Table 1 Characteristics of the participants by tertiles of total dairy intakes at metabolomics assessmenta

Associations between dairy groups and AMH decline rate-related metabolites

Least-squares means of 8 metabolites related to AMH fast decline rate across tertile categories of dairy groups are shown in Table 2. After adjusting for potential covariates including age, BMI, energy intake, grains and legumes intakes, the relative abundance of phosphate, valine, leucine, isoleucine, and proline was significantly lower with higher intakes of total dairy intakes. The findings were similar for fermented dairy intakes except for isoleucine that was non-significant (p-trend = 0.09); the associations remained significant after correction for multiple testing. However, the concentrations of the metabolites did not differ significantly across the tertiles of milk intakes.

Table 2 Relative abundance of the 8 metabolites related to the AMH decline rate by tertiles of dairy intakea

Associations between dairy groups and AMH decline rate–related metabolites are shown in Table 3. All metabolites were inversely associated with total dairy intakes in unadjusted models. The associations remained significant after accounting for age and BMI of the participants and correction for multiplicity (FDR-adjusted q-value < 0.2). However, after entering other dietary variables as covariates, the associations remained significant for phosphate, valine, leucine, proline, and urea. Higher intakes of fermented dairy were associated with a lower abundance of valine, leucine, proline, and urea after adjusting for all potential covariates (model 2). In the age and BMI-adjusted model, milk intake was inversely associated with pyroglutamic acid. However, the association became non-significant after controlling for other dietary variables.

Table 3 Standardized β for associations between dairy intakes and AMH decline rate-related metabolitesa,b

Mediating analysis

The mediating effects of phosphate, BCAAs, proline, and urea, the metabolites associated with total dairy intakes based on linear regression, on the association between total dairy intakes and AMH fast decline rate are reported in Table 4. After adjusting for age, BMI, energy, grains, and legumes intakes, the indirect effects of phosphate, proline, and BCAAs were significant. The direct effects of total dairy intakes were not significant in any of the simple mediation models. The sum of phosphate, proline, and BCAAs was calculated and used as a mediator to estimate the total indirect effects of the 5 metabolites. Total dairy intake was inversely associated with the pooled estimation of the five metabolites [β = −0.93 (95% CIs −1.58, − 0.28; p = 0.005)], and the metabolites were positively related to fast decline in AMH [β = 0.13 (95% CIs 0.05, 0.22; p = 0.002)] after adjusting for covariates. The metabolites' total indirect effects were significant, while the direct effect of total dairy intake on the fast decline in AMH was non-significant [β = −0.28 (95% CIs − 0.67, 0.10); p = 0.152].

Table 4 Mediation effects of dairy intakes on AMH fast decline rate through metabolitesa,b


In this study, we examined cross-sectional associations between dairy products and metabolites previously found to be related to the higher odds of AMH fast decline rate. After adjusting for age, BMI, intakes of energy, grains, and legumes the intensities of five AMH decline rate-related metabolites including phosphate, valine, leucine, isoleucine, and proline reduced significantly from the first to the third tertile of total dairy intake. Total dairy intakes as a continuous variable were also inversely associated with phosphate, valine, leucine, proline, and urea. Using mediating analyses, we identified that phosphate, BCAAs, and proline fully mediated the prospective inverse association of total dairy intakes and odds of AMH fast decline rate.

Reduced fertility, end of fertility, menstrual cycle irregularity, and menopause are events that women experience during their reproductive life span due to the gradual loss of numbers of ovarian follicles with aging [1]. Though the ovarian aging process naturally occurs in all women, the rates of follicular loss and the age at which the reproductive events happen considerably vary [1,2,3,4]. As fast ovarian reserve decline rate and early menopause increases the risk of cardiovascular diseases, diabetes, neurological diseases, and mortality [19,20,21], understanding modifiable determinants of ovarian aging are clinically important for managing fertility and other health risks in women. Few studies investigated the association between dairy intake and timing of menopause [7, 8, 22, 23]; among them, two studies showed the later timing of menopause with higher intakes of low-fat dairy and skimmed milk [7, 8]. We previously examined the prospective association between dairy intake and the annual decline rate of AMH as the best biomarkers of ovarian reserve. Findings of the study suggested inverse associations of total dairy, fermented dairy, and milk intakes and odds of AMH fast decline rate [9]. To study molecular biology underlying the association, we evaluated the associations between dairy intake and AMH decline rate-related metabolites in the current study.

Finding objective markers to assess dairy consumption accurately is of interest to better understand this heterogeneous food group's health effects. Some metabolomics studies have been conducted to find biomarkers reflecting habitual dairy consumption in free-living individuals [24,25,26,27]. These studies have proposed numerous metabolites, the majority of which were lipids, related to total dairy or specific dairy foods intakes. However, no reliable and validated biomarkers for dairy intake have been identified to date [28]. It is unclear that the dairy-related metabolites are indicative of dietary intakes or dairy-induced changes in metabolism. Our findings suggest that higher dairy intake are associated with lower circulatory phosphate, BCAAs, and proline. These findings reflect the metabolic pathways through which dairy may associate with AMH decline rate which is beyond markers of its dietary intake.

Approximately 21% of the amino acid residue of dairy protein is BCAAs, and therefore the food group is an important dietary source of the essential amino acids [29]. Proline, a non-essential amino-acids, is also the second most abundant amino acid in a dairy protein with a proportion of 9%. The inverse association between dairy and serum BCAAs and proline was observed in the study, unlike that the food groups are one of the rich sources of these amino-acids. However, it is unclear whether dietary intakes of the amino-acids influence their fasting circulatory levels as the number of the studies is limited, and their findings are inconsistent [30,31,32,33]. A prospective study showed higher plasma levels of BCAAs with their higher dietary intakes [30]. A cross-sectional study also reported positive associations for total animal protein and red meat with plasma BCAAs, but dairy protein intake was unrelated to the blood levels of BCAAs [33]. No significant association was observed between dietary and fasting plasma levels of BCAAs in another cross-sectional study [31]. Short term changes in dairy intakes by 1.5–2 serving of daily milk intake for one month could not affect plasma BCAAs, proline, and other free amino acids in an intervention study of 102 healthy adults [32]. Whereas the BCAAs catabolism is tightly regulated, increased dietary intakes of BCAAs have been proposed not to contribute a remarkable elevation in their circulatory levels in prandial state [34, 35]. In addition to dietary intake, protein turnover and gene expression of activity of the branched-chain-keto acid dehydrogenase (BCKAD), the rate-limiting enzyme of BCAAs catabolism, also influence their circulatory levels [34, 36]. Although the cause and effect relationship between insulin resistance and elevated BCAA is unknown, it is believed that insulin activates BCKD, so insulin resistance leads to the decreased catabolism of the amino-acids [35]. Associations of dairy intakes with gene expression and activity of the enzyme have not been studied, but dairy intakes through improved insulin sensitivity [37] can be related to the lower circulatory BCAAs.

Type of milk, the microorganism, and process increase diversity of dairy products. The health effects of the diverse food groups may be different [38]. In this study, we categorized dairy products into two groups of fermented dairy and milk. Considering the subtypes, the inverse association between milk intake and AMH decline rate metabolites were not significant. In our participants, fermented dairy is a dominant type of dairy food consumed, and milk constitutes only 36.2% of total dairy intakes. This amount of milk consumption may be inadequate to find any significant relationship with the metabolites. Since the milk and fermented dairy could reduce the odds of AMH fast decline rate to almost the same extent [9], the fermentation process might not promote these dairy-metabolite associations.

In our study, total dairy intake could explain the variance of AMH decline rate-related metabolites modestly, ranging from 2.6% in phosphate to 6.7% in leucine, so dairy intake might not be the main variable influencing these metabolites. However, the mediating analyses showed that phosphate, proline, and BCAAs mediated the prospective association between total dairy intake and odds of AMH fast decline rate. In other words, higher dairy intake through the lower abundance of phosphate, proline, and BCAAs was associated with lower odds of AMH fast decline.

Our findings provided some biological insights supporting the hypothesis that dairy intake may reduce the AMH decline rate. However, the limitations of the study should be considered. First, habitual dairy intake was indicated using FFQ, which is prone to recall bias. Second, metabolomics analysis was performed using serum samples collected at 2005–2008 and stored at − 80. Though the stability of metabolites at this temperature is different, in general, the prolonged storage of samples has been shown to change the concentrations of the metabolites, which can be a potential source of bias in a data set [39, 40]. Third, we cannot conclude any causal effects of dairy intakes on AMH decline rates-related metabolites based on the cross-sectional analyses. Fourth, the metabolites were highly correlated, and therefore we could not examine their independent mediating roles using multiple mediation analysis. Fifth the annual AMH decline rate was calculated based on AMH values measured at the baseline and the fifth follow-up and assuming a fixed pattern for decline rate; while an overtime intra-individual variations may be proposed.


In conclusion, after accounting for age, BMI, and dietary intakes, higher total dairy consumption was associated with lower intensities of serum phosphate, proline, and BCAAs, and urea, metabolites that positively related to the AMH decline rate. Lower concentrations of phosphate, proline, and BCAAs indirectly mediated the prospective inverse association of dairy intakes and odds of AMH fast decline rate. The study's findings provide some knowledge on the metabolites that may be involved in the association between dairy intake and ovarian aging, which can be confirmed using a targeted metabolomics approach. Further investigation in other populations using a combination of metabolomics platforms increases the numbers of metabolites related to the dairy-ovarian aging association.

Availability of data and materials

The dataset analyzed during the current study are available from the corresponding author on upon request.



Analysis of variance


Anti-Müllerian hormone


Branched-chain amino acids


Electron impact


Enzyme immunoassay


Fasting blood glucose


False discovery rate


Food frequency questionnaire


Gas chromatography–mass spectrometry


General linear model


Golm Metabolome Database


High-density lipoprotein cholesterol


Human Metabolome Database


National Institute of Standards and Technology


Nurses' Health Study


Mass Bank of North America




Tehran Lipid and Glucose Study


  1. 1.

    Broekmans F, Soules M, Fauser B. Ovarian aging: mechanisms and clinical consequences. Endocr Rev. 2009;30(5):465–93.

    CAS  Article  Google Scholar 

  2. 2.

    Freeman EW, Sammel MD, Lin H, Boorman DW, Gracia CR. Contribution of the rate of change of antimullerian hormone in estimating time to menopause for late reproductive-age women. Fertil Steril. 2012;98(5):1254–9 (e1-2).

    CAS  Article  Google Scholar 

  3. 3.

    Gohari MR, Ramezani Tehrani F, Chenouri S, Solaymani-Dodaran M, Azizi F. Individualized predictions of time to menopause using multiple measurements of antimullerian hormone. Menopause. 2016;23(8):839–45.

    Article  Google Scholar 

  4. 4.

    de Kat AC, van der Schouw YT, Eijkemans MJ, Herber-Gast GC, Visser JA, Verschuren WM, et al. Back to the basics of ovarian aging: a population-based study on longitudinal anti-Mullerian hormone decline. BMC Med. 2016;14(1):151.

    Article  Google Scholar 

  5. 5.

    Richardson M, Guo M, Fauser B, Macklon N. Environmental and developmental origins of ovarian reserve. Hum Reprod Update. 2014;20(3):353–69.

    CAS  Article  Google Scholar 

  6. 6.

    Moslehi N, Mirmiran P, Tehrani FR, Azizi F. Current evidence on associations of nutritional factors with ovarian reserve and timing of menopause: a systematic review. Adv Nutr. 2017;8(4):597–612.

    PubMed  PubMed Central  Google Scholar 

  7. 7.

    Carwile JL, Willett WC, Michels KB. Consumption of low-fat dairy products may delay natural menopause. J Nutr. 2013;143(10):1642–50.

    CAS  Article  Google Scholar 

  8. 8.

    Purdue-Smithe AC, Whitcomb BW, Manson JE, Hankinson SE, Rosner BA, Troy LM, et al. A prospective study of dairy-food intake and early menopause. Am J Epidemiol. 2019;188(1):188–96.

    Article  Google Scholar 

  9. 9.

    Moslehi N, Mirmiran P, Azizi F, Tehrani FR. Do dietary intakes influence the rate of decline in anti-Mullerian hormone among eumenorrheic women? A population-based prospective investigation. Nutr J. 2019;18(1):83.

    CAS  Article  Google Scholar 

  10. 10.

    Jamil Z, Fatima SS, Ahmed K, Malik R. Anti-mullerian hormone: above and beyond conventional ovarian reserve markers. Dis Markers. 2016;2016:5246217.

    Article  Google Scholar 

  11. 11.

    Moolhuijsen LM, Visser JA. Anti-Müllerian hormone and ovarian reserve: update on assessing ovarian function. J Clin Endocrinol Metab. 2020;105(11):3361–73.

    Article  Google Scholar 

  12. 12.

    Peng B, Li H, Peng X-X. Functional metabolomics: from biomarker discovery to metabolome reprogramming. Protein Cell. 2015;6(9):628–37.

    CAS  Article  Google Scholar 

  13. 13.

    O’Gorman A, Brennan L. The role of metabolomics in determination of new dietary biomarkers. Proc Nutr Soc. 2017;76(3):295–302.

    CAS  Article  Google Scholar 

  14. 14.

    Moslehi N, Mirmiran P, Marzbani R, Rezadoost H, Mirzaie M, Azizi F, et al. Serum metabolomics study of women with different annual decline rates of anti-Müllerian hormone: an untargeted gas chromatography–mass spectrometry-based study. Hum Reprod. 2020;36(3):721–33.

    Article  Google Scholar 

  15. 15.

    Azizi F, Ghanbarian A, Momenan AA, Hadaegh F, Mirmiran P, Hedayati M, et al. Prevention of non-communicable disease in a population in nutrition transition: Tehran Lipid and Glucose Study phase II. Trials. 2009;10(1):1–15.

    Article  Google Scholar 

  16. 16.

    Azizi F, Zadeh-Vakili A, Takyar M. Review of rationale, design, and initial findings: Tehran Lipid and Glucose Study. Int J Endocrinol Metab. 2018;16(4 Suppl):e84777.

    PubMed  PubMed Central  Google Scholar 

  17. 17.

    Esfahani FH, Asghari G, Mirmiran P, Azizi F. Reproducibility and relative validity of food group intake in a food frequency questionnaire developed for the Tehran Lipid and Glucose Study. J Epidemiol. 2010;20(2):150–8.

    Article  Google Scholar 

  18. 18.

    Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol. 1995;57(1):289–300.

    Google Scholar 

  19. 19.

    Shuster LT, Rhodes DJ, Gostout BS, Grossardt BR, Rocca WA. Premature menopause or early menopause: long-term health consequences. Maturitas. 2010;65(2):161–6.

    Article  Google Scholar 

  20. 20.

    Brand JS, Van Der Schouw YT, Onland-Moret NC, Sharp SJ, Ong KK, Khaw K-T, et al. Age at menopause, reproductive life span, and type 2 diabetes risk: results from the EPIC-InterAct study. Diabetes Care. 2013;36(4):1012–9.

    Article  Google Scholar 

  21. 21.

    De Kat AC, Verschuren WM, Eijkemans MJ, Broekmans FJ, Van Der Schouw YT. Anti-Müllerian hormone trajectories are associated with cardiovascular disease in women: results from the Doetinchem cohort study. Circulation. 2017;135(6):556–65.

    Article  Google Scholar 

  22. 22.

    Nagel G, Altenburg H-P, Nieters A, Boffetta P, Linseisen J. Reproductive and dietary determinants of the age at menopause in EPIC-Heidelberg. Maturitas. 2005;52(3–4):337–47.

    Article  Google Scholar 

  23. 23.

    Dunneram Y, Greenwood DC, Burley VJ, Cade JE. Dietary intake and age at natural menopause: results from the UK Women’s Cohort Study. J Epidemiol Community Health. 2018;72(8):733–40.

    Article  Google Scholar 

  24. 24.

    Zheng Y, Yu B, Alexander D, Steffen LM, Boerwinkle E. Human metabolome associates with dietary intake habits among African Americans in the atherosclerosis risk in communities study. Am J Epidemiol. 2014;179(12):1424–33.

    Article  Google Scholar 

  25. 25.

    Pallister T, Jennings A, Mohney RP, Yarand D, Mangino M, Cassidy A, et al. Characterizing blood metabolomics profiles associated with self-reported food intakes in female twins. PLoS ONE. 2016;11(6):e0158568.

    Article  Google Scholar 

  26. 26.

    Wang Y, Gapstur SM, Carter BD, Hartman TJ, Stevens VL, Gaudet MM, et al. Untargeted metabolomics identifies novel potential biomarkers of habitual food intake in a cross-sectional study of postmenopausal women. J Nutr. 2018;148(6):932–43.

    Article  Google Scholar 

  27. 27.

    Hruby A, Dennis C, Jacques PF. Dairy intake in 2 American adult cohorts associates with novel and known targeted and nontargeted circulating metabolites. J Nutr. 2020;150(5):1272–83.

    Article  Google Scholar 

  28. 28.

    Münger LH, Garcia-Aloy M, Vázquez-Fresno R, Gille D, Rosana ARR, Passerini A, et al. Biomarker of food intake for assessing the consumption of dairy and egg products. Genes Nutr. 2018;13(1):26.

    Article  Google Scholar 

  29. 29.

    Chartrand D, Da Silva MS, Julien P, Rudkowska I. Influence of amino acids in dairy products on glucose homeostasis: the clinical evidence. Can J Diabetes. 2017;41(3):329–37.

    Article  Google Scholar 

  30. 30.

    Zheng Y, Li Y, Qi Q, Hruby A, Manson JE, Willett WC, et al. Cumulative consumption of branched-chain amino acids and incidence of type 2 diabetes. Int J Epidemiol. 2016;45(5):1482–92.

    Article  Google Scholar 

  31. 31.

    Jennings A, MacGregor A, Pallister T, Spector T, Cassidy A. Associations between branched chain amino acid intake and biomarkers of adiposity and cardiometabolic health independent of genetic factors: a twin study. Int J Cardiol. 2016;223:992–8.

    Article  Google Scholar 

  32. 32.

    Prodhan UK, Milan AM, Thorstensen EB, Barnett MP, Stewart RA, Benatar JR, et al. Altered dairy protein intake does not alter circulatory branched chain amino acids in healthy adults: a randomized controlled trial. Nutrients. 2018;10(10):1510.

    Article  Google Scholar 

  33. 33.

    Rousseau M, Guénard F, Garneau V, Allam-Ndoul B, Lemieux S, Pérusse L, et al. Associations between dietary protein sources, plasma BCAA and short-chain acylcarnitine levels in adults. Nutrients. 2019;11(1):173.

    CAS  Article  Google Scholar 

  34. 34.

    Siomkajło M, Daroszewski J. Branched chain amino acids: Passive biomarkers or the key to the pathogenesis of cardiometabolic diseases? Adv Clin Exp Med. 2019;28(9):1263–9.

    Article  Google Scholar 

  35. 35.

    Le Couteur DG, Solon-Biet SM, Cogger VC, Ribeiro R, de Cabo R, Raubenheimer D, et al. Branched chain amino acids, aging and age-related health. Ageing Res Rev. 2020;64:101198.

    Article  Google Scholar 

  36. 36.

    Adeva-Andany MM, López-Maside L, Donapetry-García C, Fernández-Fernández C, Sixto-Leal C. Enzymes involved in branched-chain amino acid metabolism in humans. Amino Acids. 2017;49(6):1005–28.

    CAS  Article  Google Scholar 

  37. 37.

    Sochol KM, Johns TS, Buttar RS, Randhawa L, Sanchez E, Gal M, et al. The effects of dairy intake on insulin resistance: a systematic review and meta-analysis of randomized clinical trials. Nutrients. 2019;11(9):2237.

    CAS  Article  Google Scholar 

  38. 38.

    Fernández M, Hudson JA, Korpela R, de los Reyes-Gavilán CG. Impact on human health of microorganisms present in fermented dairy products: an overview. BioMed Res Int. 2015;2015:412714.

    PubMed  PubMed Central  Google Scholar 

  39. 39.

    Haid M, Muschet C, Wahl S, Römisch-Margl W, Prehn C, Möller G, et al. Long-term stability of human plasma metabolites during storage at− 80 C. J Proteome Res. 2018;17(1):203–11.

    CAS  Article  Google Scholar 

  40. 40.

    Wagner-Golbs A, Neuber S, Kamlage B, Christiansen N, Bethan B, Rennefahrt U, et al. Effects of long-term storage at − 80 C on the human plasma metabolome. Metabolites. 2019;9(5):99.

    CAS  Article  Google Scholar 

Download references


Not applicable.


This study was supported by Shahid Beheshti University of Medical Sciences (Grant No. 13366-4).

Author information




NM: conceptual design of the study, performed experiments, metabolomics data processing, and statistical analyses, and wrote the manuscript. RM: performed experiments and metabolomics data processing. HR: provided experimental supervisory support and contributed to the critical revision of the manuscript. FRT: conceptual design of the study, helping in statistical analyses and interpretation, and contributed to the critical revision of the manuscript. PM and FA: conceptual design of the study, supervised the analyses and contributed to the critical revision of the manuscript. All authors read and approved the final version of manuscript.

Corresponding authors

Correspondence to Parvin Mirmiran or Fahimeh Ramezani Tehrani.

Ethics declarations

Ethics approval and consent to participate

The Ethics committee of the Research Institute for Endocrine Sciences at Shahid Beheshti University of Medical Sciences (IR.SBMU.ENDOCRINE.REC.A396.534) approved the study. Written informed consent was obtained from all the participants.

Consent for publication

Not applicable.

Competing interests

The authors have no conflicts of interest to declare.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Moslehi, N., Marzbani, R., Rezadoost, H. et al. Serum metabolomics study of the association between dairy intake and the anti-müllerian hormone annual decline rate. Nutr Metab (Lond) 18, 66 (2021).

Download citation


  • Anti-müllerian hormone
  • Branched-chain amino acids
  • Dairy
  • Mediation
  • Metabolomics
  • Ovarian reserve