Application of Volcano Plots in Analyses of mRNA Differential Expressions with Microarrays

Volcano plot displays unstandardized signal (e.g. log-fold-change) against noise-adjusted/standardized signal (e.g. t-statistic or -log10(p-value) from the t test). We review the basic and an interactive use of the volcano plot, and its crucial role in understanding the regularized t-statistic. The joint filtering gene selection criterion based on regularized statistics has a curved discriminant line in the volcano plot, as compared to the two perpendicular lines for the "double filtering" criterion. This review attempts to provide an unifying framework for discussions on alternative measures of differential expression, improved methods for estimating variance, and visual display of a microarray analysis result. We also discuss the possibility to apply volcano plots to other fields beyond microarray.Comment: 8 figure

Differential expression of skeletal muscle genes following administration of clenbuterol to exercised horses.

BackgroundClenbuterol, a beta2-adrenergic receptor agonist, is used therapeutically to treat respiratory conditions in the horse. However, by virtue of its mechanism of action it has been suggested that clenbuterol may also have repartitioning affects in horses and as such the potential to affect performance. Clenbuterol decreases the percent fat and increases fat-free mass following high dose administration in combination with intense exercise in horses. In the current study, microarray analysis and real-time PCR were used to study the temporal effects of low and high dose chronic clenbuterol administration on differential gene expression of several skeletal muscle myosin heavy chains, genes involved in lipid metabolism and the β2-adrenergic receptor. The effect of clenbuterol administration on differential gene expression has not been previously reported in the horse, therefore the primary objective of the current study was to describe clenbuterol-induced temporal changes in gene expression following chronic oral administration of clenbuterol at both high and low doses.ResultsSteady state clenbuterol concentrations were achieved at approximately 50 h post administration of the first dose for the low dose regimen and at approximately 18-19 days (10 days post administration of 3.2 μg/kg) for the escalating dosing regimen. Following chronic administration of the low dose (0.8 μg/kg BID) of clenbuterol, a total of 114 genes were differentially expressed, however, none of these changes were found to be significant following FDR adjustment of the p-values. A total of 7,093 genes were differentially expressed with 3,623 genes up regulated and 3,470 genes down regulated following chronic high dose administration. Of the genes selected for further study by real-time PCR, down-regulation of genes encoding myosin heavy chains 2 and 7, steroyl CoA desaturase and the β2-adrenergic receptor were noted. For most genes, expression levels returned towards baseline levels following cessation of drug administration.ConclusionThis study showed no evidence of modified gene expression following chronic low dose administration of clenbuterol to horses. However, following chronic administration of high doses of clenbuterol alterations were noted in transcripts encoding various myosin heavy chains, lipid metabolizing enzymes and the β2-adrenergic receptor

Transcriptomic effects of the non-steroidal anti-inflammatory drug Ibuprofen in the marine bivalve Mytilus galloprovincialis Lam

The transcriptomic effects of Ibuprofen (IBU) in the digestive gland tissue of Mytilus galloprovincialis Lam. specimens exposed at low environmental concentrations (250 ng L-1) are presented. Using a 1.7 K feature cDNA microarray along with linear models and empirical Bayes statistical methods 225 differentially expressed genes were identified in mussels treated with IBU across a 15-day period. Transcriptional dynamics were typical of an adaptive response with a peak of gene expression change at day 7 (177 features, representing about 11% of sequences available for analysis) and an almost full recovery at the end of the exposure period. Functional genomics by means of Gene Ontology term analysis unraveled typical mussel stress responses i.e. aminoglycan (chitin) metabolic processes but also more specific effects such as the regulation of NF-kappa B transcription factor activity. (C) 2016 Elsevier Ltd. All rights reserved

Differential gene expression graphs: A data structure for classification in DNA microarrays

This paper proposes an innovative data structure to be used as a backbone in designing microarray phenotype sample classifiers. The data structure is based on graphs and it is built from a differential analysis of the expression levels of healthy and diseased tissue samples in a microarray dataset. The proposed data structure is built in such a way that, by construction, it shows a number of properties that are perfectly suited to address several problems like feature extraction, clustering, and classificatio

Cytological and transcript analyses reveal fat and lazy persister-like bacilli in tuberculous sputum

As nonreplicating tubercle bacilli are tolerant to the cidal action of antibiotics and resistant to multiple stresses, identification of this persister-like population of tubercle bacilli in sputum presents exciting and tractable new opportunities to investigate both responses to chemotherapy and the transmission of tuberculosis

