PTEN Regulates PI(3,4)P2 Signaling Downstream of Class I PI3K.

The PI3K signaling pathway regulates cell growth and movement and is heavily mutated in cancer. Class I PI3Ks synthesize the lipid messenger PI(3,4,5)P3. PI(3,4,5)P3 can be dephosphorylated by 3- or 5-phosphatases, the latter producing PI(3,4)P2. The PTEN tumor suppressor is thought to function primarily as a PI(3,4,5)P3 3-phosphatase, limiting activation of this pathway. Here we show that PTEN also functions as a PI(3,4)P2 3-phosphatase, both in vitro and in vivo. PTEN is a major PI(3,4)P2 phosphatase in Mcf10a cytosol, and loss of PTEN and INPP4B, a known PI(3,4)P2 4-phosphatase, leads to synergistic accumulation of PI(3,4)P2, which correlated with increased invadopodia in epidermal growth factor (EGF)-stimulated cells. PTEN deletion increased PI(3,4)P2 levels in a mouse model of prostate cancer, and it inversely correlated with PI(3,4)P2 levels across several EGF-stimulated prostate and breast cancer lines. These results point to a role for PI(3,4)P2 in the phenotype caused by loss-of-function mutations or deletions in PTEN.

MODOMICS: a database of RNA modification pathways. 2021 update.

The MODOMICS database has been, since 2006, a manually curated and centralized resource, storing and distributing comprehensive information about modified ribonucleosides. Originally, it only contained data on the chemical structures of modified ribonucleosides, their biosynthetic pathways, the location of modified residues in RNA sequences, and RNA-modifying enzymes. Over the years, prompted by the accumulation of new knowledge and new types of data, it has been updated with new information and functionalities. In this new release, we have created a catalog of RNA modifications linked to human diseases, e.g., due to mutations in genes encoding modification enzymes. MODOMICS has been linked extensively to RCSB Protein Data Bank, and sequences of experimentally determined RNA structures with modified residues have been added. This expansion was accompanied by including nucleotide 5'-monophosphate residues. We redesigned the web interface and upgraded the database backend. In addition, a search engine for chemically similar modified residues has been included that can be queried by SMILES codes or by drawing chemical molecules. Finally, previously available datasets of modified residues, biosynthetic pathways, and RNA-modifying enzymes have been updated. Overall, we provide users with a new, enhanced, and restyled tool for research on RNA modification. MODOMICS is available at https://iimcb.genesilico.pl/modomics/.

A mark of disease: how mRNA modifications shape genetic and acquired pathologies.

RNA modifications have recently emerged as a widespread and complex facet of gene expression regulation. Counting more than 170 distinct chemical modifications with far-reaching implications for RNA fate, they are collectively referred to as the epitranscriptome. These modifications can occur in all RNA species, including messenger RNAs (mRNAs) and noncoding RNAs (ncRNAs). In mRNAs the deposition, removal, and recognition of chemical marks by writers, erasers and readers influence their structure, localization, stability, and translation. In turn, this modulates key molecular and cellular processes such as RNA metabolism, cell cycle, apoptosis, and others. Unsurprisingly, given their relevance for cellular and organismal functions, alterations of epitranscriptomic marks have been observed in a broad range of human diseases, including cancer, neurological and metabolic disorders. Here, we will review the major types of mRNA modifications and editing processes in conjunction with the enzymes involved in their metabolism and describe their impact on human diseases. We present the current knowledge in an updated catalog. We will also discuss the emerging evidence on the crosstalk of epitranscriptomic marks and what this interplay could imply for the dynamics of mRNA modifications. Understanding how this complex regulatory layer can affect the course of human pathologies will ultimately lead to its exploitation toward novel epitranscriptomic therapeutic strategies.

Positioning Europe for the EPITRANSCRIPTOMICS challenge.

The genetic alphabet consists of the four letters: C, A, G, and T in DNA and C,A,G, and U in RNA. Triplets of these four letters jointly encode 20 different amino acids out of which proteins of all organisms are built. This system is universal and is found in all kingdoms of life. However, bases in DNA and RNA can be chemically modified. In DNA, around 10 different modifications are known, and those have been studied intensively over the past 20 years. Scientific studies on DNA modifications and proteins that recognize them gave rise to the large field of epigenetic and epigenomic research. The outcome of this intense research field is the discovery that development, ageing, and stem-cell dependent regeneration but also several diseases including cancer are largely controlled by the epigenetic state of cells. Consequently, this research has already led to the first FDA approved drugs that exploit the gained knowledge to combat disease. In recent years, the ~150 modifications found in RNA have come to the focus of intense research. Here we provide a perspective on necessary and expected developments in the fast expanding area of RNA modifications, termed epitranscriptomics.

