Anti-Inflammatory and Anti-asthmatic Effects of TMDCT Decoction in Eosinophilic Asthma Through Treg/Th17 Balance

Tuo-Min-Ding-Chuan decoction (TMDCT) is a Traditional Chinese Medicine (TCM) formula consisting of twelve herbs that can relieve the symptoms and treat allergic asthma. Yet, the underlying mechanism of action is still unclear. In this study, we investigated the effect of TMDCT in regulating Treg/Th17 cells immune balance and explored potential metabolic and gut biomarkers associated with Treg and Th17 cells in eosinophilic asthma mice treated by TMDCT. We found that TMDCT increases Treg cells percentage and decreases Th17 cells percentage in the ovalbumin (OVA) -induced eosinophilic asthma mice model. Furthermore, Imidazoleacetic acid, dL-glutamine, L-pyroglutamic acid, 2-deoxy-d-glucose were preliminary identified as biomarkers in plasma metabolites treated by TMDCT, meanwhile genus Desulfovibrio, genus Butyricimonas and genus Prevotella 9 were preliminary identified as gut microbiota biomarkers after TMDCT treatment. These results provide an experimental foundation for the treatment of allergic asthma with Chinese herbal compounds.


INTRODUCTION
Asthma is a chronic inflammatory disease which is difficult to control affects more than 300 million people worldwide. With current rising trends, it is expected that this number will reach 400 million people by 2025. According to some reports, there are nearly 250,000 asthma-related deaths each year, many of which are avoidable (Christiansen and Zuraw, 2019).
Inflammatory disorder in allergic asthma is characterized by broncho-construction, bronchial hyperresponsiveness, and even tissue damage and the most common type of asthma (National Asthma Education and Prevention Program, 2007). If not properly treated, allergic asthma can progress to chronic obstructive lung disease or other disease associated with airways and lung tissue. The burden of allergic diseases has increased over recent years, as evidenced by a high incidence of patients suffering from these diseases and incurring high financial costs (Bousquet et al., 2016). Over the years, anti-inflammatory drugs have been developed, such as inhaled steroids and bronchodilators that can relieve symptoms. Yet, existing drugs cannot completely cure the patient, and have an elevated recurrence rate. Moreover, conventional drugs have been associated with certain side effects. E.g., an inhaled corticosteroid (ICS), the dominant treatment in type 2-high asthmatic inflammation, can suppress the endocrine system, which assists body's immune system fight against infection, while long-term steroid treatments can lead to anxiety and depression. Allergic asthma is induced by an abnormal type 2 immune response to inhaled allergens, such as house dust mites (HDM), grass pollen, animal dander, and mold (Caminati et al., 2018). In the airway, the inflammatory response is mainly caused by Th2 type inflammation, inducing cytokines IL-4, IL-5, and IL-13 overexpression, which consequently activate the expression of IgE and the infiltration of eosinophils and mast cells in the airway.
Previous studies have shown that the imbalance of Th1/Th2 cells may be involved in the pathogenesis of airway inflammation in asthma (Berker et al., 2017). Besides that, some recent studies suggested that insufficient Th1 cell differentiation is not the only cause leading to over-differentiation and activation of Th2 cells. It has been discovered that Treg cells have an essential role in regulating the body's immune balance. Dysfunctional Treg cells, unable to effectively suppress excessive Th2 response, may lead to asthma and other allergic diseases.
Th2-type inflammation is involved in all types of asthma (mild, moderate, and severe). Also, inflammation caused by Th17 cells has been associated with a progression of asthma (Israel and Reddel, 2017). The imbalance between Treg cells and Th cells such as Th1, Th2 and Th17 cell responses in allergic asthma is associated with the development of allergic asthma (Noval Rivas and Chatila, 2016). Genetic and immunological evidence suggested that Treg cells can promote tolerance to allergens and prevent allergic disorders (Palomares et al., 2014). Thus, it is believed that the inflammation of allergic asthma inflammation and increasing Treg cells number or stimulating the proliferation of immunosuppressive cells such as Treg cells is a new strategy for the treatment of allergic asthma.
Recently, some new drugs with the ability to exert Treg cells in the body's immune system or stimulate differentiation of Treg cells have been tested. Although great development has been achieved using allergen immunotherapy (AIT) and other immunotherapy, treatment of allergic asthma and its recurrence remains challenging. New treatment options that can regulate the immune balance of Treg/Th17 are currently emerging.
TMDCT is a TCM formula consisting of twelve herbs that can relieve the symptoms and treat allergic asthma. Nevertheless, the underlying mechanism of action is still unclear. TMDCT contain 12 kinds of traditional Chinese medicine, it can suppresses inflammation of allergic asthma and inhibit the degranulation of mast cell (Qin et al., 2021). In our study, we investigated the effect of TMDCT in regulating Treg/Th17 cells balance and explored potential metabolic and gut biomarkers associated with Treg and Th17 cells in eosinophilic asthma mice treated with TMDCT using nontargeted metabolome and 16S rDNA technology.

