affy--analysis of Affymetrix GeneChip data at the probe level

Genome-wide effects of acute progressive feed restriction in liver and white adipose tissue

Acute progressive feed restriction (APFR) represents a specific form of caloric restriction in which feed availability is increasingly curtailed over a period of a few days to a few weeks. It is often used for control animals in toxicological and pharmacological studies on compounds causing body weight loss to equalize weight changes between experimental and control groups and thereby, intuitively, to also set their metabolic states to the same phase. However, scientific justification for this procedure is lacking. In the present study, we analyzed by microarrays the impact on hepatic gene expression in rats of two APFR regimens that caused identical diminution of body weight (19%) but differed slightly in duration (4 vs. 10 days). In addition, white adipose tissue (WAT) was also subjected to the transcriptomic analysis on day-4. The data revealed that the two regimens led to distinct patterns of differentially expressed genes in liver, albeit some major pathways of energy metabolism were similarly affected (particularly fatty acid and amino acid catabolism). The reason for the divergence appeared to be entrainment by the longer APFR protocol of peripheral oscillator genes, which resulted in derailment of circadian rhythms and consequent interaction of altered diurnal fluctuations with metabolic adjustments in gene expression activities. WAT proved to be highly unresponsive to the 4-day APFR as only 17 mRNA levels were influenced by the treatment. This study demonstrates that body weight is a poor proxy of metabolic state and that the customary protocols of feed restriction can lead to rhythm entrainment.

Correction of technical bias in clinical microarray data improves concordance with known biological information

A method is presented to correct for technical bias in clinical microarray data, which increases concordance with known biological relationships.

The performance of gene expression microarrays has been well characterized using controlled reference samples, but the performance on clinical samples remains less clear. We identified sources of technical bias affecting many genes in concert, thus causing spurious correlations in clinical data sets and false associations between genes and clinical variables. We developed a method to correct for technical bias in clinical microarray data, which increased concordance with known biological relationships in multiple data sets.

Dynamics of the yeast transcriptome during wine fermentation reveals a novel fermentation stress response

In this study, genome-wide expression analyses were used to study the response of Saccharomyces cerevisiae to stress throughout a 15-day wine fermentation. Forty per cent of the yeast genome significantly changed expression levels to mediate long-term adaptation to fermenting grape must. Among the genes that changed expression levels, a group of 223 genes was identified, which was designated as fermentation stress response (FSR) genes that were dramatically induced at various points during fermentation. FSR genes sustain high levels of induction up to the final time point and exhibited changes in expression levels ranging from four- to 80-fold. The FSR is novel; 62% of the genes involved have not been implicated in global stress responses and 28% of the FSR genes have no functional annotation. Genes involved in respiratory metabolism and gluconeogenesis were expressed during fermentation despite the presence of high concentrations of glucose. Ethanol, rather than nutrient depletion, seems to be responsible for entry of yeast cells into the stationary phase.

Human and chimpanzee gene expression differences replicated in mice fed different diets

Although the human diet is markedly different from the diets of closely related primate species, the influence of diet on phenotypic and genetic differences between humans and other primates is unknown. In this study, we analyzed gene expression in laboratory mice fed diets typical of humans and of chimpanzees. The effects of human diets were found to be significantly different from that of a chimpanzee diet in the mouse liver, but not in the brain. Importantly, 10% of the genes that differ in their expression between humans and chimpanzee livers differed also between the livers of mice fed the human and chimpanzee diets. Furthermore, both the promoter sequences and the amino acid sequences of these diet-related genes carry more differences between humans and chimpanzees than random genes. Our results suggest that the mouse can be used to study at least some aspects of human-specific traits.

Variations in the progranulin gene affect global gene expression in frontotemporal lobar degeneration

Frontotemporal lobar degeneration is a fatal neurodegenerative disease that results in progressive decline in behavior, executive function and sometimes language. Disease mechanisms remain poorly understood. Recently, however, the DNA- and RNA-binding protein TDP-43 has been identified as the major protein present in the hallmark inclusion bodies of frontotemporal lobar degeneration with ubiquitinated inclusions (FTLD-U), suggesting a role for transcriptional dysregulation in FTLD-U pathophysiology. Using the Affymetrix U133A microarray platform, we profiled global gene expression in both histopathologically affected and unaffected areas of human FTLD-U brains. We then characterized differential gene expression with biological pathway analyses, cluster and principal component analyses, and subgroup analyses based on brain region and progranulin (GRN) gene status. Comparing 17 FTLD-U brains to 11 controls, we identified 414 upregulated and 210 downregulated genes in frontal cortex (P-value < 0.001). Moreover, cluster and principal component analyses revealed that samples with mutations or possibly pathogenic variations in the GRN gene (GRN+, 7/17) had an expression signature that was distinct from both normal controls and FTLD-U samples lacking GRN gene variations (GRN-, 10/17). Within the subgroup of GRN+ FTLD-U, we found >1300 dysregulated genes in frontal cortex (P-value < 0.001), many participating in pathways uniquely dysregulated in the GRN+ cases. Our findings demonstrate a distinct molecular phenotype for GRN+ FTLD-U, not readily apparent on clinical or histopathological examination, suggesting distinct pathophysiological mechanisms for GRN+ and GRN- subtypes of FTLD-U. In addition, these data from a large number of human brains provide a valuable resource for future testing of disease hypotheses.

Genome-wide microarray expression analysis of CD4+ T Cells from nonobese diabetic congenic mice identifies Cd55 (Daf1) and Acadl as candidate genes for type 1 diabetes

NOD.Idd3/5 congenic mice have insulin-dependent diabetes (Idd) regions on chromosomes 1 (Idd5) and 3 (Idd3) derived from the nondiabetic strains B10 and B6, respectively. NOD.Idd3/5 mice are almost completely protected from type 1 diabetes (T1D) but the genes within Idd3 and Idd5 responsible for the disease-altering phenotype have been only partially characterized. To test the hypothesis that candidate Idd genes can be identified by differential gene expression between activated CD4+ T cells from the diabetes-susceptible NOD strain and the diabetes-resistant NOD.Idd3/5 congenic strain, genome-wide microarray expression analysis was performed using an empirical Bayes method. Remarkably, 16 of the 20 most differentially expressed genes were located in the introgressed regions on chromosomes 1 and 3, validating our initial hypothesis. The two genes with the greatest differential RNA expression on chromosome 1 were those encoding decay-accelerating factor (DAF, also known as CD55) and acyl-coenzyme A dehydrogenase, long chain, which are located in the Idd5.4 and Idd5.3 regions, respectively. Neither gene has been implicated previously in the pathogenesis of T1D. In the case of DAF, differential expression of mRNA was extended to the protein level; NOD CD4+ T cells expressed higher levels of cell surface DAF compared with NOD.Idd3/5 CD4+ T cells following activation with anti-CD3 and -CD28. DAF up-regulation was IL-4 dependent and blocked under Th1 conditions. These results validate the approach of using congenic mice together with genome-wide analysis of tissue-specific gene expression to identify novel candidate genes in T1D.

The Rcs phosphorelay is a cell envelope stress response activated by peptidoglycan stress and contributes to intrinsic antibiotic resistance

Gram-negative bacteria possess stress responses to maintain the integrity of the cell envelope. Stress sensors monitor outer membrane permeability, envelope protein folding, and energization of the inner membrane. The systems used by gram-negative bacteria to sense and combat stress resulting from disruption of the peptidoglycan layer are not well characterized. The peptidoglycan layer is a single molecule that completely surrounds the cell and ensures its structural integrity. During cell growth, new peptidoglycan subunits are incorporated into the peptidoglycan layer by a series of enzymes called the penicillin-binding proteins (PBPs). To explore how gram-negative bacteria respond to peptidoglycan stress, global gene expression analysis was used to identify Escherichia coli stress responses activated following inhibition of specific PBPs by the beta-lactam antibiotics amdinocillin (mecillinam) and cefsulodin. Inhibition of PBPs with different roles in peptidoglycan synthesis has different consequences for cell morphology and viability, suggesting that not all perturbations to the peptidoglycan layer generate equivalent stresses. We demonstrate that inhibition of different PBPs resulted in both shared and unique stress responses. The regulation of capsular synthesis (Rcs) phosphorelay was activated by inhibition of all PBPs tested. Furthermore, we show that activation of the Rcs phosphorelay increased survival in the presence of these antibiotics, independently of capsule synthesis. Both activation of the phosphorelay and survival required signal transduction via the outer membrane lipoprotein RcsF and the response regulator RcsB. We propose that the Rcs pathway responds to peptidoglycan damage and contributes to the intrinsic resistance of E. coli to beta-lactam antibiotics.

