- Open Access
Increased richness and diversity of the vaginal microbiota and spontaneous preterm birth
Microbiomevolume 6, Article number: 117 (2018)
The bacterial community present in the female lower genital tract plays an important role in maternal and neonatal health. Imbalances in this microbiota have been associated with negative reproductive outcomes, such as spontaneous preterm birth (sPTB), but the mechanisms underlying the association between a disturbed microbiota and sPTB remain poorly understood. An intrauterine infection ascending from the vagina is thought to be an important contributor to the onset of preterm labour. Our objective was to characterize the vaginal microbiota of pregnant women who had sPTB (n = 46) and compare to those of pregnant women who delivered at term (n = 170). Vaginal swabs were collected from women at 11–16 weeks of gestational age. Microbiota profiles were created by PCR amplification and pyrosequencing of the cpn60 universal target region.
Profiles clustered into seven community state types: I (Lactobacillus crispatus dominated), II (Lactobacillus gasseri dominated), III (Lactobacillus iners dominated), IVA (Gardnerella vaginalis subgroup B or mix of species), IVC (G. vaginalis subgroup A dominated), IVD (G. vaginalis subgroup C dominated) and V (Lactobacillus jensenii dominated). The microbiota of women who experienced preterm birth (< 37 weeks gestation) had higher richness and diversity and higher Mollicutes prevalence when compared to those of women who delivered at term. The two groups did not cluster according to CST, likely because CST assignment is driven in most cases by the dominance of one particular species, overwhelming the contributions of more rare taxa. In conclusion, we did not identify a specific microbial community structure that predicts sPTB, but differences in microbiota richness, diversity and Mollicutes prevalence were observed between groups.
Although a causal relationship remains to be determined, our results confirm previous reports of an association between Mollicutes and sPTB and further suggest that a more diverse microbiome may be important in the pathogenesis of some cases.
Preterm birth is defined as delivery before 37 completed weeks of gestational age  and can be further sub-categorized in extremely preterm (≤ 27+6 weeks+days), very preterm (28 to 31+6) and late preterm (32 to 36+6) . Preterm birth comprises 11% of all livebirths worldwide, and its complications are estimated to cause 35% of world’s neonatal deaths, which represents 3.1 million deaths annually . Children who are born prematurely also have higher rates of cardiovascular disorders, respiratory distress syndrome, neurodevelopmental disabilities and learning difficulties compared with children born at term .
Preterm birth is a complex multi-factorial condition with several known risk factors, such as low and high maternal ages [5,6,7], low BMI , black ethnicity , tobacco use, heavy alcohol intake, illicit drug use , close temporal proximity to a previous delivery , and multiple gestation . Although studied extensively, some preterm cases remain unexplained for women with no known risk factors. Intrauterine infection with organisms ascending from the vagina has been hypothesized as an important contributor to preterm birth since many organisms isolated from the amniotic fluid/membranes of women who experienced preterm birth are also found in the lower genital tract of pregnant women [12,13,14,15]. A large number of studies support this hypothesis based on the strong association between intra-amniotic bacterial infection and preterm birth [12, 14,15,16,17,18,19,20].
The microbiological diagnosis of a ‘normal’ or disturbed vaginal microbiota has historically been based on the Nugent score, the current gold standard diagnostic method that relies on Gram stain of vaginal smears . The ‘normal’ vaginal microbiota in non-pregnant reproductive aged women is understood to be dominated by Lactobacillus species, while an abnormal microbiota (defined as bacterial vaginosis) is characterized by low abundance of lactobacilli and an overgrowth of anaerobic bacteria, such as Gardnerella vaginalis, Prevotella spp., Bacteroides spp., Mobiluncus spp. and Mycoplasma hominis . In low-risk pregnant women, it has been shown that the vaginal microbiota has reduced richness and diversity and increased abundance of lactobacilli compared to those of non-pregnant women [23,24,25,26,27]. An abnormal microbiota has been previously associated with preterm birth , but only a few in depth culture-independent studies of the vaginal microbiota of women who had preterm birth have been published, with inconsistent conclusions [29,30,31,32].
The objective of this study was to assess whether there are differences in the vaginal microbiota composition, early in gestation, of women who had spontaneous preterm birth (sPTB) and term delivery that could be further investigated as diagnostic indicators of preterm birth risk. Microbiome profiling was based on sequencing of the cpn60 universal target, which provides higher resolution than 16S rRNA variable regions  and allows the resolution of Gardnerella vaginalis subgroups, a hallmark bacteria in the disturbed microbiota .
Study population and sampling
This retrospective cohort study analysed the vaginal microbiota of women who experienced spontaneous preterm birth (sPTB) and compared the resulting microbial profiles to those of pregnant women who delivered at term. The bacterial profiles of pregnant Canadian women at low risk of sPTB who had term deliveries (n = 170) were previously generated by our research group . The vaginal microbial profiles of Canadian women who had preterm birth originated from samples of this previous study (n = 7)  and from the Ontario Birth Study (n = 39), resulting in 46 samples. The Ontario Birth Study (ontariobirthstudy.com) is an open longitudinal pregnancy cohort at Mount Sinai Hospital, Toronto, Canada. It is a platform for studies of both pregnancy complications as well as Developmental Origins of Health and Disease related research. The PTB rates for the low-risk cohort and OBS cohorts were 4 and 6.2%, respectively. All biospecimens, including maternal vaginal swabs and maternal and infant blood, are collected concurrently with routine clinical specimens to reduce the burden on study participants. Detailed demographic and lifestyle characteristics are obtained from women during pregnancy and postpartum, and clinical information is extracted from the health records. For the purposes of this report, self-administered vaginal swabs were taken at 16 weeks gestation and placed in dry tubes prior to being placed in − 80 °C for storage in the Lunenfeld Tanenbaum Research Institute Biospecimen Storage and Processing Laboratory. Specimens from all cohorts were processed similarly in terms of sample collection, storage, DNA extraction, library preparation and sequencing.
Clinical and behavioural questionnaire data (pregnancy history, family and personal medical history, psychosocial health, demographic factors and other lifestyle and environmental exposures) were transferred to the Research Electronic Data Capture (REDCap) database protected by a secure server . For the PTB group, eligible participants for this study were women who had undergone preterm delivery at greater than 20 weeks but less than 37 weeks gestational age, where onset of labour occurred spontaneously or in association with cervical incompetence or preterm premature rupture of membranes (PPROM). Vaginal swabs collected from pregnant women (both PTB and term groups) at 11–16 weeks of gestational age were used for bacterial genomic analysis.
Total nucleic acid was extracted from swabs using the MagMAX™ Total Nucleic Acid Isolation Kit (Life Technologies, Burlington, ON, Canada) as per manufacturer’s instructions. Kit reagents are aliquoted to eliminate repeated accessing of open reagents, and samples are processed in small batches using filter tips to prevent cross-contamination. Pipettes and other lab surfaces are regularly treated with DNA surface decontaminant (DNA Away, Thermo Fisher Scientific, Waltham, MA). Samples from both cohorts were processed in exactly the same way in terms of swab type, storage temperature (no stabilizer was used), DNA extraction, library preparation and sequencing.
Total bacterial DNA (qPCR) and detection of Mollicutes (PCR)
Quantitative PCR (qPCR)
Total bacterial DNA quantity in each sample was estimated using a SYBR Green assay based on amplification of the V3 region of the 16S rRNA gene. Primer sequences were as follows: SRV3-1 (5′-CGGYCCAGACTCCTAC-3′), SRV3-2 (5′-TTACCGCGGCTGCTGGCAC-3′) . Reactions were run on a MyiQ thermocycler using the following cycling parameters: 95 °C for 3 min, followed by 30 cycles of 95 °C for 15 s, 62 °C for 15 s and 72 °C for 15 s, with a final extension at 72 °C for 5 min .
Some Mollicutes (Mycoplasma and Ureaplasma) species lack a cpn60 gene . Thus, we performed a family-specific semi-nested PCR targeting the 16S rRNA gene to detect Mollicutes , and a PCR targeting the multiple-banded antigen gene to detect Ureaplasma spp. PCR products from U. parvum and U. urealyticum can be differentiated by size .
cpn60 universal target (UT) PCR and pyrosequencing
Universal primer PCR targeting the 549–567 bp cpn60 UT region was performed using a mixture of cpn60 primers consisting of a 1:3 M ratio of primers H279/H280:H1612/H1613, as described previously [41,42,43]. To allow multiplexing of samples in a single sequencing run, primers were modified at the 5′ end with one of 24 unique decamer multiplexing identification (MID) sequences, as per the manufacturer’s recommendations (Roche, Brandford, CT, USA). Amplicons were pooled in equimolar amounts for sequencing on the Roche GS Junior sequencing platform. The sequencing libraries were prepared using the GS DNA library preparation kit, and emulsion PCR (emPCR) was performed with a GS emPCR kit (Roche Diagnostics, Laval, Canada).
Samples were handled in small batches to avoid cross-contamination, and experimental controls were included at several steps in the study. Regular monitoring of DNA extraction controls in our lab by universal PCR confirms that these procedures are sufficient to eliminate detectable template contamination of study samples. A no template control was also included in each set of PCR reaction as negative controls. Experimental controls were not sequenced as they did not yield any amplification.
Analysis of operational taxonomic units (OTU)
Raw sequence data was processed by using the default on-rig procedures from 454/Roche. Filter-passing reads were used in the subsequent analyses for each of the pyrosequencing libraries. MID-partitioned sequences were mapped using Bowtie 2 (http://bowtie-bio.sourceforge.net/bowtie2/) on to a manually curated reference set of 1561 OTU sequences representing the human vaginal microbiota. Bowtie 2 was run using the default end-to-end alignment mode.
The OTU reference set was generated originally by de novo assembly of cpn60 sequence reads from 546 vaginal microbiomes using the microbial Profiling Using Metagenomic Assembly pipeline (mPUMA, http://mpuma.sourceforge.net)  with Trinity as the assembly tool  (Additional file 1). OTU were labeled according to their nearest reference sequence determined by watered-Blast comparison  to the cpn60 reference database, cpnDB_nr (downloaded from http://www.cpndb.ca ). This reference assembly approach allows us to compare the microbial profiles from various cohorts under investigation, including the 46 pregnant women who had sPTB described in this study.
The result of mapping is an OTU frequency table (Additional file 2) that was used for microbiome data analysis. Some analyses were also performed at species level, i.e. combined OTU that have the same nearest neighbour.
Comparisons of socio-demographic characteristics of cohorts and participants were based on analysis of variance (ANOVA), t test and chi-square, performed in IBM SPSS (Statistical Package for the Social Sciences, version 21) at 5% level of significance. For analysis of associations between socio-demographic characteristics and microbiota profiles (CST), a false discovery rate (FDR) correction for multiple comparisons was applied .
Alpha (Shannon diversity and Chao1 estimated species richness) and beta diversity (jackknifed Bray–Curtis dissimilarity matrices) were calculated as the mean of 100 subsamplings of 1000 reads (or all reads available when less than 1000) in QIIME (Quantitative Insights Into Microbial Ecology) . Plots of alpha diversity measures against bootstrap sample number were generated in R and visually inspected to ensure that an adequate sampling depth for each sample was achieved.
For community state type (CST) analysis, a Jensen–Shannon distance matrix was calculated using the ‘vegdist’ function in the vegan package  with a custom distance function that calculates the square root of the Jensen–Shannon divergence . This distance matrix was used for hierarchical clustering using the ‘hclust’ function in R, with Ward linkage.
The function aldex.clr from the ALDEx2 package in R was used to compare the differential relative abundance of individual taxa in term and preterm groups . Significant differences were determined based on the false discovery rate (FDR), which is the result of a Benjamini–Hochberg corrected p value from a Welch’s t test calculated within ALDEx2.
Description of the study population and pregnancy outcomes
Socio-demographic characteristics of women who had spontaneous preterm birth (n = 46) and women who had term deliveries (n = 170) are summarized in Table 1. There were no significant differences in maternal age, BMI, ethnicity, smoking status, consumption of alcohol or use of probiotics between term and preterm groups (all p > 0.05). Average maternal age was 33 for participants in both cohorts. Average body mass index (BMI) was 22.9 and 24.2 for women in the term and preterm groups, respectively. Most women in both cohorts identified themselves as white ethnicity, followed by East Asian and South/Southeast Asian (Table 1). Consumption of tobacco (term 2.3%; preterm 0%), alcohol (term 5.9%; preterm 4.3%) or probiotic supplements (term 4.1%; preterm 6.5%) was low among women in both groups (chi-square, all p > 0.05).
Most women in the preterm group had a Bachelor/graduate degree (29/46) and an average house income higher than CAD 100,000 per year (25/46). A minority of women who had preterm birth (5/46) reported consumption of substances without prescription prior pregnancy, of which 3/46 women consumed marijuana/hashish, 1/46 woman consumed tranquilizers/nerve pills and 1/46 woman consumed cocaine/crack. Approximately 74% of the participants in the preterm group reported a pre-existing condition. A total of 12/46 women had some condition related to mental health, such depression or anxiety. Seventeen percent (8/46) had a neurological condition, including migraine headaches, and 24% (11/46) had a genitourinary condition (corpus luteal cyst, bicornuate uterus, cervical polyp, cervical dysplasia (2), uterine polyp, ovarian cyst, polycystic ovarian syndrome, urinary tract infections with and without kidney stones (3).
Characteristics regarding pregnancy and neonatal outcomes are described in Table 2. Pregnancy outcome information was not available for one woman in the preterm group as she was lost to follow-up. There were no significant differences in gestational age at enrolment, mode of conception or fetal sex between groups (all p > 0.05). Average gestational age at delivery was 39+3 weeks for the women who delivered at term and 34+2 weeks for women who had preterm birth, most of which were considered late preterm, i.e. delivery between 32 and 36+6 weeks of gestational age. Women in the preterm group were more likely to have experienced preterm birth or miscarriage in their previous pregnancy (chi-square, p < 0.001). They also had higher percentage of caesarean sections than women who delivered at term. Number of previous gestations also differed between groups; women who had preterm birth were more likely to be primigravida (22/46) in comparison with women who had term deliveries (45/170). There was a significant difference between term and preterm groups regarding birth weight and number of infants admitted to level 3 neonatal intensive care unit (NICU) (Table 2). Apgar score at 1 (term 8.75 ± 0.6; preterm 8.38 ± 1.1) and 5 min (term 8.97 ± 0.17; preterm 8.76 ± 0.7) between groups also differed (t test, all p < 0.001). One preterm infant (1/46) died shortly after birth (20 weeks of gestational age).
Among women who delivered preterm, 63% (29/46) had premature rupture of membranes (PPROM), 10.8% (5/46) had gestational diabetes and 4.3% (2/46) had anemia unresponsive to therapy. Twenty-four percent of women (11/46) presented one of the following conditions: maternal elevated liver enzymes, short cervix and incompetent cervix; fetal ascites, fetal distress and large foetus for gestational age; and placental findings of marginal cord insertion, two-vessel umbilical cord, placenta previa and low-lying placenta.
Sequencing results and OTU analysis
Raw sequence data files for the samples described in this study were deposited to the NCBI Sequence Read Archive (Accession SRP073152, BioProject PRJNA317763; BioProject PRJNA403856). Total dataset contained 1,635,072 cpn60 reads; median and average read count per sample were 4936 and 7569 (range 402–37,378), respectively. Average read length was 424 bp. Results of Bowtie2 mapping showed that these reads corresponded to 728 OTUs from the reference assembly (Additional file 1).
Microbiota profiles were created by PCR amplification and pyrosequencing of the cpn60 universal target region. Hierarchical clustering of vaginal microbiota profiles resulted in seven community state types (CST): I (Lactobacillus crispatus dominated), II (Lactobacillus gasseri dominated), III (Lactobacillus iners dominated), IVA (Gardnerella vaginalis subgroup B or mix of different species), IVC (G. vaginalis subgroup A dominated), IVD (G. vaginalis subgroup C dominated) and V (Lactobacillus jensenii dominated) (Fig. 1). Each CST is defined by the dominance of one species of Lactobacillus (I, II, III, V), Gardnerella vaginalis (IVC, IVD) or a mixture of bacteria species (IVA), as previously described [52, 53].
Overall microbiota profiles did not cluster together based on gestational age at delivery (Figs. 1 and 2). Most microbial profiles from the preterm group (80.5%) were assigned to Latobacillus-dominated CST: CST I (37% of profiles), CST III (17.4%), CST V (15.2%) and CST II (10.9%). The remaining profiles (19.5%) were assigned to CST IVA, IVC or IVD (Table 1). The CST IVA was the most heterogeneous group, represented by the dominance of Lactobacillus delbrueckii, Bifidobacterium dentium, Bifidobacterium infantis, Atopobium vaginae, Bifidobacterium breve or a mixture of different bacteria species. The CST IVC was dominated by G. vaginalis subgroup A and Megasphaera spp., and CST IVD was dominated by G. vaginalis subgroup C (Fig. 1).
Ecological analysis and total bacterial load
Assessment of alpha diversity revealed that microbiomes of women who delivered preterm were richer (Chao1 richness 46.3 ± 24.1) and more diverse (Shannon diversity index 1.8 ± 1.1) when compared to those of women in the term group (36.2 ± 14.8; 1.2 ± 0.8) (t test, p < 0.01) (Table 1). Total bacterial load was estimated based on qPCR targeting the 16S rRNA gene, and it was expressed as log 16S rRNA gene copy number per swab. Higher bacterial loads were detected in samples from the preterm group (7.7 ± 0.9) compared to term group (8.0 ± 0.7) (t test, p = 0.049) (Table 1).
Bacteria species relative abundance and prevalence
To investigate whether there was an association between individual taxa and sPTB, the abundance and prevalence of each species was evaluated. The ALDEx2 analysis assessed the relative abundance of each taxa (at the OTU and species level) in term and preterm groups. Eight OTU/species were more abundant in the term group in comparison with preterm, all of which were considered rare members of the bacterial community (Fig. 3). L. acidophilus represented 1% of the total reads in the dataset and had a low relative abundance average of 1.98% (range 0–69%) and 0.18% (range 0–0.87%) in samples from term and preterm groups respectively. All the other seven bacteria together represented only 0.4% of the total reads in the dataset.
Bacteria prevalence (presence/absence) was also assessed (only species with at least 10 total reads were included). A total of 60 taxa had significant differences in prevalence between term and preterm groups; 11 species had greater prevalence in the term cohort and 49 species were more prevalent in the preterm cohort (Table 3). Bifidobacterium infantis, for example, was two times more prevalent in the term group in comparison with preterm, and Prevotella timonensis was 1.58 times more prevalent in the preterm group (Table 3). Several Prevotella spp. were associated with both term and preterm. Prevotella amnii and P. tannerae had greater prevalence in the term cohort, whereas P. timonensis, P. bivia, P. corporis and P. bucalis were more prevalent in the preterm group (Table 3). It is important to note that read depth distribution did not differ between term and preterm cohorts (t test, p > 0.05); therefore, the differences observed here in bacteria prevalence were unlikely to be driven by cohort sequencing bias.
Mollicutes (Mycoplasma and/or Ureaplasma) were detected by family-specific conventional PCR in 28/46 (60%) of pregnant women who delivered preterm (Table 1). Ureaplasma species were detected by genus-specific PCR in samples of 14/46 (30%) women who had PTB, with all women testing positive for U. parvum and none for U. urealyticum. Women who delivered at term were less likely to be PCR positive for Mollicutes compared to women who had PTB (Table 1). No significant differences were observed in Ureaplasma prevalence between the two groups (Table 1). An association between Mollicutes/Ureaplasma detection and the composition of the vaginal microbiota, represented as CST, was also investigated. Detection of Mollicutes and Ureaplasma was not associated with any CST in particular when investigated in the term cohort, preterm cohort or both groups together (chi-square, p > 0.05).
Relationships between microbiological and socio-demographic characteristics within the preterm group
The association between CST (I, II, III, IVC, IVD, V) from profiles of women who delivered preterm and several microbiologic-socio-demographic characteristics was investigated. Only two associations were significant: CST and microbiota richness, and CST and microbiota diversity (ANOVA, p < 0.001). There was no significant association between CST and the remaining following metadata: microbiota richness and diversity (continuous variable), presence of Mollicutes and Ureaplasma (yes/no), log 16S rRNA gene copies (continuous variable), maternal age (continuous variable; 18–25, 26–35, 36–45), BMI category (underweight, normal, overweight, obese; < 25, ≥ 25), ethnicity (White, East Asian, South Asian, Black, Hispanic, Other), natural conception (yes/no), parity (0, > 1), gestational age (continuous variable), mode of delivery (vaginal, C-section), pre-existing condition (yes/no), folic acid intake before or during pregnancy (yes/no), drinking alcohol (yes/no), neonate in high level care (yes/no), birth weight (continuous variable) and Apgar score at 1 and 5 min (1–9).
In this study, we determined the composition of the vaginal microbiota of women who had spontaneous preterm birth and compared these profiles to those of women who delivered at term, previously reported by our research group . The availability of foundational data on women who delivered at term and the infeasibility of collecting large numbers of samples at 11–16 weeks gestation from women who would go on to deliver pre-term, our study design included comparison of samples collected in a previously published study  and from the OBS. To minimize any batch effects, we were rigorous in implementation of consistent sample processing and did extensive analysis of the clinical and demographic characteristics to ensure they were well matched (Table 1). The cohorts were comparable in terms of maternal age, BMI, ethnicity, consumption of tobacco, alcohol and probiotics, which is of interest given that several of these characteristics have been previously associated with preterm delivery. In particular, previous described factors included low and high maternal ages [5,6,7], low BMI , black ethnicity , high levels of tobacco, alcohol and illicit drugs consumption , close temporal proximity to a previous delivery  and multiple gestation . This cohort is unique in that it did offer the opportunity to have gestational age at delivery as the main characteristic distinguishing these two groups recognizing that the majority of preterm births occurred beyond 32 weeks gestation.
A difference in number of previous gestations was observed between groups, with women who experienced preterm birth more likely to be primigravida in comparison with women who had term deliveries. It has been recently demonstrated that women with a prior conception, regardless of whether or not this proceeded to a birth, have a decrease in the relative abundance of L. crispatus and a concomitant increase in the abundance of other Lactobacillus species as well as Gardnerella .
Other known risk factors for sPTB include maternal medical disorders like hypertension, asthma, diabetes and thyroid disease . Although some women in both cohorts reported these conditions, there were not enough participants to stratify the data based on the individual disorder and therefore was not possible to investigate the interaction between those medical conditions and gestation outcome. We were, however, able to confirm previous reports of history of prematurity as a risk factor for preterm birth .
Since many organisms isolated from the amniotic cavity of women who experienced preterm birth are also found in the genital tract [12,13,14,15], an intrauterine infection ascending from the vagina is one of the currently hypothesized triggers of PTB . In this study, however, we did not identify a signature microbiota composition (CST) associated with preterm birth. This observation is consistent with the results presented by others [29, 30]. CST assignments are largely driven by the dominance of a single species, which may mask differences in rare taxa that would differentiate term and preterm groups, and indeed, further analysis revealed that the vaginal microbiota of women who experienced preterm birth was richer and more diverse than those of women who delivered at term. Also, most women (84.8%) in our study were considered late preterm and although we cannot address this question, it is possible that sPTB driven by an ascending infection would be more evident in a high-risk cohort or extreme preterm cases. A recent study of a high-risk pregnant cohort has reported that L. iners was strongly associated with short cervix and preterm birth, as L. crispatus was associated with term deliveries . Those differences in study outcomes indicate that the pathogenesis of sPTB in low- and high-risk groups might be different. Identifying differences in the causes of early and late sPTB and the role of the vaginal microbiota in those processes will require further study.
One controversy that challenges the current hypothesis of preterm caused by an ascending infection is that antibiotic administration to pregnant women with a disturbed vaginal microbiota does not improve outcome in most cases, as demonstrated by study trials [58, 59] and systematic reviews [60,61,62]. One explanation for the inefficacy of antibiotic treatment in the prevention of preterm birth relies is the high rates of antibiotic resistance among bacterial-vaginosis-associated bacteria [63, 64]. In this case, antibiotics not only do not kill the targeted bacteria, but might also reduce the vaginal Lactobacillus population leading to an even more disturbed microbiota, as recently demonstrated .
In addition to differences in richness and diversity, differences in the microbiota between the two cohorts regarding bacterial abundance and prevalence were also identified. The ALDEx2 analysis indicated that eight rare taxa were more abundant in the term group, which does not necessarily mean they are associated with a ‘healthier’ state or implicated in preventing sPTB. Since these bacteria are detected at very low abundance within the microbiota profiles, their biological significance in the vaginal microbiome is questionable. Differences in the prevalence of several other taxa between groups were also observed. For example, more women in the term group had Prevotella amnii and P. tannerae detected in their vaginal samples, whereas P. timonensis, P. bivia, P. corporis and P. bucalis were more frequently detected in samples from women in the preterm group (Table 3). Prevotella spp. have been previously associated with bacterial vaginosis and preterm labour [22, 66, 67], and our results indicate that different Prevotella species might have different roles in sPTB. Several of the taxa that were significantly different in their prevalence among women in the two groups also had low sequence read counts (Table 3). Further investigation would be required to determine if these rare members of the microbial community play a yet unknown role in sPTB.
It is also important to note that the number of bacterial species with greater prevalence in the preterm (49/60) was higher than in the term (11/60) cohort (Table 3), which is consistent with our results of increased microbial richness and diversity in the samples from women who experienced preterm birth. This might indicate that increased richness, rather than the presence of specific taxa, might be associated with sPTB. Those differences might also be an indicator of physiological/biochemical dissimilarities in the vaginal microbiomes of women who deliver at term or preterm. In other words, the physiological state that leads to sPTB might also create an environment that supports a richer/more diverse microbiota.
Our results also confirmed previous reports of an association between Mycoplasma and preterm birth . Mollicutes were detected significantly more often in women in the preterm group compared to women in the term group, but no differences were observed in Ureaplasma prevalence between groups indicating that the difference in Mollicutes prevalence is primarily driven by the presence of Mycoplasma spp. Although individual Mycoplasma species could not be discerned based on assay used in our study, both Mycoplasma genitalium [69,70,71] and Mycoplasma hominis [72,73,74,75] have been previously associated with negative reproductive outcomes including PTB.
Collectively, our overall findings were similar to other two studies, which provided us the opportunity to compare different study designs (based on different cohorts and barcode gene) that addressed the same research question. Hyman and colleagues  described the vaginal microbiota of 83 pregnant women (term n = 66, preterm n = 17) based on Sanger sequencing of cloned 16S rRNA genes. Samples were collected at each trimester and preterm was defined as delivery before 37 weeks of gestation. There was no correlation between preterm and absence/low abundance of Lactobacillus in the microbiota; in other words, preterm outcome could not be predicted based on CST. Similar to our results, they found an association between increased microbiota diversity and preterm delivery among women of white ethnicity (n = 40) (data from women of others ethnicities was not included in the analysis because of small sample sizes). Although there was no association between CST and ethnicity, it is important to note that most women enrolled in this study described themselves as being white, and it is possible that an increased sample size of participants of other ethnicities could result in a different conclusion.
Romero and colleagues  also investigated the vaginal microbiota of pregnant women who experienced preterm, defined as delivery before 34 weeks of gestation (term n = 72, preterm n = 18). The profiles were created by 16S rRNA amplicon sequencing, and samples were collected every 4 weeks until 24 weeks of gestation and then every 2 weeks. They found no differences in the frequency of different CST between women who had term and preterm deliveries. Likewise, no differences in bacteria relative abundance were observed between the two cohorts, although only bacteria that were present in at least 25% of samples were included in the analysis. These results are consistent with our findings of bacterial abundance based on the ALDEx analysis since we only found significant differences in relative abundance for eight rare bacteria. Unlike Hyman et al.  and our results, Romero et al.  did not find differences in microbiota diversity between women who delivered preterm and at term. One possible explanation for this contradictory result might be related to differences in participant ethnicity among these studies. While most women in our study and the Hyman et al. study described themselves as white, the majority of participants in the Romero et al. study described themselves as African American. It has been reported that the composition of the vaginal microbiota is strongly associated with a woman’s ethnicity [52, 76]. Other studies have also demonstrated that black ethnicity is associated with an increased microbiota diversity in comparison with white ethnicity , which could have masked differences in bacterial diversity between term and preterm cohorts in the Romero study.
Contrary to our overall findings, DiGiulio and colleagues  found a strong association between the non-Lactobacillus-dominated CST IV and preterm birth in a case–control study based on the 16S rRNA amplicon sequencing. Pregnant women (preterm n = 34, term n = 15), mostly of white ethnicity, were sampled weekly throughout gestation. Interestingly, the authors pointed out that if samples had been collected less frequently, short-term ‘excursions’ to CST IV would have been missed and probably the association between CST IV and preterm birth would have been less obvious. The detection of a temporary microbiota disturbance represented by a change from a Lactobacillus-dominated CST to CST IV may have been missed in our study since samples were not collected longitudinally. Moreover, a recent study has demonstrated that PTB–microbiota associations are population-dependent ; lower Lactobacillus and higher Gardnerella abundances were associated with PTB in a low-risk predominantly Caucasian cohort, but not in a high-risk predominantly African American cohort. These population-dependent associations might contribute to explain contradictory conclusions among different studies and emphasize the importance of investigating the vaginal microbiota of different populations with varying ethnic backgrounds and from different geographical locations.
Most samples in the preterm group were dominated by Lactobacillus, yet, they collectively had higher richness and diversity compared to samples from the term group. The increased microbiota richness/diversity might indicate a transient state between Lactobacillus-dominated CST and non-Lactobacillus-dominated, i.e., CST IV (A, C or D). In other words, the increased richness and diversity we observed might be a remnant characteristic of the previous disturbed microbiota. In summary, although we did not “detect” a specific microbial community structure that is associated with preterm birth, the increased microbiota richness/diversity was associated with preterm birth. In addition, the association with differences in Prevotella species and Mycoplasma presence may point to signature species associated with preterm birth.
Taken together, our results suggest that the differences in the microbiota of women who had preterm deliveries, such as increased microbiota richness and diversity and greater prevalence of Mollicutes and other bacteria, may have a role in sPTB. Other differences between cohorts might have been masked by the presence of highly dominant bacteria like Lactobacillus. At the overall level, we did not identify a specific vaginal microbial community structure at 11–16 weeks gestation age that predicts sPTB. Also, differences in relative abundance of bacterial species between term and preterm groups were only significant for a few low abundance species. Although a causal relationship remains to be determined, our results confirm previous reports of an association between Mollicutes and preterm birth, and further suggest that a diverse bacterial community may contribute to the microbiome’s role in sPTB. Alternatively, the more rich and diverse microbiotas of the preterm group may reflect physiological differences between the groups that affect selection of bacteria. This study provides valuable evidence of subtle alterations in the microbiome associated with preterm birth that requires further study utilizing sequencing methodology. In addition, future study should include evaluation of the microbial metabolite production and host response to further elucidate factors leading to sPTB and identify women at risk early in pregnancy.
WHO, March of Dimes, PMNCH, Save the Children. In: Howson CP, Kinney MV, Lawn JE, editors. Born too soon: the global action report on preterm birth. Geneva: World Health Organization; 2012.
Quinn J, Munoz FM, Gonik B, Frau L, Cutland C, Mallett-Moore T, et al. Preterm birth: case definition & guidelines for data collection, analysis, and presentation of immunisation safety data. Vaccine. 2016;34:6047–56.
Blencowe H, Cousens S, Oestergaard MZ, Chou D, Moller AB, Narwal R, et al. National, regional, and worldwide estimates of preterm birth rates in the year 2010 with time trends since 1990 for selected countries: a systematic analysis and implications. Lancet. Elsevier Ltd. 2012;379:2162–72.
Behrman RE, Butler AS. In: Behrman RE, Butler AS, editors. Preterm birth: causes, consequences, and prevention. Washington, D.C.: National Academies Press; 2007.
Fraser AM, Brockert JE, Ward RH. Association of young maternal age with adverse reproductive outcomes. N Engl J Med. 1995;332:1113–8.
Jacobsson B, Ladfors L, Milsom I. Advanced maternal age and adverse perinatal outcome. Obstet Gynecol. 2004;104:727–33.
Schempf AH, Branum AM, Lukacs SL, Schoendorf KC. Maternal age and parity-associated risks of preterm birth: differences by race/ethnicity. Paediatr Perinat Epidemiol. 2007;21:34–43.
Han Z, Mulla S, Beyene J, Liao G, McDonald SD. Maternal underweight and the risk of preterm birth and low birth weight: a systematic review and meta-analyses. Int J Epidemiol. 2011;40:65–101.
Kistka ZAF, Palomar L, Lee KA, Boslaugh SE, Wangler MF, Cole FS, et al. Racial disparity in the frequency of recurrence of preterm birth. Am J Obstet Gynecol. 2007;196:1–6.
Conde-Agudelo A, Rosas-Bermúdez A, Kafury-Goeta AC. Birth spacing and risk of adverse perinatal outcomes. JAMA. 2006;295:1809–23.
Stock S, Norman J. Preterm and term labour in multiple pregnancies. Semin Fetal Neonatal Med Elsevier Ltd. 2010;15:336–41.
Gardella C, Riley DE, Hitti J, Agnew K, Krieger JN, Eschenbach DA. Identification and sequencing of bacterial rDNAs in culture-negative amniotic fluid from women in premature labor. Am J Perinatol. 2004;21:319–23.
Krohn MA, Hillier SL, Nugent RP, Cotch MF, Carey JC, Gibbs RS, et al. The genital flora of women with intraamniotic infection. J Infect Dis. 1995;171:1475–80.
Hillier SL, Martius J, Krohn M, Kiviat N, Holmes KK, Eschenbach DA. A case–control study of chorioamnionic infection and histologic chorioamnionitis in prematurity. N Engl J Med. 1988;319:972–8.
Romero R, Sirtori M, Oyarzun E, Avila C, Mazor M, Callahan R, et al. Infection and labor V. Prevalence, microbiology, and clinical significance of intraamniotic infection in women with preterm labor and intact membranes. Am J Obstet Gynecol. 1989;161:817–24.
Watts HD, Krohn MA, Hillier SL, Eschenbach DA. The association of occult amniotic fluid infection with gestational age and neonatal outcome among women in preterm labor. Obstet Gynecol. 1992;79:351–7.
Yoon BH, Romero R, Bin MJ, Shim S-S, Kim M, Kim G, et al. Clinical significance of intra-amniotic inflammation in patients with preterm labor and intact membranes. Am J Obstet Gynecol. 2001;185:1130–6.
DiGiulio DB, Romero R, Amogan HP, Kusanovic JP, Bik EM, Gotsch F, et al. Microbial prevalence, diversity and abundance in amniotic fluid during preterm labor: a molecular and culture-based investigation. PLoS One. 2008;3:e3056.
Hillier SL, Witkin SS, Krohn MA, Watts DH, Kiviat NB, Eschenbach DA. The relationship of amniotic fluid cytokines and preterm delivery, amniotic fluid infection, histologic chorioamnionitis, and chorioamnion infection. Obstet Gynecol. 1993;81:941–8.
Jacobsson B, Mattsby-Baltzer I, Hagberg H. Interleukin-6 and interleukin-8 in cervical and amniotic fluid: relationship to microbial invasion of the chorioamniotic membranes. BJOG An Int. J. Obstet. Gynaecol. 2005;112:719–24.
Nugent RP, Krohn MA, Hillier SL. Reability of diagnosing bacterial vaginosis is improved by standardized method of gram stain interpretation. J Clin Microbiol. 1991;29:297–301.
Hill GB. The microbiology of bacterial vaginosis. Am J Obstet Gynecol. 1993;169:450–4.
Romero R, Hassan SS, Gajer P, Tarca AL, Fadrosh DW, Nikita L, et al. The composition and stability of the vaginal microbiota of normal pregnant women is different from that of non-pregnant women. Microbiome. 2014;2:1–19.
Freitas AC, Chaban B, Bocking A, Rocco M, Yang S, Hill JE, et al. The vaginal microbiome of pregnant women is less rich and diverse, with lower prevalence of Mollicutes, compared to non-pregnant women. Sci Rep. 2017;7:1–16.
Aagaard KM, Riehle K, Ma J, Segata N, Mistretta TA, Coarfa C, et al. A metagenomic approach to characterization of the vaginal microbiome signature in pregnancy. PLoS One. 2012;7:e36466.
Walther-António MRS, Jeraldo P, Berg Miller ME, Yeoman CJ, Nelson KE, Wilson BA, et al. Pregnancy’s stronghold on the vaginal microbiome. PLoS One. 2014;9:e98514.
MacIntyre DA, Chandiramani M, Lee YS, Kindinger L, Smith A, Angelopoulos N, et al. The vaginal microbiome during pregnancy and the postpartum period in a European population. Sci Rep. 2015;5:1–9.
Leitich H, Bodner-Adler B, Brunbauer M, Kaider A, Egarter C, Husslein P. Bacterial vaginosis as a risk factor for preterm delivery: a meta-analysis. Am J Obstet Gynecol. 2003;189:139–47.
Romero R, Hassan SS, Gajer P, Tarca AL, Fadrosh DW, Bieda J, et al. The vaginal microbiota of pregnant women who subsequently have spontaneous preterm labor and delivery and those with a normal delivery at term. Microbiome. 2014;2:1–15.
Hyman RW, Fukushima M, Jiang H, Fung E, Rand L, Johnson B, et al. Diversity of the vaginal microbiome correlates with preterm birth. Reprod Sci. 2014;21:32–40.
DiGiulio DB, Callahan BJ, Mcmurdie PJ, Costello EK, Lyell DJ, Robaczewska A, et al. Temporal and spatial variation of the human microbiota during pregnancy. Proc Natl Acad Sci U S A. 2015;112:11060–5.
Callahan BJ, DiGiulio DB, Goltsman DSA, Sun CL, Costello EK, Jeganathan P, et al. Replication and refinement of a vaginal microbial signature of preterm birth in two racially distinct cohorts of US women. Proc Natl Acad Sci U S A. 2017;114:9966–71.
Links MG, Dumonceaux TJ, Hemmingsen SM, Hill JE. The chaperonin-60 universal target is a barcode for bacteria that enables de novo assembly of metagenomic sequence data. PLoS One. 2012;7:e49755.
Paramel Jayaprakash T, Schellenberg JJ, Hill JE. Resolution and characterization of distinct cpn60-based subgroups of Gardnerella vaginalis in the vaginal microbiota. PLoS One. 2012;7:e43009.
Harris PA, Taylor R, Thielke R, Payne J, Gonzalez N, Conde JG. Research electronic data capture (REDCap)—a metadata-driven methodology and workflow process for providing translational research informatics support. J Biomed Inform Elsevier Inc. 2009;42:377–81.
Lee DH, Zo YG, Kim SJ. Nonradioactive method to study genetic profiles of natural bacterial communities by PCR—single-strand-conformation polymorphism. Appl Environ Microbiol. 1996;62:3112–20.
Chaban B, Albert AYK, Links MG, Gardy J, Tang P, Hill JE. Characterization of the upper respiratory tract microbiomes of patients with pandemic H1N1 influenza. PLoS One. 2013;8:e69559.
Hill JE, Penny SL, Crowell KG, Goh SH, Hemmingsen SM. cpnDB: a chaperonin sequence database. Genome Res. 2004;14:1669–75.
van Kuppeveld FJM, van der Logt JTM, Angulo AF, van Zoest MJ, Quint WGV, Niesters HGM, et al. Genus- and species-specific identification of mycoplasmas by 16S rRNA amplification. Appl Environ Microbiol. 1992;58:2606–15.
Watson HL, Blalock DK, Cassell GH. Variable antigens of Ureaplasma urealyticum containing both serovar-specific and serovar-cross-reactive epitopes. Infect Immun. 1990;58:3679–88.
Chaban B, Links MG, Jayaprakash TP, Wagner EC, Bourque DK, Lohn Z, et al. Characterization of the vaginal microbiota of healthy Canadian women through the menstrual cycle. Microbiome. 2014;2:1–12.
Hill JE, Town JR, Hemmingsen SM. Improved template representation in cpn60 polymerase chain reaction (PCR) product libraries generated from complex templates by application of a specific mixture of PCR primers. Environ Microbiol. 2006;8:741–6.
Schellenberg JJ, Links MG, Hill JE, Dumonceaux TJ, Kimani J, Jaoko W, et al. Molecular definition of vaginal microbiota in East African commercial sex workers. Appl Environ Microbiol. 2011;77:4066–74.
Links MG, Chaban B, Hemmingsen SM, Muirhead K, Hill JE. mPUMA: a computational approach to microbiota analysis by de novo assembly of operational taxonomic units based on protein-coding barcode sequences. Microbiome. 2013;1:1–7.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–52.
Schellenberg JJ, Links MG, Hill JE, Dumonceaux TJ, Peters GA, Tyler SD, et al. Pyrosequencing of the chaperonin-60 universal target as a tool for determining microbial community composition. Appl Environ Microbiol. 2009;75:2889–98.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc. 1995;57:289–300.
Caporaso JG, Kuczynski J, Stombaugh JI, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.
Oksanen AJ, Blanchet FG, Kindt R, Legendre P, Minchin PR, Hara RBO, et al. Community Ecology Package “vegan”. R version 2.3-1. Available: https://cran.r-project.org/web/packages/vegan/. 2015.
Endres DM, Schindelin JE. A new metric for probability distributions. IEEE Trans Inf Theory. 2003;49:1858–60.
Fernandes AD, Macklaim JM, Linn TG, Reid G, Gloor GB. ANOVA-like differential expression (ALDEx) analysis for mixed population RNA-Seq. PLoS One. 2013;8:e67019.
Ravel J, Gajer P, Abdo Z, Schneider GM, Koenig SSK, McCulle SL, et al. Vaginal microbiome of reproductive-age women. Proc Natl Acad Sci U S A. 2011;108:4680–7.
Albert AYK, Chaban B, Wagner EC, Schellenberg JJ, Links MG, van Schalkwyk J, et al. A study of the vaginal microbiome in healthy Canadian women utilizing cpn60-based molecular profiling reveals distinct Gardnerella subgroup community state types. PLoS One. 2015;10:e0135620.
Nasioudis D, Forney LJ, Schneider GM, Gliniewicz K, France M, Boester A, et al. Influence of pregnancy history on the vaginal microbiome of pregnant women in their first trimester. Sci. Rep. Springer US; 2017;7:1–6.
Mercer BM, Goldenberg RL, Moawad AH, Meis PJ, Iams JD, Das AF, et al. The preterm prediction study: effect of gestational age and cause of preterm birth on subsequent obstetric outcome. Am J Obstet Gynecol. 1999;181:1216–21.
Goldenberg RL, Hauth JC, Andrews WW. Intrauterine infection and preterm delivery. N Engl J Med. 2000;342:1500–7.
Kindinger LM, Bennett PR, Lee YS, Marchesi JR, Smith A, Cacciatore S, et al. The interaction between vaginal microbiota, cervical length, and vaginal progesterone treatment for preterm birth risk. Microbiome. 2017;5:1–14.
Kekki M, Kurki T, Pelkonen J, Kukinen-Raty M, Cacciatore B, Paavonen J. Vaginal clindamycin in preventing preterm birth and peripartal infections in asymptomatic women with bacterial vaginosis: a randomized, controlled trial. Obstet Gynecol. 2001;97:643–8.
Andrews WW, Goldenberg RL, Hauth JC, Cliver SP, Copper R, Conner M. Interconceptional antibiotics to prevent spontaneous preterm birth: a randomized clinical trial. Am J Obstet Gynecol. 2006;194:617–23.
Morency A, Bujold E. The effect of second-trimester antibiotic therapy on the rate of preterm birth. J Obstet Gynaecol Canada. 2007;29:35–44.
Okun N, Gronau KA, Hannah ME. Antibiotics for bacterial vaginosis or Trichomonas vaginalis in pregnancy: a systematic review. Obstet Gynecol. 2005;105:857–68.
Thinkhamrop J, Hofmeyr GJ, Adetoro O, Lumbiganon P, Ota E. Antibiotic prophylaxis during the second and third trimester to reduce adverse pregnancy outcomes and morbidity. Cochrane Database Syst Rev. 2015;1:CD002250.
Beigi RH, Austin MN, Meyn LA, Krohn MA, Hillier SL. Antimicrobial resistance associated with the treatment of bacterial vaginosis. Am J Obstet Gynecol. 2004;191:1124–9.
Budd W, Bostwick DG, Woody J, Hunt C. Antimicrobial resistance genes and modelling of treatment failure in bacterial vaginosis: clinical study of 289 symptomatic women. J Med Microbiol. 2016;65:377–86.
Brown RG, Marchesi JR, Lee YS, Smith A, Lehne B, Kindinger LM, Terzidou V, Holmes E, Nicholson JK, Bennett PR, et al. Vaginal dysbiosis increases risk of preterm fetal membrane rupture, neonatal sepsis and is exacerbated by erythromycin. BMC Med. 2018;16(1):9.
Paramel Jayaprakash T, Wagner EC, van Schalkwyk J, Albert AYK, Hill JE, Money DM. High diversity and variability in the vaginal microbiome in women following preterm premature rupture of membranes (PPROM): a prospective cohort study. PLoS One. 2016;11:e0166794.
Ling Z, Kong J, Liu F, Zhu H, Chen X, Wang Y, et al. Molecular analysis of the diversity of vaginal microbiota associated with bacterial vaginosis. BMC Genomics. 2010;11:1–16.
Foxman B, Wen A, Srinivasan U, Goldberg D, Marrs CF, Owen J, et al. Mycoplasma, bacterial vaginosis-associated bacteria BVAB3, race, and risk of preterm birth in a high-risk cohort. Am J Obstet Gynecol. 2014;210:226.e1–7.
Hitti J, Garcia P, Totten P, Paul K, Astete S, Holmes KK. Correlates of cervical Mycoplasma genitalium and risk of preterm birth among Peruvian women. Sex Transm Dis. 2010;37:81–5.
Edwards RK, Ferguson RJ, Reyes L, Brown M, Theriaque DW, Duff P. Assessing the relationship between preterm delivery and various microorganisms recovered from the lower genital tract. J. Matern. Neonatal Med. 2006;19:357–63.
Lis R, Rowhani-Rahbar A, Manhart LE. Mycoplasma genitalium infection and female reproductive tract disease: a meta-analysis. Clin Infect Dis. 2015;61:418–26.
Donders GG, Van Calsteren K, Bellen G, Reybrouck R, Van Den Bosch T, Riphagen I, et al. Predictive value for preterm birth of abnormal vaginal flora, bacterial vaginosis and aerobic vaginitis during the first trimester of pregnancy. BJOG An Int J Obstet Gynaecol. 2009;116:1315–24.
Kwak D-W, Hwang H-S, Kwon J-Y, Park Y-W, Kim Y-H. Co-infection with vaginal Ureaplasma urealyticum and Mycoplasma hominis increases adverse pregnancy outcomes in patients with preterm labor or preterm premature rupture of membranes. J Matern Neonatal Med. 2014;27:333–7.
Hillier SL, Nugent RP, Eschenbach DA, Krohn MA, Gibbs RS, Martin DH, et al. Association between bacterial vaginosis and preterm delivery of a low-birth-weight infant. N Engl J Med. 1995;333:1737–42.
Lamont RF, Taylor-Robinson D, Wigglesworth JS, Furr PM, Evans RT, Elder MG. The role of mycoplasmas, ureaplasmas and chlamydiae in the genital tract of women presenting in spontaneous early preterm labour. J Med Microbiol. 1987;24:253–7.
Zhou X, Brown CJ, Abdo Z, Davis CC, Hansmann MA, Joyce P, et al. Differences in the composition of vaginal microbial communities found in healthy Caucasian and black women. ISME J. 2007;112:121–33.
Fettweis JM, Brooks JP, Serrano MG, Sheth NU, Girerd PH, Edwards DJ, et al. Differences in vaginal microbiome in African American women versus women of European ancestry. Microbiology. 2014;160:2272–82.
The authors are grateful to the women who participated in the study and acknowledge the contribution and support of Ontario Birth Study Team members and the participants. The VOGUE Research Group is Deborah Money, Alan Bocking, Sean Hemmingsen, Janet Hill, Gregor Reid, Tim Dumonceaux, Gregory Gloor, Matthew Links, Kieran O’Doherty, Patrick Tang, Julianne van Schalkwyk and Mark Yudin.
Financial support was provided by a joint Canadian Institutes of Health Research (CIHR) Emerging Team Grant and a Genome British Columbia (GBC) grant awarded to DMM, AB and JEH (grant reference #108030) as well as a CIHR grant MOP-82799 to AB. ACF was supported by a University of Saskatchewan graduate scholarship. Funding for the Ontario Birth Study was provided by Mount Sinai Hospital Foundation, Lunenfeld-Tanenbaum Research Institute and the Department of Obstetrics and Gynecology at Mount Sinai Hospital.
Availability of data and materials
The dataset supporting the results of this article is available in the NCBI SRA repository (Accession SRP073152, BioProject PRJNA317763; BioProject PRJNA403856).
Ethics approval and consent to participate
This study received ethical approval from the University of British Columbia Children’s and Women’s Research Ethics Board (Approval Number H14-01954), and Mount Sinai Hospital Research Ethics Board (Approval Number 15-0184-E). All participants provided written consent at enrolment.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.