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

Affiliation: Institut de recherche en immunologie et en cancérologie, Université de Montréal

Regulation of gene expression in the dinoflagellate Lingulodinium polyedrum

Les dinoflagellés sont des eucaryotes unicellulaires que l’on retrouve autant en eau douce qu’en milieu marin. Ils sont particulièrement connus pour causer des fleurs d’algues toxiques nommées ‘marée-rouge’, ainsi que pour leur symbiose avec les coraux et pour leur importante contribution à la fixation du carbone dans les océans. Au point de vue moléculaire, ils sont aussi connus pour leur caractéristiques nucléaires uniques, car on retrouve généralement une quantité immense d’ADN dans leurs chromosomes et ceux-ci sont empaquetés et condensés sous une forme cristalline liquide au lieu de nucléosomes. Les gènes encodés par le noyau sont souvent présents en multiples copies et arrangés en tandem et aucun élément de régulation transcriptionnelle, y compris la boite TATA, n’a encore été observé. L’organisation unique de la chromatine des dinoflagellés suggère que différentes stratégies sont nécessaires pour contrôler l’expression des gènes de ces organismes. Dans cette étude, j’ai abordé ce problème en utilisant le dinoflagellé photosynthétique Lingulodinium polyedrum comme modèle. L. polyedrum est d’un intérêt particulier, car il a plusieurs rythmes circadiens (journalier). À ce jour, toutes les études sur l’expression des gènes lors des changements circadiens ont démontrées une régulation à un niveau traductionnel. Pour mes recherches, j’ai utilisé les approches transcriptomique, protéomique et phosphoprotéomique ainsi que des études biochimiques pour donner un aperçu de la mécanique de la régulation des gènes des dinoflagellés, ceci en mettant l’accent sur l’importance de la phosphorylation du système circadien de L. polyedrum. L’absence des protéines histones et des nucléosomes est une particularité des dinoflagellés. En utilisant la technologie RNA-Seq, j’ai trouvé des séquences complètes encodant des histones et des enzymes modifiant les histones. L polyedrum exprime donc des séquences conservées codantes pour les histones, mais le niveau d’expression protéique est plus faible que les limites de détection par immunodétection de type Western. Les données de séquençage RNA-Seq ont également été utilisées pour générer un transcriptome, qui est une liste des gènes exprimés par L. polyedrum. Une recherche par homologie de séquences a d’abord été effectuée pour classifier les transcrits en diverses catégories (Gene Ontology; GO). Cette analyse a révélé une faible abondance des facteurs de transcription et une surprenante prédominance, parmi ceux-ci, des séquences à domaine Cold Shock. Chez L. polyedrum, plusieurs gènes sont répétés en tandem. Un alignement des séquences obtenues par RNA-Seq avec les copies génomiques de gènes organisés en tandem a été réalisé pour examiner la présence de transcrits polycistroniques, une hypothèse formulée pour expliquer le manque d’élément promoteur dans la région intergénique de la séquence de ces gènes. Cette analyse a également démontré une très haute conservation des séquences codantes des gènes organisés en tandem. Le transcriptome a également été utilisé pour aider à l’identification de protéines après leur séquençage par spectrométrie de masse, et une fraction enrichie en phosphoprotéines a été déterminée comme particulièrement bien adapté aux approches d’analyse à haut débit. La comparaison des phosphoprotéomes provenant de deux périodes différentes de la journée a révélée qu’une grande partie des protéines pour lesquelles l’état de phosphorylation varie avec le temps est reliées aux catégories de liaison à l’ARN et de la traduction. Le transcriptome a aussi été utilisé pour définir le spectre des kinases présentes chez L. polyedrum, qui a ensuite été utilisé pour classifier les différents peptides phosphorylés qui sont potentiellement les cibles de ces kinases. Plusieurs peptides identifiés comme étant phosphorylés par la Casein Kinase 2 (CK2), une kinase connue pour être impliquée dans l’horloge circadienne des eucaryotes, proviennent de diverses protéines de liaison à l’ARN. Pour évaluer la possibilité que quelques-unes des multiples protéines à domaine Cold Shock identifiées dans le transcriptome puissent moduler l’expression des gènes de L. polyedrum, tel qu’observé chez plusieurs autres systèmes procaryotiques et eucaryotiques, la réponse des cellules à des températures froides a été examinée. Les températures froides ont permis d’induire rapidement un enkystement, condition dans laquelle ces cellules deviennent métaboliquement inactives afin de résister aux conditions environnementales défavorables. Les changements dans le profil des phosphoprotéines seraient le facteur majeur causant la formation de kystes. Les phosphosites prédits pour être phosphorylés par la CK2 sont la classe la plus fortement réduite dans les kystes, une découverte intéressante, car le rythme de la bioluminescence confirme que l’horloge a été arrêtée dans le kyste.

