Integrative transcriptome and metabolome analysis identifies flavonoid biosynthesis as well as ascorbate and aldarate metabolism driving natural cold tolerance variation in dormant peach (Prunus persica)

기술 요약: 휴면기 복숭아의 내한성에 관한 통합 전사체 및 대사체 분석

문제 제기
겨울철 저온 스트레스는 복숭아(Prunus persica)의 성장, 월동 생존 및 지리적 분포를 제한하는 주요 요인이다. CBF 의존적 경로는 모델 식물의 저온 반응에 대해 잘 규명된 메커니즘이지만, 다양한 복숭아 품종 간의 내한성 자연 변이를 조절하는 분자적 메커니즘은 여전히 명확하지 않다. 기존 연구들은 주로 과실의 냉해(chilling injury)나 생육기 저온 내성에 초점을 맞추었으며, 일반적으로 제한된 유전자형 세트(통상적으로 내한성 품종 하나와 민감성 품종 하나)를 활용하였기에 내한성 메커니즘의 일반성을 포착하는 데 한계가 있었다. 따라서 잘 특성화된 표현형을 가진 광범위한 품종들을 사용하여, 내한성 차이의 근저에 있는 조절 네트워크를 밝힐 필요가 있다.

연구 방법
본 연구는 전사체학과 광범위 표적 대사체학을 결합한 통합적 "멀티 오믹스(multi-omics)" 접근법을 채택하였다.

1. 식물 재료 선정: 73종의 복숭아 품종 패널로부터 내한성 구배를 스크리닝하여 20종을 선별하였다. 현장 동해 조사 및 조직의 50% 사멸 온도(LT₅₀)의 실험실 측정을 바탕으로, 4종의 내한성 그룹(T 그룹: 'Hangongzhu', 'Huichun', 'Donghe No1', 'Qiuyue')과 4종의 내한성 민감 그룹(S 그룹: '21st Century', 'Zhonghuashoutao', 'Jinhuangjin', 'Zhongtao 13')을 포함한 총 8종을 선정하였다.
2. 실험 설계: T 그룹과 S 그룹의 1년생 휴면 가지를 두 가지 조건에 노출시켰다: 대조군 처리(4°C에서 12시간) 및 저온 스트레스 처리(-25°C에서 2시간 유지 후 4°C로 서서히 복귀).
3. 오믹스 분석:
   • 전사체 분석 (RNA-seq): 16개 샘플(4개 품종 × 2개 그룹 × 2개 처리) 전체에서 총 RNA를 추출하였다. Illumina 플랫폼을 통해 시퀀싱을 수행하였으며, Prunus persica v2.0.a1 참조 유전체에 정렬하였다. 차등 발현 유전자(DEG) 분석은 FDR < 0.05, |log₂FC| ≥ 1.5를 기준으로 수행되었다.
   • 대사체 분석: 저온 스트레스 조건의 샘플을 대상으로 UPLC-ESI-MS/MS를 이용한 광범위 표적 대사체 분석을 실시하였다. 차등 축적 대사물질(DAM)은 VIP > 1 및 |log₂FC| ≥ 1.0을 기준으로 선별하였다.
4. 데이터 통합: DEG와 DAM에 대해 KEGG 및 GO 농축 분석을 수행하였다. Cytoscape를 사용하여 DEG와 DAM 간의 상관관계 네트워크를 구축하였다. 핵심 후보 유전자는 qRT-PCR을 통해 검증하였다.

주요 결과

