DYSREGULATION OF MEOX2 FOLLOWING WT1 MUTATION IN KIDNEY DEVELOPMENT AND WILMS TUMORIGENESIS

Wilms tumor (WT) is a childhood tumor of the kidney and a productive model for understanding the role of genetic alteration and interactions in tumorigenesis. The Wilms tumor gene 1 (WT1) is a transcriptional factor and one of the few genes known to have genetic alterations in WT and has been shown be inactivated in 20% of WTs. However, the mechanisms of how WT1 mutations lead to Wilms tumorigenesis and its influence on downstream genes are unknown. Since it has been established that WT1 is a transcriptional regulator, it has been hypothesized that the loss of WT1 leads to the dysregulation of downstream genes, in turn result in the formation of WTs. To identify the dysregulated downstream genes following WT1 mutations, an Affymetrix GeneChip Human Genome Array was previously conducted to assess the differentially expressed genes in the WT1-wildtype human and WT1-mutant human WTs. Approximately 700 genes were identified as being significantly dysregulated. These genes were further prioritized based on their statistical significance, fold change, chromosomal region, spatial pattern of gene expression and known or putative cellular functions. Mesenchyme homeobox 2 (MEOX2) was one of the most significantly upregulated genes in WT1-mutant WT. MEOX2 is known to play a role in cell proliferation, apoptosis, and differentiation. In addition to its biological roles, it is expressed during early kidney development in the condensed mesenchyme similar to WT1. Furthermore, the use of the Match® web-based tool from the BIOBASE Biological Data base identified a significant predicted WT1 binding site within the first intron of MEOX2. The similarity in spatial gene expression in the developing kidney and the significant predicted WT1 binding site found in the first intron of MEOX2 lead to the development of my hypothesis that MEOX2 is upregulated via a WT1-dependent manner. Here as a part of my master’s work, I have validated the Affymetrix GeneChip Human Genome Array data using an independent set of Wilms tumors. MEOX2 remained upregulated in the mutant WT1 Wilms tumor by 41-fold. Wt1 and Meox2 gene expression were assessed in murine newborn kidney; both Wt1 and Meox2 were expressed in the condensed, undifferentiated metanephric mesenchyme. I have shown that the in vivo ablation of Wt1 during embryonic development at embryonic day (E) 13.5 resulted in the slight increase of Meox2 gene expression by two fold. In order to functionally demonstrate the effect of the loss of Wt1 on Meox2 gene expression in undifferentiated metanephric mesenchyme, I have generated a kidney mesenchymal cell line to genetically ablate Wt1 in vitro by adenoviral infection. The ablation of Wt1 in the kidney mesenchymal cell line resulted in the upregulation of Meox2 by 61-fold. Moreover, the upregulation of Meox2 resulted in the significant induction of p21 and Itgb5. In addition to the dysregulation of these genes the ablation of Wt1 in the kidney mesenchymal cells resulted in decrease in cell growth and loss of cellular adherence. However, it is uncertain whether the upregulation of Meox2 caused this particular cellular phenotype. Overall, I have demonstrated that the upregulation of Meox2 is Wt1-dependent during early kidney development.

Transcriptional analysis of left-sided colitis, pancolitis, and ulcerative colitis-associated dysplasia.

BACKGROUND It is unknown why patients with extensive ulcerative colitis (UC) have a higher risk of colorectal cancer compared with patients with left-sided UC. This study characterizes the inflammatory processes in left-sided UC, pancolitis, and UC-associated dysplasia at the transcriptional level to identify potential biomarkers and transcripts of importance for the carcinogenic behavior of chronic inflammation. METHODS The Affymetrix GeneChip Human Genome U133 Plus 2.0 was applied on colonic biopsies from UC patients with left-sided UC, pancolitis, dysplasia, and controls. Reverse transcription polymerase chain reaction and immunohistochemistry were performed for validating selected transcripts in the initial cohort and in 2 independent cohorts of patients with UC. Microarray data were analyzed by principal component analysis, and reverse transcription polymerase chain reaction and immunohistochemistry data by the Wilcoxon's rank-sum test. RESULTS The principal component analysis results revealed separate clusters for left-sided UC, pancolitis, dysplasia, and controls. Close clustering of dysplastic and pancolitic samples indicated similarities in gene expression. Indeed, 101 and 656 parallel upregulated and downregulated transcripts, respectively, were identified in specimens from dysplasia and pancolitis. Validation of selected transcripts hereof identified insulin receptor alpha (INSRA) and MAP kinase interacting serine/threonine kinase 2 (MKNK2) with an enhanced expression in dysplasia compared with left-sided UC and controls, whereas laminin γ2 (LAMC2) was found with a lower expression in dysplasia compared with the remaining 3 groups. CONCLUSIONS This study demonstrates pancolitis and left-sided UC as distinct inflammatory processes at the transcriptional level, and identifies INSRA, MKNK2, and LAMC2 as potential critical transcripts in the inflammation-driven preneoplastic process of UC.

Identification of diabetic retinopathy genes through a genome-wide association study among Mexican-Americans from Starr County, Texas

To identify genetic susceptibility loci for severe diabetic retinopathy, 286 Mexican-Americans with type 2 diabetes from Starr County, Texas completed detailed physical and ophthalmologic examinations including fundus photography for diabetic retinopathy grading. 103 individuals with moderate-to-severe non-proliferative diabetic retinopathy or proliferative diabetic retinopathy were defined as cases for this study. DNA samples extracted from study subjects were genotyped using the Affymetrix GeneChip® Human Mapping 100K Set, which includes 116,204 single nucleotide polymorphisms (SNPs) across the whole genome. Single-marker allelic tests and 2- to 8-SNP sliding-window Haplotype Trend Regression implemented in HelixTreeTM were first performed with these direct genotypes to identify genes/regions contributing to the risk of severe diabetic retinopathy. An additional 1,885,781 HapMap Phase II SNPs were imputed from the direct genotypes to expand the genomic coverage for a more detailed exploration of genetic susceptibility to diabetic retinopathy. The average estimated allelic dosage and imputed genotypes with the highest posterior probabilities were subsequently analyzed for associations using logistic regression and Fisher's Exact allelic tests, respectively. To move beyond these SNP-based approaches, 104,572 directly genotyped and 333,375 well-imputed SNPs were used to construct genetic distance matrices based on 262 retinopathy candidate genes and their 112 related biological pathways. Multivariate distance matrix regression was then used to test hypotheses with genes and pathways as the units of inference in the context of susceptibility to diabetic retinopathy. This study provides a framework for genome-wide association analyses, and implicated several genes involved in the regulation of oxidative stress, inflammatory processes, histidine metabolism, and pancreatic cancer pathways associated with severe diabetic retinopathy. Many of these loci have not previously been implicated in either diabetic retinopathy or diabetes. In summary, CDC73, IL12RB2, and SULF1 had the best evidence as candidates to influence diabetic retinopathy, possibly through novel biological mechanisms related to VEGF-mediated signaling pathway or inflammatory processes. While this study uncovered some genes for diabetic retinopathy, a comprehensive picture of the genetic architecture of diabetic retinopathy has not yet been achieved. Once fully understood, the genetics and biology of diabetic retinopathy will contribute to better strategies for diagnosis, treatment and prevention of this disease.^

Methionine adenosyltransferase 1A knockout mice are predisposed to liver injury and exhibit increased expression of genes involved in proliferation

Liver-specific and nonliver-specific methionine adenosyltransferases (MATs) are products of two genes, MAT1A and MAT2A, respectively, that catalyze the formation of S-adenosylmethionine (AdoMet), the principal biological methyl donor. Mature liver expresses MAT1A, whereas MAT2A is expressed in extrahepatic tissues and is induced during liver growth and dedifferentiation. To examine the influence of MAT1A on hepatic growth, we studied the effects of a targeted disruption of the murine MAT1A gene. MAT1A mRNA and protein levels were absent in homozygous knockout mice. At 3 months, plasma methionine level increased 776% in knockouts. Hepatic AdoMet and glutathione levels were reduced by 74 and 40%, respectively, whereas S-adenosylhomocysteine, methylthioadenosine, and global DNA methylation were unchanged. The body weight of 3-month-old knockout mice was unchanged from wild-type littermates, but the liver weight was increased 40%. The Affymetrix genechip system and Northern and Western blot analyses were used to analyze differential expression of genes. The expression of many acute phase-response and inflammatory markers, including orosomucoid, amyloid, metallothionein, Fas antigen, and growth-related genes, including early growth response 1 and proliferating cell nuclear antigen, is increased in the knockout animal. At 3 months, knockout mice are more susceptible to choline-deficient diet-induced fatty liver. At 8 months, knockout mice developed spontaneous macrovesicular steatosis and predominantly periportal mononuclear cell infiltration. Thus, absence of MAT1A resulted in a liver that is more susceptible to injury, expresses markers of an acute phase response, and displays increased proliferation.

