Large-scale genetic perturbations reveal regulatory networks and an abundance of gene-specific repressors.

To understand regulatory systems, it would be useful to uniformly determine how different components contribute to the expression of all other genes. We therefore monitored mRNA expression genome-wide, for individual deletions of one-quarter of yeast genes, focusing on (putative) regulators. The resulting genetic perturbation signatures reflect many different properties. These include the architecture of protein complexes and pathways, identification of expression changes compatible with viability, and the varying responsiveness to genetic perturbation. The data are assembled into a genetic perturbation network that shows different connectivities for different classes of regulators. Four feed-forward loop (FFL) types are overrepresented, including incoherent type 2 FFLs that likely represent feedback. Systematic transcription factor classification shows a surprisingly high abundance of gene-specific repressors, suggesting that yeast chromatin is not as generally restrictive to transcription as is often assumed. The data set is useful for studying individual genes and for discovering properties of an entire regulatory system.

Shortening of membrane lipid acyl chains compensates for phosphatidylcholine deficiency in choline-auxotroph yeast.

Phosphatidylcholine (PC) is an abundant membrane lipid component in most eukaryotes, including yeast, and has been assigned multiple functions in addition to acting as building block of the lipid bilayer. Here, by isolating S. cerevisiae suppressor mutants that exhibit robust growth in the absence of PC, we show that PC essentiality is subject to cellular evolvability in yeast. The requirement for PC is suppressed by monosomy of chromosome XV or by a point mutation in the ACC1 gene encoding acetyl-CoA carboxylase. Although these two genetic adaptations rewire lipid biosynthesis in different ways, both decrease Acc1 activity, thereby reducing average acyl chain length. Consistently, soraphen A, a specific inhibitor of Acc1, rescues a yeast mutant with deficient PC synthesis. In the aneuploid suppressor, feedback inhibition of Acc1 through acyl-CoA produced by fatty acid synthase (FAS) results from upregulation of lipid synthesis. The results show that budding yeast regulates acyl chain length by fine-tuning the activities of Acc1 and FAS and indicate that PC evolved by benefitting the maintenance of membrane fluidity.

Functional overlap and regulatory links shape genetic interactions between signaling pathways.

To understand relationships between phosphorylation-based signaling pathways, we analyzed 150 deletion mutants of protein kinases and phosphatases in S. cerevisiae using DNA microarrays. Downstream changes in gene expression were treated as a phenotypic readout. Double mutants with synthetic genetic interactions were included to investigate genetic buffering relationships such as redundancy. Three types of genetic buffering relationships are identified: mixed epistasis, complete redundancy, and quantitative redundancy. In mixed epistasis, the most common buffering relationship, different gene sets respond in different epistatic ways. Mixed epistasis arises from pairs of regulators that have only partial overlap in function and that are coupled by additional regulatory links such as repression of one by the other. Such regulatory modules confer the ability to control different combinations of processes depending on condition or context. These properties likely contribute to the evolutionary maintenance of paralogs and indicate a way in which signaling pathways connect for multiprocess control.

EBV MicroRNA BART16 Suppresses Type I IFN Signaling.

Type I IFNs play critical roles in orchestrating the antiviral defense by inducing direct antiviral activities and shaping the adaptive immune response. Viruses have evolved numerous strategies to specifically interfere with IFN production or its downstream mediators, thereby allowing successful infection of the host to occur. The prototypic human gammaherpesvirus EBV, which is associated with infectious mononucleosis and malignant tumors, harbors many immune-evasion proteins that manipulate the adaptive and innate immune systems. In addition to proteins, the virus encodes >40 mature microRNAs for which the functions remain largely unknown. In this article, we identify EBV-encoded miR-BART16 as a novel viral immune-evasion factor that interferes with the type I IFN signaling pathway. miR-BART16 directly targets CREB-binding protein, a key transcriptional coactivator in IFN signaling, thereby inducing CREB-binding protein downregulation in EBV-transformed B cells and gastric carcinoma cells. miR-BART16 abrogates the production of IFN-stimulated genes in response to IFN-α stimulation and it inhibits the antiproliferative effect of IFN-α on latently infected BL cells. By obstructing the type I IFN-induced antiviral response, miR-BART16 provides a means to facilitate the establishment of latent EBV infection and enhance viral replication.

The specificity and topology of chromatin interaction pathways in yeast.

