|
Products Related to West Nile, Dengue, Malaria, T.B, Chikungunya, Sars, Zika |
|
Product# 15007: Recombinant Treponema pallidum 47kD Membrane Protein (E.coli) |
|
Product# 15009: Recombinant Treponema pallidum TmpA Membrane Protein (E.coli) |
|
Product# 15005: Recombinant Treponema pallidum 17kD Membrane Protein (E.coli) |
Abstract
We measured the microbial community structure of genital ulcers in women. Swabs from clinically detected ulcers were tested for HSV-2 and Treponema pallidum by polymerase chain reaction (PCR). HSV-2 and T. pallidum were detected by serum antibody testing. Microbial community structure was characterized by high-throughput 16 s rRNA gene amplicon sequencing. Multiple group testing and Elastic net and Lasso regressions identified taxa associated with differences in factors of interest. Among 49 ulcer specimens from 49 HSV-2 seropositive women, by PCR HSV-2 was recovered from 28 (57%) specimens and T. pallidum from none; one woman showed serologic evidence of syphilis. Overall, 63% of women were HIV-positive and 49% had an uncircumcised male sex partner. By both multiple group testing and regression, Porphyromonas (FDR p-value = 0.02), Prevotella (FDR p-value = 0.03), Anaerococcus (FDR p-value = 0.07), and Dialister (FDR p-value = 0.09) were detected at higher relative abundance in HSV-2 PCR-positive than negative ulcers. The presence of HSV-2 in a lesion was associated with presumed bacterial agents of Bacterial vaginosis. Differences in bacterial communities may contribute to HSV-2 ulcer pathogenesis, severity, or prolonged healing. If these results are confirmed, future studies may consider the influence of BV treatment on women’s GUD and HSV-2 incidence and recurrence.
Introduction
Medical male circumcision (MMC) reduces the risk of heterosexual HIV acquisition in African men by ~60%1–3. Among Ugandan MMC trial participants, penile anaerobic bacteria were reduced following circumcision4. The reduction in penile anaerobic bacteria appears to be transmissible, with a 40–60% reduced risk of bacterial vaginosis (BV) among female sex partners of circumcised men5. MMC reduced the risk of genital ulcer disease (GUD) in all three MMC trials by 37–52%2,3,6, and reduced GUD by 22% in female sex partners of circumcised men in the Ugandan trial5. In our MMC trial in Kenya, despite reduction in GUD, we found no reduction in HSV-2 incidence6,7. To identify the cause of GUD reduction, we conducted high-throughput amplicon sequencing of microbial 16 S rRNA genes on clinically detected genital ulcer specimens. Fusobacteria and Sneathia had nearly 6-fold greater odds of recovery from clinically detected non-herpetic and non-syphilitic GUD compared to STI-associated GUD8. These two BV-associated anaerobes are more common among uncircumcised men4,8. Bacteria from these genera have inflammatory and tissue destructive properties causing oral mucosal ulcers9–11. In light of this novel yet plausible mechanism by which anaerobic bacteria might contribute to genital ulcer pathogenesis, and the relationships between male circumcision status, penile microbiome, GUD, and HIV, we sought to describe the bacterial community structure of GUD in women and whether that differed by ulcer recovery of HSV-2, HIV status, male partner’s circumcision status, and clinical presentation.
Go to:
Methods
This study received ethical approval from the Institutional Review Boards of the University of Illinois at Chicago and University of Nairobi. All methods were performed in accordance with relevant guidelines and regulations. Women aged 18 years and older were recruited from sexually transmitted infection (STI) clinic visits where genital ulcers were clinically detected. Consenting women were compensated 200 Kenyan Shillings (~$2.00 USD) for their time. HIV testing was offered to all women and conducted according to Kenyan national guidelines.
Physical Examination Findings of Genital Ulcer Disease
All patients underwent standardized medical examination and history as part of routine clinical care. External genital and speculum examination by trained clinicians recorded the presence or absence of genital ulcers, the location and number of ulcers. Clinicians swabbed genital ulcers using a Dacron or polyester tipped swab. Women presenting with genital ulcer complaint were told about the study prior to history and examination and if interested, informed consent was obtained at that time. If genital ulcer was detected during the examination, only women who provided informed consent were invited to participate in the study after the completion of the history and examination.
Socio Demographic and Behavioral Data Collection
Sociodemographic and behavioral risk information obtained by personal interview included age, marital status, educational attainment, employment status, number of circumcised and uncircumcised male sex partners in the last 30 days and last 6 months, male sex partner’s circumcision status at the last sexual intercourse, condom use at last sex and frequency in the past 6 months, and other behaviors expected to be associated with genital microflora exposure.
Detection of HSV-2 and Treponema pallidum
Serum specimens were tested at the study clinic for HSV-2 antibody (Kalon HSV-2 IgG ELISA, Kalon Biological Limited, Aldershot, United Kingdom) according to the manufacturer’s instructions using recommended cutoffs. Syphilis infection was detected using the rapid plasma reagin test (RPR) (Macro-Vue™, Becton Dickinson, New Jersey, United States), confirmed by the T. pallidum hemagglutination (TPHA) assay (Randox Laboratories Ltd., Ardmore, United Kingdom).
Genital ulcer specimens were frozen at −20 °C, and shipped to Johns Hopkins University STD Research Laboratory for detection of HSV-1, HSV-2 and T. pallidum by multiplex polymerase chain reaction (PCR). Swabs were treated with 350 µl of 1X TRIS-EDTA (TE) buffer and vortexed for 15 seconds. The swab was removed from the TE buffer with a pair of tweezers, and the swab was expressed on the side of the tube ~3 times by rotating and pressing it into the tube. The swab was added back to the tube, and 200 µl of the 1X TE buffer was pipetted into a sample processing tray for extraction with the Roche MagNA Pure LC robot, leaving 150 µl for amplicon sequencing. All samples were extracted utilizing the Roche MagNA Pure LC robot with DNA Isolation Kit I, following the DNA I Blood/Cells High Performance extraction protocol and PCR was performed for HSV-1, HSV-2, and T. pallidum 12.
Analysis of microbial community structure
Genomic DNA (gDNA) from ulcers was PCR amplified and prepared for next-generation sequencing (NGS) using a modified two-step targeted amplicon sequencing approach, similar to that described previously13,14, with primers 515 F and 806 R15, targeting the V4 variable region of Bacterial and Archaeal small subunit (SSU) ribosomal RNA (rRNA) gene. The primers contained 5′ common sequence tags (common sequence 1 and 2 [CS1 and CS2]) as described previously16. The forward primer, CS1–515F (ACACTGACGACATGGTTCTACAGTGCCAGCMGCCGCGGTAA) and reverse primer, CS2-806R (TACGGTAGCAGAGACTTGGTCTGGACTACHVGGGTWTCTAAT) were synthesized by Integrated DNA Technologies (IDT; Coralville, Iowa) as standard oligonucleotides. The PCR amplifications were performed in 10 microliter reactions in 96-well plates. A mastermix for the entire plate was made using the 2X AccuPrime SuperMix II (Life Technologies, Gaithersburg, MD). Final concentration of CS1_515 F and CS2_806 R primers was 500 nM. Cycling conditions were as follows: 95 °C for 5 minutes, followed by 28 cycles of 95 °C for 30”, 55 °C for 45” and 68 °C for 30”. A final, 7 minute elongation step was performed at 68 °C. A single microliter of the first amplification reaction was transferred to a second amplification reaction with the same mastermix conditions, but with different primers. No purification was performed. Each well received a unique primer pair obtained from the Access Array Barcode Library for Illumina sequencers (Fluidigm, South San Francisco, CA; Item# 100–4876), which contain Illumina sequencing adapters, a sample-specific barcode (reverse primer), and CS1 or CS2 sequences. Cycling conditions were as follows: 95 °C for 5 minutes, followed by eight cycles of 95 °C for 30”, 60 °C for 30” and 68 °C for 30”. A final, 7 minute elongation step was performed at 68 °C. Samples were pooled in equal volume and purified using solid phase reversible immobilization (SPRI) cleanup, implemented with AMPure XP beads (Beckman Coulter, Brea, CA, USA) at a ratio of 0.6X (v:v) SPRI solution to sample volume.
The library pool was spiked with 15% non-indexed PhiX control library provided by Illumina and loaded onto a MiSeq v2 flow cell at a concentration of 8 pM for cluster formation and sequencing. Sequencing was performed at the W.M. Keck Center for Comparative and Functional Genomics at the University of Illinois at Urbana-Champaign (UIUC), and data were analyzed using the Casava1.8 pipeline.
Raw sequence data were imported into the software package CLC genomics workbench (v7.0; CLC Bio, Qiagen, Boston, MA). Sequences were merged, quality trimmed (Q20), and filtered to remove residual phiX contamination, and exported as FASTA files. Subsequently, sequence data were processed using the software package QIIME (v1.8.0;17). Briefly, sequences were screened for chimeras using the usearch61 algorithm18 using de novo and reference-based detection methods, and putative chimeric sequences were removed from the dataset. All the reads were then clustered into operational taxonomic units (OTU), using a threshold similarity of 97%. Representative sequences from each OTU were extracted, and these sequences were classified using the “assign_taxonomy” algorithm implementing the RDP classifier, with the Greengenes reference OTU build (v13_8;19). A taxon-by-sample abundance matrix (biological observation matrix, BIOM;20) was generated at taxonomic levels from phylum to genus using the “make_OTU_table” algorithm. The BIOM was rarefied to 1,250 sequences per sample within QIIME to avoid analytical issues associated with variable sequence number between samples21.
Statistical Analysis
We used three analytic approaches to compare the bacterial community of genital ulcers by detection of HSV-2, male sex partner circumcision status, and HIV status: (1) test for global differences in bacterial communities by factors of interest; (2) compare diversity of bacterial communities by factors of interest; (3) identify specific bacterial taxa differing by factors of interest. Analyses were conducted at the taxonomic level of genus. Additionally, we compared clinical characteristics by recovery of HSV-2 from ulcers. Differences between categorical explanatory variables and factors were assessed by the chi-square test, or Fisher’s exact test when cell size was < 5. Inferential analyses were conducted using STATA/SE 13.0 for Windows (Stata Corp., College Station, TX).
First, we tested global bacterial communities by factors using analysis of similarity (ANOSIM). Rarefied sequence data were transformed [log(X + 1)] prior to generation of the Bray-Curtis resemblance matrix. Second, diversity at the genus level was measured by the Shannon diversity index (log base e) and was compared by factors with Wilcoxon rank sum test. ANOSIM and Shannon diversity index analyses were conducted in the software package Primer6 (Primer-E, version 6.1.13, United Kingdom). Third, to identify specific bacterial taxa that differed by HSV-2 PCR result, male partner circumcision status, and HIV status, we compared the number of sequence reads of each genus using the Kruskal-Wallis test. We report test statistic p-values and false discovery rate (FDR)-corrected p-values. Group significance testing was performed within QIIME, employing the Kruskal-Wallis non-parametric test on untransformed number of sequence reads, with calculation of false discovery rate (FDR) corrected p-values17. To verify the robustness of these results and considering the dependence between bacterial taxa, we also applied regression to identify specific bacterial bacteria that differed by factors. Elastic net regression was chosen as our primary regression approach as its method of variable selection is not bounded by the number of samples. It uses an additional ridge regression penalty to nullify this constraint22. For the validation process, k-fold cross validation was used as the sample size was not sufficient to be divided into sizable training and validation sets. We applied numerous values of k in order to find out the change in the parameters of our equation. We tested numerous values for the ridge regression parameter to identify the optimum model for appropriate variable selection. To assess the robustness of the results of Elastic net, LASSO regression was implemented23 and we identified 7 bacteria that were selected by both Elastic net and LASSO regression as having an association with any of the outcomes (HSV-2 status, HIV status, male sex partner circumcision status) at the p < 0.10 level. We applied linear regression to these 7 bacteria for each of the three outcomes, and report coefficients, p-values and adjusted R-squared values, which was multivariable adjusted when a bacterium was found to be significantly associated (p < 0.10) with more than one dependent variable. Before the analysis by linear regression, the data was standardized with µ = 0 and σ = 1. Regression analyses were executed in SAS software, version 9.4 (SAS Institute Inc., Cary, NC, USA).
Principal coordinates analysis was used to visualize bacterial community differences by the three dependent variables (HSV-2 PCR status, HIV status, and male sex partner circumcision status). We calculated both within and between group weighted Unifrac distances to find out the phylogenetic distances between the different groups of samples according to the dependent variables. The Weighted Unifrac distances and principal coordinate analysis plots were operationalized using the phyloseq, ggplot2, plyr, ape and wesanderson packages in the R 3.3.2 environment (R Development Core Team).
Data Access
The 16 S rRNA gene amplicon sequence data from this study have been submitted to the NCBI Sequence ReadArchive (http://www.ncbi.nlm.nih.gov/Traces/sra/sra.cgi) under accession number PRJNA281808.
Go to:
Results
For this pilot study, the target sample size was 60 subjects, and 58 subjects were enrolled in the allotted timeline. Data were rarified to 1,250 sequences, which resulted in maintaining 56 subjects. Overall, 49 (88%) women were HSV-2 seropositive and 7 (12%) were HSV-2 seronegative. Among the 49 HSV-2 seropositive women, by PCR HSV-2 was detected in 28 (57%) ulcers and T. pallidum in none. Serologically, syphilis was detected in 1 (2.0%) woman (titer 1:4; TPHA-positive); there were no sequence reads from the genus Treponema detected in this specimen (0/34,176 sequences before rarifaction). One serologically negative subject had an OTU annotated as Treponema (30/47,102 sequences before rarefaction) with HSV-2 recovered from the ulcer. Among the 7 HSV-2 seronegative women, there was no serologic evidence of syphilis. HSV-2 was PCR-detected in three (38%) of seven ulcers from HSV-2 seronegative women while T. pallidum was not detected; amplicon sequencing did not identify any Treponema spp. sequence reads in HSV-2 seronegative women. All GUD specimens were negative for HSV-1. Therefore, no STI etiology was identified in four (7%) of 56 women with clinically detected genital ulcers. Given the small number of HSV-2 seronegative subjects, inferential analyses were restricted to the 49 HSV-2 seropositive women. From these 49 specimens we obtained 1,664,532 sequences with an average read-length of 250 base pairs. An average of 33,970 sequences were obtained from each sample (min = 1,261, max = 64,939).
Socio-Demographic and Behavioral Characteristics
HSV-2 seropositive women were median age 26 years (range 18–59) and 76% married or cohabiting (Table 1). Overall, 60% of women reported having one sex partner in the past six months, with 17% reporting two or more sex partners, and 49% reported having any uncircumcised sex partner in the past six months. Male sex partner’s circumcision status at the last sex, past 30 days, and past six months were highly correlated (Pearson correlation coefficient = 0.68–0.79); subsequent analyses proceeded with circumcision status in the past six months as most women reported one sex partner in the past six months. Having an uncircumcised male sex partner did not differ between women with HSV-2 PCR positive (50%) or negative (47%) ulcers (p = 0.86). Overall, 63% of women were HIV positive, and this did not differ significantly between women with HSV-2 PCR positive or negative ulcers (71% vs. 50%, p = 0.13). There were no statistically significant differences in sociodemographic, behavioral characteristics, or the joint distribution of HIV status and male sex partner circumcision status between women with HSV-2 PCR positive or negative ulcers (Table 1).
Table 1
Subject Characteristics by HSV-2 mPCR Results.
|
Characteristic* |
HSV-2 mPCR Results |
P-Value |
|
|
Positive, N = 28 n (%) |
Negative, N = 21 n (%) |
||
|
Median age in years |
26.5 |
26 |
0.87 |
|
HIV status |
|||
|
Negative |
8 (29) |
10 (50) |
0.13 |
|
Positive |
20 (71) |
10 (50) |
|
|
Antiretroviral use reported |
8 (40) |
5 (50) |
|
|
Married or cohabiting |
|||
|
No |
8 (29) |
4 (19) |
0.52 |
|
Yes |
20 (71) |
17 (81) |
|
|
Educational attainment |
|||
|
Primary or less |
17 (61) |
16 (76) |
0.25 |
|
Some secondary or higher |
11 (39) |
5 (24) |
|
|
Used antibiotics, past 1 month |
|||
|
No |
12 (43) |
8 (38) |
0.74 |
|
Yes |
16 (57) |
13 (62) |
|
|
^Uncircumcised male sex partner last 6 months |
|||
|
No |
14 (50) |
10 (53) |
0.86 |
|
Yes |
14 (50) |
9 (47) |
|
|
HIV status and male partner circumcision status |
|||
|
HIV negative, circumcised partner |
6 (33) |
3 (11) |
0.25 |
|
HIV positive, circumcised partner |
3 (17) |
10 (37) |
|
|
HIV negative, uncircumcised partner |
3 (17) |
4 (15) |
|
|
HIV positive, uncircumcised partner |
6 (33) |
10 (37) |
|
|
Number of sex partners |
|||
|
None |
8 (29) |
3 (15) |
0.55 |
|
One |
15 (53) |
14 (70) |
|
|
Two or more |
5 (18) |
3 (15) |
|
|
Condom used at last sex |
|||
|
No |
12 (43) |
14 (67) |
0.1 |
|
Yes |
16 (57) |
7 (33) |
|
|
Median Shannon diversity index [95% CI] |
2.1 [1.8–2.3] |
1.9 [1.3–2.1] |
0.05 |
|
Genus-level, log base e |
|||
|
Median number of genus-level bacterial taxa [95% CI] |
26.5 [22–35] |
17 [13–23] |
0.02 |
Open in a separate window
^Uncircumcised represents having any uncircumcised sex partner in the past 6 months; Circumcised represents all sex partners in the past 6 months were circumcised.
*Not all cells sum to N due to missing responses. +Recall period for sexual behaviors is past 6 months, unless otherwise specified.
CI = Confidence Interval.
Clinical Characteristics by HSV-2 PCR Results
Compared to HSV-2 positive ulcers, the treating clinician was more likely to describe HSV-2 PCR negative ulcers as papular (32% vs. 4%, p = 0.02), and less likely to describe HSV-2 PCR negative ulcers as “patch” (0% vs. 21%, p = 0.03), and hyperpigmented (18% vs. 47%, p = 0.08; Table 2). The clinical impression was associated with PCR detection of HSV-2: 54% of ulcers in which HSV-2 was recovered were diagnosed as “GUD of herpetic appearance” compared to 24% of those that were HSV-2 negative (p = 0.05).
Table 2
Clinical Characteristics by HSV-2 mPCR Results.
|
Characteristic* |
HSV-2 mPCR Results |
P-value |
|
|
Positive, N = 28 |
Negative, N = 21 |
||
|
n (%) |
n (%) |
||
|
Symptoms^ |
|||
|
Vaginal discharge (n = 48) |
10 (37) |
10 (48) |
0.46 |
|
Genital itching |
13 (46) |
12 (57) |
0.46 |
|
Painful urination |
17 (61) |
10 (48) |
0.36 |
|
Lower abdominal pain (n = 48) |
4 (15) |
5 (24) |
0.48 |
|
Number of ulcers |
|||
|
1 |
14 (50) |
9 (43) |
|
|
2–3 |
11 (39) |
5 (24) |
|
|
4–10 |
3 (11) |
7 (33) |
0.14 |
|
Median duration of ulcer (IQR) (n = 43) |
7 (4–25.5) |
6.5 (4–10.5) |
0.24 |
|
Ulcer location^ |
|||
|
Labia majora |
23 (82) |
15 (71) |
0.37 |
|
Labia minora |
5 (18) |
4 (19) |
1 |
|
Introitus |
3 (11) |
2 (10) |
1 |
|
Vault |
2 (7) |
2 (10) |
1 |
|
Ulcer description^ |
|||
|
Multiple (vs. Singular) |
15 (54) |
12 (57) |
0.8 |
|
Macular (n = 47) |
22 (79) |
11 (58) |
0.13 |
|
Papular (n = 46) |
1 (4) |
6 (32) |
0.02 |
|
Pustular (n = 46) |
0 (0) |
2 (11) |
0.17 |
|
Patch (n = 48) |
6 (21) |
0 (0) |
0.03 |
|
Hyperpigmentation (n = 36) |
9 (47) |
3 (18) |
0.08 |
|
Hypopigmentation (n = 36) |
5 (26) |
8 (47) |
0.2 |
|
Other Examination findings^ |
|||
|
Discharge from os (n = 45) |
10 (40) |
7 (35) |
0.73 |
|
Cervical friability/bleeding (n = 45) |
2 (8) |
1 (5) |
1 |
|
Clinical impression^ |
|||
|
GUD, herpetic |
15 (54) |
5 (25) |
0.05 |
|
GUD, non-herpetic |
15 (54) |
15 (75) |
0.13 |
|
Vaginitis |
3 (11) |
6 (32) |
0.13 |
|
Medications dispensed+ |
|||
|
Acyclovir only |
5 (18) |
5 (24) |
|
|
Acyclovir, penicillin, erythromycin |
10 (36) |
3 (14) |
|
|
Penicillin, erythromycin |
7 (25) |
12 (57) |
|
|
Erythromycin only |
3 (11) |
1 (5) |
|
|
Penicillin only |
2 (7) |
0 (0) |
|
|
Acyclovir, erythromycin |
1 (4) |
0 (0) |
0.13 |
Open in a separate window
*Not all cells sum to N due to missing data; the sample size for variables with missing responses are indicated.
^Not mutually exclusive.
+Acyclovir dispensed as cream or tabs; Penicillin dispensed as Penicillin G 2.4 million units.
IQR = Interquartile range.
Global Comparison of Microbial Community Structure by mPCR Results
Analysis of similarity (Table 3) showed that the difference in bacterial community composition between HSV-2 PCR-positive and PCR-negative women was statistically significant (p = 0.002). There were no statistically significant differences in community composition by HIV status or male sex partner circumcision status. The principal components plots are complementary visualization to these analytic results: the community centroid for samples from women with HSV-2 PCR positive ulcers are distinct from the centroid of samples from women whose ulcers were HSV-2 PCR negative (Fig. 1), while centroids representing bacterial communities by HIV status (Fig. 2) and male sex partner circumcision status (Fig. 3) show much greater overlap and shorter distance between centroid.
Table 3
Results of Analysis of Similarity Measures: Comparison of Microbial Community Structure by Patient Factors.
|
Variable |
Global R statistic |
P-value |
|
All samples (N = 49) |
||
|
HSV-2 PCR results |
0.148 |
0.002 |
|
Circumcision status of male sex partner |
0.010 |
0.277 |
|
HIV status |
0.011 |
0.315 |
|
HSV-2 PCR Positive (N = 28) |
||
|
Circumcision status of male sex partner |
−0.028 |
0.766 |
|
HIV status |
−0.050 |
0.693 |
|
HSV-2 PCR Negative (N = 21) |
||
|
Circumcision status of male sex partner |
0.165 |
0.053 |
|
HIV status |
0.273 |
0.011 |
Figure 1
Principal coordinates analysis of Weighted UniFrac values between HSV-2 PCR positive (blue) and HSV-2 PCR negative (red) samples with the axes scaled by the percentage of the variance that they contain to summarize the microbial community compositional differences between samples. Each point represents a single sample, and the distance between points represents how compositionally different the samples are from one another.
Figure 2
Principal coordinates analysis of Weighted UniFrac values between HIV positive positive (blue) and HIV negative (red) samples.
Figure 3
Principal coordinates analysis of Weighted UniFrac values between samples from women with uncircumcised male sex partners (blue) and circumcised male sex partners (red).
As reflected by the median Shannon index (Table 1), there was greater diversity in microbial communities from ulcers in which HSV-2 was detected compared to ulcers in which HSV-2 was not detected (1.9 vs. 2.1, p = 0.05). The median number of taxa recovered was also higher for ulcers in which HSV-2 was detected (26.5 vs. 17, p = 0.02).
Specific Bacterial Taxa Differing by HSV-2 Status, HIV Status, and Male Sex Partner Circumcision Status
In Kruskal-Wallis analysis, compared to HSV-2 PCR negative ulcers (Table 4), HSV-2 PCR positive ulcers had higher mean sequence count for bacteria from the genera Porphyromonas (FDR p = 0.02) and Prevotella (FDR p = 0.03), with strong trends for Anaerococcus (FDR p = 0.07) and Dialister (FDR p = 0.09) (Table 3). While Finegoldia, Gardnerella, Corynebacterium, and Anaerococcus demonstrated higher mean sequence count for ulcers from HIV-positive women compared to ulcers from HIV-negative women (p < 0.05 each), only the difference for Finegoldia trended toward FDR significance (FDR p = 0.09). This was echoed in regression analyses, with a statistically significant 0.81 standard deviation reduction in Finegoldia for HIV-negative women (Table 5). Only the mean sequence of Gardnerella differed by male sex partner circumcision status (p = 0.01), and the FDR p-value was non-significant (p = 0.20). Results of regression (Table 5) also identified Porphyromonas as a significant predictor of HSV-2 positivity and Finegoldia with HIV positivity.
Table 4
Mean sequence reads of bacterial taxa by HSV-2 PCR status, HIV status, and male sex partner circumcision status with Kruskall-Wallis group tests of significance
|
Bacteria |
HSV-2 PCR status |
P-value |
FDR P-value |
HIV Status |
P-value |
FDR P-value |
Male partner circumcision status |
P-value |
FDR P-value |
|||
|
Positive |
Negative |
Positive |
Negative |
Uncircum |
Circum |
|||||||
|
Porphyromonas |
22 |
1 |
<0.01 |
0.02 |
13 |
13 |
0.79 |
0.79 |
10 |
17 |
0.72 |
0.95 |
|
Prevotella |
156 |
80 |
<0.01 |
0.03 |
145 |
95 |
0.85 |
0.90 |
141 |
113 |
0.09 |
0.41 |
|
Anaerococcus |
111 |
52 |
0.01 |
0.07 |
81 |
99 |
0.04 |
0.15 |
90 |
88 |
0.87 |
0.94 |
|
Dialister |
18 |
18 |
0.02 |
0.09 |
15 |
24 |
0.79 |
0.90 |
13 |
24 |
0.37 |
0.72 |
|
Sneathia |
100 |
22 |
0.03 |
0.11 |
93 |
27 |
0.10 |
0.29 |
112 |
29 |
0.12 |
0.41 |
|
Finegoldia |
62 |
46 |
0.17 |
0.44 |
33 |
95 |
0.01 |
0.09 |
51 |
63 |
0.83 |
0.94 |
|
Peptoniphilus |
93 |
93 |
0.25 |
0.57 |
84 |
113 |
0.48 |
0.59 |
123 |
73 |
0.45 |
0.72 |
|
Acinetobacter |
27 |
87 |
0.41 |
0.81 |
46 |
66 |
0.11 |
0.29 |
27 |
65 |
0.13 |
0.41 |
|
Corynebacterium |
41 |
44 |
0.53 |
0.88 |
24 |
75 |
0.03 |
0.15 |
45 |
43 |
0.97 |
0.97 |
|
Staphylococcus |
19 |
101 |
0.56 |
0.88 |
21 |
111 |
0.15 |
0.29 |
13 |
88 |
0.11 |
0.41 |
|
Gardnerella |
88 |
95 |
0.70 |
0.88 |
125 |
41 |
0.01 |
0.11 |
125 |
50 |
0.01 |
0.20 |
|
Lactobacillus |
165 |
254 |
0.75 |
0.88 |
176 |
190 |
0.32 |
0.51 |
164 |
226 |
0.50 |
0.72 |
|
Streptococcus |
21 |
21 |
0.77 |
0.88 |
13 |
58 |
0.14 |
0.29 |
35 |
26 |
0.47 |
0.72 |
|
Aerococcus |
2 |
13 |
0.96 |
0.96 |
11 |
0 |
0.15 |
0.40 |
14 |
1 |
0.29 |
0.77 |
Open in a separate window
FDR P-value = False discovery rate corrected P-value
“Uncircum” = Uncircumcised; “Circum” = Circumcised.
We present results from the 14 bacteria representing the top 10 bacteria of significance (ranked by P-value) for each factor. Because some bacteria were in the top 10 for some outcomes but not others (e.g., Gardnerella for HIV status and male sex partner circumcision status, but not for HSV-2 status), this leads to presentation of more than 10 bacteria total.
Table 5
Results of linear regression: Bacteria associated with HSV-2, HIV, and circumcision status.
|
Genus |
HSV-2 positivity |
HIV positivity |
Uncircumcised male partner |
P-value |
R-squared Adjusted |
|
Porphyromonas |
0.91 (0.26) |
— |
— |
0.0010 |
0.1901 |
|
Peptostreptococcus |
0.83 (0.27) |
— |
— |
0.0029 |
0.1556 |
|
Adlercreutzia |
0.80 (0.27) |
— |
— |
0.0043 |
0.1430 |
|
Finegoldia |
— |
−0.81 (0.28) |
— |
0.0049 |
0.1411 |
|
WAL_1855D |
0.70 (0.27) |
— |
— |
0.0135 |
0.1044 |
|
Sneathia |
0.66 (0.28) |
— |
0.44 (0.28) |
0.0237 |
0.1180 |
|
Gardnerella |
— |
— |
0.59 (0.28) |
0.0391 |
0.0710 |
|
Phascolarctobacterium |
— |
— |
0.53 (0.29) |
0.0699 |
0.0505 |
Interpretation: These results represent 7 different models for the 7 bacteria identified in Elasticnet and Lasso regression as associated with HSV-2 ulcer status, HIV status, and/or male sex partner’s circumcision status at the p < 0.10 level. As Sneathia was associated with both HSV-2 status and male sex partner circumcision status, the model controlled for both factors. Data are standardized (µ = 0 and σ = 1) and we report the coefficient (standard error), model p-value and adjusted R-squared.
Go to:
Discussion
Among HSV-2 seropositive women, HSV-2 was detected by PCR in 57% of clinically identified ulcers, and the bacterial community of these ulcers differed by whether HSV-2 was detected with PCR. All women were HSV-2 seropositive; the extent to which it matters that GUD of putatively similar etiology has variation in bacterial community composition relates to the potential association between these different bacterial community structures with co-infections or impact on severity of disease or treatment.
This study was conducted as response to our previous finding that 55% of men with clinically detected genital ulcers were serologically negative for HSV-2 and syphilis, and 39% also had no HSV-2 or T. pallidum recovered by PCR8. This high proportion of non-etiologically defined GUD is in keeping with studies throughout Sub-Saharan Africa24–29 and India30,31. Using amplicon sequencing in our study of men, we found that anaerobic bacteria associated with BV were more commonly recovered from genital ulcers in these men without HSV-2 or T. pallidum. Our current findings present two striking differences to these previous results. First, a minority of women (9%) with GUD were without STI-associated etiology. Subjects in the current study self-presented with symptoms, whereas in the previous study, samples were collected from men taking part in a randomized clinical trial of MMC to reduce HIV incidence3. Men in the trial underwent physical examination every 6 months, and were encouraged to come to the clinic for symptoms at interim visits. HIV was rare (4%) in our prior analysis as men had to be HIV negative for trial entry, compared to 63% HIV-positivity among women in the current study. HIV can alter the natural history of genital herpes, producing more severe and frequent outbreaks and delayed healing32. Therefore, in the trial we increased our chance of detecting non-STI related ulcers and ulcers that may not have been severe enough to lead to health care seeking otherwise. Secondly, our present study found that ulcers in which HSV-2 was detected by PCR had more BV-associated bacteria – specifically, Porphyromonas, Prevotella, Peptostreptococcus, Sneathia, and possibly Anaerococcus and Dialister – than ulcers from which HSV-2 was not detected. Ulcers in which HSV-2 was not detected by PCR may have been at a different point in healing, and therefore HSV-2 was no longer detectable. However, the median duration of ulcer symptoms did not differ by whether or not HSV-2 was recovered from the ulcer. In a cohort of Kenyan women, Kaul et al. found that baseline HSV-2 infection was associated with increased incidence of BV33. At the same time, prospective epidemiologic and clinical investigations support the increased risk of HSV-2 acquisition among women with BV34. Whether HSV-2 increases risk of BV and/or BV increases risk of HSV-2, such association would increase the likelihood of co-detection of HSV-2 and BV-associated bacteria in cross-sectional study design. Though published data regarding the effect of BV on the natural history of HSV-2 are limited, Ursell et al. demonstrate that BV-associated bacteria (G. vaginalis, A. vaginae) have been associated with increased duration and number of HSV-2 related outbreaks35.
A primary public health concern is that among HSV-2 seropositive individuals there is increased risk of HIV acquisition36 and transmission37, highlighted by the 63% HIV prevalence among our sample of HSV-2 infected women. We detected significant global bacterial community dissimilarity between HIV-positive and HIV-negative women only among ulcers that were HSV-2 PCR negative. We observed that Finegoldia was significantly reduced among women with HIV. Prospective studies of the effect of HIV on the vaginal microbiome are limited, but our recent analysis of a cohort of U.S. women showed that HIV-status did not alter the change in the vaginal microbiome over an 8- to 10-year period38.
The microbial community of ulcers from women with an uncircumcised partner reflected more BV-associated bacteria. This is consistent with what we expected based on our previous findings and those of others showing that female partners of uncircumcised men are more likely to have BV, and that this stems from the greater abundance of BV-associated anaerobic bacteria among uncircumcised men compared to circumcised men39. However, in our study this difference was only observed in ulcers from which HSV-2 was not recovered. We are unable to determine why the differences by HIV status and male partner circumcision status were detected only among HSV-2 PCR negative ulcers. Although these women were all HSV-2 seropositive, the GUD detected could have had different or multiple etiologic factors. This may be supported by findings of the differences in ulcer characteristics between HSV-2 PCR positive and negative ulcers, and overall, the clinical impression was less often “GUD, herpetic” for PCR negative (25%) than positive ulcers (54%).
Limitations
We recruited women seeking care at an STI clinic, which likely biased our sample towards those with more severe and HSV-2 related GUD. Data on male partner circumcision status were self-reported by women, and there may have been misclassification. In a study of women’s knowledge of partners’ male circumcision status in Uganda, 1.2% of women misreported the status of confirmed circumcised men and 8.2% of women misreported the status of uncircumcised men40, though study of this in Zambia and Swaziland find misclassification rates 11–15% for either status41. We did not conduct species-level analysis, which will be necessary in subsequent studies with more specific hypotheses. This pilot study had low power to detect modest associations for individual bacteria, especially in the context of multiple testing correction, which was necessary for hypothesis exploration. However, we used multiple analytic approaches to increase confidence in the robustness of results: detection of differences is supported by results of ANOSIM, Bray-Curtis dissimilarity measure, Kruskal-Wallis group testing, and regression analyses, with convergence on specific bacteria. We did not collect data on frequency of previous GUD, occurrence of prodrome, or clinician characterized stage of healing. We did not diagnose BV by Amsel’s criteria or Nugent score, as the clinic at which the study took place provides only syndromic diagnosis. This information would have aided interpretation of the observed differences in bacterial community structure. These data emphasize the need for longitudinal study to understand the relationships between the vaginal microbiome and women’s risk of HSV-2, HIV, BV, and the role of male sex partner’s circumcision status in these outcomes.
Go to:
Conclusions
Among HSV-2 seropositive women, the bacterial community of clinically detected GUD was more diverse and had higher abundance of BV-associated bacteria when HSV-2 was detected in the ulcer. This is of public health and clinical importance as differences in bacterial communities may contribute to HSV-2 related lesion induction, ulcer pathogenesis, severity, or prolonged healing. If these results are confirmed, future studies may consider the influence of BV treatment on women’s GUD and HSV-2 incidence and recurrence.
Go to:
Acknowledgements
This research was supported by grant number P01AI082971 (PI Landay) from the Division of AIDS, National Institute of Allergies and Infectious Disease of the United States National Institutes of Health.
Go to:
Author Contributions
S.D.M. Conceived study hypotheses and design; analyzed data; wrote first draft of manuscript A.K.P. Selected analysis approach, conducted regression analyses, generated plots; wrote analytic methods section of manuscript S.J.G. Oversaw amplicon sequencing procedures and selection and application of methods; interpreted analysis results; wrote laboratory methods section of manuscript; contributed to presentation of results; critical review and revision of manuscript A.N. Executed bioinformatics pipeline; conducted analyses in QIIME, contributed to presentation of results. E.O.J. Assisted with design of data collection methods; developed study protocols and ensured integrity of implementation; regulatory maintenance and compliance; critical review and revision of manuscript C.A.G. Oversaw PCR procedures and selection and application of methods; contributed to interpretation of results; critical review and revision of manuscript S.B. Contributed to study protocol develop and implementation; trained study staff; contributed to interpretation of results A.L. Contributed to interpretation of results; critical review and revision of manuscript R.C.B. Assisted with development of study design and protocol; contributed to interpretation of results; critical review and revision of manuscript.
Go to:
Notes
Go to:
Competing Interests
The authors declare that they have no competing interests.
Go to:
Footnotes
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Go to:
References
- Auvert B, et al. Randomized, controlled intervention trial of male circumcision for reduction of HIV infection risk: the ANRS 1265 Trial. PLoS Med. 2005;2:e298. doi: 10.1371/journal.pmed.0020298. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Gray RH, et al. Male circumcision for HIV prevention in men in Rakai, Uganda: a randomised trial. Lancet. 2007;369:657–66. doi: 10.1016/S0140-6736(07)60313-4. [PubMed] [CrossRef] [Google Scholar]
- Bailey RC, et al. Male circumcision for HIV prevention in young men in Kisumu, Kenya: a randomised controlled trial. Lancet. 2007;369:643–56. doi: 10.1016/S0140-6736(07)60312-2. [PubMed] [CrossRef] [Google Scholar]
- Price LB, et al. The effects of circumcision on the penis microbiome. PLoS One. 2010;5:e8422. doi: 10.1371/journal.pone.0008422. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Gray RH, 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–7. [PMC free article] [PubMed] [Google Scholar]
- Mehta SD, et al. Circumcision status and incident HSV-2 infection, genital ulcer disease, and HIV infection. AIDS. 2012;26:1141–49. doi: 10.1097/QAD.0b013e328352d116. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Mehta SD, et al. Medical male circumcision and herpes simplex virus 2 acquisition: posttrial surveillance in Kisumu, Kenya. J. Infect. Dis. 2013;208:1869–76. doi: 10.1093/infdis/jit371. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Mehta SD, et al. Microbial diversity of genital ulcer disease in men enrolled in a randomized trial of male circumcision in Kisumu, Kenya. PLoS One. 2012;7:e38991. doi: 10.1371/journal.pone.0038991. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Brook I. Microbiology and management of peritonsillar, retropharyngeal, and parapharyngeal abscesses. Oral Maxillofac. Surg. 2004;62:1545–50. doi: 10.1016/j.joms.2003.12.043. [PubMed] [CrossRef] [Google Scholar]
- Enwonwu CO, et al. Pathogenesis of cancrum oris (noma): confounding interactions of malnutrition with infection. Am. J. Trop. Med. Hyg. 1999;60:223–32. doi: 10.4269/ajtmh.1999.60.223. [PubMed] [CrossRef] [Google Scholar]
- Wilson M. Biological activities of lipopolysaccharides from oral bacteria and their relevance to the pathogenesis of chronic periodontitis. Sci. Prog. 1995;78:19–34. [PubMed] [Google Scholar]
- Suntoke TR, et al. Evaluation of multiplex real-time PCR for detection of Haemophilus ducreyi, Treponema pallidum, herpes simplex virus type 1 and 2 in the diagnosis of genital ulcer disease in the Rakai District, Uganda. Sex. Transm. Infect. 2009;8:97–101. [PMC free article] [PubMed] [Google Scholar]
- de Carcer DA, Denman SE, McSweeney C, Morrison M. Strategy for modular tagged high-throughput amplicon sequencing. Appl. Environ. Microbiol. 2011;77:6310–12. doi: 10.1128/AEM.05146-11. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Bybee SM, et al. Targeted amplicon sequencing (TAS): A scalable next-gen approach to multilocus, multitaxa phylogenetics. Genome Biol. Evol. 2011;3:1312–23. doi: 10.1093/gbe/evr106. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Caporaso JG, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc. Natl. Acad. Sci. USA. 2011;108(Suppl 1):4516–22. doi: 10.1073/pnas.1000080107. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Moonsamy PV, et al. High throughput HLA genotyping using 454 sequencing and the Fluidigm Access Array™ System for simplified amplicon library preparation. Tissue Antigens. 2013;81:141–49. doi: 10.1111/tan.12071. [PubMed] [CrossRef] [Google Scholar]
- Caporaso JG, et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods. 2010;7:335–36. doi: 10.1038/nmeth.f.303. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Haas BJ, et al. Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res. 2011;21:494–504. doi: 10.1101/gr.112730.110. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- McDonald D, et al. An improved greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2012;6:610–18. doi: 10.1038/ismej.2011.139. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- McDonald D., et al. The biological observation matrix (BIOM) format or: How I learned to stop worrying and love the ome-ome. Gigascience. 1, 7-217X-1-7 (2012). [PMC free article] [PubMed]
- Gihring TM, Green SJ, Schadt CW. Massively parallel rRNA gene sequencing exacerbates the potential for biased community diversity comparisons due to variable library sizes. Environ. Microbiol. 2012;14:285–90. doi: 10.1111/j.1462-2920.2011.02550.x. [PubMed] [CrossRef] [Google Scholar]
- Zou H, Hastie T. Regularization and variable selection via the elastic net. J. R. Statist. Soc. B. 2005;67:301–20. doi: 10.1111/j.1467-9868.2005.00503.x. [CrossRef] [Google Scholar]
- Tibshirani R. Regression shrinkage and selection via the lasso. J. R. Statist. Soc. B. 1996;58:267–88. [Google Scholar]
- Plummer FA, et al. Clinical and microbiologic studies of genital ulcers in Kenyan women. Sex. Transm. Dis. 1985;12:193–97. doi: 10.1097/00007435-198510000-00005. [PubMed] [CrossRef] [Google Scholar]
- Moodley P, et al. Association between HIV-1 infection, the etiology of genital ulcer disease, and response to syndromic management. Sex. Transm. Dis. 2003;30:241–45. doi: 10.1097/00007435-200303000-00013. [PubMed] [CrossRef] [Google Scholar]
- O’Farrell N, Hoosen AA, Coetzee KD, van den Ende J. Sexual behaviour in Zulu men and women with genital ulcer disease. Genitourin. Med. 1992;68:245–48. [PMC free article] [PubMed] [Google Scholar]
- Paz-Bailey G, et al. Changes in the etiology of sexually transmitted diseases in Botswana between 1993 and 2002: implications for the clinical management of genital ulcer disease. Clin. Infect. Dis. 2005;41:1304–12. doi: 10.1086/496979. [PubMed] [CrossRef] [Google Scholar]
- Behets FM, et al. Chancroid, primary syphilis, genital herpes, and lymphogranuloma venereum in Antananarivo, Madagascar. J. Infect. Dis. 1999;180:1382–85. doi: 10.1086/315005. [PubMed] [CrossRef] [Google Scholar]
- Ahmed HJ, et al. Etiology of genital ulcer disease and association with human immunodeficiency virus infection in two Tanzanian cities. Sex. Transm. Dis. 2003;30:114–19. doi: 10.1097/00007435-200302000-00004. [PubMed] [CrossRef] [Google Scholar]
- Prabhakar P, et al. Genital ulcer disease in India: etiologies and performance of current syndrome guidelines. Sex. Transm. Dis. 2012;39:906–10. doi: 10.1097/OLQ.0b013e3182663e22. [PubMed] [CrossRef] [Google Scholar]
- Risbud A, et al. The etiology of genital ulcer disease by multiplex polymerase chain reaction and relationship to HIV infection among patients attending sexually transmitted disease clinics in Pune, India. Sex. Transm. Dis. 1999;26:55–62. doi: 10.1097/00007435-199901000-00009. [PubMed] [CrossRef] [Google Scholar]
- Gupta R, Warren T, Wald A. Genital herpes. Lancet. 2007;370:2127–37. doi: 10.1016/S0140-6736(07)61908-4. [PubMed] [CrossRef] [Google Scholar]
- Kaul R, et al. Prevalent herpes simplex virus type 2 infection is associated with altered vaginal flora and an increased susceptibility to multiple sexually transmitted infections. J. Infect. Dis. 2007;196:1692–97. doi: 10.1086/522006. [PubMed] [CrossRef] [Google Scholar]
- Cherpes TL, et al. Genital tract shedding of herpes simplex virus type 2 in women: effects of hormonal contraception, bacterial vaginosis, and vaginal group B Streptococcus colonization. Clin. Infect. Dis. 2005;40:1422–28. doi: 10.1086/429622. [PubMed] [CrossRef] [Google Scholar]
- Ursell LK, et al. Comparison of the vaginal microbial communities in women with recurrent genital HSV receiving acyclovir intravaginal rings. Antiviral. Res. 2014;102:87–94. doi: 10.1016/j.antiviral.2013.12.004. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Freeman EE, et al. Herpes simplex virus 2 infection increases HIV acquisition in men and women: systematic review and meta-analysis of longitudinal studies. AIDS. 2006;20:73–83. doi: 10.1097/01.aids.0000198081.09337.a7. [PubMed] [CrossRef] [Google Scholar]
- Gray RH, et al. Probability of HIV-1 transmission per coital act in monogamous, heterosexual, HIV-1–discordant couples in Rakai, Uganda. Lancet. 2001;357:1149–53. doi: 10.1016/S0140-6736(00)04331-2. [PubMed] [CrossRef] [Google Scholar]
- Mehta SD, et al. The vaginal microbiota over an 8- to 10-year period in a cohort of HIV-infected and HIV-uninfected women. PLoS One. 2015;10:e0116894. doi: 10.1371/journal.pone.0116894. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Liu CM, et al. Penile microbiota and female partner bacterial vaginosis in Rakai, Uganda. MBio. 2015;6:e00589. [PMC free article] [PubMed] [Google Scholar]
- Kong X, et al. The accuracy of women’s reports of their partner’s male circumcision status in Rakai, Uganda. AIDS. 2013;27:662–64. doi: 10.1097/QAD.0b013e32835c557c. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
- Hewett PC, Haberland N, Apicella L, Mensch BS. The (mis)reporting of male circumcision status among men and women in Zambia and Swaziland: A randomized evaluation of interview methods. PLoS One. 2012;7:e36251. doi: 10.1371/journal.pone.0036251. [PMC free article] [PubMed] [CrossRef] [Google Scholar]
Articles from Scientific Reports are provided here courtesy of Nature Publishing Group
