Cells balance glycolysis with respiration to support their metabolic needs in different environmental or physiological contexts. With abundant glucose, many cells prefer to grow by aerobic glycolysis or fermentation. Using 161 natural isolates of fission yeast, we investigated the genetic basis and phenotypic effects of the fermentation-respiration balance. The laboratory and a few other strains depended more on respiration. This trait was associated with a single nucleotide polymorphism in a conserved region of Pyk1, the sole pyruvate kinase in fission yeast. This variant reduced Pyk1 activity and glycolytic flux. Replacing the "low-activity" pyk1 allele in the laboratory strain with the "high-activity" allele was sufficient to increase fermentation and decrease respiration. This metabolic rebalancing triggered systems-level adjustments in the transcriptome and proteome and in cellular traits, including increased growth and chronological lifespan but decreased resistance to oxidative stress. Thus, low Pyk1 activity does not lead to a growth advantage but to stress tolerance. The genetic tuning of glycolytic flux may reflect an adaptive trade-off in a species lacking pyruvate kinase isoforms.

Single-cell proteomics aims to characterize biological function and heterogeneity at the level of proteins in an unbiased manner. It is currently limited in proteomic depth, throughput, and robustness, which we address here by a streamlined multiplexed workflow using data-independent acquisition (mDIA). We demonstrate automated and complete dimethyl labeling of bulk or single-cell samples, without losing proteomic depth. Lys-N digestion enables five-plex quantification at MS1 and MS2 level. Because the multiplexed channels are quantitatively isolated from each other, mDIA accommodates a reference channel that does not interfere with the target channels. Our algorithm RefQuant takes advantage of this and confidently quantifies twice as many proteins per single cell compared to our previous work (Brunner et al, PMID 35226415), while our workflow currently allows routine analysis of 80 single cells per day. Finally, we combined mDIA with spatial proteomics to increase the throughput of Deep Visual Proteomics seven-fold for microdissection and four-fold for MS analysis. Applying this to primary cutaneous melanoma, we discovered proteomic signatures of cells within distinct tumor microenvironments, showcasing its potential for precision oncology.

Genomic variation impacts on cellular networks by affecting the abundance (e.g., protein levels) and the functional states (e.g., protein phosphorylation) of their components. Previous work has focused on the former, while in this context, the functional states of proteins have largely remained neglected. Here, we generated high-quality transcriptome, proteome, and phosphoproteome data for a panel of 112 genomically well-defined yeast strains. Genetic effects on transcripts were generally transmitted to the protein layer, but specific gene groups, such as ribosomal proteins, showed diverging effects on protein levels compared with RNA levels. Phosphorylation states proved crucial to unravel genetic effects on signaling networks. Correspondingly, genetic variants that cause phosphorylation changes were mostly different from those causing abundance changes in the respective proteins. Underscoring their relevance for cell physiology, phosphorylation traits were more strongly correlated with cell physiological traits such as chemical compound resistance or cell morphology, compared with transcript or protein abundance. This study demonstrates how molecular networks mediate the effects of genomic variants to cellular traits and highlights the particular importance of protein phosphorylation.

The complexity of many cellular and organismal traits results from the integration of genetic and environmental factors via molecular networks. Network structure and effect propagation are best understood at the level of functional modules, but so far, no concept has been established to include the global network state. Here, we show when and how genetic perturbations lead to molecular changes that are confined to small parts of a network versus when they lead to modulation of network states. Integrating multi-omics profiling of genetically heterogeneous budding and fission yeast strains with an array of cellular traits identified a central state transition of the yeast molecular network that is related to PKA and TOR (PT) signaling. Genetic variants affecting this PT state globally shifted the molecular network along a single-dimensional axis, thereby modulating processes including energy and amino acid metabolism, transcription, translation, cell cycle control, and cellular stress response. We propose that genetic effects can propagate through large parts of molecular networks because of the functional requirement to centrally coordinate the activity of fundamental cellular processes.

Proteins are major effectors and regulators of biological processes that can elicit multiple functions depending on their interaction with other proteins. The organization of proteins into macromolecular complexes and their quantitative distribution across these complexes is, therefore, of great biological and clinical significance. In this paper, we describe an integrated experimental and computational technique to quantify hundreds of protein complexes in a single operation. The method consists of size exclusion chromatography (SEC) to fractionate native protein complexes, SWATH/DIA mass spectrometry to precisely quantify the proteins in each SEC fraction, and the computational framework CCprofiler to detect and quantify protein complexes by error-controlled, complex-centric analysis using prior information from generic protein interaction maps. Our analysis of the HEK293 cell line proteome delineates 462 complexes composed of 2,127 protein subunits. The technique identifies novel sub-complexes and assembly intermediates of central regulatory complexes while assessing the quantitative subunit distribution across them. We make the toolset CCprofiler freely accessible and provide a web platform, SECexplorer, for custom exploration of the HEK293 proteome modularity.

Our current understanding of how natural genetic variation affects gene expression beyond well-annotated coding genes is still limited. The use of deep sequencing technologies for the study of expression quantitative trait loci (eQTLs) has the potential to close this gap. Here, we generated the first recombinant strain library for fission yeast and conducted an RNA-seq-based QTL study of the coding, non-coding, and antisense transcriptomes. We show that the frequency of distal effects (trans-eQTLs) greatly exceeds the number of local effects (cis-eQTLs) and that non-coding RNAs are as likely to be affected by eQTLs as protein-coding RNAs. We identified a genetic variation of swc5 that modifies the levels of 871 RNAs, with effects on both sense and antisense transcription, and show that this effect most likely goes through a compromised deposition of the histone variant H2A.Z. The strains, methods, and datasets generated here provide a rich resource for future studies.

Human erythropoiesis is an exquisitely controlled multistep developmental process, and its dysregulation leads to numerous human diseases. Transcriptome and epigenome studies provided insights into system-wide regulation, but we currently lack a global mechanistic view on the dynamics of proteome and post-translational regulation coordinating erythroid maturation. We established a mass spectrometry (MS)-based proteomics workflow to quantify and dynamically track 7,400 proteins and 27,000 phosphorylation sites of five distinct maturation stages of in vitro reconstituted erythropoiesis of CD34+ HSPCs. Our data reveal developmental regulation through drastic proteome remodeling across stages of erythroid maturation encompassing most protein classes. This includes various orchestrated changes in solute carriers indicating adjustments to altered metabolic requirements. To define the distinct proteome of each maturation stage, we developed a computational deconvolution approach which revealed stage-specific marker proteins. The dynamic phosphoproteomes combined with a kinome-targeted CRISPR/Cas9 screen uncovered coordinated networks of erythropoietic kinases and pinpointed downregulation of c-Kit/MAPK signaling axis as key driver of maturation. Our system-wide view establishes the functional dynamic of complex phosphosignaling networks and regulation through proteome remodeling in erythropoiesis.