Alteration of tumor suppressor BMP5 in sporadic colorectal cancer: a genomic and transcriptomic profiling based study

A new TGF-β risk score predicts clinical and immune landscape in colorectal cancer patients

Background

Aberrant TGF-β signaling pathway can lead to invasive phenotype of colorectal cancer (CRC), resulting in poor prognosis. It is pivotal to develop an effective prognostic factor on the basis of TGF-β-related genes to accurately identify risk of CRC patients.

Methods

We performed differential analysis of TGF-β-related genes in CRC patients from databases and previous literature to obtain TGF-β-related differentially expressed genes (TRDEGs). LASSO-Cox regression was utilized to build a CRC prognostic feature model based on TRDEGs. The model was validated using two GEO validation sets. Wilcoxon rank-sum test was utilized to test correlation of model with clinical factors. ESTIMATE algorithm and ssGSEA and tumor mutation burden (TMB) analysis were used to analyze immune landscape and mutation burden of high-risk (HR) and low-risk (LR) groups. CellMiner database was utilized to identify therapeutic drugs with high sensitivity to the feature genes.

Results

We established a six-gene risk prognostic model with good predictive accuracy, which independently predicted CRC patients' prognoses. The HR group was more likely to experience immunotherapy benefits due to higher immune infiltration and TMB. The feature gene TGFB2 could inhibit the efficacy of drugs such as XAV-939, Staurosporine, and Dasatinib, but promote the efficacy of drugs such as CUDC-305 and by-product of CUDC-305. Similarly, RBL1 could inhibit the drug action of Fluphenazine and Imiquimod but promote that of Irofulven.

Conclusion

A CRC risk prognostic signature was developed on basis of TGF-β-related genes, which provides a reference for risk and further therapeutic selection of CRC patients.

Epithelial-mesenchymal interaction protects normal colonocytes from 4-HNE-induced phenotypic transformation

INTRODUCTION

Recent studies have shown that epithelial-stromal interactions could play a role in the development of colorectal cancer. Here, we investigated the role of fibroblasts in the transformation of normal colonocytes induced by 4-HNE.

METHODS

Normal Co colonocytes and nF fibroblasts from the same mouse colon were exposed, in monoculture (m) or coculture (c), to 4-HNE (5 μM) twice weekly for 3 weeks. Gene expression was then analysed and the ability of Co colonocytes to grow in anchorage-independent conditions was tested in soft agar. Fibroblasts previously treated or not with 4-HNE were also seeded in culture inserts positioned above the agar layers to allow paracrine exchanges with colonocytes.

RESULTS

First, 60% of the genes studied were modulated by coculture in Co colonocytes, with notably increased expression of BMP receptors. Furthermore, while 4-HNE increased the ability of monoculture-treated Co colonocytes to form colonies, this effect was not observed in coculture-treated Co colonocytes. Adding a selective BMPR1 inhibitor during the treatment phase abolished the protective effect of coculture. Conversely, addition of a BMP4 agonist to the medium of monoculture-treated Co colonocytes prevented phenotypic transformation by 4-HNE. Second, the presence of nF(m)-HNE fibroblasts during the soft agar assay increased the number and size of Co(m) colonocyte colonies, regardless of whether these cells had been previously treated with 4-HNE in monoculture. For soft agar assays performed with nF(c) and Co(c) cells initially treated in coculture, only the reassociation between Co(c)-HNE and nF(c)-HNE resulted in a small increase in the number of colonies.

CONCLUSIONS

During the exposure phase, the epithelial-mesenchymal interaction protected colonocytes from 4-HNE-induced phenotypic transformation via activation of the BMP pathway. This intercellular dialogue also limited the ability of fibroblasts to subsequently promote colonocyte-anchorage-independent growth. In contrast, fibroblasts pre-exposed to 4-HNE in monoculture strongly increased the ability of Co(m) colonocytes to form colonies.

The Construction of a Multi-Gene Risk Model for Colon Cancer Prognosis and Drug Treatments Prediction

In clinical practice, colon cancer is a prevalent malignant tumor of the digestive system, characterized by a complex and progressive process involving multiple genes and molecular pathways. Historically, research efforts have primarily focused on investigating individual genes; however, our current study aims to explore the collective impact of multiple genes on colon cancer and to identify potential therapeutic targets associated with these genes. For this research, we acquired the gene expression profiles and RNA sequencing data of colon cancer from TCGA. Subsequently, we conducted differential gene expression analysis using R, followed by GO and KEGG pathway enrichment analyses. To construct a protein-protein interaction (PPI) network, we selected survival-related genes using the log-rank test and single-factor Cox regression analysis. Additionally, we performed LASSO regression analysis, immune infiltration analysis, mutation analysis, and cMAP analysis, as well as an investigation into ferroptosis. Our differential expression and survival analyses identified 47 hub genes, and subsequent LASSO regression analysis refined the focus to 23 key genes. These genes are closely linked to cancer metastasis, proliferation, apoptosis, cell cycle regulation, signal transduction, cancer microenvironment, immunotherapy, and neurodevelopment. Overall, the hub genes discovered in our study are pivotal in colon cancer and are anticipated to serve as important biological markers for the diagnosis and treatment of the disease.

PEX26 Functions as a Metastasis Suppressor in Colorectal Cancer

BACKGROUND/AIMS

Aberrant Peroxisomal Biogenesis Factor 26 (PEX26) occurs in multiple cell process. However, the role of PEX26 in colorectal cancer (CRC) development remains unknown. We aimed to study PEX26 expression, regulation, and function in CRC cells.

METHODS

Using the bioinformatic analysis, real-time quantitative PCR, and immunohistochemistry staining, we detected the expression of PEX26 in CRC and normal tissues. We performed functional experiments in vitro to elucidate the effect of PEX26 on CRC cells. We analyzed the RNA-seq data to reveal the downstream regulating network of PEX26.

RESULTS

PEX26 is significantly down-regulated in CRC and its low expression correlates with the poor overall survival of CRC patients. We further demonstrated that PEX26 over-expression inhibits the ability of CRC cell migration, invasion, and epithelial-mesenchymal transition (EMT), while PEX26 knockdown promotes the malignant phenotypes of migration, invasion, and EMT via activating the Wnt pathway.

CONCLUSION

Overall, our results showed that the loss of PEX26 contributes to the malignant phenotype of CRC. PEX26 may serve as a novel metastasis repressor for CRC.

Exosomal-miR-522-3p derived from cancer-associated fibroblasts accelerates tumor metastasis and angiogenesis via repression bone morphogenetic protein 5 in colorectal cancer

BACKGROUND

Colorectal cancer (CRC) is a gastrointestinal tract malignancy. Exosomes secreted by cancer-associated fibroblasts (CAFs) are reported to participate in tumor progression by delivering noncoding RNA or small proteins. However, the function of exosomal miR-522-3p in CRC remains unclear.

METHODS

CAFs were derived from tumor tissues, and exosomes were identified by western blot or TEM/NTA and originated from CAFs/NFs. The viability, invasion, and migration of HUVECs and CRC cells was examined using MTT, Transwell, and wound healing assays, respectively. The molecular interactions were validated using dual luciferase reporter assay and RIP. Xenograft and lung metastasis mouse models were generated to assess tumor growth and metastasis.

