Candidate Homologs of REST-Associated Chromatin Regulatory Components in Choanoflagellates

기술 요약: 깃털세포류(Choanoflagellates) 내 REST 관련 염색질 조절 구성 요소의 후보 상동성

문제 제기
REST/NRSF는 C2H2 아연 집게(zinc finger)를 통해 RE1/NRSE DNA 모티프에 결합하고 코어레스(co-repressor) 복합체(CoREST/RCOR, SIN3, HDACs, LSD1 포함)를 모집하여 신경 세포 유전자 프로그램을 억제하는, 잘 규명된 동물 전사 억제 인자이다. 이 조절 논리의 진화적 기원을 이해하기 위해서는 이 경로의 구성 요소 중 어떤 것이 동물 계통 이전에 존재했는지 결정하는 것이 필수적이다. 동물의 가장 가까운 친척인 깃털세포류는 이 탐구를 위한 핵심 외집단(outgroup) 역할을 한다. 본 연구는 두 가지 별개의 질문을 다룬다: (1) 깃털세포류가 REST/NRSF의 직접적인 직계 상동 유전자(ortholog)를 암호화하는가? (2) 그들은 동물의 REST 관련 억제 기구(DNA 결합 단백질, 공억제 인자, 염색질 변형 효소 등 더 넓은 단백질 패밀리)의 후보 상동체들을 암호화하는가?

방법론
본 연구는 세 개의 완전한 게놈 프로테옴(Monosiga brevicollis, Salpingoeca rosetta, Choanoeca flexa)과 20개의 전사체 유래 데이터 세트를 포함하여 총 702,196개의 예측 단백질로 구성된 23개의 깃털세포류 단백질 데이터 세트에 대해 계층적 계산 조사를 수행하였다.

1. HMM 기반 스크리닝: REST 기능과 관련된 도메인을 탐지하기 위해 Pfam/InterPro에서 큐레이션된 HMM 패널을 구축하였다. 여기에는 다음이 포함된다:

   • DNA 결합: C2H2 아연 집게, bHLH, homeobox, POU, Myb/SANT, forkhead, MADS, GATA, HMG box, T-box.
   • 억제 구조: BTB/POZ, KRAB, SCAN.
   • 공억제 복합체: ELM2, SIN3 PAH, 히스톤 탈아세틸화 효소(HDAC), LSD1 유형 아민 산화효소, SWIRM.
   • 염색질 리더/변형 효소: PHD, bromodomain, chromodomain, BAH, SET, SIR2.
     패널은 gathering threshold를 적용하여 모든 단백질 세트에 대해 HMMER로 검색되었다.

2. 후보 우선순위 지정: 5,123개의 후보 단백질이 식별되었다. 이들은 특정 도메인(예: C2H2, BTB/POZ, corepressor 도메인)의 존재 여부와 최적의 Swiss-Prot 히트 주석(예: "REST", "HDAC", "transcriptional repressor")에 가중치를 두는 점수 공식에 따라 순위가 매겨졌다. 후보들은 Low, Moderate, High, Priority 티어로 계층화되었다.

3. 서열 검증: DIAMOND 검색(E-value < 1e-5)과 상위 50개 우선순위 후보에 대한 2차 BLASTP 확인을 통해 후보들을 검증하였다. HH-suite는 데이터베이스 제약으로 인해 사용되지 않았으나 향로 확장 가능성으로 언급되었다.

4. 구조 준비: 40개의 우선순위 후보에 대해 AlphaFold 준비용 FASTA 파일을 생성하였다. 향후 구조 비교를 용이하게 하기 위해 최적의 Swiss-Prot 상동체로부터 참조 구조를 AlphaFold DB에서 검색하였다(12개의 고유한 UniProt 참조 구조 확보).

5. 모티프 분석: 세 개의 완전한 게놈 어셈블리에 대해 정형적인 RE1/NRSE 모티프(TTCAGCACCACGGACAGCGCC)를 대상으로 프로모터 스캔(유전자 상류 최대 1 kb까지)을 수행하였으며, 최대 5개의 미스매치를 허용하였다. 농축도를 평가하기 위해 분석적 GC-매칭 null 모델을 계산하였고, Poisson tail 확률과 Benjamini-Hochberg FDR 값을 산출하였다.

주요 결과

