Skip to main content
  • Research article
  • Open access
  • Published:

Number of parity/live birth(s) and cardiovascular disease among Iranian women and men: results of over 15 years of follow-up

Abstract

Background

Most previous studies conducted in non-Middle Eastern populations have suggested that an increase in the number of parity/live birth(s) leads to cardiovascular disease (CVD) development, although their findings were inconclusive on this issue for both sexes. Biologic and socioeconomic pathways were suggested to explain this association. We studied this issue among urban Iranian men and women.

Methods

In this population-based cohort study, which included 3929 women and 2571 men aged ≥30 years, data for the number of parity/live birth(s) were obtained by a standard questionnaire. Participants were then annually followed for CVD events. Multivariable Cox proportional hazard models were used to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) for the number of parity/live birth(s) and other traditional CVD risk factors.

Results

During more than 15 years of follow-up, 456 and 524 CVD events have occurred among women and men, respectively. Among women, a J-shaped association was found between the number of live births and incident CVD with the lowest risk for women with two live births. Among women in multivariable analyses, each unit increase in parity had a HR of 1.05 (CI: 1.01–1.10), and having ≥4 parity was associated with a HR of 1.86 (0.97–3.56, p-value = 0.061). Among men, in comparison with participants who had 1 child, multivariable HRs of having 2, 3, and ≥ 4 children were 1.97 (1.24–3.12), 2.08 (1.31–3.31), and 2.08 (1.30–3.34), respectively.

Conclusion

To the best of our knowledge, the current study is the first report on this issue in the Middle East and North Africa region, a region with a high burden of CVD. It can now be suggested that the number of parity/live birth(s) is linked to CVD among the Iranian population, with this issue being more prominent among men. Further research is needed to support our results and clarify the pathways between the number of parity/live birth(s) and CVD development among Iranian populations by considering potential risk factors, especially psycho-socio-economic risk factors.

Peer Review reports

Background

According to the latest GBD (Global Burden of Disease) data, cardiovascular disease (CVD) is the most common non-communicable disease, responsible for 17.8 million deaths in 2017 [1]. Likewise, two consequences of CVD (stroke and coronary heart disease (CHD)) are the leading causes of mortality around the world [2, 3]. Statistics show that the effects of this phenomenon in the Middle East and North Africa (MENA) region are more prominent [2, 4]. In this region, Iran, while reporting modest mortality from CVD, is among the countries with the highest prevalence of CVD, with an approximate prevalence of 11% (7.7% for CHD) [2, 5].

In addition to traditional modifiable CVD risk factors such as obesity, hypertension, hypercholesterolemia, diabetes, and smoking, it was shown that some reproductive factors could have a role in CVD development [6]. Pregnancy, as one of the most important events during a women’s life, is associated with several cardiometabolic changes, including weight gain, dyslipidemia, increased plasma glucose, and insulin resistance. While these kinds of changes are beneficial to both mother and fetus, they can also increase the prevalence of potential CVD risk factors. Additional stressors during pregnancy such as endothelial dysfunction, inflammation, and hemostasis process may also exacerbate this situation [7,8,9]. Results of previous studies on the association of parity (i.e., number of pregnancies reaching viable gestational age) and gravidity (i.e., number of times a woman is or has been pregnant, regardless of the pregnancy outcome) with the risk of CVD are inconclusive [9, 10]. Peters et al. showed that among 10 European countries, parous women had a 19% higher risk of CHD than nulliparous women, with this risk reaching up to 95% among women who had ≥5 children [10]. On the other hand, some others reported no significant association between parity and CVD [11,12,13]. Moreover, a J-shaped relationship was reported between parity and CVD among a Swedish female population, in which women with 2 births had the lowest risk of CVD, as compared to nulliparous and other parous women [8].

In addition to biomedical factors, some socioeconomic and lifestyle factors, which have a two-way causal relationship with the number of children (family size), can play a role in CVD development [14, 15]. These include social and economic characteristics, family and social supports, educational level, employment, income, lifestyle etc. This socioeconomic and lifestyle pathway is applicable not only to mothers but also to fathers. A few previous studies examined the effect of the number of live births on CVD development among men; however, this issue is controversial among men as well [16, 17]. A large study conducted on a Chinese population reported that having more children was associated with CVD outcomes among men [17]. On the other hand, Eisenberg et al. found a negative association between the number of children and cardiovascular mortality among American men [16].

To the best of our knowledge, most previous studies related to these issues have been conducted in developed and Western countries, and there is a lack of information on the effect of the number of parity/live birth(s) on CVD development in the MENA region, a region with a high burden of CVD [2, 4]. The aim of the current study is to determine whether the number of parity/live birth(s) is an independent risk factor for incident CVD among Iranian women and men aged ≥30 years, over 15 years of follow-up. Data was collected from the oldest cohort of the MENA region, namely the Tehran Lipid and Glucose Study (TLGS).