RESULTS

Exosomes extracted from CAFs/NFs showed high expression of CD63, CD81, and TSG101. CAF-derived exosomes significantly increased the viability, angiogenesis, invasion, and migration of HUVECs and CRC cells, thereby aggravating tumor growth, invasion, and angiogenesis in vivo. miR-522-3p was upregulated in CAF-derived exosomes and CRC tissues. Depletion of miR-522-3p reversed the effect of exosomes derived from CAFs in migration, angiogenesis, and invasion of HUVECs and CRC cells. Furthermore, bone morphogenetic protein 5 (BMP5) was identified as a target gene of miR-522-3p, and upregulation of BMP5 reversed the promoting effect of miR-522-3p mimics or CAF-derived exosomes on cell invasion, migration, and angiogenesis of HUVECs and CRC cells.

CONCLUSION

Exosomal miR-522-3p from CAFs promoted the growth and metastasis of CRC through downregulating BMP5, which might provide new strategies for the treatment of CRC.

Genomic alterations affecting tumor-infiltrating lymphocytes and PD-L1 expression patterns in triple-negative breast cancer

BACKGROUND

Tumor-infiltrating lymphocytes (TILs) and programmed cell death 1 ligand 1 (PD-L1) remain imperfect in predicting clinical outcomes of triple-negative breast cancer because outcomes do not always correlate with the expression of these biomarkers. Genomic and transcriptomic alterations that may contribute to the expression of these biomarkers remain incompletely uncovered.

METHODS

We evaluated PD-L1 immunohistochemistry scores (SP142 and 28-8 assays) and TILs in our triple-negative breast cancer multiomics dataset and 2 immunotherapy clinical trial cohorts. Then, we analyzed genomic and transcriptomic alterations correlated with TILs, PD-L1 expression, and patient outcomes.

RESULTS

Despite TILs serving as a decent predictor for triple-negative breast cancer clinical outcomes, exceptions remained. Our study revealed that several genomic alterations were correlated with unexpected events. In particular, PD-L1 expression may cause a paradoxical relationship between TILs and prognosis in certain patients. Consequently, we classified triple-negative breast cancers into 4 groups based on PD-L1 and TIL levels. The TIL-negative PD-L1-positive and TIL-positive PD-L1-negative groups were not typical "hot" tumors; both were associated with worse prognoses and lower immunotherapy efficacy than TIL-positive PD-L1-positive tumors. Copy number variation of PD-L1 and oncogenic signaling activation were correlated with PD-L1 expression in the TIL-negative PD-L1-positive group, whereas GSK3B-induced degradation may cause undetectable PD-L1 expression in the TIL-positive PD-L1-negative group. These factors have the potential to affect the predictive function of both PD-L1 and TILs.

CONCLUSIONS

Several genomic and transcriptomic alterations may cause paradoxical effects among TILs, PD-L1 expression, and prognosis in triple-negative breast cancer. Investigating and targeting these factors will advance precision immunotherapy for patients with this disease.

BMP signaling in cancer stemness and differentiation

The BMP (Bone morphogenetic protein) signaling pathway plays a central role in metazoan biology, intricately shaping embryonic development, maintaining tissue homeostasis, and influencing disease progression. In the context of cancer, BMP signaling exhibits context-dependent dynamics, spanning from tumor suppression to promotion. Cancer stem cells (CSCs), a modest subset of neoplastic cells with stem-like attributes, exert substantial influence by steering tumor growth, orchestrating therapy resistance, and contributing to relapse. A comprehensive grasp of the intricate interplay between CSCs and their microenvironment is pivotal for effective therapeutic strategies. Among the web of signaling pathways orchestrating cellular dynamics within CSCs, BMP signaling emerges as a vital conductor, overseeing CSC self-renewal, differentiation dynamics, and the intricate symphony within the tumor microenvironment. Moreover, BMP signaling's influence in cancer extends beyond CSCs, intricately regulating cellular migration, invasion, and metastasis. This multifaceted role underscores the imperative of comprehending BMP signaling's contributions to cancer, serving as the foundation for crafting precise therapies to navigate multifaceted challenges posed not only by CSCs but also by various dimensions of cancer progression. This article succinctly encapsulates the diverse roles of the BMP signaling pathway across different cancers, spanning glioblastoma multiforme (GBM), diffuse intrinsic pontine glioma (DIPG), colorectal cancer, acute myeloid leukemia (AML), lung cancer, prostate cancer, and osteosarcoma. It underscores the necessity of unraveling underlying mechanisms and molecular interactions. By delving into the intricate tapestry of BMP signaling's engagement in cancers, researchers pave the way for meticulously tailored therapies, adroitly leveraging its dualistic aspects-whether as a suppressor or promoter-to effectively counter the relentless march of tumor progression.

HEATR3 involved in the cell proliferation, metastasis and cell cycle development of bladder cancer acts as a tumor suppressor

The study was designed to detect the expression and clinical significance of the HEATR3 gene in bladder cancer (BCa) and to preliminarily explore whether this gene can affect the occurrence and development of BCa through the AKT/ERK signaling pathway. The expression and prognostic value of HEATR3 were explored based on The Cancer Genome Atlas (TCGA) and Genotypic Tissue Expression (GTEx) databases. Microarray immunohistochemical analysis was performed in 30 BCa cases to investigate the level of HEATR3 protein and to explore the relationship between HEATR3 and the clinicopathological features of BCa. Western Blot and qRT-PCR were used to detect HEATR3 protein and mRNA in BCa cell lines (5637, TCCSUP, SW780) and fallopian tube epithelial cell (SV-HUC-1). CCK8 method was employed to study the proliferation of BCa cells after heat treatment. Transwell assay was conducted to analyze the effect of HEATR3 on cell migration and invasion. And cell cycle and apoptosis were detected by flow cytometry. Furthermore, Western Blot assay was used to probe the effects of down-regulation of HEATR3 expression on the expression and phosphorylation levels of AKT and ERK proteins in BCa cells. Bioinformatics analysis showed that HEATR3 was significantly up-regulated in BCa, and high HEATR3 expression was associated with poor prognosis of BCa patients. In vitro experiments demonstrated that HEATR3 expression was up-regulated in BCa tissues compared with that in adjacent tissues. HEATR3 protein was also up-regulated in malignant cell lines. HEATR3 knockdown in BCa cells could inhibit cell proliferation, invasion and migration, block cell cycle and promote cell apoptosis. At the same time, HEATR3 knockdowns reduced the expression levels of p-AKT and p-ERK proteins. HEATR3 knockdown inhibits the development of BCa cells through the AKT/ERK signaling pathway. and it may become one of the most promising molecular targets for BCa treatment.

Loss of Tumor Suppressor C9orf9 Promotes Metastasis in Colorectal Cancer