Galectin-10 mRNA is overexpressed in peripheral blood of aspirin-induced asthma

BACKGROUND

In sensitive patients, aspirin is associated with nasal and bronchial inflammation, eliciting local symptoms. Although the disease is clinically well characterized, its physiopathology is incompletely understood and noninvasive procedures, allowing an effective distinction between aspirin-induced asthma (AIA) and aspirin-tolerant asthma (ATA) are missing.

OBJECTIVES

The aims of the study were to compare AIA and ATA cohorts for clinical characteristics and to screen peripheral blood for differential mRNA expression.

METHODS

Patients experiencing symptoms following aspirin ingestion were considered as aspirin sensitive. Peripheral blood was collected to quantify mRNA expression, using microarray technology and quantitative RT-PCR.

RESULTS

Data indicated that AIA and ATA share large number of similarities for clinical phenotype. Screening of mRNA expression using microarray showed an overexpression of galectin-10 mRNA in AIA (AIA/ATA ratio = 1.9, P < 0.05). Results were confirmed using qRT-PCR. A positive correlation was established between microarray and qRT-PCR results for galectin-10 mRNA expression (r = 0.92, P < 0.0001). Finally, qRT-PCR results were validated on a subset of asthmatics and controls, showing an increased expression of galectin-10 mRNA in AIA vs ATA (P < 0.001) and vs controls (P < 0.01).

CONCLUSIONS

Our results demonstrate that AIA and ATA remain difficult to distinguish using clinical criteria. Employing two molecular biological methods, we demonstrate that galectin-10 mRNA is overexpressed in AIA, suggesting a novel candidate gene and a potentially innovative pathway for mucosal inflammation in aspirin intolerance.

Pre-processing of microarray data and analysis of differential expression

Microarrays have become a widely used technology in molecular biology research. One of their main uses is to measure gene expression. Compared to older expression measuring assays such as Northern blotting, analyzing gene expression data from microarrays is inherently more complex due to the massive amounts of data they produce. The analysis of microarray data requires biologists to collaborate with bioinformaticians or learn the basics of statistics and programming. Many software tools for microarray data analysis are available. Currently one of the most popular and freely available software tools is Bioconductor. This chapter uses Bioconductor to preprocess microarray data, detect differentially expressed genes, and annotate the gene lists of interest.

Large-scale mRNA expression profiling in the common ice plant, Mesembryanthemum crystallinum, performing C3 photosynthesis and Crassulacean acid metabolism (CAM)

The common ice plant (Mesembryanthemum crystallinum L.) has emerged as a useful model for molecular genetic studies of Crassulacean acid metabolism (CAM) because CAM can be induced in this species by water deficit or salinity stress. Non-redundant sequence information from expressed sequence tag data was used to fabricate a custom oligonucleotide microarray to compare large-scale mRNA expression patterns in M. crystallinum plants conducting C(3) photosynthesis versus CAM. Samples were collected every 4 h over a 24 h time period at the start of the subjective second day from plants grown under constant light and temperature conditions in order to capture variation in mRNA expression due to salinity stress and circadian clock control. Of 8455 genes, a total of 2343 genes (approximately 28%) showed a significant change as judged by analysis of variance (ANOVA) in steady-state mRNA abundance at one or more time points over the 24 h period. Of these, 858 (10%) and 599 (7%) exhibited a greater than two-fold ratio (TFR) increase or decrease in mRNA abundance, respectively. Functional categorization of these TFR genes revealed that many genes encoding products that function in CAM-related C(4) acid carboxylation/decarboxylation, glycolysis/gluconeogenesis, polysaccharide, polyol, and starch biosynthesis/degradation, protein degradation, transcriptional activation, signalling, stress response, and transport facilitation, and novel, unclassified proteins exhibited stress-induced increases in mRNA abundance. In contrast, salt stress resulted in a significant decrease in transcript abundance for genes encoding photosynthetic functions, protein synthesis, and cellular biogenesis functions. Many genes with CAM-related functions exhibited phase shifts in their putative circadian expression patterns following CAM induction. This report establishes an extensive catalogue of gene expression patterns for future investigations aimed at understanding the complex, transcriptional hierarchies that govern CAM-specific expression patterns. A novel graph-theoretic approach called 'Max Clique Builder' is introduced that identifies and organizes sets of coordinately regulated genes, such as those encoding subunits of the vacuolar H(+)-ATPase complex, into tighter functionally related clusters with more similar expression patterns compared with standard hierarchical clustering methods.

Normalization of microarray expression data using within-pedigree pool and its effect on linkage analysis

"Genetical genomics", the study of natural genetic variation combining data from genetic marker-based studies with gene expression analyses, has exploded with the recent development of advanced microarray technologies. To account for systematic variation known to exist in microarray data, it is critical to properly normalize gene expression traits before performing genetic linkage analyses. However, imposing equal means and variances across pedigrees can over-correct for the true biological variation by ignoring familial correlations in expression values. We applied the robust multiarray average (RMA) method to gene expression trait data from 14 Centre d'Etude du Polymorphisme Humain (CEPH) Utah pedigrees provided by GAW15 (Genetic Analysis Workshop 15). We compared the RMA normalization method using within-pedigree pools to RMA normalization using all individuals in a single pool, which ignores pedigree membership, and investigated the effects of these different methods on 18 gene expression traits previously found to be linked to regions containing the corresponding structural locus. Familial correlation coefficients of the expressed traits were stronger when traits were normalized within pedigrees. Surprisingly, the linkage plots for these traits were similar, suggesting that although heritability increases when traits are normalized within pedigrees, the strength of linkage evidence does not necessarily change substantially.

Map of differential transcript expression in the normal human large intestine

While there is considerable research related to using differential gene expression to predict disease phenotype classification, e.g., neoplastic tissue from nonneoplastic controls, there is little understanding of the range of expression in normal tissues. Understanding patterns of gene expression in nonneoplastic tissue, including regional anatomic expression changes within an organ, is vital to understanding gene expression changes in diseased tissue. To explore the gene expression change along the proximal-distal axis of the large intestine, we analyzed microarray data in 184 normal human specimens using univariate and multivariate techniques. We found 219 probe sets that were differentially expressed between the proximal and distal colorectal regions and 115 probe sets that were differentially expressed between the terminal segments, i.e., the cecum and rectum. We did not observe any probe sets that were statistically different between any two contiguous colorectal segments. The dominant expression pattern (65 probe sets) follows a dichotomous expression pattern consistent with the midgut-hindgut embryonic origins of the gut while a second pattern (50 probe sets) depicts a gradual change in transcript levels from the cecum to the rectum. While the dichotomous pattern includes roughly equal numbers of probe sets that are elevated proximally and distally, nearly all probe sets that show a gradual change demonstrate increasing expression levels moving from proximal to distal segments. These patterns describe an expression map of individual transcript variation as well as multigene expression patterns along the large intestine. This is the first gene expression map of an entire human organ.

Integration of expression profiles and genetic mapping data to identify candidate genes in intracranial aneurysm

Intracranial aneurysm (IA) is a complex genetic disease for which, to date, 10 loci have been identified by linkage. Identification of the risk-conferring genes in the loci has proven difficult, since the regions often contain several hundreds of genes. An approach to prioritize positional candidate genes for further studies is to use gene expression data from diseased and nondiseased tissue. Genes that are not expressed, either in diseased or nondiseased tissue, are ranked as unlikely to contribute to the disease. We demonstrate an approach for integrating expression and genetic mapping data to identify likely pathways involved in the pathogenesis of a disease. We used expression profiles for IAs and nonaneurysmal intracranial arteries (IVs) together with the 10 reported linkage intervals for IA. Expressed genes were analyzed for membership in Kyoto Encyclopedia of Genes and Genomes (KEGG) biological pathways. The 10 IA loci harbor 1,858 candidate genes, of which 1,561 (84%) were represented on the microarrays. We identified 810 positional candidate genes for IA that were expressed in IVs or IAs. Pathway information was available for 294 of these genes and involved 32 KEGG biological function pathways represented on at least 2 loci. A likelihood-based score was calculated to rank pathways for involvement in the pathogenesis of IA. Adherens junction, MAPK, and Notch signaling pathways ranked high. Integration of gene expression profiles with genetic mapping data for IA provides an approach to identify candidate genes that are more likely to function in the pathology of IA.

Hematopoietic progenitor cells (HPC) from mobilized peripheral blood display enhanced migration and marrow homing compared to steady-state bone marrow HPC

OBJECTIVE

Faster engraftment of G-CSF-mobilized peripheral blood (MPB) transplants compared to steady-state bone marrow (ssBM) is well documented and clinically relevant. A number of different factors likely contribute to this outcome. In the present study we explored whether independent of cell number there are intrinsic differences in the efficiency of progenitor cell homing to marrow between MPB and ssBM.

METHODS

Mobilization was achieved by continuous infusion of G-CSF alone or in combination with other mobilizing agents. In vivo homing assays, in vitro migration assays, gene expression analysis, and flow cytometry were utilized to compare homing-related in vivo and in vitro properties of MPB and ssBM HPC.

RESULTS

Marrow homing of murine MPB HPC, generated by different mobilizing schemes, was reproducibly significantly superior to that of ssBM, in lethally irradiated as well as in nonirradiated hosts. This phenotype was independent of MMP9, selectins, and beta2- and alpha4-integrins. Superior homing was also observed for human MPB HPC transplanted into NOD/SCIDbeta2microglobulin(-/-) recipients. Inhibition of HPC migration abrogated the homing advantage of MPB but did not affect homing of ssBM HPC, whereas enhancement of motility by CD26 inhibition improved marrow homing only of ssBM HPC. Enhanced SDF-1-dependent chemotaxis and low CD26 expression on MPB HPC were identified as potential contributing factors. Significant contributions of the putative alternative SDF-1 receptor, RDC1, were unlikely based on gene expression data.

CONCLUSION

The data suggest increased motility as a converging endpoint of complex changes seen in MPB HPC which is likely responsible for their favorable homing.

Oocyte regulation of metabolic cooperativity between mouse cumulus cells and oocytes: BMP15 and GDF9 control cholesterol biosynthesis in cumulus cells

Oocyte-derived bone morphogenetic protein 15 (BMP15) and growth differentiation factor 9 (GDF9) are key regulators of follicular development. Here we show that these factors control cumulus cell metabolism, particularly glycolysis and cholesterol biosynthesis before the preovulatory surge of luteinizing hormone. Transcripts encoding enzymes for cholesterol biosynthesis were downregulated in both Bmp15(-/-) and Bmp15(-/-) Gdf9(+/-) double mutant cumulus cells, and in wild-type cumulus cells after removal of oocytes from cumulus-cell-oocyte complexes. Similarly, cholesterol synthesized de novo was reduced in these cumulus cells. This indicates that oocytes regulate cumulus cell cholesterol biosynthesis by promoting the expression of relevant transcripts. Furthermore, in wild-type mice, Mvk, Pmvk, Fdps, Sqle, Cyp51, Sc4mol and Ebp, which encode enzymes required for cholesterol synthesis, were highly expressed in cumulus cells compared with oocytes; and oocytes, in the absence of the surrounding cumulus cells, synthesized barely detectable levels of cholesterol. Furthermore, coincident with reduced cholesterol synthesis in double mutant cumulus cells, lower levels were also detected in cumulus-cell-enclosed double mutant oocytes compared with wild-type oocytes. Levels of cholesterol synthesis in double mutant cumulus cells and oocytes were partially restored by co-culturing with wild-type oocytes. Together, these results indicate that mouse oocytes are deficient in synthesizing cholesterol and require cumulus cells to provide products of the cholesterol biosynthetic pathway. Therefore, oocyte-derived paracrine factors, particularly, BMP15 and GDF9, promote cholesterol biosynthesis in cumulus cells, probably as compensation for oocyte deficiencies in cholesterol production.

RReportGenerator: automatic reports from routine statistical analysis using R

With the establishment of high-throughput (HT) screening methods there is an increasing need for automatic analysis methods. Here we present RReportGenerator, a user-friendly portal for automatic routine analysis using the statistical platform R and Bioconductor. RReportGenerator is designed to analyze data using predefined analysis scenarios via a graphical user interface (GUI). A report in pdf format combining text, figures and tables is automatically generated and results may be exported. To demonstrate suitable analysis tasks we provide direct web access to a collection of analysis scenarios for summarizing data from transfected cell arrays (TCA), segmentation of CGH data, and microarray quality control and normalization.

AVAILABILITY

RReportGenerator, a user manual and a collection of analysis scenarios are available under a GNU public license on http://www-bio3d-igbmc.u-strasbg.fr/~wraff

Analysis of gene expression during neurite outgrowth and regeneration

Background

The ability of a neuron to regenerate functional connections after injury is influenced by both its intrinsic state and also by extrinsic cues in its surroundings. Investigations of the transcriptional changes undergone by neurons during in vivo models of injury and regeneration have revealed many transcripts associated with these processes. Because of the complex milieu of interactions in vivo, these results include not only expression changes directly related to regenerative outgrowth and but also unrelated responses to surrounding cells and signals. In vitro models of neurite outgrowth provide a means to study the intrinsic transcriptional patterns of neurite outgrowth in the absence of extensive extrinsic cues from nearby cells and tissues.

Results

We have undertaken a genome-wide study of transcriptional activity in embryonic superior cervical ganglia (SCG) and dorsal root ganglia (DRG) during a time course of neurite outgrowth in vitro. Gene expression observed in these models likely includes both developmental gene expression patterns and regenerative responses to axotomy, which occurs as the result of tissue dissection. Comparison across both models revealed many genes with similar gene expression patterns during neurite outgrowth. These patterns were minimally affected by exposure to the potent inhibitory cue Semaphorin3A, indicating that this extrinsic cue does not exert major effects at the level of nuclear transcription. We also compared our data to several published studies of DRG and SCG gene expression in animal models of regeneration, and found the expression of a large number of genes in common between neurite outgrowth in vitro and regeneration in vivo.

Conclusion

Many gene expression changes undergone by SCG and DRG during in vitro outgrowth are shared between these two tissue types and in common with in vivo regeneration models. This suggests that the genes identified in this in vitro study may represent new candidates worthy of further study for potential roles in the therapeutic regrowth of neuronal connections.

Transcriptomic and metabolite analyses of Cabernet Sauvignon grape berry development

BACKGROUND

Grape berry development is a dynamic process that involves a complex series of molecular genetic and biochemical changes divided into three major phases. During initial berry growth (Phase I), berry size increases along a sigmoidal growth curve due to cell division and subsequent cell expansion, and organic acids (mainly malate and tartrate), tannins, and hydroxycinnamates accumulate to peak levels. The second major phase (Phase II) is defined as a lag phase in which cell expansion ceases and sugars begin to accumulate. Véraison (the onset of ripening) marks the beginning of the third major phase (Phase III) in which berries undergo a second period of sigmoidal growth due to additional mesocarp cell expansion, accumulation of anthocyanin pigments for berry color, accumulation of volatile compounds for aroma, softening, peak accumulation of sugars (mainly glucose and fructose), and a decline in organic acid accumulation. In order to understand the transcriptional network responsible for controlling berry development, mRNA expression profiling was conducted on berries of V. vinifera Cabernet Sauvignon using the Affymetrix GeneChip Vitis oligonucleotide microarray ver. 1.0 spanning seven stages of berry development from small pea size berries (E-L stages 31 to 33 as defined by the modified E-L system), through véraison (E-L stages 34 and 35), to mature berries (E-L stages 36 and 38). Selected metabolites were profiled in parallel with mRNA expression profiling to understand the effect of transcriptional regulatory processes on specific metabolite production that ultimately influence the organoleptic properties of wine.

