Front. Immunol. Frontiers in Immunology Front. Immunol. 1664-3224 Frontiers Media S.A. 10.3389/fimmu.2023.1121273 Immunology Original Research The causality between intestinal flora and allergic diseases: Insights from a bi-directional two-sample Mendelian randomization analysis Jin Qiubai 1 † Ren Feihong 1 2 † Dai Dan 1 Sun Nan 1 Qian Yiyun 3 Song Ping 1 * 1 Guang’anmen Hospital, China Academy of Chinese Medical Sciences, Beijing, China 2 Graduate School, Beijing University of Chinese Medicine, Beijing, China 3 Peking Union Medical College Hospital, Chinese Academy of Medical Science, Beijing, China

Edited by: Yongbo Kang, Shanxi Medical University, China

Reviewed by: Xinwei Huang, Tongji University, China; Niki Papapostolou, University General Hospital Attikon, Greece

*Correspondence: Ping Song, songping@vip.126.com

†These authors have contributed equally to this work and share first authorship

This article was submitted to Nutritional Immunology, a section of the journal Frontiers in Immunology

 09 03 2023 2023 14 1121273 11 12 2022 28 02 2023 Copyright © 2023 Jin, Ren, Dai, Sun, Qian and Song 2023 Jin, Ren, Dai, Sun, Qian and Song

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

 Background

Growing evidence shows a significant association between intestinal flora and allergic diseases, specifically atopic dermatitis (AD), allergic rhinitis (AR), and allergic asthma (AA). However, the causality has not yet been clarified.

 Objective

We conducted a bidirectional two-sample Mendelian randomization (TSMR) analysis to study the causal relationships between intestinal flora classification and AD, AR, or AA.

 Materials and methods

We obtained summary data of intestinal flora, AD, AR, and AA from a genome-wide association research. The inverse-variance weighted method is the primary method for analyzing causality in the TSMR analysis. Several sensitivity analyses were conducted to examine the stability of TSMR results. Reverse TSMR analysis was also performed to assess whether there was a reverse causality.

 Results

A total of 7 bacterial taxa associated with AD, AR, and AA were identified by the current TSMR analysis. Specifically, the genus Dialister(P=0.034)and genus Prevotella(P=0.047)were associated with a higher risk of AD, whereas class Coriobacteriia (P=0.034) and its child taxon, order Coriobacteriales (P=0.034) and family Coriobacteriaceae (P=0.034), all had a protective effect on AR. In addition, the family Victivallaceae (P=0.019) was identified as a risk factor for AR. We also noticed a positive association between the genus Holdemanella (P=0.046) and AA. The reverse TSMR analysis didn’t suggest any evidence of reverse causality from allergic diseases to the intestinal flora.

 Conclusion

We confirmed the causal relationship between intestinal flora and allergic diseases and provided an innovative perspective for research on allergic diseases: targeted regulation of dysregulation of specific bacterial taxa to prevent and treat AD, AR, and AA.

 Mendelian randomization intestinal flora atopic dermatitis allergic rhinitis allergic asthma causality Introduction

The incidence of atopic dermatitis (AD), allergic rhinitis (AR), and allergic asthma (AA) have dramatically increased over the past three decades, resulting in a considerable burden to society (1). AD is a common inflammatory skin disease characterized by recurrent eczematous lesions. AR is a non-infectious inflammatory disease of the nasal mucosa, with typical symptoms of paroxysmal sneezing, watery nasal discharge, nasal itching and congestion. AA is an airway inflammatory disease that results in recurring wheezing, chest tightness, shortness of breath, and mucus production. Currently, there is no effective radical treatment for AD, AR, or AA. Glucocorticoids and antihistamines are commonly used for treating these three allergic diseases; but symptoms tend to rebound after drug withdrawal (2). AD, AR and AA are associated with genetic, dietary, and environmental factors (such as air pollution and exposure to allergens); however, the underlying causes remain unclear (3). Therefore, there is an urgent need to identify potential causal risk factors for AD, AR and AA.

The intestinal flora is a dynamic ecosystem known as the ‘second genome’ (4). Intestinal dysbiosis can cause metabolic disorders of intestinal microorganisms and further lead to immune dysfunction (5). With an in-depth study of the gut-skin and gut-lung axes, more attention has been paid to the effect of intestinal flora on the skin and respiratory tract (6–8). Several observational studies have reported that the abundance of certain intestinal flora changes significantly in patients with AD, AR, and AA compared to healthy individuals, indicating a potential correlation between intestinal flora and these three allergic diseases (9–11). Intestinal flora can regulate adaptive immunity by maintaining the balance between effector T cells (Th1, Th2, and Th17) and regulatory T cells, which may explain the effect of intestinal flora on AD, AR and AA (12, 13).

