A central role for TRPS1 in the control of cell cycle and cancer development

Lineage-specific proteome remodeling of diverse lung cancer cells by targeted epigenetic inhibitors

Epigenetic inhibitors exhibit powerful antiproliferative and anticancer activities. However, cellular responses to small-molecule epigenetic inhibition are heterogenous and dependent on factors such as the genetic background, metabolic state, and on-/off-target engagement of individual small-molecule drugs. To determine the mechanisms that drive these heterogeneous cellular responses, we quantified chromatin, proteome, and transcriptome remodeling due to histone deacetylase inhibitor (HDACi) -treated cells derived from diverse genetic backgrounds. We utilized high-throughput sample multiplexed proteomics and integrated intelligent data acquisition methods to map proteomes of cancer cell lines in response to HDACi. We determined cell type-specific and ubiquitous cellular responses based on the quantification of 10,621 total proteins. We then established how coordinated remodeling of the proteome, transcriptome and chromatin state of HDACi treated cancer cells revealed convergent (JUN, MAP2K3, CDKN1A) and divergent (CCND3, ASF1B, BRD7) molecular phenotypes. HDACi-regulated proteins differ greatly across cell lines owing to heterogeneous molecular states of these cell lines. Finally, we demonstrated that HDACi treatment drove a highly cell-type specific response that may in part be explained by cell line-specific off-target drug engagement.

TRPS1 function beyond breast: A retrospective immunohistochemical study on non-breast cytology specimens

INTRODUCTION

Trichorhinophalangeal syndrome type 1 (TRPS1) has emerged as a reliable immunohistochemistry (IHC) marker for identifying breast origin in metastatic carcinomas. This study investigates the utility of TRPS1 IHC in non-breast cytology specimens.

MATERIALS AND METHODS

A retrospective search of our pathology database for the year 2021 identified fluids (pleural and peritoneal) and liver, lung and bone fine needle aspirations (FNAs) with surgical follow-up confirming non-breast metastatic carcinomas. Cell blocks from cases with sufficient neoplastic cells underwent immunostaining using a rabbit polyclonal antibody against human TRPS1. Cases lacking tumor on deeper levels after the original work-up were excluded from the study. Two pathologists independently interpreted the TRPS1 staining.

RESULTS

Of 136 cases assessed, 31 (22.79%) exhibited positive TRPS1 staining, while 105 (77.21%) were nonreactive. Positivity rates were observed in tumors of Mullerian cell origin, gastrointestinal tract (GIT), and lung origin at 28.85%, 25%, and 21.57%, respectively. Of the tumors of Mullerian cell origin 10 (66.67%) were serous carcinomas, 4 (26.67%) were endometrioid carcinomas, and one (6.67%) was a clear cell carcinoma. Lung tumors comprised seven (63.64%) squamous cell carcinomas and four (36.36%) adenocarcinomas, while the gastrointestinal tumors consisted of 14 (80%) adenocarcinomas and one (20%) squamous cell carcinoma.

CONCLUSIONS

Although recognized as a sensitive marker for mammary carcinomas, TRPS1 immunostaining was also detected in Mullerian, lung, and GIT carcinomas. This highlights the significance of being cautious when depending solely on TRPS1 immunostaining to distinguish metastatic breast tumors.

A Comparative Evaluation of TRPS1 and GATA3 in adenoid cystic, secretory, and acinic cell carcinomas of the breast and salivary gland

GATA3 is the most used marker to determine tumors' breast origin, but its diagnostic value in triple-negative breast cancer (TNBC) is limited. The newly identified TRPS1 is highly sensitive and specific for breast carcinoma, especially TNBC. Here, we compared the utility of TRPS1 and GATA3 expression in a subset of salivary gland-type breast tumors (including adenoid cystic, acinic cell, and secretory carcinomas [AdCC, ACC, and SC, respectively]), and we compared TRPS1 and GATA3 expression of such tumors with head and neck (H&N) and AdCC of upper respiratory tumors. TRPS1 was strongly expressed in basaloid TNBC and AdCCs with solid components, including 100 % of mixed and solid breast AdCCs. However, TRPS1 was positive in only 50 % cribriform AdCCs. Expression patterns of TRPS1 in H&N and upper respiratory AdCC were similar. TRPS1 was positive in 30 % of H&N cribriform AdCCs but was strongly expressed in mixed AdCC (67 %) and solid AdCC (100 %). In the upper respiratory AdCCs, TRPS1 was positive in 58.4 % of cribriform AdCCs and positive in 100 % of AdCCs with solid components. On the contrary, GATA3 was negative in predominant AdCCs of the breast, H&N, and upper respiratory tract. These data show that GATA3 and TRPS1 expression varies AdCCs. In addition, TRPS1 and GATA3 expression patterns were similar SC and ACC of breast and H&N. Both markers were positive in SC and negative in ACC. Therefore, TRPS1 and GATA3 cannot be used to differentiate salivary gland-type carcinomas of breast origin from those of upper respiratory or H&N origin.

Peripheral Ischemia Imprints Epigenetic Changes in Hematopoietic Stem Cells to Propagate Inflammation and Atherosclerosis

BACKGROUND

Peripheral ischemia caused by peripheral artery disease is associated with systemic inflammation, which may aggravate underlying comorbidities such as atherosclerosis and heart failure. However, the mechanisms of increased inflammation and inflammatory cell production in patients with peripheral artery disease remain poorly understood.

METHODS

We used peripheral blood collected from patients with peripheral artery disease and performed hind limb ischemia (HI) in Apoe-/- mice fed a Western diet and C57BL/6J mice with a standard laboratory diet. Bulk and single-cell RNA sequencing analysis, whole-mount microscopy, and flow cytometry were performed to analyze hematopoietic stem and progenitor cell (HSPC) proliferation, differentiation, and relocation.

RESULTS