Packaging of DNA into chromatin has a profound impact on gene expression. To understand how changes in chromatin influence transcription, we analyzed 165 mutants of chromatin machinery components in Saccharomyces cerevisiae. mRNA expression patterns change in 80% of mutants, always with specific effects, even for loss of widespread histone marks. The data are assembled into a network of chromatin interaction pathways. The network is function based, has a branched, interconnected topology, and lacks strict one-to-one relationships between complexes. Chromatin pathways are not separate entities for different gene sets, but share many components. The study evaluates which interactions are important for which genes and predicts additional interactions, for example between Paf1C and Set3C, as well as a role for Mediator in subtelomeric silencing. The results indicate the presence of gene-dependent effects that go beyond context-dependent binding of chromatin factors and provide a framework for understanding how specificity is achieved through regulating chromatin.

FoxO6 affects Plxna4-mediated neuronal migration during mouse cortical development.

The forkhead transcription factor FoxO6 is prominently expressed during development of the murine neocortex. However, its function in cortical development is as yet unknown. We now demonstrate that cortical development is altered in FoxO6+/- and FoxO6-/- mice, showing migrating neurons halted in the intermediate zone. Using a FoxO6-directed siRNA approach, we substantiate the requirement of FoxO6 for a correct radial migration in the developing neocortex. Subsequent genome-wide transcriptome analysis reveals altered expression of genes involved in cell adhesion, axon guidance, and gliogenesis upon silencing of FoxO6 We then show that FoxO6 binds to DAF-16-binding elements in the Plexin A4 (Plxna4) promoter region and affects Plxna4 expression. Finally, ectopic Plxna4 expression restores radial migration in FoxO6+/- and siRNA-mediated knockdown models. In conclusion, the presented data provide insights into the molecular mechanisms whereby transcriptional programs drive cortical development.

The genomic landscape of compensatory evolution.

Adaptive evolution is generally assumed to progress through the accumulation of beneficial mutations. However, as deleterious mutations are common in natural populations, they generate a strong selection pressure to mitigate their detrimental effects through compensatory genetic changes. This process can potentially influence directions of adaptive evolution by enabling evolutionary routes that are otherwise inaccessible. Therefore, the extent to which compensatory mutations shape genomic evolution is of central importance. Here, we studied the capacity of the baker's yeast genome to compensate the complete loss of genes during evolution, and explored the long-term consequences of this process. We initiated laboratory evolutionary experiments with over 180 haploid baker's yeast genotypes, all of which initially displayed slow growth owing to the deletion of a single gene. Compensatory evolution following gene loss was rapid and pervasive: 68% of the genotypes reached near wild-type fitness through accumulation of adaptive mutations elsewhere in the genome. As compensatory mutations have associated fitness costs, genotypes with especially low fitnesses were more likely to be subjects of compensatory evolution. Genomic analysis revealed that as compensatory mutations were generally specific to the functional defect incurred, convergent evolution at the molecular level was extremely rare. Moreover, the majority of the gene expression changes due to gene deletion remained unrestored. Accordingly, compensatory evolution promoted genomic divergence of parallel evolving populations. However, these different evolutionary outcomes are not phenotypically equivalent, as they generated diverse growth phenotypes across environments. Taken together, these results indicate that gene loss initiates adaptive genomic changes that rapidly restores fitness, but this process has substantial pleiotropic effects on cellular physiology and evolvability upon environmental change. Our work also implies that gene content variation across species could be partly due to the action of compensatory evolution rather than the passive loss of genes.

Mapping of a FEB3 homologous febrile seizure locus on mouse chromosome 2 containing candidate genes Scn1a and Scn3a.

