Improvement of HAART in Brazil, 1998–2008: a nationwide assessment of survival times after AIDS diagnosis among men who have sex with men
BMC Public Health volume 15, Article number: 226 (2015)
In 1996, Brazil became the first developing country to provide free, universal access to HAART, laboratory monitoring, and clinical care to any eligible patient. As of June 2014, approximately 400,000 patients were under treatment, making it the most comprehensive HIV treatment initiative implemented thus far in a middle-income country, worldwide. The Brazilian epidemic is highly concentrated among men who have sex with men (MSM).
Four national information systems were combined and Cox regression was used to conduct retrospective cohort analysis of HAART availability/access on all-cause mortality among MSM diagnosed with AIDS reported to the information systems between 1998–2008, adjusting for demographic, clinical, and behavioral factors and controlling for spatially-correlated survival data by including a frailty effect. Multiple imputation by chained equations was used to handle missing data.
Among 50,683 patients, 10,326 died during the 10 year of period. All-cause mortality rates declined following introduction of HAART, and were higher among non-white patients and those starting HAART with higher viral load and lower CD4 counts. In multivariable analysis adjusted for race, age at AIDS diagnosis, and baseline CD4 cell count, MSM diagnosed in latter periods had almost a 50% reduction in the risk of death, compared to those diagnosed between 1998–2001 (2002–2005 adjHR: 0.54, 95% CI:0.51-0.57; 2006–2008 adjHR: 0.51, 95% CI:0.48-0.55). After controlling for spatially correlated survival data, mortality remained higher among those diagnosed in the earliest diagnostic cohort and lower among non-white patients and those starting HAART with higher viral load and lower CD4 lymphocyte counts.
Universal and free access to HAART has helped achieve impressive declines in AIDS mortality in Brazil. However, after a 10-years follow-up, differential AIDS-related mortality continue to exist. Efforts are needed to identify and eliminate these health disparities, therefore improving the Brazilian response towards HIV/AIDS epidemic.
By the end of 2013, 35 million [33.2 million-37.2 million] people were living with HIV worldwide and 1.5 million [1.4 million-1.7 million] AIDS-related deaths were reported in the latest Epidemiological Bulletin released by UNAIDS. The introduction of combination antiretroviral therapies has led to dramatic improvements in AIDS-related morbidity and mortality in countries where this treatment is available .
Highly active antiretroviral therapy (HAART) significantly improves the prognosis of HIV-infected persons by reducing HIV viral load, increasing CD4+ lymphocyte counts and delaying progression to AIDS, ultimately reducing mortality rates [2,3].
Pre-exposure prophylaxis (PrEP) for HIV infection has greatly improved the field of HIV prevention , and a recent systematic review suggests its efficacy among high-risk groups such as men who have sex with men – MSM .
Early detection and management of HIV infection are key challenges worldwide, but a growing body of evidence supports the immediate use of HAART to maintain CD4 count and functionality, limit the size of the HIV reservoir, and reduce the risk of onward viral transmission . Following such groundbreaking results, Brazil has recently announced plans to adopt the “Treatment as Prevention” strategy to curb HIV and AIDS .
However, AIDS-related deaths continue to occur, even in high-income countries, and the emergence of drug-resistant HIV variants and drug toxicities remain major barriers and challenges to successful long-term antiretroviral efficacy [8,9].
Although stable, on average, the epidemic in Central and South America is currently concentrated in specific groups, mainly among MSM [1,10]. In Brazil MSM accounted for more than 32% of the cumulative reported cases of AIDS from 1980 to 2010 among men older than 13 years of age . Data from a Brazilian meta-analysis indicated a pooled HIV prevalence of 13.6% (95% CI:8.2-20.2) among MSM . National data from sexually transmitted infections’ clinics in Brazil demonstrated a 1.7% overall prevalence for HIV among clinic attendees . However, when stratified for subgroups defined by sexual exposure, prevalence among MSM was found to be 5.7%, compared to 1.1% among heterosexual men. Data from three voluntary counseling and testing services reported 24.8% prevalence among MSM and 4.3% among heterosexual men . HIV incidence at VCT units in Rio de Janeiro, Brazil, indicated incidence up to 11 times higher among MSM when compared to heterosexual men [13,14].
Brazil was the first middle-income country to provide full access to HAART, laboratory monitoring and clinical care at no cost at the point of health care delivery to any eligible patient, since 1996 . As of June 2014, approximately 400,000 patients were receiving HAART in Brazil , making it the most comprehensive HIV treatment initiative implemented thus far in a middle-income country, worldwide .
Recent Brazilian observational studies have brought essential information about mortality patterns among HIV-positive patients under follow-up in Brazil’s referral centers [17,18]. However, clinical cohorts from selected referral centers do not assess real-life conditions found in busy, understaffed public health clinics scattered all over Brazil, a country where social and regional heterogeneities are very relevant in the field of AIDS care [19,20], as well as health care in a broad sense .
Herein, we report differences in survival from AIDS diagnosis according to period of AIDS diagnosis, within the unique Brazilian setting. The study evaluated all MSM receiving treatment in the country between 1998 and 2008, therefore avoiding selection bias present in analyses from specific subsets (such as referral centers from major southeastern metropolitan areas, particularly Rio de Janeiro and Sao Paulo) of the country’s population living with HIV/AIDS.
This study utilized four databases comprising different longitudinal information of all people living with HIV/AIDS (PLWHA) under treatment and care through the Brazilian public health system. These databases comprise all patients receiving any antiretroviral medicine in Brazil since such medicines are not procured, purchased, or delivered by private health facilities or pharmacies for patients who are not hospitalized.
These databases contain the core information of Brazil’s surveillance system. Their basic characteristics are briefly summarized as follows:
Socio-demographic and basic clinical (for the sake of diagnosis) information on AIDS cases (SINAN-AIDS - Information System for Notifiable Diseases/AIDS);
Data from exams conducted in the network of accredited public laboratories, particularly TCD4+/TCD8 lymphocyte counting and HIV viral load (SISCEL - National Database for Laboratory Tests);
Information about monthly ARV refills and therapeutic regimen changes over time (SICLOM - Logistics Control System of ARV Medicines); and
The date and cause/s of death (SIM - National Mortality Database).
Whereas patients may seek private medical care and perform laboratory exams in private facilities, all PLWHA receive their antiretroviral drugs through SUS (Brazil’s Unified Public Health System), the single accredited source of ARVs in Brazil, and should be included in SICLOM database. At the time we conducted our analysis, SICLOM data have been updated to improve the comprehensiveness of collected data, which tend to be limited to the core information which is essential for medicines dispensing at the several health services belonging to the network of accredited public health units. Thus, to utilize data from all regions, we could not include additional information such as measures of adherence, side effects, drug regimen shifts, etc. as these data were not yet available for all regions and services.
Data from the entire population of PLWHA with AIDS diagnosis between 1998 and 2008, whose transmission category was homo/bisexual contact were extracted from those databanks and merged. The databank linkage procedures have been described in detail elsewhere .
For this analysis, participants were included if they were diagnosed as having AIDS between January 1, 1998 and December 31, 2008, and if their most probable route of HIV infection was male-male unprotected anal intercourse. Those patients were included because we aimed to analyze differences in survival from AIDS diagnosis among the MSM population, for which higher rates of initiation and better adherence to HAART during three different AIDS diagnosis periods: 1998–2001, 2002–2005, and 2006–2008 have been documented in the literature and by our own exploratory analyses.
The final merged database comprised the following covariates: date of birth, race, date of AIDS diagnosis, dates of initiation of HAART, CD4 counts and viral load at diagnosis, cause(s) and date of death. CD4 counts at initial diagnosis were classified into three categories: < 200 cells/mm3, 200–350 cells/mm3, and >350 cells/mm3. Viral load at initial diagnosis was classified as <4.00 log10 copies, 4.00-5.00 log10 copies, and >5.00 log10 copies.
The primary endpoint of interest in the present study was all-cause mortality during 1998–2008. Survival was calculated as the time elapsed from date of AIDS diagnosis until the date of death. Individuals who remained alive were censored at the end of data collection period (December 31, 2008). Time from AIDS diagnosis was transformed from days to years since diagnosis. Since analyses were based on secondary data, we did not have access to information on clinical visits; therefore we were not able to calculate losses to follow-up.
Kaplan-Meier curves were fitted to compare probabilities of survival according to selected period of diagnosis and CD4 count at diagnosis using the log-rank test. Cox proportional hazard regression modeling was used to examine covariates associated with survival among the following: age at AIDS diagnosis, race, period of AIDS diagnosis (1998–2001; 2002–2005; 2006–2008), first CD4+ lymphocyte count, HAART uptake, and first HIV-RNA viral load exam. The likelihood ratio test was used to select variables to be entered in the Cox models, adjusting for potential confounders. The assumption of proportional hazards was evaluated through the examination of Schoenfeld residuals .
Hazard rates are strongly influenced by selection effects operating in the population (e.g. higher availability of referral centers in a given locality), and individuals surviving up to a certain time will on average be less frail than the original population, therefore unobserved individual heterogeneity (or frailty) was also taken into account .
With 8.5 million square kilometers and over 202 million inhabitants according to the National Census, 2010, updated as of July 2013  and regular updates issued by the Brazilian Institute for Geography and Statistics (IBGE), Brazil is the largest and most populous country in Latin America, and faces deep health and socio-economic differences between and within each region. In an attempt to account for the heterogeneity between and within different Brazilian regions and localities, Cox regression analyses were repeated with the introduction of a random effect term to assess the potential association between survival times within a cluster – state or municipality of residency – using gamma-distributed frailty . A frailty indicator >1 denotes a faster event (death) rate, while groups will frailty <1 survival tend to be longer .
To estimate confidence intervals for Cox's proportional hazards models, robust methods were used. Maximum penalized likelihood estimation was used to fit frailty models. To handle missing data, we implemented multiple imputation by chained equations using classification and regression trees (MICE-CART) as the conditional models for imputation . This method can result in more plausible imputations, and hence more reliable inferences, in complex settings than the naive application of standard sequential regression imputation techniques.
Creating multiple imputations, as opposed to single imputations, accounts for the statistical uncertainty in the imputations. In addition, the chained equations approach is very flexible and can handle variables of varying types (e.g., continuous or binary), and has been reported as very reliable when used in different datasets [27,28]. All analyses were performed using R software version 3.0.2 (library survival) and Stata 13.0 (College Station, TX) .
No identifiable personal data were used for this study. The dataset used in the study is not openly available. Permission to use non-identifiable, individual data extracted from HIV/AIDS administrative databases was granted by the Brazilian Ministry of Health, Department of STDs, AIDS and Viral Hepatitis. Ethical approval for use of encrypted and aggregated data was also obtained from the Oswaldo Cruz Foundation Research Ethics Committee (Protocol CEP/ENSP 179/08) and the Brazilian National Research Ethics Committee (CONEP 15308).
Characteristics of the patients
Between January 1 1998 and December 31 2008, 50,683 MSM were diagnosed with AIDS in Brazil. Of these, 20,702 (40.9%) were diagnosed between 1998–2001, 19,280 (38.0%) in 2002–2005, and 10,701 (21.1%) in the last period under analysis, 2006–08. The overall mean age at AIDS diagnosis was 35.4 years (SD: 9.5). Data on race from 39.0% of our sample was missing. Although all patients had a confirmed AIDS (clinical) diagnosis, there were no information respecting 29.5% of the CD4 counts and 42.5% HIV-1 RNA viral load determinations in the merged database (Table 1). Due to the significant amount of missing data, MICE-CART multiple imputation was conducted before running survival analysis.
There were 10,326 deaths during the 10 year period. The median follow-up times for patients in three diagnosis periods are: period 1998–2001: 8.36 years; period 2002–2005: 4.68 years; and period 2006–2008: 1.44 years.
In Kaplan-Meier analyses stratified by diagnostic period, those diagnosed in the first period had significantly lower survival since AIDS diagnosis (p < 0.0001), with visible improvements among patients from the second and third cohort period (p < 0.0001). Taking into consideration the short follow-up time for participants in diagnostic period 2006–2008, there appeared to be little difference in survival between the last two diagnostic periods (Figure 1).
In the unadjusted analysis, the risk of death was lower for both periods 2002–2005 (HR: 0.74; 95% CI:0.71-0.77; p < 0.001) and 2006–2008 (HR: 0.81; 95% CI:0.76-0.86; p < 0.001), compared to the first period of diagnosis (1998–2001). The unadjusted HR for the period 2006–2008 was higher (HR: 1.09; 95% CI:1.02-1.17; p < 0.007) when compared to the reference of 2002–2005 to assess difference between latter periods (data not shown).
Predictors of survival
Controlling for other variables, the risk of death among latter cohorts of diagnosis was lower when compared to the cohort diagnosed between 1998–2001 (Period 2002–2005 adjHR: 0.54, 95CI:0.51-0.57; Period 2006–2008 adjHR: 0.51, 95%CI:0.48-0.55). There was no difference between the cohort diagnosed in 2006–2008 compared to the cohort diagnosed in 2002–2005, after adjustment in the Cox model (adjHR: 0.95; 95% CI:0.89-1.02). Additional predictors of longer survival included younger age at diagnosis and higher CD4 counts at diagnosis. Patients with an HIV-RNA viral load >5log10 copies/ml at diagnosis had a higher risk of death, as well as non-White patients (Table 2).
Upon multivariate Cox frailty regression that accounted for a random effect associated with the Brazilian state where patients lived (Table 2), the increased risk of death observed in patients diagnosed in the first period of time persisted. Both higher CD4 cells counts and lower HIV-1 RNA viral load remained as significant predictors of survival. Non-White patients remained with a higher risk of death; however, there was no longer a difference between white participants and those with unspecified race.
The estimated frailty variance was significantly different from zero, suggesting that the risk of death after AIDS diagnosis was heterogeneous among subjects living in different states (Table 2). We adjusted the final model with and without MICE-CART imputation procedure; however there were no relevant influence in the results (data not shown).
Greater frailty effects (and therefore greater mortality) were found in larger states located in the Tropical Rain Forest region and in states located in the Northeast, Brazil’s poorest region. On the other hand, industrialized areas from the Southern and Southeastern regions accounted for the smallest frailty effect (Figure 2).
This is the first Brazilian nationwide analysis of survival after AIDS diagnosis among MSM. Our study compared different time periods after the full introduction of HAART over the whole country, following the approval of the federal legislation regulating the right to the universal access to antiretroviral medicines for people living with AIDS, as of the end of 1996 . Improvements in survival were observed for latter periods of diagnosis, with almost 50% reduction in risk of death when compared to the earliest cohort. Protective effects were also observed with lower age and higher CD4 classifications at initial diagnosis.
Previous papers documented the dramatic improvement of survival before and after HAART , and identified lower mortality rates among PLWHA from Brazil, when compared with patients from other Latin American countries . Changes in the main causes of death from AIDS to non-AIDS related conditions have also been identified [18,32], however, those studies were conducted in a pool of referral centers or in a single referral center, and none of them assessed nationwide data, taking in consideration their pronounced heterogeneity.
The observed improvements in the latter periods of diagnosis compared to those diagnosed between 1998–2001 are favorable. No difference, however, has yet been detected between those diagnosed in 2006–2008 compared to those diagnosed in 2002–2005. This finding may be attributed to barriers and problems in the current treatment program and the concrete operations of the respective network of public health services. Among the most relevant of such challenges, some are highlighted as follows: i) the increasing pressures over costs , challenging the very sustainability of universal access over time; ii) the deep heterogeneity of the Brazilian network of health services, either in the specific field of AIDS management and care  or from the perspective of the health system in a broad sense ; and, iii) last but not the least, the late presentation of a substantial fraction of patients . Nonetheless, it is important to note that the median follow-up time for patients diagnosed in 2006–2008 is substantially shorter and future research may observe improvements after longer durations of follow-up of this cohort.
Other problems to be tackled by clinicians and health managers include the progressive, albeit moderate, increase of secondary resistance  and the permanent challenge posed by side effects of different antiretroviral medications, such as lipodystrophy and metabolic syndromes . However, despite all such caveats and difficulties, the present paper documents a progressive increase of survival times among gay men over successive periods of the decade 1998–2008. Similar findings were identified by a multi-site study comprising data from several North American cohorts. According to the study conducted by Samji and colleagues , who analyzed data from almost 23,000 treatment-naive ART patients from the MSM population, the life-expectancy at age 20 years increased with calendar time, from 53.3 in 2000–2002 to 69.3 in 2006–2007. Future Brazilian studies should include more recent timeframes, as well as additional information about HAART uptake and treatment adherence.
Our statistical modeling incorporated frailty effects and demonstrated a pronounced effect of geographic, social and/or health care quality upon survival times. It is not possible to disentangle the independent effects of each one of these factors; however, our study highlighted the negative impact on survival among patients living in less industrialized and most deprived regions of Brazil. Patients living in such areas are likely to suffer from an adverse and synergic effect of poverty and lack of social support, receiving less than optimal care in settings where health services are scarce, overburdened, less accessible and of lower quality .
As in any analysis of large secondary datasets, there is a trade-off between the augmented comprehensiveness of patient data under analysis and the deficiencies in terms of coverage and quality of databases. In this sense, the missing information identified in the current study constitutes serious limitations, not only in terms of the resulting analyses, but also from the broader perspective of monitoring and policymaking.
Such barriers and caveats must be progressively overcome by policymakers, managers, health professionals and the civil society in a continental-sized, deeply heterogeneous country, housing the most comprehensive network of ARV dispensing units worldwide. With increasing costs and the pressing need to switch therapeutic regimens in order to minimize both HIV resistance and different adverse effects of ARVs, substantial improvement of information systems is mandatory and urgently needed, as already highlighted by Grangeiro and colleagues .
The improvement in survival from the first period of HAART is undeniable and should be commended. However, results should be viewed in light of some study limitations. First, the process of merging four national databases resulted in missing data among key variables, such as CD4+ lymphocyte counts and viral load measures. Such limitations forced us to conduct multiple imputation (MICE-CART) of several variables before running the survival analysis. We did not conduct MICE-CART for race, a key variable with a significant proportion of missing information. This was a methodological and operational decision: patients with private health insurance usually have scarce information available on national datasets and our MICE-CART could not be accurately conducted with the available information. This problem calls for an urgent need to provide better training of professionals and clerks in charge of monitoring the Brazilian AIDS epidemic and to increase the accountability and consistency of such systems. A second limitation lies in the comparison of survival across time periods of initial diagnosis. Because our dataset is limited to 1998 to 2008, length of follow-up is shorter for those diagnosed in 2002–2005 and 2006–2008, limiting comparisons in survival after six to ten years following diagnosis. Nonetheless, these data are sufficient to allow for comparisons of trends in survival across different periods of diagnosis and CD4 count to AIDS diagnosis during the early years of HAART.
Brazil has been a worldwide leader in the global efforts to curb AIDS and to provide a standard of care and treatment for PLWHA at no cost at the point of delivery. Notwithstanding, the concerted effort from across different sectors of government, community leadership, and the civil society at large may be seriously compromised in the absence of a comprehensive and accurate monitoring system. Although our analyses documented some progress since the early years of HAART, current levels of missing information remain unacceptably high - a critical challenge as we enter in the third decade of the HIV/AIDS epidemic, with more complex and costly treatments, more patients under follow-up and more data being collected.
UNAIDS. The GAP Report. Geneva: Joint United Nations Programme on HIV/AIDS; 2014.
Van Sighem AI, Gras LA, Reiss P, Brinkman K, de Wolf F. ATHENA national observational cohort study. Life expectancy of recently diagnosed asymptomatic HIV-infected patients approaches that of uninfected individuals. AIDS. 2010;24(10):1527–35.
Kimmel AD, Charles M, Deschamps MM, Severe P, Edwards AM, Johnson WD, et al. Lives saved by expanding HIV treatment availability in resource-limited settings: the example of Haiti. J Acquir Immune Defic Syndr. 2013;63(2):e40–8.
Person AK, Hicks CB. Pre-exposure prophylaxis–one more tool for HIV prevention. Curr HIV Res. 2012;10(2):117–22.
Okwundu CI, Uthman OA, Okoromah CA. Antiretroviral pre-exposure prophylaxis (PrEP) for preventing HIV in high-risk individuals. Cochrane Database Syst Rev. 2012;7, CD007189.
Smith MK, Rutstein SE, Powers KA, Fidler S, Miller WC, Eron Jr JJ, et al. The detection and management of early HIV infection: a clinical and public health emergency. J Acquir Immune Defic Syndr. 2013;63 Suppl 2:S187–99.
Brazil to adopt B.C.’s Treatment as Prevention strategy as country’s national HIV/AIDS policy. http://cfenet.ubc.ca/news/releases/brazil-adopt-bc’s-treatment-prevention-strategy-country’s-national-hivaids-policy
Sullivan PS, Jones JS, Baral SD. The global north: HIV epidemiology in high-income countries. Curr Opin HIV AIDS. 2014;9(2):199–205.
Hosseinipour MC, Gupta RK, Van Zyl G, Eron JJ, Nachega JB. Emergence of HIV drug resistance during first- and second-line antiretroviral therapy in resource-limited settings. J Infect Dis. 2013;207 Suppl 2:S49–56.
Bastos FI, Cáceres C, Galvão J, Veras MA, Castilho EA. AIDS in Latin America: assessing the current status of the epidemic and the ongoing response. Int J Epidemiol. 2008;37(4):729–37.
Brazilian Ministry of Health. HIV/AIDS and Viral Hepatitis Department. HIV/AIDS and STI Epidemiological Bulletin. Available at http://www.aids.gov.br/publicacao/2013/boletim-epidemiologico-aids-e-dst-2013
Malta M, Magnanini MM, Mello MB, Pascom AR, Linhares Y, Bastos FI. HIV prevalence among female sex workers, drug users and men who have sex with men in Brazil: a systematic review and meta-analysis. BMC Public Health. 2010;10:317.
Velasco-de-Castro CA, Grinsztejn B, Veloso VG, Bastos FI, Pilotto JH, Fernandes N, et al. HIV-1 diversity and drug resistance mutations among people seeking HIV diagnosis in voluntary counseling and testing sites in Rio de Janeiro. Brazil PLoS One. 2014;9(1):e87622.
Barroso PF, Harrison LH, de Fatima Melo M, Batista SM, da Silva Bastos M, da Rosa Faulhaber JC, et al. Identification of a high-risk heterosexual cohort for HIV vaccine efficacy trials in Rio de Janeiro, Brazil, using a sensitive/less-sensitive assay: an update. J Acquir Immune Defic Syndr. 2004;36(3):880–1.
Malta M, Bastos FI, da Silva CM, Pereira GF, Lucena FF, Fonseca MG, et al. Differential survival benefit of universal HAART access in Brazil: a nation-wide comparison of injecting drug users versus men who have sex with men. J Acquir Immune Defic Syndr. 2009;52(5):629–35.
Greco DB, Simão M. Brazilian policy of universal access to AIDS treatment: sustainability challenges and perspectives. AIDS. 2007;21 Suppl 4:S37–45.
Oliveira MT, Latorre Mdo R, Greco DB. The impact of late diagnosis on the survival of patients following their first AIDS-related hospitalization in Belo Horizonte. Brazil AIDS Care. 2012;24(5):635–41.
Grinsztejn B, Luz PM, Pacheco AG, Santos DV, Velasque L, Moreira RI, et al. Changing mortality profile among HIV-infected patients in Rio de Janeiro, Brazil: shifting from AIDS to non-AIDS related conditions in the HAART era. PLoS One. 2013;8(4):e59768.
Grangeiro A, Escuder MM, Pereira JC. Late entry into HIV care: lessons from Brazil, 2003 to 2006. BMC Infect Dis. 2012;12:99.
Nemes MI, Alencar TM, Basso CR, Castanheira ER, Melchior R, Alves MT, et al. Assessment of outpatient services for AIDS patients, Brazil: comparative study 2001/2007. Rev Saude Publica. 2013;47(1):137–46.
Victora CG, Barreto ML, Do Carmo Leal M, Monteiro CA, Schmidt MI, Paim J, et al. Lancet brazil series working group. Health conditions and health-policy innovations in brazil: the way forward. Lancet. 2011;377(9782):2042–53.
Schoenfeld D. Partial residuals for the proportional hazards regression model. Biometrika. 1982;69:239–41.
Aalen OO. Effects of frailty in survival analysis. Stat Methods Med Res. 1994;3:227–43.
Brazilian Institute of Geography and Statistics - IBGE. National Population Estimates, July 2013. Available at http://www.ibge.gov.br/home/estatistica/populacao/estimativa2013/estimativa_dou.shtm
Dias SS, Andreozzi V, Martins MO, Torgal J. Predictors of mortality in HIV-associated hospitalizations in Portugal: a hierarchical survival model. BMC Health Serv Res. 2009;9:125.
White IR, Royston P, Wood AM. Multiple imputation using chained equations: Issues and guidance for practice. Stat Med. 2011;30(4):377–99.
Eisemann N, Waldmann A, Katalinic A. Imputation of missing values of tumour stage in population-based cancer registration. BMC Med Res Methodol. 2011;11:129.
R Core Development Team [computer software] version 2.13.0. Vienna: R Foundation for Statistical Computing. Available at http://www.r-project.org/
Samji H, Cescon A, Hogg RS, Modur SP, Althoff KN, Buchacz K, et al. Cohort collaboration on research and design (NA-ACCORD) of IeDEA. Closing the Gap: increases in life expectancy among treated HIV-positive individuals in the united states and Canada. PLoS One. 2013;8(12):e81355.
Marins JR, Jamal LF, Chen SY, Barros MB, Hudes ES, Barbosa AA, et al. Dramatic improvement in survival among adult Brazilian AIDS patients. AIDS. 2003;17(11):1675–82.
Tuboi SH, Schechter M, McGowan CC, Cesar C, Krolewiecki A, Cahn P, et al. Mortality during the first year of potent antiretroviral therapy in HIV-1-infected patients in 7 sites throughout Latin America and the Caribbean. J Acquir Immune Defic Syndr. 2009;51(5):615–23.
Pacheco AG, Tuboi SH, May SB, Moreira LF, Ramadas L, Nunes EP, et al. Temporal changes in causes of death among HIV-infected patients in the HAART era in Rio de Janeiro. Brazil J Acquir Immune Defic Syndr. 2009;51(5):624–30.
Nunn AS, Fonseca EM, Bastos FI, Gruskin S, Salomon JA. Evolution of antiretroviral drug costs in Brazil in the context of free and universal access to AIDS treatment. PLoS Med. 2007;4:e305.
Westin MR, Biscione FM, Fonseca M, Ordones M, Rodrigues M, Greco DB, et al. Resistance-associated mutation prevalence according to subtypes B and Non-B of HIV type 1 in antiretroviral-experienced patients in minas gerais. Brazil AIDS Res Hum Retroviruses. 2011;27(9):981–7.
Signorini DJ, Netto AM, Gabbay S, Monteiro MC, Signorini DH, Andrade MD, et al. A comparison of sonographic assessments and clinical questionnaire in the diagnosis of HIV-associated lipodystrophy. J Int Assoc Physicians AIDS Care (Chic). 2011;10(6):351–6.
Grangeiro A, Escuder MM, Menezes PR, Alencar R. Ayres De Castilho E. Late entry into HIV care: estimated impact on AIDS mortality rates in Brazil, 2003–2006. PLos One. 2011;6(1):14585.
The authors would like to thank the Brazilian National STD/AIDS and Viral Hepatitis Program for authorizing our access to their Nation-wide datasets
Dr. Strathdee acknowledges support from the Fogarty International Center, grant R25-TW007500.
Dr. Malta acknowledges support from INOVA/ENSP FIOCRUZ, University of California Center for AIDS Research (NIAID 5 P30 AI 036214), PAPES-VI/CNPq and CAPES (Proc. No BEX 2582/13-9).
Dr. Beyrer acknowledges support of the Johns Hopkins University Center for AIDS Research (P30AI094189).
The authors declare no conflict of interest. All analyses are based on public available databases, stripped of personal identification.
MM designed the study, organized the dataset, coordinated the statistical analysis and drafted the manuscript. CMFPS and MMM participated in the design of the study and performed the statistical analysis. AW contributed to statistical analysis and to the draft of the methods and results. AP reviewed the statistical analysis and helped to draft the manuscript. FIB, SAS and CB conceived of the study, and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Malta, M., da Silva, C.M.F.P., Magnanini, M.M. et al. Improvement of HAART in Brazil, 1998–2008: a nationwide assessment of survival times after AIDS diagnosis among men who have sex with men. BMC Public Health 15, 226 (2015). https://doi.org/10.1186/s12889-015-1530-y