QiangGuYin Modulates the OPG/RANKL/RANK Pathway by Increasing Secretin Levels during Treatment of Primary Type I Osteoporosis

QiangGuYin (QGY) is a common Traditional Chinese medicine prescription for the treatment of osteoporosis. Previous clinical studies have found that QGY effectively improves bone mineral density (BMD) in postmenopausal women, but its underlying mechanism remains unclear. The osteoprotegerin (OPG)/receptor activator of nuclear factor kappa B ligand (RANKL)/receptor activator of nuclear factor kappa B (RANK) pathway is a classic pathway involved in osteoporosis. Secretin levels are a serum marker of osteoporosis, but their effect on the OPG/RANKL/RANK pathway has not been reported. Hence, we investigated the relationship between the OPG/RANKL/RANK pathway and secretin and further revealed the mechanism underlying the effect of QGY in the treatment of osteoporosis. Mice were divided into secretin knockdown, secretin overexpression, and corresponding control groups. Micro-computed tomography was used to detect BMD in different groups, and the results show that QGY significantly improved BMD in mice of the secretin knockdown group. To further verify this, the serum levels of OPG, RANKL, RANK, and secretin were measured by enzyme-linked immunosorbent assays, and femur levels of OPG, RANKL, RANK, and secretin were evaluated by real-time quantitative PCR and western blotting. The results show that the expression of OPG was inhibited and that of RANKL and RANK was increased in mice from the secretin knockdown group, whereas the expression of OPG was upregulated and that of RANKL and RANK was downregulated after QGY intervention. Therefore, QGY inhibited bone resorption by promoting the expression of secretin and modulating the OPG/RANKL/RANK pathway. In addition to the effect of QGY, we also revealed the general regulatory effect of secretin on the OPG/RANKL/RANK pathway. We conclude that QGY modulates the OPG/RANKL/RANK pathway by increasing secretin levels during treatment of primary type I osteoporosis. This work provides a theoretical basis for the clinical use of QGY in the treatment of osteoporosis.


Introduction
Osteoporosis has become a public health issue worldwide and is considered a major health problem, second only to coronary heart disease, by the World Health Organization [1]. Osteoporosis is a skeletal disease characterized by an imbalance between bone resorption and bone formation, which leads to low bone mineral density (BMD) and microarchitectural deterioration of bone tissue, accompanied by a consequently higher bone fragility and susceptibility to fracture [2][3][4]. Furthermore, it was estimated that the prevalence of osteoporosis in the Chinese population (greater than 50 years of age) will increase from 5.04% (1990) to 7.46% (2050) in men and from 26.28% (1990) to 39.19% (2050) in women [5]. e drugs used to treat osteoporosis include vitamin D analogs, calcitonin, bisphosphonates, estrogen, selective estrogen receptor modulators, and parathyroid hormone [6,7]. In general, drug intervention for osteoporosis may be divided into two categories, bone anabolic agents and antibone resorption agents. Unfortunately, there are several side effects, such as fever, vomiting, atypical femoral fractures, and jaw necrosis, associated with these drugs, thereby compromising quality of life [8]. Traditional Chinese medicine (TCM) has a long history in the treatment of osteoporosis [9]. Developing a new antiosteoporosis TCM prescription might provide a novel strategy for the treatment of osteoporosis.
It is well known that osteoporosis is closely associated with diseases of the digestive system [10][11][12][13]. Secretin, a hormone, is secreted by duodenal S cells [14] and plays an important role in the pathological processes of digestive system diseases and osteoporosis. Previous studies have revealed that secretin is a serum protein marker of type I osteoporosis in postmenopausal women [15,16]. Structurally, secretin is closely associated with calcitonin, a hormone that inhibits osteoclast function [10,13,17]. However, secretin was also known to decrease gastric H + secretion in the duodenum [15], thus inhibiting calcium absorption, resulting in low bone mass [18]. However, the mechanism of action of secretin in osteoporosis remains unclear. Meanwhile, it has been reported that the osteoprotegerin (OPG)/ receptor activator of nuclear factor kappa B ligand (RANKL)/receptor activator of nuclear factor kappa B (RANK) signaling pathway is a major regulator of osteoporosis [4,19]. erefore, we speculated that secretin might play a crucial role in regulating the OPG/RANKL/RANK signaling pathway in bone remodeling in primary type I osteoporosis.
QiangGuYin (QGY) is a TCM prescription that has been widely used in the prevention and treatment of osteoporosis clinically, and its effectiveness and safety have been previously established [1,3]. In TCM theory, QGY has been reported to possess tonifying "Qi" and warming meridians. Moreover, it was reported that QGY relieves clinical symptoms of pain, deformation of the spine, and fractures by improving BMD and reducing fracture risk. QGY consists of 10 [20]. Hedysarum multijugum Maxim. [21] increases bone mineral density (BMD) content by improving the level of OPG and reducing the level of RANKL in ovariectomy-induced mice. Drynariae Rhizoma enhances the strength and hardness of bone [22]. Meanwhile, our previous studies reported that Drynaria flavonoids improve osteoporosis by promoting the expression of OPG and inhibiting RANKL and RANK in the sera of osteoporosis model rats [23]. erefore, it is important to elucidate the relationship among the QGY prescription, osteoporosis, secretin, and the OPG/RANKL/RANK pathway. We hypothesized that QGY regulates OPG/RANKL/RANK pathway by improving the level of secretin to treat the primary type I osteoporosis. e purpose of this study was to investigate the therapeutic effects of QGY on bone loss and its possible mechanisms. In this study, we analyzed the mRNA and protein levels of secretin, OPG, RANKL, and RANK by enzyme-linked immunosorbent assays (ELISA), western blotting, and real-time quantitative PCR (RT-PCR). Moreover, we confirmed that the release of secretin modulated the OPG/RANKL/RANK pathway using secretin knockdown and secretin overexpression models. We believe that the results of this study may help develop an effective treatment for osteoporosis. . According to the TCM decoction method, all herbs were mixed and boiled for 1.5 h, and the process was repeated three times. e decoction was concentrated to obtain a 1.5 g/mL concentration of crude drug and preserved at 4°C.

