DeepWAS: Multivariate genotype-phenotype associations by directly integrating regulatory information using deep learning.

Genome-wide association studies (GWAS) identify genetic variants associated with traits or diseases. GWAS never directly link variants to regulatory mechanisms. Instead, the functional annotation of variants is typically inferred by post hoc analyses. A specific class of deep learning-based methods allows for the prediction of regulatory effects per variant on several cell type-specific chromatin features. We here describe "DeepWAS", a new approach that integrates these regulatory effect predictions of single variants into a multivariate GWAS setting. Thereby, single variants associated with a trait or disease are directly coupled to their impact on a chromatin feature in a cell type. Up to 61 regulatory SNPs, called dSNPs, were associated with multiple sclerosis (MS, 4,888 cases and 10,395 controls), major depressive disorder (MDD, 1,475 cases and 2,144 controls), and height (5,974 individuals). These variants were mainly non-coding and reached at least nominal significance in classical GWAS. The prediction accuracy was higher for DeepWAS than for classical GWAS models for 91% of the genome-wide significant, MS-specific dSNPs. DSNPs were enriched in public or cohort-matched expression and methylation quantitative trait loci and we demonstrated the potential of DeepWAS to generate testable functional hypotheses based on genotype data alone. DeepWAS is available at https://github.com/cellmapslab/DeepWAS.

KLB is associated with alcohol drinking, and its gene product ß-Klotho is necessary for FGF21 regulation of alcohol preference.

Excessive alcohol consumption is a major public health problem worldwide. Although drinking habits are known to be inherited, few genes have been identified that are robustly linked to alcohol drinking. We conducted a genome-wide association metaanalysis and replication study among >105,000 individuals of European ancestry and identified ß-Klotho (KLB) as a locus associated with alcohol consumption (rs11940694; P = 9.2 × 10-12). ß-Klotho is an obligate coreceptor for the hormone FGF21, which is secreted from the liver and implicated in macronutrient preference in humans. We show that brain-specific ß-Klotho KO mice have an increased alcohol preference and that FGF21 inhibits alcohol drinking by acting on the brain. These data suggest that a liver-brain endocrine axis may play an important role in the regulation of alcohol drinking behavior and provide a unique pharmacologic target for reducing alcohol consumption.

Integrative pathway genomics of lung function and airflow obstruction.

Chronic respiratory disorders are important contributors to the global burden of disease. Genome-wide association studies (GWASs) of lung function measures have identified several trait-associated loci, but explain only a modest portion of the phenotypic variability. We postulated that integrating pathway-based methods with GWASs of pulmonary function and airflow obstruction would identify a broader repertoire of genes and processes influencing these traits. We performed two independent GWASs of lung function and applied gene set enrichment analysis to one of the studies and validated the results using the second GWAS. We identified 131 significantly enriched gene sets associated with lung function and clustered them into larger biological modules involved in diverse processes including development, immunity, cell signaling, proliferation and arachidonic acid. We found that enrichment of gene sets was not driven by GWAS-significant variants or loci, but instead by those with less stringent association P-values. Next, we applied pathway enrichment analysis to a meta-analyzed GWAS of airflow obstruction. We identified several biologic modules that functionally overlapped with those associated with pulmonary function. However, differences were also noted, including enrichment of extracellular matrix (ECM) processes specifically in the airflow obstruction study. Network analysis of the ECM module implicated a candidate gene, matrix metalloproteinase 10 (MMP10), as a putative disease target. We used a knockout mouse model to functionally validate MMP10's role in influencing lung's susceptibility to cigarette smoke-induced emphysema. By integrating pathway analysis with population-based genomics, we unraveled biologic processes underlying pulmonary function traits and identified a candidate gene for obstructive lung disease.

Identification of novel genetic Loci associated with thyroid peroxidase antibodies and clinical thyroid disease.

