Expression and functional profiling reveal distinct gene classes involved in fatty acid metabolism

Peroxisome prognostications: Exploring the birth, life, and death of an organelle

Peroxisomes play a central role in human health and have biochemical properties that promote their use in many biotechnology settings. With a primary role in lipid metabolism, peroxisomes share a niche with lipid droplets within the endomembrane-secretory system. Notably, factors in the ER required for the biogenesis of peroxisomes also impact the formation of lipid droplets. The dynamic interface between peroxisomes and lipid droplets, and also between these organelles and the ER and mitochondria, controls their metabolic flux and their dynamics. Here, we review our understanding of peroxisome biogenesis to propose and reframe models for understanding how peroxisomes are formed in cells. To more fully understand the roles of peroxisomes and to take advantage of their many properties that may prove useful in novel therapeutics or biotechnology applications, we recast mechanisms controlling peroxisome biogenesis in a framework that integrates inference from these models with experimental data.

A decade of advances in transposon-insertion sequencing

It has been 10 years since the introduction of modern transposon-insertion sequencing (TIS) methods, which combine genome-wide transposon mutagenesis with high-throughput sequencing to estimate the fitness contribution or essentiality of each genetic component in a bacterial genome. Four TIS variations were published in 2009: transposon sequencing (Tn-Seq), transposon-directed insertion site sequencing (TraDIS), insertion sequencing (INSeq) and high-throughput insertion tracking by deep sequencing (HITS). TIS has since become an important tool for molecular microbiologists, being one of the few genome-wide techniques that directly links phenotype to genotype and ultimately can assign gene function. In this Review, we discuss the recent applications of TIS to answer overarching biological questions. We explore emerging and multidisciplinary methods that build on TIS, with an eye towards future applications.

In this Review, several experts discuss progress in the decade since the development of transposon-based approaches for bacterial genetic screens. They describe how advances in both experimental technologies and analytical strategies are resulting in insights into diverse biological processes.

Genetic determinants of Pseudomonas aeruginosa fitness during biofilm growth

Pseudomonas aeruginosa is an environmental bacterium and an opportunistic human pathogen. It is also a well-established model organism to study bacterial adaptation to stressful conditions, such as those encountered during an infection process in the human host. Advancing knowledge on P. aeruginosa adaptation to biofilm growth conditions is bound to reveal novel strategies and targets for the treatment of chronic biofilm-associated infections. Here, we generated transposon insertion libraries in three P. aeruginosa strain backgrounds and determined the relative frequency of each insertion following biofilm growth using transposon sequencing. We demonstrate that in general the SOS response, several tRNA modifying enzymes as well as adaptation to microaerophilic growth conditions play a key role in bacterial survival under biofilm growth conditions. On the other hand, presence of genes involved in motility and PQS signaling were less important during biofilm growth. Several mutants exhibiting transposon insertions in genes detected in our screen were validated for their biofilm growth capabilities and biofilm specific transcriptional responses using independently generated transposon mutants. Our results provide new insights into P. aeruginosa adaptation to biofilm growth conditions. The detection of previously unknown determinants of biofilm survival supports the use of transposon insertion sequencing as a global genomic technology for understanding the establishment of difficult to treat biofilm-associated infections.

ESCRT-III is required for scissioning new peroxisomes from the endoplasmic reticulum

Dynamic control of peroxisome proliferation is integral to the peroxisome's many functions. The endoplasmic reticulum (ER) serves as a source of preperoxisomal vesicles (PPVs) that mature into peroxisomes during de novo peroxisome biogenesis and support growth and division of existing peroxisomes. However, the mechanism of PPV formation and release from the ER remains poorly understood. In this study, we show that endosomal sorting complexes required for transport (ESCRT)-III are required to release PPVs budding from the ER into the cytosol. Absence of ESCRT-III proteins impedes de novo peroxisome formation and results in an aberrant peroxisome population in vivo. Using a cell-free PPV budding assay, we show that ESCRT-III proteins Vps20 and Snf7 are necessary to release PPVs from the ER. ESCRT-III is therefore a positive effector of membrane scission for vesicles budding both away from and toward the cytosol. These findings have important implications for the evolutionary timing of emergence of peroxisomes and the rest of the internal membrane architecture of the eukaryotic cell.

Functional genomics of lipid metabolism in the oleaginous yeast Rhodosporidium toruloides

The basidiomycete yeast Rhodosporidium toruloides (also known as Rhodotorula toruloides) accumulates high concentrations of lipids and carotenoids from diverse carbon sources. It has great potential as a model for the cellular biology of lipid droplets and for sustainable chemical production. We developed a method for high-throughput genetics (RB-TDNAseq), using sequence-barcoded Agrobacterium tumefaciens T-DNA insertions. We identified 1,337 putative essential genes with low T-DNA insertion rates. We functionally profiled genes required for fatty acid catabolism and lipid accumulation, validating results with 35 targeted deletion strains. We identified a high-confidence set of 150 genes affecting lipid accumulation, including genes with predicted function in signaling cascades, gene expression, protein modification and vesicular trafficking, autophagy, amino acid synthesis and tRNA modification, and genes of unknown function. These results greatly advance our understanding of lipid metabolism in this oleaginous species and demonstrate a general approach for barcoded mutagenesis that should enable functional genomics in diverse fungi.

The fungus Rhodosporidium toruloides can grow on substances extracted from plant matter that is inedible to humans such as corn stalks, wood pulp, and grasses. Under some growth conditions, the fungus can accumulate massive stores of hydrocarbon-rich fats and pigments. A community of scientists and engineers has begun genetically modifying R. toruloides to convert these naturally produced fats and pigments into fuels, chemicals and medicines. These could form sustainable replacements for products made from petroleum or harvested from threatened animal and plant species.

Fungi, plants, animals and other eukaryotes store fat in specialized compartments called lipid droplets. The genes that control the metabolism – the production, use and storage – of fat in lipid bodies have been studied in certain eukaryotes, including species of yeast. However, R. toruloides is only distantly related to the most well-studied of these species. This means that we cannot be certain that a gene will play the same role in R. toruloides as in those species.

