Machine learning identifies prognostic subtypes of the tumor microenvironment of NSCLC.

The tumor microenvironment (TME) plays a fundamental role in tumorigenesis, tumor progression, and anti-cancer immunity potential of emerging cancer therapeutics. Understanding inter-patient TME heterogeneity, however, remains a challenge to efficient drug development. This article applies recent advances in machine learning (ML) for survival analysis to a retrospective study of NSCLC patients who received definitive surgical resection and immune pathology following surgery. ML methods are compared for their effectiveness in identifying prognostic subtypes. Six survival models, including Cox regression and five survival machine learning methods, were calibrated and applied to predict survival for NSCLC patients based on PD-L1 expression, CD3 expression, and ten baseline patient characteristics. Prognostic subregions of the biomarker space are delineated for each method using synthetic patient data augmentation and compared between models for overall survival concordance. A total of 423 NSCLC patients (46% female; median age [inter quantile range]: 67 [60-73]) treated with definite surgical resection were included in the study. And 219 (52%) patients experienced events during the observation period consisting of a maximum follow-up of 10 years and median follow up 78 months. The random survival forest (RSF) achieved the highest predictive accuracy, with a C-index of 0.84. The resultant biomarker subtypes demonstrate that patients with high PD-L1 expression combined with low CD3 counts experience higher risk of death within five-years of surgical resection.

Unsupervised representation learning on high-dimensional clinical data improves genomic discovery and prediction.

Although high-dimensional clinical data (HDCD) are increasingly available in biobank-scale datasets, their use for genetic discovery remains challenging. Here we introduce an unsupervised deep learning model, Representation Learning for Genetic Discovery on Low-Dimensional Embeddings (REGLE), for discovering associations between genetic variants and HDCD. REGLE leverages variational autoencoders to compute nonlinear disentangled embeddings of HDCD, which become the inputs to genome-wide association studies (GWAS). REGLE can uncover features not captured by existing expert-defined features and enables the creation of accurate disease-specific polygenic risk scores (PRSs) in datasets with very few labeled data. We apply REGLE to perform GWAS on respiratory and circulatory HDCD-spirograms measuring lung function and photoplethysmograms measuring blood volume changes. REGLE replicates known loci while identifying others not previously detected. REGLE are predictive of overall survival, and PRSs constructed from REGLE loci improve disease prediction across multiple biobanks. Overall, REGLE contain clinically relevant information beyond that captured by existing expert-defined features, leading to improved genetic discovery and disease prediction.

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.

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.

Comparing 90-Day Postoperative Mortality After Neoadjuvant Proton-Based Versus Photon-Based Chemoradiotherapy for Esophageal Cancer.

Purpose: Evidence suggests that proton-beam therapy (PBT) results in less toxicity and postoperative complications compared to photon-based radiotherapy in patients who receive chemoradiotherapy followed by esophagectomy for cancer. Ninety-day mortality (90DM) is an important measure of the postoperative (nononcologic) outcome as proxy of quality-of-care. We hypothesize that PBT could reduce 90DM compared to photon-based radiotherapy. Materials and Methods: From a single-center retrospective database patients treated with chemoradiotherapy before esophagectomy for cancer were selected (1998-2022). Univariable logistic regression was used to study the association of radiotherapy modality with 90DM. Three separate methods were applied to adjust for confounding bias, including multivariable logistic regression, propensity score matching, and inverse probability of treatment weighting. Stratified analysis for the age threshold that maximized the difference in 90DM (ie, ≥67 vs <67 years) was performed. Results: A total of 894 eligible patients were included and 90DM was 5/202 (2.5%) in the PBT versus 29/692 (4.2%) in the photon-based radiotherapy group (P = .262). After adjustment for age and tumor location, PBT versus photon-based radiotherapy was not significantly associated with 90DM (P = .491). The 90DM was not significantly different for PBT versus photon-based radiotherapy in the propensity score matching (P = .379) and inverse probability of treatment weighting cohort (P = .426). The stratified analysis revealed that in patients aged ≥67 years, PBT was associated with decreased 90DM (1.3% vs 8.8%; P = .026). Higher age significantly increased 90DM risk within the photon-based radiotherapy (8.8% vs 2.7%; P = .001), but not within the PBT group (1.3% vs 3.2%; P = .651). Conclusion: No statistically significant difference was observed in postoperative 90DM after esophagectomy for cancer between PBT and photon-based neoadjuvant chemoradiotherapy. However, among older patients a signal was observed that PBT may reduce 90DM risk.

Personalized Composite Dosimetric Score-Based Machine Learning Model of Severe Radiation-Induced Lymphopenia Among Patients With Esophageal Cancer.

PURPOSE: Radiation-induced lymphopenia (RIL) is common among patients undergoing radiation therapy (RT)' Severe RIL has been linked to adverse outcomes. The severity and risk of RIL can be predicted from baseline clinical characteristics and dosimetric parameters. However, dosimetric parameters, e.g. dose-volume (DV) indices, are highly correlated with one another and are only weakly associated with RIL. Here we introduce the novel concept of "composite dosimetric score" (CDS) as the index that condenses the dose distribution in immune tissues of interest to study the dosimetric dependence of RIL. We derived an improved multivariate classification scheme for risk of grade 4 RIL (G4RIL), based on this novel RT dosimetric feature, for patients receiving chemo RT for esophageal cancer. METHODS AND MATERIALS: DV indices were extracted for 734 patients who received chemo RT for biopsy-proven esophageal cancer. Nonnegative matrix factorization was used to project the DV indices of lung, heart, and spleen into a single CDS; XGBoost was employed to explore significant interactions among predictors; and logistic regression was applied to combine the resultant CDS with baseline clinical factors and interaction terms to facilitate individualized prediction of immunotoxicity. Five-fold cross-validation was applied to evaluate the model performance. RESULTS: The CDS for selected immune organs at risk (ie, heart, lung, and spleen) (OR 1.791; 95 CI [1.350, 2.377]) was a statistically significant risk determinant for G4RIL. Pearson correlation coefficients for CDS versus G4RIL risk for individual immune organs at risk were greater than any single DV indicx. Personalized prediction of G4RIL based on CDS and 4 clinical risk factors yielded an area under the curve value of 0.78. Interaction between age and CDS revealed that G4RIL risk increased more sharply with increasing CDS for patients aged ≥65 years. CONCLUSIONS: Risk of immunotoxicity for patients undergoing chemo RT for esophageal cancer can be predicted by CDS. The CDS concept can be extended to immunotoxicity in other cancer types and in dose-response models currently based on DV indices. Personalized treatment planning should leverage composite dosimetric scoring methods rather than using individual or subsets of DV indices.

Clinical factors associated with racial differences in the prevalence of occult hypoxemia: a retrospective case-control study.

Background: Recent studies showed that Black patients more often have falsely normal oxygen saturation on pulse oximetry compared to White patients. However, whether the racial differences in occult hypoxemia are mediated by other clinical differences is unknown. Methods: We conducted a retrospective case-control study utilizing two large ICU databases (eICU and MIMIC-IV). We defined occult hypoxemia as oxygen saturation on pulse oximetry within 92-98% despite oxygen saturation on arterial blood gas below 90%. We assessed associations of commonly measured clinical factors with occult hypoxemia using multivariable logistic regression and conducted mediation analysis of the racial effect. Results: Among 24,641 patients, there were 1,855 occult hypoxemia cases and 23,786 controls. In both datasets, Black patients were more likely to have occult hypoxemia (unadjusted odds ratio 1.66 [95%-CI: 1.41-1.95] in eICU and 2.00 [95%-CI: 1.22-3.14] in MIMIC-IV). In multivariable models, higher respiratory rate, PaCO2 and creatinine as well as lower hemoglobin were associated with increased odds of occult hypoxemia. Differences in the commonly measured clinical markers accounted for 9.2% and 44.4% of the racial effect on occult hypoxemia in eICU and MIMIC-IV, respectively. Conclusion: Clinical differences, in addition to skin tone, might mediate some of the racial differences in occult hypoxemia.

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.

