Genomic analysis reveals distinct mechanisms and functional classes of SOX10-regulated genes in melanocytes

p300 KAT Regulates SOX10 Stability and Function in Human Melanoma

SOX10 is a lineage-specific transcription factor critical for melanoma tumor growth; on the other hand, SOX10 loss-of-function drives the emergence of therapy-resistant, invasive melanoma phenotypes. A major challenge has been developing therapeutic strategies targeting SOX10's role in melanoma proliferation while preventing a concomitant increase in tumor cell invasion. In this study, we report that the lysine acetyltransferase (KAT) EP300 and SOX10 gene loci on chromosome 22 are frequently co-amplified in melanomas, including UV-associated and acral tumors. We further show that p300 KAT activity mediates SOX10 protein stability and that the p300 inhibitor A-485 downregulates SOX10 protein levels in melanoma cells via proteasome-mediated degradation. Additionally, A-485 potently inhibits proliferation of SOX10+ melanoma cells while decreasing invasion in AXLhigh/MITFlow melanoma cells through downregulation of metastasis-related genes. We conclude that the SOX10/p300 axis is critical to melanoma growth and invasion and that inhibition of p300 KAT activity through A-485 may be a worthwhile therapeutic approach for SOX10-reliant tumors.

SIGNIFICANCE

The p300 KAT inhibitor A-485 blocks SOX10-dependent proliferation and SOX10-independent invasion in hard-to-treat melanoma cells.

Gene signature reveals decreased SOX10-dependent transcripts in malignant cells from immune checkpoint inhibitor-resistant cutaneous melanomas

Evidence is mounting for cross-resistance between immune checkpoint and targeted kinase inhibitor therapies in cutaneous melanoma patients. Since the loss of the transcription factor, SOX10, causes tolerance to MAPK pathway inhibitors, we used bioinformatic techniques to determine if reduced SOX10 expression/activity is associated with immune checkpoint inhibitor resistance. We integrated SOX10 ChIP-seq, knockout RNA-seq, and knockdown ATAC-seq data from melanoma cell models to develop a robust SOX10 gene signature. We used computational methods to validate this signature as a measure of SOX10-dependent activity in independent single-cell and bulk RNA-seq SOX10 knockdown, cell line panel, and MAPK inhibitor drug-resistant datasets. Evaluation of patient single-cell RNA-seq data revealed lower levels of SOX10-dependent transcripts in immune checkpoint inhibitor-resistant tumors. Our results suggest that SOX10-deficient melanoma cells are associated with cross-resistance between targeted and immune checkpoint inhibitors and highlight the need to identify therapeutic strategies that target this subpopulation.

Recent advances in the regulation mechanism of SOX10

Neural crest (NC) is the primitive neural structure in embryonic stage, which develops from ectodermal neural plate cells and epithelial cells. When the neural fold forms into neural tube, neural crest also forms a cord like structure above the neural tube and below the ectoderm. Neural crest cells (NCC) have strong migration and proliferation abilities. A number of tissue cells differentiate from neural crest cells, such as melanocytes, central and peripheral neurons, glial cells, craniofacial cells, osteoblasts, chondrocytes and smooth muscle cells. The migration and differentiation of neural crest cells are regulated by a gene network where a variety of genes, transcriptional factors, signal pathways and growth factors are involved.

Convergent biological pathways underlying the Kallmann syndrome-linked genes Hs6st1 and Fgfr1

Kallmann syndrome (KS) is a congenital disorder characterized by idiopathic hypogonadotropic hypogonadism and olfactory dysfunction. KS is linked to variants in >34 genes, which are scattered across the human genome and show disparate biological functions. Although the genetic basis of KS is well studied, the mechanisms by which disruptions of these diverse genes cause the same outcome of KS are not fully understood. Here we show that disruptions of KS-linked genes affect the same biological processes, indicating convergent molecular mechanisms underlying KS. We carried out machine learning-based predictions and found that KS-linked mutations in heparan sulfate 6-O-sulfotransferase 1 (HS6ST1) are likely loss-of-function mutations. We next disrupted Hs6st1 and another KS-linked gene, fibroblast growth factor receptor 1 (Fgfr1), in mouse neuronal cells and measured transcriptome changes using RNA sequencing. We found that disruptions of Hs6st1 and Fgfr1 altered genes in the same biological processes, including the upregulation of genes in extracellular pathways and the downregulation of genes in chromatin pathways. Moreover, we performed genomics and bioinformatics analyses and found that Hs6st1 and Fgfr1 regulate gene transcription likely via the transcription factor Sox9/Sox10 and the chromatin regulator Chd7, which are also associated with KS. Together, our results demonstrate how different KS-linked genes work coordinately in a convergent signaling pathway to regulate the same biological processes, thus providing new insights into KS.

The DNA-Binding High-Mobility Group Box Domain of Sox Family Proteins Directly Interacts with RNA In Vitro

