I-TRAP: A method to identify transcriptional regulator activated promoters

BACKGROUND: The differential expression of virulence genes is often used by microbial pathogens in adapting to the environment of their host. The differential expression of such sets of genes can be regulated by RNA polymerase sigma factors. Some sigma factors are differentially expressed, which can provide a means to identifying other differentially expressed genes such as those whose expression are controlled by the sigma factor. METHODS: To identify sigma factor-regulated genes, we developed a method, termed I-TRAP, for the identification of transcriptional regulator activated promoters. The I-TRAP method is based on the fact that some genes will be differentially expressed in the presence and absence of a transcriptional regulator. I-TRAP uses a DNA library in a promoter-trap vector that contains two reporter genes, one to allow the selection of active promoters in the presence of the transcriptional regulator and a second to allow screening for promoter activity in the absence of the transcriptional regulator. RESULTS: To illustrate the development and use of the I-TRAP approach, the construction of the vectors, host strains, and library necessary to identify SigmaE-regulated genes of Mycobacterium tuberculosis is described. CONCLUSION: The I-TRAP method should be a versatile and useful method for identifying and characterizing promoter activity under a variety of conditions and in response to various regulatory proteins. In our study, we isolated 360 clones that may contain plasmids carrying SigmaE-regulated promoters genes of M. tuberculosis

Screening hub genes in coronary artery disease based on integrated analysis

Background: Coronary artery disease (CAD) is the leading cause of mortality worldwide. Identifying key pathogenic genes benefits the understanding molecular mechanism of CAD. Methods: In this study, 5 microarray data sets from the blood sample of 312 CADs and 277 healthy controls were downloaded. Limma and metaMA packages were used to identify differentially expressed genes. The functional enrichment analysis of differentially expressed genes was further performed by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes. Additionally, protein–protein interac­tion and transcript factors-target networks were performed based on top 10 up- and down-regulated differentially expressed genes to further study the biological function. Last, real-time quantitative poly­merase chain reaction (RT-qPCR) was used to validate the integrated analysis result. Results: A total of 528 differentially expressed genes were obtained. All differentially expressed genes were significantly involved in signal transduction and the MAPK signaling pathway. Among MAPK signaling pathway, IL1R2, ARRB2 and PRKX were associated with CAD. Furthermore, there were 4 common differentially expressed genes including PLAUR, HSPH1, ZMYND11 and S100A8 in the protein–protein interaction and transcript factors-target networks, which played crucial roles in the development of CAD. In quantitative RT-qPCR, the expression of PRKX, HSPH1 and ZMYND11 was down-regulated and consistent with the integrated analysis. Conclusions: Identified 7 differentially expressed genes (IL1R2, ARRB2, PRKX, PLAUR, HSPH1, ZMYND11 and S100A8) may play crucial roles in the development of CAD

Transcriptional changes in trichothiodystrophy cells

Mutations in three of the genes encoding the XPB, XPD and TTDA components of transcription factor TFIIH can result in the clinical phenotype of trichothiodystrophy (TTD). Different mutations in XPB and XPD can instead cause xeroderma pigmentosum (XP). The completely different features of these disorders have been attributed to TTD being a transcription syndrome. In order to detect transcriptional differences between TTD and XP cells from the XP-D complementation group, we have compared gene expression profiles in cultured fibroblasts from normal, XP and TTD donors. Although we detected transcriptional differences between individual cell strains, using an algorithm of moderate stringency, we did not identify any genes whose expression was reproducibly different in proliferating fibroblasts from each type of donor. Following UV-irradiation, many genes were up- and down-regulated in all three cell types. The microarray analysis indicated some apparent differences between the different donor types, but on more detailed inspection, these turned out to be false positives. We conclude that there are minimal differences in gene expression in proliferating fibroblasts from TTD, XP-D and normal donors

Identification of transcriptional regulatory networks specific to pilocytic astrocytoma.

