Host cell invasion by Toxoplasma gondii and other apicomplexan parasites requires transmembrane adhesins that mediate binding to receptors on the substrate and host cell to facilitate motility and invasion. Rhomboid proteases (ROMs) are thought to cleave adhesins within their transmembrane segments, thus allowing the parasite to disengage from receptors and completely enter the host cell. To examine the specific roles of individual ROMs during invasion, we generated single, double, and triple knockouts for the three ROMs expressed in T. gondii tachyzoites. Analysis of these mutants demonstrated that ROM4 is the primary protease involved in adhesin processing and host cell invasion, whereas ROM1 or ROM5 plays negligible roles in these processes. Deletion of ROM4 blocked the shedding of adhesins such as MIC2 (microneme protein 2), causing them to accumulate on the surface of extracellular parasites. Increased surface adhesins led to nonproductive attachment, altered gliding motility, impaired moving junction formation, and reduced invasion efficiency. Despite the importance of ROM4 for efficient invasion, mutants lacking all three ROMs were viable and MIC2 was still efficiently removed from the surface of invaded mutant parasites, implying the existence of ROM-independent mechanisms for adhesin removal during invasion. Collectively, these results suggest that although ROM processing of adhesins is not absolutely essential, it is important for efficient host cell invasion by T. gondii. Importance: Apicomplexan parasites such as Toxoplasma gondii express surface proteins that bind host cell receptors to aid invasion. Many of these adhesins are subject to cleavage by rhomboid proteases (ROMs) within their transmembrane segments during invasion. Previous studies have demonstrated the importance of adhesin cleavage for parasite invasion and proposed that the ROMs responsible for processing would be essential for parasite survival. In T. gondii, ROM5 was thought to be the critical ROM for adhesin shedding due to its robust protease activity in vitro and posterior localization on the parasite surface. Here, we knocked out all three ROMs in T. gondii tachyzoites and found that ROM4, but not ROM5, was key for adhesin cleavage. However, none of the ROMs individually or in combination was essential for cell entry, further emphasizing that essential pathways such as invasion typically rely on redundant pathways to ensure survival.

Apicomplexan parasites actively invade host cells using a mechanism predicted to be powered by a parasite actin-dependent myosin motor. In the model apicomplexan Toxoplasma gondii, inducible knockout of the actin gene, ACT1, was recently demonstrated to limit but not completely abolish invasion. This observation has led to the provocative suggestion that T. gondii possesses alternative, ACT1-independent invasion pathways. Here, we dissected the residual invasive ability of Δact1 parasites. Surprisingly, we were able to detect residual ACT1 protein in inducible Δact1 parasites as long as 5 days after ACT1 deletion. We further found that the longer Δact1 parasites were propagated after ACT1 deletion, the more severe an invasion defect was observed. Both findings are consistent with the quantity of residual ACT1 retained in Δact1 parasites being responsible for their invasive ability. Furthermore, invasion by the Δact1 parasites was also sensitive to the actin polymerization inhibitor cytochalasin D. Finally, there was no clear defect in attachment to host cells or moving junction formation by Δact1 parasites. However, Δact1 parasites often exhibited delayed entry into host cells, suggesting a defect specific to the penetration stage of invasion. Overall, our results support a model where residual ACT1 protein retained in inducible Δact1 parasites facilitates their limited invasive ability and confirm that parasite actin is essential for efficient penetration into host cells during invasion.

