Asthma is a chronic respiratory disease usually characterized by inflammatory changes in the lower airways that requires long-term management.
- severe asthma
- genetics
- pharmacogenetics
- precision medicine
1. Introduction
Asthma is a chronic respiratory disease usually characterized by inflammatory changes in the lower airways that requires long-term management. It has emerged as a major public health condition estimated to affect over than 339 million people worldwide and influenced by a puzzling interaction of environmental and genetic factors [1][2]. The rising prevalence of asthma in industrialized and developing countries over the past few decades has highlighted the importance of understanding the biological mechanisms that contribute to asthma development [3].
Asthma is a chronic respiratory disease usually characterized by inflammatory changes in the lower airways that requires long-term management. It has emerged as a major public health condition estimated to affect over than 339 million people worldwide and influenced by a puzzling interaction of environmental and genetic factors [1,2]. The rising prevalence of asthma in industrialized and developing countries over the past few decades has highlighted the importance of understanding the biological mechanisms that contribute to asthma development [3].
Asthma is a complex syndrome with a variable natural history and clinical development influenced by multiple genes and environmental exposures [2][4][5]. As illustrated in
Asthma is a complex syndrome with a variable natural history and clinical development influenced by multiple genes and environmental exposures [2,4,5]. As illustrated in
, the severe asthma phenotype is expressed by diverse genetic predispositions that may favor abnormal pulmonary growth with lower values of lung function, lower respiratory tract infections in early life, bronchial hyperresponsiveness and airway inflammation, and delayed immune system maturation and allergic sensitization in childhood [4]. Moreover, phenotypic expression and the epigenetics pressure in some genes might have a significant influence on specific clusters of patients, for instance, in those with occupational asthma, smoking asthmatics, or those who were exposed to allergens in the first few years of life [6].
Figure 1.
Severe asthma patients require close monitoring with typically challenging clinical management and a considerable financial burden. Despite pharmacologic treatment with anti-inflammatory and bronchodilators, a cluster of patients remain poorly controlled, and at increased risk of airway remodeling and disease progression [4][7].
Severe asthma patients require close monitoring with typically challenging clinical management and a considerable financial burden. Despite pharmacologic treatment with anti-inflammatory and bronchodilators, a cluster of patients remain poorly controlled, and at increased risk of airway remodeling and disease progression [4,7].
The goal of understanding the genetic and epigenetic basis that affect the therapeutical response in severe asthma is to assist the optimization of pharmacological treatment as well as the development of potential new target-guided therapy [4]. Pharmacogenomics may play a key role in the development of therapeutic options that better suit each individual’s unique needs.
Treatment response for asthma can be assessed through different clinical parameters, such as lung function measurements, quality of life, asthma symptoms and exacerbations, but also by other indicators of failure of the treatment, such as: (1) hospital admissions for asthma, and (2) a course of oral corticosteroids or (3) emergency room visits because of asthma. Moreover, several factors can influence poor outcomes—among them, age, continued exposure to allergens and pollutants, and genetic and epigenetic factors.
2. Pharmacogenomics of the Therapeutic Response of Asthma
It is well-known that the therapeutic response in asthma is heterogeneous, which can be at least partially explained by interindividual variability. Over the past few years, researchers have evaluated these associations through hypothesis-free strategies, such as the Genome-Wide Association Study (GWAS), which has uncovered novel pharmacogenetic variants related to inhaled corticosteroids (ICS) and bronchodilator responses. Samedy-Bates and collaborators [5] evaluated the association between ICS use and the bronchodilator response (BDR) in three pediatric populations with persistent asthma (African American, Puerto Rican, and Mexican American children). ICS treatment was significantly associated with increased BDR among Mexican Americans. However, no difference was found in the other two populations, indicating that the benefit of using ICS may vary between racial/ethnic groups, as a consequence of the genetic difference between populations.A systematic review published in 2019 analyzed a subset of seven GWAS on the inhaled corticosteroid response [8]. It reports that few genes’ effects have been replicated in GWAS or even in candidate gene studies. Such an observation can be explained by lack of harmonization between the analyzed phenotypes or even the genetic architecture of the studied populations. According to the systematic review, the most replicated variants were rs242941 and rs1876828 in CRHR1 [9], the rs37973 in GLCCI1 [9][10][11][12], the rs28364072 in FCER2 [12][13][14][15][16], the rs2240017 in TBX21 [9][12][13][14][15][16], the rs41423247 in NR3C1 [9], and variants in the 17q21 locus [17].
A systematic review published in 2019 analyzed a subset of seven GWAS on the inhaled corticosteroid response [8]. It reports that few genes’ effects have been replicated in GWAS or even in candidate gene studies. Such an observation can be explained by lack of harmonization between the analyzed phenotypes or even the genetic architecture of the studied populations. According to the systematic review, the most replicated variants were rs242941 and rs1876828 in CRHR1 [9], the rs37973 in GLCCI1 [9,10,11,12], the rs28364072 in FCER2 [12,13,14,15,16], the rs2240017 in TBX21 [9,12,13,14,15,16], the rs41423247 in NR3C1 [9], and variants in the 17q21 locus [17]. In the period considered for the present review, only three GWAS were published and may contribute to new loci associated with a poor response to inhaled corticosteroids. In the first, Hernandez–Pacheco et al. [18] carried out a meta-analysis involving two GWAS of asthma exacerbations in Hispanic and African American children treated with ICS. No variant was associated with the level of statistical significance required for a GWAS. However, 11 variants were suggestively associated with asthma exacerbations. Of these, only the variant rs5995653, at the intergenic region of APOBEC3B and APOBEC3C, was replicated in another population studied. A second GWAS carried out by the same group evaluated asthma exacerbations in a different population, including children of European descent treated with ICS from eight studies [19]. Ten variants were associated with asthma exacerbations in patients treated with ICS. Moreover, one variant, rs67026078, located within the intergenic region ofCACNA2D3
andWNT5A,
have been replicated among Europeans, but not validated in non-European populations. The same group performed a combined analysis of transcriptomic and genetic data in relation to treatment responses among individuals with asthma. The LTBP1 was the most consistently associated gene in a favorable ICS response among asthma patients. Analysis of this gene revealed the rs11681246 was negatively associated with asthma exacerbations regardless of ICS use in Europeans but not in non-Europeans. Conversely, the rs76390075, also located in LTBP1, was negatively associated with the same outcome in the non-European population (Table 1) [20]. Similar to the response to ICS, the bronchodilator drug response (BDR) also varies among population groups as a consequence of the genetic variability. As demonstrated by Spear and colleagues [21] in African American children with asthma, the rs73650726, located in the non-coding RNA LOC105376110, was negatively associated to BDR with statistical significance at the GWAS (Table 1). Additionally, the other three single-nucleotide polymorphisms (SNPs) in the PRKG1 (Protein Kinase, CGMP-Dependent, Type I) were associated with the BDR in the meta-analysis across African Americans and Latin Americans living in the USA (rs7903366, rs7070958, and rs7081864). Recently, a meta-analysis of Genome-Wide Association Studies (meta-GWAS) of asthma exacerbations on regular long-acting beta 2 agonist (LABA) use was performed in children and young adults with asthma from six studies. No variant was associated with the outcome with the GWAS significance level. However, eight variants were suggestively associated with a poor response to regular LABA treatment (p-value threshold ≤ 5 × 10−6
), including variants in two loci (TBX3
andEPHA7)
previously implicated in the response to short-acting beta2-agonists (SABA) [22]. Further studies in larger populations with a diverse genetic structure will be required to validate these results.3. Candidate Genes Studies Linked to Therapeutic Response of Asthma
3.1. Influence of Genetic Variant in ADRB2 on Response to the Bronchodilator
A systematic review has reported that the ADRB2 encoding the beta2 adrenergic receptor is the most-studied gene on the pharmacogenetics of the bronchodilator response [23]. Several studies have found associations with three non-synonymous genetic variants, rs1042713 (Arg16Gly), rs1042714 (Gln27Glu), and rs180888 (Thr164Ile) with the BDR in children and adults with asthma [24][25][26][27][28][29] although these results are still controversial [30][31][32][33] (
Table 1.
| SNP | Score PharmGKB | Gene | Outcome | Reference |
|---|---|---|---|---|
| rs11681246 | - | LTBP1 | Negatively associated with asthma exacerbations regardless of ICS use | [20] |
| rs76390075 | - | LTBP1 | [20] | |
| rs73650726 | - | LOC105376110 | Negatively associated with BDR | [21] |
| rs7903366 | - | PRKG1 | Positively associated with the BDR | [21] |
| rs7081864 | - | PRKG1 | Positively associated with the BDR | [21] |
| rs7070958 | - | PRKG1 | Negatively associated with the BDR | [21] |
| rs1042713 | 2A | ADRB2 | Better response to treatment with A allele | [24][30][24,30] |
| Worst response to treatment with A allele | [29] | |||
| No association with therapeutic response | [25][32][34][25,32,34] | |||
| rs1042714 | - | ADRB2 | Better response to treatment | [23][23[29],29] |
| No association with therapeutic response | [24][27][30][34][24,27,30,34] | |||
| rs1042713/rs1042714 | - | ADRB2 | Better response to treatment | [35] |
| rs180888 | - | ADRB2 | Uncontrolled asthma during LABA treatment | [28] |
| rs37972 | - | GLCCI1 | Attenuation of the response to ICS treatment | [36][41] |
| rs37973 | 3 | GLCCI1 | Worst response to treatment with G Allele | [37][38][36][38,40,41] |
| Better response to treatment with G allele | [39][42] | |||
| rs2228570 | - | VDR | Risk of resistance to inhaled glucocorticoids | [40][46] |
| Higher exacerbation severity scores and poorer β2-agonist treatment response | [41][47] |
1
1
A study carried out by Bhosale et al. [42] investigated the molecular impact of the rs1042714 variant on the structure of the β2AR and binding of the beta agonists albuterol, isoproterenol, and terbutaline. It was observed that the G allele (Glu variant) alters electrostatic interactions in the region of the binding site which results in better bonding with the three ligands tested compared to the wild type, which can explain the variable therapeutic response to adrenergic agonists.
®
https://www.pharmgkb.org/ accessed 13 Mar 2021) was consulted. This is a NIH-funded resource managed at Stanford University that provides information about how genetic variants affect the therapeutic response, assigning levels from 1–4, with level 1 meeting the highest criteria [43]. Among the pharmacogenetic variants for the therapeutic response of asthma, the variant rs1042713 has the highest score on the PharmGKB
®
3.2. Variants in GLCCI1 Associated with Response to Corticosteroids
Current guidelines recommend ICS as first-line therapy in persistent asthma as its anti-inflammatory effect reduces airway inflammation and risk of exacerbations, and improves disease control [37]. Despite the proven clinical efficacy of ICS, the intraindividual response is variable and likely influenced by genetic background and environmental exposures [44]. GLCCI1 is involved in inflammatory cell apoptosis and its expression is induced by glucocorticoids [38].
A GWAS carried out by Tantisira and colleagues [36] has shown that genetic variations in glucocorticoid-induced transcript 1 genes (GLCCI1) are linked with decreased lung function in subjects treated with ICS. The single nucleotide variant (SNV) rs37972 was associated with changes in FEV
1
1 decline of 30 mL/year or greater [38]. Conversely, however, a Slovenian study reported that rs37973 SNP was positively associated with ICS responsiveness, even though it was considerably modified by smoking and atopy [39] (
The association between genetic variants and exacerbation risk was recently addressed in the largest pharmacogenetic candidate gene study to date in the Dutch population-based Rotterdam study with replication in the American GERA cohort [37]. Single-nucleotide polymorphisms for GLCCI1 (rs37973) significantly increased exacerbation risk in patients treated with ICS in the Dutch cohort, but an opposite effect was seen in the American replication cohort. The rs37973 is level 3 on PharmGKB
®
3.3. Genetic Variants in the Vitamin D Pathway and the Therapeutic Response on Asthma
The active metabolite of vitamin D, 1α-25-dihydroxyvitamin D, has a substantial immunomodulatory influence on the expression of proinflammatory transcriptional regulators, the inflammatory cellular response, cytokines, and IgE regulation. Variants in the vitamin D pathway gene (VDR, CYP2R1, CYP24A1, CYP27B1, DBP) are related, typically, to molecular pathways that affect the phenotype of asthma [45].
A study with steroid-resistant asthma patients demonstrated that the oral treatment with calcitriol, the active form of vitamin D, augments IL-10 production and reduces IL-17A production from peripheral blood mononuclear cells after exposure to dexamethasone, when compared to patients with steroid-sensitive asthma [46]. Additionally, budesonide and calcitriol up-regulate the expression of VDR and Smad7 and inhibit the production of proteins involved in airway remodeling. These effects are enhanced when both (budesonide and calcitriol) are used in combination [47]. Together, these findings highlight the likely importance of the vitamin D pathway in the response to corticosteroids.
Mohamed and collaborators [40] evaluated the association between variants in the Vitamin D receptor (VDR) gene, and the response to glucocorticoids revealed a three times greater risk of individuals with the T allele of the rs2228570 (FokI) variant to be resistant to glucocorticoids. In another study evaluating the same variant, children with asthma carrying the polymorphic allele had higher exacerbation severity scores and a poorer β
2-agonist treatment response [41] (
4. Interaction between Genetics and Age in Response to Asthma Treatment
An aspect seldom addressed in genetic studies is the impact of age on genetic predisposition. In this context, a study by Dahlin and colleagues [48] performed a genome-wide interaction study (GWIS) to evaluate, for the first time, the interaction of genetic variation and age on ICS response (measured by the occurrence of exacerbations) on 1321 adults and children of European ancestry with asthma. The top SNP was the rs34631960 in THSD4 (thrombospondin type 1 domain containing protein 4) that was associated with risk of poor ICS response with increasing age (P < 5 × 10−8
). Moreover, the rs2328386 in HIVEP2 (human immunodeficiency virus type I enhancer binding protein 2) was a protective factor for the lack of response to ICS with increasing age. Another 107 age-by-genotype interaction met genome-wide suggestive significance. A similar analysis was carried out between the interaction of genetic factors and age on the response to the bronchodilator from childhood to adulthood with asthma. Sordillo et al. [49] demonstrated for the first time that the rs295137 and the rs2626393 had their association with the BDR modified by age. In early childhood through to adolescence, the rs295137 variant (T allele) near the SPATS2L gene was associated with an increased BDR response, while the rs2626393 variant (C allele) near the ASB3 gene was associated with lower BDR with a significant decrease in the magnitude of effect from adolescence to adulthood. These findings can be explained by the impact of the accumulation of environmental exposures on the genetic predisposition, determined by epigenetic changes.5. Epigenetic Mechanisms Involved in the Lack of Therapeutic Control of Asthma
Epigenetic mechanisms mediate nonstructural changes in gene expression that are regulated by DNA methylation patterns, microRNA expression, histone modifications, or a variety of other mechanisms [50]. These changes in gene expression patterns might be related to environmental exposure over time, since increasing genetic divergence in monozygotic twins has been described with advancing age [4].
MicroRNAs (miRNAs) are small noncoding RNA molecules that act as a regulator of gene transcription. Recent studies have demonstrated the role of miRNAs in the therapeutic response of patients with asthma which reinforces the potential for their future use as a predictor biomarker of treatment outcomes. Li et al. [51] identified seven miRNAs (hsa-miR-155-5p, hsa-miR4433b-5p, hsa-miR-532-5p, hsa-miR-345-5p, hsa-miR-652-3p, hsa-miR-126-3p, and hsa-miR-335-5p) associated with lung function changes over 4 years on ICS (budesonide) treatment. Improvement of lung function linked to a better ICS response was significantly associated with hsa-miR-155-5p, hsa-miR4433b-5p, hsa-miR-126-3p, and hsa-miR-335-5p miRNAs, while hsa-miR-532-5p, hsa-miR-345-5p, and hsa-miR-652-3p were linked to relative decrease in FEV1%. Additionally, the hsa-miR-155-5p and hsa-miR-532-5p were functionally validated on dexamethasone-related NF-κB transrepression, and together were predictive of ICS response [41]. Another microRNA with potential for therapeutic prediction is the miR-16, which inhibits the expression of ADRB2, has been negatively correlated with FEV
MicroRNAs (miRNAs) are small noncoding RNA molecules that act as a regulator of gene transcription. Recent studies have demonstrated the role of miRNAs in the therapeutic response of patients with asthma which reinforces the potential for their future use as a predictor biomarker of treatment outcomes. Li et al. [51] identified seven miRNAs (hsa-miR-155-5p, hsa-miR4433b-5p, hsa-miR-532-5p, hsa-miR-345-5p, hsa-miR-652-3p, hsa-miR-126-3p, and hsa-miR-335-5p) associated with lung function changes over 4 years on ICS (budesonide) treatment. Improvement of lung function linked to a better ICS response was significantly associated with hsa-miR-155-5p, hsa-miR4433b-5p, hsa-miR-126-3p, and hsa-miR-335-5p miRNAs, while hsa-miR-532-5p, hsa-miR-345-5p, and hsa-miR-652-3p were linked to relative decrease in FEV1%. Additionally, the hsa-miR-155-5p and hsa-miR-532-5p were functionally validated on dexamethasone-related NF-κB transrepression, and together were predictive of ICS response [47]. Another microRNA with potential for therapeutic prediction is the miR-16, which inhibits the expression of ADRB2, has been negatively correlated with FEV
1
and may function as a biomarker of the response to salmeterol therapy [52].
The DNA methylation is an epigenetic process related to the addition of a methyl group to a cytosine residue in a CpG site and has also been shown to interfere with the therapeutic outcome of asthma. Wang and colleagues [53] have evaluated the association between blood DNA methylation and ICS response in children with persistent asthma measured by the percentage change in FEV
1
eight weeks after treatment initiation. In those individuals, hypermethylation of cg27254601 was associated with both a response to ICS and BOLA2 expression. Additionally, 551 CpG sites were also differentially methylated. The blood ADRB2 5′ untranslated region (5′-UTR) methylation level was also positive associated with a risk of uncontrolled asthma in children. This methylation level was positively correlated with whole blood aluminum concentration in asthmatic children, highlighting the potential that environmental exposures have to modify epigenetic factors [54].
Together, these findings indicate that the lack of response to asthma treatment is a consequence not only of genetic predisposition, but also of environmental/epigenetic factors. Future studies should be carried out to investigate whether environmental exposures are efficient to reprogram genetic predisposition through changes in epigenetic patterns.