Pooling samples within microarray studies: a comparative analysis of rat liver transcription response to prototypical toxicants

Systematic transcriptome-wide meta-analysis across endocrine disrupting chemicals reveals shared and unique liver pathways, gene networks, and disease associations

Cardiometabolic disorders (CMD) are a growing public health problem across the world. Among the known cardiometabolic risk factors are compounds that induce endocrine and metabolic dysfunctions, such as endocrine disrupting chemicals (EDCs). To date, how EDCs influence molecular programs and cardiometabolic risks has yet to be fully elucidated, especially considering the complexity contributed by species-, chemical-, and dose-specific effects. Moreover, different experimental and analytical methodologies employed by different studies pose challenges when comparing findings across studies. To explore the molecular mechanisms of EDCs in a systematic manner, we established a data-driven computational approach to meta-analyze 30 human, mouse, and rat liver transcriptomic datasets for 4 EDCs, namely bisphenol A (BPA), bis(2-ethylhexyl) phthalate (DEHP), tributyltin (TBT), and perfluorooctanoic acid (PFOA). Our computational pipeline uniformly re-analyzed pre-processed quality-controlled microarray data and raw RNAseq data, derived differentially expressed genes (DEGs) and biological pathways, modeled gene regulatory networks and regulators, and determined CMD associations based on gene overlap analysis. Our approach revealed that DEHP and PFOA shared stable transcriptomic signatures that are enriched for genes associated with CMDs, suggesting similar mechanisms of action such as perturbations of peroxisome proliferator-activated receptor gamma (PPARγ) signaling and liver gene network regulators VNN1 and ACOT2. In contrast, TBT exhibited highly divergent gene signatures, pathways, network regulators, and disease associations from the other EDCs. In addition, we found that the rat, mouse, and human BPA studies showed highly variable transcriptomic patterns, providing molecular support for the variability in BPA responses. Our work offers insights into the commonality and differences in the molecular mechanisms of various EDCs and establishes a streamlined data-driven workflow to compare molecular mechanisms of environmental substances to elucidate the underlying connections between chemical exposure and disease risks.

Overnight Corticosterone and Gene Expression in Mouse Hippocampus: Time Course during Resting Period

The aim of the experiment was to test the effect of an elevated level of glucocorticoids on the mouse hippocampal transcriptome after 12 h of treatment with corticosterone that was administered during an active phase of the circadian cycle. Additionally, we also tested the circadian changes in gene expression and the decay time of transcriptomic response to corticosterone. Gene expression was analyzed using microarrays. Obtained results show that transcriptomic responses to glucocorticoids are heterogeneous in terms of the decay time with some genes displaying persistent changes in expression during 9 h of rest. We have also found a considerable overlap between genes regulated by corticosterone and genes implicated previously in stress response. The examples of such genes are Acer2, Agt, Apod, Aqp4, Etnppl, Fabp7, Fam107a, Fjx1, Fmo2, Galnt15, Gjc2, Heph, Hes5, Htra1, Jdp2, Kif5a, Lfng, Lrg1, Mgp, Mt1, Pglyrp1, Pla2g3, Plin4, Pllp, Ptgds, Ptn, Slc2a1, Slco1c1, Sult1a1, Thbd and Txnip. This indicates that the applied model is a useful tool for the investigation of mechanisms underlying the stress response.

Altered hippocampal gene expression, glial cell population, and neuronal excitability in aminopeptidase P1 deficiency

Inborn errors of metabolism are often associated with neurodevelopmental disorders and brain injury. A deficiency of aminopeptidase P1, a proline-specific endopeptidase encoded by the Xpnpep1 gene, causes neurological complications in both humans and mice. In addition, aminopeptidase P1-deficient mice exhibit hippocampal neurodegeneration and impaired hippocampus-dependent learning and memory. However, the molecular and cellular changes associated with hippocampal pathology in aminopeptidase P1 deficiency are unclear. We show here that a deficiency of aminopeptidase P1 modifies the glial population and neuronal excitability in the hippocampus. Microarray and real-time quantitative reverse transcription-polymerase chain reaction analyses identified 14 differentially expressed genes (Casp1, Ccnd1, Myoc, Opalin, Aldh1a2, Aspa, Spp1, Gstm6, Serpinb1a, Pdlim1, Dsp, Tnfaip6, Slc6a20a, Slc22a2) in the Xpnpep1^−/− hippocampus. In the hippocampus, aminopeptidase P1-expression signals were mainly detected in neurons. However, deficiency of aminopeptidase P1 resulted in fewer hippocampal astrocytes and increased density of microglia in the hippocampal CA3 area. In addition, Xpnpep1^−/− CA3b pyramidal neurons were more excitable than wild-type neurons. These results indicate that insufficient astrocytic neuroprotection and enhanced neuronal excitability may underlie neurodegeneration and hippocampal dysfunction in aminopeptidase P1 deficiency.

Risk to human health related to the presence of perfluoroalkyl substances in food

The European Commission asked EFSA for a scientific evaluation on the risks to human health related to the presence of perfluoroalkyl substances (PFASs) in food. Based on several similar effects in animals, toxicokinetics and observed concentrations in human blood, the CONTAM Panel decided to perform the assessment for the sum of four PFASs: PFOA, PFNA, PFHxS and PFOS. These made up half of the lower bound (LB) exposure to those PFASs with available occurrence data, the remaining contribution being primarily from PFASs with short half-lives. Equal potencies were assumed for the four PFASs included in the assessment. The mean LB exposure in adolescents and adult age groups ranged from 3 to 22, the 95th percentile from 9 to 70 ng/kg body weight (bw) per week. Toddlers and 'other children' showed a twofold higher exposure. Upper bound exposure was 4- to 49-fold higher than LB levels, but the latter were considered more reliable. 'Fish meat', 'Fruit and fruit products' and 'Eggs and egg products' contributed most to the exposure. Based on available studies in animals and humans, effects on the immune system were considered the most critical for the risk assessment. From a human study, a lowest BMDL 10 of 17.5 ng/mL for the sum of the four PFASs in serum was identified for 1-year-old children. Using PBPK modelling, this serum level of 17.5 ng/mL in children was estimated to correspond to long-term maternal exposure of 0.63 ng/kg bw per day. Since accumulation over time is important, a tolerable weekly intake (TWI) of 4.4 ng/kg bw per week was established. This TWI also protects against other potential adverse effects observed in humans. Based on the estimated LB exposure, but also reported serum levels, the CONTAM Panel concluded that parts of the European population exceed this TWI, which is of concern.

Questioning the Utility of Pooling Samples in Microarray Experiments with Cell Lines

We describe a microarray experiment using the MCF-7 breast cancer cell line in two different experimental conditions for which the same number of independent pools as the number of individual samples was hybridized on Affymetrix GeneChips. Unexpectedly, when using individual samples, the number of probe sets found to be differentially expressed between treated and untreated cells was about three times greater than that found using pools. These findings indicate that pooling samples in microarray experiments where the biological variability is expected to be small might not be helpful and could even decrease one's ability to identify differentially expressed genes.

The PPARα-dependent rodent liver tumor response is not relevant to humans: addressing misconceptions

