Name: 작은 데이터 세트에서 예측 작업 및 다크 바이오마커 검출을 위한 전사체 특징의 전사체 조절 보기 생성 (동영상) | JoVE Journal
Uploaded: 2024-03-01T13:50:00.000+00:00
Duration: 3 min 37 s
Description: 학습하세요 mqTrans 전사체 조절 보기를 활용한 다크 바이오마커 검출 방법을 Python 가상 환경에서 구현하는 방법을. 이 연구는 소규모 데이터 세트에서 221개의 다크 바이오마커를 성공적으로 식별하여 기존 방법의 한계를 극복하고 정밀 진단을 위한 새로운 바이오마커 스크리닝 기법을 제시합니다.

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

이 작업은 시니어 및 주니어 기술 혁신팀(20210509055RQ), 구이저우성 과학기술프로젝트(ZK2023-297), 구이저우성 위생건강위원회 과학기술재단(gzwkj2023-565), 길림성 교육부 과학기술프로젝트(JJKH20220245KJ, JJKH20220226SK), 중국 국가자연과학재단(U19A2061), 길림성 빅데이터 지능컴퓨팅 중점연구소의 지원을 받았다 (20180622002JC) 및 중앙 대학을 위한 기초 연구 기금, JLU. 우리는 이 프로토콜의 엄격함과 명확성을 실질적으로 향상시키는 데 중요한 역할을 한 건설적인 비평에 대해 리뷰 편집자와 세 명의 익명 리뷰어에게 진심으로 감사를 표합니다.

참고: 다음 프로토콜은 주요 모듈의 정보학 분석 절차 및 Python 명령에 대한 세부 정보를 설명합니다. 그림 2는 이 프로토콜에서 사용되는 예제 명령과 함께 세 가지 주요 단계를 설명하며 자세한 기술 세부 사항은 이전에 발표된 저작물26,38을 참조하십시오. 컴퓨터 시스템의 일반 사용자 계정으로 다음 프로토콜을 수행하고 관리자 또는 루트 계정을 사용하지 마십시오. 이것은 계산 프로토콜이며 생물 의학적 위험 요소가 없습니다.
1. Python 환경 준비
<ol><li>가상 환경을 만듭니다.<ol><li>이 연구는 Python 3.7에서 Python 프로그래밍 언어와 Python 가상 환경(VE)을 사용했습니다. 다음 단계를 수행합니다(그림 3A). conda create -n healthmodel 파이썬=3.7 conda create 는 새 VE를 만드는 명령입니다. 매개 변수 -n은 새 환경의 이름(이 경우 healthmodel)을 지정합니다. 그리고 python=3.7 은 설치할 Python 버전을 지정합니다. 위의 명령을 지원하는 원하는 이름과 Python 버전을 선택합니다.</li>
		<li>명령을 실행한 후 출력은 그림 3B와 유사합니다. y 를 입력하고 프로세스가 완료될 때까지 기다립니다.</li>
	</ol></li>
	<li>가상 환경 활성화<ol><li>대부분의 경우 다음 명령을 사용하여 생성된 VE를 활성화합니다(그림 3C). conda 활성화 healthmodel</li>
		<li>일부 플랫폼에서 사용자가 활성화를 위해 플랫폼별 구성 파일을 업로드해야 하는 경우 VE 활성화에 대한 플랫폼별 지침을 따릅니다.</li>
	</ol></li>
	<li>PyTorch 1.13.1 설치<ol><li>PyTorch는 인공 지능(AI) 알고리즘을 위한 인기 있는 Python 패키지입니다. CUDA 11.7 GPU 프로그래밍 플랫폼을 기반으로 하는 PyTorch 1.13.1을 예로 사용합니다. https://pytorch.org/get-started/previous-versions/ 에서 다른 버전을 찾으십시오. 다음 명령을 사용합니다(그림 3D). pip3 토치 설치 torchvision torchaudio 참고: PyTorch 버전 1.12 이상을 사용하는 것이 좋습니다. 그렇지 않으면 공식 torch_geometric 웹 사이트 https://pytorch-geometric.readthedocs.io/en/latest/install/installation.html 에 명시된 것처럼 필요한 패키지를 torch_geometric 설치하는 것이 어려울 수 있습니다.</li>
	</ol></li>
	<li>torch-geometric을 위한 추가 패키지 설치
	<ol><li>https://pytorch-geometric.readthedocs.io/en/latest/install/installation.html 의 지침에 따라 다음 명령을 사용하여 torch_scatter, torch_sparse, torch_cluster 및 torch_spline_conv 패키지를 설치합니다(그림 3E). pip 설치 pyg_lib torch_scatter torch_sparse torch_cluster torch_spline_conv -f https://data.pyg.org/whl/torch-1.13.0+cu117.html</li>
	</ol></li>
	<li>torch-geometric 패키지를 설치합니다.<ol><li>이 연구에는 torch-geometric 패키지의 특정 버전 2.2.0이 필요합니다. 명령을 실행합니다(그림 3F). pip 설치 torch_geometric==2.2.0</li>
	</ol></li>
	<li>다른 패키지를 설치합니다.<ol><li>pandas와 같은 패키지는 일반적으로 기본적으로 사용할 수 있습니다. 그렇지 않은 경우 pip 명령을 사용하여 설치합니다. 예를 들어 pandas 및 xgboost를 설치하려면 다음을 실행합니다. 핍 설치 팬더 pip 설치 xgboost</li>
	</ol></li>
</ol>2. 사전 훈련된 HealthModel을 사용하여 mqTrans 기능 생성
<ol><li>코드 및 미리 학습된 모델을 다운로드합니다.<ol><li>http://www.healthinformaticslab.org/supp/resources.php(HealthModel-mqTrans-v1-00.tar.gz)라는 웹 사이트에서 코드와 사전 훈련된 HealthModel을 다운로드합니다(그림 4A). 다운로드한 파일은 사용자가 지정한 경로로 압축을 풀 수 있습니다. 구현된 프로토콜의 자세한 공식화 및 지원 데이터는26에서 찾을 수 있습니다.</li>
	</ol></li>
	<li>HealthModel을 실행하기 위한 파라미터를 소개합니다.<ol><li>먼저 명령행에서 작업 디렉토리를 HealthModel-mqTrans 폴더로 변경하십시오. 다음 구문을 사용하여 코드를 실행합니다. 파이썬 main.py <데이터 폴더> <모델 폴더> <출력 폴더> 각 매개 변수와 데이터, 모델 및 출력 폴더에 대한 세부 정보는 다음과 같습니다. 데이터 폴더: 원본 데이터 폴더이며 각 데이터 파일은 csv 형식입니다. 이 데이터 폴더에는 두 개의 파일이 있습니다(2.3단계 및 2.4단계의 자세한 설명 참조). 이러한 파일은 개인 데이터로 바꿔야 합니다. data.csv: 전사체 매트릭스 파일입니다. 첫 번째 행에는 특징(또는 유전자) ID가 나열되고 첫 번째 열에는 샘플 ID가 표시됩니다. 유전자 목록에는 조절 인자(TF 및 lincRNA)와 조절된 mRNA 유전자가 포함됩니다. label.csv: 샘플 레이블 파일입니다. 첫 번째 열에는 샘플 ID가 나열되고 이름이 "label"인 열에는 샘플 레이블이 제공됩니다. model folder: 모델에 대한 정보를 저장할 폴더입니다. HealthModel.pth: 미리 학습된 HealthModel입니다. regulatory_geneIDs.csv: 이 연구에 사용된 조절 유전자 ID입니다. target_geneIDs.csv: 본 연구에서 사용된 표적 유전자. adjacent_matrix.csv: 조절 유전자의 인접 매트릭스. output folder: 출력 파일은 코드에 의해 생성된 이 폴더에 기록됩니다. test_target.csv: Z-Normalization 및 imputation 후 표적 유전자의 유전자 발현 값. pred_target.csv: 표적 유전자의 예측된 유전자 발현 값. mq_target.csv: 표적 유전자의 예측된 유전자 발현 값.</li>
	</ol></li>
	<li>csv 형식의 전사체 매트릭스 파일을 준비합니다.<ol><li>각 행은 샘플을 나타내고 각 열은 유전자를 나타냅니다(그림 4B). 트랜스크립토믹 데이터 매트릭스 파일의 이름을 데이터 폴더에 data.csv으로 지정합니다. 알림: 이 파일은 Microsoft Excel과 같은 소프트웨어에서 데이터 매트릭스를 .csv 형식으로 수동으로 저장하여 생성할 수 있습니다. 전사체 매트릭스는 또한 컴퓨터 프로그래밍에 의해 생성될 수 있다.</li>
	</ol></li>
	<li>csv 형식의 레이블 파일을 준비합니다.<ol><li>전사체 매트릭스 파일과 마찬가지로 레이블 파일의 이름을 데이터 폴더의 label.csv로 지정합니다(그림 4C). 참고: 첫 번째 열은 샘플 이름을 제공하고 각 샘플의 클래스 레이블은 레이블이라는 열에 제공됩니다. 레이블 열의 0 값은 이 샘플이 음수임을 의미하고 1은 양성 샘플을 의미합니다.</li>
	</ol></li>
	<li>mqTrans 기능을 생성하십시오.<ol><li>다음 명령을 실행하여 mqTrans 기능을 생성하고 그림 4D에 표시된 출력을 가져오십시오. mqTrans 기능은 ./output/mq_targets.csv 파일로 생성되고 레이블 파일은 ./output/label.csv 파일로 다시 저장됩니다. 추가 분석의 편의를 위해 mRNA 유전자의 원래 발현 값도 ./output/ test_target.csv 파일로 추출됩니다. 파이썬 ./Get_mqTrans/code/main.py ./data ./Get_mqTrans/model ./output</li>
	</ol></li>