There is a growing body of evidence that a substantial number of protein domains identified as DNA-binding also interact with RNA to regulate biological processes. Several recent studies have revealed that the Sox2 transcription factor binds RNA through its high-mobility group box (HMGB) domain in vitro and in vivo. A high degree of conservation of this domain among members of the Sox family of transcription factors suggests that RNA-binding activity may be a general feature of these proteins. To address this hypothesis, we examined a subset of HMGB domains from human Sox family of proteins for their ability to bind both DNA and RNA in vitro. We observed selective, high-affinity interactions between Sox family HMGB domains and various model RNA elements, including a four-way junction RNA, a hairpin RNA with an internal bulge, G-quadruplex RNA, and a fragment of long noncoding RNA ES2, which is known to directly interact with Sox2. Importantly, the HMGB domains bind these RNA ligands significantly tighter than nonconsensus dsDNA and in some cases with affinities rivaling those of their consensus dsDNA sequences. These data suggest that RNA binding is a conserved feature of the Sox family of transcription factors with the potential to modulate unappreciated biological functions.

Three novel variants in SOX10 gene: Waardenburg and PCWH syndromes

Background

Waardenburg syndrome (WS) is a rare genetic disorder characterized by musculoskeletal abnormalities, deafness and hypopigmentation of hair and skin. This article’s aim is to investigate clinical and genetic characteristics of WS in three unrelated Caucasian individuals.

Case presentation

The first patient was a 25-year-old female with congenital bilateral hearing loss, bright-blue-eyes, hypopigmentation of hair and skin, megacolon, language retardation, tenosynovitis and neuromas. The second case was an infant symptomatic from birth, with dysphagia, Hirschsprung disease and neurological abnormalities. The third patient was a 14-year-old boy with congenital bilateral hearing loss and ileocolic Hirschsprung disease. In order to identify variants in potentially causal genes of the patients’ phenotype, genetical testing was conducted: targeted clinical exome, targeted exome and trio exome, respectively. We identified three novel variants spread throughout the coding sequence of SOX10. The c.395C>G variant identified de novo in patient 1 was a single nucleotide substitution in exon 2. The c.850G>T variant identified as heterozygous in patient 2 was a loss-of-function variant that generated a premature stop codon. The c.966dupT variant identified in patient 3 was a duplication that generated a premature stop codon. It had been identified in his father, arising a possible germinal mosaicism. According to in silico predictors the variant identified in patient 1 was considered as pathogenic, whereas the other two were classified as likely pathogenic.

Conclusions

An exact description of the mutations responsible for WS provides useful information to explain clinical features of WS and contributes to better genetic counselling of WS patients.

SOX10 requirement for melanoma tumor growth is due, in part, to immune-mediated effects

Developmental factors may regulate the expression of immune modulatory proteins in cancer, linking embryonic development and cancer cell immune evasion. This is particularly relevant in melanoma because immune checkpoint inhibitors are commonly used in the clinic. SRY-box transcription factor 10 (SOX10) mediates neural crest development and is required for melanoma cell growth. In this study, we investigate immune-related targets of SOX10 and observe positive regulation of herpesvirus entry mediator (HVEM) and carcinoembryonic-antigen cell-adhesion molecule 1 (CEACAM1). Sox10 knockout reduces tumor growth in vivo, and this effect is exacerbated in immune-competent models. Modulation of CEACAM1 expression but not HVEM elicits modest effects on tumor growth. Importantly, Sox10 knockout effects on tumor growth are dependent, in part, on CD8+ T cells. Extending this analysis to samples from patients with cutaneous melanoma, we observe a negative correlation with SOX10 and immune-related pathways. These data demonstrate a role for SOX10 in regulating immune checkpoint protein expression and anti-tumor immunity in melanoma.

SOX10: 20 years of phenotypic plurality and current understanding of its developmental function

SOX10 belongs to a family of 20 SRY (sex-determining region Y)-related high mobility group box-containing (SOX) proteins, most of which contribute to cell type specification and differentiation of various lineages. The first clue that SOX10 is essential for development, especially in the neural crest, came with the discovery that heterozygous mutations occurring within and around SOX10 cause Waardenburg syndrome type 4. Since then, heterozygous mutations have been reported in Waardenburg syndrome type 2 (Waardenburg syndrome type without Hirschsprung disease), PCWH or PCW (peripheral demyelinating neuropathy, central dysmyelination, Waardenburg syndrome, with or without Hirschsprung disease), intestinal manifestations beyond Hirschsprung (ie, chronic intestinal pseudo-obstruction), Kallmann syndrome and cancer. All of these diseases are consistent with the regulatory role of SOX10 in various neural crest derivatives (melanocytes, the enteric nervous system, Schwann cells and olfactory ensheathing cells) and extraneural crest tissues (inner ear, oligodendrocytes). The recent evolution of medical practice in constitutional genetics has led to the identification of SOX10 variants in atypical contexts, such as isolated hearing loss or neurodevelopmental disorders, making them more difficult to classify in the absence of both a typical phenotype and specific expertise. Here, we report novel mutations and review those that have already been published and their functional consequences, along with current understanding of SOX10 function in the affected cell types identified through in vivo and in vitro models. We also discuss research options to increase our understanding of the origin of the observed phenotypic variability and improve the diagnosis and medical care of affected patients.

Therapeutic combination silencing VEGF and SOX10 increases the antiangiogenic effect in the mouse melanoma model B16-F10 - in vitro and in vivo studies

INTRODUCTION

Gene therapy is an innovative form of treatment of genetic diseases, in which psiRNA molecules silencing specific genes are applied.

