Running RNA-seq analyses pipelines in KBase
KBase offers a powerful suite of expression analysis tools. Starting with short reads, you can use the tool suite to assemble, quantify long transcripts, and identify differentially expressed genes. You can also compare the expression data with the flux when studying metabolic models in KBase and identify pathways where expression and flux agree or conflict.
KBase requires a reference genome to guide the analysis of short reads.
- 2.Import Short Reads: The reads must be a set of single-end, paired-end, or interleaved paired-end reads in FASTA, FASTQ, or SRA format.
- 3.Create a SampleSet: Run the Create RNA-seq Sample Set App to group together your reads into an RNA-seq sample set with associated experimental metadata to run RNA-seq Apps in batch mode wherever appropriate.
- 4.QC SampleSet: Run FastQC to assess the read quality of the reads set from the previous step and if needed, run Trimmomatic, Cutadapt, or PRINSEQ to pre-process or filter the reads before starting RNA-seq analysis.
The RNA-seq pipeline in KBase is modular and consists of three steps. You can pick any of the multiple Apps available for a given step depending on your preference or individual characteristics of the App.
- 1.Read Alignment: Align reads to map short reads to the reference genome. The output is a set of BAM alignments and Qualimap report. You can download the alignment output object generated by aligner Apps for further analysis.
- 2.Transcriptome Assembly and Quantification: Assemble aligned reads to generate full-length transcripts and quantify transcripts and genes as appropriate. You can view downloadable normalized full expression matrices in FPKM (fragments per kilobase of exon model per million mapped reads) and TPM (transcripts per million).
- 3.Differential Gene Expression: Generate gene- or transcript-level differential expression based on the quantification. Run Create Feature Set/Filtered Expression Matrix From Differential Expression after selecting appropriate q-value and fold change cutoffs as input parameters for the filtering of the differential gene expression.
- 1.Filtering: You can create a filtered expression matrix and associated feature set based on fold-change or adjusted p-value. You can also filter an expression matrix based on LOR or ANOVA.
- 2.Clustering: Depending on preference, run the Hierarchical, K-Means or WGCNA clustering App to group features into clusters based on gene expression. You can also visualize the clusters as an interactive heatmap.
- 3.Functional Enrichment: Assess the functional enrichment in plant genomes for a set of features using associated GO terms.
- 4.Integration into Metabolic Models: Assimilate the expression data from RNA-seq into the metabolic models to compare reaction fluxes with gene expression and thus identify pathways where expression and flux agree or conflict.
- E. coli RNA-seq Analysis Tutorial – bacterium-based example of an RNAseq workflow using a HISAT2/StringTie/DESeq2 pipeline
- Arabidopsis RNA-seq Analysis Tutorial – plant-based example of an RNAseq workflow using a HISAT2/StringTie/DESeq2 pipeline
- Case Study: Genome-wide Transcriptomics and Plant Primary Metabolism in response to Drought Stress in Sorghum – plant-based example of using KBase to integrate full genome RNA-seq analysis and a metabolic model to generate a reaction matrix
- Kumari et al. (2021) A KBase case study on genome-wide transcriptomics and plant primary metabolism in response to drought stress in Sorghum. Current Plant Biology 28. https://doi.org/10.1016/j.cpb.2021.100229