Autophagy-Related MicroRNA: Tumor miR-125b and Thyroid Cancers

Background: Autophagy is a stress response mechanism that causes cellular components to degrade. Its defects were associated with multiple pathologies, including cancers. Thyroid cancer is known to be the most prevalent form of malignant neoplasm among endocrine tumors. The aim of the study was to seek and comprehensively explore the role of autophagy related genes and proteins play in thyroid cancers through bioinformatics analysis with their detection in the tissue samples. Methods: Bioinformatics analysis was performed to investigate autophagy related proteins and genes involvement in thyroid cancer progression. The experimental verification was done in cancer samples of one hundred and three patients with thyroid pathology included in the study. The miR-125blevel was detected by PCR in real time. Results and discussion: The bioinformatics analysis verified the miR-125b as a regulatory mechanism in autophagy. Its expression in patients with PTC was reduced by 6.75 times in cancer patients compared to the patients with benign tumors. The BRAFV600E mutations were associated with a decrease in hsa-miR-125b expression by 12.67 times compared to tumors with the wild-type gene. Conclusions: Our findings revealed involvement of the autophagy related proteins in cancer progression. The significant mechanisms of regulation are non-coding RNA sequences implicated in a variety of oncogenic processes. We found that miR-125b is a potential maker in thyroid cancer invasion, BRAV600E mutational status and risk of recurrence.


Introduction
Autophagy is a stress response mechanism that leads to the degradation of cellular components, the defects of which are associated with multiple conditions [1]. Mounting evidence reveals that autophagy is a factor in thyroid cancer. Modulating autophagy has become a prominent target for improving the monitoring of cancer patients [2]. The principal triggers that determine the induction of autophagy are starvation, reactive oxygen species, endoplasmic stress and hypoxia. There are multiple channels that modulate autophagy, such as AMPK and PI3K orchestrating the protein cascade (ULK1, ATG16L1, MAP1LC3B, PRKAA1, ATG12, ATG14, ATG5, mTOR and BECN1) [2].
Papillary thyroid cancer (PTC) is the most frequent differentiated tumor. The malignant neoplasms rate in the thyroid gland is increasing annually worldwide [3]. The molecular mechanisms underlying the progress of PTC are complex and diversified. A number of key mutations are associated with PTC progression, including tyrosine kinase RET1 translocation, RAF protein mutation (BRAFV600E) and telomerase mutation (TERT1) [4]. However, their role in PTC is relatively low compared to the multiplication of oncogenic signals such as epigenetic factors.
We aimed to seek and comprehensively explore the role of autophagy related genes and proteins play in thyroid cancers through bioinformatics analysis with their detection in the tissue samples.

Protein-Protein Interaction (PPI) Analysis
The Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) used to review all known and predicted associations between proteins. The interacting network of autophagy-related proteins with other associated proteins were visualized through the STRING database. The cutoff value of interaction score was set to 0.15. A Venn diagram was utilized to display the interactions.

