Frequent MYC coamplification and DNA hypomethylation of multiple genes on 8q in 8p11-p12-amplified breast carcinomas

Pan-cancer analysis of the prognostic and immunological roles of SHP-1/ptpn6

SHP-1, a nonreceptor protein tyrosine phosphatase encoded by ptpn6, has been regarded as a regulatory protein of hematopoietic cell biology for years. However, there is now increasing evidence to support its role in tumors. Thus, the role of ptpn6 for prognosis and immune regulation across 33 tumors was investigated, aiming to explore its functional heterogeneity and clinical significance in pan-cancer. Differential expression of ptpn6 was found between cancer and adjacent normal tissues, and its expression was significantly correlated with the prognosis of tumor patients. In most cancers, ptpn6 expression was significantly associated with immune infiltration. This was further confirmed by ptpn6-related genes/proteins enrichment analysis. Additionally, genetic alterations in ptpn6 was observed in most cancers. As for epigenetic changes, it's phosphorylation levels significantly altered in 6 tumors, while methylation levels significantly altered in 12 tumors. Notably, the methylation levels of ptpn6 were significantly decreased in 11 tumors, accompanied by its increased expression in 8 of them, suggesting that the hypomethylation may be related to its increased expression. Our results show that ptpn6 plays a specific role in tumor immunity and exerts a pleiotropic effect in a variety of tumors. It can serve as a prognostic factor for some cancers. Especially in LGG, KIRC, UCS and TGCT, the increased expression of ptpn6 is associated with poor prognosis and high immune infiltration. This aids in understanding the role of ptpn6 in tumor biology, and can provide insight into presenting a potential biomarker for poor prognosis and immune infiltration in cancers.

Characterization of the Rat Osteosarcoma Cell Line UMR-106 by Long-Read Technologies Identifies a Large Block of Amplified Genes Associated with Human Disease

BACKGROUND/OBJECTIVES

The rat osteosarcoma cell line UMR-106 is widely used for the study of bone cancer biology but it has not been well characterized with modern genomic methods.

METHODS

To better understand the biology of UMR-106 cells we used a combination of optical genome mapping (OGM), long-read sequencing nanopore sequencing and RNA sequencing.The UMR-106 genome was compared to a strain-matched Sprague-Dawley rat for variants associated with human osteosarcoma while expression data were contrasted with a public osteoblast dataset.

RESULTS

Using the COSMIC database to identify the most affected genes in human osteosarcomas we found somatic mutations in Tp53 and H3f3a. OGM identified a relatively small number of differences between the cell line and a strain-matched control animal but did detect a ~45 Mb block of amplification that included Myc on chromosome 7 which was confirmed by long-read sequencing. The amplified region showed several blocks of non-contiguous rearranged sequence implying complex rearrangements during their formation and included 14 genes reported as biomarkers in human osteosarcoma, many of which also showed increased transcription. A comparison of 5mC methylation from the nanopore reads of tumor and control samples identified genes with distinct differences including the OS marker Cdkn2a.

CONCLUSIONS

This dataset illustrates the value of long DNA methods for the characterization of cell lines and how inter-species analysis can inform us about the genetic nature underlying mutations that underpin specific tumor types. The data should be a valuable resource for investigators studying osteosarcoma, in general, and specifically the UMR-106 model.

Identification and validation of genes associated with prognosis of cisplatin-resistant ovarian cancer

PURPOSE

To investigate the role of prognostic genes related to cisplatin resistance in ovarian cancer during disease progression.

METHOD

The gene expression profile of the NCI-60 cell line was acquired through comprehensive analysis of the GEO database accession GSE116439. We performed a thorough analysis of gene expression differences in samples from seven individuals exposed to cisplatin concentrations of 0 nM compared to seven samples exposed to 15000 nM over a 24-h period. Key genes were initially identified through LASSO regression, followed by their enrichment through differential gene function analysis (GO) and pathway enrichment analysis (KEGG). Subsequently, a prognostic risk model was established for these key genes. The prognostic model's performance was assessed through K-M survival curves and ROC curves. To examine the variance in immune cell infiltration between the high and low-risk groups, CIBERSORTx analysis was employed. Finally, validation of prognostic gene expression in cisplatin-resistant ovarian cancer was carried out using clinical samples, employing RT-qPCR and Western Blot techniques.

RESULTS

A total of 132 differential genes were found between cisplatin resistance and control group, and 8 key prognostic genes were selected by analysis, namely VPS13B, PLGRKT, CDKAL1, TBC1D22A, TAP1, PPP3CA, CUX1 and PPP1R15A. The efficacy of the risk assessment model derived from prognostic biomarkers, as indicated by favorable performance on both Kaplan-Meier survival curves and ROC curves. Significant variations in the abundance of Macrophages M1, T cells CD4 memory resting, T cells follicular helper, and T cells gamma delta were observed between the high and low-risk groups. To further validate our findings, RT-qPCR and Western Blot analyses were employed, confirming differential expression of the identified eight key genes between the two groups.

CONCLUSION

VPS13B, TBC1D22A, PPP3CA, CUX1 and PPP1R15A were identified as poor prognostic genes of cisplatin resistance in ovarian cancer, while PLGRKT, CDKAL1 and TAP1 were identified as good prognostic genes. This offers a novel perspective for future advancements in ovarian cancer treatment, suggesting potential avenues for the development of new therapeutic targets.

NDRGs in Breast Cancer: A Review and In Silico Analysis

The N-myc downstream regulated gene family (NDRGs) includes four members: NDRG1, NDRG2, NDRG3, and NDRG4. These members exhibit 53-65% amino acid identity. The role of NDRGs in tumor growth and metastasis appears to be tumor- and context-dependent. While many studies have reported that these family members have tumor suppressive roles, recent studies have demonstrated that NDRGs, particularly NDRG1 and NDRG2, function as oncogenes, promoting tumor growth and metastasis. Additionally, NDRGs are involved in regulating different signaling pathways and exhibit diverse cellular functions in breast cancers. In this review, we comprehensively outline the oncogenic and tumor suppressor roles of the NDRG family members in breast cancer, examining evidence from in vitro and in vivo breast cancer models as well as tumor tissues from breast cancer patients. We also present analyses of publicly available genomic and transcriptomic data from multiple independent cohorts of breast cancer patients.

Characterization of MYBL1 Gene in Triple-Negative Breast Cancers and the Genes' Relationship to Alterations Identified at the Chromosome 8q Loci

The MYBL1 gene is a strong transcriptional activator involved in events associated with cancer progression. Previous data show MYBL1 overexpressed in triple-negative breast cancer (TNBC). There are two parts to this study related to further characterizing the MYBL1 gene. We start by characterizing MYBL1 reference sequence variants and isoforms. The results of this study will help in future experiments in the event there is a need to characterize functional variants and isoforms of the gene. In part two, we identify and validate expression and gene-related alterations of MYBL1, VCIP1, MYC and BOP1 genes in TNBC cell lines and patient samples selected from the Breast Invasive Carcinoma TCGA 2015 dataset available at cBioPortal.org. The four genes are located at chromosomal regions 8q13.1 to 8q.24.3 loci, regions previously identified as demonstrating a high percentage of alterations in breast cancer. We identify alterations, including changes in expression, deletions, amplifications and fusions in MYBL1, VCPIP1, BOP1 and MYC genes in many of the same patients, suggesting the panel of genes is involved in coordinated activity in patients. We propose that MYBL1, VCPIP1, MYC and BOP1 collectively be considered as genes associated with the chromosome 8q loci that potentially play a role in TNBC pathogenesis.

NSD3 in Cancer: Unraveling Methyltransferase-Dependent and Isoform-Specific Functions

