Numerical methods for handling uncertainty in microarray data: an example analyzing perturbed mitochondrial function in yeast

Multipurpose instantaneous microarray detection of acute encephalitis causing viruses and their expression profiles

Detection of multiple viruses is important for global analysis of gene or protein content and expression, opening up new prospects in terms of molecular and physiological systems for pathogenic diagnosis. Early diagnosis is crucial for disease treatment and control as it reduces inappropriate use of antiviral therapy and focuses surveillance activity. This requires the ability to detect and accurately diagnose infection at or close to the source/outbreak with minimum delay and the need for specific, accessible point-of-care diagnosis able to distinguish causative viruses and their subtypes. None of the available viral diagnostic assays combine a point-of-care format with the complex capability to identify a large range of human and animal viruses. Microarray detection provides a useful, labor-saving tool for detection of multiple viruses with several advantages, such as convenience and prevention of cross-contamination of polymerase chain reaction (PCR) products, which is of foremost importance in such applications. Recently, real-time PCR assays with the ability to confirm the amplification product and quantitate the target concentration have been developed. Furthermore, nucleotide sequence analysis of amplification products has facilitated epidemiological studies of infectious disease outbreaks and monitoring of treatment outcomes for infections, in particular for viruses that mutate at high frequency. This review discusses applications of microarray technology as a potential new tool for detection and identification of acute encephalitis-causing viruses in human serum, plasma, and cell cultures.

Evidence for a novel mitochondria-to-nucleus signalling pathway in respiring cells lacking i-AAA protease and the ABC-transporter Mdl1

Peptides generated upon degradation of mitochondrial proteins by various ATP-dependent proteases are continuously released from mitochondria raising the intriguing possibility of a role of these peptides in interorganellar communication. Here, we have determined genome-wide transcript profiles of mutant yeast cells defective in mitochondrial peptide export. Deletion of YME1, coding for the i-AAA protease in the inner membrane, abolished peptide generation in the intermembrane space and led to the induction of nuclear genes with functions in mitochondrial gene expression and the biogenesis of the respiratory chain. On the other hand, deletion of MDL1, coding for an ABC-transporter involved in peptide export from the matrix space, only had minor effects on nuclear gene expression. It strengthened, however, the response in Deltayme1 cells suggesting a link between mitochondrial peptide export and nuclear gene expression. The response in Yme1-deficient cells depended on respiratory growth and was not observed in fermenting yeast cells. Inhibition of the F1FO-ATP synthase induced Deltayme1 responsive genes whereas inhibition of the respiratory chain or dissipation of the mitochondrial membrane potential resulted in their repression. These findings suggest the existence of a novel mitochondria-to-nucleus signalling pathway in respiring cells which allows the re-adjustment of the biogenesis of the respiratory chain in response to an altered activity of the F1FO-ATP synthase.

A genomic and proteomic study of the spectrum of nonalcoholic fatty liver disease

Nonalcoholic fatty liver disease (NAFLD) is the most common chronic liver disease, and some of its forms are progressive. This study describes the profiling of hepatic gene expression and serum protein content in patients with different subtypes of NAFLD. Liver biopsy specimens from 98 bariatric surgery patients were classified as normal, steatosis alone, steatosis with nonspecific inflammation, and nonalcoholic steatohepatitis (NASH). Microarray hybridizations were performed in triplicate and the microarray expression levels of a selected group of genes were confirmed using real-time quantitative reverse-transcriptase polymerase chain reaction. Serum protein profiles of the same patients were determined by SELDI-TOF mass spectrometry. Of 98 obese patients, 91 were diagnosed with NAFLD (12 steatosis alone, 52 steatosis with nonspecific inflammation, and 27 NASH), and 7 patients without NAFLD served as obese controls. Each group of NAFLD patients was compared with the obese controls, and 22 genes with more than twofold differences in expression levels were revealed. Proteomics analyses were performed for the same group comparisons and revealed twelve significantly different protein peaks. In conclusion, this genomic/proteomic analysis suggests differential expression of several genes and protein peaks in patients within and across the forms of NAFLD. These findings may help clarify the pathogenesis of NAFLD and identify potential targets for therapeutic intervention.

Hepatic gene expression in patients with obesity-related non-alcoholic steatohepatitis

BACKGROUND

Non-alcoholic fatty liver disease (NAFLD) is among the most common causes of chronic liver disease. NAFLD includes a spectrum of clinicopathologic syndromes that includes non-alcoholic steatohepatitis (NASH) that has potential for progression. The pathogenesis of NASH is poorly characterized.

AIM

