Optimising the analysis of transcript data using high density oligonucleotide arrays and genomic DNA-based probe selection

Background: Affymetrix GeneChip arrays are widely used for transcriptomic studies in a diverse range of species. Each gene is represented on a GeneChip array by a probe- set, consisting of up to 16 probe-pairs. Signal intensities across probe- pairs within a probe-set vary in part due to different physical hybridisation characteristics of individual probes with their target labelled transcripts. We have previously developed a technique to study the transcriptomes of heterologous species based on hybridising genomic DNA (gDNA) to a GeneChip array designed for a different species, and subsequently using only those probes with good homology. Results: Here we have investigated the effects of hybridising homologous species gDNA to study the transcriptomes of species for which the arrays have been designed. Genomic DNA from Arabidopsis thaliana and rice (Oryza sativa) were hybridised to the Affymetrix Arabidopsis ATH1 and Rice Genome GeneChip arrays respectively. Probe selection based on gDNA hybridisation intensity increased the number of genes identified as significantly differentially expressed in two published studies of Arabidopsis development, and optimised the analysis of technical replicates obtained from pooled samples of RNA from rice. Conclusion: This mixed physical and bioinformatics approach can be used to optimise estimates of gene expression when using GeneChip arrays.

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.

Cardiac protein kinases: the cardiomyocyte kinome and differential kinase expression in human failing hearts

Aims. Protein kinases are potential therapeutic targets for heart failure, but most studies of cardiac protein kinases derive from other systems, an approach that fails to account for specific kinases expressed in the heart and the contractile cardiomyocytes. We aimed to define the cardiomyocyte kinome (i.e. the protein kinases expressed in cardiomyocytes) and identify kinases with altered expression in human failing hearts. Methods and Results. Expression profiling (Affymetrix microarrays) detected >400 protein kinase mRNAs in rat neonatal ventricular myocytes (NVMs) and/or adult ventricular myocytes (AVMs), 32 and 93 of which were significantly upregulated or downregulated (>2-fold), respectively, in AVMs. Data for AGC family members were validated by qPCR. Proteomics analysis identified >180 cardiomyocyte protein kinases, with high relative expression of mitogen-activated protein kinase cascades and other known cardiomyocyte kinases (e.g. CAMKs, cAMP-dependent protein kinase). Other kinases are poorly-investigated (e.g. Slk, Stk24, Oxsr1). Expression of Akt1/2/3, BRaf, ERK1/2, Map2k1, Map3k8, Map4k4, MST1/3, p38-MAPK, PKCδ, Pkn2, Ripk1/2, Tnni3k and Zak was confirmed by immunoblotting. Relative to total protein, Map3k8 and Tnni3k were upregulated in AVMs vs NVMs. Microarray data for human hearts demonstrated variation in kinome expression that may influence responses to kinase inhibitor therapies. Furthermore, some kinases were upregulated (e.g. NRK, JAK2, STK38L) or downregulated (e.g. MAP2K1, IRAK1, STK40) in human failing hearts. Conclusions. This characterization of the spectrum of kinases expressed in cardiomyocytes and the heart (cardiomyocyte and cardiac kinomes) identified novel kinases, some of which are differentially expressed in failing human hearts and could serve as potential therapeutic targets.

Expression of the receptor tyrosine kinase ephb2 on dendritic cells is modulated by toll-like receptor ligation but is not required for t cell activation

The Eph receptor tyrosine kinases interact with their ephrin ligands on adjacent cells to facilitate contact-dependent cell communication. Ephrin B ligands are expressed on T cells and have been suggested to act as co-stimulatory molecules during T cell activation. There are no detailed reports of the expression and modulation of EphB receptors on dendritic cells, the main antigen presenting cells that interact with T cells. Here we show that mouse splenic dendritic cells (DC) and bone-marrow derived DCs (BMDC) express EphB2, a member of the EphB family. EphB2 expression is modulated by ligation of TLR4 and TLR9 and also by interaction with ephrin B ligands. Co-localization of EphB2 with MHC-II is also consistent with a potential role in T cell activation. However, BMDCs derived from EphB2 deficient mice were able to present antigen in the context of MHC-II and produce T cell activating cytokines to the same extent as intact DCs. Collectively our data suggest that EphB2 may contribute to DC responses, but that EphB2 is not required for T cell activation. This result may have arisen because DCs express other members of the EphB receptor family, EphB3, EphB4 and EphB6, all of which can interact with ephrin B ligands, or because EphB2 may be playing a role in another aspect of DC biology such as migration.

Regulation of Bcl-xL expression by H2O2 in cardiac myocytes.

Oxidative stress promotes cardiac myocyte apoptosis through the mitochondrial death pathway. Since Bcl-2 family proteins are key regulators of apoptosis, we examined the effects of H2O2 on the expression of principal Bcl-2 family proteins (Bcl-2, Bcl-xL, Bax, Bad) in neonatal rat cardiac myocytes. Protein expression was assessed by immunoblotting. Bcl-2, Bax, and Bad were all down-regulated in myocytes exposed to 0.2 mm H2O2, a concentration that induces apoptosis. In contrast, although Bcl-xL levels initially declined, the protein was re-expressed from 4-6 h. Bcl-xL mRNA was up-regulated from 2 to 4 h in neonatal rat or mouse cardiac myocytes exposed to H2O2, consistent with the re-expression of protein. Four different untranslated first exons have been identified for the Bcl-x gene (exons 1, 1B, 1C, and 1D, where exon 1 is the most proximal and exon 1D the most distal to the coding region). All were detected in mouse or rat neonatal cardiac myocytes, but exon 1D was not expressed in adult mouse hearts. In neonatal mouse or rat cardiac myocytes, H2O2 induced the expression of exons 1B, 1C, and 1D, but not exon 1. These data demonstrate that the Bcl-x gene is selectively responsive to oxidative stress, and the response is mediated through distal promoter regions.

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.

Retinoid-Related Receptor (ROR) alpha mRNA Expression Is Altered in the Brain of Male Mice Lacking All Ligand-Binding Thyroid Hormone Receptor (TR) Isoforms

In the vertebrate brain, the thalamus serves as a relay and integration station for diverse neuronal information en route from the periphery to the cortex. Deficiency of TH during development results in severe cerebral abnormalities similar to those seen in the mouse when the retinoic acid receptor (ROR)α gene is disrupted. To investigate the effect of the thyroid hormone recep-tors (TRs) on RORalpha gene expression, we used intact male mice, in which the genes encoding the α and beta TRs have been deleted. In situ hybridization for RORalpha mRNA revealed that this gene is expressed in specific areas of the brain including the thalamus, pons, cerebellum, cortex, and hippocampus. Our quantitative data showed differences in RORalpha mRNA expression in different subthalamic nuclei between wild-type and knock-out mice. For example, the centromedial nucleus of the thalamus, which plays a role in mediating nociceptive and visceral information from the brainstem to the basal ganglia and cortical regions, has less expression of RORalpha mRNA in the knockout mice (-37%) compared to the wild-type controls. Also, in the dorsal geniculate (+72%) and lateral posterior nuclei (+58%) we found more RORalpha mRNA in dKO as compared to dWT animals. Such differences in RORalpha mRNA expression may play a role in the behavioral alterations resulting from congenital hypothyroidism.