Abstract
Background The 2017 and 2018 wildfire seasons in British Columbia (BC), Canada were unprecedented. Among all the pollutants in wildfire smoke, fine particulate matter (PM2.5) poses the most significant risk to human health. There is limited research on prenatal wildfire-related PM2.5 exposure and its impacts on infant health.
Objectives We assessed the effect of prenatal wildfire-related PM2.5 exposure on respiratory infections in the first year of life. We explored exposure windows based on prenatal respiratory tract development and gestational trimester.
Methods All live singleton infants in utero during the wildfire seasons (July to September) from 2016 to 2019 were identified using the BC Perinatal Data Registry (BCPDR). The Canadian Optimized Statistical Smoke Exposure Model (CanOSSEM) was used to estimate daily PM2.5 exposures throughout pregnancy and the first year of life. Infants who were diagnosed with otitis media or a lower respiratory infection by age 1 year were identified, including infections associated with dispensations of the antibiotic amoxicillin. The odds ratios (OR) for exposures during specific developmental windows and gestational trimesters were estimated with logistic regression, conditioned on epidemiologic week of birth to control for seasonal trends. Analyses were adjusted for multiple environmental and individual-level covariates.
Results There were 134,475 infants included in the cohort, of which 22,093 were diagnosed with at least one outcome of interest. For every 10 μg/m3 increase in wildfire-related PM2.5, the OR [95% confidence interval] for otitis media associated with the fourth window of eustachian tube development (weeks 19-28) was 1.08 [1.02, 1.16]. Results were similar for infections associated with amoxicillin dispensations and during the second trimester (weeks 14-27) of development. Similarly, the canalicular stage of lower respiratory tract development (weeks 18-27) and the second trimester were associated with increased odds of lower respiratory tract infections.
Conclusion Specific windows of prenatal exposure to wildfire-related PM2.5 during the second trimester were associated with increased odds of otitis media and lower respiratory infections by age 1 year.
Study questions Is prenatal exposure to wildfire-related PM2.5 associated with risk of respiratory infections in the first year of life? If so, what are the critical prenatal exposure windows?
What’s already known Few studies have focused on prenatal exposure to wildfire-related PM2.5 and respiratory outcomes in the very early stages of life, nor the critical windows for gestational exposures. Evidence suggests that wildfire-related PM2.5 may impact the foetus directly or indirectly, potentially leading to adverse consequences for foetal respiratory tract development and vulnerability of the respiratory system in early life.
What this study adds Prenatal exposure to wildfire-related PM2.5 during specific periods of respiratory tract development were associated with increased risk of otitis media and lower respiratory infections by age 1 year. Results were similar for all infections and those associated with antibiotic dispensations. This study suggests that wildfire-related PM2.5 exposure during specific developmental windows can affect respiratory health in early life. Public health practitioners and healthcare providers should work to protect pregnant people and their children from the detrimental effects of wildfire smoke.
Background
Exposure to air pollution during pregnancy has been linked to negative effects in early life, such as cardiovascular and metabolic outcomes, respiratory and allergic responses, and neurodevelopmental impairments 1. Air pollution can directly impact pregnancy through dispersion across the tissue barriers and permeation across the cellular membrane of the placenta2. Indirectly, oxidative stress and inflammatory reactions in mothers can reduce foetal nutrient and oxygen supply 2,3. Interactions between air pollution and the epigenome can lead to adverse birth outcomes, including decreased respiratory function, though the mechanisms are not completely understood3.
Wildfire smoke is a complex form of air pollution, for which PM2.5 is generally accepted to be the most health-relevant marker. Compared with other PM2.5 sources such as traffic and industry, wildfire-related PM2.5 is far more episodic and the magnitudes of the exposure can be very high. Studies have shown that PM2.5 from wildfires may be more harmful to human respiratory health than ambient PM2.5 from other sources4,5. Although some studies have examined the relationship between wildfire smoke exposure and respiratory outcomes in adults and school-age children6–8, very few have looked at the effects of prenatal exposure on respiratory outcomes in early life.
A recent study in Australia found that a small cohort of infants prenatally exposed to smoke from a coal mine had more respiratory symptoms and diagnoses of respiratory tract infections in early life than an unexposed group9. Research from a birth cohort in Poland suggested that exposure to wildfire-related PM2.5 during pregnancy may increase the frequency of wheezing attacks in the first two years of life10. However, neither of these studies have assessed the relationship between prenatal exposure to wildfire-related PM2.5 during specific developmental windows.
Acute respiratory infections are the leading cause of morbidity and mortality in early childhood, even in high-income countries11. According to one study conducted in Canada, respiratory infections can impose substantial strain on acute care infrastructure and the medical system, and the burden is expected to increase with population growth12. Respiratory infections are also the primary reason for antibiotic use in early life, which can have long-lasting impacts on the infant microbiome and childhood respiratory health13. The most prevalent infections involving antibiotic treatment are acute otitis media, bronchitis, and upper respiratory tract infections14.
Exposure to wildfire smoke during pregnancy, or certain windows of pregnancy, may increase the risk of infant respiratory infections in early life. It follows that limiting exposure to wildfire smoke during pregnancy could protect infant health. To address this important public health question, our study focusses on the critical stages of gestational development for parts of the upper and lower respiratory tract associated with the two most diagnosed respiratory infections: otitis media (i.e., middle ear infections) and lower respiratory tract infections.
Otitis media is an upper respiratory infection that frequently requires extensive paediatric care in the early stages of life. It is characterized by acute or chronic inflammation of the middle ear, resulting from Eustachian tube dysfunction15,16. Lower respiratory tract infections include bronchitis, bronchiolitis, and pneumonia, all of which can lead to serious complications. Most lower respiratory infections in children are caused by bacteria, but they can also be caused by viruses and fungi17. The complex and protracted process of lower respiratory tract development, commonly referred to as foetal lung development, begins during the embryonic stage. Damage and changes to the upper and lower respiratory tracts during foetal development may lead to higher risk of respiratory diseases, including acute and chronic infections 18–20.
Our study aimed to investigate the prenatal impact of wildfire-related PM2.5 on respiratory infections during the first year of life using a population-based case-control design. We focused on specific stages of prenatal respiratory tract development and gestational trimesters to examine their association with respiratory infections in early life. Evidence from this study will be used to articulate preventive public health strategies to protect pregnant people and their infants.
Methods
Study setting and population
This study is set in the province of British Columbia (BC), on the west coast of Canada. Like other parts of western North America, BC is heavily forested and prone to seasonal wildfires21, especially from July through September22. The total area of BC is almost 1 million square kilometres. The 2016 population was approximately 5 million residents, with most living in large urban centres in the southern coastal areas and the rest living in smaller urban, rural, or remote areas throughout the interior and the north (Figure 1A).
In this study we focussed on the wildfire seasons from 2016-2019 because of the contrasts between these four years. The 2017 and 2018 wildfire seasons were both record-setting, with more than 1.2 million hectares burned in 2017 and more than 1.3 million burned in 2018 22,23. In both summers, the smoke impacts were widespread and prolonged (Figure 1B). By contrast, the 2016 and 2019 wildfire seasons were well below the 10-average, with less than 0.13 million hectares burned in both years, and the smoke impacts were limited23. The study included all singleton births that were prenatally exposed to the peak wildfire seasons (July – September) of 2016 to 2019 during any stage of pregnancy. Given an average gestation period of 38 weeks, this includes most infants born in BC between June 2016 and February 2020.
Health data sources
All health data were accessed through the Population Data BC research platform. The population-based cohort was identified using the BC Perinatal Data Registry (BCPDR)24. The BCPDR compiles extensive data on mothers, foetuses, and newborns, covering almost all births in the province. Most infants enter the registry at approximately 15 weeks of gestation. The BCPDR provides comprehensive information on delivery, the postpartum period, and includes many maternal covariates, such as BMI and smoking status, which are described in more detail in the covariates section to follow. The BCPDR also includes infant birthdates and estimated dates of conception, which facilitates precise exposure assessment for critical developmental windows25. Finally, the BCPDR includes the residential forward sortation area (FSA) of the mother, which is the first three digits of the 6-digit postal code.
BC has a single-payer healthcare insurance program called the Medical Services Plan (MSP), which covers most people in the province26,27. The MSP client roster was used to create a weekly residential history for mothers during their pregnancy and for infants during their first year of life28. The residential histories were created using 6-digit postal codes. In urban areas, each 6-digit postal code typically covers one building or one city block, but in rural settings they can cover much larger areas. In all cases, we used the geographic coordinates of the single link indicator, which indicates the location of most of the population within the postal code.
We used the MSP outpatient physician billings file to identify all infants who were diagnosed with otitis media or a lower respiratory infection by age 1 year. All prescription medications dispensed by community pharmacies in BC are recorded in the PharmaNet database29. By linking the MSP billings with PharmaNet records for each infant, we were able to identify infections that were associated with a subsequent antibiotic dispensation, as described in further detail in the outcomes section to follow.
Exposure
The Canadian Optimized Statistical Smoke Exposure Model (CanOSSEM) is an empirical framework that estimates daily population exposure to PM2.5 from biomass burning across Canada from 2010 to 2022 at a resolution of 5 km × 5 km for all inhabited areas 30. The model was built using PM2.5 values from over 300 air quality monitoring stations in the National Air Pollution Surveillance (NAPS) program, which encompasses both rural and urban areas. The predictive variables include daily smoke plumes from remote sensing images, measurements of Aerosol Optical Depth, measurements of fire radiative power, and meteorological parameters such as temperature, humidity, and wind speed from satellite products 30. CanOSSEM demonstrates excellent performance, with approximately 98% of estimated concentrations within 5 µg/m3 of the observed PM2.5 value for the validation dataset30.
Daily estimated PM2.5 concentrations for each 6-digit postal code were derived from CanOSSEM. The estimates were then matched with the residential histories created from the MSP client roster and BCPDR records, resulting in daily maternal (from conception to the birthdate) and infant (from the birthdate to age 1 year) PM2.5 exposure history. The maternal and infant residential histories were created by identifying a continuous set of 6-digit postal codes covering every day of pregnancy or the first year of life. If there were duplicated records for any date, we always chose the 6-digit postal code that matched the 3-digit FSA in the BCPDR for the main analysis. As a sensitivity analysis, we always chose the 6-digit postal code that did not match the 3-digit postal code in the FSA, which resulted in different exposure estimates.
Outcomes
Diagnoses of otitis media and lower respiratory infections were identified from the MSP outpatient physician billings filed using the International Classification of Diseases, 9th Revision (ICD-9) codes. We used MSP billings for specific codes to identify otitis media (ICD-9 codes 381-382) and lower respiratory infections (ICD-9 codes 466, 480-487) in the first year after birth. Subsequent dispensations of amoxicillin were identified from PharmaNet. We only included dispensations of amoxicillin, which was the most prescribed antibiotic in the first year of life and it is also the preferred medication for treating infant otitis media and lower respiratory infections in Canada based on clinical guidelines 31,32. We included any amoxicillin dispensation within seven days of a physician visit for a respiratory infection.
Covariates
Analyses were adjusted for several potentially confounding variables from the BCPDR and MSP, reflecting maternal and infant health. Maternal covariates included: age at delivery; BMI category; maternal smoking status; and neighbourhood socioeconomic status (SES) quintile, all of which have been identified as potential confounders in prior research 33–35. Infant variables included indicators for: biological sex; preterm birth; and low birthweight34,36,37. Preterm birth refers to a newborn with a gestational age of less than 37 weeks37. Low birthweight is defined as an absolute weight less than 2500 g, regardless of gestational age 36. Missing values of categorical variables were addressed by including an unknown category.
Analyses were also adjusted for temperature and postnatal exposure to wildfire-related PM2.5 during the first year of life, which are also associated with an increased risk of respiratory illness 38–40. Exposure over the whole year was used for all infants, regardless of outcome status. We used the same temperature data included in CanOSSEM, which were derived from the NASA Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2) 30.
Exclusion criteria
We excluded infants missing 25% or more days of prenatal wildfire-related PM2.5 or temperature exposure data due to incomplete residential histories or missing CanOSSEM estimates. The proportion of missing data was calculated for the whole pregnancy, and for each critical exposure window. Similarly, infants missing more than 25% of postnatal exposure in the first year of life were also excluded. We also excluded infants who could not be matched with the MSP and PharmaNet datasets, and those missing values for maternal age.
Critical developmental windows
For otitis media, the critical windows of Eustachian tube development were defined as: 0-9 weeks (stage 1); 10-12 weeks (stage 2); 13-18 weeks (stage 3); 19-28 weeks (stage 4); and >28 weeks (stage 5). For the lower respiratory tract the developmental windows were defined as: 0-7 weeks (embryonic stage); 8-17 weeks (pseudoglandular stage); 18-27 weeks (canalicular stage); 28-36 weeks (saccular stage); and >37 (alveolar stage)41,42,43. We also examined the effects of each gestational trimester 0-13 weeks (Trimester 1); 14-27 weeks (Trimester 2) ; >27 weeks (Trimester 3)44,45. These critical time windows were identified from the existing literature based on significant physiological changes in both pregnant individuals and their foetuses (eTable 1).
Statistical analysis
All data preparation and statistical analysis were done using R 4.0.5 46. We used conditional logistic regression to quantify the association between prenatal wildfire-related PM2.5 exposure and respiratory infections by age 1 year. Models were conditioned on the epidemiologic week of birth to address confounding by seasonality, which is associated with both wildfire smoke exposures (higher in summer months) and respiratory infections (higher in winter months). Epidemiologic weeks allow comparison of population health trends between years, with the first epidemiologic week of the year beginning on the first Sunday in January. Some studies suggest matching by conception date rather than birthdate to address seasonality 47, and we did this as a sensitivity analysis.
We used the daily average estimates of wildfire-related PM2.5 concentration during each prenatal exposure window as the main predictor variable. The dependent outcome variables were all otitis media, all lower respiratory tract infections, and the otitis media or lower respiratory tract infections associated with amoxicillin dispensations. All estimates were adjusted for exposure during all stages of development, and the covariates described in the previous section (Eq 1).
Where Infection is a binary indicator of infection status by age 1 year; PM2.5-first and PM2.5-last are the average exposure in the first and last developmental windows or trimesters, respectively, and … indicates all windows between; postnatal PM2.5 is the average exposure during the first year of life; Temp-first and Temp-last are the average temperature the first and last developmental windows or trimesters, respectively, and … indicates all windows between; sex categorizes the biological sex of the infant as male, female or unknown; preterm is a binary indicator of birth with gestational week less than 37; low birthweight is a binary indicator of birth weight less than 2500 g; age indicates maternal age at delivery; smoking categorizes maternal smoking status as never smoker, former smoker, or current smoker; BMI indicates whether the mother was underweight, normal, overweight or obese; and SES represents maternal neighbourhood income quintiles (1 = most deprived to 5 = least deprived) during pregnancy.
Ethical approval
This work received approval from the UBC Research Ethics Board (Certificate number: H20-01077) and Population Data BC (approval number: 20-197).
Results
Descriptive summary
There were 164,415 singleton births in BC that were prenatally exposed to the peak wildfire seasons (July – September) of 2016 to 2019. Of these, 134,475 were matched to 121,280 unique mothers and included in the study (Figure 2). Maternal characteristics (Table 1) were calculated based on each pregnancy. The following 29,940 infants were excluded from the study (Figure 2): stillbirths (n=698); infants with more than 25% missing PM2.5 exposure data or more than 25% missing temperature data, and infants with missing PM2.5 exposure and temperature data at any critical developmental windows (n=27,848); infants that could not be matched to the MSP outpatient physician billings or PharmaNet files (n=1,393); Infants with missing maternal age (n=1).
Among all infants in the study, 16.4% (n=22,093) were diagnosed with one or more of the outcomes of interest (i.e., cases) and 83.6% (n=112,382) did not have any outcomes of interest (i.e., controls). There were 9,508 cases of otitis media, of which 6,395 were associated with amoxicillin dispensations. There were 14,464 cases of lower respiratory infections, of which 3,170 were associated with amoxicillin dispensations.
Significant differences (p-value <0.05) were observed in the distributions of all covariates and between cases and controls (Table 1). The biggest differences were male sex (57% cases, 50% controls), preterm births (3.7% cases, 2.7% controls), an ever smokers (13.4% cases, 11.6% controls). Small but statistically significant differences between groups were also observed in the PM2.5 and temperature exposures. The median (IQR) daily wildfire-related PM2.5 exposure during pregnancy was 7.5 (2.0) μg/m3 for cases and 7.4 (2.0) μg/m3 for controls (Table 1).
Wildfire-related PM2.5 exposure, temperature, and respiratory infections all showed clear seasonal trends by epidemiologic week across the study period (Figure 3). The wildfire-related PM2.5 and temperature peak in mid-summer, while the number of infections peak in the winter months. Furthermore, there were clear differences in wildfire-related PM2.5 exposure and outcomes among the infants prenatally exposed during the severe 2017 and 2018 wildfire seasons compared with those exposed during the 2016 and 2019 seasons (Figure 3).
Prenatal exposure to wildfire-related PM2.5 and otitis media
Exposure to wildfire-related PM2.5 during the fourth stage (19-28 weeks) of Eustachian tube development was significantly associated with an increased odds of otitis media diagnosis and associated amoxicillin dispensations. The odds ratio (OR) [95% confidence interval] for a 10 μg/m3 increase in average PM2.5 exposure during the fourth window was 1.08 [1.02, 1.16] for all diagnoses and 1.11 [1.02,1.20] for those associated with amoxicillin dispensations. A protective effect of wildfire-related PM2.5 was observed in both the first (0-9 weeks) and second (10-12 weeks) stages of Eustachian tube development (Figure 4).
Similarly, the second trimester (14-27 weeks) was significantly associated with increased odds of otitis media diagnosis and associated amoxicillin dispensations by age 1 year. The OR for a 10 μg/m3 increase in wildfire-related PM2.5 was 1.10 [1.03, 1.20] for otitis media diagnosis and 1.12 [1.03, 1.21] for diagnoses with associated amoxicillin dispensations. Once again, exposure during the first trimester (0 – 13 weeks) was associated with a protective effect (Figure 4). The sensitivity analyses using exposure estimates from alternate residential histories and conditioned on epidemiologic week of conception produced very similar results (eFigure 2).
Prenatal exposure to wildfire-related PM2.5 and lower respiratory infections
Exposure to wildfire-related PM2.5 during the canalicular stage (18-27 weeks) was significantly associated with an increased odds of lower respiratory infection diagnosis and associated amoxicillin dispensations. The odds ratio (OR) [95% confidence interval] for a 10 μg/m3 increase in average PM2.5 exposure during the canalicular stage was 1.08 [1.03, 1.14] for all diagnoses and 1.12 [1.01,1.25] for those associated with amoxicillin dispensations. Exposure during the pseudoglandular stage (8-17 weeks) was associated with a protective effect for lower respiratory infections associated with an amoxicillin dispensation. The second trimester (14-27 weeks) was significantly associated with increased odds of lower respiratory tract infection diagnosis by age 1 year with an OR of 1.06 [1.00, 1.12]. There were no significant associations observed for amoxicillin dispensations related to lower respiratory tract infection diagnosis for the analyses by trimester (Figure 5). Once again, the sensitivity analyses showed very similar results (eFigure 3).
Comments
Principle finding
This population-based case-control study examined wildfire-related PM2.5 exposure during critical windows of respiratory prenatal development and associations with respiratory infections by age 1 year. We found that wildfire-related PM2.5 exposure during the second trimester was associated with increased odds of otitis media and lower respiratory infections, including those treated with amoxicillin. We also consistently observed that wildfire-related PM2.5 exposure during the first trimester was associated with protective effects. Analyses based on critical windows of organ system development may provide more information than those based on trimester when examining the effects of environmental exposures during pregnancy.
Strengths
Our study has several strengths. First, we used the rich, population-based BCPDR to include most infants in BC who were prenatally exposed to wildfire smoke from 2016-2019. Because wildfire smoke was highly variable across these four years, we have excellent exposure contrasts within the cohort. Second, we were able to link BCPDR data with population-based MSP and PharmaNet data to identify all infections and those with associated amoxicillin dispensations. Third, we matched the analyses on epidemiologic week of birth to address confounding by seasonality, which is associated with both the exposure and the outcome. Fourth, we used a high-quality model of wildfire-related PM2.5 at high spatial resolution to assign daily residential exposures during pregnancy and the first year of life. Finally, we adjusted analyses for postnatal exposure to wildfire-related PM2.5, which can also be associated with the outcomes of interest.
Limitations
Our study used residential history as an indicator for exposure measurements, but we cannot account for mobility or errors in the residential history data, which leads to exposure misclassification and potential attenuation of the true effects. Second, we did not have individual-level SES information, and we used neighbourhood-level SES as approximation. Neighbourhood-level SES can even out the distribution of outcomes across income decile, and cause the estimates to be biased towards the null48. Third, we were unable to incorporate critical covariates identified in literature including maternal education and exposure to second-hand smoke because this information is not reliably captured in the BCPDR. Finally, all analyses were based on administrative data that are not collected for health research purposes. We lacked clinical details of each case, and we assume that physician billings for respiratory infections are a robust indicator of a true infection, especially when followed by an antibiotic dispensation.
Interpretation
Our findings are consistent with Hsu. et al., who reported that the mid-pregnancy period (16–25 weeks) is a critical window for the relationship between PM2.5 exposure and diagnosis with asthma by the age of six 49. Hazlehurst et al. also reported that higher PM2.5 exposure during the canalicular stage of lower respiratory tract development was linked to an increased risk of asthma 50. Another study identified the pseudoglandular and canalicular stages as the most critical periods for prenatal PM2.5 exposure and the subsequent development of asthma and wheezing after birth51. The critical developmental windows identified for respiratory outcomes are consistent with the sensitive exposure windows we report.
One other observational study specific to wildfire smoke found that the second and third trimesters were the most sensitive windows for use of upper respiratory medications before age 8 years52. We also observed the increased risk of otitis media associated with exposure during the second trimester. Another study using a mouse model found that maternal exposure to PM2.5 had a greater impact in the later stages of pregnancy, which could be related to the duration of exposure and accumulation of PM2.5.
Different from earlier studies, we consistently observed that wildfire-related PM2.5 exposure during the first stages of pregnancy was associated with protective effects for respiratory infections. We do not have an explanation for this observation, but we hypothesize that it could be related to a healthy survivor effect. Infants do not enter the BCPDR until 15 weeks’ gestation, and thus we have no information regarding the loss of earlier pregnancies. However, increased prenatal exposure to PM2.5 has been associated with increased risk of pregnancy loss53, and population exposure to smoke has been associated with a significant reductions in the size of birth cohorts54. There are also several studies indicating that elevated exposure to air pollution during preconception and early pregnancy is associated with spontaneous abortion prior to 20 weeks 55–57. If wildfire smoke exposure in early pregnancy leads to pregnancy loss, the foetuses that survive may develop into healthier infants and children 58.
Our results suggest that prenatal exposure to wildfire-related PM2.5 may interfere with foetal development. During development of the Eustachian tubes, the fourth stage (19-28 weeks) emerges as a pivotal period characterized by differentiation of certain seromucous glandular extensions? within the lumens, which is associated with maturation of submucosal gland 41. The glandular mucus secreted by the submucosal gland aids in the removal of inhaled particles and provides protection to the respiratory airway by combating bacterial threats59.
We hypothesize that prenatal exposure to wildfire-related PM2.5 during the fourth stage of Eustachian tube development may be associated with inadequate submucosal gland development. The compromise in the protective functioning of the Eustachian tube potentially increases susceptibility to upper respiratory infections, such as otitis media, during early life. Similarly, the canalicular stage of foetal lower respiratory tract development is associated with the differentiation of type I and type II cells, as well as formation of the alveolar-capillary barrier. Insufficient maturation during the canalicular stage may result in underdeveloped peripheral airways, hindering efficient gas exchange and respiratory balance. This, in turn, can contribute to the proliferation of alveoli populated by inflammatory cells 19,60, which might make infants more prone to lower respiratory infections.
We identified a significant association between prenatal exposure to wildfire-related PM2.5 and respiratory infections in infants, as evidenced by increased odds of infection diagnoses and related antibiotic dispensations. To mitigate the detrimental impact of wildfire smoke on infant respiratory health, clear communication regarding the risks and protective measures is needed for those who are pregnant. For example, use of indoor air cleaning technologies should be strongly recommended for this population 61. Future studies should focus on the biological mechanisms underlying the effects of prenatal exposure to wildfire smoke, including epigenetics and markers of inflammation. Additionally, extending the follow-up period of the cohort will enable us to examine the longer-term effects associated with prenatal exposure to wildfire smoke.
Data Availability
Access to data provided by the Data Stewards is subject to approval but can be requested for research projects through the Data Stewards or their designated service providers. All inferences, opinions, and conclusions drawn in this publication are those of the author(s), and do not reflect the opinions or policies of the Data Stewards (the British Columbia Ministry of Health, the Perinatal Services BC).
Author Contributions
Qingyi Lan: Conceptualization; formal analysis; writing-original draft preparation; writing-review and editing. Kate Weinberger: Conceptualization; writing-review and editing. Sabrina Luke: Conceptualization; writing-review and editing. Eric Lavigne: Conceptualization; writing-review and editing. Scott Weichenthal: Conceptualization; writing-review and editing. Sarah B. Henderson: Supervision; conceptualization; writing-reviewing and editing.
Conflict of Interest disclosure
None.
Disclaimer
Access to data provided by the Data Stewards is subject to approval but can be requested for research projects through the Data Stewards or their designated service providers. All inferences, opinions, and conclusions drawn in this publication are those of the author(s), and do not reflect the opinions or policies of the Data Stewards (the British Columbia Ministry of Health, the Perinatal Services BC).
Acknowledgements
This study was funded by Canadian Institutes of Health Research Grant (CIHR). Qingyi Lan received Graduate Research Stipend from the CIHR, International Tuition Award, President’s Academic Excellence Initiative PhD Award, and Faculty of Medicine Graduate Award.