Transcriptome organization of white blood cells through gene co-expression network analysis in a large RNA-seq dataset.

Gene co-expression network analysis enables identification of biologically meaningful clusters of co-regulated genes (modules) in an unsupervised manner. We present here the largest study conducted thus far of co-expression networks in white blood cells (WBC) based on RNA-seq data from 624 individuals. We identify 41 modules, 13 of them related to specific immune-related functions and cell types (e.g. neutrophils, B and T cells, NK cells, and plasmacytoid dendritic cells); we highlight biologically relevant lncRNAs for each annotated module of co-expressed genes. We further characterize with unprecedented resolution the modules in T cell sub-types, through the availability of 95 immune phenotypes obtained by flow cytometry in the same individuals. This study provides novel insights into the transcriptional architecture of human leukocytes, showing how network analysis can advance our understanding of coding and non-coding gene interactions in immune system cells.

The intestinal microbiome of inflammatory bowel disease across the pediatric age range.

Dysbiosis is associated with pediatric and adult-onset inflammatory bowel disease (IBD), but the role of dysbiosis and the microbiome in very early onset IBD (VEO-IBD) has not yet been described. Here, we aimed to demonstrate the impact of age and inflammation on microbial community structure using shotgun metagenomic sequencing in children with VEO-IBD, pediatric-onset IBD, and age-matched pediatric healthy controls (HC) observed longitudinally over the course of 8 weeks. We found disease-related differences in alpha and beta diversity between HC and children with IBD or VEO-IBD. Using a healthy microbial maturity index modeled from HC across the age range to characterize their gut microbiota, we found that children with pediatric-onset IBD and VEO-IBD had lower maturity than their age-matched HC groups, suggesting a disease effect on the microbial community. In addition, patients with pediatric IBD had significantly lower maturity than those with VEO-IBD, who had more heterogeneity at the youngest ages, highlighting differences in these two cohorts that were not captured in standard comparisons of alpha and beta diversity. These results demonstrate that young age and inflammation independently impact microbial community structure. However, the effect is not additive in the youngest patients, likely because of the heterogeneous and dynamic stool microbiome in this population.

Effect of Genetic Factors, Age and Sex on Levels of Circulating Extracellular Vesicles and Platelets.

Extracellular vesicles (EVs) mediate cell interactions in biological processes, such as receptor activation or molecule transfer. Estimates of variation by age and sex have been limited by small sample size, and no report has assessed the contribution of genetic factors to levels of EVs. Here, we evaluated blood levels of 25 EV and 3 platelet traits in 974 individuals (933 genotyped) and reported the first genome-wide association study (GWAS) on levels of these traits. EV levels all decreased with age, whereas the trend for their surface markers was more heterogeneous. Platelets and CD31dim platelet EVs significantly increased in females compared to males, although CD31 expression on both platelets and platelet EVs decreased in females. Levels of the other EV subsets were similar between sexes. GWAS revealed three statistically significant genetic signals associated with EV levels in the F10 and GBP1 genes and in the intergenic region between LRIG1 and KBTBD8. These add to a signal in the 3'UTR of RHOF associated with CD31 expression on platelets that was previously found to be associated with other platelet traits. These findings suggest that EV formation is not a simple, constant adjunct of metabolism but is under both age-related and genetic control that can be independent of the regulation of the levels of the cells from which the EVs derive.

Mutations in MYO9B are associated with Charcot-Marie-Tooth disease type 2 neuropathies and isolated optic atrophy.