AIM

The study evaluated the anti-tumour effect of psiRNA silencing preparations of the vascular endothelial growth factor (VEGF) and Sry-related HMG-Box gene 10 (SOX10) on melanoma (B16-F10) by inhibiting angiogenesis.

MATERIAL AND METHODS

The preparations based on plasmid vectors psiRNA silencing the gene SOX10 and VEGF that form complexes with cationic lipid (psiRNA/carrier) have been developed. psiRNA preparations were tested on the mouse melanoma cell line B16-F10, both in vitro and in vivo. The silencing activity of transfected melanoma cells with the obtained psiRNA preparations was examined using the qPCR and Western blot methods. The anti-tumour activity of psiRNA preparations on melanoma tumour cells was then evaluated in a mouse in vivo model.

RESULTS

In vitro studies have shown that the B16-F10 cells efficiently transfect non-viral preparations - psiRNA: Lyovec (74-89%). Worth mentioning is the fact that silencing SOX10 in B16-F10 melanoma cells increases the expression of the COL18A1 gene (compared to the preparation inhibiting only VEGF), which codes the endostatin to stop angiogenesis. In vivo results show that the level of haemoglobin in tumours of mice treated with psiRNA formulations was over 6 times lower than controls and tumour mass was 60-80% lower.

CONCLUSIONS

The novel study proves that simultaneous inhibition of SOX10 and VEGF enhances the antiangiogenic action and thus contributes to a significant halt of disease development. In addition, these data expand knowledge about SOX10 regulation and functions.

Epigenetic regulation during melanocyte development and homeostasis

Melanocytes originate in the neural crest as precursor cells which then migrate and proliferate to reach their destination where they differentiate into pigment-producing cells. Melanocytes not only determine the colour of hair, skin and eyes but also protect against the harmful effects of UV irradiation. The establishment of the melanocyte lineage is regulated by a defined set of transcription factors and signalling pathways that direct the specific gene expression programmes underpinning melanoblast specification, survival, migration, proliferation and differentiation. In addition, epigenetic modifiers and replacement histones play key roles in regulating gene expression and its timing during the different steps of this process. Here, we discuss the evidence for the role of epigenetic regulators in melanocyte development and function and how they interact with transcription factors and signalling pathways to establish and maintain this important cell lineage.

Epithelial-to-mesenchymal-like transition events in melanoma

Epithelial-to-mesenchymal transition (EMT), a process through which epithelial tumor cells acquire mesenchymal phenotypic properties, contributes to both metastatic dissemination and therapy resistance in cancer. Accumulating evidence indicates that nonepithelial tumors, including melanoma, can also gain mesenchymal-like properties that increase their metastatic propensity and decrease their sensitivity to therapy. In this review, we discuss recent findings, illustrating the striking similarities-but also knowledge gaps-between the biology of mesenchymal-like state(s) in melanoma and mesenchymal state(s) from epithelial cancers. Based on this comparative analysis, we suggest hypothesis-driven experimental approaches to further deepen our understanding of the EMT-like process in melanoma and how such investigations may pave the way towards the identification of clinically relevant biomarkers for prognosis and new therapeutic strategies.

Glioblastoma epigenome profiling identifies SOX10 as a master regulator of molecular tumour subtype

Glioblastoma frequently exhibits therapy-associated subtype transitions to mesenchymal phenotypes with adverse prognosis. Here, we perform multi-omic profiling of 60 glioblastoma primary tumours and use orthogonal analysis of chromatin and RNA-derived gene regulatory networks to identify 38 subtype master regulators, whose cell population-specific activities we further map in published single-cell RNA sequencing data. These analyses identify the oligodendrocyte precursor marker and chromatin modifier SOX10 as a master regulator in RTK I-subtype tumours. In vitro functional studies demonstrate that SOX10 loss causes a subtype switch analogous to the proneural-mesenchymal transition observed in patients at the transcriptomic, epigenetic and phenotypic levels. SOX10 repression in an in vivo syngeneic graft glioblastoma mouse model results in increased tumour invasion, immune cell infiltration and significantly reduced survival, reminiscent of progressive human glioblastoma. These results identify SOX10 as a bona fide master regulator of the RTK I subtype, with both tumour cell-intrinsic and microenvironmental effects.

Novel SOX10 Mutations in Waardenburg Syndrome: Functional Characterization and Genotype-Phenotype Analysis

