Regulation of T cells in airway disease by beta-agonist

Circulating microRNAs associated with bronchodilator response in childhood asthma

BACKGROUND

Bronchodilator response (BDR) is a measure of improvement in airway smooth muscle tone, inhibition of liquid accumulation and mucus section into the lumen in response to short-acting beta-2 agonists that varies among asthmatic patients. MicroRNAs (miRNAs) are well-known post-translational regulators. Identifying miRNAs associated with BDR could lead to a better understanding of the underlying complex pathophysiology.

OBJECTIVE

The purpose of this study is to identify circulating miRNAs associated with bronchodilator response in asthma and decipher possible mechanism of bronchodilator response variation.

METHODS

We used available small RNA sequencing on blood serum from 1,134 asthmatic children aged 6 to 14 years who participated in the Genetics of Asthma in Costa Rica Study (GACRS). We filtered the participants into the highest and lowest bronchodilator response (BDR) quartiles and used DeSeq2 to identify miRNAs with differential expression (DE) in high (N = 277) vs. low (N = 278) BDR group. Replication was carried out in the Leukotriene modifier Or Corticosteroids or Corticosteroid-Salmeterol trial (LOCCS), an adult asthma cohort. The putative target genes of DE miRNAs were identified, and pathway enrichment analysis was performed.

RESULTS

We identified 10 down-regulated miRNAs having odds ratios (OR) between 0.37 and 0.76 for a doubling of miRNA counts and one up-regulated miRNA (OR = 2.26) between high and low BDR group. These were assessed for replication in the LOCCS cohort, where two miRNAs (miR-200b-3p and miR-1246) were associated. Further, functional annotation of 11 DE miRNAs were performed as well as of two replicated miRs. Target genes of these miRs were enriched in regulation of cholesterol biosynthesis by SREBPs, ESR-mediated signaling, G1/S transition, RHO GTPase cycle, and signaling by TGFB family pathways.

CONCLUSION

MiRNAs miR-1246 and miR-200b-3p are associated with both childhood and adult asthma BDR. Our findings add to the growing body of evidence that miRNAs play a significant role in the difference of asthma treatment response among patients as it points to genomic regulatory machinery underlying difference in bronchodilator response among patients.

TRIAL REGISTRATION

LOCCS cohort [ClinicalTrials.gov number NCT00156819, Registration date 20050912], GACRS cohort [ClinicalTrials.gov number NCT00021840].

Circulating MicroRNAs associated with Bronchodilator Response in Childhood Asthma

Rationale

Bronchodilator response (BDR) is a measure of improvement in airway smooth muscle tone, inhibition of liquid accumulation and mucus section into the lumen in response to short-acting beta-2 agonists that varies among asthmatic patients. MicroRNAs (miRNAs) are well-known post-translational regulators. Identifying miRNAs associated with BDR could lead to a better understanding of the underlying complex pathophysiology.

Objective

The purpose of this study is to identify circulating miRNAs associated with bronchodilator response in asthma and decipher possible mechanism of bronchodilator response variation.

Methods

We used available small RNA sequencing on blood serum from 1,134 asthmatic children aged 6 to 14 years who participated in the Genetics of Asthma in Costa Rica Study (GACRS). We filtered the participants into high and low bronchodilator response (BDR) quartiles and used DeSeq2 to identify miRNAs with differential expression (DE) in high (N= 277) vs low (N= 278) BDR group. Replication was carried out in the Leukotriene modifier Or Corticosteroids or Corticosteroid-Salmeterol trial (LOCCS), an adult asthma cohort. The putative target genes of DE miRNAs were identified, and pathway enrichment analysis was performed.

Results

We identified 10 down-regulated miRNAs having odds ratios (OR) between 0.37 and 0.76 for a doubling of miRNA counts and one up-regulated miRNA (OR=2.26) between high and low BDR group. These were assessed for replication in the LOCCS cohort, where two miRNAs (miR-200b-3p and miR-1246) were associated. Further, functional annotation of 11 DE miRNAs were performed as well as of two replicated miRs. Target genes of these miRs were enriched in regulation of cholesterol biosynthesis by SREBPs, ESR-mediated signaling, G1/S transition, RHO GTPase cycle, and signaling by TGFB family pathways.

Conclusion