The whole genome sequencing of tumor samples identifies thousands of somatic mutations. However, the function of these genes or mutations in regulating cancer progression remains unclear. We previously performed exome sequencing in patients with colorectal cancer, and identified one splicing mutation in C9orf9. The subsequent target sequencing of C9orf9 gene based on a validation cohort of 50 samples also found two function mutations, indicating that the loss of wild-type C9orf9 may participate in the tumorigenesis of colorectal cancer. In this research, we aimed to further confirm the function of C9orf9 in the CRC phenotype. Our Q-PCR analysis of the tumor and matched normal samples found that C9orf9 was downregulated in the CRC samples. Function assays revealed that C9orf9 exerts its tumor suppressor role mainly on cancer cell migration and invasion, and its loss was essential for certain tumor-microenvironment signals to induce EMT and metastasis in vivo. RNA-sequencing showed that stable-expressing C9orf9 can inhibit the expression of several metastasis-related genes and pathways, including vascular endothelial growth factor A (VEGFA), one of the essential endothelial cell mitogens which plays a critical role in normal physiological and tumor angiogenesis. Overall, our results showed that the loss of C9orf9 contributes to the malignant phenotype of CRC. C9orf9 may serve as a novel metastasis repressor for CRC.

DDI2 promotes tumor metastasis and resists antineoplastic drugs-induced apoptosis in colorectal cancer

The normal colorectal mucosa undergoes precancerous lesions that can develop over time into colorectal cancer (CRC). In the stage of precancerous lesions, DNA replication stress may lead to genome instability. We have performed whole-exome sequencing on genomic DNA obtained from three cases of CRC tissues and identified a novel frameshift mutation of DNA damage inducible 1 homolog 2 gene (DDI2, c. 854 del T). To date, there is no direct evidence that DDI2 is involved in the carcinogenesis of CRC. In this study, we demonstrated that DDI2 is upregulated in the early stage of CRC based on clinical samples and public databases. We also found that 5FU, a standard chemotherapeutic agent for CRC treatment, increased DDI2 mRNA levels in a dose-dependent manner. Depression of DDI2 inhibited CRC cell proliferation, migration and invasion both in vitro and in vivo. Transcriptome sequencing revealed that DDI2 was involved in the mitogen-activated protein kinase (MAPK) pathway. Furthermore, DDI2 resists a MAPK kinase (MEK) inhibitor (trametinib) and a PolyADP-ribose polymerase 1 (PARP1) inhibitor (talazoparib) induced apoptosis in CRC cells. Thus, our results indicate that DDI2 may play a vital role in the carcinogenesis of CRC and could serve as a promising therapeutic target for CRC.

Gene losses in the common vampire bat illuminate molecular adaptations to blood feeding

Vampire bats are the only mammals that feed exclusively on blood. To uncover genomic changes associated with this dietary adaptation, we generated a haplotype-resolved genome of the common vampire bat and screened 27 bat species for genes that were specifically lost in the vampire bat lineage. We found previously unknown gene losses that relate to reduced insulin secretion (FFAR1 and SLC30A8), limited glycogen stores (PPP1R3E), and a unique gastric physiology (CTSE). Other gene losses likely reflect the biased nutrient composition (ERN2 and CTRL) and distinct pathogen diversity of blood (RNASE7) and predict the complete lack of cone-based vision in these strictly nocturnal bats (PDE6H and PDE6C). Notably, REP15 loss likely helped vampire bats adapt to high dietary iron levels by enhancing iron excretion, and the loss of CYP39A1 could have contributed to their exceptional cognitive abilities. These findings enhance our understanding of vampire bat biology and the genomic underpinnings of adaptations to blood feeding.

Diagnostic and therapeutic biomarkers in colorectal cancer: a review

Colorectal cancer (CRC) is a public health concern and the second most common type of cancer among men and women causing a significant mortality. Biomarkers closely linked to the disease morbidity could holds potential as diagnostic and/or prognostic biomarker for the disease. This review provides an overview of recent advances in the search for colorectal cancer biomarkers through genomics and proteomics according to clinical function and application. Specifically, a number of biomarkers were identified and discussed. Emphasis was placed on their clinical applications relative to the diagnosis and prognosis of CRC. The discovery of more sensitive and specific markers for CRC is an urgent need, and the study of molecular targets is extremely important in this process, as they will allow for a better understanding of colorectal carcinogenesis, identification and validation of potential genetic signatures.

G9a Knockdown Suppresses Cancer Aggressiveness by Facilitating Smad Protein Phosphorylation through Increasing BMP5 Expression in Luminal A Type Breast Cancer

Epigenetic abnormalities affect tumor progression, as well as gene expression and function. Among the diverse epigenetic modulators, the histone methyltransferase G9a has been focused on due to its role in accelerating tumorigenesis and metastasis. Although epigenetic dysregulation is closely related to tumor progression, reports regarding the relationship between G9a and its possible downstream factors regulating breast tumor growth are scarce. Therefore, we aimed to verify the role of G9a and its presumable downstream regulators during malignant progression of breast cancer. G9a-depleted MCF7 and T47D breast cancer cells exhibited suppressed motility, including migration and invasion, and an improved response to ionizing radiation. To identify the possible key factors underlying these effects, microarray analysis was performed, and a TGF-β superfamily member, BMP5, was selected as a prominent target gene. It was found that BMP5 expression was markedly increased by G9a knockdown. Moreover, reduction in the migration/invasion ability of MCF7 and T47D breast cancer cells was induced by BMP5. Interestingly, a G9a-depletion-mediated increase in BMP5 expression induced the phosphorylation of Smad proteins, which are the intracellular signaling mediators of BMP5. Accordingly, we concluded that the observed antitumor effects may be based on the G9a-depletion-mediated increase in BMP5 expression and the consequent facilitation of Smad protein phosphorylation.

Upregulated bone morphogenetic protein 5 enhances proliferation and epithelial-mesenchymal transition process in benign prostatic hyperplasia via BMP/Smad signaling pathway

BACKGROUND

Benign prostatic hyperplasia (BPH) is one of the most common illnesses in aging men. Recent studies found that bone morphogenetic protein 5 (BMP5) is upregulated in BPH tissues, however, the role of BMP5 in the development of BPH has not been examined. The current study aims to elucidate the potential roles of BMP5 and related signaling pathways in BPH.

METHODS

Human prostate cell lines (BPH-1, WPMY-1) and human/rat hyperplastic prostate tissues were utilized. Western blot, quantitative real-time polymerase chain reaction, immunofluorescent staining, and immunohistochemical staining were performed. BMP5-silenced and -overexpressed cell models were generated and then cell cycle progression, apoptosis, and proliferation were determined. The epithelial-mesenchymal transition (EMT) was also quantitated. And rescue experiments by BMP/Smad signaling pathway agonist or antagonist were accomplished. Moreover, BPH-related tissue microarray analysis was performed and associations between clinical parameters and expression of BMP5 were analyzed.

RESULTS

Our study demonstrated that BMP5 was upregulated in human and rat hyperplastic tissues and localized both in the epithelial and stromal compartments of the prostate tissues. E-cadherin was downregulated in hyperplastic tissues, while N-cadherin and vimentin were upregulated. Overexpression of BMP5 enhanced cell proliferation and the EMT process via phosphorylation of Smad1/5/8, while knockdown of BMP5 induced cell cycle arrest at G0/G1 phase and blocked the EMT process. Moreover, a BMP/Smad signaling pathway agonist and antagonist reversed the effects of BMP5 silencing and overexpression, respectively. In addition, BMP5 expression positively correlated with prostate volume and total prostate-specific antigen.

CONCLUSION

