ALL Metrics
-
Views
-
Downloads
Get PDF
Get XML
Cite
Export
Track
Research Article
Revised

HIV prevalence correlated with circumcision prevalence and high-risk sexual behavior in India's states: an ecological study

[version 2; peer review: 2 approved]
PUBLISHED 17 Jun 2019
Author details Author details
OPEN PEER REVIEW
REVIEWER STATUS

Abstract

Background: HIV prevalence varies between 0% and 1.6% in India's states. The factors underpinning this variation are poorly defined.
Methods: We evaluated the relationship between HIV prevalence by state and a range of risk factors in the Indian 2015 National Family Health Survey. Pearson’s correlation was used to assess the relationship between HIV prevalence and each variable. The prevalence of each risk factor was compared between five high-HIV-prevalence states (>1% prevalence) and a large low-HIV-prevalence state (Uttar Pradesh; HIV prevalence, 0.06%).
Results: There was an association between HIV prevalence and men's mean lifetime number of partners (r = 0.55; P = 0.001) and men reporting sex with a non-married, non-cohabiting partner (r = 0.40; P = 0.014). In general, men in high-prevalence states were less likely to be circumcised and (with the exception of Chandigarh) use condoms at last sex. In two high prevalence states (Mizoram and Nagaland), men reported a higher number of lifetime partners and a higher prevalence of multiple partners and high-risk sex in the past year.
Conclusions: Variation in circumcision prevalence and sexual behavior may contribute to the large variations in HIV prevalence by state in India.

Keywords

India, HIV prevalence, high-risk sex, sexual network, ecological, circumcision

Revised Amendments from Version 1

The queries raised by the peer reviewers have been addressed as follows. Tables 2 and 3 have been merged as suggested. The legend for Figure 1 has been amended as suggested. More background information has been provided on the topics requested. Other small edits have been made to make each table and figure interpretable on its own.

See the author's detailed response to the review by John McBride
See the author's detailed response to the review by Curtis Jolly and Pauline E Jolly

Introduction

There is little consensus as to whether or not differences in sexual behavior play an important role in determining the large differences in HIV prevalence between populations. Certain authors have claimed that non-behavioral factors such as differences in prevalence of herpes simplex virus 2 infection, circumcision rates and STI treatment efficacy are responsible for differences in HIV prevalence1,2. Other authors have found that HIV prevalence is associated with high-risk sexual behavior35. These latter studies have included ecological-level analyses. The rationale behind ecological studies is that the prevalence of sexually transmitted infections (STIs) is, to an important extent, determined by the structure of the local sex network6. This is a population-level characteristic and thus ecological studies are required to assess if correlates of this network structure are associated with HIV prevalence7. A number of studies3,5,8,9, but not all studies10, have found that correlates of network connectivity, such as rate of partner change and the prevalence of sexual partner concurrency, are positively associated with HIV prevalence. These associations have been found to be strongest when comparing ethnic groups and regions within countries3,5,8,9,11,12. Nationally representative, HIV-serolinked Demographic Health Surveys (DHS) have been a particularly useful resource for these studies3,8,11.

India is an interesting country to test the network connectivity thesis. HIV was first detected in India in 1986 and since then a range of sources have documented higher HIV prevalence in certain North East States and to a lesser extent Karnataka, Andhra Pradesh and surrounding areas1317. Individual-level analyses have found that well-established risk factors, such as partner number, contact with sex workers, intravenous drug usage and lack of circumcision, were associated with HIV infection13,1720. These analyses have, however, not explained the differential HIV prevalence by state. A number of authors have argued that the higher HIV prevalence in Northeast India may be due to patterns of intravenous drug use (IVDU) in this region. Differences in contact with sex workers, patterns of migration, gender inequality and non-traditional forms of sex work, such as devadasi in Karnataka (where women are 'married' to temple deities and have sex with temple attendees), are some of the other factors that have been advanced as reasons for differences in HIV spread in India, often with little supporting evidence13,1620. In this analysis we address this issue by assessing whether there is an ecological-state-level association between various HIV risk factors and HIV prevalence.

Methods

