Next Article in Journal
Association of HLA-B*51:01, HLA-B*55:01, CYP2C9*3, and Phenytoin-Induced Cutaneous Adverse Drug Reactions in the South Indian Tamil Population
Next Article in Special Issue
α4β1 and αMβ2 Integrin Expression and Pro-Proliferative Properties of Eosinophil Subtypes in Asthma
Previous Article in Journal
Of Screening, Stratification, and Scores
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JPM
Volume 11
Issue 8
10.3390/jpm11080733
jpm-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of ROBO2 as a Potential Locus Associated with Inhaled Corticosteroid Response in Childhood Asthma

by
Natalia Hernandez-Pacheco
1,2,3,*,†,
Mario Gorenjak
4,
Jiang Li
5,
Katja Repnik
4,6,
Susanne J. Vijverberg
7,8,9,
Vojko Berce
4,10,
Andrea Jorgensen
11,
Leila Karimi
12,
Maximilian Schieck
13,14,
Lesly-Anne Samedy-Bates
15,16,
Roger Tavendale
17,
Jesús Villar
3,18,19,
Somnath Mukhopadhyay
17,20,
Munir Pirmohamed
21,
Katia M. C. Verhamme
12,
Michael Kabesch
13,
Daniel B. Hawcutt
22,23,
Steve Turner
24,
Colin N. Palmer
17,
Kelan G. Tantisira
5,25,
Esteban G. Burchard
15,16,
Anke H. Maitland-van der Zee
7,8,9,
Carlos Flores
1,3,26,27,
Uroš Potočnik
4,6,*,‡ and
Maria Pino-Yanes
2,3,27,‡add Show full author list remove Hide full author list
1
Research Unit, Hospital Universitario N.S. de Candelaria, Universidad de La Laguna, Carretera General del Rosario 145, 38010 Santa Cruz de Tenerife, Spain
2
Genomics and Health Group, Department of Biochemistry, Microbiology, Cell Biology and Genetics, Universidad de La Laguna, Avenida Astrofísico Francisco Sánchez s/n, Faculty of Science, Apartado 456, 38200 San Cristóbal de La Laguna, Spain
3
CIBER de Enfermedades Respiratorias, Instituto de Salud Carlos III, Avenida de Monforte de Lemos, 5, 28029 Madrid, Spain
4
Center for Human Molecular Genetics and Pharmacogenomics, Faculty of Medicine, University of Maribor, Taborska Ulica 8, 2000 Maribor, Slovenia
5
The Channing Division of Network Medicine, Department of Medicine, Brigham & Women’s Hospital and Harvard Medical School, 75 Francis St, Boston, MA 02115, USA
6
Laboratory for Biochemistry, Molecular Biology, and Genomics, Faculty of Chemistry and Chemical Engineering, University of Maribor, Smetanova ulica 17, 2000 Maribor, Slovenia
7
Department of Respiratory Medicine, Amsterdam UMC, University of Amsterdam, Meibergdreef 9, 1105 AZ Amsterdam, The Netherlands
8
Division of Pharmacoepidemiology and Clinical Pharmacology, Faculty of Science, Utrecht University, Princetonplein 5, 3584 CC Utrecht, The Netherlands
9
Department of Pediatric Respiratory Medicine and Allergy, Emma’s Children Hospital, Amsterdam UMC, University of Amsterdam, Meibergdreef 9, 1105 AZ Amsterdam, The Netherlands
10
Department of Pediatrics, University Medical Centre Maribor, Ljubljanska Ulica 5, 2000 Maribor, Slovenia
11
Department of Biostatistics, University of Liverpool, Crown Street, Liverpool L69 3BX, UK
12
Department of Medical Informatics, Erasmus University Medical Center, Dr. Molewaterplein 40, 3015 GD Rotterdam, The Netherlands
13
Department of Pediatric Pneumology and Allergy, University Children’s Hospital Regensburg (KUNO), Franz-Josef-Strauß-Allee 11, 93053 Regensburg, Germany
14
Department of Human Genetics, Hannover Medical School, Carl-Neuberg-Straße 1, 30625 Hannover, Germany
15
Department of Medicine, University of California, San Francisco, CA 94143, USA
16
Department of Bioengineering and Therapeutic Sciences, University of California, 533 Parnassus Ave, San Francisco, CA 94143, USA
17
Population Pharmacogenetics Group, Biomedical Research Institute, Ninewells Hospital, and Medical School, University of Dundee, Dundee DD1 9SY, UK
18
Multidisciplinary Organ Dysfunction Evaluation Research Network, Research Unit, Hospital Universitario Dr. Negrín, Calle Barranco de la Ballena s/n, 35019 Las Palmas de Gran Canaria, Spain
19
Keenan Research Center for Biomedical Science, Li Ka Shing Knowledge Institute, St Michael’s Hospital, 30 Bond St, Toronto, ON M5B 1W8, Canada
20
Academic Department of Paediatrics, Brighton and Sussex Medical School, Royal Alexandra Children’s Hospital, 94 N-S Rd, Falmer, Brighton BN2 5BE, UK
21
Department of Molecular and Clinical Pharmacology, Institute of Translational Medicine, University of Liverpool, 200 London Rd, Liverpool L3 9TA, UK
22
Department of Women’s and Children’s Health, University of Liverpool, Liverpool L69 3BX, UK
23
Alder Hey Children’s Hospital, E Prescot Rd, Liverpool L14 5AB, UK
24
Child Health, University of Aberdeen, King’s College, Aberdeen AB24 3FX, UK
25
Division of Pulmonary and Critical Care Medicine, Department of Medicine, Brigham and Women’s Hospital, and Harvard Medical School, 75 Francis St, Boston, MA 02115, USA
26
Genomics Division, Instituto Tecnológico y de Energías Renovables (ITER), Polígono Industrial de Granadilla, 38600 Granadilla, Spain
27
Instituto de Tecnologías Biomédicas (ITB), Universidad de La Laguna, Faculty of Health Sciences, Apartado 456, 38200 San Cristóbal de La Laguna, Spain
 
add Show full affiliation list
 
remove Hide full affiliation list
*
Authors to whom correspondence should be addressed.
†
Current affiliation: Department of Clinical Science and Education, Södersjukhuset, Karolinska Institutet, Sjukhusbacken 10, 118 83 Stockholm, Sweden.
‡
Authors contributed equally to this work.
J. Pers. Med. 2021, 11(8), 733; https://doi.org/10.3390/jpm11080733
Submission received: 18 July 2021 / Accepted: 26 July 2021 / Published: 28 July 2021
(This article belongs to the Special Issue Personalized Therapy of Asthma)
Browse Figures
Versions Notes

Abstract

:
Inhaled corticosteroids (ICS) are the most common asthma controller medication. An important contribution of genetic factors in ICS response has been evidenced. Here, we aimed to identify novel genetic markers involved in ICS response in asthma. A genome-wide association study (GWAS) of the change in lung function after 6 weeks of ICS treatment was performed in 166 asthma patients from the SLOVENIA study. Patients with an improvement in lung function ≥8% were considered as ICS responders. Suggestively associated variants (p-value ≤ 5 × 10−6) were evaluated in an independent study (n = 175). Validation of the association with asthma exacerbations despite ICS use was attempted in European (n = 2681) and admixed (n = 1347) populations. Variants previously associated with ICS response were also assessed for replication. As a result, the SNP rs1166980 from the ROBO2 gene was suggestively associated with the change in lung function (OR for G allele: 7.01, 95% CI: 3.29–14.93, p = 4.61 × 10−7), although this was not validated in CAMP. ROBO2 showed gene-level evidence of replication with asthma exacerbations despite ICS use in Europeans (minimum p-value = 1.44 × 10−5), but not in admixed individuals. The association of PDE10A-T with ICS response described by a previous study was validated. This study suggests that ROBO2 could be a potential novel locus for ICS response in Europeans.

Graphical Abstract

1. Introduction

Asthma is the most common chronic disease in childhood and causes a high burden on the quality of life of the patients and their families. In addition, asthma has an impact in economic terms on the healthcare system, school, and/or work absenteeism [1,2]. This is a complex respiratory disorder characterized by inflammation and reversible obstruction of airways [3], and with diverse manifestations of symptoms, such as wheeze, breathlessness, chest tightness, and cough [1].
Inhaled corticosteroids (ICS) are the most effective and widely prescribed asthma preventive medication [1,2,4,5]. Patients with asthma benefit from ICS therapy through decreased airway inflammation, improved pulmonary capacity, and reduced asthma-related symptoms and exacerbations [6]. Although ICS has demonstrated efficacy in improving symptoms in most children with asthma, between 30% and 40% do not have a complete response to ICS treatment. Furthermore, 10–15% of the children treated with ICS may experience asthma exacerbations worsening or even suffer severe adverse effects [4,7]. Furthermore, not only interindividual differences in ICS response have been described, but also differences among different populations and ethnic groups [4,8]. These differences have been evidenced to be the result of the interaction of several factors, including comorbidities, environmental exposures, and the clinical characteristics of the disease, among others [9]. For instance, the type of asthma (e.g., atopic or non-atopic) could partly contribute to the responsiveness to ICS therapy in some patients. Therefore, subjects with atopic asthma and high levels of blood eosinophils could experience a greater benefit from ICS in contrast to those with a neutrophilic pattern [10,11]. Nonetheless, the important contribution of the individual’s genetic composition has also been suggested [12,13].
Different clinical markers that have been commonly used to evaluate ICS response include the asthma control test [14], asthma symptoms scores [15,16], information about exacerbations [17,18,19], and change in lung function after therapy [1,20]. Among these, performing serial measurements of lung function after a short period of therapy is the most commonly used marker for assessing treatment response [6,18]. The difference between forced expiratory volume in one second (FEV1) values measured at the beginning of treatment and a few weeks [21,22] or months [23] later provides substantial information about ICS response [23]. Importantly, the change in FEV1 after 6 weeks of treatment with ICS has been proposed to be a good predictor of long-term asthma control [21,22]. Although some limitations have been attributed to the evaluation of the lung function (e.g., measurement variability during the day, experience and potential errors driven by the operator, type and calibration of the equipment, and the interpretative algorithm), this approach provides a quantitative and objective measure of the response to asthma treatment [24,25]. Some authors have suggested that variability in ICS response may be explained by the interaction of several factors, including the individual’s genetic composition by means of heritability estimates [9,26]. It has been suggested that approximately 60–80% of the total variation in asthma treatment response might be explained by genetic factors [17,27,28].
Pharmacogenetic studies of ICS have been recently carried out mostly using the genome-wide association study (GWAS) approach [27]. To date, a total of twelve published GWASs have explored the association with ICS response mostly in European populations [28,29,30,31,32,33,34,35,36,37,38,39]. These have identified the association of 28 genetic variants located within or near 17 genes with different measurements of ICS response, being the most common definition of the change in FEV1 after a short period of treatment with ICS. Nonetheless, the validation of some of these associations has suggested that the assessment of the history of recent asthma exacerbations despite ICS treatment can also be used as a proxy of asthma treatment response in different populations [38,39]. Despite the effort of these studies, the genes identified do not explain the response to ICS treatment. Thus, these have not yet provided real improvements in the clinical strategies of asthma management [40], and further genetic variants are expected to be involved in ICS responsiveness [12].
Here, we conducted a pilot study to identify novel genetic variants associated with the response to ICS treatment by means of a GWAS of the change in FEV1 after initiating ICS therapy in asthma patients of European descent. Association with asthma exacerbations of the markers identified was attempted in children and youths treated with ICS from different populations.

2. Results

2.1. Characteristics of the Study Populations

One hundred sixty-six children and young adult asthma patients from the SLOVENIA study [38,39,41] with reported use of any ICS in the last 12 months were included in the discovery phase (Table 1). Of these, 94 were ICS non-responders (cases) and 72 were responders (controls), based on a threshold of 8% FEV1 improvement after 6 weeks of ICS therapy. The individuals included were 10.9 ± 3.4 years old on average, showing a similar mean age in both groups (cases: 10.7 ± 3.2 years, controls: 11.2 ± 3.5 years). ICS responders showed a substantial improvement in pulmonary capacity after 6 weeks of treatment with ICS (16.9% ± 8.7%).
Patients from the European and non-European studies included in the replication of results assessing the association with asthma exacerbations despite ICS treatment showed a similar mean age to those in the SLOVENIA study. However, followMAGICS included older participants (17.2 ± 3.0 years). Since asthma exacerbations were differentially defined among studies, there was variation in the exacerbation rates, ranging from 11.0% in PACMAN to 66.4% in GALA II. Further details about the clinical and demographic characteristics of the populations included in the replication with asthma exacerbations despite ICS use can be found in previous publications [38,39,42].

2.2. Association Results of the Change in FEV1 after ICS Treatment

A total of 7.5 million common single-nucleotide polymorphisms (SNPs) with minor allele frequency (MAF) ≥ 1% with an imputation quality (Rsq) ≥ 0.3 were tested for association with the binary outcome related to the change in FEV1 after ICS treatment in asthma patients from the SLOVENIA study. No evidence of genomic inflation due to population stratification effects was revealed by the value of λGC = 1.00 (Figure S1). No associations were found at genome-wide significance level (p-value ≤ 5 × 10−8), but the SNP rs1166980 located in the ROBO2 gene was found to be suggestively (p-value ≤ 5 × 10−6) associated with ICS responsiveness in asthma patients (odds ratio (OR) for G allele: 7.01, 95% confidence interval (CI): 3.29–14.93, p = 4.61 × 10−7) (Figure 1 and Figure 2).
The association of rs1166980 was also found when assessing the quantitative measurement of the change in FEV1 in the same SLOVENIA participants, since the risk allele for non-response to ICS was also associated with lower lung function improvement (β for G allele: −6.54, 95% CI: −9.74–−3.34, p = 9.41 × 10−5). However, no evidence of association was found for the SNP rs1166980 when the binary (OR for G allele: 1.23, 95% CI: 0.69–1.77, p = 0.453) or the quantitative variables (β for G allele: −0.03, 95% CI: −0.07–0.01, p = 0.144) of the change in FEV1 after 2 months of ICS use were tested for association in the independent study Childhood Asthma Management Program (CAMP) (Table S1).

2.3. Validation of the Association with Asthma Exacerbations despite ICS Use

The association of the SNP rs1166980 was not replicated in Europeans, Latinos/Hispanics, and African Americans from ten independent studies participating in the Pharmacogenomics in Childhood Asthma (PiCA) Consortium [43] when assessing asthma exacerbations despite ICS use as an outcome. At genomic region level, a total of 5919 variants within 100 kilobases (kb) upstream and downstream from ROBO2 were assessed in Europeans. From these, eleven SNPs were significantly associated with asthma exacerbations despite ICS use after accounting for the 164 independent variants located within this region (Bonferroni-like correction significance threshold of p ≤ 3.04 × 10−4). The SNP rs72891545 was the most significant association with ICS response using asthma exacerbations as the outcome (OR for A allele: 4.79, 95% CI: 2.36–9.73, p = 1.44 × 10−5) (Table S2, Figure S2). A total of 6453 variants within a +/−100 kb window from ROBO2 were evaluated in admixed populations. However, no significant associations with asthma exacerbations despite ICS were found after applying a Bonferroni-like correction (p ≤ 1.22 × 10−4 for 411 independent variants).

2.4. Sensitivity Analyses Accounting for Asthma Severity

The association of the SNP rs1166980 with a binary variable of the change in FEV1 after ICS treatment remained statistically significant after also including basal FEV1 measurements as a covariate in the regression models (OR for G allele: 7.21, 95% CI: 3.15–16.50, p = 2.95 × 10−6). Similar results were found when a quantitative outcome of the change in FEV1 was tested in association (β for G allele: −5.58, 95% CI: −8.72–−2.44, p = 6.42 × 10−4).
Further sensitivity analyses were performed in 2282 individuals from six of the eight studies from PiCA from European descent populations with classification of asthma severity based on treatment steps [44]. Specifically, asthma severity was included as a covariate in the regression models for the variant rs72891545. No major differences were found (OR for A allele: 2.66, 95% CI: 1.58–4.45, p = 1.71 × 10−3) compared to the original association models performed for the subset of patients with available information about treatment steps (OR for A allele: 3.66, 95% CI: 2.08–6.43, p = 1.32 × 10−4).

2.5. In Silico Evaluation of the Variants Associated with Different Definitions of ICS Response

According to publicly available functional evidence from the Encyclopedia of DNA Elements (ENCODE) [45], the ROBO2 intronic variant, rs1166980, associated with the change in FEV1 in SLOVENIA, could be involved in the regulation of the gene expression, since this nucleotide change is related to the modification of the binding site for different transcription factors (e.g., CTCF, Sox, and p300) [45]. Moreover, rs72891545, the most significant association signal with asthma exacerbations, has also been predicted to involve the modification of several transcription factor binding sites (ELF1, Myc, Sp4, YY1, Zfx). Although no previous evidence of significant expression quantitative trait loci (eQTL) or methylation quantitative trait loci (meQTL) was found for any of these variants [45], both have been significantly associated as protein quantitative trait loci (pQTL) (p ≤ 0.01) [46,47]. Specifically, rs1166980 has been associated with the expression of 38 different proteins in plasma (Table S3) [48,49,50], including some involved in the immune response (e.g., PI3, COLEC11, and TNFRSF21), inflammatory processes (HPGD), and drug metabolism (TBXAS1). Moreover, it has been evidenced to be associated with protein expression levels of PMEPA1, a negative regulator of the activity of transforming growth factor β (TGF-β), which has been widely proposed to play a key role in allergy and asthma pathophysiology [51]. On the other hand, rs72891545 had been previously associated with the expression levels of 19 proteins [49] implicated in different functions, including some related to asthma pathophysiology, such as the maintenance of the innate immunity (ERAP2, FCN1, and SIGLEC9), collagen deposition (COL8A1), and the TGF-β signaling pathway (MSTN) (Table S4). Additionally, several proteins, whose expression could be regulated by these variants, had been previously associated with asthma-related traits and/or allergic diseases (e.g., OLFML3, PCSK3, ZNF180, GALNT16, and KLK14) (Tables S3 and S4) [52].

2.6. Validation of Previous Associations with ICS Response

Among the 26 SNPs previously associated with ICS response through GWAS approaches in studies independent from the PiCA consortium (Table S5), two variants were found to be nominally associated with the binary outcome related to the change in FEV1 after 6 weeks of ICS treatment in the SLOVENIA study: rs2395672 at CMTR1 (OR for G allele: 1.78, 95% CI: 1.03–3.05, p = 0.037) and rs3827907 at EDDM3B (OR for C allele: 0.52, 95% CI: 0.32–0.84, p = 7.40 × 10−3) (Table S6). However, these did not remain significant after adjusting for the total number of variants assessed (p ≤ 1.92 × 10−3). At the genomic-region level, a total of 33,617 variants located within a 100 kb window from genes previously associated with ICS response were assessed. This resulted in evidence of suggestive replication for two variants located in the intergenic region of PDE10A and T after Bonferroni-like correction of the significance threshold within each genomic region: rs9365939 (OR for G allele: 0.41, 95% CI: 0.26–0.65, p = 1.92 × 10−4) and rs2118353 (OR for T allele: 0.41, 95% CI: 0.26–0.65, p = 1.92 × 10−4) (Table S7). However, these associations would not be considered significant after correcting for the total number of independent SNPs tested across all the genomic regions (p ≤ 2.89 × 10−5 for 1728 independent variants).

3. Discussion

