Seasonal variability does not impact in vitro fertilization success

It is unknown whether seasonal variation influences the outcome of in vitro fertilization (IVF). Previous studies related to seasonal variation of IVF were all small sample size, and the results were conflicting. We performed a retrospective cohort study evaluating the relationship between seasonal variability and live birth rate in the year of 2014–2017. Patients were grouped into four seasons (Winter (December-February), Spring (March-May), Summer (June-August), and Autumn (September-November)) according to the day of oocyte pick-up (OPU). Multivariate logistic regression analysis was performed to evaluate association between seasonal variation and live birth. Models were adjusted for covariates including temperature, sunshine hour, infertility type, infertility duration, infertility factor and BMI. In total 38,476 women were enrolled, of which 25,097 underwent fresh cycles, 13,379 were frozen embryo transfer. Live birth rates of fresh embryo transfer were 50.36%, 53.14%, 51.94% and 51.33% for spring, summer, autumn and winter, respectively. Clinical pregnancy rate between the calendar months varied between 55.1% and 63.4% in fresh embryo transfer (ET) and between 58.8% and 65.1% in frozen embryo transfer (FET) (P-values 0.073 and 0.220). In the unadjusted model and adjust model, seasonal variation was not associated with live birth. In conclusion, there was no significant difference of seasonal variations in the outcome of IVF with fresh embryo transfer and frozen embryo transfer.

mark of fertilization. Embryos were cultured in G5-medium (Vitrolife, Sweden). Embryo scoring was performed by the combination of blastomere number, size, and fragmentation 7 . Three to five days after oocyte retrieval 1 or 2 good quality embryos were transferred into the uterus.
Our IVF center is in Xi'an, Shaanxi province of China. It is dominated by monsoon climate of medium latitudes, with high temperature and rain in summer and cold dry climate in winter, which gives us a perfect model to study the relationship of seasonal variation and IVF outcome. Monthly weather conditions (°C) and sunshine hours (Hrs) were extracted from the website of Shaanxi Provincial Bureau of Statistics. Meteorological data of Shaanxi province (China) during the year of 2014-2017 were shown in Supplementary Table 1.
The inclusion criteria were women who completed embryo transfer either in fresh cycle or frozen cycle. The exclusion criteria were: (1). history of recurrent pregnancy loss, (2). uterine pathology, (3) cycles that were cancelled due to no embryo available or failure of embryo thaw survival.
Among the 25097 IVF/ICSI stimulation cycles, 13,223(52.69%) cycles reached fresh embryo transfer. The other cycles embryos were frozen due to high risk of ovarian hyperstimulation syndrome (OHSS) (n = 2969), serum progesterone level >5.72 nmol/L (n = 832), no embryo available (n = 119), PGT (n = 1188), abnormal endometrium (n = 1782), patients with upper respiratory tract infection (n = 232) or diarrhea (n = 118) and other conditions (n = 4634) on the day of ET. Intramuscular injection of 60 mg Progesterone or 600 mg of vaginal micronized progesterone per day was given from the day of embryo transfer. The calendar months of embryo transfer were grouped to seasons, and each season lasting 3 months: spring (March-May), summer (June-August), fall (September-November) and winter (December-February).
Definition of clinical outcomes. Clinical pregnancy was defined as one or more gestational sacs confirmed by ultrasound. Miscarriage was defined as spontaneous loss of a clinical pregnancy before 22 completed weeks of gestational age. Ectopic pregnancy was defined as gestational sac observed by ultrasound outside the uterine cavity. Live birth was defined as the number of deliveries that resulted in at least one live birth >24 weeks.
The descriptive data on participants' characteristics were summarized using the mean and standard deviation for continuous variables. Counts and proportions were used for the categorical variables. Chi-squared tests were performed to compare the categorical variables. For multivariate analyses, logistic regression was utilized. Crude odds ratios (OR) and adjusted ORs (AOR) with 95% confidence interval (CI) were calculated to assess the association between covariates and live birth. In the adjusted Model I: we adjusted for temperature and sunshine hour. In the adjusted Model II: we adjusted for temperature, sunshine hour, infertility type, infertility duration, infertility factor and BMI. All of the analyses were performed with SPSS Version 13.0 (Statistical Package for Social Science, Inc., Chicago, IL, USA). The level of significance was set at p < 0.05.