Autoimmune thyroid diseases (AITD) are common, affecting 2-5% of the general population. Individuals with positive thyroid peroxidase antibodies (TPOAbs) have an increased risk of autoimmune hypothyroidism (Hashimoto's thyroiditis), as well as autoimmune hyperthyroidism (Graves' disease). As the possible causative genes of TPOAbs and AITD remain largely unknown, we performed GWAS meta-analyses in 18,297 individuals for TPOAb-positivity (1769 TPOAb-positives and 16,528 TPOAb-negatives) and in 12,353 individuals for TPOAb serum levels, with replication in 8,990 individuals. Significant associations (P<5×10(-8)) were detected at TPO-rs11675434, ATXN2-rs653178, and BACH2-rs10944479 for TPOAb-positivity, and at TPO-rs11675434, MAGI3-rs1230666, and KALRN-rs2010099 for TPOAb levels. Individual and combined effects (genetic risk scores) of these variants on (subclinical) hypo- and hyperthyroidism, goiter and thyroid cancer were studied. Individuals with a high genetic risk score had, besides an increased risk of TPOAb-positivity (OR: 2.18, 95% CI 1.68-2.81, P = 8.1×10(-8)), a higher risk of increased thyroid-stimulating hormone levels (OR: 1.51, 95% CI 1.26-1.82, P = 2.9×10(-6)), as well as a decreased risk of goiter (OR: 0.77, 95% CI 0.66-0.89, P = 6.5×10(-4)). The MAGI3 and BACH2 variants were associated with an increased risk of hyperthyroidism, which was replicated in an independent cohort of patients with Graves' disease (OR: 1.37, 95% CI 1.22-1.54, P = 1.2×10(-7) and OR: 1.25, 95% CI 1.12-1.39, P = 6.2×10(-5)). The MAGI3 variant was also associated with an increased risk of hypothyroidism (OR: 1.57, 95% CI 1.18-2.10, P = 1.9×10(-3)). This first GWAS meta-analysis for TPOAbs identified five newly associated loci, three of which were also associated with clinical thyroid disease. With these markers we identified a large subgroup in the general population with a substantially increased risk of TPOAbs. The results provide insight into why individuals with thyroid autoimmunity do or do not eventually develop thyroid disease, and these markers may therefore predict which TPOAb-positives are particularly at risk of developing clinical thyroid dysfunction.

Splenic Vein Thrombosis as a Rare Complication of Disseminated Tuberculosis - Imaging Diagnosis and Case Report.

BACKGROUND: Venous thromboembolism is a known but rare complication associated with Mycobacterium tuberculosis infection. The reported incidence of venous thromboembolism is 1.5-3.4% of infected patients, and it occurs due to a hypercoagulable state induced by the associated inflammation. CASE REPORT: A young woman with pulmonary tuberculosis was found to have disseminated tuberculosis and a clinically unsuspected partial thrombus in the splenic vein on imaging. Ultrasound demonstrated hepato-splenomegaly with multiple granulomas as well as ascites and a left-sided pleural effusion. An increased calibre of the splenic vein with a hyperechogenicity within it raised the suspicion of a thrombus, which was confirmed on a contrast-enhanced CT examination. CECT of the abdomen also showed a small peripheral splenic infarct, while CECT of the chest revealed bilateral miliary lesions in the lungs along with necrotic mediastinal lymphadenopathy. The final imaging diagnosis was disseminated tuberculosis complicated by splenic vein thrombosis. A timely institution of anti-coagulant and anti-tubercular treatment led to a complete resolution of the splenic vein thrombosis. CONCLUSIONS: Contrast-enhanced CT serves as a useful imaging tool for the detection of venous thrombosis and for the estimation of a complete burden of the disease. This condition should be kept in mind by both clinicians and radiologists and looked for in order to prevent life-threatening complications.

Genome-wide meta-analysis identifies six novel loci associated with habitual coffee consumption.