This is a pilot study describing the results of one of the few GWASs of the response to ICS in asthma carried out in ICS-naïve patients to date. A variant located in the ROBO2 gene was the most significantly associated with a binary variable of ICS responsiveness based on FEV1 change after 6 weeks of ICS treatment. Consistent with this finding, the same result was also found with the quantitative variable of this change, although no evidence of replication with the same outcome was found in an independent study considering a wider timeframe. The association of ROBO2 was validated at genomic-region level, by analyzing asthma exacerbations despite ICS use in Europeans. However, no evidence of replication was found in Latinos/Hispanics and African Americans. Moreover, the effects of the variants at ROBO2 might be driven by the response to asthma treatment rather than disease severity since significance remained after including covariates as proxies of asthma severity.
ROBO2 encodes one of the members of the roundabout guidance receptor’s family, which are immunoglobulins highly conserved across species. Four ROBO proteins have been identified in humans [53]. These are transmembrane receptors binding Slit guidance ligands (SLIT) [53,54,55] with functions initially linked to the development of the nervous system, including the modulation of axon guidance and cell migration [53,54,55]. However, they have also been demonstrated to be involved in different processes, including some that are part of the asthma pathophysiology. Specifically, ROBO/SLIT has been related to cell adhesion, migration, growth, and survival [56], in addition to the morphogenesis of normal and aberrant pulmonary tissues [56]. There is also some evidence suggesting the implication of the ROBO signaling pathway in the regulation of innate immunity [57,58]. Previous studies have related the increased production of chemoattractants with the inhibition of the ROBO expression, promoting the migration of immune cells from the blood circulatory system to the pulmonary tissue and airway remodeling processes, decreasing the number of alveoli [59]. Specifically, ROBO2 is involved in the signal transduction of SLIT2 [55], which has been shown to have an important function in pulmonary diseases [55,58,60]. SLIT2 is implicated in the regulation of chemotaxis and migration of several types of immune cells into the lung (e.g., leukocytes, T lymphocytes, dendritic cells, macrophages, and neutrophils), preventing inflammatory processes [57,61,62]. Nonetheless, members of ROBO/SLIT have also been found to be involved in the eosinophil chemotaxis induced by eotaxin, enhancing allergic airway inflammation [63]. Although little is known about its role in allergic processes outside the respiratory system, these pieces of evidence suggest its potential implication in different types of allergic responses. In addition, SLIT2 has been found to inhibit the migration of monocytes from the systemic circulation and their differentiation into fibrocytes in the lung, proliferation of fibroblasts, and collagen production through TGF-β activity [64]. SLIT2 and ROBO2 have been evidenced to prevent fibrotic processes in several diseases, including pulmonary fibrosis [64]. Moreover, Lin et al. detected decreased levels of ROBO2 and SLIT2 in chronic obstructive pulmonary disease (COPD) patients [58], a disease with underlying mechanisms shared with asthma [65,66].
These findings suggest that SLIT ligands could be involved in anti-inflammatory and antifibrotic processes. Thus, ROBO proteins could also be involved in processes occurring later in life with important implications in pulmonary disorders [57,64]. Interestingly, Ning et al. found that ROBO2/SLIT2 could be involved in the migration of rat airway smooth muscle (ASM) cells within the airway wall [67], which has been proposed to be a key feature of the structural changes taking place in asthma pathophysiology [68]. The migration of ASM induced by platelet-derived growth factor was detected to be inhibited by the administration of SLIT2 [67]. Furthermore, ROBO2 has been evidenced to be an important factor in triggering airways constriction in asthma and COPD [69]. The expression levels of ROBO2 and SLIT2 have been negatively correlated with COPD progression [58]. Specifically, it has been hypothesized that the downregulation of both genes could activate Cdc42 and Rac2 GTPases, promoting the migration of neutrophils and T lymphocytes into the lung, triggering inflammation [58], a mechanism that could also take place in asthma.
Additionally, ROBO2 has been suggestively associated with post-bronchodilator spirometric measures in African Americans. Specifically, three intronic ROBO2 variants were associated with FEV1 and the ratio between FEV1 and the forced vital capacity (FVC) measured after the administration of short-acting β2 agonists (SABA) [70]. Furthermore, this gene was identified to be a shared genetic factor for asthma susceptibility among European, African American, and Latino/Hispanic populations [71]. Ding et al. also suggested that ROBO2 may be part of biological networks related to inflammatory diseases and disorders of the immune system [71]. Altogether, this evidence suggests that ROBO2 could play an important role in asthma phenotypes, including the response to ICS in asthma.
The findings of this study suggest the potential implication of ROBO2 in the response to asthma treatment with ICS specifically in European populations, with a lack of evidence of replication in Latinos/Hispanics and African Americans. This could be explained by the fact that Latinos/Hispanics and African Americans widely differ from homogenous populations in terms of linkage disequilibrium patterns, allele frequency, gene–gene, and gene–environment interactions as a result of the recent admixture processes as the origin of these populations [72]. In fact, previous studies have suggested the potential existence of ancestry-specific genetic variation in the response to this medication [38,39].
As part of this study, we also assessed the replication of SNPs and genes that have previously been implicated in ICS response, providing evidence of an association of the intergenic region of the PDE10A and T genes with FEV1 change after ICS treatment. However, this association was demonstrated for different SNPs from those described in the study reporting the association of this genomic region with ICS response [28]. Interestingly, evidence of replication had also been found for this locus with ICS response measured as the occurrence or absence of asthma exacerbations in European populations [39]. These findings suggest that PDE10A-T could play an important role in ICS response in asthma patients, at least those of European descent. Furthermore, the association of this locus with ICS response was first identified in children and adults, which has been now replicated in children and young adults, suggesting the potential existence of some genetic factors shared between childhood and adulthood asthma. However, further investigation is needed to be able to better understand the genetic contribution to asthma treatment response in different age groups.
The preliminary results of this study demonstrate the importance of omics approaches to provide insights into asthma-related traits. Nonetheless, further investigation of the genetic factors underlying ICS response in combination with the assessment of the molecular modification of the individual’s genetic material in response to environmental exposures is crucial to better understand the mechanisms underlying ICS responsiveness in asthma patients [73,74].
This study has several limitations that need to be acknowledged. First, the sample size of the discovery phase was limited, which could cause only one variant to be suggestively associated with ICS responsiveness, with a lack of genome-wide significant associations. Nonetheless, the fact that FEV1 was measured in ICS-naïve patients with asthma, an approach that is rare among studies, explains part of the difficulty in achieving a larger sample size. Second, information related to spirometry recordings before and after a short period of ICS treatment was available only in one independent study, and no evidence of replication using the same measurement was found, which could be explained by the reduced number of individuals included in the analyses and the different timeframe considered. However, the association of ROBO2 with ICS response among asthma patients was also found when evaluating the association with asthma exacerbations despite ICS treatment. Third, different definitions of asthma exacerbations were used based on retrospective information from European and admixed asthmatic children treated with ICS, which could not be fully informative about the response to asthma treatment. Fourth, information about the specific type of ICS used and the doses administered, or indices of treatment adherence were not available for any of the studies evaluated. Fifth, covariates related to the type of asthma, exposure to environmental triggers, or the presence of comorbidities were not included in the analyses given the lack of such information.
In conclusion, this study suggests the association with a variant in ROBO2 and the change in FEV1 after ICS treatment in European children and young adults with asthma. This association was validated using asthma exacerbations despite ICS use as an alternative outcome in independent European and admixed populations. Taken together with the biologically plausible role of ROBO2 in pulmonary and immune functions, ROBO2 potentially represents a novel locus determining the response to ICS in patients with asthma. Larger studies coupled with functional evaluation are required to fully understand the role of ROBO2 in responsiveness to ICS in patients with asthma.

4. Materials and Methods

4.1. Study Population Analyzed in the Discovery Phase

Patients with mild-to-moderate persistent asthma from the SLOVENIA study were included in the GWAS of change in lung function after ICS treatment. Children and young adults (5–18 years old) of Slovenian origin were included in this study. Co-existence of other chronic inflammatory disorders, except for asthma and atopic diseases, was considered as an exclusion criterion [38,41].
A subset of patients with reports of at least one use of any type of ICS and/or combination with long-acting β2 agonists (LABA) in the 12 months preceding the study enrolment was analyzed. Availability of genome-wide genotypes, data on the change in FEV1 after 6 weeks of ICS therapy, and information regarding asthma exacerbations were considered as inclusion criteria in the GWAS analyses. FEV1, expressed as the percentage of the predicted value based on sex, age, and height of the patients, was measured before the beginning of ICS treatment (when the patients were ICS-naïve) and 6 weeks after the start of the treatment using a Vitalograph 2150 spirometer (Compact, Buckingham, UK), according to the standard guidelines [41,75]. On the basis of these measurements, the percent change in FEV1 (ΔFEV1) was calculated as (post-FEV1–pre-FEV1)/(pre-FEV1) × 100. Based on a threshold of 8% FEV1 improvement, which has been shown to be a good predictor of asthma treatment response in children [38,76], participants were classified as ICS responders (ΔFEV1 ≥ 8%) or non-responders (ΔFEV1 < 8%).

4.2. Genotyping and Imputation of Genetic Variants in SLOVENIA

The SLOVENIA samples were genotyped using the Illumina Global Screening Array-24 v1.0 BeadChip (Illumina Inc., San Diego, CA, USA). Quality control (QC) analyses were carried out with PLINK 1.9 [77,78] and genetic variants across the whole genome were imputed following the procedures described elsewhere [38,39].

4.3. Association Testing with the Change in FEV1 Defined as a Binary Variable

The association of genetic variants with the binary variable of ICS response was tested using logistic regression models with the binary Wald test implemented in EPACTS 3.2.6 [79], including age and sex as covariates. Association analyses were also adjusted by the first two principal components (PCs) of genetic ancestry estimated by means of EIGENSOFT [80]. This model was selected since it showed the best fit with the expected values in terms of significance, assuming no association as null hypothesis attending to λGC values, estimated through the R package gap [81], and quantile–quantile plots.
Results were filtered to retain common SNPs (MAF ≥ 1%), and imputation quality Rsq ≥ 0.3 and variants that reached a significance threshold of p-value ≤ 5 × 10−6 were deemed suggestively associated and followed up for replication in additional studies. This threshold was set arbitrarily, following what has been commonly adopted by previous GWAS [38,39].

4.4. Association with the Quantitative Change in FEV1 after ICS Treatment

SNPs suggestively associated with the binary outcome related to the change in FEV1 after 6 weeks of ICS treatment were assessed in the same group of asthma patients from the SLOVENIA study, but evaluating the association with the quantitative form of this outcome. Linear regression models were performed through linear Wald tests in EPACTS 3.2.6 [79], adjusted by the same aforementioned covariates.
Validation was also attempted with the change in FEV1 after a short period under ICS therapy in an independent study of European ancestry. Asthma patients with reported use of ICS in the previous 12 months from the CAMP study were included in the analyses (Table S1). ICS responders and non-responders were classified based on the change in FEV1 after 2 months of ICS use using the same criteria considered in the discovery phase. Genetic variants were imputed using phase 3 of the 1000 Genomes Project (1KGP) [82] through the Michigan Imputation Server [83]. Further information is described in the Supplementary Material. Binary and quantitative variables of the ICS response, measured as the change in FEV1, were tested in association through logistic and linear regressions using PLINK 1.9 [77,78], respectively. Validation of the association was considered for nominally significant variants (p-value ≤ 0.05) with the same direction of the association effect as in the discovery phase.

4.5. Replication of Results Analyzing the Association with Asthma Exacerbations despite ICS Use in Additional Studies

The genetic markers found to be associated with the change in FEV1 after ICS treatment were attempted for validation with the absence or presence of asthma exacerbations despite the use of ICS. This was done in ten independent studies included in the PiCA consortium [43]. Association testing was undertaken in asthma patients (2–25 years old) treated with ICS in the previous year, separately performed in two groups of studies on the basis of their ancestry.
On one hand, eight independent European studies were analyzed: BREATHE; Effectiveness and Safety of Treatment with Asthma Therapy in Children (ESTATe); the follow-up stage of the Multicentre Asthma Genetics in Childhood Study (followMAGICS); Genetics of the Scottish Health Research Register (GoSHARE); the Pharmacogenetics of Asthma medication in Children: Medication with Anti-Inflammatory effects (PACMAN); the Pediatric Asthma Gene Environment Study (PAGES); the Pharmacogenetics of Adrenal Suppression with Inhaled Steroids (PASS) study; SLOVENIA. Additionally, two recently admixed populations with African ancestry were also included in association analyses: Latinos/Hispanics and African Americans from the Genes–Environment and Admixture in Latino Americans (GALA II) study, and African Americans included in the Study of African Americans, Asthma, Genes, and Environments (SAGE), respectively.
Severe asthma exacerbations were defined as the need for emergency care, hospitalizations, or systemically administered corticosteroids because of asthma in the previous 6 or 12 months depending on the study. Alternatively, moderated exacerbations evidenced by unscheduled general practitioner or respiratory system specialist visits and school absences were considered for BREATHE–PAGES, BREATHE, and followMAGICS given the lack of information. ICS use was defined using the same criteria described for participants in the SLOVENIA study. Further description of the characteristics of the study populations, genotyping, imputation, and association analyses are available in the Supplementary Material and elsewhere [38,39,42].
Association with the presence/absence of asthma exacerbations despite ICS use was assessed in each study using the same methodology described in previous publications [38,39,42]. The association of the SNPs identified in the discovery phase was separately evaluated in each ancestry group of studies. Replication was carried out at the SNP level, but also genomic regions were considered, including variants located within a 100 kb window upstream and downstream from the limits of the genes where the variants were located. Only common SNPs with MAF ≥ 1% and Rsq ≥ 0.3 shared among the populations included in each group were included. Replication results were considered significant for those SNPs that reached the Bonferroni-corrected significance threshold, estimated as α = 0.05/number of independent signals within each genomic region, an approach considered to provide the closest approximation to permutation-based methods [84]. For this, independent variants were separately estimated for Europeans and non-Europeans through empirical autocorrelations based on the –log10 p-value of each SNP analyzed using the R package coda [77,78,85].

4.6. Sensitivity Analyses of Asthma Treatment Response

Sensitivity analyses were carried out for the variants identified to ascertain whether the association effect detected was driven by disease severity rather than a response to asthma medications. First, regression analyses evaluating the association with the binary variable of the change in FEV1 after ICS therapy were also adjusted by the FEV1 measured at the beginning of ICS therapy. On the other hand, a modified classification into treatment steps [44] as a proxy of asthma severity was included as a covariate. Nonetheless, this could not be carried out in SLOVENIA, since this study had incomplete information about the different medications included in the definition of treatment steps. Therefore, sensitivity analyses accounting for treatment step classification assessing the association with asthma exacerbations despite ICS use. The same methodology previously described was applied [39,42].

4.7. In Silico Functional Evaluation of Variants Associated with ICS Response

The potential functional implications of the variants suggestively associated with the change in FEV1 after 6 weeks of ICS treatment in SLOVENIA and those with evidence of replication with asthma exacerbations in ICS users were assessed using several publicly available databases. The potential role as eQTL, DNase hypersensitivity sites, and histone marks was assessed using HaploReg v4.1 [86], based on the functional evidence from ENCODE [45]. Significant evidence as pQTL or meQTL previously reported by publicly available studies was also explored using the PhenoScanner v2 tool [46,47].

4.8. Validation of Previous Associations with ICS Response

Previous GWASs, apart from those carried out as part of the PiCA consortium, had identified a total of 26 SNPs near or within 15 genes associated with ICS response in several populations (Table S5). Validation of these associations was attempted at SNP level using the results of the GWAS of the binary variable of the change in FEV1 after ICS treatment performed in the SLOVENIA study. Evidence of replication was considered for significant variants at the nominal level (p ≤ 0.05) with the same direction of the association effect as in the discovery phase. Replication was also evaluated at the genomic-region level, including a 100 kb window from the limits of the genes of ICS response previously identified. A Bonferroni-like correction was applied accounting for the number of independent variants analyzed within each genomic region. To avoid being too stringent and conservative, evidence of replication was considered for those association signals reaching the Bonferroni-corrected significance threshold estimated as α = 0.05/number of independent variants [38,39] instead of considering all the variants within the region, which could be correlated by LD.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/jpm11080733/s1, Supplementary methods, Figure S1: Quantile-quantile plot of association results of ICS response measured as the binary outcome related to the change in FEV1 after ICS treatment, Figure S2: Regional plot of association results with asthma exacerbations despite ICS use in European children and young adults, Table S1: Clinical and demographic characteristics of the asthma patients treated with ICS from CAMP included in the evaluation of the association with the change in FEV1 after ICS therapy, Table S2: Genomic-region replication of ROBO2 with asthma exacerbations despite ICS use in European populations. Evidence for significant variants after Bonferroni-like correction, Table S3: Proteins with expression levels in plasma associated with rs1166980, Table S4: Proteins with expression levels affected by the SNP rs72891545, Table S5: Genes associated with ICS response by genome-wide association studies published to date, Table S6: Results of SNP-level replication of previous associations of ICS response in the results of the GWAS of the change in FEV1 after ICS treatment performed in the SLOVENIA study, Table S7: Genomic-region replication of previous associations of ICS response. Evidence of association with the change in FEV1 after ICS treatment in asthma patients.

Author Contributions

Conceptualization, U.P. and M.P.-Y.; data curation, N.H.-P., M.G., J.L., K.R., S.J.V., V.B., A.J., L.K. and L.-A.S.-B.; formal analysis, N.H.-P., M.G., J.L., K.R., L.K., L.-A.S.-B. and S.T.; funding acquisition, S.J.V., M.S., R.T., S.M., M.P., K.M.C.V., M.K., D.B.H., S.T., C.N.P., K.G.T., E.G.B., A.H.M.-v.d.Z., U.P. and M.P.-Y.; investigation, N.H.-P., M.G., K.R., S.J.V., V.B., L.K., M.S., L.-A.S.-B., R.T., J.V., S.M., M.P., K.M.C.V., M.K., D.B.H., S.T., C.N.P., K.G.T., E.G.B., A.H.M.-v.d.Z., C.F., U.P. and M.P.-Y.; methodology, N.H.-P., M.G., K.R., S.J.V., L.K., L.-A.S.-B., J.V., U.P. and M.P.-Y.; supervision, C.F., U.P. and M.P.-Y.; validation, N.H.-P., M.G., J.L., S.J.V., A.J. and L.K.; writing—original draft, N.H.-P., C.F., U.P. and M.P.-Y.; writing—review and editing, N.H.-P., M.G., J.L., K.R., S.J.V., V.B., A.J., L.K., M.S., L.-A.S.-B., R.T., J.V., S.M., M.P., K.M.C.V., M.K., D.B.H., S.T., C.N.P., K.G.T., E.G.B., A.H.M.-v.d.Z., C.F., U.P. and M.P.-Y. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the award AC15/00015 funded by the Instituto de Salud Carlos III (ISCIII) through Strategic Action for Health Research (AES) and European Community (EC) within the Active and Assisted Living (AAL) Programme framework (MP-Y) and the SysPharmPedia grant from the ERACoSysMed 1st Joint Transnational Call from the European Union under the Horizon 2020. This study was also funded by the Spanish Ministry of Science, Innovation, and Universities, the State Research Agency, and the European Regional Development Fund from the European Union (MICIU/AEI/FEDER, UE, SAF2017-83417R). N.H-P was supported by a fellowship (FI16/00136) from Instituto de Salud Carlos III (ISCIII) and co-funded by the European Social Funds from the European Union (ESF) “ESF invests in your future”. MP-Y was funded by the Ramón y Cajal Program (RYC-2015-17205) by the Spanish Ministry of Science, Innovation and Universities. The PACMAN study was funded by a strategic alliance between GlaxoSmithKline and Utrecht Institute for Pharmaceutical Sciences. The SLOVENIA study was financially supported by the Slovenian Research Agency (research core funding No. P3-0067) and from SysPharmPedia grant, co-financed by Ministry of Education, Science and Sport Slovenia (MIZS) (contract number C3330-16-500106). GALA II was supported by the National Heart, Lung, and Blood Institute of the National Institute of Health (NIH) grants R01HL117004 and X01HL134589; study enrolment supported by the Sandler Family Foundation, the American Asthma Foundation, the RWJF Amos Medical Faculty Development Program, Harry Wm. and Diana V. Hind Distinguished in Pharmaceutical Sciences II and the National Institute of Environmental Health Sciences grant R01ES015794. SAGE was funded by the National Heart, Lung, and Blood Institute of the National Institute of Health (NIH) grants R01HL117004 and X01HL134589; study enrolment supported by the Sandler Family Foundation, the American Asthma Foundation, the RWJF Amos Medical Faculty Development Program, Harry Wm. and Diana V. Hind Distinguished in Pharmaceutical Sciences II. L.S.B. was supported by National Institute of General Medical Sciences of the National Institutes of Health [5T32GM007546-42]. The SHARE Bioresource (GoSHARE) and SHARE have ongoing funding from NHS Research Scotland and established by funding from The Wellcome Trust Biomedical Resource [Grant No. 099177/Z/12/Z]. The PASS study was funded by the NHS Chair of Pharmacogenetics via the UK Department of Health. MP is Emeritus NIHR Senior Investigator. ESTATe was funded by an independent research grant by ZonMw project (113201006). Genotyping of samples from BREATHE-PAGES, GoSHARE, and SCSGES was carried out at CeGen-PRB3-ISCIII; supported by ISCIII and European Regional Development Fund (ERDF) (PT17/0019).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the following Institutional Review Boards. The Slovenian National Medical Ethics Committee (Ljubljana, Slovenia) approved SLOVENIA (reference number: 0120-569/2017/4). The clinic’s institutional review board (IRB) approved CAMP (Boston, USA) (ethics approval number: 1999-P-001549). The Medical Ethics Committee of the University Medical Centre Utrecht (Utrecht, the Netherlands) approved PACMAN (protocol number: 08/023). PAGES was approved by the Cornwall and Plymouth Research Ethics Committee (Plymouth, UK). GoSHARE and BREATHE were approved by the Tayside Committee on Medical Research Ethics (Dundee, UK). The Liverpool Paediatric Research Ethics Committee (Liverpool, UK) (reference number: 08/H1002/56) approved PASS. followMAGICS was approved by the Ethik-Kommission der Bayerischen Landesärztekammer (Munich, Germany) (ethics reference number: 01218). The Medische Ethische Toetsings Commissie, Erasmus Medical Centre (Rotterdam, the Netherlands) (ethics approval number: MEC-2011-474) approved ESTATe. The Human Research Protection Program Institutional Review Board of the University of California, San Francisco (San Francisco, USA) approved GALA II and SAGE (ethics approval numbers: 10-00889 and 10-02877, respectively).