Role of histone methylation in the regulation of COX-2, iNOS, and mPGES-1 gene expression in human chondrocytes: Implication for Osteoarthritis

L'arthrose (OA) est une maladie articulaire dégénérative, classée comme la forme la plus fréquente au monde. Elle est caractérisée par la dégénérescence du cartilage articulaire, l’inflammation de la membrane synoviale, et le remodelage de l’os sous-chondral. Ces changements structurels et fonctionnels sont dues à de nombreux facteurs. Les cytokines, les prostaglandines (PG), et les espèces réactives de l'oxygène sont les principaux médiateurs impliqués dans la pathophysiologie de l'OA. L'interleukine-1β (IL-1β) est une cytokine pro-inflammatoire majeure qui joue un rôle crucial dans l'OA. L'IL-1β induit l'expression de la cyclooxygénase-2 (COX-2), la microsomale prostaglandine E synthase-1 (mPGES-1), la synthase inductible de l'oxyde nitrique (iNOS), ainsi que leurs produits la prostaglandine E2 (PGE2) et l'oxyde nitrique (NO). Ce sont des médiateurs essentiels de la réponse inflammatoire au cours de l'OA qui contribuent aux mécanismes des douleurs, de gonflement, et de destruction des tissus articulaires. Les modifications épigénétiques jouent un rôle très important dans la régulation de l’expression de ces gènes pro-inflammatoires. Parmi ces modifications, la méthylation/ déméthylation des histones joue un rôle critique dans la régulation des gènes. La méthylation/ déméthylation des histones est médiée par deux types d'enzymes: les histones méthyltransférases (HMT) et les histones déméthylases (HDM) qui favorisent l’activation et/ou la répression de la transcription. Il est donc nécessaire de comprendre les mécanismes moléculaires qui contrôlent l’expression des gènes de la COX-2, la mPGES-1, et l’iNOS. L'objectif de cette étude est de déterminer si la méthylation/déméthylation des histones contribute à la régulation de l’expression des gènes COX-2, mPGES-1, et iNOS dans des chondrocytes OA humains induits par l'IL-1β. Nous avons montré que la méthylation de la lysine K4 de l'histone H3 (H3K4) par SET-1A contribue à l’activation des gènes COX-2 et iNOS dans les chondrocytes humains OA induite par l'IL-1β. Nous avons également montré que la lysine K9 de l’histone H3 (H3K9) est déméthylée par LSD1, et que cette déméthylation contribue à l’expression de la mPGES-1 induite par IL-1β dans les chondrocytes humains OA. Nous avons aussi trouvé que les niveaux d'expression des enzymes SET-1A et LSD1 sont élevés au niveau du cartilage OA. Nos résultats montrent, pour la première fois, l'implication de la méthylation/ déméthylation des histones dans la régulation de l’expression des gènes COX-2, mPGES-1, et iNOS. Ces données suggèrent que ces mécanismes pourraient être une cible potentielle pour une intervention pharmacologique dans le traitement de la physiopathologie de l'OA.

GABAA and GABAB Receptor Gene Expression and Functional Regulation During Pancreatic Regeneration and Insulin Secretion in Rats

The present study demonstrate the functional alterations of the GABAA and GABAB receptors and the gene expression during the regeneration of pancreas following partial pancreatectomy. The role of these receptors in insulin secretion and pancreatic DNA synthesis using the specific agonists and antagonists also are studied in vitro. The alterations of GABAA and GABAR receptor function and gene expression in the brain stem, crebellum and hypothalamus play an important role in the sympathetic regulation of insulin secretion during pancreatic regeneration. Previous studies have given much information linking functional interaction between GABA and the peripheral nervous system. The involvement of specific receptor subtypes functional regulation during pancreatic regeneration has not given emphasis and research in this area seems to be scarce. We have observed a decreased GABA content, down regulation of GABAA receptors and an up regulation of GABAB receptors in the cerebral cortex, brain stem and hypothalamus. Real Time-PCR analysis confirmed the receptor data in the brain regions. These alterations in the GABAA and GABAB receptors of the brain are suggested to govern the regenerative response and growth regulation of the pancreas through sympathetic innervation. In addition, receptor binding studies and Real Time-PCR analysis revealed that during pancreatic regeneration GABAA receptors were down regulated and GABAB receptors were up regulated in pancreatic islets. This suggests an inhibitory role for GABAA receptors in islet cell proliferation i.e., the down regulation of this receptor facilitates proliferation. Insulin secretion study during 1 hour showed GABA has inhibited the insulin secretion in a dose dependent manner in normal and hyperglycaemic conditions. Bicuculline did not antagonize this effect. GABAA agonist, muscimol inhibited glucose stimulated insulin secretion from pancreatic islets except in the lowest concentration of 1O-9M in presence of 4mM glucose.Musclmol enhanced insulin secretion at 10-7 and 10-4M muscimol in presence of 20mM glucose- 4mM glucose represents normal and 20mM represent hyperglycaemic conditions. GABAB agonist, baclofen also inhibited glucose induced insulin secretion and enhanced at the concentration of 1O-5M at 4mM glucose and at 10-9M baclofen in presence of 20mM glucose. This shows a differential control of the GABAA and GABAB receptors over insulin release from the pancreatic islets. During 24 hours in vitro insulin secretion study it showed that low concentration of GABA has inhibited glucose stimulated insulin secretion from pancreatic islets. Muscimol, the GABAA agonist, inhibited the insulin secretion but, gave an enhanced secretion of insulin in presence of 4mM glucose at 10-7 , 10-5 and 1O-4M muscimol. But in presence of 20mM glucose muscimol significantly inhibited the insulin secretion. GABAB agonist, baclofen also inhibited glucose induced insulin secretion in presence of both 4mM and 20mM glucose. This shows the inhibitory role of GABA and its specific receptor subtypes over insulin synthesis from pancreatic bete-islets. In vitro DNA synthesis studies showed that activation of GABAA receptor by adding muscimol, a specific agonist, inhibited islet DNA synthesis. Also, the addition of baclofen, a specific agonist of GABAB receptor resulted in the stimulation of DNA synthesis.Thus the brain and pancreatic GABAA and GABAB receptor gene expression differentially regulates pancreatic insulin secretion and islet cell proliferation during pancreatic regeneration. This will have immense clinical significance in therapeutic applications in the management of Diabetes mellitus.

Juvenile-mature wood transition in pine: correlation between wood properties and candidate gene expression profiles

There is a strong desire to exploit transcriptomics data from model species for the genetic improvement of non-model crops. Here, we use gene expression profiles from the commercial model Pinus taeda to identify candidate genes implicated in juvenile-mature wood transition in the non-model relative, P. sylvestris. Re-analysis of 'public domain' SAGE data from xylem tissues of P. taeda revealed 283 mature-abundant and 396 juvenile-abundant tags (P < 0.01), of which 70 and 137, respectively matched to genes with known function. Based on sequence similarity, we then isolated 16 putative homologues of genes that in P. taeda exhibited widest divergence in expression between juvenile and mature samples. Candidate expression levels in P. sylvestris were almost invariably differential between juvenile and mature woody tissue samples among two cohorts of five trees collected from the same seed source and selected for genetic uniformity by genetic distance analysis. However, the direction of differential expression was not always consistent with that described in the original P. taeda SAGE data. Correlation was observed between gene expression and juvenile-mature wood anatomical characteristics by OPLS analysis. Four candidates (alpha-tubulin, porin MIP1, lipid transfer protein and aquaporin like protein) apparently had greatest influence on the wood traits measured. Speculative function of these genes in relation to juvenile-mature wood transition is briefly explored. Thus, we demonstrate the feasibility of exploiting SAGE data from a model species to identify consistently differentially expressed candidates in a related non-model species.

Dietary alpha-tocopherol affects differential gene expression in rat testes

Gene-chip technology was employed to study the effect of dietary vitamin E (VE) on gene expression in rat testes. Male albino rats were fed with either a diet deficient in VE or a standard diet containing VE. Differential gene expression was monitored at five individual time-points over a period of 14 months with all animals individually pro. led. Low VE intake resulted in the consistent upregulation of 7-dehydrocholesterol reductase and GATA binding protein 4, both involved in testosterone synthesis. Cyclin D3, important in cell cycle progression and Wilms tumor 1, related to cancer development, were also up-regulated in the vitamin E deficient animals. This study demonstrates that low dietary VE intake has long-term effects on gene expression in the testes. Our data provides insights into the possible molecular mechanisms underlying the beneficial effects of vitamin E on the male reproductive organ.

