BIOCELL DOI:10.32604/biocell.2022.019724
Article

A potential impact of A Disintegrin and Metalloproteinase Domain-Like Protein Decysin-1 (ADAMDEC1) on clear cell renal cell carcinoma propagation

MAGDALENA RUDZIŃSKA-RADECKA*

Foundation of Research and Science Development, Warsaw, 01-793, Poland
*Address correspondence to: Magdalena Rudzińska-Radecka, magdalena.radecka@fundacjabirn.pl; magdda.rudzinska@gmail.com
Received: 11 October 2021; Accepted: 27 December 2021

Abstract: Clear cell renal cell carcinoma (KIRC) is the most common and aggressive malignancy subtype of renal neoplasm that arises from proximal convoluted tubules. It is characterized by poor clinical outcomes and high mortality of patients due to the lack of specific biomarkers for varying stages of the disease and no effective treatment. Proteases are associated with the development of several malignant tumors in humans by their ability to degrade extracellular matrices, facilitating metastasis. Herein, differentially expressed genes in KIRC cases compared to healthy kidneys were screened out from the Gene Expression Profiling Interactive Analysis (GEPIA) database. This data was applied to determine the most elevated protease in KIRC and as a result, A Disintegrin and Metalloproteinase Domain-Like Protein Decysin-1 (ADAMDEC1) was selected. This expression pattern was exclusive for KIRC and not observed for papillary and chromophobe renal cell carcinomas, in which ADAMDEC1 was at the same level in tumors and non-cancer specimens. Furthermore, the ADAMDEC1 significant increase was detected in the fourteen other human malignancies compared to healthy samples, which suggested its strong involvement in cancer development. Next, GEPIA and Pathology Atlas correlated ADAMDEC1 high expression with more advanced tumor grade and shorter survival of KIRC patients. Xena Functional Genomics Explorer presented that ADAMDEC1 could be hypermethylated in some tumor cases and one somatic mutation in the gene sequence was detected. Finally, a Search Tool for the Retrieval of Interacting Genes/Proteins; STRING base was utilized to predict the interactions of ADAMDEC1 with other molecules and construct the signaling network. In summary, ADAMDEC1 showed the tremendous potential to be the predictive marker for the KIRC and its development. Therefore, this review with data analysis can be a good base for further in vitro and in vivo research that experimentally can confirm the ADAMDEC1 as prognostic biomarkers and therapeutic target of KIRC.

Keywords: A Disintegrin and Metalloproteinase Domain-Like Protein Decysin-1; ADAMDEC1; Clear cell renal cell carcinoma; Patient specimens; Databases

Abbreviations

Introduction

Proteases, enzymes that break the protein peptide bonds, are located in the cytoplasm, mitochondria, lysosome, and extracellular matrix (King et al., 2014). By their action mechanism, they can be classified into six groups: serine, cysteine, threonine, aspartate, glutamic acid proteases and metalloproteases (Eatemadi et al., 2017).

They are responsible for non-specific degradative functions but also can catalyze specific proteolytic processing. Thus proteases are relevant in the control of multiple biological processes in living organisms. Hence, proteases contribute to (i) the fate, localization, and activity of many proteins, (ii) protein-protein interactions, (iii) generation of new bioactive molecules, (iv) cellular information processing (v) molecular signal pathways (López-Otín and Bond, 2008). As a result, proteases influence transcription, cell proliferation and differentiation, tissue morphogenesis/remodeling, neurogenesis, proteostasis, angiogenesis, ovulation, fertilization, wound healing, inflammation, autophagy, apoptosis and others (Cabral-Pacheco et al., 2020; Rudzińska et al., 2020; Rudzińska et al., 2021; Bond, 2019). Consistent with these essential roles in cell behavior and survival, alterations in proteolytic balance underlie multiple pathological conditions such as neurodegenerative disorders, inflammatory and cardiovascular diseases as well as cancer (Martinelli and Rugarli, 2010; Bühling et al., 2006; Klingler and Hardt, 2012; Yang et al., 2009). Many human tumors show the up-regulated protease level at an early stage of development, so proteases are thought to be crucial in tumor angiogenesis, invasion, and infiltration (Koblinski et al., 2000).

