Majority of differentially expressed genes are down-regulated during malignant transformation in a four-stage model

Personalized differential expression analysis in triple-negative breast cancer

Identification of individual-level differentially expressed genes (DEGs) is a pre-step for the analysis of disease-specific biological mechanisms and precision medicine. Previous algorithms cannot balance accuracy and sufficient statistical power. Herein, RankCompV2, designed for identifying population-level DEGs based on relative expression orderings, was adjusted to identify individual-level DEGs. Furthermore, an optimized version of individual-level RankCompV2, named as RankCompV2.1, was designed based on the assumption that the rank positions of genes and relative rank differences of gene pairs would influence the identification of individual-level DEGs. In comparison to other individualized analysis algorithms, RankCompV2.1 performed better on statistical power, computational efficiency, and acquired coequal accuracy in both simulation and real paired cancer-normal data from ten cancer types. Besides, single sample GSEA and Gene Set Variation Analysis analysis showed that pathways enriched with up-regulated and down-regulated genes presented higher and lower enrichment scores, respectively. Furthermore, we identified 16 genes that were universally deregulated in 966 triple-negative breast cancer (TNBC) samples and interacted with Food and Drug Administration (FDA)-approved drugs or antineoplastic agents, indicating notable therapeutic targets for TNBC. In addition, we also identified genes with highly variable deregulation status and used these genes to cluster TNBC samples into three subgroups with different prognoses. The subgroup with the poorest outcome was characterized by down-regulated immune-regulated pathways, signal transduction pathways, and apoptosis-related pathways. Protein-protein interaction network analysis revealed that OAS family genes may be promising drug targets to activate tumor immunity in this subgroup. In conclusion, RankCompV2.1 is capable of identifying individual-level DEGs with high accuracy and statistical power, analyzing mechanisms of carcinogenesis and exploring therapeutic strategy.

Disentangling the splicing factor programs underlying complex molecular phenotypes

The regulation of exon inclusion through alternative splicing tunes the cell's behavior by increasing the functional diversity of the transcriptome and the proteome. Splicing factors work in concert to generate gene isoform pools that contribute to cell phenotypes yet their activity is controlled by multiple regulatory and signaling layers. This hinders identification of functional, phenotype-specific splicing factors using traditional single-omic measurements, such as their mutational state or expression. To address this challenge, we propose repurposing the virtual inference of protein activity by enriched regulon analysis (VIPER) to measure splicing factor activity solely from their downstream exon transcriptomic inclusion signatures. This approach is effective in assessing the effect of co-occurring splicing factor perturbations, as well as their post-translational regulation. As proof of concept, we dissect recurrent splicing factor programs underlying tumorigenesis including aberrantly activated factors acting as oncogenes and inactivated ones acting as tumor suppressors, which are undetectable by more conventional methodologies. Activation and inactivation of these cancer splicing programs effectively stratifies overall survival, as well as cancer hallmarks such as proliferation and immune evasion. Altogether, repurposing network-based inference of protein activity for splicing factor networks distills common, functionally relevant splicing programs in otherwise heterogeneous molecular contexts.

TRPS1 is a Highly Sensitive Marker for Breast Cancer: A Tissue Microarray Study Evaluating More Than 19,000 Tumors From 152 Different Tumor Entities

Trichorhinophalangeal syndrome 1 (TRPS1) is a nuclear protein highly expressed in breast epithelial cells. TRPS1 immunohistochemistry (IHC) has been suggested as a breast cancer marker. To determine the diagnostic and prognostic utility of TRPS1 IHC, tissue microarrays containing 19,201 samples from 152 different tumor types and subtypes were analyzed. GATA3 IHC was performed in a previous study. TRPS1 staining was seen in 86 of 152 tumor categories with 36 containing at least one strongly positive case. TRPS1 staining predominated in various types of breast carcinomas (51%-100%), soft tissue tumors (up to 100%), salivary gland tumors (up to 46%), squamous cell carcinomas (up to 35%), and gynecological cancers (up to 40%). TRPS1 positivity occurred in 1.8% of 1083 urothelial neoplasms. In invasive breast carcinoma of no special type, low TRPS1 expression was linked to high grade ( P = 0.0547), high pT ( P < 0.0001), nodal metastasis ( P = 0.0571), loss of estrogen receptor and progesterone receptor expression ( P < 0.0001 each), and triple-negative status ( P < 0.0001) but was unrelated to patient survival ( P = 0.8016). In squamous cell carcinomas from 11 different sites, low TRPS1 expression was unrelated to tumor phenotype. Positivity for both TRPS1 and GATA3 occurred in 47.4% to 100% of breast cancers, up to 30% of salivary gland tumors, and 29 (0.3%) of 9835 tumors from 134 other cancer entities. TRPS1 IHC has high utility for the identification of cancers of breast (or salivary gland) origin, especially in combination with GATA3. The virtual absence of TRPS1 positivity in urothelial neoplasms is useful for the distinction of GATA3-positive urothelial carcinoma from breast cancer.

The breast cancer coagulome in the tumor microenvironment and its role in prognosis and treatment response to chemotherapy

BACKGROUND

The procoagulant phenotype in cancer is linked to thrombosis, cancer progression, and immune response. A novel treatment that reduces the risk of both thrombosis and cancer progression without excess bleeding risk remains to be identified.

OBJECTIVES

Here, we aimed to broadly investigate the breast tumor coagulome and its relation to prognosis, treatment response to chemotherapy, and the tumor microenvironment.

METHODS

Key coagulation-related genes (n = 35) were studied in a Norwegian cohort with tumor (n = 134) and normal (n = 189) tissue and in the Cancer Genome Atlas (n = 1052) data set. We performed gene set variation analysis in the Norwegian cohort, and in the Cancer Genome Atlas cohort, associations with the tumor microenvironment and prognosis were evaluated. Analyses were performed with cBioPortal, Estimation of Stromal and Immune cells in Malignant Tumors Using Expression Data, Tumor Immune Estimation Resource, the integrated repository portal for tumor-immune system interactions, Tumor Immune Single-cell Hub 2, and the receiver operating characteristic plotter. Six independent breast cancer cohorts were used to study the tumor coagulome and treatment response to chemotherapy.

RESULTS

Twenty-two differentially expressed coagulation-related genes were identified in breast tumors. Several coagulome factors were correlated with tumor microenvironment characteristics and were expressed by nonmalignant cells in the tumor microenvironment. PLAT and F8 were independent predictors of better overall survival and progression-free survival, respectively. F12 and PLAU were predictors of worse progression-free survival. The PROCR-THBD-PLAT signature showed a promising predictive value (area under the curve, 0.75; 95% CI, 0.69-0.81; P = 3.6 × 10-17) for combination chemotherapy with fluorouracil, epirubicin, and cyclophosphamide.

CONCLUSION

The breast tumor coagulome showed potential in prediction of prognosis and chemotherapy response. Cells within the tumor microenvironment are sources of coagulome factors and may serve as therapeutic targets of coagulation factors.

Roles of Histone H2A Variants in Cancer Development, Prognosis, and Treatment

