Precision Approaches to Chronic Obstructive Pulmonary Disease Management.

Chronic obstructive pulmonary disease (COPD) is a leading cause of morbidity and mortality worldwide. COPD heterogeneity has hampered progress in developing pharmacotherapies that affect disease progression. This issue can be addressed by precision medicine approaches, which focus on understanding an individual's disease risk, and tailoring management based on pathobiology, environmental exposures, and psychosocial issues. There is an urgent need to identify COPD patients at high risk for poor outcomes and to understand at a mechanistic level why certain individuals are at high risk. Genetics, omics, and network analytic techniques have started to dissect COPD heterogeneity and identify patients with specific pathobiology. Drug repurposing approaches based on biomarkers of specific inflammatory processes (i.e., type 2 inflammation) are promising. As larger data sets, additional omics, and new analytical approaches become available, there will be enormous opportunities to identify high-risk individuals and treat COPD patients based on their specific pathophysiological derangements. These approaches show great promise for risk stratification, early intervention, drug repurposing, and developing novel therapeutic approaches for COPD.

Polygenic transcriptome risk scores for COPD and lung function improve cross-ethnic portability of prediction in the NHLBI TOPMed program.

While polygenic risk scores (PRSs) enable early identification of genetic risk for chronic obstructive pulmonary disease (COPD), predictive performance is limited when the discovery and target populations are not well matched. Hypothesizing that the biological mechanisms of disease are shared across ancestry groups, we introduce a PrediXcan-derived polygenic transcriptome risk score (PTRS) to improve cross-ethnic portability of risk prediction. We constructed the PTRS using summary statistics from application of PrediXcan on large-scale GWASs of lung function (forced expiratory volume in 1 s [FEV1] and its ratio to forced vital capacity [FEV1/FVC]) in the UK Biobank. We examined prediction performance and cross-ethnic portability of PTRS through smoking-stratified analyses both on 29,381 multi-ethnic participants from TOPMed population/family-based cohorts and on 11,771 multi-ethnic participants from TOPMed COPD-enriched studies. Analyses were carried out for two dichotomous COPD traits (moderate-to-severe and severe COPD) and two quantitative lung function traits (FEV1 and FEV1/FVC). While the proposed PTRS showed weaker associations with disease than PRS for European ancestry, the PTRS showed stronger association with COPD than PRS for African Americans (e.g., odds ratio [OR] = 1.24 [95% confidence interval [CI]: 1.08-1.43] for PTRS versus 1.10 [0.96-1.26] for PRS among heavy smokers with ≥ 40 pack-years of smoking) for moderate-to-severe COPD. Cross-ethnic portability of the PTRS was significantly higher than the PRS (paired t test p < 2.2 × 10-16 with portability gains ranging from 5% to 28%) for both dichotomous COPD traits and across all smoking strata. Our study demonstrates the value of PTRS for improved cross-ethnic portability compared to PRS in predicting COPD risk.

Blood-based Transcriptomic and Proteomic Biomarkers of Emphysema.

Rationale: Emphysema is a chronic obstructive pulmonary disease phenotype with important prognostic implications. Identifying blood-based biomarkers of emphysema will facilitate early diagnosis and development of targeted therapies. Objectives: To discover blood omics biomarkers for chest computed tomography-quantified emphysema and develop predictive biomarker panels. Methods: Emphysema blood biomarker discovery was performed using differential gene expression, alternative splicing, and protein association analyses in a training sample of 2,370 COPDGene participants with available blood RNA sequencing, plasma proteomics, and clinical data. Internal validation was conducted in a COPDGene testing sample (n = 1,016), and external validation was done in the ECLIPSE study (n = 526). Because low body mass index (BMI) and emphysema often co-occur, we performed a mediation analysis to quantify the effect of BMI on gene and protein associations with emphysema. Elastic net models with bootstrapping were also developed in the training sample sequentially using clinical, blood cell proportions, RNA-sequencing, and proteomic biomarkers to predict quantitative emphysema. Model accuracy was assessed by the area under the receiver operating characteristic curves for subjects stratified into tertiles of emphysema severity. Measurements and Main Results: Totals of 3,829 genes, 942 isoforms, 260 exons, and 714 proteins were significantly associated with emphysema (false discovery rate, 5%) and yielded 11 biological pathways. Seventy-four percent of these genes and 62% of these proteins showed mediation by BMI. Our prediction models demonstrated reasonable predictive performance in both COPDGene and ECLIPSE. The highest-performing model used clinical, blood cell, and protein data (area under the receiver operating characteristic curve in COPDGene testing, 0.90; 95% confidence interval, 0.85-0.90). Conclusions: Blood transcriptome and proteome-wide analyses revealed key biological pathways of emphysema and enhanced the prediction of emphysema.

