scholarly journals Development of a bioinformatics platform for analysis of quantitative transcriptomics and proteomics data: the OMnalysis

PeerJ ◽  
2021 ◽  
Vol 9 ◽  
pp. e12415
Author(s):  
Punit Tyagi ◽  
Mangesh Bhide

Background In the past decade, RNA sequencing and mass spectrometry based quantitative approaches are being used commonly to identify the differentially expressed biomarkers in different biological conditions. Data generated from these approaches come in different sizes (e.g., count matrix, normalized list of differentially expressed biomarkers, etc.) and shapes (e.g., sequences, spectral data, etc.). The list of differentially expressed biomarkers is used for functional interpretation and retrieve biological meaning, however, it requires moderate computational skills. Thus, researchers with no programming expertise find difficulty in data interpretation. Several bioinformatics tools are available to analyze such data; however, they are less flexible for performing the multiple steps of visualization and functional interpretation. Implementation We developed an easy-to-use Shiny based web application (named as OMnalysis) that provides users with a single platform to analyze and visualize the differentially expressed data. The OMnalysis accepts the data in tabular form from edgeR, DESeq2, MaxQuant Perseus, R packages, and other similar software, which typically contains the list of differentially expressed genes or proteins, log of the fold change, log of the count per million, the P value, q-value, etc. The key features of the OMnalysis are multiple image type visualization and their dimension customization options, seven multiple hypothesis testing correction methods to get more significant gene ontology, network topology-based pathway analysis, and multiple databases support (KEGG, Reactome, PANTHER, biocarta, NCI-Nature Pathway Interaction Database PharmGKB and STRINGdb) for extensive pathway enrichment analysis. OMnalysis also fetches the literature information from PubMed to provide supportive evidence to the biomarkers identified in the analysis. In a nutshell, we present the OMnalysis as a well-organized user interface, supported by peer-reviewed R packages with updated databases for quick interpretation of the differential transcriptomics and proteomics data to biological meaning. Availability The OMnalysis codes are entirely written in R language and freely available at https://github.com/Punit201016/OMnalysis. OMnalysis can also be accessed from - http://lbmi.uvlf.sk/omnalysis.html. OMnalysis is hosted on a Shiny server at https://omnalysis.shinyapps.io/OMnalysis/. The minimum system requirements are: 4 gigabytes of RAM, i3 processor (or equivalent). It is compatible with any operating system (windows, Linux or Mac). The OMnalysis is heavily tested on Chrome web browsers; thus, Chrome is the preferred browser. OMnalysis works on Firefox and Safari.

2020 ◽  
Author(s):  
Kumari Sonal Choudhary ◽  
Eoin Fahy ◽  
Kevin Coakley ◽  
Manish Sud ◽  
Mano R Maurya ◽  
...  

