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

Parity and bladder cancer risk: a dose-response meta-analysis



Multiple studies have reported evidence of an inverse association between parity and bladder cancer risk. However, a comprehensive and quantitative assessment of this association has never been conducted. We conducted this study to clarify this issue.


Systematic search of PubMed and Embase was performed to identify all the studies. Studies were selected based on strict screening with inclusion and exclusion criteria. Summary relative risks (RR) with 95% confidence intervals (CI) were calculated by using a fixed-effect model, and the generalized least squares trend estimation was employed to compute study-specific RR and 95% CI per live birth increase. Heterogeneity and publication bias were also evaluated.


Twelve studies (6,214 cases and 2,693,350 non-cases) were eligible in this meta-analysis. The pooled RR of bladder cancer for parous versus nulliparous women was 0.76 (95% CI: 0.70–0.82). Results were similar in the studies that adjusted for BMI(RR = 0.66; 95% CI: 0.53–0.81), cigarette smoking (RR = 0.67; 95% CI: 0.57–0.79), and age (RR = 0.77; 95% CI: 0.71–0.84). The dose-response meta-analysis showed a lower bladder cancer risk (RR = 0.95; 95% CI: 0.92–0.98) for each live birth increase in parous women. No evidence of publication bias or significant heterogeneity was detected in the above-mentioned analyses.


The finding from current meta-analysis suggest that parity may be related to decreased risk of bladder cancer.

Peer Review reports


Bladder cancer (BC) is one of the most common genitourinary malignancies in the United States, with an estimated 76, 960 newly diagnosed cases and 16 390 deaths in 2016 [1]. Incidence rate of BC is notably gender-specific with men having a 3~4 fold greater risk than women for developing the disease [1]. Given the evidence from current knowledge, the gender disparity of incidence rate for BC cannot be completely explained by known risk factors, such as cigarette smoking and occupation [2].

Although the exact biological mechanisms for the gender disparity in BC development and progression are not well established, some studies have suggested that sex hormones and their receptors may contribute to this difference [35]. During pregnancy, changes in maternal hormones, as well as physical structure changes in the pelvic or lower abdominal region following childbirth, may lead to etiological changes that decrease BC risk. Till date, numerous observational studies [611] have investigated the roles of reproductive factors, such as age at menarche and menopause, parity, age at first birth, and number of children, in the development of BC, but results have been inconsistent, probably because of limited sample size included each individual study. Although published reviews and meta-analyses have focused on this topic [12, 13], up to now, a comprehensive and quantitative assessment of the relationship between parity number and BC risk has not been reported.

Therefore, we find it necessary to further assess the relationship between parity and BC risk by conducting a meta-analysis of current epidemiologic studies and providing a quantitative dose-response analysis.


This study protocol was approved by the institutional review board at the West China School of Medicine before initiation and don’t need ethical standard statement.

Literature search and eligibility criteria

We searched online databases MEDLINE and EMBASE for the studies assessing the relationship of parity and BC risk before Sep 2016, using the following terms with every possible combination considered: parity, pregnancy, reproductivity, fertility, bladder, genitourinary tract, cancer, carcinoma, tumor. We searched the references of all eligible publications to identify additional relevant publications. Studies were included if they met the following criteria: (1) the study had a prospective or case-control design; (2) the study evaluated the relationship between parity and BC risk; (3) data were provided or would allow the calculation of relative risk (RR) or odds ratio (OR) and a 95% confidence intervals (CI). Exclusion criteria were as follows: (1) was not involved with the associations between parity and BC risk; (2) non-human studies, reviews, and comments; (3) studies not reporting primary outcomes; (4) studies based on overlapping patients. When duplicate publications were identified or data in separate publications originated from the same cohort, we used the publication with the largest number of cases and the most applicable data.

Data extraction and quality assessment

