To further resolve the genetic architecture of the inflammatory bowel diseases ulcerative colitis and Crohn's disease, we sequenced the whole genomes of 4,280 patients at low coverage and compared them to 3,652 previously sequenced population controls across 73.5 million variants. We then imputed from these sequences into new and existing genome-wide association study cohorts and tested for association at ∼12 million variants in a total of 16,432 cases and 18,843 controls. We discovered a 0.6% frequency missense variant in ADCY7 that doubles the risk of ulcerative colitis. Despite good statistical power, we did not identify any other new low-frequency risk variants and found that such variants explained little heritability. We detected a burden of very rare, damaging missense variants in known Crohn's disease risk genes, suggesting that more comprehensive sequencing studies will continue to improve understanding of the biology of complex diseases.
In genome-wide association studies (GWAS) for thousands of phenotypes in large biobanks, most binary traits have substantially fewer cases than controls. Both of the widely used approaches, the linear mixed model and the recently proposed logistic mixed model, perform poorly; they produce large type I error rates when used to analyze unbalanced case-control phenotypes. Here we propose a scalable and accurate generalized mixed model association test that uses the saddlepoint approximation to calibrate the distribution of score test statistics. This method, SAIGE (Scalable and Accurate Implementation of GEneralized mixed model), provides accurate P values even when case-control ratios are extremely unbalanced. SAIGE uses state-of-art optimization strategies to reduce computational costs; hence, it is applicable to GWAS for thousands of phenotypes by large biobanks. Through the analysis of UK Biobank data of 408,961 samples from white British participants with European ancestry for > 1,400 binary phenotypes, we show that SAIGE can efficiently analyze large sample data, controlling for unbalanced case-control ratios and sample relatedness.
Epidemiological studies demonstrate an association between asthma and mental health disorders, although little is known about the shared genetics and causality of this association. Thus, we aim to investigate shared genetic and the causal link between asthma and mental health disorders.We conducted a large-scale genome-wide cross-trait association study to investigate genetic overlap between asthma from UK Biobank and 8 mental health disorders from Psychiatric Genomics Consortium, including: attention deficit hyperactivity disorder (ADHD), anxiety disorder (ANX), autism spectrum disorder, bipolar disorder, eating disorder, major depressive disorder (MDD), posttraumatic stress disorder, and schizophrenia, with a sample size of 7556 to 446 032.In the single trait genome-wide association analysis, we replicated 130 and discovered 31 novel independent loci that are associated with asthma. We identified that ADHD, ANX and MDD have strong genetic correlation with asthma at the genome-wide level. Cross-trait meta-analysis identified 7 loci jointly associated with asthma and ADHD, 1 loci with asthma and ANX and 10 loci with asthma and MDD. Functional analysis revealed that the identified variants regulated gene expression in major tissues belonging to exocrine/endocrine, digestive, respiratory and hemic/immune system. Mendelian randomisation analyses suggested that ADHD and MDD (including 6.7% samples overlap with asthma) might increase the risk of asthma.This large-scale genome-wide cross-trait analysis identified shared genetics and potential causal links between asthma and three mental health disorders (ADHD, ANX, and MDD). Such shared genetics implicate potential new biological functions that are in common among them.
Asthma
Asthma and attention deficit hyperactivity disorder
Cigarette smoking is a leading cause of preventable mortality worldwide. Nicotine dependence, which reduces the likelihood of quitting smoking, is a heritable trait with firmly established associations with sequence variants in nicotine acetylcholine receptor genes and at other loci. To search for additional loci, we conducted a genome-wide association study (GWAS) meta-analysis of nicotine dependence, totaling 38,602 smokers (28,677 Europeans/European Americans and 9925 African Americans) across 15 studies. In this largest-ever GWAS meta-analysis for nicotine dependence and the largest-ever cross-ancestry GWAS meta-analysis for any smoking phenotype, we reconfirmed the well-known CHRNA5-CHRNA3-CHRNB4 genes and further yielded a novel association in the DNA methyltransferase gene DNMT3B. The intronic DNMT3B rs910083-C allele (frequency=44-77%) was associated with increased risk of nicotine dependence at P=3.7 × 10-8 (odds ratio (OR)=1.06 and 95% confidence interval (CI)=1.04-1.07 for severe vs mild dependence). The association was independently confirmed in the UK Biobank (N=48,931) using heavy vs never smoking as a proxy phenotype (P=3.6 × 10-4, OR=1.05, and 95% CI=1.02-1.08). Rs910083-C is also associated with increased risk of squamous cell lung carcinoma in the International Lung Cancer Consortium (N=60,586, meta-analysis P=0.0095, OR=1.05, and 95% CI=1.01-1.09). Moreover, rs910083-C was implicated as a cis-methylation quantitative trait locus (QTL) variant associated with higher DNMT3B methylation in fetal brain (N=166, P=2.3 × 10-26) and a cis-expression QTL variant associated with higher DNMT3B expression in adult cerebellum from the Genotype-Tissue Expression project (N=103, P=3.0 × 10-6) and the independent Brain eQTL Almanac (N=134, P=0.028). This novel DNMT3B cis-acting QTL variant highlights the importance of genetically influenced regulation in brain on the risks of nicotine dependence, heavy smoking and consequent lung cancer.
The production of antibodies to anti-tumor necrosis factor alpha (TNF) agents is one of the main causes of treatment failure in Crohn's disease (CD). To date, however, the contribution of genetics to anti-TNF immunogenicity in CD is still unknown. The objective of the present study was to identify genetic variation associated with anti-TNF immunogenicity in CD. We performed a two-stage genome-wide association study in a cohort of 96 and 123 adalimumab-treated patients, respectively. In the discovery stage, we identified a genome-wide significant association between the CD96 locus and the production of antibodies to anti-TNF treatment (P = 1.88e-09). This association was validated in the replication stage (P < 0.05). The risk allele for anti-TNF immunogenicity was found to be also associated with a lack of response to anti-TNF therapy (P = 0.019). These findings represent an important step toward the understanding of the immunogenicity-based mechanisms that underlie anti-TNF response in CD.
Anti-adalimumab antibody production in response to treatment in Crohn's disease
Age-related hearing impairment (ARHI) is the most common sensory impairment in the aging population; a third of individuals are affected by disabling hearing loss by the age of 65. It causes social isolation and depression and has recently been identified as a risk factor for dementia. The genetic risk factors and underlying pathology of ARHI are largely unknown, meaning that targets for new therapies remain elusive, yet heritability estimates range between 35% and 55%. We performed genome-wide association studies (GWASs) for two self-reported hearing phenotypes, using more than 250,000 UK Biobank (UKBB) volunteers aged between 40 and 69 years. Forty-four independent genome-wide significant loci (p < 5E-08) were identified, considerably increasing the number of established trait loci. Thirty-four loci are novel associations with hearing loss of any form, and only one of the ten known hearing loci has a previously reported association with an ARHI-related trait. Gene sets from these loci are enriched in auditory processes such as synaptic activities, nervous system processes, inner ear morphology, and cognition, while genetic correlation analysis revealed strong positive correlations with multiple personality and psychological traits for the first time. Immunohistochemistry for protein localization in adult mouse cochlea implicate metabolic, sensory, and neuronal functions for NID2, CLRN2, and ARHGEF28. These results provide insight into the genetic landscape underlying ARHI, opening up novel therapeutic targets for further investigation. In a wider context, our study also highlights the viability of using self-report phenotypes for genetic discovery in very large samples when deep phenotyping is unavailable.
BACKGROUND: Adjuvant radiotherapy (RT) can increase the risk of developing pain; however, the molecular mechanisms of RT-related pain remain unclear. The current study aimed to identify susceptibility loci and enriched pathways for clinically relevant acute post-RT pain, defined as having moderate to severe pain (pain score ≥ 4) at the completion of RT. METHODS: We conducted a genome-wide association study (GWAS) with 1,344,832 single-nucleotide polymorphisms (SNPs), a gene-based analysis using PLINK set-based tests of 19,621 genes, and a functional enrichment analysis of a gene list of 875 genes with p < 0.05 using NIH DAVID functional annotation module with KEGG pathways and GO terms (n = 380) among 1112 breast cancer patients. RESULTS: About 29% of patients reported acute post-RT pain. None of SNPs nor genes reached genome-wide significant level. Four SNPs showed suggestive associations with post-RT pain; rs16970540 in RFFL or near the LIG3 gene (p = 1.7 × 10-6), rs4584690, and rs7335912 in ABCC4/MPR4 gene (p = 5.5 × 10-6 and p = 7.8 × 10-6, respectively), and rs73633565 in EGFL6 gene (p = 8.1 × 10-6). Gene-based analysis suggested the potential involvement of neurotransmitters, olfactory receptors, and cytochrome P450 in post-RT pain, whereas functional analysis showed glucuronidation (FDR-adjusted p value = 9.46 × 10-7) and olfactory receptor activities (FDR-adjusted p value = 0.032) as the most significantly enriched biological features. CONCLUSIONS: This is the first GWAS suggesting that post-RT pain is a complex polygenic trait influenced by many biological processes and functions such as glucuronidation and olfactory receptor activities. If validated in larger populations, the results can provide biological targets for pain management to improve cancer patients' quality of life. Additionally, these genes can be further tested as predictive biomarkers for personalized pain management.
Genome-wide association studies (GWAS) have found hundreds of novel loci associated with full blood count (FBC) phenotypes. However, most of these studies were performed in a single phenotype framework without putting into consideration the clinical relatedness among traits. In this work, in addition to the standard univariate GWAS, we also use two different multivariate methods to perform the first multiple traits GWAS of FBC traits in ?7000 individuals from the Ugandan General Population Cohort (GPC). We started by performing the standard univariate GWAS approach. We then performed our first multivariate method, in this approach, we tested for marker associations with 15 FBC traits simultaneously in a multivariate mixed model implemented in GEMMA while accounting for the relatedness of individuals and pedigree structures, as well as population substructure. In this analysis, we provide a framework for the combination of multiple phenotypes in multivariate GWAS analysis and show evidence of multi-collinearity whenever the correlation between traits exceeds the correlation coefficient threshold of r2 >=0.75. This approach identifies two known and one novel loci. In the second multivariate method, we applied principal component analysis (PCA) to the same 15 correlated FBC traits. We then tested for marker associations with each PC in univariate linear mixed models implemented in GEMMA. We show that the FBC composite phenotype as assessed by each PC expresses information that is not completely encapsulated by the individual FBC traits, as this approach identifies three known and five novel loci that were not identified using both the standard univariate and multivariate GWAS methods. Across both multivariate methods, we identified six novel loci. As a proof of concept, both multivariate methods also identified known loci, HBB and ITFG3. The two multivariate methods show that multivariate genotype-phenotype methods increase power and identify novel genotype-phenotype associations not found with the standard univariate GWAS in the same dataset.
Immunoglobulins are the effector molecules of the adaptive humoral immune system. In a genome-wide association study of 19,219 individuals, we found 38 new variants and replicated 5 known variants associating with IgA, IgG or IgM levels or with composite immunoglobulin traits, accounted for by 32 loci. Variants at these loci also affect the risk of autoimmune diseases and blood malignancies and influence blood cell development. Notable associations include a rare variant at RUNX3 decreasing IgA levels by shifting isoform proportions (rs188468174[C>T]: P = 8.3 × 10-55, β = -0.90 s.d.), a rare in-frame deletion in FCGR2B abolishing IgG binding to the encoded receptor (p.Asn106del: P = 4.2 × 10-8, β = 1.03 s.d.), four IGH locus variants influencing class switching, and ten new associations with the HLA region. Our results provide new insight into the regulation of humoral immunity.
(1) Background: In a previous study, we found that two phenotypes related to platelet reactivity, measured with the PFA-100 system, were highly heritable. The aim of the present study was to identify genetic determinants that influence the variability of these phenotypes: closure time of collagen-ADP (Col-ADP) and of collagen-epinephrine (Col-Epi). (2) Methods: As part of the GAIT-2 (Genetic Analysis of Idiopathic Thrombophilia (2) Project, 935 individuals from 35 large Spanish families were studied. A genome-wide association study (GWAS) with ? 10 M single nucleotide polymorphisms (SNPs) was carried out with Col-ADP and Col-Epi phenotypes. (3) Results: The study yielded significant genetic signals that mapped to the ABO locus. After adjusting both phenotypes for the ABO genotype, these signals disappeared. After adjusting for von Willebrand factor (VWF) or for coagulation factor VIII (FVIII), the significant signals disappeared totally for Col-Epi phenotype but only partially for Col-ADP phenotype. (4) Conclusion: Our results suggest that the ABO locus exerts the main genetic influence on PFA-100 phenotypes. However, while the effect of the ABO locus on Col-Epi phenotype is mediated through VWF and/or FVIII, the effect of the ABO locus on Col-ADP phenotype is partly produced through VWF and/or FVIII, and partly through other mechanisms.
Lung function measures are used in the diagnosis of chronic obstructive pulmonary disease. In 38,199 European ancestry individuals, we studied genome-wide association of forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC) and FEV1/FVC with 1000 Genomes Project (phase 1)-imputed genotypes and followed up top associations in 54,550 Europeans. We identify 14 novel loci (P<5 × 10(-8)) in or near ENSA, RNU5F-1, KCNS3, AK097794, ASTN2, LHX3, CCDC91, TBX3, TRIP11, RIN3, TEKT5, LTBP4, MN1 and AP1S2, and two novel signals at known loci NPNT and GPR126, providing a basis for new understanding of the genetic determinants of these traits and pulmonary diseases in which they are altered.
The concentrations of high- and low-density-lipoprotein cholesterol and triglycerides are influenced by smoking, but it is unknown whether genetic associations with lipids may be modified by smoking. We conducted a multi-ancestry genome-wide gene-smoking interaction study in 133,805 individuals with follow-up in an additional 253,467 individuals. Combined meta-analyses identified 13 new loci associated with lipids, some of which were detected only because association differed by smoking status. Additionally, we demonstrate the importance of including diverse populations, particularly in studies of interactions with lifestyle factors, where genomic and lifestyle differences by ancestry may contribute to novel findings.
HDL cholesterol levels in never smokers
LDL cholesterol levels in current smokers
Triglyceride levels x smoking (ever vs never smokers) interaction (2df test)
HDL cholesterol levels x smoking (current vs non smokers) interaction (2df test)
LDL cholesterol levels in never smokers
LDL cholesterol levels x smoking (ever vs never smokers) interaction (2df test)
LDL cholesterol levels x smoking (current vs non smokers) interaction (2df test)
Triglyceride levels in ever smokers
LDL cholesterol levels x smoking (current vs non smokers) interaction (1df test)
Triglyceride levels x smoking (current vs non smokers) interaction (1df test)
LDL cholesterol levels
HDL cholesterol levels x smoking (ever vs never smokers) interaction (2df test)
LDL cholesterol levels x smoking (ever vs never smokers) interaction (1df test)
Triglyceride levels x smoking (current vs non smokers) interaction (2df test)
HDL cholesterol levels in non smokers
HDL cholesterol levels x smoking (ever vs never smokers) interaction (1df test)
Triglyceride levels
HDL cholesterol levels
HDL cholesterol levels x smoking (current vs non smokers) interaction (1df test)
Triglyceride levels x smoking (ever vs never smokers) interaction (1df test)
Primary open angle glaucoma (POAG), a major cause of blindness worldwide, is a complex disease with a significant genetic contribution. We performed Exome Array (Illumina) analysis on 3504 POAG cases and 9746 controls with replication of the most significant findings in 9173 POAG cases and 26 780 controls across 18 collections of Asian, African and European descent. Apart from confirming strong evidence of association at CDKN2B-AS1 (rs2157719 [G], odds ratio [OR] = 0.71, P = 2.81 × 10(-33)), we observed one SNP showing significant association to POAG (CDC7-TGFBR3 rs1192415, ORG-allele = 1.13, Pmeta = 1.60 × 10(-8)). This particular SNP has previously been shown to be strongly associated with optic disc area and vertical cup-to-disc ratio, which are regarded as glaucoma-related quantitative traits. Our study now extends this by directly implicating it in POAG disease pathogenesis.
The United States is experiencing the worst opioid overdose (OpOD) crisis in its history. We carried out a genome-wide association study on OpOD severity among 3 477 opioid-exposed individuals, 1 019 of whom experienced OpODs, including 2 032 European Americans (EAs) (653 overdose cases), and 1 445 African Americans (AAs) (366 overdose cases). Participants were scored 1 to 4 based on their reported overdose status and the number of times that medical treatment was required. Genome-wide association study (GWAS) of EAs and AAs separately resulted in two genome-wide significant (GWS) signals in AAs but none in EAs. The first signal was represented by three closely mapped variants (rs115208233, rs116181528, and rs114077267) located near mucolipin 1 (MCOLN1) and patatin-like phospholipase domain containing 6 (PNPLA6), and the other signal was represented by rs369098800 near dead-box helicase 18 (DDX18). There were no additional GWS signals in the trans-population meta-analysis, so that post-GWAS analysis focused on these loci. In network analysis, MCOLN1 was coexpressed with PNPLA6, but only MCOLN1-associated genes were enriched in functional categories relevant to OpOD, including calcium and cation channel activities; no enrichment was observed for PNPLA6-associated genes. Drug repositioning analysis was carried out in the connectivity map (CMap) database for MCOLN1 (PNPLA6 was not available in CMap) and showed that the opioid agonist drug-induced expression profile is similar to that of MCOLN1 overexpression and yielded the highest-ranked expression profile of 83 drug classes. Thus, MCOLN1 may be a risk gene for OpOD, but replication is needed. This knowledge could be helpful in the identification of drug targets for preventing OpOD.
Importance: Anxiety and stress-related disorders are among the most common mental disorders. Although family and twin studies indicate that both genetic and environmental factors play an important role underlying their etiology, the genetic underpinnings of anxiety and stress-related disorders are poorly understood. Objectives: To estimate the single-nucleotide polymorphism-based heritability of anxiety and stress-related disorders; to identify novel genetic risk variants, genes, or biological pathways; to test for pleiotropic associations with other psychiatric traits; and to evaluate the association of psychiatric comorbidities with genetic findings. Design, Setting, Participants: This genome-wide association study included individuals with various anxiety and stress-related diagnoses and controls derived from the population-based Lundbeck Foundation Initiative for Integrative Psychiatric Research (iPSYCH) study. Lifetime diagnoses of anxiety and stress-related disorders were obtained through the national Danish registers. Genes of interest were further evaluated in mice exposed to chronic social defeat. The study was conducted between June 2016 and November 2018. Main Outcomes and Measures: Diagnoses of a relatively broad diagnostic spectrum of anxiety and stress-related disorders. Results: The study sample included 12 655 individuals with various anxiety and stress-related diagnoses and 19 225 controls. Overall, 17 740 study participants (55.6%) were women. A total of 7308 participants (22.9%) were born between 1981-1985, 8840 (27.7%) between 1986-1990, 8157 (25.6%) between 1991-1995, 5918 (18.6%) between 1996-2000, and 1657 (5.2%) between 2001-2005. Standard association analysis revealed variants in PDE4B to be associated with anxiety and stress-related disorder (rs7528604; P = 5.39 × 10-11; odds ratio = 0.89; 95% CI, 0.86-0.92). A framework of sensitivity analyses adjusting for mental comorbidity supported this result showing consistent association of PDE4B variants with anxiety and stress-related disorder across analytical scenarios. In mouse models, alterations in Pde4b expression were observed in those mice displaying anxiety-like behavior after exposure to chronic stress in the prefrontal cortex (P = .002; t = -3.33) and the hippocampus (P = .001; t = -3.72). We also found a single-nucleotide polymorphism heritability of 28% (standard error = 0.027) and that the genetic signature of anxiety and stress-related overlapped with psychiatric traits, educational outcomes, obesity-related phenotypes, smoking, and reproductive success. Conclusions and Relevance: This study highlights anxiety and stress-related disorders as complex heritable phenotypes with intriguing genetic correlations not only with psychiatric traits, but also with educational outcomes and multiple obesity-related phenotypes. Furthermore, we highlight the candidate gene PDE4B as a robust risk locus pointing to the potential of PDE4B inhibitors in treatment of these disorders.
Anxiety and stress-related disorders (excluding co-morbid schizophrenia)
Stress-related disorders
Anxiety and stress-related disorders (psychiatric co-morbidity weighting design)
Anxiety and stress-related disorders (psychiatric co-morbidity covariate model)
Anxiety disorders
Anxiety and stress-related disorders
Anxiety and stress-related disorders (including obsessive-compulsive disorder)
Anxiety and stress-related disorders (excluding adjustment disorder)
Understanding genetic architecture of plasma lipidome could provide better insights into lipid metabolism and its link to cardiovascular diseases (CVDs). Here, we perform genome-wide association analyses of 141 lipid species (n = 2,181 individuals), followed by phenome-wide scans with 25 CVD related phenotypes (n = 511,700 individuals). We identify 35 lipid-species-associated loci (P <5 ×10-8), 10 of which associate with CVD risk including five new loci-COL5A1, GLTPD2, SPTLC3, MBOAT7 and GALNT16 (false discovery rate<0.05). We identify loci for lipid species that are shown to predict CVD e.g., SPTLC3 for CER(d18:1/24:1). We show that lipoprotein lipase (LPL) may more efficiently hydrolyze medium length triacylglycerides (TAGs) than others. Polyunsaturated lipids have highest heritability and genetic correlations, suggesting considerable genetic regulation at fatty acids levels. We find low genetic correlations between traditional lipids and lipid species. Our results show that lipidomic profiles capture information beyond traditional lipids and identify genetic variants modifying lipid levels and risk of CVD.
Recently, HLA class II loci, including HLA-DPB1, have been reported to be associated with interindividual variance in the hepatitis B (HB) vaccine response. In this study, we investigated significant single nucleotide polymorphisms (SNPs) for anti-HBs antibody levels in 6867 healthy Koreans using a genome-wide association study (GWAS). In GWAS, the top 20 SNPs that showed significant association with anti-HBs levels (P < 1.0 × 10-29 ) all resided in HLA-DPB1. Utilizing PCR sequencing, we verified the relationship of the top 3 most significant SNPs (rs1042169, rs9277355 and rs9277356) from the GWAS and genotypes of HLA-DPB1 with the HB vaccine response in Korean infants who received a scheduled vaccination. The DPB1*04:02 allele has G, C and A nucleotides for the 3SNP sites, and was significantly more frequent in responders than in nonresponders (10.9% vs 1.0%, Pc = 0.018). DPB1*05:01 was significantly more frequent in nonresponders than in responders (49.0% vs 31.1%, Pc = 0.018). In multivariate logistic regression, DPB1*04:02 showed a significant association with both vaccine response (P = 0.037, OR = 8.465) and high-titre response (P = 0.027, OR = 9.860). The haplotypes rs1042169 G - rs9277355 C - rs9277356 A showed a significant association with a high-titre response only (P = 0.002, OR = 2.941). In conclusion, DPB1*04:02 possessing rs1042169 G - rs9277355 C - rs9277356 A is an independent predictor of the HB vaccine response in Koreans.
Uterine fibroids affect up to 77% of women by menopause and account for up to $34 billion in healthcare costs each year. Although fibroid risk is heritable, genetic risk for fibroids is not well understood. We conducted a two-stage case-control meta-analysis of genetic variants in European and African ancestry women with and without fibroids classified by a previously published algorithm requiring pelvic imaging or confirmed diagnosis. Women from seven electronic Medical Records and Genomics (eMERGE) network sites (3,704 imaging-confirmed cases and 5,591 imaging-confirmed controls) and women of African and European ancestry from UK Biobank (UKB, 5,772 cases and 61,457 controls) were included in the discovery genome-wide association study (GWAS) meta-analysis. Variants showing evidence of association in Stage I GWAS (P < 1 × 10-5) were targeted in an independent replication sample of African and European ancestry individuals from the UKB (Stage II) (12,358 cases and 138,477 controls). Logistic regression models were fit with genetic markers imputed to a 1000 Genomes reference and adjusted for principal components for each race- and site-specific dataset, followed by fixed-effects meta-analysis. Final analysis with 21,804 cases and 205,525 controls identified 326 genome-wide significant variants in 11 loci, with three novel loci at chromosome 1q24 (sentinel-SNP rs14361789; P = 4.7 × 10-8), chromosome 16q12.1 (sentinel-SNP rs4785384; P = 1.5 × 10-9) and chromosome 20q13.1 (sentinel-SNP rs6094982; P = 2.6 × 10-8). Our statistically significant findings further support previously reported loci including SNPs near WT1, TNRC6B, SYNE1, BET1L, and CDC42/WNT4. We report evidence of ancestry-specific findings for sentinel-SNP rs10917151 in the CDC42/WNT4 locus (P = 1.76 × 10-24). Ancestry-specific effect-estimates for rs10917151 were in opposite directions (P-Het-between-groups = 0.04) for predominantly African (OR = 0.84) and predominantly European women (OR = 1.16). Genetically-predicted gene expression of several genes including LUZP1 in vagina (P = 4.6 × 10-8), OBFC1 in esophageal mucosa (P = 8.7 × 10-8), NUDT13 in multiple tissues including subcutaneous adipose tissue (P = 3.3 × 10-6), and HEATR3 in skeletal muscle tissue (P = 5.8 × 10-6) were associated with fibroids. The finding for HEATR3 was supported by SNP-based summary Mendelian randomization analysis. Our study suggests that fibroid risk variants act through regulatory mechanisms affecting gene expression and are comprised of alleles that are both ancestry-specific and shared across continental ancestries.
INTRODUCTION: Inherited susceptibility to lung cancer risk in never-smokers is poorly understood. The major reason for this gap in knowledge is that this disease is relatively uncommon (except in Asians), making it difficult to assemble an adequate study sample. In this study we conducted a genome-wide association study on the largest, to date, set of European-descent never-smokers with lung cancer. METHODS: We conducted a two-phase (discovery and replication) genome-wide association study in never-smokers of European descent. We further augmented the sample by performing a meta-analysis with never-smokers from the recent OncoArray study, which resulted in a total of 3636 cases and 6295 controls. We also compare our findings with those in smokers with lung cancer. RESULTS: We detected three genome-wide statistically significant single nucleotide polymorphisms rs31490 (odds ratio [OR]: 0.769, 95% confidence interval [CI]: 0.722-0.820; p value 5.31 × 10-16), rs380286 (OR: 0.770, 95% CI: 0.723-0.820; p value 4.32 × 10-16), and rs4975616 (OR: 0.778, 95% CI: 0.730-0.829; p value 1.04 × 10-14). All three mapped to Chromosome 5 CLPTM1L-TERT region, previously shown to be associated with lung cancer risk in smokers and in never-smoker Asian women, and risk of other cancers including breast, ovarian, colorectal, and prostate. CONCLUSIONS: We found that genetic susceptibility to lung cancer in never-smokers is associated to genetic variants with pan-cancer risk effects. The comparison with smokers shows that top variants previously shown to be associated with lung cancer risk only confer risk in the presence of tobacco exposure, underscoring the importance of gene-environment interactions in the etiology of this disease.
Ninety per cent of the human population has been right-handed since the Paleolithic, yet the brain signature and genetic basis of handedness remain poorly characterized. Here, we correlated brain imaging phenotypes from ?9000 UK Biobank participants with handedness, and with loci found significantly associated with handedness after we performed genome-wide association studies (GWAS) in ?400 000 of these participants. Our imaging-handedness analysis revealed an increase in functional connectivity between left and right language networks in left-handers. GWAS of handedness uncovered four significant loci (rs199512, rs45608532, rs13017199, and rs3094128), three of which are in-or expression quantitative trait loci of-genes encoding proteins involved in brain development and patterning. These included microtubule-related MAP2 and MAPT, as well as WNT3 and MICB, all implicated in the pathogenesis of diseases such as Parkinson's, Alzheimer's and schizophrenia. In particular, with rs199512, we identified a common genetic influence on handedness, psychiatric phenotypes, Parkinson's disease, and the integrity of white matter tracts connecting the same language-related regions identified in the handedness-imaging analysis. This study has identified in the general population genome-wide significant loci for human handedness in, and expression quantitative trait loci of, genes associated with brain development, microtubules and patterning. We suggest that these genetic variants contribute to neurodevelopmental lateralization of brain organization, which in turn influences both the handedness phenotype and the predisposition to develop certain neurological and psychiatric diseases.