However, owing to the lack of evidence from randomized controlled trials, it is unclear whether there is a definite causality between intestinal flora and these three allergic diseases. As the gold standard for causal inference in epidemiological studies, randomized controlled trials are sometimes difficult to conduct because of ethical limitations and high costs. Two sample Mendelian randomization (TSMR) is an effective alternative (14). Genome-wide association studies (GWAS) have made great contributions to the identification of genetic variants related to diseases, mainly single nucleotide polymorphism (SNP), which can increase our understanding of the genetic basis of many complex traits in human diseases (15). TSMR uses genetic variation related to exposure as an alternative indicator of exposure to study the causality between exposure and outcome (16). The selected SNPs are also called instrumental variances (IVs). TSMR simulates randomization based on the random distribution of genetic variants during gametogenesis, conceptually similar to randomized controlled trials (17). Since these genetic variants precede diseases progression and are independent of lifestyle and environmental factors after birth, TSMR can minimize the influence of confounding factors and reverse causality (16).

In this study, we used the latest available GWAS database published in 2021 (18) for TSMR analysis to investigate the possible causality between intestinal flora and AD, AR, and AA, to provide an innovative perspective for the research of allergic diseases: targeted regulation of specific bacterial taxa to prevent and treat AD, AR, and AA.

 Materials and methods Study design

TSMR was used to analyze the causal relationship between intestinal bacterial taxa and allergic diseases (AD, AR, and AA). The overall design of this study is shown in Figure 1 . To obtain reliable results, three hypotheses need to be satisfied when performing TSMR analysis (1): there is a strong correlation between genetic variants and exposure factors (2); there is no correlation between genetic variants and confounders; and (3) genetic variants can only affect the outcome through exposure factors, but not through other methods, that is, horizontal pleiotropy is not allowed ( Figure 1 ). Genetic variants that satisfy these three hypotheses can be included in TSMR analysis as instrumental variables (16).

Overview of present MR analyses and assumptions. AD, atopic dermatitis; AR, allergic rhinitis; AA, allergic asthma; SNP, single nucleotide polymorphism.

 Data sources and selection of instrumental variables

