A segmental nearest neighbor normalization and gene identification method gives superior results for DNA-array analysis

Transcriptomic studies of solventogenic clostridia, Clostridium acetobutylicum and Clostridium beijerinckii

Solventogenic clostridia are not a strictly defined group within the genus Clostridium but its representatives share some common features, i.e. they are anaerobic, non-pathogenic, non-toxinogenic and endospore forming bacteria. Their main metabolite is typically 1-butanol but depending on species and culture conditions, they can form other metabolites such as acetone, isopropanol, ethanol, butyric, lactic and acetic acids, and hydrogen. Although these organisms were previously used for the industrial production of solvents, they later fell into disuse, being replaced by more efficient chemical production. A return to a more biological production of solvents therefore requires a thorough understanding of clostridial metabolism. Transcriptome analysis, which reflects the involvement of individual genes in all cellular processes within a population, at any given (sampling) moment, is a valuable tool for gaining a deeper insight into clostridial life. In this review, we describe techniques to study transcription, summarize the evolution of these techniques and compare methods for data processing and visualization of solventogenic clostridia, particularly the species Clostridium acetobutylicum and Clostridium beijerinckii. Individual approaches for evaluating transcriptomic data are compared and their contributions to advancements in the field are assessed. Moreover, utilization of transcriptomic data for reconstruction of computational clostridial metabolic models is considered and particular models are described. Transcriptional changes in glucose transport, central carbon metabolism, the sporulation cycle, butanol and butyrate stress responses, the influence of lignocellulose-derived inhibitors on growth and solvent production, and other respective topics, are addressed and common trends are highlighted.

Megakaryocytic Maturation in Response to Shear Flow Is Mediated by the Activator Protein 1 (AP-1) Transcription Factor via Mitogen-activated Protein Kinase (MAPK) Mechanotransduction

Megakaryocytes (MKs) are exposed to shear flow as they migrate from the bone marrow hematopoietic compartment into circulation to release pro/preplatelets into circulating blood. Shear forces promote DNA synthesis, polyploidization, and maturation in MKs, and platelet biogenesis. To investigate mechanisms underlying these MK responses to shear, we carried out transcriptional analysis on immature and mature stem cell-derived MKs exposed to physiological shear. In immature (day (d)9) MKs, shear exposure up-regulated genes related to growth and MK maturation, whereas in mature (d12) MKs, it up-regulated genes involved in apoptosis and intracellular transport. Following shear-flow exposure, six activator protein 1 (AP-1) transcripts (ATF4,JUNB,JUN,FOSB,FOS, andJUND) were up-regulated at d9 and two AP-1 proteins (JunD and c-Fos) were up-regulated both at d9 and d12. We show that mitogen-activated protein kinase (MAPK) signaling is linked to both the shear stress response and AP-1 up-regulation. c-Jun N-terminal kinase (JNK) phosphorylation increased significantly following shear stimulation, whereas JNK inhibition reduced shear-induced JunD expression. Although p38 phosphorylation did not increase following shear flow, its inhibition reduced shear-induced JunD and c-Fos expression. JNK inhibition reduced fibrinogen binding and P-selectin expression of d12 platelet-like particles (PLPs), whereas p38 inhibition reduced fibrinogen binding of d12 PLPs. AP-1 expression correlated with increased MK DNA synthesis and polyploidization, which might explain the observed impact of shear on MKs. To summarize, we show that MK exposure to shear forces results in JNK activation, AP-1 up-regulation, and downstream transcriptional changes that promote maturation of immature MKs and platelet biogenesis in mature MKs.

Transcription factors and genetic circuits orchestrating the complex, multilayered response of Clostridium acetobutylicum to butanol and butyrate stress

Background

Organisms of the genus Clostridium are Gram-positive endospore formers of great importance to the carbon cycle, human normo- and pathophysiology, but also in biofuel and biorefinery applications. Exposure of Clostridium organisms to chemical and in particular toxic metabolite stress is ubiquitous in both natural (such as in the human microbiome) and engineered environments, engaging both the general stress response as well as specialized programs. Yet, despite its fundamental and applied significance, it remains largely unexplored at the systems level.

Results

We generated a total of 96 individual sets of microarray data examining the transcriptional changes in C. acetobutylicum, a model Clostridium organism, in response to three levels of chemical stress from the native metabolites, butanol and butyrate. We identified 164 significantly differentially expressed transcriptional regulators and detailed the cellular programs associated with general and stressor-specific responses, many previously unexplored. Pattern-based, comparative genomic analyses enabled us, for the first time, to construct a detailed picture of the genetic circuitry underlying the stress response. Notably, a list of the regulons and DNA binding motifs of the stress-related transcription factors were identified: two heat-shock response regulators, HrcA and CtsR; the SOS response regulator LexA; the redox sensor Rex; and the peroxide sensor PerR. Moreover, several transcriptional regulators controlling stress-responsive amino acid and purine metabolism and their regulons were also identified, including ArgR (arginine biosynthesis and catabolism regulator), HisR (histidine biosynthesis regulator), CymR (cysteine metabolism repressor) and PurR (purine metabolism repressor).

Conclusions

Using an exceptionally large set of temporal transcriptional data and regulon analyses, we successfully built a STRING-based stress response network model integrating important players for the general and specialized metabolite stress response in C. acetobutylicum. Since the majority of the transcription factors and their target genes are highly conserved in other organisms of the Clostridium genus, this network would be largely applicable to other Clostridium organisms. The network informs the molecular basis of Clostridium responses to toxic metabolites in natural ecosystems and the microbiome, and will facilitate the construction of genome-scale models with added regulatory-network dimensions to guide the development of tolerant strains.

Proposed megakaryocytic regulon of p53: the genes engaged to control cell cycle and apoptosis during megakaryocytic differentiation

During endomitosis, megakaryocytes undergo several rounds of DNA synthesis without division leading to polyploidization. In primary megakaryocytes and in the megakaryocytic cell line CHRF, loss or knock-down of p53 enhances cell cycling and inhibits apoptosis, leading to increased polyploidization. To support the hypothesis that p53 suppresses megakaryocytic polyploidization, we show that stable expression of wild-type p53 in K562 cells (a p53-null cell line) attenuates the cells' ability to undergo polyploidization during megakaryocytic differentiation due to diminished DNA synthesis and greater apoptosis. This suggested that p53's effects during megakaryopoiesis are mediated through cell cycle- and apoptosis-related target genes, possibly by arresting DNA synthesis and promoting apoptosis. To identify candidate genes through which p53 mediates these effects, gene expression was compared between p53 knock-down (p53-KD) and control CHRF cells induced to undergo terminal megakaryocytic differentiation using microarray analysis. Among substantially downregulated p53 targets in p53-KD megakaryocytes were cell cycle regulators CDKN1A (p21) and PLK2, proapoptotic FAS, TNFRSF10B, CASP8, NOTCH1, TP53INP1, TP53I3, DRAM1, ZMAT3 and PHLDA3, DNA-damage-related RRM2B and SESN1, and actin component ACTA2, while antiapoptotic CKS1B, BCL2, GTSE1, and p53 family member TP63 were upregulated in p53-KD cells. Additionally, a number of cell cycle-related, proapoptotic, and cytoskeleton-related genes with known functions in megakaryocytes but not known to carry p53-responsive elements were differentially expressed between p53-KD and control CHRF cells. Our data support a model whereby p53 expression during megakaryopoiesis serves to control polyploidization and the transition from endomitosis to apoptosis by impeding cell cycling and promoting apoptosis. Furthermore, we identify a putative p53 regulon that is proposed to orchestrate these effects.

Role of tumor suppressor p53 in megakaryopoiesis and platelet function