Identification of hepatic molecular mechanisms of action of alpha-tocopherol using global gene expression profile analysis in rats

The recent discovery that vitamin E (VE) regulates gene activity at the transcriptional level indicates that VE may exert part of its biological effects by mechanisms which may be independent of its well-recognised antioxidant function. The objective of this study was the identification of hepatic vitamin E-sensitive genes and examination of the effects of VE on their corresponding biological endpoints. Two groups of male rats were randomly assigned to either a VE-sufficient diet or to a control diet deficient in VE for 290 days. High-density oligonucleotide microarrays comprising over 7000 genes were used to assess the transcriptional response of the liver. Differential gene expression was monitored over a period of 9 months, at four different time-points, and rats were individually profiled. This experimental strategy identified several VE-sensitive genes, which were chronically altered by dietary VE. VE supplementation down-regulated scavenger receptor CD36, coagulation factor IX and 5-alpha-steroid reductase type 1 mRNA levels while hepatic gamma glutamyl-cysteinyl synthetase was significantly up-regulated. Measurement of the corresponding biological endpoints such as activated partial thromboplastin time, plasma dihydrotestosterone and hepatic glutathione substantiated the gene chip data which indicated that dietary VE plays an important role in a range of metabolic processes within the liver. (C) 2004 Elsevier B.V. All rights reserved.

Gene expression profiling of human hibernating myocardium: Increased expression of B-type natriuretic peptide and proenkephalin in hypocontractile vs normally-contracting regions of the heart

A greater understanding of the molecular basis of hibernating myocardium may assist in identifying those patients who would most benefit from revascularization. Paired heart biopsies were taken from hypocontractile and normally-contracting myocardium (identified by cardiovascular magnetic resonance) from 6 patients with chronic stable angina scheduled for bypass grafting. Gene expression profiles of hypocontractile and normally-contracting samples were compared using Affymetrix microarrays. The data for patients with confirmed hibernating myocardium were analysed separately and a different, though overlapping, set (up to 380) of genes was identified which may constitute a molecular fingerprint for hibernating myocardium. The expression of B-type natriuretic peptide (BNP) was increased in hypocontractile relative to normally-contracting myocardium. The expression of BNP correlated most closely with the expression of proenkephalin and follistatin 3, which may constitute additional heart failure markers. Our data illustrate differential gene expression in hypocontractile and/hibernating myocardium relative to normally-contracting myocardium within individual human hearts. Changes in expression of these genes, including increased relative expression of natriuretic and other factors, may constitute a molecular signature for hypocontractile and/or hibernating myocardium.

Pathogenic LRRK2 mutations do not alter gene expression in cell model systems or human brain tissue

Point mutations in LRRK2 cause autosomal dominant Parkinson's disease. Despite extensive efforts to determine the mechanism of cell death in patients with LRRK2 mutations, the aetiology of LRRK2 PD is not well understood. To examine possible alterations in gene expression linked to the presence of LRRK2 mutations, we carried out a case versus control analysis of global gene expression in three systems: fibroblasts isolated from LRRK2 mutation carriers and healthy, non-mutation carrying controls; brain tissue from G2019S mutation carriers and controls; and HEK293 inducible LRRK2 wild type and mutant cell lines. No significant alteration in gene expression was found in these systems following correction for multiple testing. These data suggest that any alterations in basal gene expression in fibroblasts or cell lines containing mutations in LRRK2 are likely to be quantitatively small. This work suggests that LRRK2 is unlikely to play a direct role in modulation of gene expression, although it remains possible that this protein can influence mRNA expression under pathogenic cicumstances.

Ethylene and 1-methylcyclopropene differentially regulate gene expression during onion sprout suppression

