Search

In mammals, X chromosome dosage is balanced between sexes through the silencing of one X chromosome in females. Recent single-cell RNA sequencing analysis demonstrated that the inactivation of the X chromosome is accompanied by the upregulation of the active X chromosome (Xa) during mouse embryogenesis. Here, we have investigated if the reactivation of inactive-X (Xi) leads to the loss of Xa upregulation in different cellular or developmental contexts. We find that while Xi reactivation and loss of Xa upregulation are tightly coupled in mouse embryonic epiblast and induced pluripotent stem cells, that is not the case in germ cells. Moreover, we demonstrate that partial reactivation of Xi in mouse extra-embryonic endoderm stem cells and human B cells does not result in the loss of Xa upregulation. Finally, we have established a mathematical model for the transcriptional coordination of two X chromosomes. Together, we conclude that the reactivation of Xi is not always synchronized with the loss of Xa upregulation., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

A common feature of bacterial, fungal and cancer cell populations upon treatment is the presence of tolerant and persistent cells able to survive, and sometimes grow, even in the presence of usually inhibitory or lethal drug concentrations, driven by non-genetic differences among individual cells in a population. Here we review and compare data obtained on drug survival in bacteria, fungi and cancer cells to unravel common characteristics and cellular pathways, and to point their singularities. This comparative work also allows to cross-fertilize ideas across fields. We particularly focus on the role of gene expression variability in the emergence of cell-cell non-genetic heterogeneity because it represents a possible common basic molecular process at the origin of most persistence phenomena and could be monitored and tuned to help improve therapeutic interventions., Competing Interests: Declaration of competing interest The authors declare no competing interests., (Copyright © 2024. Published by Elsevier Inc.)

Planar cell polarity (PCP) - tissue-scale alignment of the direction of asymmetric localization of proteins at the cell-cell interface - is essential for embryonic development and physiological functions. Abnormalities in PCP can result in developmental imperfections, including neural tube closure defects and misaligned hair follicles. Decoding the mechanisms responsible for PCP establishment and maintenance remains a fundamental open question. While the roles of various molecules - broadly classified into "global" and "local" modules - have been well-studied, their necessity and sufficiency in explaining PCP and connecting their perturbations to experimentally observed patterns have not been examined. Here, we develop a minimal model that captures the proposed features of PCP establishment - a global tissue-level gradient and local asymmetric distribution of protein complexes. The proposed model suggests that while polarity can emerge without a gradient, the gradient not only acts as a global cue but also increases the robustness of PCP against stochastic perturbations. We also recapitulated and quantified the experimentally observed features of swirling patterns and domineering non-autonomy, using only three free model parameters - the rate of protein binding to membrane, the concentration of PCP proteins, and the gradient steepness. We explain how self-stabilizing asymmetric protein localizations in the presence of tissue-level gradient can lead to robust PCP patterns and reveal minimal design principles for a polarized system., Competing Interests: DS, SR, MJ, MR The authors declare that no competing interests exist., (© 2024, Singh et al.)

Elucidating the emergent dynamics of cellular differentiation networks is crucial to understanding cell-fate decisions. Toggle switch - a network of mutually repressive lineage-specific transcription factors A and B - enables two phenotypes from a common progenitor: (high A, low B) and (low A, high B). However, the dynamics of networks enabling differentiation of more than two phenotypes from a progenitor cell has not been well-studied. Here, we investigate the dynamics of a four-node network A, B, C, and D inhibiting each other, forming a toggle tetrahedron. Our simulations show that this network is multistable and predominantly allows for the co-existence of six hybrid phenotypes where two of the nodes are expressed relatively high as compared to the remaining two, for instance (high A, high B, low C, low D). Finally, we apply our results to understand naïve CD4 + T cell differentiation into Th1, Th2, Th17 and Treg subsets, suggesting Th1/Th2/Th17/Treg decision-making to be a two-step process., (© 2024. The Author(s).)

Metastasis is the leading cause of cancer-related deaths, yet its regulatory mechanisms are not fully understood. Small-cell lung cancer (SCLC) is the most metastatic form of lung cancer, with most patients presenting with widespread disease, making it an ideal model for studying metastasis. However, the lack of suitable preclinical models has limited such studies. We utilized well-annotated rapid autopsy-derived tumors to develop xenograft models that mimic key features of SCLC, including histopathology, rapid and widespread development of metastasis to the liver, brain, adrenal, bone marrow, and kidneys within weeks, and response to chemotherapy. By integrating in vivo lineage selection with comprehensive transcriptomic and epigenomic analyses, we identified critical cellular programs driving metastatic organotropism to the liver and brain, the most common sites of SCLC metastasis. Our findings reveal the key role of nuclear-cytoskeletal interactions in SCLC liver metastasis. Specifically, the loss of the nuclear envelope protein lamin A/C, encoded by the LMNA gene, increased nuclear deformability and significantly increased the incidence of liver metastasis. Human liver metastases exhibited reduced LMNA expression compared to other metastatic sites, correlating with poorer patient outcomes and increased mortality. This study introduces novel preclinical models for SCLC metastasis and highlights pathways critical for organ-specific metastasis, offering new avenues for the development of targeted therapies to prevent or treat metastatic disease., Competing Interests: Conflict of Interest disclosure: A.T. received grants to NCI from EMD Serono Research & Development, AstraZeneca, Gilead Sciences, and ProLynx during the conduct of the study.

Neuroblastoma is the most frequent extracranial pediatric tumor and leads to 15% of all cancer-related deaths in children. Tumor relapse and therapy resistance in neuroblastoma are driven by phenotypic plasticity and heterogeneity between noradrenergic (NOR) and mesenchymal (MES) cell states. Despite the importance of this phenotypic plasticity, the dynamics and molecular patterns associated with these bidirectional cell-state transitions remain relatively poorly understood. Here, we analyze multiple RNA-seq datasets at both bulk and single-cell resolution, to understand the association between NOR- and MES-specific factors. We observed that NOR-specific and MES-specific expression patterns are largely mutually exclusive, exhibiting a "teams-like" behavior among the genes involved, reminiscent of our earlier observations in lung cancer and melanoma. This antagonism between NOR and MES phenotypes was also associated with metabolic reprogramming and with immunotherapy targets PD-L1 and GD2 as well as with experimental perturbations driving the NOR-MES and/or MES-NOR transition. Further, these "teams-like" patterns were seen only among the NOR- and MES-specific genes, but not in housekeeping genes, possibly highlighting a hallmark of network topology enabling cancer cell plasticity.