Adeno-Associated Virus (AAV) Construction and
Injection. Mice were injected with AAV vectors only once after 1 week of feeding, which mediated the efficient transduction of bone tissue at 5 × 10 10 vector genomes/g body weight, followed by QGY administration. In this study, AAVs were resuspended in phosphate-buffered saline (PBS) and used. AAV carrying the secretin-targeting shRNA or scramble control shRNA was purchased from Genomeditech (Shanghai, China). e shRNA sequence was as follows: secretin AAV, 5′-AAAAGGACCAAGGTCTGCTGTT-GATTTCGATCAACAGCAGACCTTGGTCC-3′.

Animal Grouping and QGY Administration.
Forty-two healthy C57BL/6 mice (specific-pathogen-free (SPF) level, 6 months of age, female) were acquired from Peking Weitong Lihua Experimental Animal Co., Ltd. All experimental procedures strictly adhered to the Institutional Animal Care and Use Committee (IACUC) of the Zhejiang Chinese Medical University and were approved by the Committee on the Ethics of Animal Experiments of Zhejiang Chinese Medical University (11223). Animals were raised in an SPF environment for 1 week prior to the start of the study and were kept in standard environmental conditions of temperature (26 ± 2°C), humidity (70%), and light (12 h/12 h light/dark). Adeno-associated virus (AAV) was injected into the caudal vein to establish mouse models of secretin knockdown and secretin overexpression [24]. e DNA sequence for secretin overexpression was combined with the AAV vector containing a green fluorescent protein (GFP) sequence, which was named the secretin overexpression group. Similarly, an AAV-secretin knockdown mouse model was developed by injecting AAV carrying the secretin-targeting shRNA. e same AAV vector containing only GFP was used as a control and named the secretin overexpression control or secretin knockdown control groups. After 1 week of adaptive feeding, the C57BL/6 mice were randomly distributed into seven groups as follows: control group (Con), secretin overexpression group (Sec-over), secretin overexpression control group (Sec-over-Con), secretin overexpression + QGY group (Sec-over + Q), secretin knockdown group (Sec-kno), secretin knockdown control group (Sec-kno-Con), and secretin knockdown + QGY group (Sec-kno + Q). After tail vein injection of AAV, the mice were adaptively fed for 3 weeks and then divided into the relevant groups. Next, QGY was given by gavage for 8 weeks. e QGY treatment group was perfused with QGY at a dose of 350 mg/10 g/d, whereas the control groups were orally administered normal saline. Peripheral blood was collected from the carotid artery, and bilateral femurs were cryopreserved in liquid nitrogen for future use.
2.5. ELISA. Serum was obtained from peripheral blood specimens using standard methods and stored at −80°C until analysis. e levels of secretin, OPG, RANKL, and RANK in serum were tested using ELISA kits (Jingmei Biological Company, China), according to the manufacturer's instructions. Absorbance was measured using a SpectraMax Plus 384 (Molecular Devices, USA) with a test wavelength of 450 nm.

RT-PCR Analysis.
e mRNA levels of secretin, OPG, RANKL, and RANK were determined by RT-PCR. Total RNA was separated from femurs using TRIzol Reagent (15596018, ermo, USA) according to the manufacturer's protocol, and then, 1.0 μg of total RNA was reverse transcribed into cDNA using the PrimeScript RT reagent kit (RR037B, Takara, Japan). PCR amplification was performed and normalized using TB Green ® Premix Ex Taq ™ (RR420B, Takara, Japan) and a 7500 Fast Real-time PCR instrument (Applied Biosystems, Foster City, USA). Relative mRNA levels were calculated using the 2 −ΔΔCt method [25]. e following primers were used for PCR:

Statistical Analysis.
All statistical analyses were conducted using SPSS 22.0 and GraphPad 8.0 Prism software. All data are expressed as the mean ± SD, and one-way analysis of variance was used to compare multiple groups. Statistical significance was set at P < 0.05.

QGY Enhances BMD and Improves Bone Microstructure.
BMD and bone structural properties have been widely used to evaluate the severity of osteoporosis. Our previous clinical studies have found that QGY significantly improves BMD in patients with osteoporosis and that secretin is a specific protein for osteoporosis. erefore, we speculated that QGY might reduce bone loss by affecting the expression of secretin [1,16]. To test this hypothesis, we examined the BMD and bone structural properties in Sec-over and Sec-kno mice. Micro-CTanalysis of the femoral shaft revealed that BMD was decreased in the Sec-kno group and increased in the Sec-over group, as compared to that in the control groups. Importantly, BMD was increased in the Sec-kno group after QGY intervention ( Figure 1). Furthermore, the morphological parameters of bone tissue indicated that the low expression of secretin resulted in a remarkable decrease in BV/TV and Tb.
and Tb. N and an increase in Tb. Sp. In addition, morphological changes in trabecular bone suggested improvements after treatment with QGY. e results indicated that QGY increased BMD and improved bone microstructure in the Sec-kno group of mice.

QGY Increases Expression of Secretin in Serum.
We initially developed two animal models, Sec-kno and Secover, to explore the specific relationship between QGY and secretin. e expression of secretin in mice from the Secover and Sec-kno groups showed significant changes ( * P < 0.05) compared with that in mice from the Con group. Secretin levels were found to be significantly higher after QGY treatment in the Sec-over group, as compared to that before intervention. It was, therefore, verified that QGY could promote the expression of secretin (Figure 2).

