COPD Subtypes Are Differentially Associated with Cardiovascular Events and COPD Exacerbations.

BACKGROUND: The coronary artery calcium score (CACS) and ratio of the pulmonary artery to aorta diameters (PA:A ratio) measured from chest CT scans have been established as predictors of cardiovascular events and chronic obstructive pulmonary disease (COPD) exacerbations, respectively. However, little is known about the reciprocal relationship between these predictors and outcomes. Furthermore, the prognostic implications of COPD subtypes on clinical outcomes remain insufficiently characterized. RESEARCH QUESTION: How can these two chest CT-derived parameters predict subsequent cardiovascular events and COPD exacerbations in different COPD subtypes? STUDY DESIGN AND METHODS: Using COPDGene study data, we assessed prospective cardiovascular disease (CVD) and COPD exacerbation risk in COPD subjects (Global Initiative for Chronic Obstructive Lung Disease spirometric grades 2-4), focusing on CACS and PA:A ratio at study enrollment, with logistic regression models. These outcomes were analyzed in three COPD subtypes: 1,042 Non-emphysema-predominant COPD (NEPD; low attenuation area at -950 Hounsfield units [LAA-950]<5%), 1,324 Emphysema-predominant COPD (EPD; LAA-950≥10%), and 465 Intermediate Emphysema COPD (IE; 5≤LAA-950<10%). RESULTS: Our study indicated significantly higher overall risk for cardiovascular events in subjects with higher CACS (≥median; Odds Ratio (OR): 1.61, 95% Confidence Interval (CI)=1.30-2.00) and increased COPD exacerbations in those with higher PA:A ratios (≥1; OR: 1.80, 95% CI=1.46-2.23). Notably, NEPD subjects showed a stronger association between these indicators and clinical events compared to EPD (with CACS/CVD, NEPD vs. EPD, OR 2.02 vs. 1.41; with PA:A ratio/COPD exacerbation, NEPD vs. EPD, OR 2.50 vs. 1.65); the difference in odds ratios between COPD subtypes was statistically significant for CACS/CVD. INTERPRETATION: Two chest CT parameters, CACS and PA:A ratio, hold distinct predictive values for cardiovascular events and COPD exacerbations that are influenced by specific COPD subtypes. TRIAL REGISTRATION: ClinicalTrials.gov Identifier: NCT00608764.

Sex-biased Regulation of Extracellular Matrix Genes in COPD.

Compared to men, women often develop COPD at an earlier age with worse respiratory symptoms despite lower smoking exposure. However, most preventive, and therapeutic strategies ignore biological sex differences in COPD. Our goal was to better understand sex-specific gene regulatory processes in lung tissue and the molecular basis for sex differences in COPD onset and severity. We analyzed lung tissue gene expression and DNA methylation data from 747 individuals in the Lung Tissue Research Consortium (LTRC), and 85 individuals in an independent dataset. We identified sex differences in COPD-associated gene regulation using gene regulatory networks. We used linear regression to test for sex-biased associations of methylation with lung function, emphysema, smoking, and age. Analyzing gene regulatory networks in the control group, we identified that genes involved in the extracellular matrix (ECM) have higher transcriptional factor targeting in females than in males. However, this pattern is reversed in COPD, with males showing stronger regulatory targeting of ECM-related genes than females. Smoking exposure, age, lung function, and emphysema were all associated with sex-specific differential methylation of ECM-related genes. We identified sex-based gene regulatory patterns of ECM-related genes associated with lung function and emphysema. Multiple factors including epigenetics, smoking, aging, and cell heterogeneity influence sex-specific gene regulation in COPD. Our findings underscore the importance of considering sex as a key factor in disease susceptibility and severity.

Identification of factors directly linked to incident chronic obstructive pulmonary disease: A causal graph modeling study.

BACKGROUND: Beyond exposure to cigarette smoking and aging, the factors that influence lung function decline to incident chronic obstructive pulmonary disease (COPD) remain unclear. Advancements have been made in categorizing COPD into emphysema and airway predominant disease subtypes; however, predicting which healthy individuals will progress to COPD is difficult because they can exhibit profoundly different disease trajectories despite similar initial risk factors. This study aimed to identify clinical, genetic, and radiological features that are directly linked-and subsequently predict-abnormal lung function. METHODS AND FINDINGS: We employed graph modeling on 2,643 COPDGene participants (aged 45 to 80 years, 51.25% female, 35.1% African Americans; enrollment 11/2007-4/2011) with smoking history but normal spirometry at study enrollment to identify variables that are directly linked to future lung function abnormalities. We developed logistic regression and random forest predictive models for distinguishing individuals who maintain lung function from those who decline. Of the 131 variables analyzed, 6 were identified as informative to future lung function abnormalities, namely forced expiratory flow in the middle range (FEF25-75%), average lung wall thickness in a 10 mm radius (Pi10), severe emphysema, age, sex, and height. We investigated whether these features predict individuals leaving GOLD 0 status (normal spirometry according to Global Initiative for Obstructive Lung Disease (GOLD) criteria). Linear models, trained with these features, were quite predictive (area under receiver operator characteristic curve or AUROC = 0.75). Random forest predictors performed similarly to logistic regression (AUROC = 0.7), indicating that no significant nonlinear effects were present. The results were externally validated on 150 participants from Specialized Center for Clinically Oriented Research (SCCOR) cohort (aged 45 to 80 years, 52.7% female, 4.7% African Americans; enrollment: 7/2007-12/2012) (AUROC = 0.89). The main limitation of longitudinal studies with 5- and 10-year follow-up is the introduction of mortality bias that disproportionately affects the more severe cases. However, our study focused on spirometrically normal individuals, who have a lower mortality rate. Another limitation is the use of strict criteria to define spirometrically normal individuals, which was unavoidable when studying factors associated with changes in normalized forced expiratory volume in 1 s (FEV1%predicted) or the ratio of FEV1/FVC (forced vital capacity). CONCLUSIONS: This study took an agnostic approach to identify which baseline measurements differentiate and predict the early stages of lung function decline in individuals with previous smoking history. Our analysis suggests that emphysema affects obstruction onset, while airway predominant pathology may play a more important role in future FEV1 (%predicted) decline without obstruction, and FEF25-75% may affect both.

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.

Association of PHACTR1 with Coronary Artery Calcium Differs by Sex and Cigarette Smoking.

Background: Coronary artery calcium (CAC) is a marker of subclinical atherosclerosis and is a complex heritable trait with both genetic and environmental risk factors, including sex and smoking. Methods: We performed genome-wide association (GWA) analyses for CAC among all participants and stratified by sex in the COPDGene study (n = 6144 participants of European ancestry and n = 2589 participants of African ancestry) with replication in the Diabetes Heart Study (DHS). We adjusted for age, sex, current smoking status, BMI, diabetes, self-reported high blood pressure, self-reported high cholesterol, and genetic ancestry (as summarized by principal components computed within each racial group). For the significant signals from the GWA analyses, we examined the single nucleotide polymorphism (SNP) by sex interactions, stratified by smoking status (current vs. former), and tested for a SNP by smoking status interaction on CAC. Results: We identified genome-wide significant associations for CAC in the chromosome 9p21 region [CDKN2B-AS1] among all COPDGene participants (p = 7.1 × 10-14) and among males (p = 1.0 × 10-9), but the signal was not genome-wide significant among females (p = 6.4 × 10-6). For the sex stratified GWA analyses among females, the chromosome 6p24 region [PHACTR1] had a genome-wide significant association (p = 4.4 × 10-8) with CAC, but this signal was not genome-wide significant among all COPDGene participants (p = 1.7 × 10-7) or males (p = 0.03). There was a significant interaction for the SNP rs9349379 in PHACTR1 with sex (p = 0.02), but the interaction was not significant for the SNP rs10757272 in CDKN2B-AS1 with sex (p = 0.21). In addition, PHACTR1 had a stronger association with CAC among current smokers (p = 6.2 × 10-7) than former smokers (p = 7.5 × 10-3) and the SNP by smoking status interaction was marginally significant (p = 0.03). CDKN2B-AS1 had a strong association with CAC among both former (p = 7.7 × 10-8) and current smokers (p = 1.7 × 10-7) and the SNP by smoking status interaction was not significant (p = 0.40). Conclusions: Among current and former smokers of European ancestry in the COPDGene study, we identified a genome-wide significant association in the chromosome 6p24 region [PHACTR1] with CAC among females, but not among males. This region had a significant SNP by sex and SNP by smoking interaction on CAC.

Overlap between COPD genetic association results and transcriptional quantitative trait loci.

Rationale: Genome-wide association studies (GWAS) have identified multiple genetic loci associated with chronic obstructive pulmonary disease (COPD). When integrated with GWAS results, expression quantitative trait locus (eQTL) studies can provide insight into biological mechanisms involved in disease by identifying single nucleotide polymorphisms (SNPs) that contribute to whole gene expression. However, there are multiple genetically driven regulatory and isoform-specific effects which cannot be detected in traditional eQTL analyses. Here, we identify SNPs that are associated with alternative splicing (sQTL) in addition to eQTLs to identify novel functions for COPD associated genetic variants. Methods: We performed RNA sequencing on whole blood from 3743 subjects in the COPDGene Study. RNA sequencing data from lung tissue of 1241 subjects from the Lung Tissue Research Consortium (LTRC), and whole genome sequencing data on all subjects. Associations between all SNPs within 1000 kb of a gene (cis-) and splice and gene expression quantifications were tested using tensorQTL. In COPDGene a total of 11,869,333 SNPs were tested for association with 58,318 splice clusters, and 8,792,206 SNPs were tested for association with 70,094 splice clusters in LTRC. We assessed colocalization with COPD-associated SNPs from a published GWAS[1]. Results: After adjustment for multiple statistical testing, we identified 28,110 splice-sites corresponding to 3,889 unique genes that were significantly associated with genotype in COPDGene whole blood, and 58,258 splice-sites corresponding to 10,307 unique genes associated with genotype in LTRC lung tissue. We found 7,576 sQTL splice-sites corresponding to 2,110 sQTL genes were shared between whole blood and lung, while 20,534 sQTL splice-sites in 3,518 genes were unique to blood and 50,682 splice-sites in 9,677 genes were unique to lung. To determine what proportion of COPD-associated SNPs were associated with transcriptional splicing, we performed colocalization analysis between COPD GWAS and sQTL data, and found that 38 genomic windows, corresponding to 38 COPD GWAS loci had evidence of colocalization between QTLs and COPD. The top five colocalizations between COPD and lung sQTLs include NPNT , FBXO38 , HHIP , NTN4 and BTC . Conclusions: A total of 38 COPD GWAS loci contain evidence of sQTLs, suggesting that analysis of sQTLs in whole blood and lung tissue can provide novel insights into disease mechanisms.

Epigenetic and proteomic signatures associate with clonal hematopoiesis expansion rate.

Clonal hematopoiesis of indeterminate potential (CHIP), whereby somatic mutations in hematopoietic stem cells confer a selective advantage and drive clonal expansion, not only correlates with age but also confers increased risk of morbidity and mortality. Here, we leverage genetically predicted traits to identify factors that determine CHIP clonal expansion rate. We used the passenger-approximated clonal expansion rate method to quantify the clonal expansion rate for 4,370 individuals in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) cohort and calculated polygenic risk scores for DNA methylation aging, inflammation-related measures and circulating protein levels. Clonal expansion rate was significantly associated with both genetically predicted and measured epigenetic clocks. No associations were identified with inflammation-related lab values or diseases and CHIP expansion rate overall. A proteome-wide search identified predicted circulating levels of myeloid zinc finger 1 and anti-Müllerian hormone as associated with an increased CHIP clonal expansion rate and tissue inhibitor of metalloproteinase 1 and glycine N-methyltransferase as associated with decreased CHIP clonal expansion rate. Together, our findings identify epigenetic and proteomic patterns associated with the rate of hematopoietic clonal expansion.

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.

Multi-center Korean cohort study based on RNA-sequencing data targeting COPD patients.

In 2023, WHO ranked chronic obstructive pulmonary disease (COPD) as the third leading cause of death, with 3.23 million fatalities in 2019. The intricate nature of the disease, which is influenced by genetics, environment, and lifestyle, is evident. The effect of air pollution and changes in atmospheric substances because of global warming highlight the need for this research. These environmental shifts are associated with the emergence of various respiratory infections such as COVID-19. RNA sequencing is pivotal in airway diseases, including COPD, as it enables comprehensive transcriptome analysis, biomarker discovery, and uncovers novel pathways. It facilitates personalized medicine by tracking dynamic changes in gene expression in response to various triggers. However, the limited research on East Asian populations may overlook the unique nuances of COPD development and progression. Bridging this gap and using peripheral blood samples for systemic analysis are crucial for comprehensive and globally applicable COPD diagnosis and treatment.

Rheumatoid arthritis and changes on spirometry by smoking status in two prospective longitudinal cohorts.

