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.

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.

Mitochondrial DNA Copy Number Variation in Asthma Risk, Severity, and Exacerbations

Rationale

Although airway oxidative stress and inflammation are central to asthma pathogenesis, there is limited knowledge of the relationship of asthma risk, severity, or exacerbations to mitochondrial dysfunction, which is pivotal to oxidant generation and inflammation.

Objectives

We investigated whether mitochondrial DNA copy number (mtDNA-CN) as a measure of mitochondrial function is associated with asthma diagnosis, severity, oxidative stress, and exacerbations.

Methods

We measured mtDNA-CN in blood in two cohorts. In the UK Biobank (UKB), we compared mtDNA-CN in mild and moderate-severe asthmatics to non-asthmatics. In the Severe Asthma Research Program (SARP), we evaluated mtDNA-CN in relation to asthma severity, biomarkers of oxidative stress and inflammation, and exacerbations.

Measures and Main Results

In UK Biobank, asthmatics (n = 29,768) have lower mtDNA-CN compared to non-asthmatics (n = 239,158) (beta, -0.026 [95% CI, -0.038 to -0.014], P = 2.46×10-5). While lower mtDNA-CN is associated with asthma, mtDNA-CN did not differ by asthma severity in either UKB or SARP. Biomarkers of inflammation show that asthmatics have higher white blood cells (WBC), neutrophils, eosinophils, fraction exhaled nitric oxide (FENO), and lower superoxide dismutase (SOD) than non-asthmatics, confirming greater oxidative stress in asthma. In one year follow-up in SARP, higher mtDNA-CN is associated with reduced risk of three or more exacerbations in the subsequent year (OR 0.352 [95% CI, 0.164 to 0.753], P = 0.007).

Conclusions

Asthma is characterized by mitochondrial dysfunction. Higher mtDNA-CN identifies an exacerbation-resistant asthma phenotype, suggesting mitochondrial function is important in exacerbation risk.

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.

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.

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.

Genetic control of mRNA splicing as a potential mechanism for incomplete penetrance of rare coding variants

Exonic variants present some of the strongest links between genotype and phenotype. However, these variants can have significant inter-individual pathogenicity differences, known as variable penetrance. In this study, we propose a model where genetically controlled mRNA splicing modulates the pathogenicity of exonic variants. By first cataloging exonic inclusion from RNA-sequencing data in GTEx V8, we find that pathogenic alleles are depleted on highly included exons. Using a large-scale phased whole genome sequencing data from the TOPMed consortium, we observe that this effect may be driven by common splice-regulatory genetic variants, and that natural selection acts on haplotype configurations that reduce the transcript inclusion of putatively pathogenic variants, especially when limiting to haploinsufficient genes. Finally, we test if this effect may be relevant for autism risk using families from the Simons Simplex Collection, but find that splicing of pathogenic alleles has a penetrance reducing effect here as well. Overall, our results indicate that common splice-regulatory variants may play a role in reducing the damaging effects of rare exonic variants.

Impact of Bronchiectasis on COPD Severity and Alpha-1 Antitrypsin Deficiency as a Risk Factor in Individuals with a Heavy Smoking History

Rationale

Bronchiectasis is common among those with heavy smoking histories, but risk factors for bronchiectasis, including alpha-1 antitrypsin deficiency, and its implications for COPD severity are uncharacterized in such individuals.

Objectives

To characterize the impact of bronchiectasis on COPD and explore alpha-1antitrypsin as a risk factor for bronchiectasis.

Methods

SubPopulations and InteRmediate Outcome Measures In COPD Study (SPIROMICS) participants (N=914; ages 40-80 years; ≥20-pack-year smoking) had high-resolution computed tomography (CT) scans interpreted visually for bronchiectasis, based on airway dilation without fibrosis or cicatrization. We performed regression-based models of bronchiectasis with clinical outcomes and quantitative CT measures. We deeply sequenced the gene encoding -alpha-1 antitrypsin, SERPINA1, in 835 participants to test for rare variants, focusing on the PiZ genotype (Glu366Lys, rs28929474).

Measurements and Main Results

We identified bronchiectasis in 365 (40%) participants, more frequently in women (45% versus 36%, p=0.0045), older participants (mean age=66[standard deviation (SD)=8.3] versus 64[SD=9.1] years, p=0.0083), and those with lower lung function (forced expiratory volume in 1 second [FEV1 ] percentage predicted=66%[SD=27] versus 77%[SD=25], p<0.0001; FEV1 to forced vital capacity [FVC] ratio=0.54[0.17] versus 0.63[SD=0.16], p<0.0001). Participants with bronchiectasis had greater emphysema (%voxels ≤-950 Hounsfield units, 11%[SD=12] versus 6.3%[SD=9], p<0.0001) and parametric response mapping functional small airways disease (26[SD=15] versus 19[SD=15], p<0.0001). Bronchiectasis was more frequent in the combined PiZZ and PiMZ genotype groups compared to those without PiZ, PiS, or other rare pathogenic variants (N=21 of 40 [52%] versus N=283 of 707[40%], odds ratio [OR]=1.97; 95% confidence interval [CI]=1.002, 3.90, p=0.049), an association attributed to White individuals (OR=1.98; 95%CI = 0.9956, 3.9; p=0.051).

Conclusions

Bronchiectasis was common in those with heavy smoking histories and was associated with detrimental clinical and radiographic outcomes. Our findings support alpha-1antitrypsin guideline recommendations to screen for alpha-1 antitrypsin deficiency in an appropriate bronchiectasis subgroup with a significant smoking history.

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.

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.

Clinical Implications of Longitudinal Blood Eosinophil Counts in Patients With Severe Asthma

BACKGROUND

The stability and variability of blood eosinophil counts (BECs) to phenotype patients with severe asthma is not fully understood.

OBJECTIVE

This post hoc, longitudinal, pooled analysis of placebo-arm patients from 2 phase 3 studies evaluated the clinical implications of BEC stability and variability in moderate-to-severe asthma.

METHODS

This analysis included patients from SIROCCO and CALIMA who received maintenance medium- to high-dosage inhaled corticosteroids plus long-acting β2-agonists; 2:1 patients with BECs of 300 cells/μL or higher and less than 300 cells/μL were enrolled. The BECs were measured 6 times over 1 year in a centralized laboratory. Exacerbations, lung function, and Asthma Control Questionnaire 6 scores were documented across patients grouped by BEC (<300 cells/μL or ≥300 cells/μL) and variability (<80% or ≥80% BECs less than or greater than 300 cells/μL).

RESULTS

Among 718 patients, 42.2% (n = 303) had predominantly high, 30.9% (n = 222) had predominantly low, and 26.9% (n = 193) had variable BECs. Prospective exacerbation rates (mean ± SD) were significantly greater in patients with predominantly high (1.39 ± 2.20) and variable (1.41 ± 2.09) BECs versus predominantly low (1.05 ± 1.66) BECs. Similar results were observed for the number of exacerbations while on placebo.

CONCLUSIONS

Although patients with variable BECs had intermittently high and low BECs, they experienced similar exacerbation rates to the predominantly high group, which were greater than those in the predominantly low group. A high BEC supports an eosinophilic phenotype in clinical settings without additional measurements, whereas a low BEC requires repeated measurements because it could reflect intermittently high or predominantly low BECs.