This study was designed to identify differences in hepatic gene expression in patients with NASH and to relate such differences to their clinical characteristics.

DESIGN

Consecutive patients undergoing bariatric surgery were prospectively recruited. Extensive clinical data and two liver biopsy specimens were obtained at the time of enrollment. A single hepatopathologist reviewed and classified the liver biopsies. Patients with excessive alcohol use and other causes of liver disease were excluded. A group of 29 NASH patients, 12 with steatosis alone, seven obese controls and six non-obese controls were selected for further investigation. Customized cDNA microarrays containing 5220 relevant genes were designed specifically for this study. Microarray experiments were run in triplicate for each sample and a selected group of genes were confirmed using real-time PCR.

OUTCOME MEASURE

Differential hepatic gene expressions in patients with NASH as compared with controls.

RESULTS

Thirty-four genes with significant differential expression were identified in patients with NASH when compared with non-obese controls. Moreover, 19 of these genes showed no significant expression differences in obese vs. non-obese controls, suggesting a stronger association of these genes to NASH.

CONCLUSIONS

Several differentially expressed genes in patients with NASH are related to lipid metabolism and extracellular matrix remodeling. Additionally, genes related to liver regeneration, apoptosis, and the detoxification process were differentially expressed. These findings may help clarify the molecular pathogenesis of NASH and identify potential targets for therapeutic intervention.

Development of a new reference standard for microarray experiments

Often microarray studies require a reference to indirectly compare the samples under observation. References based on pooled RNA from different cell lines have already been described (here referred to as RNA-R), but they usually do not exhaustively represent the set of genes printed on a chip, thus requiring many adjustments during the analyses. A reference could also be generated in vitro transcribing the collection of cDNA clones printed on the microarray in use (here referred to as T3-R). Here we describe an alternative and simpler PCR-based methodology to construct a similar reference (Chip-R), and we extensively test and compare it to both RNA-R and T3-R. The use of both Chip-R and T3-R dramatically increases the number of signals on the slides and gives more reproducible results than RNA-R. Each reference preparation is also evaluated in a simple microarray experiment comparing two different RNA populations. Our results show that the introduction of a reference always interferes with the analysis. Indeed, the direct comparison is able to identify more up- or down-regulated genes than any reference-mediated analysis. However, if a reference has to be used, Chip-R and T3-R are able to guarantee more reliable results than RNA-R.

Comparison of hybridization behavior between double and single strands of targets and the application of asymmetric PCR targets in cDNA microarray

Double stranded targets on the cDNA microarray contain representatives of both the coding and noncoding strands, which will introduce hybridization competition with probes. Here, the effect of double and single strands of targets on the signal intensity and the ratios of Cy5/Cy3 within the same slide were compared. The results show that single stranded targets can increase the hybridization efficiency without changing the Cy5/Cy3 ratio. Based on these results, a new strategy was established by generating cDNA targets with asymmetric PCR, instead of conventional PCR, to increase the sensitivity of the cDNA microarray. Furthermore, the feasibility of this approach was validated. The results indicate that the cDNA microarray system based on asymmetric PCR is more sensitive, with no decrease in the reliability and reproducibility as compared with that based on conventional symmetric PCR.

Alterations of O-glycosylation, cell wall, and mitochondrial metabolism in Kluyveromyces lactis cells defective in KlPmr1p, the Golgi Ca2+-ATPase

In yeast the P-type Ca(2+)-ATPase of the Golgi apparatus, Pmr1p, is the most important player in calcium homeostasis. In Kluyveromyces lactis KlPMR1 inactivation leads to pleiotropic phenotypes, including reduced N-glycosylation and altered cell wall morphogenesis. To study the physiology of K. lactis when KlPMR1 was inactivated microarrays containing all Saccharomyces cerevisiae coding sequences were utilized. Alterations in O-glycosylation, consistent with the repression of KlPMT2, were found and a terminal N-acetylglucosamine in the O-glycans was identified. Klpmr1Delta cells showed increased expression of PIRs, proteins involved in cell wall maintenance, suggesting that responses to cell wall weakening take place in K. lactis. We found over-expression of KlPDA1 and KlACS2 genes involved in the Acetyl-CoA synthesis and down-regulation of KlIDP1, KlACO1, and KlSDH2 genes involved in respiratory metabolism. Increases in oxygen consumption and succinate dehydrogenase activity were also observed in mutant cells. The described approach highlighted the unexpected involvement of KlPMR1 in energy-yielding processes.

Genome-wide identification of DNaseI hypersensitive sites using active chromatin sequence libraries

