Although allele frequency data for most HLA loci provide strong evidence for balancing selection at the allele level, the DPB1 locus is a notable exception, with allele frequencies compatible with neutral evolution (genetic drift) or directional selection in most populations. This discrepancy is especially interesting as evidence for balancing selection has been seen at the nucleotide and amino acid (AA) sequence levels for DPB1. We describe methods used to examine the global distribution of DPB1 alleles and their constituent AA sequences. These methods allow investigation of the influence of natural selection in shaping DPβ diversity in a hierarchical fashion for DPB1 alleles, all polymorphic DPB1 exon 2-encoded AA positions, as well as all pairs and trios of these AA positions. In addition, we describe how asymmetric linkage disequilibrium for all DPB1 exon 2-encoded AA pairs can be used to complement other methods. Application of these methods provides strong evidence for the operation of balancing selection on AA positions 56, 85-87, 36, 55 and 84 (listed in decreasing order of the strength of selection), but no evidence for balancing selection on DPB1 alleles.

The DPB1 locus is notable among the classical HLA loci in that allele frequencies at this locus are consistent with genetic drift, whereas the frequencies of specific DPβ amino acids are consistent with the action of balancing selection. We investigated the influence of natural selection in shaping the diversity of three functional categories of DPB1 diversity defined by specific amino acid motifs, DPB1 T-cell epitopes, DPB1 supertypes and DP1-DP4 serologic categories (SCs), via Ewens-Watterson (EW) selective neutrality and asymmetric Linkage Disequilibrium (ALD) analyses in a worldwide sample of 136 populations. These EW analyses provide strong evidence for the operation of balancing selection on DP SCs, but no evidence for balancing selection on T-cell epitopes or supertypes. We further investigated the global distribution of SCs. Each SC is common in a different region of the world, with the DP1 SC most common in Southeast Asia and Oceania, the DP2 SC in North and South America, the DP3 SC in South America, and the DP4 SC in Europe. The DP2 SC is present in all populations, while 14% of populations are missing at least one DP1, DP3, or DP4 SC. We observed consistent DPA1∼DP SC haplotype associations across 10 populations from five global regions, and found that asymmetric linkage disequilibrium (LD) between the DPB1 locus and the four most-common DPA1 alleles (DPA1*01:03, *02:01, *02:02 and *03:01) is determined by variation at DPβ AA positions 85-87. These positions are in LD with both DPα positions 31 and 50. We conclude from these EW analyses that natural selection is primarily operating to maintain population-level diversity of DP SCs, rather than DPB1 alleles or other functional categories of DPB1 diversity.

Abstract Background Using microarray profiling of airway epithelial cells, we previously identified a Th2-high molecular phenotype of asthma based on expression of periostin, CLCA1 and serpinB2 and characterized by specific inflammatory, remodeling, and treatment response features. The goal of the current study was to develop a qPCR-based assay of Th2 inflammation to overcome the limitations of microarray-based methods. Methods Airway epithelial brushings were obtained by bronchoscopy from two clinical studies comprising 44 healthy controls and 62 subjects with asthma, 39 of whom were studied before and after a standardized 8 week course of inhaled corticosteroids (ICS). The qPCR-based expression of periostin, CLCA1 and serpinB2 were combined into a single metric. Results In asthma, the three-gene-mean of periostin, CLCA1 and serpinB2 correlated with FeNO (r = 0.75, p = 0.0002), blood eosinophils (r = 0.58, p = 0.003) and PC20 methacholine (r = -0.65, p = 0.0006), but not total serum IgE (r = 0.33, p = 0.1). Higher baseline three-gene-mean correlated with greater improvement in FEV1 with ICS at 2, 4 and 8 weeks (all p < 0.05). By ROC analysis, the area under the curve (AUC) of the three-gene-mean for FEV1 improvement with ICS at 4 and 8 weeks was 0.94 and 0.87, respectively, which are higher than the AUCs of FeNO, blood eosinophils, IgE or PC20. Th2 airway inflammation as measured by this three-gene-mean also had predictive capacity for an improvement in symptoms. Conclusions The three-gene-mean of periostin, CLCA1 and serpinB2 in airway epithelial brushings identifies Th2-high and low populations, is correlated with other Th2 biomarkers, and performs well for prediction of FEV1 improvement with ICS. The three-gene-mean provides a measurement of Th2 airway inflammation that is clinically relevant and that can serve as a valuable tool to evaluate non-invasive biomarkers to predict treatment responses to existing and emerging asthma therapies.

