miRNA-mRNA Regulatory Network and Factors Associated with Prediction of Prognosis in Hepatocellular Carcinoma

A protective role of nintedanib in peritoneal fibrosis through H19-EZH2-KLF2 axis via impeding mesothelial-to-mesenchymal transition

BACKGROUND

Peritoneal fibrosis (PF), a common complication of long-term peritoneal dialysis, accounts for peritoneal ultrafiltration failure to develop into increased mortality. Nintedanib has previously been shown to protect against multi-organ fibrosis, including PF. Unfortunately, the precise molecular mechanism underlying nintedanib in the pathogenesis of PF remains elusive.

METHODS

The mouse model of PF was generated by chlorhexidine gluconate (CG) injection with or without nintedanib administration, either with the simulation for the cell model of PF by constructing high-glucose (HG)-treated human peritoneal mesothelial cells (HPMCs). HE and Masson staining were applied to assess the histopathological changes of peritoneum and collagen deposition. FISH, RT-qPCR, western blot and immunofluorescence were employed to examine distribution or expression of targeted genes. Cell viability was detected using CCK-8 assay. Cell morphology was observed under a microscope. RNA immunoprecipitation (RIP) and chromatin immunoprecipitation (ChIP) assays were applied to validate the H19-EZH2-KLF2 regulatory axis.

RESULTS

Aberrantly overexpressed H19 was observed in both the mouse and cell model of PF, of which knockdown significantly blocked HG-induced mesothelial-to-mesenchymal transition (MMT) of HPMCs. Moreover, loss of H19 further strengthened nintedanib-mediated suppressive effects against MMT process in a mouse model of PF. Mechanistically, H19 could epigenetically repressed KLF2 via recruiting EZH2. Furthermore, TGF-β/Smad pathway was inactivated by nintedanib through mediating H19/KLF2 axis.

CONCLUSION

In summary, nintedanib disrupts MMT process through regulating H19/EZH2/KLF2 axis and TGF-β/Smad pathway, which laid the experimental foundation for nintedanib in the treatment of PF.

CHST4 Gene as a Potential Predictor of Clinical Outcome in Malignant Pleural Mesothelioma

Malignant pleural mesothelioma (MPM) develops primarily from asbestos exposures and has a poor prognosis. In this study, The Cancer Genome Atlas was used to perform a comprehensive survival analysis, which identified the CHST4 gene as a potential predictor of favorable overall survival for patients with MPM. An enrichment analysis of favorable prognostic genes, including CHST4, showed immune-related ontological terms, whereas an analysis of unfavorable prognostic genes indicated cell-cycle-related terms. CHST4 mRNA expression in MPM was significantly correlated with Bindea immune-gene signatures. To validate the relationship between CHST4 expression and prognosis, we performed an immunohistochemical analysis of CHST4 protein expression in 23 surgical specimens from surgically treated patients with MPM who achieved macroscopic complete resection. The score calculated from the proportion and intensity staining was used to compare the intensity of CHST4 gene expression, which showed that CHST4 expression was stronger in patients with a better postoperative prognosis. The median overall postoperative survival was 107.8 months in the high-expression-score group and 38.0 months in the low-score group (p = 0.044, log-rank test). Survival after recurrence was also significantly improved by CHST4 expression. These results suggest that CHST4 is useful as a prognostic biomarker in MPM.

Inhibition of hepatocellular carcinoma growth via modulation of the miR-221/SOX11 axis by curcumin and berberine

