Single-cell profiling of MC1R-inhibited melanocytes.

The human red hair color (RHC) trait is caused by increased pheomelanin (red-yellow) and reduced eumelanin (black-brown) pigment in skin and hair due to diminished melanocortin 1 receptor (MC1R) function. In addition, individuals harboring the RHC trait are predisposed to melanoma development. While MC1R variants have been established as causative of RHC and are a well-defined risk factor for melanoma, it remains unclear mechanistically why decreased MC1R signaling alters pigmentation and increases melanoma susceptibility. Here, we use single-cell RNA sequencing (scRNA-seq) of melanocytes isolated from RHC mouse models to define a MC1R-inhibited Gene Signature (MiGS) comprising a large set of previously unidentified genes which may be implicated in melanogenesis and oncogenic transformation. We show that one of the candidate MiGS genes, TBX3, a well-known anti-senescence transcription factor implicated in melanoma progression, binds both E-box and T-box elements to regulate genes associated with melanogenesis and senescence bypass. Our results provide key insights into further mechanisms by which melanocytes with reduced MC1R signaling may regulate pigmentation and offer new candidates of study toward understanding how individuals with the RHC phenotype are predisposed to melanoma.

The MFSD12 p.Tyr182His common variant is sufficient to alter mouse agouti coat color.

MFSD12 functions as a transmembrane protein required for import of cysteine into melanosomes and lysosomes. The MFSD12 locus has been associated with phenotypic variation in skin color across African, Latin American, and East Asian populations. The frequency of a particular MFSD12 coding variant, rs2240751 (MAF = 0.08), has been reported to correlate with solar radiation and occur at highest frequency in Peruvian (PEL MAF = 0.48) and Han Chinese (CHB MAF = 0.40) populations, suggesting it could be causative for associated phenotypic variation in skin color. We have generated a mouse knock-in allele, Mfsd12Y182H , to model the human missense p.Tyr182His human variant. We demonstrate that the variant transcript is stably expressed and that agouti mice homozygote for the variant allele are viable with an altered coat color. This in vivo data confirms that the MFSD12 p.Tyr182His variant functions as a hypomorphic allele sufficient to alter mammalian pigmentation.

Loss of MC1R signaling implicates TBX3 in pheomelanogenesis and melanoma predisposition.

The human Red Hair Color (RHC) trait is caused by increased pheomelanin (red-yellow) and reduced eumelanin (black-brown) pigment in skin and hair due to diminished melanocortin 1 receptor (MC1R) function. In addition, individuals harboring the RHC trait are predisposed to melanoma development. While MC1R variants have been established as causative of RHC and are a well-defined risk factor for melanoma, it remains unclear mechanistically why decreased MC1R signaling alters pigmentation and increases melanoma susceptibility. Here, we use single-cell RNA-sequencing (scRNA-seq) of melanocytes isolated from RHC mouse models to reveal a Pheomelanin Gene Signature (PGS) comprising genes implicated in melanogenesis and oncogenic transformation. We show that TBX3, a well-known anti-senescence transcription factor implicated in melanoma progression, is part of the PGS and binds both E-box and T-box elements to regulate genes associated with melanogenesis and senescence bypass. Our results provide key insights into mechanisms by which MC1R signaling regulates pigmentation and how individuals with the RHC phenotype are predisposed to melanoma.

Transcriptome of HPßCD-treated Niemann-Pick disease type C1 cells highlights GPNMB as a biomarker for therapeutics.

The rare, fatal neurodegenerative disorder Niemann-Pick disease type C1 (NPC1) arises from lysosomal accumulation of unesterified cholesterol and glycosphingolipids. These subcellular pathologies lead to phenotypes of hepatosplenomegaly, neurological degeneration and premature death. The timing and severity of NPC1 clinical presentation is extremely heterogeneous. This study analyzed RNA-Seq data from 42 NPC1 patient-derived, primary fibroblast cell lines to determine transcriptional changes induced by treatment with 2-hydroxypropyl-ß-cyclodextrin (HPßCD), a compound currently under investigation in clinical trials. A total of 485 HPßCD-responsive genes were identified. Pathway enrichment analysis of these genes showed significant involvement in cholesterol and lipid biosynthesis. Furthermore, immunohistochemistry of the cerebellum as well as measurements of plasma from Npc1m1N null mice treated with HPßCD and adeno-associated virus gene therapy suggests that one of the identified genes, GPNMB, may serve as a useful biomarker of treatment response in NPC1 disease. Overall, this large NPC1 patient-derived dataset provides a comprehensive foundation for understanding the genomic response to HPßCD treatment.

