- 1Department of Urology, Institute of Urology and National Clinical Research Center for Geriatrics, West China Hospital, Sichuan University, Chengdu, China
- 2Center of Biomedical Big Data, West China Hospital, Sichuan University, Chengdu, Sichuan, China
Background: Hypogonadism has become a major cause endangering men’s health and quality of life all over the world. Testosterone Therapy (TT) is a widely accepted treatment for relieving hypogonadal symptoms. However, the effect of different administrations of TT on prostate safety is still unclear.
Methods: We did a thorough search of PubMed, Embase and Cochrane Library to identify eligible studies up to January 2022. Randomized controlled trials (RCTs) and Cohort studies evaluating the impacts of using different formulations of TT on prostate parameters were included. Changes of prostate-specific antigen (PSA) level and prostate cancer (Pca) cases were used as the primary outcomes. Quality of individual studies was estimated by RoB2 (Cochrane tool for assessing the risk of bias in randomized trials) and the Newcastle-Ottawa scale (Tool for assessing non-RCTs). Certainty of evidence for each study was evaluated according to the evidence assessment criteria of the Oxford Evidence-based Medicine Center. Random-effect network meta-analysis(NMA)was performed based on the Bayesian model.
Results: Thirty-five studies (30 RCTs and 5 Cohort studies) with 7,740 participants were included. TT administration led to fewer Pca patients (RR=0.62, 95%CI [0.39,0.99], I2=0%), while little decreasing in PSA level (MD=-0.05, 95%CI [-0.08, -0.02], I2=0%). The NMA revealed that compared with other formulations, the intramuscular injection was the most likely to rank first in decreasing Pca cases. The TT also resulted in more biopsy cases (RR=2.38, 95%CI [1.01,5.60], I2=0%). As for NMA, intramuscular injection also performed relatively better in fewer prostate biopsy cases compared with transdermal group.
Conclusion: TT does not lead to abnormal PSA changes and increased risk of Pca in patients with hypogonadism or low testosterone level. Compared with other preparations of TT, intramuscular injection proved better in minimizing Pca cases and was more likely to result in fewer prostate biopsy cases.
1 Introduction
Male hypogonadism, a clinical syndrome caused by androgen deficiency, has become a major cause endangering men’s health and quality of life all over the world. It is roughly estimated that about a quarter of men are affected by androgen deficiency (1). Some studies have shown that testosterone levels in healthy people may also gradually decline with age (2, 3). Hypogonadism becomes more common with age grow and usually occurs in men over the age of 40, but not limited to this (4). Among middle-aged men, the incidence of hypogonadism ranges from 2.1% to 12.8% (5). It is reported that about half a million men are diagnosed with androgen deficiency every year in the United States, especially in people with certain diseases, such as obesity and diabetes (6). This common hypogonadism in middle-aged people is also known as late-onset hypogonadism (LOH). LOH usually induces many other diseases by reducing sex hormones, such as atherosclerosis, hyperglycemia, hypertension, and so on (7, 8).
Testosterone Therapy (TT) has become a widely accepted treatment to alleviate the symptoms of hypogonadism. However, there still remains many controversies about the safety of TT (9). One of the most controversial is the possible consequences for prostate safety risk (10). Traditional theories suggest that the presence of testosterone may affect or even promote the growth of prostate. Nevertheless, some studies have shown that testosterone treatment does not lead to histological changes in the prostate (11, 12). Marks et al (12) and Baillargeon et al (13) reported that TT does not increase the risk of Pca or lead to more aggressive prostate tumors. The EAU guidelines also pointed out that there is no evidence that testosterone therapy increases the risk of prostate cancer (4). Therefore, the association between TT and prostate safety risk remains unclear.
Previously, Pastuszak et al (14) conducted a study which compared the long-term effects of different testosterone preparations on men with hypogonadism. Inspired by their article, we try to explore the impact of TT on prostate safety risk from the perspective of different administrations of testosterone. Testosterone is currently available in gel, patch, intramuscular, oral and other preparations. They perform differently in the route of administration, pharmacokinetics and adverse events. Therefore, we can reasonably speculate that different administration methods of testosterone have different impacts on the prostate safety risk of TT. Nevertheless, there are few head-to-head RCTs to evaluate the prostate safety risk of different formulations, and network meta-analysis technology is needed to assess the comparison between preparations. In the absence of RCTs with direct pairwise comparison, network meta-analysis (NMA) is the most robust way to compare the interaction of different treatment groups in multiple compilation studies.
Therefore, the aim of this study is to make full use of the clinical data obtained from the included studies to conduct a NMA to assess the prostate safety risk of different administrations of TT and find out the hierarchical structure of them.
2 Methods
2.1 Search strategy
According to the Preferred Reporting Items for Systematic Reviews and Meta-analyses (PRISMA) reporting guideline and its extension for NMA, we conducted the systematic review and network meta-analysis (15, 16). Three main electronic databases (PubMed, Cochrane Library and Embase) were searched to find out potentially relevant research up to January 2022, with language restricted to English. We used the following MeSH terms for search: “Prostatic Neoplasms”, “Testosterone”, and “Hypogonadism”. The complete search strategy used for PubMed was shown in Supplementary Table 1. We also manually retrieved the reference lists of relevant studies and reviews. Two authors (BZ and XYX) independently reviewed the literature, and the inconsistencies were discussed and solved with the third author (SQ).
2.2 Selection criteria
The literature was reviewed by two separate authors independently (BZ and XYX). We evaluated the eligibility of studies using the population, intervention, comparator, outcome, and study (PICOS) method: (P) research involving population with hypogonadism or low testosterone level; (I) received different testosterone administrations (mentioned which administration was used); (C) group assigned with placebo or no TT was regarded as a comparator; (O) reporting one or more parameters of prostate safety; (S) RCTs or Cohort Studies. The selection process was shown in Figure 1 based on the PRISMA flowchart.
Figure 1 Preferred Reporting Items for Systematic Reviews and Meta-analyses flow diagram of study selection process.
2.3 Data extraction
The required clinical data were extracted by two separate reviewers from the included literature (BZ and XYX), including the following: author, year of publication, intervention measures, type of study, sample size, administration method, testosterone cut off level for study entry, exclusion criteria of prostate-specific antigen (PSA), median age, duration of study, dosage of the testosterone and the functional outcomes such as: Pca cases, PSA, prostate nodule, biopsy cases, abnormal PSA, International Prostate Symptom Score (IPSS) and prostate volume. Any inconsistencies were discussed and solved with the third reviewer (SQ).
2.4 Primary and secondary outcomes
As for the safety risk of prostate, we mainly evaluated it from two aspects: Pca and prostate growth. Referring to previous studies, we used the Pca cases, changes in PSA levels, prostate biopsy cases, prostate nodules and abnormal PSA to assess the risk of Pca (17). In terms of prostate growth, we used IPSS and prostate volume changes to evaluate. To ensure the reliability of the results, we used the two indicators which included most studies: the change of PSA level (27 studies) and Pca cases (19 studies) as the primary outcomes. And the rest indicators were used as the secondary outcomes in this study.
2.5 Study quality assessment
The risk of bias in each research was evaluated independently by two reviewers (BZ and XYX). According to the recommendations of the Cochrane manual, we selected RoB2, the most commonly used tool in randomized trials, as the tool to evaluate the risk of bias in RCTs (18). As for non-RCT studies, the Newcastle–Ottawa scale was used to assess the quality of included cohort studies (19). In this study, we believe that the total score ≧ 6 is high quality for cohort studies. Furthermore, consistent with previous studies, the certainty of evidence for each study was evaluated according to the evidence assessment criteria of the Oxford Evidence-based Medicine Center (20).
2.6 Consistency and sensitivity analysis
The NMA was based on the consistency assumption. We used the node-split method to verify the consistency between direct comparison and indirect comparison (21). The sensitivity analysis was conducted to examine the reliability of our findings. We excluded cohort studies which may be the source of heterogeneity in results and reanalyzed the clinical data without cohort studies. The results were compared with our previous findings to see whether the two results were consistent.
2.7 Data synthesis and statistical analyses
The continuous outcomes in our study were transformed into the same units of measurement and were pooled as mean difference (MD). And the relative risk (RR) was used to evaluate dichotomous variables. If the extracted data from the included studies were recorded in the median (interquartile range), we calculated the mean ± standard deviation (SD) according to the method described by Hozo et al. (22). We performed the pairwise meta-analysis using Review Manager v5.4 software to calculate the MD and RR of primary and secondary outcomes between different testosterone preparations and the control group. Given the possible potential heterogeneity, we used the Mantel-Haenszel (M-H) random-effects model for evaluation. The heterogeneity between included studies was estimated using chi-square (p < 0.05) and the I2 statistic (23). If I2 > 50%, we consider there may be moderate-to-high heterogeneity between studies (24).
The NMA was performed in a Bayesian hierarchical framework, and we used the random-effects model to evaluate the direct and indirect comparison between different testosterone formulations (25). In order to ensure the high level of evidence of NMA, we only used the data from RCTs for the NMA. And we used the Brooks-Gelman-Rubin method to evaluate the convergence (26). Intra-chain and inter-chain variances were compared by the potential scale reduction factor (PSRF), which was calculated in this method. If PSRF is close to 1, we consider that it has reached the approximate convergence (27). In the process of gradual convergence, there are four chains, the variance scaling factor was 2.5, simulation iterations were 100,000, tuning iterations were 20,000, and the thinning interval was 10. Rank probabilities were calculated for all different testosterone formulations and control group, and the surface under the cumulative ranking curve (SUCRA) is estimated by dividing the cumulative probability of all rankings by the number of rankings minus 1.
All statistical tests were 2-sided, and 95% confidence intervals (CIs) were reported. P < 0.05 was defined as statistically significant in this study. We used the “gemtc” “rjags” “meta” and “metafor” packages from R 4.1.2 (R project) (28) and Review Manager v5.4 software to conduct all the statistical analyses and forest plots.
3 Results
3.1 Study selection and network structure
We identified 716 articles for appraisal, including 49 studies which were searched from reference lists of relevant studies and reviews. After excluding the literature that does not meet the inclusion criteria, a total of 35 studies involving 7,740 participants were included. The detailed selection process was shown in Figure 1. In addition, the 35 included studies were consisted of 30 RCTs (11, 12, 29–56) and 5 cohort studies (14, 57–60).
When evaluating Pca cases, abnormal PSA and prostate nodule for NMA, we found that only one literature with a value of 0 was included in certain preparation, which may lead to abnormal results. In order to assure the robustness of the results, we excluded these unstable studies for evaluation in our NMA. The network structure for Pca cases and PSA level changes were displayed in Figure 2. The size of the circle represents the sample size of each arm, and the thickness of the line represents the number of head-to-head studies. If there is no connection between circles, it indicates that there are no research with a direct comparison between the two preparations at present.
Figure 2 Network structure of the comparisons of primary outcomes for the Bayesian network meta-analysis. (A) Pca Cases; (B) PSA Level Changes. The size of the circles represents the sample size of each arm (in parentheses), and the thickness of the line represents the number of head-to-head studies (beside the line).
3.2 Characteristics of included studies
The characteristics of the whole included studies were shown in Table 1. According to different definitions, the baseline testosterone levels of all included 7,740 participants were in the range of testosterone deficiency or normal low levels. The average age of all the participants ranged from 44 to 80 years old. The duration of participants’ follow-up ranged from 1 month to 12 years after the first administration. In this study, there were three ways of TT administration to compare the primary and secondary results, including oral, intramuscular and transdermal. Of the 35 articles included, 16 studies (11, 30–32, 36–46, 58) compared transdermal testosterone with the control group (placebo or No testosterone), 13 studies (12, 33–35, 47–52, 55–57) compared intramuscular testosterone with the control group, 3 studies (29, 53, 54) compared oral testosterone with the control group, and 3 studies (14, 59, 60) compared intramuscular testosterone with transdermal testosterone. Among the participants, 171 participants were included in the oral testosterone group, 1,404 participants in the intramuscular testosterone group, and 2,750 participants in the transdermal testosterone group.
3.3 Assessment of study quality
The bias risk assessment results of the studies were provided in Table 2 and Table 3. According to the recommendations of the Cochrane manual, we used the RoB2 tool to assess the risk of bias of all included RCTs. Among the 30 RCTs, six trials (11, 33, 35, 41, 49, 51) had some concerns about the bias in measurement of outcomes, most of which were caused by inappropriate measurement methods. In addition, there is a risk of missing outcome data bias in one literature (42), mainly because of the failure to assess the impact of missing data on the results. As for cohort studies, the Newcastle–Ottawa scale was used to evaluate the quality of included research. Of the five cohort studies, most were considered to be of high quality, and three studies (58–60) were regarded as low quality due to poor control of confounding factors and too short follow-up time. According to the evidence evaluation criteria of the Oxford evidence-based medicine center, we evaluated the evidence level of each included study, and the results were shown in Table 1.
3.4 Results of primary outcomes
3.4.1 Pca cases
Seventeen studies with 6,361 participants and three testosterone preparations, constituted the pairwise meta-analysis to evaluate the Pca cases. The pooled results suggested that the use of TT could lead to fewer Pca cases than the control group (RR=0.62, 95%CI [0.39,0.99], I2 = 0%). Comparing the effects of different testosterone preparations and the control group on the risk of Pca, we found that except for intramuscular testosterone group showed fewer cases of Pca (RR=0.54, 95%CI [0.32,0.90], I2 = 0%), the oral group (RR=0.19, 95%CI [0.01,4.01]) and transdermal group (RR=1.24, 95%CI [0.45,3.38], I2 = 0%) where there was no significant difference compared with the control group in Pca cases(Figure 3A).
There were 15 RCTs with 5,333 participants and two testosterone preparations included in the NMA. Since in the oral testosterone group there was only one study with a value of 0, which would lead to abnormal results in the evaluation. For the sake of robust results, we excluded the oral group for evaluation. The intramuscular injection group accounted for 12.5% of the data, the transdermal group accounted for 41.7%, and the remaining 45.8% were in the control group (Figure 2A). The NMA results were displayed in Table 4A. We used the staircase tables to evaluate the effect of intramuscular injection, transdermal and control groups on the Pca cases. Although there was no significant difference among intramuscular injection, transdermal and control groups, intramuscular injection could show the trend of leading fewer Pca cases. According to Bayesian model, we used SUCRA to rank the probability that each testosterone preparation ameliorated the risk of Pca. The results showed that the probability of intramuscular injection, transdermal and control groups were 18%, 84% and 48%, respectively. And owing to the reverse scaled outcome of Pca cases (where lower values indicate a better outcome), the intramuscular injection was the most likely to rank first in decreasing Pca cases (Supplementary Figure 1A).
3.4.2 PSA level changes
Twenty-four studies involving 4052 participants and three testosterone preparations, constituted the pairwise meta-analysis to assess changes in PSA levels. The pooled results suggested that the use of TT could result in little decrease in PSA level compared with the control group (MD=-0.05, 95%CI [-0.08, -0.02], I2 = 98%). Comparing the effects of different testosterone preparations and the control group on the changes of PSA level, we found that the transdermal group showed minute PSA increase ((MD=0.14, 95%CI [0.08,0.19], I2 = 88%), while the injection group (MD=0.04, 95%CI [-0.25,0.32], I2 = 95%) and the oral group (MD=-0.00, 95%CI [-0.00,0.00], I2 = 0%) showed that there was no significant statistical difference compared with the control group. (Figure 3B)
In the NMA, there were 22 RCTS with 3180 participants and three testosterone preparations included. The transdermal group accounted for 42.9% of the data, the injection group accounted for 8%, oral group accounted for 5.4% and the remaining 43.7% were in the control group (Figure 2B). The results of NMA were shown in the staircase Table 4B. In the absence of evidence of direct comparison between the oral, injection and transdermal groups, we used the NMA to produce an indirect comparison between them. The results indicated that the transdermal group showed more PSA changes than the control group (MD=0.22, 95%CI [0.10,0.34]) and there was no significant difference among the remaining groups. In the same way, we calculated the ranking probability of each group according to the Bayesian model. Similarity, due to the reverse scaled outcome of PSA changes, we found that the oral group had the highest probability of ranking first in ameliorating PSA level changes (Supplementary Figure 1B). The SUCRA indicated that the probability of oral, intramuscular injection, transdermal and control groups were 25%, 69%, 88% and 18%, respectively.
3.5 Results of secondary outcomes
A total of 29 studies evaluated secondary outcome indicators. Supplementary Figure 2 displayed the results of our pairwise meta-analysis for the efficacy of different administrations of testosterone compared with the control group.
In terms of Pca risk, 16 studies evaluated abnormal PSA, 10 studies assessed prostate biopsy, and 6 studies evaluated prostate nodules. Our results indicated that the prostate biopsy was the only one where there was statistical difference between TT and the control group (RR=2.38, 95%CI [1.01,5.60], I2 = 0%). The comparison of effects of different testosterone preparations and control group on prostate biopsy indicated that no significant difference was found in both transdermal group (RR=2.44, 95%CI [0.97,6.13], I2 = 0%) and intramuscular injection group (RR=2.00, 95%CI [0.19,20.61]) compared with the control group. Subsequently, in the absence of evidence of direct comparison between the injection and transdermal groups, we conducted the NMA of prostate biopsy to generate an indirect comparison between them. The results suggested that except for the comparison between the transdermal group and the control group (control vs transdermal, RR=0.16, 95%CI [0.03,0.62]), there was no significant difference between other groups (Supplementary Table 2). The control group had the highest probability of ranking first in fewer biopsy cases, and injection was more likely to perform better compared with transdermal group.
As for the risk of prostate growth, 12 studies evaluating IPSS and 10 studies evaluating prostate volume were included. In the pairwise meta-analysis, the pooled results suggested that the use of TT had no significant effect on IPSS (MD=-0.37, 95%CI [-3.45,2.71], I2 = 99%) and prostate volume (MD=3.35, 95%CI [-0.25,6.95], I2 = 95%) compared with control group.
3.6 Sensitivity analysis
We excluded cohort studies which may be the source of heterogeneity to perform the sensitivity analysis. After excluding all cohort studies, there were 30 RCTs remained. Since the studies evaluating prostate biopsy and prostate nodules are all RCTs, the results of sensitivity analysis were consistent with previous results. In pairwise meta-analysis, sensitivity analysis showed that the effect of TT on PSA level (MD=0.01, 95%CI [-0.00,0.02], I2 = 80%) was not statistically different from that of the control group, and the use of TT had no significant effect on Pca cases compared with the control group (RR=0.85, 95%CI [0.44,1.62], I2 = 0%), which was different from previous results. Comparing the effects of different testosterone preparations and control group on the changes of PSA level, we found that there was no statistical difference between oral and control group (MD=-0.00, 95%CI [-0.00,0.00], I2 = 0%), intramuscular injection and control group (MD=0.05, 95%CI [-0.12,0.21], I2 = 45%), except that the transdermal group showed little PSA increase (MD=0.22, 95%CI [0.07,0.37], I2 = 89%) (Supplementary Figure 3). As for the other outcome indicators, the results of sensitivity analysis were consistent with our previous results. In general, according to the results of sensitivity analysis, our analysis results are relatively robust.
4 Discussion
Up to now, the controversy about the safety of prostate treated with testosterone has not been resolved. EAU Guidelines on Sexual and Reproductive Health (2022) (4)states that testosterone treatment has historically raised concerns about the possibility of benign prostatic hyperplasia (BPH) associated with prostate growth, but a large number of studies have no evidence to prove any significant association between TT and abnormal prostate growth and Pca events. In addition, it is unclear how different testosterone administrations affect prostate safety, and which testosterone preparations are better. Therefore, we conducted an indirect and direct comparison of the effects of different testosterone preparations on prostate safety in a network meta-analysis. As far as we know, this study is the first NMA to assess the impact of different testosterone preparations on prostate safety risk.
Our findings suggested that the use of TT was associated with fewer Pca cases, little decreased PSA level and more prostate biopsy cases, while there was no significant difference between TT and the control group in other outcomes. However, considering that the unit of PSA is ng/ml, we believe that this small difference caused by TT is of no clinical significance, so we can hold that TT will not significantly affect the change of PSA level. No matter compared with the control group or other testosterone preparations, injection testosterone all had the greater potential to perform better in decreasing Pca cases. Although compared with other groups, the intramuscular injection group had the longest follow-up time (mean duration: 22.4 months) which indicates more likely to observe more Pca cases, this also indirectly shows the robustness of our results from another point of view. We believe that the difference between sensitivity analysis and our previous Pca results is also due to the same reason. The average follow-up time of included RCTs (mean duration:11.3 months) is even less than one year, which may lead to an unreliable conclusion that there is no significant relationship between TT and the incidence of Pca. As for the prostate biopsy, our results indicated that the control group had the highest probability to rank first in ameliorating prostate biopsy cases. Although TT may be associated with more biopsy cases, injection testosterone was relatively better than transdermal group in terms of fewer prostate biopsy cases. For why intramuscular injection seems to be more effective in protecting prostate, we speculate that it may be caused by the following reasons: 1. The pharmacokinetics of the three preparations are different. Compared with oral and transdermal testosterone administrations, intramuscular injection can maintain a stable high circulating testosterone level for a long time, so as to better protect the prostate (61). 2. Previous evidence showed that obesity is associated with an increased risk of prostate cancer (62). Intramuscular injection is more beneficial to muscles generation, which can more effectively reduce the risk of obesity and thus lower the incidence of prostate disease, compared with oral and transdermal testosterone (63).
Our results regarding the impact of TT on prostate safety are consistent with previously published studies on this topic. As for the effect of TT on prostate growth, according to traditional theory, prostate growth depends on the presence of androgens, so it has been suggested that TT may increase prostate volume (64). However, some meta-analysis studies have shown that TT has no significant effect on prostate volume and lower urinary tract symptoms (LUTS). In a meta-analysis performed by Guo et al (65), which involved 16 trials and 1,921 participants, the reviewers found that there was no significant change in IPSS levels in the participants who received the TT compared with the placebo group (MD=0.01, 95%CI [-0.37,0.39], P=0.96). Another meta-analysis involving 51 trials yielded similar results, showing no statistical differences in IPSS and lower urinary tract symptoms between participants treated with TT and the control group (66). This may be explained by Morgentaler and Traish’s “saturation hypothesis” (67), that is when androgen receptors in the human prostate are “saturated” with circulating androgens, the prostate does not grow in response to external androgen supplementation. Several large clinical trials have also validated our results (68, 69). Regarding the Pca risk of TT, some previous meta-analyses have also been consistent with our results. Kang et al (70) conducted a meta-analysis including 15 clinical trials to find out the correlation between TT and the risk of Pca. Their results indicated that using testosterone did not increase PSA levels, and TT was not significantly associated with the elevated PSA, which was in accordance with our results. Similar results have been found in some large clinical trials, which revealed that the incidence of Pca in the population receiving TT was always lower than the incidence rate reported in the general population (71, 72). This phenomenon of using testosterone may protect the prostate could be explained by the hypothesis proposed by xu et al (73). They believed that the significant decrease of testosterone level may play an important role in adverse prostate events, and stable testosterone level can significantly reduce the incidence of prostate disease. In addition, as mentioned above, testosterone can promote the formation of fat free tissue and reduce obesity rate, so as to protect the prostate.
Although no previous studies have conducted a network meta-analysis of the effects of various testosterone administrations on prostate safety yet, some systematic reviews of the currently available evidence for direct comparison between them have been conducted. Cui et al (74) conducted a systematic review of the effects of TT on prostate growth in 2013, and compared performance differences between different administrations. Their results suggested that TT did not increase the risk of prostate growth in either the short or long term, regardless of the administration method, which was consistent with our results. Additionally, Cui et al. also explored the effect of TT on Pca risk in 2014 (17). Their results of meta-analysis suggested that TT was safe in the short term and did not contribute to the development of Pca, regardless of administrations. Although the authors suggested that the study’s failure to detect an increased risk of Pca in TT group may be due to the fact that the included trials monitored the participants closely and withdrew those participants when their Pca risk was suspected. Nevertheless, we believe that this is inevitable, and we included more studies (our study vs Cui et al:23 vs 8) evaluating the impact of TT on Pca events, leading to different results compared with them.
A total of 30 RCTs and 5 cohort studies were included in this network meta-analysis. According to the quality assessment scale we used, most studies were of high quality. From a clinical perspective, the results of this study provide a new idea for further exploring the relationship between TT and prostate safety risk. The results of sensitivity analysis are basically consistent with those of our study, indicating that our results are robust and reliable.
Nevertheless, there are some limitations in our research. Firstly, there was some heterogeneity in testosterone dose, baseline PSA levels, and so on among the participants included in the study. For example, the dosage and timing of testosterone administration varied from study to study, which may affect the final outcomes of our analysis. Although subgroup and sensitivity analyses were conducted to assess the quality of the results, heterogeneity cannot be completely avoided. Secondly, most of the participants in our analysis were the elderly, and only in one study (11) the average age of the population was 44 years old. More research on middle-aged people should be included to make the analysis more representative. Thirdly, we have three included studies (44, 49, 51) with a total number of participants less than 20 in our analysis. Too few participants may affect the robustness of the analysis results. Finally, although studies included in our current analysis were followed for a relatively longer period than some previous studies, since Pca often takes 7-10 years to shift from inert to aggressive (75), more long-term follow-up data are still needed to better evaluate the impact of TT on Pca risk.
5 Conclusions
TT is an effective treatment for patients with hypogonadism. As for the fear that testosterone treatment may result in abnormal prostate growth and Pca risk, however, our research shows that TT does not cause abnormal prostate growth, and may even reduce the risk of Pca. In terms of the preparation type of TT, intramuscular seems to be more prominent in reducing the risk of Pca, and TT will not increase the risk of abnormal prostate growth regardless of the administration method. Further high-quality clinical research are required to confirm this observation.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
LY and QW were responsible for the idea and design of the study. BZ, XX and SQ did the analysis and interpreted the results of the study. XS and ZZ searched the databases and extracted clinical data from the included studies. BZ and XX wrote the first draft of this manuscript. SQ, LY and QW revised the manuscript and wrote the final version. All authors contributed to the article and approved the submitted version.
Funding
The work was supported by National Natural Science Foundation of China (Grant No.81902578, 81974098), China Post-docotral Science Foundation (2017M612971), Post-doctoral Science Research Foundation of Sichuan University (2020SCU12041), Post-Doctor Research Project, West China Hospital, Sichuan University (2018HYBH085), National Clinical Research Center for Geriatrics, West China Hospital, Sichuan University (Z2018C01).
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/fendo.2022.1009900/full#supplementary-material
References
1. Nackeeran S, Kohn T, Gonzalez D, White J, Ory J, Ramasamy R. The effect of route of testosterone on changes in hematocrit: A systematic review and Bayesian network meta-analysis of randomized trials. J Urol (2022) 207(1):44–51. doi: 10.1097/JU.0000000000002188
2. Wu FC, Tajar A, Pye SR, Silman AJ, Finn JD, O’Neill TW, et al. Hypothalamic-Pituitary-Testicular axis disruptions in older men are differentially linked to age and modifiable risk factors: The European Male aging study. J Clin Endocrinol Metab (2008) 93(7):2737–45. doi: 10.1210/jc.2007-1972
3. Kaufman JM, Vermeulen A. The decline of androgen levels in elderly men and its clinical and therapeutic implications. Endocr Rev (2005) 26(6):833–76. doi: 10.1210/er.2004-0013
4. Eau Guidelines on Sexual and Reproductive Health. Presented at the Eau Annual Congress Amsterdam 2022. Isbn 978-94-92671-16-5. Available at: Https://D56bochluxqnz.Cloudfront.Net/Documents/Full-Guideline/Eau-Guidelines-on-Sexual-and-Reproductive-Health-2022_2022-03-29-084141_Megw.Pdf.
5. Hall SA, Esche GR, Araujo AB, Travison TG, Clark RV, Williams RE, et al. Correlates of low testosterone and symptomatic androgen deficiency in a population-based sample. J Clin Endocrinol Metab (2008) 93(10):3870–7. doi: 10.1210/jc.2008-0021
6. Basaria S. Male Hypogonadism. Lancet (2014) 383(9924):1250–63. doi: 10.1016/s0140-6736(13)61126-5
7. Kelly DM, Jones TH. Testosterone and obesity. Obes Rev (2015) 16(7):581–606. doi: 10.1111/obr.12282
8. Muller M, Grobbee DE, den Tonkelaar I, Lamberts SW, van der Schouw YT. Endogenous sex hormones and metabolic syndrome in aging men. J Clin Endocrinol Metab (2005) 90(5):2618–23. doi: 10.1210/jc.2004-1158
9. Hudson J, Cruickshank M, Quinton R, Aucott L, Aceves-Martins M, Gillies K, et al. Adverse cardiovascular events and mortality in men during testosterone treatment: An individual patient and aggregate data meta-analysis. Lancet Healthy Longevity (2022) 3(6):e381–e93. doi: 10.1016/s2666-7568(22)00096-4
10. Barone B, Napolitano L, Abate M, Cirillo L, Reccia P, Passaro F, et al. The role of testosterone in the elderly: What do we know? Int J Mol Sci (2022) 23(7):3535. doi: 10.3390/ijms23073535
11. Thirumalai A, Cooper LA, Rubinow KB, Amory JK, Lin DW, Wright JL, et al. Stable intraprostatic dihydrotestosterone in healthy medically castrate men treated with exogenous testosterone. J Clin Endocrinol Metab (2016) 101(7):2937–44. doi: 10.1210/jc.2016-1483
12. Marks LS, Mazer NA, Mostaghel E, Hess DL, Dorey FJ, Epstein JI, et al. Effect of testosterone replacement therapy on prostate tissue in men with late-onset hypogonadism: A randomized controlled trial. JAMA (2006) 296(19):2351–61. doi: 10.1001/jama.296.19.2351
13. Baillargeon J, Kuo YF, Fang X, Shahinian VB. Long-term exposure to testosterone therapy and the risk of high grade prostate cancer. J Urol (2015) 194(6):1612–6. doi: 10.1016/j.juro.2015.05.099
14. Pastuszak AW, Gomez LP, Scovell JM, Khera M, Lamb DJ, Lipshultz LI. Comparison of the effects of testosterone gels, injections, and pellets on serum hormones, erythrocytosis, lipids, and prostate-specific antigen. Sex Med (2015) 3(3):165–73. doi: 10.1002/sm2.76
15. Page MJ, McKenzie JE, Bossuyt PM, Boutron I, Hoffmann TC, Mulrow CD, et al. The prisma 2020 statement: An updated guideline for reporting systematic reviews. BMJ (2021) 372:n71. doi: 10.1136/bmj.n71
16. Hutton B, Salanti G, Caldwell DM, Chaimani A, Schmid CH, Cameron C, et al. The prisma extension statement for reporting of systematic reviews incorporating network meta-analyses of health care interventions: Checklist and explanations. Ann Internal Med (2015) 162(11):777–84. doi: 10.7326/M14-2385
17. Cui Y, Zong H, Yan H, Zhang Y. The effect of testosterone replacement therapy on prostate cancer: A systematic review and meta-analysis. Prostate Cancer Prostatic Dis (2014) 17(2):132–43. doi: 10.1038/pcan.2013.60
18. Sterne JAC, Savovic J, Page MJ, Elbers RG, Blencowe NS, Boutron I, et al. Rob 2: A revised tool for assessing risk of bias in randomised trials. BMJ (2019) 366:l4898. doi: 10.1136/bmj.l4898
19. Stang A. Critical evaluation of the Newcastle-Ottawa scale for the assessment of the quality of nonrandomized studies in meta-analyses. Eur J Epidemiol (2010) 25(9):603–5. doi: 10.1007/s10654-010-9491-z
20. Brozek JL, Akl EA, Alonso-Coello P, Lang D, Jaeschke R, Williams JW, et al. Grading quality of evidence and strength of recommendations in clinical practice guidelines. part 1 of 3. an overview of the grade approach and grading quality of evidence about interventions. Allergy (2009) 64(5):669–77. doi: 10.1111/j.1398-9995.2009.01973.x
21. Dias S, Welton NJ, Caldwell DM, Ades AE. Checking consistency in mixed treatment comparison meta-analysis. Stat Med (2010) 29(7-8):932–44. doi: 10.1002/sim.3767
22. Hozo SP, Djulbegovic B, Hozo I. Estimating the mean and variance from the median, range, and the size of a sample. BMC Med Res Method (2005) 5:13. doi: 10.1186/1471-2288-5-13
23. Higgins JP, Thompson SG. Quantifying heterogeneity in a meta-analysis. Stat Med (2002) 21(11):1539–58. doi: 10.1002/sim.1186
24. Eng C, Kramer CK, Zinman B, Retnakaran R. Glucagon-like peptide-1 receptor agonist and basal insulin combination treatment for the management of type 2 diabetes: A systematic review and meta-analysis. Lancet (2014) 384(9961):2228–34. doi: 10.1016/s0140-6736(14)61335-0
25. Greco T, Landoni G, Biondi-Zoccai G, D’Ascenzo F, Zangrillo A. A Bayesian network meta-analysis for binary outcome: How to do it. Stat Methods Med Res (2016) 25(5):1757–73. doi: 10.1177/0962280213500185
26. Gelman A, Rubin DB. Inference from iterative simulation using multiple sequences. Stat Sci (1992) 7(4):457–72. doi: 10.1214/ss/1177011136
27. Brooks SP, Gelman A. General methods for monitoring convergence of iterative simulations. J Comput Graphical Stat (1998) 7(4):434–55. doi: 10.1080/10618600.1998.10474787
28. van Valkenhoef G, Lu G, de Brock B, Hillege H, Ades AE, Welton NJ. Automating network meta-analysis. Res Synth Methods (2012) 3(4):285–99. doi: 10.1002/jrsm.1054
29. Rhee H, Navaratnam A, Oleinikova I, Gilroy D, Scuderi Y, Heathcote P, et al. A novel liver-targeted testosterone therapy for sarcopenia in androgen deprived men with prostate cancer. J Endocr Soc (2021) 5(9):bvab116. doi: 10.1210/jendso/bvab116
30. Cunningham GR, Ellenberg SS, Bhasin S, Matsumoto AM, Parsons JK, Preston P, et al. Prostate-specific antigen levels during testosterone treatment of hypogonadal older men: Data from a controlled trial. J Clin Endocrinol Metab (2019) 104(12):6238–46. doi: 10.1210/jc.2019-00806
31. Brock G, Heiselman D, Maggi M, Kim SW, Rodriguez Vallejo JM, Behre HM, et al. Effect of testosterone solution 2% on testosterone concentration, sex drive and energy in hypogonadal men: Results of a placebo controlled study. J Urol (2016) 195(3):699–705. doi: 10.1016/j.juro.2015.10.083
32. Snyder PJ, Bhasin S, Cunningham GR, Matsumoto AM, Stephens-Shields AJ, Cauley JA, et al. Effects of testosterone treatment in older men. N Engl J Med (2016) 374(7):611–24. doi: 10.1056/NEJMoa1506119
33. Hackett G, Cole N, Bhartia M, Kennedy D, Raju J, Wilkinson P. Testosterone replacement therapy improves metabolic parameters in hypogonadal men with type 2 diabetes but not in men with coexisting depression: The blast study. J Sex Med (2014) 11(3):840–56. doi: 10.1111/jsm.12404
34. Borst SE, Yarrow JF, Conover CF, Nseyo U, Meuleman JR, Lipinska JA, et al. Musculoskeletal and prostate effects of combined testosterone and finasteride administration in older hypogonadal men: A randomized, controlled trial. Am J Physiol Endocrinol Metab (2014) 306(4):E433–42. doi: 10.1152/ajpendo.00592.2013
35. Del Fabbro E, Garcia JM, Dev R, Hui D, Williams J, Engineer D, et al. Testosterone replacement for fatigue in hypogonadal ambulatory males with advanced cancer: A preliminary double-blind placebo-controlled trial. Supportive Care Cancer (2013) 21(9):2599–607. doi: 10.1007/s00520-013-1832-5
36. Hildreth KL, Barry DW, Moreau KL, Vande Griend J, Meacham RB, Nakamura T, et al. Effects of testosterone and progressive resistance exercise in healthy, highly functioning older men with low-normal testosterone levels. J Clin Endocrinol Metab (2013) 98(5):1891–900. doi: 10.1210/jc.2013-2227. 10.1210/jc.2012-3695.
37. Behre HM, Tammela TL, Arver S, Tolra JR, Bonifacio V, Lamche M, et al. A randomized, double-blind, placebo-controlled trial of testosterone gel on body composition and health-related quality-of-Life in men with hypogonadal to low-normal levels of serum testosterone and symptoms of androgen deficiency over 6 months with 12 months open-label follow-up. Aging Male (2012) 15(4):198–207. doi: 10.3109/13685538.2012.699562
38. Jones TH, Arver S, Behre HM, Buvat J, Meuleman E, Moncada I, et al. Testosterone replacement in hypogonadal men with type 2 diabetes and/or metabolic syndrome (the Times2 study). Diabetes Care (2011) 34(4):828–37. doi: 10.2337/dc10-1233
39. Basaria S, Coviello AD, Travison TG, Storer TW, Farwell WR, Jette AM, et al. Adverse events associated with testosterone administration. N Engl J Med (2010) 363(2):109–22. doi: 10.1056/NEJMoa1000485
40. Kaufman JM, Miller MG, Garwin JL, Fitzpatrick S, McWhirter C, Brennan JJ. Efficacy and safety study of 1.62% testosterone gel for the treatment of hypogonadal men. J Sex Med (2011) 8(7):2079–89. doi: 10.1111/j.1743-6109.2011.02265.x
41. Kenny AM, Prestwood KM, Gruman CA, Marcello KM, Raisz LG. Effects of transdermal testosterone on bone and muscle in older men with low bioavailable testosterone levels. J Gerontol A Biol Sci Med Sci (2001) 56(5):M266–72. doi: 10.1093/gerona/56.5.m266
42. Kenny AM, Kleppinger A, Annis K, Rathier M, Browner B, Judge JO, et al. Effects of transdermal testosterone on bone and muscle in older men with low bioavailable testosterone levels, low bone mass, and physical frailty. J Am Geriatr Soc (2010) 58(6):1134–43. doi: 10.1111/j.1532-5415.2010.02865.x
43. Nair KS, Rizza RA, O’Brien P, Dhatariya K, Short KR, Nehra A, et al. Dhea in elderly women and dhea or testosterone in elderly men. N Engl J Med (2006) 355(16):1647–59. doi: 10.1056/NEJMoa054629
44. Simon D, Charles MA, Lahlou N, Nahoul K, Oppert JM, Gouault-Heilmann M, et al. Androgen therapy improves insulin sensitivity and decreases leptin level in healthy adult men with low plasma total testosterone: A 3-month randomized placebo-controlled trial. Diabetes Care (2001) 24(12):2149–51. doi: 10.2337/diacare.24.12.2149
45. Srinivas-Shankar U, Roberts SA, Connolly MJ, O’Connell MD, Adams JE, Oldham JA, et al. Effects of testosterone on muscle strength, physical function, body composition, and quality of life in intermediate-frail and frail elderly men: A randomized, double-blind, placebo-controlled study. J Clin Endocrinol Metab (2010) 95(2):639–50. doi: 10.1210/jc.2009-1251
46. Steidle C, Schwartz S, Jacoby K, Sebree T, Smith T, Bachand R. Aa2500 testosterone gel normalizes androgen levels in aging males with improvements in body composition and sexual function. J Clin Endocrinol Metab (2003) 88(6):2673–81. doi: 10.1210/jc.2002-021058
47. Amory JK, Watts NB, Easley KA, Sutton PR, Anawalt BD, Matsumoto AM, et al. Exogenous testosterone or testosterone with finasteride increases bone mineral density in older men with low serum testosterone. J Clin Endocrinol Metab (2004) 89(2):503–10. doi: 10.1210/jc.2003-031110
48. Aversa A, Bruzziches R, Francomano D, Rosano G, Isidori AM, Lenzi A, et al. Effects of testosterone undecanoate on cardiovascular risk factors and atherosclerosis in middle-aged men with late-onset hypogonadism and metabolic syndrome: Results from a 24-month, randomized, double-blind, placebo-controlled study. J Sexual Med (2010) 7(10):3495–503. doi: 10.1111/j.1743-6109.2010.01931.x
49. Ferrando AA, Sheffield-Moore M, Yeckel CW, Gilkison C, Jiang J, Achacosa A, et al. Testosterone administration to older men improves muscle function: Molecular and physiological mechanisms. Am J Physiol Endocrinol Metab (2002) 282(3):E601–7. doi: 10.1152/ajpendo.00362.2001
50. Kalinchenko SY, Tishova YA, Mskhalaya GJ, Gooren LJ, Giltay EJ, Saad F. Effects of testosterone supplementation on markers of the metabolic syndrome and inflammation in hypogonadal men with the metabolic syndrome: The double-blinded placebo-controlled Moscow study. Clin Endocrinol (Oxf) (2010) 73(5):602–12. doi: 10.1111/j.1365-2265.2010.03845.x
51. Kenny AM, Fabregas G, Song C, Biskup B, Bellantonio S. Effects of testosterone on behavior, depression, and cognitive function in older men with mild cognitive loss. J Gerontol A Biol Sci Med Sci (2004) 59(1):75–8. doi: 10.1093/gerona/59.1.M75
52. Shigehara K, Sugimoto K, Konaka H, Iijima M, Fukushima M, Maeda Y, et al. Androgen replacement therapy contributes to improving lower urinary tract symptoms in patients with hypogonadism and benign prostate hypertrophy: A randomised controlled study. Aging Male (2011) 14(1):53–8. doi: 10.3109/13685538.2010.518178
53. Emmelot-Vonk MH, Verhaar HJ, Nakhai PHR, Aleman A, Lock TM, Bosch JL, et al. Effect of testosterone supplementation on functional mobility, cognition, and other parameters in older men: A randomized controlled trial. JAMA (2008) 299(1):39–52. doi: 10.1001/jama.2007.51
54. Wittert GA, Chapman IM, Haren MT, Mackintosh S, Coates P, Morley JE. Oral testosterone supplementation increases muscle and decreases fat mass in healthy elderly males with low-normal gonadal status. J Gerontol A Biol Sci Med Sci (2003) 58(7):618–25. doi: 10.1093/gerona/58.7.M618
55. Wittert G, Bracken K, Robledo KP, Grossmann M, Yeap BB, Handelsman DJ, et al. Testosterone treatment to prevent or revert type 2 diabetes in men enrolled in a lifestyle programme (T4dm): A randomised, double-blind, placebo-controlled, 2-year, phase 3b trial. Lancet Diabetes Endocrinol (2021) 9(1):32–45. doi: 10.1016/s2213-8587(20)30367-3
56. Groti Antonič K, Antonič B, Pfeifer M. Effects of testosterone therapy on erythrocytosis and prostate adverse events in obese males with functional hypogonadism and type 2 diabetes in a 2-year clinical trial. Androgens: Clin Res Ther (2020) 1(1):85–93. doi: 10.1089/andro.2020.0008
57. Saad F, Caliber M, Doros G, Haider KS, Haider A. Long-term treatment with testosterone undecanoate injections in men with hypogonadism alleviates erectile dysfunction and reduces risk of major adverse cardiovascular events, prostate cancer, and mortality. Aging Male (2020) 23(1):81–92. doi: 10.1080/13685538.2019.1575354
58. Pastuszak AW, Pearlman AM, Lai WS, Godoy G, Sathyamoorthy K, Liu JS, et al. Testosterone replacement therapy in patients with prostate cancer after radical prostatectomy. J Urol (2013) 190(2):639–44. doi: 10.1016/j.juro.2013.02.002
59. Rhoden EL, Morgentaler A. Influence of demographic factors and biochemical characteristics on the prostate-specific antigen (Psa) response to testosterone replacement therapy. Int J Impot Res (2006) 18(2):201–5. doi: 10.1038/sj.ijir.3901394
60. Guay AT, Perez JB, Fitaihi WA, Vereb M. Testosterone treatment in hypogonadal men: Prostate-specific antigen level and risk of prostate cancer. Endocr Pract (2000) 6(2):132–8. doi: 10.4158/EP.6.2.132
61. Zitzmann M, Cremers JF, Krallmann C, Kliesch S. The heat-registry (Hematopoietic affection by testosterone): Comparison of a transdermal gel vs long-acting intramuscular testosterone undecanoate in hypogonadal men. Aging Male (2022) 25(1):134–44. doi: 10.1080/13685538.2022.2063830
62. Allott EH, Masko Em, Freedland SJ. Obesity and Prostate Cancer: Weighing the Evidence. Eur Urol (2013) 63(5):800–9. doi: 10.1016/j.eururo.2012.11.013
63. Skinner JW, Otzel DM, Bowser A, Nargi D, Agarwal S, Peterson MD, et al. Muscular responses to testosterone replacement vary by administration route: A systematic review and meta-analysis. J Cachexia Sarcopenia Muscle (2018) 9(3):465–81. doi: 10.1002/jcsm.12291
64. Rastrelli G, Vignozzi L, Corona G, Maggi M. Testosterone and benign prostatic hyperplasia. Sexual Med Rev (2019) 7(2):259–71. doi: 10.1016/j.sxmr.2018.10.006
65. Guo C, Gu W, Liu M, Peng BO, Yao X, Yang B, et al. Efficacy and safety of testosterone replacement therapy in men with hypogonadism: A meta-analysis study of placebo-controlled trials. Exp Ther Med (2016) 11(3):853–63. doi: 10.3892/etm.2015.2957
66. Fernández-Balsells MM, Murad MH, Lane M, Lampropulos JF, Albuquerque F, Mullan RJ, et al. Clinical review 1: Adverse effects of testosterone therapy in adult men: A systematic review and meta-analysis. J Clin Endocrinol Metab (2010) 95(6):2560–75. doi: 10.1210/jc.2009-2575
67. Morgentaler A, Traish AM. Shifting the paradigm of testosterone and prostate cancer: The saturation model and the limits of androgen-dependent growth. Eur Urol (2009) 55(2):310–20. doi: 10.1016/j.eururo.2008.09.024
68. Debruyne FM, Behre HM, Roehrborn CG, Maggi M, Wu FC, Schröder FH, et al. Testosterone treatment is not associated with increased risk of prostate cancer or worsening of lower urinary tract symptoms: Prostate health outcomes in the registry of hypogonadism in men. BJU Int (2017) 119(2):216–24. doi: 10.1111/bju.13578
69. Rastrelli G, Giovannini L, Calogero AE, Gianfrilli D, Serra E, Pizzocaro A, et al. Predictors and clinical consequences of starting androgen therapy in men with low testosterone: Results from the siamo-noi registry. J Endocrinological Invest (2016) 39(6):695–708. doi: 10.1007/s40618-016-0461-5
70. Kang DY, Li HJ. The effect of testosterone replacement therapy on prostate-specific antigen (Psa) levels in men being treated for hypogonadism: A systematic review and meta-analysis. Medicine (2015) 94(3):e410. doi: 10.1097/md.0000000000000410
71. Haider A, Zitzmann M, Doros G, Isbarn H, Hammerer P, Yassin A. Incidence of prostate cancer in hypogonadal men receiving testosterone therapy: Observations from 5-year median followup of 3 registries. J Urol (2015) 193(1):80–6. doi: 10.1016/j.juro.2014.06.071
72. Wallis CJ, Lo K, Lee Y, Krakowsky Y, Garbens A, Satkunasivam R, et al. Survival and cardiovascular events in men treated with testosterone replacement therapy: An intention-to-Treat observational cohort study. Lancet Diabetes Endocrinol (2016) 4(6):498–506. doi: 10.1016/s2213-8587(16)00112-1
73. Xu X, Chen X, Hu H, Dailey AB, Taylor BD. Current opinion on the role of testosterone in the development of prostate cancer: A dynamic model. BMC Cancer (2015) 15:806. doi: 10.1186/s12885-015-1833-5
74. Cui Y, Zhang Y. The effect of androgen-replacement therapy on prostate growth: A systematic review and meta-analysis. Eur Urol (2013) 64(5):811–22. doi: 10.1016/j.eururo.2013.03.042
Keywords: hypogonadism, administrations, prostate safety, meta-analysis, systematic review, testosterone therapy
Citation: Zeng B, Qiu S, Xiong X, Su X, Zhang Z, Wei Q and Yang L (2022) The effect of different administrations of testosterone therapy on adverse prostate events: A Bayesian network meta-analysis. Front. Endocrinol. 13:1009900. doi: 10.3389/fendo.2022.1009900
Received: 02 August 2022; Accepted: 17 October 2022;
Published: 07 November 2022.
Edited by:
Haolin Chen, Wenzhou Medical University, ChinaCopyright © 2022 Zeng, Qiu, Xiong, Su, Zhang, Wei and Yang. 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: Qiang Wei, weiqiang933@163.com; Lu Yang, wycleflue@scu.edu.cn
†These authors have contributed equally to this work