Signatures of co-deregulated genes and their transcriptional regulators in colorectal cancer

Abstract The deregulated genes in colorectal cancer (CRC) vary significantly across different studies. Thus, a systems biology approach is needed to identify the co-deregulated genes (co-DEGs), explore their molecular networks, and spot the major hub proteins within these networks. We reanalyzed 19 GEO gene expression profiles to identify and annotate CRC versus normal signatures, single-gene perturbation, and single-drug perturbation signatures. We identified the co-DEGs across different studies, their upstream regulating kinases and transcription factors (TFs). Connectivity Map was used to identify likely repurposing drugs against CRC within each group. The functional changes of the co-upregulated genes in the first category were mainly associated with negative regulation of transforming growth factor β production and glomerular epithelial cell differentiation; whereas the co-downregulated genes were enriched in cotranslational protein targeting to the membrane. We identified 17 hub proteins across the co-upregulated genes and 18 hub proteins across the co-downregulated genes, composed of well-known TFs (MYC, TCF3, PML) and kinases (CSNK2A1, CDK1/4, MAPK14), and validated most of them using GEPIA2 and HPA, but also through two signature gene lists composed of the co-up and co-downregulated genes. We further identified a list of repurposing drugs that can potentially target the co-DEGs in CRC, including camptothecin, neostigmine bromide, emetine, remoxipride, cephaeline, thioridazine, and omeprazole. Similar analyses were performed in the co-DEG signatures in single-gene or drug perturbation experiments in CRC. MYC, PML, CDKs, CSNK2A1, and MAPKs were common hub proteins among all studies. Overall, we identified the critical genes in CRC and we propose repurposing drugs that could be used against them.

Download Full-text

Comprehensive Clinical and Molecular Characterization of KRASG12C-Mutant Colorectal Cancer

JCO Precision Oncology ◽

10.1200/po.20.00256 ◽

2021 ◽

pp. 613-621

Author(s):

Jason T. Henry ◽

Oluwadara Coker ◽

Saikat Chowdhury ◽

John Paul Shen ◽

Van K. Morris ◽

...

Keyword(s):

Colorectal Cancer ◽

Clinical Trials ◽

De Novo ◽

Expression Profiles ◽

Single Gene ◽

Gene Expression Profiles ◽

Progression Free Survival ◽

Cancer Center ◽

Data Set ◽

External Data

PURPOSE KRAS p.G12C mutations occur in approximately 3% of metastatic colorectal cancers (mCRC). Recently, two allosteric inhibitors of KRAS p.G12C have demonstrated activity in early phase clinical trials. There are no robust studies examining the behavior of this newly targetable population. METHODS We queried the MD Anderson Cancer Center data set for patients with colorectal cancer who harbored KRAS p.G12C mutations between January 2003 and September 2019. Patients were analyzed for clinical characteristics, overall survival (OS), and progression-free survival (PFS) and compared against KRAS nonG12C. Next, we analyzed several internal and external data sets to assess immune signatures, gene expression profiles, hypermethylation, co-occurring mutations, and proteomics. RESULTS Among the 4,632 patients with comprehensive molecular profiling, 134 (2.9%) were found to have KRAS p.G12C mutations. An additional 53 patients with single gene sequencing were included in clinical data but excluded from prevalence analysis allowing for 187 total patients. Sixty-five patients had de novo metastatic disease and received a median of two lines of chemotherapy without surgical intervention. For the first three lines of chemotherapy, the median PFS was 6.4 months (n = 65; 95% CI, 5.0 to 7.4 months), 3.9 months (n = 47; 95% CI, 2.9 to 5.9 months), and 3.0 months (n = 21; 95% CI, 2.0 to 3.4 months), respectively. KRAS p.G12C demonstrated higher rates of basal EGFR activation compared with KRAS nonG12C. When compared with an internal cohort of KRAS nonG12C, KRAS p.G12C patients had worse OS. CONCLUSION PFS is poor for patients with KRAS p.G12C metastatic colorectal cancer. OS was worse in KRAS p.G12C compared with KRAS nonG12C patients. Our data highlight the innate resistance to chemotherapy for KRAS p.G12C patients and serve as a historical comparator for future clinical trials.