A curated gene list for expanding the horizons of pigmentation biology.

Over the past century, studies of human pigmentary disorders along with mouse and zebrafish models have shed light on the many cellular functions associated with visible pigment phenotypes. This has led to numerous genes annotated with the ontology term "pigmentation" in independent human, mouse, and zebrafish databases. Comparisons among these datasets revealed that each is individually incomplete in documenting all genes involved in integument-based pigmentation phenotypes. Additionally, each database contained inherent species-specific biases in data annotation, and the term "pigmentation" did not solely reflect integument pigmentation phenotypes. This review presents a comprehensive, cross-species list of 650 genes involved in pigmentation phenotypes that was compiled with extensive manual curation of genes annotated in OMIM, MGI, ZFIN, and GO. The resulting cross-species list of genes both intrinsic and extrinsic to integument pigment cells provides a valuable tool that can be used to expand our knowledge of complex, pigmentation-associated pathways.

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

SOX10 is required for melanocyte development and maintenance, and has been linked to melanoma initiation and progression. However, the molecular mechanisms by which SOX10 guides the appropriate gene expression programs necessary to promote the melanocyte lineage are not fully understood. Here we employ genetic and epigenomic analysis approaches to uncover novel genomic targets and previously unappreciated molecular roles of SOX10 in melanocytes. Through global analysis of SOX10-binding sites and epigenetic characteristics of chromatin states, we uncover an extensive catalog of SOX10 targets genome-wide. Our findings reveal that SOX10 predominantly engages 'open' chromatin regions and binds to distal regulatory elements, including novel and previously known melanocyte enhancers. Integrated chromatin occupancy and transcriptome analysis suggest a role for SOX10 in both transcriptional activation and repression to regulate functionally distinct classes of genes. We demonstrate that distinct epigenetic signatures and cis-regulatory sequence motifs predicted to bind putative co-regulatory transcription factors define SOX10-activated and SOX10-repressed target genes. Collectively, these findings uncover a central role of SOX10 as a global regulator of gene expression in the melanocyte lineage by targeting diverse regulatory pathways.

The transcription factors Ets1 and Sox10 interact during murine melanocyte development.

Melanocytes, the pigment-producing cells, arise from multipotent neural crest (NC) cells during embryogenesis. Many genes required for melanocyte development were identified using mouse pigmentation mutants. The variable spotting mouse pigmentation mutant arose spontaneously at the Jackson Laboratory. We identified a G-to-A nucleotide transition in exon 3 of the Ets1 gene in variable spotting, which results in a missense G102E mutation. Homozygous variable spotting mice exhibit sporadic white spotting. Similarly, mice carrying a targeted deletion of Ets1 exhibit hypopigmentation; nevertheless, the function of Ets1 in melanocyte development is unknown. The transcription factor Ets1 is widely expressed in developing organs and tissues, including the NC. In the chick, Ets1 is required for the expression of Sox10, a transcription factor critical for the development of various NC derivatives, including melanocytes. We show that Ets1 is required early for murine NC cell and melanocyte precursor survival in vivo. Given the importance of Ets1 for Sox10 expression in the chick, we investigated a potential genetic interaction between these genes by comparing the hypopigmentation phenotypes of single and double heterozygous mice. The incidence of hypopigmentation in double heterozygotes was significantly greater than in single heterozygotes. The area of hypopigmentation in double heterozygotes was significantly larger than would be expected from the addition of the areas of hypopigmentation of single heterozygotes, suggesting that Ets1 and Sox10 interact synergistically in melanocyte development. Since Sox10 is also essential for enteric ganglia development, we examined the distal colons of Ets1 null mutants and found a significant decrease in enteric innervation, which was exacerbated by Sox10 heterozygosity. At the molecular level, Ets1 was found to activate an enhancer critical for Sox10 expression in NC-derived structures. Furthermore, enhancer activation was significantly inhibited by the variable spotting mutation. Together, these results suggest that Ets1 and Sox10 interact to promote proper melanocyte and enteric ganglia development from the NC.

Ectopic differentiation of melanocyte stem cells is influenced by genetic background.

Hair graying in mouse is attributed to the loss of melanocyte stem cell function and the progressive depletion of the follicular melanocyte population. Single-gene, hair graying mouse models have pointed to a number of critical pathways involved in melanocyte stem cell biology; however, the broad range of phenotypic variation observed in human hair graying suggests that additional genetic variants involved in this process may yet be discovered. Using a sensitized approach, we ask here whether natural genetic variation influences a predominant cellular mechanism of hair graying in mouse, melanocyte stem cell differentiation. We developed an innovative method to quantify melanocyte stem cell differentiation by measuring ectopically pigmented melanocyte stem cells in response to the melanocyte-specific transgene Tg(Dct-Sox10). We make the novel observation that the production of ectopically pigmented melanocyte stem cells varies considerably across strains. The success of sensitizing for melanocyte stem cell differentiation by Tg(Dct-Sox10) sets the stage for future investigations into the genetic basis of strain-specific contributions to melanocyte stem cell biology.