A blood and bronchoalveolar lavage protein signature of rapid FEV1 decline in smoking-associated COPD

Accelerated progression of chronic obstructive pulmonary disease (COPD) is associated with increased risks of hospitalization and death. Prognostic insights into mechanisms and markers of progression could facilitate development of disease-modifying therapies. Although individual biomarkers exhibit some predictive value, performance is modest and their univariate nature limits network-level insights. To overcome these limitations and gain insights into early pathways associated with rapid progression, we measured 1305 peripheral blood and 48 bronchoalveolar lavage proteins in individuals with COPD [n = 45, mean initial forced expiratory volume in one second (FEV1) 75.6 ± 17.4% predicted]. We applied a data-driven analysis pipeline, which enabled identification of protein signatures that predicted individuals at-risk for accelerated lung function decline (FEV1 decline ≥ 70 mL/year) ~ 6 years later, with high accuracy. Progression signatures suggested that early dysregulation in elements of the complement cascade is associated with accelerated decline. Our results propose potential biomarkers and early aberrant signaling mechanisms driving rapid progression in COPD.

Impact of Bronchiectasis on COPD Severity and Alpha-1 Antitrypsin Deficiency as a Risk Factor in Individuals with a Heavy Smoking History

Rationale

Bronchiectasis is common among those with heavy smoking histories, but risk factors for bronchiectasis, including α1-antitrypsin deficiency and its implications for COPD severity are uncharacterized in such individuals.

Objectives

To characterize the impact of bronchiectasis on COPD and explore α1-antitrypsin as a risk factor for bronchiectasis.

Methods

SPIROMICS participants (N=914; ages 40-80 years; ≥20 pack-years smoking) had HRCT scans interpreted visually for bronchiectasis, based on airway dilation without fibrosis or cicatrization. We performed regression-based models of bronchiectasis with clinical outcomes and quantitative CT measures. We deeply sequenced the gene encoding α1-antritrypsin, SERPINA1, in 835 participants to test for rare variants, focusing on PiZ (Glu³⁶⁶Lys, rs28929474).

Measurements and Main Results

We identified bronchiectasis in 365 (40%), more frequently in women (45% versus 36%, p=0.0045), older participants (mean age=66[SD=8.3] versus 64[SD=9.1] years, p=0.0083), and those with lower lung function (FEV₁%predicted=66%[SD=27] versus 77%[SD=25], p<0.0001; FEV₁/FVC=0.54[0.17] versus 0.63[SD=0.16], p<0.0001]. Participants with bronchiectasis had greater emphysema (%voxels ≤-950HFU, 11%[SD=12] versus 6.3%[SD=9], p<0.0001) and PRM^fSAD (26[SD=15] versus 19[SD=15], p<0.0001). Bronchiectasis was more frequent in the combined PiZZ and PiMZ genotype groups compared to those without PiZ, PiS, or other rare pathogenic variants (N=21 of 40[52%] versus N=283 of 707[40%], OR=1.97; 95%CI=1.002, 3.90, p=0.049), an association attributed to whites (OR=1.98; 95%CI = 0.9956, 3.9; p=0.051).

Conclusions

Bronchiectasis was common in those with heavy smoking histories and was associated with detrimental clinical and radiographic outcomes. Our findings support α1-antitrypsin guideline recommendations to screen for α1-antitrypsin deficiency in an appropriate bronchiectasis subgroup with a significant smoking history.

The genetic determinants of recurrent somatic mutations in 43,693 blood genomes

Nononcogenic somatic mutations are thought to be uncommon and inconsequential. To test this, we analyzed 43,693 National Heart, Lung and Blood Institute Trans-Omics for Precision Medicine blood whole genomes from 37 cohorts and identified 7131 non-missense somatic mutations that are recurrently mutated in at least 50 individuals. These recurrent non-missense somatic mutations (RNMSMs) are not clearly explained by other clonal phenomena such as clonal hematopoiesis. RNMSM prevalence increased with age, with an average 50-year-old having 27 RNMSMs. Inherited germline variation associated with RNMSM acquisition. These variants were found in genes involved in adaptive immune function, proinflammatory cytokine production, and lymphoid lineage commitment. In addition, the presence of eight specific RNMSMs associated with blood cell traits at effect sizes comparable to Mendelian genetic mutations. Overall, we found that somatic mutations in blood are an unexpectedly common phenomenon with ancestry-specific determinants and human health consequences.

Multi-ancestry genome-wide association analyses improve resolution of genes and pathways influencing lung function and chronic obstructive pulmonary disease risk

Lung-function impairment underlies chronic obstructive pulmonary disease (COPD) and predicts mortality. In the largest multi-ancestry genome-wide association meta-analysis of lung function to date, comprising 580,869 participants, we identified 1,020 independent association signals implicating 559 genes supported by ≥2 criteria from a systematic variant-to-gene mapping framework. These genes were enriched in 29 pathways. Individual variants showed heterogeneity across ancestries, age and smoking groups, and collectively as a genetic risk score showed strong association with COPD across ancestry groups. We undertook phenome-wide association studies for selected associated variants as well as trait and pathway-specific genetic risk scores to infer possible consequences of intervening in pathways underlying lung function. We highlight new putative causal variants, genes, proteins and pathways, including those targeted by existing drugs. These findings bring us closer to understanding the mechanisms underlying lung function and COPD, and should inform functional genomics experiments and potentially future COPD therapies.

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.

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.

Genetic control of mRNA splicing as a potential mechanism for incomplete penetrance of rare coding variants

Exonic variants present some of the strongest links between genotype and phenotype. However, these variants can have significant inter-individual pathogenicity differences, known as variable penetrance. In this study, we propose a model where genetically controlled mRNA splicing modulates the pathogenicity of exonic variants. By first cataloging exonic inclusion from RNA-seq data in GTEx v8, we find that pathogenic alleles are depleted on highly included exons. Using a large-scale phased WGS data from the TOPMed consortium, we observe that this effect may be driven by common splice-regulatory genetic variants, and that natural selection acts on haplotype configurations that reduce the transcript inclusion of putatively pathogenic variants, especially when limiting to haploinsufficient genes. Finally, we test if this effect may be relevant for autism risk using families from the Simons Simplex Collection, but find that splicing of pathogenic alleles has a penetrance reducing effect here as well. Overall, our results indicate that common splice-regulatory variants may play a role in reducing the damaging effects of rare exonic variants.

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.

The Impact of Insulin Resistance on Loss of Lung Function and Response to Treatment in Asthma

Rationale: The role of obesity-associated insulin resistance (IR) in airflow limitation in asthma is uncertain. Objectives: Using data in the Severe Asthma Research Program 3 (SARP-3), we evaluated relationships between homeostatic measure of IR (HOMA-IR), lung function (cross-sectional and longitudinal analyses), and treatment responses to bronchodilators and corticosteroids. Methods: HOMA-IR values were categorized as without (<3.0), moderate (3.0-5.0), or severe (>5.0). Lung function included FEV1 and FVC measured before and after treatment with inhaled albuterol and intramuscular triamcinolone acetonide and yearly for 5 years. Measurements and Main Results: Among 307 participants in SARP-3, 170 (55%) were obese and 140 (46%) had IR. Compared with patients without IR, those with IR had significantly lower values for FEV1 and FVC, and these lower values were not attributable to obesity effects. Compared with patients without IR, those with IR had lower FEV1 responses to β-adrenergic agonists and systemic corticosteroids. The annualized decline in FEV1 was significantly greater in patients with moderate IR (-41 ml/year) and severe IR (-32 ml/year,) than in patients without IR (-13 ml/year, P < 0.001 for both comparisons). Conclusions: IR is common in asthma and is associated with lower lung function, accelerated loss of lung function, and suboptimal lung function responses to bronchodilator and corticosteroid treatments. Clinical trials in patients with asthma and IR are needed to determine if improving IR might also improve lung function.