BackgroundPilocytic Astrocytomas (PAs) are common low-grade central nervous system malignancies for which few recurrent and specific genetic alterations have been identified. In an effort to better understand the molecular biology underlying the pathogenesis of these pediatric brain tumors, we performed higher-order transcriptional network analysis of a large gene expression dataset to identify gene regulatory pathways that are specific to this tumor type, relative to other, more aggressive glial or histologically distinct brain tumours.MethodsRNA derived from frozen human PA tumours was subjected to microarray-based gene expression profiling, using Affymetrix U133Plus2 GeneChip microarrays. This data set was compared to similar data sets previously generated from non-malignant human brain tissue and other brain tumour types, after appropriate normalization.ResultsIn this study, we examined gene expression in 66 PA tumors compared to 15 non-malignant cortical brain tissues, and identified 792 genes that demonstrated consistent differential expression between independent sets of PA and non-malignant specimens. From this entire 792 gene set, we used the previously described PAP tool to assemble a core transcriptional regulatory network composed of 6 transcription factor genes (TFs) and 24 target genes, for a total of 55 interactions. A similar analysis of oligodendroglioma and glioblastoma multiforme (GBM) gene expression data sets identified distinct, but overlapping, networks. Most importantly, comparison of each of the brain tumor type-specific networks revealed a network unique to PA that included repressed expression of ONECUT2, a gene frequently methylated in other tumor types, and 13 other uniquely predicted TF-gene interactions.ConclusionsThese results suggest specific transcriptional pathways that may operate to create the unique molecular phenotype of PA and thus opportunities for corresponding targeted therapeutic intervention. Moreover, this study also demonstrates how integration of gene expression data with TF-gene and TF-TF interaction data is a powerful approach to generating testable hypotheses to better understand cell-type specific genetic programs relevant to cancer

Identifying mRNA targets of microRNA dysregulated in cancer: with application to clear cell Renal Cell Carcinoma

BACKGROUND. MicroRNA regulate mRNA levels in a tissue specific way, either by inducing degradation of the transcript or by inhibiting translation or transcription. Putative mRNA targets of microRNA identified from seed sequence matches are available in many databases. However, such matches have a high false positive rate and cannot identify tissue specificity of regulation. RESULTS. We describe a simple method to identify direct mRNA targets of microRNA dysregulated in cancers from expression level measurements in patient matched tumor/normal samples. The word "direct" is used here in a strict sense to: a) represent mRNA which have an exact seed sequence match to the microRNA in their 3'UTR, b) the seed sequence match is strictly conserved across mouse, human, rat and dog genomes, c) the mRNA and microRNA expression levels can distinguish tumor from normal with high significance and d) the microRNA/mRNA expression levels are strongly and significantly anti-correlated in tumor and/or normal samples. We apply and validate the method using clear cell Renal Cell Carcinoma (ccRCC) and matched normal kidney samples, limiting our analysis to mRNA targets which undergo degradation of the mRNA transcript because of a perfect seed sequence match. Dysregulated microRNA and mRNA are first identified by comparing their expression levels in tumor vs normal samples. Putative dysregulated microRNA/mRNA pairs are identified from these using seed sequence matches, requiring that the seed sequence be conserved in human/dog/rat/mouse genomes. These are further pruned by requiring a strong anti-correlation signature in tumor and/or normal samples. The method revealed many new regulations in ccRCC. For instance, loss of miR-149, miR-200c and mir-141 causes gain of function of oncogenes (KCNMA1, LOX), VEGFA and SEMA6A respectively and increased levels of miR-142-3p, miR-185, mir-34a, miR-224, miR-21 cause loss of function of tumor suppressors LRRC2, PTPN13, SFRP1, ERBB4, and (SLC12A1, TCF21) respectively. We also found strong anti-correlation between VEGFA and the miR-200 family of microRNA: miR-200a*, 200b, 200c and miR-141. Several identified microRNA/mRNA pairs were validated on an independent set of matched ccRCC/normal samples. The regulation of SEMA6A by miR-141 was verified by a transfection assay. CONCLUSIONS. We describe a simple and reliable method to identify direct gene targets of microRNA in any cancer. The constraints we impose (strong dysregulation signature for microRNA and mRNA levels between tumor/normal samples, evolutionary conservation of seed sequence and strong anti-correlation of expression levels) remove spurious matches and identify a subset of robust, tissue specific, functional mRNA targets of dysregulated microRNA.Cancer Institute of New Jersy; New Jersey Commission for Cacner Research; Lineberger Comprehensive Cancer Center Tissue Procurement and Genomics Core Facility; Crawford Fun

Comparing gene discovery from Affymetrix GeneChip microarrays and Clontech PCR-select cDNA subtraction: a case study