We observed augmented numbers of leukocytes in the blood of patients with peripheral artery disease and Apoe-/- mice with HI. RNA sequencing and whole-mount imaging of the bone marrow revealed HSPC migration into the vascular niche from the osteoblastic niche and their exaggerated proliferation and differentiation. Single-cell RNA sequencing demonstrated alterations in the genes responsible for inflammation, myeloid cell mobilization, and HSPC differentiation after HI. Heightened inflammation in Apoe-/- mice after HI aggravated atherosclerosis. Surprisingly, bone marrow HSPCs expressed higher amounts of the receptors for IL (interleukin)-1 and IL-3 after HI. Concomitantly, the promoters of Il1r1 and Il3rb had augmented H3K4me3 and H3K27ac marks after HI. Genetic and pharmacological inhibition of these receptors resulted in suppressed HSPC proliferation, reduced leukocyte production, and ameliorated atherosclerosis.

CONCLUSIONS

Our findings demonstrate increased inflammation, HSPC abundance in the vascular niches of the bone marrow, and elevated IL-3Rb and IL-1R1 (IL-1 receptor 1) expression in HSPC following HI. Furthermore, the IL-3Rb and IL-1R1 signaling plays a pivotal role in HSPC proliferation, leukocyte abundance, and atherosclerosis aggravation after HI.

A Genome-Wide Association Study Identified Novel Genetic Susceptibility Loci for Oral Cancer in Taiwan

Taiwan has the highest incidence rate of oral cancer in the world. Although oral cancer is mostly an environmentally induced cancer, genetic factors also play an important role in its etiology. Genome-wide association studies (GWAS) have identified nine susceptibility regions for oral cancers in populations of European descent. In this study, we performed the first GWAS of oral cancer in Taiwan with 1529 cases and 44,572 controls. We confirmed two previously reported loci on the 6p21.33 (HLA-B) and 6p21.32 (HLA-DQ gene cluster) loci, highlighting the importance of the human leukocyte antigen and, hence, the immunologic mechanisms in oral carcinogenesis. The TERT-CLMPT1L locus on 5p15.33, the 4q23 ADH1B locus, and the LAMC3 locus on 9q34.12 were also consistent in the Taiwanese. We found two new independent loci on 6p21.32, rs401775 in SKIV2L gene and rs9267798 in TNXB gene. We also found two suggestive novel Taiwanese-specific loci near the TPRS1 gene on 8q23.3 and in the TMED3 gene on 15q25.1. This study identified both common and unique oral cancer susceptibility loci in the Taiwanese as compared to populations of European descent and shed significant light on the etiology of oral cancer in Taiwan.

Functional mechanisms of TRPS1 in disease progression and its potential role in personalized medicine

The gene of transcriptional repressor GATA binding 1 (TRPS1), as an atypical GATA transcription factor, has received considerable attention in a plethora of physiological and pathological processes, and may become a promising biomarker for targeted therapies in diseases and tumors. However, there still lacks a comprehensive exploration of its functions and promising clinical applications. Herein, relevant researches published in English from 2000 to 2022 were retrieved from PubMed, Google Scholar and MEDLINE, concerning the roles of TRPS1 in organ differentiation and tumorigenesis. This systematic review predominantly focused on summarizing the structural characteristics and biological mechanisms of TRPS1, its involvement in tricho-rhino-phalangeal syndrome (TRPS), its participation in the development of multiple tissues, the recent advances of its vital features in metabolic disorders as well as malignant tumors, in order to prospect its potential applications in disease detection and cancer targeted therapy. From the clinical perspective, the deeply and thoroughly understanding of the complicated context-dependent and cell-lineage-specific mechanisms of TRPS1 would not only gain novel insights into the complex etiology of diseases, but also provide the fundamental basis for the development of therapeutic drugs targeting both TRPS1 and its critical cofactors, which would facilitate individualized treatment.

Some Pleomorphic Adenomas of the Breast Share PLAG1 Rearrangements with the Analogous Tumor of the Salivary Glands

AIMS

Pleomorphic adenoma (PA) of the breast, and especially its malignant transformation, is extremely rare and represents a diagnostic pitfall. Molecular alterations in this entity have not been investigated. We aimed to examine the clinicopathologic features of our breast PAs and perform molecular analysis.

METHODS

Seven cases of breast PA including two cases of carcinoma ex PA were analyzed. PLAG1 and HMGA2 gene rearrangements were assayed by FISH and RNA-Seq, respectively. RT-PCR and Sanger sequencing were used to verify RNA sequencing results.

RESULTS

All seven cases of breast PA occurred in women. The histological features were similar to the analogous tumor in salivary glands, including a dual epithelial-myoepithelial component and negativity of ER, PR, and HER2 by immunohistochemistry. Of the two cases with carcinoma ex PA, one demonstrated minimal invasion and one was extensively invasive. PLAG1 rearrangements were identified in two cases (28.6%), but no rearrangements of HMG2A were found. A novel fusion product in PAs, TRPS1-PLAG1, was identified in one case. No patients had recurrence or metastasis with a follow-up period of 6 to 158 months.

CONCLUSIONS

Breast PA is rare, but it is an important differential diagnosis of breast pathology with the potential to develop carcinoma ex PA. We reported a novel TRPS1-PLAG1 fusion gene in breast PA.

Clinical significance and biological function of transcriptional repressor GATA binding 1 in gastric cancer: a study based on data mining, RT-qPCR, immunochemistry, and vitro experiment

Transcriptional repressor GATA binding 1 (TRPS1) is a newly discovered transcription factor, which has been reported in many tumors, except for gastric cancer (GC). In this study, we aimed to grope for clinical significance and biological function of TRPS1 in GC. TRPS1 expression in GC and its relationship with clinicopathological features were analyzed based on public databases, and verified by immunohistochemistry and RT-qPCR. Kaplan-Meier survival curve and Cox regression model were used to estimate the influence of TRPS1 on the univariate prognosis and multivariate survival risk factors of GC. The effects of TRPS1 on malignant biological behaviors of GC cells were studied by CCK8 cell proliferation, scratch test, and Transwell assay. The function of TRPS1 was further analyzed by signaling pathway analysis. TRPS1 mRNA expression in GC tissues was up-regulated and was of great significance in some prognostic factors. Protein expression of TRPS1 in tumor tissues was significantly higher than that in paracancerous tissues. Over-expression of TRPS1 was a poor prognostic indicator for GC patients. TRPS1 knockdown could inhibit the proliferation, migration, and invasion of GC cells. The important role of TRPS1 was in the extracellular matrix, and it was involved in actin binding and proteoglycan in cancer. The hub genes of TRPS1 (FN1, ITGB1) were defined. TRPS1 may be a tumor promoter and promote the development of GC by influencing the malignant biological behaviors of GC. TRPS1 is expected to be a key diagnostic and prognostic indicator for GC patients.

