Distinct genomic traits of acral and mucosal melanomas revealed by targeted mutational profiling

The incidence of melanoma is rising globally including China. Comparing to Caucasians, the incidence of non‐cutaneous melanomas is significantly higher in Chinese. Herein, we performed genomic profiling of 89 Chinese surgically resected primary melanomas, including acral (n = 54), cutaneous (n = 22), and mucosal (n = 13), by hybrid capture‐based next‐generation sequencing. We show that mucosal melanomas tended to harbor more pathogenic mutations than other types of melanoma, though the biological significance of this finding remains uncertain. Chromosomal arm‐level alterations including 6q, 9p, and 10p/q loss were highly recurrent in all subtypes, but mucosal melanoma was significantly associated with increased genomic instability. Importantly, 7p gain significantly correlated with unfavorable clinical outcomes in non‐cutaneous melanomas, representing an intriguing prognostic biomarker of those subtypes. Furthermore, focal amplification of 4q12 (KIT, KDR, and PDGFRα) and RAD51 deletion were more abundant in mucosal melanoma, while NOTCH2 amplification was enriched in acral melanoma. Additionally, cutaneous melanomas had higher mutation load than acral melanomas, while mucosal melanomas did not differ from other subtypes in mutation burden. Together, our data revealed important features of acral and mucosal melanomas in Chinese including distinctive driver mutation pattern and increased genomic instability. These findings highlight the possibilities of combination therapies in the clinical management of melanoma.     

ERK-regulated differential expression of the Mitf 6a/b splicing isoforms in melanoma

The master regulator of the melanocyte lineage Mitf is intimately involved in development as well as melanoma, controlling cell survival, differentiation, proliferation and metastasis/migration. Consistent with its central role, Mitf expression and Mitf post-translational modifications are tightly regulated. An additional potential level of regulation is afforded by differential splicing of Mitf exon-6 leading to the generation of two isoforms that differ by the presence of six amino-acids in the Mitf (+) isoform and which have differential effects on cell cycle progression. However, whether the ratio of the two isoforms is regulated and whether their expression correlates with melanoma progression is not known. Here, we show that the differential expression of the Mitf 6a/b isoforms is dependent on the MAPKinase signalling, being linked to the activation of MEK1-ERK2, but not to N-RAS/B-RAF mutation status. In addition, quantification of Mitf 6a/b splicing forms in 86 melanoma samples revealed substantially increased levels of the Mitf (−) form in a subset of metastatic melanomas. The results suggest that differential expression of the Mitf 6a/b isoforms may represent an additional mechanism for regulating Mitf function and melanoma biology.     

The Thr300Ala variant of ATG16L1 is associated with decreased risk of brain metastasis in patients with non-small cell lung cancer

Non-small cell lung cancer (NSCLC) often metastasizes to the brain, but identifying which patients will develop brain metastases (BM) is difficult. Macroautophagy/autophagy is critical for cancer initiation and progression. We hypothesized that genetic variants of autophagy-related genes may affect brain metastases (BM) in NSCLC patients. We genotyped 16 single nucleotide polymorphisms (SNPs) in 7 autophagy-related (ATG) genes (ATG3, ATG5, ATG7, ATG10, ATG12, ATG16L1, and MAP1LC3/LC3) by using DNA from blood samples of 323 NSCLC patients. Further, we evaluated the potential associations of these genes with subsequent BM development. Lung cancer cell lines stably transfected with ATG16L1: rs2241880 (T300A) were established. Mouse models of brain metastasis were developed using cells transfected with ATG16L1–300T or ATG16L1–300A. ATG10: rs10036653 and ATG16L1: rs2241880 were significantly associated with a decreased risk of BM (respective hazard ratios [HRs]=0.596, 95% confidence interval [CI] 0.398–0.894, P = 0.012; and HR = 0. 655, 95% CI 0.438–0.978, P = 0.039, respectively). ATG12: rs26532 was significantly associated with an increased risk of BM (HR=1.644, 95% CI 1.049–2.576, P = 0.030). Invasion and migration assays indicated that transfection with ATG16L1–300T (vs. 300A) stimulated the migration of A549 cells. An in vivo metastasis assay revealed that transfection with ATG16L1–300T (vs. 300A) significantly increased brain metastasis. Our results indicate that genetic variations in autophagy-related genes can predict BM and that genome analysis would facilitate stratification of patients for BM prevention trials.     

Autophagy suppresses melanoma tumorigenesis by inducing senescence

Whether and how autophagy is involved in tumorigenesis is poorly understood. We approached this question by investigating a relatively large cohort of patients with mostly early primary melanoma for their expression of 2 markers for autophagy, the protein ATG5 (autophagy-related 5) and MAP1LC3B/LC3 (microtubule-associated protein 1 light chain 3B). Surprisingly, we discovered that both ATG5 and LC3 levels are decreased in patients with melanomas as compared with those with benign nevi. We wondered why reduced autophagy should facilitate early tumor development. Using an in vitro model of melanoma tumorigenesis, in which a mutated oncogene, BRAF (v-raf murine sarcoma viral oncogene homolog B), had been introduced into normal human melanocytes, we were able to show that downregulation of ATG5 promoted the proliferation of melanocytes because it facilitated bypassing oncogene-induced senescence (OIS). Our work supports previous reports that had argued that autophagy actually suppresses tumorigenesis and explains the possible mechanism. Furthermore, our findings suggest that the status of ATG5 and autophagy could serve as a diagnostic marker for distinguishing benign from malignant tumors of melanocytes.     

Commentary: Autophagy suppresses melanoma tumorigenesis by inducing senescence

Whether and how autophagy is involved in tumorigenesis is poorly understood. We approached this question by investigating a relatively large cohort of patients with mostly early primary melanoma for their expression of 2 markers for autophagy, the protein ATG5 (autophagy-related 5) and MAP1LC3B/LC3 (microtubule-associated protein 1 light chain 3B). Surprisingly, we discovered that both ATG5 and LC3 levels are decreased in patients with melanomas as compared with those with benign nevi. We wondered why reduced autophagy should facilitate early tumor development. Using an in vitro model of melanoma tumorigenesis, in which a mutated oncogene, BRAF (v-raf murine sarcoma viral oncogene homolog B), had been introduced into normal human melanocytes, we were able to show that downregulation of ATG5 promoted the proliferation of melanocytes because it facilitated bypassing oncogene-induced senescence (OIS). Our work supports previous reports that had argued that autophagy actually suppresses tumorigenesis and explains the possible mechanism. Furthermore, our findings suggest that the status of ATG5 and autophagy could serve as a diagnostic marker for distinguishing benign from malignant tumors of melanocytes.     

Novel primate-specific genes, RMEL 1, 2 and 3, with highly restricted expression in melanoma, assessed by new data mining tool

Melanoma is a highly aggressive and therapy resistant tumor for which the identification of specific markers and therapeutic targets is highly desirable. We describe here the development and use of a bioinformatic pipeline tool, made publicly available under the name of EST2TSE, for the in silico detection of candidate genes with tissue-specific expression. Using this tool we mined the human EST (Expressed Sequence Tag) database for sequences derived exclusively from melanoma. We found 29 UniGene clusters of multiple ESTs with the potential to predict novel genes with melanoma-specific expression. Using a diverse panel of human tissues and cell lines, we validated the expression of a subset of three previously uncharacterized genes (clusters Hs.295012, Hs.518391, and Hs.559350) to be highly restricted to melanoma/melanocytes and named them RMEL1, 2 and 3, respectively. Expression analysis in nevi, primary melanomas, and metastatic melanomas revealed RMEL1 as a novel melanocytic lineage-specific gene up-regulated during melanoma development. RMEL2 expression was restricted to melanoma tissues and glioblastoma. RMEL3 showed strong up-regulation in nevi and was lost in metastatic tumors. Interestingly, we found correlations of RMEL2 and RMEL3 expression with improved patient outcome, suggesting tumor and/or metastasis suppressor functions for these genes. The three genes are composed of multiple exons and map to 2q12.2, 1q25.3, and 5q11.2, respectively. They are well conserved throughout primates, but not other genomes, and were predicted as having no coding potential, although primate-conserved and human-specific short ORFs could be found. Hairpin RNA secondary structures were also predicted. Concluding, this work offers new melanoma-specific genes for future validation as prognostic markers or as targets for the development of therapeutic strategies to treat melanoma.     

Biological processes, properties and molecular wiring diagrams of candidate low-penetrance breast cancer susceptibility genes

Background

The process of malignant transformation, progression and metastasis of melanoma is poorly understood. Gene expression profiling of human cancer has allowed for a unique insight into the genes that are involved in these processes. Thus, we have attempted to utilize this approach through the analysis of a series of primary, non-metastatic cutaneous tumors and metastatic melanoma samples.

Methods

We have utilized gene microarray analysis and a variety of molecular techniques to compare 40 metastatic melanoma (MM) samples, composed of 22 bulky, macroscopic (replaced) lymph node metastases, 16 subcutaneous and 2 distant metastases (adrenal and brain), to 42 primary cutaneous cancers, comprised of 16 melanoma, 11 squamous cell, 15 basal cell skin cancers. A Human Genome U133 Plus 2.0 array from Affymetrix, Inc. was utilized for each sample. A variety of statistical software, including the Affymetrix MAS 5.0 analysis software, was utilized to compare primary cancers to metastatic melanomas. Separate analyses were performed to directly compare only primary melanoma to metastatic melanoma samples. The expression levels of putative oncogenes and tumor suppressor genes were analyzed by semi- and real-time quantitative RT-PCR (qPCR) and Western blot analysis was performed on select genes.

Results

We find that primary basal cell carcinomas, squamous cell carcinomas and thin melanomas express dramatically higher levels of many genes, including SPRR1A/B, KRT16/17, CD24, LOR, GATA3, MUC15, and TMPRSS4, than metastatic melanoma. In contrast, the metastatic melanomas express higher levels of genes such as MAGE, GPR19, BCL2A1, MMP14, SOX5, BUB1, RGS20, and more. The transition from non-metastatic expression levels to metastatic expression levels occurs as melanoma tumors thicken. We further evaluated primary melanomas of varying Breslow's tumor thickness to determine that the transition in expression occurs at different thicknesses for different genes suggesting that the "transition zone" represents a critical time for the emergence of the metastatic phenotype. Several putative tumor oncogenes (SPP-1, MITF, CITED-1, GDF-15, c-Met, HOX loci) and suppressor genes (PITX-1, CST-6, PDGFRL, DSC-3, POU2F3, CLCA2, ST7L), were identified and validated by quantitative PCR as changing expression during this transition period. These are strong candidates for genes involved in the progression or suppression of the metastatic phenotype.

Conclusion

The gene expression profiling of primary, non-metastatic cutaneous tumors and metastatic melanoma has resulted in the identification of several genes that may be centrally involved in the progression and metastatic potential of melanoma. This has very important implications as we continue to develop an improved understanding of the metastatic process, allowing us to identify specific genes for prognostic markers and possibly for targeted therapeutic approaches.     

Characterization of unusual families of ATG8-like proteins and ATG12 in the protozoan parasite Leishmania major

