The STARD9/Kif16a kinesin associates with mitotic microtubules and regulates spindle pole assembly.

During cell division, cells form the microtubule-based mitotic spindle, a highly specialized and dynamic structure that mediates proper chromosome transmission to daughter cells. Cancer cells can show perturbed mitotic spindles and an approach in cancer treatment has been to trigger cell killing by targeting microtubule dynamics or spindle assembly. To identify and characterize proteins necessary for spindle assembly, and potential antimitotic targets, we performed a proteomic and genetic analysis of 592 mitotic microtubule copurifying proteins (MMCPs). Screening for regulators that affect both mitosis and apoptosis, we report the identification and characterization of STARD9, a kinesin-3 family member, which localizes to centrosomes and stabilizes the pericentriolar material (PCM). STARD9-depleted cells have fragmented PCM, form multipolar spindles, activate the spindle assembly checkpoint (SAC), arrest in mitosis, and undergo apoptosis. Interestingly, STARD9-depletion synergizes with the chemotherapeutic agent taxol to increase mitotic death, demonstrating that STARD9 is a mitotic kinesin and a potential antimitotic target.

Mapping the NPHP-JBTS-MKS protein network reveals ciliopathy disease genes and pathways.

Nephronophthisis (NPHP), Joubert (JBTS), and Meckel-Gruber (MKS) syndromes are autosomal-recessive ciliopathies presenting with cystic kidneys, retinal degeneration, and cerebellar/neural tube malformation. Whether defects in kidney, retinal, or neural disease primarily involve ciliary, Hedgehog, or cell polarity pathways remains unclear. Using high-confidence proteomics, we identified 850 interactors copurifying with nine NPHP/JBTS/MKS proteins and discovered three connected modules: "NPHP1-4-8" functioning at the apical surface, "NPHP5-6" at centrosomes, and "MKS" linked to Hedgehog signaling. Assays for ciliogenesis and epithelial morphogenesis in 3D renal cultures link renal cystic disease to apical organization defects, whereas ciliary and Hedgehog pathway defects lead to retinal or neural deficits. Using 38 interactors as candidates, linkage and sequencing analysis of 250 patients identified ATXN10 and TCTN2 as new NPHP-JBTS genes, and our Tctn2 mouse knockout shows neural tube and Hedgehog signaling defects. Our study further illustrates the power of linking proteomic networks and human genetics to uncover critical disease pathways.

In silico tools for accurate HLA and KIR inference from clinical sequencing data empower immunogenetics on individual-patient and population scales.

Immunogenetic variation in humans is important in research, clinical diagnosis and increasingly a target for therapeutic intervention. Two highly polymorphic loci play critical roles, namely the human leukocyte antigen (HLA) system, which is the human version of the major histocompatibility complex (MHC), and the Killer-cell immunoglobulin-like receptors (KIR) that are relevant for responses of natural killer (NK) and some subsets of T cells. Their accurate classification has typically required the use of dedicated biological specimens and a combination of in vitro and in silico efforts. Increased availability of next generation sequencing data has led to the development of ancillary computational solutions. Here, we report an evaluation of recently published algorithms to computationally infer complex immunogenetic variation in the form of HLA alleles and KIR haplotypes from whole-genome or whole-exome sequencing data. For both HLA allele and KIR gene typing, we identified tools that yielded >97% overall accuracy for four-digit HLA types, and >99% overall accuracy for KIR gene presence, suggesting the readiness of in silico solutions for use in clinical and high-throughput research settings.

Previously reported placebo-response-associated variants do not predict patient outcomes in inflammatory disease Phase III trial placebo arms.