LncRNA PVT1 Regulates TRPS1 Expression in Breast Cancer by Sponging miR-543

BACKGROUND

Breast cancer is the most common female malignancy with high invasion and metastasis abilities. Studies have shown that long non-coding RNA (lncRNA) plasmacytoma variant translocation 1 gene (PVT1) is an oncogene and is positively correlated with progression and metastasis of breast tumors. However, the detailed mechanism of PVT1 in breast cancer tumorigenesis is not fully understood.

METHODS

Real-time polymerase quantitative chain reaction (RT-qPCR) was performed to identify the expression levels of PVT1, miR-543 and trichorhinophalangeal syndrome-1 gene (TRPS1) in breast cancer tissues and cells. Cell proliferation was measured by plate clone formation and 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl-2H-tetrazol-3-ium bromide (MTT) assay. Apoptosis and motility of MCF-7 and MDA-MB-436 cells were assessed with flow cytometry assay and transwell migration and invasion analyses, respectively. In addition, a model was established to probe the function of PVT1 silencing in vivo. The target relationship among PVT1, miR-543 or TRPS1 was confirmed by dual-luciferase reporter analysis, RNA immunoprecipitation (RIP) and RNA pull down assays. The protein expression level of TRPS1 was evaluated with Western blot assay.

RESULTS

PVT1 expression was upregulated in breast cancer tissues and cell lines. In addition, PVT1 silencing inhibited breast cancer cell growth and motility, while increased apoptosis. Meanwhile, the effects of PVT1 or miR-543 could be reversed by introducing overexpressed plasmid of miR-543 or TRPS1 in breast cancer cell lines, respectively.

CONCLUSION

Knockdown of PVT1 repressed breast cancer cell growth and motility, and induced apoptosis in vitro and reduced tumor volume and weight in vivo. Mechanically, the overexpression of PVT1 enhanced TRPS1 level by negatively targeted miR-543 in breast cancer.

TRPS1 Confers Multidrug Resistance of Breast Cancer Cells by Regulating BCRP Expression

Multidrug resistance (MDR) is the major obstruction in the successful treatment of breast cancer (BCa). The elucidation of molecular events that confer chemoresistance of BCa is of major therapeutic importance. Several studies have elucidated the correlation of TRPS1 and BCa. Here we focused on the role of TRPS1 in acquisition of chemoresistance, and reported a unique role of TRPS1 renders BCa cells resistant to chemotherapeutic drugs via the regulation of BCRP expression. Bioinformation analysis based on publicly available BCa data suggested that TRPS1 overexpression linked to chemoresistance. Mechanistically, TRPS1 regulated BCRP expression and efflux transportation. Overexpression of TRPS1 led to upregulation of BCRP while its inhibition resulted in repression of BCRP. The correlation of TRPS1 and BCRP was further confirmed by immunohistochemistry in 180 BCa samples. MTT assay demonstrated that manipulation of TRPS1 expression affects the chemosensitivity of BCa cells. In total, high expression of TRPS1 confers MDR of BCa which is mediated by BCRP. Our data demonstrated a new insight into mechanisms and strategies to overcome chemoresistance in BCa.

ENGINE: a Phase III randomized placebo controlled study of enzastaurin/R-CHOP as frontline therapy in high-risk diffuse large B-cell lymphoma patients with the genomic biomarker DGM1

While combination of rituximab, cyclophosphamide, doxorubicin, vincristine and prednisone (R-CHOP) cures most patients with diffuse large B-cell lymphoma (DLBCL), those with high-risk international prognostic index disease have inferior survival. Enzastaurin as a potent inhibitor of PKC-β and PI3K/AKT pathway suppressor has been tested in many clinical trials including two key studies in DLBCL: Phase III maintenance study (Preventing Relapse in Lymphoma Using Daily Enzastaurin [PRELUDE]) and a first-line Phase II study (S028). DNA extracted from PRELUDE patients' blood samples was retrospectively genotyped identifying a novel genetic biomarker, DGM1 that showed high correlation with response to enzastaurin. A similar finding observed in the S028 study suggested that addition of enzastaurin to R-CHOP may significantly improve outcomes as frontline therapy for high-risk DGM1 positive DLBCL patients. ENGINE is a global, multicenter, placebo-controlled and randomized study to compare the effect of R-CHOP/enzastaurin as frontline treatment in high-risk DLBCL patients. The primary end point for this study is overall survival in patients who are DGM1 positive. Clinical Trial Registration Identifier: NCT03263026.

An Effective Graph Clustering Method to Identify Cancer Driver Modules

Identifying the molecular modules that drive cancer progression can greatly deepen the understanding of cancer mechanisms and provide useful information for targeted therapies. Most methods currently addressing this issue primarily use mutual exclusivity without making full use of the extra layer of module property. In this paper, we propose MCLCluster to identity cancer driver modules, which use somatic mutation data, Cancer Cell Fraction (CCF) data, gene functional interaction network and protein-protein interaction (PPI) network to derive the module property on mutual exclusivity, connectivity in PPI network and functionally similarity of genes. We have taken three effective measures to ensure the effectiveness of our algorithm. First, we use CCF data to choose stronger signals and more confident mutations. Second, the weighted gene functional interaction network is used to quantify the gene functional similarity in PPI. The third, graph clustering method based on Markov is exploited to extract the candidate module. MCLCluster is tested in the two TCGA datasets (GBM and BRCA), and identifies several well-known oncogenes driver modules and some modules with functionally associated driver genes. Besides, we compare it with Multi-Dendrix, FSME Cluster and RME in simulated dataset with background noise and passenger rate, MCLCluster outperforming all of these methods.

Novel PLAG1 Gene Rearrangement Distinguishes a Subset of Uterine Myxoid Leiomyosarcoma From Other Uterine Myxoid Mesenchymal Tumors

