Primary angle closure glaucoma (PACG) is a major cause of blindness worldwide. We conducted a genome-wide association study (GWAS) followed by replication in a combined total of 10,503 PACG cases and 29,567 controls drawn from 24 countries across Asia, Australia, Europe, North America, and South America. We observed significant evidence of disease association at five new genetic loci upon meta-analysis of all patient collections. These loci are at EPDR1 rs3816415 (odds ratio (OR) = 1.24, P = 5.94 × 10(-15)), CHAT rs1258267 (OR = 1.22, P = 2.85 × 10(-16)), GLIS3 rs736893 (OR = 1.18, P = 1.43 × 10(-14)), FERMT2 rs7494379 (OR = 1.14, P = 3.43 × 10(-11)), and DPM2-FAM102A rs3739821 (OR = 1.15, P = 8.32 × 10(-12)). We also confirmed significant association at three previously described loci (P < 5 × 10(-8) for each sentinel SNP at PLEKHA7, COL11A1, and PCMTD1-ST18), providing new insights into the biology of PACG.
Primary angle closure glaucoma (PACG) is a major cause of blindness worldwide. We conducted a genome-wide association study including 1,854 PACG cases and 9,608 controls across 5 sample collections in Asia. Replication experiments were conducted in 1,917 PACG cases and 8,943 controls collected from a further 6 sample collections. We report significant associations at three new loci: rs11024102 in PLEKHA7 (per-allele odds ratio (OR)=1.22; P=5.33×10(-12)), rs3753841 in COL11A1 (per-allele OR=1.20; P=9.22×10(-10)) and rs1015213 located between PCMTD1 and ST18 on chromosome 8q (per-allele OR=1.50; P=3.29×10(-9)). Our findings, accumulated across these independent worldwide collections, suggest possible mechanisms explaining the pathogenesis of PACG.
We conducted a genome-wide association study of oral cavity and pharyngeal cancer in 6,034 cases and 6,585 controls from Europe, North America and South America. We detected eight significantly associated loci (P < 5 × 10-8), seven of which are new for these cancer sites. Oral and pharyngeal cancers combined were associated with loci at 6p21.32 (rs3828805, HLA-DQB1), 10q26.13 (rs201982221, LHPP) and 11p15.4 (rs1453414, OR52N2-TRIM5). Oral cancer was associated with two new regions, 2p23.3 (rs6547741, GPN1) and 9q34.12 (rs928674, LAMC3), and with known cancer-related loci-9p21.3 (rs8181047, CDKN2B-AS1) and 5p15.33 (rs10462706, CLPTM1L). Oropharyngeal cancer associations were limited to the human leukocyte antigen (HLA) region, and classical HLA allele imputation showed a protective association with the class II haplotype HLA-DRB1*1301-HLA-DQA1*0103-HLA-DQB1*0603 (odds ratio (OR) = 0.59, P = 2.7 × 10-9). Stratified analyses on a subgroup of oropharyngeal cases with information available on human papillomavirus (HPV) status indicated that this association was considerably stronger in HPV-positive (OR = 0.23, P = 1.6 × 10-6) than in HPV-negative (OR = 0.75, P = 0.16) cancers.
Aldosterone-producing adenomas (APAs) are the commonest curable cause of hypertension. Most have gain-of-function somatic mutations of ion channels or transporters. Herein we report the discovery, replication and phenotype of mutations in the neuronal cell adhesion gene CADM1. Independent whole exome sequencing of 40 and 81 APAs found intramembranous p.Val380Asp or p.Gly379Asp variants in two patients whose hypertension and periodic primary aldosteronism were cured by adrenalectomy. Replication identified two more APAs with each variant (total, n = 6). The most upregulated gene (10- to 25-fold) in human adrenocortical H295R cells transduced with the mutations (compared to wildtype) was CYP11B2 (aldosterone synthase), and biological rhythms were the most differentially expressed process. CADM1 knockdown or mutation inhibited gap junction (GJ)-permeable dye transfer. GJ blockade by Gap27 increased CYP11B2 similarly to CADM1 mutation. Human adrenal zona glomerulosa (ZG) expression of GJA1 (the main GJ protein) was patchy, and annular GJs (sequelae of GJ communication) were less prominent in CYP11B2-positive micronodules than adjacent ZG. Somatic mutations of CADM1 cause reversible hypertension and reveal a role for GJ communication in suppressing physiological aldosterone production.
Most aldosterone-producing adenomas (APAs) have gain-of-function somatic mutations of ion channels or transporters. However, their frequency in aldosterone-producing cell clusters of normal adrenal gland suggests a requirement for codriver mutations in APAs. Here we identified gain-of-function mutations in both CTNNB1 and GNA11 by whole-exome sequencing of 3/41 APAs. Further sequencing of known CTNNB1-mutant APAs led to a total of 16 of 27 (59%) with a somatic p.Gln209His, p.Gln209Pro or p.Gln209Leu mutation of GNA11 or GNAQ. Solitary GNA11 mutations were found in hyperplastic zona glomerulosa adjacent to double-mutant APAs. Nine of ten patients in our UK/Irish cohort presented in puberty, pregnancy or menopause. Among multiple transcripts upregulated more than tenfold in double-mutant APAs was LHCGR, the receptor for luteinizing or pregnancy hormone (human chorionic gonadotropin). Transfections of adrenocortical cells demonstrated additive effects of GNA11 and CTNNB1 mutations on aldosterone secretion and expression of genes upregulated in double-mutant APAs. In adrenal cortex, GNA11/Q mutations appear clinically silent without a codriver mutation of CTNNB1.
Autism spectrum disorder (ASD) is a highly heritable and heterogeneous group of neurodevelopmental phenotypes diagnosed in more than 1% of children. Common genetic variants contribute substantially to ASD susceptibility, but to date no individual variants have been robustly associated with ASD. With a marked sample-size increase from a unique Danish population resource, we report a genome-wide association meta-analysis of 18,381 individuals with ASD and 27,969 controls that identified five genome-wide-significant loci. Leveraging GWAS results from three phenotypes with significantly overlapping genetic architectures (schizophrenia, major depression, and educational attainment), we identified seven additional loci shared with other traits at equally strict significance levels. Dissecting the polygenic architecture, we found both quantitative and qualitative polygenic heterogeneity across ASD subtypes. These results highlight biological insights, particularly relating to neuronal function and corticogenesis, and establish that GWAS performed at scale will be much more productive in the near term in ASD.
Genome-wide association studies (GWAS) have identified 12 epithelial ovarian cancer (EOC) susceptibility alleles. The pattern of association at these loci is consistent in BRCA1 and BRCA2 mutation carriers who are at high risk of EOC. After imputation to 1000 Genomes Project data, we assessed associations of 11 million genetic variants with EOC risk from 15,437 cases unselected for family history and 30,845 controls and from 15,252 BRCA1 mutation carriers and 8,211 BRCA2 mutation carriers (3,096 with ovarian cancer), and we combined the results in a meta-analysis. This new study design yielded increased statistical power, leading to the discovery of six new EOC susceptibility loci. Variants at 1p36 (nearest gene, WNT4), 4q26 (SYNPO2), 9q34.2 (ABO) and 17q11.2 (ATAD5) were associated with EOC risk, and at 1p34.3 (RSPO1) and 6p22.1 (GPX6) variants were specifically associated with the serous EOC subtype, all with P < 5 × 10(-8). Incorporating these variants into risk assessment tools will improve clinical risk predictions for BRCA1 and BRCA2 mutation carriers.
Allelic heterogeneity in disease-causing genes presents a substantial challenge to the translation of genomic variation into clinical practice. Few of the almost 2,000 variants in the cystic fibrosis transmembrane conductance regulator gene CFTR have empirical evidence that they cause cystic fibrosis. To address this gap, we collected both genotype and phenotype data for 39,696 individuals with cystic fibrosis in registries and clinics in North America and Europe. In these individuals, 159 CFTR variants had an allele frequency of ł0.01%. These variants were evaluated for both clinical severity and functional consequence, with 127 (80%) meeting both clinical and functional criteria consistent with disease. Assessment of disease penetrance in 2,188 fathers of individuals with cystic fibrosis enabled assignment of 12 of the remaining 32 variants as neutral, whereas the other 20 variants remained of indeterminate effect. This study illustrates that sourcing data directly from well-phenotyped subjects can address the gap in our ability to interpret clinically relevant genomic variation.
Genome-wide association studies have identified breast cancer risk variants in over 150 genomic regions, but the mechanisms underlying risk remain largely unknown. These regions were explored by combining association analysis with in silico genomic feature annotations. We defined 205 independent risk-associated signals with the set of credible causal variants in each one. In parallel, we used a Bayesian approach (PAINTOR) that combines genetic association, linkage disequilibrium and enriched genomic features to determine variants with high posterior probabilities of being causal. Potentially causal variants were significantly over-represented in active gene regulatory regions and transcription factor binding sites. We applied our INQUSIT pipeline for prioritizing genes as targets of those potentially causal variants, using gene expression (expression quantitative trait loci), chromatin interaction and functional annotations. Known cancer drivers, transcription factors and genes in the developmental, apoptosis, immune system and DNA integrity checkpoint gene ontology pathways were over-represented among the highest-confidence target genes.
Genome-wide association studies (GWAS) and large-scale replication studies have identified common variants in 79 loci associated with breast cancer, explaining ∼14% of the familial risk of the disease. To identify new susceptibility loci, we performed a meta-analysis of 11 GWAS, comprising 15,748 breast cancer cases and 18,084 controls together with 46,785 cases and 42,892 controls from 41 studies genotyped on a 211,155-marker custom array (iCOGS). Analyses were restricted to women of European ancestry. We generated genotypes for more than 11 million SNPs by imputation using the 1000 Genomes Project reference panel, and we identified 15 new loci associated with breast cancer at P < 5 × 10(-8). Combining association analysis with ChIP-seq chromatin binding data in mammary cell lines and ChIA-PET chromatin interaction data from ENCODE, we identified likely target genes in two regions: SETBP1 at 18q12.3 and RNF115 and PDZK1 at 1q21.1. One association appears to be driven by an amino acid substitution encoded in EXO1.
We analyzed 3,872 common genetic variants across the ESR1 locus (encoding estrogen receptor α) in 118,816 subjects from three international consortia. We found evidence for at least five independent causal variants, each associated with different phenotype sets, including estrogen receptor (ER(+) or ER(-)) and human ERBB2 (HER2(+) or HER2(-)) tumor subtypes, mammographic density and tumor grade. The best candidate causal variants for ER(-) tumors lie in four separate enhancer elements, and their risk alleles reduce expression of ESR1, RMND1 and CCDC170, whereas the risk alleles of the strongest candidates for the remaining independent causal variant disrupt a silencer element and putatively increase ESR1 and RMND1 expression.
Breast cancer susceptibility variants frequently show heterogeneity in associations by tumor subtype1-3. To identify novel loci, we performed a genome-wide association study including 133,384 breast cancer cases and 113,789 controls, plus 18,908 BRCA1 mutation carriers (9,414 with breast cancer) of European ancestry, using both standard and novel methodologies that account for underlying tumor heterogeneity by estrogen receptor, progesterone receptor and human epidermal growth factor receptor 2 status and tumor grade. We identified 32 novel susceptibility loci (P
Genome-wide association studies (GWAS) and fine-mapping efforts to date have identified more than 100 prostate cancer (PrCa)-susceptibility loci. We meta-analyzed genotype data from a custom high-density array of 46,939 PrCa cases and 27,910 controls of European ancestry with previously genotyped data of 32,255 PrCa cases and 33,202 controls of European ancestry. Our analysis identified 62 novel loci associated (P C, p.Pro1054Arg) in ATM and rs2066827 (OR = 1.06; P = 2.3 × 10-9; T>G, p.Val109Gly) in CDKN1B. The combination of all loci captured 28.4% of the PrCa familial relative risk, and a polygenic risk score conferred an elevated PrCa risk for men in the ninetieth to ninety-ninth percentiles (relative risk = 2.69; 95% confidence interval (CI): 2.55-2.82) and first percentile (relative risk = 5.71; 95% CI: 5.04-6.48) risk stratum compared with the population average. These findings improve risk prediction, enhance fine-mapping, and provide insight into the underlying biology of PrCa1.
In the version of this article initially published, the name of author Manuela Gago-Dominguez was misspelled as Manuela Gago Dominguez. The error has been corrected in the HTML and PDF version of the article.
Prostate cancer is a highly heritable disease with large disparities in incidence rates across ancestry populations. We conducted a multiancestry meta-analysis of prostate cancer genome-wide association studies (107,247 cases and 127,006 controls) and identified 86 new genetic risk variants independently associated with prostate cancer risk, bringing the total to 269 known risk variants. The top genetic risk score (GRS) decile was associated with odds ratios that ranged from 5.06 (95% confidence interval (CI), 4.84-5.29) for men of European ancestry to 3.74 (95% CI, 3.36-4.17) for men of African ancestry. Men of African ancestry were estimated to have a mean GRS that was 2.18-times higher (95% CI, 2.14-2.22), and men of East Asian ancestry 0.73-times lower (95% CI, 0.71-0.76), than men of European ancestry. These findings support the role of germline variation contributing to population differences in prostate cancer risk, with the GRS offering an approach for personalized risk prediction.
The transferability and clinical value of genetic risk scores (GRSs) across populations remain limited due to an imbalance in genetic studies across ancestrally diverse populations. Here we conducted a multi-ancestry genome-wide association study of 156,319 prostate cancer cases and 788,443 controls of European, African, Asian and Hispanic men, reflecting a 57% increase in the number of non-European cases over previous prostate cancer genome-wide association studies. We identified 187 novel risk variants for prostate cancer, increasing the total number of risk variants to 451. An externally replicated multi-ancestry GRS was associated with risk that ranged from 1.8 (per standard deviation) in African ancestry men to 2.2 in European ancestry men. The GRS was associated with a greater risk of aggressive versus non-aggressive disease in men of African ancestry (P = 0.03). Our study presents novel prostate cancer susceptibility loci and a GRS with effective risk stratification across ancestry groups.
Galloway-Mowat syndrome (GAMOS) is an autosomal-recessive disease characterized by the combination of early-onset nephrotic syndrome (SRNS) and microcephaly with brain anomalies. Here we identified recessive mutations in OSGEP, TP53RK, TPRKB, and LAGE3, genes encoding the four subunits of the KEOPS complex, in 37 individuals from 32 families with GAMOS. CRISPR-Cas9 knockout in zebrafish and mice recapitulated the human phenotype of primary microcephaly and resulted in early lethality. Knockdown of OSGEP, TP53RK, or TPRKB inhibited cell proliferation, which human mutations did not rescue. Furthermore, knockdown of these genes impaired protein translation, caused endoplasmic reticulum stress, activated DNA-damage-response signaling, and ultimately induced apoptosis. Knockdown of OSGEP or TP53RK induced defects in the actin cytoskeleton and decreased the migration rate of human podocytes, an established intermediate phenotype of SRNS. We thus identified four new monogenic causes of GAMOS, describe a link between KEOPS function and human disease, and delineate potential pathogenic mechanisms.
Pancreatic cancer is the fourth leading cause of cancer death in the developed world. Both inherited high-penetrance mutations in BRCA2 (ref. 2), ATM, PALB2 (ref. 4), BRCA1 (ref. 5), STK11 (ref. 6), CDKN2A and mismatch-repair genes and low-penetrance loci are associated with increased risk. To identify new risk loci, we performed a genome-wide association study on 9,925 pancreatic cancer cases and 11,569 controls, including 4,164 newly genotyped cases and 3,792 controls in 9 studies from North America, Central Europe and Australia. We identified three newly associated regions: 17q25.1 (LINC00673, rs11655237, odds ratio (OR) = 1.26, 95% confidence interval (CI) = 1.19-1.34, P = 1.42 × 10(-14)), 7p13 (SUGCT, rs17688601, OR = 0.88, 95% CI = 0.84-0.92, P = 1.41 × 10(-8)) and 3q29 (TP63, rs9854771, OR = 0.89, 95% CI = 0.85-0.93, P = 2.35 × 10(-8)). We detected significant association at 2p13.3 (ETAA1, rs1486134, OR = 1.14, 95% CI = 1.09-1.19, P = 3.36 × 10(-9)), a region with previous suggestive evidence in Han Chinese. We replicated previously reported associations at 9q34.2 (ABO), 13q22.1 (KLF5), 5p15.33 (TERT and CLPTM1), 13q12.2 (PDX1), 1q32.1 (NR5A2), 7q32.3 (LINC-PINT), 16q23.1 (BCAR1) and 22q12.1 (ZNRF3). Our study identifies new loci associated with pancreatic cancer risk.
Transforming growth factor (TGF)-β1 (encoded by TGFB1) is the prototypic member of the TGF-β family of 33 proteins that orchestrate embryogenesis, development and tissue homeostasis1,2. Following its discovery 3 , enormous interest and numerous controversies have emerged about the role of TGF-β in coordinating the balance of pro- and anti-oncogenic properties4,5, pro- and anti-inflammatory effects 6 , or pro- and anti-fibrinogenic characteristics 7 . Here we describe three individuals from two pedigrees with biallelic loss-of-function mutations in the TGFB1 gene who presented with severe infantile inflammatory bowel disease (IBD) and central nervous system (CNS) disease associated with epilepsy, brain atrophy and posterior leukoencephalopathy. The proteins encoded by the mutated TGFB1 alleles were characterized by impaired secretion, function or stability of the TGF-β1-LAP complex, which is suggestive of perturbed bioavailability of TGF-β1. Our study shows that TGF-β1 has a critical and nonredundant role in the development and homeostasis of intestinal immunity and the CNS in humans.