NSD3 (nuclear receptor-binding SET domain protein 3) is a member of the NSD histone methyltransferase family of proteins. In recent years, it has been identified as a potential oncogene in certain types of cancer. The NSD3 gene encodes three isoforms, the long version (NSD3L), a short version (NSD3S) and the WHISTLE isoforms. Importantly, the NSD3S isoform corresponds to the N-terminal region of the full-length protein, lacking the methyltransferase domain. The chromosomal location of NSD3 is frequently amplified across cancer types, such as breast, lung, and colon, among others. Recently, this amplification has been correlated to a chromothripsis event, that could explain the different NSD3 alterations found in cancer. The fusion proteins containing NSD3 have also been reported in leukemia (NSD3-NUP98), and in NUT (nuclear protein of the testis) midline carcinoma (NSD3-NUT). Its role as an oncogene has been described by modulating different cancer pathways through its methyltransferase activity, or the short isoform of the protein, through protein interactions. Specifically, in this review we will focus on the functions that have been characterized as methyltransferase dependent, and those that have been correlated with the expression of the NSD3S isoform. There is evidence that both the NSD3L and NSD3S isoforms are relevant for cancer progression, establishing NSD3 as a therapeutic target. However, further functional studies are needed to differentiate NSD3 oncogenic activity as dependent or independent of the catalytic domain of the protein, as well as the contribution of each isoform and its clinical significance in cancer progression.

The emerging role of fatty acid binding protein 5 (FABP5) in cancers

Fatty acid binding protein 5 (FABP5, or epidermal FABP) is an intracellular chaperone of fatty acid molecules that regulates lipid metabolism and cell growth. In patient-derived tumours, FABP5 expression is increased up to tenfold, often co-expressed with other cancer-related proteins. High tumoral FABP5 expression is associated with poor prognosis. FABP5 activates transcription factors (TFs) leading to increased expression of proteins involved in tumorigenesis. Genetic and pharmacological preclinical studies show that inhibiting FABP5 reduces protumoral markers, whereas elevation of FABP5 promotes tumour growth and spread. Thus, FABP5 might be a valid target for novel therapeutics. The evidence base is currently strongest for liver, prostate, breast, and brain cancers, and squamous cell carcinoma (SCC), which could represent relevant patient populations for any drug discovery programme. Teaser: This review presents the growing evidence that upregulated fatty acid binding protein 5 (FABP5) plays a role in the progression of multiple cancer types, and may represent a novel therapeutic target.

LncRNA BC promotes lung adenocarcinoma progression by modulating IMPAD1 alternative splicing

BACKGROUND

The therapeutic value of targeted therapies in patients with lung cancer is reduced when tumours acquire secondary resistance after an initial period of successful treatment. However, the molecular events behind the resistance to targeted therapies in lung cancer remain largely unknown.

AIMS

To discover the important role and mechanism of lncRNA BC in promoting tumor metastasis and influencing clinical prognosis of LUAD.

MATERIALS & METHODS

Microarrays were used to screen a comprehensive set of lncRNAs with differential expression profiles in lung cancer cells. The functional role and mechanism of lncRNA were further investigated by gain- and loss-of-function assays. RNA pull-down, protein assays, and mass spectrometry were used to identify proteins that interacted with lncRNA. TaqMan PCR was used to measure lncRNA in lung adenocarcinoma and adjacent nontumor tissues from 428 patients. The clinical significance of lncRNA identified was statistically confirmed in this cohort of patients.

RESULTS

In this study, we show that the long non-coding RNA BC009639 (BC) is involved in acquired resistance to EGFR-targeted therapies. Among the 235 long non-coding RNAs that were differentially expressed in lung cancer cell lines, with different metastatic potentials, BC promoted growth, invasion, metastasis, and resistance to EGFR-tyrosine kinase inhibitors (EGFR-TKIs), both in vitro and in vivo. BC was highly expressed in 428 patients with lung adenocarcinoma (LUAD) and high BC expression correlated with reduced efficacy of EGFR-TKI therapy. To uncover the molecular mechanism of BC-mediated EGFR-TKI resistance in lung cancer, we screened and identified nucleolin and hnRNPK that interact with BC. BC formed the splicing complex with nucleolin and hnRNPK to facilitate the production of a non-protein-coding inositol monophosphatase domain containing 1 (IMPAD1) splice variant, instead of the protein-coding variant. The BC-mediated alternative splicing (AS) of IMPAD1 resulted in the induction of the epithelial-mesenchymal transition and resistance to EGFR-TKI in lung cancer. High BC expression correlated with clinical progress and poor survival among 402 patients with LUAD.

DISSCUSSION

Through alternative splicing, BC boosted the non-coding IMPAD1-203 transcript variant while suppressing the IMPAD1-201 variant. In order to control the processing of pre-mRNA, BC not only attracted RNA binding proteins (NCL, IGF2BP1) or splicing factors (hnRNPK), but also controlled the formation of the splicing-regulator complex by creating RNA-RNA-duplexes.

CONCLUSION

Our results reveal an important role for BC in mediating resistance to EGFR-targeted therapy in LUAD through IMPAD1 AS and in implication for the targeted therapy resistance.

SQLE inhibition suppresses the development of pancreatic ductal adenocarcinoma and enhances its sensitivity to chemotherapeutic agents in vitro

PURPOSE

In this study, we sought to explore the function of seven important enzymes (MSMO1, EBP, HMGCS1, IDI2, DHCR7, FDFT1, and SQLE) involved in cholesterol biosynthesis especially SQLE in PDAC therapy.

METHODS AND RESULTS

The TCGA and Oncomine dataset were used to explore the expression of the seven enzymes in normal and cancerous pancreatic individual, and their anti-proliferation efficiency against PDAC cells was measured by cell viability assay. Expression level and prognostic values of SQLE were evaluated by western blot and Kaplan-Meier analysis. The influence of SQLE knockdown by shRNA in PDAC cells was assessed by transwell, colony formation and cell cycle analysis. RNA-seq and GSEA were utilized to investigate the potential mechanisms. The synergistic effect of SQLE inhibitor, terbinafine, combined with six chemotherapeutic drugs in PDAC cells was tested by CCK-8 method. We demonstrated that downregulation of those enzymes especially SQLE significantly suppressed PDAC cells survival. SQLE was upregulated in PDAC cell lines, and the elevated level of SQLE was correlated with poor prognosis in pancreatic cancer samples. SQLE knockdown could significantly inhibit the proliferation and migration of PDAC cells. Cell cycle was blocked in S phase after SQLE silencing. Mechanistically, GSEA analysis with RNA-seq data revealed that SQLE silencing negatively mediated mTORC1 and TNFα/NF-κB signaling pathways. Besides, SQLE inhibitor terbinafine enhanced chemotherapeutic sensitivity of the six compounds.

CONCLUSIONS

This study demonstrated that SQLE was a novel target for PDAC therapy. The synergism role of SQLE inhibition and chemotherapy may be potential therapeutic strategy for pancreatic cancer treatment.

ENPP2 Promoter Methylation Correlates with Decreased Gene Expression in Breast Cancer: Implementation as a Liquid Biopsy Biomarker

Autotaxin (ATX), encoded by the ctonucleotide pyrophosphatase/phosphodiesterase 2 (ENPP2) gene, is a key enzyme in lysophosphatidic acid (LPA) synthesis. We have recently described ENPP2 methylation profiles in health and multiple malignancies and demonstrated correlation to its aberrant expression. Here we focus on breast cancer (BrCa), analyzing in silico publicly available BrCa methylome datasets, to identify differentially methylated CpGs (DMCs) and correlate them with expression. Numerous DMCs were identified between BrCa and healthy breast tissues in the gene body and promoter-associated regions (PA). PA DMCs were upregulated in BrCa tissues in relation to normal, in metastatic BrCa in relation to primary, and in stage I BrCa in relation to normal, and this was correlated to decreased mRNA expression. The first exon DMC was also investigated in circulating cell free DNA (ccfDNA) isolated by BrCa patients; methylation was increased in BrCa in relation to ccfDNA from healthy individuals, confirming in silico results. It also differed between patient groups and was correlated to the presence of multiple metastatic sites. Our data indicate that promoter methylation of ENPP2 arrests its transcription in BrCa and introduce first exon methylation as a putative biomarker for diagnosis and monitoring which can be assessed in liquid biopsy.

Progesterone receptor isoform ratio dictates antiprogestin/progestin effects on breast cancer growth and metastases: A role for NDRG1