India is a federal union comprising 29 states and 7 union territories (Figure 1). We refer to these 36 entities as states. The states are typically large and vary in surface area from Rajastan’s 342,239km2 to Goa’s 3,702km2 (median 88,752km2). The Union Territories are far smaller (median size 491km2) and for this reason are indicated with circles on the map in Figure 1. India’s vast size incorporates a wide range of ethnic and liguistic groups. It has an estimated 2000 ethnic groups and classifies 23 languages as official15. There are considerable differences between states in ethnic composition as well as in education al attainment, life expectancy and poverty rates14,15.

354c2362-29f3-4393-8a33-9ef84690aed6_figure1.gif

Figure 1. Map of India's 29 states and 7 union territories.

Reprinted from 21 under a CC BY license, with permission from Dörrbecker, original copyright 2005.

Data source

We used the 2015 National Family Health Survey (NFHS-4) for this study. The NHFS-4 received ethical committee clearance for data analyses such as the one performed here. As a result, no specific ethics committee approval was necessary for this study.

The NFHS-4 used a household-based, two-stage stratified sampling approach. The first stage selected primary sampling units (PSUs) from the 2011 National Census. PSUs were villages in rural areas and census enumeration blocks in urban areas. In the second stage of the survey, 22 households were randomly selected with systematic sampling from every selected rural and urban cluster.

A total of 98% of selected households were successfully interviewed. In the interviewed households, women aged 15–49 and men aged 15–54 were eligible to participate. The individual response rate for women was 97% and 92% for men. The response rates were high in all states except for Delhi and Chandigarh (Table 1). The behavioral questionnaire was administered in 17 local languages using computer-assisted personal interviewing. Further details of the survey are provided in Table 1 and the NFHS-4 report14.

Table 1. Prevalence of HIV and associated risk factors in 15–49-year-old women and 15–54-year-old men by state in India in 2015 (arranged according to descending HIV prevalence).