Rare genetic variants explain missing heritability in smoking

Common genetic variants explain less variation in complex phenotypes than inferred from family-based studies, and there is a debate on the source of this 'missing heritability'. We investigated the contribution of rare genetic variants to tobacco use with whole-genome sequences from up to 26,257 unrelated individuals of European ancestries and 11,743 individuals of African ancestries. Across four smoking traits, single-nucleotide-polymorphism-based heritability ([Formula: see text]) was estimated from 0.13 to 0.28 (s.e., 0.10-0.13) in European ancestries, with 35-74% of it attributable to rare variants with minor allele frequencies between 0.01% and 1%. These heritability estimates are 1.5-4 times higher than past estimates based on common variants alone and accounted for 60% to 100% of our pedigree-based estimates of narrow-sense heritability ([Formula: see text], 0.18-0.34). In the African ancestry samples, [Formula: see text] was estimated from 0.03 to 0.33 (s.e., 0.09-0.14) across the four smoking traits. These results suggest that rare variants are important contributors to the heritability of smoking.

Clinical and molecular implications of RGS2 promoter genetic variation in severe asthma

BACKGROUND

Regulator of G protein signaling (RGS) 2 terminates bronchoconstrictive Gαq signaling; murine RGS2 knockout demonstrate airway hyperresponsiveness. While RGS2 promoter variants rs2746071 and rs2746072 associate with a clinical mild asthma phenotype, their impact on human airway smooth muscle (HASM) contractility and asthma severity outcomes is unknown.

OBJECTIVE

We sought to determine whether reductions in RGS2 expression seen with these 2 RGS2 promoter variants augment HASM contractility and associate with an asthma severity phenotype.

METHODS

We transfected HASM with a range of RGS2-specific small interfering RNA (siRNA) concentrations and determined RGS2 protein expression by Western blot analysis and intracellular calcium flux induced by histamine (a Gαq-coupled H1 receptor bronchoconstrictive agonist). We conducted regression-based genotype association analyses of RGS2 variants from 611 patients from the National Heart, Lung, and Blood Institute Severe Asthma Research Program 3.

RESULTS

RGS2-specific siRNA caused dose-dependent increases in histamine-stimulated bronchoconstrictive intracellular calcium signaling (2-way ANOVA, P < .0001) with a concomitant decrease in RGS2 protein expression. RGS2-specific siRNA did not affect Gαq-independent ionomycin-induced intracellular calcium signaling (P = .42). The minor allele frequency of rs2746071 and rs2746072 was 0.46 and 0.28 among African American/non-Hispanic Black patients and was 0.28 and 0.27 among non-Hispanic White patients, among whom these single nucleotide polymorphisms were in stronger linkage disequilibrium (r2 = 0.97). Among non-Hispanic White patients, risk allele homozygotes for rs2746072 and rs2746071 each had nearly 2-fold greater asthma exacerbation rates relative to alternative genotypes with wild-type alleles (Padditive = 2.86 × 10-5/Precessive = 5.22 × 10-6 and Padditive = 3.46 × 10-6/Precessive = 6.74 × 10-7, respectively) at baseline, which was confirmed by prospective longitudinal exacerbation data.

CONCLUSION

RGS2 promoter variation associates with a molecular and clinical phenotype characterized by enhanced bronchoconstrictive stimulation in vitro and higher asthma exacerbations rates in non-Hispanic White patients.

DNA sequencing analysis of cystic fibrosis transmembrane conductance regulator gene identifies cystic fibrosis-associated variants in the Severe Asthma Research Program

BACKGROUND

Heterozygote carriers of potentially pathogenic variants in the cystic fibrosis transmembrane conductance regulator (CFTR) gene have increased asthma risk. However, the frequency and impact of CFTR variation among individuals with asthma is unknown.

OBJECTIVE

To determine whether potentially pathogenic CFTR variants associate with disease severity and whether individuals with two potentially pathogenic variants exist in a severe asthma-enriched cohort.

METHODS

We analyzed sequencing data spanning a 190.5Kb region of CFTR in participants from the Severe Asthma Research Program (SARP1-3). Potentially pathogenic, rare CFTR variants (frequency < 0.05) were classified as CF-causing or of varying clinical consequences (VVCC) (CFTR2. org). Regression-based models tested for association between CFTR genotypes (0-2 potentially pathogenic variants) and severity outcomes.

RESULTS

Of 1401 participants, 9.5% (134) had one potentially pathogenic variant, occurring more frequently in non-Hispanic white (NHW, 10.1% [84 of 831]) compared to African American individuals (AA, 5.2% [22 of 426]). We found ≥2 potentially pathogenic CFTR variants in 1.4% (19); 0.5% (4) of NHW and 2.8% (12) of AA. Potentially pathogenic CFTR variant genotypes (≥1 or ≥2 variants) were not cumulatively associated with lung function or exacerbations. In NHW, we found three F508del compound heterozygotes with F508del and a VVCC (two 5 T; TG12[c.1210-11 T > G] and one Arg1070Trp) and a homozygote for the VVCC, 5 T; TG12.

CONCLUSIONS

We found potentially pathogenic CFTR variants within a severe asthma-enriched cohort, including three compound heterozygote genotypes variably associated with CF in NHW individuals. These findings provide the rationale for CFTR sequencing and phenotyping of CF-related traits in individuals with severe asthma.

Lung function, airway and peripheral basophils and eosinophils are associated with molecular pharmacogenomic endotypes of steroid response in severe asthma

INTRODUCTION

Asthma is a complex disease with heterogeneous expression/severity. There is growing interest in defining asthma endotypes consistently associated with different responses to therapy, focusing on type 2 inflammation (Th2) as a key pathological mechanism. Current asthma endotypes are defined primarily by clinical/laboratory criteria. Each endotype is likely characterised by distinct molecular mechanisms that identify optimal therapies.

METHODS

We applied unsupervised (without a priori clinical criteria) principal component analysis on sputum airway cells RNA-sequencing transcriptomic data from 19 asthmatics from the Severe Asthma Research Program at baseline and 6-8 weeks follow-up after a 40 mg dose of intramuscular corticosteroids. We investigated principal components PC1, PC3 for association with 55 clinical variables.

RESULTS

