TGFβ1-Induced EMT in the MCF10A Mammary Epithelial Cell Line Model Is Executed Independently of SNAIL1 and ZEB1 but Relies on JUNB-Coordinated Transcriptional Regulation

Correction: Antón-García et al. TGFβ1-Induced EMT in the MCF10A Mammary Epithelial Cell Line Model Is Executed Independently of SNAIL1 and ZEB1 but Relies on JUNB-Coordinated Transcriptional Regulation. Cancers 2023, 15, 588

In the original publication [...].

Single-cell morphodynamical trajectories enable prediction of gene expression accompanying cell state change

Extracellular signals induce changes to molecular programs that modulate multiple cellular phenotypes, including proliferation, motility, and differentiation status. The connection between dynamically adapting phenotypic states and the molecular programs that define them is not well understood. Here we develop data-driven models of single-cell phenotypic responses to extracellular stimuli by linking gene transcription levels to "morphodynamics" - changes in cell morphology and motility observable in time-lapse image data. We adopt a dynamics-first view of cell state by grouping single-cell trajectories into states with shared morphodynamic responses. The single-cell trajectories enable development of a first-of-its-kind computational approach to map live-cell dynamics to snapshot gene transcript levels, which we term MMIST, Molecular and Morphodynamics-Integrated Single-cell Trajectories. The key conceptual advance of MMIST is that cell behavior can be quantified based on dynamically defined states and that extracellular signals alter the overall distribution of cell states by altering rates of switching between states. We find a cell state landscape that is bound by epithelial and mesenchymal endpoints, with distinct sequences of epithelial to mesenchymal transition (EMT) and mesenchymal to epithelial transition (MET) intermediates. The analysis yields predictions for gene expression changes consistent with curated EMT gene sets and provides a prediction of thousands of RNA transcripts through extracellular signal-induced EMT and MET with near-continuous time resolution. The MMIST framework leverages true single-cell dynamical behavior to generate molecular-level omics inferences and is broadly applicable to other biological domains, time-lapse imaging approaches and molecular snapshot data.

Transforming growth factor β1 upregulates 6-phosphofructo-2-kinase/fructose 2,6-bisphosphatase-4 expression in A549 and MCF-10A cells

Transforming growth factor β1 (TGFβ1) induces a cellular process known as epithelial-mesenchymal transition (EMT) associated with metabolic reprogramming, including enhanced glycolysis. Given the involvement of 6-phosphofructo-2-kinase/fructose 2,6-bisphosphatase (PFKFB) enzymes in glycolysis, we aimed to investigate whether TGFβ1 regulates expressions of PFKFB genes and if PFKFBs are required for TGFβ1-driven phenotypes. A549 and MCF-10A cell lines were used as TGFβ1-driven EMT models. Messenger RNA expressions of PFKFB and EMT genes were determined by real-time quantitative polymerase chain reaction. A small interfering RNA approach was used to deplete PFKFB4 expression. A Matrigel invasion assay was conducted to assess the effect of PFKFB4 silencing on the TGFβ1-enhanced invasion of A549 cells. F2,6BP levels were analyzed using an enzyme-coupled assay. Glucose and lactate concentrations were determined using colorimetric assays. TGFβ1 robustly induced expression of the fourth isoform of PFKFBs, PFKFB4, in both cell lines. PFKFB4 depletion partially inhibits mesenchymal transdifferentiation caused by TGFβ1 in A549 cells, as assessed by microscopy. Inductions of Snail in MCF-10A cells and Fibronectin in A549 cells and repressions of E-cadherin in both cell lines by TGFβ1 are attenuated by PFKFB4 silencing. PFKFB4 silencing reduces F2,6BP and glycolytic activity, although TGFβ1 alone does not affect these parameters. Finally, PFKFB4 depletion suppresses the TGFβ1-driven invasion of A549 cells through Matrigel. Presented data suggest that TGFβ1 induces the expression of PFKFB4 in A549 and MCF-10 cells, and PFKFB4 may be required for TGFβ1-driven phenotypes such as EMT and invasion in these models.

Population dynamics of EMT elucidates the timing and distribution of phenotypic intra-tumoral heterogeneity

The Epithelial-to-Mesenchymal Transition (EMT) is a hallmark of cancer metastasis and morbidity. EMT is a non-binary process, and cells can be stably arrested en route to EMT in an intermediate hybrid state associated with enhanced tumor aggressiveness and worse patient outcomes. Understanding EMT progression in detail will provide fundamental insights into the mechanisms underlying metastasis. Despite increasingly available single-cell RNA sequencing (scRNA-seq) data that enable in-depth analyses of EMT at the single-cell resolution, current inferential approaches are limited to bulk microarray data. There is thus a great need for computational frameworks to systematically infer and predict the timing and distribution of EMT-related states at single-cell resolution. Here, we develop a computational framework for reliable inference and prediction of EMT-related trajectories from scRNA-seq data. Our model can be utilized across a variety of applications to predict the timing and distribution of EMT from single-cell sequencing data.

Regulation of Cdh2 by the AP-1 family transcription factor Junb in TM4 Sertoli cells

Cadherins are transmembrane proteins that mediate cell-to-cell adhesion and various cellular processes. In Sertoli cells of the testis, Cdh2 contributes to the development of the testis and the formation of the blood-testis barrier, being essential for germ cells' protection. Analyses of chromatin accessibility and epigenetic marks in adult mouse testis have shown that the region from -800 to +900 bp respective to Cdh2 transcription start site (TSS) is likely the active regulatory region of this gene. In addition, the JASPAR 2022 matrix has predicted an AP-1 binding element at about -600 bp. Transcription factors of the activator protein 1 (AP-1) family have been implicated in the regulation of the expression of genes encoding cell-to-cell interaction proteins such as Gja1, Nectin2 and Cdh3. To test the potential regulation of Cdh2 by members of the AP-1 family, siRNAs were transfected into TM4 Sertoli cells. The knockdown of Junb led to a decrease in Cdh2 expression. ChIP-qPCR and luciferase reporter assays with site-directed mutagenesis confirmed the recruitment of Junb to several AP-1 regulatory elements in the proximal region of the Cdh2 promoter in TM4 cells. Further investigation with luciferase reporter assays showed that other AP-1 members can also activate the Cdh2 promoter albeit to a lesser extent than Junb. Taken together, these data suggest that in TM4 Sertoli cells, Junb is responsible for the regulation of Cdh2 expression which requires its recruitment to the proximal region of the Cdh2 promoter.