NS398 as a potential drug for autosomal-dominant polycystic kidney disease: Analysis using bioinformatics,and zebrafish and mouse models

Autosomal-dominant polycystic kidney disease (ADPKD) is characterized by uncontrolled renal cyst formation, and few treatment options are available. There are many parallels between ADPKD and clear-cell renal cell carcinoma (ccRCC); however, few studies have addressed the mechanisms linking them. In this study, we aimed to investigate their convergences and divergences based on bioinformatics and explore the potential of compounds commonly used in cancer research to be repurposed for ADPKD. We analysed gene expression datasets of ADPKD and ccRCC to identify the common and disease-specific differentially expressed genes (DEGs). We then mapped them to the Connectivity Map database to identify small molecular compounds with therapeutic potential. A total of 117 significant DEGs were identified, and enrichment analyses results revealed that they are mainly enriched in arachidonic acid metabolism, p53 signalling pathway and metabolic pathways. In addition, 127 ccRCC-specific up-regulated genes were identified as related to the survival of patients with cancer. We focused on the compound NS398 as it targeted DEGs and found that it inhibited the proliferation of Pkd1^−/− and 786-0 cells. Furthermore, its administration curbed cystogenesis in Pkd2 zebrafish and early-onset Pkd1-deficient mouse models. In conclusion, NS398 is a potential therapeutic agent for ADPKD.     

基于转录组的瑶药半枫荷根和叶的差异分析

为了探究半枫荷(Semiliquidambar cathayensis)根和叶的基因表达差异和关键活性成分合成通路中关键基因的表达规律,本研究对半枫荷的根和叶进行转录组测序和生物信息分析。59378个差异表达基因归类到在GO分类的3个大类中,主要与生物学过程有关(50.59%)。626个差异表达基因注释KOG数据库的24个分类中。81个差异表达基因参与苯丙烷类化合物的生物合成,110个差异表达基因参与黄酮类化合物的生物合成,211个差异表达基因参与萜类化合物的生物合成。本研究获得了半枫荷根和叶的转录组信息特征,为今后半枫荷基因功能鉴定、次生代谢途径解析及调控机制的研究提供依据。     

Ossification of the posterior longitudinal ligament related genes identification using microarray gene expression profiling and bioinformatics analysis

Ossification of the posterior longitudinal ligament (OPLL) is a kind of disease with physical barriers and neurological disorders. The objective of this study was to explore the differentially expressed genes (DEGs) in OPLL patient ligament cells and identify the target sites for the prevention and treatment of OPLL in clinic. Gene expression data GSE5464 was downloaded from Gene Expression Omnibus; then DEGs were screened by limma package in R language, and changed functions and pathways of OPLL cells compared to normal cells were identified by DAVID (The Database for Annotation, Visualization and Integrated Discovery); finally, an interaction network of DEGs was constructed by string. A total of 1536 DEGs were screened, with 31 down-regulated and 1505 up-regulated genes. Response to wounding function and Toll-like receptor signaling pathway may involve in the development of OPLL. Genes, such as PDGFB, PRDX2 may involve in OPLL through response to wounding function. Toll-like receptor signaling pathway enriched genes such as TLR1, TLR5, and TLR7 may involve in spine cord injury in OPLL. PIK3R1 was the hub gene in the network of DEGs with the highest degree; INSR was one of the most closely related genes of it. OPLL related genes screened by microarray gene expression profiling and bioinformatics analysis may be helpful for elucidating the mechanism of OPLL.     

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.

     

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.     

多球壳菌素介导线粒体功能调控酵母细胞衰老

多球壳菌素是前期鉴定出的一种抗衰老活性小分子,为了较全面理解多球壳菌素的抗衰机制,本研究利用模式生物芽殖酵母为材料,借助DNA Microarray技术、生物信息学手段并结合相关功能实验分析了多球壳菌素对对数期芽殖酵母细胞基因表达谱、基因本体聚类及相关信号通路的影响。结果显示,多球壳菌素对细胞转录组产生了显著影响,共导致1 648个基因的差异表达(FDR0.05,Fold Change1.5),其中843个基因显著上调,805个基因显著下调。进一步对基因本体聚类、信号通路分析及功能实验验证显示,线粒体相关功能和信号通路是多球壳菌素作用的一个主要靶点,本研究对多球壳菌素的进一步开发利用奠定了重要的理论基础。     

Exploring the correlation between deltamethrin stress and Keap1-Nrf2-ARE pathway from Drosophila melanogaster RNASeq data

Deltamethrin (DM) is widely used in a variety of pest control, resulting in serious drug resistance. Keap1-Nrf2-ARE is the antioxidant stress pathway. We identified 268 genes differentially expressed (DEGs) in Drosophila Kc cells treated with DM, including up-regulated 180 genes and down-regulated 88 genes compared with the control group (fold-change≥2, qValue≤0.001) by RNA-seq, they are mainly linked to metabolic process, stimulation response, immune system process. When the cells are treated with DM in the case of overexpression of the Keap1 gene, the cytochrome P450 family genes were significantly down-regulated, and some diseases-related genes and non-coding genes also changed. Our data shown that Keap1-Nrf2-ARE pathway may play an important role in DM stress, which will provide a new direction for studying the mechanism of insect resistance.     