Progesterone receptors (PR) ligands are being tested in luminal breast cancer. There are mainly two PR isoforms, PRA and PRB, and their ratio (PRA/PRB) may be predictive of antiprogestin response. Our aim was to investigate: the impact of the PR isoform ratio on metastatic behavior, the PR isoform ratio in paired primary tumors and lymph node metastases (LNM) and, the effect of antiprogestin/progestins on metastatic growth. Using murine and human metastatic models, we demonstrated that tumors with PRB > PRA (PRB-H) have a higher proliferation index but less metastatic ability than those with PRA > PRB (PRA-H). Antiprogestins and progestins inhibited metastatic burden in PRA-H and PRB-H models, respectively. In breast cancer samples, LNM retained the same PRA/PRB ratio as their matched primary tumors. Moreover, PRA-H LNM expressed higher total PR levels than the primary tumors. The expression of NDRG1, a metastasis suppressor protein, was higher in PRB-H compared with PRA-H tumors and was inversely regulated by antiprogestins/progestins. The binding of the corepressor SMRT at the progesterone responsive elements of the NDRG1 regulatory sequences, together with PRA, impeded its expression in PRA-H cells. Antiprogestins modulate the interplay between SMRT and AIB1 recruitment in PRA-H or PRB-H contexts regulating NDRG1 expression and thus, metastasis. In conclusion, we provide a mechanistic interpretation to explain the differential role of PR isoforms in metastatic growth and highlight the therapeutic benefit of using antiprogestins in PRA-H tumors. The therapeutic effect of progestins in PRB-H tumors is suggested. This article is protected by copyright. All rights reserved.

ENPP2 Methylation in Health and Cancer

Autotaxin (ATX) encoded by Ectonucleotide Pyrophosphatase/Phosphodiesterase 2 (ENPP2) is a key enzyme in Lysophosphatidic Acid (LPA) synthesis implicated in cancer. Although its aberrant expression has been reported, ENPP2 methylation profiles in health and malignancy are not described. We examined in silico the methylation of ENPP2 analyzing publicly available methylome datasets, to identify Differentially Methylated CpGs (DMCs) which were then correlated with expression at gene and isoform levels. Significance indication was set to be FDR corrected p-value < 0.05. Healthy tissues presented methylation in all gene body CGs and lower levels in Promoter Associated (PA) regions, whereas in the majority of the tumors examined (HCC, melanoma, CRC, LC and PC) the methylation pattern was reversed. DMCs identified in the promoter were located in sites recognized by multiple transcription factors, suggesting involvement in gene expression. Alterations in methylation were correlated to an aggressive phenotype in cancer cell lines. In prostate and lung adenocarcinomas, increased methylation of PA CGs was correlated to decreased ENPP2 mRNA expression and to poor prognosis parameters. Collectively, our results corroborate that methylation is an active level of ATX expression regulation in cancer. Our study provides an extended description of the methylation status of ENPP2 in health and cancer and points out specific DMCs of value as prognostic biomarkers.

Cholesterol Metabolic Reprogramming in Cancer and Its Pharmacological Modulation as Therapeutic Strategy

Cholesterol is a ubiquitous sterol with many biological functions, which are crucial for proper cellular signaling and physiology. Indeed, cholesterol is essential in maintaining membrane physical properties, while its metabolism is involved in bile acid production and steroid hormone biosynthesis. Additionally, isoprenoids metabolites of the mevalonate pathway support protein-prenylation and dolichol, ubiquinone and the heme a biosynthesis. Cancer cells rely on cholesterol to satisfy their increased nutrient demands and to support their uncontrolled growth, thus promoting tumor development and progression. Indeed, transformed cells reprogram cholesterol metabolism either by increasing its uptake and de novo biosynthesis, or deregulating the efflux. Alternatively, tumor can efficiently accumulate cholesterol into lipid droplets and deeply modify the activity of key cholesterol homeostasis regulators. In light of these considerations, altered pathways of cholesterol metabolism might represent intriguing pharmacological targets for the development of exploitable strategies in the context of cancer therapy. Thus, this work aims to discuss the emerging evidence of in vitro and in vivo studies, as well as clinical trials, on the role of cholesterol pathways in the treatment of cancer, starting from already available cholesterol-lowering drugs (statins or fibrates), and moving towards novel potential pharmacological inhibitors or selective target modulators.

RRM2B Is Frequently Amplified Across Multiple Tumor Types: Implications for DNA Repair, Cellular Survival, and Cancer Therapy

RRM2B plays a crucial role in DNA replication, repair and oxidative stress. While germline RRM2B mutations have been implicated in mitochondrial disorders, its relevance to cancer has not been established. Here, using TCGA studies, we investigated RRM2B alterations in cancer. We found that RRM2B is highly amplified in multiple tumor types, particularly in MYC-amplified tumors, and is associated with increased RRM2B mRNA expression. We also observed that the chromosomal region 8q22.3–8q24, is amplified in multiple tumors, and includes RRM2B, MYC along with several other cancer-associated genes. An analysis of genes within this 8q-amplicon showed that cancers that have both RRM2B-amplified along with MYC have a distinct pattern of amplification compared to cancers that are unaltered or those that have amplifications in RRM2B or MYC only. Investigation of curated biological interactions revealed that gene products of the amplified 8q22.3–8q24 region have important roles in DNA repair, DNA damage response, oxygen sensing, and apoptosis pathways and interact functionally. Notably, RRM2B-amplified cancers are characterized by mutation signatures of defective DNA repair and oxidative stress, and at least RRM2B-amplified breast cancers are associated with poor clinical outcome. These data suggest alterations in RR2MB and possibly the interacting 8q-proteins could have a profound effect on regulatory pathways such as DNA repair and cellular survival, highlighting therapeutic opportunities in these cancers.

MYC Enhances Cholesterol Biosynthesis and Supports Cell Proliferation Through SQLE

Oncogene c-Myc (referred in this report as MYC) promotes tumorigenesis in multiple human cancers. MYC regulates numerous cellular programs involved in cell growth and cell metabolism. Tumor cells exhibit obligatory dependence on cholesterol metabolism, which provides essential membrane components and metabolites to support cell growth. To date, how cholesterol biosynthesis is delicately regulated to promote tumorigenesis remains unclear. Here, we show that MYC enhances cholesterol biosynthesis and promotes cell proliferation. Through transcriptional upregulation of SQLE, a rate-limiting enzyme in cholesterol synthesis pathway, MYC increases cholesterol production and promotes tumor cell growth. SQLE overexpression restores the cellular cholesterol levels in MYC-knockdown cells. More importantly, in SQLE-depleted cells, enforced expression of MYC has no effect on cholesterol levels. Therefore, our findings reveal that SQLE is critical for MYC-mediated cholesterol synthesis, and further demonstrate that SQLE may be a potential therapeutic target in MYC-amplified cancers.

Squalene epoxidase expression is associated with breast tumor progression and with a poor prognosis in breast cancer

Differentially expressed genes (DEGs) have been previously identified using massive parallel RNA sequencing in matched normal, breast cancer (BC) and nodal metastatic tissues. Squalene epoxidase (SQLE), one of these DEGs, is a key enzyme in cholesterol synthesis. The aim of the present study was to investigate the potential involvement of SQLE in the tumorigenic process of BC and to determine its association with the clinical outcome of BC. SQLE mRNA expression was measured using reverse transcription-quantitative PCR in 10 pairs of ductal carcinoma in situ (DCIS) and BC tissues and their adjacent normal tissues. Immunohistochemical staining of SQLE on tissue microarray was performed in 26 normal breast, 79 DCIS and 198 BC samples. The role of SQLE as a prognostic biomarker in patients with BC has been verified using BreastMark. SQLE mRNA expression was significantly increased in DCIS and BC tissues compared with that in their adjacent normal tissues. High SQLE expression was detected in 0, 48.1 and 40.4% of normal breast, DCIS and BC tissues, respectively. SQLE expression in DCIS and BC tissues was significantly higher than that in normal breast tissues. High SQLE expression was observed in DCIS with higher nuclear grade, comedo-type necrosis and HER2 positivity. High SQLE expression in BC was associated with larger tumor size, nodal metastases, higher stage, HER2 subtype and distant metastatic relapse. High SQLE expression was associated with poor disease-free and overall survival, and independently predicted poor disease-free survival in patients with BC. Following BreastMark analysis, high SQLE mRNA expression in BC was significantly associated with a poor prognosis in the ‘all’, lymph node negative, lymph node positive, luminal A subtype and luminal B subtype groups. Therefore, SQLE expression may be upregulated during the tumorigenic process of BC, and high SQLE expression may be a useful biomarker for predicting a poor prognosis in patients with BC.

