Aberrant activation of TCL1A promotes stem cell expansion in clonal haematopoiesis.

Mutations in a diverse set of driver genes increase the fitness of haematopoietic stem cells (HSCs), leading to clonal haematopoiesis1. These lesions are precursors for blood cancers2-6, but the basis of their fitness advantage remains largely unknown, partly owing to a paucity of large cohorts in which the clonal expansion rate has been assessed by longitudinal sampling. Here, to circumvent this limitation, we developed a method to infer the expansion rate from data from a single time point. We applied this method to 5,071 people with clonal haematopoiesis. A genome-wide association study revealed that a common inherited polymorphism in the TCL1A promoter was associated with a slower expansion rate in clonal haematopoiesis overall, but the effect varied by driver gene. Those carrying this protective allele exhibited markedly reduced growth rates or prevalence of clones with driver mutations in TET2, ASXL1, SF3B1 and SRSF2, but this effect was not seen in clones with driver mutations in DNMT3A. TCL1A was not expressed in normal or DNMT3A-mutated HSCs, but the introduction of mutations in TET2 or ASXL1 led to the expression of TCL1A protein and the expansion of HSCs in vitro. The protective allele restricted TCL1A expression and expansion of mutant HSCs, as did experimental knockdown of TCL1A expression. Forced expression of TCL1A promoted the expansion of human HSCs in vitro and mouse HSCs in vivo. Our results indicate that the fitness advantage of several commonly mutated driver genes in clonal haematopoiesis may be mediated by TCL1A activation.

Implications of Race Adjustment in Lung-Function Equations.

BACKGROUND: Adjustment for race is discouraged in lung-function testing, but the implications of adopting race-neutral equations have not been comprehensively quantified. METHODS: We obtained longitudinal data from 369,077 participants in the National Health and Nutrition Examination Survey, U.K. Biobank, the Multi-Ethnic Study of Atherosclerosis, and the Organ Procurement and Transplantation Network. Using these data, we compared the race-based 2012 Global Lung Function Initiative (GLI-2012) equations with race-neutral equations introduced in 2022 (GLI-Global). Evaluated outcomes included national projections of clinical, occupational, and financial reclassifications; individual lung-allocation scores for transplantation priority; and concordance statistics (C statistics) for clinical prediction tasks. RESULTS: Among the 249 million persons in the United States between 6 and 79 years of age who are able to produce high-quality spirometric results, the use of GLI-Global equations may reclassify ventilatory impairment for 12.5 million persons, medical impairment ratings for 8.16 million, occupational eligibility for 2.28 million, grading of chronic obstructive pulmonary disease for 2.05 million, and military disability compensation for 413,000. These potential changes differed according to race; for example, classifications of nonobstructive ventilatory impairment may change dramatically, increasing 141% (95% confidence interval [CI], 113 to 169) among Black persons and decreasing 69% (95% CI, 63 to 74) among White persons. Annual disability payments may increase by more than $1 billion among Black veterans and decrease by $0.5 billion among White veterans. GLI-2012 and GLI-Global equations had similar discriminative accuracy with regard to respiratory symptoms, health care utilization, new-onset disease, death from any cause, death related to respiratory disease, and death among persons on a transplant waiting list, with differences in C statistics ranging from -0.008 to 0.011. CONCLUSIONS: The use of race-based and race-neutral equations generated similarly accurate predictions of respiratory outcomes but assigned different disease classifications, occupational eligibility, and disability compensation for millions of persons, with effects diverging according to race. (Funded by the National Heart Lung and Blood Institute and the National Institute of Environmental Health Sciences.).

Whole genome sequencing based analysis of inflammation biomarkers in the Trans-Omics for Precision Medicine (TOPMed) consortium.

