Spatio-temporal variation in the uptake of the Human Papilloma Virus (HPV) vaccine among Malawian girls between 2019 and 2024

논문 요약: 말라위 HPV 백신 접종률의 시공간적 변동 분석 (2019-2024)

1. 연구 배경 및 문제 제기 (Problem)

• 문제의 심각성: 말라위는 전 세계적으로 자궁경부암 발병률 (ASIR 10 만 명당 67.9 명) 과 사망률이 가장 높은 국가 중 하나입니다. 자궁경부암은 인간 유두종 바이러스 (HPV) 감염이 주요 원인입니다.
• 목표와 현실의 괴리: 세계보건기구 (WHO) 는 2030 년 자궁경부암 퇴치를 위해 15 세 미만 여아의 HPV 백신 완전 접종률을 90% 로 달성할 것을 권고하고 있습니다. 말라위는 2019 년 1 월 국가 HPV 백신 프로그램을 시작했으나, 접종률이 목표에 크게 미치지 못하고 있습니다.
• 연구의 필요성: 기존 연구들은 특정 지역 (Dowa, Salima 등) 에 국한된 소규모 분석에 그쳤으며, 전국 28 개 구 (District) 를 대상으로 한 체계적인 시공간적 (Spatio-temporal) 분석과 접종률에 영향을 미치는 환경적, 사회적 요인을 규명한 연구는 부재했습니다.

2. 연구 방법론 (Methodology)

• 데이터 소스: 2019 년부터 2024 년까지 말라위 보건부 (MOH) 로부터 수집된 28 개 구의 행정적 백신 접종 데이터 (1 차 및 2 차 접종) 를 사용했습니다. 대상은 9~14 세 여아입니다.
  • 데이터 전처리: 보고된 접종률이 100% 를 초과하는 오류 (1 차 6%, 2 차 19%) 가 있어, 이를 99% 로 보정 (Capping) 했습니다.
• 공간적 공변량 (Covariates) 선정:
  • 사회경제적/접근성: 야간 조명 강도 (NTL, 도시화 및 개발 수준 대리 변수), 도로까지의 거리.
  • 기후/환경: 연평균 강수량, 최고 기온, 식생 지수 (NDVI), 고도 (고도와 기온의 다중공선성으로 인해 고도 제외).
  • 역학/행동: HIV 유병률, 십대 임신 비율 (성적 행동 및 건강 서비스 접근성 대리 변수).
  • 시간적 변수: 초기 (2019-2020), 중기 (2021-2022, 팬데믹 기간), 후기 (2023-2024) 로 구분.
• 통계 모델링:
  • 모델: R-INLA (Integrated Nested Laplace Approximation) 패키지를 활용한 시공간 베이지안 모델을 사용했습니다.
  • 구조: 이항 분포 (Binomial distribution) 를 가정하고, 고정 효과 (covariates) 와 무작위 효과 (구 수준, 구 - 연도 상호작용) 를 포함하는 로짓 (Logit) 모델을 구축했습니다.
  • 시간적 추세: 비선형 시간 추세를 포착하기 위해 다양한 함수 (Broken-stick, Penalised spline, Cubic smoothing spline) 를 비교하여 WAIC (Widely Applicable Information Criterion) 를 기준으로 최적 모델을 선택했습니다.
  • 공간 상관관계: Empirical Bayes Index (EBI) 를 통해 잔차 공간 상관관계를 검증한 결과, 유의미한 공간 상관관계가 없어 구조적 공간 무작위 효과 (ICAR 등) 는 최종 모델에서 제외했습니다.

3. 주요 결과 (Key Results)

• 접종률 현황:
  • 연구 기간 동안 1 차 접종률 평균은 46.83%, **2 차 접종률 평균은 32.44%**로 WHO 목표 (90%) 에 크게 미치지 못했습니다.
  • 2019-2020 년에는 높은 접종률을 보였으나, 2021-2022 년 (팬데믹 기간) 에 급격히 하락했으며, 2023-2024 년에는 부분적으로 회복되었으나 초기 수준에는 미치지 못했습니다.
• 시공간적 이질성:
  • 모든 28 개 구에서 연구 기간 내내 WHO 90% 목표를 달성한 지역은 없었습니다.
  • Balaka, Karonga, Nsanje 등 일부 지역은 두 접종 모두 45% 미만의 낮은 접종률을 지속적으로 보였습니다.
• 접종률에 영향을 미치는 요인 (Covariates):
  • 야간 조명 (NTL, 도시화): 음의 상관관계가 유의미했습니다. 즉, 도시화 수준이 높은 지역일수록 접종률이 낮았습니다 (1 차: β=-0.599, 2 차: β=-2.164). 이는 도시 지역에서의 학교 기반 프로그램 참여 저조나 백신 기피 현상과 관련이 있을 것으로 추정됩니다.
  • 도로 거리: 2 차 접종에서 양의 상관관계가 있었습니다. 도로에서 먼 지역 (외곽/시골) 일수록 오히려 접종률이 높았으며, 이는 NTL 결과와 일관되게 농촌 지역이 상대적으로 더 높은 접종률을 보임을 시사합니다.
  • HIV 유병률: 2 차 접종과 양의 상관관계가 있었습니다 (β=2.440). HIV 부담이 높은 지역은 관련 건강 서비스 (예: HIV 관리 프로그램 내 HPV 백신 통합) 를 통해 접종 접근성이 더 좋았을 가능성이 큽니다.
  • 십대 임신 비율: 2 차 접종과 음의 상관관계가 있었습니다 (β=-1.589). 십대 임신 비율이 높은 지역은 성적으로 활동적이거나 임신 중인 소녀들이 백신 접종 대상에서 제외되거나 접근성이 낮아 접종률이 떨어지는 것으로 분석됩니다.
  • 기후 요인: 강수량, 기온, 식생 지수의 증가가 접종률 감소와 연관되었습니다.

4. 주요 기여 및 의의 (Contributions & Significance)

• 최초의 전국적 분석: 말라위 HPV 백신 도입 이후 전국 28 개 구를 대상으로 한 최초의 시공간적 분석을 수행하여, 지역별 접종 격차 (Coldspots) 를 정량화했습니다.
• 역설적 발견: 일반적으로 도시화가 진전된 지역이 의료 접근성이 좋을 것이라는 통념과 달리, 말라위에서는 도시화 수준이 높을수록 접종률이 낮고, 오히려 외곽/농촌 지역에서 접종률이 높게 나타나는 역설적 패턴을 발견했습니다. 이는 학교 기반 접종 프로그램의 도시 지역 실패 요인을 규명하는 데 중요한 시사점을 줍니다.
• 정책적 함의:
  • 분산형 계획 (Decentralised planning): 전국 평균 수치에 의존하기보다 지역별 특성에 맞는 맞춤형 전략이 필요합니다.
  • 타겟팅 개입: 도시 지역 (특히 고밀도 주거지) 에 대한 모바일 outreach(외 Outreach) 프로그램 강화, 십대 임신이 높은 지역을 위한 성교육 및 백신 접근성 개선, HIV 클리닉과의 연계 강화가 필요합니다.
  • 팬데믹 영향: 팬데믹 기간의 접종률 급락은 백신 공급망 및 학교 폐쇄의 직접적인 영향을 받았음을 확인했습니다.

5. 결론

본 연구는 말라위의 HPV 백신 접종률이 WHO 목표에 크게 미치지 못하며, 도시화, 기후 조건, 팬데믹, 그리고 지역별 역학적 요인 (HIV, 십대 임신) 에 의해 복잡한 시공간적 패턴을 보임을 밝혔습니다. 자궁경부암 퇴치를 위해서는 전국적 차원의 균일한 접근보다는 지역별 맥락 (Context-specific) 을 고려한 분산형 전략과 도시 지역을 포함한 취약 계층 대상의 표적화된 개입이 시급히 요구됩니다.