Conifer defence against insects: microarray gene expression profiling of Sitka spruce (Picea sitchensis) induced by mechanical wounding or feeding by spruce budworms (Choristoneura occidentalis) or white pine weevils (Pissodes strobi) reveals large

Conifers are resistant to attack from a large number of potential herbivores or pathogens. Previous molecular and biochemical characterization of selected conifer defence systems support a model of multigenic, constitutive and induced defences that act on invading insects via physical, chemical, biochemical or ecological (multitrophic) mechanisms. However, the genomic foundation of the complex defence and resistance mechanisms of conifers is largely unknown. As part of a genomics strategy to characterize inducible defences and possible resistance mechanisms of conifers against insect herbivory, we developed a cDNA microarray building upon a new spruce (Picea spp.) expressed sequence tag resource. This first-generation spruce cDNA microarray contains 9720 cDNA elements representing c. 5500 unique genes. We used this array to monitor gene expression in Sitka spruce (Picea sitchensis) bark in response to herbivory by white pine weevils (Pissodes strobi, Curculionidae) or wounding, and in young shoot tips in response to western spruce budworm (Choristoneura occidentalis, Lepidopterae) feeding. Weevils are stem-boring insects that feed on phloem, while budworms are foliage feeding larvae that consume needles and young shoot tips. Both insect species and wounding treatment caused substantial changes of the host plant transcriptome detected in each case by differential gene expression of several thousand array elements at 1 or 2 d after the onset of treatment. Overall, there was considerable overlap among differentially expressed gene sets from these three stress treatments. Functional classification of the induced transcripts revealed genes with roles in general plant defence, octadecanoid and ethylene signalling, transport, secondary metabolism, and transcriptional regulation. Several genes involved in primary metabolic processes such as photosynthesis were down-regulated upon insect feeding or wounding, fitting with the concept of dynamic resource allocation in plant defence. Refined expression analysis using gene-specific primers and real-time PCR for selected transcripts was in agreement with microarray results for most genes tested. This study provides the first large-scale survey of insect-induced defence transcripts in a gymnosperm and provides a platform for functional investigation of plant-insect interactions in spruce. Induction of spruce genes of octadecanoid and ethylene signalling, terpenoid biosynthesis, and phenolic secondary metabolism are discussed in more detail.

Spatial analysis of biomineralization associated gene expression from the mantle organ of the pearl oyster Pinctada maxima

Background: Biomineralization is a process encompassing all mineral containing tissues produced within an organism. One of the most dynamic examples of this process is the formation of the mollusk shell, comprising a variety of crystal phases and microstructures. The organic component incorporated within the shell is said to dictate this architecture. However general understanding of how this process is achieved remains ambiguous. The mantle is a conserved organ involved in shell formation throughout molluscs. Specifically the mantle is thought to be responsible for secreting the protein component of the shell. This study employs molecular approaches to determine the spatial expression of genes within the mantle tissue to further the elucidation of the shell biomineralization. Results: A microarray platform was custom generated (PmaxArray 1.0) from the pearl oyster Pinctada maxima. PmaxArray 1.0 consists of 4992 expressed sequence tags (ESTs) originating from mantle tissue. This microarray was used to analyze the spatial expression of ESTs throughout the mantle organ. The mantle was dissected into five discrete regions and analyzed for differential gene expression with PmaxArray 1.0. Over 2000 ESTs were determined to be differentially expressed among the tissue sections, identifying five major expression regions. In situ hybridization validated and further localized the expression for a subset of these ESTs. Comparative sequence similarity analysis of these ESTs revealed a number of the transcripts were novel while others showed significant sequence similarities to previously characterized shell related genes.

A transcription factor map as revealed by a genome-wide gene expression analysis of whole-blood mRNA transcriptome in multiple sclerosis

