Combinatorial Signal Perception in the BMP Pathway

High order expression dependencies finely resolve cryptic states and subtypes in single cell data

Single cells are typically typed by clustering into discrete locations in reduced dimensional transcriptome space. Here we introduce Stator, a data-driven method that identifies cell (sub)types and states without relying on cells' local proximity in transcriptome space. Stator labels the same single cell multiply, not just by type and subtype, but also by state such as activation, maturity or cell cycle sub-phase, through deriving higher-order gene expression dependencies from a sparse gene-by-cell expression matrix. Stator's finer resolution is clear from analyses of mouse embryonic brain, and human healthy or diseased liver. Rather than only coarse-scale labels of cell type, Stator further resolves cell types into subtypes, and these subtypes into stages of maturity and/or cell cycle phases, and yet further into portions of these phases. Among cryptically homogeneous embryonic cells, for example, Stator finds 34 distinct radial glia states whose gene expression forecasts their future GABAergic or glutamatergic neuronal fate. Further, Stator's fine resolution of liver cancer states reveals expression programmes that predict patient survival. We provide Stator as a Nextflow pipeline and Shiny App.

Bone morphogenetic protein-3 is a negative regulator of transforming growth factor beta and fibrosis

Fibrosis results in one-third of all deaths globally and is a major healthcare challenge. Fibrosis is scarring caused by the excess deposition of extracellular matrix proteins by fibroblasts. Inhibition of pathways downstream of transforming growth factor β (TGF-β) a pluripotent growth factor, has potent antifibrotic effects in different organs. Here we show that loss of bone morphogenetic protein (BMP-3) is a feature of kidney fibrosis, independent of the initiating injury, suggesting loss of this cytokine is a core fibrotic mechanism. TGF-β decreased BMP3 expression in human fibroblasts is possibly a feed-forward loop that contributes to increased and sustained TGF-β activity. Recombinant human BMP-3 reduced TGF-β induced fibroblast contraction, migration and invasion, pathways that lead to scarring and tissue stiffening. BMP-3 reduced TGF-β stimulated collagen cross-linking, and Ox-LDL receptor 1, a regulator of collagen deposition. BMP-3 inhibited TGF-β stimulated lysyl oxidase activity. Lysyl oxidase mediated collagen cross-linking is a critical process in TGF-β induced fibrosis. We propose that BMP-3 alters fibroblast responses to TGF-β, shifting the balance from fibrosis to repair. Recombinant human BMP-3 shows promise for development as a novel therapeutic for fibrosis.

B7H6 is the predominant activating ligand driving natural killer cell-mediated killing in patients with liquid tumours: evidence from clinical, in silico, in vitro, and in vivo studies

BACKGROUND

Natural killer (NK) cells are a subset of innate lymphoid cells that are inherently capable of recognizing and killing infected or tumour cells. This has positioned NK cells as a promising live drug for tumour immunotherapy, but limited success suggests incomplete knowledge of their killing mechanism. NK cell-mediated killing involves a complex decision-making process based on integrating activating and inhibitory signals from various ligand-receptor repertoires. However, the relative importance of the different activating ligand-receptor interactions in triggering NK killing remains unclear.

METHODS

We employed a systematic approach combining clinical, in silico, in vitro, and in vivo data analysis to quantify the impact of various activating ligands. Clinical data analysis was conducted using massive pan-cancer data (n = 10,595), where patients with high NK cell levels were stratified using CIBERSORT. Subsequently, multivariate Cox regression and Kaplan-Meier (KM) survival analysis were performed based on activating ligand expression. To examine the impact of ligand expression on NK killing at the cellular level, we assessed surface expression of five major activating ligands (B7H6, MICA/B, ULBP1, ULBP2/5/6, and ULBP3) of human tumour cell lines of diverse origins (n = 33) via flow cytometry (FACs) and their NK cell-mediated cytotoxicity on by calcein-AM assay using human primary NK cells and NK-92 cell lines. Based on this data, we quantified the contribution of each activating ligand to the NK killing activity using mathematical models and Bayesian statistics. To further validate the results, we performed calcein-AM assays upon ligand knockdown and overexpression, conjugation assays, and co-culture assays in activating ligand-downregulated/overexpressed in liquid tumour (LT) cell lines. Moreover, we established LT-xenograft mouse models to assess the efficacy of NK cell targeting toward tumours with dominant ligands.

FINDINGS

Through the clinical analysis, we discovered that among nearly all 18 activating ligands, only patients with LT who were NK cell-rich and specifically had higher B7H6 level exhibited a favorable survival outcome (p = 0.0069). This unexpected dominant role of B7H6 was further confirmed by the analysis of datasets encompassing multiple ligands and a variety of tumours, which showed that B7H6 exhibited the highest contribution to NK killing among five representative ligands. Furthermore, LT cell lines (acute myeloid leukemia (AML), B cell lymphoma, and T-acute lymphocytic leukemia (ALL)) with lowered B7H6 demonstrated decreased susceptibility to NK cell-mediated cytotoxicity compared to those with higher levels. Even within the same cell line, NK cells selectively targeted cells with higher B7H6 levels. Finally, LT-xenograft mouse models (n = 24) confirmed that higher B7H6 results in less tumour burden and longer survival in NK cell-treated LT mice (p = 0.0022).

INTERPRETATION

While NK cells have gained attention for their potent anti-tumour effects without causing graft-versus-host disease (GvHD), thus making them a promising off-the-shelf therapy, our limited understanding of NK killing mechanisms has hindered their clinical application. This study illuminates the crucial role of the activating ligand B7H6 in driving NK cell killing, particularly in the context of LT. Therefore, the expression level of B7H6 could serve as a prognostic marker for patients with LT. Moreover, for the development of NK cell-based immunotherapy, focusing on increasing the level of B7H6 on its cognate receptor, NKp30, could be the most effective strategy.

FUNDING

This work was supported by the National Research Council of Science & Technology (NST) grant (CAP-18-02-KRIBB, GTL24021-000), a National Research Foundation grant (2710012258, 2710004815), and an Institute for Basic Science grant (IBS-R029-C3).

Self-organized pattern formation in the developing mouse neural tube by a temporal relay of BMP signaling

Developing tissues interpret dynamic changes in morphogen activity to generate cell type diversity. To quantitatively study bone morphogenetic protein (BMP) signaling dynamics in the mouse neural tube, we developed an embryonic stem cell differentiation system tailored for growing tissues. Differentiating cells form striking self-organized patterns of dorsal neural tube cell types driven by sequential phases of BMP signaling that are observed both in vitro and in vivo. Data-driven biophysical modeling showed that these dynamics result from coupling fast negative feedback with slow positive regulation of signaling by the specification of an endogenous BMP source. Thus, in contrast to relays that propagate morphogen signaling in space, we identify a BMP signaling relay that operates in time. This mechanism allows for a rapid initial concentration-sensitive response that is robustly terminated, thereby regulating balanced sequential cell type generation. Our study provides an experimental and theoretical framework to understand how signaling dynamics are exploited in developing tissues.

EasyFlow: An open-source, user-friendly cytometry analyzer with graphic user interface (GUI)