Febrile seizures (FS) are the most common seizure type in children. Recurrent FS are a risk factor for developing temporal lobe epilepsy later in life and are known to have a strong genetic component. Experimental FS (eFS) can be elicited in mice by warm-air induced hyperthermia. We used this model to screen the chromosome substitution strain (CSS) panel derived from C57BL/6J and A/J for FS susceptibility and identified C57BL/6J-Chr2A /NaJ (CSS2), as the strain with the strongest FS susceptibility phenotype. The aim of this study was to map FS susceptibility loci and select candidate genes on mouse chromosome 2. We generated an F2 population by intercrossing the hybrids (F1 ) that were derived from CSS2 and C57BL/6J mice. All CSS2-F2 individuals were genotyped and phenotyped for eFS susceptibility, and QTL analysis was performed. Candidate gene selection was based on bioinformatics analyses and differential brain expression between CSS2 and C57BL/6J strains determined by microarray analysis. Genetic mapping of the eFS susceptibility trait identified two significant loci: FS-QTL2a (LOD-score 3.6) and FS-QTL2b (LOD-score 6.2). FS-QTL2a contained 44 genes expressed in the brain at post natal day 14. Four of these (Arl6ip6, Cytip, Fmnl2 Ifih1) contained a non-synonymous SNP comparing CSS2 and C57BL/6J, six genes (March7, Nr4a2, Gpd2, Grb14, Scn1a, Scn3a) were differentially expressed between these strains. A region within FS-QTL2a is homologous to the human FEB3 locus. The fact that we identify mouse FS-QTL2a with high FEB3 homology is strong support for the validity of the eFS mouse model to study genetics of human FS.

A high-resolution gene expression atlas of epistasis between gene-specific transcription factors exposes potential mechanisms for genetic interactions.

BACKGROUND: Genetic interactions, or non-additive effects between genes, play a crucial role in many cellular processes and disease. Which mechanisms underlie these genetic interactions has hardly been characterized. Understanding the molecular basis of genetic interactions is crucial in deciphering pathway organization and understanding the relationship between genotype, phenotype and disease. RESULTS: To investigate the nature of genetic interactions between gene-specific transcription factors (GSTFs) in Saccharomyces cerevisiae, we systematically analyzed 72 GSTF pairs by gene expression profiling double and single deletion mutants. These pairs were selected through previously published growth-based genetic interactions as well as through similarity in DNA binding properties. The result is a high-resolution atlas of gene expression-based genetic interactions that provides systems-level insight into GSTF epistasis. The atlas confirms known genetic interactions and exposes new ones. Importantly, the data can be used to investigate mechanisms that underlie individual genetic interactions. Two molecular mechanisms are proposed, "buffering by induced dependency" and "alleviation by derepression". CONCLUSIONS: These mechanisms indicate how negative genetic interactions can occur between seemingly unrelated parallel pathways and how positive genetic interactions can indirectly expose parallel rather than same-pathway relationships. The focus on GSTFs is important for understanding the transcription regulatory network of yeast as it uncovers details behind many redundancy relationships, some of which are completely new. In addition, the study provides general insight into the complex nature of epistasis and proposes mechanistic models for genetic interactions, the majority of which do not fall into easily recognizable within- or between-pathway relationships.

A mRNA landscape of bovine embryos after standard and MAPK-inhibited culture conditions: a comparative analysis.

BACKGROUND: Genes and signalling pathways involved in pluripotency have been studied extensively in mouse and human pre-implantation embryos and embryonic stem (ES) cells. The unsuccessful attempts to generate ES cell lines from other species including cattle suggests that other genes and pathways are involved in maintaining pluripotency in these species. To investigate which genes are involved in bovine pluripotency, expression profiles were generated from morula, blastocyst, trophectoderm and inner cell mass (ICM) samples using microarray analysis. As MAPK inhibition can increase the NANOG/GATA6 ratio in the inner cell mass, additionally blastocysts were cultured in the presence of a MAPK inhibitor and changes in gene expression in the inner cell mass were analysed. RESULTS: Between morula and blastocyst 3,774 genes were differentially expressed and the largest differences were found in blastocyst up-regulated genes. Gene ontology (GO) analysis shows lipid metabolic process as the term most enriched with genes expressed at higher levels in blastocysts. Genes with higher expression levels in morulae were enriched in the RNA processing GO term. Of the 497 differentially expressed genes comparing ICM and TE, the expression of NANOG, SOX2 and POU5F1 was increased in the ICM confirming their evolutionary preserved role in pluripotency. Several genes implicated to be involved in differentiation or fate determination were also expressed at higher levels in the ICM. Genes expressed at higher levels in the ICM were enriched in the RNA splicing and regulation of gene expression GO term. Although NANOG expression was elevated upon MAPK inhibition, SOX2 and POU5F1 expression showed little increase. Expression of other genes in the MAPK pathway including DUSP4 and SPRY4, or influenced by MAPK inhibition such as IFNT, was down-regulated. CONCLUSION: The data obtained from the microarray studies provide further insight in gene expression during bovine embryonic development. They show an expression profile in pluripotent cells that indicates a pluripotent, epiblast-like state. The inability to culture ICM cells as stem cells in the presence of an inhibitor of MAPK activity together with the reported data indicates that MAPK inhibition alone is not sufficient to maintain a pluripotent character in bovine cells.

