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

Discrete survival model analysis of a couple’s smoking pattern and outcomes of assisted reproduction



Cigarette smoking has been associated with worse infertility treatment outcomes, yet some studies have found null or inconsistent results.


We followed 225 couples who underwent 354 fresh non-donor assisted reproductive technology (ART) cycles between 2006 and 2014. Smoking history was self-reported at study entry. We evaluated the associations between smoking patterns and ART success using multivariable discrete time Cox proportional hazards models with six time periods: cycle initiation to egg retrieval, retrieval to fertilization, fertilization to embryo transfer (ET), ET to implantation, implantation to clinical pregnancy, and clinical pregnancy to live birth to estimate hazard ratios (HR) and 95% CIs. Time-dependent interactions between smoking intensity and ART time period were used to identify vulnerable periods.


Overall, 26% of women and 32% of men reported ever smoking. The HR of failing in the ART cycle without attaining live birth for male and female ever smokers was elevated, but non-significant, compared to never smokers regardless of intensity (HR = 1.02 and 1.30, respectively). Female ever smokers were more likely to fail prior to oocyte retrieval (HR: 3.37; 95% CI: 1.00, 12.73). Every one cigarette/day increase in smoking intensity for females was associated with a HR of 1.02 of failing ART (95% CI: 0.97, 1.08), regardless of duration or current smoking status. Women with higher smoking intensities were most likely to fail a cycle prior to oocyte retrieval (HR: 1.07; 95% CI: 1.00, 1.16). Among past smokers, every additional year since a man had quit smoking reduced the risk of failing ART by 4% (HR: 0.96; 95% CI: 0.91, 1.00) particularly between clinical pregnancy and live birth (HR: 0.86; 95% CI: 0.76, 0.96).


Female smoking intensity, regardless of current smoking status, is positively associated with the risk of failing ART cycles between initiation and oocyte retrieval. In men who ever smoked, smoking cessation may reduce the probability of failing ART, particularly between clinical pregnancy and live birth.

Trial registration

NCT00011713. Registered: 27 February 2001.


Tobacco is among the leading causes of death and disability worldwide [1], with over 6 million people dying of tobacco-related causes per year [2]. It is a common exposure in the US, with smoking reported for approximately 25% of reproductive-aged women, 12% of women during pregnancy [3], and 25% of men aged 25 to 44 years [4]. Cigarette smoke contains about 4000 compounds that have been associated with increased risk of cardiovascular disease, lung disease, cancer [5, 6] and adverse reproductive outcomes [7, 8].

Previous studies have shown that a couple’s likelihood of success when undergoing assisted reproduction varies depending on their smoking history [9]. Female smoking has been associated with poorer ovarian response markers, lower oocyte counts, and lower rates of implantation [10] and live birth [11, 12] following assisted reproductive technology (ART). Male smoking has also been associated with lower live birth rates [12], independent of female smoking. To date, most studies have been limited in their assessment of smoking history, generally comparing “smokers” to “non-smokers” without making clear distinctions among past and never smokers and rarely taking into account smoking intensity, duration or cessation. Furthermore none of the previous studies have utilized methods which allow distinction of smoking effects by stage of the ART cycle, and thus have not always been able to account for early failures, which may have resulted in contradictory findings [13].

Thus it remains unclear whether characteristics of a couple’s smoking history, such as duration and intensity of smoking, the age at smoking initiation, and years since smoking cessation could affect the probability of succeeding at ART. To address this, we evaluated the association of pre-treatment smoking patterns and infertility treatment outcomes in a prospective cohort of couples undergoing ART at the Massachusetts General Hospital (MGH) Fertility Center, using novel statistical methods to allow for a changing effect of smoking on ART failure across time periods within an ART cycle.


Study overview

Participants were couples enrolled in the Environment And Reproductive Health (EARTH) study, an ongoing prospective cohort study started in 2006 aimed at identifying determinants of fertility among couples presenting to the MGH Fertility Center (Boston, MA) [14]. The study was approved by the Institutional Review Boards of MGH and the Harvard T. H. Chan School of Public Health. All participants provided written informed consent.

Study population

Couples presenting to the MGH Fertility Center were given information about the EARTH Study by the medical team. Study personnel approached couples who met the age requirements (18–46 years for women; 18–55 years for men), and who initially planned to use their own gametes for treatment. Approximately 60% of those contacted participated in the study. Participants were included in the main analysis if they underwent at least one fresh non-donor ART cycle between 2006 and 2014 (n = 392). Due to our interest in smoking patterns among couples, we excluded 156 women whose partner had not yet enrolled or elected not to enroll in the EARTH study. Some aspects differed in the excluded couples compared to those included. On average, the female partners from the excluded couples were slightly older (36 vs. 35 years), had higher BMIs (25 vs.24 kg/m2), had more pack-years of smoking (1.4 vs. 0.7 pack-years), and were less likely to have completed at least a college education (78% vs. 86%). The cycles from the excluded couples were also more likely to fail prior to embryo transfer (15% vs. 8%). After exclusions, 236 couples who underwent a total of 354 fresh non-donors ART cycles were included in the analysis.