Inflammation biomarkers can provide valuable insight into the role of inflammatory processes in many diseases and conditions. Sequencing based analyses of such biomarkers can also serve as an exemplar of the genetic architecture of quantitative traits. To evaluate the biological insight, which can be provided by a multi-ancestry, whole-genome based association study, we performed a comprehensive analysis of 21 inflammation biomarkers from up to 38 465 individuals with whole-genome sequencing from the Trans-Omics for Precision Medicine (TOPMed) program (with varying sample size by trait, where the minimum sample size was n = 737 for MMP-1). We identified 22 distinct single-variant associations across 6 traits-E-selectin, intercellular adhesion molecule 1, interleukin-6, lipoprotein-associated phospholipase A2 activity and mass, and P-selectin-that remained significant after conditioning on previously identified associations for these inflammatory biomarkers. We further expanded upon known biomarker associations by pairing the single-variant analysis with a rare variant set-based analysis that further identified 19 significant rare variant set-based associations with 5 traits. These signals were distinct from both significant single variant association signals within TOPMed and genetic signals observed in prior studies, demonstrating the complementary value of performing both single and rare variant analyses when analyzing quantitative traits. We also confirm several previously reported signals from semi-quantitative proteomics platforms. Many of these signals demonstrate the extensive allelic heterogeneity and ancestry-differentiated variant-trait associations common for inflammation biomarkers, a characteristic we hypothesize will be increasingly observed with well-powered, large-scale analyses of complex traits.

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.

Chitinase 3-like-1 Inhibits Innate Antitumor and Tissue Remodeling Immune Responses by Regulating CD47-SIRPα- and CD24-Siglec10-Mediated Phagocytosis.

Innate immune responses such as phagocytosis are critically linked to the generation of adaptive immune responses against the neoantigens in cancer and the efferocytosis that is essential for homeostasis in diseases characterized by lung injury, inflammation, and remodeling as in chronic obstructive pulmonary disease (COPD). Chitinase 3-like-1 (CHI3L1) is induced in many cancers where it inhibits adaptive immune responses by stimulating immune checkpoint molecules (ICPs) and portends a poor prognosis. CHI3L1 is also induced in COPD where it regulates epithelial cell death. In this study, we demonstrate that pulmonary melanoma metastasis inhibits macrophage phagocytosis by stimulating the CD47-SIRPα and CD24-Siglec10 phagocytosis checkpoint pathways while inhibiting macrophage "eat me" signals from calreticulin and HMGB1. We also demonstrate that these effects on macrophage phagocytosis are associated with CHI3L1 stimulation of the SHP-1 and SHP-2 phosphatases and inhibition of the accumulation and phosphorylation of cytoskeleton-regulating nonmuscle myosin IIa. This inhibition of innate immune responses such as phagocytosis provides a mechanistic explanation for the ability of CHI3L1 to stimulate ICPs and inhibit adaptive immune responses in cancer and diseases such as COPD. The ability of CHI3L1 to simultaneously inhibit innate immune responses, stimulate ICPs, inhibit T cell costimulation, and regulate a number of other oncogenic and inflammation pathways suggests that CHI3L1-targeted therapeutics are promising interventions in cancer, COPD, and other disorders.

Inherited causes of clonal haematopoiesis in 97,691 whole genomes.

Age is the dominant risk factor for most chronic human diseases, but the mechanisms through which ageing confers this risk are largely unknown1. The age-related acquisition of somatic mutations that lead to clonal expansion in regenerating haematopoietic stem cell populations has recently been associated with both haematological cancer2-4 and coronary heart disease5-this phenomenon is termed clonal haematopoiesis of indeterminate potential (CHIP)6. Simultaneous analyses of germline and somatic whole-genome sequences provide the opportunity to identify root causes of CHIP. Here we analyse high-coverage whole-genome sequences from 97,691 participants of diverse ancestries in the National Heart, Lung, and Blood Institute Trans-omics for Precision Medicine (TOPMed) programme, and identify 4,229 individuals with CHIP. We identify associations with blood cell, lipid and inflammatory traits that are specific to different CHIP driver genes. Association of a genome-wide set of germline genetic variants enabled the identification of three genetic loci associated with CHIP status, including one locus at TET2 that was specific to individuals of African ancestry. In silico-informed in vitro evaluation of the TET2 germline locus enabled the identification of a causal variant that disrupts a TET2 distal enhancer, resulting in increased self-renewal of haematopoietic stem cells. Overall, we observe that germline genetic variation shapes haematopoietic stem cell function, leading to CHIP through mechanisms that are specific to clonal haematopoiesis as well as shared mechanisms that lead to somatic mutations across tissues.

Detecting and dissecting signaling crosstalk via the multilayer network integration of signaling and regulatory interactions.

