- Letter to the Editor
- Open Access
An association analysis to identify genetic variants linked to asthma and rhino-conjunctivitis in a cohort of Sicilian children
Italian Journal of Pediatrics volume 45, Article number: 16 (2019)
Asthma and rhino-conjunctivitis are common chronic diseases in childhood. In this cross-sectional study, we performed a gene association analysis with current asthma and rhino-conjunctivitis in a cohort of Sicilian children aged 10–15 years. Overall, our findings reveal the importance of different genetic variants at 4p14, 16p12.1, 17q12, 6p12.2 and 17q21.1, identifying possible candidate genes responsible for susceptibility to asthma and rhino-conjunctivitis.
To the Editor,
Asthma and rhino-conjunctivitis (RC) are common diseases worldwide that are frequently associated. Observed differences in prevalence of asthma and RC may be explained by genetic susceptibility, though environmental factors play a relevant role . In order to increase genomic information on Sicilian children, this research has explored some genetic variants to discover possible association with asthma and RC.
A representative sample of 1050 children within the “Palermo Junior High School” (PJHS II)  study were investigated through questionnaires, spirometry, and skin prick test (SPT) to quantify the prevalence of asthma and RC, in association with allergic sensitization and respiratory function, and to evaluate the role of environmental and host risk factors for allergic respiratory diseases. The study was approved by the local Institutional Ethical Committee. All parents of the enrolled children signed a written informed consent.
Two different phenotypes were identified: Current Asthma (CA) defined as asthma ever plus at least a wheeze episode in the last 12 months, RC defined as sneezing, or runny, or blocked nose apart from common cold or flu in the last 12 months and nose problem accompanied by itching and/or watering eyes. The concomitant presence of CA and RC was merged into the CA group; children without CA and RC (nAnRC) were used as controls.
A total of 52 Single Nucleotide Polymorphisms (SNPs), involved in the innate immune system pathways were selected for genotyping by Matrix-Assisted Laser Desorption/Ionization (MALDI-TOF-MS). Out of the 52 initially selected SNPs, 7 were complete drop-outs and the remaining SNPs were successfully genotyped. The individuals were genotyped with the Illumina Bead-Chip (Illumina Inc., San Diego, CA, USA); the PLINK v1.07 software was used to perform standard quality control. SNPs were excluded if they had low call rates (proportion of genotyped called < 90%), were not in Hardy-Weinberg equilibrium (HWE, p < 0.001) on the nAnRC subjects, or had a low minor allele frequency (MAF < 1%). A total of 22 SNPs were used for further analyses.
Mean values were compared among children with CA, RC and nAnRC using the analysis of variance (ANOVA). Differences of categorical variables were evaluated using Chi-squared test. Associations between single SNPs and CA and RC were analysed by applying the case/control model of the SNPassoc R package, adjusting for sex, age, body mass index (BMI), SPT+ (at least one positive), exposure to current environmental tobacco smoke and traffic.
The demographic and lung function characteristics of the 1050 subjects are shown in Table 1. The study sample was composed by 523 (49.8%) Female and 527 (50.2%) Male, aged 12.07 ± 0.74 years on average. Subjects were categorized into CA (n = 61), RC (n = 184) and nAnRC (n = 805). Subjects with CA and RC more frequently had SPT+; subjects with CA were younger than nAnRC and RC subjects.
Chromosome, gene, SNP name, quality control tag, alleles coding (Major/minor), minor allele frequency (MAF), test for Hardy-Weinberg Equilibrium, percentage of missing values (%) and genotyping distribution are reported in Table 2.
By applying the case/control model, no SNP reached the Bonferroni corrected significance threshold (P value < 0.002, i.e., 0.05/#tests), and only one SNP reached the Bonferroni corrected suggested significance threshold (P value < 0.005, i.e., 0.10/#tests). However, we also included those SNPs reaching the nominal significance threshold (P value < 0.05) just to highlight modest associations with the two studied phenotypes. CA was strongly associated only with rs4252665 and modestly with rs1801275 and rs17616434; RC was modestly associated with rs7741835, rs8079416, rs3859192, rs3894194 and rs7212938. Table 3 shows the genotypic frequencies of the associated SNPs in CA/RC and nAnRC groups, and the adjusted OR and 95% CI from the logistic regression model for the only SNP strongly associated. The SNP rs4252665 showed significantly different genotypic frequencies between the two groups, i.e., the CA subjects had a high frequency of CT heterozygote genotype compared with nAnRC, which were mostly homozygous. Indeed, using the overdominant genetic model, in which the baseline is the homozygous genotype (CC/TT), the CT genotype of rs4252665 showed a large increased risk of CA (ORC/T = 5.75; 95% CI = 2.03–16.29). The SNPs modestly associated with CA showed a high frequency of the major allele homozygote genotype compared with nAnRC, in which the genotypes were mostly characterized by the presence of the minor allele. Furthermore, with respect to nAnRC, SNPs modestly associated with RC showed small variations in the genotypic frequencies. In particular, the SNPs rs8079416, rs7741835 and rs3894194 showed high frequencies of heterozygote genotypes compared with nAnRC, which are frequently homozygous, whilst the SNPs rs7212938 and rs3859192 had high frequencies of the minor allele homozygote genotypes compared with nAnRC, in which genotypes are characterized by the presence of the major allele.
The detected association signals for CA were located within the Toll-like receptor (TLR1) on chromosome 4, the interleukin 4 receptor (IL4R) on chromosome 16 and the Erb-b2 receptor tyrosine kinase 2 (ERBB2) on chromosome 17. It is known that Toll-like receptors (TLRs) represent a major group of receptors for the specific recognition of pathogen-associated molecular patterns of microbes capable of activating innate and adaptive immunity that reduce the risk for asthma . The IL4R gene is known to encode a protein that regulates IgE production and it has been shown that allelic variations in this gene are associated with atopy, allergic rhinitis and asthma . Recently, some loci of ERBB2, which belong to the encoding region 17q12, have been reported to be in linkage disequilibrium with loci in the region 17q21 encoding gasdermin A (GSDMA) gene, previously associated with childhood asthma [5,6,7].
With regard to the RC, the modestly associated genes were interleukin 17 (IL17) on chromosome 6, leucine rich repeat containing 3C (LRRC3C), and GSDMA on chromosome 17. IL17 is a pro-inflammatory cytokine that targets epithelial cells48 and its expression in the nasal mucosa has been associated with allergic rhinitis and its degree of severity [8, 9]. To our knowledge, no functional studies have been published on LRRC3C, although, within the human genome, the gene LRRC32 has been associated with eczema and allergic rhinitis , and probably some similarities between the two proteins encoded by LRRC3C and LRRC32 exist. Finally, GSDMA gene has been associated with childhood asthma and allergic disease in many populations. In particular, region 17q21 has been originally identified in the first GWAS on childhood asthma , and GSDMA variants were suggested to be strong risk factors for asthma and airway inflammation .
Overall, our findings reveal the importance of different genetic variants at 4p14, 16p12.1 17q12, 6p12.2 and 17q21.1, identifying possible candidate genes responsible for CA and RC in the Sicilian child population. These results are a preliminary step in understanding the pathophysiology of asthma and rhino-conjunctivitis in a paediatric population in the Mediterranean area and need to be verified by further studies using more advanced technologies. Furthermore, novel methodologies combining genome-wide association study (GWAS; ) and expression quantitative trait locus (eQTL ) such as summary-data based Mendelian randomization (SMR; ), PrediXcan , MetaXcan , would be useful in discovering new genetic variants linked to these allergic respiratory diseases in this geographic area. Unlike traditional single-variant tests, these innovative approaches based on SNP-gene linkage will provide valuable insights on disease causality. Noteworthy, the integrative analysis of GWAS and eQTL studies, by identifying gene-trait-associated changes in the expression, would mitigate some tasks associated with a GWAS approach, allowing the discover of genetic variants which can affect gene expression . Moreover, since some gene functions are often pleiotropic, this combined approach would allow a better comprehension of the pathways through which pleiotropy can affect clinical phenotypes.
In conclusion, the present study could facilitate the application of novel therapeutics and preventive strategies arising from the genomics era of precision medicine.
Body mass index
Expression quantitative trait locus
Human epidermal growth factor receptor 2
Genome-wide association study
Interleukin 4 receptor
Leucine rich repeat containing 3C
Minor allele frequency
Not asthma and not rhino-conjunctivitis
Summary-data based Mendelian randomization
Single Nucleotide Polymorphism
Skin prick tests
Galassi C, De Sario M, Biggeri A, et al. Changes in prevalence of asthma and allergies among children and adolescents in Italy: 1994–2002. Pediatrics. 2006;117(1):34–42.
Cibella F, Cuttitta G, La Grutta S, et al. Proportional Venn diagram and determinants of allergic respiratory diseases in Italian adolescents. Pediatr Allergy Immunol. 2011;22:60–8.
Kormann MS, Depner M, Hartl D, et al. Toll-like receptor heterodimer variants protect from childhood asthma. J Allergy Clin Immunol. 2008;122(1):86–92.
Huang ZY, Cheng BJ, Cai GJ, et al. Association of the IL-4R Q576R polymorphism and asthma in the Chinese Han population: a meta-analysis. Genet Mol Res. 2015;14(1):2900–11.
Stein MM, Thompson EE, Schoettler N, et al. A decade of research on the 17q12–21 asthma locus: Piecing together the puzzle. J Allergy Clin Immunol. 2018;142(3):749–64 e3.
Verlaan DJ, Berlivet S, Hunninghake GM, et al. Allele-specific chromatin remodeling in the ZPBP2/GSDMB/ORMDL3 locus associated with the risk of asthma and autoimmune disease. Am J Hum Genet. 2009;85(3):377–93.
Ferreira MA, Matheson MC, Tang CS, et al. Genome-wide association analysis identifies 11 risk variants associated with the asthma with hay fever phenotype. J Allergy Clin Immunol. 2014;133(6):1564–71.
Albano GD, Di Sano C, Bonanno A, et al. Th17 immunity in children with allergic asthma and rhinitis: a pharmacological approach. PLoS One. 2013;8(4):e58892.
Ramasamy A, Curjuric I, Coin LJ, et al. A genome-wide meta-analysis of genetic variants associated with allergic rhinitis and grass sensitization and their interaction with birth order. J Allergy Clin Immunol. 2011;128(5):996–1005.
Resende EP, Todo-Bom A, Loureiro C, et al. Asthma and rhinitis have different genetic profiles for IL13, IL17A and GSTP1 polymorphisms. Rev Port Pneumol. 2017;23(1):10–6.
Ozaki K, Ohnishi Y, Iida A, et al. Functional SNPs in the lymphotoxin-alpha gene that are associated with susceptibility to myocardial infarction. Nat Genet. 2002;32(4):650–4.
Brem RB, Yvert G, Clinton R, et al. Genetic dissection of transcriptional regulation in budding yeast. Science. 2002;296(5568):752–5.
Zhu Z, Zhang F, Hu H, et al. Integration of summary data from GWAS and eQTL studies predicts complex trait gene targets. Nat Genet. 2016;48(5):481–7.
Gamazon ER, Wheeler HE, Shah KP, et al. A gene-based association method for mapping traits using reference transcriptome data. Nat Genet. 2015;47(9):1091–8.
Barbeira AN, Dickinson SP, Bonazzola R, et al. Exploring the phenotypic consequences of tissue specific gene expression variation inferred from GWAS summary statistics. Nat Commun. 2018;9(1):1825.
Hauberg ME, Zhang W, Giambartolomei C, et al. Large-scale identification of common trait and disease variants affecting gene expression. Am J Hum Genet. 2017;100:885–94.
We are grateful to Dr. Maximillian Schieck and to Prof. M. Kabesch for supporting gene analysis. We also thank all the school staff, children, and parents who made the study possible.
The present research was fully supported by a grant from the Regional Agency for Environment Protection (ARPA Sicilia) [DDG No. 303/2005].
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
The study was approved by the local ethics committee (A.O.U.P. “Paolo Giaccone”, Palermo, Italy), and written informed consent was provided by parents of all participants.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.