Histones are nuclear proteins essential for packaging genomic DNA and epigenetic gene regulation. Paralogs that can substitute core histones (H2A, H2B, H3, and H4), named histone variants, are constitutively expressed in a replication-independent manner throughout the cell cycle. With specific chaperones, they can be incorporated to chromatin to modify nucleosome stability by modulating interactions with nucleosomal DNA. This allows the regulation of essential fundamental cellular processes for instance, DNA damage repair, chromosomal segregation, and transcriptional regulation. Among all the histone families, histone H2A family has the largest number of histone variants reported to date. Each H2A variant has multiple functions apart from their primary role and some, even be further specialized to perform additional tasks in distinct lineages, such as testis specific shortH2A (sH2A). In the past decades, the discoveries of genetic alterations and mutations in genes encoding H2A variants in cancer had revealed variants' potentiality in driving carcinogenesis. In addition, there is growing evidence that H2A variants may act as novel prognostic indicators or biomarkers for both early cancer detection and therapeutic treatments. Nevertheless, no studies have ever concluded all identified variants in a single report. Here, in this review, we summarize the respective functions for all the 19 mammalian H2A variants and their roles in cancer biology whilst potentiality being used in clinical setting.

Loss of Upk1a and Upk1b expression is linked to stage progression in urothelial carcinoma of the bladder

BACKGROUND

Uroplakin-1a (Upk1a) and uroplakin-1b (Upk1b) have recently been identified as diagnostic markers for the distinction of urothelial carcinomas from other solid tumor entities. Both proteins play an important role in the stabilization and strengthening of epithelial cells that line the bladder.

METHODS

To evaluate the prognostic role of uroplakin expression in urothelial carcinomas, more than 2700 urothelial neoplasms were analyzed in a tissue microarray format by immunohistochemistry. To further assess the diagnostic role of uroplakin immunohistochemistry, results were compared with preexisting GATA3 data.

RESULT

The fraction of Upk1a/Upk1b positive cases decreased slightly from pTaG2 low-grade (88% positive for Upk1a/87% positive for Upk1b) and pTaG2 high-grade (92%/89%) to pTaG3 (83%/88%; p > 0.05) and was lower in muscle-invasive (pT2-4) carcinomas (42%/64%; p < 0.0001/p < 0.0001 for pTa vs. pT2-4). Within pT2-4 carcinomas, high expression of Upk1a and Upk1b was linked to nodal metastasis and lymphatic vessel infiltration (p < 0.05) but unrelated to patient outcome. There were significant associations between Upk1a, Upk1b and GATA3 immunostaining (p < 0.0001 each), but 11% of GATA3 negative cancers were Upk1a/b positive and 8% of Upk1a/b negative cancers were GATA3 positive. Absence of GATA3/Upk1a/b staining was significantly linked to poor patient survival in the subgroup of 126 pT4 carcinomas (p = 0.0004) but not in pT2 and pT3 cancers.

CONCLUSIONS

In summary, the results of our study demonstrate that Upk1a and/or Upk1b immunohistochemistry can complement GATA3 for the distinction of urothelial carcinomas. Furthermore, a progressive loss of Upk1a/b expression during stage progression and a prognostic role of the combination GATA3/Upk1a/Upk1b in pT4 carcinomas is evident.

Multi-omics integration analysis unveils heterogeneity in breast cancer at the individual level

Identifying robust breast cancer subtypes will help to reveal the cancer heterogeneity. However, previous breast cancer subtypes were based on population-level quantitative gene expression, which is affected by batch effects and cannot be applied to individuals. We detected differential gene expression, genomic, and epigenomic alterations to identify driver differential expression at the individual level. The individual driver differential expression reflected the breast cancer patients' heterogeneity and revealed four subtypes. Mesenchymal subtype as the most aggressive subtype harbored deletion and downregulated expression of genes in chromosome 11q23 region. Specifically, silencing of the SDHD gene in 11q23 promoted the invasion and migration of breast cancer cells in vitro by the epithelial-mesenchymal transition. The immunologically hot subtype displayed an immune-hot microenvironment, including high T-cell infiltration and upregulated PD-1 and CTLA4. Luminal and genomic-unstable subtypes showed opposite macrophage polarization, which may be regulated by the ligand-receptor pairs of CD99. The integration of multi-omics data at the individual level provides a powerful framework for elucidating the heterogeneity of breast cancer.

Regulation of cellular contractile force, shape and migration of fibroblasts by oncogenes and Histone deacetylase 6

The capacity of cells to adhere to, exert forces upon and migrate through their surrounding environment governs tissue regeneration and cancer metastasis. The role of the physical contractile forces that cells exert in this process, and the underlying molecular mechanisms are not fully understood. We, therefore, aimed to clarify if the extracellular forces that cells exert on their environment and/or the intracellular forces that deform the cell nucleus, and the link between these forces, are defective in transformed and invasive fibroblasts, and to indicate the underlying molecular mechanism of control. Confocal, Epifluorescence and Traction force microscopy, followed by computational analysis, showed an increased maximum contractile force that cells apply on their environment and a decreased intracellular force on the cell nucleus in the invasive fibroblasts, as compared to normal control cells. Loss of HDAC6 activity by tubacin-treatment and siRNA-mediated HDAC6 knockdown also reversed the reduced size and more circular shape and defective migration of the transformed and invasive cells to normal. However, only tubacin-mediated, and not siRNA knockdown reversed the increased force of the invasive cells on their surrounding environment to normal, with no effects on nuclear forces. We observed that the forces on the environment and the nucleus were weakly positively correlated, with the exception of HDAC6 siRNA-treated cells, in which the correlation was weakly negative. The transformed and invasive fibroblasts showed an increased number and smaller cell-matrix adhesions than control, and neither tubacin-treatment, nor HDAC6 knockdown reversed this phenotype to normal, but instead increased it further. This highlights the possibility that the control of contractile force requires separate functions of HDAC6, than the control of cell adhesions, spreading and shape. These data are consistent with the possibility that defective force-transduction from the extracellular environment to the nucleus contributes to metastasis, via a mechanism that depends upon HDAC6. To our knowledge, our findings present the first correlation between the cellular forces that deforms the surrounding environment and the nucleus in fibroblasts, and it expands our understanding of how cells generate contractile forces that contribute to cell invasion and metastasis.

ALD-R491 regulates vimentin filament stability and solubility, cell contractile force, cell migration speed and directionality

Metastasizing cells express the intermediate filament protein vimentin, which is used to diagnose invasive tumors in the clinic. However, the role of vimentin in cell motility, and if the assembly of non-filamentous variants of vimentin into filaments regulates cell migration remains unclear. We observed that the vimentin-targeting drug ALD-R491 increased the stability of vimentin filaments, by reducing filament assembly and/or disassembly. ALD-R491-treatment also resulted in more bundled and disorganized filaments and an increased pool of non-filamentous vimentin. This was accompanied by a reduction in size of cell-matrix adhesions and increased cellular contractile forces. Moreover, during cell migration, cells showed erratic formation of lamellipodia at the cell periphery, loss of coordinated cell movement, reduced cell migration speed, directionality and an elongated cell shape with long thin extensions at the rear that often detached. Taken together, these results indicate that the stability of vimentin filaments and the soluble pool of vimentin regulate the speed and directionality of cell migration and the capacity of cells to migrate in a mechanically cohesive manner. These observations suggest that the stability of vimentin filaments governs the adhesive, physical and migratory properties of cells, and expands our understanding of vimentin functions in health and disease, including cancer metastasis.