Summary statistics of the intestinal flora were obtained from a large-scale GWAS study by the MiBioGen consortium (https://mibiogen.gcc.rug.nl), involving 18340 European ethnic participants from 11 countries with 122,110 loci of variation (18). We screened the IVs of intestinal bacterial taxa at five levels (phylum, class, order, family, and genus) from this GWAS. Fifteen bacterial taxa without specific species names were eliminated. The GWAS statistics for AD, AR, and AA were obtained from the data released by FinnGen Research (https://r7.finngen.fi/) in July 2022. The diagnostic criteria of AD were based on ICD-8, ICD-9 and ICD-10 standards and the GWAS statistics contain 16,383,295 loci of variations from 10,277 cases and 278,795 controls. The diagnostic criteria of AR were based on ICD-9, and ICD-10 standards, and the GWAS statistics of AR contain 16383313 loci of variation from 8430 cases and 298829 controls. The diagnostic criteria of AA were based on ICD-10 standards and the GWAS statistics of AR contain 16383313 loci of variation from 8430 cases and 298829 controls.

To obtain more complete results, we used a genome-wide significance threshold (5 × 10-8) and a locus-wide significance threshold (1 × 10-5), respectively to screen SNPs related to exposure (19, 20). Linkage disequilibrium analysis was performed to satisfy TSMR’s hypothesis 1. The linkage disequilibrium correlation coefficient was set to r2<0.001 and clumping window >10000kb to ensure no linkage disequilibrium among the included IVs. To avoid horizontal pleiotropy, IVs associated with risk factors for allergic diseases were excluded using PhenoScanner V2 (21). Palindromic and incompatible SNPs were excluded when harmonizing the statistics of exposure and outcome, and SNPs related to exposure that could not be matched in the GWAS outcome statistics were excluded. To avoid the influence of weak instrument bias on causal inference, we used the formula F=β2 exposure/SE2 exposure to calculate the strength of the IVs (16, 20, 22) and eliminate IVs with F < 10 (23).

 Statistical analysis

TSMR was conducted to analyse the causality between bacterial taxa and AD, AR, and AA. In the absence of horizontal pleiotropy, the inverse-variance weighted (IVW) method can be the primary method for analyzing causality in TSMR analysis (24). Before that, we implemented the Cochrane’s Q test to evaluate the heterogeneity between the IVs. If heterogeneity was detected(P<0.05), the random-effects IVW model could provide a more conservative estimate; otherwise, the fixed-effect IVW model would be used (25). Other methods of TSMR analysis, including the weighted median estimator (WM) and MR-Egger regression (26), can supplement the IVW method and provide wider confidence intervals (27). These three TSMR methods for causal inference have their model assumptions. The IVW method is suitable for situations where horizontal pleiotropy does not exist (24); the WM method assumes that less than 50% of IVs have horizontal pleiotropy (28). The MR-Egger regression assumes that more than 50% of IVs are affected by horizontal pleiotropy (26).

If the result of the TSMR analysis was nominally significant (P < 0.05), we considered that there might be a causal relationship between the flora and outcome (29). If the causality between bacterial taxa and outcome is identified as significant by two or more TSMR methods, the result is considered robust (5).

The existence of horizontal pleiotropy may challenge the second TSMR hypothesis; therefore, we adopted various methods to monitor possible horizontal pleiotropy. Specifically, the p-value of the MR-Egger intercept test and MR pleiotropy residual sum and outlier (MR-PRESSO) global test can be used to assess the existence of horizontal pleiotropy, and P < 0.05 was considered statistically significant (5, 30). The MR-PRESSO outlier test can adjust horizontal pleiotropy by detecting and removing outliers (31), and the number of distributions in the MR-PRESSO analysis was set to 1000 (17).

Additionally, we conducted a leave-one-out sensitivity analysis of the identified significant results to determine whether the causal relationship of the TSMR analysis was caused by a single SNP (32). Finally, a reverse TSMR analysis was performed between allergic diseases (AD, AR and AA) and the identified significant bacterial taxa using positive TSMR analysis to examine whether a reverse causal association existed. The reverse TSMR procedure was the same as that for the above TSMR analysis. TSMR analyses were performed using the ‘TwoSampleMR’ (version 0.5.6) in R software (version 4.2.1).

 Ethical approval

The GWAS data used in this study were public de-identified data. The ethics committee approved these data; therefore, there was no need for additional ethical approval.

 Results

IVs were screened according to the conditions described above. The details of the SNPs, that were eventually included in the TSMR analysis of intestinal flora and allergic diseases, are presented in Supplementary Table 1 . After harmonization, the number of SNP involved in each pair of bacterial taxa and allergic diseases was more than three. The F-statistics of all SNPs were greater than ten, indicating that there are no weak IVs. Moreover, it should be noted that there is an inclusive relationship between intestinal flora classifications. Thus, the SNPs included in the class and their relevant order may overlap heavily. For example, SNPs of the order Coriobacteriales, class Coriobacteriia, and family Coriobacteriaceae.

Results of the TSMR analysis (locus-wide significance, <italic>P</italic><1×10<sup>-5</sup>)

The causal relationship between each pair of bacterial taxa and allergic disease was analyzed using the three TSMR methods ( Supplementary Table 2 ). Twenty-five potential causal associations between bacterial traits and allergic diseases were identified using one or more TSMR methods ( Figure 2 ). Among them, two bacterial taxa related to AD, four bacterial taxa associated with AR and one bacterial taxon related to AA were cross-validated using the IVW and WM methods ( Table 1 and Figure 3 ). We mainly focused on these seven relatively stable causal relationships.

Causal analysis of intestinal flora and allergic diseases (locus-wide significance, P<1×10-5). MR-PRESSO, Mendelian Randomization Pleiotropy Residual Sum and Outlier; IVW, inverse-variance weighted method; WM, weighted median estimator; AD, atopic dermatitis; AR, allergic rhinitis; AA, allergic asthma.

Summary of causality between intestinal flora and AD, AR or AA (P<1×10-5).

Human gut microbiota Nsnps Traits Method OR 95%CI P-value
Genus Dialister 11 AD WM 0.812 0.664-0.994 0.043
IVW 0.834 0.714-0.987 0.034
Genus Prevotella 11 AD WM 0.872 0.783-0.970 0.012
IVW 0.924 0.854-0.999 0.047
Class Coriobacteriia 13 AD WM 0.752 0.582-0.972 0.030
IVW 0.789 0.634-0.982 0.034
Family Coriobacteriaceae 13 AR WM 0.752 0.583-0.970 0.028
IVW 0.789 0.634-0.982 0.034
Family Victivallaceae 11 AR WM 1.120 1.001-1.252 0.047
IVW 1.107 1.017-1.205 0.019
Order Coriobacteriales 13 AR WM 0.752 0.574-0.985 0.039
IVW 0.789 0.634-0.982 0.034
Genus Holdemanella 11 AA WM 1.175 1.018-1.355 0.028
IVW 1.124 1.002-1.261 0.046

AD, atopic dermatitis; AR, allergic rhinitis; AA, allergic asthma; Nsnp, number of single nucleotide polymorphism; WM, weighted median estimator; IVW, inverse-variance weighted method; OR, odd ratio; CI, confidence interval.

Forest plot of the causality between cross-validated 7 bacterial taxa with the risks of AD, AR, or AA. IVW, inverse-variance weighted method; WM, weighted median estimator; AD, atopic dermatitis; AR, allergic rhinitis; AA, allergic asthma.

We also performed the leave-one-out sensitivity analysis for the identified significant bacterial taxa, and the results further validated the robustness of our results ( Supplement Figure 1 ). In the absence of heterogeneity, horizontal pleiotropy, and outliers, the results of TSMR analysis are credible.

AD

Seven causal associations from bacterial taxa to AD were identified by the IVW method. Considering the cross-validation, the results of the two bacterial taxa remained stable. In specific, our TSMR analysis found that genus Dialister (OR: 0.839, 95% confidence interval (CI): 0.714-0.987, P=0.034) and genus Prevotella (OR: 0.924, 95% CI: 0.854-0.999, P=0.047) were associated with a higher risk of AD. In the sensitivity analysis, Cochrane’s Q test did not suggest evidence of heterogeneity in the genus Dialister (P=0.214) and genus Prevotella (P=0.672) ( Supplementary Table 3 ). The MR-Egger intercept test observed no horizontal pleiotropy in the genus Dialister (P=0.188) and genus Prevotella (P=0.914) ( Supplementary Table 4 ). Similarly, MR-PRESSO global test didn’t detect any horizontal pleiotropy in the genus Dialister (P=0.231) and genus Prevotella (P=0.667) ( Supplementary Table 5 ). For the MR-PRESSO outlier test, no outlier was found in the genus Dialister and genus Prevotella ( Supplementary Table 5 ).

 AR

Eleven causal relationships from bacterial taxa to AR were identified by the IVW method. Considering the cross-validation, the results of the four bacterial taxa remained stable. Specifically, class Coriobacteriia (OR: 0.789, 95% CI: 0.634-0.982, P=0.034) and its child taxon, order Coriobacteriales and family Coriobacteriaceae, all had a protective effect on AR. On the contrary, the family Victivallaceae (OR: 1.107, 95% CI: 1.017-1.205, P=0.019) was associated with a higher risk of AR. In the sensitivity analysis, Cochrane’s Q test did not suggest evidence of heterogeneity in class Coriobacteriia (P=0.071), family Coriobacteriaceae (P=0.071), family Victivallaceae(P=0.938) and order Coriobacteriales(P=0.071) ( Supplementary Table 3 ). No horizontal pleiotropy was observed by MR-Egger intercept test in class Coriobacteriia (P=0.609), family Coriobacteriaceae (P=0.609), family Victivallaceae(P=0.507) and order Coriobacteriales(P=0.609) ( Supplementary Table 4 ). Similarly, MR-PRESSO global test didn’t detect any horizontal pleiotropy in class Coriobacteriia (P=0.081), family Coriobacteriaceae (P=0.081), family Victivallaceae(P=0.938) and order Coriobacteriales(P=0.949) ( Supplementary Table 5 ). For the MR-PRESSO outlier test, no outlier was found in class Coriobacteriia, family Coriobacteriaceae, family Victivallaceae and order Coriobacteriales( Supplementary Table 5 ).

 AA

Three causal associations from bacterial taxa to AA were identified by the IVW method. Considering the cross-validation, only one bacterial taxon remained stable. Specifically, our TSMR analysis found that genus Holdemanella (OR: 1.124, 95% CI: 1.002-1.261, P=0.046). In the sensitivity analysis, Cochrane’s Q test did not suggest evidence of heterogeneity in the genus Holdemanella (P=0.198) ( Supplementary Table 3 ). No horizontal pleiotropy was observed by the MR-Egger intercept test in the genus Holdemanella (P=0.886) ( Supplementary Table 4 ). Similarly, MR-PRESSO global test didn’t detect any horizontal pleiotropy in the genus Holdemanella (P=0.227) ( Supplementary Table 5 ). For the MR-PRESSO outlier test, no outlier was found in the genus Holdemanella( Supplementary Table 5 ).

 Results of the TSMR analysis (genome-wide significance threshold, <italic>P</italic><5 × 10<sup>-8</sup>)

In the TSMR analysis of intestinal flora as a whole and allergic diseases, the IVW, WM, and MR-Egger regression methods did not find any significant causal associations. In the sensitivity analysis, Cochrane’s Q test did not suggest evidence of heterogeneity, MR-Egger intercept test and MR-PRESSO global test didn’t detect any horizontal pleiotropy, and the MR-PRESSO outlier test didn’t find any outliers. All the results are shown in Supplementary Table 6 .

 Reverse TSMR analysis

The results of reverse TSMR analysis are presented in Supplementary Table 7 . Considering cross-validation, we did not find any reverse causal relationships between the intestinal flora classification shown in Table 1 and allergic diseases.

 Discussion

Using large-scale GWAS statistics, seven bacterial traits associated with AD, AR, and AA were identified by the current TSMR analysis ( Figure 4 ). In accordance with some prospective observational studies (10, 33) and animal experiments (34, 35), our TSMR study also revealed that the genus Dialister and genus Prevotella may be protective factors for AD. Th2-skewed and Th17-skewed immune dysregulation are some of AD’s most significant pathogenesis mechanisms (36), whereas Treg cells can inhibit Th2 and Th17 allergic inflammation and restore immune tolerance (37). Furthermore, the genus Dialister is a propionate producer in the intestinal tract (38), and genus Prevotella can break down fibers and produce propionate and butyrate (39). Several studies have discovered that the content of short-chain fatty acids (SCFAs) such as propionate and butyrate in patients with AD is significantly lower than that in healthy individuals (40, 41). As important SCFAs, propionate and butyrate can inhibit the Th2-skewed and Th17-skewed inflammation in AD. Specifically, both can inhibit histone deacetylase (42) and induce the differentiation of peripheral CD4+T cells to Treg cells, thus producing anti-inflammatory cytokine IL-10 and inhibiting the function of Th2 and Th17 cells (43). In addition, some studies have shown that butyrate can inhibit the release of histamine and other inflammatory mediators from mast cells by inhibiting the interaction between Immunoglobulin E and mast cells (44). Therefore, we speculate that the protective effect of these two bacterial taxa on AD might be related to SCFAs, especially propionate and butyrate.

Bacterial taxa associated with AD, AR or AA identified by the current MR analysis. The blue arrow indicates that the bacterial taxa is the protective factor of the outcome and the red arrow indicates that the bacterial taxa is the risk factor of the outcome. MR, Mendelian randomization; IVW, inverse-variance weighted method; WM, weighted median estimator.

The Genus Holdemanella was considered as risk factor for AA based on the TSMR analysis. Reportedly, the abundance of the genus Holdemanella is negatively correlated with the propionate content in patients with diabetes and cognitive impairment (45). Thus, we speculate that the genus Holdemanella may promote Th2 inflammation in AA by affecting SCFAs such as propionate (46). However, the exact mechanism is still unclear and it is necessary to further study the possible role of the genus Holdemanella.

Class Coriobacteriia and its child taxa, order Coriobacteriales, and family Coriobacteriaceae, all have negative effects on AR, whereas, family Victivallaceae is the risk factor for AR. The relationship between these intestinal florae and allergic diseases was cross-validated using two IVW and WM methods. However, the function of this bacteria is poorly understood. Currently, there are no studies on the relationship between these bacterial taxa and allergic diseases, and the specific flora of AR and AA were reported for the first time in this study. Therefore, our study may provide a new perspective for mechanistic research on AR and AA.

In addition to the seven stable causal associations above, the IVW method yielded several interesting results, which are supported by previous studies. We also discovered that the family Bacteroidaceae and its child taxon genus Bacteroides were both risk factors for AD. Studies have reported that the abundance of the genus Bacteroides in patients with AD is significantly higher than that in healthy people (47), and its proportion is positively correlated with the severity of AD symptoms (33). Lipopolysaccharide, the metabolite of the family Bacteroidaceae and genus Bacteroides, can promote Th2 inflammation in AD (48). Moreover, the genus Bacteroides and genus Prevotella share a common ancestor but have inhibitory effects on each other (49). Genus Bacteroides is dominant in people who consume protein and animal fat in their main diet. In contrast, genus Prevotella is predominant in people who take fruits and vegetables as their main diet (50). The research of Nosrati et al. showed that vegetable intake could improve AD symptoms (51). This suggests that the relationship between intestinal flora and AD can be studied from the perspective of diet in the future. After all, it is much easier to adjust eating habits than to change genetic or environmental factors.

Butyrate, a metabolite of gut microbiota, may be an important connector between gut microbiota and allergic diseases (52). Genus Subdoligranulum (53) and genus Collinsella (54), both butyrate producers, were also determined to be protective factors of AR and AA by the IVW method. However, there are some differences in the role of butyrate in patients with AR and AA. Specifically, Th2-skewed and Th17-skewed immune dysregulation are important pathogenic mechanisms of AR (55). Similarly, the genus Subdoligranulum may also inhibit Th2 and Th17-mediated inflammation in AR by producing butyrate. In addition, IL-4, an important cytokine in Th2 inflammation, can impair the airway epithelial barrier function in AR patients. In this case, butyrate can improve the function of the airway epithelial barrier (56), which may explain the protective effect of genus Subdoligranulum on AR. Moreover, the inflow of eosinophils into the lung parenchyma is a hallmark of AA (46). Butyrate flows into eosinophils through monocarboxylate transporters and promotes apoptosis (57). In addition, type 2 innate lymphoid cells (ILC2) can promote T2 immunity in AA (46). Butyrate can inhibit the release of type 2 immune factors such as IL-5 and IL-13 by ILC2 cells, which can reduce the inflammatory response of AA (58).

This study has several advantages. Firstly, this is the first bi-directional TSMR study to reveal the causal relationship between intestinal flora and allergic diseases (AD, AR, and AA), which is not disturbed by confounding factors or reverse causality. Second, we set strict conditions for the screening of instrumental variables, and only when more than two TSMR methods identify the causal relationship can it be considered conceivable. Thirdly, we provided evidence for research on the intestinal-skin and intestinal-lung axes from a genetic perspective. Seven bacterial taxa associated with AD, AR, and AA were identified using TSMR analysis. These identified significant bacterial taxa could serve as candidate microbiome interventions in future clinical trials of allergic diseases. Meanwhile, our findings may provide an innovative perspective for research on allergic diseases: targeted regulation of specific bacterial taxa such as supplementing beneficial bacteria and inhibiting the growth of harmful bacteria to prevent and treat AD, AR and AA.

There are also some limitations of this study. Firstly, the number of instrumental variables involved in GWAS statistics of intestinal flora is limited, and there are no data available at the species level. Secondly, we could not determine whether there were overlapping participants in the GWAS data of the exposures and outcomes involved in this study. Thirdly, demographic data were lacking in the original research; therefore, we could not perform subgroup analysis on factors such as gender. Our results need to be verified by further clinical and basic research. In future study, we will increase the sample size and more accurately explore the relationship between intestinal flora and allergic diseases at the species level.

 Data availability statement

The original contributions presented in the study are included in the article/ Supplementary Material . Further inquiries can be directed to the corresponding author.

 Author contributions

PS and QJ designed the study. NS and YQ collected the data. QJ, FR, DD, and NS performed the computations and manuscript writing. All authors contributed to the article and approved the submitted version.

 Funding

The laboratory of the authors benefits from ongoing support from National Natural Science Foundation of China (82074448) and the Scientific and technological innovation project of China Academy of Chinese Medical Sciences (CI2021A02301).

 Acknowledgments

We would like to appreciate the participants and researchers of the FinnGen study. We also thank the MiBioGen consortium for providing the GWAS data of intestinal flora. In addition, we want to acknowledge Figdraw (www.figdraw.com) for their assistance in Figures 1 and 4 .

 Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1121273/full#supplementary-material

Abbreviations

AD, atopic dermatitis; AR, allergic rhinitis; AA, allergic asthma; TSMR, two sample mendelian randomization; GWAS, genome-wide association study; IVW, inverse-variance weighted; SNP, single nucleotide polymorphism; IVs, instrumental variances; WM, weighted median estimator; MR-PRESSO, mendelian randomization pleiotropy Residual Sum and Outlier; OR, odd ratio; CI, confidence interval.

 References Nwaru BI Virtanen SM . Allergenic food introduction and childhood risk of allergic or autoimmune disease. Jama (2017) 317(1):86. doi: 10.1001/jama.2016.18329 Cook Q Argenio K Lovinsky-Desir S . The impact of environmental injustice and social determinants of health on the role of air pollution in asthma and allergic disease in the united states. J Allergy Clin Immunol (2021) 148(5):1089–101.e5. doi: 10.1016/j.jaci.2021.09.018 Boutin RCT Sbihi H Dsouza M Malhotra R Petersen C Dai D . Mining the infant gut microbiota for therapeutic targets against atopic disease. Allergy (2020) 75(8):2065–68. doi: 10.1111/all.14244 Grice E A Segre JA . The human microbiome: our second genome. Annu Rev Genomics Hum Genet (2012) 13:151–70. doi: 10.1146/annurev-genom-090711-163814 Ni JJ Xu Q Yan SS Han BX Zhang H Wei XT . Gut microbiota and psychiatric disorders: A two-sample mendelian randomization study. Front Microbiol (2021) 12:737197. doi: 10.3389/fmicb.2021.737197 Budden KF Gellatly SL Wood DL Cooper MA Morrison M Hugenholtz P . Emerging pathogenic links between microbiota and the gut-lung axis. Nat Rev Microbiol (2017) 15(1):55–63. doi: 10.1038/nrmicro.2016.142 Salem I Ramser A Isham N Ghannoum MA . The gut microbiome as a major regulator of the gut-skin axis. Front Microbiol (2018) 9:1459. doi: 10.3389/fmicb.2018.01459 Chatenoud L Bertuccio P Turati F Galeone C Naldi L Chatenoud L . Markers of microbial exposure lower the incidence of atopic dermatitis. Allergy (2020) 75(1):104–15. doi: 10.1111/all.13990 Lee E Lee SY Kang MJ Kim K Won S Kim BJ . Clostridia in the gut and onset of atopic dermatitis via eosinophilic inflammation. Ann Allergy Asthma Immunol (2016) 117(1):91–92.e1. doi: 10.1016/j.anai.2016.04.019 Galazzo G Van Best N Bervoets L Dapaah IO Savelkoul PH Hornef MW . Development of the microbiota and associations with birth mode, diet, and atopic disorders in a longitudinal analysis of stool samples, collected from infancy through early childhood. Gastroenterology (2020) 158(6):1584–96. doi: 10.1053/j.gastro.2020.01.024 Watts AM West NP Zhang P Smith PK Cripps A W Cox AJ . The gut microbiome of adults with allergic rhinitis is characterised by reduced diversity and an altered abundance of key microbial taxa compared to controls. Int Arch Allergy Immunol (2021) 182(2):94–105. doi: 10.1159/000510536 Yatsunenko T Rey FE Manary MJ Trehan I Dominguez-Bello MG Contreras M . Human gut microbiome viewed across age and geography. Nature (2012) 486(7402):222–7. doi: 10.1038/nature11053 Kosiewicz MM Dryden GW Chhabra A Alard P . Relationship between gut microbiota and development of T cell associated disease. FEBS Lett (2014) 588(22):4195–206. doi: 10.1016/j.febslet.2014.03.019 Smith G D Ebrahim S . 'Mendelian randomization': Can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol (2003) 32(1):1–22. doi: 10.1093/ije/dyg070 Li M Pezzolesi MG . Advances in understanding the genetic basis of diabetic kidney disease. Acta Diabetol (2018) 55(11):1093–104. doi: 10.1007/s00592-018-1193-0 Emdin CA Khera A V Kathiresan S . Mendelian randomization. Jama (2017) 318(19):1925–26. doi: 10.1001/jama.2017.17219 Li Y Li Q Cao Z Wu J . The causal association of polyunsaturated fatty acids with allergic disease: A two-sample mendelian randomization study. Front Nutr (2022) 9:962787. doi: 10.3389/fnut.2022.962787 Van Der Velde KJ Imhann F Charbon B Pang C Van Enckevort D Slofstra M . MOLGENIS research: Advanced bioinformatics data software for non-bioinformaticians. Bioinformatics (2019) 35(6):1076–78. doi: 10.1093/bioinformatics/bty742 Sanna S van Zuydam NR Mahajan A Kurilshikov A Vich Vila A Võsa U . Causal relationships among the gut microbiome, short-chain fatty acids and metabolic diseases. Nat Genet (2019) 51(4):600–05. doi: 10.1038/s41588-019-0350-x Zhang Y Zhang X Chen D Lu J Gong Q Fang J . Causal associations between gut microbiome and cardiovascular disease: A mendelian randomization study. Front Cardiovasc Med (2022) 9:971376. doi: 10.3389/fcvm.2022.971376 Kamat MA Blackshaw JA Young R Surendran P Burgess S Danesh J . PhenoScanner V2: An expanded tool for searching human genotype-phenotype associations. Bioinformatics (2019) 35(22):4851–53. doi: 10.1093/bioinformatics/btz469 Lawlor DA Harbord RM Sterne JA Timpson N Davey Smith G . Mendelian randomization: Using genes as instruments for making causal inferences in epidemiology. Stat Med (2008) 27(8):1133–63. doi: 10.1002/sim.3034 Burgess S Thompson SG . Avoiding bias from weak instruments in mendelian randomization studies. Int J Epidemiol (2011) 40(3):755–64. doi: 10.1093/ije/dyr036 Burgess S Butterworth A Thompson SG . Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol (2013) 37(7):658–65. doi: 10.1002/gepi.21758 Greco MF Minelli C Sheehan N A Thompson JR . Detecting pleiotropy in mendelian randomisation studies with summary data and a continuous outcome. Stat Med (2015) 34(21):2926–40. doi: 10.1002/sim.6522 Bowden J Davey Smith G Burgess S . Mendelian randomization with invalid instruments: effect estimation and bias detection through egger regression. Int J Epidemiol (2015) 44(2):512–25. doi: 10.1093/ije/dyv080 Slob EAW Burgess S . A comparison of robust mendelian randomization methods using summary data. Genet Epidemiol (2020) 44(4):313–29. doi: 10.1002/gepi.22295 Bowden J Davey Smith G Haycock P C Burgess S . Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol (2016) 40(4):304–14. doi: 10.1002/gepi.21965 Liu B Ye D Yang H Song J Sun X Mao Y . Two-sample mendelian randomization analysis investigates causal associations between gut microbial genera and inflammatory bowel disease, and specificity causal associations in ulcerative colitis or crohn's disease. Front Immunol (2022) 13:921546. doi: 10.3389/fimmu.2022.921546 Burgess S Thompson SG . Interpreting findings from mendelian randomization using the MR-egger method. Eur J Epidemiol (2017) 32(5):377–89. doi: 10.1007/s10654-017-0255-x Verbanck M Chen CY Neale B Do R . Detection of widespread horizontal pleiotropy in causal relationships inferred from mendelian randomization between complex traits and diseases. Nat Genet (2018) 50(5):693–98. doi: 10.1038/s41588-018-0099-7 Xiang K Wang P Xu Z Hu YQ He YS Chen Y . Causal effects of gut microbiome on systemic lupus erythematosus: A two-sample mendelian randomization study. Front Immunol (2021) 12:667097. doi: 10.3389/fimmu.2021.667097 Han CY Kwon SK Yeom M Hahm DH Park JW Park HJ . Exploring the differences in the gut microbiome in atopic dermatitis according to the presence of gastrointestinal symptoms. J Clin Med (2022) 11(13):3690. doi: 10.3390/jcm11133690 Larsen JM Musavian HS Butt TM Ingvorsen C Thysen A H Brix S . Chronic obstructive pulmonary disease and asthma-associated proteobacteria, but not commensal prevotella spp., promote toll-like receptor 2-independent lung inflammation and pathology. Immunology (2015) 144(2):333–42. doi: 10.1111/imm.12376 Laigaard A Krych L Zachariassen LF Ellegaard-Jensen L Nielsen DS Hansen AK . Dietary prebiotics promote intestinal prevotella in association with a low-responding phenotype in a murine oxazolone-induced model of atopic dermatitis. Sci Rep (2020) 10(1):21204. doi: 10.1038/s41598-020-78404-0 Kwon MS Lim SK Jang JY Lee J Park HK Kim N . Lactobacillus sakei WIKIM30 ameliorates atopic dermatitis-like skin lesions by inducing regulatory T cells and altering gut microbiota structure in mice. Front Immunol (2018) 9:1905. doi: 10.3389/fimmu.2018.01905 Looman KIM van Meel ER Grosserichter-Wagener C Vissers FJM Klingenberg JH de Jong NW . Associations of Th2, Th17, treg cells, and IgA (+) memory b cells with atopic disease in children: The generation r study. Allergy (2020) 75(1):178–87. doi: 10.1111/all.14010 Morotomi M Nagai F Sakon H Tanaka R . Dialister succinatiphilus sp. nov. and barnesiella intestinihominis sp. nov., isolated from human faeces. Int J Syst Evol Microbiol (2008) 58(Pt 12):2716–20. doi: 10.1099/ijs.0.2008/000810-0 Kumari M Dasriya VL Nataraj BH Nagpal R Behare PV . Lacticaseibacillus rhamnosus-derived exopolysaccharide attenuates d-Galactose-Induced oxidative stress and inflammatory brain injury and modulates gut microbiota in a mouse model. Microorganisms (2022) 10(10):2046. doi: 10.3390/microorganisms10102046 Kim HJ Kim HY Lee SY Seo JH Lee E Hong SJ . Clinical efficacy and mechanism of probiotics in allergic diseases. Korean J Pediatr (2013) 56(9):369–76. doi: 10.3345/kjp.2013.56.9.369 Song H Yoo Y Hwang J Na Y C Kim HS . Faecalibacterium prausnitzii subspecies-level dysbiosis in the human gut microbiome underlying atopic dermatitis. J Allergy Clin Immunol (2016) 137(3):852–60. doi: 10.1016/j.jaci.2015.08.021 Anania C Brindisi G Martinelli I Bonucci E D'Orsi M Ialongo S . Probiotics function in preventing atopic dermatitis in children. Int J Mol Sci (2022) 23(10):5409. doi: 10.3390/ijms23105409 Arpaia N Campbell C Fan X Dikiy S van der Veeken J deRoos P . Metabolites produced by commensal bacteria promote peripheral regulatory T-cell generation. Nature (2013) 504(7480):451–5. doi: 10.1038/nature12726 Luu M Monning H Visekruna A . Exploring the molecular mechanisms underlying the protective effects of microbial SCFAs on intestinal tolerance and food allergy. Front Immunol (2020) 11:1225. doi: 10.3389/fimmu.2020.01225 Du Y Li X An Y Song Y Lu Y . Association of gut microbiota with sort-chain fatty acids and inflammatory cytokines in diabetic patients with cognitive impairment: A cross-sectional, non-controlled study. Front Nutr (2022) 9:930626. doi: 10.3389/fnut.2022.930626 Yip W Hughes MR Li Y Cait A Hirst M Mohn WW . Butyrate shapes immune cell fate and function in allergic asthma. Front Immunol (2021) 12:628453. doi: 10.3389/fimmu.2021.628453 Mahdavinia M Rasmussen HE Botha M Binh Tran TD Van den Berg JP Sodergren E . Effects of diet on the childhood gut microbiome and its implications for atopic dermatitis. J Allergy Clin Immunol (2019) 143(4):1636–37.e5. doi: 10.1016/j.jaci.2018.11.034 Curtis MM Hu Z Klimko C Narayanan S Deberardinis R Sperandio V . The gut commensal bacteroides thetaiotaomicron exacerbates enteric infection through modification of the metabolic landscape. Cell Host Microbe (2014) 16(6):759–69. doi: 10.1016/j.chom.2014.11.005 Graessler J Qin Y Zhong H Zhang J Licinio J Wong ML . Metagenomic sequencing of the human gut microbiome before and after bariatric surgery in obese patients with type 2 diabetes: correlation with inflammatory and metabolic parameters. Pharmacogenomics J (2013) 13(6):514–22. doi: 10.1038/tpj.2012.43 Arumugam M Raes J Pelletier E Le Paslier D Yamada T Mende DR . Enterotypes of the human gut microbiome. Nature (2011) 473(7346):174–80. doi: 10.1038/nature09944 Nosrati A Afifi L Danesh MJ Lee K Yan D Beroukhim K . Dietary modifications in atopic dermatitis: patient-reported outcomes. J Dermatolog Treat (2017) 28(6):523–38. doi: 10.1080/09546634.2016.1278071 Corrêa RO Castro PR Moser R Ferreira CM Quesniaux VFJ Vinolo MAR . Butyrate: Connecting the gut-lung axis to the management of pulmonary disorders. Front Nutr (2022) 9:1011732. doi: 10.3389/fnut.2022.1011732 Gradisteanu Pircalabioru G Liaw J Gundogdu O Corcionivoschi N Ilie I Oprea L . Effects of the lipid profile, type 2 diabetes and medication on the metabolic syndrome-associated gut microbiome. Int J Mol Sci (2022) 23(14):7509. doi: 10.3390/ijms23147509 Qin P Zou Y Dai Y Luo G Zhang X Xiao L . Characterization a novel butyric acid-producing bacterium collinsella aerofaciens subsp. shenzhenensis subsp. Nov. Microorganisms (2019) 7(3):78. doi: 10.3390/microorganisms7030078 Huang X Chen Y Zhang F Yang Q Zhang G . Peripheral Th17/Treg cell-mediated immunity imbalance in allergic rhinitis patients. Braz J Otorhinolaryngol (2014) 80(2):152–5. doi: 10.5935/1808-8694.20140031 Xu LZ Yang LT Qiu SQ Yang G Luo XQ Miao BP . Combination of specific allergen and probiotics induces specific regulatory b cells and enhances specific immunotherapy effect on allergic rhinitis. Oncotarget (2016) 7(34):54360–69. doi: 10.18632/oncotarget.10946 Theiler A Bärnthaler T Platzer W Richtig G Peinhaupt M Rittchen S . Butyrate ameliorates allergic airway inflammation by limiting eosinophil trafficking and survival. J Allergy Clin Immunol (2019) 144(3):764–76. doi: 10.1016/j.jaci.2019.05.002 Thio CL Chi PY Lai A C Chang YJ . Regulation of type 2 innate lymphoid cell-dependent airway hyperreactivity by butyrate. J Allergy Clin Immunol (2018) 142(6):1867–83.e12. doi: 10.1016/j.jaci.2018.02.032