StateNAge
Men &
Womena
Age
Mena
Age
Womena
HIV
Total
N
tested
HIV
(%)
HIV
Testing
Rate
(%)
Individual
Response
Rate Men
(%)
Individual
Response
Rate
Women
(%)
Circumcision
(%)
Condom
Men (%)
Condom
Women
(%)
High
risk
sex
Men
(%)
High
risk
sex
Women
(%)
Multiple
partners
Men (%)
Multiple
partners
Women
(%)
Lifetime
sex
partners
Mene
Lifetime
sex
partners
Womene
90+
Partners
Men
(%)b
90+
Partners
Women
(%)b
% attained
no or only
primary
educationc
% in the
Poorest
Wealth
Quintiled
Mizoram266732.931.934.926071.697.195.798.32.52.99.120.00.03.80.13.41.10.10.220.25.1
Nagaland241033.132.634.121221.485.292.395.84.33.110.514.90.21.00.21.71.20.00.628.211.1
Andhra
Pradesh
252433.433.533.323301.285.585.693.810.60.61.92.80.01.10.00.91.00.00.441.13.6
Manipur293533.132.434.328891.29794.197.16.83.37.77.00.00.90.00.81.00.00.119.19.8
Chandigarh19432.931.635.51781.178.378.486.87.038.025.513.10.00.00.00.81.00.00.018.00.0
Delhi98230.029.332.66970.855.552.982.221.822.532.25.60.00.40.20.81.00.02.224.50.6
Daman and
Diu
63129.528.931.45610.786.490.994.67.29.012.94.70.00.00.00.61.10.07.520.30.0
Karnataka633332.632.832.260330.690.989.794.518.17.77.74.20.41.73.70.71.20.30.532.56.0
Telangana181332.032.031.916770.681.282.691.815.12.03.95.80.62.20.80.91.00.31.439.45.5
Maharashtra754631.931.632.371030.489.289.294.314.810.519.27.40.21.10.30.81.00.41.023.68.9
Bihar934331.330.832.191750.29796.398.416.64.47.04.10.31.50.80.91.20.31.045.242.0
Gujarat947832.532.033.690280.290.588.394.99.06.98.96.40.31.30.31.11.10.80.731.88.8
Chhattisgarh614732.131.632.960300.296.294.197.54.85.710.59.60.11.40.11.01.10.20.340.530.6
Haryana570031.330.832.155120.296.29798.97.817.418.48.70.61.90.70.91.00.20.025.11.5
Tamil Nadu855633.333.333.383380.296.496.198.516.96.26.04.40.11.50.60.81.00.20.325.22.8
Uttarakhand333331.831.133.332030.295.392.497.214.121.820.16.90.61.70.50.71.00.30.826.94.3
Punjab502432.331.633.549320.295.193.497.83.722.021.714.70.81.80.70.91.00.20.225.21.2
Madhya Pradesh1641231.931.432.7159640.19695.697.78.89.314.310.30.32.20.51.11.10.70.544.728.1
Rajasthan994931.330.632.697870.196.695.297.98.211.111.15.80.11.00.10.81.00.10.740.915.6
Assam668232.232.531.663690.192.290.196.128.26.39.11.70.20.60.20.71.00.30.534.721.9
West Bengal429732.533.131.541110.1939397.124.78.011.13.30.50.80.50.91.00.20.139.722.5
Jammu and
Kashmir
922532.531.734.189450.195.392.297.463.715.015.92.20.00.40.30.71.10.20.331.48.5
Himachal
Pradesh
371933.532.834.835480.191.784.495.44.314.821.76.10.11.10.31.11.00.10.018.61.5
Jharkhand648631.331.131.762290.192.990.195.417.45.44.84.50.31.00.50.91.10.21.440.043.7
Uttar Pradesh2100131.330.433.0205690.191.986.996.119.113.114.58.50.21.50.31.01.10.21.139.523.1
Arunachal
Pradesh
337432.632.632.731000.187.188.693.35.18.113.414.30.22.80.41.11.10.91.040.317.6
Meghalaya190431.531.032.518230.194.191.196.97.84.55.95.70.30.90.00.81.10.00.041.09.9
Sikkim137433.133.033.313640.097.497.398.13.19.318.013.10.01.00.01.21.00.10.032.90.4
Puducherry109833.733.933.310950.098.696.599.420.03.39.51.50.00.60.00.71.00.00.216.61.3
Odisha751433.133.232.971970.09391.296.84.211.413.21.80.10.70.50.91.03.81.538.532.8
Andaman
and Nicobar
islands
70534.534.834.12950.09685.394.210.47.616.29.91.72.81.51.11.00.50.023.54.3
Dadra and
Nagar Haveli
32230.529.931.86900.094.193.597.617.51.316.86.30.01.80.00.81.00.00.035.420.2
Kerala318934.033.534.82290.098.396.697.628.75.012.43.60.21.00.20.81.00.00.07.10.6
Lakshadweep24134.634.235.412820.09297.498.894.510.45.34.70.00.60.00.61.00.00.08.70.0
Goa129134.634.435.030660.09495.198.37.615.430.110.80.20.90.30.81.20.10.813.90.0
Tripura141932.432.731.913430.091.989.296.42.44.44.81.80.00.40.00.71.00.60.029.514.4

aAverage age in years by state.

b "90+ Partners" – This variable refers to the percent per state reporting more than 90 lifetime sex partners (for men and women) separately, of all respondents.

c Highest education level attained: Percent whose highest educational attainment is primary on no schooling.

d Percent in the Poorest Wealth Quintile: this variable refers to the percent of the respondents from this region that were calculated to fall in the poorest 20% (quintile) of the nationally derived wealth band. The wealth quintiles were derived from an asset index.

e Lifetime number of partners refers to the mean number of lifetime partners reported by those who reported being sexually experienced, excluding those with more than 90 partners.

Measures

HIV Testing and Prevalence estimates

In a random subsample of households, a dried blood spot from a finger-prick blood specimen was obtained from eligible women age 15–49 and men age 15–54 who consented to laboratory HIV testing. Respondents did not have access to the test results but were all referred to counseling and testing services in the local area. Dried blood samples were collected and subsequently tested for HIV using a Microlisa ELISA. Positive tests and a random 2% sample of negative samples were then tested with SD Bioline 1/2 ELISA (Standard Diagnostics Inc., Kyonggi-do, South Korea). Discordance between the two tests was resolved with Western Blot (BioRad). With the exception of Delhi, HIV testing rates were high in all states (median 95.9%, IQR 94.4-97.9%). Coverage was higher in women (95%) and men (90%) in rural areas as well as urban areas (women, 91%; men, 84%).