Methods

Study design and study population

This study was conducted within the framework of TLGS, which is a community-based cohort study on a representative sample of Tehran’s residents. TLGS pursues aims to determine the prevalence and incidence of non-communicable diseases and their risk factors and also prevent these diseases by advancing healthier lifestyles. TLGS enrollment has been conducted in two phases; phase one (January 31, 1999- July 03, 2001) and phase two (October 20, 2001- September 22, 2005). Data collection is planned to continue for at least 20 years at approximately 3-year intervals. More details of TLGS design and enrollment have been reported elsewhere [18].

A total of 5223 women (phase I=4459, phase II=764) and 4330 men (phase I=3468, phase II=862)aged ≥30 years were enrolled. Participants were excluded if they were single (194 women and 179 men), with prevalent CVD at baseline (281 women and 318 men), and with no live births (21 women and 8 men), leading to a total of 4727 women and 3825 men. Moreover, other exclusions included those with missing data on the number of parity/live birth(s) (198 women and 1045 men), missing information on covariates including body mass index (BMI), fasting plasma glucose (FPG), total cholesterol (TC), systolic/diastolic blood pressure (SBP/DBP), education level, smoking status, family history of premature CVD, history of miscarriage, and history of oral contraceptive pill (OCP) use (198 women and 177 men, considering overlap features in the number of missing data between covariates). Also, subjects with no follow-up data (402 women and 32 men) were excluded. Finally, 3929 women and 2571 men were eligible for our analysis.

Clinical and laboratory measurements

Using a structured questionnaire, a trained nurse collected data for demographic data, past medical history, drug history, family history of premature CVD, education level, marital status, smoking habits, and physical activity level. Questionnaires also included data on the history of miscarriage and the number of parity/live birth(s).

According to TLGS setting [18], measurements of weight and height were done with shoes removed and wearing light clothing. Weight was measured to the nearest 100 g. The height of individuals was measured by a tape meter in a standing position. After a 15-min rest in a sitting position, two measurements of SBP and DBP were taken on the right arm. The mean of two measurements was defined as the subject’s blood pressure (BP).

On the day of blood collection, after 12 to 14 h of overnight fasting, a blood sample was collected between 7:00 and 9:00 AM from all participants in the TLGS research laboratory. Samples were analyzed on the same day. Details for laboratory measurements, including FPG, TC, high-density lipoprotein cholesterol (HDL-C), and triglycerides (TG) are published elsewhere [18].

Outcome assessment

CVD data collection has been described in detail elsewhere [18]. Briefly, a trained nurse telephoned all participants annually and asked them about any cardiovascular events that had occurred during the past year. For all reported events, a home visit was made by a trained physician; the physician collected data from medical documents or death certifications (in case of mortality). Finally, the outcome committee, which included an internist, an endocrinologist, a cardiologist, and an epidemiologist, evaluated the outcome data and adjudicated events. A CHD event included cases of unstable angina pectoris (new cardiac symptoms or changing symptom patterns and positive ECG findings with normal biomarkers), angiographic-proven CHD, definite myocardial infarction (MI) (diagnosed by electrocardiogram (ECG) and biomarkers), probable MI (positive ECG findings and cardiac symptoms plus missing biomarkers or positive ECG findings plus equivocal biomarkers), and cardiac death (any death in the hospital due to CHD based on the above-mentioned criteria, or sudden cardiac death by cardiac disease happening ≤1 h after symptoms initiation according to verbal autopsy files). CVD events was a combination of any CHD event, fatal or non-fatal stroke (defined as a new neurological deficit lasting ≥24 h), and cerebrovascular death.

Definition of terms

Type 2 diabetes mellitus (T2DM) was defined as one of these two criteria: i) FPG ≥7 mmol/L or ii) using glucose-lowering medications. Having TC ≥5.18 mmol/L or a history of lipid-lowering medications usage considered hypercholesterolemia. Hypertension was defined as either these three criteria: i) SBP ≥140 mmHg, ii) DBP ≥90 mmHg or iii) using antihypertensive medications. Smoking habit was categorized into three groups: i) current smoker, ii)former smokers and iii) never smokers. Education levels were classified into three categories: i) illiterate/primary school (reference group), ii) below diploma/diploma, and iii) above diploma. A history of CHD/stroke in a male first-degree relative aged < 55 years or female first-degree relative aged < 65 years was defined as a positive family history of premature CVD. For participants who were enrolled in phase I, the Lipid Research Clinic (LRC) questionnaire was used [19]; being physically active less than 3 days per week was considered as low physical activity. Using the Modifiable Activity Questionnaire (MAQ) for participants who were enrolled in phase II, subjects with < 600 MET (metabolic equivalent task-minutes per week) were in the low physical active group [18, 20, 21]. Parity was defined as the number of live childbirth plus the number of stillbirth (defined as the birth of an infant which died in the mother’s uterus after 20 weeks of gestation). A loss of an embryo or fetus before the 20th week of pregnancy defined as miscarriage.

