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.

Using genomic DNA-based probe-selection to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species

High-density oligonucleotide (oligo) arrays are a powerful tool for transcript profiling. Arrays based on GeneChip® technology are amongst the most widely used, although GeneChip® arrays are currently available for only a small number of plant and animal species. Thus, we have developed a method to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species and tested the method by analysing the transcriptome of Brassica oleracea L., a species for which no GeneChip® array is available, using a GeneChip® array designed for Arabidopsis thaliana (L.) Heynh. Genomic DNA from B. oleracea was labelled and hybridised to the ATH1-121501 GeneChip® array. Arabidopsis thaliana probe-pairs that hybridised to the B. oleracea genomic DNA on the basis of the perfect-match (PM) probe signal were then selected for subsequent B. oleracea transcriptome analysis using a .cel file parser script to generate probe mask files. The transcriptional response of B. oleracea to a mineral nutrient (phosphorus; P) stress was quantified using probe mask files generated for a wide range of gDNA hybridisation intensity thresholds. An example probe mask file generated with a gDNA hybridisation intensity threshold of 400 removed > 68 % of the available PM probes from the analysis but retained >96 % of available A. thaliana probe-sets. Ninety-nine of these genes were then identified as significantly regulated under P stress in B. oleracea, including the homologues of P stress responsive genes in A. thaliana. Increasing the gDNA hybridisation intensity thresholds up to 500 for probe-selection increased the sensitivity of the GeneChip® array to detect regulation of gene expression in B. oleracea under P stress by up to 13-fold. Our open-source software to create probe mask files is freely available http://affymetrix.arabidopsis.info/xspecies/ webcite and may be used to facilitate transcriptomic analyses of a wide range of plant and animal species in the absence of custom arrays.

Cellular gene expression altered by human cytomegalovirus: Global monitoring with oligonucleotide arrays

Mechanistic insights to viral replication and pathogenesis generally have come from the analysis of viral gene products, either by studying their biochemical activities and interactions individually or by creating mutant viruses and analyzing their phenotype. Now it is possible to identify and catalog the host cell genes whose mRNA levels change in response to a pathogen. We have used DNA array technology to monitor the level of ≈6,600 human mRNAs in uninfected as compared with human cytomegalovirus-infected cells. The level of 258 mRNAs changed by a factor of 4 or more before the onset of viral DNA replication. Several of these mRNAs encode gene products that might play key roles in virus-induced pathogenesis, identifying them as intriguing targets for further study.

Identification of genes differentially regulated by interferon α, β, or γ using oligonucleotide arrays

The pleiotropic activities of interferons (IFNs) are mediated primarily through the transcriptional regulation of many downstream effector genes. The mRNA profiles from IFN-α, -β, or -γ treatments of the human fibrosarcoma cell line, HT1080, were determined by using oligonucleotide arrays with probe sets corresponding to more than 6,800 human genes. Among these were transcripts for known IFN-stimulated genes (ISGs), the expression of which were consistent with previous studies in which the particular ISG was characterized as responsive to either Type I (α, β) or Type II (γ) IFNs, or both. Importantly, many novel IFN-stimulated genes were identified that were diverse in their known biological functions. For instance, several novel ISGs were identified that are implicated in apoptosis (including RAP46/Bag-1, phospholipid scramblase, and hypoxia inducible factor-1α). Furthermore, several IFN-repressed genes also were identified. These results demonstrate the usefulness of oligonucleotide arrays in monitoring mammalian gene expression on a broad and unprecedented scale. In particular, these findings provide insights into the basic mechanisms of IFN actions and ultimately may contribute to better therapeutic uses for IFNs.

Real-time detection of DNA hybridization and melting on oligonucleotide arrays by using optical wave guides.