Cell cycle population effects in perturbation studies.

Growth condition perturbation or gene function disruption are commonly used strategies to study cellular systems. Although it is widely appreciated that such experiments may involve indirect effects, these frequently remain uncharacterized. Here, analysis of functionally unrelated Saccharyomyces cerevisiae deletion strains reveals a common gene expression signature. One property shared by these strains is slower growth, with increased presence of the signature in more slowly growing strains. The slow growth signature is highly similar to the environmental stress response (ESR), an expression response common to diverse environmental perturbations. Both environmental and genetic perturbations result in growth rate changes. These are accompanied by a change in the distribution of cells over different cell cycle phases. Rather than representing a direct expression response in single cells, both the slow growth signature and ESR mainly reflect a redistribution of cells over different cell cycle phases, primarily characterized by an increase in the G1 population. The findings have implications for any study of perturbation that is accompanied by growth rate changes. Strategies to counter these effects are presented and discussed.

E2F7 represses a network of oscillating cell cycle genes to control S-phase progression.

E2F transcription factors are known to be important for timely activation of G(1)/S and G(2)/M genes required for cell cycle progression, but transcriptional mechanisms for deactivation of cell cycle-regulated genes are unknown. Here, we show that E2F7 is highly expressed during mid to late S-phase, occupies promoters of G(1)/S-regulated genes and represses their transcription. ChIP-seq analysis revealed that E2F7 binds preferentially to genomic sites containing the TTCCCGCC motif, which closely resembles the E2F consensus site. We identified 89 target genes that carry E2F7 binding sites close to the transcriptional start site and that are directly repressed by short-term induction of E2F7. Most of these target genes are known to be activated by E2Fs and are involved in DNA replication, metabolism and DNA repair. Importantly, induction of E2F7 during G(0)-G(1)/S resulted in S-phase arrest and DNA damage, whereas expression of E2F7 during G(2)/M failed to disturb cell cycle progression. These findings provide strong evidence that E2F7 directly controls the downswing of oscillating G(1)/S genes during S-phase progression.

Single-cell transcriptomics reveals immune suppression and cell states predictive of patient outcomes in rhabdomyosarcoma.

Paediatric rhabdomyosarcoma (RMS) is a soft tissue malignancy of mesenchymal origin that is thought to arise as a consequence of derailed myogenic differentiation. Despite intensive treatment regimens, the prognosis for high-risk patients remains dismal. The cellular differentiation states underlying RMS and how these relate to patient outcomes remain largely elusive. Here, we use single-cell mRNA sequencing to generate a transcriptomic atlas of RMS. Analysis of the RMS tumour niche reveals evidence of an immunosuppressive microenvironment. We also identify a putative interaction between NECTIN3 and TIGIT, specific to the more aggressive fusion-positive (FP) RMS subtype, as a potential cause of tumour-induced T-cell dysfunction. In malignant RMS cells, we define transcriptional programs reflective of normal myogenic differentiation and show that these cellular differentiation states are predictive of patient outcomes in both FP RMS and the less aggressive fusion-negative subtype. Our study reveals the potential of therapies targeting the immune microenvironment of RMS and suggests that assessing tumour differentiation states may enable a more refined risk stratification.

Mesenchymal tumor organoid models recapitulate rhabdomyosarcoma subtypes.

Rhabdomyosarcomas (RMS) are mesenchyme-derived tumors and the most common childhood soft tissue sarcomas. Treatment is intense, with a nevertheless poor prognosis for high-risk patients. Discovery of new therapies would benefit from additional preclinical models. Here, we describe the generation of a collection of 19 pediatric RMS tumor organoid (tumoroid) models (success rate of 41%) comprising all major subtypes. For aggressive tumors, tumoroid models can often be established within 4-8 weeks, indicating the feasibility of personalized drug screening. Molecular, genetic, and histological characterization show that the models closely resemble the original tumors, with genetic stability over extended culture periods of up to 6 months. Importantly, drug screening reflects established sensitivities and the models can be modified by CRISPR/Cas9 with TP53 knockout in an embryonal RMS model resulting in replicative stress drug sensitivity. Tumors of mesenchymal origin can therefore be used to generate organoid models, relevant for a variety of preclinical and clinical research questions.