Statistical analyses

All analyses were done separately for each sex. Baseline characteristics across the number of live births (1, 2, 3, and ≥ 4) are expressed as mean ± standard deviation (SD) for continuous and number (%) for categorical variables. To compare baseline characteristics among different groups, the ANOVA test (or Kruskal–Wallis test for skewed variables) was employed for continuous variables. The Chi-square test was also applied for categorical ones.

To be able to capture a potential nonlinear association between the number of live births and incident CVD, restricted cubic splines with 4 knots, which defined the 5th, 25th, 75th, and 95th percentile, were used. This method enabled us to modify the model across the number of live births while considering a dose-response relationship [22]. Cox proportional hazard models were applied to evaluate the association of the number of parity/live birth(s) with incident CVD in two models: model 1 adjusted with baseline measurements of age; model 2 further adjusted with BMI, T2DM, hypertension, hypercholesterolemia, educational levels, smoking status, and family history of premature CVD for both sexes, as well as history of miscarriage and OCP use for women. The hazard ratios (HRs) and 95% confidence intervals (95% CI) were reported for the number of parity/live birth(s) and other CVD risk factors. The time to event was defined as either the time of censoring or the event occurring, whichever came first. We censored participants in the case of not-CVD causes of death, leaving the district or being in the study until 20 March 2016, without any event. Moreover, we rerun our data analysis using age as time axis, rather than applied as a confounder in a sensitive analysis.

For the Cox models, the proportionality was assessed with the Schoenfeld residual test. All proportionality assumptions were appropriate. Statistical analyses were performed using STATA version 14 (StataCorp LP, College Station, Texas) statistical software. For defining statistical significance, P-values < 0.05 were considered.

Results

Baseline characteristics according to the number of live births are presented in Table 1 for both sexes. Generally, among continuous variables, cardiometabolic risk profiles became worse with increases in the number of live births. Thus, having more live births was associated with older age, higher BMI (only among women), increased BP, and higher levels of FPG, TC, and TG (only among women). Similarly, among categorical variables, having more live births was associated with higher prevalence of T2DM, hypercholesterolemia, and hypertension as well as higher use of glucose-lowering, lipid-lowering (only among women) and antihypertensive medications. Furthermore, a positive history of miscarriage became more prevalent with increases in the number of live births. Finally, as the number of live births increased, the percentage of those with a higher degree of education levels decreased significantly.

Table 1 Baseline characteristics according to the number of live births among women and men: Tehran Lipid and Glucose Study, Iran, 1999–2016

Among women, incident CVD occurred in 456 (11.6%) women during a median (interquartile range: IQR) follow-up of 15.9 (11.7–16.5) years. Among men, the incident CVD occurred in 524 (20.4%) men during a median (IQR) follow-up of 15.7 (11.3–16.4) years.

As shown in Fig. 1, we rejected the null hypothesis that CVD risk was a linear function of the mean number of live births. Among women, the associations between the number of live births and the risk of incident CVD were J-shaped. The association of the number of live births with incident CVD appeared to decrease until the number reached two live births. After this point, the association appeared to rebound. Furthermore, the association was non-linear among men.

Fig. 1
figure 1

Restricted cubic spline curve for association (95% confidence interval) of the number of live births with incident cardiovascular disease (CVD) among Tehranian women and men: Tehran Lipid and Glucose Study, Iran, 1999–2016

Multivariable HRs of incident CVD per additional live birth are shown in Table 2 for both sexes. When live birth was adjusted with age in model 1, each additional live birth was associated with HRs of 1.07 [CI: 1.03–1.12] for women and 1.02 [0.97–1.08] for men, which were significant among women. However, after further adjustment in model 2, the HRs for each additional live births reached 1.04 [0.99–1.09] among women and 1.01 [0.95–1.07] among men.

Table 2 Multivariable hazard ratios (HR) and 95% confidence intervals (CI) of incident CVD per additional live birth among women and men: Tehran Lipid and Glucose Study, Iran, 1999–2016

Compared to women with 1 live birth in model 1, women with ≥4 live births had an age-adjusted HR of 2.17 [1.18–4.00]; after more adjustment in model 2, the corresponding HR reached 1.72 [0.92–3.21, P-value = 0.087]. Among men, compared to those with 1 live birth, men with 2, 3, and ≥ 4 live births had a higher risk of incident CVD; the risks were significant in the full-adjusted model 2 by HRs of 1.97 [1.24–3.12], 2.08 [1.31–3.31], 2.08 [1.30–3.34], respectively. Importantly, beside the number of live births, traditional CVD risk factors including age, T2DM, hypertension, hypercholesterolemia, family history of premature CVD, former smoking (only among men), and current smoking were significantly associated with incident CVD in model 2 (Table 3).