In clinical trials, a placebo response refers to improvement in disease symptoms arising from the psychological effect of receiving a treatment rather than the actual treatment under investigation. Previous research has reported genomic variation associated with the likelihood of observing a placebo response, but these studies have been limited in scope and have not been validated. Here, we analyzed whole-genome sequencing data from 784 patients undergoing placebo treatment in Phase III Asthma or Rheumatoid Arthritis trials to assess the impact of previously reported variation on patient outcomes in the placebo arms and to identify novel variants associated with the placebo response. Contrary to expectations based on previous reports, we did not observe any statistically significant associations between genomic variants and placebo treatment outcome. Our findings suggest that the biological origin of the placebo response is complex and likely to be variable between disease areas.

Identification of mutant K-Ras-dependent phenotypes using a panel of isogenic cell lines.

To assess the consequences of endogenous mutant K-Ras, we analyzed the signaling and biological properties of a small panel of isogenic cell lines. These include the cancer cell lines DLD1, HCT116, and Hec1A, in which either the WT or mutant K-ras allele has been disrupted, and SW48 colorectal cancer cells and human mammary epithelial cells in which a single copy of mutant K-ras was introduced at its endogenous genomic locus. We find that single copy mutant K-Ras causes surprisingly modest activation of downstream signaling to ERK and Akt. In contrast, a negative feedback signaling loop to EGFR and N-Ras occurs in some, but not all, of these cell lines. Mutant K-Ras also had relatively minor effects on cell proliferation and cell migration but more dramatic effects on cell transformation as assessed by growth in soft agar. Surprisingly, knock-out of the wild type K-ras allele consistently increased growth in soft agar, suggesting tumor-suppressive properties of this gene under these conditions. Finally, we examined the effects of single copy mutant K-Ras on global gene expression. Although transcriptional programs triggered by mutant K-Ras were generally quite distinct in the different cell lines, there was a small number of genes that were consistently overexpressed, and these could be used to monitor K-Ras inhibition in a panel of human tumor cell lines. We conclude that there are conserved components of mutant K-Ras signaling and phenotypes but that many depend on cell context and environmental cues.

Development of a robust flow cytometry-based pharmacodynamic assay to detect phospho-protein signals for phosphatidylinositol 3-kinase inhibitors in multiple myeloma.

BACKGROUND: The phosphatidylinositol 3-kinase (PI3K) pathway plays an important role in multiple myeloma (MM), a blood cancer associated with uncontrolled proliferation of bone marrow plasma cells. This study aimed to develop a robust clinical pharmacodynamic (PD) assay to measure the on-target PD effects of the selective PI3K inhibitor GDC-0941 in MM patients. METHODS: We conducted an in vitro drug wash-out study to evaluate the feasibility of biochemical approaches in measuring the phosphorylation of S6 ribosomal protein (S6), one of the commonly used PD markers for PI3K pathway inhibition. We then developed a 7-color phospho-specific flow cytometry assay, or phospho flow assay, to measure the phosphorylation state of intracellular S6 in bone marrow aspirate (BMA) and peripheral blood (PB). Integrated mean fluorescence intensity (iMFI) was used to calculate fold changes of phosphorylation. Assay sensitivity was evaluated by comparing phospho flow with Meso Scale Discovery (MSD) and immunohistochemistry (IHC) assays. Finally, a sample handling method was developed to maintain the integrity of phospho signal during sample shipping and storage to ensure clinical application. RESULTS: The phospho flow assay provided single-cell PD monitoring of S6 phosphorylation in tumor and surrogate cells using fixed BMA and PB, assessing pathway modulation in response to GDC-0941 with sensitivity similar to that of MSD assay. The one-shot sample fixation and handling protocol herein demonstrated exceptional preservation of protein phosphorylation. In contrast, the IHC assay was less sensitive in terms of signal quantification while the biochemical approach (MSD) was less suitable to assess PD activities due to the undesirable impact associated with cell isolation on the protein phosphorylation in tumor cells. CONCLUSIONS: We developed a robust PD biomarker assay for the clinical evaluation of PI3K inhibitors in MM, allowing one to decipher the PD response in a relevant cell population. To our knowledge, this is the first report of an easily implemented clinical PD assay that incorporates an unbiased one-shot sample handling protocol, all (staining)-in-one (tube) phospho flow staining protocol, and an integrated modified data analysis for PD monitoring of kinase inhibitors in relevant cell populations in BMA and PB. The methods described here ensure a real-time, reliable and reproducible PD readout, which can provide information for dose selection as well as help to identify optimal combinations of targeted agents in early clinical trials.

