Multilevel Analysis of the Predictors of HIV Prevalence among Pregnant Women Enrolled in Annual HIV Sentinel Surveillance in Four States in Southern India

Abstract

Background

Heterogeneity of the HIV epidemic across districts of south India is reflected in HIV positivity among antenatal clinic (ANC) attendees. Along with individual factors, contextual factors also need consideration for effective HIV interventions. Thus, identifying district and individual level factors that influence ANC HIV positivity assumes importance to intervene effectively.

Methods

Data on HIV sentinel surveillance among the ANC population were obtained from the National AIDS Control Organization (NACO) between years 2004 and 2007. Data from serial cross-sectional studies among female sex workers (FSWs) conducted during this time period in 24 districts were used to generate district level variables corresponding to parameters concerning this high risk population. Other district level data were obtained from various official/governmental agencies. Multilevel logistic regression was used to identify individual and district level factors associated with ANC-HIV positivity.

Results

The average ANC-HIV prevalence from 2004 to 2007 in the 24 integrated biological and behavioural assessments (IBBA) districts ranged from 0.25 to 3.25%. HIV positivity was significantly higher among ANC women with age≥25 years [adjusted odds ratio (AOR):1.49; 95% confidence interval (95%CI):1.27 to 1.76] compared to those with age<25 years; illiterate (AOR:1.62; 95%CI:1.03 to 2.54) compared to literate; employed in agriculture (AOR:1.34; 95%CI:1.11 to 1.62) or with occupations like driver/helper/industry/factory workers/hotel staff (AOR:1.59; 95%CI:1.26 to 2.01) compared to unemployed. District level HIV prevalence among FSWs (AOR:1.03; 95%CI:1.0 to 1.05) and percentage women marrying under 18 years were significantly associated with ANC-HIV positivity (AOR:1.02; 95%CI:1.00 to 1.04).

Conclusion

Illiteracy of the woman, higher HIV prevalence among FSWs and early marriage were associated with HIV positivity among pregnant women in southern India. In addition to targeted HIV preventive interventions among FSWs, studying and changing the behavior of FSW clients and addressing structural drivers of the epidemic might indirectly help reduce HIV infection among women in southern India.

Abstract

Background

Heterogeneity of the HIV epidemic across districts of south India is reflected in HIV positivity among antenatal clinic (ANC) attendees. Along with individual factors, contextual factors also need consideration for effective HIV interventions. Thus, identifying district and individual level factors that influence ANC HIV positivity assumes importance to intervene effectively.

Methods

Data on HIV sentinel surveillance among the ANC population were obtained from the National AIDS Control Organization (NACO) between years 2004 and 2007. Data from serial cross-sectional studies among female sex workers (FSWs) conducted during this time period in 24 districts were used to generate district level variables corresponding to parameters concerning this high risk population. Other district level data were obtained from various official/governmental agencies. Multilevel logistic regression was used to identify individual and district level factors associated with ANC-HIV positivity.

Results

The average ANC-HIV prevalence from 2004 to 2007 in the 24 integrated biological and behavioural assessments (IBBA) districts ranged from 0.25 to 3.25%. HIV positivity was significantly higher among ANC women with age≥25 years [adjusted odds ratio (AOR):1.49; 95% confidence interval (95%CI):1.27 to 1.76] compared to those with age<25 years; illiterate (AOR:1.62; 95%CI:1.03 to 2.54) compared to literate; employed in agriculture (AOR:1.34; 95%CI:1.11 to 1.62) or with occupations like driver/helper/industry/factory workers/hotel staff (AOR:1.59; 95%CI:1.26 to 2.01) compared to unemployed. District level HIV prevalence among FSWs (AOR:1.03; 95%CI:1.0 to 1.05) and percentage women marrying under 18 years were significantly associated with ANC-HIV positivity (AOR:1.02; 95%CI:1.00 to 1.04).

Conclusion

Illiteracy of the woman, higher HIV prevalence among FSWs and early marriage were associated with HIV positivity among pregnant women in southern India. In addition to targeted HIV preventive interventions among FSWs, studying and changing the behavior of FSW clients and addressing structural drivers of the epidemic might indirectly help reduce HIV infection among women in southern India.

