drwxrwxr-x. 3 bdash bdash 149 Mar 24 17:58 Sp_ds
-rwxrwx—. 1 bdash bdash 175846179 Mar 24 13:01 Sp.ds.1M.left.fq
-rwxrwx—. 1 bdash bdash 175846179 Mar 24 13:01 Sp.ds.1M.right.fq
drwxrwxr-x. 3 bdash bdash 149 Mar 24 18:11 Sp_hs
-rwxrwx—. 1 bdash bdash 175736042 Mar 24 13:02 Sp.hs.1M.left.fq
-rwxrwx—. 1 bdash bdash 175736042 Mar 24 13:02 Sp.hs.1M.right.fq
drwxrwxr-x. 3 bdash bdash 149 Mar 24 18:28 Sp_log
-rwxrwx—. 1 bdash bdash 175741215 Mar 24 13:03 Sp.log.1M.left.fq
-rwxrwx—. 1 bdash bdash 175741215 Mar 24 13:03 Sp.log.1M.right.fq
drwxrwxr-x. 3 bdash bdash 149 Mar 24 18:20 Sp_plat
-rwxrwx—. 1 bdash bdash 175899533 Mar 24 13:03 Sp.plat.1M.left.fq
-rwxrwx—. 1 bdash bdash 175899533 Mar 24 13:03 Sp.plat.1M.right.fq
drwxr-xr-x. 2 bdash bdash 4096 Mar 25 11:27 trinity_out_dir
Combine the RSEM-computed isoform fragment counts into a matrix file for the isoforms
-reading file: Sp_ds/RSEM.isoforms.results
-reading file: Sp_hs/RSEM.isoforms.results
-reading file: Sp_plat/RSEM.isoforms.results
-reading file: Sp_log/RSEM.isoforms.results
/home/bdash/miniconda2/opt/trinity-2.5.1/util/support_scripts/run_TMM_scale_matrix.pl –matrix isoforms.isoform.TPM.not_cross_norm > isoforms.isoform.TMM.EXPR.matrixCMD: R –vanilla -q < __tmp_runTMM.R 1>&2
library(edgeR)
Loading required package: limma
rnaseqMatrix = read.table(“isoforms.isoform.TPM.not_cross_norm”, header=T, row.names=1, com=’’, check.names=F)
rnaseqMatrix = as.matrix(rnaseqMatrix)
rnaseqMatrix = round(rnaseqMatrix)
exp_study = DGEList(counts=rnaseqMatrix, group=factor(colnames(rnaseqMatrix)))
exp_study = calcNormFactors(exp_study)
exp_study\(samples\)eff.lib.size = exp_study\(samples\)lib.size * exp_study\(samples\)norm.factors
write.table(exp_study$samples, file=“isoforms.isoform.TPM.not_cross_norm.TMM_info.txt”, quote=F, sep=", row.names=F)
/home/bdash/miniconda2/opt/trinity-2.5.1/util/support_scripts/run_TMM_scale_matrix.pl –matrix isoforms.gene.TPM.not_cross_norm > isoforms.gene.TMM.EXPR.matrixCMD: R –vanilla -q < __tmp_runTMM.R 1>&2
library(edgeR)
Loading required package: limma
rnaseqMatrix = read.table(“isoforms.gene.TPM.not_cross_norm”, header=T, row.names=1, com=’’, check.names=F)
rnaseqMatrix = as.matrix(rnaseqMatrix)
rnaseqMatrix = round(rnaseqMatrix)
exp_study = DGEList(counts=rnaseqMatrix, group=factor(colnames(rnaseqMatrix)))
exp_study = calcNormFactors(exp_study)
exp_study\(samples\)eff.lib.size = exp_study\(samples\)lib.size * exp_study\(samples\)norm.factors
write.table(exp_study$samples, file=“isoforms.gene.TPM.not_cross_norm.TMM_info.txt”, quote=F, sep=", row.names=F)
-rw-rw-r–. 1 bdash bdash 384230 Mar 25 11:28 isoforms.gene.counts.matrix
-rw-rw-r–. 1 bdash bdash 419995 Mar 25 11:28 isoforms.gene.TMM.EXPR.matrix
-rw-rw-r–. 1 bdash bdash 361306 Mar 25 11:28 isoforms.gene.TPM.not_cross_norm
-rw-rw-r–. 1 bdash bdash 234 Mar 25 11:28 isoforms.gene.TPM.not_cross_norm.TMM_info.txt
-rw-rw-r–. 1 bdash bdash 431596 Mar 25 11:28 isoforms.isoform.counts.matrix
-rw-rw-r–. 1 bdash bdash 475183 Mar 25 11:28 isoforms.isoform.TMM.EXPR.matrix
-rw-rw-r–. 1 bdash bdash 411944 Mar 25 11:28 isoforms.isoform.TPM.not_cross_norm
-rw-rw-r–. 1 bdash bdash 234 Mar 25 11:28 isoforms.isoform.TPM.not_cross_norm.TMM_info.txt
-rw-rw-r–. 1 bdash bdash 526 Mar 25 11:28 __tmp_runTMM.R
-rw-rw-r–. 1 bdash bdash 142305396 Mar 24 17:55 bowtie.bam
-rw-rw-r–. 1 bdash bdash 0 Mar 24 17:55 bowtie.bam.ok
Sp_ds Sp_hs Sp_plat Sp_log
TRINITY_DN639_c0_g1_i1 1.00 1.00 6.00 2.00
TRINITY_DN3154_c0_g1_i2 103.26 101.58 108.72 43.73
TRINITY_DN1606_c0_g1_i1 0.00 0.00 0.00 0.00
TRINITY_DN4534_c0_g1_i1 0.00 0.00 0.88 1.58
TRINITY_DN694_c0_g1_i1 317.00 235.00 377.00 506.00
TRINITY_DN3962_c0_g2_i1 0.00 0.00 0.00 1.00
TRINITY_DN3383_c0_g2_i1 2.00 0.00 3.00 10.00
TRINITY_DN2506_c0_g1_i2 0.00 0.00 1.82 4.60
TRINITY_DN1195_c0_g1_i1 1.00 0.00 0.00 4.00
trinity2$ run_DE_analysis.pl –matrix isoforms.isoform.counts.matrix –output isoforms_DE –dispersion 0.1 –method edgeR
-note, no biological replicates identified, so setting min reps = 2.
Got 4 samples, and got: 5 data fields.
Header: Sp_ds Sp_hs Sp_plat Sp_log
Next: TRINITY_DN639_c0_g1_i1 1.00 1.00 6.00 2.00
$VAR1 = { ‘Sp_ds’ => 1,
'Sp_hs' => 2,
'Sp_log' => 4,
'Sp_plat' => 3
};
$VAR1 = { ‘Sp_plat’ => [ ‘Sp_plat’ ],
'Sp_log' => [
'Sp_log'
],
'Sp_ds' => [
'Sp_ds'
],
'Sp_hs' => [
'Sp_hs'
]
};
Contrasts to perform are:
$VAR1 = [ [ ‘Sp_ds’, ‘Sp_hs’ ],
[
'Sp_ds',
'Sp_log'
],
[
'Sp_ds',
'Sp_plat'
],
[
'Sp_hs',
'Sp_log'
],
[
'Sp_hs',
'Sp_plat'
],
[
'Sp_log',
'Sp_plat'
]
];
CMD: R –vanilla -q < isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_hs.Sp_ds.vs.Sp_hs.EdgeR.Rscript
CMD: R –vanilla -q < isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_log.Sp_ds.vs.Sp_log.EdgeR.Rscript
CMD: R –vanilla -q < isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_plat.Sp_ds.vs.Sp_plat.EdgeR.Rscript
CMD: R –vanilla -q < isoforms.isoform.counts.matrix.Sp_hs_vs_Sp_log.Sp_hs.vs.Sp_log.EdgeR.Rscript
CMD: R –vanilla -q < isoforms.isoform.counts.matrix.Sp_hs_vs_Sp_plat.Sp_hs.vs.Sp_plat.EdgeR.Rscript
CMD: R –vanilla -q < isoforms.isoform.counts.matrix.Sp_log_vs_Sp_plat.Sp_log.vs.Sp_plat.EdgeR.Rscript
drwxrwxr-x. 2 bdash bdash 4096 Mar 25 11:35 isoforms_DE
trinity2]$ ls -l isoforms_DE/
-rw-rw-r–. 1 bdash bdash 232175 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_hs.edgeR.count_matrix
-rw-rw-r–. 1 bdash bdash 744943 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_hs.edgeR.DE_results
-rw-rw-r–. 1 bdash bdash 56194 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_hs.edgeR.DE_results.MA_n_Volcano.pdf
-rw-rw-r–. 1 bdash bdash 1246 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_hs.Sp_ds.vs.Sp_hs.EdgeR.Rscript
-rw-rw-r–. 1 bdash bdash 236512 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_log.edgeR.count_matrix
-rw-rw-r–. 1 bdash bdash 785584 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_log.edgeR.DE_results
-rw-rw-r–. 1 bdash bdash 58462 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_log.edgeR.DE_results.MA_n_Volcano.pdf
-rw-rw-r–. 1 bdash bdash 1252 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_log.Sp_ds.vs.Sp_log.EdgeR.Rscript
-rw-rw-r–. 1 bdash bdash 238730 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_plat.edgeR.count_matrix
-rw-rw-r–. 1 bdash bdash 750585 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_plat.edgeR.DE_results
-rw-rw-r–. 1 bdash bdash 54153 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_plat.edgeR.DE_results.MA_n_Volcano.pdf
-rw-rw-r–. 1 bdash bdash 1258 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_plat.Sp_ds.vs.Sp_plat.EdgeR.Rscript
-rw-rw-r–. 1 bdash bdash 231199 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_hs_vs_Sp_log.edgeR.count_matrix
-rw-rw-r–. 1 bdash bdash 743784 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_hs_vs_Sp_log.edgeR.DE_results
-rw-rw-r–. 1 bdash bdash 55135 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_hs_vs_Sp_log.edgeR.DE_results.MA_n_Volcano.pdf
-rw-rw-r–. 1 bdash bdash 1252 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_hs_vs_Sp_log.Sp_hs.vs.Sp_log.EdgeR.Rscript
-rw-rw-r–. 1 bdash bdash 235014 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_hs_vs_Sp_plat.edgeR.count_matrix
-rw-rw-r–. 1 bdash bdash 795142 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_hs_vs_Sp_plat.edgeR.DE_results
-rw-rw-r–. 1 bdash bdash 58927 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_hs_vs_Sp_plat.edgeR.DE_results.MA_n_Volcano.pdf
-rw-rw-r–. 1 bdash bdash 1258 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_hs_vs_Sp_plat.Sp_hs.vs.Sp_plat.EdgeR.Rscript
-rw-rw-r–. 1 bdash bdash 240471 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_log_vs_Sp_plat.edgeR.count_matrix
-rw-rw-r–. 1 bdash bdash 833522 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_log_vs_Sp_plat.edgeR.DE_results
-rw-rw-r–. 1 bdash bdash 62074 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_log_vs_Sp_plat.edgeR.DE_results.MA_n_Volcano.pdf
-rw-rw-r–. 1 bdash bdash 1264 Mar 25 11:35 isoforms.isoform.counts.matrix.Sp_log_vs_Sp_plat.Sp_log.vs.Sp_plat.EdgeR.Rscript
** Found 568 features as differentially expressed.
CMD: /home/bdash/miniconda2/opt/trinity-2.5.1/Analysis/DifferentialExpression/PtR -m DE_isoforms.matrix –log2 –heatmap –min_colSums 0 –min_rowSums 0 –gene_dist euclidean –sample_dist euclidean –sample_cor_matrix –center_rows –save
CMD: R –vanilla -q < DE_isoforms.matrix.R
library(cluster)
library(Biobase)
Loading required package: BiocGenerics
Loading required package: parallel
Attaching package: ‘BiocGenerics’
………………….
Welcome to Bioconductor
…………………..
library(qvalue)
……………….
source(“/home/bdash/miniconda2/opt/trinity-2.5.1/Analysis/DifferentialExpression/R/heatmap.3.R”)
source(“/home/bdash/miniconda2/opt/trinity-2.5.1/Analysis/DifferentialExpression/R/misc_rnaseq_funcs.R”)
source(“/home/bdash/miniconda2/opt/trinity-2.5.1/Analysis/DifferentialExpression/R/pairs3.R”)
source(“/home/bdash/miniconda2/opt/trinity-2.5.1/Analysis/DifferentialExpression/R/vioplot2.R”)
trinity2$ ls -l isoforms_DE/
-rw-rw-r–. 1 bdash bdash 107 Mar 25 11:42 DE_feature_counts.P0.001_C2.matrix
-rw-rw-r–. 1 bdash bdash 30267 Mar 25 11:42 DE_isoforms.matrix
-rw-rw-r–. 1 bdash bdash 53939 Mar 25 11:42 DE_isoforms.matrix.log2.centered.dat
-rw-rw-r–. 1 bdash bdash 36578 Mar 25 11:42 DE_isoforms.matrix.log2.centered.genes_vs_samples_heatmap.pdf
-rw-rw-r–. 1 bdash bdash 290 Mar 25 11:42 DE_isoforms.matrix.log2.centered.sample_cor.dat
-rw-rw-r–. 1 bdash bdash 5867 Mar 25 11:42 DE_isoforms.matrix.log2.centered.sample_cor_matrix.pdf
-rw-rw-r–. 1 bdash bdash 3887 Mar 25 11:42 DE_isoforms.matrix.R
-rw-rw-r–. 1 bdash bdash 1342120 Mar 25 11:42 DE_isoforms.matrix.RData
-rw-rw-r–. 1 bdash bdash 10031 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_hs.edgeR.DE_results.P0.001_C2.DE.subset
-rw-rw-r–. 1 bdash bdash 5830 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_hs.edgeR.DE_results.P0.001_C2.Sp_ds-UP.subset
-rw-rw-r–. 1 bdash bdash 4268 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_hs.edgeR.DE_results.P0.001_C2.Sp_hs-UP.subset
-rw-rw-r–. 1 bdash bdash 24 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_hs.edgeR.DE_results.samples
-rw-rw-r–. 1 bdash bdash 28526 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_log.edgeR.DE_results.P0.001_C2.DE.subset
-rw-rw-r–. 1 bdash bdash 19319 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_log.edgeR.DE_results.P0.001_C2.Sp_ds-UP.subset
-rw-rw-r–. 1 bdash bdash 9274 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_log.edgeR.DE_results.P0.001_C2.Sp_log-UP.subset
-rw-rw-r–. 1 bdash bdash 26 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_log.edgeR.DE_results.samples
-rw-rw-r–. 1 bdash bdash 6892 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_plat.edgeR.DE_results.P0.001_C2.DE.subset
-rw-rw-r–. 1 bdash bdash 2225 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_plat.edgeR.DE_results.P0.001_C2.Sp_ds-UP.subset
-rw-rw-r–. 1 bdash bdash 4734 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_plat.edgeR.DE_results.P0.001_C2.Sp_plat-UP.subset
-rw-rw-r–. 1 bdash bdash 28 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_ds_vs_Sp_plat.edgeR.DE_results.samples
-rw-rw-r–. 1 bdash bdash 21828 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_hs_vs_Sp_log.edgeR.DE_results.P0.001_C2.DE.subset
-rw-rw-r–. 1 bdash bdash 17739 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_hs_vs_Sp_log.edgeR.DE_results.P0.001_C2.Sp_hs-UP.subset
-rw-rw-r–. 1 bdash bdash 4156 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_hs_vs_Sp_log.edgeR.DE_results.P0.001_C2.Sp_log-UP.subset
-rw-rw-r–. 1 bdash bdash 26 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_hs_vs_Sp_log.edgeR.DE_results.samples
-rw-rw-r–. 1 bdash bdash 35385 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_hs_vs_Sp_plat.edgeR.DE_results.P0.001_C2.DE.subset
-rw-rw-r–. 1 bdash bdash 17568 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_hs_vs_Sp_plat.edgeR.DE_results.P0.001_C2.Sp_hs-UP.subset
-rw-rw-r–. 1 bdash bdash 17884 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_hs_vs_Sp_plat.edgeR.DE_results.P0.001_C2.Sp_plat-UP.subset
-rw-rw-r–. 1 bdash bdash 28 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_hs_vs_Sp_plat.edgeR.DE_results.samples
-rw-rw-r–. 1 bdash bdash 52861 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_log_vs_Sp_plat.edgeR.DE_results.P0.001_C2.DE.subset
-rw-rw-r–. 1 bdash bdash 19495 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_log_vs_Sp_plat.edgeR.DE_results.P0.001_C2.Sp_log-UP.subset
-rw-rw-r–. 1 bdash bdash 33433 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_log_vs_Sp_plat.edgeR.DE_results.P0.001_C2.Sp_plat-UP.subset
-rw-rw-r–. 1 bdash bdash 30 Mar 25 11:42 isoforms.isoform.counts.matrix.Sp_log_vs_Sp_plat.edgeR.DE_results.samples
library(pdftools)
pdf_convert('isoforms_DE/DE_isoforms.matrix.log2.centered.genes_vs_samples_heatmap.pdf', format = "jpeg", pages = NULL, filenames = NULL, dpi = 300, antialias = TRUE, opw = "", upw = "", verbose = TRUE)
Converting page 1 to DE_isoforms.matrix.log2.centered.genes_vs_samples_heatmap_1.jpeg... done!
[1] "DE_isoforms.matrix.log2.centered.genes_vs_samples_heatmap_1.jpeg"
DE_genes.heatmap
pdf_convert('isoforms_DE/DE_isoforms.matrix.log2.centered.sample_cor_matrix.pdf', format = "jpeg", pages = NULL, filenames = NULL, dpi = 300, antialias = TRUE, opw = "", upw = "", verbose = TRUE)
Converting page 1 to DE_isoforms.matrix.log2.centered.sample_cor_matrix_1.jpeg... done!
[1] "DE_isoforms.matrix.log2.centered.sample_cor_matrix_1.jpeg"
DE_genes.cor_matrix