To assemble the most comprehensive list possible of the genes in R. toruloides that affect the production, use, or storage of fat in lipid bodies, Coradetti, Pinel et al. constructed a population of hundreds of thousands of mutant fungal strains, each with its own unique DNA ‘barcode’. The effects that mutations in over 6,000 genes had on growth and fat accumulation in these fungi were measured simultaneously in several experiments. This general approach is not new, but technical limitations had, until now, restricted its use in fungi to a few species.

Coradetti, Pinel et al. identified hundreds of genes that affected the ability of R. toruloides to metabolise fat. Many of these genes were related to genes with known roles in fat metabolism in other eukaryotes. Other genes are involved in different cell processes, such as the recycling of waste products in the cell. Their identification adds weight to the view that the links between these cellular processes and fat metabolism are deep and widespread amongst eukaryotes. Finally, some of the genes identified by Coradetti, Pinel et al. are not closely related to any well-studied genes. Further study of these genes could help us to understand why R. toruloides can accumulate much larger amounts of fat than most other fungi.

The methods developed by Coradetti, Pinel et al. should be possible to implement in many species of fungi. As a result these techniques may eventually contribute to the development of new treatments for human fungal diseases, the protection of important food crops, and a deeper understanding of the roles various fungi play in the broader ecosystem.

The functional readthrough extension of malate dehydrogenase reveals a modification of the genetic code

Translational readthrough gives rise to C-terminally extended proteins, thereby providing the cell with new protein isoforms. These may have different properties from the parental proteins if the extensions contain functional domains. While for most genes amino acid incorporation at the stop codon is far lower than 0.1%, about 4% of malate dehydrogenase (MDH1) is physiologically extended by translational readthrough and the actual ratio of MDH1x (extended protein) to ‘normal' MDH1 is dependent on the cell type. In human cells, arginine and tryptophan are co-encoded by the MDH1x UGA stop codon. Readthrough is controlled by the 7-nucleotide high-readthrough stop codon context without contribution of the subsequent 50 nucleotides encoding the extension. All vertebrate MDH1x is directed to peroxisomes via a hidden peroxisomal targeting signal (PTS) in the readthrough extension, which is more highly conserved than the extension of lactate dehydrogenase B. The hidden PTS of non-mammalian MDH1x evolved to be more efficient than the PTS of mammalian MDH1x. These results provide insight into the genetic and functional co-evolution of these dually localized dehydrogenases.

Antibiotics Disrupt Coordination between Transcriptional and Phenotypic Stress Responses in Pathogenic Bacteria

Bacterial genes that change in expression upon environmental disturbance have commonly been seen as those that must also phenotypically matter. However, several studies suggest that differentially expressed genes are rarely phenotypically important. We demonstrate, for Gram-positive and Gram-negative bacteria, that these seemingly uncoordinated gene sets are involved in responses that can be linked through topological network analysis. However, the level of coordination is stress dependent. While a well-coordinated response is triggered in response to nutrient stress, antibiotics trigger an uncoordinated response in which transcriptionally and phenotypically important genes are neither linked spatially nor in their magnitude. Moreover, a gene expression meta-analysis reveals that genes with large fitness changes during stress have low transcriptional variation across hundreds of other conditions, and vice versa. Our work suggests that cellular responses can be understood through network models that incorporate regulatory and genetic relationships, which could aid drug target predictions and genetic network engineering.

Pex35 is a regulator of peroxisome abundance

Peroxisomes are cellular organelles with vital functions in lipid, amino acid and redox metabolism. The cellular formation and dynamics of peroxisomes are governed by PEX genes; however, the regulation of peroxisome abundance is still poorly understood. Here, we use a high-content microscopy screen in Saccharomyces cerevisiae to identify new regulators of peroxisome size and abundance. Our screen led to the identification of a previously uncharacterized gene, which we term PEX35, which affects peroxisome abundance. PEX35 encodes a peroxisomal membrane protein, a remote homolog to several curvature-generating human proteins. We systematically characterized the genetic and physical interactome as well as the metabolome of mutants in PEX35, and we found that Pex35 functionally interacts with the vesicle-budding-inducer Arf1. Our results highlight the functional interaction between peroxisomes and the secretory pathway.

Pex17p-dependent assembly of Pex14p/Dyn2p-subcomplexes of the peroxisomal protein import machinery

Peroxisomal matrix protein import is facilitated by cycling receptors that recognize their cargo proteins in the cytosol by peroxisomal targeting sequences (PTS). In the following, the assembled receptor-cargo complex is targeted to the peroxisomal membrane where it docks to the docking-complex as part of the peroxisomal translocation machinery. The docking-complex is composed of Pex13p, Pex14p and in yeast also Pex17p, whose function is still elusive. In order to characterize the function of Pex17p, we compared the composition and size of peroxisomal receptor-docking complexes from wild-type and pex17Δ cells. Our data demonstrate that the deficiency of Pex17p affects the stoichiometry of the constituents of an isolated 600kDa complex and that pex17Δ cells lack a high molecular weight complex (>900kDa) of unknown function. We identified the dynein light chain protein Dyn2p as an additional core component of the Pex14p/Pex17p-complex. Both, Pex14p and Pex17p interact directly with Dyn2p, but in vivo, Pex17p turned out to be prerequisite for an association of Dyn2p with Pex14p. Finally, like pex17Δ also dyn2Δ cells lack the high molecular weight complex. As dyn2Δ cells also display reduced peroxisomal function, our data indicate that Dyn2p-dependent formation of the high molecular weight Pex14p-complex is required to maintain peroxisomal function on wild-type level.

Effect of ω-3 supplementation on placental lipid metabolism in overweight and obese women

BACKGROUND

The placentas of obese women accumulate lipids that may alter fetal lipid exposure. The long-chain omega-3 fatty acids (n–3 FAs) docosahexaenoic acid (DHA) and eicosapentaenoic acid (EPA) alter FA metabolism in hepatocytes, although their effect on the placenta is poorly understood.

OBJECTIVE

We aimed to investigate whether n–3 supplementation during pregnancy affects lipid metabolism in the placentas of overweight and obese women at term.

DESIGN

