Genetic diversity fuels gene discovery for tobacco and alcohol use.

Tobacco and alcohol use are heritable behaviours associated with 15% and 5.3% of worldwide deaths, respectively, due largely to broad increased risk for disease and injury1-4. These substances are used across the globe, yet genome-wide association studies have focused largely on individuals of European ancestries5. Here we leveraged global genetic diversity across 3.4 million individuals from four major clines of global ancestry (approximately 21% non-European) to power the discovery and fine-mapping of genomic loci associated with tobacco and alcohol use, to inform function of these loci via ancestry-aware transcriptome-wide association studies, and to evaluate the genetic architecture and predictive power of polygenic risk within and across populations. We found that increases in sample size and genetic diversity improved locus identification and fine-mapping resolution, and that a large majority of the 3,823 associated variants (from 2,143 loci) showed consistent effect sizes across ancestry dimensions. However, polygenic risk scores developed in one ancestry performed poorly in others, highlighting the continued need to increase sample sizes of diverse ancestries to realize any potential benefit of polygenic prediction.

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.

Loss of Airway Phylogenetic Diversity Is Associated with Clinical and Pathobiological Markers of Disease Development in COPD.

RATIONALE: The airway microbiome has the potential to shape COPD pathogenesis, but its relationship to outcomes in milder disease is unestablished. OBJECTIVES: Identify sputum microbiome characteristics associated with markers of COPD in participants of the SubPopulations and InteRmediate Outcome Measures of COPD Study (SPIROMICS). METHODS: Sputum DNA from 877 participants were analyzed using 16S rRNA gene sequencing. Relationships between baseline airway microbiota composition and clinical, radiographic and muco-inflammatory markers, including longitudinal lung function trajectory, were examined. MEASUREMENTS AND MAIN RESULTS: Participant data represented predominantly milder disease (GOLD 0-2: N=732/877). Phylogenetic diversity (range of different species within a sample) correlated positively with baseline lung function, declined with higher GOLD stage, and correlated negatively with symptom burden, radiographic markers of airway disease, and total mucin concentrations (p<0.001). In co-variate adjusted regression models, organisms robustly associated with better lung function included members of Alloprevotella, Oribacterium, and Veillonella. Conversely, lower lung function, greater symptoms and radiographic measures of small airway disease associated with enrichment in members of Streptococcus, Actinobacillus, Actinomyces, and other genera. Baseline sputum microbiota features also associated with lung function trajectory during SPIROMICS follow up (stable/improved, decliner, or rapid decliner). The 'stable/improved' group (slope of FEV1 regression ≥ 66th percentile) had higher bacterial diversity at baseline, associated with enrichment in Prevotella, Leptotrichia, and Neisseria. In contrast, the 'rapid decliner' group (FEV1 slope ≤ 33rd percentile) had significantly lower baseline diversity, associated with enrichment in Streptococcus. CONCLUSIONS: In SPIROMICS baseline airway microbiota features demonstrate divergent associations with better or worse COPD-related outcomes.

Dysanapsis Genetic Risk Predicts Lung Function Across the Lifespan.

Rationale Dysanapsis refers to a mismatch between airway tree caliber and lung size arising early in life. Dysanapsis assessed by computed tomography (CT) is evident by early adulthood and associated with chronic obstructive pulmonary disease (COPD) risk later in life. Objective By examining the genetic factors associated with CT-assessed dysanapsis, we aimed to elucidate its molecular underpinnings and physiological significance across the lifespan. Methods We performed a genome-wide association study (GWAS) of CT-assessed dysanapsis in 11,951 adults, including individuals from two population-based and two COPD-enriched studies. We applied colocalization analysis to integrate GWAS and gene expression data from whole blood and lung. Genetic variants associated with dysanapsis were combined into a genetic risk score that was applied to examine association with lung function in children from a population-based birth cohort (n=1,278) and adults from the UK Biobank (n=369,157). Measurements and Main Results CT-assessed dysanapsis was associated with genetic variants from 21 independent signals in 19 gene regions, implicating HHIP, DSP, and NPNT as potential molecular targets based on colocalization of their expression. Higher dysanapsis genetic risk score was associated with obstructive spirometry among 5 year old children and among adults in the 5th, 6th and 7th decades of life. Conclusions CT-assessed dysanapsis is associated with variation in genes previously implicated in lung development and dysanapsis genetic risk is associated with obstructive lung function from early life through older adulthood. Dysanapsis may represent an endo-phenotype link between the genetic variations associated with lung function and COPD.

