cDNA targets improve whole blood gene expression profiling and enhance detection of pharmocodynamic biomarkers: a quantitative platform analysis.

BACKGROUND: Genome-wide gene expression profiling of whole blood is an attractive method for discovery of biomarkers due to its non-invasiveness, simple clinical site processing and rich biological content. Except for a few successes, this technology has not yet matured enough to reach its full potential of identifying biomarkers useful for clinical prognostic and diagnostic applications or in monitoring patient response to therapeutic intervention. A variety of technical problems have hampered efforts to utilize this technology for identification of biomarkers. One significant hurdle has been the high and variable concentrations of globin transcripts in whole blood total RNA potentially resulting in non-specific probe binding and high background. In this study, we investigated and quantified the power of three whole blood profiling approaches to detect meaningful biological expression patterns. METHODS: To compare and quantify the impact of different mitigation technologies, we used a globin transcript spike-in strategy to synthetically generate a globin-induced signature and then mitigate it with the three different technologies. Biological differences, in globin transcript spiked samples, were modeled by supplementing with either 1% of liver or 1% brain total RNA. In order to demonstrate the biological utility of a robust globin artifact mitigation strategy in biomarker discovery, we treated whole blood ex vivo with suberoylanilide hydroxamic acid (SAHA) and compared the overlap between the obtained signatures and signatures of a known biomarker derived from SAHA-treated cell lines and PBMCs of SAHA-treated patients. RESULTS: We found cDNA hybridization targets detect at least 20 times more specific differentially expressed signatures (2597) between 1% liver and 1% brain in globin-supplemented samples than the PNA (117) or no treatment (97) method at FDR = 10% and p-value < 3x10-3. In addition, we found that the ex vivo derived gene expression profile was highly concordant with that of the previously identified SAHA pharmacodynamic biomarkers. CONCLUSIONS: We conclude that an amplification method for gene expression profiling employing cDNA targets effectively mitigates the negative impact on data of abundant globin transcripts and greatly improves the ability to identify relevant gene expression based pharmacodynamic biomarkers from whole blood.

Regulation of gene expression by lithium and depletion of inositol in slices of adult rat cortex.

Lithium inhibits inositol monophosphatase at therapeutically effective concentrations, and it has been hypothesized that depletion of brain inositol levels is an important chemical alteration for lithium's therapeutic efficacy in bipolar disorder. We have employed adult rat cortical slices as a model to investigate the gene regulatory consequences of inositol depletion effected by lithium using cytidine diphosphoryl-diacylglycerol as a functionally relevant biochemical marker to define treatment conditions. Genes coding for the neuropeptide hormone pituitary adenylate cyclase activating polypeptide (PACAP) and the enzyme that processes PACAP's precursor to the mature form, peptidylglycine alpha-amidating monooxygenase, were upregulated by inositol depletion. Previous work has shown that PACAP can increase tyrosine hydroxylase (TH) activity and dopamine release, and we found that the gene for GTP cyclohydrolase, which effectively regulates TH through synthesis of tetrahydrobiopterin, was also upregulated by inositol depletion. We propose that modulation of brain PACAP signaling might represent a new opportunity in the treatment of bipolar disorder.

Mutations that affect vacuole biogenesis inhibit proliferation of the endoplasmic reticulum in Saccharomyces cerevisiae.

In yeast, increased levels of the sterol biosynthetic enzyme, 3-hydroxy-3-methylglutaryl coenzyme A (HMG-CoA) reductase isozyme, Hmg1p, induce assembly of nuclear-associated ER membranes called karmellae. To identify additional genes involved in karmellae assembly, we screened temperature-sensitive mutants for karmellae assembly defects. Two independently isolated, temperature-sensitive strains that were also defective for karmellae biogenesis carried mutations in VPS16, a gene involved in vacuolar protein sorting. Karmellae biogenesis was defective in all 13 other vacuole biogenesis mutants tested, although the severity of the karmellae assembly defect varied depending on the particular mutation. The hypersensitivity of 14 vacuole biogenesis mutants to tunicamycin was well correlated with pronounced defects in karmellae assembly, suggesting that the karmellae assembly defect reflected alteration of ER structure or function. Consistent with this hypothesis, seven of eight mutations causing defects in secretion also affected karmellae assembly. However, the vacuole biogenesis mutants were able to proliferate their ER in response to Hmg2p, indicating that the mutants did not have a global defect in the process of ER biogenesis.

A microarray platform comparison for neuroscience applications.

To address the need for high sensitivity in gene expression profiling of small neural tissue samples ( approximately 100 ng total RNA), we compared a novel RT-PCR-IVT protocol using fluor-reverse pairs on inkjet oligonucleotide microarrays and an RT-IVT protocol using 33P labeling on nylon cDNA arrays. The comparison protocol was designed to evaluate these systems for sensitivity, specificity, reproducibility, and linearity. We developed parameters, thresholds, and testing conditions that could be used to differentiate various systems that spanned detection chemistry and instrumentation; probe number and selection criteria; and sample processing protocols. We concluded that the inkjet system had better performance in sensitivity, specificity, and reproducibility than the nylon system, and similar performance in linearity. Between these two platforms, the data indicates that the inkjet system would perform better for the transcriptional profiling of 100 ng total RNA samples for neuroscience studies.

Multi-tissue coexpression networks reveal unexpected subnetworks associated with disease.