A whole genome sequencing study of moderate to severe asthma identifies a lung function locus associated with asthma risk.

Genome-wide association studies (GWAS) have identified many common variant loci associated with asthma susceptibility, but few studies investigate the genetics underlying moderate-to-severe asthma risk. Here, we present a whole-genome sequencing study comparing 3181 moderate-to-severe asthma patients to 3590 non-asthma controls. We demonstrate that asthma risk is genetically correlated with lung function measures and that this component of asthma risk is orthogonal to the eosinophil genetics that also contribute to disease susceptibility. We find that polygenic scores for reduced lung function are associated with younger asthma age of onset. Genome-wide, seven previously reported common asthma variant loci and one previously reported lung function locus, near THSD4, reach significance. We replicate association of the lung function locus in a recently published GWAS of moderate-to-severe asthma patients. We additionally replicate the association of a previously reported rare (minor allele frequency < 1%) coding variant in IL33 and show significant enrichment of rare variant burden in genes from common variant allergic disease loci. Our findings highlight the contribution of lung function genetics to moderate-to-severe asthma risk, and provide initial rare variant support for associations with moderate-to-severe asthma risk at several candidate genes from common variant loci.

Coordinated concentration changes of transcripts and metabolites in Saccharomyces cerevisiae.

Metabolite concentrations can regulate gene expression, which can in turn regulate metabolic activity. The extent to which functionally related transcripts and metabolites show similar patterns of concentration changes, however, remains unestablished. We measure and analyze the metabolomic and transcriptional responses of Saccharomyces cerevisiae to carbon and nitrogen starvation. Our analysis demonstrates that transcripts and metabolites show coordinated response dynamics. Furthermore, metabolites and gene products whose concentration profiles are alike tend to participate in related biological processes. To identify specific, functionally related genes and metabolites, we develop an approach based on Bayesian integration of the joint metabolomic and transcriptomic data. This algorithm finds interactions by evaluating transcript-metabolite correlations in light of the experimental context in which they occur and the class of metabolite involved. It effectively predicts known enzymatic and regulatory relationships, including a gene-metabolite interaction central to the glycolytic-gluconeogenetic switch. This work provides quantitative evidence that functionally related metabolites and transcripts show coherent patterns of behavior on the genome scale and lays the groundwork for building gene-metabolite interaction networks directly from systems-level data.

Homeostatic adjustment and metabolic remodeling in glucose-limited yeast cultures.

We studied the physiological response to glucose limitation in batch and steady-state (chemostat) cultures of Saccharomyces cerevisiae by following global patterns of gene expression. Glucose-limited batch cultures of yeast go through two sequential exponential growth phases, beginning with a largely fermentative phase, followed by an essentially completely aerobic use of residual glucose and evolved ethanol. Judging from the patterns of gene expression, the state of the cells growing at steady state in glucose-limited chemostats corresponds most closely with the state of cells in batch cultures just before they undergo this "diauxic shift." Essentially the same pattern was found between chemostats having a fivefold difference in steady-state growth rate (the lower rate approximating that of the second phase respiratory growth rate in batch cultures). Although in both cases the cells in the chemostat consumed most of the glucose, in neither case did they seem to be metabolizing it primarily through respiration. Although there was some indication of a modest oxidative stress response, the chemostat cultures did not exhibit the massive environmental stress response associated with starvation that also is observed, at least in part, during the diauxic shift in batch cultures. We conclude that despite the theoretical possibility of a switch to fully aerobic metabolism of glucose in the chemostat under conditions of glucose scarcity, homeostatic mechanisms are able to carry out metabolic adjustment as if fermentation of the glucose is the preferred option until the glucose is entirely depleted. These results suggest that some aspect of actual starvation, possibly a component of the stress response, may be required for triggering the metabolic remodeling associated with the diauxic shift.

