Abstract
Objective:
Current research results on the association between parity and cardiovascular disease (CVD) risk are inconclusive. The purpose of this study was to examine the relationship between parity and risk of CVD in women over 45 years of age.
Materials and Methods:
Data were from the National Health and Nutrition Examination Survey for 2007–2018. A total of 8,882 women ≥45 years of age were included. The exposure of parity referred to the number of live births, and the outcome variable was the occurrence of CVD events, including stroke, heart attack, coronary heart disease (CHD), heart failure, and angina. We performed logistic regression to calculate the unadjusted and adjusted odds ratios (ORs; 95% confidence intervals [CIs]) controlling for confounding factors.
Results:
Among 8,882 women, the mean age was 62.4 ± 10.8 years, with a range of 45 to 80 years. The weighted prevalence of CVD, stroke, heart attack, CHD, heart failure, and angina in parous women were significantly higher than those of nulliparous women (p < 0.05). After adjusting for demographic factors, CVD risk factors, reproductive factors, the ORs of CVD with parity 1–2, 3–4, and 5+ were 1.85 (95% CI: 1.29–2.64), 1.70 (95% CI: 1.15–2.50), and 1.92 (95% CI: 1.28–2.88), respectively. The odds of stroke and heart attack were also significantly positively related to increasing parity compared with nulliparity. However, compared with nulliparity, parity of 3 was associated with a slightly lower risk of CVD.
Conclusions:
Our findings indicated that parity was significantly positively associated with CVD, giving birth to three children associated with a slightly lower risk compared with nulliparity. Further cohort studies are warranted to confirm the findings.
Introduction
Most countries around the world are facing the harsh reality of significantly dropping fertility rates and declining willingness to bear children. 1 –3 Among women 15–44 years of age in the United States, the birth rate has dropped from 70‰ in 2007 to 55.8‰ in 2020, a drop of about 20% in 13 years. 4 In Thailand, the average number of births per woman during her reproductive lifespan has decreased significantly every year, and is expected to continue to decline to 1.3 by 2040. 5 The birth rate in China has dropped astonishingly from 37‰ in 1950 to 8.52‰ in 2020 with a drop of 77% among women 15–49 years of age, and the average number of births per woman during her reproductive lifespan declined from 6.0 in the 1960s to 1.3 in 2020. 6 Faced with such changes in the global fertility rate, the relationship between giving birth and women's future health is an issue worthy of attention. Whether giving birth will influence the development of chronic conditions in later life, and the number of births most associated with healthy outcomes are controversial.
Cardiovascular disease (CVD) is the leading cause of morbidity and mortality worldwide. 7 For women, in addition to traditional CVD risk factors, some reproductive factors may also affect the development of maternal CVD in later life. 8 Pregnancy is related to cardiometabolic changes such as weight gain and dyslipidemia. 9,10 Physiological changes related to pregnancy may make women susceptible to CVD through inflammation and endothelial function. 11,12 So, a broader agenda has been warranted to integrate reproductive health into women's CVD research.
The association of nulliparity and parity with CVD has been controversial. 13 –15 Several previous studies found that nulliparous women had the lowest risk of CVD, 16 –18 while others indicated that nulliparity increased women's CVD risk compared with women with parity of 1 or 2. 19 –21 Some population-based studies did not find a significant association between parity and CVD in women. 13,14,22 Therefore, the purpose of the study was to use National Health and Nutrition Examination Survey (NHANES) data to examine the association between parity and CVD in women over 45 years of age in a large and representative population in the United States. Our hypothesis was that parous women 45 years of age or older would be likely to have an increased risk of CVD.
Materials and Methods
Data source and study population
NHANES is a nationally representative cross-sectional set of health surveys conducted by the National Center for Health Statistics (NCHS) of the Centers for Disease Control and Prevention (CDC). 23 The CDC's NCHS Ethics Review Board approved NHANES data collection from the participants. It utilized a stratified, multistage probability sampling design for selection of noninstitutionalized American civilians. The survey includes assessments of diet, general health, physical examinations, and lifestyle factors. The datasets are from public-use files. We analyzed data from six 2-year NHANES datasets for the period 2007–2018. The research population was women 45 years of age or older at the time of interview. There were 10,379 women over 45 years of age in the dataset. We excluded 121 women with the response of “don't know” or “refused” or missing values for the conditions of stroke, heart attack, heart failure, coronary heart disease (CHD), and angina.
The datasets have the variable of age at these CVD events and age at the last live birth. We also excluded those participants who had these conditions of CVDs before their last live births and nulliparous women who had these CVD conditions before the age of 29 years (the average age of last live birth, n = 69). Then, we imputed missing values of parity with the sum of the number of vaginal and cesarean deliveries, and deleted cases of unknown parity, excluding the participants whose history of pregnancy was the response, “don't know” or “refused” or missing (n = 1307). Thus, a total of 8,882 participants were included in the present analyses.
Dependent and independent variables
The primary outcomes were stroke, heart attack, CHD, heart failure, and angina pectoris based on the self-reported responses to the questions “Has a health care provider or physician ever told you that you had a stroke/heart attack/CHD/heart failure/angina?”. If any of the above diseases were reported, we considered the participant to have CVD. The independent variable of parity referred to the number of live births. Parity was determined by responses to several questions, including “Have you ever been pregnant?” “Are you pregnant now?” and “How many live births have you ever had?” We categorized parity as follows: if participants were never pregnant or the number of live births was zero, we set the value of parity as “nulliparity”; if the participant had at least one live birth, they were parous women. For parous women, the number of live births was categorized as 1–2, 3–4, and 5+ of parity, respectively. 24 We also consider parity as continuous independent variable including 0, 1, 2, 3, 4, and 5+ to explore the association between parity and CVD.
Potential confounding variables
The potential covariates included age at the interview, race, marital status, education level, family monthly poverty level index (FMPLI), body mass index (BMI), smoking, alcohol drinking, physical activity, healthy eating, family history of heart attack, hypertension, diabetes mellitus, age at menarche, age at menopause, hysterectomy, bilateral oophorectomy, ever taking birth control pills, and taking female hormones. 8,10 –12 Age at the interview, BMI, physical activity, age at menarche, and age at menopause were continuous variables. Height and weight were measured in the mobile examination centers by trained health technicians. Body weight was measured using a digital weight scale, and standing height was measured using a stadiometer with a fixed vertical backplate and adjustable headpiece. BMI was weight (kilogram) divided by height (meters) squared, rounded to one decimal place. We quantified physical activity by metabolic equivalent scores (METs).
The METs was the sum of the minutes of vigorous work-related/leisure-time activity and moderate vigorous work-related/leisure-time activity each week. The weight coefficients of vigorous activity and moderate activity were 8 and 4, respectively. Then, we calculated the average METs/hour per week as the value of physical activity variable. 25 The classification of race/ethnicity was non-Hispanic White, non-Hispanic Black, Hispanic and others, including Asian and Multiracial. Education level included high school or below and college/university or above. We divided marital status into married/living with partner, widowed/divorced/separated, and never married. Smoking was based on the two questions: “Have you smoked at least 100 cigarettes in your life?” and “Do you currently smoke?”.
We grouped smoking status into three categories: now smokers referred to participants smoking at least 100 cigarettes in life, and still smoking either every day or some days; former smokers referred to those smoking at least 100 cigarettes in life but did not smoke currently; never smoking referred to smoking less than 100 cigarettes in life. Alcohol drinking was based on the responses to the question “Have you had at least 12 alcohol drinks per year?” We defined the variable of healthy diet based on the question “How healthy is the diet,” including five different responses “excellent,” “very good,” “good,” “fair,” and “poor.” Family history of heart attack was determined from the question “Have your close relative, including father, mother, sisters, or brothers, ever had heart attack?”
The systolic and diastolic blood pressures (BP) of participants were measured by physicians up to four times after sitting and resting for 5 minutes. We calculated the average of systolic and diastolic BP of the second, third, and fourth times, discarding the BP value of the first time. We defined hypertension based on the questions, “Has a health care provider or physician ever told you that you had high blood pressure?” and “Are you taking prescription for hypertension,” and the average systolic BP was ≥140 mmHg or the average diastolic BP was ≥90mmHg. We defined diabetes mellitus by the questions, “Has a health care provider or physician ever told you that you had diabetes?” and “Are you taking antidiabetic medication/insulin?” and the fasting plasma glucose ≥126 mg/dL. 26 A history of high cholesterol, hysterectomy, bilateral oophorectomy, taking female hormones, taking birth control pills, all came from self-reports by the survey participants.
Statistical analyses
We followed NHANES analysis guidelines for merging data from 2007 to 2018 and applying sample weights. 27 We performed univariate analysis for descriptive statistics, producing weighted proportions (%). We used the Rao–Scott Chi-square test to compare weighted proportions among different parity levels. We used t-tests to compare the mean values of continuous variables between two groups, applying analysis of variance for comparing the mean values in multiple groups. We conducted multivariable analyses using logistic regression accounting for survey weights to obtain unadjusted and adjusted odds ratio (OR) values and 95% confidence intervals (CIs) treating CVD as the dependent variable and parity as the independent variable. We tried all potential covariates, and then compared the models with goodness-of-fit tests and model selection criteria (Akaike Information Criterion [AIC] and log-likelihood).
When performing the goodness-of-fit test, we listed four models. Model 1 was unadjusted, and model 2 was the full model, including all the covariates of age, race, education, marital status, FMPLI, smoking, alcohol drinking, BMI, physical activity, healthy diet, high blood pressure, cholesterol, diabetes mellitus, age at menarche, age at menopause, hysterectomy, oophorectomy, taking birth control pills, and taking female hormones. Model 3 consisted of significant variables in model 2 and added interaction term parity*race. Model 4 was the final model, omitting nonsignificant covariates in model 3. To obtain the optimal model, we took advantage of AIC, log-likelihood, and R-squared (R 2 ) criteria. Model 4 indicated a better fit model with the largest log-likelihood and R 2 and the smallest AIC. All statistical analyses were performed using SAS version 9.4 (SAS Institute Inc., Cary, NC, USA). A two-tailed p-value of <0.05 was considered having statistical significance.
Results
Basic characteristics of the study population
The study population consisted of 8,882 women over the age of 45 years, including 1,035 nulliparous women, 3,667 women with parity of 1–2, 3,017 women with a parity of 3–4, and 1,163 women with a parity of 5 or above. The mean age was 62.4 ± 10.8 years, with a range from 45 to 80 years. The weighted prevalence of nulliparity in women over 45 years was 13.0%, which could be representative of ∼7,501,386 females at the same age in the United States.
Distributions of race, marital status, education, and FMPLI differed significantly by parity level with nulliparous women and women in parity level 1–2 being more likely to have never been married, have a college education or higher and have a FMPLI >1.85 than women in other parity levels. Proportions of smoking, alcohol drinking, hypertension, diabetes mellitus, ever taking birth controls, taking female hormones, and having had hysterectomy differed significantly among the different parity groups. The mean of BMI, physical activity, age at menarche, and age at menopause were also statistically different in the groups. Proportions with high cholesterol or having had a bilateral oophorectomy or family history of heart attacks did not differ significantly among parity groups (Table 1).
Characteristics of Women 45 Years of Age and Over by Parity in National Health and Nutrition Examination Survey 2007–2018
The p value came from the analysis of variance or Rao–Scott Chi-square test among different subgroups. Physical activity indicated MET score/hour per week.
Means ± SDs.
p < 0.0001.
p < 0.05.
Alcohol drink, had at least 12 alcohol drinks per year; CHD, coronary heart disease; CVD, cardiovascular disease; FMPLI, family monthly poverty level index; NH, non-Hispanic; SD, standard deviation.
Comparison of the CVD prevalence between nulliparous and parous women
The unweighted and weighted prevalence of CVD in parous women were all significantly higher than those of nulliparous women. Specifically, the weighted prevalence of CVD, stroke, heart attack, CHD, heart failure, and angina were 12.1%, 4.9%, 3.9%, 4.0%, 3.6%, and 3.0%, respectively, in women with parity 1+, while the corresponding prevalence were 6.5%, 2.3%, 1. 6%, 1.8%, 2.0%, and 1.5% in nulliparous women (Table 2).
The Cardiovascular Disease Prevalence by Parity for Women ≥45 Years of Age in the United States, National Health and Nutrition Examination Survey 2007–2018
The Rao–Scott Chi-square p value was performed for comparing the weighted prevalence between parity groups.
Association between parity and CVD
The final model 4 included the covariates of age, race, FMPLI, smoking, BMI, physical activity, family history of heart attack, healthy diet, high BP, high cholesterol, diabetes mellitus, and age at menopause, omitting the interaction of parity and race that was not significant in CVD. The adjusted association of CVD with parity 1–2, 3–4, and 5+ compared with nulliparity were 1.85 (95% CI: 1.29–2.64), 1.70 (95% CI: 1.15–2.50), and 1.92 (95% CI: 1.28–2.88), respectively. Model 4 indicated that stroke and heart attack were also significantly related to parity compared with nulliparity after controlling for all the confounders. However, no significant relation was observed between parity and CHD, heart failure, and angina after adjusting for covariates (Table 3).
The Unadjusted and Adjusted Odds Ratios of Cardiovascular Diseases by Parity, National Health and Nutrition Examination Survey 2007–2018
Model 1 was unadjusted. Model 2 adjusted for age, race, education, marital status, FMPLI, smoke, alcohol drinking, BMI, physical activity, healthy diet, high blood pressure, cholesterol, diabetes mellitus, age at menarche, age at menopause, hysterectomy, oophorectomy, taking birth control pills and female hormones. Model 3 consisted of significant variables in model 2 and added the interaction term of parity*race. Model 4 was the final model, omitting nonsignificant variables in model 3.
BMI, body mass index; CI, confidence interval; OR, odds ratio.
Figure 1 presents the unadjusted and adjusted associations of CVD from model 4 with individual parity levels instead of grouping parity 1–2, 3–4, and 5+ presented in Table 3. Compared with nulliparity, the unadjusted ORs for CVD were similar for parity of 1–3, but were higher for parity of 4 and 5+. After adjustment for all the covariates in model 4, adjusted ORs were significant for all levels of parity compared with nulliparous women: 1.86 (95% CI: 1.21–2.86), 1.85 (95% CI: 1.27–2.69), 1.56 (95% CI: 1.05–2.32), 1.97 (95% CI: 1.28–3.03), and 1.94 (95% CI: 1.29–2.90).