Pathways related to PMA-differentiated THP1 human monocytic leukemia cells revealed by RNA-Seq

Previous analyses have reported that the human monocytic cell line THP1 can be differentiated into cells with macrophage-like characteristics by phorbol 12-myristate 13-acetate(PMA). However, little is known about the mechanism responsible for regulating this differentiation process. Here, we performed high-throughput RNA-Seq analysis to investigate the genes differently expressed in THP1 cells treated with and without PMA and examined those that may be responsible for the PMA-induced differentiation of monocytes into macrophages. We found 3,000 genes to be differentially expressed after PMA treatment. Gene ontology analysis revealed that genes related to cellular processes and regulation of biological processes were significantly enriched. KEGG analysis also demonstrated that the differentially expressed genes(DEGs) were significantly enriched in the PI3K/AKT signaling pathway and phagosome pathway. Importantly, we reveal an important role of the PI3K/AKT pathway in PMA-induced THP1 cell differentiation. The identified DEGs and pathways may facilitate further study of the detailed molecular mechanisms of THP1 differentiation. Thus, our results provide numerous potential therapeutic targets for modulation of the differentiation of this disease.     

Integrative analyses of genes associated with idiopathic pulmonary fibrosis

Idiopathic pulmonary fibrosis (IPF), characterized by irreversible scarring and progressive destruction of the lung tissue, is one of the most common types of idiopathic interstitial pneumonia worldwide. However, there are no reliable candidates for curative therapies. Hence, elucidation of the mechanisms of IPF genesis and exploration of potential biomarkers and prognostic indicators are essential for accurate diagnosis and treatment of IPF. Recently, efficient microarray and bioinformatics analyses have promoted an understanding of the molecular mechanisms of disease occurrence and development, which is necessary to explore genetic alternations and identify potential diagnostic biomarkers. However, high false-positive rates results have been observed based on single microarray datasets. In the current study, we performed a comprehensive analysis of the differential expression, biological functions, and interactions of IPF-related genes. Three publicly available microarray datasets including 54 IPF samples and 34 normal samples were integrated by performing gene set enrichment analysis and analyzing differentially expressed genes (DEGs). Our results identified 350 DEGs genetically associated with IPF. Gene ontology analyses revealed that the changes in the modules were mostly enriched in the positive regulation of smooth muscle cell proliferation, positive regulation of inflammatory responses, and the extracellular space. Kyoto encyclopedia of genes and genomes enrichment analysis of DEGs revealed that IPF involves the TNF signaling pathway, NOD-like receptor signaling pathway, and PPAR signaling pathway. To identify key genes related to IPF in the protein-protein interaction network, 20 hub genes were screened out with highest scores. Our results provided a framework for developing new pathological molecular networks related to specific diseases in silico.     

Screening of hub genes and prediction of putative drugs in arsenic-related bladder carcinoma: An in silico study

BackgroundEvidence showed that inorganic arsenic (iAs) can trigger malignant transformation in cells with complex mechanisms. Thus, we aimed to investigate the possible molecules, pathways and therapeutic drugs for iAs-induced bladder cancer (BC) by using bioinformatics approaches.MethodsMicroarray-based data were analyzed to screen the differentially expressed genes (DEGs) between iAs-related BC cells and controls. Then, the roles of DEGs were annotated and the hub genes were screened out by protein-protein interaction network. The key genes were further selected from the hub genes through an assessment of the prognostic values. Afterward, potential drugs were predicted by using CMAP analysis.ResultsAnalysis of a dataset (GSE90023) generated 21 upregulated and 47 downregulated DEGs, which were enriched in various signaling pathways. Among the DEGs, four hub genes including WNT7B, SFRP1, DNAJB2, and ATF3, were identified as the key genes because they might predict poor prognosis in BC patients. Lastly, Cantharidin was predicted to be a potential drug reversing iAs-induced malignant transformation in urinary epithelium cells.ConclusionThe present study found several hub genes involved in iAs-induced malignant transformation in urinary epithelium cells, and predicted several small agents for iAs toxicity prevention or therapy.     

Screening and identification of the core immune-related genes and immune cell infiltration in severe burns and sepsis

Severe burns often have a high mortality rate due to sepsis, but the genetic and immune crosstalk between them remains unclear. In the present study, the GSE77791 and GSE95233 datasets were analysed to identify immune-related differentially expressed genes (DEGs) involved in disease progression in both burns and sepsis. Subsequently, weighted gene coexpression network analysis (WGCNA), gene enrichment analysis, protein–protein interaction (PPI) network construction, immune cell infiltration analysis, core gene identification, coexpression network analysis and clinical correlation analysis were performed. A total of 282 common DEGs associated with burns and sepsis were identified. Kyoto Encyclopedia of Genes and Genomes pathway analysis identified the following enriched pathways in burns and sepsis: metabolic pathways; complement and coagulation cascades; legionellosis; starch and sucrose metabolism; and ferroptosis. Finally, six core DEGs were identified, namely, IL10, RETN, THBS1, FGF13, LCN2 and MMP9. Correlation analysis showed that some core DEGs were significantly associated with simultaneous dysregulation of immune cells. Of these, RETN upregulation was associated with a worse prognosis. The immune-related genes and dysregulated immune cells in severe burns and sepsis provide potential research directions for diagnosis and treatment.     

