MiR-26a-5p Inhibits Cell Proliferation and Enhances Doxorubicin Sensitivity in HCC Cells via Targeting AURKA

miR-26a is a Key Therapeutic Target with Enormous Potential in the Diagnosis and Prognosis of Human Disease

MicroRNA-26a (miR-26a) belongs to small non-coding regulatory RNA molecules emerging as fundamental post-transcriptional regulators inhibiting gene expression that plays vital roles in various processes of human diseases such as depression, renal ischemia and reperfusion injury, liver injury and some refractory cancer. In this review, we expound on the results of studies about miR-26a with emphasis on its function in animal models or in vitro cell culture to simulate the most common human disease in the clinic. Furthermore, we also illustrate the underlying mechanisms of miR-26a in strengthening the antitumor activity of antineoplastic drugs. Importantly, dysregulation of miR-26a has been related to many chronic and malignant diseases, especially in neurological disorders in the brain such as depression and neurodegenerative diseases as well as cancers such as papillary thyroid carcinoma, hepatocellular carcinoma and so on. It follows that miR-26a has a strong possibility to be a potential therapeutic target for the treatment of neurological disorders and cancers. Although the research of miRNAs has made great progress in the last few decades, much is yet to be discovered, especially regarding their underlying mechanisms and roles in the complex diseases of humans. Consequently, miR-26a has been analyzed in chronic and malignant diseases, and we discuss the dysregulation of miR-26a and functional roles in the development and pathogenesis of these diseases, which is very helpful for understanding their mechanisms as new biomarkers for diagnosing and curing diseases in the near future.

Machine learning and experimental screening of chromatin regulator signatures and potential drugs in hepatitis B related hepatocellular carcinoma

Many evidences have confirmed that chromatin regulator factors (CRs) are involved in the progression of cancer, but its potential mechanism of affecting hepatitis B related hepatocellular carcinoma still needs to be studied. Our study detected the CRs that affect hepatitis B related hepatocellular carcinoma (HBV-HCC) through machine learning analysis, conducted the analysis of immune cells, constructed the relevant risk model and immune function infiltration, and predicted the potential therapeutic drugs. We found that these CRs were significantly related to the immune cells of Macrophages, B cells, CD8+T cells, etc., and PBK, AURKA, TOP2A and AURKB were the potential risk CRs of HBV-HCC. The expression levels of these four CRs increased in HepG2.2.15 cells and the liver of HBV-HCC patients, consistent with the predicted risk model. Subsequently, ten potential drugs closely related to the risk CRs were finally obtained, experimental research on resveratrol has shown that it can inhibit the proliferation of HepG2.2.15 cells and potentially inhibit the occurrence and development of HBV-HCC. Our study provides novel insights into the function of CRs in HBV-HCC and certain ideas for more accurate targeted therapy.Communicated by Ramaswamy H. Sarma.

Synthetic miR-26a mimics delivered by tumor exosomes repress hepatocellular carcinoma through downregulating lymphoid enhancer factor 1

BACKGROUND

The dysregulation of exosomal microRNAs plays an important role in the progression of hepatocarcinogenesis. In this study, we investigated the therapeutic potential of synthetic exosomal miR-26a against HCC cells and explored the feasibility of tumor-derived exosomes as drug delivery vehicles.

METHODS

Proliferation and migration assays were performed to examine the effects of miR-26a on HCC in vitro. The direct target gene of miR-26a was identified through miRecords analysis and target validation. The transferring efficiency and anti-HCC effect of exosomes with different origin were studied and the optimal miR-26a delivery method was established and verified in vitro and in vivo. In addition, the relationships between prognosis of HCC patients and miR-26a expression in HCC serum and exosomes were retrospectively analyzed.

RESULTS

Here, we found that tumor cell-derived exosomes were taken in preferentially by HCC cells and promoted HCC progression through Wnt pathway by low-density lipoprotein receptor-related protein 6 (LRP6). HCC cells with vacuolar protein sorting-associated protein 35 knocked down were adopted to generate engineered LRP6-exosomes. The engineered HCC-derived exosomes loading miR-26a inhibited HCC progression in vitro and in vivo effectively. Overexpression of miR-26a impaired the growth and migration of HCC by targeting lymphoid enhancer factor 1 (LEF1). Moreover, low expression of exosomal miR-26a was an independent prognostic factor for recurrence and survival in HCC patients.