Our novel data suggest that BMP5 modulated cell proliferation and the EMT process through the BMP/Smad signaling pathway which could contribute to the development of BPH. However, further studies are required to determine the exact mechanism. Our study also indicated that BMP/Smad signaling may be rediscovered as a promising new therapeutic target for the treatment of BPH.

The Prognostic Value of Bone Morphogenetic Proteins and Their Receptors in Lung Adenocarcinoma

Background

Bone morphogenetic proteins (BMPs) regulate tumor progression via binding to their receptors (BMPRs). However, the expression and clinical significance of BMPs/BMPRs in lung adenocarcinoma remain unclear due to a lack of systematic studies.

Methods

This study screened differentially expressed BMPs/BMPRs (deBMPs/BMPRs) in a training dataset combining TCGA-LUAD and GTEx-LUNG and verified them in four GEO datasets. Their prognostic value was evaluated via univariate and multivariate Cox regression analyses. LASSO was performed to construct an initial risk model. Subsequently, after weighted gene co-expression network analysis (WGCNA), differential expression analysis, and univariate Cox regression analysis, hub genes co-expressed with differentially expressed BMPs/BMPRs were filtered out to improve the risk model and explore potential mechanisms. The improved risk model was re-established via LASSO combining hub genes with differentially expressed BMPs/BMPRs as the core. In the testing cohort including 93 lung adenocarcinoma patients, immunohistochemistry (IHC) was performed to verify BMP5 protein expression and its association with prognosis.

Results

BMP2, BMP5, BMP6, GDF10, and ACVRL1 were verified as downregulated in lung adenocarcinoma. Survival analysis identified BMP5 as an independent protective prognostic factor. We also found that BMP5 was significantly correlated with EGFR expression and mutations, suggesting that BMP5 may play a role in targeted therapy. The initial risk model containing only BMP5 showed a significant correlation (HR: 1.71, 95% CI: 1.28−2.28, p: 3e-04) but low prognostic accuracy (AUC of 1-year survival: 0.6, 3-year survival: 0.6, 5-year survival: 0.63). Seventy-nine hub genes co-expressed with BMP5 were identified, and their functions were enriched in cell migration and tumor metastasis. The re-established risk model showed greater prognostic correlation (HR: 2.58, 95% CI: 1.92–3.46, p: 0) and value (AUC of 1-year survival: 0.72, 3-year survival: 0.69, and 5-year survival: 0.68). IHC results revealed that BMP5 protein was also downregulated in lung adenocarcinoma and higher expression was markedly associated with better prognosis (HR: 0.44, 95% CI: 0.23–0.85, p: 0.0145).

Conclusion

BMP5 is a potential crucial target for lung adenocarcinoma treatment based on significant differential expression and superior prognostic value.

Transcriptomic Characterization of Postmolar Gestational Choriocarcinoma

The human placenta shares properties with solid tumors, such as rapid growth, tissue invasion, cell migration, angiogenesis, and immune evasion. However, the mechanisms that drive the evolution from premalignant proliferative placental diseases—called hydatidiform moles—to their malignant counterparts, gestational choriocarcinoma, as well as the factors underlying the increased aggressiveness of choriocarcinoma arising after term delivery compared to those developing from hydatidiform moles, are unknown. Using a 730-gene panel covering 13 cancer-associated canonical pathways, we compared the transcriptomic profiles of complete moles to those of postmolar choriocarcinoma samples and those of postmolar to post-term delivery choriocarcinoma. We identified 33 genes differentially expressed between complete moles and postmolar choriocarcinoma, which revealed TGF-β pathway dysregulation. We found the strong expression of SALL4, an upstream regulator of TGF-β, in postmolar choriocarcinoma, compared to moles, in which its expression was almost null. Finally, there were no differentially expressed genes between postmolar and post-term delivery choriocarcinoma samples. To conclude, the TGF-β pathway appears to be a crucial step in the progression of placental malignancies. Further studies should investigate the value of TGF- β family members as biomarkers and new therapeutic targets.

Hsa_circ_0069244 acts as the sponge of miR-346 to inhibit non-small cell lung cancer progression by regulating XPC expression

Circular RNAs (circRNAs) play a significant role in the progression of diverse malignancies. Here, we aimed to probe the function and mechanism of circ_0069244 in non-small cell lung cancer (NSCLC). In the present study, circ_0069244 was selected from the circRNA microarray datasets (GSE112214). Quantitative real-time polymerase chain reaction (qRT-PCR) was applied to examine circ_0069244, miR-346 and XPC complex subunit, DNA damage recognition and repair factor (XPC) expression levels. Kaplan-Meier curve was employed to analyze the association between circ_0069244 expression and overall survival of NSCLC patients. Cell counting kit-8 (CCK-8) and 5-Bromo-2'-deoxyuridine (BrdU) experiments were utilized to examine the proliferation of NSCLC cells. Scratch healing and Transwell experiments were executed to examine the migration of NSCLC cells. Western blot was conducted to detect XPC expression at protein level in NSCLC cells. Bioinformatics analysis, dual-luciferase reporter gene and RNA immunoprecipitation (RIP) experiments predicted and validated the targeting relationships of circ_0069244 and miR-346, as well as miR-346 and 3'untranslated region (UTR) of XPC mRNA, respectively. We reported that circ_0069244 was remarkably down modulated in NSCLC and was linked to shorter survival and poor tumor histological grade in NSCLC patients. Functionally, circ_0069244 repressed NSCLC cell proliferation and migration. Furthermore, miR-346-5p was unveiled to be a downstream target of circ_0069244, and miR-346-5p specifically modulated XPC expression. Rescue experiments indicated that the inhibitory effect of circ_0069244 was abolished by co-expression of miR-346-5p mimics. Taken together, circ_0069244 restrained NSCLC progression by modulating the miR-346-5p/XPC axis.

Epithelial-Mesenchymal Transition Gene Signature Related to Prognostic in Colon Adenocarcinoma

Colon adenocarcinoma (COAD) remains an important cause of cancer-related mortality worldwide. Epithelial–mesenchymal transition (EMT) is a key mechanism, promoting not only the invasive or metastatic phenotype but also resistance to therapy. Using bioinformatics approaches, we studied the alteration on EMT related genes and its implication on COAD prognostic based on public datasets. For the EMT mechanisms, two overexpressed genes were identified (NOX4 and IGF2BP3), as well as five downregulated genes (BMP5, DACT3, EEF1A2, GCNT2 and SFRP1) that were related to prognosis in COAD. A qRT-PCR validation step was conducted in a COAD patient cohort comprising of 29 tumor tissues and 29 normal adjacent tissues, endorsing the expression level for BMP5, as well as for two of the miRNAs targeting key EMT related genes, revealing upregulation of miR-27a-5p and miR-146a-5p. The EMT signature can be used to develop a panel of biomarkers for recurrence prediction in COAD patients, which may contribute to the improvement of risk stratification for the patients.

Familial Occurrence of Adult Granulosa Cell Tumors: Analysis of Whole-Genome Germline Variants

Simple Summary