Phagocytic cells are the first line of innate defense against intracellular pathogens, and yet Toxoplasma gondii is renowned for its ability to survive in macrophages, although this paradigm is based on virulent type I parasites. Surprisingly, we find that avirulent type III parasites are preferentially cleared in naive macrophages, independent of gamma interferon (IFN-γ) activation. The ability of naive macrophages to clear type III parasites was dependent on enhanced activity of NADPH oxidase (Nox)-generated reactive oxygen species (ROS) and induction of guanylate binding protein 5 (Gbp5). Macrophages infected with type III parasites (CTG strain) showed a time-dependent increase in intracellular ROS generation that was higher than that induced by type I parasites (GT1 strain). The absence of Nox1 or Nox2, gp91 subunit isoforms of the Nox complex, reversed ROS-mediated clearance of CTG parasites. Consistent with this finding, both Nox1-/- and Nox2-/- mice showed higher susceptibility to CTG infection than wild-type mice. Additionally, Gbp5 expression was induced upon infection and the enhanced clearance of CTG strain parasites was reversed in Gbp5-/- macrophages. Expression of a type I ROP18 allele in CTG prevented clearance in naive macrophages, suggesting that it plays a role counteracting Gbp5. Although ROS and Gbp5 have been linked to activation of the NLRP3 inflammasome, clearance of CTG parasites did not rely on induction of pyroptosis. Collectively, these findings reveal that not all strains of T. gondii are adept at avoiding clearance in macrophages and define new roles for ROS and Gbps in controlling this important intracellular pathogen.IMPORTANCEToxoplasma infections in humans and other mammals are largely controlled by IFN-γ produced by the activated adaptive immune system. However, we still do not completely understand the role of cell-intrinsic functions in controlling Toxoplasma or other apicomplexan infections. The present work identifies intrinsic activities in naive macrophages in counteracting T. gondii infection. Using an avirulent strain of T. gondii, we highlight the importance of Nox complexes in conferring protection against parasite infection both in vitro and in vivo We also identify Gbp5 as a novel macrophage factor involved in limiting intracellular infection by avirulent strains of T. gondii The rarity of human infections caused by type III strains suggests that these mechanisms may also be important in controlling human toxoplasmosis. These findings further extend our understanding of host responses and defense mechanisms that act to control parasitic infections at the cellular level.

Deleterious germline variants in CDKN2A account for around 40% of familial melanoma cases, and rare variants in CDK4, BRCA2, BAP1 and the promoter of TERT have also been linked to the disease. Here we set out to identify new high-penetrance susceptibility genes by sequencing 184 melanoma cases from 105 pedigrees recruited in the UK, The Netherlands and Australia that were negative for variants in known predisposition genes. We identified families where melanoma cosegregates with loss-of-function variants in the protection of telomeres 1 gene (POT1), with a proportion of family members presenting with an early age of onset and multiple primary tumors. We show that these variants either affect POT1 mRNA splicing or alter key residues in the highly conserved oligonucleotide/oligosaccharide-binding (OB) domains of POT1, disrupting protein-telomere binding and leading to increased telomere length. These findings suggest that POT1 variants predispose to melanoma formation via a direct effect on telomeres.

At least 17 genomic regions are established as harboring melanoma susceptibility variants, in most instances with genome-wide levels of significance and replication in independent samples. Based on genome-wide single nucleotide polymorphism (SNP) data augmented by imputation to the 1,000 Genomes reference panel, we have fine mapped these regions in over 5,000 individuals with melanoma (mainly from the GenoMEL consortium) and over 7,000 ethnically matched controls. A penalized regression approach was used to discover those SNP markers that most parsimoniously explain the observed association in each genomic region. For the majority of the regions, the signal is best explained by a single SNP, which sometimes, as in the tyrosinase region, is a known functional variant. However in five regions the explanation is more complex. At the CDKN2A locus, for example, there is strong evidence that not only multiple SNPs but also multiple genes are involved. Our results illustrate the variability in the biology underlying genome-wide susceptibility loci and make steps toward accounting for some of the "missing heritability."

Apicomplexan parasites contain a conserved protein CelTOS that, in malaria parasites, is essential for traversal of cells within the mammalian host and arthropod vector. However, the molecular role of CelTOS is unknown because it lacks sequence similarity to proteins of known function. Here, we determined the crystal structure of CelTOS and discovered CelTOS resembles proteins that bind to and disrupt membranes. In contrast to known membrane disruptors, CelTOS has a distinct architecture, specifically binds phosphatidic acid commonly present within the inner leaflet of plasma membranes, and potently disrupts liposomes composed of phosphatidic acid by forming pores. Microinjection of CelTOS into cells resulted in observable membrane damage. Therefore, CelTOS is unique as it achieves nearly universal inner leaflet cellular activity to enable the exit of parasites from cells during traversal. By providing novel molecular insight into cell traversal by apicomplexan parasites, our work facilitates the design of therapeutics against global pathogens.