PC3 was associated with baseline Th2 clinical features including blood (rank-sum p=0.0082) and airway (rank-sum p=0.0024) eosinophilia, FEV₁ change (Kendall tau-b R=-0.333 (-0.592 to -0.012)) and follow-up FEV₁ albuterol response (Kendall tau-b R=0.392 (0.079 to 0.634)). PC1 with blood basophlia (rank-sum p=0.0191). The top 5% genes contributing to PC1, PC3 were enriched for distinct immune system/inflammation ontologies suggesting distinct subject-specific clusters of transcriptomic response to corticosteroids. PC3 association with FEV₁ change was reproduced in silico in a comparable independent 14-subject (baseline, 8 weeks after daily inhaled corticosteroids (ICS)) airway epithelial cells microRNAome dataset.

CONCLUSIONS

Transcriptomic PCs from this unsupervised methodology define molecular pharmacogenomic endotypes that may yield novel biology underlying different subject-specific responses to corticosteroid therapy in asthma, and optimal personalised asthma care. Top contributing genes to these PCs may suggest new therapeutic targets.

Genetic Associations and Architecture of Asthma-COPD Overlap

BACKGROUND

Some people have characteristics of both asthma and COPD (asthma-COPD overlap), and evidence suggests they experience worse outcomes than those with either condition alone.

RESEARCH QUESTION

What is the genetic architecture of asthma-COPD overlap, and do the determinants of risk for asthma-COPD overlap differ from those for COPD or asthma?

STUDY DESIGN AND METHODS

We conducted a genome-wide association study in 8,068 asthma-COPD overlap case subjects and 40,360 control subjects without asthma or COPD of European ancestry in UK Biobank (stage 1). We followed up promising signals (P < 5 × 10^-6) that remained associated in analyses comparing (1) asthma-COPD overlap vs asthma-only control subjects, and (2) asthma-COPD overlap vs COPD-only control subjects. These variants were analyzed in 12 independent cohorts (stage 2).

RESULTS

We selected 31 independent variants for further investigation in stage 2, and discovered eight novel signals (P < 5 × 10^-8) for asthma-COPD overlap (meta-analysis of stage 1 and 2 studies). These signals suggest a spectrum of shared genetic influences, some predominantly influencing asthma (FAM105A, GLB1, PHB, TSLP), others predominantly influencing fixed airflow obstruction (IL17RD, C5orf56, HLA-DQB1). One intergenic signal on chromosome 5 had not been previously associated with asthma, COPD, or lung function. Subgroup analyses suggested that associations at these eight signals were not driven by smoking or age at asthma diagnosis, and in phenome-wide scans, eosinophil counts, atopy, and asthma traits were prominent.

INTERPRETATION

We identified eight signals for asthma-COPD overlap, which may represent loci that predispose to type 2 inflammation, and serious long-term consequences of asthma.

Quantitative CT Characteristics of Cluster Phenotypes in the Severe Asthma Research Program Cohorts

Background Clustering key clinical characteristics of participants in the Severe Asthma Research Program (SARP), a large, multicenter prospective observational study of patients with asthma and healthy controls, has led to the identification of novel asthma phenotypes. Purpose To determine whether quantitative CT (qCT) could help distinguish between clinical asthma phenotypes. Materials and Methods A retrospective cross-sectional analysis was conducted with the use of qCT images (maximal bronchodilation at total lung capacity [TLC], or inspiration, and functional residual capacity [FRC], or expiration) from the cluster phenotypes of SARP participants (cluster 1: minimal disease; cluster 2: mild, reversible; cluster 3: obese asthma; cluster 4: severe, reversible; cluster 5: severe, irreversible) enrolled between September 2001 and December 2015. Airway morphometry was performed along standard paths (RB1, RB4, RB10, LB1, and LB10). Corresponding voxels from TLC and FRC images were mapped with use of deformable image registration to characterize disease probability maps (DPMs) of functional small airway disease (fSAD), voxel-level volume changes (Jacobian), and isotropy (anisotropic deformation index [ADI]). The association between cluster assignment and qCT measures was evaluated using linear mixed models. Results A total of 455 participants were evaluated with cluster assignments and CT (mean age ± SD, 42.1 years ± 14.7; 270 women). Airway morphometry had limited ability to help discern between clusters. DPM fSAD was highest in cluster 5 (cluster 1 in SARP III: 19.0% ± 20.6; cluster 2: 18.9% ± 13.3; cluster 3: 24.9% ± 13.1; cluster 4: 24.1% ± 8.4; cluster 5: 38.8% ± 14.4; P < .001). Lower whole-lung Jacobian and ADI values were associated with greater cluster severity. Compared to cluster 1, cluster 5 lung expansion was 31% smaller (Jacobian in SARP III cohort: 2.31 ± 0.6 vs 1.61 ± 0.3, respectively, P < .001) and 34% more isotropic (ADI in SARP III cohort: 0.40 ± 0.1 vs 0.61 ± 0.2, P < .001). Within-lung Jacobian and ADI SDs decreased as severity worsened (Jacobian SD in SARP III cohort: 0.90 ± 0.4 for cluster 1; 0.79 ± 0.3 for cluster 2; 0.62 ± 0.2 for cluster 3; 0.63 ± 0.2 for cluster 4; and 0.41 ± 0.2 for cluster 5; P < .001). Conclusion Quantitative CT assessments of the degree and intraindividual regional variability of lung expansion distinguished between well-established clinical phenotypes among participants with asthma from the Severe Asthma Research Program study. © RSNA, 2022 Online supplemental material is available for this article. See also the editorial by Verschakelen in this issue.

Mendelian randomization supports bidirectional causality between telomere length and clonal hematopoiesis of indeterminate potential

Human genetic studies support an inverse causal relationship between leukocyte telomere length (LTL) and coronary artery disease (CAD), but directionally mixed effects for LTL and diverse malignancies. Clonal hematopoiesis of indeterminate potential (CHIP), characterized by expansion of hematopoietic cells bearing leukemogenic mutations, predisposes both hematologic malignancy and CAD. TERT (which encodes telomerase reverse transcriptase) is the most significantly associated germline locus for CHIP in genome-wide association studies. Here, we investigated the relationship between CHIP, LTL, and CAD in the Trans-Omics for Precision Medicine (TOPMed) program (n = 63,302) and UK Biobank (n = 47,080). Bidirectional Mendelian randomization studies were consistent with longer genetically imputed LTL increasing propensity to develop CHIP, but CHIP then, in turn, hastens to shorten measured LTL (mLTL). We also demonstrated evidence of modest mediation between CHIP and CAD by mLTL. Our data promote an understanding of potential causal relationships across CHIP and LTL toward prevention of CAD.

A genome-wide association study of bronchodilator response in participants of European and African ancestry from six independent cohorts

Introduction

Bronchodilator response (BDR) is a measurement of acute bronchodilation in response to short-acting β2-agonists, with a heritability between 10 and 40%. Identifying genetic variants associated with BDR may lead to a better understanding of its complex pathophysiology.

Methods

We performed a genome-wide association study (GWAS) of BDR in six adult cohorts with participants of European ancestry (EA) and African ancestry (AA) including community cohorts and cohorts ascertained on the basis of obstructive pulmonary disease. Validation analysis was carried out in two paediatric asthma cohorts.

Results