The pathobiological role of p53 has been widely studied, however, its role in normophysiology is relatively unexplored. We previously showed that p53 knock-down increased ploidy in megakaryocytic cultures. This study aims to examine the effect of p53 loss on in vivo megakaryopoiesis, platelet production, and function, and to investigate the basis for greater ploidy in p53(-/-) megakaryocytic cultures. Here, we used flow cytometry to analyze ploidy, DNA synthesis, and apoptosis in murine cultured and bone marrow megakaryocytes following thrombopoietin administration and to analyze fibrinogen binding to platelets in vitro. Culture of p53(-/-) marrow cells for 6 days with thrombopoietin gave rise to 1.7-fold more megakaryocytes, 26.1% ± 3.6% of which reached ploidy classes ≥64 N compared to 8.2% ± 0.9% of p53(+/+) megakaryocytes. This was due to 30% greater DNA synthesis in p53(-/-) megakaryocytes and 31% greater apoptosis in p53(+/+) megakaryocytes by day 4 of culture. Although the bone marrow and spleen steady-state megakaryocytic content and ploidy were similar in p53(+/+) and p53(-/-) mice, thrombopoietin administration resulted in increased megakaryocytic polyploidization in p53(-/-) mice. Although their platelet counts were normal, p53(-/-) mice exhibited significantly longer bleeding times and p53(-/-) platelets were less sensitive than p53(+/+) platelets to agonist-induced fibrinogen binding and P-selectin secretion. In summary, our in vivo and ex vivo studies indicate that p53 loss leads to increased polyploidization during megakaryopoiesis. Our findings also suggest for the first time a direct link between p53 loss and the development of fully functional platelets resulting in hemostatic deficiencies.

A genomic-library based discovery of a novel, possibly synthetic, acid-tolerance mechanism in Clostridium acetobutylicum involving non-coding RNAs and ribosomal RNA processing

We generated a genomic library from sheared Clostridium acetobutylicum ATCC 824 DNA, whereby inserts can be expressed in both directions from the thiolase promoter, P(thl). Serial transfer of library-bearing C. acetobutylicum cultures exposed to increasing butyrate concentrations enriched for inserts containing fragments of rRNA genetic loci. The selected library inserts were placed so that antisense (to the rRNAs) non-coding RNAs (ncRNAs) would be transcribed from P(thl). Different enriched inserts imparted similar butyrate-tolerance characteristics. A minimal tolerance fragment (RDNA7) was identified as the 16S-rRNA promoter region. Expressed on plasmid pRD7 off P(thl), RDNA7 can produce putative ncRNAs termed ncRNA(RD7). C. acetobutylicum 824(pRD7) showed superior resistance to butyrate and other carboxylic acids. Transcriptional analysis of butyrate stress identified 120 differentially expressed genes between 824(pRD7) and 824(pSOS95del). The few upregulated genes included the ffh gene of the putative signal recognition particle (SRP) system. Northern analysis of ncRNA(RD7) and corresponding antisense RNAs demonstrated multiple ncRNA(RD7) molecules in 824(pRD7). Several corresponding antisense RNA molecules were identified both in 824(pRD7) and 824(pSOS95del), but at much higher levels in 824(pRD7). Northern analysis of 16S rRNA expression suggested complex RDNA7-dependent rRNA processing. Our data suggest that by hybridizing against unprocessed rRNA precursors, ncRNA(RD7) alters rRNA processing, and these alterations result in acid tolerance, possibly through a mechanism involving the Ffh protein.

Microarray-based gene expression analysis as a process characterization tool to establish comparability of complex biological products: scale-up of a whole-cell immunotherapy product

Whole-cell immunotherapies and other cellular therapies have shown promising results in clinical trials. Due to the complex nature of the whole cell product and of the sometimes limited correlation of clinical potency with the proposed mechanism of action, these cellular immunotherapy products are generally not considered well characterized. Therefore, one major challenge in the product development of whole cell therapies is the ability to demonstrate comparability of product after changes in the manufacturing process. Such changes are nearly inevitable with increase in manufacturing experience leading to improved and robust processes that may have higher commercial feasibility. In order to comprehensively assess the impact of the process changes on the final product, and thus establish comparability, a matrix of characterization assays (in addition to lot release assays) assessing the various aspects of the cellular product are required. In this study, we assessed the capability of DNA-microarray-based, gene-expression analysis as a characterization tool using GVAX cancer immunotherapy cells manufactured by Cell Genesys, Inc. The GVAX immunotherapy product consists two prostate cancer cell lines (CG1940 and CG8711) engineered to secrete human GM-CSF. To demonstrate the capability of the assay, we assessed the transcriptional changes in the product when produced in the presence or absence of fetal bovine serum, and under normal and hypoxic conditions, both changes intended to stress the cell lines. We then assessed the impact of an approximately 10-fold process scale-up on the final product at the transcriptional level. These data were used to develop comparisons and statistical analyses suitable for characterizing culture reproducibility and cellular product similarity. Use of gene-expression data for process characterization proved to be a reproducible and sensitive method for detecting differences due to small or large changes in culture conditions as might be encountered in process scale-up or unanticipated bioprocess failures. Gene expression analysis demonstrated that cell products of representative lots under the same production process and at the same production scale were statistically identical. Large process changes that resulted from the artificial stress conditions used (absence of FBS and induction of hypoxia) displayed profoundly different gene expression patterns. We propose the use of simple t-test analysis in combination with the herein introduced expression ratio with mean intensity (ERMI) analysis as useful tools for process characterization by global gene expression analysis.

Distinct gene subsets in pterygia formation and recurrence: dissecting complex biological phenomenon using genome wide expression data

BACKGROUND

Pterygium is a common ocular surface disease characterized by fibrovascular invasion of the cornea and is sight-threatening due to astigmatism, tear film disturbance, or occlusion of the visual axis. However, the mechanisms for formation and post-surgical recurrence of pterygium are not understood, and a valid animal model does not exist. Here, we investigated the possible mechanisms of pterygium pathogenesis and recurrence.

METHODS

First we performed a genome wide expression analysis (human Affymetrix Genechip, >22000 genes) with principal component analysis and clustering techniques, and validated expression of key molecules with PCR. The controls for this study were the un-involved conjunctival tissue of the same eye obtained during the surgical resection of the lesions. Interesting molecules were further investigated with immunohistochemistry, Western blots, and comparison with tear proteins from pterygium patients.

RESULTS

Principal component analysis in pterygium indicated a signature of matrix-related structural proteins, including fibronectin-1 (both splice-forms), collagen-1A2, keratin-12 and small proline rich protein-1. Immunofluorescence showed strong expression of keratin-6A in all layers, especially the superficial layers, of pterygium epithelium, but absent in the control, with up-regulation and nuclear accumulation of the cell adhesion molecule CD24 in the pterygium epithelium. Western blot shows increased protein expression of beta-microseminoprotein, a protein up-regulated in human cutaneous squamous cell carcinoma. Gene products of 22 up-regulated genes in pterygium have also been found by us in human tears using nano-electrospray-liquid chromatography/mass spectrometry after pterygium surgery. Recurrent disease was associated with up-regulation of sialophorin, a negative regulator of cell adhesion, and never in mitosis a-5, known to be involved in cell motility.

CONCLUSION

Aberrant wound healing is therefore a key process in this disease, and strategies in wound remodeling may be appropriate in halting pterygium or its recurrence. For patients demonstrating a profile of 'recurrence', it may be necessary to manage as a poorer prognostic case and perhaps, more adjunctive treatment after resection of the primary lesion.

Comparative transcriptional analysis of embryoid body versus two-dimensional differentiation of murine embryonic stem cells

