Metabolomics,Unknown,Transcriptomics,Genomics,Proteomics

Dataset Information

The Time-Series Transcriptomic Responses of THP-1 Cells to W-Beijing M. tuberculosis Strains of different sublineages

ABSTRACT: We report a whole-genome expression profiling of PMA-treated THP-1 human macrophages as a host system, infected by W-Beijing strains of different sublineages or the H37Rv reference strain. Surprisingly, we found that host transcriptional responses were irrespective of inter-sublineage variations. Based on this finding, we further showed that such common host transcriptional responses were probably under the control of a regulatory network involving STATs, IRF-1, IRF-7 and Oct-1 transcriptional factors. We also suggested that these putative regulators might act in a cooperative manner to induce the expression of host-responsive genes, particularly those involved in cytokine-cytokine receptor interactions. THP-1 cells were infected by 11 different sublineages plus one reference strain H37Rv (two replications) at three time points (i.e., 4h, 18h and 48h). For each infection as well as a non-infection control, RNA was extracted using Trizol (Life Technologies)TRIZOL® reagent (Invitrogen), and subsequently amplified and labeled with biotin according to the standard Affymetrix® protocol. Specifically, RNA was DNase-treated and purified with RNeasy columns (Qiagen). The purified RNA was quantified using the ND-1000 Spectrophotometer (NanoDrop Technologies) and Agilent 2100 Bioanalyzer (Agilent Technologies). The biotinylated cRNA was then fragmented and subjected to hybridization on Human Genome U133 Plus 2.0 Array (Affymetrix, USA). The scanned images were converted to cell intensity files (CEL) using GeneChip® Operating Software (GCOS) (Affymetrix). These CEL raw expression data were normalized using Robust Multi-array Averaging (RMA) with quantile normalization in R (Bioconductor). Detection call-based filter was applied to remove all the probe-sets whose expression values were consistently below an empirically-determined value of minimum sensitivity, which were evaluated according to the 95th percentile of all the ‘Absent’ call-flagged signals of the entire dataset. Following the normalization and filtering, a gene expression matrix (termed as TB expression matrix) was constructed. The TB expression matrix was then subjected to component plane presentation integrated self-organizing map (CPP-SOM) for visual comparisons, or subjected to Cluster3.0/TreeView-1.0.8 softwares with Euclidean distance algorithm for sample classification of 13 time-course transcriptome profiles. Also, Linear Models for Microarray Data (LIMMA) bioconductor library[18] was applied to identify those differentially expressed genes between any two successive time points. The criteria for identifying the top significant genes for the designed contrast was based on P-values (< 0.01) corrected using Benjamini and Hochberg false discovery rate (FDR) procedure. Sets of genes of interest such as differentially expressed genes between 4h and 18h infections were uploaded onto Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.7 for identifying enriched gene ontology and KEGG pathways (FDR<0.01). For TFBS enrichment analysis, we applied PRIMA (PRomoter Integration in Microarray Analysis) program to identify TFs whose binding sites are significantly over-represented in a given gene set (i.e., differentially expressed genes) compared to background. Bonferroni-corrected P-values (< 0.01) were utilized to determine the significance of the binding site enrichments of TFs. Such computationally predicted TFs were considered as candidate regulators responsible for the putative co-regulations of differentially expressed genes.

ORGANISM(S): Homo sapiens

SUBMITTER: Hai Fang

PROVIDER: E-GEOD-29628 | biostudies-arrayexpress |

REPOSITORIES: biostudies-arrayexpress

ACCESS DATA

Similar Datasets