Coffee, a major dietary source of caffeine, is among the most widely consumed beverages in the world and has received considerable attention regarding health risks and benefits. We conducted a genome-wide (GW) meta-analysis of predominately regular-type coffee consumption (cups per day) among up to 91,462 coffee consumers of European ancestry with top single-nucleotide polymorphisms (SNPs) followed-up in ~30 062 and 7964 coffee consumers of European and African-American ancestry, respectively. Studies from both stages were combined in a trans-ethnic meta-analysis. Confirmed loci were examined for putative functional and biological relevance. Eight loci, including six novel loci, met GW significance (log10Bayes factor (BF)>5.64) with per-allele effect sizes of 0.03-0.14 cups per day. Six are located in or near genes potentially involved in pharmacokinetics (ABCG2, AHR, POR and CYP1A2) and pharmacodynamics (BDNF and SLC6A4) of caffeine. Two map to GCKR and MLXIPL genes related to metabolic traits but lacking known roles in coffee consumption. Enhancer and promoter histone marks populate the regions of many confirmed loci and several potential regulatory SNPs are highly correlated with the lead SNP of each. SNP alleles near GCKR, MLXIPL, BDNF and CYP1A2 that were associated with higher coffee consumption have previously been associated with smoking initiation, higher adiposity and fasting insulin and glucose but lower blood pressure and favorable lipid, inflammatory and liver enzyme profiles (P<5 × 10(-8)).Our genetic findings among European and African-American adults reinforce the role of caffeine in mediating habitual coffee consumption and may point to molecular mechanisms underlying inter-individual variability in pharmacological and health effects of coffee.

Meta-analysis of two genome-wide association studies identifies four genetic loci associated with thyroid function.

Thyroid hormones play key roles in cellular growth, development and metabolism. Although there is a strong genetic influence on thyroid hormone levels, the genes involved are widely unknown. The levels of circulating thyroid hormones are tightly regulated by thyrotropin (TSH), which also represents the most important diagnostic marker for thyroid function. Therefore, in order to identify genetic loci associated with TSH levels, we performed a discovery meta-analysis of two genome-wide association studies including two cohorts from Germany, KORA (n = 1287) and SHIP (n = 2449), resulting in a total sample size of 3736. Four genetic loci at 5q13.3, 1p36, 16q23 and 4q31 were associated with serum TSH levels. The lead single-nucleotide polymorphisms of these four loci were located within PDE8B encoding phosphodiesterase 8B, upstream of CAPZB that encodes the ß-subunit of the barbed-end F-actin-binding protein, in a former 'gene desert' that was recently demonstrated to encode a functional gene (LOC440389) associated with thyroid volume, and upstream of NR3C2 encoding the mineralocorticoid receptor. The latter association for the first time suggests the modulation of thyroid function by mineral corticoids. All four loci were replicated in three additional cohorts: the HUNT study from Norway (n = 1487) and the two German studies CARLA (CARLA, n = 1357) and SHIP-TREND (n = 883). Together, these four quantitative trait loci accounted for ∼3.3% of the variance in TSH serum levels. These results contribute to our understanding of genetic factors and physiological mechanisms mediating thyroid function.

Genome-wide association study identifies four genetic loci associated with thyroid volume and goiter risk.

Thyroid disorders such as goiters represent important diseases, especially in iodine-deficient areas. Sibling studies have demonstrated that genetic factors substantially contribute to the interindividual variation of thyroid volume. We performed a genome-wide association study of this phenotype by analyzing a discovery cohort consisting of 3620 participants of the Study of Health in Pomerania (SHIP). Four genetic loci were associated with thyroid volume on a genome-wide level of significance. Of these, two independent loci are located upstream of and within CAPZB, which encodes the ß subunit of the barbed-end F-actin binding protein that modulates actin polymerization, a process crucial in the colloid engulfment during thyroglobulin mobilization in the thyroid. The third locus marks FGF7, which encodes fibroblast growth factor 7. Members of this protein family have been discussed as putative signal molecules involved in the regulation of thyroid development. The fourth locus represents a "gene desert" on chromosome 16q23, located directly downstream of the predicted coding sequence LOC440389, which, however, had already been removed from the NCBI database as a result of the standard genome annotation processing at the time that this study was initiated. Experimental proof of the formerly predicted mature mRNA, however, demonstrates that LOC440389 indeed represents a real gene. All four associations were replicated in an independent sample of 1290 participants of the KORA study. These results increase the knowledge about genetic factors and physiological mechanisms influencing thyroid volume.

Luxembourg Parkinson's study -comprehensive baseline analysis of Parkinson's disease and atypical parkinsonism.

Background: Deep phenotyping of Parkinson's disease (PD) is essential to investigate this fastest-growing neurodegenerative disorder. Since 2015, over 800 individuals with PD and atypical parkinsonism along with more than 800 control subjects have been recruited in the frame of the observational, monocentric, nation-wide, longitudinal-prospective Luxembourg Parkinson's study. Objective: To profile the baseline dataset and to explore risk factors, comorbidities and clinical profiles associated with PD, atypical parkinsonism and controls. Methods: Epidemiological and clinical characteristics of all 1,648 participants divided in disease and control groups were investigated. Then, a cross-sectional group comparison was performed between the three largest groups: PD, progressive supranuclear palsy (PSP) and controls. Subsequently, multiple linear and logistic regression models were fitted adjusting for confounders. Results: The mean (SD) age at onset (AAO) of PD was 62.3 (11.8) years with 15% early onset (AAO < 50 years), mean disease duration 4.90 (5.16) years, male sex 66.5% and mean MDS-UPDRS III 35.2 (16.3). For PSP, the respective values were: 67.6 (8.2) years, all PSP with AAO > 50 years, 2.80 (2.62) years, 62.7% and 53.3 (19.5). The highest frequency of hyposmia was detected in PD followed by PSP and controls (72.9%; 53.2%; 14.7%), challenging the use of hyposmia as discriminating feature in PD vs. PSP. Alcohol abstinence was significantly higher in PD than controls (17.6 vs. 12.9%, p = 0.003). Conclusion: Luxembourg Parkinson's study constitutes a valuable resource to strengthen the understanding of complex traits in the aforementioned neurodegenerative disorders. It corroborated several previously observed clinical profiles, and provided insight on frequency of hyposmia in PSP and dietary habits, such as alcohol abstinence in PD.Clinical trial registration: clinicaltrials.gov, NCT05266872.

Immune Cell Infiltration into the Brain After Ischemic Stroke in Humans Compared to Mice and Rats: a Systematic Review and Meta-Analysis.

Although several studies have suggested that anti-inflammatory strategies reduce secondary infarct growth in animal stroke models, clinical studies have not yet demonstrated a clear benefit of immune modulation in patients. Potential reasons include systematic differences of post-ischemic neuroinflammation between humans and rodents. We here performed a systematic review and meta-analysis to summarize and compare the spatial and temporal distribution of immune cell infiltration in human and rodent stroke. Data on spatiotemporal distribution of immune cells (T cells, macrophages, and neutrophils) and infarct volume were extracted. Data from all rodent studies were pooled by means of a random-effect meta-analysis. Overall, 20 human and 188 rodent stroke studies were included in our analyses. In both patients and rodents, the infiltration of macrophages and neutrophils preceded the lymphocytic influx. Macrophages and neutrophils were the predominant immune cells within 72 h after infarction. Although highly heterogeneously across studies, the temporal profile of the poststroke immune response was comparable between patients and rodents. In rodent stroke, the extent of the immune cell infiltration depended on the duration and location of vessel occlusion and on the species. The density of infiltrating immune cells correlated with the infarct volume. In summary, we provide the first systematic analysis and comparison of human and rodent post-ischemic neuroinflammation. Our data suggest that the inflammatory response in rodent stroke models is comparable to that in patients with stroke. However, the overall heterogeneity of the post-ischemic immune response might contribute to the translational failure in stroke research.

Primary Angiitis of the CNS: A Systematic Review and Meta-analysis.

BACKGROUND AND OBJECTIVES: To facilitate and improve the diagnostic and therapeutic process by systematically reviewing studies on patients with primary angiitis of the CNS (PACNS). METHODS: We searched PubMed, looking at the period between 1988 and February 2020. Studies with adult patients with PACNS were included. We extracted and pooled proportions using fixed-effects models. Main outcomes were proportions of patients with certain clinical, imaging, and laboratory characteristics and neurologic outcomes. RESULTS: We identified 46 cohort studies including a total of 911 patients (41% biopsy confirmed, 43% angiogram confirmed, and 16% without clear assignment to the diagnostic procedure). The most frequent onset symptoms were focal neurologic signs (63%), headache (51%), and cognitive impairment (41%). Biopsy- compared with angiogram-confirmed cases had higher occurrences of cognitive impairment (55% vs 39%) and seizures (36% vs 16%), whereas focal neurologic signs occurred less often (56% vs 95%). CSF abnormalities were present in 75% vs 65% and MRI abnormalities in 97% vs 98% of patients. Digital subtraction angiography was positive in 33% of biopsy confirmed, and biopsy was positive in 8% of angiogram-confirmed cases. In 2 large cohorts, mortality was 23% and 8%, and the relapse rate was 30% and 34%, during a median follow-up of 19 and 57 months, respectively. There are no randomized trials on the treatment of PACNS. The initial treatment usually includes glucocorticoids and cyclophosphamide. DISCUSSION: PACNS is associated with disabling symptoms, frequent relapses, and significant mortality. Differences in symptoms and neuroimaging results and low overlap between biopsy and angiogram suggest that biopsy- and angiogram-confirmed cases represent different histopathologic types of PACNS. The optimal treatment is unknown.

Tuberculosis, a great masquerader: A case series unveiling rare sites of musculoskeletal involvement through imaging.

Skeletal tubercular infections that do not involve the spine or large joints are rarely encountered. This case series aims to highlight the importance of imaging in diagnosing skeletal tuberculosis (TB) at uncommon sites in clinically unsuspected patients by demonstrating specific imaging findings. We present the clinical details and imaging findings of seven pathologically confirmed cases of extraspinal skeletal TB. A multimodality imaging approach including radiography, ultrasonography (USG) and computed tomography (CT) scan was used in most cases. The imaging studies revealed an infective soft tissue collection over different sites including the sternoclavicular joint, acromion process, chest wall and temporo-mandibular joint, along with destruction and erosion of the underlying or adjacent bones. In tubercular endemic countries, strong clinical suspicion should be entertained in cases presenting with a soft tissue collection, even around unusual skeletal sites.

Age-specific incidence of breast cancer in breast cancer survivors and their first-degree relatives.

BACKGROUND: The incidence rate of breast cancer in first-degree relatives of women with breast cancer has been hypothesized to become constant at a predetermined age in accordance with observations of a high, roughly constant incidence rate of contralateral breast cancer by age. We attempted to test this hypothesis in the Danish population with cancer registry data. METHODS: We determined the age-specific incidence rates of contralateral breast cancer in Danish women who had a first breast cancer before they were 50 years of age and the rates of breast cancer among their first-degree female relatives during 1943 to 1999. The observed rates were tested for trends chi test or evaluated in Cox proportional hazard models. RESULTS: A high incidence rate of contralateral breast cancer was observed in women aged 25-44 years, followed by a decreasing rate, which reached a level corresponding to the rate per breast in the general female population at age 65. At ages older than the index patients age at diagnosis, their first-degree female relatives showed significantly increasing incidence rates of breast cancer by age, with a relatively constant absolute difference of 187 breast cancers per 100,000 person-years between the observed rates and the expected rates. CONCLUSION: The rate of contralateral breast cancer is particular high at young ages but the excess ebbs as the cohort ages, perhaps due to elimination of predisposed individuals at early ages from the cohort of survivors. First-degree relatives seem to share breast cancer susceptibility genes with the family proband resulting in a constant excess rate of breast cancer throughout life.

Trends in hormone therapy and breast cancer incidence - results from the German Network of Cancer Registries.

OBJECTIVE: Menopausal hormone therapy (HT) is a proven risk factor for breast cancer. Recent reports presented declining trends in breast cancer incidence, which were attributed to parallel declining trends in HT. Therefore, we analyzed recent data on hormone therapy and breast cancer incidence in Germany. METHODS: We performed a population-based survey using breast cancer incidence (from cancer registries) and hormone prescription data (from health insurances) for the time period from 1997 to 2006. Age-standardized rates were calculated and joinpoint regression analyses for trends were performed. RESULTS: Prescription of HT started to decline in 1999; about 2 years later, a parallel decline in breast cancer incidence was observed. HT prescription decreased by 69% up to 2006, and breast cancer incidence by 6.8% for all age groups (2002-2005) and 12.8% in the age group from 50 to 69 years. The reductions in HT prescription and breast cancer incidence were markedly correlated across the federal states in Germany. CONCLUSION: The recent decline in breast cancer incidence following decreasing HT prescription suggests a causal relationship between the risk of breast cancer and HT in Germany. Even after the 'big HT drop' in Germany, significant differences in HT prescription as well as breast cancer incidence persist at the federal state level. These results should be discussed in a public health context.

Genome-wide meta-analysis of depression identifies 102 independent variants and highlights the importance of the prefrontal brain regions.

Major depression is a debilitating psychiatric illness that is typically associated with low mood and anhedonia. Depression has a heritable component that has remained difficult to elucidate with current sample sizes due to the polygenic nature of the disorder. To maximize sample size, we meta-analyzed data on 807,553 individuals (246,363 cases and 561,190 controls) from the three largest genome-wide association studies of depression. We identified 102 independent variants, 269 genes, and 15 genesets associated with depression, including both genes and gene pathways associated with synaptic structure and neurotransmission. An enrichment analysis provided further evidence of the importance of prefrontal brain regions. In an independent replication sample of 1,306,354 individuals (414,055 cases and 892,299 controls), 87 of the 102 associated variants were significant after multiple testing correction. These findings advance our understanding of the complex genetic architecture of depression and provide several future avenues for understanding etiology and developing new treatment approaches.

Lack of Association Between Shape and Volume of Subcortical Brain Structures and Restless Legs Syndrome.

OBJECTIVE: Previous studies on patients with restless legs syndrome (RLS) yielded inconclusive results in the magnetic resonance imaging (MRI)-based analyses of alterations of subcortical structures in the brain. The aim of this study was to compare volumes as well as shapes of subcortical structures and the hippocampus between RLS cases and controls. Additionally, the associations between the genetic risks for RLS and subcortical volumes were investigated. METHODS: We compared volumetric as well as shape differences assessed by 3 T MRI in the caudate nucleus, hippocampus, globus pallidus, putamen, and thalamus in 39 RLS cases versus 117 controls, nested within a population-based sample. In a subsample, we explored associations between known genetic risk markers for RLS and the volumes of the subcortical structures and the hippocampus. RESULTS: No significant differences between RLS cases and controls in subcortical and hippocampal shapes and volumes were observed. Furthermore, the genetic risk for RLS was unrelated to any alterations of subcortical and hippocampal gray matter volume. INTERPRETATION: We conclude that neither RLS nor the genetic risk for the disease give rise to changes in hippocampal and subcortical shapes and gray matter volumes.

DNA methylation as a mediator of HLA-DRB1*15:01 and a protective variant in multiple sclerosis.

The human leukocyte antigen (HLA) haplotype DRB1*15:01 is the major risk factor for multiple sclerosis (MS). Here, we find that DRB1*15:01 is hypomethylated and predominantly expressed in monocytes among carriers of DRB1*15:01. A differentially methylated region (DMR) encompassing HLA-DRB1 exon 2 is particularly affected and displays methylation-sensitive regulatory properties in vitro. Causal inference and Mendelian randomization provide evidence that HLA variants mediate risk for MS via changes in the HLA-DRB1 DMR that modify HLA-DRB1 expression. Meta-analysis of 14,259 cases and 171,347 controls confirms that these variants confer risk from DRB1*15:01 and also identifies a protective variant (rs9267649, p < 3.32 × 10-8, odds ratio = 0.86) after conditioning for all MS-associated variants in the region. rs9267649 is associated with increased DNA methylation at the HLA-DRB1 DMR and reduced expression of HLA-DRB1, suggesting a modulation of the DRB1*15:01 effect. Our integrative approach provides insights into the molecular mechanisms of MS susceptibility and suggests putative therapeutic strategies targeting a methylation-mediated regulation of the major risk gene.

Meta-analysis of exome array data identifies six novel genetic loci for lung function.

Background: Over 90 regions of the genome have been associated with lung function to date, many of which have also been implicated in chronic obstructive pulmonary disease. Methods: We carried out meta-analyses of exome array data and three lung function measures: forced expiratory volume in one second (FEV 1), forced vital capacity (FVC) and the ratio of FEV 1 to FVC (FEV 1/FVC). These analyses by the SpiroMeta and CHARGE consortia included 60,749 individuals of European ancestry from 23 studies, and 7,721 individuals of African Ancestry from 5 studies in the discovery stage, with follow-up in up to 111,556 independent individuals. Results: We identified significant (P<2·8x10 -7) associations with six SNPs: a nonsynonymous variant in RPAP1, which is predicted to be damaging, three intronic SNPs ( SEC24C, CASC17 and UQCC1) and two intergenic SNPs near to LY86 and FGF10. Expression quantitative trait loci analyses found evidence for regulation of gene expression at three signals and implicated several genes, including TYRO3 and PLAU. Conclusions: Further interrogation of these loci could provide greater understanding of the determinants of lung function and pulmonary disease.

Genome-wide association analyses for lung function and chronic obstructive pulmonary disease identify new loci and potential druggable targets.

Chronic obstructive pulmonary disease (COPD) is characterized by reduced lung function and is the third leading cause of death globally. Through genome-wide association discovery in 48,943 individuals, selected from extremes of the lung function distribution in UK Biobank, and follow-up in 95,375 individuals, we increased the yield of independent signals for lung function from 54 to 97. A genetic risk score was associated with COPD susceptibility (odds ratio per 1 s.d. of the risk score (∼6 alleles) (95% confidence interval) = 1.24 (1.20-1.27), P = 5.05 × 10-49), and we observed a 3.7-fold difference in COPD risk between individuals in the highest and lowest genetic risk score deciles in UK Biobank. The 97 signals show enrichment in genes for development, elastic fibers and epigenetic regulation pathways. We highlight targets for drugs and compounds in development for COPD and asthma (genes in the inositol phosphate metabolism pathway and CHRM3) and describe targets for potential drug repositioning from other clinical indications.

Evidence for large-scale gene-by-smoking interaction effects on pulmonary function.

Background: Smoking is the strongest environmental risk factor for reduced pulmonary function. The genetic component of various pulmonary traits has also been demonstrated, and at least 26 loci have been reproducibly associated with either FEV 1 (forced expiratory volume in 1 second) or FEV 1 /FVC (FEV 1 /forced vital capacity). Although the main effects of smoking and genetic loci are well established, the question of potential gene-by-smoking interaction effect remains unanswered. The aim of the present study was to assess, using a genetic risk score approach, whether the effect of these 26 loci on pulmonary function is influenced by smoking. Methods: We evaluated the interaction between smoking exposure, considered as either ever vs never or pack-years, and a 26-single nucleotide polymorphisms (SNPs) genetic risk score in relation to FEV 1 or FEV 1 /FVC in 50 047 participants of European ancestry from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) and SpiroMeta consortia. Results: We identified an interaction ( ßint = -0.036, 95% confidence interval, -0.040 to -0.032, P = 0.00057) between an unweighted 26 SNP genetic risk score and smoking status (ever/never) on the FEV 1 /FVC ratio. In interpreting this interaction, we showed that the genetic risk of falling below the FEV /FVC threshold used to diagnose chronic obstructive pulmonary disease is higher among ever smokers than among never smokers. A replication analysis in two independent datasets, although not statistically significant, showed a similar trend in the interaction effect. Conclusions: This study highlights the benefit of using genetic risk scores for identifying interactions missed when studying individual SNPs and shows, for the first time, that persons with the highest genetic risk for low FEV 1 /FVC may be more susceptible to the deleterious effects of smoking.