Background: Large datasets for investigating vaginal flora change at frequent, repeated intervals are limited and graphical methods for exploring such data are inadequate. We report 2-year weekly vaginal flora changes based on Gram stain using lasagna plots.
Methods: Weekly vaginal flora patterns were evaluated among 211 sexually experienced women with ≥18 months of follow-up in Rakai, Uganda. Vaginal flora swabs were self-collected weekly and categorized by Nugent Gram stain criteria (0–3, normal; 4–6, intermediate; 7–10, bacterial vaginosis [BV]). Vaginal flora patterns were analyzed as the percentage of weekly observations with BV (longitudinal prevalence) and illustrated by lasagna plots. Characteristics of women were compared across tertiles of longitudinal prevalence of BV.
Results: Ninety-five percent of women had at least 1 episode of BV over 2 years, with one-third of women spending more than half (52%–100%) of their time with BV. Vaginal pH >4.5 increased with increasing tertiles of longitudinal prevalence of BV (P < 0.001). Weekly fluctuation in vaginal flora states, as measured by a change in flora states from the before current visit, was highest in the middle (41.9%) compared with the lower (30.1%) and upper tertiles (27.8%, P < 0.001). HIV status and reported vaginal symptoms did not differ significantly across BV tertiles.
Conclusions: Women exhibited different patterns of vaginal flora changes over time, which could not be described by baseline behaviors. Lasagna plots aided in describing the natural history of BV within and across women and may be applied to future BV natural history studies.
Few factors differed between women with higher versus lower frequencies of bacterial vaginosis. Lasagna plots facilitated the visualization and interpretation of weekly vaginal flora patterns from a large sample of women.
From the *Department of Population, Family, and Reproductive Health, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD; †Epidemiology Branch, Division of Epidemiology, Statistics, and Prevention Research, Eunice Kennedy Shriver National Institute of Child Health and Human Development, Rockville, MD; ‡Makerere University School of Public Health, Kampala, Uganda; §Rakai Health Sciences Program, Kalisizo, Uganda; ¶Department of Biostatistics, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD; and ∥Makerere University College of Health Sciences, Kampala, Uganda
The authors thank Bruce Swihart for his feedback and assistance with lasagna plotting, including providing his program developed in R to generate the figures used in this paper.
Supported in part by the Intramural Research Program of the National Institutes of Health (NIH), Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD), the NIH, NICHD grant R01AI47608 (to M.J.W.); the NIH, National Institute for Allergy and Infectious Diseases (NIAID) grant T32AI050056 (to J.M.Z.), the Cooperative Agreement number R36PS001104 (to M.E.T.) from the Center for Disease Control and Prevention (CDC).
Its contents are solely the responsibility of the authors and do not represent the official views of the CDC.
Correspondence: Marie Thoma, PhD, MHS, Epidemiology Branch/DESPR, Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, 6100 Executive Blvd, Room 7B13E, Rockville, MD 20852. E-mail: firstname.lastname@example.org.
Received for publication February 14, 2011, and accepted June 1, 2011.
Bacterial vaginosis (BV) is the most common vaginal syndrome among reproductive-aged women and is marked by a shift in the normal Lactobacillus-dominant flora to a more diverse flora consisting of commensal anaerobic bacteria.1 Clinically, BV is characterized by vaginal discharge, odor, and an increased vaginal pH (>4.5); however, symptoms are nonspecific and as many as 50% of women remain asymptomatic.2,3 Women with BV are more likely to experience reproductive and obstetric complications, including preterm birth, low birth weight, and sexually transmitted infections.4 Despite its prevalence and adverse sequelae, the etiology and natural history of BV are poorly understood.
The dynamics of vaginal flora variability over time within women add to the complexity of studying this condition. A recent National Institutes of Health workshop on BV research priorities stressed the need for larger studies with more frequent, prospective sampling.1 Longitudinal studies of vaginal flora change with infrequent observations are likely to miss episodes of BV or to detect transient BV, rather than documenting the full variability of vaginal flora patterns. A limited number of studies with frequent, repeated measurements have investigated short-term fluctuations in vaginal microflora at daily or weekly intervals; however, these studies were constrained by small sample sizes and short study duration.5–12 In addition, the display of vaginal flora patterns over time was often limited to a few individuals' representative patterns rather than utilizing the entirety of the data across all individuals. This limitation is further complicated by inadequate graphical methods for large longitudinal datasets, such as spaghetti plots, which can lead to multiple, intersecting lines that can obscure trends and patterns.13 Establishing patterns and adequate display of vaginal flora trajectories within and between women over time is needed from large studies to better understand the natural history of BV.
We report data from a cohort of ever sexually active, reproductive-aged women in Rakai, Uganda, followed weekly for up to 2 years and classify women according to their longitudinal vaginal flora profile. Vaginal flora patterns within and between women are illustrated using lasagna plotting, which is a recently developed exploratory method for representing and summarizing longitudinal, categorical data across a large sample of individuals.13
MATERIALS AND METHODS
A 2-year cohort study was conducted between 2001 and 2003 among 312 consenting females aged 13 to 39 years in rural Rakai District, Uganda. Women who were sexually active and postmenarcheal at baseline (N = 255) and who remained in the study for at least 18 months of observation (N = 211, 82.7% of the eligible sample) were included in this analysis. Ethical approval for this study was provided by institutional review boards (IRBs) in Uganda (the Scientific and Ethics Committee of the Uganda Virus Research Institute) and the United States (the Johns Hopkins Bloomberg School of Public Health IRB and the Columbia University Medical Center IRB).
The study design and data collection has been described previously.14 Briefly, women were enrolled regardless of HIV status, pregnancy, or sexual history. Participants were interviewed in their homes every week for up to 2 years and provided a self-collected vaginal swab for BV and vaginal pH assessment. More detailed questionnaires were administered at baseline, monthly, and 6-monthly visits. Self-collected vaginal swabs were rolled onto slides and air dried, Gram stained, and assessed using the Nugent quantitative morphologic classification for vaginal flora.15 Vaginal flora status based on the Nugent score was defined as normal (0–3 points), intermediate (4–6 points), or bacterial vaginosis (7–10 points). Vaginal pH was determined by BAKER-pHIX pH papers (pH range: 4.0–9.0, Phillipsburg, NJ) affixed to a pediatric tongue depressor and assessed during all nonmenstruating visits.
A serologic sample was analyzed for HIV every 6 months using 2 HIV enzyme immunosorbent assays (Vironostika HIV-1, Organon Teknicka, Charlotte, NC; Cambridge Biotech, Worcester, MA) and Western blot (HIV-1 Western Blot, BioMerieuz-Vitek, St. Louis, MI) for discordant results. InPouch TV culture (BioMed Diagnostics, San Jose, CA) was used for Trichomonas vaginalis detection from vaginal swabs collected every 6 months. At baseline and monthly, urinary hCG was assessed for pregnancy in all women who were not visibly pregnant, had a last menstrual period over 30 days before the visit, were not using Depo Provera, or had lactational amenorrhoea.
Self-reported abnormal genital discharge, dysuria, and ulcers were treated syndromically, using a single, directly observed 2 g dose of metronidazole for discharge, which was one of the recommended regimens for BV during the study period.16 Persistent abnormal discharge was treated with 1 g dose of azithromycin or 400 mg of cefixime for pregnant women. Self-reported treatment was used in the analysis.
For women who completed up to 18 to 24 months of follow-up, we estimated the percentage of follow-up time spent in a given vaginal flora state, which is referred to as “longitudinal prevalence” and is often used to represent the prevalence of common recurrent infections, such as episodes of diarrhea or respiratory infections, in the same individual over time.17 For example, the longitudinal prevalence of BV was estimated by calculating the total number of visits with BV divided by the total number of observations per woman and, subsequently, divided into tertiles. This calculation was also applied to estimate the average frequency (percentage) of time women spent with the following: normal and intermediate vaginal flora states, vaginal pH greater than 4.5, self-reported vaginal symptoms and treatment, menses, pregnancy, weekly sex, new partners in the past month, and substances inserted into the vagina in the past month. The frequency of vaginal flora state transitioning within a woman was estimated by the number of times the prior weekly state differed from the current state divided by the total number of paired, consecutive visits observed.
Vaginal flora patterns within and between women over time were represented by lasagna plots, which is an exploratory graphical method proposed by Swihart et al. for use in large longitudinal epidemiologic studies and analogous to heat maps often used in the genomics literature.13,18 This method was developed to overcome limitations of the spaghetti plot for visualizing longitudinal data, which is subject to overlapping trajectories, indistinguishable data patterns, and difficulties in displaying missing data. In contrast, lasagna plots allow for visualization of transitional, categorical data, and simultaneous illustration of individual and group level (e.g., women grouped according to an external covariate, such as age) information for large groups of individuals. Briefly, the visualization of transitional data can be represented by an m x n history matrix H, where m rows are the number of individuals and n columns are the number of repeated intervals. Thus, each row represents the subject-specific data across n intervals. Graphical visualization of the data is achieved by assigning colors to each element hij based on the vaginal flora status (normal, intermediate, or BV) for that visit. In addition, this framework allows for row, column, and cluster sorting facilitating the graphical display of information and patterns in the data. This method precludes overlapping of trajectories to better reveal trends and patterns.
BV and other vaginal flora states vary over time and are subject to remission and recurrence for different lengths of time across individuals. Two matrices were generated to illustrate these data using lasagna plots. The first matrix (Fig. 1) illustrates the temporal ordering of vaginal flora states and transitions across individuals (Y-axis) over time (X-axis). The second matrix (Fig. 2) sorts vaginal flora states within each individual to display the proportion of follow-up time spent with BV (longitudinal prevalence of BV). Both figures order women by their longitudinal prevalence of BV from lowest to highest and are divided into tertiles of this distribution. Lasagna plots were generated using R 2.8.1 statistical software.19
Baseline characteristics and longitudinal frequencies were compared by tertiles of longitudinal prevalence of BV. Baseline associations were assessed by chi-squared tests for categorical variables and one-way analysis of variance for continuous variables. Medians and interquartile ranges (IQR) of 2-year longitudinal frequencies were compared by Kruskal-Wallis one-way analysis of variance with correction for ties. A 2-sided P value of 0.05 was used to assess statistical significance. Data were analyzed using Stata/SE 10.0 for Windows (Stata Corporation, College Station, TX).
There were a total of 255 ever sexually active, postmenarcheal women at baseline. The mean age of participants at baseline was 24.7 years (SD = 6.4) with 25.7% aged between 14 and 19 years. Eleven percent of them were HIV-positive and almost half of the women (47.6%) had BV at enrollment. Four women seroconverted to HIV by the end of follow-up. Participant retention rates were high with 211 (82.7%) women remaining under observation beyond 18 months and contributing an average of 83.6 (SD = 8.2) weekly visits per woman, which constitutes an 86.2% compliance with the weekly visit schedule over 24 months. Losses to follow-up were 6.3%, 5.9%, and 5.1% by 6 months, 6 to 12 months, and 12 to 18 months, respectively.
Of the 211 women observed for ≥18 months of observation, there were 71 women in the lowest tertile (0%–19.3% longitudinal prevalence of BV), 70 women in the middle tertile (20.2%–52.2% longitudinal prevalence of BV), and 70 in the upper tertile (52.3%–100% longitudinal prevalence of BV). The vaginal flora patterns are depicted by lasagna plots in Figures 1 and 2 and ordered by each woman's longitudinal prevalence of BV from lowest to highest. Figure 1 illustrates temporal transitioning between weekly vaginal flora states. In the lowest tertile, most of the women had persistent normal vaginal flora (light blue) with transient intermediate status (red) or BV (black) episodes. In the upper tertile, most of the women had persistent BV over 2 years with some transitory intermediate episodes and resolution to normal flora. The middle tertile showed the greatest fluctuation between vaginal flora states over time. Figure 1 also shows a cluster of missing observations (off-white boxes) that occurred around 36 and 84 weeks of observation, which correspond to the Christmas and New Years holiday and are similar in timing across all tertiles of longitudinal prevalence. Figure 2 is a within-woman sort on BV status (i.e., collapsed over time) to display the total time each women spent in BV, intermediate and normal flora states. Almost all women had BV for at least one visit during follow-up (n = 202, 95.7%), and women at the bottom of the figure spent almost their entire observation time in a BV state.
Comparisons of the sociodemographic, behavioral, and clinical characteristics by tertiles of longitudinal prevalence of BV are shown in Table 1. Most baseline characteristics did not differ significantly between tertiles of longitudinal prevalence with the exception of baseline vaginal flora status, in which the enrollment prevalence of BV significantly increased with increasing tertiles of longitudinal prevalence of BV (12.7%, 55.1%, and 74.3% for lower, middle, and upper tertiles, respectively, P < 0.001). Differences in mean age at first sex and sex with a circumcised partner were of borderline significance (P = 0.07 and 0.13, respectively).
Transitioning between Nugent score categories (normal, intermediate, or BV) per woman was highest for the middle tertile (median: 41.9%, IQR: 33.3%–50.8%) compared with the lower tertile (median: 30.1%, IQR: 22.0%–39.2%) and upper tertiles (median: 27.8%, IQR: 17.9%–39.0%, P < 0.001). Women in the upper tertile spent the greatest amount of time in a BV state (median: 69.8%, IQR: 61.6–85.5), the middle tertile in an intermediate state (median: 20.0%, IQR: 15.1–30.2), and the lowest tertile in a normal vaginal flora state (median: 76.1%, IQR: 67.9–86.2). A high vaginal pH (pH >4.5) was the only clinical predictor that increased significantly with increasing tertiles of longitudinal prevalence (P < 0.001). Reported frequencies of vaginal discharge or odor and treatment for vaginal symptoms were low and not significantly different across all 3 tertiles. Restricting the analysis to those who experienced any symptom of discharge, odor, or received treatment did not change inference (P = 0.79, 0.28, and 0.91, respectively). There were no differences in the frequency of reported sex in the past week, pregnancy status, substances inserted into the vagina new sex partners, and swabs collected during menstruation across tertiles (data not shown). Finally, the pattern of missing visits seems to be fairly consistent across all individuals, and there were no differences in the average number of observed visits per tertile (lower = 83.7, middle = 83, and upper = 84.2; P = 0.70).
Among women in rural Rakai, Uganda, BV was common with 95% of women having at least 1 weekly episode of BV diagnosed by Nugent Gram stain criteria over a 2-year period. Apart from a woman's starting vaginal flora state, there were few baseline characteristics that predicted differences across tertiles of longitudinal prevalence of BV. Our data also indicate that vaginal symptoms of discharge and odor did not vary across different patterns of vaginal flora, but that the frequency of high vaginal pH (pH >4.5) increased with increasing tertiles of longitudinal prevalence of BV. Lasagna plots further illustrated temporal patterns in vaginal flora and graphically illustrated the duration of onset and remission of BV, transitions between vaginal flora states, and the overall proportion of time spent in a given vaginal flora state within and between women. Incorporating measures of chronicity, such as longitudinal prevalence, into future BV studies may better describe the burden of BV in women.
Women in the upper tertile of longitudinal prevalence of BV tended to persist in a BV state. It has been suggested that women with persistent BV may represent a unique group with distinct vaginal flora that differ from women with less frequent or no BV occurrences.20 It has also been hypothesized that specific bacteria in women with persistent BV by Nugent criteria may maintain a healthy vaginal environment, as defined by lactic acid production and low vaginal pH, despite the lack of lactobacilli.21 However, this assertion was not supported by our data, which shows an increasing percentage of women with a vaginal pH >4.5 associated with increasing tertiles of longitudinal prevalence of BV (Table 1). Women in the middle tertile of longitudinal prevalence of BV had the highest frequency of transitioning between vaginal flora states. Although transient shifts in BV status may be influenced by recent sex or menstruation,14 this association was not reflected across tertiles of longitudinal prevalence of BV. Thus, other factors, such as vaginal bacterial composition, may explain the higher frequency of transitioning in this group over a 2-year period.
There was some evidence to suggest that women who had predominantly normal flora throughout follow-up (lower tertile) were more likely to initiate sex at a later age and have a circumcised partner compared to those with moderate or frequent BV occurrences. Uncircumcised men may increase their partner's exposure to anaerobic bacteria leading to more vaginal flora changes22 and greater BV prevalence.23
There were few differences in clinical factors across tertiles of longitudinal BV prevalence. The prevalence of T. vaginalis was higher in middle and upper tertiles; however, the number of women with T. vaginalis was low and was not statistically significant. Longitudinal studies evaluating BV prevalence in HIV-positive women are conflicting and show either no effect or increased risk in HIV-positive women.24,25 However, both studies supported a protective effect of treatment, higher CD4 cell counts, or lower viral load on BV among HIV-positive women. In this study, there were a total of 21 HIV-positive women who most likely represented women at an earlier stage of HIV infection. Finally, it is unknown why some women develop symptoms of BV and others do not. We did not find that the frequency of self-reported vaginal symptoms varied by greater persistence of high Nugent scores over time. While asymptomatic BV has been associated with adverse health outcomes, there is some debate regarding the utility of Nugent scores for defining a “healthy” vaginal flora environment.26 Nevertheless, the Nugent score method continues to be recommended as the gold standard until new diagnostic methods are developed.1
This study is unique in that it provides a data-rich structure for evaluating the natural history of BV at weekly intervals over 2 years. In addition, we presented methods for illustrating and exploring trends in transitional, categorical data. However, our analyses were limited to use of Nugent score morphology, which does not provide information on bacterial composition and abundance that may influence vaginal flora trajectories. In addition, findings from these rural Ugandan women may not be generalizable to other populations. We previously reported a cross-sectional prevalence of BV in this study population of 47.6%,14 which is similar to the prevalence of non-Hispanic black women in the United States (51.4%).27 Comparisons of the longitudinal prevalence of BV by race in the United States and other populations are warranted.
In summary, BV was common and tended to be persistent in this population. The use of longitudinal prevalence of BV illustrated by lasagna plots provides a new approach to the analysis of vaginal flora in large populations with frequent sampling. Future research should explore the vaginal microbiota using molecular technologies to further elucidate bacterial communities in women with different durations of time spent with BV.
1. Marrazzo JM, Martin DH, Watts DH, et al. Bacterial vaginosis: Identifying research gaps proceedings of a workshop sponsored by DHHS/NIH/NIAID. Sex Transm Dis 2010; 37:732–744.
2. Amsel R, Totten PA, Spiegel CA, et al. Nonspecific vaginitis. Diagnostic criteria and microbial and epidemiologic associations. Am J Med 1983; 74:14–22.
3. Klebanoff MA, Schwebke JR, Zhang J, et al. Vulvovaginal symptoms in women with bacterial vaginosis. Obstet Gynecol 2004; 104:267–272.
4. Koumans EH, Kendrick JS. Preventing adverse sequelae of bacterial vaginosis: A public health program and research agenda. Sex Transm Dis 2001; 28:292–297.
5. Hay PE, Ugwumadu A, Chowns J. Sex, thrush and bacterial vaginosis. Int J STD AIDS 1997; 8:603–608.
6. Keane FE, Ison CA, Taylor-Robinson D. A longitudinal study of the vaginal flora over a menstrual cycle. Int J STD AIDS 1997; 8:489–494.
7. Morison L, Ekpo G, West B, et al. Bacterial vaginosis in relation to menstrual cycle, menstrual protection method, and sexual intercourse in rural Gambian women. Sex Transm Infect 2005; 81:242–247.
8. Priestley CJ, Jones BM, Dhar J, et al. What is normal vaginal flora? Genitourin Med 1997; 73:23–28.
9. Schwebke JR, Morgan SC, Weiss HL. The use of sequential self-obtained vaginal smears for detecting changes in the vaginal flora. Sex Transm Dis 1997; 24:236–239.
10. Schwebke JR, Richey CM, Weiss HL. Correlation of behaviors with microbiological changes in vaginal flora. J Infect Dis
11. Brotman RM, Ravel J, Cone RA, et al. Rapid fluctuation of the vaginal microbiota measured by Gram stain analysis. Sex Transm Infect 2010; 86:297–302.
12. Srinivasan S, Liu C, Mitchell CM, et al. Temporal variability of human vaginal bacteria and relationship with bacterial vaginosis. PLoS One 2010; 5:e10197.
13. Swihart BJ, Caffo B, James BD, et al. Lasagna plots: A saucy alternative to spaghetti plots. Epidemiology 2010; 21:621–625.
14. Thoma ME, Gray RH, Kiwanuka N, et al. The short-term variability of bacterial vaginosis diagnosed by Nugent Gram stain criteria among sexually active women in Rakai, Uganda. Sex Transm Dis 2011; 38:111–116.
15. Nugent RP, Krohn MA, Hillier SL. Reliability of diagnosing bacterial vaginosis is improved by a standardized method of gram stain interpretation. J Clin Microbiol 1991; 29:297–301.
16. Koumans EH, Markowitz LE, Hogan V. Indications for therapy and treatment recommendations for bacterial vaginosis in nonpregnant and pregnant women: A synthesis of data. Clin Infect Dis 2002; 35:S152–S172.
17. Schmidt WP, Genser B, Barreto ML, et al. Sampling strategies to measure the prevalence of common recurrent infections in longitudinal studies. Emerg Themes Epidemiol 2010; 7:5.
18. Baumgartner R, Somorjai R. Graphical display of fMRI data: Visualizing multidimensional space. Magn Reson Imaging 2001; 19:283–286.
19. R Development Core Team. R: A language and environment for statistical computing. 2008. Vienna, Austria, R Foundation for Statistical Computing.
20. Brotman RM, Erbelding EJ, Jamshidi RM, et al. Findings associated with recurrence of bacterial vaginosis among adolescents attending sexually transmitted disease clinics. J Pediatr Adolesc Gynecol 2007; 20:225–231.
21. Witkin SS, Linhares IM, Giraldo P. Bacterial flora of the female genital tract: Function and immune regulation. Best Pract Res Clin Obstet Gynaecol 2007; 21:347–354.
22. Price LB, Liu CM, Johnson KE, et al. The effects of circumcision on the penis microbiome. PLoS One 2010; 5:e8422.
23. Gray RH, Kigozi G, Serwadda D, et al. The effects of male circumcision on female partners' genital tract symptoms and vaginal infections in a randomized trial in Rakai, Uganda. Am J Obstet Gynecol 2009; 200:42.e1–42.e7.
24. Watts DH, Springer G, Minkoff H, et al. The occurrence of vaginal infections among HIV-infected and high-risk HIV-uninfected women: Longitudinal findings of the women's interagency HIV study. J Acquir Immune Defic Syndr 2006; 43:161–168.
25. Jamieson DJ, Duerr A, Klein RS, et al. Longitudinal analysis of bacterial vaginosis: Findings from the HIV epidemiology research study. Obstet Gynecol 2001; 98:656–663.
26. Ravel J, Gajer P, Abdo Z, et al. Vaginal microbiome of reproductive-age women. Proc Natl Acad Sci USA 2011; 108(suppl 1):4680–4687.
© Copyright 2011 American Sexually Transmitted Diseases Association
27. Koumans EH, Sternberg M, Bruce C, et al. The prevalence of bacterial vaginosis in the United States, 2001–2004: Associations with symptoms, sexual behaviors, and reproductive health. Sex Transm Dis 2007; 34:864–869.