Smoking assessment

Participants were asked to report whether they had ever smoked cigarettes (defined as smoking at least 1 cigarette/day for a year or at least 20 packs in their lifetime) in a brief staff-administered questionnaire. Individuals reporting a positive smoking history were asked about their age at initiation, current smoking status, smoke inhalation practices, and the brand of cigarette they used (including size, tobacco type, and filter). Couples reporting current use reported what day they had their last cigarette, how many cigarettes they currently smoke per day, average number of cigarettes smoked per day over their years as a current smoker, and whether they had ever quit for 6 months or more and, if yes, how many years they quit. Former users reported how old they were they quit smoking, if they quit in the last year, the date, quantity of cigarettes smoked per day during the period they smoked, and if before they quit there were periods (≥6 months) when they didn’t smoke. The total smoking duration was calculated in current users by subtracting the age of smoking initiation from their current age. In former smokers we subtracted the age of smoking initiation from their age on cessation. For both we subtracted their interim cessation periods from their total smoking duration. Pack-years of smoking was calculated as: \( \frac{Average\kern0.24em number\; of\; cigarettes\; smoked\; a\; day}{20\; cigarettes\; per\; pack}\times Total\; smoking\; duration(years) \) [15]. Questions on the use of other tobacco products (i.e. cigars, pipes, chewing tobacco and snuff or dip tobacco) were included on the questionnaire, but due to their infrequent use in our sample (<2% for cigars and <1% for the rest), only cigarettes smoking was considered as an exposure in our analyses.


At enrollment, a brief nurse-administered questionnaire was used to collect data on demographic characteristics, medical history, and lifestyle factors. A trained research study staff member measured each participant’s height and weight.

ART Outcomes

Women underwent a pretreatment cycle of oral contraceptives for 2–5 weeks to suppress ovulation before their ART cycles, unless contraindicated. On day 3 of induced menses, patients began controlled ovarian stimulation. Women underwent one of three stimulation protocols as clinically indicated: 1) luteal-phase GnRH agonist protocol; 2) follicular-phase GnRH-agonist/flare protocol; or 3) GnRH-antagonist protocol. Patients were monitored during gonadotropin stimulation for serum estradiol and endometrial thickness through 2 days before egg retrieval. Human chorionic gonadotrophin was administered 36 h before the scheduled egg retrieval procedure to induce ovulation. Details of oocyte retrieval have been previously described [16]. Embryologists classified oocytes as germinal vesicle, metaphase I, metaphase II (MII), or degenerated. Following retrieval, couples underwent conventional insemination or intracytoplasmic sperm injection (ICSI) for fertilization. We defined successful implantation as a serum β-hCG level >6 mIU/mL 15–20 days after egg retrieval, clinical pregnancy as the presence of an intrauterine pregnancy in ultrasound, and live birth as the birth of a neonate at or after 24 weeks of gestation. All clinical information was abstracted from electronic medical records.

Statistical analysis

Both male and female participants were classified into three categories: never, former, and current smokers according to self-reported smoking status. We then cross-classified couples into four categories according to their ever smoking status. We evaluated the effect of smoking intensity using the variables average cigarettes/day and total years of smoking as continuous exposures. Descriptive statistics were calculated for demographic and reproductive characteristics according to current smoking status prior to ART treatment. A Chi-square test or Fisher’s exact test was used to test for differences across categories for discrete variables and the Kruskal-Wallis test for differences across categories for continuous variables.

Multivariable Cox proportional hazards models for discrete survival time were used to estimate the hazard ratios (HRs) and 95% confidence intervals (CIs) for the associations between smoking patterns and risk of failing ART prior to live birth, as proposed by Maity [17] and Missmer [18]. In all of these models, a robust sandwich covariance estimate was used to account for the multiple cycles per woman. HRs estimate the odds of failing an ART cycle at any point, conditional on not failing at an earlier moment during the same cycle. Women were considered at risk of failing ART for the duration of their initiated cycle until their specific point of failure. There were six periods during which women could fail a cycle: between (1) cycle initiation and oocyte retrieval, (2) retrieval and fertilization, (3) fertilization and embryo transfer (ET), (4) ET and implantation, (5) implantation and clinical pregnancy, and (6) between clinical pregnancy and live birth. Estimated probabilities of survival for the average women in our cohort (according to mean or reference levels for covariates other than smoking) were used to construct survival curves across these time points.