Whole-transcriptome, high-throughput RNA sequence analysis of the bovine macrophage response to Mycobacterium bovis infection in vitro

BACKGROUND: Mycobacterium bovis, the causative agent of bovine tuberculosis, is an intracellular pathogen that can persist inside host macrophages during infection via a diverse range of mechanisms that subvert the host immune response. In the current study, we have analysed and compared the transcriptomes of M. bovis-infected monocyte-derived macrophages (MDM) purified from six Holstein-Friesian females with the transcriptomes of non-infected control MDM from the same animals over a 24 h period using strand-specific RNA sequencing (RNA-seq). In addition, we compare gene expression profiles generated using RNA-seq with those previously generated by us using the high-density Affymetrix® GeneChip® Bovine Genome Array platform from the same MDM-extracted RNA. RESULTS: A mean of 7.2 million reads from each MDM sample mapped uniquely and unambiguously to single Bos taurus reference genome locations. Analysis of these mapped reads showed 2,584 genes (1,392 upregulated; 1,192 downregulated) and 757 putative natural antisense transcripts (558 upregulated; 119 downregulated) that were differentially expressed based on sense and antisense strand data, respectively (adjusted P-value ≤ 0.05). Of the differentially expressed genes, 694 were common to both the sense and antisense data sets, with the direction of expression (i.e. up- or downregulation) positively correlated for 693 genes and negatively correlated for the remaining gene. Gene ontology analysis of the differentially expressed genes revealed an enrichment of immune, apoptotic and cell signalling genes. Notably, the number of differentially expressed genes identified from RNA-seq sense strand analysis was greater than the number of differentially expressed genes detected from microarray analysis (2,584 genes versus 2,015 genes). Furthermore, our data reveal a greater dynamic range in the detection and quantification of gene transcripts for RNA-seq compared to microarray technology. CONCLUSIONS: This study highlights the value of RNA-seq in identifying novel immunomodulatory mechanisms that underlie host-mycobacterial pathogen interactions during infection, including possible complex post-transcriptional regulation of host gene expression involving antisense RNA

Do-it-yourself: construction of a custom cDNA macroarray platform with high sensitivity and linear range

Background: Research involving gene expression profiling and clinical applications, such as diagnostics and prognostics, often require a DNA array platform that is flexibly customisable and cost-effective, but at the same time is highly sensitive and capable of accurately and reproducibly quantifying the transcriptional expression of a vast number of genes over the whole transcriptome dynamic range using low amounts of RNA sample. Hereto, a set of easy-to-implement practical optimisations to the design of cDNA-based nylon macroarrays as well as sample (33)P-labeling, hybridisation protocols and phosphor screen image processing were analysed for macroarray performance. Results: The here proposed custom macroarray platform had an absolute sensitivity as low as 50,000 transcripts and a linear range of over 5 log-orders. Its quality of identifying differentially expressed genes was at least comparable to commercially available microchips. Interestingly, the quantitative accuracy was found to correlate significantly with corresponding reversed transcriptase - quantitative PCR values, the gold standard gene expression measure (Pearson's correlation test p < 0.0001). Furthermore, the assay has low cost and input RNA requirements (0.5 mu g and less) and has a sound reproducibility. Conclusions: Results presented here, demonstrate for the first time that self-made cDNA-based nylon macroarrays can produce highly reliable gene expression data with high sensitivity and covering the entire mammalian dynamic range of mRNA abundances. Starting off from minimal amounts of unamplified total RNA per sample, a reasonable amount of samples can be assayed simultaneously for the quantitative expression of hundreds of genes in an easily customisable and cost-effective manner

Association Analysis Techniques for Discovering Functional Modules from Microarray Data

An application of great interest in microarray data analysis is the identification of a group of genes that show very similar patterns of expression in a data set, and are expected to represent groups of genes that perform common/similar functions, also known as functional modules. Although clustering offers a natural solution to this problem, it suffers from the limitation that it uses all the conditions to compare two genes, whereas only a subset of them may be relevant. Association analysis offers an alternative route for finding such groups of genes that may be co-expressed only over a subset of the experimental conditions used to prepare the data set. The techniques in this field attempt to find groups of data objects that contain coherent values across a set of attributes, in an exhaustive and efficient manner. In this paper, we illustrate how a generalization of the techniques in association analysis for real-valued data can be utilized to extract coherent functional modules from large microarray data sets