An integrative single-cell and spatial transcriptomics atlas highlights candidate regulatory factors in the development of gerbera capitulum

기술 요약: Gerbera hybrida 두상화 발달의 통합적 단일 세포 및 공간 전사체 지도

문제 제기
두상화(capitulum, 꽃머리)를 특징으로 하는 Asteraceae(국과)는 가장 성공적인 진화적 데이타를 가진 속씨식물 계통 중 하나이다. 거베라(Gerbera hybrida) 두상화의 형태학적 복잡성은 잘 기록되어 있으나, 특히 화서 분화 과정 중 발생하는 전사체 수준의 조절 메커니즘은 여전히 이해가 부족하다. 몇몇 핵심 MADS-box 전사 인자들이 식별되었음에도 불구하고, 규제 지형은 불완전하며, 조밀한 화서 구조 내에서의 유전자 발현 공간 역학은 세포 해상도로 해결되지 않았다. 기존 연구들은 세포 유형 특이성과 공간적 맥락을 결结合하는 통합적 관점이 결여되어 있어, 두상화 패턴 형성 및 관다발 발달을 조절하는 유전자 조절 네트워크(GRN)를 분석하는 데 한계가 있다.

연구 방법론
저자들은 단일 세포 RNA 시퀀싱(scRNA-seq)과 고해 resolution 공간 전사체 기술(Stereo-seq)을 결합한 통합적 접근 방식을 사용하여 G. hybrida 두상화의 역동적인 지도를 구축하였다.

• 시료 수집: 화서 분화 단계를 포착하기 위해 이배체(double-haploid) 품종인 G. hybrida SH6로부터 세 단계의 연속적인 발달 단계(S1–S3)의 두상화를 선택하였다. 프로토플라스트(원형질체) 분리를 용이하게 하기 위해 포엽(involucral bracts)을 제거하였다.
• 단일 세포 시퀀싱: 10x Genomics 3' v3 플랫폼을 사용하여 세 단계에 대한 scRNA-seq 데이터를 생성하였다. 품질 관리 후 약 25,790개의 고품질 세포를 유지하였다. 데이터는 정규화, 고변이 유전자 탐지, PCA, 확산 매핑(diffusion mapping), Leiden 클러스터링을 포함하는 Scanpy를 사용하여 처리되었다.
• 공간 전사체 기술: S1–S3 범위를 포괄하는 19개의 두상화 샘플에 Stereo-seq을 적용하였다. 저품질 칩을 제외한 후 14개의 샘플을 분석하였다. 데이터의 희소성과 해상도 한계를 해결하기 위해 50×50 픽셀 도메인을 묶어 "의사 세포(pseudo-cells)"를 생성하였다.
• 데이터 통합: cell2location 알고리즘을 사용하여 단일 세포 참조 지도를 공간 전사체 데이터에 매핑함으로써, 세포 유형의 공간적 분포를 추론하였다.
• 다운스트림 분석:
  • 세포 유형 주석 달기: 마커 유전자 발현 및 차등 발현 유전자(DEG)의 유전자 기능 분석(GO enrichment)을 사용하여 클러스터를 주석 처리하였다.
  • 궤적 추론: moscot 패키지를 활용하여 단계 간의 세포 유형 전이 및 발달 궤적을 모델링하였다.
  • MADS-box 분석: 계통 발생 분석을 통해 G. hybrida 게놈에서 90개의 MADS-box 유전자를 식별하였다. 공간 상관관계 분석(convolution-based)과 계층적 클러스터링을 사용하여 이들의 발현 패턴을 분류하였다.
  • 조절 네트워크 추론: pySCENIC 파이프라인을 통해 scRNA-seq 데이터로부터 GRN을 추론하였다. 특정 유전자인 GAGL12의 경우, GRN 예측, CArG-box 모티프 스캐닝(3kb 상류 영역), 그리고 체관 특이적 DEG의 교집합을 통해 잠재적 타겟을 식별하였다.
  • 실험적 검증: GAGL12와 다른 MADS-domain 단백질 간의 단백질-단백질 상호작용을 검증하기 위해 효모 이중 잡종(Y2H) 분석을 수행하였다.

주요 결과