</ol>3. mqTrans 기능 선택
<ol><li>기능 선택 코드의 구문<ol><li>먼저 작업 디렉토리를 HealthModel-mqTrans 폴더로 변경하십시오. 다음 구문을 사용합니다. 파이썬 ./FS_classification/testMain.py <in-data-file.csv> <in-label-file.csv> <output folder> <select_feature_number> <test_size> <combine> <combine 파일> 각 매개변수의 세부 정보는 다음과 같습니다. in-data-file: 입력 데이터 파일 in-label-file: 입력 데이터 파일의 레이블 output 폴더: 이 폴더에는 Output-score.xlsx (기능 선택 방법 및 해당 분류자의 정확도) 및 Output-SelectedFeatures.xlsx (각 기능 선택 알고리즘에 대해 선택한 기능 이름)를 포함하여 두 개의 출력 파일이 저장됩니다.<ol><li>select_feature_number: 1에서 데이터 파일의 기능 수까지 기능 수를 선택합니다.</li>
			<li>test_size: 분할할 테스트 샘플의 비율을 설정합니다. 예를 들어, 0.2는 입력 데이터셋이 0.8:0.2의 비율로 train: test 부분 집합으로 임의로 분할됨을 의미합니다.</li>
			<li>combine: true인 경우, 기능 선택(예: 원래 표현식 값 및 mqTrans 기능)을 위해 두 개의 데이터 파일을 함께 결합합니다. false인 경우, 기능 선택에 하나의 데이터 파일(예: 원래 표현식 값 또는 mqTrans 기능)만 사용하십시오.</li>
			<li>combine file: comparent가 true이면 이 파일 이름을 제공하여 결합된 데이터 행렬을 저장합니다. 참고: 이 파이프라인은 생성된 mqTrans 기능이 분류 태스크에서 어떻게 수행되는지 보여주는 것을 목표로 하며, 다음 조작을 위해 섹션 2에서 생성된 파일을 직접 사용합니다.</li>
		</ol></li>
	</ol></li>
	<li>mqTrans 기능 선택에 대한 기능 선택 알고리즘을 실행하십시오.<ol><li>사용자가 mqTrans 기능 또는 원래 기능을 선택하는 경우 결합 =False 를 설정하십시오.</li>
		<li>먼저 800개의 원래 기능을 선택하고 데이터 세트를 학습(test=0.8:0.2)으로 분할합니다. 파이썬 ./FS_classification/testMain.py ./output/test_target.csv ./output/label.csv ./result 800 0.2 거짓</li>
		<li>사용자가 mqTrans 기능을 원래 표현식 값과 결합하여 기능을 선택하려는 경우 combine =True를 설정하십시오. 여기서 시연 예제는 800개의 기능을 선택하고 데이터 세트를 train으로 분할하는 것입니다(test=0.8:0.2). 파이썬 ./FS_classification/testMain.py ./output/mq_targets.csv ./output/label.csv ./result_combine 800 0.2 True ./output/test_target.csv 참고: 그림 5 는 출력 정보를 보여줍니다. 이 프로토콜에 필요한 보충 파일은 HealthModel-mqTrans-v1-00.tar 폴더(보충 코딩 파일 1)에 있습니다.</li>
	</ol></li>
</ol>

트랜스크립톰에서 어두운 바이오마커를 식별하기 위한 GAT 프레임워크