The versatility of cellular response arises from the communication, or crosstalk, of signaling pathways in a complex network of signaling and transcriptional regulatory interactions. Understanding the various mechanisms underlying crosstalk on a global scale requires untargeted computational approaches. We present a network-based statistical approach, MuXTalk, that uses high-dimensional edges called multilinks to model the unique ways in which signaling and regulatory interactions can interface. We demonstrate that the signaling-regulatory interface is located primarily in the intermediary region between signaling pathways where crosstalk occurs, and that multilinks can differentiate between distinct signaling-transcriptional mechanisms. Using statistically over-represented multilinks as proxies of crosstalk, we infer crosstalk among 60 signaling pathways, expanding currently available crosstalk databases by more than five-fold. MuXTalk surpasses existing methods in terms of model performance metrics, identifies additions to manual curation efforts, and pinpoints potential mediators of crosstalk. Moreover, it accommodates the inherent context-dependence of crosstalk, allowing future applications to cell type- and disease-specific crosstalk.

Genetic regulators of sputum mucin concentration and their associations with COPD phenotypes.

Hyper-secretion and/or hyper-concentration of mucus is a defining feature of multiple obstructive lung diseases, including chronic obstructive pulmonary disease (COPD). Mucus itself is composed of a mixture of water, ions, salt and proteins, of which the gel-forming mucins, MUC5AC and MUC5B, are the most abundant. Recent studies have linked the concentrations of these proteins in sputum to COPD phenotypes, including chronic bronchitis (CB) and acute exacerbations (AE). We sought to determine whether common genetic variants influence sputum mucin concentrations and whether these variants are also associated with COPD phenotypes, specifically CB and AE. We performed a GWAS to identify quantitative trait loci for sputum mucin protein concentration (pQTL) in the Sub-Populations and InteRmediate Outcome Measures in COPD Study (SPIROMICS, n = 708 for total mucin, n = 215 for MUC5AC, MUC5B). Subsequently, we tested for associations of mucin pQTL with CB and AE using regression modeling (n = 822-1300). Replication analysis was conducted using data from COPDGene (n = 5740) and by examining results from the UK Biobank. We identified one genome-wide significant pQTL for MUC5AC (rs75401036) and two for MUC5B (rs140324259, rs10001928). The strongest association for MUC5B, with rs140324259 on chromosome 11, explained 14% of variation in sputum MUC5B. Despite being associated with lower MUC5B, the C allele of rs140324259 conferred increased risk of CB (odds ratio (OR) = 1.42; 95% confidence interval (CI): 1.10-1.80) as well as AE ascertained over three years of follow up (OR = 1.41; 95% CI: 1.02-1.94). Associations between rs140324259 and CB or AE did not replicate in COPDGene. However, in the UK Biobank, rs140324259 was associated with phenotypes that define CB, namely chronic mucus production and cough, again with the C allele conferring increased risk. We conclude that sputum MUC5AC and MUC5B concentrations are associated with common genetic variants, and the top locus for MUC5B may influence COPD phenotypes, in particular CB.

FGF20 and PGM2 variants are associated with childhood asthma in family-based whole-genome sequencing studies.

BACKGROUND: Asthma is a heterogeneous common respiratory disease that remains poorly understood. The established genetic associations fail to explain the high estimated heritability, and the prevalence of asthma differs between populations and geographic regions. Robust association analyses incorporating different genetic ancestries and whole-genome sequencing data may identify novel genetic associations. METHODS: We performed family-based genome-wide association analyses of childhood-onset asthma based on whole-genome sequencing (WGS) data for the 'The Genetic Epidemiology of Asthma in Costa Rica' study (GACRS) and the Childhood Asthma Management Program (CAMP). Based on parent-child trios with children diagnosed with asthma, we performed a single variant analysis using an additive and a recessive genetic model and a region-based association analysis of low-frequency and rare variants. RESULTS: Based on 1180 asthmatic trios (894 GACRS trios and 286 CAMP trios, a total of 3540 samples with WGS data), we identified three novel genetic loci associated with childhood-onset asthma: rs4832738 on 4p14 ($P=1.72\ast{10}^{-9}$, recessive model), rs1581479 on 8p22 ($P=1.47\ast{10}^{-8}$, additive model) and rs73367537 on 10q26 ($P=1.21\ast{10}^{-8}$, additive model in GACRS only). Integrative analyses suggested potential novel candidate genes underlying these associations: PGM2 on 4p14 and FGF20 on 8p22. CONCLUSION: Our family-based whole-genome sequencing analysis identified three novel genetic loci for childhood-onset asthma. Gene expression data and integrative analyses point to PGM2 on 4p14 and FGF20 on 8p22 as linked genes. Furthermore, region-based analyses suggest independent potential low-frequency/rare variant associations on 8p22. Follow-up analyses are needed to understand the functional mechanisms and generalizability of these associations.

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.

A statistical physics approach for disease module detection.

Extensive evidence indicates that the pathobiological processes of a complex disease are associated with perturbation in specific neighborhoods of the human protein-protein interaction (PPI) network (also known as the interactome), often referred to as the disease module. Many computational methods have been developed to integrate the interactome and omics profiles to extract context-dependent disease modules. Yet, existing methods all have fundamental limitations in terms of rigor and/or efficiency. Here, we developed a statistical physics approach based on the random-field Ising model (RFIM) for disease module detection, which is both mathematically rigorous and computationally efficient. We applied our RFIM approach to genome-wide association studies (GWAS) of ten complex diseases to examine its performance for disease module detection. We found that our RFIM approach outperforms existing methods in terms of computational efficiency, connectivity of disease modules, and robustness to the interactome incompleteness.

Variants in autophagy genes MTMR12 and FAM134A are putative modifiers of the hepatic phenotype in α1-antitrypsin deficiency.

BACKGROUND AND AIMS: In the classical form of α1-antitrypsin deficiency, a misfolded variant α1-antitrypsin Z accumulates in the endoplasmic reticulum of liver cells and causes liver cell injury by gain-of-function proteotoxicity in a sub-group of affected homozygotes but relatively little is known about putative modifiers. Here, we carried out genomic sequencing in a uniquely affected family with an index case of liver failure and 2 homozygous siblings with minimal or no liver disease. Their sequences were compared to sequences in well-characterized cohorts of homozygotes with or without liver disease, and then candidate sequence variants were tested for changes in the kinetics of α1-antitrypsin variant Z degradation in iPS-derived hepatocyte-like cells derived from the affected siblings themselves. APPROACH AND RESULTS: Specific variants in autophagy genes MTMR12 and FAM134A could each accelerate the degradation of α1-antitrypsin variant Z in cells from the index patient, but both MTMR12 and FAM134A variants were needed to slow the degradation of α1-antitrypsin variant Z in cells from a protected sib, indicating that inheritance of both variants is needed to mediate the pathogenic effects of hepatic proteotoxicity at the cellular level. Analysis of homozygote cohorts showed that multiple patient-specific variants in proteostasis genes are likely to explain liver disease susceptibility at the population level. CONCLUSIONS: These results validate the concept that genetic variation in autophagy function can determine susceptibility to liver disease in α1-antitrypsin deficiency and provide evidence that polygenic mechanisms and multiple patient-specific variants are likely needed for proteotoxic pathology.

Partial correlation network analysis identifies coordinated gene expression within a regional cluster of COPD genome-wide association signals.

Chronic obstructive pulmonary disease (COPD) is a complex disease influenced by well-established environmental exposures (most notably, cigarette smoking) and incompletely defined genetic factors. The chromosome 4q region harbors multiple genetic risk loci for COPD, including signals near HHIP, FAM13A, GSTCD, TET2, and BTC. Leveraging RNA-Seq data from lung tissue in COPD cases and controls, we estimated the co-expression network for genes in the 4q region bounded by HHIP and BTC (~70MB), through partial correlations informed by protein-protein interactions. We identified several co-expressed gene pairs based on partial correlations, including NPNT-HHIP, BTC-NPNT and FAM13A-TET2, which were replicated in independent lung tissue cohorts. Upon clustering the co-expression network, we observed that four genes previously associated to COPD: BTC, HHIP, NPNT and PPM1K appeared in the same network community. Finally, we discovered a sub-network of genes differentially co-expressed between COPD vs controls (including FAM13A, PPA2, PPM1K and TET2). Many of these genes were previously implicated in cell-based knock-out experiments, including the knocking out of SPP1 which belongs to the same genomic region and could be a potential local key regulatory gene. These analyses identify chromosome 4q as a region enriched for COPD genetic susceptibility and differential co-expression.

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.

Lung Transcriptomics Links Emphysema to Barrier Dysfunction and Macrophage Subpopulations.