Onion (Allium cepa) is regarded as a nonclimacteric vegetable. In onions, however, ethylene can suppress sprouting while the ethylene-binding inhibitor 1-methylcyclopropene (1-MCP) can also suppress sprout growth; yet, it is unknown how ethylene and 1-MCP elicit the same response. In this study, onions were treated with 10 mu L L(-1) ethylene or 1 mu L L(-1) 1-MCP individually or in combination for 24 h at 20 degrees C before or after curing (6 weeks) at 20 degrees C or 28 degrees C and then stored at 1 degrees C. Following curing, a subset of these same onions was stored separately under continuous air or ethylene (10 mu L L(-1)) at 1 degrees C. Onions treated with ethylene and 1-MCP in combination after curing for 24 h had reduced sprout growth as compared with the control 25 weeks after harvest. Sprout growth following storage beyond 25 weeks was only reduced through continuous ethylene treatment. This observation was supported by a higher proportion of down-regulated genes characterized as being involved in photosynthesis, measured using a newly developed onion microarray. Physiological and biochemical data suggested that ethylene was being perceived in the presence of 1-MCP, since sprout growth was reduced in onions treated with 1-MCP and ethylene applied in combination but not when applied individually. A cluster of probes representing transcripts up-regulated by 1-MCP alone but down-regulated by ethylene alone or in the presence of 1-MCP support this suggestion. Ethylene and 1-MCP both down-regulated a probe tentatively annotated as an ethylene receptor as well as ethylene-insensitive 3, suggesting that both treatments down-regulate the perception and signaling events of ethylene.

Gene expression changes in phosphorus deficient potato (Solanum tuberosum L.) leaves and the potential for diagnostic gene expression markers

Background: There are compelling economic and environmental reasons to reduce our reliance on inorganic phosphate (Pi) fertilisers. Better management of Pi fertiliser applications is one option to improve the efficiency of Pi fertiliser use, whilst maintaining crop yields. Application rates of Pi fertilisers are traditionally determined from analyses of soil or plant tissues. Alternatively, diagnostic genes with altered expression under Pi limiting conditions that suggest a physiological requirement for Pi fertilisation, could be used to manage Pifertiliser applications, and might be more precise than indirect measurements of soil or tissue samples. Results: We grew potato (Solanum tuberosum L.) plants hydroponically, under glasshouse conditions, to control their nutrient status accurately. Samples of total leaf RNA taken periodically after Pi was removed from the nutrient solution were labelled and hybridised to potato oligonucleotide arrays. A total of 1,659 genes were significantly differentially expressed following Pi withdrawal. These included genes that encode proteins involved in lipid, protein, and carbohydrate metabolism, characteristic of Pi deficient leaves and included potential novel roles for genes encoding patatin like proteins in potatoes. The array data were analysed using a support vector machine algorithm to identify groups of genes that could predict the Pi status of the crop. These groups of diagnostic genes were tested using field grown potatoes that had either been fertilised or unfertilised. A group of 200 genes could correctly predict the Pi status of field grown potatoes. Conclusions: This paper provides a proof-of-concept demonstration for using microarrays and class prediction tools to predict the Pi status of a field grown potato crop. There is potential to develop this technology for other biotic and abiotic stresses in field grown crops. Ultimately, a better understanding of crop stresses may improve our management of the crop, improving the sustainability of agriculture.

A Brassica exon array for whole-transcript gene expression profiling

Affymetrix GeneChip (R) arrays are used widely to study transcriptional changes in response to developmental and environmental stimuli. GeneChip (R) arrays comprise multiple 25-mer oligonucleotide probes per gene and retain certain advantages over direct sequencing. For plants, there are several public GeneChip (R) arrays whose probes are localised primarily in 39 exons. Plant whole-transcript (WT) GeneChip (R) arrays are not yet publicly available, although WT resolution is needed to study complex crop genomes such as Brassica, which are typified by segmental duplications containing paralogous genes and/or allopolyploidy. Available sequence data were sampled from the Brassica A and C genomes, and 142,997 gene models identified. The assembled gene models were then used to establish a comprehensive public WT exon array for transcriptomics studies. The Affymetrix GeneChip (R) Brassica Exon 1.0 ST Array is a 5 mu M feature size array, containing 2.4 million 25-base oligonucleotide probes representing 135,201 gene models, with 15 probes per gene distributed among exons. Discrimination of the gene models was based on an E-value cut-off of 1E(-5), with <= 98 sequence identity. The 135 k Brassica Exon Array was validated by quantifying transcriptome differences between leaf and root tissue from a reference Brassica rapa line (R-o-18), and categorisation by Gene Ontologies (GO) based on gene orthology with Arabidopsis thaliana. Technical validation involved comparison of the exon array with a 60-mer array platform using the same starting RNA samples. The 135 k Brassica Exon Array is a robust platform. All data relating to the array design and probe identities are available in the public domain and are curated within the BrassEnsembl genome viewer at http://www.brassica.info/BrassEnsembl/index.html.