Flow cytometry enables quantitative measurements of fluorescence in single cells. The technique was widely used for immunology to identify populations with different surface protein markers. More recently, the usage of flow cytometry has been extended to additional readouts, including intracellular proteins and fluorescent protein transgenes, and is widely utilized to study developmental biology, systems biology, microbiology, and many other fields. A common file format (FCS format, defined by the International Society for Advancement of Cytometry (ISAC)) has been universally adopted, facilitating data exchange between different machines. A diverse spectrum of software packages has been developed for the analysis of flow cytometry data. However, those are either 1) costly proprietary softwares, 2) open source packages with prerequisite installation of R or Python and sometimes require users to have experience in coding, or 3) online tools that are limiting for analysis of large data sets. Here, we present EasyFlow, an open-source flow cytometry analysis graphic user interface (GUI) based on Matlab or Python, that can be installed and run locally across platforms (Windows, MacOS, and Linux) without requiring previous coding knowledge. The Python version (EasyFlowQ) is also developed on a popular plotting framework (Matplotlib) and modern user interface toolkit (Qt), allowing more advanced users to customize and keep contributing to the software, as well as its tutorials. Overall, EasyFlow serves as a simple-to-use tool for inexperienced users with little coding experience to use locally, as well as a platform for advanced users to further customize for their own needs.

Systems profiling reveals recurrently dysregulated cytokine signaling responses in ER+ breast cancer patients' blood

Cytokines operate in concert to maintain immune homeostasis and coordinate immune responses. In cases of ER+ breast cancer, peripheral immune cells exhibit altered responses to several cytokines, and these alterations are correlated strongly with patient outcomes. To develop a systems-level understanding of this dysregulation, we measured a panel of cytokine responses and receptor abundances in the peripheral blood of healthy controls and ER+ breast cancer patients across immune cell types. Using tensor factorization to model this multidimensional data, we found that breast cancer patients exhibited widespread alterations in response, including drastically reduced response to IL-10 and heightened basal levels of pSmad2/3 and pSTAT4. ER+ patients also featured upregulation of PD-L1, IL6Rα, and IL2Rα, among other receptors. Despite this, alterations in response to cytokines were not explained by changes in receptor abundances. Thus, tensor factorization helped to reveal a coordinated reprogramming of the immune system that was consistent across our cohort.

Sulforaphane Promotes Proliferation of Porcine Granulosa Cells via the H3K27ac-Mediated GDF8-ALK5-ERK Pathway

Follicle development, a crucial process in reproductive biology, hinges upon the dynamic proliferation of granulosa cells (GCs). Growth differentiation factor-8 (GDF8) is well-known as myostatin for inhibiting skeletal muscle growth, and it also exists in ovarian GCs and follicle fluid. However, the relationship between GCs proliferation and GDF8 remains elusive. Sulforaphane (SFN) is a potent bioactive compound, which in our study has been demonstrated to induce the expression of GDF8 in GCs. Meanwhile, we discover a novel role of SFN in promoting the proliferation of porcine GCs. Specifically, SFN enhances GCs proliferation by accelerating the progression of the cell cycle through the G1 phase to the S phase. By performing gene expression profiling, we showed that the promoting proliferative effects of SFN are highly correlated with the TGF-β signaling pathways and cell cycle. Among the ligand factors of TGF-β signaling, we identify GDF8 as a critical downstream effector of SFN, which acts through ALK5 to mediate SFN-induced proliferation and G1/S transition. In addition, we identify a noncanonical downstream pathway by which GDF8 induces the activation of MAPK/ERK to facilitate the cell cycle progression in GCs. Moreover, we reveal that the expression of GDF8 is regulated by SFN through epigenetic modifications of H3K27 acetylation. These findings not only provide mechanistic insights into the regulation of GCs proliferation but also establish a previously unrecognized role of GDF8 in follicle development, which have significant implications for developing strategies to improve female fertility.

Physical extraction of antigen and information

To respond and adapt, cells use surface receptors to sense environmental cues. While biochemical signal processing inside the cell is studied in depth, less is known about how physical processes during cell-cell contact impact signal acquisition. New experiments found that fast-evolving immune B cells in germinal centers (GCs) apply force to acquire antigen clusters prior to internalization, suggesting adaptive benefits of physical information extraction. We present a theory of stochastic antigen transfer and show that maximizing information gain via physical extraction can explain the dramatic phenotypic transition from naive to GC B cells-attenuated receptor signaling, enhanced force usage, and decentralized contact architecture. Our model suggests that binding-lifetime measurement and physical extraction serve as complementary modes of antigen recognition, greatly extending the dynamic range of affinity discrimination when combined. This physical-information framework further predicts that the optimal size of receptor clusters decreases as affinity improves, rationalizing the use of a multifocal synaptic pattern seen in GC B cells. By linking extraction dynamics to selection fidelity via discriminatory performance, we propose that cells may physically enhance information acquisition to sustain adaptive evolution.

Principles of Computation by Competitive Protein Dimerization Networks

Many biological signaling pathways employ proteins that competitively dimerize in diverse combinations. These dimerization networks can perform biochemical computations, in which the concentrations of monomers (inputs) determine the concentrations of dimers (outputs). Despite their prevalence, little is known about the range of input-output computations that dimerization networks can perform (their "expressivity") and how it depends on network size and connectivity. Using a systematic computational approach, we demonstrate that even small dimerization networks (3-6 monomers) are expressive, performing diverse multi-input computations. Further, dimerization networks are versatile, performing different computations when their protein components are expressed at different levels, such as in different cell types. Remarkably, individual networks with random interaction affinities, when large enough (≥8 proteins), can perform nearly all (~90%) potential one-input network computations merely by tuning their monomer expression levels. Thus, even the simple process of competitive dimerization provides a powerful architecture for multi-input, cell-type-specific signal processing.

BMP10 accelerated spinal astrocytic activation in neuropathic pain via ALK2/smad1/5/8 signaling

Background

Astrocytic activation in the spinal dorsal horn contributes to the central sensitization of neuropathic pain. Bone morphogenetic protein (BMP) 10, one of the BMPs highly expressed in the central nervous system, has been demonstrated to have an accelerated effect on astrocytic activation. This study aimed to investigate the functional effects of BMP10 on the activation of astrocytes in the spinal dorsal horn of animal model of neuropathic pain and to explore potential mechanisms involved in this process.

Methods

A neuropathic pain mice model was established using the spared nerve injury (SNI). Western blot analysis was performed to detect the expressional levels of BMP10, activin receptor-like receptor 2 (ALK2), Smad1/5/8, phosphorylated Smad1/5/8, and glial fibrillary acidic protein (GFAP). Immunofluorescence staining was used to detect BMP10, ALK2, and GFAP distribution and expression. The behavioral changes in mice were evaluated using paw withdrawal threshold (PWT), thermal withdrawal latency (TWL), and open field test (OFT). The BMP10 siRNA, Smad1 siRNA, BMP10 peptide, and ALK2-IN-2 (ALK2 inhibitor) were intrathecally administrated to mice. A model of lipopolysaccharide (LPS)-stimulated astrocytes was established to investigate the effect of Smad1. The transfection efficiency of siRNAs was detected by western blot and qRT-PCR analysis.

Results

BMP10 levels were increased in the L4-6 ipsilateral spinal dorsal horn of SNI mice and particularly elevated in astrocytes. Consistently, GFAP and phosphorylated Smad1/5/8 were upregulated in the L4-6 ipsilateral spinal dorsal horn after SNI, indicating the activation of astrocytes and Smad1/5/8 signaling. An intrathecal injection of BMP10 siRNA abrogated pain hypersensitivity and astrocytic activation in SNI mice. In addition, intrathecal administration of BMP10 peptide evoked pain hypersensitivity and astrocytic activation in normal mice, and this action was reversed by inhibiting the ALK2. Furthermore, targeting Smad1 in vitro with the help of siRNA inhibited the activation of astrocytes induced by LPS. Finally, targeting Smad1 abrogated BMP10-induced hypersensitivity and activation of astrocytes.

Conclusion

These findings indicate that the BMP10/ALK2/Smad1/5/8 axis plays a key role in pain hypersensitivity after peripheral nerve injury, which indicates its stimulative ability toward astrocytes.