Expression of HOXB2, a retinoic acid signaling target in pancreatic cancer and pancreatic intraepithelial neoplasia

Purpose: Despite significant progress in understanding the molecular pathology of pancreatic cancer and its precursor lesion: pancreatic intraepithelial neoplasia (PanIN), there remain no molecules with proven clinical utility as prognostic or therapeutic markers. Here, we used oligonucleotide microarrays to interrogate mRNA expression of pancreatic cancer tissue and normal pancreas to identify novel molecular pathways dysregulated in the development and progression of pancreatic cancer. Experimental Design: RNA was hybridized to Affymetrix Genechip HG-U133 oligonucleotide microarrays. A relational database integrating data from publicly available resources was created to identify candidate genes potentially relevant to pancreatic cancer. The protein expression of one candidate, homeobox B2 (HOXB2), in PanIN and pancreatic cancer was assessed using immunohistochemistry. Results: We identified aberrant expression of several components of the retinoic acid (RA) signaling pathway (RARa, MUC4, Id-1, MMP9, uPAR, HB-EGF, HOXB6, and HOXB2), many of which are known to be aberrantly expressed in pancreatic cancer and Pan IN. HOXB2, a downstream target of RA, was up-regulated 6.7-fold in pancreatic cancer compared with normal pancreas. Immunohistochemistry revealed ectopic expression of HOXB2 in 15% of early Pan IN lesions and 48 of 128 (38%) pancreatic cancer specimens. Expression of HOXB2 was associated with nonresectable tumors and was an independent predictor of poor survival in resected tumors. Conclusions: We identified aberrant expression of RA signaling components in pancreatic cancer, including HOXB2, which was expressed in a proportion of PanIN lesions. Ectopic expression of HOXB2 was associated with a poor prognosis for all patients with pancreatic cancer and was an independent predictor of survival in patients who underwent resection.

Microarray analysis of expression of cell death-associated genes in rat spinal cord cells exposed to cyclic tensile stresses in vitro

The application of mechanical insults to the spinal cord results in profound cellular and molecular changes, including the induction of neuronal cell death and altered gene expression profiles. Previous studies have described alterations in gene expression following spinal cord injury, but the specificity of this response to mechanical stimuli is difficult to investigate in vivo. Therefore, we have investigated the effect of cyclic tensile stresses on cultured spinal cord cells from E15 Sprague-Dawley rats, using the FX3000 Flexercell Strain Unit. We examined cell morphology and viability over a 72 hour time course. Microarray analysis of gene expression was performed using the Affymetrix GeneChip System, where categorization of identified genes was performed using the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) systems. Changes in expression of 12 genes were validated with quantitative real-time reverse transcription polymerase chain reaction (RT-PCR).

Identification and functional analysis of novel genes expressed in the Anterior Visceral Endoderm

During early vertebrate development, the correct establishment of the body axes is critical. The anterior pole of the mouse embryo is established when Distal Visceral Endoderm (DVE) cells migrate to form the Anterior Visceral Endoderm (AVE). Symmetrical expression of Lefty1, Cer1 and Dkk1 determines the direction of DVE migration and the future anterior side. In addition to the establishment of the Anterior-Posterior axis, the AVE has also been implicated in anterior neural specification. To better understand the role of the AVE in these processes, we have performed a differential screening using Affymetrix GeneChip technology with AVE cells isolated from cer1P-EGFP transgenic mouse embryos. We found 175 genes which were upregulated in the AVE and 36 genes in the Proximal-posterior sample. Using DAVID software, we characterized the AVE cell population regarding cellular component, molecular function and biological processes. Among the genes that were found to be upregulated in the AVE, several novel genes were identified. Four of these transcripts displaying high-fold change in the AVE were further characterized by in situ hybridization in early stages of development in order to validate the screening. From those four selected genes, one, denominated Adtk1, was chosen to be functionally characterized by targeted inactivation in ES cells. Adtk1 encodes for a serine/threonine kinase. Adtk1 null mutants are smaller and present short limbs due to decreased mineralization, suggesting a potential role in chondrogenesis during limb development. Taken together, these data point to the importance of reporting novel genes present in the AVE.

Identification of differentially expressed genes in the heart precursor cells of the chick embryo

Genetic evidence has implicated several genes as being critical for heart development. However, the inducers of these genes as well as their targets and pathways they are involved with, remain largely unknown. Previous studies in the avian embryo showed that at HH4 Cerberus (cCer) transcripts are detected in the anterior endomesoderm including the heart precursor cells and later in the left lateral plate mesoderm. We have identified a promoter element of chick cCer able to drive EGFP expression in a population of cells that consistently exit from the anterior primitive streak region, from as early as stage HH3+, and that later will populate the heart. Using this promoter element as a tool allowed us to identify novel genes previously not known to potentially play a role in heart development. In order to identify and study genes expressed and involved in the correct development and differentiation of the vertebrate heart precursor cell (HPC) lineages, a differential screening using Affymetrix GeneChip system technologies was performed. Remarkably, this screening led to the identification of more than 700 transcripts differentially expressed in the heart forming regions (HFR). Bioinformatic tools allowed us to filter the large amount of data generated from this approach and to select a few transcripts for in vivo validation. Whole-mount in situ hybridization and sectioning of selected genes showed heart and vascular expression patterns for these transcripts during early chick development. We have developed an effective strategy to specifically identify genes that are differentially expressed in the HPC lineages. Within this set we have identified several genes that are expressed in the heart, blood and vascular lineages, which are likely to play a role in their development. These genes are potential candidates for future functional studies on early embryonic patterning.

Genic expression of the uterine leiomyoma in reproductive-aged women after treatment with goserelin

Objective: To identify the genes presenting different expression in uterine leiomyomas after goserelin treatment. Design: Retrospective analyses of tissue obtained in a prospective clinical study. Setting: School of Medicine of the University of Sao Paulo. Patient(s): 30 nulliparous black women aged 20 to 45 years with symptoms of uterine leiomyoma, uterine volume over 300 mL, and surgical indications for myomectomy. Intervention(s): Fifteen patients were given a monthly dose of 3.6 mg of goserelin over 3 months before surgery (group A), and 15 patients underwent surgery without any previous treatment (group B). Five random samples from each group were analyzed using the microarray technique with the Affymetrix platform (GeneChip Rat Genome 230 2.0 Array). Main Outcome Measure(s): Quantification of transcript expression levels of uterine fibroids in patients treated or not treated with goserelin. Result(s): Of the total of 47,000 sequences that were analyzed, representing approximately 38,500 human genes already characterized, we found a differential expression of 174 genes. Of these, 70 were up-regulated (33 genes with known function) and 104 were down-regulated (65 genes with known function) in samples from group A (treated) when compared with group B (nontreated). Conclusion(s): The genic expression of uterine leiomyomas changes in women who have had goserelin treatment when compared with nontreated patients. (Fertil Steril (R) 2010; 94: 1072-7. (C) 2010 by American Society for Reproductive Medicine.)

IQRray, a new method for Affymetrix microarray quality control, and the homologous organ conservation score, a new benchmark method for quality control metrics.