MiRNAs miR-1246 and miR-200b-3p are associated with both childhood and adult asthma BDR. Our findings add to the growing body of evidence that miRNAs play a significant role in the difference of asthma treatment response among patients as it points to genomic regulatory machinery underlying difference in bronchodilator response among patients.

Trial registration

LOCCS cohort [ClinicalTrials.gov number: NCT00156819], GACRS cohort [ClinicalTrials.gov number: NCT00021840].

Integration of Transcriptomic Data Identifies Global and Cell-Specific Asthma-Related Gene Expression Signatures

Over 140,000 transcriptomic studies performed in healthy and diseased cell and tissue types, at baseline and after exposure to various agents, are available in public repositories. Integrating results of transcriptomic datasets has been an attractive approach to identify gene expression signatures that are more robust than those obtained for individual datasets, especially datasets with small sample size. We developed Reproducible Analysis and Validation of Expression Data (RAVED), a pipeline that facilitates the creation of R Markdown reports detailing reproducible analysis of publicly available transcriptomic data, and used it to analyze asthma and glucocorticoid response microarray and RNA-Seq datasets. Subsequently, we used three approaches to integrate summary statistics of these studies and identify cell/tissue-specific and global asthma and glucocorticoid-induced gene expression changes. Transcriptomic integration methods were incorporated into an online app called REALGAR, where end-users can specify datasets to integrate and quickly obtain results that may facilitate design of experimental studies.

A genome-wide association and admixture mapping study of bronchodilator drug response in African Americans with asthma

Short-acting β₂-adrenergic receptor agonists (SABAs) are the most commonly prescribed asthma medications worldwide. Response to SABAs is measured as bronchodilator drug response (BDR), which varies among racial/ethnic groups in the U.S^{1, 2}. However, the genetic variation that contributes to BDR is largely undefined in African Americans with asthma³. To identify genetic variants that may contribute to differences in BDR in African Americans with asthma, we performed a genome-wide association study (GWAS) of BDR in 949 African American children with asthma, genotyped with the Axiom World Array 4 (Affymetrix, Santa Clara, CA) followed by imputation using 1000 Genomes phase III genotypes. We used linear regression models adjusting for age, sex, body mass index (BMI) and genetic ancestry to test for an association between BDR and genotype at single nucleotide polymorphisms (SNPs). To increase power and distinguish between shared vs. population-specific associations with BDR in children with asthma, we performed a meta-analysis across 949 African Americans and 1,830 Latinos (Total=2,779). Lastly, we performed genome-wide admixture mapping to identify regions whereby local African or European ancestry is associated with BDR in African Americans. We identified a population-specific association with an intergenic SNP on chromosome 9q21 that was significantly associated with BDR (rs73650726, p=7.69×10⁻⁹). A trans-ethnic meta-analysis across African Americans and Latinos identified three additional SNPs within the intron of PRKG1 that were significantly associated with BDR (rs7903366, rs7070958, and rs7081864, p≤5×10⁻⁸). Our results failed to replicate in three additional populations of 416 Latinos and 1,615 African Americans. Our findings indicate that both population specific and shared genetic variation contributes to differences in BDR in minority children with asthma, and that the genetic underpinnings of BDR may differ between racial/ethnic groups.

Defective Resensitization in Human Airway Smooth Muscle Cells Evokes β-Adrenergic Receptor Dysfunction in Severe Asthma

β₂-adrenergic receptor (β₂AR) agonists (β₂-agonist) are the most commonly used therapy for acute relief in asthma, but chronic use of these bronchodilators paradoxically exacerbates airway hyper-responsiveness. Activation of βARs by β-agonist leads to desensitization (inactivation) by phosphorylation through G-protein coupled receptor kinases (GRKs) which mediate β-arrestin binding and βAR internalization. Resensitization occurs by dephosphorylation of the endosomal βARs which recycle back to the plasma membrane as agonist-ready receptors. To determine whether the loss in β-agonist response in asthma is due to altered βAR desensitization and/or resensitization, we used primary human airway smooth muscle cells (HASMCs) isolated from the lungs of non-asthmatic and fatal-asthmatic subjects. Asthmatic HASMCs have diminished adenylyl cyclase activity and cAMP response to β-agonist as compared to non-asthmatic HASMCs. Confocal microscopy showed significant accumulation of phosphorylated β₂ARs in asthmatic HASMCs. Systematic analysis of desensitization components including GRKs and β-arrestin showed no appreciable differences between asthmatic and non-asthmatic HASMCs. However, asthmatic HASMC showed significant increase in PI3Kγ activity and was associated with reduction in PP2A activity. Since reduction in PP2A activity could alter receptor resensitization, endosomal fractions were isolated to assess the agonist ready β₂ARs as a measure of resensitization. Despite significant accumulation of β₂ARs in the endosomes of asthmatic HASMCs, endosomal β₂ARs cannot robustly activate adenylyl cyclase. Furthermore, endosomes from asthmatic HASMCs are associated with significant increase in PI3Kγ and reduced PP2A activity that inhibits β₂AR resensitization. Our study shows that resensitization, a process considered to be a homeostasis maintaining passive process is inhibited in asthmatic HASMCs contributing to β₂AR dysfunction which may underlie asthma pathophysiology and loss in asthma control.

