Classification of tumor types with gene expression data

The Cancer Genome Atlas (TCGA) has generated comprehensive molecular profiles. Researchers at the National Institute of Environmental Health Sciences aimed to identify a set of genes whose expression patterns can distinguish diverse tumor types. Those features may serve as biomarkers for tumor diagnosis and drug development.

Using RNA-seq expression data, the researchers undertook a pan-cancer classification of 9,096 TCGA tumor samples representing 31 tumor types. They randomly assigned 75% of samples into training and 25% into testing, proportionally allocating samples from each tumor type.

The researchers could correctly classify more than 90% of the test set samples. Accuracies were high for all but three of the 31 tumor types, in particular, for READ (rectum adenocarcinoma) which was largely indistinguishable from COAD (colon adenocarcinoma). They also carried out pan-cancer classification, separately for males and females, on 23 sex non-specific tumor types (those unrelated to reproductive organs). Results from these gender-specific analyses largely recapitulated results when gender was ignored. Remarkably, more than 80% of the 100 most discriminative genes selected from each gender separately overlapped. Genes that were differentially expressed between genders included BNC1, FAT2, FOXA1, and HOXA11. FOXA1 has been shown to play a role for sexual dimorphism in liver cancer. The differentially discriminative genes identified might be important for the gender differences in tumor incidence and survival.

Heatmap representation of the expression patterns of the top 50 genes across all 9096 samples

rna-seq

Each row (gene) was centered by the median expression value across all samples. A hierarchical clustering analysis was carried out for both samples and genes using the Euclidean distance as the similarity metric

These researchers were able to identify many sets of 20 genes that could correctly classify more than 90% of the samples from 31 different tumor types using TCGA RNA-seq data. This accuracy is remarkable given the number of the tumor types and the total number of samples involved. They achieved similar results when they analyzed 23 non-sex-specific tumor types separately for males and females. They regard the frequency with which a gene appeared in those sets as measuring its importance for tumor classification. One third of the 50 most frequently appearing genes were pseudogenes; the degree of enrichment may be indicative of their importance in tumor classification. Lastly, the researchers identified a few genes that might play a role in sexual dimorphism in certain cancers.

Li Y, Kang K, Krahn JM, Croutwater N, Lee K, Umbach DM, Li L. (2017) A comprehensive genomic pan-cancer classification using The Cancer Genome Atlas gene expression data. BMC Genomics 18(1):508. [article]

Leave a Reply

Your email address will not be published. Required fields are marked *

*

Time limit is exhausted. Please reload CAPTCHA.