Influenza A defective viral genomes and non-infectious particles are increased by host PI3K inhibition via anti-cancer drug alpelisib.

RNA viruses produce abundant defective viral genomes during replication, setting the stage for interactions between viral genomes that alter the course of pathogenesis. Harnessing these interactions to develop antivirals has become a recent goal of intense research focus. Despite decades of research, the mechanisms that regulate the production and interactions of Influenza A defective viral genomes are still unclear. The role of the host is essentially unexplored; specifically, it remains unknown whether host metabolism can influence the formation of defective viral genomes and the particles that house them. To address this question, we manipulated host cell anabolic signaling activity and monitored the production of defective viral genomes and particles by A/H1N1 and A/H3N2 strains, using a combination of single-cell immunofluorescence quantification, third-generation long-read sequencing, and the cluster-forming assay, a method we developed to titer defective and fully-infectious particles simultaneously. Here we show that alpelisib (Piqray), a highly selective inhibitor of mammalian Class 1a phosphoinositide-3 kinase (PI3K) receptors, significantly changed the proportion of defective particles and viral genomes (specifically deletion-containing viral genomes) in a strain-specific manner, under conditions that minimize multiple cycles of replication. Alpelisib pre-treatment of cells led to an increase in defective particles in the A/H3N2 strain, while the A/H1N1 strain showed a decrease in total viral particles. In the same infections, we found that defective viral genomes of polymerase and antigenic segments increased in the A/H1N1 strain, while the total particles decreased suggesting defective interference. We also found that the average deletion size in polymerase complex viral genomes increased in both the A/H3N2 and A/H1N1 strains. The A/H1N1 strain, additionally showed a dose-dependent increase in total number of defective viral genomes. In sum, we provide evidence that host cell metabolism can increase the production of defective viral genomes and particles at an early stage of infection, shifting the makeup of the infection and potential interactions among virions. Given that Influenza A defective viral genomes can inhibit pathogenesis, our study presents a new line of investigation into metabolic states associated with less severe flu infection and the potential induction of these states with metabolic drugs.

Spatiotemporal Clusters of ERK Activity Coordinate Cytokine-induced Inflammatory Responses in Human Airway Epithelial Cells.

RATIONALE: Spatially coordinated ERK signaling events ("SPREADs") transmit radially from a central point to adjacent cells via secreted ligands for EGFR and other receptors. SPREADs maintain homeostasis in non-pulmonary epithelia, but it is unknown whether they play a role in the airway epithelium or are dysregulated in inflammatory disease. OBJECTIVES: (1) To characterize spatiotemporal ERK activity in response to pro-inflammatory ligands, and (2) to assess pharmacological and metabolic regulation of cytokine-mediated SPREADs. METHODS: SPREADs were measured by live-cell ERK biosensors in human bronchial epithelial cell lines (HBE1 and 16HBE) and primary human bronchial epithelial (pHBE) cells, in both submerged and biphasic Air-Liquid Interface (ALI) culture conditions (i.e., differentiated cells). Cells were exposed to pro-inflammatory cytokines relevant to asthma and chronic obstructive pulmonary disease (COPD), and to pharmacological treatments (gefitinib, tocilizumab, hydrocortisone) and metabolic modulators (insulin, 2-deoxyglucose) to probe the airway epithelial mechanisms of SPREADs. Phospho-STAT3 immunofluorescence was used to measure localized inflammatory responses to IL-6. RESULTS: Pro-inflammatory cytokines significantly increased the frequency of SPREADs. Notably, differentiated pHBE cells display increased SPREAD frequency that coincides with airway epithelial barrier breakdown. SPREADs correlate with IL-6 peptide secretion and localized pSTAT3. Hydrocortisone, inhibitors of receptor signaling, and suppression of metabolic function decreased SPREAD occurrence. CONCLUSIONS: Pro-inflammatory cytokines modulate SPREADs in human airway epithelial cells via both secreted EGFR and IL6R ligands. SPREADs correlate with changes in epithelial barrier permeability, implying a role for spatiotemporal ERK signaling in barrier homeostasis and dysfunction during inflammation. The involvement of SPREADs in airway inflammation suggests a novel signaling mechanism that could be exploited clinically to supplement corticosteroid treatment for asthma and COPD.

Deciphering the History of ERK Activity from Fixed-Cell Immunofluorescence Measurements.

The Ras/ERK pathway drives cell proliferation and other oncogenic behaviors, and quantifying its activity in situ is of high interest in cancer diagnosis and therapy. Pathway activation is often assayed by measuring phosphorylated ERK. However, this form of measurement overlooks dynamic aspects of signaling that can only be observed over time. In this study, we combine a live, single-cell ERK biosensor approach with multiplexed immunofluorescence staining of downstream target proteins to ask how well immunostaining captures the dynamic history of ERK activity. Combining linear regression, machine learning, and differential equation models, we develop an interpretive framework for immunostains, in which Fra-1 and pRb levels imply long term activation of ERK signaling, while Egr-1 and c-Myc indicate recent activation. We show that this framework can distinguish different classes of ERK dynamics within a heterogeneous population, providing a tool for annotating ERK dynamics within fixed tissues.

Applications of Plant-Made Fibroblast Growth Factor for Human Pluripotent Stem Cells.

Human embryonic stem cells (hESCs) and induced pluripotent stem cells (hiPSCs) hold great potential in regenerative medicine. These cells can be expanded indefinitely in theory and are able to differentiate into different types of cells for cell therapies, drug screening, and basic biology studies. The reliable and effective propagation of hESCs and hiPSCs is important for their downstream applications. Basic fibroblast growth factor (bFGF) is critical to hESCs and hiPSCs for maintaining their pluripotency. Plant-produced growth factors are safe to use without potential contamination of infectious viruses and are less expensive to produce. In this study, we used rice cell-made basic fibroblast growth factor (RbFGF) to propagate hESCs and hiPSCs for at least eight passages. Both hESCs and hiPSCs cultured with RbFGF not only maintained the morphology but also the specific expression (OCT4, SSEA4, SOX2, and TRA-1-60) of PSCs, similar to those cultured with the commercial Escherichia coli-produced bFGF. Furthermore, both gene chip-based PluriTest and TaqMan hPSC Scorecard pluripotency analysis demonstrated the pluripotent expression profile of the hESCs cultured with RbFGF. In vitro trilineage assays further showed that these hESCs and hiPSCs cultured on RbFGF were capable of giving rise to cell derivatives of ectoderm, mesoderm, and endoderm, further demonstrating their pluripotency. Finally, chromosome stability was also maintained in hESCs cultured with RbFGF as demonstrated by normal karyotypes. This study suggests broad applications for plant-made growth factors in stem cell culture and regenerative medicine.

A guide to ERK dynamics, part 1: mechanisms and models.

Extracellular signal-regulated kinase (ERK) has long been studied as a key driver of both essential cellular processes and disease. A persistent question has been how this single pathway is able to direct multiple cell behaviors, including growth, proliferation, and death. Modern biosensor studies have revealed that the temporal pattern of ERK activity is highly variable and heterogeneous, and critically, that these dynamic differences modulate cell fate. This two-part review discusses the current understanding of dynamic activity in the ERK pathway, how it regulates cellular decisions, and how these cell fates lead to tissue regulation and pathology. In part 1, we cover the optogenetic and live-cell imaging technologies that first revealed the dynamic nature of ERK, as well as current challenges in biosensor data analysis. We also discuss advances in mathematical models for the mechanisms of ERK dynamics, including receptor-level regulation, negative feedback, cooperativity, and paracrine signaling. While hurdles still remain, it is clear that higher temporal and spatial resolution provide mechanistic insights into pathway circuitry. Exciting new algorithms and advanced computational tools enable quantitative measurements of single-cell ERK activation, which in turn inform better models of pathway behavior. However, the fact that current models still cannot fully recapitulate the diversity of ERK responses calls for a deeper understanding of network structure and signal transduction in general.

A guide to ERK dynamics, part 2: downstream decoding.

Signaling by the extracellular signal-regulated kinase (ERK) pathway controls many cellular processes, including cell division, death, and differentiation. In this second installment of a two-part review, we address the question of how the ERK pathway exerts distinct and context-specific effects on multiple processes. We discuss how the dynamics of ERK activity induce selective changes in gene expression programs, with insights from both experiments and computational models. With a focus on single-cell biosensor-based studies, we summarize four major functional modes for ERK signaling in tissues: adjusting the size of cell populations, gradient-based patterning, wave propagation of morphological changes, and diversification of cellular gene expression states. These modes of operation are disrupted in cancer and other related diseases and represent potential targets for therapeutic intervention. By understanding the dynamic mechanisms involved in ERK signaling, there is potential for pharmacological strategies that not only simply inhibit ERK, but also restore functional activity patterns and improve disease outcomes.

Modeling elucidates context dependence in adipose regulation.

Single-cell data and computational simulations reveal the dynamics of the transcription factors HIF1α and PPARγ during adipocyte differentiation and maturation. Modeling feedback within this network predicts a HIF1α-mediated choice between lipid accumulation and incomplete differentiation. In vitro experiments support this model, with implications for adipose dynamics in metabolic disorders involving hypoxia.

Continuous sensing of nutrients and growth factors by the mTORC1-TFEB axis.