Figures

Table 3

Table 4

Fig 1

Table 1

Table 2

Table 3

Table 4

Fig 1

Table 1

Table 2

Citation: Thamattoor U, Thomas T, Banandur P, S R, Duchesne T, Abdous B, et al. (2015) Multilevel Analysis of the Predictors of HIV Prevalence among Pregnant Women Enrolled in Annual HIV Sentinel Surveillance in Four States in Southern India. PLoS ONE 10(7): e0131629. https://doi.org/10.1371/journal.pone.0131629

Editor: Xi Luo, Brown University, UNITED STATES

Received: January 18, 2014; Accepted: June 4, 2015; Published: July 6, 2015

Copyright: © 2015 Thamattoor et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Funding: This research was funded by the Bill & Melinda Gates Foundation. The views expressed herein are those of the authors and do not necessarily reflect the official policy or position of the Bill & Melinda Gates Foundation.

Competing interests: The authors have declared that no competing interests exist.

Background

India is the second most populous country in the world and there is an estimated 2.3 million people living with HIV/AIDS in India [1, 2].The HIV epidemic in India is heterogeneous, both within and between districts in the four high prevalence southern Indian states, namely Andhra Pradesh, Karnataka[3], Tamil Nadu and Maharashtra [4, 5]. HIV transmission in South India is mainly heterosexual. Over 80% of HIV-infected women in the general population acquire the infection from their husbands who buy sex or have sexual intimate partners other than wives [6]. During 2007, HIV sentinel surveillance was conducted at 646 antenatal clinics, and samples were collected from 245,516 pregnant women throughout the country [1]. HIV prevalence among antenatal clinic attendees (ANC) in the four southern states was found to be five times more than in the rest of the country. An ecological study on district level high-risk population variables has shown an association between HIV prevalence among female sex workers (FSWs) and ANC HIV prevalence, which was considered as a proxy for general population HIV prevalence in southern India [4]. Another independent study on south Indian pregnant women showed that individual level characteristics such as illiteracy and being employed but not in a service oriented job could also be associated with HIV risk [7]. Hence it is important to simultaneously examine the influence of district level as well as individual characteristics on HIV risk in this population. In addition, the associations previously identified in the published ecological analysis [4] could be spurious because of ecological bias and lack of appropriate control for confounding [8].

In India, since early 2004, a comprehensive HIV prevention programme, namely Avahan (a Sanskrit word meaning “a call to action”), the India AIDS Initiative of the Bill & Melinda Gates Foundation, has been operational in the six Indian states most affected by the HIV epidemic[9]. Cross-sectional studies, known as integrated biological and behavioural assessments (IBBA), were conducted over a 19 month period between November 2005 and June 2007 across 29 districts in India where the Avahan program for high risk groups had been implemented: Andhra Pradesh, Karnataka, Maharashtra, Tamil Nadu, Manipur and Nagaland (the latter two states are located in the North-East of the country where the HIV epidemic is driven by injection drug use), and among four segments of the National Highways, to collect data on the prevalence of HIV and sexually transmitted infections. Data on HIV risk behaviours and exposure to intervention programs in a total of over 25,000 female sex workers (FSWs), men who have sex with men (MSM), transgender, injection drug users and other bridge groups such as clients of FSW and truck drivers [10] were collected. Therefore, with the aim of validating the ecological association previously identified between ANC HIV positivity and the level of HIV prevalence among FSWs, we conducted a study to assess the association of individual and population level variables with ANC HIV positivity in the 24 IBBA districts of the four southern states, using a multilevel modelling approach.

Methods

Ethical considerations: The study was approved by the ethics committees of all institutes that were involved in the data collection for this study: National AIDS Control Organization, Delhi, the National AIDS Research Institute, Pune, the National Institute of Epidemiology, Chennai (Tamil Nadu), the National Institute of Nutrition, Hyderabad (Andhra Pradesh), and St. John’s Medical College, Bangalore (Karnataka), India, as well as Family Health International, Arlington, VA, USA, and the University of Manitoba, Winnipeg, Canada. Finally, regulatory approval for the conduct of the IBBA and its protocols was obtained from the Health Ministry Screening Committee of Indian Council of Medical Research, Government of India.

