Mapping hormone-regulated cell-cell interaction networks in the human breast at single-cell resolution

Single-nucleus chromatin accessibility and transcriptomic map of breast tissues of women of diverse genetic ancestry

Single-nucleus analysis allows robust cell-type classification and helps to establish relationships between chromatin accessibility and cell-type-specific gene expression. Here, using samples from 92 women of several genetic ancestries, we developed a comprehensive chromatin accessibility and gene expression atlas of the breast tissue. Integrated analysis revealed ten distinct cell types, including three major epithelial subtypes (luminal hormone sensing, luminal adaptive secretory precursor (LASP) and basal-myoepithelial), two endothelial and adipocyte subtypes, fibroblasts, T cells, and macrophages. In addition to the known cell identity genes FOXA1 (luminal hormone sensing), EHF and ELF5 (LASP), TP63 and KRT14 (basal-myoepithelial), epithelial subtypes displayed several uncharacterized markers and inferred gene regulatory networks. By integrating breast epithelial cell gene expression signatures with spatial transcriptomics, we identified gene expression and signaling differences between lobular and ductal epithelial cells and age-associated changes in signaling networks. LASP cells and fibroblasts showed genetic ancestry-dependent variability. An estrogen receptor-positive subpopulation of LASP cells with alveolar progenitor cell state was enriched in women of Indigenous American ancestry. Fibroblasts from breast tissues of women of African and European ancestry clustered differently, with accompanying gene expression differences. Collectively, these data provide a vital resource for further exploring genetic ancestry-dependent variability in healthy breast biology.

Breast organoid suspension cultures maintain long-term estrogen receptor expression and responsiveness

Organoid cultures offer a powerful technology to investigate many different aspects of development, physiology, and pathology of diverse tissues. Unlike standard tissue culture of primary breast epithelial cells, breast organoids preserve the epithelial lineages and architecture of the normal tissue. However, existing organoid culture methods are tedious, difficult to scale, and do not robustly retain estrogen receptor (ER) expression and responsiveness in long-term culture. Here, we describe a modified culture method to generate and maintain organoids as suspension cultures in reconstituted basement membrane (™Matrigel). This method improves organoid growth and uniformity compared to the conventional Matrigel dome embedding method, while maintaining the fidelity of the three major epithelial lineages. Using this adopted method, we are able to culture and passage purified hormone sensing (HS) cells that retain ER responsiveness upon estrogen stimulation in long-term culture. This culture system presents a valuable platform to study the events involved in initiation and evolution of ER-positive breast cancer.

Coordinated immune dysregulation in juvenile dermatomyositis revealed by single-cell genomics

Juvenile dermatomyositis (JDM) is one of several childhood-onset autoimmune disorders characterized by a type I IFN response and autoantibodies. Treatment options are limited due to an incomplete understanding of how the disease emerges from dysregulated cell states across the immune system. We therefore investigated the blood of patients with JDM at different stages of disease activity using single-cell transcriptomics paired with surface protein expression. By immunophenotyping peripheral blood mononuclear cells, we observed skewing of the B cell compartment toward an immature naive state as a hallmark of JDM at diagnosis. Furthermore, we find that these changes in B cells are paralleled by T cell signatures suggestive of Th2-mediated inflammation that persist despite disease quiescence. We applied network analysis to reveal that hyperactivation of the type I IFN response in all immune populations is coordinated with previously masked cell states including dysfunctional protein processing in CD4+ T cells and regulation of cell death programming in NK cells, CD8+ T cells, and γδ T cells. Together, these findings unveil the coordinated immune dysregulation underpinning JDM and provide insight into strategies for restoring balance in immune function.

A single-cell atlas enables mapping of homeostatic cellular shifts in the adult human breast

Here we use single-cell RNA sequencing to compile a human breast cell atlas assembled from 55 donors that had undergone reduction mammoplasties or risk reduction mastectomies. From more than 800,000 cells we identified 41 cell subclusters across the epithelial, immune and stromal compartments. The contribution of these different clusters varied according to the natural history of the tissue. Age, parity and germline mutations, known to modulate the risk of developing breast cancer, affected the homeostatic cellular state of the breast in different ways. We found that immune cells from BRCA1 or BRCA2 carriers had a distinct gene expression signature indicative of potential immune exhaustion, which was validated by immunohistochemistry. This suggests that immune-escape mechanisms could manifest in non-cancerous tissues very early during tumor initiation. This atlas is a rich resource that can be used to inform novel approaches for early detection and prevention of breast cancer.

Mammary Microvessels are Sensitive to Menstrual Cycle Sex Hormones