OBJECTIVE: To compare longitudinal changes in spirometric measures between patients with rheumatoid arthritis (RA) and non-RA comparators. METHODS: We analysed longitudinal data from two prospective cohorts: the UK Biobank and COPDGene. Spirometry was conducted at baseline and a second visit after 5-7 years. RA was identified based on self-report and disease-modifying antirheumatic drug use; non-RA comparators reported neither. The primary outcomes were annual changes in the per cent-predicted forced expiratory volume in 1 s (FEV1%) and per cent predicted forced vital capacity (FVC%). Statistical comparisons were performed using multivariable linear regression. The analysis was stratified based on baseline smoking status and the presence of obstructive pattern (FEV1/FVC <0.7). RESULTS: Among participants who underwent baseline and follow-up spirometry, we identified 233 patients with RA and 37 735 non-RA comparators. Among never-smoking participants without an obstructive pattern, RA was significantly associated with more FEV1% decline (ß=-0.49, p=0.04). However, in ever smokers with ≥10 pack-years, those with RA exhibited significantly less FEV1% decline than non-RA comparators (ß=0.50, p=0.02). This difference was more pronounced among those with an obstructive pattern at baseline (ß=1.12, p=0.01). Results were similar for FEV1/FVC decline. No difference was observed in the annual FVC% change in RA versus non-RA. CONCLUSIONS: Smokers with RA, especially those with baseline obstructive spirometric patterns, experienced lower FEV1% and FEV1/FVC decline than non-RA comparators. Conversely, never smokers with RA had more FEV1% decline than non-RA comparators. Future studies should investigate potential treatments and the pathogenesis of obstructive lung diseases in smokers with RA.

Polygenic and transcriptional risk scores identify chronic obstructive pulmonary disease subtypes.

Rationale: Genetic variants and gene expression predict risk of chronic obstructive pulmonary disease (COPD), but their effect on COPD heterogeneity is unclear. Objectives: Define high-risk COPD subtypes using both genetics (polygenic risk score, PRS) and blood gene expression (transcriptional risk score, TRS) and assess differences in clinical and molecular characteristics. Methods: We defined high-risk groups based on PRS and TRS quantiles by maximizing differences in protein biomarkers in a COPDGene training set and identified these groups in COPDGene and ECLIPSE test sets. We tested multivariable associations of subgroups with clinical outcomes and compared protein-protein interaction networks and drug repurposing analyses between high-risk groups. Measurements and Main Results: We examined two high-risk omics-defined groups in non-overlapping test sets (n=1,133 NHW COPDGene, n=299 African American (AA) COPDGene, n=468 ECLIPSE). We defined "High activity" (low PRS/high TRS) and "severe risk" (high PRS/high TRS) subgroups. Participants in both subgroups had lower body-mass index (BMI), lower lung function, and alterations in metabolic, growth, and immune signaling processes compared to a low-risk (low PRS, low TRS) reference subgroup. "High activity" but not "severe risk" participants had greater prospective FEV 1 decline (COPDGene: -51 mL/year; ECLIPSE: - 40 mL/year) and their proteomic profiles were enriched in gene sets perturbed by treatment with 5-lipoxygenase inhibitors and angiotensin-converting enzyme (ACE) inhibitors. Conclusions: Concomitant use of polygenic and transcriptional risk scores identified clinical and molecular heterogeneity amongst high-risk individuals. Proteomic and drug repurposing analysis identified subtype-specific enrichment for therapies and suggest prior drug repurposing failures may be explained by patient selection.

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.

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

Chronic Obstructive Pulmonary Disease Exacerbations Increase the Risk of Subsequent Cardiovascular Events: A Longitudinal Analysis of the COPDGene Study.

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.

Determinants of mosaic chromosomal alteration fitness.

Clonal hematopoiesis (CH) is characterized by the acquisition of a somatic mutation in a hematopoietic stem cell that results in a clonal expansion. These driver mutations can be single nucleotide variants in cancer driver genes or larger structural rearrangements called mosaic chromosomal alterations (mCAs). The factors that influence the variations in mCA fitness and ultimately result in different clonal expansion rates are not well understood. We used the Passenger-Approximated Clonal Expansion Rate (PACER) method to estimate clonal expansion rate as PACER scores for 6,381 individuals in the NHLBI TOPMed cohort with gain, loss, and copy-neutral loss of heterozygosity mCAs. Our mCA fitness estimates, derived by aggregating per-individual PACER scores, were correlated (R2 = 0.49) with an alternative approach that estimated fitness of mCAs in the UK Biobank using population-level distributions of clonal fraction. Among individuals with JAK2 V617F clonal hematopoiesis of indeterminate potential or mCAs affecting the JAK2 gene on chromosome 9, PACER score was strongly correlated with erythrocyte count. In a cross-sectional analysis, genome-wide association study of estimates of mCA expansion rate identified a TCL1A locus variant associated with mCA clonal expansion rate, with suggestive variants in NRIP1 and TERT.