A Novel Long Non-Coding RNA lnc030 Maintains Breast Cancer Stem Cell Stemness by Stabilizing SQLE mRNA and Increasing Cholesterol Synthesis

Cancer stem cells (CSCs) are considered the roots of cancer metastasis and recurrence (CSCs), due in part to their self-renewal and therapy resistance properties. However, the underlying mechanisms for the regulation of CSC stemness are poorly understood. Recently, increasing evidence shows that long non-coding RNAs (lncRNAs) are critical regulators for cancer cell function in various malignancies including breast cancer, but how lncRNAs regulate the function of breast cancer stem cells (BCSCs) remains to be determined. Herein, using lncRNA/mRNA microarray assays, a novel lncRNA (named lnc030) is identified, which is highly expressed in BCSCs in vitro and in vivo, as a pivotal regulator in maintaining BCSC stemness and promoting tumorigenesis. Mechanistically, lnc030 cooperates with poly(rC) binding protein 2(PCBP2) to stabilize squalene epoxidase (SQLE) mRNA, resulting in an increase of cholesterol synthesis. The increased cholesterol in turn actives PI3K/Akt signaling, which governs BCSC stemness. In summary, these findings demonstrate that a new, lnc030-based mechanism for regulating cholesterol synthesis and stemness properties of BCSCs. The lnc030-SQLE-cholesterol synthesis pathway may serve as an effective therapeutic target for BCSC elimination and breast cancer treatment.

High Squalene Epoxidase in Tumors Predicts Worse Survival in Patients With Hepatocellular Carcinoma: Integrated Bioinformatic Analysis on NAFLD and HCC

This study aimed to identify candidate biomarkers for predicting outcomes in nonalcoholic fatty liver disease (NAFLD) and hepatocellular carcinoma (HCC). Using Gene Expression Omnibus and The Cancer Genome Atlas (TCGA) databases, we identified common upregulated differential expressed genes (DEGs) in patients with NAFLD/nonalcoholic steatohepatitis (NASH) and HCC and conducted survival analysis of these upregulated DEGs with HCC outcomes. Two common upregulated DEGs including squalene epoxidase (SQLE) and EPPK1 messenger RNA (mRNA) were significantly upregulated in NAFLD, NASH, and HCC tissues, both in GSE45436 (P < .001) and TCGA profile (P < .001). Both SQLE and EPPK1 mRNA were upregulated in 15.56% and 8.06% patients with HCC in TCGA profile. Overexpression of SQLE in tumors was significantly associated with worse overall survival (OS) and disease-free survival (DFS) in patients with HCC (log-rank P = .027 and log-rank P = .048, respectively), while no statistical significances of OS and DFS were found in EPPK1 groups (both log-rank P > .05). For validation, SQLE upregulation contributed to significantly worse OS in patients wih HCC using Kaplan-Meier plotter analysis (hazard ratio = 1.43, 95% confidence interval: 1.01-2.02, log-rank P = .043). In addition, high level of SQLE significantly associated with advanced neoplasm histologic grade, advanced AJCC stage, and α-fetoprotein elevation (P = .036, .045, and .029, respectively). Squalene epoxidase is associated with OS and DFS and serves as a novel prognostic biomarker for patients with HCC.

Personalized cancer therapy prioritization based on driver alteration co-occurrence patterns

Identification of actionable genomic vulnerabilities is key to precision oncology. Utilizing a large-scale drug screening in patient-derived xenografts, we uncover driver gene alteration connections, derive driver co-occurrence (DCO) networks, and relate these to drug sensitivity. Our collection of 53 drug-response predictors attains an average balanced accuracy of 58% in a cross-validation setting, rising to 66% for a subset of high-confidence predictions. We experimentally validated 12 out of 14 predictions in mice and adapted our strategy to obtain drug-response models from patients’ progression-free survival data. Our strategy reveals links between oncogenic alterations, increasing the clinical impact of genomic profiling.

IMPAD1 and KDELR2 drive invasion and metastasis by enhancing Golgi-mediated secretion

Non-small cell lung cancer (NSCLC) is the deadliest form of cancer worldwide, due in part to its proclivity to metastasize. Identifying novel drivers of invasion and metastasis holds therapeutic potential for the disease. We conducted a gain-of-function invasion screen, which identified two separate hits, IMPAD1 and KDELR2, as robust, independent drivers of lung cancer invasion and metastasis. Given that IMPAD1 and KDELR2 are known to be localized to the ER-Golgi pathway, we studied their common mechanism of driving in vitro invasion and in vivo metastasis and demonstrated that they enhance Golgi-mediated function and secretion. Therapeutically inhibiting matrix metalloproteases (MMPs) suppressed both IMPAD1- and KDELR2-mediated invasion. The hits from this unbiased screen and the mechanistic validation highlight Golgi function as one of the key cellular features altered during invasion and metastasis.

Lipids and cancer: Emerging roles in pathogenesis, diagnosis and therapeutic intervention

With the advent of effective tools to study lipids, including mass spectrometry-based lipidomics, lipids are emerging as central players in cancer biology. Lipids function as essential building blocks for membranes, serve as fuel to drive energy-demanding processes and play a key role as signaling molecules and as regulators of numerous cellular functions. Not unexpectedly, cancer cells, as well as other cell types in the tumor microenvironment, exploit various ways to acquire lipids and extensively rewire their metabolism as part of a plastic and context-dependent metabolic reprogramming that is driven by both oncogenic and environmental cues. The resulting changes in the fate and composition of lipids help cancer cells to thrive in a changing microenvironment by supporting key oncogenic functions and cancer hallmarks, including cellular energetics, promoting feedforward oncogenic signaling, resisting oxidative and other stresses, regulating intercellular communication and immune responses. Supported by the close connection between altered lipid metabolism and the pathogenic process, specific lipid profiles are emerging as unique disease biomarkers, with diagnostic, prognostic and predictive potential. Multiple preclinical studies illustrate the translational promise of exploiting lipid metabolism in cancer, and critically, have shown context dependent actionable vulnerabilities that can be rationally targeted, particularly in combinatorial approaches. Moreover, lipids themselves can be used as membrane disrupting agents or as key components of nanocarriers of various therapeutics. With a number of preclinical compounds and strategies that are approaching clinical trials, we are at the doorstep of exploiting a hitherto underappreciated hallmark of cancer and promising target in the oncologist's strategy to combat cancer.

Aberrantly Methylated cfDNA in Body Fluids as a Promising Diagnostic Tool for Early Detection of Breast Cancer

Breast malignancies are the leading type of cancer among women. Its prevention and early detection, particularly in young women, remains challenging. To this end, cell-free DNA (cfDNA) detected in body fluids demonstrates great potential for early detection of tissue transformation and altered molecular setup, such as epigenetic profiles. Aberrantly methylated cfDNA in body fluids could therefore serve as a potential diagnostic and prognostic tool in breast cancer management. Abnormal methylation may lead to both an activation of oncogenes via hypomethylation and an inactivation of tumor suppressor genes by hypermethylation. We update the state of the art in the area of aberrant cfDNA methylation analyses as a diagnostic and prognostic tool in breast cancer, report on the main technological challenges, and provide an outlook for advancing the overall management of breast malignancies based on cfDNA as a target for diagnosis and tailored therapies.

RNA-Seq analysis of the pathogenesis of STZ-induced male diabetic mouse liver

OBJECTIVE

Diabetes mellitus (DM) is a chronic disease characterized by hyperglycemia resulting from defects in insulin secretion, insulin action, or both. The liver is a key organ involved in glucose metabolism, and the major target proteins' changes in the pathogenesis are still unknown.

METHODS

A diabetic mouse model was induced by intraperitoneal injection of streptozotocin (STZ) solution and the RNA-Seq analysis was used to evaluate the transcription differences in the livers of diabetic mice of this study. And then, the differentially expressed genes were validated between a normal mouse group (n = 6) and a diabetic mouse group (n = 6) using quantitative real-time PCR (qRT-PCR) and Western blotting analysis. In addition, we also constructed protein-protein interacting (PPI) networks of up-regulated and down-regulated genes.

RESULTS

