Skip to main content

The interrelation between microbial immunoglobulin coating, vaginal microbiota, ethnicity, and preterm birth

Abstract

Background

Vaginal microbiota composition is associated with spontaneous preterm birth (sPTB), depending on ethnicity. Host-microbiota interactions are thought to play an important underlying role in this association between ethnicity, vaginal microbiota and sPTB.

Methods

In a prospective cohort of nulliparous pregnant women, we assessed vaginal microbiota composition, vaginal immunoglobulins (Igs), and local inflammatory markers. We performed a nested case–control study with 19 sPTB cases, matched based on ethnicity and midwifery practice to 19 term controls.

Results

Of the 294 included participants, 23 pregnancies ended in sPTB. We demonstrated that Lactobacillus iners-dominated microbiota, diverse microbiota, and ethnicity were all independently associated with sPTB. Microbial Ig coating was associated with both microbiota composition and ethnicity, but a direct association with sPTB was lacking. Microbial IgA and IgG coating were lowest in diverse microbiota, especially in women of any ethnic minority. When correcting for microbiota composition, increased microbial Ig coating correlated with increased inflammation.

Conclusion

In these nulliparous pregnant women, vaginal microbiota composition is strongly associated with sPTB. Our results support that vaginal mucosal Igs might play a pivotal role in microbiota composition, microbiota-related inflammation, and vaginal community disparity within and between ethnicities. This study provides insight in host-microbe interaction, suggesting that vaginal mucosal Igs play an immunomodulatory role similar to that in the intestinal tract.

Video Abstract

Background

An estimated 15 million babies are born preterm each year worldwide [1]. Preterm birth (PTB) is defined as birth before 37 completed weeks of gestation and is a major cause for perinatal mortality and neonatal morbidity [2]. Currently, the prevalence ranges from 5 to 18% across various countries [3]. PTB is usually specified based on its onset, which is either a spontaneous onset of labor (sPTB) or induction of labor or primary caesarean section for maternal or fetal indications [2, 4]. The etiology of sPTB is multifactorial and remains poorly understood. The most important risk factor is a previous sPTB, but most sPTB occur in nulliparous women who lack obstetric history [5]. Several other risk factors have been identified including maternal characteristics such as ethnicity, socio-economic status (SES), body mass index, and maternal smoking, as well as characteristics of the pregnancy like fetal sex, a short mid-trimester cervical length, and intra-amniotic infection [6,7,8,9].

Vaginal microbiota play an important role during pregnancy as its composition and dynamics are hypothesized to have an association with sPTB [10,11,12,13]. Vaginal microbiota are either of low diversity, mainly dominated by a single Lactobacillus species, or consist of a diverse range of (facultative) anaerobic bacteria [14]. Higher diversity of the vaginal microbiota is related to bacterial vaginosis, a disease characterized by increased vaginal discharge and increased susceptibly for invading pathogens [15, 16]. Infection and accompanying inflammatory responses are important risk factors for preterm labor or preterm prelabour rupture of membranes, and about 40% of sPTB is associated with infection [17,18,19].

Ethnicity is significantly associated with the vaginal microbiota composition. While both low and high diversity vaginal microbiota are found in women of all ethnicities, Lactobacillus crispatus-dominant vaginal microbiota is present more often in White European women while diverse vaginal microbiota is present more often in women with a sub-Saharan African descent [20]. The etiology of this association is thus far unknown.

Immunoglobulins (Igs) in the mucosal tissue of the female genital tract are key mediators of mucosal immunology and are important in the defense against infections in the reproductive tract [21]. IgA, the predominant antibody in the intestinal tract, can influence gut microbiota composition [22,23,24]. Deviations in Ig coating of intestinal bacteria have been associated with inflammatory bowel diseases [25, 26]. In the vaginal tract, there is more IgG than IgA, in contrast to other mucosal surfaces [27]. The etiology of the differences in abundance and type of Igs between mucosal sites is not well understood. In a previous study, our group demonstrated increased microbial IgA coating of L. crispatus-dominant vaginal microbiota [28]. Ig coating of vaginal microbiota might play a role in the disparity in microbial community composition between ethnicities and might be associated with gynecological and obstetric diseases.

In this prospective cohort study, we collected vaginal swabs of nulliparous healthy pregnant women at antenatal booking in the first trimester to investigate vaginal microbiota composition and microbial immunoglobulin coating and studied the associations with ethnicity and sPTB. Furthermore, we performed a nested case–control study matching participants with sPTB to participants with uncomplicated term birth. For this subset, we measured unbound Igs and a broad set of inflammatory cytokines and chemokines in vaginal fluid.

Methods

Study design and participants

For this study, we used data and vaginal swab material from women included in the PROPELLOR cohort. The study protocol and methods are described in a previous publication [29]. In short, the study included nulliparous women ≥ 18 years who received antenatal care at participating midwifery practices in the Netherlands before 24 weeks of gestation and had a low-risk singleton pregnancy at their first visit. For this study, all participants of whom a vaginal swab and pregnancy outcome was available (no loss to follow-up) were included. Nulliparity was defined as never having had a pregnancy progress beyond 16 weeks of gestation [30]. At the first prenatal visit, usually between 8 and 12 weeks pregnancy, women were approached to participate in this study by their midwife. A self-administered vaginal swab (eSwab, Copan Diagnostics Inc., Murietta, USA) was collected at inclusion. Swabs with the original medium were stored at − 20 °C until transfer at a shipping temperature at − 20 °C to the central storage facility at − 80 °C, storage duration 4–6 years. Written informed consent was obtained from each participant. The study received ethical clearance through the institutional review board of the Academic Medical Center in Amsterdam, the Netherlands (registration number NL43414.018.13).

We analyzed the microbiota composition and microbial bound Igs for all participants. Due to limited funding, we did a nested case control selection of subjects to assess the role of additional inflammatory markers. The additional inflammatory markers included unbound immunoglobulins, cytokines, chemokines, and anti-microbial proteins. Because of financial constraints and missing meta-data, 19 out of 23 women with sPTB were matched based on their ethnicity and midwifery practice to 19 women with term birth. If there was not an appropriate control available within the same midwifery practice, a control was chosen from a midwifery practice in a similar socio-economic region.

Definitions

The primary outcome measure was sPTB, defined as spontaneous onset of labor or spontaneous preterm prelabor rupture of membranes. We differentiated between sPTB between 23 and 37 weeks of gestation and late second trimester loss between 16 and 22 weeks of gestation. Ethnicity was based on participant self-identification. SES was based on status scores from the Netherlands Institute for Social Research [31].

Vaginal microbiota and bioinformatics analysis

We pre-treated the vaginal swabs with lysozyme, mutanolysin (Sigma Aldrich, St. Louis, USA), lysostaphin (AMBI, New York, USA), proteinase K, and RNAse A (Thermo Fisher, Waltham, USA). DNA was extracted using the NucliSENS EasyMAG platform, according to manufacturer protocol (BioMérieux, Marcy l’Etoile, France). We used dual indexed universal primers (319F and 806R) for PCR amplification of the V3–V4 regions of the 16S rRNA genes, as described by Fadrosh et al. [32]. PCR products were normalized and pooled. We purified the samples with Agencourt AMPure XP magnetic beads (BeckmanCoulter, Fullerton, USA). Paired-end sequencing was performed on the Illumina MiSeq platform, according to manufacturer protocol (Illumina, San Diego, USA).

