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.

Transcriptomic profiling of tissue environments critical for post-embryonic patterning and morphogenesis of zebrafish skin.

Pigment patterns and skin appendages are prominent features of vertebrate skin. In zebrafish, regularly patterned pigment stripes and an array of calcified scales form simultaneously in the skin during post-embryonic development. Understanding the mechanisms that regulate stripe patterning and scale morphogenesis may lead to the discovery of fundamental mechanisms that govern the development of animal form. To learn about cell types and signaling interactions that govern skin patterning and morphogenesis, we generated and analyzed single-cell transcriptomes of skin from wild-type fish as well as fish having genetic or transgenically induced defects in squamation or pigmentation. These data reveal a previously undescribed population of epidermal cells that express transcripts encoding enamel matrix proteins, suggest hormonal control of epithelial-mesenchymal signaling, clarify the signaling network that governs scale papillae development, and identify a critical role for the hypodermis in supporting pigment cell development. Additionally, these comprehensive single-cell transcriptomic data representing skin phenotypes of biomedical relevance should provide a useful resource for accelerating the discovery of mechanisms that govern skin development and homeostasis.

Haplotype-based analysis resolves missing heritability in oculocutaneous albinism type 1B.

Oculocutaneous albinism (OCA) is a rare disorder of pigment production. Affected individuals have variably decreased global pigmentation and visual-developmental changes that lead to low vision. OCA is notable for significant missing heritability, particularly among individuals with residual pigmentation. Tyrosinase (TYR) is the rate-limiting enzyme in melanin pigment biosynthesis and mutations that decrease enzyme function are one of the most common causes of OCA. We present the analysis of high-depth short-read TYR sequencing data for a cohort of 352 OCA probands, ∼50% of whom were previously sequenced without yielding a definitive diagnostic result. Our analysis identified 66 TYR single-nucleotide variants (SNVs) and small insertion/deletions (indels), 3 structural variants, and a rare haplotype comprised of two common frequency variants (p.Ser192Tyr and p.Arg402Gln) in cis-orientation, present in 149/352 OCA probands. We further describe a detailed analysis of the disease-causing haplotype, p.[Ser192Tyr; Arg402Gln] ("cis-YQ"). Haplotype analysis suggests that the cis-YQ allele arose by recombination and that multiple cis-YQ haplotypes are segregating in OCA-affected individuals and control populations. The cis-YQ allele is the most common disease-causing allele in our cohort, representing 19.1% (57/298) of TYR pathogenic alleles in individuals with type 1 (TYR-associated) OCA. Finally, among the 66 TYR variants, we found several additional alleles defined by a cis-oriented combination of minor, potentially hypomorph-producing alleles at common variant sites plus a second, rare pathogenic variant. Together, these results suggest that identification of phased variants for the full TYR locus are required for an exhaustive assessment for potentially disease-causing alleles.

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.

Correlation of age of onset and clinical severity in Niemann-Pick disease type C1 with lysosomal abnormalities and gene expression.

Niemann-Pick disease type C1 (NPC1) is a rare, prematurely fatal lysosomal storage disorder which exhibits highly variable severity and disease progression as well as a wide-ranging age of onset, from perinatal stages to adulthood. This heterogeneity has made it difficult to obtain prompt diagnosis and to predict disease course. In addition, small NPC1 patient sample sizes have been a limiting factor in acquiring genome-wide transcriptome data. In this study, primary fibroblasts from an extensive cohort of 41 NPC1 patients were used to validate our previous findings that the lysosomal quantitative probe LysoTracker can be used as a predictor for age of onset and disease severity. We also examined the correlation between these clinical parameters and RNA expression data from primary fibroblasts and identified a set of genes that were significantly associated with lysosomal defects or age of onset, in particular neurological symptom onset. Hierarchical clustering showed that these genes exhibited distinct expression patterns among patient subgroups. This study is the first to collect transcriptomic data on such a large scale in correlation with clinical and cellular phenotypes, providing a rich genomic resource to address NPC1 clinical heterogeneity and discover potential biomarkers, disease modifiers, or therapeutic targets.