Confounding was evaluated through directed acyclic graphs based on prior knowledge and descriptive statistics in the study population. The multivariable models were adjusted for age, BMI, protocol, educational level and partner’s smoking history (never or ever smoker). We centered smoking intensity, years of smoking duration and pack-years by the whole sample’s mean for continuous analyses. For analysis of smoking history (ever vs. never smoker), smoking status (never, former and current smoker) and couple’s cross-classifications, we further adjusted by centered pack-years. We analyzed as continuous variables smoking intensity and duration in the entire population years since smoking cessation and age at smoking initiation among ever smokers (58 women and 71 men). Never smokers were assigned a value of 0 average cigarettes/day and 0 years of total smoking for the continuous analyses. Since in Massachusetts insurance companies require current smokers to have cotinine levels below the non-smoker reference value for 2 months to provide authorization for IVF coverage, we assigned a cessation time of 2 months for current smokers. All results are presented as HRs for a 1 unit increase in the continuous exposure. We further adjusted the models for smoking intensity (cigarettes/day) and smoking duration (years) using a binary variable for self-reported current smoking status. Models for smoking cessation and age at initiation were further adjusted for pack-years. The model of age at smoking initiation was also further adjusted for a binary variable for self-reported current smoking-status. To evaluate whether there were specific stages of the cycle which were more vulnerable to the effects of smoking, we fit multivariate models with time-dependent interaction between the discrete time periods and smoking exposures.

We also ran a variety of sensitivity analysis to address potential issues related to misclassification of the exposure, selection bias into the study, and model fitting. In accordance with a prior methodological paper on smoking and lung cancer [19], we conducted a sensitivity analysis re-classifying the participant’s self-reported smoking status according to their time of cessation in three discrete categories: never smokers (self-reported never smokers), former smokers (self-reported former smokers with > 2 years of smoking cessation) and “clinically active smokers” (self-reported current smokers and self-reported former smokers with ≤ 2 years of smoking cessation). Also, in order to assess the possibility of selection bias in our sample (due to the exclusion of women who did not have a partner enrolled), we conducted a sensitivity analysis including all woman regardless of whether their partner was enrolled and analyzed the association between smoking history and probability of failing ART without adjusting for the male partner’s smoking history. Finally, we performed an analysis limited to the couple’s first cycle to address the potential limitation of this model to account for the imbalanced number of cycles per couple.

To illustrate the possible effect of survival bias, we used generalized linear mixed models with random intercepts to analyze the associations between smoking patterns and oocyte yield (Poisson distribution and log link); endometrial wall thickness and peak serum estradiol at trigger day (binomial distribution and log link) among women who underwent oocyte retrieval, and probability of live birth per initiated cycle and per embryo transfer (binary distribution and logit link). These results are presented as population marginal means adjusted for covariates. All analyses were conducted using the Statistical Analysis System Software package SAS 9.4 (SAS Institute Inc., Cary, NC) and considered two-sided significance levels < 0.05 as statistically significant.


Our sample consisted mostly of non-Hispanic Caucasians (86% of women; 88% of men) with a mean age of 34.8 and 36.6 years for female and male participants, respectively. About 26% of women were self-reported ever smokers (23% former and 3% current smokers); 32% of men were ever smokers (24% former and 7% current smokers). On average, women who ever smoked started smoking when they were 17 years, smoked 6 cigarettes/day for 9 years, and if they had quit, they did so 7.5 years before cycle initiation. Similarly, men who ever smoked started smoking when they were an average of 18 years old, and on average smoked 9 cigarettes per day for 12 years, and if they had quit, they did so 6 years before cycle initiation. Men and women who were current smokers were more likely to have a partner who currently smoked and were less likely to have a college education compared to never smokers (Tables 1 and 2). Current smokers also reported smoking for longer durations compared to former smokers. Current male smoker were more likely to have undergone reproductive tract surgery prior cycle initiation. Current female smokers were more likely to have their first cycle cancelled prior to embryo transfer and had higher baseline FSH levels.

Table 1 Baseline characteristics of 225 couples from the EARTH study in relation to female smoking status
Table 2 Baseline characteristics of 225 couples from the EARTH study in relation to male smoking status

