Importance of replication in microarray gene expression studies: Statistical methods and evidence from repetitive cDNA hybridizations

We present statistical methods for analyzing replicated cDNA microarray expression data and report the results of a controlled experiment. The study was conducted to investigate inherent variability in gene expression data and the extent to which replication in an experiment produces more consistent and reliable findings. We introduce a statistical model to describe the probability that mRNA is contained in the target sample tissue, converted to probe, and ultimately detected on the slide. We also introduce a method to analyze the combined data from all replicates. Of the 288 genes considered in this controlled experiment, 32 would be expected to produce strong hybridization signals because of the known presence of repetitive sequences within them. Results based on individual replicates, however, show that there are 55, 36, and 58 highly expressed genes in replicates 1, 2, and 3, respectively. On the other hand, an analysis by using the combined data from all 3 replicates reveals that only 2 of the 288 genes are incorrectly classified as expressed. Our experiment shows that any single microarray output is subject to substantial variability. By pooling data from replicates, we can provide a more reliable analysis of gene expression data. Therefore, we conclude that designing experiments with replications will greatly reduce misclassification rates. We recommend that at least three replicates be used in designing experiments by using cDNA microarrays, particularly when gene expression data from single specimens are being analyzed.

Expression profiling and sequence diversity of novel DREB genes from common bean (Phaseolus vulgaris L.) and their association with drought-related traits

Common bean is a major dietary component in several countries, but its productivity is negatively affected by abiotic stresses. Dissecting candidate genes involved in abiotic stress tolerance is a paramount step toward the improvement of common bean performance under such constraints. Thereby, this thesis presents a systematic analysis of the DEHYDRATION RESPONSIVE ELEMENT-BINDING (DREB) gene subfamily, which encompasses genes that regulate several processes during stress responses, but with limited information for common bean. First, a series of in silico analyses with sequences retrieved from the P. vulgaris genome on Phytozome supported the categorization of 54 putative PvDREB genes distributed within six phylogenetic subgroups (A-1 to A-6), along the 11 chromosomes. Second, we cloned four novel PvDREB genes and determined their inducibility-factors, including the dehydration-, salinity- and cold-inducible genes PvDREB1F and PvDREB5A, and the dehydration- and cold-inducible genes PvDREB2A and PvDREB6B. Afterwards, nucleotide polymorphisms were searched through Sanger sequencing along those genes, revealing a high number of single nucleotide polymorphisms within PvDREB6B by the comparison of Mesoamerican and Andean genotypes. The nomenclature of PvDREB6B is discussed in details. Furthermore, we used the BARCBean6K_3 SNP platform to identify and genotype the closest SNP to each one of the 54 PvDREB genes. We selected PvDREB6B for a broader study encompassing a collection of wild common bean accessions of Mesoamerican origin. The population structure of the wild beans was accessed using sequence polymorphisms of PvDREB6B. The genetic clusters were partially associated with variation in latitude, altitude, precipitation and temperature throughout the areas such beans are distributed. With an emphasis on drought stress, an adapted tube-screening method in greenhouse conditions enabled the phenotyping of several drought-related traits in the wild collection. Interestingly, our data revealed a correlation between root depth, plant height and biomass and the environmental data of the location of the accessions. Correlation was also observed between the population structure determined through PvDREB6B and the environmental data. An association study combining data from the SNP array and DREB polymorphisms enabled the detection of SNP associated with drought-related traits through a compressed mixed linear model (CMLM) analysis. This thesis highlighted important features of DREB genes in common bean, revealing candidates for further strategies aimed at improvement of abiotic stress tolerance, with emphasis on drought tolerance

Tissue-specific gene expression in soybean (Glycine max) detected by cDNA microarray analysis

