Improving DNA array data quality by minimising 'neighbourhood' effects

Dose-dependent regulation of target gene expression and cell proliferation by c-Myc levels

The proto-oncogene c-myc encodes a basic helix-loop-helix leucine zipper transcription factor (c-Myc). c-Myc plays a crucial role in cell growth and proliferation. Here, we examined how expression of c-Myc target genes and cell proliferation depend on variation of c-Myc protein levels. We show that proliferation rates, the number of cells in S-phase, and cell size increased in a dose-dependent manner in response to increasing c-Myc levels. Likewise, the mRNA levels of c-Myc responsive genes steadily increased with rising c-Myc levels. Strikingly, steady-state mRNA levels of c-Myc target genes did not saturate even at highest c-Myc concentrations. These characteristics predestine c-Myc levels as a cellular rheostat for the control and fine-tuning of cell proliferation and growth rates.

c-Myc and Rel/NF-kappaB are the two master transcriptional systems activated in the latency III program of Epstein-Barr virus-immortalized B cells

The Epstein-Barr virus (EBV) latency III program imposed by EBNA2 and LMP1 is directly responsible for immortalization of B cells in vitro and is thought to mediate most immunodeficiency-related posttransplant lymphoproliferative diseases in vivo. To answer the question whether and how this proliferation program is related to c-Myc, we have established the transcriptome of both c-Myc and EBV latency III proliferation programs using a Lymphochip specialized microarray. In addition to EBV-positive latency I Burkitt lymphoma lines and lymphoblastoid cell lines (LCLs), we used an LCL expressing an estrogen-regulatable EBNA2 fusion protein (EREB2-5) and derivative B-cell lines expressing a constitutively active or tetracycline-regulatable c-myc gene. A total of 897 genes were found to be fourfold or more up- or downregulated in either one or both proliferation programs compared to the expression profile of resting EREB2-5 cells. A total of 661 (74%) of these were regulated similarly in both programs. Numerous repressed genes were known targets of STAT1, and most induced genes were known to be upregulated by c-Myc and to be involved in cell proliferation. In keeping with the gene expression patterns, inactivation of c-Myc by a chemical inhibitor or by conditional expression of dominant-negative c-Myc and Max mutants led to proliferation arrest of LCLs. Most genes differently regulated in both proliferation programs corresponded to genes induced by NF-kappaB in LCLs, and many of them coded for immunoregulatory and/or antiapoptotic molecules. Thus, c-Myc and NF-kappaB are the two main transcription factors responsible for the phenotype, growth pattern, and biological properties of cells driven into proliferation by EBV.

Single-cell transcript profiling of barley attacked by the powdery mildew fungus

In many plant-pathogen interactions, there are several possible outcomes for simultaneous attacks on the same leaf. For instance, an attack by the powdery mildew fungus on one barley leaf epidermal cell may succeed in infection and formation of a functional haustorium, whereas a neighboring cell attacked at the same time may resist fungal penetration. To date, the mixed cellular responses seen even in susceptible host leaves have made it difficult to relate induced changes in gene expression to resistance or susceptibility in bulk leaf samples. By microextraction of cell-specific mRNA and subsequent cDNA array analysis, we have successfully obtained separate gene expression profiles for specific mildew-resistant and -infected barley cells. Thus, for the first time, it is possible to identify genes that are specifically regulated in infected cells and, presumably, involved in fungal establishment. Further, although much is understood about the genetic basis of effective papilla resistance associated with mutant mlo barley, we provide here the first evidence for gene regulation associated with effective papilla-based nonspecific resistance expressed in nominally "susceptible" wild-type barley.

The majority of Escherichia coli mRNAs undergo post-transcriptional modification in exponentially growing cells

Polyadenylation of RNAs by poly(A) polymerase I (PAP I) in Escherichia coli plays a significant role in mRNA decay and general RNA quality control. However, many important features of this system, including the prevalence of polyadenylated mRNAs in the bacterium, are still poorly understood. By comparing the transcriptomes of wild-type and pcnB deletion strains using macroarray analysis, we demonstrate that >90% of E.coli open reading frames (ORFs) transcribed during exponential growth undergo some degree of polyadenylation by PAP I, either as full-length transcripts or decay intermediates. Detailed analysis of over 240 transcripts suggests that Rho-independent transcription terminators serve as polyadenylation signals. Conversely, mRNAs terminated in a Rho-dependent fashion are probably not substrates for PAP I, but can be modified by the addition of long polynucleotide tails through the biosynthetic activity of polynucleotide phosphorylase (PNPase). Furthermore, real-time PCR analysis indicates that the extent of polyadenylation of individual full-length transcripts such as lpp and ompA varies significantly in wild-type cells. The data presented here demonstrates that polyadenylation in E.coli occurs much more frequently than previously envisioned.

Development of a coral cDNA array to examine gene expression profiles in Montastraea faveolata exposed to environmental stress

The development of a cDNA array of coral genes and its application to investigate changes in coral gene expression associated with stressful conditions is described. The array includes both well-characterized and previously unidentified coral genes from Acropora cervicornis and Montastraea faveolata. Corals were exposed to either natural or anthropogenic stressors to elicit the expression of stress genes for isolation and incorporation onto the array. A total of 32 genes involved in protein synthesis, apoptosis, cell signaling, metabolism, cellular defense and inflammation were included on the array. Labeled cDNA from coral (Montastraea faveolata) exposed to elevated seawater temperature, salinity and ultraviolet light was tested against the microarray to determine patterns of gene expression associated with each stressor. Carbonic anhydrase, thioredoxin, a urokinase plasminogen activator receptor (uPAR) and three ribosomal genes demonstrated differential expression across all replicates on the array and between replicate colonies. Specific gene expression patterns produced in response to different stressors demonstrate the potential for gene expression profiling in characterizing the coral stress response.