Severe Lymphopenia During Chemoradiation Therapy for Esophageal Cancer: Comprehensive Analysis of Randomized Phase 2B Trial of Proton Beam Therapy Versus Intensity Modulated Radiation Therapy.

PURPOSE: Lymphocytes play an important role in antitumor immunity; however, they are also especially vulnerable to depletion during chemoradiation therapy (CRT). The purpose of this study was to compare the incidence of grade 4 lymphopenia (G4L) between proton beam therapy (PBT) and intensity modulated photon radiation therapy (IMRT) in patients with esophageal cancer treated with CRT in a completed randomized trial and to ascertain patient heterogeneity to G4L risk based on treatment and established prognostic factors. METHODS AND MATERIALS: Between April 2012 and March 2019, a single-institution, open-label, nonblinded, phase 2 randomized trial (NCT01512589) was conducted at the University of Texas MD Anderson Cancer Center. Patients were randomly assigned to IMRT or PBT, either definitively or preoperatively. This secondary analysis of the randomized trial was G4L during concurrent CRT according to Common Terminology Criteria for Adverse Events version 5.0. RESULTS: Among 105 patients evaluable for analysis, 44 patients (42%) experienced G4L at a median of 28 days after the start date of concurrent CRT. Induction chemotherapy (P = .003), baseline absolute lymphocyte count (P < .001), radiation therapy modality (P = .002), and planning treatment volume (P = .033) were found to be significantly associated with G4L. Multivariate classification analysis partitioned patients into 5 subgroups for whom the incidence of G4L was observed in 0%, 14%, 35%, 70%, and 100% of patients. The benefit of PBT over IMRT was most pronounced in patients with an intermediate baseline absolute lymphocyte count and large planning treatment volume (P = .011). CONCLUSIONS: This is the first prospective evidence that limiting dose scatter by PBT significantly reduced the incidence of G4L, especially in the intermediate-risk patients. The implication of this immune-sparing effect of PBT, especially in the context of standard adjuvant immunotherapy, needs further examination in the current phase 3 randomized trials.

Plasma metabolomics and quantitative interstitial abnormalities in ever-smokers.

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

Clonal Somatic Mutations in Chronic Lung Diseases Are Associated with Reduced Lung Function.

Rationale: Constantly exposed to the external environment and mutagens such as tobacco smoke, human lungs have one of the highest somatic mutation rates among all human organs. However, the relationship of these mutations to lung disease and function is not known. Objectives: To identify the prevalence and significance of clonal somatic mutations in chronic lung diseases. Methods: We analyzed the clonal somatic mutations from 1,251 samples of normal and diseased noncancerous lung tissue RNA sequencing with paired whole-genome sequencing from the Lung Tissue Research Consortium. We examined the associations of somatic mutations with lung function, disease status, and computationally deconvoluted cell types in two of the most common diseases represented in our dataset, chronic obstructive pulmonary disease (COPD; 29%) and idiopathic pulmonary fibrosis (IPF; 13%). Measurements and Main Results: Clonal somatic mutational burden was associated with reduced lung function in both COPD and IPF. We identified an increased prevalence of clonal somatic mutations in individuals with IPF compared with normal control subjects and individuals with COPD independent of age and smoking status. IPF clonal somatic mutations were enriched in disease-related and airway epithelial-expressed genes such as MUC5B in IPF. Patients who were MUC5B risk variant carriers had increased odds of developing somatic mutations of MUC5B that were explained by increased expression of MUC5B. Conclusions: Our identification of an increased prevalence of clonal somatic mutation in diseased lung that correlates with airway epithelial gene expression and disease severity highlights for the first time the role of somatic mutational processes in lung disease genetics.

WHOLE GENOME SEQUENCING ANALYSIS OF BODY MASS INDEX IDENTIFIES NOVEL AFRICAN ANCESTRY-SPECIFIC RISK ALLELE.