CONCLUSIONS

Our findings suggested the exosomal miR-26a could serve as a non-invasive prognostic marker for HCC patients. Genetically modified tumor-derived exosomes showed preferable transfection efficiency but reduced Wnt activity, which provides a novel therapeutic strategy for HCC.

Doxorubicin and other anthracyclines in cancers: Activity, chemoresistance and its overcoming

Anthracyclines have been important and effective treatments against a number of cancers since their discovery. However, their use in therapy has been complicated by severe side effects and toxicity that occur during or after treatment, including cardiotoxicity. The mode of action of anthracyclines is complex, with several mechanisms proposed. It is possible that their high toxicity is due to the large set of processes involved in anthracycline action. The development of resistance is a major barrier to successful treatment when using anthracyclines. This resistance is based on a series of mechanisms that have been studied and addressed in recent years. This work provides an overview of the anthracyclines used in cancer therapy. It discusses their mechanisms of activity, toxicity, and chemoresistance, as well as the approaches used to improve their activity, decrease their toxicity, and overcome resistance.

A systematic study on the treatment of hepatitis B-related hepatocellular carcinoma with drugs based on bioinformatics and key target reverse network pharmacology and experimental verification

BACKGROUND

Chronic hepatitis B virus (HBV) infection is the major etiology of hepatocellular carcinoma (HCC). However, the mechanism of hepatitis B-related hepatocellular carcinoma (HBV-related HCC) is still unclear. Therefore, understanding the pathogenesis and searching for drugs to treat HBV-related HCC was an effective strategy to treat this disease.

PURPOSE

Bioinformatics was used to predict the potential targets of HBV-related HCC. The reverse network pharmacology of key targets was used to analyze the clinical drugs, traditional Chinese medicine (TCM) and small molecules of TCM in the treatment of HBV-related HCC.

METHODS

In this study, three microarray datasets totally containing 330 tumoral samples and 297 normal samples were selected from the GEO database. These microarray datasets were used to screen DEGs. And the expression profile and survival of 6 key genes were analyzed. In addition, Comparative Toxicogenomics Database and Coremine Medical database were used to enrich clinical drugs and TCM of HBV-related HCC by the 6 key targets. Then the obtained TCM were classified based on the Chinese Pharmacopoeia. Among these top 6 key genes, CDK1 and CCNB1 had the most connection nodes and the highest degree and were the most significantly expressed. In general, CDK1 and CCNB1 tend to form a complex, which is conducive to cell mitosis. Hence, this study mainly studied CDK1 and CCNB1. HERB database was used to predict small molecules TCM. The inhibition effect of quercetin, celastrol and cantharidin on HepG2.2.15 cells and Hep3B cells was verified by CCK8 experiment. The effects of quercetin, celastrol and cantharidin on CDK1 and CCNB1 of HepG2.2.15 cells and Hep3B cells were determined by Western Blot.

RESULTS

In short, 272 DEGs (53 upregulated and 219 downregulated) were identified. Among these DEGs, 6 key genes with high degree were identified, which were AURKA, BIRC5, CCNB1, CDK1, CDKN3 and TYMS. Kaplan-Meier plotter analysis showed that higher expression levels of AURKA, BIRC5, CCNB1, CDK1, CDKN3 and TYMS were associated with poor OS. According to the first 6 key targets, a variety of drugs and TCM were identified. These results showed that clinical drugs included targeted drugs, such as sorafenib, palbociclib and Dasatinib. and chemotherapy drugs, such as cisplatin and doxorubicin. TCM, such as the TCM flavor was mainly warm and bitter, and the main meridians were liver and lung. Small molecules of TCM included flavonoids, terpenoids, alkaloids and glycosides, such as quercetin, celastrol, cantharidin, hesperidin, silymarin, casticin, berberine and ursolic acid, which have great potential in anti-HBV-related HCC. For molecular docking of chemical components, the molecules with higher scores were flavonoids, alkaloids, etc. Three representative types of TCM small molecules were verified respectively, and it was found that quercetin, celastrol and cantharidin inhibited the proliferation of HepG2.2.15 cells and Hep3B cells along concentration gradient. Quercetin, celastrol and cantharidin decreased CDK1 expression in HepG2.2.15 and Hep3B cells, but for CCNB1, only cantharidin decreased CCNB1 expression in the two strains of cells.

