Sung Huang Laurent Tsai1,2†
Ching-Wei Hu1,2†
Shih-Chieh Shao3
Eric H. Tischler4Olufunmilayo H. Obisesan5,6Dominique Vervoort5
Wei Cheng Chen2,7Jiun-Ruey Hu8
Liang-Tseng Kuo2,9*- 1Department of Orthopaedic Surgery, Keelung Chang Gung Memorial Hospital, Keelung City, Taiwan
- 2School of Medicine, Chang Gung University, Taoyuan, Taiwan
- 3Department of Pharmacy, Keelung Chang Gung Memorial Hospital, Keelung City, Taiwan
- 4Department of Orthopaedic Surgery and Rehabilitation Medicine, Downstate Medical Center, State University of New York, New York, NY, United States
- 5Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, United States
- 6Johns Hopkins School of Medicine, Baltimore, MD, United States
- 7Department of Orthopaedic Surgery, Chang Gung Memorial Hospital, Taoyuan, Taiwan
- 8Department of Medicine, Vanderbilt University Medical Center, Nashville, TN, United States
- 9Division of Sports Medicine, Department of Orthopaedic Surgery, Chang Gung Memorial Hospital, Chiayi, Taiwan
Importance: Previous studies have shown the effectiveness and safety of direct oral anticoagulants (DOACs), including lower fracture risks, compared to warfarin. However, direct or indirect comparisons between different DOACs are scarce in the literature.
Objective: This study aims to compare fracture risks among different DOACs and warfarin, including apixaban, rivaroxaban, dabigatran, and edoxaban, in patients with non-valvular atrial fibrillation (NVAF) or venous thromboembolism (VTE).
Methods: We searched PubMed/MEDLINE, Embase, Cochrane CENTRAL, and Web of Science for randomized controlled trials and cohort studies comparing the fracture risks among patients who used warfarin or DOACs, up to March 2021. Two authors extracted data and appraised the risk of bias of included studies. The primary outcome was fracture risk. We performed pairwise meta-analyses to compare differences between medications and network meta-analyses using frequentist random-effects models to compare through indirect evidence. We used surface under the cumulative ranking curve (SUCRA) and mean ranks to determine the probability of a DOAC ranking best in terms of fracture risk.
Results: Thirty-one studies were included in the final analysis. Twenty-four randomized controlled trials and seven cohort studies with 455,343 patients were included in the systematic review and network meta-analysis. Compared to warfarin, the risk of any fractures was lowest with apixaban [relative risk (RR) = 0.59; 95% confidence interval (CI): 0.48–0.73], followed by rivaroxaban (RR: 0.72; 95% CI: 0.60–0.86), edoxaban (RR: 0.88; 95% CI: 0.62–1.23), and dabigatran (RR = 0.90; 95% CI: 0.75–1.07). No substantial inconsistency between direct and indirect evidence was detected for all outcomes.
Conclusions: All DOACs were safer than warfarin concerning the risk of fracture; however, apixaban had the lowest relative risk of fracture within the class of DOACs. Further head-to-head prospective studies should confirm the comparative safety profiles of DOACs regarding fractures.
Key Points
Question: What is the comparative risk of fractures in patients using different direct oral anticoagulants (DOACs) and warfarin?
Findings: This systematic review with network meta-analysis including 31 studies with 463,495 patients found that, compared to warfarin, the risk of fracture was lowest with apixaban, followed by rivaroxaban, edoxaban, and dabigatran. Our results suggested that among DOACs, apixaban carried the lowest fracture risk.
Meaning: DOACs were safer than warfarin with regard to the risk of fracture. Among the DOACs, apixaban had the lowest relative risk of fracture. Healthcare professionals should be informed about different fracture risk profiles associated with different DOACs in order to select the most appropriate DOACs for patients.
Introduction
As society ages, the prevalence of musculoskeletal and cardiovascular comorbidities increases. Osteoporosis, increasing with age (1), can increase the risk of osteoporotic fracture and subsequent death and disability in the older population (2). The incidence of non-valvular atrial fibrillation (NVAF), another concern in the elderly, continues to increase globally (3). Oral anticoagulants, including vitamin K antagonists (VKAs) and direct oral anticoagulants (DOACs), are recommended for patients with NVAF for the treatment or prevention of stroke and thromboembolism (4).
Warfarin, a classic VKA, has been the mainstay treatment for stroke prevention in patients with AF for decades. Of note, VKA use has been associated with an increase in osteoporotic fragility fractures (5–9). Great concern has been raised by a Medicare population-based study (5), in which AF patients using warfarin for longer than one year show an elevated risk of fragility fracture, compared to those not using warfarin. Bone quality is compromised due to the inhibition of vitamin K-dependent carboxylation of bone metabolism-associated proteins such as osteopontin and matrix Gla (10–15). Despite the potential risk of fragility fracture, warfarin has remained necessary for decades due to the lack of alternatives.
DOACs, recently approved for stroke prevention in AF patients, have been introduced for use as an alternative to warfarin. Given at least equal efficacy in stroke prevention and additional advantages including lower bleeding risk and reduced monitoring requirement compared to warfarin (4), the guidelines of the American College of Cardiology/American Heart Association and the Heart Rhythm Society currently recommend DOACs over warfarin for stroke prevention in NVAF patients (16–18). Consequently, in the United States, DOACs are now more common than VKAs in cardiovascular management (4). Furthermore, DOACs have not been reported to affect bone metabolism proteins (19). Binding et al. report that among 37,350 patients receiving DOACs for over 180 days with no previous use of osteoporotic medications, DOACs are associated with a significantly lower risk of any major osteoporotic fractures, compared to VKAs (6).
As DOACs continue to be a commonplace medication among elderly patients, it is essential to assess the comparative safety profiles, most notably with regard to fractures, within this drug class. Although recent studies have compared fracture risks among the OACs (20–22), the optimal choice of DOAC remains uncertain. Therefore, we performed this systematic review and network meta-analysis to evaluate the network, direct and indirect effects of fracture risk among different DOAC users.
Methods
Research Protocol and Search Question
The PICO search protocol framework was followed to address the hypothesis: DOAC use in patients with NVAF or VTE (Population of interest), can lead to a varying reduction in the risk of fractures, depending on which individual DOAC medication is used (Comparator/Intervention). Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) statement guidelines were followed for study protocol review and the study was registered in PROSPERO (CRD42020206788).
Eligibility Criteria and Primary Outcome
Studies were eligible if they met the following criteria: (1) They included adult patients using DOACs for NVAF or VTE. (2) They were observational studies or randomized controlled trials (RCTs). (3) They compared the fracture risk between DOACs and warfarin or other DOACs. Relevant exclusion criteria included: (1) single-arm studies, case reports, small case series of <10 patients, reviews, basic science experiments and animal- or cadaver studies; (2) studies including patients with severe infection or under immunosuppression; and (3) conference abstracts without corresponding full-length papers.
Search Strategy and Study Selection
On March 27th, 2021, we systematically searched PubMed/MEDLINE, Embase, Ovid, Cochrane Central Register of Controlled Trials (CENTRAL), Web of Science and Scopus for articles using the combination of keywords and medical subject heading (MeSH), adjusted for each database, including: “atrial fibrillation,” “anticoagulant,” “direct oral anticoagulant,” “non-vitamin K antagonist oral anticoagulants,” “vitamin K antagonist oral anticoagulants,” “ warfarin,” “Dabigatran,” “Pradaxa,” “Rivaroxaban,” “Xarelto,” “Apixaban,” “Eliquis,” “Edoxaban,” “Savaysa,” “non-vitamin K antagonist oral anticoagulants,” “novel oral anticoagulants,” “new oral anticoagulants,” “factor Xa inhibitors,” “factor IIa inhibitors,” “fracture,” “osteoporosis” and “osteoporotic fractures.” We also searched the reference lists of the included studies to identify additional studies, and the trial register (clinicaltrials.gov) for any ongoing trials. In addition, we contacted specialists in the field for any ongoing trials or unpublished data. We applied no language restrictions. The detailed search strategy is presented in the Supplementary Table 1.
Two reviewers (SHLT, CWH) independently evaluated eligible studies by their titles and abstracts and then reviewed the full text of relevant articles for further qualification. All disagreements between reviewers were resolved by reaching a consensus through discussion, and a third reviewer (LTK) was consulted where necessary.
Data Collection and Quality Assessment
Two independent reviewers (SHLT, CWH) extracted all data onto a pre-planned Microsoft Excel spreadsheet (version 16.32). Data fields included study characteristics (authors, year of publication, region of study, data source, study design, period of study), study arms, sample size of overall study and, by study arms, patient age, outcome as defined above, inclusion criteria of each study, specific definition of treatment arm, and source of funding.
The quality of included studies was assessed by two independent reviewers (SHLT, CWH). We evaluated all included RCTs via the RoB (Cochrane risk-of-bias tool for randomized trials) (23), and the non-RCTs via the Newcastle-Ottawa Scale (24). Grade assessment was also performed (25). All discrepancies were resolved by discussion, and a third reviewer (LTK) was consulted where necessary.
Statistical Analysis and Quantitative Data Synthesis
All statistical analyses were undertaken using Network commands for statistical software package Stata (Version 15). A pairwise function was first used to transform raw data to a contrast-based format and generate treatment effect and standard error for each pairwise comparison. A network meta-analysis was then performed to estimate network meta-analysis models with a frequentist approach derived from graph theoretical methods. The random-effects model was incorporated by adding the estimated heterogeneity τ2, based on the Dersimonian-Laird estimator (26). Subsequently, we examined the structure of our network comparison by applying the netgraph function, with vertices demonstrating treatments and the thickness of edges corresponding to the number of studies.
As a conservative assumption, a random-effects pooled relative risk (RR) with a 95% confidence interval (CI) was calculated to summarize the efficacy of each treatment. Forest plots were constructed to display findings with VKAs as the reference group. Given the I2 value increased with the larger populations included in the meta-analysis, τ2 was used to measure heterogeneity; 0.04, 0.16, and 0.36 corresponded to a low, moderate, and high degree of heterogeneity, respectively. Subgroup analysis based on treatment comparison was conducted to evaluate heterogeneity within studies. Sensitivity analysis was also performed in the presence of publication bias or significant heterogeneity. We also estimated the probabilities of each treatment being at each rank for each outcome. We obtained a treatment hierarchy using the surface under the cumulative ranking curve (SUCRA) and mean ranks; the SUCRA value is 0 when a treatment is the worst option and 1 when a treatment is the best option (27).
Furthermore, we assessed the potential inconsistency between direct and indirect comparisons using the design-by-treatment interaction model (28), and side-splitting models (29). The design-by-treatment interaction model provides a global assessment of consistency across the entire network. The side-splitting method separates evidence into direct and indirect evidence and then evaluates differences between them (28, 29). We used the Egger's test and a funnel plot to assess small-study bias (30, 31). Symmetry around the effect estimates line indicated lower chance of publication bias or small study effects (32).
Subgroup Analyses
Where data were available, we planned to perform subgroup analyses including:
1. Fracture location: spinal fracture, hip fracture, and all fractures.
2. DOAC indications: NVAF or VTE/PE.
3. Type of study design: RCTs vs. Non-RCTs.
4. Studies with a given drug dose.
5. Studies with male predominance.
6. Studies with patients aged <65.
Results
Literature Search and Selection Process
A total of 9,332 articles were identified through the database search. After the removal of duplicates, 1,149 articles remained. An additional 13 articles were identified after checking the reference lists of eligible studies. One thousand one hundred and seventeen articles were excluded by checking the titles and abstracts. After checking the full-text of the remaining 45 articles against the inclusion and exclusion criteria, 11 articles were excluded, whereby eight had the wrong study design, two had the wrong patient population, and one had the wrong outcomes (Supplementary Table 2). Ultimately, 31 studies were included in the network meta-analysis (Figure 1).
Figure 1. PRISMA flow diagram of the study.
Study Characteristics, Cohort Description and Treatment Definition
Our network meta-analysis included 24 randomized controlled trials and seven cohort studies with a total of 455,343 patients receiving five different anticoagulants. Two hundred twenty-one thousand two hundred three patients used warfarin, 78,810 used dabigatran, 106,996 used rivaroxaban, 35,359 used apixaban and the remaining 12,975 patients were edoxaban users. The network graphs are presented in Figures 2A–C, and the main characteristics of the included studies are reported in Table 1. The included studies were conducted in Asia (six studies; 50,203 patients), the Americas (three studies; 270,202 patients), Europe (one study; 14,376 patients), and multinational settings (21 studies; 120,562). The included patients had a median age of 69.05 years (range: 54 to 89 years). A smaller proportion of participants were female (median: 38%) (Table 1). AF and VTE prophylaxis were indications for DOAC use among 93.43% (N = 214,198) and 6.57% (N = 15,058) of patients, respectively, across the 31 studies. The assumption of transitivity was accepted because no variability was identified in the study and population baselines (Supplementary Tables 3, 4). Supplementary Figure 1 summarizes the detailed risk of bias assessments.
Figure 2. Network diagrams of comparisons of different treatment outcomes in patients receiving warfarin and different direct oral anticoagulants (DOACs). (A) Comparisons of all fracture risks in patients receiving warfarin and DOACs. (B) Comparisons of spinal fracture risks in patients receiving warfarin and DOACs. (C) Comparisons of hip fracture risks in patients receiving warfarin and DOACs.
Table 1. Characteristics of included studies.
Methodological Quality and Assessment of Risk of Bias
The main sources of RoB in the included RCTs were blinding of participants, personnel, and incomplete outcome data. Connolly et al. (34), EINSTEIN Investigators et al. (38), and EINSTEIN–PE Investigators et al. (42) had a high risk of performance bias, while Gibson et al. (49), Hohnloser et al. (58), Piazza et al. (51) and Weitz et al. (37) had a risk of attrition bias (Supplementary Figures 1A,B). The quality of non-RCTs was fairly good (Supplementary Table 5). Most studies had funding from multinational pharmaceutical companies. Only Huang et al. (63) and Wang et al. (61) did not report external funding.
Fracture Risk
We summarized our random-effects network meta-analysis and pairwise comparison of fracture risks in Figure 3; Supplementary Table 6. We ranked the risk of any fractures of DOACs against warfarin and the SUCRA probability (Supplementary Figures 3, 4; Supplementary Table 7).
Figure 3. Pooled estimates of the network meta-analysis. Comparisons, column vs. row, should be read from left to right and are ordered relative to overall effectiveness. (A) Pooled risk ratios (95% confidence intervals [CI]) for all fractures. (B) Pooled risk ratios (95% CI) for spinal fractures. (C) Pooled risk ratios (95% CI) for hip fractures.
Any Fracture Risk
This outcome was reported by 31 studies with 455,343 participants. The overall structure is shown in Figure 2A. VKA users had 5,553 fractures (5,553/2,21,203, 2.51%), dabigatran users had 2,578 fractures (2,578/78,810, 3.27%), rivaroxaban users had 2,025 fractures (2,025/1,06,996, 1.89%), apixaban users had 666 fractures (666/35,359, 1.88%) and edoxaban users had 254 fractures (254/12,975, 1.96%). Comparing network estimates of fracture risk between DOACs and warfarin, apixaban users (RR: 0.59; 95% CI: 0.50 to 0.71) and rivaroxaban users (RR: 0.72; 95% CI: 0.64 to 0.84) showed a statistically significant reduction in fracture risk, compared to warfarin users. No significant fracture risk reduction was observed among edoxaban, dabigatran, and warfarin users (Figure 3A). In terms of any fracture risk, apixaban (SUCRA = 98.0%) was most likely to be ranked the best, followed by rivaroxaban (SUCRA = 71.6%) (Supplementary Figures 3A,B, 4A; Supplementary Table 7).
Spine Fracture Risk
This outcome was reported in 10 studies with 83,842 participants (34, 38, 40–46, 54). The overall structure is shown in Figure 2B. VKA users had 61 fractures (61/4,1849, 0.15%), dabigatran users had 5 fractures (5/6,839, 0.07%), rivaroxaban users had 15 fractures (15/11,920, 0.13%) and edoxaban users had 16 fractures (16/11,153, 0.14%). No spinal fracture event was reported among apixaban users. Pooled estimates revealed no significant differences among apixaban users (RR: 0.07; 95% CI: 0.01 to 0.57), rivaroxaban users (RR: 0.75; 95% CI: 0.34 to 1.69), edoxaban users (RR: 0.78; 95% CI: 0.33 to 1.81) and dabigatran users (RR: 1.72; 95% CI: 0.32 to 9.17), when compared to warfarin users (Figure 3B). Apixaban (SUCRA = 98.5%) was most likely to be ranked the best in terms of risks for spine fracture (Supplementary Figures 3A,B, 4B; Supplementary Table 7).
Hip Fracture Risk
This outcome was reported in 16 studies with 228,133 participants (34, 35, 38, 40–42, 44–49, 51, 54, 60, 61). The overall structure is shown in Figure 2C. VKA users had 453 fractures (453/1,02,133, 0.44%), dabigatran users had 195 fractures (195/42,434, 0.46%), rivaroxaban users had 158 fractures (158/44,631, 0.35%), apixaban users had 69 fractures (69/27,726, 0.25%), and edoxaban users had 36 fractures (36/11,209, 0.32%). Overall, apixaban users generated the lowest pooled fracture risk estimate (RR: 0.56; 95% CI: 0.43 to 0.74), followed by rivaroxaban users (RR: 0.73; 95% CI: 0.60 to 0.88), edoxaban users (RR: 0.73; 95% CI: 0.47 to 1.12) and dabigatran users (RR: 1.06; 95% CI: 0.89 to 1.26), compared to warfarin users (Figure 3C). Apixaban (SUCRA = 90.4%) was most likely to be ranked the best in terms of risks for hip fracture (Supplementary Figures 3A,B, 4C; Supplementary Table 7).
Subgroup Analyses
The detailed results of subgroup analyses were presented in Supplementary Tables 8–10. Of note, in 22 studies with the indication of NVAF, apixaban reported the lowest fracture risk compared to warfarin (RR: 0.59; 95% CI: 0.58 to 0.75), followed by rivaroxaban (RR: 0.70; 95% CI: 0.58 to 0.85, Supplementary Table 10). Eight studies reported the indication of VTE, none of the DOACs were statistically significant in fracture reduction compared to warfarin. The results were similar in the subgroup of patients older than 65 and male predominant studies. Advanced age and male sex are both common characteristics of the NVAF population, with both subgroups concluding the lowest fracture risk in apixaban users (Supplementary Table 10).
Exploration for Inconsistency and Publication Bias
We found no evidence of global inconsistency in any of the outcomes using the design-by-treatment interaction models (Supplementary Table 11A). Furthermore, no substantial inconsistency between direct and indirect comparisons was observed in the side-splitting models (Supplementary Table 11B). Supplementary Figure 4A shows the comparison-adjusted funnel plots of fracture risks in the included studies, which revealed no significant funnel plot asymmetry. Lastly, the Egger test revealed no evidence of small-study bias (Supplementary Figure 4A).
Grade
We incorporated the GRADE judgments for network estimates of fracture risks. The certainty of evidence for the risk between anticoagulants varied; it was moderate for most of the comparisons involving DOACs against warfarin with regards to risks for any fracture, spine fracture, and hip fracture. The certainty of evidence was mostly moderate to low for the comparisons between different DOACs (Supplementary Tables 12, 13A–C).
Discussion
The current study aimed to identify fracture risks among patients prescribed DOACs and warfarin. The principal findings of this study were that patients who were prescribed apixaban carried the lowest fracture risk, followed by rivaroxaban, edoxaban and dabigatran, compared to patients prescribed warfarin. When assessing future fracture risk, it is crucial to consider both patient medication and medical history, given that 30% of patients presenting with a proximal femoral fracture receive anticoagulation therapy (64). Owens et al. reported that dabigatran, rivaroxaban, and apixaban might be used safely in NVAF patients with specific valvular heart diseases including aortic stenosis, aortic regurgitation, and mitral regurgitation (65). By contrast, patients with moderate to severe mitral stenosis or mechanical valves should continue to receive warfarin, as these patients have routinely been excluded from NVAF clinical trials (65). Furthermore, NVAF or VTE patients may require long-term anticoagulation therapy. Previous studies have reported that long-term exposure to VKAs is associated with an increased risk of fractures (66). These findings could be an important reference for clinicians when evaluation of fracture risk is necessary for patients at high risk of fractures, such as the elderly, who need to be on anticoagulation for NVAF.
As regards DOAC use and fracture risk, the literature remains conflicted. Both Lau et al. and Lutsey et al. report that DOACs carry a lower risk of fractures in patients with NVAF in the US and Hong-Kong, respectively, compared to warfarin (21, 60). However, Lucenteforte et al. find no differences in fracture risk between DOACs and VKA in patients with NVAF in Italy (55). These discrepancies in results may be attributed to heterogeneity of the study populations and the studies' power to detect event rate differences, whereby new RCTs and cohort studies that have since appeared add to our understanding of DOACs, especially those new to us (apixaban). Our systematic review and network meta-analysis evaluated 24 RCTs and seven cohort studies and observed that DOAC use was associated with a 21% risk reduction in reported fractures, compared to patients receiving warfarin.
Physiologically, the difference in fracture risk between DOACs and VKAs may be attributed to pharmacologic bone mineral density. Extensive literature survey reveals that both hip and vertebral fractures are most common among osteoporotic patients (67). VKAs inhibit the carboxylation of vitamin K-dependent bone mineralization proteins, including osteocalcin, matrix Gla protein, and periostin, increasing fracture risk (10, 11, 14, 68). Inhibition of osteocalcin carboxylation reduces adherence to calcium and hydroxyapatite, decreasing bone mineral density (BMD) and increasing the risk of osteoporosis (69). In animal studies, Fusaro et al. determined that among rats administered warfarin, a significant decrease in histomorphometric bone volume and increase in trabecular separation was observed, compared to both Dabigatran and placebo groups (11). In human studies, Rezaieyazdi et al. observed a marked reduction in BMD (g/cm2) and T-score of the lumbar spine among 70 rheumatic valvular heart disease patients taking warfarin, compared to controls (P = 0.048) (12). Warfarin use was the only risk factor of significant importance on spinal T-score (P < 0.03) (12). These findings support the utility of DOACs in decreasing fractures, compared to VKAs.
Kuo et al. (22) queried the Taiwan National Health Insurance database and reported that among 56,795 patients prescribed DOACs, dabigatran users show a lower incidence of osteoporotic fracture and spine fracture than patients receiving standard-dose rivaroxaban and apixaban. Our findings regarding the lower fracture risk of DOACs compared to warfarin have supported the already favorable clinical efficacy and side effect profiles of DOACs, compared to warfarin. Apixaban is superior to warfarin in the prevention of stroke and systemic embolism. The rates of stroke and ICH are both significantly lower in the ARISTOTLE trial (40). In a meta-analysis of 28 RCTs comparing DOACs with warfarin all DOACs have a higher rate of major GI bleeding, except apixaban (70). Our study findings also showed a statistically significant, lowered risk of fracture for apixaban, compared to warfarin. This lends support to the safety of apixaban use in elderly patients with regard to GI bleeding profiles, especially if these patients are at high risk of fracture.
Although DOACs have been reported to decrease fracture risk with protective bone mineralization properties compared to VKAs, not all fractures pose the same risk; therefore, subgroup analysis of anatomic fracture location is critical. Concerning hip fractures, our study determined that compared to warfarin, all of the DOACs except dabigatran exhibited a decreased hip fracture risk in the following descending order: apixaban, rivaroxaban, and edoxaban. Consistent with our findings, Huang et al. report a statistically significant risk reduction in hip fractures among adult users of DOACs, compared to VKAs, with varying risk reduction rates among the DOACs (63). Further research is required to determine the pharmacological mechanism of apixaban that contributes to fracture risk reduction in comparison to other DOACs.
Unlike osteoporotic hip and spine fractures among the elderly, trauma is typically associated with a high energy mechanism in the younger population with fewer comorbidities that do not require anticoagulation. Our findings showed that when only patients below 65 were included, no significant effect was seen among the DOACs, compared to warfarin. Most current literature focuses on patients with pre-existing comorbidities requiring anticoagulation treatment (71, 72). Second, most NVAF or VTE patients may require long-term anticoagulation therapy, and previous studies have indicated that long-term exposure to VKAs is associated with an increased risk of fractures (69, 73). In our study, we found that long-term DOAC exposure of at least one year also decreased the risk of fractures by 21%, compared to warfarin. Although the fracture types, treatment duration, and patients' sex or age varied among the included studies, the resulting overall robustness was proven by the subgroup analyses. The older female population is already known to be associated with increased fracture risk (74). We found that the female- and male population achieved similar effects when using DOACs to decrease fracture risks. Apixaban had the lowest fracture risks (RR: 0.55; 95% CI: 0.46 to 0.65), compared to warfarin, in the predominantly male studies. When we included only studies with younger patients (aged < 65), no significant effect was seen among the DOACs, compared to warfarin, which may be explained by the diminished overall sample size as a result of including only these studies. Most studies evaluated patients older than 65.
To the best of our knowledge, at present, this is the most comprehensive, up-to-date network meta-analysis to analyze the fracture risk among patients receiving DOACs and VKA. However, some limitations must be addressed. It should be noted that there was cohort heterogeneity among the studies. Although over 90% of studies analyzed involved AF patients, some study cohorts included trauma patients who received oral anticoagulation for thromboprophylaxis. Reassuringly, our sensitivity analysis revealed consistent results in patients with varying indications for anticoagulation. It should also be noted that potential confounders, including age, sex, race, and comorbidities, were adjusted for, using propensity score matching to allow for robust, accurate data comparison. Additionally, most studies did not provide BMD data as it rarely was a primary or secondary outcome; therefore, further research is required to quantify the association of oral anticoagulants with measured changes in T-score. Future meta-analyses on individual-level participant data and head-to-head prospective studies will be beneficial to confirm the findings above.
Conclusion
In summary, this network meta-analysis demonstrated that apixaban had the lowest pooled fracture risk, compared to other DOACs, and that the four major DOACs had lower fracture risk than warfarin. Similar results were found in sensitivity analyses with lower heterogeneity and inconsistency. These findings might benefit clinical practice for the individualized use of anticoagulants; however, future, large head-to-head prospective studies are required to validate these findings.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author Contributions
ST, C-WH, S-CS, and L-TK conceived and designed the study, interpreted data, and contributed to the final version of this report. ST and C-WH selected the articles and extracted the data. ST and L-TK analyzed the data. ST, C-WH, and L-TK wrote the draft. ET, OO, DV, WC, and J-RH made critical revisions. All authors agreed with the results and conclusions reported. All authors contributed to the article and approved the submitted version.
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.
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/fcvm.2022.896952/full#supplementary-material
Abbreviations
CI, confidence interval; DOAC, direct oral anticoagulants; NVAF, non-valvular atrial fibrillation; PICO, Patient-Intervention-Comparison-Outcome; PRISMA, Preferred Reporting Items for Systematic Reviews and Meta-Analyses; RR, relative risk; VKA, vitamin K antagonists.
References
1. Osteoporosis. (2021). Available online at: https://www.cdc.gov/nchs/fastats/osteoporosis.htm (accessed March 15, 2022).
2. Schousboe JT. Mortality after osteoporotic fractures: what proportion is caused by fracture and is preventable? J Bone Miner Res. (2017) 32:1783–8. doi: 10.1002/jbmr.3216
3. Mou L, Norby FL, Chen LY, O'Neal WT, Lewis TT, Loehr LR, et al. Lifetime risk of atrial fibrillation by race and socioeconomic status: ARIC study (Atherosclerosis Risk in Communities). Circ Arrhythm Electrophysiol. (2018) 11:e006350. doi: 10.1161/CIRCEP.118.006350
4. Chen A, Stecker EA, Warden B. Direct oral anticoagulant use: a practical guide to common clinical challenges. J Am Heart Assoc. (2020) 9:e017559. doi: 10.1161/JAHA.120.017559
5. Gage BF. Warfarin-induced fractures in atrial fibrillation? J Am Coll Cardiol. (2019) 74:2159–61. doi: 10.1016/j.jacc.2019.08.1026
6. Binding C, Bjerring Olesen J, Abrahamsen B, Staerk L, Gislason G, Nissen Bonde A. Osteoporotic fractures in patients with atrial fibrillation treated with conventional versus direct anticoagulants. J Am Coll Cardiol. (2019) 74:2150–8. doi: 10.1016/j.jacc.2019.08.1025
7. Aryannejad A, Gandominejad A, Tabary M, Noroozi N, Abbasi A, Araghi F, et al. Protective effect of modafinil on skin flap survival in the experimental random-pattern skin flap model in rats: the role of ATP-sensitive potassium channels and nitric oxide pathway. J Plast Reconst Aesthet Surg. (2021) 74:1346–54. doi: 10.1016/j.bjps.2020.10.084
8. Sugiyama T. Osteoporotic fractures associated with dabigatran vs warfarin. JAMA. (2017) 318:90–1. doi: 10.1001/jama.2017.6908
9. Sugiyama T. An update on hip fracture risk associated with anticoagulant therapy: warfarin versus direct oral anticoagulants. Expert Opin Drug Saf. (2020) 19:1219–20. doi: 10.1080/14740338.2020.1811230
10. Pietschmann P, Woloszczuk W, Panzer S, Kyrle P, Smolen J. Decreased serum osteocalcin levels in phenprocoumon-treated patients. J Clin Endocrinol Metab. (1988) 66:1071–4. doi: 10.1210/jcem-66-5-1071
11. Fusaro M, Dalle Carbonare L, Dusso A, Arcidiacono MV, Valenti MT, Aghi A, et al. Differential effects of dabigatran and warfarin on bone volume and structure in rats with normal renal function. PLoS ONE. (2015) 10:e0133847. doi: 10.1371/journal.pone.0133847
12. Rezaieyazdi Z, Falsoleiman H, Khajehdaluee M, Saghafi M, Mokhtari-Amirmajdi E. Reduced bone density in patients on long-term warfarin. Int J Rheum Dis. (2009) 12:130–5. doi: 10.1111/j.1756-185X.2009.01395.x
13. Sugiyama T, Kugimiya F, Kono S, Kim YT, Oda H. Warfarin use and fracture risk: an evidence-based mechanistic insight. Osteoporos Int. (2015) 26:1231–2. doi: 10.1007/s00198-014-2912-1
14. Sugiyama T, Takaki T, Sakanaka K, Sadamaru H, Mori K, Kato Y, et al. Warfarin-induced impairment of cortical bone material quality and compensatory adaptation of cortical bone structure to mechanical stimuli. J Endocrinol. (2007) 194:213–22. doi: 10.1677/JOE-07-0119
15. Tufano A, Coppola A, Cerbone AM, Ruosi C, Franchini M. Preventing postsurgical venous thromboembolism: pharmacological approaches. Semin Thromb Hemost. (2011) 37:252–66. doi: 10.1055/s-0031-1273089
16. January CT, Wann LS, Alpert JS, Calkins H, Cigarroa JE, Cleveland JC, et al. 2014 AHA/ACC/HRS guideline for the management of patients with atrial fibrillation. J Am Coll Cardiol. (2014) 64:e1–e76. doi: 10.1016/j.jacc.2014.03.021
17. Kirchhof P, Benussi S, Kotecha D, Ahlsson A, Atar D, Casadei B, et al. 2016 ESC Guidelines for the management of atrial fibrillation developed in collaboration with EACTS. Kardiol Pol. (2016) 74:1359–469. doi: 10.5603/KP.2016.0172
18. January CT, Samuel Wann L, Calkins H, Chen LY, Cigarroa JE, Cleveland JC, et al. 2019 AHA/ACC/HRS Focused update of the 2014 AHA/ACC/HRS guideline for the management of patients with atrial fibrillation: a report of the American College of Cardiology/American Heart Association Task Force on Clinical Practice Guidelines and the Heart Rhythm Society in collaboration with the society of thoracic surgeons. Circulation. (2019) 140:e125–51. doi: 10.1161/CIR.0000000000000665
19. Tufano A, Contaldi P, Franchini M, Minno G, Coppola A. Oral anticoagulant drugs and the risk of osteoporosis: new anticoagulants better than old? Semin Thromb Hemost. (2015) 41:382–8. doi: 10.1055/s-0034-1543999
20. Huang H-K, Liu PP-S, Hsu J-Y, Lin S-M, Peng CC-H, Wang J-H, et al. Risk of osteoporosis in patients with atrial fibrillation using non-vitamin K antagonist oral anticoagulants or warfarin. J Am Heart Assoc. (2020) 9:e013845. doi: 10.1161/JAHA.119.013845
21. Lau WCY, Cheung C-L, Man KKC, Chan EW, Sing CW, Lip GYH, et al. Association between treatment with apixaban, dabigatran, rivaroxaban, or warfarin and risk for osteoporotic fractures among patients with atrial fibrillation: a population-based cohort study. Ann Intern Med. (2020) 173:1–9. doi: 10.7326/M19-3671
22. Kuo L-T, Lin S-J, Wu VC-C, Chang J-J, Chu P-H, Lin Y-S. Direct oral anticoagulants and the risk of osteoporotic fractures in patients with non-valvular atrial fibrillation. Ther Adv Musculoskelet Dis. (2021) 13:1759720X211011374. doi: 10.1177/1759720X211011374
23. Higgins JPT, Altman DG, Gotzsche PC, Juni P, Moher D, Oxman AD, et al. The Cochrane Collaboration's tool for assessing risk of bias in randomised trials. BMJ. (2011) 343:d5928–d5928. doi: 10.1136/bmj.d5928
24. Wells G, Wells G, Shea B, Shea B, O'Connell D, Peterson J. The Newcastle-Ottawa Scale (NOS) for Assessing the Quality of Nonrandomised Studies in Meta-Analyses. (2014). Available online at: https://www.semanticscholar.org/paper/c293fb316b6176154c3fdbb8340a107d9c8c82bf (accessed June 30, 2021).
25. Rosli RRM, Norhayati MN, Ismail SB. Effectiveness of iron polymaltose complex in treatment and prevention of iron deficiency anemia in children: a systematic review and meta-analysis. PeerJ. (2021) 9:e10527. doi: 10.7717/peerj.10527
26. DerSimonian R, Laird N. Meta-analysis in clinical trials. Control Clin Trials. (1986) 7:177–88. doi: 10.1016/0197-2456(86)90046-2
27. Salanti G, Ades AE, Ioannidis JPA. Graphical methods and numerical summaries for presenting results from multiple-treatment meta-analysis: an overview and tutorial. J Clin Epidemiol. (2011) 64:163–71. doi: 10.1016/j.jclinepi.2010.03.016
28. White IR, Barrett JK, Jackson D, Higgins JPT. Consistency and inconsistency in network meta-analysis: model estimation using multivariate meta-regression. Res Synth Methods. (2012) 3:111–25. doi: 10.1002/jrsm.1045
29. Dias S, Welton NJ, Caldwell DM, Ades AE. Checking consistency in mixed treatment comparison meta-analysis. Stat Med. (2010) 29:932–44. doi: 10.1002/sim.3767
30. Egger M, Davey Smith G, Schneider M, Minder C. Bias in meta-analysis detected by a simple, graphical test. BMJ. (1997) 315:629–34. doi: 10.1136/bmj.315.7109.629
31. Sterne JA, Egger M. Funnel plots for detecting bias in meta-analysis: guidelines on choice of axis. J Clin Epidemiol. (2001) 54:1046–55. doi: 10.1016/S0895-4356(01)00377-8
32. Sterne JA, Gavaghan D, Egger M. Publication and related bias in meta-analysis: power of statistical tests and prevalence in the literature. J Clin Epidemiol. (2000) 53:1119–29. doi: 10.1016/S0895-4356(00)00242-0
33. Ezekowitz MD, Reilly PA, Nehmiz G, Simmers TA, Nagarakanti R, Parcham-Azad K, et al. Dabigatran with or without concomitant aspirin compared with warfarin alone in patients with nonvalvular atrial fibrillation (PETRO Study). Am J Cardiol. (2007) 100:1419–26. doi: 10.1016/j.amjcard.2007.06.034
34. Connolly SJ, Ezekowitz MD, Yusuf S, Eikelboom J, Oldgren J, Parekh A, et al. Dabigatran versus warfarin in patients with atrial fibrillation. N Engl J Med. (2009) 361:1139–51. doi: 10.1056/NEJMoa0905561
35. Schulman S, Kearon C, Kakkar AK, Mismetti P, Schellong S, Eriksson H, et al. Dabigatran versus warfarin in the treatment of acute venous thromboembolism. N Engl J Med. (2009) 361:2342–52. doi: 10.1056/NEJMoa0906598
36. U.S. National Library of Medicine. A Dose Response Study of Dabigatran Etexilate(BIBR 1048) in Pharmacodynamics and Safety in Patients With Non-valvular Atrial Fibrillation in Comparison to Warfarin. (2010). Available online at: https://clinicaltrials.gov/ct2/show/NCT01136408?cond=Nct01136408&draw=2&rank=1 (accessed April 10, 2022).
37. Weitz JI, Connolly SJ, Patel I, Salazar D, Rohatagi S, Mendell J, et al. Randomised, parallel-group, multicentre, multinational phase 2 study comparing edoxaban, an oral factor Xa inhibitor, with warfarin for stroke prevention in patients with atrial fibrillation. Thromb Haemost. (2010) 104:633–41. doi: 10.1160/TH10-01-0066
38. EINSTEIN Investigators, Bauersachs R, Berkowitz SD, Brenner B, Buller HR, Decousus H, et al. Oral rivaroxaban for symptomatic venous thromboembolism. N Engl J Med. (2010) 363:2499–510. doi: 10.1056/NEJMoa1007903
39. Chung N, Jeon H-K, Lien L-M, Lai W-T, Tse H-F, Chung W-S, et al. Safety of edoxaban, an oral factor Xa inhibitor, in Asian patients with non-valvular atrial fibrillation. Thromb Haemost. (2011) 105:535–44. doi: 10.1160/TH10-07-0451
40. Granger CB, Alexander JH, McMurray JJV, Lopes RD, Hylek EM, Hanna M, et al. Apixaban versus warfarin in patients with atrial fibrillation. N Engl J Med. (2011) 365:981–92. doi: 10.1056/NEJMoa1107039
41. Patel MR, Mahaffey KW, Garg J, Pan G, Singer DE, Hacke W, et al. Rivaroxaban versus warfarin in nonvalvular atrial fibrillation. N Engl J Med. (2011) 365:883–91. doi: 10.1056/NEJMoa1009638
42. EINSTEIN–PE Investigators, Büller HR, Prins MH, Lensin AW, Decousus H, Jacobson BF, et al. Oral rivaroxaban for the treatment of symptomatic pulmonary embolism. N Engl J Med. (2012) 366:1287–97. doi: 10.1056/NEJMoa1113572
43. Hori M, Matsumoto M, Tanahashi N, Momomura S-I, Uchiyama S, Goto S, et al. Rivaroxaban vs. warfarin in Japanese patients with atrial fibrillation. Circulation. (2012) 76:2104–11. doi: 10.1253/circj.CJ-12-0454
44. Hokusai-VTE Investigators, Büller HR, Décousus H, Grosso MA, Mercuri M, Middeldorp S, et al. Edoxaban versus warfarin for the treatment of symptomatic venous thromboembolism. N Engl J Med. (2013) 369:1406–15. doi: 10.1056/NEJMoa1306638
45. Agnelli G, Buller HR, Cohen A, Curto M, Gallus AS, Johnson M, et al. Oral apixaban for the treatment of acute venous thromboembolism. N Engl J Med. (2013) 369:799–808. doi: 10.1056/NEJMoa1302507
46. Giugliano RP, Ruff CT, Braunwald E, Murphy SA, Wiviott SD, Halperin JL, et al. Edoxaban versus warfarin in patients with atrial fibrillation. N Engl J Med. (2013) 369:2093–104. doi: 10.1056/NEJMoa1310907
47. Schulman S, Kearon C, Kakkar AK, Schellong S, Eriksson H, Baanstra D, et al. Extended use of dabigatran, warfarin, or placebo in venous thromboembolism. N Engl J Med. (2013) 368:709–18. doi: 10.1056/NEJMoa1113697
48. Schulman S, Kakkar AK, Goldhaber SZ, Schellong S, Eriksson H, Mismetti P, et al. Treatment of acute venous thromboembolism with dabigatran or warfarin and pooled analysis. Circulation. (2014) 129:764–72. doi: 10.1161/CIRCULATIONAHA.113.004450
49. Gibson CM, Mehran R, Bode C, Halperin J, Verheugt FW, Wildgoose P, et al. Prevention of bleeding in patients with atrial fibrillation undergoing PCI. N Engl J Med. (2016) 375:2423–34. doi: 10.1056/NEJMoa1611594
50. Goette A, Merino JL, Ezekowitz MD, Zamoryakhin D, Melino M, Jin J, et al. Edoxaban versus enoxaparin–warfarin in patients undergoing cardioversion of atrial fibrillation (ENSURE-AF): a randomised, open-label, phase 3b trial. Lancet. (2016) 388:1995–2003. doi: 10.1016/S0140-6736(16)31474-X
51. Piazza G, Mani V, Goldhaber SZ, Grosso MA, Mercuri M, Lanz HJ, et al. Magnetic resonance venography to assess thrombus resolution with edoxaban monotherapy versus parenteral anticoagulation/warfarin for symptomatic deep vein thrombosis: a multicenter feasibility study. Vasc Med. (2016) 21:361–8. doi: 10.1177/1358863X16645853
52. Bengtson LGS, Lutsey PL, Chen LY, MacLehose RF, Alonso A. Comparative effectiveness of dabigatran and rivaroxaban versus warfarin for the treatment of non-valvular atrial fibrillation. J Cardiol. (2017) 69:868–76. doi: 10.1016/j.jjcc.2016.08.010
53. Calkins H, Willems S, Gerstenfeld EP, Verma A, Schilling R, Hohnloser SH, et al. Uninterrupted dabigatran versus warfarin for ablation in atrial fibrillation. N Engl J Med. (2017) 376:1627–36. doi: 10.1056/NEJMoa1701005
54. Cannon CP, Bhatt DL, Oldgren J, Lip GYH, Ellis SG, Kimura T, et al. Dual Antithrombotic therapy with dabigatran after PCI in atrial fibrillation. N Engl J Med. (2017) 377:1513–24. doi: 10.1056/NEJMoa1708454
55. Lucenteforte E, Bettiol A, Lombardi N, Mugelli A, Vannacci A. Risk of bone fractures among users of oral anticoagulants: an administrative database cohort study. Eur J Intern Med. (2017) 44:e30–1. doi: 10.1016/j.ejim.2017.07.022
56. Norby FL, Bengtson LGS, Lutsey PL, Chen LY, MacLehose RF, Chamberlain AM, et al. Comparative effectiveness of rivaroxaban versus warfarin or dabigatran for the treatment of patients with non-valvular atrial fibrillation. BMC Cardiovasc Disord. (2017) 17:238. doi: 10.1186/s12872-017-0672-5
57. Ezekowitz MD, Pollack CV Jr, Halperin JL, England RD, VanPelt Nguyen S, Spahr J, et al. Apixaban compared to heparin/vitamin K antagonist in patients with atrial fibrillation scheduled for cardioversion: the EMANATE trial. Eur Heart J. (2018) 39:2959–71. doi: 10.1093/eurheartj/ehy148
58. Hohnloser SH, Camm J, Cappato R, Diener H-C, Heidbüchel H, Mont L, et al. Uninterrupted edoxaban vs. vitamin K antagonists for ablation of atrial fibrillation: the ELIMINATE-AF trial. Eur Heart J. (2019) 40:3013–21. doi: 10.1093/eurheartj/ehz190
59. Ferro JM, Coutinho JM, Dentali F, Kobayashi A, Alasheev A, Canhão P, et al. Safety and efficacy of dabigatran etexilate vs dose-adjusted warfarin in patients with cerebral venous thrombosis: a randomized clinical trial. JAMA Neurol. (2019) 76:1457–65. doi: 10.1001/jamaneurol.2019.2764
60. Lutsey PL, Norby FL, Ensrud KE, MacLehose RF, Diem SJ, Chen LY, et al. Association of anticoagulant therapy with risk of fracture among patients with atrial fibrillation. JAMA Intern Med. (2020) 180:245–53. doi: 10.1001/jamainternmed.2019.5679
61. Wang S, Qiang R, Zhang H. Rivaroxaban versus warfarin in Chinese elder diabetic patients with non-valvular atrial fibrillation: a retrospective study on the effectiveness and associated risks (EAR) analysis. Trop J Pharm Res. (2020) 19:893–8. doi: 10.4314/tjpr.v19i4.31
62. Halperin JL, Levine GN, Al-Khatib SM, Birtcher KK, Bozkurt B, Brindis RG, et al. Further evolution of the ACC/AHA clinical practice guideline recommendation classification system: a report of the American college of cardiology/American heart association task force on clinical practice guidelines. Circulation. (2016) 133:1426–8. doi: 10.1161/CIR.0000000000000312
63. Huang H-K, Liu PP-S, Hsu J-Y, Lin S-M, Peng CC-H, Wang J-H, et al. Fracture risks among patients with atrial fibrillation receiving different oral anticoagulants: a real-world nationwide cohort study. Eur Heart J. (2020) 41:1100–8. doi: 10.1093/eurheartj/ehz952
64. Papachristos IV, Giannoudis PV. Proximal femur fractures in patients taking anticoagulants. EFORT Open Rev. (2020) 5:699–706. doi: 10.1302/2058-5241.5.190071
65. Owens RE, Kabra R, Oliphant CS. Direct oral anticoagulant use in nonvalvular atrial fibrillation with valvular heart disease: a systematic review. Clin Cardiol. (2017) 40:407–12. doi: 10.1002/clc.22659
66. Veronese N, Bano G, Bertozzo G, Granziera S, Solmi M, Manzato E, et al. Vitamin K antagonists' use and fracture risk: results from a systematic review and meta-analysis. J Thromb Haemost. (2015) 13:1665–75. doi: 10.1111/jth.13052
67. Wright NC, Looker AC, Saag KG, Curtis JR, Delzell ES, Randall S, et al. The recent prevalence of osteoporosis and low bone mass in the United States based on bone mineral density at the femoral neck or lumbar spine. J Bone Miner Res. (2014) 29:2520–6. doi: 10.1002/jbmr.2269
68. Morishima Y, Kamisato C, Honda Y, Furugohri T, Shibano T. The effects of warfarin and edoxaban, an oral direct factor Xa inhibitor, on gammacarboxylated (Gla-osteocalcin) and undercarboxylated osteocalcin (uc-osteocalcin) in rats. Thromb Res. (2013) 131:59–63. doi: 10.1016/j.thromres.2012.08.304
69. Dadwal G, Schulte-Huxel T, Kolb G. Effect of antithrombotic drugs on bone health. Z Gerontol Geriatr. (2020) 53:457–62. doi: 10.1007/s00391-019-01590-8
70. Radadiya D, Devani K, Brahmbhatt B, Reddy C. Major gastrointestinal bleeding risk with direct oral anticoagulants: does type and dose matter? – A systematic review and network meta-analysis. Eur J Gastroenterol Hepatol. (2021) 33:e50–8. doi: 10.1097/MEG.0000000000002035
71. Hamada SR, Espina C, Guedj T, Buaron R, Harrois A, Figueiredo S, et al. High level of venous thromboembolism in critically ill trauma patients despite early and well-driven thromboprophylaxis protocol. Ann Intensive Care. (2017) 7:97. doi: 10.1186/s13613-017-0315-0
72. Hamidi M, Zeeshan M, Sakran JV, Kulvatunyou N, O'Keeffe T, Northcutt A, et al. Direct oral anticoagulants vs low-molecular–weight heparin for thromboprophylaxis in nonoperative pelvic fractures. J Am Coll Surg. (2019) 228:89–97. doi: 10.1016/j.jamcollsurg.2018.09.023
73. Caraballo PJ, Heit JA, Atkinson EJ, Silverstein MD, O'Fallon WM, Castro MR, et al. Long-term use of oral anticoagulants and the risk of fracture. Arch Intern Med. (1999) 159:1750–6. doi: 10.1001/archinte.159.15.1750
Keywords: non-vitamin K antagonist oral anticoagulants, fracture, network meta-analysis, direct-acting oral anticoagulant (DOAC), warfarin, osteoporosis, atrial fibrillation, venous thromboembolism
Citation: Tsai SHL, Hu C-W, Shao S-C, Tischler EH, Obisesan OH, Vervoort D, Chen WC, Hu J-R and Kuo L-T (2022) Comparative Risks of Fracture Among Direct Oral Anticoagulants and Warfarin: A Systematic Review and Network Meta-Analysis. Front. Cardiovasc. Med. 9:896952. doi: 10.3389/fcvm.2022.896952
Received: 15 March 2022; Accepted: 21 April 2022;
Published: 23 May 2022.
Edited by:
Michael Wu, Warren Alpert Medical School of Brown University, United StatesReviewed by:
Daniela Poli, Careggi University Hospital, ItalyPierre Sabouret, Sorbonne University, France
Copyright © 2022 Tsai, Hu, Shao, Tischler, Obisesan, Vervoort, Chen, Hu and Kuo. 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: Liang-Tseng Kuo, light71829@gmail.com
†These authors have contributed equally to this work