RATIONALE: While many studies have examined gene expression in lung tissue, the gene regulatory processes underlying emphysema are still not well understood. Finding efficient non-imaging screening methods and disease-modifying therapies has been challenging, but knowledge of the transcriptomic features of emphysema may help in this effort. OBJECTIVES: Our goals were to identify emphysema-associated biological pathways through transcriptomic analysis of bulk lung tissue, to determine the lung cell types in which these emphysema-associated pathways are altered, and to detect unique and overlapping transcriptomic signatures in blood and lung samples. METHODS: Using RNA-sequencing data from 446 samples in the Lung Tissue Research Consortium (LTRC) and 3,606 blood samples from the COPDGene study, we examined the transcriptomic features of chest computed tomography-quantified emphysema. We also leveraged publicly available lung single-cell RNA-sequencing data to identify cell types showing COPD-associated differential expression of the emphysema pathways found in the bulk analyses. MEASUREMENTS AND MAIN RESULTS: In the bulk lung RNA-seq analysis, 1,087 differentially expressed genes and 34 dysregulated pathways were significantly associated with emphysema. We observed alternative splicing of several genes and increased activity in pluripotency and cell barrier function pathways. Lung tissue and blood samples shared differentially expressed genes and biological pathways. Multiple lung cell types displayed dysregulation of epithelial barrier function pathways, and distinct pathway activities were observed among various macrophage subpopulations. CONCLUSIONS: This study identified emphysema-related changes in gene expression and alternative splicing, cell-type specific dysregulated pathways, and instances of shared pathway dysregulation between blood and lung.

Genetically Predicted Body Mass Index and Mortality in Chronic Obstructive Pulmonary Disease.

Rationale: Body mass index (BMI) is associated with chronic obstructive pulmonary disease (COPD) mortality, but the underlying mechanisms are unclear. The effect of genetic variants aggregated into a polygenic score may elucidate the 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 (PGS) for BMI (PGSBMI) and tested for associations of the PGSBMI with all-cause, respiratory, and cardiovascular mortality in participants with COPD from the COPDGene (Genetic Epidemiology of COPD), ECLIPSE (Evaluation of COPD Longitudinally to Identify Predictive Surrogate End-points), 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 the 20th and 80th percentiles), and discordantly high (BMIdiff >80th percentile) BMI. We applied Cox models, examined potential nonlinear associations of the PGSBMI and BMIdiff with mortality, and summarized results with meta-analysis. Measurements and Main Results: We observed significant nonlinear 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 (hazard ratio [HR], 1.29; 95% confidence interval [CI], 1.12-1.49), but not for respiratory or all-cause mortality. Compared with participants with concordant measured and genetically predicted BMI, those with discordantly low BMI had higher risks for all-cause mortality (HR, 1.57; 95% CI, 1.41-1.74) and respiratory death (HR, 2.01; 95% CI, 1.61-2.51). Conclusions: In people with COPD, a higher genetically predicted BMI is associated with higher cardiovascular mortality but not respiratory mortality. Individuals with a discordantly low BMI have higher all-cause and respiratory mortality rates than those with a concordant BMI.

Development, Progression, and Mortality of Suspected Interstitial Lung Disease in COPDGene.

RATIONALE: Some with interstitial lung abnormalities (ILA) have suspected interstitial lung disease (ILD), a subgroup with adverse outcomes. Rates of development and progression of suspected ILD and their effect on mortality are unknown. OBJECTIVES: To determine rates of development and progression of suspected ILD and assess effects of individual ILD and progression criteria on mortality. METHODS: Participants from COPDGene were included. ILD was defined as ILA and fibrosis and/or FVC <80% predicted. Prevalent ILD was assessed at enrollment, incident ILD and progression at 5-year follow-up. CT progression was assessed visually and FVC decline as relative change. Multivariable Cox regression tested associations between mortality and ILD groups. RESULTS: Of 9,588 participants at enrollment, 267 (2.8%) had prevalent ILD. Those with prevalent ILD had 52% mortality after median 10.6 years, which was higher than ILA (33%; HR=2.0; p<0.001). The subgroup of prevalent ILD with fibrosis only had worse mortality (59%) than ILA (HR=2.2; p<0.001). 97 participants with prevalent ILD completed 5-year follow-up: 32% had stable CT and relative FVC decline <10%, 6% FVC decline ≥10% only, 39% CT progression only, and 22% both CT progression and FVC decline ≥10%. Mortality rates were 32%, 50%, 45%, and 46% respectively; those with CT progression only had worse mortality than ILA (HR=2.6; p=0.005). At 5-year follow-up, incident ILD occurred in 168/4,843 participants without prevalent ILD and had worse mortality than ILA (HR=2.5; p<0.001). CONCLUSION: Rates of mortality and progression are high among those with suspected ILD in COPDGene; fibrosis and radiologic progression are important predictors of mortality.