Utility of eosinophil peroxidase as a biomarker of eosinophilic inflammation in asthma.

BACKGROUND: The relative utility of eosinophil peroxidase (EPX) and blood and sputum eosinophil counts as disease biomarkers in asthma is uncertain. OBJECTIVE: We sought to determine the utility of EPX as a biomarker of systemic and airway eosinophilic inflammation in asthma. METHODS: EPX protein was measured by immunoassay in serum and sputum in 110 healthy controls to establish a normal reference range and in repeated samples of serum and sputum collected during 3 years of observation in 480 participants in the Severe Asthma Research Program 3. RESULTS: Over 3 years, EPX levels in patients with asthma were higher than normal in 27% to 31% of serum samples and 36% to 53% of sputum samples. Eosinophils and EPX correlated better in blood than in sputum (rs values of 0.74 and 0.43, respectively), and high sputum EPX levels occurred in 27% of participants with blood eosinophil counts less than 150 cells/µL and 42% of participants with blood eosinophil counts between 150 and 299 cells/µL. Patients with persistently high sputum EPX values for 3 years were characterized by severe airflow obstruction, frequent exacerbations, and high mucus plug scores. In 59 patients with asthma who started mepolizumab during observation, serum EPX levels normalized in 96% but sputum EPX normalized in only 49%. Lung function remained abnormal even when sputum EPX normalized. CONCLUSIONS: Serum EPX is a valid protein biomarker of systemic eosinophilic inflammation in asthma, and sputum EPX levels are a more sensitive biomarker of airway eosinophilic inflammation than sputum eosinophil counts. Eosinophil measures in blood frequently miss airway eosinophilic inflammation, and mepolizumab frequently fails to normalize airway eosinophilic inflammation even though it invariably normalizes systemic eosinophilic inflammation.

A novel DNase assay reveals low DNase activity in severe asthma.

Secreted deoxyribonucleases (DNases), such as DNase-I and DNase-IL3, degrade extracellular DNA, and endogenous DNases have roles in resolving airway inflammation and guarding against autoimmune responses to nucleotides. Subsets of patients with asthma have high airway DNA levels, but information about DNase activity in health and in asthma is lacking. To characterize DNase activity in health and in asthma, we developed a novel kinetic assay using a Taqman probe sequence that is quickly cleaved by DNase-I to produce a large product signal. We used this kinetic assay to measure DNase activity in sputum from participants in the Severe Asthma Research Program (SARP)-3 (n = 439) and from healthy controls (n = 89). We found that DNase activity was lower than normal in asthma [78.7 relative fluorescence units (RFU)/min vs. 120.4 RFU/min, P < 0.0001]. Compared to patients with asthma with sputum DNase activity in the upper tertile activity levels, those in the lower tertile of sputum DNase activity were characterized clinically by more severe disease and pathologically by airway eosinophilia and airway mucus plugging. Carbamylation of DNase-I, a post-translational modification that can be mediated by eosinophil peroxidase, inactivated DNase-I. In summary, a Taqman probe-based DNase activity assay uncovers low DNase activity in the asthma airway that is associated with more severe disease and airway mucus plugging and may be caused, at least in part, by eosinophil-mediated carbamylation.NEW & NOTEWORTHY We developed a new DNase assay and used it to show that DNase activity is impaired in asthma airways.

Genetic and T2 biomarkers linked to the efficacy of HDM sublingual immunotherapy in asthma.

