Abstract

Many clinical trials and meta-analyses related to acupuncture for osteoporosis (OP) have been published. However, identifying the evidence from these studies still remains a challenge for acupuncturists. We conducted a systematic search of the Chinese Biomedical Medicine (CBM), VIP Database, Wanfang Data, China National Knowledge Infrastructure (CNKI), PubMed, Springer, Cochrane Library, and Embase to identify relevant trials, systematic reviews, and/or meta-analyses up to October 31, 2018. Data were extracted to assess the methodological quality using Veritas plots and to explore potential acupuncture prescriptions using the Traditional Chinese Medicine inheritance support system (TCMISS). In addition, potential mechanisms of core acupoints identified by data mining were summarized based on published studies. A total of 218 clinical trials and ten meta-analyses were included, involving 212 acupuncture prescriptions, 102 acupoints, 13 meridians, three extra meridians, and one Ashi point. The mean Veritas score of publication year, type of study, Assessment of Multiple Systematic Reviews 2, Preferred Reporting Items for Systematic Reviews and Meta-Analyses, heterogeneity, and publication bias were 5.5, 7.2, 6, 5.6, 5.8, and 7.5, respectively. The study of Pan et al. received the highest Veritas score of 8.67 points. The most frequently used meridian was BL. Acupoint combinations BL23 and BL20, BL23 and GV4, and BL23 and ST36 were used frequently. The core acupoints association networks were acupoints BL23, BL20, ST36, GV4, SP6, CV4, and KI3. The potential mechanisms of core acupoints involved upregulated expression of members in OPG/RANKL, Wnt/β-catenin, and MAPK pathways, such as LRP5, β-catenin, Runx2, and OPG. In conclusion, our Veritas plots enable acupuncturists to evaluate key attributes of meta-analysis quality related to acupuncture for primary OP and to improve the quality of evidence-based medicine relating to acupuncture. Data mining analysis revealed an association network of meridians, acupoint combinations, core acupoints, and the underlying mechanisms of acupuncture for primary OP.

1. Introduction

Osteoporosis (OP) is a disease characterized by low bone mineral density (BMD) and a high risk of fractures [1], with clinical manifestations of pain, muscle weakness, and fractures [2]. Approximately 9 million people worldwide suffer from OP or fragility (low trauma) fractures each year [3]. Previous studies have found that acupuncture could significantly increase the BMD of lumbar vertebrae, the values of trabecular area, and trabecular bone number, as well as reducing trabecular separation in postmenopausal OP rats [4]. Similar but slightly different results occur simultaneously in patients with OP [5]. In recent years, the number of patients seeking complementary and alternative therapies, including acupuncture, has grown continuously in clinical practice [6], leading to a large number of clinical trials and meta-analyses related to acupuncture for OP being published in various journals and databases [5, 7]. However, these studies are influenced by a series of biases and increasingly include nonrandomized clinical trials [8]. Identifying quality evidence and potential treatment strategies from these studies and applying it to clinical practice still remains a challenge for clinicians [9, 10]. Therefore, identifying clinical evidence and exploring treatment strategies of acupuncture has become an urgent problem.

The Veritas plot is a graphic tool adapted from radar plots to describe multiattribute data [10]. This technique was first applied to identify the highest quality meta-analytical literature of cardiac surgery [10]. As an evidence-synthesis graphic tool, Veritas plots are now widely used to identify and interpret variability in meta-analyses and to help clinicians use the highest quality evidence of acupuncture in clinical practice [11]. However, Veritas plots have not been applied to acupuncture for primary OP. TCMISS is a data mining tool invented by the Institute of Chinese Materia Medica China Academy of Chinese Medical Sciences for exploring potential treatment strategies of acupuncture or Chinese medicine [12, 13]. It includes different functions such as text mining, association rules analysis, and complex system entropy methods, making it easy and convenient for researchers without algorithmic experience to use it for data mining [14]. It can be used to guide clinicians in selecting appropriate acupoints or Chinese herbs to treat specific diseases [14, 15]. However, data mining has not been conducted to explore potential acupuncture prescriptions for primary OP.

Given the uncertainty about the quality of meta-analyses and the treatment strategy of acupuncture for OP, we used Veritas plots to assess quality evidence of acupuncture and used TCMISS to explore potential treatment strategies of acupuncture for primary OP. In addition, we summarized the major antiosteoporotic effects and potential mechanisms of core acupoints identified by data mining. These will help clinicians establish a specific acupuncture prescription for OP and use it in clinical practice.

2. Methods

2.1. Participants

Patients with primary OP were included as participants in this study. The diagnostic criteria for primary OP were based on the 1999 World Health Organization criteria: BMD is measured by dual energy X-ray absorptiometry as a T-score below -2.5 standard deviations (SD) in the femur neck, lumbar spine, or total femur region [16]. There were no restrictions on age, sex, race, language, or region.

2.2. Selection Criteria

This study included randomized control trials (RCTs) and/or clinical control trials (CCTs) and systematic reviews and/or meta-analyses of acupuncture for OP. Cross-sectional studies, comments, cohort studies, animal experiments, and reviews were excluded.

2.3. Search Strategy

A systematic electronic search was conducted of CBM, VIP Database, Wanfang Data, CNKI, PubMed, Springer, Cochrane Library, and Embase for relevant RCTs and/or CCTs and systematic reviews and/or meta-analyses. All of the databases were searched from their inception to October 31, 2018. The search algorithm for PubMed was as follows: (“acupuncture” or “electroacupuncture” or “needle warming therapy” or “eye acupuncture” or “foot acupuncture” or “wrist-ankle acupuncture” or “hand acupuncture” or “hand-foot acupuncture” or “face acupuncture” or “nose acupuncture” or “Ashi points”) and (“age-related osteoporosis” or “senile osteoporosis” or “postmenopausal osteoporosis” or “primary osteoporosis”) and (“randomized controlled trial” or “trial” or “clinical trial” or “systematic review” or “meta-analysis”), with no restriction on subheadings. Similar but adapted search terms were used in other databases. Clinical trials published in abstract form were excluded unless sufficient information could be obtained from the abstract or authors. Reference lists of the retrieved papers were also reviewed to discover other studies, and studies not included in the databases mentioned previously were additionally searched.

2.4. Data Selection

The results from different databases were imported into NoteExpress 3.2.0 (http://www.inoteexpress.com/aegean/). Duplicate studies were then identified and deleted with NoteExpress 3.2.0. Two authors conducted data selection independently (Lin HX and Wang XT), reviewed the abstracts and full texts, and extracted the relevant information from the included studies. Any discrepancies were resolved by consensus.

2.5. Data Extraction

Two reviewers (Wang XT and Lin HX) extracted study information for each eligible study, including the first author names, publication year, data sources, time range for inclusion in the study, detailed search strategies, search process, sample sizes, number of patients, bias assessment tools, outcomes, and methods of reporting publication bias. Data was checked by Yingjie Mo.

2.6. Veritas Score Indicators

Veritas scores were determined for publication year, study type, Assessment of Multiple Systematic Reviews (AMSTAR) 2 score, Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) score, heterogeneity, and publication bias to assess the quality of included systematic reviews or meta-analyses. The year a study was published was an important factor of heterogeneity, as disease and acupuncture techniques may change over time, whereas technical expertise may have led to unfavorable results in early studies [17]. The AMSTAR 2, updated in 2017, could provide readers with a better summary assessment to determine whether the systematic review has met the methodological safeguards against bias [18]. PRISMA statement could accurately reflect the ability of researchers to write in a comprehensible manner, instead of the way they conducted the review [19, 20]. Heterogeneity in study design or intervention has a significant impact on the conduct and outcomes of a meta-analysis [21] and often affects acupuncturists in some important decisions, such as whether a treatment is applicable to all targeted people [22]. Publication bias is an inherent problem in meta-analyses because some studies are not published in indexed journals, and some negative results were not publicly available [23, 24].

2.7. Ranking of Evidence

The meta-analyses were ranked according to six categories: year of publication, study type, AMSTAR2 score, PRISMA score, heterogeneity, and publication bias. For publication year, if the publication year of multiple studies was the same, we ranked studies by the time range for inclusion in the study. In the AMSTAR 2 score, if the study meets the criteria of AMSTAR 2 item, it scores 1 point in this item; 0.5 points for those that partially satisfy the criteria; otherwise 0 points. In the PRISMA score, 2 points were given for studies that completely meet the criteria of the PRISMA item; 1 point for those that partially satisfy the criteria; otherwise 0 points. The heterogeneity score is the average of the heterogeneity scores for clinical outcomes. Heterogeneity among studies was calculated using the chi-square test and I2 statistic [25]. 0% ≤ I2 ≤ 50% indicated no heterogeneity and the score is 3 points, 50% < I2 ≤ 75% was considered slightly significant heterogeneity and scores 2 points, and I2 > 75% was regarded as significant heterogeneity and scores 1 point.

The scoring system was applied as follows [10]: in each category, the best study got the highest score of n, where n = the number of studies. The second best study got n-1 points, and so on. If two studies performed equally well, the included study with the next highest score would get n-2 points. We considered the Veritas score as a summary statistic. The Veritas score of each study is the average scores in the six dimensions of quality.

2.8. Data Mining Analysis

