Introdução: As Reações de estresse à vacinação (REV) consistem em manifestações físicas desencadeadas por fatores emocionais, como ansiedade, medo ou tensão diante da imunização que surgem não por causa do imunizante em si. Esses eventos não possuem uma causa fisiológica identificável. Embora geralmente benignas, essas reações podem impactar negativamente a confiança da população nas vacinas e gerar hesitação vacinal, prejudicando as metas de cobertura. Este estudo teve como objetivo identificar aglomerados espaço-temporais de REV registradas após a administração da vacina contra o HPV, no período de 2014 a 2024, a fim de subsidiar a tomada de decisões estratégicas no âmbito da vigilância em saúde no Brasil.
Métodos: Trata-se de um estudo observacional retrospectivo, com abordagem quantitativa, baseado em dados secundários de sistemas oficiais de vigilância (SI-PNI e e-SUS Notifica). Foram incluídos casos de REV em adolescentes de 9 a 14 anos vacinados com a vacina contra o HPV, bem como um grupo controle composto por outros ESAVI sem sintomas compatíveis com REV. Utilizou-se o software SaTScan™ para aplicar os modelos espaço-temporais de Poisson e de Bernoulli, com simulação de Monte Carlo e nível de significância de p ≤ 0,05, objetivando detectar clusters.
Resultados: No período analisado, mais de 108 milhões de doses da vacina contra o HPV foram administradas em adolescentes de 9 a 14 anos no Brasil, resultando em 3.503 notificações de Eventos Supostamente Atribuíveis à Vacinação ou Imunização (ESAVI), das quais 123 (0,11 por 100 mil doses) foram classificadas como REV. A maioria dos casos ocorreu em meninas (78,9%), concentrando-se nos primeiros anos após a introdução do imunizante, especialmente em 2014. A análise espaço-temporal realizada pelo modelo retrospectivo de Bernoulli no SaTScan™ identificou cinco aglomerados de REV, dos quais dois apresentaram significância estatística, com duração de 2 a 38 meses e ampla abrangência territorial; o cluster com maior risco relativo (RR: 8,56) localizou-se nas regiões Sudeste e Sul. Já a aplicação do modelo de Poisson identificou sete clusters, sendo três estatisticamente significativos, com duração entre 30 e 119 dias e destacando-se grandes agrupamentos nos anos de 2014 e 2015.
Conclusão: Até onde se sabe, este é o primeiro estudo a aplicar o SaTScan para detecção de aglomerados espaço-temporais de REV utilizando simultaneamente os modelos de Poisson (dados agregados) e Bernoulli (dados individuados). As REV ocorreram majoritariamente em meninas, com pico nos primeiros anos da introdução da vacina, e apresentaram manifestações leves, com evolução favorável. Esse padrão, compatível com a literatura internacional, reflete o impacto de fatores psicossociais no início das campanhas. A identificação de aglomerados espaço-temporais de REV permitiu detectar padrões de ocorrência no tempo e no território, contribuindo para o aprimoramento da resposta da vigilância. Esses achados reforçam o potencial da análise como ferramenta estratégica para decisões oportunas em saúde pública.