Identification of Gene Expression Changes from Colitis to CRC in the Mouse CAC Model

A connection between colorectal carcinogenesis and inflammation is well known, but the underlying molecular mechanisms have not been elucidated. Chemically induced colitis-associated cancer (CAC) is an outstanding mouse model for studying the link between inflammation and cancer. Additionally, the CAC model is used for examining novel diagnostic, prognostic, and predictive markers for use in clinical practice. Here, a CAC model was established in less than 100 days using azoxymethane (AOM) with dextran sulfate sodium salt (DSS) in BALB/c mice. We examined the mRNA expression profiles of three groups: control untreated mice (K), DSS-induced chronic colitis mice (D), and AOM/DSS-induced CAC (AD) mice. We identified 6301 differentially expressed genes (DEGs) among the three groups, including 93 persistently upregulated genes and 139 persistently downregulated genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses revealed that the most persistent DEGs were significantly enriched in metabolic or inflammatory components in the tumor microenvironment. Furthermore, several associated DEGs were identified as potential DEGs by protein-protein interaction (PPI) network analysis. We selected 14 key genes from the DEGs and potential DEGs for further quantitative real-time PCR (qPCR) verification. Six persistently upregulated, 3 persistently downregulated DEGs, and the other 3 genes showed results consistent with the microarray data. We demonstrated the regulation of 12 key genes specifically involved in Wnt signaling, cytokine and cytokine receptor interactions, homeostasis, and tumor-associated metabolism during colitis-associated CRC. Our results suggest that a close relationship between metabolic and inflammatory mediators of the tumor microenvironment is present in CAC.     

Identification of key genes and pathways in benign prostatic hyperplasia

Benign prostatic hyperplasia (BPH) is one of the most common causes of lower urinary tract symptoms (LUTS) in elderly man. However, the underlying molecular mechanisms of BPH have not been completely elucidated. We identified the key genes and pathways by using analysis of Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were identified using edgeR. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed for the DEGs by Database for Annotation, Visualization and Integrated Discovery (DAVID) database and ConsensusPathDB, respectively. Then, protein–protein interaction (PPI) networks were established by the Search Tool for the Retrieval of Interacting Genes (STRING) database and visualized by Cytoscape software. Finally, we identified 660 DEGs ultimately including 268 upregulated genes and 392 downregulated genes. GO analysis revealed that DEGs were mainly enriched in extracellular exosome, identical protein binding, mitochondrial adenosine triphosphate (ATP) synthesis coupled proton transport, extracelluar matrix, focal adhesion, cytosol, Golgi apparatus, cytoplasm, protein binding, and Golgi membrane. Focal adhesion pathway, FoxO signaling pathway, and autophagy pathway were selected. Ubiquitin-conjugating enzyme E2 C (UBE2C), serine/threonine kinase (AKT1), mitogen-activated protein kinase 1 (MAPK1), cyclin B1 (CCNB1), polo-like kinase 1 (PLK1) were filtrated as the hub genes according to the degree of connectivity from the PPI network. The five hub genes including UBE2C, AKT1, MAPK1, CCNB1, PLK1 may play key roles in the pathogenesis of benign prostatic hyperplasia (BPH). Focal adhesion pathway, FoxO signaling pathway, and autophagy pathway may be crucial for the progression of BPH.     

Bioinformatics analysis reveals disturbance mechanism of MAPK signaling pathway and cell cycle in Glioblastoma multiforme

Background & objectives

To analyze the reversal gene pairs and identify featured reversal genes related to mitogen-activated protein kinases (MAPK) signaling pathway and cell cycle in Glioblastoma multiforme (GBM) to reveal its pathogenetic mechanism.

Methods

We downloaded the gene expression profile GSE4290 from the Gene Expression Omnibus database, including 81 gene chips of GBM and 23 gene chips of controls. The t test was used to analyze the DEGs (differentially expressed genes) between 23 normal and 81 GBM samples. Then some perturbing metabolic pathways, including MAPK (mitogen-activated protein kinases) and cell cycle signaling pathway, were extracted from KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway database. Cancer genes were obtained from the database of Cancer Gene Census. The reversal gene pairs between DEGs and cancer genes were further analyzed in MAPK and cell cycle signaling pathway.

Results