Project description:In order to determine BCL6 target genes an EBV negative Burkitt's lymphoma cell line, DG75, was stably transfected with a tetracycline transactivator and tight doxycycline responsive expression of GFP was established. The endogenous BCL6 genes of this cell line were disrupted by homologous recombination and a BCL6 cDNA downstream of tetracycline responsive elements (TRE) was inserted to produce Bcl6-/-:tetBCL6-HA cells. Westerns demonstrated doxycycline dependent BCL6 expression.Bcl6-/-:tet. BCL6-HA cells (clone AB7) were either grown without doxycycline (control) or with 1 ug/ml doxycycline for 16, 48 or 96 hours. Total RNA was extracted using RNeasy minipreps (Qiagen) and concentration and quality were checked on the NanoDrop® ND- 1000 spectrophotometer (NanoDrop Technologies, USA) and the RNA Nano 6000 kit (Agilent Technologies) on a 2100 Bioanalyzer (Agilent Technologies). One hundred ng of total RNA was processed with the GeneChip® Eukaryotic Whole Transcript Sense Target Labelling Assay kit (Affymetrix) according to the manufacturer's details. Hybridisation and scanning of GeneChips® was carried out at the CSC/IC Microarray Centre, MRC Clinical Sciences Centre Imperial College London and data analysis by Bioinformatics Support Service, Imperial College London. Briefly, pre- processing of data was performed using GeneSpring GX 10.0.2 software (Agilent Technologies) which applied the "Exon RMA16"¸ algorhithm to the data set. Exon RMA16 performs background correction, quantile normalisation, median polish summarisation and variance stabilisation of 16. In background correction, intensity values of each individual array are corrected for non-specific binding by subtracting the average signal intensity of the area between spots from each probe set. Normalisation is required so multiple chips can be compared to each other. Quantile normalisation adjusts the distribution of probe intensity of each array analysed and so that the distribution of probe intensities for each array in a set of arrays is the same. Probe summarisation refers to the conversion of probe level values (there are approximately 26 probes per gene on each GeneChip®) to a single probe set expression value. Variance stabilisation of 16 refers to the addition of the value 16 to the expression values. By increasing the expression value, the variance of the data set is reduced and the distribution (defined by its mean and its variance) is stabilised.

Project description:A network of transcription factors (TFs) determines cell identity, but identity can be altered by overexpressing a combination of TFs. In principle, this opens the possibility of achieving one of the goals of regenerative medicine generating the desired differentiated cells from pluripotent stem cells, such as embryonic stem (ES) cells and induced pluripotent stem (iPS) cells. However, choosing and verifying combinations of TFs for specific cell differentiation has been daunting due to a large number of possible combinations of ~2,000 TFs. Here we report an unbiased systematic approach identifying individual TFs that can direct efficient and rapid specific lineage differentiation. We start from a correlation matrix of global gene expression responses to the induction of single TFs and global gene expression profiles of a variety of tissues and organs. Based on the correlation matrix, we select TFs as examples and show that their overexpression differentiates ES cells into cells of specific organs, as predicted: Sfpi1 for blood cells, Hnf4a or Foxa1 for hepatocytes, and Ascl1 for neurons. Furthermore, we show that transfection of synthetic mRNAs of Sfpi1, Hnf4a, or Ascl1 generate correspondingly specific target cells. These results demonstrate both the wide-ranging utility of this approach to identify potent TFs for cell differentiation, and also the unanticipated capacity of single TFs to directly guide differentiation to specific lineage fates. Previously, we generated the global gene expression profiles obtained by overexpressing single TFs using the NIA mouse ES cell bank, which consists of 137 mouse ES cell lines carrying a doxycycline-regulatable TF. We then generated the correlation matrix by comparing TF-induced gene expression profiles and the expression profiles of a variety of cell types. TFs predicted by the correlation matrix for specific cell differentiation were selected and subjected to the detailed differentiation assays. ES cell lines were cultured in the GMEM with 1% FBS, 10% knockout serum (invitrogen), and LIF (Millipore). To induce differentiation, ES cell lines were cultured in the differentiation medium (DM) [(alpha minimal essential medium, aMEM) supplemented with 10% fetal calf serum and 5 × 10–5 M 2-mercaptoethanol] with or without doxycycline (1 g/ml). The following organ-specific cell culture media were also used: StemPro-34 Serum Free Media (invitrogen) for blood cells, Hepatocyte culture media kit (BD Biosciences) for hepatocytes, and NeuroCult differentiation kit (StemCell Technologies Inc) for neurons.

Dataset Information

The Time-Series Transcriptomic Responses of THP-1 Cells to W-Beijing M. tuberculosis Strains of different sublineages

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure