- Research
- Open access
- Published:
Circulating miRNAs associate with historical childhood asthma hospitalization in different serum vitamin D groups
Respiratory Research volume 25, Article number: 118 (2024)
Abstract
Background
Vitamin D may help to alleviate asthma exacerbation because of its anti-inflammation effect, but the evidence is inconsistent in childhood asthma. MiRNAs are important mediators in asthma pathogenesis and also excellent non-invasive biomarkers. We hypothesized that circulating miRNAs are associated with asthma exacerbation and modified by vitamin D levels.
Methods
We sequenced baseline serum miRNAs from 461 participants in the Childhood Asthma Management Program (CAMP). Logistic regression was used to associate miRNA expression with asthma exacerbation through interaction analysis first and then stratified by vitamin D insufficient and sufficient groups. Microarray from lymphoblastoid B-cells (LCLs) treated by vitamin D or sham of 43 subjects in CAMP were used for validation in vitro. The function of miRNAs was associated with gene modules by weighted gene co-expression network analysis (WGCNA).
Results
We identified eleven miRNAs associated with asthma exacerbation with vitamin D effect modification. Of which, five were significant in vitamin D insufficient group and nine were significant in vitamin D sufficient group. Six miRNAs, including hsa-miR-143-3p, hsa-miR-192-5p, hsa-miR-151a-5p, hsa-miR-24-3p, hsa-miR-22-3p and hsa-miR-451a were significantly associated with gene modules of immune-related functions, implying miRNAs may mediate vitamin D effect on asthma exacerbation through immune pathways. In addition, hsa-miR-143-3p and hsa-miR-451a are potential predictors of childhood asthma exacerbation at different vitamin D levels.
Conclusions
miRNAs are potential mediators of asthma exacerbation and their effects are directly impacted by vitamin D levels.
Introduction
Asthma is the most common chronic disease in children characterized by airway inflammation in response to certain triggers, such as pollen, cold and other respiratory infections [1]. Asthma exacerbations are still the fourth leading cause of hospital stays among children aged 0–17 years with 100.7 stays per 100,000 in 2018 [2] and the direct costs of pediatric asthma were estimated at 6.31 billion in 2018 dollars [3].
Vitamin D is a fat-soluble secosteroid which can be synthesized by skin in the presence of ultraviolet B (UVB) radiation from sunlight or taken in from the supplementation [4]. Vitamin D is processed by the liver to form calcifediol (25(OH)D3) and then transported to kidney where 25(OH)D3 is converted to the active metabolite calcitriol (1,25(OH)2D3) [5]. Vitamin D is well-studied in regulating the level of calcium in the body and bone formation and also plays an important role in immune regulation and respiratory infections [6]. In our previous study, we associated 25-hydroxyvitamin D (25(OH)D) levels in serum at baseline with hospitalization or emergency department (ED) visits in the Childhood Asthma Management Program (CAMP). We confirmed a higher odds of hospitalization and ED visit for insufficient vitamin D status, but the mechanisms were still unclear [7].
MiRNAs are small non-coding RNAs with 18–25 nucleotides and act as post-transcriptional regulators in various biological processes and signaling pathways, such as cell apoptosis, immune and inflammation [8, 9]. MiRNAs also function in the asthma pathogenesis [10,11,12] and circulating miRNAs are good non-invasion biomarkers of asthma exacerbation, remission and ICS prognosis [13,14,15,16].
In this study, we investigated the association between miRNAs and asthma exacerbation as modified by serum vitamin D levels. We associated miRNA expression at baseline with historical asthma hospitalization in consideration of the interaction between miRNA and vitamin D first and then in vitamin D insufficient group and sufficient group separately. To investigate the potential pathways of the interaction effect that mediated by miRNAs, we performed microarray of lymphoblastoid B-cells treated by vitamin D or sham from 43 asthmatic subjects in CAMP to identified the differentially expressed genes (DEGs) and these DEGs were clustered into gene modules by WGCNA. These gene modules were enriched with functions by GO enrichment analysis and miRNAs with significant interaction effect were associated with these functional gene modules. We also demonstrate that these miRNAs can predict both historical and future hospitalizations at an excellent level.
Methods
Participant selection
We selected 461 subjects from the Childhood Asthma Management Program (CAMP) (Clinicaltrials.gov: NCT00000575 Date: October 27, 1999), which was a double-blind randomized controlled clinical trial [17, 18]. These participants provided serum samples at enrollment. The serum concentration of 25-hydroxyvitamin D3 was measured using a radioimmunoassay method in Dr. Bruce Hollis’ laboratory at the Medical University of South Carolina [19, 20] and then categorized into insufficient (≤ 30 ng/ml) and sufficient (> 30 ng/ml) [21]. The history of previous exacerbation was collected by questionnaires during the screening visits [18].
Small RNA sequencing and profiling
We isolated total RNA of 461 subjects from serum samples using Qiagen miRNeasy Serum/Plasma extraction kit and QIAcube automation. All samples were quantified using the Nanodrop spectrophotometer prior to plating with the RNA concentration 30.97 ± 20.97 ng/μl (see Additional file 1: Fig. S1). We built the small RNA-Seq libraries by Norgen Biotek Small RNA Library Prep Kit and sequenced on the Illumina NextSeq 500 platform at 51 bp single end reads. The sequencing data was deposited in the Gene Expression Omnibus with the accession number GSE134897 [13]. COMPSRA was employed to evaluate the read quality and trim adapters [22]. Sequencing reads with quality score lower than 20 were removed. The qualified reads were aligned to human genome hg38 by STAR (v2.7.10b) [23] and miRNAs were annotated by COMPSRA on the basis of miRbase [24].
Statistical analysis
Logistic regression was used to examine the association between miRNAs and the historical asthma hospitalization (Yes: Y = 1 and No: Y = 0). Age, sex, race and BMI were considered a priori as confounders in the multivariate logistic regression model. We investigated the interaction between miRNA and vitamin D groups and then examined the association in both the vitamin D insufficient group (Group = 1) and the sufficient group (Group = 0).
Validation via functional module
We isolated lymphoblastoid B-cells (LCLs) from blood of 43 asthmatic subjects enrolled in the Childhood Asthma Management Program (CAMP) [25]. The LCLs were then cultured in RPMI 1640 medium supplemented with 5% fetal bovine serum, and 1× Penicillin/Streptomycin/l-Glutamine. Afterward, the cells were treated with a sham (culture media) and 1 μM of 1,25-OH vitamin D for 72 h. RNA was extracted from the cells using the RNA Miniprep column purification system provided by Stratagene (La Jolla, CA). Total RNA samples were profiled on the Illumina HT12 V4 microarray. Quantile normalization was used to normalize the gene expression. Differentially expressed genes (DEGs) were identified with P ≤ 0.05 [adjusted by Benjamini–Hochberg (BH) Procedure]. Weighted Gene Correlation Network Analysis (WGCNA) v1.72.1 was used to identify co-modulated gene modules and these modules were also associated with miRNAs [26]. GO enrichment analysis of the gene modules were conducted using R package “clusterProfiler” [27] and p values were adjusted by Benjamini–Hochberg (BH) Procedure.
Validation via external dataset
We studied the dataset GSE106885 in Gene Expression Omnibus (GEO), which contains four asthma human bronchial epithelial cell samples and treated with calcitriol, poly I:C, both and sham (culture media) respectively [28]. The RNA-seq read files of samples treated with calcitriol and sham were downloaded from SRA (SRP124965) and taken as the external dataset. The raw reads were filtered by sequencing quality first than then aligned to human reference genome hg38 by STAR (v2.7.10b) [23]. The differentially expressed genes were calculated by edgeR (v3.40.0) [29] and p values were adjusted by Benjamini–Hochberg (BH) Procedure. The miRNA-target interactions comes from miRTarBase with strong evidence [30].
Prediction
The random forest model was used to predict the historical asthma hospitalization and future asthma hospitalization during the first year CAMP clinical trial by baseline miRNA expression level. The R package “randomForest” was employed to conduct the prediction [31] and the area under receiver operating characteristic curve (AUROC) was used for evaluating model performance.
Results
Baseline characteristics
A total of 461 samples were investigated including 139 subjects in the vitamin D insufficient (30 ng/ml) group and 322 subjects in the vitamin D sufficient (> 30 ng/ml) group. The baseline characteristics are shown in Table 1. Subjects in the vitamin D insufficient group had older age (9.12 ± 2.16 vs. 8.74 ± 2.12 yr) and higher BMI (18.67 ± 3.69 vs. 17.82 ± 3.15 kg/m2). Race was significantly different between vitamin d insufficient and sufficient groups (P < 0.01). The sex and primary outcome of child ever in hospital for asthma was not significantly different between the vitamin D insufficient and sufficient groups.
Significant miRNAs in the vitamin D interactions analysis
We identified 11 miRNAs that were associated with historical hospitalization in the interaction analysis (nominal P ≤ 0.05, FDR ≤ 0.1) (see Table 2). The effect modification by vitamin D was observed in both positive and negative directions. Hsa-miR-24-3p had the strongest protective effect (OR = 0.31; P = 0.005) and hsa-miR-22-3p had the strongest risk effect (OR = 2.08; P = 0.02). Hsa-miR-146b-5p was the most significant miRNA (OR = 0.32; P = 0.001) and had a strong protective effect modification by vitamin D.
Significant miRNAs in the stratified analysis
We examined the association between baseline miRNAs expression and dichotomized historical asthma hospitalization in a retrospective manner in both the vitamin D insufficient group and the sufficient group.
In the insufficient group, five miRNAs were nominally associated with risk of hospitalization in both univariate and multivariate analysis (nominal P ≤ 0.05, FDR ≤ 0.2) and the results were listed in Table 3. Three miRNAs (hsa-miR-24-3p, hsa-miR-423-3p and hsa-miR146b-5p) had the protective effect on asthma exacerbation, of which hsa-miR-24-3p had the smallest OR values (OR = 0.36; P = 0.005) in the multivariate models. Two miRNAs (hsa-miR-192-5p and hsa-miR-143-3p) had a higher risk effect on asthma exacerbation, of which hsa-miR-192-5p had the strongest OR values (OR = 1.52; P = 0.02) in the multivariate models. In addition, all the five miRNAs reported significant effect modification in the interaction analysis (see Table 2).
In the vitamin D sufficient group, nine miRNAs were associated with hospitalization in the univariate model and the association for each individual miRNA remains significant after adjusting for age, sex, race and BMI (nominal P ≤ 0.05, FDR ≤ 0.1) (see Table 4). Three miRNAs (hsa-miR-22-3p, hsa-miR-16-5p and hsa-miR-451a) were associated with a decreased risk of hospitalization, while six miRNAs (hsa-miR-125b-5p, hsa-miR-92b-3p, hsa-miR146b-5p, hsa-miR-423-5p, hsa-miR-125a-5p and hsa-miR-99b-5p) were associated with an increased risk of hospitalization. Hsa-miR-22-3p was associated with the lowest odds of hospitalization (OR = 0.62; P = 0.004) in the multivariate model and hsa-miR-99b-5p was associated with the highest odds of hospitalization (OR = 1.57; P = 0.05). Four miRNAs (hsa-miR-22-3p, hsa-miR-16-5p, hsa-miR-451a and hsa-miR-146b-5p) reported significant effect modification in the interaction analysis (see Table 2).
Hsa-miR-146b-5p was associated with hospitalization in both vitamin D insufficient (OR = 0.50; P = 0.02) and sufficient (OR = 1.54; P = 0.03) groups, as well as in the interaction analysis (OR = 0.32; P = 0.001), which implied potential function in asthma pathogenesis via vitamin D.
Validate the function of miRNAs in the interaction analysis in vitro
To validate whether the target genes of miRNAs that associated with hospitalization in the interaction analysis were functional at different vitamin D concentrations, we conducted microarray experiment from 43 participants in CAMP (see Additional file 1: Table S1). The lymphoblast-like B cells (LCLs) were isolated from blood and treated with 1,25-OH vitamin D or sham. We identified 811 genes that were differentially expressed between vitamin D treated LCLs and controls, of which 180 genes (22.19% of the differentially expressed genes) overlapped with the miRNA target genes (see Fig. 1A, B).
Based on the expression profiles of differentially expressed genes, we performed WGCNA analysis and identified 9 gene modules that were associated with at least one miRNA in the interaction analysis (see Fig. 1C). GO enrichment analysis was performed to annotate the function of these gene modules and we found that Module brown, Module red, Module blue and Module turquoise had significant biological processes (see Fig. 1D).
Module brown was enriched into leukocyte development related pathways (e.g. GO:1902106, negative regulation of leukocyte differentiation, FDR = 7.94 × 10–5) and negatively associated with five miRNAs, including hsa-miR-143-3p, hsa-miR-151a-5p, hsa-miR-192-5p, hsa-miR-22-3p and hsa-miR-24-3p. Of which, hsa-miR-143-3p (r = − 0.62; P = 2 × 10–10) had the largest absolute value of correlation coefficient. Module red was enriched into both vitamin D (e.g. GO: 0071295, cellular response to vitamin, FDR = 0.008) and leukocyte migration (e.g. GO:0050900, leukocyte migration, FDR = 0.009) pathways and associated with hsa-miR-22-3p, hsa-miR-24-3p and hsa-miR-451a. Of which, hsa-miR-24-3p (r = − 0.43; P = 3 × 10–5) had the strongest negative correlation. Module blue was enriched into both immune response (e.g. GO:0002366, leukocyte activation involved in immune response, FDR = 0.01) and cell adhesion (e.g. GO:0022409, positive regulation of cell–cell adhesion, FDR = 0.01) and associated with has-miR-192-5p and hsa-miR-451a. Module turquoise was also enriched into immune response (e.g. GO:0002764, immune response-regulating signaling pathway, FDR = 2.17 × 10–4) and associated with hsa-miR-143-3p and hsa-miR192-5p.
Validate the function of miRNAs in the interaction analysis through external dataset
To gain a more comprehensive understanding of the effect modification between miRNAs and vitamin D, we download the dataset GSE106885 from Gene Expression Omnibus (GEO) which contains asthma human bronchial epithelial cell samples treated with calcitriol (bioactive vitamin D) and sham (culture media). We identified 1702 genes that were differentially expressed between the calcitriol treated and sham samples, of which 454 genes (26.7% of differentially expressed genes) were overlapped with the miRNA target genes (see Additional file 1: Figs. S2, S3). These overlapped genes were enriched into the signaling pathways about cell apoptosis, cell proliferation and immune function, which suggested that the effect modification of vitamin D on miRNAs may vary in different cell types (see Additional file 1: Fig. S4).
Predict historical and future asthma hospitalization via baseline miRNAs expression
We tried to predict the past asthma hospitalization via baseline miRNA expression in vitamin D insufficient and sufficient groups first and then examined that whether these miRNAs were able to predict the risk of future asthma hospitalization during the first year CAMP clinical trial. Among the five significant miRNAs in the vitamin D insufficient group (see Table 3), hsa-miR-143-3p performed the prediction best with AUROC 0.8 for the past asthma hospitalization and AUROC 0.89 for the future asthma hospitalization (see Fig. 2A). Among the nine significant miRNAs in the vitamin D sufficient group (see Table 4), hsa-miR-451a predicted best with AUROC 0.78 for the past asthma hospitalization and AUROC 0.73 for the future asthma hospitalization (see Fig. 2B). Both hsa-miR-143-3p and hsa-miR-451a had significant effect modification on vitamin D (see Table 2).
Discussion
In our study, we identified eleven miRNAs in the interaction analysis, of which five miRNAs were associated with historical asthma hospitalization in the vitamin D insufficient group and nine miRNAs in the vitamin D sufficient group, suggesting miRNAs may be a pivotal vitamin D-related mediator of asthma severity in childhood asthma. These miRNAs were mainly associated with functional modules in immune system, such as leukocyte development, leukocyte migration, immune response, through microarray experiment in vitro. In addition, hsa-miR-143-3p and hsa-miR-451a could predict both historical and future hospitalization with good–excellent AUROC in vitamin D insufficient and sufficient groups.
Hsa-miR-146b-5p had a significant protective effect on childhood asthma in the vitamin D insufficient group and the significant risk effect in the vitamin D sufficient group. Bianchi and colleagues validated that hsa-miR-146b-5p was upregulated in T cells when treated with vitamin D [32]. Hsa-miR-146b-5p shared the same “seed” in sequence with hsa-miR-146a-5p which exhibited strong anti-inflammatory effect by down-regulating the expression of NF-κB [33, 34].
Hsa-miR-24-3p had the strongest protective effect in the vitamin D insufficient group and was associated with Module brown (leukocyte development) and Module red (leukocyte migration and vitamin D response) (see Fig. 1C, D). Lal et al. confirmed that miR-24-3p was consistently up-regulated in differentiated blood cells [35, 36]. Pua et al. reported that miR-24-3p could inhibit IL-4 production in T cells in vitro and alleviate allergic airway hypersensitivity inflammatory responses [37]. Xu et al. found that miR-24-3p could attenuate IL-1β induced chondrocyte injury associated with osteoarthritis [38]. MiR-24-3p also exerts anti-inflammatory function in LPS-stimulated macrophages by inhibiting pro-inflammatory cytokines [39].
Hsa-miR-22-3p had the strongest protective effect in the vitamin D sufficient group and was also associated with Module brown and Module red (see Fig. 1C, D). Lal et al. also reported that miR-22-3p was up-regulated when HL60 cells differentiating to monocytes using vitamin D3 [36]. Youn et al. confirmed the anti-inflammatory effect of miR-22-3p in macrophages [40]. Guo et al. validated that miR-22-3p was downregulated in the OVA-induced murine asthma model and LPS-induced bronchial epithelia cells, while overexpression of miR-22-3p ameliorated lung injury and inhibited epithelia cell injury [41]. Alvarez-Diaz and the group reported that miR-22-3p could be induced by 1,25(OH)2D3, which implied vitamin D had a strong effect on miRNA expression [42].
Hsa-miR-192-5p had the highest risk effect on asthma exacerbation in vitamin D insufficient group and was associated with Module brown, Module blue (immune response and cell adhesion) and Module turquoise (immune response) (see Fig. 1C, D). Zhang et al. showed miR-192-5p was associated with T follicular helper cell differentiation by regulating CXCR5 in childhood asthma [43] and Chu et al. found miR-192-5p could regulate the immune response by targeting Il-1RI [44].
The WGCNA reported that miRNAs interacting with vitamin D were associated with immune function in LCLs, while analysis of external dataset showed that these miRNAs were also associated with cell apoptosis and proliferation in bronchial epithelial cells. Both the validations suggest that vitamin D is a comprehensive regulator that may be involved in different functions in different cell types through interactions with miRNAs [45].
We attempted to use miRNAs to predict the past and future asthma hospitalization in both retrospective and prospective analysis through random forest method in different vitamin D groups. Hsa-miR-143-3p and hsa-miR-451a outperformed other miRNAs with the best AUROC, suggesting potential biomarkers of vitamin D related childhood asthma.
Our study had several limitations. First, we conducted our research in CAMP which is a clinical trial of childhood asthma, so the results should only be generalized to adult asthma with additional studies. Second, we lack a replication cohort for our results, but the validation for the WGCNA and enrichment analysis from microarray in vitro supports our main results.
Conclusions
In summary, miRNAs are potential mediators and effect modifiers of asthma exacerbation as directly impacted by vitamin D levels.
Availability of data and materials
The datasets generated and/or analysed during the current study are available in the Gene Expression Omnibus (GSE134897) repository.
Abbreviations
- CAMP:
-
Childhood Asthma Management Program
- UVB:
-
Ultraviolet B
- ED:
-
Emergency department
- DEGs:
-
Differentially expressed genes
- WGCNA:
-
Weighted Gene Correlation Network Analysis
- LCLs:
-
Lymphoblast-like B cells
- BH:
-
Benjamini–Hochberg
- AUROC:
-
Area under receiver operating characteristic curve
References
Martinez FD, Vercelli D. Asthma. Lancet. 2013;382:1360–72.
Healthcare cost and utilization project (HCUP). https://www.hcup-us.ahrq.gov/. Accessed 15 Dec 2022.
Perry R, Braileanu G, Palmer T, Stevens P. The economic burden of pediatric asthma in the United States: literature review of current evidence. Pharmacoeconomics. 2019;37:155–67.
Pfeffer PE, Hawrylowicz CM. Vitamin D and lung disease. Thorax. 2012;67:1018–20.
Bikle DD. Vitamin D metabolism, mechanism of action, and clinical applications. Chem Biol. 2014;21:319–29.
Ali NS, Nanji K. A review on the role of vitamin D in asthma. Cureus. 2017;9: e1288.
Brehm JM, Schuemann B, Fuhlbrigge AL, Hollis BW, Strunk RC, Zeiger RS, Weiss ST, Litonjua AA, Childhood Asthma Management Program Research G. Serum vitamin D levels and severe asthma exacerbations in the childhood asthma management program study. J Allergy Clin Immunol. 2010;126(52–58): e55.
Ha M, Kim VN. Regulation of microRNA biogenesis. Nat Rev Mol Cell Biol. 2014;15:509–24.
Gebert LFR, MacRae IJ. Regulation of microRNA function in animals. Nat Rev Mol Cell Biol. 2019;20:21–37.
Weidner J, Bartel S, Kilic A, Zissler UM, Renz H, Schwarze J, Schmidt-Weber CB, Maes T, Rebane A, Krauss-Etschmann S, Radinger M. Spotlight on microRNAs in allergy and asthma. Allergy. 2021;76:1661–78.
AlashkarAlhamwe B, Potaczek DP, Miethe S, Alhamdan F, Hintz L, Magomedov A, Garn H. Extracellular vesicles and asthma-more than just a co-existence. Int J Mol Sci. 2021;22:4984.
Grueso-Navarro E, Navarro P, Laserna-Mendieta EJ, Lucendo AJ, Arias-Gonzalez L. Blood-based biomarkers for eosinophilic esophagitis and concomitant atopic diseases: a look into the potential of extracellular vesicles. Int J Mol Sci. 2023;24:3669.
Li J, Panganiban R, Kho AT, McGeachie MJ, Farnam L, Chase RP, Weiss ST, Lu Q, Tantisira KG. Circulating microRNAs and treatment response in childhood asthma. Am J Respir Crit Care Med. 2020;202:65–72.
Kho AT, McGeachie MJ, Moore KG, Sylvia JM, Weiss ST, Tantisira KG. Circulating microRNAs and prediction of asthma exacerbation in childhood asthma. Respir Res. 2018;19:128.
Wang AL, Li J, Kho AT, McGeachie MJ, Tantisira KG. Enhancing the prediction of childhood asthma remission: Integrating clinical factors with microRNAs. J Allergy Clin Immunol. 2021;147(1093–1095): e1091.
Alhamdan F, Greulich T, Daviaud C, Marsh LM, Pedersen F, Tholken C, Pfefferle PI, Bahmer T, Potaczek DP, Tost J, Garn H. Identification of extracellular vesicle microRNA signatures specifically linked to inflammatory and metabolic mechanisms in obesity-associated low type-2 asthma. Allergy. 2023;78:2944–58.
Covar RA, Fuhlbrigge AL, Williams P, Kelly HW, The Childhood Asthma Management Program Research G. The childhood asthma management program (CAMP): contributions to the understanding of therapy and the natural history of childhood asthma. Curr Respir Care Rep. 2012;1:243–50.
Childhood Asthma Management Program Research Group. The childhood asthma management program (CAMP): design, rationale, and methods. Control Clin Trials. 1999;20:91–120.
Hollis BW, Kamerud JQ, Selvaag SR, Lorenz JD, Napoli JL. Determination of vitamin D status by radioimmunoassay with an 125I-labeled tracer. Clin Chem. 1993;39:529–33.
Bischoff-Ferrari HA, Giovannucci E, Willett WC, Dietrich T, Dawson-Hughes B. Estimation of optimal serum concentrations of 25-hydroxyvitamin D for multiple health outcomes. Am J Clin Nutr. 2006;84:18–28.
Holick MF. Vitamin D deficiency. N Engl J Med. 2007;357:266–81.
Li J, Kho AT, Chase RP, Pantano L, Farnam L, Amr SS, Tantisira KG. COMPSRA: a COMprehensive platform for small RNA-Seq data analysis. Sci Rep. 2020;10:4552.
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, Batut P, Chaisson M, Gingeras TR. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21.
Kozomara A, Griffiths-Jones S. miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res. 2011;39:D152-157.
Kho AT, Sharma S, Qiu W, Gaedigk R, Klanderman B, Niu S, Anderson C, Leeder JS, Weiss ST, Tantisira KG. Vitamin D related genes in lung development and asthma pathogenesis. BMC Med Genom. 2013;6:47.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinform. 2008;9:559.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16:284–7.
Boutaoui N, Puranik S, Zhang R, Wang T, Hui DH, Brehm J, Forno E, Chen W, Celedon JC. Epigenome-wide effects of vitamin D on asthma bronchial epithelial cells. Epigenetics. 2019;14:844–9.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.
Hsu SD, Lin FM, Wu WY, Liang C, Huang WC, Chan WL, Tsai WT, Chen GZ, Lee CJ, Chiu CM, et al. miRTarBase: a database curates experimentally validated microRNA-target interactions. Nucleic Acids Res. 2011;39:D163-169.
Liaw A, Wiener M. Classification and regression by randomForest. R news. 2002;2:18–22.
Bianchi N, Emming S, Zecca C, Monticelli S. Vitamin D and IFN-beta modulate the inflammatory gene expression program of primary human T lymphocytes. Front Immunol. 2020;11: 566781.
Jolivet G, Devinoy E, Fontaine ML, Houdebine LM. Structure of the gene encoding rabbit alpha s1-casein. Gene. 1992;113:257–62.
Barutta F, Bellini S, Mastrocola R, Bruno G, Gruden G. MicroRNA and microvascular complications of diabetes. Int J Endocrinol. 2018;2018:6890501.
Lal A, Navarro F, Maher CA, Maliszewski LE, Yan N, O’Day E, Chowdhury D, Dykxhoorn DM, Tsai P, Hofmann O, et al. miR-24 inhibits cell proliferation by targeting E2F2, MYC, and other cell-cycle genes via binding to “seedless” 3′UTR microRNA recognition elements. Mol Cell. 2009;35:610–25.
Lal A, Pan Y, Navarro F, Dykxhoorn DM, Moreau L, Meire E, Bentwich Z, Lieberman J, Chowdhury D. miR-24-mediated downregulation of H2AX suppresses DNA repair in terminally differentiated blood cells. Nat Struct Mol Biol. 2009;16:492–8.
Pua HH, Steiner DF, Patel S, Gonzalez JR, Ortiz-Carpena JF, Kageyama R, Chiou NT, Gallman A, de Kouchkovsky D, Jeker LT, et al. MicroRNAs 24 and 27 suppress allergic inflammation and target a network of regulators of T helper 2 cell-associated cytokine production. Immunity. 2016;44:821–32.
Xu J, Qian X, Ding R. MiR-24-3p attenuates IL-1beta-induced chondrocyte injury associated with osteoarthritis by targeting BCL2L12. J Orthop Surg Res. 2021;16:371.
Fordham JB, Naqvi AR, Nares S. miR-24 regulates macrophage polarization and plasticity. J Clin Cell Immunol. 2015;6:362.
Youn GS, Park JK, Lee CY, Jang JH, Yun SH, Kwon HY, Choi SY, Park J. MicroRNA-22 negatively regulates LPS-induced inflammatory responses by targeting HDAC6 in macrophages. BMB Rep. 2020;53:223–8.
Guo S, Chen R, Zhang L, Wu M, Wei Y, Dai W, Jiang Y, Kong X. microRNA-22-3p plays a protective role in a murine asthma model through the inhibition of the NLRP3-caspase-1-IL-1beta axis. Exp Physiol. 2021;106:1829–38.
Alvarez-Diaz S, Valle N, Ferrer-Mayorga G, Lombardia L, Herrera M, Dominguez O, Segura MF, Bonilla F, Hernando E, Munoz A. MicroRNA-22 is induced by vitamin D and contributes to its antiproliferative, antimigratory and gene regulatory effects in colon cancer cells. Hum Mol Genet. 2012;21:2157–65.
Zhang D, Wu Y, Sun G. miR-192 suppresses T follicular helper cell differentiation by targeting CXCR5 in childhood asthma. Scand J Clin Lab Invest. 2018;78:236–42.
Chu Q, Xu T. miR-192 targeting IL-1RI regulates the immune response in miiuy croaker after pathogen infection in vitro and in vivo. Fish Shellfish Immunol. 2016;54:537–43.
Janubova M, Zitnanova I. The effects of vitamin D on different types of cells. Steroids. 2023;202: 109350.
Acknowledgements
The authors thank all members of the Childhood Asthma Management Program (CAMP) research group and volunteers who have participated in this study.
Funding
This research was funded by the Shenzhen Science and Technology Program [Grant Number JCYJ20230807110302006]; the Shenzhen Key Laboratory of Chinese Medicine Active Substance Screening and Translational Research [Grant Number ZDSYS20220606100801003]; the Seventh Affiliated Hospital of Sun Yan-Sen University [Grant Number ZSQYBRJH0024]; NIH grants R01 HL127332, R01 HL129935, R01 HL162570, and R01 HL139634.
Author information
Authors and Affiliations
Contributions
Conceptualization: JL and KGT; methodology: XH; software: XH; validation: XH; formal analysis: XH and MJ; resources: KGT and STW; data curation: KGT and ATK; writing—original draft preparation, JL and XH; writing—review and editing, JL, XH, MJ, ATK, AT, AWL, HG, MJM, KGT and STW; visualization: XH and MJ; supervision: JL; project administration: JL and KGT; funding acquisition: JL, KGT. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Review Board (or Ethics Committee) of Mass General Brigham (Protocol #: 2017P001799 Date: June 14, 2022). Informed consent was obtained from all subjects involved in the study.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Table S1.
Demographic features of 43 subjects with asthma in LCLs microarray. Figure S1. RNA concentration of the subjects for sequencing in CAMP study. Figure S2. Volcano plots show the differentially expressed genes (DEGs) between the calcitriol treated and sham bronchial epithelial cell samples in GSE106885. Figure S3. Venn plot of target genes of 454 significant miRNAs in the interaction analysis and the differentially expressed genes in GSE106885. Figure S4. GO enrichment analysis of differentially expressed genes in dataset GSE106885.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Hong, X., Jiang, M., Kho, A.T. et al. Circulating miRNAs associate with historical childhood asthma hospitalization in different serum vitamin D groups. Respir Res 25, 118 (2024). https://0-doi-org.brum.beds.ac.uk/10.1186/s12931-024-02737-x
Received:
Accepted:
Published:
DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12931-024-02737-x