A Unifying Theory of Branching Morphogenesis.

The morphogenesis of branched organs remains a subject of abiding interest. Although much is known about the underlying signaling pathways, it remains unclear how macroscopic features of branched organs, including their size, network topology, and spatial patterning, are encoded. Here, we show that, in mouse mammary gland, kidney, and human prostate, these features can be explained quantitatively within a single unifying framework of branching and annihilating random walks. Based on quantitative analyses of large-scale organ reconstructions and proliferation kinetics measurements, we propose that morphogenesis follows from the proliferative activity of equipotent tips that stochastically branch and randomly explore their environment but compete neutrally for space, becoming proliferatively inactive when in proximity with neighboring ducts. These results show that complex branched epithelial structures develop as a self-organized process, reliant upon a strikingly simple but generic rule, without recourse to a rigid and deterministic sequence of genetically programmed events.

Mechanisms that clear mutations drive field cancerization in mammary tissue.

Oncogenic mutations are abundant in the tissues of healthy individuals, but rarely form tumours1-3. Yet, the underlying protection mechanisms are largely unknown. To resolve these mechanisms in mouse mammary tissue, we use lineage tracing to map the fate of wild-type and Brca1-/-;Trp53-/- cells, and find that both follow a similar pattern of loss and spread within ducts. Clonal analysis reveals that ducts consist of small repetitive units of self-renewing cells that give rise to short-lived descendants. This offers a first layer of protection as any descendants, including oncogenic mutant cells, are constantly lost, thereby limiting the spread of mutations to a single stem cell-descendant unit. Local tissue remodelling during consecutive oestrous cycles leads to the cooperative and stochastic loss and replacement of self-renewing cells. This process provides a second layer of protection, leading to the elimination of most mutant clones while enabling the minority that by chance survive to expand beyond the stem cell-descendant unit. This leads to fields of mutant cells spanning large parts of the epithelial network, predisposing it for transformation. Eventually, clone expansion becomes restrained by the geometry of the ducts, providing a third layer of protection. Together, these mechanisms act to eliminate most cells that acquire somatic mutations at the expense of driving the accelerated expansion of a minority of cells, which can colonize large areas, leading to field cancerization.

Metabolic adaptation towards glycolysis supports resistance to neoadjuvant chemotherapy in early triple negative breast cancers.

BACKGROUND: Neoadjuvant chemotherapy (NAC) is the standard of care for patients with early-stage triple negative breast cancers (TNBC). However, more than half of TNBC patients do not achieve a pathological complete response (pCR) after NAC, and residual cancer burden (RCB) is associated with dismal long-term prognosis. Understanding the mechanisms underlying differential treatment outcomes is therefore critical to limit RCB and improve NAC efficiency. METHODS: Human TNBC cell lines and patient-derived organoids were used in combination with real-time metabolic assays to evaluate the effect of NAC (paclitaxel and epirubicin) on tumor cell metabolism, in particular glycolysis. Diagnostic biopsies (pre-NAC) from patients with early TNBC were analyzed by bulk RNA-sequencing to evaluate the predictive value of a glycolysis-related gene signature. RESULTS: Paclitaxel induced a consistent metabolic switch to glycolysis, correlated with a reduced mitochondrial oxidative metabolism, in TNBC cells. In pre-NAC diagnostic biopsies from TNBC patients, glycolysis was found to be upregulated in non-responders. Furthermore, glycolysis inhibition greatly improved response to NAC in TNBC organoid models. CONCLUSIONS: Our study pinpoints a metabolic adaptation to glycolysis as a mechanism driving resistance to NAC in TNBC. Our data pave the way for the use of glycolysis-related genes as predictive biomarkers for NAC response, as well as the development of inhibitors to overcome this glycolysis-driven resistance to NAC in human TNBC patients.

Single-cell analysis of regions of interest (SCARI) using a photosensitive tag.

The functional activity and differentiation potential of cells are determined by their interactions with surrounding cells. Approaches that allow unbiased characterization of cell states while at the same time providing spatial information are of major value to assess this environmental influence. However, most current techniques are hampered by a tradeoff between spatial resolution and cell profiling depth. Here, we develop a photocage-based technology that allows isolation and in-depth analysis of live cells from regions of interest in complex ex vivo systems, including primary human tissues. The use of a highly sensitive 4-nitrophenyl(benzofuran) cage coupled to a set of nanobodies allows high-resolution photo-uncaging of different cell types in areas of interest. Single-cell RNA-sequencing of spatially defined CD8+ T cells is used to exemplify the feasibility of identifying location-dependent cell states. The technology described here provides a valuable tool for the analysis of spatially defined cells in diverse biological systems, including clinical samples.

