1
|
Jackson D, Sweeting M, van Aert RCM, Bujkiewicz S, Abrams KR, Viechtbauer W. A New Frequentist Implementation of the Daniels and Hughes Bivariate Meta-Analysis Model for Surrogate Endpoint Evaluation. Biom J 2025; 67:e70048. [PMID: 40105204 PMCID: PMC11921291 DOI: 10.1002/bimj.70048] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2024] [Revised: 01/23/2025] [Accepted: 02/13/2025] [Indexed: 03/20/2025]
Abstract
Surrogate endpoints are used when the primary outcome is difficult to measure accurately. Determining if a measure is suitable to use as a surrogate endpoint is a challenging task and a variety of meta-analysis models have been proposed for this purpose. The Daniels and Hughes bivariate model for trial-level surrogate endpoint evaluation is gaining traction but presents difficulties for frequentist estimation and hitherto only Bayesian solutions have been available. This is because the marginal model is not a conventional linear model and the number of unknown parameters increases at the same rate as the number of studies. This second property raises immediate concerns that the maximum likelihood estimator of the model's unknown variance component may be downwardly biased. We derive maximum likelihood estimating equations to motivate a bias adjusted estimator of this parameter. The bias correction terms in our proposed estimating equation are easily computed and have an intuitively appealing algebraic form. A simulation study is performed to illustrate how this estimator overcomes the difficulties associated with maximum likelihood estimation. We illustrate our methods using two contrasting examples from oncology.
Collapse
Affiliation(s)
- Dan Jackson
- Statistical Innovation Group, AstraZeneca, Cambridge, UK
| | | | | | - Sylwia Bujkiewicz
- Biostatistics Research Group, Department of Population Health Sciences, University of Leicester, Leicester, UK
| | - Keith R Abrams
- Department of Statistics and Warwick Medical School, University of Warwick, Coventry, UK
| | | |
Collapse
|
2
|
Yang L, Raveendran G, Meng X, Lin J, Meng Z. Predict progression free survival and overall survival using objective response rate for anti-PD1/PDL1 therapy development. BMC Cancer 2024; 24:912. [PMID: 39075397 PMCID: PMC11287896 DOI: 10.1186/s12885-024-12664-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2024] [Accepted: 07/19/2024] [Indexed: 07/31/2024] Open
Abstract
In oncology anti-PD1 / PDL1 therapy development for solid tumors, objective response rate (ORR) is commonly used clinical endpoint for early phase study decision making, while progression free survival (PFS) and overall survival (OS) are widely used for late phase study decision making. Developing predictive models to late phase outcomes such as median PFS (mPFS) and median OS (mOS) based on early phase clinical outcome ORR could inform late phase study design optimization and probability of success (POS) evaluation. In existing literature, there are ORR / mPFS / mOS association and surrogacy investigations with limited number of included clinical trials. In this paper, without establishing surrogacy, we attempt to predict late phase survival (mPFS and mOS) based on early efficacy ORR and optimize late phase trial design for anti-PD1 / PDL1 therapy development. In order to include adequate number of eligible clinical trials, we built a comprehensive quantitative clinical trial landscape database (QLD) by combining information from different sources such as clinicaltrial.gov, publications, company press releases for relevant indications and therapies. We developed a generalizable algorithm to systematically extract structured data for scientific accuracy and completeness. Finally, more than 150 late phase clinical trials were identified for ORR / mPFS (ORR / mOS) predictive model development while existing literature included at most 50 trials. A tree-based machine learning regression model has been derived to account for ORR / mPFS (ORR / mOS) relationship heterogeneity across tumor type, stage, line of therapy, treatment class and borrow strength simultaneously when homogeneity persists. The proposed method ensures that the predictive model is robust and have explicit structure for clinical interpretation. Through cross validation, the average predictive mean square error of the proposed model is competitive to random forest and extreme gradient boosting methods and outperforms commonly used additive or interaction linear regression models. An example application of the proposed ORR / mPFS (ORR / mOS) predictive model on late phase trial POS evaluation for anti-PD1 / PDL1 combination therapy was illustrated.
Collapse
Affiliation(s)
- Lei Yang
- Sanofi Bridgewater, New Jersey, US
| | | | | | - Ji Lin
- Sanofi, Cambridge, MA, US
| | | |
Collapse
|
3
|
Waser NA, Quintana M, Schweikert B, Chaft JE, Berry L, Adam A, Vo L, Penrod JR, Fiore J, Berry DA, Goring S. Pathological response in resectable non-small cell lung cancer: a systematic literature review and meta-analysis. JNCI Cancer Spectr 2024; 8:pkae021. [PMID: 38521542 PMCID: PMC11101053 DOI: 10.1093/jncics/pkae021] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2023] [Revised: 12/06/2023] [Accepted: 03/15/2024] [Indexed: 03/25/2024] Open
Abstract
BACKGROUND Surrogate endpoints for overall survival in patients with resectable non-small cell lung cancer receiving neoadjuvant therapy are needed to provide earlier treatment outcome indicators and accelerate drug approval. This study's main objectives were to investigate the association among pathological complete response, major pathological response, event-free survival and overall survival and to determine whether treatment effects on pathological complete response and event-free survival correlate with treatment effects on overall survival. METHODS A comprehensive systematic literature review was conducted to identify neoadjuvant studies in resectable non-small cell lung cancer. Analysis at the patient level using frequentist and Bayesian random effects (hazard ratio [HR] for overall survival or event-free survival by pathological complete response or major pathological response status, yes vs no) and at the trial level using weighted least squares regressions (hazard ratio for overall survival or event-free survival vs pathological complete response, by treatment arm) were performed. RESULTS In both meta-analyses, pathological complete response yielded favorable overall survival compared with no pathological complete response (frequentist, 20 studies and 6530 patients: HR = 0.49, 95% confidence interval = 0.42 to 0.57; Bayesian, 19 studies and 5988 patients: HR = 0.48, 95% probability interval = 0.43 to 0.55) and similarly for major pathological response (frequentist, 12 studies and 1193 patients: HR = 0.36, 95% confidence interval = 0.29 to 0.44; Bayesian, 11 studies and 1018 patients: HR = 0.33, 95% probability interval = 0.26 to 0.42). Across subgroups, estimates consistently showed better overall survival or event-free survival in pathological complete response or major pathological response compared with no pathological complete response or no major pathological response. Trial-level analyses showed a moderate to strong correlation between event-free survival and overall survival hazard ratios (R2 = 0.7159) but did not show a correlation between treatment effects on pathological complete response and overall survival or event-free survival. CONCLUSION There was a strong and consistent association between pathological response and survival and a moderate to strong correlation between event-free survival and overall survival following neoadjuvant therapy for patients with resectable non-small cell lung cancer.
Collapse
Affiliation(s)
| | | | | | - Jamie E Chaft
- Department of Medicine, Memorial Sloan Kettering Cancer Center, New York, NY, USA
| | | | - Ahmed Adam
- Insights, Evidence and Value, ICON plc, Burlington, ON, Canada
| | - Lien Vo
- Health Economics and Outcomes Research, Bristol Myers Squibb, Lawrenceville, NJ, USA
| | - John R Penrod
- Health Economics and Outcomes Research, Bristol Myers Squibb, Lawrenceville, NJ, USA
| | - Joseph Fiore
- Health Economics and Outcomes Research, Bristol Myers Squibb, Lawrenceville, NJ, USA
| | | | - Sarah Goring
- Insights, Evidence and Value, ICON plc, Burlington, ON, Canada
| |
Collapse
|
4
|
Liu H, Milenković‐Grišić A, Krishnan SM, Jönsson S, Friberg LE, Girard P, Venkatakrishnan K, Vugmeyster Y, Khandelwal A, Karlsson MO. A multistate modeling and simulation framework to learn dose-response of oncology drugs: Application to bintrafusp alfa in non-small cell lung cancer. CPT Pharmacometrics Syst Pharmacol 2023; 12:1738-1750. [PMID: 37165943 PMCID: PMC10681430 DOI: 10.1002/psp4.12976] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2023] [Revised: 04/06/2023] [Accepted: 04/14/2023] [Indexed: 05/12/2023] Open
Abstract
The dose/exposure-efficacy analyses are often conducted separately for oncology end points like best overall response, progression-free survival (PFS) and overall survival (OS). Multistate models offer to bridge these dose-end point relationships by describing transitions and transition times from enrollment to response, progression, and death, and evaluating transition-specific dose effects. This study aims to apply the multistate pharmacometric modeling and simulation framework in a dose optimization setting of bintrafusp alfa, a fusion protein targeting TGF-β and PD-L1. A multistate model with six states (stable disease [SD], response, progression, unknown, dropout, and death) was developed to describe the totality of endpoints data (time to response, PFS, and OS) of 80 patients with non-small cell lung cancer receiving 500 or 1200 mg of bintrafusp alfa. Besides dose, evaluated predictor of transitions include time, demographics, premedication, disease factors, individual clearance derived from a pharmacokinetic model, and tumor dynamic metrics observed or derived from tumor size model. We found that probabilities of progression and death upon progression decreased over time since enrollment. Patients with metastasis at baseline had a higher probability to progress than patients without metastasis had. Despite dose failed to be statistically significant for any individual transition, the combined effect quantified through a model with dose-specific transition estimates was still informative. Simulations predicted a 69.2% probability of at least 1 month longer, and, 55.6% probability of at least 2-months longer median OS from the 1200 mg compared to the 500 mg dose, supporting the selection of 1200 mg for future studies.
Collapse
Affiliation(s)
- Han Liu
- Department of PharmacyUppsala UniversityUppsalaSweden
| | | | | | - Siv Jönsson
- Department of PharmacyUppsala UniversityUppsalaSweden
| | | | - Pascal Girard
- Merck Institute of Pharmacometrics, an affiliate of Merck KGaALausanneSwitzerland
| | - Karthik Venkatakrishnan
- EMD Serono Research & Development Institute, Inc., an affiliate of Merck KGaABillericaMassachusettsUSA
| | - Yulia Vugmeyster
- EMD Serono Research & Development Institute, Inc., an affiliate of Merck KGaABillericaMassachusettsUSA
| | | | | |
Collapse
|
5
|
Turner DC, Wada R, Zhou H, Wang X, de Greef R, Valiathan C, Zhang L, Zhang N, Kuchimanchi M, Chen T, Ballas M, Visser SAG. Model-based meta-analysis of non-small cell lung cancer with standard of care PD-1 inhibitors and chemotherapy for early development decision making. CPT Pharmacometrics Syst Pharmacol 2023; 12:1751-1763. [PMID: 36642813 PMCID: PMC10681483 DOI: 10.1002/psp4.12917] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2022] [Revised: 12/22/2022] [Accepted: 12/27/2022] [Indexed: 01/17/2023] Open
Abstract
Single-arm cohorts/trials are often used in early phase oncology programs to support preliminary clinical activity assessments for investigational products, administered alone or in combination with standard of care (SOC) agents. Benchmarking clinical activity of those combinations against other treatments, including SOC, requires indirect comparisons against external trials, which presents challenges including cross-study differences in trial populations/other factors. To facilitate such nonrandomized comparisons, we developed a comprehensive model-based meta-analysis (MBMA) framework to quantitatively adjust for factors related to efficacy in metastatic non-small cell lung cancer (mNSCLC). Data were derived from 15 published studies assessing key programmed cell death protein-1 (PD-1) inhibitors pembrolizumab (n = 8) and nivolumab (n = 7), representing current SOC in mNSCLC. In the first stage, a mixed-effects logistic regression model for overall response rate (ORR) was developed accounting for effects of various population covariates on ORR. The ORR model results indicated an odds ratio (OR) of 1.02 for squamous versus non-squamous histology and OR of 1.20 for PD-ligand 1 tumor proportion score (TPS) per every 10% increase of TPS level. Next, a nonparametric mixed-effects model for overall survival (OS) was developed with ORR/other clinical covariates as input. Subsequently, MBMA simulations of relevant hypothetical scenarios involving single-arm trial design predicted OS hazard ratios as a function of ORR with matched patient characteristics. Findings from this MBMA and derived parameter estimates can be generally applied by the reader as a framework for interpreting efficacy data from early phase trials to support ORR-based go/no-go decisions and futility rules, illustrated through examples in this report.
Collapse
Affiliation(s)
- David C. Turner
- GSKCollegevillePennsylvaniaUSA
- Present address:
GenentechSan FranciscoCaliforniaUSA
| | - Russ Wada
- CertaraMenlo ParkCaliforniaUSA
- Present address:
QuanTx ConsultingMountain ViewCaliforniaUSA
| | | | - Xiaowei Wang
- GSKCollegevillePennsylvaniaUSA
- Present address:
ModernaCambridgeMassachusettsUSA
| | | | - Chandni Valiathan
- GSKCollegevillePennsylvaniaUSA
- Present address:
J&JNew BrunswickNew JerseyUSA
| | | | | | | | | | - Marc Ballas
- GSKCollegevillePennsylvaniaUSA
- Present address:
NovocurePotomacMarylandUSA
| | | |
Collapse
|
6
|
Saad ED, Coart E, Deltuvaite-Thomas V, Garcia-Barrado L, Burzykowski T, Buyse M. Trial Design for Cancer Immunotherapy: A Methodological Toolkit. Cancers (Basel) 2023; 15:4669. [PMID: 37760636 PMCID: PMC10527464 DOI: 10.3390/cancers15184669] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2023] [Revised: 08/12/2023] [Accepted: 09/05/2023] [Indexed: 09/29/2023] Open
Abstract
Immunotherapy with checkpoint inhibitors (CPIs) and cell-based products has revolutionized the treatment of various solid tumors and hematologic malignancies. These agents have shown unprecedented response rates and long-term benefits in various settings. These clinical advances have also pointed to the need for new or adapted approaches to trial design and assessment of efficacy and safety, both in the early and late phases of drug development. Some of the conventional statistical methods and endpoints used in other areas of oncology appear to be less appropriate in immuno-oncology. Conversely, other methods and endpoints have emerged as alternatives. In this article, we discuss issues related to trial design in the early and late phases of drug development in immuno-oncology, with a focus on CPIs. For early trials, we review the most salient issues related to dose escalation, use and limitations of tumor response and progression criteria for immunotherapy, the role of duration of response as an endpoint in and of itself, and the need to conduct randomized trials as early as possible in the development of new therapies. For late phases, we discuss the choice of primary endpoints for randomized trials, review the current status of surrogate endpoints, and discuss specific statistical issues related to immunotherapy, including non-proportional hazards in the assessment of time-to-event endpoints, alternatives to the Cox model in these settings, and the method of generalized pairwise comparisons, which can provide a patient-centric assessment of clinical benefit and be used to design randomized trials.
Collapse
Affiliation(s)
- Everardo D. Saad
- International Drug Development Institute, Louvain-la-Neuve (IDDI), 1340 Ottignies-Louvain-la-Neuve, Belgium; (E.C.); (V.D.-T.); (L.G.-B.); (T.B.); (M.B.)
| | - Elisabeth Coart
- International Drug Development Institute, Louvain-la-Neuve (IDDI), 1340 Ottignies-Louvain-la-Neuve, Belgium; (E.C.); (V.D.-T.); (L.G.-B.); (T.B.); (M.B.)
| | - Vaiva Deltuvaite-Thomas
- International Drug Development Institute, Louvain-la-Neuve (IDDI), 1340 Ottignies-Louvain-la-Neuve, Belgium; (E.C.); (V.D.-T.); (L.G.-B.); (T.B.); (M.B.)
| | - Leandro Garcia-Barrado
- International Drug Development Institute, Louvain-la-Neuve (IDDI), 1340 Ottignies-Louvain-la-Neuve, Belgium; (E.C.); (V.D.-T.); (L.G.-B.); (T.B.); (M.B.)
| | - Tomasz Burzykowski
- International Drug Development Institute, Louvain-la-Neuve (IDDI), 1340 Ottignies-Louvain-la-Neuve, Belgium; (E.C.); (V.D.-T.); (L.G.-B.); (T.B.); (M.B.)
- Interuniversity Institute for Biostatistics and Statistical Bioinformatics (I-BioStat), Hasselt University, B-3500 Hasselt, Belgium
| | - Marc Buyse
- International Drug Development Institute, Louvain-la-Neuve (IDDI), 1340 Ottignies-Louvain-la-Neuve, Belgium; (E.C.); (V.D.-T.); (L.G.-B.); (T.B.); (M.B.)
- Interuniversity Institute for Biostatistics and Statistical Bioinformatics (I-BioStat), Hasselt University, B-3500 Hasselt, Belgium
| |
Collapse
|
7
|
Zhang Z, Pan Q, Lu M, Zhao B. Intermediate endpoints as surrogates for outcomes in cancer immunotherapy: a systematic review and meta-analysis of phase 3 trials. EClinicalMedicine 2023; 63:102156. [PMID: 37600482 PMCID: PMC10432823 DOI: 10.1016/j.eclinm.2023.102156] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/20/2023] [Revised: 07/24/2023] [Accepted: 07/26/2023] [Indexed: 08/22/2023] Open
Abstract
Background Cancer immunotherapy shows unique efficacy kinetics that differs from conventional treatment. These characteristics may lead to the prolongation of trial duration, hence reliable surrogate endpoints are urgently needed. We aimed to systematically evaluate the study-level performance of commonly reported intermediate clinical endpoints for surrogacy in cancer immunotherapy. Methods We searched the Embase, PubMed, and Cochrane databases, between database inception and October 18, 2022, for phase 3 randomised trials investigating the efficacy of immunotherapy in patients with advanced solid tumours. An updated search was done on July, 15, 2023. No language restrictions were used. Eligible trials had to set overall survival (OS) as the primary or co-primary endpoint and report at least one intermediate clinical endpoint including objective response rate (ORR), disease control rate (DCR), progression-free survival (PFS), and 1-year overall survival. Other key inclusion and exclusion criteria included: (1) adult patients (>18 years old) with advanced solid tumour; (2) no immunotherapy conducted in the control arms; (3) follow-up is long enough to achieve OS; (4) data should be public available. A two-stage meta-analytic approach was conducted to evaluate the magnitude of the association between these intermediate endpoints and OS. A surrogate was identified if the coefficient of determination (R2) was 0.7 or greater. Leave-one-out cross-validation and pre-defined subgroup analysis were conducted to examine the heterogeneity. Potential publication bias was evaluated using the Egger's and Begg's tests. This trial was registered with PROSPERO, number CRD42022381648. Findings 52,342 patients with 15 types of tumours from 77 phase 3 studies were included. ORR (R2 = 0.11; 95% CI, 0.00-0.24), DCR (R2 = 0.01; 95% CI, 0.00-0.01), and PFS (R2 = 0.40; 95% CI, 0.23-0.56) showed weak associations with OS. However, a strong correlation was observed between 1-year survival and clinical outcome (R2 = 0.74; 95% CI, 0.64-0.83). These associations remained relatively consistent across pre-defined subgroups stratified based on tumour types, masking methods, line of treatments, drug targets, treatment strategies, and follow-up durations. No significant heterogeneities or publication bias were identified. Interpretation 1-year milestone survival was the only identified surrogacy endpoint for outcomes in cancer immunotherapy. Ongoing investigations and development of new endpoints and incorporation of biomarkers are needed to identify potential surrogate markers that can be more robust than 1-year survival. This work may provide important references in assisting the design and interpretation of future clinical trials, and constitute complementary information in drafting clinical practice guidelines. Funding None.
Collapse
Affiliation(s)
- Zhishan Zhang
- Quanzhou First Hospital Affiliated to Fujian Medical University, Quanzhou, China
| | - Qunxiong Pan
- Quanzhou First Hospital Affiliated to Fujian Medical University, Quanzhou, China
| | - Mingdong Lu
- The Second Affiliated Hospital & Yuying Children's Hospital, Wenzhou Medical University, Wenzhou, China
| | - Bin Zhao
- Quanzhou First Hospital Affiliated to Fujian Medical University, Quanzhou, China
| |
Collapse
|
8
|
Hua T, Gao Y, Zhang R, Wei Y, Chen F. Validating ORR and PFS as surrogate endpoints in phase II and III clinical trials for NSCLC patients: difference exists in the strength of surrogacy in various trial settings. BMC Cancer 2022; 22:1022. [PMID: 36171546 PMCID: PMC9520950 DOI: 10.1186/s12885-022-10046-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2022] [Accepted: 08/31/2022] [Indexed: 11/10/2022] Open
Abstract
OBJECTIVE This study aims to systematically validate the performance of surrogate endpoints in phase II and III clinical trials for NSCLC patients under various trial settings. METHODS A literature search retrieved all registered phase II and III trials of NSCLC patients in which OS, with at least one of ORR and PFS, were reported. Associations between surrogate and true endpoints were assessed on two levels. On the arm level, three pairs of correlations, i.e., ORR vs. median OS, ORR vs. median PFS, and median PFS vs. median OS, were analysed using Spearman's rho. On the trial level, similarly, three pairs of correlations, i.e., ΔORR vs. HR of OS, ΔORR vs. HR of PFS, and HR of PFS vs. HR of OS, were analysed using Spearman's rho and weighted linear regression model respectively. Finally, sensitivity analyses were performed to explore surrogacy under various trial settings. RESULTS At arm level, three pairs of correlations are all high (Spearman's rho = 0.700, 0.831, 0.755, respectively). At trial level, there is a low correlation between ΔORR and HR of OS, a high correlation between ΔORR and HR of PFS and a moderate correlation between HR of PFS and HR of OS (Spearman's rho = 0.462, 0.764, 0.584, respectively). In the sensitivity analysis, we find correlations between surrogate and true endpoints vary by different trial settings. It is noteworthy that the strength of surrogacy of these intermediate endpoints in targeted therapy is greater than that in immunotherapy. CONCLUSION According to the arm-level and trial level-analysis, we suggest that in phase II and III trials of targeted therapy and immunotherapy for NSCLC patients: 1) ORR lacks validity for the surrogacy of OS, excluding in first-line therapy, and 2) ORR may be an appropriate surrogate endpoint for PFS, and 3) PFS may be considered a modest surrogacy for OS, with better performance in first-line therapy trials. Moreover, to provide more convincing evidence of surrogacy of the surrogate endpoints, patient-level analyses are in desperate need.
Collapse
Affiliation(s)
- Tiantian Hua
- Present Address: Department of Epidemiology and Health Statistics, School of Public Health, Southeast University, Nanjing, China.,Department of Biostatistics, School of Public Health, Nanjing Medical University, Nanjing, China
| | - Yuan Gao
- Department of Biostatistics, School of Public Health, Nanjing Medical University, Nanjing, China
| | - Ruyang Zhang
- Department of Biostatistics, School of Public Health, Nanjing Medical University, Nanjing, China
| | - Yongyue Wei
- Department of Biostatistics, School of Public Health, Nanjing Medical University, Nanjing, China
| | - Feng Chen
- Present Address: Department of Epidemiology and Health Statistics, School of Public Health, Southeast University, Nanjing, China. .,Department of Biostatistics, School of Public Health, Nanjing Medical University, Nanjing, China.
| |
Collapse
|