mTORC1 senses nutrients and growth factors and phosphorylates downstream targets, including the transcription factor TFEB, to coordinate metabolic supply and demand. These functions position mTORC1 as a central controller of cellular homeostasis, but the behavior of this system in individual cells has not been well characterized. Here, we provide measurements necessary to refine quantitative models for mTORC1 as a metabolic controller. We developed a series of fluorescent protein-TFEB fusions and a multiplexed immunofluorescence approach to investigate how combinations of stimuli jointly regulate mTORC1 signaling at the single-cell level. Live imaging of individual MCF10A cells confirmed that mTORC1-TFEB signaling responds continuously to individual, sequential, or simultaneous treatment with amino acids and the growth factor insulin. Under physiologically relevant concentrations of amino acids, we observe correlated fluctuations in TFEB, AMPK, and AKT signaling that indicate continuous activity adjustments to nutrient availability. Using partial least squares regression modeling, we show that these continuous gradations are connected to protein synthesis rate via a distributed network of mTORC1 effectors, providing quantitative support for the qualitative model of mTORC1 as a homeostatic controller and clarifying its functional behavior within individual cells.

Cellular transformation promotes the incorporation of docosahexaenoic acid into the endolysosome-specific lipid bis(monoacylglycerol)phosphate in breast cancer.

Bis(monoacylglycero)phosphates (BMPs), a class of lipids highly enriched within endolysosomal organelles, are key components of the lysosomal intraluminal vesicles responsible for activating sphingolipid catabolic enzymes. While BMPs are understudied relative to other phospholipids, recent reports associate BMP dysregulation with a variety of pathological states including neurodegenerative diseases and lysosomal storage disorders. Since the dramatic lysosomal remodeling characteristic of cellular transformation could impact BMP abundance and function, we employed untargeted lipidomics approaches to identify and quantify BMP species in several in vitro and in vivo models of breast cancer and comparative non-transformed cells and tissues. We observed lower BMP levels within transformed cells relative to normal cells, and consistent enrichment of docosahexaenoic acid (22:6) fatty acyl chain-containing BMP species in both human- and mouse-derived mammary tumorigenesis models. Our functional analysis points to a working model whereby 22:6 BMPs serve as reactive oxygen species scavengers in tumor cells, protecting lysosomes from oxidant-induced lysosomal membrane permeabilization. Our findings suggest that breast tumor cells might divert polyunsaturated fatty acids into BMP lipids as part of an adaptive response to protect their lysosomes from elevated reactive oxygen species levels, and raise the possibility that BMP-mediated lysosomal protection is a tumor-specific vulnerability that may be exploited therapeutically.

Protocol for electrical modulation of ERK signaling dynamics in cell lines.

Here, we describe a protocol for modulating the dynamics of the extracellular signal-regulated kinases (ERK) pathway in a customized alternating current (AC) electric field stimulation chamber. We use an ERK translocation reporter that can accurately represent the intracellular ERK activity in real time without chemical agents or gene disruption. ERK activation is assessed by comparing the relative intensity of nuclear fluorescence to cytosolic fluorescence in live-cell conditions. The approach can be applied to other signaling pathways as well. For complete details on the use and execution of this protocol, please refer to Guo et al. (2021).

ERK signaling dynamics: Lights, camera, transduction.

Three-dimensional mammary epithelial acini are a model for understanding how microenvironment-driven signaling coordinates cell behavior and tissue morphogenesis. In this issue of Developmental Cell, Ender et al. use live-cell imaging to capture dynamic spatiotemporal patterns of ERK activity that instruct cell migration and survival fates in developing acini.

Electrically synchronizing and modulating the dynamics of ERK activation to regulate cell fate.

Intracellular signaling dynamics play fundamental roles in cell biology. Precise modulation of the amplitude, duration, and frequency of signaling activation will be a powerful approach to investigate molecular mechanisms as well as to engineer signaling to control cell behaviors. Here, we showed a practical approach to achieve precise amplitude modulation (AM), frequency modulation (FM), and duration modulation (DM) of MAP kinase activation. Alternating current (AC) electrical stimulation induced synchronized ERK activation. Amplitude and duration of ERK activation were controlled by varying stimulation strength and duration. ERK activation frequencies were arbitrarily modulated with trains of short AC applications with accurately defined intervals. Significantly, ERK dynamics coded by well-designed AC can rewire PC12 cell fate independent of growth factors. This technique can be used to synchronize and modulate ERK activation dynamics, thus would offer a practical way to control cell behaviors in vivo without the use of biochemical agents or genetic manipulation.

Entosis is induced by ultraviolet radiation.