MicroRNA and Genes Associations
For pathway enrichment analysis, both Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and Gene Ontology (GO) enrichment analysis was performed to visualize the cellular pathways. The cutoff Q value and P value were set at 0.05 and 0.01, respectively. The DIANA miRpath and PANTHER databases were used for bioinformatics analysis. HumanTargetScan (ver. 8.0) and miRTarBase (ver. 9.0) were used to seek target genes for miR-125. To identify the transcription factors, Transmir base (version 2.0) was used (https://www.cuilab.cn/transmir; accessed on 27 January 2023).

Clinical and Experimental Tools
The study included one hundred and three patients with thyroid pathology treated at the Cancer Research Institute of the Tomsk National Research Medical Center: sixty-seven patients had PTC and thirty-six patients had benign thyroid pathology. The T1-2N0M0 stage detected in twenty-seven patients and T3-4N0-1M0 was detected in forty patients. Metastases in regional lymph nodes were presented in twenty-seven patients. According to the histological subtype of PTC, patients were divided into a follicular subtype group (12 patients) and a classical subtype group (55 patients). The BRAFV600E mutation was identified in 18 patients. According to the criteria of the American Thyroid Association (ATA, American Thyroid Association), patients were divided into recurrence risk groups as follows: low risk of thyroid cancer recurrence was observed in 23 patients, intermediate risk in 25 patients and high risk in 19 patients. This work was approved by the Local Ethical Committee in the Cancer Research Institute of the Tomsk National Research Medical Center. All procedures were carried out in accordance with the Protocol of the Helsinki Declaration on Human Rights (1964). All patients included in the study signed an informed consent to participate in the study.
The material of the study was samples of tumor and unchanged tissue obtained during surgery. Samples were taken in points located at a distance of at least 1 cm from the tumor boundaries. The samples were frozen and stored at a temperature of −80 • C.

DNA Extraction
DNA was isolated using the FFPET DNA Extraction Kit (Biolink, Russia). To assess the amount of isolated DNA, its concentration was assessed on a NanoDrop-2000 spectrophotometer (Thermo Scientific, Waltham, MA, USA). The resulting DNA was used for real-time PCR.

BRAFV600E Mutation Analysis
The BRAFV600E mutation was determined using the Real-time-PCR-BRAF-V600E kit (Biolink, Novosibirsk, Russia), designed to detect the GTG→GGG point mutation in codon 600 of the BRAF gene. The analysis was carried out by real-time allele-specific PCR.

Quantitative PCR with Reverse Transcription in Real Time
The tumor samples were incubated in RNAlater solution (Ambion, Foster City, CA, USA) for 24 h at 4 • C and then stored at −80 • C. MicroRNA was extracted using an isolation Kit (Biolabmix, Novosibirsk, Russia). The quality of the isolated nucleic acids was carried out using capillary electrophoresis on a TapeStation device (Agilient Technologies, Santa Clara, CA, USA). The RIN ranged from 2.2-3.3.
Statistical analysis was performed using the Statistica 12.0 software package (Tulsa, OK, USA). Normal distribution check was conducted using the Kolmogorov-Smirnov test. The results of gene expression determination are presented as Me (Q1; Q3). The Mann-Whitney test was used to assess significant differences. Differences were considered significant when p < 0.05.

Investigation of MicroRNAs Regulating Proteins Associated with Autophagy
Autophagy is a universal process of cellular viability as well as cancer cells. The key autophagy associated proteins are used to study the STRING database: ULK1, ATG16L1, MAP1LC3B, PRKAA1, ATG12, ATG14, ATG5, mTOR and BECN1 ( Figure 1). The analysis showed the following associations between autophagy-associated proteins with multiple links. The edges colored purple reflect experimentally proven relationships, bluerelationships, information about which is obtained from supervised databases. Green edges reflect the neighborhood in the genome, blue-the joint occurrence, red-the fusion of genes. Light green ribs mean a joint mention of these proteins in Pub-Med Abstracts, black-co-expression, and light blue-homology. The found data highlighted the multiple interactions between the proteins.
The bioinformatics analysis in the DIANA microRNA database identified microRNAs regulating proteins associated with autophagy (Table S1). Twenty-two microRNAs were found for ATG12. ATG14 was regulated by the twenty-three microRNAs, ATG16L1 by seventy-five microRNAs, ATG5 by twenty-two microRNAs, BECN1 by forty microRNAs, MAP1LC3 by thirty microRNAs, ULK1 by one hundred six microRNAs, mTOR by thirtythree microRNAs, and PRKAA1 by two hundred sixteen microRNAs.

MicroRNAs Associated with the Thyroid Cancers
At the next step, data on the microRNAs associated with the thyroid cancer and the proteins associated with autophagy were analyzed. Twenty-three microRNAs were identified with upregulated-expression and thirty-nine had oncosuppressor activity, which regulate all stages of oncogenesis (Table 1). However, among the targets associated with autophagy, only one microRNA is shown. Thus, miR-125b are able to regulate transcription and translation of the ATG5 protein, mediating the resistance to therapy in papillary thyroid cancer [19].
Genes 2023, 14, x FOR PEER REVIEW Figure 1. Protein-protein interactions between proteins associated with autophagy (https: db.org/; accessed on 27 January 2023). The edges colored purple reflect experimentally prov tionships, blue-relationships, information about which is obtained from supervised da Green edges reflect the neighborhood in the genome, blue-the joint occurrence, red-the f genes. Light green ribs mean a joint mention of these proteins in Pub-Med Abstracts, bla expression, and light blue-homology. Note: Interactions: -from verified databases -experimentally determined interactions. Predicted interactions: -neighborhood in nome; -gene fusion; -occurrence. Other: -joint mention of these pro PubMed Abstract; -co-expression; -homology.
The bioinformatics analysis in the DIANA microRNA database identified micr regulating proteins associated with autophagy (Table S1). Twenty-two microRNA found for ATG12. ATG14 was regulated by the twenty-three microRNAs, ATG16L1 enty-five microRNAs, ATG5 by twenty-two microRNAs, BECN1 by forty micro MAP1LC3 by thirty microRNAs, ULK1 by one hundred six microRNAs, mTOR by three microRNAs, and PRKAA1 by two hundred sixteen microRNAs.

MicroRNAs Associated with the Thyroid Cancers
At the next step, data on the microRNAs associated with the thyroid cancer a proteins associated with autophagy were analyzed. Twenty-three microRNAs wer tified with upregulated-expression and thirty-nine had oncosuppressor activity, regulate all stages of oncogenesis (Table 1). However, among the targets associate autophagy, only one microRNA is shown. Thus, miR-125b are able to regulate tra tion and translation of the ATG5 protein, mediating the resistance to therapy in pa thyroid cancer [19].  The bioinformatics analysis in the DIANA microRNA database identified microRNAs regulating proteins associated with autophagy (Table S1). Twenty-two microRNAs were found for ATG12. ATG14 was regulated by the twenty-three microRNAs, ATG16L1 by seventy-five microRNAs, ATG5 by twenty-two microRNAs, BECN1 by forty microRNAs, MAP1LC3 by thirty microRNAs, ULK1 by one hundred six microRNAs, mTOR by thirtythree microRNAs, and PRKAA1 by two hundred sixteen microRNAs.

MicroRNAs Associated with the Thyroid Cancers
At the next step, data on the microRNAs associated with the thyroid cancer and the proteins associated with autophagy were analyzed. Twenty-three microRNAs were identified with upregulated-expression and thirty-nine had oncosuppressor activity, which regulate all stages of oncogenesis (Table 1). However, among the targets associated with autophagy, only one microRNA is shown. Thus, miR-125b are able to regulate transcription and translation of the ATG5 protein, mediating the resistance to therapy in papillary thyroid cancer [19]. -neighborhood in the geher: -joint mention of these proteins in mology. microRNA database identified microRNAs (Table S1). Twenty-two microRNAs were wenty-three microRNAs, ATG16L1 by sev-icroRNAs, BECN1 by forty microRNAs, hundred six microRNAs, mTOR by thirtyd sixteen microRNAs.
Cancers associated with the thyroid cancer and the zed. Twenty-three microRNAs were iden--nine had oncosuppressor activity, which wever, among the targets associated with s, miR-125b are able to regulate transcripating the resistance to therapy in papillary cancer.

MicroRNAs Associated with the Thyroid Cancers
At the next step, data on the microRNAs associated with the thyroid cancer and the proteins associated with autophagy were analyzed. Twenty-three microRNAs were identified with upregulated-expression and thirty-nine had oncosuppressor activity, which regulate all stages of oncogenesis (Table 1). However, among the targets associated with autophagy, only one microRNA is shown. Thus, miR-125b are able to regulate transcription and translation of the ATG5 protein, mediating the resistance to therapy in papillary thyroid cancer [19].

MicroRNAs Associated with the Thyroid Cancers
At the next step, data on the microRNAs associated with the thyroid cancer and the proteins associated with autophagy were analyzed. Twenty-three microRNAs were identified with upregulated-expression and thirty-nine had oncosuppressor activity, which regulate all stages of oncogenesis (Table 1). However, among the targets associated with autophagy, only one microRNA is shown. Thus, miR-125b are able to regulate transcription and translation of the ATG5 protein, mediating the resistance to therapy in papillary thyroid cancer [19].  The bioinformatics analysis in the DIANA microRNA database identified microRNAs regulating proteins associated with autophagy (Table S1). Twenty-two microRNAs were found for ATG12. ATG14 was regulated by the twenty-three microRNAs, ATG16L1 by seventy-five microRNAs, ATG5 by twenty-two microRNAs, BECN1 by forty microRNAs, MAP1LC3 by thirty microRNAs, ULK1 by one hundred six microRNAs, mTOR by thirtythree microRNAs, and PRKAA1 by two hundred sixteen microRNAs.

MicroRNAs Associated with the Thyroid Cancers
At the next step, data on the microRNAs associated with the thyroid cancer and the proteins associated with autophagy were analyzed. Twenty-three microRNAs were identified with upregulated-expression and thirty-nine had oncosuppressor activity, which regulate all stages of oncogenesis (Table 1). However, among the targets associated with autophagy, only one microRNA is shown. Thus, miR-125b are able to regulate transcription and translation of the ATG5 protein, mediating the resistance to therapy in papillary thyroid cancer [19].  The bioinformatics analysis in the DIANA microRNA database identified microRNAs regulating proteins associated with autophagy (Table S1). Twenty-two microRNAs were found for ATG12. ATG14 was regulated by the twenty-three microRNAs, ATG16L1 by seventy-five microRNAs, ATG5 by twenty-two microRNAs, BECN1 by forty microRNAs, MAP1LC3 by thirty microRNAs, ULK1 by one hundred six microRNAs, mTOR by thirtythree microRNAs, and PRKAA1 by two hundred sixteen microRNAs.

MicroRNAs Associated with the Thyroid Cancers
At the next step, data on the microRNAs associated with the thyroid cancer and the proteins associated with autophagy were analyzed. Twenty-three microRNAs were identified with upregulated-expression and thirty-nine had oncosuppressor activity, which regulate all stages of oncogenesis (Table 1). However, among the targets associated with autophagy, only one microRNA is shown. Thus, miR-125b are able to regulate transcription and translation of the ATG5 protein, mediating the resistance to therapy in papillary thyroid cancer [19].  The bioinformatics analysis in the DIANA microRNA database identified microRNAs regulating proteins associated with autophagy (Table S1). Twenty-two microRNAs were found for ATG12. ATG14 was regulated by the twenty-three microRNAs, ATG16L1 by seventy-five microRNAs, ATG5 by twenty-two microRNAs, BECN1 by forty microRNAs, MAP1LC3 by thirty microRNAs, ULK1 by one hundred six microRNAs, mTOR by thirtythree microRNAs, and PRKAA1 by two hundred sixteen microRNAs.

MicroRNAs Associated with the Thyroid Cancers
At the next step, data on the microRNAs associated with the thyroid cancer and the proteins associated with autophagy were analyzed. Twenty-three microRNAs were identified with upregulated-expression and thirty-nine had oncosuppressor activity, which regulate all stages of oncogenesis (Table 1). However, among the targets associated with autophagy, only one microRNA is shown. Thus, miR-125b are able to regulate transcription and translation of the ATG5 protein, mediating the resistance to therapy in papillary thyroid cancer [19].  The bioinformatics analysis in the DIANA microRNA database identified microRNAs regulating proteins associated with autophagy (Table S1). Twenty-two microRNAs were found for ATG12. ATG14 was regulated by the twenty-three microRNAs, ATG16L1 by seventy-five microRNAs, ATG5 by twenty-two microRNAs, BECN1 by forty microRNAs, MAP1LC3 by thirty microRNAs, ULK1 by one hundred six microRNAs, mTOR by thirtythree microRNAs, and PRKAA1 by two hundred sixteen microRNAs.

MicroRNAs Associated with the Thyroid Cancers
At the next step, data on the microRNAs associated with the thyroid cancer and the proteins associated with autophagy were analyzed. Twenty-three microRNAs were identified with upregulated-expression and thirty-nine had oncosuppressor activity, which regulate all stages of oncogenesis (Table 1). However, among the targets associated with autophagy, only one microRNA is shown. Thus, miR-125b are able to regulate transcription and translation of the ATG5 protein, mediating the resistance to therapy in papillary thyroid cancer [19]. A bioinformatic analysis was performed to look for potential target genes and proteinprotein interactions between their products ( Table 2). Forty-four genes were found to be targets for the miR-125b. The detailed review showed the participation of the miR-125b in multiple processes: cell motility, immune response, DNA methylation, growth and transcriptional factors signal transduction. The MAPK and AKT/mTOR signaling cascades are found to be the key ones in autophagy [63]. PIK3CA, the initial component of the AKT/mTOR cascade [64]. CCDC6, protein 6 containing an antenna domain, interacts with CREB1 (reaction element binding protein cAMP 1) and suppresses its transcriptional activity [65].

MAPK3
Mitogen-activated protein kinase 3; The STRING database was used for the creation of the protein-protein interactions between the proteins of miR-125b associated genes ( Figure 2). The complex scheme recovers only four target proteins that are involved in thyroid cancer progression.

UBR4
E3 ubiquitin-protein ligase UBR4; WNT11 Wingless-type mmtv integration site family, member 11 The STRING database was used for the creation of the protein-protein interactions between the proteins of miR-125b associated genes (Figure 2). The complex scheme recovers only four target proteins that are involved in thyroid cancer progression.  Protein-protein interactions between proteins associated with autophagy (https://stringccessed on 27 January 2023). The edges colored purple reflect experimentally proven relablue-relationships, information about which is obtained from supervised databases. es reflect the neighborhood in the genome, blue-the joint occurrence, red-the fusion of ht green ribs mean a joint mention of these proteins in Pub-Med Abstracts, black-con, and light blue-homology. Note: Interactions: -from verified databases; entally determined interactions. Predicted interactions: -neighborhood in the ge--gene fusion; -occurrence. Other: -joint mention of these proteins in bstract; -co-expression; -homology. bioinformatics analysis in the DIANA microRNA database identified microRNAs g proteins associated with autophagy (Table S1). Twenty-two microRNAs were ATG12. ATG14 was regulated by the twenty-three microRNAs, ATG16L1 by sev-microRNAs, ATG5 by twenty-two microRNAs, BECN1 by forty microRNAs, 3 by thirty microRNAs, ULK1 by one hundred six microRNAs, mTOR by thirty-roRNAs, and PRKAA1 by two hundred sixteen microRNAs.
croRNAs Associated with the Thyroid Cancers he next step, data on the microRNAs associated with the thyroid cancer and the associated with autophagy were analyzed. Twenty-three microRNAs were identh upregulated-expression and thirty-nine had oncosuppressor activity, which all stages of oncogenesis (Table 1). However, among the targets associated with y, only one microRNA is shown. Thus, miR-125b are able to regulate transcriptranslation of the ATG5 protein, mediating the resistance to therapy in papillary ancer [19].
ioinformatics analysis in the DIANA microRNA database identified microRNAs proteins associated with autophagy (Table S1). Twenty-two microRNAs were ATG12. ATG14 was regulated by the twenty-three microRNAs, ATG16L1 by sev-microRNAs, ATG5 by twenty-two microRNAs, BECN1 by forty microRNAs, 3 by thirty microRNAs, ULK1 by one hundred six microRNAs, mTOR by thirty-oRNAs, and PRKAA1 by two hundred sixteen microRNAs.
roRNAs Associated with the Thyroid Cancers e next step, data on the microRNAs associated with the thyroid cancer and the ssociated with autophagy were analyzed. Twenty-three microRNAs were idenh upregulated-expression and thirty-nine had oncosuppressor activity, which ll stages of oncogenesis (Table 1). However, among the targets associated with y, only one microRNA is shown. Thus, miR-125b are able to regulate transcripranslation of the ATG5 protein, mediating the resistance to therapy in papillary ncer [19].

of 18
otein-protein interactions between proteins associated with autophagy (https://stringessed on 27 January 2023). The edges colored purple reflect experimentally proven relalue-relationships, information about which is obtained from supervised databases. reflect the neighborhood in the genome, blue-the joint occurrence, red-the fusion of t green ribs mean a joint mention of these proteins in Pub-Med Abstracts, black-coand light blue-homology. Note: Interactions: -from verified databases; ntally determined interactions. Predicted interactions: -neighborhood in the ge--gene fusion; -occurrence. Other: -joint mention of these proteins in stract; -co-expression; -homology.
ioinformatics analysis in the DIANA microRNA database identified microRNAs proteins associated with autophagy (Table S1). Twenty-two microRNAs were TG12. ATG14 was regulated by the twenty-three microRNAs, ATG16L1 by sev-icroRNAs, ATG5 by twenty-two microRNAs, BECN1 by forty microRNAs, by thirty microRNAs, ULK1 by one hundred six microRNAs, mTOR by thirty-RNAs, and PRKAA1 by two hundred sixteen microRNAs.
oRNAs Associated with the Thyroid Cancers next step, data on the microRNAs associated with the thyroid cancer and the sociated with autophagy were analyzed. Twenty-three microRNAs were idenupregulated-expression and thirty-nine had oncosuppressor activity, which l stages of oncogenesis (Table 1). However, among the targets associated with , only one microRNA is shown. Thus, miR-125b are able to regulate transcripanslation of the ATG5 protein, mediating the resistance to therapy in papillary cer [19]. t of microRNAs associated with thyroid cancer.

MicroRNAs Associated with the Thyroid Cancers
At the next step, data on the microRNAs associated proteins associated with autophagy were analyzed. Twen tified with upregulated-expression and thirty-nine had regulate all stages of oncogenesis (Table 1). However, am autophagy, only one microRNA is shown. Thus, miR-125 tion and translation of the ATG5 protein, mediating the re thyroid cancer [19].  The bioinformatics analysis in the DIANA microRNA regulating proteins associated with autophagy (Table S1). found for ATG12. ATG14 was regulated by the twenty-three enty-five microRNAs, ATG5 by twenty-two microRNAs, MAP1LC3 by thirty microRNAs, ULK1 by one hundred six three microRNAs, and PRKAA1 by two hundred sixteen m

MicroRNAs Associated with the Thyroid Cancers
At the next step, data on the microRNAs associated w proteins associated with autophagy were analyzed. Twent tified with upregulated-expression and thirty-nine had o regulate all stages of oncogenesis (Table 1). However, am autophagy, only one microRNA is shown. Thus, miR-125b tion and translation of the ATG5 protein, mediating the re thyroid cancer [19]. in-protein interactions between proteins associated with autophagy (https://stringd on 27 January 2023). The edges colored purple reflect experimentally proven rela--relationships, information about which is obtained from supervised databases. flect the neighborhood in the genome, blue-the joint occurrence, red-the fusion of een ribs mean a joint mention of these proteins in Pub-Med Abstracts, black-colight blue-homology. Note: Interactions: -from verified databases; lly determined interactions. Predicted interactions: -neighborhood in the gegene fusion; -occurrence. Other: -joint mention of these proteins in act; -co-expression; -homology.
formatics analysis in the DIANA microRNA database identified microRNAs teins associated with autophagy (Table S1). Twenty-two microRNAs were 12. ATG14 was regulated by the twenty-three microRNAs, ATG16L1 by sev-roRNAs, ATG5 by twenty-two microRNAs, BECN1 by forty microRNAs, thirty microRNAs, ULK1 by one hundred six microRNAs, mTOR by thirty-As, and PRKAA1 by two hundred sixteen microRNAs.
As Associated with the Thyroid Cancers xt step, data on the microRNAs associated with the thyroid cancer and the iated with autophagy were analyzed. Twenty-three microRNAs were idenregulated-expression and thirty-nine had oncosuppressor activity, which ages of oncogenesis (Table 1). However, among the targets associated with ly one microRNA is shown. Thus, miR-125b are able to regulate transcriplation of the ATG5 protein, mediating the resistance to therapy in papillary r [19].
microRNAs associated with thyroid cancer.
The bioinformatics analysis in the DIANA microRNA dat regulating proteins associated with autophagy (Table S1). Tw found for ATG12. ATG14 was regulated by the twenty-three mi enty-five microRNAs, ATG5 by twenty-two microRNAs, BE MAP1LC3 by thirty microRNAs, ULK1 by one hundred six m three microRNAs, and PRKAA1 by two hundred sixteen micro

MicroRNAs Associated with the Thyroid Cancers
At the next step, data on the microRNAs associated with proteins associated with autophagy were analyzed. Twenty-th tified with upregulated-expression and thirty-nine had onco regulate all stages of oncogenesis (Table 1). However, among autophagy, only one microRNA is shown. Thus, miR-125b ar tion and translation of the ATG5 protein, mediating the resist thyroid cancer [19].  The bioinformatics analysis in the DIANA microRNA data regulating proteins associated with autophagy (Table S1). Tw found for ATG12. ATG14 was regulated by the twenty-three mic enty-five microRNAs, ATG5 by twenty-two microRNAs, BE MAP1LC3 by thirty microRNAs, ULK1 by one hundred six mi three microRNAs, and PRKAA1 by two hundred sixteen microR

MicroRNAs Associated with the Thyroid Cancers
At the next step, data on the microRNAs associated with proteins associated with autophagy were analyzed. Twenty-th tified with upregulated-expression and thirty-nine had onco regulate all stages of oncogenesis (Table 1). However, among autophagy, only one microRNA is shown. Thus, miR-125b are tion and translation of the ATG5 protein, mediating the resista thyroid cancer [19]. KEGG pathway enrichment analysis revealed microRNA targets. They are components of the following signaling pathways: ECM-receptor interaction, proteoglycans in cancer, lysine degradation, adherent junction, viral carcinogenesis, thyroid cancer, focal adhesion, one carbon pool by folate, TGF-beta signaling pathway, colorectal cancers pathway (APC and Lynch syndrome), central carbon metabolism in cancer and FoxO signaling pathway. Of particular interest are the components involved in the molecular mechanisms of thyroid cancer development. Its main components are listed in Table 3.
The ECM-receptor interaction is essential in the cancer heterogeneity investigation [66]. Proteoglycans are molecules, responsible for cellular behavior modulation [67]. Cadherins and catenins are the central cell-cell adhesion molecules in adherens junctions (AJs). The significance of cadherin and catenin in cancer initiation and progression is well known [68]. The additional pathogenic mechanism is Lysine-specific demethylase 1 (LSD1) that stabilizes hypoxia-inducible factor 1α (HIF1α) to promote tumor progression [69]. The associations between viruses and cancer have been conducted in several studies and verified in meta-analyses [70]. Molecular alterations underlying PTC progression include deregulation of focal adhesion kinase (FAK), and PTC progression correlates with mRNA FAK-Del33 and pY397-FAK [71].  KEGG analysis revealed that one carbon pool by folate was connected to PTC pathogenesis [72]. A study by Xie et al. suggested that TGF-β1 triggers invasion and migration of PTC cells by inhibiting the expression of lncRNA-NEF [73].
We revealed gene and signaling components responsible for PTC manifestation. Lynch syndrome describing a familial cancer syndrome with four DNA mismatch repair genes mutations was found to be associated with PTC [74]. A high risk of familial adenomatous polyposis (FAP) was found in PTC patients [75].
PTC exhibits impaired energy metabolism, abnormal cell division, and more aggressive behavior through metabolic reprogramming and changes in the central carbon metabolism [76]. The FOXo signaling cascade is the important regulator of oncogenesis, which involves the mechanism of the RBM47/SNHG5/FOXO3 axis in PTC and also affects the triggering of autophagy [77].
Recent studies revealed the EZH2 gene (histone methyltransferase) in modulating the cancer invasiveness and metastatic potential [81]. NF-κB induction is a key process in cancer progression and promotes NFKB1 [82]. The somatic mutations were identified in AKT genes and is associated with PTC risk [83].
Thyroid cancers are prevalent in women, indicating the sex hormone disorders as an essential risk factor [84]. The peroxisome proliferator receptor gamma (PPARγ) alteration associated with steroids transcriptional factors was detected only in follicular thyroid cancer and follicular adenoma [85]. STAT3 transcriptional factor regulates PTC metastasis via glycolysis modification [78]. PRDM1 is a key target in Hashimoto thyroiditis [79]. TP53 is an oncosuppressor associated with metastatic PTC [80].

MicroRNA 125b Expression in Thyroid Pathologies
Recent studies revealed the EZH2 gene (histone methyltransferase) in modulating the cancer invasiveness and metastatic potential [81]. NF-κB induction is a key process in cancer progression and promotes NFKB1 [82]. The somatic mutations were identified in AKT genes and is associated with PTC risk [83].
Thyroid cancers are prevalent in women, indicating the sex hormone disorders as an essential risk factor [84]. The peroxisome proliferator receptor gamma (PPARγ) alteration associated with steroids transcriptional factors was detected only in follicular thyroid cancer and follicular adenoma [85].

MicroRNA 125b Expression in Thyroid Pathologies
The miRNA expression was reduced by 6.75 times in patients with PTC compared with the benign tumors. miRNAs are found to play a multifaceted role in oncogenesis. It is believed that they could be both oncoproteins and oncosuppressors, guiding the cancer progression potential [11]. Table 4 presents data on the microRNA expression in PTC. Thus, tumor size and regional lymph nodes involvement were not associated with the hsa-miR-125b level. However, in patients with the thyroid gland capsule invasion, hsa-miR-125-b expression increased 2.36 times compared to patients without signs of invasion. Note: #-significance of differences compared to patients with benign tumors, p < 0.05; *-significance of differences compared to patients with signs of tumor, p < 0.05; **-significance of differences compared to patients with follicular subtype of PTC, p < 0.05; ***-significance of differences compared to patients with BRAFV600E mutation, p < 0.05; ****-significance of differences compared with patients in group with low risk of recurrence after the treatment, p < 0.05.
The histological subtype of thyroid cancers was associated with miR-125b expression. A decrease of 37.0 times was found in the follicular subtype compared to tumors with the classical subtype. Additionally, in patients with positive BRAF V600E mutational status, a significant decrease in hsa-miR-125b levels by 12.67 times was shown compared to tumors without the mutation.
To assess the miR-125b contribution to the PTC relapse prediction, we studied the miRNA expression in groups according to the criteria of the American Thyroid Association (ATA, American Thyroid Association). In patients with an intermediate risk of recurrence, the miR 125-b expression was significantly increased by 107.09 and 11.78 times compared with the low-and high-risk groups, respetively.

Discussion
Autophagy is a universal process of the cellular viability as well as cancer cells. For nine key autophagy associated proteins (ULK1, ATG16L1, MAP1LC3B, PRKAA1, ATG12, ATG 14, ATG5, mTOR and BECN1), we found five hundred sixty-two microRNAs (in Supplements). Among them only sixty-nine are found to be correlated with PTC. However, among the targets associated with autophagy, only miR-125b is able to regulate transcription and translation of the ATG5 protein, mediating the resistance to therapy in PTC [19].
Non-coding RNA sequences, such as miRNAs, are essential for the autophagic flux modification, from upstream signaling pathways to later stages of autophagy. We found the theoretical and experimental verification of hsa-miR-125b implication in thyroid cancer progression. The revealed forty-four genes were found to be targets for miR-125b. KEEG data analysis and the DIANA microRNA investigation tool show that a variety of genes and transcriptional factors guided the PTC progression. But only few of them are under consideration in PTC research. We indicated four target proteins associated with PTC. They belong to the two signaling cascades that contribute to the cancer progression. MAPK1 and MAPK3 are the components of the rout that orchestrate the rapid cancer cell division and determine the aggressive cancer cells behavior [63]. PIK3CA activates AKT/mTOR signaling and operates the protein's metabolism as a key mechanism of cells adaptation and survival [64]. The novel genetic marker is CCDC6, the components of the gene's translocation in PTC [65]. It should be noted that the role of most molecular and genetic factors is still unknown, which causes the unpredictable pattern of behavior of tumor cells and leads to the low treatment response.
Currently, the diagnostic application of microRNAs in human tumor diseases, as well as for recurrence detection after therapy, is widely discussed. Multidirectional changes in the miR-125b level in the PTC were revealed. The presented data indicate the miRNAs involvement in oncogenesis. There was a decrease in the hsa-miR-125b rate found in patients with PTC compared with benign pathology tissues. But the PTC classical subtype was associated with an increased level of expression.
It is known that the tumor invasion is the cancer progression, associated with the intracellular signaling, transcription and growth factors changes [12,37,65] and reflected in the aggressive tumor behavior. In addition, a high hsa-miR-125b expression was associated with the BRAF V600E status, which was disclosed in the results of the experiment on melanoma in cell culture [17]. The study showed an association between the cancer mutational status and hsa-miR-125b expression.
The tumor invasive potential may be associated with the autophagy, modulating the transformed cell behavior [18]. Activation of "cell self-digestion" is significant in the PTC progression, including the therapy effectiveness. It was found that miRNA expression differs depending on the disease recurrence risk. The maximum values of the expression were recorded for patients with an intermediate risk, which is probably of decisive importance, since tumors with an intermediate risk of relapses may have biological features associated with an aggressive phenotype. In previous research, the LC3B protein content, an autophagosome protein, was noted in this group of patients [86]. Some components of the autophagic machinery have important non-autophagic cellular functions. High MAP1LC3B levels could result from cancer progression and reflect the low sensitivity of anti-cancer therapy [87].

Conclusions
Five hundred and sixty-seven miRNAs associated with nine autophagy proteins (ULK1, ATG16L1, MAP1LC3B, PRKAA1, ATG12, ATG 14, ATG5, mTOR and BECN1) were identified, of which sixty-two are accompanied by thyroid cancer, and only miR-125b is associated with the autophagy modulation in this pathology.
The involvement of the non-coding RNA sequences remains unclear. It found the multiple genes and transcription factors used in oncogenesis. The concept of the paper was to find the pool of miRNAs associated with autophagy, and then identify the most powerful ones that regulate the oncogenesis in the thyroid gland. The interesting trend was again verified for the thyroid cancers, that determination of one isolated single pathway and route does not explain and predict cancer behavior. There should be found the molecular factors combination that are essential for the oncogenesis, providing new prominent targets in search for anti-cancer therapy.
Autophagy is a universal, powerful tool for thyroid cancer progression. Our findings revealed the hsa-miR-125b expression changes in thyroid cancer invasion, in BRAF V600E mutant tumors and in intermediate recurrence risk. High expression of hsa-miR-125b is a feature of the classical tumor subtype, as well as for tumors with invasion into the thyroid gland capsule. The data obtained indicate the role of miRNAs in the formation of the tumor invasive properties. The role of miRNAs is planned to be investigated further by a more comprehensive study.

Informed Consent Statement:
Written informed consent has been obtained from the patient(s) to publish this paper. Data Availability Statement: All the relevant data have been provided in the manuscript, and any supplementary datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest:
The authors affirm that they have no competing interest.