Identity and dynamics of mammary stem cells during branching morphogenesis.

During puberty, the mouse mammary gland develops into a highly branched epithelial network. Owing to the absence of exclusive stem cell markers, the location, multiplicity, dynamics and fate of mammary stem cells (MaSCs), which drive branching morphogenesis, are unknown. Here we show that morphogenesis is driven by proliferative terminal end buds that terminate or bifurcate with near equal probability, in a stochastic and time-invariant manner, leading to a heterogeneous epithelial network. We show that the majority of terminal end bud cells function as highly proliferative, lineage-committed MaSCs that are heterogeneous in their expression profile and short-term contribution to ductal extension. Yet, through cell rearrangements during terminal end bud bifurcation, each MaSC is able to contribute actively to long-term growth. Our study shows that the behaviour of MaSCs is not directly linked to a single expression profile. Instead, morphogenesis relies upon lineage-restricted heterogeneous MaSC populations that function as single equipotent pools in the long term.

Stem cell lineage survival as a noisy competition for niche access.

Understanding to what extent stem cell potential is a cell-intrinsic property or an emergent behavior coming from global tissue dynamics and geometry is a key outstanding question of systems and stem cell biology. Here, we propose a theory of stem cell dynamics as a stochastic competition for access to a spatially localized niche, giving rise to a stochastic conveyor-belt model. Cell divisions produce a steady cellular stream which advects cells away from the niche, while random rearrangements enable cells away from the niche to be favorably repositioned. Importantly, even when assuming that all cells in a tissue are molecularly equivalent, we predict a common ("universal") functional dependence of the long-term clonal survival probability on distance from the niche, as well as the emergence of a well-defined number of functional stem cells, dependent only on the rate of random movements vs. mitosis-driven advection. We test the predictions of this theory on datasets of pubertal mammary gland tips and embryonic kidney tips, as well as homeostatic intestinal crypts. Importantly, we find good agreement for the predicted functional dependency of the competition as a function of position, and thus functional stem cell number in each organ. This argues for a key role of positional fluctuations in dictating stem cell number and dynamics, and we discuss the applicability of this theory to other settings.

An Intravital Microscopy Toolbox to Study Mammary Gland Dynamics from Cellular Level to Organ Scale.

The architecture of the mouse mammary gland is highly dynamic and constantly remodeled during pubertal development and estrous cycle-driven sprouting and regression of alveolar side branches. During each of these developmental stages, turnover is driven by distinct subsets of mammary epithelial cells. Extensive previous research has shed light on the unique morphological and cell biological characteristics of each stage. However, technological shortcomings failed to capture the dynamics and single-cell contributions to mammary remodeling. Here, we developed in vivo imaging strategies to follow the same mammary ducts over time and quantify the dynamics of mammary gland growth and remodeling from single-cell level to organ scale. Using a combination of intravital microscopy and genetic reporter systems we show how proliferative heterogeneity drives ductal morphogenesis during different developmental stages. To visualize pubertal growth at the cellular level, we performed long-term time-lapse imaging of extending terminal end buds through a mammary imaging window. We show that single-cells within the terminal end buds are extremely motile and continuously exchange position whilst the duct is elongating. To visualize short-term remodeling in the adult mammary gland at the single cell level, we performed multi-day intravital imaging in photoconvertible Kikume Green-Red mice and fluorescent ubiquitination-based cell cycle indicator mice. We demonstrate that the contribution of single-cells to estrous-driven remodeling is highly variable between cells in the same micro-environment. To assess the effects of this dynamic proliferative contribution on the long-term stability of tissue architecture, we developed a repeated skin flap method to assess mammary gland morphology by intravital microscopy over extended time spans for up to six months. Interestingly, in contrast to the short-term dynamic remodeling, the long-term morphology of the mammary gland remains remarkably stable. Together, our tool box of imaging strategies allows to identify and map transient and continuing dynamics of single cells to the architecture of the mammary gland.

E2F7 and E2F8 promote angiogenesis through transcriptional activation of VEGFA in cooperation with HIF1.

The E2F family of transcription factors plays an important role in controlling cell-cycle progression. While this is their best-known function, we report here novel functions for the newest members of the E2F family, E2F7 and E2F8 (E2F7/8). We show that simultaneous deletion of E2F7/8 in zebrafish and mice leads to severe vascular defects during embryonic development. Using a panel of transgenic zebrafish with fluorescent-labelled blood vessels, we demonstrate that E2F7/8 are essential for proper formation of blood vessels. Despite their classification as transcriptional repressors, we provide evidence for a molecular mechanism through which E2F7/8 activate the transcription of the vascular endothelial growth factor A (VEGFA), a key factor in guiding angiogenesis. We show that E2F7/8 directly bind and stimulate the VEGFA promoter independent of canonical E2F binding elements. Instead, E2F7/8 form a transcriptional complex with the hypoxia inducible factor 1 (HIF1) to stimulate VEGFA promoter activity. These results uncover an unexpected link between E2F7/8 and the HIF1-VEGFA pathway providing a molecular mechanism by which E2F7/8 control angiogenesis.