Genetically Predicted Body Mass Index and Mortality in COPD.

RATIONALE: BMI is associated with COPD mortality, but the underlying mechanisms are unclear. The effect of genetic variants aggregated into a polygenic score may elucidate causal mechanisms and predict risk. OBJECTIVES: To examine the associations of genetically predicted BMI with all-cause and cause-specific mortality in COPD. METHODS: We developed a polygenic score for BMI (PGSBMI) and tested for associations of the PGSBMI with all-cause, respiratory, and cardiovascular mortality in participants with COPD from the COPDGene, ECLIPSE, and Framingham Heart studies. We calculated the difference between measured BMI and PGS-predicted BMI (BMIdiff) and categorized participants into groups of discordantly low (BMIdiff < 20th percentile), concordant (BMIdiff between 20th - 80th percentile), and discordantly high (BMIdiff > 80th percentile) BMI. We applied Cox models, examined potential non-linear associations of the PGSBMI and BMIdiff with mortality, and summarized results with meta-analysis. MEASUREMENTS AND MAIN RESULTS: We observed significant non-linear associations of measured BMI and BMIdiff, but not PGSBMI, with all-cause mortality. In meta-analyses, a one standard deviation increase in the PGSBMI was associated with an increased hazard for cardiovascular mortality (HR=1.29, 95% CI=1.12-1.49), but not with respiratory or all-cause mortality. Compared to participants with concordant measured and genetically predicted BMI, those with discordantly low BMI had higher mortality risk for all-cause (HR=1.57, CI=1.41-1.74) and respiratory death (HR=2.01, CI=1.61-2.51). CONCLUSIONS: In people with COPD, higher genetically predicted BMI is associated with higher cardiovascular mortality but not respiratory mortality. Individuals with discordantly low BMI have higher all-cause and respiratory mortality compared to those with concordant BMI.

A robust and adaptive framework for interaction testing in quantitative traits between multiple genetic loci and exposure variables.

The identification and understanding of gene-environment interactions can provide insights into the pathways and mechanisms underlying complex diseases. However, testing for gene-environment interaction remains a challenge since a.) statistical power is often limited and b.) modeling of environmental effects is nontrivial and such model misspecifications can lead to false positive interaction findings. To address the lack of statistical power, recent methods aim to identify interactions on an aggregated level using, for example, polygenic risk scores. While this strategy can increase the power to detect interactions, identifying contributing genes and pathways is difficult based on these relatively global results. Here, we propose RITSS (Robust Interaction Testing using Sample Splitting), a gene-environment interaction testing framework for quantitative traits that is based on sample splitting and robust test statistics. RITSS can incorporate sets of genetic variants and/or multiple environmental factors. Based on the user's choice of statistical/machine learning approaches, a screening step selects and combines potential interactions into scores with improved interpretability. In the testing step, the application of robust statistics minimizes the susceptibility to main effect misspecifications. Using extensive simulation studies, we demonstrate that RITSS controls the type 1 error rate in a wide range of scenarios, and we show how the screening strategy influences statistical power. In an application to lung function phenotypes and human height in the UK Biobank, RITSS identified highly significant interactions based on subcomponents of genetic risk scores. While the contributing single variant interaction signals are weak, our results indicate interaction patterns that result in strong aggregated effects, providing potential insights into underlying gene-environment interaction mechanisms.

Association of clonal hematopoiesis with chronic obstructive pulmonary disease.