Hepatocellular carcinoma (HCC) is a fatal malignancy that has limited treatment options. This study focused on the potential therapeutic effects of curcumin (CUR) and berberine (BBR) on the miR-221/SRY-box transcription factor 11 (SOX11) axis in HCC. We investigated the combined effects of CUR and BBR on HEPG2 and Huh7 cell survival and miR-221 expression using Cell Counting Kit-8 assays and RT-qPCR, respectively. Western blotting was used to detect changes in the apoptosis-related caspase-3/9 protein levels. We performed bioinformatics analysis and dual-luciferase assays and measured apoptotic protein levels to assess the role of the miR-221/SOX11 axis in mediating the effects of CUR-BBR. Both CUR and BBR suppressed HCC cell growth in a dose-dependent manner, with the most potent combined effect observed at a 2:1 ratio. CUR-BBR treatment significantly downregulated miR-221 expression, and miR-221 overexpression partially reversed the CUR-BBR-mediated decrease in cell survival. In addition, SOX11 was found to be a direct target of miR-221. CUR-BBR treatment upregulated SOX11 expression, and overexpression of SOX11 restored the inhibitory effects of CUR-BBR on cell growth, migration, and invasion and promoted apoptosis in the presence of miR-221. Furthermore, CUR-BBR activated pro-apoptotic proteins caspase-3/9 through the miR-221/SOX11 axis. The combined effect of CUR-BBR played an important role in inhibiting the growth of HCC cells. This combined effect was achieved by regulating the miR-221/SOX11 axis and activating the synthesis of pro-apoptotic proteins. Our findings highlight a promising combined therapeutic approach for HCC and underscore the importance of targeting the miR-221/SOX11 axis.

Knocking Down Gm16685 Decreases Liver Granuloma in Murine Schistosomiasis Japonica

Long noncoding RNAs (lncRNAs) can regulate key genes and pathways in liver disease development. Moreover, macrophages are speculated to play an important role in regulating granulomatous inflammation during schistosomiasis. However, the role of lncRNAs in the formation of liver granulomas by influencing the polarization of macrophages in Schistosoma japonicum infection is unclear. Our study aimed to determine whether lncRNAs can play a role in S. japonicum-induced hepatic egg granulomas and elucidate their effect on macrophages. We established S. japonicum infection models and screened the target lncRNA Gm16685 highly expressed in schistosomiasis mice using high-throughput sequencing. Hematoxylin and eosin staining revealed that the knockdown of Gm16685 reduced the area of egg granulomas. Moreover, M1 macrophage factor genes were significantly downregulated in Gm16685 knockdown livers. Meanwhile, M2 macrophage factor genes were significantly upregulated, which was consistent with the protein detection results. Hepatocytes, hepatic stellate cells, and macrophages were isolated from mouse models infected with S. japonicum, with Gm16685 being significantly upregulated in macrophages. Moreover, the knockdown of Gm16685 in RAW264.7 cells revealed similar results to in liver tissue. RNA fluorescence in situ hybridization (FISH) and nucleocytoplasmic separation experiments revealed that Gm16685 was predominantly localized in the cytoplasm of cells. We found that miR-205-5p was upregulated after Gm16685 was knocked down. After overexpression of miR-205-5p, the expression of Gm16685 and inflammatory factors was significantly downregulated. These results indicate that Gm16685 can participate in the pathogenesis of hepatic disease in schistosomiasis and promote M1 macrophage polarization by regulating miR-205-5p. Thus, our study may provide a new target for schistosomiasis japonica treatment.

Evolution and function of developmentally dynamic pseudogenes in mammals

BACKGROUND

Pseudogenes are excellent markers for genome evolution, which are emerging as crucial regulators of development and disease, especially cancer. However, systematic functional characterization and evolution of pseudogenes remain largely unexplored.

RESULTS

To systematically characterize pseudogenes, we date the origin of human and mouse pseudogenes across vertebrates and observe a burst of pseudogene gain in these two lineages. Based on a hybrid sequencing dataset combining full-length PacBio sequencing, sample-matched Illumina sequencing, and public time-course transcriptome data, we observe that abundant mammalian pseudogenes could be transcribed, which contribute to the establishment of organ identity. Our analyses reveal that developmentally dynamic pseudogenes are evolutionarily conserved and show an increasing weight during development. Besides, they are involved in complex transcriptional and post-transcriptional modulation, exhibiting the signatures of functional enrichment. Coding potential evaluation suggests that 19% of human pseudogenes could be translated, thus serving as a new way for protein innovation. Moreover, pseudogenes carry disease-associated SNPs and conduce to cancer transcriptome perturbation.

CONCLUSIONS

Our discovery reveals an unexpectedly high abundance of mammalian pseudogenes that can be transcribed and translated, and these pseudogenes represent a novel regulatory layer. Our study also prioritizes developmentally dynamic pseudogenes with signatures of functional enrichment and provides a hybrid sequencing dataset for further unraveling their biological mechanisms in organ development and carcinogenesis in the future.

