Analyzing the disease module associated with osteosarcoma via a network- and pathway-based approach

Osteosarcoma-enriched transcripts paradoxically generate osteosarcoma-suppressing extracellular proteins

Osteosarcoma (OS) is the common primary bone cancer that affects mostly children and young adults. To augment the standard-of-care chemotherapy, we examined the possibility of protein-based therapy using mesenchymal stem cells (MSCs)-derived proteomes and OS-elevated proteins. While a conditioned medium (CM), collected from MSCs, did not present tumor-suppressing ability, the activation of PKA converted MSCs into induced tumor-suppressing cells (iTSCs). In a mouse model, the direct and hydrogel-assisted administration of CM inhibited tumor-induced bone destruction, and its effect was additive with cisplatin. CM was enriched with proteins such as calreticulin, which acted as an extracellular tumor suppressor by interacting with CD47. Notably, the level of CALR transcripts was elevated in OS tissues, together with other tumor-suppressing proteins, including histone H4, and PCOLCE. PCOLCE acted as an extracellular tumor-suppressing protein by interacting with amyloid precursor protein, a prognostic OS marker with poor survival. The results supported the possibility of employing a paradoxical strategy of utilizing OS transcriptomes for the treatment of OS.

TCF-3-mediated transcription of lncRNA HNF1A-AS1 targeting oncostatin M expression inhibits epithelial-mesenchymal transition via TGFβ signaling in gastroenteropancreatic neuroendocrine neoplasms

Long noncoding RNAs play key roles in several cancers, but their potential functions in gastroenteropancreatic neuroendocrine neoplasms remain to be investigated. We performed GeneChip assay to explore differentiated lncRNAs in gastric NENs and peri-cancerous tissues. The regulation of HNF1A-AS1 on biological behavior of GEP-NENs cells and in vivo xenograft model was confirmed by CCK8, colony formation assay, transwell, western blot and qRT-PCR. We next detected the potential transcription factors and the binding sites between them with bioinformatic analysis. qRT-PCR was performed to analyze the exact relationship between them. HNF1A-AS1 expression was decreased in gastric NENs tissues (p < 0.01). Over-expression of HNF1A-AS1 suppressed cellular proliferation, migration and invasion. Knockdown of transcription factor 3 inhibited the expression of HNF1A-AS1 and promoted cellular migration and invasion. Oncostatin M was identified as the downstream target of HNF1A-AS1. Inhibition of transforming growth factor-β activity inhibited HNF1A-AS1/Oncostatin M-mediated epithelial-mesenchymal transition. Our data suggest that transcription factor 3/HNF1A-AS1/Oncostatin M axis inhibits the tumorigenesis and metastasis of gastroenteropancreatic neuroendocrine neoplasms via transforming growth factor-β signaling.

Identification of Tumor Microenvironment-Related Prognostic Biomarkers in Luminal Breast Cancer

Background: The tumor microenvironment (TME) has been reported to have significant value in the diagnosis and prognosis of cancers. This study aimed to identify key biomarkers in the TME of luminal breast cancer (BC).

Methods: We obtained immune scores (ISs) and stromal scores (SSs) for The Cancer Genome Atlas (TCGA) luminal BC cohort from the online ESTIMATE (Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data) portal. The relationships between ISs and SSs and the overall survival of luminal BC patients were assessed by the Kaplan-Meier method. The differentially expressed messenger RNAs (DEmRNAs) related to the ISs and SSs were subjected to functional enrichment analysis. Additionally, a competing endogenous RNA (ceRNA) network was constructed with differentially expressed microRNAs (DEmiRNAs) and long noncoding RNAs (DElncRNAs). Furthermore, a protein–protein interaction (PPI) network was established to analyze the DEmRNAs in the ceRNA network. Then, survival analysis of biomarkers involved in the ceRNA network was carried out to explore their prognostic value. Finally, these biomarkers were validated using the luminal BC dataset from the Gene Expression Omnibus (GEO) database.

Results: The results showed that ISs were significantly associated with longer survival times of luminal BC patients. Functional enrichment analysis showed that the DEmRNAs were mainly associated with immune response, antigen binding, and the extracellular region. In the PPI network, the top 10 DEmRNAs were identified as hub genes that affected the TME of luminal BC. Finally, two DEmiRNAs, two DElncRNAs, and 17 DEmRNAs of the ceRNA network associated with the TME were shown to have prognostic value. Subsequently, the expression of 15 prognostic biomarkers was validated in one additional dataset (GSE81002). In particular, one lncRNA (GVINP1) and five mRNAs (CCDC69, DOCK2, IKZF1, JCHAIN, and NCKAP1L) were novel biomarkers.

Conclusions: Our studies demonstrated that ISs were associated with the survival of luminal BC patients, and a set of novel biomarkers that might play a prognostic role in the TME of luminal BC was identified.

An integrative, genomic, transcriptomic and network-assisted study to identify genes associated with human cleft lip with or without cleft palate

BACKGROUND

Cleft lip with or without cleft palate (CL/P) is one of the most common congenital human birth defects. A combination of genetic and epidemiology studies has contributed to a better knowledge of CL/P-associated candidate genes and environmental risk factors. However, the etiology of CL/P remains not fully understood. In this study, to identify new CL/P-associated genes, we conducted an integrative analysis using our in-house network tools, dmGWAS [dense module search for Genome-Wide Association Studies (GWAS)] and EW_dmGWAS (Edge-Weighted dmGWAS), in a combination with GWAS data, the human protein-protein interaction (PPI) network, and differential gene expression profiles.

RESULTS

A total of 87 genes were consistently detected in both European and Asian ancestries in dmGWAS. There were 31.0% (27/87) showed nominal significance with CL/P (gene-based p < 0.05), with three genes showing strong association signals, including KIAA1598, GPR183, and ZMYND11 (p < 1 × 10^- 3). In EW_dmGWAS, we identified 253 and 245 module genes associated with CL/P for European ancestry and the Asian ancestry, respectively. Functional enrichment analysis demonstrated that these genes were involved in cell adhesion, protein localization to the plasma membrane, the regulation of the apoptotic signaling pathway, and other pathological conditions. A small proportion of genes (5.1% for European ancestry; 2.4% for Asian ancestry) had prior evidence in CL/P as annotated in CleftGeneDB database. Our analysis highlighted nine novel CL/P candidate genes (BRD1, CREBBP, CSK, DNM1L, LOR, PTPN18, SND1, TGS1, and VIM) and 17 previously reported genes in the top modules.

CONCLUSIONS

The genes identified through superimposing GWAS signals and differential gene expression profiles onto human PPI network, as well as their functional features, helped our understanding of the etiology of CL/P. Our multi-omics integrative analyses revealed nine novel candidate genes involved in CL/P.