A Disintegrin and Metalloprotease-Like Decysin 1 (ADAMDEC1) is a unique, highly conserved secreted metalloprotease with a very rare zinc-binding motif (HEXXHXXGXXD) within the metalloprotease domain (Mueller et al., 1997). It contains pro- and catalytic domains but has a truncated disintegrin domain and lacks the typical transmembrane domain and the cytoplasmic tail (Mueller et al., 1997). The recent studies point to several vital roles of ADAMDEC1 in human diseases, including atherosclerosis (Verdugo et al., 2013), pulmonary sarcoidosis (Crouser et al., 2009), osteoarthritis (Galligan et al., 2007), Crohn’s disease (Smith et al., 2015), gastrointestinal malignancies (Oh et al., 2017) and glioblastoma (Vauléon et al., 2012). However, the natural role of ADAMDEC1 remains unknown, although its characteristic features suggest that it plays a fundamental role in the physiology of mammals.

Renal cell carcinoma (RCC) accounts for about 3% of adult malignancies, with an accelerating trend (increasing 2%–3% per decade) (Padala et al., 2020). Clear cell renal cell carcinoma (KIRC) is the most common subtype characterized by diagnostic in advanced stage and high drug-resistant nature (Makhov et al., 2018). This type of cancer does not show up early warning signs, so understanding the molecular pathogenesis of RCC and refining the diagnostic and treatment options are highly desirable (Motzer et al., 1997). The cancer propagation is attributed to proteolytic degradation of the extracellular matrix initiated by all classes of proteases (Skrzydlewska et al., 2005). It was already presented that various proteases can contribute to RCC malignancies (Drendel et al., 2017; Rudzińska et al., 2020).

Here, all up-regulated proteases were screened in 623 cases of KIRC based on the Gene Expression Profiling Interactive Analysis (GEPIA) database. Analysis elucidated ADAMDEC1 as the strongest accelerated protease in KIRC compared to healthy kidney tissues. This expression pattern was exclusive for KIRC subtype and distinguished it of chromophobe (KICH) and papillary renal (KIRP) cell carcinomas, which presented a similar ADAMDEC1 level in normal and tumor specimens. Furthermore, according to Pathology Atlas and GEPIA, the ADAMDEC1 expression level increased gradually as the tumor stage and correlated with the shorter survival of KIRC patients. Interestingly, ADAMDEC1 was elevated in other 14 types of tumors, including glioblastoma, breast, cervical, lung and others.

Additionally, Xena Functional Genomics Explorer showed that ADAMDEC1 expression could be modified by hypermethylation and one somatic mutation was detected in the gene sequence. Finally, a Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) base was used to predict the interactions of ADAMDEC1 with other proteins and visualize the signaling network.

This work points out the perspective role of ADAMDEC1 in KIRC and its particular implementation of patient outcomes. In addition, ADAMDEC1 showed the potential for experimental analysis to confirm its diagnostic value and molecular target for KIRC.

Characterization of ADAMDEC1 structure and function

ADAMDEC1 mRNA sequence was discovered in dendritic cells of human tonsils and categorized into the ADAM family (Mueller et al., 1997). However, the protein is short and consists 470 aa. Its structure is untypical since ADAMDEC1 is missing the cytoplasmic tail, transmembrane and cysteine-rich domains, thus consisting of only a signaling peptide, prodomain, metalloprotease and disintegrin-like domains (Fig. 1) (Mueller et al., 1997). Additionally, analyzing the canonical structure of ADAM enzymes-histidine (H) repeat motif (HEXXHXXGXXH) within the zinc-binding site of the metalloprotease domain, in ADAMDEC1, the third histidine is replaced by aspartate (D) serving in the motif sequence of HEXXHXXGXXD (Mueller et al., 1997; Lund et al., 2013). This rare structure was connected with the reduced catalytic activity, diminished substrate specificity and enhanced resistance to the regulation by endogenous inhibitors of metalloproteinases (TIMPs) (Lund et al., 2013; Lund et al., 2015).

Figure 1: Schematic structure of ADAMDEC1. The unique structure ADAMDEC1 consists of a signaling peptide, prodomain, metalloprotease and disintegrin-like domains.

Physiologically, the highest expression of ADAMDEC1 was found in the gastrointestinal, expressed predominantly in the macrophages in the lamina propria. The lower intensity of ADAMDEC1 was detected in seen, lymph nodes, tonsils, spleen, urinary bladder and placenta (Mueller et al., 1997; Kumagai et al., 2021).