Genetic alterations in uterine myxoid leiomyosarcoma are unknown. We investigate the clinicopathologic features of 19 uterine tumors previously diagnosed as myxoid leiomyosarcomas in which tumoral RNA was subjected to targeted RNA sequencing. PLAG1, BCOR, BCORL1, HMGA2, and ALK break-apart fluorescence in situ hybridization (FISH) and BCOR, PLAG1, and ALK immunohistochemistry were performed in cases which failed or lacked fusions by sequencing. The diagnosis of myxoid leiomyosarcoma was confirmed in 15 cases after exclusion of 4 tumors with BCOR and ALK rearrangements. These 15 patients presented at a median age of 50 years with stage I (3), II (2), III (2), and IV (1) tumors, respectively; stage was unknown in 7 cases. Tumor size ranged from 10 to 24 cm. Matrix was myxoid in all tumors and also eosinophilic in 2. Cells were spindled, epithelioid, and both in 10, 2, and 3 tumors and showed mild, moderate, and severe nuclear atypia in 3, 8, and 4 tumors, respectively. Mitotic index ranged from <1 to 14/10 HPF, while tumor necrosis was present in 6 (40%). Novel TRPS1-PLAG1 or RAD51B-PLAG1 fusions were detected by sequencing in 4 tumors, 3 of which were also confirmed by FISH. Diffuse PLAG1 expression was seen in 7 tumors, including 4 with PLAG1 rearrangement. No morphologic differences were seen among PLAG1 fusion-positive and fusion-negative tumors. No PLAG1, HMGA2, ALK, BCOR, or BCORL1 rearrangements were detected by FISH in 11 tumors. On the basis of sequencing and FISH results, PLAG1 rearrangements resulting in PLAG1 expression underpin ~25% of myxoid leiomyosarcomas and may serve as a useful diagnostic biomarker. Immunohistochemistry, targeted RNA sequencing, and/or FISH may distinguish myxoid leiomyosarcoma from its morphologic mimics.

Reciprocal Regulation of TRPS1 and miR-221 in Intervertebral Disc Cells

Intervertebral disc (IVD), a moderately moving joint located between the vertebrae, has a limited capacity for self-repair, and treating injured intervertebral discs remains a major challenge. The development of innovative therapies to reverse IVD degeneration relies primarily on the discovery of key molecules that, occupying critical points of regulatory mechanisms, can be proposed as potential intradiscal injectable biological agents. This study aimed to elucidate the underlying mechanism of the reciprocal regulation of two genes differently involved in IVD homeostasis, the miR-221 microRNA and the TRPS1 transcription factor. Human lumbar IVD tissue samples and IVD primary cells were used to specifically evaluate gene expression and perform functional analysis including the luciferase gene reporter assay, chromatin immunoprecipitation, cell transfection with hTRPS1 overexpression vector and antagomiR-221. A high-level expression of TRPS1 was significantly associated with a lower pathological stage, and TRPS1 overexpression strongly decreased miR-221 expression, while increasing the chondrogenic phenotype and markers of antioxidant defense and stemness. Additionally, TRPS1 was able to repress miR-221 expression by associating with its promoter and miR-221 negatively control TRPS1 expression by targeting the TRPS1-3'UTR gene. As a whole, these results suggest that, in IVD cells, a double-negative feedback loop between a potent chondrogenic differentiation suppressor (miR-221) and a regulator of axial skeleton development (TRPS1) exists. Our hypothesis is that the hostile degenerated IVD microenvironment may be counteracted by regenerative/reparative strategies aimed at maintaining or stimulating high levels of TRPS1 expression through inhibition of one of its negative regulators such as miR-221.

Lost or Forgotten: The nuclear cathepsin protein isoforms in cancer

While research into the role of cathepsins has been progressing at an exponential pace over the years, research into their respective isoform proteins has been less frenetic. In view of the functional and biological potential of such protein isoforms in model systems for cancer during their initial discovery, much later they have offered a new direction in the field of cathepsin basic and applied research. Consequently, the analysis of such isoforms has laid strong foundations in revealing other important regulatory aspects of the cathepsin proteins in general. In this review article, we address these key aspects of cathepsin isoform proteins, with particular emphasis on how they have shaped what is now known in the context of nuclear cathepsin localization and what potential these hold as nuclear-based therapeutic targets in cancer.

A Novel Method for Identifying the Potential Cancer Driver Genes Based on Molecular Data Integration

The identification of the cancer driver genes is essential for personalized therapy. The mutation frequency of most driver genes is in the middle (2-20%) or even lower range, which makes it difficult to find the driver genes with low-frequency mutations. Other forms of genomic aberrations, such as copy number variations (CNVs) and epigenetic changes, may also reflect cancer progression. In this work, a method for identifying the potential cancer driver genes (iPDG) based on molecular data integration is proposed. DNA copy number variation, somatic mutation, and gene expression data of matched cancer samples are integrated. In combination with the method of iKEEG, the "key genes" of cancer are identified, and the change in their expression levels is used for auxiliary evaluation of whether the mutated genes are potential drivers. For a mutated gene, the concept of mutational effect is defined, which takes into account the effects of copy number variation, mutation gene itself, and its neighbor genes. The method mainly includes two steps: the first step is data preprocessing. First, DNA copy number variation and somatic mutation data are integrated. Then, the integrated data are mapped to a given interaction network, and the diffusion kernel is used to form the mutation effect matrix. The second step is to obtain the key genes by using the iKGGE method, and construct the connection matrix by means of the gene expression data of the key genes and mutation impact matrix of the mutated genes. Experiments on TCGA breast cancer and Glioblastoma multiforme datasets demonstrate that iPDG is effective not only to identify the known cancer driver genes but also to discover the rare potential driver genes. When measured by functional enrichment analysis, we find that these genes are clearly associated with these two types of cancers.

Folate Repletion after Deficiency Induces Irreversible Genomic and Transcriptional Changes in Human Papillomavirus Type 16 (HPV16)-Immortalized Human Keratinocytes