Panel of Bacillus subtilis reporter strains indicative of various modes of action

In a recent project, we collected the transcriptional profiles of Bacillus subtilis 168 after treatment with a large set of diverse antibacterial agents. One result of the data analysis was the identification of marker genes that are indicative of certain compounds or compound classes. We cloned these promoter regions in front of a luciferase reporter gene and reintroduced the constructs individually into the B. subtilis chromosome. Strains were analyzed for their responsiveness after treatment with a set of 37 antibacterials. Twelve functional reporter strains were generated that were selectively and significantly upregulated by the compounds. The selectivity of the reporter strains ranged from generic pathways like protein biosynthesis, cell wall biosynthesis, and fatty acid biosynthesis to compound classes (quinolones and glycopeptides) and individual compounds (rifampin, cycloserine, and clindamycin). Five of the strains are amenable for high-throughput applications, e.g., pathway-specific screening. In summary, we successfully generated B. subtilis reporter strains that are indicative of the mechanisms of action of various classes of antibacterials. The set of reporter strains presented herein can be used for mode-of-action analyses and for whole-cell screening of compound libraries in a mode-of-action-specific manner.

Prediction of mechanisms of action of antibacterial compounds by gene expression profiling

We have generated a database of expression profiles carrying the transcriptional responses of the model organism Bacillus subtilis following treatment with 37 well-characterized antibacterial compounds of different classes. The database was used to build a predictor for the assignment of the mechanisms of action (MoAs) of antibacterial compounds by the use of support vector machines. This predictor was able to correctly classify the MoA class for most compounds tested. Furthermore, we provide evidence that the in vivo MoA of hexachlorophene does not match the MoA predicted from in vitro data, a situation frequently faced in drug discovery. A database of this kind may facilitate the prioritization of novel antibacterial entities in drug discovery programs. Potential applications and limitations are discussed.

Validation of a novel, fully integrated and flexible microarray benchtop facility for gene expression profiling

Here we describe a novel microarray platform that integrates all functions needed to perform any array-based experiment in a compact instrument on the researcher's laboratory benchtop. Oligonucle otide probes are synthesized in situ via a light- activated process within the channels of a three-dimensional microfluidic reaction carrier. Arrays can be designed and produced within hours according to the user's requirements. They are processed in a fully automatic workflow. We have characterized this new platform with regard to dynamic range, discrimination power, reproducibility and accuracy of biological results. The instrument detects sample RNAs present at a frequency of 1:100 000. Detection is quantitative over more than two orders of magnitude. Experiments on four identical arrays with 6398 features each revealed a mean coefficient of variation (CV) value of 0.09 for the 6398 unprocessed raw intensities indicating high reproducibility. In a more elaborate experiment targeting 1125 yeast genes from an unbiased selection, a mean CV of 0.11 on the fold change level was found. Analyzing the transcriptional response of yeast to osmotic shock, we found that biological data acquired on our platform are in good agreement with data from Affymetrix GeneChips, quantitative real-time PCR and--albeit somewhat less clearly--to data from spotted cDNA arrays obtained from the literature.

Feature extraction and signal processing for nylon DNA microarrays

Background

High-density DNA microarrays require automatic feature extraction methodologies and softwares. These can be a potential source of non-reproducibility of gene expression measurements. Variation in feature location or in signal integration methodology may be a significant contribution to the observed variance in gene expression levels.

Results

We explore sources of variability in feature extraction from DNA microarrays on Nylon membrane with radioactive detection. We introduce a mathematical model of the signal emission and derive methods for correcting biases such as overshining, saturation or variation in probe amount. We also provide a quality metric which can be used qualitatively to flag weak or untrusted signals or quantitatively to modulate the weight of each experiment or gene in higher level analyses (clustering or discriminant analysis).

Conclusions

Our novel feature extraction methodology, based on a mathematical model of the radioactive emission, reduces variability due to saturation, neighbourhood effects and variable probe amount. Furthermore, we provide a fully automatic feature extraction software, BZScan, which implements the algorithms described in this paper.

Genomic analysis in Escherichia coli demonstrates differential roles for polynucleotide phosphorylase and RNase II in mRNA abundance and decay

Previous work has shown that simultaneous inactivation of polynucleotide phosphorylase (PNPase) and RNase II (both 3' 5' exonucleases) in Escherichia coli leads to the loss of cell viability and the accumulation of partially degraded mRNA species. In order to help to distinguish how these two enzymes globally affect the abundance and decay of mRNAs, we have carried out a genome-wide analysis of the steady-state levels of E. coli transcripts using deletion mutations in either rnb or pnp. The data show that, in exponentially growing cells, inactivation of PNPase leads to an increase in the steady-state level of more expressed mRNAs (17.3%) than inactivation of RNase II (7.3%). In contrast, the steady-state levels of a large number of E. coli mRNAs (31%) are decreased in the absence of RNase II, including almost all the ribosomal protein genes, suggesting that a major function of this enzyme is to protect specific mRNAs from the activity of other ribonucleases. Array data were confirmed by Northern analysis of 12 individual mRNAs. A comparison between the steady-state levels and the half-lives of individual mRNAs indicates that there may be a direct interaction between transcription and mRNA decay for some of the transcripts. In addition, results are presented to show significant phenotypic differences between the pnp-7 point mutant and the pnp delta 683 deletion allele.