CONCLUSION

In conclusion, AURKA, BIRC5, CCNB1, CDK1, CDKN3 and TYMS could be potential targets for the diagnosis and prognosis of HBV-related HCC. Clinical drugs include chemotherapeutic and targeted drug, traditional Chinese medicine is mainly bitter and warm TCM. Small molecular of TCM including flavonoids, terpenoids and glycosides and alkaloids, which have great potential in anti-HBV-related HCC. This study provides potential therapeutic targets and novel strategies for the treatment of HBV-related HCC.

Differential translation of mRNA isoforms underlies oncogenic activation of cell cycle kinase Aurora A

Aurora Kinase A (AURKA) is an oncogenic kinase with major roles in mitosis, but also exerts cell cycle- and kinase-independent functions linked to cancer. Therefore, control of its expression, as well as its activity, is crucial. A short and a long 3'UTR isoform exist for AURKA mRNA, resulting from alternative polyadenylation (APA). We initially observed that in triple-negative breast cancer, where AURKA is typically overexpressed, the short isoform is predominant and this correlates with faster relapse times of patients. The short isoform is characterized by higher translational efficiency since translation and decay rate of the long isoform are targeted by hsa-let-7a tumor-suppressor miRNA. Additionally, hsa-let-7a regulates the cell cycle periodicity of translation of the long isoform, whereas the short isoform is translated highly and constantly throughout interphase. Finally, disrupted production of the long isoform led to an increase in proliferation and migration rates of cells. In summary, we uncovered a new mechanism dependent on the cooperation between APA and miRNA targeting likely to be a route of oncogenic activation of human AURKA.

The role of AURKA/miR-199b-3p in hepatocellular carcinoma cells

BACKGROUND

Previous studies proved that AURKA functions as an oncogene in several cancers. This article aimed to probe the miRNA-induced regulatory mechanism of AURKA in hepatocellular carcinoma (HCC).

METHODS

Differentially expressed genes in TCGA-LIHC dataset were screened by bioinformatics methods. Levels of miR-199b-3p and AURKA mRNA were examined by qRT-PCR. Western blot was utilized to evaluate protein levels of AURKA, p-AKT, and AKT. Dual-luciferase assay was introduced to explore their interaction. MTT, colony formation, scratch healing, transwell, and flow cytometry assays were introduced into cell proliferation, migration, invasion, and apoptosis assessment. The impact of miR-199b-3p/AURKA axis on HCC tumor growth was determined in a tumor xenograft model.

RESULTS

We found that AURKA was highly expressed in HCC and was coupled to poor prognosis of HCC. As manifested by cellular assays, compared to the normal cells HL-7702, AURKA presented notably high expression in HCC cell lines. Overexpressed AURKA evidently impelled the proliferation, colony formation, migration, and invasion of HCC cells while suppressing apoptosis. The regulatory gene upstream of AURKA was predicted to be miR-199b-3p by bioinformatics method, and there was a markedly negative correlation between the two. Overexpressed miR-199b-3p constrained HCC cell proliferation, migration, and invasion while fostering apoptosis, which could be counteracted by upregulating AURKA. MiR-199b-3p repressed the tumor growth in vivo by targeting AURKA and affected PI3K/AKT signaling pathway.

CONCLUSION

To summarize, this study implied the regulatory mechanism of miR-199b-3p/AURKA axis in HCC, and supplied optional therapeutic targets for HCC patients.

Regulating the regulator: a survey of mechanisms from transcription to translation controlling expression of mammalian cell cycle kinase Aurora A

Aurora Kinase A (AURKA) is a positive regulator of mitosis with a strict cell cycle-dependent expression pattern. Recently, novel oncogenic roles of AURKA have been uncovered that are independent of the kinase activity and act within multiple signalling pathways, including cell proliferation, survival and cancer stem cell phenotypes. For this, cellular abundance of AURKA protein is per se crucial and must be tightly fine-tuned. Indeed, AURKA is found overexpressed in different cancers, typically as a result of gene amplification or enhanced transcription. It has however become clear that impaired processing, decay and translation of AURKA mRNA can also offer the basis for altered AURKA levels. Accordingly, the involvement of gene expression mechanisms controlling AURKA expression in human diseases is increasingly recognized and calls for much more research. Here, we explore and create an integrated view of the molecular processes regulating AURKA expression at the level of transcription, post-transcription and translation, intercalating discussion on how impaired regulation underlies disease. Given that targeting AURKA levels might affect more functions compared to inhibiting the kinase activity, deeper understanding of its gene expression may aid the design of alternative and therapeutically more successful ways of suppressing the AURKA oncogene.