A total of 10 623 EA and 3597 AA participants were included in the analyses. No single nucleotide polymorphism (SNP) was associated with BDR at the conventional genome-wide significance threshold (p<5×10-8). Performing fine mapping and using a threshold of p<5×10-6 to identify suggestive variants of interest, we identified three SNPs with possible biological relevance: rs35870000 (within FREM1), which may be involved in IgE- and IL5-induced changes in airway smooth muscle cell responsiveness; rs10426116 (within ZNF284), a zinc finger protein, which has been implicated in asthma and BDR previously; and rs4782614 (near ATP2C2), involved in calcium transmembrane transport. Validation in paediatric cohorts yielded no significant SNPs, possibly due to age-genotype interaction effects.

Conclusion

Ancestry-stratified and ancestry-combined GWAS meta-analyses of over 14 000 participants did not identify genetic variants associated with BDR at the genome-wide significance threshold, although a less stringent threshold identified three variants showing suggestive evidence of association. A common definition and protocol for measuring BDR in research may improve future efforts to identify variants associated with BDR.

Genetic determinants of telomere length from 109,122 ancestrally diverse whole-genome sequences in TOPMed

Genetic studies on telomere length are important for understanding age-related diseases. Prior GWAS for leukocyte TL have been limited to European and Asian populations. Here, we report the first sequencing-based association study for TL across ancestrally-diverse individuals (European, African, Asian and Hispanic/Latino) from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program. We used whole genome sequencing (WGS) of whole blood for variant genotype calling and the bioinformatic estimation of telomere length in n=109,122 individuals. We identified 59 sentinel variants (p-value <5×10-9) in 36 loci associated with telomere length, including 20 newly associated loci (13 were replicated in external datasets). There was little evidence of effect size heterogeneity across populations. Fine-mapping at OBFC1 indicated the independent signals colocalized with cell-type specific eQTLs for OBFC1 (STN1). Using a multi-variant gene-based approach, we identified two genes newly implicated in telomere length, DCLRE1B (SNM1B) and PARN. In PheWAS, we demonstrated our TL polygenic trait scores (PTS) were associated with increased risk of cancer-related phenotypes.

Genetic variation in genes regulating skeletal muscle regeneration and tissue remodelling associated with weight loss in chronic obstructive pulmonary disease

BACKGROUND

Chronic obstructive pulmonary disease (COPD) is the third leading cause of death globally. COPD patients with cachexia or weight loss have increased risk of death independent of body mass index (BMI) and lung function. We tested the hypothesis genetic variation is associated with weight loss in COPD using a genome-wide association study approach.

METHODS

Participants with COPD (N = 4308) from three studies (COPDGene, ECLIPSE, and SPIROMICS) were analysed. Discovery analyses were performed in COPDGene with replication in SPIROMICS and ECLIPSE. In COPDGene, weight loss was defined as self-reported unintentional weight loss > 5% in the past year or low BMI (BMI < 20 kg/m2 ). In ECLIPSE and SPIROMICS, weight loss was calculated using available longitudinal visits. Stratified analyses were performed among African American (AA) and Non-Hispanic White (NHW) participants with COPD. Single variant and gene-based analyses were performed adjusting for confounders. Fine mapping was performed using a Bayesian approach integrating genetic association results with linkage disequilibrium and functional annotation. Significant gene networks were identified by integrating genetic regions associated with weight loss with skeletal muscle protein-protein interaction (PPI) data.

RESULTS

At the single variant level, only the rs35368512 variant, intergenic to GRXCR1 and LINC02383, was associated with weight loss (odds ratio = 3.6, 95% confidence interval = 2.3-5.6, P = 3.2 × 10-8 ) among AA COPD participants in COPDGene. At the gene level in COPDGene, EFNA2 and BAIAP2 were significantly associated with weight loss in AA and NHW COPD participants, respectively. The EFNA2 association replicated among AA from SPIROMICS (P = 0.0014), whereas the BAIAP2 association replicated in NHW from ECLIPSE (P = 0.025). The EFNA2 gene encodes the membrane-bound protein ephrin-A2 involved in the regulation of developmental processes and adult tissue homeostasis such as skeletal muscle. The BAIAP2 gene encodes the insulin-responsive protein of mass 53 kD (IRSp53), a negative regulator of myogenic differentiation. Integration of the gene-based findings participants with PPI data revealed networks of genes involved in pathways such as Rho and synapse signalling.

CONCLUSIONS

The EFNA2 and BAIAP2 genes were significantly associated with weight loss in COPD participants. Collectively, the integrative network analyses indicated genetic variation associated with weight loss in COPD may influence skeletal muscle regeneration and tissue remodelling.

Genome-wide association study of asthma, total IgE, and lung function in a cohort of Peruvian children

BACKGROUND

Genetic ancestry plays a role in asthma health disparities.

OBJECTIVE

Our aim was to evaluate the impact of ancestry on and identify genetic variants associated with asthma, total serum IgE level, and lung function.

METHODS

A total of 436 Peruvian children (aged 9-19 years) with asthma and 291 without asthma were genotyped by using the Illumina Multi-Ethnic Global Array. Genome-wide proportions of indigenous ancestry populations from continental America (NAT) and European ancestry from the Iberian populations in Spain (IBS) were estimated by using ADMIXTURE. We assessed the relationship between ancestry and the phenotypes and performed a genome-wide association study.

RESULTS

The mean ancestry proportions were 84.7% NAT (case patients, 84.2%; controls, 85.4%) and 15.3% IBS (15.8%; 14.6%). With adjustment for asthma, NAT was associated with higher total serum IgE levels (P < .001) and IBS was associated with lower total serum IgE levels (P < .001). NAT was associated with higher FEV1 percent predicted values (P < .001), whereas IBS was associated with lower FEV1 values in the controls but not in the case patients. The HLA-DR/DQ region on chromosome 6 (Chr6) was strongly associated with total serum IgE (rs3135348; P = 3.438 × 10-10) and was independent of an association with the haplotype HLA-DQA1∼HLA-DQB1:04.01∼04.02 (P = 1.55 × 10-05). For lung function, we identified a locus (rs4410198; P = 5.536 × 10-11) mapping to Chr19, near a cluster of zinc finger interacting genes that colocalizes to the long noncoding RNA CTD-2537I9.5. This novel locus was replicated in an independent sample of pediatric case patients with asthma with similar admixture from Brazil (P = .005).

CONCLUSION

This study confirms the role of HLA in atopy, and identifies a novel locus mapping to a long noncoding RNA for lung function that may be specific to children with NAT.

Pharmacogenetic studies of long-acting beta agonist and inhaled corticosteroid responsiveness in randomised controlled trials of individuals of African descent with asthma

BACKGROUND

Pharmacogenetic studies in asthma cohorts, primarily made up of White people of European descent, have identified loci associated with response to inhaled beta agonists and corticosteroids (ICSs). Differences exist in how individuals from different ancestral backgrounds respond to long-acting beta agonist (LABA) and ICSs. Therefore, we sought to understand the pharmacogenetic mechanisms regulating therapeutic responsiveness in individuals of African descent.

METHODS