Reducing Uncertainty Through Mutual Information in Structural and Systems Biology

Systems biology models are useful models of complex biological systems that may require a large amount of experimental data to fit each model's parameters or to approximate a likelihood function. These models range from a few to thousands of parameters depending on the complexity of the biological system modeled, potentially making the task of fitting parameters to the model difficult - especially when new experimental data cannot be gathered. We demonstrate a method that uses structural biology predictions to augment systems biology models to improve systems biology models' predictions without having to gather more experimental data. Additionally, we show how systems biology models' predictions can help evaluate novel structural biology hypotheses, which may also be expensive or infeasible to validate.

Intercellular Communication in Shoot Meristems

The shoot meristem of land plants maintains the capacity for organ generation throughout its lifespan due to a group of undifferentiated stem cells. Most meristems are shaped like a dome with a precise spatial arrangement of functional domains, and, within and between these domains, cells interact through a network of interconnected signaling pathways. Intercellular communication in meristems is mediated by mobile transcription factors, small RNAs, hormones, and secreted peptides that are perceived by membrane-localized receptors. In recent years, we have gained deeper insight into the underlying molecular processes of the shoot meristem, and we discuss here how plants integrate internal and external inputs to control shoot meristem activities.

Protein Design Using Structure-Prediction Networks: AlphaFold and RoseTTAFold as Protein Structure Foundation Models

Designing proteins with tailored structures and functions is a long-standing goal in bioengineering. Recently, deep learning advances have enabled protein structure prediction at near-experimental accuracy, which has catalyzed progress in protein design as well. We review recent studies that use structure-prediction neural networks to design proteins, via approaches such as activation maximization, inpainting, or denoising diffusion. These methods have led to major improvements over previous methods in wet-lab success rates for designing protein binders, metalloproteins, enzymes, and oligomeric assemblies. These results show that structure-prediction models are a powerful foundation for developing protein-design tools and suggest that continued improvement of their accuracy and generality will be key to unlocking the full potential of protein design.

Silibinin is a suppressor of the metastasis-promoting transcription factor ID3

BACKGROUND

ID3 (inhibitor of DNA binding/differentiation-3) is a transcription factor that enables metastasis by promoting stem cell-like properties in endothelial and tumor cells. The milk thistle flavonolignan silibinin is a phytochemical with anti-metastatic potential through largely unknown mechanisms.

HYPOTHESIS/PURPOSE

We have mechanistically investigated the ability of silibinin to inhibit the aberrant activation of ID3 in brain endothelium and non-small cell lung cancer (NSCLC) models.

METHODS

Bioinformatic analyses were performed to investigate the co-expression correlation between ID3 and bone morphogenic protein (BMP) ligands/BMP receptors (BMPRs) genes in NSCLC patient datasets. ID3 expression was assessed by immunoblotting and qRT-PCR. Luciferase reporter assays were used to evaluate the gene sequences targeted by silibinin to regulate ID3 transcription. In silico computational modeling and LanthaScreen TR-FRET kinase assays were used to characterize and validate the BMPR inhibitory activity of silibinin. Tumor tissues from NSCLC xenograft models treated with oral silibinin were used to evaluate the in vivo anti-ID3 effects of silibinin.

RESULTS

Analysis of lung cancer patient datasets revealed a top-ranked positive association of ID3 with the BMP9 endothelial receptor ACVRL1/ALK1 and the BMP ligand BMP6. Silibinin treatment blocked the BMP9-induced activation of the ALK1-phospho-SMAD1/5-ID3 axis in brain endothelial cells. Constitutive, acquired, and adaptive expression of ID3 in NSCLC cells were all significantly downregulated in response to silibinin. Silibinin blocked ID3 transcription via BMP-responsive elements in ID3 gene enhancers. Silibinin inhibited the kinase activities of BMPRs in the micromolar range, with the lower IC50 values occurring against ACVRL1/ALK1 and BMPR2. In an in vivo NSCLC xenograft model, tumoral overexpression of ID3 was completely suppressed by systematically achievable oral doses of silibinin.

CONCLUSIONS

ID3 is a largely undruggable metastasis-promoting transcription factor. Silibinin is a novel suppressor of ID3 that may be explored as a novel therapeutic approach to interfere with the metastatic dissemination capacity of NSCLC.

Liebig's law of the minimum in the TGF-β/SMAD pathway

Cells use signaling pathways to sense and respond to their environments. The transforming growth factor-β (TGF-β) pathway produces context-specific responses. Here, we combined modeling and experimental analysis to study the dependence of the output of the TGF-β pathway on the abundance of signaling molecules in the pathway. We showed that the TGF-β pathway processes the variation of TGF-β receptor abundance using Liebig's law of the minimum, meaning that the output-modifying factor is the signaling protein that is most limited, to determine signaling responses across cell types and in single cells. We found that the abundance of either the type I (TGFBR1) or type II (TGFBR2) TGF-β receptor determined the responses of cancer cell lines, such that the receptor with relatively low abundance dictates the response. Furthermore, nuclear SMAD2 signaling correlated with the abundance of TGF-β receptor in single cells depending on the relative expression levels of TGFBR1 and TGFBR2. A similar control principle could govern the heterogeneity of signaling responses in other signaling pathways.

Cellular and Molecular Mechanisms of Heterotopic Ossification in Fibrodysplasia Ossificans Progressiva

Fibrodysplasia ossificans progressiva (FOP) is a debilitating genetic disorder characterized by recurrent episodes of heterotopic ossification (HO) formation in muscles, tendons, and ligaments. FOP is caused by a missense mutation in the ACVR1 gene (activin A receptor type I), an important signaling receptor involved in endochondral ossification. The ACVR1R206H mutation induces increased downstream canonical SMAD-signaling and drives tissue-resident progenitor cells with osteogenic potential to participate in endochondral HO formation. In this article, we review aberrant ACVR1R206H signaling and the cells that give rise to HO in FOP. FOP mouse models and lineage tracing analyses have been used to provide strong evidence for tissue-resident mesenchymal cells as cellular contributors to HO. We assess how the underlying mutation in FOP disrupts muscle-specific dynamics during homeostasis and repair, with a focus on muscle-resident mesenchymal cells known as fibro-adipogenic progenitors (FAPs). Accumulating research points to FAPs as a prominent HO progenitor population, with ACVR1R206H FAPs not only aberrantly differentiating into chondro-osteogenic lineages but creating a permissive environment for bone formation at the expense of muscle regeneration. We will further discuss the emerging role of ACVR1R206H FAPs in muscle regeneration and therapeutic targeting of these cells to reduce HO formation in FOP.

Opsin-Free Activation of Bmp Receptors by a Femtosecond Laser

Bone morphogenetic protein (BMP) signaling plays a vital role in differentiation, organogenesis, and various cell processes. As a member of TGF-β superfamily, the BMP initiation usually accompanies crosstalk with other signaling pathways and simultaneously activates some of them. It is quite challenging to solely initiate an individual pathway. In this study, an opsin-free optical method to specifically activate BMP receptors (BMPR) and subsequent pSmad1/5/8 cascades by a single-time scan of a tightly-focused femtosecond laser in the near infrared range is reported. Via transient two-photon excitation to intrinsic local flavins near the cell membrane, the photoactivation drives conformational changes of preformed BMPR complexes to enable their bonding and phosphorylation of the GS domain in BMPR-I by BMPR-II. The pSmad1/5/8 signaling is initiated by this method, while p38 and pSmad2 are rarely perturbed. Based on a microscopic system, primary adipose-derived stem cells in an area of 420  × 420 µm2 are photoactivated by a single-time laser scanning for 1.5 s and exhibit pSmad1/5/8 upregulation and osteoblastic differentiation after 21 days. Hence, an opsin-free, specific, and noninvasive optical method to initiate BMP signaling, easily accomplished by a two-photon microscope system is reported.

