tb /

Filename Size Date modified Message
78 B
55.0 KB
3.1 KB
7.8 KB
12.0 KB
825 B
1.2 KB
4.2 KB
36.7 KB
10.2 KB
211 B
10.9 KB
2.0 KB
26.0 KB
2.3 KB
8.0 KB
26.4 KB
65 B
589 B
1.0 KB
7.0 KB
8.4 KB
2.7 KB
3.1 KB
2.6 KB
6.6 KB
4.1 KB
7.5 KB
11.2 KB
4.1 KB
1.1 KB
2.7 KB

Mycobacterial infection induces a specific human innate immune response

This repo contains the code for the following manuscript:

Blischak, J. D. et al. Mycobacterial infection induces a specific human innate immune response. Sci. Rep. 5, 16882; doi: 10.1038/srep16882 (2015).

Other links related to this publication:


We have collected gene expression data from macrophages isolated from six individuals. They were infected with different strains of Mycobacterium tuberculosis and other bacteria (nine total "conditions" including infections and controls). Furthermore, we collected data at 4, 18, and 48 hours post-infection.

Main analysis files

  • differential-expression.Rmd - Performs DE analysis with limma/voom.

  • joint_analysis.Rmd - Runs Cormotif to classify genes according to their expression patterns across the infections.

  • candidate_genes.Rmd - Plots expression of immune response genes.

  • reQTL.Rmd - Compares FDRs from DE analysis of response eQTL-associated genes (identified in Barreiro, Tailleux et al. 2012) to the background of expressed genes.

Preprocessing steps

Uses R package ngspipeline.

Sort files

./organize_files.R ../data/flow-cells.csv ../samples/

Map reads

find ../samples/ -name "*fastq.gz" -exec bash -c \
'echo "./map_fastq.R {} ~/subread-indexed-genomes/hg19/hg19" | \
qsub -l h_vmem=16g,bigio=1 -V -N map-`basename {}` -cwd -o stdout -e stderr' \;

Count reads per gene

mkdir std-counts
find ../samples/ -name "*bam" -exec bash -c \
'echo "./count_reads.R -sl ../data/gene_lengths.txt {} ../data/exons.txt -o /KG/jdblischak/counts" | \
qsub -l h_vmem=8g -V -N count-`basename {}` -cwd -o std-counts -j y' \;

Combine annotation and gene counts

anno  <-  "../data/annotation.txt"
count_files <- list.files("/KG/jdblischak/counts/",
                          pattern = "counts$", full.name = TRUE)
out_name <- "../data/counts_per_run.txt"
add_anno_to_counts(anno, count_files, out_name)

Sum gene counts per sample



The code is available via the GPLv3 license (because we use multiple R/Biocodncutor packages which are GPL'd). The text is available via the CC BY 3.0 license.