Optimized network based natural language processing approach to reveal disease comorbidities in COVID-19.

A novel virus emerged from Wuhan, China, at the end of 2019 and quickly evolved into a pandemic, significantly impacting various industries, especially healthcare. One critical lesson from COVID-19 is the importance of understanding and predicting underlying comorbidities to better prioritize care and pharmacological therapies. Factors like age, race, and comorbidity history are crucial in determining disease mortality. While clinical data from hospitals and cohorts have led to the identification of these comorbidities, traditional approaches often lack a mechanistic understanding of the connections between them. In response, we utilized a deep learning approach to integrate COVID-19 data with data from other diseases, aiming to detect comorbidities with mechanistic insights. Our modified algorithm in the mpDisNet package, based on word-embedding deep learning techniques, incorporates miRNA expression profiles from SARS-CoV-2 infected cell lines and their target transcription factors. This approach is aligned with the emerging field of network medicine, which seeks to define diseases based on distinct pathomechanisms rather than just phenotypes. The main aim is discovery of possible unknown comorbidities by connecting the diseases by their miRNA mediated regulatory interactions. The algorithm can predict the majority of COVID-19's known comorbidities, as well as several diseases that have yet to be discovered to be comorbid with COVID-19. These potentially comorbid diseases should be investigated further to raise awareness and prevention, as well as informing the comorbidity research for the next possible outbreak.

Control analysis of the eukaryotic cell cycle using gene copy-number series in yeast tetraploids.

BACKGROUND: In the model eukaryote, Saccharomyces cerevisiae, previous experiments have identified those genes that exert the most significant control over cell growth rate. These genes are termed HFC for high flux control. Such genes are overrepresented within pathways controlling the mitotic cell cycle. RESULTS: We postulated that the increase/decrease in growth rate is due to a change in the rate of progression through specific cell cycle steps. We extended and further developed an existing logical model of the yeast cell cycle in order elucidate how the HFC genes modulated progress through the cycle. This model can simulate gene dosage-variation and calculate the cycle time, determine the order and relative speed at which events occur, and predict arrests and failures to correctly execute a step. To experimentally test our model's predictions, we constructed a tetraploid series of deletion mutants for a set of eight genes that control the G2/M transition. This system allowed us to vary gene copy number through more intermediate levels than previous studies and examine the impact of copy-number variation on growth, cell-cycle phenotype, and response to different cellular stresses. CONCLUSIONS: For the majority of strains, the predictions agreed with experimental observations, validating our model and its use for further predictions. Where simulation and experiment diverged, we uncovered both novel tetraploid-specific phenotypes and a switch in the determinative execution point of a key cell-cycle regulator, the Cdc28 kinase, from the G1/S to the S/G2 boundaries.

Absolute quantification of the glycolytic pathway in yeast: deployment of a complete QconCAT approach.

The availability of label-free data derived from yeast cells (based on the summed intensity of the three strongest, isoform-specific peptides) permitted a preliminary assessment of protein abundances for glycolytic proteins. Following this analysis, we demonstrate successful application of the QconCAT technology, which uses recombinant DNA techniques to generate artificial concatamers of large numbers of internal standard peptides, to the quantification of enzymes of the glycolysis pathway in the yeast Saccharomyces cerevisiae. A QconCAT of 88 kDa (59 tryptic peptides) corresponding to 27 isoenzymes was designed and built to encode two or three analyte peptides per protein, and after stable isotope labeling of the standard in vivo, protein levels were determined by LC-MS, using ultra high performance liquid chromatography-coupled mass spectrometry. We were able to determine absolute protein concentrations between 14,000 and 10 million molecules/cell. Issues such as efficiency of extraction and completeness of proteolysis are addressed, as well as generic factors such as optimal quantotypic peptide selection and expression. In addition, the same proteins were quantified by intensity-based label-free analysis, and both sets of data were compared with other quantification methods.

A comparative performance evaluation of imputation methods in spatially resolved transcriptomics data.