A secondary analysis of a double-blind randomized controlled trial was conducted in healthy overweight and obese pregnant women who were randomly assigned to DHA plus EPA (2 g/d) or placebo twice a day from early pregnancy to term. Placental FA uptake, esterification, and oxidation pathways were studied by measuring the expression of key genes in the placental tissue of women supplemented with placebo and n–3 and in vitro in isolated trophoblast cells in response to DHA and EPA treatment.

RESULTS

Total lipid content was significantly lower in the placentas of overweight and obese women supplemented with n–3 FAs than in those supplemented with placebo (14.14 ± 1.03 compared with 19.63 ± 1.45 mg lipid/g tissue; P < 0.05). The messenger RNA expression of placental FA synthase (FAS) and diacylglycerol O-acyltransferase 1 (DGAT1) was negatively correlated with maternal plasma enrichment in DHA and EPA (P < 0.05). The expression of placental peroxisome proliferator–activated receptor γ (r = −0.39, P = 0.04) and its target genes DGAT1 (r = −0.37, P = 0.02) and PLIN2 (r = −0.38, P = 0.04) significantly decreased, with an increasing maternal n–3:n–6 ratio (representing the n–3 status) near the end of pregnancy. The expression of genes that regulate FA oxidation or uptake was not changed. Birth weight and length were significantly higher in the offspring of n–3-supplemented women than in those in the placebo group (P < 0.05), but no differences in the ponderal index were observed. Supplementation of n–3 significantly decreased FA esterification in isolated trophoblasts without affecting FA oxidation.

CONCLUSION

Supplementing overweight and obese women with n–3 FAs during pregnancy inhibited the ability of the placenta to esterify and store lipids. This trial was registered at clinicaltrials.gov as NCT00957476.

Peroxisomes take shape

Peroxisomes carry out various oxidative reactions that are tightly regulated to adapt to the changing needs of the cell and varying external environments. Accordingly, they are remarkably fluid and can change dramatically in abundance, size, shape and content in response to numerous cues. These dynamics are controlled by multiple aspects of peroxisome biogenesis that are coordinately regulated with each other and with other cellular processes. Ongoing studies are deciphering the diverse molecular mechanisms that underlie biogenesis and how they cooperate to dynamically control peroxisome utility. These important challenges should lead to an understanding of peroxisome dynamics that can be capitalized upon for bioengineering and the development of therapies to improve human health.

Molecular mechanisms of system responses to novel stimuli are predictable from public data

Systems scale models provide the foundation for an effective iterative cycle between hypothesis generation, experiment and model refinement. Such models also enable predictions facilitating the understanding of biological complexity and the control of biological systems. Here, we demonstrate the reconstruction of a globally predictive gene regulatory model from public data: a model that can drive rational experiment design and reveal new regulatory mechanisms underlying responses to novel environments. Specifically, using ∼ 1500 publically available genome-wide transcriptome data sets from Saccharomyces cerevisiae, we have reconstructed an environment and gene regulatory influence network that accurately predicts regulatory mechanisms and gene expression changes on exposure of cells to completely novel environments. Focusing on transcriptional networks that induce peroxisomes biogenesis, the model-guided experiments allow us to expand a core regulatory network to include novel transcriptional influences and linkage across signaling and transcription. Thus, the approach and model provides a multi-scalar picture of gene dynamics and are powerful resources for exploiting extant data to rationally guide experimentation. The techniques outlined here are generally applicable to any biological system, which is especially important when experimental systems are challenging and samples are difficult and expensive to obtain-a common problem in laboratory animal and human studies.

Dynein light chain interaction with the peroxisomal import docking complex modulates peroxisome biogenesis in yeast

Dynein is a large macromolecular motor complex that moves cargo along microtubules. A motor-independent role for the light chain of dynein, Dyn2p, in peroxisome biology in Saccharomyces cerevisiae was suggested from its interaction with Pex14p, a component of the peroxisomal matrix protein import docking complex. Here we show that cells of the yeast Yarrowia lipolytica deleted for the gene encoding the homologue of Dyn2p are impaired in peroxisome function and biogenesis. These cells exhibit compromised growth on medium containing oleic acid as the carbon source, the metabolism of which requires functional peroxisomes. Their peroxisomes have abnormal morphology, atypical matrix protein localization, and an absence of proteolytic processing of the matrix enzyme thiolase, which normally occurs upon its import into the peroxisome. We also show physical and genetic interactions between Dyn2p and members of the docking complex, particularly Pex17p. Together, our results demonstrate a role for Dyn2p in the assembly of functional peroxisomes and provide evidence that Dyn2p acts in cooperation with the peroxisomal matrix protein import docking complex to effect optimal matrix protein import.

Indirect and suboptimal control of gene expression is widespread in bacteria

Gene regulation in bacteria is usually described as an adaptive response to an environmental change so that genes are expressed when they are required. We instead propose that most genes are under indirect control: their expression responds to signal(s) that are not directly related to the genes' function. Indirect control should perform poorly in artificial conditions, and we show that gene regulation is often maladaptive in the laboratory. In Shewanella oneidensis MR-1, 24% of genes are detrimental to fitness in some conditions, and detrimental genes tend to be highly expressed instead of being repressed when not needed. In diverse bacteria, there is little correlation between when genes are important for optimal growth or fitness and when those genes are upregulated. Two common types of indirect control are constitutive expression and regulation by growth rate; these occur for genes with diverse functions and often seem to be suboptimal. Because genes that have closely related functions can have dissimilar expression patterns, regulation may be suboptimal in the wild as well as in the laboratory.

Global analysis of condition-specific subcellular protein distribution and abundance

Cellular control of protein activities by modulation of their abundance or compartmentalization is not easily measured on a large scale. We developed and applied a method to globally interrogate these processes that is widely useful for systems-level analyses of dynamic cellular responses in many cell types. The approach involves subcellular fractionation followed by comprehensive proteomic analysis of the fractions, which is enabled by a data-independent acquisition mass spectrometry approach that samples every available mass to charge channel systematically to maximize sensitivity. Next, various fraction-enrichment ratios are measured for all detected proteins across different environmental conditions and used to group proteins into clusters reflecting changes in compartmentalization and relative conditional abundance. Application of the approach to characterize the response of yeast proteins to fatty acid exposure revealed dynamics of peroxisomes and novel dynamics of MCC/eisosomes, specialized plasma membrane domains comprised of membrane compartment occupied by Can1 (MCC) and eisosome subdomains. It also led to the identification of Fat3, a fatty acid transport protein of the plasma membrane, previously annotated as Ykl187.