Background

Using microarray profiling of airway epithelial cells, we previously identified a Th2-high molecular phenotype of asthma based on expression of periostin, CLCA1 and serpinB2 and characterized by specific inflammatory, remodeling, and treatment response features. The goal of the current study was to develop a qPCR-based assay of Th2 inflammation to overcome the limitations of microarray-based methods.

Methods

Airway epithelial brushings were obtained by bronchoscopy from two clinical studies comprising 44 healthy controls and 62 subjects with asthma, 39 of whom were studied before and after a standardized 8 week course of inhaled corticosteroids (ICS). The qPCR-based expression of periostin, CLCA1 and serpinB2 were combined into a single metric.

Results

In asthma, the three-gene-mean of periostin, CLCA1 and serpinB2 correlated with FeNO (r = 0.75, p = 0.0002), blood eosinophils (r = 0.58, p = 0.003) and PC20 methacholine (r = -0.65, p = 0.0006), but not total serum IgE (r = 0.33, p = 0.1). Higher baseline three-gene-mean correlated with greater improvement in FEV1 with ICS at 2, 4 and 8 weeks (all p < 0.05). By ROC analysis, the area under the curve (AUC) of the three-gene-mean for FEV1 improvement with ICS at 4 and 8 weeks was 0.94 and 0.87, respectively, which are higher than the AUCs of FeNO, blood eosinophils, IgE or PC20. Th2 airway inflammation as measured by this three-gene-mean also had predictive capacity for an improvement in symptoms.

Conclusions

The three-gene-mean of periostin, CLCA1 and serpinB2 in airway epithelial brushings identifies Th2-high and low populations, is correlated with other Th2 biomarkers, and performs well for prediction of FEV1 improvement with ICS. The three-gene-mean provides a measurement of Th2 airway inflammation that is clinically relevant and that can serve as a valuable tool to evaluate non-invasive biomarkers to predict treatment responses to existing and emerging asthma therapies.

Rationale

Changes in airway epithelial cell differentiation, driven in part by IL-13, are important in asthma. Micro-RNAs (miRNAs) regulate cell differentiation in many systems and could contribute to epithelial abnormalities in asthma.

Objectives

To determine whether airway epithelial miRNA expression is altered in asthma and identify IL-13-regulated miRNAs.

Methods

We used miRNA microarrays to analyze bronchial epithelial brushings from 16 steroid-naive subjects with asthma before and after inhaled corticosteroids, 19 steroid-using subjects with asthma, and 12 healthy control subjects, and the effects of IL-13 and corticosteroids on cultured bronchial epithelial cells. We used quantitative polymerase chain reaction to confirm selected microarray results.

Measurements and main results

Most (12 of 16) steroid-naive subjects with asthma had a markedly abnormal pattern of bronchial epithelial miRNA expression by microarray analysis. Compared with control subjects, 217 miRNAs were differentially expressed in steroid-naive subjects with asthma and 200 in steroid-using subjects with asthma (false discovery rate < 0.05). Treatment with inhaled corticosteroids had modest effects on miRNA expression in steroid-naive asthma, inducing a statistically significant (false discovery rate < 0.05) change for only nine miRNAs. qPCR analysis confirmed differential expression of 22 miRNAs that were highly differentially expressed by microarrays. IL-13 stimulation recapitulated changes in many differentially expressed miRNAs, including four members of the miR-34/449 family, and these changes in miR-34/449 family members were resistant to corticosteroids.

Conclusions

Dramatic alterations of airway epithelial cell miRNA levels are a common feature of asthma. These alterations are only modestly corrected by inhaled corticosteroids. IL-13 effects may account for some of these alterations, including repression of miR-34/449 family members that have established roles in airway epithelial cell differentiation. Clinical trial registered with www.clinicaltrials.gov (NCT 00595153).

Background