1. 보존된 저온 반응 경로: T 그룹과 S 그룹 모두 저온 스트레스에 반응하여 MAPK 신호 전달 경로와 식물 호르몬 신호 전달 경로의 보편적인 활성화를 보였다. 그러나 이러한 경로들은 내한성 그룹 간에 유의미한 차등 발현을 보이지 않았는데, 이는 이들이 유전형 간의 차이를 유도하는 요인이기보다는 일반적인 스트레스 반응임을 시사한다.
2. 플라보노이드 생합성:
   • 전사 수준: 대조 조건(저온 스트레스 없음) 하에서, T 그룹과 S 그룹 사이의 플라보노이드 생합성 경로에서 DEG의 유의미한 농축이 관찰되었다. 6개의 유전자가 확인되었으며, 특히 Prupe.3G185000(UDP-glycosyltransferase 71A16, UGT71A16 인코딩)은 T 그룹에서 유의하게 상향 조절된 반면, 다른 5개의 구조적 효소 유전자들은 하향 조절되었다.
   • 대사 수준: 저온 스트레스 하에서, 플라보노이드 생합성 경로의 DAM은 S 그룹에 비해 T 그룹에서 유의하게 농축되고 높은 수준으로 축적되었다.
   • 시사점: 본 연구는 T 그룹이 스트레스 전 단계에서 차별화된 기초 발현 패턴과 높은 플라보노이드 축적량을 유지하며, 이것이 우수한 활성산소(ROS) 제거 능력에 기여할 것임을 시사한다.
3. 아스코르브산 및 알다레이트 대사:
   • 전사 수준: 저온 스트레스 하에서, 이 경로는 T 그룹과 S 그룹 간의 DEG에서 유의하게 농축되었다. L-아스코르브산 산화효소(AO) 또는 AO 유사 단백질을 인코딩하는 3개의 유전자가 T 그룹에서 현저하게 하향 조절되었다. 특히, Prupe.2G077800은 S 그룹 대비 T 그룹에서 급격한 하향 조절(log₂FC = -9.84)을 보였다.
   • 대사 수준: 관련 대사물질인 L-아스코르브산-2-글루코사이드(AA2G)는 T 그룹에서 유의하게 상향 조절되었다.
   • 시사점: 내한성 품종에서의 AO 유전자 하향 조절은 아스코르브산(AsA)의 산화를 감소시켜, 결과적으로 더 큰 내인성 AsA 풀(pool)을 확보하고 항산화 능력을 강화하는 것으로 보인다.
4. 식물-병원균 상호작な: 이 경로는 대조군과 스트레스 조건 비교 모두에서 농축되었으나, 상응하는 대사체 데이터는 부족하였다. 몇몇 질병 저항성 단백질들이 차등 발현되었으며, 일부는 T 그룹에서 일관된 상향 조절을 보였다.
5. 후보 유전자: 플라보노이드 생합성 및 아스코르브산/알다레이트 대사 경로에서 9개의 유전자가 내한성 차이를 주도하는 핵심 후보로 식별되었다. qRT-PCR 검증을 통해 RNA-seq 데이터의 신뢰성을 확인하였다.

의의 및 주장
저자들은 본 연구가 휴면기 복숭아의 내한성 자연 변이에 대한 새로운 메커니즘적 통찰을 제공한다고 주장한다. 주요 기여는 보존된 스트레스 반응 메커니즘과 유전형 간의 차이를 유도하는 구성적 차이를 구분해 낸 점이다.

• 메커니즘적 통찰: 본 연구는 MAPK 및 호르몬 신호 전달이 저온에 의해 보편적으로 활성화되지만, 내한성 변이를 결정짓는 핵심 요인은 플라보노이드 생합성과 아스코르브산/알다레이트 대사의 구성적 차이라는 점을 상정한다.
• 특이적 동인: 본 연구는 내한성 품종이 플라보노이드 축적을 선호하는 독특한 기초 발현 프로파일과 L-아스코르브산 산화효소 유전자의 특이적 하향 조절을 보유하고 있으며, 이것이 ROS 제거를 위한 아스코르브산 풀을 보존한다는 것을 밝혀냈다.
• 육종 적용: 본 연구는 내한성 복숭아 품종의 분자 육종을 위한 표적이 될 수 있는 가치 있는 후보 유전자(예: Prupe.3G185000, Prupe.2G077800)와 대사물질(플라보노이드, AA2G)을 제공한다.

본 논문은 이러한 특정 경로와 유전자를 "핵심 요인" 및 "후보" 대상으로 식별하는 데 집중하며, 현재 데이터의 범위를 벗어난 검증되지 않은 응용이나 미래의 공학적 제안으로 과도하게 확장하지 않고 절제된 어조를 유지하고 있다.