Rescue of gene expression by modified REST decoy oligonucleotides in a cellular model of Huntington's disease

Transcriptional dysfunction is a prominent hallmark of Huntington's disease (HD). Several transcription factors have been implicated in the aetiology of HD progression and one of the most prominent is repressor element 1 (RE1) silencing transcription factor (REST). REST is a global repressor of neuronal gene expression and in the presence of mutant Huntingtin increased nuclear REST levels lead to elevated RE1 occupancy and a concomitant increase in target gene repression, including brain-derived neurotrophic factor. It is of great interest to devise strategies to reverse transcriptional dysregulation caused by increased nuclear REST and determine the consequences in HD. Thus far, such strategies have involved RNAi or mutant REST constructs. Decoys are double-stranded oligodeoxynucleotides corresponding to the DNA-binding element of a transcription factor and act to sequester it, thereby abrogating its transcriptional activity. Here, we report the use of a novel decoy strategy to rescue REST target gene expression in a cellular model of HD. We show that delivery of the decoy in cells expressing mutant Huntingtin leads to its specific interaction with REST, a reduction in REST occupancy of RE1s and rescue of target gene expression, including Bdnf. These data point to an alternative strategy for rebalancing the transcriptional dysregulation in HD.

Using U0126 to dissect the role of the extracellular signal-regulated kinase 1/2 (ERK1/2) cascade in the regulation of gene expression by endothelin-1 in cardiac myocytes

The hypertrophic agonist endothelin-1 rapidly but transiently activates the extracellular signal-regulated kinase 1/2 (ERK1/2) cascade (and other signalling pathways) in cardiac myocytes, but the events linking this to hypertrophy are not understood. Using Affymetrix rat U34A microarrays, we identified the short-term (2-4 h) changes in gene expression induced in neonatal myocytes by endothelin-1 alone or in combination with the ERK1/2 cascade inhibitor, U0126. Expression of 15 genes was significantly changed by U0126 alone, and expression of an additional 78 genes was significantly changed by endothelin-1. Of the genes upregulated by U0126, four are classically induced through the aryl hydrocarbon receptor (AhR) by dioxins suggesting that U0126 activates the xenobiotic response element in cardiac myocytes potentially independently of effects on ERK1/2 signalling. The 78 genes showing altered expression with endothelin-1 formed five clusters: (i) three clusters showing upregulation by endothelin-1 according to time course (4 h > 2 h; 2 h > 4 h; 2 h approximately 4 h) with at least partial inhibition by U0126; (ii) a cluster of 11 genes upregulated by endothelin-1 but unaffected by U0126 suggesting regulation through signalling pathways other than ERK1/2; (iii) a cluster of six genes downregulated by endothelin-1 with attenuation by U0126. Thus, U0126 apparently activates the AhR in cardiac myocytes (which must be taken into account in protracted studies), but careful analysis allows identification of genes potentially regulated acutely via the ERK1/2 cascade. Our data suggest that the majority of changes in gene expression induced by endothelin-1 are mediated by the ERK1/2 cascade.

Chemotherapeutic agents and gene expression in cardiac myocytes

It is becoming apparent that anti-cancer chemotherapies are increasingly associated with cardiac dysfunction or even congestive heart failure (Minotti et al., 2004; Eliott, 2006; Suter et al., 2004; Ren, 2005). Our data suggest that one of the contributing factors to the cardiotoxicitiy of these drugs may be the activation of the AhR-response (including the increased expression of Cyp1a1) and/or other detoxification program in cardiac myocytes themselves. The induction of such responses may have secondary effects (e.g. to increase the level of intracellular oxidative stress), which may influence the contractility or even survival of cardiac myocytes. Furthermore, the specific response of cardiac myocytes, both with respect to the metabolizing enzymes and the export channels, potentially differs from other cells (e.g. we failed to detect any increase in expression of other “classical” AhR-responsive genes, Ugt1a1 and Ugt1a6). This could account for, for example, the observation that doxoribicinol (the 13-hydroxy form of doxorubicin) accumulates in cardiac myocytes but not in hepatocytes (Del Tacca et al., 1985; Olson et al., 1988). Given the vulnerability of the heart and the almost irreparable damage that can be done by severe oxidative stress, further studies would seem to be merited specifically on the effects of chemotherapeutic agents on cardiac myocytes.