Notch-Delta-Jagged (NDJ) signaling among neighboring cells contributes crucially to spatiotemporal pattern formation and developmental decision-making. Despite numerous detailed mathematical models, their high-dimensionality parametric space limits analytical treatment, especially regarding local microenvironmental fluctuations. Using the low-dimensional dynamics of the recently postulated least microenvironmental uncertainty principle (LEUP) framework, we showcase how the LEUP formalism recapitulates a noisy NDJ spatial patterning. Our LEUP simulations show that local phenotypic entropy increases for lateral inhibition but decreases for lateral induction. This distinction allows us to identify a critical parameter that captures the transition from a Notch-Delta-driven lateral inhibition to a Notch-Jagged-driven lateral induction phenomenon and suggests random phenotypic patterning in the case of lack of dominance of either Notch-Delta or Notch-Jagged signaling. Our results enable an analytical treatment to map the high-dimensional dynamics of NDJ signaling on tissue-level patterning and can possibly be generalized to decode operating principles of collective cellular decision-making., Competing Interests: The authors declared no competing interests., (© 2024 The Author(s).)

The ultraviolet (UV) radiation triggers a pigmentation response in human skin, wherein, melanocytes rapidly activate divergent maturation and proliferation programs. Using single-cell sequencing, we demonstrate that these 2 programs are segregated in distinct subpopulations in melanocytes of human and zebrafish skin. The coexistence of these 2 cell states in cultured melanocytes suggests possible cell autonomy. Luria-Delbrück fluctuation test reveals that the initial establishment of these states is stochastic. Tracking of pigmenting cells ascertains that the stochastically acquired state is faithfully propagated in the progeny. A systemic approach combining single-cell multi-omics (RNA+ATAC) coupled to enhancer mapping with H3K27 acetylation successfully identified state-specific transcriptional networks. This comprehensive analysis led to the construction of a gene regulatory network (GRN) that under the influence of noise, establishes a bistable system of pigmentation and proliferation at the population level. This GRN recapitulates melanocyte behaviour in response to external cues that reinforce either of the states. Our work highlights that inherent stochasticity within melanocytes establishes dedicated states, and the mature state is sustained by selective enhancers mark through histone acetylation. While the initial cue triggers a proliferation response, the continued signal activates and maintains the pigmenting subpopulation via epigenetic imprinting. Thereby our study provides the basis of coexistence of distinct populations which ensures effective pigmentation response while preserving the self-renewal capacity., Competing Interests: RSG is the co-founder of Vyome Biosciences Pvt Ltd., a biopharmaceutical company working in the area of dermatology. Other authors declare no competing interest., (Copyright: © 2024 Aggarwal et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Exploring the complexity of the epithelial-to-mesenchymal transition (EMT) unveils a diversity of potential cell fates; however, the exact timing and mechanisms by which early cell states diverge into distinct EMT trajectories remain unclear. Studying these EMT trajectories through single-cell RNA sequencing is challenging due to the necessity of sacrificing cells for each measurement. In this study, we employed optimal-transport analysis to reconstruct the past trajectories of different cell fates during TGF-beta-induced EMT in the MCF10A cell line. Our analysis revealed three distinct trajectories leading to low EMT, partial EMT, and high EMT states. Cells along the partial EMT trajectory showed substantial variations in the EMT signature and exhibited pronounced stemness. Throughout this EMT trajectory, we observed a consistent downregulation of the EED and EZH2 genes. This finding was validated by recent inhibitor screens of EMT regulators and CRISPR screen studies. Moreover, we applied our analysis of early-phase differential gene expression to gene sets associated with stemness and proliferation, pinpointing ITGB4 , LAMA3 , and LAMB3 as genes differentially expressed in the initial stages of the partial versus high EMT trajectories. We also found that CENPF , CKS1B , and MKI67 showed significant upregulation in the high EMT trajectory. While the first group of genes aligns with findings from previous studies, our work uniquely pinpoints the precise timing of these upregulations. Finally, the identification of the latter group of genes sheds light on potential cell cycle targets for modulating EMT trajectories., Competing Interests: Competing interests statement:F.M. is a co-founder of and has equity in Harbinger Health, has equity in Zephyr AI, and serves as a consultant for both companies. She is also on the board of directors of Exscientia Plc. F.M. declares that none of these relationships are directly or indirectly related to the content of this manuscript. One of the authors (H.L.) was co-author on a review paper (in 2021) with one of the referees (Y.K.), and a different author (M.K.J.) was co-author on a different review paper (in 2023) with the other referee (Q.N.). All other authors declare no conflicts.

Background: The glucocorticoid receptor (GR) is frequently expressed in breast cancer (BC), and its prognostic implications are contingent on estrogen receptor (ER) status. To address conflicting reports and explore therapeutic potential, a GR signature (GRsig) independent of ER status was developed. We also investigated cell type-specific GR protein expression in BC tumor epithelial cells and infiltrating lymphocytes., Methods: GRsig was derived from Dexamethasone treated cell lines through a bioinformatic pipeline. Immunohistochemistry assessed GR protein expression. Associations between GRsig and tumor phenotypes (proliferation, cytolytic activity (CYT), immune cell distribution, and epithelial-to-mesenchymal transition (EMT) were explored in public datasets. Single-cell RNA sequencing data evaluated context-dependent GR roles, and a cell type-specific prognostic role was assessed in an independent BC cohort., Results: High GRsig levels were associated with a favorable prognosis across BC subtypes. Tumor-specific high GRsig correlated with lower proliferation, increased CYT, and anti-tumorigenic immune cells. Single-cell data analysis revealed higher GRsig expression in immune cells, negatively correlating with EMT while a positive correlation was observed with EMT primarily in tumor and stromal cells. Univariate and multivariate analyses demonstrated the robust and independent predictive capability of GRsig for favorable prognosis. GR protein expression on immune cells in triple-negative tumors indicated a favorable prognosis., Conclusion: This study underscores the cell type-specific role of GR, where its expression on tumor cells is associated with aggressive features like EMT, while in infiltrating lymphocytes, it predicts a better prognosis, particularly within TNBC tumors. The GRsig emerges as a promising independent prognostic indicator across diverse BC subtypes., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024. Published by Elsevier Inc.)

Cancer cell populations comprise phenotypes distributed among the epithelial-mesenchymal (E-M) spectrum. However, it remains unclear which population-level processes give rise to the observed experimental distribution and dynamical changes in E-M heterogeneity, including (1) differential growth, (2) cell-state switching, and (3) population density-dependent growth or state-transition rates. Here, we analyze the necessity of these three processes in explaining the dynamics of E-M population distributions as observed in PMC42-LA and HCC38 breast cancer cells. We find that, while cell-state transition is necessary to reproduce experimental observations of dynamical changes in E-M fractions, including density-dependent growth interactions (cooperation or suppression) better explains the data. Further, our models predict that treatment of HCC38 cells with transforming growth factor β (TGF-β) signaling and Janus kinase 2/signal transducer and activator of transcription 3 (JAK2/3) inhibitors enhances the rate of mesenchymal-epithelial transition (MET) instead of lowering that of E-M transition (EMT). Overall, our study identifies the population-level processes shaping the dynamics of spontaneous E-M heterogeneity in breast cancer cells., Competing Interests: The authors declare no competing interests., (© 2024 The Author(s).)

Small-cell lung cancer (SCLC) is the most fatal form of lung cancer. Intratumoral heterogeneity, marked by neuroendocrine (NE) and non-neuroendocrine (non-NE) cell states, defines SCLC, but the cell-extrinsic drivers of SCLC plasticity are poorly understood. To map the landscape of SCLC tumor microenvironment (TME), we apply spatially resolved transcriptomics and quantitative mass spectrometry-based proteomics to metastatic SCLC tumors obtained via rapid autopsy. The phenotype and overall composition of non-malignant cells in the TME exhibit substantial variability, closely mirroring the tumor phenotype, suggesting TME-driven reprogramming of NE cell states. We identify cancer-associated fibroblasts (CAFs) as a crucial element of SCLC TME heterogeneity, contributing to immune exclusion, and predicting exceptionally poor prognosis. Our work provides a comprehensive map of SCLC tumor and TME ecosystems, emphasizing their pivotal role in SCLC's adaptable nature, opening possibilities for reprogramming the TME-tumor communications that shape SCLC tumor states., Competing Interests: Declaration of interests A.T. received grants to NCI from EMD Serono Research & Development, AstraZeneca, Gilead Sciences, and ProLynx during the conduct of the study., (Published by Elsevier Inc.)

Collective cell invasion (CCI), a canon of most invasive solid tumors, is an emergent property of the interactions between cancer cells and their surrounding extracellular matrix (ECM). However, tumor populations invariably consist of cells expressing variable levels of adhesive proteins that mediate such interactions, disallowing an intuitive understanding of how tumor invasiveness at a multicellular scale is influenced by spatial heterogeneity of cell-cell and cell-ECM adhesion. Here, we have used a Cellular Potts model-based multiscale computational framework that is constructed on the histopathological principles of glandular cancers. In earlier efforts on homogenous cancer cell populations, this framework revealed the relative ranges of interactions, including cell-cell and cell-ECM adhesion that drove collective, dispersed, and mixed multimodal invasion. Here, we constitute a tumor core of two separate cell subsets showing distinct intra- and inter-subset cell-cell or cell-ECM adhesion strengths. These two subsets of cells are arranged to varying extents of spatial intermingling, which we call the heterogeneity index (HI). We observe that low and high inter-subset cell adhesion favors invasion of high-HI and low-HI intermingled populations with distinct intra-subset cell-cell adhesion strengths, respectively. In addition, for explored values of cell-ECM adhesion strengths, populations with high HI values collectively invade better than those with lower HI values. We then asked how spatial invasion is regulated by progressively intermingled cellular subsets that are epithelial, i.e., showed high cell-cell but poor cell-ECM adhesion, and mesenchymal, i.e., with reversed adhesion strengths to the former. Here too, inter-subset adhesion plays an important role in contextualizing the proportionate relationship between HI and invasion. An exception to this relationship is seen for cases of heterogeneous cell-ECM adhesion where sub-maximal HI patterns with higher outer localization of cells with stronger ECM adhesion collectively invade better than their relatively higher-HI counterparts. Our simulations also reveal how adhesion heterogeneity qualifies collective invasion, when either cell-cell or cell-ECM adhesion type is varied but results in an invasive dispersion when both adhesion types are simultaneously altered., Competing Interests: Declaration of interests The authors do not declare any conflicts of interest., (Copyright © 2024 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Intra-tumoral phenotypic heterogeneity promotes tumor relapse and therapeutic resistance and remains an unsolved clinical challenge. Decoding the interconnections among different biological axes of plasticity is crucial to understand the molecular origins of phenotypic heterogeneity. Here, we use multi-modal transcriptomic data-bulk, single-cell, and spatial transcriptomics-from breast cancer cell lines and primary tumor samples, to identify associations between epithelial-mesenchymal transition (EMT) and luminal-basal plasticity-two key processes that enable heterogeneity. We show that luminal breast cancer strongly associates with an epithelial cell state, but basal breast cancer is associated with hybrid epithelial/mesenchymal phenotype(s) and higher phenotypic heterogeneity. Mathematical modeling of core underlying gene regulatory networks representative of the crosstalk between the luminal-basal and epithelial-mesenchymal axes elucidate mechanistic underpinnings of the observed associations from transcriptomic data. Our systems-based approach integrating multi-modal data analysis with mechanism-based modeling offers a predictive framework to characterize intra-tumor heterogeneity and identify interventions to restrict it., Competing Interests: The authors declare no conflicts of interest., (© 2024 The Author(s).)

Epigenetic alterations that lead to differential expression of microRNAs (miRNAs/miR) are known to regulate tumour cell states, epithelial-mesenchymal transition (EMT) and the progression to metastasis in breast cancer. This study explores the key contribution of miRNA-18a in mediating a hybrid E/M cell state that is pivotal to the malignant transformation and tumour progression in the aggressive ER-negative subtype of breast cancer. The expression status and associated effects of miR-18a were evaluated in patient-derived breast tumour samples in combination with gene expression data from public datasets, and further validated in in vitro and in vivo breast cancer model systems. The clinical relevance of the study findings was corroborated against human breast tumour specimens (n = 446 patients). The down-regulated expression of miR-18a observed in ER-negative tumours was found to drive the enrichment of hybrid epithelial/mesenchymal (E/M) cells with luminal attributes, enhanced traits of migration, stemness, drug-resistance and immunosuppression. Further analysis of the miR-18a targets highlighted possible hypoxia-inducible factor 1-alpha (HIF-1α)-mediated signalling in these tumours. This is a foremost report that validates the dual role of miR-18a in breast cancer that is subtype-specific based on hormone receptor expression. The study also features a novel association of low miR-18a levels and subsequent enrichment of hybrid E/M cells, increased migration and stemness in a subgroup of ER-negative tumours that may be attributed to HIF-1α mediated signalling. The results highlight the possibility of stratifying the ER-negative disease into clinically relevant groups by analysing miRNA signatures.

Although suppressed cAMP levels have been linked to cancer for nearly five decades, the molecular basis remains uncertain. Here, we identify endosomal pH as a novel regulator of cytosolic cAMP homeostasis and a promoter of transformed phenotypic traits in colorectal cancer. Combining experiments and computational analysis, we show that the Na+/H+ exchanger NHE9 contributes to proton leak and causes luminal alkalinization, which induces resting [Ca2+], and in consequence, represses cAMP levels, creating a feedback loop that echoes nutrient deprivation or hypoxia. Higher NHE9 expression in cancer epithelia is associated with a hybrid epithelial-mesenchymal (E/M) state, poor prognosis, tumor budding, and invasive growth in vitro and in vivo. These findings point to NHE9-mediated cAMP suppression as a pseudostarvation-induced invasion state and potential therapeutic vulnerability in colorectal cancer. Our observations lay the groundwork for future research into the complexities of endosome-driven metabolic reprogramming and phenotype switching and the biology of cancer progression., Implications: Endosomal pH regulator NHE9 actively controls cytosolic Ca2+ levels to downregulate the adenylate cyclase-cAMP system, enabling colorectal cancer cells to acquire hybrid E/M characteristics and promoting metastatic progression., (©2024 American Association for Cancer Research.)

Epithelial to Mesenchymal transition (EMT) drives cancer metastasis and is governed by genetic and epigenetic alterations at multiple levels of regulation. It is well established that loss/mutation of p53 confers oncogenic function to cancer cells and promotes metastasis. Though transcription factors like ZEB1, SLUG, SNAIL and TWIST have been implied in EMT signalling, p53 mediated alterations in the epigenetic machinery accompanying EMT are not clearly understood. This work attempts to explore epigenetic signalling during EMT in colorectal cancer (CRC) cells with varying status of p53. Towards this, we have induced EMT using TGFβ on CRC cell lines with wild type, null and mutant p53 and have assayed epigenetic alterations after EMT induction. Transcriptomic profiling of the four CRC cell lines revealed that the loss of p53 confers more mesenchymal phenotype with EMT induction than its mutant counterparts. This was also accompanied by upregulation of epigenetic writer and eraser machinery suggesting an epigenetic signalling cascade triggered by TGFβ signalling in CRC. Significant agonist and antagonistic relationships observed between EMT factor SNAI1 and SNAI2 with epigenetic enzymes KDM6A/6B and the chromatin organiser SATB1 in p53 null CRC cells suggest a crosstalk between epigenetic and EMT factors. The observed epigenetic regulation of EMT factor SNAI1 correlates with poor clinical outcomes in 270 colorectal cancer patients taken from TCGA-COAD. This unique p53 dependent interplay between epigenetic enzymes and EMT factors in CRC cells may be exploited for development of synergistic therapies for CRC patients presenting to the clinic with loss of p53., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

The epithelial-to-mesenchymal transition (EMT) provides crucial insights into the metastatic process and possesses prognostic value within the cancer context. Here, we present COMET, an R package for inferring EMT trajectories and inter-state transition rates from single-cell RNA sequencing data. We describe steps for finding the optimal number of EMT genes for a specific context, estimating EMT-related trajectories, optimal fitting of continuous-time Markov chain to inferred trajectories, and estimating inter-state transition rates., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

In this Perspective, Journal of Cell Science invited researchers working on cell and tissue polarity to share their thoughts on unique, emerging or open questions relating to their field. The goal of this article is to feature 'voices' from scientists around the world and at various career stages, to bring attention to innovative and thought-provoking topics of interest to the cell biology community. These voices discuss intriguing questions that consider polarity across scales, evolution, development and disease. What can yeast and protists tell us about the evolution of cell and tissue polarity in animals? How are cell fate and development influenced by emerging dynamics in cell polarity? What can we learn from atypical and extreme polarity systems? How can we arrive at a more unified biophysical understanding of polarity? Taken together, these pieces demonstrate the broad relevance of the fascinating phenomenon of cell polarization to diverse fundamental biological questions., Competing Interests: Competing interests The authors declare no competing or financial interests., (© 2024. Published by The Company of Biologists Ltd.)

As the SIRT i analogue series (HL1-HL6) show potent antitumor activity in vitro, we synthesized their corresponding zinc(II) complexes (ZnL1-ZnL6) and investigated their potential as anticancer agents. The Zn(II) complexes showed substantially greater cytotoxicity than HL1-HL6 alone in several cancer cell-types. Notably, distinct structure-activity relationships confirmed the significance of tert -butyl (ZnL2) pharmacophore inclusion in their activity. ZnL2 complexes were found to transmetalate with copper ions inside cells, causing the formation of redox-active copper complexes that induced reactive oxygen species (ROS) production, mitochondrial membrane depolarization, ATP decay, and cell death. This is the first study to exhibit Zn(II) complexes that mediate their activity via transmetalation with copper ions to undergo paraptosis cell death pathway. To further confirm if the SIRT1/2 inhibitory property of SIRTi analogues is conserved, a docking simulation study is performed. The binding affinity and specific interactions of the Cu(II) complex obtained after transmetalation with ZnL2 were found to be higher for SIRT2 ( K i = 0.06 μM) compared to SIRT1 ( K i = 0.25 μM). Thus, the concurrent regulation of several biological targets using a single drug has been shown to have synergistic therapeutic effects, which are crucial for the effective treatment of cancer., Competing Interests: The authors declare no competing financial interest., (© 2024 The Authors. Published by American Chemical Society.)

The aggressive nature of glioblastoma (GBM) - one of the deadliest forms of brain tumors - is majorly attributed to underlying phenotypic heterogeneity. Early attempts to classify this heterogeneity at a transcriptomic level in TCGA GBM cohort proposed the existence of four distinct molecular subtypes: Proneural, Neural, Classical, and Mesenchymal. Further, a single-cell RNA sequencing (scRNA-seq) analysis of primary tumors also reported similar four subtypes mimicking neurodevelopmental lineages. However, it remains unclear whether these four subtypes identified via bulk and single-cell transcriptomics are mutually exclusive or not. Here, we perform pairwise correlations among individual genes and gene signatures corresponding to these proposed subtypes and show that the subtypes are not distinctly mutually antagonistic in either TCGA or scRNA-seq data. We observed that the proneural (or neural progenitor-like)-mesenchymal axis is the most prominent antagonistic pair, with the other two subtypes lying on this spectrum. These results are reinforced through a meta-analysis of over 100 single-cell and bulk transcriptomic datasets as well as in terms of functional association with metabolic switching, cell cycle, and immune evasion pathways. Finally, this proneural-mesenchymal antagonistic trend percolates to the association of relevant transcription factors with patient survival. These results suggest rethinking GBM phenotypic characterization for more effective therapeutic targeting efforts., Competing Interests: The authors declare no conflict of interest., (© 2024 The Authors.)

Colorectal cancer (CRC) is highly heterogeneous with variable survival outcomes and therapeutic vulnerabilities. A commonly used classification system in CRC is the Consensus Molecular Subtypes (CMS) based on gene expression patterns. However, how these CMS categories connect to axes of phenotypic plasticity and heterogeneity remains unclear. Here, in our analysis of CMS-specific TCGA data and 101 bulk transcriptomic datasets, we found the epithelial phenotype score to be consistently positively correlated with scores of glycolysis, OXPHOS and FAO pathways, while PD-L1 activity scores positively correlated with mesenchymal phenotype scoring, revealing possible interconnections among plasticity axes. Single-cell RNA-sequencing analysis of patient samples revealed that that CMS2 and CMS3 subtype samples were relatively more epithelial as compared to CMS1 and CMS4. CMS1 revealed two subpopulations: one close to CMS4 (more mesenchymal) and the other closer to CMS2 or CMS3 (more epithelial), indicating a partial EMT-like behavior. Consistent observations were made in single-cell analysis of metabolic axes and PD-L1 activity scores. Together, our results quantify the patterns of two functional interconnected axes of phenotypic heterogeneity - EMT and metabolic reprogramming - in a CMS-specific manner in CRC., Competing Interests: Declaration of Competing Interest The authors declare no conflict of interest., (Copyright © 2023. Published by Elsevier Inc.)

Ectodermal appendages, such as the mammary gland (MG), are thought to have evolved from hair-associated apocrine glands to serve the function of milk secretion. Through the directed differentiation of mouse embryonic stem cells (mESCs), here, we report the generation of multilineage ESC-derived mammary organoids (MEMOs). We adapted the skin organoid model, inducing the dermal mesenchyme to transform into mammary-specific mesenchyme via the sequential activation of Bone Morphogenetic Protein 4 (BMP4) and Parathyroid Hormone-related Protein (PTHrP) and inhibition of hedgehog (HH) signaling. Using single-cell RNA sequencing, we identified gene expression profiles that demonstrate the presence of mammary-specific epithelial cells, fibroblasts, and adipocytes. MEMOs undergo ductal morphogenesis in Matrigel and can reconstitute the MG in vivo. Further, we demonstrate that the loss of function in placode regulators LEF1 and TBX3 in mESCs results in impaired skin and MEMO generation. In summary, our MEMO model is a robust tool for studying the development of ectodermal appendages, and it provides a foundation for regenerative medicine and disease modeling., Competing Interests: Declaration of interests The authors declare no competing interests., (Published by Elsevier Inc.)

The epithelial-mesenchymal transition (EMT) program is crucial for transforming carcinoma cells into a partially mesenchymal state, enhancing their chemoresistance, migration, and metastasis. This shift in cell state is tightly regulated by cellular mechanisms that are not yet fully characterized. One intriguing EMT aspect is the rewiring of the proteoglycan landscape, particularly the induction of heparan sulfate proteoglycan (HSPG) biosynthesis. This proteoglycan functions as a co-receptor that accelerates cancer-associated signaling pathways through its negatively-charged residues. However, the precise mechanisms through which EMT governs HSPG biosynthesis and its role in cancer cell plasticity remain elusive. Here, we identified exostosin glycosyltransferase 1 (EXT1), a central enzyme in HSPG biosynthesis, to be selectively upregulated in aggressive tumor subtypes and cancer cell lines, and to function as a key player in breast cancer aggressiveness. Notably, ectopic expression of EXT1 in epithelial cells is sufficient to induce HSPG levels and the expression of known mesenchymal markers, subsequently enhancing EMT features, including cell migration, invasion, and tumor formation. Additionally, EXT1 loss in MDA-MB-231 cells inhibits their aggressiveness-associated traits such as migration, chemoresistance, tumor formation, and metastasis. Our findings reveal that EXT1, through its role in HSPG biosynthesis, governs signal transducer and activator of transcription 3 (STAT3) signaling, a known regulator of cancer cell aggressiveness. Collectively, we present the EXT1/HSPG/STAT3 axis as a central regulator of cancer cell plasticity that directly links proteoglycan synthesis to oncogenic signaling pathways., Competing Interests: Competing interests statement:The authors declare no competing interest.

High-grade serous ovarian carcinoma (HGSC) is associated with late-stage disease presentation and poor prognosis, with a limited understanding of early transformation events. Our study analyzes HGSC tumor progression and organ-specific metastatic dissemination to identify hypoxia-associated molecular, cellular, and histological alterations. Clinical characteristics of the HGSC were replicated in orthotopic xenografts, which involve metastatic dissemination and the prevalence of group B tumors (volume: >0.0625 ≤ 0.5 cm3). Enhanced hyaluronic acid (HA) deposition, expanded tumor vasculature, and increased necrosis contributed to the remodeling of tumor tissue architecture. The proliferative potential of tumor cells and the ability to form glands were also altered during tumor growth. Flow cytometry and label chase-based molecular profiling across the tumor regenerative hierarchy identified the hypoxia-vasculogenic niche and the hybrid epithelial-mesenchymal tumor-cell state as determinants of self-renewal capabilities of progenitors and cancer stem cells. A regulatory network and mathematical model based on tumor histology and molecular signatures predicted hypoxia-inducible factor 1-alpha (HIF1A) as a central node connecting HA synthesis, epithelial-mesenchymal transition, metabolic, vasculogenic, inflammatory, and necrotic pathways in HGSC tumors. Thus, our findings provide a temporal resolution of hypoxia-associated events that sculpt HGSC tumor growth; an in-depth understanding of it may aid in the early detection and treatment of HGSC., (© 2022 S. Karger AG, Basel.)

Background: Epithelial-to-mesenchymal transition (EMT) is a developmental program that consists of the loss of epithelial features concomitant with the acquisition of mesenchymal features. Activation of EMT in cancer facilitates the acquisition of aggressive traits and cancer invasion. EMT plasticity (EMP), the dynamic transition between multiple hybrid states in which cancer cells display both epithelial and mesenchymal markers, confers survival advantages for cancer cells in constantly changing environments during metastasis., Methods: RNAseq analysis was performed to assess genome-wide transcriptional changes in cancer cells depleted for histone regulators FLASH, NPAT, and SLBP. Quantitative PCR and Western blot were used for the detection of mRNA and protein levels. Computational analysis was performed on distinct sets of genes to determine the epithelial and mesenchymal score in cancer cells and to correlate FLASH expression with EMT markers in the CCLE collection., Results: We demonstrate that loss of FLASH in cancer cells gives rise to a hybrid E/M phenotype with high epithelial scores even in the presence of TGFβ, as determined by computational methods using expression of predetermined sets of epithelial and mesenchymal genes. Multiple genes involved in cell-cell junction formation are similarly specifically upregulated in FLASH-depleted cells, suggesting that FLASH acts as a repressor of the epithelial phenotype. Further, FLASH expression in cancer lines is inversely correlated with the epithelial score. Nonetheless, subsets of mesenchymal markers were distinctly up-regulated in FLASH, NPAT, or SLBP-depleted cells., Conclusions: The ZEB1 low /SNAIL high /E-cadherin high phenotype described in FLASH-depleted cancer cells is driving a hybrid E/M phenotype in which epithelial and mesenchymal markers coexist., Competing Interests: Declaration of Competing Interest The authors declare no conflict of interest., (Copyright © 2023. Published by Elsevier Inc.)

Multistability is central to biological systems. It plays a crucial role in adaptation, evolvability, and differentiation. The presence of positive feedback loops can enable multistability. The simplest of such feedback loops are (a) a mutual inhibition (MI) loop, (b) a mutual activation (MA) loop, and (c) self-activation. While it is established that all three motifs can give rise to bistability, the characteristic differences in the bistability exhibited by each of these motifs is relatively less understood. Here, we use dynamical simulations across a large ensemble of parameter sets and initial conditions to study the bistability characteristics of these motifs. Furthermore, we investigate the utility of these motifs for achieving coordinated expression through cyclic and parallel coupling amongst them. Our analysis revealed that MI-based architectures offer discrete and robust control over gene expression, multistability, and coordinated expression among multiple genes, as compared to MA-based architectures. We then devised a combination of MI and MA architectures to improve coordination and multistability. Such designs help enhance our understanding of the control structures involved in robust cell-fate decisions and provide a way to achieve controlled decision-making in synthetic systems., (© 2023 S. Karger AG, Basel.)

Intra-tumor heterogeneity contributes to treatment failure and poor survival in urothelial bladder carcinoma (UBC). Analyzing transcriptome from a UBC cohort, we report that intra-tumor transcriptomic heterogeneity indicates co-existence of tumor cells in epithelial and mesenchymal-like transcriptional states and bi-directional transition between them occurs within and between tumor subclones. We model spontaneous and reversible transition between these partially heritable states in cell lines and characterize their population dynamics. SMAD3, KLF4 and PPARG emerge as key regulatory markers of the transcriptional dynamics. Nutrient limitation, as in the core of large tumors, and radiation treatment perturb the dynamics, initially selecting for a transiently resistant phenotype and then reconstituting heterogeneity and growth potential, driving adaptive evolution. Dominance of transcriptional states with low PPARG expression indicates an aggressive phenotype in UBC patients. We propose that phenotypic plasticity and dynamic, non-genetic intra-tumor heterogeneity modulate both the trajectory of disease progression and adaptive treatment response in UBC., (© 2023. The Author(s).)

MicroRNAs (miRNAs) are small noncoding RNAs that regulate gene expression post-transcriptionally in eukaryotes by binding with target mRNAs and preventing translation. miRNA-mediated feedback motifs are ubiquitous in various genetic networks that control cellular decision making. A key question is how such a feedback mechanism may affect gene expression noise. To answer this, we have developed a mathematical model to study the effects of a miRNA-dependent negative-feedback loop on mean expression and noise in target mRNAs. Combining analytics and simulations, we show the existence of an expression threshold demarcating repressed and expressed regimes in agreement with earlier studies. The steady-state mRNA distributions are bimodal near the threshold, where copy numbers of mRNAs and miRNAs exhibit enhanced anticorrelated fluctuations. Moreover, variation of negative-feedback strength shifts the threshold locations and modulates the noise profiles. Notably, the miRNA-mRNA binding affinity and feedback strength collectively shape the bimodality. We also compare our model with a direct auto-repression motif, where a gene produces its own repressor. Auto-repression fails to produce bimodal mRNA distributions as found in miRNA-based indirect repression, suggesting the crucial role of miRNAs in creating phenotypic diversity. Together, we demonstrate how miRNA-dependent negative feedback modifies the expression threshold and leads to a broader parameter regime of bimodality compared to the no-feedback case., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Phenotypic plasticity was recently incorporated as a hallmark of cancer. This plasticity can manifest along many interconnected axes, such as stemness and differentiation, drug-sensitive and drug-resistant states, and between epithelial and mesenchymal cell-states. Despite growing acceptance for phenotypic plasticity as a hallmark of cancer, the dynamics of this process remains poorly understood. In particular, the knowledge necessary for a predictive understanding of how individual cancer cells and populations of cells dynamically switch their phenotypes in response to the intensity and/or duration of their current and past environmental stimuli remains far from complete. Here, we present recent investigations of phenotypic plasticity from a systems-level perspective using two exemplars: epithelial-mesenchymal plasticity in carcinomas and phenotypic switching in melanoma. We highlight how an integrated computational-experimental approach has helped unravel insights into specific dynamical hallmarks of phenotypic plasticity in different cancers to address the following questions: a) how many distinct cell-states or phenotypes exist?; b) how reversible are transitions among these cell-states, and what factors control the extent of reversibility?; and c) how might cell-cell communication be able to alter rates of cell-state switching and enable diverse patterns of phenotypic heterogeneity? Understanding these dynamic features of phenotypic plasticity may be a key component in shifting the paradigm of cancer treatment from reactionary to a more predictive, proactive approach., Competing Interests: Declaration of Competing Interest The authors declare no conflict of interests., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Epithelial-mesenchymal transition (EMT) is an important axis of phenotypic plasticity-a hallmark of cancer metastasis. Raf kinase-B inhibitor protein (RKIP) and BTB and CNC homology 1 (BACH1) are reported to influence EMT. In breast cancer, they act antagonistically, but the exact nature of their roles in mediating EMT and associated other axes of plasticity remains unclear. Here, analysing transcriptomic data, we reveal their antagonistic trends in a pan-cancer manner in terms of association with EMT, metabolic reprogramming and immune evasion via PD-L1. Next, we developed and simulated a mechanism-based gene regulatory network that captures how RKIP and BACH1 engage in feedback loops with drivers of EMT and stemness. We found that RKIP and BACH1 belong to two antagonistic 'teams' of players-while BACH1 belonged to the one driving pro-EMT, stem-like and therapy-resistant cell states, RKIP belonged to the one enabling pro-epithelial, less stem-like and therapy-sensitive phenotypes. Finally, we observed that low RKIP levels and upregulated BACH1 levels associated with worse clinical outcomes in many cancer types. Together, our systems-level analysis indicates that the emergent dynamics of underlying regulatory network enable the antagonistic patterns of RKIP and BACH1 with various axes of cancer cell plasticity, and with patient survival data.