Supplementation of micronutrients like folate is a double-edged sword in terms of their ambivalent role in cell metabolism. Although several epidemiological studies support a protective role of folate in carcinogenesis, there are also data arguing for an opposite effect. To address this issue in the context of human papillomavirus (HPV)-induced transformation, the molecular events of different folate availability on human keratinocytes immortalized by HPV16 E6 and E7 oncoproteins were examined. Several sublines were established: Control (4.5 µM folate), folate deficient (0.002 µM folate), and repleted cells (4.5 µM folate). Cells were analyzed in terms of oncogene expression, DNA damage and repair, karyotype changes, whole-genome sequencing, and transcriptomics. Here we show that folate depletion irreversibly induces DNA damage, impairment of DNA repair fidelity, and unique chromosomal alterations. Repleted cells additionally underwent growth advantage and enhanced clonogenicity, while the above mentioned impaired molecular properties became even more pronounced. Overall, it appears that a period of folate deficiency followed by repletion can shape immortalized cells toward an anomalous phenotype, thereby potentially contributing to carcinogenesis. These observations should elicit questions and inquiries for broader additional studies regarding folate fortification programs, especially in developing countries with micronutrient deficiencies and high HPV prevalence.

Atypical GATA protein TRPS1 plays indispensable roles in mouse two-cell embryo

Zygotic genome activation (ZGA) is one of the most critical events at the beginning of mammalian preimplantation embryo development (PED). The mechanisms underlying mouse ZGA remain unclear although it has been widely studied. In the present study, we identified that tricho-rhino-phalangeal syndrome 1 (TRPS1), an atypical GATA family member, is an important factor for ZGA in mouse PED. We found that the Trps1 mRNA level peaked at the one-cell stage while TRPS1 protein did so at the two/four-cell stage. Knockdown of Trps1 by the microinjection of Trps1 siRNA reduced the developmental rate of mouse preimplantation embryos by approximately 30%, and increased the expression of ZGA marker genes MuERV-L and Zscan4d via suppressing the expression of major histone markers H3K4me3 and H3K27me3. Furthermore, Trps1 knockdown decreased the expression of Sox2 but increased Oct4 expression. We conclude that TRPS1 may be indispensable for zygotic genome activation during mouse PED.

Atypical GATA transcription factor TRPS1 represses gene expression by recruiting CHD4/NuRD(MTA2) and suppresses cell migration and invasion by repressing TP63 expression

Transcriptional repressor GATA binding 1 (TRPS1), an atypical GATA transcription factor, functions as a transcriptional repressor and is also implicated in human cancers. However, the underlying mechanism of TRPS1 contributing to malignancy remains obscure. In the current study, we report that TRPS1 recognizes both gene proximal and distal transcription start site (TSS) sequences to repress gene expression. Co-IP mass spectrometry and biochemical studies showed that TRPS1 binds to CHD4/NuRD(MTA2). Genome-wide and molecular studies revealed that CHD4/NuRD(MTA2) is required for TRPS1 transcriptional repression. Mechanically, TRPS1 and CHD4/NuRD(MTA2) form precision-guided transcriptional repression machinery in which TRPS1 guides the machinery to specific target sites by recognizing GATA elements, and CHD4/NuRD(MTA2) represses the transcription of target genes. Furthermore, TP63 was identified and validated to be a direct target of TRPS1-CHD4/NuRD(MTA2) complex, which represses TP63 expression by involving decommission of TP63 enhancer in the described precision-guided manner, leading to a reduction of the ΔNp63 level and contributing to migration and invasion of cancer cells.

Transcriptional repressor GATA binding 1-mediated repression of SRY-box 2 expression suppresses cancer stem cell functions and tumor initiation

Cancer stem cells (CSCs) have been reported in a variety of cancers. SRY-box 2 (SOX2) is a member of the SOX family of transcription factors and has been shown to play a critical role in maintaining the functions of CSCs and promoting tumor initiation. However, the underlying mechanisms for the transcriptional regulation of the SOX2 gene in CSCs are unclear. In this study, using in silico and experimental approaches, we identified transcriptional repressor GATA binding 1 (TRPS1), an atypical GATA-type transcription factor, as a critical transcriptional regulator that represses SOX2 expression and thereby suppresses cancer stemness and tumorigenesis. Mechanistically, TRPS1 repressed SOX2 expression by directly targeting the consensus GATA-binding element in the SOX2 promoter as elucidated by ChIP and luciferase reporter assays. Of note, in vitro mammosphere formation assays in culture and in vivo xenograft tumor initiation experiments in mouse models revealed that TRPS1-mediated repression of SOX2 expression suppresses CSC functions and tumor initiation. Taken together, our study provides detailed mechanistic insights into CSC functions and tumor initiation by the TRPS1-SOX2 axis.

TRPS1 regulates oestrogen receptor binding and histone acetylation at enhancers

The chromatin state is finely tuned to regulate function and specificity for transcription factors such as oestrogen receptor alpha (ER), which contributes to cell growth in breast cancer. ER transcriptional potential is mediated, in large part, by the specific associated proteins and co-factors that interact with it. Despite the identification and characterisation of several ER coregulators, a complete and systematic view of ER-regulating chromatin modifiers is lacking. By exploiting a focused siRNA screen that investigated the requirement for a library of 330 chromatin regulators in ER-mediated cell growth, we find that the NuRD and coREST histone deacetylation complexes are critical for breast cancer cell proliferation. Further, by proteomic and genomics approaches, we discover the transcription factor TRPS1 to be a key interactor of the NuRD and coREST complexes. Interestingly, TRPS1 gene amplification occurs in 28% of human breast tumours and is associated with poor prognosis. We propose that TRPS1 is required to repress spurious binding of ER, where it contributes to the removal of histone acetylation. Our data suggest that TRPS1 is an important ER-associated transcriptional repressor that regulates cell proliferation, chromatin acetylation and ER binding at the chromatin of cis-regulatory elements.

Tricho-rhino-phalangeal syndrome 1 protein functions as a scaffold required for ubiquitin-specific protease 4-directed histone deacetylase 2 de-ubiquitination and tumor growth

Background

Although numerous studies have reported that tricho-rhino-phalangeal syndrome type I (TRPS1) protein, the only reported atypical GATA transcription factor, is overexpressed in various carcinomas, the underlying mechanism(s) by which it contributes to cancer remain unknown.

Methods

Both overexpression and knockdown of TRPS1 assays were performed to examine the effect of TRPS1 on histone deacetylase 2 (HDAC2) protein level and luminal breast cancer cell proliferation. Also, RT-qRCR, luciferase reporter assay and RNA-sequencing were used for transcription detection. Chromatin immunoprecipitation (ChIP) using H4K16ac antibody in conjunction with qPCR was used for determining H4K16ac levels in targeted genes. Furthermore, in vitro cell proliferation assay and in vivo tumor xenografts were used to detect the effect of TRPS1 on tumor growth.