BACKGROUND: Hypersensitivity to house dust mite (HDM) allergens is a common cause of allergic asthma symptoms and can be effectively treated with allergy immunotherapy (AIT). OBJECTIVE: To investigate whether genetic and type 2 (T2) inflammatory biomarkers correlate with disease severity in subjects with allergic asthma, and whether this can be modified by AIT. METHODS: MITRA (NCT01433523) was a phase III, randomised, double-blind, placebo-controlled trial of HDM sublingual immunotherapy (SLIT)-tablets in adults with HDM allergic asthma. Post hoc analyses of the study population (N=742) evaluated associations between T2 inflammatory (blood eosinophils, eosinophil cationic protein (ECP), total IgE and tryptase) and genetic (single-nucleotide polymorphisms, SNP) biomarkers (n=582) for the primary study endpoint (time to first moderate/severe asthma exacerbation). SNP associations were verified in HDM-positive subgroup from an independent 3-year Severe Asthma Research Programme (SARP3) subject cohort. RESULTS: An increased asthma exacerbation risk in subjects homozygous for SNP rs7216389 (chromosomal locus 17q12-21) was reduced (p=0.037) by treatment with HDM SLIT (HR=0.37 (95% CI 0.22 to 0.64), p<0.001). The associations between exacerbation risk and 17q12-21 SNPs were replicated in the SARP3 HDM-positive subgroup. High levels of T2 biomarkers were associated with increased risk of asthma exacerbations in the placebo group. HDM SLIT-tablet treatment reduced this risk (blood eosinophils: HR=0.50 (95% CI 0.30 to 0.85); ECP: HR=0.45 (95% CI 0.29 to 0.87); tryptase: HR=0.45 (95% CI 0.25 to 0.80)). The treatment effect was higher (p=0.006) for subjects with a higher number of elevated T2 biomarkers. CONCLUSIONS: HDM SLIT-tablet AIT is efficacious in HDM-sensitised asthma subjects with a genetic asthma predisposition and/or an underlying T2 endotype. TRIAL REGISTRATION NUMBER: NCT01433523.

Protein-Protein interactive networks identified in bronchoalveolar lavage of severe compared to nonsevere asthma.

INTRODUCTION: Previous bronchoalveolar lavage fluid (BALF) proteomic analysis has evaluated limited numbers of subjects for only a few proteins of interest, which may differ between asthma and normal controls. Our objective was to examine a more comprehensive inflammatory biomarker panel in quantitative proteomic analysis for a large asthma cohort to identify molecular phenotypes distinguishing severe from nonsevere asthma. METHODS: Bronchoalveolar lavage fluid from 48 severe and 77 nonsevere adult asthma subjects were assessed for 75 inflammatory proteins, normalized to BALF total protein concentration. Validation of BALF differences was sought through equivalent protein analysis of autologous sputum. Subjects' data, stratified by asthma severity, were analysed by standard statistical tests, principal component analysis and 5 machine learning algorithms. RESULTS: The severe group had lower lung function and greater health care utilization. Significantly increased BALF proteins for severe asthma compared to nonsevere asthma were fibroblast growth factor 2 (FGF2), TGFα, IL1Ra, IL2, IL4, CCL8, CCL13 and CXCL7 and significantly decreased were platelet-derived growth factor a-a dimer (PDGFaa), vascular endothelial growth factor (VEGF), interleukin 5 (IL5), CCL17, CCL22, CXCL9 and CXCL10. Four protein differences were replicated in sputum. FGF2, PDGFaa and CXCL7 were independently identified by 5 machine learning algorithms as the most important variables for discriminating severe and nonsevere asthma. Increased and decreased proteins identified for the severe cluster showed significant protein-protein interactions for chemokine and cytokine signalling, growth factor activity, and eosinophil and neutrophil chemotaxis differing between subjects with severe and nonsevere asthma. CONCLUSION: These inflammatory protein results confirm altered airway remodelling and cytokine/chemokine activity recruiting leukocytes into the airways of severe compared to nonsevere asthma as important processes even in stable status.