A number of industrial chemicals and therapeutic agents cause liver tumors in rats and mice by activating the nuclear receptor peroxisome proliferator-activated receptor α (PPARα). The molecular and cellular events by which PPARα activators induce rodent hepatocarcinogenesis have been extensively studied elucidating a number of consistent mechanistic changes linked to the increased incidence of liver neoplasms. The weight of evidence relevant to the hypothesized mode of action (MOA) for PPARα activator-induced rodent hepatocarcinogenesis is summarized here. Chemical-specific and mechanistic data support concordance of temporal and dose-response relationships for the key events associated with many PPARα activators. The key events (KE) identified in the MOA are PPARα activation (KE1), alteration in cell growth pathways (KE2), perturbation of hepatocyte growth and survival (KE3), and selective clonal expansion of preneoplastic foci cells (KE4), which leads to the apical event-increases in hepatocellular adenomas and carcinomas (KE5). In addition, a number of concurrent molecular and cellular events have been classified as modulating factors, because they potentially alter the ability of PPARα activators to increase rodent liver cancer while not being key events themselves. These modulating factors include increases in oxidative stress and activation of NF-kB. PPARα activators are unlikely to induce liver tumors in humans due to biological differences in the response of KEs downstream of PPARα activation. This conclusion is based on minimal or no effects observed on cell growth pathways and hepatocellular proliferation in human primary hepatocytes and absence of alteration in growth pathways, hepatocyte proliferation, and tumors in the livers of species (hamsters, guinea pigs and cynomolgus monkeys) that are more appropriate human surrogates than mice and rats at overlapping dose levels. Despite this overwhelming body of evidence and almost universal acceptance of the PPARα MOA and lack of human relevance, several reviews have selectively focused on specific studies that, as discussed, contradict the consensus opinion and suggest uncertainty. In the present review, we systematically address these most germane suggested weaknesses of the PPARα MOA.

Proteogenomics Reveals Enriched Ribosome Assembly and Protein Translation in Pectoralis major of High Feed Efficiency Pedigree Broiler Males

Background: In production animal agriculture, the cost of feed represents 60–70% of the total cost of raising an animal to market weight. Thus, development of viable biomarkers for feed efficiency (FE, g gain/g feed) to assist in genetic selection of breeding stock remains an important goal in commercial breeding programs.

Methods: Global gene (cDNA microarray, RNAseq) and protein expression (shotgun proteomics) analyses have been conducted on breast muscle samples obtained from pedigree broiler males (PedM) exhibiting high and low FE phenotypes. Using the entire datasets (i.e., no cutoffs for significance or fold difference in expression) the number of genes or proteins that were expressed numerically higher or lower in the high FE compared to the low FE phenotype for key terms or functions, e.g., ribosomal, mitochondrial ribosomal, tRNA, RNA binding motif, RNA polymerase, small nuclear ribonucleoprotein, and protein tyrosine phosphatase, were determined. Bionomial distribution analysis (exact) was then conducted on these datasets to determine significance between numerically up or down expression.

Results: Processes associated with mitochondrial proteome expression (e.g., mitochondrial ribosomal proteins, mitochondrial transcription, mitochondrial tRNA, and translation) were enriched in breast muscle from the high FE compared to the low FE pedigree male broiler phenotype. Furthermore, the high FE phenotype exhibited enrichment of ribosome assembly (e.g., RNA polymerase, mitochondrial and cytosolic ribosomes, small, and heterogeneous nuclear ribonucleoproteins), as well as nuclear transport and protein translation processes compared to the low FE phenotype. Quality control processes (proteosomes and autophagy) were also enriched in the high FE phenotype. In contrast, the low FE phenotype exhibited enrichment of cytoskeletal proteins, protein tyrosine phosphatases, and tyrosine kinases compared to the high FE phenotype. These results suggest that processes of mitochondrial and cytosolic ribosomal construction, activity, and protein translation would be enhanced in high FE breast muscle, and that phosphorylation of tyrosine moieties of proteins could be prolonged in the high compared to low FE phenotype. The results indicate the presence of a proteogenomic architecture that could enhance ribosome construction, protein translation, and quality control processes and contribute to the phenotypic expression of feed efficiency in this PedM broiler model.

Network-aided Bi-Clustering for discovering cancer subtypes

Bi-clustering is a widely used data mining technique for analyzing gene expression data. It simultaneously groups genes and samples of an input gene expression data matrix to discover bi-clusters that relevant samples exhibit similar gene expression profiles over a subset of genes. The discovered bi-clusters bring insights for categorization of cancer subtypes, gene treatments and others. Most existing bi-clustering approaches can only enumerate bi-clusters with constant values. Gene interaction networks can help to understand the pattern of cancer subtypes, but they are rarely integrated with gene expression data for exploring cancer subtypes. In this paper, we propose a novel method called Network-aided Bi-Clustering (NetBC). NetBC assigns weights to genes based on the structure of gene interaction network, and it iteratively optimizes sum-squared residue to obtain the row and column indicative matrices of bi-clusters by matrix factorization. NetBC can not only efficiently discover bi-clusters with constant values, but also bi-clusters with coherent trends. Empirical study on large-scale cancer gene expression datasets demonstrates that NetBC can more accurately discover cancer subtypes than other related algorithms.

The Effect of Acute and Chronic Social Stress on the Hippocampal Transcriptome in Mice

Psychogenic stress contributes to the formation of brain pathology. Using gene expression microarrays, we analyzed the hippocampal transcriptome of mice subjected to acute and chronic social stress of different duration. The longest period of social stress altered the expression of the highest number of genes and most of the stress-induced changes in transcription were reversible after 5 days of rest. Chronic stress affected genes involved in the functioning of the vascular system (Alas2, Hbb-b1, Hba-a2, Hba-a1), injury response (Vwf, Mgp, Cfh, Fbln5, Col3a1, Ctgf) and inflammation (S100a8, S100a9, Ctla2a, Ctla2b, Lcn2, Lrg1, Rsad2, Isg20). The results suggest that stress may affect brain functions through the stress-induced dysfunction of the vascular system. An important issue raised in our work is also the risk of the contamination of brain tissue samples with choroid plexus. Such contamination would result in a consistent up- or down-regulation of genes, such as Ttr, Igf2, Igfbp2, Prlr, Enpp2, Sostdc1, 1500015O10RIK (Ecrg4), Kl, Clic6, Kcne2, F5, Slc4a5, and Aqp1. Our study suggests that some of the previously reported, supposedly specific changes in hippocampal gene expression, may be a result of the inclusion of choroid plexus in the hippocampal samples.

The plasticizer benzyl butyl phthalate (BBP) alters the ecdysone hormone pathway, the cellular response to stress, the energy metabolism, and several detoxication mechanisms in Chironomus riparius larvae

Butyl benzyl phthalate (BBP) has been extensively used worldwide as a plasticizer in the polyvinyl chloride (PVC) industry and the manufacturing of many other products, and its presence in the aquatic environment is expected for decades. In the present study, the toxicity of BBP was investigated in Chironomus riparius aquatic larvae. The effects of acute 24-h and 48-h exposures to a wide range of BBP doses were evaluated at the molecular level by analysing changes in genes related to the stress response, the endocrine system, the energy metabolism, and detoxication pathways, as well as in the enzyme activity of glutathione S-transferase. BBP caused a dose and time-dependent toxicity in most of the selected biomarkers. 24-h exposures to high doses affected larval survival and lead to a significant response of several heat-shock genes (hsp70, hsp40, and hsp27), and to a clear endocrine disrupting effect by upregulating the ecdysone receptor gene (EcR). Longer treatments with low doses triggered a general repression of transcription and GST activity. Furthermore, delayed toxicity studies were specially relevant, since they allowed us to detect unpredictable toxic effects, not immediately manifested after contact with the phthalate. This study provides novel and interesting results on the toxic effects of BBP in C. riparius and highlights the suitability of this organism for ecotoxicological risk assessment, especially in aquatic ecosystems.

