Immune milieu and microbiome of the distal urethra in Ugandan men: impact of penile circumcision and implications for HIV susceptibility
Microbiome volume 10, Article number: 7 (2022)
Coronal sulcus (CS) anaerobe abundance and IL-8 levels are linked to HIV acquisition, and are dramatically reduced after penile circumcision (PC). The distal urethra may be the site of some HIV acquisition before PC, and presumably most acquisition post PC. We describe the immune milieu and microbiome of the distal urethra in uncircumcised Ugandan men, and define the impact of PC. Participants consisted of HIV-negative, genital symptom-free adult Ugandan men undergoing PC (n = 51). Urethral and coronal sulcus swabs were collected at baseline and at 6- and 12-months post-PC. Soluble immune factors were quantified by multiplex ELISA, and bacterial abundance assessed by 16S rRNA qPCR and sequencing.
At baseline, the urethra was enriched compared to the CS for most cytokines (including IL-8 and MIP-1β) and soluble E-cadherin (sE-cadherin, an epithelial disruption marker), although CS levels of IL-1α and IL-1β were higher. Baseline total bacterial abundance was ≥ 20-fold higher in the CS than the urethra (median 27,100 vs. 1200 gene copies/swab, p = 0.001), and anaerobes comprised 58% of CS bacteria vs. 42% of urethral bacteria. PC did not alter urethral IL-8 (median 806 at baseline vs. 1130 pg/ml at 12 months; p = 0.062) and urethral sE-cadherin increased (113,223 vs. 158,385 pg/ml, p = 0.009), despite five- and sevenfold drops in total bacterial and anaerobe abundance after PC, respectively. However, PC dramatically reduced CS levels of sE-cadherin (15,843 vs. 837 pg/ml, p < 0.001) and most cytokines (IL-8; 34 vs. 3 pg/ml, p < 0.001), while reducing total bacterial and anaerobe abundance by 13-fold and 60-fold, respectively (both P ≤ 0.004).
The urethra is immunologically rich with characteristics of an HIV-susceptible tissue site. However, PC had no impact on urethral immunology and may have reduced epithelial integrity, despite modest reductions in total bacteria and anaerobes, suggesting that HIV protection from PC is not mediated via immune or microbiome alterations in the urethra.
Penile circumcision (PC) reduces the risk of human immunodeficiency virus type 1 (HIV) acquisition by almost two-thirds, with protection mediated through several proposed physical and biological mechanisms [1,2,3]. This one-time, cost-effective HIV prevention tool has been rolled out on a large scale in resource-limited regions with a high HIV burden . However, in many regions of sub-Saharan Africa, up to 50% of eligible men decline the procedure [5,6,7]. Therefore, a better understanding of the biological mechanism(s) for PC-mediated protection may inform the design of alternative prevention tools focused on uncircumcised men who decline PC.
At least part of the mechanism by which PC protects against HIV infection is stochastic, mediated through the reduced surface area of penile tissues exposed to HIV-containing genital fluids during sex and the direct removal of tissue-associated HIV-susceptible target cells that include activated CD4+ T lymphocyte subsets, Langerhans cells, macrophages, and dendritic cells [8, 9]. In addition, foreskin removal obliterates the subpreputial space, where HIV-containing genital fluids may be retained between the glans penis and inner foreskin after sex, thereby prolonging virus-tissue contact. In addition to these physical alterations, surgical obliteration of the subpreputial space via PC exposes the coronal sulcus (CS) to air, profoundly altering both the penile microbiome composition and immune milieu. The total bacterial load and the proportion of anaerobes is dramatically reduced in the CS after PC, with a shift from anaerobes typically associated with bacterial vaginosis (e.g., Gardnerella, Prevotella, and Peptostreptococcus) toward facultative anaerobic bacteria generally considered to be normal skin flora (e.g., Staphylococcus and Corynebacteria). There are also substantial reductions in the levels of the chemoattractant cytokine interleukin (IL)-8 [9,10,11], a chemokine which has been directly linked to the increased density of anaerobes in the prepuce, to an increased density of foreskin HIV target cells, and to an increased risk of HIV acquisition risk . This suggests that the microenvironment of the subpreputial space in an uncircumcised man sustains an anaerobic microbiome that induces host tissue inflammation and thereby increases HIV susceptibility .
There is still residual HIV acquisition after PC, and the urethra is thought to be the major site for virus acquisition in circumcised men . Indeed, it has been hypothesized that the urethra may be the tissue site for more virus acquisition in uncircumcised men than has been appreciated, and that a major mechanism underpinning PC-mediated protection could be alterations in the microbiome and immunology of the penile urethra, instead of the CS . In keeping with this, the penile urethra contains several potential HIV target cell subsets, including CD4+ CCR5+ macrophages, dendritic cells, and activated CD4+ T cells [14, 15] while urethral macrophages may act as an HIV reservoir after HIV acquisition and treatment . Urethral innate and immune defenses include physical flushing during urination, mucus secretion to trap pathogens, antimicrobial factors that include defensins and lysozyme, and immunoglobulin production (especially secretory immunoglobulin A, IgA) . The anatomy of the urethra naturally limits aeration, and it is possible that the overhanging foreskin of an uncircumcised penis might further enhance the anaerobic environment and/or further prolong urethral contact with HIV-containing coital secretions. However, PC does not protect against other urethrally acquired genital infections such as gonorrhea and chlamydia , despite reducing the risk of “skin acquired” infections such as human papillomavirus (HPV), herpes simplex virus type 2 (HSV-2), and genital ulcer disease [19, 20]. In order to assess the impact of PC on the urethral environment and microbiome, we prospectively enrolled Ugandan men presenting for elective PC into a longitudinal clinical protocol, with repeat sampling of the distal urethra and CS prior to and after surgical PC.
A total of 51 participants were enrolled prior to elective PC, with 46 (90%) and 35 (69%) of participants attending scheduled 6- and 12-month post-PC follow-up visits. The mean age of participants was 22 years (Table 1); 45 men (88%) reported having ever had vaginal intercourse, with 87% of these reporting multiple (> 1) lifetime sexual partners, and just over half reporting multiple female partners within the previous 6 months. A minority of men (4%) reported antibiotic exposure during the preceding 3 months, and almost all (94%) reported regular retraction of the foreskin during washing.
Immunology and bacteriology of the urethra and coronal sulcus in uncircumcised men
At baseline (pre-PC), the immune milieu of the urethra and coronal sulcus was distinct (Table 2). Specifically, the distal urethra was significantly enriched for multiple cytokines and other immune factors compared to the CS, including IL-8, matrix metallopeptidase 9 (MMP-9), macrophage inflammatory protein (MIP)-1β, resistin and tissue inhibitor of metalloproteinases (TIMP)-1 (all P < 0.001). Vascular endothelial growth factor (VEGF) levels did not differ between tissue sites (P > 0.10), and IL-1α and IL-1β levels were relatively enriched in the CS compared to the distal urethra (P < 0.001) (Table 2). Given the role of epithelial barrier integrity in impairing HIV entry and the potential epithelial disruption induced by inflammatory cytokines , we also assayed levels of soluble E-cadherin (sE-cadherin), a tight junction protein, and found that these were significantly elevated in the distal urethra compared to CS (P < 0.001) (Table 2). In the CS, baseline levels of sE-cadherin correlated strongly with several proinflammatory biomarkers, including increased IL-8 levels (r = 0.624; P < 0.0001), while no correlation was observed between levels of sE-cadherin and inflammatory biomarkers in the distal urethra.
In contrast to most immune factors, baseline (pre-PC) total bacterial abundance, defined as the 16S ribosomal RNA copy number/swab, was 17-fold higher in the CS than the distal urethra (median 25,980,000 vs. 1,468,000 gene copies/swab, p = 0.001). In addition, anaerobes constituted a higher proportion of the CS microbiome, where they comprised 58% of the microbiome, compared to 42% in the distal urethra (Fig. 1A, Table 3). The composition of the urethral microbiome was also quite distinct from the CS microbiome, with Streptococcus being the most abundant and comprising an average of 16% of urethral bacteria, followed by Veillonella, Prevotella, unclassified Clostridiales, Corynebacterium, and Anaerococcus at 10%, 10%, 8%, 6%, and 5% respectively. By contrast, Prevotella was the most abundant CS bacteria (mean 15%), followed by Corynebacterium, Finegoldia, unclassified Clostridiales, and Peptoniphilus at 11%, 11%, 9%, and 8% respectively (Fig. 1A, Table 3). Notably, key female genital tract bacteria, particularly Lactobacillus and Gardnerella, were rare at both penile sites.
Differential impact of penile circumcision on the coronal sulcus versus distal urethra immune milieu
Levels of the proinflammatory chemoattractant IL-8 at the CS fell significantly 6 months after PC (median 5.39 vs. 33.87 pg/ml, p < 0.001), with further reductions at 12 months (3.41 vs. 33.87 pg/ml, p < 0.001; Fig. 2A). Similar falls after PC were observed for IL-1β, resistin, sE-cadherin, and VEGF (Fig. 2B-C, E-F), although there was no impact of PC on CS levels of MIP-1β or MMP-9 (data not shown). Unexpectedly, levels of the proinflammatory cytokine IL-1α in the coronal sulcus actually increased significantly at both 6 months (788.48 vs. 358.77 pg/ml; p < 0.001) and 12 months (538.02 vs. 358.77 pg/ml, p = 0.018; Fig. 2D) after PC.
In contrast to the profound immune alterations seen at the coronal sulcus, there was unexpectedly little impact of PC on the urethral immune milieu. Urethral IL-8 levels were unaltered at 6 months post PC (754.96 vs. 805.59 pg/ml; P = 0.906) and actually tended to increase at 12 months (1130.21 vs. 805.59 pg/ml; P = 0.062; Fig. 3A). Other immune biomarkers remained unchanged, with the exception of a transient increase in urethral IL-1α levels at 6 months (67.81 vs. 45.90 pg/ml; P = 0.032; Fig. 3D) and a sustained increase in levels of sE-cadherin at both 6 months (179,590.44 vs. 113,223.06 pg/ml; P = 0.032) and 12 months (158,384.77 vs. 113,223.06 pg/ml; P = 0.009; Fig. 3E).
Penile circumcision and the coronal sulcus and urethral microbiome
Total bacterial load remained higher at the CS than the urethra at follow-up visits (Fig. 1B), although there were significant reductions at both sites post-PC, with a disproportionate reduction in anaerobic taxa (Fig. 1A). At the CS, there was a >13-fold reduction in total bacterial load at both 6 months and 12 months, while reductions in the urethra were less marked (sixfold) and only significant at 6 months (Fig. 1B).
Broadly, post-PC, there was a shift toward “skin-associated” gram-positive aerobes at the CS (from 58% anaerobe dominance at baseline to 90% aerobe dominance at 12 months), largely composed of Corynebacterium (50%) and Staphylococcus (20%). At the distal urethra, the aerobe Streptococcus dominated the microbiome both pre- and post-PC (16% and 26% respectively). PC increased the proportion of Corynebacterium both at the CS and in the urethra (Fig. 1A).
Circumcision substantially reduced the absolute abundance of anaerobic taxa at both penile tissue sites. This was most profound in the CS, where there was a 60-fold drop (from 12,733,000 to 204,098 gene copies per swab; P < 0.001) but was also significant in the urethra (14-fold drop). Specifically, at the CS, PC significantly reduced the absolute abundances of several specific bacterial genera linked to BV and HIV risk in both men and women. PC profoundly reduced the abundance of most BV-associated bacterial genera in the CS at 12 months, particularly Prevotella (median 0 vs. 2,159,328 gene copies per swab, P < 0.0001; Fig. 4A). While there had been a marked proportionate increase in gram positive aerobes after PC, their absolute abundance did not significantly change (e.g., Corynebacterium, 894,487 vs. 272,834, P = 0.512; and Staphylococcus, 370,878 vs. 21,426, P = 0.120; Fig. 4E, G). However, there were similar sustained drops in the absolute abundance of several BV-associated taxa, including Prevotella and Peptoniphilus (0 vs. 40,604 gene copies/swab, P = 0.02; and 0 vs. 34,433 gene copies/swab, P < 0.0001 respectively; Fig. 4B, D). Similar to the CS, despite the proportionate change, there was no increase in the absolute abundance of skin commensals Staphylococcus and Corynebacterium at the urethra (P > 0.5; Fig. 4F, H).
Immune associations of specific penile bacterial genera
The bacterial vaginosis (BV)-associated taxa Prevotella and Dialister in the CS have been most strongly associated with HIV seroconversion in uncircumcised men, while a CS microbiome enriched for skin-associated taxa such as Corynebacterium and Staphylococcus is associated with HIV protection after PC . Therefore, we next assessed how the absolute abundances of these taxa correlated with immune parameters at both the CS and urethra prior to PC, specifically with IL-8 (a chemoattractant cytokine linked to penile HIV acquisition) and E-cadherin (a tight junction protein linked to epithelial disruption).
At the CS, Dialister abundance correlated positively with levels of IL-8 and sE-cadherin (Spearman coefficients 0.400, p = 0.006 and 0.659, p < 0.0001; Fig. 5), while Corynebacterium abundance was negatively correlated with both parameters (Spearman coefficient −0.326, p = 0.029 and −0.307, p = 0.037; Fig. 5). Likewise, CS Prevotella abundance tended to be positively associated with these immune parameters (Spearman coefficients 0.183, p = 0.229 and 0.600, p < 0.0001) and Staphylococcus abundance negatively (Spearman coefficients −0.214, p = 0.158 and −0.420, p = 0.004). Bacterial-immune associations were weaker in the distal urethra, although Prevotella abundance was again positively associated with IL-8 levels (Spearman coefficients 0.350, p = 0.018), and both Corynebacterium and Staphylococcus abundance with reduced sE-cadherin (Spearman coefficient −0.304, p = 0.043, Fig. 5 and −0.306, p = 0.041 respectively).
Overall, at the CS of uncircumcised men, skin-associated bacterial taxa were associated with less inflammation and enhanced epithelial integrity, while seroconversion-associated taxa were associated with more inflammation and reduced integrity; similar but weaker associations were seen in the distal urethra.
The HIV protection that is afforded by PC [1,2,3] may be mediated via several mechanisms, including direct removal of HIV-susceptible foreskin tissue [22, 23] and reductions in inflammatory cytokines and pro-inflammatory anaerobes in the coronal sulcus . The urethra is assumed to be the site of residual penile HIV acquisition after PC, but it has been hypothesized that this tissue site may play a greater role than appreciated prior to PC, due to urethral immune and microbiome parameters induced by an overhanging foreskin, so that HIV protection post-PC might relate to an altered urethral immune or microbiological milieu . We addressed this research question in a cohort of adult Ugandan men presenting for elective PC. As expected from prior studies [9, 10, 24], we demonstrated that PC dramatically reduced most pro-inflammatory anaerobes and inflammatory cytokines/chemokines in the coronal sulcus, with simultaneous reductions in soluble E-cadherin suggesting enhanced epithelial integrity. Prior to PC, the microbiome and immune milieu of the distal urethra were quite distinct to the coronal sulcus, and the impact of PC also differed at this tissue site. The urethra of the uncircumcised penis had a much lower bacterial load than the CS, with a smaller proportion of anaerobes, and there were higher levels of IL-8, sE-cadherin, and most other soluble immune parameters (with the exception of IL1α and IL1β). However, while PC substantially reduced the bacterial load and proportionately enriched the microbiome of both tissue sites for skin-associated gram-positive aerobes, this was much less marked in the urethra than the coronal sulcus, and PC had a minimal impact on urethral cytokines and chemokines while increasing urethral sE-cadherin levels. Overall, these results suggest that HIV protection following PC is not likely to be mediated through alterations in the urethral microbiome or urethral immune milieu.
Differences at the urethra and coronal sulcus in levels of proinflammatory/chemoattractant biomarkers such as IL-8, and the contrasting effects of PC on the local immune milieu, may have important implications for HIV susceptibility. IL-8 is produced by epithelial and antigen presenting cells, and serves as both a chemoattractant and proinflammatory cytokine [25, 26]. Elevated genital IL-8 levels are associated with enhanced HIV acquisition in both women [27, 28] and uncircumcised men , where they correlate with an increased density of CD4+ T cell targets in the endocervix and foreskin, respectively [9, 29]. Although baseline urethral IL-8 levels were over 20-fold higher in the urethra than the coronal sulcus, urethral levels of IL-8 and other cytokines/chemokines were unaltered by PC, despite modest reductions in the urethral bacterial load and a microbiome shift toward skin-associated gram-positive aerobes. In addition, no association was seen between urethral anaerobes and local inflammation, even in those men with the highest anaerobic bacterial load at baseline. It remains unclear whether urethral anaerobes simply do not induce inflammation at this site, whether additional factors such as increased local epithelial disruption (evidenced by high E-cadherin levels) have a greater contribution to inflammation, or whether cytokines induced by inflammation are not able to accumulate due to intermittent “flushing” during urination.
Soluble E-cadherin levels appear to be a biomarker of epithelial disruption at other mucosal sites [30,31,32]. This transmembrane glycoprotein connects epithelial cells at adherens junctions and is integral in mediating cell adhesion and contact inhibition of proliferation ; therefore, elevated levels indicate disrupted cell-cell junctions and a lack of epithelial integrity. Prior to PC, urethral sE-cadherin levels were almost tenfold higher than in the coronal sulcus, despite a lower total bacterial load and reduced anaerobic burden. While G. vaginalis induced higher soluble sE-cadherin levels and the expression of IL-8 and IL-6 in both the murine genital tract  and a human ex vivo cervicovaginal model , we only observed a correlation between inflammatory anaerobes, IL-8 levels and sE-cadherin in the coronal sulcus, but not the urethra. In addition, there was an unexpected increase in urethral sE-cadherin levels post-PC, which we hypothesize, may relate to minor post-surgical irritation of the now “unshielded” distal urethra (e.g., through rubbing on clothing).
Finally, we observed an unexpected and substantial increase IL-1α levels after PC at the CS, despite profound reductions in pro-inflammatory anaerobes and a reduction in other cytokines and chemokines (including IL-8). Although perhaps not intuitive, this may make biological sense: while IL-1α is predominantly thought of as a proinflammatory cytokine, it also mediates several epidermal barrier functions that include stimulation of keratinocyte proliferation and differentiation, upregulation of genes associated with cell adhesion and the synthesis of lipids needed for the formation and maturation of the stratum corneum [34,35,36]. As a result, all epithelial cells (especially epidermal cells) constitutively express relatively high levels of IL-1α, and it is possible that elevated levels may relate to healing of the epidermal barrier [37, 38].
There are some limitations to our study that merit discussion. We were not able to directly link immune/microbiome changes seen at either penile tissue site with actual HIV acquisition, since this would have required a much larger sample size. However, profound immune and bacterial changes were seen at both sites, and the lack of alteration in urethral immunology post-PC, coupled with elevated sE-cadherin, strongly suggests that the protection mediated by PC is unlikely to be mediated urethrally. Nonetheless, additional studies will be needed to confirm and expand our observations, potentially with the examination of additional immune markers and an assessment of PC impact on other microbes such as fungi and viruses. Furthermore, sE-cadherin is an indirect marker of epithelial disruption, and ideally future studies would correlate this with direct microscopic evidence of epithelial damage. We were unable to assess immune cell populations in the urethra, such as T cells, dendritic cell subsets, neutrophils, and other innate immune cells, due to practical difficulties in tissue sampling at this site. However, previous studies using cadaveric or surgical specimens have demonstrated that many of these cell subsets are abundant in the urethra [15, 39], can serve as a reservoir for HIV in infected men taking ART , and would be expected to be profoundly altered by the local microbiome and immune milieu.
In summary, PC does not alter urethral immunology and may reduce urethral epithelial integrity, despite reductions in urethral bacterial load and a modest reduction in the relative proportion of anaerobes. In contrast, PC dramatically reduced the CS bacterial load and anaerobe proportion, enhanced epithelial integrity in the CS, and reduced most CS inflammatory chemokines/cytokines. This suggests that HIV protection post-PC is mediated by the removal of inflamed, HIV-susceptible foreskin tissues rather than by immunologic or microbiome alterations in the urethra.
Study enrollments follow-up visits and sample processing
Study participants consisted of uncircumcised HIV-negative Ugandan men who did not have any genital STI symptoms, aged at least 18 years electively presenting at the Rakai Health Sciences Program for voluntary male medical PC to reduce their HIV risk. Penile samples were collected at baseline (pre-PC) and at 6- and 12-months following PC. A social-behavioral questionnaire was also administered at each study visit. The study clinician used sterile polyester tips (Puritan Medical Products, ME, USA) premoistened in phosphate-buffered saline (PBS) to swab the inner foreskin (pre-PC) or coronal sulcus (post-PC), and nylon flocked urethral swabs (Hardy Diagnostics, CA, USA) were used to swab the distal urethra. Swabs were immediately placed into 500 μL of PBS and transferred to the laboratory on ice. In the laboratory, swabs were vigorously vortexed for 60 s and then the swab head was inverted prior to a quick spin to dry out the swab, which was then discarded. Each sample tube was pulse-vortexed and two aliquots each containing 250 μL of supernatant were logged and frozen at −80 °C.
Multiplex chemiluminescent ELISA
Levels of 9 soluble immune biomarkers were assayed using a multiplex electro-chemiluminescence ELISA platform (Meso Scale Discovery, Rockville, MD) by research personnel fully blinded to PC status. The biomarker panel included the prototypic proinflammatory cytokines interleukin 1 alpha (IL-1α), and interleukin 1 beta (IL-1β); the chemoattractant chemokines interleukin 8 (IL-8) and macrophage inflammatory protein 1 beta (MIP-1β); and a biomarker of epithelial integrity/breakdown (E-cadherin). Other exploratory analytes on the panel included resistin, an atypical proinflammatory biomarker, and the novel biomarkers tissue inhibitor of metalloproteases 1 (TIMP-1), vascular endothelial growth factor (VEGF), and matrix metalloproteinase 9 (MMP-9). Analyte concentrations within each sample were calculated from a standard curve generated using serial dilutions of stock analyte from the manufacturer and this was prepared for each plate run. Lower limits of detection derived from each plate run were applied to samples flagged as below detection for both duplicates irrespective of the coefficient of variation. All samples belonging to a given individual were run on the same plate to limit any potential impact of plate-plate variation. On each plate, a frozen control media aliquot was plated to monitor inter plate/run variability. Our primary immune endpoint was IL-8, since coronal sulcus levels of IL-8 were the strongest immune predictor of HIV acquisition in uncircumcised Ugandan men.
Microbiome laboratory analysis and controls used
DNA was extracted from 80 μL of diluted swab eluent using a combination of enzymatic and chemical lysis. Briefly, each sample was treated with an enzymatic cocktail containing 122 μL Tris-EDTA, 50 μL 10 mg/mL lysozyme (L6876-1G, Sigma-Aldrich), 4 μL 25 KU/mL mutanolysin (M4782-5KU, Sigma-Aldrich), and 3 μL 4 U/μL lysostaphin (SAE0091-2MG, Sigma-Aldrich) at 37 °C for 1 h, followed by extraction using MagMax DNA Multi-Sample Ultra 2.0 Kit (including Proteinase K treatment) with 80 μL final elution volume. Penile microbiome analysis was characterized by 16S rRNA gene-based broad-range real-time PCR  and sequencing. The sequencing analysis was performed using a modified protocol from Fadrosh et al.  with forward (341F) and reverse (786R) primers from Liu et al. . Sequencing was performed on MiSeq platform using MiSeq Reagent Kit v3 (600 cycle).
During processing, primer sequences were removed using cutadapt v2.4  and the resultant sequences were quality trimmed using Trimmomatic v0.39 . DADA2 v1.10  modules were used for reads-filtering, chimera check, and inferred error models to identify amplicon sequence variants (ASVs). The ASVs were classified at each taxonomic level at 80% bootstrap confidence level using the Naïve Bayesian Classifier (v.2.12) . Classification results for each sample were enumerated to generate an abundance matrix for analysis. Additional details can be found at https://github.com/araclab/mb_analysis.
Using the resultant qPCR and sequencing outputs, absolute abundance of each penile bacterial genus and species was calculated as: Absolute abundance of a taxon per swab = Total bacterial load per swab (measured by qPCR as total copies of 16S rRNA gene per swab) × proportional abundance of the given taxon (measured by sequencing as the number of sequences assigned to a taxon in a given sample, divided by the total number of sequences obtained for the sample). Sequence data for this study can be accessed at SRA project number PRJNA738496.
Negative extraction control (NEC) is included with each batch of extraction and analyzed by qPCR, where Cp > 33-35 are considered acceptable. NECs are further included in sequencing analysis. No template control (NTC) and positive template controls (PTC) are included for each amplicon PCR plate and analyzed to assess cross-contamination during PCR and verify PCR performance.
Microbiome statistical analysis
Intra-individual paired analysis was performed using the Wilcoxon-matched pairs signed-rank test in SPSS version 24 (Armonk, New York, USA). Intra-individual baseline cytokine comparisons were performed between the two sites (coronal sulcus vs. urethra) as well as cytokine levels at baseline compared to months 6 and 12. Non-parametric Spearman rank-order correlation was run between specific cytokines/parameters with probable biological interactions. Analyses were visualized using GraphPad Prism version 6 (La Jolla, CA, USA). Missing data (e.g., arising from a missed visit) was excluded from all paired intra-individual analysis.
Availability of data and materials
Sequence data for this study can be accessed at SRA project number PRJNA738496. Details for the bioinformatics analyses can be found at https://github.com/araclab/mb_analysis. Please contact author for additional requests.
Auvert B, Taljaard D, Lagarde E, Sobngwi-Tambekou J, Sitta R, Puren A. Randomized, controlled intervention trial of male circumcision for reduction of HIV infection risk: the ANRS 1265 Trial. PLoS Med. 2005;2(11):e298.
Bailey RC, Moses S, Parker CB, Agot K, Maclean I, Krieger JN, et al. Male circumcision for HIV prevention in young men in Kisumu, Kenya: a randomised controlled trial. Lancet. 2007;369(9562):643–56.
Gray RH, Kigozi G, Serwadda D, Makumbi F, Watya S, Nalugoda F, et al. Male circumcision for HIV prevention in men in Rakai, Uganda: a randomised trial. Lancet. 2007;369(9562):657–66.
Njeuhmeli E, Forsythe S, Reed J, Opuni M, Bollinger L, Heard N, et al. Voluntary medical male circumcision: modeling the impact and cost of expanding male circumcision for HIV prevention in eastern and southern Africa. PLoS Medicine. 2011;8(11):e1001132.
Ledikwe JH, Nyanga RO, Hagon J, Grignon JS, Mpofu M, Semo BW. Scaling-up voluntary medical male circumcision - what have we learned? HIV AIDS (Auckl). 2014;6:139–46.
Moyo S, Mhloyi M, Chevo T, Rusinga O: Men’s attitudes: a hindrance to the demand for voluntary medical male circumcision--a qualitative study in rural Mhondoro-Ngezi, Zimbabwe. Glob Public Health 2015, 10(5-6):708-720.
Sgaier SK, Reed JB, Thomas A, Njeuhmeli E. Achieving the HIV prevention impact of voluntary medical male circumcision: lessons and challenges for managing programs. PLoS Med. 2014;11(5):e1001641.
Prodger JL, Gray R, Kigozi G, Nalugoda F, Galiwango R, Hirbod T, et al. Foreskin T-cell subsets differ substantially from blood with respect to HIV co-receptor expression, inflammatory profile, and memory status. Mucosal Immunol. 2012;5(2):121–8.
Prodger JL, Gray RH, Shannon B, Shahabi K, Kong X, Grabowski K, et al. Chemokine levels in the penile coronal sulcus correlate with HIV-1 acquisition and are reduced by male circumcision in Rakai, Uganda. PLoS Pathog. 2016;12(11):e1006025.
Liu CM, Hungate BA, Tobian AA, Serwadda D, Ravel J, Lester R, Kigozi G, Aziz M, Galiwango RM, Nalugoda F et al: Male circumcision significantly reduces prevalence and load of genital anaerobic bacteria. mBio 2013, 4(2):e00076.
Liu CM, Prodger JL, Tobian AAR, Abraham AG, Kigozi G, Hungate BA, Aziz M, Nalugoda F, Sariya S, Serwadda D et al: Penile anaerobic dysbiosis as a risk factor for HIV infection. MBio 2017, 8(4).
Prodger JL, Kaul R. The biology of how circumcision reduces HIV susceptibility: broader implications for the prevention field. AIDS Res Ther. 2017;14(1):49.
Anderson D, Politch JA, Pudney J. HIV infection and immune defense of the penis. Am J Reprod Immunol. 2011;65(3):220–9.
Pudney J, Anderson DJ. Immunobiology of the human penile urethra. Am J Pathol. 1995;147(1):155–65.
Ganor Y, Zhou Z, Bodo J, Tudor D, Leibowitch J, Mathez D, et al. The adult penile urethra is a novel entry site for HIV-1 that preferentially targets resident urethral macrophages. Mucosal Immunol. 2013;6(4):776–86.
Ganor Y, Real F, Sennepin A, Dutertre CA, Prevedel L, Xu L, et al. HIV-1 reservoirs in urethral macrophages of patients under suppressive antiretroviral therapy. Nature Microbiol. 2019;4(4):633–44.
Pudney J, Anderson D. Innate and acquired immunity in the human penile urethra. J Reprod Immunol. 2011;88(2):219–27.
Mehta SD, Moses S, Agot K, Parker C, Ndinya-Achola JO, Maclean I, et al. Adult male circumcision does not reduce the risk of incident Neisseria gonorrhoeae, Chlamydia trachomatis, or Trichomonas vaginalis infection: results from a randomized, controlled trial in Kenya. J Infect Dis. 2009;200(3):370–8.
Tobian AA, Serwadda D, Quinn TC, Kigozi G, Gravitt PE, Laeyendecker O, et al. Male circumcision for the prevention of HSV-2 and HPV infections and syphilis. N Engl J Med. 2009;360(13):1298–309.
Weiss HA, Thomas SL, Munabi SK, Hayes RJ. Male circumcision and risk of syphilis, chancroid, and genital herpes: a systematic review and meta-analysis. Sexually Transmitted Infections. 2006;82(2):101–9 discussion 110.
Arnold KB, Burgener A, Birse K, Romas L, Dunphy LJ, Shahabi K, et al. Increased levels of inflammatory cytokines in the female reproductive tract are associated with altered expression of proteases, mucosal barrier proteins, and an influx of HIV-susceptible target cells. Mucosal Immunol. 2016;9(1):194–205.
Kigozi G, Wawer M, Ssettuba A, Kagaayi J, Nalugoda F, Watya S, et al. Foreskin surface area and HIV acquisition in Rakai, Uganda (size matters). AIDS. 2009;23(16):2209–13.
Kigozi G, Liu CM, Park D, Packman ZR, Gray RH, Kaul R, et al. Foreskin surface area is not associated with sub-preputial microbiome composition or penile cytokines. PLoS One. 2020;15(6):e0234256.
Price LB, Liu CM, Johnson KE, Aziz M, Lau MK, Bowers J, et al. The effects of circumcision on the penis microbiome. PLoS One. 2010;5(1):e8422.
Baggiolini M, Clark-Lewis I. Interleukin-8, a chemotactic and inflammatory cytokine. FEBS Letters. 1992;307(1):97–101.
Harada A, Sekido N, Akahoshi T, Wada T, Mukaida N, Matsushima K. Essential involvement of interleukin-8 (IL-8) in acute inflammation. J Leukocyte Biol. 1994;56(5):559–64.
Liebenberg LJ, Masson L, Arnold KB, McKinnon LR, Werner L, Proctor E, et al. Genital-systemic chemokine gradients and the risk of HIV acquisition in women. J Acquired Immune Deficiency Syndromes. 2017;74(3):318–25.
Masson L, Passmore JA, Liebenberg LJ, Werner L, Baxter C, Arnold KB, et al. Genital inflammation and the risk of HIV acquisition in women. Clin Infect Dis. 2015;61(2):260–9.
Shannon B, Gajer P, Yi TJ, Ma B, Humphrys MS, Thomas-Pavanel J, et al. Distinct effects of the cervicovaginal microbiota and herpes simplex type 2 infection on female genital tract immunology. J Infect Dis. 2017;215(9):1366–75.
Anton L, Sierra LJ, DeVine A, Barila G, Heiser L, Brown AG, et al. Common cervicovaginal microbial supernatants alter cervical epithelial function: mechanisms by which Lactobacillus crispatus contributes to cervical health. Front Microbiol. 2018;9:2181.
Furukawa F, Fujii K, Horiguchi Y, Matsuyoshi N, Fujita M, Toda K, et al. Roles of E- and P-cadherin in the human skin. Microscopy Research Technique. 1997;38(4):343–52.
Sierra LJ, Brown AG, Barila GO, Anton L, Barnum CE, Shetye SS, et al. Colonization of the cervicovaginal space with Gardnerella vaginalis leads to local inflammation and cervical remodeling in pregnant mice. PloS one. 2018;13(1):e0191524.
Garcia MA, Nelson WJ, Chavez N: Cell-cell junctions organize structural and signaling networks. Cold Spring Harbor Perspectives Biol 2018, 10(4).
Jiang YJ, Lu B, Crumrine D, Man MQ, Elias PM, Feingold KR. IL-1alpha accelerates stratum corneum formation and improves permeability barrier homeostasis during murine fetal development. J Dermatol Sci. 2009;54(2):88–98.
Jung YJ, Jung M, Kim M, Hong SP, Choi EH. IL-1alpha stimulation restores epidermal permeability and antimicrobial barriers compromised by topical tacrolimus. J Investigative Dermatol. 2011;131(3):698–705.
Barland CO, Zettersten E, Brown BS, Ye J, Elias PM, Ghadially R. Imiquimod-induced interleukin-1 alpha stimulation improves barrier homeostasis in aged murine epidermis. J Investigative Dermatol. 2004;122(2):330–6.
Hauser C, Saurat JH, Schmitt A, Jaunin F, Dayer JM. Interleukin 1 is present in normal human epidermis. J Immunol. 1986;136(9):3317–23.
Gahring LC, Buckley A, Daynes RA. Presence of epidermal-derived thymocyte activating factor/interleukin 1 in normal human stratum corneum. J Clin Investigation. 1985;76(4):1585–91.
Fischetti L, Barry SM, Hope TJ, Shattock RJ. HIV-1 infection of human penile explant tissue and protection by candidate microbicides. AIDS. 2009;23(3):319–28.
Liu CM, Aziz M, Kachur S, Hsueh PR, Huang YT, Keim P, et al. BactQuant: an enhanced broad-coverage bacterial quantitative real-time PCR assay. BMC Microbiol. 2012;12:56.
Fadrosh DW, Ma B, Gajer P, Sengamalay N, Ott S, Brotman RM, et al. An improved dual-indexing approach for multiplexed 16S rRNA gene sequencing on the Illumina MiSeq platform. Microbiome. 2014;2(1):6.
Martin M: Cutadapt removes adapter sequences from high-throughput sequencing reads. 2011 2011, 17(1):3.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.
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.
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.
We thank all the men who kindly participated in the study and staff at the Rakai Health Sciences Program for their help in study implementation.
This work was supported by the Canadian Institutes for Health Research (TMI-138656 to RK), the Fogarty HIV Research Training Program of the National Institutes of Health (4D43TW009578-04 to RMG), and RK is supported by the University of Toronto/OHTN Endowed Chair in HIV Research. This work was also supported in part by extramural funding from the National Institutes of Health (1R01AI128779 to AART) and R01A123002-01A1 (to CL). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Ethics approval and consent to participate
The Research and Ethics Board (REB) of the Uganda Virus Research Institute in Entebbe, Uganda, and the Institutional Review Board at the University of Toronto in Canada provided review, ethical oversight, and approval of the study proposal. All study participants gave written informed consent in accordance with the Declaration of Helsinki.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Galiwango, R.M., Park, D.E., Huibner, S. et al. Immune milieu and microbiome of the distal urethra in Ugandan men: impact of penile circumcision and implications for HIV susceptibility. Microbiome 10, 7 (2022). https://doi.org/10.1186/s40168-021-01185-9