ABSTRACTWith the advent of high throughput mass spectrometric methods, metabolomics has emerged as an essential area of research in biomedicine with the potential to provide deep biological insights into normal and diseased functions in physiology. However, to achieve the potential offered by metabolomics measures, there is a need for biologist-friendly integrative analysis tools that can transform data into mechanisms that relate to phenotypes. Here, we describe MetENP, an R package, and a user-friendly web application deployed at the Metabolomics Workbench site extending the metabolomics enrichment analysis to include species-specific pathway analysis, pathway enrichment scores, gene-enzyme information, and enzymatic activities of the significantly altered metabolites. MetENP provides a highly customizable workflow through various user-specified options and includes support for all metabolite species with available KEGG pathways. MetENPweb is a web application for calculating metabolite and pathway enrichment analysis.Availability and ImplementationThe MetENP package is freely available from Metabolomics Workbench GitHub: (https://github.com/metabolomicsworkbench/MetENP), the web application, is freely available at (https://www.metabolomicsworkbench.org/data/analyze.php)


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Xin Yuan ◽  
Shenqiang Hu ◽  
Liang Li ◽  
Chunchun Han ◽  
Hehe Liu ◽  
...  

Abstract Background Despite their important functions and nearly ubiquitous presence in cells, an understanding of the biology of intracellular lipid droplets (LDs) in goose follicle development remains limited. An integrated study of lipidomic and transcriptomic analyses was performed in a cellular model of stearoyl-CoA desaturase (SCD) function, to determine the effects of intracellular LDs on follicle development in geese. Results Numerous internalized LDs, which were generally spherical in shape, were dispersed throughout the cytoplasm of granulosa cells (GCs), as determined using confocal microscopy analysis, with altered SCD expression affecting LD content. GC lipidomic profiling showed that the majority of the differentially abundant lipid classes were glycerophospholipids, including PA, PC, PE, PG, PI, and PS, and glycerolipids, including DG and TG, which enriched glycerophospholipid, sphingolipid, and glycerolipid metabolisms. Furthermore, transcriptomics identified differentially expressed genes (DEGs), some of which were assigned to lipid-related Gene Ontology slim terms. More DEGs were assigned in the SCD-knockdown group than in the SCD-overexpression group. Integration of the significant differentially expressed genes and lipids based on pathway enrichment analysis identified potentially targetable pathways related to glycerolipid/glycerophospholipid metabolism. Conclusions This study demonstrated the importance of lipids in understanding follicle development, thus providing a potential foundation to decipher the underlying mechanisms of lipid-mediated follicle development.


2021 ◽  
Vol 80 (Suppl 1) ◽  
pp. 1037.2-1038
Author(s):  
X. Sun ◽  
S. X. Zhang ◽  
S. Song ◽  
T. Kong ◽  
C. Zheng ◽  
...  

Background:Psoriasis is an immune-mediated, genetic disease manifesting in the skin or joints or both, and also has a strong genetic predisposition and autoimmune pathogenic traits1. The hallmark of psoriasis is sustained inflammation that leads to uncontrolled keratinocyte proliferation and dysfunctional differentiation. And it’s also a chronic relapsing disease, which often necessitates a long-term therapy2.Objectives:To investigate the molecular mechanisms of psoriasis and find the potential gene targets for diagnosis and treating psoriasis.Methods:Total 334 gene expression data of patients with psoriasis research (GSE13355 GSE14905 and GSE30999) were obtained from the Gene Expression Omnibus database. After data preprocessing and screening of differentially expressed genes (DEGs) by R software. Online toll Metascape3 was used to analyze Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of DEGs. Interactions of proteins encoded by DEGs were discovered by Protein-protein interaction network (PPI) using STRING online software. Cytoscape software was utilized to visualize PPI and the degree of each DEGs was obtained by analyzing the topological structure of the PPI network.Results:A total of 611 DEGs were found to be differentially expressed in psoriasis. GO analysis revealed that up-regulated DEGs were mostly associated with defense and response to external stimulus while down-regulated DEGs were mostly associated with metabolism and synthesis of lipids. KEGG enrichment analysis suggested they were mainly enriched in IL-17 signaling, Toll-like receptor signaling and PPAR signaling pathways, Cytokine-cytokine receptor interaction and lipid metabolism. In addition, top 9 key genes (CXCL10, OASL, IFIT1, IFIT3, RSAD2, MX1, OAS1, IFI44 and OAS2) were identified through Cytoscape.Conclusion:DEGs of psoriasis may play an essential role in disease development and may be potential pathogeneses of psoriasis.References:[1]Boehncke WH, Schon MP. Psoriasis. Lancet 2015;386(9997):983-94. doi: 10.1016/S0140-6736(14)61909-7 [published Online First: 2015/05/31].[2]Zhang YJ, Sun YZ, Gao XH, et al. Integrated bioinformatic analysis of differentially expressed genes and signaling pathways in plaque psoriasis. Mol Med Rep 2019;20(1):225-35. doi: 10.3892/mmr.2019.10241 [published Online First: 2019/05/23].[3]Zhou Y, Zhou B, Pache L, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun 2019;10(1):1523. doi: 10.1038/s41467-019-09234-6 [published Online First: 2019/04/05].Acknowledgements:This project was supported by National Science Foundation of China (82001740), Open Fund from the Key Laboratory of Cellular Physiology (Shanxi Medical University) (KLCP2019) and Innovation Plan for Postgraduate Education in Shanxi Province (2020BY078).Disclosure of Interests:None declared


2018 ◽  
Vol 2018 ◽  
pp. 1-10 ◽  
Author(s):  
Binbin Xie ◽  
Yiran Li ◽  
Rongjie Zhao ◽  
Yuzi Xu ◽  
Yuhui Wu ◽  
...  

Chemoresistance is a significant factor associated with poor outcomes of osteosarcoma patients. The present study aims to identify Chemoresistance-regulated gene signatures and microRNAs (miRNAs) in Gene Expression Omnibus (GEO) database. The results of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) included positive regulation of transcription, DNA-templated, tryptophan metabolism, and the like. Then differentially expressed genes (DEGs) were uploaded to Search Tool for the Retrieval of Interacting Genes (STRING) to construct protein-protein interaction (PPI) networks, and 9 hub genes were screened, such as fucosyltransferase 3 (Lewis blood group) (FUT3) whose expression in chemoresistant samples was high, but with a better prognosis in osteosarcoma patients. Furthermore, the connection between DEGs and differentially expressed miRNAs (DEMs) was explored. GEO2R was utilized to screen out DEGs and DEMs. A total of 668 DEGs and 5 DEMs were extracted from GSE7437 and GSE30934 differentiating samples of poor and good chemotherapy reaction patients. The Database for Annotation, Visualization, and Integrated Discovery (DAVID) was used to perform GO and KEGG pathway enrichment analysis to identify potential pathways and functional annotations linked with osteosarcoma chemoresistance. The present study may provide a deeper understanding about regulatory genes of osteosarcoma chemoresistance and identify potential therapeutic targets for osteosarcoma.