A total 8523 DEGs were obtained including 4090 up-regulated and 4433 down-regulated genes. Among them, ras-related protein rab-13(RAB13), neuroblastoma breakpoint family member 10 (NBPF10) and disks large homologue 4 (DLG4) were found to be involved in GBM for the first time. We obtained MAPK and cell cycle signaling pathways from KEGG database. By analyzing perturbing mechanism in these two pathways, we identified several reversal gene pairs, including NRAS (neuroblastoma RAS) and CDK2 (cyclin-dependent kinase 2), CCND1 (cyclin D1) and FGFR (fibroblast growth factor receptor). Further analysis showed that NRAS and CDK2 were positively related with GBM. However, FGFR2 and CCND1 were negatively related with GBM.

Interpretation & conclusions

These findings suggest that newly identified DEGs and featured reversal gene pairs participated in MAPK and cell cycle signaling pathway may provide a new therapeutic line of approach to GBM.     

Transcriptome analysis of the response of silkworm to drastic changes in ambient temperature

Bombyx mori is a poikilothermic insect and is economically important for silk production. Drastic changes in the ambient temperature have a negative impact on sericulture. However, the reason as to why high temperature is associated with the occurrence of diseases in silkworm and the response of silkworm to low temperature remain unclear and were the focus of the present study. Dazao silkworm exposed to 13 °C (DZ-13), 25 °C (DZ-25), and 37 °C (DZ-37) were used for RNA-seq analysis. There were 478 and 194 upregulated differentially expressed genes (DEGs) in DZ-13 and DZ-37 while 49 and 273 downregulated DEGs in DZ-13 and DZ-37, respectively. Eight DEGs were co-upregulated, in which seven genes were for heat shock proteins (Hsps), implying that Hsps play important roles in the tolerance of silkworm to high and low temperature. Gene ontology analysis revealed that the developmental process was downregulated in DZ-13. All the DEGs in the oxidative phosphorylation and insulin signaling pathways were upregulated in DZ-13. Several cuticular proteins and ATP synthesis-related genes were upregulated in DZ-13, suggesting that thickening of the cuticle and increase in the ATPase expression would help silkworms to protect themselves from low temperature-induced stress. Several immune-related genes, such as BmRel and BmSerpin-2, were downregulated in DZ-37, revealing that the resistance of silkworm is decreased under high temperature shock resulting in susceptibility to pathogens. Thus, the increase in the thermo-tolerance of silkworm should be related to the enhancement in the pathogen resistance.

     

Identification of key candidate genes and pathways in multiple myeloma by integrated bioinformatics analysis

Multiple myeloma (MM) is a common hematologic malignancy for which the underlying molecular mechanisms remain largely unclear. This study aimed to elucidate key candidate genes and pathways in MM by integrated bioinformatics analysis. Expression profiles GSE6477 and GSE47552 were obtained from the Gene Expression Omnibus database, and differentially expressed genes (DEGs) with p < .05 and [logFC] > 1 were identified. Functional enrichment, protein–protein interaction network construction and survival analyses were then performed. First, 51 upregulated and 78 downregulated DEGs shared between the two GSE datasets were identified. Second, functional enrichment analysis showed that these DEGs are mainly involved in the B cell receptor signaling pathway, hematopoietic cell lineage, and NF-kappa B pathway. Moreover, interrelation analysis of immune system processes showed enrichment of the downregulated DEGs mainly in B cell differentiation, positive regulation of monocyte chemotaxis and positive regulation of T cell proliferation. Finally, the correlation between DEG expression and survival in MM was evaluated using the PrognoScan database. In conclusion, we identified key candidate genes that affect the outcomes of patients with MM, and these genes might serve as potential therapeutic targets.     

Prognostic value of metabolic genes in lung adenocarcinoma via integrative analyses

BackgroundLung adenocarcinoma (LUAD) is the most common malignant lung tumor. Metabolic pathway reprogramming is an important hallmark of physiologic changes in cancers. However, the mechanisms through which these metabolic genes and pathways function in LUAD as well as their prognostic values have not been fully established.MethodsFour publicly available datasets from GEO and TCGA were used to identify differentially expressed genes (DEGs) in LUAD, which were then subjected to GO and KEGG pathway enrichment analysis. Associations between metabolic gene expressions with overall survival, tumor stage, TP53 mutation status, and infiltrated immune cells were investigated. Protein-protein interactions were evaluated using GeneMANIA and Metascape.ResultsBy integrating four public datasets, 247 DEGs were identified in LUAD. These DEGs were significantly enriched in regulation of chromosome segregation, centromeric region, and histone kinase activity GO terms, as well as in cell cycle, p53 signaling pathway, metabolic pathways, and other KEGG pathways. Elevated expressions of ten metabolic genes in LUAD were significantly associated with poor survival outcomes. These metabolic genes were highly expressed in more advanced tumor stage and TP53 mutated patients. Moreover, expression levels were significantly correlated with tumor-infiltrating immune cells. PPI interaction analysis revealed that the top 20 genes interacting with each metabolic gene were significantly enriched in DNA replication, response to radiation, and central carbon metabolism in cancer.ConclusionThis study elucidates on molecular changes in metabolic genes in LUAD, which may inform the development of genetically oriented diagnostic approaches and effective treatment options.     

