Comprehensive transcriptomic analysis of BjGL1-knockout Brassica juncea: novel insights into leaf trichome formation and defense-associated biological processes

*기술 요약: BjGL1 결실 Brassica juncea의 포괄적 전사체 분석*

문제 제lemnt (Problem Statement)
Brassica juncea(갓)는 기름, 향신료 및 의약품을 위한 중요한 작물이지만, 생물적 스트레스, 특히 진딧물 침해로 인해 상당한 어려움을 겪고 있다. 잎의 모상체(trichome)는 이러한 해충에 대한 일차적인 물리적 및 화학적 장벽 역할을 한다. 이전 연구를 통해 두 개의 상동 유전자인 BjA06.GL1과 BjB02.GL1을 동시에 결실시키면 잎의 모상체가 완전히 소실되고 진딧물 저항성이 감소한다는 사실이 밝혀졌으나, 이 유전자들이 영향을 미치는 하위 분자 조절 네트워크와 광범위한 생물학적 과정은 여전히 명확히 이해되지 않고 있다. 구체적으로, 유전자 교정된 겨자에서 이들 유전자의 부재가 모상체 형성, 이차 대사산물 생합성 및 방어 반응과 관련된 유전자 발현을 어떻게 변화시키는지 규명하는 전사체 데이터가 부족한 실정이다.

방법론 (Methodology)
본 연구는 CRISPR/Cas9 유전자 편집과 RNA 시퀀싱(RNA-seq) 및 정량적 실시간 PCR(qRT-PCR) 검증을 결합한 멀티 오믹스 접근법을 사용하였다:

1. 식물 재료 및 유전자 편집: 연구진은 CRISPR/Cas9을 사용하여 B. juncea 내의 BjA06.GL1 및 BjB02.GL1의 첫 번째와 세 번째 엑손을 동시에 표적하였다. 모상체 결핍 표현형(BG로 명명)을 보이는 이중 결실 돌연변이 T₁ 세대 식물을 정상적인 모상체 밀도를 가진 야생형 대조군(BT)과 비교하였다.
2. 유전자형 및 오프 타겟(Off-target) 검증: 돌연변이 확인을 위해 DNA 추출 후 PCR 증폭 및 생거 시퀀싱(Sanger sequencing)을 수행하였다. 오프 타겟 효과는 B. juncea 참조 게놈에 대한 BLAST 정렬을 통해 엄격하게 스크리닝하였으며, 이후 잠재적 오프 타겟 로커스에 대한 PCR 및 시퀀싱을 실시하여 높은 편집 특이성을 확인하였고 감지된 오프 타겟 돌연변이는 없었다.
3. 전사체 프로파일링: BT와 BG 라인(각각 3개의 생물학적 반복구)의 21일 된 잎에서 RNA를 추출하였다. 라이브러리는 Illumina NovaSeq 6000 플랫폼에서 시퀀싱되었다. 클린 리드(Clean reads)는 HIS2를 사용하여 B. juncea (Xuecai) 참조 게놈에 정렬되었으며, StringTie을 통해 전사체 조립이 수행되었다.
4. 차등 발현 분석: DEGs(차등 발현 유전자)는 |log₂ fold change| ≥ 1 및 조정된 P-value < 0.05의 임계값을 사용하여 DESeq2로 식별되었다. 유전자 발현은 FPKM으로 정량화되었다.
5. 기능적 주석 및 검증: DEGs는 GO(Gene Ontology) 및 KEGG 경로 농축 분석의 대상이 되었다. 또한, 발현 패턴에 따라 유전자를 그룹화하기 위해 K-means 클러스터링이 사용되었다. 전사 인자 및 대사 경로 유전자를 포함한 주요 DEGs는 qRT-PCR을 통해 검증되었다.

주요 기여 및 결과 (Key Contributions and Results)

