Trans-ethnic and Ancestry-Specific Blood-Cell Genetics in 746,667 Individuals from 5 Global Populations.

Most loci identified by GWASs have been found in populations of European ancestry (EUR). In trans-ethnic meta-analyses for 15 hematological traits in 746,667 participants, including 184,535 non-EUR individuals, we identified 5,552 trait-variant associations at p < 5 × 10-9, including 71 novel associations not found in EUR populations. We also identified 28 additional novel variants in ancestry-specific, non-EUR meta-analyses, including an IL7 missense variant in South Asians associated with lymphocyte count in vivo and IL-7 secretion levels in vitro. Fine-mapping prioritized variants annotated as functional and generated 95% credible sets that were 30% smaller when using the trans-ethnic as opposed to the EUR-only results. We explored the clinical significance and predictive value of trans-ethnic variants in multiple populations and compared genetic architecture and the effect of natural selection on these blood phenotypes between populations. Altogether, our results for hematological traits highlight the value of a more global representation of populations in genetic studies.

The Polygenic and Monogenic Basis of Blood Traits and Diseases.

Blood cells play essential roles in human health, underpinning physiological processes such as immunity, oxygen transport, and clotting, which when perturbed cause a significant global health burden. Here we integrate data from UK Biobank and a large-scale international collaborative effort, including data for 563,085 European ancestry participants, and discover 5,106 new genetic variants independently associated with 29 blood cell phenotypes covering a range of variation impacting hematopoiesis. We holistically characterize the genetic architecture of hematopoiesis, assess the relevance of the omnigenic model to blood cell phenotypes, delineate relevant hematopoietic cell states influenced by regulatory genetic variants and gene networks, identify novel splice-altering variants mediating the associations, and assess the polygenic prediction potential for blood traits and clinical disorders at the interface of complex and Mendelian genetics. These results show the power of large-scale blood cell trait GWAS to interrogate clinically meaningful variants across a wide allelic spectrum of human variation.

Genetic diversity fuels gene discovery for tobacco and alcohol use.

Tobacco and alcohol use are heritable behaviours associated with 15% and 5.3% of worldwide deaths, respectively, due largely to broad increased risk for disease and injury1-4. These substances are used across the globe, yet genome-wide association studies have focused largely on individuals of European ancestries5. Here we leveraged global genetic diversity across 3.4 million individuals from four major clines of global ancestry (approximately 21% non-European) to power the discovery and fine-mapping of genomic loci associated with tobacco and alcohol use, to inform function of these loci via ancestry-aware transcriptome-wide association studies, and to evaluate the genetic architecture and predictive power of polygenic risk within and across populations. We found that increases in sample size and genetic diversity improved locus identification and fine-mapping resolution, and that a large majority of the 3,823 associated variants (from 2,143 loci) showed consistent effect sizes across ancestry dimensions. However, polygenic risk scores developed in one ancestry performed poorly in others, highlighting the continued need to increase sample sizes of diverse ancestries to realize any potential benefit of polygenic prediction.

Ancestry- and sex-specific effects underlying inguinal hernia susceptibility identified in a multiethnic genome-wide association study meta-analysis.

Inguinal hernias are some of the most frequently diagnosed conditions in clinical practice and inguinal hernia repair is the most common procedure performed by general surgeons. Studies of inguinal hernias in non-European populations are lacking, though it is expected that such studies could identify novel loci. Further, the cumulative lifetime incidence of inguinal hernia is nine times greater in men than women, however, it is not clear why this difference exists. We conducted a genome-wide association meta-analysis of inguinal hernia risk across 513 120 individuals (35 774 cases and 477 346 controls) of Hispanic/Latino, African, Asian and European descent, with replication in 728 418 participants (33 491 cases and 694 927 controls) from the 23andMe, Inc dataset. We identified 63 genome-wide significant loci (P < 5 × 10-8), including 41 novel. Ancestry-specific analyses identified two loci (LYPLAL1-AS1/SLC30A10 and STXBP6-NOVA1) in African ancestry individuals. Sex-stratified analyses identified two loci (MYO1D and ZBTB7C) that are specific to women, and four (EBF2, EMX2/RAB11FIP2, VCL and FAM9A/FAM9B) that are specific to men. Functional experiments demonstrated that several of the associated regions (EFEMP1 and LYPLAL1-SLC30A10) function as enhancers and show differential activity between risk and reference alleles. Our study highlights the importance of large-scale genomic studies in ancestrally diverse populations for identifying ancestry-specific inguinal hernia susceptibility loci and provides novel biological insights into inguinal hernia etiology.