Table 3 Multivariable hazard ratios (HR) and 95% confidence intervals (CI) of incident CVD according to the number of live births among women and men: Tehran Lipid and Glucose Study, Iran, 1999–2016

Moreover, we conducted a similar analysis among men aged ≥45 years, who had a high probability of having completed family size (the average number of live births was 3.11 in men aged ≥30 years vs. 4.07 for men aged ≥45 years). Similar to the men aged ≥30 years, among men older than 45 years, having ≥2 live births was associated with an increased CVD risk, with similar effect sizes, although the results did not reach a significant level (Supplementary Table 1). On the other hand, there was no difference between the average number of live birth in women aged ≥30 years (2.46) vs. aged ≥45 years (2.43).

Each unit increase in parity had an age-adjusted HR of 1.08 [1.04–1.13] among women in model 1; the HR remained at significant levels (1.05 [1.01–1.10]), even after further adjustment in model 2 (Table 4). Compared to women with 1 parity, the age-adjusted HRs were 0.97 [0.47–2.00] and 1.41 [0.72–2.76] among those with 2 and 3 parity, respectively. Moreover, having ≥4 parity had an age-adjusted HR of 2.34 [1.24–4.43], and after further adjustment in model 2, it also remained marginally significant (1.86 [0.97–3.56], P-value: 0.061). Importantly, beside the number of parity, traditional CVD risk factors including age, T2DM, hypertension, hypercholesterolemia, family history of premature CVD, and current smoking were significantly associated with incident CVD in model 2 (Table 5).

Table 4 Multivariable hazard ratios (HR) and 95% confidence intervals (CI) of incident CVD per additional parity among women: Tehran Lipid and Glucose Study, Iran, 1999–2016
Table 5 Multivariable hazard ratios (HR) and 95% confidence intervals (CI) of incident CVD according to the number of parity among women: Tehran Lipid and Glucose Study, Iran, 1999–2016

Compared to our main analysis (using follow-up time as the time axis), using age as the time axis yielded slightly lower effect sizes, but the results were overall similar and the 95% CIs overlapped; however, among women, the effect of parity as a continuous variable on incident CVD lost its significance in the sensitive analysis (Supplementary Table 2).

Discussion

This is the first population-based study conducted in the MENA region which examines the impact of the number of parity/live birth(s) on incident CVD events among both sexes during more than 15 years of follow-up.

Among women, the number of live births has a J-shaped association with incident CVD, with the lowest risk for those with 2 births. After adjustment for a wide series of important risk factors, including age, T2DM, hypertension, hypercholesterolemia, smoking status, and positive history of premature CVD -all of which remained significant risk factors in our analysis - each unit increase in parity was associated with a 5% higher risk of CVD events among women. Moreover, women with ≥4 parity had a more than 80% higher risk for incident CVD (marginally significant). Focusing on men, in comparison with participants who had 1 live birth, those with 2, 3, and ≥ 4 live births had about 100% higher risk for incident CVD in the presence of important traditional risk factors.

The findings of the current study about the association of parity number with incident CVD among women are consistent with a meta-analysis study on this issue by Li et al. [9] They found a non-linear J-shaped dose-response relationship between the number of parity and CVD among women. They also reported that each unit increase in the number of live birth led to a 4% increased risk of incident CVD; however, the authors showed significant high heterogeneity between studies (I2 = 89.6%). In another meta-analysis [23], Lv et al. also found a similar J-shaped association between parity number and CVD mortality with the lowest risk for women with 4 live births. They also found that each live birth was associated with a 1% non-significant increased risk for CVD mortality (I2 = 86.4%). It should be noted that the significant heterogeneity in these two meta-analyses [9, 23] could be related to different study populations, sample sizes, and other epidemiologic aspects of studies. Marginally significant increased CVD risk for women with ≥4 parity in our results agrees with the findings of other studies [24, 25]. Furthermore, our findings are in line with some previous cohort studies on CVD morbidity and mortality risk assessment, which show a J- or U-shaped association with the lowest risk for women with 2 live births or parity [8, 26]. However, in some others, there is no relation among women across the number of parity/live birth(s) [11, 13, 27], or their significance of association was lost after adjustment for other factors [12, 28, 29].