The NFHS-4 was designed to provide HIV prevalence estimates that are representative for the women (15–49 years) and men (15–54 years) for the whole country and rural and urban areas. It also oversampled certain states so as to be able to generate representative samples for 11 groups of states (Table 2). These were states that had been found to have higher HIV prevalence rates in previous surveys, as well as the Uttar Pradesh-centered group of states, which was chosen as a low-HIV-prevalence comparator. These 11 groups of states consisted of composites of vast populations with considerable heterogeneity in sexual behaviors and circumcision prevalence (Table 1). Our primary research question involved assessing if HIV-related risk factors differed between high- and low-HIV-prevalence states. To avoid problems related to averaging out large differences in prevalence of risk factors in the groups of states, in our primary analysis we compared risk factors with estimated HIV prevalence by state. In a sensitivity analysis this analysis was repeated using the 11 groups of states.

High HIV prevalence states. UNAIDS defines populations with a 15-49-year-old HIV prevalence of greater than 1% as a generalized HIV epidemic22. We therefore defined states with a 15-49-year-old HIV prevalence of above 1% as high-HIV-prevalence states.

Low HIV prevalence comparator state. Uttar Pradesh was chosen as the low-HIV-prevalence comparator state for a number of reasons. The previous HIV-serolinked NHFS, as well as the current NFHS survey, found an adult HIV prevalence of less than 0.1% in Uttar Pradesh15. It also has the largest population of any state in India and as a result it had the largest sample size of all states in NFHS-4. In sensitivity analyses we repeated the analyses using all Indian states with a HIV prevalence below 0.2% as the low HIV prevalence comparator population.

All HIV-prevalences were calculated for 15–49-year-olds.

Independent variables

Each of following predictor variables were calculated separately for men and women and were limited to those between the ages of 15–49 years for women and 15–54 years for men.

  • Condom use at last sex: The percentage of respondents who reported using a condom at last sex, amongst those who have had sex in the past 12 months.

  • Male circumcision: The percentage of men who reported being circumcised.

  • High-risk sex (sex with a non-married, non-cohabiting partner): The percentage of respondents who reported sex with a non-marital, non-cohabiting partner in the past 12 months, amongst all respondents who reported sex in the past 12 months.

  • Multiple partners in the past year: The percentage with two or more sexual partners in the past 12 months amongst all respondents who reported sex in the past 12 months.

  • Lifetime sex partners: The mean number of reported lifetime sex partners amongst all respondents who reported having had sex, excluding those with greater than 90 partners.

The reason for excluding those with 90 or more partners was that in a small number of states these individuals exerted a large effect on the mean number of lifetime partners. Our primary research question was whether or not there were population differences in lifetime number of partners for the majority of the population. The percent of the population with 90 partners or more was small in all states: men, median 0.17 (IQR 0-0.29); women, median 0.39 (IQR 0.2-0.89) (Table 1). As a result we found it more informative to calculate mean lifetime number of partners excluding those with more than 90 partners.

Two socioeconomic control variables were assessed:

  • Education attained: Percent of state respondents whose highest educational attainment is primary or no schooling.

  • Poverty: Percentage of the respondents from this state that were calculated to fall in the poorest 20% (quintile) of the nationally derived wealth band. The wealth quintiles were derived from an asset index.

Statistical analysis

All analyses are ecological in nature and conducted with HIV prevalence by state or group of states as the outcome variable. The analyses were conducted using STATA 13.0 (College Station, TX) and were all adjusted to account for the complex sampling strategies of the survey using the survey (SVY) command. The analyses were stratified by gender. The two-sample Wilcoxon rank-sum test was used to assess if there was a difference in the median number of lifetime partners between the high-HIV-prevalence states and Uttar Pradesh/low-HIV-prevalence states. Cohen's d was used to assess effect size. Histograms were used to depict the distribution of the number of lifetime sexual partners by gender and state. Chi-squared tests were used to assess differences in categorical variables. Pearson’s correlation was used to assess the relationship between state-level HIV prevalence and the prevalence of each risk factor.

Results

Demographic variables and HIV prevalence

An overview of the sample size, mean ages and other demographic characteristics of men and women by state are provided in Table 1. The average age of respondents in the regions ranged from 28.9 to 34.8 years in men (median 32.0 [IQR 31.1-33.0]) and 31.4 to 35.5 years in women (median 33.0 [IQR 31.4-35.5]).

The overall 2015 HIV prevalence in 15–49-year-olds was 0.25% in men and 0.23% in women. The median HIV prevalence by state was 0.13% (IQR 0.05-0.48%) and the prevalence in Uttar Pradesh was 0.06% (Table 1).

The prevalence of HIV exceeded 1% in 5 states. Three of these were in the North East (Manipur, 1.2%; Nagaland, 1.4%; and Mizoram, 1.6%) and two elsewhere (Andhra Pradesh, 1.2%; and Chandigarh, 1.1%). India’s prior HIV prevalence survey (NFHS-3) in 2005–2006 was designed to provide HIV prevalence estimates from six states15. Of these, Manipur had the highest HIV prevalence (1.13%), followed by Andra Pradesh (0.97%), Karnataka (0.69%), Maharashtra (0.62%), Tamil Nadu (0.34%) and Uttar Pradesh (0.07%) in 15–49 year old men and women.

Comparison of risk factors between high and low HIV prevalence states

Circumcision. The prevalence of circumcision was low in all states, including Uttar Pradesh (19.1%), but was significantly lower in all the high-HIV-prevalence states (2.5–10.6%; all P<0.0005).

Condoms at last sex. Men/women in Andhra Pradesh (1.9/0.6%), Manipur (7.7/3.3%), Mizoram (9.1/2.9%) and Nagaland (10.5/3.1) were less likely, but those from Chandigarh (25.5/38.0%) were more likely to report using a condom at last sex than those from Uttar Pradesh (13.5/13.1%; all P-values <0.05).

High-risk sex. The prevalence of reported high-risk sex in men was 8.3% in Uttar Pradesh. A higher proportion of men in Mizoram (20%) and Nagaland (15%) but a lower proportion in Andhra Pradesh (2.8%) reported higher risk sex (all P<0.0005). Reported proportions in women were low and did not differ between states.

Multiple partners in past year. Men in Mizoram were more likely to report multiple partners (3.8%) in the past 12 months than Uttar Pradesh (1.5%) (P<0.005).