Chronic obstructive pulmonary disease (COPD) is associated with age and smoking, but other determinants of the disease are incompletely understood. Clonal hematopoiesis of indeterminate potential (CHIP) is a common, age-related state in which somatic mutations in clonal blood populations induce aberrant inflammatory responses. Patients with CHIP have an elevated risk for cardiovascular disease, but the association of CHIP with COPD remains unclear. We analyzed whole-genome sequencing and whole-exome sequencing data to detect CHIP in 48 835 patients, of whom 8444 had moderate to very severe COPD, from four separate cohorts with COPD phenotyping and smoking history. We measured emphysema in murine models in which Tet2 was deleted in hematopoietic cells. In the COPDGene cohort, individuals with CHIP had risks of moderate-to-severe, severe, or very severe COPD that were 1.6 (adjusted 95% confidence interval [CI], 1.1-2.2) and 2.2 (adjusted 95% CI, 1.5-3.2) times greater than those for noncarriers. These findings were consistently observed in three additional cohorts and meta-analyses of all patients. CHIP was also associated with decreased FEV1% predicted in the COPDGene cohort (mean between-group differences, -5.7%; adjusted 95% CI, -8.8% to -2.6%), a finding replicated in additional cohorts. Smoke exposure was associated with a small but significant increased risk of having CHIP (odds ratio, 1.03 per 10 pack-years; 95% CI, 1.01-1.05 per 10 pack-years) in the meta-analysis of all patients. Inactivation of Tet2 in mouse hematopoietic cells exacerbated the development of emphysema and inflammation in models of cigarette smoke exposure. Somatic mutations in blood cells are associated with the development and severity of COPD, independent of age and cumulative smoke exposure.

A Polygenic Risk Score for Idiopathic Pulmonary Fibrosis and Interstitial Lung Abnormalities.

Rationale: In addition to rare genetic variants and the MUC5B locus, common genetic variants contribute to idiopathic pulmonary fibrosis (IPF) risk. The predictive power of common variants outside the MUC5B locus for IPF and interstitial lung abnormalities (ILAs) is unknown. Objectives: We tested the predictive value of IPF polygenic risk scores (PRSs) with and without the MUC5B region on IPF, ILA, and ILA progression. Methods: We developed PRSs that included (PRS-M5B) and excluded (PRS-NO-M5B) the MUC5B region (500-kb window around rs35705950-T) using an IPF genome-wide association study. We assessed PRS associations with area under the receiver operating characteristic curve (AUC) metrics for IPF, ILA, and ILA progression. Measurements and Main Results: We included 14,650 participants (1,970 IPF; 1,068 ILA) from six multi-ancestry population-based and case-control cohorts. In cases excluded from genome-wide association study, the PRS-M5B (odds ratio [OR] per SD of the score, 3.1; P = 7.1 × 10-95) and PRS-NO-M5B (OR per SD, 2.8; P = 2.5 × 10-87) were associated with IPF. Participants in the top PRS-NO-M5B quintile had ∼sevenfold odds for IPF compared with those in the first quintile. A clinical model predicted IPF (AUC, 0.61); rs35705950-T and PRS-NO-M5B demonstrated higher AUCs (0.73 and 0.7, respectively), and adding both genetic predictors to a clinical model yielded the highest performance (AUC, 0.81). The PRS-NO-M5B was associated with ILA (OR, 1.25) and ILA progression (OR, 1.16) in European ancestry participants. Conclusions: A common genetic variant risk score complements the MUC5B variant to identify individuals at high risk of interstitial lung abnormalities and pulmonary fibrosis.

Polygenic risk scores identify heterogeneity in asthma and chronic obstructive pulmonary disease.

