Tumor-specific gene expression patterns with gene expression profiles

Gene expression profiles of 14 common tumors and their counterpart normal tissues were analyzed with machine learning methods to address the problem of selection of tumor-specific genes and analysis of their differential expressions in tumor tissues. First, a variation of the Relief algorithm, “RFE_Relief algorithm” was proposed to learn the relations between genes and tissue types. Then, a support vector machine was employed to find the gene subset with the best classification performance for distinguishing cancerous tissues and their counterparts. After tissue-specific genes were removed, cross validation experiments were employed to demonstrate the common deregulated expressions of the selected gene in tumor tissues. The results indicate the existence of a specific expression fingerprint of these genes that is shared in different tumor tissues, and the hallmarks of the expression patterns of these genes in cancerous tissues are summarized at the end of this paper.     

Tumor-specific gene expression patterns with gene expression profiles

Gene expression profiles of 14 common tumors and their counterpart normal tissues were analyzed with machine learning methods to address the problem of selection of tumor-specific genes and analysis of their differential expressions in tumor tissues. First, a variation of the Relief algorithm, "RFE_Relief algorithm" was proposed to learn the relations between genes and tissue types. Then, a support vector machine was employed to find the gene subset with the best classification performance for distinguishing cancerous tissues and their counterparts. After tissue-specific genes were removed, cross validation experiments were employed to demonstrate the common deregulated expressions of the selected gene in tumor tissues. The results indicate the existence of a specific expression fingerprint of these genes that is shared in different tumor tissues, and the hallmarks of the expression patterns of these genes in cancerous tissues are summarized at the end of this paper.     

Aberrant CpG methylation of the TFAP2A gene constitutes a mechanism for loss of TFAP2A expression in human metastatic melanoma

Metastatic melanoma is a deadly treatment-resistant form of skin cancer whose global incidence is on the rise. During melanocyte transformation and melanoma progression the expression profile of many genes changes. Among these, a gene implicated in several steps of melanocyte development, TFAP2A, is frequently silenced; however, the molecular mechanism of TFAP2A silencing in human melanoma remains unknown. In this study, we measured TFAP2A mRNA expression in primary human melanocytes compared to 11 human melanoma samples by quantitative real-time RT-PCR. In addition, we assessed CpG DNA methylation of the TFAP2A promoter in these samples using bisulfite sequencing. Compared to primary melanocytes, which showed high TFAP2A mRNA expression and no promoter methylation, human melanoma samples showed decreased TFAP2A mRNA expression and increased promoter methylation. We further show that increased CpG methylation correlates with decreased TFAP2A mRNA expression. Using The Cancer Genome Atlas, we further identified TFAP2A as a gene displaying among the most decreased expression in stage 4 melanomas vs. non-stage 4 melanomas, and whose CpG methylation was frequently associated with lack of mRNA expression. Based on our data, we conclude that TFAP2A expression in human melanomas can be silenced by aberrant CpG methylation of the TFAP2A promoter. We have identified aberrant CpG DNA methylation as an epigenetic mark associated with TFAP2A silencing in human melanoma that could have significant implications for the therapy of human melanoma using epigenetic modifying drugs.     

Aberrant CpG methylation of the TFAP2A gene constitutes a mechanism for loss of TFAP2A expression in human metastatic melanoma

Metastatic melanoma is a deadly treatment-resistant form of skin cancer whose global incidence is on the rise. During melanocyte transformation and melanoma progression the expression profile of many genes changes. Among these, a gene implicated in several steps of melanocyte development, TFAP2A, is frequently silenced; however, the molecular mechanism of TFAP2A silencing in human melanoma remains unknown. In this study, we measured TFAP2A mRNA expression in primary human melanocytes compared to 11 human melanoma samples by quantitative real-time RT-PCR. In addition, we assessed CpG DNA methylation of the TFAP2A promoter in these samples using bisulfite sequencing. Compared to primary melanocytes, which showed high TFAP2A mRNA expression and no promoter methylation, human melanoma samples showed decreased TFAP2A mRNA expression and increased promoter methylation. We further show that increased CpG methylation correlates with decreased TFAP2A mRNA expression. Using The Cancer Genome Atlas, we further identified TFAP2A as a gene displaying among the most decreased expression in stage 4 melanomas vs. non-stage 4 melanomas, and whose CpG methylation was frequently associated with lack of mRNA expression. Based on our data, we conclude that TFAP2A expression in human melanomas can be silenced by aberrant CpG methylation of the TFAP2A promoter. We have identified aberrant CpG DNA methylation as an epigenetic mark associated with TFAP2A silencing in human melanoma that could have significant implications for the therapy of human melanoma using epigenetic modifying drugs.     