Download Full-text

Use of gene set analysis (GSA) for molecular classification of responders and nonresponders to FOLFOX therapy in colorectal cancer (CRC).

Journal of Clinical Oncology ◽

10.1200/jco.2012.30.15_suppl.3628 ◽

2012 ◽

Vol 30 (15_suppl) ◽

pp. 3628-3628

Author(s):

Stephen Yip ◽

Maggie Chon U. Cheang ◽

Hagen F. Kennecke

Keyword(s):

Gene Expression ◽

Colorectal Cancer ◽

Metal Ion ◽

Expression Profiles ◽

Single Gene ◽

Chromosome 17 ◽

Aminopeptidase Activity ◽

P Value ◽

Hematopoietic Stem ◽

Gene Sets

3628 Background: Folinic acid (FOL) fluorouracil (F) and oxaliplatin (OX) chemotherapy is a commonly used therapy for CRC. Lacking in current literature are clinically relevant classifiers for potential responders. GSA technique is statistical method that detects significance of sets of genes, instead of examining a gene-by-gene basis. The objective of this study was to identify differential functionally annotated gene expression profiles associated with response to FOLFOX therapy in CRC tumors using gene-by-gene and GSA approaches. Methods: Genome wide expression profile data were collected on pre-treatment tumor tissues from patients with unresectable CRC receiving FOLFOX therapy (n = 83, Affymetrix HG U133A, GSE28702, Tsuji et al. BJC 2012). Gene expression was compared between responders (n = 42) and Non-responders (n = 41). GSA was conducted on 3272 curated gene sets from the Molecular Signatures Database (Subramanian, Tamayo et al. 2005, PNAS 102, 15545-15550) annotated by biological pathway, biochemical function and clinical behavior. Significant analysis of Microarray (SAM) and GSA (Tibshirani et al.) was done to identify gene sets associated with FOLFOX response. Results: Differential expressions of 23 genes were significantly associated with response, based on a single gene approach (p-value < 0.05). 13 of these were located on Chromosome 17 (p < 0.001). Among these, the top 5 ranked genes included NPEPPS, MBTD1, CEP44, LTA4H and CPNE4 which are involved in metal ion binding and aminopeptidase activity. GSA revealed only 44 out of 3272 gene sets were significantly associated with response, with a false discovery rate less than 25%. Increased expression of B-lymphocyte differentiation and Ras-signalling-related gene sets was associated with responders while mTOR signaling and hematopoietic stem cell-related genes set were associated with non-responders. Conclusions: Our data showed that differential biological pathways could be identified to predict response to FOLFOX therapy for CRC patients. Analysis may be useful to help define clinically relevant biologic subtypes among patients with metastatic colorectal cancer.

Download Full-text

Immunological and hematological outcomes following protracted low dose/low dose rate ionizing radiation and simulated microgravity

Scientific Reports ◽

10.1038/s41598-021-90439-5 ◽

2021 ◽

Vol 11 (1) ◽

Author(s):

Amber M. Paul ◽

Eliah G. Overbey ◽

Willian A. da Silveira ◽

Nathaniel Szewczyk ◽

Nina C. Nishiyama ◽

...

Keyword(s):

Low Dose ◽

Transforming Growth Factor ◽

Expression Profiles ◽

Transforming Growth Factor Β ◽

Hindlimb Unloading ◽

B Cell Differentiation ◽

Specific Regulation ◽

Β Receptor ◽

Gene Orthologs ◽

Whole Transcriptome

