Multitrait genetic association analysis identifies 50 new risk loci for gastro-oesophageal reflux, seven new loci for Barrett's oesophagus and provides insights into clinical heterogeneity in reflux diagnosis.

OBJECTIVE: Gastro-oesophageal reflux disease (GERD) has heterogeneous aetiology primarily attributable to its symptom-based definitions. GERD genome-wide association studies (GWASs) have shown strong genetic overlaps with established risk factors such as obesity and depression. We hypothesised that the shared genetic architecture between GERD and these risk factors can be leveraged to (1) identify new GERD and Barrett's oesophagus (BE) risk loci and (2) explore potentially heterogeneous pathways leading to GERD and oesophageal complications. DESIGN: We applied multitrait GWAS models combining GERD (78 707 cases; 288 734 controls) and genetically correlated traits including education attainment, depression and body mass index. We also used multitrait analysis to identify BE risk loci. Top hits were replicated in 23andMe (462 753 GERD cases, 24 099 BE cases, 1 484 025 controls). We additionally dissected the GERD loci into obesity-driven and depression-driven subgroups. These subgroups were investigated to determine how they relate to tissue-specific gene expression and to risk of serious oesophageal disease (BE and/or oesophageal adenocarcinoma, EA). RESULTS: We identified 88 loci associated with GERD, with 59 replicating in 23andMe after multiple testing corrections. Our BE analysis identified seven novel loci. Additionally we showed that only the obesity-driven GERD loci (but not the depression-driven loci) were associated with genes enriched in oesophageal tissues and successfully predicted BE/EA. CONCLUSION: Our multitrait model identified many novel risk loci for GERD and BE. We present strong evidence for a genetic underpinning of disease heterogeneity in GERD and show that GERD loci associated with depressive symptoms are not strong predictors of BE/EA relative to obesity-driven GERD loci.

Identification of novel risk loci, causal insights, and heritable risk for Parkinson's disease: a meta-analysis of genome-wide association studies.

BACKGROUND: Genome-wide association studies (GWAS) in Parkinson's disease have increased the scope of biological knowledge about the disease over the past decade. We aimed to use the largest aggregate of GWAS data to identify novel risk loci and gain further insight into the causes of Parkinson's disease. METHODS: We did a meta-analysis of 17 datasets from Parkinson's disease GWAS available from European ancestry samples to nominate novel loci for disease risk. These datasets incorporated all available data. We then used these data to estimate heritable risk and develop predictive models of this heritability. We also used large gene expression and methylation resources to examine possible functional consequences as well as tissue, cell type, and biological pathway enrichments for the identified risk factors. Additionally, we examined shared genetic risk between Parkinson's disease and other phenotypes of interest via genetic correlations followed by Mendelian randomisation. FINDINGS: Between Oct 1, 2017, and Aug 9, 2018, we analysed 7·8 million single nucleotide polymorphisms in 37 688 cases, 18 618 UK Biobank proxy-cases (ie, individuals who do not have Parkinson's disease but have a first degree relative that does), and 1·4 million controls. We identified 90 independent genome-wide significant risk signals across 78 genomic regions, including 38 novel independent risk signals in 37 loci. These 90 variants explained 16-36% of the heritable risk of Parkinson's disease depending on prevalence. Integrating methylation and expression data within a Mendelian randomisation framework identified putatively associated genes at 70 risk signals underlying GWAS loci for follow-up functional studies. Tissue-specific expression enrichment analyses suggested Parkinson's disease loci were heavily brain-enriched, with specific neuronal cell types being implicated from single cell data. We found significant genetic correlations with brain volumes (false discovery rate-adjusted p=0·0035 for intracranial volume, p=0·024 for putamen volume), smoking status (p=0·024), and educational attainment (p=0·038). Mendelian randomisation between cognitive performance and Parkinson's disease risk showed a robust association (p=8·00 × 10-7). INTERPRETATION: These data provide the most comprehensive survey of genetic risk within Parkinson's disease to date, to the best of our knowledge, by revealing many additional Parkinson's disease risk loci, providing a biological context for these risk factors, and showing that a considerable genetic component of this disease remains unidentified. These associations derived from European ancestry datasets will need to be followed-up with more diverse data. FUNDING: The National Institute on Aging at the National Institutes of Health (USA), The Michael J Fox Foundation, and The Parkinson's Foundation (see appendix for full list of funding sources).

The Parkinson's Disease Mendelian Randomization Research Portal.

BACKGROUND: Mendelian randomization is a method for exploring observational associations to find evidence of causality. OBJECTIVE: To apply Mendelian randomization between risk factors/phenotypic traits (exposures) and PD in a large, unbiased manner, and to create a public resource for research. METHODS: We used two-sample Mendelian randomization in which the summary statistics relating to single-nucleotide polymorphisms from 5,839 genome-wide association studies of exposures were used to assess causal relationships with PD. We selected the highest-quality exposure genome-wide association studies for this report (n = 401). For the disease outcome, summary statistics from the largest published PD genome-wide association studies were used. For each exposure, the causal effect on PD was assessed using the inverse variance weighted method, followed by a range of sensitivity analyses. We used a false discovery rate of 5% from the inverse variance weighted analysis to prioritize exposures of interest. RESULTS: We observed evidence for causal associations between 12 exposures and risk of PD. Of these, nine were effects related to increasing adiposity and decreasing risk of PD. The remaining top three exposures that affected PD risk were tea drinking, time spent watching television, and forced vital capacity, but these may have been biased and were less convincing. Other exposures at nominal statistical significance included inverse effects of smoking and alcohol. CONCLUSIONS: We present a new platform which offers Mendelian randomization analyses for a total of 5,839 genome-wide association studies versus the largest PD genome-wide association studies available (https://pdgenetics.shinyapps.io/MRportal/). Alongside, we report further evidence to support a causal role for adiposity on lowering the risk of PD. © 2019 The Authors. Movement Disorders published by Wiley Periodicals, Inc. on behalf of International Parkinson and Movement Disorder Society.

Anti-IL-13Rα2 therapy promotes recovery in a murine model of inflammatory bowel disease.

There continues to be a major need for more effective inflammatory bowel disease (IBD) therapies. IL-13Rα2 is a decoy receptor that binds the cytokine IL-13 with high affinity and diminishes its STAT6-mediated effector functions. Previously, we found that IL-13Rα2 was necessary for IBD in mice deficient in the anti-inflammatory cytokine IL-10. Here, we tested for the first time a therapeutic antibody specifically targeting IL-13Rα2. We also used the antibody and Il13ra2-/- mice to dissect the role of IL-13Rα2 in IBD pathogenesis and recovery. Il13ra2-/- mice were modestly protected from induction of dextran sodium sulfate (DSS)-induced colitis. Following a 7-day recovery period, Il13ra2-/- mice or wild-type mice administered the IL-13Rα2-neutralizing antibody had significantly improved colon health compared to control mice. Neutralizing IL-13Rα2 to increase IL-13 bioavailability promoted resolution of IBD even if neutralization occurred only during recovery. To link our observations in mice to a large human cohort, we conducted a phenome-wide association study of a more active variant of IL-13 (R130Q) that has reduced affinity for IL-13Rα2. Human subjects carrying R130Q reported a lower risk for Crohn's disease. Our findings endorse moving anti-IL-13Rα2 into preclinical drug development with the goal of accelerating recovery and maintaining remission in Crohn's disease patients.

Genetic Architectures of Childhood- and Adult-Onset Asthma Are Partly Distinct.

The extent to which genetic risk factors are shared between childhood-onset (COA) and adult-onset (AOA) asthma has not been estimated. On the basis of data from the UK Biobank study (n = 447,628), we found that the variance in disease liability explained by common variants is higher for COA (onset at ages between 0 and 19 years; h2g = 25.6%) than for AOA (onset at ages between 20 and 60 years; h2g = 10.6%). The genetic correlation (rg) between COA and AOA was 0.67. Variation in age of onset among COA-affected individuals had a low heritability (h2g = 5%), which we confirmed in independent studies and also among AOA-affected individuals. To identify subtype-specific genetic associations, we performed a genome-wide association study (GWAS) in the UK Biobank for COA (13,962 affected individuals) and a separate GWAS for AOA (26,582 affected individuals) by using a common set of 300,671 controls for both studies. We identified 123 independent associations for COA and 56 for AOA (37 overlapped); of these, 98 and 34, respectively, were reproducible in an independent study (n = 262,767). Collectively, 28 associations were not previously reported. For 96 COA-associated variants, including five variants that represent COA-specific risk factors, the risk allele was more common in COA- than in AOA-affected individuals. Conversely, we identified three variants that are stronger risk factors for AOA. Variants associated with obesity and smoking had a stronger contribution to the risk of AOA than to the risk of COA. Lastly, we identified 109 likely target genes of the associated variants, primarily on the basis of correlated expression quantitative trait loci (up to n = 31,684). GWAS informed by age of onset can identify subtype-specific risk variants, which can help us understand differences in pathophysiology between COA and AOA and so can be informative for drug development.

Genome-wide association analyses of chronotype in 697,828 individuals provides insights into circadian rhythms.

Being a morning person is a behavioural indicator of a person's underlying circadian rhythm. Using genome-wide data from 697,828 UK Biobank and 23andMe participants we increase the number of genetic loci associated with being a morning person from 24 to 351. Using data from 85,760 individuals with activity-monitor derived measures of sleep timing we find that the chronotype loci associate with sleep timing: the mean sleep timing of the 5% of individuals carrying the most morningness alleles is 25 min earlier than the 5% carrying the fewest. The loci are enriched for genes involved in circadian regulation, cAMP, glutamate and insulin signalling pathways, and those expressed in the retina, hindbrain, hypothalamus, and pituitary. Using Mendelian Randomisation, we show that being a morning person is causally associated with better mental health but does not affect BMI or risk of Type 2 diabetes. This study offers insights into circadian biology and its links to disease in humans.

Phenome-wide association studies across large population cohorts support drug target validation.

Phenome-wide association studies (PheWAS) have been proposed as a possible aid in drug development through elucidating mechanisms of action, identifying alternative indications, or predicting adverse drug events (ADEs). Here, we select 25 single nucleotide polymorphisms (SNPs) linked through genome-wide association studies (GWAS) to 19 candidate drug targets for common disease indications. We interrogate these SNPs by PheWAS in four large cohorts with extensive health information (23andMe, UK Biobank, FINRISK, CHOP) for association with 1683 binary endpoints in up to 697,815 individuals and conduct meta-analyses for 145 mapped disease endpoints. Our analyses replicate 75% of known GWAS associations (P < 0.05) and identify nine study-wide significant novel associations (of 71 with FDR < 0.1). We describe associations that may predict ADEs, e.g., acne, high cholesterol, gout, and gallstones with rs738409 (p.I148M) in PNPLA3 and asthma with rs1990760 (p.T946A) in IFIH1. Our results demonstrate PheWAS as a powerful addition to the toolkit for drug discovery.

Gene expression imputation identifies candidate genes and susceptibility loci associated with cutaneous squamous cell carcinoma.

Cutaneous squamous cell carcinoma (cSCC) is a common skin cancer with genetic susceptibility loci identified in recent genome-wide association studies (GWAS). Transcriptome-wide association studies (TWAS) using imputed gene expression levels can identify additional gene-level associations. Here we impute gene expression levels in 6891 cSCC cases and 54,566 controls in the Kaiser Permanente Genetic Epidemiology Research in Adult Health and Aging (GERA) cohort and 25,558 self-reported cSCC cases and 673,788 controls from 23andMe. In a discovery-validation study, we identify 19 loci containing 33 genes whose imputed expression levels are associated with cSCC at false discovery rate < 10% in the GERA cohort and validate 15 of these candidate genes at Bonferroni significance in the 23andMe dataset, including eight genes in five novel susceptibility loci and seven genes in four previously associated loci. These results suggest genetic mechanisms contributing to cSCC risk and illustrate advantages and disadvantages of TWAS as a supplement to traditional GWAS analyses.

Whole exome sequencing reveals HSPA1L as a genetic risk factor for spontaneous preterm birth.

Preterm birth is a leading cause of morbidity and mortality in infants. Genetic and environmental factors play a role in the susceptibility to preterm birth, but despite many investigations, the genetic basis for preterm birth remain largely unknown. Our objective was to identify rare, possibly damaging, nucleotide variants in mothers from families with recurrent spontaneous preterm births (SPTB). DNA samples from 17 Finnish mothers who delivered at least one infant preterm were subjected to whole exome sequencing. All mothers were of northern Finnish origin and were from seven multiplex families. Additional replication samples of European origin consisted of 93 Danish sister pairs (and two sister triads), all with a history of a preterm delivery. Rare exonic variants (frequency <1%) were analyzed to identify genes and pathways likely to affect SPTB susceptibility. We identified rare, possibly damaging, variants in genes that were common to multiple affected individuals. The glucocorticoid receptor signaling pathway was the most significant (p<1.7e-8) with genes containing these variants in a subgroup of ten Finnish mothers, each having had 2-4 SPTBs. This pathway was replicated among the Danish sister pairs. A gene in this pathway, heat shock protein family A (Hsp70) member 1 like (HSPA1L), contains two likely damaging missense alleles that were found in four different Finnish families. One of the variants (rs34620296) had a higher frequency in cases compared to controls (0.0025 vs. 0.0010, p = 0.002) in a large preterm birth genome-wide association study (GWAS) consisting of mothers of general European ancestry. Sister pairs in replication samples also shared rare, likely damaging HSPA1L variants. Furthermore, in silico analysis predicted an additional phosphorylation site generated by rs34620296 that could potentially affect chaperone activity or HSPA1L protein stability. Finally, in vitro functional experiment showed a link between HSPA1L activity and decidualization. In conclusion, rare, likely damaging, variants in HSPA1L were observed in multiple families with recurrent SPTB.

Assessment of rosacea symptom severity by genome-wide association study and expression analysis highlights immuno-inflammatory and skin pigmentation genes.

Rosacea is a common, chronic skin disease of variable severity with limited treatment options. The cause of rosacea is unknown, but it is believed to be due to a combination of hereditary and environmental factors. Little is known about the genetics of the disease. We performed a genome-wide association study (GWAS) of rosacea symptom severity with data from 73 265 research participants of European ancestry from the 23andMe customer base. Seven loci had variants associated with rosacea at the genome-wide significance level (P < 5 × 10-8). Further analyses highlighted likely gene regions or effector genes including IRF4 (P = 1.5 × 10-17), a human leukocyte antigen (HLA) region flanked by PSMB9 and HLA-DMB (P = 2.2 × 10-15), HERC2-OCA2 (P = 4.2 × 10-12), SLC45A2 (P = 1.7 × 10-10), IL13 (P = 2.8 × 10-9), a region flanked by NRXN3 and DIO2 (P = 4.1 × 10-9), and a region flanked by OVOL1and SNX32 (P = 1.2 × 10-8). All associations with rosacea were novel except for the HLA locus. Two of these loci (HERC-OCA2 and SLC45A2) and another precedented variant (rs1805007 in melanocortin 1 receptor) with an association P value just below the significance threshold (P = 1.3 × 10-7) have been previously associated with skin phenotypes and pigmentation, two of these loci are linked to immuno-inflammation phenotypes (IL13 and PSMB9-HLA-DMA) and one has been associated with both categories (IRF4). Genes within three loci (PSMB9-HLA-DMA, HERC-OCA2 and NRX3-DIO2) were differentially expressed in a previously published clinical rosacea transcriptomics study that compared lesional to non-lesional samples. The identified loci provide specificity of inflammatory mechanisms in rosacea, and identify potential pathways for therapeutic intervention.

Multiethnic GWAS Reveals Polygenic Architecture of Earlobe Attachment.

The genetic basis of earlobe attachment has been a matter of debate since the early 20th century, such that geneticists argue both for and against polygenic inheritance. Recent genetic studies have identified a few loci associated with the trait, but large-scale analyses are still lacking. Here, we performed a genome-wide association study of lobe attachment in a multiethnic sample of 74,660 individuals from four cohorts (three with the trait scored by an expert rater and one with the trait self-reported). Meta-analysis of the three expert-rater-scored cohorts revealed six associated loci harboring numerous candidate genes, including EDAR, SP5, MRPS22, ADGRG6 (GPR126), KIAA1217, and PAX9. The large self-reported 23andMe cohort recapitulated each of these six loci. Moreover, meta-analysis across all four cohorts revealed a total of 49 significant (p < 5 × 10-8) loci. Annotation and enrichment analyses of these 49 loci showed strong evidence of genes involved in ear development and syndromes with auricular phenotypes. RNA sequencing data from both human fetal ear and mouse second branchial arch tissue confirmed that genes located among associated loci showed evidence of expression. These results provide strong evidence for the polygenic nature of earlobe attachment and offer insights into the biological basis of normal and abnormal ear development.

Erratum to: A multi-stage genome-wide association study of uterine fibroids in African Americans.

The article "A multi-stage genome-wide association study of uterine fibroids in African Americans", written by Jacklyn N. Hellwege, was originally published Online First without open access. After publication in volume 136, issue 10, page 1363-1373 the author decided to opt for Open Choice and to make the article an open access publication. Therefore, the copyright of the article has been changed to

Genetic Associations with Gestational Duration and Spontaneous Preterm Birth.

BACKGROUND: Despite evidence that genetic factors contribute to the duration of gestation and the risk of preterm birth, robust associations with genetic variants have not been identified. We used large data sets that included the gestational duration to determine possible genetic associations. METHODS: We performed a genomewide association study in a discovery set of samples obtained from 43,568 women of European ancestry using gestational duration as a continuous trait and term or preterm (<37 weeks) birth as a dichotomous outcome. We used samples from three Nordic data sets (involving a total of 8643 women) to test for replication of genomic loci that had significant genomewide association (P<5.0×10-8) or an association with suggestive significance (P<1.0×10-6) in the discovery set. RESULTS: In the discovery and replication data sets, four loci (EBF1, EEFSEC, AGTR2, and WNT4) were significantly associated with gestational duration. Functional analysis showed that an implicated variant in WNT4 alters the binding of the estrogen receptor. The association between variants in ADCY5 and RAP2C and gestational duration had suggestive significance in the discovery set and significant evidence of association in the replication sets; these variants also showed genomewide significance in a joint analysis. Common variants in EBF1, EEFSEC, and AGTR2 showed association with preterm birth with genomewide significance. An analysis of mother-infant dyads suggested that these variants act at the level of the maternal genome. CONCLUSIONS: In this genomewide association study, we found that variants at the EBF1, EEFSEC, AGTR2, WNT4, ADCY5, and RAP2C loci were associated with gestational duration and variants at the EBF1, EEFSEC, and AGTR2 loci with preterm birth. Previously established roles of these genes in uterine development, maternal nutrition, and vascular control support their mechanistic involvement. (Funded by the March of Dimes and others.).

A multi-stage genome-wide association study of uterine fibroids in African Americans.

Uterine fibroids are benign tumors of the uterus affecting up to 77% of women by menopause. They are the leading indication for hysterectomy, and account for $34 billion annually in the United States. Race/ethnicity and age are the strongest known risk factors. African American (AA) women have higher prevalence, earlier onset, and larger and more numerous fibroids than European American women. We conducted a multi-stage genome-wide association study (GWAS) of fibroid risk among AA women followed by in silico genetically predicted gene expression profiling of top hits. In Stage 1, cases and controls were confirmed by pelvic imaging, genotyped and imputed to 1000 Genomes. Stage 2 used self-reported fibroid and GWAS data from 23andMe, Inc. and the Black Women's Health Study. Associations with fibroid risk were modeled using logistic regression adjusted for principal components, followed by meta-analysis of results. We observed a significant association among 3399 AA cases and 4764 AA controls at rs739187 (risk-allele frequency = 0.27) in CYTH4 (OR (95% confidence interval) = 1.23 (1.16-1.30), p value = 7.82 × 10-9). Evaluation of the genetic association results with MetaXcan identified lower predicted gene expression of CYTH4 in thyroid tissue as significantly associated with fibroid risk (p value = 5.86 × 10-8). In this first multi-stage GWAS for fibroids among AA women, we identified a novel risk locus for fibroids within CYTH4 that impacts gene expression in thyroid and has potential biological relevance for fibroids.

Meta-analysis identifies novel risk loci and yields systematic insights into the biology of male-pattern baldness.

Male-pattern baldness (MPB) is a common and highly heritable trait characterized by androgen-dependent, progressive hair loss from the scalp. Here, we carry out the largest GWAS meta-analysis of MPB to date, comprising 10,846 early-onset cases and 11,672 controls from eight independent cohorts. We identify 63 MPB-associated loci (P<5 × 10-8, METAL) of which 23 have not been reported previously. The 63 loci explain ∼39% of the phenotypic variance in MPB and highlight several plausible candidate genes (FGF5, IRF4, DKK2) and pathways (melatonin signalling, adipogenesis) that are likely to be implicated in the key-pathophysiological features of MPB and may represent promising targets for the development of novel therapeutic options. The data provide molecular evidence that rather than being an isolated trait, MPB shares a substantial biological basis with numerous other human phenotypes and may deserve evaluation as an early prognostic marker, for example, for prostate cancer, sudden cardiac arrest and neurodegenerative disorders.

Two-stage genome-wide association study identifies a novel susceptibility locus associated with melanoma.

Genome-wide association studies have identified 21 susceptibility loci associated with melanoma. These loci implicate genes affecting pigmentation, nevus count, telomere maintenance, and DNA repair in melanoma risk. Here, we report the results of a two-stage genome-wide association study of melanoma. The stage 1 discovery phase consisted of 4,842 self-reported melanoma cases and 286,565 controls of European ancestry from the 23andMe research cohort and the stage 2 replication phase consisted of 1,804 melanoma cases and 1,026 controls from the University of Texas M.D. Anderson Cancer Center. We performed a combined meta-analysis totaling 6,628 melanoma cases and 287,591 controls. Our study replicates 20 of 21 previously known melanoma-loci and confirms the association of the telomerase reverse transcriptase, TERT, with melanoma susceptibility at genome-wide significance. In addition, we uncover a novel polymorphism, rs187843643 (OR = 1.96; 95% CI = [1.54, 2.48]; P = 3.53 x 10-8), associated with melanoma. The SNP rs187842643 lies within a noncoding RNA 177kb downstream of BASP1 (brain associated protein-1). We find that BASP1 expression is suppressed in melanoma as compared with benign nevi, providing additional evidence for a putative role in melanoma pathogenesis.

Gene-based analysis of regulatory variants identifies 4 putative novel asthma risk genes related to nucleotide synthesis and signaling.

BACKGROUND: Hundreds of genetic variants are thought to contribute to variation in asthma risk by modulating gene expression. Methods that increase the power of genome-wide association studies (GWASs) to identify risk-associated variants are needed. OBJECTIVE: We sought to develop a method that aggregates the evidence for association with disease risk across expression quantitative trait loci (eQTLs) of a gene and use this approach to identify asthma risk genes. METHODS: We developed a gene-based test and software package called EUGENE that (1) is applicable to GWAS summary statistics; (2) considers both cis- and trans-eQTLs; (3) incorporates eQTLs identified in different tissues; and (4) uses simulations to account for multiple testing. We applied this approach to 2 published asthma GWASs (combined n = 46,044) and used mouse studies to provide initial functional insights into 2 genes with novel genetic associations. RESULTS: We tested the association between asthma and 17,190 genes that were found to have cis- and/or trans-eQTLs across 16 published eQTL studies. At an empirical FDR of 5%, 48 genes were associated with asthma risk. Of these, for 37, the association was driven by eQTLs located in established risk loci for allergic disease, including 6 genes not previously implicated in disease cause (eg, LIMS1, TINF2, and SAFB). The remaining 11 significant genes represent potential novel genetic associations with asthma. The association with 4 of these replicated in an independent GWAS: B4GALT3, USMG5, P2RY13, and P2RY14, which are genes involved in nucleotide synthesis or nucleotide-dependent cell activation. In mouse studies, P2ry13 and P2ry14-purinergic receptors activated by adenosine 5-diphosphate and UDP-sugars, respectively-were upregulated after allergen challenge, notably in airway epithelial cells, eosinophils, and neutrophils. Intranasal exposure with receptor agonists induced the release of IL-33 and subsequent eosinophil infiltration into the lungs. CONCLUSION: We identified novel associations between asthma and eQTLs for 4 genes related to nucleotide synthesis/signaling and demonstrated the power of gene-based analyses of GWASs.

Genome-wide association study identifies 14 novel risk alleles associated with basal cell carcinoma.

Basal cell carcinoma (BCC) is the most common cancer worldwide with an annual incidence of 2.8 million cases in the United States alone. Previous studies have demonstrated an association between 21 distinct genetic loci and BCC risk. Here, we report the results of a two-stage genome-wide association study of BCC, totalling 17,187 cases and 287,054 controls. We confirm 17 previously reported loci and identify 14 new susceptibility loci reaching genome-wide significance (P<5 × 10(-8), logistic regression). These newly associated SNPs lie within predicted keratinocyte regulatory elements and in expression quantitative trait loci; furthermore, we identify candidate genes and non-coding RNAs involved in telomere maintenance, immune regulation and tumour progression, providing deeper insight into the pathogenesis of BCC.

Genome-wide association study identifies novel susceptibility loci for cutaneous squamous cell carcinoma.

Cutaneous squamous cell carcinoma represents the second most common cutaneous malignancy, affecting 7-11% of Caucasians in the United States. The genetic determinants of susceptibility to cutaneous squamous cell carcinoma remain largely unknown. Here we report the results of a two-stage genome-wide association study of cutaneous squamous cell carcinoma, totalling 7,404 cases and 292,076 controls. Eleven loci reached genome-wide significance (P<5 × 10(-8)) including seven previously confirmed pigmentation-related loci: MC1R, ASIP, TYR, SLC45A2, OCA2, IRF4 and BNC2. We identify an additional four susceptibility loci: 11q23.3 CADM1, a metastasis suppressor gene involved in modifying tumour interaction with cell-mediated immunity; 2p22.3; 7p21.1 AHR, the dioxin receptor involved in anti-apoptotic pathways and melanoma progression; and 9q34.3 SEC16A, a putative oncogene with roles in secretion and cellular proliferation. These susceptibility loci provide deeper insight into the pathogenesis of squamous cell carcinoma.

Germ line variants predispose to both JAK2 V617F clonal hematopoiesis and myeloproliferative neoplasms.

We conducted a genome-wide association study (GWAS) to identify novel predisposition alleles associated with Philadelphia chromosome-negative myeloproliferative neoplasms (MPNs) and JAK2 V617F clonal hematopoiesis in the general population. We recruited a web-based cohort of 726 individuals with polycythemia vera, essential thrombocythemia, and myelofibrosis and 252 637 population controls unselected for hematologic phenotypes. Using a single-nucleotide polymorphism (SNP) array platform with custom probes for the JAK2 V617F mutation (V617F), we identified 497 individuals (0.2%) among the population controls who were V617F carriers. We performed a combined GWAS of the MPN cases plus V617F carriers in the control population (n = 1223) vs the remaining controls who were noncarriers for V617F (n = 252 140). For these MPN cases plus V617F carriers, we replicated the germ line JAK2 46/1 haplotype (rs59384377: odds ratio [OR] = 2.4, P = 6.6 × 10(-89)), previously associated with V617F-positive MPN. We also identified genome-wide significant associations in the TERT gene (rs7705526: OR = 1.8, P = 1.1 × 10(-32)), in SH2B3 (rs7310615: OR = 1.4, P = 3.1 × 10(-14)), and upstream of TET2 (rs1548483: OR = 2.0, P = 2.0 × 10(-9)). These associations were confirmed in a separate replication cohort of 446 V617F carriers vs 169 021 noncarriers. In a joint analysis of the combined GWAS and replication results, we identified additional genome-wide significant predisposition alleles associated with CHEK2, ATM, PINT, and GFI1B All SNP ORs were similar for MPN patients and controls who were V617F carriers. These data indicate that the same germ line variants endow individuals with a predisposition not only to MPN, but also to JAK2 V617F clonal hematopoiesis, a more common phenomenon that may foreshadow the development of an overt neoplasm.