Smoking status of the female partner, the male partner, and the cross-classified couple smoking status were unrelated to the risk of cycle failure (Table 3). The HR of failing ART for female ever smokers compared to never smokers was 1.30 (95% CI: 0.84, 2.02; P-Value =0.24) and estimated HRs remained similar even after accounting for smoking intensity and duration. When evaluating the effects of smoking during specific periods of potential ART failure, ever female smokers were marginally more likely to fail prior to oocyte retrieval (HR: 3.57; 95% CI: 1.00, 12.73; P-Value = 0.05) but not during other periods of the ART cycle. In the self-reported current smoking status analysis, the HR of failing ART for former and current female smokers was 1.23 (95% CI: 0.79, 1.91; P-Value = 0.35) and 2.41 (95% CI: 0.62, 9.39; P-Value = 0.20) when compared to never smokers. Similarly, current smokers (HR: 6.0; 95% CI: 0.96, 37.07; P-Value = 0.05) and former smokers (HR: 2.05; 95% CI: 0.50, 8.34; P-Value = 0.32) were more likely to fail prior to oocyte retrieval compared to never smokers but this association was not significant for the other ART time points and only borderline significant for current smokers. Women whose male partners were former and current smokers had no increased risk of failing ART when compared to never smokers. In addition, couple classifications in which either or both of the partners ever smoked, as compared to couples in which neither partner ever smoked, did not have an increase in risk of cycle failure.

Table 3 Estimated cumulative probabilities of live birth and corresponding hazard ratios of failing at ART by discrete smoking patterns (EARTH study, N = 225 couples, 354 ART cycles)

For the intensity and duration analyses (Table 4), a 1 cigarette/day increase in female smoking intensity was associated with a 7% increased risk of failing at ART (HR = 1.07, 95% CI: 0.97, 1.08; P-Value = 0.45) after further adjustment for smoking duration and current smoking status; while a 1 year increase in smoking duration yielded a HR of 1.01 (95% CI: 0.97, 1.05 P-Value = 0.61) after further adjustment for smoking intensity and current smoking status. Greater intensity of female smoking was associated with a borderline higher risk of failing a cycle prior to oocyte retrieval (HR: 1.07; 95% CI: 1.00, 1.16; P-Value = 0.07); In contrast smoking duration was not significantly related to the probability of failing ART during any period of the ART cycle. Neither time of cessation nor age of smoking initiation was significantly associated with the risk of failing ART (Table 4). We found a borderline significant protective effect of smoking cessation in ever-smoker males; for every additional year since a male smoker quit, the risk of failing during an ART cycle decreased by 4% (HR: 0.96; 95% CI: 0.91, 1.00; P-Value <0.05) regardless of smoking duration and intensity. This association was particularly pronounced for the time period between clinical pregnancy and live birth (HR: 0.86; 95% CI: 0.76, 0.96; P-Value = 0.01).

Table 4 Hazard ratios for failing ART by continuous smoking patterns (EARTH study, N = 225 couples, 354 ART cycles)

Reclassification of smoking status by taking into account time since smoking cessation did not significantly alter our results. For instance, compared to never smokers, former smokers had a 23% increase in risk (HR = 1.23,95% CI: 0.78, 1.95; P-Value = 0.37) while “clinically active” smokers had a HR of 1.62 (95% CI: 0.70, 3.78; P-Value = 0.37) of failing at ART; similar to the main analysis, the increased risk of failing at ART was predominantly prior to oocyte retrieval (HR: 2.47; 95% CI: 1.27, 4.78; P-Value = 0.01) for “clinically active” smokers. The sensitivity analysis including the women without a male partner enrolled in the study provided consistent results to those seem in the main sample. Specifically, female ever smokers had a borderline significant HR of 1.38 (95% CI: 0.99, 1.92; P-Value = 0.06) compared to never smokers; also, this increased risk was significantly higher for failing prior to oocyte retrieval (HR: 2.45; 95% CI: 1.12, 5.36; P-Value = 0.02) and it was also significant between embryo transfer and implantation (HR: 1.63; 95% CI: 1.03, 2.58; P-Value = 0.04). An increase of 1 cigarette/day produced an overall HR of 1.01 (95% CI: 0.97, 1.06; P-Value = 0.65) for failing at ART and the HR was not significantly different according to the particular moment of the cycle. Finally, sensitivity analyses limited to the couple’s first cycle produced results that were consistent with the main analysis, although with loss of significance given the reduction in power: female ever smokers had an HR of 1.34 (95% CI: 0.76, 2.34; P-Value = 0.31) of failing ART compared to never smokers and this increased risk was more pronounced prior to oocyte retrieval (HR: 3.00; 95% CI: 0.56, 16.08; P-Value = 0.20) although these increases were not statistically significant.