The whole-genome transcriptomic analysis associated ADAMDEC1 expression deviations with Crohn’s disease (de Bruyn et al., 2014), rheumatoid arthritis (Galligan et al., 2007) and atherosclerosis (Papaspyridonos et al., 2006).

ADAMDEC1 downregulation in the inflamed ileal biopsy samples collected from patients with Crohn’s disease compared to healthy controls was detected. This correlation was independent of the degree of mucosal inflammation, thus underlying ADADEC1 role in the pathogenesis of this disease (de Bruyn et al., 2014).

About 24 times enhanced ADAMDEC1 level was detected in fibroblast-like synovial cells isolated from joint tissue taken from patients with rheumatoid arthritis (Galligan et al., 2007). This observation was confirmed by transcriptional datasets generated from the synovial tissue rheumatoid arthritis and osteoarthritis (Li et al., 2019).

Next, the enhanced presence of ADAMDEC1 in the circulating monocytes was found to have the third significant positive association with the occurrence of atherosclerotic plaques in the common carotid artery. Thus, it can suggest the potential function of ADAMDEC1 in atherosclerosis and plaque instability development (Verdugo et al., 2013).

In the context of cancerogenesis, the drastic reduction of ADADEC1 in the colorectal adenoma compared to the normal colorectal tissues was associated with poor prognosis of patients. However, it was not connected with the loss of macrophages in tissues, verified by CD68 expression (Kalmár et al., 2019).

Analysis of patients with craniopharyngioma revealed accelerated ADAMDEC1 expression in collected specimens compared to normal brain tissues. The primary craniopharyngioma-derived cells treated with an anti-estrogen drug, tamoxifen, reduced tumor cell proliferation and also the expression of ADAMDEC1 mRNA/protein levels (Xu et al., 2012).

A potential mechanistic role for ADAMDEC1 in developing oral squamous cell carcinoma (OSCC) has been shown in experiments based on EGF-receptor (EGFR) activation. An ectopic expression of EGFR is detected in approximately 90% of OSCC, and it is a factor that can drive an aggressive phenotype and contribute to decreased response to radiotherapy. The thrombin-stimulated platelets were shown to secrete ADAMDEC1, resulting in the cleavage of the platelet membrane-bound pro-EGF. The soluble EGF then resulted in a migratory and invasive phenotypic shift in the OSCC via EGFR signaling (Chen et al., 2018).

Finally, Jimenez-Pascual et al. (2019), Vauléon et al. (2012) identified ADAMDEC1 as one of the top up-regulated ADAM genes in glioblastoma multiforme (GBM), which correlates with poor patients prognoses. The study showed that ADAMDEC1 rapidly solubilizes the membrane-bound FGF2 and stimulates FGFR1 expressed on the GBM stem cells. Further, the transcription factor Zinc finger E-box-binding homeobox 1 (ZEB1) was activated and consequently induced ADAMDEC1 expression by inhibiting miR-203 (functioning as the tumor-suppressor) and creating a positive feedback loop. Targeting this signaling pathway decreased self-renewal and tumor growth, which highlighted the potential role of ADAMDEC1 in cancer propagation (Jimenez-Pascual et al., 2019).

Renal cell carcinoma

RCC is originated from the renal epithelium and accounts for >90% of all kidney malignancies. The disease encompasses >10 histological and molecular subtypes, of which KIRC is most common (70%–80%; accounts for most cancer-related deaths), next, papillary (10%–15%), chromophobe (5%–10%) and collecting duct (1%) RCC (Hsieh et al., 2017; Harbin et al., 2015; Muglia and Prando, 2015). These three subtypes are the most frequent histological groups that together represent more than 90% of all RCC (Muglia and Prando, 2015).

Approximately 2%–3% of all RCC are hereditary and several autosomal dominant syndromes (e.g., Hippel–Lindau disease) are described. Well-validated targets, including VHL, VEGFR and mTOR and pathways such as HGF/c-MET and Wnt/β-catenin, are strongly related to RCC pathogenesis (Maher, 2018).