RESULTS

Over the course of berry development whole fruit tissues were found to express an average of 74.5% of probes represented on the Vitis microarray, which has 14,470 Unigenes. Approximately 60% of the expressed transcripts were differentially expressed between at least two out of the seven stages of berry development (28% of transcripts, 4,151 Unigenes, had pronounced (> or =2 fold) differences in mRNA expression) illustrating the dynamic nature of the developmental process. The subset of 4,151 Unigenes was split into twenty well-correlated expression profiles. Expression profile patterns included those with declining or increasing mRNA expression over the course of berry development as well as transient peak or trough patterns across various developmental stages as defined by the modified E-L system. These detailed surveys revealed the expression patterns for genes that play key functional roles in phytohormone biosynthesis and response, calcium sequestration, transport and signaling, cell wall metabolism mediating expansion, ripening, and softening, flavonoid metabolism and transport, organic and amino acid metabolism, hexose sugar and triose phosphate metabolism and transport, starch metabolism, photosynthesis, circadian cycles and pathogen resistance. In particular, mRNA expression patterns of transcription factors, abscisic acid (ABA) biosynthesis, and calcium signaling genes identified candidate factors likely to participate in the progression of key developmental events such as véraison and potential candidate genes associated with such processes as auxin partitioning within berry cells, aroma compound production, and pathway regulation and sequestration of flavonoid compounds. Finally, analysis of sugar metabolism gene expression patterns indicated the existence of an alternative pathway for glucose and triose phosphate production that is invoked from véraison to mature berries.

CONCLUSION

These results reveal the first high-resolution picture of the transcriptome dynamics that occur during seven stages of grape berry development. This work also establishes an extensive catalog of gene expression patterns for future investigations aimed at the dissection of the transcriptional regulatory hierarchies that govern berry development in a widely grown cultivar of wine grape. More importantly, this analysis identified a set of previously unknown genes potentially involved in critical steps associated with fruit development that can now be subjected to functional testing.

Dynamic changes of proteases and protease inhibitors revealed by microarray analysis in CA3 and entorhinal cortex during epileptogenesis in the rat

We investigated expression of genes involved in the proteolytic process during epileptogenesis in a rat model of temporal lobe epilepsy (TLE). In a previous microarray study we found prominent activation of this process, which reached highest expression during the acute and latent phase (1 week after SE) in CA3 and entorhinal cortex (EC). Detailed analysis shows differences in dynamics of the changes of several protease genes such as cathepsins, caspases, matrix metalloproteinases, and plasminogen activators. Most genes were acutely upregulated while others were mainly activated during the latent phase. Interestingly several proteolytic genes were still elevated in the chronic epileptic phase. Various protease inhibitors followed a similar time course. The identification of changes in the activation of genes involved in proteolysis at critical phases during epileptogenesis could point to potential time specific targets for intervention. The fact that several proteolytic genes were still activated in the chronic epileptic phase makes them interesting candidates to modify and slow down seizure progression.

Gene expression profile analysis of epilepsy-associated gangliogliomas

Gangliogliomas (GG) constitute the most frequent tumor entity in young patients undergoing surgery for intractable epilepsy. The histological composition of GG, with the presence of dysplastic neurons, corroborates their maldevelopmental origin. However, their histogenesis, the pathogenetic relationship with other developmental lesions, and the molecular alterations underlying the epileptogenicity of these tumors remain largely unknown. We performed gene expression analysis using the Affymetrix Gene Chip System (U133 plus 2.0 array). We used GENMAPP and the Gene Ontology database to identify global trends in gene expression data. Our analysis has identified various interesting genes and processes that are differentially expressed in GG when compared with normal tissue. The immune and inflammatory responses were the most prominent processes expressed in GG. Several genes involved in the complement pathway displayed a high level of expression compared with control expression levels. Higher expression was also observed for genes involved in cell adhesion, extracellular matrix and proliferation processes. We observed differential expression of genes as cyclin D1 and cyclin-dependent kinases, essential for neuronal cell cycle regulation and differentiation. Synaptic transmission, including GABA receptor signaling was an under-expressed process compared with control tissue. These data provide some suggestions for the molecular pathogenesis of GG. Furthermore, they indicate possible targets that may be investigated in order to dissect the mechanisms of epileptogenesis and possibly counteract the epileptogenic process in these developmental lesions.

Transforming properties of TC-1 in human breast cancer: interaction with FGFR2 and beta-catenin signaling pathways

Breast cancer development is associated with gene amplification and over expression that is believed to have a causative role in oncogenesis. Previous studies have demonstrated that over expression of TC-1(C8orf4) mRNA occurs in approximately 50% of breast cancer cell lines and primary tumor specimens. Here, we show that TC-1 has transforming properties in human mammary epithelial (HME) cells and its expression is mechanistically linked to FGFR signaling cascades. In vitro experiments demonstrate that TC-1 over expression mediates both anchorage-independent and growth factor-independent proliferation of HME cells. TC-1 was down regulated by the FGFR inhibitor PD173074 in the breast cancer cell line SUM-52 that also has an FGFR2 gene amplification and over expression. Furthermore, forced expression of FGFR2 in HME cells increased the level of expression of endogenous TC-1 mRNA. TC-1 has been implicated as a modulator of Wnt/beta-catenin signaling in 293 cells and in gastric cancer cells. However, while we did find increased expression of a subset of beta-catenin target genes in TC-1 over expressing cells, we did not find an association of TC-1 with global expression of beta-catenin target genes in our cells. Taken together, our data suggest that TC-1 over expression is transforming and may link with the FGFR pathway in a subset of breast cancer.

Transcriptional profiling of bipotential embryonic liver cells to identify liver progenitor cell surface markers

The ability to purify to homogeneity a population of hepatic progenitor cells from adult liver is critical for their characterization prior to any therapeutic application. As a step in this direction, we have used a bipotential liver cell line from 14 days postcoitum mouse embryonic liver to compile a list of cell surface markers expressed specifically by liver progenitor cells. These cells, known as bipotential mouse embryonic liver (BMEL) cells, proliferate in an undifferentiated state and are capable of differentiating into hepatocyte-like and cholangiocyte-like cells in vitro. Upon transplantation, BMEL cells are capable of differentiating into hepatocytes and cholangiocytes in vivo. Microarray and Gene Ontology (GO) analysis of gene expression in the 9A1 and 14B3 BMEL cell lines grown under proliferating and differentiating conditions was used to identify cell surface markers preferentially expressed in the bipotential undifferentiated state. This analysis revealed that proliferating BMEL cells express many genes involved in cell cycle regulation, whereas differentiation of BMEL cells by cell aggregation causes a switch in gene expression to functions characteristic of mature hepatocytes. In addition, microarray data and protein analysis indicated that the Notch signaling pathway could be involved in maintaining BMEL cells in an undifferentiated stem cell state. Using GO annotation, a list of cell surface markers preferentially expressed on undifferentiated BMEL cells was generated. One marker, Cd24a, is specifically expressed on progenitor oval cells in livers of diethyl 1,4-dihydro-2,4,6-trimethyl-3,5-pyridinedicarboxylate-treated animals. We therefore consider Cd24a expression a candidate molecule for purification of hepatic progenitor cells. Disclosure of potential conflicts of interest is found at the end of this article.

Adaptive gene expression divergence inferred from population genomics

Detailed studies of individual genes have shown that gene expression divergence often results from adaptive evolution of regulatory sequence. Genome-wide analyses, however, have yet to unite patterns of gene expression with polymorphism and divergence to infer population genetic mechanisms underlying expression evolution. Here, we combined genomic expression data--analyzed in a phylogenetic context--with whole genome light-shotgun sequence data from six Drosophila simulans lines and reference sequences from D. melanogaster and D. yakuba. These data allowed us to use molecular population genetics to test for neutral versus adaptive gene expression divergence on a genomic scale. We identified recent and recurrent adaptive evolution along the D. simulans lineage by contrasting sequence polymorphism within D. simulans to divergence from D. melanogaster and D. yakuba. Genes that evolved higher levels of expression in D. simulans have experienced adaptive evolution of the associated 3' flanking and amino acid sequence. Concomitantly, these genes are also decelerating in their rates of protein evolution, which is in agreement with the finding that highly expressed genes evolve slowly. Interestingly, adaptive evolution in 5' cis-regulatory regions did not correspond strongly with expression evolution. Our results provide a genomic view of the intimate link between selection acting on a phenotype and associated genic evolution.