Overview of MicroRNAs as Diagnostic and Prognostic Biomarkers for High-Incidence Cancers in 2021

MicroRNAs (miRNAs) are small non-coding RNAs (ncRNAs) about 22 nucleotides in size, which play an important role in gene regulation and are involved in almost all major cellular physiological processes. In recent years, the abnormal expression of miRNAs has been shown to be associated with human diseases including cancer. In the past ten years, the link between miRNAs and various cancers has been extensively studied, and the abnormal expression of miRNAs has been reported in various malignant tumors, such as lung cancer, gastric cancer, colorectal cancer, liver cancer, breast cancer, and prostate cancer. Due to the high malignancy grade of these cancers, it is more necessary to develop the related diagnostic and prognostic methods. According to the study of miRNAs, many potential cancer biomarkers have been proposed for the diagnosis and prognosis of diseases, especially cancer, thus providing a new theoretical basis and perspective for cancer screening. The use of miRNAs as biomarkers for diagnosis or prognosis of cancer has the advantages of being less invasive to patients, with better accuracy and lower price. In view of the important clinical significance of miRNAs in human cancer research, this article reviewed the research status of miRNAs in the above-mentioned cancers in 2021, especially in terms of diagnosis and prognosis, and provided some new perspectives and theoretical basis for the diagnosis and treatment of cancers.

Identification of stage-associated exosome miRNAs in colorectal cancer by improved robust and corroborative approach embedded miRNA-target network

Background

Colorectal cancer (CRC) is a common gastrointestinal tumor with high morbidity and mortality. At the molecular level, patients at different stages present considerable heterogeneity. Although the miRNA in exosome is an effective biomarker to reveal tumor progression, studies based on stage-associated exosome miRNA regulatory network analysis still lacking. This study aims to identify CRC stage-associated exosome miRNAs and reveal their potential function in tumor progression.

Methods

In this study, serum and cellular exosome miRNA expression microarrays associated with CRC were downloaded from GEO database. Stage-common (SC) and stage-specific (SS) differentially expressed miRNAs were extracted and their targets were identified based on 11 databases. Furthermore, miRNA SC and SS regulatory function networks were built based on the CRC phenotypic relevance of miRNA targets, and the corresponding transcription factors were identified. Concurrently, the potential stage-associated miRNAs were identified by receiver-operating characteristic (ROC) curve analysis, survival analysis, drug response analysis, ceRNA analysis, pathway analysis and a comprehensive investigation of 159 publications.

Results

Ten candidate stage-associated miRNAs were identified, with three SC (miR-146a-5p, miR-22-3p, miR-23b-3p) and seven SS (I: miR-301a-3p, miR-548i; IIIA: miR-23a-3p; IV: miR-194-3p, miR-33a-3p, miR-485-3p, miR-194-5p) miRNAs. Additionally, their targets were enriched in several vital cancer-associated pathways such as TGF-beta, p53, and hippo signaling pathways. Moreover, five key hotspot target genes (CCNA2, MAPK1, PTPRD, MET, and CDKN1A) were demonstrated to associated with better overall survival in CRC patients. Finally, miR-23b-3p, miR-301a-3p and miR-194-3p were validated being the most stably expressed stage-associated miRNAs in CRC serum exosomes, cell exosomes and tissues.

Conclusions

These CRC stage-associated exosome miRNAs aid to further mechanism research of tumor progression and provide support for better clinical management in patients with different stages.

The Prognostic Model and Drug Sensitivity of LKB1-Mutant Lung Adenocarcinoma Based on Immune Landscape

Background: Lung cancer is the most common cause of cancer-related deaths worldwide. LKB1-mutant lung adenocarcinoma (LUAD) is a unique subtype of this deadly cancer. LKB1 mutations cause functional changes in a variety of cell processes, including immune functions, that affect prognosis. To date, the potential role of immunity in the prognosis of LKB1-mutant LUAD is not well understood.