QGY and Secretin Have Positive Regulatory Effects on the OPG/RANKL/RANK Pathway in Serum.
e OPG/RANKL/ RANK pathway is a classic pathway involved in osteoporosis. To explore whether QGY and secretin were related to the OPG/RANKL/RANK signaling pathway, we monitored the expression of OPG, RANKL, and RANK in mouse serum from different treatment groups (by ELISA) after the oral administration of QGY. e results revealed that OPG was increased in mice of the Sec-over group and decreased in mice of the Sec-kno group ( * P < 0.05), as compared to that in the Con group, as shown in Figure 3(a). e expression of OPG was increased after the oral administration of QGY in the Sec-kno + Q group compared to that in the Sec-kno group ( # P < 0.05). Simultaneously, the expression of RANKL and RANK decreased in the Sec-over group and increased in the Sec-kno group as compared to levels in the Con group. e expression of RANKL and RANK was significantly decreased after treatment with QGY as compared to that in mice of the Sec-kno group (Figures 3(b) and  3(c)). e results verified that QGY could regulate the OPG/ RANKL/RANK pathway by affecting the expression of secretin in blood.

QGY Regulates the Expression of Secretin, OPG, RANKL, and RANK in Femur Tissue.
It is well known that OPG, RANKL, and RANK constitute a system that regulates bone remodeling. To further explore the potential mechanisms through which QGY exerts effect on osteoporosis, the expression of secretin, OPG, RANKL, and RANK in femur tissue was studied using western blotting and PCR. e western blot results are shown in Figure 4; secretin and OPG protein expression was downregulated, whereas that of RANKL and RANK was upregulated, in femur tissue from mice in the Sec-kno group. Meanwhile, the expression of secretin and OPG was increased and that of RANK and RANKL was decreased after treatment with QGY. PCR results are shown in Figure 5, and a similar trend was observed.

Discussion
Osteoporosis is a degenerative bone disease related to aging [7]. It is well known that an imbalance between bone resorption (driven by osteoclasts) and bone formation (driven by osteoblasts) results in decreased bone mass [5,26]. erefore, studying the mechanism of osteoblasts and osteoclasts could help in designing therapeutic strategies to effectively prevent a reduction in bone mass. Based on current knowledge, the OPG/RANL/RANK signaling pathway is the main regulator of osteoporosis and modulates osteoclast differentiation, induction, activation, and maintenance [19,27]. RANKL binds to RANK on the surface of osteoclasts to form a RANKL-RANK complex, which stimulates the activation, formation, and differentiation of osteoclasts [28,29]. OPG inhibits bone resorption by preventing the binding of RANKL to the RANK receptor by blocking the RANKL binding site [30,31]. It has been reported that increased levels of RANKL/OPG directly affect osteoclasts [32]. In addition, digestive system diseases, especially inflammatory bowel disease, are significantly associated with an increased risk of osteoporosis [33]. e mechanism leading to low BMD and increased fracture risk caused by inflammatory bowel disease is closely related to OPG/RANKL/RANK pathway. In several studies, it has been shown that inflammatory bowel disease results in increased concentrations of proinflammatory cytokines, including tumor necrosis factor alpha (TNF-α), interleukin (IL)-1β, IL-6, and IL-17 [34]. It is well known that cytokines such as TNF-α, IL-1β, and IL-6 increase the production of receptor activator of nuclear factor kappa B ligand (RANKL) by preosteoblasts, thereby promoting osteoclastogenesis.
ese cytokines bind to their respective receptors, which causes increased activation of the p38 mitogen-activated protein kinase (MAPK), nuclear factor kappa B (NF-κB), and c-Jun N-terminal kinase (JNK) pathways, resulting in excessive bone loss and reduced BMD [35]. On the other hand, digestive system diseases lead to malabsorption of trace elements such as calcium, vitamin C, and vitamin D [36]. Studies suggest that insufficient calcium intake leads to excessive parathyroid hormone production and decreased BMD. Vitamin C is absorbed in the small * * * #

Con
Sec-over Sec-over-con Sec-over+Q Sec-kno Sec-kno-con Sec-kno+Q intestine, and a higher consumption of this vitamin decreases the risk of osteoporosis by 33% [37]. Vitamin D acts directly and indirectly on the bone, regulates calcium absorption in the gastrointestinal tract, determines the proper serum calcium level, influences osteoblasts mineralization and differentiation, and decreases the synthesis of parathormone and the reabsorption of phosphates from bones [38]. Secretin protein was the first animal hormone discovered by Bayliss and Starling in 1902 [39]. e cells producing secretin were named "S" cells and are mainly distributed in the duodenal mucosa, with a small number in the jejunum, ileum, and gastric antrum. It was found that secretin has similar functions with PTH and calcitonin, and there is an important relationship between secretin and RANKL of the OPG/ RANKL/RANK pathway of bone metabolism, as shown by the string protein-protein interaction network analysis tool. e OPG/RANKL/RANK pathway and its downstream factors are relatively well studied, but it is not clear whether Relative mRNA level of RANKL

Con
Sec-over Sec-over-con Sec-over+Q Sec-kno Sec-kno-con Sec-kno+Q Relative mRNA level of RANK

Con
Sec-over Sec-over-con Sec-over+Q Sec-kno Sec-kno-con Sec-kno+Q (c) Figure 5: Effect of QiangGuYin (QGY) treatment on the mRNA expression of OPG, RANK, and RANKL in the femur tissue. (a-c) e mRNA expression of OPG was decreased, whereas that of RANK and RANKL was increased, after secretin knockdown. e mRNA expression of OPG was increased, whereas that of RANK and RANKL was significantly decreased, in the Sec-kno + Q group, compared to levels in the Sec-kno group. Data are expressed as the mean ± SD (n � 6). * P < 0.05 vs. the Con group; # P < 0.05 vs. the Sec-kno group. Groups were as follows: control group (Con), secretin overexpression group (Sec-over), secretin overexpression control group (Sec-over-Con), secretin overexpression + QGY group (Sec-over + Q), secretin knockdown group (Sec-kno), secretin knockdown control group (Seckno-Con), and secretin knockdown + QGY group (Sec-kno + Q). e protein expression of secretin and OPG was decreased, whereas that of RANK and RANKL was increased, after secretin knockdown. e expression of secretin and OPG was increased, whereas that of RANK and RANKL was significantly decreased in the Sec-kno + Q group, as compared to that in the Sec-kno group. Data are expressed as the mean ± SD (n � 3). * P < 0.05 vs. the Con group; # P < 0.05 vs. the Seckno group. Groups were as follows: control group (Con), secretin overexpression group (Sec-over), secretin overexpression control group (Sec-over-Con), secretin overexpression + QGY group (Sec-over + Q), secretin knockdown group (Sec-kno), secretin knockdown control group (Sec-kno-Con), and secretin knockdown + QGY group (Sec-kno + Q). there are influencing factors or proteins upstream of this pathway. Based on our experimental results, we propose a superior subordinate relationship, the secretin/OPG/ RANKL/RANK pathway. ese results also show the close relationship between digestive system diseases and osteoporosis.
Osteoporosis has always been a research hotspot in orthopedics and other related fields. ere are many methods to obtain animal models of osteoporosis, such as surgical induction [40] (removal of bilateral ovaries/testes), drug induction [41] (adrenocortical hormone and other drugs), loss induction, diet induction [42] (low calcium diet), gene knockout [43] (knockout of the OPG gene), and spontaneous degeneration. e molding time is generally about 2 months. In our previous study, we found that secretin is a specific serum protein marker for patients with type I osteoporosis [2]. erefore, we used the AAV-mediated knockdown of the secretin gene to obtain an osteoporosis mouse model. After 3 weeks, we evaluated whether the model was successful by determining BMD, an important index to evaluate bone strength [44]. According to the micro-CT analysis of femoral shafts, it was revealed that BMD was decreased in mice from the Sec-kno group and increased in mice from the Sec-over group as compared to that in mice from Con groups. Moreover, morphological parameters of bone tissue indicated that the low expression of secretin resulted in a significant decrease in BV/TV, Tb. , and Tb. N and an increase in Tb. Sp. is experiment showed that our method of obtaining an osteoporosis model with secretin knockdown was successful, and the required time was much shorter compared with previous methods. e morphological changes of trabecular bone were also improved after treatment with QGY. It was found that QGY improved BMD and bone microstructure in Sec-kno-group mice. us, it was important to elucidate the mechanistic details of the QGY function that led to an effective improvement in BMD in mice, as well as the role that secretin played in this process. ELISA and western blot assays indicated that QGY had a positive effect on bone loss in mice from the Sec-kno group. e results also indicated that BMD in mice from the Sec-kno group was significantly decreased but that the level of BMD was significantly improved after QGY treatment. It has been suggested that QGY has a positive effect on BMD by increasing the level of secretin. Moreover, RANKL, RANK, and OPG constitute a three-in-one system regulating bone resorption and bone formation [45]. In our experiments, we detected the protein expression of RANKL, RANK, and OPG in both serum and bone tissue by ELISA and western blotting, as well as the mRNA levels of RANKL, RANK, and OPG in bone tissue by RT-qPCR. e results demonstrated that knockdown of secretin inhibited OPG expression and promoted RANKL and RANK expression. After the oral administration of QGY, the expression of RANKL and RANK was inhibited, whereas OPG expression was upregulated. e experimental results support our hypothesis and revealed the antiosteoporosis mechanism of QGY. Our previous clinical trials found that QGY significantly improves BMD and reduces clinical symptoms such as low back pain in patients with osteoporosis but failed to reveal the underlying mechanism [1]. e present work used the new modeling method of secretin knockdown, which shows that QGY modulates the OPG/RANKL/RANK pathway by increasing secretin levels during treatment of primary type I osteoporosis.
Recent studies have shown that QGY is widely used in the prevention and treatment of osteoporosis [1]. Pharmacological studies have suggested that the effective components of natural Chinese medicine can attenuate the imbalance in bone resorption and reduce bone microstructural degradation [9,46,47]. Interestingly, the active ingredients in the prescription of QGY play an important role in preventing low bone mass [7,[48][49][50][51][52]. Polysaccharides and total flavonoids were determined to be the active ingredients of Hedysarum multijugum Maxim., which could reduce bone resorption and inhibit osteoclasts by increasing OPG levels and reducing RANKL and RANK levels considerably, thereby lowering the ratio of RANKL/OPG. Moreover, polysaccharides effectively inhibit osteoclastogenesis by reducing tumor necrosis factor-(TNF-) α [7,46]. Antler polypeptide is the active ingredient of Colla Cornus Cervi, which increases osteoblast differentiation and inhibits TNF-α-mediated inhibition of osteoclastogenesis via regulation of the NF-κB/p65 pathway, as well as suppressing the NF-κB signaling pathway in osteoclasts in vitro [48]. Flavonoids are the active ingredient of Drynariae Rhizoma, which decrease the levels of TNF-α and interleukin-6 (IL-6), involved in the induction of bone resorption, and suppress osteoclastogenesis by modulating the OPG/RANKL system in osteoblastic cells [48]. Eucommiae Cortex, rich in polyphenols such as phenolic acids, flavonoids, and lignans, is a kidney-tonifying herbal medicine that has been shown to improve bone histomorphometric parameters in lead-exposed rats. e possible mechanism underlying the effects of Eucommia Cortex extract is the promotion of bone formation by increasing the expression of OPG and the suppression of bone resorption by decreasing the level of RANKL in serum [22]. Gentiopicroside is the main active ingredient of Gentiana macrophylla Pall., which significantly suppresses RANKL-induced osteoclast formation in bone marrow mesenchymal stem cells (BMMs), with the possible mechanism being the inhibition of RANKL-induced activation of NF-κB and Jun N-terminal kinase (JNK) signaling pathways in BMMs [50]. Previous studies have also found that Saposhnikoviae divaricata chromone extract reduces the levels of NF-κB by suppressing p-p38, p-mitogen-activated protein kinase 1 (p-ERK), and p-JNK expression, which might be involved in the transcription of proinflammatory factors and the phosphorylation of the transcription factor NF-κB. Moreover, it could suppress the production of IL-6, TNF-α, and IL-1β in tissues [52].

Conclusion
Our findings demonstrate that QGY modulates the OPG/ RANKL/RANK pathway by increasing secretin levels during treatment of primary type I osteoporosis. In addition to the effect of QGY, we also revealed the general regulatory effect of secretin on the OPG/RANK/RANKL pathway, which has not been reported. erefore, QGY might be a promising drug candidate for treating primary type I osteoporosis. 8 Evidence-Based Complementary and Alternative Medicine Abbreviations QGY: QiangGuYin BMD: Bone mineral density TCM: Traditional Chinese medicine AAV: Adeno-associated virus OPG/RANK/ RANKL: Osteoprotegerin/receptor activator of the nuclear factor kappa B (NF-κB)/receptor activator of the NF-κB ligand ELISA: Enzyme-linked immunosorbent assay ELISA: Enzyme-linked immunosorbent assay RT-PCR: Real-time quantitative PCR SPF: Specific-pathogen-free GFP: Green fluorescent protein Con: Control group Sec-over: Secretin overexpression group Sec-over-Con: Secretin overexpression control group Sec-over + Q: Secretin overexpression + QGY group Sec-kno: Secretin knockdown group Sec-kno-Con: Secretin knockdown control group Sec-kno + Q: Secretin knockdown + QGY group PBS: Phosphate-buffered saline Tb. : Trabecular thickness Tb. N: Trabecular number BV/TV: Bone volume/tissue volume Tb. Sp: Trabecular separation PVDF: Polyvinylidene fluoride TNF-α: Tumor necrosis factor IL-6: Interleukin-6 BMMs: Bone marrow mesenchymal stem cells JNK: Jun N-terminal kinase.
Data Availability e data are reflected in the manuscript; for more information, contact the corresponding author: Dr. Wu Lianguo, email: zcmu@vip.163.com.

Conflicts of Interest
e authors have no conflicts of interest to declare.

Authors' Contributions
Xuefei Li, Longkang Cui, and Wenhua Chen contributed equally to this work.