We have constructed cDNA microarrays for soybean (Glycine max L. Merrill), containing approximately 4,100 Unigene ESTs derived from axenic roots, to evaluate their application and utility for functional genomics of organ differentiation in legumes. We assessed microarray technology by conducting studies to evaluate the accuracy of microarray data and have found them to be both reliable and reproducible in repeat hybridisations. Several ESTs showed high levels (>50 fold) of differential expression in either root or shoot tissue of soybean. A small number of physiologically interesting, and differentially expressed sequences found by microarray analysis were verified by both quantitative real-time RT-PCR and Northern blot analysis. There was a linear correlation (r(2) = 0.99, over 5 orders of magnitude) between microarray and quantitative real-time RT-PCR data. Microarray analysis of soybean has enormous potential not only for the discovery of new genes involved in tissue differentiation and function, but also to study the expression of previously characterised genes, gene networks and gene interactions in wild-type, mutant or transgenic; plants.

Effects of A1 adenosine receptor overexpression on normoxic and post-ischemic gene expression

Objectives: To identify potential molecular genetic determinants of cardiovascular ischemic tolerance in wild-type and transgenic hearts overexpressing A(1) adenosine receptors (A(1)ARs). Methods: cDNA microarrays were used to explore expression of 1824 genes ill wild-type hearts and ischemia-tolerant mouse hearts overexpressing A(1)ARs. Results: Overexpression of A(1)ARs reduced post-ischemic contractile dysfunction, limited arrhythmogenesis, and reduced necrosis by similar to80% in hearts subjected to 30 min global ischemia 60 mill reperfusion. Cardioprotection was abrogated by acute A(1)AR antagonism, and only a small number (19) of genes were modified by A(1)AR overexpression in normoxic hearts. Ischemia-reperfusion significantly altered expression of 75 genes in wild-type hearts (14 induced, 61 down-regulated), including genes for metabolic enzymes, structural/motility proteins, cell signaling proteins, defense/growth proteins, and regulators of transcription and translation. A(1)AR overexpression reversed the majority of gene down-regulation whereas gene induction was generally unaltered. Additionally, genes involved in cell defence, signaling and gene expression were selectively modified by ischemia in transgenic hearts (33 induced, 10 down-regulated), possibly contributing to the protected phenotype. Real-time PCR verified changes in nine selected genes, revealing concordance with array data. Transcription of the A(1)AR gene was also modestly reduced post-ischemia, consistent with impaired functional sensitivity to A(1)AR stimulation Conclusions: Data are presented regarding the early post-ischemic gene profile of intact heart. Reduced A(1)AR transcription is observed which may contribute to poor outcome from ischemia. A(1)AR overexpression selectively modifies post-ischemic gene expression, potentially contributing to ischemic-tolerance. (C) 2003 European Society of Cardiology. Published by Elsevier Science B.V. All rights reserved.

Effector ExoU from the type III secretion system is an important modulator of gene expression in lung epithelial cells in response to Pseudomonas aeruginosa infection

Pseudomonas aeruginosa is an important pathogen in immunocompromised patients and secretes a diverse set of virulence factors that aid colonization and influence host cell defenses. An important early step in the establishment of infection is the production of type III-secreted effectors translocated into host cells by the bacteria. We used cDNA microarrays to compare the transcriptomic response of lung epithelial cells to P. aeruginosa mutants defective in type IV pili, the type III secretion apparatus, or in the production of specific type III-secreted effectors. Of the 18,000 cDNA clones analyzed, 55 were induced or repressed after 4 It of infection and could be classified into four different expression patterns. These include (i) host genes that are induced or repressed in a type III secretion-independent manner (32 clones), (ii) host genes induced specifically by ExoU (20 clones), and (iii) host genes induced in an ExoU-independent but type III secretion dependent manner (3 clones). In particular, ExoU was essential for the expression of immediate-early response genes, including the transcription factor c-Fos. ExoU-dependent gene expression was mediated in part by early and transient activation of the AN transcription factor complex. In conclusion, the present study provides a detailed insight into the response of epithelial cells to infection and indicates the significant role played by the type III virulence mechanism in the initial host response.

Microarray expression profiling in melanoma reveals a BRAF mutation signature

We have used microarray gene expression pro. ling and machine learning to predict the presence of BRAF mutations in a panel of 61 melanoma cell lines. The BRAF gene was found to be mutated in 42 samples (69%) and intragenic mutations of the NRAS gene were detected in seven samples (11%). No cell line carried mutations of both genes. Using support vector machines, we have built a classifier that differentiates between melanoma cell lines based on BRAF mutation status. As few as 83 genes are able to discriminate between BRAF mutant and BRAF wild-type samples with clear separation observed using hierarchical clustering. Multidimensional scaling was used to visualize the relationship between a BRAF mutation signature and that of a generalized mitogen-activated protein kinase ( MAPK) activation ( either BRAF or NRAS mutation) in the context of the discriminating gene list. We observed that samples carrying NRAS mutations lie somewhere between those with or without BRAF mutations. These observations suggest that there are gene-specific mutation signals in addition to a common MAPK activation that result from the pleiotropic effects of either BRAF or NRAS on other signaling pathways, leading to measurably different transcriptional changes.

Expression profiling in alcoholism research

This article represents the proceedings of a symposium at the 2004 International Society for Biomedical Research on Alcoholism in Mannheim, Germany, organized and co-chaired by Susan E. Bergeson and Wolfgang Sommer. The presentations and presenter were (1) Gene Expression in Brains of AlcoholPreferring and Non-Preferring Rats, by Howard J. Edenberg (2) Candidate Treatment Targets for Alcoholism: Leads from Functional Genomics Approaches, by Wolfgang Sommer (3) Microarray Analysis of Acute and Chronic Alcohol Response in Brain, by Susan E. Bergeson (4) On the Integration of QTL and Gene Expression Analysis, by Robert J. Hitzemann (5) Microarray and Proteomic Analysis of the Human Alcoholic Brain, by Peter R. Dodd.

Differential gene expression in the developing mouse ureter

In many instances, kidney dysgenesis results as a secondary consequence to defects in the development of the ureter. Through the use of mouse genetics a number of genes associated with such malformations have been identified, however, the cause of many other abnormalities remain unknown. In order to identify novel genes involved in ureter development we compared gene expression in embryonic day (E) 12.5, E15.5 and postnatal day (P) 75 ureters using the Compugen mouse long oligo microarrays. A total of 248 genes were dynamically upregulated and 208 downregulated between E12.5 and P75. At E12.5, when the mouse ureter is comprised of a simple cuboidal epithelium surrounded by ureteric mesenchyme, genes previously reported to be expressed in the ureteric mesenchyme, foxC1 and foxC2 were upregulated. By E15.5 the epithelial layer develops into urothelium, impermeable to urine, and smooth muscle develops for the peristaltic movement of urine towards the bladder. The development of these two cell types coincided with the upregulation of UPIIIa, RAB27b and PPAR gamma reported to be expressed in the urothelium, and several muscle genes, Acta1, Tnnt2, Myocd, and Tpm2. In situ hybridization identified several novel genes with spatial expression within the smooth muscle, Acta1; ureteric mesenchyme and smooth muscle, Thbs2 and Co15a2; and urothelium, Kcnj8 and Adh1. This study marks the first known report defining global gene expression of the developing mouse ureter and will provide insight into the molecular mechanisms underlying kidney and lower urinary tract malformations. (c) 2005 Elsevier B.V. All rights reserved.

Comparative gene expression in brain regions of human alcoholics

The mesocorticolimbic system is the reward centre of the brain and the major target for drugs of abuse including alcohol. Neuroadaptive changes in this region are thought to underlie the process of tolerance and dependence. Recently, several research groups have searched for alcohol-responsive genes using high-throughput microarrays and well-characterized human post-mortem material. Comparison of data from these studies of cortical regions highlights the differences in experimental approach and selection of cases. However, alcohol-responsive gene sets associated with transcription, oxidative stress and energy production were common to these studies. In marked contrast, alcohol-responsive genes in the nucleus accumbens and the ventral tegmental area are primarily associated with changes in neurotransmission and signal transduction. These data support the concept that, within cortical regions, changes in gene expression are associated with alcoholism-related pathology. In the dopaminergic tract of the mesocorticolimbic system, alcohol-responsive gene sets suggest long-term neuroplastic changes in synaptic transmission.

Expression profiling correlates with treatment response in women with advanced serous epithelial ovarian cancer

The majority of epithelial ovarian carcinomas are of serous subtype, with most women presenting at an advanced stage. Approximately 70% respond to initial chemotherapy but eventually relapse. We aimed to find markers of treatment response that might be suitable for routine use, using the gene expression profile of tumor tissue. Thirty one women with histologically-confirmed late-stage serous ovarian cancer were classified into 3 groups based on response to treatment (nonresponders, responders with relapse less than 12 months and responders with no relapse within 12 months). Gene expression profiles of these specimens were analyzed with respect to treatment response and survival (minimum 36 months follow-up). Patients' clinical features did not correlate with prognosis, or with specific gene expression patterns of their tumors. However women who did not respond to treatment could be distinguished from those who responded with no relapse within 12 months based on 34 gene transcripts (p < 0.02). Poor prognosis was associated with high expression of inhibitor of differentiation-2 (ID2) (p = 0.001). High expression of decorin (DCN) and ID2 together was strongly associated with reduced survival (p = 0.003), with an estimated 7-fold increased risk of dying (95% CI 1.9-29.6; 14 months survival) compared with low expression (44 months). Immunohistochemical analysis revealed both nuclear and cytoplasmic distribution of ID2 in ovarian tumors. High percentage of nuclear staining vas associated with poor survival, although not statistically significantly. In conclusion, elevated expression of ID2 and DCN was significantly associated with poor prognosis in a homogeneous group of ovarian cancer patients for whom survival could not be predicted from clinical factors. (c) 2006 Wiley-Liss, Inc.

Expression profiling of purified mouse gonadal somatic cells during the critical time window of sex determination reveals novel candidate genes for human sexual dysgenesis syndromes

Despite the identification of SRY as the testis-determining gene in mammals, the genetic interactions controlling the earliest steps of male sex determination remain poorly understood. In particular, the molecular lesions underlying a high proportion of human XY gonadal dysgenesis, XX maleness and XX true hermaphroditism remain undiscovered. A number of screens have identified candidate genes whose expression is modulated during testis or ovary differentiation in mice, but these screens have used whole gonads, consisting of multiple cell types, or stages of gonadal development well beyond the time of sex determination. We describe here a novel reporter mouse line that expresses enhanced green fluorescent protein under the control of an Sf1 promoter fragment, marking Sertoli and granulosa cell precursors during the critical period of sex determination. These cells were purified from gonads of male and female transgenic embryos at 10.5 dpc (shortly after Sry transcription is activated) and 11.5 dpc (when Sox9 transcription begins), and their transcriptomes analysed using Affymetrix genome arrays. We identified 266 genes, including Dhh, Fgf9 and Ptgds, that were upregulated and 50 genes that were downregulated in 11.5 dpc male somatic gonad cells only, and 242 genes, including Fst, that were upregulated in 11.5 dpc female somatic gonad cells only. The majority of these genes are novel genes that lack identifiable homology, and several human orthologues were found to map to chromosomal loci implicated in disorders of sexual development. These genes represent an important resource with which to piece together the earliest steps of sex determination and gonad development, and provide new candidates for mutation searching in human sexual dysgenesis syndromes.

Spatial gene expression in the T-stage mouse metanephros

The E11.5 mouse metanephros is comprised of a T-stage ureteric epithelial tubule sub-divided into tip and trunk cells surrounded by metanephric mesenchyme (MM). Tip cells are induced to undergo branching morphogenesis by the MM. In contrast, signals within the mesenchyme surrounding the trunk prevent ectopic branching of this region. In order to identify novel genes involved in the molecular regulation of branching morphogenesis we compared the gene expression profiles of isolated tip, trunk and MM cells using Compugen mouse long oligo microarrays. We identified genes enriched in the tip epithelium, sim-1, Arg2, Tacstd1, Crlf-1 and BMP7; genes enriched in the trunk epithelium, Innp1, Itm2b, Mkrn1, SPARC, Emu2 and Gsta3 and genes spatially restricted to the mesenchyme surrounding the trunk, CSPG2 and CV-2, with overlapping and complimentary expression to BMP4, respectively. This study has identified genes spatially expressed in regions of the developing kidney involved in branching morphogenesis, nephrogenesis and the development of the collecting duct system, calyces, renal pelvis and ureter. (c) 2006 Elsevier B.V. All rights reserved.

Gene expression in human alcoholism: Microarray analysis of frontal cortex