IFN-γ activates cells to restrict intracellular pathogens by upregulating cellular effectors including the p65 family of guanylate-binding proteins (GBPs). Here we test the role of Gbp1 in the IFN-γ-dependent control of T. gondii in the mouse model. Virulent strains of T. gondii avoided recruitment of Gbp1 to the parasitophorous vacuole in a strain-dependent manner that was mediated by the parasite virulence factors ROP18, an active serine/threonine kinase, and the pseudokinase ROP5. Increased recruitment of Gbp1 to Δrop18 or Δrop5 parasites was associated with clearance in IFN-γ-activated macrophages in vitro, a process dependent on the autophagy protein Atg5. The increased susceptibility of Δrop18 mutants in IFN-γ-activated macrophages was reverted in Gbp1(-/-) cells, and decreased virulence of this mutant was compensated in Gbp1(-/-) mice, which were also more susceptible to challenge with type II strain parasites of intermediate virulence. These findings demonstrate that Gbp1 plays an important role in the IFN-γ-dependent, cell-autonomous control of toxoplasmosis and predict a broader role for this protein in host defense.

Virulence of complex pathogens in mammals is generally determined by multiple components of the pathogen interacting with the functional complexity and multiple layering of the mammalian immune system. It is most unusual for the resistance of a mammalian host to be overcome by the defeat of a single defence mechanism. In this study we uncover and analyse just such a case at the molecular level, involving the widespread intracellular protozoan pathogen Toxoplasma gondii and one of its most important natural hosts, the house mouse (Mus musculus). Natural polymorphism in virulence of Eurasian T. gondii strains for mice has been correlated in genetic screens with the expression of polymorphic rhoptry kinases (ROP kinases) secreted into the host cell during infection. We show that the molecular targets of the virulent allelic form of ROP18 kinase are members of a family of cellular GTPases, the interferon-inducible IRG (immunity-related GTPase) proteins, known from earlier work to be essential resistance factors in mice against avirulent strains of T. gondii. Virulent T. gondii strain ROP18 kinase phosphorylates several mouse IRG proteins. We show that the parasite kinase phosphorylates host Irga6 at two threonines in the nucleotide-binding domain, biochemically inactivating the GTPase and inhibiting its accumulation and action at the T. gondii parasitophorous vacuole membrane. Our analysis identifies the conformationally active switch I region of the GTP-binding site as an Achilles' heel of the IRG protein pathogen-resistance mechanism. The polymorphism of ROP18 in natural T. gondii populations indicates the existence of a dynamic, rapidly evolving ecological relationship between parasite virulence factors and host resistance factors. This system should be unusually fruitful for analysis at both ecological and molecular levels since both T. gondii and the mouse are widespread and abundant in the wild and are well-established model species with excellent analytical tools available.

Truncating germline mutations in the tumor suppressor gene BRCA-1 associated protein-1 (BAP1) have been reported in families predisposed to developing a wide range of different cancer types including uveal melanoma and cutaneous melanoma. There has also been an association between amelanotic tumor development and germline BAP1 mutation suggesting a possible phenotypic characteristic of BAP1 mutation carriers. Though there have been many types of cancer associated with germline BAP1 mutation, the full spectrum of disease association is yet to be ascertained. Here we describe a Danish family with predominantly uveal melanoma but also a range of other tumor types including lung, neuroendocrine, stomach, and breast cancer; as well as pigmented skin lesions. Whole-exome sequencing identified a BAP1 splice mutation located at c.581-2A>G, which leads to a premature truncation of BAP1 in an individual with uveal melanoma. This mutation was carried by several other family members with melanoma or various cancers. The finding expands on the growing profile of BAP1 as an important uveal and cutaneous melanoma tumor suppressor gene and implicates its involvement in the development of lung, and stomach cancer.