AbstractUsing a ground-based model to simulate spaceflight [21-days of single-housed, hindlimb unloading (HLU) combined with continuous low-dose gamma irradiation (LDR, total dose of 0.04 Gy)], an in-depth survey of the immune and hematological systems of mice at 7-days post-exposure was performed. Collected blood was profiled with a hematology analyzer and spleens were analyzed by whole transcriptome shotgun sequencing (RNA-sequencing). The results revealed negligible differences in immune differentials. However, hematological system analyses of whole blood indicated large disparities in red blood cell differentials and morphology, suggestive of anemia. Murine Reactome networks indicated majority of spleen cells displayed differentially expressed genes (DEG) involved in signal transduction, metabolism, cell cycle, chromatin organization, and DNA repair. Although immune differentials were not changed, DEG analysis of the spleen revealed expression profiles associated with inflammation and dysregulated immune function persist to 1-week post-simulated spaceflight. Additionally, specific regulation pathways associated with human blood disease gene orthologs, such as blood pressure regulation, transforming growth factor-β receptor signaling, and B cell differentiation were noted. Collectively, this study revealed differential immune and hematological outcomes 1-week post-simulated spaceflight conditions, suggesting recovery from spaceflight is an unremitting process.

Download Full-text

Identification and validation of a miRNA-related expression signature for tumor mutational burden in colorectal cancer

World Journal of Surgical Oncology ◽

10.1186/s12957-021-02137-1 ◽

2021 ◽

Vol 19 (1) ◽

Author(s):

Lijun Xu ◽

Qing Zheng

Keyword(s):

Colorectal Cancer ◽

Immune Response ◽

Validation Cohort ◽

Expression Profiles ◽

Functional Enrichment ◽

Immune Checkpoints ◽

Positive Correlation ◽

Expression Signature ◽

Mutational Burden ◽

Tumor Mutational Burden

Abstract Background Tumor mutational burden (TMB) is a promising predictor, which could stratify colorectal cancer (CRC) patients based on the response to immune checkpoint inhibitors (ICIs). MicroRNAs (miRNAs) act as the key regulators of anti-cancer immune response. However, the relationship between TMB and miRNA expression profiles is not elucidated in CRC. Methods Differentially expressed miRNAs (DE miRNAs) between the TMBhigh group and the TMBlow group were identified for the CRC cohort of the TCGA database. In the training cohort, a miRNA-related expression signature for predicting TMB level was developed by the least absolute shrinkage and selection operator (LASSO) method and tested with reference to its discrimination, calibration, and decision curve analysis (DCA) in the validation cohort. Functional enrichment analysis of these TMB-related miRNAs was performed. The correlation between this miRNA-related expression signature and three immune checkpoints was analyzed. Results Twenty-one out of 43 DE miRNAs were identified as TMB-related miRNAs, which were used to develop a miRNA-related expression signature. This TMB-related miRNA signature demonstrated great discrimination (AUCtest set = 0.970), satisfactory calibration (P > 0.05), and clinical utility in the validation cohort. Functional enrichment results revealed that these TMB-related miRNAs were mainly involved in biological processes associated with immune response and signaling pathways related with cancer. This miRNA-related expression signature showed a median positive correlation with PD-L1 (R = 0.47, P < 0.05) and CTLA4 (R = 0.39, P < 0.05) and a low positive correlation with PD-1 (R = 0.16, P < 0.05). Conclusion This study presents a miRNA-related expression signature which could stratify CRC patients with different TMB levels.

Download Full-text

Genome-wide analysis identifies critical DNA methylations within NTRKs genes in colorectal cancer

Journal of Translational Medicine ◽

10.1186/s12967-021-02740-6 ◽

2021 ◽

Vol 19 (1) ◽

Author(s):

Zijian Chen ◽

Zenghong Huang ◽

Yanxin Luo ◽

Qi Zou ◽

Liangliang Bai ◽

...

Keyword(s):

Colorectal Cancer ◽

Dna Methylation ◽

Promoter Methylation ◽

Expression Profiles ◽

Methylation Status ◽

Gene Promoter ◽

Normal Mucosa ◽

Suppressor Gene ◽

Survival Outcome ◽

Prognostic Models