BACKGROUND: Asthma and chronic obstructive pulmonary disease (COPD) have distinct and overlapping genetic and clinical features. OBJECTIVE: We sought to test the hypothesis that polygenic risk scores (PRSs) for asthma (PRSAsthma) and spirometry (FEV1 and FEV1/forced vital capacity; PRSspiro) would demonstrate differential associations with asthma, COPD, and asthma-COPD overlap (ACO). METHODS: We developed and tested 2 asthma PRSs and applied the higher performing PRSAsthma and a previously published PRSspiro to research (Genetic Epidemiology of COPD study and Childhood Asthma Management Program, with spirometry) and electronic health record-based (Mass General Brigham Biobank and Genetic Epidemiology Research on Adult Health and Aging [GERA]) studies. We assessed the association of PRSs with COPD and asthma using modified random-effects and binary-effects meta-analyses, and ACO and asthma exacerbations in specific cohorts. Models were adjusted for confounders and genetic ancestry. RESULTS: In meta-analyses of 102,477 participants, the PRSAsthma (odds ratio [OR] per SD, 1.16 [95% CI, 1.14-1.19]) and PRSspiro (OR per SD, 1.19 [95% CI, 1.17-1.22]) both predicted asthma, whereas the PRSspiro predicted COPD (OR per SD, 1.25 [95% CI, 1.21-1.30]). However, results differed by cohort. The PRSspiro was not associated with COPD in GERA and Mass General Brigham Biobank. In the Genetic Epidemiology of COPD study, the PRSAsthma (OR per SD: Whites, 1.3; African Americans, 1.2) and PRSspiro (OR per SD: Whites, 2.2; African Americans, 1.6) were both associated with ACO. In GERA, the PRSAsthma was associated with asthma exacerbations (OR, 1.18) in Whites; the PRSspiro was associated with asthma exacerbations in White, LatinX, and East Asian participants. CONCLUSIONS: PRSs for asthma and spirometry are both associated with ACO and asthma exacerbations. Genetic prediction performance differs in research versus electronic health record-based cohorts.

Prevalence and mortality associations of interstitial lung abnormalities in rheumatoid arthritis within a multicentre prospective cohort of smokers.

OBJECTIVE: To investigate the prevalence and mortality impact of interstitial lung abnormalities (ILAs) in RA and non-RA comparators. METHODS: We analysed associations between ILAs, RA, and mortality in COPDGene, a multicentre prospective cohort study of current and past smokers, excluding known interstitial lung disease (ILD) or bronchiectasis. All participants had research chest high-resolution CT (HRCT) reviewed by a sequential reading method to classify ILA as present, indeterminate or absent. RA cases were identified by self-report RA and DMARD use; non-RA comparators had neither an RA diagnosis nor used DMARDs. We examined the association and mortality risk of RA and ILA using multivariable logistic regression and Cox regression. RESULTS: We identified 83 RA cases and 8725 non-RA comparators with HRCT performed for research purposes. ILA prevalence was 16.9% in RA cases and 5.0% in non-RA comparators. After adjusting for potential confounders, including genetics, current/past smoking and other lifestyle factors, ILAs were more common among those with RA compared with non-RA [odds ratio 4.76 (95% CI 2.54, 8.92)]. RA with ILAs or indeterminate for ILAs was associated with higher all-cause mortality compared with non-RA without ILAs [hazard ratio (HR) 3.16 (95% CI 2.11, 4.74)] and RA cases without ILA [HR 3.02 (95% CI 1.36, 6.75)]. CONCLUSIONS: In this cohort of smokers, RA was associated with ILAs and this persisted after adjustment for current/past smoking and genetic/lifestyle risk factors. RA with ILAs in smokers had a 3-fold increased all-cause mortality, emphasizing the importance of further screening and treatment strategies for preclinical ILD in RA.

Rheumatoid arthritis, quantitative parenchymal lung features, and mortality among smokers.