The genomic proliferation of transposable elements in colonizing populations: Schistosoma mansoni in the new world

Transposable elements (TEs) are mobile genes with an inherent ability to move within and among genomes. Theory predicts that TEs proliferate extensively during physiological stress due to the breakdown of TE repression systems. We tested this hypothesis in Schistosoma mansoni, a widespread trematode parasite that causes the human disease schistosomiasis. According to phylogenetic analysis, S. mansoni invaded the new world during the last 500 years. We hypothesized that new world strains of S. mansoni would have more copies of TEs than old world strains due to the physiological stress associated with invasion of the new world. We quantified the copy number of six TEs (Saci-1, Saci-2 and Saci-3, Perere-1, Merlin-sm1, and SmTRC1) in the genome and the transcriptome of old world and new world strains of S. mansoni, using qPCR relative quantification. As predicted, the genomes of new world parasites contain significantly more copies of class I and class II TEs in both laboratory and field strains. However, such differences are not observed in the transcriptome suggesting that either TE silencing mechanisms have reactivated to control the expression of these elements or the presence of inactive truncated copies of TEs.

Development and application of a comparative fatty acid analysis method to investigate voriconazole-induced hepatotoxicity

BACKGROUND

As fatty acids play an important role in biological regulation, the profiling of fatty acid expression has been used to discover various disease markers and to understand disease mechanisms. This study developed an effective and accurate comparative fatty acid analysis method using differential labeling to speed up the metabolic profiling of fatty acids.

METHODS

Fatty acids were derivatized with unlabeled (D0) or deuterated (D3) methanol, followed by GC-MS analysis. The comparative fatty acid analysis method was validated using a series of samples with different ratios of D0/D3-labeled fatty acid standards and with mouse liver extracts.

RESULTS

Using a lipopolysaccharide (LPS)-treated mouse model, we found that the fatty acid profiles after LPS treatment were similar between the conventional single-sample analysis approach and the proposed comparative approach, with a Pearson's correlation coefficient of approximately 0.96. We applied the comparative method to investigate voriconazole-induced hepatotoxicity and revealed the toxicity mechanism as well as the potential of using fatty acids as toxicity markers.

CONCLUSIONS

In conclusion, the comparative fatty acid profiling technique was determined to be fast and accurate and allowed the discovery of potential fatty acid biomarkers in a more economical and efficient manner.

Microarray RNA transcriptional profiling: Part I. Platforms, experimental design and standardization

This review summarizes, in a balanced and comprehensive manner, the various components of microarrays and their types, substrate architecture, platforms for microarray probe implementation, standardizations and confounders. The review is intended to familiarize the beginner with the principles of experimental design and the selection of an appropriate microarray platform. This parallel technology has revolutionized transcriptomic approaches to data profiling and has a major role in the identification of expressed genes, classification and diagnosis studies. The technology is still evolving and guidelines for standardization and reporting have been developed and are being improved.

Identification of serum insulin-like growth factor binding protein 1 as diagnostic biomarker for early-stage alcohol-induced liver disease

Background

Alcohol consumption is a major cause of liver disease in humans. The use and monitoring of biomarkers associated with early, pre-clinical stages of alcohol-induced liver disease (pre-ALD) could facilitate diagnosis and treatment, leading to improved outcomes.

Methods

We investigated the pathological, transcriptomic and protein changes in early stages of pre-ALD in mice fed the Lieber-Decarli liquid diet with or without alcohol for four months to identify biomarkers for the early stage of alcohol induced liver injury. Mice were sampled after 1, 2 and 4 months treatment.

Results

Pathological examination revealed a modest increase in fatty liver changes in alcohol-treated mice. Transcriptomics revealed gene alterations at all time points. Most notably, the Igfbp1 (Insulin-Like Growth Factor Binding Protein 1) was selected as the best candidate gene for early detection of liver damage since it showed early and continuously enhanced induction during the treatment course. Consistent with the microarray data, both Igfbp1mRNA expression in the liver tissue and the IGFBP1 serum protein levels showed progressive and significant increases over the course of pre-ALD development.

Conclusions

The results suggest that in conjunction with other tests, serum IGFBPI protein could provide an easily measured biomarker for early detection of alcohol-induced liver injury in humans.

Gene expression in breast muscle associated with feed efficiency in a single male broiler line using a chicken 44K microarray. II. Differentially expressed focus genes

Global RNA expression in breast muscle obtained from a male broiler line phenotyped for high or low feed efficiency (FE) was investigated using microarray analysis. Microarray procedures and validation were reported previously. By using an overlay function of a software program (Ingenuity Pathway Analysis, IPA) in which canonical pathways are projected onto a set of genes, a subset of 27 differentially expressed focus genes were identified. Focus genes that were upregulated in the high FE phenotype were associated with important signal transduction pathways (Jnk, G-coupled, and retinoic acid) or in sensing cell energy status and stimulating energy production that would likely enhance growth and development of muscle tissue. In contrast, focus genes that were upregulated in the low FE muscle phenotype were associated with cytoskeletal architecture (e.g., actin-myosin filaments), fatty acid oxidation, growth factors, or ones that would likely be induced in response to oxidative stress. The results of this study provide additional information on gene expression and the cellular basis of feed efficiency in broilers.

Effects of pooling samples on the performance of classification algorithms: a comparative study

A pooling design can be used as a powerful strategy to compensate for limited amounts of samples or high biological variation. In this paper, we perform a comparative study to model and quantify the effects of virtual pooling on the performance of the widely applied classifiers, support vector machines (SVMs), random forest (RF), k-nearest neighbors (k-NN), penalized logistic regression (PLR), and prediction analysis for microarrays (PAMs). We evaluate a variety of experimental designs using mock omics datasets with varying levels of pool sizes and considering effects from feature selection. Our results show that feature selection significantly improves classifier performance for non-pooled and pooled data. All investigated classifiers yield lower misclassification rates with smaller pool sizes. RF mainly outperforms other investigated algorithms, while accuracy levels are comparable among all the remaining ones. Guidelines are derived to identify an optimal pooling scheme for obtaining adequate predictive power and, hence, to motivate a study design that meets best experimental objectives and budgetary conditions, including time constraints.

RNA isolation for transcriptomics of human and mouse small skin biopsies

Background

Isolation of RNA from skin biopsies presents a challenge, due to the tough nature of skin tissue and a high presence of RNases. As we lacked the dedicated equipment, i.e. homogenizer or bead-beater, needed for the available RNA from skin isolation methods, we adapted and tested our zebrafish single-embryo RNA-isolation protocol for RNA isolation from skin punch biopsies.

Findings