Analysis of protein-altering variants in telomerase genes and their association with MUC5B common variant status in patients with idiopathic pulmonary fibrosis: a candidate gene sequencing study.

BACKGROUND: Idiopathic pulmonary fibrosis (IPF) risk has a strong genetic component. Studies have implicated variations at several loci, including TERT, surfactant genes, and a single nucleotide polymorphism at chr11p15 (rs35705950) in the intergenic region between TOLLIP and MUC5B. Patients with IPF who have risk alleles at rs35705950 have longer survival from the time of IPF diagnosis than do patients homozygous for the non-risk allele, whereas patients with shorter telomeres have shorter survival times. We aimed to assess whether rare protein-altering variants in genes regulating telomere length are enriched in patients with IPF homozygous for the non-risk alleles at rs35705950. METHODS: Between Nov 1, 2014, and Nov 1, 2016, we assessed blood samples from patients aged 40 years or older and of European ancestry with sporadic IPF from three international phase 3 clinical trials (INSPIRE, CAPACITY, ASCEND), one phase 2 study (RIFF), and US-based observational studies (Vanderbilt Clinical Interstitial Lung Disease Registry and the UCSF Interstitial Lung Disease Clinic registry cohorts) at the Broad Institute (Cambridge, MA, USA) and Human Longevity (San Diego, CA, USA). We also assessed blood samples from non-IPF controls in several clinical trials. We did whole-genome sequencing to assess telomere length and identify rare protein-altering variants, stratified by rs35705950 genotype. We also assessed rare functional variation in TERT exons and compared telomere length and disease progression across genotypes. FINDINGS: We assessed samples from 1510 patients with IPF and 1874 non-IPF controls. 30 (3%) of 1046 patients with an rs35705950 risk allele had a rare protein-altering variant in TERT compared with 34 (7%) of 464 non-risk allele carriers (odds ratio 0·40 [95% CI 0·24-0·66], p=0·00039). Subsequent analyses identified enrichment of rare protein-altering variants in PARN and RTEL1, and rare variation in TERC in patients with IPF compared with controls. We expanded our study population to provide a more accurate estimation of rare variant frequency in these four loci, and to calculate telomere length. The proportion of patients with at least one rare variant in TERT, PARN, TERC, or RTEL1 was higher in patients with IPF than in controls (149 [9%] of 1739 patients vs 205 [2%] of 8645 controls, p=2·44 × 10-8). Patients with IPF who had a variant in any of the four identified telomerase component genes had telomeres that were 3·69-16·10% shorter than patients without a variant in any of the four genes and had an earlier mean age of disease onset than patients without one or more variants (65·1 years [SD 7·8] vs 67·1 years [7·9], p=0·004). In the placebo arms of clinical trials, shorter telomeres were significantly associated with faster disease progression (1·7% predicted forced vital capacity per kb per year, p=0·002). Pirfenidone had treatment benefit regardless of telomere length (p=4·24 × 10-8 for telomere length lower than the median, p=0·0044 for telomere length greater than the median). INTERPRETATION: Rare protein-altering variants in TERT, PARN, TERC, and RTEL1 are enriched in patients with IPF compared with controls, and, in the case of TERT, particularly in individuals without a risk allele at the rs35705950 locus. This suggests that multiple genetic factors contribute to sporadic IPF, which might implicate distinct mechanisms of pathogenesis and disease progression. FUNDING: Genentech, National Institutes of Health, Francis Family Foundation, Pulmonary Fibrosis Foundation, Nina Ireland Program for Lung Health, US Department of Veterans Affairs.