Transcriptome sequencing analysis revealed 370 up-regulated differentially expressed genes and 281 down-regulated differentially expressed genes in the diabetes model. The gene ontology (GO) analysis and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis results showed that the differentially expressed genes were mainly involved in immunity, enzyme activity, metabolism, and steroid synthesis. PPI analysis results indicated that the main 15 core differential proteins (Cyp51a1, Acsl4, Ugt1a1, Stat1, Gsta2, Cbr1, Aldh1a1, Fasn, Ces1, Camk2b, Tap1, Egr1, Sqle, Lpin1, Fabp5) were involved in the pathogenesis of diabetes. The qRT-PCR results showed that expression changes of four genes (Acsl4, Stat1, Gsta2, Fabp5) were in different directions from those of RNA-Seq. Western blotting results indicated that Sqle expression change at the protein level was in opposition direction from qRT-PCR, and we speculated that Sqle may be involved in the post-transcriptional modification process.

CONCLUSIONS

Our data speculated that the pathogenesis of diabetes may be mediated mainly through steroid biosynthesis, metabolic processes, and immune responses. Further researches on these pathways may provide new targets for the prevention and treatment of diabetes.

Transcriptome analysis reveals that cyclophosphamide induces premature ovarian failure by blocking cholesterol biosynthesis pathway

AIMS

The present study aimed to investigate the effects of cyclophosphamide (Cytoxan, CTX) on premature ovarian failure (POF) in mice and its regulatory mechanisms by transcriptome analysis.

MAIN METHODS

Female C57BL/6 mice were treated with a single intraperitoneal injection of 70 mg/kg CTX. Serum levels of estradiol (E₂) and follicle stimulating hormone (FSH) were measured by enzyme-linked immunosorbent assay (ELISA), and follicular structure differences were observed by hematoxylin and eosin (H&E) staining. The main mechanism of POF was investigated by RNA-seq data, protein-protein interaction (PPI) networks and qPCR analysis.

KEY FINDINGS

The serum levels of E₂ were significantly decreased and those of FSH were significantly increased compared to the control group. The ovarian weights of the mice in the CTX group were reduced, and abnormal follicular structures were also observed in the CTX group. The RNA-seq data show that the downregulated genes were related to the cholesterol biosynthesis pathway. The PPI network and qPCR analyses further confirm that the PPAR signaling pathway and the ovarian infertility genes were also involved in blocking the cholesterol biosynthesis pathway. The differences were statistically significant.

SIGNIFICANCE

Our results indicate that CTX may exert its anti-tumor effects by inactivating the cholesterol biosynthesis pathway, and simultaneously reducing the supply of estrogen precursor materials, ultimately leading to the occurrence of POF. Our data provided a preliminary theoretical basis for resolving the clinical toxicity and side effects of CTX.

Sulfation pathways from red to green

Sulfur is present in the amino acids cysteine and methionine and in a large range of essential coenzymes and cofactors and is therefore essential for all organisms. It is also a constituent of sulfate esters in proteins, carbohydrates, and numerous cellular metabolites. The sulfation and desulfation reactions modifying a variety of different substrates are commonly known as sulfation pathways. Although relatively little is known about the function of most sulfated metabolites, the synthesis of activated sulfate used in sulfation pathways is essential in both animal and plant kingdoms. In humans, mutations in the genes encoding the sulfation pathway enzymes underlie a number of developmental aberrations, and in flies and worms, their loss-of-function is fatal. In plants, a lower capacity for synthesizing activated sulfate for sulfation reactions results in dwarfism, and a complete loss of activated sulfate synthesis is also lethal. Here, we review the similarities and differences in sulfation pathways and associated processes in animals and plants, and we point out how they diverge from bacteria and yeast. We highlight the open questions concerning localization, regulation, and importance of sulfation pathways in both kingdoms and the ways in which findings from these "red" and "green" experimental systems may help reciprocally address questions specific to each of the systems.

A 17-marker panel for global genomic instability in breast cancer

Genomic instability is a hallmark of cancer that plays a pivotal role in breast cancer development and evolution. A number of existing prognostic gene expression signatures for breast cancer are based on proliferation-related genes. Here, we identified a 17-marker panel associated with genome stability. A total of 136 primary breast carcinomas were stratified by genome stability. Matched gene expression profiles showed an innate segregation based on genome stability. We identified a 17-marker panel stratifying the training and validation cohorts into high- and low-risk patients. The 17 genes associated with genomic instability strongly impacted clinical outcome in breast cancer. Pathway analyses determined chromosome organisation, cell cycle regulation, and RNA processing as the underlying biological processes, thereby offering options for drug development and treatment tailoring. Our work supports the applicability of the 17-marker panel to improve clinical outcome prediction for breast cancer patients based on a signature accounting for genomic instability.

High-Resolution Bisulfite-Sequencing of Peripheral Blood DNA Methylation in Early-Onset and Familial Risk Breast Cancer Patients

PURPOSE

Understanding and explaining hereditary predisposition to cancer has focused on the genetic etiology of the disease. However, mutations in known genes associated with breast cancer, such as BRCA1 and BRCA2, account for less than 25% of familial cases of breast cancer. Recently, specific epigenetic modifications at BRCA1 have been shown to promote hereditary breast cancer, but the broader potential for epigenetic contribution to hereditary breast cancer is not yet well understood.

EXPERIMENTAL DESIGN

We examined DNA methylation through deep bisulfite sequencing of CpG islands and known promoter or regulatory regions in peripheral blood DNA from 99 patients with familial or early-onset breast or ovarian cancer, 6 unaffected BRCA mutation carriers, and 49 unaffected controls.

RESULTS

In 9% of patients, we observed altered methylation in the promoter regions of genes known to be involved in cancer, including hypermethylation at the tumor suppressor PTEN and hypomethylation at the proto-oncogene TEX14. These alterations occur in the form of allelic methylation that span up to hundreds of base pairs in length.

CONCLUSIONS

Our observations suggest a broader role for DNA methylation in early-onset, familial risk breast cancer. Further studies are warranted to clarify these mechanisms and the benefits of DNA methylation screening for early risk prediction of familial cancers.

Importance of Copy Number Alterations of FGFR1 and C-MYC Genes in Triple Negative Breast Cancer

Background

Triple negative breast cancer (TNBC) is characterized by aggressive clinical course and is unresponsive to anti-HER2 and endocrine therapy. TNBC is difficult to treat and is often lethal. Given the need to find new targets for therapy we explored clinicopathological significance of copy number gain of FGFR1 and c-MYC. Our aim was to determine the impact of FGFR1 and c-MYC copy number gain on clinical course and outcome of TNBC.

Methods

FGFR1 and c-MYC gene copy number alterations were evaluated in 78 archive TNBC samples using TaqMan based quantitative real time PCR assays.

Results

50% of samples had increased c-MYC copy number. c-MYC copy number gain was associated with TNBC in contrast to ER positive cancers. Our results showed significant correlation between c-MYC copy number gain and high grade of TNBCs. This suggests that c-MYC copy number could be an useful prognostic marker for TNBC patients. c-MYC copy number gain was associated with high pTNM stage as well as lobular and medullary tumor subtypes. 43% of samples had increased FGFR1 copy number. No correlations between FGFR1 copy number gain and clinicopathological variables were observed.

Conclusions

We identified c-MYC copy number gain as a prognostic marker for TNBC. Our results indicate that c- MYC may contribute to TNBC progression. We observed no significant association between c-MYC and/or FGFR1 copy number status and patient survival.

Aberrant promoter 2 methylation‑mediated downregulation of protein tyrosine phosphatase, non‑receptor type 6, is associated with progression of esophageal squamous cell carcinoma

