Pulmonary function decline is a major contributor to morbidity and mortality among smokers. Post bronchodilator FEV1 and FEV1/FVC ratio are considered the standard assessment of airflow obstruction. We performed a genome-wide association study (GWAS) in 9919 current and former smokers in the COPDGene study (6659 non-Hispanic Whites [NHW] and 3260 African Americans [AA]) to identify associations with spirometric measures (post-bronchodilator FEV1 and FEV1/FVC). We also conducted meta-analysis of FEV1 and FEV1/FVC GWAS in the COPDGene, ECLIPSE, and GenKOLS cohorts (total n = 13,532).

Chronic obstructive pulmonary disease (COPD) is the leading cause of respiratory mortality worldwide. Genetic risk loci provide new insights into disease pathogenesis. We performed a genome-wide association study in 35,735 cases and 222,076 controls from the UK Biobank and additional studies from the International COPD Genetics Consortium. We identified 82 loci associated with P < 5 × 10-8; 47 of these were previously described in association with either COPD or population-based measures of lung function. Of the remaining 35 new loci, 13 were associated with lung function in 79,055 individuals from the SpiroMeta consortium. Using gene expression and regulation data, we identified functional enrichment of COPD risk loci in lung tissue, smooth muscle, and several lung cell types. We found 14 COPD loci shared with either asthma or pulmonary fibrosis. COPD genetic risk loci clustered into groups based on associations with quantitative imaging features and comorbidities. Our analyses provide further support for the genetic susceptibility and heterogeneity of COPD.

A growing number of studies clearly demonstrate a substantial association between chronic obstructive pulmonary disease (COPD) and cardiovascular diseases (CVD), although little is known about the shared genetics that contribute to this association.

While whole genome prediction (WGP) methods have recently demonstrated successes in the prediction of complex genetic diseases, they have not yet been applied to asthma and related phenotypes. Longitudinal patterns of lung function differ between asthmatics, but these phenotypes have not been assessed for heritability or predictive ability. Herein, we assess the heritability and genetic predictability of asthma-related phenotypes.

Due to the pleiotropic effects of nitric oxide (NO) within the lungs, it is likely that NO is a significant factor in the pathogenesis of chronic obstructive pulmonary disease (COPD). The aim of this study was to test for association between single nucleotide polymorphisms (SNPs) in three NO synthase (NOS) genes and lung function, as well as to examine gene expression and protein levels in relation to the genetic variation.

We performed a genome-wide association study for chronic obstructive pulmonary disease (COPD) in three population cohorts, including 2,940 cases and 1,380 controls who were current or former smokers with normal lung function. We identified a new susceptibility locus at 4q22.1 in FAM13A and replicated this association in one case-control group (n = 1,006) and two family-based cohorts (n = 3,808) (rs7671167, combined P = 1.2 x 10(-11), combined odds ratio in case-control studies 0.76, 95% confidence interval 0.69-0.83).

Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) significantly associated with chronic obstructive pulmonary disease (COPD). However, many genetic variants show suggestive evidence for association but do not meet the strict threshold for genome-wide significance. Integrative analysis of multiple omics datasets has the potential to identify novel genes involved in disease pathogenesis by leveraging these variants in a functional, regulatory context.

Many well-powered genome-wide association studies have identified genetic determinants of self-reported smoking behaviors and measures of nicotine dependence, but most have not considered the role of structural variants, such as copy number variation (CNVs), influencing these phenotypes. Here, we included 2,889 African American and 6,187 non-Hispanic White subjects from the COPDGene cohort (http://www.copdgene.org) to carefully investigate the role of polymorphic CNVs across the genome on various measures of smoking behavior. We identified a CNV component (a hemizygous deletion) on chromosome 3p26.1 associated with two quantitative phenotypes related to smoking behavior among African Americans. This polymorphic hemizygous deletion is significantly associated with pack-years and cigarettes smoked per day among African American subjects in the COPDGene study. We sought evidence of replication in African Americans from the population based Atherosclerosis Risk in Communities (ARIC) study. While we observed similar CNV counts, the extent of exposure to cigarette smoking among ARIC subjects was quite different and the smaller sample size of heavy smokers in ARIC severely limited statistical power, so we were unable to replicate our findings from the COPDGene cohort. But meta-analyses of COPDGene and ARIC study subjects strengthened our association signal. However, a few linkage studies have reported suggestive linkage to the 3p26.1 region, and a few genome-wide association studies (GWAS) have reported markers in the gene (GRM7) nearest to this 3p26.1 area of polymorphic deletions are associated with measures of nicotine dependence among subjects of European ancestry.

Genetic factors influence chronic obstructive pulmonary disease (COPD) risk, but the individual variants that have been identified have small effects. We hypothesised that a polygenic risk score using additional variants would predict COPD and associated phenotypes.

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

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes)1. In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%.

Emphysema is an important feature of chronic obstructive pulmonary disease (COPD). Genetic factors likely affect emphysema pathogenesis, but this question has predominantly been studied in those of European ancestry. In this study, we sought to determine genetic components of emphysema severity and characterize the potential function of the associated loci in Korean population. We performed a genome-wide association study (GWAS) on quantitative emphysema in subjects with or without COPD from two Korean COPD cohorts. We investigated the functional consequences of the loci using epigenetic annotation and gene expression data. We also compared our GWAS results with an epigenome-wide association study and previous differential gene expression analysis. In total, 548 subjects (476 [86.9%] male) including 514 COPD patients were evaluated. We identified one genome-wide significant SNP (P < 5.0 × 10-8), rs117084279, near PIBF1. We identified an additional 57 SNPs (P < 5.0 × 10-6) associated with emphysema in all subjects, and 106 SNPs (P < 5.0 × 10-6) in COPD patients. Of these candidate SNPs, 2 (rs12459249, rs11667314) near CYP2A6 were expression quantitative trait loci in lung tissue and a SNP (rs11214944) near NNMT was an expression quantitative trait locus in whole blood. Of note, rs11214944 was in linkage disequilibrium with variants in enhancer histone marks in lung tissue. Several genes near additional SNPs were identified in our previous EWAS study with nominal level of significance. We identified a novel SNP associated with quantitative emphysema on CT. Including the novel SNP, several candidate SNPs in our study may provide clues to the genetic etiology of emphysema in Asian populations. Further research and validation of the loci will help determine the genetic factors for the development of emphysema.

Ever larger Structural Variant (SV) catalogs highlighting the diversity within and between populations help researchers better understand the links between SVs and disease. The identification of SVs from DNA sequence data is non-trivial and requires a balance between comprehensiveness and precision. Here we present a catalog of 355,667 SVs (59.34% novel) across autosomes and the X chromosome (50bp+) from 138,134 individuals in the diverse TOPMed consortium. We describe our methodologies for SV inference resulting in high variant quality and >90% allele concordance compared to long-read de-novo assemblies of well-characterized control samples. We demonstrate utility through significant associations between SVs and important various cardio-metabolic and hemotologic traits. We have identified 690 SV hotspots and deserts and those that potentially impact the regulation of medically relevant genes. This catalog characterizes SVs across multiple populations and will serve as a valuable tool to understand the impact of SV on disease development and progression.

The association between genetic variants on the X chromosome to risk of COPD has not been fully explored. We hypothesize that the X chromosome harbors variants important in determining risk of COPD related phenotypes and may drive sex differences in COPD manifestations.

Chronic obstructive pulmonary disease (COPD) is the third leading cause of death globally. COPD patients with cachexia or weight loss have increased risk of death independent of body mass index (BMI) and lung function. We tested the hypothesis genetic variation is associated with weight loss in COPD using a genome-wide association study approach.

The serpin family A member 1 (SERPINA1) Z allele is present in approximately one in 25 individuals of European ancestry. Z allele homozygosity (Pi*ZZ) is the most common cause of alpha 1-antitrypsin deficiency and is a proven risk factor for cirrhosis. We examined whether heterozygous Z allele (Pi*Z) carriers in United Kingdom (UK) Biobank, a population-based cohort, are at increased risk of liver disease. We replicated findings in Massachusetts General Brigham Biobank, a hospital-based cohort. We also examined variants associated with liver disease and assessed for gene-gene and gene-environment interactions. In UK Biobank, we identified 1,493 cases of cirrhosis, 12,603 Z allele heterozygotes, and 129 Z allele homozygotes among 312,671 unrelated white British participants. Heterozygous carriage of the Z allele was associated with cirrhosis compared to noncarriage (odds ratio [OR], 1.53; P = 1.1×10-04); homozygosity of the Z allele also increased the risk of cirrhosis (OR, 11.8; P = 1.8 × 10-09). The OR for cirrhosis of the Z allele was comparable to that of well-established genetic variants, including patatin-like phospholipase domain containing 3 (PNPLA3) I148M (OR, 1.48; P = 1.1 × 10-22) and transmembrane 6 superfamily member 2 (TM6SF2) E167K (OR, 1.34; P = 2.6 × 10-06). In heterozygotes compared to noncarriers, the Z allele was associated with higher alanine aminotransferase (ALT; P = = 4.6 × 10-46), aspartate aminotransferase (AST; P = 2.2 × 10-27), alkaline phosphatase (P = 3.3 × 10-43), gamma-glutamyltransferase (P = 1.2 × 10-05), and total bilirubin (P = 6.4 × 10-06); Z allele homozygotes had even greater elevations in liver biochemistries. Body mass index (BMI) amplified the association of the Z allele for ALT (P interaction = 0.021) and AST (P interaction = 0.0040), suggesting a gene-environment interaction. Finally, we demonstrated genetic interactions between variants in PNPLA3, TM6SF2, and hydroxysteroid 17-beta dehydrogenase 13 (HSD17B13); there was no evidence of epistasis between the Z allele and these variants. Conclusion: SERPINA1 Z allele heterozygosity is an important risk factor for liver disease; this risk is amplified by increasing BMI.

The analysis of whole-genome sequencing studies is challenging due to the large number of rare variants in noncoding regions and the lack of natural units for testing. We propose a statistical method to detect and localize rare and common risk variants in whole-genome sequencing studies based on a recently developed knockoff framework. It can (1) prioritize causal variants over associations due to linkage disequilibrium thereby improving interpretability; (2) help distinguish the signal due to rare variants from shadow effects of significant common variants nearby; (3) integrate multiple knockoffs for improved power, stability, and reproducibility; and (4) flexibly incorporate state-of-the-art and future association tests to achieve the benefits proposed here. In applications to whole-genome sequencing data from the Alzheimer's Disease Sequencing Project (ADSP) and COPDGene samples from NHLBI Trans-Omics for Precision Medicine (TOPMed) Program we show that our method compared with conventional association tests can lead to substantially more discoveries.

Chronic obstructive pulmonary disease (COPD) is a highly morbid and heterogenous disease. While COPD is defined by spirometry, many COPD characteristics are seen in cigarette smokers with normal spirometry. The extent to which COPD and COPD heterogeneity is captured in omics of lung tissue is not known.

The determinants and health outcomes of lung function trajectories in adults among the general population are poorly understood. We aimed to identify and characterise clusters of lung function trajectories in adults aged ≥45 years.

Although chronic low-grade inflammation does not cause immediate clinical symptoms, over the longer term, it can enhance other insults or age-dependent damage to organ systems and thereby contribute to age-related disorders, such as respiratory disorders, heart disease, metabolic disorders, autoimmunity, and cancer. However, the molecular mechanisms governing low-level inflammation are largely unknown. We discovered that Bcl-2-interacting killer (Bik) deficiency causes low-level inflammation even at baseline and the development of spontaneous emphysema in female but not male mice. Similarly, a single nucleotide polymorphism that reduced Bik levels was associated with increased inflammation and enhanced decline in lung function in humans. Transgenic expression of Bik in the airways of Bik-deficient mice inhibited allergen- or LPS-induced lung inflammation and reversed emphysema in female mice. Bik deficiency increased nuclear but not cytosolic p65 levels because Bik, by modifying the BH4 domain of Bcl-2, interacted with regulatory particle non-ATPase 1 (RPN1) and RPN2 and enhanced proteasomal degradation of nuclear proteins. Bik deficiency increased inflammation primarily in females because Bcl-2 and Bik levels were reduced in lung tissues and airway cells of female compared with male mice. Therefore, controlling low-grade inflammation by modifying the unappreciated role of Bik and Bcl-2 in facilitating proteasomal degradation of nuclear proteins may be crucial in treating chronic age-related diseases.