Barnard, J

Abstract Introduction Interstitial lung disease-associated pulmonary hypertension (ILD-PH) is a clinically challenging entity that lies on a pathophysiologic spectrum between pulmonary arterial hypertension (PAH) and interstitial lung disease (ILD). We compared the metabolomes of patients with ILD-PH to those with PAH and ILD to characterize the PAH-ILD axis and determine whether ILD-PH patients can be phenotyped based on their metabolomic similarity to PAH. Methods Using untargeted metabolomics profiles from patients with PAH, ILD, and ILD-PH in the large, U.S.-based, NIH-funded, multi-center PVDOMICS cohort, we trained an anchor-based classifier on PAH and ILD metabolomes using supervised PLS-DA with a 75/25 train/test split. Primary and secondary PH or risk-group classifications in PVDOMICS were based on expert adjudication of clinical data including right heart catheterization hemodynamics, pulmonary function testing, chest imaging, and laboratory results. We then applied the trained model to each ILD-PH patient to assign a PAH similarity score. We classified patients as having a PAH-like or ILD-like metabolomic phenotype using an optimal PAH similarity score cutoff of 0.5 then compared PAH metabolomic similarity to expert-adjudicated primary and secondary WSPH classifications in PVDOMICS. We implemented FELLA, a network-based metabolomics enrichment method, to holistically characterize the biologic underpinnings of PAH, ILD-PH, and ILD metabolomes. Results We analyzed 244 patients with PAH, 54 with ILD, and 47 with ILD-PH. The anchor-based classifier demonstrated high train-test performance with a c-statistic of 0.81 and 78% agreement with expert-adjudicated diagnoses. Among the ILD-PH patients, 23 were classified as PAH-like and 24 as ILD-like (74% agreement with expert adjudication, Cohen’s kappa 0.49, Hedge’s g estimate of 0.86). FELLA enrichment revealed core glycolytic and bioenergetic pathways were enriched in both pure PAH and ILD-PH. Comparison of diffusion Z-scores between groups demonstrated subtle but biologically coherent differences in several reaction-level nodes, including lactate, NADPH oxidoreductase, and phosphofructokinase, which were enriched in ILD-PH. Conclusion Patients with ILD-PH demonstrate clinically relevant metabolomic phenotypes. Examination of the network biology underlying ILD-PH may point to potential therapeutic targets in this vulnerable subgroup with lagging clinical outcomes. This abstract is funded by: Pulmonary Fibrosis Foundation PFF Scholars Program

Abstract Rationale Pulmonary hypertension (PH) is a heterogeneous condition classified into 5 WSPH groups but lacking detailed molecular distinction. The PVDOMICS cohort includes all PH groups alongside healthy and disease comparators, is deeply clinically phenotyped and has RNA sequencing data. Pseudogenes are conserved genome elements that due to sequence-defects cannot be translated into functional proteins. Some pseudogenes never-the-less produce RNA-transcripts which can regulate the expression of other, protein coding genes. Therefore, our study was focused on identifying differentially expressed pseudogene-transcripts within the dataset from the PVDOMICS study and gaining insights into their potential function. Methods Peripheral venous whole blood was collected from PVDOMICS participants [n = 525 PH patients; n = 312 controls without PH who were comparators (n = 256) or healthy (n = 56)]. RNA was isolated, RNA sequencing was performed. After quality control, differential transcript expression between PH and control groups was analyzed using the limma-voom pipeline, rSEA with adjustment for age, sex, race, smoking, WSPH group, latent factors and false discovery rate. We focused on a) transcripts expressed at increased levels (n = 218): the first 200 DEG (differentially-expressed-genes)-transcripts based on log2 difference plus increased DEG within the first 200 most significant p-values; and b) transcripts expressed at decreased levels (n = 292): last 200 DEG based on log2 difference plus decreased DEG within the first 200 most significant p-values. Results Six pseudogene transcripts were among the increased and 35 pseudogene transcripts were among the decreased DEG in the PH group compared to controls. Some of the differentially expressed pseudogenes (DEP) belonged to gene families, like the ribosomal protein family (7 DEP), mitochondria expressed transcripts (3 DEP), zinc finger proteins (2 DEP). Associations of pseudogenes with biological function or disease were extracted via GeneCards(R) and information from GWAS analysis. The DUSP8P5 pseudogene was among the DEP and is known for its association with left heart disease and hypertension; 2 additional DEP have known association with hypertension. Further known associations were with inhibitory or inflammatory immune processes (8 DEP), bone density (6 DEP), metabolism and diabetes (6 DEP), respiratory function or disease (4 DEP), erythrocyte biology (3 DEP), cancer (5 DEP), kidney disease (1 DEP), and demographic parameters like eye or hair color (2 DEP). Ten DEP had no known association with function or disease. Conclusions We identified pseudogene RNA-transcripts associated with PH across WSPH groups. The local vs. systemic effects of these pseudogene transcripts on PH remain to be determined and regulators of pseudogene-transcript expression (e.g. environmental exposures) identified. This abstract is funded by: 1 X01 HL148611 (TopMed), PVDOMICS NHLBI, PVDOMICS Pulmonary Hypertension Association, HT9425-25-1-0533 (GG), 1R21OH012955-01 (GG)