Although most of the previous studies on the current issue focused only on women, the limited data on men was controversial as well [16]. Similar to the current report, Peters et al. reported that in comparison with Chinese men who had 1 child, men with ≥2 children and men with ≥3 children had a higher risk of incident CHD and stroke, respectively. Among Chinese men, they also showed that each additional child significantly increased the risk of CHD and stroke by 3 and 2%, respectively [17]; however, we did not find a linear relationship between the number of live births and incident CVD among Tehranian men. On the other hand, Eisenberg et al. found a negative association between the number of offspring and CVD mortality. Based on their report, each additional child decreased the risk of CVD mortality by 2% among American men [16]. Furthermore, some previous studies could not detect any significant association between the number of live births and CVD mortality and morbidity among men [30, 31]. However, some others reported J- or U-shaped associations among men, which were similar to their female population study [26, 32].

A possible explanation for the association between parity and incident CVD in women is the biologic pathway. During pregnancy, some physiologic changes can have adverse effects for incident CVD which remained even after delivery, including weight gain, dyslipidemia, increased plasma glucose and insulin resistance as well as endothelial dysfunction and inflammatory and hemostatic processes [7,8,9]. In the light of multiple pregnancies, exposure time to these changes increased [9, 33]. This accumulative effect of repeating parity on traditional risk factors might be an important neglected residual confounding factor in the current study and similar researches in this field.

Beside the impact of the biologic pathway, psycho-socio-economic factors were reported to have a potentially important role in the pathway between the number of children and CVD development in both parents [14, 15]. In a previous study on a Swedish population, Barclay et al. compared the effect of the number of live births on CVD mortality between adoptive and biologic parents [14]. They showed that CVD mortality is higher in biological parents than adoptive parents, which means that the biologic pathway had inevitable effects. On the other hand, finding a U-shaped pattern in this study among adoptive male parents suggested that the biological pathway couldn’t be the only explanation for their findings, and the socioeconomic and lifestyle pathway should be considered as another explanation. Moreover, some previous studies [15, 17, 26] have found a similar pattern for men and women in their population study, which confirmed the strong role of socioeconomic and lifestyle pathway. It should be considered that during our study recruitment period, the policy of the Iranian government was based on reducing population growth, so there was minimal economic support for Iranian parents [34]. Furthermore, according to the data of the statistical center of Iran in 2017, household income per capita reduced with increased family size [35]. It means parents with more children might be under financial pressure. These economic problems could have adverse effects on socioeconomic status, leading to CVD development through psychosocial factors (anxiety and depression development), using unhealthier diet, smoking initiation, limiting leisure-time physical activity, poor access to health care, and little knowledge about diseases [36,37,38]. Since head of the Iranian family is traditionally the father [39], the stronger effect the of number of children among men is expectable through the socioeconomic and life style pathway; the issue needs serious and deep investigations in future studies.

Although we couldn’t find a strong protective role for childbearing in our data set, previous studies suggested that adult children may provide important social and financial support to their parents [14, 40]. These supports might be more important at older ages, where parents often need help for basic daily needs, accessing routine medical care, and prevention of social isolation [40]. Moreover, it was suggested that childbearing could cause developing a healthier lifestyle in parents [14, 40].

In addition to the CVD events, in previous studies, all-cause mortality [40, 41], diabetes [42], different types of cancer [43,44,45] were also found as medical outcomes influenced by childbearing. Therefore, it is impossible for us to accurately judge the ideal number of children in Tehran society. Moreover, in addition to the physical health of parents and families, other factors such as mental health, socioeconomic status, supports, employment, income, and life satisfaction should also be considered in the determination of the ideal children’s number. Hence, it seems that due to the large differences in the mentioned factors in the community, instead of talking about a fixed value of the ideal number of children for all families, it is better to determine the ideal number of children for each family individually based on a predictive model including all the potential effective factors.

The strengths of the current study consist of addressing the effect of the number of parity/live birth(s) on CVD development among both sexes in the MENA region for the first time, with standardized measurements for assessment of traditional risk factors rather than relying on self-reported data and adjudicated CVD events. The current study was limited in several ways. First, potential risk factors were considered at the time of baseline phases, and possible changes in risk factors as well as the number of parity/live birth(s), were not considered during the follow-up period. Second, the number of participants who hadn’t had children was too low to permit us to compare the impact of nulliparity with ever parity. Third, we did not have access to valid data on participant job status, income, and diet, which can clarify the socioeconomic and lifestyle pathway. Although we considered education status as a socioeconomic determinant, there was no strong relationship between education levels and economic status among the Iranian population [46]. Fourth, some other reproductive factors of female participants such as breastfeeding, age at first and last live birth, age at menarche, and age at natural menopause were found to have roles in CVD development. However, due to lack of data at the baseline recruitment, these reproductive factors were not accessible to us. Finally, our population study included only residents of a metropolitan city and cannot be generalized to rural populations.

Conclusion

