Structural characterisation of stress-responsive miRNA precursors of rice (Oryza sativa L.).

기술 요약: 벼(Oryza sativa L.)의 스트레스 반응형 miRNA 전구체의 구조적 특성 규명

문제 제기
벼(Oryza sativa L.)는 다양한 생물적(예: 도열병, 세균성 벼 흰잎마름병, 잎집마름병) 및 비생물적(예: 가뭄, 염분, 중금속) 스트레스로 인해 상당한 생산성 손실을 겪고 있는 핵심적인 글로벌 주식 작물이다. 마이크로RNA(miRNA)는 유전자 발현을 조절함으로써 스트레스 반응을 조절하는 핵심적인 사후 전사 조절자로 확립되어 있으나, 벼의 전구체 분자(pre-miRNA)에 대한 구조적 다양성에 관한 지식은 여전히 불충분하다. miRNA의 생물학적 활성, 생성(biogenesis), 안정성 및 가공 효율은 전구체의 이차 구조 형태에 크게 의존한다. 그러나 스트레스 반응형 벼 miRNA의 조절 기능에 대한 스템(stem), 루프(loop), 벌지(bulge)와 같은 구조적 요소의 구체적인 기여도는 아직 충분히 규명되지 않았다. 본 연구는 약 150개의 스트레스 반응형 벼 pre-miRNA의 이차 구조적 특성에 대한 이해의 격차를 다룬다.

연구 방법
본 연구는 생물정보학적 접근 방식을 사용하여 스트레스 반응형 벼 pre-miRNA의 이차 구조를 분석하였다.

• 데이터 검색: miRBase(Release V22)에서 스트레스 반응과 관련된 약 150개의 전구체 서열(osa-mir)을 검색하였다.
• 구조 예측: 서열을 UNAFold 서버에 제출하여 이차 구조를 예측하였다. 파라미터에는 37°C에서의 최소 자유 에너지(MFE) 및 분할 함수 계산이 포함되었으며, 매달린 에너지(dangling energy) 고려 사항이 적용되었다.
• 큐레이션: 확립된 규칙(Zhang et al., 2005)에 따라 예측된 구조를 수동으로 큐레이션하였다.
• 통계 분석:
  • 클러스터링: 스템, 내부 루프, 벌지, 터미널 루프의 개수를 기준으로 샘플을 그룹화하기 위해 R Studio를 사용하여 이원 계층적 클러스터링 히트맵을 생성하였다.
  • 주성분 분석(PCA): 차원을 축소하고 구조적 구성 요소(스템, 벌지, 루프) 간의 분산을 식별하기 위해 PCA를 수행하여 이들이 구조적 안정성에 미치는 영향을 결정하였다.
  • 열역학 분석: 열역학적 안정성을 평가하기 위해 자유 에너지($\Delta G$) 값을 계산하였다.

주요 결과

• 구조적 패턴: 스템, 내부 루프, 벌지, 터미널 루프로 구성된 네 가지 뚜렷한 구조적 패턴이 확인되었다.
  • 스템: 스템의 개수는 1개에서 6개 사이였다. 약 80%의 전구체가 1~2개의 스템을 포함하고 있었으며, osa-miR399d가 최대 6개를 나타냈다.
  • 내부 루프: 개수는 0개에서 6개 사이였다. 특히, 34%의 전구체가 내부 루프가 완전히 결여되어 있었으며, osa-miR408과 osa-miR167e가 각각 6개의 루프를 가져 가장 높은 복잡성을 보였다.
  • 벌지: 개수는 0개에서 7개 사이였다. 약 23%의 전구체는 벌지가 없었으며, osa-miR156b가 최대 7개를 나타냈다.
  • 터미널 루프: 대부분(85%)의 전구체는 정형적인 헤어핀 구조와 일치하게 1~2개의 터미널 루프를 보유하고 있었다.
• 열역학적 안정성: $\Delta G$ 값이 –28.70 kcal/mol에서 –342.00 kcal/mol까지 분포함에 따라 구조적 안정성에서 상당한 차이가 관찰되었다.
  • osa-miR7695가 가장 안정적인 구조($\Delta G = -342.00$ kcal/mol)로 확인되었다.
  • osa-miR395(특히 osa-miR395u)는 가장 불안정한 구조($\Delta G = -28.70$ kcal/mol)였다.
  • 이는 분석된 패밀리들 사이에서 약 12배의 안정성 차이를 나타낸다.
• 다변량 분석:
  • PCA: 첫 두 주성분이 전체 분산의 100%를 설명하였다(PC1: 65.2%, PC2: 34.8%). PC1은 벌지와 스템 루프를 분리하였고, PC2는 다른 구조적 차이를 포착하였다. 이 분석은 구조적 그룹 간의 뚜렷한 분리를 보여주었으며, 이는 벌지, 스템, 루프의 고유한 특성을 시사한다.
  • 히트맵 클러스터링: 비지도 클러스터링 결과, 스템과 루프 영역이 벌지보다 서로 더 밀접하게 연관되어 있음을 보여주었다. 특정 miRNA 패밀리(예: osa-miR156, osa-miR166, osa-miR169, osa-miR395, osa-miR399, osa-miR171)는 뚜렷한 클러스터를 형성하였으며, 이는 공유된 구조적 선호도를 나타낸다.

의의 및 주장
저자들은 이차 구조 클러스터링과 PCA 결과 사이의 강력한 일치성은 pre-miRNA의 이차 구조 변화가 그 생성 및 조절 기능에 직접적인 영향을 미칠 수 있음을 시사한다고 주장한다.

• 구조가 기능에 미치는 영향: 연구는 내부 루프 및 벌지와 같은 구조적 요소가 Dicer 가공 효율 및 표적 인식에 결정적이라고 상정한다. 예를 들어, 내부 루프의 존재는 성숙한 miRNA의 길이에 영향을 줄 수 있으며, 벌지는 절단의 정밀도에 영향을 줄 수 있다.
• 안정성 및 반감기: 자유 에너지 값의 넓은 범위는 구조적 안정성의 상당한 다양성을 나타낸다. 저자들은 낮은 $\Delta G$를 갖는 매우 안정적인 구조가 뉴클레아제 분해에 대한 저항성이 높고 세포 내 반감기가 길다는 점을 언급하였는데, 이는 헤어핀 miRNA에 대한 기존 연구와 일치하는 결과이다.
• 응용: 이 결과는 스트레스 내성을 강화하기 위한 miRNA 기반 전략을 설계하는 데 유용한 구조적 통찰력을 제공한다. 저자들은 이러한 구조-기능 상관관계를 이해하는 것이 비생물적 및 생물적 스트레스에 대한 회복력이 향상된 벼 품종을 개발하기 위한 필수적인 단계라고 제언한다.

본 논문은 구조적 다양성이 유의미하지만, osa-miRNA의 구체적인 구조-기능 상관관계를 완전히 규명하고 새로운 스트레스 반응형 전구체를 식별하기 위해서는 향후 연구가 필요하다고 결론짓는다.