Background: Androgen receptor (AR) is considered a marker of better prognosis in hormone receptor positive breast cancers (BC), however, its role in triple negative breast cancer (TNBC) is controversial. This may be attributed to intrinsic molecular differences or scoring methods for AR positivity. We derived AR regulated gene score and examined its utility in BC subtypes., Methods: AR regulated genes were derived by applying a bioinformatic pipeline on publicly available microarray data sets of AR+ BC cell lines and gene score was calculated as average expression of six AR regulated genes. Tumors were divided into AR high and low based on gene score and associations with clinical parameters, circulating androgens, survival and epithelial to mesenchymal transition (EMT) markers were examined, further evaluated in invitro models and public datasets., Results: 53% (133/249) tumors were classified as AR gene score high and were associated with significantly better clinical parameters, disease-free survival (86.13 vs 72.69 months, log rank p = 0.032) when compared to AR low tumors. 36% of TNBC (N = 66) were AR gene score high with higher expression of EMT markers (p = 0.024) and had high intratumoral levels of 5α-reductase, enzyme involved in intracrine androgen metabolism. In MDA-MB-453 treated with dihydrotestosterone, SLUG expression increased, E-cadherin decreased with increase in migration and these changes were reversed with bicalutamide. Similar results were obtained in public datasets., Conclusion: Deciphering the role of AR in BC is difficult based on AR protein levels alone. Our results support the context dependent function of AR in driving better prognosis in ER positive tumors and EMT features in TNBC tumors., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023. Published by Elsevier Inc.)