Results
We included 25097 ovarian stimulation and OPU cycles of which 13223 cycles resulted in fresh embryo transfer and 13379 frozen embryo transfer cycles. The comparison of the general information among seasonal groups is shown in Table 1. In summer, the mean age of the patients was youngest (31.10 ± 5.07). BMI (22.49 ± 3.27) was highest and duration of infertility (3.96 ± 3.03) was longest in spring. Infertility type did not differ among seasonal groups. Tubal factor accounts for the largest proportion in summer (41.52%). In summer, patients were allocated agonist protocol most (81.80%). Male age was oldest in spring (33.13 ± 5.82), while sperm concentration was highest in autumn (54.97 ± 274.44) with significant difference.
We harvested a mean number of 10.47 ± 6.66 oocytes (Supplementary Table 2). In February, the number of OPU reached highest while in December, it dropped to the lowest. However, clinical pregnancy rate did not differ among 12 months. Likewise, analysis of FET cycles did not show any significant differences in the clinical pregnancy rate. In March, live birth rate seems to be lower. Thus, we compared live birth rate from the year of 2014 to 2017 and found live birth rate did not differ between different months.
IVF outcome among seasonal groups were shown in Table 2. The number of oocytes retrieved per OPU cycle and number of embryos transferred did not significantly differ between seasons. Meanwhile, gonadotropin (Gn) dosage and duration were largest in autumn. In spring and winter, the number of top-quality embryos was largest among seasons. In summer, clinical pregnancy rate reached highest, however, no difference was found in live birth rate and birth weight. Miscarriage rate and ectopic pregnancy rate also did not differ among seasons. A logistic regression model was then used to assess the association between seasonal variation and live birth while adjusting for potential confounders presented in Table 3. In the unadjusted model, seasonal variation was not associated with live birth. After adjusting for covariates, the association partially disappear, which indicated seasonal variation may impact live birth rate in IVF. In addition, the effect of winter is independent of sunshine and temperature.

Discussion
The seasonal variability of fertility has been studied in animals. However, seasonal influence on outcome of IVF in human is still under discussion. The effect of seasonal change on natural conception mainly in two ways: weather condition and temperature affect male sperm quality 8 , sunshine exposure influence female's ovulation. The optimal spermatogenesis in man requires 2-6°C lower testicular temperature than body and any rise above physiological temperature of the testes has adverse effect on spermatogenesis 9,10 . IVF procedure exert an extensive control over temperature, moisture, light and other condition, and huge effort has been made to standardize protocols and laboratory procedures and minimize any external influence on IVF outcome. Nevertheless, the effect of seasonal variability on IVF has been reported by some articles 6,11 , while others have reached different conclusions 12 . The inconsistent results are attributed to different inclusion criteria, widely diverse environment with different temperature and sunlight exposure when studies were conducted. (2019) 9:17185 | https://doi.org/10.1038/s41598-019-53919-3 www.nature.com/scientificreports www.nature.com/scientificreports/ In this study, we analyzed the seasonal variation of different variables related to IVF outcome. We use 3-month seasonal periods to minimize the cycle uneven distribution. Our results showed no significant differences in the number of oocytes retrieved, live birth rate or birth weight. Our findings were consistent with some studies 13 . One hypothesis is that in summer, male sperm quality is affected by temperature. Our data confirmed this hypothesis showing that male sperm concentration was lowest in summer with significant difference. However, no difference was found in clinical pregnancy rate among months and live birth rate among seasons. One explanation is that sperm concentration could be compensated by utilizing a concentrated amount of motile sperm in IVF or changing to ICSI for fertilization. However, our data from logistic regression showed seasonal variation may affect live birth rate especially in winter compared with spring. One cross-sectional study in China found lower serum vitamin D levels were associated with clinical pregnancy and live birth rate following IVF 14 . Serum 25(OH)-D levels were lowest in spring compared with other seasons, which could in turn affect pregnancy outcomes. Studies have demonstrated the fact that seasonal variation in serum AMH correlated with changes in seasonal 25(OH)-D (being 18% lower in winter than in summer) 15 .
There are other methods for controlling for time collinearity. A time-series study could be performed to assess the association of seasonal variation and live birth rate and evaluating the degree of the effects on estimation. In addition, sensitivity analyses can be conducted to examine whether there were nonlinear association between months and live birth rate. This study showed there was a trend of decreasing clinical pregnancy in February and March, but it did not reach significant difference. One explanation of monthly fluctuations is the climate change with sunlight exposure, a direct melatonin or neurotransmitter effect on the end-organ 16,17 . One study from Finland suggested melatonin secretion is increased during the dark season (November-January) when the amount of daylight is about 3-4 h, while ovarian activity is decreased 12 . Data from animal studies indicate the length of the nocturnal melatonin pulse is of central importance in gonadal regulation 18 .
In conclusion, by the analysis of the largest sample size to date, we confirm there was no significant difference of seasonal variations in the outcome of IVF with fresh embryo transfer and frozen embryo transfer. A routine IVF treatment should not be changed by season.