Background Several lines of evidence suggests that transcription factors are involved in the pathogenesis of Multiple Sclerosis (MS) but a complete mapping the whole network has been elusive. One of the reasons is that there are several clinical subtypes of MS and transcription factors which may be involved in one subtype may not be in others. We investigated the possibility that this network could be mapped using microarray technologies and modern bioinformatics methods on a dataset from whole blood in 99 untreated MS patients (36 Relapse Remitting MS, 43 Primary Progressive MS, and 20 Secondary Progressive MS) and 45 age-matched healthy controls, Methodology/Principal Findings We have used two different analytical methodologies: a differential expression analysis and a differential co-expression analysis, which have converged on a significant number of regulatory motifs that seem to be statistically overrepresented in genes which are either differentially expressed (or differentially co-expressed) in cases and controls (e.g. V$KROX_Q6, p-value < 3.31E-6; V$CREBP1_Q2, p-value < 9.93E-6, V$YY1_02, p-value < 1.65E-5). Conclusions/significance: Our analysis uncovered a network of transcription factors that potentially dysregulate several genes in MS or one or more of its disease subtypes. Analysing the published literature we have found that these transcription factors are involved in the early T-lymphocyte specification and commitment as well as in oligodendrocytes dedifferentiation and development. The most significant transcription factors motifs were for the Early Growth response EGR/KROX family, ATF2, YY1 (Yin and Yang 1), E2F-1/DP-1 and E2F-4/DP-2 heterodimers, SOX5, and CREB and ATF families.

Gene expression analysis of diagnostic biopsies predicts pathological response to neoadjuvant chemoradiotherapy of esophageal cancer

OBJECTIVE: This study explored gene expression differences in predicting response to chemoradiotherapy in esophageal cancer. PURPOSE:: A major pathological response to neoadjuvant chemoradiation is observed in about 40% of esophageal cancer patients and is associated with favorable outcomes. However, patients with tumors of similar histology, differentiation, and stage can have vastly different responses to the same neoadjuvant therapy. This dichotomy may be due to differences in the molecular genetic environment of the tumor cells. BACKGROUND DATA: Diagnostic biopsies were obtained from a training cohort of esophageal cancer patients (13), and extracted RNA was hybridized to genome expression microarrays. The resulting gene expression data was verified by qRT-PCR. In a larger, independent validation cohort (27), we examined differential gene expression by qRT-PCR. The ability of differentially-regulated genes to predict response to therapy was assessed in a multivariate leave-one-out cross-validation model. RESULTS: Although 411 genes were differentially expressed between normal and tumor tissue, only 103 genes were altered between responder and non-responder tumor; and 67 genes differentially expressed >2-fold. These included genes previously reported in esophageal cancer and a number of novel genes. In the validation cohort, 8 of 12 selected genes were significantly different between the response groups. In the predictive model, 5 of 8 genes could predict response to therapy with 95% accuracy in a subset (74%) of patients. CONCLUSIONS: This study has identified a gene microarray pattern and a set of genes associated with response to neoadjuvant chemoradiation in esophageal cancer. The potential of these genes as biomarkers of response to treatment warrants further investigation. Copyright © 2009 by Lippincott Williams & Wilkins.

G9a histone methyltransferase contributes to imprinting in the mouse placenta

Whereas DNA methylation is essential for genomic imprinting, the importance of histone methylation in the allelic expression of imprinted genes is unclear. Imprinting control regions (ICRs), however, are marked by histone H3-K9 methylation on their DNA-methylated allele. In the placenta, the paternal silencing along the Kcnq1 domain on distal chromosome 7 also correlates with the presence of H3-K9 methylation, but imprinted repression at these genes is maintained independently of DNA methylation. To explore which histone methyltransferase (HMT) could mediate the allelic H3-K9 methylation on distal chromosome 7, and at ICRs, we generated mouse conceptuses deficient for the SET domain protein G9a. We found that in the embryo and placenta, the differential DNA methylation at ICRs and imprinted genes is maintained in the absence of G9a. Accordingly, in embryos, imprinted gene expression was unchanged at the domains analyzed, in spite of a global loss of H3-K9 dimethylation (H3K9me2). In contrast, the placenta-specific imprinting of genes on distal chromosome 7 is impaired in the absence of G9a, and this correlates with reduced levels of H3K9me2 and H3K9me3. These findings provide the first evidence for the involvement of an HMT and suggest that histone methylation contributes to imprinted gene repression in the trophoblast.

Gene expression profiling of cuticular proteins across the moult cycle of the crab Portunus pelagicus