Asthma control in the United States: relationships between short-acting ß2-agonist and systemic corticosteroid use.

BACKGROUND: Asthma control assessment is based on impairment (current symptoms) and risk (exacerbation history). OBJECTIVE: To understand the extent of uncontrolled asthma, we assessed relationships between prescription fills for systemic corticosteroids (SCS) and short-acting ß2-agonists (SABA) as risk and impairment markers, respectively. METHODS: Annual SCS and SABA fills among US patients with asthma were evaluated by a retrospective analysis of IQVIA Longitudinal Access and Adjudication Data. Patient severity was assigned based on GINA step-therapy level. Exacerbations were evaluated by SCS fills within 12 months of a first asthma prescription fill. Uncontrolled asthma was defined as ≥2 SCS and/or ≥3 SABA fills annually. Individual patient relationships between SCS and SABA fills were assessed by Pearson's correlation coefficient. RESULTS: 4,506,527 patients were included: 15% had ≥2 SCS fills, 29% had ≥3 SABA fills, 37% fulfilled either or both criteria. If only SCS were assessed, 22% treated as mild-to-moderate and 27% as severe asthma would have been misclassified as controlled. If only SABA use was evaluated, 8% treated as mild-to-moderate and 11% as severe asthma would have been misclassified. Overall, 81% of uncontrolled asthma occurred in patients treated for mild-to-moderate disease. Among patients with ≥2 SCS fills, mean SABA fills were 2.9; the correlation between SCS and SABA fills per patient was significant but weak (r=0.18; p<0.001). CONCLUSION: High symptom burden and SCS exposures are not limited to severe asthma but are also characteristic in patients treated for mild-to-moderate disease. Both impairment and risk assessments are required to understand the full extent of uncontrolled asthma across disease severities.

Real-world effectiveness of dupilumab in patients with asthma: Findings from the US ADVANTAGE study.

BACKGROUND: Dupilumab is approved as an add-on maintenance therapy for patients (≥6 years) with moderate-to-severe asthma. Better understanding of real-world effectiveness is needed. OBJECTIVE: To characterize the real-world effectiveness of dupilumab in asthma management. METHODS: This retrospective study included patients (≥12 years of age) diagnosed with asthma, initiating dupilumab between November 2018 and September 2020. The study used a US electronic medical record database (TriNetX Dataworks, Cambridge, Massachusetts). Asthma exacerbation rates before and after the initiation of dupilumab were analyzed using generalized estimating equations models with Poisson probabilistic link to estimate incidence rate ratios (IRRs). Sensitivity analyses were conducted based on previous exacerbation data, eosinophil levels, history of atopic dermatitis or chronic rhinosinusitis with nasal polyps, previous use of biologics, and presence of SARS-CoV-2 (COVID-19). RESULTS: A total of 2400 patients initiating dupilumab met all study criteria. After initiation of dupilumab, risk of asthma exacerbation was reduced by 44% (IRR, 0.56; 95% CI, 0.47-0.57; P = <0.0001) and systemic corticosteroid prescriptions by 48% (IRR, 0.52; 95% CI, 0.48, 0.56; P = <0.0001) compared with those before initiation of dupilumab. Adjustment for COVID-19 showed a greater reduction in asthma exacerbations (IRR, 0.50; 95% CI, 0.45-0.55; P = <0.0001). CONCLUSION: Current real-world efficacy evidence indicates that dupilumab reduces asthma exacerbations and total systemic corticosteroid prescriptions in clinical practice. The effectiveness of dupilumab was observed independent of exacerbation history, eosinophil levels, or COVID-19 impact.

Low CC16 mRNA Expression Levels in Bronchial Epithelial Cells Are Associated with Asthma Severity.

