Genome-wide characterization of vegetable soybean core collection reveals genetic diversity for breeding efforts in Sub-Saharan Africa

기술 요약: 채소용 대두 핵심 컬렉션의 전장 유전체 특성 분석

문제 제기
채소용 대두(Glycine max (L.) Merr)는 사하라 이남 아프리카(SSA)에서 영양가를 높이고 작물 체계를 다양화할 수 있는 잠재력이 큰, 아직 충분히 활용되지 못한 콩과 작물이다. 이러한 가치에도 불구하고, 이 지역의 육종 노력은 포괄적인 유전적 특성 분석의 부재로 인해 제한되어 왔다. SSA 내 곡물용 대두 생산은 개량 품종을 통해 확대되었으나, 채소용 대두의 도입은 부적응, 낮은 수확 안정성, 도입 유전자원의 불충분한 농업적 성능으로 인해 제약을 받고 있다. 대두 다양성에 관한 기존 연구들은 종종 저밀도 분자 마커(예: SSRs, AFLPs)에 의존해 왔으며, 이는 현대 육종에 필요한 미세 규모의 집단 구조나 유전적 관계를 정의하기에는 해상도가 부족하다. 또한, 채소용 대로의 개선을 위한 핵심 자원인 세계 채소 센터(World Vegetable Center)의 대립형질 대립유전자(large-seeded) 대두 핵심 컬렉션 내 유전적 다양성, 집단 구조 및 진화적 관계에 대한 포괄적인 전장 유전체 기반 특성 분석이 이루어지지 않았다.

연구 방법론
본 연구는 다양한 기원을 가진 191개의 대두 계통 패널을 활용하였다:

• 세계 채소 센터의 채소용 대두 핵심 계통 109종.
• Seed Savers Exchange(미국)의 채소용 대두 26종.
• 국립농업연구기관(일본)의 채소용 대두 3종.
• SSA의 우량 곡물용 대두 53종 (IITA, CSIR, Makerere University 및 Seed Co Limited로부터 확보).

유전체 워크플로우:

1. 시퀀싱: Illumina NovaSeq 6000 플랫폼을 사용하여 제한 효소 결합 DNA(Restriction site-associated DNA, RAD) 시퀀싱을 수행하였다 (150 bp paired-end reads).
2. 변이 호출(Variant Calling): 원시 리드(Raw reads)는 Stacks (v2.68)를 사용하여 처리되었으며, BWA를 사용하여 Williams 82.a5 참조 유전체에 정렬되었다. 변이는 Best Practices 워크플로우에 따라 GATK (v4.2.6.1)를 사용하여 호출되었다.
3. 필터링 및 임퓨테이션(Imputation): 고신뢰도 SNP를 품질 임계값(QUAL > 30, DP > 10)을 기준으로 추출하였다. 누락된 유전형은 LD KNNi 방법을 사용하여 임퓨테이션하였다. 최종 필터링 결과, 유전형 호출율(GCR) ≥ 90% 및 미세 대립유전자 빈도(MAF) ≥ 5%를 충족하는 21,970개의 고품질 SNP 데이터셋을 확보하였다.
4. 통계 분석:
   • 집단 구조: 주성분 분석(PCA), 판별 분석(DAPC), 그리고 최적의 클러스터 수($K$)를 결정하기 위한 교차 엔트로피 최소화(LEA 패키지)를 통해 분석하였다.
   • 계통 발생: phangorn 패키지를 사용하여 근접 이웃법(Neighbor-joining) 트리를 구축하였다.
   • 차별성: 분자 분산 분석(AMOFA), 쌍별 고정 지수($F_{ST}$), 그리고 Nei의 유전적 거리를 계산하였다.
   • 다양성 지표: 추론된 하위 집단에 대해 다형성 정보 함량(PIC), 기대 이형접합도($H_e$), 유효 대립유전자 수($N_e$)를 산출하였다.

주요 결과

