- 1Affiliated Hospital of Nanjing University of Chinese Medicine, Nanjing, China
- 2Pediatric Institution of Nanjing University of Chinese Medicine, Nanjing, China
- 3Jiangsu Key Laboratory of Pediatric Respiratory Disease, Nanjing, China
Airway remodeling is a primary pathological feature of asthma. The current therapy for asthma mainly targets reducing inflammation but not particularly airway remodeling. Therefore, it is worthwhile to develop alternative and more effective therapies to attenuate remodeling. Gu-Ben-Fang-Xiao Decoction (GBFXD) has been used to effectively and safely treat asthma for decades. In this study, GBFXD regulated airway inflammation, collagen deposition, and the molecules relevant to airway remodeling such as Vimentin, α-SMA, hydroxyproline, and E-cadherin in chronic remission asthma (CRA) murine model. Proteomic analysis indicated that the overlapping differentially expressed proteins (DEPs) (Model/Control and GBFXD/Model) were mainly collagens and laminins, which were extracellular matrix (ECM) proteins. In addition, the KEGG analysis showed that GBFXD could regulate pathways related to airway remodeling including ECM-receptor interactions, focal adhesion, and the PI3K/AKT signaling pathway, which were the top three significantly enriched pathways containing the most DEPs for both Model/Control and GBFXD/Model. Further validation research showed that GBFXD regulated reticulon-4 (RTN4) and suppressed the activation of the PI3K/AKT pathway to alleviate ECM proteins deposition. In conclusion, our findings indicate that GBFXD possibly regulate the PI3K/AKT pathway via RTN4 to improve airway remodeling, which provides a new insight into the molecular mechanism of GBFXD for the treatment of CRA.
Introduction
Asthma is a heterogeneous disease of high prevalence worldwide, characterized by respiratory symptoms such as recurrent wheeze, coughing, and shortness of breath caused by lumen narrowing and airflow obstruction (Camoretti-Mercado and Lockey, 2020). Epidemiological studies predict that the number of patients with asthma will reach 400 million by 2025 (Voskamp et al., 2020). Asthma is the most prevalent chronic respiratory disease in children, affecting approximately 10% of the general population and up to 20% of children, with approximately a 50% increase in prevalence each decade (Croisant, 2014). Inhaled corticosteroids (ICS) has been the most effective anti-inflammatory medication in controlling symptoms and improving the quality of life for patients with asthma (Costello and Cushen, 2020). However, the effect of ICS on attenuating airway remodeling remains controversial (Halwani et al., 2016; Yap et al., 2019) and the current therapies for asthma do not particularly target airway remodeling (Prakash et al., 2017). The underlying mechanism of airway remodeling in asthma remains unclear, but chronic inflammation and/or mechanical stresses from repeated bronchoconstriction independent of inflammation play vital roles (Grainge et al., 2011; Prakash et al., 2017). Importantly, airway remodeling may occur early during asthma and is also related to airway hyperresponsiveness (AHR) (Jendzjowsky and Kelly, 2019). Therefore, developing alternative and more effective therapies to attenuate airway inflammation and remodeling is essential and urgent.
Traditional Chinese medicine (TCM) formulae are usually composed of one or two medicinal herbs as the principal component and the others act as adjuvants to enhance the effects (Tao et al., 2017). The therapeutic efficacy of TCM formulae is widely accepted to treat disease from a holistic point of view. High-throughput proteomic analysis of the functional proteins can help clarify the complex mechanisms of TCM. Gu-Ben-Fang-Xiao Decoction (GBFXD), an empirical formula prescribed by a famous Chinese professor Jiang Yuren, contains 11 herbs (Table 1) and has been used to treat pediatric asthma for decades. Previous clinical studies in which the decoction obtained by the above mentioned herbs were administrated for children orally once a day have confirmed the effectiveness and safety in reducing asthma recurrence (Yuan et al., 2010). Our previous research demonstrated that GBFXD could modulate the microbiota-acetate-regulatory T cell (Treg) axis (Dong et al., 2020), downregulate the asthma susceptibility gene ORMDL3 (Huang et al., 2016), and suppress endoplasmic reticulum (ER) stress responses (Lu et al., 2016) in a murine model of chronic remission asthma (CRA). In addition, GBFXD has been shown to regulate cholesterol transport, activate complement factors (Xing et al., 2019), and inhibit macrophage polarization (Liu et al., 2019) in a murine model of chronic persistent asthma (CPA); however, little is known regarding the mechanism of GBFXD on airway remodeling in CRA model.
The structural changes of the airway wall known as airway remodeling, which may explain persistent airflow obstructions and recurrent attacks in some asthmatic patients, can occur early in childhood and increase the risk of advancing into clinical asthma (Bonato et al., 2018). Recent studies reported that airway remodeling in asthmatic children could also arise independently on inflammation (Baraldo et al., 2011), thus much attention should be paid to the prevention and treatment of airway remodeling. Increased deposition of ECM proteins in the reticular basement membrane region, lamina propria, and submucosa is a feature of asthmatic airways and facilitates the airway wall thickening (Hough et al., 2020). Meantime, the molecules such as collagen, fibronectin related to ECM are used to confirm the airway modeling. Collagen is the most abundant component of ECM in the lung and a potential marker of remodeling in asthma (Ojiaku et al., 2017). There is an increase in the deposition of collagen I and collagen III in asthma by bronchial biopsies (Benayoun et al., 2003; Chakir et al., 2003). In addition, fibronectin is the ECM adhesion protein and increased fibronectin level of the airway wall smooth muscle in fatal asthma was reported (Araujo et al., 2008). Emerging studies reported that basement membrane thickening is still present even in clinical and complete remissions (van den Toorn et al., 2001; Broekema et al., 2011). Our previous study reported GBFXD reduced airway remodeling by Masson staining in a murine model of CRA (Huang et al., 2016), but the mechanism remains unknown. Carpaij et al. declared consistent findings that inflammatory markers were elevated in asthma remission compared to those in the healthy controls, and were lower compared to those in persistent asthma cases (Carpaij et al., 2019). Similarly, in previous study we established the murine models of CPA and CRA, airway inflammation, epithelial goblet cells, and AHR were more severe in the CPA model than those in the CRA model (Liu et al., 2019). Considering the differences between CPA and CRA as well as multiple TCM targets, we intend to confirm differentially expressed proteins (DEPs) after GBFXD treatment and further elucidate the possible mechanisms of GBFXD in CRA by iTRAQ-based proteomics.
Materials and Methods
GBFXD Preparation
All herbs were purchased from the Affiliated Hospital of Nanjing University of Chinese Medicine, identified by Professor Shengjin Liu from Nanjing University of Chinese Medicine, and deposited into the Herbarium of Traditional Chinese Medicine, Nanjing University of Chinese Medicine (Table 1). The herbs of GBFXD were mixed, decocted as previously described (Dong et al., 2020), and stored with a final concentration of 3 g/ml of crude herbs for further use.
Components of GBFXD Determined by UPLC-ESI/LTQ-Orbitrap-MS
The quality control of herbs was performed based on previous established methods (Dong et al., 2020) with minor modifications. Briefly, a Waters AQUITY H-Class ultra-high-performance liquid chromatograph (UPLC) was used and chromatographic separation was achieved by the AQUITY UPLC HSS T3 chromatographic column (2.1 × 100 mm, 1.8 μm). The column oven temperature was set at 25°C. The 0.1% formic acid in water (A) and 0.1% formic acid in acetonitrile (B) were chosen as the mobile phases at a flow rate of 0.2 ml/min. The elution conditions were as follows: 0 min: 10% B, 3 min: 10% B, 26 min: 90% B, 30 min: 10% B. The LTQ ORBITRAP XL tandem mass spectrometer was used with the following parameters: ion source: HESI, scan ranges: 100–1,000 m/z, and collision energy: 35. The positive ionization mode was used with the following parameters: heater temperature at 300°C, sheath gas flow rate at 45 arb, aux gas flow rate at 10 arb, spray voltage at 4 kV, capillary temperature at 300°C, capillary voltage at 35 V, and tube lens voltage at 110 V. In addition, the negative ionization mode was used the following parameters: heater temperature at 300°C, sheath gas flow rate at 45 arb, aux gas flow rate at 10 arb, spray voltage at 4 kV, capillary temperature at 300°C, capillary voltage at −35 V, tube lens voltage at −110 V. The reference standards included Prim-o-glucosylcimifugin (catalog number: 141021), Liquiritin (catalog number: 141029), Lobetyolin (catalog number: 140909), Astragaloside IV (catalog number: 140321), Magnolin (catalog number: 140,820), Schisandrol A (catalog number: 141211), Schisandrin A (catalog number: 19070202) (Chengdu Pufei De Biotech Co. Ltd., China), Calycosin-7-O-glucoside (catalog number: 120819), Isoliquiritin (catalog number: 120903), Glycyrrhizic acid (catalog number: 121213), Astragaloside I (catalog number: 120525), Astragaloside II (catalog number: 120528), Astragaloside III (catalog number: 120329), and Calycosin (catalog number: 120615) (Chengdu Herbpurity Co. Ltd., China). The detailed chemical information of GBFXD was shown in Figure 1 and Table 2.
FIGURE 1. UPLC chromatograms analysis of GBFXD with the representative active ingredients. (A) The total chromatograms in positive ion mode (black, GBFXD; Red, standards); (B) The total chromatograms in negative ion mode (black, GBFXD; Red, standards).
Animal Experiment and Sample Preparation
Four-week-old female specific-pathogen-free BALB/c mice (weight: 16–18 g, 8 mice per group) were purchased from Charles River Laboratories (Beijing, China) and the murine model of CRA was established as previously described (Dong et al., 2020). Briefly, the mice were subjected to adaptive feeding for one week prior to treatment. On Days 1 and 8, mice were sensitized with an intraperitoneal injection of 0.2 ml normal saline (NS) containing 100 μg OVA mixed with 1 mg Al(OH)3. The mice were challenged with 2.5% erosolized OVA for 30 min/day from Days 15–28, and once every 3 days from Days 32–55. On Days 29, 42, and 55, the mice were intranasally administered with 50 μL (1.0 × 105 TCID50/mL) respiratory syncytial virus (RSV). From Days 56–85, the mice in the treatment groups were intragastrically administered with different doses of GBFXD (12, 24, 36 g/kg) or Montelukast (2.6 mg/kg) per day. Montelukast, the leukotriene receptor antagonist and clinically most widely used for asthma remission (Lazarinis et al., 2018), was used as a positive control drug for GBFXD. In addition, the dosage conversion from human to mouse was calculated with reference to the book of “Research Methods in Pharmacology of Chinese Materia Medica” edited by Chen Qi (Chen, 2011). The dosage of GBFXD (36 g/kg) for mouse was converted from the clinical daily dosage for human. Due to the further increase in the concentration will enhance the difficulty of intragastric gavage, so the doses (12, 24, and 36 g/kg) are set and also consistent with our previous study (Huang et al., 2016; Lu et al., 2016; Dong et al., 2020). Mice in the control and model group received the same volume of NS. The procedure is shown in Figure 2A. All experimental procedures were performed and approved by the animal ethics committee of Nanjing University of Chinese medicine (No. 201809A019, No. 201810A021).
FIGURE 2. Effect of GBFXD on airway inflammation by hematoxylin and eosin (HE) staining, remodeling by Sirius red staining and HYP levels. (A) Experimental protocol to establish the murine OVA-RSV induced CRA murine model and the treatment strategy. (B) HE staining examined inflammatory infiltration and Sirius red staining identified collagen deposition (×200), Scale bar: 100 μm. (C–D) Total inflammation scores and the quantification of Sirius red staining by AOD in each group mice. (E) The HYP levels in lung of each group mice measured by ELISA. Values are presented as mean ± SD (n = 3 or 4 mice per group). #p < 0.05 compared to control group, and *p < 0.05 compared to model group.
Histopathology and Immunohistochemistry
The middle sections of the right lungs in each group were fixed with 4% paraformaldehyde. The samples were embedded in paraffin, cut into sections of 3–4 μm in thickness, and stained with hematoxylin-eosin (HE) or Sirius red. The inflammation score was determined as previously described (Xing et al., 2019). Other sections from each group were analyzed by immunohistochemistry to explore the expressions of E-cadherin, Vimentin (Servicebio, Wuhan, China), and α-SMA (Proteintech, Wuhan, China). The images were captured by a fluorescence microscope (Olympus BX43, Tokyo, Japan) at ×200 magnification. The results of staining and immunohistochemistry were analyzed by Image-Pro Plus 6.0 and quantified as average optical density (AOD).
Enzyme-Linked Immunosorbent Assay
Hydroxyproline (HYP) levels in the lungs from each group were determined by ELISA (Westang, Shanghai, China), according to manufacturer’s protocol.
Proteomics
Protein Extraction and Digestion
The frozen lung tissues were powdered in liquid nitrogen and proteins were extracted by the RIPA buffer (Beyotime, Haimen, China) with 10 mM dithiothreitol and protease inhibitor cocktails (Roche, Basel, Switzerland). The mixture was sonicated and centrifuged at 4°C and 30,000× g for 15 min. The supernatant was added to five volumes of cold acetone containing 10% (v/v) trichloroacetic acid, well mixed, and incubated at −20°C overnight. Afterward, the mixture was centrifuged again at 4°C, 30,000× g and the suspension was discarded. The sediment was washed three times with cold acetone, air-dried, and dissolved in RIPA buffer. The extracted proteins were quantified by a BCA kit (Thermo Fisher Scientific, United States); 300 µg proteins from each sample were mixed with sequencing-grade trypsin (Promega, Madison, WI, United States) at an enzyme to protein ratio of 1:50 and incubated for 16 h at 37°C. The peptides acquired from enzymatic digestion were dried by vacuum centrifugation.
iTRAQ Labeling and High-pH Reversed-phase Fractionation
The obtained peptides were processed by 8-plex iTRAQ reagent (AB Sciex, Framingham, MA, United States) per manufacturer’s instructions. The control samples were labeled with 115 iTRAQ tags, the model samples were labeled with 116 iTRAQ tags, and the GBFXD samples were labeled with 113 iTRAQ tags. High pH RP fractionation was performed with the U3000 HPLC chromatography system (Thermo Fisher Scientific, United States). The iTRAQ-labeled peptides were reconstituted with 100 µL high pH RP buffer A (98% H2O, 2% acetonitrile, pH 10.0) and loaded onto a C18 column (250 × 3 mm, 5 µm particle size, Phenomenex, United States). The column was eluted (flow rate: 0.2 ml/min) with the following parameters: 3–18% buffer B (2% H2O, 98% acetonitrile, pH 10.0) for 30 min, 18–32% B for 15 min, 32–98% B for 6 min, and 98% B for 15 min. The elution was monitored by detecting the absorbance at 214 nm.
LC-MS/MS Analysis
The peptides were re-dissolved in buffer A (2% acetonitrile, 0.1% formate) and centrifuged (4°C, 20,000× g, 10 min). Using the Nano LC System auto sampler (Thermo Fisher Scientific, United States), 10 µL peptides were loaded onto a 2 cm C18 trap column and eluted onto a 15 cm analytical C18 column with an inner diameter of 75 µm using a 3–55% buffer B (84% acetonitrile, 0.1% formate). The elution process lasted 112 min at a flow rate of 300 nL/min. With 2.2 kV voltage, the peptides were ionized by nano-electrospray ionization and data-dependent MS/MS was performed by the LTQ-Orbitrap XL mass spectrometer (Thermo Fisher Scientific). With the MS1 full scan (resolution: 60,000), the five most abundant precursor ions above the 5,000 counts threshold were chosen for MS/MS analysis along with a dynamic exclusion duration of 60 s. Normalized collision energy for high-energy collision dissociation was set to 40.0 and then the ion fragments were tested in the Orbitrap at a resolution of 7,500 with scan ranges of 350–1800 Da (MS1) and 100–1800 Da (MS2).
Identification and Quantification of Proteins
Raw data was analyzed by Proteome Discoverer (version: 2.0, Thermo Fisher Scientific, United States). The SEQUEST search engine was used to identify the proteins, and a mass tolerance of 10 ppm or 0.02 Da was permitted for intact peptides or fragmented ions with an allowance of no more than two missed cleavages with trypsin digestion. Except for iTRAQ labels, the fixed modification was carbamidomethylation (C), whereas the variable modifications were oxidation (M) and deamidation (NQ). With SEQUEST probability analysis, peptides at the 99% confidence interval were counted as identified, which may reduce the probability of false identifications. At least two unique peptides were required to quantify the proteins. To weigh and normalize the quantitative protein ratios, the median ratio in SEQUEST was applied. The DEPs were identified by fold change >1.5 or <0.67 with FDR < 1%.
Bioinformatic Analysis of DEPs
OmicsBean (http://www.omicsbean.cn), including gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, was used to analyze DEPs. In additon, the Heat map was produced by MetaboAnalyst 5.0 (https://www.metaboanalyst.ca).
Real-Time PCR
Total RNA in the lung was extracted by the Fast Pure Cell/Tissue Total RNA Isolation Kit (Vazyme Biotech Co., Ltd., Nanjing, China), from which cDNA was synthesized with the HiScript® II Q RT SuperMix for qPCR (Vazyme Biotech Co., Ltd., Nanjing, China) and quantified using the ChamQ universal SYBR qPCR Master Mix (Vazyme Biotech Co., Ltd., Nanjing, China). The results were analyzed using the QuantStudio 7 Flex detection system (Applied Biosystems Co., United States) and the 2−ΔΔCT method to assess the levels of mRNAs encoding targeted genes normalized to GAPDH. The primers are shown in Table 3.
Western Blotting
The proteins were extracted from the lungs in each group and measured by the BCA kit. Equal amounts of proteins were separated by SDS-PAGE, transferred to polyvinylidene difluoride membranes (0.45 μm, Millipore, United States), and incubated with primary antibodies including AKT (1:2,000 dilution), P-AKT (1:5,000 dilution), RTN4 (1:500 dilution) (Proteintech, Wuhan, China), Col1a2 (1:1,000 diluted), Col4a1 (1:1,000 dilution), PI3K (1:1,000 dilution) (Immunoway, Plano, TX, United States), P-PI3K (1:1,000 dilution), GAPDH (1:10,000 dilution (Affinity Biosciences, Cincinnati, OH, United States), PTEN (1:1,000 dilution) (Cell Signaling Technology Inc., United States) and Vimentin (1:1,000 dilution) (Arigo, Shanghai, China). The ChemiDoc™ MP Imaging system (Bio-Rad Co., United States) was used to quantify the protein bands with the Image Lab software (Bio-Rad Co., United States). The values were expressed as intensities normalized to GAPDH or their respective total proteins.
Statistical Analysis
Data were expressed as mean and standard deviation (SD). One-way analysis of variance (ANOVA) followed by the Dunnett’s post hoc test was used to determine statistical significances between groups. Differences with p < 0.05 were considered significant. All statistical analyses were performed in GraphPad Prism 6.0 (GraphPad Software Inc., San Diego, CA, United States).
Results
Identification of GBFXD Components
A total of 14 chemical compounds were identified in GBFXD, namely Prim-o-glucosylcimifugin, Calycosin-7-O-glucoside, Liquiritin, Lobetyolin, Isoliquiritin, Astragaloside II, Calycosin, Astragaloside IV, Astragaloside III, Glycyrrhizic acid, Magnolin, Astragaloside I, Schisandrol A, and Schisandrin A (Figures 1A,B and Table 2).
GBFXD Ameliorated Airway Inflammation and Remodeling in a Murine Model of CRA
HE and Sirius red staining were analyzed to determine airway inflammation and remodeling. The mice in the model group exhibited significantly increased inflammation and collagen deposition in the respiratory tract (Figures 2B–D). These pulmonary pathological symptoms were obviously alleviated in the mice with GBFXD treatment (24 g/kg or 36 g/kg) compared with those in the model group, although the administration of GBFXD (12 g/kg) and Montelukast did not improve collagen deposition (Figures 2B,D). For molecules relevant to airway remodeling, the expressions of Vimentin, α-SMA and HYP were increased, while the E-cadherin was reduced in CRA mice. More importantly, GBFXD (24 g/kg or 36 g/kg) drastically restored the expressions of those proteins (Figures 2E, 3A–D). In addition, the changes in Vimentin were further confirmed by Western blotting (Supplementary Figure S1). The above results confirmed that GBFXD could relieve airway inflammation and remodeling in the murine model of CRA.
FIGURE 3. GBFXD decreased epithelial-to-mesenchymal transition in CRA mice. (A) Markers of lung EMT were performed by immunohistochemical staining (×200), Scale bar: 100 μm. (B–D) Relative protein expressions of E-cadherin (epithelial mark), Vimentin and α-SMA (mesenchymal marks) were quantified by AOD. Values are presented as mean ± SD (n = 3 mice per group). #p < 0.05 compared to control group, and *p < 0.05 compared to model group.
DEPs Identification by iTRAQ-Based Proteomics
To explore the underlying mechanisms of GBFXD in CRA, iTRAQ-based proteomics was used to identify the DEPs in the control, model and GBFXD (36 g/kg) groups. In addition, the dose of GBFXD is the same as before (Liu et al., 2019; Xing et al., 2019). A total of 2,319 proteins were identified, among which 1,472 were quantified. As shown in Figures 4A,B 79 DEPs in the model/control were identified (Table 4, upregulated: 23, downregulated: 56), whereas 24 DEPs were identified in the GBFXD/model group (Table 5, upregulated: 2, downregulated: 22). There were 18 overlapping DEPs, which were upregulated in the model/control group, and downregulated by GBXFD treatment (Figure 4C). Further analysis found that these proteins were mostly collagen and laminin, which were extracellular matrix (ECM) proteins and closely related to airway remodeling. Downregulated expression of fibronectin (Fn1), a key marker of airway remodeling, was also observed after GBFXD treatment.
FIGURE 4. The differentially expressed proteins (DEPs) in lung between Model/Control and GBFXD/Model. (A) Venn-diagram of DAPs. (B) The left was DAPs involved in Model/Control and the right was GBFXD/Model (the numbers of upregulated, downregulated proteins characterized by Red and Green, respectively). (C) Heat map of overlap DEPs.
Bioinformatic Analysis of DEPs
GO annotations including biological process (BP), cellular component (CC), and molecular function (MF) were used to obtain the functional information of DEPs. As shown in Figures 5A,B, the DEPs of the Model/control group indicated significant enrichment of 700 BP terms, 106 CC terms, and 97 MF terms. After GBFXD treatment, the DEPs of the GBFXD/Model group showed significant enrichment of 621 BP terms, 70 CC terms, and 32 MF terms. The top 10 significantly enriched GO terms were shown in Figures 5C,D.
FIGURE 5. Functional annotation and categories of DAPs. (A–B) Bioinformatics analysis of DEPs by categories: BP, CC, MF, and KEGG pathway, respectively (A: Model/Control, B: GBFXD/Model). Counts for each category represent the total associated terms in the database with the query protein list. Terms are statistically significant with a p-value <0.05. (C–D) The ten most significantly enriched terms of the GO hierarchy, information of percentage, and number of involved proteins in a term are shown in the left and right y-axes, respectively (C: Model/Control, D: GBFXD/Model). (E–F) Enriched KEGG pathways of DAPs (E: Model/Control, F: GBFXD/Model). The number of involved proteins and corresponding p-values are shown on the right side of the column.
The top three terms of DEPs for BP from the model/control group were muscle contraction, ECM organization, and cell adhesion. Most proteins were involved in contractile fiber, ECM component, collagen trimer for CC, and receptor binding, ECM structural constituent, cell adhesion molecule binding for MF (Figure 5C). The most enriched BP terms of DEPs for BP from the GBFXD/model group include cell adhesion, ECM organization, and extracellular structure organization (Figure 5D). In addition, the ECM, ECM component, basement membrane for CC and receptor binding, cell adhesion molecule binding, structural molecule activity for MF were the most enriched. Consequently, the GO analysis indicated that GBFXD may impact ECM related to airway remodeling.
According to KEGG pathway enrichment analysis, the DEPs from the model/control and GBFXD/model groups were shown to be simultaneously involved in ECM-receptor interactions, the PI3K/AKT signaling pathway, and focal adhesion, which were the top three significantly enriched pathways (Figures 5E,F) and have been reported to participate in airway remodeling in asthma (Athari, 2019; Yap et al., 2019; Zhou et al., 2019; Ali et al., 2020).
Validation of DEPs
The expression of mRNAs encoding collagen (Cola1, Col1a2, Col4a1, Col4a2), laminin (LAMA3, LAMA5, LAMB2, LAMB3), and Fn1 in the lungs of the model group were significantly higher than those in the control group (Figures 6A–C). Meanwhile, the expression levels of these mRNAs were decreased in mice treated with GBFXD or Montelukast. Data on the protein levels of Col1a2, Col4a1, and RTN4 also showed consistent results (Figures 6D–G) with iTRAQ proteomic analysis and the trends in the mRNA levels.
FIGURE 6. Validation of DEPs identified in proteomics. (A–C) The mRNA levels of Collagen (Col1a1, Col1a2, Col4a1, Col4a2), Laminin (LAMA3, LAMA5, LAMB2, LAMB3) and Fn1 in lung were measured by Real-time PCR. (D) The proteins of Col1a2, Col4a1, and RTN4 were analyzed by Western Blot. (E–G) Quantification of Col1a2, Col4a1, and RTN4. Values are presented as mean ± SD (n = 3 or 4 mice per group). #p < 0.05 compared to control group, and *p < 0.05 compared to model group.
GBFXD Probably Inhibited the Activation of the PI3K/AKT Signaling Pathway Mediated by RTN4
We further confirmed the pharmacological mechanisms of GBFXD. Recent publications suggested that the RTN4 promoted angiogenesis in proliferative diabetic retinopathy via PI3K/Akt signal pathway and epithelial-mesenchymal transition in non-small cell lung cancer cells (Zhang et al., 2017; Wu et al., 2018). Based on our proteomics results, GBFXD could regulate PI3K/AKT pathway and RTN4 was also an overlapping important DEP. In addition, our verification experiments also proved that GBFXD could regulate the expression of RTN4. Therefore, we hypothesized that GBFXD could alleviate airway remodeling by inhibiting the pathway of PI3K/AKT via RTN4. The expression of PTEN was lower, and the phosphorylation of PI3K and AKT were remarkably higher in the model group (Figures 7A–D) compared to the control group. More importantly, the dysregulation was reversed after GBFXD treatment. Consequently, the data indicated that GBFXD ameliorated airway remodeling by inhibiting the activation of the PI3K/AKT signaling pathway via RTN4 to reduce ECM deposition possibly.
FIGURE 7. Effect of GBFXD on PI3K/AKT signaling pathway. (A) Western blotting was used to detect the protein expressions of PTEN, P-PI3K, PI3K, P-AKT, and AKT. (B–D) Quantification of PTEN, P-PI3K/PI3K and P-AKT/AKT. Values are presented as mean ± SD (n = 3 or 4 mice per group). #p < 0.05 compared to control group, and *p < 0.05 compared to model group.
Discussion
Asthma is a common chronic inflammatory respiratory disease of highly prevalence worldwide. There is currently no treatment to cure this disease, but some patients may spontaneously enter asthma remission later in life (Carpaij et al., 2017). Asthma remission often refers to the state in which patients are no longer burdened by the disease to some extent, and no longer require any asthma medication. However, most of the patients might still have airway remodeling, hyperresponsiveness, or declined lung function, and very few can go into complete asthma remission with no respiratory symptoms (Broekema et al., 2011; Carpaij et al., 2017). Even though some therapies such as ICS or biological treatment can fully control the symptoms, they have not been associated with real remission. Recently, immunotherapy has been shown to have a vital effect on asthma remission and cumulating evidence indicates that TCM has unique advantages in treating asthma primarily by regulating immunity (Bao et al., 2019; Cui et al., 2019; Yao et al., 2019). GBFXD was shown to be effective in improving the life of asthmatic children and reducing the rate of relapses. Moreover, findings from our previous study have indicated that GBFXD could regulate Th1/Th2, macrophage polarization (Liu et al., 2019), CD4+/CD8+ (Lu et al., 2016) and Treg cells (Dong et al., 2020) to enhance the immune activity, which has triggered our increasing interest in further research on GBFXD.
Asshown in the results of the HE and Sirius red staining, airway inflammatory infiltration and collagen deposition were effectively relieved after GBFXD treatment. Airway remodeling is closely related to the dysregulation of epithelial-mesenchymal transition (EMT) by which epithelial cells transform into the mesenchymal phenotype after epithelial damage (Guan et al., 2020). EMT is mainly associated with decreased levels of epithelial phenotypic proteins including E-cadherin, ZO-1, and increased levels of mesenchymal phenotypic proteins such as N-cadherin, Vimentin, and α-SMA (Bartis et al., 2014). Subsequently, we showed the pharmacological effects of GBFXD on typical markers of EMT (Vimentin, α-SMA, and E-cadherin). HYP, an indicator of the collagen content in the lung tissue, is associated with the occurrence of lung fibrosis (Sun et al., 2020). In this study, HYP levels were significantly reduced in mice treated with GBFXD, which further confirmed the effect of GBFXD on airway remodeling.
However, the detailed mechanisms of GBFXD on alleviating airway remodeling in CRA remained unclear, which were further explored by proteomic analyses. The identified overlapping DEPs were mostly collagen and laminin. Reticular basement membrane with ECM protein deposition is a characteristic of airway remodeling. The ECM proteins include structural proteins, adhesion proteins, glycosaminoglycans, and proteoglycans, and examples include collagen, laminin, and Fn1 (Olczyk et al., 2014). More importantly, the deposition of ECM molecules is increased in asthmatic airways (Royce et al., 2012; Mostaço-Guidolin et al., 2019). In this study, increased mRNA or protein levels of DEPs were observed for collagen (Col1a1, Col1a2, Col4a1, Col4a2), laminin (LAMA3, LAMA5, LAMB2, LAMB3), and Fn1 in CRA mice, while the expressions of these molecules decreased after GBFXD treatment. Collectively, our validation of DEPs was consistent with proteomic findings.
The KEGG analysis showed that GBFXD could regulate pathways related to airway remodeling, such as ECM-receptor interactions, focal adhesion, and the PI3K/AKT signaling pathway, which were the top three significantly enriched pathways containing the most DEPs for both Model/Control and GBFXD/Model. RTN4 was an important common DEP and downregulated after GBFXD treatment. A previous study showed that RTN4 was primarily localized in the ER and was commonly known as Nogo-B in the peripheral tissues including the lung (Pathak et al., 2018). Accumulating research also reported the importance of RTN4 in tissue repairing and hepatic fibrosis (Yu et al., 2009; Tashiro et al., 2013). In addition, RTN4 is necessary for the migration and contraction of airway smooth muscle cells in airway remodeling in chronic asthma (Xu et al., 2011), promoting angiogenesis in proliferative diabetic retinopathy via the PI3K/Akt signal pathway (Zhang et al., 2017), and facilitating EMT in non-small cell lung cancer (Wu et al., 2018). However, the relationship of RTN4 and PI3K/AKT pathway with airway remodeling in CRA was still unclear. According to our present findings (Figures 6D–G, 7A–D), GBFXD regulated RTN4 and suppressed the activation of the PI3K/AKT pathway. Recent research has demonstrated an association between the PI3K/AKT pathway and airway remodeling (Hong et al., 2017; Shao et al., 2019) and the KEGG analysis of our proteomics also confirmed this finding. Consequently, we hypothesize that potential mechanisms of GBFXD on airway remodeling in CRA may be related to inhibiting RTN4 and subsequently suppressing the activation of the PI3K/AKT signaling pathway to reduce the deposition of ECM related molecules. More in-depth research will be performed to prove the conclusion. In addition, GBFXD in CPA mainly regulated the pathways of mitochondrial energy metabolism and macrophage polarization which were also related to airway inflammation and remodeling as reported in our previous study (Liu et al., 2019). Herein, there may be different pathological mechanisms between CPA and CRA, which may explain the differences in HE, PAS staining and Penh between CPA and CRA mice in previous studies (Liu et al., 2019) and warrant further exploration. The proteomics data have also confirmed that GBFXD could regulate airway remodeling in CPA or CRA even though through the different pathways, suggesting multiple targets for the treatment of asthma.
In conclusion, our data demonstrated the effect of GBFXD on airway remodeling by regulating the ECM proteins in CRA mice. The potential mechanisms may be related to inhibiting RTN4 and subsequently suppressing the activation of the PI3K/AKT signaling pathway to reduce the deposition of ECM related molecules. This study provides a new insight into the mechanism of GBFXD on improving airway remodeling for the treatment of CRA.
Data Availability Statement
The mass spectrometry proteomics data have been deposited to the ProteomeX change Consortium via the PRIDE partner repository with the dataset identifier PXD020711.
Ethics Statement
All experimental procedures were performed and approved by the animal ethics committee of Nanjing University of Chinese medicine (No.201809A019, No.201810A021).
Author Contributions
QQ and YN performed the data analysis and wrote the manuscript, JJ provided information for research design, YH, YM, LS, YY, and ST provided experimental support, ZX designed this study and edited the final manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (Grant No. 81774367 and Grant No. 82074489), the leading academics training program of Chinese medicine in Jiangsu Province, China (No. SLJ0224), Postgraduate Research and Practice Innovation Program of Jiangsu Province, China (Grant No. KYCX20_1468), and the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD), the Open Projects of the Discipline of Chinese Medicine of Nanjing University of Chinese Medicine Supported by the Subject of Academic priority discipline of Jiangsu Higher Education Institutions (No. ZYX03KF).
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2021.588588/full#supplementary-material.
References
Ali, M. K., Kim, R. Y., Brown, A. C., Mayall, J. R., Karim, R., Pinkerton, J. W., et al. (2020). Crucial role for lung iron level and regulation in the pathogenesis and severity of asthma. Eur. Respir. J. 55, 1901340. doi:10.1183/13993003.01340-2019
Araujo, B. B., Dolhnikoff, M., Silva, L. F. F., Elliot, J., Lindeman, J. H. N., Ferreira, D. S., et al. (2008). Extracellular matrix components and regulators in the airway smooth muscle in asthma. Eur. Respir. J. 32, 61–69. doi:10.1183/09031936.00147807
Athari, S. S. (2019). Targeting cell signaling in allergic asthma. Sig Transduct Target. Ther. 4, 45. doi:10.1038/s41392-019-0079-0
Bao, K., Yuan, W., Zhou, Y., Chen, Y., Yu, X., Wang, X., et al. (2019). A Chinese prescription yu-ping-feng-san administered in remission restores bronchial epithelial barrier to inhibit house dust mite-induced asthma recurrence. Front. Pharmacol. 10, 1698. doi:10.3389/fphar.2019.01698
Baraldo, S., Turato, G., Bazzan, E., Ballarin, A., Damin, M., Balestro, E., et al. (2011). Noneosinophilic asthma in children: relation with airway remodelling. Eur. Respir. J. 38, 575–583. doi:10.1183/09031936.00168210
Bartis, D., Mise, N., Mahida, R. Y., Eickelberg, O., and Thickett, D. R. (2014). Epithelial-mesenchymal transition in lung development and disease: does it exist and is it important? Thorax 69, 760–765. doi:10.1136/thoraxjnl-2013-204608
Benayoun, L., Druilhe, A., Dombret, M.-C., Aubier, M., and Pretolani, M. (2003). Airway structural alterations selectively associated with severe asthma. Am. J. Respir. Crit. Care Med. 167, 1360–1368. doi:10.1164/rccm.200209-1030OC
Bonato, M., Bazzan, E., Snijders, D., Tinè, M., Biondini, D., Turato, G., et al. (2018). Clinical and pathologic factors predicting future asthma in wheezing children. A longitudinal study. Am. J. Respir. Cel Mol Biol. 59, 458–466. doi:10.1165/rcmb.2018-0009OC
Broekema, M., Timens, W., Vonk, J. M., Volbeda, F., Lodewijk, M. E., Hylkema, M. N., et al. (2011). Persisting remodeling and less airway wall eosinophil activation in complete remission of asthma. Am. J. Respir. Crit. Care Med. 183, 310–316. doi:10.1164/rccm.201003-0494OC
Camoretti-Mercado, B., and Lockey, R. F. (2020). Bitter taste receptors in the treatment of asthma: opportunities and challenges. J. Allergy Clin. Immunol. 146, 776. doi:10.1016/j.jaci.2020.04.036
Carpaij, O. A., Burgess, J. K., Kerstjens, H. A. M., Nawijn, M. C., and van den Berge, M. (2019). A review on the pathophysiology of asthma remission. Pharmacol. Ther. 201, 8–24. doi:10.1016/j.pharmthera.2019.05.002
Carpaij, O. A., Nieuwenhuis, M. A. E., Koppelman, G. H., van den Berge, M., Postma, D. S., and Vonk, J. M. (2017). Childhood factors associated with complete and clinical asthma remission at 25 and 49 years. Eur. Respir. J. 49, 1601974. doi:10.1183/13993003.01974-2016
Chakir, J., Shannon, J., Molet, S., Fukakusa, M., Elias, J., Laviolette, M., et al. (2003). Airway remodeling-associated mediators in moderate to severe asthma: effect of steroids on TGF-β, IL-11, IL-17, and type I and type III collagen expression. J. Allergy Clin. Immunol. 111, 1293–1298. doi:10.1067/mai.2003.1557
Chen, Q. (2011). Research methods in Pharmacology of Chinese Materia Medica (version 3). Beijing, China: People's Medical Publishing House
Costello, R. W., and Cushen, B. (2020). Looking back to go forward: adherence to inhaled therapy before biologic therapy in severe asthma. Eur. Respir. J. 55, 2000954. doi:10.1183/13993003.00954-2020
Croisant, S. (2014). Epidemiology of asthma: prevalence and burden of disease. Adv. Exp. Med. Biol. 795, 17–29. doi:10.1007/978-1-4614-8603-9_2
Cui, J., Xu, F., Tang, Z., Wang, W., Hu, L. L., Yan, C., et al. (2019). Bu-Shen-Yi-Qi formula ameliorates airway remodeling in murine chronic asthma by modulating airway inflammation and oxidative stress in the lung. Biomed. Pharmacother. 112, 108694. doi:10.1016/j.biopha.2019.108694
Dong, Y., Yan, H., Zhao, X., Lin, R., Lin, L., Ding, Y., et al. (2020). Gu-ben-fang-xiao decoction ameliorated murine asthma in remission stage by modulating microbiota-acetate-tregs Axis. Front. Pharmacol. 11, 549. doi:10.3389/fphar.2020.00549
Grainge, C. L., Lau, L. C. K., Ward, J. A., Dulay, V., Lahiff, G., Wilson, S., et al. (2011). Effect of bronchoconstriction on airway remodeling in asthma. N. Engl. J. Med. 364, 2006–2015. doi:10.1056/NEJMoa1014350
Guan, R., Wang, J., Cai, Z., Li, Z., Wang, L., Li, Y., et al. (2020). Hydrogen sulfide attenuates cigarette smoke-induced airway remodeling by upregulating SIRT1 signaling pathway. Redox Biol. 28, 101356. doi:10.1016/j.redox.2019.101356
Halwani, R., Sultana, A., Al-Kufaidy, R., Jamhawi, A., Vazquez-Tello, A., and Al-Muhsen, S. (2016). Th-17 regulatory cytokines inhibit corticosteroid induced airway structural cells apoptosis. Respir. Res. 17, 6. doi:10.1186/s12931-015-0307-2
Hong, W., Peng, G., Hao, B., Liao, B., Zhao, Z., Zhou, Y., et al. (2017). Nicotine-induced airway smooth muscle cell proliferation involves TRPC6-dependent calcium influx via α7 nAChR. Cell. Physiol Biochem. 43, 986–1002. doi:10.1159/000481651
Hough, K. P., Curtiss, M. L., Blain, T. J., Liu, R.-M., Trevor, J., Deshane, J. S., et al. (2020). Airway remodeling in asthma. Front. Med. 7, 191. doi:10.3389/fmed.2020.00191
Huang, Z., Gao, L., Zhao, X., Ling, H., and Chen, W. (2016). Effect of Gubenfangxiao decoction on respiratory syncytial virus-induced asthma and expression of asthma susceptibility gene orosomucoid 1-like protein 3 in mice. J. Tradit Chin. Med. 36, 101–106. doi:10.1016/s0254-6272(16)30015-2
Jendzjowsky, N. G., and Kelly, M. M. (2019). The role of airway myofibroblasts in asthma. Chest 156, 1254–1267. doi:10.1016/j.chest.2019.08.1917
Lazarinis, N., Bood, J., Gomez, C., Kolmert, J., Lantz, A.-S., Gyllfors, P., et al. (2018). Leukotriene E4 induces airflow obstruction and mast cell activation through the cysteinyl leukotriene type 1 receptor. J. Allergy Clin. Immunol. 142, 1080–1089. doi:10.1016/j.jaci.2018.02.024
Liu, L.-w., Xing, Q.-q., Zhao, X., Tan, M., Lu, Y., Dong, Y.-m., et al. (2019). Proteomic analysis provides insights into the therapeutic effect of GU-BEN-FANG-XIAO decoction on a persistent asthmatic mouse model. Front. Pharmacol. 10, 441. doi:10.3389/fphar.2019.00441
Lu, Y., Xu, J.-Y., Zhang, X.-H., and Zhao, X. (2016). Gu-Ben-Fang-Xiao decoction attenuates sustained airway inflammation by suppressing ER stress response in a murine asthma remission model of respiratory syncytial virus infection. J. Ethnopharmacology 192, 496–509. doi:10.1016/j.jep.2016.09.039
Mostaço-Guidolin, L. B., Osei, E. T., Ullah, J., Hajimohammadi, S., Fouadi, M., Li, X., et al. (2019). Defective fibrillar collagen organization by fibroblasts contributes to airway remodeling in asthma. Am. J. Respir. Crit. Care Med. 200, 431–443. doi:10.1164/rccm.201810-1855OC
Ojiaku, C. A., Yoo, E. J., and Panettieri, R. A. (2017). Transforming growth factor β1 function in airway remodeling and hyperresponsiveness. The missing link? Am. J. Respir. Cel Mol Biol 56, 432–442. doi:10.1165/rcmb.2016-0307TR
Olczyk, P., Mencner, Ł., and Komosinska-Vassev, K. (2014). The role of the extracellular matrix components in cutaneous wound healing. Biomed. Res. Int. 2014, 1. doi:10.1155/2014/747584
Pathak, G. P., Shah, R., Kennedy, B. E., Murphy, J. P., Clements, D., Konda, P., et al. (2018). RTN4 knockdown dysregulates the AKT pathway, destabilizes the cytoskeleton, and enhances paclitaxel-induced cytotoxicity in cancers. Mol. Ther. 26, 2019–2033. doi:10.1016/j.ymthe.2018.05.026
Prakash, Y. S., Halayko, A. J., Gosens, R., Panettieri, R. A., Camoretti-Mercado, B., and Penn, R. B. (2017). An official American thoracic society research statement: current challenges facing research and therapeutic advances in airway remodeling. Am. J. Respir. Crit. Care Med. 195, e4–e19. doi:10.1164/rccm.201611-2248ST
Royce, S. G., Cheng, V., Samuel, C. S., and Tang, M. L. K. (2012). The regulation of fibrosis in airway remodeling in asthma. Mol. Cell Endocrinol. 351, 167–175. doi:10.1016/j.mce.2012.01.007
Shao, Y., Chong, L., Lin, P., Li, H., Zhu, L., Wu, Q., et al. (2019). MicroRNA-133a alleviates airway remodeling in asthtama through PI3K/AKT/mTOR signaling pathway by targeting IGF1R. J. Cel Physiol 234, 4068–4080. doi:10.1002/jcp.27201
Sun, B., Shi, Y., Li, Y., Jiang, J., Liang, S., Duan, J., et al. (2020). Short-term PM2.5 exposure induces sustained pulmonary fibrosis development during post-exposure period in rats. J. Hazard. Mater. 385, 121566. doi:10.1016/j.jhazmat.2019.121566
Tao, Z., Meng, X., Han, Y.-q., Xue, M.-m., Wu, S., Wu, P., et al. (2017). Therapeutic mechanistic studies of ShuFengJieDu capsule in an acute lung injury animal model using quantitative proteomics Technology. J. Proteome Res. 16, 4009–4019. doi:10.1021/acs.jproteome.7b00409
Tashiro, K., Satoh, A., Utsumi, T., Chung, C., and Iwakiri, Y. (2013). Absence of Nogo-B (reticulon 4B) facilitates hepatic stellate cell apoptosis and diminishes hepatic fibrosis in mice. Am. J. Pathol. 182, 786–795. doi:10.1016/j.ajpath.2012.11.032
van den Toorn, L. M., Overbeek, S. E., de Jongste, J. C., Leman, K., Hoogsteden, H. C., and Prins, J. B. (2001). Airway inflammation is present during clinical remission of atopic asthma. Am. J. Respir. Crit. Care Med. 164, 2107–2113. doi:10.1164/ajrccm.164.11.2006165
Voskamp, A. L., Kormelink, T. G., van Wijk, R. G., Hiemstra, P. S., Taube, C., de Jong, E. C., et al. (2020). Modulating local airway immune responses to treat allergic asthma: lessons from experimental models and human studies. Semin. Immunopathol. 42, 95–110. doi:10.1007/s00281-020-00782-4
Wu, D., Zhao, B., Qi, X., Peng, F., Fu, H., Chi, X., et al. (2018). Nogo-B receptor promotes epithelial-mesenchymal transition in non-small cell lung cancer cells through the Ras/ERK/Snail1 pathway. Cancer Lett. 418, 135–146. doi:10.1016/j.canlet.2018.01.030
Xing, Q.-Q., Liu, L.-W., Zhao, X., Lu, Y., Dong, Y.-M., and Liang, Z.-Q. (2019). Serum proteomics analysis based on label-free revealed the protective effect of Chinese herbal formula Gu-Ben-Fang-Xiao. Biomed. Pharmacother. 119, 109390. doi:10.1016/j.biopha.2019.109390
Xu, W., Hong, W., Shao, Y., Ning, Y., Cai, Z., and Li, Q. (2011). Nogo-B regulates migration and contraction of airway smooth muscle cells by decreasing ARPC 2/3 and increasing MYL-9 expression. Respir. Res. 12, 14. doi:10.1186/1465-9921-12-14
Yao, L., Wang, S., Wei, P., Bao, K., Yuan, W., Wang, X., et al. (2019). Huangqi-Fangfeng protects against allergic airway remodeling through inhibiting epithelial-mesenchymal transition process in mice via regulating epithelial derived TGF-β1. Phytomedicine 64, 153076. doi:10.1016/j.phymed.2019.153076
Yap, H. M., Israf, D. A., Harith, H. H., Tham, C. L., and Sulaiman, M. R. (2019). Crosstalk between signaling pathways involved in the regulation of airway smooth muscle cell hyperplasia. Front. Pharmacol. 10, 1148. doi:10.3389/fphar.2019.01148
Yu, J., Fernández-Hernando, C., Suarez, Y., Schleicher, M., Hao, Z., Wright, P. L., et al. (2009). Reticulon 4B (Nogo-B) is necessary for macrophage infiltration and tissue repair. Proc. Natl. Acad. Sci. 106, 17511–17516. doi:10.1073/pnas.0907359106
Yuan, X. J., Sun, Y. Q., Wang, S. M., Li, Y. N., Wang, M. Q., and Liu, X. R. (2010). Clinical research of Gubenfangxiao decoction combined with point application on chronic asthmatic children in 100 cases. Zhong Hua Zhong Yi Yao Za Zhi 21, 2306–2309.
Zhang, Y., Wang, L., Zhang, Y., Wang, M., Sun, Q., Xia, F., et al. (2017). Nogo-B promotes angiogenesis in proliferative diabetic retinopathy via VEGF/PI3K/akt pathway in an autocrine manner. Cell Physiol Biochem 43, 1742–1754. doi:10.1159/000484061
Keywords: GU-BEN-FANG-XIAO decoction, airway remodeling, extracellular matrix deposition, chronic remission, asthma
Citation: Xing Q, You Y, Zhao X, Ji J, Yan H, Dong Y, Ren L, Ding Y and Hou S (2021) iTRAQ-Based Proteomics Reveals Gu-Ben-Fang-Xiao Decoction Alleviates Airway Remodeling via Reducing Extracellular Matrix Deposition in a Murine Model of Chronic Remission Asthma. Front. Pharmacol. 12:588588. doi: 10.3389/fphar.2021.588588
Received: 29 July 2020; Accepted: 01 March 2021;
Published: 14 June 2021.
Edited by:
Adolfo Andrade-Cetto, National Autonomous University of Mexico, MexicoReviewed by:
Yu Chiang Hung, Kaohsiung Chang Gung Memorial Hospital, TaiwanChia-Hung Yen, Kaohsiung Medical University, Taiwan
Copyright © 2021 Xing, You, Zhao, Ji, Yan, Dong, Ren, Ding and Hou. 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.
*Correspondence: Xia Zhao, emhhb3hpYWh5QG5qdWNtLmVkdS5jbg==
†These authors have contributed equally to the study