BxGenomics enables biologists to easily analyze single-cell and bulk RNA-Seq data, identify changed genes and enriched pathways, and further visualize the results and compare across projects. The integrated (sc)RNA-Seq data solution of BxGenomics includes three key parts.
- Efficient data management system for public and private datasets
- Powerful analytical pipelines
- Interactive data visualization tools
Explore Public Single-Cell and Bulk RNA-Seq Data
We developed an integrated web database including single-cell RNA-Seq data from over 400 published datasets from multiple species and many tissues. It provides an efficient data management system for single cell datasets. Users can search the database by multiple attributes such as publication, key words, species and cell types. Advanced filter and across-dataset gene expression tools enable users to quickly find the datasets of interest and compare the expression of gene(s) across these.
Importantly, each scRNA dataset can be loaded into the cellxgene VIP tool, which enables fast exploration of individual dataset in the manner of interactivity and scalability to gain more refined insights such as cell composition, gene expression profiles, and differentially expressed genes among cell types by leveraging more than 20 frequently applied plotting functions and high-level analysis methods in single cell research.
Figure legend: Cellxgene VIP includes analytical modules that provide essential functions for interactive visualization and generation of publication-ready plots. (a) Multi-tSNE/UMAP plot visually highlights which cells expressing cell markers on selected embedding (UMAP based on harmony batch correction in this example). (b) Dual-gene plot highlights cells express SYT1 and GAD1 (green SYT1 only, red GAD1 only, yellow co-expression of STY1 and GAD1), expression cutoff 2.2. (c) Stacked barplot demonstrates the fraction of each major cell type across each sample (C are Control and MS are MS patients). (d) Trackplot shows expression of lineage marker genes across individual cells in annotated clusters. (e) Violin plot shows the AQP4 gene expression across cell types. (f) Sankey diagram (a.k.a. Riverplot) provides quick and easy way to explore the inter-dependent relationship of variables in the MS snRNAseq dataset8. (g) Density plots shows expression of marker genes across annotated clusters and split across cell types. (h) Stacked violin and Dot plot are the key visualizations of selected cell markers across cell types. They highlight their selective expression and validates the scRNAseq approach and visualization method. (i) Command Line Interface (CLI) exposed by mini Jupyter Notebook to provide maximal flexibility of doing various analytics on the whole or sliced single cell dataset.
About 700,000 human samples from over 12,000 projects in GEO and SRA are included in the BxGenomics Bulk RNA-Seq View platform. For each project, the sample meta data is extracted and can be used to annotate samples and specify conditions for comparisons. For each sample, gene count data is used to perform statistical analysis, and gene TMP data is used for visualization of gene expression.
Explore Public RNA-Seq DataBxGenomics system includes data processing capabilities to enable users to perform differentially expression gene (DEG) analysis and pass the results to the downstream data mining application QuickOmics.
Quickomics is a feature-rich R Shiny-powered tool that enable biologists to fully explore RNA-Seq results and perform advanced analysis in an easy-to-use interface. It covers a broad range of secondary and tertiary analytical tasks. Each functional module is equipped with customizable options and generates both interactive and publication-ready plots to uncover biological insights from data.
Figure legend: Selected Quickomics functions applied to a dataset of microglial RNA-seq gene expression from three mouse genotypes over time. (A) PCA based on full dataset highlights primary sample separation by mouse age at which the cells were isolated. (B) Volcano plot visualizes differentially expressed genes, most of which show reduced expression in 2mo KO compared to 2mo_WT microglia. For spacing purpose, absolute log2FC (Fold Change) and negative log10 adjusted p-value are capped at 1.5 and 15, respectively. (C) Correlation analysis between two comparisons shows that aging and Cx3cr1-KO have a similar effect on gene expression. (D) Pattern clustering identifies subsets of genes with similar expression over the samples. The clustering is mostly driven by age, with the KO genotype having a similar, but smaller effect. (E) Heatmap of all samples allows the identification of gene clusters with expression regulated by age and/or genotype. Key genes and the pathways they belong to are highlighted on the right. (F) After pathway enrichment analysis, KEGG pathways (Kanehisa and Goto, 2000) of interest can be displayed in a cellular context. The color bars with each stripe representing one comparison show log2 fold changes in various comparisons, allowing project-wide insights for patterns of expression. (G) Correlation network shows potential links between genes of interest.
Analyze Your Own Data
BxGenomics enables biologists to easily analyze bulk and single-cell RNA-Seq data, identify changed genes and enriched pathways, and further visualize the results and compare across projects using interactive data mining tools.
Why Use BxGenomics and Related Services?
Biologically Meaningful Answers
Immediately find changed genes and functional pathways in instant reports.
Full Potential of Data
Explore your results in an interactive online database to delve into details.
Expert Bioinformatics Support
Rest assured that your data are in good hands and help is always available.
Data Analysis Service Overview
| Type of Analysis | Bulk RNA-Seq | scRNA-Seq |
|---|---|---|
| Input | Raw fastq files or gene count and gene TPM files | Count matrix file, or raw fastq files |
| Gene Level Analysis | QC (PCA, covariates); gene expression plot, clustering and correlation; heatmap with functional gene sets. | Multi-UMAP/tSNE; violin and dot plot for genes of interest; stacked barplot for cell types; Sankey diagram for variable relationship. |
| Advanced Analysis | Differential gene expressed between conditions; Venn diagram of comparisons; functional enrichment and pathway overlay. | Identify marker genes; assign cell types based on reference data; differential gene expression and gene set enrichment. |
| Visualization Tool | QuickOmics | Cellxgene VIP |
RNA-Seq Analysis service is easy to use. Just provide raw data and sample description, and a professional report and personal online database will be delivered to you in 1-2 weeks.
An expert-designed analysis pipeline will work for you behind the scenes, including comprehensive QC of raw data and gene counts, robust statistical analysis for differentially expressed genes, advanced functional pathway analysis, and more. All results are reviewed by experts before final delivery.
Once the RNA-Seq data is upload to the BxGenomics data mining platform, authorized users can access the data with a browser anytime from anywhere with internet connection.
Different kinds of gene IDs are automatically recognized and converted so gene expression data are easily integrated between different projects analyzed with different types of gene IDs.
The BxGenomics platform serves dual purposes, as a data mining system of all RNA-Seq results in the laboratory to enable continuous discovery, and as a data archive system to enable data longevity and long-term access.
The BxGenomics RNA-Seq service team consists of highly educated scientists and customer support managers with years of experience in biological studies and genomics.
The team works collaboratively with customers to provide the utmost service for the success of every RNA-Seq project. Customers benefit from continuous high-quality technical support by a friendly and responsive team.
Methods and References
Methods
Bulk RNA-Seq Analysis
- Generate gene count, one of the methods will be used based on the project.
- Subread + Featurecount
- STAR + RSEM
- Kallisto or Salmon
- Differential Expression, limma or DESeq2
- Data Visualization, QuickOmics
Bulk RNA-Seq Analysis
- Count Matrix from 10X raw data: Cellranger
- Data analysis, Seurat, Scanpy, Harmony, LIGER
- Data Visualization, Cellxgene VIP
- Single-Cell database, CellDepot
References
- Liao Y, Duan B, Zhang Y, Zhang X, Xia B. Excessive ER-phagy mediated by the autophagy receptor FAM134B results in ER stress, the unfolded protein response, and cell death in HeLa cells. J Biol Chem. 2019 Dec 27;294(52):20009-20023. doi: 10.1074/jbc.RA119.008709. Epub 2019 Nov 20. PMID: 31748416, PMCID: PMC6937584
- Kong G, You X, Wen Z, Chang YI, Qian S, Ranheim EA, Letson C, Zhang X, Zhou Y, Liu Y, Rajagopalan A, Zhang J, Stieglitz E, Loh M, Hofmann I, Yang D, Zhong X, Padron E, Zhou L, Pear WS, Zhang J. Downregulating Notch counteracts Kras(G12D)-induced ERK activation and oxidative phosphorylation in myeloproliferative neoplasm. Leukemia. 2019 Mar;33(3):671-685. doi: 10.1038/s41375-018-0248-0. Epub 2018 Sep 11. PMID: 30206308, PMCID: PMC6405304
- Zhang J, Kong G, Rajagopalan A, Lu L, Song J, Hussaini M, Zhang X, Ranheim EA, Liu Y, Wang J, Gao X, Chang YI, Johnson KD, Zhou Y, Yang D, Bhatnagar B, Lucas DM, Bresnick EH, Zhong X, Padron E, Zhang J. p53-/- synergizes with enhanced NrasG12D signaling to transform megakaryocyte-erythroid progenitors in acute myeloid leukemia. Blood. 2017 Jan 19;129(3):358-370. doi: 10.1182/blood-2016-06-719237. Epub 2016 Nov 4. PMID: 27815262, PMCID: PMC5248933
- Suter B, Zhang X, Pesce CG, Mendelsohn AR, Dinesh-Kumar SP, Mao JH. Next-Generation Sequencing for Binary Protein-Protein Interactions. Front Genet. 2015 Dec 17;6:346. doi: 10.3389/fgene.2015.00346. eCollection 2015. PubMed [citation] PMID: 26734059, PMCID: PMC4681833
- Combined MEK and JAK inhibition abrogates murine myeloproliferative neoplasm. Kong G, Wunderlich M, Yang D, Ranheim EA, Young KH, Wang J, Chang YI, Du J, Liu Y, Tey SR, Zhang X, Juckett M, Mattison R, Damnernsawad A, Zhang J, Mulloy JC, Zhang J. The Journal of Clinical Investigation. 2014 May 8; 124(6): 2762-2773 PMC [article] PMCID: PMC4038579, PMID: 24812670, DOI: 10.1172/JCI74182
- Still AJ, Floyd BJ, Hebert AS, Bingman CA, Carson JJ, Gunderson DR, Dolan BK, Grimsrud PA, Dittenhafer-Reed KE, Stapleton DS, Keller MP, Westphall MS, Denu JM, Attie AD, Coon JJ, Pagliarini DJ. Quantification of mitochondrial acetylation dynamics highlights prominent sites of metabolic regulation. J Biol Chem. 2013 Sep 6;288(36):26209-19. doi: 10.1074/jbc.M113.483396. Epub 2013 Jul 17. PubMed [citation] PMID: 23864654, PMCID: PMC3764825
- Wang J, Kong G, Liu Y, Du J, Chang YI, Tey SR, Zhang X, Ranheim EA, Saba-El-Leil MK, Meloche S, Damnernsawad A, Zhang J, Zhang J. Nras(G12D/+) promotes leukemogenesis by aberrantly regulating hematopoietic stem cell functions. Blood. 2013 Jun 27;121(26):5203-7. doi: 10.1182/blood-2012-12-475863. Epub 2013 May 17. PMID: 23687087, PMCID: PMC3695364
- Grimsrud PA, Carson JJ, Hebert AS, Hubler SL, Niemi NM, Bailey DJ, Jochem A, Stapleton DS, Keller MP, Westphall MS, Yandell BS, Attie AD, Coon JJ, Pagliarini DJ. A quantitative map of the liver mitochondrial phosphoproteome reveals posttranslational control of ketogenesis. Cell Metab. 2012 Nov 7;16(5):672-83. doi: 10.1016/j.cmet.2012.10.004. PMID: 23140645, PMCID: PMC3506251
- Mougeot JL, Li Z, Price AE, Wright FA, Brooks BR. Microarray analysis of peripheral blood lymphocytes from ALS patients and the SAFE detection of the KEGG ALS pathway. BMC Med Genomics. 2011 Oct 25;4:74. doi: 10.1186/1755-8794-4-74. PMID: 22027401, PMCID: PMC3219589
- Galliher-Beckley AJ, Williams JG, Cidlowski JA. Ligand-independent phosphorylation of the glucocorticoid receptor integrates cellular stress pathways with nuclear receptor signaling. Mol Cell Biol. 2011 Dec;31(23):4663-75. doi: 10.1128/MCB.05866-11. Epub 2011 Sep 19. PMID: 21930780, PMCID: PMC3232926
- Al-Dhaheri M, Wu J, Skliris GP, Li J, Higashimato K, Wang Y, White KP, Lambert P, Zhu Y, Murphy L, Xu W. CARM1 is an important determinant of ERα-dependent breast cancer cell differentiation and proliferation in breast cancer cells. Cancer Res. 2011 Mar 15;71(6):2118-28. doi: 10.1158/0008-5472.CAN-10-2426. Epub 2011 Jan 31. PMID: 21282336, PMCID: PMC3076802
- Grimsrud PA, den Os D, Wenger CD, Swaney DL, Schwartz D, Sussman MR, Ané JM, Coon JJ. Large-scale phosphoprotein analysis in Medicago truncatula roots provides insight into in vivo kinase activity in legumes. Plant Physiol. 2010 Jan;152(1):19-28. doi: 10.1104/pp.109.149625. Epub 2009 Nov 18. PMID: 19923235, PMCID: PMC2799343
- Zhu Y, Davis S, Stephens R, Meltzer PS, Chen Y. GEOmetadb: powerful alternative search engine for the Gene Expression Omnibus. Bioinformatics. 2008 Dec 1;24(23):2798-800. doi: 10.1093/bioinformatics/btn520. Epub 2008 Oct 7. PubMed [citation] PMID: 18842599, PMCID: PMC2639278
- Keller MP, Choi Y, Wang P, Davis DB, Rabaglia ME, Oler AT, Stapleton DS, Argmann C, Schueler KL, Edwards S, Steinberg HA, Chaibub Neto E, Kleinhanz R, Turner S, Hellerstein MK, Schadt EE, Yandell BS, Kendziorski C, Attie AD. A gene expression network model of type 2 diabetes links cell cycle regulation in islets with diabetes susceptibility. Genome Res. 2008 May;18(5):706-16. doi: 10.1101/gr.074914.107. Epub 2008 Mar 17. PMID: 18347327, PMCID: PMC2336811
- Zhu Y, Zhu Y, Xu W. EzArray: a web-based highly automated Affymetrix expression array data management and analysis system. BMC Bioinformatics. 2008 Jan 24;9:46. doi: 10.1186/1471-2105-9-46. PMID: 18218103, PMCID: PMC2265266