Background: Crustaceans represent an attractive model to study biomineralization and cuticle matrix formation, as these events are precisely timed to occur at certain stages of the moult cycle. Moulting, the process by which crustaceans shed their exoskeleton, involves the partial breakdown of the old exoskeleton and the synthesis of a new cuticle. This cuticle is subdivided into layers, some of which become calcified while others remain uncalcified. The cuticle matrix consists of many different proteins that confer the physical properties, such as pliability, of the exoskeleton. Results: We have used a custom cDNA microarray chip, developed for the blue swimmer crab Portunus pelagicus, to generate expression profiles of genes involved in exoskeletal formation across the moult cycle. A total of 21 distinct moult-cycle related differentially expressed transcripts representing crustacean cuticular proteins were isolated. Of these, 13 contained copies of the cuticle_1 domain previously isolated from calcified regions of the crustacean exoskeleton, four transcripts contained a chitin_bind_4 domain (RR consensus sequence) associated with both the calcified and un-calcified cuticle of crustaceans, and four transcripts contained an unannotated domain (PfamB_109992) previously isolated from C. pagurus. Additionally, cryptocyanin, a hemolymph protein involved in cuticle synthesis and structural integrity, also displays differential expression related to the moult cycle. Moult stage-specific expression analysis of these transcripts revealed that differential gene expression occurs both among transcripts containing the same domain and among transcripts containing different domains. Conclusion: The large variety of genes associated with cuticle formation, and their differential expression across the crustacean moult cycle, point to the complexity of the processes associated with cuticle formation and hardening. This study provides a molecular entry path into the investigation of the gene networks associated with cuticle formation.

Gene expression profile of epithelial cells and mesenchymal cells derived from limbal explant culture

Purpose: Limbal stem cell deficiency is a challenging clinical problem and the current treatment involves replenishing the depleted limbal stem cell (LSC) pool by either limbal tissue transplantation or use of cultivated limbal epithelial cells (LEC). Our experience of cultivating the LEC on denuded human amniotic membrane using a feeder cell free method, led to identification of mesenchymal cells of limbus (MC-L), which showed phenotypic resemblance to bone marrow derived mesenchymal stem cells (MSC-BM). To understand the transcriptional profile of these cells, microarray experiments were carried out.Methods: RNA was isolated from cultured LEC, MC-L and MSC-BM and microarray experiments were carried out by using Agilent chip (4x44 k). The microarray data was validated by using Realtime and semiquntitative reverse transcription polymerase chain reaction. Results: The microarray analysis revealed specific gene signature of LEC and MC-L, and also their complementary role related to cytokine and growth factor profile, thus supporting the nurturing roles of the MC-L. We have also observed similar and differential gene expression between MC-L and MSC-BM.Conclusions: This study represents the first extensive gene expression analysis of limbal explant culture derived epithelial and mesenchymal cells and as such reveals new insight into the biology, ontogeny, and in vivo function of these cells.

Association of gene expression with sequential proliferation, differentiation and tumor formation in murine skin

Differential gene expression in two established initiation and promotion skin carcinogenesis models during promotion and tumor formation was determined by microarray technology with the purpose of distinguishing the genes more associated with neoplastic transformation from those linked with proliferation and differentiation. The first model utilized dimethylbenz[a]anthracene initiation and 12-O-tetradecanoylphorbol 13-acetate (TPA) promotion in the FVB/N mouse, and the second TPA promotion of the Tg.Ac mouse, which is endogenously initiated by virtue of an activated Ha-ras transgene. Comparison of gene expression profiles across the two models identified genes whose altered expression was associated with papilloma formation rather than TPA-induced proliferation and differentiation. DMBA suppressed TPA-induced differentiation which allowed identification of those genes associated more specifically with differentiation rather than proliferation. EASE (Expression Analysis Systemic Explorer) indicated a correlation between muscle-associated genes and skin differentiation, whereas genes involved with protein biosynthesis were strongly correlated with proliferation. For verification the altered expression of selected genes were confirmed by RT-PCR; Carbonic anhydrase 2, Thioredoxin 1 and Glutathione S-transferase omega 1 associated with papilloma formation and Enolase 3, Cystatin 6 and Filaggrin associated with TPA-induced proliferation and differentiation. In situ analysis located the papillomas Glutathione S-transferase omega 1 expression to the proliferating areas of the papillomas. Thus we have identified profiles of differential gene expression associated with the tumorigenesis and promotion stages for skin carcinogenesis in the mouse.

cudaMap: a GPU accelerated program for gene expression connectivity mapping

