A Predictive Model of the Oxygen and Heme Regulatory Network in Yeast

Deciphering gene regulatory mechanisms through the analysis of high-throughput expression data is a challenging computational problem. Previous computational studies have used large expression datasets in order to resolve fine patterns of coexpression, producing clusters or modules of potentially coregulated genes. These methods typically examine promoter sequence information, such as DNA motifs or transcription factor occupancy data, in a separate step after clustering. We needed an alternative and more integrative approach to study the oxygen regulatory network in Saccharomyces cerevisiae using a small dataset of perturbation experiments. Mechanisms of oxygen sensing and regulation underlie many physiological and pathological processes, and only a handful of oxygen regulators have been identified in previous studies. We used a new machine learning algorithm called MEDUSA to uncover detailed information about the oxygen regulatory network using genome-wide expression changes in response to perturbations in the levels of oxygen, heme, Hap1, and Co2+. MEDUSA integrates mRNA expression, promoter sequence, and ChIP-chip occupancy data to learn a model that accurately predicts the differential expression of target genes in held-out data. We used a novel margin-based score to extract significant condition-specific regulators and assemble a global map of the oxygen sensing and regulatory network. This network includes both known oxygen and heme regulators, such as Hap1, Mga2, Hap4, and Upc2, as well as many new candidate regulators. MEDUSA also identified many DNA motifs that are consistent with previous experimentally identified transcription factor binding sites. Because MEDUSA's regulatory program associates regulators to target genes through their promoter sequences, we directly tested the predicted regulators for OLE1, a gene specifically induced under hypoxia, by experimental analysis of the activity of its promoter. In each case, deletion of the candidate regulator resulted in the predicted effect on promoter activity, confirming that several novel regulators identified by MEDUSA are indeed involved in oxygen regulation. MEDUSA can reveal important information from a small dataset and generate testable hypotheses for further experimental analysis. Supplemental data are included

A sequence that directs transcriptional initiation in yeast

While RNA polymerase II of the yeast Saccharomyces cerevisiae initiates transcription at discrete sites, these sites are located over a wide range of distances from the TATA box for different genes. This variability has led to a number of proposals for consensus sequences located at the initiation site which, in conjunction with the TATA box, would direct initiation. We tested this hypothesis via oligonucleotide-directed mutagenesis, by placing the sequence CAAG, a member of one of these consensus sequences, upstream of the coding sequence of the CYC7 gene at a site at which initiation does not occur. The distance between the TATA sequence and this putative initiation site was varied by inserting it into the wild-type gene and three deletion mutants. The results demonstrated that this sequence can serve as an initiation site when located 49, 77, or 106 bp from the TATA sequence, but not when located 30 bp away.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/46971/1/294_2004_Article_BF00312597.pd

Genome-Wide Fitness and Expression Profiling Implicate Mga2 in Adaptation to Hydrogen Peroxide

Caloric restriction extends lifespan, an effect once thought to involve attenuation of reactive oxygen species (ROS) generated by aerobic metabolism. However, recent evidence suggests that caloric restriction may in fact raise ROS levels, which in turn provides protection from acute doses of oxidant through a process called adaptation. To shed light on the molecular mechanisms of adaptation, we designed a series of genome-wide deletion fitness and mRNA expression screens to identify genes involved in adaptation to hydrogen peroxide. Combined with known transcriptional interactions, the integrated data implicate Yap1 and Skn7 as central transcription factors of both the adaptive and acute oxidative responses. They also identify the transcription factors Mga2 and Rox1 as active exclusively in the adaptive response and show that Mga2 is essential for adaptation. These findings are striking because Mga2 and Rox1 have been thought to control the response to hypoxic, not oxidative, conditions. Expression profiling of mga2Δ and rox1Δ knockouts shows that these factors most strongly regulate targets in ergosterol, fatty-acid, and zinc metabolic pathways. Direct quantitation of ergosterol reveals that its basal concentration indeed depends on Mga2, but that Mga2 is not required for the decrease in ergosterol observed during adaptation

Extent of Structural Asymmetry in Homodimeric Proteins: Prevalence and Relevance

Most homodimeric proteins have symmetric structure. Although symmetry is known to confer structural and functional advantage, asymmetric organization is also observed. Using a non-redundant dataset of 223 high-resolution crystal structures of biologically relevant homodimers, we address questions on the prevalence and significance of asymmetry. We used two measures to quantify global and interface asymmetry, and assess the correlation of several molecular and structural parameters with asymmetry. We have identified rare cases (11/223) of biologically relevant homodimers with pronounced global asymmetry. Asymmetry serves as a means to bring about 2∶1 binding between the homodimer and another molecule; it also enables cellular signalling arising from asymmetric macromolecular ligands such as DNA. Analysis of these cases reveals two possible mechanisms by which possible infinite array formation is prevented. In case of homodimers associating via non-topologically equivalent surfaces in their tertiary structures, ligand-dependent mechanisms are used. For stable dimers binding via large surfaces, ligand-dependent structural change regulates polymerisation/depolymerisation; for unstable dimers binding via smaller surfaces that are not evolutionarily well conserved, dimerisation occurs only in the presence of the ligand. In case of homodimers associating via interaction surfaces with parts of the surfaces topologically equivalent in the tertiary structures, steric hindrance serves as the preventive mechanism of infinite array. We also find that homodimers exhibiting grossly symmetric organization rarely exhibit either perfect local symmetry or high local asymmetry. Binding of small ligands at the interface does not cause any significant variation in interface asymmetry. However, identification of biologically relevant interface asymmetry in grossly symmetric homodimers is confounded by the presence of similar small magnitude changes caused due to artefacts of crystallisation. Our study provides new insights regarding accommodation of asymmetry in homodimers