Exploring the role of non-coding RNAs in autophagy

As a self-degradative mechanism, macroautophagy/autophagy has a role in the maintenance of energy homeostasis during critical periods in the development of cells. It also controls cellular damage through the eradication of damaged proteins and organelles. This process is accomplished by tens of ATG (autophagy-related) proteins. Recent studies have shown the involvement of non-coding RNAs in the regulation of autophagy. These transcripts mostly modulate the expression of ATG genes. Both long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) have been shown to modulate the autophagy mechanism. Levels of several lncRNAs and miRNAs are altered in this process. In the present review, we discuss the role of lncRNAs and miRNAs in the regulation of autophagy in diverse contexts such as cancer, deep vein thrombosis, spinal cord injury, diabetes and its complications, acute myocardial infarction, osteoarthritis, pre-eclampsia and epilepsy.Abbreviations: AMI: acute myocardial infarction; ATG: autophagy-related; lncRNA: long non-coding RNA; miRNA: microRNA.

miR-26a-5p suppresses nasopharyngeal carcinoma progression by inhibiting PTGS2 expression

Nasopharyngeal carcinoma (NPC) has a low five-year survival rate, and its pathogenesis remains unclear. There is an urgent need to improve our understanding of the genetic regulation of NPC tumorigenesis and development. The role of miR-26a-5p in NPC growth regulation and the expression of its target, PTGS2, was analyzed. Quantitative Real-time PCR assay was used to detect miR-26a-5p and PTGS2 expression in human NPC tissues and cell lines. The RNA pull-down dual-luciferase reporter assay was used to determine the association between miR-26a-5p and PTGS2. The effects of miR-26a-5p and PTGS2 on NPC cell viability, proliferation, migration, and invasion were measured by CCK-8, BrdU, and Transwell assays. miR-26a-5p expression in NPC tissues and cell lines was significantly decreased. The overexpression of miR-26a-5p inhibited the viability, proliferation, migration, and invasion of NPC cells. miR-26a-5p bound to the 3-'untranslated region of PTGS2, thus reducing PTGS2 protein levels. miR-26a-5p inhibited NPC development by reducing the expression of its target PTGS2.

The related miRNAs involved in doxorubicin resistance or sensitivity of various cancers: an update

Doxorubicin (DOX) is an effective chemotherapy agent against a wide variety of tumors. However, intrinsic or acquired resistance diminishes the sensitivity of cancer cells to DOX, which leads to a cancer relapse and treatment failure. Resolutions to this challenge includes identification of the molecular pathways underlying DOX sensitivity/resistance and the development of innovative techniques to boost DOX sensitivity. DOX is classified as a Topoisomerase II poison, which is cytotoxic to rapidly dividing tumor cells. Molecular mechanisms responsible for DOX resistance include effective DNA repair and resumption of cell proliferation, deregulated development of cancer stem cell and epithelial to mesenchymal transition, and modulation of programmed cell death. MicroRNAs (miRNAs) have been shown to potentiate the reversal of DOX resistance as they have gene-specific regulatory functions in DOX-responsive molecular pathways. Identifying the dysregulation patterns of miRNAs for specific tumors following treatment with DOX facilitates the development of novel combination therapies, such as nanoparticles harboring miRNA or miRNA inhibitors to eventually prevent DOX-induced chemoresistance. In this article, we summarize recent findings on the role of miRNAs underlying DOX sensitivity/resistance molecular pathways. Also, we provide latest strategies for utilizing deregulated miRNA patterns as biomarkers or miRNAs as tools to overcome chemoresistance and enhance patient's response to DOX treatment.

Ascites in ovarian cancer: MicroRNA deregulations and their potential roles in ovarian carcinogenesis

