Corrigendum: Total flavonoids of Rhizoma Drynariae ameliorate bone growth in experimentally induced tibial dyschondroplasia in chickens via regulation of OPG/RANKL axis
- 1College of Veterinary Medicine, Henan Agricultural University, Zhengzhou, China
- 2Sanquan College of Xinxiang Medical University, Xinxiang, China
- 3National Center for International Research on Animal Genetics, Breeding and Reproduction (NCIRAGBR), Huazhong Agricultural University, Wuhan, China
- 4Livestock and Dairy Development Department, Pishin, Pakistan
Background: Rhizoma Drynariae, traditional Chinese herb, is widely used to treat and prevent bone disorders. However, experimental evidence on the use of Rhizoma Drynariae extract, total flavonoids of Rhizoma Drynariae (TFRD) to treat tibial dyschondroplasia (TD) in chickens and its underlying mechanisms have not been investigated.
Purpose: To evaluate the therapeutic effect of TFRD on leg disease caused by TD and elucidate its mechanisms in modulating the bone status.
Methods: Thiram-induced chicken TD model has been established. The tibia status was evaluated by analyzing tibia-related parameters including tibial weight, tibial length and its growth plate width and by performing histopathological examination. The expression of tibial bone development-related genes and proteins was confirmed by western blotting and qRT-PCR.
Results: The results showed that administration of TFRD mitigated lameness, increased body weight, recuperated growth plate width in broilers affected by TD and the increase of tibia weight and tibia length is significantly positively correlated with body weight. Compared with the TD group broilers, 500 mg/kg TFRD evidently reduced the damage width of the growth plate and improved its blood vessel distribution by elevating the gene expression levels of BMP-2 and Runx2 and OPG/RANKL ratio. Furthermore, correlation analysis found that the damage width of the growth plate was negatively correlated with the expression levels of BMP-2 and OPG.
Conclusion: The present study revealed that TFRD could promote the bone growth via upregulating OPG/RANKL ratio, suggesting that TFRD might be a potential novel drug in the treatment of TD in chickens.
Introduction
Tibial dyschondroplasia (TD) is a commonly occurring leg problem in fast growing birds caused by vascular injury or ischemia resulting in abnormal chondrocyte differentiation in the proximal tibial growth plate, in which leads to the claudication and brittle fracture (Huang et al., 2017, 2018; Zhang et al., 2018). The development of growth plate (GP) includes integrated processes of chondrocyte differentiation, proliferation, mineralization and osteogenesis (Lui et al., 2014). Among them, the vascularization of the GP is essential for the proliferation and differentiation of chondrocytes to promote bone mineralization (Huang et al., 2018). However, it is noted that non-vascularized, non-mineralized white opaque masses are generally observed during GP development in TD broilers (Zhang et al., 2018). TD broilers are usually characterized by difficulty in standing, inflexible gait, limited intake of food and water, retarded growth and slow weight gain, which cause a decrease in production performance (Huang et al., 2018).
The broiler chickens with TD can also reduce the meat yield and quality, thereby causing huge economic losses to the poultry industry (Huang et al., 2019a; Huang et al., 2021a). Although there is a lot of evidence that angiogenesis dysfunction (Huang et al., 2018; Huang et al., 2019b), genetic selection (Kapell et al., 2017), manganese deficiency (Dong et al., 2021), calcium and phosphorus metabolism disorders (Ledwaba and Roberson 2003), and dietary dithiocarbamates (Rath et al., 2007; Huang et al., 2019a) plays a role in the establishment of TD lesions. However, the cause and pathogenesis of TD lesion formation in chickens and how to prevent and minimize its harmful effects on the poultry industry remains the challenge for future research.
Previous studies demonstrated that the occurrence of TD in broilers is related to the abnormal development of tibial bone (Huang et al., 2017). Bone morphogenetic protein-2 (BMP-2) is an important factor regulating bone formation that promotes the differentiation of bone marrow mesenchymal stem cells into osteoblasts (Song et al., 2014). Moreover, many studies have revealed that most members of BMP can promote the growth and development of osteoblasts (Feng and Derynck, 2005; Bruderer et al., 2014). BMP-2 binds to cell surface receptors during bone formation and activates downstream signaling molecules called Smads, Smads (including Smad1, Smad5 and Smad8), which are activated and combined with Smad4 to form Smads complex. It is then transferred into the nucleus as transcriptional enhancers by activating runt-related transcription factor-2 (Runx2) and promoting osteoblast differentiation and maturation (Nuschke et al., 2014; Pierrefite-Carle et al., 2015). In recent years, osteoprotegerin (OPG)/receptor activator of nuclear factor kappa-B ligand (RANKL) axis has been extensively studied that the high expression of OPG and low expression of RANKL can effectively promote the expression of osteoblasts and accelerate the recovery of bone in through (Ibáñez, et al., 2019; Amin, et al., 2020; Ma, et al., 2020), which may be through the combination of OPG and the key factor RANKL in the process of osteoclast differentiation to inhibit the combination of RANKL and RANK, thereby blocking the differentiation of osteoclast precursors into osteoclasts (Udagawa et al., 2021). However, the underlying mechanism of OPG/RANKL axis to regulate tibial development in TD broilers remains largely unknown.
Recently, the incidence of broiler leg disease can be reduced by immunization, reasonably increasing stocking density, nutrition management and light time (Garner et al., 2012; Kierończyk et al., 2017), which is not expected to increase production costs and labor. In addition, it is essential to select low-residue drugs for prevention and treatment due to the short growth cycle and long drug metabolism time of broilers. Traditional Chinese medicine (TCM) has a long history, less harmful ingredients, fewer side effects, lesser residues, and is relatively safe (Zhang et al., 2016, 2018; Sun et al., 2020). Therefore, TCM may become the first choice for the prevention and treatment of various bone disorders including TD (Yao et al., 2018). Total flavonoids from Rhizoma Drynariae (TFRD), a Chinese herbal medicine product was extracted from the root of Drynaria roosii Nakaike (Chen et al., 2021). TFRD can replenish the kidney, strengthen the bones, promote the healing of fracture, and relieve pain. It has been proved to be useful in the treatment and prevention of bone disorders by changing the function of articular cartilage cells, increasing the proliferation and differentiation of osteoblasts, and promoting bone strength (Song et al., 2016; Zhang et al., 2017; Yao et al., 2018). Therefore, the present study aimed to investigate whether TFRD has the therapeutic and protective effects on tibial development of TD broiler chickens via regulating the OPG/RANKL axis.
Materials and Methods
Animal Ethics and Experimental Materials
All experimental procedures were conducted strictly following the Guidelines of the Care and Use of Laboratory Animals in China. This study was approved by the ethics committee of Henan Agricultural University, Zhengzhou, China (Permit No: 17-0126).
All Arbor Acres (AA) chickens (regardless of gender; 1-day-old; 48.25 ± 1.15 g) were obtained from a commercial hatchery (Xingda Poultry Industry Co., Ltd., Kaifeng, China). Thiram (AR, purity ≥98%, #C10036460) was acquired from Macklin Biochemical Co., Ltd. (Shanghai, China). TFRD (#K198752, purity ≥95%) was detected by HPLC and purchased from Xi’an Kailai Biological Engineering Co., Ltd. (Xi’an, China).
Experimental Design
All broiler chickens were divided randomly into six separate groups (n = 40/group, 10 chicks per replicate and four replicates per group): the control group (Con), thiram-induced TD group (TD), TFRD prevention group (LTFRD), dose gradient TFRD-treated group (TD + LTFRD, TD + MTFRD and TD + HTFRD). The management protocols used for the broilers in each group was as follows: Broilers in the Con group were fed a normal diet, the LTFRD group were fed a normal diet with 125 mg/kg TFRD, and the experimental groups (TD, TD + LTFRD, TD + MTFRD, and TD + HTFRD) were fed the same diet as the Con group supplemented with 100 mg/kg thiram from day 4 to day 7. For TFRD-treated groups including TD + LTFRD, TD + MTFRD, and TD + HTFRD group, all broiler chickens were administered with 125 mg/kg, 250 mg/kg, 500 mg/kg of the TFRD in fed during the entire experimental period. The dosage of TFRD were based on previous reports (Yao et al., 2018; Lv et al., 2021). All broilers were housed in a controlled environment maintained at the recommended standard room temperature (23–35°C), relative humidity (60–70%), 12:12 h dark or light cycle, ventilation, and hygienic conditions as required for broilers.
Sample Collection and Analysis of Tibia-Related Parameters
All groups were reared for 21 days, and daily body weight gain in every group was calculated. Ten chicks per group were randomly selected and euthanized on days 7, 14, and 21. Then, the removal of the shank muscles and the exposure of the tibia growth plate region was performed. In addition, the bone histomorphometric parameters, such as the weight, length, and the tibia growth plate (TGP) width were determined by an electronic balance with a sensitivity of 0.001 g (#FA1204N, Jinghai Instrument Co., Ltd., Shanghai, China) and digital calipers (#SATA91511, TATA Company, Shanghai, China). Then the tibia bone sections were fixed in paraformaldehyde (4%) for hematoxylin and eosin (H&E) staining. The TGPs were frozen in liquid nitrogen and stored at −80°C for further analysis, such as quantitative real-time polymerase chain reaction (RT-qPCR) and western blotting.
Flavone Content Measurement
The flavone was analyzed by HPLC at BioNovoGene (Suzhou, China) (http://www.bionovogene.com/). Analysis of the samples was performed using liquid chromatography-mass spectrometry (LC-MS), carried out on a Waters ACQUITY UPLC equipped with an AB 4000 Triple Quadrupole Mass Spectrometer (AB 4000). The chromatographic conditions were as follows: ACQU ITY UPLC® BEH C18 column (2.1 × 100 mm, 1.7 μm, Waters, United States), injection volume 5 μl, column temperature 40°C, mobile phase A1-0.1% formic acid water, mobile phase B- Methanol, flow rate 0.25 ml/min. The mass spectrometry conditions were as follows: Electrospray ionization (ESI) source, negative ion ionization mode. The ion source temperature was 500°C, the ion source voltage was −4,500 V, the collision gas was 6 psi, the curtain gas was 30 psi, and the atomizing gas and auxiliary gas were both 50 psi. Scans were performed using multiple reaction monitoring (MRM). The samples were quantified according to the standard curve (Supplementary Table S1).
Hematoxylin and Eosin Staining
The tibia samples from each group were fixed in 4% paraformaldehyde. The paraffin sections were prepared by conventional methods, and then dewaxed in xylene, dehydrated in graded ethanol solutions and embedded in paraffin. Tibia section (5 μm) were stained with hematoxylin and eosin (H&E), and examined by optical microscopy using a Motic BA600-4 microscope (MoticⓇ, Xiamen, China) to assess tibial histopathology.
RNA Isolation and qRT-PCR
Total RNA of the tibia sample was extracted using Trizol reagent (China Biological Engineering Co., Ltd., China). The integrity of total RNA was checked on denaturing formaldehyde agarose gel and the concentration of RNA was quantified by Nanodrop 2000 (Thermo scientific) spectrophotometer. Primers specific for GAPDH (Glyceraldehyde-3 phosphate dehydrogenase), OPG, RNAKL, BMP-2 and Runx2 were designed by Premier software (#Version 5.0., Canada) and synthesized by Tsingke Biotechnology Co., Ltd. (Zhengzhou, China). All the primer sequences in this study are shown in Table 1. GAPDH was used as a loading control. The qRT-PCR was performed according to the SYBR Green I PCR Master Mix (Vazyme Biotech Co., Ltd., China). The PCR amplification thermal cycling consisted of predeformation at 95°C for 30 s, followed by 40 cycles consisting of denaturation at 95°C for 10 s, extension at 72°C for 30 s, and annealing at 60°C for 30 s. Calculated according to the expression of mRNA using the 2−ΔΔCT method (Lin et al., 2022). All the assays were performed three times to ensure accuracy.
Western Blot Analysis
Tibial growth plates from all the groups were initially homogenized in PBS (freezing-cold) solution and kept at 4°C for 2 h. After that, the homogenate was centrifuged at 14,000 rpm for 10 min to collect the total protein. Coomassie Brilliant Blue G-250 method was used to calculate the content of total protein (Zhang et al., 2018). Moreover, protein from the samples was separated with the help of SDS polyacrylamide gel and then moved to polyvinylidene fluoride (PVDF) membranes, which were incubated for 2 h in 5% skimmed milk. Next, membranes were incubated overnight at 4°C with primary antibodies against OPG (#0431R, 1:1,000, Bioss, Beijing, China), RANKL (#0747R, 1:1,000, Bioss, Beijing, China), and β-actin (#WL01372, 1:1,000, Wanleibio, Shenyang, China). Subsequently, membranes were washed extensively using Tris-buffered saline, containing 0.1% Tween 20 (TBST) for 30 min and then incubated for 2 h at room temperature with secondary antibody (HRP goat anti-rabbit IgG, #WLA023a, 1:5,000 dilution; Wanleibio, Shenyang, China). Later, with the help of TBST, membranes were washed 30 min again rigorously, and images were taken with an imaging system (AI600, CE, United States). At the end, gray level of the exposed strip was analyzed.
Statistical Analysis
All of the statistical analyses were performed with SPSS. Software (#version 26.0, SPSS Inc., Chicago, IL, United States) by using one-way ANOVA followed by the LSD test as the post-hoc test. The correlation analysis between body weight and tibial parameters, tibia growth parameters and bone formation, and bone resorption-related indices were performed by Prism statistical analysis software (#Version 8.00, GraphPad Software, Inc.). The PCA analysis of bone resorption and bone formation-related genes through genetic cloud tools (a free online data analysis platform, https://www.omicstudio.cn). All experiments were performed at least in triplicates with data expressed as mean ± standard deviation (SD), and it was considered that p < 0.05 represents statistically significant.
Results
Effects of Total Flavonoids of Rhizoma Drynariae on Tibia Growth in Broilers With Tibial Dyschondroplasia
The tibia growth indexes of chickens are presented in Figure 1. The body weight and tibia weight of the TD group were significantly lower than those of the Con group (p < 0.001 and p < 0.001, respectively). In addition, no statistically difference in body weight and tibia weight have been observed between LTFRD and Con groups on day 7 (p = 0.172 and p = 0.857, respectively), day 14 (p = 0.556 and p = 0.941, respectively), day 21 (p = 0.828 and p = 0.668, respectively) (Figures 1B,C). The TFRD treatment group displayed a good therapeutic effect in body weight and tibia weight, significantly different from the TD group, especially in the high-dose TFRD group. The increase of tibia weight index in TD group also reflected the growth retardation of TD broilers, whereas TFRD supplementation could alleviate the increase of tibia weight index caused by TD (Figure 1D).
FIGURE 1. Effects of TFRD on tibial growth in TD broilers. (A) The tibial of chickens. (B) Body weight of chickens. (C) Tibia weight. (D) Tibia weight index. (E) Tibia length. (F) Tibia growth plate width. (G) Tibia growth plate index. (H) The pearson correlation analysis between tibia weight and body weight. (I) The pearson correlation analysis between tibia length and body weight. (J) The pearson correlation analysis between tibia growth plate width and body weight. The results are represented as mean ± SD. a, b, c, and d represent significant differences between groups (p < 0.05).
As shown in Figures 1E–G, the tibia length was significantly decreased and the TGP width was markedly increased in the TD group compared with the Con group on day 7 (p = 0.042 and p = 0.004, respectively), day 14 (p = 0.008 and p < 0.001, respectively) and day 21 (p < 0.001 and p < 0.001, respectively). TFRD treatment can ameliorate TD-induced abnormality of these two indicators and has the best effect on day 21. The higher the TGP index for TD broilers, the more severe the tibial damage. In the present study, the TGP index in the TD group was higher than that in the Con group, indicating that the tibia in TD broilers was seriously damaged. Conversely, TFRD supplementation can effectively reverse the increase in the TGP index. These results showed that TFRD supplementation alleviated tibia length, TGP width, and TGP index in TD broilers.
Then we conducted a correlation analysis between tibia-related parameters and bodyweight of broilers. As shown in Figure 1, the bodyweight of broilers was significantly positively correlated with tibia length (p < 0.01, r = 0.997) and tibia weight (p < 0.01, r = 0.966). The body weight was weakly correlated with TGP width (p = 0.038, r = 0.804) (Figures 1H–J). The above results indicated that TFRD could promote the tibia growth of TD broilers.
Effects of Total Flavonoids of Rhizoma Drynariae on Morphology and Histology of the Tibia in Broilers With Tibial Dyschondroplasia
The results of tibia morphology and histology on the tibia of broilers revealed that the TGP width of LTFRD group was not abnormal, and the arrangement of bone trabecula was normal as compared to the Con group (Figures 2A,B). Additionally, the growth plate of the tibia in the TD group became thicker, and the structure of the trabecular bone was destroyed and broken as compared to the Con group. Especially on day 21, the growth plate width of tibia in the TD group became abnormally thick. After treatment with different doses of TFRD, it was found that the growth plate width of TD + LTFRD group, TD + MTFRD group and TD + HTFRD group is gradually restored to normal when compared with TD group, and the arrangement of bone trabeculae is also becoming orderly, especially on day 21. These results showed that TFRD could restore TD injury of broiler chickens by improving the structure of the TGP.
FIGURE 2. Effects of TFRD on growth plate of tibia bone in TD broilers. (A) The morphology result of tibial growth plate. GP, growth plate; TDL, tibial dyschondroplasias lesion. (B) The HE staining observation result of tibial growth plate. Scar bar = 1,000 µm.
Effects of Total Flavonoids of Rhizoma Drynariae on Vascular Distribution of the Tibia Growth Plate in Broilers With Tibial Dyschondroplasia
Bone formation is dependent upon the richly vascular system. In the present study, the blood vessels in the TGP area of the LTFRD group were more abundant than those in the Con group (Figure 3). It is noted that there is no invasion of blood vessels in the corresponding area of the tibia in the TD group, suggesting that the growth of the tibia in TD broilers may be inhibited (Figure 3). Compared with the TD group, after 2 weeks of intervention with TFRD, the blood vessels in the TGP area of broilers in the TD + LTFRD, TD + MTFRD and TD + HTFRD groups were obviously increased (Figure 3). Collectively, the above findings indicated that TFRD treatment could enhance vascular infiltration in TGP area in TD broilers to a certain extent.
FIGURE 3. Effects of TFRD on vascular invasion of tibial growth plate in TD broilers. HE staining was used to observe the vascular infiltration of growth plates in each group. The black arrow points to the blood vessel. Scar bar = 100 µm.
Effects of Total Flavonoids of Rhizoma Drynariae on Gene Expression Levels of BMP-2/Runx2 and OPG/RANKL Pathways of TGP in Broilers With Tibial Dyschondroplasia
PCA analysis was performed as an unsupervised pattern recognition method to explore the expression of bone-related genes in TD broilers and distinguish the potential effects of different doses of TFRD on bone-related genes (Figures 4A–D). PCA scores of BMP-2, Runx2, OPG, and RANKL genes showed that the distance between the TD and Con groups was greater than that between LTFRD and Con groups. Moreover, there was a dose-dependent effect in the TFRD-treated chicken groups. The TD + HTFRD group was farther away from the TD group, which suggests that supplementation of 500 mg/kg TFRD in drinking water has the best protective effect on protecting the bone injury of TD broilers.
FIGURE 4. Effects of TFRD on gene expression levels of BMP-2/Runx2 and OPG/RANKL in tibial growth plate in TD broilers. (A–D) Principal component analysis (PCA) of the BMP-2, Runx2, OPG and RANKL gene expressions with an unsupervised pattern recognition method. (E) The mRNA expression of genes BMP-2, Runx2, OPG, RANKL and the ratio of OPG/RANKL at day 7, 14 and 21. (F) The heat map shows the mRNA levels of the bone-related genes. The results are represented as mean ± SD. a, b, c, and d represent significant differences between groups (p < 0.05).
The expression levels of BMP-2/Runx2 and OPG/RANKL genes in TGP are displayed in Figure 4E. Compared with the Con group, the mRNA expression levels of Runx2 on day 7 (p = 0.002) and day 21 (p = 0.020) were dramatically decreased in the TD group. Similarly, the expression levels of BMP-2, OPG and OPG/RANKL ratio were significantly decreased in the TD group during the experiment. Conversely, the mRNA expression of RANKL gene was significantly up-regulated throughout the experimental period in the TD group compared to the Con group (p < 0.001). Compared with the TD group, the three different doses of TFRD treatment clearly up-regulated the expression of BMP-2 on day 14 (p = 0.710, p = 0.108, and p = 0.001, respectively) and day 21 (p = 0.357, p = 0.069, and p = 0.032, respectively), the expression of Runx2 on day 14 (p = 0.466, p = 0.008, and p < 0.001, respectively) and day 21 (p = 0.017, p < 0.001, and p < 0.001, respectively), the expression of OPG on day 14 (p = 0.026, p = 0.010, and p = 0.002, respectively) and day 21 (p < 0.001, p < 0.001, and p < 0.001, respectively), and the ratio of OPG/RANKL on day 14 (p < 0.001, p < 0.001, and p = 0.001, respectively) and day 21 (p < 0.001, p < 0.001, and p < 0.001, respectively; Figure 4E). For the expression level of the RANKL gene, the three different doses of TFRD treatment could significantly reduce its expression on the TGP on day 14 and day 21 when compared with the TD group (p < 0.001, Figure 4E).
A heat map in Figure 4F further presented the expression levels of osteogenesis-related genes. The expression of bone resorption gene RANKL in TD broilers (red) was higher than that of Con group (blue), especially on days 14 and 21. Additionally, the expression of bone formation-related genes OPG, BMP-2, Runx2 and OPG/RANKL ratio in TD broilers was evidently decreased. A high dose of TFRD had a better effect on restoring the expression of osteogenesis-related genes. Collectively, these results indicated that TFRD could ameliorate the expression levels of osteogenesis-related genes on the TGP of broiler chickens with TD.
Effects of Total Flavonoids of Rhizoma Drynariae on the Protein Expression Levels of OPG/RANKL of Tibia Growth Plate in Broilers With Tibial Dyschondroplasia
The OPG and RANKL protein levels were assessed by western blotting analysis in the tibia of broiler chickens (Figure 5). The results demonstrated that the expression levels of OPG protein and OPG/RANKL ratio were significantly upregulated in the LTFRD group as compared to the Con group on day 7 (p < 0.001 and p = 0.129, respectively), day 14 (p < 0.001 and p < 0.001, respectively), and day 21 (p = 0.013 and p = 0.016, respectively). Compared with the Con group, the OPG protein and OPG/RANKL ratio were significantly down-regulated in the TD group on day 7 (p < 0.001 and p = 0.26, respectively), day 14 (p < 0.001 and p = 0.001, respectively), and day 21 (p = 0.005 and p = 0.008, respectively; Figures 5A,B,D). Moreover, it is found that TD + HTFRD group apparently increased the expression levels of OPG protein and the ratio of OPG/RANKL on days 14 (p = 0.002, and p = 0.008, respectively) and 21 (p = 0.007, and p = 0.011, respectively) as compared to TD group (Figures 5A,B). Additionally, the RANKL protein level was significantly up-regulated in the TGP of TD chicken during the experiment as compared to the Con group (p < 0.001). A high dose of TFRD treatment decreased the protein level of RANKL significantly as compared to TD group on day 14 and day 21 (p < 0.001 and p < 0.001, respectively). In addition, there was no difference in the expression of RANKL protein between the LTFRD and Con groups on day 14 and day 21 (p = 0.124 and p = 0.761, respectively, Figure 5C). These findings indicated that TFRD could activate the OPG/RANKL signaling pathway in broiler chickens with TD.
FIGURE 5. Effects of LTFRD on OPG/RANKL protein level of tibial growth plate in TD broilers. (A) The gray scale analysis of OPG, RANKL and β-ACTIN. (B) The protein level of OPG. (C) The protein level of RANKL. (D) The protein level of OPG/RANKL The results are represented as mean ± SD. a, b, c, and d represent significant differences between groups (p < 0.05).
Correlation Analysis Between Tibia Growth Parameters and Bone Formation-Related Indices in Broilers
To further analyze the potential relationship between osteogenesis-related biological indicators and tibial growth parameters, Pearson correlation analysis was performed (Figure 6). The results showed that body weight, tibia weight and tibia length were positively correlated with osteogenesis-related genes OPG, BMP-2 and Runx2, but not significant; they were negatively correlated with RANKL protein (r = −0.291, −0.302, and −0.296; p = 0.052, 0.044, and 0.048, respectively). Interestingly, the TGP width and its index were negatively correlated with the gene expression of BMP-2 (r = −0.582 and r = −0.634; p < 0.001 and p < 0.00, respectively) and OPG (r = −0.525 and r = −0.600; p < 0.001 and p < 0.001, respectively). Moreover, the protein expression of OPG and the protein expression of OPG/RANKL ratio were also negatively correlated with the TGP width (r = −0.669 and r = −0.627; p < 0.001 and p < 0.001, respectively) and its index (r = −0.669 and r = −0.675; p < 0.001 and p < 0.001, respectively). These results indicate that BMP-2/Runx2 and OPG/RANKL pathways may be involved in the regulation of tibia formation in TD broilers by affecting the function of tibia growth plate.
FIGURE 6. Correlation heat map analysis of tibial growth parameters and osteogenesis-related biological indicators in broilers. Note: The different rectangles are colored based on the Pearson correlation coefficients between tibia growth performance expression levels and bone formation-related indices. The intensity of color represents the degree of correlation, red represents positive correlation, blue represents negative correlation. *indicates significant correlation (p < 0.05), ***indicates highly significant correlation (p < 0.001).
Discussion
Previous studies have indicated that TFRD, as a Rhizoma Drynariae extract, has been widely used in the prevention and treatment of many bone disorders including femoral head necrosis (Lv et al., 2021), osteoporosis (Wang et al., 2011), large tibial defects (Sun et al., 2021), by promoting osteogenesis, increasing bone mineral density, and preventing fracture in animal experiments (Shen et al., 2020). In vitro experiments have also shown that TFRD can increase the activity of osteoblasts, enhance the osteogenic effect of the membrane and improve the resistance of osteoblasts (Yu et al., 2021). In fast-growing poultry, tibial dyschondroplasia is an intractable tibiotarsal bone disorder that affects the proximal growth plate of the tibia bone (Huang et al., 2019a). In the present study, our results revealed that TFRD supplementation reduces lameness in TD broilers by increasing tibia growth parameters, reducing growth plate damage, improving growth plate vascular distribution, and regulating the expression of OPG/RANKL ratio, suggesting that TFRD had the effect on preventing and treating broiler chickens with TD (Figure 7).
FIGURE 7. The schematic diagram of the protective mechanism of TFRD on TD broilers by regulating the OPG/RANKL expression.
Tibia growth performances have been widely used to evaluate the overall bone growth and development of broiler. The reduction of tibia weight and tibia indexes length indicates that the development of broiler bone is impeded. In the current study, the body weight gain of broilers in the TD group was significantly lower than those in the Con group. Additionally, the tibia weight and length were significantly decreased, while the tibia growth plate was markedly wider in TD broilers. TFRD as a TCM extract with bone strength (Zhang et al., 2017) confirm increased tibial volume fraction and accelerated fracture healing in fractured mice (Sun et al., 2021). In the present study, we demonstrated similar results that TFRD can increase body weight, tibia weight and tibia length, and decrease growth plate width, indicating TFRD plays a role in bone strengthening in TD broilers. The previous finding of Yao and coworkers (Yao et al., 2018) are consistent with the results of this study: TFRD also reduced mortality and accelerated recovery from lameness in TD broilers. These findings suggested that TFRD has a significant beneficial effect on the tibia growth performance of TD broilers.
The most direct reflection of the low tibia growth performance of the broiler is the transparent uncalcified embolism in the growth plate. In the present study, we investigated irregular transparent cartilage emboli with thickening in TD broilers (Figure 2A). These results are consistent with previous studies on thiram damage to growth plates of broilers (Zhang et al., 2018). The occurrence of transparent uncalcified embolism is mainly due to decreased angiogenesis in the tibial cartilage region, resulting in insufficient nutritional supply of osteogenic related cells, blocked metabolism and accumulation in the growth plate, resulting in calcification termination. TFRD has been shown to repair the properties of chondrocytes and promote the metabolism of osteogenic related cells, raise the proliferation of osteoblasts while inhibiting the activity of osteoclasts, thereby enhancing bone mass (Wei et al., 2017; Yao et al., 2018).
The TFRD was previously found to promote osteoblast activity, improve the accumulation of growth plate calcification and reduce growth plate width in TD broilers (Yao et al., 2018). In this experiment, the tibia growth plate of broilers in the LTFRD group showed good growth characteristics, and there was no difference in GP area between the LTFRD group and the normal group. Furthermore, TD lesion area of the growth plate of TD broilers treated with TFRD was significantly reduced, especially on day 21, and TD + HTFRD group had the best effect. These findings indicated that TFRD could be used as a long-term additive and has a good effect on improving osteoblast metabolism, reducing growth plate calcification in TD broilers.
To further study the occurrence and development of TD growth plate, we observed the vessels of the growth plate. There was almost no vascular infiltration in the growth plate in TD group, indicating that the growth plate of TD broilers was seriously damaged. Restoration of bone blood and nutrient supply is very important for tibial growth plate restoration of TD broilers. TCM can promote vascular regeneration and reconstruction (Huang et al., 2021b; Lv et al., 2021). TFRD, as a TCM widely used in bone diseases, has been proven that have the potential to promote vascular remodeling. Lv et al. (2021). treated avascular necrosis of the femoral head mice with TFRD and found that TFRD could improve the supply of femoral blood and nutrition, promote blood vessel regeneration, and strengthen the recovery of femoral osteonecrosis. Besides, TFRD stimulated vascular endothelial cell production and alleviated decreased vascular abundance in fractured mice (Shen et al., 2020). The previous results findings are similar to our study on TD broilers blood vessels, in the current study, vascular infiltration was enriched in Con and LTFRD groups, the vascular enrichment of TD broilers was gradually improved after TFRD treatment, and the effect was most obvious on day 21, which was the same as that of Con group. Therefore, TFRD can ameliorate the impaired growth plate development of TD broilers by promoting vascular growth.
The obstructed metabolism of osteoblasts and osteoclasts and accumulation in growth plate are the direct causes of transparent uncalcified embolism in TD broilers (Huang et al., 2019b). Promoting the maturation of osteoblasts and inhibiting differentiation of osteoclasts are important factors in alleviating TD broilers. BMP-2/Runx2 is a key signaling pathway that regulates osteoblast activity and plays a decisive role in bone formation. BMP-2 has strong osteoinductive potential in vivo and in vitro. It is an important stimulus for bone development (Soysa and Alles, 2016).
The Runx2 is a transcription factor important for the maturation of chondrocytes. It is also an important gene for osteoblast differentiation. TFRD has been proven to promote osteoblast activity and treat osteoporosis, bone fragility fractures, and bone loss (Lian and Stein., 2003; Wei et al., 2017). Previous studies have shown that TFRD can treat osteoporosis in rabbits and rats through BMP-2 related signaling pathway (Wong and Rabie, 2006; Wei et al., 2017). In TD broilers, TFRD improves the growth performance and restore normal vigor and significantly increases the expression of BMP-2 and Runx2 genes (Yao et al., 2018). We found that TFRD ameliorates TD-induced down-regulation of BMP-2 and Runx2 genes, suggesting that TFRD could enhance osteoblast activity and promote its expression during bone development in TD broilers. In addition, research has shown that TFRD can promote bone formation by activating the BMP-2/Smad signaling pathway and significantly improving the expression level of BMP-2/Runx2 in tibial vascular injury rats (Sun et al., 2021). These findings suggest that BMP-2 and Runx2 are important genes controlling tibial achondroplasia, and TFRD promotes osteoblast activity by regulating the expression of BMP-2 and Runx2, thereby ameliorating bone injury in TD broilers.
The OPG/RANKL is an important bone resorption regulator that inhibits osteoclast formation and activation. The OPG binds to the RANKL to form blocking competition for its receptor RANK, which is close to the differentiation and fusion of osteoclast precursors induced by osteoblasts. It regulates osteoclast formation, differentiation, maturation and apoptosis (Piemontese et al., 2016). Therefore, we detected the expression of OPG/RANKL relative content in TD broilers. The results exhibited that OPG expression was significantly down-regulated and RANKL expression was significantly upregulated in TD broilers compared with the Con group. This indicated that TD broiler osteoclast differentiation was strong, and tibia was seriously damaged. Several studies have shown that TFRD induces potential anti-osteoporosis activity by regulating targets in bone metabolic signaling pathways OPG/RANKL, and prevents and treats osteoporosis by inhibiting bone resorption or stimulating bone (Ibáñez, et al., 2019; Amin, et al., 2020; Ma, et al., 2020). Similar results were obtained in our study with TD broilers: TFRD treatment increased OPG expression level and decreased RANKL expression level of TD broilers, thereby improving tibial dysplasia in TD broilers. The high-dose group of TFRD had the best treatment effect on day 21. The ratio of OPG/RANKL can reflect the recovery rate of damaged bone. Our study found that TFRD supplementation significantly increased the OPG/RANKL ratio of TD broilers. Taken together, these results indicate that TFRD can inhibit the activity of osteoclasts by regulating the expression of OPG and RANKL, thus promoting bone development.
Conclusion
The TFRD has a potential therapeutic effect on improving growth performance and promoting the vascular distribution of the tibial growth plate in TD broilers. Importantly, TFRD regulates bone formation via OPG/RANKL expression and thus diminishes lameness and recovers TD in broilers. Which provide a valuable reference for further exploring the effects of TFRD on poultry leg disease. And TFRD can be used as a potential target medicine to treat thiram-induced tibial dyschondroplasia.
Data Availability Statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Ethics Statement
All broiler chicks experiments were followed the guidelines of institutional Animal Welfare and Research Ethics Committee of the College of Animal Science and Veterinary Medicine, Henan Agricultural University Zhengzhou, China (Permit No: 17-0126).
Author Contributions
TX and JZ analyzed the data and wrote the paper. SH, AS, and MA contributed to the revision and editing of the manuscript. CZ, LXL, and LL performed the experiments. WJ and QC provided technical and/or conceptual assistance. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the China Postdoctoral Science Foundation (No. 2020M672234), the Outstanding Talents of Henan Agricultural University (No.30500421), and Key Scientific Research Project of Henan Higher Education Institutions of China (No. 21A230013).
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.
The reviewer AL declared a shared affiliation with the author AS to the handling editor at the time of review.
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/fphar.2022.881057/full#supplementary-material
Abbreviations
BMP-2, bone morphogenetic protein-2; GP, growth plate; OPG, osteoprotegerin; RANKL, receptor activator of nuclear factor kappa-B ligand; Runx2, runt-related transcription factor-2; TCM, traditional Chinese medicine; TD, tibial dyschondroplasia; TDL, tibial dyschondroplasias lesion; TFRD, total flavonoids of Rhizoma Drynariae; TGP, tibia growth plate.
References
Amin, N., Boccardi, V., Taghizadeh, M., and Jafarnejad, S. (2020). Probiotics and Bone Disorders: the Role of RANKL/RANK/OPG Pathway. Aging Clin. Exp. Res. 32 (3), 363–371. doi:10.1007/s40520-019-01223-5
Bruderer, M., Richards, R. G., Alini, M., and Stoddart, M. J. (2014). Role and Regulation of RUNX2 in Osteogenesis. Eur. Cell. Mater 28, 269–286. doi:10.22203/ecm.v028a19
Chen, G. Y., Chen, J. Q., Liu, X. Y., Xu, Y., Luo, J., Wang, Y. F., et al. (2021). Total Flavonoids of Rhizoma Drynariae Restore the MMP/TIMP Balance in Models of Osteoarthritis by Inhibiting the Activation of the NF-Κb and PI3K/AKT Pathways. Evid. Based Complement. Altern. Med. 2021, 6634837. doi:10.1155/2021/6634837
Dong, P.-F., Jin, C., Lian, C.-Y., Wang, L., and Wang, Z.-Y. (2021). Enhanced Extracellular Matrix Degradation in Growth Plate Contributes to Manganese Deficiency-Induced Tibial Dyschondroplasia in Broiler Chicks. Biol. Trace Elem. Res 2021, 1–10. doi:10.1007/s12011-021-02921-w
Feng, X. H., and Derynck, R. (2005). Specificity and Versatility in Tgf-Beta Signaling through Smads. Annu. Rev. Cell. Dev. Biol. 21, 659–693. doi:10.1146/annurev.cellbio.21.022404.142018
Garner, J. P., Kiess, A. S., Mench, J. A., Newberry, R. C., and Hester, P. Y. (2012). The Effect of Cage and House Design on Egg Production and Egg Weight of White Leghorn Hens: an Epidemiological Study. Poult. Sci. 91 (7), 1522–1535. doi:10.3382/ps.2011-01969
Huang, N. C., Dai, L. G., Kang, L. Y., Huang, N. C., Fu, K. Y., Hsieh, P. S., et al. (2021b). Beneficial Effects of Astragaloside IV-Treated and 3-Dimensional-Cultured Endothelial Progenitor Cells on Angiogenesis and Wound Healing. Ann. Plast. Surg. 86 (2S Suppl. 1), S3–S12. doi:10.1097/SAP.0000000000002655
Huang, S., Kong, A., Cao, Q., Tong, Z., and Wang, X. (2019b). The Role of Blood Vessels in Broiler Chickens with Tibial Dyschondroplasia. Poult. Sci. 98 (12), 6527–6532. doi:10.3382/ps/pez497
Huang, S. C., Cao, Q. Q., Cao, Y. B., Yang, Y. R., Xu, T. T., Yue, K., et al. (2021a). Morinda Officinalis Polysaccharides Improve Meat Quality by Reducing Oxidative Damage in Chickens Suffering from Tibial Dyschondroplasia. Food Chem. 344, 128688. doi:10.1016/j.foodchem.2020.128688
Huang, S. C., Li, L., Rehman, M. U., Gao, J. D., Zhang, L. H., Tong, X. L., et al. (2019a). Tibial Growth Plate Vascularization Is Inhibited by the Dithiocarbamate Pesticide Thiram in Chickens: Potential Relationship to Peripheral Platelet Counts Alteration. Environ. Sci. Pollut. Res. Int. 26 (36), 36322–36332. doi:10.1007/s11356-019-06664-z
Huang, S. C., Rehman, M. U., Lan, Y. F., Qiu, G., Zhang, H., Iqbal, M. K., et al. (2017). Tibial Dyschondroplasia Is Highly Associated with Suppression of Tibial Angiogenesis through Regulating the HIF-1α/VEGF/VEGFR Signaling Pathway in Chickens. Sci. Rep. 7 (1), 9089. doi:10.1038/s41598-017-09664-6
Huang, S. C., Zhang, L. H., Zhang, J. L., Rehman, M. U., Tong, X. L., Qiu, G., et al. (2018). Role and Regulation of Growth Plate Vascularization during Coupling with Osteogenesis in Tibial Dyschondroplasia of Chickens. Sci. Rep. 8 (1), 3680. doi:10.1038/s41598-018-22109-y
Ibáñez, L., Rouleau, M., Wakkach, A., and Blin-Wakkach, C. (2019). Gut Microbiome and Bone. Jt. Bone Spine 86 (1), 43–47. doi:10.1016/j.jbspin.2018.02.008
Kapell, D. N. R. G., Hocking, P. M., Glover, P. K., Kremer, V. D., and Avendaño, S. (2017). Genetic Basis of Leg Health and its Relationship with Body Weight in Purebred turkey Lines. Poult. Sci. 96 (6), 1553–1562. doi:10.3382/ps/pew479
Kierończyk, B., Rawski, M., Józefiak, D., and Świątkiewicz, S. (2017). Infectious and Non-infectious Factors Associated with Leg Disorders in Poultry - a Review. Ann. Anim. Sci. 17, 645–669. doi:10.1515/aoas-2016-0098
Ledwaba, M. F., and Roberson, K. D. (2003). Effectiveness of Twenty-Five-Hydroxycholecalciferol in the Prevention of Tibial Dyschondroplasia in Ross Cockerels Depends on Dietary Calcium Level. Poult. Sci. 82 (11), 1769–1777. doi:10.1093/ps/82.11.1769
Lian, J. B., and Stein, G. S. (2003). Runx2/Cbfa1: a Multifunctional Regulator of Bone Formation. Curr. Pharm. Des. 9 (32), 2677–2685. doi:10.2174/1381612033453659
Lin, L., Fu, P., Chen, N., Gao, N., Cao, Q., Yue, K., et al. (2022). Total Flavonoids of Rhizoma Drynariae Protect Hepatocytes against Aflatoxin B1-Induced Oxidative Stress and Apoptosis in Broiler Chickens. Ecotoxicol. Environ. Saf. 230, 113148. doi:10.1016/j.ecoenv.2021.113148
Lui, J. C., Nilsson, O., and Baron, J. (2014). Recent Research on the Growth Plate: Recent Insights into the Regulation of the Growth Plate. J. Mol. Endocrinol. 53 (1), T1–T9. doi:10.1530/JME-14-0022
Lv, W., Yu, M., Yang, Q., Kong, P., and Yan, B. (2021). Total Flavonoids of Rhizoma Drynariae Ameliorate Steroid-induced A-vascular N-ecrosis of the F-emoral H-ead via the PI3K/AKT P-athway. Mol. Med. Rep. 23 (5), 345. doi:10.3892/mmr.2021.11984
Ma, Z., Li, S., and Sun, Y. (2020). Effect of Enhanced Masticatory Force on OPG, RANKL and MGF in Alveolar Bone of Ovariectomized Rats. J. Appl. Oral Sci. 28, e20190409. doi:10.1590/1678-7757-2019-0409
Nuschke, A., Rodrigues, M., Stolz, D. B., Chu, C. T., Griffith, L., and Wells, A. (2014). Human Mesenchymal Stem Cells/multipotent Stromal Cells Consume Accumulated Autophagosomes Early in Differentiation. Stem Cell. Res. Ther. 5 (6), 140. doi:10.1186/scrt530
Piemontese, M., Xiong, J., Fujiwara, Y., Thostenson, J. D., and O'Brien, C. A. (2016). Cortical Bone Loss Caused by Glucocorticoid Excess Requires RANKL Production by Osteocytes and Is Associated with Reduced OPG Expression in Mice. Am. J. Physiol. Endocrinol. Metab. 311 (3), E587–E593. doi:10.1152/ajpendo.00219.2016
Pierrefite-Carle, V., Santucci-Darmanin, S., Breuil, V., Camuzard, O., and Carle, G. F. (2015). Autophagy in Bone: Self-Eating to Stay in Balance. Ageing Res. Rev. 24 (Pt B), 206–217. doi:10.1016/j.arr.2015.08.004
Rath, N. C., Huff, W. E., and Huff, G. R. (2007). Thiram-induced Changes in the Expression of Genes Relating to Vascularization and Tibial Dyschondroplasia. Poult. Sci. 86, 2390–2395. doi:10.3382/ps.2007-00219
Shen, Z., Chen, Z., Li, Z., Zhang, Y., Jiang, T., Lin, H., et al. (2020). Total Flavonoids of Rhizoma Drynariae Enhances Angiogenic-Osteogenic Coupling during Distraction Osteogenesis by Promoting Type H Vessel Formation through PDGF-BB/PDGFR-β Instead of HIF-1α/ VEGF Axis. Front. Pharmacol. 11, 503524. doi:10.3389/fphar.2020.503524
Song, C., Song, C., and Tong, F. (2014). Autophagy Induction Is a Survival Response against Oxidative Stress in Bone Marrow-Derived Mesenchymal Stromal Cells. Cytotherapy 16 (10), 1361–1370. doi:10.3389/fphar.2020.50352410.1016/j.jcyt.2014.04.006
Song, S. H., Zhai, Y. K., Li, C. Q., Yu, Q., Lu, Y., Zhang, Y., et al. (2016). Effects of Total Flavonoids from Drynariae Rhizoma Prevent Bone Loss In Vivo and In Vitro. Bone Rep. 5, 262–273. doi:10.1016/j.bonr.2016.09.001
Soysa, N. S., and Alles, N. (2016). Osteoclast Function and Bone-Resorbing Activity: An Overview. Biochem. Biophys. Res. Commun. 476 (3), 115–120. doi:10.1016/j.bbrc.2016.05.019
Sun, W., Li, M., Zhang, Y., Huang, Y., Zhan, Q., Ren, Y., et al. (2021). Total Flavonoids of Rhizoma Drynariae Ameliorates Bone Formation and Mineralization in BMP-Smad Signaling Pathway Induced Large Tibial Defect Rats. Biomed. Pharmacother. 138, 111480. doi:10.1016/j.biopha.2021.111480
Sun, Y., Huang, Y., Chen, T., Li, X., Chen, J., Wang, Z., et al. (2020). Effect of Downregulation of Serum MMP-3 Levels by Traditional Chinese Medicine Ingredients Combined with Methotrexate on the Progression of Bone Injury in Patients with Rheumatoid Arthritis: A Protocol for a Systematic Review and Meta-Analysis. Med. Baltim. 99 (43), e22841. doi:10.1097/MD.0000000000022841
Udagawa, N., Koide, M., Nakamura, M., Nakamichi, Y., Yamashita, T., Uehara, S., et al. (2021). Osteoclast Differentiation by RANKL and OPG Signaling Pathways. J. Bone Min. Metab. 39 (1), 19–26. doi:10.1007/s00774-020-01162-6
Wang, X., Zhen, L., Zhang, G., Wong, M. S., Qin, L., and Yao, X. (2011). Osteogenic Effects of Flavonoid Aglycones from an Osteoprotective Fraction of Drynaria Fortunei-Aan In Vitro Efficacy Study. Phytomedicine 18 (10), 868–872. doi:10.1016/j.phymed.2011.01.022
Wei, X., Xu, A., Shen, H., and Xie, Y. (2017). Qianggu Capsule for the Treatment of Primary Osteoporosis: Evidence from a Chinese Patent Medicine. BMC Complement. Altern. Med. 17 (1), 108. doi:10.1186/s12906-017-1617-3
Wong, R. W., and Rabie, A. B. (2006). Systemic Effect of Crude Extract from Rhizome of Drynaria Fortunei on Bone Formation in Mice. Phytother. Res. 20 (4), 313–315. doi:10.1002/ptr.1842
Yao, W., Zhang, H., Jiang, X., Mehmood, K., Iqbal, M., Li, A., et al. (2018). Effect of Total Flavonoids of Rhizoma Drynariae on Tibial Dyschondroplasia by Regulating BMP-2 and Runx2 Expression in Chickens. Front. Pharmacol. 9, 1251. doi:10.3389/fphar.2018.01251
Yu, F., Lv, Q., Tong, Z., Song, W., Zhao, Z., Yuan, B., et al. (2021). Effect of Total Flavonoids from Drynaria Rhizome on Bone Loss in Ovariectomized Rats. Trop. J. Pharm. Res. 18 (6), 1285–1289. doi:10.4314/tjpr.v18i6.20
Zhang, H., Mehmood, K., Li, K., Rehman, M. U., Jiang, X., Huang, S., et al. (2018). Icariin Ameliorate Thiram-Induced Tibial Dyschondroplasia via Regulation of WNT4 and VEGF Expression in Broiler Chickens. Front. Pharmacol. 9, 123. doi:10.3389/fphar.2018.00123
Zhang, N. D., Han, T., Huang, B. K., Rahman, K., Jiang, Y. P., Xu, H. T., et al. (2016). Traditional Chinese Medicine Formulas for the Treatment of Osteoporosis: Implication for Antiosteoporotic Drug Discovery. J. Ethnopharmacol. 189, 61–80. doi:10.1016/j.jep.2016.05.025
Keywords: bone development, Chinese herbal medicine, leg disease, tibial dyschondroplasia, total flavonoids of Rhizoma Drynariae
Citation: Xu T, Zheng J, Jin W, Li L, Lin L, Shaukat A, Zhang C, Cao Q, Ashraf M and Huang S (2022) Total Flavonoids of Rhizoma Drynariae Ameliorate Bone Growth in Experimentally Induced Tibial Dyschondroplasia in Chickens via Regulation of OPG/RANKL Axis. Front. Pharmacol. 13:881057. doi: 10.3389/fphar.2022.881057
Received: 22 February 2022; Accepted: 16 May 2022;
Published: 26 May 2022.
Edited by:
Karl Tsim, Hong Kong University of Science and Technology, Hong Kong SAR, ChinaReviewed by:
Aoyun Li, Huazhong Agricultural University, ChinaFazul Nabi, Lasbela University of Agriculture, Water and Marine Sciences, Pakistan
Kun Li, Nanjing Agricultural University, China
Copyright © 2022 Xu, Zheng, Jin, Li, Lin, Shaukat, Zhang, Cao, Ashraf and Huang. 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: Shucheng Huang, aHVhbmcuc2NAaGVuYXUuZWR1LmNu
†These authors have contributed equally to this work