1. 세포 지도 및 역학: 본 연구는 G2/M 및 S기 세포, 표피, 화서/초기 화서 분생 조직(IM/early-FM), 화서 분생 조직(FM), 내피, 체관, 전재/물관, 광합성 유조직, 포엽 프리모디아를 포함한 10개의 주요 추정 세포 유형으로 집계된 33개의 클러스터를 식별하였다.

   • 발달적 전이: 세포 구성의 유의미한 변화를 확인하였다: IM/early-FM 세포(S1에서 풍부함)는 두상화가 성숙함에 따라 감소한 반의, FM 세포는 증가하였으며, 이는 분생 조직이 화기(floral organs)로 전환되는 것과 일치한다.
   • 체관 아형: 체관 세포는 세 가지 하위 집단(phloem_a, b, c)으로 구분되었다. RNA 속도(RNA velocity) 및 세포 전이 분석은 phloem_b와 phloem_c가 종말 상태를 향해 진행되는 반면, phloem_a는 안정적으로 유지되는 등 뚜렷한 발달 궤적을 보여주었다.

2. MADS-box 유전자의 공간적 발현 패턴: 저자들은 90개의 MADS-box 유전자를 네 가지 뚜렷한 공간적 발현 패턴으로 분류하였다:

   • 패턴 1 (두상화 전체): 광범위하게 발현되는 GSQUA 유전자(1–9)와 GRCD2/7을 포함하며, 이는 일반적인 두상화 및 화기 발달에서의 역할을 시사한다.
   • 패턴 2 (상단): 화기 및 기관에 농축되어 있으며, 기관 정체성(organ identity)에 관여하는 유전자들(예: GAGA1/2, GGLO1, GDEF2, GRCD1/4/5/8)을 포함한다.
   • 패턴 3 (중단): 조직 확장 및 관다발 지지와 관련이 있다. 이 패턴은 GAGL12와 GANR1L10이 체관 세포에 독특하게 농축되어 있음을 강조하였다.
   • 패턴 4 (하단): 광합성 및 내열성과 관련된 GSOC2 및 GSOC4(TM3-clade)만을 포함한다.

3. 관다발 조절자로서의 GAGL12 발견:

   • 공간적 특이성: GAGL12는 추정 체관 세포에 특이적으로 농축되어 있으며, 두상화 내의 추론된 체관 위치와 일치하는 공간적 분포를 보였다.
   • 단백질 상호작용: Y2H 분석은 GAGL12 단백질이 GRCD2/6/7, GSQUA1, GFLC6 및 MIKC*-유형 단백질(Ger_023266)을 포함한 여러 MADS-domain 단백질과 상호작용함을 확인하였다.
   • 타겟 예측: 통합 분석을 통해 플라스모데스마타(원형질 연락사) 기능 및 세포 간 수송에 관여하는 PDCB3(plasmodesmata callose-binding protein 3)와 PDLP6를 포함한 GAGL12의 잠재적 타겟 유전자를 식별하였다.

의의 및 주장
본 논문은 화기 분화 과정 중 G. hybrida 두상화의 최초의 통합적, 역동적 전사체 지도를 제공한다고 주장한다. scRNA-seq와 Stereo-seq을 결합함으로써, 본 연구는 벌크 시퀀싱이 달성할 수 없는 세포 이질성과 유전자 발현의 공간적 맥락을 해결하였다.

저자들은 이 연구가 복잡한 생식 구조의 조절 제어를 위한 기초를 마련한다고 강조한다. 특히, 체관 특이적 발현, 알려진 관다발 조절자와의 상호작용, 그리고 플라스모데스마타와 관련된 예측된 타겟들을 근거로 GAGL12를 두상화 관다발의 후보 조절자로 식별한 것은 Asteraceae의 관다발 발달에 대한 분자적 제어에 대한 새로운 가설을 제시한다. 또한, 본 연구는 특이적 공간 패턴을 가진 다수의 미규명 MADS-box 유전자를 강조하며, 이는 향후 기능적 특성 규명을 위한 로드맵을 제공한다. 저자들은 세포 유형 주석이 마커 유전자와 GO 농축에 의해 뒷받침되지만 여전히 "추정적(putative)"이며, 새로 식별된 유전자들의 정확한 생물학적 기능은 아직 밝혀져야 한다는 점을 언급하며 신중한 어조를 유지하였다.