Distinct PKA and Epac compartmentalization in airway function and plasticity

Asthma and chronic obstructive pulmonary disease (COPD) are obstructive lung diseases characterized by airway obstruction, airway inflammation and airway remodelling. Next to inflammatory cells and airway epithelial cells, airway mesenchymal cells, including airway smooth muscle cells and (myo)fibroblasts, substantially contribute to disease features by the release of inflammatory mediators, smooth muscle contraction, extracellular matrix deposition and structural changes in the airways. Current pharmacological treatment of both diseases intends to target the dynamic features of the endogenous intracellular suppressor cyclic AMP (cAMP). This review will summarize our current knowledge on cAMP and will emphasize on key discoveries and paradigm shifts reflecting the complex spatio-temporal nature of compartmentalized cAMP signalling networks in health and disease. As airway fibroblasts and airway smooth muscle cells are recognized as central players in the development and progression of asthma and COPD, we will focus on the role of cAMP signalling in their function in relation to airway function and plasticity. We will recapture on the recent identification of cAMP-sensing multi-protein complexes maintained by cAMP effectors, including A-kinase anchoring proteins (AKAPs), proteins kinase A (PKA), exchange protein directly activated by cAMP (Epac), cAMP-elevating seven-transmembrane (7TM) receptors and phosphodiesterases (PDEs) and we will report on findings indicating that the pertubation of compartmentalized cAMP signalling correlates with the pathopysiology of obstructive lung diseases. Future challenges include studies on cAMP dynamics and compartmentalization in the lung and the development of novel drugs targeting these systems for therapeutic interventions in chronic obstructive inflammatory diseases.

Modulation of T cell immune functions by the prostaglandin E(2) - cAMP pathway in chronic inflammatory states

Cyclic AMP is the intracellular second messenger for a variety of immunoregulatory inflammatory mediators such as prostaglandin E2, adenosine and histamine that signal to effector T cells from monocytes, macrophages and regulatory T cells. Protein kinase A (PKA) type I localizes to lipid rafts in effector T cells during T cell activation and directly modulates proximal signal events including phosphorylation of C-terminal Src kinase (Csk), which initiates a negative signal pathway that fine-tunes the T cell activation process. The PKA-Csk immunoregulatory pathway is scaffolded by the A kinase anchoring protein ezrin, the Csk binding protein phosphoprotein associated with glycosphingolipid-enriched membrane microdomains and the linker protein ezrin/radixin/moesin binding protein of 50 kDa. This pathway is hyperactivated in chronic infections with an inflammatory component such as HIV, other immunodeficiencies and around solid tumours as a consequence of local inflammation leading to inhibition of anti-tumour immunity. LINKED ARTICLES This article is part of a themed section on Novel cAMP Signalling Paradigms. To view the other articles in this section visit http://dx.doi.org/10.1111/bph.2012.166.issue-2.

Genome-wide association analysis in asthma subjects identifies SPATS2L as a novel bronchodilator response gene