A unique missense allele of BAF155, a core BAF chromatin remodeling complex protein, causes neural tube closure defects in mice.

Failure of embryonic neural tube closure results in the second most common class of birth defects known as neural tube defects (NTDs). While NTDs are likely the result of complex multigenic dysfunction, it is not known whether polymorphisms in epigenetic regulators may be risk factors for NTDs. Here we characterized Baf155(msp3) , a unique ENU-induced allele in mice. Homozygous Baf155(mps3) embryos exhibit highly penetrant exencephaly, allowing us to investigate the roles of an assembled, but malfunctional BAF chromatin remodeling complex in vivo at the time of neural tube closure. Evidence of defects in proliferation and apoptosis were found within the neural tube. RNA-Seq analysis revealed that surprisingly few genes showed altered expression in Baf155 mutant neural tissue, given the broad epigenetic role of the BAF complex, but included genes involved in neural development and cell survival. Moreover, gene expression changes between individual mutants were variable even though the NTD was consistently observed. This suggests that inconsistent gene regulation contributes to failed neural tube closure. These results shed light on the role of the BAF complex in the process of neural tube closure and highlight the importance of studying missense alleles to understand epigenetic regulation during critical phases of development.

SOX10 ablation arrests cell cycle, induces senescence, and suppresses melanomagenesis.

The transcription factor SOX10 is essential for survival and proper differentiation of neural crest cell lineages, where it plays an important role in the generation and maintenance of melanocytes. SOX10 is also highly expressed in melanoma tumors, but a role in disease progression has not been established. Here, we report that melanoma tumor cell lines require wild-type SOX10 expression for proliferation and SOX10 haploinsufficiency reduces melanoma initiation in the metabotropic glutamate receptor 1 (Grm1(Tg)) transgenic mouse model. Stable SOX10 knockdown in human melanoma cells arrested cell growth, altered cellular morphology, and induced senescence. Melanoma cells with stable loss of SOX10 were arrested in the G1 phase of the cell cycle, with reduced expression of the melanocyte determining factor microphthalmia-associated transcription factor, elevated expression of p21WAF1 and p27KIP2, hypophosphorylated RB, and reduced levels of its binding partner E2F1. As cell-cycle dysregulation is a core event in neoplastic transformation, the role for SOX10 in maintaining cell-cycle control in melanocytes suggests a rational new direction for targeted treatment or prevention of melanoma.

Efficacy of N-acetylcysteine in phenotypic suppression of mouse models of Niemann-Pick disease, type C1.

Niemann-Pick disease, type C1 (NPC1), which arises from a mutation in the NPC1 gene, is characterized by abnormal cellular storage and transport of cholesterol and other lipids that leads to hepatic disease and progressive neurological impairment. Oxidative stress has been hypothesized to contribute to the NPC1 disease pathological cascade. To determine whether treatments reducing oxidative stress could alleviate NPC1 disease phenotypes, the in vivo effects of the antioxidant N-acetylcysteine (NAC) on two mouse models for NPC1 disease were studied. NAC was able to partially suppress phenotypes in both antisense-induced (NPC1ASO) and germline (Npc1-/-) knockout genetic mouse models, confirming the presence of an oxidative stress-related mechanism in progression of NPC1 phenotypes and suggesting NAC as a potential molecule for treatment. Gene expression analyses of NAC-treated NPC1ASO mice suggested NAC affects pathways distinct from those initially altered by Npc1 knockdown, data consistent with NAC achieving partial disease phenotype suppression. In a therapeutic trial of short-term NAC administration to NPC1 patients, no significant effects on oxidative stress in these patients were identified other than moderate improvement of the fraction of reduced CoQ10, suggesting limited efficacy of NAC monotherapy. However, the mouse model data suggest that the distinct antioxidant effects of NAC could provide potential treatment of NPC1 disease, possibly in concert with other therapeutic molecules at earlier stages of disease progression. These data also validated the NPC1ASO mouse as an efficient model for candidate NPC1 drug screening, and demonstrated similarities in hepatic phenotypes and genome-wide transcript expression patterns between the NPC1ASO and Npc1-/- models.

Comparison of melanoblast expression patterns identifies distinct classes of genes.