The challenge of the Human Genome Project is to increase the rate of DNA sequence acquisition by two orders of magnitude to complete sequencing of the human genome by the year 2000. The present work describes a rapid detection method using a two-dimensional optical wave guide that allows measurement of real-time binding or melting of a light-scattering label on a DNA array. A particulate label on the target DNA acts as a light-scattering source when illuminated by the evanescent wave of the wave guide and only the label bound to the surface generates a signal. Imaging/visual examination of the scattered light permits interrogation of the entire array simultaneously. Hybridization specificity is equivalent to that obtained with a conventional system using autoradiography. Wave guide melting curves are consistent with those obtained in the liquid phase and single-base discrimination is facile. Dilution experiments showed an apparent lower limit of detection at 0.4 nM oligonucleotide. This performance is comparable to the best currently known fluorescence-based systems. In addition, wave guide detection allows manipulation of hybridization stringency during detection and thereby reduces DNA chip complexity. It is anticipated that this methodology will provide a powerful tool for diagnostic applications that require rapid cost-effective detection of variations from known sequences.

An integrated genomic analysis of human glioblastoma Multiforme

Glioblastoma multiforme ( GBM) is the most common and lethal type of brain cancer. To identify the genetic alterations in GBMs, we sequenced 20,661 protein coding genes, determined the presence of amplifications and deletions using high- density oligonucleotide arrays, and performed gene expression analyses using next- generation sequencing technologies in 22 human tumor samples. This comprehensive analysis led to the discovery of a variety of genes that were not known to be altered in GBMs. Most notably, we found recurrent mutations in the active site of isocitrate dehydrogenase 1 ( IDH1) in 12% of GBM patients. Mutations in IDH1 occurred in a large fraction of young patients and in most patients with secondary GBMs and were associated with an increase in overall survival. These studies demonstrate the value of unbiased genomic analyses in the characterization of human brain cancer and identify a potentially useful genetic alteration for the classification and targeted therapy of GBMs.

Gene Expression in Brain: A Window on Ethanol Dependence, Neuroadaptation, and Preference

This article represents the proceedings of a symposium at the 2002 joint RSA/ISBRA Conference in San Francisco, California. The organizer was Paula L. Hoffman and the co-chairs were Paula L. Hoffman and Michael Miles. The presentations were (1) Introduction and overview of the use of DNA microarrays, by Michael Miles; (2) DNA microarray analysis of gene expression in brains of P and NP rats, by Howard J. Edenberg; (3) Gene expression patterns in brain regions of AA and ANA rats, by Wolfgang Sommer; (4) Patterns of gene expression in brains of selected lines of mice that differ in ethanol tolerance, by Boris Tabakoff; (5) Gene expression profiling related to initial sensitivity and tolerance in gamma-protein kinase C mutants, by Jeanne Wehner; and (6) Gene expression patterns in human alcoholic brain: from microarrays to protein profiles, by Joanne Lewohl.

Pregnancy-induced chromatin remodeling in the breast of postmenopausal women.

Early pregnancy and multiparity are known to reduce the risk of women to develop breast cancer at menopause. Based on the knowledge that the differentiation of the breast induced by the hormones of pregnancy plays a major role in this protection, this work was performed with the purpose of identifying what differentiation-associated molecular changes persist in the breast until menopause. Core needle biopsies (CNB) obtained from the breast of 42 nulliparous (NP) and 71 parous (P) postmenopausal women were analyzed in morphology, immunocytochemistry and gene expression. Whereas in the NP breast, nuclei of epithelial cells were large and euchromatic, in the P breast they were small and hyperchromatic, showing strong methylation of histone 3 at lysine 9 and 27. Transcriptomic analysis performed using Affymetrix HG_U133 oligonucleotide arrays revealed that in CNB of the P breast, there were 267 upregulated probesets that comprised genes controlling chromatin organization, transcription regulation, splicing machinery, mRNA processing and noncoding elements including XIST. We concluded that the differentiation process induced by pregnancy is centered in chromatin remodeling and in the mRNA processing reactome, both of which emerge as important regulatory pathways. These are indicative of a safeguard step that maintains the fidelity of the transcription process, becoming the ultimate mechanism mediating the protection of the breast conferred by full-term pregnancy.

DNA methylation profiles of human active and inactive X chromosomes.

X-chromosome inactivation (XCI) is a dosage compensation mechanism that silences the majority of genes on one X chromosome in each female cell. To characterize epigenetic changes that accompany this process, we measured DNA methylation levels in 45,X patients carrying a single active X chromosome (X(a)), and in normal females, who carry one X(a) and one inactive X (X(i)). Methylated DNA was immunoprecipitated and hybridized to high-density oligonucleotide arrays covering the X chromosome, generating epigenetic profiles of active and inactive X chromosomes. We observed that XCI is accompanied by changes in DNA methylation specifically at CpG islands (CGIs). While the majority of CGIs show increased methylation levels on the X(i), XCI actually results in significant reductions in methylation at 7% of CGIs. Both intra- and inter-genic CGIs undergo epigenetic modification, with the biggest increase in methylation occurring at the promoters of genes silenced by XCI. In contrast, genes escaping XCI generally have low levels of promoter methylation, while genes that show inter-individual variation in silencing show intermediate increases in methylation. Thus, promoter methylation and susceptibility to XCI are correlated. We also observed a global correlation between CGI methylation and the evolutionary age of X-chromosome strata, and that genes escaping XCI show increased methylation within gene bodies. We used our epigenetic map to predict 26 novel genes escaping XCI, and searched for parent-of-origin-specific methylation differences, but found no evidence to support imprinting on the human X chromosome. Our study provides a detailed analysis of the epigenetic profile of active and inactive X chromosomes.

Characterization of a genomic signature of pregnancy identified in the breast.

The objective of this study was to comprehensively compare the genomic profiles in the breast of parous and nulliparous postmenopausal women to identify genes that permanently change their expression following pregnancy. The study was designed as a two-phase approach. In the discovery phase, we compared breast genomic profiles of 37 parous with 18 nulliparous postmenopausal women. In the validation phase, confirmation of the genomic patterns observed in the discovery phase was sought in an independent set of 30 parous and 22 nulliparous postmenopausal women. RNA was hybridized to Affymetrix HG_U133 Plus 2.0 oligonucleotide arrays containing probes to 54,675 transcripts, scanned and the images analyzed using Affymetrix GCOS software. Surrogate variable analysis, logistic regression, and significance analysis of microarrays were used to identify statistically significant differences in expression of genes. The false discovery rate (FDR) approach was used to control for multiple comparisons. We found that 208 genes (305 probe sets) were differentially expressed between parous and nulliparous women in both discovery and validation phases of the study at an FDR of 10% and with at least a 1.25-fold change. These genes are involved in regulation of transcription, centrosome organization, RNA splicing, cell-cycle control, adhesion, and differentiation. The results provide initial evidence that full-term pregnancy induces long-term genomic changes in the breast. The genomic signature of pregnancy could be used as an intermediate marker to assess potential chemopreventive interventions with hormones mimicking the effects of pregnancy for prevention of breast cancer.

Immunopathogenesis of Asthma and Atopic Diseases – The specific Role of a selected Panel of Genes in human T helper Cell Differentiation