We tested our new RNA-isolation protocol in two experiments: a large-scale study with 97 human skin samples, and a small study with 16 mouse skin samples. Human skin was sampled with 4.0 mm biopsy punches and for the mouse skin different punch diameter sizes were tested; 1.0, 1.5, 2.0, and 2.5 mm. The average RNA yield in human samples was 1.5 μg with an average RNA quality RIN value of 8.1. For the mouse biopsies, the average RNA yield was 2.4 μg with an average RIN value of 7.5. For 96% of the human biopsies and 100% of the mouse biopsies we obtained enough high-quality RNA. The RNA samples were successfully tested in a transcriptomics analysis using the Affymetrix and Roche NimbleGen platforms.

Conclusions

Using our new RNA-isolation protocol, we were able to consistently isolate high-quality RNA, which is apt for further transcriptomics analysis. Furthermore, this method is already useable on biopsy material obtained with a punch diameter as small as 1.5 mm.

Pathways to immunity: temporal dynamics of the bumblebee (Bombus terrestris) immune response against a trypanosomal gut parasite

Immune response dynamics in insects from natural host-parasite associations are poorly understood, despite accumulating evidence of ecological immune phenomena in these systems. Using a gene discovery approach, we have identified genes relating to signalling, enzymatic processes and respiration that were up-regulated in the bumblebee, Bombus terrestris, during infection with the trypanosomatid parasite, Crithidia bombi. In addition, we have mapped dynamic changes in the temporal expression of these genes and three candidate antimicrobial peptide (AMP) immune genes, Abaecin, Defensin and Hymenoptaecin, from 1 to 24 h after C. bombi infection. We show that dynamic changes in expression occur for individual genes at distinct phases of the immune response to C. bombi that correspond to early, intermediate and late stages of infection.

Regulation of Proteome Maintenance Gene Expression by Activators of Peroxisome Proliferator-Activated Receptor α

The nuclear receptor peroxisome proliferator-activated receptor α (PPARα) is activated by a large number of xenobiotic and hypolipidemic compounds called peroxisome proliferator chemicals (PPCs). One agonist of PPARα (WY-14,643) regulates responses in the mouse liver to chemical stress in part by altering expression of genes involved in proteome maintenance (PM) including protein chaperones in the heat shock protein (Hsp) family and proteasomal genes (Psm) involved in proteolysis. We hypothesized that other PPARα activators including diverse hypolipidemic and xenobiotic compounds also regulate PM genes in the rat and mouse liver. We examined the expression of PM genes in rat and mouse liver after exposure to 7 different PPCs (WY-14,643, clofibrate, fenofibrate, valproic acid, di-(2-ethylhexyl) phthalate, perfluorooctanoic acid, and perfluorooctane sulfonate) using Affymetrix microarrays. In rats and mice, 174 or 380 PM genes, respectively, were regulated by at least one PPC. The transcriptional changes were, for the most part, dependent on PPARα, as most changes were not observed in similarly treated PPARα-null mice and the changes were not consistently observed in rats treated with activators of the nuclear receptors CAR or PXR. In rats and mice, PM gene expression exhibited differences compared to typical direct targets of PPARα (e.g., Cyp4a family members). PM gene expression was usually delayed and in some cases, it was transient. Dose-response characterization of protein expression showed that Hsp86 and Hsp110 proteins were induced only at higher doses. These studies demonstrate that PPARα, activated by diverse PPC, regulates the expression of a large number of genes involved in protein folding and degradation and support an expanded role for PPARα in the regulation of genes that protect the proteome.

Establishing reference genes for use in real-time quantitative PCR analysis of early equine embryos

Real-time quantitative PCR (qPCR) is invaluable for investigating changes in gene expression during early development, since it can be performed on the limited quantities of mRNA contained in individual embryos. However, the reliability of this method depends on the use of validated stably expressed reference genes for accurate data normalisation. The aim of the present study was to identify and validate a set of reference genes suitable for studying gene expression during equine embryo development. The stable expression of four carefully selected reference genes and one developmentally regulated gene was examined by qPCR in equine in vivo embryos from morula to expanded blastocyst stage. SRP14, RPL4 and PGK1 were identified by geNorm analysis as stably expressed reference genes suitable for data normalisation. RPL13A expression was less stable and changed significantly during the period of development examined, rendering it unsuitable as a reference gene. As anticipated, CDX2 expression increased significantly during embryo development, supporting its possible role in trophectoderm specification in the horse. In summary, it was demonstrated that evidence-based selection of potential reference genes can reduce the number needed to validate stable expression in an experimental system; this is particularly useful when dealing with tissues that yield small amounts of mRNA. SRP14, RPL4 and PGK1 are stable reference genes suitable for normalising expression for genes of interest during in vivo morula to expanded blastocyst development of horse embryos.

Urine peptidomics for clinical biomarker discovery

Urine-based proteomic profiling is a novel approach that may result in the discovery of noninvasive biomarkers for diagnosing patients with different diseases, with the aim to ultimately improve clinical outcomes. Given new and emerging analytical technologies and data mining algorithms, the urine peptidome has become a rich resource to uncover naturally occurring peptide biomarkers for both systemic and renal diseases. However, significant analytical hurdles remain in sample collection and storage, experimental design, data analysis, and statistical inference. This study summarizes, focusing on our experiences and perspectives, the progress in addressing these challenges to enable high-throughput urine peptidomics-based biomarker discovery.

RNA isolation method for single embryo transcriptome analysis in zebrafish

Background

Transcriptome analysis during embryogenesis usually requires pooling of embryos to obtain sufficient RNA. Hence, the measured levels of gene-expression represent the average mRNA levels of pooled samples and the biological variation among individuals is confounded. This can irreversibly reduce the robustness, resolution, or expressiveness of the experiment. Therefore, we developed a robust method to isolate abundant high-quality RNA from individual embryos to perform single embryo transcriptome analyses using zebrafish as a model organism. Available methods for embryonic zebrafish RNA isolation minimally utilize ten embryos. Further downscaling of these methods to one embryo is practically not feasible.

Findings

We developed a single embryo RNA extraction method based on sample homogenization in liquid nitrogen, RNA extraction with phenol and column purification. Evaluation of this method showed that: the quality of the RNA was very good with an average RIN value of 8.3-8.9; the yield was always ≥ 200 ng RNA per embryo; the method was applicable to all stages of zebrafish embryogenesis; the success rate was almost 100%; and the extracted RNA performed excellent in microarray experiments in that the technical variation was much lower than the biological variation.

Conclusions

Presented is a high-quality, robust RNA isolation method. Obtaining sufficient RNA from single embryos eliminates the necessity of sample pooling and its associated drawbacks. Although our RNA isolation method has been setup for transcriptome analysis in zebrafish, it can also be used for other model systems and other applications like (q)PCR and transcriptome sequencing.

Characterization of peroxisome proliferator-activated receptor alpha--independent effects of PPARalpha activators in the rodent liver: di-(2-ethylhexyl) phthalate also activates the constitutive-activated receptor