Results

We found that TRPS1 scaffolding recruits and enhances interaction between USP4 and HDAC2 leading to HDAC2 de-ubiquitination and H4K16 deacetylation. We detected repression of a set of cellular growth-related genes by the TRPS1-USP4-HDAC2 axis indicating it is essential in tumor growth. In vitro and in vivo experiments confirmed that silencing TRPS1 reduced tumor growth, whereas overexpression of HDAC2 restored tumor growth.

Conclusion

Our study deciphered the TRPS1-USP4-HDAC2 axis as a novel mechanism that contributes to tumor growth. Significantly, our results revealed the scaffolding function of TPRS1 in USP4-directed HDAC2 de-ubiquitination and provided new mechanistic insights into the crosstalk between TRPS1, ubiquitin, and histone modification systems leading to tumor growth.

Electronic supplementary material

The online version of this article (10.1186/s13058-018-1018-7) contains supplementary material, which is available to authorized users.

Network analysis of SRC-1 reveals a novel transcription factor hub which regulates endocrine resistant breast cancer

Steroid receptor coactivator 1 (SRC-1) interacts with nuclear receptors and other transcription factors (TFs) to initiate transcriptional networks and regulate downstream genes which enable the cancer cell to evade therapy and metastasise. Here we took a top-down discovery approach to map out the SRC-1 transcriptional network in endocrine resistant breast cancer. First, rapid immunoprecipitation mass spectrometry of endogenous proteins (RIME) was employed to uncover new SRC-1 TF partners. Next, RNA sequencing (RNAseq) was undertaken to investigate SRC-1 TF target genes. Molecular and patient-derived xenograft studies confirmed STAT1 as a new SRC-1 TF partner, important in the regulation of a cadre of four SRC-1 transcription targets, NFIA, SMAD2, E2F7 and ASCL1. Extended network analysis identified a downstream 79 gene network, the clinical relevance of which was investigated in RNAseq studies from matched primary and local-recurrence tumours from endocrine resistant patients. We propose that SRC-1 can partner with STAT1 independently of the estrogen receptor to initiate a transcriptional cascade and control regulation of key endocrine resistant genes.

Trps1 is associated with the multidrug resistance of lung cancer cell by regulating MGMT gene expression

Multidrug resistance (MDR) often leads to chemotherapy failure of lung cancer and has been linking to the cellular expression of several DNA transcription- and repair-related genes such as Trps1 and MGMT. However, their roles in the formation of MDR are largely unknown. In this study, overexpression/knockdown, luciferase assay and ChIP assay were performed to study the relationship between Trps1 and MGMT, as well as their roles in MDR formation. Our results demonstrated that Trps1 and MGMT expression both increased in drug-resistant lung cancer cell line (H446/CDDP). Silencing of Trps1 resulted in downregulation of MGMT expression and decrease in the multidrug sensitivity of H446/CDDP cells, while Trps1 overexpression exhibited the opposite effects in H446 cells. Ectopic expression of MGMT had no effect on Trps1 expression, but enhanced the IC50 values of H446 cells or rescued the IC50 values of Trps1-silenced H446/CDDP cells in treatment of multidrug. Our data further showed that, mechanistically, Trps1 acted as a transcription activator that directly induced MGMT transcription by binding to the MGMT promoter. Taken together, we consider that upregulation of Trps1 induces MGMT transcription contributing to the formation of MDR in lung cancer cells. Our findings proved potential targets for reversing MDR in clinical chemotherapy of lung cancer.

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.

Quionolone carboxylic acid derivatives as HIV-1 integrase inhibitors: Docking-based HQSAR and topomer CoMFA analyses

Quionolone carboxylic acid derivatives as inhibitors of HIV-1 integrase were investigated as a potential class of drugs for the treatment of acquired immunodeficiency syndrome (AIDS). Hologram quantitative structure-activity relationships (HQSAR) and translocation comparative molecular field vector analysis (topomer CoMFA) were applied to a series of 48 quionolone carboxylic acid derivatives. The most effective HQSAR model was obtained using atoms and bonds as fragment distinctions: cross-validation q² = 0.796, standard error of prediction SD_CV = 0.36, the non-cross-validated r² = 0.967, non-cross validated standard error SD = 0.17, the correlation coefficient of external validation Q_ext² = 0.955, and the best hologram length HL = 180. topomer CoMFA models were built based on different fragment cutting models, with the most effective model of q² = 0.775, SD_CV = 0.37, r² = 0.967, SD = 0.15, Q_ext² = 0.915, and F = 163.255. These results show that the models generated form HQSAR and topomer CoMFA were able to effectively predict the inhibitory potency of this class of compounds. The molecular docking method was also used to study the interactions of these drugs by docking the ligands into the HIV-1 integrase active site, which revealed the likely bioactive conformations. This study showed that there are extensive interactions between the quionolone carboxylic acid derivatives and THR80, VAL82, GLY27, ASP29, and ARG8 residues in the active site of HIV-1 integrase. These results provide useful insights for the design of potent new inhibitors of HIV-1 integrase.

Expression of miRNA-26b-5p and its target TRPS1 is associated with radiation exposure in post-Chernobyl breast cancer

Ionizing radiation is a well-recognized risk factor for the development of breast cancer. However, it is unknown whether radiation-specific molecular oncogenic mechanisms exist. We investigated post-Chernobyl breast cancers from radiation-exposed female clean-up workers and nonexposed controls for molecular changes. Radiation-associated alterations identified in the discovery cohort (n = 38) were subsequently validated in a second cohort (n = 39). Increased expression of hsa-miR-26b-5p was associated with radiation exposure in both of the cohorts. Moreover, downregulation of the TRPS1 protein, which is a transcriptional target of hsa-miR-26b-5p, was associated with radiation exposure. As TRPS1 overexpression is common in sporadic breast cancer, its observed downregulation in radiation-associated breast cancer warrants clarification of the specific functional role of TRPS1 in the radiation context. For this purpose, the impact of TRPS1 on the transcriptome was characterized in two radiation-transformed breast cell culture models after siRNA-knockdown. Deregulated genes upon TRPS1 knockdown were associated with DNA-repair, cell cycle, mitosis, cell migration, angiogenesis and EMT pathways. Furthermore, we identified the interaction partners of TRPS1 from the transcriptomic correlation networks derived from gene expression data on radiation-transformed breast cell culture models and sporadic breast cancer tissues provided by the TCGA database. The genes correlating with TRPS1 in the radiation-transformed breast cell lines were primarily linked to DNA damage response and chromosome segregation, while the transcriptional interaction partners in the sporadic breast cancers were mostly associated with apoptosis. Thus, upregulation of hsa-miR-26b-5p and downregulation of TRPS1 in radiation-associated breast cancer tissue samples suggests these molecules representing radiation markers in breast cancer.