Whole genome sequence analysis of platelet traits in the NHLBI Trans-Omics for Precision Medicine (TOPMed) initiative.

Platelets play a key role in thrombosis and hemostasis. Platelet count (PLT) and mean platelet volume (MPV) are highly heritable quantitative traits, with hundreds of genetic signals previously identified, mostly in European ancestry populations. We here utilize whole genome sequencing (WGS) from NHLBI's Trans-Omics for Precision Medicine initiative (TOPMed) in a large multi-ethnic sample to further explore common and rare variation contributing to PLT (n = 61 200) and MPV (n = 23 485). We identified and replicated secondary signals at MPL (rs532784633) and PECAM1 (rs73345162), both more common in African ancestry populations. We also observed rare variation in Mendelian platelet-related disorder genes influencing variation in platelet traits in TOPMed cohorts (not enriched for blood disorders). For example, association of GP9 with lower PLT and higher MPV was partly driven by a pathogenic Bernard-Soulier syndrome variant (rs5030764, p.Asn61Ser), and the signals at TUBB1 and CD36 were partly driven by loss of function variants not annotated as pathogenic in ClinVar (rs199948010 and rs571975065). However, residual signal remained for these gene-based signals after adjusting for lead variants, suggesting that additional variants in Mendelian genes with impacts in general population cohorts remain to be identified. Gene-based signals were also identified at several genome-wide association study identified loci for genes not annotated for Mendelian platelet disorders (PTPRH, TET2, CHEK2), with somatic variation driving the result at TET2. These results highlight the value of WGS in populations of diverse genetic ancestry to identify novel regulatory and coding signals, even for well-studied traits like platelet traits.

Whole-genome sequencing association analysis of quantitative red blood cell phenotypes: The NHLBI TOPMed program.

Whole-genome sequencing (WGS), a powerful tool for detecting novel coding and non-coding disease-causing variants, has largely been applied to clinical diagnosis of inherited disorders. Here we leveraged WGS data in up to 62,653 ethnically diverse participants from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program and assessed statistical association of variants with seven red blood cell (RBC) quantitative traits. We discovered 14 single variant-RBC trait associations at 12 genomic loci, which have not been reported previously. Several of the RBC trait-variant associations (RPN1, ELL2, MIDN, HBB, HBA1, PIEZO1, and G6PD) were replicated in independent GWAS datasets imputed to the TOPMed reference panel. Most of these discovered variants are rare/low frequency, and several are observed disproportionately among non-European Ancestry (African, Hispanic/Latino, or East Asian) populations. We identified a 3 bp indel p.Lys2169del (g.88717175_88717177TCT[4]) (common only in the Ashkenazi Jewish population) of PIEZO1, a gene responsible for the Mendelian red cell disorder hereditary xerocytosis (MIM: 194380), associated with higher mean corpuscular hemoglobin concentration (MCHC). In stepwise conditional analysis and in gene-based rare variant aggregated association analysis, we identified several of the variants in HBB, HBA1, TMPRSS6, and G6PD that represent the carrier state for known coding, promoter, or splice site loss-of-function variants that cause inherited RBC disorders. Finally, we applied base and nuclease editing to demonstrate that the sentinel variant rs112097551 (nearest gene RPN1) acts through a cis-regulatory element that exerts long-range control of the gene RUVBL1 which is essential for hematopoiesis. Together, these results demonstrate the utility of WGS in ethnically diverse population-based samples and gene editing for expanding knowledge of the genetic architecture of quantitative hematologic traits and suggest a continuum between complex trait and Mendelian red cell disorders.

Whole-genome sequencing in diverse subjects identifies genetic correlates of leukocyte traits: The NHLBI TOPMed program.

