2024年3月1日
Here, we introduce a protocol for converting transcriptomic data into a mqTrans view, enabling the identification of dark biomarkers. While not differentially expressed in conventional transcriptomic analyses, these biomarkers exhibit differential expression in the mqTrans view. The approach serves as a complementary technique to traditional methods, unveiling previously overlooked biomarkers.
Our research focuses on the intergenetic interactions of disease. We found that many genes have complex intertwined relationships with each other. By exploring these relationships, we aim in making more precise diagnosis in disease treatment and management.
We found that many dark biomarkers with traditional combinational methods ignored, are supported by many medical literatures. This indicates that approach can further assist biologists in reducing time cycle for biomarker screening. The biggest challenge during the experiment was addressing the issue of small sample sizes of different disease types.
To overcome this difficulty, we collected as many healthy samples as possible to train a reference model, which we called Health Model. To begin, create a new virtual environment named Health Model with Python version 3.7. in the Slurm Cluster supercomputing platform, execute the module load anaconda command.
Once the command is executed, a confirmation prompt appears on the screen. Enter Y to proceed and wait for the process to complete. Then activate the virtual environment following platform specific instructions.
Next, run the command to install PyTorch 1.13.1. Install additional packages for torch geometric, such as torch underscore scatter, torch underscore sparse, torch underscore cluster, and torch underscore spline underscore convulsion, following the installation guidelines. Then install the torch geometric package version 2.2.0.
Download the code and the pre-trained Health Model from the Health Informatics Lab website. Decompress the file to a desired path. Then change the working directory in the command line to the Health Model MQ trans folder.
Execute the command to generate MQ trans features and obtain the outputs. The MQ trans features will be generated as output MQ targets CSV, and the label file will be received as output label CSV. Additionally, original expression values of the MRNA genes will be extracted as file output test targets CSV.
Next, use the feature selection algorithm for selecting MQ trans features. If selecting MQ trans features or original features without combining them, set combine to false. Select 800 original features and split the dataset into 0.8 to 0.2 for training and testing.
To combine MQ trans features with the original expression values for feature selection, set combine to true. Dark biomarkers with differential MQ trans values, but undifferential mRNA expression were identified. Among 3, 062 features, 221 dark biomarkers were detected.
A general scarcity of dark biomarkers was observed in comparison to traditional biomarkers across most cancer types except BRCA, MESO, and TGCT.
本研究探讨了与疾病相关的基因之间复杂的相互作用,重点在于识别传统方法常忽略的“暗”生物标志物。所提出的mqTrans视角有助于深入理解这些生物标志物,其表达水平相较于传统转录组学分析表现出差异性表达。
将多任务图注意力网络与转录组学数据整合,能够检测出传统差异表达分析所忽略的调控特征和隐性生物标志物。该方法通过揭示隐藏的调控信号,提高了早期发现阶段的预测置信度,有助于更稳健地进行靶点验证以及基于风险调整的研发组合决策。HealthModel 框架解决了小规模数据集的挑战,拓展了转录组谱在生物制药研发管线中的应用价值。
HealthModel 和 mqTrans 工作流程通过从转录组数据生成调控特征视图,连接了早期发现、靶点验证和临床前研究。