Data on individual level factors were collected from annual HIV sentinel surveillance among the ANC population (ANC HSS) conducted by the National AIDS Control Organization (NACO) [1], Government of India, for the years 2004 to 2007. Data on individual level factors included age of the respondent at the time of interview, education (illiterate, up to grade 5, grade 5–12, graduation and above), migrant status (yes, no), locality (rural, urban) and occupation (agriculture, business owner, service, truck/auto/taxi driver/helper/industry/factory workers/hotel staff, unemployed). These data were obtained on request from NACO.

In the four southern states, the IBBA was carried out among FSWs in 24 districts. The response rate of IBBA among FSWs ranged from 44% to 90% across different districts[11]. The IBBA data gave rise to numerous peer-reviewed publications [4, 11, 12]. Avahan developed a computerized management information system (CMIS) data (2005–2009) during the course of implementation of its program and data on program inputs, infrastructure, outreach, and clinical service utilization was developed. Several indicators on FSW in CMIS data were validated with IBBA data for Maharashtra and Tamil Nadu [13, 14].

Data on district level variables on high risk groups were obtained from the first round of FSW IBBA conducted between 2004 and 2007. The IBBA data were used to compute HIV prevalence in FSWs and mean number of clients of FSWs for each of the districts. These two variables were the only ones included from the FSW IBBA data, because in a previous study, they were identified as the only significant predictors of ANC HIV prevalence out of a large number of data extracted from the IBBA (including the IBBA carried out among MSM and clients of FSWs) [4]. Since Avahan is essentially an urban intervention, we only used the ANC HIV prevalence data in urban areas, except for the district of Belgaum, Karnataka state, where the Avahan intervention covered both urban and rural areas. These data were obtained on request from the National AIDS Research Institute [10]. For all the 24 IBBA districts of the four southern Indian states, data on population level variables were collected from different sources such as the census of India, IndiaStat.Com website, district level household and facility website, gateway to districts of India website etc. (see Additional file 1: S1 Table.).

A total of 49 district level variables were hypothesised to be affecting HIV prevalence at the district level including the two high risk variables, mean number of clients reported by FSWs and HIV prevalence in FSWs. All these data used in this study are either publicly available or on request and the relevant links are provided in Additional file 1: S1 Table.

Statistical Analysis

Multilevel logistic regression analysis was performed with HIV positivity from ANC HSS data as the outcome variable, and individual (level 1) and district (level 2) variables included as independent variables for the multilevel modelling analysis. The algorithm for inclusion of variables in multilevel modelling is shown in Fig 1.

[Figure omitted. See PDF.]

Fig 1. Schematic representation for selection of variables for multilevel modelling analysis.

https://doi.org/10.1371/journal.pone.0131629.g001

Selection of Individual variables

Individual level variables included in the final multilevel model were identified using logistic regression analysis. Individual level variables significantly associated with ANC HIV positivity (p < 0.01) in the univariate analysis were considered for the multivariate model. Only those variables that were significant at p < 0.05 level in the multivariate model were kept in the final multilevel logistic regression model.

Selection of district level variables

All 49 district level variables were individually checked for their association with the mean of the yearly prevalence of ANC HIV from 2004 to 2007 using simple linear regression. The district level variables significant at p<0.05 were included in the multiple linear regression model. The variables which were significant at p<0.05 in the multiple linear regression analysis were eligible to be considered in the multilevel logistic regression model (Fig 1). In addition, the high risk group variables (mean number of clients in the last week reported by FSWs and HIV prevalence among FSWs), were included in the multiple linear regression a priori. The selection of individual and district level variables was performed using SAS 9.2 for Windows version 7 (SAS Institute Inc., Cary, NC, USA).

Analyses using multilevel modelling techniques

The multilevel modelling analysis was performed using STATA IC 11.1 for Windows (Statacorp LP, Texas, USA). A random intercept logistic model was used to determine factors that were associated with inter-district variations in HIV positivity in pregnant women by fitting a two level model (individuals at level 1 nested within district at level 2). Different models were constructed and compared. A null or unconditional model without any exposure variables was specified to decompose the amount of variance that existed between districts. In the next model, individual-level variables (age, educational and employment status) were included. The model was further extended to include the district level variables that were statistically significant in the multiple linear regression.