A full understanding of transcriptional regulation requires integration of information obtained from multiple experimental datasets. These include datasets annotating gene expression within the context of an entire organism under normal and genetically perturbed conditions. Here we describe an expression dataset annotating pigment cell-expressed genes of the developing melanocyte and retinal pigmented epithelium lineages. Expression images are annotated and available at http://research.nhgri.nih.gov/manuscripts/Loftus/March2009/. Data are also summarized in a standardized manner using a universal melanoblast scoring scale that accounts for the embryonic location of cells and regional cell density. This approach allowed us to classify 14 pigment genes into four groupings classified by cell lineage expression, temporal-spatial context, and differential alteration in response to altered MITF and SOX10 status. Significant differences in regional populations were also observed across inbred strain backgrounds, highlighting the value of this approach to identify modifier allele influences on melanoblast number and distributions. This analysis revealed novel features of in vivo expression patterns that are not measurable by in vitro-based assays, providing data that in combination with genomic analyses will allow modeling of pigment cell gene expression in development and disease.

A Sox10 expression screen identifies an amino acid essential for Erbb3 function.

The neural crest (NC) is a population of embryonic stem cells that gives rise to numerous cell types, including the glia and neurons of the peripheral and enteric nervous systems and the melanocytes of the skin and hair. Mutations in genes and genetic pathways regulating NC development lead to a wide spectrum of human developmental disorders collectively called neurocristopathies. To identify molecular pathways regulating NC development and to understand how alterations in these processes lead to disease, we established an N-ethyl-N-nitrosourea (ENU) mutagenesis screen utilizing a mouse model sensitized for NC defects, Sox10(LacZ/+). Out of 71 pedigrees analyzed, we identified and mapped four heritable loci, called modifier of Sox10 expression pattern 1-4 (msp1-4), which show altered NC patterning. In homozygous msp1 embryos, Sox10(LacZ) expression is absent in cranial ganglia, cranial nerves, and the sympathetic chain; however, the development of other Sox10-expressing cells appears unaffected by the mutation. Linkage analysis, sequencing, and complementation testing confirmed that msp1 is a new allele of the receptor tyrosine kinase Erbb3, Erbb3(msp1), that carries a single amino acid substitution in the extracellular region of the protein. The ENU-induced mutation does not alter protein expression, however, it is sufficient to impair ERBB3 signaling such that the embryonic defects observed in msp1 resemble those of Erbb3 null alleles. Biochemical analysis of the mutant protein showed that ERBB3 is expressed on the cell surface, but its ligand-induced phosphorylation is dramatically reduced by the msp1 mutation. These findings highlight the importance of the mutated residue for ERBB3 receptor function and activation. This study underscores the utility of using an ENU mutagenesis to identify genetic pathways regulating NC development and to dissect the roles of discrete protein domains, both of which contribute to a better understanding of gene function in a cellular and developmental setting.

A sensitized mutagenesis screen identifies Gli3 as a modifier of Sox10 neurocristopathy.

Haploinsufficiency for the transcription factor SOX10 is associated with the pigmentary deficiencies of Waardenburg syndrome (WS) and is modeled in Sox10 haploinsufficient mice (Sox10(LacZ/+)). As genetic background affects WS severity in both humans and mice, we established an N-ethyl-N-nitrosourea (ENU) mutagenesis screen to identify modifiers that increase the phenotypic severity of Sox10(LacZ/+) mice. Analysis of 230 pedigrees identified three modifiers, named modifier of Sox10 neurocristopathies (Mos1, Mos2 and Mos3). Linkage analysis confirmed their locations on mouse chromosomes 13, 4 and 3, respectively, within regions distinct from previously identified WS loci. Positional candidate analysis of Mos1 identified a truncation mutation in a hedgehog(HH)-signaling mediator, GLI-Kruppel family member 3 (Gli3). Complementation tests using a second allele of Gli3 (Gli3(Xt-J)) confirmed that a null mutation of Gli3 causes the increased hypopigmentation in Sox10(LacZ/+);Gli3(Mos1/)(+) double heterozygotes. Early melanoblast markers (Mitf, Sox10, Dct, and Si) are reduced in Gli3(Mos1/)(Mos1) embryos, indicating that loss of GLI3 signaling disrupts melanoblast specification. In contrast, mice expressing only the GLI3 repressor have normal melanoblast specification, indicating that the full-length GLI3 activator is not required for specification of neural crest to the melanocyte lineage. This study demonstrates the feasibility of sensitized screens to identify disease modifier loci and implicates GLI3 and other HH signaling components as modifiers of human neurocristopathies.