Ovarian cancer comprises the most lethal gynecologic malignancy and is accompanied by the high potential for the incidence of metastasis, recurrence and chemotherapy resistance, often associated with a formation of ascitic fluid. The differentially expressed ascites-derived microRNAs may be linked to ovarian carcinogenesis. The article focuses on a number of miRNAs that share a common expression pattern as determined by independent studies using ascites samples and with regard to their functions and outcomes in experimental and clinical investigations.Let-7b and miR-143 have featured as tumor suppressors in ovarian cancer, which is in line with data on other types of cancer. Although two miRNAs, i.e. miR-26a-5p and miR-145-5p, act principally as tumor suppressor miRNAs, they occasionally exhibit oncogenic roles. The performance of miR-95-3p, upregulated in ascites, is open to debate given the current lack of supportive data on ovarian cancer; however, data on other cancers indicates its probable oncogenic role. Different findings have been reported for miR-182-5p and miR-200c-3p; in addition to their presumed oncogenic roles, contrasting findings have indicated their ambivalent functions. Further research is required for the identification and evaluation of the potential of specific miRNAs in the diagnosis, prediction, treatment and outcomes of ovarian cancer patients.

The role of miRNAs and lncRNAs in conferring resistance to doxorubicin

Doxorubicin is a chemotherapeutic agent that inhibits topoisomerase II, intercalates within DNA base pairs and results in oxidative DNA damage, thus inducing cell apoptosis. Although it is effective in the treatment of a wide range of human cancers, the emergence of resistance to this drug can increase tumour growth and impact patients' survival. Numerous molecular mechanisms and signalling pathways have been identified that induce resistance to doxorubicin via stimulation of cell proliferation, cell cycle switch and preclusion of apoptosis. A number of microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) have also been identified that alter sensitivity to doxorubicin. Understanding the particular impact of these non-coding RNAs in conferring resistance to doxorubicin has considerable potential to improve selection of chemotherapeutic regimens for cancer patients. Moreover, modulation of expression of these transcripts is a putative strategy for combating resistance. In the current paper, the influence of miRNAs and lncRNAs in the modification of resistance to doxorubicin is discussed.

The regulatory effect of has-circ-0001146/miR-26a-5p/MNAT1 network on the proliferation and invasion of osteosarcoma

Osteosarcoma is a malignant bone tumour with the lowest survival rates out of all paediatric cancers and is primarily diagnosed in children and adolescents. MNAT1 is a subunit in the cyclin-dependent kinase-activating kinase complex. Abnormal up-regulation of MNAT1 has been associated with the poor prognosis of multiple cancers. Bioinformatics analysis showed that has-circ-0001146 and miR-26a-5p were involved in the regulation of MNAT1 in osteosarcoma. The present study investigated the regulatory effects of has-circ-0001146 and miR-26a-5p on MNAT1 expression using luciferase reporter and RNA-pull down assays. The effects of the has-circ-0001146/miR26a-5p/Mnat1 network on the proliferation and invasion of osteosarcoma were evaluated by cell viability, apoptosis, migration, and invasion assays. Osteosarcoma tissues showed higher MNAT1 and has-circ-0001146 expression than adjacent normal tissues, although the expression of MNAT1 was not significantly up-regulated in sarcomas according to TCGA databases. As indicated by luciferase reporter and RNA-pull down assays, miR-26a-5p was able to bind to both has-circ-0001146 and MNAT1 mRNA. The depletion of has-circ-0001146 as well as the increase of miR-26a-5p decreased MNAT1 expression in osteosarcoma cells, while the reduction of miR-26a-5p was associated with increased MNAT1 expression. These data suggested that has-circ-0001146 promoted MNAT1 expression by competitively binding to miR-26a-5p with MNAT1 mRNA. The depletion of has-circ-0001146 or MNAT1 or the increase of miR-26a-5p inhibited osteosarcoma cell viability and invasion, and increased apoptosis. Reduction of miR-26a-5p conversely promoted osteosarcoma cell viability and invasion. The present study confirmed that has-circ-0001146 blocked miR-26a-5p targeting MNAT1 in osteosarcoma cells, thereby promoting the malignant behaviours of osteosarcoma cells.

MicroRNA signature in liver cancer

Liver cancer is the 7^th utmost frequent neoplasm and the 4^th principal source of cancer deaths. This malignancy is linked with several environmental and lifestyle-related factors emphasizing the role of epigenetics in its pathogenesis. MicroRNAs (miRNAs) have been regarded as potent epigenetic mechanisms partaking in the pathogenesis of liver cancer. Dysregulation of miRNAs has been related with poor outcome of patients with liver cancer. In the current manuscript, we provide a concise review of the results of recent studies about the role of miRNAs in the progression of liver cancer and their diagnostic and prognostic utility.

