ABSTRACTObjectivesLow bone density (LBD) in the postmenopausal period has long been a pervasive public health concern; however, the association between parity and LBD has yet to be fully elucidated. Thus, we investigated the association between parity and LBD in postmenopausal Korean women.
MethodsThis study used baseline data from 1287 Korean postmenopausal women aged 40 years or older enrolled in the Cardiovascular and Metabolic Diseases Etiology Research Center community-based cohort study conducted in Korea from 2013 to 2017. The main exposure was parity (nullipara, 1, 2, 3+). The main outcome was LBD, including osteopenia and osteoporosis, based on bone mineral density measured using quantitative computed tomography of the lumbar spine (L1-2).
ResultsThe mean age of participants was 57.1 years, and the median parity was 2. Of the 1287 participants, 594 (46.2%) had osteopenia and 147 (11.4%) had osteoporosis. No significant difference in the prevalence of LBD was found between nullipara and parous women, whereas higher parity was associated with a higher risk of LBD among parous women; the adjusted odds ratio (95% confidence interval) for the presence of LBD was 1.40 (0.97 to 2.02) for a parity of 2 and 1.95 (1.23 to 3.09) for a parity of 3 relative to a parity of 1.
INTRODUCTIONPostmenopausal women are considered to be at high risk of developing osteoporosis, as estrogen deficiency accelerates bone turnover with net bone loss [1]. The prevalence of osteopenia and osteoporosis among Korean women aged 50 years and older is 48.7% and 38.0%, respectively [2]. Reducing the population burden of fractures requires paying more attention to women with osteopenia, because over half of fragility fractures in the general population arise in these individuals [3]. The bone mineral density (BMD) at every skeletal site in Korean men and women was found to be significantly lower than that of United States participants when comparing the results of the 2008-2011 Korea National Health and Nutrition Examination Survey (KNHANES) and the 2007-2009 United States National Health and Nutrition Examination Survey [2]. It has been estimated that the residual lifetime probability of osteoporosis-related fractures in 50-year-old Korean women is 59.5% [4]. Therefore, with the aging Korean population, the assessment of postmenopausal osteopenia and osteoporosis is particularly important.
Among the factors related to low bone density (LBD), parity and lactation remain controversial. Some studies have suggested that in postmenopausal women, a history of prolonged breastfeeding is an independent risk for osteoporosis [5,6]. Conversely, other studies have reported that the effect of breastfeeding on LBD in the postmenopausal period is not significant. It has been reported that BMD almost fully recovers after the restart of menstruation and weaning of lactation [7,8], and that spinal BMD significantly increased over time even when lactation continued for more than 6 months [7].
However, the effect of parity on the risk of LBD is not clear in the literature. Pregnancy and childbirth are factors that can lead to hormonal changes in the postpartum period, as well as throughout pregnancy. Parity also exerts a non-hormonal effect on BMD because lifestyle (e.g., nutrition, sleep, smoking, drinking, and physical activity) changes that can affect BMD occur during the parenting process following childbirth. Previous studies have reported diverse results regarding the association between parity and LBD. Some studies have reported a positive association between parity and the risk of osteoporosis or osteoporotic fracture [1], but others have not [9,10]. Another study reported that parous women had significantly higher BMD in the total hip than nullipara women, but there was no difference in the femoral neck or lumbar spine [11].
Thus, we investigated whether parity was associated with LBD among middle-aged postmenopausal women.
METHODSParticipants and Data SourceThis cross-sectional study was conducted using baseline data from the Cardiovascular and Metabolic Diseases Etiology Research Center (CMERC) study, a community-based cohort study conducted in Korea. The CMERC dataset includes data obtained from questionnaire interviews, physical examinations, neuropsychological testing, electrocardiography, a laboratory blood workup, and quantitative computed tomography (QCT). The CMERC baseline data included a sample of 2634 middle-aged (30-64 years) women with an obstetric history and BMD measured using QCT to examine associations between parity and LBD, including osteopenia and osteoporosis. Of these, 367 participants younger than 40 years and 500 premenopausal women were excluded. Participants who were on hormone replacement (n=395) or steroid therapy (n=4), which could affect BMD, were also excluded. Four participants with missing QCT results and 77 participants with missing covariate data were excluded. Thus, 1287 participants were included in the final analysis.
MeasurementsDual X-ray absorptiometry (DXA) and QCT are the most common tools for measuring BMD [12]. Compared to DXA, QCT is less affected by surrounding tissues and may avoid the overestimation of BMD on DXA that is associated with spinal degeneration, abdominal aortic calcification, and other sclerotic lesions [13]. When using QCT, spinal volumetric BMD (vBMD) is known to be more sensitive than BMD measurements at other sites [12]. Thus, we analyzed the QCT-measured spinal vBMD of L1 to L2 trabecular bone as the outcome.
The participants were scanned using a Somatom Definition AS+ 128-channel computed tomography (CT) scanner (Siemens Healthcare, Forchheim, Germany), a Somatom Sensation 64-channel CT scanner (Siemens Healthcare), or a General Electric LightSpeed VCT scanner (GE Medical Systems, Milwaukee, WI, USA) at 120 kVp and 150 mAs, using a 50-cm scan field of view. A liquid dipotassium phosphate (K2HPO4) phantom (Model 3; Mindways Software, Austin, TX, USA) was included in each scan to prevent scanner drift and to calculate the vBMD. Images obtained of the lumbar spine (L1-2) were reconstructed with a standard body reconstruction algorithm at 3-mm slice thickness, with an in-plane pixel size of 512×512 and a display field-of-view of 250 mm. All scanned data were analyzed using QCTPro software (Mindways Software). Lumbar trabecular vBMD was calculated as the average vBMD of the L1 and L2 trabecular bones.
We defined LBD as including osteopenia and osteoporosis. According to the official positions of the International Society for Clinical Densitometry 2019 and American College of Radiology 2018 revised guidelines [14], we defined LBD as BMD <120 mg/cm³. Osteoporosis was defined as a lumbar trabecular BMD ≤80 mg/cm³ or as taking anti-osteoporosis medication.
Height was measured within 0.1 cm using a stadiometer while the participant was standing in an upright position (Jenix, Seoul, Korea). Body weight was determined within 0.1 kg using an electronic scale (CAS, Seongnam, Korea). Body mass index (BMI) was calculated as body weight in kilograms divided by height in meters squared.
Trained examiners administered a questionnaire on obstetric history, demographic factors, lifestyle factors, and medical history in face-to-face interviews.
Information on reproductive health, including menopausal status, defined as cessation of menstruation for a minimum of 1 consecutive year, and the number of pregnancies and their outcomes (childbirth, stillbirth, or abortion) were obtained based on participants’ self-reports in face-to-face interviews. The postmenopausal period was calculated as the difference between the age at the time of examination and age at the time of menopause. Working status was classified as working or non-working. The non-working group included housewives, students, and unemployed individuals. Physical activity was assessed using the Korean version of the International Physical Activity Questionnaire standard to yield the metabolic equivalent of task (MET). Participants were classified into 3 groups according to the MET score: inactive, minimally active, and health-enhancing physical activity. Participants were categorized according to their cigarette smoking status as ever-smokers (current or ex-smokers) and non-smokers (never smokers). Drinking status was categorized into “yes” (current alcohol drinkers and former drinkers) and “no” (non-drinkers). Marital status was divided into ”unmarried” (never married, divorced, widowed, separated) and “married” (married and living with spouse). Education level was divided into middle school or less, high school, and university or beyond. Household income was classified into tertiles. A history of comorbidities, such as diabetes, hypertension, thyroid disease, and chronic kidney disease, was identified based on self-reported physician’s diagnoses. Current information about antidepressant and antiosteoporotic medication use was also obtained from participants’ self-reports.
The history of parity was obtained from participants who had been pregnant using the question, “How many total birth experiences have you had?” Participants were divided into four groups (nullipara, 1, 2, 3+) based on parity. Because there were only 17 participants with a parity of 4 or higher, they were merged into the 3+ category.
Statistical AnalysisRegression modeling was used to assess the associations between the predictors and LBD. First, we examined the difference in LBD between non-parous women (nullipara) and parous women using unadjusted and adjusted logistic regression models. We then examined the association between parity and LBD among parous women using unadjusted and adjusted logistic regression models. Age, BMI, postmenopausal period, working status, educational level, household income, marital status, smoking and alcohol intake, physical activity, and medical history were included as covariates in the adjusted models. We also performed a sensitivity analysis that restricted the outcomes only to osteoporosis using the fully adjusted model, as well as a subgroup analysis of the risk of LBD according to parity stratified by age group in tertiles. We additionally fitted generalized linear regression models to assess the association between parity and lumbar trabecular BMD (mg/cm³) as a continuous variable. Data analysis was performed from January to December 2020 using SAS version 9.4 (SAS Institute Inc., Cary, NC, USA).
RESULTSThe characteristics according to parity group are shown in Table 1. The mean age of the participants was 57.1 years. Only 40 (3.1%) women were nullipara, 153 (11.9%) were primipara, and 1094 (85.0%) were multipara.
We first evaluated the difference in LBD risk between nulliparous and parous women, but there was no significant difference (p=0.502) (Supplemental Material 1). Since non-parous women may have a history of pregnancy and abortion, we excluded nulliparous women from the following analysis to estimate LBD risk according to parity.
Women in the higher parity group had lower educational levels (p<0.001), a higher proportion of non-working status (p=0.002), and higher proportion with non-drinking status (p=0.032). These results could be explained by the fact that a higher number of older women were included in the higher parity group. The prevalence of LBD was 53% (10% osteoporosis+43% osteopenia) among non-parous women; the prevalence increased gradually according to parity among parous women: 48% (4% osteoporosis+44% osteopenia) for a parity of 1, 56% (8% osteoporosis+48% osteopenia) for a parity of 2, and 69% (13% osteoporosis+56% osteopenia) for a parity of 3+ (Figure 1).
The factors associated with LBD in the parous groups are shown in Table 2. The multivariable analysis revealed that compared to a parity of 1, the risk of LBD increased by 1.40 times and 1.95 times for parities of 2 and 3+, respectively. Older age, longer postmenopausal period, lower education level, and non-drinking status were also significantly associated with LBD in both the unadjusted and adjusted models (Table 2).
When we limited the outcome to osteoporosis, similar but stronger associations were observed between parity number and LBD; the odds ratio (95% confidence interval) was 2.46 (1.10 to 5.50) for a parity of 2 and 3.08 (1.30 to 7.30) for a parity of 3+ relative to a parity of 1 (Table 3).
We examined the association between parity and LBD in 3 subgroups stratified by age. In women aged 40–55 years and 56–59 years, similar to the trend for participants of all ages before stratification, women with higher parity had a higher risk of LBD. However, this trend was less significant in the older group (age 60–64 years), although we did not find a significant interaction between age group and parity (p=0.797) (Table 4). Since some studies reported that the risk of LBD increases exponentially with age, we conducted another sensitivity analysis with a model including age [2] but did not find different results (Supplemental Materials 2 and 3).
DISCUSSIONWe found no significant difference in LBD between non-parous and parous women, but observed that higher parity was associated with a higher risk of LBD among parous women.
Previous studies have also suggested that women with higher parity have a higher risk of LBD. A previous study evaluated the influence of parity on BMD and the risk of postmenopausal fracture and found that women with higher parity had significantly lower lumbar and femoral BMD after adjusting for age and BMI [15]. In another study, it was revealed that femoral BMD significantly decreased as parity increased, although no significant differences were found in the frequency of lumbar osteoporosis in the study group [16].
Some explanations have been proposed for the dose-relationship between parity and LBD. The recovery of bone density loss accelerates with the resumption of menstruation that was interrupted during pregnancy and breastfeeding. Women who experience repeated pregnancy, as well as women who have a short interval between pregnancies, have insufficient time to recover BMD as the number of pregnancies increases. In fact, pregnancy interval and number of pregnancies are independent factors associated with bone recovery, not breastfeeding [1]. This could explain the result that the risk of LBD increased depending on parity among parous women in our study.
We hypothesized that women with high parity may neglect aspects of lifestyle management, such as physical activity and nutrition, during parenting, which may increase the risk of LBD. However, in our study, having a high number of live births was not significantly associated with lower physical activity or higher BMI.
Since age was significantly and independently associated with decreased BMD [17], it is very important to consider age; however, an age-stratified association was not conducted in previous studies. In the fully adjusted model of the subgroup analysis following stratification into age tertiles, the intensity of the association between LBD and parity in the older age group was slightly diminished compared to the intensity in the other younger age groups. One possible explanation is a cohort effect. Each age group consisted of Korean women born between the 1950s and 1970s. During this period, Korea experienced sudden social changes immediately after the Korean war, such that women in each age group had quite different early life experiences; this could have affected their BMD according to their time of birth.
Nulliparous women have a heterogeneous background of not giving birth. Since various characteristics in nulliparous women are different from those in parous women, nulliparous women were excluded from the main analysis to estimate the risk of LBD according to parity.
The present study has some limitations. First, there was no breastfeeding history information in the data. Pregnancy and lactation affect bone metabolism and calcium homeostasis; during this period, increased intestinal calcium absorption and breast milk secretion cause calcium loss from the maternal skeleton, but bone loss could be restored within 6-12 months after weaning [1,18,19]. During 6 months following exclusive lactation and despite continued lactation, spinal BMD was found to increase significantly in a previous study [7]. Several recent studies have concluded that there is no significant association between decreased BMD and lactation, and that lactation does not increase the risk of fracture [19]. Therefore, it is believed that the significance of the effect of parity on LBD would not change if breastfeeding history had been included. Second, misclassification of parity was possible because women’s reproductive history was collected based on self-reports. Moreover, women’s reproductive information, such as that regarding abortions and pregnancies, may not be accurate because of a tendency for low rates of obstetric and gynecologic visits during pregnancy. Third, the CMERC cohort comprised a middle-aged (30-64 years) study population with no participants aged ≥65 years. Therefore, in this study, we could not analyze the risk of LBD according to parity number in older women, who have a high prevalence of osteoporosis and osteopenia. Fourth, there is a lack of generalizability when we compare the prevalence of LBD in our sample with a nationwide study sample. The prevalence of LBD in this study was lower than that in a 2008-2011 nationwide study for the 50-year to 59-year age group in each sample (CMERC: women aged 50-59 years: osteopenia 44.3%+osteoporosis 8.1%; KNHANES: women aged 50-59 years: osteopenia 58.6%+osteoporosis 15.4%) [2]. The lower prevalence of LBD in our study can be explained by the fact that BMD measured using DXA tends to be overestimated compared to the BMD measured using QCT, as in our study, and that our study sample comprised relatively healthy people. Lastly, calcium intake was not adjusted for the association between LBD and parity in this study. Calcium is an essential element in the diet, but its role in the pathogenesis of osteoporosis remains controversial [20]. In this study, daily calcium intake was not significantly associated with LBD (age-adjusted p-value=0.453). Based on previous literature and our analysis results, daily calcium intake was not included in the final model.
The main result of our analysis is that higher parity was associated with a higher risk of LBD in postmenopausal Korean women who had given birth. In particular, LBD of weight-bearing bones can lead to compression fractures. Accordingly, BMD monitoring is needed to detect LBD at an early stage and perform pharmacological and non-pharmacological interventions among postmenopausal women in Korea who have experienced multiple births.
In conclusion, we found a dose-related association between parity and LBD in Korean postmenopausal parous women. This association remained after adjusting for potential confounding covariates. Future studies should examine the longitudinal association between parity and the risk of developing osteoporosis and osteoporotic fractures.
Conflict of InterestCONFLICT OF INTEREST
The authors have no conflicts of interest associated with the material presented in this paper.
NotesFUNDING
This study was supported by a grant (HI13C0715) from the Korea Health Technology R&D Project, funded by the Ministry of Health and Welfare of the Republic of Korea.
NotesAUTHOR CONTRIBUTIONS
Conceptualization: ES, HCK. Data curation: ES, HCK, YL. Formal analysis: ES, HCK. Funding acquisition: HCK. Methodology: ES, HCK. Visualization: ES, HCK, YL. Writing – original draft: ES, HCK. Writing – review & editing: ES, HCK, YL.
Figure. 1.Frequency of osteopenia and osteoporosis of the lumbar spine in each parity group. Values are presented as number (%). ![]() Table 1.General characteristics of the study population according to parity Table 2.Factors associated with LBD in parous postmenopausal women
Table 3.Factors associated with osteoporosis in parous postmenopausal women
Table 4.Associations between parity and low bone density (LBD) according to age group
SUPPLEMENTAL MATERIALSSupplemental materials are available at https://doi.org/10.3961/jpmph.21.162.
Supplemental Material 1.Risk of low bone density in parous versus nulliparous women Supplemental Material 2.Associations between parity and low bone density according to age and age² controlled in the fully adjusted model Supplemental Material 3.Association between parity and osteoporosis according to age and age² controlled in the fully adjusted model REFERENCES1. Yan G, Huang Y, Cao H, Wu J, Jiang N, Cao X. Association of breastfeeding and postmenopausal osteoporosis in Chinese women: a community-based retrospective study. BMC Womens Health 2019;19(1):110.
![]() ![]() ![]() 2. Park EJ, Joo IW, Jang MJ, Kim YT, Oh K, Oh HJ. Prevalence of osteoporosis in the Korean population based on Korea National Health and Nutrition Examination Survey (KNHANES), 2008-2011. Yonsei Med J 2014;55(4):1049-1057.
![]() ![]() ![]() 3. Pasco JA, Seeman E, Henry MJ, Merriman EN, Nicholson GC, Kotowicz MA. The population burden of fractures originates in women with osteopenia, not osteoporosis. Osteoporos Int 2006;17(9):1404-1409.
![]() ![]() 4. Park C, Ha YC, Jang S, Jang S, Yoon HK, Lee YK. The incidence and residual lifetime risk of osteoporosis-related fractures in Korea. J Bone Miner Metab 2011;29(6):744-751.
![]() ![]() 5. Hwang IR, Choi YK, Lee WK, Kim JG, Lee IK, Kim SW, et al. Association between prolonged breastfeeding and bone mineral density and osteoporosis in postmenopausal women: KNHANES 2010-2011. Osteoporos Int 2016;27(1):257-265.
![]() ![]() 6. Lee EN. Effects of parity and breastfeeding duration on bone density in postmenopausal women. Asian Nurs Res (Korean Soc Nurs Sci) 2019;13(2):161-167.
![]() ![]() 7. Cooke-Hubley S, Kirby BJ, Valcour JE, Mugford G, Adachi JD, Kovacs CS. Spine bone mineral density increases after 6 months of exclusive lactation, even in women who keep breastfeeding. Arch Osteoporos 2017;12(1):73.
![]() 8. More C, Bettembuk P, Bhattoa HP, Balogh A. The effects of pregnancy and lactation on bone mineral density. Osteoporos Int 2001;12(9):732-737.
![]() ![]() 9. Crandall CJ, Liu J, Cauley J, Newcomb PA, Manson JE, Vitolins MZ, et al. Associations of parity, breastfeeding, and fractures in the Women’s Health Observational Study. Obstet Gynecol 2017;130(1):171-180.
![]() ![]() ![]() 10. Wang Q, Huang Q, Zeng Y, Liang JJ, Liu SY, Gu X, et al. Parity and osteoporotic fracture risk in postmenopausal women: a dose-response meta-analysis of prospective studies. Osteoporos Int 2016;27(1):319-330.
![]() ![]() 11. Song SY, Kim Y, Park H, Kim YJ, Kang W, Kim EY. Effect of parity on bone mineral density: a systematic review and meta-analysis. Bone 2017;101: 70-76.
![]() ![]() 12. Ma XH, Zhang W, Wang Y, Xue P, Li YK. Comparison of the spine and hip BMD assessments derived from quantitative computed tomography. Int J Endocrinol 2015;2015: 675340.
![]() ![]() ![]() 13. Li N, Li XM, Xu L, Sun WJ, Cheng XG, Tian W. Comparison of QCT and DXA: osteoporosis detection rates in postmenopausal women. Int J Endocrinol 2013;2013: 895474.
![]() ![]() ![]() 14. American College of Radiology. ACR–SPR–SSR practice parameter for the performance of musculoskeletal quantitative computed tomography (QCT). 2018 [cited 2021 Mar 1]. Available from: https://www.acr.org/-/media/ACR/Files/Practice-Parameters/qct.pdf.
15. Allali F, Maaroufi H, Aichaoui SE, Khazani H, Saoud B, Benyahya B, et al. Influence of parity on bone mineral density and peripheral fracture risk in Moroccan postmenopausal women. Maturitas 2007;57(4):392-398.
![]() ![]() 16. Peker N, Tosun ÖÇ. Is grand multiparity a risk factor for the development of postmenopausal osteoporosis? Clin Interv Aging 2018;13: 505-508.
![]() ![]() ![]() 17. Kritz-Silverstein D, Barrett-Connor E, Hollenbach KA. Pregnancy and lactation as determinants of bone mineral density in postmenopausal women. Am J Epidemiol 1992;136(9):1052-1059.
![]() ![]() 18. Kalkwarf HJ, Specker BL. Bone mineral loss during lactation and recovery after weaning. Obstet Gynecol 1995;86(1):26-32.
![]() ![]() |
|