OBJECTIVES: There have been limited investigations of the prevalence and mortality impact of quantitative computed tomography (QCT) parenchymal lung features in rheumatoid arthritis (RA). We examined the cross-sectional prevalence and mortality associations of QCT features, comparing RA and non-RA participants. METHODS: We identified participants with and without RA in COPDGene, a multicentre cohort study of current or former smokers. Using a k-nearest neighbor quantifier, high resolution CT chest scans were scored for percentage of normal lung, interstitial changes, and emphysema. We examined associations between QCT features and RA using multivariable linear regression. After dichotomizing participants at the 75th percentile for each QCT feature among non-RA participants, we investigated mortality associations by RA/non-RA status and quartile 4 vs quartiles 1-3 of QCT features using Cox regression. We assessed for statistical interactions between RA and QCT features. RESULTS: We identified 82 RA cases and 8820 non-RA comparators. In multivariable linear regression, RA was associated with higher percentage of interstitial changes (ß = 1.7 ± 0.5, p= 0.0008) but not emphysema (ß = 1.3 ± 1.7, p= 0.44). Participants with RA and >75th percentile of emphysema had significantly higher mortality than non-RA participants (HR 5.86, 95%CI 3.75-9.13) as well as RA participants (HR 5.56, 95%CI 2.71-11.38) with ≤75th percentile of emphysema. There were statistical interactions between RA and emphysema for mortality (multiplicative p= 0.014; attributable proportion 0.53, 95%CI 0.30-0.70). CONCLUSIONS: Using machine learning-derived QCT data in a cohort of smokers, RA was associated with higher percentage of interstitial changes. The combination of RA and emphysema conferred >5-fold higher mortality.

Lower myostatin and higher MUC1 levels are associated with better response to mepolizumab and omalizumab in asthma: a protein-protein interaction analyses.

INTRODUCTION: Biomarkers are needed to inform the choice of biologic therapy in patients with asthma given the increasing number of biologics. We aimed to identify proteins associated with response to omalizumab and mepolizumab. METHODS: Aptamer-based proteomic profiling (SomaScan) was used to assess 1437 proteins from 51 patients with moderate to severe asthma who received omalizumab (n = 29) or mepolizumab (n = 22). Response was defined as the change in asthma-related exacerbations in the 12 months following therapy initiation. All models were adjusted for age, sex, and pre-treatment exacerbation rate. Additionally, body mass index was included in the omalizumab model and eosinophil count in the mepolizumab model. We evaluated the association between molecular signatures and response using negative binomial regression correcting for the false discovery rate (FDR) and gene set enrichment analyses (GSEA) to identify associated pathways. RESULTS: Over two-thirds of patients were female. The average age for omalizumab patients was 42 years and 57 years for mepolizumab. At baseline, the average exacerbation rate was 1.5/year for omalizumab and 2.4/year for mepolizumab. Lower levels of LOXL2 (unadjusted p: 1.93 × 10E-05, FDR-corrected: 0.028) and myostatin (unadjusted: 3.87 × 10E-05, FDR-corrected: 0.028) were associated with better response to mepolizumab. Higher levels of CD9 antigen (unadjusted: 5.30 × 10E-07, FDR-corrected: 0.0006) and MUC1 (unadjusted: 1.15 × 10E-06, FDR-corrected: 0.0006) were associated with better response to omalizumab, and LTB4R (unadjusted: 1.12 × 10E-06, FDR-corrected: 0.0006) with worse response. Protein-protein interaction network modeling showed an enrichment of the TNF- and NF-kB signaling pathways for patients treated with mepolizumab and multiple pathways involving MAPK, including the FcER1 pathway, for patients treated with omalizumab. CONCLUSIONS: This study provides novel fundamental data on proteins associated with response to mepolizumab or omalizumab in severe asthma and warrants further validation as potential biomarkers for therapy selection.

Alpha-1 Antitrypsin MZ Heterozygosity Is an Endotype of Chronic Obstructive Pulmonary Disease.