Gene expression control based on clustered regularly interspaced short palindromic repeats (CRISPR) has emerged as a powerful approach for constructing synthetic gene circuits. While the use of CRISPR interference (CRISPRi) is already well-established in prokaryotic circuits, CRISPR activation (CRISPRa) is less mature, and a combination of the two in the same circuits is only just emerging. Here, we report that combining CRISPRi with SoxS-based CRISPRa in Escherichia coli can lead to context-dependent effects due to different affinities in the formation of CRISPRa and CRISPRi complexes, resulting in loss of predictable behavior. We show that this effect can be avoided by using the same scaffold guide RNA structure for both complexes.

In the epithelium, cell density and cell proliferation are closely connected to each other through contact inhibition of proliferation (CIP). Depending on cell density, CIP proceeds through three distinct stages: the free-growing stage at low density, the pre-epithelial transition stage at medium density, and the post-epithelial transition stage at high density. Previous studies have elucidated how cell morphology, motion, and mechanics vary in these stages. However, it remains unknown whether cellular metabolism also has a density-dependent behavior. By measuring the mitochondrial membrane potential at different cell densities, here we reveal a heterogeneous landscape of metabolism in the epithelium, which appears qualitatively distinct in three stages of CIP and did not follow the trend of other CIP-associated parameters, which increases or decreases monotonically with increasing cell density. Importantly, epithelial cells established a collective metabolic heterogeneity exclusively in the pre-epithelial transition stage, where the multicellular clusters of high- and low-potential cells emerged. However, in the post-epithelial transition stage, the metabolic potential field became relatively homogeneous. Next, to study the underlying dynamics, we constructed a system biology model, which predicted the role of cell proliferation in metabolic potential toward establishing collective heterogeneity. Further experiments indeed revealed that the metabolic pattern spatially correlated with the proliferation capacity of cells, as measured by the nuclear localization of a pro-proliferation protein, YAP. Finally, experiments perturbing the actomyosin contractility revealed that, while metabolic heterogeneity was maintained in the absence of actomyosin contractility, its ab initio emergence depended on the latter. Taken together, our results revealed a density-dependent collective heterogeneity in the metabolic field of a pre-epithelial transition-stage epithelial monolayer, which may have significant implications for epithelial form and function., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