Comprehensive identification of sequences that regulate transcription is one of the major goals of genome biology. Focal alteration in chromatin structure in vivo, detectable through hypersensitivity to DNaseI and other nucleases, is the sine qua non of a diverse cast of transcriptional regulatory elements including enhancers, promoters, insulators, and locus control regions. We developed an approach for genome-scale identification of DNaseI hypersensitive sites (HSs) via isolation and cloning of in vivo DNaseI cleavage sites to create libraries of active chromatin sequences (ACSs). Here, we describe analysis of >61,000 ACSs derived from erythroid cells. We observed peaks in the density of ACSs at the transcriptional start sites of known genes at non-gene-associated CpG islands, and, to a lesser degree, at evolutionarily conserved noncoding sequences. Peaks in ACS density paralleled the distribution of DNaseI HSs. ACSs and DNaseI HSs were distributed between both expressed and nonexpressed genes, suggesting that a large proportion of genes reside within open chromatin domains. The results permit a quantitative approximation of the distribution of HSs and classical cis-regulatory sequences in the human genome.

Determining a detectable threshold of signal intensity in cDNA microarray based on accumulated distribution

In microarray data mining, one of the key problems is how to handle weak signals. Based on a bent piecewise linear accumulated distribution generally found in the microarray data, a new detectable threshold finding method is proposed to filter genes with unreliable information in this paper. More reliable and reproducible data is produced for the subsequent data mining.

Listeria monocytogenes sigma B regulates stress response and virulence functions

While the stress-responsive alternative sigma factor sigma(B) has been identified in different species of Bacillus, Listeria, and Staphylococcus, the sigma(B) regulon has been extensively characterized only in B. subtilis. We combined biocomputing and microarray-based strategies to identify sigma(B)-dependent genes in the facultative intracellular pathogen Listeria monocytogenes. Hidden Markov model (HMM)-based searches identified 170 candidate sigma(B)-dependent promoter sequences in the strain EGD-e genome sequence. These data were used to develop a specialized, 208-gene microarray, which included 166 genes downstream of HMM-predicted sigma(B)-dependent promoters as well as selected virulence and stress response genes. RNA for the microarray experiments was isolated from both wild-type and Delta sigB null mutant L. monocytogenes cells grown to stationary phase or exposed to osmotic stress (0.5 M KCl). Microarray analyses identified a total of 55 genes with statistically significant sigma(B)-dependent expression under the conditions used in these experiments, with at least 1.5-fold-higher expression in the wild type over the sigB mutant under either stress condition (51 genes showed at least 2.0-fold-higher expression in the wild type). Of the 55 genes exhibiting sigma(B)-dependent expression, 54 were preceded by a sequence resembling the sigma(B) promoter consensus sequence. Rapid amplification of cDNA ends-PCR was used to confirm the sigma(B)-dependent nature of a subset of eight selected promoter regions. Notably, the sigma(B)-dependent L. monocytogenes genes identified through this HMM/microarray strategy included both stress response genes (e.g., gadB, ctc, and the glutathione reductase gene lmo1433) and virulence genes (e.g., inlA, inlB, and bsh). Our data demonstrate that, in addition to regulating expression of genes important for survival under environmental stress conditions, sigma(B) also contributes to regulation of virulence gene expression in L. monocytogenes. These findings strongly suggest that sigma(B) contributes to L. monocytogenes gene expression during infection.

Global transcription analysis of Krebs tricarboxylic acid cycle mutants reveals an alternating pattern of gene expression and effects on hypoxic and oxidative genes

To understand the many roles of the Krebs tricarboxylic acid (TCA) cycle in cell function, we used DNA microarrays to examine gene expression in response to TCA cycle dysfunction. mRNA was analyzed from yeast strains harboring defects in each of 15 genes that encode subunits of the eight TCA cycle enzymes. The expression of >400 genes changed at least threefold in response to TCA cycle dysfunction. Many genes displayed a common response to TCA cycle dysfunction indicative of a shift away from oxidative metabolism. Another set of genes displayed a pairwise, alternating pattern of expression in response to contiguous TCA cycle enzyme defects: expression was elevated in aconitase and isocitrate dehydrogenase mutants, diminished in alpha-ketoglutarate dehydrogenase and succinyl-CoA ligase mutants, elevated again in succinate dehydrogenase and fumarase mutants, and diminished again in malate dehydrogenase and citrate synthase mutants. This pattern correlated with previously defined TCA cycle growth-enhancing mutations and suggested a novel metabolic signaling pathway monitoring TCA cycle function. Expression of hypoxic/anaerobic genes was elevated in alpha-ketoglutarate dehydrogenase mutants, whereas expression of oxidative genes was diminished, consistent with a heme signaling defect caused by inadequate levels of the heme precursor, succinyl-CoA. These studies have revealed extensive responses to changes in TCA cycle function and have uncovered new and unexpected metabolic networks that are wired into the TCA cycle.

