A multi-cohort genome-wide association study in African ancestry individuals reveals risk loci for primary open-angle glaucoma.

Primary open-angle glaucoma (POAG), the leading cause of irreversible blindness worldwide, disproportionately affects individuals of African ancestry. We conducted a genome-wide association study (GWAS) for POAG in 11,275 individuals of African ancestry (6,003 cases; 5,272 controls). We detected 46 risk loci associated with POAG at genome-wide significance. Replication and post-GWAS analyses, including functionally informed fine-mapping, multiple trait co-localization, and in silico validation, implicated two previously undescribed variants (rs1666698 mapping to DBF4P2; rs34957764 mapping to ROCK1P1) and one previously associated variant (rs11824032 mapping to ARHGEF12) as likely causal. For individuals of African ancestry, a polygenic risk score (PRS) for POAG from our mega-analysis (African ancestry individuals) outperformed a PRS from summary statistics of a much larger GWAS derived from European ancestry individuals. This study quantifies the genetic architecture similarities and differences between African and non-African ancestry populations for this blinding disease.

Complex effects of sequence variants on lipid levels and coronary artery disease.

Many sequence variants have additive effects on blood lipid levels and, through that, on the risk of coronary artery disease (CAD). We show that variants also have non-additive effects and interact to affect lipid levels as well as affecting variance and correlations. Variance and correlation effects are often signatures of epistasis or gene-environmental interactions. These complex effects can translate into CAD risk. For example, Trp154Ter in FUT2 protects against CAD among subjects with the A1 blood group, whereas it associates with greater risk of CAD in others. His48Arg in ADH1B interacts with alcohol consumption to affect lipid levels and CAD. The effect of variants in TM6SF2 on blood lipids is greatest among those who never eat oily fish but absent from those who often do. This work demonstrates that variants that affect variance of quantitative traits can allow for the discovery of epistasis and interactions of variants with the environment.

Limb development genes underlie variation in human fingerprint patterns.

Fingerprints are of long-standing practical and cultural interest, but little is known about the mechanisms that underlie their variation. Using genome-wide scans in Han Chinese cohorts, we identified 18 loci associated with fingerprint type across the digits, including a genetic basis for the long-recognized "pattern-block" correlations among the middle three digits. In particular, we identified a variant near EVI1 that alters regulatory activity and established a role for EVI1 in dermatoglyph patterning in mice. Dynamic EVI1 expression during human development supports its role in shaping the limbs and digits, rather than influencing skin patterning directly. Trans-ethnic meta-analysis identified 43 fingerprint-associated loci, with nearby genes being strongly enriched for general limb development pathways. We also found that fingerprint patterns were genetically correlated with hand proportions. Taken together, these findings support the key role of limb development genes in influencing the outcome of fingerprint patterning.

Root System Depth in Arabidopsis Is Shaped by EXOCYST70A3 via the Dynamic Modulation of Auxin Transport.

Root system architecture (RSA), the distribution of roots in soil, plays a major role in plant survival. RSA is shaped by multiple developmental processes that are largely governed by the phytohormone auxin, suggesting that auxin regulates responses of roots that are important for local adaptation. However, auxin has a central role in numerous processes, and it is unclear which molecular mechanisms contribute to the variation in RSA for environmental adaptation. Using natural variation in Arabidopsis, we identify EXOCYST70A3 as a modulator of the auxin system that causes variation in RSA by acting on PIN4 protein distribution. Allelic variation and genetic perturbation of EXOCYST70A3 lead to alteration of root gravitropic responses, resulting in a different RSA depth profile and drought resistance. Overall our findings suggest that the local modulation of the pleiotropic auxin pathway can gives rise to distinct RSAs that can be adaptive in specific environments.

Genomic Analyses from Non-invasive Prenatal Testing Reveal Genetic Associations, Patterns of Viral Infections, and Chinese Population History.

We analyze whole-genome sequencing data from 141,431 Chinese women generated for non-invasive prenatal testing (NIPT). We use these data to characterize the population genetic structure and to investigate genetic associations with maternal and infectious traits. We show that the present day distribution of alleles is a function of both ancient migration and very recent population movements. We reveal novel phenotype-genotype associations, including several replicated associations with height and BMI, an association between maternal age and EMB, and between twin pregnancy and NRG1. Finally, we identify a unique pattern of circulating viral DNA in plasma with high prevalence of hepatitis B and other clinically relevant maternal infections. A GWAS for viral infections identifies an exceptionally strong association between integrated herpesvirus 6 and MOV10L1, which affects piwi-interacting RNA (piRNA) processing and PIWI protein function. These findings demonstrate the great value and potential of accumulating NIPT data for worldwide medical and genetic analyses.

