Association of Polymorphic Variants of HHIP, ADRB2 and IL-33 Genes with Clinical Manifestations of Bronchial Asthma in Children =============================================================================================================================== * Yu.S. Alieva * E.G. Furman * E.I. Kondratyeva * E.V. Loshkova * V.S. Sheludko * V.S. Sokolovsky * M.S. Ponomareva * E.A. Khuzina * R.S. Aushova ## Abstract Studying the contribution of genetic mechanisms to the development of bronchial asthma (BA) is to look for associations of the disease and its phenotypes with polymorphic markers of candidate genes. **Objective** To investigate the association of polymorphic variants of HHIP, ADRB2 and IL-33 genes with the phenotypes of clinical course of BA in children and the effectiveness of therapy of the disease **Patients and methods** A cohort single-center study of 90 bronchial asthma patients aged 5 to 17 years with an established diagnosis of bronchial asthma of varying degree of severity and control was conducted. The allele and genotype frequencies of polymorphic loci of the following genes were studied: rs12551256-A and rs146597587-G of IL-33 gene; rs12504628 of HHIP gene and ARG16GLY rs1042713 of ADRB2 gene in 90 BA patients with regard to the severity and control of the disease. In children with severe BA, as well as in children with poorly controlled/uncontrolled asthma (n=26), sequencing of the entire coding sequence of the IL-33 gene located on the 9th chromosome in the 9p24.1 region was additionally performed (search for mutations in 9 exons). **Results** Comparison of genetic markers in patients with severe BA (tBA) and non-severe BA (nBA) revealed a reduced risk of severe disease realization among those carrying the TT genotype (OR=0.221 (95% CI: 0.059-0.828; χ2=5.759; p=0.056)) and the T allele (OR=0.491 (95% CI: 0.190-1.269; χ2=4.270; p=0.039)) of the studied genetic variant rs12504628 (T>C) of the HHIP gene, the frequency of the CC genotype in severe BA was 64%, versus 28% in nonsevere BA, and the C allele 77% versus 52%.Comparison of genetic markers in patients with a combination of atopic dermatitis (AtD) and bronchial asthma (BA+AtD) and BA without AtD (BA without AtD) revealed an increased risk of combining asthma and dermatitis among individuals carrying the TT genotype (OR=2.875 (95% CI: 1.130-7.316; χ2=5.751; p=0.056)) of genetic variant rs12504628 (T>C) of the HHIP gene. Sequencing and exome analysis of the IL-33 gene showed a statistically significant positive association between the frequency of lesions in exons 4 (r=0.417; p=0.034) and 6 (r=0.593; p=0.001) on the one hand and the severity of BA on the other. Nucleotide substitutions in these exons were found to be more frequently associated with the severe course of bronchial asthma. **Conclusion** It was shown that TT genotype of genetic variant rs12504628 (T>C) of HHIP gene reduces the risk of severe BA, but increases the risk of atopic dermatitis combined with BA by 2.8 times. The CC+ST genotype of the HHIP gene increases the risk of drug allergy against the background of BA by 2.9 times. Polymorphic variants in exons 4 and 6 of the IL-33 gene are more often combined with moderate and severe asthma, and nucleotide substitutions in exons 4 and 6 are associated with a severe course of BA. ## Introduction Bronchial asthma (BA) is one of the most important problems of modern theoretical and practical medicine. According to the report of The Global Asthma Network, currently, about 348 million people suffer from this disease, at least 14% of them are children [1]. It is estimated that in 2019, the number of patients with BA was 262 million and 461,000 deaths from the disease were reported [2]. BA is a heterogeneous disease based on chronic airway inflammation characterized by reversible bronchoobstructive syndrome, recurrent episodes of wheezing, dyspnoea, chest congestion and cough. The symptoms vary in time and intensity [3, 4]. The main goals of treatment of patients with BA are to achieve and maintain optimal control of the disease and to prevent exacerbations [5]. Despite a wide BA availability of inhaled glucocorticosteroids (iGCS) and standardized guidelines for asthma management, disease control remains suboptimal in most children. More than 50% of all children with BA experience at least one exacerbation each year, including children with non-severe asthma. In Russia, the problem of asthma control is also very urgent, as only 23% of patients achieve complete control of the disease. Among the reasons for insufficient control of BA are low adherence to therapy (43%), lack of elimination of triggers (29%), presence of comorbidities (15%), smoking (15%) and others [6]. Single nucleotide substitutions in the genome account for the impact of genetic polymorphism on the disease phenotype, and predetermine differences in the clinical manifestations of the disease, including symptom control. A characteristic feature of molecular medicine as a science based on data on the molecular structure of the human genome is its focus on the correction of the pathological process in a particular person, taking into account his or her unique genetic features. Another peculiarity is preventive orientation, when information about genome obtained long before obvious manifestations of a disease can prevent its development. Genetic predisposition can manifest itself in interaction with environmental factors, that forms a pathological phenotype. The most common method of studying the contribution of genetic mechanisms to the development of BA is to look for associations of the disease and its phenotypes with polymorphic markers of candidate genes. Over the last decade, genetic studies have identified numerous candidate genes, which predispose to BA. However, the results reported in the literature are often contradictory, that confirms the need to further investigate the associations of BA with polymorphic markers of candidate genes. Some polymorphic markers of candidate genes may contribute to the development of airway obstruction due to loss of lung elasticity, while others contribute to chronic inflammation leading to airway obstruction or poor response to medications such as β2-agonists or (iGCS) [7]. Genetic aspects of BA control in children continue to be studied. In particular, it is known that polymorphism of the β2-adrenergic receptor gene (ADRB2 gene) is associated with the therapeutic response of patients to bronchodilators - β2-agonists. Stimulation of β2-adrenergic receptors leads to bronchodilation and improvement of bronchial conduction, and also affects T-cell function, eosinophilic inflammation and causes a decrease in the secretion of proinflammatory mediators from mast cells [8]. The product of the HHIP gene (Hedgehog family interacting protein gene) is an evolutionarily conserved signaling protein, playing an important role in a wide range of processes. There is evidence of an association of the single-nucleotide polymorphism rs1828591 of the HHIP gene with predisposition to the development of bronchial obstruction [9]. The association of the single nucleotide polymorphism rs1512288 of the HHIP gene with reversibility of bronchial obstruction has been reported, while no such association has been found with bronchial hyperresponsiveness [10]. Among a significant number of BA candidate genes, cytokine-alarmin genes, which play a key role in all stages of allergic reactions, occupy a special place. In the direction of molecular genetic methods, the interleukin-33 (IL-33) cytokine-alarmin gene is promising for research. IL-33 is one of the central signaling molecules of immune reactions in BA [11]. The role of IL-33, has been confirmed in the pathogenesis of BA in children [12,13]. IL-33 levels have been found to be elevated in the sputum and bronchial biopsy specimens of asthmatic patients [14]. IL-33 is a tissue-derived cytokine, which induces and enhances eosinophilic inflammation and has emerged as a promising new target for the treatment of asthma and allergic diseases [15]. Binding of IL-33 to its receptor (ST2) increases the expression of several pro-inflammatory mediators (IL-5, IL-4 and IL-13), and thereby affects Th2-cell-mediated eosinophilic airway inflammation. Increased IL-33 was associated with the presence of the rs1342326 allele of the IL-33 gene. The IL-33 gene rs1342326 polymorphism was associated with a lower risk of asthma in children in Tunisian population and higher expression of IL-33 cytokine [16]. Single nucleotide polymorphism rs992969 of IL-33 gene, was reported to be associated with blood eosinophil levels, asthma and eosinophilic asthma. The polymorphism rs4008366 of IL-33 gene, showed a weak association with eosinophilic asthma [17]. A whole-genome sequencing in the Icelandic population revealed a rare variant of the IL-33 gene (NM_001199640:exon7:c.487-1G>C (rs146597587-C), allele frequency = 0.65%), which disrupts the canonical acceptor splice site upstream of the last coding exon of the gene [18]. This variant also occurs with low frequency in European population and is associated with lower eosinophil counts and reduced risk of asthma in Europeans (OR = 0.47; 95%). Heterozygotes have approximately 40% lower total IL-33 mRNA expression than non-carriers. This polymorphism results in the formation of a shortened form of IL-33 protein. The shortened variant does not form an IL-33R/ST2 complex and does not activate ST2-expressing cells. These data demonstrate that rs146597587-C is a loss-of-function cytokine mutation [18]. Despite a number of reported studies of HHIP, ADRB2 and IL-33 gene polymorphisms in BA, the significance of the association of polymorphisms of the above genes with the clinical course of BA in children has not yet been definitively determined. In connection with the above-mentioned, there is a need to improve the comprehensive assessment of the degree of BA control and to determine the influence of clinical, laboratory, functional and genetic characteristics on it. ### Aim of the study To investigate the association of polymorphic variants of HHIP, ADRB2 and IL-33 genes with the phenotypes of clinical course of BA in children and the effectiveness of therapy of the disease ## Materials and Methods A cohort single-center study of 90 bronchial asthma patients aged 5 to 17 years with an established diagnosis of bronchial asthma of varying degree of severity and control between November 2019 and March 2021 was conducted. The diagnosis of BA was established on the basis of current clinical guidelines. The study was conducted on the basis of the Regional Children’s Clinical Hospital of Perm and polyclinics of the Perm region within the framework of the RFBR research grant. Inclusion criteria: children with diagnosed BA aged from 5 to 17 years; Exclusion criteria: any acute respiratory infections at the time of examination, age less than 5 years (due to the impossibility of spirography in this age group). All participants underwent a comprehensive examination of clinical condition and external respiratory function. Subsequently, a set of diagnostic procedures was implemented for them, including the study of genetic polymorphism of HHIP, ADRB2 and IL-33 genes to establish the relationship with clinical phenotypes, laboratory and instrumental examination parameters determining the course of bronchial asthma and the degree of disease control. All patients were collected anamnesis of the disease, including allergological anamnesis, general blood analysis, rhinocytogram, level of total and specific IgE, immunogram as indicated. Spirographic examination, pulse oximetry, chest X-ray examination, peak expiratory flow rate indices were performed. Based on the clinical parameters and respiratory function, the severity degree of exacerbation by clinical guidelines was identified. The severity of BA exacerbations was determined according to the following clinical criteria - clinical symptoms, PSV, respiratory rate, pulse rate, frequency of emergency medication use, and nocturnal awakenings. The level of BA control was determined according to Asthma Control Test - ACT and C-ACT and the Composite Asthma Severity Index (CASI). The material for the molecular genetic study was DNA isolated from dried capillary blood spots from 90 children. The allele and genotype frequencies of polymorphic loci of the following genes were studied: rs12551256-A and rs146597587-G of IL-33 gene in 70 children; rs12504628 of HHIP gene and ARG16GLY rs1042713 of ADRB2 gene in 90 BA patients with regard to the severity and control of the disease. The polymerase chain reaction (PCR) method was used to detect mutant alleles of genes. In children with severe BA, as well as in children with poorly controlled/uncontrolled asthma (n=26), sequencing of the entire coding sequence of the IL-33 gene located on the 9th chromosome in the 9p24.1 region was additionally performed (search for mutations in 9 exons). The selection of cases was carried out by the method of continuous sampling. The methodology for determining the sample size was based on the use of a specialized formula, when the size of the general population is unknown: ![Formula][1] where n is the sample size, t is the coefficient depending on the confidence level chosen by the researcher, p is the proportion of respondents with the presence of the trait under study, q = 1-p - the proportion of respondents who do not have the characteristic under study, Δ - marginal sampling error. Statistical processing of the results was carried out using Microsoft Excel 2010 statistical software packages. The hypothesis of normality of distribution of the studied indicators was tested using the Shapiro-Wilk criterion. In order to extend the conclusions to the general population (95% reliability), some indicators were presented as M ±2m (% ±2m). The two-sample Student’s t-test or its nonparametric analogue, the Mann-Whitney U-test was used to compare the dependent and independent groups of features characterizing the level of control and/or severity of the disease (depending on the type of distributions of the analyzed parameters). Pearson’s χ2 test was used to analyze the conjugation tables. The relationship between the variables was studied using correlation analysis. The significance level for the tested hypotheses was assumed to be 0.05. The expected distribution of genotypes was estimated using the Hardy-Weinberg formula: ![Formula][2] where q is the frequency of occurrence of recessive gene, p is the frequency of occurrence of dominant gene, q² is the frequency of occurrence of genotype aa, p² is the frequency of occurrence of genotype AA; 2pq is the frequency of occurrence of genotype Aa. The χ2 test was used to assess whether the observed genotype distribution matched the expected one based on Hardy-Weinberg equilibrium. For statistically significant parameters, relative risk (RR) was calculated using the formula: ![Formula][3] The confidence (95%) interval for the OR was calculated using the formulas: ![Formula][4] A, B, C, D in the formulas are the number of observations in the cells of the four-field conjugation table. ### Ethical review This study was conducted in accordance with the principles of the World Medical Association Declaration of Helsinki. Informed consent was obtained in writing from the child’s legal representative in accordance with local laws and regulations before the child was included in the study. Conclusion of the local ethical committee at PSMU No. 5/20 dated 4 August 2020. ## Results and Discussion All children of the study group were on pulmonologist’s dispensary registration with the diagnosis of bronchial asthma. The median age of the patients was 13 years [Q1-Q3: 9; 15 years]. A total of 100 children participated in the study, of whom 72 (72%) were boys and 28 (28%) were girls. In 62 patients BA had a mild course, 27 children were observed with moderate BA and 11 children with severe asthma. The largest number of patients in the study cohort (67 patients) had incomplete control of bronchial asthma, 7 children were diagnosed with no control and 26 patients had complete control. When collecting anamnestic data, 87% of children were found to have hypoallergenic factors (presence of pets, carpets, flowering indoor plants, mould, passive smoking, gassing of the place of residence with exhaust fumes from cars or nearby industrial plants). Outpatient medical records of children under dispensary supervision by a district paediatrician were analyzed to assess the comorbid background. Having studied the comorbidities of the examined children, allergic rhinitis - 80.0%, atopic dermatitis - 32.0% and pollinosis - 40.4% accounted for the highest percentage (Table 1). Most of the children had two or more comorbidities. View this table: [Table 1:](http://medrxiv.org/content/early/2023/12/15/2023.12.14.23299853/T1) Table 1: Comorbidities in children with BA In the structure of complaints of the examined children, there dominated dyspnoea during physical activity (70.0 ± 9.0), dyspnoea outside in spring and summer (34.0 ± 9.3), respiratory infection (42.0 ± 9.7), labored breathing with wheezing (45.0 ± 9.8) and dry attack-like cough (48.0 ± 9.8). The examination of external respiratory function using spirography revealed obstructive type disorders in 29 children, and 18% of children showed significant deterioration of parameters and their shift to the “red” zone, that once again confirms the hypothesis about the lack of disease control in BA patients on the background of baseline therapy. ### Characteristics of the studied genotypes The distribution of genotypes in BA patients corresponded to the Hardy-Weinberg equilibrium, except for the genetic variant rs146597587 of IL-33 gene (G>C), where only carriers of one GG genotype were present (Table 2). View this table: [Table 2:](http://medrxiv.org/content/early/2023/12/15/2023.12.14.23299853/T2) Table 2: Distribution of genotype frequencies of the studied polymorphisms (rs12504628 (T>C), rs1042713 (G>A), rs12551256, rs146597587 (G>C)) in the group of bronchial asthma patients ### Search for associations of genetic markers of HHIP, ADRB2 and IL-33 genes with the clinical course of bronchial asthma Comparison of genetic markers in patients with severe BA (tBA) and non-severe BA (nBA) revealed a reduced risk of severe disease realization among those carrying the TT genotype (OR=0.221 (95% CI: 0.059-0.828; χ2=5.759; p=0.056)) and the T allele (OR=0.491 (95% CI: 0.190-1.269; χ2=4.270; p=0.039)) of the studied genetic variant rs12504628 (T>C) of the HHIP gene, the frequency of the CC genotype in severe BA was 64%, versus 28% in nonsevere BA, and the C allele 77% versus 52% (Table. 3). View this table: [Table 3:](http://medrxiv.org/content/early/2023/12/15/2023.12.14.23299853/T3) Table 3: Case-control analysis of the studied genetic variants in severe and nonsevere BA Comparison of genetic markers in patients with a combination of atopic dermatitis (AtD) and bronchial asthma (BA+AtD) and BA without AtD (BA without AtD) revealed an increased risk of combining asthma and dermatitis among individuals carrying the TT genotype (OR=2.875 (95% CI: 1.130-7.316; χ2=5.751; p=0.056)) of genetic variant rs12504628 (T>C) of the HHIP gene. Analysis of genetic markers in patients with a combination of congenital malformations (CM) of the tracheobronchial tree with asthma (BA+CM) and bronchial asthma without malformations (BA without CM) revealed associations with the combination of asthma and CM with the AA genotype (OR=0.182 (95% CI: 0.051-0.646; χ2=8.567; p=0.014)) genetic variant rs1042713 (G>A) of the ADRB2 gene, which has a protective value against CM, it was shown that AA genotype carriers had bronchial tree CM in 40% of cases, compared to 11% of AA genotype carriers without CM. Analysis of genetic markers in patients with BA and aggravated allergy history among 1st degree relatives (BA+ hereditary history) and bronchial asthma without aggravated hereditary allergy history (BA without hereditary history) revealed that AA genotype carriers (OR=0.112 (95% CI: 0.013-0.932; χ2=5.554; p=0.062)) and carriers of the A allele (OR=0.453 (95% CI: 0.213-0.964; χ2=3.537; p=0.059)) genetic variant rs1042713 (G>A) of the ADRB2 gene, have a lower incidence of an aggravated allergy history among 1st degree relatives (4% vs 26% for genotype AA and 26: vs 44% for allele A). Analysis of genetic markers in patients with BA and drug allergy (BA+DA) and bronchial asthma without drug allergy (BA without DA) revealed that carriers of the CC+TS genotype were more likely to have a combination of asthma and drug allergy (OR=2.917 (95% CI: 1.009-8.427; χ2=4.984; p=0.083)) genetic variant rs12504628 (T>C) of the HHIP gene, a protective role was shown for the T allele in relation to the realization of drug allergy (31% vs 52%) (OR=0.416 (95% CI: 0.190-0.909; χ2=4.204; p=0.040)). Analysis of genetic markers in patients with uncontrolled BA and bronchial asthma with partial and complete control of disease symptoms showed no associations with the genetic variants studied. Carrying the TT genotype of genetic variant rs12504628 (T>C) of the HHIP gene reduces the risk of severe BA, but increases the risk of concomitant atopic dermatitis by 2.8 times. Carriage of CC+CT genotypes increases the risk of drug allergy on the background of BA by 2.9 times (Table 4). View this table: [Table 4:](http://medrxiv.org/content/early/2023/12/15/2023.12.14.23299853/T4) Table 4: Associations of bronchial asthma with genetic variants studied Carrying the AA genotype of the ADRB2 gene is associated with a decreased risk of allergic history among 1st degree relatives and a decreased risk of congenital malformations of the tracheobronchial tree against the background of BA (Table 4). Sequencing and exome analysis of the IL-33 gene showed a statistically significant positive association between the frequency of lesions in exons 4 (r=0.417; p=0.034) and 6 (r=0.593; p=0.001) on the one hand and the severity of BA on the other. Nucleotide substitutions in these exons were found to be more frequently associated with the severe course of bronchial asthma. The study [14] shows that the homozygous variant on Arg16 of the ADRB2 gene is associated with an increased risk of developing a more severe form of bronchial asthma, as well as its nocturnal form compared to homozygotes on the Gly16 allele. An analysis of 28 published studies on the association of beta2-adrenoreceptor gene polymorphisms with BA phenotypes confirmed an association between the Gly16 polymorphism and nocturnal asthma, but no association was found between the Arg16Gly variant and bronchial hyperresponsiveness [20, 21]. In our study, we were able to confirm that the AA genotype of the ADRB2 gene was associated with a decreased risk of aggravated allergoanamnesis among 1st degree relatives and a decreased risk of realization of congenital malformations of the tracheobronchial tree against the background of BA. It should be noted that children with the indicated AA genotype of the ADRB2 gene in our sample did not receive long-acting b2-agonists as monotherapies of baseline therapy and did not use monotherapy with short-acting b2-agonists as emergency drugs during the last 6 months. Consequently, it is not possible to assess the risk of disease exacerbation in Arg16 genotype carriers who resort to the use of b2-agonists. The HHIP gene has now been found to affect both small and large airways [22]. The presence of allele A of the rs13118928 polymorphism of the HHIP gene may be associated with the emphysema-hyperinflation phenotype in patients with chronic obstructive pulmonary disease [23]. It is known that the state of external respiratory function is the most important criterion of bronchial asthma severity. We found that carrying the TT genotype of the genetic variant rs12504628 (T>C) of the HHIP gene reduces the risk of severe BA, but increases the risk of atopic dermatitis combined with BA by 2.8 times. The presence of CC+CT genotypes increases the risk of drug allergy against the background of BA by 2.9 times. On the one hand, we did not detect any association of genetic variants of the IL-33 gene rs12551256 polymorphism with the clinical course of BA. When analyzing the genotypes of the rs146597587 (G>C) polymorphism of the IL-33 gene, all children were carriers of the same genotype. On the other hand, nucleotide substitutions in exons 4 and 6 of the IL-33 gene were found to be associated with a severe course of bronchial asthma. This justifies the expediency of further study of IL-33 gene exon polymorphism and its associations with the clinical course of bronchial asthma in children with an increased patients’ sample size. ## Study limitations The association of HHIP, ADRB2 and IL-33 gene polymorphisms in children cannot be extrapolated to the entire population of Russian children due to the small size of the study sample. It is possible that when the sample size is increased, the distributions of genotypes and alleles of these genes will differ from those presented in this article. We did not perform multivariate analysis with correction for the detected gene associations, taking into account the carriage of polymorphic variants of other genes and environmental factors, that may affect the results of the evaluation of the effect of the studied genes. ## Conclusion The study of HHIP, ADRB2 and IL-33 gene polymorphisms in BA children with different disease phenotypes revealed an association between gene polymorphisms and disease severity, as well as with comorbidities. It was shown that TT genotype of genetic variant rs12504628 (T>C) of HHIP gene reduces the risk of severe BA, but increases the risk of atopic dermatitis combined with BA by 2.8 times. The CC+ST genotype of the HHIP gene increases the risk of drug allergy against the background of BA by 2.9 times. The AA genotype of the ADRB2 gene is associated with the absence of an aggravated allergoanamnesis among 1st degree relatives and a decreased risk of congenital malformations of the tracheobronchial tree against the background of BA. Polymorphic variants in exons 4 and 6 of the IL-33 gene are more often combined with moderate and severe asthma, and nucleotide substitutions in exons 4 and 6 are associated with a severe course of BA. This suggests that exons 4 and 6 should be paid attention to in clinical practice to predict the course of the disease and timely correction of baseline therapy. Thus, this study has established associations of polymorphic variants of HHIP, ADRB2 and IL-33 genes with clinical manifestations of bronchial asthma in children, which can be taken into account in the personalized monitoring of these patients and can help to achieve complete control of the disease. This research was supported by a joint grant from the Israel Ministry of Science and Technology (MOST, 3-16500), the Russian Centre for Scientific Information (RFBR) (Joint Research Project 19-515-06001) ## Data Availability All data produced in the present study are available upon reasonable request to the authors * Received December 14, 2023. * Revision received December 14, 2023. * Accepted December 15, 2023. * © 2023, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.Global Initiative for Asthma. Global Strategy for Asthma Management and Prevention, 2020. Available from: [http://www.ginasthma.org](http://www.ginasthma.org). 2. 2.Global burden of 369 diseases and injuries in 204 countries and territories, 1990-2019: a systematic analysis for the Global Burden of Disease Study 2019. Lancet. 2020; 396 (10258): 1204–1222. doi:10.1016/S0140-6736(20)30925-9 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30925-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F15%2F2023.12.14.23299853.atom) 3. 3.Zhmurov D.V., Parfenteva M.A., Semyonova Y.V. Bronchial asthma. Colloquium-journal. 2020; 14: 66–72. doi: 10.24411/2520-6990-2020-11894 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.24411/2520-6990-2020-11894&link_type=DOI) 4. 4.Ovsyannikov, D. Yu. Bronchial asthma in children / D. Yu. Ovsyannikov, E. G. Furman T. I. Eliseeva. - Moscow : Peoples’ Friendship University of Russia (PFUR), 2019. -211 c. - ISBN 978-5-209-09256-8. 5. 5.Bateman ED, Hurd SS, Barnes PJ, et al. Global strategy for asthma management and prevention: GINA executive summary. Eur. Respir. J. 2008 Jan; 31 (1): 143–178. doi: 10.1183l of Medicine, 2000. 342(19): p. 1378–1384. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiZXJqIjtzOjU6InJlc2lkIjtzOjg6IjMxLzEvMTQzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMTIvMTUvMjAyMy4xMi4xNC4yMzI5OTg1My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 6. 6.Dudenbostel, T. and D. C Allegra L, Cremonesi G, Girbino G, et al. Real-life prospective study on asthma control in Italy: cross-sectional phase results. Respir. Med. 2012; 106 (2): 205–214. doi:10.1016/j.rmed.2011.10.001 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.rmed.2011.10.001&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22035853&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F15%2F2023.12.14.23299853.atom) 7. 7.Cherkashina I., Nikulina S., Logvinenko N., Maksimov V., Liberdovskaya E. Clinical and genetic analysis of patients with bronchial asthma. Pulmonology. 2009; 2: 77–81. doi:10.18093/0869-0189-2009-2-77-81. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18093/0869-0189-2009-2-77-81&link_type=DOI) 8. 8.Ponomaryova, M. S., Furman, E. G., Khuzina, E. A., Yarulina, A. M., Zhdanovich, E. A. (2015). Familial polymorphism of ADRB2 gene in bronchial asthma in childhood. Perm Medical Journal, 32(5), 30–36. 9. 9.Cherkashina, I. I., Razvodovskaya, A. V., Nikulina, S. Yu. Polymorphisms of some genes in bronchial asthma patients of Krasnoyarsk residents // Pulmonology. - 2016. - T. 26. - No. 3. - c. 293–302. 10. 10.Li X, Howard TD, Moore WC, Ampleford EJ, Li H, Busse WW, Calhoun WJ, Castro M, Chung KF, Erzurum SC, Fitzpatrick AM, Gaston B, Israel E, Jarjour NN, Teague WG, Wenzel SE, Peters SP, Hawkins GA, Bleecker ER, Meyers DA. Importance of hedgehog interacting protein and other lung function genes in asthma. J Allergy Clin Immunol. 2011 Jun;127(6):1457–65. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jaci.2011.01.056&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21397937&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F15%2F2023.12.14.23299853.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000291048500022&link_type=ISI) 11. 11.Korppi M, Teräsjärvi J, Lauhkonen E, et al. IL33 rs1342326 gene variation is associated with allergic rhinitis at school age after infant bronchiolitis. Acta Paediatr. 2020 Oct; 109 (10): 2112–2116. doi: 10.1111/apa.15175. PMID: 31955459. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/apa.15175&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31955459&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F15%2F2023.12.14.23299853.atom) 12. 12.Wang Y, Wang L, Hua S. Interleukin-33 in children with asthma: A systematic review and meta-analysis. Allergol. Immunopathol. (Madr.). 2017 Jul-Aug; 45 (4): 387–392. doi: 10.1016/j.aller.2016.12.007. PMID: 28410870. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.aller.2016.12.007&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28410870&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F15%2F2023.12.14.23299853.atom) 13. 13.Saikumar Jayalatha AK, Hesse L, Ketelaar ME, et al. The central role of IL-33/IL-1RL1 pathway in asthma: From pathogenesis to intervention. Pharmacol. Ther. 2021 Sep; 225: 107847. doi: 10.1016/j.pharmthera.2021.107847. PMID: 33819560. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pharmthera.2021.107847&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33819560&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F15%2F2023.12.14.23299853.atom) 14. 14.Hamzaoui A, Berraies A, Kaabachi W, et al. Induced sputum levels of IL-33 and soluble ST2 in young asthmatic children. J. Asthma. 2013 Oct; 50 (8): 803–809. doi: 10.3109/02770903.2013.816317. PMID: 23855553. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3109/02770903.2013.816317&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23855553&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F15%2F2023.12.14.23299853.atom) 15. 15.Cayrol C, Girard JP. Interleukin-33 (IL-33): A nuclear cytokine from the IL-1 family. Immunol. Rev. 2018 Jan; 281 (1): 154–168. doi: 10.1111/imr.12619. PMID: 29247993. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/imr.12619&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29247993&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F15%2F2023.12.14.23299853.atom) 16. 16.Charrad R, Kaabachi W, Berraies A, Hamzaoui K, Hamzaoui A. IL-33 gene variants and protein expression in pediatric Tunisian asthmatic patients. Cytokine. 2018 Apr;104:85–91. 17. 17.Ketelaar ME, Portelli MA, Dijk FN, Shrine N, Faiz A, Vermeulen CJ, Xu CJ, Hankinson J, Bhaker S, Henry AP, Billington CK, Shaw DE, Johnson SR, Benest AV, Pang V, Bates DO, Pogson ZEK, Fogarty A, McKeever TM, Singapuri A, Heaney LG, Mansur AH, Chaudhuri R, Thomson NC, Holloway JW, Lockett GA, Howarth PH, Niven R, Simpson A, Tobin MD, Hall IP, Wain LV, Blakey JD, Brightling CE, Obeidat M, Sin DD, Nickle DC, Bossé Y, Vonk JM, van den Berge M, Koppelman GH, Sayers I, Nawijn MC. Phenotypic and functional translation of IL33 genetics in asthma. J Allergy Clin Immunol. 2021 Jan;147(1):144–157 18. 18.Smith D, Helgason H, Sulem P, Bjornsdottir US, Lim AC, Sveinbjornsson G, Hasegawa H, Brown M, Ketchem RR, Gavala M, Garrett L, Jonasdottir A, Jonasdottir A, Sigurdsson A, Magnusson OT, Eyjolfsson GI, Olafsson I, Onundarson PT, Sigurdardottir O, Gislason D, Gislason T, Ludviksson BR, Ludviksdottir D, Boezen HM, Heinzmann A, Krueger M, Porsbjerg C, Ahluwalia TS, Waage J, Backer V, Deichmann KA, Koppelman GH, Bønnelykke K, Bisgaard H, Masson G, Thorsteinsdottir U, Gudbjartsson DF, Johnston JA, Jonsdottir I, Stefansson K. A rare IL33 loss-of-function mutation reduces blood eosinophil counts and protects from asthma. PLoS Genet. 2017 Mar 8;13(3):e1006659 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1006659&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28273074&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F15%2F2023.12.14.23299853.atom) 19. 19.Alieva Yu.S., Furman E.G., Khuzina E.A., Abdrakhmanova S.T., Ponomareva M.S., Sheludko V.S., Sokolovsky V.L. Biomarkers of inflammation and control of bronchial asthma in children. Vopr. prakt. pediatr. (Clinical Practice in Pediatrics). 2023; 18(5): 13–20 20. 20.Contopoulos-Ioannidis D. G., Manoli E. N., Ioannidis J. P. A. Meta-analysis of the association of β2-adrenergic receptor polymorphisms with asthma phenotypes //Journal of Allergy and Clinical Immunology. – 2005. – T. 115. – No. 5. – C. 963–972. 21. 21.Furman, E. G. Bronchial asthma in children: inflammatory markers and the state of external respiratory function / E. G. Furman, I. P. Koryukina. - Perm : Book format, 2010. - 183 C. 22. 22.van der Plaat DA, de Jong K, Lahousse L, Faiz A, Vonk JM, van Diemen CC, Nedeljkovic I, Amin N, Obeidat M, van Duijn CM, Boezen HM, Postma DS. The Well-Known Gene HHIP and Novel Gene MECR Are Implicated in Small Airway Obstruction. Am J Respir Crit Care Med. 2016 Nov 15;194(10):1299–1302. 23. 23.Bártholo TP, Porto LC, Pozzan R, Nascimento A, Da Costa CH. Evaluation Of HHIP Polymorphisms And Their Relationship With Chronic Obstructive Pulmonary Disease Phenotypes. Int J Chron Obstruct Pulmon Dis. 2019 Oct 3;14:2267–2272 [1]: /embed/graphic-1.gif [2]: /embed/graphic-2.gif [3]: /embed/graphic-3.gif [4]: /embed/graphic-4.gif