An update to the critical vaccination threshold for measles herd immunity

기술 요약: 홍역 집단 면역을 위한 임계 백신 접종 문턱값의 업데이트

문제 제기
홍역 역학의 고전적 가정은 백신에 의한 면역이 평생 지속된다는 것이며, 이에 따라 집단 면역을 달성하기 위한 임계 백신 접종 문턱값($p_c$)은 $1 - 1/R_0$가 된다. 그러나 최근의 증거들은 백신 유도 면역이 시간이 지남에 따라 약화될 수 있음을 시사하며, 이는 백신 접종 후 수십 년이 지난 개인을 다시 감수성 있는 상태로 만들 수 있다. 이러한 면역 약화는 기존의 공중 보건 지표, 특히 제어 재생산 수($R_c$)와 임계 백신 접종 문턱값($p_c$)에 대한 도전을 제기한다. 본 논문은 백신 면역이 영구적이지 않은 인구 집단에서 이러한 지표들을 재평가할 필요성을 다루며, 2024-2025년 캐나다 남서부 온타리오(Southwestern Ontario)의 홍역 발생 사례를 케이스 스터디로 활용한다.

방법론
저자들은 표준 SIR(Susceptible-Infected-Recovered) 프레임워크를 확장하여 SIRVWS 구조를 갖춘 결정론적 수학 모델을 개발하였다. 인구는 다음과 같이 구분된다:

• S: 완전한 감수성을 가진 개인.
• I: 감수성 계층으로부터 기원한 감염된 개인.
• R: 평생 면역을 가진 회복된 개인.
• V: 완전한 보호를 받는 백신 접종 개인.
• W: 면역이 약화된 단계(부분적 면역 상태)의 백신 접종 개인.
• $I_W$: 약화된 면역 계층으로부터 기원한 감염된 개인.

모델은 인구 역학(출생 및 사망)과 면역 약화($V \to W \to S$) 및 면역 강화($W \to V$, 노출 시)를 위한 특정 전이를 포함한다. 주요 파라미터로는 약화율($k_1, k_2$), 부분적 면역자의 감염성 감소율($G$), 면역 강화 비례 상수($\nu$), 그리고 노출 시 감염 대 면역 강화의 확률($q$)이 포함된다.

저자들은 차세대 행렬(Next Generation Matrix)법을 사용하여 다음의 분석적 식을 도출하였다:

1. 제어 재생산 수 ($R_c$).
2. 집단 면역을 위해 필요한 임계 백신 접종 문턱값 ($p_c$).
3. 발생 확률 및 최종 유행 규모.

이러한 이론적 결과들은 남서부 온타리오 공중 보건 지역(SOPH) 발생 사례에 적용되었다. 저자들은 라틴 하이퍼큐브 샘플링(LHS)을 사용하여 타당한 범위 내에서 100,000개의 파라미터 세트를 생성하여 유행 역학을 시뮬레이션하였으며, 보고된 사례 데이터(2025년 8월 기준 확진 사례 769건)에 모델을 적합시켜 $R_0$, 약화율, 전파에 관여하는 유효 인구 규모 등의 파라미터를 추정하였다.

주요 기여 및 결과

• 수정된 지표: 본 연구는 면역 약화를 명시적으로 고려하여 $R_c$와 $p_c$에 대한 업데이트된 방정식을 도출한다. 저자들은 고전적 공식($R_c = R_0(1-p)$ 및 $p_c = 1 - 1/R_0$)이 하한선 역할을 하며, 면역 약화를 고려할 경우 실제 문턱값은 이보다 더 높다는 것을 입증한다.
• 집단 면역의 불가능성: 주요 발견 중 하나는 매우 느린 면역 약화(면역 반감기 77~99년에 해당하는 수준)조차도 집단 면역 역학을 근본적으로 변화시킬 수 있다는 것이다. 분석 결과, 특정 파라미터 영역에서는 높은 백신 접종률에도 불구하고 $p_c$가 1을 초과하게 되어 백신 유도 집단 면역을 달성하는 것이 불가능함을 확인하였다.
• 유행 확률 vs. 집단 면역: 면역 약화 조건하에서 집단 면역 달성이 불가능할 가능성이 있음에도 불구하고, 모델은 유형 재생산 수가 충분히 감소한다면(예: 약화된 계층의 낮은 감수성 또는 감소된 감염성) 유행 발생 확률은 여전히 매우 낮게 유지될 수 있음을 보여준다.
• SOPH 발생 사례 적용:
  • 모델은 SOPH 발생에 대한 기초 재생산 수($R_0$)를 약 12.3으로, 제어 재생산 수($R_c$)를 2.7로 추정하였다.
  • 추정된 백신 접종률은 약 86%였다.
  • 약화율은 $k_1 \approx 0.009$ year$^{-1}$ ($V \to W$), $k_2 \approx 0.007$ year$^{-1}$ ($W \to S$)로 추정되었으며, 이는 각각 약 77년과 99년의 면역 반감기에 해당한다. 이 수치는 감염 가능성이 백신 접종 후 약 51년 후에 돌아올 수 있음을 시사하는 최근의 혈청학적 연구들과 일치한다.
  • 전체 감염의 약 77%가 완전한 감수성 계층($S$)에서 발생하였으며, 이는 대다수의 사례가 미접종자에게서 발생했다는 보고와 일치한다.
  • 전파에 관여하는 유효 인구($N_0$)는 4,400명에서 6,500명 사이인 것으로 추정되었으며, 이는 전체 지역 인구의 극히 일부이다.
• 접종률에 대한 민감도: 78%와 92%의 백신 접종률을 비교한 시뮬레이션 결과, 접종률이 높을수록 (모델 적합 제약으로 인해) 더 높은 $R_0$ 추정치를 나타내지만, 전체 유행 규모는 작아지고 정점 도달 시간은 빨라졌다. 반대로, 낮은 접종률(78%)은 더 큰 최종 유행 규모와 약 6개월 늦은 정점 발생을 초래했다.

의의 및 주장
본 논문은 백신 유도 면역이 약화되는 고접종 인구 집단에서 널리 사용되는 홍역 집단 면역 문턱값을 재평가해야 한다고 주장한다. 저자들은 자신들의 업데이트된 수학적 정의가 평생 면역을 가정하는 고전적 모델보다 공중 보건 의사결정에 더 정확한 근거를 제공한다고 단언한다.

연구는 역설적인 역학 관계를 강조한다: 면역 약화로 인해 이론적인 집단 면역 문턱값을 충족할 수 없는 상황에서도, 부분적 면역자의 감소된 감수성과 감염성 덕분에 실제 대규모 유행의 위험은 낮게 유지될 수 있다. 저자들은 자신들의 모델이 위험 평가 도구로서 필요한 업데이트를 제공하지만, 균질 혼합(homogeneous mixing) 가정에 의존하며 연령 구조 및 확률적 효과를 무시하고 있다는 점을 언급하며 이를 향후 연구를 위한 한계점으로 제시하였다. 본 연구는 새로운 백신 전략을 제안하는 것이 아니라, 백신 보호의 장기적 지속 가능성을 판단하고 유행 데이터를 해석하기 위한 정교한 수학적 프레임워크를 제공하는 데 목적이 있다.