Entosis is a cell death mechanism that is executed through neighbor cell ingestion and killing that occurs in cancer tissues and during development. Here, we identify JNK and p38 stress-activated kinase signaling as an inducer of entosis in cells exposed to ultraviolet (UV) radiation. Cells with high levels of stress signaling are ingested and killed by those with low levels, a result of heterogeneity arising within cell populations over time. In stressed cells, entosis occurs as part of mixed-cell death response with parallel induction of apoptosis and necrosis, and we find that inhibition of one form of cell death leads to increased rates of another. Together, these findings identify stress-activated kinase signaling as a new inducer of entosis and demonstrate cross talk between different forms of cell death that can occur in parallel in response to UV radiation.

Transient phases of OXPHOS inhibitor resistance reveal underlying metabolic heterogeneity in single cells.

Cell-to-cell heterogeneity in metabolism plays an unknown role in physiology and pharmacology. To functionally characterize cellular variability in metabolism, we treated cells with inhibitors of oxidative phosphorylation (OXPHOS) and monitored their responses with live-cell reporters for ATP, ADP/ATP, or activity of the energy-sensing kinase AMPK. Across multiple OXPHOS inhibitors and cell types, we identified a subpopulation of cells resistant to activation of AMPK and reduction of ADP/ATP ratio. This resistant state persists transiently for at least several hours and can be inherited during cell divisions. OXPHOS inhibition suppresses the mTORC1 and ERK growth signaling pathways in sensitive cells, but not in resistant cells. Resistance is linked to a multi-factorial combination of increased glucose uptake, reduced protein biosynthesis, and G0/G1 cell-cycle status. Our results reveal dynamic fluctuations in cellular energetic balance and provide a basis for measuring and predicting the distribution of cellular responses to OXPHOS inhibition.

Combining Microbial Culturing With Mathematical Modeling in an Introductory Course-Based Undergraduate Research Experience.

Quantitative techniques are a critical part of contemporary biology research, but students interested in biology enter college with widely varying quantitative skills and attitudes toward mathematics. Course-based undergraduate research experiences (CUREs) may be an early way to build student competency and positive attitudes. Here we describe the design, implementation, and assessment of an introductory quantitative CURE focused on halophilic microbes. In this CURE, students culture and isolate halophilic microbes from environmental and food samples, perform growth assays, then use mathematical modeling to quantify the growth rate of strains in different salinities. To assess how the course may impact students' future academic plans and attitudes toward the use of math in biology, we used pre- and post-quarter surveys. Students who completed the course showed more positive attitudes toward science learning and an increased interest in pursuing additional quantitative biology experiences. We argue that the classroom application of microbiology methods, combined with mathematical modeling using student-generated data, provides a degree of student ownership, collaboration, iteration, and discovery that makes quantitative learning both relevant and exciting to students.

Oncogenic mutant RAS signaling activity is rescaled by the ERK/MAPK pathway.

Activating mutations in RAS are present in ~ 30% of human tumors, and the resulting aberrations in ERK/MAPK signaling play a central role in oncogenesis. However, the form of these signaling changes is uncertain, with activating RAS mutants linked to both increased and decreased ERK activation in vivo. Rationally targeting the kinase activity of this pathway requires clarification of the quantitative effects of RAS mutations. Here, we use live-cell imaging in cells expressing only one RAS isoform to quantify ERK activity with a new level of accuracy. We find that despite large differences in their biochemical activity, mutant KRAS isoforms within cells have similar ranges of ERK output. We identify roles for pathway-level effects, including variation in feedback strength and feedforward modulation of phosphatase activity, that act to rescale pathway sensitivity, ultimately resisting changes in the dynamic range of ERK activity while preserving responsiveness to growth factor stimuli. Our results reconcile seemingly inconsistent reports within the literature and imply that the signaling changes induced by RAS mutations early in oncogenesis are subtle.

Systems-Level Properties of EGFR-RAS-ERK Signaling Amplify Local Signals to Generate Dynamic Gene Expression Heterogeneity.

Intratumoral heterogeneity is associated with aggressive tumor behavior, therapy resistance, and poor patient outcomes. Such heterogeneity is thought to be dynamic, shifting over periods of minutes to hours in response to signaling inputs from the tumor microenvironment. However, models of this process have been inferred from indirect or post-hoc measurements of cell state, leaving the temporal details of signaling-driven heterogeneity undefined. Here, we developed a live-cell model system in which microenvironment-driven signaling dynamics can be directly observed and linked to variation in gene expression. Our analysis reveals that paracrine signaling between two cell types is sufficient to drive continual diversification of gene expression programs. This diversification emerges from systems-level properties of the EGFR-RAS-ERK signaling cascade, including intracellular amplification of amphiregulin-mediated paracrine signals and differential kinetic filtering by target genes including Fra-1, c-Myc, and Egr1. Our data enable more precise modeling of paracrine-driven transcriptional variation as a generator of gene expression heterogeneity. A record of this paper's transparent peer review process is included in the Supplemental Information.