Enrichment on steps, not genes, improves inference of differentially expressed pathways

Enrichment analysis is frequently used in combination with differential expression data to investigate potential commonalities amongst lists of genes and generate hypotheses for further experiments. However, current enrichment analysis approaches on pathways ignore the functional relationships between genes in a pathway, particularly OR logic that occurs when a set of proteins can each individually perform the same step in a pathway. As a result, these approaches miss pathways with large or multiple sets because of an inflation of pathway size (when measured as the total gene count) relative to the number of steps. We address this problem by enriching on step-enabling entities in pathways. We treat sets of protein-coding genes as single entities, and we also weight sets to account for the number of genes in them using the multivariate Fisher's noncentral hypergeometric distribution. We then show three examples of pathways that are recovered with this method and find that the results have significant proportions of pathways not found in gene list enrichment analysis.

Time-integrated BMP signaling determines fate in a stem cell model for early human development

How paracrine signals are interpreted to yield multiple cell fate decisions in a dynamic context during human development in vivo and in vitro remains poorly understood. Here we report an automated tracking method to follow signaling histories linked to cell fate in large numbers of human pluripotent stem cells (hPSCs). Using an unbiased statistical approach, we discover that measured BMP signaling history correlates strongly with fate in individual cells. We find that BMP response in hPSCs varies more strongly in the duration of signaling than the level. However, both the level and duration of signaling activity control cell fate choices only by changing the time integral. Therefore, signaling duration and level are interchangeable in this context. In a stem cell model for patterning of the human embryo, we show that signaling histories predict the fate pattern and that the integral model correctly predicts changes in cell fate domains when signaling is perturbed. Our data suggest that mechanistically, BMP signaling is integrated by SOX2.

Logical design of synthetic cis-regulatory DNA for genetic tracing of cell identities and state changes

Descriptive data are rapidly expanding in biomedical research. Instead, functional validation methods with sufficient complexity remain underdeveloped. Transcriptional reporters allow experimental characterization and manipulation of developmental and disease cell states, but their design lacks flexibility. Here, we report logical design of synthetic cis-regulatory DNA (LSD), a computational framework leveraging phenotypic biomarkers and trans-regulatory networks as input to design reporters marking the activity of selected cellular states and pathways. LSD uses bulk or single-cell biomarkers and a reference genome or custom cis-regulatory DNA datasets with user-defined boundary regions. By benchmarking validated reporters, we integrate LSD with a computational ranking of phenotypic specificity of putative cis-regulatory DNA. Experimentally, LSD-designed reporters targeting a wide range of cell states are functional without minimal promoters. Applied to broadly expressed genes from human and mouse tissues, LSD generates functional housekeeper-like sLCRs compatible with size constraints of AAV vectors for gene therapy applications. A mesenchymal glioblastoma reporter designed by LSD outperforms previously validated ones and canonical cell surface markers. In genome-scale CRISPRa screens, LSD facilitates the discovery of known and novel bona fide cell-state drivers. Thus, LSD captures core principles of cis-regulation and is broadly applicable to studying complex cell states and mechanisms of transcriptional regulation.

The roles and regulatory mechanisms of TGF-β and BMP signaling in bone and cartilage development, homeostasis and disease

Transforming growth factor-βs (TGF-βs) and bone morphometric proteins (BMPs) belong to the TGF-β superfamily and perform essential functions during osteoblast and chondrocyte lineage commitment and differentiation, skeletal development, and homeostasis. TGF-βs and BMPs transduce signals through SMAD-dependent and -independent pathways; specifically, they recruit different receptor heterotetramers and R-Smad complexes, resulting in unique biological readouts. BMPs promote osteogenesis, osteoclastogenesis, and chondrogenesis at all differentiation stages, while TGF-βs play different roles in a stage-dependent manner. BMPs and TGF-β have opposite functions in articular cartilage homeostasis. Moreover, TGF-β has a specific role in maintaining the osteocyte network. The precise activation of BMP and TGF-β signaling requires regulatory machinery at multiple levels, including latency control in the matrix, extracellular antagonists, ubiquitination and phosphorylation in the cytoplasm, nucleus-cytoplasm transportation, and transcriptional co-regulation in the nuclei. This review weaves the background information with the latest advances in the signaling facilitated by TGF-βs and BMPs, and the advanced understanding of their diverse physiological functions and regulations. This review also summarizes the human diseases and mouse models associated with disordered TGF-β and BMP signaling. A more precise understanding of the BMP and TGF-β signaling could facilitate the development of bona fide clinical applications in treating bone and cartilage disorders.

Combinatorial expression motifs in signaling pathways

In animal cells, molecular pathways often comprise families of variant components, such as ligands or receptors. These pathway components are differentially expressed by different cell types, potentially tailoring pathway function to cell context. However, it has remained unclear how pathway expression profiles are distributed across cell types and whether similar profiles can occur in dissimilar cell types. Here, using single-cell gene expression datasets, we identified pathway expression motifs, defined as recurrent expression profiles that are broadly distributed across diverse cell types. Motifs appeared in core pathways, including TGF-β, Notch, Wnt, and the SRSF splice factors, and involved combinatorial co-expression of multiple components. Motif usage was weakly correlated between pathways in adult cell types and during dynamic developmental transitions. Together, these results suggest a mosaic view of cell type organization, in which different cell types operate many of the same pathways in distinct modes.

Self-organized BMP signaling dynamics underlie the development and evolution of digit segmentation patterns in birds and mammals

Repeating patterns of synovial joints are a highly conserved feature of articulated digits, with variations in joint number and location resulting in diverse digit morphologies and limb functions across the tetrapod clade. During the development of the amniote limb, joints form iteratively within the growing digit ray, as a population of distal progenitors alternately specifies joint and phalanx cell fates to segment the digit into distinct elements. While numerous molecular pathways have been implicated in this fate choice, it remains unclear how they give rise to a repeating pattern. Here, using single-cell RNA sequencing and spatial gene expression profiling, we investigate the transcriptional dynamics of interphalangeal joint specification in vivo. Combined with mathematical modeling, we predict that interactions within the BMP signaling pathway-between the ligand GDF5, the inhibitor NOGGIN, and the intracellular effector pSMAD-result in a self-organizing Turing system that forms periodic joint patterns. Our model is able to recapitulate the spatiotemporal gene expression dynamics observed in vivo, as well as phenocopy digit malformations caused by BMP pathway perturbations. By contrasting in silico simulations with in vivo morphometrics of two morphologically distinct digits, we show how changes in signaling parameters and growth dynamics can result in variations in the size and number of phalanges. Together, our results reveal a self-organizing mechanism that underpins amniote digit segmentation and its evolvability and, more broadly, illustrate how Turing systems based on a single molecular pathway may generate complex repetitive patterns in a wide variety of organisms.

Ischemia-reperfusion injury: molecular mechanisms and therapeutic targets