The triple negative breast cancer (TNBC) subtype is one of the most aggressive forms of breast cancer that has poor clinical outcome and is an unmet clinical challenge. Accumulating evidence suggests that intratumoral heterogeneity or the presence of phenotypically distinct cell populations within a tumor play a crucial role in chemoresistance, tumor progression and metastasis. An increased understanding of the molecular regulators of intratumoral heterogeneity is crucial to the development of effective therapeutic strategies in TNBC. To this end, we used an unbiased approach to identify a molecular mediator of intratumoral heterogeneity in breast cancer by isolating two tumor cell populations (T1 and T2) from the 4T1 TNBC model. Phenotypic characterization revealed that the cells are different in terms of their morphology, proliferation and self-renewal ability in vitro as well as primary tumor formation and metastatic potential in vivo . Bioinformatic analysis followed by Kaplan Meier survival analysis in TNBC patients identified Metastasis associated colon cancer 1 (Macc1) as one of the top candidate genes mediating the aggressive phenotype in the T1 tumor cells. The role of Macc1 in regulating the proliferative phenotype was validated and taken forward in a therapeutic context with Lovastatin, a small molecule transcriptional inhibitor of Macc1 to target the T1 cell population. This study increases our understanding of the molecular underpinnings of intratumoral heterogeneity in breast cancer that is critical to improve the treatment of women currently living with the highly aggressive TNBC subtype., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Thankamony, Ramkomuth, Ramesh, Murali, Chakraborty, Karthikeyan, Varghese, Jaikumar, Jolly, Swarbrick and Nair.)