Rationale: Multiple studies have demonstrated an increased risk of chronic obstructive pulmonary disease (COPD) in heterozygous carriers of the AAT (alpha-1 antitrypsin) Z allele. However, it is not known if MZ subjects with COPD are phenotypically different from noncarriers (MM genotype) with COPD. Objectives: To assess if MZ subjects with COPD have different clinical features compared with MM subjects with COPD. Methods: Genotypes of SERPINA1 were ascertained by using whole-genome sequencing data in three independent studies. We compared outcomes between MM subjects with COPD and MZ subjects with COPD in each study and combined the results in a meta-analysis. We performed longitudinal and survival analyses to compare outcomes in MM and MZ subjects with COPD over time. Measurements and Main Results: We included 290 MZ subjects with COPD and 6,184 MM subjects with COPD across the three studies. MZ subjects had a lower FEV1% predicted and greater quantitative emphysema on chest computed tomography scans compared with MM subjects. In a meta-analysis, the FEV1 was 3.9% lower (95% confidence interval [CI], -6.55% to -1.26%) and emphysema (the percentage of lung attenuation areas <-950 HU) was 4.14% greater (95% CI, 1.44% to 6.84%) in MZ subjects. We found one gene, PGF (placental growth factor), to be differentially expressed in lung tissue from one study between MZ subjects and MM subjects. Conclusions: Carriers of the AAT Z allele (those who were MZ heterozygous) with COPD had lower lung function and more emphysema than MM subjects with COPD. Taken with the subtle differences in gene expression between the two groups, our findings suggest that MZ subjects represent an endotype of COPD.

Development of a Blood-based Transcriptional Risk Score for Chronic Obstructive Pulmonary Disease.

Rationale: The ability of peripheral blood biomarkers to assess chronic obstructive pulmonary disease (COPD) risk and progression is unknown. Genetics and gene expression may capture important aspects of COPD-related biology that predict disease activity. Objectives: Develop a transcriptional risk score (TRS) for COPD and assess the contribution of the TRS and a polygenic risk score (PRS) for disease susceptibility and progression. Methods: We randomly split 2,569 COPDGene (Genetic Epidemiology of COPD) participants with whole-blood RNA sequencing into training (n = 1,945) and testing (n = 624) samples and used 468 ECLIPSE (Evaluation of COPD Longitudinally to Identify Predictive Surrogate End-points) COPD cases with microarray data for replication. We developed a TRS using penalized regression (least absolute shrinkage and selection operator) to model FEV1/FVC and studied the predictive value of TRS for COPD (Global Initiative for Chronic Obstructive Lung Disease 2-4), prospective FEV1 change (ml/yr), and additional COPD-related traits. We adjusted for potential confounders, including age and smoking. We evaluated the predictive performance of the TRS in the context of a previously derived PRS and clinical factors. Measurements and Main Results: The TRS included 147 transcripts and was associated with COPD (odds ratio, 3.3; 95% confidence interval [CI], 2.4-4.5; P < 0.001), FEV1 change (ß, -17 ml/yr; 95% CI, -28 to -6.6; P = 0.002), and other COPD-related traits. In ECLIPSE cases, we replicated the association with FEV1 change (ß, -8.2; 95% CI, -15 to -1; P = 0.025) and the majority of other COPD-related traits. Models including PRS, TRS, and clinical factors were more predictive of COPD (area under the receiver operator characteristic curve, 0.84) and annualized FEV1 change compared with models with one risk score or clinical factors alone. Conclusions: Blood transcriptomics can improve prediction of COPD and lung function decline when added to a PRS and clinical risk factors.

Associations of Monocyte Count and Other Immune Cell Types with Interstitial Lung Abnormalities.

Rationale: Higher blood monocyte counts are associated with worse survival in adults with clinically diagnosed pulmonary fibrosis. Their association with the development and progression of interstitial lung abnormalities (ILA) in humans is unknown. Objectives: We evaluated the associations of blood monocyte count, and other immune cell types, with ILA, high-attenuation areas, and FVC in four independent cohorts. Methods: We included participants with measured monocyte counts and computed tomographic (CT) imaging enrolled in MESA (Multi-Ethnic Study of Atherosclerosis, n = 484), AGES-Reykjavik (Age/Gene Environment Susceptibility Study, n = 3,547), COPDGene (Genetic Epidemiology of COPD, n = 2,719), and the ECLIPSE (Evaluation of COPD Longitudinally to Identify Predictive Surrogate End-points, n = 646). Measurements and Main Results: After adjustment for covariates, a 1-SD increment in blood monocyte count was associated with ILA in MESA (odds ratio [OR], 1.3; 95% confidence interval [CI], 1.0-1.8), AGES-Reykjavik (OR, 1.2; 95% CI, 1.1-1.3), COPDGene (OR, 1.3; 95% CI, 1.2-1.4), and ECLIPSE (OR, 1.2; 95% CI, 1.0-1.4). A higher monocyte count was associated with ILA progression over 5 years in AGES-Reykjavik (OR, 1.2; 95% CI, 1.0-1.3). Compared with participants without ILA, there was a higher percentage of activated monocytes among those with ILA in MESA. Higher monocyte count was associated with greater high-attenuation areas in MESA and lower FVC in MESA and COPDGene. Associations of other immune cell types were less consistent. Conclusions: Higher blood monocyte counts were associated with the presence and progression of interstitial lung abnormalities and lower FVC.