Correlation-based network integration of lung RNA sequencing and DNA methylation data in chronic obstructive pulmonary disease.

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.

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.

HHIP protein interactions in lung cells provide insight into COPD pathogenesis.

Chronic obstructive pulmonary disease (COPD) is the third leading cause of death worldwide. The primary causes of COPD are environmental, including cigarette smoking; however, genetic susceptibility also contributes to COPD risk. Genome-Wide Association Studies (GWASes) have revealed more than 80 genetic loci associated with COPD, leading to the identification of multiple COPD GWAS genes. However, the biological relationships between the identified COPD susceptibility genes are largely unknown. Genes associated with a complex disease are often in close network proximity, i.e. their protein products often interact directly with each other and/or similar proteins. In this study, we use affinity purification mass spectrometry (AP-MS) to identify protein interactions with HHIP , a well-established COPD GWAS gene which is part of the sonic hedgehog pathway, in two disease-relevant lung cell lines (IMR90 and 16HBE). To better understand the network neighborhood of HHIP , its proximity to the protein products of other COPD GWAS genes, and its functional role in COPD pathogenesis, we create HUBRIS, a protein-protein interaction network compiled from 8 publicly available databases. We identified both common and cell type-specific protein-protein interactors of HHIP. We find that our newly identified interactions shorten the network distance between HHIP and the protein products of several COPD GWAS genes, including DSP, MFAP2, TET2 , and FBLN5 . These new shorter paths include proteins that are encoded by genes involved in extracellular matrix and tissue organization. We found and validated interactions to proteins that provide new insights into COPD pathobiology, including CAVIN1 (IMR90) and TP53 (16HBE). The newly discovered HHIP interactions with CAVIN1 and TP53 implicate HHIP in response to oxidative stress.

Computational Deconvolution of Cell Type-Specific Gene Expression in COPD and IPF Lungs Reveals Disease Severity Associations.

RATIONALE: Chronic obstructive pulmonary disease (COPD) and idiopathic pulmonary fibrosis (IPF) are debilitating diseases associated with divergent histopathological changes in the lungs. At present, due to cost and technical limitations, profiling cell types is not practical in large epidemiology cohorts (n>1000). Here, we used computational deconvolution to identify cell types in COPD and IPF lungs whose abundances and cell type-specific gene expression are associated with disease diagnosis and severity. METHODS: We analyzed lung tissue RNA-seq data from 1026 subjects (COPD, n=465; IPF, n=213; control, n=348) from the Lung Tissue Research Consortium. We performed RNA-seq deconvolution, querying thirty-eight discrete cell-type varieties in the lungs. We tested whether deconvoluted cell-type abundance and cell type-specific gene expression were associated with disease severity. RESULTS: The abundance score of twenty cell types significantly differed between IPF and control lungs. In IPF subjects, eleven and nine cell types were significantly associated with forced vital capacity (FVC) and diffusing capacity for carbon monoxide (DLCO), respectively. Aberrant basaloid cells, a rare cells found in fibrotic lungs, were associated with worse FVC and DLCO in IPF subjects, indicating that this aberrant epithelial population increased with disease severity. Alveolar type 1 and vascular endothelial (VE) capillary A were decreased in COPD lungs compared to controls. An increase in macrophages and classical monocytes was associated with lower DLCO in IPF and COPD subjects. In both diseases, lower non-classical monocytes and VE capillary A cells were associated with increased disease severity. Alveolar type 2 cells and alveolar macrophages had the highest number of genes with cell type-specific differential expression by disease severity in COPD and IPF. In IPF, genes implicated in the pathogenesis of IPF, such as matrix metallopeptidase 7, growth differentiation factor 15, and eph receptor B2, were associated with disease severity in a cell type-specific manner. CONCLUSION: Utilization of RNA-seq deconvolution enabled us to pinpoint cell types present in the lungs that are associated with the severity of COPD and IPF. This knowledge offers valuable insight into the alterations within tissues in more advanced illness, ultimately providing a better understanding of the underlying pathological processes that drive disease progression.

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.