Introduction
The introduction of immune checkpoint inhibitors (ICIs) as cancer treatment has markedly transformed the prognosis of patients in different subtypes of metastatic cancer by enhancing cytotoxic T‐cells activity though the programmed cell death 1 (PD-1)–programmed cell death 1 ligand 1 (PD-L1) axis [1] , [2] , [3] , [4] . However, ICIs are correlated with a variety of toxicities as results of an excessively activated immune system and, therefore, considered as immune‐related adverse events (irAEs). Despite these toxicities may potentially affect all host tissues [5] , [6] , [7] , [8] , [9] , the effect on the skeleton is poorly investigated due to the few and small case series reported in literature [10] , [11] . These studies described skeletal adverse effects in patients treated with ICIs including serious vertebral fractures that, in some cases, affected multiple sites [10] , [11] .
Osteoimmunology is recently emerging as a new interdisciplinary field to investigate the bi-directional interplay between bone and immune systems. Indeed, systemic activation of T cells or alterations in pro-inflammatory cytokines production stimulates osteoclastogenesis [12] , [13] , [14] , [15] , [16] , [17] , similarly, bone cells are able to influence immune cells mainly through the receptor activator of nuclear factor kappa-B ligand (RANKL)–receptor activator of NF-κB (RANK)–osteoprotegerin (OPG) system [18] . Skeletal complications affect patients’ survival and quality of life [19] , thus, bone health represents a clinically relevant issue for the management of the cancer patients. Moreover, the improvement of life expectancy in patients treated with ICIs makes the attention for skeletal integrity worthy of clinical interest. All these evidence prompted us to explore the effect of immune activation, induced by ICIs, on bone remodeling. We hypothesized that ICIs treatment enhanced bone resorption through the activation of the T cells, thus increasing the risk of fractures and bone loss.
In the present study, we evaluated the changes in plasma levels of reference markers of bone turnover such as type I collagen C-terminal telopeptide (CTX-I) and N -terminal propeptide of type I procollagen (PINP) in patients affected by advanced non-small cell lung cancer (NSCLC) or renal cell carcinoma (RCC) treated with ICIs therapy. In addition, we correlated the levels of these bone turnover markers with treatment response and over-all survival (OS).
Material and methods
Study design and patient characteristics
A consecutive prospective series of 44 patients affected by advanced NSCLC and RCC treated with ICIs as monotherapy were enrolled from 2017 to 2019 and followed up until 06–2022 at Fondazione Policlinico Universitario Campus Bio Medico. Sample size estimation has been designed to achieved a power of 80 % and a level of significance 5 % (one-side) for detecting an effect size of 0.4 between CTX-I and/or PINP pairs (Time 0 (T0) vs 3 months (T1)) requiring a minimum number of 39 patients. The study was conducted in accordance with the principles of the Helsinki Declaration. All experimental protocols were approved by the Internal Review and Ethics Boards of the Campus Bio Medico University Hospital of Rome (Prot. N. 48.17OSS) and all patients provided informed consent. The inclusion criteria were patients who were at least 18 years old, with a performance status of 0–1, no signs of active autoimmune disease, without presence of bone metastases and/or previous osteoporotic fractures. We also excluded patients with conditions affecting bone, vitamin D and/or calcium metabolism (chronic liver disease, paget's disease of bone, renal failure, malabsorption, hypercortisolism); medications altering bone metabolism (e.g. denosumab, bisphosphonates, teriparatide, glucocorticoids, aromatase inhibitors, estrogen). Plasma samples were collected at the day of the first cycle of treatment be-fore the infusion and after 3 months of treatment. The patients' disease had to be measurable per Response Evaluation Criteria in Solid Tumors (RECIST) 1.1 at baseline and had to be periodically evaluated for response to treatment by Computed Tomography Scan (CT scan). Time to Treatment Failure (TTF) was defined as the interval from initiation of ICIs to its discontinuation for any reason, including cancer progression, adverse events, patient choice, or death; Overall Survival (OS) was defined as the time from the first infusion of ICIs to death or last news.
CTX-I and PINP assessment
All samples were collected in the morning hours in the fasted state. Samples were centrifugated 1500 rpm for 15 min at room temperature and serum was isolated within two hours, aliquoted and stored at −80 °C. CTX-I and PINP were quantified by Human CTX-I ELISA Kit (Cusabio Technology LLC, Huston, TX) and Human PINP ELISA Kit (Cusabio Technology LLC, Huston, TX), respectively, according to the manufacturers’ instructions.
CTX-I and PINP assessment
Descriptive statistics were reported and compared using Wilcoxon’s rank to assess CTX-I and PINP changes. CTX-I and PINP values are presented as mean ± standard error (SE). Non-parametric Spearman test was used for correlation analyses. Survival curves were estimated by the Kaplan–Meier method and compared with the Log-rank test. All statistical tests were two-sided with a P value of ≤ 0.05 considered statistically significant. All statistical analyses were performed using SPSS version 20.0 (SPSS, Chicago, IL) and R version 3.5.0 (R Institute for Statistical Computing, Vienna, Austria).
Results
Clinic-pathological findings of the patient population
From 2017 to 2019, a consecutive series of 44 patients treated with anti–PD1 therapy as monotherapy was prospectively enrolled including thirty-six patients (82 %) affected by advanced NSCLC and eight patients (18 %) suffering from advanced RCC. Twenty-three patients (52 %) received nivolumab, sixteen patients (36 %) were treated with pembrolizumab, and five patients (12 %) undergone to atezolizumab. Seventeen patients (39 %) received an anti-PD1 therapy as first-line therapy and 27 (61 %) as second or third line. The median follow up was 45 months (95 % CI 42– non-reached). The median age was 70 years (range 24–86), and female patients, all in post-menopausal state, were 19 (43 %), Patient characteristics at the study entry are summarized in Table 1 . Table 1 Clinic-pathological variables of study population. Overall (N = 44) Sex Female 19 (43 %) Male 25 (57 %) Age ≥ 70 22 (50 %) < 70 22 (50 %) ECOG 0 29 (66 %) 1 15 (34 %) Tumor type NSCLC 36 (82 %) RCC 8 (18 %) Treatment Nivolumab 23 (52 %) Pembrolizumab 16 (36 %) Atezolizumab 5 (12 %) Treatment Line First line 17 (39 %) Second or Third 27 (61 %) Primary Tumor Non resected 19 (43 %) Resected 25 (57 %) N. metastatic sites 1 15 (34 %) 2 17 (39 %) ˃2 12 (27 %) ECOG (Eastern Cooperative Oncology Group); NSCLC (Non-Small Cell Lung Cancer); RCC (Renal Cell Carcinoma); N. (number).
Table 1
Clinic-pathological variables of study population.
ECOG (Eastern Cooperative Oncology Group); NSCLC (Non-Small Cell Lung Cancer); RCC (Renal Cell Carcinoma); N. (number).
| Overall (N = 44) | |
| Sex | |
| Female | 19 (43 %) |
| Male | 25 (57 %) |
| Age | |
| ≥ 70 | 22 (50 %) |
| < 70 | 22 (50 %) |
| ECOG | |
| 0 | 29 (66 %) |
| 1 | 15 (34 %) |
| Tumor type | |
| NSCLC | 36 (82 %) |
| RCC | 8 (18 %) |
| Treatment | |
| Nivolumab | 23 (52 %) |
| Pembrolizumab | 16 (36 %) |
| Atezolizumab | 5 (12 %) |
| Treatment Line | |
| First line | 17 (39 %) |
| Second or Third | 27 (61 %) |
| Primary Tumor | |
| Non resected | 19 (43 %) |
| Resected | 25 (57 %) |
| N. metastatic sites | |
| 1 | 15 (34 %) |
| 2 | 17 (39 %) |
| ˃2 | 12 (27 %) |
CTX-I and PINP changes during ICIs
PINP and CTX-I were evaluated at both baseline (T0) and after 3 months (T1). We found a significant increase of CTX-I levels after 3 months of ICIs treatment (mean T0 0.388 ± 0.037 µg/L vs T1 0.455 ± 0.042 µg/L; p = 0.045), while we observed a decreasing trend towards the reduction of PINP levels (mean T0 3432 ± 123 pg/ml vs T1 3256 ± 169 pg/ml; p = 0.073) ( Fig. 1 A and B). In particular, twenty-three patients (52 %) showed a CTX-I beyond least significant increase and seven patients (16 %) experienced a CTX-I above least significant reduction [20] . Fig. 1 CTX-I and PINP changes. Upper panel . CTX-I (A) and PINP (B) values at baseline and after 3 months of ICIs treatment. Lower panel . Rain-cloud difference plot showing ΔCTX-I (A) and ΔPINP (B) calculated for each patients.

Correlation of CTX-I and PINP changes with clinical variables
Correlation with baseline CTX-I and PINP levels with clinically relevant variables (sex, age, ECOG PS, line of treatment, tumor and treatment types) showed no significant association (p > 0.05 for all).
As expected, low baseline CTX-I levels were observed in patients with the overweight range of Body Mass Index (BMI) (data not shown). Similarly, no significant correlations were observed when CTX-I and PINP changes (T0 vs T1) were associated with all clinic-pathological parameters including BMI.
Interestingly, CTX-I increase was associated with poor prognosis in terms of TTF (p = 0.029) and OS (p = 0.030) ( Fig. 2 ). In patients with or without increased CTX-I, median TTF was 5 months (95 % CI 4.3 – 17.1) vs 12 months (95 % CI 8.0 – non-reached), respectively ( Fig. 2 A). Similarly, median OS was lower in patients with increased CTX-I (11 months (95 % CI 6.0 – 27.0) vs 37 months (95 % CI 13.0 – non-reached) ( Fig. 2 B). Moreover, ΔCTX-I (CTX-I T1 – CTX-I T0) positively correlated with tumor size changes at best response (r = − 0.36, p = 0.018) showing that a CTX-I increase was significantly associated with an increased tumor size during ICIs ( Fig. 3 ). An additional survival analysis showed that patients with a CTX-I above “least significant increase” have a greater TTF compared to patients with a CTX-I above “least significant decrease” ( Supplementary Fig. 1 ). No significant association was found between PINP and any clinical parameter of outcome. Fig. 2 Survival analysis . Kaplan-Meier curves reporting the TTF (A) and OS (B) of patients who showed an increased CTX-I after 3 months of ICIs treatment. Fig. 3 CTX-I and tumor size changes. Scatter Plot representing the correlation between ΔCTX-I and tumor size change at best response.


Discussion
To our knowledge, this is the first study evaluating the bone remodeling status during ICIs using CTX-I and PINP, useful non-invasive biomarkers for monitoring bone health and the efficacy of anti-resorptive and anabolic therapies [21] , [22] , [23] . In our cohort, we found a significant increase of CTX-I levels, after 3 months of treatment, and a concomitantly PINP reduction with a trend towards statistical significance. Interestingly, four patients (9 %) developed new lumbar fractures after 3 months of therapy suggesting a potential relationship between ICIs, increased osteoclast activity and fracture risk. To date, how ICIs influence bone metabolism is poorly investigated, even if the established role of activated T cells in skeletal remodeling makes plausible a direct relationship between ICIs therapy and bone health. Studies from autoimmune/inflammatory diseases associated to bone loss, such as rheumatoid arthritis and postmenopausal osteoporosis showed that the increased levels of pro-inflammatory cytokines including interferon-γ (IFN-γ), tumor necrosis factor-α (TNF-α) and interleukin-6 (IL-6) levels stimulated osteoclast differentiation and activity [12] , [24] , [25] . Moreover, in these inflammatory diseases, T -cell activation induced RANK-L production leading to a disruption of the RANK-L/OPG balance, increasing bone turnover and promoting bone loss [26] , [27] , [28] . Similarly, ICIs therapy induce the secretion of pro-inflammatory cytokines that are implicated in tumor destruction but, at the same time, could adversely affect bone remodeling stimulating osteoclast function, thus resulting in an increased risk of fractures. These evidences lead us to sup-pose that there is potential relationship between ICIs treatment, increased CTX-I and the development of fractures observed in our cohort.
Additional research is warranted to elucidate the significant correlation between increased CTX-I and poor prognosis, in terms of survival and treatment response. One potential explanation is that osteoclasts could create a systemic immunosuppressive microenvironment that promotes tumor growth. Evidence from multiple myeloma support the immunosuppressive role of osteoclasts that.
directly inhibit proliferation of activated CD4 + and CD8 + effector T cells. In addition, the secretion of some molecules by osteoclasts, such as galectin-9 and a proliferation-inducing ligand (APRIL), induced T cells apoptosis and enhanced tumor growth [29] .
It has widely described that Th17 cells stimulate bone resorption directly through the expression of RANKL on their surface or, indirectly, by the secretion of IL-17 that enhances the expression of RANKL in osteoblasts [30] , [31] . A recent paper demonstrated that an increased infiltration of T helper 17 (Th17) cells caused primary resistance to ICIs in lung cancer patients [32] . Thus, another possible scenario could be that specific tumor-infiltrated T -cells in tumors resistant to ICIs secrete pro-osteoclastogenic cytokines. As results of an excessive osteoclast activation, we observed a CTX-I increase in primary refractory patients. Based on all these evidences, it reasonable to hypothesize a bidirectional interaction between osteoclasts and specific populations of tumor-infiltrated T -cell, such as Th17 cells, which creates a vicious cycle that promote cancer cell progression and bone osteolysis. Finally, we cannot exclude that the bone turnover increase in patients who experienced disease progression is associated to a general physical health deterioration that affect also bone tissue.
The results of this study are in accordance with our previous paper that demonstrated the correlation between the high levels of circulating RANKL and poor prognosis in patients with metastatic RCC treated with ICIs [33] .
Our exploratory study includes a small, but adequately powered number of patients with stringent exclusion criteria that makes the data reliable and unaffected by potential biases or confounding factors. The major limitation of the study is the lack of a control group of patients who were not treated with ICIs. The adequate control group would have been patients mainly treated with chemotherapy, but it is known that anti-neoplastic agents and steroids used to prevent emesis both affect bone homeostasis creating potentially confounding elements in the analyses. Although we can not rule out that CTX-I increase is independent from ICIs therapy, the association between CTX-I changes and ICIs response make less plausible this second hypothesis. The study is limited also by immortal time biases since patients who died before 3 months of ICIs treatment were not included. Furthermore, bone turnover markers were evaluated only at baseline and after 3 months without a longitudinal assessment of their dynamics over time in long-responder patients.
Overall, the key message of this study is that ICIs may act as triggers for skeletal events, but larger prospective studies are needed to confirm if this treatment really affects the bone. In this regard, laboratory and imaging studies should be performed in patients who are candidates to ICIs treatment for risk stratification and early detection of porotic fractures. In addition, the presence of risk factors for impairments of bone metabolism/fracture risk factors such as pre-existing osteoporosis, fragility fractures, concomitant inflammatory arthritis or other autoimmune disease, genetic or environmental factors should be taken in consideration before starting treatment. As future clinical perspective, the use of anti-osteoclastic drugs such as bisphosphonate and denosumab could represent a novel combination strategy with ICIs not only to prevent skeletal lesions but also to improve patients’ clinical outcome. In this regards, recent studies evaluating the combination efficacy of anti-osteoclastic agents and ICIs have already shown promising results [34] , [35] , [36] , [37] .
Ethics approval and consent to participate .
The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board of the Fondazione Policlinico Universitario Campus Bio Medico (Prot. N. 48.17OSS).
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
References
1 Emens L.A. Ascierto P.A. Darcy P.K. Demaria S. Eggermont A.M.M. Redmond W.L. Seliger B. Marincola FM Cancer immunotherapy: Opportunities and challenges in the rapidly evolving clinical landscape Eur J Cancer. 81 2017 116 129 28623775 10.1016/j.ejca.2017.01.035
2 Gandhi L. Rodríguez-Abreu D. Gadgeel S. Esteban E. Felip E. De Angelis F. Pembrolizumab plus Chemotherapy in Metastatic Non-Small-Cell Lung Cancer N Engl J Med 378 2018 2078 2092 29658856 10.1056/NEJMoa1801005
3 Motzer R.J. Escudier B. McDermott D.F. George S. Hammers H.J. Srinivas S et al Nivolumab versus Everolimus in Advanced Renal-Cell Carcinoma N Engl J Med 373 2015 1803 1813 26406148 10.1056/NEJMoa1510665 PMC5719487
4 Hamid O. Robert C. Daud A. Five-year survival outcomes for patients with advanced melanoma treated with pembrolizumab in KEYNOTE-001 Ann Oncol 30 2019 582 588 30715153 10.1093/annonc/mdz011 PMC6503622
5 Postow M.A. Sidlow R. Hellmann MD Immune-related adverse events associated with immune checkpoint blockade N Engl J Med 378 2018 158 168 29320654 10.1056/NEJMra1703481
6 Vozy A. De Martin E. Johnson D.B. Increased reporting of fatal hepatitis associated with immune checkpoint inhibitors Eur J Cancer 123 2019 112 115 31678768 10.1016/j.ejca.2019.09.022
7 Guerrero E. Johnson D.B. Bachelot A. Immune checkpoint inhibitor-associated hypophysitis World Health Organisation VigiBase report analysis Eur J Cancer 113 2019 10 13 30954881 10.1016/j.ejca.2019.03.002
8 Raschi E. Mazzarella A. Antonazzo I.C. Toxicities with Immune Checkpoint Inhibitors: Emerging Priorities From Disproportionality Analysis of the FDA Ad-verse Event Reporting System Target Oncol 14 2019 205 221 30927173 10.1007/s11523-019-00632-w
9 Raschi E. Antonazzo I.C. La Placa M. Ardizzoni A. Poluzzi E. De Ponti F Serious cutaneous toxicities with immune checkpoint inhibitors in the U.S. Food and Drug Administration adverse event reporting system Oncologist 11 2019 e1228 e1231 10.1634/theoncologist.2019-0250 PMC6853099 31387950
10 Moseley K.F. Naidoo J. Bingham C.O. Immune-related adverse events with immune checkpoint inhibitors affecting the skeleton: a seminal case series J Immunother Cancer 6 2018 104 30305172 10.1186/s40425-018-0417-8 PMC6180387
11 Filippini D.M. Gatti M. Di Martino V. Bone fracture as a novel immune-related adverse event with immune checkpoint inhibitors: Case series and large-scale pharmacovigilance analysis Int J Cancer 149 2021 675 683 33844854 10.1002/ijc.33592 PMC8251715
12 Pietschmann P. Mechtcheriakova D. Meshcheryakova A. Immunology of osteoporosis: a mini-review Gerontology 62 2016 128 137 26088283 10.1159/000431091 PMC4821368
13 D'Amelio P. Grimaldi A. Di Bella S. Estrogen deficiency increases osteoclastogenesis up-regulating T cells activity: a key mechanism in osteoporosis Bone 43 2008 92 100 18407820 10.1016/j.bone.2008.02.017
14 Kawai T. Matsuyama T. Hosokawa Y. B and T lymphocytes are the primary sources of RANKL in the bone resorptive lesion of periodontal disease Am J Pathol 169 2006 987 998 16936272 10.2353/ajpath.2006.060180 PMC1698808
15 Kotake S. Nanke Y. Mogi M. IFN-gamma-producing human T cells directly induce osteoclastogenesis from human monocytes via the expression of RANKL Eur J Immunol 35 2005 3353 3363 16220542 10.1002/eji.200526141
16 Li Y. Toraldo G. Li A. B cells and T cells are critical for the preservation of bone homeostasis and attainment of peak bone mass in vivo Blood 109 2007 3839 3848 17202317 10.1182/blood-2006-07-037994 PMC1874582
17 Sato K. Suematsu A. Okamoto K. Th17 functions as an osteoclastogenic helper T cell subset that links T cell activation and bone destruction J Exp Med 203 2006 2673 2682 17088434 10.1084/jem.20061775 PMC2118166
18 Ponzetti M. Rucci N Updates on Osteoimmunology: What's New on the Cross-Talk Between Bone and Immune Sys-tem Front Endocrinol (Lausanne) 10 2019 236 31057482 10.3389/fendo.2019.00236 PMC6482259
19 von Moos R. Costa L. Gonzalez-Suarez E. Management of bone health in solid tumours: From bisphosphonates to a monoclonal antibody Cancer Treat Rev 76 2019 57 67 31136850 10.1016/j.ctrv.2019.05.003
20 Claudon A. Vergnaud P. Valverde C. New automated multiplex assay for bone turnover markers in osteoporosis Clin Chem 54 2008 1554 1563 18676587 10.1373/clinchem.2008.105866
21 Vilaca T. Gossiel F. Eastell R Bone Turnover Markers: Use in Fracture Prediction J Clin Densitom 20 2017 346 352 28716498 10.1016/j.jocd.2017.06.020
22 Bauer D. Lane N. Leary E. National Bone Health Alliance Bone Turnover Marker Project: current practices and the need for US harmonization, standardization, and common reference ranges Osteoporos Int. 23 2012 2425 2433 22797491 10.1007/s00198-012-2049-z PMC4011662
23 Dai Z. Wang R. Ang L.W. Bone turnover biomarkers and risk of osteoporotic hip fracture in an Asian population Bone 83 2016 171 177 26555636 10.1016/j.bone.2015.11.005 PMC4724247
24 Cauley J.A. Barbour K.E. Harrison S.L. Inflammatory markers and the risk of hip and vertebral fractures in men: the osteoporotic fractures in men (MrOS) J Bone Miner Res 31 2016 2129 2138 27371811 10.1002/jbmr.2905 PMC5240475
25 Walsh M.C. Takegahara N. Kim H. Updating osteoimmunology: regulation of bone cells by innate and adaptive immunity Nat Rev Rheumatol 14 2018 146 156 29323344 10.1038/nrrheum.2017.213 PMC5821527
26 Weitzmann M.N. The role of inflammatory cytokines, the RANKL/OPG Axis, and the Immunoskeletal Interface in physiological bone turnover and osteoporosis Scientifica (Cairo) 2013 2013 125705 10.1155/2013/125705 PMC3820310 24278766
27 Nakashima T. Takayanagi H RANKL signal and osteoimmunology Clin Calcium 21 2011 1131 1140 21814017
28 Showalter A. Limaye A. Oyer J.L. Cytokines in immunogenic cell death: applications for cancer immunotherapy Cytokine 97 2017 123 132 28648866 10.1016/j.cyto.2017.05.024 PMC5572581
29 Tai Y.T. Cho S.F. Anderson KC Osteoclast Immunosuppressive Effects in Multiple Myeloma: Role of Programmed Cell Death Ligand 1 Front Immunol. 9 2018 1822 30147691 10.3389/fimmu.2018.01822 PMC6095980
30 Ono T. Takayanagi H Osteoimmunology in bone fracture healing Curr Osteoporos Rep 15 2017 367 375 28647888 10.1007/s11914-017-0381-0
31 Zhu L. Hua F. Ding W. The correlation between the Th17/Treg cell balance and bone health Immun Ageing 17 2020 30 33072163 10.1186/s12979-020-00202-z PMC7557094
32 Peng D.H. Rodriguez B.L. Diao L. Th17 cells contribute to combination MEK inhibitor and anti-PD-L1 therapy resistance in KRAS/p53 mutant lung cancers Nat Commun 12 2021 2606 33972557 10.1038/s41467-021-22875-w PMC8110980
33 Simonetti S. Iuliani M. Stellato M. Extensive plasma proteomic profiling revealed receptor activator of nuclear factor kappa-Β ligand (RANKL) as emerging biomarker of nivolumab clinical benefit in patients with metastatic renal cell carcinoma J Immunother Cancer 10 2022 e005136 36104102 10.1136/jitc-2022-005136 PMC9476128
34 Dougall W.C. Roman Aguilera A. Smyth MJ Dual targeting of RANKL and PD-1 with a bispecific antibody improves anti-tumor immunity Clin Transl Immunol 8 2019 e01081 10.1002/cti2.1081 PMC6763724 31572609
35 Ahern E. Harjunpää H. O'Donnell JS, et al RANKL blockade improves efficacy of PD1-PD-L1 blockade or dual PD1-PD-L1 and CTLA4 blockade in mouse models of cancer Oncoimmunology 7 2018 e1431088 29872559 10.1080/2162402X.2018.1431088 PMC5980384
36 Smyth M. Yagita H. McArthur GA Combination Anti-CTLA-4 and Anti-RANKL in Metastatic Melanoma J Clin Oncol 34 2016 e104 e106 25311221 10.1200/JCO.2013.51.3572
37 Bostwick A.D. Salama A.K. Hanks BA Rapid complete response of metastatic melanoma in a patient undergoing ipilimumab immunotherapy in the setting of active ulcerative colitis J Immunother Cancer 3 2015 19 25992290 10.1186/s40425-015-0064-2 PMC4437559
Supplementary data
The following are the Supplementary data to this article: Supplementary data 1
Supplementary data to this article can be found online at https://doi.org/10.1016/j.jbo.2022.100459 .