High-dose dexamethasone injection disordered metabolism and multiple protein kinases expression in the mouse kidney

Glucocorticoids (GCs) have been widely used in clinical treatment as anti-inflammatory, anti-shock and immunosuppressive medicines. However, the effect of excessive GCs on immune response and metabolism of kidney remains unclear. Here, we profiled the gene expression of kidney from mice with high-dose dexamethasone (DEX) treatment. A total of 1193 differentially expressed genes (DEGs) were screened in DEX treatment group compared with the saline group, including 715 down- regulated and 478 up-regulated. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of these DEGs showed extracellular matrix (ECM)–receptor interaction, cell adhesion molecules signaling pathway were significantly enriched, and that the vast majority of DEGs were involved in monocarboxylic acid metabolism, leukocyte cell–cell adhesion and fatty acid metabolism. Gene set enrichment analysis (GSEA) revealed that DEGs were strongly associated with immune-response and cell adhesion gene sets, such as Fc γ R-mediated phagocytosis, leukocyte transendothelial migration, T-cell receptor signaling pathway, cell adhesion, ECM–receptor interaction and focal adhesion-associated pathways. KEGG pathway analysis of differentially expressed kinases (DEKs) showed T-cell receptor and forkhead box class O signaling pathway were enriched. Furthermore, we found multiple protein kinases expression were dysregulated greatly after dexamethasone treatment, including classical effector of GCs stimulation-serum and GC-regulated kinase. These protein kinases are involved in multiple signaling pathways in mice kidney, such as mitogen-activated protein kinase (MAPK) and phosphoinositide 3-kinase (PI3K)/Akt signaling pathway. We profiled the gene expression of the kidney from high-dose dexamethasone-treated mice and provided important information for further study the mechanism of side effects of GCs in clinical therapy.     

Knockdown of CTRP6 reduces the deposition of intramuscular and subcutaneous fat in pigs via different signaling pathways

The regulation of porcine subcutaneous (SC) and intramuscular (IM) fat deposition significantly affects pork quality and the lean meat percentage of the carcass, respectively. The adipokine C1q/tumor necrosis factor-related protein 6 (CTRP6), plays a significant role in regulating animal fat deposition. The purpose of this study was to understand the effects of CTRP6 gene knockdown in IM and SC adipocytes by RNA-seq analysis. A total of 1830 and 2936 differentially expressed genes (DEGs) were identified in SC and IM adipocytes, respectively. 844 were down- and 2092 were upregulated in SC adipocytes, while 648 were down- and 1182 were upregulated in IM adipocytes. Furthermore, 1778 DEGs were detected only in SC adipocytes, 672 DEGs only in IM adipocytes, and 1158 DEGs in both types of adipocytes. GO analysis indicated that DEGs involved in adipocyte differentiation were significantly enriched in both SC and IM adipocytes following treatment with CTRP6-siRNA. Moreover, KEGG pathway enrichment analysis revealed differences of metabolic regulation between IM and SC adipocytes. With CTRP6-silencing, the signaling pathways related to Ras and arachidonic acid metabolism were significantly enriched in IM adipocytes, while four other signaling pathways, encompassing the TNF, MAPK, p53 and adipokine pathway were specifically enriched in SC adipocytes. Interestingly, the effect of CTRP6-siRNA treatment was attenuated by the specific Ras activator ML-097 in IM adipocytes, while the specific p53 activator SJ-172550 had the corresponding effect in SC adipocytes. Altogether, we suggest that CTRP6 may be a differential regulator of the development and metabolism of IM and SC adipose tissues.     

Differences in potential key genes and pathways between primary and radiation-associated angiosarcoma of the breast

