Next Article in Journal
Polysomnographically Defined Restless Sleep Disorder and Periodic Limb Movements during Sleep in Children Born Prematurely
Next Article in Special Issue
Mepolizumab in Severe Pediatric Asthma: Certainties and Doubts through a Single-Center Experience and Review of the Literature
Previous Article in Journal
Evaluation of Adolescents’ Awareness of Seat Belt Use and the Relationship with Risky Behaviors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Interleukin-17 Receptor A Gene Variants on Asthma Susceptibility and Clinical Manifestations in Children and Adolescents

1
Department of Medical Biochemistry and Molecular Biology, Faculty of Medicine, Suez Canal University, Ismailia 41522, Egypt
2
Center of Excellence in Molecular and Cellular Medicine (CEMCM), Faculty of Medicine, Suez Canal University, Ismailia 41522, Egypt
3
Department of Pediatrics, Faculty of Medicine, Suez Canal University, Ismailia 41522, Egypt
4
Department of Basic Medical Sciences, Faculty of Medicine, Jazan University, Jazan 45141, Saudi Arabia
5
Department of Chest Diseases and Tuberculosis, Faculty of Medicine, Suez Canal University, Ismailia 41522, Egypt
6
Department of Medical Laboratory, College of Applied Medical Sciences, Taibah University, Yanbu 46423, Saudi Arabia
7
Department of Clinical Pathology, Faculty of Medicine, Suez Canal University, Ismailia 41522, Egypt
8
Department of Surgery, School of Medicine, Tulane University, New Orleans, LA 70112, USA
9
Medical Genetics Unit, Department of Histology and Cell Biology, Suez Canal University, Ismailia 41522, Egypt
10
Biochemistry Department, Animal Health Research Institute, Mansoura Branch, Giza 12618, Egypt
11
Department of Biochemistry, Faculty of Medicine, Northern Border University, Arar P.O. Box 1321, Saudi Arabia
*
Author to whom correspondence should be addressed.
Children 2024, 11(6), 657; https://doi.org/10.3390/children11060657
Submission received: 16 April 2024 / Revised: 24 May 2024 / Accepted: 26 May 2024 / Published: 28 May 2024

Abstract

:
Several single nucleotide polymorphisms (SNPs) in multiple interleukin receptor genes could be associated with asthma risk and/or phenotype. Interleukin-17 (IL-17) has been implicated in tissue inflammation and autoimmune diseases. As no previous studies have uncovered the potential role of IL17 receptor A (RA) gene variants in asthma risk, we aimed to explore the association of four IL17RA SNPs (i.e., rs4819554A/G, rs879577C/T, rs41323645G/A, and rs4819555C/T) with asthma susceptibility/phenotype in our region. TaqMan allelic discrimination analysis was used to genotype 192 individuals. We found that the rs4819554 G/G genotype significantly reduced disease risk in the codominant (OR = 0.15, 95%CI = 0.05–0.45, p < 0.001), dominant (OR = 0.49, 95%CI = 0.26–0.93, p = 0.028), and recessive (OR = 0.18, 95%CI = 0.07–0.52, p < 0.001) models. Similarly, rs879577 showed reduced disease risk associated with the T allele across all genetic models. However, the A allele of rs41323645 was associated with increased disease risk in all models. The G/A and A/A genotypes have higher ORs of 2.47 (95%CI = 1.19–5.14) and 3.86 (95%CI = 1.62–9.18), respectively. Similar trends are observed in the dominant 2.89 (95%CI = 1.47–5.68, p = 0.002) and recessive 2.34 (95%CI = 1.10–4.98, p = 0.025) models. For the rs4819555 variant, although there was no significant association identified under any models, carriers of the rs4819554*A demonstrated an association with a positive family history of asthma (71.4% in carriers vs. 27% in non-carriers; p = 0.025) and the use of relievers for >2 weeks (52.2% of carriers vs. 28.8% of non-carriers; p = 0.047). Meanwhile, the rs4819555*C carriers displayed a significant divergence in the asthma phenotype, specifically atopic asthma (83.3% vs. 61.1%; p = 0.007), showed a higher prevalence of chest tightness (88.9% vs. 61.5%; p = 0.029), and were more likely to report comorbidities (57.7% vs. 16.7%, p = 0.003). The most frequent haplotype in the asthma group was ACAC, with a frequency of 22.87% vs. 1.36% in the controls (p < 0.001). In conclusion, the studied IL17RA variants could be essential in asthma susceptibility and phenotype in children and adolescents.

Graphical Abstract

1. Introduction

Asthma is a prevalent chronic respiratory illness in pediatrics, characterized by mucus hypersecretion, persistent airway inflammation, and bronchial hyper-responsiveness, leading to the narrowing of the bronchi and increased resistance to airflow, often reversible naturally or with treatment [1]. It is prevalent worldwide, especially in developed countries where the prevalence is increasing to epidemic proportions. According to the statistical survey of the “Global Asthma Network’s Phase 1”, the prevalence of asthma symptoms was 9.1%, 11.0%, and 6.6% for children, adolescents, and adults, respectively [2].
Asthma is a multifaceted disease influenced by various factors, including genetics. Discovering genes associated with asthma, understanding genetic variations, and examining cellular signaling pathways could be crucial for predicting disease outcomes and developing targeted treatments [3]. It is considered a classic T-cell-mediated condition driven by various cytokines orchestrating the inflammatory response. Predominantly, these include the Th2-secreted cytokines such as interleukin-4 (IL-4), IL-5, and IL-13. Additionally, IL-17 plays a notable role in the development and severity of asthma [4].
The IL-17 cytokine family, which significantly influences both acute and chronic immune responses [5], is composed of six members: IL-17A, B, C, D, E, and F. Within this family, IL-17A and IL-17F are linked explicitly to T helper type 17 (Th17) cells [6]. These cytokines signal through five receptors named IL-17 receptors (IL-17R) A–E. Among these, IL-17A primarily signals via the IL-17R complex, which includes the IL-17RA and IL-17RC subunits, with IL-17RA displaying a higher binding affinity [7,8] (Figure 1) [9,10].
Interleukin (IL)-17A, the primary cytokine produced by T helper 17 (Th17) cells, plays a crucial role in defending host mucosal surfaces [11] and acts as a significant inflammatory mediator in various autoimmune diseases, such as rheumatoid arthritis (RA) [12,13]. The IL-17A/IL-17RA signaling pathway has been associated with the development of several acute and chronic respiratory conditions, including asthma [14,15,16,17,18,19,20]. Research has shown that IL-17A promotes the growth and survival of airway smooth muscle cells in both asthmatic and non-asthmatic individuals through IL-17RA activation [21]. Additionally, IL-17A can trigger the production of inflammatory cytokines by human bronchial fibroblasts, including IL-6, IL-8, and IL-11, potentially leading to further inflammation and tissue remodeling [16].
In patients with asthma, there is an upregulation of IL-17 mRNA and proteins found within the lungs, sputum, and bronchoalveolar lavage fluids [4,22]. Research has shown that respiratory infections caused by viruses, bacteria, and fungi can impact the Th17/IL-17A pathway, consequently exacerbating both the onset and progression of asthma [23,24].
Previous studies have indicated that various single nucleotide polymorphisms (SNPs) in several interleukin and interleukin receptor genes, including IL-4 and IL-4R [25,26,27,28], IL-9 [29], IL-10 [30], IL-5 and IL5R [31], IL13 [32], IL-33 and its receptor IL-1R-like 1 [33,34], IL-18 [35], IL-13R [36], and IL-23R [37], may be linked to an increased risk of developing asthma.
Although several SNPs related to IL-17A have been linked to susceptibility of a variety of immune-mediated/inflammatory disorders [38,39,40,41,42,43,44], including asthma [45,46,47,48,49], with potential modification of related gene expression and Th17 cells activities [40,50], no previous studies have explored the potential role of IL17RA gene variants in asthma susceptibility, in particular in our population. In this sense, the current study was designed to investigate for the first time the association between four SNPs (i.e., rs4819554, rs879577, rs41323645, and rs4819555) in IL-17RA and asthma susceptibility/phenotype in our region. These SNPs were chosen according to the criteria outlined in the Methods section.

2. Materials and Methods

2.1. Population Characteristics

Pediatric patients, ranging in age from 6 to 15 years, diagnosed with asthma were included from the “Pediatrics, Allergy, and Immunology, as well as Chest Diseases and Tuberculosis” clinics at the Suez Canal University Hospital in Ismailia, Egypt. The research was carried out among two distinct groups, with 96 unrelated patients diagnosed with bronchial asthma, comprising the asthma group and 96 age-similar healthy individuals (i.e., the control group). Diagnoses were formulated following the Global Initiative for Asthma (GINA) regulations. Notable factors like co-existing health conditions, the extent of control over asthma symptoms over the previous month, disease seriousness, treatment level and responsiveness, and consistency of therapy adherence were evaluated by the treating physicians based on GINA’s guidelines [51]. To ensure accuracy and standardization in line with GINA recommendations, we implemented a secondary review process. Prior to the inclusion of participants in the study, the initial evaluations made by the treating physicians underwent an additional layer of verification by a separate team of respiratory specialists. This team was not involved in the initial clinical care of the patients, providing an independent confirmation of the asthma classification.
Chest X-rays were taken to rule out alternate causes of airflow restrictions. Asthma patients with chronic comorbidities, those who exhibited signs of respiratory infection in the previous three weeks, those treated with systemic steroids in the past fortnight, or those non-adherent to therapy were not included in the study for several reasons: (1) impact on test results: steroid therapy can suppress inflammation and modify airway responsiveness, potentially altering the results of bronchial challenge tests. Excluding patients recently receiving steroid therapy helps ensure the test accurately reflects their underlying respiratory condition. (2) Consistency of results: patients not adherent to their medication regimen may have variable airway inflammation and responsiveness, making it challenging to interpret the results of bronchial challenge tests reliably. Excluding non-adherent patients helps maintain consistency in the study population and improves the validity of the findings. (3) Safety concerns: in some cases, patients not adherent to their medication regimen may have poorly controlled asthma or other respiratory conditions, putting them at higher risk of adverse events during bronchial challenge testing. Excluding these patients helps mitigate potential safety concerns and ensures the ethical conduct of the study.
Control participants were children accompanying their siblings to the Pediatrics Clinic, with no precedent of atopic disorders, history of medications, or any health problems, including chronic respiratory diseases. The “body mass index (BMI)” was determined by dividing the weight of the individual (in kilograms) by their height (in meters) raised to the power of two [52]. Pubertal development was evaluated using the Tanner stages of growth [53].
Adherence to the “Declaration of Helsinki” outlines was ensured. The research gained approval from the “Medical Research Ethics Committee” at the Faculty of Medicine, Suez Canal University, Ismailia, Egypt, and parents or guardians of the study participants provided their written consent.

2.2. Pulmonary Function and Methacholine Challenge Tests

Baseline pulmonary function tests were conducted using the ExpAir spirometry system from Medisoft (Viasys Healthcare, Conshohocken, PA, USA), following the guidelines of the American Thoracic Society (ATS) and European Respiratory Society (ERS) [54]. The parameters measured included “peak expiratory flow rate (PEFR)”, “forced vital capacity (FVC)”, and the “forced expiratory volume in one second (FEV1)”. Participants were administered 400 µg of salbutamol (known commercially as Ventolin by GlaxoSmithKline) using a metered-dose inhaler and a spacer. Forced spirometry was repeated 15 min after administering the bronchodilator salbutamol. The response to salbutamol was calculated using the following formula: “[BDRBASE = ((post-salbutamol FEV1 − pre-salbutamol FEV1)/pre-salbutamol FEV1) × 100]” [55].
Airway hyperresponsiveness was evaluated using a methacholine challenge test, following the clinical practice guidelines [56]. The amount of methacholine required to induce a 20% reduction in FEV1 was determined and recorded in mg/mL. A PC20 value of less than 8 mg/mL indicated positive airway hyperresponsiveness [57].

2.3. Laboratory Investigations

Venous blood samples were collected in EDTA and plain tubes. The EDTA tubes were used for absolute peripheral blood eosinophil count using the “Coulter counting method (Beckman Coulter Ltd., Brea, CA, USA)” and confirmed by manual microscopic differential count. Eosinophil counts exceeding 400 cells/mm3 indicated absolute eosinophilia [58]. Total serum IgE levels were measured from the plain tubes using an enzyme-linked immunosorbent assay (ELISA) with the “Accu-Bind® ELISA kit (Monobind Inc., Lake Forest, CA, USA)”. Total IgE levels of ≥100 IU/mL were classified as an elevated level.

2.4. Genomic DNA Extraction and IL17RA Variants Allelic Discrimination Analysis