Rationale: CC16 is a protein mainly produced by nonciliated bronchial epithelial cells (BECs) that participates in host defense. Reduced CC16 protein concentrations in BAL and serum are associated with asthma susceptibility. Objectives: Few studies have investigated the relationship between CC16 and asthma progression, and none has focused on BECs. In this study, we sought to determine if CC16 mRNA expression levels in BECs are associated with asthma severity. Methods: Association analyses between CC16 mRNA expression levels in BECs (242 asthmatics and 69 control subjects) and asthma-related phenotypes in Severe Asthma Research Program were performed using a generalized linear model. Measurements and Main Results: Low CC16 mRNA expression levels in BECs were significantly associated with asthma susceptibility and asthma severity, high systemic corticosteroids use, high retrospective and prospective asthma exacerbations, and low pulmonary function. Low CC16 mRNA expression levels were significantly associated with high T2 inflammation biomarkers (fractional exhaled nitric oxide and sputum eosinophils). CC16 mRNA expression levels were negatively correlated with expression levels of Th2 genes (IL1RL1, POSTN, SERPINB2, CLCA1, NOS2, and MUC5AC) and positively correlated with expression levels of Th1 and inflammation genes (IL12A and MUC5B). A combination of two nontraditional T2 biomarkers (CC16 and IL-6) revealed four asthma endotypes with different characteristics of T2 inflammation, obesity, and asthma severity. Conclusions: Our findings indicate that low CC16 mRNA expression levels in BECs are associated with asthma susceptibility, severity, and exacerbations, partially through immunomodulation of T2 inflammation. CC16 is a potential nontraditional T2 biomarker for asthma development and progression.

Skeletal Muscle Adiposity and Lung Function Trajectory in the Severe Asthma Research Program.

Rationale: Extrapulmonary manifestations of asthma, including fatty infiltration in tissues, may reflect systemic inflammation and influence lung function and disease severity. Objectives: To determine if skeletal muscle adiposity predicts lung function trajectory in asthma. Methods: Adult SARP III (Severe Asthma Research Program III) participants with baseline computed tomography imaging and longitudinal postbronchodilator FEV1% predicted (median follow-up 5 years [1,132 person-years]) were evaluated. The mean of left and right paraspinous muscle density (PSMD) at the 12th thoracic vertebral body was calculated (Hounsfield units [HU]). Lower PSMD reflects higher muscle adiposity. We derived PSMD reference ranges from healthy control subjects without asthma. A linear multivariable mixed-effects model was constructed to evaluate associations of baseline PSMD and lung function trajectory stratified by sex. Measurements and Main Results: Participants included 219 with asthma (67% women; mean [SD] body mass index, 32.3 [8.8] kg/m2) and 37 control subjects (51% women; mean [SD] body mass index, 26.3 [4.7] kg/m2). Participants with asthma had lower adjusted PSMD than control subjects (42.2 vs. 55.8 HU; P < 0.001). In adjusted models, PSMD predicted lung function trajectory in women with asthma (ß = -0.47 Δ slope per 10-HU decrease; P = 0.03) but not men (ß = 0.11 Δ slope per 10-HU decrease; P = 0.77). The highest PSMD tertile predicted a 2.9% improvement whereas the lowest tertile predicted a 1.8% decline in FEV1% predicted among women with asthma over 5 years. Conclusions: Participants with asthma have lower PSMD, reflecting greater muscle fat infiltration. Baseline PSMD predicted lung function decline among women with asthma but not men. These data support an important role of metabolic dysfunction in lung function decline.

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.

Real-world effectiveness of dupilumab versus benralizumab and mepolizumab.