Apicomplexan parasites such as Toxoplasma gondii rely on actin-based motility to cross biological barriers and invade host cells. Key structural and biochemical differences in host and parasite actins make this an attractive target for small-molecule inhibitors. Here we took advantage of recent advances in the synthesis of cyclic depsipeptide compounds that stabilize filamentous actin to test the ability of chondramides to disrupt growth of T. gondii in vitro. Structural modeling of chondramide A (2) binding to an actin filament model revealed variations in the binding site between host and parasite actins. A series of 10 previously synthesized analogues (2b-k) with substitutions in the β-tyrosine moiety blocked parasite growth on host cell monolayers with EC₅₀ values that ranged from 0.3 to 1.3 μM. In vitro polymerization assays using highly purified recombinant actin from T. gondii verified that synthetic and natural product chondramides target the actin cytoskeleton. Consistent with this, chondramide treatment blocked parasite invasion into host cells and was more rapidly effective than pyrimethamine, a standard therapeutic agent. Although the current compounds lack specificity for parasite vs host actin, these studies provide a platform for the future design and synthesis of synthetic cyclic peptide inhibitors that selectively disrupt actin dynamics in parasites.

The volatile anesthetic isoflurane protects against renal ischemia and reperfusion injury by releasing renal tubular TGF-β1. As adenosine is a powerful cytoprotective molecule, we tested whether TGF-β1 generated by isoflurane induces renal tubular ecto-5'-nucleotidase (CD73) and adenosine to protect against renal ischemia and reperfusion injury. Isoflurane induced new CD73 synthesis and increased adenosine generation in cultured kidney proximal tubule cells and in mouse kidney. Moreover, a TGF-β1-neutralizing antibody prevented isoflurane-mediated induction of CD73 activity. Mice anesthetized with isoflurane after renal ischemia and reperfusion had significantly reduced plasma creatinine and decreased renal tubular necrosis, neutrophil infiltration, and apoptosis compared with pentobarbital-anesthetized mice. Isoflurane failed to protect against renal ischemia and reperfusion injury in CD73-deficient mice, in mice pretreated with a selective CD73 inhibitor, or in mice treated with an adenosine receptor antagonist. The TGF-β1-neutralizing antibody or the CD73 inhibitor attenuated isoflurane-mediated protection against HK-2 cell apoptosis. Thus, isoflurane causes TGF-β1-dependent induction of renal tubular CD73 and adenosine generation to protect against renal ischemia and reperfusion injury. Modulation of this pathway may have important therapeutic implications to reduce morbidity and mortality arising from ischemic acute kidney injury.

Macrophages are specialized to detect and destroy intracellular microbes and yet a number of pathogens have evolved to exploit this hostile niche. Here we demonstrate that the obligate intracellular parasite Toxoplasma gondii disarms macrophage innate clearance mechanisms by secreting a serine threonine kinase called ROP18, which binds to and phosphorylates immunity-related GTPases (IRGs). Substrate profiling of ROP18 revealed a preference for a conserved motif within switch region I of the GTPase domain, a modification predicted to disrupt IRG function. Consistent with this, expression of ROP18 was both necessary and sufficient to block recruitment of Irgb6, which was in turn required for parasite destruction. ROP18 phosphorylation of IRGs prevented clearance within inflammatory monocytes and IFN-γ-activated macrophages, conferring parasite survival in vivo and promoting virulence. IRGs are implicated in clearance of a variety of intracellular pathogens, suggesting that other virulence factors may similarly thwart this innate cellular defense mechanism.

Liver failure due to ischemia and reperfusion (IR) and subsequent acute kidney injury are significant clinical problems. We showed previously that liver IR selectively reduced plasma sphinganine-1-phosphate levels without affecting sphingosine-1-phosphate (S1P) levels. Furthermore, exogenous sphinganine-1-phosphate protected against both liver and kidney injury induced by liver IR. In this study, we elucidated the signaling mechanisms of sphinganine-1-phosphate-mediated renal and hepatic protection. A selective S1P(1) receptor antagonist blocked the hepatic and renal protective effects of sphinganine-1-phosphate, whereas a selective S1P(2) or S1P(3) receptor antagonist was without effect. Moreover, a selective S1P(1) receptor agonist, SEW-2871, provided similar degree of liver and kidney protection compared with sphinganine-1-phosphate. Furthermore, in vivo gene knockdown of S1P(1) receptors with small interfering RNA abolished the hepatic and renal protective effects of sphinganine-1-phosphate. In contrast to sphinganine-1-phosphate, S1P's hepatic protection was enhanced with an S1P(3) receptor antagonist. Inhibition of extracellular signal-regulated kinase, Akt or pertussis toxin-sensitive G-proteins blocked sphinganine-1-phosphate-mediated liver and kidney protection in vivo. Taken together, our results show that sphinganine-1-phosphate provided renal and hepatic protection after liver IR injury in mice through selective activation of S1P(1) receptors and pertussis toxin-sensitive G-proteins with subsequent activation of ERK and Akt.

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.