We de-multiplexed raw sequences and removed adapters, barcodes, and heterogeneity spacers using Cutadapt 3.5 [33]. Processing of de-multiplexed sequence data and taxonomic classification was performed using the software QIIME 2 version 2021.8 [34]. Forward and reverse reads were trimmed to 260 and 210 basepairs respectively based on a visible drop in average quality beyond these points (visualization performed with https://view.qiime2.org/). We included two unsampled swabs as negative controls. We deemed samples with a read count < 200 reads too similar to the controls and excluded these from analysis. Amplicon sequence variants (ASVs) were generated with DADA2 [35]. We used a pre-trained Naive Bayes classifier of the Silva v138 reference database to assign genus and species (if possible) to each ASV [36]. Because the Silva v138 reference database does not include Lactobacillus crispatus, one of the most abundant vaginal species, all Lactobacillus sequences without species assignment were further refined using the Nucleotide BLAST (BLASTn) function on the National Center for Biotechnology Information NCBI website. We manually identified sequences belonging to Candidatus Lachnocurva vaginae (formerly Bacterial Vaginosis Associated Bacterium (BVAB) 1), BVAB2 and TM7-H1. The data were not normalized or rarefied. The code for bioinformatics analyses is available in Supplementary material. Raw sequence reads were dehumanized for publication reasons. We grouped samples based on microbiota profile using VALENCIA centroid classification tool [37]. This tool divides vaginal microbiota into community state types (CSTs) based on their taxonomic composition, by calculating their similarity to a set of reference centroids. It identifies seven CSTS, 4 dominated by lactobacilli (CST I by L. crispatus, CST II by L. gasseri, CST III by L. iners, CST V by L. jensenii), and three depleted of lactobacilli (CST IV-A with majority Candidatus Lachnocurva vaginae and Gardnerella vaginalis, CST IV-B with majority G. vaginalis and Atopobium vaginae, CST IV-C with low abundance of G. vaginalis and Candidatus Lachnocurva vaginae). For statistical power, we reduced these to 4 groups, combining CST II and CST V into one group and combining CST IV-A, CST IV-B, and CST IV-C into one group.

Microbial immunoglobulin coating and inflammatory markers

We determined microbial Ig coating as described in a previous publication using flow cytometry [28]. We calculated coating index by multiplying the percentage of bacteria with bound immunoglobulin with the median fluorescence intensity (MFI). We determined total IgA, IgA1, IgA2, secretory IgA (SIgA), and IgG levels in vaginal swabs by enzyme-linked immunosorbent, as described in a previous publication [28]. We measured vaginal cytokines, chemokines, and anti-microbial peptides with a Multiplex assay using a Bio-Plex 200 according to the manufacturer’s instructions, and human beta defensin-2 (HBD-2) in a separate assay according to manufacturer’s protocol with minor adaptions. We measured total vaginal protein concentration to correct for inter-participant variation according to manufacturer’s protocol. Additional methods and manufacturers can be found in Supplementary methods.

Statistical analysis

If necessary, data were log-transformed to derive normal distributions. Missing data were handled using multiple imputation creating 10 imputation datasets, except for the main determinants’ vaginal microbiota profile and immunoglobulin coating. All variables, including microbiota composition, microbial bound Igs, and midwifery practice, were used for imputation. Numerical results are based on pooled estimates over 10 imputation sets using Rubin’s rules [38].

For analyses on ethnicity, the variable was dichotomised between White European and non-White European. SES was divided into low and middle/high. We performed Firth’s correction logistic regression analysis for sPTB, calculating ORs and 95% CIs. We defined a priori which associations to estimate and confounders to use and, due to the limited number of events, accounted for the single most important confounder in multivariable logistic regression. Even so, overfitting (overestimations of associations due to small number of events) could be an issue. To reduce overfitting because of the low incidence of sPTB, we applied Firth’s correction for all ORs. Firth’s correction uses penalized likelihood which aims to shrink estimated associations that are overly optimistic [39]. For the association between individual taxa, the linear discriminant analysis effect size (LEfSe) algorithm was used [40]. For this analysis, only bacterial taxa were included with > 1% of total read count. This algorithm calculates the median relative abundance of all taxa and compared this between participants with and without sPTB. It uses factorial Kruskal–Wallis rank-sum test to detect differential abundances of bacterial taxa between these groups. The estimated effect size of the differentially abundant taxa was calculated using linear discriminant analysis, with a minimum threshold of 2.0. For the remaining analyses of the entire cohort, we used Student’s t-test and ANOVA with post hoc Bonferroni correction for continuous variables and the Chi-square test for categorical variables. For the nested case–control study, we calculated standardized β-coefficient by linear regression with adjustment for microbiota composition and post hoc Bonferroni correction.

Statistical analyses were performed using IBM SPSS statistics (version 28) and R version 3.3.2 (R Core Team (2016)) with the mice, miceadds, and logistf packages. A p-value of < 0.05 was considered statistically significant. Data were visualized using GraphPad Prism (version 9).

Role of the funding source

The study sponsors had no role in the study design, collection, analysis, and interpretation of the data, the writing of the report, and decision to submit this paper for publication.

Results

Study population

A total of 294 participants were included in this study. Key demographics are shown in Table 1. Of the participants, 189 (73.8%) self-identified as White European and 92 (31.3%) had a low SES. sPTB occurred in 23 (7.8%) participants, of which 18 (6.1%) between 23 and 37 weeks of gestation and in five (1.7%) between 16 and 22 weeks of gestation. A miscarriage < 16 weeks of gestation occurred in one (0.3%) participant and in eight (2.7%) participants birth was induced < 37 weeks of gestation for maternal or fetal indications. Key demographics of participants included in the nested case–control study are shown in Table S1.

Table 1 Cohort characteristics (N=154). SD: standard deviation. IQR: interquartile range. *including miscarriages < 16+0 weeks of gestation and pregnancy terminations. **including miscarriages < 16+0 weeks of gestation and pregnancy terminations and participants with cervical cerclage or pessary in situ before start routine cervical length measurements

Microbiota

Of the 294 samples, three samples were excluded because of a read count below the threshold. The remaining samples had an average of 28,195 reads per sample. 16 s rDNA sequence analysis identified community state types (CSTs), grouped into four clusters: dominated by L. crispatus (CST I, n = 139), dominated by L. gasseri or L. jensenii (CST II/V, n = 19), dominated by L. iners (CST III, n = 70), and diverse microbiota (CST IV, n = 63) (Fig. 1A and Table S2). L. iners-dominated (CST III) and diverse microbiota communities (CST IV) were the most common vaginal microbiota communities found in women experiencing sPTB, and logistic regression revealed an increased odds ratio (OR) for L. iners-dominated (CST III) and diverse microbiota (CST IV) compared to women with microbiota dominated by L. crispatus (CST I) (OR 5.2, 95% confidence interval (CI) 1.6–16.5 and OR 5.2, 95% CI 1.6–16.9, respectively, Fig. 1B and Table 2). No sPTB occurred in participants with L. gasseri/L. jensenii-dominated microbiota (CST II/V), resulting in an inaccurate OR. With linear discriminant analysis effect size (LEfSe), the individual taxa L. iners, Finegoldia, and Prevotella amnii were associated with sPTB (Fig. 1C).

Fig. 1
figure 1

Vaginal microbiota and their association with spontaneous preterm birth (sPTB). A Distribution of vaginal community state types (CSTs) of all participants. Numbers represent total participants per group. B Distribution of community state types of participants with and without sPTB. Numbers represent total participants per group. C The linear discriminant analysis (LDA) score, calculated with linear discriminant analysis effect size (LEfSe) algorithm of the association of individual taxa with sPTB. The bar represents the effect size of the taxa associated with sPTB

Table 2 Associations with spontaneous preterm birth. All logistic regression analyses with Firth’ss correction adjusted for gestational age. Only including participants with a pregnancy that progressed beyond 16 weeks of gestation. a Multivariate analyses are adjusted for ethnicity (White European/non-White European), body mass index (continuous) and smoking (no/yes or stopped during pregnancy). Statistically significant data (p<0.05) are presented in bold

Bacteria bound immunoglobulins

We measured immunoglobulin coating levels using flow cytometry, calculated coating index of IgA and IgG, and divided these in quartiles. IgA and IgG coating indices were not associated with sPTB (Table 2). However, immunoglobulin coating was associated with microbiota composition (IgA p < 0.001, IgG p < 0.001) (Fig. 2). IgA and IgG coating index was statistically significantly lower in diverse microbiota (CST IV) compared to L. crispatus (CST I) and L. iners (CST III)-dominated microbiota (IgA L. crispatus/CST I p < 0.001 and L. iners/CST III p = 0.005, IgG L. crispatus/CST I p < 0.001 and L. iners/CST III p < 0.001). Two samples with the lowest IgA and IgG coating had both diverse microbiota composition (CST IV), as depicted in Fig. 2. In one sample, G. vaginalis was the predominant species and the other was dominated by Enterobacteriaceae. One of the samples with low IgA and IgG coating, with G. vaginalis as predominant species, also had lower than average read count (323 reads). Analysis without this sample showed attenuated, but still statistically significant results for IgA coating index (diverse microbiota/CST IV compared to L. crispatus/CST I p < 0.001 and L. iners/CST III p = 0.013). For IgG coating, the additional analyses showed the same results as analysis with all participants. We further investigated these results dividing the samples into sub-CSTs (Figure S1, Table S3). Within the group of diverse microbiota, CST IV-A with high to moderate levels of Candidatus Lachnocurva vaginae and G. vaginalis had IgA and IgG levels similar to Lactobacillus-dominated CSTs. CST IV-C2 dominated by Enterococcus spp. showed the lowest IgA and IgG coating. Statistical tests were not possible due to small groups.

Fig. 2
figure 2

Microbial immunoglobulin coating in different community state types. Bars represent mean with standard deviation. ** p < 0.01, *** p < 0.001

Ethnicity

Self-identified ethnicity was associated with sPTB, with an increased risk of sPTB for non-White European women (OR 3.8, 95% CI 1.5–9.4, Table 2). When combining both ethnicity and microbiota composition in a regression model, having L. iners-dominated (CST III) or diverse vaginal microbiota (CST IV) and having a non-White European ethnicity remained statistically significantly associated with sPTB, with slightly attenuated adjusted ORs (aORs) (aOR 3.9, 95% CI 1.2–12.9, aOR 3.9, 95% CI 1.2–13.2; and aOR 2.6, 95% CI 1.0–6.5 respectively, Table 2). Several baseline and pregnancy characteristics were associated with sPTB (Table 2). After adjusting for ethnicity and microbiota profile, urinary tract infection during pregnancy and vaginal blood loss in the 1st or 2nd trimester remained associated with sPTB (aOR 4.0, 95% CI 1.3–12.9, and aOR 3.2, 95% CI 1.2–8.7, respectively).

Ethnicity was also associated with microbiota composition and microbial IgA coating. White European women most often had L. crispatus-dominated (CST I) microbiota (n = 121, 57.1%) and non-White European women had most often L. iners-dominated (CST III) microbiota (n = 30, 38.0%) (p < 0.001, Fig. 3A). Non-White European participants had lower IgA coating compared to White European women when diverse microbiota was present (p < 0.001, Fig. 3B). Analysis without the sample with low read count showed the same results.

Fig. 3
figure 3

Associations with ethnicity. A Distribution of community state types (CSTs) in White European and non-White European. Numbers represent total participants per group. B Microbial IgA coating in White European and non-White European participants within various vaginal CSTs. Bars represent mean with standard deviation. *** p < 0.001

Cytokines, chemokines, and anti-inflammatory peptides

In the nested case–control study, we determined the association between the measured unbound immunoglobulins, cytokines, chemokines, or peptides and sPTB, while adjusting for vaginal microbiota composition (Table S4). None was statistically significantly associated with sPTB. We investigated whether microbial Ig coating and unbound Igs were associated with inflammation. Inflammatory cytokines and chemokines IL-1α, IL-1β, IL-2, IL-6, IL-8, CCL4, and CCL5 were positively associated with one or more of the microbial bound and unbound Igs (Table 3). Also, microbial bound and unbound Igs showed positive correlation (Table S5).

Table 3 Correlation between immunoglobulins and inflammatory analytes

Discussion

Our study recapitulated well-known associations between ethnicity and vaginal microbiota composition, with non-White European women having more L. iners-dominated (CST III) and diverse microbiota (CST IV), and the association between diverse microbiota and sPTB [11, 14, 41, 42]. The association between L. iners and sPTB is previously described, but also an association in the opposite direction with term birth is described [43,44,45]. Our study adds to the suspicion that L. iners is more foe that friend during pregnancy in nulliparous women [46]. Our study describes that diverse microbiota have decreased microbial IgA and IgG coating compared to Lactobacillus-dominated microbiota, and that non-White European women with diverse microbiota had lower microbial IgA coating compared to White European women with the same vaginal microbiota profile. With this triad of associations, we anticipated finding an association between decreased IgA coating and sPTB, as diverse microbiota profiles and non-White European participants were over-represented in sPTB cases. The absence of this association was therefore remarkable, and could be reconciled by our finding that lower microbial Ig coating was also associated with lower local inflammation (cytokines and chemokines). Taken together, our findings show that vaginal microbiota and the local immunomodulatory properties of immunoglobulins each play a part in the pathophysiology of preterm birth.

Several studies on sPTB and vaginal microbiota showed similar associations with diverse vaginal microbiota, L. iners, and Prevotella spp. [10,11,12, 42, 47,48,49]. Unique in our study is that the study population only comprises nulliparous women with low risk for sPTB. Parity harbors associations with microbiota composition and sPTB. Vaginal L. iners and Gardnerella dominance is associated with previous birth and the risk for recurrent sPTB is 30% [50,51,52]. Therefore, associations between vaginal microbiota and sPTB might be different for nulliparous and multiparous women. This is corroborated by absence of such associations in two recent studies investigating only women at risk for recurrent sPTB [53, 54]. What makes this especially interesting is that the association between vaginal microbiota composition is present in the first trimester. Risk stratification for sPTB early in pregnancy is limited, especially for nulliparous women. Because obstetric history is a strong prognostic factor, prediction models have limited effect in nulliparous women, while this is the largest group at risk [5, 55]. Several treatments can reduce the risk for sPTB, but these are mainly available to women with an increased risk based on obstetric history. Previous studies illustrated the additional value of vaginal microbiota markers to the prediction of sPTB [12, 56, 57]. Our study confirms this, especially identifying a very low risk for sPTB in White European women with L. crispatus-dominated (CST I) vaginal microbiota. While this requires further research, determining vaginal microbiota composition in pregnancy could help to provide treatment to nulliparous women with an increased risk for sPTB.

Another strength is that our cohort is large enough to include samples of various vaginal microbiota profiles. This allowed us to further investigate the association between immunoglobulin coating and microbiota profile, compared to our previous study [28]. Also, our study population is ethnically diverse, reflecting the population in large Dutch urban areas and provided us the possibility to investigate ethnic disparities.

A limitation of our study is that despite the size, there were only few sPTB cases, limiting our statistical power. Another limitation is that we did not sample over time and thus longitudinal research was not possible. Also, we found no differences in pro-inflammatory mediators in women with sPTB compared to women delivered at term.

In contrast to recent studies, that showed increased levels of pro-inflammatory mediators in vaginal fluid of women with sPTB, including IL-1β, IL-2, IL-6, IL-8, eotaxin, CCL4, and CCL5 [12, 53, 58, 59]. The most likely explanation for this discrepancy could be the gestational age at collection of the vaginal fluid. In our study, material was collected in early pregnancy (8–12 weeks), while other studies only found statistical differences in samples collected later in pregnancy (> 20 weeks). In one study with multiple sampling moments, no differences were found in samples collected at a first time point between 12 and 16 weeks of gestation, while an increase in several pro-inflammatory mediators between the first and second time point (between 20 and 24 weeks of gestation) was associated with sPTB [53]. These results suggest that inflammatory markers are increasing during pregnancy and clear deviations are not yet found early in pregnancy.

Our study revealed remarkable results concerning microbial bound and unbound Igs. A previous study from our group demonstrated higher IgA coating in L. crispatus-dominated vaginal microbiota in non-pregnant healthy women. Due to the larger sample size of the current study, we were able to further elucidate the association between microbial Ig coating. We demonstrated that both IgA and IgG coating are increased in not only L. crispatus but all Lactobacillus dominated microbiota. The previous longitudinal study focussed on the changes over time and showed higher microbial IgA and IgG coating during menses. As pregnancy is hormonally very different without regular vaginal bleeding, we were interested in the vaginal microbial immunoglobulin coating during pregnancy. Unfortunately, we could not study this longitudinally during pregnancy in the current study. But we did confirm that microbial IgA and IgG coating is also high during pregnancy. Also, the increased the sample size made it possible to study differences in microbial Ig coating in women from different ethnicities.

The association between high microbial coating and high inflammation seems to contradict the association between diverse microbiota and low microbial coating. In previous studies, diverse Lactobacillus-depleted microbiota are associated with increased inflammation [12, 53]. Therefore, one would expect to find an association between low microbial coating and high inflammation. However, these results are similar to earlier data from the gut, with IgA both associated with healthy, diversified microbiota, and with inflammatory diseases [25, 60,61,62,63]. It has been suggested that low-affinity IgA contributes to healthy gut microbiota, while high-affinity IgA is involved in pathogen clearance [64]. Based on our results, we hypothesize that a similar regulation can take place in the female genital tract. The role of microbial IgG coating remains unclear, as IgG levels are very low in the intestinal tract, and research on vaginal microbial IgG coating especially in relation to Lactobacillus spp. is limited. A recent study demonstrated that unbound vaginal IgG levels were highest in diverse microbiota and increasing IgG levels during pregnancy were associated with sPTB [53]. It remains to be elucidated what the exact role of microbial bound and unbound Igs is in vaginal mucosa and genital tract related health outcomes.

Vaginal microbiota and related sPTB risk is associated with ethnicity [10, 11]. Our results imply that ethnicity is also associated with immunoglobulin levels and microbial immunoglobulin binding in the vaginal mucosa. In serum, differential immunoglobulin levels between ethnicity have been identified in studies performed several decades ago. Increased levels of IgG and IgA have been found in Black compared to White populations [65,66,67]. Also, vaginal cytokine levels in White and Black women have been reported to differ and to be differentially influenced by vaginal microbiota composition [12, 68, 69]. The underlying mechanisms in different immunoglobulin levels and their relation to vaginal microbiota between ethnically diverse women has been understudied and remains open for further investigation both in the circulation and at mucosal surfaces.

Conclusions

In conclusion, while microbial immunoglobulin coating is associated with vaginal microbiota composition and ethnicity, it is not associated with sPTB. We did find a strong association between L. iners-dominated and diverse vaginal microbiota and sPTB in nulliparous women. In addition, we further explored the association between microbial Ig coating and vaginal microbiota composition, showing that diverse vaginal microbiota have lower IgA and IgG coating than Lactobacillus-dominated microbiota. Further research should investigate whether microbial immunoglobulin coating plays a role in maintaining a Lactobacillus-dominated microbiota profile and whether it is involved in the ethnic disparities of vaginal microbiota composition.

Availability of data and materials

Due to enhanced privacy legislation regarding the presence of human DNA sequences in publicly available datasets, we cannot make raw sequencing data as used in the analyses in this study publicly available. For publication purposes, human DNA reads were removed from the sequence files using the software HoCoRT and the human genome assembly GRCh38.p14 as reference [70]. As note, the amount of human DNA detected in each file was lower than 1% and should not affect the outcome of any subsequent analysis. The cleaned sequencing data is available under study accession number PRJEB71956, sample accession numbers ERS17760025 to ERS17760318. The read count of individual ASVs per sample is available in Table S2. The data that support the findings of this study, including the raw sequencing data, are available from the corresponding author on reasonable request. Extensive data and material availability is described in a previous publication  [29].

References

  1. Howson CP, Kinney MV, McDougall L, Lawn JE, Born Too Soon Preterm Birth Action G. Born too soon: preterm birth matters. Reprod Health. 2013;10 Suppl 1:S1.

    Article  PubMed  Google Scholar 

  2. Goldenberg RL, Culhane JF, Iams JD, Romero R. Epidemiology and causes of preterm birth. Lancet. 2008;371(9606):75–84.

    Article  PubMed  PubMed Central  Google Scholar 

  3. 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. 2012;379(9832):2162–72.

    Article  PubMed  Google Scholar 

  4. Stout MJ, Busam R, Macones GA, Tuuli MG. Spontaneous and indicated preterm birth subtypes: interobserver agreement and accuracy of classification. Am J Obstet Gynecol. 2014;211(5):530 e1-4.

    Article  PubMed  Google Scholar 

  5. Schaaf JM, Ravelli AC, Mol BW, Abu-Hanna A. Development of a prognostic model for predicting spontaneous singleton preterm birth. Eur J Obstet Gynecol Reprod Biol. 2012;164(2):150–5.

    Article  PubMed  Google Scholar 

  6. Iams JD, Goldenberg RL, Meis PJ, Mercer BM, Moawad A, Das A, et al. The length of the cervix and the risk of spontaneous premature delivery. National Institute of Child Health and Human Development Maternal Fetal Medicine Unit Network. N Engl J Med. 1996;334(9):567–72.

    Article  CAS  PubMed  Google Scholar 

  7. Liu P, Xu L, Wang Y, Zhang Y, Du Y, Sun Y, Wang Z. Association between perinatal outcomes and maternal pre-pregnancy body mass index. Obes Rev. 2016;17(11):1091–102.

    Article  CAS  PubMed  Google Scholar 

  8. Soneji S, Beltran-Sanchez H. Association of maternal cigarette smoking and smoking cessation with preterm birth. JAMA Netw Open. 2019;2(4):e192514.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Peelen M, Kazemier BM, Ravelli ACJ, de Groot CJM, van der Post JAM, Mol BWJ, et al. Ethnic differences in the impact of male fetal gender on the risk of spontaneous preterm birth. J Perinatol. 2021;41(9):2165–72.

    Article  CAS  PubMed  Google Scholar 

  10. 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(37):9966–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Elovitz MA, Gajer P, Riis V, Brown AG, Humphrys MS, Holm JB, Ravel J. Cervicovaginal microbiota and local immune response modulate the risk of spontaneous preterm delivery. Nat Commun. 2019;10(1):1305.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Fettweis JM, Serrano MG, Brooks JP, Edwards DJ, Girerd PH, Parikh HI, et al. The vaginal microbiome and preterm birth. Nat Med. 2019;25(6):1012–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Peelen MJ, Luef BM, Lamont RF, de Milliano I, Jensen JS, Limpens J, et al. The influence of the vaginal microbiota on preterm birth: a systematic review and recommendations for a minimum dataset for future research. Placenta. 2019;79:30–9.

    Article  PubMed  Google Scholar 

  14. Ravel J, Gajer P, Abdo Z, Schneider GM, Koenig SS, McCulle SL, et al. Vaginal microbiome of reproductive-age women. Proc Natl Acad Sci U S A. 2011;108(Suppl 1):4680–7.

    Article  CAS  PubMed  Google Scholar 

  15. Bayigga L, Kateete DP, Anderson DJ, Sekikubo M, Nakanjako D. Diversity of vaginal microbiota in sub-Saharan Africa and its effects on HIV transmission and prevention. Am J Obstet Gynecol. 2019;220(2):155–66.

    Article  PubMed  Google Scholar 

  16. Wiesenfeld HC, Hillier SL, Krohn MA, Landers DV, Sweet RL. Bacterial vaginosis is a strong predictor of Neisseria gonorrhoeae and Chlamydia trachomatis infection. Clin Infect Dis. 2003;36(5):663–8.

    Article  PubMed  Google Scholar 

  17. Agrawal V, Hirsch E. Intrauterine infection and preterm labor. Semin Fetal Neonatal Med. 2012;17(1):12–9.

    Article  PubMed  Google Scholar 

  18. Green ES, Arck PC. Pathogenesis of preterm birth: bidirectional inflammation in mother and fetus. Semin Immunopathol. 2020;42(4):413–29.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Romero R, Espinoza J, Goncalves LF, Kusanovic JP, Friel LA, Nien JK. Inflammation in preterm and term labour and delivery. Semin Fetal Neonatal Med. 2006;11(5):317–26.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Borgdorff H, van der Veer C, van Houdt R, Alberts CJ, de Vries HJ, Bruisten SM, et al. The association between ethnicity and vaginal microbiota composition in Amsterdam, the Netherlands. PLoS One. 2017;12(7):e0181135.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Chen A, McKinley SA, Wang S, Shi F, Mucha PJ, Forest MG, Lai SK. Transient antibody-mucin interactions produce a dynamic molecular shield against viral invasion. Biophys J. 2014;106(9):2028–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Breedveld A, van Egmond M. IgA and FcαRI: pathological roles and therapeutic opportunities. Front Immunol. 2019;10(553).

  23. Bunker JJ, Erickson SA, Flynn TM, Henry C, Koval JC, Meisel M, et al. Natural polyreactive IgA antibodies coat the intestinal microbiota. Science. 2017;358(6361).

  24. Okai S, Usui F, Yokota S, Hori IY, Hasegawa M, Nakamura T, et al. High-affinity monoclonal IgA regulates gut microbiota and prevents colitis in mice. Nat Microbiol. 2016;1(9):16103.

    Article  CAS  PubMed  Google Scholar 

  25. Palm NW, de Zoete MR, Cullen TW, Barry NA, Stefanowski J, Hao L, et al. Immunoglobulin A coating identifies colitogenic bacteria in inflammatory bowel disease. Cell. 2014;158(5):1000–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Harmsen HJ, Pouwels SD, Funke A, Bos NA, Dijkstra G. Crohn’s disease patients have more IgG-binding fecal bacteria than controls. Clin Vaccine Immunol. 2012;19(4):515–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Usala SJ, Usala FO, Haciski R, Holt JA, Schumacher GF. IgG and IgA content of vaginal fluid during the menstrual cycle. J Reprod Med. 1989;34(4):292–4.

    CAS  PubMed  Google Scholar 

  28. Breedveld AC, Schuster HJ, van Houdt R, Painter RC, Mebius RE, van der Veer C, et al. Enhanced IgA coating of bacteria in women with Lactobacillus crispatus-dominated vaginal microbiota. Microbiome. 2022;10(1):15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Schuster HJ, Peelen M, Hajenius PJ, van Beukering MDM, van Eekelen R, Schonewille M, et al. Risk factors for spontaneous preterm birth among healthy nulliparous pregnant women in the Netherlands, a prospective cohort study. Health Sci Rep. 2022;5(3):e585.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Heijneman MJ, Evers JLH, Massuger LFAG, Steegers EAP. Obstetrie en gynaecologie. De voortplanting van de mens. Maarssen: Elsevier; 2008.

    Google Scholar 

  31. Knol FA. Van hoog naar laag: van laag naar hoog. Sociaal en Cultureel planbureau. 1998.

  32. Fadrosh DW, Ma B, Gajer P, Sengamalay N, Ott S, Brotman RM, Ravel J. An improved dual-indexing approach for multiplexed 16S rRNA gene sequencing on the Illumina MiSeq platform. Microbiome. 2014;2(1):6.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17(1):10–2.

    Article  Google Scholar 

  34. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37(8):852–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJ, Holmes SP. DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13(7):581–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73(16):5261–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. France MT, Ma B, Gajer P, Brown S, Humphrys MS, Holm JB, et al. VALENCIA: a nearest centroid classification method for vaginal microbial communities based on composition. Microbiome. 2020;8(1):166.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Rubin D. Multiple imputation for nonresponse in surveys. New York: Wiley; 2004.

    Google Scholar 

  39. Wang X. Firth logistic regression for rare variant association tests. Front Genet. 2014;5:187.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, Huttenhower C. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12(6):R60.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Chang DH, Shin J, Rhee MS, Park KR, Cho BK, Lee SK, Kim BC. Vaginal microbiota profiles of native Korean women and associations with high-risk pregnancy. J Microbiol Biotechnol. 2020;30(2):248–58.

    Article  CAS  PubMed  Google Scholar 

  42. Tabatabaei N, Eren AM, Barreiro LB, Yotova V, Dumaine A, Allard C, Fraser WD. Vaginal microbiome in early pregnancy and subsequent risk of spontaneous preterm birth: a case-control study. BJOG. 2019;126(3):349–58.

    Article  CAS  PubMed  Google Scholar 

  43. 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):6.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Kumar M, Murugesan S, Singh P, Saadaoui M, Elhag DA, Terranegra A, et al. Vaginal microbiota and cytokine levels predict preterm delivery in Asian women. Front Cell Infect Microbiol. 2021;11:639665.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Payne MS, Newnham JP, Doherty DA, Furfaro LL, Pendal NL, Loh DE, Keelan JA. A specific bacterial DNA signature in the vagina of Australian women in midpregnancy predicts high risk of spontaneous preterm birth (the Predict1000 study). Am J Obstet Gynecol. 2021;224(2):206 e1-e23.

    Article  PubMed  Google Scholar 

  46. Petrova MI, Reid G, Vaneechoutte M, Lebeer S. Lactobacillus iners: Friend or Foe? Trends Microbiol. 2017;25(3):182–91.

    Article  CAS  PubMed  Google Scholar 

  47. Brown RG, Al-Memar M, Marchesi JR, Lee YS, Smith A, Chan D, et al. Establishment of vaginal microbiota composition in early pregnancy and its association with subsequent preterm prelabor rupture of the fetal membranes. Transl Res. 2019;207:30–43.

    Article  PubMed  PubMed Central  Google Scholar 

  48. 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(35):11060–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Freitas AC, Bocking A, Hill JE, Money DM, Group VR. Increased richness and diversity of the vaginal microbiota and spontaneous preterm birth. Microbiome. 2018;6(1):117.

    Article  PubMed  PubMed Central  Google Scholar 

  50. 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. 2017;7(1):10201.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Kervinen K, Holster T, Saqib S, Virtanen S, Stefanovic V, Rahkonen L, et al. Parity and gestational age are associated with vaginal microbiota composition in term and late term pregnancies. 2021.

  52. Phillips C, Velji Z, Hanly C, Metcalfe A. Risk of recurrent spontaneous preterm birth: a systematic review and meta-analysis. BMJ Open. 2017;7(6):e015402.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Chan D, Bennett PR, Lee YS, Kundu S, Teoh TG, Adan M, et al. Microbial-driven preterm labour involves crosstalk between the innate and adaptive immune response. Nat Commun. 2022;13(1):975.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Goodfellow L, Verwijs MC, Care A, Sharp A, Ivandic J, Poljak B, et al. Vaginal bacterial load in the second trimester is associated with early preterm birth recurrence: a nested case-control study. BJOG. 2021;128(13):2061–72.

    Article  CAS  PubMed  Google Scholar 

  55. Meertens LJE, van Montfort P, Scheepers HCJ, van Kuijk SMJ, Aardenburg R, Langenveld J, et al. Prediction models for the risk of spontaneous preterm birth based on maternal characteristics: a systematic review and independent external validation. Acta Obstet Gynecol Scand. 2018;97(8):907–20.

    Article  PubMed  Google Scholar 

  56. Flaviani F, Hezelgrave NL, Kanno T, Prosdocimi EM, Chin-Smith E, Ridout AE, et al. Cervicovaginal microbiota and metabolome predict preterm birth risk in an ethnically diverse cohort. JCI Insight. 2021;6(16).

  57. Park S, Oh D, Heo H, Lee G, Kim SM, Ansari A, et al. Prediction of preterm birth based on machine learning using bacterial risk score in cervicovaginal fluid. Am J Reprod Immunol. 2021;86(3):e13435.

    Article  CAS  PubMed  Google Scholar 

  58. Amabebe E, Reynolds S, He X, Wood R, Stern V, Anumba DOC. Infection/inflammation-associated preterm delivery within 14 days of presentation with symptoms of preterm labour: a multivariate predictive model. PLoS ONE. 2019;14(9):e0222455.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Ashford K, Chavan NR, Wiggins AT, Sayre MM, McCubbin A, Critchfield AS, O’Brien J. Comparison of serum and cervical cytokine levels throughout pregnancy between preterm and term births. AJP Rep. 2018;8(2):e113–20.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Huus KE, Bauer KC, Brown EM, Bozorgmehr T, Woodward SE, Serapio-Palacios A, et al. Commensal bacteria modulate immunoglobulin A binding in response to host nutrition. Cell Host Microbe. 2020;27(6):909-21 e5.

    Article  CAS  PubMed  Google Scholar 

  61. Kau AL, Planer JD, Liu J, Rao S, Yatsunenko T, Trehan I, et al. Functional characterization of IgA-targeted bacterial taxa from undernourished Malawian children that produce diet-dependent enteropathy. Sci Transl Med. 2015;7(276):276ra24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Hansen IS, Hoepel W, Zaat SAJ, Baeten DLP, den Dunnen J. Serum IgA immune complexes promote proinflammatory cytokine production by human macrophages, monocytes, and Kupffer cells through FcalphaRI-TLR cross-talk. J Immunol. 2017;199(12):4124–31.

    Article  CAS  PubMed  Google Scholar 

  63. Sterlin D, Fadlallah J, Adams O, Fieschi C, Parizot C, Dorgham K, et al. Human IgA binds a diverse array of commensal bacteria. J Exp Med. 2020;217(3).

  64. Jackson MA, Pearson C, Ilott NE, Huus KE, Hegazy AN, Webber J, et al. Accurate identification and quantification of commensal microbiota bound by host immunoglobulins. Microbiome. 2021;9(1):33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Maddison SE, Stewart CC, Farshy CE, Reimer CB. The relationship of race, sex, and age to concentrations of serum immunoglobulins expressed in international units in healthy adults in the USA. Bull World Health Organ. 1975;52(2):179–85.

    CAS  PubMed  PubMed Central  Google Scholar 

  66. Shulman G, Gilich GC, Andrew MJ. Serum immunoglobulins G, A and M in White and Black adults on the Witwatersrand. S Afr Med J. 1975;49(29):1160–4.

    CAS  PubMed  Google Scholar 

  67. Grundbacher FJ. Heritability estimates and genetic and environmental correlations for the human immunoglobulins G, M, and A. Am J Hum Genet. 1974;26(1):1–12.

    CAS  PubMed  PubMed Central  Google Scholar 

  68. Ryckman KK, Williams SM, Krohn MA, Simhan HN. Racial differences in cervical cytokine concentrations between pregnant women with and without bacterial vaginosis. J Reprod Immunol. 2008;78(2):166–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Lennard K, Dabee S, Barnabas SL, Havyarimana E, Blakney A, Jaumdally SZ, et al. Microbial composition predicts genital tract inflammation and persistent bacterial vaginosis in South African adolescent females. Infect Immun. 2018;86(1).

  70. Rumbavicius I, Rounge TB, Rognes T. HoCoRT: host contamination removal tool. BMC Bioinformatics. 2023;24(1):371.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Funding

