Comprehensive gene expression analysis reveals multiple signal pathways associated with prostate cancer

Prostate cancer (PC) depends on androgenic signaling for growth and survival. To data, the exact molecular mechanism of hormone controlling proliferation and tumorigenesis in the PC remains unclear. Therefore, in this study, we explored the differentially expressed genes (DEGs) and identified featured genes related to hormone stimulus from PC. Two sets of gene expression data, including PC and normal control sample, were downloaded from Gene Expression Omnibus (GEO) database. The t-test was used to identify DEGs between PC and controls. Gene ontology (GO) functional annotation was applied to analyze the function of DEGs and screen hormone-related DEGs. Then these hormone-related DEGs were further analyzed in constructed cancer network and Human Protein Reference Database to screen important signaling pathways they participated in. A total of 912 DEGs were obtained which included 326 up-regulated genes and 586 down-regulated genes. GO functional enrichment analysis identified 50 hormone-related DEGs associated with PC. After pathway and PPI network analysis, we found these hormone-related DEGs participated in several important signaling pathways including TGF-β (TGFB2, TGFB3 and TGFBR2), MAPK (TGFB2, TGFB3 and TGFBR2), insulin (PIK3R3, SHC1 and EIF4EBP1), and p53 signaling pathways (CCND2 and CDKN1A). In addition, a total of five hormone-related DEGs (SHC1, CAV1, RXRA, CDKN1A and SRF) were located in the center of PPI network and 12 hormone-related DEGs formed six protein modules. These important signal pathways and hormone-related DEGs may provide potential therapeutic targets for PC.     

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.     

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.     

应用GSEA和WGCNA方法分析细菌性败血症宿主关键差异基因及意义

【目的】采用生物信息学方法分析公共数据库来源的细菌性败血症患者全血转录组学表达谱,探讨细菌败血症相关的宿主关键差异基因及意义。【方法】基于GEO数据库中GSE80496和GSE72829全血转录组基因数据集,采用GEO2R、基因集富集分析(GSEA)联用加权基因共表达网络分析(WGCNA)筛选细菌性败血症患者相比健康人群显著改变的差异基因,通过R软件对交集基因进行GO功能分析和KEGG富集分析。同时,通过String 11.0和Cytoscape分析枢纽基因,验证枢纽基因在数据集GSE72809(Health组52例,Definedsepsis组52例)全血标本中的表达情况,并探讨婴儿性别、月(胎)龄、出生体重、是否接触抗生素等因素与靶基因表达谱间的关系。【结果】分析GSE80496和GSE72829数据集分别筛选得到932个基因和319个基因,联合WGCNA枢纽模块交集得到与细菌性败血症发病相关的10个枢纽基因(MMP9、ITGAM、CSTD、GAPDH、PGLYRP1、FOLR3、OSCAR、TLR5、IL1RN和TIMP1);GSEA分析获得关键通路(氨基酸糖类-核糖代谢、PPAR信号通路、聚糖生物合成通路、自噬调控通路、补体、凝血因子级联反应、尼古丁和烟酰胺代谢、不饱和脂肪酸生物合成和阿尔兹海默症通路)及生物学过程(类固醇激素分泌、腺苷酸环化酶的激活、细胞外基质降解和金属离子运输)。【结论】本项研究通过GEO2R、GSEA联用WGCNA分析,筛选出与细菌性败血症发病相关的2个枢纽模块、10个枢纽基因以及一些关键信号通路和生物学过程,可为后续深入研究细菌性败血症致病机制奠定理论依据。     

smc5基因敲除斑马鱼肝脏转录组学分析