Although granulosa cell tumors can occur in rare syndromes and one familial case of a granulosa cell tumor has been described, a genetic predisposition for granulosa cell tumors has not been identified. Through our collaborations with patients, we identified four families in which two women of each family were diagnosed with an adult granulosa cell tumor. Although predicted deleterious variants in PIK3C2G, BMP5, and LRP2 were found, we could not identify an overlapping genetic variant or affected locus that may explain a genetic predisposition for granulosa cell tumors. The age of onset in the familial patients was significantly lower (median 38 years, range from 17 to 60) than in sporadic patients (median between 50 and 55 years). Furthermore, breast cancer, polycystic ovary syndrome, and subfertility were seen in these families.

Abstract

Adult granulosa cell tumor (AGCT) is a rare ovarian cancer subtype, with a peak incidence around 50–55 years. Although AGCT can occur in specific syndromes, a genetic predisposition for AGCT has not been identified. The aim of this study is to identify a genetic variant in families with AGCT patients, potentially contributing to tumor evolution. We identified four families, each including two women diagnosed with AGCT. Whole-genome sequencing was performed to identify overlapping germline variants or affected genes. Familial relationship was evaluated using genealogy and genomic analyses. Patient characteristics, medical (family) history, and pedigrees were collected. Findings were compared to a reference group of 33 unrelated AGCT patients. Mean age at diagnosis was 38 years (range from 17 to 60) versus 51 years in the reference group, and seven of eight patients were premenopausal. In two families, three first degree relatives were diagnosed with breast cancer. Furthermore, polycystic ovary syndrome (PCOS) and subfertility was reported in three families. Predicted deleterious variants in PIK3C2G, BMP5, and LRP2 were identified. In conclusion, AGCTs occur in families and could potentially be hereditary. In these families, the age of AGCT diagnosis is lower and cases of breast cancer, PCOS, and subfertility are present. We could not identify an overlapping genetic variant or affected locus that may explain a genetic predisposition for AGCT.

Decreased level of miR-1301 promotes colorectal cancer progression via activation of STAT3 pathway

The molecular pathogenesis of colorectal cancer (CRC) has been widely investigated in recent years. Accumulating evidence has indicated that microRNA (miRNA) dysregulation participates in the processes of driving CRC initiation and progression. Aberrant expression of miR-1301 has been found in various tumor types. However, its role in CRC remains to be elucidated. In the present study, we identified miR-1301 was enriched in normal colorectal tissues and significantly down-regulated in CRC. Decreased level of miR-1301 strongly correlated with aggressive pathological characteristics, including advanced stage and metastasis. Bioinformatics and dual luciferase assay demonstrated that STAT3 is a direct target of miR-1301. Gain and loss-of-function assays showed that miR-1301 had no effect on cell proliferation. Overexpression of miR-1301 suppressed cell migration and invasion capacity of pSTAT3-positive LoVo cells, but not pSTAT3-negative SW480 cells, while inhibition of miR-1301 consistently promoted cell migration and invasion in both cell lines. Additionally, miR-1301 inhibition restored the suppressed migration and invasion of STAT3-knockdown LoVo cells. MiR-1301 functioned as a tumor suppressor to modulate the IL6/STAT3 signaling pathway. In summary, this study highlights the significant role of miR-1301/STAT3 axis in CRC metastasis.

Paving the way for small-molecule drug discovery

Small-molecule drugs are organic compounds affecting molecular pathways by targeting important proteins, which have a low molecular weight, making them penetrate cells easily. Small-molecule drugs can be developed from leads derived from rational drug design or isolated from natural resources. As commonly used medications, small-molecule drugs can be taken orally, which enter cells to act on intracellular targets. These characteristics make small-molecule drugs promising candidates for drug development, and they are increasingly favored in the pharmaceutical market. Despite the advancements in molecular genetics and effective new processes in drug development, the drugs currently used in clinical practice are inadequate due to their poor efficacy or severe side effects. Therefore, developing new safe and efficient drugs is a top priority for disease control and curing.

Prognostic Signatures Based on Thirteen Immune-Related Genes in Colorectal Cancer

Background

The immunosuppressive microenvironment is closely related to tumorigenesis and cancer development, including colorectal cancer (CRC). The aim of the current study was to identify new immune biomarkers for the diagnosis and treatment of CRC.

Materials and Methods

CRC data were downloaded from the Gene Expression Omnibus and The Cancer Genome Atlas databases. Sequences of immune-related genes (IRGs) were obtained from the ImmPort and InnateDB databases. Gene set enrichment analysis (GSEA) and transcription factor regulation analysis were used to explore potential mechanisms. An immune-related classifier for CRC prognosis was conducted using weighted gene co-expression network analysis (WGCNA), Cox regression analysis, and least absolute shrinkage and selection operator (LASSO) analysis. ESTIMATE and CIBERSORT algorithms were used to explore the tumor microenvironment and immune infiltration in the high-risk CRC group and the low-risk CRC group.

Results

By analyzing the IRGs that were significantly associated with CRC in the module, a set of 13 genes (CXCL1, F2RL1, LTB4R, GPR44, ANGPTL5, BMP5, RETNLB, MC1R, PPARGC1A, PRKDC, CEBPB, SYP, and GAB1) related to the prognosis of CRC were identified. An IRG-based prognostic signature that can be used as an independent potentially prognostic indicator was generated. The ROC curve analysis showed acceptable discrimination with AUCs of 0.68, 0.68, and 0.74 at 1-, 3-, and 5- year follow-up respectively. The predictive performance was validated in the train set. The potential mechanisms and functions of prognostic IRGs were analyzed, i.e., NOD-like receptor signaling, and transforming growth factor beta (TGFβ) signaling. Besides, the stromal score and immune score were significantly different in high-risk group and low-risk group (p=4.6982e-07, p=0.0107). Besides, the proportions of resting memory CD4⁺ T cells was significantly higher in the high-risk groups.

Conclusions

The IRG-based classifier exhibited strong predictive capacity with regard to CRC. The survival difference between the high-risk and low-risk groups was associated with tumor microenvironment and immune infiltration of CRC. Innovative biomarkers for the prediction of CRC prognosis and response to immunological therapy were identified in the present study.

Systematic analysis using a bioinformatics strategy identifies SFTA1P and LINC00519 as potential prognostic biomarkers for lung squamous cell carcinoma

Lung cancer has high incidence and mortality rates, in which lung squamous cell carcinoma (LUSC) is a primary type of non-small cell lung carcinoma (NSCLC). The aim of our study was to discover long non-coding RNAs (lncRNAs) associated with diagnose and prognosis for LUSC. RNA sequencing data obtained from LUSC samples were extracted from The Cancer Genome Atlas database (TCGA). Two prognosis-associated lncRNAs (including SFTA1P and LINC00519) were selected from LUSC samples, and the expression levels were also verified to be associated abnormal in LUSC clinical samples. Our findings demonstrate that lncRNAs SFTA1P and LINC00519 exert important functions in human LUSC and may serve as new targets for LUSC diagnosis and therapy.

C8orf48 inhibits the tumorigenesis of colorectal cancer by regulating the MAPK signaling pathway

AIMS

Colorectal cancer (CRC) is a leading cause of cancer-related death globally. Thus, in this study, we aimed to investigate chromosome 8 open reading frame 48 (C8orf48) as a biomarker for early detection of CRC.

MAIN METHODS