Spatially resolved transcriptomics technologies have drawn enormous attention by providing RNA expression patterns together with their spatial information. Even though improved techniques are being developed rapidly, the technologies which give spatially whole transcriptome level profiles suffer from dropout problems because of the low capture rate. Imputation of missing data is one strategy to eliminate this technical problem. We evaluated the imputation performance of five available methods (SpaGE, stPlus, gimVI, Tangram and stLearn) which were indicated as capable of making predictions for the dropouts in spatially resolved transcriptomics datasets. The evaluation was performed qualitatively via visualization of the predictions against the original values and quantitatively with Pearson's correlation coefficient, cosine similarity, root mean squared log-error, Silhouette Index and Calinski Harabasz Index. We found that stPlus and gimVI outperform the other three. However, the performance of all methods was lower than expected which indicates that there is still a gap for imputation tools dealing with dropout events in spatially resolved transcriptomics.

An integrated study to decipher immunosuppressive cellular communication in the PDAC environment.

Pancreatic ductal adenocarcinoma (PDAC) is one the most aggressive cancers and characterized by a highly rigid and immunosuppressive tumor microenvironment (TME). The extensive cellular interactions are known to play key roles in the immune evasion, chemoresistance, and poor prognosis. Here, we used the spatial transcriptomics, scRNA-seq, and bulk RNA-seq datasets to enhance the insights obtained from each to decipher the cellular communication in the TME. The complex crosstalk in PDAC samples was revealed by the single-cell and spatial transcriptomics profiles of the samples. We show that tumor-associated macrophages (TAMs) are the central cell types in the regulation of microenvironment in PDAC. They colocalize with the cancer cells and tumor-suppressor immune cells and take roles to provide an immunosuppressive environment. LGALS9 gene which is upregulated in PDAC tumor samples in comparison to healthy samples was also found to be upregulated in TAMs compared to tumor-suppressor immune cells in cancer samples. Additionally, LGALS9 was found to be the primary component in the crosstalk between TAMs and the other cells. The widespread expression of P4HB gene and its interaction with LGALS9 was also notable. Our findings point to a profound role of TAMs via LGALS9 and its interaction with P4HB that should be considered for further elucidation as target in the combinatory immunotherapies for PDAC.

SUMA: a lightweight machine learning model-powered shared nearest neighbour-based clustering application interface for scRNA-Seq data.