Role of the repressor Oaf3p in the recruitment of transcription factors and chromatin dynamics during the oleate response

Cellular responses to environmental stimuli are mediated by the co-ordinated activity of multiple control mechanisms, which result in the dynamics of cell function. Communication between different levels of regulation is central for this adaptability. The present study focuses on the interplay between transcriptional regulators and chromatin modifiers to co-operatively regulate transcription in response to a fatty acid stimulus. The genes involved in the β-oxidation of fatty acids are highly induced in response to fatty acid exposure by four gene-specific transcriptional regulators, Oaf (oleate-activated transcription factor) 1p, Pip2p (peroxisome induction pathway 2), Oaf3p and Adr1p (alcohol dehydrogenase regulator 1). In the present study, we examine the interplay of these factors with Htz1p (histone variant H2A.Z) in regulating POT1 (peroxisomal oxoacyl thiolase 1) encoding peroxisomal thiolase and PIP2 encoding the autoregulatory oleate-specific transcriptional activator. Temporal resolution of ChIP (chromatin immunoprecipitation) data indicates that Htz1p is required for the timely removal of the transcriptional repressor Oaf3p during oleate induction. Adr1p plays an important role in the assembly of Htz1p-containing nucleosomes on the POT1 and PIP2 promoters. We also investigated the function of the uncharacterized transcriptional inhibitor Oaf3p. Deletion of OAF3 led to faster POT1 mRNA accumulation than in the wild-type. Most impressively, a highly protected nucleosome structure on the POT1 promoter during activation was observed in the OAF3 mutant cells in response to oleate induction.

Exploring the yeast acetylome using functional genomics

Lysine acetylation is a dynamic posttranslational modification with a well-defined role in regulating histones. The impact of acetylation on other cellular functions remains relatively uncharacterized. We explored the budding yeast acetylome with a functional genomics approach, assessing the effects of gene overexpression in the absence of lysine deacetylases (KDACs). We generated a network of 463 synthetic dosage lethal (SDL) interactions involving class I and II KDACs, revealing many cellular pathways regulated by different KDACs. A biochemical survey of genes interacting with the KDAC RPD3 identified 72 proteins acetylated in vivo. In-depth analysis of one of these proteins, Swi4, revealed a role for acetylation in G1-specific gene expression. Acetylation of Swi4 regulates interaction with its partner Swi6, both components of the SBF transcription factor. This study expands our view of the yeast acetylome, demonstrates the utility of functional genomic screens for exploring enzymatic pathways, and provides functional information that can be mined for future studies.

Asymmetric positive feedback loops reliably control biological responses

Positive feedback is a common mechanism enabling biological systems to respond to stimuli in a switch-like manner. Such systems are often characterized by the requisite formation of a heterodimer where only one of the pair is subject to feedback. This ASymmetric Self-UpREgulation (ASSURE) motif is central to many biological systems, including cholesterol homeostasis (LXRα/RXRα), adipocyte differentiation (PPARγ/RXRα), development and differentiation (RAR/RXR), myogenesis (MyoD/E12) and cellular antiviral defense (IRF3/IRF7). To understand why this motif is so prevalent, we examined its properties in an evolutionarily conserved transcriptional regulatory network in yeast (Oaf1p/Pip2p). We demonstrate that the asymmetry in positive feedback confers a competitive advantage and allows the system to robustly increase its responsiveness while precisely tuning the response to a consistent level in the presence of varying stimuli. This study reveals evolutionary advantages for the ASSURE motif, and mechanisms for control, that are relevant to pharmacologic intervention and synthetic biology applications.

A regression model approach to enable cell morphology correction in high-throughput flow cytometry

Cells exposed to stimuli exhibit a wide range of responses ensuring phenotypic variability across the population. Such single cell behavior is often examined by flow cytometry; however, gating procedures typically employed to select a small subpopulation of cells with similar morphological characteristics make it difficult, even impossible, to quantitatively compare cells across a large variety of experimental conditions because these conditions can lead to profound morphological variations. To overcome these limitations, we developed a regression approach to correct for variability in fluorescence intensity due to differences in cell size and granularity without discarding any of the cells, which gating ipso facto does. This approach enables quantitative studies of cellular heterogeneity and transcriptional noise in high-throughput experiments involving thousands of samples. We used this approach to analyze a library of yeast knockout strains and reveal genes required for the population to establish a bimodal response to oleic acid induction. We identify a group of epigenetic regulators and nucleoporins that, by maintaining an 'unresponsive population,' may provide the population with the advantage of diversified bet hedging.

Extensive divergence of yeast stress responses through transitions between induced and constitutive activation

Closely related species show a high degree of differences in gene expression, but the functional significance of these differences remains unclear. Similarly, stress responses in yeast typically involve differential expression of numerous genes, and it is unclear how many of these are functionally significant. To address these issues, we compared the expression programs of four yeast species under different growth conditions, and found that the response of these species to stress has diverged extensively. On an individual gene basis, most transcriptional responses are not conserved in any pair of species, and there are very limited common responses among all four species. We present evidence that many evolutionary changes in stress responses are compensated either (i) by the response of related genes or (ii) by changes in the basal expression levels of the genes whose responses have diverged. Thus, stress-related genes are often induced upon stress in some species but maintain high levels even in the absence of stress at other species, indicating a transition between induced and constitutive activation. In addition, ~15% of the stress responses are specific to only one of the four species, with no evidence for compensating effects or stress-related annotations, and these may reflect fortuitous regulation that is unimportant for the stress response (i.e., biological noise). Frequent compensatory changes and biological noise may explain how diverged expression responses support similar physiological responses.

Trade-off between responsiveness and noise suppression in biomolecular system responses to environmental cues