The mammary gland is a highly vascularized organ influenced by sex hormones including estrogen (E2) and progesterone (P4). Beyond whole-organism studies in rodents or cell monocultures, hormonal effects on the breast microvasculature remain largely understudied. Recent methods to generate 3D microvessels on-chip have enabled direct observation of complex vascular processes; however, these models often use non-tissue-specific cell types, such as human umbilical vein endothelial cells (HUVECs) and fibroblasts from various sources. Here, novel mammary-specific microvessels are generated by coculturing primary breast endothelial cells and fibroblasts under optimized culture conditions. These microvessels are mechanosensitive (to interstitial flow) and require endothelial-stromal interactions to develop fully perfusable vessels. These mammary-specific microvessels are also responsive to exogenous stimulation by sex hormones. When treated with combined E2 and P4, corresponding to the four phases of the menstrual cycle (period, follicular, ovular, and luteal), vascular remodeling and barrier function are altered in a phase-dependent manner. The presence of high E2 (ovulation) promotes vascular growth and remodeling, corresponding to high depletion of proangiogenic factors, whereas high P4 concentrations (luteal) promote vascular regression. The effects of combined E2 and P4 hormones are not only dose-dependent but also tissue-specific, as are shown by similarly treating non-tissue-specific HUVEC microvessels.

Coordinated immune dysregulation in Juvenile Dermatomyositis revealed by single-cell genomics

Juvenile Dermatomyositis (JDM) is one of several childhood-onset autoimmune disorders characterized by a type I interferon response and autoantibodies. Treatment options are limited due to incomplete understanding of how the disease emerges from dysregulated cell states across the immune system. We therefore investigated the blood of JDM patients at different stages of disease activity using single-cell transcriptomics paired with surface protein expression. By immunophenotyping peripheral blood mononuclear cells, we observed skewing of the B cell compartment towards an immature naive state as a hallmark of JDM. Furthermore, we find that these changes in B cells are paralleled by signatures of Th2-mediated inflammation. Additionally, our work identified SIGLEC-1 expression in monocytes as a composite measure of heterogeneous type I interferon activity in disease. We applied network analysis to reveal that hyperactivation of the type I interferon response in all immune populations is coordinated with dysfunctional protein processing and regulation of cell death programming. This analysis separated the ubiquitously expressed type I interferon response into a central hub and revealed previously masked cell states. Together, these findings reveal the coordinated immune dysregulation underpinning JDM and provide novel insight into strategies for restoring balance in immune function.

Microfluidics-free single-cell genomics with templated emulsification

Current single-cell RNA-sequencing approaches have limitations that stem from the microfluidic devices or fluid handling steps required for sample processing. We develop a method that does not require specialized microfluidic devices, expertise or hardware. Our approach is based on particle-templated emulsification, which allows single-cell encapsulation and barcoding of cDNA in uniform droplet emulsions with only a vortexer. Particle-templated instant partition sequencing (PIP-seq) accommodates a wide range of emulsification formats, including microwell plates and large-volume conical tubes, enabling thousands of samples or millions of cells to be processed in minutes. We demonstrate that PIP-seq produces high-purity transcriptomes in mouse-human mixing studies, is compatible with multiomics measurements and can accurately characterize cell types in human breast tissue compared to a commercial microfluidic platform. Single-cell transcriptional profiling of mixed phenotype acute leukemia using PIP-seq reveals the emergence of heterogeneity within chemotherapy-resistant cell subsets that were hidden by standard immunophenotyping. PIP-seq is a simple, flexible and scalable next-generation workflow that extends single-cell sequencing to new applications.

Single-cell and spatial analyses reveal a tradeoff between murine mammary proliferation and lineage programs associated with endocrine cues

Although distinct epithelial cell types have been distinguished in glandular tissues such as the mammary gland, the extent of heterogeneity within each cell type and the degree of endocrine control of this diversity across development are incompletely understood. By combining mass cytometry and cyclic immunofluorescence, we define a rich array of murine mammary epithelial cell subtypes associated with puberty, the estrous cycle, and sex. These subtypes are differentially proliferative and spatially segregate distinctly in adult versus pubescent glands. Further, we identify systematic suppression of lineage programs at the protein and RNA levels as a common feature of mammary epithelial expansion during puberty, the estrous cycle, and gestation and uncover a pervasive enrichment of ribosomal protein genes in luminal cells elicited specifically during progesterone-dominant expansionary periods. Collectively, these data expand our knowledge of murine mammary epithelial heterogeneity and connect endocrine-driven epithelial expansion with lineage suppression.

Loss of PPARγ activity characterizes early protumorigenic stromal reprogramming and dictates the therapeutic window of opportunity