Decreased expression of Apaf‐1 with progression of melanoma

Defects in apoptotic system may contribute in the pathogenesis and resistance of malignant melanoma cells to chemotherapy. Apoptotic protease‐activating factor‐1 (Apaf‐1) is a cell death effector that acts with cytochrome c and caspase‐9 to mediate apoptosis. Recently it was shown that metastatic melanomas often lose Apaf‐1 and are concomitantly resistant to apoptosis. It is not known, however, whether Apaf‐1 protein is lost during melanoma progression from localized to metastatic tumor. To this end, we evaluated Apaf‐1 protein expression by immunohistochemistry in 10 cases of human nevi, 11 melanomas in situ, 26 primary melanomas and 15 metastases. Significant decreases in Apaf‐1 expression was observed when comparing nevi and melanomas (chi‐square = 33.719; P < 0.0001). Moreover, primary melanomas with greater tumor thickness showed lesser expression of Apaf‐1 (chi‐square = 16.182; P < 0.003). Intriguingly, we were unable to detect Apaf‐1 expression in lesions of metastatic melanomas. These data demonstrated that there is an inverse correlation between Apaf‐1 expression and pathologic stage of melanoma. This suggests that the decreased expression of Apaf‐1 seen in correlation with melanoma progression renders melanoma more resistant to chemotherapy.     

Systematic classification of melanoma cells by phenotype-specific gene expression mapping

There is growing evidence that the metastatic spread of melanoma is driven not by a linear increase in tumorigenic aggressiveness, but rather by switching back and forth between two different phenotypes of metastatic potential. In vitro these phenotypes are respectively defined by the characteristics of strong proliferation/weak invasiveness and weak proliferation/strong invasiveness. Melanoma cell phenotype is tightly linked to gene expression. Taking advantage of this, we have developed a gene expression-based tool for predicting phenotype called Heuristic Online Phenotype Prediction. We demonstrate the predictive utility of this tool by comparing phenotype-specific signatures with measurements of characteristics of melanoma phenotype-specific biology in different melanoma cell lines and short-term cultures. We further show that 86% of 536 tested melanoma lines and short-term cultures are significantly associated with the phenotypes we describe. These findings reinforce the concept that a two-state system, as described by the phenotype switching model, underlies melanoma progression.     

Melanomas of unknown primary have a mutation profile consistent with cutaneous sun‐exposed melanoma

Melanoma of unknown primary (MUP) is an uncommon phenomenon whereby patients present with metastatic disease without an evident primary site. To determine their likely site of origin, we combined exome sequencing from 33 MUPs to assess the total rate of somatic mutations and degree of UV mutagenesis. An independent cohort of 91 archival MUPs was also screened for 46 hot spot mutations highly prevalent in melanoma including BRAF, NRAS, KIT, GNAQ, and GNA11. Results showed that the majority of MUPs exhibited high somatic mutation rates, high ratios of C>T/G>A transitions, and a high rate of BRAF (45 of 101, 45%) and NRAS (32 of 101, 32%) mutations, collectively indicating a mutation profile consistent with cutaneous sun‐exposed melanomas. These data suggest that a significant proportion of MUPs arise from regressed or unrecognized primary cutaneous melanomas or arise de novo in lymph nodes from nevus cells that have migrated from the skin.     

Nonrandom RNAseq gene expression associated with RNAlater and flash freezing storage methods