When living systems detect changes in their external environment their response must be measured to balance the need to react appropriately with the need to remain stable, ignoring insignificant signals. Because this is a fundamental challenge of all biological systems that execute programs in response to stimuli, we developed a generalized time-frequency analysis (TFA) framework to systematically explore the dynamical properties of biomolecular networks. Using TFA, we focused on two well-characterized yeast gene regulatory networks responsive to carbon-source shifts and a mammalian innate immune regulatory network responsive to lipopolysaccharides (LPS). The networks are comprised of two different basic architectures. Dual positive and negative feedback loops make up the yeast galactose network; whereas overlapping positive and negative feed-forward loops are common to the yeast fatty-acid response network and the LPS-induced network of macrophages. TFA revealed remarkably distinct network behaviors in terms of trade-offs in responsiveness and noise suppression that are appropriately tuned to each biological response. The wild type galactose network was found to be highly responsive while the oleate network has greater noise suppression ability. The LPS network appeared more balanced, exhibiting less bias toward noise suppression or responsiveness. Exploration of the network parameter space exposed dramatic differences in system behaviors for each network. These studies highlight fundamental structural and dynamical principles that underlie each network, reveal constrained parameters of positive and negative feedback and feed-forward strengths that tune the networks appropriately for their respective biological roles, and demonstrate the general utility of the TFA approach for systems and synthetic biology.

Biological systems constantly balance noise suppression with responsiveness. In a fluctuating environment, some changes are insignificant to living cells while others represent cues to which they must respond. These stimuli are interpreted by molecular circuits that enable the cell to strike an appropriate balance between responsiveness and noise suppression. This trade-off is governed by the structure and kinetic parameters of molecular networks, which have been tuned by evolutionary selection for different stimuli and responses. We consider three regulatory circuits (two from yeast and one from mammalian cells), which respond to different environments and involve very different physiological processes. To investigate the responses to a time varying signal, we developed a generalized time-frequency analysis framework for studying such trade-offs using mathematical models of regulatory circuits and explore how the structure and parameters of the circuit affect the trade-offs between noise suppression and responsiveness. The generalized TFA approach represents an effective tool for exploring and analyzing different systems-level dynamical properties. Making use of such properties can facilitate prediction and network control for systems- and synthetic biology applications.

ResponseNet: revealing signaling and regulatory networks linking genetic and transcriptomic screening data

Cellular response to stimuli is typically complex and involves both regulatory and metabolic processes. Large-scale experimental efforts to identify components of these processes often comprise of genetic screening and transcriptomic profiling assays. We previously established that in yeast genetic screens tend to identify response regulators, while transcriptomic profiling assays tend to identify components of metabolic processes. ResponseNet is a network-optimization approach that integrates the results from these assays with data of known molecular interactions. Specifically, ResponseNet identifies a high-probability sub-network, composed of signaling and regulatory molecular interaction paths, through which putative response regulators may lead to the measured transcriptomic changes. Computationally, this is achieved by formulating a minimum-cost flow optimization problem and solving it efficiently using linear programming tools. The ResponseNet web server offers a simple interface for applying ResponseNet. Users can upload weighted lists of proteins and genes and obtain a sparse, weighted, molecular interaction sub-network connecting their data. The predicted sub-network and its gene ontology enrichment analysis are presented graphically or as text. Consequently, the ResponseNet web server enables researchers that were previously limited to separate analysis of their distinct, large-scale experiments, to meaningfully integrate their data and substantially expand their understanding of the underlying cellular response. ResponseNet is available at http://bioinfo.bgu.ac.il/respnet.

Environment-responsive transcription factors bind subtelomeric elements and regulate gene silencing

Chromosome position analysis of ChIP-chip data revealed that several carbon source and stress-responsive yeast transcription factors conditionally bind subtelomeric X elements.

Integration of several microarray gene expression data sets showed that, in this context, the factors conditionally control the boundaries and strength of subtelomeric silencing.

Regulation of silencing by a fatty acid-responsive factor was found to be dependent on Sir2p and independent of Hda1p.

These findings provide a critical link for establishing the mechanisms by which telomere biology is coordinated with other cellular processes including responses to environmental stimuli, aging and adaptation.

It is well established that environmental conditions modulate gene expression through local binding of a variety of conditionally active transcription factors, each responsive to specific environmental cues. However, another prevalent mechanism of gene regulation in eukaryotic cells is the long-range control of groups of genes by chromatin modifications or other position-dependent mechanisms. One such phenomenon, gene silencing, is an important and evolutionarily conserved mode of regulation that controls expression of subtelomeric genes. These genes are enriched for stress response and metabolic genes and their regulation is controlled by the spreading of silencing molecules from chromosome ends (telomeres) into subtelomeric regions. Levels of subtelomeric silencing have been linked to cellular lifespan, and study of the regulation of silencing is fundamental to our understanding of human aging. The spread of silencing in subtelomeric regions is discontinuous, and is controlled by various genomic elements that can either relay and enhance silencing from telomeres (proto-silencing) or create boundaries that protect some genomic regions from silencing. In yeast, every subtelomeric region contains an X element that proto-silences centromere-proximal genes, and also insulates telomere-proximal genes from silencing.

In this paper, we identify a regulatory mechanism to control X element-mediated proto-silencing and insulating activities in response to environmental cues. The mechanism was identified using chromosome position analysis of microarray-based chromatin immunoprecipitation (ChIP-chip) data for environment-responsive TFs and genome-wide gene expression data under the same conditions. The mechanism involves the conditional association of environment-responsive transcription factors to X elements. The binding at X elements results in regulation of proto-silencing of centromere-proximal genes, or insulation of telomere-proximal genes (depending on the factor) in response to environmental stimuli related to stress response and metabolism. One example is shown below (Figure 4B). Transcription factor, Oaf1p, conditionally binds X elements in the presence of fatty acids and enhances proto-silencing specifically under this condition. Oaf1p and several other factors implicated here are known to control adjacent genes at intrachromosomal positions, suggesting their dual functionality in both gene-specific transcriptional regulation, and long-range position-dependent mechanism. Investigation of this mechanism during the response to fatty acid exposure showed that conditional proto-silencing activity is dependent on Sir2p, a molecule known to be involved in subtelomeric silencing related to aging. This study reveals a path cells can use to coordinate subtelomeric silencing related to aging with cellular environment, and with the activities of other cellular processes.