2020 ◽  
Author(s):  
Lun Wu ◽  
Ying Wei ◽  
Wen-Bo Zhou ◽  
Jiao Zhou ◽  
Li-Hua Yang ◽  
...  

Abstract Background Borax, a boron compound, which is becoming widely recognized for its biological effects, including antioxidant activity, cytotoxicity, and potential therapeutic benefits. However, the specific molecular mechanisms underlying borax-induced anti-tumor effect still remain to be to further elucidated. MicroRNAs (miRNAs) may play key roles in cellular processes including tumor progression, cell apoptosis and cytotoxicity. Thus, this study aimed to investigate, whether miRNAs were involved in the borax-mediated anti-tumor effect using miRNA profiling of a human liver cancer cell line (HepG2) using gene-chip analysis.Methods Total RNA was extracted and purified from HepG2 cells that were treated with 4 mM borax for either 2 or 24 h. The samples underwent microarray analysis using an Agilent Human miRNA Array. Differentially expressed miRNAs were analysed by volcano plot and heatmap, and were validated using real-time fluorescent quantitative PCR (qPCR).ResultsAmong this, 2- or 24-h exposure to borax significantly altered the expression level of miRNAs in HepG2 cells, 4 or 14 were upregulated and 3 were downregulated compared with the control group, respectively (≥2-fold; P<0.05). GO enrichment analysis and KEGG pathway enrichment analysis revealed that target genes of differentially expressed miRNAs in HepG2 cells predominantly participated in MAPK signaling pathway, TGF-beta signaling pathway, NF-kappa B signaling pathway, etc; in 2-h borax treatment group, while Ras signaling pathway, FoxO signaling pathway, Cellular senescence, etc; involved in 24-h treatment group.Conclusions Result indicates that borax-induced anti-tumor effect may be associated with alterations in miRNAs.


2021 ◽  
Author(s):  
Li Guoquan ◽  
Du Junwei ◽  
He Qi ◽  
Fu Xinghao ◽  
Ji Feihong ◽  
...  

Abstract BackgroundHashimoto's thyroiditis (HT), also known as chronic lymphocytic thyroiditis, is a common autoimmune disease, which mainly occurs in women. The early manifestation was hyperthyroidism, however, hypothyroidism may occur if HT was not controlled for a long time. Numerous studies have shown that multiple factors, including genetic, environmental, and autoimmune factors, were involved in the pathogenesis of the disease, but the exact mechanisms were not yet clear. The aim of this study was to identify differentially expressed genes (DEGs) by comprehensive analysis and to provide specific insights into HT. MethodsTwo gene expression profiles (GSE6339, GSE138198) about HT were downloaded from the Gene Expression Omnibus (GEO) database. The DEGs were assessed between the HT and normal groups using the GEO2R. The DEGs were then sent to the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. The hub genes were discovered using Cytoscape and CytoHubba. Finally, NetworkAnalyst was utilized to create the hub genes' targeted microRNAs (miRNAs). ResultsA total of 62 DEGs were discovered, including 60 up-regulated and 2 down-regulated DEGs. The signaling pathways were mainly engaged in cytokine interaction and cytotoxicity, and the DEGs were mostly enriched in immunological and inflammatory responses. IL2RA, CXCL9, IL10RA, CCL3, CCL4, CCL2, STAT1, CD4, CSF1R, and ITGAX were chosen as hub genes based on the results of the protein-protein interaction (PPI) network and CytoHubba. Five miRNAs, including mir-24-3p, mir-223-3p, mir-155-5p, mir-34a-5p, mir-26b-5p, and mir-6499-3p, were suggested as likely important miRNAs in HT. ConclusionsThese hub genes, pathways and miRNAs contribute to a better understanding of the pathophysiology of HT and offer potential treatment options for HT.


PeerJ ◽  
2020 ◽  
Vol 8 ◽  
pp. e8831 ◽  
Author(s):  
Xiaojiao Guan ◽  
Yao Yao ◽  
Guangyao Bao ◽  
Yue Wang ◽  
Aimeng Zhang ◽  
...  

Esophageal cancer is a common malignant tumor in the world, and the aim of this study was to screen key genes related to the development of esophageal cancer using a variety of bioinformatics analysis tools and analyze their biological functions. The data of esophageal squamous cell carcinoma from the Gene Expression Omnibus (GEO) were selected as the research object, processed and analyzed to screen differentially expressed microRNAs (miRNAs) and differential methylation genes. The competing endogenous RNAs (ceRNAs) interaction network of differentially expressed genes was constructed by bioinformatics tools DAVID, String, and Cytoscape. Biofunctional enrichment analysis was performed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). The expression of the screened genes and the survival of the patients were verified. By analyzing GSE59973 and GSE114110, we found three down-regulated and nine up-regulated miRNAs. The gene expression matrix of GSE120356 was calculated by Pearson correlation coefficient, and the 11696 pairs of ceRNA relation were determined. In the ceRNA network, 643 lncRNAs and 147 mRNAs showed methylation difference. Functional enrichment analysis showed that these differentially expressed genes were mainly concentrated in the FoxO signaling pathway and were involved in the corresponding cascade of calcineurin. By analyzing the clinical data in The Cancer Genome Atlas (TCGA) database, it was found that four lncRNAs had an important impact on the survival and prognosis of esophageal carcinoma patients. QRT-PCR was also conducted to identify the expression of the key lncRNAs (RNF217-AS1, HCP5, ZFPM2-AS1 and HCG22) in ESCC samples. The selected key genes can provide theoretical guidance for further research on the molecular mechanism of esophageal carcinoma and the screening of molecular markers.


2020 ◽  
Author(s):  
Bolin Wu ◽  
Haitao Shang ◽  
Xitian Liang ◽  
Huajing Yang Huajing Yang ◽  
Hui Jing ◽  
...  

Abstract Background: Hepatocellular carcinoma (HCC) poses a severe threat to human health. The NET-1 protein has been proved to be strongly associated with HCC proliferation and metastasis in our previous study. Methods: Here, we developed a label-free proteome mass spectrometry workflow to analyze formalin-fixed and paraffin-embedded HCC xenograft samples collected in our previous study. Results: The result showed that 78 proteins were differentially expressed after NET-1 protein inhibited. Among them, the expression of 61 proteins up-regulated and the expression of 17 proteins were significantly down-regulated. Of the differentially expressed proteins, the vast majority of Gene Ontology enrichment terms belong to the biological process. The KEGG pathway enrichment analysis showed that the 78 differentially expressed proteins significantly enriched in 45 pathways. We concluded that the function of the NET-1 gene is not only to regulate HCC but also to participate in a variety of biochemical metabolic pathways in the human body. Furthermore, the protein-protein interaction analysis indicated that the interactions of differentially expressed proteins are incredibly sophisticated. All the protein-protein interactions happened after the NET-1 gene has been silenced. Conclusions: Finally, our study also provides a useful proposal for targeted therapy based on tetraspanin proteins to treat HCC, and further mechanism investigations are needed to reveal a more detailed mechanism of action for NET-1 protein regulation of HCC.


2020 ◽  
Vol 10 (1) ◽  
Author(s):  
Guoqing Shen ◽  
Xiao Zhang ◽  
Jie Gong ◽  
Yang Wang ◽  
Pengdan Huang ◽  
...  

AbstractEach year from April to May, high mortality rates are reported in red swamp crayfish (Procambarus clarkii) cultured in Jiangsu and other regions, in China, and this phenomenon has come to be known as “Black May” disease (BMD). Therefore, in order to investigate the possible causes of this disease, this study gathered BMD-affected P. clarkii samples and performed transcriptome analysis on hepatopancreas, gill, and muscle tissues. A total of 19,995,164, 149,212,804, and 222,053,848 clean reads were respectively obtained from the gills, muscle, and hepatopancreas of BMD-affected P. clarkii, and 114,024 unigenes were identified. The number of differentially expressed genes (DEGs) in gill, muscle, and hepatopancreas was 1703, 964, and 476, respectively. GO and KEGG enrichment analyses of the DEGs were then conducted. Based on KEGG pathway enrichment analysis, the most significantly differentially expressed pathways were mainly those involved with metabolism, human disease, and cellular processes. Further analysis of the significantly DEGs revealed that they were mainly related to the mitochondrial-mediated apoptosis pathway and that the expression of these DEGs was mostly down-regulated. Moreover, the expression of genes related to immune and metabolism-related pathways was also significantly down-regulated, and these significantly-inhibited pathways were the likely causes of P. clarkii death. Therefore, our results provide a basis for the identification of BMD causes.


2019 ◽  
Vol 244 (18) ◽  
pp. 1648-1657
Author(s):  
Yuan Li ◽  
Dan Yang ◽  
Bo Sun ◽  
Xu Zhang ◽  
Fangda Li ◽  
...  

As a common disease, abdominal aortic aneurysm (AAA) features permanently progressively dilated abdominal aorta. Various cytokines are implicated in AAA pathogenesis. Clarification of involved cytokines combined with functional analysis may provide new insights into AAA pathogenesis. Using a mouse model, this study analyzed the cytokine profiles in AAA. Cytokines were measured in AAA tissues of saline control or angiotensin II-treated ApoE−/− mice using an antibody array of 200 cytokines, cytokine receptors, and related proteins. Statistical analysis revealed that 21 of 200 proteins were differentially expressed in AAA. These differentially expressed proteins were subjected to function and pathway enrichment analysis, which revealed that leukocyte migration and positive regulation of cell adhesion were the most significant biological processes. Specific signaling pathways, including Janus kinase/signal transducers and activators of transcription and cytokine–cytokine receptor interaction, were prominent in Kyoto encyclopedia of genes and genomes pathway enrichment analysis. Importantly, our data identified cytokines which had not previously been illustrated in AAA pathogenic pathways. Bivariate correlation analysis between these cytokines and protease activity showed that granulocyte colony-stimulating factor (G-CSF), macrophage inflammatory protein 1 g, cardiotrophin 1, milk fat globule-EGF factor 8 protein, interleukin 33, and periostin were positively correlated with matrix metalloprotease 1 (MMP-1), MMP-9, cathepsin B, and cathepsin L. G-CSF was positively correlated with cathepsin L. In conclusion, these results demonstrate that cytokine profile is significantly altered in AAA, and that the newly identified crucial cytokines may function potentially in AAA pathogenesis. Impact statement Various cytokines are known contributors to abdominal aortic aneurysm (AAA) pathologic processes, but the mechanisms underlying the pathogenesis remains unclear. We illustrated the altered cytokine profiles in AAA by high throughput antibody array of 200 cytokines, cytokine receptors and related proteins, as well as bioinformatics analysis of differentially expressed proteins in lesion tissues from AAA mice infused with angiotensin II. Functional analyses of differentially expressed cytokines showed clustering on cell migration and adhesion processes. More importantly, crucial cytokines whose association with AAA formation had not been established were identified. Significant correlations were found between these cytokines and protease activity. This study identifies several crucial markers for further researches on the molecular basis of AAA.


Sign in / Sign up

Export Citation Format

Share Document