Given the observed relation between female smoking and early ART failures, we analyzed the associations between smoking patterns and oocytes counts among women with successful egg retrieval using conventional methods to illustrate the possible effect of survival bias. In these models, self-reported current smokers appeared to have a higher number of total and MII phased oocytes retrieved than former and never smokers even after multivariable adjustment, although these differences did not attain statistical significance (Additional file 1: Table S1); however, 16.7% of cycles from self-reported current smokers could not be included in this analysis because they failed prior to retrieval compared to only 2.0% of cycles from never smokers and 6.5% from former smokers (Additional file 1: Table S1). We did not observe associations with endometrial thickness or peak estradiol levels. When we analyzed the associations between smoking patterns and live birth following ART among the cycles of women undergoing an embryo transfer and all initiated cycles, the predicted probability of live birth was the lowest for current smokers in both cases. However, the differences by smoking status for women undergoing an embryo transfer were attenuated compared to those observed taking into consideration all initiated cycles (Additional file 2: Table S2). This was to be expected as 25% of cycles from self-reported current smokers failed prior to ET compared to only 6.8% of cycles from never smokers.


In a prospective cohort of couples undergoing infertility treatment with ART, we found that the female smoking intensity was related to a higher risk of cycle failure prior to oocyte retrieval, even after accounting for duration of smoking and current smoking status. Our results suggest that the effects of smoking on ovarian response to controlled ovarian stimulation may persist after smoking cessation. We also found a borderline protective effect for male smoking cessation and cycle failure between clinical pregnancy and live birth in former smokers. Our conventional analysis highlighted the possibility that studies ignoring failures prior to oocyte retrieval for intermediate outcomes and failures prior to embryo transfer for clinical outcomes may be susceptible to survival bias, possibly explaining some inconsistencies in the literature of the effects of smoking on human fertility, such as a higher number of oocytes or no differences in clinical outcomes comparing current to never smokers [13].

Despite the low prevalence of current smoking in our cohort, our results do suggest that female smoking is detrimental to overall ART success. These results are in accordance with previous cohort studies and meta-analyses which reported lower probabilities of implantation, clinical pregnancy and live birth following ART in smokers compared to never smokers [9, 12, 20, 21]. Similar to our findings, many of these previous studies also found that female smokers were more likely to fail before an embryo transfer [12, 21, 22]. The higher odds of failing between cycle initiation and oocyte retrieval could imply that cigarette smoke disrupts the response to controlled ovarian hyper stimulation protocols or disrupts the ovary itself. This mechanism could be related to findings in animal models, in which cigarette smoke has been found to alter the serum levels of sex hormones and gametogenesis [23, 24]. Although this mechanism accounts only for contemporaneous exposures there seems to be a relation between lower fecundity and cumulative smoking exposure in epidemiological studies [25]. Other epidemiological studies have shown that current smokers were more likely to have worse ovarian stimulation outcomes compared to never smokers [26, 27]. Our findings also suggest an increased risk of failure at implantation that may have been underpowered in our study, as further evidenced by the results of the second sensitivity analysis which included a larger sample of women.

Not all epidemiological studies have found negative associations between smoking and ART success [13, 28]. Although differences in study design and power may have influenced these results, these studies also tended to analyze clinical outcomes only among cycles reaching the point of embryo transfer. Due to the effects of smoking on early ART outcomes (e.g. response to controlled ovarian stimulation), restricting analyses to only cycles with embryo transfer creates a selection bias in which only the women who survived until then are included, biasing the observed burden of smoking downwards [29] even despite the well-recognized associations between smoking and miscarriage, stillbirth and neonatal death [30]. When we analyzed the association between smoking patterns and live birth per embryo transfer (as opposed to per initiated cycle; Additional file 2: Table S2), we illustrated how this issue may be at play. We also found paradoxical results when we evaluated the association between smoking patterns and oocyte counts among women with successful egg retrieval, as current smokers appeared to yield as many oocytes compared to former and never smokers. However, this association appeared to be largely driven by the 12% of cycles from smokers that failed prior to egg retrieval which was more than double the percentage of failures from never (2.0%) and former smokers (6.3%). Future observational studies should use the appropriate statistics to account for these effects, particularly when exposures may exert effects early in gestation.

Other than a possible beneficial effect of smoking cessation among male ever-smokers, we did not observe any significant associations when studying male smoking patterns and ART outcomes. Nevertheless, male smoking has been associated with lower semen quality parameters [31] and live birth rates [22] in other epidemiological studies. While our analyses were restricted by power, the fact that former male smokers seemed to more quickly benefit from smoking cessation could be explained by the biology of sperm division and differentiation. The mechanism for the observed reduction in cycle failure between clinical pregnancy and live birth due to paternal smoking cessation could be due to lower maternal exposure to second hand smoking and/or a direct paternal exposure [32], although the exact mechanisms remain obscure.

