1 Department of Pediatrics, Kursk State Medical University, 305041 Kursk, Russia
2 Department of Biology, Medical Genetics and Ecology, Kursk State Medical University, 305041 Kursk, Russia
3 Laboratory of Genomic Research, Research Institute for Genetic and Molecular Epidemiology, Kursk State Medical University, 305041 Kursk, Russia
4 Department of Twin Research and Genetic Epidemiology, King’s College London, SE1 7EH London, UK
5 Laboratory of Statistical Genetics and Bioinformatics, Research Institute for Genetic and Molecular Epidemiology, Kursk State Medical University, 305041 Kursk, Russia
Abstract
Despite significant advances in identifying susceptibility genes for bronchial asthma (BA), the immunogenetic and other molecular mechanisms underlying disease risk and progression remain incompletely understood. This study aimed to investigate the association between two common single-nucleotide polymorphisms (SNPs), rs17293632 and rs2033784, located within the SMAD3 (SMAD family member 3) gene—previously identified in genome-wide association studies—and susceptibility to atopic asthma and allergen sensitivity in Russian children.
This study included DNA samples from 999 unrelated children from the Kursk region, comprising 526 patients with BA (316 boys and 210 girls) and 473 healthy children. SNP genotyping was performed using the MassARRAY-4 system.
Both SNPs were associated with an increased risk of BA under the additive genetic model (p ≤ 0.001). In the sex-stratified analysis, the rs17293632 (odds ratio (OR) = 1.79, 95% confidence interval (CI): 1.15–2.77; p = 9.0 × 10-4) and rs2033784 (OR = 1.80, 95% CI: 1.28–2.52; p = 6.0 × 10-4) variants were exclusively associated with asthma risk among girls, whereas no significant associations were observed in boys. The rs17293632T-rs2033784G haplotype was associated with an increased risk of asthma in girls (OR = 2.01, 95% CI: 1.38–2.92, p = 3.0 × 10-4). Notably, SMAD3 gene polymorphisms were associated with an increased risk of asthma exclusively among children living in urban environments (p ≤ 0.0001). Polymorphisms in the SMAD3 gene have also been linked to specific allergies, including those to horsehair, chicken eggs, beef, wheat, oatmeal, barley, and foxtail.
The present study of Russian children confirmed that the rs17293632 and rs2033784 polymorphisms of the SMAD3 gene are significant genetic markers associated with BA and allergies. However, an association between these SMAD3 variants and asthma was observed only among girls, marking the first time this association has been reported. Further research is necessary to determine whether the SMAD3 gene constitutes a viable therapeutic target for treating asthma and allergic diseases.
Keywords
- asthma
- hypersensitivity
- etiology
- Smad proteins
- genetic predisposition to disease
- single nucleotide polymorphism
Bronchial asthma (BA) is a common, clinically diverse, multifactorial disease resulting from complex interactions among multiple genetic and environmental factors that collectively contribute to its pathogenesis [1, 2, 3, 4]. The pathogenesis of BA involves a complex interplay of multiple immune cells—including mast cells, eosinophils, T lymphocytes, macrophages, neutrophils, and dendritic cells—and their inflammatory mediators, such as cytokines, chemokines, cysteinyl leukotrienes, histamine, nitric oxide, and prostaglandin D2 [5, 6]. Dysregulated T helper type 2 (Th2) inflammation, predominantly driven by cytokines IL-4, IL-5, and IL-13, underlies the core pathophysiological features of asthma: airway inflammation, hyperresponsiveness, and remodeling [6].
A causal role in asthma development within the general population has been
unequivocally established for several environmental factors, including
vaccination practices, reduced immune stimulation by bacterial products,
occupational exposures, and contact with particulate matter [7, 8]. Bronchial
asthma is a genetically heterogeneous disorder with a substantial heritable
component, as unequivocally demonstrated by familial aggregation, higher
concordance rates in monozygotic compared to dizygotic twins, and findings from
large-scale population-based genetic studies [1, 3, 4]. Despite significant
advances in identifying asthma susceptibility genes through genome-wide
association studies (GWAS) and candidate gene approaches, the immunogenetic
mechanisms underlying disease risk and progression remain incompletely understood
[2, 3, 9, 10, 11]. Multiple signaling pathways have been identified as contributing
factors in the pathogenesis and progression of asthma. The Smad signaling
pathway, in which the SMAD3 protein functions as an intracellular signal
transducer and transcriptional modulator, is activated by transforming growth
factor beta (TGF-
Numerous multi-ancestral GWAS have identified specific single nucleotide polymorphisms (SNPs) within the SMAD3 gene that are associated with the development and progression of bronchial asthma by modulating airway structure, immune responses, and inflammatory pathways [15, 16, 17]. Although GWAS have produced significant findings across diverse global populations, many SNPs identified within the SMAD3 gene require further validation in independent cohorts to confirm their clinical relevance and practical applicability. The present study aimed to investigate the association between two common SNPs, rs17293632 and rs2033784, located within the SMAD3 gene—previously identified through GWAS [11, 17]—and susceptibility to atopic asthma and sensitivity to allergens in pediatric populations from the Russian Federation.
The study included DNA samples from 999 unrelated children aged 1 to 16 years.
This cohort comprised 526 patients with allergic bronchial asthma who were
receiving inpatient treatment at the Kursk Regional Children’s Clinical Hospital,
and 473 unrelated children without any signs of bronchopulmonary or allergic
pathology at the time of sample collection, all residing in the Kursk region. The
group of patients with BA included children with intermittent, mild persistent,
moderate persistent, and severe persistent asthma. The diagnosis of bronchial
asthma was confirmed by qualified allergologists and pulmonologists in accordance
with national guidelines [18], based on a typical clinical presentation,
exclusion of other causes of bronchial obstruction, and findings from laboratory
and instrumental examinations, as described previously [19, 20, 21]. Clinical
examination was conducted following the standard protocol and included inspection
(skin color), palpation, percussion of the lung fields and heart, and
auscultation to detect abnormal breath sounds such as harsh or diminished
breathing and the presence of dry rales. The severity of bronchial obstruction
was clinically evaluated based on the number of choking attacks during the day,
the frequency of nocturnal symptoms, and the daily number of
Genetic studies were conducted at Research Institute of Genetic and Molecular
Epidemiology at Kursk State Medical University (Kursk, Russia). Genomic DNA was
isolated from peripheral venous blood samples of the study participants using the
classical phenol-chloroform extraction method. The GWAS Catalogue
(https://www.ebi.ac.uk/gwas/home) was used to identify SNPs of the SMAD3
gene associated with the isolated asthma phenotype in European populations.
Associations identified in mixed populations were excluded when selecting
polymorphisms. We selected SNPs associated with the following disease phenotypes
(trait names): “Asthma”, “Asthma (moderate or severe)”, “Asthma (childhood
onset)”, and “Atopic asthma”. As a result, we identified nine SNPs
(rs10152593, rs12592283, rs12912010, rs17293632, rs28617673, rs56062135,
rs56375023, rs72743461, and rs744910) with a minor allele frequency (MAF)
The sample size of patients for the study was determined through statistical
power calculations performed for each phase of the association analysis.
Statistical power calculations were performed using the Genetic Association Study
(GAS) power calculator
(https://csg.sph.umich.edu/abecasis/gas_power_calculator/). Assuming a
significance level of
The baseline and clinical characteristics of the study patients are summarized in Table 1. The group of patients with bronchial asthma was matched with the control group by gender (p = 0.72). The mean age of the BA patients was 8.1 years, compared to 8.2 years in the control group. The mean age at disease onset was 5.1 years. Concomitant allergic conditions—including food allergy, atopic dermatitis, allergic rhinitis, pollinosis, and acute exogenous urticaria—were observed in 402 children within the study group. The distribution of asthma severity was as follows: mild BA in 264 patients (50.2%), moderate BA in 184 patients (35%), and severe BA in 78 patients (14.8%). Among the study group, 214 children lived in urban areas, while 312 resided in rural areas.
| Characteristics | Children with BA | Control group | p-value |
| n = 526 | n = 474 | ||
| Age |
8.1 |
8.2 |
0.35 |
| Boys, n (%) | 316 (60) | 279 (59) | 0.72 |
| Girls, n (%) | 210 (40) | 195 (41) | |
| Place of residence (Urban/Countryside), n (%) | 214 (41) | 188 (39) | 0.03 |
| 312 (59) | 286 (61) | ||
| Concomitant pathology | |||
| Food allergy, n (%) | 306 (58.2) | - | - |
| Atopic dermatitis, n (%) | 12 (2.28) | - | - |
| Allergic rhinitis, n (%) | 123 (23.38) | - | - |
| Pollinosis, n (%) | 12 (2.28) | - | - |
| Insect allergy, n (%) | 64 (12.16) | - | - |
| Urticaria, n (%) | 76 (14.45) | - | - |
BA, bronchial asthma. Bold indicates statistically significant values
(p
The genotype frequencies of rs17293632 and rs2033784 polymorphisms in the
SMAD3 gene were in Hardy-Weinberg equilibrium in both patient groups
(p
| SNP ID | Genotype/Allele | Gene frequency (%) | OR (95% CI)1 | |
| Control group | Children with BA | p2 | ||
| Overall groups | ||||
| rs17293632 | C/C | 298 (63) | 255 (48.5) | 1.48 (1.18–1.85) |
| C/T | 153 (32.3) | 189 (36) | 6.0 × 10–4 | |
| T/T | 14 (4.7) | 28 (12.5) | ||
| rs2033784 | A/A | 254 (53.7) | 216 (41.1) | 1.39 (1.13–1.70) |
| A/G | 183 (38.7) | 214 (40.7) | 0.001 | |
| G/G | 26 (7.6) | 43 (18.2) | ||
| Boys | ||||
| rs17293632 | C/C | 162 (59) | 153 (53) | 1.23 (0.92–1.64) |
| C/T | 100 (37) | 118 (41) | 0.15 | |
| T/T | 11 (4) | 15 (6) | ||
| rs2033784 | A/A | 134 (49) | 129 (45) | 1.17 (0.90–1.52) |
| A/G | 119 (44) | 134 (47) | 0.25 | |
| G/G | 19 (7) | 25 (8) | ||
| Girls | ||||
| rs17293632 | C/C | 136 (70) | 102 (55) | 1.79 (1.15–2.77) |
| C/T | 53 (28) | 71 (38) | 9.0 × 10–4 | |
| T/T | 3 (2) | 13 (7) | ||
| rs2033784 | A/A | 120 (63) | 87 (47) | 1.80 (1.28–2.52) |
| A/G | 64 (33) | 80 (43) | 6.0 × 10–4 | |
| G/G | 7 (4) | 18 (10) | ||
1 Odds ratio and 95% confidence intervals.
2 p-level for SNP-BA association. Bold indicates statistically
significant values (p
SNP, single nucleotide polymorphism; OR, odds ratio; CI, confidence interval.
Notably, we also found that polymorphisms of the SMAD3 gene were
significantly associated with an increased risk of asthma in children living in
urban areas, whereas no such association was observed in those from rural regions
(Table 3). Specifically, the rs17293632 variant showed an odds ratio (OR) of 2.00
with a 95% confidence interval (CI) of 1.32–3.03 (p
| SNP ID | Genotype/Allele | Gene frequency (%) | OR (95% CI)1 | |
| Control group | Children with BA | p2 | ||
| Children living in rural regions | ||||
| rs17293632 | C/C | 177 (63.4) | 163 (58) | 1.29 (0.91–1.82) |
| C/T | 91 (32.6) | 108 (38.4) | 0.15 | |
| T/T | 11 (4) | 10 (3.6) | ||
| rs2033784 | A/A | 150 (54.1) | 136 (48.4) | 1.26 (0.90–1.76) |
| A/G | 110 (39.7) | 127 (45.2) | 0.17 | |
| G/G | 17 (6.2) | 18 (6.4) | ||
| Children living in urban areas | ||||
| rs17293632 | C/C | 121 (65) | 92 (48.2) | 2.00 (1.32–3.03) |
| C/T | 62 (33.3) | 81 (42.4) | p | |
| T/T | 3 (1.7) | 18 (9.4) | ||
| rs2033784 | A/A | 104 (55.9) | 80 (41.7) | 1.73 (1.25–2.39) |
| A/G | 73 (39.2) | 87 (45.3) | 7.0 × 10-4 | |
| G/G | 9 (4.9) | 25 (13) | ||
1 Odds ratio and 95% confidence intervals.
2 p-level for SNP-BA association.
* Data on the place of residence were available for 937 individuals. Bold
indicates statistically significant values (p
The results of the SMAD3 haplotype association analysis are
shown in Table 4. All three SMAD3 haplotypes had frequencies greater
than 5%. The rs17293632 and rs2033784 polymorphisms exhibited strong linkage
disequilibrium (D’
| Haplotype | SNP | Control group | Children with BA | OR (95% CI)1 | p2 | |
| rs17293632 | rs2033784 | |||||
| Overall groups | ||||||
| H1 | C | A | 0.746 | 0.682 | 1.00 | - |
| H2 | T | G | 0.195 | 0.259 | 1.48 (1.18–1.86) | 8.0 × 10-4 |
| H3 | C | G | 0.059 | 0.058 | 1.10 (0.74–1.62) | 0.64 |
| Boys | ||||||
| H1 | C | A | 0.711 | 0.681 | 1.00 | - |
| H2 | T | G | 0.224 | 0.258 | 1.22 (0.92–1.63) | 0.17 |
| H3 | C | G | 0.064 | 0.061 | 0.99 (0.60–1.63) | 0.97 |
| Girls | ||||||
| H1 | C | A | 0.797 | 0.685 | 1.00 | - |
| H2 | T | G | 0.153 | 0.261 | 2.01 (1.38–2.92) | 3.0 × 10-4 |
| H3 | C | G | 0.045 | 0.054 | 1.28 (0.69–2.39) | 0.44 |
1 Odds ratio and 95% confidence intervals.
2 p-level for haplotype-BA association. Bold indicates
statistically significant values (p
Data from skin prick tests were used to assess the association between
SMAD3 gene polymorphisms and allergy types in children with asthma. To
address the issue of multiple comparisons, when analyzing associations with
multiple allergens, we adjusted for multiple testing using the false discovery
rate (FDR) procedure via an online calculator
(https://www.sdmproject.com/utilities/?show=FDR). Table 5 presents a summary of
the observed associations between SMAD3 gene polymorphisms and various
allergy types in children with bronchial asthma. Both polymorphisms of
SMAD3 were associated with various types of allergies; however, the
nature and direction of these associations varied depending on the specific
allergic condition. In particular, the SNP rs2033784 was associated with
increased sensitivity to horse hair (FDR = 0.016), chicken eggs (FDR = 0.03),
wheat (FDR = 0.03), oatmeal (FDR = 0.04), and foxtail (FDR = 0.03). Nonetheless,
this variant demonstrated a protective effect against certain allergy types,
including barley (FDR = 0.05) and beef (FDR = 0.01). The polymorphism rs17293632
showed a significant association with increased sensitivity in children to wheat
(FDR = 0.03), timothy grass (FDR = 0.03), oatmeal (FDR = 0.03), and foxtail (FDR
= 0.03). Meanwhile, this variant was protective against beef allergy (FDR =
0.004). We found no significant associations between the rs17293632 and rs2033784
polymorphisms in the SMAD3 gene and the severity of bronchial asthma,
total IgE levels, or pulmonary function parameters (p
| SNP | Genotype | ОR (95% СI)1 | p2 | FDR |
| Horse hair | ||||
| rs17293632 | C/C-T/C | 0.87 (0.47–1.12) | 0.09 | 0.10 |
| T/T | ||||
| rs2033784 | A/A | 0.63 (0.47–0.85) | 0.002 | 0.02 |
| G/G | ||||
| Chicken eggs | ||||
| rs17293632 | C/C-T/C | 1.94 (0.89–3.16) | 0.12 | 0.12 |
| T/T | ||||
| rs2033784 | A/A-G/A | 2.26 (1.19–4.30) | 0.01 | 0.03 |
| G/G | ||||
| Wheat | ||||
| rs17293632 | C/C-T/C | 2.92 (1.35–6.31) | 0.007 | 0.03 |
| T/T | ||||
| rs2033784 | A/A-G/A | 2.14 (1.12–4.07) | 0.02 | 0.03 |
| G/G | ||||
| Barley | ||||
| rs17293632 | C/C-T/C | 0.69 (0.46–1.09) | 0.08 | 0.10 |
| T/T | ||||
| rs2033784 | A/A | 0.65 (0.43–0.98) | 0.04 | 0.05 |
| G/A-G/G | ||||
| Beef | ||||
| rs17293632 | C/C-T/T | 0.56 (0.37–0.83) | 0.004 | 0.02 |
| T/C | ||||
| rs2033784 | A/A | 0.52 (0.35–0.76) | 8.0 × 10–4 | 0.01 |
| G/A-G/G | ||||
| Timothy-grass | ||||
| rs17293632 | C/C-C/T | 2.41 (1.07–5.44) | 0.02 | 0.03 |
| T/T | ||||
| rs2033784 | A/A | 2.17 (0.95–6.04) | 0.09 | 0.10 |
| G/A-G/G | ||||
| Oatmeal | ||||
| rs17293632 | C/C-T/T | 1.65 (1.45–1.94) | 0.02 | 0.03 |
| C/T | ||||
| rs2033784 | A/A-G/G | 1.67 (1.47–1.97) | 0.03 | 0.04 |
| A/G | ||||
| Foxtail | ||||
| rs17293632 | C/C-C/T | 2.81 (1.24–6.35) | 0.009 | 0.03 |
| T/T | ||||
| rs2033784 | A/A-A/G | 2.04 (1.07–3.86) | 0.02 | 0.03 |
| G/G | ||||
1 OR – odds ratio and 95% confidence intervals;
2 p-level for SNP-allergy type association.
FDR, false discovery rate. Bold indicates statistically significant values
(p
The present study found, for the first time, that the rs17293632 and rs2033784 polymorphisms of the SMAD3 gene are significant markers of increased genetic susceptibility to bronchial asthma in Russian children in a sex-specific manner. However, both polymorphisms were found to be associated with asthma risk exclusively in girls. Furthermore, the observed associations were exclusively present among children residing in urban environments, while in children from rural areas, these genetic variants showed no significant influence on disease risk. Furthermore, we identified significant correlations between polymorphisms in the SMAD3 gene and multiple allergy phenotypes, demonstrating increased susceptibility to certain allergens while concurrently exhibiting reduced sensitivity to others.
Multiple polymorphisms within the SMAD3 gene have been identified as being associated with an increased risk of asthma and other allergic disorders, including atopic dermatitis, allergic rhinitis, and hay fever [15, 16, 17, 24, 25]. SMAD3 gene polymorphisms, including those examined in this study, are also linked to the development of various immune-mediated diseases such as ankylosing spondylitis, Crohn’s disease, psoriasis, ulcerative colitis, autoimmune thyroid disease, systemic lupus erythematosus, and type 1 diabetes mellitus [26, 27].
The availability of bioinformatics tools and genomic-transcriptomic databases, such as the GTEx Portal and the eQTLGen Consortium, which are derived from large population samples, enables the assessment of the functionality of genetic polymorphisms in relation to their impact on gene expression. This, in turn, facilitates the interpretation of genotypic correlations. The studied polymorphisms are located in intronic regions of the SMAD3 gene. The functional effects of the polymorphisms were evaluated using genome-transcriptome data available at the GTEx Portal (https://gtexportal.org). The asthma-associated alleles rs17293632-T and rs2033784-G are associated to decreased SMAD3 gene expression in thyroid tissue and esophageal mucosa, and, to a lesser extent, in the lung (only for the rs17293632-T allele). Additionally, a decreased SMAD3 gene expression in the carriers of the rs17293632-T and rs2033784-G alleles was observed in blood, according to data from the eQTLGen Consortium portal (https://www.eqtlgen.org/cis-eqtls.html). According to Ensembl data (https://www.ensembl.org/) on genes and regulation obtained from the ENCODE database, SNP rs17293632 participates in epigenetic regulatory functions within an enhancer region located in the lung and airway tissues. This activity is particularly prominent in bronchial epithelial cells, which are characterized by the presence of H3K27me3, H3K4me3, CTCF, and DNase I markers. It is also notable in lung fibroblasts, which exhibit markers such as H3K27ac, H3K27me3, H3K4me1, H3K4me3, CTCF, and DNase I.
We hypothesize that the association between SMAD3 gene polymorphisms
and susceptibility to asthma and sensitization to food allergens may result from
a disruption of TGF-
Our study has several limitations. The observed sex-specific associations between SMAD3 polymorphisms and asthma in our study are likely attributable to the relatively limited sample size, which is considerably smaller compared to the genome-wide association studies performed by large international consortia. Although the phenomenon of sexual dimorphism in the relationship between SNPs and predisposition to bronchial asthma has been identified in various studies [10, 37, 38], future genetic association studies investigating the SMAD3 gene should be conducted with substantially larger sample sizes. This will provide more robust and reliable estimates of the gene’s association with asthma development, particularly considering the observed sexual dimorphism in disease susceptibility. Since the associations between SMAD3 polymorphisms and various types of allergies were weak, further studies with larger sample sizes are necessary to confirm these findings. Furthermore, the lack of direct functional experimental data in our study limits the mechanistic explanations of the observed associations between SMAD3 polymorphisms and asthma or allergy. An interesting finding was that polymorphisms in the SMAD3 gene are associated with an increased risk of asthma in urban children, suggesting that specific environmental factors—particularly air pollution—may underlie this relationship. Toxicogenomic studies in animal models can help elucidate the underlying mechanisms and assess whether exposure to air pollutants modulates SMAD3 gene expression, thereby contributing to the development of asthma and allergies.
The present study of Russian children confirmed that the rs17293632 and
rs2033784 polymorphisms of the SMAD3 gene are significant genetic
markers associated with susceptibility to bronchial asthma. However, the
association of SMAD3 variants with asthma was observed only among girls,
marking the first time this finding has been reported. The associations between
polymorphisms and asthma, as well as various types of allergies that we
identified, confirm the important role of SMAD3 in the development of
immunopathological disorders underlying asthma and sensitization to multiple
allergens. Based on available genomic and transcriptomic data, there is strong
evidence that the asthma-associated loss-of-function alleles rs17293632 and
rs2033784 reduce SMAD3 gene expression, thereby disrupting
Smad3-mediated TGF-
BA, bronchial asthma; GWAS, genome-wide association study; Th2, T helper type 2; OR, odds ratio; CI, confidence intervals.
The data presented in this study are available upon reasonable request from the corresponding author.
AS and AP designed the research study. AP conceptualized and supervised the research. AS and AB examined the patients. AS, AP, and MF analyzed the data. AS, OB, and VP performed the laboratory genetic studies. AS drafted the manuscript. AP reviewed and edited the final manuscript. All authors contributed to the critical revision of the manuscript for important intellectual content. All authors read and approved the final manuscript. All authors have participated sufficiently in the work and agreed to be accountable for all aspects of the work.
The study was conducted in accordance with the guidelines of the Declaration of Helsinki. Informed consent for children’s participation in the study was obtained from their parents or legal guardians. The study protocol was approved by the regional ethics committee of Kursk State Medical University (protocol No. 4, dated 09.04.2018).
We thank all the children and their parents who participated in this study.
This research was conducted with the financial support of Kursk State Medical University.
The authors declare no conflict of interest. Given his role as the Editorial Board member, Alexey V Polonikov had no involvement in the peer-review of this article and has no access to information regarding its peer review. Full responsibility for the editorial process for this article was delegated to Gustavo Caetano-Anollés.
References
Publisher’s Note: IMR Press stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