MOTIVATION: Microarray results accumulated in public repositories are widely reused in meta-analytical studies and secondary databases. The quality of the data obtained with this technology varies from experiment to experiment, and an efficient method for quality assessment is necessary to ensure their reliability. RESULTS: The lack of a good benchmark has hampered evaluation of existing methods for quality control. In this study, we propose a new independent quality metric that is based on evolutionary conservation of expression profiles. We show, using 11 large organ-specific datasets, that IQRray, a new quality metrics developed by us, exhibits the highest correlation with this reference metric, among 14 metrics tested. IQRray outperforms other methods in identification of poor quality arrays in datasets composed of arrays from many independent experiments. In contrast, the performance of methods designed for detecting outliers in a single experiment like Normalized Unscaled Standard Error and Relative Log Expression was low because of the inability of these methods to detect datasets containing only low-quality arrays and because the scores cannot be directly compared between experiments. AVAILABILITY AND IMPLEMENTATION: The R implementation of IQRray is available at: ftp://lausanne.isb-sib.ch/pub/databases/Bgee/general/IQRray.R. CONTACT: Marta.Rosikiewicz@unil.ch SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Bone anabolic potential of soy isoflavones and plant extracts and genomic changes during differentiation of hPOB-tert osteoblast-like cells

Summary For the nutritional management of bone health and the prevention of osteoporosis it is important to identify nutrients that positively influence the bone remodeling process at the cellular level. Soy isoflavones show promising osteoprotective effects in animals and humans but their mechanism of action in bone cells is yet poorly understood. Firstly, soy tissue cultures were characterized as a new and optimized source of isoflavones. A large variability in the isoflavone content was observed and high-producing strains (46.3 mg/g dry wt isoflavones) were identified. In the Ishikawa cells bioassay, the estrogenicity of isoflavones was confirmed to be 1000 to 10000 less than 17Mestradiol and that of the malonyl forms was shown for the first time (EC50 of 350 nM and 1880 nM for malonylgenistin and malonyldaidzin, respectively). The estrogenic activity of soya tissue culture extracts correlated to their isoflavone content. Secondly, the effects of phytonutrients on BMP-2 gene expression and on the mevalonate synthesis pathway, as key mediators of bone formation, were investigated. Dietary achievable concentrations of genistein and daidzein (10vM), and statins (4xM) but not 17M estradiol (10nM), induced BMP-2 gene expression (by up to 3-fold) and inhibited the cholesterol biosynthetic pathway (by up to 50%) in the human osteoblastic cell line hP0B¬tert. In addition, several plant extracts (Cyperus rotundus, Lindera benzoin and Cnidium monnieri) induced BMP-2 gene expression but this induction was not restricted to the inhibition of the cholesterol synthesis pathway neither to the estrogenicity. Finally, the gene expression profiles during hP0B-tert differentiation induced by vitamin D and dexamethasone were analyzed with the Affymetrix human GeneChip. 1665 different genes and 98 ESTs were significantly regulated. The expression profiles of bone-related genes was largely in agreement with previously documented patterns, supporting the physiological relevance of the genomic results and the hP0B-tert cell line as a valid model of human osteoblast differentiation. The expression of alternative differentiation markers during the osteogenic treatment of hP0B-tert cells indicated that the adipocyte and myoblast differentiation pathways were repressed, confirming that these culture conditions allowed only osteoblast differentiation. The gene ontology analysis identified further sub-groups of genes that may be involved in the bone formation process. Aims of the thesis In order to define new strategies for the nutritional management of bone health and for the prevention of osteoporosis the major goal of the present work was to investigate the potential of phytonutrients to positively modulate the bone formation process at the cellular level and, in particular: 1.To select and optimise alternative plant sources containing high levels of isoflavones with estrogenic activity (Chapter 3). 2.To compare the effects of statins and phytonutrients on BMP-2 gene expression and on the mevalonate synthesis pathway and to select new plant extracts with a bone anabolic potential (Chapter 4). 3.To further characterize the new human periosteal cell line, hP0B-tert, as a bone- formation model, by elucidating its gene expression profile during differentiation induced by vitamin D and dexamethasone (Chapter 5).

Affymetrix Whole-Transcript Human Gene 1.0 ST array is highly concordant with standard 3' expression arrays.