Peroxisome proliferator chemicals (PPC) are thought to mediate their effects in rodents on hepatocyte growth and liver cancer through the nuclear receptor peroxisome proliferator-activated receptor (PPAR) alpha. Recent studies indicate that the plasticizer di-(2-ethylhexyl) phthalate (DEHP) increased the incidence of liver tumors in PPARalpha-null mice. We hypothesized that some PPC, including DEHP, induce transcriptional changes independent of PPARalpha but dependent on other nuclear receptors, including the constitutive-activated receptor (CAR) that mediates phenobarbital (PB) effects on hepatocyte growth and liver tumor induction. To determine the potential role of CAR in mediating effects of PPC, a meta-analysis was performed on transcript profiles from published studies in which rats and mice were exposed to PPC and compared the profiles to those produced by exposure to PB. Valproic acid, clofibrate, and DEHP in rat liver and DEHP in mouse liver induced genes, including Cyp2b family members that are known to be regulated by CAR. Examination of transcript changes by Affymetrix ST 1.0 arrays and reverse transcription-PCR in the livers of DEHP-treated wild-type, PPARalpha-null, and CAR-null mice demonstrated that (1) most (approximately 94%) of the transcriptional changes induced by DEHP were PPARalpha-dependent, (2) many PPARalpha-independent genes overlapped with those regulated by PB, (3) induction of genes Cyp2b10, Cyp3a11, and metallothionine-1 by DEHP was CAR dependent but PPARalpha-independent, and (4) induction of a number of genes (Cyp8b1, Gstm4, and Gstm7) was independent of both CAR and PPARalpha. Our results indicate that exposure to PPARalpha activators including DEHP leads to activation of multiple nuclear receptors in the rodent liver.

Toxicogenomics: transcription profiling for toxicology assessment

Toxicogenomics, the application of transcription profiling to toxicology, has been widely used for elucidating the molecular and cellular actions of chemicals and other environmental stressors on biological systems, predicting toxicity before any functional damages, and classification of known or new toxicants based on signatures of gene expression. The success of a toxicogenomics study depends upon close collaboration among experts in different fields, including a toxicologist or biologist, a bioinformatician, statistician, physician and, sometimes, mathematician. This review is focused on toxicogenomics studies, including transcription profiling technology, experimental design, significant gene extraction, toxicological results interpretation, potential pathway identification, database input and the applications of toxicogenomics in various fields of toxicological study.

Expression profiles of the pluripotency marker gene POU5F1 and validation of reference genes in rabbit oocytes and preimplantation stage embryos

BACKGROUND

The surge in the number of gene expression studies and tendencies to increase the quality of analysis have necessitated the identification of stable reference genes. Although rabbits are classical experimental model animals, stable reference genes have not been identified for normalization. The aims of this study were to compare the expression profiles of the widely used reference genes in rabbit oocytes and preimplantation stage embryos, and to select and validate stable ones to use as reference.

RESULTS

Quantitative real time PCR method was used to evaluate 13 commonly used references (Actb, Gapdh, Hprt1, H2afz, Ubc, Ppia, Eef1e1, Polr2a, Tbp, G6pdx, B2m, Pgk1, and Ywhaz) and POU5F1 (Oct4) genes. Expressions of these genes were examined in multiple individual embryos of seven different preimplantation developmental stages and embryo types (in vivo and in vitro). Initial analysis identified three genes (Ubc, Tbp, and B2m) close to the detection limit with irregular expression between the different stages. As variability impedes the selection of stable genes, these were excluded from further analysis. The expression levels of the remaining ten genes, varied according to developmental stage and embryo types. These genes were ranked using the geNorm software and finally the three most stable references (H2afz, Hprt1, and Ywhaz) were selected. Normalization factor was calculated (from the geometric averages of the three selected genes) and used to normalize the expressions of POU5F1 gene. The results showed the expected expression patterns of the POU5F1 during development.

CONCLUSION

Compared to the earlier studies with similar objectives, the comparison of large number of genes, the use of multiple individual embryos as compared to pools, and simultaneous analyses of in vitro and in vivo derived embryo samples were unique approaches in our study. Based on quantification, pattern and geNorm analyses, we found the three genes (H2afz, Hprt1, and Ywhaz) to be the most stable across developmental stages and embryo types, and the geometric averages of these genes can be used for appropriate normalization.

Gene expression in human peripheral blood mononuclear cells upon acute ischemic stroke

BACKGROUND AND PURPOSE

Ischemic stroke provokes a systemic inflammatory response. The purpose of this study was to characterize this response on the gene expression level in circulating mononuclear leukocytes from acute ischemic stroke (AIS) patients.

METHODS

RNA from peripheral blood mononuclear cells (PBMCs) of AIS patients (24 + 2 hours after onset of symptoms) was analyzed with Affymetrix U133A GeneChips using a pooled design. We compared the gene expression signature from AIS patients (n = 20), stroke survivors (n = 15), patients with acute traumatic brain injury (ATBI, n = 15) and healthy control subjects without vascular risk factors (n = 15).

RESULTS

Expression levels of 9682 probe sets with present calls on each GeneChip were compared. Between AIS patients and stroke survivors or between AIS patients and ATBI patients there were no significant differences in expression values of single genes after correction for multiple testing. However, comparison of the PBMC expression profiles from AIS patients and healthy subjects revealed significantly different expression (p = 0.012) of a single probe set, specific for phosphodiesterase 4 D (PDE4D). In order to detect modest expression differences in multiple genes with a presumed cumulative effect we studied the gene expression of functional groups of genes by global statistical tests. Analysis of 11 gene groups revealed differential expression between AIS patients and healthy subjects for genes involved in the inflammatory response (GeneOntology GO:0006954). Genes encoding the N-formyl peptide receptor-like 1 (FPRL1), interleukin-1 receptor antagonist (IL1RN) and complement component 3a receptor 1 (C3AR1) contributed most to the observed difference.

CONCLUSIONS

This transcriptome analysis did not identify significant changes between circulating mononuclear cells from AIS patients 24 hours after stroke and closely matched stroke survivors. However, comparing AIS patients with healthy control subjects revealed measurable differences in PDE4D and in inflammatory response genes when considered as a set.

Public consortium efforts in toxicogenomics

Public consortia provide a forum for addressing questions requiring more resources than one organization alone could bring to bear and engaging many sectors of the scientific community. They are particular well suited for tackling some of the questions encountered in the field of toxicogenomics, where the number of studies and microarray analyses would be prohibitively expensive for a single organization to carry out. Five consortia that stand out in the field of toxicogenomics are the Institutional Life Sciences Institute (ILSI) Health and Environmental Sciences Institute (HESI) Committee on the Application of Genomics to Mechanism Based Risk Assessment, the Toxicogenomics Research Consortium, the MicroArray Quality Control (MAQC) Consortium, the InnoMed PredTox effort, and the Predictive Safety Testing Consortium. Collectively, these consortia efforts have addressed issues such as reproducibility of microarray results, standard practice for assays and analysis, relevance of microarray results to conventional end points, and robustness of statistical models on diverse data sets. Their results demonstrate the impact that the pooling of resources, experience, expertise, and insight found in consortia can have.

Proteomic applications in ecotoxicology

Within the growing body of proteomics studies, issues addressing problems of ecotoxicology are on the rise. Generally speaking, ecotoxicology uses quantitative expression changes of distinct proteins known to be involved in toxicological responses as biomarkers. Unlike these directed approaches, proteomics examines how multiple expression changes are associated with a contamination that is suspected to be detrimental. Consequently, proteins involved in toxicological responses that have not been described previously may be revealed. Following identification of key proteins indicating exposure or effect, proteomics can potentially be employed in environmental risk assessment. To this end, bioinformatics may unveil protein patterns specific to an environmental stress that would constitute a classifier able to distinguish an exposure from a control state. The combined use of sets of marker proteins associated with a given pollution impact may prove to be more reliable, as they are based not only on a few unique markers which are measured independently, but reflect the complexity of a toxicological response. Such a proteomic pattern might also integrate some of the already established biomarkers of environmental toxicity. Proteomics applications in ecotoxicology may also comprise functional examination of known classes of proteins, such as glutathione transferases or metallothioneins, to elucidate their toxicological responses.