Bronchodilator response (BDR) is an important asthma phenotype that measures reversibility of airway obstruction by comparing lung function (i.e. FEV₁) before and after the administration of a short-acting β₂-agonist, the most common rescue medications used for the treatment of asthma. BDR also serves as a test of β₂-agonist efficacy. BDR is a complex trait that is partly under genetic control. A genome-wide association study (GWAS) of BDR, quantified as percent change in baseline FEV₁ after administration of a β₂-agonist, was performed with 1,644 non-Hispanic white asthmatic subjects from six drug clinical trials: CAMP, LOCCS, LODO, a medication trial conducted by Sepracor, CARE, and ACRN. Data for 469,884 single-nucleotide polymorphisms (SNPs) were used to measure the association of SNPs with BDR using a linear regression model, while adjusting for age, sex, and height. Replication of primary P-values was attempted in 501 white subjects from SARP and 550 white subjects from DAG. Experimental evidence supporting the top gene was obtained via siRNA knockdown and Western blotting analyses. The lowest overall combined P-value was 9.7E-07 for SNP rs295137, near the SPATS2L gene. Among subjects in the primary analysis, those with rs295137 TT genotype had a median BDR of 16.0 (IQR = [6.2, 32.4]), while those with CC or TC genotypes had a median BDR of 10.9 (IQR = [5.0, 22.2]). SPATS2L mRNA knockdown resulted in increased β₂-adrenergic receptor levels. Our results suggest that SPATS2L may be an important regulator of β₂-adrenergic receptor down-regulation and that there is promise in gaining a better understanding of the biological mechanisms of differential response to β₂-agonists through GWAS.

Bronchodilator response (BDR) is an important asthma phenotype that measures reversibility of airway obstruction by comparing lung function before and after the administration of short-acting β₂-agonists, common medications used for asthma treatment. We performed a genome-wide association study of BDR with 1,644 white asthmatic subjects from six drug clinical trials and attempted to replicate these findings in 1,051 white subjects from two independent cohorts. The most significant associated variant was near the SPATS2L gene. We knocked down SPATS2L mRNA in human airway smooth muscle cells and found that β₂-adrenergic receptor levels increased, suggesting that SPATS2L may be a regulator of BDR. Our results highlight the promise of pursuing GWAS results that do not necessarily reach genome-wide significance and are an example of how results from pharmacogenetic GWAS can be studied functionally.

Norepinephrine suppresses IFN-γ and TNF-α production by murine intestinal intraepithelial lymphocytes via the β₁ adrenoceptor

We examined whether norepinephrine (NE) had direct effects on cytokine production by murine intestinal intraepithelial lymphocytes (IELs), compared with splenocytes. CD3⁺ IELs and CD3⁺ splenocytes expressed α(1B), α(1D), α(2C), β₁, β₂, and β₃ adrenoceptors (ARs). NE significantly suppressed IFN-γ and TNF-α production by IELs and splenocytes ex vivo. The suppressive effects of NE in IELs were reversed by β₁ AR antagonist CGP-20712A, whereas those in splenocytes were reversed by β₂ AR antagonist ICI118,551. In IELs, β₁ AR agonist xamoterol mimicked the suppressive effects of NE. These results indicated NE regulates intestinal mucosal immune responses mediated by IELs via β₁ AR.

Combined anti-inflammatory effects of β2-adrenergic agonists and PDE4 inhibitors on astrocytes by upregulation of intracellular cAMP

Inflammation is an important hallmark of all neurodegenerative diseases and activation of different glial populations may be involved in the progression of some of these disorders. Especially, the activation of astroglia can lead to long-term detrimental morphological changes, such as scar formation. Therefore, improved strategies to modulate inflammation in these cells are currently being investigated. We investigated the interaction of phosphodiesterase (PDE) 4 inhibitors, such as rolipram, with other agents raising cellular cAMP levels. When used alone, none of the PDE4 inhibitors increased cAMP levels. The adenylate cyclase activator forskolin, the β(2)-adrenergic agonist clenbuterol and the mixed β(1)/β(2)-adrenergic agonist isoproterenol increased intracellular cAMP levels of cortical murine astrocytes. This increase was synergistically elevated by rolipram or the PDE4 inhibitor RO-201724, but not by inhibition of PDE3. Inflammatory stimulation of the cells with the cytokines TNF-α, IL-1β and IFN-γ strongly induced PDE4B and augmented overall PDE4 activity, while PDE3 activity was low. Clenbuterol and forskolin caused downregulation of cytokines and chemokines such as IL-6 and MCP-1. This effect was further enhanced by rolipram, but not by the PDE3 inhibitor milrinone. The cAMP-raising drug combinations attenuated the upregulation of TNF-α and IL-6 mRNA and the secretion of IL-6, but did not affect initial NF-κB signalling triggered by the stimulating cytokines. These results indicate that PDE4 may be a valuable anti-inflammatory target in brain diseases, especially under conditions associated with stimulation of cAMP-augmenting astrocyte receptors as is observed by clenbuterol treatment.