p53 (also known as TP53) mutation and amyloid formation are long associated with cancer pathogenesis; however, the direct demonstration of the link between p53 amyloid load and cancer progression is lacking. Using multi-disciplinary techniques and 59 tissues (53 oral and stomach cancer tumor tissue samples from Indian individuals with cancer and six non-cancer oral and stomach tissue samples), we showed that p53 amyloid load and cancer grades are highly correlated. Furthermore, next-generation sequencing (NGS) data suggest that not only mutant p53 (e.g. single-nucleotide variants, deletions, and insertions) but wild-type p53 also formed amyloids either in the nucleus (50%) and/or in the cytoplasm in most cancer tissues. Interestingly, in all these cancer tissues, p53 displays a loss of DNA-binding and transcriptional activities, suggesting that the level of amyloid load correlates with the degree of loss and an increase in cancer grades. The p53 amyloids also sequester higher amounts of the related p63 and p73 (also known as TP63 and TP73, respectively) protein in higher-grade tumor tissues. The data suggest p53 misfolding and/or aggregation, and subsequent amyloid formation, lead to loss of the tumor-suppressive function and the gain of oncogenic function, aggravation of which might determine the cancer grade., Competing Interests: Competing interests The authors declare no competing or financial interests., (© 2023. Published by The Company of Biologists Ltd.)

