Are gene expression microarray analyses reliable? A review of studies of retinoic acid responsive genes.
Ontology highlight
ABSTRACT:
SUBMITTER: van der Spek PJ
PROVIDER: S-EPMC5171945 | biostudies-literature | 2003 Feb
REPOSITORIES: biostudies-literature
Publications
Are gene expression microarray analyses reliable? A review of studies of retinoic acid responsive genes.
van der Spek Peter J PJ Kremer Andreas A Murry Lynn L Walker Michael G MG
Genomics, proteomics & bioinformatics 20030201 1
Microarray analyses of gene expression are widely used, but reports of the same analyses by different groups give widely divergent results, and raise questions regarding reproducibility and reliability. We take as an example recent published reports on microarray experiments that were designed to identify retinoic acid responsive genes. These reports show substantial differences in their results. In this article, we review the methodology, results, and potential causes of differences in these ap ...[more]
Similar Datasets
Project description:The mechanism of action of retinoic acid (RA) is of broad relevance to cell and developmental biology, nutrition, and cancer chemotherapy. RA is known to induce expression of the Burkitt's lymphoma receptor 1 (BLR1) gene which propels RA-induced cell cycle arrest and differentiation of HL-60 human myeloblastic leukemia cells, motivating the present analysis of transcriptional regulation of blr1 expression by RA. The RA-treated HL-60 cells used here expressed all RA receptor (RAR) and retinoid X receptor (RXR) subtypes (as detected by Northern analysis) except RXRgamma. Treatment with RAR- and RXR-selective ligands showed that RARalpha synergized with RXRalpha to transcriptionally activate blr1 expression. A 5'-flanking region capable of supporting RA-induced blr1 activation in HL-60 cells was found to contain a 205-bp sequence in the distal portion that was necessary for transcriptional activation by RA. Within this sequence DNase I footprinting revealed that RA induced binding of a nuclear protein complex to an element containing two GT boxes. Electromobility shift assays (EMSAs) and supershift assays showed that this element bound recombinant RARalpha and RXRalpha. Without RA there was neither complex binding nor transcriptional activation. Both GT boxes were needed for binding the complex, and mutation of either GT box caused the loss of transcriptional activation by RA. The ability of this cis-acting RAR-RXR binding element to activate transcription in response to RA also depended on downstream sequences where an octamer transcription factor 1 (Oct1) site and a nuclear factor of activated T cells (NFATc) site between this element and the transcriptional start, as well as a cyclic AMP response element binding factor (CREB) site between the transcriptional start and first exon of the blr1 gene, were necessary. Each of these sites bound its corresponding transcription factor. A transcription factor-transcription factor binding array analysis of nuclear lysate from RA-treated cells indicated several prominent RARalpha binding partners; among these, Oct1, NFATc3, and CREB2 were identified by competition EMSA and supershift and chromatin immunoprecipitation assays as components of the complex. RA upregulated expression of these three factors. In sum the results of the present study indicate that RA-induced expression of blr1 expression depends on a novel RA response element. This cis-acting element approximately 1 kb upstream of the transcriptional start consists of two GT boxes that bind RAR and RXR in a nuclear protein complex that also contains Oct1, NFATc3, and CREB2 bound to their cognate downstream consensus binding sites.
| S-EPMC355834 | biostudies-literature
Project description:limma is an R/Bioconductor software package that provides an integrated solution for analysing data from gene expression experiments. It contains rich features for handling complex experimental designs and for information borrowing to overcome the problem of small sample sizes. Over the past decade, limma has been a popular choice for gene discovery through differential expression analyses of microarray and high-throughput PCR data. The package contains particularly strong facilities for reading, normalizing and exploring such data. Recently, the capabilities of limma have been significantly expanded in two important directions. First, the package can now perform both differential expression and differential splicing analyses of RNA sequencing (RNA-seq) data. All the downstream analysis tools previously restricted to microarray data are now available for RNA-seq as well. These capabilities allow users to analyse both RNA-seq and microarray data with very similar pipelines. Second, the package is now able to go past the traditional gene-wise expression analyses in a variety of ways, analysing expression profiles in terms of co-regulated sets of genes or in terms of higher-order expression signatures. This provides enhanced possibilities for biological interpretation of gene expression differences. This article reviews the philosophy and design of the limma package, summarizing both new and historical features, with an emphasis on recent enhancements and features that have not been previously described.
| S-EPMC4402510 | biostudies-literature
Project description:Retinoic acid (RA) signaling regulates a variety of developmental processes through controlling the expression of numerous genes. Here, we have identified and characterized RA-responsive genes in mouse kidney development. Analysis of isolated embryonic kidneys cultured in the presence and absence of RA identified 33 candidates of RA-responsive genes. Most of these candidate genes were down-regulated by treatment with the RA receptor antagonist. Many of them have potential binding sites for Elf5, one of the RA-responsive genes, in their promoter region. Whole-mount in situ hybridization demonstrated that specific expression of Elf5 in the ureteric trunk depends on RA. RA-dependent expression in the ureteric trunk was also demonstrated for the sodium channel subunit Scnn1b, which has been shown to be the marker gene of the collecting duct. In contrast, the expression of Ecm1, Tnfsf13b and IL-33 was detected in the stromal mesenchymal cells. Both Tnfsf13b and IL-33 were previously shown to cause NF-κB activation. We have demonstrated that the inhibition of NF-κB signaling with specific inhibitors suppresses branching morphogenesis of the ureteric bud. Our study thus identifies and characterizes RA-dependent upregulated genes in kidney development, and suggests an involvement of NF-κB signaling in the branching morphogenesis.
2014-11-14 | GSE57656 | GEO
Project description:Retinoic acid (RA) signaling regulates a variety of developmental processes through controlling the expression of numerous genes. Here, we have identified and characterized RA-responsive genes in mouse kidney development. Analysis of isolated embryonic kidneys cultured in the presence and absence of RA identified 33 candidates of RA-responsive genes. Most of these candidate genes were down-regulated by treatment with the RA receptor antagonist. Many of them have potential binding sites for Elf5, one of the RA-responsive genes, in their promoter region. Whole-mount in situ hybridization demonstrated that specific expression of Elf5 in the ureteric trunk depends on RA. RA-dependent expression in the ureteric trunk was also demonstrated for the sodium channel subunit Scnn1b, which has been shown to be the marker gene of the collecting duct. In contrast, the expression of Ecm1, Tnfsf13b and IL-33 was detected in the stromal mesenchymal cells. Both Tnfsf13b and IL-33 were previously shown to cause NF-?B activation. We have demonstrated that the inhibition of NF-?B signaling with specific inhibitors suppresses branching morphogenesis of the ureteric bud. Our study thus identifies and characterizes RA-dependent upregulated genes in kidney development, and suggests an involvement of NF-?B signaling in the branching morphogenesis. To obtain genome-wide gene expression profiles of kidneys cultured in 10% FBS DMEM, serum-free DMEM, and serum-free DMEM with RA, we performed microarray analysis on RNA prepared from these three samples. Total RNA was prepared using the RNeasy Mini Kit (QIAGEN) according to the manufacturer's instruction. Synthesis of cDNA, in vitro transcription and biotin labeling cRNA, and hybridization to GeneChip Mouse Gene 1.0 ST array (Affymetrix) were performed according to Affymetrix protocols.
2014-11-14 | E-GEOD-57656 | biostudies-arrayexpress
Project description:The vertebrate nuclear receptor subfamily 2, group f (nr2f) genes encode orphan receptors that have the capacity to act as negative regulators of retinoic acid (RA) signaling.We describe embryonic and larval expression of four of the six zebrafish nr2f genes, nr2f1a, nr2f1b, nr2f2, and nr2f5. These genes show highly regulated patterns of expression within the central nervous system, including in the developing hindbrain, as well as in the mesoderm and endoderm. We also investigated the role of RA and fibroblast growth factor (Fgf) signaling in regulating early nr2f gene expression. RA is not required for nr2f expression in the hindbrain; however, exogenous RA can repress this expression. Conversely, we find that RA positively regulates nr2f1a expression in trunk endoderm and mesoderm. Fgf signaling is not required for nr2f expression onset in the hindbrain; however, it may play a role in maintaining rhombomere-specific expression.We report detailed expression analysis of four nr2f genes in all three germ layers. The onset of nr2f expression in the hindbrain does not require RA or Fgf signals. Our finding that RA positively regulates nr2f1a expression in the trunk supports the possibility that Nr2fs function in a negative feedback loop to modulate RA signaling in this region.
| S-EPMC3459307 | biostudies-literature
Project description:The quantitative real-time polymerase chain reaction (qRT-PCR) approach has become a widely used method to analyze expression patterns of target genes. The selection of an optimal reference gene is a prerequisite for the accurate normalization of gene expression in qRT-PCR. The present study constitutes the first systematic evaluation of potential reference genes in Rhododendron molle G. Don. Eleven candidate reference genes in different tissues and flowers at different developmental stages of R. molle were assessed using the following three software packages: GeNorm, NormFinder, and BestKeeper. The results showed that EF1-? (elongation factor 1-alpha), 18S (18s ribosomal RNA), and RPL3 (ribosomal protein L3) were the most stable reference genes in developing rhododendron flowers and, thus, in all of the tested samples, while tublin (TUB) was the least stable. ACT5 (actin), RPL3, 18S, and EF1-? were found to be the top four choices for different tissues, whereas TUB was not found to favor qRT-PCR normalization in these tissues. Three stable reference genes are recommended for the normalization of qRT-PCR data in R. molle. Furthermore, the expression profiles of RmPSY (phytoene synthase) and RmPDS (phytoene dehydrogenase) were assessed using EF1-?, 18S, ACT5, RPL3, and their combination as internals. Similar trends were found, but these trends varied when the least stable reference gene TUB was used. The results further prove that it is necessary to validate the stability of reference genes prior to their use for normalization under different experimental conditions. This study provides useful information for reliable qRT-PCR data normalization in gene studies of R. molle.
| S-EPMC5067439 | biostudies-literature
Project description:A common observation in the analysis of gene expression data is that many genes display similarity in their expression patterns and therefore appear to be co-regulated. However, the variation associated with microarray data and the complexity of the experimental designs make the acquisition of co-expressed genes a challenge. We developed a novel method for Extracting microarray gene expression Patterns and Identifying co-expressed Genes, designated as EPIG. The approach utilizes the underlying structure of gene expression data to extract patterns and identify co-expressed genes that are responsive to experimental conditions.Through evaluation of the correlations among profiles, the magnitude of variation in gene expression profiles, and profile signal-to-noise ratio's, EPIG extracts a set of patterns representing co-expressed genes. The method is shown to work well with a simulated data set and microarray data obtained from time-series studies of dauer recovery and L1 starvation in C. elegans and after ultraviolet (UV) or ionizing radiation (IR)-induced DNA damage in diploid human fibroblasts. With the simulated data set, EPIG extracted the appropriate number of patterns which were more stable and homogeneous than the set of patterns that were determined using the CLICK or CAST clustering algorithms. However, CLICK performed better than EPIG and CAST with respect to the average correlation between clusters/patterns of the simulated data. With real biological data, EPIG extracted more dauer-specific patterns than CLICK. Furthermore, analysis of the IR/UV data revealed 18 unique patterns and 2661 genes out of approximately 17,000 that were identified as significantly expressed and categorized to the patterns by EPIG. The time-dependent patterns displayed similar and dissimilar responses between IR and UV treatments. Gene Ontology analysis applied to each pattern-related subset of co-expressed genes revealed underlying biological processes affected by IR- and/or UV- induced DNA damage.EPIG competed with CLICK and performed better than CAST in extracting patterns from simulated data. EPIG extracted more biological informative patterns and co-expressed genes from both C. elegans and IR/UV-treated human fibroblasts. Using Gene Ontology analysis of the genes in the patterns extracted by EPIG, several key biological categories related to p53-dependent cell cycle control were revealed from the IR/UV data. Among them were mitotic cell cycle, DNA replication, DNA repair, cell cycle checkpoint, and G0-like status transition. EPIG can be applied to data sets from a variety of experimental designs.
| S-EPMC2194742 | biostudies-other
Project description:Accurately predicting and testing the types of Pulmonary arterial hypertension (PAH) of each patient using cost-effective microarray-based expression data and machine learning algorithms could greatly help either identifying the most targeting medicine or adopting other therapeutic measures that could correct/restore defective genetic signaling at the early stage. Furthermore, the prediction model construction processes can also help identifying highly informative genes controlling PAH, leading to enhanced understanding of the disease etiology and molecular pathways. In this study, we used several different gene filtering methods based on microarray expression data obtained from a high-quality patient PAH dataset. Following that, we proposed a novel feature selection and refinement algorithm in conjunction with well-known machine learning methods to identify a small set of highly informative genes. Results indicated that clusters of small-expression genes could be extremely informative at predicting and differentiating different forms of PAH. Additionally, our proposed novel feature refinement algorithm could lead to significant enhancement in model performance. To summarize, integrated with state-of-the-art machine learning and novel feature refining algorithms, the most accurate models could provide near-perfect classification accuracies using very few (close to ten) low-expression genes.
| S-EPMC6705875 | biostudies-literature
Project description:BACKGROUND: RT-qPCR is a preferred method for rapid and reliable quantification of gene expression studies. Appropriate application of RT-qPCR in such studies requires the use of reference gene(s) as an internal control to normalize mRNA levels between different samples for an exact comparison of gene expression level. However, recent studies have shown that no single reference gene is universal for all experiments. Thus, the identification of high quality reference gene(s) is of paramount importance for the interpretation of data generated by RT-qPCR. Only a few studies on reference genes have been done in plants and none in peach (Prunus persica L. Batsch). Therefore, the present study was conducted to identify suitable reference gene(s) for normalization of gene expression in peach. RESULTS: In this work, eleven reference genes were investigated in different peach samples using RT-qPCR with SYBR green. These genes are: actin 2/7 (ACT), cyclophilin (CYP2), RNA polymerase II (RP II), phospholipase A2 (PLA2), ribosomal protein L13 (RPL13), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), 18S ribosomal RNA (18S rRNA), tubblin beta (TUB), tubblin alpha (TUA), translation elongation factor 2 (TEF2) and ubiquitin 10 (UBQ10). All eleven reference genes displayed a wide range of Cq values in all samples, indicating that they expressed variably. The stability of these genes except for RPL13 was determined by three different descriptive statistics, geNorm, NormFinder and BestKeeper, which produced highly comparable results. CONCLUSION: Our study demonstrates that expression stability varied greatly between genes studied in peach. Based on the results from geNorm, NormFinder and BestKeeper analyses, for all the sample pools analyzed, TEF2, UBQ10 and RP II were found to be the most suitable reference genes with a very high statistical reliability, and TEF2 and RP II for the other sample series, while 18S rRNA, RPL13 and PLA2 were unsuitable as internal controls. GAPDH and ACT also performed poorly and were less stable in our analysis. To achieve accurate comparison of levels of gene expression, two or more reference genes must be used for data normalization. The combinations of TEF2/UBQ10/RP II and TEF2/RP II were suggested for use in all samples and subsets, respectively.
| S-EPMC3224724 | biostudies-literature
Project description:Neural tube defects (NTDs) occur in as many as 0.5-2 per 1000 live births in the United States. One of the most common and severe neural tube defects is meningomyelocele (MM) resulting from failed closure of the caudal end of the neural tube. MM has been induced by retinoic acid teratogenicity in rodent models. We hypothesized that genetic variants influencing retinoic acid (RA) induction via retinoic acid receptors (RARs) may be associated with risk for MM.We analyzed 47 single nucleotide polymorphisms (SNPs) that span across the three retinoic acid receptor genes using the SNPlex genotyping platform. Our cohort consisted of 610 MM families.One variant in the RARA gene (rs12051734), three variants in the RARB gene (rs6799734, rs12630816, rs17016462), and a single variant in the RARG gene (rs3741434) were found to be statistically significant at p < 0.05.RAR genes were associated with risk for MM. For all associated SNPs, the rare allele conferred a protective effect for MM susceptibility.
| S-EPMC3044482 | biostudies-literature