Waardenburg syndrome (WS) is a prevalent hearing loss syndrome, concomitant with focal skin pigmentation abnormalities, blue iris, and other abnormalities of neural crest-derived cells, including Hirschsprung’s disease. WS is clinically and genetically heterogeneous and it is classified into four major types WS type I, II, III, and IV (WS1, WS2, WS3, and WS4). WS1 and WS3 have the presence of dystopia canthorum, while WS3 also has upper limb anomalies. WS2 and WS4 do not have the dystopia canthorum, but the presence of Hirschsprung’s disease indicates WS4. There is a more severe subtype of WS4 with peripheral nerve and/or central nervous system involvement, namely peripheral demyelinating neuropathy, central dysmyelinating leukodystrophy, WS, and Hirschsprung’s disease or PCW/PCWH. We characterized the genetic defects underlying WS2, WS4, and the WS4-PCW/PCWH) using Sanger and whole-exome sequencing and cytogenomic microarray in seven patients from six unrelated families, including two with WS2 and five with WS4. We also performed multiple functional studies and analyzed genotype–phenotype correlations. The cohort included a relatively high frequency (80%) of individuals with neurological variants of WS4. Six novel SOX10 mutations were identified, including c.89C > A (p.Ser30^∗), c.207_8 delCG (p.Cys71Hisfs^∗62), c.479T > C (p.Leu160Pro), c.1379 delA (p.Tyr460Leufs^∗42), c.425G > C (p.Trp142Ser), and a 20-nucleotide insertion, c.1155_1174dupGCCCCACTATGGCTCAGCCT (p.Phe392Cysfs^∗117). All pathogenic variants were de novo. The results of reporter assays, western blotting, immunofluorescence, and molecular modeling supported the deleterious effects of the identified mutations and their correlations with phenotypic severity. The prediction of genotype–phenotype correlation and functional pathology, and dominant negative effect vs. haploinsufficiency in SOX10-related WS were influenced not only by site (first two vs. last coding exons) and type of mutation (missense vs. truncation/frameshift), but also by the protein expression level, molecular weight, and amino acid content of the altered protein. This in vitro analysis of SOX10 mutations thus provides a deeper understanding of the mechanisms resulting in specific WS subtypes and allows better prediction of the phenotypic manifestations, though it may not be always applicable to in vivo findings without further investigations.

Sorting Sox: Diverse Roles for Sox Transcription Factors During Neural Crest and Craniofacial Development

Sox transcription factors play many diverse roles during development, including regulating stem cell states, directing differentiation, and influencing the local chromatin landscape. Of the twenty vertebrate Sox factors, several play critical roles in the development the neural crest, a key vertebrate innovation, and the subsequent formation of neural crest-derived structures, including the craniofacial complex. Herein, we review the specific roles for individual Sox factors during neural crest cell formation and discuss how some factors may have been essential for the evolution of the neural crest. Additionally, we describe how Sox factors direct neural crest cell differentiation into diverse lineages such as melanocytes, glia, and cartilage and detail their involvement in the development of specific craniofacial structures. Finally, we highlight several SOXopathies associated with craniofacial phenotypes.

Cross-species analysis of enhancer logic using deep learning

Deciphering the genomic regulatory code of enhancers is a key challenge in biology because this code underlies cellular identity. A better understanding of how enhancers work will improve the interpretation of noncoding genome variation and empower the generation of cell type-specific drivers for gene therapy. Here, we explore the combination of deep learning and cross-species chromatin accessibility profiling to build explainable enhancer models. We apply this strategy to decipher the enhancer code in melanoma, a relevant case study owing to the presence of distinct melanoma cell states. We trained and validated a deep learning model, called DeepMEL, using chromatin accessibility data of 26 melanoma samples across six different species. We show the accuracy of DeepMEL predictions on the CAGI5 challenge, where it significantly outperforms existing models on the melanoma enhancer of IRF4 Next, we exploit DeepMEL to analyze enhancer architectures and identify accurate transcription factor binding sites for the core regulatory complexes in the two different melanoma states, with distinct roles for each transcription factor, in terms of nucleosome displacement or enhancer activation. Finally, DeepMEL identifies orthologous enhancers across distantly related species, where sequence alignment fails, and the model highlights specific nucleotide substitutions that underlie enhancer turnover. DeepMEL can be used from the Kipoi database to predict and optimize candidate enhancers and to prioritize enhancer mutations. In addition, our computational strategy can be applied to other cancer or normal cell types.

Boundary sequences flanking the mouse tyrosinase locus ensure faithful pattern of gene expression

Control of gene expression is dictated by cell-type specific regulatory sequences that physically organize the structure of chromatin, including promoters, enhancers and insulators. While promoters and enhancers convey cell-type specific activating signals, insulators prevent the cross-talk of regulatory elements within adjacent loci and safeguard the specificity of action of promoters and enhancers towards their targets in a tissue specific manner. Using the mouse tyrosinase (Tyr) locus as an experimental model, a gene whose mutations are associated with albinism, we described the chromatin structure in cells at two distinct transcriptional states. Guided by chromatin structure, through the use of Chromosome Conformation Capture (3C), we identified sequences at the 5' and 3' boundaries of this mammalian gene that function as enhancers and insulators. By CRISPR/Cas9-mediated chromosomal deletion, we dissected the functions of these two regulatory elements in vivo in the mouse, at the endogenous chromosomal context, and proved their mechanistic role as genomic insulators, shielding the Tyr locus from the expression patterns of adjacent genes.

Melanoma plasticity and phenotypic diversity: therapeutic barriers and opportunities

An incomplete view of the mechanisms that drive metastasis, the primary cause of cancer-related death, has been a major barrier to development of effective therapeutics and prognostic diagnostics. Increasing evidence indicates that the interplay between microenvironment, genetic lesions, and cellular plasticity drives the metastatic cascade and resistance to therapies. Here, using melanoma as a model, we outline the diversity and trajectories of cell states during metastatic dissemination and therapy exposure, and highlight how understanding the magnitude and dynamics of nongenetic reprogramming in space and time at single-cell resolution can be exploited to develop therapeutic strategies that capitalize on nongenetic tumor evolution.