Data were independently extracted from the included studies by two investigators (YJB and XMW), and disagreement was settled by discussion with the third investigator. The following data were extracted from each study: first author’ name, publication year, country of study, participant demographics, sample size, duration of follow-up or study period, parity number categories, corresponding OR/RR (with their 95% CIs) for each category, and variables adjusted for in the analysis.

The methodological quality assessment was conducted for each of the included studies using the Newcastle-Ottawa Scale quality assessment [14] on three broad perspectives: selection, comparability and exposure or outcome. Two investigators read each included study and scored them independently. Disagreements were resolved by consensus.

Statistical analysis

Given that bladder cancer is a rare disease, the OR was assumed to be the same as the RR and we therefore report all effect sizes as RR for simplicity. We calculated summary RRs and 95% CIs associated with parity and BC risk, if the study considered nulliparous as a reference. Between-study heterogeneity was evaluated with Q and I2 statistics [15]. For the Q-test, we used P < 0.10 as evidence of heterogeneity. An I2 score exceeding 50% is considered to indicate the presence of heterogeneity. The random effect model was used to provide summary estimations if there was heterogeneity, or else the fixed effect model was used. Heterogeneity was explored by stratified analysis. Publication bias was explored using Begg’s and Egger’s test [16, 17] and funnel plots. P < 0.05 indicated the existence of publication bias.

We carried out a dose-response meta-analysis using the method proposed by Greenland et al. [18] and Orsini et al. [19] to estimate study-specific slopes (linear trends) and 95% CIs from the log-RRs and CIs across the categories of parity number. This analysis requires the number of cases and controls for case-control studies or the number of cases and person-years for cohort studies, and the RR with 95% CI for at least three quantitative categories of parity number are presented. The value assigned to each category was the midpoint of the upper and lower boundaries and was adjusted for half range of the neighborhood categories when categories were open-ended. Random-effects model was applied in our study. This was done by modeling parity number using restricted cubic splines with three knots at percentiles 5, 50, and 95% of the distribution. A P value for nonlinearity was computed by testing the null hypothesis that the coefficient of the second spline was equal to zero. Statistical analyses were performed with Stata version 12.0 software (Stata Corp, College Station, TX).


Study characteristics and quality assessment

The article selection process is shown in Fig. 1. Thirteen articles [613, 2024] met the selection criteria. Two studies reported data on the same cohort [9, 22] and we included the more recently one which reported the largest number of cases [22]. Thus, 12 publications were included in this meta-analysis [68, 1013, 2024], and their characteristics are presented in Table 1. These studies included a total of 6,214 cases and 2,693,350 non-cases and were published from 1992 to 2013. Among these studies, one study reported three separate outcomes: the Los Angeles Bladder Cancer Study, the Shanghai Bladder Cancer Study, and the California Teachers Study [13]. Thus, the present study used seven outcomes from case-control studies and seven outcomes from cohort studies to produce a meta-analysis. Of the included studies, 13 were conducted in Europe, nine in America and one in China. All studies included met quality criteria ranging from 5 to 7 stars.

Fig. 1
figure 1

Flow diagram of the study selection process

Table 1 Characteristics of eligible studies included in this meta-analysis of parity and bladder cancer risk

Ever versus never parity

All included studies in present study investigated the association between ever parity and BC risk. The summary multivariable-adjusted RR of BC for ever versus nulliparous was 0.76 (95% CI, 0.70–0.82), without heterogeneity (I2 = 0%, p = 0.795, Fig. 2), indicating an inverse association between parity and BC risk. For women who have never smoked, the summary RR of BC for ever versus nulliparous was 0.47 (95% CI, 0.35–0.63), without heterogeneity (I2 = 0%, p = 0.562). For women who have ever smoked, the summary RR was 0.90 (95% CI, 0.67–1.21) for ever versus nulliparous, without heterogeneity (I2 = 0%, p = 0.531).

Fig. 2
figure 2

Forest plot (fixed-effects model) of ever parity and bladder cancer risk