Abstract Rationale Deficiency of alpha-1 antitrypsin (AAT), the protein encoded by SERPINA1, remains the strongest genetic risk factor for COPD. Individuals heterozygote for the type Z variant, termed the PI*MZ genotype (estimated global prevalence n = 35,000,000), experience varying respiratory morbidity beyond expected from tobacco smoking alone. AAT is a pleiotropic protein which inhibits not only neutrophil elastase but also neutrophil chemotaxis mediated by leukotriene B4 (LTB4) [by antagonizing the LTB4 receptor (LTB4R)] and interleukin-8 (IL-8). Although LTB-4 and IL-8 gene variants and protein levels have independently been implicated in COPD, their potential role in PI*MZ is unknown. Methods COPDGene and SPIROMICS are observational tobacco-smoke and COPD-enriched cohort studies. Participants were phenotyped by spirometry and CT imaging, with whole genome sequencing data available via TOPMed. We analyzed participants with PI*MM (COPDGene n = 9,732; SPIROMICS n = 2,105) and PI*MZ (COPDGene n = 300; SPIROMICS n = 83). LTB-4 (ALOX5AP, LTA4H, LTB4R) and IL-8 (CXCL8, CXCR1, CXCR2) associated gene variants (exonic and intronic) with a minor allele frequency of ≥ 0.05 were extracted. We assessed statistical interactions between MZ status and chest CT outcomes or FEV1% predicted using the SPAGxECCT package with adjustment for relatedness (genetic related matrix) and ancestry (top ten PCs). Baseline protein LTB4R (SPIROMICS only) and IL-8 (both cohorts) protein levels were measured via SomaScan, and statistical interactions between MZ status and AAT concentration were explored. Models were adjusted for age, sex, BMI, smoking status, and pack-years. Given the candidate gene discovery approach, significance was defined as nominal p-value<0.05 Results We identified the following number of variants in genes of the LTB-4 or IL-8 pathways: 87 shared between cohorts; 133 unique to COPDGene; 1 unique to SPIROMICS. LTB-4-related polymorphisms were associated (p < 0.05) with any radiographic outcomes in 22 variants (COPDGene) vs. 23 variants (SPIROMICS), with a single overlap (rs17074984 in ALOX5AP). Two LTB-4-pathway variants were associated with FEV1% predicted in SPIROMICS (both in ALOX5AP), but none in COPDGene. AAT plasma concentrations displayed a statistically significant interaction with FEV1% predicted (p = 0.046) and PRMemphysema (p = 0.013) for LTBR4, but not for IL-8. There was not an interaction by MM/MZ status and protein levels. Conclusions Predominantly LTB-4, and not IL-8, pathway-related genetic variants and protein levels may contribute to disease heterogeneity among PI*MZ individuals. The association of circulating LTB4R protein with reduced lung function, particularly among those with low AAT levels, further suggest interactions between AAT and the leukotriene pathway in PI*MZ individuals. Thus, PI*MZ individuals may uniquely derive benefit from leukotriene-targeting pharmacotherapy. This abstract is funded by: NIH NHLBI K23HL173570