MEK inhibition remodels the active chromatin landscape and induces SOX10 genomic recruitment in BRAF(V600E) mutant melanoma cells

Background

The MAPK/ERK signaling pathway is an essential regulator of numerous cell processes that are crucial for normal development as well as cancer progression. While much is known regarding MAPK/ERK signal conveyance from the cell membrane to the nucleus, the transcriptional and epigenetic mechanisms that govern gene expression downstream of MAPK signaling are not fully elucidated.

Results

This study employed an integrated epigenome analysis approach to interrogate the effects of MAPK/ERK pathway inhibition on the global transcriptome, the active chromatin landscape, and protein–DNA interactions in 501mel melanoma cells. Treatment of these cells with the small-molecule MEK inhibitor AZD6244 induces hyperpigmentation, widespread gene expression changes including alteration of genes linked to pigmentation, and extensive epigenomic reprogramming of transcriptionally distinct regulatory regions associated with the active chromatin mark H3K27ac. Regulatory regions with differentially acetylated H3K27ac regions following AZD6244 treatment are enriched in transcription factor binding motifs of ETV/ETS and ATF family members as well as the lineage-determining factors MITF and SOX10. H3K27ac-dense enhancer clusters known as super-enhancers show similar transcription factor motif enrichment, and furthermore, these super-enhancers are associated with genes encoding MITF, SOX10, and ETV/ETS proteins. Along with genome-wide resetting of the active enhancer landscape, MEK inhibition also results in widespread SOX10 recruitment throughout the genome, including increased SOX10 binding density at H3K27ac-marked enhancers. Importantly, these MEK inhibitor-responsive enhancers marked by H3K27ac and occupied by SOX10 are located near melanocyte lineage-specific and pigmentation genes and overlap numerous human SNPs associated with pigmentation and melanoma phenotypes, highlighting the variants located within these regions for prioritization in future studies.

Conclusions

These results reveal the epigenetic reprogramming underlying the re-activation of melanocyte pigmentation and developmental transcriptional programs in 501mel cells in response to MEK inhibition and suggest extensive involvement of a MEK-SOX10 axis in the regulation of these processes. The dynamic chromatin changes identified here provide a rich genomic resource for further analyses of the molecular mechanisms governing the MAPK pathway in pigmentation- and melanocyte-associated diseases.

Electronic supplementary material

The online version of this article (10.1186/s13072-019-0297-2) contains supplementary material, which is available to authorized users.

Chromatin remodellers Brg1 and Bptf are required for normal gene expression and progression of oncogenic Braf-driven mouse melanoma

Somatic oncogenic mutation of BRAF coupled with inactivation of PTEN constitute a frequent combination of genomic alterations driving the development of human melanoma. Mice genetically engineered to conditionally express oncogenic Braf^V600E and inactivate Pten in melanocytes following tamoxifen treatment rapidly develop melanoma. While early-stage melanomas comprised melanin-pigmented Mitf and Dct-expressing cells, expression of these and other melanocyte identity genes was lost in later stage tumours that showed histological and molecular characteristics of de-differentiated neural crest type cells. Melanocyte identity genes displayed loss of active chromatin marks and RNA polymerase II and gain of heterochromatin marks, indicating epigenetic reprogramming during tumour progression. Nevertheless, late-stage tumour cells grown in culture re-expressed Mitf, and melanocyte markers and Mitf together with Sox10 coregulated a large number of genes essential for their growth. In this melanoma model, somatic inactivation that the catalytic Brg1 (Smarca4) subunit of the SWI/SNF complex and the scaffolding Bptf subunit of the NuRF complex delayed tumour formation and deregulated large and overlapping gene expression programs essential for normal tumour cell growth. Moreover, we show that Brg1 and Bptf coregulated many genes together with Mitf and Sox10. Together these transcription factors and chromatin remodelling complexes orchestrate essential gene expression programs in mouse melanoma cells.

SOX10 regulates multiple genes to direct eumelanin versus pheomelanin production in domestic rock pigeon

The domesticated rock pigeon (Columba livia) has been bred for hundreds of years to display an immense variety of ornamental attributes such as feather color and color patterns. Color is influenced by multiple loci that impact the type and amount of melanin deposited on the feathers. Pigeons homozygous for the "recessive red" mutation, which causes downregulation of Sox10, display brilliant red feathers instead of blue/black feathers. Sox10 encodes a transcription factor important for melanocyte differentiation and function, but the genes that mediate its promotion of black versus red pigment are unknown. Here, we present a transcriptomic comparison of regenerating feathers from wild-type and recessive red pigeons to identify candidate SOX10 targets. Our results identify both known and novel targets, including many genes not previously implicated in pigmentation. These data highlight the value of using novel, emerging model organisms to gain insight into the genetic basis of pigment variation.

Genome-wide association studies for tick resistance in Bos taurus × Bos indicus crossbred cattle: A deeper look into this intricate mechanism