Abstract Background Neurotrophic tropomyosin receptor kinases (NTRKs) are a gene family function as oncogene or tumor suppressor gene in distinct cancers. We aimed to investigate the methylation and expression profiles and prognostic value of NTRKs gene in colorectal cancer (CRC). Methods An analysis of DNA methylation and expression profiles in CRC patients was performed to explore the critical methylations within NTRKs genes. The methylation marker was validated in a retrospectively collected cohort of 229 CRC patients and tested in other tumor types from TCGA. DNA methylation status was determined by quantitative methylation-specific PCR (QMSP). Results The profiles in six CRC cohorts showed that NTRKs gene promoter was more frequently methylated in CRC compared to normal mucosa, which was associated with suppressed gene expression. We identified a specific methylated region within NTRK3 promoter targeted by cg27034819 and cg11525479 that best predicted survival outcome in CRC. NTRK3 promoter methylation showed independently predictive value for survival outcome in the validation cohort (P = 0.004, HR 2.688, 95% CI [1.355, 5.333]). Based on this, a nomogram predicting survival outcome was developed with a C-index of 0.705. Furthermore, the addition of NTRK3 promoter methylation improved the performance of currently-used prognostic model (AIC: 516.49 vs 513.91; LR: 39.06 vs 43.64, P = 0.032). Finally, NTRK3 promoter methylation also predicted survival in other tumors, including pancreatic cancer, glioblastoma and stomach adenocarcinoma. Conclusions This study highlights the essential value of NTRK3 methylation in prognostic evaluation and the potential to improve current prognostic models in CRC and other tumors.

Download Full-text

Identification of Key Candidate Genes Involved in the Progression of Idiopathic Pulmonary Fibrosis

Molecules ◽

10.3390/molecules26041123 ◽

2021 ◽

Vol 26 (4) ◽

pp. 1123

Author(s):

Yu Cui ◽

Jie Ji ◽

Jiwei Hou ◽

Yi Tan ◽

Xiaodong Han

Keyword(s):

Idiopathic Pulmonary Fibrosis ◽

Pulmonary Fibrosis ◽

Candidate Genes ◽

Expression Profiles ◽

Biological Processes ◽

Protein Protein Interaction ◽

Novel Treatments ◽

Kegg Pathway Analysis ◽

Cell Components ◽

Upregulated Genes

Idiopathic pulmonary fibrosis (IPF) is a lethal, agnogenic interstitial lung disease with limited therapeutic options. To investigate vital genes involved in the development of IPF, we integrated and compared four expression profiles (GSE110147, GSE53845, GSE24206, and GSE10667), including 87 IPF samples and 40 normal samples. By reanalyzing these datasets, we managed to identify 62 upregulated genes and 20 downregulated genes in IPF samples compared with normal samples. Differentially expressed genes (DEGs) were analyzed by gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis to illustrate relevant pathways of IPF, biological processes, molecular function, and cell components. The DEGs were then subjected to protein–protein interaction (PPI) for network analysis, serving to find 11 key candidate genes (ANXA3, STX11, THBS2, MMP1, MMP9, MMP7, MMP10, SPP1, COL1A1, ITGB8, IGF1). The result of RT-qPCR and immunohistochemical staining verified our finding as well. In summary, we identified 11 key candidate genes related to the process of IPF, which may contribute to novel treatments of IPF.

Download Full-text

Dynamic Changes of the Anthocyanin Biosynthesis Mechanism During the Development of Heading Chinese Cabbage (Brassica rapa L.) and Arabidopsis Under the Control of BrMYB2

Frontiers in Plant Science ◽

10.3389/fpls.2020.593766 ◽

2020 ◽

Vol 11 ◽

Author(s):

Qiong He ◽

Qianqian Lu ◽

Yuting He ◽

Yaxiu Wang ◽

Ninan Zhang ◽

...

Keyword(s):

Chinese Cabbage ◽

Anthocyanin Biosynthesis ◽

Expression Profiles ◽

Expression Patterns ◽

Anthocyanin Accumulation ◽

Total Anthocyanin ◽

Head Development ◽

Whole Plant ◽

Heading Chinese Cabbage ◽

Upregulated Genes