Our study has certain limitations. Due to our design, it may not be possible to extrapolate our results to the general population. However, our cohort of men and women are similar to the population of infertility patients nationwide [33] suggesting that our results might be representative of other couples undergoing ART. While we lacked specific information on the women who declined to participate in the EARTH study, we had a similar smoking prevalence as other ART cohorts in the US [21]. Furthermore, in a previous study on semen quality from MGH, we found that men who refused to participate were comparable to the men who did participate in terms of baseline characteristics [34]. This suggests that the likelihood of selection bias in our study (e.g. smoking men and women being less likely to participate) is low. The women excluded had different baseline characteristics than the included sample; however, our results were similar after performing a sensitivity analysis among all EARTH women and, in fact, some associations gained significance with this increase in power. Another generalizability concern is that in Massachusetts, insurance companies will generally not allow current smokers to have authorization for IVF coverage until cotinine levels are below the non-smoker reference value for 2 months. Therefore self-reported current smokers might not have been smoking for the past two months, which might have attenuated the magnitude of associations towards the null. To address this we assumed a conservative cessation time of at least 2 months for current smokers when modeling time of smoking cessation a continuous exposure. Smoking history was assessed through self-report in our population as opposed to more objective markers like cotinine. While it is possible there was some misclassification of current smokers due to its subjectivity, a person’s self-report of smoking is considered reliable [35]. To minimize the possible effect of passive smoking, we adjusted for the partner’s smoking history. Although our survival models can account for multiple cycles per women, they are unable to account for the imbalance in number of cycles per couple (e.g. women with more cycles are more likely to fail) [17], yet even when we restricted the analysis to just the first cycle per couple the results were similar. Residual and unmeasured confounding is always possible as this was an observational cohort study; however, we attempted to minimize this by adjusting for number of demographic and lifestyle factors and focusing on a couple analyses. Our main limitation was the small numbers of former and current smokers, reducing our power to detect associations. Also, the subjects who did smoke reported being mostly “light” or “moderate” smokers in clinical terms. Our study strengths include our prospective design, our detailed assessment of smoking history, and our ability to account for many confounders including the partner’s smoking.


Female smoking intensity was positively associated with the risk of failing in ART cycles between cycle initiation and oocyte retrieval. Among former male smokers, every additional year of smoking cessation was associated with a slightly decreased risk of failing at ART between clinical pregnancy and live birth. Future studies should further investigate the long-term effects of smoking intensity on reproductive potential and take care to use the appropriate analytical approaches to account for failures at early stages of reproduction to reduce the likelihood of survivor bias.


95% CI:

95% confidence interval


Assisted reproductive technology


Body mass index


Environment And Reproductive Health study


Embryo transfer


Hazard ratio


Intracytoplasmic sperm injection


In-vitro fertilization


Massachusetts General Hospital