This study was financed by ZonMw, the Netherlands Organization of Health Research and Development, project number VICI 91814650, and Amsterdam Reproduction and Development (AR&D 2016). The study sponsors had no role in the study design, collection, analysis, and interpretation of the data, the writing of the report, and decision to submit this paper for publication.

Author information

Authors and Affiliations

Authors

Contributions

M.K., P.J.H., and R.C.P. conceived and led the clinical study. M.v.E. and R.v.H. led immunology and microbiome data generation. H.J.S. and A.C.B. generated immunology data. H.J.S. and R.v.H. generated microbiome data. H.J.S., A.C.B., S.M., and R.v.E. performed data processing and formal analyses. H.J.S., A.C.B., R.C.P, P.H.M.S., M.v.E., and R.v.H. performed data interpretation. H.J.S. and A.C.B. wrote the first draft of the manuscript. All authors critically reviewed, read, and approved the final manuscript.

Corresponding author

Correspondence to H. J. Schuster.

Ethics declarations

Ethics approval and consent to participate

The study received ethical clearance through the institutional review board of the Academic Medical Center in Amsterdam, the Netherlands (registration number NL43414.018.13). All participants provided written informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Supplementary Material 1.

Supplementary methods.

Supplementary Material 2.

Table S1. Characteristics of participants with short cervix (n=136). SD: standard deviation. IQR: interquartile range. *including miscarriage.

Supplementary Material 3.

Table S2. Read count and community state type per sample. 

Supplementary Material 4.

Table S3. Sexually transmitted infections. Multiple samples are available per participant and pathogens can be detected in some, but not all, samples from the participant. Therefore, results are presented as sexually transmitted infections (STIs) per participant in total and per sample individually.

Supplementary Material 5.

Table S4. Associations with spontaneous preterm birth.

Supplementary Material 6.

Table S5. Correlation between microbial bound and unbound immunoglobulins.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Schuster, H.J., Breedveld, A.C., Matamoros, S.P.F. et al. The interrelation between microbial immunoglobulin coating, vaginal microbiota, ethnicity, and preterm birth. Microbiome 12, 99 (2024). https://doi.org/10.1186/s40168-024-01787-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40168-024-01787-z

Keywords