Cellular heterogeneity in TNF/TNFR1 signalling: live cell imaging of cell fate decisions in single cells.

Cellular responses to TNF are inherently heterogeneous within an isogenic cell population and across different cell types. TNF promotes cell survival by activating pro-inflammatory NF-κB and MAPK signalling pathways but may also trigger apoptosis and necroptosis. Following TNF stimulation, the fate of individual cells is governed by the balance of pro-survival and pro-apoptotic signalling pathways. To elucidate the molecular mechanisms driving heterogenous responses to TNF, quantifying TNF/TNFR1 signalling at the single-cell level is crucial. Fluorescence live-cell imaging techniques offer real-time, dynamic insights into molecular processes in single cells, allowing for detection of rapid and transient changes, as well as identification of subpopulations, that are likely to be missed with traditional endpoint assays. Whilst fluorescence live-cell imaging has been employed extensively to investigate TNF-induced inflammation and TNF-induced cell death, it has been underutilised in studying the role of TNF/TNFR1 signalling pathway crosstalk in guiding cell-fate decisions in single cells. Here, we outline the various opportunities for pathway crosstalk during TNF/TNFR1 signalling and how these interactions may govern heterogenous responses to TNF. We also advocate for the use of live-cell imaging techniques to elucidate the molecular processes driving cell-to-cell variability in single cells. Understanding and overcoming cellular heterogeneity in response to TNF and modulators of the TNF/TNFR1 signalling pathway could lead to the development of targeted therapies for various diseases associated with aberrant TNF/TNFR1 signalling, such as rheumatoid arthritis, metabolic syndrome, and cancer.

Defining the condensate landscape of fusion oncoproteins.

Fusion oncoproteins (FOs) arise from chromosomal translocations in ~17% of cancers and are often oncogenic drivers. Although some FOs can promote oncogenesis by undergoing liquid-liquid phase separation (LLPS) to form aberrant biomolecular condensates, the generality of this phenomenon is unknown. We explored this question by testing 166 FOs in HeLa cells and found that 58% formed condensates. The condensate-forming FOs displayed physicochemical features distinct from those of condensate-negative FOs and segregated into distinct feature-based groups that aligned with their sub-cellular localization and biological function. Using Machine Learning, we developed a predictor of FO condensation behavior, and discovered that 67% of ~3000 additional FOs likely form condensates, with 35% of those predicted to function by altering gene expression. 47% of the predicted condensate-negative FOs were associated with cell signaling functions, suggesting a functional dichotomy between condensate-positive and -negative FOs. Our Datasets and reagents are rich resources to interrogate FO condensation in the future.

Overexpression of IκB⍺ modulates NF-κB activation of inflammatory target gene expression.

Cells respond to inflammatory stimuli such as cytokines by activation of the nuclear factor-κB (NF-κB) signalling pathway, resulting in oscillatory translocation of the transcription factor p65 between nucleus and cytoplasm in some cell types. We investigate the relationship between p65 and inhibitor-κB⍺ (IκBα) protein levels and dynamic properties of the system, and how this interaction impacts on the expression of key inflammatory genes. Using bacterial artificial chromosomes, we developed new cell models of IκB⍺-eGFP protein overexpression in a pseudo-native genomic context. We find that cells with high levels of the negative regulator IκBα remain responsive to inflammatory stimuli and maintain dynamics for both p65 and IκBα. In contrast, canonical target gene expression is dramatically reduced by overexpression of IκBα, but can be partially rescued by overexpression of p65. Treatment with leptomycin B to promote nuclear accumulation of IκB⍺ also suppresses canonical target gene expression, suggesting a mechanism in which nuclear IκB⍺ accumulation prevents productive p65 interaction with promoter binding sites. This causes reduced target promoter binding and gene transcription, which we validate by chromatin immunoprecipitation and in primary cells. Overall, we show how inflammatory gene transcription is modulated by the expression levels of both IκB⍺ and p65. This results in an anti-inflammatory effect on transcription, demonstrating a broad mechanism to modulate the strength of inflammatory response.