Results

District level factors associated with ANC HIV Prevalence

ANC HIV prevalence ranged from 0.25% in Chennai district to 3.25% in Belgaum district (Table 1). The percentage of women marrying under 18 years was highly variable between the districts, with the lowest percentages in Coimbatore (4%) and Hyderabad (5%), and the highest in Belgaum (43%). Overall, the average male and female literacy rates were 79% (range: 67% to 91%) and 61% (range: 43% to 81%), respectively. A total of four out of 49 district level variables were significantly associated with ANC HIV prevalence in simple linear regression models at p<0.05 level. For the non-IBBA variables, these were mean age at marriage for girls which was negatively associated with HIV prevalence (p = 0.014), and percentage of women marrying under 18 years which was positively associated with HIV prevalence (p = 0.014) (see Additional file 1: S1 Table). Inter-district variation of HIV in FSWs ranged from 2% to 38%. The mean number of clients in the last week reported by FSWs (p = 0.036) and HIV prevalence among FSWs (p = 0.002) were also positively associated with ANC HIV prevalence in the univariate linear regression model. In the multiple linear regression model, HIV prevalence among FSWs and percentage of women marrying under 18 years were significant at p<0.05 such that for a percentage increase in HIV prevalence among FSWs, HIV ANC prevalence was higher by about 0.01 percentage (β = 0.01, p = 0.004) and a percentage increase in the women marrying under 18 years, HIV ANC prevalence was higher by about 0.05 percentage (β = 0.05, p = 0.027).

[Figure omitted. See PDF.]

Table 1. Selected characteristics of districts where integrated biological and behavioral assessment (IBBA) was conducted.

https://doi.org/10.1371/journal.pone.0131629.t001

Individual factors affecting ANC HIV positivity

The overall average percentage of ANC women aged <25 years was 70% (range: 59, 87). Table 2 shows that 26% of ANC women were illiterate, 35% of the women were employed in the agriculture sector and 41% were unemployed. The results of univariate logistic regression analyses of HIV positivity are presented in Table 3. All the individual level variables except locality and migrant status were considered in the multiple logistic regression model.

[Figure omitted. See PDF.]

Table 2. Individual level characteristics of ANC* women in IBBA districts; Annual HIV Sentinel Surveillance 2004–2007.

https://doi.org/10.1371/journal.pone.0131629.t002

[Figure omitted. See PDF.]

Table 3. Association of individual level variables with HIV positivity in antenatal women.

https://doi.org/10.1371/journal.pone.0131629.t003

In the multiple logistic regression model (Table 3), significantly higher odds of HIV positivity was seen among ANC women aged 25 years or more (AOR: 1.38;95% CI:1.17 to 1.61) compared to those aged below 25 years, women employed in agriculture (AOR: 1.4;95% CI:1.16 to 1.69) and women employed in sectors truck/auto/taxi driver/helper/industry/factory workers/hotel staff (AOR: 1.61;95% CI:1.28 to 2.03) compared to those unemployed. HIV positivity was also significantly higher among illiterate women (AOR: 2.71; 95% CI: 1.76 to 4.18), women with education below grade 5 (AOR: 2.06; 95% CI: 1.33 to3.19) and those with only grade 5–12 education (AOR: 1.84; 95% CI: 1.21 to 2.81) compared to women who had graduated from school.

In the final multilevel logistic regression model, ANC women aged 25 years or more (AOR: 1.49;95% CI:1.27 to 1.76) compared to those aged below 25 years, illiterate women (AOR: 1.62;95% CI:1.03 to 2.54) compared to women who had graduated from school and being employed in agriculture (AOR: 1.34;95% CI:1.11 to 1.62), and in occupations truck/auto/taxi driver/helper/industry/factory workers/hotel staff (AOR: 1.59;95% CI:1.26 to 2.01) compared to those unemployed, were significantly associated with ANC HIV positivity. Concerning the district level variables, the odds of HIV infection was 3% higher for ANC women for every percent increase in HIV prevalence among FSWs (AOR: 1.03; 95% CI: 1.01 to 1.05). Furthermore, the odds of HIV infection was 2% higher for ANC women for each percent increase in women marrying less than 18 years (AOR: 1.02; 95% CI: 1.00 to 1.04) in the district (Table 4).