Background: Modern cancer research often involves large datasets and the use of sophisticated statistical techniques. Together these add a heavy computational load to the analysis, which is often coupled with issues surrounding data accessibility. Connectivity mapping is an advanced bioinformatic and computational technique dedicated to therapeutics discovery and drug re-purposing around differential gene expression analysis. On a normal desktop PC, it is common for the connectivity mapping task with a single gene signature to take >2h to complete using sscMap, a popular Java application that runs on standard CPUs (Central Processing Units). Here, we describe new software, cudaMap, which has been implemented using CUDA C/C++ to harness the computational power of NVIDIA GPUs (Graphics Processing Units) to greatly reduce processing times for connectivity mapping.

Results: cudaMap can identify candidate therapeutics from the same signature in just over thirty seconds when using an NVIDIA Tesla C2050 GPU. Results from the analysis of multiple gene signatures, which would previously have taken several days, can now be obtained in as little as 10 minutes, greatly facilitating candidate therapeutics discovery with high throughput. We are able to demonstrate dramatic speed differentials between GPU assisted performance and CPU executions as the computational load increases for high accuracy evaluation of statistical significance.

Conclusion: Emerging 'omics' technologies are constantly increasing the volume of data and information to be processed in all areas of biomedical research. Embracing the multicore functionality of GPUs represents a major avenue of local accelerated computing. cudaMap will make a strong contribution in the discovery of candidate therapeutics by enabling speedy execution of heavy duty connectivity mapping tasks, which are increasingly required in modern cancer research. cudaMap is open source and can be freely downloaded from http://purl.oclc.org/NET/cudaMap.

Gene expression in understanding mechanisms of toxicity in amphibians

De uma forma geral os anfíbios são conhecidos como organismos que apresentam uma grande sensibilidade a vários tipos de contaminantes. Contudo existem casos, como o de Pelophylax perezi (rã-verde), em que estes organismos habitam áreas extremamente contaminadas. Este facto verifica-se na mina de urânio desactivada, da Cunha Baixa (Viseu, centro de Portugal), em que uma população destas rãs habita na lagoa de efluente ácido mineiro (M). Estudos ecotoxicológicos anteriores com estes organismos revelaram apenas efeitos de toxicidade ténues levantando algumas questões. Com o objectivo de elucidar quais os mecanismos que permitem a P. perezi permanecer neste local, sem sofrer aparentemente efeitos perniciosos, encetamos este trabalho. Numa primeira abordagem, avaliámos o sistema de defesa antioxidante de rãs adultas, bem como o conteúdo em metais de alguns órgãos. Desta forma verificámos alterações enzimáticas, principalmente no pulmão e acumulação de metais nos vários órgãos. Posteriormente foi realizado um estudo de expressão genética diferencial, também em organismos adultos e desta feita foram sugeridos alguns mecanismos de protecção basal que estarão por detrás da capacidade de suportar este ambiente extremamente contaminado. Numa etapa seguinte abordámos os efeitos em fases larvares, fazendo inicialmente uma exposição in situ, a vários efluentes, caracteristicamente diferentes, do complexo mineiro. Avaliámos o crescimento, a acumulação de metais e a actividade de alguns biomarcadores de stress oxidativo. Como resultado pudemos constatar que nas fases larvares para além de alguma mortalidade existe acumulação de metais bem como algumas alterações a nível de biomarcadores de stress oxidativo. Numa última abordagem realizamos uma exposição crónica dos girinos a efluente da mina com diversos níveis de pH para distinguir os efeitos da toxicidade do pH, dos efeitos da toxicidade pelo conteúdo de metais. Para tal avaliámos novamente biomarcadores de stress oxidativo, crescimento, acumulação de metais e efectuamos ainda um estudo de expressão genética diferencial. Esta última aproximação permitiu verificar que a toxicidade do efluente resulta primariamente do pH ácido, assumindo a contaminação por metais um papel secundário. Contudo o crescimento dos girinos de P. perezi apresenta-se estimulado por pHs mais baixos. São apontados ainda alguns mecanismos, em girinos, para lidar com o stress causado pela contaminação por metais.De uma forma geral pôde-se constatar que quer anfíbios adultos quer girinos expostos ao efluente apresentam valores altos de metais acumulados. Os biomarcadores de stress oxidativo na sua maioria não apresentaram respostas coerentes mediante as várias exposições. Este trabalho apresentase como um contributo importante para a ecotoxicologia de anfíbios, aumentando os níveis actuais de conhecimento sobre o efeito de contaminação proveniente de efluentes mineiros, sugerindo ainda mecanismos de resistência quer em larvas, quer para adultos.

