쥐에 묻힌 진드기 집합 패턴: 파워‑법칙 꼬리가 질병 전파를 촉진한다

본 논문은 유럽 전역에서 인간·동물 보건에 위협이 되는 진드기 매개 질병, 특히 틱-뇌염(TBE) 바이러스의 전파 메커니즘을 이해하기 위해, 이탈리아 트렌티노 북부 알프스 지역에서 9년간 수집된 장기간 데이터셋을 이용해 진드기 집합 패턴을 정량적으로 분석하고, 이를 역학 모델에 적용한 연구이다. 1. **데이터 수집 및 기본 통계** - 2000년부터 2008년까지 9년간 4,722마리의 Apodemus flavicollis(노란목쥐)를 포획·표시·재포획(CMR) 방식으로 조사하였다. - 각 쥐에 부착된 Ixodes ricinus 진드기의 수와 발달 단계(주로 유충, 소수의 nymph 및 성체)를 기록했으며, 총 55,411마리의 진드기가 수집되었다. - 진드기 부하는 평균 8마리, 중앙값 8마리, 최대 111마리(유충)·15마리(nymph)까지 관찰되었다. 2. **분포 모델링** - 진드기 부하의 꼬리를 설명하기 위해 부정이항(NB), 포아송‑로그정규(PoiLN), 그리고 이산형 파워‑법칙(PL) 세 가지 확률분포를 적용하였다. - 각 모델의 파라미터는 최대우도추정(MLE)으로 추정하고, Kolmogorov‑Smirnov(KS) 통계와 부트스트랩(p‑value)으로 적합도를 평가하였다. - 특히 ‘꼬리’ 영역(k ≥ k_min)에서 모델을 피팅하고, Vuong 비정규화 로그우도비(LLR) 검정을 통해 PL이 가장 높은 적합도를 보였으며, NB와 PoiLN은 고부하 숙주를 과소추정하는 경향을 보였다. 3. **계절적 변동 분석** - 쥐 개체수의 연간 변동을 2차 다항식으로 스무딩하고 정규화한 뒤, 정규화 값이 특정 임계치(θ) 이상인 기간을 ‘피크’ 시기로 정의하였다. - 피크와 비피크 시기의 진드기 부하 분포를 KS‑통계량으로 비교한 결과, 피크 시기에 꼬리가 더욱 두꺼워져 PL의 스케일 파라미터 α가 감소함을 확인하였다. 이는 숙주 밀도가 높을수록 진드기 부하가 더욱 불균등하게 집중된다는 의미이다. 4. **공동흡혈 전파 모델** - TBE와 같은 비체계적 전파(공동흡혈) 메커니즘을 고려한 확률적 전염 모델을 구축하였다. 모델은 숙주‑진드기 이분 그래프를 기반으로 하며, 각 숙주에 할당되는 진드기 수를 NB, PoiLN, PL, 혹은 실제 관측 분포(EMP) 중 하나로 샘플링한다. - 전염 파라미터는 감염된 진드기와 비감염 진드기가 동일 숙주에서 동시에 흡혈할 확률을 반영했으며, 시뮬레이션은 10,000번 반복 수행하였다. - 결과는 PL과 EMP이 거의 동일한 역학적 임계값(R₀)과 평형 감염율을 보이며, NB와 PoiLN은 상대적으로 높은 임계값과 낮은 감염율을 나타냈다. 특히 계절적 변동을 반영한 PL(피크 시)에서는 임계값이 더욱 낮아져 작은 초기 감염이라도 전염병이 확산될 가능성이 커졌다. 5. **생태·역학적 함의** - 진드기 부하의 ‘두꺼운 꼬리’(고부하 숙주) 존재가 공동흡혈 전파에 결정적 역할을 함을 확인하였다. 기존에 널리 사용되던 부정이항 모델은 평균 부하와 분산을 설명하는 데는 충분하지만, 고부하 숙주의 비율을 과소평가한다. - 따라서 질병 역학 모델링 시 PL과 같은 ‘두꺼운 꼬리’ 분포를 채택해야 실제 전파 위험을 정확히 추정할 수 있다. 이는 TBE와 같이 비체계적 전파가 주요 경로인 질병의 지속 가능성을 이해하고, 예방·통제 전략을 설계하는 데 중요한 통찰을 제공한다. 6. **결론** - 진드기·들쥐 간 부하 분포는 전통적인 NB보다 PL이 더 적합하며, 계절에 따라 꼬리 형태가 변한다. - PL 기반 모델은 비체계적 전파 메커니즘을 포함한 역학 시뮬레이션에서 낮은 전염 임계값과 높은 감염 지속성을 예측한다. - 이러한 결과는 진드기 매개 질병, 특히 TBE와 같은 질병의 위험 평가와 관리에 있어, ‘두꺼운 꼬리’ 분포를 고려한 정밀 모델링이 필요함을 강조한다.