Gene master regulators of papillary and anaplastic thyroid cancers

We hypothesize that distinct cell phenotypes are governed by different sets of gene master regulators (GMRs) whose strongly protected (by the homeostatic mechanisms) abundance modulates most cell processes by coordinating the expression of numerous genes from the corresponding functional pathways. Gene Commanding Height (GCH), a composite measure of gene expression control and coordination, is introduced to establish the gene hierarchy in each phenotype. If the hypothesis is true, than one can selectively destroy cancer nodules from a heterogeneous tissue by altering the expression of genes whose GCHs are high in cancer but low in normal cell phenotype. Here, we test the hypothesis and show its utility for the thyroid cancer (TC) gene therapy. First, we prove that malignant and cancer free surrounding areas of a surgically removed papillary TC (PTC) tumor are governed by different GMRs. Second, we show that stable transfection of a gene induces larger transcriptomic alterations in the cells where it has higher GCH than in other cells. For this, we profiled the transcriptomes of the papillary BCPAP and anaplastic 8505C TC cell lines before and after stable transfection with NEMP1, DDX19B, PANK2 or UBALD1. The four genes were selected to have similar expression levels but significantly different GCH scores in the two cell lines before transfection. Indeed, each of the four genes triggered larger alterations in the cells where they had larger GCH. Our results prove the feasibility of a personalized gene therapy approach that selectively targets the cancer cells from a tissue.