Understanding the process of ex vivo embryonic stem (ES) cell differentiation is important for generating higher yields of desirable cell types or lineages and for understanding fundamental aspects of ES differentiation. We used DNA microarray analysis to investigate the differentiation of mouse ES cells cultured under three differentiation conditions. Embryoid body (EB) formation was compared to differentiation on surfaces coated with either gelatin (GEL) or matrigel (MAT). Based on the transcriptional patterns of a list of literature-based "stemness" genes, ES cell differentiation on the two coated surfaces appeared similar but not identical to EB differentiation. A notable difference was the GEL and MAT upregulation but EB downregulation of nine such stemness genes, which are related to cell adhesion and epithelial differentiation. Further, GEL and MAT differentiation showed higher expression of bone formation-related genes (Spp1, Csf1, Gsn, Bmp8b, Crlf1). Gene ontology analysis shows an increase in the expression of genes related to migration and cell structure in all three conditions. Overall, GEL and MAT conditions resulted in a more similar to each other transcriptional profile than to the EB condition, and such differences are apparently related to higher nutrient and metabolite gradients and limitations in the EB versus the GEL or MAT cultures.

The transcriptional program underlying the physiology of clostridial sporulation

A detailed microarray analysis of transcription during sporulation of the strict anaerobe and endospore former Clostridium acetobutylicum is presented.

Background

Clostridia are ancient soil organisms of major importance to human and animal health and physiology, cellulose degradation, and the production of biofuels from renewable resources. Elucidation of their sporulation program is critical for understanding important clostridial programs pertaining to their physiology and their industrial or environmental applications.

Results

Using a sensitive DNA-microarray platform and 25 sampling timepoints, we reveal the genome-scale transcriptional basis of the Clostridium acetobutylicum sporulation program carried deep into stationary phase. A significant fraction of the genes displayed temporal expression in six distinct clusters of expression, which were analyzed with assistance from ontological classifications in order to illuminate all known physiological observations and differentiation stages of this industrial organism. The dynamic orchestration of all known sporulation sigma factors was investigated, whereby in addition to their transcriptional profiles, both in terms of intensity and differential expression, their activity was assessed by the average transcriptional patterns of putative canonical genes of their regulon. All sigma factors of unknown function were investigated by combining transcriptional data with predicted promoter binding motifs and antisense-RNA downregulation to provide a preliminary assessment of their roles in sporulation. Downregulation of two of these sigma factors, CAC1766 and CAP0167, affected the developmental process of sporulation and are apparently novel sporulation-related sigma factors.

Conclusion

This is the first detailed roadmap of clostridial sporulation, the most detailed transcriptional study ever reported for a strict anaerobe and endospore former, and the first reported holistic effort to illuminate cellular physiology and differentiation of a lesser known organism.

The Zur of Xanthomonas campestris functions as a repressor and an activator of putative zinc homeostasis genes via recognizing two distinct sequences within its target promoters

It has been long considered that zinc homeostasis in bacteria is maintained by export systems and uptake systems, which are separately controlled by their own regulators and the uptake systems are negatively regulated by Zur which binds to an about 30-bp AT-rich sequence known as Zur-box present in its target promoters to block the entry of RNA polymerase. Here, we demonstrated in vivo and in vitro that in addition to act as a repressor of putative Zn²⁺-uptake systems, the Zur of the bacterial phytopathogen Xanthomonas campestris pathovar campestris (Xcc) acts as an activator of a Zn²⁺ efflux pump. The Xcc Zur binds to a similar Zur-box with ∼30-bp AT-rich sequence in the promoters of the genes encoding putative Zn²⁺-uptake systems but a 59-bp GC-rich sequence with a 20-bp inverted repeat overlapping the promoter's −35 to −10 sequence of the gene encoding a Zn²⁺-export system. Mutagenesis of the inverted repeat sequence resulted in abolishment of the in vitro binding and the in vivo and in vitro activation of the export gene's promoter by Zur. These results reveal that the Xcc Zur functions as a repressor and an activator of putative zinc homeostasis genes via recognizing two distinct sequences within its target promoters.

Comparative analysis of transcriptional profiling of CD3+, CD4+ and CD8+ T cells identifies novel immune response players in T-cell activation

Background

T-cell activation is an essential step of the immune response and relies on the tightly controlled orchestration of hundreds of genes/proteins, yet the cellular and molecular events underlying this complex process are not fully understood, especially at the genome-scale. Significantly, a comparative genome-scale transcriptional analysis of two T-cell subsets (CD4+ and CD8+) against each other and against the naturally mixed population (CD3+ cells) remains unexplored.

Results