Leishmania major possesses, apparently uniquely, four families of ATG8-like genes, designated ATG8, ATG8A, ATG8B and ATG8C, and 25 genes in total.  L. major ATG8 and examples from the ATG8A, ATG8B and ATG8C families are able to complement a Saccharomyces cerevisiae ATG8-deficient strain, indicating functional conservation. Whereas ATG8 has been shown to form putative autophagosomes during differentiation and starvation of L. major, ATG8A primarily form puncta in response to starvation - suggesting a role for ATG8A in starvation-induced autophagy. Recombinant ATG8A was processed at the scissile glycine by recombinant ATG4.2 but not ATG4.1 cysteine peptidases of L. major and, consistent with this, ATG4.2-deficient L. major mutants were unable to process ATG8A and were less able to withstand starvation than wild type cells. GFP-ATG8-containing puncta were less abundant in ATG4.2 over-expression lines, in which unlipidated ATG8 predominated, which is consistent with ATG4.2 being an ATG8-deconjugating enzyme as well as an ATG8A-processing enzyme. In contrast, recombinant ATG8, ATG8B and ATG8C were all processed by ATG4.1, but not by ATG4.2. ATG8B and ATG8C both have a distinct subcellular location close to the flagellar pocket, but the occurrence of the GFP-labelled puncta suggest that they do not have a role in autophagy. L. major genes encoding possible ATG5, ATG10 and ATG12 homologues were found to complement their respective S. cerevisiae mutants, and ATG12 localised in part to ATG8-containing puncta, suggestive of a functional ATG5-ATG12 conjugation pathway in the parasite. L. major ATG12 is unusual as it requires C-terminal processing by an as yet unidentified peptidase.     

ATG8 lipidation and ATG8‐mediated autophagy in Arabidopsis require ATG12 expressed from the differentially controlled ATG12A AND ATG12B loci

Autophagic recycling of intracellular plant constituents is maintained at a basal level under normal growth conditions but can be induced in response to nutritional demand, biotic stress, and senescence. One route requires the ubiquitin‐fold proteins Autophagy‐related (ATG)‐8 and ATG12, which become attached to the lipid phosphatidylethanolamine (PE) and the ATG5 protein, respectively, during formation of the engulfing vesicle and delivery of its cargo to the vacuole for breakdown. Here, we genetically analyzed the conjugation machinery required for ATG8/12 modification in Arabidopsis thaliana with a focus on the two loci encoding ATG12. Whereas single atg12a and atg12b mutants lack phenotypic consequences, atg12a atg12b double mutants senesce prematurely, are hypersensitive to nitrogen and fixed carbon starvation, and fail to accumulate autophagic bodies in the vacuole. By combining mutants eliminating ATG12a/b, ATG5, or the ATG10 E2 required for their condensation with a method that unequivocally detects the ATG8‐PE adduct, we also show that ATG8 lipidation requires the ATG12–ATG5 conjugate. Unlike ATG8, ATG12 does not associate with autophagic bodies, implying that its role(s) during autophagy is restricted to events before the vacuolar deposition of vesicles. The expression patterns of the ATG12a and ATG12b genes and the effects of single atg12a and atg12b mutants on forming the ATG12–ATG5 conjugate reveal that the ATG12b locus is more important during basal autophagy while the ATG12a locus is more important during induced autophagy. Taken together, we conclude that the formation of the ATG12–ATG5 adduct is essential for ATG8‐mediated autophagy in plants by promoting ATG8 lipidation.     

A novel Atg5-shRNA mouse model enables temporal control of Autophagy in vivo

Macroautophagy/autophagy is an evolutionarily conserved catabolic pathway whose modulation has been linked to diverse disease states, including age-associated disorders. Conventional and conditional whole-body knockout mouse models of key autophagy genes display perinatal death and lethal neurotoxicity, respectively, limiting their applications for in vivo studies. Here, we have developed an inducible shRNA mouse model targeting Atg5, allowing us to dynamically inhibit autophagy in vivo, termed ATG5i mice. The lack of brain-associated shRNA expression in this model circumvents the lethal phenotypes associated with complete autophagy knockouts. We show that ATG5i mice recapitulate many of the previously described phenotypes of tissue-specific knockouts. While restoration of autophagy in the liver rescues hepatomegaly and other pathologies associated with autophagy deficiency, this coincides with the development of hepatic fibrosis. These results highlight the need to consider the potential side effects of systemic anti-autophagy therapies.     

miR-216b enhances the efficacy of vemurafenib by targeting Beclin-1, UVRAG and ATG5 in melanoma