Background: Changes in brain gene expression are thought to be responsible for the tolerance, dependence, and neurotoxicity produced by chronic alcohol abuse, but there has been no large scale study of gene expression in human alcoholism. Methods: RNA was extracted from postmortem samples of superior frontal cortex of alcoholics and nonalcoholics. Relative levels of RNA were determined by array techniques. We used both cDNA and oligonucleotide microarrays to provide coverage of a large number of genes and to allow cross-validation for those genes represented on both types of arrays. Results: Expression levels were determined for over 4000 genes and 163 of these were found to differ by 40% or more between alcoholics and nonalcoholics. Analysis of these changes revealed a selective reprogramming of gene expression in this brain region, particularly for myelin-related genes which were downregulated in the alcoholic samples. In addition, cell cycle genes and several neuronal genes were changed in expression. Conclusions: These gene expression changes suggest a mechanism for the loss of cerebral white matter in alcoholics as well as alterations that may lead to the neurotoxic actions of ethanol.

Retrotransposon expression profiling: Unveiling the hidden SINE transcriptome through Next-Generation Sequencing data analysis

Of the ~1.7 million SINE elements in the human genome, only a tiny number are estimated to be active in transcription by RNA polymerase (Pol) III. Tracing the individual loci from which SINE transcripts originate is complicated by their highly repetitive nature. By exploiting RNA-Seq datasets and unique SINE DNA sequences, we devised a bioinformatic pipeline allowing us to identify Pol III-dependent transcripts of individual SINE elements. When applied to ENCODE transcriptomes of seven human cell lines, this search strategy identified ~1300 Alu loci and ~1100 MIR loci corresponding to detectable transcripts, with ~120 and ~60 respectively Alu and MIR loci expressed in at least three cell lines. In vitro transcription of selected SINEs did not reflect their in vivo expression properties, and required the native 5’-flanking region in addition to internal promoter. We also identified a cluster of expressed AluYa5-derived transcription units, juxtaposed to snaR genes on chromosome 19, formed by a promoter-containing left monomer fused to an Alu-unrelated downstream moiety. Autonomous Pol III transcription was also revealed for SINEs nested within Pol II-transcribed genes raising the possibility of an underlying mechanism for Pol II gene regulation by SINE transcriptional units. Moreover the application of our bioinformatic pipeline to both RNA-seq data of cells subjected to an in vitro pro-oncogenic stimulus and of in vivo matched tumor and non-tumor samples allowed us to detect increased Alu RNA expression as well as the source loci of such deregulation. The ability to investigate SINE transcriptomes at single-locus resolution will facilitate both the identification of novel biologically relevant SINE RNAs and the assessment of SINE expression alteration under pathological conditions.

Host gene expression classifiers diagnose acute respiratory illness etiology.

Acute respiratory infections caused by bacterial or viral pathogens are among the most common reasons for seeking medical care. Despite improvements in pathogen-based diagnostics, most patients receive inappropriate antibiotics. Host response biomarkers offer an alternative diagnostic approach to direct antimicrobial use. This observational cohort study determined whether host gene expression patterns discriminate noninfectious from infectious illness and bacterial from viral causes of acute respiratory infection in the acute care setting. Peripheral whole blood gene expression from 273 subjects with community-onset acute respiratory infection (ARI) or noninfectious illness, as well as 44 healthy controls, was measured using microarrays. Sparse logistic regression was used to develop classifiers for bacterial ARI (71 probes), viral ARI (33 probes), or a noninfectious cause of illness (26 probes). Overall accuracy was 87% (238 of 273 concordant with clinical adjudication), which was more accurate than procalcitonin (78%, P < 0.03) and three published classifiers of bacterial versus viral infection (78 to 83%). The classifiers developed here externally validated in five publicly available data sets (AUC, 0.90 to 0.99). A sixth publicly available data set included 25 patients with co-identification of bacterial and viral pathogens. Applying the ARI classifiers defined four distinct groups: a host response to bacterial ARI, viral ARI, coinfection, and neither a bacterial nor a viral response. These findings create an opportunity to develop and use host gene expression classifiers as diagnostic platforms to combat inappropriate antibiotic use and emerging antibiotic resistance.