BACKGROUND AND PURPOSE: Charcot-Marie-Tooth disease (CMT) is a heterogeneous group of disorders caused by mutations in at least 100 genes. However, approximately 60% of cases with axonal neuropathies (CMT2) still remain without a genetic diagnosis. We aimed at identifying novel disease genes responsible for CMT2. METHODS: We performed whole exome sequencing and targeted next generation sequencing panel analyses on a cohort of CMT2 families with evidence for autosomal recessive inheritance. We also performed functional studies to explore the pathogenetic role of selected variants. RESULTS: We identified rare, recessive variants in the MYO9B (myosin IX) gene in two families with CMT2. MYO9B has not yet been associated with a human disease. MYO9B is an unconventional single-headed processive myosin motor protein with signaling properties, and, consistent with this, our results indicate that a variant occurring in the MYO9B motor domain impairs protein expression level and motor activity. Interestingly, a Myo9b-null mouse has degenerating axons in sciatic nerves and optic nerves, indicating that MYO9B plays an essential role in both peripheral nervous system and central nervous system axons, respectively. The degeneration observed in the optic nerve prompted us to screen for MYO9B mutations in a cohort of patients with optic atrophy (OA). Consistent with this, we found compound heterozygous variants in one case with isolated OA. CONCLUSIONS: Novel or very rare variants in MYO9B are associated with CMT2 and isolated OA.

Cross-sectional analysis of the humoral response after SARS-CoV-2 vaccination in Sardinian multiple sclerosis patients, a follow-up study.

Monitoring immune responses to SARS-CoV-2 vaccination and its clinical efficacy over time in Multiple Sclerosis (MS) patients treated with disease-modifying therapies (DMTs) help to establish the optimal strategies to ensure adequate COVID-19 protection without compromising disease control offered by DMTs. Following our previous observations on the humoral response one month after two doses of BNT162b2 vaccine (T1) in MS patients differently treated, here we present a cross-sectional and longitudinal follow-up analysis six months following vaccination (T2, n=662) and one month following the first booster (T3, n=185). Consistent with results at T1, humoral responses were decreased in MS patients treated with fingolimod and anti-CD20 therapies compared with untreated patients also at the time points considered here (T2 and T3). Interestingly, a strong upregulation one month after the booster was observed in patients under every DMTs analyzed, including those treated with fingolimod and anti-CD20 therapies. Although patients taking these latter therapies had a higher rate of COVID-19 infection five months after the first booster, only mild symptoms that did not require hospitalization were reported for all the DMTs analyzed here. Based on these findings we anticipate that additional vaccine booster shots will likely further improve immune responses and COVID-19 protection in MS patients treated with any DMT.

Large-scale sequencing identifies multiple genes and rare variants associated with Crohn's disease susceptibility.

Genome-wide association studies (GWASs) have identified hundreds of loci associated with Crohn's disease (CD). However, as with all complex diseases, robust identification of the genes dysregulated by noncoding variants typically driving GWAS discoveries has been challenging. Here, to complement GWASs and better define actionable biological targets, we analyzed sequence data from more than 30,000 patients with CD and 80,000 population controls. We directly implicate ten genes in general onset CD for the first time to our knowledge via association to coding variation, four of which lie within established CD GWAS loci. In nine instances, a single coding variant is significantly associated, and in the tenth, ATG4C, we see additionally a significantly increased burden of very rare coding variants in CD cases. In addition to reiterating the central role of innate and adaptive immune cells as well as autophagy in CD pathogenesis, these newly associated genes highlight the emerging role of mesenchymal cells in the development and maintenance of intestinal inflammation.

Clinical and laboratory predictors of monogenic very early onset inflammatory bowel disease.

BACKGROUND: Inflammatory bowel disease (IBD) is a chronic inflammatory disease of the gastrointestinal tract. Treatment for patients who have a monogenic cause of their IBD, often the youngest children, known as very early onset IBD (VEO-IBD), can be different from standard treatment for polygenic cases. Yet, ascertainment of these patients is difficult. METHODS: We analyzed cases of VEO-IBD to understand the breadth of monogenic etiology and to identify clinical, laboratory, and flow cytometric correlates of this subpopulation. RESULTS: Genetic causes of very early onset inflammatory bowel disease are highly diverse ranging from pure epithelial defects to classic T cell defects. Flow cytometry, other than testing for chronic granulomatous disease, has a low sensitivity for monogenic etiologies. Poor growth was a clinical feature associated with monogenic causality. CONCLUSIONS: Genetic testing is, at this moment, the most robust method for the identification of monogenic cases of very early onset IBD.

Genetic Analysis in African American Children Supports Ancestry-Specific Neuroblastoma Susceptibility.

BACKGROUND: Neuroblastoma is rarer in African American (AA) children compared with American children of European descent. AA children affected with neuroblastoma, however, more frequently develop the high-risk form of the disease. METHODS: We have genotyped an AA cohort of 629 neuroblastoma cases (254 high-risk) and 2,990 controls to investigate genetic susceptibility to neuroblastoma in AAs. RESULTS: We confirmed the known neuroblastoma susceptibility gene BARD1 at genome-wide significance in the subset of high-risk cases. We also estimated local admixture across the autosomal genome in the AA cases and controls and detected a signal at 4q31.22 where cases show an increase in European ancestry. A region at 17p13.1 showed increased African ancestry in the subgroup of high-risk cases with respect to intermediate- and low-risk cases. Using results from our published European American (EA) genome-wide association study (GWAS), we found that a polygenic score that included all independent SNPs showed a highly significant association (P value = 1.8 × 10-73) and explained 19% of disease risk variance in an independent EA cohort. In contrast, the best fit polygenic score (P value = 3.2 × 10-11) in AAs included only 22 independent SNPs with association P value < 2.75 × 10-6 in the EA GWAS, and explained 2% of neuroblastoma risk variance. The significance of the polygenic score dropped rapidly with inclusion of additional SNPs. CONCLUSIONS: These findings suggest that several common variants contribute to risk of neuroblastoma in an ancestry-specific fashion. IMPACT: This work supports the need for GWAS to be performed in populations of all races and ethnicities.

Molecular Diagnostic Outcomes from 700 Cases: What Can We Learn from a Retrospective Analysis of Clinical Exome Sequencing?

Clinical exome sequencing (CES) aids in the diagnosis of rare genetic disorders. Herein, we report the molecular diagnostic yield and spectrum of genetic alterations contributing to disease in 700 pediatric cases analyzed at the Children's Hospital of Philadelphia. The overall diagnostic yield was 23%, with three cases having more than one molecular diagnosis and 2.6% having secondary/additional findings. A candidate gene finding was reported in another 8.4% of cases. The clinical indications with the highest diagnostic yield were neurodevelopmental disorders (including seizures), whereas immune- and oncology-related indications were negatively associated with molecular diagnosis. The rapid expansion of knowledge regarding the genome's role in human disease necessitates reanalysis of CES samples. To capture these new discoveries, a subset of cases (n = 240) underwent reanalysis, with an increase in diagnostic yield. We describe our experience reporting CES results in a pediatric setting, including reporting of secondary findings, reporting newly discovered genetic conditions, and revisiting negative test results. Finally, we highlight the challenges associated with implementing critical updates to the CES workflow. Although these updates are necessary, they demand an investment of time and resources from the laboratory. In summary, these data demonstrate the clinical utility of exome sequencing and reanalysis, while highlighting the critical considerations for continuous improvement of a CES test in a clinical laboratory.

Effect of Different Disease-Modifying Therapies on Humoral Response to BNT162b2 Vaccine in Sardinian Multiple Sclerosis Patients.

Objectives: Vaccination against COVID-19 is highly recommended to patients affected by multiple sclerosis (MS); however, the impact of MS disease-modifying therapies (DMTs) on the immune response following vaccination has been only partially investigated. Here, we aimed to elucidate the effect of DMTs on the humoral immune response to mRNA-based anti-SARS-CoV-2 vaccines in MS patients. Methods: We obtained sera from 912 Sardinian MS patients and 63 healthy controls 30 days after the second dose of BNT162b2 vaccine and tested them for SARS-CoV-2 response using anti-Spike (S) protein-based serology. Previous SARS-CoV-2 infection was assessed by anti-Nucleocapsid (N) serology. Patients were either untreated or undergoing treatment with a total of 13 different DMTs. Differences between treatment groups comprised of at least 10 patients were assessed by generalized linear mixed-effects model. Demographic and clinical data and smoking status were analyzed as additional factors potentially influencing humoral immunity from COVID-19 vaccine. Results: MS patients treated with natalizumab, teriflunomide, azathioprine, fingolimod, ocrelizumab, and rituximab showed significantly lower humoral responses compared to untreated patients. We did not observe a statistically significant difference in response between patients treated with the other drugs (dimethyl fumarate, interferon, alemtuzumab and glatiramer acetate) and untreated patients. In addition, older age, male sex and active smoking were significantly associated with lower antibody titers against SARS-CoV-2. MS patients previously infected with SARS-CoV-2 had significantly higher humoral responses to vaccine than uninfected patients. Conclusion: Humoral response to BNT162b2 is significantly influenced by the specific DMTs followed by patients, as well as by other factors such as previous SARS-CoV-2 infection, age, sex, and smoking status. These results are important to inform targeted strategies to prevent clinically relevant COVID-19 in MS patients.

Immune Dysregulation in Human ITCH Deficiency Successfully Treated with Hematopoietic Cell Transplantation.

BACKGROUND: Mutations in ITCH, which encodes an E3 ubiquitin-protein ligase, can result in systemic autoimmunity and immunodeficiency. The clinical phenotype and mechanism of disease have not been fully characterized, resulting in a paucity of therapeutic options for this potentially fatal disease. OBJECTIVE: We aimed to (1) expand the understanding about the phenotype of human ITCH deficiency (2) further characterize the associated immune dysregulation, and (3) report the first successful hematopoietic cell transplant (HCT) in a patient with ITCH deficiency. METHODS: Disease profiling was performed in a patient with multisystem immune dysregulation. Whole exome sequencing with trio analysis and functional validation of candidate disease variants were performed, including mRNA and protein expression. Analyses to further delineate the immunophenotype included quantitative evaluation of lymphoid and myeloid subsets with flow cytometry and mass cytometry. RESULTS: A patient with multisystem immune dysregulation presenting with growth failure, very-early-onset inflammatory bowel disease, arthritis, uveitis, psoriasis, and type 1 diabetes mellitus underwent whole exome sequencing, which identified novel compound heterozygous mutations in ITCH. Reduced expression of ITCH mRNA and absent ITCH protein were found. Abnormalities in both lymphoid and myeloid lineages were identified. The patient underwent HCT. He demonstrated excellent immune reconstitution and resolution of many manifestations of his systemic disease. CONCLUSIONS: Here we report ITCH deficiency with unique clinical features of colonic very-early-onset inflammatory bowel disease, arthritis, and uveitis in the setting of immune dysregulation and further characterize the underlying immune dysregulation. We demonstrate that HCT can be an effective, and potentially curative, therapy for ITCH deficiency.

Variants in STXBP3 are Associated with Very Early Onset Inflammatory Bowel Disease, Bilateral Sensorineural Hearing Loss and Immune Dysregulation.

BACKGROUND AND AIMS: Very early onset inflammatory bowel disease [VEOIBD] is characterized by intestinal inflammation affecting infants and children less than 6 years of age. To date, over 60 monogenic aetiologies of VEOIBD have been identified, many characterized by highly penetrant recessive or dominant variants in underlying immune and/or epithelial pathways. We sought to identify the genetic cause of VEOIBD in a subset of patients with a unique clinical presentation. METHODS: Whole exome sequencing was performed on five families with ten patients who presented with a similar constellation of symptoms including medically refractory infantile-onset IBD, bilateral sensorineural hearing loss and, in the majority, recurrent infections. Genetic aetiologies of VEOIBD were assessed and Sanger sequencing was performed to confirm novel genetic findings. Western analysis on peripheral blood mononuclear cells and functional studies with epithelial cell lines were employed. RESULTS: In each of the ten patients, we identified damaging heterozygous or biallelic variants in the Syntaxin-Binding Protein 3 gene [STXBP3], a protein known to regulate intracellular vesicular trafficking in the syntaxin-binding protein family of molecules, but not associated to date with either VEOIBD or sensorineural hearing loss. These mutations interfere with either intron splicing or protein stability and lead to reduced STXBP3 protein expression. Knock-down of STXBP3 in CaCo2 cells resulted in defects in cell polarity. CONCLUSION: Overall, we describe a novel genetic syndrome and identify a critical role for STXBP3 in VEOIBD, sensorineural hearing loss and immune dysregulation.

Colonoids From Patients With Pediatric Inflammatory Bowel Disease Exhibit Decreased Growth Associated With Inflammation Severity and Durable Upregulation of Antigen Presentation Genes.

BACKGROUND: Defining epithelial cell contributions to inflammatory bowel disease (IBD) is essential for the development of much needed therapies for barrier repair. Children with very early onset (VEO)-IBD have more extensive, severe, and refractory disease than older children and adults with IBD and, in some cases, have defective barrier function. We therefore evaluated functional and transcriptomic differences between pediatric IBD (VEO and older onset) and non-IBD epithelium using 3-dimensional, biopsy-derived organoids. METHODS: We measured growth efficiency relative to histopathological and clinical parameters in patient enteroid (ileum) and colonoid (colon) lines. We performed RNA-sequencing on patient colonoids and subsequent flow cytometry after multiple passages to evaluate changes that persisted in culture. RESULTS: Enteroids and colonoids from pediatric patients with IBD exhibited decreased growth associated with histological inflammation compared with non-IBD controls. We observed increased LYZ expression in colonoids from pediatric IBD patients, which has been reported previously in adult patients with IBD. We also observed upregulation of antigen presentation genes HLA-DRB1 and HLA-DRA, which persisted after prolonged passaging in patients with pediatric IBD. CONCLUSIONS: We present the first functional evaluation of enteroids and colonoids from patients with VEO-IBD and older onset pediatric IBD, a subset of which exhibits poor growth. Enhanced, persistent epithelial antigen presentation gene expression in patient colonoids supports the notion that epithelial cell-intrinsic differences may contribute to IBD pathogenesis.

Author Correction: Complex genetic signatures in immune cells underlie autoimmunity and inform therapy.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Complex genetic signatures in immune cells underlie autoimmunity and inform therapy.

We report on the influence of ~22 million variants on 731 immune cell traits in a cohort of 3,757 Sardinians. We detected 122 significant (P < 1.28 × 10-11) independent association signals for 459 cell traits at 70 loci (53 of them novel) identifying several molecules and mechanisms involved in cell regulation. Furthermore, 53 signals at 36 loci overlapped with previously reported disease-associated signals, predominantly for autoimmune disorders, highlighting intermediate phenotypes in pathogenesis. Collectively, our findings illustrate complex genetic regulation of immune cells with highly selective effects on autoimmune disease risk at the cell-subtype level. These results identify drug-targetable pathways informing the design of more specific treatments for autoimmune diseases.

Impaired Redox and Protein Homeostasis as Risk Factors and Therapeutic Targets in Toxin-Induced Biliary Atresia.

BACKGROUND & AIMS: Extrahepatic biliary atresia (BA) is a pediatric liver disease with no approved medical therapy. Recent studies using human samples and experimental modeling suggest that glutathione redox metabolism and heterogeneity play a role in disease pathogenesis. We sought to dissect the mechanistic basis of liver redox variation and explore how other stress responses affect cholangiocyte injury in BA. METHODS: We performed quantitative in situ hepatic glutathione redox mapping in zebrafish larvae carrying targeted mutations in glutathione metabolism genes and correlated these findings with sensitivity to the plant-derived BA-linked toxin biliatresone. We also determined whether genetic disruption of HSP90 protein quality control pathway genes implicated in human BA altered biliatresone toxicity in zebrafish and human cholangiocytes. An in vivo screening of a known drug library was performed to identify novel modifiers of cholangiocyte injury in the zebrafish experimental BA model, with subsequent validation. RESULTS: Glutathione metabolism gene mutations caused regionally distinct changes in the redox potential of cholangiocytes that differentially sensitized them to biliatresone. Disruption of human BA-implicated HSP90 pathway genes sensitized zebrafish and human cholangiocytes to biliatresone-induced injury independent of glutathione. Phosphodiesterase-5 inhibitors and other cyclic guanosine monophosphate signaling activators worked synergistically with the glutathione precursor N-acetylcysteine in preventing biliatresone-induced injury in zebrafish and human cholangiocytes. Phosphodiesterase-5 inhibitors enhanced proteasomal degradation and required intact HSP90 chaperone. CONCLUSION: Regional variation in glutathione metabolism underlies sensitivity to the biliary toxin biliatresone and may account for the reported association between BA transplant-free survival and glutathione metabolism gene expression. Human BA can be causatively linked to genetic modulation of protein quality control. Combined treatment with N-acetylcysteine and cyclic guanosine monophosphate signaling enhancers warrants further investigation as therapy for BA.

Exome Sequencing in Individuals with Isolated Biliary Atresia.

Biliary atresia (BA) is a severe pediatric liver disease resulting in necroinflammatory obliteration of the extrahepatic biliary tree. BA presents within the first few months of life as either an isolated finding or with additional syndromic features. The etiology of isolated BA is unknown, with evidence for infectious, environmental, and genetic risk factors described. However, to date, there are no definitive causal genes identified for isolated BA in humans, and the question of whether single gene defects play a major role remains open. We performed exome-sequencing in 101 North American patients of European descent with isolated BA (including 30 parent-child trios) and considered several experimental designs to identify potentially deleterious protein-altering variants that may be involved in the disease. In a case-only analysis, we did not identify genes with variants shared among more than two probands, and burden tests of rare variants using a case-case control design did not yield significant results. In the trio analysis of 30 simplex families (patient and parent trios), we identified 66 de novo variants in 66 genes including potentially deleterious variants in STIP1 and REV1. STIP1 is a co-chaperone for the heat-shock protein, HSP90, and has been shown to have diverse functions in yeast, flies and mammals, including stress-responses. REV1 is known to be a key player in DNA repair pathway and to interact with HSP90. In conclusion, our results do not support the hypothesis that a simple genetic model is responsible for the majority of cases of isolated BA. Our finding of de novo variants in genes linked to evolutionarily conserved stress responses (STIP1 and REV1) suggests that exploration of how genetic susceptibility and environmental exposure may interact to cause BA is warranted.

Rare copy number variants in over 100,000 European ancestry subjects reveal multiple disease associations.

Copy number variants (CNVs) are suggested to have a widespread impact on the human genome and phenotypes. To understand the role of CNVs across human diseases, we examine the CNV genomic landscape of 100,028 unrelated individuals of European ancestry, using SNP and CGH array datasets. We observe an average CNV burden of ~650 kb, identifying a total of 11,314 deletion, 5625 duplication, and 2746 homozygous deletion CNV regions (CNVRs). In all, 13.7% are unreported, 58.6% overlap with at least one gene, and 32.8% interrupt coding exons. These CNVRs are significantly more likely to overlap OMIM genes (2.94-fold), GWAS loci (1.52-fold), and non-coding RNAs (1.44-fold), compared with random distribution (P < 1 × 10-3). We uncover CNV associations with four major disease categories, including autoimmune, cardio-metabolic, oncologic, and neurological/psychiatric diseases, and identify several drug-repurposing opportunities. Our results demonstrate robust frequency definition for large-scale rare variant association studies, identify CNVs associated with major disease categories, and illustrate the pleiotropic impact of CNVs in human disease.

The Unique Disease Course of Children with Very Early onset-Inflammatory Bowel Disease.

BACKGROUND: Insight into the pathogenesis of very early onset-inflammatory bowel disease (VEO-IBD) has expanded through the identification of causative monogenic defects detected in a subset of patients. However, the clinical course of this population remains uncertain. The study objective is to determine whether VEO-IBD is associated with more severe disease, defined as increased surgical intervention and growth failure, than older pediatric IBD. Secondary outcomes included therapeutic response and hospitalizations. METHODS: Subjects with IBD diagnosed younger than 6 years old (VEO-IBD) were compared with children diagnosed 6 to 10 (intermediate-onset) and older than 10 years of age (older-onset IBD). Metadata obtained from the medical record included age of onset, disease phenotype and location, surgeries, medical therapy, and comorbid conditions. Length of follow-up was at least 1 year from diagnosis. RESULTS: There were 229, 221, and 521 subjects with VEO, intermediate-onset, and older-onset IBD, respectively. Very early onset-inflammatory bowel disease subjects underwent more diverting ileostomies (P < 0.001) and colectomies (P < 0.001) than the older children. There was less improvement in weight- and height-for-age Z scores during the follow-up period in subjects with VEO-IBD. Additionally, subjects with VEO-IBD had higher rates of medication failure at 1 year and were more frequently readmitted to the hospital. Targeted therapy was successfully used almost exclusively in VEO-IBD. CONCLUSION: Patients with VEO-IBD can have a more severe disease course with increased surgical interventions and poor growth as compared with older-onset IBD patients. Further, VEO-IBD patients are more likely to be refractory to conventional therapies. Strategies using targeted therapy in these children can improve outcome and, in some cases, be curative.

Neural crest-derived tumor neuroblastoma and melanoma share 1p13.2 as susceptibility locus that shows a long-range interaction with the SLC16A1 gene.

Neuroblastoma (NB) and malignant cutaneous melanoma (CMM) are neural crest cells (NCC)-derived tumors and may have a shared genetic basis, but this has not been investigated systematically by genome-wide association studies (GWAS). We took a three-staged approach to conduct cross-disease meta-analysis of GWAS for NB and CMM (2101 NB cases and 4202 controls; 12 874 CMM cases and 23 203 controls) to identify shared loci. Findings were replicated in 1403 NB cases and 1403 controls of European ancestry and in 636 NB, 508 CMM cases and 2066 controls of Italian origin. We found a cross-association at locus 1p13.2 (rs2153977, odds ratio = 0.91, P = 5.36 × 10-8). We also detected a suggestive (P < 10-7) NB-CMM cross-association at 2q37.1 with opposite effect on cancer risk. Pathway analysis of 110 NB-CMM risk loci with P < 10-4 demonstrated enrichment of biological processes such as cell migration, cell cycle, metabolism and immune response, which are essential of human NCC development, underlying both tumors. In vitro and in silico analyses indicated that the rs2153977-T protective allele, located in an NB and CMM enhancer, decreased expression of SLC16A1 via long-range loop formation and altered a T-box protein binding site. Upon depletion of SLC16A1, we observed a decrease of cellular proliferation and invasion in both NB and CMM cell lines, suggesting its role as oncogene. This is the largest study to date examining pleiotropy across two NC cell-derived tumors identifying 1p13.2 as common susceptibility locus for NB and CMM risk. We demonstrate that combining genome-wide association studies results across cancers with same origins can identify new loci common to neuroblastoma and melanoma arising from tissues which originate from neural crest cells. Our results also show 1p13.2 confer risk to neuroblastoma and melanoma by regulating SLC16A1.