Association of the functional ovarian reserve with serum metabolomic profiling by nuclear magnetic resonance spectroscopy: a cross-sectional study of ~ 400 women
BMC Medicine volume 18, Article number: 247 (2020)
Women with diminished ovarian reserve are known to have increased cardiovascular risk, whether there is a continuous association between the ovarian reserve biomarkers; anti-Müllerian hormone (AMH), antral follicle count (AFC) and cardio-metabolic risk factors are unknown.
A cross-sectional study of 398 women intending to undergo IVF with pre-treatment early follicular AMH and AFC measurements. Serum lipids, lipoprotein subclasses and low-molecular-weight metabolites were quantified by NMR spectroscopy (155 metabolic measures). Associations were analysed using multivariable regression.
Participants were mean 35.5 (SD 4.43) years old and had a median AMH of 16 pmol/l (IQR 8.8, 28.0 pmol/l) and a median AFC of 12 (IQR 7.16). AMH showed positive associations with HDL, omega-6 and polyunsaturated fatty acids and the amino acids isoleucine, leucine and tyrosine, with effects ranging from 0.11 (95%CI 0.004 to 0.21) for total lipids in small HDL to 0.16 (0.06 to 0.26) for isoleucine, for a mean difference of one SD of metabolite per one SD increment in AMH, and negatively with acetate: − 0.31(− 0.22, − 0.004) SD per 1 SD AMH. AFC was positively associated with alanine, glutamine and glycine. Results were consistent, though less precisely estimated, when restricted to those women who were preparing for treatment because of their partner’s infertility.
In women intending to have IVF, AMH and AFC were not associated with traditional lipid measured but were associated with a number of novel cardiovascular risk factors. Prospective studies will be required for replication, determination of causality and confirmation that ovarian reserve is impacting on metabolism rather than variation in metabolism is influencing ovarian reserve.
Female reproductive ageing is the result of a gradual decrease in both the quantity and quality of oocytes . Genetic, environmental and lifestyle factors are all known to contribute to the timing and depletion of the ovarian reserve. Markers of diminished ovarian reserve, such as low anti-Müllerian hormone (AMH) and low antral follicle count (AFC), have been shown to associate with earlier menopause [2, 3]. Several observational studies have investigated the association between these markers of diminished ovarian reserve and cardiovascular risk factors, with several [4,5,6,7,8] but not all studies [9,10,11], suggesting that a diminished ovarian reserve may be associated with an unfavourable circulating cardiometabolic risk profile and cardiovascular events. However, these studies have been limited to a restricted number of established cardiovascular risk factors, including evaluation of total and LDL cholesterol , homeostatic model assessment as a surrogate for insulin resistance (HOMA-IR)  or have considered these established risk factors together as a composite outcome of cardiometabolic risk [4, 5].
Detailed metabolic profiling or metabolomics has been applied successfully to identify novel biomarkers for the development of cardiovascular disease [13, 14] and all-cause mortality , with improved prediction as compared to models containing conventional risk factors . Serum nuclear magnetic resonance (NMR) metabolomics which enables reproducible quantification of circulating lipids and abundant metabolites  has facilitated its use in the assessment of the changes in metabolites with adiposity , glycaemia , pregnancy  and menopausal status .
The aim of the current study was to assess the association of ovarian reserve (as measured by AMH and AFC) with 155 circulating metabolic measures. These measures were profiled by a high-throughput NMR metabolomics platform, covering a wide range of metabolic pathways including lipoprotein lipids, fatty acids, amino acids, ketone bodies and glycaemic traits, which are highly relevant to cardiometabolic risk and overall health.
Study design and participants
This is a cross-sectional study of women aged 18 to 45 who presented at Glasgow Royal Infirmary, UK, for assessment prior to assisted conception between 1 April 2017 and 31 March 2019. Exclusion criteria were a documented positive pregnancy test at the time of presentation, body mass index (BMI) ≥ 35 kg/m2 and/or requiring oocyte or embryo donation. A total of 400 women were recruited, and of these, 398 (99%) had complete data on AMH, AFC and at least one NMR metabolite and were included in the analyses presented in this paper.
The study was conducted according to the ICH Guideline for good clinical practice, the Declaration of Helsinki and the Convention of the Council of Europe. All women provided written informed consent. The study protocol was approved prior to study initiation by the relevant institutional review boards (see Supplementary Material).
Demographic, lifestyle, fertility and medical history was obtained by self-reported questionnaire and clinical data by linkage to electronic medical records.
AFC was determined by two- or three-dimensional transvaginal ultrasound (AFC was defined as the total number of antral follicles with a size of 2–10 mm in both ovaries) on menstrual cycle days 2–4. Follicle counts greater than 20 follicles per ovary were classed as ≥ 20, consistent with the diagnostic threshold for PCOS, and that ovary was not counted further . Due to the potential for inter-sonographer variability, sonographers were provided with training and were asked to follow published practical recommendations for accurate transvaginal ultrasound ; each sonographer used standard equipment (Acuson Sequoia, Siemens Germany).
Non-fasted blood samples were collected during the same visit that transvaginal ultrasound was performed. AMH was measured using the Beckman Coulter AMH automated method on a clinically validated immunoassay platform (Access 2, Beckman Coulter, USA). The assay was calibrated and quality controlled using the manufacturer’s reagents and is known to have a measuring range of 0.08–24 ng/ml (0.57–171 pmol/l). The limit of quantitation (LoQ) was 0.02 ng/ml (0.014pmo/l), with the 20% CV LoQ 0.08 ng/ml (0.57 pmol/l). The coefficient of variation between runs for two levels of control ran at < 4.4%.
Additional blood samples were taken for NMR analyses and immediately spun and frozen at − 80 °C, and all NMR assays completed for this study were undertaken within 1 year of storage and with no previous freeze/thaw cycles.
Profiling of 155 lipid and metabolite measures was performed by a high-throughput targeted NMR platform [Nightingale Health© (Helsinki, Finland)] at the University of Bristol. The platform applies a single experimental setup, which allows for the simultaneous quantification of routine lipids, 14 lipoprotein subclasses and individual lipids transported by these particles, multiple fatty acids, glucose, the glycolysis precursors lactate and pyruvate, ketone bodies and amino acids in absolute concentration units (mostly mmol/l). The NMR-based metabolite quantification is achieved through measurements of three molecular windows from each sample. Two of the spectra (LIPO and LMWM windows) are acquired from native serum and one spectrum from serum lipid extracts (LIPID window). The NMR spectra were measured using Bruker AVANCE III spectrometer operating at 600 MHz. Measurements of native serum samples and serum lipid extracts are conducted at 37 °C and 22 °C, respectively. Details of this platform have been published previously [16, 23], and it has been widely applied in genetic and observational epidemiological studies [13,14,15, 17, 19, 20, 24,25,26]. Further details of the platform are provided in Additional file 1 (Supplemental text, Table S1 and Fig. S1 ).
Metabolite quantification and quality control
The NMR spectra were analysed for absolute metabolite quantification (molar concentration) in an automated fashion. For each metabolite, a ridge regression model was applied for quantification in order to overcome the problems of heavily overlapping spectral data. In the case of the lipoprotein lipid data, quantification models were calibrated using high-performance liquid chromatography methods and individually cross-validated against NMR-independent lipid data. Low-molecular-weight metabolites, as well as lipid extract measures, were quantified as millimoles per litre based on regression modelling calibrated against a set of manually fitted metabolite measures. The calibration data are quantified based on iterative line-shape fitting analysis using the PERCH NMR software (PERCH Solutions Ltd., Kuopio, Finland). Absolute quantification cannot be directly established for the lipid extract measures due to experimental variation in the lipid extraction protocol. Therefore, serum extract metabolites are scaled via the total cholesterol as quantified from the native serum LIPO spectrum.
Assessment of potential confounders
In relation to our analyses and control for confounders, we used currently recommended practice of defining confounders a priori (before undertaking analyses) using directed acyclic graphs (DAGs) [27, 28]. This approach defines confounders as any characteristic that is known to cause variation in the exposure (here ovarian biomarkers AMH and AFC) and outcome (NMR metabolites) or is plausibly a cause of exposure and outcome. Using a priori knowledge and relevant literature, our selected confounders were age, BMI, educational attainment, ethnicity, family history of cardiovascular disease (defined as first-degree relative affected) physical activity, alcohol intake, smoking status, duration and cause of infertility and whether infertility was primary or secondary (with secondary defined as a woman unable to establish a clinical pregnancy but who has previously been diagnosed with a clinical pregnancy) . These are known to or plausibly influence both ovarian reserve and the NMR metabolites . Figure 1a shows the DAG for this, our main, confounder-adjusted analyses. Following reviewer comments, we also considered the extent to which PCOS might be part of a confounding path and whether we should also adjust for it. There is evidence ovarian reserve and PCOS may share underlying common causes, including genetic variation and intrauterine exposures. If these, specific genetic variants and exposures are also related to cardiometabolic health, then PCOS may be on a confounding path between ovarian reserve and cardiometabolic health. Three possibilities are considered in Fig. 1b–d. Whether or not to adjust for PCOS is in part related to the fact that we do not have genetic or intrauterine data on the women included in this study, and so, these would be potential unmeasured confounders. In 1B, we assume that there is no causal relationship between PCOS and cardiometabolic health. Whilst observational studies have shown associations of PCOS with type 2 diabetes, and some also show associations with CHD, a recent Mendelian randomisation study has suggested no causal effect of PCOS on type 2 diabetes, coronary heart disease or stroke , making this scenario plausible. In this scenario, PCOS is not on a confounding path, and adjusting for it should not alter the results. Figure 1c and d suggest PCOS does causally influence cardiometabolic health. In 1C, PCOS is on a confounding path, and we would want to adjust for it. In 1D, it is on a mediating path, and we would not want to adjust for it. Given we cannot be certain which of these scenarios is correct in additional analyses we, repeat all of our main analyses with those with known PCOS removed. We chose to remove them rather than add to a multivariable model because there were only 24 cases. Weight and height [used to calculate the body mass index (BMI)] were measured in light clothing and unshod. Weight was measured to the nearest 0.1 kg using Tanita scales; height was measured to the nearest 0.1 cm using a Harpenden stadiometer. Smoking status was categorised as ever versus never (to be considered for state-funded assisted conception, women had to have not smoked for at least 3 months, and this was confirmed by a negative cotinine breath test). All other confounders, including an existing diagnosis of PCOS, were obtained by a questionnaire or from medical notes when the women were originally recruited. Full details of the assessment of these confounders are provided in the supplementary material.
All analyses were conducted using R version 3.4.2 (R Foundation for Statistical Computing, Vienna, Austria). An analysis plan was written in June 2019. Characteristics were summarised as n, total range, mean, standard deviation, median and 25th and 75th quantiles (IQR) as appropriate. Multivariable linear regression was used to examine the associations of functional ovarian reserve markers (treated as exposures) with serum metabolic profiles (treated as outcomes). Robust standard errors were estimated for all associations, as some metabolite concentrations had skewed distributions. The metabolic measures were scaled to standard deviation (SD) units (by subtracting the mean and dividing by the standard deviation of all women included in the analyses). This scaling allows easy comparison of multiple metabolic measures with different units or with large differences in their concentration distributions. AMH and AFC were also scaled to SD units in the same way as the metabolites. This was done for ease of comparison of results between the two biomarkers. Associations were adjusted for all a priori selected confounders (age, ethnicity, education, family history of cardiovascular disease, BMI, physical activity, alcohol, smoking, duration of infertility and whether the woman had primary or secondary infertility).
In addition to presenting our main results as the difference in mean metabolite in SD units per SD of AMH or AFC, we also present the full results (confounder adjusted) in the metabolite, AMH and AFC original units in the supplementary material. One woman was taking lipid-lowering medication, and removal of her from the analyses did not alter any of the findings from the main results (results available from the corresponding author on request). We repeated our main analyses only on those women with a known partner cause of infertility (N = 87 (22% of the cohort). We knew a priori that this sample size would lack statistical power for reliable results but wanted to compare the point estimates in this group to the results of the whole cohort to provide some indication as to whether our results might be driven by the cause of infertility in the women or potentially generalisable to women of reproductive age without infertility. We compared the magnitudes of the results in women who were going to undergo assisted conception because of partner infertility with those of the whole cohort using a scatterplot. As 99% of the women recruited had full data on ovarian reserve and at least one metabolite, with < 1% of these having missing covariable data (see Table 1), we did not need to undertake any additional analyses to explore potential biases due to missing data. To explore the departure from linearity where there was evidence of an association, AMH and AFC were split into quarters of their distribution and regression models were run with these quarters as a continuous score and as a 4 level categorical variable (with 3 indicators). A likelihood ratio test was used to compare these two models. Statistical support that the second (3 indicators) model was a better fit of the data would suggest a possible non-linear association.
Accounting for multiple testing
Due to the correlated nature of the metabolic biomarkers, over 95% of the variation in the 155 metabolic biomarkers was explained by 14 principal components. Therefore, multiple testing correction, accounting for 14 independent tests using the Bonferroni method, resulted in p < 0.0036 (0.05/14) being denoted as statistically significant.
Three hundred and ninety-eight women (99% of the 400 recruited) with available AMH and AFC levels and data on at least one metabolite were included in the study. The characteristics of the participants are shown in Table 1. Mean (SD) age of the women was 35.5 (4.4) years, and the majority (92%) were White European, with over 50% having a university degree, 26% being ever smokers, 68% exercising more than 3–4 times per week, 52% having a family history of cardiovascular disease and median alcohol consumption being 4 units per week. Mean (SD) BMI was 25 (3) kg/m2, and 95% were due to start their first treatment cycle, with 49% having an unknown cause of infertility and 24% a cause related to their partner. The median AMH was 16.1 pmol/l (IQR 8.8, 28.0 pmol/l) and median AFC of 12 (IQR 7, 16). AMH and AFC were positively correlated (Spearman’s correlation coefficient = 0.55, p < 0.001).
The unadjusted associations of AMH and AFC with selected confounders are shown in Additional file 1 Tables S2 and S3, respectively. Age was negatively associated with AMH (difference in mean − 0.09 SD per 1-year older age, 95% CI − 0.11, − 0.07) and AFC (− 0.06 SD per 1-year older age, 95% CI − 0.08, − 0.04). An ovulatory cause of infertility compared with an unexplained cause was associated with a higher AMH (difference in mean 0.93 SD, 95% CI 0.48, 1.37) and AFC (0.65, 95% CI 0.20, 1.11). Secondary infertility was strongly negatively associated with AFC (− 0.34SD (95% CI − 0.55, − 0.13), with a weaker negative association with AMH (− 0.11 (95%CI − 0.32, 0.11)). Family history of cardiometabolic diseases was negatively associated, and non-White European ethnicity positively associated with AMH and AFC. There was no strong statistical evidence of an association between BMI or other potential confounders and either measure of ovarian reserve.
The adjusted (for all potential confounders) associations between AMH and AFC and the respective metabolomics measures are SD units per SD of AMH/AFC and are presented in Figs. 2, 3 and 4. With the unadjusted associations presented Additional file 1 Figures S2 - S4. Confounder-adjusted results in the original units of AMH/AFC and the metabolites are shown in Additional file 2 Tables S4 and S5.
In confounder-adjusted analyses, AMH levels were positively correlated with measures of components of medium HDL and small HDL and with concentrations of omega-6 fatty acids and polyunsaturated fatty acids (PUFA). AMH was also positively associated with the amino acids isoleucine, leucine and tyrosine, and it was negatively associated with acetate concentrations. Although there were positive associations with lipoprotein subclasses, other cholesterol subtypes, glycerides and apolipoproteins, several of which had point estimates of ~ 0.1 SD or larger differences; these had wide confidence intervals that included the null and did not reach our multiple testing threshold for statistical significance.
Overall, the associations for AFC with the metabolic measures were in the same direction as those observed for AMH; however, they were weaker in magnitude and for many of them close to unity (Figs. 2, 3 and 4). The notable exceptions were the positive associations with the amino acids alanine, glutamine and glycine.
There was no strong evidence of departure from linearity for any of the observed associations (all likelihood ratio p values comparing a model with four categories as 3 indicator variables to the simpler model as a 4 level score ≥ 0.05). When analyses were repeated only in those undergoing assisted conception due to partner infertility (n = 87), most of the 155 associations of AMH were stronger than those in the main analyses including all women (Fig. 5a). This included the strengthening of the previously noted associations of AMH with the metabolomic markers. Positive associations were also observed for the various lipid concentrations within the medium and small very-low-density lipoproteins, intermediate-density lipoprotein and low-density lipoproteins, cholesterol, the mono-unsaturated fatty acids and alanine, with a negative association with glucose emerging; the negative association with acetate attenuated (Additional file 1 Figures S5-S7). In this subgroup, some associations of AFC also strengthened (Fig. 5b), with positive associations with various lipid concentrations across the size range of the very-low-density lipoproteins, triglycerides, apolipoprotein-B, mono-unsaturated fatty acids and the amino acids alanine and glutamine (Additional file 1 Figures S8-S10). When analyses were compared between the whole cohort and after exclusion of women with PCOS (n = 24), the associations were similar for AMH and AFC with correlations of 0.97 and 0.98, respectively (Additional file 1 Figure S11).
In this cross-sectional study of women attending a fertility clinic for evaluation, we identify several novel associations between the biomarkers of ovarian reserve, as measured by AMH and AFC, and circulating metabolites. Specifically, AMH showed positive associations with the amino acids isoleucine, leucine and tyrosine; medium HDL and small HDL; omega-6 fatty acids; and polyunsaturated fatty acids and a negative association with acetate. AFC had directionally consistent associations to those seen for AMH, but the magnitudes of association (measured on the same scale as SD difference in metabolite for a 1 SD increase in AMH or AFC) were weaker overall. Statistical evidence for positive associations with alanine, glutamine and glycine were observed. In our main analyses, including all women, we did not find strong statistical support for either biomarker being associated with an extensive lipid profile. This was still the case after the exclusion of women with PCOS. When we restricted analyses to the subgroup (22%) of women with male partner infertility, observed associations strengthened and we observed positive associations of AMH and AFC with a range of lipids, suggesting that the women infertility phenotypes may mask associations with lipid profiles.
Previous studies have suggested an association of lower ovarian reserve with vascular health [7, 33, 34]. Whether alteration in serum lipids contributes to this association is unclear. One small study (n = 50) reported a positive correlation between AMH and total cholesterol and LDL-C, and inverse correlation with HDL-C , whilst a larger study (N = 252) of women with PCOS found a weak positive correlation with HDL-C, which was attenuated when adjusted for BMI . In contrast, larger studies of Chinese women did not observe an association between AMH and any of total cholesterol, LDL-C, HDL-C or triglycerides in either women with PCOS (N = 304), infertile women (N = 1896)  or in a general cohort of women (N = 6763) . This is consistent with our own findings, with limited statistical evidence of an association between AMH and lipoprotein subclasses, lipoprotein particle size, cholesterol, LDL or HDL subtypes or triglycerides, phosphoglycerides, or apolipoproteins when the whole population was considered. Collectively, the evidence to date (including from our study) would support recent observational studies suggesting that the association of AMH with atherosclerosis may be independent of lipid levels [7, 33, 34]. As an alternative mechanism, it has been proposed that AMH has direct effects on cardiovascular tissue , making it more prone to injury and atherosclerosis . Though specific evidence for this is lacking. We did not demonstrate associations between AFC and lipids, which have to our knowledge not been previously explored, further supporting suggestions that alternative mechanisms to dyslipidaemia may underlie the relationships between ovarian reserve and cardiovascular disease. Though as acknowledged below, the greater variability in AFC compared with AMH may have attenuated some of these results towards the null.
We report a positive association between AMH and higher circulating levels of PUFA and omega-6 fatty acids. However, a systematic review and meta-analysis found no evidence of either observational associations between dietary intake or measured circulating concentrations of omega-6 on cardiovascular diseases or any effect of dietary supplementation with omega-6 in randomised controlled trials , suggesting that even if the association of AMH with omega-6 is causal, this is unlikely to be a mechanism for preventing cardiovascular diseases.
Branch chain amino acids (BCAA), such as isoleucine, leucine and valine, have been found to be positively associated with a number of cardiometabolic risk factors, including adiposity, fasting glucose, insulin resistance, blood pressure, dyslipidaemia and indicators of coronary artery disease, in cross-sectional studies . There are also positive associations of BCAA and aromatic amino acids with incident cardiovascular events in several large prospective studies [14, 42,43,44]. However, whether these are causally related is at present unknown. Our observed AMH-BCAA association is of interest, but the replication of our findings and confirmation of biological plausibility of causality would be required to assess whether the association of AMH with CVD may be mediated to some extent by BCAA.
The same NMR platform has previously been used to identify 14 metabolites associated with all-cause mortality in a meta-analysis of 12 cohorts and 44,168 participants , with subgroup analyses of 7603 participants identifying seven metabolites (XXL-VLDL-L, PUFA, lactate, histidine, leucine, phenylalanine and albumin) inversely associated with cardiovascular mortality and three (glucose, lactate and glycoprotein acetyls) positively associated with cardiovascular mortality . Overall, our observed direction of associations was consistent with markers of low ovarian reserve relating to these metabolites that have previously been shown to robustly associate with cardiovascular disease mortality. Associations of AMH with PUFA and leucine met our adjusted statistical threshold. These findings, if replicated in larger cohorts, would support the overall concept that a diminished ovarian reserve may be associated with an unfavourable circulating cardiometabolic risk profile .
Our study has several strengths. To our knowledge, we are not aware of any study with a similar or larger sample size with detailed phenotypic and metabolite measurements. All women attended during the early follicular phase for the measurement of AMH and AFC. We included all women across the range of ovarian reserve including some women with extremely low and high AMH and AFC. We wrote, and have used, a prespecified analysis plan which incorporated adjustment for a wide range of a priori specified potential confounders. Subgroup analyses in those women who were awaiting assisted conception because of male partner infertility were also assessed, as these women may reflect a general population of women of reproductive age.
We do however acknowledge several limitations. Our analyses are cross-sectional and therefore could be explained by variation in metabolism (e.g. of amino acids) influencing ovarian reserve rather than the other way around, as we have assumed. Furthermore, we cannot assume that the small number of associations identified is causal. Residual confounding may have resulted from crude questionnaire measurements of physical activity, alcohol intake and family history of cardiovascular disease and the lack of any data on dietary intake. However, we adjusted for the measures we had of alcohol, physical activity and family history of cardiovascular disease, as well as education and BMI, which influence diet and physical activity, or are influenced by it, and thus may have captured key confounding paths. Of note, our adjustment for confounders is more extensive than previous studies of the association of AMH with lipids and other cardiovascular risk factors. Whilst this is one of the larger studies to explore these associations, our results were imprecisely estimated, with wide confidence intervals, and we are not aware of any independent study that has measures of ovarian reserve and multiple metabolites (or even the amino acids and fatty acids that we observed associations with) in which to attempt to replicate our findings.
Women were awaiting IVF and related to that were confirmed (through cotinine breath test) non-smokers and were of a relatively restricted BMI range. The women were also largely White European and educated to degree level. This homogeneous relatively healthy population may have resulted in some selection bias and may mean that our results do not generalise to a general population of women of reproductive age or other infertile populations. AFC was measured by several operators, and the known intra- and inter-observer variability of AFC  and application of a threshold for counting per ovary, with the highest value for AFC of 40 (range 0–40), as compared to AMH of 171 (range 1–171 pmol/l), may explain why we did not see the same strength of associations of AFC with the various metabolites compared to what we observed with AMH. AMH was measured in a single laboratory on an automated analyser with a wide analytical range and low coefficient of variation.
Analyses were undertaken on non-fasting samples. This was necessary to align with clinical processes for a population who are undergoing assisted conception, where caloric restraint may be detrimental. In collaborations of several studies using this same NMR analysis platform results have not differed notably between studies in which the analyses were undertaken in participants who had been advised to fast and those who had not, including in analyses exploring the associations of these metabolites with cardiovascular diseases .
The NMR platform used in these analyses covers considerably more of the lipidome than conventional clinical chemistry measures (total cholesterol, LDL-C, HDL-C and triglycerides) that have previously been explored in relation to ovarian reserve, and in addition includes fatty acids, amino acids, glycolysis metabolites, ketone bodies and an inflammatory marker. We acknowledge it misses a high proportion of the currently quantifiable metabolites in human serum/plasma, including markers of energy balance, microbiota metabolism, vitamins, co-factors and xenobiotics, that may be influenced by ovarian reserve. High-throughput analyses of a wider range of metabolites measured by mass spectrometry are possible in epidemiological studies but considerably more expensive than the NMR platform used here, resulting in their use being frequently restricted to subsamples of cohorts .
This study provides novel insight into the association of the ovarian biomarkers, AMH and AFC, with metabolic profiles. Taken together with other recent studies, our results suggest that dyslipidaemia may have a limited role to play in the relationship between ovarian reserve and cardiovascular diseases. The novel associations we find with some fatty acids and amino acids may have a role in mediating any effect of ovarian reserve on cardiometabolic diseases, but these require replication in large prospective studies and we cannot assume that the small number of associations identified are causal.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Antral follicle count
Branched-chain amino acid
Body mass index
Directed acyclic graph
High-density lipoprotein cholesterol
Homeostatic Model Assessment of Insulin Resistance
In vitro fertilisation
Low-density lipoprotein cholesterol
Nuclear magnetic resonance
Polycystic ovarian syndrome
Polyunsaturated fatty acid
Very low-density lipoprotein
te Velde ER, Pearson PL. The variability of female reproductive ageing. Hum Reprod Update. 2002;8(2):141–54.
Broekmans FJ, Faddy MJ, Scheffer G, te Velde ER. Antral follicle counts are related to age at natural fertility loss and age at menopause. Menopause. 2004;11(6 Part 1 of 2):607–14.
Dolleman M, Faddy MJ, van Disseldorp J, van der Schouw YT, Messow CM, Leader B, Peeters PH, McConnachie A, Nelson SM, Broekmans FJ. The relationship between anti-Mullerian hormone in women receiving fertility assessments and age at menopause in subfertile women: evidence from large population studies. J Clin Endocrinol Metab. 2013;98(5):1946–53.
Bleil ME, Gregorich SE, McConnell D, Rosen MP, Cedars MI. Does accelerated reproductive aging underlie premenopausal risk for cardiovascular disease? Menopause. 2013;20(11):1139–46.
de Kat AC, Verschuren WM, Eijkemans MJ, van der Schouw YT, Broekmans FJ. The association of low ovarian reserve with cardiovascular disease risk: a cross-sectional population-based study. Hum Reprod. 2016;31(8):1866–74.
Tehrani FR, Erfani H, Cheraghi L, Tohidi M, Azizi F. Lipid profiles and ovarian reserve status: a longitudinal study. Hum Reprod. 2014;29(11):2522–9.
de Kat AC, Verschuren WM, Eijkemans MJC, Broekmans FJM, van der Schouw YT. Anti-Mullerian hormone trajectories are associated with cardiovascular disease in women. Circulation. 2017;135(6):556–65.
Dam V, van der Schouw YT, Onland-Moret NC, Groenwold RHH, Peters SAE, Burgess S, Wood AM, Chirlaque M-D, Moons KGM, Oliver-Williams C, et al. Association of menopausal characteristics and risk of coronary heart disease: a pan-European case–cohort analysis. Int J Epidemiol. 2019;48(4):1275–85.
Yarde F, Group tOS, Spiering W, Group tOS, Franx A, Group tOS, Visseren FLJ, Group tOS, Eijkemans MJC, Group tOS et al: Association between vascular health and ovarian ageing in type 1 diabetes mellitus. Hum Reprod 2016, 31(6):1354–1362.
Anderson EL, Fraser A, McNally W, Sattar N, Lashen H, Fleming R, Nelson SM, Lawlor DA. Anti-mullerian hormone is not associated with cardiometabolic risk factors in adolescent females. PLoS One. 2013;8(5):e64510.
Ramezani Tehrani F, Montazeri SA, Khalili D, Cheraghi L, Broekmans FJ, Momenan AA, de Kat AC, Azizi F. Age-specific anti-Müllerian hormone and electrocardiographic silent coronary artery disease. Climacteric. 2016;19(4):344–8.
Park HT, Cho GJ, Ahn KH, Shin JH, Kim YT, Hur JY, Kim SH, Lee KW, Kim T. Association of insulin resistance with anti-Mullerian hormone levels in women without polycystic ovary syndrome (PCOS). Clin Endocrinol. 2010;72(1):26–31.
Akbaraly T, Wurtz P, Singh-Manoux A, Shipley MJ, Haapakoski R, Lehto M, Desrumaux C, Kahonen M, Lehtimaki T, Mikkila V, et al. Association of circulating metabolites with healthy diet and risk of cardiovascular disease: analysis of two cohort studies. Sci Rep. 2018;8(1):8620.
Wurtz P, Havulinna AS, Soininen P, Tynkkynen T, Prieto-Merino D, Tillin T, Ghorbani A, Artati A, Wang Q, Tiainen M, et al. Metabolite profiling and cardiovascular event risk: a prospective study of 3 population-based cohorts. Circulation. 2015;131(9):774–85.
Deelen J, Kettunen J, Fischer K, van der Spek A, Trompet S, Kastenmuller G, Boyd A, Zierer J, van den Akker EB, Ala-Korpela M, et al. A metabolic profile of all-cause mortality risk identified in an observational study of 44,168 individuals. Nat Commun. 2019;10(1):3346.
Wurtz P, Kangas AJ, Soininen P, Lawlor DA, Davey Smith G, Ala-Korpela M. Quantitative serum nuclear magnetic resonance metabolomics in large-scale epidemiology: a primer on -Omic Technologies. Am J Epidemiol. 2017;186(9):1084–96.
Wurtz P, Wang Q, Kangas AJ, Richmond RC, Skarp J, Tiainen M, Tynkkynen T, Soininen P, Havulinna AS, Kaakinen M, et al. Metabolic signatures of adiposity in young adults: Mendelian randomization analysis and effects of weight change. PLoS Med. 2014;11(12):e1001765.
Wurtz P, Tiainen M, Makinen VP, Kangas AJ, Soininen P, Saltevo J, Keinanen-Kiukaanniemi S, Mantyselka P, Lehtimaki T, Laakso M, et al. Circulating metabolite predictors of glycemia in middle-aged men and women. Diabetes Care. 2012;35(8):1749–56.
Wang Q, Wurtz P, Auro K, Makinen VP, Kangas AJ, Soininen P, Tiainen M, Tynkkynen T, Jokelainen J, Santalahti K, et al. Metabolic profiling of pregnancy: cross-sectional and longitudinal evidence. BMC Med. 2016;14(1):205.
Wang Q, Ferreira DLS, Nelson SM, Sattar N, Ala-Korpela M, Lawlor DA. Metabolic characterization of menopause: cross-sectional and longitudinal evidence. BMC Med. 2018;16(1):17.
Teede HJ, Misso ML, Costello MF, Dokras A, Laven J, Moran L, Piltonen T, Norman RJ. Recommendations from the international evidence-based guideline for the assessment and management of polycystic ovary syndrome. Hum Reprod. 2018;33(9):1602–18.
Broekmans FJ, de Ziegler D, Howles CM, Gougeon A, Trew G, Olivennes F. The antral follicle count: practical recommendations for better standardization. Fertil Steril. 2010;94(3):1044–51.
Soininen P, Kangas AJ, Wurtz P, Tukiainen T, Tynkkynen T, Laatikainen R, Jarvelin MR, Kahonen M, Lehtimaki T, Viikari J, et al. High-throughput serum NMR metabonomics for cost-effective holistic studies on systemic metabolism. Analyst. 2009;134(9):1781–5.
Fischer K, Kettunen J, Wurtz P, Haller T, Havulinna AS, Kangas AJ, Soininen P, Esko T, Tammesoo ML, Magi R, et al. Biomarker profiling by nuclear magnetic resonance spectroscopy for the prediction of all-cause mortality: an observational study of 17,345 persons. PLoS Med. 2014;11(2):e1001606.
White SL, Pasupathy D, Sattar N, Nelson SM, Lawlor DA, Briley AL, Seed PT, Welsh P, Poston L. Metabolic profiling of gestational diabetes in obese women during pregnancy. Diabetologia. 2017;60(10):1903–12.
Mills HL, Patel N, White SL, Pasupathy D, Briley AL, Santos Ferreira DL, Seed PT, Nelson SM, Sattar N, Tilling K, et al. The effect of a lifestyle intervention in obese pregnant women on gestational metabolic profiles: findings from the UK Pregnancies Better Eating and Activity Trial (UPBEAT) randomised controlled trial. BMC Med. 2019;17(1):15.
Pearce N, Lawlor DA. Causal inference-so much more than statistics. Int J Epidemiol. 2016;45(6):1895–903.
Williamson EJ, Aitken Z, Lawrie J, Dharmage SC, Burgess JA, Forbes AB. Introduction to causal diagrams for confounder selection. Respirology (Carlton, Vic). 2014;19(3):303–11.
Zegers-Hochschild F, Adamson GD, Dyer S, Racowsky C, de Mouzon J, Sokol R, Rienzi L, Sunde A, Schmidt L, Cooke ID, et al. The international glossary on infertility and fertility care, 2017. Human Reproduction (Oxford, England). 2017;32(9):1786–801.
Matthews KA, Crawford SL, Chae CU, Everson-Rose SA, Sowers MF, Sternfeld B, Sutton-Tyrrell K. Are changes in cardiovascular disease risk factors in midlife women due to chronological aging or to the menopausal transition? J Am Coll Cardiol. 2009;54(25):2366–73.
Zhu T, Goodarzi MO. SUN-016 Mendelian randomization reveals that polycystic ovary syndrome is not a causal risk factor for type 2 diabetes, coronary heart disease, or stroke. J Endocr Soc. 2020;4(Suppl 1):SUN-016.
Pearce N, Lawlor DA. Causal inference—so much more than statistics. Int J Epidemiol. 2017;45(6):1895–903.
Appt SE, Chen H, Clarkson TB, Kaplan JR. Premenopausal antimullerian hormone concentration is associated with subsequent atherosclerosis. Menopause. 2012;19(12):1353–9.
Looby SE, Fitch KV, Srinivasa S, Lo J, Rafferty D, Martin A, Currier JC, Grinspoon S, Zanni MV. Reduced ovarian reserve relates to monocyte activation and subclinical coronary atherosclerotic plaque in women with HIV. AIDS (London, England). 2016;30(3):383–93.
Skałba P, Cygal A, Madej P, Dąbkowska-Huć A, Sikora J, Martirosian G, Romanik M, Olszanecka-Glinianowicz M. Is the plasma anti-Müllerian hormone (AMH) level associated with body weight and metabolic, and hormonal disturbances in women with and without polycystic ovary syndrome? Eur J Obstetrics Gynecol Reproductive Biol. 2011;158(2):254–9.
Feldman RA, O’Neill K, Butts SF, Dokras A. Antimullerian hormone levels and cardiometabolic risk in young women with polycystic ovary syndrome. Fertil Steril. 2017;107(1):276–81.
Cui Y, Shi Y, Cui L, Han T, Gao X, Chen Z-J. Age-specific serum antimüllerian hormone levels in women with and without polycystic ovary syndrome. Fertility Sterility. 2014;102(1):230–6 e232.
Cui L, Qin Y, Gao X, Lu J, Geng L, Ding L, Qu Z, Zhang X, Chen Z-J. Antimüllerian hormone: correlation with age and androgenic and metabolic factors in women from birth to postmenopause. Fertility Sterility. 2016;105(2):481–5 e481.
McLennan IS, Pankhurst MW. Anti-Mullerian hormone is a gonadal cytokine with two circulating forms and cryptic actions. J Endocrinol. 2015;226(3):R45–57.
Chowdhury R, Warnakula S, Kunutsor S, Crowe F, Ward HA, Johnson L, Franco OH, Butterworth AS, Forouhi NG, Thompson SG, et al. Association of dietary, circulating, and supplement fatty acids with coronary risk: a systematic review and meta-analysis. Ann Intern Med. 2014;160(6):398–406.
Tobias DK, Mora S, Verma S, Lawler PR. Altered branched chain amino acid metabolism: toward a unifying cardiometabolic hypothesis. Curr Opin Cardiol. 2018;33(5):558–64.
Shah SH, Bain JR, Muehlbauer MJ, Stevens RD, Crosslin DR, Haynes C, Dungan J, Newby LK, Hauser ER, Ginsburg GS, et al. Association of a peripheral blood metabolic profile with coronary artery disease and risk of subsequent cardiovascular events. Circ Cardiovasc Genet. 2010;3(2):207–14.
Magnusson M, Lewis GD, Ericson U, Orho-Melander M, Hedblad B, Engstrom G, Ostling G, Clish C, Wang TJ, Gerszten RE, et al. A diabetes-predictive amino acid score and future cardiovascular disease. Eur Heart J. 2013;34(26):1982–9.
Ruiz-Canela M, Toledo E, Clish CB, Hruby A, Liang L, Salas-Salvado J, Razquin C, Corella D, Estruch R, Ros E, et al. Plasma branched-chain amino acids and incident cardiovascular disease in the PREDIMED Trial. Clin Chem. 2016;62(4):582–92.
Iliodromiti S, Anderson RA, Nelson SM: Technical and performance characteristics of anti-Mullerian hormone and antral follicle count as biomarkers of ovarian response. Hum Reprod Update 2015;21(6):698–710. https://doi.org/10.1093/humupd/dmu062. Epub 2014 Dec 8.
Sovio U, Goulding N, McBride N, Cook E, Gaccioli F, Charnock-Jones DS, Lawlor DA, Smith GCS. A maternal serum metabolite ratio predicts fetal growth restriction at term. Nat Med. 2020;26(3):348–53.
We thank the women who participated in the study.
This work was supported by a Kuwait Government Fellowship to KA, the National Institute for Health Research Biomedical Centre at the University Hospitals Bristol NHS Foundation Trust and the University of Bristol (AT, DAL, and SMN), a European Research Council grant (DevelopObese; 669545 to DAL), European Union’s Horizon 2020 research and innovation programme (LifeCycle; 733206 to DA), a British Heart Foundation Grant (AA/18/7/34219 to DAL) and a National Institute for Health Research Senior Investigator award (NF-0616-10102 to DAL). AT, NG and DAL work in a Unit that receives support from the University of Bristol and the MRC (MC_UU_00011/6). The views expressed in this publication are those of the authors and not necessarily those of the NHS, the National Institute for Health Research or the Department of Health and Social Care, or any other funders mentioned here.
Ethics approval and consent to participate
NHS Health Research Authority provided ethical approval for the study. REC reference 16/WM/0308. IRAS project ID:202216. All participants provided signed consent.
Consent for publication
No funding bodies had any role in the study design, data collection and analysis, decision to publish or preparation of the manuscript. SMN has participated in Advisory Boards and received speaker’s fees from Beckman Coulter, Ferring, Finox, Merck, MSD and Roche Diagnostics. DAL has received grant funding for other studies, not related to this one, from government, charity and industry funders, including Roche Diagnostics and Medtronic. DAL is a member of the Editorial Board for BMC Medicine.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
List of NMR metabolomic measures assessed in this study, with their units of quantification. Table S2. Associations of possible confounders with AMH. Table S3. Associations of possible confounders with AFC. Fig. S1. Summary of NMR platform used in this paper. Fig. S2. Associations of lipoprotein classes with AMH/AFC (unadjusted). Fig. S3. Associations of lipoprotein classes and fatty acids with AMH/AFC (unadjusted). Fig. S4. Associations of metabolic traits with AMH/AFC (unadjusted). Fig. S5. Association of lipoprotein classes with AMH in women with male factor infertility or no male partner. Fig. S6. Association of lipoprotein classes and fatty acids with AMH in women with male factor infertility or no male partner. Fig. S7. Association of metabolic traits with AMH in women with male factor infertility or no male partner. Fig. S8. Association of lipoprotein classes with AFC in women with male factor infertility or no male partner. Fig. S9. Association of lipoprotein classes and fatty acids with AFC in women with male factor infertility or no male partner. Fig. S10 Association of metabolic traits with AFC in women with male factor infertility or no male partner. Fig. S11. Comparison of associations for AMH and AFC for whole cohort and with PCOS cases (n = 24) removed.
Adjusted associations of metabolites with AMH in original units. Table S5. Adjusted associations of metabolites with AFC in original units.
About this article
Cite this article
Al Rashid, K., Taylor, A., Lumsden, M.A. et al. Association of the functional ovarian reserve with serum metabolomic profiling by nuclear magnetic resonance spectroscopy: a cross-sectional study of ~ 400 women. BMC Med 18, 247 (2020). https://doi.org/10.1186/s12916-020-01700-z