Hyaluronan nanoscale clustering and Hyaluronan synthase 2 expression are linked to the invasion of child fibroblasts and infantile fibrosarcoma in vitro and in vivo

Infantile fibrosarcoma is a rare childhood tumour that originates in the fibrous connective tissue of the long bones for which there is an urgent need to identify novel therapeutic targets. This study aims to clarify the role of the extracellular matrix component hyaluronan in the invasion of child fibroblasts and Infantile fibrosarcoma into the surrounding environment. Using nanoscale super-resolution STED (Stimulated emission depletion) microscopy followed by computational image analysis, we observed, for the first time, that invasive child fibroblasts showed increased nanoscale clustering of hyaluronan at the cell periphery, as compared to control cells. Hyaluronan was not observed within focal adhesions. Bioinformatic analyses further revealed that the increased nanoscale hyaluronan clustering was accompanied by increased gene expression of Hyaluronan synthase 2, reduced expression of Hyaluronidase 2 and CD44, and no change of Hyaluronan synthase 1 and Hyaluronidases 1, 3, 4 or 5. We further observed that the expression of the Hyaluronan synthase 1, 2 and 3, and the Hyaluronidase 3 and 5 genes was linked to reduced life expectancy of fibrosarcoma patients. The invasive front of infantile fibrosarcoma tumours further showed increased levels of hyaluronan, as compared to the tumour centre. Taken together, our findings are consistent with the possibility that while Hyaluronan synthase 2 increases the levels, the Hyaluronidases 3 and 5 reduce the weight of hyaluronan, resulting in the nanoscale clustering of hyaluronan at the leading edge of cells, cell invasion and the spread of Infantile fibrosarcoma.

Development and validation of a prognostic 9-gene signature for colorectal cancer

Introduction

Colorectal cancer (CRC) is one of the most prevalent cancers globally with a high mortality rate. Predicting prognosis using disease progression and cancer pathologic stage is insufficient, and a prognostic factor that can accurately evaluate patient prognosis needs to be developed. In this study, we aimed to infer a prognostic gene signature to identify a functional signature associated with the prognosis of CRC patients.

Methods

First, we used univariate Cox regression, least absolute shrinkage and selection operator (lasso) regression, and multivariate Cox regression analyses to screen genes significantly associated with CRC patient prognosis, from colorectal cancer RNA sequencing data in The Cancer Genome Atlas (TCGA) database. We then calculated the risk score (RS) for each patient based on the expression of the nine candidate genes and developed a prognostic signature.

Results

Based on the optimal cut-off on the receiver operating characteristic (ROC) curve, patients were separated into high- and low-risk groups, and the difference in overall survival between the two groups was examined. Patients in the low-risk group had a better overall survival rate than those in the high-risk group. The results were validated using the GSE72970, GSE39582, and GSE17536 Gene Expression Omnibus (GEO) datasets, and the same conclusions were reached. ROC curve test of the RS signature also indicated that it had excellent accuracy. The RS signature was then compared with traditional clinical factors as a prognostic indicator, and we discovered that the RS signature had superior predictive ability.

Conclusion

The RS signature developed in this study has excellent predictive power for the prognosis of patients with CRC and broad applicability as a prognostic indicator for patients.

Metastasising Fibroblasts Show an HDAC6-Dependent Increase in Migration Speed and Loss of Directionality Linked to Major Changes in the Vimentin Interactome

Metastasising cells express the intermediate filament protein vimentin, which is used to diagnose invasive tumours in the clinic. We aimed to clarify how vimentin regulates the motility of metastasising fibroblasts. STED super-resolution microscopy, live-cell imaging and quantitative proteomics revealed that oncogene-expressing and metastasising fibroblasts show a less-elongated cell shape, reduced cell spreading, increased cell migration speed, reduced directionality, and stronger coupling between these migration parameters compared to normal control cells. In total, we identified and compared 555 proteins in the vimentin interactome. In metastasising cells, the levels of keratin 18 and Rab5C were increased, while those of actin and collagen were decreased. Inhibition of HDAC6 reversed the shape, spreading and migration phenotypes of metastasising cells back to normal. Inhibition of HDAC6 also decreased the levels of talin 1, tropomyosin, Rab GDI β, collagen and emilin 1 in the vimentin interactome, and partially reversed the nanoscale vimentin organisation in oncogene-expressing cells. These findings describe the changes in the vimentin interactome and nanoscale distribution that accompany the defective cell shape, spreading and migration of metastasising cells. These results support the hypothesis that oncogenes can act through HDAC6 to regulate the vimentin binding of the cytoskeletal and cell–extracellular matrix adhesion components that contribute to the defective motility of metastasising cells.

Comparative transcriptome analysis between patient and endometrial cancer cell lines to determine common signaling pathways and markers linked to cancer progression

The rising incidence and mortality of endometrial cancer (EC) in the United States calls for an improved understanding of the disease's progression. Current methodologies for diagnosis and treatment rely on the use of cell lines as models for tumor biology. However, due to inherent heterogeneity and differential growing environments between cell lines and tumors, these comparative studies have found little parallels in molecular signatures. As a consequence, the development and discovery of preclinical models and reliable drug targets are delayed. In this study, we established transcriptome parallels between cell lines and tumors from The Cancer Genome Atlas (TCGA) with the use of optimized normalization methods. We identified genes and signaling pathways associated with regulating the transformation and progression of EC. Specifically, the LXR/RXR activation, neuroprotective role for THOP1 in Alzheimer's disease, and glutamate receptor signaling pathways were observed to be mostly downregulated in advanced cancer stage. While some of these highlighted markers and signaling pathways are commonly found in the central nervous system (CNS), our results suggest a novel function of these genes in the periphery. Finally, our study underscores the value of implementing appropriate normalization methods in comparative studies to improve the identification of accurate and reliable markers.

Transcriptome Analysis Identifies Doublesex and Mab-3 Related Transcription Factor (DMRT3) in Nasal Polyp Epithelial Cells of Patients Suffering from Non-Steroidal Anti-Inflammatory Drug-Exacerbated Respiratory Disease (AERD)

Background: Aspirin-exacerbated respiratory disease (AERD) is a syndrome characterised by chronic rhinosinusitis, nasal polyps, asthma and aspirin intolerance. An imbalance of eicosanoid metabolism with anover-production of cysteinyl leukotrienes (CysLTs) has been associated with AERD. However, the precise mechanisms underlying AERD are unknown. Objective: To establish the transcriptome of the nasal polyp airway epithelial cells derived from AERD patients to discover gene expression patterns in this disease. Methods: Nasal airway epithelial cells were isolated from 12 AERD polyps and 8 AERD non-polyp nasal mucosa samples as controls from the same subjects. Utilising the Illumina HiSeq 2500 platform, RNA samples were sequenced. Potential gene candidate DMRT3 was selected from the differentially-expressed genes for validation. Results: Comparative transcriptome profiling of nasal epithelial cells was accomplished in AERD. A total of 20 genes had twofold mean regulation expression differences or greater. In addition, 8 genes were upregulated, including doublesex and mab-3 related transcription factor 3 (DMRT3), and 12 genes were downregulated. Differentially regulated genes comprised roles in inflammation, defence and immunity. Metabolic process and embryonic development pathways were significantly enriched. Enzyme-linked immune sorbent assay (ELISA) results of DMRT3 in AERD patients were significantly upregulated compared to controls (p = 0.03). Immunohistochemistry (IHC) of AERD nasal polyps localised DMRT3 and was predominantly released in the airway epithelia. Conclusion: Findings suggest that DMRT3 could be potentially involved in nasal polyp development in AERD patients. Furthermore, several genes are downregulated, hinting at the dedifferentiation phenomenon in AERD polyps. However, further studies are imperative to confirm the exact mechanism of polyp formation in AERD patients.