INTRODUCTION: In the United States, this real-world study compared the effectiveness of dupilumab, benralizumab, andmepolizumab in reducing exacerbations and systemic corticosteroid (SCS) prescriptions among patients with asthma. METHODS: Patients (12 years old) who initiated dupilumab, benralizumab, or mepolizumab (index) between November 2018 and September 2020 were identified by using electronic medical record data. Subjects were included if they had greater than or equal to 12 months of data before and after the index date and two or more severe asthma-related exacerbations before the index date. Differences in baseline characteristics were addressed by using inverse probability treatment weighting (IPTW). Pairwisecomparisons between dupilumab and benralizumab, or mepolizumab were conducted by using negative binomial regression, adjusting for baseline rates and unbalance characteristics (greater than or equal to 10% standardized differences) after IPTW. RESULTS: Overall, a total of 1737 subjects met all criteria: 825 dupilumab, 461 benralizumab, and 451 mepolizumab initiators.In the postindex period, dupilumab was associated with a 24% and 28% significant reduction in the risk of severe asthmaexacerbations versus benralizumab (incidence rate ratio [IRR] 0.76 [95% confidence interval {CI}, 0.67-0.86)] and mepolizumab(IRR 0.72 [95% CI, 0.63-0.82]), respectively. In addition, dupilumab treatment significantly reduced SCS prescriptionsby 16% and 25% versus benralizumab and mepolizumab, respectively (p < 0.05). CONCLUSION: This study represents one of the largest real-world comparisons of biologics (dupilumab, benralizumab, and mepolizumab) for asthma in the United States to date. This analysis shows that the use of dupilumab was associated with a significantly greater reduction in both severe asthma exacerbations and SCS prescriptions compared with benralizumab and mepolizumab.

Bronchial epithelial cell transcriptional responses to inhaled corticosteroids dictate severe asthmatic outcomes.

BACKGROUND: Inhaled corticosteroids (CSs) are the backbone of asthma treatment, improving quality of life, exacerbation rates, and mortality. Although effective for most, a subset of patients with asthma experience CS-resistant disease despite receiving high-dose medication. OBJECTIVE: We sought to investigate the transcriptomic response of bronchial epithelial cells (BECs) to inhaled CSs. METHODS: Independent component analysis was performed on datasets, detailing the transcriptional response of BECs to CS treatment. The expression of these CS-response components was examined in 2 patient cohorts and investigated in relation to clinical parameters. Supervised learning was used to predict BEC CS responses using peripheral blood gene expression. RESULTS: We identified a signature of CS response that was closely correlated with CS use in patients with asthma. Participants could be separated on the basis of CS-response genes into groups with high and low signature expression. Patients with low expression of CS-response genes, particularly those with a severe asthma diagnosis, showed worse lung function and quality of life. These individuals demonstrated enrichment for T-lymphocyte infiltration in endobronchial brushings. Supervised machine learning identified a 7-gene signature from peripheral blood that reliably identified patients with poor CS-response expression in BECs. CONCLUSIONS: Loss of CS transcriptional responses within bronchial epithelium was related to impaired lung function and poor quality of life, particularly in patients with severe asthma. These individuals were identified using minimally invasive blood sampling, suggesting these findings may enable earlier triage to alternative treatments.

CCL5 is a potential bridge between type 1 and type 2 inflammation in asthma.

BACKGROUND: Type 1 (T1) inflammation (marked by IFN-γ expression) is now consistently identified in subsets of asthma cohorts, but how it contributes to disease remains unclear. OBJECTIVE: We sought to understand the role of CCL5 in asthmatic T1 inflammation and how it interacts with both T1 and type 2 (T2) inflammation. METHODS: CCL5, CXCL9, and CXCL10 messenger RNA expression from sputum bulk RNA sequencing, as well as clinical and inflammatory data were obtained from the Severe Asthma Research Program III (SARP III). CCL5 and IFNG expression from bronchoalveolar lavage cell bulk RNA sequencing was obtained from the Immune Mechanisms in Severe Asthma (IMSA) cohort and expression related to previously identified immune cell profiles. The role of CCL5 in tissue-resident memory T-cell (TRM) reactivation was evaluated in a T1high murine severe asthma model. RESULTS: Sputum CCL5 expression strongly correlated with T1 chemokines (P < .001 for CXCL9 and CXCL10), consistent with a role in T1 inflammation. CCL5high participants had greater fractional exhaled nitric oxide (P = .009), blood eosinophils (P < .001), and sputum eosinophils (P = .001) in addition to sputum neutrophils (P = .001). Increased CCL5 bronchoalveolar lavage expression was unique to a previously described T1high/T2variable/lymphocytic patient group in the IMSA cohort, with IFNG trending with worsening lung obstruction only in this group (P = .083). In a murine model, high expression of the CCL5 receptor CCR5 was observed in TRMs and was consistent with a T1 signature. A role for CCL5 in TRM activation was supported by the ability of the CCR5 inhibitor maraviroc to blunt reactivation. CONCLUSION: CCL5 appears to contribute to TRM-related T1 neutrophilic inflammation in asthma while paradoxically also correlating with T2 inflammation and with sputum eosinophilia.