Literature-based compound profiling: application to toxicogenomics

Introduction: To reduce continuously increasing costs in drug development, adverse effects of drugs need to be detected as early as possible in the process. In recent years, compound-induced gene expression profiling methodologies have been developed to assess compound toxicity, including Gene Ontology term and pathway over-representation analyses. The objective of this study was to introduce an additional approach, in which literature information is used for compound profiling to evaluate compound toxicity and mode of toxicity. Methods: Gene annotations were built by text mining in Medline abstracts for retrieval of co-publications between genes, pathology terms, biological processes and pathways. This literature information was used to generate compound-specific keyword fingerprints, representing over-represented keywords calculated in a set of regulated genes after compound administration. To see whether keyword fingerprints can be used for assessment of compound toxicity, we analyzed microarray data sets of rat liver treated with 11 hepatotoxicants. Results: Analysis of keyword fingerprints of two genotoxic carcinogens, two nongenotoxic carcinogens, two peroxisome proliferators and two randomly generated gene sets, showed that each compound produced a specific keyword fingerprint that correlated with the experimentally observed histopathological events induced by the individual compounds. By contrast, the random sets produced a flat aspecific keyword profile, indicating that the fingerprints induced by the compounds reflect biological events rather than random noise. A more detailed analysis of the keyword profiles of diethylhexylphthalate, dimethylnitrosamine and methapyrilene (MPy) showed that the differences in the keyword fingerprints of these three compounds are based upon known distinct modes of action. Visualization of MPy-linked keywords and MPy-induced genes in a literature network enabled us to construct a mode of toxicity proposal for MPy, which is in agreement with known effects of MPy in literature. Conclusion: Compound keyword fingerprinting based on information retrieved from literature is a powerful approach for compound profiling, allowing evaluation of compound toxicity and analysis of the mode of action.

Technical, experimental, and biological variations in isobaric tags for relative and absolute quantitation (iTRAQ)

We assess the reliability of isobaric-tags for relative and absolute quantitation (iTRAQ), based on different types of replicate analyses taking into account technical, experimental, and biological variations. In total, 10 iTRAQ experiments were analyzed across three domains of life involving Saccharomyces cerevisiae KAY446, Sulfolobus solfataricus P2, and Synechocystis sp. PCC 6803. The coverage of protein expression of iTRAQ analysis increases as the variation tolerance increases. In brief, a cutoff point at +/-50% variation (+/-0.50) would yield 88% coverage in quantification based on an analysis of biological replicates. Technical replicate analysis produces a higher coverage level of 95% at a lower cutoff point of +/-30% variation. Experimental or iTRAQ variations exhibit similar behavior as biological variations, which suggest that most of the measurable deviations come from biological variations. These findings underline the importance of replicate analysis as a validation tool and benchmarking technique in protein expression analysis.

Quantitative evaluation and selection of reference genes in mouse oocytes and embryos cultured in vivo and in vitro

Background

Real-time PCR is an efficient tool to measure transcripts and provide valuable quantitative information on gene expression of preimplantation stage embryos. Finding valid reference genes for normalization is essential to interpret the real-time PCR results accurately, and understand the biological dynamics during early development. The use of reference genes also known as housekeeping genes is the most widely applied approach. However, the different genes are not systematically compared, and as a result there is no uniformity between studies in selecting the reference gene. The goals of this study were to compare a wide selection of the most commonly used housekeeping genes in mouse oocytes and preimplantation stage embryos produced under different culture conditions, and select the best stable genes for normalization of gene expression data.

Results

Quantitative real time PCR method was used to evaluate 12 commonly used housekeeping genes (Actb, Gapdh, H2afz, Hprt, Ppia, Ubc, Eef1e1, Tubb4, Hist2h2aa1, Tbp, Bmp7, Polr2a) in multiple individual embryos representing six different developmental stages. The results were analysed, and stable genes were selected using the geNorm software. The expression pattern was almost similar despite differences in the culture system; however, the transcript levels were affected by culture conditions. The genes have showed various stabilities, and have been ranked accordingly.

Conclusion

Compared to earlier studies with similar objectives, we used a unique approach in analysing larger number of genes, comparing embryo samples derived in vivo or in vitro, analysing the expression in the early and late maternal to zygote transition periods separately, and using multiple individual embryos. Based on detailed quantification, pattern analyses and using the geNorm application, we found Ppia, H2afz and Hprt1 genes to be the most stable across the different stages and culture conditions, while Actb, the classical housekeeping gene, showed the least stability. We recommend the use of the geometric averages of those three genes for normalization in mouse preimplantation-stage gene expression studies.

Underlying mechanisms of pharmacology and toxicity of a novel PPAR agonist revealed using rodent and canine hepatocytes

Marked species-specific responses to agonists of the peroxisome proliferator-activated alpha receptor (PPAR alpha) have been observed in rats and dogs, two species typically used to assess the potential human risk of pharmaceuticals in development. In this study, we used primary cultured rat and dog hepatocytes to investigate the underlying mechanisms of a novel PPAR alpha and -gamma coagonist, LY465608, relative to fenofibrate, a prototypical PPAR alpha agonist. As expected, rat hepatocytes incubated with these two agonists demonstrated an increase in peroxisome number as evaluated by electron microscopy, whereas the peroxisome number remained unchanged in dog hepatocytes. Biochemical analysis showed that rat hepatocytes responded to PPAR agonists with an induction of both peroxisomal and mitochondrial beta-oxidation (PBox and MBox) activities. Dog hepatocytes treated with both PPAR agonists, however, did not show increased PBox activity but did demonstrate increased MBox activity. Analysis of peroxisomal beta-oxidation gene expression markers by quantitative real-time PCR confirmed that PPAR agonists induced the peroxisomal enzymes, acyl-coenzyme A (CoA) oxidase (Acox), enoyl-CoA hydratase/L-3-hydroxyacyl-CoA dehydrogenase (Ehhadh), and 3-ketoacyl-CoA thiolase (Acaa1) at the transcriptional level in rat hepatocytes, but not dog hepatocytes. Expression of mRNA for the mitochondrial beta-oxidation gene hydroxyacyl-CoA dehydrogenase/3-ketoacyl-CoA thiolase (Hadhb), however, increased in both rat and dog hepatocytes, consistent with biochemical measurements of peroxisomal and mitochondrial beta-oxidation. Repeat-dose nonclinical safety studies of LY465608 revealed abnormities in mitochondrial morphology and evidence of single-cell necrosis following 30 days of dosing exclusively in dogs, but not in rats. Microarray analysis indicated that dog hepatocytes, but not rat hepatocytes, treated with LY465608 had an expression profile consistent with abnormalities in the regulation of cell renewal and death, oxidative stress, and mitochondrial bioenergetics, which may explain the canine-specific toxicity observed in vivo with this compound. This increased sensitivity to mitochondrial toxicity of canine hepatocytes relative to rat hepatocytes identified using gene expression was confirmed using the fluorescent indicator tetramethylrhodamine ethyl ester (TMRE) and flow cytometry. At doses of 0.1 microM LY465608, canine hepatocytes showed a greater shift in fluorescence indicative of mitochondrial damage than observed with rat hepatocytes treated at 10 microM. In summary, using rat and dog primary hepatocytes, we replicated the pharmacologic and toxicologic effects of LY465608 observed in vivo during preclinical development and propose an underlying mechanism for these species-specific effects.