The human protein tyrosine phosphatase, non‑receptor type 6 (PTPN6) gene is located on chromosome 12p13 and encodes an Mr 68,000 non‑receptor type protein‑tyrosine phosphatase. The PTPN6 gene has been considered as a candidate tumor suppressor in hematological and solid malignancies, and promoter methylation may be an epigenetic modification silencing its expression. However, the detailed role of PTPN6 and its promoter methylation status in the pathogenesis of esophageal squamous cell carcinoma (ESCC) has not been fully elucidated. The aim of the present study was to investigate PTPN6 expression in ESCC tissues and esophageal cancer cell lines, detect the effect of CpG hypermethylation on the activity of PTPN6, and additionally elucidate the role and prognostic significance of PTPN6 in ESCC tumorigenesis and progression. The expression of PTPN6 was identified to be significantly downregulated in esophageal cancer cell lines and ESCC tissues. Marked upregulation of PTPN6 was detected in 5‑aza‑2'‑deoxycytidine‑treated esophageal cancer cells, and frequent hypermethylation of the CpG sites within the P2 promoter (P2) was detected in ESCC tissues and esophageal cancer cell lines. The expression and methylation status of PTPN6 was associated with tumor node metastasis stage, pathological differentiation and lymph node metastasis in patients with ESCC. Aberrant hypermethylation of the P2 exhibited marked tumor specificity and was identified to be associated with the expression level of PTPN6. Downregulation and hypermethylation of PTPN6 were identified to be associated with poor ESCC patient survival. Furthermore, upregulation of PTPN6 inhibited the proliferation and invasion of esophageal cancer cells in vitro. The results of the present study suggest that PTPN6 may serve as a tumor suppressor in ESCC, and it may serve as a potential target for antitumor therapy.

Bottom-up, integrated -omics analysis identifies broadly dosage-sensitive genes in breast cancer samples from TCGA

The massive genomic data from The Cancer Genome Atlas (TCGA), including proteomics data from Clinical Proteomic Tumor Analysis Consortium (CPTAC), provides a unique opportunity to study cancer systematically. While most observations are made from a single type of genomics data, we apply big data analytics and systems biology approaches by simultaneously analyzing DNA amplification, mRNA and protein abundance. Using multiple genomic profiles, we have discovered widespread dosage compensation for the extensive aneuploidy observed in TCGA breast cancer samples. We do identify 11 genes that show strong correlation across all features (DNA/mRNA/protein) analogous to that of the well-known oncogene HER2 (ERBB2). These genes are generally less well-characterized regarding their role in cancer and we advocate their further study. We also discover that shRNA knockdown of these genes has an impact on cancer cell growth, suggesting a vulnerability that could be used for cancer therapy. Our study shows the advantages of systematic big data methodologies and also provides future research directions.

A multifunctional toolkit for target-directed cancer therapy

Here we present 2shRNA, a shRNA-based nanobinder, which can simultaneously attack two therapeutic targets involved in drug resistance pathways and can additionally bind accessory molecules such as cell targeting peptides or fluorophores. We create 2shRNAs designed to specifically kill HER2+ breast cancer cells in the absence of a transfecting agent.

Autotaxin in Pathophysiology and Pulmonary Fibrosis

Lysophospholipid signaling is emerging as a druggable regulator of pathophysiological responses, and especially fibrosis, exemplified by the relative ongoing clinical trials in idiopathic pulmonary fibrosis (IPF) patients. In this review, we focus on ectonucleotide pyrophosphatase-phosphodiesterase 2 (ENPP2), or as more widely known Autotaxin (ATX), a secreted lysophospholipase D (lysoPLD) largely responsible for extracellular lysophosphatidic acid (LPA) production. In turn, LPA is a bioactive phospholipid autacoid, forming locally upon increased ATX levels and acting also locally through its receptors, likely guided by ATX's structural conformation and cell surface associations. Increased ATX activity levels have been detected in many inflammatory and fibroproliferative conditions, while genetic and pharmacologic studies have confirmed a pleiotropic participation of ATX/LPA in different processes and disorders. In pulmonary fibrosis, ATX levels rise in the broncheoalveolar fluid (BALF) and stimulate LPA production. LPA engagement of its receptors activate multiple G-protein mediated signal transduction pathways leading to different responses from pulmonary cells including the production of pro-inflammatory signals from stressed epithelial cells, the modulation of endothelial physiology, the activation of TGF signaling and the stimulation of fibroblast accumulation. Genetic or pharmacologic targeting of the ATX/LPA axis attenuated disease development in animal models, thus providing the proof of principle for therapeutic interventions.

Genome-wide multi-omics profiling of the 8p11-p12 amplicon in breast carcinoma

Genomic instability contributes to the neoplastic phenotype by deregulating key cancer-related genes, which in turn can have a detrimental effect on patient outcome. DNA amplification of the 8p11-p12 genomic region has clinical and biological implications in multiple malignancies, including breast carcinoma where the amplicon has been associated with tumor progression and poor prognosis. However, oncogenes driving increased cancer-related death and recurrent genetic features associated with the 8p11-p12 amplicon remain to be identified. In this study, DNA copy number and transcriptome profiling data for 229 primary invasive breast carcinomas (corresponding to 185 patients) were evaluated in conjunction with clinicopathological features to identify putative oncogenes in 8p11-p12 amplified samples. Illumina paired-end whole transcriptome sequencing and whole-genome SNP genotyping were subsequently performed on 23 samples showing high-level regional 8p11-p12 amplification to characterize recurrent genetic variants (SNPs and indels), expressed gene fusions, gene expression profiles and allelic imbalances. We now show previously undescribed chromothripsis-like patterns spanning the 8p11-p12 genomic region and allele-specific DNA amplification events. In addition, recurrent amplification-specific genetic features were identified, including genetic variants in the HIST1H1E and UQCRHL genes and fusion transcripts containing MALAT1 non-coding RNA, which is known to be a prognostic indicator for breast cancer and stimulated by estrogen. In summary, these findings highlight novel candidate targets for improved treatment of 8p11-p12 amplified breast carcinomas.

Cigarette smoke and chewing tobacco alter expression of different sets of miRNAs in oral keratinocytes

Carcinogenic effect of tobacco in oral cancer is through chewing and/or smoking. Significant differences exist in development of oral cancer between tobacco users and non-users. However, molecular alterations induced by different forms of tobacco are yet to be fully elucidated. We developed cellular models of chronic exposure to chewing tobacco and cigarette smoke using immortalized oral keratinocytes. Chronic exposure to tobacco resulted in increased cell scattering and invasiveness in immortalized oral keratinocytes. miRNA sequencing using Illumina HiSeq 2500 resulted in the identification of 10 significantly dysregulated miRNAs (4 fold; p ≤ 0.05) in chewing tobacco treated cells and 6 in cigarette smoke exposed cells. We integrated this data with global proteomic data and identified 36 protein targets that showed inverse expression pattern in chewing tobacco treated cells and 16 protein targets that showed inverse expression in smoke exposed cells. In addition, we identified 6 novel miRNAs in chewing tobacco treated cells and 18 novel miRNAs in smoke exposed cells. Integrative analysis of dysregulated miRNAs and their targets indicates that signaling mechanisms leading to oncogenic transformation are distinct between both forms of tobacco. Our study demonstrates alterations in miRNA expression in oral cells in response to two frequently used forms of tobacco.

The Autotaxin-Lysophosphatidic Acid Axis Promotes Lung Carcinogenesis

Pathogenesis and progression of lung cancer are governed by complex interactions between the environment and host genetic susceptibility, which is further modulated by genetic and epigenetic changes. Autotaxin (ATX, ENPP2) is a secreted glycoprotein that catalyzes the extracellular production of lysophosphatidic acid (LPA), a growth-factor-like phospholipid that is further regulated by phospholipid phosphatases (PLPP). LPA's pleiotropic effects in almost all cell types are mediated through at least six G-protein coupled LPA receptors (LPAR) that exhibit overlapping specificities, widespread distribution, and differential expression profiles. Here we use both preclinical models of lung cancer and clinical samples (from patients and healthy controls) to investigate the expression levels, activity, and biological role of the above components of the ATX/LPA axis in lung cancer. ENPP2 was genetically altered in 8% of patients with lung cancer, whereas increased ATX staining and activity were detected in patient biopsies and sera, respectively. Moreover, PLPP3 expression was consistently downregulated in patients with lung cancer. Comparable observations were made in the two most widely used animal models of lung cancer, the carcinogen urethane-induced and the genetically engineered K-ras^G12D -driven models, where genetic deletion of Enpp2 or Lpar1 resulted in disease attenuation, thus confirming a procarcinogenic role of LPA signaling in the lung. Expression profiling data analysis suggested that metabolic rewiring may be implicated in the procarcinogenic effects of the ATX/LPA axis in K-ras- ^G12D -driven lung cancer pathogenesis.Significance: These findings establish the role of ATX/LPA in lung carcinogenesis, thus expanding the mechanistic links between pulmonary fibrosis and cancer. Cancer Res; 78(13); 3634-44. ©2018 AACR.