摘要 目的：探讨smc5基因敲除对斑马鱼肝脏基因表达谱的影响,进一步明确smc5突变对斑马鱼代谢的影响。方法：用CRISPR/Cas9技术构建smc5基因敲除斑马鱼模型,取3个月的smc5^-/-和野生型斑马鱼肝脏进行转录组测序,创建基因表达谱文库,观察smc5基因敲除后斑马鱼肝脏基因表达谱的变化,将筛选出的差异表达基因进行功能富集,并运用荧光定量PCR对KEGG通路中显著的差异表达基因进行验证。结果：成功构建出7号外显子上2碱基缺失造成移码突变的smc5基因敲除斑马鱼模型。RNA-seq发现smc5^-/-斑马鱼的肝脏基因表达谱变化显著,包含p53的多个通路激活,如细胞周期和凋亡。糖酵解、脂肪酸降解与代谢、丙酮酸代谢等相关通路显著下调。荧光定量PCR结果与RNA-seq结果一致。结论：smc5基因敲除下调斑马鱼肝脏糖脂代谢。本研究结果为进一步研究SMC5基因在糖脂代谢调控中的潜在机制奠定基础。     

沙葱萤叶甲卵滞育的转录组学分析

【目的】建立沙葱萤叶甲Galeruca daurica滞育卵转录组数据库,挖掘卵滞育相关的基因以及代谢和信号通路,在转录组水平探讨卵滞育的分子机制。【方法】采用Illumina NovaSeq6000高通量测序平台对沙葱萤叶甲滞育卵与解除滞育卵进行转录组测序,并进行生物信息学分析;利用DESeq软件分析沙葱萤叶甲滞育卵与解除滞育卵中的差异表达基因,对差异表达基因进行KEGG通路富集分析;利用qRT PCR技术对10个差异表达基因的表达模式进行验证。【结果】基于沙葱萤叶甲滞育卵与解除滞育卵转录组测序结果,共获得53 389个unigene,其中差异表达基因2 145个,24个差异表达基因与保幼激素信号及脂肪酸生物合成和降解相关。与解除滞育卵相比,滞育卵转录组中1 297个基因上调表达,富集于124条KEGG通路,其中核糖体通路显著富集;848个基因下调表达,富集于73条KEGG通路,其中MAPK信号通路和糖胺聚糖生物合成通路显著富集。qRT-PCR结果表明,随机选取的10个差异表达基因的表达趋势与RNA-Seq转录组测序结果完全一致。【结论】保幼激素,脂肪酸生物合成和降解,核糖体,MAPK信号及糖胺聚糖生物合成等通路可能在沙葱萤叶甲卵滞育调控中起着重要的作用。     

Identification of Key Genes and Pathways Involved in the Heterogeneity of Intrinsic Growth Ability Between Neurons After Spinal Cord Injury in Adult Zebrafish

In the adult central nervous system (CNS), axon regeneration is a major hurdle for functional recovery after trauma. The intrinsic growth potential of an injured axon varies widely between neurons. The underlying molecular mechanisms of such heterogeneity are largely unclear. In the present study, the adult zebrafish dataset GSE56842 were downloaded. Differentially expressed genes (DEGs) were sorted and deeply analyzed by bioinformatics methods. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of DEGs were performed with the DAVID. A DEGs-associated protein–protein interaction network was constructed from the STRING database and visualized with Cytoscape software. In total, 621 DEGs were identified. GO analysis showed that the biological processes of DEGs focused mainly on the Notch signaling pathway, cell differentiation and positive regulation of neuron differentiation. The molecular functions mainly included calcium-transporting ATPase activity and calcium ion binding and structural constituents of the cytoskeleton. The cellular components included the plasma membrane, spectrin, and cytoplasmic and membrane-bound vesicles. KEGG pathway analysis showed that these DEGs were mainly involved in the metabolic pathway and Notch signaling pathway, and subnetworks revealed that genes within modules were involved in the metabolic pathway, Wnt signaling pathway, and calcium signaling pathway. This study identified DEG candidate genes and pathways involved in the heterogeneity of the intrinsic growth ability between neurons after spinal cord injury in adult zebrafish, which could facilitate our understanding of the molecular mechanisms underlying axon regeneration, and these candidate genes and pathways could be therapeutic targets for the treatment of CNS injury.

     