Systematic order-dependent effect in expression values, variance, detection calls and differential expression in Affymetrix GeneChips

MOTIVATION

Affymetrix GeneChips are common 3' profiling platforms for quantifying gene expression. Using publicly available datasets of expression profiles from human and mouse experiments, we sought to characterize features of GeneChip data to better compare and evaluate analyses for differential expression, regulation and clustering. We uncovered an unexpected order dependence in expression data that holds across a variety of chips in both human and mouse data.

RESULTS

Order dependence among GeneChips affected relative expression measures pre-processed and normalized with the Affymetrix MAS5.0 algorithm and the robust multi-array average summarization method. The effect strongly influenced detection calls and tests for differential expression and can potentially significantly bias experimental results based on GeneChip profiling.

Trans-10, cis-12 conjugated linoleic acid activates the integrated stress response pathway in adipocytes

Trans-10, cis-12 conjugated linoleic acid (t10c12 CLA) causes fat loss in mouse white adipose tissue (WAT) and adipocytes in culture. The early transcriptome changes in treated WAT and 3T3-L1 adipocytes were analyzed using high-density microarrays to better characterize the signaling pathways responding to t10c12 CLA. Gene expression responses between 4 and 24 h after treatment showed a common set of early gene expression changes indicative of an integrated stress response (ISR). The responses of 3T3-L1 preadipocytes treated with t10c12 CLA or adipocytes treated with the cis-9, trans-11 isomer of CLA did not show the ISR, indicating the effect is specific to adipocytes responding to t10c12 CLA. Western blot analysis found increased phosphorylation of eIF2 alpha and increased production of ATF4 confirming at least part of the response to t10c12 CLA is mediated through the ISR pathway. Immunofluorescence microscopy found that the cell type expressing ATF3, an indicator of the ISR, was early stage adipocytes containing oil droplets but lacking the abundant levels of fatty acid binding protein-4 (FABP4) (AP2) found in mature adipocytes. Our data suggests that the ISR precedes and is possibly the cause of the later induction of proinflammatory cytokines observed in t10c12 CLA treated adipocytes. The release of proinflammatory cytokines may explain how the ISR in early stage adipocytes causes lipid loss in mature adipocytes.

Speed-mapping quantitative trait loci using microarrays

We developed a rapid, economical method for high-resolution quantitative trait locus (QTL) mapping using microarrays for selective genotyping of pooled DNA samples. We generated 21,207 F2 flies from two inbred Drosophila melanogaster strains with known QTLs affecting lifespan, and hybridized DNA pools of young and old flies to microarrays. We used changes of gene frequency of 2,326 single-feature polymorphisms (SFPs) to map previously identified and additional QTLs affecting lifespan.

CD154 tone sets the signaling pathways and transcriptome generated in model CD40-pluricompetent L3055 Burkitt's lymphoma cells

Activated B cells reacting to small amounts of CD40L (CD154) maintain homeostasis by suppressing default apoptosis. Additional outcomes, particularly differentiation, demand higher CD40 occupancy. Here, focusing on survival, we compared changes in the transcriptome of pleiotropically competent, early passage L3055 Burkitt's lymphoma cells confronted with low (picomolar) and high (nanomolar) concentrations of CD154 to gain insight into how a single receptor sets these distinct phenotypes. Of 267 genes altering transcriptional activity in response to strong CD154 tone, only 25 changed coordinately on low receptor occupancy. Seven of the top nine common up-regulated genes were targets of NF-kappaB. Direct measurement and functional inhibition of the NF-kappaB pathway revealed it to be central to a CD40-dependent survival signature. Although the canonical NF-kappaB axis was engaged by both signaling strengths equally, robust alternative pathway activation was a feature selective to a strong CD40 signal. Discriminatory exploitation of the two separate arms of NF-kappaB activation may indicate a principle whereby a cell senses and reacts differentially to shifting ligand availability. Identifying components selectively coupling CD40 to each axis could indicate targets for disruption in B cell pathologies underpinned by ectopic and/or hyper-CD154 activity such as neoplasia and some autoimmunities.

Comparative analysis of microarray normalization procedures: effects on reverse engineering gene networks

MOTIVATION

An increasingly common application of gene expression profile data is the reverse engineering of cellular networks. However, common procedures to normalize expression profiles generated using the Affymetrix GeneChips technology were originally developed for a rather different purpose, namely the accurate measure of differential gene expression between two or more phenotypes. As a result, current evaluation strategies lack comprehensive metrics to assess the suitability of available normalization procedures for reverse engineering and, in general, for measuring correlation between the expression profiles of a gene pair.

RESULTS

We benchmark four commonly used normalization procedures (MAS5, RMA, GCRMA and Li-Wong) in the context of established algorithms for the reverse engineering of protein-protein and protein-DNA interactions. Replicate sample, randomized and human B-cell data sets are used as an input. Surprisingly, our study suggests that MAS5 provides the most faithful cellular network reconstruction. Furthermore, we identify a crucial step in GCRMA responsible for introducing severe artifacts in the data leading to a systematic overestimate of pairwise correlation. This has key implications not only for reverse engineering but also for other methods, such as hierarchical clustering, relying on accurate measurements of pairwise expression profile correlation. We propose an alternative implementation to eliminate such side effect.

Analysis of potential transcriptomic biomarkers for Huntington's disease in peripheral blood

Highly quantitative biomarkers of neurodegenerative disease remain an important need in the urgent quest for disease-modifying therapies. For Huntington's disease (HD), a genetic test is available (trait marker), but necessary state markers are still in development. In this report, we describe a large battery of transcriptomic tests explored as state biomarker candidates. In an attempt to exploit the known neuroinflammatory and transcriptional perturbations of disease, we measured relevant mRNAs in peripheral blood cells. The performance of these potential markers was weak overall, with only one mRNA, immediate early response 3 (IER3), showing a modest but significant increase of 32% in HD samples compared with controls. No statistically significant differences were found for any other mRNAs tested, including a panel of 12 RNA biomarkers identified in a previous report [Borovecki F, Lovrecic L, Zhou J, Jeong H, Then F, Rosas HD, Hersch SM, Hogarth P, Bouzou B, Jensen RV, et al. (2005) Proc Natl Acad Sci USA 102:11023-11028]. The present results may nonetheless inform the future design and testing of HD biomarker strategies.

Genome-wide expression profiling of the Arabidopsis female gametophyte identifies families of small, secreted proteins

The female gametophyte of flowering plants, the embryo sac, develops within the diploid (sporophytic) tissue of the ovule. While embryo sac–expressed genes are known to be required at multiple stages of the fertilization process, the set of embryo sac–expressed genes has remained poorly defined. In particular, the set of genes responsible for mediating intracellular communication between the embryo sac and the male gametophyte, the pollen grain, is unknown. We used high-throughput cDNA sequencing and whole-genome tiling arrays to compare gene expression in wild-type ovules to that in dif1 ovules, which entirely lack embryo sacs, and myb98 ovules, which are impaired in pollen tube attraction. We identified nearly 400 genes that are downregulated in dif1 ovules. Seventy-eight percent of these embryo sac–dependent genes were predicted to encode for secreted proteins, and 60% belonged to multigenic families. Our results define a large number of candidate extracellular signaling molecules that may act during embryo sac development or fertilization; less than half of these are represented on the widely used ATH1 expression array. In particular, we found that 37 out of 40 genes encoding Domain of Unknown Function 784 (DUF784) domains require the synergid-specific transcription factor MYB98 for expression. Several DUF784 genes were transcribed in synergid cells of the embryo sac, implicating the DUF784 gene family in mediating late stages of embryo sac development or interactions with pollen tubes. The coexpression of highly similar proteins suggests a high degree of functional redundancy among embryo sac genes.

During the sexual reproduction of flowering plants, a pollen tube delivers sperm cells to a specialized group of cells known as the embryo sac, which contains the egg cell. It is known that embryo sacs are active participants in guiding the growth of pollen tubes, in facilitating fertilization, and in initiating seed development. However, the genes responsible for the complex biology of embryo sacs are poorly understood. The authors use two recently developed technologies, whole-genome tiling microarrays and high-throughput cDNA sequencing, to identify hundreds of genes expressed in embryo sacs of Arabidopsis thaliana. Most embryo sac–dependent genes have no known function, and include entire families of related genes that are only expressed in embryo sacs. Furthermore, most embryo sac–dependent genes encode small proteins that are potentially secreted from their cells of origin, suggesting that they may act as intracellular signals or to modify the extracellular matrix during fertilization or embryo sac development. These results illustrate the extent to which our understanding of plant sexual reproduction is limited and identifies hundreds of candidate genes for future studies investigating the molecular biology of the embryo sac.

Macromolecule biosynthesis: a key function of sleep

The function(s) of sleep remains a major unanswered question in biology. We assessed changes in gene expression in the mouse cerebral cortex and hypothalamus following different durations of sleep and periods of sleep deprivation. There were significant differences in gene expression between behavioral states; we identified 3,988 genes in the cerebral cortex and 823 genes in the hypothalamus with altered expression patterns between sleep and sleep deprivation. Changes in the steady-state level of transcripts for various genes are remarkably common during sleep, as 2,090 genes in the cerebral cortex and 409 genes in the hypothalamus were defined as sleep specific and changed (increased or decreased) their expression during sleep. The largest categories of overrepresented genes increasing expression with sleep were those involved in biosynthesis and transport. In both the cerebral cortex and hypothalamus, during sleep there was upregulation of multiple genes encoding various enzymes involved in cholesterol synthesis, as well as proteins for lipid transport. There was also upregulation during sleep of genes involved in synthesis of proteins, heme, and maintenance of vesicle pools, as well as antioxidant enzymes and genes encoding proteins of energy-regulating pathways. We postulate that during sleep there is a rebuilding of multiple key cellular components in preparation for subsequent wakefulness.

Circulating blood leukocyte gene expression profiles: effects of the Ames dwarf mutation on pathways related to immunity and inflammation

Aging is associated with a decline of immune competence and an increase in markers of inflammation. There is considerable evidence that inflammatory processes play a role in aging and the determination of lifespan. Hypopituitary Ames dwarf mice have extended longevity and exhibit many symptoms of delayed aging, although various aspects of immune function are suppressed in the mutants. In the present study, the expression of genes related to immunity and inflammation was compared in peripheral blood leukocytes (PBL) from Ames dwarf and normal mice using Affymetrix GeneChip arrays. Among the more than 3000 probe sets that were differentially expressed, 273 were identified as being associated with immunity and/or inflammation. Pathway analysis revealed interactions among 91 of these probe sets, centered on casp3, bcl2, il4, prkca, mapk14 and TGFbeta1. Ames dwarf mice had reduced leukocyte expression of casp3 and TGFbeta and increased expression of Bcl2. Alterations in the expression of these genes suggest likely functional changes in apoptosis, B and T cell homeostasis, prostaglandin synthesis, humoral immunity, chemokine activity, complement activation, hemostasis and wound healing pathways. Collectively, these results suggest that activation of both anti-inflammatory pathways and an anti-clotting mechanism combined with reduced turnover of leukocytes may contribute to delayed aging and extended longevity of Ames dwarf mice. We are also aware that alterations in gene expression in PBLs can be due to different composition of PBL populations when comparing Ames dwarf to WT animals, and it will be interesting to investigate these genes in particular PBL populations in the future. However, whole leukocytes population represents the function of immune system in these organisms.

Mutant huntingtin's effects on striatal gene expression in mice recapitulate changes observed in human Huntington's disease brain and do not differ with mutant huntingtin length or wild-type huntingtin dosage

To test the hypotheses that mutant huntingtin protein length and wild-type huntingtin dosage have important effects on disease-related transcriptional dysfunction, we compared the changes in mRNA in seven genetic mouse models of Huntington's disease (HD) and postmortem human HD caudate. Transgenic models expressing short N-terminal fragments of mutant huntingtin (R6/1 and R6/2 mice) exhibited the most rapid effects on gene expression, consistent with previous studies. Although changes in the brains of knock-in and full-length transgenic models of HD took longer to appear, 15- and 22-month CHL2(Q150/Q150), 18-month Hdh(Q92/Q92) and 2-year-old YAC128 animals also exhibited significant HD-like mRNA signatures. Whereas it was expected that the expression of full-length huntingtin transprotein might result in unique gene expression changes compared with those caused by the expression of an N-terminal huntingtin fragment, no discernable differences between full-length and fragment models were detected. In addition, very high correlations between the signatures of mice expressing normal levels of wild-type huntingtin and mice in which the wild-type protein is absent suggest a limited effect of the wild-type protein to change basal gene expression or to influence the qualitative disease-related effect of mutant huntingtin. The combined analysis of mouse and human HD transcriptomes provides important temporal and mechanistic insights into the process by which mutant huntingtin kills striatal neurons. In addition, the discovery that several available lines of HD mice faithfully recapitulate the gene expression signature of the human disorder provides a novel aspect of validation with respect to their use in preclinical therapeutic trials.

The utility of MAS5 expression summary and detection call algorithms

BACKGROUND

Used alone, the MAS5.0 algorithm for generating expression summaries has been criticized for high False Positive rates resulting from exaggerated variance at low intensities.

RESULTS

Here we show, with replicated cell line data, that, when used alongside detection calls, MAS5 can be both selective and sensitive. A set of differentially expressed transcripts were identified that were found to be changing by MAS5, but unchanging by RMA and GCRMA. Subsequent analysis by real time PCR confirmed these changes. In addition, with the Latin square datasets often used to assess expression summary algorithms, filtered MAS5.0 was found to have performance approaching that of its peers.

CONCLUSION

When used alongside detection calls, MAS5 is a sensitive and selective algorithm for identifying differentially expressed genes.

Methodological considerations for gene expression profiling of human brain

Gene expression profiles of postmortem brain tissue represent important resources for understanding neuropsychiatric illnesses. The impact(s) of quality covariables on the analysis and results of gene expression studies are important questions. This paper addressed critical variables which might affect gene expression in two brain regions. Four broad groups of quality indicators in gene expression profiling studies (clinical, tissue, RNA, and microarray quality) were identified. These quality control indicators were significantly correlated, however one quality variable did not account for the total variance in microarray gene expression. The data showed that agonal factors and low pH correlated with decreased integrity of extracted RNA in two brain regions. These three parameters also modulated the significance of alterations in mitochondrial-related genes. The average F-ratio summaries across all transcripts showed that RNA degradation from the AffyRNAdeg program accounted for higher variation than all other quality factors. Taken together, these findings confirmed prior studies, which indicated that quality parameters including RNA integrity, agonal factors, and pH are related to differences in gene expression profiles in postmortem brain. Individual candidate genes can be evaluated with these quality parameters in post hoc analysis to help strengthen the relevance to psychiatric disorders. We find that clinical, tissue, RNA, and microarray quality are all useful variables for collection and consideration in study design, analysis, and interpretation of gene expression results in human postmortem studies.

Functional Role for IκBNS in T Cell Cytokine Regulation As Revealed by Targeted Gene Disruption

Triggering of the TCR by cognate peptide/MHC ligands induces expression of I kappa BNS, a member of the I kappa B family of NF-kappaB inhibitors whose expression is associated with apoptosis of immature thymocytes. To understand the role of I kappa BNS in TCR triggering, we created a targeted disruption of the I kappa BNS gene. Surprisingly, mice lacking I kappa BNS show normal thymic progression but both thymocytes and T cells manifest reduced TCR-stimulated proliferation. Moreover, I kappa BNS knockout thymocytes and T cells produce significantly less IL-2 and IFN-gamma than wild-type cells. Transfection analysis demonstrates that I kappa BNS and c-Rel individually increase IL-2 promoter activity. The effect of I kappa BNS on the IL-2 promoter, unlike c-Rel, is dependent on the NF-kappaB rather than the CD28RE site; mutation of the NF-kappaB site extinguishes the induction of transcription by I kappa BNS in transfectants and prevents association of I kappa BNS with IL-2 promoter DNA. Microarray analyses confirm the reduction in IL-2 production and some IFN-gamma-linked transcripts in I kappa BNS knockout T cells. Collectively, our findings demonstrate that I kappa BNS regulates production of IL-2 and other cytokines induced via "strong" TCR ligation.