Tight control of MEK-ERK activation is essential in regulating proliferation, survival, and cytokine production of CD34+-derived neutrophil progenitors.

A plethora of extracellular stimuli regulate growth, survival, and differentiation responses through activation of the MEK-ERK MAPK signaling module. Using CD34+ hematopoietic progenitor cells, we describe a novel role for the MEK-ERK signaling module in the regulation of proliferation, survival, and cytokine production during neutrophil differentiation. Addition of the specific MEK1/2 inhibitor U0126 resulted in decreased proliferation of neutrophil progenitors. Conversely, transient activation of a conditionally active MEK1 mutant resulted in the expansion of progenitor cells, which thereafter differentiated normally into mature neutrophils. In contrast, chronic MEK1 activation was found to induce cell death of CD34+ neutrophil progenitors. Microarray analysis of CD34+ progenitor cells showed that activation of MEK1 resulted in changes in expression of a variety of cell-cycle modulating genes. Furthermore, conditional activation of MEK1 resulted in a dramatic increase in the expression of mRNA transcripts encoding a large number of hematopoietic cytokines, chemokines, and growth factors. These findings identify a novel role for MEK-ERK signaling in regulating the balance between proliferation and apoptosis during neutrophil differentiation, and they suggest the need for tight control of MEK-ERK activation to prevent the development of bone marrow failure.

Type I interferon receptor-independent and -dependent host transcriptional responses to mouse hepatitis coronavirus infection in vivo.

BACKGROUND: The role of type I IFNs in protecting against coronavirus (CoV) infections is not fully understood. While CoVs are poor inducers of type I IFNs in tissue culture, several studies have demonstrated the importance of the type I IFN response in controlling MHV infection in animals. The protective effectors against MHV infection are, however, still unknown. RESULTS: In order to get more insight into the antiviral gene expression induced in the brains of MHV-infected mice, we performed whole-genome expression profiling. Three different mouse strains, differing in their susceptibility to infection with MHV, were used. In BALB/c mice, which display high viral loads but are able to control the infection, 57 and 121 genes were significantly differentially expressed (> or = 1.5 fold change) upon infection at 2 and 5 days post infection, respectively. Functional association network analyses demonstrated a strong type I IFN response, with Irf1 and Irf7 as the central players. At 5 days post infection, a type II IFN response also becomes apparent. Both the type I and II IFN response, which were more pronounced in mice with a higher viral load, were not observed in 129SvEv mice, which are much less susceptible to infection with MHV. 129SvEv mice lacking the type I interferon receptor (IFNAR-/-), however, were not able to control the infection. Gene expression profiling of these mice identified type I IFN-independent responses to infection, with IFN-gamma as the central player. As the BALB/c and the IFNAR-/- 129SvEv mice demonstrated very similar viral loads in their brains, we also compared their gene expression profiles upon infection with MHV in order to identify type I IFN-dependent transcriptional responses. Many known IFN-inducible genes were detected, several of which have previously been shown to play an important protective role against virus infections. We speculate that the additional type I IFN-dependent genes that we discovered may also be important for protection against MHV infection. CONCLUSION: Transcriptional profiling of mice infected with MHV demonstrated the induction of a robust IFN response, which correlated with the viral load. Profiling of IFNAR-/- mice allowed us to identify type I IFN-independent and -dependent responses. Overall, this study broadens our present knowledge of the type I and II IFN-mediated effector responses during CoV infection in vivo.

Maintenance of head and neck tumor gene expression profiles upon lymph node metastasis.

Spread of cancer and development of solid metastases at distant sites is the main cause of cancer-related deaths. To understand and treat metastases, it is important to determine at which stages the most pivotal steps for development of metastases occur. In head and neck squamous cell carcinoma (HNSCC), metastasis nearly always occurs first in local lymph nodes before development of distant metastasis. Here, we have investigated gene expression patterns in HNSCC lymph node metastases using DNA microarrays. Several types of analyses show that the gene expression patterns in lymph node metastases are most similar to the corresponding primary tumors from which they arose, as long as samples contain sufficient proportions of tumor cells. Strikingly, gene expression patterns of metastatic primary HNSCC are largely maintained upon spread to the lymph node. Only a single gene, metastasis-associated gene 1 (MTA1), was found to show consistently changed expression between a large number of matched primary tumor-lymph node metastasis pairs. The maintained expression pattern includes the predictive signature for HNSCC lymph node metastasis. These results underscore the importance of the primary tumor gene expression profile for development and treatment of metastasis. The findings also agree with the concept that disseminated cancer cells alter the surrounding tissue into a metastatic environment that resembles the primary tumor microenvironment.