Differences with parity number

Six cohort and six case-control studies investigated possible association between parity and BC risk. The results of the effects of different parity number on BC risk are presented in Table 2. To clarify the effects, we divided the parity number into three groups. The first group (1–2 births) contained seven reports, compared with nulliparous,, the pooled RR of BC associated with giving birth to two children was 0.82 (95% CI 0.71–0.94) with I2 = 17.7% (P = 0.295). The second group (3–4 births) contained seven reports, compared with nulliparous, the overall RR of BC was 0.79 (95% CI 0.68–0.91) with I2 = 0% (P = 0.997). The third group (≥5 births) contained six reports, compared with nulliparous, the overall RR of BC was 0.76 (95% CI 0.66–0.88) with I2 = 0% (P = 0.994).

Table 2 Summary risk estimates of the association between parity and bladder cancer

Dose-response meta-analysis

Dose-response from nine studies [6, 8, 1113, 2023] showed a decreased in BC risk of 0.95 (95% CI: 0.92–0.98) per live birth increase, without heterogeneity (I2 = 12.79%, P = 0.7832). Compared with nulliparous individuals, the pooled RRs (95% CI) of BC were 0.84 (0.76–0.94), 0.76 (0.67–0.87), 0.74 (0.65–0.85), 0.75 (0.66–0.85), and 0.76 (0.65–0.88) for 1, 2, 3, 4, and 5 live births, respectively (Fig. 3).

Fig. 3
figure 3

Nonlinear dose-response relationship between parity and bladder cancer risk

Publication bias

No evidence of a significant publication bias was observed in our analyses as assessed using the Egger’s test or Begg’s test and no asymmetry was seen in the funnel plots of included studies.

Subgroup analyses

To avoid the influence of cigarette smoking, we evaluated the influence of smoking status by adjustment in the eight reports that considered smoking. However, when data were stratified by adjustment for smoking, we did not find a significant difference between summary RRs that were adjusted and those that were not adjusted for smoking (Table 2). Moreover, significant inverse associations were persisted after the subgroup analyses stratified by whether the study included adjustment for specific potential confounders, such as age, Body Mass Index (BMI), type of control subjects or others potential confounders (Table 2). An additional file shows this in more detail (see Additional file 1).


Numerous epidemiological studies have focused on BC risk factors in order to explain the gender disparity in BC incidence, while there is also a great interest in identifying factors that affect the risk of BC among women. Although previous studies have reported on the associations between parity and BC, to the best of our knowledge, no other earlier studies have clearly documented a dose-response pattern between parity number and BC risk. In present study, we conducted a dose-response meta-analysis to characterize the association between parity and BC risk. Findings from the present study indicated that parous was significantly inversely associated with BC risk compared with nulliparous, especially for non-smokers. Overall, the risk of BC decreased by 24% for women who had given birth and 53% for those who had never smoked and were ever parous.

Although the exact mechanisms underlying this inverse association between parity and BC risk are not completely established, several potential mechanisms have been proposed. First, during pregnancy, estrogen and progesterone levels increase drastically and induce substantial alterations to the structure, function, and histology of the bladder [25], possibly affecting the development and progression of BC. Both estrogen and progesterone receptors are found in bladder tissues suggesting that endocrine regulation could directly influence BC development [26]. Shen et al. [27] found that antiestrogens have an inhibitory effect on the growth of BC cells in vitro, indicating that estrogens may increase the development and growth of bladder malignancies. Some research has reported that progesterone suppresses the activity of the estrogen receptor during pregnancy [28]. Therefore, estrogen and progesterone have an antagonistic effects on human cells. The substantial increased in their levels during pregnancy may responsible for the decreased risk of BC for women. Moreover, a hormone-related protective effect based on the prity number may mirror the lifetime period of lactation and prolonged oxytocin inhibition of steroid hormones [6].