Apoptotic priming is defined by the dynamic exchange of Bcl-2 proteins between mitochondria and cytosol.

Apoptosis is regulated by interactions between the BH3-only and multi-domain Bcl-2 family proteins. These interactions are integrated on the outer mitochondrial membrane (OMM) where they set the threshold for apoptosis, known as mitochondrial priming. However, how mitochondrial priming is controlled at the level of single cells remains unclear. Retrotranslocation of Bcl-XL has been proposed as one mechanism, removing pro-apoptotic Bcl-2 proteins from the OMM, thus reducing priming. Contrary to this view, we now show that Bcl-XL retrotranslocation is inhibited by binding to its BH3-only partners, resulting in accumulation of these protein complexes on mitochondria. We find that Bcl-XL retrotranslocation dynamics are tightly coupled to mitochondrial priming. Quantifying these dynamics indicates the heterogeneity in priming between cells within a population and predicts how they subsequently respond to a pro-apoptotic signal.

Phase Separation Mediates NUP98 Fusion Oncoprotein Leukemic Transformation.

NUP98 fusion oncoproteins (FO) are drivers in pediatric leukemias and many transform hematopoietic cells. Most NUP98 FOs harbor an intrinsically disordered region from NUP98 that is prone to liquid-liquid phase separation (LLPS) in vitro. A predominant class of NUP98 FOs, including NUP98-HOXA9 (NHA9), retains a DNA-binding homeodomain, whereas others harbor other types of DNA- or chromatin-binding domains. NUP98 FOs have long been known to form puncta, but long-standing questions are how nuclear puncta form and how they drive leukemogenesis. Here we studied NHA9 condensates and show that homotypic interactions and different types of heterotypic interactions are required to form nuclear puncta, which are associated with aberrant transcriptional activity and transformation of hematopoietic stem and progenitor cells. We also show that three additional leukemia-associated NUP98 FOs (NUP98-PRRX1, NUP98-KDM5A, and NUP98-LNP1) form nuclear puncta and transform hematopoietic cells. These findings indicate that LLPS is critical for leukemogenesis by NUP98 FOs. SIGNIFICANCE: We show that homotypic and heterotypic mechanisms of LLPS control NUP98-HOXA9 puncta formation, modulating transcriptional activity and transforming hematopoietic cells. Importantly, these mechanisms are generalizable to other NUP98 FOs that share similar domain structures. These findings address long-standing questions on how nuclear puncta form and their link to leukemogenesis. This article is highlighted in the In This Issue feature, p. 873.

Natural killer cell immune synapse formation and cytotoxicity are controlled by tension of the target interface.

Natural killer (NK) cells can kill infected or transformed cells via a lytic immune synapse. Diseased cells may exhibit altered mechanical properties but how this impacts NK cell responsiveness is unknown. We report that human NK cells were stimulated more effectively to secrete granzymes A and B, FasL (also known as FasLG), granulysin and IFNγ, by stiff (142 kPa) compared to soft (1 kPa) planar substrates. To create surrogate spherical targets of defined stiffness, sodium alginate was used to synthesise soft (9 kPa), medium (34 kPa) or stiff (254 kPa) cell-sized beads, coated with antibodies against activating receptor NKp30 (also known as NCR3) and the integrin LFA-1 (also known as ITGAL). Against stiff beads, NK cells showed increased degranulation. Polarisation of the microtubule-organising centre and lytic granules were impaired against soft targets, which instead resulted in the formation of unstable kinapses. Thus, by varying target stiffness to characterise the mechanosensitivity of immune synapses, we identify soft targets as a blind spot in NK cell recognition. This article has an associated First Person interview with the co-first authors of the paper.

Transcription Factor Pit-1 Affects Transcriptional Timing in the Dual-Promoter Human Prolactin Gene.