Rhipicephalus (Boophilus) microplus is the main cattle ectoparasite in tropical areas. Gir × Holstein crossbred cows are well adapted to different production systems in Brazil. In this context, we performed genome-wide association study (GWAS) and post-GWAS analyses for R. microplus resistance in an experimental Gir × Holstein F₂ population. Single nucleotide polymorphisms (SNP) identified in GWAS were used to build gene networks and to investigate the breed of origin for its alleles. Tick artificial infestations were performed during the dry and rainy seasons. Illumina BovineSNP50 BeadChip (Illumina Inc., San Diego, CA) and single-step BLUP procedure was used for GWAS. Post-GWAS analyses were performed by gene ontology terms enrichment and gene transcription factors networks, generated from enriched transcription factors, identified from the promoter sequences of selected gene sets. The genetic origin of marker alleles in the F₂ population was assigned using the breed of origin of alleles approach. Heritability estimates for tick counts were 0.40 ± 0.11 in the rainy season and 0.54 ± 0.11 in the dry season. The top ten 0.5-Mbp windows with the highest percentage of genetic variance explained by SNP markers were found in chromosomes 10 and 23 for both the dry and rainy seasons. Gene network analyses allowed the identification of genes involved with biological processes relevant to immune system functions (TREM1, TREM2, and CD83). Gene-transcription factors network allowed the identification of genes involved with immune functions (MYO5A, TREML1, and PRSS16). In resistant animals, the average proportion of animals showing significant SNPs with paternal and maternal alleles originated from Gir breed was 44.8% whereas the proportion of animals with both paternal and maternal alleles originated from Holstein breed was 11.3%. Susceptible animals showing both paternal and maternal alleles originated from Holstein breed represented 44.6% on average, whereas both paternal and maternal alleles originated from Gir breed animals represented 9.3%. This study allowed us to identify candidate genes for tick resistance in Gir × Holstein crossbreds in both rainy and dry seasons. According to the origin of alleles analysis, we found that most animals classified as resistant showed 2 alleles from Gir breed, while the susceptible ones showed alleles from Holstein. Based on these results, the identified genes may be thoroughly investigated in additional experiments aiming to validate their effects on tick resistance phenotype in cattle.

A direct link between MITF, innate immunity, and hair graying

Melanocyte stem cells (McSCs) and mouse models of hair graying serve as useful systems to uncover mechanisms involved in stem cell self-renewal and the maintenance of regenerating tissues. Interested in assessing genetic variants that influence McSC maintenance, we found previously that heterozygosity for the melanogenesis associated transcription factor, Mitf, exacerbates McSC differentiation and hair graying in mice that are predisposed for this phenotype. Based on transcriptome and molecular analyses of Mitf^mi-vga9/+ mice, we report a novel role for MITF in the regulation of systemic innate immune gene expression. We also demonstrate that the viral mimic poly(I:C) is sufficient to expose genetic susceptibility to hair graying. These observations point to a critical suppressor of innate immunity, the consequences of innate immune dysregulation on pigmentation, both of which may have implications in the autoimmune, depigmenting disease, vitiligo.

Hair pigmentation over the course of a lifetime depends on melanocyte stem cells that reside in the hair follicle. As old hairs fall out and new hairs grow in, melanocyte stem cells serve as a reservoir for the melanocytes that produce the pigment that gives hair its visible color. The loss of these stem cells leads to the growth of nonpigmented, or gray, hairs. Evaluating mouse models of hair graying can reveal key aspects of melanocyte stem cell biology. Using this approach, we discovered a novel role for the melanogenesis associated transcription factor, MITF, in repressing the expression of innate immune genes within cells of the melanocyte lineage. The importance of this repression is revealed in animals that have a predisposition for hair graying. In these animals, artificial elevation of the innate immune response, either through a genetic mechanism or via exposure to viral mimic, results in significant melanocyte and melanocyte stem cell loss and leads to the production of an increased number of gray hairs. These observations highlight the negative effects of innate immune activation on melanocyte and melanocyte stem cell physiology and suggest a connection between viral infection and hair graying.

Beyond MITF: Multiple transcription factors directly regulate the cellular phenotype in melanocytes and melanoma

MITF governs multiple steps in the development of melanocytes, including specification from neural crest, growth, survival, and terminal differentiation. In addition, the level of MITF activity determines the phenotype adopted by melanoma cells, whether invasive, proliferative, or differentiated. However, MITF does not act alone. Here, we review literature on the transcription factors that co-regulate MITF-dependent genes. ChIP-seq studies have indicated that the transcription factors SOX10, YY1, and TFAP2A co-occupy subsets of regulatory elements bound by MITF in melanocytes. Analyses at single loci also support roles for LEF1, RB1, IRF4, and PAX3 acting in combination with MITF, while sequence motif analyses suggest that additional transcription factors colocalize with MITF at many melanocyte-specific regulatory elements. However, the precise biochemical functions of each of these MITF collaborators and their contributions to gene expression remain to be elucidated. Analogous to the transcriptional networks in morphogen-patterned tissues during embryogenesis, we anticipate that the level of MITF activity is controlled not only by the concentration of activated MITF, but also by additional transcription factors that either quantitatively or qualitatively influence the expression of MITF-target genes.

Identification and functional analysis of SOX10 phosphorylation sites in melanoma

