RNA-seq identifies a distinct response in macrophages infected with Mycobacterium tuberculosis for a serial time points.

Purpose: The goals of this study are to obtain the NGS-derived transcriptome profiling (RNA-seq) for THP-1 macrophages response to Mycobacterium tuberculosis (H37Rv and H37Ra) Methods: mRNA and long noncoding RNA profiles of THP-1 macrophages infected with H37Rv and H37Ra for 1, 4, 12, 24, 48 hours were generated by deep sequencing, using Illumina Hiseq3000. The sequence reads that passed quality filters were first mapped to the latest UCSC transcript set using Bowtie2 (version 2.1.0). Then the gene expression level was estimated using RSEM (RNA-Seq by Expectation Maximization, v1.2.15), for lncRNA analysis, reads were mapped to lncRNA transcript set from LNCipedia.org. The sequence reads were normalized with TMM (trimmed mean of M-values) to identify differentially expressed genes (DEGs) using the edgeR package edgeR. qRT–PCR validation was performed using SYBR Green assays. Results: Using an optimized data analysis workflow, we mapped about 20 million sequence reads per sample to the human genome (GRCh38/hg38) and identified 25,343 mRNA and 47877 long non-coding RNA transcripts. Conclusions: Our study represents the detailed analysis of transcriptomes for THP-1 macrophages response to H37Rv and H37Ra, generated by RNA-seq technology. NGS-derived transcriptome profiling (RNA-seq) for THP-1 macrophages response to Mycobacterium tuberculosis H37Rv and H37Ra.