RCC is a disease with complex etiologies, which the combined effect of multiple genes may cause. In this context, VHL (von Hippel-Lindau), p53, p16, p21 and p27 were shown as the primary tumor suppressor genes in RCC, in which VHL and p53 were certified to result in the development of RCC (Liu et al., 2015). Dysfunction of VHL will lead to constitutively aberrant activation of the hypoxic response, like upregulation of vascular endothelial growth factor (VEGF), which is considered to play essential roles in tumor development and angiogenesis (Zhang and Zhang, 2018). p53 is shown to suppress tumor growth and induce cell apoptosis in KIRC (Zhao et al., 2016).

Recently the genetic profiling of human RCC has increasingly been used to identify the potential of epigenetic regulatory mechanisms. It was shown that single-stranded molecules (lncRNA, miRNA) (Trilla-Fuertes et al., 2020; Rudzinska et al., 2021) and DNA methylation (Lasseigne and Brooks, 2018) are deeply involved in renal pathogenesis.

DNA methylation involves the covalent transfer of a methyl group to the C-5 position of the cytosine DNA ring by DNA methyltransferases and causes gene silencing. The genome-wide DNA methylation study in RCC identified increased global methylation in more aggressive cancer types and a potential risk factor associated with malignant transformation. The variations in sequence hypermethylation were dependents on medical aspects, such as aggressive cancer, tumor size, etc.

The potential therapeutic targets, such as small-molecule multikinase inhibitors that target VEGF receptors (sunitinib and sorafenib), the anti-VEGF antibody bevacizumab, and a mammalian target of rapamycin inhibitor temsirolimus, still are not efficient enough in the KIRC currency and chemioresistance is often observed (Larkin et al., 2009; Rudzińska et al., 2021; Penticuff and Kyprianou, 2015). It is connected with the fact that molecular mechanisms underlying this disease have not been fully elucidated. In addition, more than 50% of KIRC are detected incidentally because of a lack of effective diagnostics in the early stages of the disease; it leads to metastasis and increases the mortality rate (Czarnecka et al., 2014). This scenario emphasizes the need to discover new biomarkers for RCC subtypes with the tumor grade-dependent manner. Furthermore, exploring molecular biology of KIRC can deepen the knowledge about its mechanism and identify the novel molecular targets.

ADAMDEC1 Expression Level in Renal Cell Carcinomas and Other Human Tumors

Material and methods

ADAMDEC1 expression level analysis