Gene transcription occurs in short bursts interspersed with silent periods, and these kinetics can be altered by promoter structure. The effect of alternate promoter architecture on transcription bursting is not known. We studied the human prolactin (hPRL) gene that contains 2 promoters, a pituitary-specific promoter that requires the transcription factor Pit-1 and displays dramatic transcriptional bursting activity and an alternate upstream promoter that is active in nonpituitary tissues. We studied large hPRL genomic fragments with luciferase reporters, and used bacterial artificial chromosome recombineering to manipulate critical promoter regions. Stochastic switch mathematical modelling of single-cell time-lapse luminescence image data revealed that the Pit-1-dependent promoter showed longer, higher-amplitude transcriptional bursts. Knockdown studies confirmed that the presence of Pit-1 stabilized and prolonged periods of active transcription. Pit-1 therefore plays an active role in establishing the timing of transcription cycles, in addition to its cell-specific functions.

Calcium dynamics and chromatin remodelling underlie heterogeneity in prolactin transcription.

Pituitary cells have been reported to show spontaneous calcium oscillations and dynamic transcription cycles. To study both processes in the same living cell in real time, we used rat pituitary GH3 cells stably expressing human prolactin-luciferase or prolactin-EGFP reporter gene constructs loaded with a fluorescent calcium indicator and measured activity using single-cell time-lapse microscopy. We observed heterogeneity between clonal cells in the calcium activity and prolactin transcription in unstimulated conditions. There was a significant correlation between cells displaying spontaneous calcium spikes and cells showing spontaneous bursts in prolactin expression. Notably, cells showing no basal calcium activity showed low prolactin expression but elicited a significantly greater transcriptional response to BayK8644 compared to cells showing basal calcium activity. This suggested the presence of two subsets of cells within the population at any one time. Fluorescence-activated cell sorting was used to sort cells into two populations based on the expression level of prolactin-EGFP however, the bimodal pattern of expression was restored within 26 h. Chromatin immunoprecipitation showed that these sorted populations were distinct due to the extent of histone acetylation. We suggest that maintenance of a heterogeneous bimodal population is a fundamental characteristic of this cell type and that calcium activation and histone acetylation, at least in part, drive prolactin transcriptional competence.

Multiplexing information flow through dynamic signalling systems.

We consider how a signalling system can act as an information hub by multiplexing information arising from multiple signals. We formally define multiplexing, mathematically characterise which systems can multiplex and how well they can do it. While the results of this paper are theoretical, to motivate the idea of multiplexing, we provide experimental evidence that tentatively suggests that the NF-κB transcription factor can multiplex information about changes in multiple signals. We believe that our theoretical results may resolve the apparent paradox of how a system like NF-κB that regulates cell fate and inflammatory signalling in response to diverse stimuli can appear to have the low information carrying capacity suggested by recent studies on scalar signals. In carrying out our study, we introduce new methods for the analysis of large, nonlinear stochastic dynamic models, and develop computational algorithms that facilitate the calculation of fundamental constructs of information theory such as Kullback-Leibler divergences and sensitivity matrices, and link these methods to a new theory about multiplexing information. We show that many current models such as those of the NF-κB system cannot multiplex effectively and provide models that overcome this limitation using post-transcriptional modifications.

ER stress-linked autophagy stabilizes apoptosis effector PERP and triggers its co-localization with SERCA2b at ER-plasma membrane junctions.

Specific molecular interactions that underpin the switch between ER stress-triggered autophagy-mediated cellular repair and cellular death by apoptosis are not characterized. This study reports the unexpected interaction elicited by ER stress between the plasma membrane (PM)-localized apoptosis effector PERP and the ER Ca2+ pump SERCA2b. We show that the p53 effector PERP, which specifically induces apoptosis when expressed above a threshold level, has a heterogeneous distribution across the PM of un-stressed cells and is actively turned over by the lysosome. PERP is upregulated following sustained starvation-induced autophagy, which precedes the onset of apoptosis indicating that PERP protein levels are controlled by a lysosomal pathway that is sensitive to cellular physiological state. Furthermore, ER stress stabilizes PERP at the PM and induces its increasing co-localization with SERCA2b at ER-PM junctions. The findings highlight a novel crosstalk between pro-survival autophagy and pro-death apoptosis pathways and identify, for the first time, accumulation of an apoptosis effector to ER-PM junctions in response to ER stress.

Correction: Genotype-by-environment interactions affecting heterosis in maize.

[This corrects the article DOI: 10.1371/journal.pone.0191321.].

Disentangling juxtacrine from paracrine signalling in dynamic tissue.

Prolactin is a major hormone product of the pituitary gland, the central endocrine regulator. Despite its physiological importance, the cell-level mechanisms of prolactin production are not well understood. Having significantly improved the resolution of real-time-single-cell-GFP-imaging, the authors recently revealed that prolactin gene transcription is highly dynamic and stochastic yet shows space-time coordination in an intact tissue slice. However, it still remains an open question as to what kind of cellular communication mediates the observed space-time organization. To determine the type of interaction between cells we developed a statistical model. The degree of similarity between two expression time series was studied in terms of two distance measures, Euclidean and geodesic, the latter being a network-theoretic distance defined to be the minimal number of edges between nodes, and this was used to discriminate between juxtacrine from paracrine signalling. The analysis presented here suggests that juxtacrine signalling dominates. To further determine whether the coupling is coordinating transcription or post-transcriptional activities we used stochastic switch modelling to infer the transcriptional profiles of cells and estimated their similarity measures to deduce that their spatial cellular coordination involves coupling of transcription via juxtacrine signalling. We developed a computational model that involves an inter-cell juxtacrine coupling, yielding simulation results that show space-time coordination in the transcription level that is in agreement with the above analysis. The developed model is expected to serve as the prototype for the further study of tissue-level organised gene expression for epigenetically regulated genes, such as prolactin.

C9orf72 Poly(PR) Dipeptide Repeats Disturb Biomolecular Phase Separation and Disrupt Nucleolar Function.

Repeat expansion in the C9orf72 gene is the most common cause of the neurodegenerative disorder amyotrophic lateral sclerosis (C9-ALS) and is linked to the unconventional translation of five dipeptide-repeat polypeptides (DPRs). The two enriched in arginine, poly(GR) and poly(PR), infiltrate liquid-like nucleoli, co-localize with the nucleolar protein nucleophosmin (NPM1), and alter the phase separation behavior of NPM1 in vitro. Here, we show that poly(PR) DPRs bind tightly to a long acidic tract within the intrinsically disordered region of NPM1, altering its phase separation with nucleolar partners to the extreme of forming large, soluble complexes that cause droplet dissolution in vitro. In cells, poly(PR) DPRs disperse NPM1 from nucleoli and entrap rRNA in static condensates in a DPR-length-dependent manner. We propose that R-rich DPR toxicity involves disrupting the role of phase separation by NPM1 in organizing ribosomal proteins and RNAs within the nucleolus.

Selection Signatures Underlying Dramatic Male Inflorescence Transformation During Modern Hybrid Maize Breeding.

Inflorescence capacity plays a crucial role in reproductive fitness in plants, and in production of hybrid crops. Maize is a monoecious species bearing separate male and female flowers (tassel and ear, respectively). The switch from open-pollinated populations of maize to hybrid-based breeding schemes in the early 20th century was accompanied by a dramatic reduction in tassel size, and the trend has continued with modern breeding over the recent decades. The goal of this study was to identify selection signatures in genes that may underlie this dramatic transformation. Using a population of 942 diverse inbred maize accessions and a nested association mapping population comprising three 200-line biparental populations, we measured 15 tassel morphological characteristics by manual and image-based methods. Genome-wide association studies identified 242 single nucleotide polymorphisms significantly associated with measured traits. We compared 41 unselected lines from the Iowa Stiff Stalk Synthetic (BSSS) population to 21 highly selected lines developed by modern commercial breeding programs, and found that tassel size and weight were reduced significantly. We assayed genetic differences between the two groups using three selection statistics: cross population extended haplotype homozogysity, cross-population composite likelihood ratio, and fixation index. All three statistics show evidence of selection at genomic regions associated with tassel morphology relative to genome-wide null distributions. These results support the tremendous effect, both phenotypic and genotypic, that selection has had on maize male inflorescence morphology.