BACKGROUND: Obesity is a particularly complex disease that at least partially involves genetic and environmental perturbations to gene-networks connecting the hypothalamus and several metabolic tissues, resulting in an energy imbalance at the systems level. RESULTS: To provide an inter-tissue view of obesity with respect to molecular states that are associated with physiological states, we developed a framework for constructing tissue-to-tissue coexpression networks between genes in the hypothalamus, liver or adipose tissue. These networks have a scale-free architecture and are strikingly independent of gene-gene coexpression networks that are constructed from more standard analyses of single tissues. This is the first systematic effort to study inter-tissue relationships and highlights genes in the hypothalamus that act as information relays in the control of peripheral tissues in obese mice. The subnetworks identified as specific to tissue-to-tissue interactions are enriched in genes that have obesity-relevant biological functions such as circadian rhythm, energy balance, stress response, or immune response. CONCLUSIONS: Tissue-to-tissue networks enable the identification of disease-specific genes that respond to changes induced by different tissues and they also provide unique details regarding candidate genes for obesity that are identified in genome-wide association studies. Identifying such genes from single tissue analyses would be difficult or impossible.

Endoplasmic reticulum-associated degradation is required for cold adaptation and regulation of sterol biosynthesis in the yeast Saccharomyces cerevisiae.

Endoplasmic reticulum-associated degradation (ERAD) mediates the turnover of short-lived and misfolded proteins in the ER membrane or lumen. In spite of its important role, only subtle growth phenotypes have been associated with defects in ERAD. We have discovered that the ERAD proteins Ubc7 (Qri8), Cue1, and Doa10 (Ssm4) are required for growth of yeast that express high levels of the sterol biosynthetic enzyme, 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR). Interestingly, the observed growth defect was exacerbated at low temperatures, producing an HMGR-dependent cold sensitivity. Yeast strains lacking UBC7, CUE1, or DOA10 also assembled aberrant karmellae (ordered arrays of membranes surrounding the nucleus that assemble when HMGR is expressed at high levels). However, rather than reflecting the accumulation of abnormal karmellae, the cold sensitivity of these ERAD mutants was due to increased HMGR catalytic activity. Mutations that compromise proteasomal function also resulted in cold-sensitive growth of yeast with elevated HMGR, suggesting that improper degradation of ERAD targets might be responsible for the observed cold-sensitive phenotype. However, the essential ERAD targets were not the yeast HMGR enzymes themselves. The sterol metabolite profile of ubc7Delta cells was altered relative to that of wild-type cells. Since sterol levels are known to regulate membrane fluidity, the viability of ERAD mutants expressing normal levels of HMGR was examined at low temperatures. Cells lacking UBC7, CUE1, or DOA10 were cold sensitive, suggesting that these ERAD proteins have a role in cold adaptation, perhaps through effects on sterol biosynthesis.

Proliferation of the endoplasmic reticulum occurs normally in cells that lack a functional unfolded protein response.

Increased expression of certain ER membrane proteins leads to biogenesis of novel ER membrane arrays. These structures provide models in which to explore the mechanisms by which cells control the size and organization of organelles in response to changing physiological demands. In yeast, elevated levels of HMG-CoA reductase induce ER arrays known as karmellae. Cox and co-workers (1997) discovered that karmellae assembly is toxic to ire1 mutants. These mutants are unable to initiate the unfolded protein response, which enables cells to adjust levels of ER chaperones in response to stresses. We sought to determine whether the karmellae-dependent death of ire1 mutants was due to karmellae assembly or to increased levels of HMG-CoA reductase activity. Unexpectedly, we found that ire1 cells could assemble normal levels of karmellae that were structurally identical to those of wild-type cells. In addition, karmellae assembly did not itself induce the unfolded protein response. Certain ire1 strains produced significant numbers of transformants that were unable to utilize galactose as sole carbon source. These results suggest that the karmellae-dependent death of certain ire1 strains may simply reflect their inability to grow on galactose.

Parallel analysis of tagged deletion mutants efficiently identifies genes involved in endoplasmic reticulum biogenesis.

Increased levels of HMG-CoA reductase induce cell type- and isozyme-specific proliferation of the endoplasmic reticulum. In yeast, the ER proliferations induced by Hmg1p consist of nuclear-associated stacks of smooth ER membranes known as karmellae. To identify genes required for karmellae assembly, we compared the composition of populations of homozygous diploid S. cerevisiae deletion mutants following 20 generations of growth with and without karmellae. Using an initial population of 1,557 deletion mutants, 120 potential mutants were identified as a result of three independent experiments. Each experiment produced a largely non-overlapping set of potential mutants, suggesting that differences in specific growth conditions could be used to maximize the comprehensiveness of similar parallel analysis screens. Only two genes, UBC7 and YAL011W, were identified in all three experiments. Subsequent analysis of individual mutant strains confirmed that each experiment was identifying valid mutations, based on the mutant's sensitivity to elevated HMG-CoA reductase and inability to assemble normal karmellae. The largest class of HMG-CoA reductase-sensitive mutations was a subset of genes that are involved in chromatin structure and transcriptional regulation, suggesting that karmellae assembly requires changes in transcription or that the presence of karmellae may interfere with normal transcriptional regulation.

Optimization of oligonucleotide arrays and RNA amplification protocols for analysis of transcript structure and alternative splicing.

Microarrays offer a high-resolution means for monitoring pre-mRNA splicing on a genomic scale. We have developed a novel, unbiased amplification protocol that permits labeling of entire transcripts. Also, hybridization conditions, probe characteristics, and analysis algorithms were optimized for detection of exons, exon-intron edges, and exon junctions. These optimized protocols can be used to detect small variations and isoform mixtures, map the tissue specificity of known human alternative isoforms, and provide a robust, scalable platform for high-throughput discovery of alternative splicing.