Risco espacial de óbito de pacientes com aids em Campinas, São Paulo, Brasil

Thiago Santos Mota Maria Rita Donalisio Liciana Vaz de Arruda Silveira Sobre os autores

RESUMO:

Objetivo:

O estudo teve como objetivo analisar o risco espacial de mortalidade por aids em setores censitários de município brasileiro de grande porte.

Método:

Foram estudadas três coortes retrospectivas de indivíduos notificados por aids e residentes no município de Campinas, São Paulo, nos períodos de 1980 a 1990 (coorte 1), 1996 a 2000 (coorte 2) e 2001 a 2005 (coorte 3), a fim de identificar a ocorrência do óbito. Foi ajustado um modelo semiparamétrico Bayesiano, empregando o método de Aproximação de Laplace Aninhada e Integrada (INLA), que permitiu obter mapas de risco de óbito nas três coortes estudadas.

Resultados:

A incidência de óbito nas coortes 1, 2 e 3 foram, respectivamente, 72,73; 32,21 e 13,11%. Os mapas de risco de óbito apontaram tendência de queda da mortalidade e evidenciaram os setores que apresentaram maior e menor risco em cada período. O modelo mostrou que, nas três coortes, os fatores associados ao maior risco de óbito foram: ser do sexo masculino e possuir idade de diagnóstico superior a 49 anos. A orientação homossexual/bissexual foi associada ao menor risco de morrer.

Conclusões:

Mesmo diante da diminuição do risco de morte por aids nos últimos anos, após o acesso à terapia antirretroviral altamente ativa, a maior incidência se concentra entre os mais pobres nas três coortes estudadas. Os mapas e o fatores de risco obtidos orientam possíveis ações e monitoramento da doença no município.

Palavras-chave:
Aids; Mortalidade; Terapia antirretroviral; Análise de sobrevida

INTRODUÇÃO

A epidemia da aids ao longo das últimas três décadas vem afetando vários segmentos da população, principalmente adultos jovens, com um grande impacto na economia e na qualidade de vida dos individuos11. Fonseca MGP, Bastos FI. Twenty-five years of the AIDS epidemic in Brazil: principal epidemiological findings, 1980-2005. Cad Saúde Pública. 2007. DOI: 10.1590/S0102-311X2007001500002
https://doi.org/10.1590/S0102-311X200700...
.

Vários estudos apontam uma redução drástica da mortalidade por aids em várias partes do mundo, em decorrência da utilização da terapia antirretroviral altamente ativa (HAART) 22. Fordyce EJ, Singh TP, Nash D, Gallagher B, Forlenza S. Survival rates in NYC in the era of combination ART. J Acquir Immune Defic Syndr. 2002; 30(1): 111-8.,33. Lima VD, Hogg RS, Harrigan PR, Moore D, Yip B, Wood E, et al. Continued improvement in survival among HIV-infected individuals with newer forms of highly active antiretroviral therapy. AIDS. 2007; 21(6): 685-92. DOI: 10.1097/QAD.0b013e32802ef30c
https://doi.org/10.1097/QAD.0b013e32802e...
,44. Schneider MF, Gange SJ, Williams CM, Anastos K, Greenblatt RM, Kingsley L, et al. Patterns of the hazard of death after AIDS through the evolution of antiretroviral therapy: 1984-2004. AIDS. 2005. DOI: 10.1097/01.aids.0000189864.90053.22
https://doi.org/10.1097/01.aids.00001898...
. No Brasil, o acesso à medicação antirretroviral combinada tem sido garantido pelo Ministério da Saúde desde o ano de 1996. O seu impacto no aumento da sobrevida dos pacientes, na redução de internações, de infecções oportunistas e no aumento de doenças crônicas tem sido registrado em vários estudos no país33. Lima VD, Hogg RS, Harrigan PR, Moore D, Yip B, Wood E, et al. Continued improvement in survival among HIV-infected individuals with newer forms of highly active antiretroviral therapy. AIDS. 2007; 21(6): 685-92. DOI: 10.1097/QAD.0b013e32802ef30c
https://doi.org/10.1097/QAD.0b013e32802e...
,55. Chequer P, Hearst N, Hudes ES, Castilho E, Rutherford G, Loures L, et al. Determinants of survival in adult Brazilian AIDS patients, 1982-1989. The Brazilian State AIDS Program Co-Ordinators. AIDS. 1992.,66. Guibu IA, Barros MB de A, Donalísio MR, Tayra Â, Alves MCGP. Survival of AIDS patients in the Southeast and South of Brazil: analysis of the 1998-1999 cohort. Cad Saúde Pública. 2011; 27(Suppl 1): S79-92. http://dx.doi.org/10.1590/S0102-311X2011001300009
http://dx.doi.org/10.1590/S0102-311X2011...
,77. Tancredi MV, Waldman EA. Survival of AIDS patients in Sao Paulo-Brazil in the pre- and post-HAART eras: a cohort study. BMC Infect Dis. 2014. DOI: 10.1186/s12879-014-0599-8
https://doi.org/10.1186/s12879-014-0599-...
.