Squalene epoxidase as a promising metabolic target in cancer treatment

Oncogenic alteration of the cholesterol synthesis pathway is a recognized mechanism of metabolic adaptation. In the present review, we focus on squalene epoxidase (SE), one of the two rate-limiting enzymes in cholesterol synthesis, retracing its history since its discovery as an antimycotic target to its description as an emerging metabolic oncogene by amplification with clinical relevance in cancer. We review the published literature assessing the association between SE over-expression and poor prognosis in this disease. We assess the works demonstrating how SE promotes tumor cell proliferation and migration, and displaying evidence of cancer cell demise in presence of human SE inhibitors in in vitro and in vivo models. Taken together, robust scientific evidence has by now accumulated pointing out SE as a promising novel therapeutic target in cancer treatment.

Diagnostic and prognostic implications of ribosomal protein transcript expression patterns in human cancers

Background

Ribosomes, the organelles responsible for the translation of mRNA, are comprised of four rRNAs and ~ 80 ribosomal proteins (RPs). Although canonically assumed to be maintained in equivalent proportions, some RPs have been shown to possess differential expression across tissue types. Dysregulation of RP expression occurs in a variety of human diseases, notably in many cancers, and altered expression of some RPs correlates with different tumor phenotypes and patient survival. Little work has been done, however, to characterize overall patterns of RP transcript (RPT) expression in human cancers.

Methods

To investigate the impact of global RPT expression patterns on tumor phenotypes, we analyzed RPT expression of ~ 10,000 human tumors and over 700 normal tissues from The Cancer Genome Atlas (TCGA) using t-distributed stochastic neighbor embedding (t-SNE). Clusters of tumors identified by t-SNE were then analyzed with chi-squared and t-tests to compare phenotypic data, ANOVA to compare individual RPT expression, and Kaplan-Meier curves to assess survival differences.

Results

Normal tissues and cancers possess distinct and readily discernible RPT expression patterns that are independent of their absolute levels of expression. In tumors, RPT patterning is distinct from that of normal tissues, identifies heretofore unrecognized tumor subtypes, and in many cases correlates with molecular, pathological, and clinical features, including survival.

Conclusions

RPT expression patterns are both tissue-specific and tumor-specific. These could be used as a powerful and novel method of tumor classification, offering a potential clinical tool for prognosis and therapeutic stratification.

Electronic supplementary material

The online version of this article (10.1186/s12885-018-4178-z) contains supplementary material, which is available to authorized users.

TRPS1 Suppresses Breast Cancer Epithelial-mesenchymal Transition Program as a Negative Regulator of SUZ12

Breast cancer (BC) is among the most common malignant diseases and metastasis is the handcuff of treatment. Cancer metastasis is a multistep process associated with the epithelial-mesenchymal transition (EMT) program. Several studies have demonstrated that transcriptional repressor GATA binding 1 (TRPS1) played important roles in development and progression of primary BC. In this study we sought to identify the mechanisms responsible for this function of TRPS1 in the continuum of the metastatic cascade. Here we described that TRPS1 was significantly associated with BC metastasis using public assessable datasets. Clinically, loss of TRPS1 expression in BC was related to higher histological grade. In vitro functional study and bioinformatics analysis revealed that TRPS1 inhibited cell migration and EMT in BC. Importantly, we identified SUZ12 as a novel target of TRPS1 related to EMT program. ChIP assay demonstrated TRPS1 directly inhibited SUZ12 transcription by binding to the SUZ12 promoter. Loss of TRPS1 resulted in increased SUZ12 binding and H3K27 tri-methylation at the CDH1 promoter and repression of E-cadherin. In all, our data indicated that TRPS1 maintained the expression of E-cadherin by inhibiting SUZ12, which might provide novel insight into how loss of TRPS1 contributed to BC progression.

The exomic landscape of t(14;18)-negative diffuse follicular lymphoma with 1p36 deletion

Predominantly diffuse t(14;18) negative follicular lymphoma (FL) with 1p36 deletion shows distinctive clinical, morphological and molecular features that distinguish it from classical FL. In order to investigate whether it possesses a unique mutation profile, we performed whole exome sequencing of six well-characterised cases. Our analysis showed that the mutational landscape of this subtype is largely distinct from classical FL. It appears to harbour several recurrent mutations, affecting STAT6, CREBBP and basal membrane protein genes with high frequency. Our data support the view that this FL subtype should be considered a separate entity from classical FL.

A Novel 18-Marker Panel Predicting Clinical Outcome in Breast Cancer

Background: Gene expression profiling has made considerable contributions to our understanding of cancer biology and clinical care. This study describes a novel gene expression signature for breast cancer-specific survival that was validated using external datasets.Methods: Gene expression signatures for invasive breast carcinomas (mainly luminal B subtype) corresponding to 136 patients were analyzed using Cox regression, and the effect of each gene on disease-specific survival (DSS) was estimated. Iterative Bayesian model averaging was applied on multivariable Cox regression models resulting in an 18-marker panel, which was validated using three external validation datasets. The 18 genes were analyzed for common pathways and functions using the Ingenuity Pathway Analysis software. This study complied with the REMARK criteria.Results: The 18-gene multivariable model showed a high predictive power for DSS in the training and validation cohort and a clear stratification between high- and low-risk patients. The differentially expressed genes were predominantly involved in biological processes such as cell cycle, DNA replication, recombination, and repair. Furthermore, the majority of the 18 genes were found to play a pivotal role in cancer.Conclusions: Our findings demonstrated that the 18 molecular markers were strong predictors of breast cancer-specific mortality. The stable time-dependent area under the ROC curve function (AUC(t)) and high C-indices in the training and validation cohorts were further improved by fitting a combined model consisting of the 18-marker panel and established clinical markers.Impact: Our work supports the applicability of this 18-marker panel to improve clinical outcome prediction for breast cancer patients. Cancer Epidemiol Biomarkers Prev; 26(11); 1619-28. ©2017 AACR.

The role of c-Myc-RBM38 loop in the growth suppression in breast cancer

Background

RNA-binding protein 38 (RBM38) is a member of the RNA recognition motif (RRM) family of RNA-binding proteins (RBPs). RBM38 often exerts its function by forming regulatory loops with relevant genes. c-Myc is an oncogenic transcription factor that is upregulated in one-third of breast cancers and involved in many cellular processes in this malignancy. In our previous study, RBM38 was identified as a tumor suppressor in breast cancer. In the present study, we investigated the mechanisms underlying the regulation of this tumor suppressor.

Methods

Lentivirus transfections, Western blotting analysis, qRT-PCR and immunohistochemistry were employed to study the expression of c-Myc and RBM38. Chromatin immunoprecipitation and dual-luciferase reporter assays were performed to investigate the direct relationship between c-Myc protein and the RBM38 gene. RNA immunoprecipitation combined with dual-luciferase reporter assays was conducted to confirm the direct relationship between the RBM38 protein and the c-Myc transcript.

Results

Knockdown of c-Myc increased RBM38 expression by binding directly to specific DNA sequences (5′-CACGTG-3′), known as the E-box motif, in the promoter region of RBM38 gene. Additionally, RBM38 destabilized the c-Myc transcript by directly targeting AU-rich elements (AREs) in the 3′-untranslated region (3′-UTR) of c-Myc mRNA to suppress c-Myc expression. Moreover, specific inhibitors of c-Myc transcriptional activity inhibited RBM38-induced suppression of growth, implying that RBM38 acts as a tumor suppressor via a mechanism that depends, at least partially, on the reduction of c-Myc expression in breast cancer.

Conclusions

RBM38 and c-Myc form a unique mutually antagonistic RBM38-c-Myc loop in breast cancer.

Electronic supplementary material

The online version of this article (doi:10.1186/s13046-017-0521-5) contains supplementary material, which is available to authorized users.

DNA methylation signatures in peripheral blood strongly predict all-cause mortality