RNA expression and methylation profiles were downloaded from The Cancer Genome Atlas (TCGA) database. Cell proliferation, migration and invasion assays were performed to confirm the function of C8orf48 in CRC cells. Dual-luciferase reporter assay was used to identify that C8orf48 was the direct target of miR-556. Genomics of Drug Sensitivity in Cancer (GDSC) database, gene set enrichment analysis (GSEA) and western blot analysis were performed to explore the mechanism of C8orf48.

KEY FINDINGS

we found that C8orf48 is down-regulated in clinical samples of CRC tissues. Enrichment analysis showed that C8orf48 is associated with methylation biomarkers in CRC, and TCGA database confirmed that the methylation of C8orf48 is up-regulated in the early stage of CRC. We further revealed that the overexpression of C8orf48 decreased CRC cell proliferation, migration and invasion. Luciferase reporter indicated that C8orf48 was the direct target of the oncogene miR-556. Additionally, we used GDSC database, GSEA database and western blot analysis to demonstrate that C8orf48 plays a suppressor role in CRC by inhibiting MAPK signaling pathway.

SIGNIFICANCE

C8orf48 was identified as a biomarker for early detection of CRC for the first time, and might provide novel information for CRC prediction and therapy.

Pathologic evolution-related Gene Analysis based on both single-cell and bulk transcriptomics in Colorectal Cancer

Purpose: The patients diagnosed with colorectal cancer (CRC) are likely to undergo differential outcomes in clinical survival owing to different pathologic stages. However, signatures in association with pathologic evolution and CRC prognosis are not clearly defined. This study aimed to identify pathologic evolution-related genes in CRC based on both single-cell and bulk transcriptomics.

Patients and methods: The CRC single-cell transcriptomic dataset (GSE81861, n=590) with clinical information and tumor microenvironmental tissues was collected to identify the pathologic evolution-related genes. The colonic adenocarcinoma and rectum adenocarcinoma transcriptomics from The Cancer Genome Atlas were obtained as the training dataset (n=363) and 5 other CRC transcriptomics cohorts from Gene Expression Omnibus (n=1031) were acquired as validation data. Graph-based clustering analysis algorithm was applied to identify pathologic evolution-related cell populations. Pseudotime analysis was performed to construct the trajectory plot of pathologic evolution and to define hub genes in the evolution process. Cell-type identification by estimating relative subsets of RNA transcripts was then executed to build a novel cell infiltration classifier. The prediction efficacy of this classifier was validated in bulk transcriptomic datasets.

Results: Epithelial and T cells were elucidated to be related to the pathologic stages in CRC tissues. Pseudotime analysis and survival analysis indicated that HOXC5, HOXC8 and BMP5 were the marker genes in pathologic evolution process. Our cell infiltration classifier exhibited excellent forecast efficacy in predicting pathologic stages and prognosis of CRC patients.

Conclusion: We identified pathologic evolution-related genes in single-cell transcriptomic and proposed a novel specific cell infiltration classifier to forecast the prognosis of CRC patients based on pathologic stage-related hub genes HOXC6, HOXC8 and BMP5.

Uncovering epigenetic landscape: a new path for biomarkers identification and drug development

Scientific advances in recent decades have revealed an incredible degree of plasticity in gene expression in response to various environmental, nutritional, physiological, pathological, and behavioral conditions. Epigenetics emerges in this sense, as the link between the internal (genetic) and external (environmental) factors underlying the expression of the phenotype. Methylation of DNA and histone post-translationa modifications are canonical epigenetic events. Additionally, noncoding RNAs molecules (microRNAs and lncRNAs) have also been proposed as another layer of epigenetic regulation. Together, these events are responsible for regulating gene expression throughout life, controlling cellular fate in both normal and pathological development. Despite being a relatively recent science, epigenetics has been arousing the interest of researchers from different segments of the life sciences and the general public. This review highlights the recent advances in the characterization of the epigenetic events and points promising use of these brands for the diagnosis, prognosis, and therapy of diseases. We also present several classes of epigenetic modifying compounds with therapeutic applications (so-call epidrugs) and their current status in clinical trials and approved by the FDA. In summary, hopefully, we provide the reader with theoretical bases for a better understanding of the epigenetic mechanisms and of the promising application of these marks and events in the medical clinic.

Detection of circulating BMP5 as a risk factor for Barrett's esophagus

Barrett's esophagus (BE) predisposes for the malignant condition of esophageal adenocarcinoma (EAC). Since BE patients have few or no symptoms, most of these patients are not identified and not included in surveillance programs. These BE patients are at risk of developing advanced-stage EAC. At present, non-invasive tests to identify BE patients from the general population are lacking. We and others showed that Bone Morphogenetic Protein 4 (BMP4), and other BMPs are upregulated in BE. We aimed to determine if circulating BMPs can be identified and used as blood biomarkers to identify BE patients at high risk in the general population. In this study, we could detect the different BMPs in the blood of 112 BE patients and 134 age- and sex-matched controls. Concentration levels of BMP2, BMP4, and BMP5 were elevated in BE patients, with BMP2 and BMP5 significantly increased. BMP5 remained significant after multivariate analysis and was associated with an increased risk for BE with an OR of 1.49 (p value 0.01). Per log (pg/mL) of BMP5, the odds of having BE increased by 50%. Future optimization and validation studies might be needed to prove its utility as a non-invasive method for the detection of BE in high-risk populations and screening programs.

MiR-581/SMAD7 Axis Contributes to Colorectal Cancer Metastasis: A Bioinformatic and Experimental Validation-Based Study

Metastasis is a well-known poor prognostic factor and primary cause of mortality in patients with colorectal cancer (CRC). Recently, with the progress of high through-put sequencing, aberrantly expressed non-coding RNAs (ncRNAs) were found to participate in the initiation and development of cancer. However, the mechanisms of ncRNA-mediated regulation of metastasis in CRC remain largely unknown. In this study, we systematically analyzed the expression network of microRNAs (miRNAs) and genes in CRC metastasis using bioinformatics, and discovered that the miR-581/SMAD7 axis could be a potential factor that drives CRC metastasis. A dual luciferase report assay and protein analysis confirmed the binding relationship between miR-581 and SMAD7. Further functional assays revealed that miR-581 inhibition could suppress cell proliferation and induce apoptosis in SW480 cells. Up-regulation or down-regulation of miR-581 could both affect cell invasion capacity and modulate epithelial to mesenchymal transition (EMT) via a SMAD7/TGFβ signaling pathway. In conclusion, our findings elucidated that miR-581/SMAD7 could be essential for CRC metastasis, and may serve as a potential therapeutic target for CRC patients.

Targeting Colorectal Cancer Stem Cells as an Effective Treatment for Colorectal Cancer

As one of the common cancers that threaten human life, the recurrence and metastasis of colorectal cancer seriously affect the prognosis of patients. Although new drugs and comprehensive treatments have been adopted, the current treatment effect on this tumor, especially in advanced colorectal cancer, is still not satisfactory. More and more evidence shows that tumors are likely to be a stem cell disease. In recent years, the rise of cancer stem cell theory has provided a new way for cancer treatment. Studies have found that a small number of special cells in colorectal cancer tissues that induce tumorigenesis, proliferation, and promote tumor migration and metastasis, namely, colorectal cancer stem cells. Colorectal cancer stem cells are defined with a group of cell-surface markers, such as CD44, CD133, CD24, epithelial cell adhesion factor molecule, LGR5, and acetaldehyde dehydrogenase. They are highly tumorigenic, aggressive, and chemoresistant and thus are critical in the metastasis and recurrence of colorectal cancer. Therefore, targeting colorectal cancer stem cells may become an important research direction for the future cure of colorectal cancer.