Obesity is a major public health crisis associated with high mortality rates. Previous genome-wide association studies (GWAS) investigating body mass index (BMI) have largely relied on imputed data from European individuals. This study leveraged whole-genome sequencing (WGS) data from 88,873 participants from the Trans-Omics for Precision Medicine (TOPMed) Program, of which 51% were of non-European population groups. We discovered 18 BMI-associated signals (P < 5 × 10-9). Notably, we identified and replicated a novel low frequency single nucleotide polymorphism (SNP) in MTMR3 that was common in individuals of African descent. Using a diverse study population, we further identified two novel secondary signals in known BMI loci and pinpointed two likely causal variants in the POC5 and DMD loci. Our work demonstrates the benefits of combining WGS and diverse cohorts in expanding current catalog of variants and genes confer risk for obesity, bringing us one step closer to personalized medicine.

Augmentation Therapy for Severe Alpha-1 Antitrypsin Deficiency Improves Survival and Is Decoupled from Spirometric Decline-A Multinational Registry Analysis.

Rationale: Intravenous plasma-purified alpha-1 antitrypsin (IV-AAT) has been used as therapy for alpha-1 antitrypsin deficiency (AATD) since 1987. Previous trials (RAPID and RAPID-OLE) demonstrated efficacy in preserving computed tomography of lung density but no effect on FEV1. This observational study evaluated 615 people with severe AATD from three countries with socialized health care (Ireland, Switzerland, and Austria), where access to standard medical care was equal but access to IV-AAT was not. Objectives: To assess the real-world longitudinal effects of IV-AAT. Methods: Pulmonary function and mortality data were utilized to perform longitudinal analyses on registry participants with severe AATD. Measurements and Main Results: IV-AAT confers a survival benefit in severe AATD (P < 0.001). We uncovered two distinct AATD phenotypes based on an initial respiratory diagnosis: lung index and non-lung index. Lung indexes demonstrated a more rapid FEV1 decline between the ages of 20 and 50 and subsequently entered a plateau phase of minimal decline from 50 onward. Consequentially, IV-AAT had no effect on FEV1 decline, except in patients with a Global Initiative for Chronic Obstructive Lung Disease (GOLD) stage 2 lung index. Conclusions: This real-world study demonstrates a survival advantage from IV-AAT. This improved survival is largely decoupled from FEV1 decline. The observation that patients with severe AATD fall into two major phenotypes has implications for clinical trial design where FEV1 is a primary endpoint. Recruits into trials are typically older lung indexes entering the plateau phase and, therefore, unlikely to show spirometric benefits. IV-AAT attenuates spirometric decline in lung indexes in GOLD stage 2, a spirometric group commonly outside current IV-AAT commencement recommendations.

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.

Declarations of Independence: How Embedded Multicollinearity Errors Affect Dosimetric and Other Complex Analyses in Radiation Oncology.

The statistical technique of multiple regression, commonly referred to as "multivariable regression," is often used in clinical research to quantify the relationships between multiple predictor variables and a single outcome variable of interest. The foundational theory underpinning multivariable regression assumes that all predictor variables are independent of one another. In other words, the effect of each independent variable is measured by its contribution to the regression equation while all other variables remain unchanged. In the presence of correlations between two or more variables, however, it is impossible to change one variable without a consequent change in the variable(s) it is linked to. This condition, known as "multicollinearity," can introduce errors into multivariable regression models by affecting estimates of the regression coefficients that quantify the relationship between individual predictor variables and the outcome variable. Errors that arise due to violations of the multicollinearity assumption are of special interest to radiation oncology researchers. Because of high levels of correlation among variables derived from points along individual organ dose-volume histogram (DVH) curves, as well as strong intercorrelations among dose-volume parameters in neighboring organs, dosimetric analyses are particularly subject to multicollinearity errors. For example, dose-volume parameters for the heart are strongly correlated not only with other points along the heart DVH curve but are likely also correlated with dose-volume parameters in neighboring organs such as the lung. In this paper, we describe the problem of multicollinearity in accessible terms and discuss examples of violations of the multicollinearity assumption within the radiation oncology literature. Finally, we provide recommendations regarding best practices for identifying and managing multicollinearity in complex data sets.

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.