Growth plate expression profiling: Large and small breed dogs provide new insights in endochondral bone formation.

The difference in the adult height of mammals, and hence in endochondral bone formation, is not yet fully understood and may serve to identify targets for bone and cartilage regeneration. In line with this hypothesis, the intra-species disparity between the adult height of Great Danes and Miniature Poodles was investigated at a transcriptional level. Microarray analysis of the growth plate of five Great Danes and five Miniature Poodles revealed 2,981 unique genes that were differentially expressed, including many genes with an unknown role in skeletal development. A signaling pathway impact analysis indicated activation of the cell cycle, extracellular matrix receptor interaction and the tight junction pathway, and inhibition of pathways associated with inflammation and the complement cascade. In additional validation steps, the gene expression profile of the separate growth plate zones for both dog breeds were determined. Given that the BMP signaling is known for its crucial role in skeletal development and fracture healing, and BMP-2 is used in orthopaedic and spine procedures for bone augmentation, further investigations concentrated on the BMP pathway.The canonical BMP-2 and BMP-6 signaling pathway was activated in the Great Danes compared to Miniature Poodles. In conclusion, investigating the differential expression of genes involved in endochondral bone formation in small and large breed dogs, could be a game changing strategy to provide new insights in growth plate development and identify new targets for bone and cartilage regeneration. © 2017 The Authors. Journal of Orthopaedic Research® published by Wiley Periodicals, Inc. on behalf of the Orthopaedic Research Society. J Orthop Res 36:138-148, 2018.

Genome-wide based model predicting recovery from portosystemic shunting after liver shunt attenuation in dogs.

BACKGROUND: In dogs with congenital portosystemic shunt (CPSS), recovery after surgical CPSS attenuation is difficult to predict. OBJECTIVES: Our aim was to build a model with plasma albumin concentration and mRNA expression levels of hepatic gene products as predictors of recovery from portosystemic shunting after surgery. ANIMALS: Seventy-three client-owned dogs referred for surgical attenuation of CPSS. METHODS: A prediction model was constructed using 2 case-control studies of recovered and nonrecovered dogs after surgical CPSS attenuation. In the 1st study, a dog-specific gene expression microarray analysis was used to compare mRNA expression in intraoperatively collected liver tissue between 23 recovered and 23 nonrecovered dogs. In the 2nd study, preoperative plasma albumin concentration and the expression of microarray-selected genes were confirmed by RT-qPCR in intraoperatively collected liver samples of 31 recovered and 31 nonrecovered dogs, including 35 dogs from the 1st study. RESULTS: In the 1st study, 43 genes were differently expressed in recovered and nonrecovered dogs. The mean preoperative plasma albumin concentration in recovered dogs was higher compared to nonrecovered dogs (23 and 19 g/L, respectively; P = .004). The best fitting prediction model in the 2nd study included preoperative plasma albumin concentration and intraoperative DHDH, ERLEC1, and LYSMD2 gene expression levels. CONCLUSION AND CLINICAL IMPORTANCE: A preclinical model was constructed using preoperative plasma albumin concentration and intraoperative hepatic mRNA expression of 3 genes that were unbiasedly selected from the genome to predict recovery from portosystemic shunting after shunt ligation. Further development is essential for clinical application.

FOXO Transcription Factors Both Suppress and Support Breast Cancer Progression.

FOXO transcription factors are regulators of cellular homeostasis and putative tumor suppressors, yet the role of FOXO in cancer progression remains to be determined. The data on FOXO function, particularly for epithelial cancers, are fragmentary and come from studies that focused on isolated aspects of cancer. To clarify the role of FOXO in epithelial cancer progression, we characterized the effects of inducible FOXO activation and loss in a mouse model of metastatic invasive lobular carcinoma. Strikingly, either activation or loss of FOXO function suppressed tumor growth and metastasis. We show that the multitude of cellular processes critically affected by FOXO function include proliferation, survival, redox homeostasis, and PI3K signaling, all of which must be carefully balanced for tumor cells to thrive.Significance: FOXO proteins are not solely tumor suppressors, but also support tumor growth and metastasis by regulating a multitude of cellular processes essential for tumorigenesis. Cancer Res; 78(9); 2356-69. ©2018 AACR.