BackgroundAngiosarcoma of the breast is a high-grade malignant soft tissue tumor, it can be divided into primary and radiation-associated angiosarcoma(secondary). However, the differences between primary and secondary angiosarcomas in terms of pathogenesis, clinical behavior, early diagnosis biomarkers, genetic abnormalities, and therapeutic targets remain to be fully elucidated. At the same time, due to its rarity, most of current information relating to angiosarcoma is provided by case reports. Therefore, exploring the mechanisms of primary and secondary breast angiosarcoma have important value for the discovery of new biomarkers and research into potential therapeutic targets.MethodsThe differentially expressed genes (DEGs) between 36 cases of primary angiosarcoma and 54 cases of secondary angiosarcoma were screened. Then, the DEGs were used to gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Then, a protein-protein interaction (PPI) network was constructed using the STRING database.ResultsA total of 18 DEGs were identified, of which 13 were upregulated and 5 were downregulated in secondary breast angiosarcoma. The GO enrichment analysis showed that the DEGs were most enriched in metabolism, energy pathways, and protein metabolism in biological processes. The enriched signaling pathways of DEGs were the transforming growth factor-β (TGF-β), Wnt, Hippo and PI3K-Akt signaling pathways. Then, the PPI network was conducted and hub genes were identified and they were involved in thyroid hormone, Hippo and other signaling pathways.ConclusionThis study lay the foundation for the discovery of effective and reliable molecular biomarkers and essential therapeutic targets for these malignancies.     

Transcriptome Profiling of Barley Cultivar Hua 30 MDEC in Response to <i>Agrobacterium</i> Infection

Agrobacterium-mediated transformation has been widely used in plants. However, the mechanism in plant cells’ response to Agrobacterium infection was very complex. The mechanism of the determinants in host cell remains obscure, especially in barley, which is recalcitrant for Agrobacterium-mediated transformation. In the present study, microspore-derived embryogenic calli (MDEC) from barley elite cultivar were employed as unique subjects to characterize the mechanisms during the Agrobacterium infection process. Hua 30 MDEC can be successfully infected by Agrobacterium. RNA-sequencing at different infection points (0, 2, 6, 12, 24 hpi) was performed. The average expressional intensity of the whole genomics increased from 0 to 2 hpi, and then decreased subsequently. More upregulated than downregulated differentially expressed genes (DEGs) were counted at the same time. GO enrichment analysis showed that protein modification was significantly overrepresented in upregulated DEGs. Chromosome-related biological processes, gene expression and cellular metabolic processes were significantly overrepresented in downregulated DEGs. KEGG analysis showed that plant defense responses, phenylpropanoid biosynthesis and biosynthesis of amino acids were significantly enriched across the infection time course. Nine DEGs related to defense responses were identified. All DEGs were upregulated from 2 to 24 hpi. We speculate that these genes are possibly related to Agrobacterium infection. These findings will provide deep insights into the molecular events occurring during the process of Agrobacterium-mediated transformation.     

Identification and validation of key genes associated with non-small-cell lung cancer

Non-small-cell lung cancer (NSCLC) is one of the main causes of death induced by cancer globally. However, the molecular aberrations in NSCLC patients remain unclearly. In the present study, four messenger RNA microarray datasets (GSE18842, GSE40275, GSE43458, and GSE102287) were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) between NSCLC tissues and adjacent lung tissues were obtained from GEO2R and the overlapping DEGs were identified. Moreover, functional and pathway enrichment were performed by Funrich, while the protein–protein interaction (PPI) network construction were obtained from STRING and hub genes were visualized and identified by Cytoscape software. Furthermore, validation, overall survival (OS) and tumor staging analysis of selected hub genes were performed by GEPIA. A total of 367 DEGs (95 upregulated and 272 downregulated) were obtained through gene integration analysis. The PPI network consisted of 94 nodes and 1036 edges in the upregulated DEGs and 272 nodes and 464 edges in the downregulated DEGs, respectively. The PPI network identified 46 upregulated and 27 downregulated hub genes among the DEGs, and six (such as CENPE, NCAPH, MYH11, LRRK2, HSD17B6, and A2M) of that have not been identified to be associated with NSCLC so far. Moreover, the expression differences of the mentioned hub genes were consistent with that in lung adenocarcinoma and lung squamous cell carcinoma in the TCGA database. Further analysis showed that all the six hub genes were associated with tumor staging except MYH11, while only the upregulated DEG CENPE was associated with the worse OS of patients with NSCLC. In conclusion, the current study showed that CENPE, NCAPH, MYH11, LRRK2, HSD17B6, and A2M might be the key genes contributed to tumorigenesis or tumor progression in NSCLC, further functional study is needed to explore the involved mechanisms.