Altered phenotype and gene transcription in endothelial cells, induced by Plasmodium falciparum-infected red blood cells: pathogenic or protective?

Severe malaria is associated with sequestration of Plasmodium falciparum-infected red blood cells (PRBC) in the microvasculature and elevation of intercellular adhesion molecule-1 (ICAM-1) and TNF. In vitro co-culture of human umbilical vein endothelial cells (HUVEC), with either PRBC or uninfected RBC, required the presence of low level TNF (5pg/ml) for significant up-regulation of ICAM-1, which may contribute to increased cytoadhesion in vivo. These effects were independent of P. falciparum erythrocyte membrane protein-1 (PfEMP-1)-mediated adhesion but critically dependent on cell-cell contact. Further changes included increases in IL8 release and soluble TNF receptor shedding. Microarray analysis of HUVEC transcriptome following co-culture, using a human Affymetrix microarray chip, showed significant differential regulation of genes which defined gene ontologies such as cell communication, cell adhesion, signal transduction and immune response. Our data demonstrate that endothelial cells have the ability to mobilise immune and pro-adhesive responses when exposed to both PRBC and TNF. In addition, there is also a previously un-described positive regulation by RBC and TNF and a concurrent negative regulation of a range of genes involved in inflammation and cell-death, by PRBC and TNF. We propose that the balance between positive and negative regulation demonstrated in our study will determine endothelial pathology during a malaria infection.

Growth hormone regulation of metabolic gene expression in muscle: a microarray study in hypopituitary men

Muscle is a target of growth hormone (GH) action and a major contributor to whole body metabolism. Little is known about how GH regulates metabolic processes in muscle or the extent to which muscle contributes to changes in whole body substrate metabolism during GH treatment. To identify GH-responsive genes that regulate substrate metabolism in muscle, we studied six hypopituitary men who underwent whole body metabolic measurement and skeletal muscle biopsies before and after 2 wk of GH treatment (0.5 mg/day). Transcript profiles of four subjects were analyzed using Affymetrix GeneChips. Serum insulin-like growth factor I (IGF-I) and procollagens I and III were measured by RIA. GH increased serum IGF-I and procollagens I and III, enhanced whole body lipid oxidation, reduced carbohydrate oxidation, and stimulated protein synthesis. It induced gene expression of IGF-I and collagens in muscle. GH reduced expression of several enzymes regulating lipid oxidation and energy production. It reduced calpain 3, increased ribosomal protein L38 expression, and displayed mixed effects on genes encoding myofibrillar proteins. It increased expression of circadian gene CLOCK, and reduced that of PERIOD. In summary, GH exerted concordant effects on muscle expression and blood levels of IGF-I and collagens. It induced changes in genes regulating protein metabolism in parallel with a whole body anabolic effect. The discordance between muscle gene expression profiles and metabolic responses suggests that muscle is unlikely to contribute to GH-induced stimulation of whole body energy and lipid metabolism. GH may regulate circadian function in skeletal muscle by modulating circadian gene expression with possible metabolic consequences.

Genome-enabled hitchhiking mapping identifies QTLs for stress resistance in natural Drosophila

Identification of genes underlying complex traits is an important problem. Quantitative trait loci (QTL) are mapped using marker-trait co-segregation in large panels of recombinant genotypes. Most frequently, recombinant inbred lines derived from two isogenic parents are used. Segregation patterns are also studied in pedigrees from multiple families. Great advances have been made through creative use of these techniques, but narrow sampling and inadequate power represent strong limitations. Here, we propose an approach combining the strengths of both techniques. We established a mapping population from a sample of natural genotypes, and applied artificial selection for a complex character. Selection changed the frequencies of alleles in QTLs contributing to the selection response. We infer QTLs with dense genotyping microarrays by identifying blocks of linked markers undergoing selective changes in allele frequency. We demonstrated this approach with an experimental population composed from 20 isogenic strains. Selection for starvation survival was executed in three replicated populations with three control non-selected populations. Three individuals per population were genotyped using Affymetrix GeneChips. Two regions of the genome, one each on the left arms of the second and third chromosomes, showed significant divergence between control and selected populations. For the former region, we inferred allele frequencies in selected and control populations by pyrosequencing. We conclude that the allele frequency difference, averaging approximately 40% between selected and control lines, contributed to selection response. Our approach can contribute to the fine scale decomposition of the genetics of direct and indirect selection responses, and genotype by environment interactions.

beadarray: R classes and methods for Illumina bead-based data

The R/Bioconductor package beadarray allows raw data from Illumina experiments to be read and stored in convenient R classes. Users are free to choose between various methods of image processing, background correction and normalization in their analysis rather than using the defaults in Illumina's; proprietary software. The package also allows quality assessment to be carried out on the raw data. The data can then be summarized and stored in a format which can be used by other R/Bioconductor packages to perform downstream analyses. Summarized data processed by Illumina's; BeadStudio software can also be read and analysed in the same manner.

AVAILABILITY

The beadarray package is available from the Bioconductor web page at www.bioconductor.org. A user's guide and example data sets are provided with the package.

Tissue-specific mRNA expression profiling in grape berry tissues

Background

Berries of grape (Vitis vinifera) contain three major tissue types (skin, pulp and seed) all of which contribute to the aroma, color, and flavor characters of wine. The pericarp, which is composed of the exocarp (skin) and mesocarp (pulp), not only functions to protect and feed the developing seed, but also to assist in the dispersal of the mature seed by avian and mammalian vectors. The skin provides volatile and nonvolatile aroma and color compounds, the pulp contributes organic acids and sugars, and the seeds provide condensed tannins, all of which are important to the formation of organoleptic characteristics of wine. In order to understand the transcriptional network responsible for controlling tissue-specific mRNA expression patterns, mRNA expression profiling was conducted on each tissue of mature berries of V. vinifera Cabernet Sauvignon using the Affymetrix GeneChip^® Vitis oligonucleotide microarray ver. 1.0. In order to monitor the influence of water-deficit stress on tissue-specific expression patterns, mRNA expression profiles were also compared from mature berries harvested from vines subjected to well-watered or water-deficit conditions.

Results

Overall, berry tissues were found to express approximately 76% of genes represented on the Vitis microarray. Approximately 60% of these genes exhibited significant differential expression in one or more of the three major tissue types with more than 28% of genes showing pronounced (2-fold or greater) differences in mRNA expression. The largest difference in tissue-specific expression was observed between the seed and pulp/skin. Exocarp tissue, which is involved in pathogen defense and pigment production, showed higher mRNA abundance relative to other berry tissues for genes involved with flavonoid biosynthesis, pathogen resistance, and cell wall modification. Mesocarp tissue, which is considered a nutritive tissue, exhibited a higher mRNA abundance of genes involved in cell wall function and transport processes. Seeds, which supply essential resources for embryo development, showed higher mRNA abundance of genes encoding phenylpropanoid biosynthetic enzymes, seed storage proteins, and late embryogenesis abundant proteins. Water-deficit stress affected the mRNA abundance of 13% of the genes with differential expression patterns occurring mainly in the pulp and skin. In pulp and seed tissues transcript abundance in most functional categories declined in water-deficit stressed vines relative to well-watered vines with transcripts for storage proteins and novel (no-hit) functional assignments being over represented. In the skin of berries from water-deficit stressed vines, however, transcripts from several functional categories including general phenypropanoid and ethylene metabolism, pathogenesis-related responses, energy, and interaction with the environment were significantly over-represented.

Conclusion