Background: Phenotypic heterogeneity of melanoma cells contributes to drug tolerance, increased metastasis, and immune evasion in patients with progressive disease. Diverse mechanisms have been individually reported to shape extensive intra-tumor and inter-tumor phenotypic heterogeneity, such as IFNγ signaling and proliferative to invasive transition, but how their crosstalk impacts tumor progression remains largely elusive., Methods: Here, we integrate dynamical systems modeling with transcriptomic data analysis at bulk and single-cell levels to investigate underlying mechanisms behind phenotypic heterogeneity in melanoma and its impact on adaptation to targeted therapy and immune checkpoint inhibitors. We construct a minimal core regulatory network involving transcription factors implicated in this process and identify the multiple 'attractors' in the phenotypic landscape enabled by this network. Our model predictions about synergistic control of PD-L1 by IFNγ signaling and proliferative to invasive transition were validated experimentally in three melanoma cell lines-MALME3, SK-MEL-5 and A375., Results: We demonstrate that the emergent dynamics of our regulatory network comprising MITF, SOX10, SOX9, JUN and ZEB1 can recapitulate experimental observations about the co-existence of diverse phenotypes (proliferative, neural crest-like, invasive) and reversible cell-state transitions among them, including in response to targeted therapy and immune checkpoint inhibitors. These phenotypes have varied levels of PD-L1, driving heterogeneity in immunosuppression. This heterogeneity in PD-L1 can be aggravated by combinatorial dynamics of these regulators with IFNγ signaling. Our model predictions about changes in proliferative to invasive transition and PD-L1 levels as melanoma cells evade targeted therapy and immune checkpoint inhibitors were validated in multiple RNA-seq data sets from in vitro and in vivo experiments., Conclusion: Our calibrated dynamical model offers a platform to test combinatorial therapies and provide rational avenues for the treatment of metastatic melanoma. This improved understanding of crosstalk among PD-L1 expression, proliferative to invasive transition and IFNγ signaling can be leveraged to improve the clinical management of therapy-resistant and metastatic melanoma., Competing Interests: Competing interests: None declared., (© Author(s) (or their employer(s)) 2023. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.)

The last decade has witnessed a surge of theoretical and computational models to describe the dynamics of complex gene regulatory networks, and how these interactions can give rise to multistable and heterogeneous cell populations. As the use of theoretical modeling to describe genetic and biochemical circuits becomes more widespread, theoreticians with mathematical and physical backgrounds routinely apply concepts from statistical physics, non-linear dynamics, and network theory to biological systems. This review aims at providing a clear overview of the most important methodologies applied in the field while highlighting current and future challenges. It also includes hands-on tutorials to solve and simulate some of the archetypical biological system models used in the field. Furthermore, we provide concrete examples from the existing literature for theoreticians that wish to explore this fast-developing field. Whenever possible, we highlight the similarities and differences between biochemical and regulatory networks and 'classical' systems typically studied in non-equilibrium statistical and quantum mechanics., (© 2023 IOP Publishing Ltd.)

Esophageal squamous cell carcinoma (ESCC) is the most prevalent malignant gastrointestinal tumor. Ion channels contribute to tumor growth and progression through interactions with their neighboring molecules including lipids. The dysregulation of membrane ion channels and lipid metabolism may contribute to the epithelial-mesenchymal transition (EMT), leading to metastatic progression. Herein, transcriptome profiles of patients with ESCC were analyzed by performing differential gene expression and weighted gene co-expression network analysis to identify the altered ion channels, lipid metabolism- and EMT-related genes in ESCC. A total of 1,081 differentially expressed genes, including 113 ion channels, 487 lipid metabolism-related, and 537 EMT-related genes, were identified in patients with ESCC. Thereafter, EMT scores were correlated with altered co-expressed genes. The altered co-expressed genes indicated a correlation with EMT signatures. Interactions among 22 ion channels with 3 hub lipid metabolism- and 13 hub EMT-related proteins were determined using protein-protein interaction networks. A pathway map was generated to depict deregulated signaling pathways including insulin resistance and the estrogen receptor-Ca 2+ signaling pathway in ESCC. The relationship between potential ion channels and 5-year survival rates in ESCC was determined using Kaplan-Meier plots and Cox proportional hazard regression analysis. Inositol 1,4,5-trisphosphate receptor type 3 ( ITPR3 ) was found to be associated with poor prognosis of patients with ESCC. Additionally, drugs interacting with potential ion channels, including GJA1 and ITPR3 , were identified. Understanding alterations in ion channels with lipid metabolism and EMT in ESCC pathophysiology would most likely provide potential targets for the better treatment of patients with ESCC., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Parthasarathi, Mandal, George, Gaikwad, Sasidharan, Gundimeda, Jolly, Pandey and Sharma.)

Mathematical modeling of the emergent dynamics of gene regulatory networks (GRN) faces a double challenge of (a) dependence of model dynamics on parameters, and (b) lack of reliable experimentally determined parameters. In this paper we compare two complementary approaches for describing GRN dynamics across unknown parameters: (1) parameter sampling and resulting ensemble statistics used by RACIPE (RAndom CIrcuit PErturbation), and (2) use of rigorous analysis of combinatorial approximation of the ODE models by DSGRN (Dynamic Signatures Generated by Regulatory Networks). We find a very good agreement between RACIPE simulation and DSGRN predictions for four different 2- and 3-node networks typically observed in cellular decision making. This observation is remarkable since the DSGRN approach assumes that the Hill coefficients of the models are very high while RACIPE assumes the values in the range 1-6. Thus DSGRN parameter domains, explicitly defined by inequalities between systems parameters, are highly predictive of ODE model dynamics within a biologically reasonable range of parameters., (© 2023. The Author(s).)

Consensus decision-making in social groups strongly depends on communication links that determine to whom individuals send, and from whom they receive, information. Here, we ask how consensus decisions are affected by strategic updating of links and how this effect varies with the direction of communication. We quantified the coevolution of link and opinion dynamics in a large population with binary opinions using mean-field numerical simulations of two voter-like models of opinion dynamics: an incoming model (IM) (where individuals choose who to receive opinions from) and an outgoing model (OM) (where individuals choose who to send opinions to). We show that individuals can bias group-level outcomes in their favour by breaking disagreeing links while receiving opinions (IM) and retaining disagreeing links while sending opinions (OM). Importantly, these biases can help the population avoid stalemates and achieve consensus. However, the role of disagreement avoidance is diluted in the presence of strong preferences-highly stubborn individuals can shape decisions to favour their preferences, giving rise to non-consensus outcomes. We conclude that collectively changing communication structures can bias consensus decisions, as a function of the strength of preferences and the direction of communication., Competing Interests: We declare we have no competing interests., (© 2023 The Authors.)