Abordagens geoespaciais têm sido utilizadas para estudar tendências da epidemia de aids no Brasil, apontando disparidades regionais nos riscos de adoecer e morrer88. Teixeira TR de A, Gracie R, Malta MS, Bastos FI. Social geography of AIDS in Brazil: identifying patterns of regional inequalities. Cad Saúde Pública. 2014; 30(2): 259-71. DOI: 10.1590/0102-311X00051313
https://doi.org/10.1590/0102-311X0005131...
,99. Sousa AIA de, Pinto VL. Análise espacial e temporal dos casos de aids no Brasil em 1996-2011: áreas de risco aumentado ao longo do tempo. Epidemiol Serv Saúde. 2016; 25(3): 467-76. DOI: 10.5123/S1679-49742016000300003
https://doi.org/10.5123/S1679-4974201600...
.

No estado de São Paulo, na década de 1980, a doença acometia principalmente homossexuais, usuários de drogas injetáveis e pacientes transfundidos. Nas décadas seguintes observou-se que a epidemia atingiu vários grupos da sociedade, sofrendo um processo de heterossexualização, feminização e pauperização. Atualmente a doença tem se expandido em adultos jovens com baixo nível socioeconômico1010. Santos NJS, Tayra A, Silva SR, Buchalla CM, Laurenti R. A aids no Estado de São Paulo: as mudanças no perfil da epidemia e perspectivas da vigilância epidemiológica. Rev Bras Epidemiol. 2002; 5(3): 286-310. DOI: 10.1590/S1415-790X2002000300007
https://doi.org/10.1590/S1415-790X200200...
,1111. Pellini ACG, Bellini JH, Cavalin RF, Chiaravalloti Neto F, Zanetta DMT. Distribuição espacial das notificações de aids em mulheres no município de São Paulo nos períodos de 1999-2001 e 2009-2011: uma análise sob a ótica da desigualdade socioespacial. Rev Espaço Geogr [Internet]. 2015; 18(2). Disponível em: http://www.lsie.unb.br/espacoegeografia/index.php/espacoegeografia/article/view/441 (Acessado em 05 de janeiro de 2016).
http://www.lsie.unb.br/espacoegeografia/...
.

Nessa perspectiva, com intuito de entender a dinâmica da epidemia de aids no município de Campinas, São Paulo, Brasil, foram georreferenciados os casos notificados pelo Sistema de Informação de Agravos de Notificação (SINAN) e dos óbitos que constam no banco de dados do Sistema de Informação sobre Mortalidade (SIM) de moradores da cidade, com o objetivo de analisar o risco espacial de óbito entre pacientes notificados por aids nos setores censitários desse município, ao longo do tempo.

MÉTODO

O estudo envolveu três coortes retrospectivas de pacientes maiores de 13 anos, notificados por aids pelo SINAN, moradores do município de Campinas, a nordeste do estado de São Paulo (aproximadamente 1 milhão de habitantes). Foram excluídos 1.091 (26,13%) pacientes cuja notificação coincidiu com a data do óbito para reduzir o impacto dos tempos desses pacientes na sobrevida, centrando os resultados dos estimadores nos tempos dos pacientes com acesso ao seguimento e tratamento.66. Guibu IA, Barros MB de A, Donalísio MR, Tayra Â, Alves MCGP. Survival of AIDS patients in the Southeast and South of Brazil: analysis of the 1998-1999 cohort. Cad Saúde Pública. 2011; 27(Suppl 1): S79-92. http://dx.doi.org/10.1590/S0102-311X2011001300009
http://dx.doi.org/10.1590/S0102-311X2011...
Uma coorte foi composta por indivíduos notificados nos primeiros anos da epidemia para identificar a ocorrência de óbito, de 1980 a 1990 (coorte 1 - seguimento de 10 anos, com um tempo médio de sobrevida dos indivíduos de 5,7 anos, com um desvio padrão de 6,5 anos), outra com as notificações de 1996 a 2000 (coorte 2 - seguimento de 10 anos, com um tempo médio de sobrevida de 9,2 anos, com um desvio padrão de 5,6 anos) e uma terceira, no período de 2001 a 2005 (coorte 3 - seguimento de 7 a 9 anos, com um tempo de sobrevida médio de 8,9 anos, com um desvio padrão de 3 anos). Além do SINAN, foram obtidas informações no SIM sobre os óbitos dos pacientes notificados. Não foram incluídos pacientes que não constavam na lista de notificação do SINAN. Após checagem de duplicações e inconsistência dos bancos, as variáveis utilizadas no estudo foram: sexo, idade com que o indivíduo foi diagnosticado com o HIV, escolaridade, categoria de transmissão (homossexual/bissexual ou heterossexual, uso de drogas injetáveis), contagem de células linfócitos T CD4 categorizada em contagens menores ou iguais a 350 células por mm3 de sangue como critério de aids e analisada no modelo como variável contínua; além disso também foi obtida a localização das residências georreferenciadas nos setores censitários. Na coorte 1 havia um total de 286 indivíduos, sendo 245 (85,66%) homens e 41 (14,34%) mulheres. Na coorte 2, os dados correspondiam a 1.456 indivíduos, consistindo em 1.006 (69,09%) homens e 450 (30,91%) mulheres. Na coorte 3, observou-se 1.342 indivíduos, sendo 831 (61,92%) homens e 511 (38,08%) mulheres.