The gastrointestinal tract hosts the largest collection of commensal microbes in the body. Infections at this site can cause significant perturbations in the microbiota, known as dysbiosis, that facilitate the expansion of pathobionts, and can elicit inappropriate immune responses that impair the intestinal barrier function. Dysbiosis typically occurs during intestinal infection with Toxoplasma gondii. Host resistance to T. gondii depends on a potent Th1 response. In addition, a Th17 response is also elicited. How Th17 cells contribute to the host response to T. gondii remains unclear. Here we show that class I-restricted T cell-associated molecule (CRTAM) expression on T cells is required for an optimal IL-17 production during T. gondii infection. Moreover, that the lack of IL-17, results in increased immunopathology caused by an impaired antimicrobial peptide production and bacterial translocation from the intestinal lumen to the mesenteric lymph nodes and spleen.

The protozoan parasite Toxoplasma gondii is thought to exploit monocyte trafficking to facilitate dissemination across endothelial barriers such as the blood-brain barrier. Here, we analysed the migration of parasitized monocytes in model endothelial and interstitial environments. We report that infection enhanced monocyte locomotion on the surface of endothelial cells, but profoundly inhibited monocyte transmigration across endothelial barriers. By contrast, infection robustly increased monocyte and macrophage migration through collagen-rich tissues in a Rho-ROCK-dependent manner consistent with integrin-independent interstitial migration. We further demonstrated that the secreted T. gondii protein kinase ROP17 was required for enhanced tissue migration. In vivo, ROP17-deficient parasites failed to upregulate monocyte tissue migration and exhibited an early dissemination delay, leading to prolonged mouse survival. Our findings indicate that the parasite-induced changes in monocyte motility primarily facilitate the transport of T. gondii through tissues and promote systemic dissemination, rather than shuttle parasites across the blood-brain barrier via extravasation.

Despite being a frequent cause of severe diarrheal disease in infants and an opportunistic infection in immunocompromised patients, Cryptosporidium research has lagged due to a lack of facile experimental methods. Here, we describe a platform for complete life cycle development and long-term growth of C. parvum in vitro using "air-liquid interface" (ALI) cultures derived from intestinal epithelial stem cells. Transcriptomic profiling revealed that differentiating epithelial cells grown under ALI conditions undergo profound changes in metabolism and development that enable completion of the parasite life cycle in vitro. ALI cultures support parasite expansion > 100-fold and generate viable oocysts that are transmissible in vitro and to mice, causing infection and animal death. Transgenic parasite lines created using CRISPR/Cas9 were used to complete a genetic cross in vitro, demonstrating Mendelian segregation of chromosomes during meiosis. ALI culture provides an accessible model that will enable innovative studies into Cryptosporidium biology and host interactions.