Methods for Physical Characterization of Phase-Separated Bodies and Membrane-less Organelles.

Membrane-less organelles are cellular structures which arise through the phenomenon of phase separation. This process enables compartmentalization of specific sets of macromolecules (e.g., proteins, nucleic acids), thereby regulating cellular processes by increasing local concentration, and modulating the structure and dynamics of their constituents. Understanding the connection between structure, material properties and function of membrane-less organelles requires inter-disciplinary approaches, which address length and timescales that span several orders of magnitude (e.g., Ångstroms to micrometer, picoseconds to hours). In this review, we discuss the wide variety of methods that have been applied to characterize the morphology, rheology, structure and dynamics of membrane-less organelles and their components, in vitro and in live cells.

Quantitative analysis of competitive cytokine signaling predicts tissue thresholds for the propagation of macrophage activation.

Toll-like receptor (TLR) signaling regulates macrophage activation and effector cytokine propagation in the constrained environment of a tissue. In macrophage populations, TLR4 stimulates the dose-dependent transcription of nuclear factor κB (NF-κB) target genes. However, using single-RNA counting, we found that individual cells exhibited a wide range (three orders of magnitude) of expression of the gene encoding the proinflammatory cytokine tumor necrosis factor-α (TNF-α). The TLR4-induced TNFA transcriptional response correlated with the extent of NF-κB signaling in the cells and their size. We compared the rates of TNF-α production and uptake in macrophages and mouse embryonic fibroblasts and generated a mathematical model to explore the heterogeneity in the response of macrophages to TLR4 stimulation and the propagation of the TNF-α signal in the tissue. The model predicts that the local propagation of the TLR4-dependent TNF-α response and cellular NF-κB signaling are limited to small distances of a few cell diameters between neighboring tissue-resident macrophages. In our predictive model, TNF-α propagation was constrained by competitive uptake of TNF-α from the environment, rather than by heterogeneous production of the cytokine. We propose that the highly constrained architecture of tissues enables effective localized propagation of inflammatory cues while avoiding out-of-context responses at longer distances.

Quantitative analysis reveals crosstalk mechanisms of heat shock-induced attenuation of NF-κB signaling at the single cell level.

Elevated temperature induces the heat shock (HS) response, which modulates cell proliferation, apoptosis, the immune and inflammatory responses. However, specific mechanisms linking the HS response pathways to major cellular signaling systems are not fully understood. Here we used integrated computational and experimental approaches to quantitatively analyze the crosstalk mechanisms between the HS-response and a master regulator of inflammation, cell proliferation, and apoptosis the Nuclear Factor κB (NF-κB) system. We found that populations of human osteosarcoma cells, exposed to a clinically relevant 43°C HS had an attenuated NF-κB p65 response to Tumor Necrosis Factor α (TNFα) treatment. The degree of inhibition of the NF-κB response depended on the HS exposure time. Mathematical modeling of single cells indicated that individual crosstalk mechanisms differentially encode HS-mediated NF-κB responses while being consistent with the observed population-level responses. In particular "all-or-nothing" encoding mechanisms were involved in the HS-dependent regulation of the IKK activity and IκBα phosphorylation, while others involving transport were "analogue". In order to discriminate between these mechanisms, we used live-cell imaging of nuclear translocations of the NF-κB p65 subunit. The single cell responses exhibited "all-or-nothing" encoding. While most cells did not respond to TNFα stimulation after a 60 min HS, 27% showed responses similar to those not receiving HS. We further demonstrated experimentally and theoretically that the predicted inhibition of IKK activity was consistent with the observed HS-dependent depletion of the IKKα and IKKß subunits in whole cell lysates. However, a combination of "all-or-nothing" crosstalk mechanisms was required to completely recapitulate the single cell data. We postulate therefore that the heterogeneity of the single cell responses might be explained by the cell-intrinsic variability of HS-modulated IKK signaling. In summary, we show that high temperature modulates NF-κB responses in single cells in a complex and unintuitive manner, which needs to be considered in hyperthermia-based treatment strategies.