Small-sized extracellular vesicles (EVs) derived from acute myeloid leukemia bone marrow mesenchymal stem cells transfer miR-26a-5p to promote acute myeloid leukemia cell proliferation, migration, and invasion

Bone marrow mesenchymal stem cells (BMSCs) in acute myeloid leukemia (AML) microenvironment undergo modification that includes expression of contents in the small-sized extracellular vesicles (EVs) they secrete. This study aims to investigate whether small-sized EVs from BMSCs of AML patients regulate AML progression by modifying the expression of miR-26a-5p. Small-sized EVs from BMSCs of AML patients (AML-BMSC-EVs) or healthy controls (HC-BMSC-EVs) were isolated by ultra-centrifugation and administered to AML cells (OCI/AML-2 and THP-1). Cell proliferation, migration, and invasion were evaluated by CCK-8 assay, Transwell migration and invasion assays, respectively. Compared with HC-BMSC-EVs, AML-BMSC-EVs contained higher expression of miR-26a-5p and promoted AML cell proliferation, migration, and invasion. Inhibition of miR-26a-5p expression in AML-BMSC-EVs could abrogate the promoting effects of AML-BMSC-EVs on AML cell proliferation, migration, and invasion. Furthermore, GSK3β was a direct target of miR-26a-5p. Moreover, AML-BMSC-EVs inhibited GSK3β expression and activated Wnt/β-catenin signaling in AML cells. Additionally, GSK3β overexpression in THP-1 cells counteracted the promoting effects of AML-BMSCs-EVs on THP-1 cell proliferation, migration, and invasion. AML-BMSC-EVs promoted AML progression by transferring miR-26a-5p to AML cells and subsequently activating the Wnt/β-catenin pathway.

Comparison of MicroRNA Transcriptomes Reveals the Association between MiR-148a-3p Expression and Rumen Development in Goats

Simple Summary

In ruminants, the rumen epithelium plays an important role in nutrient absorption, metabolism and transport. MicroRNAs (miRNAs) have been reported to regulate the proliferation of diverse epithelial cells. In this study, we profiled the miRNA transcriptomes of goat rumens at four development stages and screened for candidate miRNAs related to rumen development. MiR-148a-3p was found to be highly expressed in the rumen tissues and induced the proliferation of GES-1 cells by targeting QKI. Our findings provide some insights into the functional roles of miRNAs in rumen growth and functional development in ruminants.

Abstract

The rumen is an important digestive organ of ruminants. From the fetal to adult stage, the morphology, structure and function of the rumen change significantly. However, the knowledge of the intrinsic genetic regulation of these changes is still limited. We previously reported a genome-wide expression profile of miRNAs in pre-natal goat rumens. In this study, we combined and analyzed the transcriptomes of rumen miRNAs during pre-natal (E60 and E135) and post-natal (D30 and D150) stages. A total of 66 differentially expressed miRNAs (DEMs) were identified in the rumen tissues from D30 and D150 goats. Of these, 17 DEMs were consistently highly expressed in the rumens at the pre-weaning stages (E60, E135 and D30), while down-regulated at D150. Noteworthy, annotation analysis revealed that the target genes regulated by the DEMs were mainly enriched in MAPK signaling pathway, Jak-STAT signaling pathway and Ras signaling pathway. Interestingly, the expression of miR-148a-3p was significantly high in the embryonic stage and down-regulated at D150. The potential binding sites of miR-148a-3p in the 3′-UTR of QKI were predicted by the TargetScan and verified by the dual luciferase report assay. The co-localization of miR-148a-3p and QKI through in situ hybridization was observed in the rumen tissues but not in the intestinal tracts. Moreover, the expression of miR-148a-3p in the epithelium was significantly higher than that in the other layers of the rumen, suggesting that miR-148a-3p is involved in the development of the rumen epithelial cells by targeting QKI. Subsequently, miR-148a-3p inhibitor was found to induce the proliferation of GES-1 cells. Taken together, our study identified DEMs involved in the development of the rumen and provides insights into the regulation mechanism of rumen development in goats.