Integrated GWAS and linkage mapping in a multi-parent population prioritizes candidate genes for kernel number per row in maize

기술 요약: 다계통 집단을 이용한 GWAS 및 연관 매핑 통합을 통한 옥수수 이삭당 알 수(Kernel Number per Row)의 후보 유전자 우선순위 선정

문제 제기
이삭당 알 수(KNPR)는 옥수수(Zea mays L.)의 단일 개체 수량을 결정하는 핵심 요소이지만, 특히 temperate non-Reid 계통 내에서 나타나는 자연적 변이를 조절하는 분자적 메커니즘과 유리한 대립유전자 변이에 대해서는 아직 명확히 규명되지 않았다. 전통적인 양친 집단 매핑 집단은 유전적 다양성과 매핑 해상도(통상적으로 메가베이스 단위)가 부족하여 특정 후보 유전자를 식별하는 데 한계가 있는 반면, 자연 집단에서의 전전장 유전체 연관 분석(GWAS)은 집단 구조로 인한 위양성(false positives) 문제에 직면할 수 있다. 따라서 복잡한 양적 형질인 KNPR을 해독하기 위해 통제된 교배의 통계적 검정력과 연관 분석의 높은 해상도를 결합하는 전략이 필요하다.

연구 방법론
본 연구는 이친 부모 연관 매핑(biparental linkage mapping)과 GWAS 사이의 간극을 메우기 위해 개발된 다계통 집단(Multi-Parent Population, MPP)을 사용하여 통합 매핑 전략을 채택하였다.

• 집단 구축: 공통 Reid 근계 부모(Ye107)와 유전적으로 구별되는 4개의 엘리트 계통(Q11, R-2-1-1, YML1218, Chang7-2)을 교배하여 단일 종자 하위 세대(F2:8) 방식으로 530개의 재조합 근교계(RILs)를 생성하였다. 이 설계는 Reid × Reid 및 Reid × non-Reid 헤테로틱 조합을 모두 포함하는 4개의 하위 집단(Pop1–Pop4)을 형성하였다.
• 표현형 분석: MPP를 세 가지 환경(2021년 Jinghong; 2022년 및 2023년 Yanshan)에서 평가하였다. KNPR을 측정하였으며, 환경 노이즈를 최소化하기 위해 최적 선형 불편 예측값(BLUPs)을 산출하였다. 광의적 유전력($H^2$)은 경험적 상한값으로서 추정되었다.
• 유전형 분석: 530개 RIL 전체에 대해 전전장 유전체 시퀀싱(WGS)을 수행하여 630만 개 이상의 고품질 SNP를 확보하였다. 결측 데이터와 미세 대립유전자 빈도(MAF)를 기준으로 필터링한 후 903,693개의 SNP를 유지하였다.
• 유전 분석:
  • GWAS: 집단 구조를 제어하기 위해 상위 3개의 주성분(PC)과 친족 계수(kinship matrix)를 포함하는 혼합 선형 모델(MLM)을 EMMAX에서 사용하였다. 두 가지 임계값 전략을 적용하였다: 엄격한 전전장 유전체 임계값($-\log_{10}(P) = 5.66$)과 미세 효과 로커스를 포착하기 위한 암시적 임계값($-\log_{10}(P) \geq 5.0$).
  • 연관 매핑: JoinMap을 사용하여 각 하위 집단에 대한 고밀도 유전 지도를 구축하였다. QTL Cartographer의 합성 구간 매핑(CIM)을 사용하여 LOD 임계값 2.5를 기준으로 QTL을 식별하였다.
  • 통합: GWAS 신호(암시적 임계값을 충족하는)가 QTL 구간과 공통으로 존재하는 게놈 영역을 식별함으로써 후보 유전자의 우선순위를 정하였다.
• 기능적 특성 분석: 우선순위가 선정된 후보들은 하플로타입 분석, 서열 정렬(비동의 치환 변이에 집중), 프로모터 cis-element 예측, 단백질 구조 모델링(homology modeling), 그리고 세 가지 이삭 발달 단계에서의 부모 계통 간 qRT-PCR을 통한 시공간적 발현 분석을 거쳤다. 또한, GBLUP을 사용하여 후보 하플로타입의 가산적 유전 효과를 평가하기 위해 유전적 추정 육종가(GEBVs)를 추정하였다.

주요 결과

• 표현형 변이: MPP는 여러 환경에서 KNPR에 대해 연속적이고 정규 분포에 가까운 변이를 보였으며, 광의적 유전력 추정치는 79.4%에서 93.3% 사이였다.
• 매핑 해상도: 연관 불평형(LD) 붕괴가 빠르게 진행되어 약 20 kb에서 평탄화되었으며, 이는 고해상도 후보 유전자 스크리닝을 가능하게 하였다.
• 공통 국소화(Colocalization): 통합 접근법을 통해 GWAS 신호와 QTL 구간이 일치하는 염색체 1, 3, 5, 7의 네 개 중첩 로커스를 식별하였다.
• 후보 유전자: 기능적 주석 및 서열 변이에 기반하여 네 개의 유전자를 우선순위화하였다:
  1. ZmNMD3 (Chr 1): 60S 리보솜 핵 수출 단백질을 인코딩한다. 비동의 치환 변이는 없었으나 프로모터에 NAC 결합 모티프를 보유하고 있었다. ZmNMD3의 Hap3 하플로타입은 non-Reid 하위 집단(Pop3)에서만 독점적으로 발견되었으며, 높은 KNPR 및 높은 적합 GEBV와 연관되어 있었다.
  2. ZmPAPS1 (Chr 5): 핵 poly(A) 폴리머라제를 인코딩한다. 높은 KNPR을 가진 부모인 R-2-1-1은 엑손 11에서 세 개의 고유한 미스센스 변이(p.D540E, p.S588C, p.C598S)를 보유하고 있었다.
  3. ZmNAC70 (Chr 3): NAC 전사 인자를 인코딩한다. 다수의 비동의 SNP를 포함하고 있었으나, 특정 하플로타입이 대조적인 KNPR 표현형을 가진 부모 간에 공유되어 복잡하거나 배경 의존적인 효과를 시사하였다.
  4. ZmF-box1 (Chr 7): 유비퀴틴 매개 분해에 관여하는 F-box 단백질을 인코딩한다. Ye107과 Chang7-2가 공유하는 미스센스 변이(p.L279V, p.A444E, p.K466Q)를 보유하고 있었다.
• 발현 패턴: qRT-PCR 분석 결과, ZmNMD3, ZmPAPS1, ZmNAC70, ZmF-box1은 이삭 원기 발달 단계 동안 단계 특이적 발현 패턴을 보였으며, 일반적으로 부모 계통의 KNPR 수준과 상관관}있었다.
• 육종가: GEBV 분석은 non-Reid 하위 집단(Pop3)이 유리한 ZmNMD3 Hap3의 존재와 일치하게 KNPR에 대해 더 높은 전전장 가산적 유전 값을 보유하고 있음을 나타냈다.

의의 및 주장
본 논문은 다계통 집단의 사용이 특히 양친 집단에서 놓치기 쉬운 temperate non-Reid 유전체에 숨겨진 귀중한 대립유전자 변이를 발굴하는 데 효과적임을 주장한다. GWAS와 연관 매핑을 통합함으로써, 본 연구는 넓은 QTL 구간을 특정 후보 유전자로 성공적으로 좁혔으며, ZmNMD3, ZmPAPS1, ZmNAC70, ZmF-box1을 KNPR과 관련된 새로운 로커스로 우선순위화하였다.

저자들은 이 후보들이 인과적 유전자가 아닌 통계적 연관성, 공통 국소화 및 기능적 주석에 기반하여 우선순위가 정해진 가설적 후보임을 명시적으로 밝히고 있다. 본 연구는 발견된 결과들이 정밀 육종을 위한 프레임워크를 제공한다는 점을 강조한다. 특히, non-Reid 유전체에서 ZmNMD3 Hap3를 식별한 것은 단일 개체 수량 개선을 위한 방향성 유전자 도입 및 다중 유전자 피라미딩을 위한 잠재적 경로를 제시한다. 그러나 저자들은 이러한 자원들이 상업적 옥수수 개량에 배치되기 전, 인과 기전과 육종 적용성을 확인하기 위해 기능적 검증(예: CRISPR/Cas9, NILs 및 반복 야지 시험)이 필요하다는 점을 언급하며 신중한 태도를 유지하고 있다.