Figures
Abstract
Background
Monitoring the HIV epidemic in a defined population is critical for planning treatment and preventive strategies. This is especially important in sub-Saharan Africa, which harbours the highest burden of the disease.
Objective
To estimate HIV incidence in adults aged 18-47 years old and to investigate spatial variations of HIV prevalence in Manhiça, a semi-rural area of southern Mozambique.
Methods
Two cross-sectional community-based surveys were conducted in 2010 and 2012 to determine HIV prevalence. Individual participants were randomly selected from the demographic surveillance system in place in the area and voluntary HIV counselling and testing was offered at the household level. HIV incidence was calculated using prevalence estimates from the two sero-surveys. Each participant’s household was geocoded using a global information system. The Spatial Scan Statistics programme was used to identify areas with disproportionate excess in HIV prevalence.
Results
A total of 1511 adults were tested. The estimated HIV prevalence in the community was 39.9% in 2010 and 39.7% in 2012. The overall HIV incidence was 3.6 new infections per 100 person-years at risk (PYAR) [95CI 1.56; 7.88], assuming stable epidemic conditions, and tended to be higher in women (4.9/100 PYAR [95CI 1.74; 11.85]) than in men (3.2/PYAR [95CI 1.36; 9.92]). One cluster with significant excess HIV prevalence was identified at the same geographic location in both surveys. This cluster had an HIV prevalence of 79.0% in 2010 and 52.3% in 2012.
Citation: González R, Augusto OJ, Munguambe K, Pierrat C, Pedro EN, Sacoor C, et al. (2015) HIV Incidence and Spatial Clustering in a Rural Area of Southern Mozambique. PLoS ONE 10(7): e0132053. https://doi.org/10.1371/journal.pone.0132053
Editor: Jean K. Carr, St. James School of Medicine, ANGUILLA
Received: January 7, 2015; Accepted: June 9, 2015; Published: July 6, 2015
Copyright: © 2015 González 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
Data Availability: Due to ethical restrictions related to participant consent, all relevant data are available upon request to Denise Naniche (dsuzanne@clinic.ub.es).
Funding: This work was supported by the European and Developing Countries Clinical Trials Partnership (EDCTP; http://www.edctp.org/) as part of the AfrEVacc consortium (grant: CT.2006.33111.001) and co-funded by the Fondo de Investigación Sanitaria from the Spanish Ministry of Health. The CISM receives core funding from the Spanish Agency for International Cooperation (AECI) and the HIV VCT units were in part supported by the Agència de Cooperació Catalana. RG was partially supported by an EDCTP grant (IP.07.31080.002).
Competing interests: The authors have declared that no competing interests exist.
Introduction
Surveillance of the HIV epidemic is essential in sub-Saharan African (SSA) countries where the highest burden of the disease concentrates, accounting for 69% of the people living with HIV worldwide[1]. HIV prevalence estimates have been based on data from women attending the antenatal clinics and more recently on population-based surveys[2]. However, for monitoring the HIV epidemic prevalence measures may not be as useful as it is measuring the appearance of new infections. The accepted gold standard method for measuring population-level HIV incidence is a prospective cohort study that assess the occurrence of new infections in a well-defined HIV-negative population followed over time and tested at regular intervals for HIV infection[3, 4]. However, this method is logistically complex and costly. Over the last years, methods for estimating HIV incidence using data from single-round cross-sectional surveys have been developed and validated in several countries[3–7].
While the HIV epidemic has been extensively studied through prevalence and incidence estimates, geographic patterns of HIV distribution at local community level have been underexplored in SSA. The analysis of HIV epidemics at a micro geographic scale including identification of hot spots is necessary to understand the driving forces of HIV transmission patterns [8]. Spatial analysis at the local level may also enable policy makers to design effective and culturally acceptable preventive measures or outreach programs for increasing antiretroviral therapy (ART) uptake in clusters of high HIV prevalence.
Mozambique is one of the ten countries with the highest HIV prevalence in the world with 1.4 million (95% Confidence Interval (95CI) 1.2–1.5) people living with HIV according to UNAIDS estimates [9, 10]. The first population-based survey conducted in 2009 showed that national HIV prevalence was 15% (95CI 13.9–16) in individuals aged 15–49 years old reaching over 19% in the southern region of the country[11]. Recent community-based estimates from central and southern Mozambique indicate that HIV prevalence may be up to 32.6% and 39.9%, respectively in some populations[12, 13]. No data on HIV incidence at community-level are currently available.
The present study reports the results of two community based cross-sectionals surveys to estimate HIV prevalence and incidence in adults aged 18–47 years in a semi-rural district of southern Mozambique. In addition, a geographical analysis was performed in order to investigate spatial variations and clusters of HIV prevalence in the area under study.
Methods
Ethics statement
The study protocol and informed consent form were reviewed and approved by the National Committee on Health Bioethics of Mozambique and the Hospital Clínic of Barcelona Ethics Committee (Spain). Written informed consent was obtained from all study participants.
Study area and population
The study was carried out in Manhiça District, a semi-rural area in Maputo Province, in southern Mozambique. The characteristics of the study population and site have been described elsewhere [12, 14]. Briefly, since 1996 the Centro de Investigação em Saúde de Manhiça (CISM) runs a continuous demographic surveillance system (DSS) for vital events and migrations. In 2007, there were 160 000 inhabitants in the district and currently the DSS covers over 92 000 inhabitants [15]. The preferred destinations for emigrants from Manhiça are the Maputo city followed by South Africa and other districts within Maputo province[14]. Since 2003 the CISM collaborates with the Mozambican HIV/AIDS control program through the establishment and continuous support of voluntary counseling and testing centres at health facilities, provision of antiretroviral drugs and diagnosis tests, and contributing to the clinical management of patients among other activities. Estimates from 2010 showed that the HIV prevalence at the antenatal clinics of the Manhiça District Hospital was 29.4%, and of 39.9% at the community level [12].
Study design
Two cross sectional community-based studies to determine age and sex-specific HIV prevalence in individuals aged 18–47 years old were conducted in 2010 [12] and 2012. The sample size was calculated with a precision of 0.05 assuming a 20% HIV prevalence and a 95% confidence interval and it was deemed that 232 subjects per age group would be needed to determine age and sex specific prevalence. The width of the defined age groups was designed to be equal among the three age groups (18–27, 28–37 and 38–47 years of age). A secondary objective was to calculate HIV incidence in the community from cross-sectional prevalence data at two pre-defined time points [3, 5–7]. In order to ensure HIV incidence calculations for all age groups, the second cross sectional survey included an upper age group of 47–50 since the model takes into account the time interval between surveys[6].
Study procedures
Selection of participants has been described in detail elsewhere [12]. Briefly, random lists of adults living in the study area stratified by age group and sex were generated and organized by neighborhood. The study inclusion criteria were: age 18–47 or 18–50 years in the first and second surveys respectively, being resident of the main study area, and willing to participate in the study after signing a written informed consent. The individuals were visited at home by a study field worker who explained the objectives of the study. If the candidate agreed, another home visit was made to provide more information about the study and HIV testing. Voluntary HIV counselling and testing was offered at the household level but only the randomly selected individuals contributed data to prevalence estimates. Study procedures of the household visits have been described in detail elsewhere[12]. Recruitment was stopped once the minimum sample size for each age and sex group was reached. Basic sociodemographic information was recorded onto a specifically designed case report form using mobile devices. Rapid HIV testing was performed by fingerprick following national recommendations using two rapid tests; the Determine HIV 1/2 test (Abbott laboratories, North Chicago, IL, USA), and the UniGold HIV test for confirmation of an HIV-positive result (Trinity Biotech, Bray, Ireland). HIV positivity was defined only if both tests were positive. Participants with a positive result were offered medical follow up at the Manhiça outpatient clinic, which included CD4 counts, clinical management and provision of ART according to national guidelines [16].
Data management and statistical analysis
Data from the mobile devices was transferred to the CISM centralized database. The Open Data Kit software (http://opendatakit.org/) was used for data management. Information on participant’s religion and migration history was obtained from the CISM DSS [17]. Individuals under the DSS are considered immigrants if they were not born in the Manhiça DSS area and had lived there for at least three months [14]. One-way and two-way contingency tables were generated for description of the categorical variables and calculation of proportions and p-values, taking into account the probability of sampling and stratification to extrapolate the data from the survey to the community[18]. The statistical analysis was performed using Stata statistical software version 12 (Stata Corp., College Station, Texas, USA).
HIV incidence estimation
HIV incidence was derived from the HIV prevalence estimates at the two cross-sectional surveys from 2010 and 2012 according to the model proposed by Hallet et al.[6, 19, 20] and validated in other settings [4, 5]. The method has been described in detail elsewhere [5, 6]. Briefly, it is based on the synthetic cohort principle and relies on the decomposition of prevalence changes by age group of width r (usually 5 years) between two cross-sectional surveys separated by T years of time. It assumes that individuals of age a years in the first survey will be represented by individuals aged a + T years in the second survey. Thus, the HIV prevalence in the second of two surveys represents the sum of the incident HIV infections and the surviving HIV-infected individuals from the first survey, even if the survey does not include the same individuals. The model adjusts for mortality rates using assumptions for three HIV epidemic situations, early, stable and declining [5, 6, 21]. Confidence intervals of incidence estimates were generated by bootstrap [19]. HIV incidence was estimated by age groups of 5 years width.
Spatial clustering detection
For the spatial analysis, “clusters of HIV infection” were defined as “geographical areas with a disproportionate excess in HIV prevalence compared to the surrounding areas”[22]. In a specific scale, the importance of HIV infection of an area defined as a cluster differs from that of a neighborhood area. Spatial Scan Statistics (SaTScan) programme (http://www.satscan.org/) and Kulldorf’s methodology were used to identify clusters of HIV infection [22–24]. For detecting clusters, SaTScan gradually scans a window across space and compares the number of observed and expected cases inside the window at each location. The radius of the window varies continuously in size from zero to an upper limit of 50% of the entire study area. This method creates an infinite number of circular windows [24]. The window with the maximum likelihood ratio is identified as the “most likely cluster”. For each window, the likelihood function is computed as follows:
Where C is the total number of cases, c is the observed number of cases within the window and E[c] is the covariate adjusted expected number of cases within the window under the null hypothesis. SaTScan takes into account and adjusts for population density in the analysis. Finally, Monte Carlo hypothesis testing was performed to calculate p-values for detected clusters. The number of DSS individuals that were eligible for participation at the first and second survey (background population at risk) was taken into account for cluster identification (N = 13618 and N = 20309, respectively).
Results
HIV prevalence in 2010 and 2012 surveys
Detailed results of the 2010 sero-survey are described in Gonzalez et al.[12]. In 2012, 888 adults were approached and given an appointment card for a later mobile team visit, following the same methodology as the previous survey [12]. Of the 888 adults invited, 789 accepted to participate and were recruited with an acceptance rate of 88.9%, similar to that of 2010 of 86.1%[12]. There were no significant differences in sex, age or education of participants between surveys (Table 1). The community HIV prevalence in 2012 was 39.7% (95CI 36.0;43.5), similar to the 39.9% (95CI 35.9;43.8) observed in 2010[12]. The age and sex-specific HIV prevalence were also similar in both surveys (Fig 1). HIV prevalence was significantly lower in the 18–27 year age group for both men and women as compared to the older age groups in both survey years (p<0.001).
HIV incidence estimates
Indirect HIV incidence estimates derived from HIV prevalence approximated 2.9 new infections/100 person-year (PYAR) (95CI 1.40; 7.14) for early epidemic conditions, 3.6 infections/100 PYAR (95CI 1.56; 7.88) for a stable HIV epidemic, and 4.7 infections/100 PYAR (95CI 2.37; 8.99) for declining epidemic conditions. HIV incidence point estimates in the three epidemic conditions tended to be higher in women (stable epidemic: 4.9/100 PYAR [95CI 1.74; 11.85]) than in men (stable epidemic: 3.2/PYAR [95CI 1.36; 9.92]; Table 2). Age stratification by 5 years age-group revealed a peak in HIV incidence trends at 35–39 years of age which was not statistically significant (S1 Fig).
Spatial analysis of HIV prevalence
The spatial scan analysis identified a small cluster of significant excess of HIV prevalence in 2010 (<0.05) with 19 individuals out of the 722 included in the analysis, and an overall area of 0.2 kilometres2 (km2) (Fig 2). HIV prevalence within this cluster was 79.0% (versus 36.5% outside). When using data from the 2012 sero-survey, the spatial scan identified a cluster of higher HIV prevalence centred in the same location as in 2010 but spanning a larger geographical area with 109 individuals out of the 789 included in the analysis (Fig 2). The 2012 cluster was approximately 1.4 km from the Maragra sugar mill, had an overall perimeter of 7.6 km and spanned 2.7 km2. The HIV prevalence inside the 2012 cluster was 52.3% (95CI 42.9; 61.7) as compared to 34.7% (95CI 31.1; 38.3) outside (p = 0.002; OR = 2.05 [95CI 1.36; 3.08]; p = 0.001). S2 Fig shows a map of the Manhiça DSS population density in 2012 and location of the cluster, which is independent of population density. There were no significant differences in sex, age, marital status education orliteracy between individuals from the high HIV prevalence cluster and those from the lower prevalence surrounding area (Table 3). There was however a significant difference in the proportion of immigrants within the cluster (56.9%) compared to that outside (39.6%; p = 0.001).
Discussion
This study has shown that in the DSS area of Manhiça, in southern Mozambique, the average adult community HIV prevalence in 2012 was very high, nearly 40%, with little change over that found two years earlier. Prevalence was highest in the older population as compared to the 18–27 age group. The estimated number of new infections was 3.6/ 100PYAR, with a slightly higher rate in women at 4.9 /100 PYAR. The findings also show that there is spatial heterogeneity in the observed HIV prevalence within the Manhiça DSS, with a small cluster (of 0.2 km2 area) of significant excess of HIV prevalence identified in the 2010 survey which, increased to a 2.7 km2 area in the 2012 survey.
The consistent results showing high prevalence of HIV both in 2010 and 2012 in Manhiça confirms the magnitude of the epidemic in this southern region of the country. A proportion of the high prevalence is likely to be due to increased survival in HIV-positive individuals stemming from the scale up of ART coverage. High HIV prevalence has been described in young women from the centre of Mozambique and in neighbouring countries such as South Africa and it has been associated with both continued new infections, as well as increased survival of the HIV-positive population receiving ART [13, 25–28]. The impact of ARTs rollout on contributing to increased HIV prevalence has recently been shown in neighbouring KwaZulu Natal, South Africa. This was found to be the greatest at the initial phases of the ART scale up, from zero to 40% coverage (between 2004 and 2010), stabilizing afterwards despite continued increases in ART coverage [26, 29]. Mozambique, like KwaZulu Natal, began ART scale-up in 2003 and recent reports estimate that ART coverage of adults in need of treatment in Mozambique has hovered at an estimated 50% since 2010 [30, 31]. This level of ART coverage is indeed associated with decreased risk of HIV acquisition [32], however there still remains a substantial pool of HIV-infected untreated individuals potentially contributing to the propagation of new infections, with the incidence of new infections suggested to be higher in the older population [29, 33].
In this study, an overall crude incidence of 3.6 new HIV infections/100 PYAR was estimated from prevalence sero-surveys in Manhiça, assuming stable epidemic conditions. Sex-specific incidence was estimated at 4.9 infections/100 PYAR for women and 3.2 infections/100 PYAR for men Previous estimates of HIV incidence in Manhiça from prevalence data in women of reproductive age from 1999 to 2008 ranged from 11.6 to 12.2 infections/100 PYAR[19]. Although study populations are not comparable, the current data points to a possible decrease in incidence of new infections in this area of the country. In neighbouring KwaZulu Natal estimated HIV incidence has been reported from longitudinal cohorts with a range of 6.3 to 14.8 new infections per 100 PYAR in 2009 and timid decreases in new infections since 2011 [25, 34]. Two recent cohort studies conducted among women at higher risk for HIV acquisition in neighbouring Mozambican areas report similar incidences to that found in Manhiça, of 4.6 infections/100 PYAR (95CI 2.7; 7.3) and 6.5 infections/ 100 PYAR (95CI 4.1; 9.9)[28, 35]. It has been described that HIV incidence may be lower in younger population[29]. However, the sample size in this study did not allow for precise estimations of incidence by age groups. The age-specific incidence estimates from the Manhiça population suggests trends of higher incidence at 35–40 years of age, in agreement with analyses using household-based prevalence data from Zambia and Niger that showed a peak of HIV incidence in men of the same age group[5]. In Manhiça, this peak was observed for the estimated incidence and it was not seen for age-specific HIV prevalence. It has been suggested that HIV incidence at older ages is due to widowhood exposing individuals to the risk of infection as they form new partnerships[36]. Apart from that, it has been described that peak HIV incidence in men may occur 5–7 years later than in women in a mature epidemic [5, 27].
In addition to population-based prevalence and incidence measures, the estimation of geographic distribution of HIV infection is essential in developing strategies for tackling the epidemic. These results suggest a significant geospatial cluster of excess HIV prevalence in a 7.6 km perimeter area within the Manhiça DSS, which is near to the Maragra sugar mill. Indeed, in 2012 HIV prevalence was 52.3% within what it can be considered a “hot spot”, as compared to 34.7% outside (p = 0.002). This cluster near the sugar mill was not associated with demographic risk factors for HIV except a higher rate of immigration. Southern Mozambique, similar to KwaZulu Natal in South Africa, is a stronghold of sugar production. Mozambique has four sugar mills, two of which, including Maragra, are in the Maputo province. Although sugar mills employ permanent workers, much of the workforce in the sugar industry rests on hired seasonal cane cutters. In Maragra, of a total workforce of approximately 5500 in 2012, close to 3000 were seasonal workers [37]. Studies in KwaZulu Natal and Malawi have described that migrants and their partners are at greater risk for HIV and other STIs [38]. This may be associated with increased frequency of sexual relationships outside of the couple and the effect of concurrent relationships with split households [39, 40]. The mobility of sugar mill workers is likely to be associated with a higher cluster of HIV prevalence. The occurrence of micro geographic clustering of HIV prevalence was recently described in a study across SSA using national demographic data[8]. In Mozambique, several high HIV prevalence clusters were found in Maputo, Gaza, Sofala and Zambezia provinces using data from 2009 [8]. Indeed, comparable to what observed in the Manhiça DSS area, micro clusters of high prevalence have also been observed at the level of neighbourhoods in a Tanzanian study conducted with HIV AIDS indicator survey data from 2003–2004[41]. It is challenging to assess the contribution of different risk factors to micro geographic clusters in HIV due to the sheer number of potential risk factors and the small population sizes in clusters. Small numbers also preclude accurate estimates of incidence within the clusters. Cuadros et al. have proposed that clustering may reflect differences in specific behavioural and biological factors which become amplified as a consequence of the epidemic being close to its epidemic threshold in SSA, thus leading to pockets of higher or lower HIV prevalence [8].
A possible limitation of the study could be that the HIV-specific mortality rates used in the estimates were derived from the model described by Hallett et al [5] and not specific to Manhiça. As an example, if the mortality rate used in our methodology was lower than that of Manhiça, HIV incidence would have been overestimated. However, small changes in mortality rates in the methodology do not significantly affect the incidence estimates. Furthermore, this methodology has been validated and employed successfully in many other SSA epidemiological settings [4, 19].
Most HIV prevention programs are guided by HIV prevalence data from administrative districts and HIV incidence data from specific risk groups. Results from Manhiça and other areas in SSA suggest the importance of taking into account “place” in terms of local geospatial differences in HIV prevalence in order to better understand the epidemic[38, 41, 42]. Geospatial prevalence data cannot be used alone since it does not estimate the rate of new infections. It is thus important to combine incidence estimates of new infections with geospatial estimates of prevalence clusters in order to best tailor local interventions and outreach programs for increasing ART uptake and retention in the complex SSA setting.
In conclusion, this study has shown that in this area of southern Mozambique the elevated HIV adult community prevalence remains highest in the older population and is accompanied by a continued high number of new infections especially among women. In addition, a geographic “hot spot” has been identified possibly associated with immigration to the area. These findings reinforce the need for continuous monitoring of the HIV epidemic and for further research on the geographical patterns of the infection, in order to guide and consolidate effective prevention strategies.
Supporting Information
S1 Fig. HIV incidence by age group and HIV epidemic condition.
https://doi.org/10.1371/journal.pone.0132053.s001
(TIF)
S2 Fig. Population density of the study area in 2010 and 2012 with the HIV clusters designated.
https://doi.org/10.1371/journal.pone.0132053.s002
(TIF)
S1 Table. Mortality rates and HIV prevalence estimates by age group and epidemic condition.
https://doi.org/10.1371/journal.pone.0132053.s003
(DOCX)
Acknowledgments
We are grateful to the study participants, field workers, HIV counsellors and all the staff from the Demography and Social Sciences Departments at the CISM, especially to Talita Matavel, Helena Boene, Rui Anselmo and Carolina Mindu.
Author Contributions
Conceived and designed the experiments: RG CP JJA CM DN. Performed the experiments: RG OJA CS ENP DN. Analyzed the data: RG OJA CP EdL. Contributed reagents/materials/analysis tools: KM CS JJA EM PLA CM. Wrote the paper: RG OJA KM CP ENP CS EdL JJA EM PLA CM DN.
References
- 1.
UNAIDS/WHO. Global report: UNAIDS report on the global AIDS Epidemic 2012. WHO Library Cataloguing-in-Publication Data 2012;UNAIDS(JC2417E).
- 2.
WHO, UNAIDS. Reconciling Antenatal Clinic-based surveillance and Population-based survey estimates of HIV prevalence in Sub-Saharan Africa. WHO. 2003.
- 3. Gregson S, Machekano R, Donnelly CA, Mbizvo MT, Anderson RM, Katzenstein DA. Estimating HIV incidence from age-specific prevalence data: comparison with concurrent cohort estimates in a study of male factory workers, Harare, Zimbabwe. Aids. 1998;12(15):2049–58. pmid:9814874.
- 4. Kim AA, Hallett T, Stover J, Gouws E, Musinguzi J, Mureithi PK, et al. Estimating HIV incidence among adults in Kenya and Uganda: a systematic comparison of multiple methods. PLoS One. 2011;6(3):e17535. Epub 2011/03/17. pmid:21408182; PubMed Central PMCID: PMCPMC3049787.
- 5. Hallett TB, Stover J, Mishra V, Ghys PD, Gregson S, Boerma T. Estimates of HIV incidence from household-based prevalence surveys. AIDS. 2010;24(1):147–52. Epub 2009/11/17. pmid:19915447; PubMed Central PMCID: PMCPMC3535828.
- 6. Hallett TB, Zaba B, Todd J, Lopman B, Mwita W, Biraro S, et al. Estimating incidence from prevalence in generalised HIV epidemics: methods and validation. PLoS Med. 2008;5(4):e80. pmid:18590346.
- 7. Wong KF, Tsai WY, Kuhn L. Estimating HIV hazard rates from cross-sectional HIV prevalence data. Stat Med. 2006;25(14):2441–9. pmid:16220498.
- 8. Cuadros DF, Awad SF, Abu-Raddad LJ. Mapping HIV clustering: a strategy for identifying populations at high risk of HIV infection in sub-Saharan Africa. International journal of health geographics. 2013;12:28. Epub 2013/05/23. pmid:23692994; PubMed Central PMCID: PMCPMC3669110.
- 9.
UNAIDS/WHO. Global report: UNAIDS report on the global AIDS Epidemic 2010. Joint United Nations Programme on HIV/AIDS 2010;UNAIDS/10.11E (JC1958E).
- 10. Merson MH, O'Malley J, Serwadda D, Apisuk C. The history and challenge of HIV prevention. Lancet. 2008;372(9637):475–88. pmid:18687461.
- 11.
Instituto-Nacional-de-Saúde., Instituto-Nacional-de-Estatística., ICF-Macro. Inquérito Nacional de Prevalência, Riscos Comportamentais e Informaçao sobre o HIV e SIDA em Moçambique 2009. Calverton, Maryland, EUA:INS,INE e ICF Macro. 2010.
- 12.
Gonzalez R, Munguambe K, Aponte J, Bavo C, Nhalungo D, Macete E, et al. High HIV prevalence in a southern semi-rural area of Mozambique: a community-based survey. HIV Med. 2012. pmid:22500780.
- 13. Zango A, Dube K, Kelbert S, Meque I, Cumbe F, Chen PL, et al. Determinants of prevalent HIV infection and late HIV diagnosis among young women with two or more sexual partners in Beira, Mozambique. PLoS One. 2013;8(5):e63427. Epub 2013/05/22. pmid:23691046; PubMed Central PMCID: PMCPMC3656941.
- 14. Nhacolo AQ, Nhalungo DA, Sacoor CN, Aponte JJ, Thompson R, Alonso P. Levels and trends of demographic indices in southern rural Mozambique: evidence from demographic surveillance in Manhica district. BMC Public Health. 2006;6:291. pmid:17137494.
- 15.
Instituto-Nacional-de-Estatística. Tercer Recenseamento Geral da População e Habitação 2007. 2007.
- 16.
Ministério da Saúde-República de Moçambique. Introduçao de Novas Normas para o seguimiento do paciente HIV positivo. Direcçao Nacional de Assistencia Médica Ministério da Saúde de Moçambique. 2013.
- 17. Sacoor C, Nhacolo A, Nhalungo D, Aponte JJ, Bassat Q, Augusto O, et al. Profile: Manhica Health Research Centre (Manhica HDSS). Int J Epidemiol. 2013;42(5):1309–18. Epub 2013/10/26. pmid:24159076.
- 18.
Stata-Corp-LP. STATA Survey Data Reference Manual. 2009.
- 19.
Perez-Hoyos S, Naniche D, Macete E, Aponte J, Sacarlal J, Sigauque B, et al. Stabilization of HIV incidence in women of reproductive age in southern Mozambique. HIV Med. 2011; https://doi.org/10.1111/j.1468-1293.2010.00908.x
- 20. Hallett TB, Zaba B, Stover J, Brown T, Slaymaker E, Gregson S, et al. Embracing different approaches to estimating HIV incidence, prevalence and mortality. AIDS. 2014;28 Suppl 4:S523–32. Epub 2014/11/20. pmid:25406755.
- 21. Hallett T, Garnett G. Estimating incidence of HIV infection in Uganda. Jama. 2009;301(2):159; author reply 60–1. pmid:19141760.
- 22. Wand H, Ramjee G. Targeting the hotspots: investigating spatial and demographic variations in HIV infection in small communities in South Africa. J Int AIDS Soc. 2010;13:41. Epub 2010/10/29. pmid:20979651; PubMed Central PMCID: PMCPMC2984578.
- 23. Sartorius B, Kahn K, Vounatsou P, Collinson MA, Tollman SM. Space and time clustering of mortality in rural South Africa (Agincourt HDSS), 1992–2007. Global health action. 2010;3. Epub 2010/09/15. pmid:20838482; PubMed Central PMCID: PMCPMC2938122.
- 24. Kulldorff M, Nagarwalla N. Spatial disease clusters: detection and inference. Stat Med. 1995;14(8):799–810. Epub 1995/04/30. pmid:7644860.
- 25. Nel A, Mabude Z, Smit J, Kotze P, Arbuckle D, Wu J, et al. HIV incidence remains high in KwaZulu-Natal, South Africa: evidence from three districts. PLoS One. 2012;7(4):e35278. Epub 2012/04/27. pmid:22536364; PubMed Central PMCID: PMCPMC3334938.
- 26. Zaidi J, Grapsa E, Tanser F, Newell ML, Barnighausen T. Dramatic increase in HIV prevalence after scale-up of antiretroviral treatment. AIDS. 2013;27(14):2301–5. Epub 2013/05/15. pmid:23669155.
- 27. Mossong J, Grapsa E, Tanser F, Barnighausen T, Newell ML. Modelling HIV incidence and survival from age-specific seroprevalence after antiretroviral treatment scale-up in rural South Africa. AIDS. 2013;27(15):2471–9. Epub 2013/07/12. pmid:23842131; PubMed Central PMCID: PMCPMC3815011.
- 28. Feldblum PJ, Enosse S, Dube K, Arnaldo P, Muluana C, Banze R, et al. HIV Prevalence and Incidence in a Cohort of Women at Higher Risk for HIV Acquisition in Chokwe, Southern Mozambique. PLoS One. 2014;9(5):e97547. Epub 2014/05/21. pmid:24842811; PubMed Central PMCID: PMCPMC4026130.
- 29. Rehle TM, Hallett TB, Shisana O, Pillay-van Wyk V, Zuma K, Carrara H, et al. A decline in new HIV infections in South Africa: estimating HIV incidence from three national HIV surveys in 2002, 2005 and 2008. PLoS One. 2010;5(6):e11094. Epub 2010/06/19. pmid:20559425; PubMed Central PMCID: PMCPMC2885415.
- 30.
UNGASS. United Nations General Assambly Special Session on HIV and AIDS. Mozambique Progress Report 2008–2009. UNAIDS. 2010.
- 31.
National AIDS Council-Republic of Mozambique Global AIDS Response Progress Report. For the Period 2010–2011. 2012; Available: http://wwwunaidsorg/en/dataanalysis/knowyourresponse/countryprogressreports/2012countries/ce_MZ_Narrative_Report[1]pdf. Accessed February 2014.
- 32. Tanser F, Barnighausen T, Grapsa E, Zaidi J, Newell ML. High coverage of ART associated with decline in risk of HIV acquisition in rural KwaZulu-Natal, South Africa. Science (New York, NY). 2013;339(6122):966–71. Epub 2013/02/23. pmid:23430656.
- 33. Eaton JW, Rehle TM, Jooste S, Nkambule R, Kim AA, Mahy M, et al. Recent HIV prevalence trends among pregnant women and all women in sub-Saharan Africa: implications for HIV estimates. AIDS. 2014;28 Suppl 4:S507–14. Epub 2014/11/20. pmid:25406753.
- 34.
Abdool Karim Q, Kharsany AB, Frohlich JA, Werner L, Mashego M, Mlotshwa M, et al. Stabilizing HIV prevalence masks high HIV incidence rates amongst rural and urban women in KwaZulu-Natal, South Africa. Int J Epidemiol. 2010. pmid:21047913.
- 35. Dube K, Zango A, van de Wijgert J, Meque I, Ferro JJ, Cumbe F, et al. HIV incidence in a cohort of women at higher risk in Beira, Mozambique: prospective study 2009–2012. PLoS One. 2014;9(1):e84979. Epub 2014/01/30. pmid:24475035; PubMed Central PMCID: PMCPMC3903474.
- 36. Lopman BA, Nyamukapa C, Hallett TB, Mushati P, Spark-du Preez N, Kurwa F, et al. Role of widows in the heterosexual transmission of HIV in Manicaland, Zimbabwe, 1998–2003. Sex Transm Infect. 2009;85 Suppl 1:i41–8. Epub 2009/04/03. pmid:19307340; PubMed Central PMCID: PMCPMC2654148.
- 37.
Jeque J. Geraçao de emprego de larga escala na indústria açucareira. Conferência Nórdica-Moçambicana. 2012; Available: http://www.dneapmpd.gov.mz/index.php?option=com_docman&task=cat_view&gid=61&Itemid=54. Accessed February 2014.
- 38. Feldacker C, Emch M, Ennett S. The who and where of HIV in rural Malawi: Exploring the effects of person and place on individual HIV status. Health & place. 2010;16(5):996–1006. Epub 2010/07/06. pmid:20598623; PubMed Central PMCID: PMCPMC3454472.
- 39. Coffee M, Lurie MN, Garnett GP. Modelling the impact of migration on the HIV epidemic in South Africa. AIDS. 2007;21(3):343–50. Epub 2007/01/27. pmid:17255741.
- 40. Lurie MN, Williams BG, Zuma K, Mkaya-Mwamburi D, Garnett GP, Sweat MD, et al. Who infects whom? HIV-1 concordance and discordance among migrant and non-migrant couples in South Africa. AIDS. 2003;17(15):2245–52. Epub 2003/10/03. pmid:14523282.
- 41. Msisha WM, Kapiga SH, Earls FJ, Subramanian SV. Place matters: multilevel investigation of HIV distribution in Tanzania. AIDS. 2008;22(6):741–8. Epub 2008/03/22. pmid:18356604; PubMed Central PMCID: PMCPMC2789284.
- 42. Feldacker C, Ennett ST, Speizer I. It's not just who you are but where you live: an exploration of community influences on individual HIV status in rural Malawi. Social science & medicine (1982). 2011;72(5):717–25. Epub 2011/02/15. pmid:21316134.