Using RNA-seq to determine the transcriptional landscape and the hypoxic response of the pathogenic yeast Candida parapsilosis

<p>Abstract</p> <p>Background</p> <p><it>Candida parapsilosis </it>is one of the most common causes of <it>Candida </it>infection worldwide. However, the genome sequence annotation was made without experimental validation and little is known about the transcriptional landscape. The transcriptional response of <it>C. parapsilosis </it>to hypoxic (low oxygen) conditions, such as those encountered in the host, is also relatively unexplored.</p> <p>Results</p> <p>We used next generation sequencing (RNA-seq) to determine the transcriptional profile of <it>C. parapsilosis </it>growing in several conditions including different media, temperatures and oxygen concentrations. We identified 395 novel protein-coding sequences that had not previously been annotated. We removed > 300 unsupported gene models, and corrected approximately 900. We mapped the 5' and 3' UTR for thousands of genes. We also identified 422 introns, including two introns in the 3' UTR of one gene. This is the first report of 3' UTR introns in the Saccharomycotina. Comparing the introns in coding sequences with other species shows that small numbers have been gained and lost throughout evolution. Our analysis also identified a number of novel transcriptional active regions (nTARs). We used both RNA-seq and microarray analysis to determine the transcriptional profile of cells grown in normoxic and hypoxic conditions in rich media, and we showed that there was a high correlation between the approaches. We also generated a knockout of the <it>UPC2 </it>transcriptional regulator, and we found that similar to <it>C. albicans</it>, Upc2 is required for conferring resistance to azole drugs, and for regulation of expression of the ergosterol pathway in hypoxia.</p> <p>Conclusion</p> <p>We provide the first detailed annotation of the <it>C. parapsilosis </it>genome, based on gene predictions and transcriptional analysis. We identified a number of novel ORFs and other transcribed regions, and detected transcripts from approximately 90% of the annotated protein coding genes. We found that the transcription factor Upc2 role has a conserved role as a major regulator of the hypoxic response in <it>C. parapsilosis </it>and <it>C. albicans</it>.</p

Oxygen dependence of metabolic fluxes and energy generation of Saccharomyces cerevisiae CEN.PK113-1A

<p>Abstract</p> <p>Background</p> <p>The yeast <it>Saccharomyces cerevisiae </it>is able to adjust to external oxygen availability by utilizing both respirative and fermentative metabolic modes. Adjusting the metabolic mode involves alteration of the intracellular metabolic fluxes that are determined by the cell's multilevel regulatory network. Oxygen is a major determinant of the physiology of <it>S. cerevisiae </it>but understanding of the oxygen dependence of intracellular flux distributions is still scarce.</p> <p>Results</p> <p>Metabolic flux distributions of <it>S. cerevisiae </it>CEN.PK113-1A growing in glucose-limited chemostat cultures at a dilution rate of 0.1 h^-1with 20.9%, 2.8%, 1.0%, 0.5% or 0.0% O₂in the inlet gas were quantified by ¹³C-MFA. Metabolic flux ratios from fractional [U-¹³C]glucose labelling experiments were used to solve the underdetermined MFA system of central carbon metabolism of <it>S. cerevisiae</it>.</p> <p>While ethanol production was observed already in 2.8% oxygen, only minor differences in the flux distribution were observed, compared to fully aerobic conditions. However, in 1.0% and 0.5% oxygen the respiratory rate was severely restricted, resulting in progressively reduced fluxes through the TCA cycle and the direction of major fluxes to the fermentative pathway. A redistribution of fluxes was observed in all branching points of central carbon metabolism. Yet only when oxygen provision was reduced to 0.5%, was the biomass yield exceeded by the yields of ethanol and CO₂. Respirative ATP generation provided 59% of the ATP demand in fully aerobic conditions and still a substantial 25% in 0.5% oxygenation. An extensive redistribution of fluxes was observed in anaerobic conditions compared to all the aerobic conditions. Positive correlation between the transcriptional levels of metabolic enzymes and the corresponding fluxes in the different oxygenation conditions was found only in the respirative pathway.</p> <p>Conclusion</p> <p>¹³C-constrained MFA enabled quantitative determination of intracellular fluxes in conditions of different redox challenges without including redox cofactors in metabolite mass balances. A redistribution of fluxes was observed not only for respirative, respiro-fermentative and fermentative metabolisms, but also for cells grown with 2.8%, 1.0% and 0.5% oxygen. Although the cellular metabolism was respiro-fermentative in each of these low oxygen conditions, the actual amount of oxygen available resulted in different contributions through respirative and fermentative pathways.</p