Lifetime partners. Men/women in Mizoram (3.39/1.12) and Nagaland (1.74/1.15) reported a higher number of lifetime partners than Uttar Pradesh (0.97/1.07); women in Andhra Pradesh reported fewer partners (1.01) (all P<0.0005). The magnitude of the right shifting in the distribution of lifetime number of partners was large for men in Mizoram and Nagaland (Cohen's d = 0.68 and 0.43, respectively) and small elsewhere (Table 2; Figure 2).

354c2362-29f3-4393-8a33-9ef84690aed6_figure2.gif

Figure 2. Histograms of the reported number of lifetime sexual partners by state.

(A) men aged 15–54 years and (B) women aged 15–49 years. Data from five high HIV prevalence states and Uttar Pradesh. Distributions are truncated at 10/6 partners for men/women.

Two sensitivity analyses were performed. Firstly, repeating the analyses using all Indian states with HIV prevalences below 0.2% as the low-HIV-prevalence comparator population (instead of Uttar Pradesh) had little effect on the results (Table 2). Secondly, repeating the analyses using the 11 groups of states in place of the 36 states produced comparable results (Table 3). Low circumcision and condom-usage rates together with lifetime number of partners and high-risk sex (men only) remained significant risk factors in group 7 states (Mizoram, Manipur, and Nagaland) compared to the Uttar Pradesh group (Uttar Pradesh, Madhya Pradesh, Uttarakhand, and Rajasthan).

Table 2. Comparison of prevalence of HIV-related risk factors in high HIV prevalence states (>1%; (Andhra Pradesh, Chandigarh, (Manipur, Mizoram and Nagaland) with Uttar Pradesh (A) and all low HIV prevalence states (≤0.2%; B).

HIV
prevalence
(%)
Circumcision
(%)
High-risk sex (%)Lifetime number
of sex partners
(mean)
Multiple
partners past
year (%)
Condom use at
last sex (%)
MenWomenMenWomenMenWomenMenWomenMen
High HIV Prevalence
States
Andhra
Pradesh
1.2110.6***0.02.8***1.01**0.920.01.10.6***1.9***
Chandigarh1.17.0**0.013.1*1.000.800.00.038.0***25.5*
Manipur1.176.8***0.07.01.030.840.00.93.3***7.7***
Mizoram1.622.5***0.020.0***1.12***3.39***0.013.8**2.9***9.1***
Nagaland1.394.3***0.215.0***1.15***1.74***0.021.03.1***10.5
Low HIV Prevalence
States
A) Uttar Pradesh
(Low HIV ref)
0.0619.10.28.31.070.97Ref0.031.5
B) All states with HIV
prevalence ≤0.2%
0.1216.20.256.21.070.920.41.39.111.5

*P<0.05 **P<0.005 ***P<0.0005

Table 3. Comparison of HIV-related risk factors in 11 groups of states in NFHS-4.

P-values are for comparisons between group 11 (low HIV prevalence group) and groups 1 and 7 (high-HIV-prevalence groups).

GroupsHIV (%)Lifetime partners
(mean)
High-risk sex (%)Multiple
partners in
past year (%)
Circumcision
(%)
Condom (%)
Women
and Men
WomenMenWomenMenWomenMenMenWomenMen
Group 1:
Andhra
Pradesh and
Telangana
0.911.01***0.92**0.24.0**0.31.512.4*1.1***2.7***
Group 2: Bihar,
Jharkhand,
West Bengal,
and Andaman
& Nicobar
Islands
0.151.110.880.43.80.61.120.56.38.7
Group 3:
Gujarat, Dadra
& Nagar Haveli,
and Daman &
Diu
0.181.091.140.36.40.31.39.16.99.0
Group 4:
Himachal
Pradesh and
Jammu &
Kashmir
0.081.060.830.03.70.30.642.415.018.2
Group 5:
Karnataka
0.641.180.690.44.23.71.718.17.77.7
Group 6:
Maharashtra
and Goa
0.391.040.850.27.60.31.114.610.719.6
Group 7:
Mizoram,
Manipur, and
Nagaland
1.491.08***1.65***0.112.1***0.11.65.0***3.1***8.8***
Group 8:
Odisha and
Chhattisgarh
0.131.050.930.15.50.31.04.48.911.9
Group 9:
Punjab,
Haryana,
Delhi, and
Chandigarh
0.251.030.870.610.70.61.59.920.322.3
Group 10: Tamil
Nadu, Kerala,
Puducherry,
and
Lakshadweep
0.171.040.780.24.20.51.320.25.97.8
Group 11:
Uttar Pradesh,
Madhya Pradesh,
Uttarakhand,
and Rajasthan
0.101.060.970.28.40.31.514.311.913.9

*P<0.05 **P<0.005 ***P<0.0005.

State-level correlates of HIV prevalence

Controlling for education and poverty levels, the prevalence of higher-risk sex in men (r=0.40; P=0.017) and mean number of lifetime partners in men (r=0.55; P=0.001) were positively associated with HIV prevalence (Table 4).

Table 4. Pearson’s correlation coefficients for the association between HIV prevalence and various risk factors by region in India's 36 States in NFHS-4.

VariableMenWomen
Circumcision-0.23NA
Condom usage-0.007-0.07
High-risk sex0.40*-0.22
Multiple partners past year0.11-0.03
Mean no. of lifetime sex partners0.55**0.12
Education attained-0.14-0.14
Poverty-0.24-0.24

* P<0.05, ** P<0.005. † These two variables were calculated with the combined men and women's education/poverty data and hence are the same for men and women. NA, not applicable.

Discussion

As in other countries, the spread of HIV has been heterogenous in India. Whilst the absolute differences in HIV prevalence by state are not large, the relative differences are. A state-level HIV prevalence above 1% is four times the national average prevalence. These differences in HIV prevalence are also fairly stable based on both data from the 2005 NFHS survey15 and other data sources such as antenatal surveys from the 2000s until the present13,16. Previous studies have argued that the higher prevalence rates in certain states in Northeast India is predominantly due to patterns of IVDU13,16. We were unable to investigate the role of IVDU, but we found population-level differences in the prevalence of circumcision, condom usage and sexual behaviors, which could explain differential HIV spread.

The most consistent association we found was lower circumcision prevalence rates in the high HIV prevalence states. Circumcision rates were, however, as low or lower in a number of low HIV prevalence states as in the high HIV prevalence states. Although differences in condom usage were statistically significant, the absolute differences were small and one high-HIV-prevalence state had higher condom utilization than the low-HIV-prevalence comparators.

The relationship between the sexual behavioral risk factors and HIV prevalence was not uniform. Whilst sexual behaviors were in general riskier in the North East states, this was not the case in Andhra Pradesh and Chandigarh, where a number of the behavioural risk factors were actually less prevalent than Uttar Pradesh. In a country as vast and diverse as India, it should not be too surprising to find different combinations of risk factors to be responsible for HIV spread in different regions. As noted above, previous studies have suggested that various factors in Andhra Pradesh and Karnataka may play a role in the local HIV epidemics here13,20. An important finding of our study was the higher number of lifetime partners, multiple partners in the prior year and high-risk sex in Mizoram and Nagaland compared to Uttar Pradesh. Our study is thus compatible with the thesis that these behaviors would translate into a more connected sexual network which could play a role in the generation of higher HIV prevalence rates in these states.

These same three risk factors (number of lifetime partners, multiple partners in the prior year and high-risk sex) have previously been found to be associated with differences in HIV prevalence by ethnic group in Kenya3, South Africa5 and elsewhere8,23. Previous reports from DHS data found that, at an individual level, there was a stepwise increase in HIV prevalence with increasing lifetime sex partners in all 15 countries with available data. This was true for both men and women in all cases4,24. This association was also present in the Indian NFHS-3 in 2005 and in this survey14,15. That this association was present at both individual and population levels increases the probability that the association between lifetime sex partners and HIV prevalence is real25,26. In our study, this association was, however, only statistically significant for men. High-risk sex has also been found to be an individual level risk factor for HIV infection in a number of a number of countries, including India4,14. As in other areas, it is likely that these patterns of higher risk sexual behavior interact with other risk factors such as IVDU to produce the observed differences in HIV prevalence16,19.

This study has a number of limitations. The study is ecological in nature and is thus susceptible to the ecological inference fallacy. DHS surveys are not optimal for collecting sensitive sexual information27,28. The use of computer-assisted interviewing would be expected to have reduced but not eliminated this problem27. In addition, although the response rates for participation in the survey and HIV testing were high, this varied somewhat between states. The survey was designed to provide HIV prevalence estimates for 11 groups of states and not individual states. The data is thus susceptible to a large number of biases such as misclassification, nonresponse, recall and social desirability biases. In particular, other work has found evidence of culture-specific heterogeneity in answering questionnaires29,30. We cannot exclude the possibility that respondents from states where lower-risk sexual behavior was reported were subject to a greater social desirability bias, which could invalidate our findings. The available evidence, however, suggests that only minor differences in sexual behavior exist between those who do and do not answer sexual behavior questionnaires31. Finally, the analyses do not control for the influence of other variables.

Conclusion

We found a range of risk factors to be more prevalent in high-HIV-prevalence states in India. There was no clear single risk factor (or combination thereof) which appeared capable of explaining the heterogeneity of HIV spread in India. In the case of Mizoram and Nagaland, however, a higher prevalence of sexual risk behaviors may be contributing to higher HIV prevalence rates in this region. More detailed comparative studies between these populations and lower prevalence populations elsewhere in India may confirm or refute this finding.

Studies of other higher-HIV-prevalence-populations that have managed to reduce HIV incidence, including Uganda, Zimbabwe, Kenya and Thailand, have pointed to the importance of reductions in multiple partnering in effecting this decline3235. If other studies confirm our findings, then similar campaigns could be considered in Mizoram and Nagaland.

Data availability

The NFHS-4 survey is freely available from www.measureDHS.com as part of the India: Standard DHS, 2015–16. Access to the dataset requires registration, and is granted to those that wish to use the data for legitimate research purposes. A guide for how to apply for dataset access is available at: https://dhsprogram.com/data/Access-Instructions.cfm.

Comments on this article Comments (0)

Version 2
VERSION 2 PUBLISHED 15 Jan 2019
Comment
Author details Author details
Competing interests
Grant information
Copyright
Download
 
Export To
metrics
Views Downloads
F1000Research - -
PubMed Central
Data from PMC are received and updated monthly.
- -
Citations
CITE
how to cite this article
Kenyon CR. HIV prevalence correlated with circumcision prevalence and high-risk sexual behavior in India's states: an ecological study [version 2; peer review: 2 approved] F1000Research 2019, 8:60 (https://doi.org/10.12688/f1000research.17807.2)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
track
receive updates on this article
Track an article to receive email alerts on any updates to this article.

Open Peer Review

Current Reviewer Status: ?
Key to Reviewer Statuses VIEW
ApprovedThe paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approvedFundamental flaws in the paper seriously undermine the findings and conclusions
Version 2
VERSION 2
PUBLISHED 17 Jun 2019
Revised
Views
1
Cite
Reviewer Report 04 Jul 2019
Pauline E Jolly, Department of Epidemiology, University of Alabama at Birmingham, Birmingham, AL, USA 
Approved
VIEWS 1
I have gone over the paper and the revisions have been made. The paper is now acceptable for indexing. ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Jolly PE. Reviewer Report For: HIV prevalence correlated with circumcision prevalence and high-risk sexual behavior in India's states: an ecological study [version 2; peer review: 2 approved]. F1000Research 2019, 8:60 (https://doi.org/10.5256/f1000research.21562.r49983)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
Version 1
VERSION 1
PUBLISHED 15 Jan 2019
Views
9
Cite
Reviewer Report 10 Jun 2019
John McBride, College of Medicine and Dentistry, James Cook University, Cairns, Australia 
Approved
VIEWS 9
The paper provides some valuable insights into the transmission of HIV in India. The findings are consistent with other international studies in showing that an increased number of lifetime sexual partners and lower condom use is associated with higher HIV ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
McBride J. Reviewer Report For: HIV prevalence correlated with circumcision prevalence and high-risk sexual behavior in India's states: an ecological study [version 2; peer review: 2 approved]. F1000Research 2019, 8:60 (https://doi.org/10.5256/f1000research.19468.r48166)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
  • Author Response 17 Jun 2019
    Chris Kenyon, University of Cape Town, Cape Town, 7925, South Africa
    17 Jun 2019
    Author Response
    1. The wording as to the ethical approval process has been amended as suggested.
    2. The following text has been added to the methods section to assist the reader
    ... Continue reading
COMMENTS ON THIS REPORT
  • Author Response 17 Jun 2019
    Chris Kenyon, University of Cape Town, Cape Town, 7925, South Africa
    17 Jun 2019
    Author Response
    1. The wording as to the ethical approval process has been amended as suggested.
    2. The following text has been added to the methods section to assist the reader
    ... Continue reading
Views
10
Cite
Reviewer Report 12 Mar 2019
Curtis Jolly, Department of Agricultural Economics, Auburn University, Auburn, AL, USA 
Pauline E Jolly, Department of Epidemiology, University of Alabama at Birmingham, Birmingham, AL, USA 
Approved with Reservations
VIEWS 10
  1. The study is simple and is a good start for a more complex study with more in-depth analyses.
  2. The author should describe the ecological zones in more detail. What makes zones homogenous or heterogeneous (poverty,
... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Jolly C and Jolly PE. Reviewer Report For: HIV prevalence correlated with circumcision prevalence and high-risk sexual behavior in India's states: an ecological study [version 2; peer review: 2 approved]. F1000Research 2019, 8:60 (https://doi.org/10.5256/f1000research.19468.r45246)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
  • Author Response 17 Jun 2019
    Chris Kenyon, University of Cape Town, Cape Town, 7925, South Africa
    17 Jun 2019
    Author Response
    1. Thank you for this comment.
    2. We have provided more background information as the heterogeneity of peoples in India’s states (Page 3):
      India is a federal union comprising 29
    ... Continue reading
COMMENTS ON THIS REPORT
  • Author Response 17 Jun 2019
    Chris Kenyon, University of Cape Town, Cape Town, 7925, South Africa
    17 Jun 2019
    Author Response
    1. Thank you for this comment.
    2. We have provided more background information as the heterogeneity of peoples in India’s states (Page 3):
      India is a federal union comprising 29
    ... Continue reading

Comments on this article Comments (0)

Version 2
VERSION 2 PUBLISHED 15 Jan 2019
Comment
Alongside their report, reviewers assign a status to the article:
Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations - A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions
Sign In
If you've forgotten your password, please enter your email address below and we'll send you instructions on how to reset your password.

The email address should be the one you originally registered with F1000.

Email address not valid, please try again

You registered with F1000 via Google, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Google account password, please click here.

You registered with F1000 via Facebook, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Facebook account password, please click here.

Code not correct, please try again
Email us for further assistance.
Server error, please try again.