How can I analyze my miRVEL Small RNA sequencing data ?
miRVEL small RNA data can be analyzed as a Service project with Lexogen! For more information, please visit our Bioinformatic Services webpage and contact us at services@lexogen.com.
Alternatively, you can analyze the data on your own. As reference, you can use the following publicly available data analysis workflows:
Small RNA analysis pipeline built using modified Nextflow workflow tool: https://nf-co.re/smrnaseq/1.0.0.
Workflow adapted from Tam et al., (2015) Optimization of miRNA-seq data preprocessing. Briefings in Bioinformatics 16(6). This paper also provides an additional comparison of programs for read mapping (alignment), normalization methods, and differential expression analysis.
The general workflow involves trimming of adapters, size filtering of trimmed reads, read mapping:
Trimming, adapter-removal, and size filtering (BBDuk, https://archive.jgi.doe.gov/data-and-tools/software-tools/bbtools/bb-tools-user-guide/bbduk-guide/)
1.1. Use bbduk.sh to remove the adapter sequence from the 3' end of the reads and trim homopolymers.
1.2. Split up the reads into short read and long read fractions.
1.3. Short reads read into the adapter sequence AND are shorter than 31 bp and are processed for miRNA alignment.
1.4. Long reads are kept separate for alignment to the reference genome.
NOTE: Please see the following FAQs for the adapter sequences trimming for miRVEL Discovery and miRVEL Profiling data.
Mapping* (BWA, http://bio-bwa.sourceforge.net/, BBMap, https://sourceforge.net/projects/bbmap/)
2.1. Short read fraction
2.2. Filter out reads <= 14 bp -> outputs: "filtered reads: too short" statistic
2.3. Map the remaining 15 – 31 bp reads to the ribosomal RNA sequences using BWA. -> outputs: "mapped reads" and "unmapped reads", mapping statistics
2.4. Use the unmapped reads from the BWA alignment and map these against the mirBase annotation (http://www.mirbase.org; mature miRNAs, species-specific) -> output "mapped reads", read count statistics, and miRNA-count table.
2.5. Long read fraction
2.6. Map long reads against the reference genome using STAR (https://github.com/alexdobin/STAR) -> output: "mapped reads", read count statistics, gene_biotype statistics (including rRNA sequences), and gene-read table.