The online database Gene Expression Profiling Interactive Analysis (GEPIA; http://gepia.cancer-pku.cn/index.html; Tang et al., 2017) is a communal network server that includes a transcriptome sequencing dataset of 9736 tumors together with 8587 adjacent normal tissues from TCGA and Genotype-Tissue Expression (GTEx) data sets. Here, the elevated expression of proteases was analyzed in 623 KIRC cases, compared to 100 non-tumoral kidney samples. Correlation analysis was performed on the selected gene pairwise, and ANOVA test was used to calculate the statistical significance (p < 0.05 was considered as significant).

ADAMDEC1 expression level was analyzed in 16 other human tumors at the same way.

Protease list was downloaded from MEROPS The Peptidase Database (https://www.ebi.ac.uk/merops/; Rawlings et al., 2021).

Correlation with tumor stages

ADAMDEC1 expression level with tumor grade information was collected from Pathology Atlas (https://www.proteinatlas.org/humanproteome/pathology; Uhlen et al., 2017). Data were analyzed using a non-parametric U-Mann–Whitney test (GraphPad, Prism 6.00 for Windows, Graf Pad software, San Diego, CA, USA). Data are reported as violin plots and p-value < 0.05 was considered statistically significant.

Correlation with survival of patients

By GEPIA the ADAMDEC1 ectopic expression was correlated with survival of KIRC patients. The Kaplan-Meier plots log-rank allowed for identification of predictive value of analyzed gene.

Hypermethylation and somatic mutation analysis

By USC Xena Browser (https://xenabrowser.net/; Goldman et al., 2020), the hypermethylation status of ADAMDEC1 and somatic mutations were analyzed. DNA methylation profile was detected using Infinium Human Methylation 450 and Human Methylation 27 techniques, whereas the somatic mutations were detected using the MuTect2 tool.

Protein-protein (PPI) network

The potential interactions of ADAMDEC1 with other molecules was visualized by The Search Tool for the Retrieval of Interacting Genes - STRING database (https://string-db.org/; Szklarczyk et al., 2021). PPI map was constructed, with the active interaction sources, such as experiments, databases, textminig, co-expression, gene fusion, and co-occurrence.

Results

ADAMDEC1 is the strongest overexpressed protease in KIRC and up-regulated transcript in 14 other human tumors compared to non-tumor controls

Using the GEPIA database, all significantly up-regulated genes in 623 KIRC cases compared to the healthy cohort (100 samples) were listed. Next, by MEROPS The Peptidase Database, all changed proteases were selected and a shortlist of twentieth strongest changed proteases is included below (Table 1).

The strongest changed ADADEC1 was taken for further analysis and firstly, the significant (p < 0.05) evaluation of ADAMDEC1 expression in KIRC samples compared to normal kidney tissues was confirmed (Fig. 2A). Following, the higher ADAMDEC1 expression was detected in 14 other human tumors (Fig. 2B), including: invasive breast carcinoma (BRCA), cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC), lymphoid neoplasm diffuse large B-cell lymphoma (DLBC), esophageal carcinoma (ESCA), glioblastoma multiforme (GBM), head and neck squamous cell carcinoma (HNSC), lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC), ovarian serous cystadenocarcinoma (OV), pancreatic adenocarcinoma (PAAD), cutaneous skin melanoma (SKCM), stomach adenocarcinoma (STAD), testicular germ cell tumors (TGCT) and uterine corpus endometrial carcinoma (UCEC). Interestingly, ADAMDEC1 was not elevated in kidney chromophobe (KICH) and kidney renal papillary cell carcinoma (KIRP) samples to serve as distinguishing features of the KIRC subtype.

Figure 2: (A) ADAMDEC1 expression level in clear cell renal cell carcinoma samples (KIRC; tumor-red) compared to control kidney tissues (non-tumor-grey); ANOVA test; * p < 0.05. (B) ADAMDEC1 expression level in other human tumors: (BRCA) beast invasive carcinoma, (CESC) cervical squamous cell carcinoma and endocervical adenocarcinoma, (DLBC) lymphoid neoplasm diffuse large B-cell lymphoma, (ESCA) esophageal carcinoma, (GBM) glioblastoma multiforme, (HNSC) head and neck squamous cell carcinoma, (LUAD) lung adenocarcinoma, (LUSC) lung squamous cell carcinoma, (OV) ovarian serous cystadenocarcinoma, (PAAD) pancreatic adenocarcinoma, (SKCM) skin cutaneous melanoma, (STAD) stomach adenocarcinoma, (TGCT) testicular germ cell tumors and (UCEC) uterine corpus endometrial carcinoma–colored red. The lack of significant difference in (KICH) kidney chromophobe and (KIRP) kidney renal papillary cell carcinoma colored black. Analyzed by GEPIA database.

ADAMDEC1 enhanced presence correlates with high-grade KIRC tumor and worse prognosis of patients

Further, ADADEC1 expression was correlated with clinical outcomes of KIRC cases. As a result, the accelerating ADAMDEC1 expression level tendency corresponded to higher tumor stages; from I to IV stages (Fig. 3A). Significant different (*p < 0.05) was noticed between I and III and IV grades. In addition, its enhanced expression notably (p = 0.046) corresponded to a worse survival prognosis of patients (Fig. 3B).

Figure 3: The ADAMDEC1 expression pattern in the clear renal cancer cell (KIRC) and its relation to (A) tumor stages of KIRC patients: I, II, III, and IV staging. According to the American Joint Committee on Cancer Staging (AJCC), Stage I (n = 259) is characterized by the tumor only in the kidney with a max. 7 cm size. Stage II (n = 31): the tumor is found only in the kidney and is larger than 7 cm. Stage III (n = 112): the tumor has grown into major veins within the kidney or perinephric tissue, which is the connective, fatty tissue around the kidneys. However, it has not grown into the adrenal gland on the same side of the body as the tumor. Stage IV (n = 88): the tumor has spread to areas beyond Gerota’s fascia and extends into the adrenal gland on the same body side as the tumor. *p < 0.05 U-Mann-Whitney test. (B) The ADAMDEC1 mRNA expression relation to survival of patients with KIRC; GEPIA database.

ADAMDEC1 expression can be control by hypermethylation and mutations in KIRC cases

To investigate the possible relationship between hypermethylation and ADAMDEC1 expression, an analysis of methylation among KIRC cases was performed using USC Xena browser (Fig. 4C). The hypermethylation was detected in part of both tumors and non-tumor controls.

Figure 4: ADAMDEC1 hypermethylation and somatic mutations in clear cell renal cell carcinoma (KIRC) cases. (A) Hierarchical clustering of KIRC tissues–red–solid tissue normal (non-cancer tissue surrounding), blue–primary tumor, (B) age analyzed patients (C) ADAMDEC1 sequence hypermethylation (red); grey–no analyzed samples; hypermethylated sequences localization: (red line: hg38 chr8:24,383,084-24,383,084; green line: hg38 chr8:24,384,181-24,384,182; blue line: hg38 chr8:24,384,596-24,384,597; violet line: hg38 chr8:24,394,941-24,394942). (D) Somatic mutation analyzing. Blue point (a red arrow) indicates missense variant p.W23R/substitution T to C (hg38 ch8:24,384,571) in one KIRC case. Analyzed by USC Xena Browser.

Additionally, the somatic mutations were checked and missense variant p.W23R/substitution T to C (hg38 ch8:24,384,571) in one case was detected (Fig. 4D).

ADAMDEC1 potentially interact with other protein creating a signaling cluster

Proteins and protein-protein interactions form the backbone of the cellular complex machinery. Thus, for the full understanding of biological phenomena, the full network needs to be considered. The Search Tool for the Retrieval of Interacting Genes-STRING database collected scored and integrated all publicly available sources of protein-protein interaction information, and to complemented these with computational predictions. The STRING database uses protein-protein co-occurrence as one type of text mining evidence for protein interactions. Text mining (text data mining) transforms unstructured text into a structured format to identify meaningful patterns and new insights. Text mining systems are designed to extract information from the text in a domain-oriented manner (Franceschini et al., 2013).

Here, STRING was used to seek potential interactions between ADAMDEC1 with other molecules (Fig. 5A). As a result, PPI network was constructed, with the active interaction sources, including experiments, databases, textminig, co-expression, gene fusion, and co-occurrence (Fig. 5B). In the matrix, the nodes correspond to the proteins, and the edges represent the interactions. Using an interaction score threshold of 0.7 (high confidence), the string PPI analysis yielded a highly clustered network (clustering coefficient: 0.73) containing 11 nodes with 19 edges (enrichment p-value < 0.008), false discovery rate = 1.42E-0.5.

Figure 5: STRING Protein-protein interaction (PPI) analyses. (A) Predicted functional partners of ADAMDEC1 (B) PPI network connectivity of human ADAMDEC1 interactions cluster. The network contains 11 nodes with 19 edges. The confidence score threshold was set at 0.8 (high) for analyses. The gene interaction data used to build the network was based on direct physical interactions that are either experimentally derived or computationally predicted.

As the outcome, the co-expression was found for CXCL9 and SLAMF9. Text mining systems showed the interactions with proteins, including DUOXA1, ZNF511, CXCL9, KNDC1, GPR123, DUOX2A, SMPDL3A, SLAMF8, ERIN1, TACSTD2, with a different score value.

Conclusions and further perspectives

Recently, ADAMDEC1 has been connected with the development of varied human diseases, including cancers. Herein, the presented review with database analysis demonstrates the first evidence for positive relations of ADAMDEC1 expression with KIRC cases and its worse outcomes. It can implicate the potential role of ADAMDEC1 in the pathogenesis of KIRC. The presented findings support the further experimental analysis of ADAMDEC1 as a KIRC biomarker and suggest its essential function in diagnosing and promoting KIRC. If validated, this biomarker can significantly facilitate the introduction of new therapies based on protease inhibition that profoundly affects the efficient treatment of KIRC patients. Additionally, this knowledge can be helpful in research based on specific protease inhibitors against KIRC development.

Author Contribution: The author confirms sole responsibility for the following: study conception and design, data collection, analysis and interpretation of results, and manuscript preparation.

Funding Statement: The author received no specific funding for this study.

Conflicts of Interest: The author declares that they have no conflicts of interest to report regarding the present study.

This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.