Ischemia-reperfusion (I/R) injury paradoxically occurs during reperfusion following ischemia, exacerbating the initial tissue damage. The limited understanding of the intricate mechanisms underlying I/R injury hinders the development of effective therapeutic interventions. The Wnt signaling pathway exhibits extensive crosstalk with various other pathways, forming a network system of signaling pathways involved in I/R injury. This review article elucidates the underlying mechanisms involved in Wnt signaling, as well as the complex interplay between Wnt and other pathways, including Notch, phosphatidylinositol 3-kinase/protein kinase B, transforming growth factor-β, nuclear factor kappa, bone morphogenetic protein, N-methyl-D-aspartic acid receptor-Ca2+-Activin A, Hippo-Yes-associated protein, toll-like receptor 4/toll-interleukine-1 receptor domain-containing adapter-inducing interferon-β, and hepatocyte growth factor/mesenchymal-epithelial transition factor. In particular, we delve into their respective contributions to key pathological processes, including apoptosis, the inflammatory response, oxidative stress, extracellular matrix remodeling, angiogenesis, cell hypertrophy, fibrosis, ferroptosis, neurogenesis, and blood-brain barrier damage during I/R injury. Our comprehensive analysis of the mechanisms involved in Wnt signaling during I/R reveals that activation of the canonical Wnt pathway promotes organ recovery, while activation of the non-canonical Wnt pathways exacerbates injury. Moreover, we explore novel therapeutic approaches based on these mechanistic findings, incorporating evidence from animal experiments, current standards, and clinical trials. The objective of this review is to provide deeper insights into the roles of Wnt and its crosstalk signaling pathways in I/R-mediated processes and organ dysfunction, to facilitate the development of innovative therapeutic agents for I/R injury.

Context-dependent TGFβ family signalling in cell fate regulation

The transforming growth factor-β (TGFβ) family are a large group of evolutionarily conserved cytokines whose signalling modulates cell fate decision-making across varying cellular contexts at different stages of life. Here we discuss new findings in early embryos that reveal how, in contrast to our original understanding of morphogen interpretation, robust cell fate specification can originate from a noisy combination of signalling inputs and a broad range of signalling levels. We compare this evidence with novel findings on the roles of TGFβ family signalling in tissue maintenance and homeostasis during juvenile and adult life, spanning the skeletal, haemopoietic and immune systems. From these comparisons, it emerges that in contrast to robust developing systems, relatively small perturbations in TGFβ family signalling have detrimental effects at later stages in life, leading to aberrant cell fate specification and disease, for example in cancer or congenital disorders. Finally, we highlight novel strategies to target and amend dysfunction in signalling and discuss how gleaning knowledge from different fields of biology can help in the development of therapeutics for aberrant TGFβ family signalling in disease.

Systems profiling reveals recurrently dysregulated cytokine signaling responses in ER+ breast cancer patients' blood

Cytokines mediate cell-to-cell communication across the immune system and therefore are critical to immunosurveillance in cancer and other diseases. Several cytokines show dysregulated abundance or signaling responses in breast cancer, associated with the disease and differences in survival and progression. Cytokines operate in a coordinated manner to affect immune surveillance and regulate one another, necessitating a systems approach for a complete picture of this dysregulation. Here, we profiled cytokine signaling responses of peripheral immune cells from breast cancer patients as compared to healthy controls in a multidimensional manner across ligands, cell populations, and responsive pathways. We find alterations in cytokine responsiveness across pathways and cell types that are best defined by integrated signatures across dimensions. Alterations in the abundance of a cytokine's cognate receptor do not explain differences in responsiveness. Rather, alterations in baseline signaling and receptor abundance suggesting immune cell reprogramming are associated with altered responses. These integrated features suggest a global reprogramming of immune cell communication in breast cancer.

Multivalent, asymmetric IL-2-Fc fusions show enhanced selectivity for regulatory T cells

The cytokine interleukin-2 (IL-2) has the potential to treat autoimmune disease but is limited by its modest specificity toward immunosuppressive regulatory T (Treg) cells. IL-2 receptors consist of combinations of α, β, and γ chains of variable affinity and cell specificity. Engineering IL-2 to treat autoimmunity has primarily focused on retaining binding to the relatively Treg-selective, high-affinity receptor while reducing binding to the less selective, low-affinity receptor. However, we found that refining the designs to focus on targeting the high-affinity receptor through avidity effects is key to optimizing Treg selectivity. We profiled the dynamics and dose dependency of signaling responses in primary human immune cells induced by engineered fusions composed of either wild-type IL-2 or mutant forms with altered affinity, valency, and fusion to the antibody Fc region for stability. Treg selectivity and signaling response variations were explained by a model of multivalent binding and dimer-enhanced avidity-a combined measure of the strength, number, and conformation of interaction sites-from which we designed tetravalent IL-2-Fc fusions that had greater Treg selectivity in culture than do current designs. Biasing avidity toward IL2Rα with an asymmetrical multivalent design consisting of one α/β chain-binding and one α chain-binding mutant further enhanced Treg selectivity. Comparative analysis revealed that IL2Rα was the optimal cell surface target for Treg selectivity, indicating that avidity for IL2Rα may be the optimal route to producing IL-2 variants that selectively target Tregs.

Factorization and pseudofactorization of weighted graphs

For unweighted graphs, finding isometric embeddings of a graph G is closely related to decompositions of G into Cartesian products of smaller graphs. When G is isomorphic to a Cartesian graph product, we call the factors of this product a factorization of G . When G is isomorphic to an isometric subgraph of a Cartesian graph product, we call those factors a pseudofactorization of G . Prior work has shown that an unweighted graph's pseudofactorization can be used to generate a canonical isometric embedding into a product of the smallest possible pseudofactors. However, for arbitrary weighted graphs, which represent a richer variety of metric spaces, methods for finding isometric embeddings or determining their existence remain elusive, and indeed pseudofactorization and factorization have not previously been extended to this context. In this work, we address the problem of finding the factorization and pseudofactorization of a weighted graph G , where G satisfies the property that every edge constitutes a shortest path between its endpoints. We term such graphs minimal graphs, noting that every graph can be made minimal by removing edges not affecting its path metric. We generalize pseudofactorization and factorization to minimal graphs and develop new proof techniques that extend the previously proposed algorithms due to Graham and Winkler [Graham and Winkler, '85] and Feder [Feder, '92] for pseudofactorization and factorization of unweighted graphs. We show that any n -vertex, m -edge graph with positive integer edge weights can be factored in O ( m 2 ) time, plus the time to find all pairs shortest paths (APSP) distances in a weighted graph, resulting in an overall running time of O ( m 2 + n 2  log log  n ) time. We also show that a pseudofactorization for such a graph can be computed in O ( m n ) time, plus the time to solve APSP, resulting in an O ( m n + n 2  log log  n ) running time.

TGF-β signaling in health and disease

The TGF-β regulatory system plays crucial roles in the preservation of organismal integrity. TGF-β signaling controls metazoan embryo development, tissue homeostasis, and injury repair through coordinated effects on cell proliferation, phenotypic plasticity, migration, metabolic adaptation, and immune surveillance of multiple cell types in shared ecosystems. Defects of TGF-β signaling, particularly in epithelial cells, tissue fibroblasts, and immune cells, disrupt immune tolerance, promote inflammation, underlie the pathogenesis of fibrosis and cancer, and contribute to the resistance of these diseases to treatment. Here, we review how TGF-β coordinates multicellular response programs in health and disease and how this knowledge can be leveraged to develop treatments for diseases of the TGF-β system.

Terminal differentiation of villus tip enterocytes is governed by distinct Tgfβ superfamily members

The protective and absorptive functions of the intestinal epithelium rely on differentiated enterocytes in the villi. The differentiation of enterocytes is orchestrated by sub-epithelial mesenchymal cells producing distinct ligands along the villus axis, in particular Bmps and Tgfβ. Here, we show that individual Bmp ligands and Tgfβ drive distinct enterocytic programs specific to villus zonation. Bmp4 is expressed from the centre to the upper part of the villus and activates preferentially genes connected to lipid uptake and metabolism. In contrast, Bmp2 is produced by villus tip mesenchymal cells and it influences the adhesive properties of villus tip epithelial cells and the expression of immunomodulators. Additionally, Tgfβ induces epithelial gene expression programs similar to those triggered by Bmp2. Bmp2-driven villus tip program is activated by a canonical Bmp receptor type I/Smad-dependent mechanism. Finally, we establish an organoid cultivation system that enriches villus tip enterocytes and thereby better mimics the cellular composition of the intestinal epithelium. Our data suggest that not only a Bmp gradient but also the activity of individual Bmp drives specific enterocytic programs.

Reduced GS Domain Serine/Threonine Requirements of Fibrodysplasia Ossificans Progressiva Mutant Type I BMP Receptor ACVR1 in the Zebrafish

Fibrodysplasia ossificans progressiva (FOP) is a rare human genetic condition characterized by altered skeletal development and extraskeletal bone formation. All cases of FOP are caused by mutations in the type I bone morphogenetic protein (BMP) receptor gene ACVR1 that result in overactivation of the BMP signaling pathway. Activation of the wild-type ACVR1 kinase requires assembly of a tetrameric type I and II BMP receptor complex followed by phosphorylation of the ACVR1 GS domain by type II BMP receptors. Previous studies showed that the FOP-mutant ACVR1-R206H required type II BMP receptors and presumptive glycine/serine-rich (GS) domain phosphorylation for overactive signaling. Structural modeling of the ACVR1-R206H mutant kinase domain supports the idea that FOP mutations alter the conformation of the GS domain, but it is unclear how this leads to overactive signaling. Here we show, using a developing zebrafish embryo BMP signaling assay, that the FOP-mutant receptors ACVR1-R206H and -G328R have reduced requirements for GS domain phosphorylatable sites to signal compared to wild-type ACVR1. Further, ligand-independent and ligand-dependent signaling through the FOP-mutant ACVR1 receptors have distinct GS domain phosphorylatable site requirements. ACVR1-G328R showed increased GS domain serine/threonine requirements for ligand-independent signaling compared to ACVR1-R206H, whereas it exhibited reduced serine/threonine requirements for ligand-dependent signaling. Remarkably, while ACVR1-R206H does not require the type I BMP receptor partner, Bmpr1, to signal, a ligand-dependent GS domain mutant of ACVR1-R206H could signal independently of Bmpr1 only when Bmp7 ligand was overexpressed. Of note, unlike human ACVR1-R206H, the zebrafish paralog Acvr1l-R203H does not show increased signaling activity. However, in domain-swapping studies, the human kinase domain, but not the human GS domain, was sufficient to confer overactive signaling to the Acvr1l-R203H receptor. Together these results reflect the importance of GS domain activation and kinase domain functions in regulating ACVR1 signaling and identify mechanisms of reduced regulatory constraints conferred by FOP mutations. © 2023 American Society for Bone and Mineral Research (ASBMR).

Future medicine: from molecular pathways to the collective intelligence of the body

The remarkable anatomical homeostasis exhibited by complex living organisms suggests that they are inherently reprogrammable information-processing systems that offer numerous interfaces to their physiological and anatomical problem-solving capacities. We briefly review data suggesting that the multiscale competency of living forms affords a new path for biomedicine that exploits the innate collective intelligence of tissues and organs. The concept of tissue-level allostatic goal-directedness is already bearing fruit in clinical practice. We sketch a roadmap towards 'somatic psychiatry' by using advances in bioelectricity and behavioral neuroscience to design methods that induce self-repair of structure and function. Relaxing the assumption that cellular control mechanisms are static, exploiting powerful concepts from cybernetics, behavioral science, and developmental biology may spark definitive solutions to current biomedical challenges.

Reconstitution of morphogen shuttling circuits

Developing tissues form spatial patterns by establishing concentration gradients of diffusible signaling proteins called morphogens. The bone morphogenetic protein (BMP) morphogen pathway uses a family of extracellular modulators to reshape signaling gradients by actively "shuttling" ligands to different locations. It has remained unclear what circuits are sufficient to enable shuttling, what other patterns they can generate, and whether shuttling is evolutionarily conserved. Here, using a synthetic, bottom-up approach, we compared the spatiotemporal dynamics of different extracellular circuits. Three proteins-Chordin, Twsg, and the BMP-1 protease-successfully displaced gradients by shuttling ligands away from the site of production. A mathematical model explained the different spatial dynamics of this and other circuits. Last, combining mammalian and Drosophila components in the same system suggests that shuttling is a conserved capability. Together, these results reveal principles through which extracellular circuits control the spatiotemporal dynamics of morphogen signaling.

The computational capabilities of many-to-many protein interaction networks

Many biological circuits comprise sets of protein variants that interact with one another in a many-to-many, or promiscuous, fashion. These architectures can provide powerful computational capabilities that are especially critical in multicellular organisms. Understanding the principles of biochemical computations in these circuits could allow more precise control of cellular behaviors. However, these systems are inherently difficult to analyze, due to their large number of interacting molecular components, partial redundancies, and cell context dependence. Here, we discuss recent experimental and theoretical advances that are beginning to reveal how promiscuous circuits compute, what roles those computations play in natural biological contexts, and how promiscuous architectures can be applied for the design of synthetic multicellular behaviors.

Isometric Hamming embeddings of weighted graphs

A mapping α : V ( G ) → V ( H ) from the vertex set of one graph G to another graph H is an isometric embedding if the shortest path distance between any two vertices in G equals the distance between their images in H . Here, we consider isometric embeddings of a weighted graph G into unweighted Hamming graphs, called Hamming embeddings, when G satisfies the property that every edge is a shortest path between its endpoints. Using a Cartesian product decomposition of G called its canonical isometric representation, we show that every Hamming embedding of G may be partitioned into a canonical partition, whose parts provide Hamming embeddings for each factor of the canonical isometric representation of G . This implies that G permits a Hamming embedding if and only if each factor of its canonical isometric representation is Hamming embeddable. This result extends prior work on unweighted graphs that showed that an unweighted graph permits a Hamming embedding if and only if each factor is a complete graph. When a graph G has nontrivial isometric representation, determining whether G has a Hamming embedding can be simplified to checking embeddability of two or more smaller graphs.

Subfunctionalized expression drives evolutionary retention of ribosomal protein paralogs Rps27 and Rps27l in vertebrates

The formation of paralogs through gene duplication is a core evolutionary process. For paralogs that encode components of protein complexes such as the ribosome, a central question is whether they encode functionally distinct proteins or whether they exist to maintain appropriate total expression of equivalent proteins. Here, we systematically tested evolutionary models of paralog function using the ribosomal protein paralogs Rps27 (eS27) and Rps27l (eS27L) as a case study. Evolutionary analysis suggests that Rps27 and Rps27l likely arose during whole-genome duplication(s) in a common vertebrate ancestor. We show that Rps27 and Rps27l have inversely correlated mRNA abundance across mouse cell types, with the highest Rps27 in lymphocytes and the highest Rps27l in mammary alveolar cells and hepatocytes. By endogenously tagging the Rps27 and Rps27l proteins, we demonstrate that Rps27- and Rps27l-ribosomes associate preferentially with different transcripts. Furthermore, murine Rps27 and Rps27l loss-of-function alleles are homozygous lethal at different developmental stages. However, strikingly, expressing Rps27 protein from the endogenous Rps27l locus or vice versa completely rescues loss-of-function lethality and yields mice with no detectable deficits. Together, these findings suggest that Rps27 and Rps27l are evolutionarily retained because their subfunctionalized expression patterns render both genes necessary to achieve the requisite total expression of two equivalent proteins across cell types. Our work represents the most in-depth characterization of a mammalian ribosomal protein paralog to date and highlights the importance of considering both protein function and expression when investigating paralogs.

Structural insights into regulation of CCN protein activities and functions

CCN proteins play important functions during development, in repair mechanisms following tissue injury, as well as in pathophysiologic mechanisms of metastasis of cancer. CCNs are secreted proteins that have a multimodular structure and are categorized as matricellular proteins. Although the prevailing view is that CCN proteins regulate biologic processes by interacting with a wide array of other proteins in the microenvironment of the extracellular matrix, the molecular mechanisms of action of CCN proteins are still poorly understood. Not dissuading the current view, however, the recent appreciation that these proteins are signaling proteins in their own right and may even be considered preproproteins controlled by endopeptidases to release a C-terminal bioactive peptide has opened new avenues of research. Also, the recent resolution of the crystal structure of two of the domains of CCN3 have provided new knowledge with implications for the entire CCN family. These resolved structures in combination with structural predictions based upon the AlphaFold artificial intelligence tool provide means to shed new light on CCN functions in context of the notable literature in the field. CCN proteins have emerged as important therapeutic targets in several disease conditions, and clinical trials are currently ongoing. Thus, a review that critically discusses structure - function relationship of CCN proteins, in particular as it relates to interactions with other proteins in the extracellular milieu and on the cell surface, as well as to cell signaling activities of these proteins, is very timely. Suggested mechanism for activation and inhibition of signaling by the CCN protein family (graphics generated with BioRender.com ).

pYtags enable spatiotemporal measurements of receptor tyrosine kinase signaling in living cells

Receptor tyrosine kinases (RTKs) are major signaling hubs in metazoans, playing crucial roles in cell proliferation, migration, and differentiation. However, few tools are available to measure the activity of a specific RTK in individual living cells. Here, we present pYtags, a modular approach for monitoring the activity of a user-defined RTK by live-cell microscopy. pYtags consist of an RTK modified with a tyrosine activation motif that, when phosphorylated, recruits a fluorescently labeled tandem SH2 domain with high specificity. We show that pYtags enable the monitoring of a specific RTK on seconds-to-minutes time scales and across subcellular and multicellular length scales. Using a pYtag biosensor for epidermal growth factor receptor (EGFR), we quantitatively characterize how signaling dynamics vary with the identity and dose of activating ligand. We show that orthogonal pYtags can be used to monitor the dynamics of EGFR and ErbB2 activity in the same cell, revealing distinct phases of activation for each RTK. The specificity and modularity of pYtags open the door to robust biosensors of multiple tyrosine kinases and may enable engineering of synthetic receptors with orthogonal response programs.

High-resolution spatial multi-omics reveals cell-type specific nuclear compartments

The mammalian nucleus is compartmentalized by diverse subnuclear structures. These subnuclear structures, marked by nuclear bodies and histone modifications, are often cell-type specific and affect gene regulation and 3D genome organization1-3. Understanding nuclear organization requires identifying the molecular constituents of subnuclear structures and mapping their associations with specific genomic loci in individual cells, within complex tissues. Here, we introduce two-layer DNA seqFISH+, which allows simultaneous mapping of 100,049 genomic loci, together with nascent transcriptome for 17,856 genes and a diverse set of immunofluorescently labeled subnuclear structures all in single cells in cell lines and adult mouse cerebellum. Using these multi-omics datasets, we showed that repressive chromatin compartments are more variable by cell type than active compartments. We also discovered a single exception to this rule: an RNA polymerase II (RNAPII)-enriched compartment was associated with long, cell-type specific genes (> 200kb), in a manner distinct from nuclear speckles. Further, our analysis revealed that cell-type specific facultative and constitutive heterochromatin compartments marked by H3K27me3 and H4K20me3 are enriched at specific genes and gene clusters, respectively, and shape radial chromosomal positioning and inter-chromosomal interactions in neurons and glial cells. Together, our results provide a single-cell high-resolution multi-omics view of subnuclear compartments, associated genomic loci, and their impacts on gene regulation, directly within complex tissues.

A deep dive into future therapies for microcytic anemias and clinical considerations

INTRODUCTION

Microcytic anemias (MA) have frequent or rare etiologies. New discoveries in understanding and treatment of microcytic anemias need to be reviewed.

AREAS COVERED

Microcytic anemias with a focus on the most frequent causes and on monogenic diseases that are relevant for understanding biocellular mechanisms of MA. All treatments except gene therapy, with a focus on recent advances. PubMed search with references selected by expert opinion.

EXPERT OPINION

As the genetic and cellular backgrounds of dyserythropoiesis will continue to be clarified, collaboration with bioengineering of treatments acting specifically at the protein domain level will continue to provide new therapies in hematology as well as oncology and neurology.

Heterodimerization-dependent secretion of bone morphogenetic proteins in Drosophila

Combinatorial signaling is key to instruct context-dependent cell behaviors. During embryonic development, adult homeostasis, and disease, bone morphogenetic proteins (BMPs) act as dimers to instruct specific cellular responses. BMP ligands can form both homodimers or heterodimers; however, obtaining direct evidence of the endogenous localization and function of each form has proven challenging. Here, we make use of precise genome editing and direct protein manipulation via protein binders to dissect the existence and functional relevance of BMP homodimers and heterodimers in the Drosophila wing imaginal disc. This approach identified in situ the existence of Dpp (BMP2/4)/Gbb (BMP5/6/7/8) heterodimers. We found that Gbb is secreted in a Dpp-dependent manner in the wing imaginal disc. Dpp and Gbb form a gradient of heterodimers, whereas neither Dpp nor Gbb homodimers are evident under endogenous physiological conditions. We find that the formation of heterodimers is critical for obtaining optimal signaling and long-range BMP distribution.

The beauty of seeing the real thing: BMP heterodimer detection in vivo reveals dimer composition and activity

BMP family ligands can direct cells to divide, differentiate, or die, depending on cell context and specific hetero- or homodimer combinations. In this issue of Developmental Cell, Bauer et al. detect endogenous Drosophila ligand dimers in situ and show that BMP dimer composition affects both signaling range and activity.

The time integral of BMP signaling determines fate in a stem cell model for early human development

How paracrine signals are interpreted to yield multiple cell fate decisions in a dynamic context during human development in vivo and in vitro remains poorly understood. Here we report an automated tracking method to follow signaling histories linked to cell fate in large numbers of human pluripotent stem cells (hPSCs). Using an unbiased statistical approach, we discovered that measured BMP signaling history correlates strongly with fate in individual cells. We found that BMP response in hPSCs varies more strongly in the duration of signaling than the level. However, we discovered that both the level and duration of signaling activity control cell fate choices only by changing the time integral of signaling and that duration and level are therefore interchangeable in this context. In a stem cell model for patterning of the human embryo, we showed that signaling histories predict the fate pattern and that the integral model correctly predicts changes in cell fate domains when signaling is perturbed. Using an RNA-seq screen we then found that mechanistically, BMP signaling is integrated by SOX2.

State- and stimulus-specific dynamics of SMAD signaling determine fate decisions in individual cells

SMAD-mediated signaling regulates apoptosis, cell cycle arrest, and epithelial-to-mesenchymal transition to safeguard tissue homeostasis. However, it remains elusive how the relatively simple pathway can determine such a broad range of cell fate decisions and how it differentiates between varying ligands. Here, we systematically investigate how SMAD-mediated responses are modulated by various ligands of the transforming growth factor β (TGFβ) family and compare these ligand responses in quiescent and proliferating MCF10A cells. We find that the nature of the phenotypic response is mainly determined by the proliferation status, with migration and cell cycle arrest being dominant in proliferating cells for all tested TGFβ family ligands, whereas cell death is the major outcome in quiescent cells. In both quiescent and proliferating cells, the identity of the ligand modulates the strength of the phenotypic response proportional to the dynamics of induced SMAD nuclear-to-cytoplasmic translocation and, as a consequence, the corresponding gene expression changes. Interestingly, the proliferation state of a cell has little impact on the set of genes induced by SMAD signaling; instead, it modulates the relative cellular sensitivity to TGFβ superfamily members. Taken together, diversity of SMAD-mediated responses is mediated by differing cellular states, which determine ligand sensitivity and phenotypic effects, while the pathway itself merely serves as a quantitative relay from the cell membrane to the nucleus.

Molecular Computation for Molecular Classification

DNA as an informational polymer has, for the past 30 years, progressively become an essential molecule to rationally build chemical reaction networks endowed with powerful signal-processing capabilities. Whether influenced by the silicon world or inspired by natural computation, molecular programming has gained attention for diagnosis applications. Of particular interest for this review, molecular classifiers have shown promising results for disease pattern recognition and sample classification. Because both input integration and computation are performed in a single tube, at the molecular level, this low-cost approach may come as a complementary tool to molecular profiling strategies, where all biomarkers are quantified independently using high-tech instrumentation. After introducing the elementary components of molecular classifiers, some of their experimental implementations are discussed either using digital Boolean logic or analog neural network architectures.

Improving cooperativity of transcription activators by oligomerization domains in mammalian cells

Cooperative activation is critical for the applications of synthetic biology in mammalian cells. In this study, we have developed cooperative transcription factor by fusing oligomerization domain in mammalian cells. Firstly, we demonstrated that two oligomerized domains (CI434 and CI) successfully improved transcription factor cooperativity in bacterial cells but failed to increase cooperativity in mammalian cells, possibly because the additional mammalian activation domain disrupted their oligomerization capability. Therefore, we chose a different type of oligomerized domain (CarHC), whose ability to oligomerize is not dependent on its C-terminal domains, to fuse with a transcription factor (RpaR) and activation domain (VTR3), forming a potential cooperative transcription activator RpaR-CarH-VTR3 for mammalian regulatory systems. Compared with RpaR-VTR3, the cooperativity of RpaR-CarH-VTR3 was significantly improved with higher Hill coefficient and a narrower input range in the inducible switch system in mammalian cells. Moreover, a mathematical model based on statistical mechanics model was developed and the simulation results supported the hypothesis that the tetramer of the CarH domain in mammalian cells was the reason for the cooperative capacity of RpaR-CarH-VTR3.

How cells sense and integrate information from different sources

Cell signaling is a fundamental cellular process that enables cells to sense and respond to information in their surroundings. At the molecular level, signaling is primarily carried out by transmembrane protein receptors that can initiate complex downstream signal transduction cascades to alter cellular behavior. In the human body, different cells can be exposed to a wide variety of environmental conditions, and cells express diverse classes of receptors capable of sensing and integrating different signals. Furthermore, different receptors and signaling pathways can crosstalk with each other to calibrate the cellular response. Crosstalk occurs through multiple mechanisms at different levels of signaling pathways. In this review, we discuss how cells sense and integrate different chemical, mechanical, and spatial signals as well as the mechanisms of crosstalk between pathways. To illustrate these concepts, we use a few well-studied signaling pathways, including receptor tyrosine kinases and integrin receptors. Finally, we discuss the implications of dysregulated cellular sensing on driving diseases such as cancer. This article is categorized under: Cancer > Molecular and Cellular Physiology Metabolic Diseases > Molecular and Cellular Physiology.

FKBP12 is a major regulator of ALK2 activity in multiple myeloma cells

BACKGROUND

The immunophilin FKBP12 binds to TGF-β family type I receptors, including the BMP type I receptor ALK2. FKBP12 keeps the type I receptor in an inactive state and controls signaling activity. Removal of FKBP12 with drugs such as the FKBP-ligand FK506 enhances BMP activity in various cell types. In multiple myeloma cells, activation of SMAD1/5/8 leads to apoptosis. We hypothesized that removing FKBP12 from ALK2 in myeloma cells would potentiate BMP-induced ALK2-SMAD1/5/8 activity and in consequence cell death.

METHODS

Multiple myeloma cell lines were treated with FK506, or other FKBP-binding compounds, combined with different BMPs before analyzing SMAD1/5/8 activity and cell viability. SMAD1/5/8 activity was also investigated using a reporter cell line, INA-6 BRE-luc. To characterize the functional signaling receptor complex, we genetically manipulated receptor expression by siRNA, shRNA and CRISPR/Cas9 technology.

RESULTS

FK506 potentiated BMP-induced SMAD1/5/8 activation and apoptosis in multiple myeloma cell lines. By using FKBP-binding compounds with different affinity profiles, and siRNA targeting FKBP12, we show that the FK506 effect is mediated by binding to FKBP12. Ligands that typically signal via ALK3 in myeloma cells, BMP2, BMP4, and BMP10, did not induce apoptosis in cells lacking ALK3. Notably, BMP10 competed with BMP6 and BMP9 and antagonized their activity via ALK2. However, upon addition of FK506, we saw a surprising shift in specificity, as the ALK3 ligands gained the ability to signal via ALK2 and induce apoptosis. This indicates that the receptor complex can switch from an inactive non-signaling complex (NSC) to an active one by adding FK506. This gain of activity was also seen in other cell types, indicating that the observed effects have broader relevance. BMP2, BMP4 and BMP10 depended on BMPR2 as type II receptor to signal, which contrasts with BMP6 and BMP9, that activate ALK2 more potently when BMPR2 is knocked down.

CONCLUSIONS

In summary, our data suggest that FKBP12 is a major regulator of ALK2 activity in multiple myeloma cells, partly by switching an NSC into an active signaling complex. FKBP12 targeting compounds devoid of immunosuppressing activity could have potential in novel treatment strategies aiming at reducing multiple myeloma tumor load. Video Abstract.

The Topological Characteristics of Biological Ratio-Sensing Networks

Ratio sensing is a fundamental biological function observed in signal transduction and decision making. In the synthetic biology context, ratio sensing presents one of the elementary functions for cellular multi-signal computation. To investigate the mechanism of the ratio-sensing behavior, we explored the topological characteristics of biological ratio-sensing networks. With exhaustive enumeration of three-node enzymatic and transcriptional regulatory networks, we found that robust ratio sensing was highly dependent on network structure rather than network complexity. Specifically, a set of seven minimal core topological structures and four motifs were deduced to be capable of robust ratio sensing. Further investigations on the evolutionary space of robust ratio-sensing networks revealed highly clustered domains surrounding the core motifs which suggested their evolutionary plausibility. Our study revealed the network topological design principles of ratio-sensing behavior and provided a design scheme for constructing regulatory circuits with ratio-sensing behavior in synthetic biology.

Molecular versatility during pluripotency progression

A challenge during development is to ensure lineage segregation while preserving plasticity. Using pluripotency progression as a paradigm, we review how developmental transitions are coordinated by redeployments, rather than global resettings, of cellular components. We highlight how changes in response to extrinsic cues (FGF, WNT, Activin/Nodal, Netrin-1), context- and stoichiometry-dependent action of transcription factors (Oct4, Nanog) and reconfigurations of epigenetic regulators (enhancers, promoters, TrxG, PRC) may confer robustness to naïve to primed pluripotency transition. We propose the notion of Molecular Versatility to regroup mechanisms by which molecules are repurposed to exert different, sometimes opposite, functions in close stem cell configurations.