A novel application of data-consistent inversion to overcome spurious inference in genome-wide association studies

The genome-wide association studies (GWAS) typically use linear or logistic regression models to identify associations between phenotypes (traits) and genotypes (genetic variants) of interest. However, the use of regression with the additive assumption has potential limitations. First, the normality assumption of residuals is the one that is rarely seen in practice, and deviation from normality increases the Type-I error rate. Second, building a model based on such an assumption ignores genetic structures, like, dominant, recessive, and protective-risk cases. Ignoring genetic variants may result in spurious conclusions about the associations between a variant and a trait. We propose an assumption-free model built upon data-consistent inversion (DCI), which is a recently developed measure-theoretic framework utilized for uncertainty quantification. This proposed DCI-derived model builds a nonparametric distribution on model inputs that propagates to the distribution of observed data without the required normality assumption of residuals in the regression model. This characteristic enables the proposed DCI-derived model to cover all genetic variants without emphasizing on additivity of the classic-GWAS model. Simulations and a replication GWAS with data from the COPDGene demonstrate the ability of this model to control the Type-I error rate at least as well as the classic-GWAS (additive linear model) approach while having similar or greater power to discover variants in different genetic modes of transmission.

Deep Learning Integration of Chest Computed Tomography Imaging and Gene Expression Identifies Novel Aspects of COPD

Rationale

Chronic obstructive pulmonary disease (COPD) is characterized by pathologic changes in the airways, lung parenchyma, and persistent inflammation, but the links between lung structural changes and blood transcriptome patterns have not been fully described.

Objections

The objective of this study was to identify novel relationships between lung structural changes measured by chest computed tomography (CT) and blood transcriptome patterns measured by blood RNA sequencing (RNA-seq).

Methods

CT scan images and blood RNA-seq gene expression from 1223 participants in the COPD Genetic Epidemiology (COPDGene®) study were jointly analyzed using deep learning to identify shared aspects of inflammation and lung structural changes that we labeled image-expression axes (IEAs). We related IEAs to COPD-related measurements and prospective health outcomes through regression and Cox proportional hazards models and tested them for biological pathway enrichment.

Results

We identified 2 distinct IEAs: IEAemph which captures an emphysema-predominant process with a strong positive correlation to CT emphysema and a negative correlation to forced expiratory volume in 1 second and body mass index (BMI); and IEAairway which captures an airway-predominant process with a positive correlation to BMI and airway wall thickness and a negative correlation to emphysema. Pathway enrichment analysis identified 29 and 13 pathways significantly associated with IEAemph and IEAairway, respectively (adjusted p<0.001).

Conclusions

Integration of CT scans and blood RNA-seq data identified 2 IEAs that capture distinct inflammatory processes associated with emphysema and airway-predominant COPD.

Heterogeneity and Progression of Chronic Obstructive Pulmonary Disease: Emphysema-Predominant and Non-Emphysema-Predominant Disease

While variation in emphysema severity between patients with chronic obstructive pulmonary disease (COPD) is well-recognized, clinically applicable definitions of the emphysema-predominant disease (EPD) and non-emphysema-predominant disease (NEPD) subtypes have not been established. To study the clinical relevance of the EPD and NEPD subtypes, we tested the association of these subtypes with prospective decline in forced expiratory volume in 1 second (FEV1) and mortality among 3,427 subjects with Global Initiative for Chronic Obstructive Lung Disease (GOLD) spirometric grade 2-4 COPD at baseline in the Genetic Epidemiology of COPD (COPDGene) Study, an ongoing national multicenter study that started in 2007. NEPD was defined as airflow obstruction with less than 5% computed tomography (CT) quantitative densitometric emphysema at -950 Hounsfield units, and EPD was defined as airflow obstruction with 10% or greater CT emphysema. Mixed-effects models for FEV1 demonstrated larger average annual FEV1 loss in EPD subjects than in NEPD subjects (-10.2 mL/year; P < 0.001), and subtype-specific associations with FEV1 decline were identified. Cox proportional hazards models showed higher risk of mortality among EPD patients versus NEPD patients (hazard ratio = 1.46, 95% confidence interval: 1.34, 1.60; P < 0.001). To determine whether the NEPD/EPD dichotomy is captured by previously described COPDGene subtypes, we used logistic regression and receiver operating characteristic (ROC) curve analysis to predict NEPD/EPD membership using these previous subtype definitions. The analysis generally showed excellent discrimination, with areas under the ROC curve greater than 0.9. The NEPD and EPD COPD subtypes capture important aspects of COPD heterogeneity and are associated with different rates of disease progression and mortality.

Mendelian randomization study of diabetes and dementia in the Million Veteran Program

INTRODUCTION

Diabetes and dementia are diseases of high health-care burden worldwide. Individuals with diabetes have 1.4 to 2.2 times higher risk of dementia. Our objective was to evaluate evidence of causality between these two common diseases.

METHODS

We conducted a one-sample Mendelian randomization (MR) analysis in the US Department of Veterans Affairs Million Veteran program. The study included 334,672 participants ≥65 years of age with type 2 diabetes and dementia case-control status and genotype data.

RESULTS

For each standard deviation increase in genetically predicted diabetes, we found increased odds of three dementia diagnoses in non-Hispanic White participants (all-cause: odds ratio [OR] = 1.07 [1.05-1.08], P = 3.40E-18; vascular: OR = 1.11 [1.07-1.15], P = 3.63E-09, Alzheimer's disease [AD]: OR = 1.06 [1.02-1.09], P = 6.84E-04) and non-Hispanic Black participants (all-cause: OR = 1.06 [1.02-1.10], P = 3.66E-03, vascular: OR = 1.11 [1.04-1.19], P = 2.20E-03, AD: OR = 1.12 [1.02-1.23], P = 1.60E-02) but not in Hispanic participants (all P > 0.05).

DISCUSSION

We found evidence of causality between diabetes and dementia using a one-sample MR study, with access to individual level data, overcoming limitations of prior studies using two-sample MR techniques.

Multi-ancestry genome-wide study identifies effector genes and druggable pathways for coronary artery calcification

Coronary artery calcification (CAC), a measure of subclinical atherosclerosis, predicts future symptomatic coronary artery disease (CAD). Identifying genetic risk factors for CAC may point to new therapeutic avenues for prevention. Currently, there are only four known risk loci for CAC identified from genome-wide association studies (GWAS) in the general population. Here we conducted the largest multi-ancestry GWAS meta-analysis of CAC to date, which comprised 26,909 individuals of European ancestry and 8,867 individuals of African ancestry. We identified 11 independent risk loci, of which eight were new for CAC and five had not been reported for CAD. These new CAC loci are related to bone mineralization, phosphate catabolism and hormone metabolic pathways. Several new loci harbor candidate causal genes supported by multiple lines of functional evidence and are regulators of smooth muscle cell-mediated calcification ex vivo and in vitro. Together, these findings help refine the genetic architecture of CAC and extend our understanding of the biological and potential druggable pathways underlying CAC.

Genome-wide association meta-analysis identifies 17 loci associated with nonalcoholic fatty liver disease

Nonalcoholic fatty liver disease (NAFLD) is common and partially heritable and has no effective treatments. We carried out a genome-wide association study (GWAS) meta-analysis of imaging (n = 66,814) and diagnostic code (3,584 cases versus 621,081 controls) measured NAFLD across diverse ancestries. We identified NAFLD-associated variants at torsin family 1 member B (TOR1B), fat mass and obesity associated (FTO), cordon-bleu WH2 repeat protein like 1 (COBLL1)/growth factor receptor-bound protein 14 (GRB14), insulin receptor (INSR), sterol regulatory element-binding transcription factor 1 (SREBF1) and patatin-like phospholipase domain-containing protein 2 (PNPLA2), as well as validated NAFLD-associated variants at patatin-like phospholipase domain-containing protein 3 (PNPLA3), transmembrane 6 superfamily 2 (TM6SF2), apolipoprotein E (APOE), glucokinase regulator (GCKR), tribbles homolog 1 (TRIB1), glycerol-3-phosphate acyltransferase (GPAM), mitochondrial amidoxime-reducing component 1 (MARC1), microsomal triglyceride transfer protein large subunit (MTTP), alcohol dehydrogenase 1B (ADH1B), transmembrane channel like 4 (TMC4)/membrane-bound O-acyltransferase domain containing 7 (MBOAT7) and receptor-type tyrosine-protein phosphatase δ (PTPRD). Implicated genes highlight mitochondrial, cholesterol and de novo lipogenesis as causally contributing to NAFLD predisposition. Phenome-wide association study (PheWAS) analyses suggest at least seven subtypes of NAFLD. Individuals in the top 10% and 1% of genetic risk have a 2.5-fold to 6-fold increased risk of NAFLD, cirrhosis and hepatocellular carcinoma. These genetic variants identify subtypes of NAFLD, improve estimates of disease risk and can guide the development of targeted therapeutics.

Causes of and Clinical Features Associated with Death in Tobacco Cigarette Users by Lung Function Impairment

Rationale: Cigarette smoking contributes to the risk of death through different mechanisms. Objectives: To determine how causes of and clinical features associated with death vary in tobacco cigarette users by lung function impairment. Methods: We stratified current and former tobacco cigarette users enrolled in Genetic Epidemiology of Chronic Obstructive Pulmonary Disease (COPDGene) into normal spirometry, PRISm (Preserved Ratio Impaired Spirometry), Global Initiative for Chronic Obstructive Lung Disease (GOLD) 1-2 COPD, and GOLD 3-4 COPD. Deaths were identified via longitudinal follow-up and Social Security Death Index search. Causes of death were adjudicated after a review of death certificates, medical records, and next-of-kin interviews. We tested associations between baseline clinical variables and all-cause mortality using multivariable Cox proportional hazards models. Measurements and Main Results: Over a 10.1-year median follow-up, 2,200 deaths occurred among 10,132 participants (age 59.5 ± 9.0 yr; 46.6% women). Death from cardiovascular disease was most frequent in PRISm (31% of deaths). Lung cancer deaths were most frequent in GOLD 1-2 (18% of deaths vs. 9-11% in other groups). Respiratory deaths outpaced competing causes of death in GOLD 3-4, particularly when BODE index ⩾7. St. George's Respiratory Questionnaire score ⩾25 was associated with higher mortality in all groups: Hazard ratio (HR), 1.48 (1.20-1.84) normal spirometry; HR, 1.40 (1.05-1.87) PRISm; HR, 1.80 (1.49-2.17) GOLD 1-2; HR, 1.65 (1.26-2.17) GOLD 3-4. History of respiratory exacerbations was associated with higher mortality in GOLD 1-2 and GOLD 3-4, quantitative emphysema in GOLD 1-2, and airway wall thickness in PRISm and GOLD 3-4. Conclusions: Leading causes of death vary by lung function impairment in tobacco cigarette users. Worse respiratory-related quality of life is associated with all-cause mortality regardless of lung function.

Changes in Lung Volumes with Spirometric Disease Progression in COPD

Background

Abnormal lung volumes representing air trapping identify the subset of smokers with preserved spirometry who develop spirometric chronic obstructive pulmonary disease (COPD) and adverse outcomes. However, how lung volumes evolve in early COPD as airflow obstruction develops remains unclear.

Methods

To establish how lung volumes change with the development of spirometric COPD, we examined lung volumes from the pulmonary function data (seated posture) available in the U.S. Department of Veterans Affairs electronic health records (n=71,356) and lung volumes measured by computed tomography (supine posture) available from the COPD Genetic Epidemiology (COPDGene®) study (n=7969) and the SubPopulations and InterMediate Outcome Measures In COPD Study (SPIROMICS) (n=2552) cohorts, and studied their cross-sectional distributions and longitudinal changes across the airflow obstruction spectrum. Patients with preserved ratio-impaired spirometry (PRISm) were excluded from this analysis.

Results

Lung volumes from all 3 cohorts showed similar patterns of distributions and longitudinal changes with worsening airflow obstruction. The distributions for total lung capacity (TLC), vital capacity (VC), and inspiratory capacity (IC) and their patterns of change were nonlinear and included different phases. When stratified by airflow obstruction using Global initiative for chronic Obstructive Lung Disease (GOLD) stages, patients with GOLD 1 (mild) COPD had larger lung volumes (TLC, VC, IC) compared to patients with GOLD 0 (smokers with preserved spirometry) or GOLD 2 (moderate) disease. In longitudinal follow-up of baseline GOLD 0 patients who progressed to spirometric COPD, those with an initially higher TLC and VC developed mild obstruction (GOLD 1) while those with an initially lower TLC and VC developed moderate obstruction (GOLD 2).

Conclusions

In COPD, TLC, and VC have biphasic distributions, change in nonlinear fashions as obstruction worsens, and could differentiate those GOLD 0 patients at risk for more rapid spirometric disease progression.

Longitudinal Assessment of Multimorbidity Medication Patterns among Smokers in the COPDGene Cohort

Background and objectives: Chronic obstructive pulmonary disease (COPD) is usually comorbid with other chronic diseases. We aimed to assess the multimorbidity medication patterns and explore if the patterns are similar for phase 1 (P1) and 5-year follow-up phase 2 (P2) in the COPDGene cohort. Materials and Methods: A total of 5564 out of 10,198 smokers from the COPDGene cohort who completed 2 visits, P1 and P2 visits, with complete medication use history were included in the study. We conducted latent class analysis (LCA) among the 27 categories of chronic disease medications, excluding COPD treatments and cancer medications at P1 and P2 separately. The best number of LCA classes was determined through both statistical fit and interpretation of the patterns. Results: We found four classes of medication patterns at both phases. LCA showed that both phases shared similar characteristics in their medication patterns: LC0: low medication; LC1: hypertension (HTN) or cardiovascular disease (CVD)+high cholesterol (Hychol) medication predominant; LC2: HTN/CVD+type 2 diabetes (T2D) +Hychol medication predominant; LC3: Hychol medication predominant. Conclusions: We found similar multimorbidity medication patterns among smokers at P1 and P2 in the COPDGene cohort, which provides an understanding of how multimorbidity medication clustered and how different chronic diseases combine in smokers.

Mendelian randomization study of diabetes and dementia in the Million Veteran Program

INTRODUCTION

Diabetes and dementia are diseases of high healthcare burden worldwide. Individuals with diabetes have 1.4 to 2.2 times higher risk of dementia. Our objective was to evaluate evidence of causality between these two common diseases.

METHODS

We conducted a one-sample Mendelian randomization (MR) analysis in the U.S. Department of Veterans Affairs Million Veteran program. The study included 334,672 participants ≥65 years of age with type 2 diabetes and dementia case-control status and genotype data.

RESULTS

For each standard deviation increase in genetically-predicted diabetes, we found increased odds of three dementia diagnoses in non-Hispanic White participants (all-cause: OR=1.07[1.05-1.08], P =3.40E-18; vascular: OR=1.11[1.07-1.15], P =3.63E-09, Alzheimer's: OR=1.06[1.02-1.09], P =6.84E-04) and non-Hispanic Black participants (all-cause: OR=1.06[1.02-1.10], P =3.66E-03, vascular: OR=1.11[1.04-1.19], P =2.20E-03, Alzheimer's: OR=1.12 [1.02-1.23], P =1.60E-02) but not in Hispanic participants (all P >.05).

DISCUSSION

We found evidence of causality between diabetes and dementia using a one-sample MR study, with access to individual level data, overcoming limitations of prior studies utilizing two-sample MR techniques.

Multi-ancestry transcriptome-wide association analyses yield insights into tobacco use biology and drug repurposing

Most transcriptome-wide association studies (TWASs) so far focus on European ancestry and lack diversity. To overcome this limitation, we aggregated genome-wide association study (GWAS) summary statistics, whole-genome sequences and expression quantitative trait locus (eQTL) data from diverse ancestries. We developed a new approach, TESLA (multi-ancestry integrative study using an optimal linear combination of association statistics), to integrate an eQTL dataset with a multi-ancestry GWAS. By exploiting shared phenotypic effects between ancestries and accommodating potential effect heterogeneities, TESLA improves power over other TWAS methods. When applied to tobacco use phenotypes, TESLA identified 273 new genes, up to 55% more compared with alternative TWAS methods. These hits and subsequent fine mapping using TESLA point to target genes with biological relevance. In silico drug-repurposing analyses highlight several drugs with known efficacy, including dextromethorphan and galantamine, and new drugs such as muscle relaxants that may be repurposed for treating nicotine addiction.

Genetic diversity fuels gene discovery for tobacco and alcohol use

Tobacco and alcohol use are heritable behaviours associated with 15% and 5.3% of worldwide deaths, respectively, due largely to broad increased risk for disease and injury1-4. These substances are used across the globe, yet genome-wide association studies have focused largely on individuals of European ancestries5. Here we leveraged global genetic diversity across 3.4 million individuals from four major clines of global ancestry (approximately 21% non-European) to power the discovery and fine-mapping of genomic loci associated with tobacco and alcohol use, to inform function of these loci via ancestry-aware transcriptome-wide association studies, and to evaluate the genetic architecture and predictive power of polygenic risk within and across populations. We found that increases in sample size and genetic diversity improved locus identification and fine-mapping resolution, and that a large majority of the 3,823 associated variants (from 2,143 loci) showed consistent effect sizes across ancestry dimensions. However, polygenic risk scores developed in one ancestry performed poorly in others, highlighting the continued need to increase sample sizes of diverse ancestries to realize any potential benefit of polygenic prediction.

Increased chest CT derived bone and muscle measures capture markers of improved morbidity and mortality in COPD

BACKGROUND

Chronic obstructive pulmonary disease (COPD) is a disease of accelerated aging and is associated with comorbid conditions including osteoporosis and sarcopenia. These extrapulmonary conditions are highly prevalent yet frequently underdiagnosed and overlooked by pulmonologists in COPD treatment and management. There is evidence supporting a role for bone-muscle crosstalk which may compound osteoporosis and sarcopenia risk in COPD. Chest CT is commonly utilized in COPD management, and we evaluated its utility to identify low bone mineral density (BMD) and reduced pectoralis muscle area (PMA) as surrogates for osteoporosis and sarcopenia. We then tested whether BMD and PMA were associated with morbidity and mortality in COPD.

METHODS

BMD and PMA were analyzed from chest CT scans of 8468 COPDGene participants with COPD and controls (smoking and non-smoking). Multivariable regression models tested the relationship of BMD and PMA with measures of function (6-min walk distance (6MWD), handgrip strength) and disease severity (percent emphysema and lung function). Multivariable Cox proportional hazards models were used to evaluate the relationship between sex-specific quartiles of BMD and/or PMA derived from non-smoking controls with all-cause mortality.

RESULTS

COPD subjects had significantly lower BMD and PMA compared with controls. Higher BMD and PMA were associated with increased physical function and less disease severity. Participants with the highest BMD and PMA quartiles had a significantly reduced mortality risk (36% and 46%) compared to the lowest quartiles.

CONCLUSIONS

These findings highlight the potential for CT-derived BMD and PMA to characterize osteoporosis and sarcopenia using equipment available in the pulmonary setting.

Rare genetic variants explain missing heritability in smoking

Common genetic variants explain less variation in complex phenotypes than inferred from family-based studies, and there is a debate on the source of this 'missing heritability'. We investigated the contribution of rare genetic variants to tobacco use with whole-genome sequences from up to 26,257 unrelated individuals of European ancestries and 11,743 individuals of African ancestries. Across four smoking traits, single-nucleotide-polymorphism-based heritability ([Formula: see text]) was estimated from 0.13 to 0.28 (s.e., 0.10-0.13) in European ancestries, with 35-74% of it attributable to rare variants with minor allele frequencies between 0.01% and 1%. These heritability estimates are 1.5-4 times higher than past estimates based on common variants alone and accounted for 60% to 100% of our pedigree-based estimates of narrow-sense heritability ([Formula: see text], 0.18-0.34). In the African ancestry samples, [Formula: see text] was estimated from 0.03 to 0.33 (s.e., 0.09-0.14) across the four smoking traits. These results suggest that rare variants are important contributors to the heritability of smoking.

A Diabetes Genetic Risk Score Is Associated With All-Cause Dementia and Clinically Diagnosed Vascular Dementia in the Million Veteran Program

OBJECTIVE

Diabetes and dementia are diseases of high health care burden worldwide, and studies have shown that diabetes is associated with an increased relative risk of dementia. We set out to examine whether type 2 diabetes-associated genetic variants were associated with dementia and whether they differed by race/ethnicity or clinical dementia diagnosis.

RESEARCH DESIGN AND METHODS

We evaluated associations of two type 2 diabetes genetic risk scores (GRS and GRS-nonAPOE: a score without rs429358, a variant associated with Alzheimer disease [AD]) with three classifications of clinical dementia diagnoses in the Million Veteran Program (MVP): all-cause dementia, vascular dementia (VaD), and AD. We conducted our analysis stratified by European (EUR), African (AFR), and Hispanic (HIS) races/ethnicities.

RESULTS

In EUR, we found associations of the GRS with all-cause dementia (odds ratio [OR] 1.06, P = 1.60e-07) and clinically diagnosed VaD (OR 1.12, P = 5.2e-05) but not with clinically diagnosed AD (OR 1.02, P = 0.43). The GRS was not associated with any dementia outcome in AFR or HIS. When testing with GRS-nonAPOE, we found that effect size estimates in EUR increased and P values decreased for all-cause dementia (OR 1.08, P = 2.6e-12), for VaD (OR 1.14, P = 7.2e-07), and for AD (OR 1.06, P = 0.018). For AFR, the association of GRS-nonAPOE and clinically diagnosed VaD (OR 1.15, P = 0.016) was statistically significant. There were no significant findings for HIS.

CONCLUSIONS

We found evidence suggesting shared genetic pathogenesis of diabetes with all-cause dementia and clinically diagnosed VaD.

GAWMerge expands GWAS sample size and diversity by combining array-based genotyping and whole-genome sequencing

Genome-wide association studies (GWAS) have made impactful discoveries for complex diseases, often by amassing very large sample sizes. Yet, GWAS of many diseases remain underpowered, especially for non-European ancestries. One cost-effective approach to increase sample size is to combine existing cohorts, which may have limited sample size or be case-only, with public controls, but this approach is limited by the need for a large overlap in variants across genotyping arrays and the scarcity of non-European controls. We developed and validated a protocol, Genotyping Array-WGS Merge (GAWMerge), for combining genotypes from arrays and whole-genome sequencing, ensuring complete variant overlap, and allowing for diverse samples like Trans-Omics for Precision Medicine to be used. Our protocol involves phasing, imputation, and filtering. We illustrated its ability to control technology driven artifacts and type-I error, as well as recover known disease-associated signals across technologies, independent datasets, and ancestries in smoking-related cohorts. GAWMerge enables genetic studies to leverage existing cohorts to validly increase sample size and enhance discovery for understudied traits and ancestries.

Weight gain trajectories and obesity rates in intensive and conventional treatments of type 1 diabetes from the DCCT compared with a control population without diabetes

AIM

Obesity is a significant health issue for participants with type 1 diabetes undergoing intensive diabetes management. The temporal pattern and factors associated with weight gain after treatment initiation remain poorly understood including how weight gain in participants with and without type I diabetes compare. Our aim was to compare weight gain in those receiving intensive (INT) and conventional (CONV) type 1 diabetes treatment to a population without diabetes.

METHODS

Participants included men and women of 18 years and older in the Diabetes Control and Complications Trial (DCCT) randomized to INT (n = 562) or CONV (n = 568) and a prospective, observational cohort without diabetes from the Coronary Artery Development in Young Adults (CARDIA, controls) study (n = 2446). Body mass index (BMI) trajectories and obesity prevalence were compared between groups and candidate metabolic and therapeutic moderators investigated.

RESULTS

Annual weight gain with INT peaked 1.3 years after initiation and was greater than both CONV and controls before and after this peak. Obesity prevalence with INT was lower than controls at baseline, was similar to controls at 2 years and surpassed controls by 5 years. Obesity rates with CONV remained below controls at all time points. Greater annual weight gain in the DCCT was associated with lower haemoglobin A_1c , higher insulin dose and family history of type 2 diabetes.

CONCLUSIONS

Greater weight gain accompanying INT therapy occurs in two stages, leads to similar or greater obesity rates than controls after 2 years and is primarily modified by glucose control and family history, supportive of a therapeutic-genetic influence on weight trajectories.

Significant Spirometric Transitions and Preserved Ratio Impaired Spirometry Among Ever Smokers

BACKGROUND

Emerging data from longitudinal studies suggest that preserved ratio impaired spirometry (PRISm), defined by proportionate reductions in FEV₁ and FVC, is a heterogeneous population with frequent transitions to other lung function categories relative to individuals with normal and obstructive spirometry. Controversy regarding the clinical significance of these transitions exists (eg, whether transitions merely reflect measurement variability or noise).

RESEARCH QUESTION

Are individuals with PRISm enriched for transitions associated with substantial changes in lung function?

STUDY DESIGN AND METHODS

Current and former smokers enrolled in the Genetic Epidemiology of COPD (COPDGene) study with spirometry available in phases 1 through 3 (enrollment, 5-year follow-up, and 10-year follow-up) were analyzed. Postbronchodilator lung function categories were as follows: PRISm (FEV₁ < 80% predicted with FEV₁/FVC ratio ≥ 0.7), Global Initiative for Chronic Obstructive Lung Disease grade 0 (FEV₁ ≥ 80% predicted and FEV₁/FVC ≥ 0.7), and obstruction (FEV₁/FVC < 0.7). Significant transition status was affirmative if a subject belonged to two or more spirometric categories and had > 10% change in FEV₁ % predicted and/or FVC % predicted between consecutive visits. Ever-PRISm was present if a subject had PRISm at any visit. Logistic regression examined the association between significant transitions and ever-PRISm status, adjusted for age, sex, race, FEV₁ % predicted, current smoking, pack-years, BMI, and ever-positive bronchodilator response.

RESULTS

Among subjects with complete data (N = 1,775) over 10.1 ± 0.4 years of follow-up, the prevalence of PRISm remained consistent (10.4%-11.3%) between phases 1 through 3, but nearly one-half of subjects with PRISm transitioned into or out of PRISm at each visit. Among all subjects, 19.7% had a significant transition; ever-PRISm was a significant predictor of significant transitions (unadjusted OR, 10.3; 95% CI, 7.9-13.5; adjusted OR, 14.9; 95% CI, 10.9-20.7). Results were similar with additional adjustment for radiographic emphysema and gas trapping, when lower limit of normal criteria were used to define lung function categories, and when FEV₁ alone (regardless of change in FVC % predicted) was used to define significant transitions.

INTERPRETATION

PRISm is an unstable group, with frequent significant transitions to both obstruction and normal spirometry over time.

CLINICAL TRIAL REGISTRATION

ClinicalTrials.gov; No.: NCT000608764; URL: www.

CLINICALTRIALS

gov.

Longitudinal association between muscle loss and mortality in ever-smokers

BACKGROUND

Body composition measures, specifically low weight or reduced muscle mass, are associated with mortality in chronic obstructive pulmonary disease (COPD), but the effect of longitudinal body composition changes is undefined.

RESEARCH QUESTION

Is the longitudinal loss of fat-free mass (FFM) associated with increased mortality including in those with initially normal or elevated body composition metrics?

STUDY DESIGN AND METHODS

Participants with complete data for at least one visit in the COPDGene (n=9,268) and ECLIPSE studies (1,760) were included and followed for 12 and 8 years, respectively. Pectoralis muscle area (PMA) was derived from thoracic CT scans and used as a proxy for FFM. A longitudinal mixed sub-model for PMA and a Cox proportional hazards sub-model for survival were fitted on a joint distribution using a shared random intercept parameter and Markov chain Monte Carlo parameter estimation.

RESULTS

Both cohorts demonstrated a left shifted distribution of baseline FFM, not reflected in BMI, and an increase in all-cause mortality risk associated with longitudinal loss of PMA. For each one cm2 PMA loss, mortality increased 3.1% (95% CI 2.4, 3.7, p<0.001) in COPDGene, and 2.4% (95% CI 0.9, 4.0, p<0.001) in ECLIPSE. Increased mortality risk was independent of enrollment values for BMI and disease severity (BODE index quartiles) and was significant even in participants with initially greater than average PMA.

INTERPRETATION

Longitudinal loss of PMA is associated with increased all-cause mortality, regardless of BMI or initial muscle mass. Consideration of novel screening tests and further research into mechanisms contributing to muscle decline may improve risk stratification and identify novel therapeutic targets in ever-smokers.

Increased mortality associated with frequent exacerbations in COPD patients with mild-to-moderate lung function impairment, and smokers with normal spirometry

Background:

Methods:

Results:

Conclusions:

Pulmonary Arterial Pruning and Longitudinal Change in Percent Emphysema and Lung Function: The Genetic Epidemiology of COPD Study

BACKGROUND

Pulmonary endothelial damage has been shown to precede the development of emphysema in animals, and vascular changes in humans have been observed in COPD and emphysema.

RESEARCH QUESTION

Is intraparenchymal vascular pruning associated with longitudinal progression of emphysema on CT imaging or decline in lung function over 5 years?

STUDY DESIGN AND METHODS

The Genetic Epidemiology of COPD Study enrolled ever smokers with and without COPD from 2008 through 2011. The percentage of emphysema-like lung, or "percent emphysema," was assessed at baseline and after 5 years on noncontrast CT imaging as the percentage of lung voxels < -950 Hounsfield units. An automated CT imaging-based tool assessed and classified intrapulmonary arteries and veins. Spirometry measures are postbronchodilator. Pulmonary arterial pruning was defined as a lower ratio of small artery volume (< 5 mm2 cross-sectional area) to total lung artery volume. Mixed linear models included demographics, anthropomorphics, smoking, and COPD, with emphysema models also adjusting for CT imaging scanner and lung function models adjusting for clinical center and baseline percent emphysema.

RESULTS

At baseline, the 4,227 participants were 60 ± 9 years of age, 50% were women, 28% were Black, 47% were current smokers, and 41% had COPD. Median percent emphysema was 2.1 (interquartile range, 0.6-6.3) and progressed 0.24 percentage points/y (95% CI, 0.22-0.26 percentage points/y) over 5.6 years. Mean FEV1 to FVC ratio was 68.5 ± 14.2% and declined 0.26%/y (95% CI, -0.30 to -0.23%/y). Greater pulmonary arterial pruning was associated with more rapid progression of percent emphysema (0.11 percentage points/y per 1-SD increase in arterial pruning; 95% CI, 0.09-0.16 percentage points/y), including after adjusting for baseline percent emphysema and FEV1. Arterial pruning also was associated with a faster decline in FEV1 to FVC ratio (-0.04%/y per 1-SD increase in arterial pruning; 95% CI, -0.008 to -0.001%/y).

INTERPRETATION

Pulmonary arterial pruning was associated with faster progression of percent emphysema and more rapid decline in FEV1 to FVC ratio over 5 years in ever smokers, suggesting that pulmonary vascular differences may be relevant in disease progression.

TRIAL REGISTRY

ClinicalTrials.gov; No.: NCT00608764; URL: www.clinicaltrials.gov.

Caution against examining the role of reverse causality in Mendelian Randomization

Recently, Mendelian Randomization (MR) has gained in popularity as a concept to assess the causal relationship between phenotypes in genetic association studies. An extension of standard MR methodology, the MR Steiger approach, has recently been developed to infer the causal direction between two phenotypes in prospective studies. Through simulation studies, we examined and quantified the ability of the MR Steiger approach to determine the causal direction between two phenotypes (i.e., effect direction). Through simulation studies, our results show that the MR Steiger approach may fail to correctly identify the direction of causality. This is true, especially in the presence of pleiotropy. We also applied the MR Steiger method to the COPDGene study, a case-control study of chronic obstructive pulmonary disease (COPD) in current and former smokers, to examine the role of smoking on lung function. We have created an R package on Github called reverseDirection which runs simulations for user-specified scenarios to examine when the MR Steiger approach can correctly determine the causal direction between two phenotypes in any user specified scenario. In summary, our results emphasize the importance of caution when the MR Steiger approach is used in to infer the direction of causality.

Respiratory exacerbations are associated with muscle loss in current and former smokers

OBJECTIVES

Muscle wasting is a recognised extra-pulmonary complication in chronic obstructive pulmonary disease and has been associated with increased risk of death. Acute respiratory exacerbations are associated with reduction of muscle function, but there is a paucity of data on their long-term effect. This study explores the relationship between acute respiratory exacerbations and long-term muscle loss using serial measurements of CT derived pectoralis muscle area (PMA).

DESIGN AND SETTING

Participants were included from two prospective, longitudinal, observational, multicentre cohorts of ever-smokers with at least 10 pack-year history.

PARTICIPANTS

The primary analysis included 1332 (of 2501) participants from Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints (ECLIPSE) and 4384 (of 10 198) participants from Genetic Epidemiology of COPD (COPDGene) who had complete data from their baseline and follow-up visits.

INTERVENTIONS

PMA was measured on chest CT scans at two timepoints. Self-reported exacerbation data were collected from participants in both studies through the use of periodic longitudinal surveys.

MAIN OUTCOME MEASURES

Age-related and excess muscle loss over time.

RESULTS

Age, sex, race and body mass index were associated with baseline PMA. Participants experienced age-related decline at the upper end of reported normal ranges. In ECLIPSE, the exacerbation rate over time was associated with an excess muscle area loss of 1.3% (95% CI 0.6 to 1.9, p<0.001) over 3 years and in COPDGene with an excess muscle area loss of 2.1% (95% CI 1.2 to 2.8, p<0.001) over 5 years. Excess muscle area decline was absent in 273 individuals who participated in pulmonary rehabilitation.

CONCLUSIONS

Exacerbations are associated with accelerated skeletal muscle loss. Each annual exacerbation was associated with the equivalent of 6 months of age-expected decline in muscle mass. Ameliorating exacerbation-associated muscle loss represents an important therapeutic target.

Expanding the genetic architecture of nicotine dependence and its shared genetics with multiple traits

Cigarette smoking is the leading cause of preventable morbidity and mortality. Genetic variation contributes to initiation, regular smoking, nicotine dependence, and cessation. We present a Fagerström Test for Nicotine Dependence (FTND)-based genome-wide association study in 58,000 European or African ancestry smokers. We observe five genome-wide significant loci, including previously unreported loci MAGI2/GNAI1 (rs2714700) and TENM2 (rs1862416), and extend loci reported for other smoking traits to nicotine dependence. Using the heaviness of smoking index from UK Biobank (N = 33,791), rs2714700 is consistently associated; rs1862416 is not associated, likely reflecting nicotine dependence features not captured by the heaviness of smoking index. Both variants influence nearby gene expression (rs2714700/MAGI2-AS3 in hippocampus; rs1862416/TENM2 in lung), and expression of genes spanning nicotine dependence-associated variants is enriched in cerebellum. Nicotine dependence (SNP-based heritability = 8.6%) is genetically correlated with 18 other smoking traits (rg = 0.40-1.09) and co-morbidities. Our results highlight nicotine dependence-specific loci, emphasizing the FTND as a composite phenotype that expands genetic knowledge of smoking.

Co-Morbidity Patterns Identified Using Latent Class Analysis of Medications Predict All-Cause Mortality Independent of Other Known Risk Factors: The COPDGene® Study

PURPOSE

Medication patterns include all medications in an individual's clinical profile. We aimed to identify chronic co-morbidity treatment patterns through medication use among COPDGene participants and determine whether these patterns were associated with mortality, acute exacerbations of chronic obstructive pulmonary disease (AECOPD) and quality of life.

MATERIALS AND METHODS

Participants analyzed here completed Phase 1 (P1) and/or Phase 2 (P2) of COPDGene. Latent class analysis (LCA) was used to identify medication patterns and assign individuals into unobserved LCA classes. Mortality, AECOPD, and the St. George's Respiratory Questionnaire (SGRQ) health status were compared in different LCA classes through survival analysis, logistic regression, and Kruskal-Wallis test, respectively.

RESULTS

LCA identified 8 medication patterns from 32 classes of chronic comorbid medications. A total of 8110 out of 10,127 participants with complete covariate information were included. Survival analysis adjusted for covariates showed, compared to a low medication use class, mortality was highest in participants with hypertension+diabetes+statin+antiplatelet medication group. Participants in hypertension+SSRI+statin medication group had the highest odds of AECOPD and the highest SGRQ score at both P1 and P2.

CONCLUSION

Medication pattern can serve as a good indicator of an individual's comorbidities profile and improves models predicting clinical outcomes.

Genome-Wide Gene-by-Smoking Interaction Study of Chronic Obstructive Pulmonary Disease

Risk of chronic obstructive pulmonary disease (COPD) is determined by both cigarette smoking and genetic susceptibility, but little is known about gene-by-smoking interactions. We performed a genome-wide association analysis of 179,689 controls and 21,077 COPD cases from UK Biobank subjects of European ancestry recruited from 2006 to 2010, considering genetic main effects and gene-by-smoking interaction effects simultaneously (2-degrees-of-freedom (df) test) as well as interaction effects alone (1-df interaction test). We sought to replicate significant results in COPDGene (United States, 2008-2010) and SpiroMeta Consortium (multiple countries, 1947-2015) data. We considered 2 smoking variables: 1) ever/never and 2) current/noncurrent. In the 1-df test, we identified 1 genome-wide significant locus on 15q25.1 (cholinergic receptor nicotinic β4 subunit, or CHRNB4) for ever- and current smoking and identified PI*Z allele (rs28929474) of serpin family A member 1 (SERPINA1) for ever-smoking and 3q26.2 (MDS1 and EVI1 complex locus, or MECOM) for current smoking in an analysis of previously reported COPD loci. In the 2-df test, most of the significant signals were also significant for genetic marginal effects, aside from 16q22.1 (sphingomyelin phosphodiesterase 3, or SMPD3) and 19q13.2 (Egl-9 family hypoxia inducible factor 2, or EGLN2). The significant effects at 15q25.1 and 19q13.2 loci, both previously described in prior genome-wide association studies of COPD or smoking, were replicated in COPDGene and SpiroMeta. We identified interaction effects at previously reported COPD loci; however, we failed to identify novel susceptibility loci.

Genome-Wide Meta-Analyses of FTND and TTFC Phenotypes

INTRODUCTION

FTND (Fagerstrӧm test for nicotine dependence) and TTFC (time to smoke first cigarette in the morning) are common measures of nicotine dependence (ND). However, genome-wide meta-analysis for these phenotypes has not been reported.

METHODS

Genome-wide meta-analyses for FTND (N = 19,431) and TTFC (N = 18,567) phenotypes were conducted for adult smokers of European ancestry from 14 independent cohorts.

RESULTS

We found that SORBS2 on 4q35 (p = 4.05 × 10-8), BG182718 on 11q22 (p = 1.02 × 10-8), and AA333164 on 14q21 (p = 4.11 × 10-9) were associated with TTFC phenotype. We attempted replication of leading candidates with independent samples (FTND, N = 7010 and TTFC, N = 10 061), however, due to limited power of the replication samples, the replication of these new loci did not reach significance. In gene-based analyses, COPB2 was found associated with FTND phenotype, and TFCP2L1, RELN, and INO80C were associated with TTFC phenotype. In pathway and network analyses, we found that the interconnected interactions among the endocytosis, regulation of actin cytoskeleton, axon guidance, MAPK signaling, and chemokine signaling pathways were involved in ND.

CONCLUSIONS

Our analyses identified several promising candidates for both FTND and TTFC phenotypes, and further verification of these candidates was necessary. Candidates supported by both FTND and TTFC (CHRNA4, THSD7B, RBFOX1, and ZNF804A) were associated with addiction to alcohol, cocaine, and heroin, and were associated with autism and schizophrenia. We also identified novel pathways involved in cigarette smoking. The pathway interactions highlighted the importance of receptor recycling and internalization in ND.

IMPLICATIONS

Understanding the genetic architecture of cigarette smoking and ND is critical to develop effective prevention and treatment. Our study identified novel candidates and biological pathways involved in FTND and TTFC phenotypes, and this will facilitate further investigation of these candidates and pathways.

Dissecting the genetic overlap of smoking behaviors, lung cancer, and chronic obstructive pulmonary disease: A focus on nicotinic receptors and nicotine metabolizing enzyme

Smoking is a major contributor to lung cancer and chronic obstructive pulmonary disease (COPD). Two of the strongest genetic associations of smoking-related phenotypes are the chromosomal regions 15q25.1, encompassing the nicotinic acetylcholine receptor subunit genes CHRNA5-CHRNA3-CHRNB4, and 19q13.2, encompassing the nicotine metabolizing gene CYP2A6. In this study, we examined genetic relations between cigarettes smoked per day, smoking cessation, lung cancer, and COPD. Data consisted of genome-wide association study summary results. Genetic correlations were estimated using linkage disequilibrium score regression software. For each pair of outcomes, z-score-z-score (ZZ) plots were generated. Overall, heavier smoking and decreased smoking cessation showed positive genetic associations with increased lung cancer and COPD risk. The chromosomal region 19q13.2, however, showed a different correlational pattern. For example, the effect allele-C of the sentinel SNP (rs56113850) within CYP2A6 was associated with an increased risk of heavier smoking (z-score = 19.2; p = 1.10 × 10-81 ), lung cancer (z-score = 8.91; p = 5.02 × 10-19 ), and COPD (z-score = 4.04; p = 5.40 × 10-5 ). Surprisingly, this allele-C (rs56113850) was associated with increased smoking cessation (z-score = -8.17; p = 2.52 × 10-26 ). This inverse relationship highlights the need for additional investigation to determine how CYP2A6 variation could increase smoking cessation while also increasing the risk of lung cancer and COPD likely through increased cigarettes smoked per day.

A Risk Prediction Model for Mortality Among Smokers in the COPDGene® Study

BACKGROUND

Risk factor identification is a proven strategy in advancing treatments and preventive therapy for many chronic conditions. Quantifying the impact of those risk factors on health outcomes can consolidate and focus efforts on individuals with specific high-risk profiles. Using multiple risk factors and longitudinal outcomes in 2 independent cohorts, we developed and validated a risk score model to predict mortality in current and former cigarette smokers.

METHODS

We obtained extensive data on current and former smokers from the COPD Genetic Epidemiology (COPDGene®) study at enrollment. Based on physician input and model goodness-of-fit measures, a subset of variables was selected to fit final Weibull survival models separately for men and women. Coefficients and predictors were translated into a point system, allowing for easy computation of mortality risk scores and probabilities. We then used the SubPopulations and InteRmediate Outcome Measures In COPD Study (SPIROMICS) cohort for external validation of our model.

RESULTS

Of 9867 COPDGene participants with standard baseline data, 17.6% died over 10 years of follow-up, and 9074 of these participants had the full set of baseline predictors (standard plus 6-minute walk distance and computed tomography variables) available for full model fits. The average age of participants in the cohort was 60 for both men and women, and the average predicted 10-year mortality risk was 18% for women and 25% for men. Model time-integrated area under the receiver operating characteristic curve statistics demonstrated good predictive model accuracy (0.797 average), validated in the external cohort (0.756 average). Risk of mortality was impacted most by 6-minute walk distance, forced expiratory volume in 1 second and age, for both men and women.

CONCLUSIONS

Current and former smokers exhibited a wide range of mortality risk over a 10- year period. Our models can identify higher risk individuals who can be targeted for interventions to reduce risk of mortality, for participants with or without chronic obstructive pulmonary disease (COPD) using current Global initiative for obstructive Lung Disease (GOLD) criteria.

Chronic obstructive pulmonary disease and related phenotypes: polygenic risk scores in population-based and case-control cohorts

BACKGROUND

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.

METHODS

We constructed a polygenic risk score using a genome-wide association study of lung function (FEV1 and FEV1/forced vital capacity [FVC]) from the UK Biobank and SpiroMeta. We tested this polygenic risk score in nine cohorts of multiple ethnicities for an association with moderate-to-severe COPD (defined as FEV1/FVC <0·7 and FEV1 <80% of predicted). Associations were tested using logistic regression models, adjusting for age, sex, height, smoking pack-years, and principal components of genetic ancestry. We assessed predictive performance of models by area under the curve. In a subset of studies, we also studied quantitative and qualitative CT imaging phenotypes that reflect parenchymal and airway pathology, and patterns of reduced lung growth.

FINDINGS

The polygenic risk score was associated with COPD in European (odds ratio [OR] per SD 1·81 [95% CI 1·74-1·88] and non-European (1·42 [1·34-1·51]) populations. Compared with the first decile, the tenth decile of the polygenic risk score was associated with COPD, with an OR of 7·99 (6·56-9·72) in European ancestry and 4·83 (3·45-6·77) in non-European ancestry cohorts. The polygenic risk score was superior to previously described genetic risk scores and, when combined with clinical risk factors (ie, age, sex, and smoking pack-years), showed improved prediction for COPD compared with a model comprising clinical risk factors alone (AUC 0·80 [0·79-0·81] vs 0·76 [0·75-0·76]). The polygenic risk score was associated with CT imaging phenotypes, including wall area percent, quantitative and qualitative measures of emphysema, local histogram emphysema patterns, and destructive emphysema subtypes. The polygenic risk score was associated with a reduced lung growth pattern.

INTERPRETATION

A risk score comprised of genetic variants can identify a small subset of individuals at markedly increased risk for moderate-to-severe COPD, emphysema subtypes associated with cigarette smoking, and patterns of reduced lung growth.

FUNDING

US National Institutes of Health, Wellcome Trust.

Association Analysis and Meta-Analysis of Multi-Allelic Variants for Large-Scale Sequence Data

There is great interest in understanding the impact of rare variants in human diseases using large sequence datasets. In deep sequence datasets of >10,000 samples, ~10% of the variant sites are observed to be multi-allelic. Many of the multi-allelic variants have been shown to be functional and disease-relevant. Proper analysis of multi-allelic variants is critical to the success of a sequencing study, but existing methods do not properly handle multi-allelic variants and can produce highly misleading association results. We discuss practical issues and methods to encode multi-allelic sites, conduct single-variant and gene-level association analyses, and perform meta-analysis for multi-allelic variants. We evaluated these methods through extensive simulations and the study of a large meta-analysis of ~18,000 samples on the cigarettes-per-day phenotype. We showed that our joint modeling approach provided an unbiased estimate of genetic effects, greatly improved the power of single-variant association tests among methods that can properly estimate allele effects, and enhanced gene-level tests over existing approaches. Software packages implementing these methods are available online.

Smaller Left Ventricle Size at Noncontrast CT Is Associated with Lower Mortality in COPDGene Participants

Background Smokers with chronic obstructive pulmonary disease (COPD) have smaller left ventricles (LVs) due to reduced preload. Skeletal muscle wasting is also common in COPD, but less is known about its contribution to LV size. Purpose To explore the relationships between CT metrics of emphysema, venous vascular volume, and sarcopenia with the LV epicardial volume (LV_EV) (myocardium and chamber) estimated from chest CT images in participants with COPD and then to determine the clinical relevance of the LV_EV in multivariable models, including sex and anthropomorphic metrics. Materials and Methods The COPDGene study (ClinicalTrials.gov identifier: NCT00608764) is an ongoing prospective longitudinal observational investigation that began in 2006. LV_EV, distal pulmonary venous blood volume for vessels smaller than 5 mm² in cross section (BV5), CT emphysema, and pectoralis muscle area were retrospectively extracted from 3318 nongated, unenhanced COPDGene CT scans. Multivariable linear and Cox regression models were used to explore the association between emphysema, venous BV5, pectoralis muscle area, and LV_EV as well as the association of LV_EV with health status using the St George's Respiratory Questionnaire, 6-minute walk distance, and all-cause mortality. Results The median age of the cohort was 64 years (interquartile range, 57-70 years). Of the 2423 participants, 1806 were men and 617 were African American. The median LV_EV between Global Initiative for Chronic Obstructive Lung Disease (GOLD) 1 and GOLD 4 COPD was reduced by 13.9% in women and 17.7% in men (P < .001 for both). In fully adjusted models, higher emphysema percentage (β = -4.2; 95% confidence interval [CI]: -5.0, -3.4; P < .001), venous BV5 (β = 7.0; 95% CI: 5.7, 8.2; P < .001), and pectoralis muscle area (β = 2.7; 95% CI: 1.2, 4.1; P < .001) were independently associated with reduced LV_EV. Reductions in LV_EV were associated with improved health status (β = 0.3; 95% CI: 0.1, 0.4) and 6-minute walk distance (β = -12.2; 95% CI: -15.2, -9.3). These effects were greater in women than in men. The effect of reduced LV_EV on mortality (hazard ratio: 1.07; 95% CI: 1.05, 1.09) did not vary by sex. Conclusion In women more than men with chronic obstructive pulmonary disease, a reduction in the estimated left ventricle epicardial volume correlated with a loss of pulmonary venous vasculature, greater pectoralis muscle sarcopenia, and lower all-cause mortality. © RSNA, 2020 Online supplemental material is available for this article.

Common and Rare Variants Genetic Association Analysis of Cigarettes per Day Among Ever-Smokers in Chronic Obstructive Pulmonary Disease Cases and Controls

INTRODUCTION

Cigarette smoking is a major environmental risk factor for many diseases, including chronic obstructive pulmonary disease (COPD). There are shared genetic influences on cigarette smoking and COPD. Genetic risk factors for cigarette smoking in cohorts enriched for COPD are largely unknown.

METHODS

We performed genome-wide association analyses for average cigarettes per day (CPD) across the Genetic Epidemiology of COPD (COPDGene) non-Hispanic white (NHW) (n = 6659) and African American (AA) (n = 3260), GenKOLS (the Genetics of Chronic Obstructive Lung Disease) (n = 1671), and ECLIPSE (the Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints) (n = 1942) cohorts. In addition, we performed exome array association analyses across the COPDGene NHW and AA cohorts. We considered analyses across the entire cohort and stratified by COPD case-control status.

RESULTS

We identified genome-wide significant associations for CPD on chromosome 15q25 across all cohorts (lowest p = 1.78 × 10-15), except in the COPDGene AA cohort alone. Previously reported associations on chromosome 19 had suggestive and directionally consistent associations (RAB4, p = 1.95 × 10-6; CYP2A7, p = 7.50 × 10-5; CYP2B6, p = 4.04 × 10-4). When we stratified by COPD case-control status, single nucleotide polymorphisms on chromosome 15q25 were nominally associated with both NHW COPD cases (β = 0.11, p = 5.58 × 10-4) and controls (β = 0.12, p = 3.86 × 10-5) For the gene-based exome array association analysis of rare variants, there were no exome-wide significant associations. For these previously replicated associations, the most significant results were among COPDGene NHW subjects for CYP2A7 (p = 5.2 × 10-4).

CONCLUSIONS

In a large genome-wide association study of both common variants and a gene-based association of rare coding variants in ever-smokers, we found genome-wide significant associations on chromosome 15q25 with CPD for common variants, but not for rare coding variants. These results were directionally consistent among COPD cases and controls.

IMPLICATIONS

We examined both common and rare coding variants associated with CPD in a large population of heavy smokers with and without COPD of NHW and AA descent. We replicated genome-wide significant associations on chromosome 15q25 with CPD for common variants among NHW subjects, but not for rare variants. We demonstrated for the first time that common variants on chromosome 15q25 associated with CPD are similar among COPD cases and controls. Previously reported associations on chromosome 19 showed suggestive and directionally consistent associations among common variants (RAB4, CYP2A7, and CYP2B6) and for rare variants (CYP2A7) among COPDGene NHW subjects. Although the genetic effect sizes for these single nucleotide polymorphisms on chromosome 15q25 are modest, we show that this creates a substantial smoking burden over the lifetime of a smoker.

Genetic Advances in Chronic Obstructive Pulmonary Disease. Insights from COPDGene

Chronic obstructive pulmonary disease (COPD) is a common and progressive disease that is influenced by both genetic and environmental factors. For many years, knowledge of the genetic basis of COPD was limited to Mendelian syndromes, such as alpha-1 antitrypsin deficiency and cutis laxa, caused by rare genetic variants. Over the past decade, the proliferation of genome-wide association studies, the accessibility of whole-genome sequencing, and the development of novel methods for analyzing genetic variation data have led to a substantial increase in the understanding of genetic variants that play a role in COPD susceptibility and COPD-related phenotypes. COPDGene (Genetic Epidemiology of COPD), a multicenter, longitudinal study of over 10,000 current and former cigarette smokers, has been pivotal to these breakthroughs in understanding the genetic basis of COPD. To date, over 20 genetic loci have been convincingly associated with COPD affection status, with additional loci demonstrating association with COPD-related phenotypes such as emphysema, chronic bronchitis, and hypoxemia. In this review, we discuss the contributions of the COPDGene study to the discovery of these genetic associations as well as the ongoing genetic investigations of COPD subtypes, protein biomarkers, and post-genome-wide association study analysis.

Five-year Progression of Emphysema and Air Trapping at CT in Smokers with and Those without Chronic Obstructive Pulmonary Disease: Results from the COPDGene Study

Background CT is used to quantify abnormal changes in the lung parenchyma of smokers that might overlap chronic obstructive pulmonary disease (COPD), but studies on the progression of expiratory air trapping in smokers are scarce. Purpose To evaluate the relationship between longitudinal changes in forced expiratory volume in 1 second (FEV₁) and CT-quantified emphysema and air trapping in smokers. Materials and Methods Cigarette smokers with and those without COPD participating in the multicenter observational COPDGene study were evaluated. Subjects underwent inspiratory and expiratory chest CT and spirometry at baseline and 5-year follow-up. Emphysema was quantified by using adjusted lung density (ALD). Air trapping was quantified by using mean lung density at expiratory CT and CT-measured functional residual capacity-to-total lung volume ratio. Linear models were used to regress quantitative CT measurements taken 5 years apart, and models were fit with and without adding FEV₁ as a predictor. Analyses were stratified by Global Initiative for Chronic Obstructive Lung Disease (GOLD) stage (GOLD 0, no COPD; GOLD 1, mild COPD; GOLD 2, moderate COPD; GOLD 3, severe COPD; GOLD 4, very severe COPD). Subjects with preserved FEV₁-to-forced vital capacity ratio and reduced FEV₁ percentage predicted were categorized as having preserved ratio impaired spirometry (PRISm). Results A total of 4211 subjects (503 with PRISm; 2034 with GOLD 0, 388 with GOLD 1, 816 with GOLD 2, 381 with GOLD 3, 89 with GOLD 4) were evaluated. ALD decreased by 1.7 g/L (95% confidence interval [CI]: -2.5, -0.9) in subjects with GOLD 0 at baseline and by 5.3 g/L (95% CI: -6.2, -4.4) in those with GOLD 1-4 (P < .001 for both). When adjusted for changes in FEV₁, corresponding numbers were -2.2 (95% CI: -3.0, -1.3) and -4.6 g/L (95% CI: -5.6, -3.4) (P < .001 for both). Progression in air trapping was identified only in GOLD stage 2-4. Approximately 33%-50% of changes in air trapping in GOLD stages 2-4 were accounted for by changes in FEV₁. Conclusion CT measures of emphysema and air trapping increased over 5 years in smokers. Forced expiratory volume in one second accounted for less than 10% of emphysema progression and less than 50% of air trapping progression detected at CT. © RSNA, 2020 Online supplemental material is available for this article.

The Association of Multiparity with Lung Function and Chronic Obstructive Pulmonary Disease-Related Phenotypes

BACKGROUND

Apparent increased female susceptibility to chronic obstructive pulmonary disease (COPD) suggests sex hormones modulate disease pathogenesis. Little is known about associations between multiparity and lung function in smokers.

RESEARCH QUESTION

We hypothesized that multiparity is associated with lung function and measures of emphysema and airway disease.

STUDY DESIGN AND METHODS

Utilizing female participants from the 5-year follow up of the COPD Genetic Epidemiology (COPDGene^®) study we performed multivariable linear regressions to assess the effect of multiparity and number of pregnancies on forced expiratory volume in 1 second (FEV₁) percentage of predicted (% predicted), FEV₁/forced vital capacity (FVC), percent emphysema on computed tomography (CT) scans, and Pi10, a measure of airway thickening. We sampled never smokers and those with lower smoking exposure from the National Health and Nutrition Examination Survey (NHANES) 2011-2012 dataset.

RESULTS

We included 1820 participants from COPDGene^® and 418 participants from NHANES (321 never smokers, 97 ever smokers). In COPDGene^®, multiparity (beta coefficient [β] = -3.8, 95% confidence interval [CI]: [-6.5, -1.1], p = 0.005) and higher number of pregnancies were associated with lower FEV₁ % predicted. Multiparity was not associated with percent emphysema or Pi10. In individuals with no or mild obstruction, multiparity was associated with lower FEV₁ % predicted. There was an interaction with multiparity and age on FEV₁ % predicted (p = 0.025). In NHANES, there was no association between multiparity and FEV₁ % predicted in never smokers or the lower smoking exposure group.

INTERPRETATION

Multiparity was associated with lower FEV₁ % predicted in current and former smokers in COPDGene^® study participants. These preliminary results emphasize the importance of smoking abstinence in women of child-bearing age.

Machine Learning Characterization of COPD Subtypes: Insights From the COPDGene Study

COPD is a heterogeneous syndrome. Many COPD subtypes have been proposed, but there is not yet consensus on how many COPD subtypes there are and how they should be defined. The COPD Genetic Epidemiology Study (COPDGene), which has generated 10-year longitudinal chest imaging, spirometry, and molecular data, is a rich resource for relating COPD phenotypes to underlying genetic and molecular mechanisms. In this article, we place COPDGene clustering studies in context with other highly cited COPD clustering studies, and summarize the main COPD subtype findings from COPDGene. First, most manifestations of COPD occur along a continuum, which explains why continuous aspects of COPD or disease axes may be more accurate and reproducible than subtypes identified through clustering methods. Second, continuous COPD-related measures can be used to create subgroups through the use of predictive models to define cut-points, and we review COPDGene research on blood eosinophil count thresholds as a specific example. Third, COPD phenotypes identified or prioritized through machine learning methods have led to novel biological discoveries, including novel emphysema genetic risk variants and systemic inflammatory subtypes of COPD. Fourth, trajectory-based COPD subtyping captures differences in the longitudinal evolution of COPD, addressing a major limitation of clustering analyses that are confounded by disease severity. Ongoing longitudinal characterization of subjects in COPDGene will provide useful insights about the relationship between lung imaging parameters, molecular markers, and COPD progression that will enable the identification of subtypes based on underlying disease processes and distinct patterns of disease progression, with the potential to improve the clinical relevance and reproducibility of COPD subtypes.

COPDGene® 2019: Redefining the Diagnosis of Chronic Obstructive Pulmonary Disease

BACKGROUND

Chronic obstructive pulmonary disease (COPD) remains a major cause of morbidity and mortality. Present-day diagnostic criteria are largely based solely on spirometric criteria. Accumulating evidence has identified a substantial number of individuals without spirometric evidence of COPD who suffer from respiratory symptoms and/or increased morbidity and mortality. There is a clear need for an expanded definition of COPD that is linked to physiologic, structural (computed tomography [CT]) and clinical evidence of disease. Using data from the COPD Genetic Epidemiology study (COPDGene®), we hypothesized that an integrated approach that includes environmental exposure, clinical symptoms, chest CT imaging and spirometry better defines disease and captures the likelihood of progression of respiratory obstruction and mortality.

METHODS

Four key disease characteristics - environmental exposure (cigarette smoking), clinical symptoms (dyspnea and/or chronic bronchitis), chest CT imaging abnormalities (emphysema, gas trapping and/or airway wall thickening), and abnormal spirometry - were evaluated in a group of 8784 current and former smokers who were participants in COPDGene® Phase 1. Using these 4 disease characteristics, 8 categories of participants were identified and evaluated for odds of spirometric disease progression (FEV1 > 350 ml loss over 5 years), and the hazard ratio for all-cause mortality was examined.

RESULTS

Using smokers without symptoms, CT imaging abnormalities or airflow obstruction as the reference population, individuals were classified as Possible COPD, Probable COPD and Definite COPD. Current Global initiative for obstructive Lung Disease (GOLD) criteria would diagnose 4062 (46%) of the 8784 study participants with COPD. The proposed COPDGene® 2019 diagnostic criteria would add an additional 3144 participants. Under the new criteria, 82% of the 8784 study participants would be diagnosed with Possible, Probable or Definite COPD. These COPD groups showed increased risk of disease progression and mortality. Mortality increased in patients as the number of their COPD characteristics increased, with a maximum hazard ratio for all cause-mortality of 5.18 (95% confidence interval [CI]: 4.15-6.48) in those with all 4 disease characteristics.

CONCLUSIONS

A substantial portion of smokers with respiratory symptoms and imaging abnormalities do not manifest spirometric obstruction as defined by population normals. These individuals are at significant risk of death and spirometric disease progression. We propose to redefine the diagnosis of COPD through an integrated approach using environmental exposure, clinical symptoms, CT imaging and spirometric criteria. These expanded criteria offer the potential to stimulate both current and future interventions that could slow or halt disease progression in patients before disability or irreversible lung structural changes develop.

Arterial Vascular Pruning, Right Ventricular Size, and Clinical Outcomes in Chronic Obstructive Pulmonary Disease. A Longitudinal Observational Study

Rationale: Cor pulmonale (right ventricular [RV] dilation) and cor pulmonale parvus (RV shrinkage) are both described in chronic obstructive pulmonary disease (COPD). The identification of emphysema as a shared risk factor suggests that additional disease characterization is needed to understand these widely divergent cardiac processes.Objectives: To explore the relationship between computed tomography measures of emphysema and distal pulmonary arterial morphology with RV volume, and their association with exercise capacity and mortality in ever-smokers with COPD enrolled in the COPDGene Study.Methods: Epicardial (myocardium and chamber) RV volume (RVEV), distal pulmonary arterial blood vessel volume (arterial BV5: vessels <5 mm2 in cross-section), and objective measures of emphysema were extracted from 3,506 COPDGene computed tomography scans. Multivariable linear and Cox regression models and the log-rank test were used to explore the association between emphysema, arterial BV5, and RVEV with exercise capacity (6-min-walk distance) and all-cause mortality.Measurements and Main Results: The RVEV was approximately 10% smaller in Global Initiative for Chronic Obstructive Lung Disease stage 4 versus stage 1 COPD (P < 0.0001). In multivariable modeling, a 10-ml decrease in arterial BV5 (pruning) was associated with a 1-ml increase in RVEV. For a given amount of emphysema, relative preservation of the arterial BV5 was associated with a smaller RVEV. An increased RVEV was associated with reduced 6-minute-walk distance and in those with arterial pruning an increased mortality.Conclusions: Pulmonary arterial pruning is associated with clinically significant increases in RV volume in smokers with COPD and is related to exercise capacity and mortality in COPD.Clinical trial registered with www.clinicaltrials.gov (NCT00608764).

Increased apolipoprotein C3 drives cardiovascular risk in type 1 diabetes

Type 1 diabetes mellitus (T1DM) increases the risk of atherosclerotic cardiovascular disease (CVD) in humans by poorly understood mechanisms. Using mouse models of T1DM-accelerated atherosclerosis, we found that relative insulin deficiency rather than hyperglycemia elevated levels of apolipoprotein C3 (APOC3), an apolipoprotein that prevents clearance of triglyceride-rich lipoproteins (TRLs) and their remnants. We then showed that serum APOC3 levels predict incident CVD events in subjects with T1DM in the Coronary Artery Calcification in Type 1 Diabetes (CACTI) study. To explore underlying mechanisms, we investigated the impact of Apoc3 antisense oligonucleotides (ASOs) on lipoprotein metabolism and atherosclerosis in a mouse model of T1DM. Apoc3 ASO treatment abolished the increased hepatic Apoc3 expression in diabetic mice - resulting in lower levels of TRLs - without improving glycemic control. APOC3 suppression also prevented arterial accumulation of APOC3-containing lipoprotein particles, macrophage foam cell formation, and the accelerated atherosclerosis in diabetic mice. Our observations demonstrate that relative insulin deficiency increases APOC3 and that this results in elevated levels of TRLs and accelerated atherosclerosis in a mouse model of T1DM. Because serum levels of APOC3 predicted incident CVD events in the CACTI study, inhibiting APOC3 might reduce CVD risk in T1DM patients.

Turning subtypes into disease axes to improve prediction of COPD progression

Chronic obstructive pulmonary disease (COPD) is an umbrella definition encompassing multiple disease processes. COPD heterogeneity has been described as distinct subgroups of individuals (subtypes) or as continuous measures of COPD variability (disease axes). There is little consensus on whether subtypes or disease axes are preferred, and the relative value of disease axes and subtypes for predicting COPD progression is unknown. Using a propensity score approach to learn disease axes from pairs of subtypes, we demonstrate that these disease axes predict prospective forced expiratory volume in 1 s decline and emphysema progression more accurately than the subtype pairs from which they were derived.

Author Correction: New genetic signals for lung function highlight pathways and chronic obstructive pulmonary disease associations across multiple ancestries

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Apolipoprotein E-C1-C4-C2 gene cluster region and inter-individual variation in plasma lipoprotein levels: a comprehensive genetic association study in two ethnic groups

The apolipoprotein E-C1-C4-C2 gene cluster at 19q13.32 encodes four amphipathic apolipoproteins. The influence of APOE common polymorphisms on plasma lipid/lipoprotein profile, especially on LDL-related traits, is well recognized; however, little is known about the role of other genes/variants in this gene cluster. In this study, we evaluated the role of common and uncommon/rare genetic variation in this gene region on inter-individual variation in plasma lipoprotein levels in non-Hispanic Whites (NHWs) and African blacks (ABs). In the variant discovery step, the APOE, APOC1, APOC4, APOC2 genes were sequenced along with their flanking and hepatic control regions (HCR1 and HCR2) in 190 subjects with extreme HDL-C/TG levels. The next step involved the genotyping of 623 NHWs and 788 ABs for the identified uncommon/rare variants and common tagSNPs along with additional relevant SNPs selected from public resources, followed by association analyses with lipid traits. A total of 230 sequence variants, including 15 indels were identified, of which 65 were novel. A total of 70 QC-passed variants in NHWs and 108 QC-passed variants in ABs were included in the final association analyses. Single-site association analysis of SNPs with MAF>1% revealed 20 variants in NHWs and 24 variants in ABs showing evidence of association with at least one lipid trait, including several variants exhibiting independent associations from the established APOE polymorphism even after multiple-testing correction. Overall, our study has confirmed known associations and also identified novel associations in this genomic region with various lipid traits. Our data also support the contribution of both common and uncommon/rare variation in this gene region in affecting plasma lipid profile in the general population.