Comparison of the microarray-based gene expression patterns between CD3+ T cells, and the CD4+ and CD8+ subsets revealed largely conserved, but not identical, transcriptional patterns. We employed a Gene-Ontology-driven transcriptional analysis coupled with protein abundance assays in order to identify novel T-cell activation genes and cell-type-specific genes associated with the immune response. We identified potential genes involved in the communication between the two subsets (including IL23A, NR4A2, CD83, PSMB2, -8, MIF, IFI16, TNFAIP1, POU2AF1, and OTUB1) and would-be effector-function-specific genes (XCL2, SLAMF7, TNFSF4, -5, -9, CSF3, CD48 and CD244). Chemokines induced during T-cell activation, but not previously identified in T cells, include CCL20, CXCL9, -10, -11 (in all three populations), and XCL2 (preferentially in CD8+ T cells). Increased expression of other unexpected cytokines (GPI, OSM and MIF) suggests their involvement in T-cell activation with their functions yet to be examined. Differential expression of many receptors, not previously reported in the context of T-cell activation, including CCR5, CCR7, IL1R2, IL1RAP, IL6R, TNFRSF25 and TNFRSF1A, suggests their role in this immune process. Several receptors involved in TCR activation (CD3D, CD3G, TRAT1, ITGAL, ITGB1, ITGB2, CD8A and B (CD8+ T-cell specific) along with LCK, ZAP70 and TYROBP were synchronously downregulated. Members of cell-surface receptors (HLA-Ds and KLRs), none previously identified in the context of T-cell activation, were also downregulated.

Conclusion

This comparative genome-scale, transcriptional analysis of T-cell activation in the CD4+ and CD8+ subsets and the mixed CD3+ populations made possible the identification of many immune-response genes not previously identified in the context of T-cell activation. Significantly, it made possible to identify the temporal patterns of many previously known T-cell activation genes, and also identify genes implicated in effector functions of and communication between CD4+ and CD8+ T cells.

The rsmA-like gene rsmA(Xcc) of Xanthomonas campestris pv. campestris is involved in the control of various cellular processes, including pathogenesis

RsmA is an RNA-binding protein functioning as a global post-transcriptional regulator of various cellular processes in bacteria and has been demonstrated to be an important virulence regulator in many animal bacterial pathogens. However, its function in other phytopathogenic bacteria is unclear, except for the Erwinia carotovora RsmA, which acts as a negative virulence regulator. In this work, we investigated the function of the rsmA-like gene, named rsmA(Xcc), of the phytopathogen Xanthomonas campestris pv. campestris. Deletion of rsmA(Xcc) resulted in complete loss of virulence on the host plant Chinese radish, hypersensitive response on the nonhost plant pepper ECW-10R, and motility on the surface of an agar plate. The rsmA(Xcc) mutant displayed a significant reduction in the production of extracellular amylase, endoglucanase, and polysaccharide, but a significant increase in intracellular glycogen accumulation and an enhanced bacterial aggregation and cell adhesion. Microarray hybridization and semiquantitative reverse-transcription polymerase chain reaction analysis showed that deletion of rsmA(Xcc) led to significantly reduced expression of genes encoding the type III secretion system (T3SS), T3SS-effectors, and the bacterial aggregate dispersing enzyme endo-beta-1,4-mannanase. These results suggest that rsmA(Xcc) is involved in the control of various cellular processes, including pathogenesis of X. campestris pv. campestris.

Genomewide expression analysis in zebrafish mind bomb alleles with pancreas defects of different severity identifies putative Notch responsive genes

Background

Notch signaling is an evolutionarily conserved developmental pathway. Zebrafish mind bomb (mib) mutants carry mutations on mib gene, which encodes a RING E3 ligase required for Notch activation via Delta/Jagged ubiquitylation and internalization.

Methodology/Principal Findings

We examined the mib mutants for defects in pancreas development using in situ hybridization and GFP expression analysis of pancreas-specific GFP lines, carried out the global gene expression profile analysis of three different mib mutant alleles and validated the microarray data using real-time PCR and fluorescent double in situ hybridization. Our study showed that the mib mutants have diminished exocrine pancreas and this defect was most severe in mib^ta52b followed by mib^m132 and then mib^tfi91, which is consistent with the compromised Notch activity found in corresponding mib mutant alleles. Global expression profile analysis of mib mutants showed that there is a significant difference in gene expression profile of wt and three mib mutant alleles. There are 91 differentially expressed genes that are common to all three mib alleles. Through detailed analysis of microarray data, we have identified several previously characterized genes and some putative Notch-responsive genes involved in pancreas development. Moreover, results from real-time PCR and fluorescent double in situ hybridization were largely consistent with microarray data.

Conclusions/Significance

This study provides, for the first time, a global gene expression profile in mib mutants generating useful genomic resources and providing an opportunity to identify the function of novel genes involved in Notch signaling and Notch-regulated developmental processes.

A systems-biology analysis of isogenic megakaryocytic and granulocytic cultures identifies new molecular components of megakaryocytic apoptosis

Background

The differentiation of hematopoietic stem cells into platelet-forming megakaryocytes is of fundamental importance to hemostasis. Constitutive apoptosis is an integral, yet poorly understood, facet of megakaryocytic (Mk) differentiation. Understanding Mk apoptosis could lead to advances in the treatment of Mk and platelet disorders.

Results

We used a Gene-ontology-driven microarray-based transcriptional analysis coupled with protein-level and activity assays to identify genes and pathways involved in Mk apoptosis. Peripheral blood CD34⁺ hematopoietic progenitor cells were induced to either Mk differentiation or, as a negative control without observable apoptosis, granulocytic differentiation. Temporal gene-expression data were analyzed by a combination of intra- and inter-culture comparisons in order to identify Mk-associated genes. This novel approach was first applied to a curated set of general Mk-related genes in order to assess their dynamic transcriptional regulation. When applied to all apoptosis associated genes, it revealed a decrease in NF-κB signaling, which was explored using phosphorylation assays for IκBα and p65 (RELA). Up-regulation was noted among several pro-apoptotic genes not previously associated with Mk apoptosis such as components of the p53 regulon and TNF signaling. Protein-level analyses probed the involvement of the p53-regulated GADD45A, and the apoptosis signal-regulating kinase 1 (ASK1). Down-regulation of anti-apoptotic genes, including several of the Bcl-2 family, was also detected.

Conclusion

Our comparative approach to analyzing dynamic large-scale transcriptional data, which was validated using a known set of Mk genes, robustly identified candidate Mk apoptosis genes. This led to novel insights into the molecular mechanisms regulating apoptosis in Mk cells.

A model-based optimization framework for the inference of regulatory interactions using time-course DNA microarray expression data

Background

Proteins are the primary regulatory agents of transcription even though mRNA expression data alone, from systems like DNA microarrays, are widely used. In addition, the regulation process in genetic systems is inherently non-linear in nature, and most studies employ a time-course analysis of mRNA expression. These considerations should be taken into account in the development of methods for the inference of regulatory interactions in genetic networks.

Results

We use an S-system based model for the transcription and translation process. We propose an optimization-based regulatory network inference approach that uses time-varying data from DNA microarray analysis. Currently, this seems to be the only model-based method that can be used for the analysis of time-course "relative" expressions (expression ratios). We perform an analysis of the dynamic behavior of the system when the number of experimental samples available is varied, when there are different levels of noise in the data and when there are genes that are not considered by the experimenter. Our studies show that the principal factor affecting the ability of a method to infer interactions correctly is the similarity in the time profiles of some or all the genes. The less similar the profiles are to each other the easier it is to infer the interactions. We propose a heuristic method for resolving networks and show that it displays reasonable performance on a synthetic network. Finally, we validate our approach using real experimental data for a chosen subset of genes involved in the sporulation cascade of Bacillus anthracis. We show that the method captures most of the important known interactions between the chosen genes.

Conclusion

The performance of any inference method for regulatory interactions between genes depends on the noise in the data, the existence of unknown genes affecting the network genes, and the similarity in the time profiles of some or all genes. Though subject to these issues, the inference method proposed in this paper would be useful because of its ability to infer important interactions, the fact that it can be used with time-course DNA microarray data and because it is based on a non-linear model of the process that explicitly accounts for the regulatory role of proteins.

Gene expression analysis illuminates the transcriptional programs underlying the functional activity of ex vivo-expanded granulocytes

Global gene expression analysis established the temporal expression patterns and programs underlying the development of functional activity of ex vivo-expanded (EXE) human granulocytes, as well as differences compared with peripheral blood (PB) granulocytes. CD34(+) progenitor cells were cultured for 3 wk to induce rapid expansion and granulocytic differentiation, with 40% CD15(+) cells by day 3 and 90% by day 12. Phagocytic and respiratory burst activity increased with the fraction of CD15(++)CD11b(+) cells (myelocytes to segmented) and peaked by day 17. However, only 25% of CD15(++)CD11b(+) cells were phagocytic, and respiratory burst activity was one-third that of PB granulocytes. EXE granulocytes from later days and PB granulocytes showed similar expression of Fc gamma receptors (-1A, -2A, -2C, -3A) and complement receptors (-1, -3, -4). Later downregulation of CD36 (expressed by macrophages) suggests lineage plasticity early in granulocytic differentiation. Expression in mature EXE and PB granulocytes was similar for most Fc gamma receptor-mediated phagocytosis signaling proteins, including high-level expression of Hck, Fgr, and the actin-related protein 2/3 complex. Lower expression of Lyn, Cdc42, pleckstrin, and PKC beta(I) by EXE granulocytes may explain decreased phagocytosis. PB and mature EXE granulocytes expressed similar levels of NADPH oxidase complex genes and receptors for fMLP-mediated respiratory burst. Lower burst activity by EXE granulocytes may result from lower expression of Raf1 and PKC zeta. Elevated expression of toll-like receptor (TLR)2, TLR1, and CD14 in mature EXE and PB granulocytes supports a role for the TLR2 and CD14 pathway in zymosan-mediated respiratory burst activity. Lower activity in EXE granulocytes may be due to greater expression of IRAK3, which inhibits TLR-mediated signaling.

Deregulated CDC25A expression promotes mammary tumorigenesis with genomic instability

Checkpoint pathways help cells maintain genomic integrity, delaying cell cycle progression in response to various risks of fidelity, such as genotoxic stresses, compromised DNA replication, and impaired spindle control. Cancer cells frequently exhibit genomic instability, and recent studies showed that checkpoint pathways are likely to serve as a tumor-suppressive barrier in vivo. The cell cycle-promoting phosphatase CDC25A is an activator of cyclin-dependent kinases and one of the downstream targets for the CHK1-mediated checkpoint pathway. Whereas CDC25A overexpression is observed in various human cancer tissues, it has not been determined whether deregulated CDC25A expression triggers or promotes tumorigenesis in vivo. Here, we show that transgenic expression of CDC25A cooperates markedly with oncogenic ras or neu in murine mammary tumorigenesis. MMTV-CDC25A transgenic mice exhibit alveolar hyperplasia in the mammary tissue but do not develop spontaneous mammary tumors. The MMTV-CDC25A transgene markedly shortens latency of tumorigenesis in MMTV-ras mice. The MMTV-CDC25A transgene also accelerates tumor growth in MMTV-neu mice with apparent cell cycle miscoordination. CDC25A-overexpressing tumors, which invade more aggressively, exhibit various chromosomal aberrations on fragile regions, including the mouse counterpart of human 1p31-36, according to array-based comparative genomic hybridization and karyotyping. The chromosomal aberrations account for substantial changes in gene expression profile rendered by transgenic expression of CDC25A, including down-regulation of Trp73. These data indicate that deregulated control of cellular CDC25A levels leads to in vivo genomic instability, which cooperates with the neu-ras oncogenic pathway in mammary tumorigenesis.

Global transcriptional analysis delineates the differential inflammatory response interleukin-15 elicits from cultured human T cells

OBJECTIVE

Interleukin 15 (IL)-15 controls proliferation and survival of T cells, but its effects and the underlying cellular regulation are not well understood. Previous studies have focused on its effects on short-term T-cell cultures. In view of the potential problems associated with using IL-2 alone in adoptive immunotherapy protocols, we investigated the impact of IL-15 on T-cell cultures and the global transcriptional effects it elicits in such cultures.

MATERIALS AND METHODS

DNA microarrays and flow cytometry were used to examine the differential effect of 20 ng/mL IL-15 on primary serum-free T-cell cultures activated and cultured in the presence of IL-2. Quantitative reverse transcriptase polymerase chain reaction confirmed select microarray data.

RESULTS

IL-15 significantly increased ex vivo expansion of primary human T cells over the entire 11-day expansions without affecting viability. The 1133 genes were consistently differentially expressed among three donor samples. Ontological analysis demonstrated that IL-15 increases expression of genes involved in inflammatory response (e.g., tumor necrosis factor [TNF]-alpha, Oncostatin M, CD40L, and CD33) and apoptosis (e.g., TNF-related apoptosis-inducing ligand). IL-15 also induced expression of four suppressors of cytokine signaling (SOCS) family genes (SOCS1-3, cytokine-inducible SH2-containing protein), which are classical negative regulators of cytokine signaling. IL-15 strongly suppressed the expression of inhibitory natural killer cell receptor genes, including three C-type lectins (KLRB1, KLRC1, and KLRD1), as well as IL-7Ra and Granzyme H. Finally, IL-15 induced differential expression of TNF receptor superfamily members (CD27 and CD30).

CONCLUSION

These findings suggest that exogenous IL-15 may have a potential role in adoptive immunotherapy by both enhancing proliferation and modulating functionality during ex vivo T-cell expansion.

Functional genomics via multiscale analysis: application to gene expression and ChIP-on-chip data

We present a fast, versatile and adaptive-multiscale algorithm for analyzing a wide-variety of DNA microarray data. Its primary application is in normalization of array data as well as subsequent identification of 'enriched targets', e.g. differentially expressed genes in expression profiling arrays and enriched sites in ChIP-on-chip experimental data. We show how to accommodate the unique characteristics of ChIP-on-chip data, where the set of 'enriched targets' is large, asymmetric and whose proportion to the whole data varies locally.

SUPPLEMENTARY INFORMATION

Supplementary figures, related preprint, free software as well as our raw DNA microarray data with PCR validations are available at http://www.math.umn.edu/~lerman/supp/bioinfo06 as well as Bioinformatics online.

Genomewide expression profiling in the zebrafish embryo identifies target genes regulated by Hedgehog signaling during vertebrate development

Hedgehog proteins play critical roles in organizing the embryonic development of animals, largely through modulation of target gene expression. Little is currently known, however, about the kinds and numbers of genes whose expression is controlled, directly or indirectly, by Hedgehog activity. Using techniques to globally repress or activate Hedgehog signaling in zebrafish embryos followed by microarray-based expression profiling, we have discovered a cohort of genes whose expression responds significantly to loss or gain of Hedgehog function. We have confirmed the Hedgehog responsiveness of a representative set of these genes with whole-mount in situ hybridization as well as real time PCR. In addition, we show that the consensus Gli-binding motif is enriched within the putative regulatory elements of a sizeable proportion of genes that showed positive regulation in our assay, indicating that their expression is directly induced by Hedgehog. Finally, we provide evidence that the Hedgehog-dependent spatially restricted transcription of one such gene, nkx2.9, is indeed mediated by Gli1 through a single Gli recognition site located within an evolutionarily conserved enhancer fragment. Taken together, this study represents the first comprehensive survey of target genes regulated by the Hedgehog pathway during vertebrate development. Our data also demonstrate for the first time the functionality of the Gli-binding motif in the control of Hedgehog signaling-induced gene expression in the zebrafish embryo.

Genome scale analysis of diffusible signal factor regulon in Xanthomonas campestris pv. campestris: identification of novel cell-cell communication-dependent genes and functions

The bacterial pathogen Xanthomonas campestris pv. campestris (Xcc) recruits a diffusible signal factor (DSF), which has recently been structurally characterized as cis-11-methyl-2-dodecenoic acid, as a cell-cell communication signal to synchronize virulence gene expression and biofilm dispersal. In this study, we showed that despite the existance of phenotype variations in different Xcc isolates, the DSF-mediated functions were in general conserved. To investigate the genomic profiles of DSF regulation, we designed and conducted oligomicroarray analysis by comparison of the gene expression patterns of wild-type strain XC1 and its DSF-deficient mutant XC1dF, as well as those of XC1dF in the presence or absence of DSF signals. The analyses led to identification of 165 genes, whose expression was significantly influenced by DSF signals. These genes encode proteins and enzymes belonging to at least 12 functional groups. In addition to those previously known DSF-dependent activities such as production of extracellular enzymes and extracellular polysaccharides, microarray analyses also revealed new functions mediated by DSF, such as flagellum synthesis, resistance to toxins and oxidative stress, and aerobic respiration. Phenotype analyses confirmed that DSF signalling contributed to resistance to toxin acriflavin and hydrogen peroxide, and to the survival of bacterial cells at different temperatures. We conclude that DSF cell-cell signalling is not only essential for co-ordinating the expression of virulence genes but also plays a vital role in keeping up the general competence of the pathogen in ecosystems.

A novel normalization method for effective removal of systematic variation in microarray data

Normalization of cDNA and oligonucleotide microarray data has become a standard procedure to offset non-biological differences between two samples for accurate identification of differentially expressed genes. Although there are many normalization techniques available, their ability to accurately remove systematic variation has not been sufficiently evaluated. In this study, we performed experimental validation of various normalization methods in order to assess their ability to accurately offset non-biological differences (systematic variation). The limitations of many existing normalization methods become apparent when there are unbalanced shifts in transcript levels. To overcome this limitation, we have proposed a novel normalization method that uses a matching algorithm for the distribution peaks of the expression log ratio. The robustness and effectiveness of this method was evaluated using both experimental and simulated data.

Expression of insulin-like growth factor binding protein-3 in pterygium tissue

BACKGROUND

Pterygium is a fibrovascular overgrowth from the conjunctiva onto the cornea. The pathogenesis of this common ocular surface disorder is not well understood and the only treatment is surgical removal.

METHODS

DNA microarray analysis of primary pterygium tissue was carried out using uninvolved conjunctiva tissue as a comparison for gene expression levels. Real time polymerase chain reaction (PCR) was used to verify the mRNA level of expression for genes changed in pterygium. Western blot analysis and immunohistochemistry showed protein expression levels and the tissue distribution.

RESULTS

Microarray analysis revealed that mRNA levels of a number of genes were changed in primary pterygium. In particular the gene for insulin-like growth factor binding protein-3, (IGFBP3), which modulates the effects of insulin-like growth factor on cells was significantly decreased. Both the message and protein expression of IGFBP3 in pterygium were decreased compared to normal, uninvolved conjunctiva.

CONCLUSION

Decreased levels of IGFBP3 protein have been strongly correlated with the presence of cancer. Identification of the low level of expression of IGFBP3 in pterygium suggests that the pathway controlling cell proliferation has lost an important control mechanism, which may explain the continued growth of pterygium.

Breast cancer prognosis by combinatorial analysis of gene expression data

INTRODUCTION

The potential of applying data analysis tools to microarray data for diagnosis and prognosis is illustrated on the recent breast cancer dataset of van 't Veer and coworkers. We re-examine that dataset using the novel technique of logical analysis of data (LAD), with the double objective of discovering patterns characteristic for cases with good or poor outcome, using them for accurate and justifiable predictions; and deriving novel information about the role of genes, the existence of special classes of cases, and other factors.

METHOD

Data were analyzed using the combinatorics and optimization-based method of LAD, recently shown to provide highly accurate diagnostic and prognostic systems in cardiology, cancer proteomics, hematology, pulmonology, and other disciplines.

RESULTS

LAD identified a subset of 17 of the 25,000 genes, capable of fully distinguishing between patients with poor, respectively good prognoses. An extensive list of 'patterns' or 'combinatorial biomarkers' (that is, combinations of genes and limitations on their expression levels) was generated, and 40 patterns were used to create a prognostic system, shown to have 100% and 92.9% weighted accuracy on the training and test sets, respectively. The prognostic system uses fewer genes than other methods, and has similar or better accuracy than those reported in other studies. Out of the 17 genes identified by LAD, three (respectively, five) were shown to play a significant role in determining poor (respectively, good) prognosis. Two new classes of patients (described by similar sets of covering patterns, gene expression ranges, and clinical features) were discovered. As a by-product of the study, it is shown that the training and the test sets of van 't Veer have differing characteristics.

CONCLUSION

The study shows that LAD provides an accurate and fully explanatory prognostic system for breast cancer using genomic data (that is, a system that, in addition to predicting good or poor prognosis, provides an individualized explanation of the reasons for that prognosis for each patient). Moreover, the LAD model provides valuable insights into the roles of individual and combinatorial biomarkers, allows the discovery of new classes of patients, and generates a vast library of biomedical research hypotheses.

Coordinated patterns of cytochrome bd and lactate dehydrogenase expression in Bacillus subtilis

A variety of pathways for electron and carbon flow in the soil bacterium Bacillus subtilis are differentially expressed depending on whether oxygen is present in the cell environment. This study characterizes the regulation of the respiratory oxidase cytochrome bd and the NADH-linked fermentative lactate dehydrogenase (LDH). Transcription of the cydABCD operon, encoding cytochrome bd, is highly regulated and only becomes activated at low oxygen availability. This induction is not dependent on the gene encoding the redox regulator Fnr or the genes encoding the ResDE two-component regulatory system. The DNA-binding protein YdiH was found to be a principal regulator that controls cydABCD expression. Transcription from the cyd promoter is stimulated 15-fold by a region located upstream of the core promoter. The upstream region may constitute a binding site for an unidentified transcription activator that is likely to influence the level of transcription but not its timing, which is negatively controlled by YdiH. This report provides evidence that YdiH also functions as a repressor of the ldh gene encoding LDH and of a gene, ywcJ, which encodes a putative formate-nitrite transporter. Based on the similarity between YdiH and the Rex protein of Streptomyces coelicolor, it is proposed that YdiH serves as a redox sensor, the activity of which is regulated by cellular differences in the free levels of NAD+ and NADH. It is suggested that ydiH be renamed as rex.

Transcriptional program of early sporulation and stationary-phase events in Clostridium acetobutylicum

DNA microarray analysis of Clostridium acetobutylicum was used to examine the genomic-scale gene expression changes during the shift from exponential-phase growth and acidogenesis to stationary phase and solventogenesis. Self-organizing maps were used to identify novel expression patterns of functional gene classes, including aromatic and branched-chain amino acid synthesis, ribosomal proteins, cobalt and iron transporters, cobalamin biosynthesis, and lipid biosynthesis. The majority of pSOL1 megaplasmid genes (in addition to the solventogenic genes aad-ctfA-ctfB and adc) had increased expression at the onset of solventogenesis, suggesting that other megaplasmid genes may play a role in stationary-phase phenomena. Analysis of sporulation genes and comparison with published Bacillus subtilis results indicated conserved expression patterns of early sporulation genes, including spo0A, the sigF operon, and putative canonical genes of the sigma(H) and sigma(F) regulons. However, sigE expression could not be detected within 7.5 h of initial spo0A expression, consistent with the observed extended time between the appearance of clostridial forms and endospore formation. The results were compared with microarray comparisons of the wild-type strain and the nonsolventogenic, asporogenous M5 strain, which lacks the pSOL1 megaplasmid. While some results were similar, the expression of primary metabolism genes and heat shock proteins was higher in M5, suggesting a difference in metabolic regulation or a butyrate stress response in M5. The results of this microarray platform and analysis were further validated by comparing gene expression patterns to previously published Northern analyses, reporter assays, and two-dimensional protein electrophoresis data of metabolic genes (including all major solventogenesis genes), sporulation genes, heat shock proteins, and other solventogenesis-induced gene expression.

Transcriptome analysis of zebrafish embryogenesis using microarrays

Zebrafish (Danio rerio) is a well-recognized model for the study of vertebrate developmental genetics, yet at the same time little is known about the transcriptional events that underlie zebrafish embryogenesis. Here we have employed microarray analysis to study the temporal activity of developmentally regulated genes during zebrafish embryogenesis. Transcriptome analysis at 12 different embryonic time points covering five different developmental stages (maternal, blastula, gastrula, segmentation, and pharyngula) revealed a highly dynamic transcriptional profile. Hierarchical clustering, stage-specific clustering, and algorithms to detect onset and peak of gene expression revealed clearly demarcated transcript clusters with maximum gene activity at distinct developmental stages as well as co-regulated expression of gene groups involved in dedicated functions such as organogenesis. Our study also revealed a previously unidentified cohort of genes that are transcribed prior to the mid-blastula transition, a time point earlier than when the zygotic genome was traditionally thought to become active. Here we provide, for the first time to our knowledge, a comprehensive list of developmentally regulated zebrafish genes and their expression profiles during embryogenesis, including novel information on the temporal expression of several thousand previously uncharacterized genes. The expression data generated from this study are accessible to all interested scientists from our institute resource database (http://giscompute.gis.a-star.edu.sg/~govind/zebrafish/data_download.html).

Intracellular butyryl phosphate and acetyl phosphate concentrations in Clostridium acetobutylicum and their implications for solvent formation

It has been suggested (L. H. Harris, R. P. Desai, N. E. Welker, and E. T. Papoutsakis, Biotechnol. Bioeng. 67:1-11, 2000) that butyryl phosphate (BuP) is a regulator of solventogenesis in Clostridium acetobutylicum. Here, we determined BuP and acetyl phosphate (AcP) levels in fermentations of C. acetobutylicum wild type (WT), degenerate strain M5, a butyrate kinase (buk) mutant, and a phosphotransacetylase (pta) mutant. A sensitive method was developed to measure BuP and AcP in the same sample. Compared to the WT, the buk mutant had higher levels of BuP and AcP; the BuP levels were high during the early exponential phase, and there was a peak corresponding to solvent production. Consistent with this, solvent formation was initiated significantly earlier and was much stronger in the buk mutant than in all other strains. For all strains, initiation of butanol formation corresponded to a BuP peak concentration that was more than 60 to 70 pmol/g (dry weight), and higher and sustained levels corresponded to higher butanol formation fluxes. The BuP levels never exceeded 40 to 50 pmol/g (dry weight) in strain M5, which produces no solvents. The BuP profiles were bimodal, and there was a second peak midway through solventogenesis that corresponded to carboxylic acid reutilization. AcP showed a delayed single peak during late solventogenesis corresponding to acetate reutilization. As expected, in the pta mutant the AcP levels were very low, yet this strain exhibited strong butanol production. These data suggest that BuP is a regulatory molecule that may act as a phosphodonor of transcriptional factors. DNA array-based transcriptional analysis of the buk and M5 mutants demonstrated that high BuP levels corresponded to downregulation of flagellar genes and upregulation of solvent formation and stress genes.

Molecular understanding of oxygen-tension and patient-variability effects on ex vivo expanded T cells

Immunotherapy with ex vivo cultured T cells depends on a large supply of biologically active cells. Understanding the effects of culture parameters is essential for improving the proliferation and efficacy of the expanded cells. Low oxygen tension (5% pO(2)) was previously reported to improve T-cell expansion and alter cellular phenotypic characteristics compared to T cells cultured at 20% pO(2). Here we report the use of DNA-array based transcriptional analysis coupled with protein-level analysis to provide molecular insights into pO(2) and patient-variability effects on expanded primary human T cells. Analysis of seven blood samples showed that reduced pO(2) results in higher expression of genes important in lymphocyte biology, immune function, and cell-cycle progression. 20% pO(2) resulted in higher expression of genes involved in stress response, cell death, and cellular repair. Expression of granzyme A (gzmA) was found to be significantly regulated by oxygen tension with cells at 5% pO(2) having greater gzmA expression than at 20% pO(2). Protein-level analysis of gzmA was consistent with transcriptional analysis. Granzyme K (gzmK) was coexpressed with gzmA, whereas Granzyme B (gzmB) expression was found to precede the expression of both gzmA and gzmK in 15-day cultures. Temporal gene expression patterns for seven blood samples demonstrate that most genes are expressed by all patient samples in similar temporal patterns. However, several patient-specific gene clusters were identified, and one cluster was found to correlate well with cell proliferation and may potentially be used to predict patient-specific T-cell expansion.

Molecular insights into the pleiotropic effects of plasma on ex vivo-expanded T cells using DNA-microarray analysis

OBJECTIVE

Immunotherapy with ex vivo-expanded T cells depends on a large supply of biologically active cells. Understanding the effects of culture parameters is essential for improving cell expansion and efficacy. We used DNA-microarray and flow-cytometric analysis coupled with functional assays to investigate mechanistic aspects of plasma supplementation in ex vivo T-cell expansion.

METHODS

The effect of plasma supplementation on 18 primary T-cell cultures over a 15-day expansion was investigated. Transcriptional analysis of 5 samples was done with time points every 2 to 3 days throughout the 15-day expansion. Quantitative RT-PCR analysis was used to confirm selected microarray data. The expression of granzyme A and vimentin were analyzed using intracellular flow cytometry. T-cell functionality was assessed using a mixed leukocyte reaction (MLR).

RESULTS

We show that the increased expansion of plasma-supplemented cultures of primary human T cells is mostly due to increased cell survival. T cells from plasma-supplemented cultures show higher expression of immunoglobulin genes, integrins, and genes of cytotoxic granules, suggesting a possible enhanced immune function. This was confirmed using a mixed leukocyte reaction and intracellular granzyme-A measurements. A distinct gene expression pattern was correlated to viability differences between plasma-supplemented and serum-free cultures. Ontological analysis of genes in this pattern suggests that the decreased viability of serum-free cultures correlates with higher expression of actin-cytoskeleton and lipid-metabolism genes. Vimentin was found to be expressed higher in serum-free cultures.

CONCLUSIONS

These results indicate that the observed decreased cytotoxicity of T cells cultured in serum-free media may be due to increased oxidative stress and cytoskeleton degradation.

Rank Difference Analysis of Microarrays (RDAM), a novel approach to statistical analysis of microarray expression profiling data

Background

A key step in the analysis of microarray expression profiling data is the identification of genes that display statistically significant changes in expression signals between two biological conditions.

Results

We describe a new method, Rank Difference Analysis of Microarrays (RDAM), which estimates the total number of truly varying genes and assigns a p-value to each signal variation. Information on a group of differentially expressed genes includes the sensitivity and the false discovery rate. We demonstrate the feasibility and efficiency of our approach by applying it to a large synthetic expression data set and to a biological data set obtained by comparing vegetatively-growing wild type and tor2-mutant yeast strains. In both cases we observed a significant improvement of the power of analysis when our method is compared to another popular nonparametric method.

Conclusions

This study provided a valuable new statistical method to analyze microarray data. We conclude that the good quality of the results obtained by RDAM is mainly due to the quasi-perfect equalization of variation distribution, which is related to the standardization procedure used and to the measurement of variation by rank difference.

Transcriptional analysis of product-concentration driven changes in cellular programs of recombinant Clostridium acetobutylicumstrains

Antisense RNA (asRNA) downregulation alters protein expression without changing the regulation of gene expression. Downregulation of primary metabolic enzymes possibly combined with overexpression of other metabolic enzymes may result in profound changes in product formation, and this may alter the large-scale transcriptional program of the cells. DNA-array based large-scale transcriptional analysis has the potential to elucidate factors that control cellular fluxes even in the absence of proteome data. These themes are explored in the study of large-scale transcriptional analysis programs and the in vivo primary-metabolism fluxes of several related recombinant C. acetobutylicum strains: C. acetobutylicum ATCC 824(pSOS95del) (plasmid control; produces high levels of butanol snd acetone), 824(pCTFB1AS) (expresses antisense RNA against CoA transferase (ctfb1-asRNA); produces very low levels of butanol and acetone), and 824(pAADB1) (expresses ctfb1-asRNA and the alcohol-aldehyde dahydrogenase gene (aad); produce high alcohol and low acetone levels). DNA-array based transcriptional analysis revealed that the large changes in product concentrations (snd notably butanol concentration) due to ctfb1-asRNA expression alone and in combination with aad overexpression resulted in dramatic changes of the cellular transcriptome. Cluster analysis and gene expression patterns of established and putative operons involved in stress response, motility, sporulation, and fatty-acid biosynthesis indicate that these simple genetic changes dramatically alter the cellular programs of C. acetobutylicum. Comparison of gene expression and flux analysis data may point to possible flux-controling steps and suggest unknown regulatory mechanisms.

Transcriptional analysis of spo0A overexpression in Clostridium acetobutylicum and its effect on the cell's response to butanol stress

Spo0A is the regulator of stationary-phase events and is required for transcription of solvent formation genes in Clostridium acetobutylicum. In order to elucidate the role of spo0A in differentiation, we performed transcriptional analysis of 824(pMSPOA) (a spo0A-overexpressing C. acetobutylicum strain with enhanced sporulation) against a plasmid control strain. DNA microarray data were contrasted to data from a spo0A knockout strain (SKO1) that neither sporulates nor produces solvents. Transcripts of fatty acid metabolism genes, motility and chemotaxis genes, heat shock protein genes, and genes encoding the Fts family of cell division proteins were differentially expressed in the two strains, suggesting that these genes play roles in sporulation and the solvent stress response. 824(pMSPOA) alone showed significant downregulation of many glycolytic genes in stationary phase, which is consistent with metabolic flux analysis data. Surprisingly, spo0A overexpression resulted in only nominal transcriptional changes of regulatory genes (abrB and sigF) whose expression was significantly altered in SKO1. Overexpression of spo0A imparted increased tolerance and prolonged metabolism in response to butanol stress. While most of the differentially expressed genes appear to be part of a general stress response (similar to patterns in two plasmid control strains and a groESL-overexpressing strain), several genes were expressed at higher levels at early time points after butanol challenge only in 824(pMSPOA). Most of these genes were related to butyryl coenzyme A and butyrate formation and/or assimilation, but they also included the cell division gene ftsX, the gyrase subunit-encoding genes gyrB and gyrA, DNA synthesis and repair genes, and fatty acid synthesis genes, all of which might play a role in the immediate butanol stress response, and thus in enhanced butanol tolerance.

Transcriptional analysis of butanol stress and tolerance in Clostridium acetobutylicum

The effects of challenges with low (0.25%, vol/vol) and high (0.75%) concentrations of butanol on the growth, glucose metabolism, product formation, and transcriptional program of the solvent-tolerant Clostridium acetobutylicum strain 824(pGROE1) and the plasmid control strain 824(pSOS95del) were used to study solvent tolerance and stress response. Strain 824(pGROE1) was generated by groESL overexpression. The growth of 824(pGROE1) was less inhibited than that of 824(pSOS95del), and 824(pGROE1) was able to metabolize glucose over the entire course of the culture (60 h postchallenge) while glucose metabolism in 824(pSOS95del) lasted 24 h. A comparison of their respective DNA array-based transcriptional profiles identified genes with similar expression patterns (these genes are likely to be part of a general butanol stress response) and genes with opposite expression patterns (these genes are likely to be associated with increased tolerance to butanol). Both strains exhibited a butanol dose-dependent increase in expression of all major stress protein genes, including groES, dnaKJ, hsp18, and hsp90; all major solvent formation genes, including aad, ctfA and -B, adc, and bdhA and -B (an unexpected and counterintuitive finding); the butyrate formation genes (ptb and buk); the butyryl coenzyme A biosynthesis operon genes; fructose bisphosphate aldolase; and a gene with homology to Bacillus subtilis kinA. A dose-dependent decrease in expression was observed for the genes of the major fatty acid synthesis operon (also an unexpected and counterintuitive finding), several glycolytic genes, and a few sporulation genes. Genes with opposite expression kinetics included rlpA, artP, and a gene encoding a hemin permease. Taken together, these data suggest that stress, even when it derives from the solvent product itself, triggers the induction of the solvent formation genes.

Process for the integrated extraction, identification and quantification of metabolites, proteins and RNA to reveal their co-regulation in biochemical networks

A novel extraction protocol is described with which metabolites, proteins and RNA are sequentially extracted from the same sample, thereby providing a convenient procedure for the analysis of replicates as well as exploiting the inherent biological variation of independent samples for multivariate data analysis. A detection of 652 metabolites, 297 proteins and clear RNA bands in a single Arabidopsis thaliana leaf sample was validated by analysis with gas chromatography coupled to a time of flight mass spectrometer for metabolites, two-dimensional liquid chromatography coupled to mass spectrometry for proteins, and Northern blot analysis for RNA. A subset of the most abundant proteins and metabolites from replicate analysis of different Arabidopsis accessions was merged to form an integrative dataset allowing both classification of different genotypes and the unbiased analysis of the hierarchical organization of proteins and metabolites within a real biochemical network.

Overexpression of groESL in Clostridium acetobutylicum results in increased solvent production and tolerance, prolonged metabolism, and changes in the cell's transcriptional program

DNA array and Western analyses were used to examine the effects of groESL overexpression and host-plasmid interactions on solvent production in Clostridium acetobutylicum ATCC 824. Strain 824(pGROE1) was created to overexpress the groESL operon genes from a clostridial thiolase promoter. The growth of 824(pGROE1) was inhibited up to 85% less by a butanol challenge than that of the control strain, 824(pSOS95del). Overexpression of groESL resulted in increased final solvent titers 40% and 33% higher than those of the wild type and plasmid control strains, respectively. Active metabolism lasted two and one half times longer in 824(pGROE1) than in the wild type. Transcriptional analysis of 824(pGROE1) revealed increased expression of motility and chemotaxis genes and a decrease in the expression of the other major stress response genes. Decreased expression of the dnaKJ operon upon overexpression of groESL suggests that groESL functions as a modulator of the CIRCE regulon, which is shown here to include the hsp90 gene. Analysis of the plasmid control strain 824(pSOS95del) revealed complex host-plasmid interactions relative to the wild-type strain, resulting in prolonged biphasic growth and metabolism. Decreased expression of four DNA gyrases resulted in differential expression of many key primary metabolism genes. The ftsA and ftsZ genes were expressed at higher levels in 824(pSOS95del), revealing an altered cell division and sporulation pattern. Both transcriptional and Western analyses revealed elevated stress protein expression in the plasmid-carrying strain.

Spurious spatial periodicity of co-expression in microarray data due to printing design

Global transcriptome data is increasingly combined with sophisticated mathematical analyses to extract information about the functional state of a cell. Yet the extent to which the results reflect experimental bias at the expense of true biological information remains largely unknown. Here we show that the spatial arrangement of probes on microarrays and the particulars of the printing procedure significantly affect the log-ratio data of mRNA expression levels measured during the Saccharomyces cerevisiae cell cycle. We present a numerical method that filters out these technology-derived contributions from the existing transcriptome data, leading to improved functional predictions. The example presented here underlines the need to routinely search and compensate for inherent experimental bias when analyzing systematically collected, internally consistent biological data sets.

DNA array-based transcriptional analysis of asporogenous, nonsolventogenic Clostridium acetobutylicum strains SKO1 and M5

The large-scale transcriptional program of two Clostridium acetobutylicum strains (SKO1 and M5) relative to that of the parent strain (wild type [WT]) was examined by using DNA microarrays. Glass DNA arrays containing a selected set of 1,019 genes (including all 178 pSOL1 genes) covering more than 25% of the whole genome were designed, constructed, and validated for data reliability. Strain SKO1, with an inactivated spo0A gene, displays an asporogenous, filamentous, and largely deficient solventogenic phenotype. SKO1 displays downregulation of all solvent formation genes, sigF, and carbohydrate metabolism genes (similar to genes expressed as part of the stationary-phase response in Bacillus subtilis) but also several electron transport genes. A major cluster of genes upregulated in SKO1 includes abrB, the genes from the major chemotaxis and motility operons, and glycosylation genes. Strain M5 displays an asporogenous and nonsolventogenic phenotype due to loss of the megaplasmid pSOL1, which contains all genes necessary for solvent formation. Therefore, M5 displays downregulation of all pSOL1 genes expressed in the WT. Notable among other genes expressed more highly in WT than in M5 were sigF, several two-component histidine kinases, spo0A, cheA, cheC, many stress response genes, fts family genes, DNA topoisomerase genes, and central-carbon metabolism genes. Genes expressed more highly in M5 include electron transport genes (but different from those downregulated in SKO1) and several motility and chemotaxis genes. Most of these expression patterns were consistent with phenotypic characteristics. Several of these expression patterns are new or different from what is known in B. subtilis and can be used to test a number of functional-genomic hypotheses.

Antisense RNA downregulation of coenzyme A transferase combined with alcohol-aldehyde dehydrogenase overexpression leads to predominantly alcohologenic Clostridium acetobutylicum fermentations

Plasmid pAADB1 for the overexpression of the alcohol-aldehyde dehydrogenase (aad) gene and downregulation of the coenzyme A transferase (CoAT) using antisense RNA (asRNA) against ctfB (the second CoAT gene on the polycistronic aad-ctfA-ctfB message) was used in order to increase the butanol/acetone ratio of Clostridium acetobutylicum ATCC 824 fermentations. Acetone and butanol levels were drastically reduced in 824(pCTFB1AS) (expresses only an asRNA against ctfB) compared to 824(pSOS95del) (plasmid control). Compared to strain 824(pCTFB1AS), 824(pAADB1) fermentations exhibited two profound differences. First, butanol levels were ca. 2.8-fold higher in 824(pAADB1) and restored back to plasmid control levels, thus supporting the hypothesis that asRNA downregulation of ctfB leads to degradation of the whole aad-ctfA-ctfB transcript. Second, ethanol titers in 824(pAADB1) were ca. 23-fold higher and the highest (ca. 200 mM) ever reported in C. acetobutylicum. Western blot analysis confirmed that CoAT was downregulated in 824(pAADB1) at nearly the same levels as in strain 824(pCTFB1AS). Butyrate depletion in 824(pAADB1) fermentations suggested that butyryl-CoA was limiting butanol production in 824(pAADB1). This was confirmed by exogenously adding butyric acid to 824(pAADB1) fermentations to increase the butanol/ethanol ratio. DNA microarray analysis showed that aad overexpression profoundly affects the large-scale transcriptional program of the cells. Several classes of genes were differentially expressed [strain 824(pAADB1) versus strain 824(pCTFB1AS)], including genes of the stress response, sporulation, and chemotaxis. The expression patterns of the CoAT genes (ctfA and ctfB) and aad were consistent with the overexpression of aad and asRNA downregulation of ctfB.