Large scale proteomic studies create novel privacy considerations.

Privacy protection is a core principle of genomic but not proteomic research. We identified independent single nucleotide polymorphism (SNP) quantitative trait loci (pQTL) from COPDGene and Jackson Heart Study (JHS), calculated continuous protein level genotype probabilities, and then applied a naïve Bayesian approach to link SomaScan 1.3K proteomes to genomes for 2812 independent subjects from COPDGene, JHS, SubPopulations and InteRmediate Outcome Measures In COPD Study (SPIROMICS) and Multi-Ethnic Study of Atherosclerosis (MESA). We correctly linked 90-95% of proteomes to their correct genome and for 95-99% we identify the 1% most likely links. The linking accuracy in subjects with African ancestry was lower (~ 60%) unless training included diverse subjects. With larger profiling (SomaScan 5K) in the Atherosclerosis Risk Communities (ARIC) correct identification was > 99% even in mixed ancestry populations. We also linked proteomes-to-proteomes and used the proteome only to determine features such as sex, ancestry, and first-degree relatives. When serial proteomes are available, the linking algorithm can be used to identify and correct mislabeled samples. This work also demonstrates the importance of including diverse populations in omics research and that large proteomic datasets (> 1000 proteins) can be accurately linked to a specific genome through pQTL knowledge and should not be considered unidentifiable.

An ERK5-NRF2 Axis Mediates Senescence-Associated Stemness and Atherosclerosis.

BACKGROUND: ERK5 (extracellular signal-regulated kinase 5) is a dual kinase transcription factor containing an N-terminal kinase domain and a C-terminal transcriptional activation domain. Many ERK5 kinase inhibitors have been developed and tested to treat cancer and inflammatory diseases. However, recent data have raised questions about the role of the catalytic activity of ERK5 in proliferation and inflammation. We aimed to investigate how ERK5 reprograms myeloid cells to the proinflammatory senescent phenotype, subsequently leading to atherosclerosis. METHODS: A ERK5 S496A (dephosphorylation mimic) knock in (KI) mouse model was generated using CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats/clustered regularly interspaced short palindromic repeat-associated 9), and atherosclerosis was characterized by hypercholesterolemia induction. The plaque phenotyping in homozygous ERK5 S496A KI and wild type (WT) mice was studied using imaging mass cytometry. Bone marrow-derived macrophages were isolated from hypercholesterolemic mice and characterized using RNA sequencing and functional in vitro approaches, including senescence, mitochondria reactive oxygen species, and inflammation assays, as well as by metabolic extracellular flux analysis. RESULTS: We show that atherosclerosis was inhibited in ERK5 S496A KI mice. Furthermore, ERK5 S496 phosphorylation mediates both senescence-associated secretory phenotype and senescence-associated stemness by upregulating AHR (aryl hydrocarbon receptor) in plaque and bone marrow-derived macrophages isolated from hypercholesterolemic mice. We also discovered that ERK5 S496 phosphorylation could induce NRF2 (NFE2-related factor 2) SUMOylation at a novel K518 site to inhibit NRF2 transcriptional activity without altering ERK5 catalytic activity and mediates oxidized LDL (low-density lipoprotein)-induced senescence-associated secretory phenotype. Specific ERK5 kinase inhibitors (AX15836 and XMD8-92) also inhibited ERK5 S496 phosphorylation, suggesting the involvement of ERK5 S496 phosphorylation in the anti-inflammatory effects of these ERK5 kinase inhibitors. CONCLUSIONS: We discovered a novel mechanism by which the macrophage ERK5-NRF2 axis develops a unique senescence-associated secretory phenotype/stemness phenotype by upregulating AHR to engender atherogenesis. The finding of senescence-associated stemness phenotype provides a molecular explanation to resolve the paradox of senescence in proliferative plaque by permitting myeloid cells to escape the senescence-induced cell cycle arrest during atherosclerosis formation.