Although robustly expressed in the disease-free (DF) breast stroma, CD36 is consistently absent from the stroma surrounding invasive breast cancers (IBCs). In this study, we primarily observed CD36 expression in adipocytes and intralobular capillaries within the DF breast. Larger vessels concentrated in interlobular regions lacked CD36 and were instead marked by the expression of CD31. When evaluated in perilesional capillaries surrounding ductal carcinoma in situ, a nonobligate IBC precursor, CD36 loss was more commonly observed in lesions associated with subsequent IBC. Peroxisome proliferator-activated receptor γ (PPARγ) governs the expression of CD36 and genes involved in differentiation, metabolism, angiogenesis, and inflammation. Coincident with CD36 loss, we observed a dramatic suppression of PPARγ and its target genes in capillary endothelial cells (ECs) and pericytes, which typically surround and support the stability of the capillary endothelium. Factors present in conditioned media from malignant cells repressed PPARγ and its target genes not only in cultured ECs and pericytes but also in adipocytes, which require PPARγ for proper differentiation. In addition, we identified a role for PPARγ in opposing the transition of pericytes toward a tumor-supportive myofibroblast phenotype. In mouse xenograft models, early intervention with rosiglitazone, a PPARγ agonist, demonstrated significant antitumor effects; however, following the development of a palpable tumor, the antitumor effects of rosiglitazone were negated by the repression of PPARγ in the mouse stroma. In summary, PPARγ activity in healthy tissues places several stromal cell types in an antitumorigenic state, directly inhibiting EC proliferation, maintaining adipocyte differentiation, and suppressing the transition of pericytes into tumor-supportive myofibroblasts.

A spatially resolved single-cell genomic atlas of the adult human breast

The adult human breast is comprised of an intricate network of epithelial ducts and lobules that are embedded in connective and adipose tissue1-3. Although most previous studies have focused on the breast epithelial system4-6, many of the non-epithelial cell types remain understudied. Here we constructed the comprehensive Human Breast Cell Atlas (HBCA) at single-cell and spatial resolution. Our single-cell transcriptomics study profiled 714,331 cells from 126 women, and 117,346 nuclei from 20 women, identifying 12 major cell types and 58 biological cell states. These data reveal abundant perivascular, endothelial and immune cell populations, and highly diverse luminal epithelial cell states. Spatial mapping using four different technologies revealed an unexpectedly rich ecosystem of tissue-resident immune cells, as well as distinct molecular differences between ductal and lobular regions. Collectively, these data provide a reference of the adult normal breast tissue for studying mammary biology and diseases such as breast cancer.

Configurational entropy is an intrinsic driver of tissue structural heterogeneity

Tissues comprise ordered arrangements of cells that can be surprisingly disordered in their details. How the properties of single cells and their microenvironment contribute to the balance between order and disorder at the tissue-scale remains poorly understood. Here, we address this question using the self-organization of human mammary organoids as a model. We find that organoids behave like a dynamic structural ensemble at the steady state. We apply a maximum entropy formalism to derive the ensemble distribution from three measurable parameters - the degeneracy of structural states, interfacial energy, and tissue activity (the energy associated with positional fluctuations). We link these parameters with the molecular and microenvironmental factors that control them to precisely engineer the ensemble across multiple conditions. Our analysis reveals that the entropy associated with structural degeneracy sets a theoretical limit to tissue order and provides new insight for tissue engineering, development, and our understanding of disease progression.

Preneoplastic stromal cells promote BRCA1-mediated breast tumorigenesis

Women with germline BRCA1 mutations (BRCA1+/mut) have increased risk for hereditary breast cancer. Cancer initiation in BRCA1+/mut is associated with premalignant changes in breast epithelium; however, the role of the epithelium-associated stromal niche during BRCA1-driven tumor initiation remains unclear. Here we show that the premalignant stromal niche promotes epithelial proliferation and mutant BRCA1-driven tumorigenesis in trans. Using single-cell RNA sequencing analysis of human preneoplastic BRCA1+/mut and noncarrier breast tissues, we show distinct changes in epithelial homeostasis including increased proliferation and expansion of basal-luminal intermediate progenitor cells. Additionally, BRCA1+/mut stromal cells show increased expression of pro-proliferative paracrine signals. In particular, we identify pre-cancer-associated fibroblasts (pre-CAFs) that produce protumorigenic factors including matrix metalloproteinase 3 (MMP3), which promotes BRCA1-driven tumorigenesis in vivo. Together, our findings demonstrate that precancerous stroma in BRCA1+/mut may elevate breast cancer risk through the promotion of epithelial proliferation and an accumulation of luminal progenitor cells with altered differentiation.