Blood RNA sequencing shows overlapping gene expression across COPD phenotype domains.

RATIONALE: COPD can be assessed using multidimensional grading systems with components from three domains: pulmonary function tests, symptoms and systemic features. Clinically, measures may be used interchangeably, though it is not known if they share similar pathobiology. OBJECTIVE: To use RNA sequencing (RNA-seq) to determine if there is an overlap in the underlying biological mechanisms and consequences driving different components of the multidimensional grading systems. METHODS: Whole blood was collected for RNA-seq from current and former smokers in the Genetic Epidemiology of COPD study. We tested the overlap in gene expression and biological pathways associated with case-control status and quantitative COPD phenotypes within and between the three domains. RESULTS: In 2647 subjects, there were 3030 genes differentially expressed in any of the three domains or case-control status. There were five genes that overlapped between the three domains and case-control status, including G protein-coupled receptor 15(GPR15), sestrin 1 (SESN1) and interferon-induced guanylate-binding protein 1 (GBP1), which were associated with longitudinal decline in FEV1. The overlap between the three domains was enriched for pathways related to cellular components. CONCLUSIONS: We identified gene sets and pathways that overlap between 12 COPD-related phenotypes and case-control status. There were no pathways represented in the overlap between the three domains and case-control status, but we identified multiple genes that demonstrated a consistent pattern of expression across several of the phenotypes. Patterns of gene expression correlation were generally similar to the correlation of clinical phenotypes in the PFT and symptom domains but not the systemic features.

A polygenic risk score and age of diagnosis of COPD.

BACKGROUND: Genetic susceptibility may be associated with earlier onset of chronic obstructive pulmonary disease (COPD). We hypothesised that a polygenic risk score (PRS) for COPD would be associated with earlier age of diagnosis of COPD. METHODS: In 6647 non-Hispanic White (NHW) and 2464 African American (AA) participants from COPDGene, and 6812 participants from the Framingham Heart Study (FHS), we tested the relationship of the PRS and age of COPD diagnosis. Age at diagnosis was determined by: 1) self-reported age at COPD diagnosis or 2) age at visits when moderate-to-severe airflow limitation (Global Initiative for Chronic Obstructive Lung Disease (GOLD) grade 2-4) was observed on spirometry. We used Cox regression to examine the overall and time-dependent effects of the PRS on incident COPD. In the COPDGene study, we also examined the PRS's predictive value for COPD at age <50 years (COPD50) using logistic regression and area under the curve (AUC) analyses, with and without the addition of other risk factors present at early life (e.g. childhood asthma). RESULTS: In Cox models, the PRS demonstrated age-dependent associations with incident COPD, with larger effects at younger ages in both cohorts. The PRS was associated with COPD50 (OR 1.55 (95% CI 1.41-1.71) for NHW, OR 1.23 (95% CI 1.05-1.43) for AA and OR 2.47 (95% CI 2.12-2.88) for FHS participants). In COPDGene, adding the PRS to known early-life risk factors improved prediction of COPD50 in NHW (AUC 0.69 versus 0.74; p<0.0001) and AA (AUC 0.61 versus 0.64; p=0.04) participants. CONCLUSIONS: A COPD PRS is associated with earlier age of diagnosis of COPD and retains predictive value when added to known early-life risk factors.

Genetic Variation in the Mitochondrial Glycerol-3-Phosphate Acyltransferase Is Associated With Liver Injury.

