Objectives Herpes simplex virus-2 (HSV-2) and syphilis are associated with increased risk of HIV, highlighting the importance of understanding their transmission dynamics. In India, most studies of HSV-2 and syphilis incidence are in high-risk populations and may not accurately reflect infectious activity. In this study, we aim to define HSV-2/syphilis incidence and risk factors in a population sample.
Methods We conducted a longitudinal population-based survey in Andhra Pradesh, India, in two rounds: 2004–2005 and 2010–2011. Sociodemographic and behavioural data were collected, and dried blood spots tested for HSV-2 and Treponema pallidum IgG. After calculating sexually transmitted infection (STI) incidence, associated factors were assessed using modified Poisson regression and within-couple transmission rates modelled using seroconcordance/discordance data.
Results 12 617 adults participated at baseline with 8494 at follow-up. Incidence of HSV-2 and syphilis per 1000 person-years was 25.6 (95% CI 24.1 to 27.2) and 3.00 (95% CI 2.52 to 3.54). Incidence of HSV-2 was higher in women vs men (31.1 vs 20.2) and in rural vs urban residents (31.1 vs 19.0) (p<0.05 for both). STI seroincidence increased in a step-wise fashion with age and was associated with spousal seropositivity for both sexes (incidence rate ratio (IRR) 2.59 to 6.78). Within couples the rate of transmission per 1000 couple-years from men to women vs women to men was higher for HSV-2 (193.3 vs 119.0) compared with syphilis (27.6 vs 198.8), p<0.05 for both.
Conclusions HSV-2 has higher incidence among subpopulations such as women, rural residents and older-aged individuals, suggesting a need for more generalised STI prevention approaches among populations traditionally considered low risk.
- EPIDEMIOLOGY (GENERAL)
This is an Open Access article distributed in accordance with the terms of the Creative Commons Attribution (CC BY 4.0) license, which permits others to distribute, remix, adapt and build upon this work, for commercial use, provided the original work is properly cited. See: http://creativecommons.org/licenses/by/4.0/
Statistics from Altmetric.com
In India, previous work has indicated a high burden of sexually transmitted infections (STIs) and a rising prevalence of genital ulcerative diseases including herpes simplex virus-2 (HSV-2) and syphilis.1 ,2 In addition to the significant morbidity resulting directly from these infections, both HSV-2 and syphilis are associated with an increased risk of HIV, and specifically a high population attributable risk for HIV in India.3–6 This highlights the need to better understand the epidemiology of these infections in India.
Existing Indian prevalence data for HSV-2 and syphilis have shown marked geographic and gendered heterogeneity. Estimates from population-based serostudies have ranged from 5.87% to 18.9% for HSV-2 and 0.4–1.7% for syphilis.7–10 These estimates have suggested higher HSV-2 prevalence in women and have shown inconsistent urban versus rural patterning. Incidence studies of these STIs in India have come from high-risk populations and may not accurately reflect disease transmission. HSV-2 incidence estimates have ranged from an annualised 8% to 15.1 per 100 person-years,4 ,11 ,12 which are relatively higher rates than estimates from Western countries.13 ,14 Syphilis incidence estimates have ranged from 5.4 per 100 person-years to an annualised incidence of 1% and 5% for men and women, respectively.5 ,12
Population-based studies have been proposed as providing more reliable estimates of true population disease rates.15 ,16 Furthermore, examinations of transmission within couples can aid in our understanding of transmission dynamics.17 Taken together such studies can help to provide a better foundation for effective HIV/STI prevention strategies.
Here we report data from a longitudinal population-based study of HIV, HSV-2 and syphilis incidence from the south Indian state of Andhra Pradesh.18 Andhra Pradesh is thought to have the highest HIV burden in India and is also reported to have high rates of HSV-2.10 ,15 In addition to calculating incidence, we also examined associated risk factors and seroconcordance/discordance patterns among stably married couples.
Study population and setting
We conducted a longitudinal population-based survey in two rounds over 5.6 years: round 1 was conducted between September 2004 and September 2005 and round 2 between May 2010 and July 2011.18 The original study design, baseline findings and HIV incidence have been described in detail elsewhere.9 ,15 ,18 ,19 In brief, at baseline, persons 15–49 years old were selected using a stratified random sampling method to construct a representative sample of the adult population of Guntur district in Andhra Pradesh. Of the 13 838 eligible people identified at baseline, 12 617 (91.2%) participated.
For round 2 data collection, interviewers returned to each household and invited subjects to participate in the follow-up survey.18 At least five attempts were made to contact each eligible participant, including tracking and contacting those that had moved. Of the 12 557 non-homeless participants at baseline, 7222 (57.5%) had remained at the same residence (including 229 who had died) and 5335 (42.5%) had moved. Of those who had moved, 3562 (66.8%) were traced and contacted.
Trained interviewers conducted confidential interviews to assess sociodemographic and behavioural characteristics.18 Sociodemographic information from the interviews included age, sex, residence (rural vs urban), education, caste, religion, marital status, occupation and standard of living index (SLI).15 Information on risk behaviour included number of lifetime sexual partners, condom use, alcohol and tobacco use, transfusion history and same-sex activity (for men).
A finger prick method was used to obtain a dried blood spot (DBS) on Whatman No. 3 paper (Whatman International, Maidstone, UK). Lifetime prevalence of HSV-2 and syphilis was evaluated by testing for the presence of HSV-2 IgG (HerpeSelect 2 ELISA IgG, Focus Diagnostics, California, USA) using an index value of >1.1,20 and Treponema pallidum IgG (Treponostika sandwich ELISA, bioMérieux, France).21 For quality control, DBSs were compared with standard serum testing of 75 samples from an STD clinic and demonstrated 100% specificity and sensitivity. Additionally, all positive DBSs were retested and 10% of negative samples were retested, yielding 100% concordance with the original results.
HSV-2 and syphilis incidence analyses
Factors associated with participation status at follow-up were analysed using χ2 and t tests. Incident cases of HSV-2/syphilis were defined as seroconversion from negative to positive in an analytic sample of those who were HSV-2/syphilis seronegative at baseline. Total person-time was calculated as the sum of the follow-up periods of each participant and was used as the denominator for the incidence calculation. Since the exact time of infection was unknown, infection was estimated to have occurred at the midway point of the study, and those with incident STIs contributed half of their total follow-up time to total person-years.22
Factors associated with STI incidence were assessed using modified Poisson regressions with robust error estimation in a multistep process.23 First, we used univariate regression to examine the association of each variable with incident infection. Then we included all variables significant at the p<0.10 level from the univariate models in a multivariable regression that was then followed by backwards elimination to select variables remaining significant at the p<0.05 level. Sociodemographic and risk multivariable models were built and examined separately in order not to underestimate the effect of factors more distally related to infection (ie, rural residence vs number of lifetime sex partners) and then combined into a final model.9 ,10 ,19 ,24 Analyses were done using STATA V.12.0 (StataCorp, Austin, Texas, USA).
HSV-2 and syphilis dyadic transmission analysis
Raw transmission rates between couples were compared using tests of proportions. In order to examine transmission dynamics, we built a Markov model with four states corresponding to the possible serostatuses of couples: (1) negatively concordant, male (M)−, female (F)−; (2) male positive discordant, M+, F−; (3) female positive discordant, M−, F+; and (4) positively concordant, M+, F+. We let α be the rate at which the man is infected outside the couple; β the rate that the woman is infected outside the couple; δ the rate of male-to-female transmission within the couple; and γ the rate of female-to-male transmission within the couple. Thus the transition rate from state 1–>2 is α; 1–>3 is β; 2–>4 is β+δ; and 3–>4 is α+γ. Then over a T=5-year period, the transition probabilities between states are P11=e−(α+β)T, P12=(e−(α+β)T−e−(β+δ)T)α/(δ−α), P13=(e−(α+β)T−e−(α+γ)T)β/(γ−β), P14=1−P11−P12−P13, P22=e−(β+δ)T, P24=1−e−(β+δ)T, P33=e−(α+γ)T, P34=1−e−(α+γ)T. Let nij be the number of couples observed in state i at baseline and state j at follow-up. The log likelihood is then ∑ij nij log(Pij). For both HSV-2 and syphilis, we then estimated the parameters α, β, δ and γ using maximum likelihood estimation. CIs were calculated using the likelihood ratio test.
Of 12 617 baseline participants, 8362 (66.3%) provided a blood spot and follow-up interview. Of this sample, 4593 (54.1%) were rural residents and 4358 (51.3%) women. Lower participation at follow-up was associated with baseline HIV/HSV-2 but not syphilis seropositivity (HSV-2(+) 63.7 vs HSV-2(−) 67.6%, p=0.01; HIV(+) 43.1 vs HIV(−) 67.8%, p<0.001; syphilis(+) 61.4 vs syphilis(–) 66.4%, p=0.12). Women and rural residents were more likely to participate (68.3% women vs 66.3% men, p=0.02) (72.7% rural vs 62.5% urban, p<0.001). Women that participated were older than women who did not (mean (SD) 30.2 (9.5) years at baseline vs 28.8 (10.0), p<0.001).
Of the 4358 women in the final sample, 389 (8.9%) were HSV-2 seropositive and 54 (1.2%) syphilis seropositive at baseline. This yielded final female analytic samples of 3969 for HSV-2 and 4304 for syphilis. Of the 4136 men in the final sample, 213 (5.1%) were HSV-2 and 84 (2.0%) syphilis seropositive at baseline and were excluded from further analysis. This yielded final male analytic samples of 3923 for HSV-2 and 4052 for syphilis.
Over a mean follow-up time of 5.63 years (range 5.17–6.33), the incidence per 1000 person-years for HSV-2 and syphilis was 25.6 (95% CI 24.1 to 27.2) and 3.00 (95% CI 2.52 to 3.54) (table 1).
The incidence of HSV-2 was higher in women than men (31.1 vs 20.2 per 1000 person-years, p<0.001) and higher in rural than urban residents (31.1 vs 19.0 per 1000 person-years, p<0.001). The incidence of syphilis in men vs women was not statistically significant (3.22 vs 2.79 per 1000 person-years, p=0.39), while rural vs urban syphilis incidence was borderline (3.44 vs 2.48 per 1000 person-years, p=0.06). Both HSV-2 and syphilis incidence increased with age, and this was consistent for both men and women (p-trend <0.05, see online supplementary figure S1).
Factors associated with infection
For HSV-2 and syphilis, the univariate and final combined multivariate analyses are presented in table 2 for men and in table 3 for women. Only factors that were included in the sociodemographic and health/behavioural multivariate models (significant at the p<0.05 level in these models) are shown here. Please refer to online supplementary tables S1 and S2 for a full reporting of all factors that were considered.
For men, sociodemographic factors associated with higher HSV-2 incidence in univariate analyses were increasing age, rural residence, transport-related occupation, lower caste, Christian or Hindu religion, lower education and marital status. In the final combined model, only age, education and Muslim religion remained significantly associated. Of the health and behavioural associations found in the univariate analyses—circumcision alcohol use, lifetime sex partners and syphilis-seropositive spouse—only spousal HSV-2 seropositivity remained statistically significant in the final combined model. Circumcision was associated with a lower incidence of HSV-2 in the risk model (incidence rate ratio (IRR) 0.74, 95% CI 0.55 to 0.97) but was excluded from the final model because of collinearity with Muslim religion.
The model for syphilis in men showed slightly different characteristics. Although age was significant in the initial sociodemographic model, once the behavioural risk factors were included in the final combined model it was no longer statistically significant at any individual strata, but remained significant taken as a whole.
When compared with the male HSV-2 data, women had a few additional factors of importance. Similar to men, older age, lower education and HSV-2-positive spouse were significant in the univariate and final models. Additionally, in both univariate and final models, lower SLI score, field jobs and previously or currently married women were shown to be at higher risk. Unlike the male data, rural residence remained significant in the final model. Although syphilis seropositivity at baseline was associated with incident HSV-2 in the univariate and initial risk behaviour analyses, this association did not remain significant in the final combined model.
For women, the only sociodemographic variables that remained significantly associated with incident syphilis were increasing age and having no education. Of the behavioural factors in the risk model—tobacco use, multiple lifetime male partners and spousal syphilis seropositivity—only spousal HSV-2 seropositivity was significant.
For both STIs, baseline HIV positivity was examined as a risk factor for subsequent STI. Although higher rates were observed, the finding did not reach statistical significance (see online supplementary tables S1 and S2). However, incident HIV was significantly more likely to occur in participants with incident HSV-2 (IRR 3.98, 95% CI 2.13 to 7.41) (data not shown).
Of the 2716 stably married couples, 2499 (HSV-2) and 2520 (syphilis) had complete data on the infection status of each partner (see online supplementary table S3). Estimating transmission probability we found that for HSV-2 infections in previously seronegative couples there was no difference between male and female acquisition of infection (19.4 vs 20.9 per 1000 couple-years, p>0.05) (table 4).
There was a slight male predominance of outside infections for syphilis (3.1 vs 1.0 per 1000 couple-years, p<0.05). Couples that were serodiscordant at baseline had markedly different transmission dynamics for HSV-2 and syphilis with estimated male-to-female vs female-to-male transmission of 193.3 vs 119.0 per 1000 couple-years for HSV-2 and 27.6 vs 198.8 for syphilis (p<0.05 for both).
To our knowledge, this is the first population-based report of HSV-2 and syphilis seroincidence from India. Overall we found seroincidence of HSV-2 and syphilis to be high relative to baseline prevalence. STI incidence was highest with HSV-2 among rural residents and specifically rural women. Two striking findings from this analysis are the consistent association of higher STI incidence with older age and rapid STI transmission among rural populations.
The pattern of higher incidence as age increases has also been reported in studies of HSV-2 seroincidence done in the USA and New Zealand, although these studies were in younger cohorts (32 and below).13 ,14 A few possible explanations of this phenomenon are as follows: (1) there is a higher seroprevalence as groups age that could confer a higher risk of partnering with someone with an active infection. One caveat to this hypothesis is that seroprevalence is not a clear marker of active infection. However, this cohort effect is made more likely by a greater lifetime risk of having multiple partners as one ages. (2) Younger couples that get married and live in joint family households (common practice in India) may gain more autonomy in their daily activities with age and seniority. They then may be able to pursue outside relationships that were not possible beforehand. (3) Public health messages on safe sex and STI prevention may have less penetrance in older populations. (4) There may be an increased biological susceptibility for HSV-2 and syphilis infection in older individuals.
Over a short 5.6-year period, a higher rural STI incidence has brought the previously lower rural vs urban prevalence of these STIs to equivalent levels. Higher prevalence in urban areas may precede a spread to the rural population. This may be mediated in part through migrant workers as they travel from cities and couple with rural partners. This hypothesis is supported by the association of incident STI with transport-related occupation that was found here and in previous work showing increased HIV risk behaviours and visits to sex workers among migrant workers.25 ,26 Now that the prevalence has ‘evened out’ in rural and urban populations, it may be that the incidence rates will also begin to equalise. Other drivers of this association may be related to sociodemographic factors such as lower levels of education and religion. A lower level of education was the largest confounder of the rural–HSV-2 association and explained some of the increases seen here. Rural residents are also more likely to be Hindu or Christian, both groups with higher STI incidence compared with Muslims. Muslims are far more likely to be circumcised, which has been shown to be protective against STIs such as HSV-2, and may explain some of the effect seen here.
Although baseline HIV was not found to be a significant risk factor for either HSV-2 or syphilis incidence, our sample size was small (just 26 men and 22 women with HIV at baseline) and we were likely underpowered to detect true differences. In an exploratory analysis, we did find a roughly 4× higher rate of incident HIV in those with incident HSV-2. This supports multiple lines of evidence showing that active HSV-2 is a significant risk factor for HIV and highlights the importance of broader STI prevention.3–6 The rates of syphilis incidence were too low to allow an examination of its association with incident HIV.
In this study, previously married women had significantly higher rates of HSV-2 and higher absolute rates of syphilis compared with women that were currently or never married. While previous work consistently found higher prevalence in this group,19 ,27 our results provide the first temporal information about this association. It suggests that for many of these women STIs are contracted after a marriage ends and not before. This indicates that previously married women may occupy a place in society that leaves them more vulnerable to unsafe sexual practices and subsequent STI acquisition.
Our examination of the serodiscordance/concordance status of couples aids in our understanding of the dynamics of STI transmission. Previous work by Arora et al examined seroconcordance among Indian couples using prevalence data and found that in contrast to studies from South Africa, men were more likely than their spouses to be the first infected partner.17 ,28 Our estimates suggest an equivalent rate of male vs female HSV-2 acquisition in seronegative couples, but higher rates of syphilis acquisition by the male partner. Similar to the Arora paper, our data support a higher risk of HSV-2 transmission from men-to-women vs women-to-men. The transmission dynamics were markedly different for syphilis, with a higher rate of transmission from women to men within couples. Past examinations of syphilis epidemiology have shown higher rates of incident cases in men, suggesting a component of differing transmission probabilities, but conclusive estimates have been elusive.29 Our data support a higher transmission probability of syphilis from women to men versus men to women. This may reflect a greater biological susceptibility in women and/or differences in sexual practices when the female versus male partner is infected.
There were several limitations to this study. The sample was limited to one district in the state of Andhra Pradesh, and especially given the geographic heterogeneity of STI prevalence in India, our findings may not be generalisable to India as a whole. A 67% follow-up rate makes it possible that the sample that remained was subject to compositional bias. In particular, those with HIV at baseline were less likely to participate in the follow-up study, indicating that some of the highest risk participants may have dropped out. However, in other cases, there were lower participation rates in individuals with factors associated with decreased STI incidence (ie, higher education).
The HSV-2 assay used has been shown to have varied performance in different populations,20 ,30 ,31 although this was taken into account when determining the index cut-offs, and quality control testing against gold standard western blot was performed. The validity of these data is also supported by the similar patterns seen in both syphilis and HSV-2 results. Lastly, HSV-2 may also be contracted orolabially leading to seroconversion. This type of transmission may be particularly prevalent among women participating in childcare. Although orolabial HSV-2 transmission would reflect true infection, this is largely clinically irrelevant in terms of HIV risk.
In conclusion, these data show a changing epidemic, in which older persons, rural populations and women are relatively more vulnerable to incident HSV-2 and syphilis. These trends suggest a need to reorient STI prevention strategies to include messaging and programmes targeted to populations previously thought to be at lower risk of infection.
Incidence of herpes simplex virus-2 (HSV-2) and syphilis over a 5-year period was high relative to baseline prevalence.
Risk of incident HSV-2 and syphilis increased with age.
Rural women and previously married women were at the highest risk of incident HSV-2.
HSV-2 transmission was more likely from men to women while syphilis transmission was more likely from women to men.
Handling editor Jackie A Cassell
Contributors LD and RD conceived and led this population-based cohort study. JAS, VL, GAK and MAlary contributed to the design. MAkbar, GMMA and SPR supervised the field data collection. VL led the laboratory analysis and TS contributed to this. CHH, JAS and BA did the statistical analysis. CHH and JAS wrote the initial draft of the paper. All authors contributed to interpretation of findings and approved the final version of the paper.
Funding This work was supported by the US National Institutes of Health under award numbers R01AI084315 and R21HD068352, and by a grant from the Bill & Melinda Gates Foundation. The content in this publication is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health or the Bill & Melinda Gates Foundation.
None of the authors have any relevant conflicts of interest to report.
Competing interests None.
Ethics approval The protocols, procedures and analyses for this study were approved by institutional ethics committees of the Public Health Foundation of India, Nizam's Institute of Medical Sciences, the Centre hospitalier affilié universitaire de Québec and the University of Chicago.
Provenance and peer review Not commissioned; externally peer reviewed.
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.