琥珀酸对酵母细胞转录组的影响

寻找抗衰老活性小分子并研究其作用机制是衰老药物学研究的重点和热点。本文报道了一种新的抗衰老活性小分子琥珀酸,发现琥珀酸可以显著延缓芽殖酵母细胞的衰老并增强细胞的压力抗性。随后,利用DNA Microarray技术及生物信息学手段较系统分析了琥珀酸处理对基因表达谱、基因本体聚类及相关信号通路的影响。结果显示,琥珀酸处理对细胞转录组产生了显著影响,共导致3 485个基因的差异表达(P 0.05),其中1 335个基因显著上调,2 150个基因显著下调。进一步对基因本体聚类及信号通路分析显示,线粒体及核糖体生物合成相关的分子功能、细胞组分、生物学过程和信号通路可能是琥珀酸作用的主要靶点,其他可能的作用靶点还包括蛋白酶体、细胞内吞、过氧化物酶体代谢及细胞自噬等。本研究为进一步阐明琥珀酸介导的寿命及压力调控机制提供了理论参考和研究线索。     

小儿脓毒症休克相关基因的筛选及网络构建

为探究脓毒症休克与SIRS的差异表达基因及网络的构建,筛选潜在的核心基因,从GEO数据库下载相关基因表达谱GSE26378,数据分为脓毒症休克与SIRS各29个样本,通过在线软件GCBI对其进行标准化及差异基因筛选;对差异基因进行GO分析;基于KEGG进行功能通路分析以及基因信号网络分析;差异基因共表达网络分析。结果表明:两组中总共有1 456个基因被识别为差异基因(P0.05),与SIRS组相比,脓毒症休克组中有条859条下调基因,597条上调基因。GO功能富集分析显示差异基因主要参与了细胞周期、细胞免疫、细胞代谢。KEGG功能通路分析显示差异基因主要参与了MAPK信号通路、P53信号通路、wnt信号通路、细胞凋亡信号通路,细胞周期受体信号通路等。共表达分析发现基因CCNB1、NUSAP1、OIP5、SHCBP1、ZWINT、TOP2A、DLGAP5等位于网络中央部位,而基因信号网络分析发现基因PLCB1、PIK3CA、STAT3、CAMK2D、PRKCB、CREB1位于网络核心。基因芯片分析有助于发现脓毒症休克与SIRS患儿外周血单核细胞在转录组学上的改变,而生物信息学网络分析有助于发现潜在的靶点。     

基于转录组的肝豆状核变性调控网络的构建和分析

本研究旨在利用生物信息学方法构建经铜诱导的ATP7B基因敲除HepG2细胞系的转录调控网络。探讨关键转录因子在肝豆状核变性发生、发展中的潜在作用机制。收集公共基因表达数据库(gene expression omnibus, GEO)中包含野生型、ATP7B基因敲除型、铜诱导的野生型和铜诱导的ATP7B基因敲除型HepG2细胞系数据。筛选由铜诱导产生的差异表达基因(differentially expressed genes,DEGs)后进行基因本体论(gene ontology,GO)、京都基因和基因组百科全书(Kyoto encyclopedia of genes and genomes, KEGG)富集分析。基于蛋白相互作用网络,识别疾病关键基因和功能模块,并对关键功能模块中的基因进行富集分析。最后,构建转录调控网络,筛选核心转录因子。共筛选出1 034个差异表达基因,其中上调525个,下调509个。上、下调关键功能模块分别包括了3785个和3931个基因。关键功能模块中的基因主要定位于细胞-基质连接、染色体、剪接复合体、核糖体等区域,共同参与了mRNA加工、组蛋白修饰、RNA剪切...     

萱草叶片响应低温胁迫的转录组分析

低温胁迫是萱草（Hemerocallis fulva）生长过程中经常会遭遇的一种非生物胁迫。比较了萱草叶片在低温处理（10、5、0 ℃）下转录组与对照（15 ℃）数据的差异,共筛选出差异表达基因2 457个,其中上调基因1 253个,下调基因1 204个。差异表达基因主要富集在细胞过程、代谢过程和催化活性等49个GO过程,代谢途径、次生代谢产物的生物合成、植物激素信号转导等42条KEGG代谢途径中。其中参与植物激素信号转导通路的差异表达基因发生了不同程度的变化,GH3.10基因上调至对照组的13.624倍,IAA1基因下调0.120倍;参与可溶性糖合成通路的差异基因发生了0.076~28.114倍不同程度的变化。随后对3个低温处理组共有的29个差异表达基因进行热图和网络调控分析,基于基因在网络调控中的位置,对ABCF5、OFPs和SWEETs等基因在冷应答的作用进行了分析。本研究结果为进一步挖掘萱草低温响应的关键基因及耐寒萱草种质开发、分子育种提供了一定的理论支撑。