Autophagy maintains cells survival in many stressful conditions including starvation, growth factor deprivation and misfolded protein accumulation. Additionally, autophagic survival mechanisms are used by transformed tumor cells to inhibit cell death, limit drug effectiveness and possibly generate drug resistance. However, the mechanism of how cells utilize autophagy during drug resistance is not fully understood. Here, we demonstrate that miR-216b plays an important role in alleviating drug resistance by regulating autophagy in melanoma. We show that miR-216b attenuates autophagy by directly targeting three key autophagy genes Beclin-1, UVRAG and ATG5. Overexpression of these genes from miRNA immune cDNA constructs rescue autophagic activity in the presence of miR-216b. Antagomir-mediated inactivation of endogenous miR-216b led to an increase of Beclin-1, UVRAG, ATG5, and subsequent autophagic activity. More importantly, we have discovered that BRAF(V600E) inhibitor vemurafenib suppresses miR-216b activity, which in turn activates autophagy to generate drug resistance in both BRAFi-sensitive and -resistant cells. Strikingly, ectopic expression of miR-216b increases the efficacy of vemurafenib both in vitro and in vivo. Taken together, these data indicate that miR-216b regulates autophagy by suppressing three key autophagy genes, and enhances the antitumor activity of vemurafenib in BRAF(V600E) melanoma cells.     

Autophagy inhibition suppresses pulmonary metastasis of HCC in mice via impairing anoikis resistance and colonization of HCC cells

Metastasis is one of the main causes of poor prognosis for hepatocellular carcinoma (HCC), which has been linked to cell-death resistance. Autophagy is an important survival mechanism under conditions of cell stress. We hypothesized that autophagy may play a role in HCC metastasis due to its prosurvival effect. Highly metastatic HCC cell lines with stable autophagy inhibition were established via lentivirus-mediated silencing of BECN1 and ATG5 genes. Mouse models of pulmonary metastasis were then developed using the cells with or without autophagy inhibition. The analysis of lung metastasis by histopathological examination and small animal imaging showed that autophagy inhibition significantly decreased the incidence of pulmonary metastases in vivo. Further invasion, migration, detachment, lung colonization, and epithelial-mesenchymal transition (EMT) assays indicated that autophagy inhibition did not affect cell invasiveness, migration or EMT but attenuated the anoikis-resistance and lung colonization of HCC cells. Investigation of the molecular mechanisms underlying showed that the autophagy-inhibition-mediated anoikis-resistance attenuation was associated with the regulation of apoptotic signaling. As autophagy inhibition was shown to be able to suppress HCC metastasis, an autophagy-based HCC tissue-specific target therapy system (AFP-Cre/LoxP-shRNA) was constructed. In vitro and in vivo analyses showed that the system was able to efficiently inhibit autophagy of HCC cells and tissue in a tissue-specific manner. Further in vivo metastasis assay showed that intratumoral administration of the system could significantly suppress lung metastasis. Together, our findings suggest that autophagy may be involved in HCC metastasis through facilitating anoikis resistance and lung colonization of HCC cells. Autophagy-based HCC tissue-specific target therapy may be a new strategy for the management of HCC metastasis.     

Molecular pathogenesis of malignant melanoma: a different perspective from the studies of melanocytic nevus and acral melanoma

The Clark model for melanoma progression emphasizes a series of histopathological changes beginning from benign melanocytic nevus to melanoma via dysplastic nevus. Several models of the genetic basis of melanoma development and progression are based on this Clark’s multi-step model, and predict that the acquisition of a BRAF mutation can be a founder event in melanocytic neoplasia. However, our recent investigations have challenged this view, showing the polyclonality of BRAF mutations in melanocytic nevi. Furthermore, it is suggested that many melanomas, including acral and mucosal melanomas, arise de novo, not from melanocytic nevus. While mutations of the BRAF gene are frequent in melanomas on non-chronic sun damaged skin which are prevalent in Caucasians, acral and mucosal melanomas harbor mutations of the KIT gene as well as the amplifications of cyclin D1 or cyclin-dependent kinase 4 gene. Amplifications of the cyclin D1 gene are detected in normal-looking ‘field melanocytes’, which represent a latent progression phase of acral melanoma that precedes the stage of atypical melanocyte proliferation in the epidermis. Based on these observations, we propose an alternative genetic progression model for melanoma.