[Figure omitted. See PDF.]

Table 4. Multilevel analysis of factors associated with HIV in antenatal women: ANC HSS 2004–2007.

https://doi.org/10.1371/journal.pone.0131629.t004

Discussion

In this study, individual and district level variables that could characterize HIV positivity among pregnant women across 24 districts in four southern Indian states were examined using a multilevel statistical modelling approach. Among individual level characteristics, older women, illiterate women, those employed in agriculture and occupations truck/auto/taxi driver/helper/industry/factory workers/hotel staff had significantly higher odds of HIV positivity. Among district level variables, HIV prevalence among FSWs and the percentage of women marrying below 18 years were significantly associated with a higher HIV positivity among pregnant women.

This study strengthens the evidence that HIV in the FSW population is a determinant of ANC HIV positivity and thus confirms the results of a previous ecological analysis [4]. This suggests that in India, women in the general population are probably getting infected with HIV from their partners who are clients of FSWs, as strongly suggested by previous mathematical modelling based on general and high risk population data from India [6]. Studies have shown that among all HIV cases, heterosexual transmission contributed to 87% and mother-to-child transmission contributed to 5.4%. Furthermore, IDU and homosexual transmission accounted for 1.6 and 1.5% of all HIV cases, respectively. Thus, primarily heterosexual transmission from the core group (FSWs) to the general population mostly happens through the clients of FSWs to their spouses [15]. Higher age has been shown to be associated with a higher risk of HIV infection in other studies [16]. The greater risk of HIV among older women could be attributed to longer exposure to sexual activity [17–20]. This also could be due to the chronic nature of HIV infection such that women who were infected with HIV at younger age could contribute to greater prevalence in older women, as their age progresses. The higher risk of HIV among illiterate women, those employed in agriculture, and in occupations such as truck/auto/taxi driver/helper/industry/factory workers/hotel staff, suggests that women with low socio-economic status are at higher risk of HIV. The lower HIV levels found in unemployed women may be related to the fact that, in the Indian context, most of them are housewives, with more stable sexual relationships than employed women.

Early sexual activity, of which age at marriage is a proxy in the Indian context, is a known risk for HIV both in India and globally. This is partly due to the biological vulnerability of young women due to the sensitive nature of their genital tract[21]. In a study by Bhattacharya, 97% of women surveyed in India in 1992–1993 did not use any contraception before their first child was born [22]. Research conducted in Kenya and Zambia shows that young married girls are more likely to be HIV positive than their unmarried peers [23–26]. While child marriage has decreased globally over the last 30 years, it remains common in rural areas [24]. Child marriage is most common in sub-Saharan Africa and South Asia (where 42 and 48 per cent of girls, respectively, marry before age 18) [24, 27], and we observed a similar high percentage of women marrying under 18 years in some districts such as Belgaum (43%) where ANC HIV was also the highest among all IBBA districts (3.3%). Poverty, security, and family status are some of the reasons for the continued practice of child marriage in India.

A limitation of this study is that the various district level variables considered for the model were collected from the latest available reliable data sources, which encompassed different time points (2001 to 2010). However, as the outcome variable was obtained from four years of ANC HIV surveillance (2004–2007), the time interval for district level data was not unreasonable. Another limitation of the study is that except Belgaum district,

all other districts are urban localities. As multiple statistical tests were used to screen the variables (at both district and individual levels) to be included in the final model, another potential limitation is an increased likelihood of having included in the final model variables that are not truly associated with HIV in the whole population of women attending antenatal clinics.In conclusion, HIV prevalence among FSWs is associated with HIV positivity among pregnant women in southern India. Illiteracy and lower age at marriage of women were also associated with HIV positivity. These structural factors may increase women’s vulnerability to HIV infection and interventions to increase literacy and increase age at marriage could have an indirect positive impact on HIV among women. In addition to targeted HIV preventive interventions among FSWs, studying and changing the behaviour of FSW clients could help reduce HIV infection among women in southern India.