BACKGROUND AND AIMS: Most of the genetic basis of chronic liver disease remains undiscovered. APPROACH AND RESULTS: To identify genetic loci that modulate the risk of liver injury, we performed genome-wide association studies on circulating levels of alanine aminotransferase (ALT), aspartate aminotransferase (AST), alkaline phosphatase (ALP), and total bilirubin across 312,671 White British participants in the UK Biobank. We focused on variants associated with elevations in all four liver biochemistries at genome-wide significance (P < 5 × 10-8 ) and that replicated using Mass General Brigham Biobank in 19,323 European ancestry individuals. We identified a genetic locus in mitochondrial glycerol-3-phosphate acyltransferase (GPAM rs10787429) associated with increased levels of ALT (P = 1.4 × 10-30 ), AST (P = 3.6 × 10-10 ), ALP (P = 9.5 × 10-30 ), and total bilirubin (P = 2.9 × 10-12 ). This common genetic variant was also associated with an allele dose-dependent risk of alcohol-associated liver disease (odd ratio [OR] = 1.34, P = 2.6 × 10-5 ) and fatty liver disease (OR = 1.18, P = 5.8 × 10-4 ) by International Classification of Diseases, 10th Revision codes. We identified significant interactions between GPAM rs10787429 and elevated body mass index in association with ALT and AST (P = 7.1 × 10-9 and 3.95 × 10-8 , respectively), as well as between GPAM rs10787429 and weekly alcohol consumption in association with ALT, AST, and alcohol-associated liver disease (P = 4.0 × 10-2 , 1.6 × 10-2 , and 1.3 × 10-2 , respectively). Unlike previously described genetic variants that are associated with an increased risk of liver injury but confer a protective effect on circulating lipids, GPAM rs10787429 was associated with an increase in total cholesterol (P = 2.0 × 10-17 ), LDL cholesterol (P = 2.0 × 10-10 ), and HDL cholesterol (P = 6.6 × 10-37 ). Single-cell RNA-sequencing data demonstrated hepatocyte-predominant expression of GPAM in cells that co-express genes related to VLDL production (P = 9.4 × 10-103 ). CONCLUSIONS: Genetic variation in GPAM is associated with susceptibility to liver injury. GPAM may represent a therapeutic target in chronic liver disease.

Associations of the MUC5B promoter variant with timing of interstitial lung disease and rheumatoid arthritis onset.

OBJECTIVES: To investigate the associations of the common MUC5B promoter variant with timing of RA-associated interstitial lung disease (RA-ILD) and RA onset. METHODS: We identified patients with RA meeting 2010 ACR/EULAR criteria and available genotype information in the Mass General Brigham Biobank, a multihospital biospecimen and clinical data collection research study. We determined RA-ILD presence by reviewing all RA patients who had CT imaging, lung biopsy or autopsy results. We determined the dates of RA and RA-ILD diagnoses by manual records review. We examined the associations of the MUC5B promoter variant (G>T at rs35705950) with RA-ILD, RA-ILD occurring before or within 2 years of RA diagnosis and RA diagnosis at age >55 years. We used multivariable logistic regression to estimate odds ratios (ORs) for each outcome by MUC5B promoter variant status, adjusting for potential confounders including genetic ancestry and smoking. RESULTS: We identified 1005 RA patients with available genotype data for rs35705950 (mean age 45 years, 79% female, 81% European ancestry). The MUC5B promoter variant was present in 155 (15.4%) and was associated with RA-ILD [multivariable OR 3.34 (95% CI 1.97, 5.60)], RA-ILD before or within 2 years of RA diagnosis [OR 4.01 (95% CI 1.78, 8.80)] and RA onset after age 55 years [OR 1.52 (95% CI 1.08, 2.12)]. CONCLUSIONS: The common MUC5B promoter variant was associated with RA-ILD onset earlier in the RA disease course and older age of RA onset. These findings suggest that the MUC5B promoter variant may impact RA-ILD risk early in the RA disease course, particularly in patients with older-onset RA.

Blood gene expression risk profiles and interstitial lung abnormalities: COPDGene and ECLIPSE cohort studies.

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.

Lung tissue shows divergent gene expression between chronic obstructive pulmonary disease and idiopathic pulmonary fibrosis.

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.