Differential analysis between somatic mutation and germline variation profiles reveals cancer-related genes

A major aim of cancer genomics is to pinpoint which somatically mutated genes are involved in tumor initiation and progression. We introduce a new framework for uncovering cancer genes, differential mutation analysis, which compares the mutational profiles of genes across cancer genomes with their natural germline variation across healthy individuals. We present DiffMut, a fast and simple approach for differential mutational analysis, and demonstrate that it is more effective in discovering cancer genes than considerably more sophisticated approaches. We conclude that germline variation across healthy human genomes provides a powerful means for characterizing somatic mutation frequency and identifying cancer driver genes. DiffMut is available at https://github.com/Singh-Lab/Differential-Mutation-Analysis.

GATA3 and TRPS1 are distinct biomarkers and prognostic factors in breast cancer: database mining for GATA family members in malignancies

GATA transcription factors are zinc finger DNA binding proteins that activate transcription during development and cell differentiation. To date, 7 members of GATA family have been reported. However, the expression patterns and the exact roles of distinct GATA family members contributing to tumorigenesis and progression of breast cancer (BC) remain to be elucidated. Here, we studied the expression of GATA transcripts in a variety of tumor types compared with the normal controls using the ONCOMINE and GOBO databases, along with their corresponding expression profiles in an array of cancer cell lines through CCLE analysis. Based on Kaplan-Meier plotter, we further investigated the prognostic values of GATA members specifically high expressed in BC patients. It was found that, when compared with normal tissues, GATA3 and TRPS1 were distinctly high expressed in BC patients among all GATA members. GATA3 expression was significantly associated with ESR1, while TRPS1 was correlated with ERBB2. In survival analysis, GATA3 and TRPS1 mRNA high expressions were correlated to better survival in BC patients, and TRPS1 high expression was significantly associated with longer RFS in patients who have received chemotherapy. These results suggest that GATA3 and TRPS1 are distinct biomarkers and essential prognostic factors for breast cancer.

A pilot study on the prevalence of DNA palindromes in breast cancer genomes

Background

DNA palindromes are a unique pattern of repeat sequences that are present in the human genome. It consists of a sequence of nucleotides in which the second half is the complement of the first half but appearing in reverse order. These palindromic sequences may have a significant role in DNA replication, transcription and gene regulation processes. They occur frequently in human cancers by clustering at specific locations of the genome that undergo gene amplification and tumorigenesis. Moreover, some studies showed that palindromes are clustered in amplified regions of breast cancer genomes especially in chromosomes (chr) 8 and 11. With the large number of personal genomes and cancer genomes becoming available, it is now possible to study their association to diseases using computational methods. Here, we conducted a pilot study on chromosomes 8 and 11 of cancer genomes to identify computationally the differentially occurring palindromes.

Methods

We processed 69 breast cancer genomes from The Cancer Genome Atlas including serum-normal and tumor genomes, and 1000 Genomes to serve as control group. The Biological Language Modelling Toolkit (BLMT) computes palindromes in whole genomes. We developed a computational pipeline integrating BLMT to compute and compare prevalence of palindromes in personal genomes.

Results

We carried out a pilot study on chr 8 and chr 11 taking into account single nucleotide polymorphisms, insertions and deletions. Of all the palindromes that showed any variation in cancer genomes, 38% of what were near breast cancer genes happened to be the most differentiated palindromes in tumor (i.e. they ranked among the top 25% by our heuristic measure).

Conclusions

These results will shed light on the prevalence of palindromes in oncogenes and the mutations that are present in the palindromic regions that could contribute to genomic rearrangements, and breast cancer progression.

The role of HDAC2 in chromatin remodelling and response to chemotherapy in ovarian cancer

Chromatin undergoes structural changes in response to extracellular and environmental signals. We observed changes in nuclear morphology in cancer tissue biopsied after chemotherapy and hypothesised that these DNA damage-induced changes are mediated by histone deacetylases (HDACs). Nuclear morphological changes in cell lines (PE01 and PE04 models) and a xenograft model (OV1002) were measured in response to platinum chemotherapy by image analysis of nuclear texture. HDAC2 expression increased in PEO1 cells treated with cisplatin at 24h, which was accompanied by increased expression of heterochromatin protein 1 (HP1). HDAC2 and HP1 expression were also increased after carboplatin treatment in the OV1002 carboplatin-sensitive xenograft model but not in the insensitive HOX424 model. Expression of DNA damage response pathways (pBRCA1, γH2AX, pATM, pATR) showed time-dependent changes after cisplatin treatment. HDAC2 knockdown by siRNA reduced HP1 expression, induced DNA double strand breaks (DSB) measured by γH2AX, and interfered with the activation of DNA damage response induced by cisplatin. Furthermore, HDAC2 depletion affected γH2AX foci formation, cell cycle distribution, and apoptosis triggered by cisplatin, and was additive to the inhibitory effect of cisplatin in cell lines. By inhibiting expression of HDAC2, reversible alterations in chromatin patterns during cisplatin treatment were observed. These results demonstrate quantifiable alterations in nuclear morphology after chemotherapy, and implicate HDAC2 in higher order chromatin changes and cellular DNA damage responses in ovarian cancer cells in vitro and in vivo.

Nuclear cathepsin D enhances TRPS1 transcriptional repressor function to regulate cell cycle progression and transformation in human breast cancer cells