A Decision Tree to Guide Human and Mouse Mammary Organoid Model Selection.

Over the past 50 years, researchers from the mammary gland field have launched a collection of distinctive 3D cell culture systems to study multiple aspects of mammary gland physiology and disease. As our knowledge about the mammary gland evolves, more sophisticated 3D cell culture systems are required to answer more and more complex questions. Nowadays, morphologically complex mammary organoids can be generated in distinct 3D settings, along with reproduction of multiple aspects of the gland microenvironment. Yet, each 3D culture protocol comes with its advantages and limitations, where some culture systems are best suited to study stemness potential, whereas others are tailored towards the study of mammary gland morphogenesis. Therefore, prior to starting a 3D mammary culture experiment, it is important to consider and select the ideal culture model to address the biological question of interest. The number and technical requirements of novel 3D cell culture methods vastly increased over the past decades, making it currently challenging and time consuming to identify the best experimental testing. In this chapter, we provide a summary of the most promising murine and human 3D organoid models that are currently used in mammary gland biology research. For each model, we will provide a brief description of the protocol and an overview of the expected morphological outcome, the advantages of the model, and the potential pitfalls, to guide the reader to the best model of choice for specific applications.

Smad1/5 is acetylated in the dorsal aortae of the mouse embryo before the onset of blood flow, driving early arterial gene expression.

AIMS: During embryonic development, arteriovenous (AV) differentiation ensures proper blood vessel formation and maturation. Defects in arterial or venous identity cause inappropriate fusion of vessels, resulting in atypical shunts, so-called arteriovenous malformations (AVM). Currently, the mechanism behind AVM formation remains unclear and treatment options are fairly limited. Mammalian AV differentiation is initiated before the onset of blood flow in the embryo; however, this pre-flow mechanism is poorly understood. Here, we aimed to unravel the role of Smad1/5 signalling in pre-flow arterial identity, and in the process uncovered an unexpected control mechanism of Smad1/5 signalling. METHODS AND RESULTS: We establish that despite Notch1 being expressed in the pre-flow mouse embryo, it is not activated, nor is it necessary for the expression of the earliest arterial genes in the dorsal aortae (i.e., Hey1 and Gja4). Furthermore, interrupting blood flow by using the Ncx1 KO model completely prevents the activation of Notch1 signalling, suggesting a strong role of shear stress in maintaining arterial identity. We demonstrate that early expression of Hey1 and Gja4 requires SMAD1/5 signalling. Using embryo cultures, we show that Smad1/5 signalling is activated through the Alk1/Alk5/TGFßR2 receptor complex, with TGFß1 as a necessary ligand. Furthermore, our findings demonstrate that early arterial gene expression requires the acetylation of Smad1/5 proteins, rendering them more sensitive to TGFß1 stimulation. Blocking acetyl-CoA production prevents pre-flow arterial expression of Hey1 and Gja4, while stabilizing acetylation rescues their expression. CONCLUSIONS: Our findings highlight the importance of the acetyl-CoA production in the cell and provide a novel control mechanism of Smad1/5 signalling involving protein acetylation. As disturbed canonical Smad1/5 signalling is involved in several vascular conditions, our results offer new insights in treatment options for circumventing canonical Smad1/5 signalling.

Volumetric analysis of the terminal ductal lobular unit architecture and cell phenotypes in the human breast.

The major lactiferous ducts of the human breast branch out and end at terminal ductal lobular units (TDLUs). Despite their functional and clinical importance, the three-dimensional (3D) architecture of TDLUs has remained undetermined. Our quantitative and volumetric imaging of healthy human breast tissue demonstrates that highly branched TDLUs, which exhibit increased proliferation, are uncommon in the resting tissue regardless of donor age, parity, or hormonal contraception. Overall, TDLUs have a consistent shape and branch parameters, and they contain a main subtree that dominates in bifurcation events and exhibits a more duct-like keratin expression pattern. Simulation of TDLU branching morphogenesis in three dimensions suggests that evolutionarily conserved mechanisms regulate mammary gland branching in humans and mice despite their anatomical differences. In all, our data provide structural insight into 3D anatomy and branching of the human breast and exemplify the power of volumetric imaging in gaining a deeper understanding of breast biology.

An autophagy program that promotes T cell egress from the lymph node controls responses to immune checkpoint blockade.

Lymphatic endothelial cells (LECs) of the lymph node (LN) parenchyma orchestrate leukocyte trafficking and peripheral T cell dynamics. T cell responses to immunotherapy largely rely on peripheral T cell recruitment in tumors. Yet, a systematic and molecular understanding of how LECs within the LNs control T cell dynamics under steady-state and tumor-bearing conditions is lacking. Intravital imaging combined with immune phenotyping shows that LEC-specific deletion of the essential autophagy gene Atg5 alters intranodal positioning of lymphocytes and accrues their persistence in the LNs by increasing the availability of the main egress signal sphingosine-1-phosphate. Single-cell RNA sequencing of tumor-draining LNs shows that loss of ATG5 remodels niche-specific LEC phenotypes involved in molecular pathways regulating lymphocyte trafficking and LEC-T cell interactions. Functionally, loss of LEC autophagy prevents recruitment of tumor-infiltrating T and natural killer cells and abrogates response to immunotherapy. Thus, an LEC-autophagy program boosts immune-checkpoint responses by guiding systemic T cell dynamics.

Rapid autopsies to enhance metastatic research: the UPTIDER post-mortem tissue donation program.

Research on metastatic cancer has been hampered by limited sample availability. Here we present the breast cancer post-mortem tissue donation program UPTIDER and show how it enabled sampling of a median of 31 (range: 5-90) metastases and 5-8 liquids per patient from its first 20 patients. In a dedicated experiment, we show the mild impact of increasing time after death on RNA quality, transcriptional profiles and immunohistochemical staining in tumor tissue samples. We show that this impact can be counteracted by organ cooling. We successfully generated ex vivo models from tissue and liquid biopsies from distinct histological subtypes of breast cancer. We anticipate these and future findings of UPTIDER to elucidate mechanisms of disease progression and treatment resistance and to provide tools for the exploration of precision medicine strategies in the metastatic setting.

A Guide Toward Multi-scale and Quantitative Branching Analysis in the Mammary Gland.

The mammary gland consists of a bilayered epithelial structure with an extensively branched morphology. The majority of this epithelial tree is laid down during puberty, during which actively proliferating terminal end buds repeatedly elongate and bifurcate to form the basic structure of the ductal tree. Mammary ducts consist of a basal and luminal cell layer with a multitude of identified sub-lineages within both layers. The understanding of how these different cell lineages are cooperatively driving branching morphogenesis is a problem of crossing multiple scales, as this requires information on the macroscopic branched structure of the gland, as well as data on single-cell dynamics driving the morphogenic program. Here we describe a method to combine genetic lineage tracing with whole-gland branching analysis. Quantitative data on the global organ structure can be used to derive a model for mammary gland branching morphogenesis and provide a backbone on which the dynamics of individual cell lineages can be simulated and compared to lineage-tracing approaches. Eventually, these quantitative models and experiments allow to understand the couplings between the macroscopic shape of the mammary gland and the underlying single-cell dynamics driving branching morphogenesis.

Mouse intraductal modeling of primary ductal carcinoma in situ.

Mouse intraductal modeling enables efficient in vivo propagation of pre-invasive breast cancer lesions and provides a suitable micro-environment for creating patient-derived tumor xenograft models of estrogen-receptor-positive breast cancer. Here, we present a protocol for mouse intraductal modeling of primary ductal carcinoma in situ (DCIS). We describe steps for processing primary DCIS tissues and performing intraductal injections. We then detail procedures for processing intraductal lesions for 3D whole-mount imaging or serial transplantation using magnetic bead sorting. For complete details on the use and execution of this protocol, please refer to Hutten et al. (2023).1.

Multidimensional Imaging of Breast Cancer.

Breast cancer is a pathological condition characterized by high morphological and molecular heterogeneity. Not only the breast cancer cells, but also their tumor micro-environment consists of a multitude of cell types and states, which continuously evolve throughout progression of the disease. To understand breast cancer evolution within this complex environment, in situ analysis of breast cancer and their co-evolving cells and structures in space and time are essential. In this review, recent technical advances in three-dimensional (3D) and intravital imaging of breast cancer are discussed. Moreover, we highlight the resulting new knowledge on breast cancer biology obtained through these innovative imaging technologies. Finally, we discuss how multidimensional imaging technologies can be integrated with molecular profiling to understand the full complexity of breast cancer and the tumor micro-environment during tumor progression and treatment response.

Annexin A1 is a polarity cue that directs mitotic spindle orientation during mammalian epithelial morphogenesis.

Oriented cell divisions are critical for the formation and maintenance of structured epithelia. Proper mitotic spindle orientation relies on polarised anchoring of force generators to the cell cortex by the evolutionarily conserved protein complex formed by the Gαi subunit of heterotrimeric G proteins, the Leucine-Glycine-Asparagine repeat protein (LGN) and the nuclear mitotic apparatus protein. However, the polarity cues that control cortical patterning of this ternary complex remain largely unknown in mammalian epithelia. Here we identify the membrane-associated protein Annexin A1 (ANXA1) as an interactor of LGN in mammary epithelial cells. Annexin A1 acts independently of Gαi to instruct the accumulation of LGN and nuclear mitotic apparatus protein at the lateral cortex to ensure cortical anchoring of Dynein-Dynactin and astral microtubules and thereby planar alignment of the mitotic spindle. Loss of Annexin A1 randomises mitotic spindle orientation, which in turn disrupts epithelial architecture and luminogenesis in three-dimensional cultures of primary mammary epithelial cells. Our findings establish Annexin A1 as an upstream cortical cue that regulates LGN to direct planar cell divisions during mammalian epithelial morphogenesis.

CD26-negative and CD26-positive tissue-resident fibroblasts contribute to functionally distinct CAF subpopulations in breast cancer.

Cancer-associated fibroblasts (CAFs) are abundantly present in the microenvironment of virtually all tumors and strongly impact tumor progression. Despite increasing insight into their function and heterogeneity, little is known regarding the origin of CAFs. Understanding the origin of CAF heterogeneity is needed to develop successful CAF-based targeted therapies. Through various transplantation studies in mice, we show that CAFs in both invasive lobular breast cancer and triple-negative breast cancer originate from mammary tissue-resident normal fibroblasts (NFs). Single-cell transcriptomics, in vivo and in vitro studies reveal the transition of CD26+ and CD26- NF populations into inflammatory CAFs (iCAFs) and myofibroblastic CAFs (myCAFs), respectively. Functional co-culture experiments show that CD26+ NFs transition into pro-tumorigenic iCAFs which recruit myeloid cells in a CXCL12-dependent manner and enhance tumor cell invasion via matrix-metalloproteinase (MMP) activity. Together, our data suggest that CD26+ and CD26- NFs transform into distinct CAF subpopulations in mouse models of breast cancer.

Breast tumors interfere with endothelial TRAIL at the premetastatic niche to promote cancer cell seeding.

Endothelial cells (ECs) grant access of disseminated cancer cells to distant organs. However, the molecular players regulating the activation of quiescent ECs at the premetastatic niche (PMN) remain elusive. Here, we find that ECs at the PMN coexpress tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) and its cognate death receptor 5 (DR5). Unexpectedly, endothelial TRAIL interacts intracellularly with DR5 to prevent its signaling and preserve a quiescent vascular phenotype. In absence of endothelial TRAIL, DR5 activation induces EC death and nuclear factor κB/p38-dependent EC stickiness, compromising vascular integrity and promoting myeloid cell infiltration, breast cancer cell adhesion, and metastasis. Consistently, both down-regulation of endothelial TRAIL at the PMN by proangiogenic tumor-secreted factors and the presence of the endogenous TRAIL inhibitors decoy receptor 1 (DcR1) and DcR2 favor metastasis. This study discloses an intracrine mechanism whereby TRAIL blocks DR5 signaling in quiescent endothelia, acting as gatekeeper of the vascular barrier that is corrupted by the tumor during cancer cell dissemination.

A living biobank of patient-derived ductal carcinoma in situ mouse-intraductal xenografts identifies risk factors for invasive progression.

Ductal carcinoma in situ (DCIS) is a non-obligate precursor of invasive breast cancer (IBC). Due to a lack of biomarkers able to distinguish high- from low-risk cases, DCIS is treated similar to early IBC even though the minority of untreated cases eventually become invasive. Here, we characterized 115 patient-derived mouse-intraductal (MIND) DCIS models reflecting the full spectrum of DCIS observed in patients. Utilizing the possibility to follow the natural progression of DCIS combined with omics and imaging data, we reveal multiple prognostic factors for high-risk DCIS including high grade, HER2 amplification, expansive 3D growth, and high burden of copy number aberrations. In addition, sequential transplantation of xenografts showed minimal phenotypic and genotypic changes over time, indicating that invasive behavior is an intrinsic phenotype of DCIS and supporting a multiclonal evolution model. Moreover, this study provides a collection of 19 distributable DCIS-MIND models spanning all molecular subtypes.