This is a work-in-progress repository, automating the khmer protocols over a large-scale RNAseq data set:
https://khmer-protocols.readthedocs.org/en/ctb/mrnaseq/
The Marine Microbial Eukaryotic Transcriptome Sequencing Project (MMETSP) data set contains 678 cultured samples of 306 pelagic and endosymbiotic marine eukaryotic species representing more than 40 phyla (Keeling et al. 2014).
Automated scripts running the protocols:
- getdata.py, download data from NCBI and organize into individual directories for each sample/accession ID
- trim_qc.py, trim reads for quality, interleave reads
- diginorm_mmetsp.py, normalize-by-median and filter-abund from khmer, rename, combined orphans,
- assembly.py, runs Trinity de novo transcriptome assembly software
- salmon.py, runs salmon reference-free transcript quantification
References:
Keeling et al. 2014: http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.1001889
Supporting information with methods description:
http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.1001889#s6
Preliminary assembly protocol run by NCGR: https://github.com/ncgr/rbpa
NCBI Bioproject accession: PRJNA231566
http://www.ncbi.nlm.nih.gov/bioproject/PRJNA231566/
MMETSP website: http://marinemicroeukaryotes.org/
iMicrobe project with data and combined assembly downloads: http://data.imicrobe.us/project/view/104