Comprehensive Expression Analyses of Neural Cell-Type-Specific miRNAs Identify New Determinants of the Specification and Maintenance of Neuronal Phenotypes.

MicroRNAs (miRNAs) have been shown to play important roles in both brain development and the regulation of adult neural cell functions. However, a systematic analysis of brain miRNA functions has been hindered by a lack of comprehensive information regarding the distribution of miRNAs in neuronal versus glial cells. To address this issue, we performed microarray analyses of miRNA expression in the four principal cell types of the CNS (neurons, astrocytes, oligodendrocytes, and microglia) using primary cultures from postnatal d 1 rat cortex. These analyses revealed that neural miRNA expression is highly cell-type specific, with 116 of the 351 miRNAs examined being differentially expressed fivefold or more across the four cell types. We also demonstrate that individual neuron-enriched or neuron-diminished RNAs had a significant impact on the specification of neuronal phenotype: overexpression of the neuron-enriched miRNAs miR-376a and miR-434 increased the differentiation of neural stem cells into neurons, whereas the opposite effect was observed for the glia-enriched miRNAs miR-223, miR-146a, miR-19, and miR-32. In addition, glia-enriched miRNAs were shown to inhibit aberrant glial expression of neuronal proteins and phenotypes, as exemplified by miR-146a, which inhibited neuroligin 1-dependent synaptogenesis. This study identifies new nervous system functions of specific miRNAs, reveals the global extent to which the brain may use differential miRNA expression to regulate neural cell-type-specific phenotypes, and provides an important data resource that defines the compartmentalization of brain miRNAs across different cell types.

Expression profiles in the progression of ductal carcinoma in the breast

An understanding of the multi-step nature of cancer as it is in the breast, as a series of pivotal genetic/epigenetic modifications is irrefutably a milestone in diagnostics, prognostics and eventually providing a cure. Here we have utilised a variant of analysis of variance (ANOVA) as a model for the identification and tracking of specific mRNA species whose transcription has been significantly altered at each grade in the progression of ductal carcinoma, making it possible to correlate histological progression with the genetic events underlying breast cancer. We show that in the progression of ductal carcinomas, from grade 1 to 3, there is a reduction in the actual number of mRNA species, which are significantly over or under expressed. We also show that this technique can be employed to generate differential gene expression patterns, whereby the combined expression profile of the tailored spectra of genes in the comparison of each ductal grade is sufficient to render them on clearly separate arms of an array-wise hierarchical cluster dendrogram.

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.

Changes in gene expression in Arabidopsis shoots during phosphate starvation and the potential for developing smart plants

Our aim was to generate and prove the concept of "smart" plants to monitor plant phosphorus (P) status in Arabidopsis. Smart plants can be genetically engineered by transformation with a construct containing the promoter of a gene up-regulated specifically by P starvation in an accessible tissue upstream of a marker gene such as beta-glucuronidase (GUS). First, using microarrays, we identified genes whose expression changed more than 2.5-fold in shoots of plants growing hydroponically when P, but not N or K, was withheld from the nutrient solution. The transient changes in gene expression occurring immediately (4 h) after P withdrawal were highly variable, and many nonspecific, shock-induced genes were up-regulated during this period. However, two common putative cis-regulatory elements (a PHO-like element and a TATA box-like element) were present significantly more often in the promoters of genes whose expression increased 4 h after the withdrawal of P compared with their general occurrence in the promoters of all genes represented on the microarray. Surprisingly, the expression of only four genes differed between shoots of P-starved and -replete plants 28 h after P was withdrawn. This lull in differential gene expression preceded the differential expression of a new group of 61 genes 100 h after withdrawing P. A literature survey indicated that the expression of many of these "late" genes responded specifically to P starvation. Shoots had reduced P after 100 h, but growth was unaffected. The expression of SQD1, a gene involved in the synthesis of sulfolipids, responded specifically to P starvation and was increased 100 h after withdrawing P. Leaves of Arabidopsis bearing a SQD1::GUS construct showed increased GUS activity after P withdrawal, which was detectable before P starvation limited growth. Hence, smart plants can monitor plant P status. Transferring this technology to crops would allow precision management of P fertilization, thereby maintaining yields while reducing costs, conserving natural resources, and preventing pollution.