RNA sequencing is a popular next‐generation sequencing technique for assaying genome‐wide gene expression profiles. Nonetheless, it is susceptible to biases that are introduced by sample handling prior gene expression measurements. Two of the most common methods for preserving samples in both field‐based and laboratory conditions are submersion in RNAlater and flash freezing in liquid nitrogen. Flash freezing in liquid nitrogen can be impractical, particularly for field collections. RNAlater is a solution for stabilizing tissue for longer‐term storage as it rapidly permeates tissue to protect cellular RNA. In this study, we assessed genome‐wide expression patterns in 30‐day‐old fry collected from the same brood at the same time point that were flash‐frozen in liquid nitrogen and stored at ?80°C or submerged and stored in RNAlater at room temperature, simulating conditions of fieldwork. We show that sample storage is a significant factor influencing observed differential gene expression. In particular, genes with elevated GC content exhibit higher observed expression levels in liquid nitrogen flash‐freezing relative to RNAlater storage. Further, genes with higher expression in RNAlater relative to liquid nitrogen experience disproportionate enrichment for functional categories, many of which are involved in RNA processing. This suggests that RNAlater may elicit a physiological response that has the potential to bias biological interpretations of expression studies. The biases introduced to observed gene expression arising from mimicking many field‐based studies are substantial and should not be ignored.     

人睾丸与脑基因表达谱的相似性

人类的物种形成与进化问题一直是研究的一个焦点。近年来,对于人和灵长类以及果蝇等其他一些动物多种组织基因表达谱的研究表明,在人的进化过程中脑基因表达的改变最为显著,并且脑中许多基因的表达呈显著上调。信息学分析显示,在多种组织当中,人的脑与睾丸可能存在最为相似的基因表达谱。这些结果提示睾丸可能与脑类似,也在人的物种形成和进化历程中起着重要作用。本文对人睾丸和脑基因表达谱的研究进行了回顾,并提出了该研究方向今后的一些研究设想。     

Clinicopathological characteristics associated with BRAFK601E and BRAFL597 mutations in melanoma

BRAF mutations at codons L597 and K601 occur uncommonly in melanoma. Clinical and pathological associations of these mutations were investigated in a cohort of 1119 patients with known BRAF mutation status. A BRAF mutation was identified in 435 patients; Mutations at L597 and the K601E mutation were seen in 3.4 and 3.2% of these, respectively. K601E melanomas tended to occur in male patients, a median age of 58 yr, were generally found on the trunk (64%) and uncommonly associated with chronically sun‐damaged (CSD) skin. BRAF L597 melanomas occurred in older patients (median 66 yr), but were associated with CSD skin (extremities or head and neck location – 73.3%, P = 0.001). Twenty‐three percent of patients with V600E‐ and 43% of patients with K601E‐mutant melanomas presented with nodal disease at diagnosis compared to just 14% of patients with BRAF wild‐type tumors (P = 0.001 and 0.006, respectively). Overall, these mutations represent a significant minority of BRAF mutations, but have distinct clinicopathological phenotypes and clinical behaviors.     

Immune signatures in cutaneous melanoma correlate with survival independently of immunotherapy treatment

Immune checkpoint inhibitors (ICIs) have fundamentally improved survival from advanced cutaneous melanoma. Significant efforts have been made to understand the ICI response to identify ways to further improve outcomes. One such approach has been to investigate gene expression associated with response to ICI, which has identified various immune-related mRNA signatures, including a six-gene IFN-γ signature (IFN-γ⁶), an expanded immune signature (IFN-γ¹⁸), an effector T-cell gene signature (T_eff), and a T_eff-associated and IFN-γ-associated gene signature (T_eff + IFN-γ). Given that these signatures appear to reflect expression from T cells and the level of tumour-infiltrating immune cells has been associated with survival, we hypothesised that the prognostic value of the signatures is not limited to ICI treatment and investigated if they were associated with survival also in patients who never received ICI. The signatures were not present in melanoma cell lines when compared with tumour samples, confirming that the signatures were likely derived from the samples' non-tumour (immune) components. We acquired expression and survival data from five melanoma cohorts with a wide range of disease stages, treatments and metrics for survival, and correlated the expression signatures with survival. All four signatures were significantly associated (p < .05) with survival in four of five cohorts, with hazard ratios ranging from 0.69 to 0.92. We conclude that these immune signatures' association with survival is not specific to ICI-treated patients, but present in a number of settings.