The transcription factor SOX10 plays an important role in vertebrate neural crest development, including the establishment and maintenance of the melanocyte lineage. SOX10 is also highly expressed in melanoma tumors, and SOX10 expression increases with tumor progression. The suppression of SOX10 in melanoma cells activates TGF-β signaling and can promote resistance to BRAF and MEK inhibitors. Since resistance to BRAF/MEK inhibitors is seen in the majority of melanoma patients, there is an immediate need to assess the underlying biology that mediates resistance and to identify new targets for combinatorial therapeutic approaches. Previously, we demonstrated that SOX10 protein is required for tumor initiation, maintenance and survival. Here, we present data that support phosphorylation as a mechanism employed by melanoma cells to tightly regulate SOX10 expression. Mass spectrometry identified eight phosphorylation sites contained within SOX10, three of which (S24, S45 and T240) were selected for further analysis based on their location within predicted MAPK/CDK binding motifs. SOX10 mutations were generated at these phosphorylation sites to assess their impact on SOX10 protein function in melanoma cells, including transcriptional activation on target promoters, subcellular localization, and stability. These data further our understanding of SOX10 protein regulation and provide critical information for identification of molecular pathways that modulate SOX10 protein levels in melanoma, with the ultimate goal of discovering novel targets for more effective combinatorial therapeutic approaches for melanoma patients.

ERK-mediated phosphorylation regulates SOX10 sumoylation and targets expression in mutant BRAF melanoma

In human mutant BRAF melanoma cells, the stemness transcription factor FOXD3 is rapidly induced by inhibition of ERK1/2 signaling and mediates adaptive resistance to RAF inhibitors. However, the mechanism underlying ERK signaling control of FOXD3 expression remains unknown. Here we show that SOX10 is both necessary and sufficient for RAF inhibitor-induced expression of FOXD3 in mutant BRAF melanoma cells. SOX10 activates the transcription of FOXD3 by binding to a regulatory element in FOXD3 promoter. Phosphorylation of SOX10 by ERK inhibits its transcription activity toward multiple target genes by interfering with the sumoylation of SOX10 at K55, which is essential for its transcription activity. Finally, depletion of SOX10 sensitizes mutant BRAF melanoma cells to RAF inhibitors in vitro and in vivo. Thus, our work discovers a novel phosphorylation-dependent regulatory mechanism of SOX10 transcription activity and completes an ERK1/2/SOX10/FOXD3/ERBB3 axis that mediates adaptive resistance to RAF inhibitors in mutant BRAF melanoma cells.

Hypoxia-induced HIF1α targets in melanocytes reveal a molecular profile associated with poor melanoma prognosis

Hypoxia and HIF1α signaling direct tissue-specific gene responses regulating tumor progression, invasion, and metastasis. By integrating HIF1α knockdown and hypoxia-induced gene expression changes, this study identifies a melanocyte-specific, HIF1α-dependent/hypoxia-responsive gene expression signature. Integration of these gene expression changes with HIF1α ChIP-Seq analysis identifies 81 HIF1α direct target genes in melanocytes. The expression levels for 10 of the HIF1α direct targets - GAPDH, PKM, PPAT, DARS, DTWD1, SEH1L, ZNF292, RLF, AGTRAP, and GPC6 - are significantly correlated with reduced time of disease-free status in melanoma by logistic regression (P-value = 0.0013) and ROC curve analysis (AUC = 0.826, P-value < 0.0001). This HIF1α-regulated profile defines a melanocyte-specific response under hypoxia, and demonstrates the role of HIF1α as an invasive cell state gatekeeper in regulating cellular metabolism, chromatin and transcriptional regulation, vascularization, and invasion.

Sox10 regulates skin melanocyte proliferation by activating the DNA replication licensing factor MCM5

BACKGROUND

The control of cell proliferation is a fundamental aspect of tissue formation in development and regeneration. A cell type that illustrates this point particularly well is the neural crest-derived melanocyte, the pigment cell of vertebrates, as melanocytes can be followed easily during development and their pigment is directly visible in the integument of the adult. In mammals, melanocytes undergo physiological cycles of loss and proliferative regeneration during the hair cycle, and their proliferation is also critical during wound healing, repigmentation of depigmented lesions, and in melanoma formation and progression. Hence, a thorough analysis of the molecular parameters controlling melanocyte proliferation is crucial for our understanding of the physiology of this cell type both in health and disease.

OBJECTIVE

SOX10 is a critical regulator in melanocytes and melanoma cells, but its specific role in their proliferation is far from clear. In this study we analyze the role of SOX10 in regulating mammalian melanocyte proliferation in a mouse model.

METHODS

The role of SOX10 in melanoblast proliferation was analyzed in Sox10/+ mice by co-staining for melanocyte-specific markers and cell proliferation. In vitro, the role of SOX10 was studied by manipulating its levels using RNAi and analyzing the effects on DNA synthesis and cell growth and on gene expression at the RNA and protein levels.

RESULTS

Reduction of Sox10 gene dose led to a reduction in the number of melanoblasts. Knockdown of Sox10 in melanocytes led to inhibition of cell proliferation and a decrease in the expression of the minichromosome maintenance complex component 5 (MCM5). In fact, SOX10 directly activated MCM5 transcription by binding to conserved SOX10 consensus DNA sequences in the MCM5 promoter. Furthermore, the defect in cell proliferation could be rescued partially by overexpression of MCM5 in Sox10 knockdown melanocytes.