<ol><li>Mutz, K. -O., Heilkenbrinker, A., Lönne, M., Walter, J. -G., Stahl, F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Transcriptome+analysis+using+next-generation+sequencing%5BTitle%5D)+AND+%22Curr+Opin+in+Biotechnol%22%5BJournal%5D)">Transcriptome analysis using next-generation sequencing.</a> Curr Opin in Biotechnol. 24 (1), 22-30 (2013).</li>
<li>Meng, G., Tang, W., Huang, E., Li, Z., Feng, H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+comprehensive+assessment+of+cell+type-specific+differential+expression+methods+in+bulk+data%5BTitle%5D)+AND+%22Brief+Bioinform%22%5BJournal%5D)">A comprehensive assessment of cell type-specific differential expression methods in bulk data.</a> Brief Bioinform. 24 (1), 516(2023).</li>
<li>Iqbal, N., Kumar, P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Integrated+COVID-19+Predictor%3A+Differential+expression+analysis+to+reveal+potential+biomarkers+and+prediction+of+coronavirus+using+RNA-Seq+profile+data%5BTitle%5D)+AND+%22Comput+Biol+Med%22%5BJournal%5D)">Integrated COVID-19 Predictor: Differential expression analysis to reveal potential biomarkers and prediction of coronavirus using RNA-Seq profile data.</a> Comput Biol Med. 147, 105684(2022).</li>
<li>Ravichandran, S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((VB%2810%29%2C+a+new+blood+biomarker+for+differential+diagnosis+and+recovery+monitoring+of+acute+viral+and+bacterial+infections%5BTitle%5D)+AND+%22EBioMedicine%22%5BJournal%5D)">VB(10), a new blood biomarker for differential diagnosis and recovery monitoring of acute viral and bacterial infections.</a> EBioMedicine. 67, 103352(2021).</li>
<li>Lv, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Targeting+FABP4+in+elderly+mice+rejuvenates+liver+metabolism+and+ameliorates+aging-associated+metabolic+disorders%5BTitle%5D)+AND+%22Metabolism%22%5BJournal%5D)">Targeting FABP4 in elderly mice rejuvenates liver metabolism and ameliorates aging-associated metabolic disorders.</a> Metabolism. 142, 155528(2023).</li>
<li>Cruz, J. A., Wishart, D. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Applications+of+machine+learning+in+cancer+prediction+and+prognosis%5BTitle%5D)+AND+%22Cancer+Inform%22%5BJournal%5D)">Applications of machine learning in cancer prediction and prognosis.</a> Cancer Inform. 2, 59-77 (2007).</li>
<li>Cox, D. R. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aCox%2C+DR+%22Analysis+of+Survival+Data%22">Analysis of Survival Data</a>. , Chapman and Hall/CRC. London. (2018).</li>
<li>Newman, A. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Robust+enumeration+of+cell+subsets+from+tissue+expression+profiles%5BTitle%5D)+AND+%22Nat+Methods%22%5BJournal%5D)">Robust enumeration of cell subsets from tissue expression profiles.</a> Nat Methods. 12 (5), 453-457 (2015).</li>
<li>Ramilowski, J. A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+draft+network+of+ligand-receptor-mediated+multicellular+signalling+in+human%5BTitle%5D)+AND+%22Nat+Commun%22%5BJournal%5D)">A draft network of ligand-receptor-mediated multicellular signalling in human.</a> Nat Commun. 6 (1), 7866(2015).</li>
<li>Xu, Y., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((MiR-145+detection+in+urinary+extracellular+vesicles+increase+diagnostic+efficiency+of+prostate+cancer+based+on+hydrostatic+filtration+dialysis+method%5BTitle%5D)+AND+%22Prostate%22%5BJournal%5D)">MiR-145 detection in urinary extracellular vesicles increase diagnostic efficiency of prostate cancer based on hydrostatic filtration dialysis method.</a> Prostate. 77 (10), 1167-1175 (2017).</li>
<li>Wang, Y., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Profiles+of+differential+expression+of+circulating+microRNAs+in+hepatitis+B+virus-positive+small+hepatocellular+carcinoma%5BTitle%5D)+AND+%22Cancer+Biomark%22%5BJournal%5D)">Profiles of differential expression of circulating microRNAs in hepatitis B virus-positive small hepatocellular carcinoma.</a> Cancer Biomark. 15 (2), 171-180 (2015).</li>
<li>Hu, S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Transcriptional+response+profiles+of+paired+tumor-normal+samples+offer+novel%5BTitle%5D)+AND+%22Oncotarget%22%5BJournal%5D)">Transcriptional response profiles of paired tumor-normal samples offer novel.</a> Oncotarget. 8 (25), 41334-41347 (2017).</li>
<li>Xu, H., Luo, D., Zhang, F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((DcWRKY75+promotes+ethylene+induced+petal+senescence+in+carnation+%28Dianthus+caryophyllus+L%29%5BTitle%5D)+AND+%22Plant+J%22%5BJournal%5D)">DcWRKY75 promotes ethylene induced petal senescence in carnation (Dianthus caryophyllus L).</a> Plant J. 108 (5), 1473-1492 (2021).</li>
<li>Niu, H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Dynamic+role+of+Scd1+gene+during+mouse+oocyte+growth+and+maturation%5BTitle%5D)+AND+%22Int+J+Biol+Macromol%22%5BJournal%5D)">Dynamic role of Scd1 gene during mouse oocyte growth and maturation.</a> Int J Biol Macromol. 247, 125307(2023).</li>
<li>Aznaourova, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Single-cell+RNA+sequencing+uncovers+the+nuclear+decoy+lincRNA+PIRAT+as+a+regulator+of+systemic+monocyte+immunity+during+COVID-19%5BTitle%5D)+AND+%22Proc+Natl+Acad+Sci+U+S+A%22%5BJournal%5D)">Single-cell RNA sequencing uncovers the nuclear decoy lincRNA PIRAT as a regulator of systemic monocyte immunity during COVID-19.</a> Proc Natl Acad Sci U S A. 119 (36), 2120680119(2022).</li>
<li>Prakash, A., Banerjee, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((An+interpretable+block-attention+network+for+identifying+regulatory+feature+interactions%5BTitle%5D)+AND+%22Brief+Bioinform%22%5BJournal%5D)">An interpretable block-attention network for identifying regulatory feature interactions.</a> Brief Bioinform. 24 (4), (2023).</li>
<li>Zhai, Y., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Single-cell+RNA+sequencing+integrated+with+bulk+RNA+sequencing+analysis+reveals+diagnostic+and+prognostic+signatures+and+immunoinfiltration+in+gastric+cancer%5BTitle%5D)+AND+%22Comput+Biol+Med%22%5BJournal%5D)">Single-cell RNA sequencing integrated with bulk RNA sequencing analysis reveals diagnostic and prognostic signatures and immunoinfiltration in gastric cancer.</a> Comput Biol Med. 163, 107239(2023).</li>
<li>Duan, L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Dynamic+changes+in+spatiotemporal+transcriptome+reveal+maternal+immune+dysregulation+of+autism+spectrum+disorder%5BTitle%5D)+AND+%22Comput+Biol+Med%22%5BJournal%5D)">Dynamic changes in spatiotemporal transcriptome reveal maternal immune dysregulation of autism spectrum disorder.</a> Comput Biol Med. 151, 106334(2022).</li>
<li>Zolotareva, O., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Flimma%3A+a+federated+and+privacy-aware+tool+for+differential+gene+expression+analysis%29%5BTitle%5D)+AND+%22Genome+Biol%22%5BJournal%5D)">Flimma: a federated and privacy-aware tool for differential gene expression analysis).</a> Genome Biol. 22 (1), 338(2021).</li>
<li>Su, R., Zhu, Y., Zou, Q., Wei, L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Distant+metastasis+identification+based+on+optimized+graph+representation+of+gene%5BTitle%5D)+AND+%22Brief+Bioinform%22%5BJournal%5D)">Distant metastasis identification based on optimized graph representation of gene.</a> Brief Bioinform. 23 (1), (2022).</li>
<li>Xing, X., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Multi-level+attention+graph+neural+network+based+on+co-expression+gene+modules+for+disease+diagnosis+and+prognosis%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">Multi-level attention graph neural network based on co-expression gene modules for disease diagnosis and prognosis.</a> Bioinformatics. 38 (8), 2178-2186 (2022).</li>
<li>Bongini, P., Pancino, N., Scarselli, F., Bianchini, M. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aBongini%2C+P+%22BioGNN%3A+How+Graph+Neural+Networks+Can+Solve+Biological+Problems.+Artificial+Intelligence+and+Machine+Learning+for+Healthcare%3A+Vol.+1%3A+Image+and+Data+Analytics%22">BioGNN: How Graph Neural Networks Can Solve Biological Problems. Artificial Intelligence and Machine Learning for Healthcare: Vol. 1: Image and Data Analytics</a>. , Springer. Cham. (2022).</li>
<li>Muzio, G., O'Bray, L., Borgwardt, K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Biological+network+analysis+with+deep+learning%5BTitle%5D)+AND+%22Brief+Bioinform%22%5BJournal%5D)">Biological network analysis with deep learning.</a> Brief Bioinform. 22 (2), 1515-1530 (2021).</li>
<li>Luo, H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Multi-omics+integration+for+disease+prediction+via+multi-level+graph+attention+network+and+adaptive+fusion%5BTitle%5D)+AND+%22bioRxiv%22%5BJournal%5D)">Multi-omics integration for disease prediction via multi-level graph attention network and adaptive fusion.</a> bioRxiv. , (2023).</li>
<li>Feng, X., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Selecting+multiple+biomarker+subsets+with+similarly+effective+binary+classification+performances%5BTitle%5D)+AND+%22J+Vis+Exp%22%5BJournal%5D)">Selecting multiple biomarker subsets with similarly effective binary classification performances.</a> J Vis Exp. (140), e57738(2018).</li>
<li>Duan, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Orchestrating+information+across+tissues+via+a+novel+multitask+GAT+framework+to+improve+quantitative+gene+regulation+relation+modeling+for+survival+analysis%5BTitle%5D)+AND+%22Brief+Bioinform%22%5BJournal%5D)">Orchestrating information across tissues via a novel multitask GAT framework to improve quantitative gene regulation relation modeling for survival analysis.</a> Brief Bioinform. 24 (4), (2023).</li>
<li>Chicco, D., Starovoitov, V., Jurman, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+benefits+of+the+Matthews+correlation+Coefficient+%28MCC%29+over+the+diagnostic+odds+ratio+%28DOR%29+in+binary+classification+assessment%5BTitle%5D)+AND+%22IEEE+Access%22%5BJournal%5D)">The benefits of the Matthews correlation Coefficient (MCC) over the diagnostic odds ratio (DOR) in binary classification assessment.</a> IEEE Access. 9, 47112-47124 (2021).</li>
<li>Goldman, M. J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Visualizing+and+interpreting+cancer+genomics+data+via+the+Xena+platform%5BTitle%5D)+AND+%22Nat+Biotechnol%22%5BJournal%5D)">Visualizing and interpreting cancer genomics data via the Xena platform.</a> Nat Biotechnol. 38 (6), 675-678 (2020).</li>
<li>Liu, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((An+integrated+TCGA+pan-cancer+clinical+data+resource+to+drive+high-quality+survival+outcome+analytics%5BTitle%5D)+AND+%22Cell%22%5BJournal%5D)">An integrated TCGA pan-cancer clinical data resource to drive high-quality survival outcome analytics.</a> Cell. 173 (2), 400-416 (2018).</li>
<li>Uhlen, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Towards+a+knowledge-based+human+protein+atlas%5BTitle%5D)+AND+%22Nat+Biotechnol%22%5BJournal%5D)">Towards a knowledge-based human protein atlas.</a> Nat Biotechnol. 28 (12), 1248-1250 (2010).</li>
<li>Hernaez, M., Blatti, C., Gevaert, O. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Comparison+of+single+and+module-based+methods+for+modeling+gene+regulatory%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">Comparison of single and module-based methods for modeling gene regulatory.</a> Bioinformatics. 36 (2), 558-567 (2020).</li>
<li>Consortium, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+genotype-tissue+expression+%28GTEx%29+project%5BTitle%5D)+AND+%22Nat+Genet%22%5BJournal%5D)">The genotype-tissue expression (GTEx) project.</a> Nat Genet. 45 (6), 580-585 (2013).</li>
<li>Kanehisa, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((KEGG+for+taxonomy-based+analysis+of+pathways+and+genomes%5BTitle%5D)+AND+%22Nucleic+Acids+Res%22%5BJournal%5D)">KEGG for taxonomy-based analysis of pathways and genomes.</a> Nucleic Acids Res. 51, D587-D592 (2023).</li>
<li>Han, H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((TRRUST+v2%3A+an+expanded+reference+database+of+human+and+mouse+transcriptional%5BTitle%5D)+AND+%22Nucleic+Acids+Res%22%5BJournal%5D)">TRRUST v2: an expanded reference database of human and mouse transcriptional.</a> Nucleic Acids Res. 46, D380-D386 (2018).</li>
<li>Langfelder, P., Horvath, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((WGCNA%3A+an+R+package+for+weighted+correlation+network+analysis%5BTitle%5D)+AND+%22BMC+Bioinformatics%22%5BJournal%5D)">WGCNA: an R package for weighted correlation network analysis.</a> BMC Bioinformatics. 9, 559(2008).</li>
<li>Sulaimanov, N., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Inferring+gene+expression+networks+with+hubs+using+a+degree+weighted+Lasso%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">Inferring gene expression networks with hubs using a degree weighted Lasso.</a> Bioinformatics. 35 (6), 987-994 (2019).</li>
<li>Kogelman, L. J. A., Kadarmideen, H. N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Weighted+Interaction+SNP+Hub+%28WISH%29+network+method+for+building+genetic+networks%5BTitle%5D)+AND+%22BMC+Syst+Biol%22%5BJournal%5D)">Weighted Interaction SNP Hub (WISH) network method for building genetic networks.</a> BMC Syst Biol. 8, 5(2014).</li>
<li>Duan, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Pan-cancer+identification+of+the+relationship+of+metabolism-related+differentially+expressed+transcription+regulation+with+non-differentially+expressed+target+genes+via+a+gated+recurrent+unit+network%5BTitle%5D)+AND+%22Comput+Biol+Med%22%5BJournal%5D)">Pan-cancer identification of the relationship of metabolism-related differentially expressed transcription regulation with non-differentially expressed target genes via a gated recurrent unit network.</a> Comput Biol Med. 148, 105883(2022).</li>
<li>Duan, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Detection+and+independent+validation+of+model-based+quantitative+transcriptional+regulation+relationships+altered+in+lung+cancers%5BTitle%5D)+AND+%22Front+Bioeng+Biotechnol%22%5BJournal%5D)">Detection and independent validation of model-based quantitative transcriptional regulation relationships altered in lung cancers.</a> Front Bioeng Biotechnol. 8, 582(2020).</li>
<li>Fiorini, N., Lipman, D. J., Lu, Z. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Towards+PubMed+2.0%5BTitle%5D)+AND+%22eLife%22%5BJournal%5D)">Towards PubMed 2.0.</a> eLife. 6, 28801(2017).</li>
<li>Liu, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Maternal+microbiome+regulation+prevents+early+allergic+airway+diseases+in+mouse+offspring%5BTitle%5D)+AND+%22Pediatr+Allergy+Immunol%22%5BJournal%5D)">Maternal microbiome regulation prevents early allergic airway diseases in mouse offspring.</a> Pediatr Allergy Immunol. 31 (8), 962-973 (2020).</li>
<li>Childs, E. J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Association+of+common+susceptibility+variants+of+pancreatic+cancer+in+higher-risk+patients%3A+A+PACGENE+study%5BTitle%5D)+AND+%22Cancer+Epidemiol+Biomarkers+Prev%22%5BJournal%5D)">Association of common susceptibility variants of pancreatic cancer in higher-risk patients: A PACGENE study.</a> Cancer Epidemiol Biomarkers Prev. 25 (7), 1185-1191 (2016).</li>
<li>Wang, C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Thailandepsins%3A+bacterial+products+with+potent+histone+deacetylase+inhibitory+activities+and+broad-spectrum+antiproliferative+activities%5BTitle%5D)+AND+%22J+Nat+Prod%22%5BJournal%5D)">Thailandepsins: bacterial products with potent histone deacetylase inhibitory activities and broad-spectrum antiproliferative activities.</a> J Nat Prod. 74 (10), 2031-2038 (2011).</li>
<li>Lv, X., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Transcriptional+dysregulations+of+seven+non-differentially+expressed+genes+as+biomarkers+of+metastatic+colon+cancer%5BTitle%5D)+AND+%22Genes+%28Basel%29%22%5BJournal%5D)">Transcriptional dysregulations of seven non-differentially expressed genes as biomarkers of metastatic colon cancer.</a> Genes (Basel). 14 (6), 1138(2023).</li>
<li>Li, X., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Undifferentially+expressed+CXXC5+as+a+transcriptionally+regulatory+biomarker+of+breast+cancer%5BTitle%5D)+AND+%22Advanced+Biology%22%5BJournal%5D)">Undifferentially expressed CXXC5 as a transcriptionally regulatory biomarker of breast cancer.</a> Advanced Biology. , (2023).</li>
<li>Yuan, W., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+N6-methyladenosine+reader+protein+YTHDC2+promotes+gastric+cancer+progression+via+enhancing+YAP+mRNA+translation%5BTitle%5D)+AND+%22Transl+Oncol%22%5BJournal%5D)">The N6-methyladenosine reader protein YTHDC2 promotes gastric cancer progression via enhancing YAP mRNA translation.</a> Transl Oncol. 16, 101308(2022).</li>
<li>Tanabe, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((RNA+helicase+YTHDC2+promotes+cancer+metastasis+via+the+enhancement+of+the+efficiency+by+which+HIF-1%CE%B1+mRNA+is+translated%5BTitle%5D)+AND+%22Cancer+Lett%22%5BJournal%5D)">RNA helicase YTHDC2 promotes cancer metastasis via the enhancement of the efficiency by which HIF-1α mRNA is translated.</a> Cancer Lett. 376 (1), 34-42 (2016).</li>
</ol>

저자는 공개할 것이 없습니다.

프로토콜의 섹션 2(사전 학습된 HealthModel을 사용하여 mqTrans 기능 생성)는 이 프로토콜 내에서 가장 중요한 단계입니다. 섹션 1에서 계산 작업 환경을 준비한 후 섹션 2에서는 사전 훈련된 대형 참조 모델을 기반으로 전사체 데이터 세트의 mqTrans 보기를 생성합니다. 섹션 3은 바이오마커 검출 및 예측 작업을 위해 생성된 mqTrans 기능을 선택하는 데모 예제입니다. 사용자는 자신의 도구 또는 코드를 사용하여 이 mqTrans 데이터 세트에 대해 다른 전사체 분석을 수행할 수 있습니다.
원래 HealthModel 프레임워크는26에 설명된 대로 멀티태스킹 아키텍처를 사용하여 사전 훈련된 HealthModel을 더욱 구체화할 수 있습니다. 이 프로토콜은 전사체 데이터 세트의 mqTrans 보기를 생성하기 위해 사전 훈련된 참조 모델의 활용에 중점을 둡니다.
기본 사전 훈련된 참조 모델은 건강한 샘플에 대해 설정되었으며 일부 특정 작업(예: 원발암과 전이성 암 간의 조사)에는 적합하지 않을 수 있습니다. 또한 대규모 전사체 데이터 세트의 경우 계산 속도가 느립니다.
이 프로토콜의 중요성은 가장 풍부하게 사용 가능한 OMIC 데이터 유형(예: 트랜스크립톰)에 대한 보완적인 mqTrans 보기를 제공하는 것입니다. 다크 바이오마커는 기존의 전사체 분석에서 무시된 미분화 발현 유전자에서 밝힐 수 있습니다. 최근 연구에서는 총44개 샘플로 구성된 805개 샘플로 구성된 3개의 독립적인 코호트를 기반으로 전이성 대장암(mCC)의 7가지 어두운 바이오마커를 검출했습니다. 다크 바이오마커는 비차별적 발현으로 인해 제한된 습식 실험실 조사를 받았습니다. 그러나, 검출된 mCC 다크 바이오마커 YTHDC2 중 하나는 2를 포함하는 단백질 YTH 도메인을 암호화하며, 이의 단백질 수준은 인간 위암 세포(46 ) 및 결장암(47)의 전이 상태와 양의 상관관계가 있는 것으로 관찰되었다. 다크 바이오마커에 대한 새로운 생물학적 통찰력은 in vitro 및 in vivo 기술을 통해 해결되어야 합니다.
이 프로토콜은 완전히 모듈식으로 설계되었습니다. 원발암과 같은 다른 대규모 데이터 세트에서 사전 훈련된 참조 모델은 종양 전이 조사를 용이하게 합니다. 이 프로토콜은 식물, 균류 및 미생물을 포함한 다른 생명 영역의 응용 분야에서도 탐구될 것입니다.
이 프로토콜의 계산 효율성은 병렬화 및 알고리즘 최적화를 통해 향상될 계획입니다.
이 프로토콜은 전사체 데이터 세트를 새로운 mqTrans 보기로 변환하는 절차를 설명하며, 유전자의 변환된 mqTrans 값은 참조 샘플과 비교하여 전사 조절 변화를 정량적으로 측정합니다. 기본 모델은 건강한 전사체에 대해 사전 학습되었으며 참조 HealthModel로 릴리스되었습니다.
생물의학 연구자가 이 프로토콜을 쉽게 활용할 수 있도록 두 가지 다운스트림 작업의 소스 코드가 제공됩니다. 실험 데이터는 변환된 mqTrans 기능이 원래 표현식 레벨만 사용하여 예측 태스크를 개선할 수 있음을 보여줍니다. mqTrans 보기는 또한 원본 전사체 데이터에서 차등 발현 없이 일부 어두운 바이오마커의 잠재 표현형 연결을 밝힐 수 있습니다.

전사체는 샘플에 있는 모든 유전자의 발현으로 구성되며 마이크로어레이 및 RNA-seq1과 같은 고처리량 기술로 프로파일링할 수 있습니다. 데이터 세트에서 한 유전자의 발현 수준을 전사체 특징이라고 하며, 표현형과 대조군 간의 전사체 특징의 차등 표현은 이 유전자를 이 표현형 2,3의 바이오마커로 정의합니다. 전사체 바이오마커는 질병 진단4, 생물학적 메커니즘5, 생존 분석6,7 등의 연구에 광범위하게 활용되어 왔다.
건강한 조직의 유전자 활동 패턴은 생명에 대한 중요한 정보를 담고 있다 8,9. 이러한 패턴은 귀중한 통찰력을 제공하며, 양성 질환(benign disorders)10,11 및 치명적 질환(lethal diseases)12의 복잡한 발달 궤적을 이해하기 위한 이상적인 참고 자료 역할을 한다. 유전자는 서로 상호 작용하며, 전사체는 복잡한 상호 작용 후 최종 발현 수준을 나타냅니다. 이러한 패턴은 전사 조절 네트워크(transcriptional regulation network)13 및 대사 네트워크(metabolism network)14 등으로 공식화된다. 메신저 RNA(mRNA)의 발현은 전사 인자(TF) 및 긴 유전자 간 비코딩 RNA(lincRNA)에 의해 전사적으로 조절될 수 있습니다15,16,17. 종래의 차등 발현 분석은 특징간 독립성(inter-feature independence)을 가정한 이러한 복잡한 유전자 상호작용을 무시하였다18,19.
그래프 신경망(GNN)의 최근 발전은 암 연구(20)를 위한 OMIC 기반 데이터로부터 중요한 정보를 추출하는 데 있어 비범한 잠재력을 보여주는데, 예를 들어, 공발현 모듈(21)을 식별한다. GNN의 타고난 능력은 유전자22,23 사이의 복잡한 관계와 의존성을 모델링하는 데 이상적입니다.
생물의학 연구는 종종 대조군에 대한 표현형을 정확하게 예측하는 데 중점을 둡니다. 이러한 작업은 일반적으로 이진 분류(binary classifications)24,25,26으로 공식화된다. 여기서 두 클래스 레이블은 일반적으로 1과 0, true 및 false 또는 양수 및 음수27로 인코딩됩니다.
이 연구는 사전 훈련된 GAT(Graph-Attention Network) 참조 모델을 기반으로 전사체 데이터 세트의 전사 조절(mqTrans) 보기를 생성하기 위한 사용하기 쉬운 프로토콜을 제공하는 것을 목표로 했습니다. 이전에 공개된 작업(26 )으로부터의 멀티태스킹 GAT 프레임워크를 사용하여 전사체 특징을 mqTrans 특징들로 변환하였다. UCSC(University of California, Santa Cruz) Xena 플랫폼28 의 건강한 전사체에 대한 대규모 데이터 세트를 사용하여 조절 인자(TF 및 lincRNA)에서 표적 mRNA까지의 전사 조절을 정량적으로 측정한 참조 모델(HealthModel)을 사전 훈련했습니다. 생성된 mqTrans 뷰는 예측 모델을 구축하고 다크 바이오마커를 감지하는 데 사용할 수 있습니다. 이 프로토콜은 TCGA(The Cancer Genome Atlas) 데이터베이스(29 )의 결장 선암종(COAD) 환자 데이터 세트를 예시로 사용합니다. 이러한 맥락에서 I기 또는 II기의 환자는 음성 샘플로 분류되고 III 또는 IV기의 환자는 양성 샘플로 간주됩니다. 26개의 TCGA 암 유형에 대한 어두운 바이오마커와 전통적인 바이오마커의 분포도 비교됩니다.
HealthModel 파이프라인에 대한 설명 이 프로토콜에 사용된 방법론은 그림 1에 요약된 바와 같이 이전에 발표된 프레임워크26을 기반으로 합니다. 시작하려면 사용자가 입력 데이터 세트를 준비하고, 제안된 HealthModel 파이프라인에 공급하고, mqTrans 기능을 가져와야 합니다. 자세한 데이터 준비 지침은 프로토콜 섹션의 섹션 2에 나와 있습니다. 그 후, 사용자는 mqTrans 기능을 원래 전사체 기능과 결합하거나 생성된 mqTrans 기능만 진행할 수 있습니다. 생성된 데이터 세트는 기능 선택 프로세스를 거치며, 사용자는 분류를 위한 k-폴드 교차 검증에서 k에 대해 선호하는 값을 유연하게 선택할 수 있습니다. 이 프로토콜에서 사용되는 기본 평가 메트릭은 정확도입니다.
HealthModel26은 전사체 특징을 TF(전사 인자), lincRNA(긴 유전자 간 비코딩 RNA) 및 mRNA(메신저 RNA)의 세 가지 그룹으로 분류합니다. TF 특징은 Human Protein Atlas30,31에서 사용할 수 있는 주석을 기반으로 정의됩니다. 이 작업은 GTEx 데이터 세트32의 lincRNA 주석을 활용합니다. KEGG 데이터베이스(33)에서 3단계 경로에 속하는 유전자는 mRNA 특징으로서 고려된다. mRNA 특징이 TRRUST 데이터베이스(34)에 문서화된 바와 같이 표적 유전자에 대한 조절 역할을 나타내면, TF 클래스로 재분류된다는 점에 주목할 필요가 있다.
이 프로토콜은 또한 조절 인자(regulatory_geneIDs.csv)와 표적 mRNA(target_geneIDs.csv)의 유전자 ID에 대한 두 개의 예제 파일을 수동으로 생성합니다. 조절 특징(TF 및 lincRNA) 간의 쌍별 거리 매트릭스는 Pearson 상관 계수에 의해 계산되고 널리 사용되는 도구인 WGCNA(Weighted Gene Co-Expression Network Analysis)36 (adjacent_matrix.csv)에 의해 클러스터링됩니다. 사용자는 이러한 예제 구성 파일과 함께 HealthModel 파이프라인을 직접 활용하여 전사체 데이터 세트의 mqTrans 보기를 생성할 수 있습니다.
HealthModel의 기술적 세부 사항 HealthModel은 TF와 lincRNA 간의 복잡한 관계를 그래프로 나타내며, 입력 특징은 V로 표시된 꼭짓점과 E로 지정된 꼭짓점 간 가장자리 행렬 역할을 합니다. 각 샘플은 VK×1로 기호화된 K 조절 특징이 특징입니다. 특히, 데이터 세트에는 425개의 TF와 375개의 lincRNA가 포함되어 있어 K = 425 + 375 = 800의 샘플 차원을 얻을 수 있었습니다. 에지 매트릭스 E를 확립하기 위해 이 작업에서는 널리 사용되는 도구인 WGCNA35를 사용했습니다. 및 <img alt="Equation 2" src="/files/ftp_upload/66030/66030eq02.jpg" style="margin: auto;">로 표시된 <img alt="Equation 1" src="/files/ftp_upload/66030/66030eq01.jpg" style="margin: auto;"> 두 꼭짓점을 연결하는 쌍별 가중치는 Pearson 상관 계수에 의해 결정됩니다. 유전자 조절 네트워크(gene regulatory network)는 중추적인 기능적 역할을 하는 허브 유전자(hub gene)의 존재를 특징으로 하는 비늘 없는 토폴로지(scale-free topology)(36)를 나타낸다. 두 피처 또는 꼭짓점 간의 상관 관계를 계산하고<img alt="Equation 2" src="/files/ftp_upload/66030/66030eq02.jpg" style="margin: auto;">, <img alt="Equation 1" src="/files/ftp_upload/66030/66030eq01.jpg" style="margin: auto;"> 다음과 같이 TOM(Topological Overlap Measure)을 사용합니다.
<img alt="Equation 3" src="/files/ftp_upload/66030/66030eq03.jpg" style="margin: auto;">(1)
<img alt="Equation 4" src="/files/ftp_upload/66030/66030eq04.jpg" style="margin: auto;">(2)
소프트 임계값 β는 WGCNA 패키지의 'pickSoft Threshold' 함수를 사용하여 계산됩니다. 거듭제곱 지수 함수 a ij가 적용되며, 여기서 <img alt="Equation 5" src="/files/ftp_upload/66030/66030eq05.jpg" style="margin: auto;"> 는 i와 j를 제외한 유전자를 나타내고 <img alt="Equation 6" src="/files/ftp_upload/66030/66030eq06.jpg" style="margin: auto;"> 꼭짓점 연결을 나타냅니다. WGCNA는 일반적으로 사용되는 비유사성 측정(dissimilarity measure)을 사용하여 전사체 특징의 발현 프로파일을 여러 모듈로 클러스터링합니다(<img alt="Equation 7" src="/files/ftp_upload/66030/66030eq07.jpg" style="margin: auto;">37.
HealthModel 프레임워크는 원래 멀티태스킹 학습 아키텍처(multitask learning architecture)로 설계되었다(26). 이 프로토콜은 트랜스크립토믹 mqTrans 보기의 구성을 위해 모델 사전 학습 작업만 사용합니다. 사용자는 추가 작업별 전사체 샘플을 사용하여 멀티태스킹 그래프 주의 네트워크 하에서 사전 훈련된 HealthModel을 더욱 구체화하도록 선택할 수 있습니다.
기능 선택 및 분류에 대한 기술적 세부 사항 기능 선택 풀은 11개의 FS(기능 선택) 알고리즘을 구현합니다. 그 중 SK_mic(Maximal Information Coefficient)를 사용하여 K 최상의 특징을 선택하고, MIC의 FPR을 기반으로 K 특징을 선택하고(SK_fpr), MIC의 오탐지율이 가장 높은 K 특징을 선택(SK_fdr)하는 세 가지 필터 기반 FS 알고리즘입니다. 또한 세 가지 트리 기반 FS 알고리즘은 지니 지수(DT_gini), 적응형 부스트 의사 결정 트리(AdaBoost) 및 랜덤 포레스트(RF_fs)가 있는 의사 결정 트리를 사용하여 개별 기능을 평가합니다. 또한 풀은 선형 서포트 벡터 분류자(RFE_SVC)를 사용한 재귀적 기능 제거와 로지스틱 회귀 분류자(RFE_LR)를 사용한 재귀적 기능 제거라는 두 가지 래퍼 방법을 통합합니다. 마지막으로, 최상위 L1 기능 중요도 값(lSVC_L1)이 있는 선형 SVC 분류자와 최상위 L1 기능 중요도 값(LR_L1)이 있는 로지스틱 회귀 분류자의 두 가지 임베딩 알고리즘이 포함됩니다.
분류자 풀은 7개의 서로 다른 분류자를 사용하여 분류 모델을 빌드합니다. 이러한 분류기는 선형 서포트 벡터 머신(SVC), 가우스 나이브 베이즈(GNB), 로지스틱 회귀 분류기(LR), k-최근접이웃(k가 기본적으로 5로 설정됨(KNN), XGBoost, 랜덤 포레스트(RF) 및 의사 결정 트리(DT)로 구성됩니다.
데이터셋을 train: test 서브셋으로 임의로 분할하는 것은 명령줄에서 설정할 수 있습니다. 시연된 예제에서는 train: test = 8:2의 비율을 사용합니다.

<table><tbody><tr><td>아나콘다</td><td>아나콘다</td><td>버전 2020.11</td><td>Python 프로그래밍 플랫폼</td></tr><tr><td>컴퓨터</td><td>N/A</td><td>N/A</td><td>요구 사항을 충족하는 모든 범용 컴퓨터</td></tr><tr>GPU<td>카드</td><td>N/A</td> N<td>/A</td><td>CUDA 컴퓨팅 라이브러리가 있는 모든 범용 GPU 카드</td></tr><tr><td>pytorch</td><td>Pytorch</td><td>버전 1.13.1</td><td>소프트웨어</td></tr><tr><td>토치 기하학</td><td>Pytorch</td><td>버전 2.2.0</td><td>소프트웨어</td></tr></tbody></table>

작은 데이터 세트에서 예측 작업 및 다크 바이오마커 검출을 위한 전사체 특징의 전사체 조절 보기 생성

Generating the Transcriptional Regulation View of Transcriptomic Features for Prediction Task and Dark Biomarker Detection on Small Datasets

전사체는 샘플에서 많은 유전자의 발현 수준을 나타내며 생물학 연구 및 임상 실습에서 널리 사용되었습니다. 연구자들은 일반적으로 표현형 그룹과 샘플의 대조군 사이에 차이 표현을 가진 전사체 바이오마커에 초점을 맞췄습니다. 이 연구는 참조 샘플의 복잡한 유전자 간 상호 작용을 학습하기 위한 멀티태스킹 그래프-주의 네트워크(GAT) 학습 프레임워크를 제시했습니다. 시범적 참조 모델은 건강한 샘플(HealthModel)에 대해 사전 학습되었으며, 이는 독립적인 테스트 전사체의 모델 기반 정량적 전사 조절(mqTrans) 보기를 생성하는 데 직접 사용할 수 있습니다. 생성된 전사체의 mqTrans 보기는 예측 작업 및 다크 바이오마커 검출에 의해 입증되었습니다. "다크 바이오마커(dark biomarker)"라는 용어는 다크 바이오마커가 mqTrans 뷰에서 차등 표현을 나타내지만 원래 발현 수준에서는 차등 발현을 나타내지 않는다는 정의에서 비롯되었습니다. 어두운 바이오마커는 차등 발현이 없기 때문에 전통적인 바이오마커 검출 연구에서 항상 간과되었습니다. 파이프라인 HealthModelPipe의 소스 코드 및 매뉴얼은 http://www.healthinformaticslab.org/supp/resources.php 에서 다운로드할 수 있습니다.

전사체 데이터 세트의 mqTrans 보기 평가 테스트 코드는 11개의 기능 선택(FS) 알고리즘과 7개의 분류자를 사용하여 전사체 데이터 세트의 생성된 mqTrans 보기가 분류 작업에 어떻게 기여하는지 평가합니다(그림 6). 테스트 데이터 세트는 TCGA(The Cancer Genome Atlas) 데이터베이스29의 317개 결장 선암종(COAD)으로 구성됩니다. I기 또는 II기의 COAD 환자는 음성 샘플로 간주되고 III 또는 IV기의 환자는 양성 샘플로 간주됩니다.
11개의 FS 알고리즘이 테스트 코드에 구현됩니다. MIC에 의한 K 최상의 특징 선택(SK_mic), MIC의 FPR에 의한 K 특징 선택(SK_fpr), MIC의 가장 높은 FDR에 의한 K 특징 선택(SK_fpr)을 포함한 세 가지 필터 기반 FS 알고리즘이 있습니다. 세 가지 트리 기반 FS 알고리즘은 각각 지니 지수(DT_gini), 적응형 부스트 의사결정 트리(AdaBoost) 및 랜덤 포레스트(RF_fs)를 사용하여 개별 기능을 평가합니다. 또한 테스트 코드의 FS 풀은 선형 서포트 벡터 분류자(SVC)(RFE_SVC)를 사용하는 RFE(재귀 기능 제거)와 로지스틱 회귀 분류자(RFE_LR)를 사용하는 RFE, 그리고 최상위 L1 기능 중요도 값(lSVC_L1)을 사용하는 선형 SVC 분류자 및 최상위 L1 기능 중요도 값(LR_L1)을 사용하는 로지스틱 회귀 분류자 등 두 개의 래퍼를 평가합니다.
테스트 코드는 선형 서포트 벡터 머신(SVC), 가우스 나이브 베이즈(GNB), 로지스틱 회귀 분류기(LR), k-최근접 이웃, k-5 기본값(KNN), XGBoost, 랜덤 포레스트(RF) 및 의사 결정 트리(DT)를 포함한 7개의 분류기를 사용하여 분류 모델을 빌드합니다.
그림 6은 mqTrans 특징, 원래 mRNA 특징, 각 FS 알고리즘에서 권장하는 mRNA 및 mqTrans 특징의 결합된 하위 집합의 최대 테스트 정확도를 보여줍니다.
결합된 특징 서브세트(mRNA+mqTrans)는 "SK_fpr" FS 분석법에서 가장 높은 정확도인 0.7656을 달성했으며, 이는 개별 특징 유형 mqTrans(0.7188) 및 오리지널 mRNA(0.7188)보다 우수합니다. 다른 FS 알고리즘에서도 유사한 패턴을 관찰할 수 있습니다. 사용자는 출력 파일 Output-SelectedFeatures.csv에서 선택한 기능을 확인할 수 있습니다.
다크 바이오마커 검출 이전 연구에서는 표현형 및 대조군 간에 현저하게 다르게 표현된 mqTrans 값을 갖는 미분화 발현 유전자의 존재를 보여주었습니다 26,38,39. 이러한 유전자는 다크 바이오마커(dark biomarker)라고 불리는데, 이는 전통적인 바이오마커 검출 연구가 비차별적 발현으로 이를 무시하기 때문입니다. Microsoft Excel의 통계 분석 함수 t.test를 사용하여 통계적 p-값이 0.05보다 작은 경우 차등적으로 표현되는 특성을 정의할 수 있습니다.
생성된 mqTrans 값을 가진 3062개의 특징 중 221개의 다크 바이오마커가 검출되었습니다(그림 7). 3위를 차지한 유전자 ENSG00000163697(APBB2, Amyloid Beta Precursor Protein Binding Family B Member 2)은 유의미하게 다르게 표현된 mqTrans 값(mqTrans.P = 2.03 x 10-4)을 보인 반면, 원래 발현 수준은 차등 발현을 나타내지 않았습니다(mRNA.P = 3.80 x 10-1). 키워드 APBB2는 PubMed 데이터베이스40개에서 27개의 간행물에 적중했으나 결장이나 장과의 연관성은 발견되지 않았다.
또 다른 유전자 ENSG00000048052(HDAC9, 히스톤 탈아세틸화효소 9)은 표현형과 대조군 간에 실질적으로 동일한 정규 분포(mRNA.P = 9.62 x 10-1)를 유지하면서 다르게 표현된 mqTrans 값(mqTrans.P = 6.09 x 10-3)을 갖습니다. 키워드 HDAC9은 PubMed 데이터베이스에서 417개의 출판물을 기록했습니다. 3건의 연구는 초록41, 42, 43에서 "결장(colon)" 또는 "장(intestine)"이라는 키워드를 언급했다. 그러나 그들 중 누구도 대장암에서 HDAC9의 역할을 조사하지 않았다.
이 데이터는 전사 후 활동, 예를 들어 번역된 단백질 수준44,45에서 이러한 어두운 바이오마커의 추가 평가의 필요성을 시사했습니다.
대사 관련 다크 및 전통적 바이오마커의 범암 분포 대사 관련 전통적 바이오마커를 스크리닝하여 TCGA 데이터 세트38의 26개 암 유형에 걸쳐 어두운 바이오마커와 비교했습니다. 두 범주의 바이오마커 모두 초기(I기 및 II) 및 후기(III기 및 IV) 암 단계에서 유의성 수준을 식별하기 위해 통계적 평가를 거쳤습니다. 이 평가는 p-값에 대해 스튜던트 t-검정을 사용했으며, 이후 FDR(False Discovery Rate)을 사용하여 여러 테스트에 대해 수정되었습니다. 26가지 암 유형 각각에 대한 자세한 데이터는 그림 8에 나와 있습니다.
FDR 보정된 p-값이 0.05 미만인 유전자는 전통적인 바이오마커로 분류되었습니다. 대조적으로, 다크 바이오마커는 mqTrans 보기에서 FDR 보정된 p-값이 0.05 미만인 것으로 정의되었으며, 동시에 발현 수준에서 통계적으로 유의한 차이를 나타내지 않았습니다.
그림 9는 대부분의 암 유형에서 전통적인 바이오마커와 비교하여 다크 바이오마커가 일반적으로 부족하다는 것을 보여줍니다. 주목할 만한 예외로는 BRCA, MESO 및 TGCT가 있으며, 이들은 다크 바이오마커의 더 큰 유병률을 나타냅니다. 전사 인자, 메틸화 패턴, 유전자 돌연변이 및 환경 조건을 포함한 다양한 요인이 이러한 어두운 바이오마커의 전사 조절 장애를 조절할 수 있음이 밝혀졌습니다. 다크 바이오마커의 발현 수준을 혼동할 수 있는 겹치는 비코딩 RNA 전사체로 인해 더 많은 복잡성이 발생할 수 있습니다. 일부 다크 바이오마커의 전사 조절 장애는 차등 단백질 수준44,45에 의해 뒷받침되었습니다. 다크 바이오마커는 전통적인 연구에서 종종 간과되며 미래의 기계론적 연구를 위한 흥미로운 길을 제시합니다.
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/66030/66030fig01.jpg"> 그림 1: 이 프로토콜의 HealthModel 및 기능 선택 모듈에 대한 개요입니다. 사용자가 Python 프로그래밍에 익숙한 경우 기능 선택 풀과 분류자 풀의 특정 알고리즘을 바꿉니다. <a href="https://www.jove.com/files/ftp_upload/66030/66030fig01large.jpg" target="_blank">이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.</a>
<img alt="Figure 2" class="xfigimg" src="/files/ftp_upload/66030/66030fig02.jpg"> 그림 2: 이 프로토콜에 대한 전체 코드 흐름. (a) Python 환경을 준비합니다. 시작하려면 가상 환경을 만들고 필수 패키지를 설치합니다. 자세한 지침은 섹션 1을 참조하십시오. (b) mqTrans 기능을 생성합니다. 제공된 코드를 단계별로 실행하여 mqTrans 기능을 확보하십시오. 자세한 설명은 섹션 2에서 찾을 수 있습니다. (c) mqTrans 기능을 선택합니다. 이 절에서는 mqTrans 기능 평가에 초점을 맞춥니다. 자세한 내용은 섹션 3을 참조하십시오. <a href="https://www.jove.com/files/ftp_upload/66030/66030fig02large.jpg" target="_blank">이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.</a>
<img alt="Figure 3" class="xfigimg" src="/files/ftp_upload/66030/66030fig03.jpg"> 그림 3: Python을 위한 환경을 준비합니다. (A) healthmodel을 만드는 명령입니다. (B) VE 생성 과정에서 y 를 입력합니다. (C) VE를 활성화하기 위한 가장 일반적인 명령입니다. (D) 토치 설치 명령 1.13.1. (E) torch-geometric 패키지에 대한 추가 라이브러리를 설치합니다. (F) torch-geometric 패키지를 설치합니다. <a href="https://www.jove.com/files/ftp_upload/66030/66030fig03large.jpg" target="_blank">이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.</a>
<img alt="Figure 4" class="xfigimg" src="/files/ftp_upload/66030/66030fig04.jpg"> 그림 4: HealthModel을 실행하여 mqTrans 기능을 가져옵니다. (ᅡ) 코드를 다운로드합니다. (B) 데이터 파일의 예. 각 열에는 조절 인자의 모든 값이 있으며 첫 번째 항목은 유전자 ID입니다. 각 행은 지정된 샘플의 값을 제공하며 첫 번째 항목은 샘플 이름입니다. (C) 라벨 파일의 예. 첫 번째 열은 샘플 이름을 제공하고 각 샘플의 클래스 레이블은 레이블이라는 열에 제공됩니다. 레이블 열의 0 값은 이 샘플이 활성 상태임을 의미하고 1은 비활성을 의미합니다. (D) mqTrans의 출력. <a href="https://www.jove.com/files/ftp_upload/66030/66030fig04large.jpg" target="_blank">이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.</a>
<img alt="Figure 5" class="xfigimg" src="/files/ftp_upload/66030/66030fig05.jpg"> 그림 5: mqTrans 기능에 대한 기능 선택 알고리즘을 실행합니다. 기능 선택 알고리즘의 결과가 사용자에게 표시됩니다. <a href="https://www.jove.com/files/ftp_upload/66030/66030fig05large.jpg" target="_blank">이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.</a>
<img alt="Figure 6" class="xfigimg" src="/files/ftp_upload/66030/66030fig06.jpg"> 그림 6: 각 기능 선택 알고리즘의 최대 테스트 세트 정확도. 가로축은 기능 선택 알고리즘을 나열하고 세로축은 정확도 값을 제공합니다. 히스토그램은 세 가지 설정, 즉 mqTrans, mRNA, mRNA+mqTrans의 실험 데이터를 보여줍니다. <a href="https://www.jove.com/files/ftp_upload/66030/66030fig06large.jpg" target="_blank">이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.</a>
<img alt="Figure 7" class="xfigimg" src="/files/ftp_upload/66030/66030fig07.jpg"> 그림 7: mqTrans 보기에서 p-값이 가장 작은 상위 50개의 어두운 바이오마커. "Dark Biomarker" 열은 다크 바이오마커 이름을 제공합니다. "mRNA.P" 및 "mqTrans.P" 열은 표현형과 대조군 간의 통계적 t-검정 p-값입니다. p-값의 배경색은 p-값 1.00(파란색)과 0.00(빨간색) 사이에 있고 흰색은 p-값 = 0.05를 나타냅니다. <a href="https://www.jove.com/files/ftp_upload/66030/66030fig07large.jpg" target="_blank">이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.</a>
<img alt="Figure 8" class="xfigimg" src="/files/ftp_upload/66030/66030fig08.jpg"> 그림 8: TCGA(The Cancer Genome Atlas)에 있는 26개 암의 세부 정보를 서로 다른 단계에서 보여줍니다. "Cohort" 및 "Disease Tissue" 열은 각 데이터 세트에 대한 환자 그룹 및 질병이 있는 조직을 설명합니다. 마지막 4개의 열은 각각 발달 단계 I, II, III 및 IV의 샘플 수를 제공합니다. <a href="https://www.jove.com/files/ftp_upload/66030/66030fig08large.jpg" target="_blank">이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.</a>
<img alt="Figure 9" class="xfigimg" src="/files/ftp_upload/66030/66030fig09.jpg"> 그림 9: 26개 암에서 다크 바이오마커와 전통적인 바이오마커의 수. 가로축에는 26가지 암 유형이 나열되어 있습니다. 세로축은 이러한 암 유형에 대한 어두운 바이오마커와 전통적인 바이오마커의 수를 제공합니다. <a href="https://www.jove.com/files/ftp_upload/66030/66030fig09large.jpg" target="_blank">이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.</a>
보충 코딩 파일 1: HealthModel-mqTrans-v1-00.tar <a href="https://www.jove.com/files/ftp_upload/66030/HealthModel-mqTrans-v1-00.tar.gz" target="_blank">이 파일을 다운로드하려면 여기를 클릭하십시오.</a>

여기에서는 전사체 데이터를 mqTrans 보기로 변환하여 다크 바이오마커를 식별할 수 있는 프로토콜을 소개합니다. 기존의 전사체 분석에서는 차등적으로 발현되지 않지만, 이러한 바이오마커는 mqTrans 보기에서 차등 발현을 나타냅니다. 이 접근법은 전통적인 방법을 보완하는 기술 역할을 하며, 이전에 간과되었던 바이오마커를 밝혀냅니다.

Transcriptome represents the expression levels of many genes in a sample and has been widely used in biological research and clinical practice. ...

우리의 연구는 질병에 대한 유전 간 상호 작용에 중점을 둡니다. 우리는 많은 유전자들이 서로 복잡하고 얽힌 관계를 가지고 있다는 것을 발견했다. 이러한 관계를 탐구함으로써 질병 치료 및 관리에서보다 정확한 진단을 내리는 것을 목표로합니다.우리는 전통적이고 전통적인 방법이 많은 의학 문헌에 의해 뒷받침되지 않는 많은 어두운 바이오마커를 발견했습니다. 이는 접근법이 생물학자가 바이오마커 스크리닝을 위한 시간 주기를 줄이는 데 더 도움이 될 수 있음을 나타냅니다. 실험 중 가장 큰 과제는 다양한 질병 유형의 작은 샘플 크기 문제를 해결하는 것이었습니다.이러한 어려움을 극복하기 위해 가능한 한 많은 건강 샘플을 수집하여 건강 모델이라고 하는 참조 모델을 훈련했습니다.

generating-transcriptional-regulation-view-transcriptomic-features

Research

연구

JoVE Journal

Biology

생물학

RNA 시퀀싱과 생물 정보학 분석은 마우스 EMLcells의 린-CD34 + 린-CD34- 부분 집단에서 유의 및 발현 된 전사 인자를 식별하는 데 사용되었다. 이들 전사 인자는 자기 갱신 린-CD34 +를 부분적으로 분화 된 린 - CD34- 세포 사이의 스위치를 결정하는데 중요한 역할을 할 수도있다.

RNA 시퀀싱 분석하여 EML 조혈 전구체 세포의자가 재생 (self-renewal)과 분화를 조절하는 중요한 요소의 식별

여기에서는 세포 기반 리포터 유전자 분석을 3' UTR을 통해 발휘되는 전사 후 제어 메커니즘을 조절하는 화합물에 대한 화학 라이브러리를 스크리닝하는 유용한 도구로 설명합니다.

후 전사적으로 규제 유전자의 화학 변조기에 대한 높은 처리량 검사

본 연구에서는 관심 있는 전사 인자를 표적으로 하는 작은 헤어핀(sh)RNA로 세포를 transfection한 후 공개적으로 사용 가능한 데이터, 계산 리소스 및 microRNA 어레이의 높은 처리량 데이터를 사용하여 microRNA 전사체에 대한 전사 인자의 영향을 연구하는 전략을 제안합니다.

마이크로 RNA 전 사체 (Transcriptome)의 전사 인자 종속 규정을 설명

이 작업에서는 활성 인핸서를 식별하고 실험적으로 검증하는 방법에 대한 실험적 워크플로를 제공합니다.

예측 및 유전자 규제 요소의 검증은 레티노 산 유도 배아 동안 활성화 된 세포 분화 줄기

기존 알고리즘 바이오 마커 검출 데이터 집합에 대 한 하나의 솔루션을 생성합니다. 이 프로토콜의 여러 마찬가지로 효과적인 솔루션의 존재 및 고 생물 의학 연구원 제안된 과제에 대 한 데이터 세트를 조사 수 있도록 사용자 친화적인 소프트웨어를 제공. 컴퓨터 과학자 들은 팬 들은이 기능에 그들의 바이오 마커 검출 알고리즘의 정보를 제공할 수도 있습니다.

마찬가지로 효과적인 이진 분류 공연으로 여러 바이오 마커 하위 집합 선택

전사 인자의 새로운 조절자를 확인하기 위해, 우리는 이중 luciferase 기반 의 전사 기자 분석을 사용하여 배열 된 렌티 바이러스 또는 레트로 바이러스 RNAi 라이브러리를 선별하는 접근 법을 개발했습니다. 이 방법은 단일 실험에서 수백 개의 후보를 선별하는 빠르고 비교적 저렴한 방법을 제공합니다.

배열된 라이브러리의 중간 처리량 스크리닝및 이중 루시퍼라제 기반 리포터를 이용한 전사 계수 조절기 식별

IR-TEx는 아노펠레스 감비아종에서 살충제 저항 성 전사 프로파일을 탐구한다. 여기에 는 응용 프로그램 사용, 여러 전사 데이터 집합을 탐색하기 위한 수정 및 프레임워크를 사용하여 모든 플랫폼에서 생성된 모든 유기체의 전사 데이터 수집을 위한 대화형 데이터베이스를 빌드하기 위한 전체 지침이 제공됩니다.

IR-TEx: 말라리아 벡터 아노펠레스 감비아를 위해 설계된 빅 데이터 전사학을 위한 오픈 소스 데이터 통합 도구

본질적으로 무질서한 도메인은 온코겐융합 전사 인자 기능에 중요합니다. 이러한 단백질을 치료적으로 표적으로 하기 위해서는 이러한 도메인에 의해 사용되는 규제 메커니즘에 대한 보다 상세한 이해가 필요합니다. 여기서는 전사학을 사용하여 유잉 육종에서 본질적으로 무질서한 EWS 도메인의 중요한 구조적 특징을 매핑합니다.

전사 분석을 사용하여 무질서한 온코겐 전사 요인의 구조 기능 관계 매핑

어둠 / 빛으로 전환 시험은 마우스의 동작과 같은 불안을 측정하기 위하여 가장 널리 사용되는 검사 중 하나입니다. 여기, 우리가 테스트를 실시하는 방법에 대한 자세한 절차를 보여주는 동영상을 제시한다.

마우스에 대한 진한 / 라이트 전환 테스트

우리는에서 고속 광학 녹음을 분석을위한 세미 자동 광학 도면 해석 방법 (하)를...

작은 하트의 세미 - 자동 광학 도면 분석

학습하세요 mqTrans 전사체 조절 보기를 활용한 다크 바이오마커 검출 방법을 Python 가상 환경에서 구현하는 방법을. 이 연구는 소규모 데이터 세트에서 221개의 다크 바이오마커를 성공적으로 식별하여 기존 방법의 한계를 극복하고 정밀 진단을 위한 새로운 바이오마커 스크리닝 기법을 제시합니다.

작은 데이터 세트에서 예측 작업 및 다크 바이오마커 검출을 위한 전사체 특징의 전사체 조절 보기 생성

더 많은 동영상 탐색