Airway Mucus Plugs on Chest Computed Tomography Are Associated with Exacerbations in COPD.

Rationale/Objective: Acute exacerbations (AEs) of chronic obstructive pulmonary disease (COPD) are associated with significant morbidity and mortality. Whether mucus plugs are associated with prospective exacerbations has not been examined extensively. METHODS: Mucus plugs were visually-identified on baseline chest computed tomography (CT) scans from smokers with Global Initiative for Chronic Obstructive Lung Disease (GOLD) grades 2-4 COPD enrolled in two multicenter cohort studies: Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints (ECLIPSE) and COPDGene. Associations between ordinal mucus plug score categories (0/1-2/≥3) and prospectively-ascertained AEs, defined as worsening respiratory symptoms requiring systemic steroids and/or antibiotics (moderate-to-severe) and/or ER/hospitalization (severe), were assessed using multivariable-adjusted zero-inflated Poisson regression; subjects were exacerbation-free at enrollment. RESULTS: Among 3,250 participants in COPDGene (mean±SD age 63.7±8.4 years, FEV1 50.6%±17.8% predicted, 45.1% female) and 1,716 participants in ECLIPSE (age 63.3±7.1 years, FEV1 48.3%±15.8% predicted, 36.2% female), 44.4% and 46.0% had mucus plugs, respectively. The incidence rates of AEs were 61.0 (COPDGene) and 125.7 (ECLIPSE) per 100 person-years. Relative to those without mucus plugs, the presence of 1-2 and ≥3 mucus plugs was associated with increased risk (adjusted rate ratio, aRR [95%CI]=1.07[1.05-1.09] and 1.15[1.1-1.2] in COPDGene; aRR=1.06[1.02-1.09] and 1.12[1.04-1.2] in ECLIPSE, respectively) for prospective moderate-to-severe AEs. The presence of 1-2 and ≥3 mucus plugs was also associated with increased risk for severe AEs during follow-up (aRR=1.05[1.01-1.08] and 1.09[1.02-1.18] in COPDGene; aRR=1.17[1.07-1.27] and 1.37[1.15-1.62] in ECLIPSE, respectively). CONCLUSION: CT-based mucus plugs are associated with an increased risk for future COPD AEs.

Early Evidence of Chronic Obstructive Pulmonary Disease Obscured by Race-Specific Prediction Equations.

Rationale: The identification of early chronic obstructive pulmonary disease (COPD) is essential to appropriately counsel patients regarding smoking cessation, provide symptomatic treatment, and eventually develop disease-modifying treatments. Disease severity in COPD is defined using race-specific spirometry equations. These may disadvantage non-White individuals in diagnosis and care. Objectives: Determine the impact of race-specific equations on African American (AA) versus non-Hispanic White individuals. Methods: Cross-sectional analyses of the COPDGene (Genetic Epidemiology of Chronic Obstructive Pulmonary Disease) cohort were conducted, comparing non-Hispanic White (n = 6,766) and AA (n = 3,366) participants for COPD manifestations. Measurements and Main Results: Spirometric classifications using race-specific, multiethnic, and "race-reversed" prediction equations (NHANES [National Health and Nutrition Examination Survey] and Global Lung Function Initiative "Other" and "Global") were compared, as were respiratory symptoms, 6-minute-walk distance, computed tomography imaging, respiratory exacerbations, and St. George's Respiratory Questionnaire. Application of different prediction equations to the cohort resulted in different classifications by stage, with NHANES and Global Lung Function Initiative race-specific equations being minimally different, but race-reversed equations moving AA participants to more severe stages and especially between the Global Initiative for Chronic Obstructive Lung Disease (GOLD) stage 0 and preserved ratio impaired spirometry groups. Classification using the established NHANES race-specific equations demonstrated that for each of GOLD stages 1-4, AA participants were younger, had fewer pack-years and more current smoking, but had more exacerbations, shorter 6-minute-walk distance, greater dyspnea, and worse BODE (body mass index, airway obstruction, dyspnea, and exercise capacity) scores and St. George's Respiratory Questionnaire scores. Differences were greatest in GOLD stages 1 and 2. Race-reversed equations reclassified 774 AA participants (43%) from GOLD stage 0 to preserved ratio impaired spirometry. Conclusions: Race-specific equations underestimated disease severity among AA participants. These effects were particularly evident in early disease and may result in late detection of COPD.

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.