Search

Most expression quantitative trait locus (eQTL) studies to date have been performed in heterogeneous tissues as opposed to specific cell types. To better understand the cell-type-specific regulatory landscape of human melanocytes, which give rise to melanoma but account for <5% of typical human skin biopsies, we performed an eQTL analysis in primary melanocyte cultures from 106 newborn males. We identified 597,335 cis -eQTL SNPs prior to linkage disequilibrium (LD) pruning and 4997 eGenes (FDR < 0.05). Melanocyte eQTLs differed considerably from those identified in the 44 GTEx tissue types, including skin. Over a third of melanocyte eGenes, including key genes in melanin synthesis pathways, were unique to melanocytes compared to those of GTEx skin tissues or TCGA melanomas. The melanocyte data set also identified trans -eQTLs, including those connecting a pigmentation-associated functional SNP with four genes, likely through cis -regulation of IRF4 Melanocyte eQTLs are enriched in cis -regulatory signatures found in melanocytes as well as in melanoma-associated variants identified through genome-wide association studies. Melanocyte eQTLs also colocalized with melanoma GWAS variants in five known loci. Finally, a transcriptome-wide association study using melanocyte eQTLs uncovered four novel susceptibility loci, where imputed expression levels of five genes ( ZFP90 , HEBP1 , MSC , CBWD1 , and RP11-383H13.1 ) were associated with melanoma at genome-wide significant P -values. Our data highlight the utility of lineage-specific eQTL resources for annotating GWAS findings, and present a robust database for genomic research of melanoma risk and melanocyte biology., (© 2018 Zhang et al.; Published by Cold Spring Harbor Laboratory Press.)

We take a comprehensive approach to the study of regulatory control of gene expression in melanocytes that proceeds from large-scale enhancer discovery facilitated by ChIP-seq; to rigorous validation in silico, in vitro, and in vivo; and finally to the use of machine learning to elucidate a regulatory vocabulary with genome-wide predictive power. We identify 2489 putative melanocyte enhancer loci in the mouse genome by ChIP-seq for EP300 and H3K4me1. We demonstrate that these putative enhancers are evolutionarily constrained, enriched for sequence motifs predicted to bind key melanocyte transcription factors, located near genes relevant to melanocyte biology, and capable of driving reporter gene expression in melanocytes in culture (86%; 43/50) and in transgenic zebrafish (70%; 7/10). Next, using the sequences of these putative enhancers as a training set for a supervised machine learning algorithm, we develop a vocabulary of 6-mers predictive of melanocyte enhancer function. Lastly, we demonstrate that this vocabulary has genome-wide predictive power in both the mouse and human genomes. This study provides deep insight into the regulation of gene expression in melanocytes and demonstrates a powerful approach to the investigation of regulatory sequences that can be applied to other cell types.

Treacher Collins syndrome (TCOF1) is a dominant disorder of craniofacial development that has been linked previosuly to a region of chromosome 5q31.3-32. Identification of recombination events in affected individuals has reduced the candidate gene region to a 0.5-Mb area between the loci RPS14 (proximal) and ANX6 [distal]. A transcriptional map of this candidate gene region, generated by analysis of exon amplification clones, has identified the genomic location of four genes, heparan sulfate-N-sulfotransferase-N-deacetylase, glutathione peroxidase, as well as two novel, previously uncharacterized genes. Each of these genes, based on their location, must be considered candidates for TCOF1 locus.