Genomic DNA isolation was performed from peripheral blood leukocytes utilizing the “Wizard® Genomic DNA Purification Kit (Promega Co., Madison, WI, USA)”. Its quantity was validated through absorptance measurement at 260/280 nm using a “NanoDrop-1000 spectrophotometer (NanoDrop Tech., Wilmington, NC, USA)” to ensure all the included samples in the subsequent analyses have a ratio of 1.8 to 2.0. The DNA samples that did not meet this criterion were either re-purified or excluded from further analysis. Then, the extracted DNA samples were maintained at −20 °C [59].
Based on searching the IL17RA gene variants in the “National Center of Biotechnology Information Single nucleotide polymorphism database (dbSNP; www.ncbi.nlm.nih.gov)”, (accessed on 20 September 2022) for a minor allele frequency (MAF) ≥ 0.1, documented relationships to other similar diseases, and lack of previous studies exploring the association of these variants with asthma, in particular in our population, the authors selected the studied SNPs. The genotyping of these SNPs was executed using “real-time polymerase chain reaction (PCR)” allelic discrimination TaqMan assays (cat# 4351379). The assay IDs and the context sequence of each variant with the VIC/FAM-identified polymorphisms are presented in Table 1.
The PCR was executed in a 25 µL reaction volume that includes 20 ng of the extracted DNA, 12.5 µL of TaqMan Universal PCR Master Mix (comprising AmpErase UNG, AmpliTaq Gold enzyme, dNTPs, and reaction buffer), and 1.25 µL of TaqMan SNP genotyping assay (Table 1). All plates had negative controls to exclude cross-contamination. PCR amplification was conducted using the “StepOne™ real-time PCR system (Applied Biosystems, Waltham, MA, USA)”, based on two initial holds (50 °C for 2 min and 95 °C for 10 min) followed by a 40-cycle two-step PCR (95 °C denaturation for 15 s and annealing/extension 60 °C for 1 min) [60]. A randomly chosen 10% of the samples were re-genotyped across different runs to preclude false genotype calls, maintaining a 100% concordance rate. Post-PCR data analysis was carried out using “SDS software (v1.3.1., Applied Biosystems)”.

2.5. Statistical Analysis

Statistical analysis was conducted using SPSS version 27.0 (IBM Corp., Armonk, NY, USA). Study power was figured out employing G*Power software (version 3.0.10), revealing an anticipated power of 91% given a total sample size of 192, a medium effect size of 0.3, and an error probability (α) of 0.05. The calculation of allele/genotype and carriage rate frequencies of IL17RA adopted the method delineated previously [29]. The “Hardy–Weinberg equilibrium (HWE)” was calculated for each studied variant. Computed odds ratios (OR) with a 95%CI (confidence interval) across varied genetic association models were accomplished through logistics regression models [30]. The potential confounding effects were adjusted using binary logistic regression analysis. The relationship between IL17RA SNPs and clinic-laboratory features was assessed using two-sided statistical tests, including the “Chi-square test, t-test, Mann–Whitney U test, or Kruskal–Walli’s test”. A significance level of p < 0.05 (two-tailed) was considered statistically significant to determine associations.

2.6. In Silico Analysis

The molecular structure of the IL17RA gene/protein was determined by searching several available online tools such as the “GeneCards (https://www.genecards.org/), the National Center of Biotechnology Institute (https://www.ncbi.nlm.nih.gov/gene/23765), the Ensembl genome browser 110 (https://asia.ensembl.org/Homo_sapiens/Transcript/), the Protter tool for visualization of proteoforms and interactive integration of annotated/predicted sequence features of the protein (https://wlab.ethz.ch/protter/)” [61], the “UniProt protein database for the resource of protein sequence and functional information (https://www.uniprot.org/)”, “the subcellular/molecular localization of the expressed protein of interest (https://compartments.jensenlab.org/)”, and the “known/predicted protein–protein interactions network generation (https://string-db.org/)” (all previous databases were last accessed on 20 July 2023).
The functional effect of the studied SNPs was checked using the following bioinformatics tools: “MutationTaster (https://www.mutationtaster.org/, accessed on 20 July 2023)” [62]; “SIFT (Sorting Intolerant from Tolerant, https://sift.bii.a-star.edu.sg/, accessed on 20 July 2023)”, which depends on sequence homology/physical properties of amino acids to predict the effect of missense variants on protein function [63]; PANTHER, that provides “comprehensive information about the evolution of protein-coding gene families depending on protein phylogeny, function and genetic variation impacting that function” [64]; Varsome, which generates recommended pathogenicity based on the data available from multiple genomic databases [65]; “I-Mutant v2.0 predictor of protein Stability changes upon base substitution (https://folding.biofold.org/cgi-bin/i-mutant2.0.cgi, accessed on 20 July 2023)” [66]; and “HaploReg v4.2 to predict the impact of variants on regulatory motifs, and the expression from expression quantitative trait loci (eQTL) studies” [67].

3. Results

3.1. Participants Characteristics

The study consisted of 96 controls and 96 asthma patients, with an average age of around 9.7 years for both groups (p = 0.76). The sex distribution was balanced, with females accounting for 58.3% and 54.2% in the control and asthma groups, respectively (p = 0.66). A higher proportion of rural residents was observed in the asthma group (66.7% vs. 53.1%, p = 0.08). Notably, a higher percentage of individuals in the asthma group reported a family history of asthma (30.2% vs. 7.3%, p < 0.001). BMI percentiles were similar across both groups, and pubertal status and Tanner stages were comparable, indicating similar physical development levels (p = 0.88 and p = 0.65, respectively) (Table 2).

3.2. Association of IL17RA Polymorphisms with Disease Risk

Table 3 presents the relationship between IL17RA polymorphisms and disease risk across four genetic models—codominant, dominant, recessive, and over-dominant—for four different polymorphisms (rs4819554, rs879577, rs41323645, and rs4819555). The distribution of genotypes in the control group followed the HWE (p-values = 0.42, 1.0, 0.27, and 0.82), respectively.
The G/G genotype of rs4819554 significantly reduced disease risk in the codominant model (OR = 0.15, 95% = 0.05–0.45, p < 0.001) compared to the reference A/A genotype. This protective effect of the G allele was also observed in dominant (OR = 0.49, 95%CI = 0.26–0.93, p = 0.028) and recessive (OR = 0.18, 95%CI = 0.07–0.52, p < 0.001) models. Similarly, rs879577 showed a reduction in disease risk associated with the T allele, but it was consistent across all genetic models. In the codominant model, the T/T genotype (OR = 0.13, 95%CI = 0.05–0.36, p < 0.001) was associated with a lower disease risk than the reference C/C genotype. The dominant model showed a similar trend (OR = 0.19, 95%CI = 0.10–0.37, p < 0.001). The recessive model indicated a lower disease risk for the T/T genotype (OR = 0.28, 95%CI = 0.11–0.70, p = 0.004) compared to C/C-C/T genotypes. The A allele of rs41323645 was associated with increased disease risk in all models. The G/A and A/A genotypes had higher ORs of 2.47 (95%CI = 1.19–5.14) and 3.86 (95%CI = 1.62–9.18), respectively. Similar trends were observed in the dominant 2.89 (95%CI = 1.47–5.68, p = 0.002) and recessive 2.34 (95%CI = 1.10–4.98, p = 0.025) models, indicating an increased disease risk associated with the A allele. For rs4819555, no significant associations were identified in any models.
Table 4 reports the carriage rates of alleles, meaning the “proportion of individuals carrying at least one copy of the specific allele in question”. The alleles of four different polymorphisms (rs4819554, rs879577, rs4819554, and rs4819555) are compared between the control/patient subgroups. For rs4819554, carriers of the A variant were associated with higher disease risk (OR = 4.01, p = 0.002), while those of the G allele appeared protected (OR = 0.45, p = 0.012). Similarly, the rs879577*C allele was associated with a higher disease risk (OR = 3.87, p = 0.002), and the T allele was associated with a lower disease risk (OR = 0.22, p < 0.001). For rs41323645, the G allele appeared protective (OR = 0.41, p = 0.016), while the A allele was associated with higher disease risk (OR = 2.47, p = 0.006). No significant associations were observed for rs4819555 (p = 0.32 and p = 0.24, respectively).
Subgroup analysis was performed based on puberty status under a dominant model, where the presence of one or more copies of the minor allele is considered. The genotype distributions between patients and controls for four SNPs (rs4819554, rs879577, rs41323645, and rs4819555) were compared for both prepubertal and pubertal groups (Table 5). For rs4819554, in the prepuberty group, having one or more G alleles (A/G or G/G genotypes) was associated with a reduced risk of disease (OR = 0.49, 95%CI: 0.26–0.93) compared to the A/A genotype. This association, however, was not significant in the pubertal group (OR = 0.70, 95%CI: 0.27–1.77). For rs879577, individuals with one or more T allele (C/T or T/T genotypes) had a significantly lower risk of disease in both prepubertal (OR = 0.11, 95%CI: 0.04–0.30) and pubertal (OR = 0.33, 95%CI: 0.13–0.87) groups, relative to those with the C/C genotype. For rs41323645, individuals with one or more A allele (G/A or A/A genotypes) had a higher risk of disease in the prepubertal group (OR = 3.35, 95%CI: 1.30–8.62). This trend was also observed in the pubertal group, but the association was insignificant (OR = 2.47, 95%CI: 0.95–6.43). For rs4819555, no significant associations were detected in either the prepubertal or pubertal groups. However, at least one T allele (C/T or T/T genotypes) was associated with a trend towards higher risk compared to the C/C genotype. Overall, there was no significant interaction between these SNPs and puberty status (p > 0.05), suggesting that puberty does not significantly modify the effect of these SNPs on disease risk.

3.3. Haplotype Analysis and Disease Risk

Haplotype frequency analysis was performed to identify potential combinations of alleles across the SNPs that might be associated with asthma risk. The most frequent haplotype in the asthma group was ACAC, with a frequency of 22.87%, significantly higher than in the control group (1.36%, p < 0.001). This indicates that the ACAC haplotype could be linked to an increased risk of asthma (Table 6). Other haplotypes, including ACGC, ATGC, and GCAC, were also significantly associated with asthma, albeit with lower frequencies than ACAC. Notably, the ACGC haplotype was associated with a reduced risk of asthma (OR = 0.00, 95%CI: 0.00–0.11, p = 0.002). The ATGC/GCAC haplotypes exhibited comparable trends. These haplotypes were less frequently observed in the asthma group than in the control group, and their presence was linked to a notably reduced likelihood of developing asthma (OR = 0.00, p < 0.001 for ATGC; OR = 0.00, p = 0.002 for GCAC). The haplotypes ACGT, ATGT, ACAT, GTAC, GCAT, GTAT, and GTGT were present at different frequencies in asthma and control groups and were significantly associated with lower odds of having asthma (p < 0.05). In contrast, haplotype GCGC, GTGC, and ATAC showed no significant association with risk of asthma (p > 0.05).

3.4. Association of Gene Variants with Disease Severity

In exploring gene variants’ association with disease severity in pediatric asthma, distinct correlations were noted among the carriers of specific IL17RA gene variants. While many patient characteristics and disease outcomes exhibited no significant differences between carriers and non-carriers, several exceptions were identified (Table 7).
Notably, carriers of the rs4819554*A allele demonstrated a significant association with a positive family history of asthma, contrasting 71.4% in carriers to 27% in non-carriers (p = 0.025). Similarly, rs41323645*A carriers showed a significant correlation with the use of relievers for more than two weeks, with 52.2% of carriers as opposed to 28.8% of non-carriers exhibiting this usage pattern (p = 0.047).
Carriers of the rs4819555*C allele displayed a significant divergence in the asthma phenotype, specifically atopic asthma, identified in 83.3% of carriers compared to 61.1% of non-carriers (p = 0.007). The same group showed a higher prevalence of chest tightness (88.9% in carriers versus 61.5% in non-carriers, p = 0.029). Furthermore, carriers of the rs4819555*C allele were significantly more likely to report comorbidities (57.7% versus 16.7%, p = 0.003), with allergic rhinitis and eczema being more common in this group (44.9% and 38.5% respectively, compared to 16.7% in non-carriers, p = 0.033 for allergic rhinitis and p = 0.001 for eczema). Lastly, a significant difference in eosinophil count was observed between rs4819555*C carriers and non-carriers, with a median count of 50 (IQR: 20–207.5) × 106/L in carriers compared to 230 (IQR: 86–509.8) × 106/L in non-carriers (p = 0.009).
Supplementary Table S1 delineates the relationship between the carriage of IL17RA alleles and various allergies or aggravating factors contributing to asthma. Notably, emotional stress triggers were significantly associated with the rs41323645*A carriers (p = 0.047). Aside from this, our analysis indicates that the studied gene variants do not significantly affect asthma development in the context of these risk factors.

3.5. Association of Gene Variants with Therapeutic Response

As depicted in Table 8 on the carriage of risk alleles, the total number of patients varied for each risk allele, with 162 carrying rs4819554*A, 159 carrying rs879577*C, 127 carrying rs41323645*A, and 162 carrying rs4819555*C. The table compares the management of airway hyperresponsiveness (AHR) and therapy levels for each group. Across the spectrum of normal to severe AHR and therapy levels from Step 1 to Step 5, there were no significant differences in frequencies between carriers and non-carriers of the risk alleles, indicated by p-values above 0.05. The therapeutic response was assessed using several parameters, including FVC1, Pre- and Post-FEV1, PEFR1, “provocative concentration causing a 20% drop in FEV1 (PC20)”, and “bronchodilator response (BDRBASE)”. Notably, carriers of rs4819554*A marginally significantly lower values for Post FEV1 (% predicted) and PEFR1 (% predicted) than non-carriers, with p-values of 0.045 and 0.046, respectively.

3.6. Association of IL17RA Haplotypes with Asthma Severity and Therapeutic Response

Table 9 details a comparison of demographic characteristics, disease severity and control, management strategies, and therapeutic responses between two groups differentiated by the presence or absence of the risk haplotype block (rs4819554*A/rs879577*C/rs41323645*A). The last SNP rs4819555 (C/T) was excluded, as it did not show a significant association. This haplotype block is henceforth referred to as ACA (C/T). The total sample comprised 96 participants, 36 without the risk haplotype and 60 with the risk haplotype. The analysis demonstrated insignificant differences (p-value > 0.05) in these characteristics between the two groups.

3.7. Molecular Features of IL17RA and In Silico Analysis of the Impact of the Studied IL17RA Variants

The IL-17 RA gene (ID: 23765, ENSG00000177663) maps to the long arm of the human chr22 at the locus 22q11.1 spanning 17,084,954 to 17,115,693 on the plus strand (Figure 2A). It has 13 exons and eight transcripts (splice variants), with two codes for proteins (i.e., IL17RA-201/204) (Figure 2B). The cell membrane attached protein consists of 866aa, and the other secretory one consists of 832aa. The cell membrane-attached receptor comprises a signal peptide and the extracellular, transmembrane, and cytoplasmic domains (Figure 2C). The predominant subcellular localization is the cell membrane as a single pass type I protein and the extracellular secretory isoform (Figure 2D). IL17RA exhibits substantial co-expression with IL17F, which acts as a ligand for both IL17RA and IL17RC. The heterodimer formed by IL17A/IL17F is a ligand for the heterodimeric receptor complex composed of IL17RA/IL17RC. This interaction is integral in promoting the production of various cytokines, including IL6, IL8, and CSF2 [68]. Furthermore, IL17RA is involved in the activation of peripheral blood mononuclear cells, T-cell proliferation, and the inhibition of angiogenesis. It also plays a crucial role in inducing lung neutrophilia and exacerbating antigen-induced pulmonary allergic inflammation. Additionally, IL17RA is co-expressed with IL17RE, a specific receptor for IL17C, and is implicated in signal transduction via the NF-kappa-B and MAPK pathways, likely requiring TRAF3IP2/ACT1 for this signaling. This receptor may be essential in the innate immune defense against bacterial infections. The second and fourth isoforms of IL17RA may reside in inactive states within the cytoplasm, while the third and fifth isoforms might function as soluble decoy receptors (Figure 2E,F) (https://string-db.org/) (last accessed on 20 July 2023). The IL17RA gene ontology (i.e., cellular component, molecular functions, and biological process) is represented in Supplementary Table S2. The impact of the studied IL17RA variants on gene expression and/or function is illustrated in Table 10.

4. Discussion

Asthma is a chronic condition characterized by immune-mediated inflammation of the airways, leading to airway remodeling, hyperreactivity, increased mucus production, and airflow obstruction. These changes can cause symptoms that range from mild respiratory discomfort to severe, life-threatening respiratory distress [69]. Various cells and cytokines play critical roles in this process. However, the exact contributions of individual cytokines and their genetic variants to the development and severity of pediatric asthma are not yet fully understood [70,71,72].
Extensive research on genetic susceptibility to asthma has primarily targeted genes encoding polymorphic proteins, particularly those regulating immune responses [73]. Over 60 genetic loci have been linked to asthma, and genome-wide association studies (GWAS) in children and adults with asthma have identified several genes associated with immune responses [73,74,75]. These studies also found that some of these genes encode cytokines that act as immune modulators, and variations in these genes may affect one’s susceptibility to asthma [76].
IL17A is one of these cytokines and has been shown to induce structural changes in airway smooth muscle cells, which may affect susceptibility to airway hyperreactivity (AHR) and the development of asthma [77]. Specific genetic variations in IL17A, such as rs2275913 G/A, rs3819024 A/G, and rs8193036 C/T, have been linked to asthma risk. The rs2275913 variant is associated with an increased risk of developing asthma, whereas the rs8193036 variant may have protective effects [77].
The IL17 receptor A-coding gene (IL17RA; ENSG00000177663) is highly polymorphic. Several IL-17RA variations have been discovered as risky in several disorders’ pathogenesis, most of which are immune-mediated, including autoimmune diabetes [78], psoriasis [79,80,81], chronic spontaneous urticaria (CSU) [82], alopecia areata [83], ankylosing spondylitis [84]. The impact of genetic variation can vary from no effect to alteration of protein expression, processing/maturation, and activity [85].
Evidence suggests a potential link between IL-17RA SNPs and respiratory diseases. Specifically, genetic variants of IL-17R (rs882643 and rs2241049) have been associated with primary graft dysfunction in lung transplantation [86] and aspirin-exacerbated respiratory disease [87,88].
To the authors’ knowledge, this study is the first to interrogate the associations between four IL17RA SNPs: rs4819554A/G, rs879577C/T, rs41323645G/A, and rs4819555C/T with asthma susceptibility/phenotype in an Egyptian population.
Currently, the IL17RA rs4819554 shows that the G/G genotype reduces disease risk significantly, and the G allele is protective. This association was observed in the dominant/codominant and recessive models compared to the reference A/A genotype. Carriers of the A allele for this variant had a four-fold higher disease risk and demonstrated a significant association with a positive asthma FH and using the reliever for more than two weeks. Analysis of risk allele carriage and association with patient characteristics/disease outcome and therapeutic response refile that risk allele was associated with lower values of two parameters to assess therapeutic response, PEFR1 and post FEV1.
The rs4819554 SNP is a common IL17RA variant discovered upstream in the gene (c.-947G > A), and because it is located in the IL17RA promoter region, it may have a functional effect by affecting gene transcription. Compared to the other genotypes, the rs4819554 AA genotype enhanced protein production. This genotype may regulate inflammatory responses by boosting the production of TNF and other proinflammatory cytokines [89]. Interestingly, this SNP is in linkage disequilibrium with other promoter polymorphisms that influence transcription factor binding, such as the binding sites for Lyf-1 and Ik-2. These transcription factors are part of the Ikaros (IK) family, which has been associated with Th17 cell differentiation by upregulating the expression of Th17 lineage-determining genes, including IL-17A, IL-17F, IL-21, IL-22, and IL-23R, while also suppressing the expression of genes that inhibit Th17 cell development [90].
The rs4819554 has been related to enhanced vulnerability or protection against various inflammatory and immunological disorders, as well as some types of cancer. P Park et al. observed that this SNP may predispose individuals to asthma by upregulating IL17RA expression. Asthmatics exhibit higher IL-17RA mRNA and protein levels in CD14+ monocytes and mononuclear cells. This upregulation was significantly more pronounced in individuals with the AA genotype than those with the GG genotype, which aligns with the decreased frequency of the minor allele G in asthmatic patients [88].
In a study conducted by Lauhkonen and colleagues to examine the relationship between three IL17RA polymorphisms, including rs4819554, and various outcomes such as asthma, asthma medication use, allergic rhinitis, and lung function decline at 5–7 years and 11–13 years after hospitalization for bronchiolitis in the first six months of infancy, no correlation was found between the IL17RA gene variations studied and susceptibility to severe bronchiolitis in infancy or subsequent post-bronchiolitis asthma [91]. Additionally, autoimmune type 1 diabetes and other inflammatory and autoimmune disorders have been linked to this variation. The GG genotype was discovered to be linked with a lower risk than the AA and AG genotypes and to be negatively correlated with the presence of serum “glutamic acid decarboxylase antibodies” [78].
The rs4819554 SNP was significantly associated with an increased risk of psoriasis. Psoriatic patients exhibited higher frequencies of the G allele and AG genotype than controls. Additionally, patients with AA or AG genotypes had significantly elevated gene transcript and protein levels relative to those in normal controls. This suggests that the AA genotype of rs4819554 may be linked to a reduced likelihood of acquiring psoriasis [79]. The rs48419554 can impact protein biological activity and thus modify immunological aspects such as response to treatment. IL-17RA rs48419554 was previously related to etanercept response in psoriatic arthritis [92]. Using the pharma GKB database to explore clinical annotation, which provides information about variant-drug connects based basically on variant annotations for rs4819554, it was revealed that patients with the AA genotype and psoriasis might have a better response to treatment with anti-TNF drugs than patients with the AG or GG genotype and carriers of the GG genotype have the poorer response (level 3 efficacy) (http://www.pharmgkb.org) (last accessed on 30 July 2023)
Sandip et al. identified a significant association between the rs4819554 variant and cardiovascular mortality in patients with congestive heart failure, with a hazard ratio (HR) of 1.28 (95% CI = 1.02–1.59, adjusted p = 0.03) [93]. Another study reported that, in patients with coronary artery disease (CAD), the IL17RA transcript levels were elevated in total leukocytes, and rs4819554*G carriers exhibited higher transcript levels compared to controls despite no significant association being found between the rs4819554 variant and CAD risk [94]. The GG genotype of rs4819554 has also been reported as a risk factor for preeclampsia [95]. Additionally, the G allele and the combined AG + GG genotypes were strongly associated with Hashimoto’s thyroiditis, particularly in cases with a positive family history and increased illness severity [96].
The IL-17RA rs48419554 G allele may be a potential marker of disease severity in Polish ankylosing spondylitis patients [97]. Additionally, patients with chronic spontaneous urticaria (CSU) had significantly more GG genotypes than controls, and the frequency of the A allele was notably lower in urticaria patients. The G allele was considered a risk factor. Despite this, patients with the rs4819554 AA genotype were associated with longer disease duration, concurrent angioedema, positive autologous serum skin test status, advanced treatment stages, and the lowest quality of life scores [82]. This may be partially explained by an increase in the level of the IL17RA protein with the A allele, which causes the condition to become more severe [88,98].
Moreover, the rs4819554 variant was significantly associated with worse progression-free survival and overall survival in patients with non-small cell lung cancer (NSCLC) [99]. Similarly, homozygosity for the AA genotype was linked to more significant deterioration of renal function [98] and increased primary graft dysfunction after lung transplantation [86]. Conversely, the rs4819554 AG genotype was associated with smaller tumor size, earlier tumor stage at diagnosis, and a better response to initial therapy [100].
According to 1000 Genomes Project Phase 3, the minor allele (G) frequency of this variant is reported to be “0.008 in Africans, 0.41 in Americans, 0.42 in East Asians, 0.24 in South Asians, and 0.21 in Europe”. A bioinformatic analysis using haploreg4.2 illustrated that this variant was associated with an alternation of six motifs, including AP-4, ASCL2, HEN1, LBP-1, TCF2, VDR theses motifs were related to COPD, inflammation and immunity [101,102,103,104,105,106].
The C to T substitution (c.1100C > T) in exon 13 causes an amino acid shift from alanine to valine at codon 367, resulting in rs879577, another frequent missense variation. This SNP overlaps three IL-17RA transcripts. This variation was associated with CSU; CC genotype and C allele were predominant among patients and associated with poor clinical indicators, concurrent angioedema, leg lesions, and low quality of life [82]. This result may have been caused by substituting the amino acid valine (p.A367V) for alanine (p.A367A) in the IL-17RA cytoplasmic domain. Because the two amino acids have different structural and biochemical prosperities, this substitution may have altered the protein’s conformation and had functional effects.
In the Chinese population, Li et al. found an insignificant relationship between this variation and autoimmune diabetes type 1 [78]. According to Batalla et al., no conclusive evidence exists that rs879577 increases the likelihood of developing psoriasis [80]. The rs879577 variant was highly enriched in Korean alopecia areata patients [107]. The bioinformatics online tool “I-Mutant v.2 (available at https://folding.biofold.org/i-mutant//pages/I-Mutant2.0_Details.html),” accessed on 20 July 2023, shows that this variation results in higher protein stability, and Polyphen-2 v2.2.2, which predicts the functional effect of the variant on proteins, shows that this variant is benign and tolerated.
The minor allele (A) frequency among different ethnicities was found to be “0.47 in Africans, 0.28 in Americans, 0.09 in East Asians, 0.19 in South Asians, and 0.26 in Europe (www.ensembl.org)”, last accessed on 20 July 2023. QLT bioinformatic analysis revealed this polymorphism showed linkage disequilibrium with other 4 SNP on IL17RA rs34484815, rs879575, rs4819555, rs1468488 and was shown to be associated with four altered motifs including AP-2, ER alpha-a, Rad 21, and sin3AK-20. Some have been reported previously to be associated with oxidant responses, COPD, and asthma [108,109,110].
In our study population, this variant was found to reduce disease risk, and the minor allele (T) was shown to decrease risk across all models, similar to rs4819554, the reference allele of rs879577 C was associated with higher disease risk, and T allele was associated with a protective effect, further analysis regarding carrying risk allele and disease severity, outcome, allergies, and aggravating factors, therapeutic response reveal non-significant association.
The third investigated variant in this study was IL17RA_rs41323645G/A. This missense variant overlaps three transcripts, resulting in a different amino acid sequence but where the length is preserved A (Ala) > T (Thr). Ala691Thr. I-Mutant v.2 online tool, this variation results in decreased protein stability, and Polyphen-2 v2.2.2 shows that this variant is predicted to be probably damaging. The minor allele (A) frequency is reported to be 0.24 in Africans, 0.8 in Americans, not present in East Asians, 0.10 in South Asians, and 0.11 in Europe. QLT analysis reveals rs41323645 is in linkage disequilibrium with rs73381911, rs879576, rs57380532, rs11702918, and rs11703539 on IL-17RA and was associated with changes in neuron-restrictive silencer factor (NRSF) motif which shown to exert inhibitory effects on cell proliferation/migration of lung squamous cell carcinoma [111]. To our knowledge, this variation has never been studied in asthma. This mutation, along with variants in the NOD2, EPHA2, and KALRN genes, may play an important role in the emergence of sarcoidosis by sustaining a persistent proinflammatory condition in macrophages [112]. Our findings show that the A allele was more common in patients than controls and was associated with increased disease risk in all models. In contrast, the G allele appears to be protective, and further analysis of the association between the A allele and disease severity, outcome, therapeutic response, and allergy reveals non-significance except for the use of asthma reliever for more than two weeks and risk factors as emotional stress triggers and smoking exposure.
IL17RA_rs4819555C/T (c.2160C > T) (p.Pro720) synonymous variant has a frequency of the minor allele (T) reported to be “0.38 in Africans, 0.27 in Americans, 0.10 in East Asians, 0.18 in South Asians, and 0.25 in Europe (www.ensemble.org),” last accessed on 20 July 2023. In the current study, an eQTL bioinformatic analysis showed that rs4819555 was found to be in linkage disequilibrium with nearby SNPs, such as rs879575 and rs879577, which could also be potential candidates associated with asthma risk. Furthermore, the rs4819555 variant could influence binding with the transcription factors BDP1, EBF, HNF4, NRSF, PLAG1, STAT, and ZNF263, which are linked to carcinogenesis, inflammation, and immunity [113,114,115,116,117]. Although a trend of increased disease risk with the minor (T) allele was observed, there was an insignificant association with disease risk in our study. However, carriers of the T allele significantly correlate with atopic asthma and comorbidities such as allergic rhinitis and eczema. There was also a significant difference in eosinophilic count between carriers and non-carriers.
It is worth noting that the discrepancy between SNPs, either risk or protection in the literature and the current study is caused by the recognized genetic heterogeneity, different populations, ethnicity, various environmental factors, geographic distribution, variations in sample size/study methodology, and the underlying effects that genetic variants cause in a cell/context-specific manner.
Increasing understanding of genetic disease susceptibility, it has been found that considering combined genotypes, or haplotypes, in association studies is far more insightful and compelling for identifying risk related to asthma than relying solely on a single polymorphism [118,119]. The reason for this is the complexity of human genetic composition, in which many diseases, including asthma, present as multifactorial disorders involving more than a single gene variant [120,121]. Consequentially, a risk variant in isolation can provide limited information and be misleading, providing an inadequate understanding of disease susceptibility. Accordingly, in the current study, we found that the most frequent haplotype in the asthma group was ACAC, with a frequency of 22.87% in the patient group vs. 1.36% in controls, and this haplotype was associated with increased asthma risk (Table 5). Also, other haplotypes, including ACGC, ATGC, and GCAC, were significantly associated with asthma, albeit with lower frequencies than ACAC. Although these identified haplotypes were not associated with asthma severity or treatment response, they allow for a more comprehensive and accurate understanding of asthma risk. It embodies the concept of cumulative effect from multiple gene variations. Therefore, the panoramic perspective obtained by studying the interactions of multiple genetic variants confers a more sophisticated understanding of the relationships between gene combinations and disease risk, ultimately leading to enhanced predictive ability and improved personalized medicine strategies.
Although this is the first to highlight the association between IL17RA SNPs rs4819554, rs879577, rs41323645, and rs4819555 and asthma, various limitations should be addressed. The sample was relatively small, and there was an unavoidable selection bias as both patients and controls were enrolled from the same hospital. Functional analyses for the studied variants are required to elucidate the precise mechanism by which these variants confer asthma susceptibility/protection; thus, further research is needed to investigate variant effects in large populations with diverse ethnic backgrounds supported by including more IL17RA-related SNPs through rigorous genetic matching and functional analysis. Furthermore, future studies investigating the potential influence of the identified SNPs on serum biomarkers are highly recommended.

5. Conclusions

In conclusion, the present results indicate a significant association of IL17RA polymorphisms rs4819554, rs879577, and rs41323645 with disease risk under different genetic models and the predominance of the risk haplotype ACAC in the patient group. Meanwhile, the rs4819555 variant showed associations with different phenotypic features and comorbidities of the studied population. These findings need to be validated in larger and more diverse populations.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/children11060657/s1, Table S1: Association of studied gene variants with allergies and aggravating factors in asthma. Table S2: IL17RA gene ontology.

Author Contributions

Conceptualization, S.A.M., N.B.A., E.A.A., E.R., E.A.T., M.S.F. and A.I.; data curation, S.W.K., S.A., W.A., G.A.I. and A.I.; formal analysis, E.A.T.; investigation, N.B.A., E.R., S.A., W.A. and A.I.; methodology, S.A.M., G.A.I. and M.S.F.; resources, N.B.A., E.A.A., E.R., S.W.K., S.A., W.A., G.A.I., G.A.A. and A.I.; software, E.A.T.; supervision, M.S.F.; validation, E.A.A., E.R. and S.W.K.; visualization, E.A.T.; writing—original draft, S.A.M. and E.A.T.; writing—review and editing, N.B.A., E.A.A., E.R., S.W.K., S.A., W.A., G.A.I., G.A.A., M.S.F. and A.I. All authors have read and agreed to the published version of the manuscript.

Funding

The Deanship of Scientific Research at Northern Border University, Arar, KSA, funded this research work through the project number “NBU-FFR-2024–1442-03”.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committee of the Faculty of Medicine, Suez Canal University (approval no. 5101#, 22/11/2022).

Informed Consent Statement

Informed consent was taken from all parents or guardians before participating in the current study.

Data Availability Statement

All generated data in this study are included in the article and Supplement Materials.

Acknowledgments

The authors extend their gratitude to all participants who consented to partake in this study. Further, appreciation is also directed towards the “Oncology Diagnostic Unit” and the “Center of Excellence in Molecular and Cellular Medicine” at Suez Canal University, Egypt. Their valuable provision of facilities made the execution of molecular work possible.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Patel, S.J.; Teach, S.J. Asthma. Pediatr. Rev. 2019, 40, 549–567. [Google Scholar] [CrossRef]
  2. Mortimer, K.; Lesosky, M.; Garcia-Marcos, L.; Asher, M.I.; Pearce, N.; Ellwood, E.; Bissell, K.; El Sony, A.; Ellwood, P.; Marks, G.B.; et al. The burden of asthma, hay fever and eczema in adults in 17 countries: GAN Phase I study. Eur. Respir. J. 2022, 60, 2102865. [Google Scholar] [CrossRef] [PubMed]
  3. Siroux, V.; Bouzigon, E. Asthma heterogeneity: The increasing genetic evidence. Lancet Respir. Med. 2019, 7, 469–471. [Google Scholar] [CrossRef] [PubMed]
  4. Ramakrishnan, R.K.; Al Heialy, S.; Hamid, Q. Role of IL-17 in asthma pathogenesis and its implications for the clinic. Expert Rev. Respir. Med. 2019, 13, 1057–1068. [Google Scholar] [CrossRef] [PubMed]
  5. Gu, C.; Wu, L.; Li, X. IL-17 family: Cytokines, receptors and signaling. Cytokine 2013, 64, 477–485. [Google Scholar] [CrossRef] [PubMed]
  6. Mills, K.H.G. IL-17 and IL-17-producing cells in protection versus pathology. Nat. Rev. Immunol. 2023, 23, 38–54. [Google Scholar] [CrossRef] [PubMed]
  7. Moseley, T.A.; Haudenschild, D.R.; Rose, L.; Reddi, A.H. Interleukin-17 family and IL-17 receptors. Cytokine Growth Factor Rev. 2003, 14, 155–174. [Google Scholar] [CrossRef] [PubMed]
  8. Li, X.; Bechara, R.; Zhao, J.; McGeachy, M.J.; Gaffen, S.L. IL-17 receptor-based signaling and implications for disease. Nat. Immunol. 2019, 20, 1594–1602. [Google Scholar] [CrossRef]
  9. Kanehisa, M.; Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef]
  10. Kanehisa, M.; Sato, Y.; Kawashima, M.; Furumichi, M.; Tanabe, M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 2016, 44, D457–D462. [Google Scholar] [CrossRef]
  11. Chen, K.; Kolls, J.K. Interluekin-17A (IL17A). Gene 2017, 614, 8–14. [Google Scholar] [CrossRef] [PubMed]
  12. Zhang, J.; Wang, D.; Wang, L.; Wang, S.; Roden, A.C.; Zhao, H.; Li, X.; Prakash, Y.S.; Matteson, E.L.; Tschumperlin, D.J.; et al. Profibrotic effect of IL-17A and elevated IL-17RA in idiopathic pulmonary fibrosis and rheumatoid arthritis-associated lung disease support a direct role for IL-17A/IL-17RA in human fibrotic interstitial lung disease. Am. J. Physiol. Lung Cell. Mol. Physiol. 2019, 316, L487–L497. [Google Scholar] [CrossRef] [PubMed]
  13. Amatya, N.; Garg, A.V.; Gaffen, S.L. IL-17 Signaling: The Yin and the Yang. Trends Immunol. 2017, 38, 310–322. [Google Scholar] [CrossRef] [PubMed]
  14. Confalonieri, M.; Salton, F.; Fabiano, F. Acute respiratory distress syndrome. Eur. Respir. Rev. Off. J. Eur. Respir. Soc. 2017, 26, 160116. [Google Scholar] [CrossRef] [PubMed]
  15. Kudo, M.; Melton, A.C.; Chen, C.; Engler, M.B.; Huang, K.E.; Ren, X.; Wang, Y.; Bernstein, X.; Li, J.T.; Atabai, K.; et al. IL-17A produced by αβ T cells drives airway hyper-responsiveness in mice and enhances mouse and human airway smooth muscle contraction. Nat. Med. 2012, 18, 547–554. [Google Scholar] [CrossRef] [PubMed]
  16. Manni, M.L.; Robinson, K.M.; Alcorn, J.F. A tale of two cytokines: IL-17 and IL-22 in asthma and infection. Expert Rev. Respir. Med. 2014, 8, 25–42. [Google Scholar] [CrossRef]
  17. Roos, A.B.; Sethi, S.; Nikota, J.; Wrona, C.T.; Dorrington, M.G.; Sanden, C.; Bauer, C.M.; Shen, P.; Bowdish, D.; Stevenson, C.S.; et al. IL-17A and the Promotion of Neutrophilia in Acute Exacerbation of Chronic Obstructive Pulmonary Disease. Am. J. Respir. Crit. Care Med. 2015, 192, 428–437. [Google Scholar] [CrossRef] [PubMed]
  18. Tiringer, K.; Treis, A.; Fucik, P.; Gona, M.; Gruber, S.; Renner, S.; Dehlink, E.; Nachbaur, E.; Horak, F.; Jaksch, P.; et al. A Th17- and Th2-skewed cytokine profile in cystic fibrosis lungs represents a potential risk factor for Pseudomonas aeruginosa infection. Am. J. Respir. Crit. Care Med. 2013, 187, 621–629. [Google Scholar] [CrossRef]
  19. Tsai, H.C.; Velichko, S.; Hung, L.Y.; Wu, R. IL-17A and Th17 cells in lung inflammation: An update on the role of Th17 cell differentiation and IL-17R signaling in host defense against infection. Clin. Dev. Immunol. 2013, 2013, 267971. [Google Scholar] [CrossRef]
  20. Gurczynski, S.J.; Moore, B.B. IL-17 in the lung: The good, the bad, and the ugly. Am. J. Physiol. Lung Cell. Mol. Physiol. 2018, 314, L6–L16. [Google Scholar] [CrossRef]
  21. Chang, Y.; Al-Alwan, L.; Risse, P.A.; Halayko, A.J.; Martin, J.G.; Baglole, C.J.; Eidelman, D.H.; Hamid, Q. Th17-associated cytokines promote human airway smooth muscle cell proliferation. FASEB J. Off. Publ. Fed. Am. Soc. Exp. Biol. 2012, 26, 5152–5160. [Google Scholar] [CrossRef] [PubMed]
  22. Hynes, G.M.; Hinks, T.S.C. The role of interleukin-17 in asthma: A protective response? ERJ Open Res. 2020, 6, 00364–02019. [Google Scholar] [CrossRef] [PubMed]
  23. Liu, Z.; Niu, C.; Ying, L.; Zhang, Q.; Long, M.; Fu, Z. Exploration of the Serum Interleukin-17 and Interleukin-27 Expression Levels in Children with Bronchial Asthma and Their Correlation with Indicators of Lung Function. Genet. Test. Mol. Biomark. 2020, 24, 10–16. [Google Scholar] [CrossRef] [PubMed]
  24. Liu, D.; Tan, Y.; Bajinka, O.; Wang, L.; Tang, Z. Th17/IL-17 Axis Regulated by Airway Microbes Get Involved in the Development of Asthma. Curr. Allergy Asthma Rep. 2020, 20, 11. [Google Scholar] [CrossRef] [PubMed]
  25. Kousha, A.; Mahdavi Gorabi, A.; Forouzesh, M.; Hosseini, M.; Alexander, M.; Imani, D.; Razi, B.; Mousavi, M.J.; Aslani, S.; Mikaeili, H. Interleukin 4 gene polymorphism (−589C/T) and the risk of asthma: A meta-analysis and met-regression based on 55 studies. BMC Immunol. 2020, 21, 55. [Google Scholar] [CrossRef] [PubMed]
  26. Zhu, L.; Liu, T.; Wang, L.; Li, Q.; Wu, Y.; Liu, B. Polymorphisms in the interleukin 4 promoter −589C/T gene and the risk of asthma: A systematic review and meta-analysis. Transl. Pediatr. 2021, 10, 2355–2365. [Google Scholar] [CrossRef] [PubMed]
  27. Liu, T.; Yin, W.; Luo, L.; Wu, Y.; Qin, S.; Qin, X. Association between Interleukin-4-590C>T Polymorphism and the Susceptibility to Asthma: A Meta-Analysis of Case-Control Study. J. Healthc. Eng. 2022, 2022, 1712715. [Google Scholar] [CrossRef] [PubMed]
  28. Qurashi, T.A.; Bhat, G.A.; Khan, M.S.; Rasool, R.; Sameen, F.; Hassan, I.; Mudassar, S. Interleukin 4 and Interleukin 4 receptor alpha gene variants and risk of atopy—A case control study based assessment. Clin. Immunol. 2021, 229, 108783. [Google Scholar] [CrossRef]
  29. Mahdaviani, S.A.; Eskian, M.; Khorasanizadeh, M.; Bashardoost, B.; Tashayoie Nejad, S.; Jamaati, H.R.; Rezaei, A.; Sadr, M.; Aryan, Z.; Rezaei, N. Interleukin 9 serum level and single nucleotide polymorphism in patients with asthma. Acta Bio-Medica Atenei Parm. 2021, 92, e2021206. [Google Scholar] [CrossRef]
  30. Mocellin, M.; de Azeredo Leitao, L.A.; de Araujo, P.D.; Jones, M.H.; Stein, R.T.; Pitrez, P.M.; de Souza, A.P.D.; Pinto, L.A. Association between interleukin-10 polymorphisms and CD4(+)CD25(+)FOXP3(+) T cells in asthmatic children. J. Pediatr. 2021, 97, 546–551. [Google Scholar] [CrossRef]
  31. Raghdah Maytham, H.; Haidar Abdul Amir Najim, A.; Mohanad Mohsin, A. Interleukin-5 and Interleukin-5 Receptor Polymorphism in Asthma. In Chemokines Updates; Murat, Ş., Ed.; IntechOpen: Rijeka, Croatia, 2022; p. Ch. 7. [Google Scholar]
  32. Hafez, R.A.; Hassan, M.E.; Haggag, M.G.; Atef, N.; Abdallah, A.L.; Gerges, M.A. Association of Interleukin 13 rs20541 Gene Polymorphism and Serum Periostin with Asthma and Allergic Conjunctivitis Among Egyptian Patients. J. Asthma Allergy 2022, 15, 971–982. [Google Scholar] [CrossRef] [PubMed]
  33. Riikonen, R.; Terasjarvi, J.; Lauhkonen, E.; Nuolivirta, K.; He, Q.; Korppi, M. Interleukin 1 receptor-like 1 rs13408661/13431828 polymorphism is associated with persistent post-bronchiolitis asthma at school age. Acta Paediatr. 2022, 111, 628–635. [Google Scholar] [CrossRef] [PubMed]
  34. Wu, M.; Zheng, X.; Huang, J.; Hu, X. Association of IL33, IL1RL1, IL1RAP Polymorphisms and Asthma in Chinese Han Children. Front. Cell Dev. Biol. 2021, 9, 759542. [Google Scholar] [CrossRef]
  35. Ezzat, D.A.; Morgan, D.S.; Mohamed, R.A.; Mohamed, A.F. Genetic association of interleukin 18 (−607C/A, rs1946518) single nucleotide polymorphism with asthmatic children, disease severity and total IgE serum level. Cent.-Eur. J. Immunol. 2019, 44, 285–291. [Google Scholar] [CrossRef]
  36. Konstantinidis, A.K.; Barton, S.J.; Sayers, I.; Yang, I.A.; Lordan, J.L.; Rorke, S.; Clough, J.B.; Holgate, S.T.; Holloway, J.W. Genetic association studies of interleukin-13 receptor α1 subunit gene polymorphisms in asthma and atopy. Eur. Respir. J. 2007, 30, 40–47. [Google Scholar] [CrossRef]
  37. Mosayebian, A.; Ganjalikhani Hakemi, M.; Meshkat, R.; Ghasemi, R.; Khan Ahmad, H.; Samadi, M. Association between Interleukin-23 Receptor R381Q Gene. Iran. J. Allergy Asthma Immunol. 2015, 14, 386–391. [Google Scholar] [PubMed]
  38. Arikan, S.; Ozturk, O.; Duygulu, S.; Atalay, E.O.; Atalay, A. Associations of IL-17 and IL-17 receptor polymorphisms with Behcet’s disease in Denizli Province of Turkey. Immunol. Res. 2023, 71, 600–608. [Google Scholar] [CrossRef]
  39. Cakmak Genc, G.; Karakas Celik, S.; Yilmaz, B.; Piskin, N.; Altinsoy, B.; Dursun, A. Investigation of the relationship between IL17A, IL17F and ILR1N polymorphisms and COVID-19 severity: The predictive role of IL17A rs2275913 polymorphism in the clinical course of COVID-19. Int. J. Immunogenet. 2023, 50, 117–126. [Google Scholar] [CrossRef]
  40. Liberalesso, V.Y.S.W.; Azevedo, M.L.V.; Malaquias, M.A.S.; de Paula, C.B.V.; Nagashima, S.; de Souza, D.G.; Neto, P.C.; Gouveia, K.O.; Biscaro, L.C.; Giamberardino, A.L.G.; et al. The role of IL17 and IL17RA polymorphisms in lethal pandemic acute viral pneumonia (Influenza A virus H1N1 subtype). Surg. Exp. Pathol. 2023, 6, 1. [Google Scholar] [CrossRef]
  41. Ruiz de Morales, J.M.G.; Puig, L.; Dauden, E.; Canete, J.D.; Pablos, J.L.; Martin, A.O.; Juanatey, C.G.; Adan, A.; Montalban, X.; Borruel, N.; et al. Critical role of interleukin (IL)-17 in inflammatory and immune disorders: An updated review of the evidence focusing in controversies. Autoimmun. Rev. 2020, 19, 102429. [Google Scholar] [CrossRef]
  42. Shao, M.; Xu, S.; Yang, H.; Xu, W.; Deng, J.; Chen, Y.; Gao, X.; Guan, S.; Xu, S.; Shuai, Z.; et al. Association between IL-17A and IL-17F gene polymorphism and susceptibility in inflammatory arthritis: A meta-analysis. Clin. Immunol. 2020, 213, 108374. [Google Scholar] [CrossRef]
  43. Jiang, L.; Zhou, X.; Xiong, Y.; Bao, J.; Xu, K.; Wu, L. Association between interleukin-17A/F single nucleotide polymorphisms and susceptibility to osteoarthritis in a Chinese population. Medicine 2019, 98, e14944. [Google Scholar] [CrossRef]
  44. Pasha, H.F.; Tantawy, E.A.; Youssef, M.A. Osteopontin and interleukin-17A genes polymorphisms in Egyptian systemic lupus erythematosus patients: A relation to disease activity and severity. Gene 2019, 702, 107–113. [Google Scholar] [CrossRef] [PubMed]
  45. Lee, Y.H.; Song, G.G. Associations between interleukin 17A and 17F polymorphisms and asthma susceptibility: A meta-analysis. Int. J. Immunogenet. 2023, 50, 53–62. [Google Scholar] [CrossRef]
  46. Zhai, C.; Li, S.; Feng, W.; Shi, W.; Wang, J.; Wang, Q.; Chai, L.; Zhang, Q.; Yan, X.; Li, M. Association of interleukin-17a rs2275913 gene polymorphism and asthma risk: A meta-analysis. Arch. Med. Sci. AMS 2018, 14, 1204–1211. [Google Scholar] [CrossRef] [PubMed]
  47. Ganta, S.; Komaravalli, P.L.; Ahmad, S.; Gaddam, S.L. Influence of genetic variants and mRNA expression of interleukin IL17A gene in asthma susceptibility. Gene 2023, 854, 147119. [Google Scholar] [CrossRef]
  48. Du, J.W.; Xu, Z.L.; Xu, Q.X. Interaction of Interleukin 7 Receptor (IL7R) and IL6 Gene Polymorphisms with Smoking Associated with Susceptibility to Asthma in Chinese Han Adults. Immunol. Investig. 2022, 51, 1364–1371. [Google Scholar] [CrossRef] [PubMed]
  49. Yuce, G.D.; Erdogan, T.; Bozkurt, B.; Toprak, U.; Ceylan, G.G. Relationship between asthma and IL-17 gene polymorphism in a Turkish population. Ir. J. Med. Sci. 2023, 192, 269–275. [Google Scholar] [CrossRef]
  50. Huang, H.T.; Lu, Y.L.; Wang, R.; Qin, H.M.; Wang, C.F.; Wang, J.L.; Xiang, Y.; Guo, J.; Lan, Y.; Wei, Y.S. The association of IL-17A polymorphisms with IL-17A serum levels and risk of ischemic stroke. Oncotarget 2017, 8, 103499–103508. [Google Scholar] [CrossRef]
  51. El-Qutob, D.; Maillo, M. New GINA Guidelines: Controversy Still Exists. Clin. Pulm. Med. 2020, 27, 61–63. [Google Scholar] [CrossRef]
  52. Cole, T.J.; Bellizzi, M.C.; Flegal, K.M.; Dietz, W.H. Establishing a standard definition for child overweight and obesity worldwide: International survey. Bmj 2000, 320, 1240–1243. [Google Scholar] [CrossRef]
  53. Marshall, W.A.; Tanner, J.M. Variations in pattern of pubertal changes in girls. Arch. Dis. Child. 1969, 44, 291–303. [Google Scholar] [CrossRef] [PubMed]
  54. Graham, B.L.; Steenbruggen, I.; Miller, M.R.; Barjaktarevic, I.Z.; Cooper, B.G.; Hall, G.L.; Hallstrand, T.S.; Kaminsky, D.A.; McCarthy, K.; McCormack, M.C.; et al. Standardization of Spirometry 2019 Update. An Official American Thoracic Society and European Respiratory Society Technical Statement. Am. J. Respir. Crit. Care Med. 2019, 200, e70–e88. [Google Scholar] [CrossRef] [PubMed]
  55. Hussein, M.H.; Sobhy, K.E.; Sabry, I.M.; El Serafi, A.T.; Toraih, E.A. Beta2-adrenergic receptor gene haplotypes and bronchodilator response in Egyptian patients with chronic obstructive pulmonary disease. Adv. Med. Sci. 2017, 62, 193–201. [Google Scholar] [CrossRef] [PubMed]
  56. Coates, A.L.; Wanger, J.; Cockcroft, D.W.; Culver, B.H.; Bronchoprovocation Testing Task Force: Kai-Hakon, C.; Diamant, Z.; Gauvreau, G.; Hall, G.L.; Hallstrand, T.S.; Horvath, I.; et al. ERS technical standard on bronchial challenge testing: General considerations and performance of methacholine challenge tests. Eur. Respir. J. 2017, 49, 1601526. [Google Scholar] [CrossRef] [PubMed]
  57. Crapo, R.O.; Casaburi, R.; Coates, A.L.; Enright, P.L.; Hankinson, J.L.; Irvin, C.G.; MacIntyre, N.R.; McKay, R.T.; Wanger, J.S.; Anderson, S.D.; et al. Guidelines for methacholine and exercise challenge testing-1999. This official statement of the American Thoracic Society was adopted by the ATS Board of Directors, July 1999. Am. J. Respir. Crit. Care Med. 2000, 161, 309–329. [Google Scholar] [CrossRef]
  58. Price, D.B.; Rigazio, A.; Campbell, J.D.; Bleecker, E.R.; Corrigan, C.J.; Thomas, M.; Wenzel, S.E.; Wilson, A.M.; Small, M.B.; Gopalan, G.; et al. Blood eosinophil count and prospective annual asthma disease burden: A UK cohort study. Lancet Respir. Med. 2015, 3, 849–858. [Google Scholar] [CrossRef]
  59. Toraih, E.A.; Hussein, M.H.; Al Ageeli, E.; Riad, E.; AbdAllah, N.B.; Helal, G.M.; Fawzy, M.S. Structure and functional impact of seed region variant in MIR-499 gene family in bronchial asthma. Respir. Res. 2017, 18, 169. [Google Scholar] [CrossRef] [PubMed]
  60. Toraih, E.A.; Hussein, M.H.; Ibrahim, A.; AbdAllah, N.B.; Mohammad, E.; Kishk, A.M.; Fawzy, M.S. Beta2-adrenergic receptor variants in children and adolescents with bronchial asthma. Front. Biosci. 2019, 11, 61–78. [Google Scholar] [CrossRef]
  61. Omasits, U.; Ahrens, C.H.; Muller, S.; Wollscheid, B. Protter: Interactive protein feature visualization and integration with experimental proteomic data. Bioinformatics 2014, 30, 884–886. [Google Scholar] [CrossRef]
  62. Steinhaus, R.; Proft, S.; Schuelke, M.; Cooper, D.N.; Schwarz, J.M.; Seelow, D. MutationTaster2021. Nucleic Acids Res. 2021, 49, W446–W451. [Google Scholar] [CrossRef] [PubMed]
  63. Sim, N.L.; Kumar, P.; Hu, J.; Henikoff, S.; Schneider, G.; Ng, P.C. SIFT web server: Predicting effects of amino acid substitutions on proteins. Nucleic Acids Res. 2012, 40, W452–W457. [Google Scholar] [CrossRef] [PubMed]
  64. Thomas, P.D.; Ebert, D.; Muruganujan, A.; Mushayahama, T.; Albou, L.P.; Mi, H. PANTHER: Making genome-scale phylogenetics accessible to all. Protein Sci. A Publ. Protein Soc. 2022, 31, 8–22. [Google Scholar] [CrossRef] [PubMed]
  65. Kopanos, C.; Tsiolkas, V.; Kouris, A.; Chapple, C.E.; Albarca Aguilera, M.; Meyer, R.; Massouras, A. VarSome: The human genomic variant search engine. Bioinformatics 2019, 35, 1978–1980. [Google Scholar] [CrossRef] [PubMed]
  66. Capriotti, E.; Fariselli, P.; Casadio, R. I-Mutant2.0: Predicting stability changes upon mutation from the protein sequence or structure. Nucleic Acids Res. 2005, 33, W306–W310. [Google Scholar] [CrossRef] [PubMed]
  67. Ward, L.D.; Kellis, M. HaploReg: A resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants. Nucleic Acids Res. 2012, 40, D930–D934. [Google Scholar] [CrossRef] [PubMed]
  68. Gaffen, S.L. Structure and signalling in the IL-17 receptor family. Nat. Rev. Immunol. 2009, 9, 556–567. [Google Scholar] [CrossRef] [PubMed]
  69. Lee, E.; Lee, S.Y.; Park, M.J.; Hong, S.J. TNF-α (rs1800629) polymorphism modifies the effect of sensitization to house dust mite on asthma and bronchial hyperresponsiveness in children. Exp. Mol. Pathol. 2020, 115, 104467. [Google Scholar] [CrossRef]
  70. El-Shal, A.S.; Shalaby, S.M.; Abdel-Nour, H.M.; Sarhan, W.M.; Hamed Gehad, M.; Mohamed Yousif, Y. Impact of cytokines genes polymorphisms and their serum levels on childhood asthma in Egyptian population. Cytokine 2022, 157, 155933. [Google Scholar] [CrossRef]
  71. Jin, X.; Zheng, J. IL-4-C-590T locus polymorphism and susceptibility to asthma in children: A meta-analysis. J. Pediatr. 2021, 97, 264–272. [Google Scholar] [CrossRef]
  72. Elsaid, A.; Shoaib, R.M.S.; Badr, S.S.; Wahba, Y.; Ayyad, S.N. Polymorphisms of interleukin 4 and interleukin 4 receptor genes and bronchial asthma risk among Egyptian children. Clin. Biochem. 2021, 93, 66–72. [Google Scholar] [CrossRef] [PubMed]
  73. Demenais, F.; Margaritte-Jeannin, P.; Barnes, K.C.; Cookson, W.O.C.; Altmuller, J.; Ang, W.; Barr, R.G.; Beaty, T.H.; Becker, A.B.; Beilby, J.; et al. Multiancestry association study identifies new asthma risk loci that colocalize with immune-cell enhancer marks. Nat. Genet. 2018, 50, 42–53. [Google Scholar] [CrossRef] [PubMed]
  74. Schoettler, N.; Strek, M.E. Recent Advances in Severe Asthma: From Phenotypes to Personalized Medicine. Chest 2020, 157, 516–528. [Google Scholar] [CrossRef] [PubMed]
  75. Zhu, Z.; Lee, P.H.; Chaffin, M.D.; Chung, W.; Loh, P.R.; Lu, Q.; Christiani, D.C.; Liang, L. A genome-wide cross-trait analysis from UK Biobank highlights the shared genetic architecture of asthma and allergic diseases. Nat. Genet. 2018, 50, 857–864. [Google Scholar] [CrossRef] [PubMed]
  76. Laulajainen-Hongisto, A.; Lyly, A.; Hanif, T.; Dhaygude, K.; Kankainen, M.; Renkonen, R.; Donner, K.; Mattila, P.; Jartti, T.; Bousquet, J.; et al. Genomics of asthma, allergy and chronic rhinosinusitis: Novel concepts and relevance in airway mucosa. Clin. Transl. Allergy 2020, 10, 45. [Google Scholar] [CrossRef] [PubMed]
  77. Thompson, E.E.; Dang, Q.; Mitchell-Handley, B.; Rajendran, K.; Ram-Mohan, S.; Solway, J.; Ober, C.; Krishnan, R. Cytokine-induced molecular responses in airway smooth muscle cells inform genome-wide association studies of asthma. Genome Med. 2020, 12, 64. [Google Scholar] [CrossRef]
  78. Li, J.; Xu, L.; Zhao, W.; Pan, J.; Lu, J.; Lu, H.; Yan, J.; Weng, J.; Liu, F. Serum IL-17A concentration and a IL17RA single nucleotide polymorphism contribute to the risk of autoimmune type 1 diabetes. Diabetes/Metab. Res. Rev. 2022, 38, e3547. [Google Scholar] [CrossRef]
  79. Sabry, D.; Aboraia, N.; Samir, M. A potential association between psoriasin to rs4819554 of IL-17RA gene polymorphism in psoriasis Egyptian patients. Arch. Dermatol. Res. 2020, 312, 273–281. [Google Scholar] [CrossRef]
  80. Batalla, A.; Coto, E.; Gonzalez-Lara, L.; Gonzalez-Fernandez, D.; Gomez, J.; Aranguren, T.F.; Queiro, R.; Santos-Juanes, J.; Lopez-Larrea, C.; Coto-Segura, P. Association between single nucleotide polymorphisms IL17RA rs4819554 and IL17E rs79877597 and Psoriasis in a Spanish cohort. J. Dermatol. Sci. 2015, 80, 111–115. [Google Scholar] [CrossRef]
  81. Batalla, A.; Coto, E.; Gomez, J.; Eiris, N.; Gonzalez-Fernandez, D.; Gomez-De Castro, C.; Dauden, E.; Llamas-Velasco, M.; Prieto-Perez, R.; Abad-Santos, F.; et al. IL17RA gene variants and anti-TNF response among psoriasis patients. Pharmacogenomics J. 2018, 18, 76–80. [Google Scholar] [CrossRef]
  82. Nada, H.; Hassan, R.; Ibrahim, R.A.E.; Abdelsalam, O.E.; Fathy, A.; Toraih, E.A.; Atwa, M.A. Interleukin 17 receptor A haplotype analysis in chronic spontaneous urticaria: A preliminary study. J. Cosmet. Dermatol. 2021, 20, 1331–1342. [Google Scholar] [CrossRef] [PubMed]
  83. Sanad, E.M.; Salem, R.M.; Ali, H.M.; Elshafee, S.N. IL17RA Gene polymorphism in Alopecia Areata in an Egyption sample. Benha J. Appl. Sci. 2020, 5, 281–284. [Google Scholar] [CrossRef]
  84. Vidal-Castineira, J.R.; Lopez-Vazquez, A.; Diaz-Pena, R.; Diaz-Bulnes, P.; Martinez-Camblor, P.; Coto, E.; Coto-Segura, P.; Bruges-Armas, J.; Pinto, J.A.; Blanco, F.J.; et al. A Single Nucleotide Polymorphism in the Il17ra Promoter Is Associated with Functional Severity of Ankylosing Spondylitis. PLoS ONE 2016, 11, e0158905. [Google Scholar] [CrossRef] [PubMed]
  85. Lappalainen, T.; MacArthur, D.G. From variant to function in human disease genetics. Science 2021, 373, 1464–1468. [Google Scholar] [CrossRef] [PubMed]
  86. Somers, J.; Ruttens, D.; Verleden, S.E.; Vandermeulen, E.; Piloni, D.; Wauters, E.; Lambrechts, D.; Vos, R.; Verleden, G.M.; Vanaudenaerde, B.; et al. Interleukin-17 receptor polymorphism predisposes to primary graft dysfunction after lung transplantation. J. Heart Lung Transplant. Off. Publ. Int. Soc. Heart Transplant. 2015, 34, 941–949. [Google Scholar] [CrossRef] [PubMed]
  87. Fathollahpour, A.; Abyaneh, F.A.; Darabi, B.; Ebrahimi, M.; Kooti, W.; Kalmarzi, R.N. Aspirin-Exacerbated Respiratory Disease Polymorphisms; a review study. Gene 2023, 870, 147326. [Google Scholar] [CrossRef] [PubMed]
  88. Park, J.S.; Park, B.L.; Kim, M.O.; Heo, J.S.; Jung, J.S.; Bae, D.J.; Uh, S.T.; Kim, M.K.; Choi, I.S.; Cho, S.H.; et al. Association of single nucleotide polymorphisms on Interleukin 17 receptor A (IL17RA) gene with aspirin hypersensitivity in asthmatics. Hum. Immunol. 2013, 74, 598–606. [Google Scholar] [CrossRef] [PubMed]
  89. Johnston, A.; Guzman, A.M.; Swindell, W.R.; Wang, F.; Kang, S.; Gudjonsson, J.E. Early tissue responses in psoriasis to the antitumour necrosis factor-α biologic etanercept suggest reduced interleukin-17 receptor expression and signalling. Br. J. Dermatol. 2014, 171, 97–107. [Google Scholar] [CrossRef] [PubMed]
  90. Wong, L.Y.; Hatfield, J.K.; Brown, M.A. Ikaros sets the potential for Th17 lineage gene expression through effects on chromatin state in early T cell development. J. Biol. Chem. 2013, 288, 35170–35179. [Google Scholar] [CrossRef]
  91. Lauhkonen, E.; Holster, A.; Terasjarvi, J.; Nuolivirta, K.; He, Q.; Korppi, M. IL17RA variations showed no associations with post-bronchiolitis asthma or lung function. Pediatr. Int. Off. J. Jpn. Pediatr. Soc. 2021, 63, 196–201. [Google Scholar] [CrossRef]
  92. Murdaca, G.; Negrini, S.; Magnani, O.; Penza, E.; Pellecchio, M.; Puppo, F. Impact of pharmacogenomics upon the therapeutic response to etanercept in psoriasis and psoriatic arthritis. Expert Opin. Drug Saf. 2017, 16, 1173–1179. [Google Scholar] [CrossRef] [PubMed]
  93. Sandip, C.; Tan, L.; Huang, J.; Li, Q.; Ni, L.; Cianflone, K.; Wang, D.W. Common variants in IL-17A/IL-17RA axis contribute to predisposition to and progression of congestive heart failure. Medicine 2016, 95, e4105. [Google Scholar] [CrossRef] [PubMed]
  94. Coto, E.; Pascual, I.; Avanzas, P.; Cuesta-Lavona, E.; Lorca, R.; Martin, M.; Vazquez-Coto, D.; Diaz-Corte, C.; Moris, C.; Rodriguez-Reguero, J.; et al. IL17RA in early-onset coronary artery disease: Total leukocyte transcript analysis and promoter polymorphism (rs4819554) association. Cytokine 2020, 136, 155285. [Google Scholar] [CrossRef] [PubMed]
  95. Chen, A.; Zhao, H.; Wang, J.; Zhang, R.; Liu, J.; Zhao, X.; Li, C.; Jia, X.; Li, X.; Lin, Y.; et al. Haplotype Analysis of Candidate Genes Involved in Inflammation and Oxidative Stress and the Susceptibility to Preeclampsia. J. Immunol. Res. 2020, 2020, 4683798. [Google Scholar] [CrossRef] [PubMed]
  96. Cai, T.; Wang, G.; Yang, Y.; Mu, K.; Zhang, J.; Jiang, Y.; Zhang, J.A. Association Between Polymorphisms of IL-23/IL-17 Pathway and Clinical Phenotypes of Autoimmune Thyroid Diseases. Iran. J. Immunol. IJI 2022, 19, 139–149. [Google Scholar] [CrossRef] [PubMed]
  97. Wielinska, J.; Swierkot, J.; Kolossa, K.; Bugaj, B.; Chaszczewska-Markowska, M.; Jeka, S.; Bogunia-Kubik, K. Polymorphisms within Genes Coding for IL-17A and F and Their Receptor as Clinical Hallmarks in Ankylosing Spondylitis. Mediat. Inflamm. 2021, 2021, 3125922. [Google Scholar] [CrossRef] [PubMed]
  98. Coto, E.; Gomez, J.; Suarez, B.; Tranche, S.; Diaz-Corte, C.; Ortiz, A.; Ruiz-Ortega, M.; Coto-Segura, P.; Batalla, A.; Lopez-Larrea, C. Association between the IL17RA rs4819554 polymorphism and reduced renal filtration rate in the Spanish RENASTUR cohort. Hum. Immunol. 2015, 76, 75–78. [Google Scholar] [CrossRef]
  99. Zhang, X.; Ye, T.; Li, M.; Yan, H.; Lin, H.; Lu, H.; Qi, Z.; Sheng, H.; He, C. Association of Polymorphisms in Inflammation Genes with the Prognosis of Advanced Non-Small Cell Lung Cancer Patients Receiving Epidermal Growth Factor Receptor Tyrosine Kinase Inhibitors. Front. Oncol. 2022, 12, 836117. [Google Scholar] [CrossRef] [PubMed]
  100. Bertol, B.C.; Goes de Araujo, J.N.; de Carvalho, K.T.C.; Gonzaga Dos Santos, M.; Martelli-Palomino, G.; de Figueiredo Feitosa, N.L.; Maciel, L.M.Z.; Silbiger, V.N.; Lucena-Silva, N.; Freitas, J.; et al. Polymorphisms at the IL17A and IL17RA Genes are Associated with Prognosis of Papillary Thyroid Carcinoma. Arch. Med. Res. 2022, 53, 163–169. [Google Scholar] [CrossRef]
  101. Huang, Y.; Zhou, F.; Shen, C.; Wang, H.; Xiao, Y. LBP reduces theinflammatory injuryof kidney in septic rat and regulates the Keap1-Nrf2/ARE signaling pathway1. Acta Cir. Bras. 2019, 34, e20190010000003. [Google Scholar] [CrossRef]
  102. Kavanagh, K.; Hsu, F.C.; Davis, A.T.; Kritchevsky, S.B.; Rejeski, W.J.; Kim, S. Biomarkers of leaky gut are related to inflammation and reduced physical function in older adults with cardiometabolic disease and mobility limitations. GeroScience 2019, 41, 923–933. [Google Scholar] [CrossRef]
  103. Barshad, G.; Webb, L.M.; Ting, H.A.; Oyesola, O.O.; Onyekwere, O.G.; Lewis, J.J.; Rice, E.J.; Matheson, M.K.; Sun, X.H.; von Moltke, J.; et al. E-Protein Inhibition in ILC2 Development Shapes the Function of Mature ILC2s during Allergic Airway Inflammation. J. Immunol. 2022, 208, 1007–1020. [Google Scholar] [CrossRef]
  104. Huo, R.; Tian, X.; Chang, Q.; Liu, D.; Wang, C.; Bai, J.; Wang, R.; Zheng, G.; Tian, X. Targeted inhibition of β-catenin alleviates airway inflammation and remodeling in asthma via modulating the profibrotic and anti-inflammatory actions of transforming growth factor-β1. Ther. Adv. Respir. Dis. 2021, 15, 1753466620981858. [Google Scholar] [CrossRef] [PubMed]
  105. Hu, G.; Dong, T.; Wang, S.; Jing, H.; Chen, J. Vitamin D(3)-vitamin D receptor axis suppresses pulmonary emphysema by maintaining alveolar macrophage homeostasis and function. EBioMedicine 2019, 45, 563–577. [Google Scholar] [CrossRef] [PubMed]
  106. Aggeletopoulou, I.; Thomopoulos, K.; Mouzaki, A.; Triantos, C. Vitamin D-VDR Novel Anti-Inflammatory Molecules-New Insights into Their Effects on Liver Diseases. Int. J. Mol. Sci. 2022, 23, 8465. [Google Scholar] [CrossRef] [PubMed]
  107. Lew, B.L.; Cho, H.R.; Haw, S.; Kim, H.J.; Chung, J.H.; Sim, W.Y. Association between IL17A/IL17RA Gene Polymorphisms and Susceptibility to Alopecia Areata in the Korean Population. Ann. Dermatol. 2012, 24, 61–65. [Google Scholar] [CrossRef]
  108. Aneas, I.; Decker, D.C.; Howard, C.L.; Sobreira, D.R.; Sakabe, N.J.; Blaine, K.M.; Stein, M.M.; Hrusch, C.L.; Montefiori, L.E.; Tena, J.; et al. Asthma-associated genetic variants induce IL33 differential expression through an enhancer-blocking regulatory region. Nat. Commun. 2021, 12, 6115. [Google Scholar] [CrossRef]
  109. Pierrou, S.; Broberg, P.; O’Donnell, R.A.; Pawlowski, K.; Virtala, R.; Lindqvist, E.; Richter, A.; Wilson, S.J.; Angco, G.; Moller, S.; et al. Expression of genes involved in oxidative stress responses in airway epithelial cells of smokers with chronic obstructive pulmonary disease. Am. J. Respir. Crit. Care Med. 2007, 175, 577–586. [Google Scholar] [CrossRef]
  110. Costanzo, L.; Soto, B.; Meier, R.; Geraghty, P. The Biology and Function of Tissue Inhibitor of Metalloproteinase 2 in the Lungs. Pulm. Med. 2022, 2022, 3632764. [Google Scholar] [CrossRef]
  111. Wang, K.; Xu, K.; Leng, X.; Han, Y.; Fang, Q. miRNA-9 Inhibits Proliferation and Migration of Lung Squamous Cell Carcinoma Cells by Regulating NRSF/EGFR. Technol. Cancer Res. Treat. 2020, 19, 1533033820945807. [Google Scholar] [CrossRef]
  112. Besnard, V.; Calender, A.; Bouvry, D.; Pacheco, Y.; Chapelon-Abric, C.; Jeny, F.; Nunes, H.; Planes, C.; Valeyre, D. G908R NOD2 variant in a family with sarcoidosis. Respir. Res. 2018, 19, 44. [Google Scholar] [CrossRef] [PubMed]
  113. Collins, P.E.; Grassia, G.; Colleran, A.; Kiely, P.A.; Ialenti, A.; Maffia, P.; Carmody, R.J. Mapping the Interaction of B Cell Leukemia 3 (BCL-3) and Nuclear Factor κB (NF-κB) p50 Identifies a BCL-3-mimetic Anti-inflammatory Peptide. J. Biol. Chem. 2015, 290, 15687–15696. [Google Scholar] [CrossRef] [PubMed]
  114. Cabarcas-Petroski, S.; Schramm, L. BDP1 Alterations Correlate with Clinical Outcomes in Breast Cancer. Cancers 2022, 14, 1658. [Google Scholar] [CrossRef] [PubMed]
  115. Alunno, A.; Padjen, I.; Fanouriakis, A.; Boumpas, D.T. Pathogenic and Therapeutic Relevance of JAK/STAT Signaling in Systemic Lupus Erythematosus: Integration of Distinct Inflammatory Pathways and the Prospect of Their Inhibition with an Oral Agent. Cells 2019, 8, 898. [Google Scholar] [CrossRef] [PubMed]
  116. Tuli, H.S.; Sak, K.; Iqubal, A.; Garg, V.K.; Varol, M.; Sharma, U.; Chauhan, A.; Yerer, M.B.; Dhama, K.; Jain, M.; et al. STAT signaling as a target for intervention: From cancer inflammation and angiogenesis to non-coding RNAs modulation. Mol. Biol. Rep. 2022, 49, 8987–8999. [Google Scholar] [CrossRef] [PubMed]
  117. Sun, T.; Cao, Y.; Huang, T.; Sang, Y.; Dai, Y.; Tao, Z. Comprehensive analysis of fifteen hub genes to identify a promising diagnostic model, regulated networks, and immune cell infiltration in acute kidney injury. J. Clin. Lab. Anal. 2022, 36, e24709. [Google Scholar] [CrossRef] [PubMed]
  118. Shi, M.; Umbach, D.M.; Weinberg, C.R. Identification of risk-related haplotypes with the use of multiple SNPs from nuclear families. Am. J. Hum. Genet. 2007, 81, 53–66. [Google Scholar] [CrossRef] [PubMed]
  119. Vishweswaraiah, S.; Ramachandra, N.B.; Jayaraj, B.S.; Holla, A.D.; Chakraborty, S.; Agrawal, A.; Mahesh, P.A. Haplotype analysis of ADAM33 polymorphisms in asthma: A pilot study. Indian J. Med. Res. 2019, 150, 272–281. [Google Scholar] [CrossRef]
  120. Thomsen, S.F. Genetics of asthma: An introduction for the clinician. Eur. Clin. Respir. J. 2015, 2, 24643. [Google Scholar] [CrossRef]
  121. Ahmadi, K.R.; Goldstein, D.B. Multifactorial diseases: Asthma genetics point the way. Curr. Biol. CB 2002, 12, R702–R704. [Google Scholar] [CrossRef]
Figure 1. Interleukin-17 signaling pathway (Kyoto Encyclopedia of Genes and Genomes (KEGG) ID: Human (Homo sapiens), hsa04657). The IL-17 cytokine family (IL-17A-F) plays a role in acute and chronic inflammation. IL-17A, predominantly secreted by T helper 17 (TH17) cells, is essential for protecting against external pathogens and is implicated in autoimmune inflammation. IL-17F primarily protects mucosal surfaces, while IL-17E promotes Th2 immune responses. IL-17C shares similar functions with IL-17A. The IL-17 family communicates via specific receptors, such as IL17RA, which several family members use. Activation of downstream signaling pathways, including NF-kappa B, mitogen-activated protein kinases (MAPKs), and CCAAT-enhancer-binding proteins (C/EBPs), results in the production of antimicrobial peptides and the expression of cytokines and chemokines. Adopted from (https://www.genome.jp/pathway/hsa04657+23765) (last acceded 7 July 2023) [9,10].
Figure 1. Interleukin-17 signaling pathway (Kyoto Encyclopedia of Genes and Genomes (KEGG) ID: Human (Homo sapiens), hsa04657). The IL-17 cytokine family (IL-17A-F) plays a role in acute and chronic inflammation. IL-17A, predominantly secreted by T helper 17 (TH17) cells, is essential for protecting against external pathogens and is implicated in autoimmune inflammation. IL-17F primarily protects mucosal surfaces, while IL-17E promotes Th2 immune responses. IL-17C shares similar functions with IL-17A. The IL-17 family communicates via specific receptors, such as IL17RA, which several family members use. Activation of downstream signaling pathways, including NF-kappa B, mitogen-activated protein kinases (MAPKs), and CCAAT-enhancer-binding proteins (C/EBPs), results in the production of antimicrobial peptides and the expression of cytokines and chemokines. Adopted from (https://www.genome.jp/pathway/hsa04657+23765) (last acceded 7 July 2023) [9,10].
Children 11 00657 g001
Figure 2. Molecular features of interleukin-17RA (IL17RA) gene and protein. (A) The IL17RA gene (Gene ID: 23765) is located on the long arm of chromosome 22-NC_000022.11:17,084,954-17,115,693 according to Homo sapiens assembly GRCh38.p14, with 30,740 bp long, on the forward strand. It consists of 13 exons. (B) The gene has eight transcripts: the major one (IL17RA-201 transcript; ENST00000319363.11) has 8566 nucleotides, which codes for 866aa protein; the IL17RA-204 (2499 nucleotides) codes for the 832aa protein, the IL17RA-207 (739 nucleotides) codes for the129aa nonsense-mediated decay transcript, the IL17RA-203 (651 nucleotides) codes for the 59aa nonsense-mediated decay transcript, and the IL17RA-206, 205, 208, and 202 (675, 611, 847, and 621 nucleotides, respectively) are non-protein coding transcripts. (C) The encoded protein is a single polypeptide chain of 866 amino acid residues, starting with the signal peptide (1–32 aa, red color) in the amino terminus, the extracellular domain (33–320aa), the transmembrane domain (321–341aa), and the cytoplasmic domain (342–866aa). The studied missense variants-related amino acid substitutions were labeled in the cytoplasmic domain and explained in the legend with the specified colors. (D) The subcellular localization of IL17RA can be observed, with the intensity of green coloration indicating its relative abundance. (E) The interacting protein–protein network demonstrates the functional partners of the IL17RA. The generated network consists of 11 nodes (each node corresponds to a protein generated from a unique gene locus) and 45 edges (each edge signifies an association between proteins, which could imply sharing functionality but does not necessarily mean the proteins are physically bound to each other), with average local clustering coefficient = 0.889, and the protein–protein interaction enrichment p-value = 4.22 × 10−15. (F) IL17RA gene coexpression triangle matrix. The confidence level of association is demonstrated by the color intensity, given the overall expression data in humans. Data sources: “https://www.genecards.org/, https://www.ncbi.nlm.nih.gov/gene/23765, https://asia.ensembl.org/Homo_sapiens/Tran-script/, https://wlab.ethz.ch/protter/, https://www.uniprot.org/, https://compartments.jensenlab.org/, and https://string-db.org/” (last accessed on 20 July 2023).
Figure 2. Molecular features of interleukin-17RA (IL17RA) gene and protein. (A) The IL17RA gene (Gene ID: 23765) is located on the long arm of chromosome 22-NC_000022.11:17,084,954-17,115,693 according to Homo sapiens assembly GRCh38.p14, with 30,740 bp long, on the forward strand. It consists of 13 exons. (B) The gene has eight transcripts: the major one (IL17RA-201 transcript; ENST00000319363.11) has 8566 nucleotides, which codes for 866aa protein; the IL17RA-204 (2499 nucleotides) codes for the 832aa protein, the IL17RA-207 (739 nucleotides) codes for the129aa nonsense-mediated decay transcript, the IL17RA-203 (651 nucleotides) codes for the 59aa nonsense-mediated decay transcript, and the IL17RA-206, 205, 208, and 202 (675, 611, 847, and 621 nucleotides, respectively) are non-protein coding transcripts. (C) The encoded protein is a single polypeptide chain of 866 amino acid residues, starting with the signal peptide (1–32 aa, red color) in the amino terminus, the extracellular domain (33–320aa), the transmembrane domain (321–341aa), and the cytoplasmic domain (342–866aa). The studied missense variants-related amino acid substitutions were labeled in the cytoplasmic domain and explained in the legend with the specified colors. (D) The subcellular localization of IL17RA can be observed, with the intensity of green coloration indicating its relative abundance. (E) The interacting protein–protein network demonstrates the functional partners of the IL17RA. The generated network consists of 11 nodes (each node corresponds to a protein generated from a unique gene locus) and 45 edges (each edge signifies an association between proteins, which could imply sharing functionality but does not necessarily mean the proteins are physically bound to each other), with average local clustering coefficient = 0.889, and the protein–protein interaction enrichment p-value = 4.22 × 10−15. (F) IL17RA gene coexpression triangle matrix. The confidence level of association is demonstrated by the color intensity, given the overall expression data in humans. Data sources: “https://www.genecards.org/, https://www.ncbi.nlm.nih.gov/gene/23765, https://asia.ensembl.org/Homo_sapiens/Tran-script/, https://wlab.ethz.ch/protter/, https://www.uniprot.org/, https://compartments.jensenlab.org/, and https://string-db.org/” (last accessed on 20 July 2023).
Children 11 00657 g002
Table 1. TaqMan assay-related data applied in the present study.
Table 1. TaqMan assay-related data applied in the present study.
SNP Reference SequenceAssay IDLocation on Ch.22 According to GRCh38 AssemblyContext Sequence [VIC/FAM]
rs4819554C__337392_3017084145GGGAAGTAACGACTCTCTTAGGTGC[A/G]GC TGGGACACAGTCTCACAGACCAG
rs879577C__2666446_20 17108319 CTTGTTTCCTTAGATGGCCTGCCTG[C/T]GGC TGACCTGATCCCCCCACCGCTG
rs41323645 C__86401686_10 17109290 GCCTGGGCCCCTGGCTGACGGTGCC[A/G]CA GTCCGGCTGGCACTGGCGGGGG
rs4819555 C__28000962_10 17109379 GGCGAAATAGCGTCCTCTTCCTCCC[C/T]GTG GACCCCGAGGACTCGCCCCTTG
All variants are of the “transition substitution” type. All assays have the catalog number 4351379 (Applied Biosystems, ThermoFisher Scientific, Waltham, MA, USA) (https://www.thermofisher.com/), last accessed on 20 July 2023. The ancestral/minor alleles of each studied variant labeled by VIC/FAM probes are blue-colored in the context gene sequence.
Table 2. Baseline characteristics of the study population.
Table 2. Baseline characteristics of the study population.
CharacteristicsLevelsControl (n = 96)Asthma (n = 96)p-Value
AgeMean ± SD9.6 ± 3.19.7 ± 3.00.76
SexFemale56 (58.3%)52 (54.2%)0.66
Male40 (41.7%)44 (45.8%)
ResidencyRural 51 (53.1%)64 (66.7%)0.08
Urban45 (46.9%)32 (33.3%)
FH of asthmaNegative 89 (92.7%)67 (69.8%)<0.001
Positive7 (7.3%)29 (30.2%)
BMI %<85th51 (53.1%)53 (55.2%)0.95
<95th34 (35.4%)32 (33.3%)
≥95th11 (11.5%)11 (11.5%)
Pubertal statusNegative52 (54.2%)51 (53.1%)0.88
Positive44 (45.8%)45 (46.9%)
Tanner stage *Stage 152 (54.2%)51 (53.1%)0.65
Stage 214 (14.6%)17 (17.7%)
Stage 318 (18.8%)12 (12.5%)
Stage 46 (6.3%)10 (10.4%)
Stage 56 (6.3%)6 (6.3%)
Data (except age) are reported as frequency/percentage. A Chi-Square test (two-sided) was employed. The bold value indicates p-value < 0.05. * According to “Marshall and Tanner” [53]. FH: family history, BMI%: body mass index percentile.
Table 3. Genotype frequency and genetic association model.
Table 3. Genotype frequency and genetic association model.
Genetic ModelsGenotypeControls (n= 96)Patients (n = 96)OR (95%CI)p-ValueAIC
rs4819554
Codominant ModelA/A29 (30.2%)47 (49%)10.001250.4
A/G44 (45.8%)42 (43.8%)0.69 (0.35–1.36)
G/G23 (24%)7 (7.3%)0.15 (0.05–0.45)
Dominant ModelA/A29 (30.2%)47 (49%)10.028257.4
A/G-G/G67 (69.8%)49 (51%)0.49 (0.26–0.93)
Recessive ModelA/A-A/G73 (76%)89 (92.7%)1<0.001249.6
G/G23 (24%)7 (7.3%)0.18 (0.07–0.52)
Over dominant A/A-G/G52 (54.2%)54 (56.2%)10.73262
A/G44 (45.8%)42 (43.8%)1.12 (0.60–2.07)
rs879577
Codominant ModelC/C23 (24%)57 (59.4%)1<0.001237.3
C/T48 (50%)31 (32.3%)0.22 (0.10–0.45)
T/T25 (26%)8 (8.3%)0.13 (0.05–0.36)
Dominant ModelC/C23 (24%)57 (59.4%)1<0.001236.2
C/T-T/T73 (76%)39 (40.6%)0.19 (0.10–0.37)
Recessive ModelC/C-C/T71 (74%)88 (91.7%)10.004253.8
T/T25 (26%)8 (8.3%)0.28 (0.11–0.70)
Over dominant C/C-T/T48 (50%)65 (67.7%)10.003253.1
C/T48 (50%)31 (32.3%)0.37 (0.19–0.72)
rs41323645
Codominant ModelG/G42 (43.8%)23 (24%)10.004253
G/A39 (40.6%)43 (44.8%)2.47 (1.19–5.14)
A/A15 (15.6%)30 (31.2%)3.86 (1.62–9.18)
Dominant ModelG/G42 (43.8%)23 (24%)10.002252.1
G/A-A/A54 (56.2%)73 (76%)2.89 (1.47–5.68)
Recessive ModelG/G-G/A81 (84.4%)66 (68.8%)10.025257.1
A/A15 (15.6%)30 (31.2%)2.34 (1.10–4.98)
Over dominant G/G-A/A57 (59.4%)53 (55.2%)10.24260.8
G/A39 (40.6%)43 (44.8%)1.45 (0.78–2.72)
rs4819555
Codominant ModelC/C42 (43.8%)33 (34.4%)10.19260.9
C/T42 (43.8%)45 (46.9%)1.43 (0.72–2.82)
T/T12 (12.5%)18 (18.8%)2.32 (0.91–5.92)
Dominant ModelC/C42 (43.8%)33 (34.4%)10.14260
C/T-T/T54 (56.2%)63 (65.6%)1.61 (0.85–3.07)
Recessive ModelC/C-C/T84 (87.5%)78 (81.2%)10.13259.9
T/T12 (12.5%)18 (18.8%)1.91 (0.81–4.47)
Over dominant C/C-T/T54 (56.2%)51 (53.1%)10.75262.1
C/T42 (43.8%)45 (46.9%)1.11 (0.60–2.05)
Data are reported as frequency/percentage. A Chi-Square test (two-sided) was employed. Logistic binary regression was carried out, adjusted by age, sex, residency (rural vs. urban), family history, BMI percentile, and puberty defined according to “Marshall and Tanner” [53]. The table also provides the odds ratio (OR) with its 95% confidence interval (CI), the p-value, and the Akaike information criterion (AIC) for each model. Bold values indicate p-value < 0.05.
Table 4. The carriage rate of alleles and disease risk.
Table 4. The carriage rate of alleles and disease risk.
AlleleOverallControlsPatientsOR (95%CI)p-Value
rs4819554
A162 (84.4%)73 (76%)89 (92.7%)4.01 (1.63–9.86)0.002
G116 (60.4%)67 (69.8%)49 (51%)0.45 (0.25–0.82)0.012
rs879577
C159 (82.8%)71 (74%)88 (91.7%)3.87 (1.65–9.11)0.002
T112 (58.3%)73 (76%)39 (40.6%)0.22 (0.12–0.4)<0.001
rs41323645
G147 (76.6%)81 (84.4%)66 (68.8%)0.41 (0.2–0.82)0.016
A127 (66.1%)54 (56.3%)73 (76%)2.47 (1.33–4.58)0.006
rs4819555
C162 (84.4%)84 (87.5%)78 (81.3%)0.62 (0.28–1.37)0.32
T117 (60.9%)54 (56.3%)63 (65.6%)1.48 (0.83–2.66)0.24
Data are presented as frequency and percentage. A Chi-Square test (two-sided) was employed. Logistic binary regression was carried out. Unadjusted odds ratio (OR) with its 95% confidence interval (CI) is calculated for each allele in the patient group, using the control group as a reference. Bold values indicate p-values < 0.05.
Table 5. Subgroup analysis by puberty. Haplotype frequency analysis.
Table 5. Subgroup analysis by puberty. Haplotype frequency analysis.
SNP IDGenotypePrepubertyPubertyp-Value *
ControlsPatientsOR (95%CI)ControlsPatientsOR (95%CI)
rs4819554A/A29471141910.33
A/G-G/G67490.49 (0.26–0.93)30260.70 (0.27–1.77)
rs879577C/C11321122510.13
C/T-T/T41190.11 (0.04–0.30)32200.33 (0.13–0.87)
rs41323645G/G22101201310.65
G/A-A/A30413.35 (1.30–8.62)24322.47 (0.95–6.43)
rs4819555C/C26191161410.91
C/T-T/T26321.67 (0.72–3.84)28311.55 (0.58–4.11)
* The p-values presented in the last column test for interaction in the trend, which indicates if the effect of the genotype on disease risk changes with puberty status. Bold values indicate a statistical significance at p-value < 0.05.
Table 6. Subgroup analysis by puberty. Haplotype frequency analysis.
Table 6. Subgroup analysis by puberty. Haplotype frequency analysis.
HaplotypeTotalControlsPatientsCumulative FrequencyOR (95%CI)p-Value
1ACAC0.13910.01360.22870.13911---
2ACGC0.13320.15170.14720.27230.00 (0.00–0.11)0.002
3ATGC0.10160.17770.02180.37390.00 (0.00–0.05)<0.001
4GCAC0.08820.07450.06830.46210.00 (0.00–0.08)0.002
5ACGT0.08790.10640.09230.550.02 (0.00–0.36)0.009
6ATGT0.07010.07840.07250.62010.00 (0.00–0.07)0.001
7GCGC0.05030.04990.05020.67040.06 (0.00–1.46)0.08
8GTGC0.04990.05280.02940.72030.00 (-Inf–Inf)1.00
9ACAT0.04820.0920.01620.76850.02 (0.00–0.48)0.018
10GTAC0.04460.05840.00340.81310.00 (0.00–0.04)<0.001
11ATAT0.04260.02880.07250.8556NA<0.001
12GCAT0.03530.059300.89090.00 (0.00–0.14)0.003
13GCGT0.03290.03810.03360.92370.03 (0.00–1.01)0.05
14GTAT0.0327NA0.0790.95640.00 (0.00–0.12)0.004
15GTGT0.03130.0120.06650.98770.00 (0.00–0.16)0.004
16ATAC0.01230.00660.018210.00 (0.00–1.36)0.06
Bold values indicate a statistical significance of the observed differences at p-value < 0.05.
Table 7. Association of carriage of risky allele on patient characteristics and disease outcomes.
Table 7. Association of carriage of risky allele on patient characteristics and disease outcomes.
Characteristicsrs4819554*A Carrierp-Valuers879577*C Carrierp-Valuers41323645*A Carrierp-Valuers4819555*C Carrierp-Value
NoYesNoYesNoYesNoYes
Total number 789 888 2373 1878
Demographics
AgeMean ± SD9.6 ± 2.99.7 ± 3.10.849.5 ± 3.39.7 ± 30.709.9 ± 3.19.6 ± 3.10.4810.4 ± 3.39.6 ± 30.16
SexMale3 (42.9%)41 (46.1%)0.874 (50%)40 (45.5%)0.8012 (52.2%)32 (43.8%)0.636 (33.3%)38 (48.7%)0.30
ResidencyUrban4 (57.1%)28 (31.5%)0.225 (62.5%)27 (30.7%)0.115 (21.7%)27 (37%)0.218 (44.4%)24 (30.8%)0.28
FH of asthmaPositive5 (71.4%)24 (27%)0.0252 (25%)27 (30.7%)0.747 (30.4%)22 (30.1%)0.984 (22.2%)25 (32.1%)0.57
Pubertal statusPositive5 (71.4%)40 (44.9%)0.254 (50%)41 (46.6%)0.8513 (56.5%)32 (43.8%)0.347 (38.9%)38 (48.7%)0.60
BMI %<85th2 (28.6%)51 (57.3%)0.213 (37.5%)50 (56.8%)0.3812 (52.2%)41 (56.2%)0.7610 (55.6%)43 (55.1%)1.00
<95th3 (42.9%)29 (32.6%) 3 (37.5%)29 (33%) 9 (39.1%)23 (31.5%) 6 (33.3%)26 (33.3%)
≥95th2 (28.6%)9 (10.1%) 2 (25%)9 (10.2%) 2 (8.7%)9 (12.3%) 2 (11.1%)9 (11.5%)
Presentation
OnsetEarly (≤3 y)2 (28.6%)43 (48.3%)0.445 (62.5%)40 (45.5%)0.4710 (43.5%)35 (47.9%)0.8110 (55.6%)35 (44.9%)0.44
Late (>3 y)5 (71.4%)46 (51.7%) 3 (37.5%)48 (54.5%) 13 (56.5%)38 (52.1%) 8 (44.4%)43 (55.1%)
Asthma phenotypeAtopic6 (85.7%)70 (78.7%)0.674 (50%)72 (81.8%)0.1019 (82.6%)57 (78.1%)0.5311 (61.1%)65 (83.3%)0.007
Non-atopic 1 (14.3%)6 (6.7%) 2 (25%)5 (5.7%) 2 (8.7%)5 (6.8%) 3 (16.7%)4 (5.1%)
Exercise-induced0 (0%)11 (12.4%) 2 (25%)9 (10.2%) 1 (4.3%)10 (13.7%) 2 (11.1%)9 (11.5%)
Aspirin-sensitive0 (0%)2 (2.2%) 0 (0%)2 (2.3%) 1 (4.3%)1 (1.4%) 2 (11.1%)0 (0%)
Symptoms Cough7 (100%)87 (97.8%)0.698 (100%)86 (97.7%)0.6722 (95.7%)72 (98.6%)0.4218 (100%)76 (97.4%)0.49
Dyspnea3 (42.9%)53 (59.6%)0.454 (50%)52 (59.1%)0.7217 (73.9%)39 (53.4%)0.0912 (66.7%)44 (56.4%)0.60
Sputum3 (42.9%)51 (57.3%)0.703 (37.5%)51 (58%)0.2916 (69.6%)38 (52.1%)0.169 (50%)45 (57.7%)0.60
Tightness4 (57.1%)60 (67.4%)0.686 (75%)58 (65.9%)0.7116 (69.6%)48 (65.8%)0.8016 (88.9%)48 (61.5%)0.029
Wheezes 4 (57.1%)76 (85.4%)0.096 (75%)74 (84.1%)0.6220 (87%)60 (82.2%)0.7514 (77.8%)66 (84.6%)0.49
Daytime symptoms > 2 wk.2 (28.6%)46 (51.7%)0.445 (62.5%)43 (48.9%)0.7113 (56.5%)35 (47.9%)0.6312 (66.7%)36 (46.2%)0.19
Night awakening0 (0%)16 (18%)0.600 (0%)16 (18.2%)0.344 (17.4%)12 (16.4%)0.912 (11.1%)14 (17.9%)0.73
Reliever use > 2 wks.1 (14.3%)32 (36%)0.421 (12.5%)32 (36.4%)0.2612 (52.2%)21 (28.8%)0.0476 (33.3%)27 (34.6%)0.92
Activity limitations2 (28.6%)27 (30.3%)0.923 (37.5%)26 (29.5%)0.6910 (43.5%)19 (26%)0.137 (38.9%)22 (28.2%)0.40
ComorbiditiesPositive 2 (28.6%)46 (51.7%)0.443 (37.5%)45 (51.1%)0.7215 (65.2%)33 (45.2%)0.153 (16.7%)45 (57.7%)0.003
Allergic rhinitis2 (28.6%)36 (40.4%)0.702 (25%)36 (40.9%)0.4711 (47.8%)27 (37%)0.463 (16.7%)35 (44.9%)0.033
Eczema0 (0%)30 (33.7%)0.092 (25%)28 (31.8%)0.699 (39.1%)21 (28.8%)0.440 (0%)30 (38.5%)0.001
Asthma Severity Mild 5 (71.4%)35 (39.3%)0.205 (62.5%)35 (39.8%)0.4510 (43.5%)30 (41.1%)0.667 (38.9%)33 (42.3%)0.96
Moderate2 (28.6%)38 (42.7%) 2 (25%)38 (43.2%) 8 (34.8%)32 (43.8%) 8 (44.4%)32 (41%)
Severe0 (0%)16 (18%) 1 (12.5%)15 (17%) 5 (21.7%)11 (15.1%) 3 (16.7%)13 (16.7%)
Duration, years Median (IQR)6 (4–9)5 (5–8)0.885.5 (5–7.8)5 (5–8)0.846 (5–9)5 (4–8)0.166 (5–9)5 (4–8)0.08
Asthma ControlWell-controlled5 (71.4%)30 (34.9%)0.133 (37.5%)32 (37.6%)0.385 (22.7%)30 (42.3%)0.134 (23.5%)31 (40.8%)0.24
Partly-controlled2 (28.6%)41 (47.7%) 5 (62.5%)38 (44.7%) 11 (50%)32 (45.1%) 11 (64.7%)32 (42.1%)
Uncontrolled0 (0%)15 (17.4%) 0 (0%)15 (17.6%) 6 (27.3%)9 (12.7%) 2 (11.8%)13 (17.1%)
Investigations
Lab testsHigh total IgE3 (42.9%)33 (37.1%)0.762 (25%)34 (38.6%)0.7111 (47.8%)25 (34.2%)0.329 (50%)27 (34.6%)0.28
Eosinophilia 3 (42.9%)11 (12.4%)0.061 (12.5%)13 (14.8%)0.864 (17.4%)10 (13.7%)0.745 (27.8%)9 (11.5%)0.08
Total IgE (IU/mL)75 (25–245)50 (15–120)0.3343.5 (17.5–105)60 (15–120)0.6980 (15–180)45 (15–120)0.11100 (25–155.8)47.5 (15–120)0.07
Eosinophil (×106/L)240 (20–689)78 (25–235)0.43111 (21.3–243)78 (25–240)0.9678 (25–340)98 (25–240)0.85230 (86–509.8)50 (20–207.5)0.009
Data (except age) are reported as frequency or median (interquartile range). Chi-Square (two-sided), T-test, or Mann–Whitney U-tests were applied. SD: standard deviation, FH: family history, BMI%: body mass index percentile. Bold values indicate a statistical significance of the observed differences at p-value < 0.05.
Table 8. Association of carriage of risky allele on patient characteristics and disease outcomes.
Table 8. Association of carriage of risky allele on patient characteristics and disease outcomes.
Characteristicsrs4819554*A Carrierp-Valuers879577*C Carrierp-Valuers41323645*A Carrierp-Valuers4819555*C Carrierp-Value
NoYesNoYesNoYesNoYes
Total number 30162 33159 65127 30162
Management
AHRNormal3 (42.9)37 (41.6)0.774 (50)36 (40.9)0.747 (30.4)33 (45.2)0.664 (22.2)36 (46.2)0.14
Borderline3 (42.9)25 (28.1) 3 (37.5)25 (28.4) 8 (34.8)20 (27.4) 9 (50)19 (24.4)
Mild/moderate1 (14.3)23 (25.8) 1 (12.5)23 (26.1) 7 (30.4)17 (23.3) 4 (22.2)20 (25.6)
Severe0 (0)4 (4.5) 0 (0)4 (4.5) 1 (4.3)3 (4.1) 1 (5.6)3 (3.8)
Therapy LevelStep 13 (42.9)18 (20.2)0.232 (25)19 (21.6)0.594 (17.4)17 (23.3)0.913 (16.7)18 (23.1)0.92
Step 22 (28.6)17 (19.1) 3 (37.5)16 (18.2) 6 (26.1)13 (17.8) 4 (22.2)15 (19.2)
Step 32 (28.6)14 (15.7) 1 (12.5)15 (17) 4 (17.4)12 (16.4) 4 (22.2)12 (15.4)
Step 40 (0)32 (36) 1 (12.5)31 (35.2) 7 (30.4)25 (34.2) 6 (33.3)26 (33.3)
Step 50 (0)8 (9) 1 (12.5)7 (8) 2 (8.7)6 (8.2) 1 (5.6)7 (9)
Therapeutic responsePFTFVC1 (% predicted)83 (78–89)80 (72–82.5)0.0780 (70.5–88.5)80 (72–83)0.6680 (72–84)80 (72–82.5)0.8281 (77.5–84)78 (72–82)0.16
Pre FEV1 (% predicted)66 (62–80)62 (53–72)0.1170 (48–85.5)62 (55–71.5)0.5662 (50–76)62 (59–74)0.4161.5 (52–72)62 (58–76)0.73
Post FEV1 (% predicted)82 (80–88)80 (70–82)0.04581.5 (71–89.5)80 (70.5–83.5)0.4578 (70–84)80 (71–83)0.9681 (74–86)80 (70–82)0.17
PEFR1 (%predicted)86 (80–93)78 (68–86)0.04685 (66.5–89)79 (68–86)0.3980 (68–86)80 (68–86)0.9085 (71–88)79 (68–84.5)0.22
PC208 (5–8)4 (1–8)0.558 (2.8–8)4 (1.8–8)0.826 (1–8)6 (4–8)0.698 (1–16)4 (3.3–8)0.36
BDRBASE (ml)32.3 (5.1–33.3)32.3 (5.4–39.7)0.3419 (4.7–47.6)32.3 (5.6–38.6)0.8734.5 (7–42.3)31.3 (5.2–37.9)0.2432.6 (22.1–41.1)29.2 (5.2–38.5)0.26
Data are reported as frequency or median (interquartile range). Chi-Square (two-sided) or Mann–Whitney U-tests were applied. “AHR: airway hyperresponsiveness, PFT: pulmonary function test, FVC1: forced vital capacity in 1 s, FEV1: forced expiratory volume in 1 s, PEFR1: peak expiratory flow rate, PC20: provocative concentration causing a 20% drop in FEV1, BDRBASE: bronchodilator response”. Bold values indicate significance at p-value < 0.05.
Table 9. Impact of the risky haplotype block (rs4819554*A/ rs879577*C/ rs41323645*A) on disease outcomes.
Table 9. Impact of the risky haplotype block (rs4819554*A/ rs879577*C/ rs41323645*A) on disease outcomes.
CharacteristicsACA(C/T) Haplotypep-Value
NoYes
Total number 3660
Demographics
AgeMean ± SD9.7 ± 3.19.6 ± 3.00.79
SexMale17 (47.2%)27 (45%)0.84
ResidencyUrban12 (33.3%)20 (33.3%)1.00
FH of asthmaPositive12 (33.3%)17 (28.3%)0.65
Pubertal status *Positive20 (55.6%)25 (41.7%)0.21
BMI %<85th17 (47.2%)36 (60%)0.34
<95th13 (36.1%)19 (31.7%)
≥95th6 (16.7%)5 (8.3%)
Presentation
OnsetEarly (≤3 y)17 (47.2%)28 (46.7%)0.96
Late (>3 y)19 (52.8%)32 (53.3%)
Asthma phenotypeAtopic27 (75%)49 (81.7%)0.24
Non-atopic 5 (13.9%)2 (3.3%)
Exercise-induced3 (8.3%)8 (13.3%)
Aspirin-sensitive1 (2.8%)1 (1.7%)
Symptoms Cough35 (97.2%)59 (98.3%)0.71
Dyspnea22 (61.1%)34 (56.7%)0.83
Sputum20 (55.6%)34 (56.7%)0.92
Tightness24 (66.7%)40 (66.7%)1.00
Wheezes 28 (77.8%)52 (86.7%)0.27
Daytime Symptoms > 2 wks.18 (50%)30 (50%)1.00
Night awakening4 (11.1%)12 (20%)0.40
Reliever use > 2 wks.12 (33.3%)21 (35%)0.87
Activity limitations13 (36.1%)16 (26.7%)0.36
ComorbiditiesPositive 20 (55.6%)28 (46.7%)0.53
Allergic rhinitis15 (41.7%)23 (38.3%)0.83
Eczema11 (30.6%)19 (31.7%)0.91
Asthma severity Mild 18 (50%)22 (36.7%)0.38
Moderate12 (33.3%)28 (46.7%)
Severe6 (16.7%)10 (16.7%)
Asthma controlWell-controlled13 (37.1%)22 (37.9%)0.98
Partly-controlled16 (45.7%)27 (46.6%)
Uncontrolled6 (17.1%)9 (15.5%)
Management
AHRNormal12 (33.3%)28 (46.7%)0.38
Borderline14 (38.9%)14 (23.3%)
Mild/moderate9 (25%)15 (25%)
Severe1 (2.8%)3 (5%)
Therapy levelStep 17 (19.4%)14 (23.3%)0.21
Step 211 (30.6%)8 (13.3%)
Step 37 (19.4%)9 (15%)
Step 48 (22.2%)24 (40%)
Step 53 (8.3%)5 (8.3%)
Lab testsHigh IgE16 (44.4%)20 (33.3%)0.29
Eosinophilia 8 (22.2%)6 (10%)0.14
Therapeutic response
PFTFVC1 (% predicted)80 (72–84)79 (72–82)0.32
Pre FEV1 (% predicted)62 (52–77.5)62 (58–71.5)1.0
Post FEV1 (% predicted)80.5 (72.5–86)79 (70–82)0.23
PEFR1 (% predicted)83 (68.5–87.5)78 (68–84)0.16
PC208.0 (1.0–8.0)4.0 (4.0–8.0)0.91
BDRBASE (mL)32.3 (6.0–41.6)29.1 (5.2–38.3)0.39
Data (except age) are reported as frequency or median (interquartile range). Chi-Square (two-sided), T-test, or Mann–Whitney U-tests were applied. “SD: standard deviation, FH: family history, BMI: body mass index, AHR: airway hyperresponsiveness, IgE: immunoglobulin E, PFT: pulmonary function test, FVC1: forced vital capacity in 1 s, FEV1: forced expiratory volume in 1 s, PEFR1: peak expiratory flow rate, PC20: provocative concentration causing a 20% drop in FEV1, BDRBASE: bronchodilator response, SD: standard deviation”. * According to “Marshall and Tanner” [53].
Table 10. In silico analysis of the impact of the studied IL17RA variants.
Table 10. In silico analysis of the impact of the studied IL17RA variants.
SNP ID and PositionAmino Acid
change
RegionMutation TasterSIFT
Prediction (Score)
PANTHER3D Protein Viewer with a Specified Amino Acid Substitution
rs4819554
c.-947A > G
NAPromoter 808 bp 5′ of the geneNANANAChildren 11 00657 i001
rs879577
c.1100C > A
g.23366C > A
A367VCDS (Exon 13)Likely benignTolerated
0.285
Probably benign
rs41323645
c.2071G > A
g.24337G > A
A691TCDS (Exon 13)BenignTolerated
0.362
Probably benign
rs4819555
c.2160C > G
P720PCDS (Exon 13)Likely benignTolerated
1.0
Possibly damaging
SNP: single nucleotide polymorphism, SIFT: sorting intolerant from tolerant, PANTHER: Protein ANalysis THrough Evolutionary Relationships, A (amino acid): alanin, V: valine, T: threonine, P: prolin, CDS: coding sequences, NA: not applicable. Data sources: “National Center of Biotechnology Institute, SNP database (https://www.ncbi.nlm.nih.gov/snp/), Mutation Taster (https://www.genecascade.org/MutationTaster2021/#chrpos), SIFTonline tool (https://sift.bii.a-star.edu.sg/www/SIFT_dbSNP.html), PANTHER (https://www.pantherdb.org/tools/csnpScoreForm.jsp), VarSome (https://varsome.com/variant/hg38/) according to SwissModel 4nux 342–556. The VarSome 3D protein viewer for rs41323645 and rs4819555-related protein products was unavailable on the applicable model”. All databases were last accessed on 30 July 2023.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Maher, S.A.; AbdAllah, N.B.; Ageeli, E.A.; Riad, E.; Kattan, S.W.; Abdelaal, S.; Abdelfatah, W.; Ibrahim, G.A.; Toraih, E.A.; Awadalla, G.A.; et al. Impact of Interleukin-17 Receptor A Gene Variants on Asthma Susceptibility and Clinical Manifestations in Children and Adolescents. Children 2024, 11, 657. https://doi.org/10.3390/children11060657

AMA Style

Maher SA, AbdAllah NB, Ageeli EA, Riad E, Kattan SW, Abdelaal S, Abdelfatah W, Ibrahim GA, Toraih EA, Awadalla GA, et al. Impact of Interleukin-17 Receptor A Gene Variants on Asthma Susceptibility and Clinical Manifestations in Children and Adolescents. Children. 2024; 11(6):657. https://doi.org/10.3390/children11060657

Chicago/Turabian Style

Maher, Shymaa Ahmed, Nouran B. AbdAllah, Essam Al Ageeli, Eman Riad, Shahad W. Kattan, Sherouk Abdelaal, Wagdy Abdelfatah, Gehan A. Ibrahim, Eman A. Toraih, Ghada A. Awadalla, and et al. 2024. "Impact of Interleukin-17 Receptor A Gene Variants on Asthma Susceptibility and Clinical Manifestations in Children and Adolescents" Children 11, no. 6: 657. https://doi.org/10.3390/children11060657

APA Style

Maher, S. A., AbdAllah, N. B., Ageeli, E. A., Riad, E., Kattan, S. W., Abdelaal, S., Abdelfatah, W., Ibrahim, G. A., Toraih, E. A., Awadalla, G. A., Fawzy, M. S., & Ibrahim, A. (2024). Impact of Interleukin-17 Receptor A Gene Variants on Asthma Susceptibility and Clinical Manifestations in Children and Adolescents. Children, 11(6), 657. https://doi.org/10.3390/children11060657

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