Colorectal Cancer and Bone Tissue: Fantastic Relations and Where to Find Them

Colorectal cancer (CRC) is the third most common cancer worldwide. There is a need for the early diagnosis of CRC for a better prognostic outcome. It is, therefore, crucial to understand the CRC pathogenesis in all its aspects. In many cases, one of the main causes of cancer-related deaths is the presence of metastases. In this context, an often overlooked aspect is the metastatic tropism, since CRC, like other cancers, is more prone to metastasize some organs rather than others. Beyond the liver and lung, and differently from other types of cancers, a not usual site of CRC metastases is the bone. However, it may assume a crucial role in the development and the outcome of the disease. Therefore, this review aims to discuss the complex relations between bone markers and CRC pathogenesis, suggesting the use of these molecules as potential targets for therapeutic purposes. Different osteogenic molecules, some of whom are growth factors and are implicated in the different osteogenic pathways, have been proved to also be involved in CRC progression. Some of them are oncogenes, while others oncosuppressors, and in a future perspective, some of them may represent new potential CRC biomarkers.

Intracellular Fate and Impact on Gene Expression of Doxorubicin/Cyclodextrin-Graphene Nanomaterials at Sub-Toxic Concentration

The graphene road in nanomedicine still seems very long and winding because the current knowledge about graphene/cell interactions and the safety issues are not yet sufficiently clarified. Specifically, the impact of graphene exposure on gene expression is a largely unexplored concern. Herein, we investigated the intracellular fate of graphene (G) decorated with cyclodextrins (CD) and loaded with doxorubicin (DOX) and the modulation of genes involved in cancer-associated canonical pathways. Intracellular fate of GCD@DOX, tracked by FLIM, Raman mapping and fluorescence microscopy, evidenced the efficient cellular uptake of GCD@DOX and the presence of DOX in the nucleus, without graphene carrier. The NanoString nCounter™ platform provided evidence for 34 (out of 700) differentially expressed cancer-related genes in HEp-2 cells treated with GCD@DOX (25 µg/mL) compared with untreated cells. Cells treated with GCD alone (25 µg/mL) showed modification for 16 genes. Overall, 14 common genes were differentially expressed in both GCD and GCD@DOX treated cells and 4 of these genes with an opposite trend. The modification of cancer related genes also at sub-cytotoxic G concentration should be taken in consideration for the rational design of safe and effective G-based drug/gene delivery systems. The reliable advantages provided by NanoString^® technology, such as sensibility and the direct RNA measurements, could be the cornerstone in this field.

contamDE-lm: linear model-based differential gene expression analysis using next-generation RNA-seq data from contaminated tumor samples

MOTIVATION

Tumor and adjacent normal RNA samples are commonly used to screen differentially expressed genes between normal and tumor samples or among tumor subtypes. Such paired-sample design could avoid numerous confounders in differential expression (DE) analysis, but the cellular contamination of tumor samples can be an important noise and confounding factor, which can both inflate false-positive rate and deflate true-positive rate. The existing DE tools that use next-generation RNA-seq data either do not account for cellular contamination or are computationally extensive with increasingly large sample size.

RESULTS

A novel linear model was proposed to avoid the problem that could arise from tumor-normal correlation for paired samples. A statistically robust and computationally very fast DE analysis procedure, contamDE-lm, was developed based on the novel model to account for cellular contamination, boosting DE analysis power through the reduction in individual residual variances using gene-wise information. The desired advantages of contamDE-lm over some state-of-the-art methods (limma and DESeq2) were evaluated through the applications to simulated data, TCGA database and Gene Expression Omnibus (GEO) database.

AVAILABILITY AND IMPLEMENTATION

The proposed method contamDE-lm was implemented in an updated R package contamDE (version 2.0), which is freely available at https://github.com/zhanghfd/contamDE.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Utilizing Experimental Mouse Model to Identify Effectors of Hepatocellular Carcinoma Induced by HBx Antigen

Hepatocellular carcinoma (HCC) is among the ten most commonly diagnosed cancers and the fourth leading cause of cancer-related death. Patients with hepatitis B virus (HBV) infection are prone to developing chronic liver diseases (i.e., fibrosis and cirrhosis), and the HBV X antigen plays an important role in the development of HCC. The difficulty in detecting HCC at the early stages is one of the main reasons that the death rate approximates the incidence rate. The regulators controlling the downstream liver protein expression from HBV infection are unclear. Mass spectrometric techniques and customized programs were used to identify differentially expressed proteins which may be involved in the development of liver fibrosis and HCC progression in hepatitis B virus X protein transgenic mice (HBx mice). FSTL1, CTSB, and TGF-β enhanced the signaling pathway proteins during the pathogenesis of HBx. Missing proteins can be essential in cell growth, differentiation, apoptosis, migration, metastasis or angiogenesis. We found that LHX2, BMP-5 and GDF11 had complex interactions with other missing proteins and BMP-5 had both tumor suppressing and tumorigenic roles. BMP-5 may be involved in fibrosis and tumorigenic processes in the liver. These results provide us an understanding of the mechanism of HBx-induced disorders, and may serve as molecular targets for liver treatment.

Construction of a novel gene-based model for prognosis prediction of clear cell renal cell carcinoma

BACKGROUND

Clear cell renal cell carcinoma (ccRCC) comprises the majority of kidney cancer death worldwide, whose incidence and mortality are not promising. Identifying ideal biomarkers to construct a more accurate prognostic model than conventional clinical parameters is crucial.

METHODS

Raw count of RNA-sequencing data and clinicopathological data were acquired from The Cancer Genome Atlas (TCGA). Tumor samples were divided into two sets. Differentially expressed genes (DEGs) were screened in the whole set and prognosis-related genes were identified from the training set. Their common genes were used in LASSO and best subset regression which were performed to identify the best prognostic 5 genes. The gene-based risk score was developed based on the Cox coefficient of the individual gene. Time-dependent receiver operating characteristic (ROC) and Kaplan-Meier (KM) survival analysis were used to assess its prognostic power. GSE29609 dataset from GEO (Gene Expression Omnibus) database was used to validate the signature. Univariate and multivariate Cox regression were performed to screen independent prognostic parameters to construct a nomogram. The predictive power of the nomogram was revealed by time-dependent ROC curves and the calibration plot and verified in the validation set. Finally, Functional enrichment analysis of DEGs and 5 novel genes were performed to suggest the potential biological pathways.

RESULTS

PADI1, ATP6V0D2, DPP6, C9orf135 and PLG were screened to be significantly related to the prognosis of ccRCC patients. The risk score effectively stratified the patients into high-risk group with poor overall survival (OS) based on survival analysis. AJCC-stage, age, recurrence and risk score were regarded as independent prognostic parameters by Cox regression analysis and were used to construct a nomogram. Time-dependent ROC curves showed the nomogram performed best in 1-, 3- and 5-year survival predictions compared with AJCC-stage and risk score in validation sets. The calibration plot showed good agreement of the nomogram between predicted and observed outcomes. Functional enrichment analysis suggested several enriched biological pathways related to cancer.