CONCLUSION

The results suggest that the SOX10-MCM5 axis plays an important role in controlling melanocyte proliferation. Our findings provide novel insights into the regulatory mechanisms of melanocyte proliferation and may have implications for our understanding of the roles of SOX10 and MCM5 in abnormal melanocyte proliferation disorders such as cutaneous melanoma.

TFAP2 paralogs regulate melanocyte differentiation in parallel with MITF

Mutations in the gene encoding transcription factor TFAP2A result in pigmentation anomalies in model organisms and premature hair graying in humans. However, the pleiotropic functions of TFAP2A and its redundantly-acting paralogs have made the precise contribution of TFAP2-type activity to melanocyte differentiation unclear. Defining this contribution may help to explain why TFAP2A expression is reduced in advanced-stage melanoma compared to benign nevi. To identify genes with TFAP2A-dependent expression in melanocytes, we profile zebrafish tissue and mouse melanocytes deficient in Tfap2a, and find that expression of a small subset of genes underlying pigmentation phenotypes is TFAP2A-dependent, including Dct, Mc1r, Mlph, and Pmel. We then conduct TFAP2A ChIP-seq in mouse and human melanocytes and find that a much larger subset of pigmentation genes is associated with active regulatory elements bound by TFAP2A. These elements are also frequently bound by MITF, which is considered the "master regulator" of melanocyte development. For example, the promoter of TRPM1 is bound by both TFAP2A and MITF, and we show that the activity of a minimal TRPM1 promoter is lost upon deletion of the TFAP2A binding sites. However, the expression of Trpm1 is not TFAP2A-dependent, implying that additional TFAP2 paralogs function redundantly to drive melanocyte differentiation, which is consistent with previous results from zebrafish. Paralogs Tfap2a and Tfap2b are both expressed in mouse melanocytes, and we show that mouse embryos with Wnt1-Cre-mediated deletion of Tfap2a and Tfap2b in the neural crest almost completely lack melanocytes but retain neural crest-derived sensory ganglia. These results suggest that TFAP2 paralogs, like MITF, are also necessary for induction of the melanocyte lineage. Finally, we observe a genetic interaction between tfap2a and mitfa in zebrafish, but find that artificially elevating expression of tfap2a does not increase levels of melanin in mitfa hypomorphic or loss-of-function mutants. Collectively, these results show that TFAP2 paralogs, operating alongside lineage-specific transcription factors such as MITF, directly regulate effectors of terminal differentiation in melanocytes. In addition, they suggest that TFAP2A activity, like MITF activity, has the potential to modulate the phenotype of melanoma cells.

Systems analysis identifies miR-29b regulation of invasiveness in melanoma

Background

In many cancers, microRNAs (miRs) contribute to metastatic progression by modulating phenotypic reprogramming processes such as epithelial-mesenchymal plasticity. This can be driven by miRs targeting multiple mRNA transcripts, inducing regulated changes across large sets of genes. The miR-target databases TargetScan and DIANA-microT predict putative relationships by examining sequence complementarity between miRs and mRNAs. However, it remains a challenge to identify which miR-mRNA interactions are active at endogenous expression levels, and of biological consequence.

Methods

We developed a workflow to integrate TargetScan and DIANA-microT predictions into the analysis of data-driven associations calculated from transcript abundance (RNASeq) data, specifically the mutual information and Pearson’s correlation metrics. We use this workflow to identify putative relationships of miR-mediated mRNA repression with strong support from both lines of evidence. Applying this approach systematically to a large, published collection of unique melanoma cell lines – the Ludwig Melbourne melanoma (LM-MEL) cell line panel – we identified putative miR-mRNA interactions that may contribute to invasiveness. This guided the selection of interactions of interest for further in vitro validation studies.

Results

Several miR-mRNA regulatory relationships supported by TargetScan and DIANA-microT demonstrated differential activity across cell lines of varying matrigel invasiveness. Strong negative statistical associations for these putative regulatory relationships were consistent with target mRNA inhibition by the miR, and suggest that differential activity of such miR-mRNA relationships contribute to differences in melanoma invasiveness. Many of these relationships were reflected across the skin cutaneous melanoma TCGA dataset, indicating that these observations also show graded activity across clinical samples. Several of these miRs are implicated in cancer progression (miR-211, -340, -125b, −221, and -29b). The specific role for miR-29b-3p in melanoma has not been well studied. We experimentally validated the predicted miR-29b-3p regulation of LAMC1 and PPIC and LASP1, and show that dysregulation of miR-29b-3p or these mRNA targets can influence cellular invasiveness in vitro.

Conclusions

This analytic strategy provides a comprehensive, systems-level approach to identify miR-mRNA regulation in high-throughput cancer data, identifies novel putative interactions with functional phenotypic relevance, and can be used to direct experimental resources for subsequent experimental validation.

Computational scripts are available: http://github.com/uomsystemsbiology/LMMEL-miR-miner

Electronic supplementary material

The online version of this article (doi:10.1186/s12943-016-0554-y) contains supplementary material, which is available to authorized users.