We did ancestry-based pharmacogenetic studies of children (aged 5-11 years) and adolescents and adults (aged 12-69 years) from the Best African Response to Drug (BARD) trials, in which participants with asthma uncontrolled with low-dose ICS (fluticasone propionate 50 μg in children, 100 μg in adolescents and adults) received different step-up combination therapies. The hierarchal composite outcome of pairwise superior responsiveness in BARD was based on asthma exacerbations, a 31-day difference in annualised asthma-control days, or a 5% difference in percentage predicted FEV₁. We did whole-genome admixture mapping of 15 159 ancestral segments within 312 independent regions, stratified by the two age groups. The two co-primary outcome comparisons were the step up from low-dose ICS to the quintuple dose of ICS (5 × ICS: 250 μg twice daily in children and 500 μg twice daily in adolescents and adults) versus double dose (2-2·5 × ICS: 100 μg twice daily in children, 250 μg twice daily in adolescents and adults), and 5 × ICS versus 100 μg fluticasone plus a LABA (salmeterol 50 μg twice daily). We used a genome-wide significance threshold of p<1·6 × 10^-4, and tested for replication using independent cohorts of individuals of African descent with asthma.

FINDINGS

We included 249 unrelated children and 267 unrelated adolescents and adults in the BARD pharmacogenetic analysis. In children, we identified a significant admixture mapping peak for superior responsiveness to 5 × ICS versus 100 μg fluticasone plus salmeterol on chromosome 12 (odds ratio [OR_{local African}] 3·95, 95% CI 2·02-7·72, p=6·1 × 10^-5) fine mapped to a locus adjacent to RNFT2 and NOS1 (rs73399224, OR_{allele dose} 0·17, 95% CI 0·07-0·42, p=8·4 × 10^-5). In adolescents and adults, we identified a peak for superior responsiveness to 5 × ICS versus 2·5 × ICS on chromosome 22 (OR_{local African} 3·35, 1·98-5·67, p=6·8 × 10^-6) containing a locus adjacent to TPST2 (rs5752429, OR_{allele dose} 0·21, 0·09-0·52, p=5·7 × 10^-4). We replicated rs5752429 and nominally replicated rs73399224 in independent African American cohorts.

INTERPRETATION

BARD is the first genome-wide pharmacogenetic study of LABA and ICS response in clinical trials of individuals of African descent to detect and replicate genome-wide significant loci. Admixture mapping of the composite BARD trial outcome enabled the identification of novel pharmacogenetic variation accounting for differential therapeutic responses in people of African descent with asthma.

FUNDING

National Institutes of Health, National Heart, Lung, and Blood Institute.

Whole-genome sequencing in diverse subjects identifies genetic correlates of leukocyte traits: The NHLBI TOPMed program

Many common and rare variants associated with hematologic traits have been discovered through imputation on large-scale reference panels. However, the majority of genome-wide association studies (GWASs) have been conducted in Europeans, and determining causal variants has proved challenging. We performed a GWAS of total leukocyte, neutrophil, lymphocyte, monocyte, eosinophil, and basophil counts generated from 109,563,748 variants in the autosomes and the X chromosome in the Trans-Omics for Precision Medicine (TOPMed) program, which included data from 61,802 individuals of diverse ancestry. We discovered and replicated 7 leukocyte trait associations, including (1) the association between a chromosome X, pseudo-autosomal region (PAR), noncoding variant located between cytokine receptor genes (CSF2RA and CLRF2) and lower eosinophil count; and (2) associations between single variants found predominantly among African Americans at the S1PR3 (9q22.1) and HBB (11p15.4) loci and monocyte and lymphocyte counts, respectively. We further provide evidence indicating that the newly discovered eosinophil-lowering chromosome X PAR variant might be associated with reduced susceptibility to common allergic diseases such as atopic dermatitis and asthma. Additionally, we found a burden of very rare FLT3 (13q12.2) variants associated with monocyte counts. Together, these results emphasize the utility of whole-genome sequencing in diverse samples in identifying associations missed by European-ancestry-driven GWASs.

Benefits of Airway Androgen Receptor Expression in Human Asthma

Rationale: Androgens are potentially beneficial in asthma, but AR (androgen receptor) has not been studied in human airways.Objectives: To measure whether AR and its ligands are associated with human asthma outcomes.Methods: We compared the effects of AR expression on lung function, symptom scores, and fractional exhaled nitric oxide (Fe_NO) in adults enrolled in SARP (Severe Asthma Research Program). The impact of sex and of androgens on asthma outcomes was also evaluated in the SARP with validation studies in the Cleveland Clinic Health System and the NHANES (U.S. National Health and Nutrition Examination Survey).Measurements and Main Results: In SARP (n = 128), AR gene expression from bronchoscopic epithelial brushings was positively associated with both FEV₁/FVC ratio (R² = 0.135, P = 0.0002) and the total Asthma Quality of Life Questionnaire score (R² = 0.056, P = 0.016) and was negatively associated with Fe_NO (R² = 0.178, P = 9.8 × 10^-6) and NOS2 (nitric oxide synthase gene) expression (R² = 0.281, P = 1.2 × 10^-10). In SARP (n = 1,659), the Cleveland Clinic Health System (n = 32,527), and the NHANES (n = 2,629), women had more asthma exacerbations and emergency department visits than men. The levels of the AR ligand precursor dehydroepiandrosterone sulfate correlated positively with the FEV₁ in both women and men.Conclusions: Higher bronchial AR expression and higher androgen levels are associated with better lung function, fewer symptoms, and a lower Fe_NO in human asthma. The role of androgens should be considered in asthma management.

Mixed Sputum Granulocyte Longitudinal Impact on Lung Function in the Severe Asthma Research Program

Rationale: Some reports indicate longitudinal variability in sputum differential cell counts, whereas others describe stability. Highly variable sputum eosinophil percentages are associated with greater lung function loss than persistently elevated eosinophil percentages, but elevated neutrophils are linked to more severe asthma.Objectives: To examine sputum granulocyte stability or variability longitudinally and associations with important clinical characteristics.Methods: The SARP III (Severe Asthma Research Program III) cohort underwent comprehensive phenotype characterization at baseline and annually over 3 years. Adult subjects with acceptable sputum levels were assigned to one of three longitudinal sputum groups: eosinophils predominantly <2%, eosinophils predominantly ≥2%, or highly variable eosinophil percentages (>2 SDs determined from independent, repeated baseline eosinophil percentages). Subjects were similarly assigned to one of three longitudinal neutrophil groups with a 50% cut point.Measurements and Main Results: The group with predominantly <2% sputum eosinophils had the highest lung function (prebronchodilator FEV₁% predicted, P < 0.01; FEV₁/FVC ratio, P < 0.001) at baseline and throughout 3 years compared with other eosinophil groups. Healthcare use did not differ, although the highly variable eosinophil group reported more asthma exacerbations at Year 3. Longitudinal neutrophil groups showed few differences. However, a combination of predominantly ≥2% eosinophil and ≥50% neutrophil groups resulted in the lowest prebronchodilator FEV₁% predicted (P = 0.049) compared with the combination with predominantly <2% eosinophils and<50% neutrophils.Conclusions: Subjects with predominantly ≥2% sputum eosinophils in combination with predominantly ≥50% neutrophils showed greater loss of lung function, whereas those with highly variable sputum eosinophils had greater healthcare use.

Responsiveness to Parenteral Corticosteroids and Lung Function Trajectory in Adults with Moderate-to-Severe Asthma