Expression analysis of the acute phase response in channel catfish (Ictalurus punctatus) after infection with a Gram-negative bacterium

The acute phase response (APR) is a set of metabolic and physiological reactions occurring in the host in response to tissue infection or injury and is a crucial component of the larger innate immune response. The APR is best characterized by dramatic changes in the concentration of a group of plasma proteins known as acute phase proteins (APPs) which are synthesized in the liver and function in a wide range of immunity-related activities. Utilizing a new high-density in situ oligonucleotide microarray, we have evaluated the APR in channel catfish liver following infection with Edwardsiella ictaluri, a bacterial pathogen that causes enteric septicemia of catfish. Our catfish microarray design (28K) builds upon a previous 19K channel catfish array by adding recently sequenced immune transcripts from channel catfish along with 7159 unique sequences from closely related blue catfish. The analysis of microarray results using a traditional 2-fold change in gene expression cutoff and a 10% false-discovery rate revealed a well-developed APR in catfish, with particularly high upregulation (>50-fold) of genes involved in iron homeostasis (i.e. intelectin, hemopexin, haptoglobin, ferritin, and transferrin). Other classical APP genes upregulated greater than 2-fold included coagulation factors, proteinase inhibitors, transport proteins, and complement components. Upregulation of the majority of the complement cascade was observed including the membrane attack complex components and complement inhibitors. A number of pathogen recognition receptors (PRRs) and chemokines were also differentially expressed in the liver following infection. Independent testing of a selection of differentially expressed genes with real-time RT-PCR confirmed microarray results.

Gene expression profiles of murine fatty liver induced by the administration of valproic acid

Valproic acid (VPA) has been used as anticonvulsants, however, it induces hepatotoxicity such as microvesicular steatosis and necrosis in the liver. To explore the mechanisms of VPA-induced steatosis, we profiled the gene expression patterns of the mouse liver that were altered by treatment with VPA using microarray analysis. VPA was orally administered as a single dose of 100 mg/kg (low-dose) or 1000 mg/kg (high-dose) to ICR mice and the animals were killed at 6, 24, or 72 h after treatment. Serum alanine aminotransferase and aspartate aminotransferase levels were not significantly altered in the experimental animals. However, symptoms of steatosis were observed at 72 h with low-dose and at 24 h and 72 h with high-dose. After microarray data analysis, 1910 genes were selected by two-way ANOVA (P<0.05) as VPA-responsive genes. Hierarchical clustering revealed that gene expression changes depended on the time rather than the dose of VPA treatment. Gene profiling data showed striking changes in the expression of genes associated with lipid, fatty acid, and steroid metabolism, oncogenesis, signal transduction, and development. Functional categorization of 1156 characteristically up- and down-regulated genes (cutoff >1.5-fold) revealed that 60 genes were involved in lipid metabolism that was interconnected with biological pathways for biosynthesis of triglyceride and cholesterol, catabolism of fatty acid, and lipid transport. This gene expression profile may be associated with the known steatogenic hepatotoxicity of VPA and it may provide useful information for prediction of hepatotoxicity of unknown chemicals or new drug candidates through pattern recognition.

Gene expression microarrays and respiratory muscles

The routine measurement of the expression of tens of thousands of gene transcripts, simultaneously, is a defining advance of the last decade which has been made possible by microarray technology. Using this very powerful approach, a pattern has emerged from a number of studies that suggest a molecular niche for the diaphragm which is quite different from that occupied by limb muscle. All indications are that this is true not only in regard to differential gene transcription patterns in healthy muscles but also in the changes in transcription occurring in association with different diseases. Furthermore, respiratory muscle mounts a rich gene expression response to a number of disturbances, be they primary genetic defects (e.g. various types of muscular dystrophies) or non-genetic perturbations (e.g. controlled mechanical ventilation). Large numbers of genes undergo altered levels of transcription, ranging from tens to hundreds (typical) to thousands. These genes are involved in diverse cellular processes, such as contraction, intermediate metabolism, oxidative stress, apoptosis and cellular adhesion. Functional groups of genes identified as having changed expression differ in many respects from one disease to another. Previously identified pathways of muscle injury and repair are often perturbed to greater extents than previously anticipated, and processes not previously suspected of having important roles in the pathophysiology of specific disorders have been identified. Elucidation of these under-appreciated molecular events may lead to novel therapeutic interventions based on disrupting the downstream adverse consequences of the primary event or facilitating events which ameliorate the injury and/or promote muscle healing.

The hepatic transcriptome in human liver disease

The transcriptome is the mRNA transcript pool in a cell, organ or tissue with the liver transcriptome being amongst the most complex of any organ. Functional genomics methodologies are now being widely utilized to study transcriptomes including the hepatic transcriptome. This review outlines commonly used methods of transcriptome analysis, especially gene array analysis, focusing on publications utilizing these methods to understand human liver disease. Additionally, we have outlined the relationship between transcript and protein expressions as well as summarizing what is known about the variability of the transcriptome in non-diseased liver tissue. The approaches covered include gene array analysis, serial analysis of gene expression, subtractive hybridization and differential display. The discussion focuses on primate whole organ studies and in-vitro cell culture systems utilized. It is now clear that there are a vast number research opportunities for transcriptome analysis of human liver disease as we attempt to better understand both non-diseased and disease hepatic mRNA expression. We conclude that hepatic transcriptome analysis has already made significant contributions to the understanding of human liver pathobiology.

Probe-level linear model fitting and mixture modeling results in high accuracy detection of differential gene expression

Background

The identification of differentially expressed genes (DEGs) from Affymetrix GeneChips arrays is currently done by first computing expression levels from the low-level probe intensities, then deriving significance by comparing these expression levels between conditions. The proposed PL-LM (Probe-Level Linear Model) method implements a linear model applied on the probe-level data to directly estimate the treatment effect. A finite mixture of Gaussian components is then used to identify DEGs using the coefficients estimated by the linear model. This approach can readily be applied to experimental design with or without replication.

Results

On a wholly defined dataset, the PL-LM method was able to identify 75% of the differentially expressed genes within 10% of false positives. This accuracy was achieved both using the three replicates per conditions available in the dataset and using only one replicate per condition.

Conclusion

The method achieves, on this dataset, a higher accuracy than the best set of tools identified by the authors of the dataset, and does so using only one replicate per condition.

Gene Expression Profiles in Fathead Minnow Exposed to 2,4-DNT: Correlation with Toxicity in Mammals

