gene_x 0 like s 410 view s
Tags: processing, bash, pipeline, RNA-seq
YopM is a systematic antagonist of YopP in suppressing inflammatory gene expression in Yersinia-infected human macrophages
goal
# 0th heatmap see manuscript: The
# 1.th heatmap mock WAC WAP 1.5 + 6 hours: WAC, WA314 comparing to mock
# 2.th heatmap deltaYopM, deltaYopP, deltaYopQ Vergleich zu WAP (6 hours). The sample deltaYopMP is completed deleted: deltaM, deltaP, deltaQ to WA314
download the batch 1 data
#-- PRJEB10086 --
mock_6h_a
mock_6h_b
mock_90min_a
mock_90min_b
WA314_6h_a
WA314_6h_b
WA314_90min_a
WA314_90min_b
WA314dYopM_6h_a
WA314dYopM_6h_b
WA314dYopM_90min_a
WA314dYopM_90min_b
#-- PRJEB45780 --
wget ftp.sra.ebi.ac.uk/vol1/run/ERR592/ERR5924425/mock_DoA.fastq.gz
wget ftp.sra.ebi.ac.uk/vol1/run/ERR592/ERR5924426/mock_DoII.fastq.gz
#mv mock_DoA.fastq.gz mock_6h_DoA.fastq.gz
#mv mock_DoII.fastq.gz mock_6h_DoII.fastq.gz
ftp.sra.ebi.ac.uk/vol1/run/ERR592/ERR5924427/WAC1_5h_DoA.fastq.gz
ftp.sra.ebi.ac.uk/vol1/run/ERR592/ERR5924428/WAC1_5h_DoII.fastq.gz
ftp.sra.ebi.ac.uk/vol1/run/ERR592/ERR5924429/WAC6h_DoA.fastq.gz
ftp.sra.ebi.ac.uk/vol1/run/ERR592/ERR5924430/WAC6h_DoII.fastq.gz
ftp.sra.ebi.ac.uk/vol1/run/ERR592/ERR5924431/WAP1_5h_DoA.fastq.gz
ftp.sra.ebi.ac.uk/vol1/run/ERR592/ERR5924421/WAP1_5h_DoII.fastq.gz
ftp.sra.ebi.ac.uk/vol1/run/ERR592/ERR5924422/WAP6h_DoA.fastq.gz
ftp.sra.ebi.ac.uk/vol1/run/ERR592/ERR5924423/WAP6h_DoII.fastq.gz
ftp.sra.ebi.ac.uk/vol1/run/ERR592/ERR5924420/dP6h_DoII.fastq.gz
ftp.sra.ebi.ac.uk/vol1/run/ERR592/ERR5924424/dP6h_DoIII.fastq.gz
#-- PRJEB44958 --
ftp.sra.ebi.ac.uk/vol1/run/ERR608/ERR6087955/dYopP_1.5h_DonorII.fastq.gz
ftp.sra.ebi.ac.uk/vol1/run/ERR608/ERR6087956/dYopP_1.5h_DonorIII.fastq.gz
wget ftp.sra.ebi.ac.uk/vol1/run/ERR608/ERR6087951/dYopMP_1.5h_DonorII.fastq.gz
wget ftp.sra.ebi.ac.uk/vol1/run/ERR608/ERR6087952/dYopMP_1.5h_DonorIII.fastq.gz
wget ftp.sra.ebi.ac.uk/vol1/run/ERR608/ERR6087953/dYopMP_6h_DonorII.fastq.gz
wget ftp.sra.ebi.ac.uk/vol1/run/ERR608/ERR6087954/dYopMP_6h_DonorIII.fastq.gz
# ---- 1st. batch: 1.5h + 6h ----
mock
WAC
WA314 == WAP
dYopM,dYopP
E-MTAB-10602
E-MTAB-10473
Project: PRJEB45780
Pathogenic bacteria Yersinia enterocolitica injects virulence plasmid-encoded effectors through the type three secretion system into macrophages to modulate gene expression. Here we analyzed the effect on gene expression in primary human macrophages of Y. enterocolitica strains lacking effector YopP (1.5 h infection) or effectors YopP and YopM (1.5 h or 6 h infection) simultaneously using RNA-seq. This is part of a larger sequencing experiment for which other samples can be found in EMBL-EBI (www.ebi.ac.uk/arrayexpress) under accession number E-MTAB-10473 and European Nucleotide Archive (ENA) at http://www.ebi.ac.uk/ena/data/view/PRJEB10086.
./ena-file-download-read_run-PRJEB10086-submitted_ftp-20231220-1619_sorted.sh
#? check WA314_90min_a.fastq.gz == WAP1_5h_DoA.fastq.gz
#? check WA314_90min_b.fastq.gz == WAP1_5h_DoII.fastq.gz
#? check WA314_6h_a.fastq.gz == WAP6h_DoA.fastq.gz
#? check WA314_6h_b.fastq.gz == WAP6h_DoII.fastq.gz
#TODO: if they are not the same, take WA314_* as the WAP-samples in the 1st batch!
cat mock_90min_a_L004_R1_001.fastq.gz mock_90min_a_L005_R1_001.fastq.gz mock_90min_a_L006_R1_001.fastq.gz mock_90min_a_L007_R1_001.fastq.gz mock_90min_a_L008_R1_001.fastq.gz > mock_90min_a.fastq.gz
cat mock_90min_b_L004_R1_001.fastq.gz mock_90min_b_L005_R1_001.fastq.gz mock_90min_b_L006_R1_001.fastq.gz mock_90min_b_L007_R1_001.fastq.gz mock_90min_b_L008_R1_001.fastq.gz > mock_90min_b.fastq.gz
cat mock_6h_a_L004_R1_001.fastq.gz mock_6h_a_L005_R1_001.fastq.gz mock_6h_a_L006_R1_001.fastq.gz mock_6h_a_L007_R1_001.fastq.gz mock_6h_a_L008_R1_001.fastq.gz > mock_6h_a.fastq.gz
cat mock_6h_b_L004_R1_001.fastq.gz mock_6h_b_L005_R1_001.fastq.gz mock_6h_b_L006_R1_001.fastq.gz mock_6h_b_L007_R1_001.fastq.gz mock_6h_b_L008_R1_001.fastq.gz > mock_6h_b.fastq.gz
cat WA314_90min_a_L004_R1_001.fastq.gz WA314_90min_a_L005_R1_001.fastq.gz WA314_90min_a_L006_R1_001.fastq.gz WA314_90min_a_L007_R1_001.fastq.gz WA314_90min_a_L008_R1_001.fastq.gz > WA314_90min_a.fastq.gz
cat WA314_90min_b_L004_R1_001.fastq.gz WA314_90min_b_L005_R1_001.fastq.gz WA314_90min_b_L006_R1_001.fastq.gz WA314_90min_b_L007_R1_001.fastq.gz WA314_90min_b_L008_R1_001.fastq.gz > WA314_90min_b.fastq.gz
cat WA314_6h_a_L004_R1_001.fastq.gz WA314_6h_a_L005_R1_001.fastq.gz WA314_6h_a_L006_R1_001.fastq.gz WA314_6h_a_L007_R1_001.fastq.gz WA314_6h_a_L008_R1_001.fastq.gz > WA314_6h_a.fastq.gz
cat WA314_6h_b_L004_R1_001.fastq.gz WA314_6h_b_L005_R1_001.fastq.gz WA314_6h_b_L006_R1_001.fastq.gz WA314_6h_b_L007_R1_001.fastq.gz WA314_6h_b_L008_R1_001.fastq.gz > WA314_6h_b.fastq.gz
cat WA314dYopM_90min_a_L004_R1_001.fastq.gz WA314dYopM_90min_a_L005_R1_001.fastq.gz WA314dYopM_90min_a_L006_R1_001.fastq.gz WA314dYopM_90min_a_L007_R1_001.fastq.gz WA314dYopM_90min_a_L008_R1_001.fastq.gz > WA314dYopM_90min_a.fastq.gz
cat WA314dYopM_90min_b_L004_R1_001.fastq.gz WA314dYopM_90min_b_L005_R1_001.fastq.gz WA314dYopM_90min_b_L006_R1_001.fastq.gz WA314dYopM_90min_b_L007_R1_001.fastq.gz WA314dYopM_90min_b_L008_R1_001.fastq.gz > WA314dYopM_90min_b.fastq.gz
cat WA314dYopM_6h_a_L004_R1_001.fastq.gz WA314dYopM_6h_a_L005_R1_001.fastq.gz WA314dYopM_6h_a_L006_R1_001.fastq.gz WA314dYopM_6h_a_L007_R1_001.fastq.gz WA314dYopM_6h_a_L008_R1_001.fastq.gz > WA314dYopM_6h_a.fastq.gz
cat WA314dYopM_6h_b_L004_R1_001.fastq.gz WA314dYopM_6h_b_L005_R1_001.fastq.gz WA314dYopM_6h_b_L006_R1_001.fastq.gz WA314dYopM_6h_b_L007_R1_001.fastq.gz WA314dYopM_6h_b_L008_R1_001.fastq.gz > WA314dYopM_6h_b.fastq.gz
#END
./WAC1_5h_DoA.fastq.gz
./WAC1_5h_DoII.fastq.gz
./WAC6h_DoA.fastq.gz
./WAC6h_DoII.fastq.gz
./WAP1_5h_DoA.fastq.gz
./WAP1_5h_DoII.fastq.gz
./WAP6h_DoA.fastq.gz
./WAP6h_DoII.fastq.gz
./dYopP_1.5h_DonorII.fastq.gz
./dYopP_1.5h_DonorIII.fastq.gz
./dP6h_DoII.fastq.gz
./dP6h_DoIII.fastq.gz
mv dP6h_DoII.fastq.gz dYopP_6h_DoII.fastq.gz
mv dP6h_DoIII.fastq.gz dYopP_6h_DoIII.fastq.gz
mv WAC1_5h_DoA.fastq.gz WAC_1.5h_DoA.fastq.gz
mv WAC1_5h_DoII.fastq.gz WAC_1.5h_DoII.fastq.gz
mv WAC6h_DoA.fastq.gz WAC_6h_DoA.fastq.gz
mv WAC6h_DoII.fastq.gz WAC_6h_DoII.fastq.gz
mv WAP1_5h_DoA.fastq.gz WAP_1.5h_DoA.fastq.gz
mv WAP1_5h_DoII.fastq.gz WAP_1.5h_DoII.fastq.gz
mv WAP6h_DoA.fastq.gz WAP_6h_DoA.fastq.gz
mv WAP6h_DoII.fastq.gz WAP_6h_DoII.fastq.gz
mv WA314dYopM_90min_a.fastq.gz dYopM_90min_a.fastq.gz
mv WA314dYopM_90min_b.fastq.gz dYopM_90min_b.fastq.gz
mv WA314dYopM_6h_a.fastq.gz dYopM_6h_a.fastq.gz
mv WA314dYopM_6h_b.fastq.gz dYopM_6h_b.fastq.gz
#Desperated!
#wget -nc ftp://ftp.sra.ebi.ac.uk/vol1/run/ERR592/ERR5924425/mock_DoA.fastq.gz
#wget -nc ftp://ftp.sra.ebi.ac.uk/vol1/run/ERR592/ERR5924426/mock_DoII.fastq.gz
#wget -nc ftp://ftp.sra.ebi.ac.uk/vol1/run/ERR608/ERR6087953/dYopMP_6h_DonorII.fastq.gz
#wget -nc ftp://ftp.sra.ebi.ac.uk/vol1/run/ERR608/ERR6087954/dYopMP_6h_DonorIII.fastq.gz
#wget -nc ftp://ftp.sra.ebi.ac.uk/vol1/run/ERR608/ERR6087951/dYopMP_1.5h_DonorII.fastq.gz
#wget -nc ftp://ftp.sra.ebi.ac.uk/vol1/run/ERR608/ERR6087952/dYopMP_1.5h_DonorIII.fastq.gz
download the public data and prepare the batch 2 data
# ---- data public (naive and LPS is public data)----
jhuang@hamburg:/media/jhuang/Elements2/Data_Indra_RNASeq_GSM2262901/raw_data
#Dataset name GEO ID Series Experiment
#- naive macrophages GSM2679941 GSE100382 RNA-seq (RNA_N_rep1)
# SRR5746425 4,000,000 204M 136Mb 2017-06-25
# SRR5746426 4,000,000 204M 137Mb 2017-06-25
# SRR5746427 4,000,000 204M 137.2Mb 2017-06-25
# SRR5746428 4,000,000 204M 135.6Mb 2017-06-25
# SRR5746429 4,000,000 204M 134.5Mb 2017-06-25
# SRR5746430 4,000,000 204M 136.2Mb 2017-06-25
# SRR5746431 903,661 46.1M 32.7Mb 2017-06-25
cp RNA_N_rep1.fastq.gz ../raw_data/naive_rep1.fastq.gz
#- naive macrophages GSM2679942 GSE100382 RNA-seq (RNA_N_rep2)
# SRR5746432 14,176,210 723M 350.2Mb 2017-06-25
# SRR5746433 14,157,173 722M 353.4Mb 2017-06-25
cp RNA_N_rep2.fastq.gz ../raw_data/naive_rep2.fastq.gz
#- LPS stimulated macrophages GSM2679944 GSE100382 RNA-seq (RNA_L_rep1)
# SRR5746436 4,000,000 204M 135.8Mb 2017-06-25
# SRR5746437 4,000,000 204M 137.1Mb 2017-06-25
# SRR5746438 4,000,000 204M 136.6Mb 2017-06-25
# SRR5746439 4,000,000 204M 135.4Mb 2017-06-25
# SRR5746440 4,000,000 204M 135.3Mb 2017-06-25
# SRR5746441 2,178,850 111.1M 76.7Mb 2017-06-25
cp RNA_L_rep1.fastq.gz ../raw_data/LPS_rep1.fastq.gz
#- LPS stimulated macrophages GSM2679945 GSE100382 RNA-seq (RNA_L_rep2)
# SRR5746442 17,351,249 884.9M 428Mb 2017-06-25
# SRR5746443 17,318,750 883.3M 431.8Mb 2017-06-25
cp RNA_L_rep2.fastq.gz ../raw_data/LPS_rep2.fastq.gz
#- naive macrophages GSM2262901 GSE85243 RNA-seq (RPMI_d6_6095)
# SRR4004025 11,559,434 485.5M 338.3Mb 2016-11-21
cp RPMI_d6_6095.fastq.gz ../raw_data/naive_rep3.fastq.gz
#- naive macrophages GSM2262902 GSE85243 RNA-seq (RPMI_d6_6718)
# SRR4004026 22,911,696 962.3M 635Mb 2016-11-21
cp RPMI_d6_6718.fastq.gz ../raw_data/naive_rep4.fastq.gz
#- LPS stimulated macrophages GSM2262906 GSE85243 RNA-seq (RPMI_d6_Restim_8749)
# SRR4004030 26,566,129 1.1G 681.3Mb 2016-11-21
cp RPMI_d6_Restim_8749.fastq.gz ../raw_data/LPS_rep3.fastq.gz
#- LPS stimulated macrophages GSM2262907 GSE85243 RNA-seq (RPMI_d6_Restim_8754)
# SRR4004031 28,719,690 1.2G 731.2Mb 2016-11-21
cp RPMI_d6_Restim_8754.fastq.gz ../raw_data/LPS_rep4.fastq.gz
#END
# ---- 2nd. batch: 6h jhuang@hamburg:~/DATA/Data_Soeren_RNA-seq_2022/Raw_Data ----
mock
WAC
WAP
deltaYopQ -->
Vergleich zu WAP
mv MOCK_A.fastq.gz ../raw_data/
mv MOCK_B.fastq.gz ../raw_data/
mv WAC_A.fastq.gz ../raw_data/
mv WAC_B.fastq.gz ../raw_data/
mv WAP_A.fastq.gz ../raw_data/
mv WAP_B.fastq.gz ../raw_data/
mv deltaQ_A.fastq.gz ../raw_data/
mv deltaQ_B.fastq.gz ../raw_data/
run nextflow
(rnaseq) [jhuang@sage DATA]$ pwd
/home/jhuang/DATA
/home/jhuang/REFs/Homo_sapiens/Ensembl
#"/home/jhuang/REFs/Homo_sapiens/UCSC/hg38/blacklists/hg38-blacklist.bed"
nextflow run rnaseq/main.nf --input samplesheet.csv --outdir results_GRCh38 --genome GRCh38 -profile test_full -resume --max_memory 300.GB --max_time 2400.h --save_reference --aligner star_salmon --skip_deseq2_qc --skip_fastqc
ln -s ~/Tools/rnaseq/assets/multiqc_config.yaml multiqc_config.yaml
multiqc -f --config multiqc_config.yaml . 2>&1
rm multiqc_config.yaml
summary:
#ENSG00000000419 DPM1 385.616
#grep "ENSG00000000419" quant.genes.sf
#ENSG00000000419 1159.45 909.453 15.8259 385.616
Attached, you will find the results of the RNASeq analysis discussed in our last online meeting.
The raw counts are located in Yersinia_RNAseq_results_GRCh38/star_salmon/salmon.merged.gene_counts.tsv. Note that this table was generated using the STAR-Salmon strategy, with STAR aligning the reads and Salmon quantifying transcript abundances. Salmon's advanced probabilistic model assigns fractional counts to transcripts, particularly when reads map to multiple similar transcripts. This is the reason non-integer numbers appear in this table.
The samples included in the analysis are as follows:
batch 1:
./mock_90min_a.fastq.gz
./mock_90min_b.fastq.gz
./mock_6h_a.fastq.gz
./mock_6h_b.fastq.gz
./WAC_1.5h_DoA.fastq.gz
./WAC_1.5h_DoII.fastq.gz
./WAC_6h_DoA.fastq.gz
./WAC_6h_DoII.fastq.gz
./WAP_1.5h_DoA.fastq.gz
./WAP_1.5h_DoII.fastq.gz
./WAP_6h_DoA.fastq.gz
./WAP_6h_DoII.fastq.gz
./WA314_90min_a.fastq.gz
./WA314_90min_b.fastq.gz
./WA314_6h_a.fastq.gz
./WA314_6h_b.fastq.gz
./dYopM_90min_a.fastq.gz
./dYopM_90min_b.fastq.gz
./dYopM_6h_a.fastq.gz
./dYopM_6h_b.fastq.gz
./dYopP_1.5h_DonorII.fastq.gz
./dYopP_1.5h_DonorIII.fastq.gz
./dYopP_6h_DoII.fastq.gz
./dYopP_6h_DoIII.fastq.gz
./mock_6h_DoA.fastq.gz
./mock_6h_DoII.fastq.gz
./dYopMP_1.5h_DonorII.fastq.gz
./dYopMP_1.5h_DonorIII.fastq.gz
./dYopMP_6h_DonorII.fastq.gz
./dYopMP_6h_DonorIII.fastq.gz
batch public:
naive_rep1.fastq.gz
naive_rep2.fastq.gz
LPS_rep1.fastq.gz
LPS_rep2.fastq.gz
naive_rep3.fastq.gz
naive_rep4.fastq.gz
LPS_rep3.fastq.gz
LPS_rep4.fastq.gz
batch 2:
MOCK_A.fastq.gz (6h)
MOCK_B.fastq.gz (6h)
WAC_A.fastq.gz (6h)
WAC_B.fastq.gz (6h)
WAP_A.fastq.gz (6h)
WAP_B.fastq.gz (6h)
deltaQ_A.fastq.gz (6h)
deltaQ_B.fastq.gz (6h)
点赞本文的读者
还没有人对此文章表态
没有评论
RNA-seq 2024 Ute from raw counts
© 2023 XGenes.com Impressum