Travelling waves of invasion in microbial communities with phenotypic switching

문제 제기
미생물 군집은 종종 클론성 세균 개체군이 서로 다른 표현형으로 분화되는 복잡한 공간 구조를 보입니다. 잘 혼합된 환경에서의 종 공존과 관련하여 이러한 하위 개체군 구조의 생태학적 중요성은 인식되고 있지만, 공간적 경쟁 역학에 미치는 영향은 여전히 덜 이해되어 있습니다. 구체적으로 본 논문은 한 세균 종이 공간적으로 확장된 서식지에서 다른 세균 종을 침입할 때 발생하는 침입 전파파(traveling waves)에 표현형 전환, 즉 확률적 전환과 경쟁자에 대한 반응적 전환이 어떻게 영향을 미치는지 조사합니다. 본 연구는 두 경쟁 종을 포함하는 최소 모델에 초점을 맞추는데, 이는 단일 표현형을 가진 종 A와 성장 표현형 (B) 과 경쟁에 견고한 지속자 (P) 표현형 사이를 전환할 수 있는 종 B 로 구성됩니다.

방법론
저자들은 두 종 간의 공간적 경쟁을 설명하기 위한 최소 반응 - 확산 모델을 개발했습니다. 이 모델은 성장 표현형 ($B$), 지속자 표현형 ($P$), 그리고 경쟁자 ($A$) 의 풍부도를 추적하는 세 개의 결합된 편미분 방정식으로 구성됩니다. 역동성은 다음을 포함합니다:

• 로지스틱 성장 및 Lotka-Volterra 경쟁 항.
• $B$와 $P$ 사이의 확률적 전환율.
• 경쟁자 $A$의 존재에 의해 유도된 $B$에서 $P$로의 반응적 전환.
• 세균의 운동성을 나타내는 확산 항.

분석은 두 가지 상호 보완적인 접근 방식을 통해 진행됩니다:

1. 해석적 유도: 저자들은 이동 좌표 $z = x - ct$를 사용하여 시스템을 일련의 상미분 방정식으로 변환합니다. 그들은 (한 종만 존재하는) 균일한 정상 상태 근처의 선형화된 역학을 분석하여 이종 연결 (traveling waves) 의 존재에 필요한 대수적 조건을 유도합니다. 이는 파동의 앞쪽과 뒤쪽 가장자리에서 야코비안 행렬의 고윳값을 계산하여 음이 아닌 해에 필요한 최소 파동 속도 ($c_{min}$) 를 결정하는 것을 포함합니다.
2. 수치 시뮬레이션: 전체 반응 - 확산 시스템과 해당 "평균화" 모델 (표현형 하위 개체군을 단일 유효 개체군으로 축소) 을 MATLAB 에서 수치적으로 적분합니다. 저자들은 매개변수 공간을 무작위로 샘플링하여 역학에 의해 선택된 실제 파동 속도 ($c^*$) 를 해석적으로 유도된 최소 속도와 평균화 모델이 예측한 속도와 비교합니다.

주요 기여 및 결과
이 논문은 표현형 전환이 침입 속도에 미치는 영향과 관련하여 두 가지 주요하고 다소 반직관적인 결과를 도출합니다:

1. 경쟁자에 의한 침입 ($c > 0$): 경쟁자 종 (A) 이 전환 종 (B) 을 침입할 때, 표현형 전환의 존재는 이동 파동의 속도에 영향을 미치지 않습니다. 전환이 포함된 전체 모델의 최소 파동 속도는 전환이 없는 평균화 모델의 속도와 동일합니다. 이 결과는 반응적 전환율 ($\beta$) 에 관계없이 성립합니다.

2. 경쟁자에 대한 침입 ($c < 0$): 전환 종 (B) 이 경쟁자 (A) 를 침입할 때, 표현형 전환은 침입 파동의 속도를 증가시킬 수 있습니다. 구체적으로, 반응적 전환율이 너무 높지 않다면 전체 모델의 파동 속도는 평균화 모델이 예측한 속도를 초과할 수 있습니다. 이는 견고한 표현형으로 전환할 수 있는 능력이 경쟁자의 영토를 장악하는 속도를 가속화하는 공격적 메커니즘으로 작용할 수 있음을 시사합니다.

해석적 결과는 최소 파동 속도가 파동의 앞쪽 가장자리에서 선형화된 역학에 의해 결정됨을 확립합니다. $c > 0$인 경우, 전환 항은 최소 속도에 대한 최고차항 표현에 나타나지 않습니다. 반면 $c < 0$인 경우, 전환 항은 유효 성장 및 확산 매개변수를 수정하여 특정 조건에서 더 빠른 전파를 가능하게 합니다.

의의 및 주장
저자들은 종종 스트레스나 항생제 생존을 위한 방어 전략으로 간주되는 세균의 지속성이 공간적 경쟁에서 공격적 생태 전략으로 기능할 수 있다고 결론지었습니다. 지속자 표현형으로 전환함으로써 개체군은 역설적으로 경쟁자의 침입 속도를 가속화할 수 있습니다. 이 발견은 견고한 표현형의 "덩어리"가 버퍼로 작용하여 침입 파동을 늦출 것이라는 직관적 기대에 도전합니다.

이 논문은 잘 혼합된 개체군에서 표현형 변이가 공존을 안정화한다는 더 넓은 맥락에서 이 결과를 위치시키며, 동일한 메커니즘 (희귀 표현형으로 전환) 이 공간적 확장을 주도할 수도 있음을 시사합니다. 저자들은 그들의 모델이 1 차원이지만 토양과 같은 다공성 매체에서의 세균 이동과 관련이 있음을 지적합니다. 그들은 표현형 전환이 있는 시스템에서 파동 속도 선택 메커니즘을 미래 연구의 열린 수학 문제로 식별했는데, 특히 $\dot{X}=0$인 고전적 개체군 역학 모델에서 전환으로 인한 이주가 파동 속도 선택에 어떻게 영향을 미치는지에 관한 것입니다.