2型猪链球菌脑膜炎小鼠模型的构建及其转录组学分析

【背景】2型猪链球菌(Streptococcus suis serotype 2, S. suis 2)可感染宿主引起严重的脑膜炎,对养猪业和人类公共卫生安全构成重大威胁。【目的】构建S. suis 2感染小鼠脑膜炎模型,并对其脑组织进行转录组学分析,为揭示S.suis2感染宿主后引起脑膜炎的分子机制和发现潜在的治疗靶点提供理论依据。【方法】采用S. suis 2感染小鼠,并对其脑组织进行病理组织学分析确认构建脑膜炎小鼠后,对其脑组织进行转录组学分析,对比S.suis2感染和未感染小鼠的差异表达基因,并对差异表达基因进行基因本体论(geneontology,GO)功能、京都基因和基因组百科全书(Kyoto encyclopedia of genes and genomes, KEGG)通路富集和韦恩分析。【结果】脑病理组织学分析结果显示,S. suis 2感染的小鼠脑膜中有大量的炎症细胞浸润,并且血管周围出现“袖套”现象,并能从感染小鼠的组织器官中再分离出攻毒的S. suis 2菌株,结果证明构建了S. suis 2感染脑膜炎小鼠模型。转录组学分析结果表明,感染S.suis2与未感染的...     

The Imprinted Gene PEG3 Inhibits Wnt Signaling and Regulates Glioma Growth

The imprinted gene PEG3 confers parenting and sexual behaviors, alters growth and development, and regulates apoptosis. However, a molecular mechanism that can account for the diverse functions of Peg3/Pw1 is not known. To elucidate Peg3-regulated pathways, we performed a functional screen in zebrafish. Enforced overexpression of PEG3 mRNA during zebrafish embryogenesis decreased β-catenin protein expression and inhibited Wnt-dependent tail development. Peg3/Pw1 also inhibited Wnt signaling in human cells by binding to β-catenin and promoting its degradation via a p53/Siah1-dependent, GSK3β-independent proteasomal pathway. The inhibition of the Wnt pathway by Peg3/Pw1 suggested a role in tumor suppression. Hypermethylation of the PEG3 promoter in primary human gliomas led to a loss of imprinting and decreased PEG3 mRNA expression that correlated with tumor grade. The decrease in Peg3/Pw1 protein expression increased β-catenin, promoted proliferation, and inhibited p53-dependent apoptosis in human CD133⁺ glioma stem cells. Thus, mammalian imprinting utilizes Peg3/Pw1 to co-opt the Wnt pathway, thereby regulating development and glioma growth.     

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.     

Bioinformatics identification of the methylerythritol phosphate pathway associated genes in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis> with <Emphasis Type="Italic">ceh1</Emphasis> mutant