Improved systemic AAV gene therapy with a neurotrophic capsid in Niemann-Pick disease type C1 mice.

Niemann-Pick C1 disease (NPC1) is a rare, fatal neurodegenerative disease caused by mutations in NPC1, which encodes the lysosomal cholesterol transport protein NPC1. Disease pathology involves lysosomal accumulation of cholesterol and lipids, leading to neurological and visceral complications. Targeting the central nervous system (CNS) from systemic circulation complicates treatment of neurological diseases with gene transfer techniques. Selected and engineered capsids, for example, adeno-associated virus (AAV)-PHP.B facilitate peripheral-to-CNS transfer and hence greater CNS transduction than parental predecessors. We report that systemic delivery to Npc1 m1N/m1N mice using an AAV-PHP.B vector ubiquitously expressing NPC1 led to greater disease amelioration than an otherwise identical AAV9 vector. In addition, viral copy number and biodistribution of GFP-expressing reporters showed that AAV-PHP.B achieved more efficient, albeit variable, CNS transduction than AAV9 in Npc1 m1N/m1N mice. This variability was associated with segregation of two alleles of the putative AAV-PHP.B receptor Ly6a in Npc1 m1N/m1N mice. Our data suggest that robust improvements in NPC1 disease phenotypes occur even with modest CNS transduction and that improved neurotrophic capsids have the potential for superior NPC1 AAV gene therapy vectors.

A custom capture sequence approach for oculocutaneous albinism identifies structural variant alleles at the OCA2 locus.

Oculocutaneous albinism (OCA) is a heritable disorder of pigment production that manifests as hypopigmentation and altered eye development. Exon sequencing of known OCA genes is unsuccessful in producing a complete molecular diagnosis for a significant number of affected individuals. We sequenced the DNA of individuals with OCA using short-read custom capture sequencing that targeted coding, intronic, and noncoding regulatory regions of known OCA genes, and genome-wide association study-associated pigmentation loci. We identified an OCA2 complex structural variant (CxSV), defined by a 143 kb inverted segment reintroduced in intron 1, upstream of the native location. The corresponding CxSV junctions were observed in 11/390 probands screened. The 143 kb CxSV presents in one family as a copy number variant duplication for the 143 kb region. In the remaining 10/11 families, the 143 kb CxSV acquired an additional 184 kb deletion across the same region, restoring exons 3-19 of OCA2 to a copy-number neutral state. Allele-associated haplotype analysis found rare SNVs rs374519281 and rs139696407 are linked with the 143 kb CxSV in both OCA2 alleles. For individuals in which customary molecular evaluation does not reveal a biallelic OCA diagnosis, we recommend preliminary screening for these haplotype-associated rare variants, followed by junction-specific validation for the OCA2 143 kb CxSV.

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.

Strategic vision for improving human health at The Forefront of Genomics.

Starting with the launch of the Human Genome Project three decades ago, and continuing after its completion in 2003, genomics has progressively come to have a central and catalytic role in basic and translational research. In addition, studies increasingly demonstrate how genomic information can be effectively used in clinical care. In the future, the anticipated advances in technology development, biological insights, and clinical applications (among others) will lead to more widespread integration of genomics into almost all areas of biomedical research, the adoption of genomics into mainstream medical and public-health practices, and an increasing relevance of genomics for everyday life. On behalf of the research community, the National Human Genome Research Institute recently completed a multi-year process of strategic engagement to identify future research priorities and opportunities in human genomics, with an emphasis on health applications. Here we describe the highest-priority elements envisioned for the cutting-edge of human genomics going forward-that is, at 'The Forefront of Genomics'.

Genetic background modifies phenotypic severity and longevity in a mouse model of Niemann-Pick disease type C1.