Mice
BALB/c mice (female, 6-8 weeks old, 17-20 g) were obtained from Beijing Vital River Laboratory Animal Technology Co., Ltd in China. Keeping mice in specific pathogen-free conditions in Beijing University of Chinese Medicine. All mice were kept at a controlled room (25 ± 1°C, 45-60% humidity). All animal studies were conducted in accordance with the institutional animal care regulations of Beijing University of Chinese Medicine and were conducted in accordance with AAALAC and IACUC guidelines.

OVA-Induced Allergic Asthma Mice Model and Grouping
The OVA-induced eosinophilic asthma BALB/c mice model was constructed following a previously described approach (Dumas et al., 2018). Briefly, mice were intraperitoneally injected with 2 mg of OVA (Sigma-Aldrich, Cat#A5503) mixed with 2 mg Imject ™ Alum Adjuvant (Invitrogen, Cat#77161) and PBS on day 0 and day 14. From the 21st to 25th days post-injection, in the challenge phase, mice are continuously nebulized with 1% OVA for 30 min as shown in Figure 1A.
OVA-induced allergic asthma mice were then divided in the following groups: High-dose TMDCT group (40.56 g/kg/d TMDCT (H group)), Middle-dose TMDCT treatment group (20.28 g/kg/d TMDCT (M group)), Low-dose TMDCT treatment group (10.14 g/kg/d TMDCT (L group)), and positive group (dexamethasone 1 mg/kg/d). TMDCT, containing 12 kinds of Chinese medicine, was prepared as previously described (Qin et al., 2021). TMDCT and dexamethasone (1 mg/kg/d, a positive drug for suppressing eosinophilic asthma) were given daily at the same time from the 21st to 25th days 1 hour before challenge.

Histological Sections and Pathology Scoring
Lungs were fixed with 4% paraformaldehyde (PFA) at 4°C, after paraffin-embedded sections and stained them with hematoxylin and eosin (HE) (to examine cell infiltration detection) or periodic acid-schiff stain (PAS) (to examine mucus production) in lung tissues. Scoring of inflammatory cells and goblet cells was performed in at least three different fields for each lung section as described (Duan et al., 2004). The mean scores were calculated using five animals.

Bronchoalveolar Lavage Fluid, Serum Cytokines, and the Culture Supernatant of Treg and Th17 Cells Cytokines Analysis
Twenty-four hours after the last aerosol challenge, bronchoalveolar lavage (BAL) fluid was collected by syringe three times with 1 ml PBS containing 1% BSA. Eosinophils in BAL fluid were counted by cell sorting and counting instrument. Cytokine levels of IL-4, IL-5, IL-13, IL-10, IL-17A, IL-6 in BALF were analyzed using a premixed AimPlex ™ multiplex-assay kit (Cat#T2C0710709), TGF-β was analyzed using a premixed AimPlex ™ multiplex-assay kit (Cat#B111206); OVA-specific IgE in serum was detected by ELISA kit (Cayman, Cat#500840).

Untargeted Metabolomics Detection of Plasma
Plasma was collected with 1.5 ml Eppendorf Tubes containing EDTA (ethylene diamine tetraacetic acid), which were centrifuged at 4°C by 1,500 g, 15 min. Rremove the protein by Methanol/acetonitrile (1:1, v/v) and centrifuged at 14000g, 4°C for 15 min. All LC-MS analyses were performed at Shanghai Applied Protein Technology Co., Ltd.
R package (ropls) was used to analyzed the processed data. Unsupervised principal component analysis (PCA) and supervised orthogonal partial least squares discriminant analysis (OPLS-DA) were used to evaluate sample stability. The VIP (variables in the projection) value of the OPLS-DA model was calculated, metabolites with VIP ˃ 1, p-value ˂ 0.05 is considered to be the significant changed metabolites.
After multivariate and univariable analysis, searched the significant metabolites in the Human Metabolome database (http://www.hmdb.ca), METLIN (https://metlin. scripps.edu), KEGG (http://www.kegg.com) and Chemspider (http://www. chemspider.com/). Cytoscape software was used to perform enrichment analysis and visualization. Pearson correlation analysis was used to perform and determine the correlation between two variables. KEGG enrichment analysis was performed using MetaboAnalyst (www.metaboanalyst.ca). The volcano plot and clustering analysis were performed using R.

Gut Microbiota Analysis
Using CTAB/SDS method to extracte the total genome DNA. Then, the V3-V4 regions in 16S rDNA were amplified using a specific primer included in the barcode. Next, an Illumina Miseq/ HiSeq2500 platform was used to build a library. We used PCoA (Principal Co-ordinates Analysis) to study the similarity or difference of sample community composition, and LEfSe (LDA Effect Size) to quantify the biomarkers in different groups.

Statistical Analysis
In this study, t-test, one-way ANOVA test (Turkey or Dunnett), Wilcoxon rank-sum test, Kruskal-Wallis test or Wilcoxon ranksum test was used. p value ˂ 0.05 was considered as statistically significant. Pearson analysis with R3.4.2 Heatmap was used to analyze the correlation between metabolites and other index, the correlation between gut microbiota and other index.

TMDCT Alleviates Airway Inflammation in BALB/C Mice With Eosinophilic Asthma
To evaluate the treatment effect of TMDCT in treating eosinophilic asthma, an OVA-induced eosinophilic asthma model was constructed as described in Figure 1A. The most significant effect of TMDCT(H) group in reducing cell infiltration ( Figure 1C) and decreasing goblet cell hyperplasia was observed in the TMDCT(H) group ( Figure 1D), further indicating that TMDCT can alleviate airway inflammation. Moreover, TMDCT reduced airway resistance (Penh ratio), which indicated that lung function in TMDCT with a high concentration group had the best effect ( Figure 1B).
To investigate the role of TMDCT in allergic airway inflammation, we detected OVA-specific IgE levels in serum. TMDCT significantly decreased OVA-specific IgE concentration in serum ( Figure 1E). Additionally, the eosinophils number in BALF was reduced, and IL-5 levels were decreased in BALF of the TMDCT group, especially in the H group ( Figure 1E). Also, TMDCT significantly reduced IL-4 and IL-13 cytokines in BALF. These data suggested that TMDCT could alleviate allergic airway inflammation.

TMDCT can Increase Treg Cells Percentage and Decrease Th17 Cells Percentage in the Spleen
Th17 and Treg cells have important role in driving and restraining airway inflammation in patients with asthma (Seumois et al., 2020). To further investigate whether TMDCT can regulate immune response, we detected Treg and Th17 cells percentage in spleen tissues. Surprisingly, TMDCT increased Treg cells percentage and decreased Th17 cells percentage; the most significant effect was seen in the TMDCT(H) group ( Figures 2B,C).
Next, we analyzed TGFβ, IL-10, IL-6, and IL17A cytokines in BALF. As shown in Figure 2A, TMDCT stimulated TGFβ and IL-10 and suppressed IL-6 and IL17A (p value were all <0.05), which is consistent with the Treg and Th17 percentage detection.
At the mRNA level, TMDCT increased the Foxp3 (transcription factor of Treg cells) and decreased RORγt (transcription factor of Th17 cells) relative expression ( Figure 2D). These data suggested that TMDCT could decrease Th17 cells percentage, and increase Treg cells percentage in spleen tissues, among them TMDCT(H) has the best effect.

Identification Plasma Metabolites of TMDCT
To further explore the mechanism of TMDCT in regulating the immune balance in eosinophilic asthma, we examined potential plasma metabolite biomarkers of TMDCT. The TMDCT with a Frontiers in Pharmacology | www.frontiersin.org February 2022 | Volume 13 | Article 819728 high concentration was used, and the untargeted metabolomics technology was used in this experiment. As seen in Figure 3A, OPLA-DA analysis found the following: 0.3 < Q2<0.384 (in POS mode) < 0.5, 0.3 < Q2=0.329 (in NEG mode) < 0.5, suggesting that the detected model was stable. The upregulated and downregulated metabolites in positive and negative mode can be shown in Figure 3B, Hierarchical cluster of significant difference metabolites (VIP >1, p-value < 0.05) analysis ( Figure 3C) further suggested that metabolites in TMDCT or model group had a similar function; they participated in similar metabolite pathway or cell pathway. Finally, six differential metabolites in positive mode and 16 differential metabolites in the negative mode were detected (Table 1).
Next, we performed a KEGG enrichment analysis. As shown in Figure 3D, several pathways were identified. The significant metabolism pathway were our interest.

Metabolic Biomarkers in Plasma Associated With Treg and Th17 Cells in Eosinophilic Asthma Mice Treated With TMDCT
To understand the mechanism of TMDCT in regulating Treg/ Th17 cells balance, the different metabolites between TMDCT and model group were used for correlation analysis, including Penh value (Mch:12.5, 25, 50 mg/ml), OVA-IgE, Eosinophils, Treg and Th17 cells percentage. In our study, L-Arabinono-1,4lactone, Myo-inositol, DL-lactate, Uracil and 1-stearoyl-2docosahexaenoyl-sn-glycero-3-phosphocholine were associated with Penh value (Mch: 12.5, 25, 50 mg/ml), OVA-IgE, Eosinophils, Treg and Th17 cells percentage ( Figure 4A). There were obviously different (p < 0.05) in 1-stearoyl-2docosahexaenoyl-sn-glycero-3-phosphocholine and Ethylenediaminetetraacetic acid either between the model and control group, or the model group and the TMDCT group ( Figure 4B), they were all associated with OVA-IgE, Penh value (Mch: 12.5 mg/ml, 25 mg/ml, 50 mg/ml) and Eosinophil numbers while they has no correaltion with Treg or Th17 cells percentage ( Figure 4A). Moreover, Imidazoleacetic acid, DL-glutamine, L-pyroglutamic acid were all increased in model group compared with control group, while they were all decreased in TMDCT treatment group (p value were all <0.05). Besides that, Imidazoleacetic acid, DL-glutamine, L-pyroglutamic acid were all associated with OVA-IgE, Penh value (Mch: 12.5 mg/ml, 25 mg/ ml, 50 mg/ml), Eosinophil numbers and Th17 cells percentage (p value all <0.05). As shown in Figure 4B, 2-deoxy-D-glucose was decreased in model group when compared with control group (p < 0.05), while it increased in TMDCT treatment group (p < 0.05), it is also associated with OVA-IgE, Penh value (Mch: 12.5 mg/ml, 25 mg/ ml, 50 mg/ml), Eosinophil number and Th17 cells percentage (p value all <0.05). The above four metabolites were not correlated with Treg cells, which indicated that they could also affect the Treg/Th17 balance in eosinophilic asthma. The changed metabolites were enriched to phenylalanine metabolism, pyrimidine metabolism, ascorbate and aldarate metabolism pathway, histidine metabolism, glutathione metabolism and metabolic pathways ( Figure 3D).

Variation of Gut Microbiota and Key Phylotypes of Gut Microbiota in Response to OVA-Induced Eosinophilic Asthma Mice Treated by TMDCT
Gut bacteria have an essential role in the action of drugs. Thus, we detected gut microbiota in the Control group and TMDCT (H) group. As shown in Figures 5A,B, the differences in gut microbiota were found between the control and the model group. LEfSe (LDA Effect Size) software was used to discover high-dimensional biomarkers and reveal genome characteristics. LDA ˃ two was used as selection criteria in the study. Through STAMP analysis, we found that genus Rikenellaceae_ RC9_gut_group (p = 0.0497), genus Bifidobacterium (p = 0.0462), genus Rikenella (p = 0.0496), genus mouse gut metagenome (p = 0.0196), genus Butyricimonas (p = 0.0323), genus Prevotella (p = 0.0339), genus Enterococcus (p = 0.0029), We used Pearson correlation analysis with R 3.4.2 Heatmap software to analyze the correlation. As shown in Figure 6, genus Desulfovibrio, genus Muribaculum and genus Prevotella 9 were all associated with Penh value (Mch:12.5, 25, 50 mg/ml), OVA-IgE, Eosinophils number in BALF and Th17 cells percentage. These results indicated a correlation with Treg/Th17 balance in eosinophilic asthma treated by TMDCT.

DISCUSSION
TMDCT is a prescription of professor Qi Wang used to treat allergic asthma which has significant clinical effect. Yet, its mechanism of action still remains unclear. We further investigated the role of TMDCT in regulating Treg/Th17 cells balance in eosinophilic asthma.
The generation and maintenance of Treg cells and their suppressive cytokines (TGF-β, IL-10) are essential for the induction of allergen tolerance in allergic disease (Palomares et al., 2017). In addition, FOXP3+ Treg cells play a key role in maintaining peripheral tolerance by inhibiting the responsiveness to allergens (Noval Rivas and Chatila, 2016). It is also suggested that allergenspecific immunotherapy regimens increase the numbers of Treg cells (Radulovic et al., 2008). Moreover, some studies found that Treg cells have an irreplaceable role in preventing airway inflammation of allergic asthma (Tortola et al., 2019;Sun et al., 2020). While increasing immunosuppressive function of immunosuppressive cells such as Treg or Breg cells may be essential in a complete cure for allergic asthma.
IL-17A is upregulated in allergic asthma patients, and the upregulation of IL-17A is correlated with the severity of allergic  asthma. TGF-β induces the synthesis of the transcription factor of RORγt, which is specific for the Th17 cells (Ivanov et al., 2006). In this study, IL17A and TGF-β, transcription factor RORγt of Th17 cells, and transcription factor FOXP3 of Treg cells were used to analyze the balance of Treg/Th17 cells after TMDCT treatment. Our data suggested that TMDCT can suppress the symptoms of eosinophilic asthma and regulate the balance of Treg/Th17 cells. Interestingly, dexamethasone has similar anti-inflammatory and antiasthma effects when compared with TMDCT in our study. TMDCT can also treat eosinophilic asthma by regulating Treg/Th17 immune balance the same as dexamethasone does (Kianmehr et al., 2017). While there are too many drawbacks of dexamethasone, a new drug should be explored.
Next, we used non-targeted metabolome to explore potential metabolic biomarkers in the plasma that were associated with Treg and Th17 cells in eosinophilic asthma mice treated with TMDCT. Imidazoleacetic acid, DL-glutamine, 2-deoxy-D-glucose, L-pyroglutamic acid were identified as plasma biomarkers of TMDCT and were all associated with Treg/Th17 cells balance. They were enriched in phenylalanine metabolism, pyrimidine metabolism, ascorbate and aldarate metabolism pathway, histidine metabolism, glutathione metabolism and metabolic pathways, etc. All of the plasma metabolites preliminary proved that TMDCT could regulate Treg/Th17 balance in eosinophilic asthma, which is also the advantage of multi-target effects of traditional Chinese medicine TMDCT. Previous studies found that glutamine can distinct regulate Th17 differentiation (Johnson et al., 2018), which can better explain the reasons for the change in Th17 cells percentage which is consistent with our results. While, until now there are no study found imidazoleacetic acid, 2-deoxy-D-glucose, L-pyroglutamic acid neithor can regulate Treg, Th17 cells nor regulate the differentiation of Treg or Th17 cells, which are not found associated with Treg/Th17 balance in allergic asthma.
Gut, lung, and skin microbiome exposures can also influence the occurrence and the development of allergy disease (Kemter and Nagler, 2019). The gut microbiota is essential in systemic immune regulatory network of allergic disease (Sommer and Bäckhed, 2013) and can regulate drug pharmacokinetics, such as absorption and distribution (Sousa et al., 2008;Tremaroli and Bäckhed, 2012;El Aidy et al., 2015;Zhang et al., 2018). Besides that, the gut-lung axis transfers metabolites and immunomodulatory signals between the gut and lungs. Many studies have shown an increase in the number of respiratory diseases due to deviations in gut ecology (Reddy et al., 2012;Budden et al., 2017). Omics technologies, including metabolomics, can be used to begin to understand relevant molecular changes. In this study, we used omics technology to detect the effect on of TMDCT in regulating Treg/Th17 cells, and then used Non-targeted metabolome and 16S rRNA technology to analyze the biomarkers after TMDCT treatment and relative factors correlated with Treg/Th17 balance. We found that gut microbiota genus Rikenellaceae_ RC9_gut_group, genus Bifidobacterium, genus Rikenella, genus mouse gut metagenome, genus Butyricimonas, genus Prevotella, genus Enterococcus, genus Peptoniphilus, genus Dialister, genus Corynebacterium, genus Dermabacter, genus Varibaculum were significantly different after TMDCT treatment.
Bifidobacterium is anti-inflammatory bacteria that can relieve allergic asthma in mice by regulating Th1/Th2 balance (Wang et al., 2020). Allergic patients with chronic asthma usually have low levels of Bifidobacterium (Hevia et al., 2016).
Oral administration of Enterococcus faecalis can suppress allergic asthmatic response associated with attenuation of Th17 cell development (Zhang et al., 2012). Fecal transplantation containing gut Rikenellaceae bacteria can alleviate acute liver injury in mice through regulating Treg/ Th17 balance (Liu et al., 2021), while it cannot alleviate eosinophilic asthma. Corynebacterium may regulate Th17/Treg cells in the intestinal mucosal immunity (Zhang et al., 2019); However, there are still no studies in allergic asthma. So far, we have not find relevant research about genus Desulfovibrio, genus Muribaculum or genus Prevotella 9 are related with Th17 cells percentage in eosinophilic asthma disease, they were the gut microbiota biomarker after TMDCT treatment.
This study has some limitations. First, only the BALB/c eosinophilic asthma model was used in this study. The effect of TMDCT should be further examined in other asthma models. Also, clinical research and other allergic asthma models are required to prove the mechanism of TMDCT in treating allergic asthma. Much more study of the biomarkers of TMDCT treatment should be developed to further explain the mechanism of TMDCT in treating allergic asthma, which will be a new strategy to treat allergic asthma.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm. nih.gov/, PRJNA799688.

ETHICS STATEMENT
The animal study was reviewed and approved by Beijing University of Chinese Medicine.

AUTHOR CONTRIBUTIONS
YZ, JW, XZ and QW conceived and designed the experiments; YZ, HZ and TW performed the experiments; YZ analyzed the data and wrote the paper. All authors have read and agreed to the published version of the manuscript.