BACKGROUND: Several high throughput technologies have been employed to identify differentially regulated genes that may be molecular targets for drug discovery. Here we compared the sets of differentially regulated genes discovered using two experimental approaches: a subtracted suppressive hybridization (SSH) cDNA library methodology and Affymetrix GeneChip(® )technology. In this "case study" we explored the transcriptional pattern changes during the in vitro differentiation of human monocytes to myeloid dendritic cells (DC), and evaluated the potential for novel gene discovery using the SSH methodology. RESULTS: The same RNA samples isolated from peripheral blood monocyte precursors and immature DC (iDC) were used for GeneChip microarray probing and SSH cDNA library construction. 10,000 clones from each of the two-way SSH libraries (iDC-monocytes and monocytes-iDC) were picked for sequencing. About 2000 transcripts were identified for each library from 8000 successful sequences. Only 70% to 75% of these transcripts were represented on the U95 series GeneChip microarrays, implying that 25% to 30% of these transcripts might not have been identified in a study based only on GeneChip microarrays. In addition, about 10% of these transcripts appeared to be "novel", although these have not yet been closely examined. Among the transcripts that are also represented on the chips, about a third were concordantly discovered as differentially regulated between iDC and monocytes by GeneChip microarray transcript profiling. The remaining two thirds were either not inferred as differentially regulated from GeneChip microarray data, or were called differentially regulated but in the opposite direction. This underscores the importance both of generating reciprocal pairs of SSH libraries, and of real-time RT-PCR confirmation of the results. CONCLUSIONS: This study suggests that SSH could be used as an alternative and complementary transcript profiling tool to GeneChip microarrays, especially in identifying novel genes and transcripts of low abundance

The Reproducibility of Lists of Differentially Expressed Genes in Microarray Studies

Reproducibility is a fundamental requirement in scientific experiments and clinical contexts. Recent publications raise concerns about the reliability of microarray technology because of the apparent lack of agreement between lists of differentially expressed genes (DEGs). In this study we demonstrate that (1) such discordance may stem from ranking and selecting DEGs solely by statistical significance (P) derived from widely used simple t-tests; (2) when fold change (FC) is used as the ranking criterion, the lists become much more reproducible, especially when fewer genes are selected; and (3) the instability of short DEG lists based on P cutoffs is an expected mathematical consequence of the high variability of the t-values. We recommend the use of FC ranking plus a non-stringent P cutoff as a baseline practice in order to generate more reproducible DEG lists. The FC criterion enhances reproducibility while the P criterion balances sensitivity and specificity

BNP-Seq: Bayesian Nonparametric Differential Expression Analysis of Sequencing Count Data

We perform differential expression analysis of high-throughput sequencing count data under a Bayesian nonparametric framework, removing sophisticated ad-hoc pre-processing steps commonly required in existing algorithms. We propose to use the gamma (beta) negative binomial process, which takes into account different sequencing depths using sample-specific negative binomial probability (dispersion) parameters, to detect differentially expressed genes by comparing the posterior distributions of gene-specific negative binomial dispersion (probability) parameters. These model parameters are inferred by borrowing statistical strength across both the genes and samples. Extensive experiments on both simulated and real-world RNA sequencing count data show that the proposed differential expression analysis algorithms clearly outperform previously proposed ones in terms of the areas under both the receiver operating characteristic and precision-recall curves.Comment: To appear in Journal of the American Statistical Associatio

The Functional Consequences of Variation in Transcription Factor Binding

One goal of human genetics is to understand how the information for precise and dynamic gene expression programs is encoded in the genome. The interactions of transcription factors (TFs) with DNA regulatory elements clearly play an important role in determining gene expression outputs, yet the regulatory logic underlying functional transcription factor binding is poorly understood. Many studies have focused on characterizing the genomic locations of TF binding, yet it is unclear to what extent TF binding at any specific locus has functional consequences with respect to gene expression output. To evaluate the context of functional TF binding we knocked down 59 TFs and chromatin modifiers in one HapMap lymphoblastoid cell line. We then identified genes whose expression was affected by the knockdowns. We intersected the gene expression data with transcription factor binding data (based on ChIP-seq and DNase-seq) within 10 kb of the transcription start sites of expressed genes. This combination of data allowed us to infer functional TF binding. On average, 14.7% of genes bound by a factor were differentially expressed following the knockdown of that factor, suggesting that most interactions between TF and chromatin do not result in measurable changes in gene expression levels of putative target genes. We found that functional TF binding is enriched in regulatory elements that harbor a large number of TF binding sites, at sites with predicted higher binding affinity, and at sites that are enriched in genomic regions annotated as active enhancers.Comment: 30 pages, 6 figures (7 supplemental figures and 6 supplemental tables available upon request to [email protected]). Submitted to PLoS Genetic