Toxoplasma motility is both activated and suppressed by 3',5'-cyclic nucleotide signaling. Cyclic GMP (cGMP) signaling through Toxoplasma gondii protein kinase G (TgPKG) activates motility, whereas cyclic AMP (cAMP) signaling through TgPKAc1 inhibits motility. Despite their importance, it remains unclear how cGMP and cAMP levels are maintained in Toxoplasma. Phosphodiesterases (PDEs) are known to inactivate cyclic nucleotides and are highly expanded in the Toxoplasma genome. Here, we analyzed the expression and function of the 18-member TgPDE family in tachyzoites, the virulent life stage of Toxoplasma. We detected the expression of 11 of 18 TgPDEs, confirming prior expression studies. A knockdown screen of the TgPDE family revealed four TgPDEs that contribute to lytic Toxoplasma growth (TgPDE1, TgPDE2, TgPDE5, and TgPDE9). Depletion of TgPDE1 or TgPDE2 caused severe growth defects, prompting further investigation. While TgPDE1 was important for extracellular motility, TgPDE2 was important for host cell invasion, parasite replication, host cell egress, and extracellular motility. TgPDE1 displayed a plasma membrane/cytomembranous distribution, whereas TgPDE2 displayed an endoplasmic reticulum/cytomembranous distribution. Biochemical analysis of TgPDE1 and TgPDE2 purified from Toxoplasma lysates revealed that TgPDE1 hydrolyzes both cGMP and cAMP, whereas TgPDE2 was cAMP specific. Interactome studies of TgPDE1 and TgPDE2 indicated that they do not physically interact with each other or other TgPDEs but may be regulated by kinases and proteases. Our studies have identified TgPDE1 and TgPDE2 as central regulators of tachyzoite cyclic nucleotide levels and enable future studies aimed at determining how these enzymes are regulated and cooperate to control Toxoplasma motility and growth. IMPORTANCE Apicomplexan parasites require motility to actively infect host cells and cause disease. Cyclic nucleotide signaling governs apicomplexan motility, but it is unclear how cyclic nucleotide levels are maintained in these parasites. In search of novel regulators of cyclic nucleotides in the model apicomplexan Toxoplasma, we identified and characterized two catalytically active phosphodiesterases, TgPDE1 and TgPDE2, that are important for Toxoplasma's virulent tachyzoite life cycle. Enzymes that generate, sense, or degrade cyclic nucleotides make attractive targets for therapies aimed at paralyzing and killing apicomplexan parasites.

Hair colour is a polygenic phenotype that results from differences in the amount and ratio of melanins located in the hair bulb. Genome-wide association studies (GWAS) have identified many loci involved in the pigmentation pathway affecting hair colour. However, most of the associated loci overlap non-protein coding regions and many of the molecular mechanisms underlying pigmentation variation are still not understood. Here, we conduct GWAS meta-analyses of hair colour in a Canadian cohort of 12,741 individuals of European ancestry. By performing fine-mapping analyses we identify candidate causal variants in pigmentation loci associated with blonde, red and brown hair colour. Additionally, we observe colocalization of several GWAS hits with expression and methylation quantitative trait loci (QTLs) of cultured melanocytes. Finally, transcriptome-wide association studies (TWAS) further nominate the expression of EDNRB and CDK10 as significantly associated with hair colour. Our results provide insights on the mechanisms regulating pigmentation biology in humans.

Although expression quantitative trait loci (eQTLs) have been powerful in identifying susceptibility genes from genome-wide association study (GWAS) findings, most trait-associated loci are not explained by eQTLs alone. Alternative QTLs, including DNA methylation QTLs (meQTLs), are emerging, but cell-type-specific meQTLs using cells of disease origin have been lacking. Here, we established an meQTL dataset by using primary melanocytes from 106 individuals and identified 1,497,502 significant cis-meQTLs. Multi-QTL colocalization with meQTLs, eQTLs, and mRNA splice-junction QTLs from the same individuals together with imputed methylome-wide and transcriptome-wide association studies identified candidate susceptibility genes at 63% of melanoma GWAS loci. Among the three molecular QTLs, meQTLs were the single largest contributor. To compare melanocyte meQTLs with those from malignant melanomas, we performed meQTL analysis on skin cutaneous melanomas from The Cancer Genome Atlas (n = 444). A substantial proportion of meQTL probes (45.9%) in primary melanocytes is preserved in melanomas, while a smaller fraction of eQTL genes is preserved (12.7%). Integration of melanocyte multi-QTLs and melanoma meQTLs identified candidate susceptibility genes at 72% of melanoma GWAS loci. Beyond GWAS annotation, meQTL-eQTL colocalization in melanocytes suggested that 841 unique genes potentially share a causal variant with a nearby methylation probe in melanocytes. Finally, melanocyte trans-meQTLs identified a hotspot for rs12203592, a cis-eQTL of a transcription factor, IRF4, with 131 candidate target CpGs. Motif enrichment and IRF4 ChIP-seq analysis demonstrated that these target CpGs are enriched in IRF4 binding sites, suggesting an IRF4-mediated regulatory network. Our study highlights the utility of cell-type-specific meQTLs.