Mapping novel traits by array-assisted bulk segregant analysis in Saccharomyces cerevisiae.

We demonstrate a new method, microarray-assisted bulk segregant analysis, for mapping traits in yeast by genotyping pooled segregants. We apply a probabilistic model to the progeny of a single cross and as little as two microarray hybridizations to reliably map an auxotrophic marker, a Mendelian trait, and a major-effect quantitative trait locus.

Nutritional homeostasis in batch and steady-state culture of yeast.

We studied the physiological response to limitation by diverse nutrients in batch and steady-state (chemostat) cultures of S. cerevisiae. We found that the global pattern of transcription in steady-state cultures in limiting phosphate or sulfate is essentially identical to that of batch cultures growing in the same medium just before the limiting nutrient is completely exhausted. The massive stress response and complete arrest of the cell cycle that occurs when nutrients are fully exhausted in batch cultures is not observed in the chemostat, indicating that the cells in the chemostat are "poor, not starving." Similar comparisons using leucine or uracil auxotrophs limited on leucine or uracil again showed patterns of gene expression in steady-state closely resembling those of corresponding batch cultures just before they exhaust the nutrient. Although there is also a strong stress response in the auxotrophic batch cultures, cell cycle arrest, if it occurs at all, is much less uniform. Many of the differences among the patterns of gene expression between the four nutrient limitations are interpretable in light of known involvement of the genes in stress responses or in the regulation or execution of particular metabolic pathways appropriate to the limiting nutrient. We conclude that cells adjust their growth rate to nutrient availability and maintain homeostasis in the same way in batch and steady state conditions; cells in steady-state cultures are in a physiological condition normally encountered in batch cultures.

Constitutive NOTCH3 Signaling Promotes the Growth of Basal Breast Cancers.

Notch ligands signal through one of four receptors on neighboring cells to mediate cell-cell communication and control cell fate, proliferation, and survival. Although aberrant Notch activation has been implicated in numerous malignancies, including breast cancer, the importance of individual receptors in distinct breast cancer subtypes and the mechanisms of receptor activation remain unclear. Using a novel antibody to detect active NOTCH3, we report here that NOTCH3 signals constitutively in a panel of basal breast cancer cell lines and in more than one third of basal tumors. Selective inhibition of individual ligands revealed that this signal does not require canonical ligand induction. A NOTCH3 antagonist antibody inhibited growth of basal lines, whereas a NOTCH3 agonist antibody enhanced the transformed phenotype in vitro and in tumor xenografts. Transcriptomic analyses generated a Notch gene signature that included Notch pathway components, the oncogene c-Myc, and the mammary stem cell regulator Id4 This signature drove clustering of breast cancer cell lines and tumors into the common subtypes and correlated with the basal classification. Our results highlight an unexpected ligand-independent induction mechanism and suggest that constitutive NOTCH3 signaling can drive an oncogenic program in a subset of basal breast cancers. Cancer Res; 77(6); 1439-52. ©2017 AACR.

GMAP and GSNAP for Genomic Sequence Alignment: Enhancements to Speed, Accuracy, and Functionality.

The programs GMAP and GSNAP, for aligning RNA-Seq and DNA-Seq datasets to genomes, have evolved along with advances in biological methodology to handle longer reads, larger volumes of data, and new types of biological assays. The genomic representation has been improved to include linear genomes that can compare sequences using single-instruction multiple-data (SIMD) instructions, compressed genomic hash tables with fast access using SIMD instructions, handling of large genomes with more than four billion bp, and enhanced suffix arrays (ESAs) with novel data structures for fast access. Improvements to the algorithms have included a greedy match-and-extend algorithm using suffix arrays, segment chaining using genomic hash tables, diagonalization using segmental hash tables, and nucleotide-level dynamic programming procedures that use SIMD instructions and eliminate the need for F-loop calculations. Enhancements to the functionality of the programs include standardization of indel positions, handling of ambiguous splicing, clipping and merging of overlapping paired-end reads, and alignments to circular chromosomes and alternate scaffolds. The programs have been adapted for use in pipelines by integrating their usage into R/Bioconductor packages such as gmapR and HTSeqGenie, and these pipelines have facilitated the discovery of numerous biological phenomena.