• 도메인 출현 빈도: 조사는 5,123개의 후보에서 8,299개의 도메인 히트를 식별하였다. C2H2 아연 집게, BT-POZ, HDAC, SIN3 PAH 및 다양한 염색질 리더와 같은 패밀리가 23개 샘플 모두에서 검출되었으며, 이는 깃털세포류에 광범위한 조절 툴킷이 존재함을 나타낸다.
• REST/NRSF 직계 상동 유전자의 부재: 어떤 후보 단백질도 설득력 있는 최적 히트 주석(NRSF/REST)을 반환하지 않았다. 조사된 데이터 세트에서 직접적인 REST/NRSF 직계 상동 유전자의 증거는 발견되지 않았다.
• 공억제 후보 존재: 가장 강력한 신호는 REST 자체보다는 REST에 의해 모집되는 *기구(machinery)*의 후보 상동체를 식별한 데서 나타났다.
  • HDAC/SIN3: 히스톤 탈아세틸화 효소 유사 및 SIN3 유사 PAH 도메인 단백질이 상위 우선순위 티어를 지배하였다.
  • RCOR/CoREST: 5개의 후보가 ELM2 및/또는 Myb-SANT 도메인을 포함하며, 최적의 Swiss-Prot 히트는 RCOR1/2/3 또는 "REST corepressor" 단백질이었다. ELM2/SANT 도메인은 널리 퍼져 있으며 상동성 확립을 위해 추가적인 계통 발생학적 검증이 필요하므로, 이들은 잠정적 상동체로 취급된다.
  • YY1 유사 억제 인자: 10개의 C2H2 아연 집게 후보가 인간 YY1 유사 전사 억제자와 일치하였으며, 이는 REST와 구별되는 서열 특이적 억제제의 존재를 시사한다.
• 모티프 농축: 717개의 느슨하게 일치하는 RE1 유사 서열이 발견되었다. GC-매칭 null 모델은 완화된 임계값(4 및 5 미스매치)에서 S. rosetta와 C. flexa의 농축을 나타냈으나, 가장 엄격한 임계값(3 미스매치)에서는 미미한 증거(예: S. rosetta에서 2개 히트, FDR = 0.057)만을 보여, 이들이 확정된 조절 요소라기보다 후보 사이트임을 시사한다.
• 구조 데이터: 향후 예측된 깃털세포류 후보 구조와의 비교를 위해 12개의 고유한 참조 구조가 준비되었다.

의의 및 주장
본 논문은 보수적인 해석으로 결론을 맺는다: 깃털세포류는 동물 REST 관련 억제에 참여하는 여러 단백질 패밀리(특히 HDAC, SIN3 유사 단백질, 잠정적인 RCOR/CoREST 유사 단백질, C2H2 억제 인자)의 후보 상동체를 암호화하고 있다. 그러나 본 연구는 이 구성 요소들이 깃털세포류에서 통합된 REST 유사 경로로서 함께 기능한다는 것을 명시적으로 주장하지 않는다.

주된 의의는 개별 구성 요소의 존재와 기능적 경로의 존재를 구분하는 데 있다. REST 관련 억제를 위한 "부품 목록(parts list)"은 대부분 선(pre)-동물 단계에 존재하는 것으로 보이지만, 본 연구는 깃털세포류에서 직접적인 REST/NRSF 직계 상동 유전자나 확정된 통합 복합체의 증거를 찾지 못했다. 잠정적인 RCOR/CoREST 유사 후보의 식별은 동물 REST 복합체와의 가장 직접적인 연결 고리로 강조되며, 이는 향후 계통 발생 및 구조적 조사를 요한다. 모티프 결과는 기능적 결합의 증거라기보다 향후 조절 분석을 위한 탐색적 단서로 제시된다.

저자가 언급한 한계점

• 서열 유사성이 생화학적 기능을 의미하지는 않는다.
• 구조 비교가 준비되었으나 실행되지는 않았다 (로컬 AlphaFold 예측은 생성되지 않음).
• 모티프 분석은 단순한 GC-매칭 null 모델을 사용하였으며, 계통 발생적 보존성이나 인터제닉(intergenic) 배경에 대한 통제가 부족하다.
• 전사체 유래 세트는 복제 수 해석을 복잡하게 만드는 아티팩트나 이소형(isoform)을 포함할 수 있다.
• RCOR/CoREST 상동성은 상호 유사성 및 계통 발생 재구성을 통한 검증이 필요한 잠정적인 상태이다.