Feasibility of RNA Sequencing from Human Skin Microbiopsies in Actinic Keratosis and Basal Cell Carcinoma: A Paired Pilot Comparison with Shave Biopsy

기술 요약: 광선 각화증 및 기저 세포 암종에서의 인간 피부 미세 생검을 통한 RNA 시퀀싱의 타당성

문제 제기
전통적인 쉐이브 생검(shave biopsy)은 피부과적 진단 및 분자 연구의 표준이지만 침습성 측면에서 한계가 있다. 이러한 침습성은 미용적으로 민감한 부위에서의 샘플링을 제한하고, 반복적인 종단적 평가를 방해하며, 탐색적 연구를 위한 피험자 모집을 저해할 수 있다. 서브 밀리미터(sub-millimeter) 단위의 피부 미세 생검은 분자 샘핑을 위한 최소 침습적이고 봉합이 필요 없는 대안을 제공하지만, 케라티노사이트 유래 병변에 대해 RNA 시퀀싱(RNA-seq) 용도로써 기존 생검과 비교하여 그 효능을 입증하는 쌍을 이룬 근거는 부족한 실정이다.

방법론
본 연구는 두 개의 인간 병변(하나의 광선 각화증(AK, 샘플 RH001)과 하나의 기저 세포 암종(BCC, 샘플 RW002))을 대상으로 탐색적 쌍 비교 파일럿 연구를 수행하였다. 각 병변은 전통적인 쉐브 생검과 미세 생검을 모두 사용하여 샘플링되었다. 워크플로우는 다음과 같다:

• RNA 처리: RNA 무결성(integrity) 및 농도 평가 후, 가닥 구분 RNA-seq(stranded RNA-seq) 라이브러리 제작을 진행하였다.
• 시퀀싱: 라이브러리는 Illumina HiSeq 2000 플랫폼을 사용하여 100-bp 페어드 엔드(paired-end) 리드로 시퀀싱되었다.
• 생물정보학: 품질 관리는 FastQC 및 Trim Galore을 통해 수행되었다. 리드는 GENCODE v19에 의해 가이드된 TopHat v2.1을 사용하여 인간 게놈(hg19)에 정렬되었다. 유전자 수준의 카운트는 HTSeq로 생성되었으며 edgeR을 사용하여 정규화되었다.
• 분석: 샘플링 방법 간의 일치도는 log10 변환된 FPKM 산점도를 통해 평가되었다. 마커 발현은 소스 마커 세트와 Human Protein Atlas의 피부 층 주석을 사용하여 해석되었다.

주요 결과

• RNA 품질 및 수율: 조직 부피가 작음에 따라 예상된 바와 같이, 미세 생검 샘플은 현저히 낮은 초기 RNA 농도를 나타냈다 (예: AK 미세 생검 2.85 ng/µL vs. AK 쉐브 생검 59 ng/µL). 그러나 RNA 무결성 지수(RIN)는 라이브러리 제작에 적합한 수준을 유지했다 (AK 미세 생검 RIN 5.9 vs. 쉐브 6.7; BCC 미세 생검 RIN 7.1 vs. 쉐브 7.7).
• 시퀀싱 성능: 낮은 입력량에도 불구하고, 미세 생검 라이브러리는 높은 매핑률을 달ien했다. AK 쌍의 경우, 미세 생검의 매핑된 리드 비율은 94.8%였으며 쉐브 생검은 92.5%였다. BCC 쌍의 경우, 미세 생검은 91.5%, 쉐브 생검은 94.7%를 기록했다.
• 전사체 일치도: 전역적 발현 프로파일은 두 샘플링 방법 간에 높은 일치도를 보였으며, log10 FPKM 산점도의 $r^2$ 값은 두 병변 모두 0.93이었다.
• 생물학적 해석 가능성:
  • 광선 각화증: 미세 생검 샘플은 표피 상층부 분화 마커(예: FLG, FLG2, CDSN, KRT1, KRT10, CASP14)의 농축을 보였다. 예를 들어, FLG 발현은 미세 생검에서 713.61 FPKM이었던 반면 쉐브 생검에서는 76.90 FPKM이었다. 반대로, 콜라겐 전사체(COL1A1, COL3A1)는 미세 생검에서 더 낮게 나타났는데, 이는 얕은 샘플링 및 기질(stroma) 포집 감소와 일치한다.
  • 기저 세포 암종: 미세 생검은 KRT17(미세 생검 3822.90 FPKM vs. 쉐브 생검 5090.87 FPKM)을 포함한 종양 관련 마커의 발현을 유지했다.

의의 및 주장
저자들은 인간 피부 미세 생검이 시퀀싱 가능한 RNA-seq 라이브러리를 생성하고 AK 및 BCC로부터 해석 가능한 전사체 프로파일을 회복하는 데 타당하다고 결론짓는다. 데이터는 낮은 RNA 입력량에도 불구하고, 이 방법이 기존 생검과 높은 일치도를 가진 고품질 데이터를 생성함을 입증한다.

그러나 본 논문은 이러한 발견을 진단적 동등성이나 임상적 유용성에 대한 증거가 아닌 **타당성 데이터(feasibility data)**로 명시적으로 규정한다. 저자들은 다음과 같은 몇 가지 제한 사항을 강조한다:

• 본 연구는 매우 작은 샘플 크기(병변 유형당 한 쌍)를 가진 파일럿 연구이며, 진단 내에서의 생물학적 복제가 결여되어 있다.
• 일치하는 조직학적 깊이 평가가 이루어지지 않았다.
• 소스 자료에 불완전한 방법론 메타데이터가 존재한다.

따라서 저자들은 분산(dispersion) 및 볼케이노 플롯(volcano plots)을 질병 수준의 차별 발현 증거로 해석해서는 안 된다고 명시한다. 저자들은 질병 수준의 결론을 도출하기 전에 표준화된 샘플링 깊이, 조직학적 상관관계, 그리고 재현 가능한 데이터 배포를 포함한 더 큰 규모의 쌍 비교 연구가 필요하다고 주장한다.