Sustained Brown Fat Stimulation and Insulin Sensitization by a Humanized Bispecific Antibody Agonist for Fibroblast Growth Factor Receptor 1/ßKlotho Complex.

Dissipating excess calories as heat through therapeutic stimulation of brown adipose tissues (BAT) has been proposed as a potential treatment for obesity-linked disorders. Here, we describe the generation of a humanized effector-less bispecific antibody that activates fibroblast growth factor receptor (FGFR) 1/ßKlotho complex, a common receptor for FGF21 and FGF19. Using this molecule, we show that antibody-mediated activation of FGFR1/ßKlotho complex in mice induces sustained energy expenditure in BAT, browning of white adipose tissue, weight loss, and improvements in obesity-associated metabolic derangements including insulin resistance, hyperglycemia, dyslipidemia and hepatosteatosis. In mice and cynomolgus monkeys, FGFR1/ßKlotho activation increased serum high-molecular-weight adiponectin, which appears to contribute over time by enhancing the amplitude of the metabolic benefits. At the same time, insulin sensitization by FGFR1/ßKlotho activation occurs even before the onset of weight loss in a manner that is independent of adiponectin. Together, selective activation of FGFR1/ßKlotho complex with a long acting therapeutic antibody represents an attractive approach for the treatment of type 2 diabetes and other obesity-linked disorders through enhanced energy expenditure, insulin sensitization and induction of high-molecular-weight adiponectin.

A comprehensive transcriptional portrait of human cancer cell lines.

Tumor-derived cell lines have served as vital models to advance our understanding of oncogene function and therapeutic responses. Although substantial effort has been made to define the genomic constitution of cancer cell line panels, the transcriptome remains understudied. Here we describe RNA sequencing and single-nucleotide polymorphism (SNP) array analysis of 675 human cancer cell lines. We report comprehensive analyses of transcriptome features including gene expression, mutations, gene fusions and expression of non-human sequences. Of the 2,200 gene fusions catalogued, 1,435 consist of genes not previously found in fusions, providing many leads for further investigation. We combine multiple genome and transcriptome features in a pathway-based approach to enhance prediction of response to targeted therapeutics. Our results provide a valuable resource for studies that use cancer cell lines.

A high-content cellular senescence screen identifies candidate tumor suppressors, including EPHA3.

Activation of a cellular senescence program is a common response to prolonged oncogene activation or tumor suppressor loss, providing a physiological mechanism for tumor suppression in premalignant cells. The link between senescence and tumor suppression supports the hypothesis that a loss-of-function screen measuring bona fide senescence marker activation should identify candidate tumor suppressors. Using a high-content siRNA screening assay for cell morphology and proliferation measures, we identify 12 senescence-regulating kinases and determine their senescence marker signatures, including elevation of senescence-associated ß-galactosidase, DNA damage and p53 or p16 (INK4a) expression. Consistent with our hypothesis, SNP array CGH data supports loss of gene copy number of five senescence-suppressing genes across multiple tumor samples. One such candidate is the EPHA3 receptor tyrosine kinase, a gene commonly mutated in human cancer. We demonstrate that selected intracellular EPHA3 tumor-associated point mutations decrease receptor expression level and/or receptor tyrosine kinase (RTK) activity. Our study therefore describes a new strategy to mine for novel candidate tumor suppressors and provides compelling evidence that EPHA3 mutations may promote tumorigenesis only when key senescence-inducing pathways have been inactivated.

Identification and analysis of in vivo VEGF downstream markers link VEGF pathway activity with efficacy of anti-VEGF therapies.