Second, smoking habits may have influenced the incidence of BC in women with parity. It is well known that prenatal cigarette smoking exposure may cause detrimental effects on reproductive health [29]. Therefore, most smokers who prepare for pregnancy are more likely to quit smoking. Thus, compared with nulliparous female smokers, parous smokers who quit smoking may have a shorter smoking exposure history and reduction of pack years, and this may affect their BC risk. In addition, cumulative effect of abandon smoking following multiple pregnancies may also associate to decrease the subsequent risk of BC. The findings of our study suggest that the negative correlation between parity and BC appears more remarkable among women who were never exposed to cigarettes, with these parous women experiencing at least a 53% decreased risk of BC compared with nulliparous women. However, such a dramatic risk reduction was not apparent among smokers. The reasons for this results of subgroup analysis may be that the effect of smoking on BC is so huge that for the development of the disease minor influence factors such as parity are less relevant. In case of never smokers, these minor influence factors play indeed a role.

Third, structural changes in the pelvic organs and pelvic floor organization following pregnancy may induce a decreased risk for BC in parous women. Childbirth is an established risk factor for lower urinary tract dysfunction, such as increased urinary frequency, stress incontinence, and overactive bladder [30], but the possible relationship with BC remains unclear. It is well known that BC risk is related to the contact time of the urothelium with carcinogens in urine. Some study reported that increased water intake and urination frequency may lessen contact time of urinary carcinogens with urothelium of the bladder, thus, diminish the bladder cancer risk [31, 32]. Zhang et al. [31] found that increasing urination frequency during daytime was associated with decreased risk of bladder cancer. Silverman et al. [32] suggested that nocturia may be a powerful factor in reducing BC risk. Based on the results of above studies, we can be reasonable to assume that lower urinary tract dysfunction after childbirth has a protective effect on BC risk, and this may be associated with shorter duration of exposure to urine carcinogens in the bladder. Further research should verify this hypothesis.

Meta-analyses of observational studies are prone to confounding, selection, and information bias as the included observational studies and can present particular challenges that may distort the results. However, our study has some strengths. First, we searched 2 electronic databases, which ensured a broad scope, and we included seven cohorts and seven case-control studies that involved a total of 6,214 cases and 2,693,350 non-cases. This large sample size should have provided sufficient statistical power to detect this potential relationship. Second, to control misclassification, we conducted a dose-response meta-analysis. This is the most comprehensive meta-analysis evaluating the associations between parity and BC. Third, no publication biases were detected in present study and outcome may be unlikely to remain unreported as studies with complex assessment of exposure.


This current dose-response meta-analysis indicates that women with parity have an inversely associated BC risk compared with nulliparous women. The exact mechanism underlying this protective effect requires further investigation.



Bladder cancer


Body mass index


Confidence interval


Odds ratio