Para modelar o risco de óbito foi necessária uma abordagem que leva em conta as variáveis de interesse e a informação espacial dos casos. Utilizou-se o modelo de regressão semiparamétrico em um cenário Bayesiano, que permite modelar o risco de óbito em relação aos setores censitários.

O modelo semiparamétrico bayesiano utilizado está implementando no programa R, pacote INLA1212. Martino S, Akerkar R, Havard R. Approximate Bayesian Inference for Survival Models. Scand J Stat. 2011; 38(3): 514-28. DOI: 10.1111/j.1467-9469.2010.00715.x
https://doi.org/10.1111/j.1467-9469.2010...
,1313. R Development Core Team. R: A language and environment for statistical computing [Internet]. Viena: R Found Stat Comput Vienna; 2017. Disponível em: http://www.r-project.org/ (Acessado em 10 de março de 2017).
http://www.r-project.org/...
.

A suposição foi que os tempos até a ocorrência óbito sejam particionados em um intervalo finito, 0 = S 0 < S 1 < ⋅⋅⋅ < S K , t i < S K , para ∀t i com i = 1 ⋅⋅⋅, n, e que o risco basal seja constante no intervalo de tempo ]S k-1 , S k ], com k = 1 ⋅⋅⋅ , K, ou seja, h 0(t) = h k , para cada t ∈ ]S k-1 , S k ]. O modelo para o i-ésimo indivíduo no k-ésimo intervalo de tempo, pertencente ao j-ésimo setor censitário é dado por (Equação 1):

hikj(t) = h0(t) exp{Xi'β+f(zj)} = exp{Xi'β+f(zj)+bk}=exp{ηikj},(1)

Em que:

  • hikj(t):  função taxa de risco;
  • bk =log{h 0(t)} ;
  • β=(β1,…,βp )':  vetor de dimensão 1×p de parâmetros associados às covariáveis;
  • Xi:   vetor de dimensão p×1 de covariáveis observadas para o i-ésimo indivíduo; e
  • f(zj ):  efeito espacial no j-ésimo setor censitário da cidade de Campinas, São Paulo.

Para o componente f(z j ), adotou-se um modelo ICAR1414. Besag J, Kooperberg C. On conditional and intrinsic autoregressions. Biometrika. 1995; 82(4): 733-46. https://doi.org/10.1093/biomet/82.4.733
https://doi.org/10.1093/biomet/82.4.733...
, que pode ser expresso por (Equação 2):

fzj|fzw~N1ljj~wfZw, 1ljτz(2)

Em que:

  • lj:   número de setores que compartilham a mesma fronteira com o setor j; e
  • j~w:   dois setores censitários j e w são vizinhos;

Para o vetor de parâmetros são assumidas prioris gaussianas, β∼N(0 ;104I) para o vetor b=(b 1,⋅⋅⋅,b k ), prioris gaussianas, b k+1 -b kN(0,τb -1), de um passeio aleatório de primeira ordem (RW1) e para os hiperparâmetros, τb e τz , prioris independentes com distribuição gama, Γ(1;0,0001).

No cenário Bayesiano, para realizar as inferências desse modelo, a abordagem de estimação Bayesiana Aproximação de Laplace Aninhada e Integrada (INLA) é uma alternativa para o Markov Chain Monte Carlo (MCMC)1515. Rue H, Martino S, Chopin N. Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations. J R Stat Soc Ser B Stat Methodol. 2009; 71(2): 319-92..

Basicamente, a metodologia INLA foi criada para uma classes de modelos de regressão aditivos estruturados.

A ideia central da metodologia INLA é fazer integrações numéricas em um contexto Bayesiano, em que apenas aproximações determinísticas são usadas para obter as distribuições a posteriori, obtendo assim um ganho computacional expressivo1515. Rue H, Martino S, Chopin N. Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations. J R Stat Soc Ser B Stat Methodol. 2009; 71(2): 319-92..

Para agregar as localizações das residências dos indivíduos que pertencem ao mesmo setor censitário, utilizamos o software QGIS 2.8.2 e o mapa da cidade de Campinas, São Paulo, obtido no site do Instituto Brasileiro de Geografia e Estatística (IBGE), utilizando a Malha de Setor Censitário Urbano do ano 2000. Todas as análises estatísticas foram realizadas no programa R1313. R Development Core Team. R: A language and environment for statistical computing [Internet]. Viena: R Found Stat Comput Vienna; 2017. Disponível em: http://www.r-project.org/ (Acessado em 10 de março de 2017).
http://www.r-project.org/...
.