Galaxy and MEAN Stack to Create a User-Friendly Workflow for the Rational Optimization of Cancer Chemotherapy

One aspect of personalized medicine is aiming at identifying specific targets for therapy considering the gene expression profile of each patient individually. The real-world implementation of this approach is better achieved by user-friendly bioinformatics systems for healthcare professionals. In this report, we present an online platform that endows users with an interface designed using MEAN stack supported by a Galaxy pipeline. This pipeline targets connection hubs in the subnetworks formed by the interactions between the proteins of genes that are up-regulated in tumors. This strategy has been proved to be suitable for the inhibition of tumor growth and metastasis in vitro. Therefore, Perl and Python scripts were enclosed in Galaxy for translating RNA-seq data into protein targets suitable for the chemotherapy of solid tumors. Consequently, we validated the process of target diagnosis by (i) reference to subnetwork entropy, (ii) the critical value of density probability of differential gene expression, and (iii) the inhibition of the most relevant targets according to TCGA and GDC data. Finally, the most relevant targets identified by the pipeline are stored in MongoDB and can be accessed through the aforementioned internet portal designed to be compatible with mobile or small devices through Angular libraries.

Anillin is an emerging regulator of tumorigenesis, acting as a cortical cytoskeletal scaffold and a nuclear modulator of cancer cell differentiation

Remodeling of the intracellular cytoskeleton plays a key role in accelerating tumor growth and metastasis. Targeting different cytoskeletal elements is important for existing and future anticancer therapies. Anillin is a unique scaffolding protein that interacts with major cytoskeletal structures, e.g., actin filaments, microtubules and septin polymers. A well-studied function of this scaffolding protein is the regulation of cytokinesis at the completion of cell division. Emerging evidence suggest that anillin has other important activities in non-dividing cells, including control of intercellular adhesions and cell motility. Anillin is markedly overexpressed in different solid cancers and its high expression is commonly associated with poor prognosis of patient survival. This review article summarizes rapidly accumulating evidence that implicates anillin in the regulation of tumor growth and metastasis. We focus on molecular and cellular mechanisms of anillin-dependent tumorigenesis that include both canonical control of cytokinesis and novel poorly understood functions as a nuclear regulator of the transcriptional reprogramming and phenotypic plasticity of cancer cells.

Computational identification of biomarker genes for lung cancer considering treatment and non-treatment studies

BACKGROUND

Lung cancer is the number one cancer killer in the world with more than 142,670 deaths estimated in the United States alone in the year 2019. Consequently, there is an overreaching need to identify the key biomarkers for lung cancer. The aim of this study is to computationally identify biomarker genes for lung cancer that can aid in its diagnosis and treatment. The gene expression profiles of two different types of studies, namely non-treatment and treatment, are considered for discovering biomarker genes. In non-treatment studies healthy samples are control and cancer samples are cases. Whereas, in treatment studies, controls are cancer cell lines without treatment and cases are cancer cell lines with treatment.

RESULTS

The Differentially Expressed Genes (DEGs) for lung cancer were isolated from Gene Expression Omnibus (GEO) database using R software tool GEO2R. A total of 407 DEGs (254 upregulated and 153 downregulated) from non-treatment studies and 547 DEGs (133 upregulated and 414 downregulated) from treatment studies were isolated. Two Cytoscape apps, namely, CytoHubba and MCODE, were used for identifying biomarker genes from functional networks developed using DEG genes. This study discovered two distinct sets of biomarker genes - one from non-treatment studies and the other from treatment studies, each set containing 16 genes. Survival analysis results show that most non-treatment biomarker genes have prognostic capability by indicating low-expression groups have higher chance of survival compare to high-expression groups. Whereas, most treatment biomarkers have prognostic capability by indicating high-expression groups have higher chance of survival compare to low-expression groups.

CONCLUSION

A computational framework is developed to identify biomarker genes for lung cancer using gene expression profiles. Two different types of studies - non-treatment and treatment - are considered for experiment. Most of the biomarker genes from non-treatment studies are part of mitosis and play vital role in DNA repair and cell-cycle regulation. Whereas, most of the biomarker genes from treatment studies are associated to ubiquitination and cellular response to stress. This study discovered a list of biomarkers, which would help experimental scientists to design a lab experiment for further exploration of detail dynamics of lung cancer development.

Size-Exclusion Chromatography as a Technique for the Investigation of Novel Extracellular Vesicles in Cancer

Simple Summary

Extracellular vesicles (EVs) are small particles that are released by cancer cells, and they may hold vital information for researchers looking for early markers for diagnosis. Size-exclusion chromatography (SEC) is a classical technique that has become increasingly popular and can be used for rapid isolation and investigation of both their cargo and functionality. This systematic review highlights its main technical aspects, the type of materials involved and by covering the findings of the identified papers hopes to demonstrate the utility of this method in cancer research to date.

Abstract

Cancer cells release extracellular vesicles, which are a rich target for biomarker discovery and provide a promising mechanism for liquid biopsy. Size-exclusion chromatography (SEC) is an increasingly popular technique, which has been rediscovered for the purposes of extracellular vesicle (EV) isolation and purification from diverse biofluids. A systematic review was undertaken to identify all papers that described size exclusion as their primary EV isolation method in cancer research. In all, 37 papers were identified and discussed, which showcases the breadth of applications in which EVs can be utilised, from proteomics, to RNA, and through to functionality. A range of different methods are highlighted, with Sepharose-based techniques predominating. EVs isolated using SEC are able to identify cancer cells, highlight active pathways in tumourigenesis, clinically distinguish cohorts, and remain functionally active for further experiments.

Transcriptome profiling of Ewing sarcomas - treatment resistance pathways and IGF-dependency

Ewing sarcomas (ESs) are aggressive sarcomas driven by EWS fusion genes. We sought to investigate whether whole-transcriptome sequencing (RNA-seq) could be used to detect patterns associated with chemotherapy response or tumor progression after first-line treatment. Transcriptome sequencing (RNA-seq) of 13 ES cases was performed. Among the differentially expressed pathways, we identified IGF2 expression as a potential driver of chemotherapy response and progression. We investigated the effect of IGF2 on proliferation, radioresistance, apoptosis, and the transcriptome pattern in four ES cell lines and the effect of IGF2 expression in a validation series of 14 patients. Transcriptome analysis identified differentially expressed genes (adj. P < 0.005) and pathways associated with chemotherapy response (285 genes), short overall survival (662 genes), and progression after treatment (447 genes). Imprinting independent promoter P3-mediated IGF2 expression was identified in a subset of cases with aggressive clinical course. In ES cell lines, IGF2 induced proliferation, but promoted radioresistance only in CADO cells. High IGF2 expression was also significantly associated with shorter overall survival in patients with ES. Transcriptome analysis of the clinical samples and the cell lines revealed an IGF-dependent signature, potentially related to a stem cell-like phenotype. Transcriptome analysis is a potentially powerful complementary tool to predict the clinical behavior of ES and may be utilized for clinical trial stratification strategies and personalized oncology. Certain gene signatures, for example, IGF-related pathways, are coupled to biological functions that could be of clinical importance. Finally, our results indicate that IGF inhibition may be successful as a first-line therapy in conjunction with conventional radiochemotherapy for a subset of patients.