Supporting Information

[Figure omitted. See PDF.]

S1 Table. List of variables used in the univariate linear regression.

https://doi.org/10.1371/journal.pone.0131629.s001

(DOC)

Acknowledgments

The authors would like to thank the Avahan lead implementing partners, Family Health International, Hindustan Latex Family Planning Promotion Trust (HLFPPT), the International HIV/AIDS Alliance, the Karnataka Health Promotion Trust (KHPT), Pathfinder International, and Voluntary Health Services, for providing IBBA-related data. We also thank the National AIDS Control Organisation (NACO) in India, and the Karnataka State AIDS Prevention Society (KSAPS) for their support. This research was funded by the Bill & Melinda Gates Foundation. The views expressed herein are those of the authors and do not necessarily reflect the official policy or position of the Bill and Melinda Gates Foundation.

Author Contributions

Conceived and designed the experiments: MA RW RBM SM. Analyzed the data: UT TT PB RS TD BA MA. Wrote the paper: UT TT PB RS TD BA MA RW RBM SM.

Citation: Thamattoor U, Thomas T, Banandur P, S R, Duchesne T, Abdous B, et al. (2015) Multilevel Analysis of the Predictors of HIV Prevalence among Pregnant Women Enrolled in Annual HIV Sentinel Surveillance in Four States in Southern India. PLoS ONE 10(7): e0131629. https://doi.org/10.1371/journal.pone.0131629

References

1. National Control Organisation (NACO).HIV Sentinel Surveillance and HIV Estimation in India 2007. A technical Brief 2007. Available: http://www.nacoonline.org/upload/Publication/M&E%20Surveillance,%20Research/HIV%20Sentinel%20Surveillance%20and%20HIV%20Estimation%202007_A%20Technical%20Brief.pdf.

2. UNAIDS.Asia: AIDS epidemic summary: regional summary 2008. Available: http://data.unaids.org/pub/Report/2008/jc1527_epibriefs_asia_en.pdf.

3. Banandur P, Rajaram SP, Mahagaonkar SB, Bradley J, Ramesh BM, Washington RG, et al. Heterogeneity of the HIV epidemic in the general population of Karnataka state, south India. BMC Public Health. 2011;11 Suppl 6:S13. Epub 2012/03/02. doi: 1471-2458-11-S6-S13 [pii] pmid:22376218; PubMed Central PMCID: PMC3287551.

4. Alary M, Jayachandran AA, Lowndes CM, Bradley J, Demers E, Adhikary R, et al. Ecological analysis of the association between high-risk population parameters and HIV prevalence among pregnant women enrolled in sentinel surveillance in four Southern India states. Sex Transm Infect. 2010;86 Suppl 1:i10–6. Epub 2010/03/05. doi: 86/Suppl_1/i10 [pii] pmid:20167724; PubMed Central PMCID: PMC3252620.

5. Ng M, Gakidou E, Levin-Rector A, Khera A, Murray CJL, Dandona L. Assessment of population-level effect of Avahan, an HIV-prevention initiative in India. The Lancet. 378(9803):1643–52.

6. Vickerman P, Foss AM, Pickles M, Deering K, Verma S, Eric D, et al. To what extent is the HIV epidemic in southern India driven by commercial sex? A modelling analysis. AIDS. 2010;24(16):2563–72. Epub 2010/09/21. pmid:20852405.

7. Solomon E, Visnegarwala F, Philip P, Alexander G. Predictors of HIV positivity among pregnant women presenting for obstetric care in South India—a case-control study. AIDS Care. 2011.

8. Morgenstern H. Ecologic studies. In: Rothman KJ GS, Lash TL, eds. Modern epidemiology. Philadelphia: Lippincott-Williams & Wilkins. 2008;(3rd edn):511–31.

9. AVAHAN-THE INDIA AIDS INITIATIVE: The Business of HIV Prevention at Scale New Delhi, India: Bill & Melinda Gates Foundation; 2008. Available: http://www.gatesfoundation.org/avahan/Documents/Avahan_HIVPrevention.pdf.