Para identificação de regiões da cidade de maior carência socioambiental, foi utilizado o mapa de distribuição espacial do Índice de Carência Socioambiental (ICSA), segundo classes socioeconômicas e ambientais por área de cobertura dos Centros de Saúde em Campinas, desenvolvido por Costa et al.1616. Costa JV, Donalisio MR, Silveira LV de A. Spatial distribution of dengue incidence and socio-environmental conditions in Campinas, São Paulo State, Brazil, 2007. Cad Saúde Pública. 2013; 29(8): 1522-32. DOI: 10.1590/0102-311X00110912
https://doi.org/10.1590/0102-311X0011091...
. Esse mapa possibilita a localização de três estratos de baixa, moderada e alta carência socioambiental a partir de análise estatística de componentes principais de variáveis demográficas e socioambientais do Censo 2000 da Fundação Instituto Brasileiro de Geografia e Estatística (FIBGE).

Resultados

Inicialmente, realizou-se a descrição das variáveis nas três coortes estudadas, apresentada na Tabela 1. Comparando a coorte 1 com as coortes 2 e 3, nota-se um aumento nas porcentagens dos diagnósticos de aids: em mulheres, indivíduos heterossexuais e nas faixas etárias maiores que 30 anos. Em relação à escolaridade, a maioria dos diagnósticos ocorreu em indivíduos com escolaridade inferior a 8 anos. Houve uma redução da categoria de transmissão por uso de drogas injetáveis entre os pacientes notificados, o que pode ser visualizado nas coortes 2 e 3. Quando se avalia a variável óbito, observa-se uma queda drástica no número de mortes na coorte 3, quando comparada com as outras duas coortes. Nas três coortes, a maioria dos indivíduos notificados apresentava contagem de células CD4 menor que 350 células por mm3 no momento do diagnóstico.

Tabela 1:
Número e proporções de indivíduos diagnosticados com aids de acordo com as variáveis nas três coortes estudadas, Campinas, São Paulo.

Observa-se grande percentual de perdas de informações nos três períodos de estudo, como mostra a Tabela 1. Devido à coorte 1 apresentar 96% de falta de informação da variável CD4, compatível com a realidade dos serviços na primeira década da epidemia, optou-se por não incluí-la no ajuste do modelo da coorte 1, e também considerá-la como uma variável contínua nos ajustes dos modelos para maximizar a informação.

A Tabela 2 apresenta os resultados do ajuste do modelo semiparamétrico (Equações 1 e 2) com as estimativas posteriores (média e desvio padrão). A Tabela 3 mostra as razões de riscos associadas às variáveis do estudo e seus respectivos intervalos de credibilidade de 95% (IC95%). Na coorte 1, os efeitos das variáveis sexo, escolaridade e categoria de transmissão por uso de drogas injetáveis não foram significativos, pois a média posteriori igual a zero se encontra no IC95%. Nas coortes 2 e 3, a variável categoria de transmissão por uso drogas injetáveis não foi significativamente associada ao risco de morrer. Na coorte 3, também se observa que os efeitos das variáveis idade, na faixa de 30 a 49 anos, escolaridade e contagem de CD4 não foram significativos.

Tabela 2:
Estimativas posteriores (média, desvio padrão e intervalo de credibilidade de 95%) dos parâmetros do modelo semiparamétrico para as variáveis, nas coortes 1, 2 e 3.
Tabela 3:
Razão de risco de óbito e seus respectivos intervalos de credibilidade de 95% obtidos do ajuste do modelo semiparamétrico, para os indivíduos diagnosticados com aids, nas coortes 1, 2 e 3, Campinas, SP.

O risco de óbito, na coorte 3, para homens foi 1,91 vezes maior (IC95% 1,36 - 2,73) do que para as mulheres. Para os maiores de 49 anos, o risco de óbito foi 2,59 (IC95% 1,58 - 4,24) vezes o risco em menores que 30 anos. Observa-se que em todas as coortes, ter idade maior que 49 anos foi fator de risco de óbito. Para aqueles com categoria de transmissão homossexual/bissexual, o risco de óbito foi 0,48 (IC95% 0,31 - 0,72) vezes (ou 52% menor) que o risco entre os da categoria de transmissão heterossexual.

Os mapas apresentados na Figura 1 (A, B e C) mostram o risco de óbito por setor censitário na cidade de Campinas, indicando que quanto maior o valor na escala, maior foi o risco de óbito dos indivíduos notificados por aids, nos períodos de estudo.

Figura 1:
Mapas dos efeitos espaciais de área para os setores censitários (médias posteriori) para a mortalidade por aids nas coortes estudadas (A, B, C) e mapa de classes de carência socioambiental (D, obtida do artigo Costa et al.1616. Costa JV, Donalisio MR, Silveira LV de A. Spatial distribution of dengue incidence and socio-environmental conditions in Campinas, São Paulo State, Brazil, 2007. Cad Saúde Pública. 2013; 29(8): 1522-32. DOI: 10.1590/0102-311X00110912
https://doi.org/10.1590/0102-311X0011091...
) da cidade de Campinas, São Paulo.