Pulmonary emphysema subtypes defined by unsupervised machine learning on CT scans.

BACKGROUND: Treatment and preventative advances for chronic obstructive pulmonary disease (COPD) have been slow due, in part, to limited subphenotypes. We tested if unsupervised machine learning on CT images would discover CT emphysema subtypes with distinct characteristics, prognoses and genetic associations. METHODS: New CT emphysema subtypes were identified by unsupervised machine learning on only the texture and location of emphysematous regions on CT scans from 2853 participants in the Subpopulations and Intermediate Outcome Measures in COPD Study (SPIROMICS), a COPD case-control study, followed by data reduction. Subtypes were compared with symptoms and physiology among 2949 participants in the population-based Multi-Ethnic Study of Atherosclerosis (MESA) Lung Study and with prognosis among 6658 MESA participants. Associations with genome-wide single-nucleotide-polymorphisms were examined. RESULTS: The algorithm discovered six reproducible (interlearner intraclass correlation coefficient, 0.91-1.00) CT emphysema subtypes. The most common subtype in SPIROMICS, the combined bronchitis-apical subtype, was associated with chronic bronchitis, accelerated lung function decline, hospitalisations, deaths, incident airflow limitation and a gene variant near DRD1, which is implicated in mucin hypersecretion (p=1.1 ×10-8). The second, the diffuse subtype was associated with lower weight, respiratory hospitalisations and deaths, and incident airflow limitation. The third was associated with age only. The fourth and fifth visually resembled combined pulmonary fibrosis emphysema and had distinct symptoms, physiology, prognosis and genetic associations. The sixth visually resembled vanishing lung syndrome. CONCLUSION: Large-scale unsupervised machine learning on CT scans defined six reproducible, familiar CT emphysema subtypes that suggest paths to specific diagnosis and personalised therapies in COPD and pre-COPD.

Genetic analyses of chr11p15.5 region identify MUC5AC-MUC5B associated with asthma-related phenotypes.

OBJECTIVE: Genome-wide association studies (GWASs) have identified single nucleotide polymorphisms (SNPs) in chr11p15.5 region associated with asthma and idiopathic interstitial pneumonias (IIPs). We sought to identify functional genes for asthma by combining SNPs and mRNA expression in bronchial epithelial cells (BEC) in the Severe Asthma Research Program (SARP). METHODS: Correlation analyses of mRNA expression of six candidate genes (AP2A2, MUC6, MUC2, MUC5AC, MUC5B, and TOLLIP) and asthma phenotypes were performed in the longitudinal cohort (n = 156) with RNAseq in BEC, and replicated in the cross-sectional cohort (n = 155). eQTL (n = 114) and genetic association analysis of asthma severity (426 severe vs. 531 non-severe asthma) were performed, and compared with previously published GWASs of IIPs and asthma. RESULTS: Higher expression of AP2A2 and MUC5AC and lower expression of MUC5B in BEC were correlated with asthma, asthma exacerbations, and T2 biomarkers (P < 0.01). SNPs associated with asthma and IIPs in previous GWASs were eQTL SNPs for MUC5AC, MUC5B, or TOLLIP, however, they were not in strong linkage disequilibrium. The risk alleles for asthma or protective alleles for IIPs were associated with higher expression of MUC5AC and lower expression of MUC5B. rs11603634, rs12788104, and rs28415845 associated with moderate-to-severe asthma or adult onset asthma in previous GWASs were not associated with asthma severity (P > 0.8). CONCLUSIONS: SNPs associated with asthma in chr11p15.5 region are not associated with asthma severity neither with IIPs. Higher expression of MUC5AC and lower expression of MUC5B are risk for asthma but protective for IIPs.