Many common and rare variants associated with hematologic traits have been discovered through imputation on large-scale reference panels. However, the majority of genome-wide association studies (GWASs) have been conducted in Europeans, and determining causal variants has proved challenging. We performed a GWAS of total leukocyte, neutrophil, lymphocyte, monocyte, eosinophil, and basophil counts generated from 109,563,748 variants in the autosomes and the X chromosome in the Trans-Omics for Precision Medicine (TOPMed) program, which included data from 61,802 individuals of diverse ancestry. We discovered and replicated 7 leukocyte trait associations, including (1) the association between a chromosome X, pseudo-autosomal region (PAR), noncoding variant located between cytokine receptor genes (CSF2RA and CLRF2) and lower eosinophil count; and (2) associations between single variants found predominantly among African Americans at the S1PR3 (9q22.1) and HBB (11p15.4) loci and monocyte and lymphocyte counts, respectively. We further provide evidence indicating that the newly discovered eosinophil-lowering chromosome X PAR variant might be associated with reduced susceptibility to common allergic diseases such as atopic dermatitis and asthma. Additionally, we found a burden of very rare FLT3 (13q12.2) variants associated with monocyte counts. Together, these results emphasize the utility of whole-genome sequencing in diverse samples in identifying associations missed by European-ancestry-driven GWASs.

Genetic variation affects morphological retinal phenotypes extracted from UK Biobank optical coherence tomography images.

Optical Coherence Tomography (OCT) enables non-invasive imaging of the retina and is used to diagnose and manage ophthalmic diseases including glaucoma. We present the first large-scale genome-wide association study of inner retinal morphology using phenotypes derived from OCT images of 31,434 UK Biobank participants. We identify 46 loci associated with thickness of the retinal nerve fibre layer or ganglion cell inner plexiform layer. Only one of these loci has been associated with glaucoma, and despite its clear role as a biomarker for the disease, Mendelian randomisation does not support inner retinal thickness being on the same genetic causal pathway as glaucoma. We extracted overall retinal thickness at the fovea, representative of foveal hypoplasia, with which three of the 46 SNPs were associated. We additionally associate these three loci with visual acuity. In contrast to the Mendelian causes of severe foveal hypoplasia, our results suggest a spectrum of foveal hypoplasia, in part genetically determined, with consequences on visual function.

Correction: Genetic variation affects morphological retinal phenotypes extracted from UK Biobank optical coherence tomography images.

[This corrects the article DOI: 10.1371/journal.pgen.1009497.].

A large-scale transcriptome-wide association study (TWAS) of 10 blood cell phenotypes reveals complexities of TWAS fine-mapping.

Hematological measures are important intermediate clinical phenotypes for many acute and chronic diseases and are highly heritable. Although genome-wide association studies (GWAS) have identified thousands of loci containing trait-associated variants, the causal genes underlying these associations are often uncertain. To better understand the underlying genetic regulatory mechanisms, we performed a transcriptome-wide association study (TWAS) to systematically investigate the association between genetically predicted gene expression and hematological measures in 54,542 Europeans from the Genetic Epidemiology Research on Aging cohort. We found 239 significant gene-trait associations with hematological measures; we replicated 71 associations at p < 0.05 in a TWAS meta-analysis consisting of up to 35,900 Europeans from the Women's Health Initiative, Atherosclerosis Risk in Communities Study, and BioMe Biobank. Additionally, we attempted to refine this list of candidate genes by performing conditional analyses, adjusting for individual variants previously associated with hematological measures, and performed further fine-mapping of TWAS loci. To facilitate interpretation of our findings, we designed an R Shiny application to interactively visualize our TWAS results by integrating them with additional genetic data sources (GWAS, TWAS from multiple reference panels, conditional analyses, known GWAS variants, etc.). Our results and application highlight frequently overlooked TWAS challenges and illustrate the complexity of TWAS fine-mapping.

Genome-Wide Association Study Identifies Two Common Loci Associated with Pigment Dispersion Syndrome/Pigmentary Glaucoma and Implicates Myopia in its Development.

PURPOSE: To identify genetic variants associated with pigment dispersion syndrome (PDS) and pigmentary glaucoma (PG) in unrelated patients and to further understand the genetic and potentially causal relationships between PDS and associated risk factors. DESIGN: A 2-stage genome-wide association meta-analysis with replication and subsequent in silico analyses including Mendelian randomization. PARTICIPANTS: A total of 574 cases with PG or PDS and 52 627 controls of European descent. METHODS: Genome-wide association analyses were performed in 4 cohorts and meta-analyzed in 3 stages: (1) a discovery meta-analysis was performed in 3 cohorts, (2) replication was performed in the fourth cohort, and (3) all 4 cohorts were meta-analyzed to increase statistical power. Two-sample Mendelian randomization was used to determine whether refractive error and intraocular pressure exert causal effects over PDS. MAIN OUTCOME MEASURES: The association of genetic variants with PDS and whether myopia exerts causal effects over PDS. RESULTS: Significant association was present at 2 novel loci for PDS/PG. These loci and follow-up analyses implicate the genes gamma secretase activator protein (GSAP) (lead single nucleotide polymorphism [SNP]: rs9641220, P = 6.0×10-10) and glutamate metabotropic receptor 5 (GRM5)/TYR (lead SNP: rs661177, P = 3.9×10-9) as important factors in disease risk. Mendelian randomization showed significant evidence that negative refractive error (myopia) exerts a direct causal effect over PDS (P = 8.86×10-7). CONCLUSIONS: Common SNPs relating to the GSAP and GRM5/TYR genes are associated risk factors for the development of PDS and PG. Although myopia is a known risk factor, this study uses genetic data to demonstrate that myopia is, in part, a cause of PDS and PG.

Use of >100,000 NHLBI Trans-Omics for Precision Medicine (TOPMed) Consortium whole genome sequences improves imputation quality and detection of rare variant associations in admixed African and Hispanic/Latino populations.

Most genome-wide association and fine-mapping studies to date have been conducted in individuals of European descent, and genetic studies of populations of Hispanic/Latino and African ancestry are limited. In addition, these populations have more complex linkage disequilibrium structure. In order to better define the genetic architecture of these understudied populations, we leveraged >100,000 phased sequences available from deep-coverage whole genome sequencing through the multi-ethnic NHLBI Trans-Omics for Precision Medicine (TOPMed) program to impute genotypes into admixed African and Hispanic/Latino samples with genome-wide genotyping array data. We demonstrated that using TOPMed sequencing data as the imputation reference panel improves genotype imputation quality in these populations, which subsequently enhanced gene-mapping power for complex traits. For rare variants with minor allele frequency (MAF) < 0.5%, we observed a 2.3- to 6.1-fold increase in the number of well-imputed variants, with 11-34% improvement in average imputation quality, compared to the state-of-the-art 1000 Genomes Project Phase 3 and Haplotype Reference Consortium reference panels. Impressively, even for extremely rare variants with minor allele count <10 (including singletons) in the imputation target samples, average information content rescued was >86%. Subsequent association analyses of TOPMed reference panel-imputed genotype data with hematological traits (hemoglobin (HGB), hematocrit (HCT), and white blood cell count (WBC)) in ~21,600 African-ancestry and ~21,700 Hispanic/Latino individuals identified associations with two rare variants in the HBB gene (rs33930165 with higher WBC [p = 8.8x10-15] in African populations, rs11549407 with lower HGB [p = 1.5x10-12] and HCT [p = 8.8x10-10] in Hispanics/Latinos). By comparison, neither variant would have been genome-wide significant if either 1000 Genomes Project Phase 3 or Haplotype Reference Consortium reference panels had been used for imputation. Our findings highlight the utility of the TOPMed imputation reference panel for identification of novel rare variant associations not previously detected in similarly sized genome-wide studies of under-represented African and Hispanic/Latino populations.

Dysfunction of lipid sensor GPR120 leads to obesity in both mouse and human.

Free fatty acids provide an important energy source as nutrients, and act as signalling molecules in various cellular processes. Several G-protein-coupled receptors have been identified as free-fatty-acid receptors important in physiology as well as in several diseases. GPR120 (also known as O3FAR1) functions as a receptor for unsaturated long-chain free fatty acids and has a critical role in various physiological homeostasis mechanisms such as adipogenesis, regulation of appetite and food preference. Here we show that GPR120-deficient mice fed a high-fat diet develop obesity, glucose intolerance and fatty liver with decreased adipocyte differentiation and lipogenesis and enhanced hepatic lipogenesis. Insulin resistance in such mice is associated with reduced insulin signalling and enhanced inflammation in adipose tissue. In human, we show that GPR120 expression in adipose tissue is significantly higher in obese individuals than in lean controls. GPR120 exon sequencing in obese subjects reveals a deleterious non-synonymous mutation (p.R270H) that inhibits GPR120 signalling activity. Furthermore, the p.R270H variant increases the risk of obesity in European populations. Overall, this study demonstrates that the lipid sensor GPR120 has a key role in sensing dietary fat and, therefore, in the control of energy balance in both humans and rodents.

Association of cardiovascular risk factors with disease severity in cerebral cavernous malformation type 1 subjects with the common Hispanic mutation.

BACKGROUND: Cerebral cavernous malformations (CCM) are enlarged vascular lesions affecting 0.1-0.5% of the population worldwide and causing hemorrhagic strokes, seizures, and neurological deficits. Familial CCM type 1 (CCM1) is an autosomal dominant disease caused by mutations in the Krev Interaction Trapped 1 (KRIT1/CCM1) gene, and is characterized by multiple brain lesions whose number and size increase with age. The number of lesions varies widely for unknown reasons, even among carriers of similar ages with the same mutation. The purpose of this study was to investigate whether cardiovascular (CV) risk factors influence potential markers of familial CCM1 disease severity, such as lesion count and history of intracerebral hemorrhage. METHODS: We analyzed baseline data from 185 Hispanic subjects, enrolled in the Brain Vascular Malformation Consortium study between June 2010 and March 2013. All subjects were carriers of the founder Q455X 'Common Hispanic Mutation' (CHM) in the KRIT1 gene, and had a clinical diagnosis of CCM or had an affected first- or second-degree relative with CCM. We performed a cross-sectional study, collecting detailed clinical information of CCM1-CHM subjects and cerebral susceptibility-weighted magnetic resonance imaging to assess lesion count. Linear or logistic regression analysis of log-lesion count or history of intracerebral hemorrhage and CV risk factors (age, gender, obesity, diabetes, hypertension, hyperlipidemia and smoking status) and related quantitative traits (body mass index, glycosylated hemoglobin levels, blood pressure, lipids levels and pack-years of cigarette smoking) was performed accommodating familial clustering. RESULTS: CCM1-CHM subjects were mainly female (63.8%) and symptomatic at presentation (63.2%). Lesion count was highly variable (mean ± SD: 57.7 ± 110.6; range: 0-713); 90% of CCM1-CHM subjects had multiple lesions at enrollment. Age (p < 0.001) was positively correlated with lesion count and male gender (p = 0.035) was associated with a greater number of lesions. Obesity (p = 0.001) and higher body mass index (p = 0.002) were associated with fewer lesions. No association with hypertension was detected, however, systolic blood pressure (p = 0.002) was associated with fewer lesions. No significant association with lesion count was observed for diabetes, hyperlipidemia, smoking status or for related quantitative traits. History of intracerebral hemorrhage was not significantly associated with any CV risk factors, however, we found borderline associations of hemorrhage with obesity (p = 0.062), systolic blood pressure (p = 0.083) and pack-years of cigarette smoking (p = 0.055). After correction for multiple testing, age and obesity remained significantly associated with lesion count in CCM1-CHM subjects. CONCLUSIONS: These results suggest that several CV risk factors explain some of the variability in lesion count in Hispanic CCM1-CHM subjects. Although age, gender, obesity, body mass index and systolic blood pressure may influence familial CCM1 disease severity, further longitudinal studies in larger sample sizes are essential to confirm these findings.

Polymorphisms in inflammatory and immune response genes associated with cerebral cavernous malformation type 1 severity.

BACKGROUND: Familial cerebral cavernous malformation type 1 (CCM1) is an autosomal dominant disease caused by mutations in the Krev Interaction Trapped 1 (KRIT1/CCM1) gene, and characterized by multiple brain lesions that often result in intracerebral hemorrhage (ICH), seizures, and neurological deficits. Carriers of the same genetic mutation can present with variable symptoms and severity of disease, suggesting the influence of modifier factors. Evidence is emerging that inflammation and immune response play a role in the pathogenesis of CCM. The purpose of this study was to investigate whether common variants in inflammatory and immune response genes influence the severity of familial CCM1 disease, as manifested by ICH and greater brain lesion count. METHODS: Hispanic CCM1 patients (n=188) harboring the founder Q455X 'common Hispanic mutation' (CHM) in the KRIT1 gene were analyzed at baseline. Participants were enrolled between June 2010 and March 2014 either through the Brain Vascular Malformation Consortium (BVMC) study or through the Angioma Alliance organization. Clinical assessment and cerebral susceptibility-weighted magnetic resonance imaging were performed to determine ICH as well as total and large (≥5 mm in diameter) lesion counts. Samples were genotyped on the Affymetrix Axiom Genome-Wide LAT1 Human Array. We analyzed 830 variants in 56 inflammatory and immune response genes for association with ICH and residuals of log-transformed total or large lesion count adjusted for age at enrollment and gender. Variants were analyzed individually or grouped by sub-pathways or whole pathways. RESULTS: At baseline, 30.3% of CCM1-CHM subjects had ICH, with a mean ± standard deviation (SD) of 60.1±115.0 (range 0-713) for total lesions and 4.9±8.7 (range 0-104) for large lesions. The heritability estimates explained by all autosomal variants were 0.20 (SE=0.31), 0.81 (SE=0.17), and 0.48 (SE=0.19), for ICH, total lesion count, and large lesion count, respectively. TGFBR2 rs9823731 was significantly associated with ICH as well as with the total and large lesion counts (p≤0.017). Further, IL-4 rs9327638, CD14 rs778588, IL-6R rs114660934 and MSR1 rs62489577 were associated with two markers of disease severity. Finally, the whole pathway was associated with total lesion count (p=0.005) with TLR-4 rs10759930, CD14 rs778588, IL-6R rs114660934 and IGH rs57767447 mainly bearing this association. Eicosanoid signaling, extracellular pattern recognition, and immune response sub-pathways were also associated with the total lesion count. CONCLUSIONS: These results suggest that polymorphisms in inflammatory and immune response pathways contribute to variability in CCM1 disease severity and might be used as predictors of disease severity. In particular, TGFBR2 rs9823731 was associated with all three markers of CCM1 disease severity tested, suggesting that TGFBR2 might be a key participant in the mechanism underlying CCM1 disease severity and phenotype variability. However, further longitudinal studies in larger sample sizes are needed to confirm these findings.

Socio-Demographic and Preoperative Clinical Factors Associated With 5-Year Weight Trajectories After Roux-en-Y Gastric Bypass and Sleeve Gastrectomy.

Objective: To determine whether socio-demographic and preoperative clinical factors contribute to the percent total body weight loss (%TBWL) after bariatric surgery (BS). Background: BS is the most effective long-term treatment for medically complicated obesity. More information is needed about the factors that contribute to postoperative %TBWL in large and ethnically diverse cohorts. Methods: This retrospective study conducted in the Kaiser Permanente Northern California region included 7698 patients who underwent Roux-en-Y gastric bypass (RYGB) or sleeve gastrectomy (SG) between January 2009 and March 2015. Trajectory analyses were conducted from 5-year follow-up data to assign patients to "low," "average," or "high" postoperative %TBWL groups. We then evaluated whether age, sex, race/ethnicity, neighborhood deprivation index and preoperative body mass index (BMI)/weight loss, diabetes, hypertension, and sleep apnea contributed to postoperative %TBWL using logistic regression models. Results: Of 7698 patients (83.2% women), 48.6% underwent a RYGB and 51.4% underwent a SG. Postoperative %TBWL trajectories over 5 years were obtained in 6229 (81%) of 7698 eligible patients. About 27.8% and 29.3% of patients followed the "low" postoperative %TBWL trajectory, for RYGB and SG, respectively. Men, older patients, and Asian, Black, and Hispanic/Latino patients were more likely to be classified in the low postoperative %TBWL group. Patients showing lower postoperative %TBWL had a lower preoperative BMI (but lost less weight before surgery) and were more likely to have preoperative comorbidities. Conclusions: This study confirms and extends prior findings of the effects of several demographic and preoperative clinical factors on postoperative weight loss. Findings could improve the support of patients to achieve desired surgical outcomes.

Association between lifetime smoking and cutaneous squamous cell carcinoma: A 2-sample Mendelian randomization study.

Background/Purpose: Cutaneous squamous cell carcinoma (cSCC) is one of the most common malignancies worldwide. While several environmental risk factors for cSCC are well established, there is conflicting evidence on cigarette smoking (and its potential causal effect) and cSCC risk. Furthermore, it is unclear if these potential associations represent causal, modifiable risk factors for cSCC development. This study aims to assess the nature of the associations between cigarette smoking traits (smoking initiation, amount smoked, and lifetime smoking exposure) and cSCC risk using two-sample Mendelian randomization analyses. Methods: Genetic instruments, based on common genetic variants associated with cigarette smoking traits (P < 5 × 10-8), were derived from published genome-wide association studies (GWASs). For cSCC, we used GWAS summary statistics from the Kaiser Permanente GERA cohort (7701 cSCC cases and 60,167 controls; all non-Hispanic Whites). Results: We found modest evidence that genetically determined lifetime smoking was associated with cSCC (inverse-variance weighted method: OR[95% CI] = 1.47[1.09-1.98]; P = .012), suggesting it may be a causal risk factor for cSCC. We did not detect any evidence of association between genetically determined smoking initiation or amount smoked and cSCC risk. Conclusion: Study findings highlight the importance of smoking prevention and may support risk-stratified cSCC screening strategies based on carcinogen exposure and other genetic and clinical information.

Multi-ancestry genome-wide meta-analysis identifies novel basal cell carcinoma loci and shared genetic effects with squamous cell carcinoma.

Basal cell carcinoma (BCC) is one of the most common malignancies worldwide, yet its genetic determinants are incompletely defined. We perform a European ancestry genome-wide association (GWA) meta-analysis and a Hispanic/Latino ancestry GWA meta-analysis and meta-analyze both in a multi-ancestry GWAS meta-analysis of BCC, totaling 50,531 BCC cases and 762,234 controls from four cohorts (GERA, Mass-General Brigham Biobank, UK Biobank, and 23andMe research cohort). Here we identify 122 BCC-associated loci, of which 36 were novel, and subsequently fine-mapped these associations. We also identify an association of the well-known pigment gene SLC45A2 as well as associations at RCC2 and CLPTM1L with BCC in Hispanic/Latinos. We examine these BCC loci for association with cutaneous squamous cell carcinoma (cSCC) in 16,407 SCC cases and 762,486 controls of European ancestry, and 33 SNPs show evidence of association. Our study findings provide important insights into the genetic basis of BCC and cSCC susceptibility.

Multi-tissue transcriptome-wide association study identifies novel candidate susceptibility genes for cataract.

Introduction: Cataract is the leading cause of blindness among the elderly worldwide. Twin and family studies support an important role for genetic factors in cataract susceptibility with heritability estimates up to 58%. To date, 55 loci for cataract have been identified by genome-wide association studies (GWAS), however, much work remains to identify the causal genes. Here, we conducted a transcriptome-wide association study (TWAS) of cataract to prioritize causal genes and identify novel ones, and examine the impact of their expression. Methods: We performed tissue-specific and multi-tissue TWAS analyses to assess associations between imputed gene expression from 54 tissues (including 49 from the Genotype Tissue Expression (GTEx) Project v8) with cataract using FUSION software. Meta-analyzed GWAS summary statistics from 59,944 cataract cases and 478,571 controls, all of European ancestry and from two cohorts (GERA and UK Biobank) were used. We then examined the expression of the novel genes in the lens tissue using the iSyTE database. Results: Across tissue-specific and multi-tissue analyses, we identified 99 genes for which genetically predicted gene expression was associated with cataract after correcting for multiple testing. Of these 99 genes, 20 (AC007773.1, ANKH, ASIP, ATP13A2, CAPZB, CEP95, COQ6, CREB1, CROCC, DDX5, EFEMP1, EIF2S2, ESRRB, GOSR2, HERC4, INSRR, NIPSNAP2, PICALM, SENP3, and SH3YL1) did not overlap with previously reported cataract-associated loci. Tissue-specific analysis identified 202 significant gene-tissue associations for cataract, of which 166 (82.2%), representing 9 unique genes, were attributed to the previously reported 11q13.3 locus. Tissue-enrichment analysis revealed that gastrointestinal tissues represented one of the highest proportions of the Bonferroni-significant gene-tissue associations (21.3%). Moreover, this gastrointestinal tissue type was the only anatomical category significantly enriched in our results, after correcting for the number of tissue donors and imputable genes for each reference panel. Finally, most of the novel cataract genes (e.g., Capzb) were robustly expressed in iSyTE lens data. Discussion: Our results provide evidence of the utility of imputation-based TWAS approaches to characterize known GWAS risk loci and identify novel candidate genes that may increase our understanding of cataract etiology. Our findings also highlight the fact that expression of genes associated with cataract susceptibility is not necessarily restricted to lens tissue.