Observa-se na coorte 1, que as regiões sudoeste, sul e central da cidade apresentaram maiores riscos (Figura 1A). Na coorte 2 (Figura 1B), as regiões sudoeste, noroeste, sul e central exibiram maiores riscos de óbito. Por fim, na coorte 3 (Figura 1C), o risco de óbito foi maior entre as regiões norte, noroeste, parte da sudoeste e na região central. Os mapas apresentam diferenças na distribuição espacial do risco de morrer, porém predomina o maior risco em áreas com indicadores socioambientais precários, como mostra a Figura 1D. A Figura 2 apresenta os gráficos de risco basal do modelo ajustado, em cada umas das coortes. Destaca-se que o risco de óbito diminui ao longo tempo, nas coortes 2 e 3, e se torna próximo de 0 em torno de 2 mil dias de seguimento.

Figura 2:
Média posteriori (linha sólida) e intervalo de credibilidade de 95% (linha tracejada) para o logaritmo do risco basal do modelo ajustado para a mortalidade.

DISCUSSÃO

Este estudo evidencia a queda crescente do risco de morrer entre os pacientes notificados com aids e moradores do município de Campinas em três períodos. A mortalidade entre os pacientes diminuiu ao longo do tempo, como também tem sido apontado em outros estudos no Brasil e em outras partes do mundo33. Lima VD, Hogg RS, Harrigan PR, Moore D, Yip B, Wood E, et al. Continued improvement in survival among HIV-infected individuals with newer forms of highly active antiretroviral therapy. AIDS. 2007; 21(6): 685-92. DOI: 10.1097/QAD.0b013e32802ef30c
https://doi.org/10.1097/QAD.0b013e32802e...
,44. Schneider MF, Gange SJ, Williams CM, Anastos K, Greenblatt RM, Kingsley L, et al. Patterns of the hazard of death after AIDS through the evolution of antiretroviral therapy: 1984-2004. AIDS. 2005. DOI: 10.1097/01.aids.0000189864.90053.22
https://doi.org/10.1097/01.aids.00001898...
,66. Guibu IA, Barros MB de A, Donalísio MR, Tayra Â, Alves MCGP. Survival of AIDS patients in the Southeast and South of Brazil: analysis of the 1998-1999 cohort. Cad Saúde Pública. 2011; 27(Suppl 1): S79-92. http://dx.doi.org/10.1590/S0102-311X2011001300009
http://dx.doi.org/10.1590/S0102-311X2011...
,1717. Burch LS, Smith CJ, Phillips AN, Johnson MA, Lampe FC. Socioeconomic status and response to antiretroviral therapy in high-income countries. AIDS. 2016; 30(8): 1147-62. DOI: 10.1097/QAD.0000000000001068
https://doi.org/10.1097/QAD.000000000000...
. Com o acesso à terapia antirretroviral altamente ativa (TARV) a partir de 1996 devido à disponibilização pelo Ministério da Saúde, obteve-se o controle mais efetivo da proliferação viral e da recuperação do sistema imune dos pacientes.

Além do impacto na mortalidade, a morbidade dos pacientes com infecção pelo HIV tem diminuído progressivamente, elevando a sobrevida e a qualidade de vida dos pacientes infectados. Contudo, a aids continua sendo uma doença com grande impacto na morbidade e mortalidade e exige ações de prevenção e assistência contínuas.

Os estimadores das análises nas coortes estudadas identificam os homens, os maiores de 50 anos e os com categoria de transmissão heterossexual como os de maior risco de morrer nos períodos de estudo em Campinas.

Vários autores têm associado a maior sobrevida dos pacientes com aids com a idade ao diagnóstico, sexo, situação socioeconômica e categoria de exposição, como observado neste estudo em Campinas88. Teixeira TR de A, Gracie R, Malta MS, Bastos FI. Social geography of AIDS in Brazil: identifying patterns of regional inequalities. Cad Saúde Pública. 2014; 30(2): 259-71. DOI: 10.1590/0102-311X00051313
https://doi.org/10.1590/0102-311X0005131...
,1717. Burch LS, Smith CJ, Phillips AN, Johnson MA, Lampe FC. Socioeconomic status and response to antiretroviral therapy in high-income countries. AIDS. 2016; 30(8): 1147-62. DOI: 10.1097/QAD.0000000000001068
https://doi.org/10.1097/QAD.000000000000...
,1818. Arnold M, Hsu L, Pipkin S, McFarland W, Rutherford GW. Race, place and AIDS: The role of socioeconomic context on racial disparities in treatment and survival in San Francisco. Soc Sci Med. 2009. DOI: 10.1016/j.socscimed.2009.04.019
https://doi.org/10.1016/j.socscimed.2009...
,1919. Barboza R. Homens idosos e o HIV/Aids no campo da Saúde Coletiva: vulnerabilidades e desafios na quarta década da epidemia. BIS Bol Inst Saúde. 2012; 14: 81-9.. As condições socioeconômicas dos pacientes podem influenciar o risco de morrer devido a diferenças no acesso ao seguimento clínico, prevenção de complicações e na adesão aos antirretrovirais, que são menores entre as populações mais pobres44. Schneider MF, Gange SJ, Williams CM, Anastos K, Greenblatt RM, Kingsley L, et al. Patterns of the hazard of death after AIDS through the evolution of antiretroviral therapy: 1984-2004. AIDS. 2005. DOI: 10.1097/01.aids.0000189864.90053.22
https://doi.org/10.1097/01.aids.00001898...
,1111. Pellini ACG, Bellini JH, Cavalin RF, Chiaravalloti Neto F, Zanetta DMT. Distribuição espacial das notificações de aids em mulheres no município de São Paulo nos períodos de 1999-2001 e 2009-2011: uma análise sob a ótica da desigualdade socioespacial. Rev Espaço Geogr [Internet]. 2015; 18(2). Disponível em: http://www.lsie.unb.br/espacoegeografia/index.php/espacoegeografia/article/view/441 (Acessado em 05 de janeiro de 2016).
http://www.lsie.unb.br/espacoegeografia/...
,1717. Burch LS, Smith CJ, Phillips AN, Johnson MA, Lampe FC. Socioeconomic status and response to antiretroviral therapy in high-income countries. AIDS. 2016; 30(8): 1147-62. DOI: 10.1097/QAD.0000000000001068
https://doi.org/10.1097/QAD.000000000000...
,1818. Arnold M, Hsu L, Pipkin S, McFarland W, Rutherford GW. Race, place and AIDS: The role of socioeconomic context on racial disparities in treatment and survival in San Francisco. Soc Sci Med. 2009. DOI: 10.1016/j.socscimed.2009.04.019
https://doi.org/10.1016/j.socscimed.2009...
. Os mapas apresentados mostram piores estimativas de sobrevida nas regiões sul e sudoeste do município principalmente nas primeiras coortes, restando a região sudoeste como a de maior risco na coorte 2001-2005 em todos os mapas, regiões coincidentes com áreas de elevada carência socioambiental1616. Costa JV, Donalisio MR, Silveira LV de A. Spatial distribution of dengue incidence and socio-environmental conditions in Campinas, São Paulo State, Brazil, 2007. Cad Saúde Pública. 2013; 29(8): 1522-32. DOI: 10.1590/0102-311X00110912
https://doi.org/10.1590/0102-311X0011091...
.

A distribuição espacial do risco entre os setores censitários em Campinas mostra diferenças entre as coortes. Embora em escalas decrescentes de risco, os mapas indicam que as regiões mais pobres da cidade (noroeste, sudoeste) sempre foram as com maior mortalidade associada à aids. As possíveis justificativas são as precárias condições de vida, de acesso aos serviços ou o maior percentual de uso de drogas nas primeiras décadas da epidemia.

O “envelhecimento” da epidemia e o pior prognóstico em maiores de 50 anos têm sido também identificado por outros autores2020. Sousa ACA, Suassuna DSB, Costa SML. Perfil Clínico Epidemiológico de Idosos com AIDS. J Bras Doenças Sex Transm. 2009; 21(1): 22-6.. Ressalta-se a necessidade de reconhecer a vulnerabilidade à infecção nessa faixa de idade, considerando-se a menor percepção de risco, a falta de conhecimento e de utilização de métodos de prevenção, além do diagnóstico tardio nessa faixa etária1919. Barboza R. Homens idosos e o HIV/Aids no campo da Saúde Coletiva: vulnerabilidades e desafios na quarta década da epidemia. BIS Bol Inst Saúde. 2012; 14: 81-9..

Algumas variáveis clínicas, não abordadas neste estudo, também têm sido apontadas como preditores de sobrevida de pacientes com aids, como as doenças oportunistas, comorbidades, acesso e adesão aos esquemas terapêuticos e a serviços especializados33. Lima VD, Hogg RS, Harrigan PR, Moore D, Yip B, Wood E, et al. Continued improvement in survival among HIV-infected individuals with newer forms of highly active antiretroviral therapy. AIDS. 2007; 21(6): 685-92. DOI: 10.1097/QAD.0b013e32802ef30c
https://doi.org/10.1097/QAD.0b013e32802e...
,44. Schneider MF, Gange SJ, Williams CM, Anastos K, Greenblatt RM, Kingsley L, et al. Patterns of the hazard of death after AIDS through the evolution of antiretroviral therapy: 1984-2004. AIDS. 2005. DOI: 10.1097/01.aids.0000189864.90053.22
https://doi.org/10.1097/01.aids.00001898...
,66. Guibu IA, Barros MB de A, Donalísio MR, Tayra Â, Alves MCGP. Survival of AIDS patients in the Southeast and South of Brazil: analysis of the 1998-1999 cohort. Cad Saúde Pública. 2011; 27(Suppl 1): S79-92. http://dx.doi.org/10.1590/S0102-311X2011001300009
http://dx.doi.org/10.1590/S0102-311X2011...
,1717. Burch LS, Smith CJ, Phillips AN, Johnson MA, Lampe FC. Socioeconomic status and response to antiretroviral therapy in high-income countries. AIDS. 2016; 30(8): 1147-62. DOI: 10.1097/QAD.0000000000001068
https://doi.org/10.1097/QAD.000000000000...
. A maioria dos indivíduos com contagem CD4 ≤ 350 células por mm3 no momento do diagnóstico sugere o diagnóstico tardio da doença. A associação negativa da contagem de CD4 na primeira consulta (variável contínua) com a mortalidade na coorte 2 reforça esse parâmetro laboratorial como preditor de evolução clínica dos casos. Entretanto, na coorte 3, outras variáveis deslocaram o CD4 para explicar o risco de morrer, entre elas a idade e o sexo. O acesso à medicação antirretroviral combinada nos serviços de referência na cidade possivelmente diminuiu a influência desse indicador laboratorial como preditor de morte. Campinas conta com dois centros municipais de referência para DST/aids, um centro de orientação e apoio sorológicos, além de ambulatórios de referência do Hospital das Clínicas da Universidade Estadual de Campinas (Unicamp) e no Hospital Celso Pierro da Pontifícia Universidade Católica (PUC) de Campinas, com cobertura ampla de indivíduos com HIV/aids do município e região metropolitana.

Algumas limitações do estudo podem ser apontadas, como as decorrentes da utilização de dados secundários para a análise. A notificação do caso nem sempre ocorre no momento do diagnóstico da aids; muitos casos são detectados tardiamente, outros no momento de intercorrências clínicas ou do início do tratamento. A exclusão de pacientes notificados como aids na data do óbito procurou reduzir essas distorções. A subnotificação de casos de aids e óbitos pode ocorrer devido à falta de diagnóstico, bem como ao mau preenchimento de atestados de óbito. Além disso, as fichas de notificação nem sempre são bem preenchidas, restando informações às vezes imprecisas. Porém, trata-se de um banco de dados com informações valiosas sobre as tendências da doença no território e no tempo. Por outro lado, o longo período de seguimento das três coortes implica na comparação de coortes em diferentes contextos de risco, opções terapêuticas e de prevenção, com impacto no prognóstico dos pacientes. Os períodos de estudo certamente são heterogêneos, pois resultam de conjunturas epidemiológicas distintas.

CONCLUSÕES

O estudo contribui para possíveis ações e monitoramento dos serviços de saúde do município de Campinas em relação aos grupos e às regiões da cidade de maior risco de morrer, entre elas as de maior carência e piores condições socioambientais. Ações de prevenção e controle da doença e particularmente medidas para aumentar a adesão dos pacientes ao seguimento clínico e ao uso correto da medicação disponível podem ser melhor direcionadas. Mesmo diante da diminuição do risco de morrer por aids nos últimos anos, após o acesso à TARV, a maior incidência se concentra entre os homens, mais idosos, na categoria de transmissão heterossexual e entre os mais pobres nas três coortes estudadas. Os mapas obtidos e os fatores de risco orientam possíveis ações e monitoramento da doença no município.

REFERÊNCIAS

  • 1
    Fonseca MGP, Bastos FI. Twenty-five years of the AIDS epidemic in Brazil: principal epidemiological findings, 1980-2005. Cad Saúde Pública. 2007. DOI: 10.1590/S0102-311X2007001500002
    » https://doi.org/10.1590/S0102-311X2007001500002
  • 2
    Fordyce EJ, Singh TP, Nash D, Gallagher B, Forlenza S. Survival rates in NYC in the era of combination ART. J Acquir Immune Defic Syndr. 2002; 30(1): 111-8.
  • 3
    Lima VD, Hogg RS, Harrigan PR, Moore D, Yip B, Wood E, et al. Continued improvement in survival among HIV-infected individuals with newer forms of highly active antiretroviral therapy. AIDS. 2007; 21(6): 685-92. DOI: 10.1097/QAD.0b013e32802ef30c
    » https://doi.org/10.1097/QAD.0b013e32802ef30c
  • 4
    Schneider MF, Gange SJ, Williams CM, Anastos K, Greenblatt RM, Kingsley L, et al. Patterns of the hazard of death after AIDS through the evolution of antiretroviral therapy: 1984-2004. AIDS. 2005. DOI: 10.1097/01.aids.0000189864.90053.22
    » https://doi.org/10.1097/01.aids.0000189864.90053.22
  • 5
    Chequer P, Hearst N, Hudes ES, Castilho E, Rutherford G, Loures L, et al. Determinants of survival in adult Brazilian AIDS patients, 1982-1989. The Brazilian State AIDS Program Co-Ordinators. AIDS. 1992.
  • 6
    Guibu IA, Barros MB de A, Donalísio MR, Tayra Â, Alves MCGP. Survival of AIDS patients in the Southeast and South of Brazil: analysis of the 1998-1999 cohort. Cad Saúde Pública. 2011; 27(Suppl 1): S79-92. http://dx.doi.org/10.1590/S0102-311X2011001300009
    » http://dx.doi.org/10.1590/S0102-311X2011001300009
  • 7
    Tancredi MV, Waldman EA. Survival of AIDS patients in Sao Paulo-Brazil in the pre- and post-HAART eras: a cohort study. BMC Infect Dis. 2014. DOI: 10.1186/s12879-014-0599-8
    » https://doi.org/10.1186/s12879-014-0599-8
  • 8
    Teixeira TR de A, Gracie R, Malta MS, Bastos FI. Social geography of AIDS in Brazil: identifying patterns of regional inequalities. Cad Saúde Pública. 2014; 30(2): 259-71. DOI: 10.1590/0102-311X00051313
    » https://doi.org/10.1590/0102-311X00051313
  • 9
    Sousa AIA de, Pinto VL. Análise espacial e temporal dos casos de aids no Brasil em 1996-2011: áreas de risco aumentado ao longo do tempo. Epidemiol Serv Saúde. 2016; 25(3): 467-76. DOI: 10.5123/S1679-49742016000300003
    » https://doi.org/10.5123/S1679-49742016000300003
  • 10
    Santos NJS, Tayra A, Silva SR, Buchalla CM, Laurenti R. A aids no Estado de São Paulo: as mudanças no perfil da epidemia e perspectivas da vigilância epidemiológica. Rev Bras Epidemiol. 2002; 5(3): 286-310. DOI: 10.1590/S1415-790X2002000300007
    » https://doi.org/10.1590/S1415-790X2002000300007
  • 11
    Pellini ACG, Bellini JH, Cavalin RF, Chiaravalloti Neto F, Zanetta DMT. Distribuição espacial das notificações de aids em mulheres no município de São Paulo nos períodos de 1999-2001 e 2009-2011: uma análise sob a ótica da desigualdade socioespacial. Rev Espaço Geogr [Internet]. 2015; 18(2). Disponível em: http://www.lsie.unb.br/espacoegeografia/index.php/espacoegeografia/article/view/441 (Acessado em 05 de janeiro de 2016).
    » http://www.lsie.unb.br/espacoegeografia/index.php/espacoegeografia/article/view/441
  • 12
    Martino S, Akerkar R, Havard R. Approximate Bayesian Inference for Survival Models. Scand J Stat. 2011; 38(3): 514-28. DOI: 10.1111/j.1467-9469.2010.00715.x
    » https://doi.org/10.1111/j.1467-9469.2010.00715.x
  • 13
    R Development Core Team. R: A language and environment for statistical computing [Internet]. Viena: R Found Stat Comput Vienna; 2017. Disponível em: http://www.r-project.org/ (Acessado em 10 de março de 2017).
    » http://www.r-project.org/
  • 14
    Besag J, Kooperberg C. On conditional and intrinsic autoregressions. Biometrika. 1995; 82(4): 733-46. https://doi.org/10.1093/biomet/82.4.733
    » https://doi.org/10.1093/biomet/82.4.733
  • 15
    Rue H, Martino S, Chopin N. Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations. J R Stat Soc Ser B Stat Methodol. 2009; 71(2): 319-92.
  • 16
    Costa JV, Donalisio MR, Silveira LV de A. Spatial distribution of dengue incidence and socio-environmental conditions in Campinas, São Paulo State, Brazil, 2007. Cad Saúde Pública. 2013; 29(8): 1522-32. DOI: 10.1590/0102-311X00110912
    » https://doi.org/10.1590/0102-311X00110912
  • 17
    Burch LS, Smith CJ, Phillips AN, Johnson MA, Lampe FC. Socioeconomic status and response to antiretroviral therapy in high-income countries. AIDS. 2016; 30(8): 1147-62. DOI: 10.1097/QAD.0000000000001068
    » https://doi.org/10.1097/QAD.0000000000001068
  • 18
    Arnold M, Hsu L, Pipkin S, McFarland W, Rutherford GW. Race, place and AIDS: The role of socioeconomic context on racial disparities in treatment and survival in San Francisco. Soc Sci Med. 2009. DOI: 10.1016/j.socscimed.2009.04.019
    » https://doi.org/10.1016/j.socscimed.2009.04.019
  • 19
    Barboza R. Homens idosos e o HIV/Aids no campo da Saúde Coletiva: vulnerabilidades e desafios na quarta década da epidemia. BIS Bol Inst Saúde. 2012; 14: 81-9.
  • 20
    Sousa ACA, Suassuna DSB, Costa SML. Perfil Clínico Epidemiológico de Idosos com AIDS. J Bras Doenças Sex Transm. 2009; 21(1): 22-6.

  • Fonte de financiamento: Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), processo nº 2014/22856-7

Datas de Publicação

  • Publicação nesta coleção
    17 Set 2018

Histórico

  • Recebido
    20 Out 2016
  • Revisado
    02 Ago 2017
  • Aceito
    10 Ago 2017
Associação Brasileira de Pós -Graduação em Saúde Coletiva São Paulo - SP - Brazil
E-mail: revbrepi@usp.br