Feasibility of a liver transcriptomics approach to assess bovine treatment with the prohormone dehydroepiandrosterone (DHEA)
© Rijk et al; licensee BioMed Central Ltd. 2010
Received: 15 March 2010
Accepted: 16 September 2010
Published: 16 September 2010
Within the European Union the use of growth promoting agents in animal production is prohibited. Illegal use of natural prohormones like dehydroepiandrosterone (DHEA) is hard to prove since prohormones are strongly metabolized in vivo. In the present study, we investigated the feasibility of a novel effect-based approach for monitoring abuse of DHEA. Changes in gene expression profiles were studied in livers of bull calves treated orally (PO) or intramuscularly (IM) with 1000 mg DHEA versus two control groups, using bovine 44K DNA microarrays. In contrast to controlled genomics studies, this work involved bovines purchased at the local market on three different occasions with ages ranging from 6 to 14 months, thereby reflecting the real life inter-animal variability due to differences in age, individual physiology, season and diet.
As determined by principal component analysis (PCA), large differences in liver gene expression profiles were observed between treated and control animals as well as between the two control groups. When comparing the gene expression profiles of PO and IM treated animals to that of all control animals, the number of significantly regulated genes (p-value <0.05 and a fold change >1.5) was 23 and 37 respectively. For IM and PO treated calves, gene sets were generated of genes that were significantly regulated compared to one control group and validated versus the other control group using Gene Set Enrichment Analysis (GSEA). This cross validation, showed that 6 out of the 8 gene sets were significantly enriched in DHEA treated animals when compared to an 'independent' control group.
This study showed that identification and application of genomic biomarkers for screening of (pro)hormone abuse in livestock production is substantially hampered by biological variation. On the other hand, it is demonstrated that comparison of pre-defined gene sets versus the whole genome expression profile of an animal allows to distinguish DHEA treatment effects from variations in gene expression due to inherent biological variation. Therefore, DNA-microarray expression profiling together with statistical tools like GSEA represent a promising approach to screen for (pro)hormone abuse in livestock production. However, a better insight in the genomic variability of the control population is a prerequisite in order to define growth promoter specific gene sets that can be used as robust biomarkers in daily practice.
In the European Union the use of growth promoting substances in livestock production is prohibited according EC directive 96/22 . To ensure compliance with this legislation, requirements for monitoring are described in EC directive 96/23 . At national level, legislations are implemented in residue monitoring programs regulating sampling of animal matrices and residue analysis therein to guarantee fair trade, food safety and public health. Residue analysis in livestock production is in general based on chemical , immunochemical or biological [4, 5] screening methods followed by mass spectrometry based confirmation methods. Although this strategy seems to work for synthetic anabolic steroids, problems arise when compounds that also occur naturally are used.
Abuse of naturally occurring (pro)hormones is hard to prove since most of these substances are strongly metabolized in vivo. Moreover, metabolites are not always known or are present in levels not significantly different from highly fluctuating endogenous levels. This makes it difficult to prove fraudulent use based on quantification of natural occurring compounds. Nowadays, it is observed that misuse of growth promoters in cattle fattening moves towards these natural steroids and steroid esters. Moreover, inspections of livestock farms in The Netherlands occasionally result in the finding of feed or herbal additives and preparations containing so-called prohormones. Prohormones are compounds that exhibit limited or no hormonal action by themselves, however they are direct precursors of active hormones and indirectly affect natural hormone levels. Dehydroepiandrosterone (DHEA) is such a prohormone and is the most abundant occurring precursor of both androgens and estrogens in humans [6, 7]. It is claimed that orally taken DHEA improves muscle strength and is therefore illicitly used in sports to enhance performance and appearance [8, 9].
Looking for alternatives to support evidence of illegal use of growth promoting substances, gene expression analysis can be an attractive new approach. Several studies demonstrated changes in mRNA expression in bovine tissues upon treatment with growth promoters after performing real-time RT-PCR analysis on a limited number of preselected genes [10–14]. Untargeted transcriptomics approaches using microarrays allow gene expression analysis of thousands of genes simultaneously as well as identification of (new) biomarkers for screening [15, 16]. Moreover, microarray data can provide mechanistic insights in cellular processes and pathways and can be used for classification of compounds with the same mode of action (gene expression finger prints) [17, 18]. Comparative microarray analysis is therefore in potential a promising screening tool for growth promoter abuse and in particular for prohormones of which the mode of action in cattle is sometimes unclear.
In recent work we used a metabolomics approach to compare urine profiles of control and DHEA exposed bovines . This revealed several urinary steroid phase I and phase II metabolites which are potential biomarkers for DHEA treatment. In the present study we investigated the feasibility of monitoring prohormone abuse at the mRNA level using liver tissue from the same animal experiment. Gene expression profiles of control and DHEA treated animals were compared to determine differentially expressed genes and to identify biomarkers for DHEA treatment.
Animals and treatment
Experimental setup, age and weights of bovines included in the DHEA animal treatment experiment.
Intra muscular (IM 1)
Oral (PO 1)
Control (C 1-1)
Control (C 1-2)
Control (C 1-3)
Intra muscular (IM 2)
Oral (PO 2)
Intra muscular (IM 3)
Oral (PO 3)
Control (C 3-1)
Control (C 3-2)
Control (C 3-3)
Control (C 3-4)
Total RNA was extracted from tissues by homogenization in Trizol (Invitrogen Life Technologies, Breda, The Netherlands) and mixed with chloroform. The lysate was centrifuged at 12000 × g for 15 minutes at 4°C and the aqueous phase was transferred to be mixed with isopropanol which precipitates total RNA. After centrifuging (10 minutes, 12000 xg at 4°C) the pellet was washed with 75% ethanol and resuspended in RNase free water. Upon extraction the RNA was purified according to the RNeasy mini kit protocol (Qiagen, Westburg bv, Leusden, The Netherlands). After purification, RNA integrity was determined spectroscopically (Nanodrop technologies) and by gel electrophoresis. Only RNA with A260/280 and A260/230 ratios above 1.8 was used for amplification. To generate fluorescently-labeled cRNA, the Agilent Low RNA Input Fluorescent Linear Amplification Kit (Agilent Technologies, Palo Alto, CA, USA) was used according to the manufacturer's protocol. In short, 1 μg of total RNA was reverse transcribed using T7 tagged oligo-dT primer and labeled with Cy3 or Cy5 (Perkin Elmer/NEN Life Sciences, Boston, MA, USA). Liver RNAs of the treated and control animals were individually labeled with Cy5 and RNA of all 7 control animals was pooled and labeled with Cy3. After purification with the RNeasy mini kit (Qiagen), label efficiency and yield were determined using a Nanodrop spectrophotometer (Nanodrop technologies). A mixture of 1 μg of Cy3-labeled cRNA and 1 μg of Cy5-labeled cRNA was hybridized onto a 44k bovine oligo microarray (Agilent Technologies), using Agilent's gene expression hybridization kit. Hybridization was performed at 65°C for 17 hours in a hybridization oven with rotation function (Agilent Technologies). Upon hybridization, microarrays were washed and dried according to the Agilent's instructions and fluorescence measurements were performed using an Agilent Technologies G2565B microarray scanner.
Fluorescence intensities were quantified using Feature Extraction 8.5 software (Agilent Technologies). Data were imported in GeneMaths XT 1.6 (Applied Maths, St. Martens-Latem, Belgium) and signals below two times background were excluded from further analysis. Subsequently, the data was normalized as described by Pellis et al. . This normalization included correction for the random error, with the median Cy3 signal for each individual spot. Secondly, correction for the systematic error was performed with the median value of the overall Cy5 signal. After normalization, principal component analysis (PCA) was performed to visualize differences between groups and t-test statistics were performed to test for differential expression. Microarray data was floored by adjusting low intensity spots to a threshold value of 130, hereby reducing the number of less reliable genes. Next, spot intensities were 2log transformed and each gene was mean centered versus all samples. Based on these 2log transformed data differentially regulated genes were selected with a p-value <0.05 and a fold change >1.5 (> 2log 0.6) in each of the three treatment replicates versus the average from the control animals. Hierarchical clustering of the differentially regulated genes was performed using Cluster and Treeview software . Raw microarray data of the present study have been submitted to ArrayExpress (available at: http://www.ebi.ac.uk) and are stored under experiment accession number A-MEXP-1810.
Gene set enrichment analysis (GSEA)
Gene set enrichment analysis (GSEA) is a tool to identify and analyse the differential expression of biologically relevant sets of genes that share common biological functions . Using GSEA, the differentially regulated genes observed in DHEA treated animals versus one control group (e.g. controls of experiment 1) were validated by evaluating this gene set by comparing the same DHEA treated animals versus the other control group (e.g. controls of experiment 3). Therefore, separate gene sets were generated of the differentially expressed genes of respectively IM treated animals (n = 3) and OS treated animals (n = 3) versus the controls of experiment 1 as well as the controls of experiment 3. For example, the transcripts found significantly up-regulated when comparing DHEA IM versus control group 1 were included in the gene set "DHEA_IM_vs_CTR1_UP". In a similar way other gene sets were created for up- as well as down-regulated genes. Next, GSEA ranks all genes on the microarray on differential expression between DHEA exposed and controls using signal to noise statistics, resulting in a list with up-regulated genes at the top and down-regulated genes at the lower end of the list. Each of the pre-defined gene sets was tested against this list and GSEA calculated whether the genes in the gene set are randomly distributed, enriched at the top or at the lower end of the ranked list. Permutations were performed on gene sets and gene sets were considered significantly affected when the p-value was below 0.05 and the false discovery rate (FDR) below 0.25, accordingly to GSEA recommendations .
Results and discussion
Principal component analysis (PCA) and selection of differentially regulated genes
In the present study the potential strengths as well as the pitfalls of microarray experiments using calves from real-life practice were investigated. Three small animal experiments were performed independently using bull calves purchased at the local market. In this way the experimental setup was taking into account the inherent variability needed to investigate the usefulness of bovine-specific microarrays as a screening tool for prohormone abuse in veterinary control. For obvious ethical reasons larger numbers of bovines treated with banned substances could not be justified.
Of the 37 differentially expressed genes in response to the IM DHEA treatment, 4 represented unidentified transcripts and 7 were encoding for proteins with poorly known or unknown function (LOC617652, LOC527553, LOC515784, LOC617667, LOC515640, LOC497203 and LOC521795). Among the 26 transcripts that encode for known proteins, 10 transcripts are involved in immune response and inflammatory processes. Of these latter transcripts XBP1, MX1, LBP, SERPINA3, CCL24, CARD14 and PIGR were found up-regulated and ANKRD1, LYZ and DMBT were down-regulated. The remaining transcripts are involved in various processes like cell growth and proliferation (INHBE), formation of tight-junctions (CLDN7), tumor suppression (DIRAS3), cell proliferation and cell adhesion (CYR61) intra-cellular signalling (JAG2) and cell-cell interactions (LGALS4). Regarding metabolism, the GSTP1 gene was found to be down regulated >1.6 fold in all IM treated animals. GSTP1 mediates glutathione conjugation and plays an important role in detoxification of xenobiotics as well as in uptake and transport of numerous hydrophobic endogenous compounds like steroids  Moreover, it has been observed in mouse that the GSTP1 gene contains androgen receptor binding sites which regulate GSTP1 activity in response to androgens .
Comparison of PO treated animals versus all control animals revealed a total of 23 differentially expressed genes of which 7 represent unidentified transcripts or encoded for proteins with an unknown function. Again a substantial number of the differentially regulated genes are involved in immune response of which LILRA5, THBS, CLEC6A and FUT1 were found up-regulated and CCL14 and DBMT were down-regulated. Other differentially regulated genes are involved in peptidase inhibition (SERPINB8, WFDC1), G-protein signalling (RGS5) and amino acid transport (SLC6A14). Also regulated is the short/branched chain acyl-CoA dehydrogenase (ACADSB) gene, a member of the acyl-CoA dehydrogenase enzyme family which is involved in fatty acid metabolism. This could point towards regulation of fatty acid metabolism and is supported by a study in which DHEA administration to rats showed significant regulation of genes involved in fatty acid metabolism, including the very long chain acyl-CoA gene which is also a member of the acyl-CoA dehydrogenase enzyme family . Overall it can be stated that the majority of regulated genes are involved in immune response for both PO as well as IM treated animals which is in line with numerous studies reporting the significant immune modulatory properties of DHEA .
In principle the above listed genes are potential biomarkers for DHEA treatment. On the one hand, we are aware of the small number of animals used in this study which hampers proper statistics and substantially increases the chance of missing DHEA-responsive genes or detecting false-positive genes. On the other hand, combining and comparing the data of three independently performed experiments will limit the risk of false-positive genes considerably and results in identification of only the most robustly regulated genes. Therefore, we assessed whether the genes differentially expressed in animals treated with DHEA via one administration route versus animals of one control group would also be affected when compared 1) to other control animals and 2) by the other administration route. To deal with these issues we applied the statistics of gene set enrichment analysis (GSEA). In this way statistical power could be improved and regulated gene sets were tested for their robustness.
Gene set enrichment analysis (GSEA)
Significance of gene set regulation after GSEA analysis.
DHEA IM group compared to CTR1
DHEA IM group compared to CTR3
Name gene set
DHEA PO group compared to CTR1
DHEA PO group compared to CTR3
Name gene set
Although in vivo transcriptomics data of DHEA in liver is limited [25, 27], Depreter et al. identified 13 genes which were found to be up-regulated in rat liver . GSEA analysis showed significant enrichment of this gene set in DHEA IM and PO treated animals when these were compared with the controls of experiment 1 (Table 2). These results illustrate that GSEA is a powerful approach for comparative analysis of gene expression data obtained in different settings.
Controlled experiments with bovines have resulted in the identification of genomics based biomarkers which potentially can be used for screening for hormones [10–14, 28]. However, when examining bovines from real-life practice, one is dealing with biological variation like age, genetic background, environment, nutrition and disease history. In the current study, this biological variation was deliberately included and was mainly reflected by the large differences in gene expression profiles of the control populations tested. The two control groups in this study showed substantial age differences i.e. the animals in the CTR1, CTR3 and DHEA-treated groups are 6 months, 13-14 months and 8-13 months in age, respectively. Nevertheless, for the DHEA IM and PO treated animals, sets of respectively 37 and 23 genes were found differentially expressed when compared to all controls using standard statistics. These two groups of genes are specific for IM and PO treatment, respectively, and independent from biological factors like age. However, GSEA results showed a correlation between gene expression profiles of IM and PO treated animals, suggesting that there are also effects independent from the route of administration. This is in line with our earlier performed metabolomics study showing large similarities in urine metabolite profiles of IM and PO treated animals as well as metabolites specific for the route of administration .
Hence, for application of transcriptomics based screening of bovines for (pro)hormones in practice, the treatment effect should be filtered out from differences in gene expression due to inherent biological variation. Here it was shown that microarray gene expression profiling in combination with statistical methods like GSEA are able to distinguish gene expression profiles of DHEA-treated animals from non-treated control animals. It should be noted that this experiment comprised small numbers of animals and follow up experiments are required to gain statistical power and to obtain a better description of DHEA specific gene sets. Furthermore, the behaviour of such a growth promoter specific gene set should be studied in a broad spectrum of untreated control animals from daily practice, to assure the robustness of the gene set. This underlines the need to obtain more liver gene expression profiles of control animals from slaughterhouses.
The present study showed that identification of genomic biomarkers for DHEA treatment in cattle is hampered by the large biological variability as compared to genomics experiments with inbred strains of rodents under well-defined laboratory conditions. However, gene expression profiling using whole genome microarrays in combination with predefined gene sets and statistical methods like GSEA showed to be a promising approach to screen for (pro)hormone abuse in livestock production. For application in practice however, a better genomic description of the control population as well as growth promoter specific gene set are needed.
This project was financially supported by the Dutch Ministry of Agriculture, Nature and Food Quality.
- EEC Directive 96/22. Off J Eur Commun. 1996, L125: 3-9.Google Scholar
- EEC Directive 96/23. Off J Eur Commun. 1996, L125: 10-32.Google Scholar
- Stolker AAM, Zuidema T, Nielen MWF: Residue analysis of veterinary drugs and growth-promoting agents. Trends Anal Chem. 2007, 26: 967-979. 10.1016/j.trac.2007.09.008.View ArticleGoogle Scholar
- Simontacchi C, Perez De Altamirano T, Marinelli L, Angeletti R, Gabai G: Plasma steroid variations in bull calves repeatedly treated with testosterone, nortestosterone and oestradiol administered alone or in combination. Vet Res Comm. 2004, 28: 467-477. 10.1023/B:VERC.0000040244.27933.f1.View ArticleGoogle Scholar
- Bovee TFH, Pikkemaat MG: Bioactivity-based screening of antibiotics and hormones. J Chromatogr A. 2009, 1216: 8035-8050. 10.1016/j.chroma.2009.03.045.PubMedView ArticleGoogle Scholar
- Labrie F: Intracrinology. Mol Cell Endocrinol. 1991, 78: 113-118. 10.1016/0303-7207(91)90116-A.View ArticleGoogle Scholar
- Labrie F, Belanger A, Luu-The V, Labrie C, Simard J, Cusan L, Gomez J-L, Candas B: DHEA and the intracrine formation of androgens and estrogens in peripheral target tissues: Its role during aging. Steroids. 1998, 63: 322-328. 10.1016/S0039-128X(98)00007-5.PubMedView ArticleGoogle Scholar
- Dehennin L, Ferry M, Lafarge P, Peres G, Lafarge J-P: Oral administration of dehydroepiandrosterone to healthy men: Alteration of the urinary androgen profile and consequences for the detection of abuse in sport by gas chromatography-mass spectrometry. Steroids. 1998, 63: 80-87. 10.1016/S0039-128X(97)00138-4.PubMedView ArticleGoogle Scholar
- Bahrke MS, Yesalis CE: Abuse of anabolic androgenic steroids and related substances in sport and exercise. Curr Opin Pharmacol. 2004, 4: 614-620. 10.1016/j.coph.2004.05.006.PubMedView ArticleGoogle Scholar
- Reiter M, Walf VM, Christians A, Pfaffl MW, Meyer HHD: Modification of mRNA expression after treatment with anabolic agents and the usefulness for gene expression-biomarkers. Anal Chim Acta. 2007, 586: 73-81. 10.1016/j.aca.2006.10.049.PubMedView ArticleGoogle Scholar
- Toffolatti L, Rosa Gastaldo L, Patarnello T, Romualdi C, Merlanti R, Montesissa C, Poppi L, Castagnaro M, Bargelloni L: Expression analysis of androgen-responsive genes in the prostate of veal calves treated with anabolic hormones. Dom Anim Endocrinol. 2005, 30: 38-55. 10.1016/j.domaniend.2005.05.008.View ArticleGoogle Scholar
- Riedmaier I, Becker C, Pfaffl MW, Meyer HHD: The use of omic technologies for biomarker development to trace functions of anabolic agents. J Chromatogr A. 2009, 1216: 8192-8199. 10.1016/j.chroma.2009.01.094.PubMedView ArticleGoogle Scholar
- Pfaffl MW, Daxenberger A, Hageleit M, Meyer HHD: Effects of synthetic progestagens on the mRNA expression of androgen receptor, progesterone receptor, oestrogen receptor α and ß, insulin-like growth factor-1 (IGF-1) and IGF-1 receptor in heifer tissues. J Vet Med A Physiol Pathol Clin Med. 2002, 49: 57-64.PubMedView ArticleGoogle Scholar
- White ME, Johnson BJ, Hathaway MR, Dayton WR: Growth factor messenger RNA levels in muscle and liver of steroid-implanted and nonimplanted steers. J Anim Sci. 2003, 81: 965-972.PubMedGoogle Scholar
- Buterin T, Koch C, Naegeli H: Potential application of gene expression fingerprinting for food safety screening. Anal Chim Acta. 2005, 529: 33-39. 10.1016/j.aca.2004.07.012.View ArticleGoogle Scholar
- Carraro L, Ferraresso S, Cardazzo B, Romualdi C, Montesissa C, Gottardo F, Patarnello T, Castagnaro M, Bargelloni L: Expression profiling of skeletal muscle in young bulls treated with steroidal growth promoters. Physiol Genomics. 2009, 38: 138-148. 10.1152/physiolgenomics.00014.2009.PubMedView ArticleGoogle Scholar
- Labrie F, Luu-The V, Calvo E, Martel C, Cloutier J, Gauthier S, Belleau P, Morisette J, Levesque M-H, Labrie C: Tetrahydrogestrinone induces a genomic signature typical of a potent anabolic steroid. J Endocrinol. 2005, 184: 427-433. 10.1677/joe.1.05997.PubMedView ArticleGoogle Scholar
- Nuwaysir EF, Bittner M, Trent J, Barret JC, Afshari CA: Microarrays and toxicology: the advent of toxicogenomics. Mol Carcinog. 1999, 24: 153-159. 10.1002/(SICI)1098-2744(199903)24:3<153::AID-MC1>3.0.CO;2-P.PubMedView ArticleGoogle Scholar
- Rijk JCW, Lommen A, Essers ML, Groot MJ, Van Hende JM, Doeswijk TG, Nielen MWF: Metabolomics approach to anabolic steroid urine profiling of bovines treated with prohormones. Anal Chem. 2009, 81: 6879-6888. 10.1021/ac900874m.PubMedView ArticleGoogle Scholar
- Pellis L, Franssen-van Hal NLW, Burema J, Keijer J: The interclass correlation coefficient applied for evaluation of data correction, labeling methods, and rectal biopsy sampling in DNA microarray experiments. Physiol Genomics. 2003, 16: 99-106. 10.1152/physiolgenomics.00111.2003.PubMedView ArticleGoogle Scholar
- Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci. 1998, 95: 14863-14868. 10.1073/pnas.95.25.14863.PubMed CentralPubMedView ArticleGoogle Scholar
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gilette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: A knowledge-based approach for interpreting genome wide expression profiles. PNAS. 2005, 102: 15545-15550. 10.1073/pnas.0506580102.PubMed CentralPubMedView ArticleGoogle Scholar
- Coecke S, Vanhaecke T, Foriers A, Phillips IR, Vercruysse A, Shepard EA, Rogiers V: Hormonal regulation of glutathione S-transferase expression in co-cultured adult rat hepatocytes. J Endocrinol. 2000, 166: 363-371. 10.1677/joe.0.1660363.PubMedView ArticleGoogle Scholar
- Ikeda H, Serria MS, Kakizaki I, Hatayama I, Satoh K, Tsuchida S, Muramatsu M, Nishi S, Sakai M: Activation of mouse Pi-class glutathione S-transferase gene by Nrf2 (NF-E2-related factor 2) and androgen. Biochem J. 2002, 364: 563-570. 10.1042/BJ20011756.PubMed CentralPubMedView ArticleGoogle Scholar
- Depreter M, Vandesompele J, Espeel M, Speleman F, Roels F: Modulation of the peroxisomal gene expression pattern by dehydroepiandrosterone and vitamin D: therapeutic implications. J Endocrinol. 2002, 175: 779-792. 10.1677/joe.0.1750779.PubMedView ArticleGoogle Scholar
- Hazeldine J, Arlt W, Lord JM: Dehydroepiandrosterone as a regulator of immune cell function. J Steroid Biochem Mol Biol.Google Scholar
- Gu S, Ripp SL, Prough RA, Geoghegan TE: Dehydroepiandrosterone affects the expression of multiple genes in rat liver including 11β-Hydroxysteroid Dehydrogenase Type 1: A cDNA Array Analysis. Mol Pharmacol. 2003, 63: 722-731. 10.1124/mol.63.3.722.PubMedView ArticleGoogle Scholar
- Riedmaier I, Tichopad A, Reiter M, Pfaffl M, Meyer HHD: Identification of potential gene expression biomarkers for the surveillance of anabolic agents in bovine blood cells. Anal Chim Acta. 2009, 638: 106-113. 10.1016/j.aca.2009.02.014.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.