• 집단 구조: 전장 유전체 분석 결과, 컬렉션이 세 개의 주요 유전 그룹($K=3$)으로 일관되게 구분되었다.
  • 클러스터 1 & 2: 주로 채소용 대두로 구성됨.
  • 클러스터 3: 주로 SSA 우량 곡물용 대두로 구성됨.
  • $K=3$에서 121개 계통(63.35%)이 조상 계수($Q$) ≥ 0.7인 순수 클러스터로 할당되었다.
• 유전적 차별성:
  • AMOFA 분석 결과, 전체 유전 변이의 83.78%가 집단 내에 존재하며, 16.22%가 집단 간에 발생함이 밝혀졌다 ($\Phi_{ST} = 0.162, p=0.001$).
  • 쌍별 $F_{ST}$ 값은 중등도에서 강한 차별성을 나타냈으며, 클러스터 2와 클러스터 3 사이에서 가장 높은 격차($F_{ST} = 0.335$)를 보였고, 클러스터 1과 클러스터 3 사이에서 가장 낮은 격차($F_{ST} = 0.128$)를 보였다.
  • Nei의 유전적 거리는 클러스터 1과 3이 가장 가까운 진화적 관계를 공유함을 확인하였다 ($D = 0.043$).
• 유전적 다양성:
  • 클러스터 1은 가장 높은 유전적 다양성을 나타냈다 (PIC = 0.253; $H_e$ = 0.262).
  • 클러스터 3(SSA 곡물용 대두)은 가장 낮은 다양성을 보였다 (PIC = 0.188; $H_e$ = 0.178). 이는 제한된 창시자 계통으로부터의 반복적인 선발로 인한 상대적으로 좁은 유전적 기반을 시사한다.
  • 클러스터 2는 규모는 작지만(5개 계통), 중간 정도의 다양성과 높은 가산 유전 분산을 나타냈다.
  • 관찰된 이형접합도($H_o$)는 모든 클러스터에서 균일하게 낮았으며(0.009–0.019), 이는 대두의 자가 수분 특성과 일치한다.

의의 및 주장
본 논문은 세계 채소 센터의 대립형질 대(large-seeded soybean) 핵심 컬렉션이 사하라 이남 아프리카의 채소용 대두 개선을 위한 귀중한 유전적 다양성 저장고 역할을 한다고 주장한다. 본 연구의 주장은 다음과 같다:

1. 자원 식별: 식별된 집단 구조와 다양성 패턴은 유전자원 보존 및 부모 계통 선발을 위한 기초적인 프레임워크를 제공한다.
2. 육종 전략: 채소용 대두와 곡물용 대두 사이의 명확한 분리, 그리고 순계(pure lines)의 식별은 표적 부모 계통 선발의 기회를 제공한다. 유전적으로 차이가 큰 클러스터 간(예: 클러스터 1/2와 클러스터 3)의 교배는 대립유전자 재조합을 증가시키고 농업적 및 품질 형질을 위한 새로운 대립유전자의 도입을 용이하게 할 것으로 제안된다.
3. 유전체 기반: 생성된 SNP 데이터셋은 향후 마커 보조 선발, 유전체 선택, 그리고 사하라 이남 아프리카에 특화된 핵심 육종 집단 개발을 위한 토대를 마련한다.
4. 좁은 유전적 기반 해결: 본 연구 결과는 현재 SSA 곡물용 대두 품종(클러스터 3)에서 관찰되는 좁은 유전적 기반과 대조되는, 클러스터 1과 2에서 발견되는 다양한 유전자원을 전략적으로 활용함으로써 아프리카 육종 집단의 유전적 기반을 넓혀야 할 필요성을 강조한다.

저자들은 이 연구가 해당 특정 채소용 대두 자원에 대한 최초의 포괄적인 전장 유전체 특성 분석이며, 이 지역에 적응된 영양가가 높고 알이 큰 품종을 개발하기 위한 로드맵을 제공한다고 결론짓는다.