The unadjusted and adjusted odds ratios between parity and CVD in women ≥45 years of age in the United States, NHANES 2007–2018. The adjusted confounders were model 4 covariates. CVD, cardiovascular disease; NHANES, National Health and Nutrition Examination Survey.
Discussion
We found that among women over 45 years of age in the United States, the weighted prevalence of CVD, stroke, heart attack, CHD, heart failure, and angina were lower in nulliparous women than those of parous women (Table 2). Women with parity 1–2, 3–4, and 5+ were more likely to suffer from stroke, heart attack, and CVD than nulliparous women in later life (Table 3). Compared with nulliparous women, the adjusted OR between CVD and parity 3 was slightly lower than other parous women after controlling for the potential confounders (Fig. 1).
Research on the association of CVD with nulliparity has yielded conflicting findings. 11 The results of Healthy Aging and Body Composition Study in Pittsburgh were consistent with our findings, reporting nulliparous older women with the average of 74 years had a lower CVD prevalence than parous women (18.0% vs. 30.2%). 16 In contrast, Parikh et al. reported that nulliparity conferred a moderate risk of CVD for women. 21 Hannaford et al. demonstrated that nulliparous women (mean age 56 years) were at higher risk of hypertension or stroke than parous women. 20
Our findings agree with those of multiple previous studies. Catov et al. reported that compared with nulliparity, the CVD prevalence of women with parity 1+ was 1.95-fold higher (95% CI: 1.03–3.7), parity 5+ having 2.27-fold (95% CI: 1.00 5.15) after adjustment for confounders. 16 One study integrating the data of the Framingham Heart Study and NHANES National Epidemic Follow-up Study showed women with six or more pregnancies had a slightly increased risk of CVD after controlling for age, CVD risk factors, and education level. 28
However, some research found no association of parity with CVD risk factors. A longitudinal study from the United Kingdom showed that there was no statistical association between parity and CVD risk factors (hypertension, type 2 diabetes, and BMI) after controlling for behavioral and lifestyle variables. 14 Chung et al. also reported that the association of parity with CVD risk factors (systolic and diastolic pressure, fasting plasma glucose, and insulin) had no significant association after controlling for covariates. 13 This inconsistency may be due to the differences in the average age of study populations in the studies.
Compared with nulliparous women, the CVD risk of parity 3 (OR:1.56; 95% CI: 1.05–2.32) was slightly lower than that for other parous women (Fig. 1). This was different from other research findings. The cohort studies from Sweden and Iran indicated a J-shaped relationship between parity and CVD, with parity of 2 having the lowest risk of CVD in parous women. 10,21 Another study from Britain also reported a similar J-shaped association between parity and CHD, women with parity 2 having the lowest prevalence of CHD. For those women with more than two children, each additional parity was associated with 30% greater age-adjusted odds of developing CHD (OR:1.30; 95% CI: 1.17–1.44). 15 So, further studies about U.S. women are needed to confirm the findings.
The specific risk factors for women to suffer from CVD might include pregnancy and reproductive-related endocrine disorders . 29 Female pregnancy can cause changes in the hormone environment, glucose and fat metabolism, and mild inflammation. 30,31 These changes may promote the development of CVD in later life. 32,33 In the normal third trimester of pregnancy, diastolic and systolic functions of the heart may be slightly impaired. Additionally, the insulin level of pregnant women could increase due to assisting in fetal development. 34,35 In the case of increased number of pregnancies and parity, all these changes may elevate the risk of CVD in women in the future.
In addition to biological mechanisms related to giving birth, the interpretation for the association between parity and CVD also involved factors related to raising children. 36 Some researchers have explored the differences in the association of parity with CVD between mothers and fathers. For example, Ness et al. reported the number of children was not associated with father's CVD in later life, although there was significant relationship between parity and maternal CVD. 37 This implied that the association between parity and maternal CVD may be from the biological mechanisms involved in pregnancy and giving birth. We were unable to examine the association between parity and paternal CVD because the National Health and Nutrition Examination Survey data do not include data on the number of children for fathers.
If we could examine the relationship between the number of children and paternal CVD, we may be able to better explain whether the association between parity and CVD was derived from parity-related physiological mechanisms or residual effects of socioeconomic factors. However, no matter what the cause of the relationship between parity and CVD, parous women may be high-risk populations of CVD compared with nulliparous women.
The strengths of this study included its representativeness of U.S. women and large sample size. The outcome variable was based on the study design where we only included participants in which stroke, heart attack, CHD, heart failure, and angina occurred after the last live birth, which could support the temporal relationship between parity and CVD. This could offset the shortcomings of cross-sectional data regarding causal identification to a certain extent. Furthermore, we performed the weighted analysis of the data, which made the results more applicable to the U.S. women. 38 Regarding the shortcomings of this study, information concerning family history of CVD only involved heart attack, not including stroke, heart failure, and other CVD conditions. Recall bias may exist because the data were self-reported from the questionnaire. In addition, due to the small sample size of nulliparous women with a history of pregnancy, we did not extract them as a separate subgroup during the analysis. This restricted us from comparing them to participants without any history of pregnancy.
Our findings suggest that women who gave births of more than four children may have an increased risk of developing CVD in later life in the United States. This may imply that preventive screening measures for CVD are necessary to favor the subpopulation. Of course, more research is imperative to further corroborate our findings, as well as more research to examine the underlying theoretical mechanism for the association. In addition, we would like to emphasize that although the results show that nulliparous women have the lowest risk of CVD, we discourage women from childbearing because nulliparity may have other adverse physical and mental health effects on women in later life.
Conclusion
Our study used a representative U.S. population-based sample to explore the association between parity and CVD in women over 45 years of age. There were strong positive associations between increased parity and CVD, stroke, heart attack, CHD, heart failure, and angina. After controlling for potential confounding variables, increased parity was statistically positively associated with CVD, stroke, and heart attack. Among all participants, nulliparous women had the lowest risk of CVD. In parous women, giving birth to three children seemed to have a slightly lower risk of CVD. Further cohort studies are warranted to confirm the findings.
Footnotes
Acknowledgment
The authors would like to thank Centers for Disease Control and Prevention, United States, for NHANES data.
Author Disclosure Statement
No competing financial interests exist.
Funding Information
No funding was received for this article.