Chinese cabbage is an important vegetable mainly planted in Asian countries, and mining the molecular mechanism responsible for purple coloration in Brassica crops is fast becoming a research hotspot. In particular, the anthocyanin accumulation characteristic of purple heading Chinese cabbage, along with the plant’s growth and head developing, is still largely unknown. To elucidate the dynamic anthocyanin biosynthesis mechanism of Chinese cabbage during its development processes, here we investigated the expression profiles of 86 anthocyanin biosynthesis genes and corresponding anthocyanin accumulation characteristics of plants as they grew and their heads developed, between purple heading Chinese cabbage 11S91 and its breeding parents. Anthocyanin accumulation of 11S91 increased from the early head formation period onward, whereas the purple trait donor 95T2-5 constantly accumulated anthocyanin throughout its whole plant development. Increasing expression levels of BrMYB2 and BrTT8 together with the downregulation of BrMYBL2.1, BrMYBL2.2, and BrLBD39.1 occurred in both 11S91 and 95T2-5 plants during their growth, accompanied by the significantly continuous upregulation of a phenylpropanoid metabolic gene, BrPAL3.1; a series of early biosynthesis genes, such as BrCHSs, BrCHIs, BrF3Hs, and BrF3’H; as well as some key late biosynthesis genes, such as BrDFR1, BrANS1, BrUF3GT2, BrUF5GT, Br5MAT, and Brp-Cout; in addition to the transport genes BrGST1 and BrGST2. Dynamic expression profiles of these upregulated genes correlated well with the total anthocyanin contents during the processes of plant growth and leaf head development, and results supported by similar evidence for structural genes were also found in the BrMYB2 transgenic Arabidopsis. After intersubspecific hybridization breeding, the purple interior heading leaves of 11S91 inherited the partial purple phenotypes from 95T2-5 while the phenotypes of seedlings and heads were mainly acquired from white 94S17; comparatively in expression patterns of investigated anthocyanin biosynthesis genes, cotyledons of 11S91 might inherit the majority of genetic information from the white type parent, whereas the growth seedlings and developing heading tissues of 11S91 featured expression patterns of these genes more similar to 95T2-5. This comprehensive set of results provides new evidence for a better understanding of the anthocyanin biosynthesis mechanism and future breeding of new purple Brassica vegetables.

Download Full-text

Gene expression profiling of bovine endometrium during the oestrous cycle: detection of molecular pathways involved in functional changes

Journal of Molecular Endocrinology ◽

10.1677/jme.1.01799 ◽

2005 ◽

Vol 34 (3) ◽

pp. 889-908 ◽

Cited By ~ 92

Author(s):

S Bauersachs ◽

S E Ulbrich ◽

K Gross ◽

S E M Schmidt ◽

H H D Meyer ◽

...

Keyword(s):

Molecular Mechanisms ◽

Transforming Growth Factor ◽

Cdna Array ◽

Transforming Growth Factor Β ◽

Early Embryo ◽

Oestrous Cycle ◽

Cdna Libraries ◽

Functional Changes ◽

Retinoic Acid Signalling ◽

Subtracted Cdna Libraries

The endometrium plays a central role among the reproductive tissues in the context of early embryo–maternal communication and pregnancy. It undergoes typical changes during the sexual/oestrous cycle, which are regulated by the ovarian hormones progesterone and oestrogen. To identify the underlying molecular mechanisms we have performed the first holistic screen of transcriptome changes in bovine intercaruncular endometrium at two stages of the cycle – end of day 0 (late oestrus, low progesterone) and day 12 (dioestrus, high progesterone). A combination of subtracted cDNA libraries and cDNA array hybridisation revealed 133 genes showing at least a 2-fold change of their mRNA abundance, 65 with higher levels at oestrus and 68 at dioestrus. Interestingly, genes were identified which showed differential expression between different uterine sections as well. The most prominent example was the UTMP (uterine milk protein) mRNA, which was markedly upregulated in the cranial part of the ipsilateral uterine horn at oestrus. A Gene Ontology classification of the genes with known function characterised the oestrus time by elevated expression of genes, for example related to cell adhesion, cell motility and extracellular matrix and the dioestrus time by higher expression of mRNAs encoding for a variety of enzymes and transport proteins, in particular ion channels. Searching in pathway databases and literature data-mining revealed physiological processes and signalling cascades, e.g. the transforming growth factor-β signalling pathway and retinoic acid signalling, which are potentially involved in the regulation of changes of the endometrium during the oestrous cycle.

Download Full-text