To sum up, a J-shaped association was found between the number of live births and incident CVD among women. Moreover, for CVD events, besides the effects of traditional risk factors, each additional pregnancy increased the risk by 5%. We also found that women with ≥4 parity showed a marginally significant increased risk. Among men, those with only 1 child had the lowest risk of CVD events. It can be suggested that the number of parity/live birth(s) is linked to CVD among the Iranian population, with this issue being more prominent among men. Further research is needed to support our results and clarify the pathways between the number of parity/live birth(s) and CVD development among the Iranian population by considering potential factors, particularly psycho-socio-economic factors.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Abbreviations

GBD:

Global Burden of Disease

CVD:

Cardiovascular Disease

CHD:

Coronary Heart Disease

MENA:

Middle East and North Africa

TLGS:

Tehran Lipid and Glucose Study

BMI:

Body Mass Index

FPG:

Fasting Plasma Glucose

TC:

Total Cholesterol

SBP:

Systolic Blood Pressure

DBP:

Diastolic Blood Pressure

OCP:

Oral Contraceptive Pill

BP:

Blood Pressure

HDL-C:

High-density Lipoprotein Cholesterol

TG:

Triglycerides

MI:

Myocardial Infarction

ECG:

Electrocardiogram

T2DM:

Type 2 Diabetes Mellitus

LRC:

Lipid Research Clinic

MAQ:

Modifiable Activity Questionnaire

MET:

Metabolic Equivalent Task-minutes per Week

HR:

Hazard ratio

95% CI:

95% Confidence Interval

IQR:

Interquartile Range