Background/aim: Single-cell transcriptomics (scRNA-Seq) explores cellular diversity at the gene expression level. Due to the inherent sparsity and noise in scRNA-Seq data and the uncertainty on the types of sequenced cells, effective clustering and cell type annotation are essential. The graph-based clustering of scRNA-Seq data is a simple yet powerful approach that presents data as a "shared nearest neighbour" graph and clusters the cells using graph clustering algorithms. These algorithms are dependent on several user-defined parameters.Here we present SUMA, a lightweight tool that uses a random forest model to predict the optimum number of neighbours to obtain the optimum clustering results. Moreover, we integrated our method with other commonly used methods in an RShiny application. SUMA can be used in a local environment (https://github.com/hkarakurt8742/SUMA) or as a browser tool (https://hkarakurt.shinyapps.io/suma/). Materials and methods: Publicly available scRNA-Seq datasets and 3 different graph-based clustering algorithms were used to develop SUMA, and a large range for number of neighbours and variant genes was taken into consideration. The quality of clustering was assessed using the adjusted Rand index (ARI) and true labels of each dataset. The data were split into training and test datasets, and the model was built and optimised using Scikit-learn (Python) and randomForest (R) libraries. Results: The accuracy of our machine learning model was 0.96, while the AUC of the ROC curve was 0.98. The model indicated that the number of cells in scRNA-Seq data is the most important feature when deciding the number of neighbours. Conclusion: We developed and evaluated the SUMA model and implemented the method in the SUMAShiny app, which integrates SUMA with different clustering methods and enables nonbioinformatician users to cluster and visualise their scRNA data easily. The SUMAShiny app is available both for desktop and browser use.

Haploinsufficiency and the sex chromosomes from yeasts to humans.

BACKGROUND: Haploinsufficient (HI) genes are those for which a reduction in copy number in a diploid from two to one results in significantly reduced fitness. Haploinsufficiency is increasingly implicated in human disease, and so predicting this phenotype could provide insights into the genetic mechanisms behind many human diseases, including some cancers. RESULTS: In the present work we show that orthologues of Saccharomyces cerevisiae HI genes are preferentially retained across the kingdom Fungi, and that the HI genes of S. cerevisiae can be used to predict haploinsufficiency in humans. Our HI gene predictions confirm known associations between haploinsufficiency and genetic disease, and predict several further disorders in which the phenotype may be relevant. Haploinsufficiency is also clearly relevant to the gene-dosage imbalances inherent in eukaryotic sex-determination systems. In S. cerevisiae, HI genes are over-represented on chromosome III, the chromosome that determines yeast's mating type. This may be a device to select against the loss of one copy of chromosome III from a diploid. We found that orthologues of S. cerevisiae HI genes are also over-represented on the mating-type chromosomes of other yeasts and filamentous fungi. In animals with heterogametic sex determination, accumulation of HI genes on the sex chromosomes would compromise fitness in both sexes, given X chromosome inactivation in females. We found that orthologues of S. cerevisiae HI genes are significantly under-represented on the X chromosomes of mammals and of Caenorhabditis elegans. There is no X inactivation in Drosophila melanogaster (increased expression of X in the male is used instead) and, in this species, we found no depletion of orthologues to yeast HI genes on the sex chromosomes. CONCLUSION: A special relationship between HI genes and the sex/mating-type chromosome extends from S. cerevisiae to Homo sapiens, with the microbe being a useful model for species throughout the evolutionary range. Furthermore, haploinsufficiency in yeast can predict the phenotype in higher organisms.

Genome-wide assessment of the carriers involved in the cellular uptake of drugs: a model system in yeast.

BACKGROUND: The uptake of drugs into cells has traditionally been considered to be predominantly via passive diffusion through the bilayer portion of the cell membrane. The recent recognition that drug uptake is mostly carrier-mediated raises the question of which drugs use which carriers. RESULTS: To answer this, we have constructed a chemical genomics platform built upon the yeast gene deletion collection, using competition experiments in batch fermenters and robotic automation of cytotoxicity screens, including protection by 'natural' substrates. Using these, we tested 26 different drugs and identified the carriers required for 18 of the drugs to gain entry into yeast cells. CONCLUSIONS: As well as providing a useful platform technology, these results further substantiate the notion that the cellular uptake of pharmaceutical drugs normally occurs via carrier-mediated transport and indicates that establishing the identity and tissue distribution of such carriers should be a major consideration in the design of safe and effective drugs.

In silico analysis of metabolic effects of bipolar disorder on prefrontal cortex identified altered GABA, glutamate-glutamine cycle, energy metabolism and amino acid synthesis pathways.

Bipolar disorder (BP) is a lifelong psychiatric condition, which often disrupts the daily life of the patients. It is characterized by unstable and periodic mood changes, which cause patients to display unusual shifts in mood, energy, activity levels, concentration and the ability to carry out day-to-day tasks. BP is a major psychiatric condition, and it is still undertreated. The causes and neural mechanisms of bipolar disorder are unclear, and diagnosis is still mostly based on psychiatric examination, furthermore the unstable character of the disorder makes diagnosis challenging. Identification of the molecular mechanisms underlying the disease may improve the diagnosis and treatment rates. Single nucleotide polymorphisms (SNP) and transcriptome profiles of patients were studied along with signalling pathways that are thought to be associated with bipolar disorder. Here, we present a computational approach that uses publicly available transcriptome data from bipolar disorder patients and healthy controls. Along with statistical analyses, data are integrated with a genome-scale metabolic model and protein-protein interaction network. Healthy individuals and bipolar disorder patients are compared based on their metabolic profiles. We hypothesize that energy metabolism alterations in bipolar disorder relate to perturbations in amino-acid metabolism and neuron-astrocyte exchange reactions. Due to changes in amino acid metabolism, neurotransmitters and their secretion from neurons and metabolic exchange pathways between neurons and astrocytes such as the glutamine-glutamate cycle are also altered. Changes in negatively charged (-1) KIV and KMV molecules are also observed, and it indicates that charge balance in the brain is highly altered in bipolar disorder. Due to this fact, we also hypothesize that positively charged lithium ions may stabilize the disturbed charge balance in neurons in addition to its effects on neurotransmission. To the best of our knowledge, our approach is unique as it is the first study using genome-scale metabolic models in neuropsychiatric research.

Genome-Wide Analysis of Yeast Metabolic Cycle through Metabolic Network Models Reveals Superiority of Integrated ATAC-seq Data over RNA-seq Data.

Saccharomyces cerevisiae undergoes robust oscillations to regulate its physiology for adaptation and survival under nutrient-limited conditions. Environmental cues can induce rhythmic metabolic alterations in order to facilitate the coordination of dynamic metabolic behaviors. Of such metabolic processes, the yeast metabolic cycle enables adaptation of the cells to varying nutritional status through oscillations in gene expression and metabolite production levels. In this process, yeast metabolism is altered between diverse cellular states based on changing oxygen consumption levels: quiescent (reductive charging [RC]), growth (oxidative [OX]), and proliferation (reductive building [RB]) phases. We characterized metabolic alterations during the yeast metabolic cycle using a variety of approaches. Gene expression levels are widely used for condition-specific metabolic simulations, whereas the use of epigenetic information in metabolic modeling is still limited despite the clear relationship between epigenetics and metabolism. This prompted us to investigate the contribution of epigenomic information to metabolic predictions for progression of the yeast metabolic cycle. In this regard, we determined altered pathways through the prediction of regulated reactions and corresponding model genes relying on differential chromatin accessibility levels. The predicted metabolic alterations were confirmed via data analysis and literature. We subsequently utilized RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin using sequencing (ATAC-seq) data sets in the contextualization of the yeast model. The use of ATAC-seq data considerably enhanced the predictive capability of the model. To the best of our knowledge, this is the first attempt to use genome-wide chromatin accessibility data in metabolic modeling. The preliminary results showed that epigenomic data sets can pave the way for more accurate metabolic simulations. IMPORTANCE Dynamic chromatin organization mediates the emergence of condition-specific phenotypes in eukaryotic organisms. Saccharomyces cerevisiae can alter its metabolic profile via regulation of genome accessibility and robust transcriptional oscillations under nutrient-limited conditions. Thus, both epigenetic information and transcriptomic information are crucial in the understanding of condition-specific metabolic behavior in this organism. Based on genome-wide alterations in chromatin accessibility and transcription, we investigated the yeast metabolic cycle, which is a remarkable example of coordinated and dynamic yeast behavior. In this regard, we assessed the use of ATAC-seq and RNA-seq data sets in condition-specific metabolic modeling. To our knowledge, this is the first attempt to use chromatin accessibility data in the reconstruction of context-specific metabolic models, despite the extensive use of transcriptomic data. As a result of comparative analyses, we propose that the incorporation of epigenetic information is a promising approach in the accurate prediction of metabolic dynamics.

Nutrient control of eukaryote cell growth: a systems biology study in yeast.

BACKGROUND: To elucidate the biological processes affected by changes in growth rate and nutrient availability, we have performed a comprehensive analysis of the transcriptome, proteome and metabolome responses of chemostat cultures of the yeast, Saccharomyces cerevisiae, growing at a range of growth rates and in four different nutrient-limiting conditions. RESULTS: We find significant changes in expression for many genes in each of the four nutrient-limited conditions tested. We also observe several processes that respond differently to changes in growth rate and are specific to each nutrient-limiting condition. These include carbohydrate storage, mitochondrial function, ribosome synthesis, and phosphate transport. Integrating transcriptome data with proteome measurements allows us to identify previously unrecognized examples of post-transcriptional regulation in response to both nutrient and growth-rate signals. CONCLUSIONS: Our results emphasize the unique properties of carbon metabolism and the carbon substrate, the limitation of which induces significant changes in gene regulation at the transcriptional and post-transcriptional level, as well as altering how many genes respond to growth rate. By comparison, the responses to growth limitation by other nutrients involve a smaller set of genes that participate in specific pathways. See associated commentary http://www.biomedcentral.com/1741-7007/8/62.

Transcriptional landscape of cellular networks reveal interactions driving the dormancy mechanisms in cancer.

Primary cancer cells exert unique capacity to disseminate and nestle in distant organs. Once seeded in secondary sites, cancer cells may enter a dormant state, becoming resistant to current treatment approaches, and they remain silent until they reactivate and cause overt metastases. To illuminate the complex mechanisms of cancer dormancy, 10 transcriptomic datasets from the literature enabling 21 dormancy-cancer comparisons were mapped on protein-protein interaction networks and gene-regulatory networks to extract subnetworks that are enriched in significantly deregulated genes. The genes appearing in the subnetworks and significantly upregulated in dormancy with respect to proliferative state were scored and filtered across all comparisons, leading to a dormancy-interaction network for the first time in the literature, which includes 139 genes and 1974 interactions. The dormancy interaction network will contribute to the elucidation of cellular mechanisms orchestrating cancer dormancy, paving the way for improvements in the diagnosis and treatment of metastatic cancer.

Integration of transcriptomic profile of SARS-CoV-2 infected normal human bronchial epithelial cells with metabolic and protein-protein interaction networks.

A novel coronavirus (SARS-CoV-2, formerly known as nCoV-2019) that causes an acute respiratory disease has emerged in Wuhan, China and spread globally in early 2020. On January the 30th, the World Health Organization (WHO) declared spread of this virus as an epidemic and a public health emergency. With its highly contagious characteristic and long incubation time, confinement of SARS-CoV-2 requires drastic lock-down measures to be taken and therefore early diagnosis is crucial. We analysed transcriptome of SARS-CoV-2 infected human lung epithelial cells, compared it with mock-infected cells, used network-based reporter metabolite approach and integrated the transcriptome data with protein-protein interaction network to elucidate the early cellular response. Significantly affected metabolites have the potential to be used in diagnostics while pathways of protein clusters have the potential to be used as targets for supportive or novel therapeutic approaches. Our results are in accordance with the literature on response of IL6 family of cytokines and their importance, in addition, we find that matrix metalloproteinase 2 (MMP2) and matrix metalloproteinase 9 (MMP9) with keratan sulfate synthesis pathway may play a key role in the infection. We hypothesize that MMP9 inhibitors have potential to prevent "cytokine storm" in severely affected patients.

Comparison of visualization tools for single-cell RNAseq data.

In the last decade, single cell RNAseq (scRNAseq) datasets have grown in size from a single cell to millions of cells. Due to its high dimensionality, it is not always feasible to visualize scRNAseq data and share it in a scientific report or an article publication format. Recently, many interactive analysis and visualization tools have been developed to address this issue and facilitate knowledge transfer in the scientific community. In this study, we review several of the currently available scRNAseq visualization tools and benchmark the subset that allows to visualize the data on the web and share it with others. We consider the memory and time required to prepare datasets for sharing as the number of cells increases, and additionally review the user experience and features available in the web interface. To address the problem of format compatibility we have also developed a user-friendly R package, sceasy, which allows users to convert their own scRNAseq datasets into a specific data format for visualization.

Next-Generation Genome-Scale Models Incorporating Multilevel 'Omics Data: From Yeast to Human.

Genome-scale modelling in eukaryotes has been pioneered by the yeast Saccharomyces cerevisiae. Early metabolic networks have been reconstructed based on genome sequence and information accumulated in the literature on biochemical reactions. Protein-protein interaction networks have been constructed based on experimental observations such as yeast-2-hybrid method. Gene regulatory networks were based on a variety of data types, including information on TF-promoter binding and gene coexpression. The aforementioned networks have been improved gradually, and methods for their integration were developed. Incorporation of omics data including genomics, metabolomics, transcriptomics, fluxome, and phosphoproteome led to next-generation genome-scale models. The methods tested on yeast have later been implemented in human, further, cellular components found to be important in yeast physiology under (ab)normal conditions, and (dis)regulation mechanisms in yeast shed light to the healthy and disease states in human. This chapter provides a historical perspective on next-generation genome-scale models incorporating multilevel 'omics data, from yeast to human.

Exometabolic and transcriptional response in relation to phenotype and gene copy number in respiration-related deletion mutants of S. cerevisiae.

The transcriptional and metabolic impact of deleting one or both copies of a respiration-related gene has been studied in glucose-limited chemostats. Integration of literature information on phenotype with our exometabolome and transcriptome data enabled the identification of novel relationships between gene copy number, transcriptional regulation and phenotype. We found that the effect of complete respiratory deficiency on transcription was limited to downregulation of genes involved in oxidoreductase activity and iron assimilation. Partial respiratory deficiency had no significant impact on gene transcription. Changes in the copy number of two transcription-factor genes, HAP4 and MIG1, had a major impact on genes involved in mitochondrial function. Regulation of respiratory chain components encoded in the nucleus and mitochondria appears to be divided between Hap4p and Oxa1p, respectively. Similarly, repression of respiration may be imposed by the action of Mig1p and Mba1p on nuclear and mitochondrial gene expression, respectively. However, it is not clear whether Oxa1p and Mba1p regulate mitochondrial gene expression via their interaction with mitochondrial ribosomes or by some indirect means. The phenotype of nuclear petite mutants may not simply be due to the absence of respiration; e.g. Oxa1p or Bcs1p may play a role in the regulation of ribosome assembly in the nucleolus. Integration between respiration and cell growth may also result from the action of a single transcription factor. Thus, Hap4p targets genes that are required for respiration and for fitness in nutrient-limited conditions. This suggests that Hap4p may enable cells to adapt to nutrient limitation as well as diauxy.