Relative risk


  1. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2016. CA Cancer J Clin. 2016;66:7–30.

    Article  PubMed  Google Scholar 

  2. Hartge P, Harvey EB, Linehan WM, Silverman DT, Sullivan JW, Hoover RN, Fraumeni Jr JF. Unexplained excess risk of bladder cancer in men. J Natl Cancer Inst. 1990;82:1636–40.

    Article  CAS  PubMed  Google Scholar 

  3. Bolenz C, Lotan Y, Ashfaq R, Shariat SF. Estrogen and progesterone hormonal receptor expression in urothelial carcinoma of the bladder. Eur Urol. 2009;56:1093–5.

    Article  PubMed  Google Scholar 

  4. Lucca I, Fajkovic H, Klatte T. Sex steroids and gender differences in nonmuscle invasive bladder cancer. Curr Opin Urol. 2014;24:500–5.

    Article  PubMed  Google Scholar 

  5. Gakis G, Stenzl A. Gender-specific differences in muscle-invasive bladder cancer: the concept of sex steroid sensitivity. World J Urol. 2013;31:1059–64.

    Article  CAS  PubMed  Google Scholar 

  6. Weibull CE, Eloranta S, Altman D, Johansson AL, Lambe M. Childbearing and the risk of bladder cancer: a nationwide population-based cohort study. Eur Urol. 2013;63:733–8.

    Article  PubMed  Google Scholar 

  7. Huang AT, Kogevinas M, Silverman DT, Malats N, Rothman N, Tardón A, Serra C, García-Closas R, Carrato A, Cantor KP. Bladder cancer and reproductive factors among women in Spain. Cancer Causes Control. 2009;20:1907–13.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Cantwell MM, Lacey Jr JV, Schairer C, Schatzkin A, Michaud DS. Reproductive factors, exogenous hormone use and bladder cancer risk in a prospective study. Int J Cancer. 2006;119:2398–401.

    Article  CAS  PubMed  Google Scholar 

  9. Tripathi A, Folsom AR, Anderson KE, Iowa Women’s Health Study. Risk factors for urinary bladder carcinoma in postmenopausal women. The Iowa Women’s Health Study. Cancer. 2002;95:2316–23.

    Article  PubMed  Google Scholar 

  10. Pelucchi C, La Vecchia C, Negri E, Dal Maso L, Franceschi S. Smoking and other risk factors for bladder cancer in women. Prev Med. 2002;35:114–20.

    Article  PubMed  Google Scholar 

  11. McGrath M, Michaud DS, De Vivo I. Hormonal and reproductive factors and the risk of bladder cancer in women. Am J Epidemiol. 2006;163:236–44.

    Article  PubMed  Google Scholar 

  12. Dietrich K, Demidenko E, Schned A, Zens MS, Heaney J, Karagas MR. Parity, early menopause and the incidence of bladder cancer in women: a case-control study and meta-analysis. Eur J Cancer. 2011;47:592–9.

    Article  CAS  PubMed  Google Scholar 

  13. Davis-Dao CA, Henderson KD, Sullivan-Halley J, Ma H, West D, Xiang YB, Gago-Dominguez M, Stern MC, Castelao JE, Conti DV, Pike MC, Bernstein L, Cortessis VK. Lower risk in parous women suggests that hormonal factors are important in bladder cancer etiology. Cancer Epidemiol Biomarkers Prev. 2011;20:1156–70.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Stang A. Critical evaluation of the Newcastle-Ottawa scale for the assessment of the quality of nonrandomized studies in meta-analyses. Eur J Epidemiol. 2010;25:603–5.

    Article  PubMed  Google Scholar 

  15. Higgins JP, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in meta-analyses. BMJ. 2003;327:557–60.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Begg CB, Mazumdar M. Operating characteristics of a rank correlation test for publication bias. Biometrics. 1994;50:1088–101.

    Article  CAS  PubMed  Google Scholar 

  17. Egger M, Davey SG, Schneider M, Minder C. Bias in meta-analysis detected by a simple, graphical test. BMJ. 1997;315:629–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Greenland S, Longnecker MP. Methods for trend estimation from summarized dose-response data, with applications to meta-analysis. Am J Epidemiol. 1992;135:1301–9.

    CAS  PubMed  Google Scholar 

  19. Orsini N, Bellocco R, Greenland S. Meta-analysis for linear and nonlinear dose-response relations: examples, an evaluation of approximations, and software. Am J Epidemiol. 2012;175:66–73.

    Article  PubMed  Google Scholar 

  20. Cantor KP, Lynch CF, Johnson D. Bladder cancer, parity, and age at first birth. Cancer Causes Control. 1992;3:57–62.

    Article  CAS  PubMed  Google Scholar 

  21. La Vecchia C, Negri E, Franceschi S, Parazzini F. Long-term impact of reproductive factors on cancer risk. Int J Cancer. 1993;53:215–9.

    Article  PubMed  Google Scholar 

  22. Prizment AE, Anderson KE, Harlow BL, Folsom AR. Reproductive risk factors for incident bladder cancer: Iowa Women’s Health Study. Int J Cancer. 2006;120:1093–8.

    Article  Google Scholar 

  23. Daugherty SE, Lacey Jr JV, Pfeiffer RM, Park Y, Hoover RN, Silverman DT. Reproductive factors and menopausal hormone therapy and bladder cancer risk in the NIH-AARP Diet and Health Study. Int J Cancer. 2013;133:462–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Kabat GC, Kim MY, Luo J, Hou L, Cetnar J, Wactawski-Wende J, Rohan TE. Menstrual and reproductive factors and exogenous hormone use and risk of transitional cell bladder cancer in postmenopausal women. Eur J Cancer Prev. 2013;22:409–16.

    Article  CAS  PubMed  Google Scholar 

  25. Rodrìguez LV, Wang B, Shortliffe LM. Structural changes in the bladder walls of pregnant and hormone-treated rats: correlation with bladder dynamics. BJU Int. 2004;94:1366–72.

    Article  PubMed  Google Scholar 

  26. Blakeman PJ, Hilton P, Bulmer JN. Cellular proliferation in the female lower urinary tract with reference to oestrogen status. BJOG. 2001;108:813–6.

    CAS  PubMed  Google Scholar 

  27. Shen SS, Smith CL, Hsieh JT, Yu J, Kim IY, Jian W, Sonpavde G, Ayala GE, Younes M, Lerner SP. Expression of estrogen receptors-alpha and -beta in bladder cancer cell lines and human bladder tumor tissue. Cancer. 2006;106:2610–6.

    Article  CAS  PubMed  Google Scholar 

  28. Batra SC, Iosif CS. Progesterone receptors in the female lower urinary tract. J Urol. 1987;138:1301–4.

    CAS  PubMed  Google Scholar 

  29. Håkonsen LB, Ernst A, Ramlau-Hansen CH. Maternal cigarette smoking during pregnancy and reproductive health in children: a review of epidemiological studies. Asian J Androl. 2014;16:39–49.

    Article  PubMed  Google Scholar 

  30. Chaliha C, Stanton SL. Urological problems in pregnancy. BJU Int. 2002;89:469–76.

    Article  CAS  PubMed  Google Scholar 

  31. Zhang W, Xiang YB, Fang RR, Cheng JR, Yuan JM, Gao YT. Total fluid intake, urination frequency and risk of bladder cancer: a population-based case-control study in urban Shanghai. Zhonghua Liu Xing Bing Xue Za Zhi. 2010;31:1120–4.

    PubMed  Google Scholar 

  32. Silverman DT, Alguacil J, Rothman N, Real FX, Garcia-Closas M, Cantor KP, Malats N, Tardon A, Serra C, Garcia-Closas R, Carrato A, Lloreta J, Samanic C, Dosemeci M, Kogevinas M. Does increased urination frequency protect against bladder cancer? Int J Cancer. 2008;123:1644–8.

    Article  CAS  PubMed  Google Scholar 

Download references


Thanks to the help of the members of the urology department in West China Hospital of Sichuan University, especially QW, CXP and JHL.


This work was collectively supported by grant (National Natural Science Foundation of China (No. 81270841)), and Science & Technology Pillar Program from Science & Technology Department of Sichuan Province (2013SZ0034). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Availability of data and materials

All data generated or analysed during this study are included in this published article.

Authors’ contributions

Conceived and designed the experiments: PH. Analyzed the data: JYB and XMW. Contributed reagents/materials/analysis YBY, YT and JW. Wrote the manuscript: JYB, XMW and JW. All authors have read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Ping Han.

Additional file

Additional file 1:

Subgroup analyses of forest plots included in our meta-analysis. (ZIP 143 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Bai, Y., Wang, X., Yang, Y. et al. Parity and bladder cancer risk: a dose-response meta-analysis. BMC Cancer 17, 31 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: