Search

Summary: Background: Chronic Obstructive Pulmonary Disease (COPD) is a leading cause of mortality. Predicting mortality risk in patients with COPD can be important for disease management strategies. Although all-cause mortality predictors have been developed previously, limited research exists on factors directly affecting COPD-specific mortality. Methods: In a retrospective study, we used probabilistic graphs to analyse clinical cross-sectional data (COPDGene cohort), including demographics, spirometry, quantitative chest imaging, and symptom features, as well as gene expression data. COPDGene recruited current and former smokers, aged 45–80 years with >10 pack-years smoking history, from across the USA (Phase 1, 11/2007-4/2011) and invited them for a follow-up visit (Phase 2, 7/2013-7/2017). ECLIPSE cohort recruited current and former smokers (COPD patients and controls from USA and Europe), aged 45–80 with smoking history >10 pack-years (12/2005-11/2007). We applied graphical models on multi-modal data COPDGene Phase 1 participants to identify factors directly affecting all-cause and COPD-specific mortality (primary outcomes); and on Phase 2 follow-up cohort to identify additional molecular and social factors affecting mortality. We used penalized Cox regression with features selected by the causal graph to build VAPORED, a mortality risk prediction model. VAPORED was compared to existing scores (BODE: BMI, airflow obstruction, dyspnoea, exercise capacity; ADO: age, dyspnoea, airflow obstruction) on the ability to rank individuals by mortality risk, using four evaluation metrics (concordance, concordance probability estimate (CPE), cumulative/dynamic (C/D) area under the receiver operating characteristic curve (AUC), and integrated C/D AUC). The results were validated in ECLIPSE. Findings: Graphical models, applied on the COPDGene Phase 1 samples (n = 8610), identified 11 and 7 variables directly linked to all-cause and COPD-specific mortality, respectively. Although many appear in both models, non-lung comorbidities appear only in the all-cause model, while forced vital capacity (FVC %predicted) appears in COPD-specific mortality model only. Additionally, the graph model of Phase 2 data (n = 3182) identified internet access, CD4 T cells and platelets to be linked to lower mortality risk. Furthermore, using the 7 variables linked to COPD-specific mortality (forced expiratory volume in 1 s/forced vital capacity (FEV1/FVC) ration, FVC %predicted, age, history of pneumonia, oxygen saturation, 6-min walk distance, dyspnoea) we developed VAPORED mortality risk score, which we validated on the ECLIPSE cohort (3-yr all-cause mortality data, n = 2312). VAPORED performed significantly better than ADO, BODE, and updated BODE indices in predicting all-cause mortality in ECLIPSE in terms of concordance (VAPORED [0.719] vs ADO [0.693; FDR p-value 0.014], BODE [0.695; FDR p-value 0.020], and updated BODE [0.694; FDR p-value 0.021]); CPE (VAPORED [0.714] vs ADO [0.673; FDR p-value

Background Cardiovascular disease (CVD) is the most important comorbidity in patients with chronic obstructive pulmonary disease (COPD). COPD exacerbations not only contribute to COPD progression but may also elevate the risk of CVD. This study aimed to determine whether COPD exacerbations increase the risk of subsequent CVD events using up to 15 years of prospective longitudinal follow‐up data from the COPDGene (Genetic Epidemiology of Chronic Obstructive Pulmonary Disease) study. Methods and Results The COPDGene study is a large, multicenter, longitudinal investigation of COPD, including subjects at enrollment aged 45 to 80 years with a minimum of 10 pack‐years of smoking history. Cox proportional hazards models and Kaplan‐Meier survival curves were used to assess the risk of a composite end point of CVD based on the COPD exacerbation rate. Frequent exacerbators exhibited a higher cumulative incidence of composite CVD end points than infrequent exacerbators, irrespective of the presence of CVD at baseline. After adjusting for covariates, frequent exacerbators still maintained higher hazard ratios (HRs) than the infrequent exacerbator group (without CVD: HR, 1.81 [95% CI, 1.47–2.22]; with CVD: HR, 1.92 [95% CI, 1.51–2.44]). This observation remained consistently significant in moderate to severe COPD subjects and the preserved ratio impaired spirometry population. In the mild COPD population, frequent exacerbators showed a trend toward more CVD events. Conclusions COPD exacerbations are associated with an increased risk of subsequent cardiovascular events in subjects with and without preexisting CVD. Patients with COPD experiencing frequent exacerbations may necessitate careful monitoring and additional management for subsequent potential CVD. Registration URL: https://www.clinicaltrials.gov; Unique identifier: NCT00608764.

Chronic Obstructive Pulmonary Disease (COPD) is a heterogeneous, chronic inflammatory process of the lungs and, like other complex diseases, is caused by both genetic and environmental factors. Detailed understanding of the molecular mechanisms of complex diseases requires the study of the interplay among different biomolecular layers, and thus the integration of different omics data types. In this study, we investigated COPD-associated molecular mechanisms through a correlation-based network integration of lung tissue RNA-seq and DNA methylation data of COPD cases (n = 446) and controls (n = 346) derived from the Lung Tissue Research Consortium. First, we performed a SWIM-network based analysis to build separate correlation networks for RNA-seq and DNA methylation data for our case-control study population. Then, we developed a method to integrate the results into a coupled network of differentially expressed and differentially methylated genes to investigate their relationships across both molecular layers. The functional enrichment analysis of the nodes of the coupled network revealed a strikingly significant enrichment in Immune System components, both innate and adaptive, as well as immune-system component communication (interleukin and cytokine-cytokine signaling). Our analysis allowed us to reveal novel putative COPD-associated genes and to analyze their relationships, both at the transcriptomics and epigenomics levels, thus contributing to an improved understanding of COPD pathogenesis.

Abstract Background Quantitative interstitial abnormalities (QIA) are an automated computed tomography (CT) finding of early parenchymal lung disease, associated with worse lung function, reduced exercise capacity, increased respiratory symptoms, and death. The metabolomic perturbations associated with QIA are not well known. We sought to identify plasma metabolites associated with QIA in smokers. We also sought to identify shared and differentiating metabolomics features between QIA and emphysema, another smoking-related advanced radiographic abnormality. Methods In 928 former and current smokers in the Genetic Epidemiology of COPD cohort, we measured QIA and emphysema using an automated local density histogram method and generated metabolite profiles from plasma samples using liquid chromatography–mass spectrometry (Metabolon). We assessed the associations between metabolite levels and QIA using multivariable linear regression models adjusted for age, sex, body mass index, smoking status, pack-years, and inhaled corticosteroid use, at a Benjamini–Hochberg False Discovery Rate p-value of ≤ 0.05. Using multinomial regression models adjusted for these covariates, we assessed the associations between metabolite levels and the following CT phenotypes: QIA-predominant, emphysema-predominant, combined-predominant, and neither- predominant. Pathway enrichment analyses were performed using MetaboAnalyst. Results We found 85 metabolites significantly associated with QIA, with overrepresentation of the nicotinate and nicotinamide, histidine, starch and sucrose, pyrimidine, phosphatidylcholine, lysophospholipid, and sphingomyelin pathways. These included metabolites involved in inflammation and immune response, extracellular matrix remodeling, surfactant, and muscle cachexia. There were 75 metabolites significantly different between QIA-predominant and emphysema-predominant phenotypes, with overrepresentation of the phosphatidylethanolamine, nicotinate and nicotinamide, aminoacyl-tRNA, arginine, proline, alanine, aspartate, and glutamate pathways. Conclusions Metabolomic correlates may lend insight to the biologic perturbations and pathways that underlie clinically meaningful quantitative CT measurements like QIA in smokers.

Abstract The consequences of returning infectious pathogen test results identified incidentally in research studies have not been well-studied. Concerns include identification of an important health issue for individuals, accuracy of research test results, public health impact, potential emotional distress for participants, and need for IRB permissions. Blood RNA-sequencing analysis for non-human RNA in 3984 participants from the COPDGene study identified 228 participants with evidence suggestive for hepatitis C virus (HCV) infection. We hypothesized that incidentally discovered HCV results could be effectively returned to COPDGene participants with attention to the identified concerns. In conjunction with a COPDGene Participant Advisory Panel, we developed and obtained IRB approval for a process of returning HCV research results and an HCV Follow-Up Study questionnaire to capture information about previous HCV diagnosis and treatment information and participant reactions to return of HCV results. During phone calls following the initial HCV notification letter, 84 of 124 participants who could be contacted (67.7%) volunteered that they had been previously diagnosed with HCV infection. Thirty-one of these 124 COPDGene participants were enrolled in the HCV Follow-Up Study. Five of the 31 HCV Follow-Up Study participants did not report a previous diagnosis of HCV. For four of these participants, subsequent clinical HCV testing confirmed HCV infection. Thus, 30/31 Follow-Up Study participants had confirmed HCV diagnoses, supporting the accuracy of the HCV research test results. However, the limited number of participants in the Follow-Up Study precludes an accurate assessment of the false-positive and false-negative rates of the research RNA sequencing evidence for HCV. Most HCV Follow-Up Study participants (29/31) were supportive of returning HCV research results, and most participants found the process for returning HCV results to be informative and not upsetting. Newly diagnosed participants were more likely to be pleased to learn about a potentially curable infection (p = 0.027) and showed a trend toward being more frightened by the potential health risks of HCV (p = 0.11). We conclude that HCV results identified incidentally during transcriptomic research studies can be successfully returned to research study participants with a carefully designed process.

Abstract Expression quantitative trait methylation (eQTM) analysis identifies DNA CpG sites at which methylation is associated with gene expression. The present study describes an eQTM resource of CpG-transcript pairs derived from whole blood DNA methylation and RNA sequencing gene expression data in 2115 Framingham Heart Study participants. We identified 70,047 significant cis CpG-transcript pairs at p

Abstract Detection of viruses by RNA and DNA sequencing has improved the understanding of the human virome. We sought to identify blood viral signatures through secondary use of RNA-sequencing (RNA-seq) data in a large study cohort. The ability to reveal undiagnosed infections with public health implications among study subjects with available sequencing data could enable epidemiologic surveys and may lead to diagnosis and therapeutic interventions, leveraging existing research data in a clinical context. We detected viral RNA in peripheral blood RNA-seq data from a COPD-enriched population of current and former smokers. Correlation between viral detection and both reported infections and relevant disease outcomes was evaluated. We identified Hepatitis C virus RNA in 228 subjects and HIV RNA in 30 subjects. Overall, we observed 31 viral species, including Epstein-Barr virus and Cytomegalovirus. We observed an enrichment of Hepatitis C and HIV infections among subjects reporting liver disease and HIV infections, respectively. Higher interferon expression scores were observed in the subjects with Hepatitis C and HIV infections. Through secondary use of RNA-seq from a cohort of current and former smokers, we detected peripheral blood viral signatures. We identified HIV and Hepatitis C virus (HCV), highlighting potential public health implications for the approach described this study. We observed correlations with reported infections, chronic infection outcomes and the host transcriptomic response, providing evidence to support the validity of the approach.

Abstract Background Chronic obstructive pulmonary disease (COPD) varies significantly in symptomatic and physiologic presentation. Identifying disease subtypes from molecular data, collected from easily accessible blood samples, can help stratify patients and guide disease management and treatment. Methods Blood gene expression measured by RNA-sequencing in the COPDGene Study was analyzed using a network perturbation analysis method. Each COPD sample was compared against a learned reference gene network to determine the part that is deregulated. Gene deregulation values were used to cluster the disease samples. Results The discovery set included 617 former smokers from COPDGene. Four distinct gene network subtypes are identified with significant differences in symptoms, exercise capacity and mortality. These clusters do not necessarily correspond with the levels of lung function impairment and are independently validated in two external cohorts: 769 former smokers from COPDGene and 431 former smokers in the Multi-Ethnic Study of Atherosclerosis (MESA). Additionally, we identify several genes that are significantly deregulated across these subtypes, including DSP and GSTM1, which have been previously associated with COPD through genome-wide association study (GWAS). Conclusions The identified subtypes differ in mortality and in their clinical and functional characteristics, underlining the need for multi-dimensional assessment potentially supplemented by selected markers of gene expression. The subtypes were consistent across cohorts and could be used for new patient stratification and disease prognosis.

Abstract 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.

Abstract Background Interstitial lung abnormalities (ILA) are radiologic findings that may progress to idiopathic pulmonary fibrosis (IPF). Blood gene expression profiles can predict IPF mortality, but whether these same genes associate with ILA and ILA outcomes is unknown. This study evaluated if a previously described blood gene expression profile associated with IPF mortality is associated with ILA and all-cause mortality. Methods In COPDGene and ECLIPSE study participants with visual scoring of ILA and gene expression data, we evaluated the association of a previously described IPF mortality score with ILA and mortality. We also trained a new ILA score, derived using genes from the IPF score, in a subset of COPDGene. We tested the association with ILA and mortality on the remainder of COPDGene and ECLIPSE. Results In 1469 COPDGene (training n = 734; testing n = 735) and 571 ECLIPSE participants, the IPF score was not associated with ILA or mortality. However, an ILA score derived from IPF score genes was associated with ILA (meta-analysis of test datasets OR 1.4 [95% CI: 1.2–1.6]) and mortality (HR 1.25 [95% CI: 1.12–1.41]). Six of the 11 genes in the ILA score had discordant directions of effects compared to the IPF score. The ILA score partially mediated the effects of age on mortality (11.8% proportion mediated). Conclusions An ILA gene expression score, derived from IPF mortality-associated genes, identified genes with concordant and discordant effects on IPF mortality and ILA. These results suggest shared, and unique biologic processes, amongst those with ILA, IPF, aging, and death.

Abstract Background Chronic obstructive pulmonary disease (COPD) and idiopathic pulmonary fibrosis (IPF) are characterized by shared exposures and clinical features, but distinct genetic and pathologic features exist. These features have not been well-studied using large-scale gene expression datasets. We hypothesized that there are divergent gene, pathway, and cellular signatures between COPD and IPF. Methods We performed RNA-sequencing on lung tissues from individuals with IPF (n = 231) and COPD (n = 377) compared to control (n = 267), defined as individuals with normal spirometry. We grouped the overlapping differential expression gene sets based on direction of expression and examined the resultant sets for genes of interest, pathway enrichment, and cell composition. Using gene set variation analysis, we validated the overlap group gene sets in independent COPD and IPF data sets. Results We found 5010 genes differentially expressed between COPD and control, and 11,454 genes differentially expressed between IPF and control (1% false discovery rate). 3846 genes overlapped between IPF and COPD. Several pathways were enriched for genes upregulated in COPD and downregulated in IPF; however, no pathways were enriched for genes downregulated in COPD and upregulated in IPF. There were many myeloid cell genes with increased expression in COPD but decreased in IPF. We found that the genes upregulated in COPD but downregulated in IPF were associated with lower lung function in the independent validation cohorts. Conclusions We identified a divergent gene expression signature between COPD and IPF, with increased expression in COPD and decreased in IPF. This signature is associated with worse lung function in both COPD and IPF.

Abstract Background 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. Methods Participants with COPD (N = 4308) from three studies (COPDGene, ECLIPSE, and SPIROMICS) were analysed. Discovery analyses were performed in COPDGene with replication in SPIROMICS and ECLIPSE. In COPDGene, weight loss was defined as self‐reported unintentional weight loss > 5% in the past year or low BMI (BMI

Abstract The human microbiome has a role in the development of multiple diseases. Individual microbiome profiles are highly personalized, though many species are shared. Understanding the relationship between the human microbiome and disease may inform future individualized treatments. We hypothesize the blood microbiome signature may be a surrogate for some lung microbial characteristics. We sought associations between the blood microbiome signature and lung-relevant host factors. Based on reads not mapped to the human genome, we detected microbial nucleic acids through secondary use of peripheral blood RNA-sequencing from 2,590 current and former smokers with and without chronic obstructive pulmonary disease (COPD) from the COPDGene study. We used the Genome Analysis Toolkit (GATK) microbial pipeline PathSeq to infer microbial profiles. We tested associations between the inferred profiles and lung disease relevant phenotypes and examined links to host gene expression pathways. We replicated our analyses using a second independent set of blood RNA-seq data from 1,065 COPDGene study subjects and performed a meta-analysis across the two studies. The four phyla with highest abundance across all subjects were Proteobacteria, Actinobacteria, Firmicutes and Bacteroidetes. In our meta-analysis, we observed associations (q-value

Abstract Background Survivors of prematurity are at risk for abnormal childhood lung function. Few studies have addressed trajectories of lung function and risk factors for abnormal growth in childhood. This study aims to describe changes in lung function in a contemporary cohort of children born preterm followed longitudinally in pulmonary clinic for post-prematurity respiratory disease and to assess maternal and neonatal risk factors associated with decreased lung function trajectories. Methods Observational cohort of 164 children born preterm ≤ 32 weeks gestation followed in pulmonary clinic at Boston Children’s Hospital with pulmonary function testing. We collected demographics and neonatal history. We used multivariable linear regression to identify the impact of neonatal and maternal risk factors on lung function trajectories in childhood. Results We identified 264 studies from 82 subjects with acceptable longitudinal FEV1 data and 138 studies from 47 subjects with acceptable longitudinal FVC and FEV1/FVC data. FEV1% predicted and FEV1/FVC were reduced compared to childhood norms. Growth in FVC outpaced FEV1, resulting in an FEV1/FVC that declined over time. In multivariable analyses, longer duration of mechanical ventilation was associated with a lower rate of rise in FEV1% predicted and greater decline in FEV1/FVC, and postnatal steroid exposure in the NICU was associated with a lower rate of rise in FEV1 and FVC % predicted. Maternal atopy and asthma were associated with a lower rate of rise in FEV1% predicted. Conclusions Children with post-prematurity respiratory disease demonstrate worsening obstruction in lung function throughout childhood. Neonatal risk factors including exposure to mechanical ventilation and postnatal steroids, as well as maternal atopy and asthma, were associated with diminished rate of rise in lung function. These results may have implications for lung function trajectories into adulthood.

Introduction Bronchodilator response (BDR) is a measurement of acute bronchodilation in response to short-acting β2-agonists, with a heritability between 10 and 40%. Identifying genetic variants associated with BDR may lead to a better understanding of its complex pathophysiology. Methods We performed a genome-wide association study (GWAS) of BDR in six adult cohorts with participants of European ancestry (EA) and African ancestry (AA) including community cohorts and cohorts ascertained on the basis of obstructive pulmonary disease. Validation analysis was carried out in two paediatric asthma cohorts. Results A total of 10 623 EA and 3597 AA participants were included in the analyses. No single nucleotide polymorphism (SNP) was associated with BDR at the conventional genome-wide significance threshold (p

Abstract Introduction Cachexia contributes to increased mortality and reduced quality of life in Chronic Obstructive Pulmonary Disease (COPD) and may be associated with underlying gene expression changes. Our goal was to identify differential gene expression signatures associated with COPD cachexia in current and former smokers. Methods We analyzed whole-blood gene expression data from participants with COPD in a discovery cohort (COPDGene, N = 400) and assessed replication (ECLIPSE, N = 114). To approximate the consensus definition using available criteria, cachexia was defined as weight-loss > 5% in the past 12 months or low body mass index (BMI) ( 5% in the past 12 months or low BMI and 3/5 criteria: decreased muscle strength, anorexia, abnormal biochemistry (anemia or high c-reactive protein (> 5 mg/l)), fatigue, and low FFMI. Differential gene expression was assessed between cachectic and non-cachectic subjects, adjusting for age, sex, white blood cell counts, and technical covariates. Gene set enrichment analysis was performed using MSigDB. Results The prevalence of COPD cachexia was 13.7% in COPDGene and 7.9% in ECLIPSE. Fourteen genes were differentially downregulated in cachectic versus non-cachectic COPD patients in COPDGene (FDR

Abstract Background Previous studies have established a higher prevalence of vitamin D deficiency in patients with COPD, but the relationship between vitamin D levels and COPD exacerbations remains controversial. In addition, the effect of vitamin D levels on imaging characteristics remains mostly unexplored. Using cross-sectional and longitudinal follow up data from the COPDGene Study, we assessed the association between vitamin D levels on respiratory symptoms, exacerbations, and imaging characteristics. We hypothesized that vitamin D deficiency will be associated with worse respiratory-related outcomes. Methods Current and former smokers between ages 45–80 were enrolled the COPDGene Study. Subjects completed questionnaires, spirometry, six-minute walk test, and chest computed tomography scans. A subset of subjects had measurement of serum concentration of 25-hydroxyvitamin D (25(OH)D). Vitamin D deficiency was defined as serum concentration less than 20 ng/mL. Longitudinal follow up was conducted via a web-based or telephone questionnaire. Results Vitamin D levels were measured on 1544 current and former smokers, of which 981 subjects had sufficient vitamin D levels and 563 subjects had vitamin D deficiency. Subjects with vitamin D deficiency were younger with increased likelihood of being African American, being current smokers, having a lower percent predicted FEV1, and having COPD. Vitamin D deficiency was associated with worse quality of life, increased dyspnea, decreased exercise tolerance, and increased frequency of severe exacerbations. Vitamin D deficiency was also associated with increased segmental airway wall thickness on chest CT scans. Conclusion Vitamin D deficiency was associated with increased respiratory symptoms, decreased functional status, increased frequency of severe exacerbations, as well as airway wall thickening on chest CT scans. Further research is needed to determine the potential impact of vitamin D supplementation to improve disease outcomes.

Chronic obstructive pulmonary disease (COPD) is mainly caused by cigarette smoking and characterized by chronic inflammation in vulnerable individuals. However, it is unknown how genetic factors may shape chronic inflammation in COPD. To understand how hedgehog interacting protein, encoded by HHIP gene identified in the genome-wide association study in COPD, plays a role in inflammation, we utilized Hhip+/– mice that present persistent inflammation and emphysema upon aging similar to that observed in human COPD. By performing single-cell RNA sequencing of the whole lung from mice at different ages, we found that Hhip+/– mice developed a cytotoxic immune response with a specific increase in killer cell lectin-like receptor G1–positive CD8+ T cells with upregulated Ifnγ expression recapitulating human COPD. Hhip expression was restricted to a lung fibroblast subpopulation that had increased interaction with CD8+ T lymphocytes in Hhip+/– compared with Hhip+/+ during aging. Hhip-expressing lung fibroblasts had upregulated IL-18 pathway genes in Hhip+/– lung fibroblasts, which was sufficient to drive increased levels of IFN-γ in CD8+ T cells ex vivo. Our finding provides insight into how a common genetic variation contributes to the amplified lymphocytic inflammation in COPD.

In COPD, anaemia is associated with increased morbidity, but the relationship between haemoglobin over its entire observed range and morbidity is poorly understood. Such an understanding could guide future therapeutic targeting of haemoglobin in COPD management. Leveraging the COPDGene study, we conducted a cross-sectional analysis of haemoglobin from COPD participants, examining symptoms, quality of life, functional performance, and acute exacerbations of COPD (AECOPD). Haemoglobin was analysed both as a continuous variable and categorised into anaemia, normal haemoglobin, and polycythaemia groups. Fractional polynomial modelling was used for continuous analyses; categorical models were multivariable linear or negative binomial regressions. Covariates included demographics, comorbidities, emphysema, diffusing capacity, and airflow obstruction. From 2539 participants, 366 (14%) were identified as anaemic and 125 (5%) as polycythaemic. Compared with normal haemoglobin, anaemia was significantly associated with increased symptoms (COPD Assessment Test score: p=0.006, modified Medical Research Council (mMRC) Dyspnoea Score: p=0.001); worse quality of life (St. George's Respiratory Questionnaire (SGRQ) score: p

Abstract Background Cachexia is associated with increased mortality risk among chronic obstructive pulmonary disease (COPD) patients. However, low body mass index (BMI) as opposed to cachexia is often used, particularly when calculating the BODE (BMI, Obstruction, Dyspnea and Exercise) index. For this reason, we examined mortality using a consensus definition and a weight-loss definition of cachexia among COPD cases and compared two new COPD severity indices with BODE. Methods In the current report, the consensus definition for cachexia incorporated weight-loss > 5% in 12-months or low BMI in addition to 3/5 of decreased muscle strength, fatigue, anorexia, low FFMI and inflammation. The weight-loss definition incorporated weight-loss > 5% or weight-loss > 2% (if low BMI) in 12-months. The low BMI component in BODE was replaced with the consensus definition to create the CODE (Consensus cachexia, Obstruction, Dyspnea and Exercise) index and the weight-loss definition to create the WODE (Weight loss, Obstruction, Dyspnea and Exercise) index. Mortality was assessed using Kaplan-Meier survival and Cox Regression. Performance of models was compared using C-statistics. Results Among 1483 COPD cases, the prevalences of cachexia by the consensus and weight-loss definitions were 4.7 and 10.4%, respectively. Cachectic patients had a greater than three-fold increased mortality by either the consensus or the weight-loss definition of cachexia independent of BMI and lung function. The CODE index predicted mortality slightly more accurately than the BODE and WODE indices. Conclusions Cachexia is associated with increased mortality among COPD patients. Monitoring cachexia using weight-loss criteria is relatively simple and predictive of mortality among COPD cases who may be missed if only low BMI is used.

Abstract Background Multiple gene expression studies have been performed separately in peripheral blood, lung, and airway tissues to study COPD. We performed RNA-sequencing gene expression profiling of large-airway epithelium, alveolar macrophage and peripheral blood samples from the same subset of COPD cases and controls from the COPDGene study who underwent bronchoscopy at a single center. Using statistical and gene set enrichment approaches, we sought to improve the understanding of COPD by studying gene sets and pathways across these tissues, beyond the individual genomic determinants. Methods We performed differential expression analysis using RNA-seq data obtained from 63 samples from 21 COPD cases and controls (includes four non-smokers) via the R package DESeq2. We tested associations between gene expression and variables related to lung function, smoking history, and CT scan measures of emphysema and airway disease. We examined the correlation of differential gene expression across the tissues and phenotypes, hypothesizing that this would reveal preserved and private gene expression signatures. We performed gene set enrichment analyses using curated databases and findings from prior COPD studies to provide biological and disease relevance. Results The known smoking-related genes CYP1B1 and AHRR were among the top differential expression results for smoking status in the large-airway epithelium data. We observed a significant overlap of genes primarily across large-airway and macrophage results for smoking and airway disease phenotypes. We did not observe specific genes differentially expressed in all three tissues for any of the phenotypes. However, we did observe hemostasis and immune signaling pathways in the overlaps across all three tissues for emphysema, and amyloid and telomere-related pathways for smoking. In peripheral blood, the emphysema results were enriched for B cell related genes previously identified in lung tissue studies. Conclusions Our integrative analyses across COPD-relevant tissues and prior studies revealed shared and tissue-specific disease biology. These replicated and novel findings in the airway and peripheral blood have highlighted candidate genes and pathways for COPD pathogenesis.

Abstract Background Previous studies have identified lung, sputum or blood transcriptomic biomarkers associated with the severity of airflow limitation in COPD. Yet, it is not clear whether the lung pathobiology is mirrored by these surrogate tissues. The aim of this study was to explore this question. Methods We used Weighted Gene Co-expression Network Analysis (WGCNA) to identify shared pathological mechanisms across four COPD gene-expression datasets: two sets of lung tissues (L1 n = 70; L2 n = 124), and one each of induced sputum (S; n = 121) and peripheral blood (B; n = 121). Results WGCNA analysis identified twenty-one gene co-expression modules in L1. A robust module preservation between the two L datasets was observed (86%), with less preservation in S (33%) and even less in B (23%). Three modules preserved across lung tissues and sputum (not blood) were associated with the severity of airflow limitation. Ontology enrichment analysis showed that these modules included genes related to mitochondrial function, ion-homeostasis, T cells and RNA processing. These findings were largely reproduced using the consensus WGCNA network approach. Conclusions These observations indicate that major differences in lung tissue transcriptomics in patients with COPD are poorly mirrored in sputum and are unrelated to those determined in blood, suggesting that the systemic component in COPD is independently regulated. Finally, the fact that one of the preserved modules associated with FEV1 was enriched in mitochondria-related genes supports a role for mitochondrial dysfunction in the pathobiology of COPD.

MicroRNAs have been independently associated with asthma and COPD; however, it is unclear if microRNA associations will overlap when evaluating retrospective acute exacerbations. Objective: We hypothesized that peripheral blood microRNAs would be associated with retrospective acute asthma exacerbations in a pediatric asthma cohort and that such associations may also be relevant to acute COPD exacerbations. Methods: We conducted small-RNA sequencing on 374 whole-blood samples from children with asthma ages 6–14 years who participated in the Genetics of Asthma in Costa Rica Study (GACRS) and 450 current and former adult smokers with and without COPD who participated in the COPDGene study. Measurements and Main Results: After QC, we had 351 samples and 649 microRNAs for Differential Expression (DE) analysis between the frequent (n = 183) and no or infrequent exacerbation (n = 168) groups in GACRS. Fifteen upregulated miRs had odds ratios (OR) between 1.22 and 1.59 for a doubling of miR counts, while five downregulated miRs had ORs between 0.57 and 0.8. These were assessed for generalization in COPDGene, where three of the upregulated miRs (miR-532-3p, miR-296-5p, and miR-766-3p) and two of the downregulated miRs (miR-7-5p and miR-451b) replicated. Pathway enrichment analysis showed MAPK and PI3K-Akt signaling pathways were strongly enriched for target genes of DE miRNAs and miRNAs generalizing to COPD exacerbations, as well as infection response pathways to various pathogens. Conclusion: miRs (451b; 7-5p; 532-3p; 296-5p and 766-3p) associated with both childhood asthma and adult COPD exacerbations may play a vital role in airflow obstruction and exacerbations and point to shared genomic regulatory machinery underlying exacerbations in both diseases.

Abstract Background Childhood asthma is strongly influenced by genetics and is a risk factor for reduced lung function and chronic obstructive pulmonary disease (COPD) in adults. This study investigates self-reported childhood asthma in adult smokers from the COPDGene Study. We hypothesize that childhood asthma is associated with decreased lung function, increased risk for COPD, and that a genome-wide association study (GWAS) will show association with established asthma variants. Methods We evaluated current and former smokers ages 45–80 of non-Hispanic white (NHW) or African American (AA) race. Childhood asthma was defined by self-report of asthma, diagnosed by a medical professional, with onset at

Rationale Chronic obstructive pulmonary disease (COPD) comprises distinct phenotypes, all characterised by airflow limitation. Objectives We hypothesised that somatotype changes – as a surrogate of adiposity – from early adulthood follow different trajectories to reach distinct phenotypes. Methods Using the validated Stunkard's Pictogram, 356 COPD patients chose the somatotype that best reflects their current body build and those at ages 18, 30, 40 and 50 years. An unbiased group-based trajectory modelling was used to determine somatotype trajectories. We then compared the current COPD-related clinical and phenotypic characteristics of subjects belonging to each trajectory. Measurements and main results At 18 years of age, 88% of the participants described having a lean or medium somatotype (estimated body mass index (BMI) between 19 and 23 kg·m−2) while the other 12% a heavier somatotype (estimated BMI between 25 and 27 kg·m−2). From age 18 onwards, five distinct trajectories were observed. Four of them demonstrating a continuous increase in adiposity throughout adulthood with the exception of one, where the initial increase was followed by loss of adiposity after age 40. Patients with this trajectory were primarily females with low BMI and DLCO (diffusing capacity of the lung for carbon monoxide). A persistently lean trajectory was seen in 14% of the cohort. This group had significantly lower forced expiratory volume in 1 s (FEV1), DLCO, more emphysema and a worse BODE (BMI, airflow obstruction, dyspnoea and exercise capacity) score thus resembling the multiple organ loss of tissue (MOLT) phenotype. Conclusions COPD patients have distinct somatotype trajectories throughout adulthood. Those with the MOLT phenotype maintain a lean trajectory throughout life. Smoking subjects with this lean phenotype in early adulthood deserve particular attention as they seem to develop more severe COPD.

Introduction Family history is a risk factor for chronic obstructive pulmonary disease (COPD). We previously developed a COPD risk score from genome-wide genetic markers (Polygenic Risk Score, PRS). Whether the PRS and family history provide complementary or redundant information for predicting COPD and related outcomes is unknown.Methods We assessed the predictive capacity of family history and PRS on COPD and COPD-related outcomes in non-Hispanic white (NHW) and African American (AA) subjects from COPDGene and ECLIPSE studies. We also performed interaction and mediation analyses.Results In COPDGene, family history and PRS were significantly associated with COPD in a single model (PFamHx

Abstract Background Genome-wide association studies have identified several genetic risk loci for severe chronic obstructive pulmonary disease (COPD) and emphysema. However, these studies do not fully explain disease heritability and in most cases, fail to implicate specific genes. Integrative methods that combine gene expression data with GWAS can provide more power in discovering disease-associated genes and give mechanistic insight into regulated genes. Methods We applied a recently described method that imputes gene expression using reference transcriptome data to genome-wide association studies for two phenotypes (severe COPD and quantitative emphysema) and blood and lung tissue gene expression datasets. We further tested the potential causality of individual genes using multi-variant colocalization. Results We identified seven genes significantly associated with severe COPD, and five genes significantly associated with quantitative emphysema in whole blood or lung. We validated results in independent transcriptome databases and confirmed colocalization signals for PSMA4, EGLN2, WNT3, DCBLD1, and LILRA3. Three of these genes were not located within previously reported GWAS loci for either phenotype. We also identified genetically driven pathways, including those related to immune regulation. Conclusions An integrative analysis of GWAS and gene expression identified novel associations with severe COPD and quantitative emphysema, and also suggested disease-associated genes in known COPD susceptibility loci. Trial registration NCT00608764, Registry: ClinicalTrials.gov, Date of Enrollment of First Participant: November 2007, Date Registered: January 28, 2008 (retrospectively registered); NCT00292552, Registry: ClinicalTrials.gov, Date of Enrollment of First Participant: December 2005, Date Registered: February 14, 2006 (retrospectively registered).

Abstract Background 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. Results We performed expression quantitative trait locus (eQTL) analysis using genome-wide SNP genotyping and gene expression profiling of lung tissue samples from 86 COPD cases and 31 controls, testing for SNPs associated with gene expression levels. These results were integrated with a prior COPD GWAS using an ensemble statistical and network methods approach to identify relevant genes and observe them in the context of overall genetic control of gene expression to highlight co-regulated genes and disease pathways. We identified 250,312 unique SNPs and 4997 genes in the cis(local)-eQTL analysis (5% false discovery rate). The top gene from the integrative analysis was MAPT, a gene recently identified in an independent GWAS of lung function. The genes HNRNPAB and PCBP2 with RNA binding activity and the gene ACVR1B were identified in network communities with validated disease relevance. Conclusions The integration of lung tissue gene expression with genome-wide SNP genotyping and subsequent intersection with prior GWAS and omics studies highlighted candidate genes within COPD loci and in communities harboring known COPD genes. This integration also identified novel disease genes in sub-threshold regions that would otherwise have been missed through GWAS.

Abstract Background Cigarette smoking is the leading modifiable risk factor for disease and death worldwide. Previous studies quantifying gene-level expression have documented the effect of smoking on mRNA levels. Using RNA sequencing, it is possible to analyze the impact of smoking on complex regulatory phenomena (e.g. alternative splicing, differential isoform usage) leading to a more detailed understanding of the biology underlying smoking-related disease. Methods We used whole-blood RNA sequencing to describe gene and exon-level expression differences between 229 current and 286 former smokers in the COPDGene study. We performed differential gene expression and differential exon usage analyses using the voom/limma and DEXseq R packages. Samples from current and former smokers were compared while controlling for age, gender, race, lifetime smoke exposure, cell counts, and technical covariates. Results At an adjusted p-value

Abstract Background Chronic obstructive pulmonary disease (COPD) is currently the third leading cause of death and there is a huge unmet clinical need to identify disease biomarkers in peripheral blood. Compared to gene level differential expression approaches to identify gene signatures, network analyses provide a biologically intuitive approach which leverages the co-expression patterns in the transcriptome to identify modules of co-expressed genes. Methods A weighted gene co-expression network analysis (WGCNA) was applied to peripheral blood transcriptome from 238 COPD subjects to discover co-expressed gene modules. We then determined the relationship between these modules and forced expiratory volume in 1 s (FEV1). In a second, independent cohort of 381 subjects, we determined the preservation of these modules and their relationship with FEV1. For those modules that were significantly related to FEV1, we determined the biological processes as well as the blood cell-specific gene expression that were over-represented using additional external datasets. Results Using WGCNA, we identified 17 modules of co-expressed genes in the discovery cohort. Three of these modules were significantly correlated with FEV1 (FDR

Following publication of the original article [1], the authors provided an updated accession number in the “Availability of data and materials” section of the declarations.

Severe alpha-1 antitrypsin (AAT) deficiency is one of the most common serious genetic diseases in adults of European descent. Individuals with AAT deficiency have a greatly increased risk for emphysema and liver disease. Other manifestations include bronchiectasis, necrotizing panniculitis and granulomatosis with polyangiitis. Despite the frequency and potential severity, AAT deficiency remains under-recognized, and there is often a delay in diagnosis. This review will focus on three recent updates that should serve to encourage testing and diagnosis of AAT deficiency: first, the publication of a randomized clinical trial demonstrating the efficacy of intravenous augmentation therapy in slowing the progression of emphysema in AAT deficiency; second, the mounting evidence showing an increased risk of lung disease in heterozygous PI MZ genotype carriers; last, the recent publication of a clinical practice guideline, outlining diagnosis and management. Though it has been recognized for more than fifty years, AAT deficiency exemplifies the modern paradigm of precision medicine, with a diagnostic test that identifies a genetic subtype of a heterogeneous disease, leading to a targeted treatment.

Background: Electronic nicotine delivery systems (ENDS) are driving an epidemic of vaping. Identifying biomarkers of vaping and dual use (concurrent vaping and smoking) will facilitate studies of the health effects of vaping. To identify putative biomarkers of vaping and dual use, we performed association analysis in an observational cohort of 3,892 COPDGene study participants with blood transcriptomics and/or plasma proteomics data and self-reported current vaping and smoking behavior. Methods: Biomarkers of vaping and dual use were identified through differential expression analysis and related to prospective health events over six years of follow-up. To assess the predictive accuracy of multi-biomarker panels, we constructed predictive models for vaping and smoking categories and prospective health outcomes. Results: We identified three transcriptomic and three proteomic associations with vaping, and 90 transcriptomic and 100 proteomic associations to dual use. Many of these vaping or dual use biomarkers were significantly associated with prospective health outcomes, such as FEV1 decline (three transcripts and 62 proteins), overall mortality (18 transcripts and 73 proteins), respiratory mortality (two transcripts and 23 proteins), respiratory exacerbations (13 proteins) and incident cardiovascular disease (24 proteins). Multimarker models showed good performance discriminating between vaping and smoking behavior and produced informative, modestly powerful predictions of future FEV1 decline, mortality, and respiratory exacerbations. Conclusions: In summary, vaping and dual use are associated with RNA and protein blood-based biomarkers that are also associated with adverse health outcomes.