PURPOSE: The aim of this study was to identify conserved pharmacodynamic and potential predictive biomarkers of response to anti-VEGF therapy using gene expression profiling in preclinical tumor models and in patients. EXPERIMENTAL DESIGN: Surrogate markers of VEGF inhibition [VEGF-dependent genes or VEGF-dependent vasculature (VDV)] were identified by profiling gene expression changes induced in response to VEGF blockade in preclinical tumor models and in human biopsies from patients treated with anti-VEGF monoclonal antibodies. The potential value of VDV genes as candidate predictive biomarkers was tested by correlating high or low VDV gene expression levels in pretreatment clinical samples with the subsequent clinical efficacy of bevacizumab (anti-VEGF)-containing therapy. RESULTS: We show that VDV genes, including direct and more distal VEGF downstream endothelial targets, enable detection of VEGF signaling inhibition in mouse tumor models and human tumor biopsies. Retrospective analyses of clinical trial data indicate that patients with higher VDV expression in pretreatment tumor samples exhibited improved clinical outcome when treated with bevacizumab-containing therapies. CONCLUSIONS: In this work, we identified surrogate markers (VDV genes) for in vivo VEGF signaling in tumors and showed clinical data supporting a correlation between pretreatment VEGF bioactivity and the subsequent efficacy of anti-VEGF therapy. We propose that VDV genes are candidate biomarkers with the potential to aid the selection of novel indications as well as patients likely to respond to anti-VEGF therapy. The data presented here define a diagnostic biomarker hypothesis based on translational research that warrants further evaluation in additional retrospective and prospective trials.

Slow growth induces heat-shock resistance in normal and respiratory-deficient yeast.

Yeast cells respond to a variety of environmental stresses, including heat shock and growth limitation. There is considerable overlap in these responses both from the point of view of gene expression patterns and cross-protection for survival. We performed experiments in which cells growing at different steady-state growth rates in chemostats were subjected to a short heat pulse. Gene expression patterns allowed us to partition genes whose expression responds to heat shock into subsets of genes that also respond to slow growth rate and those that do not. We found also that the degree of induction and repression of genes that respond to stress is generally weaker in respiratory deficient mutants, suggesting a role for increased respiratory activity in the apparent stress response to slow growth. Consistent with our gene expression results in wild-type cells, we found that cells growing more slowly are cross-protected for heat shock, i.e., better able to survive a lethal heat challenge. Surprisingly, however, we found no difference in cross-protection between respiratory-deficient and wild-type cells, suggesting induction of heat resistance at low growth rates is independent of respiratory activity, even though many of the changes in gene expression are not.

Coordination of growth rate, cell cycle, stress response, and metabolic activity in yeast.

We studied the relationship between growth rate and genome-wide gene expression, cell cycle progression, and glucose metabolism in 36 steady-state continuous cultures limited by one of six different nutrients (glucose, ammonium, sulfate, phosphate, uracil, or leucine). The expression of more than one quarter of all yeast genes is linearly correlated with growth rate, independent of the limiting nutrient. The subset of negatively growth-correlated genes is most enriched for peroxisomal functions, whereas positively correlated genes mainly encode ribosomal functions. Many (not all) genes associated with stress response are strongly correlated with growth rate, as are genes that are periodically expressed under conditions of metabolic cycling. We confirmed a linear relationship between growth rate and the fraction of the cell population in the G0/G1 cell cycle phase, independent of limiting nutrient. Cultures limited by auxotrophic requirements wasted excess glucose, whereas those limited on phosphate, sulfate, or ammonia did not; this phenomenon (reminiscent of the "Warburg effect" in cancer cells) was confirmed in batch cultures. Using an aggregate of gene expression values, we predict (in both continuous and batch cultures) an "instantaneous growth rate." This concept is useful in interpreting the system-level connections among growth rate, metabolism, stress, and the cell cycle.