Data were analyzed by TCMISS Version 2.5 (invented by Institute of Chinese Materia Medica China Academy of Chinese Medical Sciences). We first analyzed the frequency of acupoints, meridians, and common acupoint combinations using text mining. Then we conducted an association rule analysis to determine the support value of different acupoints and the core acupoints association network of acupuncture for primary OP. We set the “the degree of support” to 60, which means the frequency of A and B appearing at the same time is 60. The “confidence level” was set to 0.6, which means that when A occurs, the probability of B appearing is 60%. We also further analyzed the cooccurrence matrix of the top 23 acupoints using Heml 1.0.3.7 (http://hemi.biocuckoo.org/down.php). The results were represented by a heat map.

2.9. Summary of Major Antiosteoporotic Effects and Potential Mechanisms of Core Acupoints

The major antiosteoporotic effects and underlying mechanisms of core acupoints identified by data mining were summarized using published studies, thereby providing evidence for the patients and acupuncturists of the applications of acupuncture for primary OP

3. Results

3.1. Eligible Studies

The research process is summarized in Figure 1. We firstly identified 2204 citations through an initial search and identified 821 potential articles after removing duplicated records. After screening the titles and abstracts, 558 studies were deleted. We included 228 studies after screening the full text of these 558 studies. In total, 218 studies were clinical trials and the acupoint characteristics were analyzed using TCMISS Version 2.5. Ten studies were meta-analyses and the evidence was ranked using Veritas scores. The baseline characteristics of the meta-analyses are presented in Table 1. The meta-analyses were published between 2014 and 2018, with data from 12 databases including Embase, Medline, the Cochrane Library, ScienceDirect, PubMed, AMED, CBM, VIP Database, Wanfang Data, CNKI, Web of Science, and The Cochrane Central Register of Controlled Trials. The databases were searched from their inception to October 31, 2018, in the latest meta-analysis. A total of 138 clinical trials involving 11,232 patients were included in 10 meta-analyses. Five meta-analyses provided a detailed search strategy. Six meta-analyses had search processes that were performed independently by two reviewers. In the bias assessment, five meta-analyses used Cochrane Collaboration’s risk of bias tool, while others used the Jadad scale. Four meta-analyses did not report publication bias.

3.2. Veritas Score

The scores of included meta-analyses are presented in Table 2. The detailed heterogeneity score, AMSTAR 2 score, and PRISMA score of included studies are shown in Supplementary Tables 1, 2, and 3, respectively. We also developed summary Veritas plots of the 10 studies (Figure 2). The mean Veritas scores of publication year, type of study, AMSTAR 2 score, PRISMA score, heterogeneity, and publication bias, were 5.5, 7.2, 6, 5.6, 5.8, and 7.5, respectively. The number of meta-analyses that exceeded the mean Veritas scores of publication year, type of study, AMSTAR 2 score, PRISMA score, heterogeneity, and publication bias was 5, 6, 5, 6, 5, and 5, respectively. The study with the highest Veritas score was by Pan et al. (2018), with 8.67 points, while the lowest Veritas score was Li et al. (2014), with 4 points. There were 5 studies with Veritas scores ≥ 6.27 points. In Figure 3, the individual Veritas scores for all of the meta-analyses are plotted.

3.3. Frequency Statistics of Acupoints

A total of 212 acupuncture prescriptions were included, which involved 102 acupoints. The top 30 acupoints of acupuncture for OP are presented in Table 3. Acupoints BL23, BL20, ST36, and GV4 were used in more than 60% of the acupuncture prescriptions, with frequency of 178, 120, 112, and 111, respectively. There were 10 acupoints with frequency ≥ 50, which were BL23, BL20, ST36, GV4, SP6, CV4, KI3, GV3, BL40, and BL11.

3.4. Frequency Statistics of Meridians

Meridians and acupoints used in acupuncture therapy for primary OP are shown in Table 4. A total of 13 meridians, three extra meridians, and one Ashi point were used in acupuncture for primary OP. There were six meridians used over 100 times, including stomach meridian of foot-yangming (ST), spleen meridian of foot-taiyin (SP), governor vessel (GV), gallbladder meridian of foot-shaoyang (GB), conception vessel (CV), and bladder meridian of foot-taiyang (BL). The most frequently used meridian was BL, which was used 640 times and involved 31 acupoints. The second most commonly used meridian was GV, which was used 269 times and involved 13 acupoints. CV was the third most frequently used meridian, with 147 times and nine points.

3.5. Common Acupoint Combinations

The common acupoint combinations of acupuncture for primary OP are presented in Table 5. A total of 15 common acupoint combinations were frequently used over 60 times. There were three acupoint combinations with frequency ≥ 100: BL23 and BL20, BL23 and GV4, and BL23 and ST36. The acupoint combination BL23 and BL20 was used 117 times in acupuncture therapy for OP, accounting for about 50% of acupuncture prescriptions. The second most frequently used acupoint combination was BL23 and GV4, with 104 times. Acupoints BL20 combined with ST36 and BL23, BL20 combined with ST36 were used as many as 76 times. These results were consistent with the cooccurrence matrix of acupoints, which are shown in Figure 4. There were 18 pairs of acupoints with a cooccurrence frequency above 50, of which four pairs were above 75.

3.6. Association Rules of Acupoints

We present the association rules for acupoints in Table 6. There were 22 acupoint combinations with confidence levels ranging from 60% to 98.7%, seven of which were above 90%, such as acupoint BL20, ST36 -> BL23; BL20 -> BL23; CV4 -> BL23; BL20, GV4 -> BL23; SP6 -> BL23; GV4 -> BL23; KI3 -> BL23. The acupoint combinations of BL20 and ST36 with BL23 have the highest confidence level of 98.7%, which means that when BL20 and ST36 were selected, the probability of selecting BL23 was 98.7%. The acupoint combination of BL20 and BL23 has the second highest confidence level of 97.5%. The core acupoints association network of acupuncture for primary OP is shown in Figure 5, including acupoints BL23, BL20, ST36, GV4, SP6, CV4, and KI3. The degree of support from strong to weak was acupoint BL23 combined with BL20, BL23 with ST36, BL23 with GV4, BL20 with ST36, BL23 with SP6, BL23 with KI3, ST36 with GV4, ST36 with SP6, BL23 with CV4, and SP6 with BL20.

3.7. Beneficial Effects of Core Acupoints on Primary OP and Potential Mechanisms

Representative examples of major antiosteoporotic effects of core acupoints and potential mechanisms are presented in Table 7. The prevention of primary OP in acupoints BL23 and BL20 was mainly to reduce the loss of bone mass, increase femoral BMD, and improve the microstructure of bone tissue. The role of acupoints ST36, SP6, CV4, and KI3 in preventing OP was mainly to increase E2 levels and BMD. The antiosteoporotic effects of acupoints GV4 and CV4 were to reduce the loss of bone mass and increase the femoral BMD. The potential antiosteoporotic mechanisms in acupoints BL23, BL20, and GV4 were to upregulate the ratio of OPG/RANKL, and LRP5, β-catenin, and Runx2 expression. In addition, acupoints BL23 and BL20 also could downregulate the expression of phosphorylated (p)-p38 and p-JNK. The potential antiosteoporosis mechanism in acupoints SP6 and KI3 was to decrease the concentration of TNFα and IgM. Acupoint ST36 could reduce OPGL expression and upregulate the expression levels of RANKL mRNA. Acupoint CV4 could upregulate the expression levels of ER-α and ERA-α mRNA in uterus, hypothalamus, and spleen.

4. Discussion

In our study, we presented a novel application of Veritas plots that acupuncturists can use to evaluate key attributes of meta-analysis quality without extensive narratives. A total of 10 meta-analyses were included. The mean Veritas score of publication year, type of study, AMSTAR 2, PRISMA, heterogeneity, and publication bias were 5.5, 7.2, 6, 5.6, 5.8, and 7.5, respectively, indicating that the included meta-analyses performed well in the following areas: timely update results, high consistency of study type, compliance with PRISMA statements, and low heterogeneity. However, the statement of the review method, the assessment of the potential impact of risk of bias in individual studies on the results, the list of excluded studies, funding sources, and conflicts of interest in AMSTAR 2 were insufficient. Some meta-analyses did not assess publication bias, indicating that we should pay more attention to the methodological guarantee of bias, and consider the clinical value of the results carefully. A similar phenomenon also occurs in acupuncture treatment of anxiety disorders [41], which indicates a direction to further improve the quality of research. With the help of the Veritas scores, we ultimately found that the study of Pan et al. [26] performed well in the year of publication, type of study, AMSTAR 2 score, PRISMA score, and publication bias and received the highest Veritas score of 8.67 points.

In this study, we used TCMISS to investigate the meridians, acupoint combinations, and core acupoints for primary OP. A total of 212 acupuncture prescriptions were included, including 102 acupoints, 13 meridians, three extra meridians, and one Ashi point. The most frequently used meridian was BL, which was used 640 times and involved 31 acupoints. The second most frequently used meridian was GV. Experiments in ovariectomized rats have demonstrated that acupuncture at BL and GV acupoints could alleviate osteoporosis by regulating the OPG/RANKL and Wnt/beta-catenin signaling pathways [38]. Acupoint combination BL23 and BL20 was used in about 50% of acupuncture prescriptions, with a frequency of 117. The second most frequent acupoint combination was BL23 and GV4, and the third was BL23 and ST36. A clinical study has demonstrated that electroacupuncture stimulation at BL20/BL23 could increase levels of osteocalcin and the BMD of lumbar vertebrae, as well as improving bone microstructure in the femur [4]. The core acupoints association network of acupuncture for primary OP included BL23, BL20, ST36, GV4, SP6, CV4, and KI3. BL23 and KI3 [6] and BL23 and GV4 [42] have been demonstrated in randomized controlled trials to improve quality of life and reduce pain in patients with primary OP. At the same time, we summarized the major antiosteoporotic effects and potential mechanisms of core acupoints based on published studies. The major antiosteoporotic effects of core acupoints were increased BMD and reduced loss of bone mass. The underlying antiosteoporotic mechanisms were to upregulate the expression of members in OPG/RANKL, Wnt/β-catenin, and MAPK pathways, such as LRP5, β-catenin, Runx2, and OPG.

5. Limitations

There are some limitations in this study. First, we did not assess the impact of population factor on the results, which is different from the Panesar et al. study [10]. Second, the number of clinical trials, especially the multicenter randomized controlled clinical trials that validated the experience of well-known Chinese medicine practitioners, is small. This may limit the applications of the results. Third, the acupoint combinations and core acupoints of acupuncture for primary OP were rarely confirmed by large-scale multicenter clinical trials or animal experiments.

6. Conclusion

Our Veritas plots enable acupuncturists to evaluate key attributes of meta-analysis quality related to acupuncture for primary OP and to improve the quality of evidence-based medicine related to acupuncture. Data mining analysis revealed an associated network involving meridians, acupoint combinations, core acupoints, and the underlying mechanisms of acupuncture for primary OP. Future large-sample multicenter RCTs and animal experiments should be conducted to confirm our findings.

Abbreviations

CBM:Chinese Biomedical Medicine
CNKI:China National Knowledge Infrastructure
AMSTAR:Assessment of Multiple Systematic Reviews tool
PRISMA:Preferred Reporting Items for Systematic Reviews and Meta-Analyse
RCT:Randomized control trial
CCT:Clinical control trial
BMD:Bone Mineral Density
BGP:Bone Gla protein
VAS:Visual analog scale
TCM:Traditional Chinese Medicine
E:Experimental group
C:Control group
OP:Osteoporosis
TCMISS:Traditional Chinese Medicine inheritance support system
AMSTAR:Assessment of Multiple Systematic Reviews
SD:Standard deviations
BL:Bladder meridian of foot-taiyang
GV:Governor vessel
CV:Conception vessel
GB:Gallbladder meridian of foot-shaoyang
ST:Stomach meridian of foot-yangming
SP:Spleen meridian of foot-taiyin
KI:Kidney meridian of foot-shaoyin
EX-B:Points of back
LR:Liver meridian of foot-jueyin
LI:Large intestine meridian of hand-yangming
LU:Lung meridian of hand-taiyin
SI:Small intestine meridian of hand-taiyang
PC:Pericardium meridian of hand-jueyin
EX-CA:Points of chest and abdomen
TE:Triple energizer meridian of hand-shaoyang
EX-UE:Points of upper extremities
EX-LE:Points of lower extremities
JNK:c-Jun N-terminal kinase
OPGL:Osteoprotegerin ligand
RANKL:Receptor activator of NF-kB ligand
Ig:Immunoglobulin
E2:Oestradiol
OPG:Osteoprotegerin
TNF:Tissue necrosis factor
ALP:Alkaline phosphatase
DPD:Deoxypyridinoline
ER-α:Estrogen receptor-α
LRP:Lipoprotein receptor-related protein
Runx:Runt-related transcription factor.

Data Availability

The detailed heterogeneity score, AMSTAR 2 score, and PRISMA score of included studies were shown in Additional Tables 1–3. The prescriptions were available from the corresponding author on reasonable request. The datasets or software generated and/or analyzed during the study were available in the Chinese Biomedical Medicine: http://www.sinomed.ac.cn/cross/; VIP Database: http://vip.hbsti.ac.cn/; Wanfang Data: http://www.wanfangdata.com.cn/index.html.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Haixiong Lin and Xiaotong Wang contributed equally to this work into All authors read and approved the final manuscript and agreed to submit it for publication in Evidence-Based Complementary and Alternative Medicine. Haixiong Lin and Xiaotong Wang conceived and designed the study. Feng Huang and Yongjun Chen revised the study. Haixiong Lin and Xiaotong Wang extracted the data, conducted statistical analysis, and wrote the manuscript. Yingjie Mo checked the data. Nenggui Xu, Feng Huang, and Yongjun Chen interpreted the results. Haixiong Lin, Xiaotong Wang, and Yingjie Mo wrote this manuscript. Chunni Lin designed picture. Nenggui Xu reviewed and proposed advice. Feng Huang and Yongjun Chen supervised the manuscript.

Acknowledgments

This work was supported by Guangdong Province Universities and Colleges Pearl River Scholar Funded Scheme (2016), Scientific Research and Innovation Team Program of Guangzhou University of Chinese Medicine (2017KYTD03), Innovation Team Program of Guangdong Provincial Department of Education (2017KCXTD006), Key Laboratory of Acupuncture and Moxibustion of Traditional Chinese Medicine in Guangdong (2012A061400017), Innovation and Strengthen School Project of Guangdong Provincial Education Department-National Major Cultivation Project (2014GKXM031), and Huang Feng Guangdong Provincial Famous Chinese Medicine Expert Heritage Studio (Index no. 006939748/2017-00583). All funding sources were used for data collection.

Supplementary Materials

Additional Table 1: heterogeneity score of included studies. Additional Table 2: the Assessment of Multiple Systematic Reviews (AMSTAR) 2 tool. Additional Table 3: Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) 2009 Checklist. (Supplementary Materials)