Methods: We systematically analyzed immune-related genes in LUAD samples from The Cancer Genome Atlas (TCGA) database. ESTIMATE and CIBERSORT algorithms were used to explore the immune microenvironment. A prognostic risk model was constructed, and prognostic, immune function, drug sensitivity, and model specificity analyses were performed to identify the effectiveness of the model.

Results: Our results showed that LKB1 mutations suppressed immune function in LUAD. A three-gene signature was constructed to stratify patients into two risk groups. The risk score was an independent predictor for overall survival (OS) in multivariate Cox regression analyses [hazard ratio (HR) > 1, p = 0.002]. Receiver operating characteristic (ROC) curve analyses confirmed that the risk score has better performance than clinicopathological characteristics. Functional analysis revealed that the immune status was different between the risk groups. ZM.447439 was an appropriate treatment for the high-risk group of patients. This risk model is only suitable for LKB1-mutant tumors; it performed poorly in LUAD patients with wild-type LKB1.

Conclusion: Our findings indicate the potential role of immunity in LKB1-mutant LUAD, providing novel insights into prognosis and guiding effective immunotherapy.

Dynamic human liver proteome atlas reveals functional insights into disease pathways

Deeper understanding of liver pathophysiology would benefit from a comprehensive quantitative proteome resource at cell type resolution to predict outcome and design therapy. Here, we quantify more than 150,000 sequence‐unique peptides aggregated into 10,000 proteins across total liver, the major liver cell types, time course of primary cell cultures, and liver disease states. Bioinformatic analysis reveals that half of hepatocyte protein mass is comprised of enzymes and 23% of mitochondrial proteins, twice the proportion of other liver cell types. Using primary cell cultures, we capture dynamic proteome remodeling from tissue states to cell line states, providing useful information for biological or pharmaceutical research. Our extensive data serve as spectral library to characterize a human cohort of non‐alcoholic steatohepatitis and cirrhosis. Dramatic proteome changes in liver tissue include signatures of hepatic stellate cell activation resembling liver cirrhosis and providing functional insights. We built a web‐based dashboard application for the interactive exploration of our resource (www.liverproteome.org).

Stanniocalcin 2 (STC2): a universal tumour biomarker and a potential therapeutical target

Stanniocalcin 2 (STC2) is a glycoprotein which is expressed in a broad spectrum of tumour cells and tumour tissues derived from human breast, colorectum, stomach, esophagus, prostate, kidney, liver, bone, ovary, lung and so forth. The expression of STC2 is regulated at both transcriptional and post-transcriptional levels; particularly, STC2 is significantly stimulated under various stress conditions like ER stress, hypoxia and nutrient deprivation. Biologically, STC2 facilitates cells dealing with stress conditions and prevents apoptosis. Importantly, STC2 also promotes the development of acquired resistance to chemo- and radio- therapies. In addition, multiple groups have reported that STC2 overexpression promotes cell proliferation, migration and immune response. Therefore, the overexpression of STC2 is positively correlated with tumour growth, invasion, metastasis and patients' prognosis, highlighting its potential as a biomarker and a therapeutic target. This review focuses on discussing the regulation, biological functions and clinical importance of STC2 in human cancers. Future perspectives in this field will also be discussed.

The Age-Specific Features and Clinical Significance of NRF2 and MAPK10 Expression in HCC Patients

Background

Nuclear factor (erythroid-derived 2)-like 2 (NRF2) functions decline with age; however, cancer cells can hijack its pathways to ensure survival and aggressiveness. Yet, the role of NRF2 in hepatocellular carcinoma (HCC) is rarely investigated in an age-specific manner. This study investigates the expression of NRF2 and its activator (MAPK10) in different age groups of HCC patients, in addition to the age-specific features of NRF2 and MAPK10 interaction and their clinical significance.

Methods

Tumor and near-tumor tissue samples of 181 HCC patients were used to complete a protein expression analysis of NRF2 and MAPK10. Patients’ survival and clinical data were collected for clinical analysis. Global databases (TCGA, ICGC) were used to collect MAPK10 genetic mutation and mRNA expression data in patients with HCC, colorectal, stomach, and pancreatic cancers.

Results