Metaphase II oocytes


  1. Havelaar AH, Kirk MD, Torgerson PR, Gibb HJ, Hald T, Lake RJ, Praet N, Bellinger DC, de Silva NR, Gargouri N, et al. World Health Organization Global Estimates and Regional Comparisons of the Burden of Foodborne Disease in 2010. PLoS Med. 2015;12:e1001923.

    Article  PubMed  PubMed Central  Google Scholar 

  2. World Health Organization (WHO). WHO report on the global tobacco epidemic, 2015: Raising taxes on tobacco. In: Tobacco Free Initiative. Geneva: WHO; 2015.

    Google Scholar 

  3. Tong VT, Dietz PM, Morrow B, D’Angelo DV, Farr SL, Rockhill KM, England LJ. Trends in smoking before, during, and after pregnancy--Pregnancy Risk Assessment Monitoring System, United States, 40 sites, 2000–2010. MMWR Surveill Summ. 2013;62:1–19.

    PubMed  Google Scholar 

  4. Agaku IT, King BA, Husten CG, Bunnell R, Ambrose BK, Hu SS, Holder-Hayes E, Day HR. Tobacco product use among adults--United States, 2012–2013. MMWR Morb Mortal Wkly Rep. 2014;63:542–7.

    PubMed  Google Scholar 

  5. Kaiserman MJ, Rickert WS. Carcinogens in tobacco smoke: benzo[a]pyrene from Canadian cigarettes and cigarette tobacco. Am J Public Health. 1992;82:1023–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Ding YS, Trommel JS, Yan XJ, Ashley D, Watson CH. Determination of 14 polycyclic aromatic hydrocarbons in mainstream smoke from domestic cigarettes. Environ Sci Technol. 2005;39:471–8.

    Article  CAS  PubMed  Google Scholar 

  7. Augood C, Duckitt K, Templeton AA. Smoking and female infertility: a systematic review and meta-analysis. Hum Reprod. 1998;13:1532–9.

    Article  CAS  PubMed  Google Scholar 

  8. Jeng HA, Chen YL, Kantaria KN. Association of cigarette smoking with reproductive hormone levels and semen quality in healthy adult men in Taiwan. J Environ Sci Health A Tox Hazard Subst Environ Eng. 2014;49:262–8.

    Article  CAS  PubMed  Google Scholar 

  9. Waylen AL, Metwally M, Jones GL, Wilkinson AJ, Ledger WL. Effects of cigarette smoking upon clinical outcomes of assisted reproduction: a meta-analysis. Hum Reprod Update. 2009;15:31–44.

    Article  CAS  PubMed  Google Scholar 

  10. Winter E, Wang J, Davies MJ, Norman R. Early pregnancy loss following assisted reproductive technology treatment. Hum Reprod. 2002;17:3220–3.

    Article  PubMed  Google Scholar 

  11. Penzias AS. Recurrent IVF failure: other factors. Fertil Steril. 2012;97:1033–8.

    Article  PubMed  Google Scholar 

  12. Klonoff-Cohen H, Natarajan L, Marrs R, Yee B. Effects of female and male smoking on success rates of IVF and gamete intra-Fallopian transfer. Hum Reprod. 2001;16:1382–90.

    Article  CAS  PubMed  Google Scholar 

  13. Cinar O, Dilbaz S, Terzioglu F, Karahalil B, Yucel C, Turk R, Taskin L, Kose SK. Does cigarette smoking really have detrimental effects on outcomes of IVF? Eur J Obstet Gynecol Reprod Biol. 2014;174:106–10.

    Article  CAS  PubMed  Google Scholar 

  14. Hauser R, Meeker JD, Duty S, Silva MJ, Calafat AM. Altered semen quality in relation to urinary concentrations of phthalate monoester and oxidative metabolites. Epidemiology. 2006;17:682–91.

    Article  PubMed  Google Scholar 

  15. Pinsky PF, Zhu CS, Kramer BS. Lung cancer risk by years since quitting in 30+ pack year smokers. J Med Screen. 2015;22:151–7.

    Article  PubMed  Google Scholar 

  16. Mok-Lin E, Ehrlich S, Williams PL, Petrozza J, Wright DL, Calafat AM, Ye X, Hauser R. Urinary bisphenol A concentrations and ovarian response among women undergoing IVF. Int J Androl. 2010;33:385–93.

    Article  CAS  PubMed  Google Scholar 

  17. Maity A, Williams PL, Ryan L, Missmer SA, Coull BA, Hauser R. Analysis of in vitro fertilization data with multiple outcomes using discrete time-to-event analysis. Stat Med. 2014;33:1738–49.

    Article  PubMed  Google Scholar 

  18. Missmer SA, Pearson KR, Ryan LM, Meeker JD, Cramer DW, Hauser R. Analysis of multiple-cycle data from couples undergoing in vitro fertilization: methodologic issues and statistical approaches. Epidemiology. 2011;22:497–504.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Leffondre K, Abrahamowicz M, Siemiatycki J, Rachet B. Modeling smoking history: a comparison of different approaches. Am J Epidemiol. 2002;156:813–23.

    Article  PubMed  Google Scholar 

  20. Neal MS, Hughes EG, Holloway AC, Foster WG. Sidestream smoking is equally as damaging as mainstream smoking on IVF outcomes. Hum Reprod. 2005;20:2531–5.

    Article  PubMed  Google Scholar 

  21. Wright KP, Trimarchi JR, Allsworth J, Keefe D. The effect of female tobacco smoking on IVF outcomes. Hum Reprod. 2006;21:2930–4.

    Article  CAS  PubMed  Google Scholar 

  22. Fuentes A, Munoz A, Barnhart K, Arguello B, Diaz M, Pommer R. Recent cigarette smoking and assisted reproductive technologies outcome. Fertil Steril. 2010;93:89–95.

    Article  CAS  PubMed  Google Scholar 

  23. Mai Z, Lei M, Yu B, Du H, Liu J. The effects of cigarette smoke extract on ovulation, oocyte morphology and ovarian gene expression in mice. PLoS One. 2014;9:e95945.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Sadeu JC, Foster WG. Cigarette smoke condensate exposure delays follicular development and function in a stage-dependent manner. Fertil Steril. 2011;95:2410–7.

    Article  CAS  PubMed  Google Scholar 

  25. Bolumar F, Olsen J, Boldsen J. Smoking reduces fecundity: a European multicenter study on infertility and subfecundity. The European Study Group on Infertility and Subfecundity. Am J Epidemiol. 1996;143:578–87.

    Article  CAS  PubMed  Google Scholar 

  26. Van Voorhis BJ, Syrop CH, Hammitt DG, Dunn MS, Snyder GD. Effects of smoking on ovulation induction for assisted reproductive techniques. Fertil Steril. 1992;58:981–5.

    Article  PubMed  Google Scholar 

  27. Freour T, Masson D, Dessolle L, Allaoua D, Dejoie T, Mirallie S, Jean M, Barriere P. Ovarian reserve and in vitro fertilization cycles outcome according to women smoking status and stimulation regimen. Arch Gynecol Obstet. 2012;285:1177–82.

    Article  PubMed  Google Scholar 

  28. Weigert M, Hofstetter G, Kaipl D, Gottlich H, Krischker U, Bichler K, Poehl M, Feichtinger W. The effect of smoking on oocyte quality and hormonal parameters of patients undergoing in vitro fertilization-embryo transfer. J Assist Reprod Genet. 1999;16:287–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Buckley JP, Keil AP, McGrath LJ, Edwards JK. Evolving methods for inference in the presence of healthy worker survivor bias. Epidemiology. 2015;26:204–12.

    Article  PubMed  Google Scholar 

  30. Hyland A, Piazza KM, Hovey KM, Ockene JK, Andrews CA, Rivard C, Wactawski-Wende J. Associations of lifetime active and passive smoking with spontaneous abortion, stillbirth and tubal ectopic pregnancy: a cross-sectional analysis of historical data from the Women’s Health Initiative. Tob Control. 2015;24:328–35.

    Article  PubMed  Google Scholar 

  31. Caserta D, Bordi G, Di Segni N, D’Ambrosio A, Mallozzi M, Moscarini M. The influence of cigarette smoking on a population of infertile men and women. Arch Gynecol Obstet. 2013;287:813–8.

    Article  PubMed  Google Scholar 

  32. Venners SA, Wang X, Chen C, Wang L, Chen D, Guang W, Huang A, Ryan L, O’Connor J, Lasley B, et al. Paternal smoking and pregnancy loss: a prospective study using a biomarker of pregnancy. Am J Epidemiol. 2004;159:993–1001.

    Article  PubMed  Google Scholar 

  33. Centers for Disease Control and Prevention (CDC), American Society for Reproductive Medicine (ASRM), Society for Assisted Reproductive Technology (SART). Assisted Reproductive Technology National Summary Report. Atlanta: US Dept of Health and Human Services; 2013.

    Google Scholar 

  34. Hauser R, Godfrey-Bailey L, Chen Z. Does the potential for selection bias in semen quality studies depend on study design? Experience from a study conducted within an infertility clinic. Hum Reprod. 2005;20:2579–83.

    Article  PubMed  Google Scholar 

  35. Soulakova JN, Hartman AM, Liu B, Willis GB, Augustine S. Reliability of adult self-reported smoking history: data from the tobacco use supplement to the current population survey 2002–2003 cohort. Nicotine Tob Res. 2012;14:952–60.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


We would like to thank all the patients and personnel involved in the EARTH study from the MGH Fertility Center and in particular thanks to Myra Keller, Ramace Dadd, Patricia Morey, and Tairmee Kangarloo, Harvard T.H. Chan School of Public Health research staff.


This work was supported by NIH grants: L50-HD085359 from NICHD, R01-ES009718 and P30-ES000002 from NIEHS.

Availability of data and materials

The datasets generated and/or analyzed during the current study are not publicly available because they contain parts of the study participant’s electronical medical records; but are available from the corresponding author on reasonable request.

Authors’ contributions

JCV: first author wrote manuscript and ran statistical analyses. JEC: Help conceived manuscript question and study design, provided manuscript review. PLW: Assisted in the creation of statistical models and manuscript review. JBF: Senior nurse in charge of coordinating participant’s enrollment, data entry and manuscript review. TLT: Senior obstetric fertility physician in charge of ART cycles and provided clinical insight and manuscript review. RH: Conceived initial study question and provided manuscript review. AJG: Senior and corresponding author, provided code and manuscript review as well a principal focus for the analysis. All authors read and approved the final manuscript.

Authors’ information

For citations and highest academic degree please use: Vanegas JC, MD; Chavarro JE, MD, ScD; Williams PL, PhD; Ford JE, RN, BSN; Toth TL, MD; Hauser R, MD, ScD, MPH; Gaskins AJ, ScD. Affiliations listed on title page.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

The study was approved by the Institutional Review Boards of MGH and the Harvard T. H. Chan School of Public Health. All participants provided written informed consent.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Audrey J. Gaskins.

Additional files

Additional file 1: Table S1.

Adjusted mean oocyte yield by smoking patterns among couples with successful egg retrieval. (EARTH Study, N = 225 couples, 354 ART cycles). (DOCX 17 kb)

Additional file 2: Table S2.

Adjusted estimated probability of cycles ending in live birth by smoking patterns (EARTH Study, N = 225 couples, 354 ART cycles). (DOCX 17 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

Vanegas, J.C., Chavarro, J.E., Williams, P.L. et al. Discrete survival model analysis of a couple’s smoking pattern and outcomes of assisted reproduction. Fertil Res and Pract 3, 5 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: