예방접종 전 백일해 장기 동역학 모델 비교 연구

초록

본 논문은 백일해의 면역 특성을 다섯 가지 가정으로 나누어, 계절적 강제와 무작위성을 포함한 SIR 모델을 구축하고, 장기 유병률 시계열의 파워 스펙트럼을 분석한다. 각 모델의 스펙트럼 차이를 실제 전염병 기록과 비교함으로써 면역 소멸·보강·재감염 메커니즘을 검증한다.

상세 분석

이 연구는 백일해와 같이 어린 시절에 강력히 면역이 형성되는 질병의 장기 역학을 이해하기 위해, 면역 반응에 대한 다섯 가지 상이한 가정을 기반으로 한 최소 모델들을 제시한다. 모든 모델은 폐쇄된 인구를 전제로 하며, 계절적 전염력 변동을 사인 함수 형태로 강제하고, 전이 과정에 포아송 잡음을 도입한 stochastic SIR 구조를 갖는다. 첫 번째 모델은 회복 후 영구 면역을 가정하여 S→I→R→S 전이 없이 고정된 R 상태를 유지한다. 두 번째 모델은 면역이 일정 평균 기간 후 완전히 소멸하는 waning immunity를 도입해 R→S 전이가 가능하도록 설계하였다. 세 번째 모델은 면역 소멸과 동시에 면역 강화(immunity boosting) 메커니즘을 포함, 재노출 시 면역이 재강화되어 소멸 속도가 감소한다는 가정을 반영한다. 네 번째 모델은 회복자가 완전 면역이 아니라 부분 면역을 유지하며, 일정 확률로 재감염이 가능한 reinfection 모델이며, 재감염 시 감염력과 회복률이 원래와 동일하게 설정된다. 마지막 다섯 번째 모델은 부분 면역이 점진적으로 감소하면서 감염 위험이 증가하는 partial waning 모델로, 재감염 시 감염력과 전염 기간이 감소된 형태를 취한다. 각 모델에 대해 동일한 기본 전염 파라미터(기본 재생산수 R0, 평균 감염 기간, 인구 규모 등)를 적용하고, 3050년 규모의 시뮬레이션을 수천 번 수행해 평균 파워 스펙트럼을 추정하였다. 스펙트럼 분석 결과, 영구 면역 모델은 단일 주기(연간)와 그 고조 파형만을 보이며, 장기 저주파 성분이 거의 사라진다. 반면 면역 소멸 모델은 23년 주기의 저주파 피크와 연간 피크가 동시에 나타나, 실제 백일해 기록에서 관찰되는 2~3년 주기의 대폭발을 재현한다. 면역 강화가 포함된 모델은 저주파 피크가 약화되고, 고주파 성분이 상대적으로 강화돼, 면역 강화가 전염병 진폭을 억제한다는 점을 시사한다. 재감염 모델은 저주파 피크가 더욱 뚜렷해지고, 스펙트럼 전반에 걸쳐 폭넓은 주파수 대역에 에너지가 분산돼, 면역이 완전하지 않을 경우 전염병이 지속적으로 재발할 가능성을 보여준다. 마지막으로 부분 면역 소멸 모델은 저주파 피크와 고주파 피크가 모두 존재하지만, 저주파 피크의 위치가 약간 이동하고 폭이 넓어, 면역 감소 속도와 재감염 위험 사이의 복합적 상호작용을 반영한다. 이러한 차별화된 스펙트럼 특성은 실제 전염 기록과 비교해 모델 선택의 근거를 제공한다. 논문은 두 개의 공개된 전염 기록(19세기 영국과 20세기 미국)을 사용해 각 모델의 스펙트럼을 실제 데이터와 정량적으로 매칭시켰으며, 면역 소멸·재감염을 포함한 모델이 가장 높은 적합도를 보였음을 보고한다. 따라서 백일해와 같은 강면역 질환에서도 면역이 완전히 지속되지 않으며, 재노출에 의한 면역 강화와 부분 면역 소멸이 전염병 패턴을 결정짓는 핵심 요인임을 제시한다.