Subtelomeric chromatin is subject to evolutionarily conserved complex epigenetic regulation and is implicated in numerous aspects of cellular function including formation of heterochromatin, regulation of stress response pathways and control of lifespan. Subtelomeric DNA is characterized by the presence of specific repeated segments that serve to propagate silencing or to protect chromosomal regions from spreading epigenetic control. In this study, analysis of genome-wide chromatin immunoprecipitation and expression data, suggests that several yeast transcription factors regulate subtelomeric silencing in response to various environmental stimuli through conditional association with proto-silencing regions called X elements. In this context, Oaf1p, Rox1p, Gzf1p and Phd1p control the propagation of silencing toward centromeres in response to stimuli affecting stress responses and metabolism, whereas others, including Adr1p, Yap5p and Msn4p, appear to influence boundaries of silencing, regulating telomere-proximal genes in Y′ elements. The factors implicated here are known to control adjacent genes at intrachromosomal positions, suggesting their dual functionality. This study reveals a path for the coordination of subtelomeric silencing with cellular environment, and with activities of other cellular processes.

The phosphoinositide 3-kinase Vps34p is required for pexophagy in Saccharomyces cerevisiae

PIds (phosphoinositides) are phosphorylated derivatives of the membrane phospholipid PtdIns that have emerged as key regulators of many aspects of cellular physiology. We have discovered a PtdIns3P-synthesizing activity in peroxisomes of Saccharomyces cerevisiae and have demonstrated that the lipid kinase Vps34p is already associated with peroxisomes during biogenesis. However, although Vps34 is required, it is not essential for optimal peroxisome biogenesis. The function of Vps34p-containing complex I as well as a subset of PtdIns3P-binding proteins proved to be mandatory for the regulated degradation of peroxisomes. This demonstrates that PtdIns3P-mediated signalling is required for pexophagy.

Comparative analysis of transcriptome and fitness profiles reveals general and condition-specific cellular functions involved in adaptation to environmental change in Saccharomyces cerevisiae

The transcriptional responses of yeast cells to a wide variety of stress conditions have been studied extensively. In addition, deletion mutant collections have been widely used to measure the combined effect of gene loss and stress on growth (fitness). Here we present a comparative analysis of 1,095 publicly available transcription and genome-wide fitness profiles in yeast, from different laboratories and experimental platforms. We analyzed these data, using T-profiler, to describe the correlation in behavior of a priori defined functional groups. Two-mode clustering analysis of the fitness T-profiles revealed that functional groups involved in regulating ribosome biogenesis and translation offer general stress resistance. These groups are closely related to growth rate and nutrient availability. General stress sensitivity was found in deletion mutant groups functioning in intracellular vesicular transport, actin cytoskeleton organization, and cell polarity, indicating that they play an key role in maintaining yeast adaptability. Analysis of the phenotypic and transcriptional variability of our a priori defined functional groups showed that the quantitative effect on fitness of both resistant and sensitive groups is highly condition-dependent. Finally, we discuss the implications of our results for combinatorial drug design.

Genome-wide analysis of effectors of peroxisome biogenesis

Peroxisomes are intracellular organelles that house a number of diverse metabolic processes, notably those required for beta-oxidation of fatty acids. Peroxisomes biogenesis can be induced by the presence of peroxisome proliferators, including fatty acids, which activate complex cellular programs that underlie the induction process. Here, we used multi-parameter quantitative phenotype analyses of an arrayed mutant collection of yeast cells induced to proliferate peroxisomes, to establish a comprehensive inventory of genes required for peroxisome induction and function. The assays employed include growth in the presence of fatty acids, and confocal imaging and flow cytometry through the induction process. In addition to the classical phenotypes associated with loss of peroxisomal functions, these studies identified 169 genes required for robust signaling, transcription, normal peroxisomal development and morphologies, and transmission of peroxisomes to daughter cells. These gene products are localized throughout the cell, and many have indirect connections to peroxisome function. By integration with extant data sets, we present a total of 211 genes linked to peroxisome biogenesis and highlight the complex networks through which information flows during peroxisome biogenesis and function.

Oleate inhibits steryl ester synthesis and causes liposensitivity in yeast

In the yeast Saccharomyces cerevisiae, neutral lipids can be synthesized by four acyltransferases, namely Dga1p and Lro1p producing triacylglycerols (TAG) and Are1p and Are2p forming steryl esters (SE). TAG and SE are stored in an organelle called lipid particles/droplet. Growth of yeast cells on oleate-supplemented media strongly induced proliferation of lipid particles and specifically the synthesis of TAG, which serve as the major pool for the excess of fatty acids. Surprisingly, SE synthesis was strongly inhibited under these conditions. Here, we show that this effect was not due to decreased expression of ARE2 encoding the major yeast SE synthase at the transcriptional level but to competitive enzymatic inhibition of Are2p by free oleate. Consequently, a triple mutant dga1Deltalro1Deltaare1DeltaARE2(+) grown on oleate did not form substantial amounts of SE and exhibited a growth phenotype similar to the dga1Deltalro1Deltaare1Deltaare2Delta quadruple mutant, including lack of lipid particles. Growth of these mutants on oleate was strongly delayed, and cell viability was decreased but rescued by adaptation. In these strains, oleate stress caused morphological changes of intracellular membranes, altered phospholipid composition and formation of an additional lipid class, ethyl esters of fatty acids. In summary, our data showed that exposure to oleate led to disturbed lipid and membrane homeostasis along with liposensitivity of the yeast.

Integrated phosphoproteomics analysis of a signaling network governing nutrient response and peroxisome induction