The methylerythritol phosphate (MEP) pathway for the production of isoprenoids is recently discovered. The current study aimed to identify MEP pathway disorder-related molecular mechanisms and potential genes in Arabidopsis thaliana. Microarray data (GSE61675) obtained from ceh1 mutant plants and corresponding parental lines were retrieved from Gene Expression Omnibus (GEO) database and were applied for differentially expressed genes (DEGs) screening. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of DEGs were performed. Protein-protein interaction (PPI) network was then constructed and displayed by Cytoscape software. Total 762 DEGs including 620 up-regulated and 142 down-regulated genes were screened. In addition, a great many of DEGs were mainly involved in biosynthesis and metabolism-related pathways, such as stilbenoid, diarylheptanoid, and gingerol biosynthesis, and biosynthesis of terpenoids and steroids. Moreover, a PPI network contained 90 down-regulated genes and 497 up-regulated genes were obtained. Up-regulated DEGs including glutaredoxin (GRX480, cytochrome BC1 synthase (BCS1, syntaxin of plants 121 (SYP121) and A. thaliana MAP kinase 11 (ATMPK11) with higher degree in this network were hub nodes. Pathways including stilbenoid, diarylheptanoid, and gingerol biosynthesis obtained in our study were consistent with previous studies. Importantly, GRX480, BCS1 and ATMPK11 could have close interactions with the MEP pathway and may play important roles in the biosynthesis of isoprenoids.     

A PTCH1 Homolog Transcriptionally Activated by p53 Suppresses Hedgehog Signaling

The p53-mediated responses to DNA damage and the Hedgehog (Hh) signaling pathway are each recurrently dysregulated in many types of human cancer. Here we describe PTCH53, a p53 target gene that is homologous to the tumor suppressor gene PTCH1 and can function as a repressor of Hh pathway activation. PTCH53 (previously designated PTCHD4) was highly responsive to p53 in vitro and was among a small number of genes that were consistently expressed at reduced levels in diverse TP53 mutant cell lines and human tumors. Increased expression of PTCH53 inhibited canonical Hh signaling by the G protein-coupled receptor SMO. PTCH53 thus delineates a novel, inducible pathway by which p53 can repress tumorigenic Hh signals.     

Analysis of the potential pathways and target genes in spinal cord injury using bioinformatics methods

Spinal cord injury (SCI) remains to be the most devastating type of trauma for patients because of long lasting disability and limited response to the acute drug administration and efforts at rehabilitation. With the purpose to identify potential targets for SCI treatment and to gain more insights into the mechanisms of SCI, the microarray data of GSE2270, including 119 raphe magnus (RM) samples and 125 sensorimotor cortex (SMTC) samples, was downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were screened in RM group and SMTC group compared with their corresponding controls, respectively. A protein–protein interaction (PPI) network was constructed based on the common DEGs identified in both RM group and SMTC group. Gene ontology (GO) and pathway enrichment analyses of the overlapping DEGs were performed. Furthermore, the common DEGs enriched in each pathway were analyzed to identify significant regulatory elements. Totally, 173 overlapping DEGs (130 up-regulated and 43 down-regulated) were identified in both RM and SMTC samples. These overlapping DEGs were enriched in different GO terms. Pathway enrichment analysis revealed that DEGs were mainly related to inflammation and immunity. CD68 molecule (CD68) was a hub protein in the PPI network. Moreover, the regulatory network showed that ras-related C3 botulinum toxin substrate 2 (RAC2), CD44 molecule (CD44), and actin related protein 2/3 complex (ARPC1B) were hub genes. RAC2, CD44, and ARPC1B may be significantly involved in the pathogenesis of SCI by participating significant pathways such as extracellular matrix-receptor signaling pathway and Toll-like receptor signaling pathway.     

miR-731 modulates the zebrafish heart morphogenesis via targeting Calcineurin/Nfatc3a pathway

BackgroundZebrafish miR-731 is orthologous of human miR-425, which has been demonstrated to have cardio-protective roles by a variety of mechanisms. The miR-731 morphants show pericardium enlargement, and many DEGs (differentially expressed genes) are enriched in ‘Cardiac muscle contraction’ and ‘Calcium signaling pathway’, implying that miR-731 plays a potential role in heart function and development. However，the in vivo physiological role of miR-731 in the heart needs to be fully defined.MethodsZebrafish miR-731 morphants were generated by morpholino knockdown, and miR-731 knockout zebrafish was generated by CRISRP/Cas9. We observed cardiac morphogenesis based on whole-mount in situ hybridization. Furthermore, RNA-seq and qRT-PCR were used to elucidate the molecular mechanism and analyze the gene expression. Double luciferase verification and Western blot were used to verify the target gene.ResultsThe depletion of miR-731 in zebrafish embryos caused the deficiency of cardiac development and function, which was associated with reduced heart rate, ventricular enlargement and heart looping disorder. In addition, mechanistic study demonstrated that Calcineurin/Nfatc3a signaling involved in miR-731 depletion induced abnormal cardiac function and developmental defects.ConclusionMiR-731 regulates cardiac function and morphogenesis through Calcineurin/Nfatc3a signaling.General significanceOur studies highlight the potential importance of miR-731 in cardiac development.