Rationale: It is unclear why select patients with moderate-to-severe asthma continue to lose lung function despite therapy. We hypothesized that participants with the smallest responses to parenteral corticosteroids have the greatest risk of undergoing a severe decline in lung function.Objectives: To evaluate corticosteroid-response phenotypes as longitudinal predictors of lung decline.Methods: Adults within the NHLBI SARP III (Severe Asthma Research Program III) who had undergone a course of intramuscular triamcinolone at baseline and at ≥2 annual follow-up visits were evaluated. Longitudinal slopes were calculated for each participant's post-bronchodilator FEV₁% predicted. Categories of participant FEV₁ slope were defined: severe decline, >2% loss/yr; mild decline, >0.5-2.0% loss/yr; no change, 0.5% loss/yr to <1% gain/yr; and improvement, ≥1% gain/yr. Regression models were used to develop predictors of severe decline.Measurements and Main Results: Of 396 participants, 78 had severe decline, 91 had mild decline, 114 had no change, and 113 showed improvement. The triamcinolone-induced difference in the post-bronchodilator FEV₁% predicted (derived by baseline subtraction) was related to the 4-year change in lung function or slope category in univariable models (P < 0.001). For each 5% decrement in the triamcinolone-induced difference the FEV1% predicted, there was a 50% increase in the odds of being in the severe decline group (odds ratio, 1.5; 95% confidence interval, 1.3-1.8), when adjusted for baseline FEV₁, exacerbation history, blood eosinophils and body mass index.Conclusions: Failure to improve the post-bronchodilator FEV₁ after a challenge with parenteral corticosteroids is an evoked biomarker for patients at risk for a severe decline in lung function.

Genetic and non-genetic factors affecting the expression of COVID-19-relevant genes in the large airway epithelium

BACKGROUND

The large airway epithelial barrier provides one of the first lines of defense against respiratory viruses, including SARS-CoV-2 that causes COVID-19. Substantial inter-individual variability in individual disease courses is hypothesized to be partially mediated by the differential regulation of the genes that interact with the SARS-CoV-2 virus or are involved in the subsequent host response. Here, we comprehensively investigated non-genetic and genetic factors influencing COVID-19-relevant bronchial epithelial gene expression.

METHODS

We analyzed RNA-sequencing data from bronchial epithelial brushings obtained from uninfected individuals. We related ACE2 gene expression to host and environmental factors in the SPIROMICS cohort of smokers with and without chronic obstructive pulmonary disease (COPD) and replicated these associations in two asthma cohorts, SARP and MAST. To identify airway biology beyond ACE2 binding that may contribute to increased susceptibility, we used gene set enrichment analyses to determine if gene expression changes indicative of a suppressed airway immune response observed early in SARS-CoV-2 infection are also observed in association with host factors. To identify host genetic variants affecting COVID-19 susceptibility in SPIROMICS, we performed expression quantitative trait (eQTL) mapping and investigated the phenotypic associations of the eQTL variants.

RESULTS

We found that ACE2 expression was higher in relation to active smoking, obesity, and hypertension that are known risk factors of COVID-19 severity, while an association with interferon-related inflammation was driven by the truncated, non-binding ACE2 isoform. We discovered that expression patterns of a suppressed airway immune response to early SARS-CoV-2 infection, compared to other viruses, are similar to patterns associated with obesity, hypertension, and cardiovascular disease, which may thus contribute to a COVID-19-susceptible airway environment. eQTL mapping identified regulatory variants for genes implicated in COVID-19, some of which had pheWAS evidence for their potential role in respiratory infections.

CONCLUSIONS

These data provide evidence that clinically relevant variation in the expression of COVID-19-related genes is associated with host factors, environmental exposures, and likely host genetic variation.

Sequencing of 53,831 diverse genomes from the NHLBI TOPMed Program

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes)1. In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%.

Comparison of whole genome sequencing and targeted sequencing for mitochondrial DNA

Mitochondrial dysfunction has emerged to be associated with a broad spectrum of diseases, and there is an increasing demand for accurate detection of mitochondrial DNA (mtDNA) variants. Whole genome sequencing (WGS) has been the dominant sequencing approach to identify genetic variants in recent decades, but most studies focus on variants on the nuclear genome. Whole genome sequencing is also costly and time consuming. Sequencing specifically targeted for mtDNA is commonly used in the diagnostic settings and has lower costs. However, there is a lack of pairwise comparisons between these two sequencing approaches for calling mtDNA variants on a population basis. In this study, we compared WGS and mtDNA-targeted sequencing (targeted-seq) in analyzing mitochondrial DNA from 1499 participants recruited into the Severe Asthma Research Program (SARP). Our study reveals that targeted-sequencing and WGS have comparable capacity to determine genotypes and to call haplogroups and homoplasmies on mtDNA. However, there exists a large variability in calling heteroplasmies, especially for low-frequency heteroplasmies, which indicates that investigators should be cautious about heteroplasmies acquired from different sequencing methods. Further research is highly desired to improve variant detection methods for mitochondrial DNA.

Loss-of-function genomic variants highlight potential therapeutic targets for cardiovascular disease

Pharmaceutical drugs targeting dyslipidemia and cardiovascular disease (CVD) may increase the risk of fatty liver disease and other metabolic disorders. To identify potential novel CVD drug targets without these adverse effects, we perform genome-wide analyses of participants in the HUNT Study in Norway (n = 69,479) to search for protein-altering variants with beneficial impact on quantitative blood traits related to cardiovascular disease, but without detrimental impact on liver function. We identify 76 (11 previously unreported) presumed causal protein-altering variants associated with one or more CVD- or liver-related blood traits. Nine of the variants are predicted to result in loss-of-function of the protein. This includes ZNF529:p.K405X, which is associated with decreased low-density-lipoprotein (LDL) cholesterol (P = 1.3 × 10-8) without being associated with liver enzymes or non-fasting blood glucose. Silencing of ZNF529 in human hepatoma cells results in upregulation of LDL receptor and increased LDL uptake in the cells. This suggests that inhibition of ZNF529 or its gene product should be prioritized as a novel candidate drug target for treating dyslipidemia and associated CVD.

The chemokine CX3CL1/fractalkine regulates immunopathogenesis during fungal-associated allergic airway inflammation

Individuals that present with difficult-to-control asthma and sensitivity to one or more fungal species are categorized as a subset of severe asthma patients belonging to a group herein referred to as severe asthma with fungal sensitization (SAFS). We have previously reported the identification of numerous cytokines and chemokines that were elevated in human asthmatics that were sensitized to fungi vs. nonfungal sensitized asthmatics. Here, we show that the unique chemokine CX3CL1 (fractalkine) is elevated in both bronchoalveolar lavage fluid and sputum from human asthmatics sensitized to fungi, implicating an association with CX3CL1 in fungal asthma severity. In an experimental model of fungal-associated allergic airway inflammation, we demonstrate that the absence of CX3CR1 signaling unexpectedly resulted in a profound impairment in lung function. Histological assessment of lung tissue revealed an unrestricted inflammatory response that was subsequently characterized by enhanced levels of neutrophils, eosinophils, and inflammatory monocytes. Neutrophilic inflammation correlated with elevated IL-17A, proinflammatory cytokines (TNF-α, IL-1α, and IL-1β), neutrophil survival factors (granulocyte colony-stimulating factor), and neutrophil-targeting chemokines (CCL3 and CCL4). Eosinophilia correlated with elevated type 2 responses (IL-5 and IL-13) whereas inflammatory monocyte levels correlated with elevated type 1 responses (IFN-γ and CXCL9) and survival factors (macrophage colony-stimulating factor). Despite enhanced inflammatory responses, the immunoregulatory cytokine IL-10 and the natural inhibitor of IL-1 signaling, IL-1RA, were significantly elevated rather than impaired. Regulatory T-cell levels were unchanged, as were levels of the anti-inflammatory cytokines IL-35 and IL-38. Taken together, the CX3CL1/CX3CR1 axis preserves lung function during fungal-associated allergic airway inflammation through a nonclassical immunoregulatory mechanism.

Inherited causes of clonal haematopoiesis in 97,691 whole genomes

Age is the dominant risk factor for most chronic human diseases; yet the mechanisms by which aging confers this risk are largely unknown.¹ Recently, the age-related acquisition of somatic mutations in regenerating hematopoietic stem cell populations leading to clonal expansion was associated with both hematologic cancer ^2–4 and coronary heart disease⁵, a phenomenon termed ‘Clonal Hematopoiesis of Indeterminate Potential’ (CHIP).⁶ Simultaneous germline and somatic whole genome sequence analysis now provides the opportunity to identify root causes of CHIP. Here, we analyze high-coverage whole genome sequences from 97,691 participants of diverse ancestries in the NHLBI TOPMed program and identify 4,229 individuals with CHIP. We identify associations with blood cell, lipid, and inflammatory traits specific to different CHIP genes. Association of a genome-wide set of germline genetic variants identified three genetic loci associated with CHIP status, including one locus at TET2 that was African ancestry specific. In silico-informed in vitro evaluation of the TET2 germline locus identified a causal variant that disrupts a TET2 distal enhancer resulting in increased hematopoietic stem cell self-renewal. Overall, we observe that germline genetic variation shapes hematopoietic stem cell function leading to CHIP through mechanisms that are both specific to clonal hematopoiesis and shared mechanisms leading to somatic mutations across tissues.

Tumour necrosis factor-related apoptosis-inducing ligand (TRAIL) and receptors in type 1, type 2 and type 17 inflammation in cross-sectional asthma study

Tumour necrosis factor-related apoptosis-inducing ligand (TRAIL) reportedly promotes, or conversely, resolves inflammation in asthma. In this study of TRAIL and cell receptors in sputum, bronchoalveolar lavage and biopsy from subjects in the Severe Asthma Research Program at Wake Forest, the high TRAIL group had significant increases in all leucocytes, and was associated with increased type 1, type 2 and type 17 cytokines, but not type 9 interleukin 9. Two variants at loci in the TRAIL gene were associated with higher sputum levels of TRAIL. Increased TRAIL decoy receptor R3/DcR1 was observed on sputum leucocytes compared with death receptor R1/DR4, suggesting reduced apoptosis and prolonged cellular inflammation.

Chronic obstructive pulmonary disease and related phenotypes: polygenic risk scores in population-based and case-control cohorts

BACKGROUND

Genetic factors influence chronic obstructive pulmonary disease (COPD) risk, but the individual variants that have been identified have small effects. We hypothesised that a polygenic risk score using additional variants would predict COPD and associated phenotypes.

METHODS

We constructed a polygenic risk score using a genome-wide association study of lung function (FEV1 and FEV1/forced vital capacity [FVC]) from the UK Biobank and SpiroMeta. We tested this polygenic risk score in nine cohorts of multiple ethnicities for an association with moderate-to-severe COPD (defined as FEV1/FVC <0·7 and FEV1 <80% of predicted). Associations were tested using logistic regression models, adjusting for age, sex, height, smoking pack-years, and principal components of genetic ancestry. We assessed predictive performance of models by area under the curve. In a subset of studies, we also studied quantitative and qualitative CT imaging phenotypes that reflect parenchymal and airway pathology, and patterns of reduced lung growth.

FINDINGS

The polygenic risk score was associated with COPD in European (odds ratio [OR] per SD 1·81 [95% CI 1·74-1·88] and non-European (1·42 [1·34-1·51]) populations. Compared with the first decile, the tenth decile of the polygenic risk score was associated with COPD, with an OR of 7·99 (6·56-9·72) in European ancestry and 4·83 (3·45-6·77) in non-European ancestry cohorts. The polygenic risk score was superior to previously described genetic risk scores and, when combined with clinical risk factors (ie, age, sex, and smoking pack-years), showed improved prediction for COPD compared with a model comprising clinical risk factors alone (AUC 0·80 [0·79-0·81] vs 0·76 [0·75-0·76]). The polygenic risk score was associated with CT imaging phenotypes, including wall area percent, quantitative and qualitative measures of emphysema, local histogram emphysema patterns, and destructive emphysema subtypes. The polygenic risk score was associated with a reduced lung growth pattern.

INTERPRETATION

A risk score comprised of genetic variants can identify a small subset of individuals at markedly increased risk for moderate-to-severe COPD, emphysema subtypes associated with cigarette smoking, and patterns of reduced lung growth.

FUNDING

US National Institutes of Health, Wellcome Trust.

Overlap of Genetic Risk between Interstitial Lung Abnormalities and Idiopathic Pulmonary Fibrosis

Rationale: Interstitial lung abnormalities (ILAs) are associated with the highest genetic risk locus for idiopathic pulmonary fibrosis (IPF); however, the extent to which there are unique associations among individuals with ILAs or additional overlap with IPF is not known.Objectives: To perform a genome-wide association study (GWAS) of ILAs.Methods: ILAs and a subpleural-predominant subtype were assessed on chest computed tomography (CT) scans in the AGES (Age Gene/Environment Susceptibility), COPDGene (Genetic Epidemiology of Chronic Obstructive Pulmonary Disease [COPD]), Framingham Heart, ECLIPSE (Evaluation of COPD Longitudinally to Identify Predictive Surrogate End-points), MESA (Multi-Ethnic Study of Atherosclerosis), and SPIROMICS (Subpopulations and Intermediate Outcome Measures in COPD Study) studies. We performed a GWAS of ILAs in each cohort and combined the results using a meta-analysis. We assessed for overlapping associations in independent GWASs of IPF.Measurements and Main Results: Genome-wide genotyping data were available for 1,699 individuals with ILAs and 10,274 control subjects. The MUC5B (mucin 5B) promoter variant rs35705950 was significantly associated with both ILAs (P = 2.6 × 10^-27) and subpleural ILAs (P = 1.6 × 10^-29). We discovered novel genome-wide associations near IPO11 (rs6886640, P = 3.8 × 10^-8) and FCF1P3 (rs73199442, P = 4.8 × 10^-8) with ILAs, and near HTRE1 (rs7744971, P = 4.2 × 10^-8) with subpleural-predominant ILAs. These novel associations were not associated with IPF. Among 12 previously reported IPF GWAS loci, five (DPP9, DSP, FAM13A, IVD, and MUC5B) were significantly associated (P < 0.05/12) with ILAs.Conclusions: In a GWAS of ILAs in six studies, we confirmed the association with a MUC5B promoter variant and found strong evidence for an effect of previously described IPF loci; however, novel ILA associations were not associated with IPF. These findings highlight common genetically driven biologic pathways between ILAs and IPF, and also suggest distinct ones.