References

  1. Roth GA, Abate D, Abate KH, Abay SM, Abbafati C, Abbasi N, et al. Global, regional, and national age-sex-specific mortality for 282 causes of death in 195 countries and territories, 1980–2017: a systematic analysis for the global burden of disease study 2017. Lancet. 2018;392(10159):1736–88.

    Article  Google Scholar 

  2. Roth GA, Johnson C, Abajobir A, Abd-Allah F, Abera SF, Abyu G, et al. Global, regional, and National Burden of cardiovascular diseases for 10 causes, 1990 to 2015. J Am Coll Cardiol. 2017;70(1):1–25.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Johnson CO, Nguyen M, Roth GA, Nichols E, Alam T, Abate D, et al. Global, regional, and national burden of stroke, 1990–2016: a systematic analysis for the global burden of disease study 2016. Lancet Neurol. 2019;18(5):439–58.

    Article  Google Scholar 

  4. Turk-Adawi K, Sarrafzadegan N, Fadhil I, Taubert K, Sadeghi M, Wenger NK, et al. Cardiovascular disease in the eastern Mediterranean region: epidemiology and risk factor burden. Nat Rev Cardiol. 2018;15(2):106–19.

    Article  PubMed  Google Scholar 

  5. Moazzeni SS, Ghafelehbashi H, Hasheminia M, Parizadeh D, Ghanbarian A, Azizi F, et al. Sex-specific prevalence of coronary heart disease among Tehranian adult population across different glycemic status: Tehran lipid and glucose study, 2008–2011. BMC Public Health. 2020;20(1):1510.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Peters SA, Woodward M. Women's reproductive factors and incident cardiovascular disease in the UK biobank. Heart. 2018;104(13):1069–75.

    Article  CAS  PubMed  Google Scholar 

  7. Ogunmoroti O, Osibogun O, Kolade OB, Ying W, Sharma G, Vaidya D, et al. Multiparity is associated with poorer cardiovascular health among women from the Multi-Ethnic Study of Atherosclerosis. Am J Obstet Gynecol. 2019;221(6):631.e631–16.

    Article  Google Scholar 

  8. Parikh NI, Cnattingius S, Dickman PW, Mittleman MA, Ludvigsson JF, Ingelsson E. Parity and risk of later-life maternal cardiovascular disease. Am Heart J. 2010;159(2):215–221.e216.

    Article  PubMed  Google Scholar 

  9. Li W, Ruan W, Lu Z, Wang D. Parity and risk of maternal cardiovascular disease: a dose–response meta-analysis of cohort studies. Eur J Prev Cardiol. 2019;26(6):592–602.

    Article  PubMed  Google Scholar 

  10. Peters SA, van der Schouw YT, Wood AM, Sweeting MJ, Moons KG, Weiderpass E, et al. Parity, breastfeeding and risk of coronary heart disease: a pan-European case-cohort study. Eur J Prev Cardiol. 2016;23(16):1755–65.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Colditz GA, Willett WC, Stampfer MJ, Rosner B, Speizer FE, Hennekens CH. A prospective study of age at menarche, parity, age at first birth, and coronary heart disease in women. Am J Epidemiol. 1987;126(5):861–70.

    Article  CAS  PubMed  Google Scholar 

  12. Jacobsen BK, Knutsen SF, Oda K, Fraser GE. Parity and total, ischemic heart disease and stroke mortality. The Adventist health study, 1976-1988. Eur J Epidemiol. 2011;26(9):711–8.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Chang HS, Odongua N, Ohrr H, Sull JW, Nam CM. Reproductive risk factors for cardiovascular disease mortality among postmenopausal women in Korea: the Kangwha cohort study, 1985-2005. Menopause. 2011;18(11):1205–12.

    Article  PubMed  Google Scholar 

  14. Barclay K, Kolk M. Parity and mortality: an examination of different explanatory mechanisms using data on biological and adoptive parents. Eur J Popul. 2019;35(1):63–85.

    Article  PubMed  Google Scholar 

  15. Magnus MC, Iliodromiti S, Lawlor DA, Catov JM, Nelson SM, Fraser A. Number of offspring and cardiovascular disease risk in men and women: the role of shared lifestyle characteristics. Epidemiology. 2017;28(6):880–8.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Eisenberg ML, Park Y, Hollenbeck AR, Lipshultz LI, Schatzkin A, Pletcher MJ. Fatherhood and the risk of cardiovascular mortality in the NIH-AARP diet and health study. Hum Reprod. 2011;26(12):3479–85.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Peters SA, Yang L, Guo Y, Chen Y, Bian Z, Millwood IY, et al. Parenthood and the risk of cardiovascular diseases among 0.5 million men and women: findings from the China Kadoorie biobank. Int J Epidemiol. 2017;46(1):180–9.

    PubMed  Google Scholar 

  18. 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:5.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Ainsworth BE, Jacobs DR Jr, Leon AS. Validity and reliability of self-reported physical activity status: the lipid research clinics questionnaire. Med Sci Sports Exerc. 1993;25(1):92–8.

    Article  CAS  PubMed  Google Scholar 

  20. Committee IR. Guidelines for data processing and analysis of the International Physical Activity Questionnaire (IPAQ)-short and long forms. http://www.ipaq.ki.se/scoring.pdf. 2005.

    Google Scholar 

  21. Momenan AA, Delshad M, Sarbazi N, Rezaei Ghaleh N, Ghanbarian A, Azizi F. Reliability and validity of the modifiable activity questionnaire (MAQ) in an Iranian urban adult population. Arch Iran Med. 2012;15(5):279–82.

    PubMed  Google Scholar 

  22. Royston P, Sauerbrei W. Multivariable modeling with cubic regression splines: a principled approach. Stata J. 2007;7(1):45–70.

    Article  Google Scholar 

  23. Lv H, Wu H, Yin J, Qian J, Ge J. Parity and cardiovascular disease mortality: a dose-response meta-analysis of cohort studies. Sci Rep. 2015;5:13411.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Zhang X, Shu X-O, Gao Y-T, Yang G, Li H, Zheng W. Pregnancy, childrearing, and risk of stroke in Chinese women. Stroke. 2009;40(8):2680–4.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Dior UP, Hochner H, Friedlander Y, Calderon-Margalit R, Jaffe D, Burger A, et al. Association between number of children and mortality of mothers: results of a 37-year follow-up study. Ann Epidemiol. 2013;23(1):13–8.

    Article  PubMed  Google Scholar 

  26. Barclay K, Keenan K, Grundy E, Kolk M, Myrskyla M. Reproductive history and post-reproductive mortality: a sibling comparison analysis using Swedish register data. Soc Sci Med. 2016;155:82–92.

    Article  PubMed  Google Scholar 

  27. Gallagher LG, Davis LB, Ray RM, Psaty BM, Gao DL, Checkoway H, et al. Reproductive history and mortality from cardiovascular disease among women textile workers in Shanghai, China. Int J Epidemiol. 2011;40(6):1510–8.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Koski-Rahikkala H, Pouta A, Pietilainen K, Hartikainen AL. Does parity affect mortality among parous women? J Epidemiol Community Health. 2006;60(11):968–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Vladutiu CJ, Meyer ML, Malek AM, Stuebe AM, Mosher A, Aggarwal S, et al. Racial differences in the association between parity and incident stroke: results from the REGARDS study. J Stroke Cerebrovasc Dis. 2017;26(4):749.

    Article  PubMed  Google Scholar 

  30. Ness RB, Cobb J, Harris T, D'Agostino RB. Does number of children increase the rate of coronary heart disease in men? Epidemiology. 1995;6:442–5.

    Article  CAS  PubMed  Google Scholar 

  31. Dekker J, Schouten E. Number of pregnancies and risk of cardiovascular disease. N Engl J Med. 1993;329(25):1893.

    Article  CAS  PubMed  Google Scholar 

  32. Jaffe DH, Eisenbach Z, Manor O. The effect of parity on cause-specific mortality among married men and women. Matern Child Health J. 2011;15(3):376–85.

    Article  PubMed  Google Scholar 

  33. Sanghavi M, Kulinski J, Ayers CR, Nelson D, Stewart R, Parikh N, et al. Association between number of live births and markers of subclinical atherosclerosis: the Dallas heart study. Eur J Prev Cardiol. 2016;23(4):391–9.

    Article  PubMed  Google Scholar 

  34. Abbasi-Shavazi MJ, McDonald P, Hosseini-Chavoshi M. The fertility transition in Iran, vol. 75: Springer; 2009.

    Google Scholar 

  35. Statistical-center-of-Iran. NATAYEGE AMARGIRI AZ HAZINE VA DARAMADE KHANEVARHAYE SHAHRI SALE 1396. 2017.

    Google Scholar 

  36. Clark AM, DesMeules M, Luo W, Duncan AS, Wielgosz A. Socioeconomic status and cardiovascular disease: risks and implications for care. Nat Rev Cardiol. 2009;6(11):712–22.

    Article  PubMed  Google Scholar 

  37. Williams J, Allen L, Wickramasinghe K, Mikkelsen B, Roberts N, Townsend N. A systematic review of associations between non-communicable diseases and socioeconomic status within low- and lower-middle-income countries. J Glob Health. 2018;8(2):020409.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Schultz WM, Kelli HM, Lisko JC, Varghese T, Shen J, Sandesara P, et al. Socioeconomic status and cardiovascular outcomes: challenges and interventions. Circulation. 2018;137(20):2166–78.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Aghajanian A. Family and family change in Iran. Fayetteville: Fayetteville State University; 2001.

    Google Scholar 

  40. Högnäs RS, Roelfs DJ, Shor E, Moore C, Reece T. J-curve? A meta-analysis and meta-regression of parity and parental mortality. Popul Res Policy Rev. 2017;36(2):273–308.

    Article  Google Scholar 

  41. Zeng Y, Ni Z-m, Liu S-y, Gu X, Huang Q, Liu J-a, et al. Parity and all-cause mortality in women and men: A dose-response meta-analysis of cohort studies. Sci Rep. 2016;6:19351.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Guo P, Zhou Q, Ren L, Chen Y, Hui Y. Higher parity is associated with increased risk of type 2 diabetes mellitus in women: a linear dose–response meta-analysis of cohort studies. J Diabetes Complicat. 2017;31(1):58–66.

    Article  Google Scholar 

  43. Wu Q-J, Li Y-Y, Tu C, Zhu J, Qian K-Q, Feng T-B, et al. Parity and endometrial cancer risk: a meta-analysis of epidemiological studies. Sci Rep. 2015;5:14243.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Lambe M, Hsieh C-c, Chan H-w, Ekbom A, Trichopoulos D, Adami H-O. Parity, age at first and last birth, and risk of breast cancer: a population-based study in Sweden. Breast Cancer Res Treat. 1996;38(3):305–11.

    Article  CAS  PubMed  Google Scholar 

  45. Zhu J, Zhu X, Tu C, Li YY, Qian KQ, Jiang C, et al. Parity and thyroid cancer risk: a meta-analysis of epidemiological studies. Cancer Med. 2016;5(4):739–52.

    Article  PubMed  Google Scholar 

  46. Falahati A, Fatahi S, Goli Y, Goli S. Education Demand and the Share of Educational in Income Inequality. Econ Model. 2017;11(37):23–44.