Phosphorylation of proteins is a key posttranslational modification in cellular signaling, regulating many aspects of cellular responses. We used a quantitative, integrated, phosphoproteomics approach to characterize the cellular responses of the yeast Saccharomyces cerevisiae to the fatty acid oleic acid, a molecule with broad human health implications and a potent inducer of peroxisomes. A combination of cryolysis and urea solubilization was used to minimize the opportunity for reorientation of the phosphoproteome, and hydrophilic interaction liquid chromatography and IMAC chemistries were used to fractionate and enrich for phosphopeptides. Using these approaches, numerous phosphorylated peptides specific to oleate-induced and glucose-repressed conditions were identified and mapped to known signaling pathways. These include several transcription factors, two of which, Pip2p and Cst6p, must be phosphorylated for the normal transcriptional response of fatty acid-responsive loci encoding peroxisomal proteins. The phosphoproteome data were integrated with results from genome-wide assays studying the effects of signaling molecule deletions and known protein-protein interactions to generate a putative fatty acid-responsive signaling network. In this network, the most highly connected nodes are those with the largest effects on cellular responses to oleic acid. These properties are consistent with a scale-free topology, demonstrating that scale-free properties are conserved in condition-specific networks.

Snf1 dependence of peroxisomal gene expression is mediated by Adr1

Eukaryotes utilize fatty acids by beta-oxidation, which occurs in the mitochondria and peroxisomes in higher organisms and in the peroxisomes in yeast. The AMP-activated protein kinase regulates this process in mammalian cells, and its homolog Snf1, together with the transcription factors Adr1, Oaf1, and Pip2, regulates peroxisome proliferation and beta-oxidation in yeast. A constitutive allele of Adr1 (Adr1(c)) lacking the glucose- and Snf1-regulated phosphorylation substrate Ser-230 was found to be Snf1-independent for regulation of peroxisomal genes. In addition, it could compensate for and even suppress the requirement for Oaf1 or Pip2 for gene induction. Peroxisomal genes were found to be regulated by oleate in the presence of glucose, as long as Adr1(c) was expressed, suggesting that the Oaf1/Pip2 heterodimer is Snf1-independent. Consistent with this observation, Oaf1 binding to promoters in the presence of oleate was not reduced in a snf1Delta strain. Exploring the mechanism by which Adr1(c) permits Snf1-independent peroxisomal gene induction, we found that strength of promoter binding did not correlate with transcription, suggesting that stable binding is not a prerequisite for enhanced transcription. Instead, enhanced transcriptional activation and suppression of Oaf1, Pip2, and Snf1 by Adr1(c) may be related to the ability of Adr1(c) to suppress the requirement for and enhance the recruitment of transcriptional coactivators in a promoter- and growth medium-dependent manner.

Dynamic changes in the subcellular distribution of Gpd1p in response to cell stress

Gpd1p is a cytosolic NAD(+)-dependent glycerol 3-phosphate dehydrogenase that also localizes to peroxisomes and plays an essential role in the cellular response to osmotic stress and a role in redox balance. Here, we show that Gpd1p is directed to peroxisomes by virtue of an N-terminal type 2 peroxisomal targeting signal (PTS2) in a Pex7p-dependent manner. Significantly, localization of Gpd1p to peroxisomes is dependent on the metabolic status of cells and the phosphorylation of aminoacyl residues adjacent to the targeting signal. Exposure of cells to osmotic stress induces changes in the subcellular distribution of Gpd1p to the cytosol and nucleus. This behavior is similar to Pnc1p, which is coordinately expressed with Gpd1p, and under conditions of cell stress changes its subcellular distribution from peroxisomes to the nucleus where it mediates chromatin silencing. Although peroxisomes are necessary for the beta-oxidation of fatty acids in yeast, the localization of Gpd1p to peroxisomes is not. Rather, shifts in the distribution of Gpd1p to different cellular compartments in response to changing cellular status suggests a role for Gpd1p in the spatial regulation of redox potential, a process critical to cell survival, especially under the complex stress conditions expected to occur in the wild.

ER stress response, peroxisome proliferation, mitochondrial unfolded protein response and Golgi stress response

The endoplasmic reticulum (ER) response has been thought a cytoprotective mechanism to cope with accumulation of unfolded proteins in the ER. Recent progress has made a quantum leap revealing that ER stress response can be regarded as an autoregulatory system adjusting the ER capacity to cellular demand. This Copernican change raised a novel fundamental question in cell biology: how do cells regulate the capacity of each organelle in accordance with cellular needs? Although this fundamental question has not been fully addressed yet, research about each organelle has been advancing. The proliferation of the peroxisome is regulated by the PPAR alpha pathway, whereas the abundance of mitochondria appears to be regulated by mitochondrial retrograde signaling and the mitochondrial unfolded protein response. The functional capacity of the Golgi apparatus may be regulated by the mechanism of the Golgi stress response. These observations strongly suggest the existence of an elaborate network of organelle autoregulation in eukaryotic cells.

Endoplasmic reticulum-associated secretory proteins Sec20p, Sec39p, and Dsl1p are involved in peroxisome biogenesis

Two pathways have been identified for peroxisome formation: (i) growth and division and (ii) de novo synthesis. Recent experiments determined that peroxisomes originate at the endoplasmic reticulum (ER). Although many proteins have been implicated in the peroxisome biogenic program, no proteins in the eukaryotic secretory pathway have been identified as having roles in peroxisome formation. Using the yeast Saccharomyces cerevisiae regulatable Tet promoter Hughes clone collection, we found that repression of the ER-associated secretory proteins Sec20p and Sec39p resulted in mislocalization of the peroxisomal matrix protein chimera Pot1p-green fluorescent protein (GFP) to the cytosol. Likewise, the peroxisomal membrane protein chimera Pex3p-GFP localized to tubular-vesicular structures in cells suppressed for Sec20p, Sec39p, and Dsl1p, which form a complex at the ER. Loss of Sec39p attenuated formation of Pex3p-derived peroxisomal structures following galactose induction of Pex3p-GFP expression from the GAL1 promoter. Expression of Sec20p, Sec39p, and Dsl1p was moderately increased in yeast grown under conditions that proliferate peroxisomes, and Sec20p, Sec39p, and Dsl1p were found to cofractionate with peroxisomes and colocalize with Pex3p-monomeric red fluorescent protein under these conditions. Our results show that SEC20, SEC39, and DSL1 are essential secretory genes involved in the early stages of peroxisome assembly, and this work is the first to identify and characterize an ER-associated secretory machinery involved in peroxisome biogenesis.