CONCLUSIONS

In our study, we constructed a gene-based model integrating clinical prognostic parameters to predict prognosis of ccRCC well, which might provide a reliable prognosis assessment tool for clinician and aid treatment decision-making in the clinic.

A Multi-Omics Analysis of Bone Morphogenetic Protein 5 (BMP5) mRNA Expression and Clinical Prognostic Outcomes in Different Cancers Using Bioinformatics Approaches

Cumulative studies have provided controversial evidence for the prognostic values of bone morphogenetic protein 5 (BMP5) in different types of cancers such as colon, breast, lung, bladder, and ovarian cancer. To address the inconsistent correlation of BMP5 expression with patient survival and molecular function of BMP5 in relation to cancer progression, we performed a systematic study to determine whether BMP5 could be used as a prognostic marker in human cancers. BMP5 expression and prognostic values were assessed using different bioinformatics tools such as ONCOMINE, GENT, TCGA, GEPIA, UALCAN, PrognoScan, PROGgene V2 server, and Kaplan-Meier Plotter. In addition, we used cBioPortal database for the identification and analysis of BMP5 mutations, copy number alterations, altered expression, and protein-protein interaction (PPI). We found that BMP5 is frequently down-regulated in our queried cancer types. Use of prognostic analysis showed negative association of BMP5 down-regulation with four types of cancer except for ovarian cancer. The highest mutation was found in the R321*/Q amino acid of BMP5 corresponding to colorectal and breast cancer whereas the alteration frequency was higher in lung squamous carcinoma datasets (>4%). In PPI analysis, we found 31 protein partners of BMP5, among which 11 showed significant co-expression (p-value < 0.001, log odds ratio > 1). Pathway analysis of differentially co-expressed genes with BMP5 in breast, lung, colon, bladder and ovarian cancers revealed the BMP5-correlated pathways. Collectively, this data-driven study demonstrates the correlation of BMP5 expression with patient survival and identifies the involvement of BMP5 pathways that may serve as targets of a novel biomarker for various types of cancers in human.

Construction and Validation of an Immune-Related Prognostic Model Based on TP53 Status in Colorectal Cancer

Growing evidence has indicated that prognostic biomarkers have a pivotal role in tumor and immunity biological processes. TP53 mutation can cause a range of changes in immune response, progression, and prognosis of colorectal cancer (CRC). Thus, we aim to build an immunoscore prognostic model that may enhance the prognosis of CRC from an immunological perspective. We estimated the proportion of immune cells in the GSE39582 public dataset using the CIBERSORT (Cell type identification by estimating relative subset of known RNA transcripts) algorithm. Prognostic genes that were used to establish the immunoscore model were generated by the LASSO (Least absolute shrinkage and selection operator) Cox regression model. We established and validated the immunoscore model in GEO (Gene Expression Omnibus) and TCGA (The Cancer Genome Atlas) cohorts, respectively; significant differences of overall survival analysis were found between the low and high immunoscore groups or TP53 subgroups. In the multivariable Cox analysis, we observed that the immunoscore was an independent prognostic factor both in the GEO cohort (HR (Hazard ratio) 1.76, 95% CI (confidence intervals): 1.26-2.46) and the TCGA cohort (HR 1.95, 95% CI: 1.20-3.18). Furthermore, we established a nomogram for clinical application, and the results suggest that the nomogram is a better predictive model for prognosis than immunoscore or TNM staging.

MiR-496 promotes migration and epithelial-mesenchymal transition by targeting RASSF6 in colorectal cancer

Aberrant loss of tumor-suppressor genes plays a crucial role in tumorigenesis and development of colorectal cancer (CRC). Extensive studies have reported tha hypermethylation of Ras association domain family member 6 (RASSF6) is common in various solid tumors. Another important mode of epigenetic regulation, microRNA (miRNA) regulation of RASSF6, is far from clear. The aim of the present work was to screen out novel miRNA regulating RASSF6, and to explore its underlying mechanism in CRC. With the use of bioinformatics, clinical sample data, and luciferase binding assay, we determined that microRNA-496 (miR-496) could be a novel oncomiR that directly binds to RASSF6. Next, a series of miR-496 mimics or inhibitor, or RASSF6 small interfering RNA (siRNA) introduced into CRC cells were applied to examine the effect of miR-496 on CRC cell viability, migration, and epithelial-mesenchymal transition (EMT). The results demonstrated that miR-496/RASSF6 could promote cell migration and EMT via Wnt signaling activation, but had no effect on cell viability. Our results confirmed that the miR-496/RASSF6 axis is involved in Wnt pathway-mediated tumor metastasis, highlighting its potential as a therapeutic target for CRC.

The Akt/FoxO/p27Kip1 axis contributes to the anti-proliferation of pentoxifylline in hypertrophic scars

Hypertrophic scars (HS) are characterized by the excessive production and deposition of extracellular matrix (ECM) proteins. Pentoxifylline (PTX), a xanthine derived antioxidant, inhibits the proliferation, inflammation and ECM accumulation of HS. In this study, we aimed to explore the effect of PTX on HS and further clarify the mechanism of PTX‐induced anti‐proliferation. We found that PTX could significantly attenuate proliferation of HS fibroblasts and fibrosis in an animal HS model. PTX inhibited the proliferation of HSFs in a dose‐ and time‐dependent manner, and this growth inhibition was mainly mediated by cell cycle arrest. Transcriptome sequencing showed that PTX affects HS formation through the PI3K/Akt/FoxO1 signalling pathway to activate p27^Kip1. PTX down‐regulated p‐Akt and up‐regulated p‐FoxO1 in TGF‐β1 stimulated fibroblasts at the protein level, and simultaneously, the expression of p27^Kip1 was activated. In a mouse model of HS, PTX treatment resulted in the ordering of collagen fibres. The results revealed that PTX regulates TGFβ1‐induced fibroblast activation and inhibits excessive scar formation. Therefore, PTX is a promising agent for the treatment of HS formation.

Alteration of the tumor suppressor SARDH in sporadic colorectal cancer: A functional and transcriptome profiling-based study

Sporadic colorectal cancer (sCRC) is one of the leading causes of cancer death worldwide. As a highly heterogeneous complex disease, the currently reported classical genetic markers for sCRC, including APC, KRAS, BRAF, and TP53 gene mutations and epigenetic alterations, can explain only some sCRC patients. Here, we first reported a deleterious c.551C>T mutation in SARDH in sCRC. SARDH was identified as a novel tumor suppressor gene and was abnormally decreased in sCRC at both the transcriptional and the translational level. SARDH mRNA levels were also down-regulated in oesophageal cancer, lung cancer, liver cancer, and pancreatic cancer in the TCGA database. SARDH overexpression inhibited the proliferation, migration, and invasion of CRC cell lines, whereas its depletion improved these processes. SARDH overexpression was down-regulated in multiple pathways, especially in the chemokine pathway. The SARDH transcript level was positively correlated with the methylation states of CXCL1 and CCL20. Therefore, we concluded that SARDH depletion is involved in the development of sCRC.