Edgetic perturbation signatures represent known and novel cancer biomarkers

Isoform switching is a recently characterized hallmark of cancer, and often translates to the loss or gain of domains mediating protein interactions and thus, the re-wiring of the interactome. Recent computational tools leverage domain-domain interaction data to resolve the condition-specific interaction networks from RNA-Seq data accounting for the domain content of the primary transcripts expressed. Here, we used The Cancer Genome Atlas RNA-Seq datasets to generate 642 patient-specific pairs of interactomes corresponding to both the tumor and the healthy tissues across 13 cancer types. The comparison of these interactomes provided a list of patient-specific edgetic perturbations of the interactomes associated with the cancerous state. We found that among the identified perturbations, select sets are robustly shared between patients at the multi-cancer, cancer-specific and cancer sub-type specific levels. Interestingly, the majority of the alterations do not directly involve significantly mutated genes, nevertheless, they strongly correlate with patient survival. The findings (available at EdgeExplorer: “http://webclu.bio.wzw.tum.de/EdgeExplorer”) are a new source of potential biomarkers for classifying cancer types and the proteins we identified are potential anti-cancer therapy targets.

Fungal canker pathogens trigger carbon starvation by inhibiting carbon metabolism in poplar stems

Carbon starvation is the current leading hypothesis of plant mortality mechanisms under drought stress; recently, it is also used to explain tree die-off in plant diseases. However, the molecular biology of the carbon starvation pathway is unclear. Here, using a punch inoculation system, we conducted transcriptome and physiological assays to investigate pathogen response in poplar stems at the early stages of Botryosphaeria and Valsa canker diseases. Transcriptome assays showed that the majority of differentially expressed genes (DEGs) in stem phloem and xylem, such as genes involved in carbon metabolism and transportation, aquaporin genes (in xylem) and genes related to the biosynthesis of secondary metabolites and the phenylpropanoid pathway (related to lignin synthesis), were downregulated at 7 days after inoculation (DAI). Results also showed that the expression of the majority of disease-resistance genes upregulated in poplar stems, which may be connected with the downregulation expression of the majority of WRKY family genes. Physiological assays showed that transpiration rate decreased but WUE (water use efficiency) increased the 3 and 7 DAI, while the net photosynthetic rate decreased at 11 DAI in Botryosphaeria infected poplars (ANOVA, P < 0.05). The NSC (non-structural carbohydrates) content assays showed that the soluble sugar content of stem phloem samples increased at 3, 7, and 11 DAI that might due to the impede of pathogen infection. However, soluble sugar content of stem xylem and root samples decreased at 11 DAI; in contrast, the starch content unchanged. Therefore, results revealed a chronological order of carbon related molecular and physiological performance: declination of genes involved in carbon and starch metabolism first (at least at 7 DAI), declination of assimilation and carbon reserve (at 11 DAI) second. Results implied a potential mechanism that affects the host carbon reserve, by directly inhibiting the expression of genes involved in carbon metabolism and transport.

Pan-cancer analysis identifies telomerase-associated signatures and cancer subtypes

Background

Cancer cells become immortalized through telomere maintenance mechanisms, such as telomerase reverse transcriptase (TERT) activation. In addition to maintaining telomere length, TERT activates manifold cell survival signaling pathways. However, telomerase-associated gene signatures in cancer remain elusive.

Methods

We performed a systematic analysis of TERT high (TERT^high) and low (TERT^low) cancers using multidimensional data from The Cancer Genome Atlas (TCGA). Multidimensional data were analyzed by propensity score matching weight algorithm. Coexpression networks were constructed by weight gene coexpression network analysis (WGCNA). Random forest classifiers were generated to identify cancer subtypes.

Results

The TERT^high-specific mRNA expression signature is associated with cell cycle-related coexpression modules across cancer types. Experimental screening of hub genes in the cell cycle module suggested TPX2 and EXO1 as potential regulators of telomerase activity and cell survival. MiRNA analysis revealed that the TERT^high-specific miR-17-92 cluster can target biological processes enriched in TERT^low cancer and that its expression is negatively correlated with the tumor/normal telomere length ratio. Intriguingly, TERT^high cancers tend to have mutations in extracellular matrix organization genes and amplify MAPK signaling. By mining the clinical actionable gene database, we uncovered a number of TERT^high-specific somatic mutations, amplifications and high expression genes containing therapeutic targets. Finally, a random forest classifier integrating telomerase-associated multi-omics signatures identifies two cancer subtypes showed profound differences in telomerase activity and patient survival.

Conclusions

In summary, our results depict a telomerase-associated molecular landscape in cancers and provide therapeutic opportunities for cancer treatment.

Electronic supplementary material

The online version of this article (10.1186/s12943-019-1035-x) contains supplementary material, which is available to authorized users.

mRNA and miRNA expression profiling of follicular variant of papillary thyroid carcinoma with and without distant metastases

Follicular Variant of Papillary Thyroid Carcinoma (FVPTC) is usually associated with a good outcome. Nevertheless, in rare cases, it develops distant metastases (1-9%). Our goal was to investigate whether mRNA and miRNA expression profiles may help distinguish between metastatic versus non-metastatic FVPTCs. Twenty-four primary FVPTCs, 12 metastatic and 12 non-metastatic, with similar clinicopathological features were selected and analyzed by nanoString nCounter technology using two distinct panels for expression analysis of 740 mRNA and 798 miRNAs. Data analysis was performed using the nanoString nSolver 3.0 software. Forty-seven mRNA and 35 miRNAs were differentially expressed between the two groups. Using these mRNA and miRNAs, metastatic and non-metastatic FVPTCs were clearly divided into two distinct clusters. Our results indicate that FVPTCs with metastatic abilities have different expression profiles compared to the non-metastatic. A prospective validation is needed to evaluate the usefulness of this molecular approach in the early identification of high-risk FVPTCs.

Decreased expression of cell adhesion genes in cancer stem-like cells isolated from primary oral squamous cell carcinomas

The goal of this study was to isolate cancer stem-like cells marked by high expression of CD44, a putative cancer stem cell marker, from primary oral squamous cell carcinomas and identify distinctive gene expression patterns in these cells. From 1 October 2013 to 4 September 2015, 76 stage III-IV primary oral squamous cell carcinoma of the gingivobuccal sulcus were resected. In all, 13 tumours were analysed by immunohistochemistry to visualise CD44-expressing cells. Expression of CD44 within The Cancer Genome Atlas-Head and Neck Squamous Cell Carcinoma RNA-sequencing data was also assessed. Seventy resected tumours were dissociated into single cells and stained with antibodies to CD44 as well as CD45 and CD31 (together referred as Lineage/Lin). From 45 of these, CD44+Lin- and CD44-Lin- subpopulations were successfully isolated using fluorescence-activated cell sorting, and good-quality RNA was obtained from 14 such sorted pairs. Libraries from five pairs were sequenced and the results analysed using bioinformatics tools. Reverse transcription quantitative polymerase chain reaction was performed to experimentally validate the differential expression of selected candidate genes identified from the transcriptome sequencing in the same 5 and an additional 9 tumours. CD44 was expressed on the surface of poorly differentiated tumour cells, and within the The Cancer Genome Atlas-Head and Neck Squamous Cell Carcinoma samples, its messenger RNA levels were higher in tumours compared to normal. Transcriptomics revealed that 102 genes were upregulated and 85 genes were downregulated in CD44+Lin- compared to CD44-Lin- cells in at least 3 of the 5 tumours sequenced. The upregulated genes included those involved in immune regulation, while the downregulated genes were enriched for genes involved in cell adhesion. Decreased expression of PCDH18, MGP, SPARCL1 and KRTDAP was confirmed by reverse transcription quantitative polymerase chain reaction. Lower expression of the cell-cell adhesion molecule PCDH18 correlated with poorer overall survival in the The Cancer Genome Atlas-Head and Neck Squamous Cell Carcinoma data highlighting it as a potential negative prognostic factor in this cancer.

Transcriptome profiling of the interconnection of pathways involved in malignant transformation and response to hypoxia

In tumor tissues, hypoxia is a commonly observed feature resulting from rapidly proliferating cancer cells outgrowing their surrounding vasculature network. Transformed cancer cells are known to exhibit phenotypic alterations, enabling continuous proliferation despite a limited oxygen supply. The four-step isogenic BJ cell model enables studies of defined steps of tumorigenesis: the normal, immortalized, transformed, and metastasizing stages. By transcriptome profiling under atmospheric and moderate hypoxic (3% O₂) conditions, we observed that despite being highly similar, the four cell lines of the BJ model responded strikingly different to hypoxia. Besides corroborating many of the known responses to hypoxia, we demonstrate that the transcriptome adaptation to moderate hypoxia resembles the process of malignant transformation. The transformed cells displayed a distinct capability of metabolic switching, reflected in reversed gene expression patterns for several genes involved in oxidative phosphorylation and glycolytic pathways. By profiling the stage-specific responses to hypoxia, we identified ASS1 as a potential prognostic marker in hypoxic tumors. This study demonstrates the usefulness of the BJ cell model for highlighting the interconnection of pathways involved in malignant transformation and hypoxic response.

GeneGini: Assessment via the Gini Coefficient of Reference "Housekeeping" Genes and Diverse Human Transporter Expression Profiles

The expression levels of SLC or ABC membrane transporter transcripts typically differ 100- to 10,000-fold between different tissues. The Gini coefficient characterizes such inequalities and here is used to describe the distribution of the expression of each transporter among different human tissues and cell lines. Many transporters exhibit extremely high Gini coefficients even for common substrates, indicating considerable specialization consistent with divergent evolution. The expression profiles of SLC transporters in different cell lines behave similarly, although Gini coefficients for ABC transporters tend to be larger in cell lines than in tissues, implying selection. Transporter genes are significantly more heterogeneously expressed than the members of most non-transporter gene classes. Transcripts with the stablest expression have a low Gini index and often differ significantly from the "housekeeping" genes commonly used for normalization in transcriptomics/qPCR studies. PCBP1 has a low Gini coefficient, is reasonably expressed, and is an excellent novel reference gene. The approach, referred to as GeneGini, provides rapid and simple characterization of expression-profile distributions and improved normalization of genome-wide expression-profiling data.

Epigenomic analysis in a cell-based model reveals the roles of H3K9me3 in breast cancer transformation

Aim: Epigenetic marks are critical regulators of chromatin and gene activity. Their roles in normal physiology and disease states, including cancer development, still remain elusive. Herein, the epigenomic change of H3K9me3, as well as its potential impacts on gene activity and genome stability, was investigated in an in vitro breast cancer transformation model. Methods: The global H3K9me3 level was studied with western blotting. The distribution of H3K9me3 on chromatin and gene expression was studied with ChIP-Seq and RNA-Seq, respectively. Results: The global H3K9me3 level decreases during transformation and its distribution on chromatin is reprogrammed. By combining with TCGA data, we identified 67 candidate oncogenes, among which five genes are totally novel. Our analysis further links H3K9me3 with transposon activity, and suggests H3K9me3 reduction increases the cell’s sensitivity to DNA damage reagents. Conclusion: H3K9me3 reduction is possibly related with breast cancer transformation by regulating gene expression and chromatin stability during transformation.

Identification of novel prostate cancer drivers using RegNetDriver: a framework for integration of genetic and epigenetic alterations with tissue-specific regulatory network

We report a novel computational method, RegNetDriver, to identify tumorigenic drivers using the combined effects of coding and non-coding single nucleotide variants, structural variants, and DNA methylation changes in the DNase I hypersensitivity based regulatory network. Integration of multi-omics data from 521 prostate tumor samples indicated a stronger regulatory impact of structural variants, as they affect more transcription factor hubs in the tissue-specific network. Moreover, crosstalk between transcription factor hub expression modulated by structural variants and methylation levels likely leads to the differential expression of target genes. We report known prostate tumor regulatory drivers and nominate novel transcription factors (ERF, CREB3L1, and POU2F2), which are supported by functional validation.

Vimentin Levels and Serine 71 Phosphorylation in the Control of Cell-Matrix Adhesions, Migration Speed, and Shape of Transformed Human Fibroblasts

Metastasizing tumor cells show increased expression of the intermediate filament (IF) protein vimentin, which has been used to diagnose invasive tumors for decades. Recent observations indicate that vimentin is not only a passive marker for carcinoma, but may also induce tumor cell invasion. To clarify how vimentin IFs control cell adhesions and migration, we analyzed the nanoscale (30-50 nm) spatial organization of vimentin IFs and cell-matrix adhesions in metastatic fibroblast cells, using three-color stimulated emission depletion (STED) microscopy. We also studied whether wild-type and phospho-deficient or -mimicking mutants of vimentin changed the size and lifetime of focal adhesions (FAs), cell shape, and cell migration, using live-cell total internal reflection imaging and confocal microscopy. We observed that vimentin exists in fragments of different lengths. Short fragments were mostly the size of a unit-length filament and were mainly localized close to small cell-matrix adhesions. Long vimentin filaments were found in the proximity of large FAs. Vimentin expression in these cells caused a reduction in FAs size and an elongated cell shape, but did not affect FA lifetime, or the speed or directionality of cell migration. Expression of a phospho-mimicking mutant (S71D) of vimentin increased the speed of cell migration. Taken together, our results suggest that in highly migratory, transformed mesenchymal cells, vimentin levels control the cell shape and FA size, but not cell migration, which instead is linked to the phosphorylation status of S71 vimentin. These observations are consistent with the possibility that not only levels, but also the assembly status of vimentin control cell migration.

Global histone modification profiling reveals the epigenomic dynamics during malignant transformation in a four-stage breast cancer model

BACKGROUND

Epigenetic regulation has emerged to be the critical steps for tumorigenesis and metastasis. Multiple histone methyltransferase and demethylase have been implicated as tumor suppressors or oncogenes recently. But the key epigenomic events in cancer cell transformation still remain poorly understood.

METHODS

A breast cancer transformation model was established via stably expressing three oncogenes in primary breast epithelial cells. Chromatin immunoprecipitation followed by the next-generation sequencing of histone methylations was performed to determine epigenetic events during transformation. Western blot, quantitative RT-PCR, and immunostaining were used to determine gene expression in cells and tissues.

RESULTS

Histones H3K9me2 and me3, two repressive marks of transcription, decrease in in vitro breast cancer cell model and in vivo clinical tissues. A survey of enzymes related with H3K9 methylation indicated that KDM3A/JMJD1A, a demethylase for H3K9me1 and me2, gradually increases during cancer transformation and is elevated in patient tissues. KDM3A/JMJD1A deficiency impairs the growth of tumors in nude mice and transformed cell lines. Genome-wide ChIP-seq analysis reveals that the boundaries of decreased H3K9me2 large organized chromatin K9 modifications (LOCKs) are enriched with cancer-related genes, such as MYC and PAX3. Further studies show that KDM3A/JMJD1A directly binds to these oncogenes and regulates their transcription by removing H3K9me2 mark.

CONCLUSIONS

Our study demonstrates reduction of histones H3K9 me2 and me3, and elevation of KDM3A/JMJD1A as important events for breast cancer, and illustrates the dynamic epigenomic mechanisms during breast cancer transformation.

NLRX1 Acts as an Epithelial-Intrinsic Tumor Suppressor through the Modulation of TNF-Mediated Proliferation

The mitochondrial Nod-like receptor protein NLRX1 protects against colorectal tumorigenesis through mechanisms that remain unclear. Using mice with an intestinal epithelial cells (IEC)-specific deletion of Nlrx1, we find that NLRX1 provides an IEC-intrinsic protection against colitis-associated carcinogenesis in the colon. These Nlrx1 mutant mice have increased expression of Tnf, Egf, and Tgfb1, three factors essential for wound healing, as well as increased epithelial proliferation during the epithelial regeneration phase following injury triggered by dextran sodium sulfate. In primary intestinal organoids lacking Nlrx1, stimulation with TNF resulted in exacerbated proliferation and expression of the intestinal stem cell markers Olfm4 and Myb. This hyper-proliferation response was associated with increased activation of Akt and NF-κB pathways in response to TNF stimulation. Together, these results identify NLRX1 as a suppressor of colonic tumorigenesis that acts by controlling epithelial proliferation in the intestine during the regeneration phase following mucosal injury.

Karyotype alteration generates the neoplastic phenotypes of SV40-infected human and rodent cells

BACKGROUND

Despite over 50 years of research, it remains unclear how the DNA tumor viruses SV40 and Polyoma cause cancers. Prevailing theories hold that virus-coded Tumor (T)-antigens cause cancer by inactivating cellular tumor suppressor genes. But these theories don't explain four characteristics of viral carcinogenesis: (1) less than one in 10,000 infected cells become cancer cells, (2) cancers have complex individual phenotypes and transcriptomes, (3) recurrent tumors without viral DNA and proteins, (4) preneoplastic aneuploidies and immortal neoplastic clones with individual karyotypes.

RESULTS

As an alternative theory we propose that viral carcinogenesis is a form of speciation, initiated by virus-induced aneuploidy. Since aneuploidy destabilizes the karyotype by unbalancing thousands of genes it catalyzes chain reactions of karyotypic and transcriptomic evolutions. Eventually rare karyotypes evolve that encode cancer-specific autonomy of growth. The low probability of forming new autonomous cancer-species by random karyotypic and transcriptomic variations predicts individual and clonal cancers. Although cancer karyotypes are congenitally aneuploid and thus variable, they are stabilized or immortalized by selections for variants with cancer-specific autonomy. Owing to these inherent variations cancer karyotypes are heterogeneous within clonal margins. To test this theory we analyzed karyotypes and phenotypes of SV40-infected human, rat and mouse cells developing into neoplastic clones. In all three systems we found (1) preneoplastic aneuploidies, (2) neoplastic clones with individual clonal but flexible karyotypes and phenotypes, which arose from less than one in 10,000 infected cells, survived over 200 generations, but were either T-antigen positive or negative, (3) spontaneous and drug-induced variations of neoplastic phenotypes correlating 1-to-1 with karyotypic variations.

CONCLUSIONS

Since all 14 virus-induced neoplastic clones tested contained individual clonal karyotypes and phenotypes, we conclude that these karyotypes have generated and since maintained these neoplastic clones. Thus SV40 causes cancer indirectly, like carcinogens, by inducing aneuploidy from which new cancer-specific karyotypes evolve automatically at low rates. This theory explains the (1) low probability of carcinogenesis per virus-infected cell, (2) the individuality and clonal flexibility of cancer karyotypes, (3) recurrence of neoplasias without viral T-antigens, and (4) the individual clonal karyotypes, transcriptomes and immortality of virus-induced neoplasias - all unexplained by current viral theories.

Integrated genome-wide analysis of genomic changes and gene regulation in human adrenocortical tissue samples

To gain insight into the pathogenesis of adrenocortical carcinoma (ACC) and whether there is progression from normal-to-adenoma-to-carcinoma, we performed genome-wide gene expression, gene methylation, microRNA expression and comparative genomic hybridization (CGH) analysis in human adrenocortical tissue (normal, adrenocortical adenomas and ACC) samples. A pairwise comparison of normal, adrenocortical adenomas and ACC gene expression profiles with more than four-fold expression differences and an adjusted P-value < 0.05 revealed no major differences in normal versus adrenocortical adenoma whereas there are 808 and 1085, respectively, dysregulated genes between ACC versus adrenocortical adenoma and ACC versus normal. The majority of the dysregulated genes in ACC were downregulated. By integrating the CGH, gene methylation and expression profiles of potential miRNAs with the gene expression of dysregulated genes, we found that there are higher alterations in ACC versus normal compared to ACC versus adrenocortical adenoma. Importantly, we identified several novel molecular pathways that are associated with dysregulated genes and further experimentally validated that oncostatin m signaling induces caspase 3 dependent apoptosis and suppresses cell proliferation. Finally, we propose that there is higher number of genomic changes from normal-to-adenoma-to-carcinoma and identified oncostatin m signaling as a plausible druggable pathway for therapeutics.

Analysis of the Human Prostate-Specific Proteome Defined by Transcriptomics and Antibody-Based Profiling Identifies TMEM79 and ACOXL as Two Putative, Diagnostic Markers in Prostate Cancer

To better understand prostate function and disease, it is important to define and explore the molecular constituents that signify the prostate gland. The aim of this study was to define the prostate specific transcriptome and proteome, in comparison to 26 other human tissues. Deep sequencing of mRNA (RNA-seq) and immunohistochemistry-based protein profiling were combined to identify prostate specific gene expression patterns and to explore tissue biomarkers for potential clinical use in prostate cancer diagnostics. We identified 203 genes with elevated expression in the prostate, 22 of which showed more than five-fold higher expression levels compared to all other tissue types. In addition to previously well-known proteins we identified two poorly characterized proteins, TMEM79 and ACOXL, with potential to differentiate between benign and cancerous prostatic glands in tissue biopsies. In conclusion, we have applied a genome-wide analysis to identify the prostate specific proteome using transcriptomics and antibody-based protein profiling to identify genes with elevated expression in the prostate. Our data provides a starting point for further functional studies to explore the molecular repertoire of normal and diseased prostate including potential prostate cancer markers such as TMEM79 and ACOXL.

The reverse evolution from multicellularity to unicellularity during carcinogenesis

Theoretical reasoning suggests that cancer may result from a knockdown of the genetic constraints that evolved for the maintenance of metazoan multicellularity. By characterizing the whole-life history of a xenograft tumour, here we show that metastasis is driven by positive selection for general loss-of-function mutations on multicellularity-related genes. Expression analyses reveal mainly downregulation of multicellularity-related genes and an evolving expression profile towards that of embryonic stem cells, the cell type resembling unicellular life in its capacity of unlimited clonal proliferation. Also, the emergence of metazoan multicellularity ~600 Myr ago is accompanied by an elevated birth rate of cancer genes, and there are more loss-of-function tumour suppressors than activated oncogenes in a typical tumour. These data collectively suggest that cancer represents a loss-of-function-driven reverse evolution back to the unicellular 'ground state'. This cancer evolution model may account for inter-/intratumoural genetic heterogeneity, could explain distant-organ metastases and hold implications for cancer therapy.

The mitochondrial protein NLRX1 controls the balance between extrinsic and intrinsic apoptosis

NLRX1 is a mitochondrial Nod-like receptor (NLR) protein whose function remains enigmatic. Here, we observed that NLRX1 expression was glucose-regulated and blunted by SV40 transformation. In transformed but not primary murine embryonic fibroblasts, NLRX1 expression mediated resistance to an extrinsic apoptotic signal, whereas conferring susceptibility to intrinsic apoptotic signals, such as glycolysis inhibition, increased cytosolic calcium and endoplasmic reticulum stress. In a murine model of colorectal cancer induced by azoxymethane, NLRX1-/- mice developed fewer tumors than wild type mice. In contrast, in a colitis-associated cancer model combining azoxymethane and dextran sulfate sodium, NLRX1-/- mice developed a more severe pathology likely due to the increased sensitivity to dextran sulfate sodium colitis. Together, these results identify NLRX1 as a critical mitochondrial protein implicated in the regulation of apoptosis in cancer cells. The unique capacity of NLRX1 to regulate the cellular sensitivity toward intrinsic versus extrinsic apoptotic signals suggests a critical role for this protein in numerous physiological processes and pathological conditions.

Development of proteome-wide binding reagents for research and diagnostics

Alongside MS, antibodies and other specific protein-binding molecules have a special place in proteomics as affinity reagents in a toolbox of applications for determining protein location, quantitative distribution and function (affinity proteomics). The realisation that the range of research antibodies available, while apparently vast is nevertheless still very incomplete and frequently of uncertain quality, has stimulated projects with an objective of raising comprehensive, proteome-wide sets of protein binders. With progress in automation and throughput, a remarkable number of recent publications refer to the practical possibility of selecting binders to every protein encoded in the genome. Here we review the requirements of a pipeline of production of protein binders for the human proteome, including target prioritisation, antigen design, 'next generation' methods, databases and the approaches taken by ongoing projects in Europe and the USA. While the task of generating affinity reagents for all human proteins is complex and demanding, the benefits of well-characterised and quality-controlled pan-proteome binder resources for biomedical research, industry and life sciences in general would be enormous and justify the effort. Given the technical, personnel and financial resources needed to fulfil this aim, expansion of current efforts may best be addressed through large-scale international collaboration.

Increasing the Content of High-Content Screening: An Overview

Target-based high-throughput screening (HTS) has recently been critiqued for its relatively poor yield compared to phenotypic screening approaches. One type of phenotypic screening, image-based high-content screening (HCS), has been seen as particularly promising.

In this article, we assess whether HCS is as high content as it can be. We analyze HCS publications and find that although the number of HCS experiments published each year continues to grow steadily, the information content lags behind. We find that a majority of high-content screens published so far (60−80%) made use of only one or two image-based features measured from each sample and disregarded the distribution of those features among each cell population. We discuss several potential explanations, focusing on the hypothesis that data analysis traditions are to blame. This includes practical problems related to managing large and multidimensional HCS data sets as well as the adoption of assay quality statistics from HTS to HCS. Both may have led to the simplification or systematic rejection of assays carrying complex and valuable phenotypic information.

We predict that advanced data analysis methods that enable full multiparametric data to be harvested for entire cell populations will enable HCS to finally reach its potential.

Chromosome 3p loss of heterozygosity is associated with a unique metabolic network in clear cell renal carcinoma

Several common oncogenic pathways have been implicated in the emergence of renowned metabolic features in cancer, which in turn are deemed essential for cancer proliferation and survival. However, the extent to which different cancers coordinate their metabolism to meet these requirements is largely unexplored. Here we show that even in the heterogeneity of metabolic regulation a distinct signature encompassed most cancers. On the other hand, clear cell renal cell carcinoma (ccRCC) strongly deviated in terms of metabolic gene expression changes, showing widespread down-regulation. We observed a metabolic shift that associates differential regulation of enzymes in one-carbon metabolism with high tumor stage and poor clinical outcome. A significant yet limited set of metabolic genes that explained the partial divergence of ccRCC metabolism correlated with loss of von Hippel-Lindau tumor suppressor (VHL) and a potential activation of signal transducer and activator of transcription 1. Further network-dependent analyses revealed unique defects in nucleotide, one-carbon, and glycerophospholipid metabolism at the transcript and protein level, which contrasts findings in other tumors. Notably, this behavior is recapitulated by recurrent loss of heterozygosity in multiple metabolic genes adjacent to VHL. This study therefore shows how loss of heterozygosity, hallmarked by VHL deletion in ccRCC, may uniquely shape tumor metabolism.

Oncogenes induce a vimentin filament collapse mediated by HDAC6 that is linked to cell stiffness

Oncogenes deregulate fundamental cellular functions, which can lead to development of tumors, tumor-cell invasion, and metastasis. As the mechanical properties of cells govern cell motility, we hypothesized that oncogenes promote cell invasion by inducing cytoskeletal changes that increase cellular stiffness. We show that the oncogenes simian virus 40 large T antigen, c-Myc, and cyclin E induce spatial reorganization of the vimentin intermediate filament network in cells. At the cellular level, this reorganization manifests as increased width of vimentin fibers and the collapse of the vimentin network. At nanoscale resolution, the organization of vimentin fibers in these oncogene-expressing cells was more entangled, with increased width of the fibers compared with control cells. Expression of these oncogenes also resulted in up-regulation of the tubulin deacetylase histone deacetylase 6 (HDAC6) and altered spatial distribution of acetylated microtubules. This oncogene expression also induced increases in cellular stiffness and promoted the invasive capacity of the cells. Importantly, HDAC6 was required and sufficient for the structural collapse of the vimentin filament network, and was required for increased cellular stiffness of the oncogene-expressing cells. Taken together, these data are consistent with the possibility that oncogenes can induce cellular stiffness via an HDAC6-induced reorganization of the vimentin intermediate filament network.