Role of the histone variant H2A.Z/Htz1p in TBP recruitment, chromatin dynamics, and regulated expression of oleate-responsive genes

The histone variant H2A.Z (Htz1p) has been implicated in transcriptional regulation in numerous organisms, including Saccharomyces cerevisiae. Genome-wide transcriptome profiling and chromatin immunoprecipitation studies identified a role for Htz1p in the rapid and robust activation of many oleate-responsive genes encoding peroxisomal proteins, in particular POT1, POX1, FOX2, and CTA1. The Swr1p-, Gcn5p-, and Chz1p-dependent association of Htz1p with these promoters in their repressed states appears to establish an epigenetic marker for the rapid and strong expression of these highly inducible promoters. Isw2p also plays a role in establishing the nucleosome state of these promoters and associates stably in the absence of Htz1p. An analysis of the nucleosome dynamics and Htz1p association with these promoters suggests a complex mechanism in which Htz1p-containing nucleosomes at fatty acid-responsive promoters are disassembled upon initial exposure to oleic acid leading to the loss of Htz1p from the promoter. These nucleosomes reassemble at later stages of gene expression. While these new nucleosomes do not incorporate Htz1p, the initial presence of Htz1p appears to mark the promoter for sustained gene expression and the recruitment of TATA-binding protein.

Bridging high-throughput genetic and transcriptional data reveals cellular responses to alpha-synuclein toxicity

Cells respond to stimuli by changes in various processes, including signaling pathways and gene expression. Efforts to identify components of these responses increasingly depend on mRNA profiling and genetic library screens. By comparing the results of these two assays across various stimuli, we found that genetic screens tend to identify response regulators, whereas mRNA profiling frequently detects metabolic responses. We developed an integrative approach that bridges the gap between these data using known molecular interactions, thus highlighting major response pathways. We used this approach to reveal cellular pathways responding to the toxicity of alpha-synuclein, a protein implicated in several neurodegenerative disorders including Parkinson's disease. For this we screened an established yeast model to identify genes that when overexpressed alter alpha-synuclein toxicity. Bridging these data and data from mRNA profiling provided functional explanations for many of these genes and identified previously unknown relations between alpha-synuclein toxicity and basic cellular pathways.

Spatiotemporal dynamics of the ER-derived peroxisomal endomembrane system

Recent studies have provided evidence that peroxisomes constitute a multicompartmental endomembrane system. The system begins to form with the targeting of certain peroxisomal membrane proteins to the ER and their exit from the ER via preperoxisomal carriers. These carriers undergo a multistep maturation into metabolically active peroxisomes containing the entire complement of peroxisomal membrane and matrix proteins. At each step, the import of a subset of proteins and the uptake of certain membrane lipids result in the formation of a distinct, more mature compartment of the peroxisomal endomembrane system. Individual peroxisomal compartments proliferate by undergoing one or several rounds of division. Herein, we discuss various strategies that evolutionarily diverse organisms use to coordinate compartment formation, maturation, and division in the peroxisomal endomembrane system. We also critically evaluate the molecular and cellular mechanisms governing these processes, outline the most important unanswered questions, and suggest directions for future research.

Cytoscape: a community-based framework for network modeling

Cytoscape is a general network visualization, data integration, and analysis software package. Its development and use has been focused on the modeling requirements of systems biology, though it has been used in other fields. Cytoscape's flexibility has encouraged many users to adopt it and adapt it to their own research by using the plugin framework offered to specialize data analysis, data integration, or visualization. Plugins represent collections of community-contributed functionality and can be used to dynamically extend Cytoscape functionality. This community of users and developers has worked together since Cytoscape's initial release to improve the basic project through contributions to the core code and public offerings of plugin modules. This chapter discusses what Cytoscape does, why it was developed, and the extensions numerous groups have made available to the public. It also describes the development of a plugin used to investigate a particular research question in systems biology and walks through an example analysis using Cytoscape.

Genome-wide analysis of signaling networks regulating fatty acid-induced gene expression and organelle biogenesis

Reversible phosphorylation is the most common posttranslational modification used in the regulation of cellular processes. This study of phosphatases and kinases required for peroxisome biogenesis is the first genome-wide analysis of phosphorylation events controlling organelle biogenesis. We evaluate signaling molecule deletion strains of the yeast Saccharomyces cerevisiae for presence of a green fluorescent protein chimera of peroxisomal thiolase, formation of peroxisomes, and peroxisome functionality. We find that distinct signaling networks involving glucose-mediated gene repression, derepression, oleate-mediated induction, and peroxisome formation promote stages of the biogenesis pathway. Additionally, separate classes of signaling proteins are responsible for the regulation of peroxisome number and size. These signaling networks specify the requirements of early and late events of peroxisome biogenesis. Among the numerous signaling proteins involved, Pho85p is exceptional, with functional involvements in both gene expression and peroxisome formation. Our study represents the first global study of signaling networks regulating the biogenesis of an organelle.

Transcriptional responses to fatty acid are coordinated by combinatorial control

In transcriptional regulatory networks, the coincident binding of a combination of factors to regulate a gene implies the existence of complex mechanisms to control both the gene expression profile and specificity of the response. Unraveling this complexity is a major challenge to biologists. Here, a novel network topology-based clustering approach was applied to condition-specific genome-wide chromatin localization and expression data to characterize a dynamic transcriptional regulatory network responsive to the fatty acid oleate. A network of four (predicted) regulators of the response (Oaf1p, Pip2p, Adr1p and Oaf3p) was investigated. By analyzing trends in the network structure, we found that two groups of multi-input motifs form in response to oleate, each controlling distinct functional classes of genes. This functionality is contributed in part by Oaf1p, which is a component of both types of multi-input motifs and has two different regulatory activities depending on its binding context. The dynamic cooperation between Oaf1p and Pip2p appears to temporally synchronize the two different responses. Together, these data suggest a network mechanism involving dynamic combinatorial control for coordinating transcriptional responses.