Niemann-Pick disease type C1 (NPC1) is a rare, fatal neurodegenerative disorder characterized by lysosomal accumulation of unesterified cholesterol and glycosphingolipids. These subcellular pathologies lead to phenotypes of hepatosplenomegaly, neurological degeneration and premature death. NPC1 is extremely heterogeneous in the timing of clinical presentation and is associated with a wide spectrum of causative NPC1 mutations. To study the genetic architecture of NPC1, we have generated a new NPC1 mouse model, Npc1em1PavNpc1em1Pav/em1Pav mutants showed notably reduced NPC1 protein compared to controls and displayed the pathological and biochemical hallmarks of NPC1. Interestingly, Npc1em1Pav/em1Pav mutants on a C57BL/6J genetic background showed more severe visceral pathology and a significantly shorter lifespan compared to Npc1em1Pav/em1Pav mutants on a BALB/cJ background, suggesting that strain-specific modifiers contribute to disease severity and survival. QTL analysis for lifespan of 202 backcross N2 mutants on a mixed C57BL/6J and BALB/cJ background detected significant linkage to markers on chromosomes 1 and 7. The discovery of these modifier regions demonstrates that mouse models are powerful tools for analyzing the genetics underlying rare human diseases, which can be used to improve understanding of the variability in NPC1 phenotypes and advance options for patient diagnosis and therapy.This article has an associated First Person interview with the first author of the paper.

Maternal immune activation modifies the course of Niemann-pick disease, type C1 in a gender specific manner.

Niemann-Pick disease, type C1 (NPC1) is a rare neurodegenerative lysosomal storage disease with a wide spectrum of clinical manifestation. Multiple genetic factors influence the NPC1 mouse phenotype, but very little attention has been given to prenatal environmental factors that might have long-term effects on the neuroinflammatory component of NPC1 pathology. Studies in other mouse models of cerebellar ataxia have shown that developmental exposures lead to Purkinje neuron degeneration later in life, suggesting that environmental exposures during development can impact cerebellar biology. Thus, we evaluated the potential effect of maternal immune activation (MIA) on disease progression in an Npc1 mouse model. The MIA paradigm used mimics viral infection using the toll like receptor 3 agonist polyinosinic-polycytidilic acid during gestation. Through phenotypic and pathologic tests, we measured motor and behavioral changes as well as cerebellar neuroinflammation and neurodegeneration. We observed a gender and genotype dependent effect of MIA on the cerebellum. While the effects of MIA have been previously shown to primarily affect male progeny, we observed increased sensitivity of female mutant progeny to prenatal exposure to treatment with polyinosinic-polycytidilic acid. Specifically, prenatal MIA resulted in female NPC1 mutant progeny with greater motor deficits and a corresponding decrease in cerebellar Purkinje neurons. Our data suggest that prenatal environmental exposures may be one factor contributing to the phenotypic variability observed in individuals with NPC1.

NPC1 Deficiency in Mice is Associated with Fetal Growth Restriction, Neonatal Lethality and Abnormal Lung Pathology.

The rare lysosomal storage disorder Niemann-Pick disease type C1 (NPC1) arises from mutation of NPC1, which encodes a lysosomal transmembrane protein essential for normal transport and trafficking of cholesterol and sphingolipids. NPC1 is highly heterogeneous in both clinical phenotypes and age of onset. Previous studies have reported sub-Mendelian survival rates for mice homozygous for various Npc1 mutant alleles but have not studied the potential mechanisms underlying this phenotype. We performed the first developmental analysis of a Npc1 mouse model, Npc1em1Pav, and discovered significant fetal growth restriction in homozygous mutants beginning at E16.5. Npc1em1Pav/em1Pav mice also exhibited cyanosis, increased respiratory effort, and over 50% lethality at birth. Analysis of neonatal lung tissues revealed lipid accumulation, notable abnormalities in surfactant, and enlarged alveolar macrophages, suggesting that lung abnormalities may be associated with neonatal lethality in Npc1em1Pav/em1Pav mice. The phenotypic severity of the Npc1em1Pav model facilitated this first analysis of perinatal lethality and lung pathology in an NPC1 model organism, and this model may serve as a useful resource for developing treatments for respiratory complications seen in NPC1 patients.

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.

Meta-analysis of GWA studies provides new insights on the genetic architecture of skin pigmentation in recently admixed populations.