10. Details Regarding IBBA Round-I & II Data. Access for Scientists & Researchers.: National Aids Research Institute. Available: http://www.nari-icmr.res.in/IBBAdataaccess.php.

11. Ramesh B, Moses S, Washington R, Isac S, Mohapatra B, Mahagaonkar S, et al. Determinants of HIV prevalence among female sex workers in four south Indian states: analysis of cross-sectional surveys in twenty-three districts. AIDS 2008;22(5):10.

12. Kumar GA, Dandona R, Alary M, Dandona L. Predictors of HIV prevalence among street-based female sex workers in Andhra Pradesh state of India: a district-level analysis. BMC Women's Health. 2014;14(1):65.

13. Mainkar MM, Pardeshi DB, Dale J, Deshpande S, Khazi S, Gautam A, et al. Targeted interventions of the Avahan program and their association with intermediate outcomes among female sex workers in Maharashtra, India. BMC public health. 2011;11(Suppl 6):S2. pmid:22375562

14. Thilakavathi S, Boopathi K, Kumar CG, Santhakumar A, Senthilkumar R, Eswaramurthy C, et al. Assessment of the scale, coverage and outcomes of the Avahan HIV prevention program for female sex workers in Tamil Nadu, India: is there evidence of an effect? BMC public health. 2011;11(Suppl 6):S3. pmid:22375609

15. UNAIDS. UNGASS Country Progress Report 2010: India. 2010.

16. Dandona R, Dandona L, Gutierrez J, Kumar A, McPherson S, Samuels F, et al. The ASCI FPP Study Team. High Risk of HIV in non-brothel based female sex workers in India. BMC Public Health. 2005;5(87).

17. Use it or lose it: How Avahan used data to shape its HIV prevention efforts in India New Delhi: Bill & Melinda Gates Foundation; 2008. Available: http://www.gatesfoundation.org/avahan/Documents/Avahan_UseItOrLooseIt.pdf.

18. Barongo L, Borgdorff M, Mosha F, Nicoll A, Grosskurth H, Senkoro K. The epidemiology of HIV-1 infection in urban areas, roadwide settlements and rural villages in Mwanza Region, Tanzania. AIDS. 1992;6:1521–8. pmid:1492935

19. Mukhopadhyay C, Nath G, Gulati A, Mohapatra S. Prevalence of HIV among low and high risk population of eastern part of northern India. J Commun Dis. 2001;33:136–42. pmid:12170934

20. Mehta S, Gupta A, Sahay S, Godbole S, Joshi S, Reynolds S. High HIV prevalence among a high-risk subgroup of women attending sexually transmitted infection clinics in Pune, India. J Acquir Immune Defic Syndr. 2006;41:75–80. pmid:16340477

21. Pettifor AE, van der Straten A, Dunbar MS, Shiboski SC, Padian NS. Early age of first sex: a risk factor for HIV infection among women in Zimbabwe. Aids. 2004;18(10):1435–42. pmid:15199320

22. Bhattacharya G. Sociocultural and Behavioral Contexts of Condom Use in Heterosexual Married Couples in India: Challenges to the HIV Prevention Program. Health Education & Behavior. 2004;31(1):101–17.

23. International Center for Research on Women: Too Young to Wed: The Lives, Rights and Health of Young Married Girls 2003 Available: http://www.icrw.org/files/publications/Too-Young-to-Wed-the-Lives-Rights-and-Health-of-Young-Married-Girls.pdf.

24. UNFPA, State of the World Population 2005. Available: http://www.unfpa.org/swp/2005/english/ch1/index.htm.

25. Early Marriage: A Harmful Traditional Practice UNICEF; 2005. Available: http://www.unicef.org/publications/files/Early_Marriage_12.lo.pdf.

26. Clark S. Early Marriage and HIV Risks in Sub-Saharan Africa. Studies in Family Planning. 2004;35(3):149–60. pmid:15511059

27. Jenson R, Thornton R. Early female marriage in the developing world, Gender and Development. JSTOR. 2003;11 (2):9–19.