Download references

Acknowledgements

The authors would like to express their appreciation to the TLGS participants and staff for their kind cooperation. The authors also wish to acknowledge Dr. Laleh Rashidian for critical editing of English grammar and syntax of the manuscript.

Funding

No funding from any source was obtained for this study.

Author information

Authors and Affiliations

Authors

Contributions

Study conception and design: S.S.M and F.H; Analysis and interpretation of data: S. A and F.H; Drafting of manuscript: S.S.M, H. T, S.A and F.H; Critical revision: S.S.M, H. T, S. A, F. A, F.R.T and F.H. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Farzad Hadaegh.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Institutional Review Board (IRB) of the Research Institute for Endocrine Sciences (RIES), Shahid Beheshti University of Medical Sciences, and all participants provided written informed consent. All methods were done in accordance with the relevant guidelines and regulations.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Table S1.

Multivariable hazard ratios (HR) and 95% confidence intervals (CI) of incident CVD according to the number of live births among men aged ≥45 years: Tehran Lipid and Glucose Study, Iran, 1999–2016.

Additional file 2: Table S2.

Multivariable hazard ratios (HR) and 95% confidence intervals (CI) of incident CVD, comparing the use of follow-up time and age as axis for the Cox model: Tehran Lipid and Glucose Study, Iran, 1999–2016.

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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) 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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Moazzeni, S.S., Toreyhi, H., Asgari, S. et al. Number of parity/live birth(s) and cardiovascular disease among Iranian women and men: results of over 15 years of follow-up. BMC Pregnancy Childbirth 21, 28 (2021). https://doi.org/10.1186/s12884-020-03499-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12884-020-03499-2

Keywords