Investigations of a combination of atopic status and age of asthma onset identify asthma subphenotypes.

OBJECTIVE: Subphenotypes of asthma may be determined by age onset and atopic status. We sought to characterize early or late onset atopic asthma with fungal or non-fungal sensitization (AAFS or AANFS) and non-atopic asthma (NAA) in children and adults in the Severe Asthma Research Program (SARP). SARP is an ongoing project involving well-phenotyped patients with mild to severe asthma. METHODS: Phenotypic comparisons were performed using Kruskal-Wallis or chi-square test. Genetic association analyses were performed using logistic or linear regression. RESULTS: Airway hyper-responsiveness, total serum IgE levels, and T2 biomarkers showed an increasing trend from NAA to AANFS and then to AAFS. Children and adults with early onset asthma had greater % of AAFS than adults with late onset asthma (46% and 40% vs. 32%; P < 0.00001). In children, AAFS and AANFS had lower % predicted FEV1 (86% and 91% vs. 97%) and greater % of patients with severe asthma than NAA (61% and 59% vs. 43%). In adults with early or late onset asthma, NAA had greater % of patients with severe asthma than AANFS and AAFS (61% vs. 40% and 37% or 56% vs. 44% and 49%). The G allele of rs2872507 in GSDMB had higher frequency in AAFS than AANFS and NAA (0.63 vs. 0.55 and 0.55), and associated with earlier age onset and asthma severity. CONCLUSIONS: Early or late onset AAFS, AANFS, and NAA have shared and distinct phenotypic characteristics in children and adults. AAFS is a complex disorder involving genetic susceptibility and environmental factors.

Reversible Airflow Obstruction Predicts Future Chronic Obstructive Pulmonary Disease Development in the SPIROMICS Cohort: An Observational Cohort Study.

Rationale: Chronic obstructive pulmonary disease (COPD) is defined by fixed spirometric ratio, FEV1/FVC < 0.70 after inhaled bronchodilators. However, the implications of variable obstruction (VO), in which the prebronchodilator FEV1/FVC ratio is less than 0.70 but increases to 0.70 or more after inhaled bronchodilators, have not been determined. Objectives: We explored differences in physiology, exacerbations, and health status in participants with VO compared with reference participants without obstruction. Methods: Data from the SPIROMICS (Subpopulations and Intermediate Outcome Measures in COPD Study) cohort were obtained. Participants with VO were compared with reference participants without obstruction. Measurements and Main Results: We assessed differences in baseline radiographic emphysema and small airway disease at study entry, baseline, and change in lung function by spirometry, functional capacity by 6-minute walk, health status using standard questionnaires, exacerbation rates, and progression to COPD between the two groups. All models were adjusted for participant characteristics, asthma history, and tobacco exposure. We assessed 175 participants with VO and 603 reference participants without obstruction. Participants with VO had 6.2 times the hazard of future development of COPD controlling for other factors (95% confidence interval, 4.6-8.3; P < 0.001). Compared with reference participants, the VO group had significantly lower baseline pre- and post-bronchodilator (BD) FEV1, and greater decline over time in post-BD FEV1, and pre- and post-BD FVC. There were no significant differences in exacerbations between groups. Conclusions: Significant risk for future COPD development exists for those with pre- but not post-BD airflow obstruction. These findings support consideration of expanding spirometric criteria defining COPD to include pre-BD obstruction. Clinical trial registered with www.clinicaltrials.gov (NCT01969344).