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.

Role of TGF-beta beta and BMP7 in the pathogenesis of oral submucous fibrosis

To understand the molecular pathogenesis of oral submucous fibrosis (OSF), which is a chronic inflammatory disease, gene expression profiling was performed in 10 OSF tissues against 8 pooled normal tissues using oligonucleotide arrays. Microarray results revealed differential expression of 5288 genes (P < a parts per thousand currency sign 0.05 and fold change >= a parts per thousand yen 1.5). Among these, 2884 are upregulated and 2404 are downregulated. Validation employing quantitative real-time PCR and immunohistochemistry confirmed upregulation of transforming growth factor-beta beta 1 (TGF-beta beta 1), TGFBIp, THBS1, SPP1, and TIG1 and downregulation of bone morphogenic protein 7 (BMP7) in OSF tissues. Furthermore, activation of TGF-beta beta pathway was evident in OSF as demonstrated by pSMAD2 strong immunoreactivity. Treatment of keratinocytes and oral fibroblasts by TGF-beta beta confirmed the regulation of few genes identified in microarray including upregulation of connective tissue growth factor, TGM2, THBS1, and downregulation of BMP7, which is a known negative modulator of fibrosis. Taken together, these data suggest activation of TGF-beta beta signaling and suppression of BMP7 expression in the manifestation of OSF.

Uncovering functionally relevant signaling pathways using microarray-based expression profiling

The introduction of microarray technology to the scientific and medical communities has fundamentally altered the way in which we now address basic biomedical questions. Microarrays technology facilitates a more complete and inclusive experimental approach where alterations in the transcript level of entire genomes can be simultaneously assayed in response to a variety of stimuli. Conceptually different approaches to the development of microarray technology have resulted in the generation of two different array formats: oligonucleotide arrays and cDNA arrays. The application of microarray and related technologies to identify specific targets of defined genes that have clearly been implicated in cancer progression requires a specific experimental approach. The objective of tiffs approach is to define changes in transcriptional profile that occur in response to modulating the expression level of the gene to be studied. The resulting altered expression profile can then be viewed as a blueprint by which that gene effects its cellular function. We have used oligonucleotide array-based expression profiling in collaboration with Affymetrix to identify downstream transcriptional targets of the BRCA1 tumor-suppressor gene as a means of defining its function. BRCA1 has been implicated in at least three functional pathways, namely, mediating the cellular response to DNA damage, as a cell cycle checkpoint protein and in the regulation of transcription. The physiological significance of these properties and their implications for the function of BRCA1 as a tumor-suppressor gene remain to be defined.

Induction of GADD45 and JNK/SAPK-dependent apoptosis following inducible expression of BRCA1

The breast cancer susceptibility gene BRCA1 encodes a protein implicated in the cellular response to DNA damage, with postulated roles in homologous recombination as well as transcriptional regulation. To identify downstream target genes, we established cell lines with tightly regulated inducible expression of BRCA1. High-density oligonucleotide arrays were used to analyze gene expression profiles at various times following BRCA1 induction. A major BRCA1 target is the DNA damage-responsive gene GADD45. Induction of BRCA1 triggers apoptosis through activation of c-Jun N-terminal kinase/stress-activated protein kinase (JNK/SAPK), a signaling pathway potentially linked to GADD45 gene family members. The p53-independent induction of GADD45 by BRCA1 and its activation of JNK/SAPK suggest a pathway for BRCA1-induced apoptosis.

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.

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.

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.

Clinical and molecular delineation of the 17q21.31 microdeletion syndrome

Background: The chromosome 17q21.31 microdeletion syndrome is a novel genomic disorder that has originally been identified using high resolution genome analyses in patients with unexplained mental retardation. Aim: We report the molecular and/or clinical characterisation of 22 individuals with the 17q21.31 microdeletion syndrome. Results: We estimate the prevalence of the syndrome to be 1 in 16 000 and show that it is highly underdiagnosed. Extensive clinical examination reveals that developmental delay, hypotonia, facial dysmorphisms including a long face, a tubular or pear-shaped nose and a bulbous nasal tip, and a friendly/amiable behaviour are the most characteristic features. Other clinically important features include epilepsy, heart defects and kidney/urologic anomalies. Using high resolution oligonucleotide arrays we narrow the 17q21.31 critical region to a 424 kb genomic segment (chr17: 41046729-41470954, hg17) encompassing at least six genes, among which is the gene encoding microtubule associated protein tau (MAPT). Mutation screening of MAPT in 122 individuals with a phenotype suggestive of 17q21.31 deletion carriers, but who do not carry the recurrent deletion, failed to identify any disease associated variants. In five deletion carriers we identify a <500 bp rearrangement hotspot at the proximal breakpoint contained within an L2 LINE motif and show that in every case examined the parent originating the deletion carries a common 900 kb 17q21.31 inversion polymorphism, indicating that this inversion is a necessary factor for deletion to occur (p< 10(25)). Conclusion: Our data establish the 17q21.31 microdeletion syndrome as a clinically and molecularly well recognisable genomic disorder.