Article Text
PDF
Abstract
Background Nivolumab is a human monoclonal antibody against programmed cell death receptor-1 (PD-1) able to rescue quiescent tumor infiltrating cytotoxic T lymphocytes (CTLs) restoring their ability to kill target cells expressing specific tumor antigen-derived epitope peptides bound to homologue human leukocyte antigen (HLA) molecules. Nivolumab is currently an active but expensive therapeutic agent for metastatic non-small cell lung cancer (mNSCLC), producing, in some cases, immune-related adverse events (irAEs). At the present, no reliable biomarkers have been validated to predict either treatment response or adverse events in treated patients.
Methods We performed a retrospective multi-institutional analysis including 119 patients with mNSCLC who received PD-1 blockade since November 2015 to investigate the predictive role of germinal class I HLA and DRB1 genotype. We investigated the correlation among patients’ outcome and irAEs frequency with specific HLA A, B, C and DRB1 alleles by reverse sequence-specific oligonucleotide (SSO) DNA typing.
Results A poor outcome in patients negative for the expression of two most frequent HLA-A alleles was detected (HLA: HLA-A*01 and or A*02; progression-free survival (PFS): 7.5 (2.8 to 12.2) vs 15.9 (0 to 39.2) months, p=0.01). In particular, HLA-A*01-positive patients showed a prolonged PFS of 22.6 (10.2 to 35.0) and overall survival (OS) of 30.8 (7.7 to 53.9) months, respectively. We also reported that HLA-A and DRB1 locus heterozygosis (het) were correlated to a worse OS if we considered het in the locus A; in reverse, long survival was correlated to het in DRB1.
Conclusions This study demonstrate that class I and II HLA allele characterization to define tumor immunogenicity has relevant implications in predicting nivolumab efficacy in mNSCLC and provide the rationale for further prospective trials of cancer immunotherapy.
- B7-H1 antigen
- antigen presentation
- lung neoplasms
- tumor biomarkers
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See http://creativecommons.org/licenses/by-nc/4.0/.
Statistics from Altmetric.com
Background
Non-small cell lung cancer (NSCLC) is the most common malignancy and the leading cause of cancer death. Chemotherapy with platinum doublets±bevacizumab, a monoclonal antibody (mAb) to the vascular endothelial growth factor (reserved for patients with non-squamous histology and low risk of bleeding only), has been for many years the standard treatment for patients in metastatic (m) and advanced stage of disease (stage IIIB-IV), a good performance status (Eastern Cooperative Oncology Group (ECOG) 0–1) and no druggable driver mutations/rearrangements of epidermal growth factor receptor (EGFR), anaplastic lymphoma kinase/echinoderm microtubule-associated protein-like 4 (ALK/EML4) and ROS proto-oncogene 1 (ROS1).1 2
In the last few years, effector immune checkpoint blockade (ICB) with mAbs to programmed cell death receptor-1 (PD-1; nivolumab and pembrolizumab) and PD-1 ligand-1 (PD-L1; atezolizumab, avelumab and durvalumab) alone or in combination with chemotherapy, radiotherapy or bevacizumab has grown as a solid treatment for these patients with remarkable results in terms of clinical benefit, progression-free survival (PFS), and overall survival (OS) which has been observed in more than 30% of the cases.3–5 On the other hand, these immune-oncological treatments may be associated to more or less severe immune-related adverse events (irAEs) and to high costs. It is, therefore, conceivable that reliable biomarkers predictive of response and irAEs are eagerly awaited.6–8
On this basis, we led a real-world multicenter retrospective investigation aimed to evaluate the predictive role of different human leukocyte antigen (HLA) alleles in patients undergone PD-1 blockade with nivolumab by taking in consideration that: (1) these mAbs promote the rescue of cytotoxic T lymphocytes (CTLs) function silenced by the PD-1/PD-L1 or PD-L2 immune checkpoint at tumor sites,8 (2) HLA molecules play a critical role in triggering CTL-mediated tumor cell killing, T cell priming and clonal expansion,9–12 and finally, (3) there is a strong correlation between expression of specific HLA alleles and development of autoimmune disease.13
Here, we correlated the germline expression of the most frequent HLA alleles in the loci A, B, C and DRB1 (reflecting a β subunit of class II HLA-DR) with outcome and risk of iraEs in patients with mNSCLC undergone to salvage therapy with nivolumab treatment.
Methods
Patients sample, treatment and monitoring
This study is part of a retrospective real-world evidence (RWE) multi-institutional study including 119 pretreated patients with mNSCLC who had received salvage therapy with nivolumab at the Medical Oncology Unit, Grand Metropolitan Hospital “Bianchi-Melacrino-Morelli”, Reggio Calabria, Medical and Translational Oncology Unit, Magna Graecia University, Catanzaro, and Section of Radiation Oncology, University of Siena, Siena, between September 2015 and December 2018. This study included patients who had been voluntary screened for A, B, C and DRB1 HLA typing at the Tissue Typing Unit at the Grand Metropolitan Hospital in Reggio Calabria.
All patients gave an informed consent to anonymous use of their examinations for research aim. All patients received PD-1 blockade with intravenous nivolumab at the dosage of 3 mg/kg in 60 min on a biweekly schedule until progression, grade 3–4 adverse events or death. All patients had received at least a previous line of platinum-based doublet±bevacizumab prior PD-1 blockade, presented with ECOG performance status ≤1, have had complete physical examination reports, histological sampling and hematological, biochemical, immune-biological, radiological and instrumental monitoring at baseline. Clinical history, physical examination and record of adverse events were evaluated prior each drug infusion. A CT scan was performed every 3 months or in any case of suspected progressive disease and evaluated according to the immune response evaluation criteria in solid tumors (iRECIST).14 Patients were monitored for blood cell counts and biochemistry prior each treatment course and were also monitored for adrenal hormone profile, adrenocorticotropic hormone (ACTH), thyroid stimulating hormone (TSH), thyroid hormones, antithyroid autoantibodies (AAbs), extractable nuclear antigen antibodies, antinucleus antibodies, antismooth cells antibodies and c/p- antineutrophil cytoplasmic antibodies on monthly bases since the beginning of treatment.
HLA molecular typing of the loci A, B, C, and DRB1 was centralized and carried out in the Tissue Typing Unit at the Grand Metropolitan Hospital in Reggio Calabria by performing low-medium-resolution reverse SSO DNA typing assays (One-Lambda Luminex Technology LABScan 100, HLA Fusion Software) on genomic DNA extracted by whole blood or peripheral blood mononuclear cells according to the kit manufactures. The HLA allele frequency of the studied loci A, B, C and DRB1 in our patients’ cohorts was obtained by comparing the results with our database bank including 3500 healthy bone marrow donors of the Calabrian National Registry.
Survival analysis
PFS and OS were analyzed with the Kaplan-Meier method and significance of survival differences were analyzed by the log-rank test. Median survival and 95% CIs were reported. HRs and their 95% CIs were estimated through the Cox regression proportional model; in the multivariate approach, a forward stepwise procedure was used and enter and remove limit set to 0.05 and 0.10, respectively. Association between irAEs and clinical and biological factors was assessed by the χ2 test. Our statistical analyses were performed by the SPSS software V.23.0.
Results
Patients’ features
Our analysis included a cohort of 119 patients with mNSCLC, 99 males and 20 females with median age of 68 years. Forty-two presented with a squamous and 77 with a non-squamous histology. Ninety-one patients had received a front-line doublet chemotherapy±bevacizumab prior immunotherapy, 27 a platinum-based metronomic chemotherapy and 1 did not receive chemotherapy. Fifty-seven patients had also received palliative radiotherapy (25–30 Gy) on symptomatic single lesions (mediastinum, nodes, soft tissue, bone or brain) prior or at beginning of nivolumab treatment.
Predictive values of clinical parameters
Patients included in the retrospective study presented a median follow-up of 22 months with a median PFS and OS of 9.7 (6.1 to 13.3) and 18.5 (14.6 to 22.4) months, respectively (figure 1A,B).
Progression-free survival (PFS) and overall survival (OS). (A, B) PSF and OS of patients with metastatic non-small cell lung cancer (mNSCLC) subjected to nivolumab treatment. (C, D) PFS and OS of the same patients with mNSCLC with and without germinal expression of human leukocyte antigen (HLA)-A*01 and or A*02. (E, F) PFS and OS of patients with mNSCLC with and without germinal heterozygosis in class I HLA-A locus. (G, H) PFS and OS of patients with mNSCLC with and without germinal heterozygosis in DRB1 locus. (#): differences statistically significant p<0.05.
IrAEs were reported in 59 patients and mainly consisted in grade 1–2 cutaneous rash, polyarthritis and thyroiditis, generally occurring after 3–4 treatment courses. Autoimmune pneumonitis occurred in five patients among the third and sixth month of treatment and additional three cases of uveitis recorded after 2 and 3 months of treatment.
HLA typing
In this patient population we detected the two most frequent germline alleles for each class I HLA and DRB1 locus and, then, we correlated their expression with both outcome and irAEs frequency. The most frequent alleles resulted A*02 (38 pts) and A*01 (25 pts) for the locus A; B*18 (20 pts) and B*35 (20 pts) for the locus B; C*07 (43 pts) and C*04 (27 pts) for the locus C and DRB1*11 (36 pts) and DRB1*03 (20 pts) for DRB1. In our analysis, we also included allele HLA-B*51 (17 pts) for its known correlation with common autoimmune diseases.15 Allele homozygosis for the loci A, B, C and DRB1 was detected in 10, 4, 11 and 13 patients, respectively. These results were in line with the HLA frequency recorded in healthy bone marrow donors in the Calabrian Regional Registry, representative of the healthy population in our geographical area.
Statistical correlations: patients’ outcome and the most frequent alleles for HLA-A, B, C and DRB1
A poor outcome treatment response was observed in patients not expressing either HLA-A*02 and/or HLA-A*01 alleles (PFS: 7.5 (2.8 to 12.2) vs 15.9 (0 to 39.2) months, p=0.01 with a trend to worse OS of 13.4 (95% CI not evaluable) months vs 20.4 (6.1 to 34.7) and months, p=0.1; figure 1C,D). In this context, HLA-A*01 positive patients showed the best outcome with a PFS and OS of 22.6 (10.2 to 35.0) and 30.8 (7.7 to 53.9) months, respectively. We also found a trend to a worse outcome in patients expressing the HLA-C*07 allele (C*07 positive vs negative; PFS=12.2 (9.3 to 15.1) vs 29.1 (7.2 to 51.0) months p=0.09). Percentage of 1-year PFS rate and 1/2-year survival rate presented in table 1 showed valuable, though not statistically significant, difference among the groups reporting outstanding 1/2-year survivals in patients showing HLA-A*01 (87.5% /55.9%), HLA-B*51 (82.4%/74.1%) and DRB1*03 (83.8%/58.7%) alleles.
Clinical features, 1-year progression-free survival rate and 1/2-year survival rate for each patient cohort expressing the most frequent alleles in class I human leukocyte antigen (HLA) A, B and C and DRB1
Statistical correlations: patients’ outcome and heterozygosis in HLA-A, B, C and DRB1 loci
Subsequently, we evaluated whether heterozygosis (het) in HLA-A, B, C and DRB1 loci could have a predictive value on patients’ outcome. First, we did not find any correlation between het and PFS (data not shown). Completely different results were instead recorded for what concerns the influence of het of HLA-A and DRB1 genes on patients’ OS. Het in the locus A was indeed correlated to a worse OS (no-het vs het=more than 40 (median not yet reached) vs 20.8 (11.4 to 30.2) months; p=0.03; figure 1E,F) while on the contrary, het in DRB1 was correlated to a prolonged survival (no het vs het=18.5 (4.6 to 32.4) vs 30.8 (20.4 to 41.2) months, p=0.05; figure 1G,H). Finally, patients showing het in the loci B and C showed a trend to prolonged survival that did not achieve statistical significance (data not shown). Due to the small patients’ sample, it was not possible to find any correlation among survival and heterozygosis in specific locus alleles.
We also performed an additional analysis on 375 patients with mNSCLC who received PD-1/PD-L1 blockade, extrapolated by a recent survey published by Chowell et al (NIHMS980063-table S1),16 and those included in our series, correlating their OS with heterozygosis in at least one among class I HLA-A, B or C loci (figure 2A) and those bearing selective heterozygosis in the locus A (figure 2B). In the first case, our analysis failed to demonstrate any statistically significant difference in both patients’ groups (from our study, Italy, ITA and from Chowell’s study, National Institutes of Health; figure 2A; p=0.39 and p=0.12). Similarly, we were unable to demonstrate a significant difference in survival correlated with HLA-A heterozygosis in the patients’ population examined in the Chowell’s study considering that patients not presenting HLA-A heterozygosis (53 out of 375 pts) showed a twice shorter follow-up compared with the other cohort and their median of survival was not achieved yet (figure 2B).
Survival of 375 patients with metastatic non-small cell lung cancer subjected to programmed cell death receptor-1 (PD-1)/PD-1 ligand-1 blockade included in the database published by Chowell et al (NIHMS980063-table s1; group National Institutes of Health (NIH)).16 (A) Patients presenting homozygosis in at least one class I human leukocyte antigen (HLA) locus (NIH-no-het) versus full heterozygosis (NIH-het). (B) Patients presenting homozygosis versus heterozygosis in HLA-A showing no statistically significant differences (p>0.1). Similar results were observed in our patients’ series (group-IT) for what concerns homozygosis in at least one class I HLA locus (IT-no-het) versus full heterozygosis (IT-het) while significant differences in survival were observed, when the overall survival of patients with homozygosis and heterozygosis in HLA-A locus was compared (p=0.03).
Multivariate analysis: patients’ outcome, gender, irAEs’ frequency, AAb rise, baseline inflammatory status, HLA A, B, C and DRB1 alleles and heterozygosis
Our univariate analysis showed that a prolonged PFS was correlated with occurrence of irAEs, low baseline erythrocyte sedimentation rate values, type of chemotherapy (metronomic chemotherapy±bevacizumab), presence of allele A*01 and, in general, haplotype A*01 and/or A*02. When independent factors among these were investigated, only occurrence of irAEs and expression of A*01 and/or A*02 retained significance (table 2). Differently, a prolonged OS was correlated with the presence of irAEs, low baseline C-reactive protein value, chemotherapy type (metronomic chemotherapy±bevacizumab), no heterozygosis in locus A and heterozygosis in DRB1. When independent factors were investigated, only no heterozygosis in locus A retained significance (table 3).
Univariate and multivariate analyses of progression-free survival (PFS) in relation to immune-related adverse events (irAEs), low baseline erythrocyte sedimentation rate (ESR) values, type of chemotherapy, presence of allele A*01 and haplotype A*01 and/or A*02
Univariate and multivariate analyses of overall survival (OS) in relation to immune-related adverse events (irAEs), low baseline erythrocyte sedimentation rate (ESR) values, type of chemotherapy, presence of allele A*01 and haplotype A*01 and/or A*02
Discussion
This retrospective study was carried out in patients with mNSCLC who received a real-world salvage therapy with nivolumab. The results of the study fulfilled the primary endpoint to investigate the possible correlations among patients’ outcome and irAEs frequency with the germinal expression of multiple class I HLA and DRB1 loci. Variability of class I HLA genes has been related to ICB response in a number of malignancies including melanoma and colorectal cancer.16 17 The association between autoimmunity and antitumoral effect induced by immunotherapy is revealed by IRAEs that more strongly emerge after the treatment with anti-PD-1 and anti-PD-L1 than in response to anti-cytotoxic T-lymphocyte protein 4(CTLA-4) antibody.18 Then, it could be useful to study HLA genetic variation in order to define biomarker of response to ICIs in terms of benefit and adverse events. The rationale of the study is based on the evidence that HLA molecules play a critical role in engaging an antigen-specific immune response acting both on primary cross-priming and peripheral T cell effector response. The results of the present study indeed showed that the PFS of these patients was strongly correlated to their germinal expression of HLA-A*01 and/or A*02 alleles, the most frequent class I HLA-A allele in the Caucasian population and Mediterranean geographical area.19 On the contrary, patients expressing less frequent class I A allele and not A*01 and A*02 showed a very poor prognosis. We also showed that the expression of homozygosis in the locus A (mainly represented by A*01 and A*02 alleles) correlated with their OS confirming the critical role for these molecules in promoting an efficient and prolonged anticancer immune response in patients subjected to PD-1 blockade. In our patient population, however, the best survival was observed in those patients presenting germinal heterozygosis in the locus DRB1 and, therefore, class II HLA-DRB1 molecules, which, differently from class I HLA, are not directly involved in the CTL-mediated killing of tumor cells enforced by PD-1 blockade, but rather to cross priming, T and B cell modulation20 and autoimmunity.21 A recent study showed that no significant correlations were observed for HLA class I zygosity and PFS or OS in three independent cohorts, and tumor genomic characteristics, such as EGFR and ALK alterations, should take in consideration to determine benefit from ICB in NSCLC. Moreover, the authors suggest that decreased expression of MHC class I and antigen presentation through MHC class II should be evaluated in NSCLC to predict response to ICB.22 In this context, our analysis confirmed the previous finding that occurrence of irAEs, AAbs rise, a low baseline inflammatory status and the kind of chemotherapy regimen received prior PD-1 blockade were all correlated to patients’ survival.23 Class II HLA molecules alike HLA-DRB1 bind peptides recognized by T cell receptors (TCRs) on CD4+ T cells with helper and/or T regulatory function; thus, they can modulate the differentiation and the proliferation of other lymphocyte subsets including CTLs and B cells.20 We believe that the occurrence of autoimmunity in these patients, in line with the results of other immunological studies, is consequential to the cross priming of antigens released from the tumor site in an immunological context of CTL mediated immune rejection6 24 25 whose cross priming may be facilitated and redirected to the B cell compartment by the germline expression of specific and multiple class II HLA alleles.
All together the results of our study are in line with the knowledge that mAbs to PD-1, such as nivolumab, exert their antitumor effects by reactivating tumor-specific CTLs26 whose existence largely precedes the diagnosis and the immunological treatment and that are consequence of a long-lasting immune response to multiple tumor-associated (TAA) and tumor-specific antigens (TSA).12 27
The ultimate effectors of PD-1 blockade are primed CTLs able to recognize target cells with a specific TCR which, in turn, binds TAA/TSA derived 9–10 mere-peptides (epitopes) bound to homologous class I HLA molecules on their surface. On the other hand, T cells priming and expansion occur within specialized lymphoid organs where T cell precursors use their TCRs to recognize peptide epitopes/class I–II HLA complexes on antigen-presenting cells able to uptake and process antigenic material released in the blood and lymphatic stream by the tumor tissues.9–12 Epitope peptide binding to HLA molecules occurs in the cytoplasm operated by the proteasome system, and it is constrained by definite HLA allele/locus specific amino acid consensus motifs. The latter have been largely described in the literature and specific network algorithms have already been defined to detect epitopic sequences and their potential binding affinity for specific HLA alleles by scanning the amino acid sequence of viral and tumor-derived protein antigens.28–30 A large number of dominant immunogenic epitope peptides from various TAA/TSA and viruses specific for the most frequent class I HLA alleles have already been identified using these specific algorhythms. Many of these peptides have been synthesized and tested in several studies worldwide.31 32 Some of these peptides derived from prostate-specific antigen, carcinoembryonic antigen, human telomerases and thymidilate synthase have also been evaluated in clinical trials of active specific immunotherapy in patients with cancer obtaining contrasting results.31–38 Moreover, the sequence divergence between alleles of HLA-I genotype is related to the increased probability of neoantigen recognition by tumor-infiltrating T cells and subsequently ability of different individual immune system response to checkpoint blockade immunotherapy.39 The main lesson from these studies was that protein antigens present one or more epitope sequence with class I/II HLA allele-specific binding motifs and that both their number and binding affinity for different HLA alleles and TCRs are not uniform. It is consequential that patients with different class I and II HLA germinal haplotypes present divergent immune reactivity to the tumor due to a completely different pool of epitopes available for immune priming, CTL recognition and tumor cell killing even in the presence of the same target antigens. Several studies have already shown that gene and/or upstream downregulation and/or one or more HLA allele and antigen-peptide processing in tumor cells are effective mechanisms of tumor immune escape and resistance to immune checkpoint in patients with cancer.40 41 Another study supports the hypothesis that an allelic variant of killer cell immunoglobulin-like receptor (KIR), KIR3DS1, mainly expressed on natural killer (NK) cells, is associated to resistance to PD-1 blockade in lung cancer.42 The recent results of a large survey by Chowell et al, including a pool of 1535 patients with advanced cancer receiving CTLA-4 and/or PD-1/PD-L1 blockade, showed that germinal homozygosis in at least one of the three class I HLA locus (A, B and C) were predictive of a poor survival.16 These results were in line with our findings concerning heterozygosis in the B, C and DRB1 loci, whereas apparently in contrast with the finding that class I HLA-A locus heterozygosis was correlated to a worse outcome in our series. This difference could be simply explained considering that they do not focus on class I HLA-A and on the other loci separately; additionally, their analysis includes patients with various malignancies with different prognosis and treatment responsiveness and finally, they pool together, patients who received anti-CTLA-4 and anti-PD-1/PD-L1 immune checkpoints inhibitors that act on different phases of the immune response and present a non-homogenous range of activity. Therefore, we performed an additional analysis on the 375 patients with mNSCLC who received PD-1 blockade and included in Chowell’s survey (supplemental NIHMS980063_table_S1)16 finding no outcome differences correlated with class I A, B and C allele heterozygosis. In the same patients’ series, we were unable to extrapolate significant data on the outcome of patients with heterozygosis in the locus A, considering the following: (1) patients with HLA-A locus homozygosis (53 out of 375) showed a twice shorter follow-up; (2) they have not achieved the median survival when the data were collected (before January 2018); (3) the two groups had not received a uniform immunological treatment.
Conclusions
All together these findings suggest that the different response to nivolumab observed in patients with less frequent class I HLA-A alleles might be a consequence of a divergent CTL response, whose entity is dependent on the different amount of high affinity epitope peptides (epitope repertoire) bound to their specific class I HLA profile on the locus A. It is, therefore, possible to speculate that patients with a germline expression of the most frequent class I HLA-A alleles in a specific geographical area (A*02 and/or A*01) also present a more appropriate CTL epitope peptide repertoire aimed to protect them in a specific environment.
If these data should be confirmed on a large series, it is reasonable to propose the combined use of ICB and allele-specific anticancer vaccines in patients who failed to rise a spontaneous tumor antigen-specific immune response. Recent findings concerning the application of PD-1/PD-L1 blockers in combination with front-line chemotherapy in patients with advanced NSCLC compared with the same chemotherapy doublet alone were partially sorrowed by the occurrence of a high risk of early deaths (within the first 60 days of treatment) due to severe immunological disorders. Even though, it is possible to speculate that immune system over-stimulation combining danger signals and massive antigen release consequent to secondary cytotoxic agents (anticancer drugs/radiotherapy/antiangiogenic agents) might be detrimental in patients who already show an inappropriate immune response.43 44 On these bases, we believe that an accurate multimodal definition of tumor immunogenicity at baseline as well as a germline class I and II HLA allele characterization deserves to be studied on a large scale in these patients in a contest of RWE in order to define not only responsive patients but also patients with high risk of lethal irAEs and to propose new treatment modalities combining immune checkpoint inhibitors and class I HLA allele-specific anticancer vaccines.
Acknowledgments
The authors wish to thank patients and their families for their participation to the study and the ARCO patients’ association for the precious fundraising support to the study.
References
Footnotes
Contributors PC, DG, MC, PT and PT prepared the manuscript and the figures. PC, RG, RA, AS, AF, PP, PT, PT and PT clinically evaluated patients enrolled for the study. RES, NI, VA and GR performed the HLA typing. DG performed statistical analysis. NS, TDG, DC, VB, VN, PP and PT contributed to the experimental design. MA, AG, AN, EA and LP improved the discussion of the results. AL, MC and PC revised the manuscript. All authors read and approved the final manuscript.
Funding This work has been supported by the NSP project nos. LO1508 and LO1309, MZ-VES project nos. 16-28637A, 16-2960A, and 17-32285A, and by funds institutional research (TA 29) of UVPS Brno.
Competing interests None declared.
Patient consent for publication Obtained.
Ethics approval All procedures were undertaken in compliance with the ethical statements of the Helsinki Declaration (1964, amended most recently in 2008) of the World Medical Association and respect of their privacy. The protocol and any amendments were reviewed and approved by an Institutional Review Board/Independent Ethics Committee prior to initiation of the study. All patients provided voluntary written informed consent.
Provenance and peer review Not commissioned; externally peer reviewed.
Data availability statement The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.