Our findings revealed an increase in NRF2 protein expression but only in younger HCC patients, along with a decline in MAPK10 ability to activate NRF2 in older patients. We also found an increased MAPK10 genetic mutation rate and decreased mRNA expression in older patients. Low MAPK10 and NRF2 expression levels were associated with shorter survival and poorer prognosis due to positive correlation with microvascular invasion, tumor thrombus, elevated AFP levels, and larger tumor size.

Conclusion

NRF2 expression and oxidative stress mechanism in HCC patients are influenced by age. This magnifies the need to consider patients’ age in treatment strategies and guidelines and re-evaluates the application of studies’ age-standardized findings in older patients who are usually excluded from relevant research.

HOTAIR Facilitates Endocrine Resistance in Breast Cancer Through ESR1/miR-130b-3p Axis: Comprehensive Analysis of mRNA-miRNA-lncRNA Network

Background

To summarize the regulatory role of mRNA-miRNA-lncRNA network associated with endocrine therapy resistance (ETR) in breast cancer.

Methods

We analyzed the differentially expressed genes (DEGs), differentially expressed lncRNAs (DELs), and differentially expressed miRNAs (DEMs) in long-term estrogen-deprived (LTED) estrogen receptor (ER)-positive breast cancer cells (LTED MCF7) (modeling relapse on endocrine therapy) and MCF7 cells in the presence of estrogen (E2) (modeling a patient at primary diagnosis) by mining GSE120929 and GSE120930 datasets. The mRNA-miRNA-lncRNA network was constructed by multiple bioinformatic tools. The prognosis of genes from the network was validated in breast cancer patients with following systemic treatment (endocrine therapy) by GEPIA, Kaplan–Meier plotter and UALCAN database.

Results

Totally, 769 DEGs, 33 DEMs, and 10 DELs were selected. The mRNA-miRNA-lncRNA network was established including 60 mRNA nodes, 6 miRNA nodes and 3 lncRNA nodes. A significant module containing 3 nodes and 3 edges was calculated based on the mRNA-miRNA-lncRNA network. The hub genes in the network are ABCG2, ESR1 and GJA1. ESR1/miR-130b-3p/HOTAIR are significantly correlated with the prognosis of breast cancer patients with endocrine therapy.

Conclusion

This study provides a novel ETR-related mRNA-miRNA-lncRNA network. Further, we suggest that ESR1/miR-130b-3p/HOTAIR may be promising targets for clinical treatment of endocrine therapy-resistant breast cancer.

Multi-omics network characterization reveals novel microRNA biomarkers and mechanisms for diagnosis and subtyping of kidney transplant rejection

BACKGROUND

Kidney transplantation is an optimal method for treatment of end-stage kidney failure. However, kidney transplant rejection (KTR) is commonly observed to have negative effects on allograft function. MicroRNAs (miRNAs) are small non-coding RNAs with regulatory role in KTR genesis, the identification of miRNA biomarkers for accurate diagnosis and subtyping of KTR is therefore of clinical significance for active intervention and personalized therapy.

METHODS

In this study, an integrative bioinformatics model was developed based on multi-omics network characterization for miRNA biomarker discovery in KTR. Compared with existed methods, the topological importance of miRNA targets was prioritized based on cross-level miRNA-mRNA and protein-protein interaction network analyses. The biomarker potential of identified miRNAs was computationally validated and explored by receiver-operating characteristic (ROC) evaluation and integrated "miRNA-gene-pathway" pathogenic survey.

RESULTS

Three miRNAs, i.e., miR-145-5p, miR-155-5p, and miR-23b-3p, were screened as putative biomarkers for KTR monitoring. Among them, miR-155-5p was a previously reported signature in KTR, whereas the remaining two were novel candidates both for KTR diagnosis and subtyping. The ROC analysis convinced the power of identified miRNAs as single and combined biomarkers for KTR prediction in kidney tissue and blood samples. Functional analyses, including the latent crosstalk among HLA-related genes, immune signaling pathways and identified miRNAs, provided new insights of these miRNAs in KTR pathogenesis.

CONCLUSIONS

A network-based bioinformatics approach was proposed and applied to identify candidate miRNA biomarkers for KTR study. Biological and clinical validations are further needed for translational applications of the findings.