The recently released Affymetrix Human Gene 1.0 ST array has two major differences compared with standard 3' based arrays: (i) it interrogates the entire mRNA transcript, and (ii) it uses DNA targets. To assess the impact of these differences on array performance, we performed a series of comparative hybridizations between the Human Gene 1.0 ST and the Affymetrix HG-U133 Plus 2.0 and the Illumina HumanRef-8 BeadChip arrays. Additionally, both RNA and DNA targets were hybridized on HG-U133 Plus 2.0 arrays. The results show that the overall reproducibility of the Gene 1.0 ST array is best. When looking only at the high intensity probes, the reproducibility of the Gene 1.0 ST array and the Illumina BeadChip array is equally good. Concordance of array results was assessed using different inter-platform mappings. Agreements are best between the two labeling protocols using HG-U133 Plus 2.0 array. The Gene 1.0 ST array is most concordant with the HG-U133 array hybridized with cDNA targets. This may reflect the impact of the target type. Overall, the high degree of correspondence provides strong evidence for the reliability of the Gene 1.0 ST array.

Regional variations in ABC transporter expression along the mouse intestinal tract.

The ATP-binding cassette (ABC) family of proteins comprise a group of membrane transporters involved in the transport of a wide variety of compounds, such as xenobiotics, vitamins, lipids, amino acids, and carbohydrates. Determining their regional expression patterns along the intestinal tract will further characterize their transport functions in the gut. The mRNA expression levels of murine ABC transporters in the duodenum, jejunum, ileum, and colon were examined using the Affymetrix MuU74v2 GeneChip set. Eight ABC transporters (Abcb2, Abcb3, Abcb9, Abcc3, Abcc6, Abcd1, Abcg5, and Abcg8) displayed significant differential gene expression along the intestinal tract, as determined by two statistical models (a global error assessment model and a classic ANOVA, both with a P < 0.01). Concordance with semiquantitative real-time PCR was high. Analyzing the promoters of the differentially expressed ABC transporters did not identify common transcriptional motifs between family members or with other genes; however, the expression profile for Abcb9 was highly correlated with fibulin-1, and both genes share a common complex promoter model involving the NFkappaB, zinc binding protein factor (ZBPF), GC-box factors SP1/GC (SP1F), and early growth response factor (EGRF) transcription binding motifs. The cellular location of another of the differentially expressed ABC transporters, Abcc3, was examined by immunohistochemistry. Staining revealed that the protein is consistently expressed in the basolateral compartment of enterocytes along the anterior-posterior axis of the intestine. Furthermore, the intensity of the staining pattern is concordant with the expression profile. This agrees with previous findings in which the mRNA, protein, and transport function of Abcc3 were increased in the rat distal intestine. These data reveal regional differences in gene expression profiles along the intestinal tract and demonstrate that a complete understanding of intestinal ABC transporter function can only be achieved by examining the physiologically distinct regions of the gut.

Segregation of regulatory polymorphisms with effects on the gluteus medius transcriptome in a purebred pig population

Background: The main goal of the present study was to analyse the genetic architecture of mRNA expression in muscle, a tissue with an outmost economic importance for pig breeders. Previous studies have used F2 crosses to detect porcine expression QTL (eQTL), so they contributed with data that mostly represents the between-breed component of eQTL variation. Herewith, we have analysed eQTL segregation in an outbred Duroc population using two groups of animals with divergent fatness profiles. This approach is particularly suitable to analyse the within-breed component of eQTL variation, with a special emphasis on loci involved in lipid metabolism. Methodology/Principal Findings: GeneChip Porcine Genome arrays (Affymetrix) were used to determine the mRNA expression levels of gluteus medius samples from 105 Duroc barrows. A whole-genome eQTL scan was carried out with a panel of 116 microsatellites. Results allowed us to detect 613 genome-wide significant eQTL unevenly distributed across the pig genome. A clear predominance of trans- over cis-eQTL, was observed. Moreover, 11 trans-regulatory hotspots affecting the expression levels of four to 16 genes were identified. A Gene Ontology study showed that regulatory polymorphisms affected the expression of muscle development and lipid metabolism genes. A number of positional concordances between eQTL and lipid trait QTL were also found, whereas limited evidence of a linear relationship between muscle fat deposition and mRNA levels of eQTL regulated genes was obtained. Conclusions/Significance: Our data provide substantial evidence that there is a remarkable amount of within-breed genetic variation affecting muscle mRNA expression. Most of this variation acts in trans and influences biological processes related with muscle development, lipid deposition and energy balance. The identification of the underlying causal mutations and the ascertainment of their effects on phenotypes would allow gaining a fundamental perspective about how complex traits are built at the molecular level.

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.