• 표현형 확인: 이중 결실 BG 라인은 주면(adaxial)과 밑면(abaxial) 모두에서 잎의 모상체가 완전히 결여된 모습을 보였으며, 이는 주사전자현미경(SEM)을 통해 확인되었다. 돌연변이는 표적 부위에서의 결실 및 삽입에 의한 프레임시프트(frameshift)를 포함하며, T₁ 세대에서 안정적으로 유전되었다.
• 전사체 경관: 본 연구는 WT(BT)와 결실체(BG) 사이에서 4,604개의 차등 발현 유전자(DEGs)를 식별하였다. 구체적으로, BG 돌연변이 대비 모상체가 있는 BT 잎에서 1,831개의 유전자가 상향 조절되었고 2,773개의 유전자가 하향 조절되었다.
• 농축된 경로: 기능적 농축 분석 결과, DEGs는 다음과 관련하여 유의미한 연관성을 보였다:
  • 사이토크롬 P450 (CYP) 대사: 내인성 전구체의 산화적 변환 및 방어용 이차 대사산물의 생합성에 결정적이다.
  • 식물-병원균 상호작용: 모상체의 존재와 면역 반응 메커니즘 사이의 직접적인 연결을 시사한다.
  • 이차 대사산물 생합성: 식물의 방어 및 스트레스 적응에 필수적인 플라보노이드 및 페닐프로파노이드 경로를 포함한다.
  • 환경 적응 및 막 수송.
• 전사 인자 네트워크: 분석 결과 357개의 차등 발현 전사 인자(TFs)가 식별되었으며, ERF, MYB, NAC 가문이 가장 풍부했다.
  • 양성 조절 인자: 핵심 유전자인 BjGL1, BjGL2, BjZFP8, BjSPL9, BjMYB106은 BT에서 유의미하게 상향 조절되었다.
  • 음성 조절 인자: 모상체 발달을 억제하는 것으로 알려진 BjTCP10, BjSPL5, 그리고 R3형 MYB 유전자(BjTRY, BjCPC, BjETC1)는 뚜렷한 발현 패턴을 보였으며, 억제제들은 일반적으로 BT의 활발한 조절 네트워크 맥락에서 하향 조절되거나 덜 활성화되어 있었다.
  • 주목할 점은, B. juncea 게놈에 4개의 GL1 상동 유전자가 존재하지만, BjA06.GL1과 BjB02.GL1만이 유의미한 차등 발현을 보였으며, 나머지 두 상동 유전자(BjuLA09G063230 및 BjuLB07G024530)는 차등 발현되지 않았다.
• 검증: qRT-PCR 결과는 RNA-seq 데이터에서 식별된 주요 조절 유전자의 발현 경향을 확인하여, 전사체 연구 결과의 신뢰성을 검증하였다.

의의 및 주장 (Significance and Claims)
저자들은 본 연구가 Brassica juncea의 모상체 매개 저항성을 지배하는 분자적 메커니즘을 이해하기 위한 중요한 이론적 토대를 제공한다고 주장한다. 본 연구의 주요 의의는 BjA06.GL1과 BjB02.GL1이 단순한 구조적 결정 요인이 아니라 광범위한 전사 네트워크를 조정하는 중심 조절자임을 입증했다는 데 있다. 이 네트워크는 단순한 형태 형성(morphogenesis)을 넘어 다음을 포함한다:

1. 방어 메커니즘: 모상체가 풍부한 식물에서 식물-병원균 상호작용 유전자와 이차 대사산물 경로가 상향 조절된다는 것은 모상체 발달이 식물의 화학적 방어 체계와 본질적으로 연결되어 있음을 시사한다.
2. 대사 조절: 본 연구 결과는 이 유전자들이 페닐프로파노이드 경로와 사이토크롬 P450 대사에 영향을 미쳐, 생물적 및 비생물적 스트레스에 대한 저항성을 높이기 위해 식물의 화학적 프로필을 변화시킬 수 있음을 나타낸다.
3. 육종 응용: 특정 조절 유전자와 경로를 밝힘으로써, 본 연구는 향에 강화된 곤충 저항성과 화학 농약 의존도 감소를 목표로 하는 B. juncea 품종 개발을 위한 분자 육종 프로그램의 가치 있는 표적을 제공한다.

본 논문은 본 연구가 BjGL1 결실과 관련된 전사체 변화를 성공적으로 지도하였으나, 하위 대사 변화와 이들 조절 네트워크의 정확한 생태적 역할을 완전히 규명하기 위해서는 대사체학(metabolomics)을 통합한 추가 연구가 필요하다고 결론짓는다.