Toxicogenomics, the genome-wide analysis of gene expression to study the effect of toxicants, has great potential for use in environmental toxicology. Applied to standard test organisms, it has possible applications in aquatic toxicology as a sensitive monitoring tool to detect the presence of contaminants while providing information on the mechanisms of action of these pollutants. We describe the use of a complementary DNA (cDNA) microarray of the fathead minnow (Pimephales promelas) a standard sentinel organism in aquatic toxicology, to better understand the mechanisms of toxicity of 2,4-dinitrotoluene (2,4-DNT) which is released in the environment through military and industrial use. We have constructed a fathead minnow microarray containing 5000 randomly picked anonymous cDNAs from a whole fish cDNA library. Expression profiles were analyzed in fish exposed to 2,4-DNT for 10 days at three concentrations (11, 22, and 44 microM, respectively) below the measured median lethal concentration (58 microM). Sequence analysis of cDNAs corresponding to differentially expressed genes affected by exposure revealed that lipid metabolism and oxygen transport genes were prominently affected in a dose-specific manner. We measured liver lipids and demonstrate that lipid metabolism is indeed perturbed following exposure. These observations correlate well with available toxicological data on 2,4-DNT. We present possible modes of action of 2,4-DNT toxicity and suggest that fathead minnow cDNA microarrays can be useful to identify mechanisms of toxicity in fish and as a predictive tool for toxicity in mammals.

Differential gene expression in the rat cochlea after exposure to impulse noise

Understanding the molecular biology of noise trauma is vital to developing effective and timely interventions. In a model of explosion-mediated impulse noise injury, differential gene expression was studied in whole rat cochlea preparations at 3 and 24 h following the exposure. We developed a technique using mRNA from a single cochlea on each oligonucleotide microarray to avoid pooling of mRNA samples. Application of a conservative statistical analysis approach resulted in the identification of 61 differentially expressed genes. Within 3 h after the exposure, there was an up-regulation of immediate early genes, mainly transcription factors and genes involved in the tissue's response to oxidative stress. No genes were found to be significantly down-regulated. At 24 h following the exposure, up-regulated genes included members of inflammatory and antioxidant pathways and one gene involved in glutathione metabolism was down-regulated. A subset of genes was confirmed by real-time reverse transcriptase-polymerase chain reaction (RT-PCR). The present study demonstrates the power of the microarray technique in providing a global view of the gene regulation following noise exposure, and in identifying genes that may be mechanistically important in hearing loss, and thereby serve as a basis for the development of therapeutic interventions.

Genomics of Hepatitis B and C Infections: Diagnostic and Therapeutic Applications of Microarray Profiling

Microarray profiling offers many potential advances in diagnostic and therapeutic intervention in human disease because of its unparalleled ability to conduct high-throughput analysis of gene expression. However, limitations of this technique relate in part to issues regarding the various methodologies and experimental designs as well as difficulties in the interpretation of results. Despite this, microarray profiling has led to a better understanding of the molecular pathogenesis of hepatitis B virus (HBV) and hepatitis C virus (HCV) infection. Key events in clearance and the development of chronicity of HCV have been identified that may prove to have a role in the development of future treatments. In addition, pharmacogenomic studies of interferon-based treatment for chronic HCV and HBV have provided mechanistic insights into the therapeutic action of interferons. These advances have implications with respect to the development of improved therapeutic agents. New biomarkers for cancer screening and gene profiles with prognostic value for survival have also been developed for hepatocellular carcinoma, which frequently complicates chronic viral hepatitis. Thus, microarray profiling offers enormous potential for improvements in antiviral therapy and our understanding of blood-borne viral hepatitis.

An Integrated Study of Acute Effects of Valproic Acid in the Liver Using Metabonomics, Proteomics, and Transcriptomics Platforms

An integrated omics approach was undertaken in order to elucidate a systems biology level understanding of the acute hepatotoxcity of valproic acid (VPA). Metabonomics, proteomics and gene expression microarray platforms were employed in this systems biology study. CD-1 female pregnant mice were injected subcutaneously with 600 mg/kg VPA or vehicle control. Urine, serum, and liver tissue were collected at 6, 12, and 24 h after dosing. Principal component analysis (PCA) of the metabonomics data showed clustering of the dosed groups away from the controls for the urine samples. Looser clustering was seen in the other sample sets investigated. However, VPA administration resulted in altered glucose concentrations in urine samples at 12 and 24 h and in aqueous liver tissue extracts at 12 h after VPA administration. Proteomics studies identified two proteins, glycogen phosphorylase and amylo-1,6-glucosidase, which were increased in dosed animals relative to control. Both of these proteins are involved in converting glycogen to glucose. Examination of the expression of 20,000 liver genes did not reveal significantly altered expression at 6, 12, or 24 h after VPA exposure. The combined studies indicated a perturbation in the glycogenolysis pathway following administration of VPA.

Gene expression profiling of diaphragm muscle in alpha2-laminin (merosin)-deficient dy/dy dystrophic mice

Deficiency of alpha2-laminin (merosin) underlies classical congenital muscular dystrophy in humans and dy/dy muscular dystrophy in mice and causes severe muscle dysfunction in both species. To gain greater insight into the biochemical and molecular events that link alpha2-laminin deficiency with muscle fiber necrosis, and the associated compensatory responses, gene expression profiles were characterized in diaphragm muscle from 8-wk-old dy/dy mice using oligonucleotide microarrays. Compared with age-matched normal muscle, dystrophic diaphragm was characterized by predominantly augmented gene expression, irrespective of the fold-change threshold. Among the 69 genes with at least plus or minus twofold significantly altered expression, 30 belonged to statistically overrepresented Gene Ontology (GO) biological process groups. These covered four specific themes: development including muscle development, cell motility with an emphasis on muscle contraction, defense/immune response, and cell adhesion. An additional 11 gene transcripts were assigned to more general overrepresented GO biological process groups (e.g., cellular process, organismal physiological process); the remaining 28 did not belong to any overrepresented groups. GO cellular constituent assignment resulted in the highest degree of overrepresentation in extracellular and muscle fiber locations, whereas GO molecular function assignment was most notable for various types of binding. RT-PCR was performed on 38 of 41 genes with at least plus or minus twofold significantly altered expression that were assigned to overrepresented GO biological process groups, with expression changes verified for 36 of 38 genes. These results indicate that several specific groups of genes have altered expression in response to genetic alpha2-laminin deficiency, with both similarities and differences compared with data reported for dystrophin-deficient muscular dystrophies.

Effect of pooling samples on the efficiency of comparative studies using microarrays

MOTIVATION

Many biomedical experiments are carried out by pooling individual biological samples. However, pooling samples can potentially hide biological variance and give false confidence concerning the data significance. In the context of microarray experiments for detecting differentially expressed genes, recent publications have addressed the problem of the efficiency of sample pooling, and some approximate formulas were provided for the power and sample size calculations. It is desirable to have exact formulas for these calculations and have the approximate results checked against the exact ones. We show that the difference between the approximate and the exact results can be large.

RESULTS

In this study, we have characterized quantitatively the effect of pooling samples on the efficiency of microarray experiments for the detection of differential gene expression between two classes. We present exact formulas for calculating the power of microarray experimental designs involving sample pooling and technical replications. The formulas can be used to determine the total number of arrays and biological subjects required in an experiment to achieve the desired power at a given significance level. The conditions under which pooled design becomes preferable to non-pooled design can then be derived given the unit cost associated with a microarray and that with a biological subject. This paper thus serves to provide guidance on sample pooling and cost-effectiveness. The formulation in this paper is outlined in the context of performing microarray comparative studies, but its applicability is not limited to microarray experiments. It is also applicable to a wide range of biomedical comparative studies where sample pooling may be involved.