T helper cell (Th) functions are crucial for proper immune defence against various intra- and extracellular pathogens. According to the specific immune responses, Th cells can be classified into subtypes, Th1 and Th2 cells being the most frequently characterized classes. Th1 and Th2 cells interact with other immune cells by regulating their functions with specific cytokine production. IFN, IL-2 and TNF- are the cytokines predominantly produced by Th1 cells whereas Th2 cells produce Th2-type cytokines, such as IL-4, IL-5 and IL-13. Upon TCR activation and in the presence of polarizing cytokines, Th cells differentiate into effector subtypes from a common precursor cell. IFN and IL-12 are the predominant Th1 polarizing cytokines whereas IL-4 directs Th2 polarization. The cytokines mediate their effects through specific receptor signalling. The differentiation process is complex, involving various signalling molecules and routes, as well as functions of the specific transcription factors. The functions of the Th1/Th2 cells are tightly regulated; however, knowledge on human Th cell differentiation is, as yet, fairly poor. The susceptibility for many immune-mediated disorders often originates from disturbed Th cell responses. Thus, research is needed for defining the molecular mechanisms involved in the differentiation and balanced functions of the Th cells. Importantly, the new information obtained will be crucial for a better understanding of the pathogenesis of immune-mediated disorders, such as asthma or autoimmune diseases. In the first subproject of this thesis, the role of genetic polymorphisms in the human STAT6, GATA3 and STAT4 genes were investigated for asthma or atopy susceptibility in Finnish asthma families by association analysis. These genes code for key transcription factors regulating Th cell differentiation. The study resulted in the identification of a GATA3 haplotype that associated with asthma and related traits (high serum IgE level). In the second subproject, an optimized method for human primary T cell transfection and enrichment was established. The method can be utilized for functional studies for the selected genes of interest. The method was also utilized in the third subproject, which aimed at the identification of novel genes involved in early human Th cell polarization (0-48h) using genome-wide oligonucleotide arrays. As a result, numerous genes and ESTs with known or unknown functions were identified in the study. Using an shRNA knockdown approach, a panel of novel IL-4/STAT6 regulated genes were identified in the functional studies of the genes. Moreover, one of the genes, NDFIP2, with a previously uncharacterized role in the human Th differentiation, was observed to promote IFN production of the differentiated Th1 cells. Taken together, the results obtained have revealed potential new relevant candidate genes serving as a basis for further studies characterizing the detailed networks involved in the human Th cell differentiation as well as in the genetic susceptibility of Th-mediated immune disorders.

Transcriptome analysis of grain development in hexaploid wheat

Background: Hexaploid wheat is one of the most important cereal crops for human nutrition. Molecular understanding of the biology of the developing grain will assist the improvement of yield and quality traits for different environments. High quality transcriptomics is a powerful method to increase this understanding. Results: The transcriptome of developing caryopses from hexaploid wheat ( Triticum aestivum, cv. Hereward) was determined using Affymetrix wheat GeneChip (R) oligonucleotide arrays which have probes for 55,052 transcripts. Of these, 14,550 showed significant differential regulation in the period between 6 and 42 days after anthesis ( daa). Large changes in transcript abundance were observed which were categorised into distinct phases of differentiation ( 6 - 10 daa), grain fill ( 12 - 21 daa) and desiccation/maturation ( 28 - 42 daa) and were associated with specific tissues and processes. A similar experiment on developing caryopses grown with dry and/or hot environmental treatments was also analysed, using the profiles established in the first experiment to show that most environmental treatment effects on transcription were due to acceleration of development, but that a few transcripts were specifically affected. Transcript abundance profiles in both experiments for nine selected known and putative wheat transcription factors were independently confirmed by real time RT-PCR. These expression profiles confirm or extend our knowledge of the roles of the known transcription factors and suggest roles for the unknown ones. Conclusion: This transcriptome data will provide a valuable resource for molecular studies on wheat grain. It has been demonstrated how it can be used to distinguish general developmental shifts from specific effects of treatments on gene expression and to diagnose the probable tissue specificity and role of transcription factors.

Effect of Ginkgo biloba (EGb 761®) on differential gene expression

Supplementation of diets with plant extracts such as ginkgo biloba extract (EGb 761®) (definition see editorial) for health and prevention of degenerative diseases is popular. However, it is often difficult to analyse the biological activities of plant extracts due to their complex nature and the possible synergistic and/or antagonistic effects of their components. Genome-wide expression monitoring with high-density oligonucleotide arrays provides one way to examine the molecular targets of plant extracts and may prove a useful tool in evaluating their therapeutic claims. Here, we will briefly describe some of our work on the effect of EGb 761® on differential gene expression in relation to its potential anti-carcinogenic, photoprotective and neuromodulatory properties.

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.