A Genetic Variant Associated with Five Vascular Diseases Is a Distal Regulator of Endothelin-1 Gene Expression.

Genome-wide association studies (GWASs) implicate the PHACTR1 locus (6p24) in risk for five vascular diseases, including coronary artery disease, migraine headache, cervical artery dissection, fibromuscular dysplasia, and hypertension. Through genetic fine mapping, we prioritized rs9349379, a common SNP in the third intron of the PHACTR1 gene, as the putative causal variant. Epigenomic data from human tissue revealed an enhancer signature at rs9349379 exclusively in aorta, suggesting a regulatory function for this SNP in the vasculature. CRISPR-edited stem cell-derived endothelial cells demonstrate rs9349379 regulates expression of endothelin 1 (EDN1), a gene located 600 kb upstream of PHACTR1. The known physiologic effects of EDN1 on the vasculature may explain the pattern of risk for the five associated diseases. Overall, these data illustrate the integration of genetic, phenotypic, and epigenetic analysis to identify the biologic mechanism by which a common, non-coding variant can distally regulate a gene and contribute to the pathogenesis of multiple vascular diseases.

A Natural Allele of a Transcription Factor in Rice Confers Broad-Spectrum Blast Resistance.

Rice feeds half the world's population, and rice blast is often a destructive disease that results in significant crop loss. Non-race-specific resistance has been more effective in controlling crop diseases than race-specific resistance because of its broad spectrum and durability. Through a genome-wide association study, we report the identification of a natural allele of a C2H2-type transcription factor in rice that confers non-race-specific resistance to blast. A survey of 3,000 sequenced rice genomes reveals that this allele exists in 10% of rice, suggesting that this favorable trait has been selected through breeding. This allele causes a single nucleotide change in the promoter of the bsr-d1 gene, which results in reduced expression of the gene through the binding of the repressive MYB transcription factor and, consequently, an inhibition of H2O2 degradation and enhanced disease resistance. Our discovery highlights this novel allele as a strategy for breeding durable resistance in rice.

Examining the shared etiology of psychopathology with genome-wide association studies.

Genome-wide association studies (GWASs) have ushered in a new era of reproducible discovery in psychiatric genetics. The field has now identified hundreds of common genetic variants that are associated with mental disorders, and many of them influence more than one disorder. By advancing the understanding of causal biology underlying psychopathology, GWAS results are poised to inform the development of novel therapeutics, stratification of at-risk patients, and perhaps even the revision of top-down classification systems in psychiatry. Here, we provide a concise review of GWAS findings with an emphasis on findings that have elucidated the shared genetic etiology of psychopathology, summarizing insights at three levels of analysis: 1) genome-wide architecture; 2) networks, pathways, and gene sets; and 3) individual variants/genes. Three themes emerge from these efforts. First, all psychiatric phenotypes are heritable, highly polygenic, and influenced by many pleiotropic variants with incomplete penetrance. Second, GWAS results highlight the broad etiological roles of neuronal biology, system-wide effects over localized effects, and early neurodevelopment as a critical period. Third, many loci that are robustly associated with multiple forms of psychopathology harbor genes that are involved in synaptic structure and function. Finally, we conclude our review by discussing the implications that GWAS results hold for the field of psychiatry, as well as expected challenges and future directions in the next stage of psychiatric genetics.

Common variants increase risk for congenital diaphragmatic hernia within the context of de novo variants.

Congenital diaphragmatic hernia (CDH) is a severe congenital anomaly often accompanied by other structural anomalies and/or neurobehavioral manifestations. Rare de novo protein-coding variants and copy-number variations contribute to CDH in the population. However, most individuals with CDH remain genetically undiagnosed. Here, we perform integrated de novo and common-variant analyses using 1,469 CDH individuals, including 1,064 child-parent trios and 6,133 ancestry-matched, unaffected controls for the genome-wide association study. We identify candidate CDH variants in 15 genes, including eight novel genes, through deleterious de novo variants. We further identify two genomic loci contributing to CDH risk through common variants with similar effect sizes among Europeans and Latinx. Both loci are in putative transcriptional regulatory regions of developmental patterning genes. Estimated heritability in common variants is ∼19%. Strikingly, there is no significant difference in estimated polygenic risk scores between isolated and complex CDH or between individuals harboring deleterious de novo variants and individuals without these variants. The data support a polygenic model as part of the CDH genetic architecture.

An integrative framework to prioritize genes in more than 500 loci associated with body mass index.

Obesity is a major risk factor for a myriad of diseases, affecting >600 million people worldwide. Genome-wide association studies (GWASs) have identified hundreds of genetic variants that influence body mass index (BMI), a commonly used metric to assess obesity risk. Most variants are non-coding and likely act through regulating genes nearby. Here, we apply multiple computational methods to prioritize the likely causal gene(s) within each of the 536 previously reported GWAS-identified BMI-associated loci. We performed summary-data-based Mendelian randomization (SMR), FINEMAP, DEPICT, MAGMA, transcriptome-wide association studies (TWASs), mutation significance cutoff (MSC), polygenic priority score (PoPS), and the nearest gene strategy. Results of each method were weighted based on their success in identifying genes known to be implicated in obesity, ranking all prioritized genes according to a confidence score (minimum: 0; max: 28). We identified 292 high-scoring genes (≥11) in 264 loci, including genes known to play a role in body weight regulation (e.g., DGKI, ANKRD26, MC4R, LEPR, BDNF, GIPR, AKT3, KAT8, MTOR) and genes related to comorbidities (e.g., FGFR1, ISL1, TFAP2B, PARK2, TCF7L2, GSK3B). For most of the high-scoring genes, however, we found limited or no evidence for a role in obesity, including the top-scoring gene BPTF. Many of the top-scoring genes seem to act through a neuronal regulation of body weight, whereas others affect peripheral pathways, including circadian rhythm, insulin secretion, and glucose and carbohydrate homeostasis. The characterization of these likely causal genes can increase our understanding of the underlying biology and offer avenues to develop therapeutics for weight loss.

Genetic determinants of IgG antibody response to COVID-19 vaccination.

Human humoral immune responses to SARS-CoV-2 vaccines exhibit substantial inter-individual variability and have been linked to vaccine efficacy. To elucidate the underlying mechanism behind this variability, we conducted a genome-wide association study (GWAS) on the anti-spike IgG serostatus of UK Biobank participants who were previously uninfected by SARS-CoV-2 and had received either the first dose (n = 54,066) or the second dose (n = 46,232) of COVID-19 vaccines. Our analysis revealed significant genome-wide associations between the IgG antibody serostatus following the initial vaccine and human leukocyte antigen (HLA) class II alleles. Specifically, the HLA-DRB1∗13:02 allele (MAF = 4.0%, OR = 0.75, p = 2.34e-16) demonstrated the most statistically significant protective effect against IgG seronegativity. This protective effect was driven by an alteration from arginine (Arg) to glutamic acid (Glu) at position 71 on HLA-DRß1 (p = 1.88e-25), leading to a change in the electrostatic potential of pocket 4 of the peptide binding groove. Notably, the impact of HLA alleles on IgG responses was cell type specific, and we observed a shared genetic predisposition between IgG status and susceptibility/severity of COVID-19. These results were replicated within independent cohorts where IgG serostatus was assayed by two different antibody serology tests. Our findings provide insights into the biological mechanism underlying individual variation in responses to COVID-19 vaccines and highlight the need to consider the influence of constitutive genetics when designing vaccination strategies for optimizing protection and control of infectious disease across diverse populations.

Genetic Factors in Mammalian Prion Diseases.

Mammalian prion diseases are a group of neurodegenerative conditions caused by infection of the central nervous system with proteinaceous agents called prions, including sporadic, variant, and iatrogenic Creutzfeldt-Jakob disease; kuru; inherited prion disease; sheep scrapie; bovine spongiform encephalopathy; and chronic wasting disease. Prions are composed of misfolded and multimeric forms of the normal cellular prion protein (PrP). Prion diseases require host expression of the prion protein gene (PRNP) and a range of other cellular functions to support their propagation and toxicity. Inherited forms of prion disease are caused by mutation of PRNP, whereas acquired and sporadically occurring mammalian prion diseases are controlled by powerful genetic risk and modifying factors. Whereas some PrP amino acid variants cause the disease, others confer protection, dramatically altered incubation times, or changes in the clinical phenotype. Multiple mechanisms, including interference with homotypic protein interactions and the selection of the permissible prion strains in a host, play a role. Several non-PRNP factors have now been uncovered that provide insights into pathways of disease susceptibility or neurotoxicity.

Whole genome sequencing of 4,787 individuals identifies gene-based rare variants in age-related macular degeneration.

Genome-wide association studies have contributed extensively to the discovery of disease-associated common variants. However, the genetic contribution to complex traits is still largely difficult to interpret. We report a genome-wide association study of 2394 cases and 2393 controls for age-related macular degeneration (AMD) via whole-genome sequencing, with 46.9 million genetic variants. Our study reveals significant single-variant association signals at four loci and independent gene-based signals in CFH, C2, C3, and NRTN. Using data from the Exome Aggregation Consortium (ExAC) for a gene-based test, we demonstrate an enrichment of predicted rare loss-of-function variants in CFH, CFI, and an as-yet unreported gene in AMD, ORMDL2. Our method of using a large variant list without individual-level genotypes as an external reference provides a flexible and convenient approach to leverage the publicly available variant datasets to augment the search for rare variant associations, which can explain additional disease risk in AMD.

The genetic dissection of fetal haemoglobin persistence in sickle cell disease in Nigeria.

The clinical severity of sickle cell disease (SCD) is strongly influenced by the level of fetal haemoglobin (HbF) persistent in each patient. Three major HbF loci (BCL11A, HBS1L-MYB, and Xmn1-HBG2) have been reported, but a considerable hidden heritability remains. We conducted a genome-wide association study for HbF levels in 1006 Nigerian patients with SCD (HbSS/HbSß0), followed by a replication and meta-analysis exercise in four independent SCD cohorts (3,582 patients). To dissect association signals at the major loci, we performed stepwise conditional and haplotype association analyses and included public functional annotation datasets. Association signals were detected for BCL11A (lead SNP rs6706648, ß = -0.39, P = 4.96 × 10-34) and HBS1L-MYB (lead SNP rs61028892, ß = 0.73, P = 1.18 × 10-9), whereas the variant allele for Xmn1-HBG2 was found to be very rare. In addition, we detected three putative new trait-associated regions. Genetically, dissecting the two major loci BCL11A and HBS1L-MYB, we defined trait-increasing haplotypes (P < 0.0001) containing so far unidentified causal variants. At BCL11A, in addition to a haplotype harbouring the putative functional variant rs1427407-'T', we identified a second haplotype, tagged by the rs7565301-'A' allele, where a yet-to-be-discovered causal DNA variant may reside. Similarly, at HBS1L-MYB, one HbF-increasing haplotype contains the likely functional small indel rs66650371, and a second tagged by rs61028892-'C' is likely to harbour a presently unknown functional allele. Together, variants at BCL11A and HBS1L-MYB SNPs explained 24.1% of the trait variance. Our findings provide a path for further investigation of the causes of variable fetal haemoglobin persistence in sickle cell disease.

A multi-ancestry genome-wide association study in type 1 diabetes.

Type 1 diabetes (T1D) is an autoimmune disease caused by destruction of the pancreatic ß-cells. Genome-wide association (GWAS) and fine mapping studies have been conducted mainly in European ancestry (EUR) populations. We performed a multi-ancestry GWAS to identify SNPs and HLA alleles associated with T1D risk and age at onset. EUR families (N = 3223), and unrelated individuals of African (AFR, N = 891) and admixed (Hispanic/Latino) ancestry (AMR, N = 308) were genotyped using the Illumina HumanCoreExome BeadArray, with imputation to the TOPMed reference panel. The Multi-Ethnic HLA reference panel was utilized to impute HLA alleles and amino acid residues. Logistic mixed models (T1D risk) and frailty models (age at onset) were used for analysis. In GWAS meta-analysis, seven loci were associated with T1D risk at genome-wide significance: PTPN22, HLA-DQA1, IL2RA, RNLS, INS, IKZF4-RPS26-ERBB3, and SH2B3, with four associated with T1D age at onset (PTPN22, HLA-DQB1, INS, and ERBB3). AFR and AMR meta-analysis revealed NRP1 as associated with T1D risk and age at onset, although NRP1 variants were not associated in EUR ancestry. In contrast, the PTPN22 variant was significantly associated with risk only in EUR ancestry. HLA alleles and haplotypes most significantly associated with T1D risk in AFR and AMR ancestry differed from that seen in EUR ancestry; in addition, the HLA-DRB1*08:02-DQA1*04:01-DQB1*04:02 haplotype was 'protective' in AMR while HLA-DRB1*08:01-DQA1*04:01-DQB1*04:02 haplotype was 'risk' in EUR ancestry, differing only at HLA-DRB1*08. These results suggest that much larger sample sizes in non-EUR populations are required to capture novel loci associated with T1D risk.

Identifying X-chromosome variants associated with age-related macular degeneration.

PURPOSE: In genome-wide association studies (GWAS), X chromosome (ChrX) variants are often not investigated. Sex-specific effects and ChrX-specific quality control (QC) are needed to examine these effects. Previous GWAS identified 52 autosomal variants associated with age-related macular degeneration (AMD) via the International AMD Genomics Consortium (IAMDGC), but did not analyze ChrX. Therefore¸ our goal was to investigate ChrX variants for association with AMD. METHODS: We genotyped 29 629 non-Hispanic White (NHW) individuals (M/F:10404/18865; AMD12,087/14723) via a custom chip and imputed after ChrX-specific QC (XWAS 3.0) using the Michigan Imputation Server. Imputation generated 1 221 623 variants on ChrX. Age, informative PCs, and subphenotypes were covariates for logistic association analyses with Fisher's correction. Gene/pathway analyses were performed with VEGAS, GSEASNP, ICSNPathway, DAVID, and mirPath. RESULTS: Logistic association on NHW individuals with sex correction identified variants in/near the genes SLITRK4, ARHGAP6, FGF13 and DMD associated with AMD (P < 1 × 10-6,Fisher's combined-corrected). Association testing of the subphenotypes of choroidal neovascularization and geographic atrophy (GA), identified variants in DMD associated with GA (P < 1 × 10-6, Fisher's combined-corrected). Via gene-based analysis with VEGAS, several genes were associated with AMD (P < 0.05, both truncated tail strength/truncated product P) including SLITRK4 and BHLHB9. Pathway analysis using GSEASNP and DAVID identified genes associated with nervous system development (FDR: P:0.02), and blood coagulation (FDR: P:0.03). Variants in the region of a microRNA (miR) were associated with AMD (P < 0.05, truncated tail strength/truncated product P). Via DIANA mirPath analysis, downstream targets of miRs showed association with brain disorders and fatty acid elongation (P < 0.05). A long noncoding RNA on ChrX near the DMD locus was also associated with AMD (P = 4 × 10-7). Epistatic analysis (t-statistic) for a quantitative trait of AMD vs control including covariates found a suggestive association in the XG gene (P = 2 × 10^-5). CONCLUSIONS: Analysis of ChrX variation identifies several potential new locifor AMD risk and these variants nominate novel AMD pathways. Further analysis is needed to refine these results and to understand their biological significance and relationship with AMD development in worldwide populations.

LDAK-GBAT: Fast and powerful gene-based association testing using summary statistics.

We present LDAK-GBAT, a tool for gene-based association testing using summary statistics from genome-wide association studies that is computationally efficient, produces well-calibrated p values, and is significantly more powerful than existing tools. LDAK-GBAT takes approximately 30 min to analyze imputed data (2.9M common, genic SNPs), requiring less than 10 Gb memory. It shows good control of type 1 error given an appropriate reference panel. Across 109 phenotypes (82 from the UK Biobank, 18 from the Million Veteran Program, and nine from the Psychiatric Genetics Consortium), LDAK-GBAT finds on average 19% (SE: 1%) more significant genes than the existing tool sumFREGAT-ACAT, with even greater gains in comparison with MAGMA, GCTA-fastBAT, sumFREGAT-SKAT-O, and sumFREGAT-PCA.

High-throughput identification of regulatory elements and functional assays to uncover susceptibility genes for nasopharyngeal carcinoma.

Large-scale genetic association studies have identified multiple susceptibility loci for nasopharyngeal carcinoma (NPC), but the underlying biological mechanisms remain to be explored. To gain insights into the genetic etiology of NPC, we conducted a follow-up study encompassing 6,907 cases and 10,472 controls and identified two additional NPC susceptibility loci, 9q22.33 (rs1867277; OR = 0.74, 95% CI = 0.68-0.81, p = 3.08 × 10-11) and 17q12 (rs226241; OR = 1.42, 95% CI = 1.26-1.60, p = 1.62 × 10-8). The two additional loci, together with two previously reported genome-wide significant loci, 5p15.33 and 9p21.3, were investigated by high-throughput sequencing for chromatin accessibility, histone modification, and promoter capture Hi-C (PCHi-C) profiling. Using luciferase reporter assays and CRISPR interference (CRISPRi) to validate the functional profiling, we identified PHF2 at locus 9q22.33 as a susceptibility gene. PHF2 encodes a histone demethylase and acts as a tumor suppressor. The risk alleles of the functional SNPs reduced the expression of the target gene PHF2 by inhibiting the enhancer activity of its long-range (4.3 Mb) cis-regulatory element, which promoted proliferation of NPC cells. In addition, we identified CDKN2B-AS1 as a susceptibility gene at locus 9p21.3, and the NPC risk allele of the functional SNP rs2069418 promoted the expression of CDKN2B-AS1 by increasing its enhancer activity. The overexpression of CDKN2B-AS1 facilitated proliferation of NPC cells. In summary, we identified functional SNPs and NPC susceptibility genes, which provides additional explanations for the genetic association signals and helps to uncover the underlying genetic etiology of NPC development.

Dissecting the polygenic basis of atherosclerosis via disease-associated cell state signatures.

Coronary artery disease (CAD) is a pandemic disease where up to half of the risk is explained by genetic factors. Advanced insights into the genetic basis of CAD require deeper understanding of the contributions of different cell types, molecular pathways, and genes to disease heritability. Here, we investigate the biological diversity of atherosclerosis-associated cell states and interrogate their contribution to the genetic risk of CAD by using single-cell and bulk RNA sequencing (RNA-seq) of mouse and human lesions. We identified 12 disease-associated cell states that we characterized further by gene set functional profiling, ligand-receptor prediction, and transcription factor inference. Importantly, Vcam1+ smooth muscle cell state genes contributed most to SNP-based heritability of CAD. In line with this, genetic variants near smooth muscle cell state genes and regulatory elements explained the largest fraction of CAD-risk variance between individuals. Using this information for variant prioritization, we derived a hybrid polygenic risk score (PRS) that demonstrated improved performance over a classical PRS. Our results provide insights into the biological mechanisms associated with CAD risk, which could make a promising contribution to precision medicine and tailored therapeutic interventions in the future.

Genome-wide association study of esophageal squamous cell cancer identifies shared and distinct risk variants in African and Chinese populations.

Esophageal squamous cell carcinoma (ESCC) has a high disease burden in sub-Saharan Africa and has a very poor prognosis. Genome-wide association studies (GWASs) of ESCC in predominantly East Asian populations indicate a substantial genetic contribution to its etiology, but no genome-wide studies have been done in populations of African ancestry. Here, we report a GWAS in 1,686 African individuals with ESCC and 3,217 population-matched control individuals to investigate its genetic etiology. We identified a genome-wide-significant risk locus on chromosome 9 upstream of FAM120A (rs12379660, p = 4.58 × 10-8, odds ratio = 1.28, 95% confidence interval = 1.22-1.34), as well as a potential African-specific risk locus on chromosome 2 (rs142741123, p = 5.49 × 10-8) within MYO1B. FAM120A is a component of oxidative stress-induced survival signals, and the associated variants at the FAM120A locus co-localized with highly significant cis-eQTLs in FAM120AOS in both esophageal mucosa and esophageal muscularis tissue. A trans-ethnic meta-analysis was then performed with the African ESCC study and a Chinese ESCC study in a combined total of 3,699 ESCC-affected individuals and 5,918 control individuals, which identified three genome-wide-significant loci on chromosome 9 at FAM120A (rs12379660, pmeta = 9.36 × 10-10), chromosome 10 at PLCE1 (rs7099485, pmeta = 1.48 × 10-8), and chromosome 22 at CHEK2 (rs1033667, pmeta = 1.47 × 10-9). This indicates the existence of both shared and distinct genetic risk loci for ESCC in African and Asian populations. Our GWAS of ESCC conducted in a population of African ancestry indicates a substantial genetic contribution to ESCC risk in Africa.