BACKGROUND: Association studies in recently admixed populations are extremely useful to identify the genetic architecture of pigmentation, due to their high genotypic and phenotypic variation. However, to date only four Genome-Wide Association Studies (GWAS) have been carried out in these populations. RESULTS: We present a GWAS of skin pigmentation in an admixed sample from Cuba (N = 762). Additionally, we conducted a meta-analysis including the Cuban sample, and admixed samples from Cape Verde, Puerto Rico and African-Americans from San Francisco. This meta-analysis is one of the largest efforts so far to characterize the genetic basis of skin pigmentation in admixed populations (N = 2,104). We identified five genome-wide significant regions in the meta-analysis, and explored if the markers observed in these regions are associated with the expression of relevant pigmentary genes in human melanocyte cultures. In three of the regions identified in the meta-analysis (SLC24A5, SLC45A2, and GRM5/TYR), the association seems to be driven by non-synonymous variants (rs1426654, rs16891982, and rs1042602, respectively). The rs16891982 polymorphism is strongly associated with the expression of the SLC45A2 gene. In the GRM5/TYR region, in addition to the rs1042602 non-synonymous SNP located on the TYR gene, variants located in the nearby GRM5 gene have an independent effect on pigmentation, possibly through regulation of gene expression of the TYR gene. We also replicated an association recently described near the MFSD12 gene on chromosome 19 (lead variant rs112332856). Additionally, our analyses support the presence of multiple signals in the OCA2/HERC2/APBA2 region on chromosome 15. A clear causal candidate is the HERC2 intronic variant rs12913832, which has a profound influence on OCA2 expression. This variant has pleiotropic effects on eye, hair, and skin pigmentation. However, conditional and haplotype-based analyses indicate the presence of other variants with independent effects on melanin levels in OCA2 and APBA2. Finally, a follow-up of genome-wide signals identified in a recent GWAS for tanning response indicates that there is a substantial overlap in the genetic factors influencing skin pigmentation and tanning response. CONCLUSIONS: Our meta-analysis of skin pigmentation GWAS in recently admixed populations provides new insights about the genetic architecture of this complex trait.

The Genetics of Human Skin and Hair Pigmentation.

Human skin and hair color are visible traits that can vary dramatically within and across ethnic populations. The genetic makeup of these traits-including polymorphisms in the enzymes and signaling proteins involved in melanogenesis, and the vital role of ion transport mechanisms operating during the maturation and distribution of the melanosome-has provided new insights into the regulation of pigmentation. A large number of novel loci involved in the process have been recently discovered through four large-scale genome-wide association studies in Europeans, two large genetic studies of skin color in Africans, one study in Latin Americans, and functional testing in animal models. The responsible polymorphisms within these pigmentation genes appear at different population frequencies, can be used as ancestry-informative markers, and provide insight into the evolutionary selective forces that have acted to create this human diversity.

Identification of Gene Variants Associated with Melanocyte Stem Cell Differentiation in Mice Predisposed for Hair Graying.

Age-related hair graying is caused by malfunction in the regenerative potential of the adult pigmentation system. The retention of hair color over the life of an organism is dependent on the ability of the melanocyte stem cells and their progeny to produce pigment each time a new hair grows. Age-related hair graying is variable in association with genetic background suggesting that quantitative trait loci influencing this trait can be identified. Identification of these quantitative trait loci may lead to the discovery of novel and interesting genes involved in stem cell biology and/or melanogenesis. With this in mind we developed previously a sensitized, mouse modifier screen and discovered that the DBA/1J background is particularly resistant to melanocyte stem cell differentiation in comparison to the C57BL/6J background. Melanocyte stem cell differentiation generally precedes hair graying and is observed in melanocyte stem cells with age. Using quantitative trait loci analysis, we have now identified three quantitative trait loci on mouse chromosomes 7, 13, and X that are associated with DBA/1J-mediated variability in melanocyte stem cell differentiation. Taking advantage of publicly-available mouse sequence and variant data, in silico protein prediction programs, and whole genome gene expression results we describe a short list of potential candidate genes that we anticipate to be involved in melanocyte stem cell biology in mice.

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.