- Open Access
Variants in the IL17 pathway genes are associated with atopic asthma and atopy makers in a South American population
Allergy, Asthma & Clinical Immunologyvolume 15, Article number: 28 (2019)
Asthma is a complex disorder with multiple phenotypes which can influence its severity and response to treatment. The TH17 lymphocytes producing IL-17A and IL17-F cytokines, may have a role on asthma inflammation. The aim of our study was to evaluate the association between genetic variants in IL17 pathway genes with asthma and atopy markers.
Materials and methods
Genotyping was performed using a commercial panel in 1245 participants of SCAALA cohort. The study included 91 SNVs in IL-17 pathway genes. Logistic regressions for asthma and atopy markers were performed using PLINK 1.9. In silico analyses were performed using rSNPbase, RegulomeDB, and Gtex portal for in silico gene expression.
Results and discussion
The T allele of rs1974226 in IL17A was positively associated with asthma (OR: 1.37; 95% CI 1.02–1.82). Also, the T allele of rs279548 was positively associated with asthma (OR: 1.30; 95% CI 1.02–1.64), atopy (OR: 1.62; 95% CI 1.05–2.50) and increased expression of the IL17RC in lung and whole blood tissues. The others genetic variants in the IL17 pathways genes were associated with both protection and risk for asthma development as well as with IgE levels.
The genetic variants in IL-17-related genes are associated with the atopic asthma phenotype and IgE production.
About 334 million people worldwide suffer from asthma, and this number tends to increase . It is expected that in 2025, 100 million new asthma cases will occur worldwide ; thus, it represents a global public health problem, especially in developing countries where the westernized life style led to an increase of such disease . Asthma is one of the most common chronic inflammatory diseases of lower airways, affecting children and young adults . It has a heterogeneous etiopathogeny which arises from various factors associated with a complex genetic basis as well as several environmental factors and individual variability. All these factors together are responsible for the wide variety of inflammatory phenotypes [10, 19, 27], leading to different response to treatment, making it difficult to establish a specific therapy for each phenotypes since most of them are still unknown .
Several genes have been associated with distinct asthma phenotypes. Classically, the allergic asthma is orchestrated by T2-type immune response which is related to the development of allergic inflammation inducing IgE release, mucus secretion and eosinophils chemotaxis to the lungs. More recent studies have showed that, part of asthma patients have basal levels of T2-type cytokines in association to no allergic sensitization, known as non-allergic asthma (non-atopic). The non-atopic asthma may be associated with TH17 T CD4+ lymphocytes . In the presence of allergens and other environmental factors, bronchial epithelial cells are stimulated to synthesize cytokines, among them, IL-6 and IL-23. Such cytokines lead the naïve T cells to activate the transcription factors RORγt and STAT3. Thus, these T cells differentiate into TH17 lymphocytes, which in turn will stimulate the release of IL-17-family cytokines, especially IL-17A and IL-17F . These cytokines stimulate the inflammatory response in bronchial epithelial cells through interaction with IL-17RA receptors that act preferably as dimers with the IL-17RC. The IL17RA/IL17RC activation can initiate a signaling cascade leading to production of downstream inflammatory mediators and cytokines (IL-6 and IL8) induced by NF-kB transcription factor  (Fig. 1). Recent studies show that IL-17A and IL-17F cytokines are overexpressed in the lung tissue and into the bloodstream of some asthmatics patients . Other evidence also have pointed out that, in asthmatic patients, IL-17A expression is increased in lung tissue, sputum and bronchoalveolar lavage fluid (BALF) .
It is well established that genetic polymorphisms play a significant role in asthma and, therefore, interest in this field has grown over the past two decades . In studies involving the Th17 pathway, SNVs in the IL17A and the IL17F were associated with allergic rhinitis and asthma . The most studied variants in the IL17F gene (rs763780 and rs2397084) were associated with asthma . Also, the SNVs rs2275913 and rs2397084 were associated with asthma in children in Tunisia  and with asthma and rhinitis in a Portuguese population . In a Asian population, the rs2275913 was associated with the development of asthma after bronchiolitis in infancy . One study also identified variants in IL17RA associated with the risk of developing aspirin-intolerant asthma . There are no previous studies associating genetic variants in the IL17RC, RORC and NFKB with asthma, however, it is plausible to believe that variants that may interfere in the regulation of these genes play an important role in asthma.
Collectively, these findings indicate that variations on IL-17 pathway genes may be involved in the pathophysiology of atopic asthma and suggest that the production of IL-17 and its accessory downstream molecules are candidates for asthma susceptibility. Therefore, the aim to this study was to evaluate how variants in IL-17-related genes can influence the development of atopy and asthma in a Latin urban population from Brazil.
The studied population comprised 1245 children born between 1994 and 2001, originally recruited for Social Changes Asthma and Allergy in Latin America (SCAALA) program in the city of Salvador, BA, Brazil, as previously described . Briefly, demographic and environmental data were collected from standardized questionnaires. Asthma diagnosis was done by the application of a Portuguese translated, ISAAC phase II questionnaire. These questionnaires were applied to each child parents or legal guardian in 2005. Blood samples for laboratory tests and isolation of genetic material were collected and skin prick test were done with several allergens in a field ambulatory. This work has been approved by the National Research Ethics Committee (reference number: 120.616) and free informed consent was properly obtained from the parents or legal guardian of each child.
Skin prick test and allergens specific IgE (sIgE)
All individuals were submitted to skin puncture tests to common inhalant allergens in our region (mites—D. pteronyssinus, D. farinae, B. tropicalis, cockroach—P. americana and B. germanica; fungi mix (ALK-ABELLO, São Paulo, Brazil). After 15 min, the wheal diameter was measured. The test was considered positive if the mean of the largest perpendicular diameters (excluding pseudopodia) was at least 3 mm greater than the negative control (Table 1).
Determination of specific IgE serum concentrations was performed for Dermatophagoides pteronyssinus, Blomia tropicalis, Blattella germanica, and Periplaneta Americana using the ImmunoCAP assay (Phadia Diagnostics AB, Uppsala Sweden). Children with ≥ 0.70 kU/L or greater of specific IgE for any allergen tested were considered to have a positive result .
Asthma symptoms and atopy definitions
Atopy was defined if the child had specific IgE ≥ 0.70 kU/L and/or Skin prick test ≥ 3 mm for at least one tested aeroallergens .
Asthma (current wheezing) was defined by having wheezing in the last 12 months, plus at least 1 of the following: (1) diagnosis of asthma ever; (2) wheezing with exercise in the last 12 months; (3) 4 or more episodes of wheezing in the last 12 months; and (4) waking up at night because of wheezing in the last 12 months. All children with asthma and atopy were classified as atopic asthmatics. The others with asthma and without atopy were considered as non-atopic asthmatics.
Genomic DNA extraction
DNA was extracted from peripheral blood samples according to the protocol Flexigene DNA Kit (Qiagen, Hilden, Germany). All samples were standardized at a gDNA concentration of 5 ng/μL and stored at − 30 °C until use.
The genotyping was performed using a commercial panel of Illumina BeadChip Human Omni2.5-8 Kit (www.illumina.com) with 2284818 SNVs, through the Consortium EPIGEN-Brazil. The information on the genes position was obtained from National Center for Biotechnology Information (NCBI) (www.ncbi.nlm.nih.gov) and used to localize genetic data corresponding to each gene as following: IL17A from NC_000006.11 (52051185.0.52055436) position at the chromosome 6, IL17F from NC_000006.11 (52101484.0.52109298) position at the chromosome 6, IL17RA from NC_000022.10 (17565849.0.17596584) position at the chromosome 12, IL17RC from NC_000003.11 (9958758.0.9975305) chromosome 3, NFKB from NC_000004.11 (103422486.0.103538459) chromosome 4, and the RORC from NC_000001.10 (151778547.0.151804348) position at the chromosome 1. A total of 79 SNVs were analyzed in this study related to genes on IL-17 pathway. Also, 277 informative markers of ancestry were genotyped and used to correct the population stratification. The selection of these markers was based on the allelic frequency difference above 50% between populations present in the HGDP-CEPH projects and 1000 Genomes.
Were included in this study, markers with genotype call rates > 98%, Hardy–Weinberg equilibrium (HWE) with p > 0.05 and population stratification and frequency of minor allele (MAF) with p > 0.05.
Logistic regressions were applied to estimate odds ratio (OR), permutational-p value, and 95% confidential interval (CI) for the association between SNVs, asthma and atopy markers (skin tests and specific IgE production), adjusted for sex, age, helminth infections and ancestry markers. The analyses were performed using the additive genetic model. A permutational p-value lower than 0.05 was considered statistically significant. The permutation test is performed to test the null hypothesis, a difference in the values is expected under the null hypothesis. Permutation procedures provide a computationally intensive approach to generating significance levels empirically. This test control the false discovery rate in order to solve the problem of multiple comparisons , preserving the correlational structure between SNVs . All the analyses above were performed using the software PLINK (version 1.9). The Linkage Disequilibrium among cases and controls was created using Haploview 4.2.
Genetic risk score for asthma symptoms
The genetic score was performed to determine the risk degree in the presence of more than one allele of certain SNVs associated with asthma symptoms in this study. Therefore, it was included in the analyses only the genotype combination of these variants. The score number was attributed according to the presence of minor allele.
All the genetic scores analyses were performed using the SNPstats web tool (www.snpstats.net/snpstats/). The forest graph was created using the GraphPad Prism 5.0 software.
In silico functional analyses
To identify the SNVs functions the National Center of Biotechnology Information (NCBI) (www.ncbi.nlm.nih.gov) was accessed, which is an open online platform that contains worldwide data about biotechnology, biomedicine and bioinformatics. It houses databases with specifics material, and among them there is the SNV database (dbSNP). The dbSNP is a helpful tool that provides extensive information about SNVs such as the SNV position, the chromosome, allele frequency for each allele variation; also, whether the SNV is intronic, exonic or near to a gene, whether it is a missense variation, which amino acids is changed. dbSNP was used to gather all the information cited above.
The statistical analyses were made using the method of Benjamini and Hochberg , and corrected by Bonferroni test. Here, the Homo sapiens organism and the option to evaluated a lists of target genes were chosen.
RegulomeDB (regulomedb.org) is a database for interpretation of regulatory variants in the human genome. It includes high-throughput, experimental data sets from Encyclopedia of DNA Elements (ENCODE) and other sources. RegulomeDB identifies putative regulatory potential and functional variants through computational predictions and manual annotations . The greater evidence for a variant to be located in a functional region the lower the score in RegulomeDB is, as follows: 1a to 1f scores, SNVs are likely to affect the binding and the expression of a gene target; 2a to 2c scores, SNVs are likely to affect binding; scores 3a and 3b are less likely to affect binding; 4, 5 and 6 scores, SNVs may have minimal binding evidence and; 7 score, there are no data about the function of a certain SNV .
rSNPBase (http://rsnp.psych.ac.cn/) is a database that provides annotations focused on regulatory SNVs involved in a wide range of regulation types, including proximal, distal and post-transcriptional regulation, to identify their potentially regulated genes , in respect of: (i) proximal regulation indicates SNVs that are involved in proximal transcriptional regulation. The (ii) distal regulation, SNVs should be involved in distal transcriptional regulation. The (iii) micro RNA regulation describes SNVs within mature miRNA and the (iv) RNA binding protein are SNVs involved in RNA binding protein-mediated post-transcriptional regulation. Through the rSNPbase platform the information about each genetic variants was accessed in UniProt database, and Ensembl genome browser (wgEncodeUwDgf), a directory containing the downloadable files associated with Encyclopedia of DNA Elements (ENCODE) data sources (Uniport, 2014).
The genotype-tissue expression project GTEx (www.gtexportal.org) of the National Institutes of Health Common Fund establishes a resource database in association with tissue bank made with the objective to study the relationship between gene expression and genetic variation, and other molecular phenotypes, in multiple reference human tissues . Expression quantitative trait loci (eQTL) mapping offers a powerful approach to elucidate the genetic component underlying altered gene expression. Genetic variation can also influence gene expression through alterations in splicing, non-coding RNAs, and RNA stability. Gene expression is differentially regulated across tissues, and many human transcripts are expressed in a limited set of cell types or during a limited developmental stage . From this project information was obtained about gene expression according to each SNV.
Characteristics of the study population
77.30% of the participants were non-asthmatic individuals, while 22.70% were asthmatics. No significant difference was found for age range and sex between the two groups. On the other hand, there was a significant statistical difference in specific IgE production for at least one allergen, specific IgE for D. pteronyssinus, B. tropicalis, B. germanica and P. americana when comparing the non-asthmatics group with asthmatic group. Also, there was a statistical difference in positive skin prick test response for at least one allergen, D. pteronyssinus, B. tropicalis and B. germanica among the two groups (Table 1).
We found 150 SNVs for the six genes in the IL-17 pathway. Of these, 26 SNV was excluded by Tests of Hardy–Weinberg equilibrium (HWE) (p ≤ 0.001) and 53 SNVs by minor allele frequency (MAF < 0.05). No sample was excluded by low genotyping (mind > 0.1) or missingness (geno > 0.1). After quality control steps, the study included 79SNVs.
Variants in IL17 pathway genes can influence asthma
In the IL17A, the T allele of rs1974226 was positively associated with asthma (OR: 1.37; 95% CI 1.02–1.82) (Table 2). In the IL17RC, the T allele of rs279548 (OR: 1.30; 95% CI 1.02–1.64), the T allele of rs11917994 (OR: 1.40; 95% CI 1.09–1.79) and the A allele of rs76234423 (OR: 1.39; 95% CI 1.05–1.83), was positively associated with asthma (Table 2). The allele C of rs6769465 was negatively associated with asthma (OR: 0.60; 95% CI 0.34–0.97) (Table 2). In the NFKB1, the A allele of rs35680095 (OR: 1.28; 95% CI 1.02–1.61) and the G allele of rs75071695 (OR: 3.25; 95% CI 1.19–8.87) were positively associated with asthma (Table 2). In the RORC gene, the allele C of rs4995918 was positively associated with asthma (OR: 1.28; 95% CI 1.04–1.57) (Table 2).
Variants in IL17 can influence atopic makers
In IL17RA, the allele C of rs10483089 was negatively associated with SPT for at least one allergen (OR: 0.76; 95% CI 0.60–0.97 (Table 3). In IL17RC, the A allele of rs115461448 was positively associated with anti-D. pteronyssinus IgE (OR: 3.09; 95% CI 1.20–8.00) and SPT for one least one allergen (OR: 2.58 and 95% CI 1.01–6.62) (Tables 3 and 4). The allele A of rs77569961 was negatively associated with production of anti-D. pteronyssinus IgE (OR: 0.69 and 95% CI 0.51–0.95), sIgE for B. tropicalis (OR: 0.74; 95% CI 0.57–0.97) and positive SPT for D. pteronyssinus (OR:0.68; 95% CI 0.47–0.98) (Tables 3 and 4). The allele A of rs279545 was positively associated with sIgE production for at least one allergen (OR: 1.35; 95% CI 1.03–1.76), sIgE for D. pteronyssinus (OR: 1.26; 95% CI 1.04–1.53) and sIgE for B. tropicalis (OR: 1.23 95% CI 1.04–1.50) (Table 4). These variants are not in linkage disequilibrium (Fig. 3d).
In NFKB1, the allele G of rs73837255 was negatively associated with anti-D. pteronyssinus and anti-B. tropicalis IgE production, SPT for D. pteronyssinus (OR: 0.69; 95% CI 0.51–0.95), for B. tropicalis (OR: 0.61; 95% CI 0.39–0.96) and for at least one tested allergen (OR: 0.63; 95% CI 0.42–0.93) (Tables 3 and 4). The A allele of rs73837241 was negatively associated with sIgE production for D. pteronyssinus (OR: 0.57; 95% CI 0.40–0.89) and for B. tropicalis (OR: 0.70; 95% CI 0.49–0.99) and SPT for D. pteronyssinus (OR: 0.62; 95% CI 0.40–0.95) (Tables 3 and 4). These SNVs were not in linkage disequilibrium (Fig. 2f). The allele T of rs4647992 was positively associated with anti-D. pteronyssinus (OR: 1.64; 95% CI 1.18–2.30), B. germanica (OR: 1.64; 95% CI 1.10–2.42) and B. tropicalis IgE production (OR: 1.61; 95% CI 1.15–2.26) (Table 4).
In RORC, the allele A of rs11578418 was positively associated with sIgE production against B. tropicalis (OR: 1.30; 95% CI 1.02–1.65) (Table 4), and positively associated with SPT for D. pteronyssinus (OR: 1.41; 95% CI 1.05–1.90) (Table 3). The allele T of rs11204894 was positively associated with sIgE anti-B. tropicalis (OR: 1.27; 95% CI 1.04–1.56) (Table 4) and positively associated with SPT for D. pteronyssinus (OR: 1.30; 95% CI 1.01–1.68) and for B. tropicalis (OR: 1.33 and 95% CI 1.06–1.67) (Table 3). Moreover, the rs11578418 and rs11204894 were in total linkage disequilibrium (Fig. 2e). The allele T of rs7540799 was positively associated with SPT for B. tropicalis (Table 3). These SNVs (rs11204894 and rs7540799) were in high linkage disequilibrium (Fig. 3e). The allele T of rs3790515 was negatively associated with SPT for D. pteronyssinus (OR: 0.62 and 95% CI 0.46–0.85) and for B. tropicalis (OR: 0.74; 95% CI 0.58–0.96) (Table 3), this variant was in total linkage disequilibrium with rs78703675 (Fig. 3).
Variants can increases the risk of asthma symptoms
Initially all genotypes combination of the rs1974226, rs11917994, rs76234423, rs279548, rs35680095, rs75071695 and rs4995918 were analyzed (data not shown). It was found that the risk of asthma susceptibility was increased only when the alleles that acted as risk factors (rs76234423, rs279548, and rs3568009) were presented together.
Genetic score analyses showed that among the individuals, the most frequent genotype was that one which has all wild-type alleles (26.52%) in the rs35680095 (G/G), rs76234423 (G/G) and rs279548 (C/C); these individuals were the reference group (score 0). Individuals were grouped according to the amount of minor allele they presented. The total of subjects included in the score “1” was 529, representing 45.5% of total population. The score “2” were present in 23.92% of population. The score “3” appeared in 5.82% of the subjects. The score “4” composed a minor part of subjects, only 0.72%. And, the least frequent score was the score “5”, with only two individuals (0.16%) (Table 5).
Statistical analyses between the gene score and asthma symptoms showed that, in the presence of one polymorphic allele (score 1), the risk of asthma symptoms was increased (OR: 1.86; 95% CI 1.09–3.15) when adding another risk allele the OR did not change too much (OR: 95% CI 1.44–4.52). However, in the presence of three minor allele (score 3), the risk of symptoms of asthma is approximately three times larger (OR: 3.25; 95% CI 1.38–7.65). All analyses were compared with the reference group (Fig. 4). Because of the low numbers of individuals with the score 4 and score 5 in our population, the gene score analysis could not be performed.
In silico IL17RC expression is increased in the presence of variants
The eQTL was examined in silico to obtain expression of the each gene according to the genotype in lung and whole blood when available in Gtex portal. All SNVs described in this study were evaluated. Nevertheless, no difference in gene expression was found for IL17A, IL17F, IL17RA, NFKB1 and RORC according to genotype for the SNVs explored in this study (data not showed). On the other hand, gene expression analysis for IL17RC (Fig. 2), the T allele of rs279548 increased the IL17RC expression in lung and whole blood (p-value ≤ 0.05); the A allele of rs279545 also increases the IL17RC expression in the whole blood (Fig. 2a–c), (p-value ≤ 0.05).
Regulatory features of SNVs
Table 6 summarizes the results of rSNPBase in which regulatory features of retrieved SNVs in four (i–iv) regulation manners are summarized as ‘Yes’ or ‘No’. All SNVs in the IL17A, IL17F, IL17RA and IL17RC interfere in proximal transcription regulation, which is related to regulatory elements associated with DNA accessibility. In NFKB and RORC the variants involved in this process were few. Most of SNVs involved in the distal regulation were present in IL17RC, NFKB and RORC; this SNVs can modify the chromatin interactions. All SNVs in NFKB are involved with miRNA regulation, only SNVs in this gene were involved with this type of regulatory feature. The miRNA are important to control the levels of RNA expression. Also, most of the SNVs operate in the post-transcriptional RNA binding protein mediated regulation.
For the significant variants found in this study, RegulomeDB was used to identify putative regulatory potential and functional SNVs (Table 6). This database has a score ranging from 1 to 7, as previously described in Methods section (Table 3). Few variants were found with high involvement in expression mechanisms (score 1–3); most of the variants were involved in transcription biding factors.
All the SNVs studied herein had some level of evidence in terms of having an impacting gene expression.
In our study, it was shown for the first time that variants in IL-17 pathway can interfere on asthma and atopy in a Brazilian population.
Studies have demonstrated the important role of IL-17 pathway in pulmonary diseases , specially in asthma, considering that the overexpression was demonstrated of IL-17A and IL17F in lungs of non-atopic asthmatic patients . Also, Molet et al.  demonstrated that in asthmatic subject, IL-17A is increased in lungs, sputum and bronchoalveolar fluid. In our work it was found that the variants in genes from IL-17 pathway can increase the risk for developing asthma and atopy; this fact is probably due to an increase of IL-17 levels gene expression in lung of asthmatics subjects. All variants that increase the risk of developing asthma are involved in proximal regulation, and this type of regulation is largely associated with gene expression .
The actions of IL-17A and IL-17F depend on the appropriate link to their receptors IL17RA and IL17RC . It was showed that variants in IL17RC are associated with increased risk of asthma (rs11917994, rs76234423 and rs279548) and atopy (rs279545, rs7627060, rs115461448, rs77569961) development. Among variants in IL17RC, the T allele of rs279548 increases the expression of IL17RC in lung and whole blood and the A allele of rs279545 increases the expression of IL17RC in whole blood. Function analyses suggest that they can be involved in proximal, distal, transcriptional regulation and post-transcriptional regulation. Moreover, this variant opens the chromatin in chr6:52055240-52055390 position, on TH17 cells lineage. Also, the post transcriptional regulation of the SNV is involved in the regulation of ELAVL1 (ELAV like RNA binding protein 1); the protein encoded by this gene is closely related with the mRNAs stabilization and gene expression .
The pro-inflammatory signalization of IL17RA and IL17RC are largely linked to NF-kB . NF-κB is a transcription factor expressed in numerous cell types, which plays a key role in the pro-inflammatory cytokines production . Variants were found in NFKB1 (rs35680095 and rs75071695) elevating the risk of asthma and atopy development; furthermore, according to our results, all variants in this gene are involved in inactivation of miRNA. Interestingly, genetic score analyses demonstrated that combination variants in IL17RC and NFKB1 together increased the risk of asthma in comparison with the presence of these variants alone. It was found that variants in NFKB1 have an important role in asthma according to studies that indicate enhanced NF-κB pathway activation in asthmatic tissues , and the NF-kB expression is increased in the airway epithelium of asthmatic humans .
Another transcription factor related with asthma is the RORγt; it has been described as a major factor that determines the phenotype of airway inflammation and steroid sensitivity in asthma . In our study a variant (rs4995918) was also found in this gene increasing the asthma risk.
In other way, some variants were also found as a protect factor against asthma (rs6769465) and atopy (rs77569961, rs73837255, rs73837241 and rs73837255). These SNVs are located in non-coding region and variants in these regions can alter the RNA splicing process and interfere with exons junctions impacting directly on the protein translation . Moreover, it was found that some of these protecting SNVs block the distal transcriptional, which can lead to a low expression of pro-inflammatory genes involved in asthma and atopy.
In conclusion, variant in genes of IL-17 pathways may influence in the development course of asthma pathology and atopy. Also, more studies are necessary to further elucidate the potential role of IL17 pathways genes on asthma and atopy, and how it could be a strategy to control this disease.
expression quantitative trait loci
genomic deoxyribonucleic acid
minor allele frequency
nuclear kappa B factor
related orphan receptor gamma
single nucleotide variant
Al-Ramli W, Prefontaine D, et al. T(H)17-associated cytokines (IL-17A and IL-17F) in severe asthma. J Allergy Clin Immunol. 2009;123(5):1185–7.
Asher I, Pearce N. Global burden of asthma among children. Int J Tuberc Lung Dis. 2014;18(11):1269–78.
Barreto ML, Cunha SS, et al. Risk factors and immunological pathways for asthma and other allergic diseases in children: background and methodology of a longitudinal study in a large urban center in Northeastern Brazil (Salvador-SCAALA study). BMC Pulm Med. 2006;6:15.
Barreto ML, et al. Prevalence of asthma symptoms among adolescents in Brazil: national adolescent school-based health survey (PeNSE 2012). Rev Bras Epidemiol. 2014;17(Suppl 1):106–15.
Bazzi MD, Sultan MA, et al. Interleukin 17A and F and asthma in Saudi Arabia: gene polymorphisms and protein levels. J Investig Allergol Clin Immunol. 2011;21(7):551–5.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. Roy. Statist. Soc. B. 1995;57:289–300.
Boyle AP, et al. Annotation of functional variation in personal genomes using RegulomeDB. Genome Res. 2012;22(9):1790–7.
Chunfang G, Wu L, and Li X. IL-17 family: cytokines, receptors and signaling. Cytokine 2013;64(2):477–485
Cooper PJ, Rodrigues LC, et al. Influence of poverty and infection on asthma in Latin America. Curr Opin Allergy Clin Immunol. 2012;12(2):171–8.
Dixon AL, Liang L, et al. A genome-wide association study of global gene expression. Nat Genet. 2007;39(10):1202–7.
Dzhugashvili M, Luengo-Gil G, et al. Role of genetic polymorphisms in NFKB-mediated inflammatory pathways in response to primary chemoradiation therapy for rectal cancer. Int J Radiat Oncol Biol Phys. 2014;90(3):595–602.
Figueiredo CA, et al. Environmental conditions, immunologic phenotypes, atopy, and asthma: new evidence of how the hygiene hypothesis operates in Latin America. J Allergy Clin Immunol. 2013;131(4):1064.e1–1068.e1.
Gaffen S. IL-17 receptor composition. Nat Rev Immunol. 2016;16:4.
Gagliardo R, et al. Persistent activation of nuclear factor-kappaB signaling pathway in severe uncontrolled asthma. Am J Respir Crit Care Med. 2003;168(10):1190–8.
Guo L, et al. rSNPBase: a database for curated regulatory SNPs. Nucleic Acids Res. 2014;42(Database issue):D1033–9.
Holster A, Teräsjärvi J, et al. IL-17A gene polymorphism rs2275913 is associated with the development of asthma after bronchiolitis in infancy. Allergol Int. 2018;67(1):109–13.
Ivanov II, McKenzie BS, et al. The orphan nuclear receptor RORgammat directs the differentiation program of proinflammatory IL-17+ T helper cells. Cell. 2006;126(6):1121–33.
Iniguez LP, Hernandez G. The evolutionary relationship between alternative splicing and gene duplication. Front Genet. 2017;8:14.
Kabesch M, Depner M, et al. Polymorphisms in eosinophil pathway genes, asthma and atopy. Allergy. 2007;62(4):423–8.
Kawaguchi M, Takahashi D, et al. IL-17F sequence variant (His161Arg) is associated with protection against asthma and antagonizes wild-type IL-17F activity. J Allergy Clin Immunol. 2006;117(4):795–801.
Ki-Hyun K, Jahan SA, and Kabir E. A review on human health perspective of air pollution with respect to allergies and asthma. Environment international 59 (2013):41–52
Korn T, Bettelli E, et al. IL-17 and Th17 cells. Annu Rev Immunol. 2009;27:485–517.
Lage-Castellanos A, Martinez-Montes E, et al. False discovery rate and permutation test: an evaluation in ERP data analysis. Stat Med. 2010;29(1):63–74.
Ma L, Xue HB, et al. The Imbalance of Th17 cells and CD4(+) CD25(high) Foxp3(+) Treg cells in patients with atopic dermatitis. J Eur Acad Dermatol Venereol. 2014;28(8):1079–86.
Maalmi H, Beraies A, et al. IL-17A and IL-17F genes variants and susceptibility to childhood asthma in Tunisia. J Asthma. 2014;51(4):348–54.
Melling N, Taskin B, Hube-Magg C, Kluth M, Minner S, Koop C, Grob T, Graefen M, Heinzer H, Tsourlakis MC, Izbicki J, Wittmer C, Huland H, Simon R, Wilczak W, Sauter G, Steurer S, Schlomm T, Krech T. Cytoplasmic accumulation of ELAVL1 is an independent predictor of biochemical recurrence associated with genomic instability in prostate cancer. Prostate. 2016;76(3):259–72. https://doi.org/10.1002/pros.23120.
Moffatt MF. Genes in asthma: new genes and new ways. Curr Opin Allergy Clin Immunol. 2008;8(5):411–7.
Molet S, Hamid Q, et al. IL-17 is increased in asthmatic airways and induces human bronchial fibroblasts to produce cytokines. J Allergy Clin Immunol. 2001;108(3):430–8.
Molet SM, Hamid QA, et al. IL-11 and IL-17 expression in nasal polyps: relationship to collagen deposition and suppression by intranasal fluticasone propionate. Laryngoscope. 2003;113(10):1803–12.
Park SM, Park JS, et al. Unraveling the genetic basis of aspirin hypersensitivity in asthma beyond arachidonate pathways. Allergy Asthma Immunol Res. 2013;5(5):258–76.
Pascual RM, Peters SP. Asthma. Med Clin North Am. 2011;95(6):1115–24.
Peters SP. Asthma phenotypes: nonallergic (intrinsic) asthma. J Allergy Clin Immunol Pract. 2014;2(6):650–2.
Purcell S, Neale B, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559–75.
Resende E, Todo-Bom A, et al. Asthma and rhinitis have different genetic profiles for IL13, IL17A and GSTP1 polymorphisms. Revista Portuguesa de Pneumologia. 2017;23(1):10–6.
Lonsdale J, Thomas J, Salvatore M, Phillips R, Lo E, Shad S, Hasz R, Walters G, Garcia F, Young N, Foster B. The genotype-tissue expression (GTEx) project. Nat Genet. 2013;45(6):580–5.
Wang M, Zhang Y, et al. Association between polymorphisms in cytokine genes IL-17A and IL-17F and development of allergic rhinitis and comorbid asthma in Chinese subjects. Hum Immunol. 2012;73(6):647–53.
Zhu M, Wang T, Chen R, Wang C, Liu S, Ji Y. Association between interleukin-17a gene polymorphisms and asthma risk: a meta-analysis. Asian Pac J Allergy Immunol. 2016;34(2):115–23.
MJS did the statistical analyzes, in silico analyzes, genetic analyzes and drafted the figures and manuscript; MBRS did the statistical analyzes; BRT did the statistical analyzes; RPE drafted the manuscript; NMAN carried out the molecular genetic studies, participated in the sequence alignment; MLB carried out the molecular genetic studies, participated in the sequence alignment; CAF carried out the molecular genetic studies, participated in the sequence alignment and drafted the manuscript; RSC: drafted the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
The authors consents publication.
Ethics approval and consent to participate
This work has been approved by the National Research Ethics Committee (reference number: 120.616) and free informed consent was properly obtained from the parents or legal guardian of each child.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.