Eosinophilic airway inflammation is heterogeneous in asthmatic patients. We recently described a distinct subtype of asthma defined by the expression of genes inducible by T(H)2 cytokines in bronchial epithelium. This gene signature, which includes periostin, is present in approximately half of asthmatic patients and correlates with eosinophilic airway inflammation. However, identification of this subtype depends on invasive airway sampling, and hence noninvasive biomarkers of this phenotype are desirable.

Objective

We sought to identify systemic biomarkers of eosinophilic airway inflammation in asthmatic patients.

Methods

We measured fraction of exhaled nitric oxide (Feno), peripheral blood eosinophil, periostin, YKL-40, and IgE levels and compared these biomarkers with airway eosinophilia in asthmatic patients.

Results

We collected sputum, performed bronchoscopy, and matched peripheral blood samples from 67 asthmatic patients who remained symptomatic despite maximal inhaled corticosteroid treatment (mean FEV(1), 60% of predicted value; mean Asthma Control Questionnaire [ACQ] score, 2.7). Serum periostin levels are significantly increased in asthmatic patients with evidence of eosinophilic airway inflammation relative to those with minimal eosinophilic airway inflammation. A logistic regression model, including sex, age, body mass index, IgE levels, blood eosinophil numbers, Feno levels, and serum periostin levels, in 59 patients with severe asthma showed that, of these indices, the serum periostin level was the single best predictor of airway eosinophilia (P = .007).

Conclusion

Periostin is a systemic biomarker of airway eosinophilia in asthmatic patients and has potential utility in patient selection for emerging asthma therapeutics targeting T(H)2 inflammation.

Mucous cell hyperplasia and airway smooth muscle (ASM) hyperresponsiveness are hallmark features of inflammatory airway diseases, including asthma. Here, we show that the recently identified calcium-activated chloride channel (CaCC) TMEM16A is expressed in the adult airway surface epithelium and ASM. The epithelial expression is increased in asthmatics, particularly in secretory cells. Based on this and the proposed functions of CaCC, we hypothesized that TMEM16A inhibitors would negatively regulate both epithelial mucin secretion and ASM contraction. We used a high-throughput screen to identify small-molecule blockers of TMEM16A-CaCC channels. We show that inhibition of TMEM16A-CaCC significantly impairs mucus secretion in primary human airway surface epithelial cells. Furthermore, inhibition of TMEM16A-CaCC significantly reduces mouse and human ASM contraction in response to cholinergic agonists. TMEM16A-CaCC blockers, including those identified here, may positively impact multiple causes of asthma symptoms.

Rationale

Quantification of type 2 inflammation provided a molecular basis for heterogeneity in asthma. Non-type 2 pathways that contribute to asthma pathogenesis are not well understood.

Objectives

To identify dysregulated pathways beyond type 2 inflammation.

Methods

We applied RNA sequencing to airway epithelial brushings obtained from subjects with stable mild asthma not on corticosteroids (n = 19) and healthy control subjects (n = 16). Sequencing reads were mapped to human and viral genomes. In the same cohort, and in a separate group with severe asthma (n = 301), we profiled blood gene expression with microarrays.

Measurements and main results

In airway brushings from mild asthma on inhaled corticosteroids, RNA sequencing yielded 1,379 differentially expressed genes (false discovery rate < 0.01). Pathway analysis revealed increased expression of type 2 markers, IFN-stimulated genes (ISGs), and endoplasmic reticulum (ER) stress-related genes. Airway epithelial ISG expression was not associated with type 2 inflammation in asthma or with viral transcripts but was associated with reduced lung function by FEV₁ (ρ = -0.72; P = 0.0004). ER stress was confirmed by an increase in XBP1 (X-box binding protein 1) splicing in mild asthma and was associated with both type 2 inflammation and ISG expression. ISGs were also the most activated genes in blood cells in asthma and were correlated with airway ISG expression (ρ = 0.55; P = 0.030). High blood ISG expression in severe asthma was similarly unrelated to type 2 inflammation.

Conclusions

ISG activation is prominent in asthma, independent of viral transcripts, orthogonal to type 2 inflammation, and associated with distinct clinical features. ER stress is associated with both type 2 inflammation and ISG expression.