These results revealed novel insights into the tissue-specific expression mRNA expression patterns of an extensive repertoire of genes expressed in berry tissues. This work also establishes an extensive catalogue of gene expression patterns for future investigations aimed at the dissection of the transcriptional regulatory hierarchies that govern tissue-specific expression patterns associated with tissue differentiation within berries. These results also confirmed that water-deficit stress has a profound effect on mRNA expression patterns particularly associated with the biosynthesis of aroma and color metabolites within skin and pulp tissues that ultimately impact wine quality.

Monitoring gene changes during antiepileptic drug therapy to widen the safety window and reduce pharmacoresistance

A novel and potentially effective way to decrease pharmacoresistance is to widen an antiepileptic drug's (AED) safety window by altering its efficacy-toxicity relationship to reduce the frequency or eliminate the occurrence of drug-induced toxicities. The first step in this process is to identify specific AED toxicities that are common, objective, reliable to assess, and exhibit both interpatient and interdrug variability. Next, one or more of the fundamental sources of biological variability (e.g., DNA, RNA, proteins, or metabolites) is selected for study. By identifying a relationship between variations in these fundamental sources and specific AED toxicities, one aims to identify pathways underlying the genesis of the toxicity (mechanisms) or biomarkers related to elevated risk for AED toxicity (predictors). These types of studies have many methodological challenges. However, if properly conducted, the resulting data could be used to design interventions that directly block the toxicity mechanisms or preselect individuals who are at high risk of developing these toxicities. Both approaches would conceivably widen the safety window of efficacious AEDs and reduce the incidence of pharmacoresistance. This article describes a general framework for this approach and provides an example of the methodology using differential gene expression measured using RNA-based microarray technology in children experiencing VPA-associated weight gain.

New candidate genes for sex-comb divergence between Drosophila mauritiana and Drosophila simulans

A large-effect QTL for divergence in sex-comb tooth number between Drosophila simulans and D. mauritiana was previously mapped to 73A-84AB. Here we identify genes that are likely contributors to this divergence. We first improved the mapping resolution in the 73A-84AB region using 12 introgression lines and 62 recombinant nearly isogenic lines. To further narrow the list of candidate genes, we assayed leg-specific expression and identified genes with transcript-level evolution consistent with a potential role in sex-comb divergence. Sex combs are formed on the prothoracic (front) legs, but not on the mesothoracic (middle) legs of Drosophila males. We extracted RNA from the prothoracic and mesothoracic pupal legs of two species to determine which of the genes expressed differently between leg types were also divergent for gene expression. Two good functional candidate genes, Scr and dsx, are located in one of our fine-scale QTL regions. In addition, three previously uncharacterized genes (CG15186, CG2016, and CG2791) emerged as new candidates. These genes are located in regions strongly associated with sex-comb tooth number differences and are expressed differently between leg tissues and between species. Further supporting the potential involvement of these genes in sex-comb divergence, we found a significant difference in sex-comb tooth number between co-isogenic D. melanogaster lines with and without P-element insertions at CG2791.

A set of stage-specific gene transcripts identified in EK stage X and HH stage 3 chick embryos

Background

The embryonic developmental process in avian species is quite different from that in mammals. The first cleavage begins 4 h after fertilization, but the first differentiation does not occur until laying of the egg (Eyal-Giladi and Kochav (EK) stage X). After 12 to 13 h of incubation (Hamburger and Hamilton (HH) stage 3), the three germ layers form and germ cell segregation in the early chick embryo are completed. Thus, to identify genes associated with early embryonic development, we compared transcript expression patterns between undifferentiated (stage X) and differentiated (HH stage 3) embryos.

Results

Microarray analysis primarily showed 40 genes indicating the significant changes in expression levels between stage X and HH stage 3, and 80% of the genes (32/40) were differentially expressed with more than a twofold change. Among those, 72% (23/32) were relatively up-regulated at stage X compared to HH stage 3, while 28% (9/32) were relatively up-regulated at HH stage 3 compared to stage X. Verification and gene expression profiling of these GeneChip expression data were performed using quantitative RT-PCR for 32 genes at developmental four points; stage X (0 h), HH stage 3 (12 h), HH stage 6 (24 h), and HH stage 9 (30 h). Additionally, we further analyzed four genes with less than twofold expression increase at HH stage 3. As a result, we identified a set of stage-specific genes during the early chick embryo development; 21 genes were relatively up-regulated in the stage X embryo and 12 genes were relatively up-regulated in the HH stage 3 embryo based on both results of microarray and quantitative RT-PCR.

Conclusion

We identified a set of genes with stage-specific expression from microarray Genechip and quantitative RT-PCR. Discovering stage-specific genes will aid in uncovering the molecular mechanisms involved the formation of the three germ layers and germ cell segregation in the early chick embryos.

Characterization of the RSL1-dependent conditional expression system in LNCaP prostate cancer cells and development of a single vector format

BACKGROUND

Conditional expression systems are useful tools for the study of gene function but the use of these systems in prostate cancer cells is limited by the undesired biological effects of the inducing ligands. The RheoSwitch system employs RheoSwitch Ligand 1 (RSL1), a non-steroidal analog of the insect hormone ecdysone, to activate a modified nuclear receptor heterodimer that controls target gene expression via GAL4 response elements. This system has not been tested in prostate cancer cells.

METHODS

We established LNCaP human prostate cancer cell lines that constitutively express RheoSwitch transcription factors to quantify RSL1-dependent expression and assess the effects of RSL1 on cell proliferation and endogenous gene expression. Potential RSL1-responsive genes were identified using Affymetrix microarrays and validated by Northern blot hybridization. A single-vector format was developed to establish cell lines that conditionally produce a recombinant protein.

RESULTS

Stable cell lines displayed tight and potent (over several orders of magnitude) RSL1-dependent regulation of a transiently transfected luciferase reporter gene. RSL1 did not affect basal or androgen-stimulated cell proliferation and exerted minimal effects on the expression of endogenous genes. Cell lines established using the single-vector system also displayed strictly RSL1-dependent production of the recombinant protein encoded by the stably integrated RSL1-responsive expression cassette.

CONCLUSIONS

The RheoSwitch system is well suited for conditional gene expression in prostate cancer cells. The single-vector format should facilitate the production of stable cell lines. This system should be useful for the study of proteins involved in prostate cancer in both cell and animal models of the disease.

Early blood profiles of virus infection in a monkey model for Lassa fever

Acute arenavirus disease in primates, like Lassa hemorrhagic fever in humans, begins with flu-like symptoms and leads to death approximately 2 weeks after infection. Our goal was to identify molecular changes in blood that are related to disease progression. Rhesus macaques (Macaca mulatta) infected intravenously with a lethal dose of lymphocytic choriomeningitis virus (LCMV) provide a model for Lassa virus infection of humans. Blood samples taken before and during the course of infection were used to monitor gene expression changes that paralleled disease onset. Changes in blood showed major disruptions in eicosanoid, immune response, and hormone response pathways. Approximately 12% of host genes alter their expression after LCMV infection, and a subset of these genes can discriminate between virulent and non-virulent LCMV infection. Major transcription changes have been given preliminary confirmation by quantitative PCR and protein studies and will be valuable candidates for future validation as biomarkers for arenavirus disease.

Analysis of cell death inducing compounds

Biomarkers for early detection of toxicity hold the promise of improving the failure rates in drug development. In the present study, gene expression levels were measured using full-genome RAE230 version 2 Affymetrix GeneChips on rat liver tissue 48 h after administration of six different compounds, three toxins (ANIT, DMN and NMF) and three non-toxins (Caeruelein, Dinitrophenol and Rosiglitazone). We identified three gene transcripts with exceptional predictive performance towards liver toxicity and/or changes in histopathology. The three genes were: glucokinase regulatory protein (GCKR), ornithine aminotransferase (OAT) and Cytochrome P450, subfamily IIC (mephenytoin 4-hydroxylase) (Cyp2C29). RT-PCR for these three genes was performed and four additional compounds were included for validation. The quantitative RT-PCR analysis confirmed the findings based on the microarray data and using the three genes a classification rate of 55 of 57 samples was achieved for the classification of not toxic versus toxic. The single most promising biomarker (OAT) alone resulted in a surprisingly 100% correctly classified samples. OAT has not previously been linked to toxicity and cell death in the literature and the novel finding represents a putative hepatotoxicity biomarker.