Genetic Diversity and Population Structure of Madhuca longifolia in Tamil Nadu, India Revealed by SSR Markers

기술 요약: 인도 타밀나두 지역 마두카 롱티폴라(Madhuca longifolia)의 유전적 다양성 및 집단 구조

문제 제기
Madhuca longifolia(마후아)는 인도에서 식량, 연료, 사료 및 의약 자원을 제공하는 생태적, 사회경제적 가치가 매우 높은 다목적 수종이다. 광범위한 분포와 유용성에도 불구하고, 이 종은 분자 수준에서의 연구가 미흡한 상태이다. 기존 연구는 주로 표현형 평가에 의존해 왔는데, 이는 속도가 느리고 환경 의존적이며, 수명이 긴 수종에서 진정한 유전적 변이를 포착하는 데 한계가 있다. 또한, 특히 남인도 지역 규모에서의 유전 자원에 대한 포괄적인 분자적 특성 규명이 제한적이다. 고해고 해상도의 유전형 분석을 위한 강력한 도구가 부족하며, 보전 및 육종 프로그램을 지원하기 위한 집단 수준의 다양성, 유전 구조 및 유전자 흐름에 대한 이해가 필요하다.

연구 방법
본 연구는 인도 타밀나두의 세 가지 지리적으로 구분된 농업 기후 지역(남부 19개 계통, 중부 19개 계통, 북부 18개 계통)에서 수집된 56개의 M. longifolia 계통의 유전적 다양성과 집단 구조를 평가하였다. 샘플링에는 M. longifolia var. longifolia와 var. latifolia가 모두 포함되었다.

• 마커 개발: 연구진은 전사체 분석을 통해 개발된 15개의 다형성 단순 반복 서열(SSR) 마커를 활용하였다. 여기에는 11개의 EST-SSR 로커스(MLO 시리즈)와 4개의 게놈 SSR 로커스(MLA 시리즈)가 포함된다. 프라이머 설계는 디(di-) 및 트리(tri-)뉴클레오타이드 반복 모티프를 대상으로 하였다.
• 유전형 분석: 변형된 CTAB 프로토콜을 사용하여 게놈 DNA를 추출하였다. PCR 증폭은 특정 열 순환 조건을 사용하여 수행되었으며, 증폭 산물은 3% MetaPhor 아가로스 겔에서 분리되었다. 오류를 최소화하기 위해 두 명의 연구자가 독립적으로 대립유전자 크기를 판독하였다.
• 데이터 분석: 유전적 파라미터는 GenAlEx v.6.5 및 POPGENE v.1.32를 사용하여 계산되었다. 주요 지표로는 대립유전자 수($N_a$), 유효 대립유전자 수($N_e$), 관찰된($H_o$) 및 기대($H_e$) 이형접합도, 섀넌 정보 지수($I$), 다형성 정보 함량(PIC), 그리고 고정 지수($F$)가 포함된다. 유전자 흐름($N_m$)은 고정 지수($F_{ST}$)를 기반으로 추정되었다.
• 구조 분석: 네이(Nei)의 유전 거리와 동일성을 사용하여 집단 분화를 평가하였다. 분산 분석(AMOVA)을 통해 집단 간, 집단 내 개체 간, 그리고 개체 내 변이를 구분하였다. 클러스터 분석은 근연 관계 분석(NJ) 덴드로그램(MEGA v.11.0)과 주좌표 분석(PCoA)을 사용하여 수행되었다.

주요 결과

• 마커 성능: 15개의 SSR 마커 모두 높은 다형성과 일관된 증폭을 나타냈다. PIC 값은 0.608에서 0.666 사이로, 모든 마커가 매우 유익한 것으로 분류되었다.
• 유전적 다양성: 집단들은 상당한 대립유전자 풍부도를 보였다. 평균 기대 이형접합도($H_e$)는 0.657이었고, 섀넌 정보 지수($I$)는 1.080이었다.
  • 남부 집단: 양(+)의 고정 지수($F = 0.096$)와 낮은 관찰 이형접합도($H_o = 0.579$)를 보였으며, 이는 가벼운 근친교배 또는 비무작위 교배를 시사한다.
  • 중부 및 북부 집단: 각각 음(-)의 고정 지수($F = -0.136$ 및 $-0.102$)를 나타냈으며, 이는 이형접합자 과잉과 주로 타가수분 행동을 나타낸다.
• 집단 구조 및 분화:
  • AMOVA: 전체 유전 변이의 75%가 개체 내에 존재하고, 18%는 집단 간 차이에 기인하며, 6%는 집단 내 개체 간 차이에 기인함을 밝혔다. 이러한 분포는 광범한 유전자 흐름과 타가수분을 뒷받침한다.
  • 유전 거리: 네이의 유전 거리는 0.553(남부와 중부 사이)에서 1.319(중부와 북부 사이)까지 나타났으며, 이는 중부와 북부 집단 사이의 가장 높은 분기도를 나타낸다.
  • 유전자 흐름: 평균 유전자 흐름($N_m$)은 1.344였으며, 로커스별 값은 0.757에서 2.801까지였다. 높은 $N_m$은 유전적 부동을 상쇄하기에 충분한 중간 내지 높은 수준의 대립유전자 교환을 시사한다.
• 클러스터링: 덴드로그램 및 PCoA 분석 결과, 계통들은 주로 지리적 기원에 따라 그룹화되었으나 일부 혼합(admixture)이 관찰되었다. 중부와 북부 집단이 가장 큰 유전적 분기를 보였다.
• 정보성 로커스: 높은 $N_m$을 가진 MLO24와 높은 PIC 및 유전자 다양성을 가진 MLO33이 분자 지문법 및 집단 연결성 평가에 특히 유용한 로커스로 식별되었다.

의의 및 주장
본 논문은 이 연구가 타밀나두 지역 M. longifolia의 보전 및 유전적 개량을 위한 견고한 토대를 제공한다고 주장한다. 핵심 기여는 다음과 같다:

1. SSR 유용성 검증: 본 연구는 전사체 유래 SSR 마커가 이 비모델 수종의 분자적 특성 규명에 매우 효과적이며 높은 식별력을 제공함을 확인하였다.
2. 보전 전략: 연구 결과는 지역 내 M. longifolia 집단이 높은 유전적 다양성과 구조화되어 있으면서도 서로 연결된 구조를 가지고 있음을 시사한다. 보전 전략은 모든 대립유전자의 폭을 포착하기 위해 세 지리적 집단을 모두 대표하도록 우선순위를 두어야 하며, 유전적으로 분기된 계통에 특별히 주의를 기울여야 한다.
3. 육종 응용: 유전적으로 크게 분기된 집단(특히 중부 대 북부)과 정보성이 높은 로커스(MLO24 및 MLO33)의 식별은 정밀 육종의 기초를 제공한다. 저자들은 이러한 발견이 종자유 생산량이나 내건성과 같은 경제적으로 중요한 형질을 개선하기 위한 마커 보조 선발(MAS) 프로그램 및 잡종 강세 육종을 지원한다고 제안한다.
4. 미래 통합: 본 연구는 향-유전체 선택(GS) 및 전장 유전체 연관 분석(GWAS)으로 전환하기 위해 필요한 집단 구조 및 친연성 매트릭스를 구축한다.

저자들은 향후 영향에 대해 신중한 어조를 유지하며, 현재의 SSR 패널이 비용 효율적이고 휴대 가능하지만, 국지적 적응에 대한 선택적 스윕(selective sweep) 분석을 수행하기 위해 게놈 전장 재서열 분석(GBS 또는 전장 유전체 재서열 분석)과 같은 차세대 염기서열 분석과의 통합이 필요할 것이라고 언급하였다.