The lysosomal protease cathepsin D (Cath-D) is overproduced in breast cancer cells (BCC) and supports tumor growth and metastasis formation. Here, we describe the mechanism whereby Cath-D is accumulated in the nucleus of ERα-positive (ER⁺) BCC. We identified TRPS1 (tricho-rhino-phalangeal-syndrome 1), a repressor of GATA-mediated transcription, and BAT3 (Scythe/BAG6), a nucleo-cytoplasmic shuttling chaperone protein, as new Cath-D-interacting nuclear proteins. Cath-D binds to BAT3 in ER⁺ BCC and they partially co-localize at the surface of lysosomes and in the nucleus. BAT3 silencing inhibits Cath-D accumulation in the nucleus, indicating that Cath-D nuclear targeting is controlled by BAT3. Fully mature Cath-D also binds to full-length TRPS1 and they co-localize in the nucleus of ER⁺ BCC where they are associated with chromatin. Using the LexA-VP16 fusion co-activator reporter assay, we then show that Cath-D acts as a transcriptional repressor, independently of its catalytic activity. Moreover, microarray analysis of BCC in which Cath-D and/or TRPS1 expression were silenced indicated that Cath-D enhances TRPS1-mediated repression of several TRPS1-regulated genes implicated in carcinogenesis, including PTHrP, a canonical TRPS1 gene target. In addition, co-silencing of TRPS1 and Cath-D in BCC affects the transcription of cell cycle, proliferation and transformation genes, and impairs cell cycle progression and soft agar colony formation. These findings indicate that Cath-D acts as a nuclear transcriptional cofactor of TRPS1 to regulate ER⁺ BCC proliferation and transformation in a non-proteolytic manner.

Transcription Factor Trps1 Promotes Tubular Cell Proliferation after Ischemia-Reperfusion Injury through cAMP-Specific 3',5'-Cyclic Phosphodiesterase 4D and AKT

Trichorhinophalangeal 1 (Trps1) is a transcription factor essential for epithelial cell morphogenesis during kidney development, but the role of Trps1 in AKI induced by ischemia-reperfusion (I/R) remains unclear. Our study investigated Trps1 expression during kidney repair after acute I/R in rats and explored the molecular mechanisms by which Trps1 promotes renal tubular epithelial cell proliferation. Trps1 expression positively associated with the extent of renal repair after I/R injury. Compared with wild-type rats, rats with knockdown of Trps1 exhibited significantly delayed renal repair in the moderate I/R model, with lower GFR levels and more severe morphologic injury, whereas rats overexpressing Trps1 exhibited significantly accelerated renal repair after severe I/R injury. Additionally, knockdown of Trps1 inhibited and overexpression of Trps1 enhanced the proliferation of renal tubular epithelial cells in rats. Chromatin immunoprecipitation sequencing assays and RT-PCR revealed that Trps1 regulated cAMP-specific 3',5'-cyclic phosphodiesterase 4D (Pde4d) expression. Knockdown of Trps1 decreased the renal protein expression of Pde4d and phosphorylated Akt in rats, and dual luciferase analysis showed that Trps1 directly activated Pde4d transcription. Furthermore, knockdown of Pde4d or treatment with the phosphatidylinositol 3 kinase inhibitor wortmannin significantly inhibited Trps1-induced tubular cell proliferation in vitro Trps1 may promote tubular cell proliferation through the Pde4d/phosphatidylinositol 3 kinase/AKT signaling pathway, suggesting Trps1 as a potential therapeutic target for kidney repair after I/R injury.

Overexpression of Trps1 contributes to tumor angiogenesis and poor prognosis of human osteosarcoma

Background

Trichorhinophalangeal syndrome 1 (Trps1) gene is a member of GATA transcription factor family and has an important function in tumorigenesis and progression. However, there are rare studies on its roles in carcinogenesis and prognostic significance in human osteosarcoma.

Methods

The expression of Trps1 was detected by immunohistochemistry, and MVD was evaluated to determine the amounts of microvessels by counting CD31-positive endothelial cells.

Results

Of the 74 cases that underwent study, Trps1-positive cases were 24. And it was associated with MVD significantly (P = 0.008). The data also exhibited more cases of remote metastasis (P = 0.013) and higher Enneking stage (P = 0.017) in Trps1-positive group compared to Trps1-negative group. Univariate analysis revealed that distant metastasis, MVD and Trps1 expression were associated with a lower 3-year overall survival rate and disease-free survival rate (P = 0.003, and P = 0.012 respectively). Furthermore, Trps1 and distant metastasis retained their significant prognostic effects on patients survival rate by multivariate analysis (P < 0.05).

Conclusions

Trps1 plays a crucial role in osteosarcoma angiogenesis, metastasis and clinical surgical stage. Trps1 can be a novel promising prognostic marker and therapeutic target, and antiangiogenic therapy which targets Trps1 molecule in patients with osteosarcoma may lead to improved prognosis and longer-term survival.

DEOD: uncovering dominant effects of cancer-driver genes based on a partial covariance selection method

MOTIVATION

The generation of a large volume of cancer genomes has allowed us to identify disease-related alterations more accurately, which is expected to enhance our understanding regarding the mechanism of cancer development. With genomic alterations detected, one challenge is to pinpoint cancer-driver genes that cause functional abnormalities.

RESULTS

Here, we propose a method for uncovering the dominant effects of cancer-driver genes (DEOD) based on a partial covariance selection approach. Inspired by a convex optimization technique, it estimates the dominant effects of candidate cancer-driver genes on the expression level changes of their target genes. It constructs a gene network as a directed-weighted graph by integrating DNA copy numbers, single nucleotide mutations and gene expressions from matched tumor samples, and estimates partial covariances between driver genes and their target genes. Then, a scoring function to measure the cancer-driver score for each gene is applied. To test the performance of DEOD, a novel scheme is designed for simulating conditional multivariate normal variables (targets and free genes) given a group of variables (driver genes). When we applied the DEOD method to both the simulated data and breast cancer data, DEOD successfully uncovered driver variables in the simulation data, and identified well-known oncogenes in breast cancer. In addition, two highly ranked genes by DEOD were related to survival time. The copy number amplifications of MYC (8q24.21) and TRPS1 (8q23.3) were closely related to the survival time with P-values = 0.00246 and 0.00092, respectively. The results demonstrate that DEOD can efficiently uncover cancer-driver genes.