Informed Consent Statement

Informed consent for inclusion was obtained from all participants or given by their parents/caregivers before they participated in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors acknowledge the patients, families, recruiters, health care providers and community clinics for their participation in all the studies included in the PiCA consortium (http://pica-consortium.org). The authors thank the contribution of Teide High-Performance Computing facilities (http://teidehpc.iter.es) provided by the Instituto Tecnológico y de Energías Renovables (ITER, S.A.) to the results of this research and also the Centro Nacional de Genotipado-Plataforma de Recursos Biomoleculares-Instituto de Salud Carlos III (CeGen-PRB3-ISCIII; www.cegen.org) for the genotyping services provided. The GALA II and SAGE study collaborators include Shannon Thyne, UCSF; Harold J. Farber, Texas Children’s Hospital; Denise Serebrisky, Jacobi Medical Center; Rajesh Kumar, Lurie Children’s Hospital of Chicago; Emerita Brigino-Buenaventura, Kaiser Permanente; Michael A. LeNoir, Bay Area Pediatrics; Kelley Meade, UCSF Benioff Children’s Hospital, Oakland; William Rodriguez-Cintron, VA Hospital, Puerto Rico; Pedro C. Avila, Northwestern University; Jose R. Rodriguez-Santana, Centro de Neumologia Pediatrica; Luisa N. Borrell, City University of New York; Adam Davis, UCSF Benioff Children’s Hospital, Oakland; Saunak Sen, University of Tennessee and Fred Lurmann, Sonoma Technologies, Inc and Michael A. LeNoir, Bay Area Pediatrics. The authors acknowledge the families and patients for their participation and thank the numerous health care providers and community clinics for their support and participation in GALA II and SAGE. In particular, the authors thank study coordinator Sandra Salazar; the recruiters who obtained the data: Duanny Alva, Gaby Ayala-Rodriguez, Lisa Caine, Elizabeth Castellanos, Jaime Colon, Denise DeJesus, Blanca Lopez, Brenda Lopez, Louis Martos, Vivian Medina, Juana Olivo, Mario Peralta, Esther Pomares, Jihan Quraishi, Johanna Rodriguez, Shahdad Saeedi, Dean Soto, Ana Taveras; and the lab researcher Celeste Eng who processed the biospecimens.

Conflicts of Interest

NH-P declares funding from Instituto de Salud Carlos III (ISCIII) and the European Social Funds. SM reports funding from The Gannochy Trust, Perth and Kinross City Council, and Scottish Enterprises Tayside. MP declares funding from the MRC Clinical Pharmacology Training Scheme, EPSRC, Astra Zeneca, and Bristol Myers Squibb. KMV reports funding from ZonMw, Yamanouchi, Pfizer/Boehringer Ingelheim, Novartis, and GSK. MK declares funding from the European Union, German Ministry of Education, Research, German Research Foundation, and other sources. KGT reports funding from the National Institutes of Health of the USA. EGB reports funding from the National Institute of Health, National Institute of Health and Environmental Health Sciences, National Institute on Minority Health and Health Disparities, National Institute of General Medical Sciences, and National Human Genome Research Institute USA, Sandler Family Foundation, American Asthma Foundation, and Robert Wood Johnson Foundation. A-HM declares funding from GlaxoSmithKline, Boehringer Ingelheim, and Astra Zeneca. UP reports funding from the Slovenian Research Agency and the Ministry of Education, Science, and Sport of Slovenia. JV declares funding from ISCIII and the European Regional Development Fund. MP-Y declares funding from ISCIII, the European Regional Development Fund, and the Spanish Ministry of Science, Innovation, and Universities, and funding from GlaxoSmithKline outside of the submitted work.

References

  1. Global Initiative for Asthma. Global Strategy for Asthma Management and Prevention. 2020. Available online: https://ginasthma.org/ (accessed on 25 October 2020).
  2. Cerasoli, F., Jr. Developing the Ideal Inhaled Corticosteroid. Chest 2006, 130 (Suppl. 1), 54S–64S. [Google Scholar]
  3. Yang, J.; Benyamin, B.; McEvoy, B.P.; Gordon, S.; Henders, A.K.; Nyholt, D.R.; Madden, P.A.; Heath, A.C.; Martin, N.G.; Montgomery, G.W.; et al. Common SNPs Explain a Large Proportion of the Heritability for Human Height. Nat. Genet. 2010, 42, 565–569. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Szefler, S.J.; Phillips, B.R.; Martinez, F.D.; Chinchilli, V.M.; Lemanske, R.F.; Strunk, R.C.; Zeiger, R.S.; Larsen, G.; Spahn, J.D.; Bacharier, L.B.; et al. Characterization of Within-Subject Responses to Fluticasone and Montelukast in Childhood Asthma. J. Allergy Clin. Immunol. 2005, 115, 233–242. [Google Scholar] [PubMed]
  5. Scelfo, C.; Galeone, C.; Bertolini, F.; Caminati, M.; Ruggiero, P.; Facciolongo, N.; Menzella, F. Towards Precision Medicine: The Application of Omics Technologies in Asthma Management. F1000Research 2018, 7, 423. [Google Scholar] [CrossRef] [PubMed]
  6. Dijkstra, A.; Vonk, J.M.; Jongepier, H.; Koppelman, G.H.; Koppelman, J.P.; ten Hacken, N.H.T.; Timens, W.; Postma, D.S. Lung Function Decline in Asthma: Association with Inhaled Corticosteroids, Smoking and Sex. Thorax 2006, 61, 105–110. [Google Scholar]
  7. Drazen, J.M.; Silverman, E.K.; Lee, T.H. Heterogeneity of Therapeutic Responses in Asthma. Br. Med. Bull. 2000, 56, 1054–1070. [Google Scholar]
  8. Mersha, T.B. Mapping Asthma-Associated Variants in Admixed Populations. Front. Genet. 2015, 6, 292. [Google Scholar]
  9. Ramadan, A.A.; Gaffin, J.M.; Israel, E.; Phipatanakul, W. Asthma and Corticosteroid Responses in Childhood and Adult Asthma. Clin. Chest Med. 2019, 40, 163–177. [Google Scholar]
  10. Fitzpatrick, A.M.; Teague, W.G.; Meyers, D.A.; Peters, S.P.; Li, X.; Li, H.; Wenzel, S.E.; Aujla, S.; Castro, M.; Bacharier, L.B.; et al. Heterogeneity of Severe Asthma in Childhood: Confirmation by Cluster Analysis of Children in the National Institutes of Health/National Heart, Lung, and Blood Institute Severe Asthma Research Program. J. Allergy Clin. Immunol. 2011, 127, 382–389.e13. [Google Scholar]
  11. Moore, W.C.; Hastie, A.T.; Li, X.; Li, H.; Busse, W.W.; Jarjour, N.N.; Wenzel, S.E.; Peters, S.P.; Meyers, D.A.; Bleecker, E.R. Sputum Neutrophil Counts Are Associated with More Severe Asthma Phenotypes Using Cluster Analysis. J. Allergy Clin. Immunol. 2014, 133, 1557–1563.e5. [Google Scholar] [CrossRef] [Green Version]
  12. Park, H.W.; Tantisira, K.G.; Weiss, S.T. Pharmacogenomics in Asthma Therapy: Where Are We and Where Do We Go? Annu. Rev. Pharmacol. Toxicol. 2015, 55, 129–147. [Google Scholar] [CrossRef]
  13. Duong-Thi-Ly, H.; Nguyen-Thi-Thu, H.; Nguyen-Hoang, L.; Nguyen-Thi-Bich, H.; Craig, T.J.; Duong-Quy, S. Effects of Genetic Factors to Inhaled Corticosteroid Response in Children with Asthma: A Literature Review. J. Int. Med. Res. 2017, 45, 1818–1830. [Google Scholar] [CrossRef] [Green Version]
  14. Al Moamary, M.S.; Al-Kordi, A.G.; Al Ghobain, M.O.; Tamim, H.M. Utilization and Responsiveness of the Asthma Control Test (ACT) at the Initiation of Therapy for Patients with Asthma: A Randomized Controlled Trial. BMC Pulm. Med. 2012, 12, 14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Gorelick, M.H.; Stevens, M.W.; Schultz, T.R.; Scribano, P.V. Performance of a Novel Clinical Score, the Pediatric Asthma Severity Score (PASS), in the Evaluation of Acute Asthma. Acad. Emerg. Med. 2004, 11, 10–18. [Google Scholar] [CrossRef] [PubMed]
  16. Keogh, K.A.; Macarthur, C.; Parkin, P.C.; Stephens, D.; Arseneault, R.; Tennis, O.; Bacal, L.; Schuh, S. Predictors of Hospitalization in Children with Acute Asthma. J. Pediatr. 2001, 139, 273–277. [Google Scholar] [CrossRef] [PubMed]
  17. Fuhlbrigge, A.; Peden, D.; Apter, A.J.; Boushey, H.A.; Camargo, C.A., Jr.; Gern, J.; Heymann, P.W.; Martinez, F.D.; Mauger, D.; Teague, W.G.; et al. Asthma Outcomes: Exacerbations. J. Allergy Clin. Immunol. 2012, 129 (Suppl. 3), S34–S48. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. National Asthma Education and Prevention Program, Third Expert Panel on the Diagnosis and Management of Asthma. Expert Panel Report 3: Guidelines for the Diagnosis and Management of Asthma. Section 2, Definition, Pathophysiology and Pathogenesis of Asthma, and Natural History of Asthma; National Heart, Lung, and Blood Institute: Bethesda, MD, USA, 2007.
  19. Aldington, S.; Beasley, R. Asthma Exacerbations. 5: Assessment and Management of Severe Asthma in Adults in Hospital. Thorax 2007, 62, 447–458. [Google Scholar] [CrossRef] [Green Version]
  20. Gorelick, M.H.; Stevens, M.W.; Schultz, T.; Scribano, P.V. Difficulty in Obtaining Peak Expiratory Flow Measurements in Children with Acute Asthma. Pediatr. Emerg. Care 2004, 20, 22–26. [Google Scholar] [CrossRef]
  21. Szefler, S.J.; Martin, R.J.; King, T.S.; Boushey, H.A.; Cherniack, R.M.; Chinchilli, V.M.; Craig, T.J.; Dolovich, M.; Drazen, J.M.; Fagan, J.K.; et al. Significant Variability in Response to Inhaled Corticosteroids for Persistent Asthma. J. Allergy Clin. Immunol. 2002, 109, 410–418. [Google Scholar] [CrossRef] [Green Version]
  22. Martin, R.J.; Szefler, S.J.; King, T.S.; Kraft, M.; Boushey, H.A.; Chinchilli, V.M.; Craig, T.J.; Dimango, E.A.; Deykin, A.; Fahy, J.V.; et al. The Predicting Response to Inhaled Corticosteroid Efficacy (PRICE) Trial. J. Allergy Clin. Immunol. 2007, 119, 73–80. [Google Scholar] [CrossRef] [Green Version]
  23. Gallucci, M.; Carbonara, P.; Pacilli, A.M.G.; di Palmo, E.; Ricci, G.; Nava, S. Use of Symptoms Scores, Spirometry, and Other Pulmonary Function Testing for Asthma Monitoring. Front. Pediatr. 2019, 7, 54. [Google Scholar] [CrossRef] [Green Version]
  24. Cooper, B.G. Limitations to Spirometry Being Performed in “the Office”. Chron. Respir. Dis. 2005, 2, 113–115. [Google Scholar] [CrossRef] [Green Version]
  25. Tepper, R.S.; Wise, R.S.; Covar, R.; Irvin, C.G.; Kercsmar, C.M.; Kraft, M.; Liu, M.C.; O’Connor, G.T.; Peters, S.P.; Sorkness, R.; et al. Asthma Outcomes: Pulmonary Physiology. J. Allergy Clin. Immunol. 2012, 129 (Suppl. 3), S65–S87. [Google Scholar] [CrossRef] [Green Version]
  26. Hernandez-Pacheco, N.; Flores, C.; Oh, S.S.; Burchard, E.G.; Pino-Yanes, M. What Ancestry Can Tell Us About the Genetic Origins of Inter-Ethnic Differences in Asthma Expression. Curr. Allergy Asthma Rep. 2016, 16, 53. [Google Scholar] [CrossRef]
  27. Hernandez-Pacheco, N.; Pino-Yanes, M.; Flores, C. Genomic Predictors of Asthma Phenotypes and Treatment Response. Front. Pediatr. 2019, 7, 6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Tantisira, K.G.; Damask, A.; Szefler, S.J.; Schuemann, B.; Markezich, A.; Su, J.; Klanderman, B.; Sylvia, J.; Wu, R.; Martinez, F.; et al. Genome-Wide Association Identifies the T Gene as a Novel Asthma Pharmacogenetic Locus. Am. J. Respir. Crit. Care Med. 2012, 185, 1286–1291. [Google Scholar] [CrossRef] [Green Version]
  29. Tantisira, K.G.; Lasky-Su, J.; Harada, M.; Murphy, A.; Litonjua, A.A.; Himes, B.E.; Lange, C.; Lazarus, R.; Sylvia, J.; Klanderman, B.; et al. Genome-wide Association between GLCCI1 and Response to Glucocorticoid Therapy in Asthma. N. Engl. J. Med. 2011, 365, 1173–1183. [Google Scholar] [CrossRef] [Green Version]
  30. Wu, A.C.; Himes, B.E.; Lasky-Su, J.; Litonjua, A.; Peters, S.P.; Lima, J.; Kubo, M.; Tamari, M.; Nakamura, Y.; Qiu, W.; et al. Inhaled Corticosteroid Treatment Modulates ZNF432 Gene Variant’s Effect on Bronchodilator Response in Asthmatics. J. Allergy Clin. Immunol. 2014, 133, 723–728.e3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Park, H.W.; Dahlin, A.; Tse, S.; Duan, Q.L.; Schuemann, B.; Martinez, F.D.; Peters, S.P.; Szefler, S.J.; Lima, J.J.; Kubo, M.; et al. Genetic Predictors Associated with Improvement of Asthma Symptoms in Response to Inhaled Corticosteroids. J. Allergy Clin. Immunol. 2014, 133, 664–669.e5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Park, T.J.; Park, J.S.; Cheong, H.S.; Park, B.L.; Kim, L.H.; Heo, J.S.; Kim, Y.K.; Kim, K.U.; Uh, S.T.; Lee, H.S.; et al. Genome-Wide Association Study Identifies ALLC Polymorphisms Correlated with FEV1 Change by Corticosteroid. Clin. Chim. Acta 2014, 436, 20–26. [Google Scholar] [CrossRef] [PubMed]
  33. Dahlin, A.; Denny, J.; Roden, D.M.; Brilliant, M.H.; Ingram, C.; Kitchner, T.E.; Linneman, J.G.; Shaffer, C.M.; Weeke, P.; Xu, H.; et al. CMTR1 Is Associated with Increased Asthma Exacerbations in Patients Taking Inhaled Corticosteroids. Immun. Inflamm. Dis. 2015, 3, 350–359. [Google Scholar] [CrossRef] [PubMed]
  34. Wang, Y.; Tong, C.; Wang, Z.; Mauger, D.; Tantisira, K.G.; Israel, E.; Szefler, S.J.; Chinchilli, V.M.; Boushey, H.A.; Lazarus, S.C.; et al. Pharmacodynamic Genome-Wide Association Study Identifies New Responsive Loci for Glucocorticoid Intervention in Asthma. Pharmacogenom. J. 2015, 15, 422–429. [Google Scholar] [CrossRef] [Green Version]
  35. Leusink, M.; Vijverberg, S.J.; Koenderman, L.; Raaijmakers, J.A.; de Jongste, J.C.; Sterk, P.J.; Duiverman, E.J.; Onland-Moret, N.C.; Postma, D.S.; de Boer, A.; et al. Genetic Variation in Uncontrolled Childhood Asthma despite ICS Treatment. Pharmacogenom. J. 2016, 16, 158–163. [Google Scholar] [CrossRef] [PubMed]
  36. Mosteller, M.; Hosking, L.; Murphy, K.; Shen, J.; Song, K.; Nelson, M.; Ghosh, S. No Evidence of Large Genetic Effects on Steroid Response in Asthma Patients. J. Allergy Clin. Immunol. 2017, 139, 797–803.e7. [Google Scholar] [CrossRef] [PubMed]
  37. Levin, A.M.; Gui, H.; Hernandez-Pacheco, N.; Yang, M.; Xiao, S.; Yang, J.J.; Hochstadt, S.; Barczak, A.J.; Eckalbar, W.L.; Rynkowski, D.; et al. Integrative Approach Identifies Corticosteroid Response Variant in Diverse Populations with Asthma. J. Allergy Clin. Immunol. 2019, 143, 1791–1802. [Google Scholar] [CrossRef]
  38. Hernandez-Pacheco, N.; Farzan, N.; Francis, B.; Karimi, L.; Repnik, K.; Vijverberg, S.J.; Soares, P.; Schieck, M.; Gorenjak, M.; Forno, E.; et al. Genome-Wide Association Study of Inhaled Corticosteroid Response in Admixed Children with Asthma. Clin. Exp. Allergy 2019, 49, 789–798. [Google Scholar] [CrossRef] [Green Version]
  39. Hernandez-Pacheco, N.; Vijverberg, S.J.; Herrera-Luis, E.; Li, J.; Sio, Y.Y.; Granell, R.; Corrales, A.; Maroteau, C.; Lethem, R.; Perez-Garcia, J.; et al. Genome-Wide Association Study of Asthma Exacerbations despite Inhaled Corticosteroids Use. Eur. Respir. J. 2020, 57, 2003388. [Google Scholar] [CrossRef]
  40. Garcia-Menaya, J.M.; Cordobes-Duran, C.; Garcia-Martin, E.; Agundez, J.A.G. Pharmacogenetic Factors Affecting Asthma Treatment Response. Potential Implications for Drug Therapy. Front. Pharmacol. 2019, 10, 520. [Google Scholar] [CrossRef]
  41. Berce, V.; Kozmus, C.E.; Potocnik, U. Association among ORMDL3 Gene Expression, 17q21 Polymorphism and Response to Treatment with Inhaled Corticosteroids in Children with Asthma. Pharmacogenom. J. 2013, 13, 523–529. [Google Scholar] [CrossRef]
  42. Hernandez-Pacheco, N.; Gorenjak, M.; Jurgec, S.; Corrales, A.; Jorgensen, A.; Karimi, L.; Vijverberg, S.J.; Berce, V.; Schieck, M.; Acosta-Herrera, M.; et al. Combined Analysis of Transcriptomic and Genetic Data for the Identification of Loci Involved in Glucocorticosteroid Response in Asthma. Allergy 2020, 76, 1238–1243. [Google Scholar] [CrossRef]
  43. Farzan, N.; Vijverberg, S.J.; Andiappan, A.K.; Arianto, L.; Berce, V.; Blanca-Lopez, N.; Bisgaard, H.; Bonnelykke, K.; Burchard, E.G.; Campo, P.; et al. Rationale and Design of the Multiethnic Pharmacogenomics in Childhood Asthma Consortium. Pharmacogenomics 2017, 18, 931–943. [Google Scholar] [CrossRef]
  44. British Thoracic Society and the Scottish Intercollegiate Guidelines Network. British guideline on the management of asthma. Thorax 2014, 69 (Suppl. 1), 1–192. [Google Scholar]
  45. The ENCODE Project Consortium. An Integrated Encyclopedia of DNA Elements in the Human Genome. Nature 2012, 489, 57–74. [Google Scholar] [CrossRef]
  46. Staley, J.R.; Blackshaw, J.; Kamat, M.A.; Ellis, S.; Surendran, P.; Sun, B.B.; Paul, D.S.; Freitag, D.; Burgess, S.; Danesh, J.; et al. PhenoScanner: A Database of Human Genotype-Phenotype Associations. Bioinformatics 2016, 32, 3207–3209. [Google Scholar] [CrossRef] [Green Version]
  47. Kamat, M.A.; Blackshaw, J.A.; Young, R.; Surendran, P.; Burgess, S.; Danesh, J.; Butterworth, A.S.; Staley, J.R. PhenoScanner V2: An Expanded Tool for Searching Human Genotype-Phenotype Associations. Bioinformatics 2019, 35, 4851–4853. [Google Scholar] [CrossRef] [Green Version]
  48. Suhre, K.; Arnold, M.; Bhagwat, A.M.; Cotton, R.J.; Engelke, R.; Raffler, J.; Sarwath, H.; Thareja, G.; Wahl, A.; DeLisle, R.K.; et al. Connecting Genetic Risk to Disease End Points through the Human Blood Plasma Proteome. Nat. Commun. 2017, 8, 14357. [Google Scholar] [CrossRef] [Green Version]
  49. Sun, B.B.; Maranville, J.C.; Peters, J.E.; Stacey, D.; Staley, J.R.; Blackshaw, J.; Burgess, S.; Jiang, T.; Paige, E.; Surendran, P.; et al. Genomic Atlas of the Human Plasma Proteome. Nature 2018, 558, 73–79. [Google Scholar] [CrossRef] [PubMed]
  50. Folkersen, L.; Fauman, E.; Sabater-Lleal, M.; Strawbridge, R.J.; Frånberg, M.; Sennblad, B.; Baldassarre, D.; Veglia, F.; Humphries, S.E.; Rauramaa, R.; et al. Mapping of 79 Loci for 83 Plasma Protein Biomarkers in Cardiovascular Disease. PLoS Genet. 2017, 13, e1006706. [Google Scholar] [CrossRef] [PubMed]
  51. Tirado-Rodriguez, B.; Ortega, E.; Segura-Medina, P.; Huerta-Yepez, S. TGF- Beta: An Important Mediator of Allergic Disease and a Molecule with Dual Activity in Cancer Development. J. Immunol. Res. 2014, 2014, 318481. [Google Scholar] [CrossRef] [Green Version]
  52. Carvalho-Silva, D.; Pierleoni, A.; Pignatelli, M.; Ong, C.; Fumis, L.; Karamanis, N.; Carmona, M.; Faulconbridge, A.; Hercules, A.; McAuley, E.; et al. Open Targets Platform: New Developments and Updates Two Years On. Nucleic Acids Res. 2019, 47, D1056–D1065. [Google Scholar] [CrossRef] [PubMed]
  53. Li, H.S.; Chen, J.H.; Wu, W.; Fagaly, T.; Zhou, L.; Yuan, W.; Dupuis, S.; Jiang, Z.H.; Nash, W.; Gick, C.; et al. Vertebrate Slit, a Secreted Ligand for the Transmembrane Protein Roundabout, Is a Repellent for Olfactory Bulb Axons. Cell 1999, 96, 807–818. [Google Scholar] [CrossRef] [Green Version]
  54. Kidd, T.; Brose, K.; Mitchell, K.J.; Fetter, R.D.; Tessier-Lavigne, M.; Goodman, C.S.; Tear, G. Roundabout Controls Axon Crossing of the CNS Midline and Defines a Novel Subfamily of Evolutionarily Conserved Guidance Receptors. Cell 1998, 92, 205–215. [Google Scholar] [CrossRef] [Green Version]
  55. Xian, J.; Clark, K.J.; Fordham, R.; Pannell, R.; Rabbitts, T.H.; Rabbitts, P.H. Inadequate Lung Development and Bronchial Hyperplasia in Mice with a Targeted Deletion in the Dutt1/Robo1 Gene. Proc. Natl. Acad. Sci. USA 2001, 98, 15062–15066. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Dickinson, R.E.; Duncan, W.C. The SLIT-ROBO Pathway: A Regulator of Cell Function with Implications for the Reproductive System. Reproduction 2010, 139, 697–704. [Google Scholar] [CrossRef] [Green Version]
  57. Tole, S.; Mukovozov, I.M.; Huang, Y.W.; Magalhaes, M.A.; Yan, M.; Crow, M.R.; Liu, G.Y.; Sun, C.X.; Durocher, Y.; Glogauer, M.; et al. The Axonal Repellent, Slit2, Inhibits Directional Migration of Circulating Neutrophils. J. Leukoc. Biol. 2009, 86, 1403–1415. [Google Scholar] [CrossRef] [Green Version]
  58. Lin, Y.Z.; Zhong, X.N.; Chen, X.; Liang, Y.; Zhang, H.; Zhu, D.L. Roundabout Signaling Pathway Involved in the Pathogenesis of COPD by Integrative Bioinformatics Analysis. Int. J. COPD 2019, 14, 2145–2162. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Branchfield, K.; Nantie, L.; Verheyden, J.M.; Sui, P.; Wienhold, M.D.; Sun, X. Pulmonary Neuroendocrine Cells Function as Airway Sensors to Control Lung Immune Response. Science 2016, 351, 707–710. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Tseng, R.C.; Lee, S.H.; Hsu, H.S.; Chen, B.H.; Tsai, W.C.; Tzao, C.; Wang, Y.C. SLIT2 Attenuation during Lung Cancer Progression Deregulates Beta-Catenin and E-Cadherin and Associates with Poor Prognosis. Cancer Res. 2010, 70, 543–551. [Google Scholar] [CrossRef] [Green Version]
  61. Wu, J.Y.; Feng, L.; Park, H.T.; Havlioglu, N.; Wen, L.; Tang, H.; Bacon, K.B.; Jiang, Z.; Zhang, X.; Rao, Y. The Neuronal Repellent Slit Inhibits Leukocyte Chemotaxis Induced by Chemotactic Factors. Nature 2001, 410, 948–952. [Google Scholar] [CrossRef]
  62. Pilling, D.; Chinea, L.E.; Consalvo, K.M.; Gomer, R.H. Different Isoforms of the Neuronal Guidance Molecule Slit2 Directly Cause Chemoattraction or Chemorepulsion of Human Neutrophils. J. Immunol. 2019, 202, 239–248. [Google Scholar] [CrossRef] [Green Version]
  63. Ye, B.Q.; Geng, Z.H.; Ma, L.; Geng, J.G. Slit2 Regulates Attractive Eosinophil and Repulsive Neutrophil Chemotaxis through Differential SrGAP1 Expression during Lung Inflammation. J. Immunol. 2010, 185, 6294–6305. [Google Scholar] [CrossRef] [Green Version]
  64. Pilling, D.; Zheng, Z.; Vakil, V.; Gomer, R.H. Fibroblasts Secrete Slit2 to Inhibit Fibrocyte Differentiation and Fibrosis. Proc. Natl. Acad. Sci. USA 2014, 111, 18291–18296. [Google Scholar] [CrossRef] [Green Version]
  65. Gaspar Marques, J.; Lobato, M.; Leiria Pinto, P.; Neuparth, N.; Carreiro Martins, P. Asthma and COPD “Overlap”: A Treatable Trait or Common Several Treatable-Traits? Eur. Ann. Allergy Clin. Immunol. 2020, 52, 148–159. [Google Scholar] [CrossRef] [Green Version]
  66. Cukic, V.; Lovre, V.; Dragisic, D.; Ustamujic, A. Asthma and Chronic Obstructive Pulmonary Disease (COPD)—Differences and Similarities. Mater. Sociomed. 2012, 24, 100–105. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Ning, Y.; Sun, Q.; Dong, Y.; Xu, W.; Zhang, W.; Huang, H.; Li, Q. Slit2-N Inhibits PDGF-Induced Migration in Rat Airway Smooth Muscle Cells: WASP and Arp2/3 Involved. Toxicology 2011, 283, 32–40. [Google Scholar] [CrossRef]
  68. Zuyderduyn, S.; Sukkar, M.B.; Fust, A.; Dhaliwal, S.; Burgess, J.K. Treating Asthma Means Treating Airway Smooth Muscle Cells. Eur. Respir. J. 2008, 32, 265–274. [Google Scholar] [CrossRef] [PubMed]
  69. Parameswaran, K.; Radford, K.; Fanat, A.; Stephen, J.; Bonnans, C.; Levy, B.D.; Janssen, L.J.; Cox, P.G. Modulation of Human Airway Smooth Muscle Migration by Lipid Mediators and Th-2 Cytokines. Am. J. Respir. Cell Mol. Biol. 2007, 37, 240–247. [Google Scholar] [CrossRef] [PubMed]
  70. Lutz, S.M.; Cho, M.H.; Young, K.; Hersh, C.P.; Castaldi, P.J.; McDonald, M.L.; Regan, E.; Mattheisen, M.; DeMeo, D.L.; Parker, M.; et al. A Genome-Wide Association Study Identifies Risk Loci for Spirometric Measures among Smokers of European and African Ancestry. BMC Genet. 2015, 16, 138. [Google Scholar] [CrossRef] [Green Version]
  71. Ding, L.; Abebe, T.; Beyene, J.; Wilke, R.A.; Goldberg, A.; Woo, J.G.; Martin, L.J.; Rothenberg, M.E.; Rao, M.; Hershey, G.K.; et al. Rank-Based Genome-Wide Analysis Reveals the Association of Ryanodine Receptor-2 Gene Variants with Childhood Asthma among Human Populations. Hum. Genom. 2013, 7, 16. [Google Scholar] [CrossRef] [Green Version]
  72. Sirugo, G.; Williams, S.M.; Tishkoff, S.A. The Missing Diversity in Human Genetic Studies. Cell 2019, 177, 26–31. [Google Scholar] [CrossRef] [Green Version]
  73. Potaczek, D.P.; Harb, H.; Michel, S.; Alhamwe, B.A.; Renz, H.; Tost, J. Epigenetics and Allergy: From Basic Mechanisms to Clinical Applications. Epigenomics 2017, 9, 539–571. [Google Scholar] [CrossRef]
  74. Alashkar Alhamwe, B.; Miethe, S.; Pogge von Strandmann, E.; Potaczek, D.P.; Garn, H. Epigenetic Regulation of Airway Epithelium Immune Functions in Asthma. Front. Immunol. 2020, 11, 1747. [Google Scholar] [CrossRef] [PubMed]
  75. Reddel, H.K.; Taylor, D.R.; Bateman, E.D.; Boulet, L.P.; Boushey, H.A.; Busse, W.W.; Casale, T.B.; Chanez, P.; Enright, P.L.; Gibson, P.G.; et al. An Official American Thoracic Society/European Respiratory Society Statement: Asthma Control and Exacerbations: Standardizing Endpoints for Clinical Asthma Trials and Clinical Practice. Am. J. Respir. Crit. Care Med. 2009, 180, 59–99. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Tse, S.M.; Gold, D.R.; Sordillo, J.E.; Hoffman, E.B.; Gillman, M.W.; Rifas-Shiman, S.L.; Fuhlbrigge, A.L.; Tantisira, K.G.; Weiss, S.T.; Litonjua, A.A. Diagnostic Accuracy of the Bronchodilator Response in Children. J. Allergy Clin. Immunol. 2013, 132, 554–559.e5. [Google Scholar] [CrossRef] [Green Version]
  77. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.; Daly, M.J.; et al. PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [Green Version]
  78. Chang, C.C.; Chow, C.C.; Tellier, L.C.; Vattikuti, S.; Purcell, S.M.; Lee, J.J. Second-Generation PLINK: Rising to the Challenge of Larger and Richer Datasets. Gigascience 2015, 4, 7. [Google Scholar] [CrossRef]
  79. Kang, H.M. EPACTS (Efficient and Parallelizable Association Container Toolbox) 2016. Available online: http://genome.sph.umich.edu/wiki/EPACTS (accessed on 11 April 2019).
  80. Price, A.L.; Patterson, N.J.; Plenge, R.M.; Weinblatt, M.E.; Shadick, N.A.; Reich, D. Principal Components Analysis Corrects for Stratification in Genome-Wide Association Studies. Nat. Genet. 2006, 38, 904–909. [Google Scholar] [CrossRef] [PubMed]
  81. Zhao, J.H. Gap: Genetic Analysis Package. R Package Version 1.2.2. 2020. Available online: https://CRAN.R-project.org/package=gap (accessed on 24 July 2019).
  82. Abecasis, G.R.; Auton, A.; Brooks, L.D.; DePristo, M.A.; Durbin, R.M.; Handsaker, R.E.; Kang, H.M.; Marth, G.T.; McVean, G.A. An Integrated Map of Genetic Variation from 1,092 Human Genomes. Nature 2012, 491, 56–65. [Google Scholar]
  83. Das, S.; Forer, L.; Schonherr, S.; Sidore, C.; Locke, A.E.; Kwong, A.; Vrieze, S.I.; Chew, E.Y.; Levy, S.; McGue, M.; et al. Next-Generation Genotype Imputation Service and Methods. Nat. Genet. 2016, 48, 1284–1287. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Kuo, K.H.M. Multiple Testing in the Context of Gene Discovery in Sickle Cell Disease Using Genome-Wide Association Studies. Genom. Insights 2017, 10, 1178. [Google Scholar] [CrossRef] [Green Version]
  85. Plummer, M.; Best, N.; Cowles, K.; Vines, K. CODA: Convergence Diagnosis and Output Analysis for MCMC. R News 2006, 6, 7–11. [Google Scholar]
  86. Ward, L.D.; Kellis, M. HaploReg v4: Systematic Mining of Putative Causal Variants, Cell Types, Regulators and Target Genes for Human Complex Traits and Disease. Nucleic Acids Res. 2016, 44, D877–D881. [Google Scholar] [CrossRef] [PubMed]
Jpm 11 00733 g001 550
Figure 1. Manhattan plot of association results of the change in FEV1 after ICS treatment in SLOVENIA. The logarithmic transformation of the p-value (−log10 p-value) is represented on the y-axis, along with the chromosome position (x-axis). The gray line indicates the suggestive significance threshold (p ≤ 5 × 10−6) considered for evidence of association with ICS response.
Figure 1. Manhattan plot of association results of the change in FEV1 after ICS treatment in SLOVENIA. The logarithmic transformation of the p-value (−log10 p-value) is represented on the y-axis, along with the chromosome position (x-axis). The gray line indicates the suggestive significance threshold (p ≤ 5 × 10−6) considered for evidence of association with ICS response.
Jpm 11 00733 g001
Jpm 11 00733 g002 550
Figure 2. Regional plot of association results with the change in FEV1 after ICS treatment expressed as a binary variable in SLOVENIA. Results are represented based on the significance of the association on the y-axis (−log10 p-value) and chromosome position (x-axis) for each SNP as a dot. The SNP rs1166980 suggestively associated with ICS response is represented by the purple diamond. The remaining SNPs are color-coded based on pairwise linkage disequilibrium (r2 values) with that SNP for European populations from 1KGP (GRCh37/hg19 build).
Figure 2. Regional plot of association results with the change in FEV1 after ICS treatment expressed as a binary variable in SLOVENIA. Results are represented based on the significance of the association on the y-axis (−log10 p-value) and chromosome position (x-axis) for each SNP as a dot. The SNP rs1166980 suggestively associated with ICS response is represented by the purple diamond. The remaining SNPs are color-coded based on pairwise linkage disequilibrium (r2 values) with that SNP for European populations from 1KGP (GRCh37/hg19 build).
Jpm 11 00733 g002
Table
Table 1. Clinical and demographic characteristics of the asthma patients from the SLOVENIA study included in the GWAS of change in FEV1 after ICS treatment.
Table 1. Clinical and demographic characteristics of the asthma patients from the SLOVENIA study included in the GWAS of change in FEV1 after ICS treatment.
TotalICS Non-
Responders a		ICS Responders bp-Value
Sample size1669472-
Gender, n (% male)98 (59.0)59 (62.8)39 (54.2)0.264 e
Mean age ± SD (years)10.9 ± 3.410.7 ± 3.211.2 ± 3.50.461 f
Lung function
Mean basal FEV1 ± SD (%) c87.1 ± 14.891.3 ± 12.781.6 ± 15.5<0.001 f
Mean post-treatment FEV1 ± SD (%) d93.7 ± 14.490.1 ± 13.698.5 ± 14.2<0.001 f
Mean ΔFEV1 ± SD (%)6.7 ± 12.1−1.2 ± 7.816.9 ± 8.7<0.001 f
a Asthma patients with ΔFEV1 < 8% after 6 weeks of ICS treatment; b Asthma patients with ΔFEV1 ≥ 8% after 6 weeks of ICS treatment; c FEV1 measured at the beginning of ICS treatment; d FEV1 measured after 6 weeks of ICS treatment; e Pearson χ2 test (df = 1; α = 0.05); f Mann–Whitney U test. FEV1: forced expiratory volume in one second; ΔFEV1: change in FEV1 after 6 weeks of ICS treatment; SD: standard deviation; NA: not available.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hernandez-Pacheco, N.; Gorenjak, M.; Li, J.; Repnik, K.; Vijverberg, S.J.; Berce, V.; Jorgensen, A.; Karimi, L.; Schieck, M.; Samedy-Bates, L.-A.; et al. Identification of ROBO2 as a Potential Locus Associated with Inhaled Corticosteroid Response in Childhood Asthma. J. Pers. Med. 2021, 11, 733. https://doi.org/10.3390/jpm11080733

AMA Style

Hernandez-Pacheco N, Gorenjak M, Li J, Repnik K, Vijverberg SJ, Berce V, Jorgensen A, Karimi L, Schieck M, Samedy-Bates L-A, et al. Identification of ROBO2 as a Potential Locus Associated with Inhaled Corticosteroid Response in Childhood Asthma. Journal of Personalized Medicine. 2021; 11(8):733. https://doi.org/10.3390/jpm11080733

Chicago/Turabian Style

Hernandez-Pacheco, Natalia, Mario Gorenjak, Jiang Li, Katja Repnik, Susanne J. Vijverberg, Vojko Berce, Andrea Jorgensen, Leila Karimi, Maximilian Schieck, Lesly-Anne Samedy-Bates, and et al. 2021. "Identification of ROBO2 as a Potential Locus Associated with Inhaled Corticosteroid Response in Childhood Asthma" Journal of Personalized Medicine 11, no. 8: 733. https://doi.org/10.3390/jpm11080733

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop