Climate and urbanization structure West Nile virus vector communities: Tools for surveillance under environmental change

문제 정의
북미 동부의 웨스트 나일 바이러스(WNV)는 주로 두 가지 잠재적(cryptic) 모기 종인 Culex pipiens와 Culex restuans에 의해 전파된다. 이 두 종은 형태학적 유사성 때문에 일상적인 감시 프로그램에서 하나의 "Cx. pipiens/restuans" 지수로 통합되어 관리되지만, 서로 구별되는 생물 계절학적, 열적, 서식지 선호도를 보인다. Culex restuans는 시즌 초반과 식생이 풍부한 서식지에서 우세한 반면, Culex pipiens는 여름이 늦게 진행될 때 정점에 도달하며 고도로 도시화된 환경에서 번성한다. 이러한 계절적 교체는 역학적으로 매우 중요하다. 왜냐 해가 갈수록 Cx. pipiens로의 전환이 인간의 WNV 위험이 가장 높은 시기와 일치하기 때문이다. 그러나 통합된 감시 지수는 환경 변화에 대한 이들의 상반된 반응을 포착하지 못한다. 기후 온난화와 도시 확장이 온도 및 토지 이용 구배를 변화시킴에 따라 매개체들의 상대적 풍부도가 변화하고 있으나, 기관들은 매개체 군집이 언제 어디서 재편될지 예측할 수 있는 종 수준의 해결된(species-resolved) 도구가 부족한 실정이다.

방법론
본 연구는 세 곳의 동부 대도시 지역의 데이터를 통합하였다: 볼티모어-워싱턴(BW) 지역의 다년도, 공간적으로 균형 잡힌 설계 연구(2019–2024)와 필라델피아 및 시카고의 기존 운영 컬렉션에서 추출한 보완적 데이터셋이다. 총 9,789마리의 Culex 표본이 분자적으로 식별되었다.

분석에는 두 가지 주요 통계적 접근 방식이 사용되었다:

1. 풍부도 모델링 (볼티모어-워싱턴 한정): 음이항 일반화 선형 혼합 모델(GLMM)을 사용하여 경관 예측 변수(예: 불투수면, 수관 밀도, 고도) 및 기후 예측 변수(예: 주간 평균 기온, 누적 적산온도)에 대한 종 특이적 풍부도를 모델링하였다. 두 종 사이의 발산적 반응을 테스트하기 위해 종 요인 모델(species-factor model)을 적합시켰다.
2. 구성적 전환 (세 지역 모두): 비선형적 구성적 전환을 정량화하기 위해 머신러닝 접근법인 그래디언트 포레스트(Gradient Forest, GF) 모델을 적용하였다. 이 모델들은 Cx. restuans에서 Cx. pipiens로 군집이 가장 급격히 변화하는 임계점을 식별하였다.
3. 공간 예측: BW 지역에서 GF 유래 생물학적 중요도 값을 주성분 분석(PCA)을 통해 변환하였고, 이를 베타-이항 일반화 선형 모델에 사용하여 초기 및 후기 계절 창(window) 동안의 예측된 종 구성을 경관 전체에 걸쳐 매핑하였다.

주요 결과

• 발산적 환경 반응: 두 종은 도시화와 온도 모두에 대해 반대 방향으로 반응하였다. Cx. pipiens의 풍부도는 불투수면이 높고 기온이 높을수록 증가한 반면, Cx. restuans의 풍부도는 두 구배를 따라 감소하였다. 두 종 모두 고도와 양의 상관관계를 보였으나, Cx. pipiens가 더 가파르게 반응하였다.
• 전환을 위한 열적 임계점: 누적 적산온도(ADD)와 주간 평균 기온은 세 지역 모두에서 종 전환의 가장 일관된 동인이었다. Cx. restuans 우세 군집에서 Cx. pipiana 우세 군집으로의 전환은 594 ~ 610 적산온도(ADD)(주간 평균 기온 약 21~22°C에 해당)라는 좁은 열적 창 내에서 가장 빠르게 발생하였다. 이 임계점은 샘플링 설계와 경관 맥락의 차이에도 불구하고 볼티모어-워싱턴, 시카고, 필라델피아에서 일관되게 나타났다.
• 경관 동인: 볼티모어-워싱턴의 경우, 경관 예측 변수(불투수면, 수관 밀도, NDVI)가 공간적 전환의 강력한 동인이었으며, 가장 급격한 변화는 교외 지역(불투수면 약 20~50%)에서 발생하였다. 필라델피아에서는 인구 밀도가 도시화 구배의 대리 변수로 작용하였고, 시카고에서는 제한된 경관 변동성으로 인해 기상 예측 변수가 전환 신호를 지배하였다.
• 공간 패턴: 공간 예측 결과, 고도로 도시화된 지역은 시즌 내내 Cx. pipiens가 풍부한 상태를 유지한 반면, 교외 경관은 군집 구성의 가장 큰 계절적 변화를 보였다.

의의 및 운영 도구
저자들은 이러한 발견이 환경 변화 하에서 대응적(reactive) 관리에서 예측적(anticipatory) 관리로 나아가기 위한 프레임워크를 제공한다고 주장한다. 본 연구는 방역 기관을 위해 세 가지 구체적인 운영 도구를 제안한다:

1. 기후 기반 예측: 식별된 적산온도 구간(594–610 ADD)은 전이 가능한 열적 임계점 역할을 한다. 기관들은 온도 데이터를 사용하여 트랩 구성에서 감지되기 전에 Cx.-pipiens 우세로의 전환을 예측함으로써, 감시 타이밍과 강도를 선제적으로 조정할 수 있다.
2. 공간적 타겟팅: 중간 정도의 불투수면을 가진 교외 경관이 가장 큰 계절적 전환 구역이라는 점은, 구성이 정체된 도시 중심부에 균일하게 집중하기보다 이 지역들을 집중 감시 및 조기 시즌 개입의 우선순위로 두어야 함을 시사한다.
3. 종 수준의 투영: 추정된 종 특이적 환경 반응을 통해, 기관들은 기존의 통합된 감시 데이터를 종 수준이 해결된 방식으로 해석할 수 있다. 이를 통해 모든 트랩에서 즉각적인 종 식별이 이루어지지 않더라도, 미래의 기후 온난화 및 토지 이용 변화에 따른 매개체 군집의 재편을 예측할 수 있다.

본 논문은 통합 지수가 수십 년간 감시 업무에 기여해 왔으나, 변화하는 환경 구배에 따른 매개체 군집의 재편을 추적하기에는 점점 더 불충분해지고 있다고 결론짓는다. 제안된 도구들은 공간적으로 타겟팅된, 미래 지향적인 WNV 위험 관리를 향한 경로를 제시한다.