Genotype-by-environment interactions affecting heterosis in maize.

The environment can influence heterosis, the phenomena in which the offspring of two inbred parents exhibits phenotypic performance beyond the inbred parents for specific traits. In this study we measured 25 traits in a set of 47 maize hybrids and their inbred parents grown in 16 different environments with varying levels of average productivity. By quantifying 25 vegetative and reproductive traits across the life cycle we were able to analyze interactions between the environment and multiple distinct instances of heterosis. The magnitude and rank among hybrids for better-parent heterosis (BPH) varied for the different traits and environments. Across the traits, a higher within plot variance was observed for inbred lines compared to hybrids. However, for most traits, variance across environments was not significantly different for inbred lines compared to hybrids. Further, for many traits the correlations of BPH to hybrid performance and BPH to better parent performance were of comparable magnitude. These results indicate that inbred lines and hybrids show similar trends in environmental response and both are contributing to observed genotype-by-environment interactions for heterosis. This study highlights the degree of heterosis is not an inherent trait of a specific hybrid, but varies depending on the trait measured and the environment where that trait is measured. Studies that attempt to correlate molecular processes with heterosis are hindered by the fact that heterosis is not a consistent attribute of a specific hybrid.

Live Imaging of Label-Free Graphene Oxide Reveals Critical Factors Causing Oxidative-Stress-Mediated Cellular Responses.

The interest in graphene and its translation into commercial products has been expanding at a high pace. Based on previously described pulmonary safety concerns for carbon nanomaterials, there is a great need to define parameters guiding interactions between graphene-based materials and the pulmonary system. The aim of the present study was to determine the importance of two critical parameters: lateral dimensions of the material and coating with proteins in relation to each other and their pulmonary impact. Endotoxin-free materials with distinct lateral dimensions, s-GO (50-200 nm) and l-GO (5-15 µm), were produced and thoroughly characterized. Exploiting intrinsic fluorescence of graphene oxide (GO) and using confocal live-cell imaging, the behavior of the cells in response to the material was visualized in real time. Although BEAS-2B cells internalized GO efficiently, l-GO was linked to higher plasma membrane interactions correlated with elevated reactive oxygen species (ROS) levels, pro-inflammatory response, and greater cytotoxicity, in agreement with the oxidative stress paradigm. For both GO types, the presence of serum alleviated lipid peroxidation of plasma membrane and decreased intracellular ROS levels. However, protein coating was not enough to entirely mitigate toxicity and inflammatory response induced by l-GO. In vitro results were validated in vivo, as l-GO was more prone to induce pulmonary granulomatous response in mice compared to s-GO. In conclusion, the lateral dimension of GO played a more important role than serum protein coating in determining biological responses to the material. It was also demonstrated that time-lapse imaging of live cells interacting with label-free GO sheets can be used as a tool to assess GO-induced cytotoxicity.

Asymmetry between Activation and Deactivation during a Transcriptional Pulse.

Transcription in eukaryotic cells occurs in gene-specific bursts or pulses of activity. Recent studies identified a spectrum of transcriptionally active "on-states," interspersed with periods of inactivity, but these "off-states" and the process of transcriptional deactivation are poorly understood. To examine what occurs during deactivation, we investigate the dynamics of switching between variable rates. We measured live single-cell expression of luciferase reporters from human growth hormone or human prolactin promoters in a pituitary cell line. Subsequently, we applied a statistical variable-rate model of transcription, validated by single-molecule FISH, to estimate switching between transcriptional rates. Under the assumption that transcription can switch to any rate at any time, we found that transcriptional activation occurs predominantly as a single switch, whereas deactivation occurs with graded, stepwise decreases in transcription rate. Experimentally altering cAMP signalling with forskolin or chromatin remodelling with histone deacetylase inhibitor modifies the duration of defined transcriptional states. Our findings reveal transcriptional activation and deactivation as mechanistically independent, asymmetrical processes.