Pipeline Overview

The ENCODE miRNA-seq pipeline can be used for libraries generated from miRNAs, size-selected from total RNA to be 30 bp or smaller. The microRNA-seq pipeline was developed by Ali Mortazavi's group at UC Irvine.

Pipeline Schematic

View the current instance of this pipeline.

Inputs:

File format	Information contained in file	File description	Notes
fastq	reads	microRNA-seq reads from single-end libraries.	Reads must meet the criteria outlined in the Uniform Processing Pipeline Restrictions.
gtf, tar	reference files	Genome annotations, adapter sequences, genome indices	See the References section below for details.

Outputs:

File format	Information contained in file	File description	Notes
bam	alignments	Produced by mapping reads to the genome.	Reads are trimmed using a proprietary version of cutAdapt, linked below under References.
bigWig	signal	Normalized RNA-seq signal	Signals are generated for both the plus and minus strands and for unique reads and unique+multimapping reads.
tsv	gene (microRNA) quantifications	Non-normalized counts.
Quality control metrics are also generated, including gene quantification level and read depth.

The mapping of the reads is done using the STAR aligner. STAR is also used to obtain counts of miRNAs (number of reads mapped to each miRNA gene in the annotations file). The miRNA counts can be normalized, for example by library size, to obtain counts-per-million for downstream analysis.

References

Adapter Trimming

An important step before aligning miRNA-seq reads is the trimming of adapters. The raw sequencing reads for samples generated for this pipeline contain 3’ and 5’ adapters. The links below contain details of the pipeline including information about the adapter sequences used and an example script to trim these adapters:

Genomic References

View the mapping assembly and genome annotation reference files used in this pipeline.

STAR Indices

This pipeline requires both assembly information for the species of interest and a gene reference. STAR creates an index for use in the mapping step. Please note that the comprehensive GENCODE references are used for the index generation step, while the miRNA subsets of the corresponding annotations are used at the alignment step to obtain miRNA counts by STAR.

The miRNA subset of GENCODE V29 for human (used at the mapping step by STAR)

The miRNA subset of GENCODE M21 for mouse (used at the mapping step by STAR)

Links and Publications

Find data generated by the pipeline here.

Pipeline schematic