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:

1. getdata.py, download data from NCBI and organize into individual directories for each sample/accession ID
2. trim_qc.py, trim reads for quality, interleave reads
3. diginorm_mmetsp.py, normalize-by-median and filter-abund from khmer, rename, combined orphans,
4. assembly.py, runs Trinity de novo transcriptome assembly software
5. 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