Despite the wide range of skin pigmentation in humans, little is known about its genetic basis in global populations. Examining ethnically diverse African genomes, we identify variants in or near SLC24A5, MFSD12, DDB1, TMEM138, OCA2 and HERC2 that are significantly associated with skin pigmentation. Genetic evidence indicates that the light pigmentation variant at SLC24A5 was introduced into East Africa by gene flow from non-Africans. At all other loci, variants associated with dark pigmentation in Africans are identical by descent in southern Asian and Australo-Melanesian populations. Functional analyses indicate that MFSD12 encodes a lysosomal protein that affects melanogenesis in zebrafish and mice, and that mutations in melanocyte-specific regulatory regions near DDB1/TMEM138 correlate with expression of UV response genes under selection in Eurasians.

Melanomas originating from mucosal surfaces have low mutation burden, genomic instability, and poor prognosis. To identify potential driver genes, we sequenced hundreds of cancer-related genes in 43 human mucosal melanomas, cataloging point mutations, amplifications, and deletions. The SPRED1 gene, which encodes a negative regulator of mitogen-activated protein kinase (MAPK) signaling, was inactivated in 37% of the tumors. Four distinct genotypes were associated with SPRED1 loss. Using a rapid, tissue-specific CRISPR technique to model these genotypes in zebrafish, we found that SPRED1 functions as a tumor suppressor, particularly in the context of KIT mutations. SPRED1 knockdown caused MAPK activation, increased cell proliferation, and conferred resistance to drugs inhibiting KIT tyrosine kinase activity. These findings provide a rationale for MAPK inhibition in SPRED1-deficient melanomas and introduce a zebrafish modeling approach that can be used more generally to dissect genetic interactions in cancer.