Transcriptional mapping and genomic analysis of the cardiac atria and ventricles

The atria and ventricles of the heart have distinct development, structure, and physiology. However, only a few of the genes that underlie the differences between these tissues are known. We used a murine cardiac cDNA microarray to identify genes differentially expressed in the atria and ventricles. The reliability of these findings is supported by highly concordant repetition of hybridization, recognition of previously known atrial and ventricular isoforms of contractile proteins, and confirmation of results by quantitative PCR and in situ hybridization. We examined the most differentially regulated genes for evolutionarily conserved noncoding sequences and found that atrial-expressed genes have more predicted myocyte enhancer factor-2 (MEF2) binding sites than ventricle-predominant genes. We confirmed that messages for MEF2 family members are more abundant in the atria, as are their protein products. Moreover, the activity of a transgenic reporter construct for MEF2 activity is preferentially upregulated in the atria in response to hypertrophic stimuli. This study provides a greater understanding of the molecular differences between atria and ventricles and establishes the framework for an anatomically detailed evaluation of cardiac transcriptional regulation.

Induction of global stress response in Saccharomyces cerevisiae cells lacking telomerase

Cellular senescence is a major intermediate step from healthy cells toward tumor cells. By using microarrays that simultaneously examine the transcription levels of 6,200 Saccharomyces cerevisiae genes, we show that 45 gene transcript levels are increased and 11 are decreased after exposure to telomere shortening and cellular senescence in a telomerase-deficient mutant. About half of the genes that showed increased expression were found induced under stress, consistent with the notion that critical short telomeres cause stress to cells. Surprisingly, the expression level of telomere recombination genes was not altered suggesting that even though recombination is a means to rescue critically short telomeres, its machinery was not controlled by telomere shortening. The expression of telomere-proximal genes was also analyzed. The possibility of induction of a program to cope with cellular senescence and active telomere-telomere recombination is discussed.

RTG-dependent mitochondria to nucleus signaling is negatively regulated by the seven WD-repeat protein Lst8p

In cells with reduced mitochondrial function, RTG1, 2 and 3 are required for expression of genes involved in glutamate synthesis. Glutamate negatively regulates RTG-dependent gene expression upstream of Rtg2p, which, in turn, acts upstream of the bHLH/Zip transcription factors, Rtg1p and Rtg3p. Here we report that some mutations [lst8-(2-5)] in LST8, an essential gene encoding a seven WD40-repeat protein required for targeting of amino acid permeases (AAPs) to the plasma membrane, bypass the requirement for Rtg2p and abolish glutamate repression of RTG-dependent gene expression. The lst8-1 mutation, however, which reduces plasma membrane expression of AAP, cannot bypass the Rtg2p requirement, but still suppresses glutamate repression of RTG target gene expression. We show that Lst8p negatively regulates RTG gene function, acting at two sites, one upstream of Rtg2p, affecting glutamate repression of RTG-dependent gene expression through Ssy1p, an AAP-like sensor of external amino acids, and the other between Rtg2p and Rtg1p-Rtg3p. These data, together with genome-wide transcription profiling, reveal pathways regulated by glutamate, and provide insight into the regulation of cellular responses to mitochondrial dysfunction.

Genome-wide responses to mitochondrial dysfunction

Mitochondrial dysfunction can lead to diverse cellular and organismal responses. We used DNA microarrays to characterize the transcriptional responses to different mitochondrial perturbations in Saccharomyces cerevisiae. We examined respiratory-deficient petite cells and respiratory-competent wild-type cells treated with the inhibitors of oxidative phosphorylation antimycin, carbonyl cyanide m-chlorophenylhydrazone, or oligomycin. We show that respiratory deficiency, but not inhibition of mitochondrial ATP synthesis per se, induces a suite of genes associated with both peroxisomal activities and metabolite-restoration (anaplerotic) pathways that would mitigate the loss of a complete tricarboxylic acid cycle. The array data suggested, and direct microscopic observation of cells expressing a derivative of green fluorescent protein with a peroxisomal matrix-targeting signal confirmed, that respiratory deficiency dramatically induces peroxisome biogenesis. Transcript profiling of cells harboring null alleles of RTG1, RTG2, or RTG3, genes known to control signaling from mitochondria to the nucleus, suggests that there are multiple pathways of cross-talk between these organelles in yeast.