DNA methylation (DNAm) has been revealed to play a role in various diseases. Here we performed epigenome-wide screening and validation to identify mortality-related DNAm signatures in a general population-based cohort with up to 14 years follow-up. In the discovery panel in a case-cohort approach, 11,063 CpGs reach genome-wide significance (FDR<0.05). 58 CpGs, mapping to 38 well-known disease-related genes and 14 intergenic regions, are confirmed in a validation panel. A mortality risk score based on ten selected CpGs exhibits strong association with all-cause mortality, showing hazard ratios (95% CI) of 2.16 (1.10–4.24), 3.42 (1.81–6.46) and 7.36 (3.69–14.68), respectively, for participants with scores of 1, 2–5 and 5+ compared with a score of 0. These associations are confirmed in an independent cohort and are independent from the ‘epigenetic clock'. In conclusion, DNAm of multiple disease-related genes are strongly linked to mortality outcomes. The DNAm-based risk score might be informative for risk assessment and stratification.

DNA methylation is modulated by environmental factors and has a role in many complex diseases. Here, the authors find that methylation at specific DNA sites is associated with all-cause mortality, and a methylation-based risk score may be informative for risk assessment and stratification.

SQLE induces epithelial-to-mesenchymal transition by regulating of miR-133b in esophageal squamous cell carcinoma

Increasing evidence suggests that microRNAs, which control gene expression at the post-transcriptional level, are aberrantly expressed in cancers and play significant roles in carcinogenesis and cancer progression. In this study, we show differential miR-133b down-expression in human esophageal squamous cell carcinoma (ESCC) cells and tissues. In addition, squalene epoxidase (SQLE), a key enzyme of cholesterol synthesis, is identified as the direct downstream target gene of miR-133b by luciferase gene reporter assay. Furthermore, ectogenic miR-133b expression and SQLE knockdown can inhibit proliferation, invasion, and metastasis, and diminish epithelial-to-mesenchymal transition (EMT) traits of ESCC in vitro, implying that miR-133b-dependent SQLE can induce tumorigenicity and that SQLE is an EMT inducer. Xenograft experiment results also proved the biological function of SQLE in vivo. Therefore, we conclude that miR-133b-dependent SQLE plays a critical role in the potential metastasis mechanisms in ESCC.

Language-Agnostic Reproducible Data Analysis Using Literate Programming

A modern biomedical research project can easily contain hundreds of analysis steps and lack of reproducibility of the analyses has been recognized as a severe issue. While thorough documentation enables reproducibility, the number of analysis programs used can be so large that in reality reproducibility cannot be easily achieved. Literate programming is an approach to present computer programs to human readers. The code is rearranged to follow the logic of the program, and to explain that logic in a natural language. The code executed by the computer is extracted from the literate source code. As such, literate programming is an ideal formalism for systematizing analysis steps in biomedical research. We have developed the reproducible computing tool Lir (literate, reproducible computing) that allows a tool-agnostic approach to biomedical data analysis. We demonstrate the utility of Lir by applying it to a case study. Our aim was to investigate the role of endosomal trafficking regulators to the progression of breast cancer. In this analysis, a variety of tools were combined to interpret the available data: a relational database, standard command-line tools, and a statistical computing environment. The analysis revealed that the lipid transport related genes LAPTM4B and NDRG1 are coamplified in breast cancer patients, and identified genes potentially cooperating with LAPTM4B in breast cancer progression. Our case study demonstrates that with Lir, an array of tools can be combined in the same data analysis to improve efficiency, reproducibility, and ease of understanding. Lir is an open-source software available at github.com/borisvassilev/lir.

Genomic alterations underlie a pan-cancer metabolic shift associated with tumour hypoxia

BACKGROUND

Altered metabolism is a hallmark of cancer. However, the role of genomic changes in metabolic genes driving the tumour metabolic shift remains to be elucidated. Here, we have investigated the genomic and transcriptomic changes underlying this shift across ten different cancer types.

RESULTS

A systematic pan-cancer analysis of 6538 tumour/normal samples covering ten major cancer types identified a core metabolic signature of 44 genes that exhibit high frequency somatic copy number gains/amplifications (>20 % cases) associated with increased mRNA expression (ρ > 0.3, q < 10(-3)). Prognostic classifiers using these genes were confirmed in independent datasets for breast and kidney cancers. Interestingly, this signature is strongly associated with hypoxia, with nine out of ten cancer types showing increased expression and five out of ten cancer types showing increased gain/amplification of these genes in hypoxic tumours (P ≤ 0.01). Further validation in breast and colorectal cancer cell lines highlighted squalene epoxidase, an oxygen-requiring enzyme in cholesterol biosynthesis, as a driver of dysregulated metabolism and a key player in maintaining cell survival under hypoxia.

CONCLUSIONS

This study reveals somatic genomic alterations underlying a pan-cancer metabolic shift and suggests genomic adaptation of these genes as a survival mechanism in hypoxic tumours.

Cholesterol biosynthesis pathway as a novel mechanism of resistance to estrogen deprivation in estrogen receptor-positive breast cancer

Background

Therapies targeting estrogenic stimulation in estrogen receptor-positive (ER+) breast cancer (BC) reduce mortality, but resistance remains a major clinical problem. Molecular studies have shown few high-frequency mutations to be associated with endocrine resistance. In contrast, expression profiling of primary ER+ BC samples has identified several promising signatures/networks for targeting.

Methods

To identify common adaptive mechanisms associated with resistance to aromatase inhibitors (AIs), we assessed changes in global gene expression during adaptation to long-term estrogen deprivation (LTED) in a panel of ER+ BC cell lines cultured in 2D on plastic (MCF7, T47D, HCC1428, SUM44 and ZR75.1) or in 3D on collagen (MCF7) to model the stromal compartment. Furthermore, dimethyl labelling followed by LC-MS/MS was used to assess global changes in protein abundance. The role of target genes/proteins on proliferation, ER-mediated transcription and recruitment of ER to target gene promoters was analysed.

Results

The cholesterol biosynthesis pathway was the common upregulated pathway in the ER+ LTED but not the ER– LTED cell lines, suggesting a potential mechanism dependent on continued ER expression. Targeting the individual genes of the cholesterol biosynthesis pathway with siRNAs caused a 30–50 % drop in proliferation. Further analysis showed increased expression of 25-hydroxycholesterol (HC) in the MCF7 LTED cells. Exogenous 25-HC or 27-HC increased ER-mediated transcription and expression of the endogenous estrogen-regulated gene TFF1 in ER+ LTED cells but not in the ER– LTED cells. Additionally, recruitment of the ER and CREB-binding protein (CBP) to the TFF1 and GREB1 promoters was increased upon treatment with 25-HC and 27-HC. In-silico analysis of two independent studies of primary ER+ BC patients treated with neoadjuvant AIs showed that increased expression of MSMO1, EBP, LBR and SQLE enzymes, required for cholesterol synthesis and increased in our in-vitro models, was significantly associated with poor response to endocrine therapy.

Conclusion

Taken together, these data provide support for the role of cholesterol biosynthesis enzymes and the cholesterol metabolites, 25-HC and 27-HC, in a novel mechanism of resistance to endocrine therapy in ER+ BC that has potential as a therapeutic target.

Electronic supplementary material

The online version of this article (doi:10.1186/s13058-016-0713-5) contains supplementary material, which is available to authorized users.

Squalene epoxidase is a bona fide oncogene by amplification with clinical relevance in breast cancer

SQLE encodes squalene epoxidase, a key enzyme in cholesterol synthesis. SQLE has sporadically been reported among copy-number driven transcripts in multi-omics cancer projects. Yet, its functional relevance has never been subjected to systematic analyses. Here, we assessed the correlation of SQLE copy number (CN) and gene expression (GE) across multiple cancer types, focusing on the clinico-pathological associations in breast cancer (BC). We then investigated whether any biological effect of SQLE inhibition could be observed in BC cell line models. Breast, ovarian, and colorectal cancers showed the highest CN driven GE among 8,783 cases from 22 cancer types, with BC presenting the strongest one. SQLE overexpression was more prevalent in aggressive BC, and was an independent prognostic factor of unfavorable outcome. Through SQLE pharmacological inhibition and silencing in a panel of BC cell lines portraying the diversity of SQLE CN and GE, we demonstrated that SQLE inhibition resulted in a copy-dosage correlated decrease in cell viability, and in a noticeable increase in replication time, only in lines with detectable SQLE transcript. Altogether, our results pinpoint SQLE as a bona fide metabolic oncogene by amplification, and as a therapeutic target in BC. These findings could have implications in other cancer types.