CellChat for systematic analysis of cell-cell communication from single-cell transcriptomics.

Recent advances in single-cell sequencing technologies offer an opportunity to explore cell-cell communication in tissues systematically and with reduced bias. A key challenge is integrating known molecular interactions and measurements into a framework to identify and analyze complex cell-cell communication networks. Previously, we developed a computational tool, named CellChat, that infers and analyzes cell-cell communication networks from single-cell transcriptomic data within an easily interpretable framework. CellChat quantifies the signaling communication probability between two cell groups using a simplified mass-action-based model, which incorporates the core interaction between ligands and receptors with multisubunit structure along with modulation by cofactors. Importantly, CellChat performs a systematic and comparative analysis of cell-cell communication using a variety of quantitative metrics and machine-learning approaches. CellChat v2 is an updated version that includes additional comparison functionalities, an expanded database of ligand-receptor pairs along with rich functional annotations, and an Interactive CellChat Explorer. Here we provide a step-by-step protocol for using CellChat v2 on single-cell transcriptomic data, including inference and analysis of cell-cell communication from one dataset and identification of altered intercellular communication, signals and cell populations from different datasets across biological conditions. The R implementation of CellChat v2 toolkit and its tutorials together with the graphic outputs are available at https://github.com/jinworks/CellChat . This protocol typically takes ~5 min depending on dataset size and requires a basic understanding of R and single-cell data analysis but no specialized bioinformatics training for its implementation.

Hydrogel crosslinking modulates macrophages, fibroblasts, and their communication, during wound healing.

Biomaterial wound dressings, such as hydrogels, interact with host cells to regulate tissue repair. This study investigates how crosslinking of gelatin-based hydrogels influences immune and stromal cell behavior and wound healing in female mice. We observe that softer, lightly crosslinked hydrogels promote greater cellular infiltration and result in smaller scars compared to stiffer, heavily crosslinked hydrogels. Using single-cell RNA sequencing, we further show that heavily crosslinked hydrogels increase inflammation and lead to the formation of a distinct macrophage subpopulation exhibiting signs of oxidative activity and cell fusion. Conversely, lightly crosslinked hydrogels are more readily taken up by macrophages and integrated within the tissue. The physical properties differentially affect macrophage and fibroblast interactions, with heavily crosslinked hydrogels promoting pro-fibrotic fibroblast activity that drives macrophage fusion through RANKL signaling. These findings suggest that tuning the physical properties of hydrogels can guide cellular responses and improve healing, offering insights for designing better biomaterials for wound treatment.

Three-dimensional genome architecture persists in a 52,000-year-old woolly mammoth skin sample.

Analyses of ancient DNA typically involve sequencing the surviving short oligonucleotides and aligning to genome assemblies from related, modern species. Here, we report that skin from a female woolly mammoth (†Mammuthus primigenius) that died 52,000 years ago retained its ancient genome architecture. We use PaleoHi-C to map chromatin contacts and assemble its genome, yielding 28 chromosome-length scaffolds. Chromosome territories, compartments, loops, Barr bodies, and inactive X chromosome (Xi) superdomains persist. The active and inactive genome compartments in mammoth skin more closely resemble Asian elephant skin than other elephant tissues. Our analyses uncover new biology. Differences in compartmentalization reveal genes whose transcription was potentially altered in mammoths vs. elephants. Mammoth Xi has a tetradic architecture, not bipartite like human and mouse. We hypothesize that, shortly after this mammoth's death, the sample spontaneously freeze-dried in the Siberian cold, leading to a glass transition that preserved subfossils of ancient chromosomes at nanometer scale.

Integrated Single-Cell Analysis Reveals Spatially and Temporally Dynamic Heterogeneity in Fibroblast States during Wound Healing.

Wound healing is a dynamic process over temporal and spatial scales. Key to repair outcomes are fibroblasts; yet, how they modulate healing across time and in different wound regions remains incompletely understood. By integrating single-cell RNA-sequencing datasets of mouse skin and wounds, we infer that fibroblasts are the most transcriptionally dynamic skin-resident cells, evolving during postnatal skin maturation and rapidly after injury toward distinct late scar states. We show that transcriptional dynamics in fibroblasts are largely driven by genes encoding extracellular matrix and signaling factors. Lineage trajectory inference and spatial gene mapping reveal that Prg4-expressing fibroblasts transiently emerge along early wound edges. Within days, they become replaced by long-lasting and likely noninterconverting fibroblast populations, including Col25a1-expressing and Pamr1-expressing fibroblasts that occupy subepidermal and deep scar regions, respectively, where they engage in reciprocal signaling with immune cells. Signaling inference shows that fibroblast-immune crosstalk repeatedly uses some signaling pathways across wound healing time, whereas use of other signaling pathways is time and space limited. Collectively, we uncovered high transcriptional plasticity by wound fibroblasts, with early states transiently forming distinct microniches along wound edges and in the fascia, followed by stable states that stratify scar tissue into molecularly dissimilar upper and lower layers.

Dermal TRPV1 innervations engage a macrophage- and fibroblast-containing pathway to activate hair growth in mice.

Pain, detected by nociceptors, is an integral part of injury, yet whether and how it can impact tissue physiology and recovery remain understudied. Here, we applied chemogenetics in mice to locally activate dermal TRPV1 innervations in naive skin and found that it triggered new regenerative cycling by dormant hair follicles (HFs). This was preceded by rapid apoptosis of dermal macrophages, mediated by the neuropeptide calcitonin gene-related peptide (CGRP). TRPV1 activation also triggered a macrophage-dependent induction of osteopontin (Spp1)-expressing dermal fibroblasts. The neuropeptide CGRP and the extracellular matrix protein Spp1 were required for the nociceptor-triggered hair growth. Finally, we showed that epidermal abrasion injury induced Spp1-expressing dermal fibroblasts and hair growth via a TRPV1 neuron and CGRP-dependent mechanism. Collectively, these data demonstrated a role for TRPV1 nociceptors in orchestrating a macrophage and fibroblast-supported mechanism to promote hair growth and enabling the efficient restoration of this mechano- and thermo-protective barrier after wounding.

Parsing patterns: Emerging roles of tissue self-organization in health and disease.

Patterned morphologies, such as segments, spirals, stripes, and spots, frequently emerge during embryogenesis through self-organized coordination between cells. Yet, complex patterns also emerge in adults, suggesting that the capacity for spontaneous self-organization is a ubiquitous property of biological tissues. We review current knowledge on the principles and mechanisms of self-organized patterning in embryonic tissues and explore how these principles and mechanisms apply to adult tissues that exhibit features of patterning. We discuss how and why spontaneous pattern generation is integral to homeostasis and healing of tissues, illustrating it with examples from regenerative biology. We examine how aberrant self-organization underlies diverse pathological states, including inflammatory skin disorders and tumors. Lastly, we posit that based on such blueprints, targeted engineering of pattern-driving molecular circuits can be leveraged for synthetic biology and the generation of organoids with intricate patterns.

Emx2 underlies the development and evolution of marsupial gliding membranes.

Phenotypic variation among species is a product of evolutionary changes to developmental programs1,2. However, how these changes generate novel morphological traits remains largely unclear. Here we studied the genomic and developmental basis of the mammalian gliding membrane, or patagium-an adaptative trait that has repeatedly evolved in different lineages, including in closely related marsupial species. Through comparative genomic analysis of 15 marsupial genomes, both from gliding and non-gliding species, we find that the Emx2 locus experienced lineage-specific patterns of accelerated cis-regulatory evolution in gliding species. By combining epigenomics, transcriptomics and in-pouch marsupial transgenics, we show that Emx2 is a critical upstream regulator of patagium development. Moreover, we identify different cis-regulatory elements that may be responsible for driving increased Emx2 expression levels in gliding species. Lastly, using mouse functional experiments, we find evidence that Emx2 expression patterns in gliders may have been modified from a pre-existing program found in all mammals. Together, our results suggest that patagia repeatedly originated through a process of convergent genomic evolution, whereby regulation of Emx2 was altered by distinct cis-regulatory elements in independently evolved species. Thus, different regulatory elements targeting the same key developmental gene may constitute an effective strategy by which natural selection has harnessed regulatory evolution in marsupial genomes to generate phenotypic novelty.

mTORC2-AKT signaling to PFKFB2 activates glycolysis that enhances stemness and tumorigenicity of intestinal epithelial cells.

Although elevated glycolysis has been widely recognized as a hallmark for highly proliferating cells like stem cells and cancer, its regulatory mechanisms are still being updated. Here, we found a previously unappreciated mechanism of mammalian target of rapamycin complex 2 (mTORC2) in regulating glycolysis in intestinal stem cell maintenance and cancer progression. mTORC2 key subunits expression levels and its kinase activity were specifically upregulated in intestinal stem cells, mouse intestinal tumors, and human colorectal cancer (CRC) tissues. Genetic ablation of its key scaffolding protein Rictor in both mouse models and cell lines revealed that mTORC2 played an important role in promoting intestinal stem cell proliferation and self-renewal. Moreover, utilizing mouse models and organoid culture, mTORC2 loss of function was shown to impair growth of gut adenoma and tumor organoids. Based on these findings, we performed RNA-seq and noticed significant metabolic reprogramming in Rictor conditional knockout mice. Among all the pathways, carbohydrate metabolism was most profoundly altered, and further studies demonstrated that mTORC2 promoted glycolysis in intestinal epithelial cells. Most importantly, we showed that a rate-limiting enzyme in regulating glycolysis, 6-phosphofructo-2-kinase (PFKFB2), was a direct target for the mTORC2-AKT signaling. PFKFB2 was phosphorylated upon mTORC2 activation, but not mTORC1, and this process was AKT-dependent. Together, this study has identified a novel mechanism underlying mTORC2 activated glycolysis, offering potential therapeutic targets for treating CRC.

Excessive nucleic acid R-loops induce mitochondria-dependent epithelial cell necroptosis and drive spontaneous intestinal inflammation.

Oxidative stress, which can be activated by a variety of environmental risk factors, has been implicated as an important pathogenic factor for inflammatory bowel disease (IBD). However, how oxidative stress drives IBD onset remains elusive. Here, we found that oxidative stress was strongly activated in inflamed tissues from both ulcerative colitis patients and Crohn's disease patients, and it caused nuclear-to-cytosolic TDP-43 transport and a reduction in the TDP-43 protein level. To investigate the function of TDP-43 in IBD, we inducibly deleted exons 2 to 3 of Tardbp (encoding Tdp-43) in mouse intestinal epithelium, which disrupted its nuclear localization and RNA-processing function. The deletion gave rise to spontaneous intestinal inflammation by inducing epithelial cell necroptosis. Suppression of the necroptotic pathway with deletion of Mlkl or the RIP1 inhibitor Nec-1 rescued colitis phenotypes. Mechanistically, disruption of nuclear TDP-43 caused excessive R-loop accumulation, which triggered DNA damage and genome instability and thereby induced PARP1 hyperactivation, leading to subsequent NAD+ depletion and ATP loss, consequently activating mitochondrion-dependent necroptosis in intestinal epithelial cells. Importantly, restoration of cellular NAD+ levels with NAD+ or NMN supplementation, as well as suppression of ALKBH7, an α-ketoglutarate dioxygenase in mitochondria, rescued TDP-43 deficiency-induced cell death and intestinal inflammation. Furthermore, TDP-43 protein levels were significantly inversely correlated with γ-H2A.X and p-MLKL levels in clinical IBD samples, suggesting the clinical relevance of TDP-43 deficiency-induced mitochondrion-dependent necroptosis. Taken together, these findings identify a unique pathogenic mechanism that links oxidative stress to intestinal inflammation and provide a potent and valid strategy for IBD intervention.

Delayed Collagen Production without Myofibroblast Formation Contributes to Reduced Scarring in Adult Skin Microwounds.

In adult mammals, wound healing predominantly follows a fibrotic pathway, culminating in scar formation. However, cutaneous microwounds generated through fractional photothermolysis, a modality that produces a constellation of microthermal zones, exhibit a markedly different healing trajectory. Our study delineates the cellular attributes of these microthermal zones, underscoring a temporally limited, subclinical inflammatory milieu concomitant with rapid re-epithelialization within 24 hours. This wound closure is facilitated by the activation of genes associated with keratinocyte migration and differentiation. In contrast to macrothermal wounds, which predominantly heal through a robust myofibroblast-mediated collagen deposition, microthermal zones are characterized by absence of wound contraction and feature delayed collagen remodeling, initiating 5-6 weeks after injury. This distinct wound healing is characterized by a rapid re-epithelialization process and a muted inflammatory response, which collectively serve to mitigate excessive myofibroblast activation. Furthermore, we identify an initial reparative phase characterized by a heterogeneous extracellular matrix protein composition, which precedes the delayed collagen remodeling. These findings extend our understanding of cutaneous wound healing and may have significant implications for the optimization of therapeutic strategies aimed at mitigating scar formation.

Studying Hair Growth Cycle and its Effects on Mouse Skin.

Researchers should be aware that hair growth cycle drives prominent molecular, cellular, and morphological changes to the entire skin. Thus, hair growth constitutes a major experimental variable that influences the interpretation of dermatological studies. Hair growth in mice is neither asynchronous nor fully synchronized; rather, it occurs in waves that dynamically propagate across the skin. In consequence, any given area of mouse skin can contain hair follicles in different stages of the cycle in close physical proximity. Furthermore, hair growth waves in mice are initiated by probabilistic events at different time points and across stochastic locations. The consequence of such stochasticity is that precise patterns of hair growth waves differ from mouse to mouse, even in littermates of the same sex. However, such physiological stochasticity is commonly misconstrued as a significant hair growth phenotype in mutant mice or in drug-treated mice. The purpose of this article is to provide a set of guidelines for designing reliably interpretable murine studies on hair growth and to highlight key experimental caveats to be avoided. It also informs on how to account for and minimize the impact of physiological hair cycle differences when designing and interpreting nonhair growth dermatological studies in mice.

A Roadmap for a Consensus Human Skin Cell Atlas and Single-Cell Data Standardization.

Single-cell technologies have become essential to driving discovery in both basic and translational investigative dermatology. Despite the multitude of available datasets, a central reference atlas of normal human skin, which can serve as a reference resource for skin cell types, cell states, and their molecular signatures, is still lacking. For any such atlas to receive broad acceptance, participation by many investigators during atlas construction is an essential prerequisite. As part of the Human Cell Atlas project, we have assembled a Skin Biological Network to build a consensus Human Skin Cell Atlas and outline a roadmap toward that goal. We define the drivers of skin diversity to be considered when selecting sequencing datasets for the atlas and list practical hurdles during skin sampling that can result in data gaps and impede comprehensive representation and technical considerations for tissue processing and computational analysis, the accounting for which should minimize biases in cell type enrichments and exclusions and decrease batch effects. By outlining our goals for Atlas 1.0, we discuss how it will uncover new aspects of skin biology.

Dynamic interplay between IL-1 and WNT pathways in regulating dermal adipocyte lineage cells during skin development and wound regeneration.

Dermal adipocyte lineage cells are highly plastic and can undergo reversible differentiation and dedifferentiation in response to various stimuli. Using single-cell RNA sequencing of developing or wounded mouse skin, we classify dermal fibroblasts (dFBs) into distinct non-adipogenic and adipogenic cell states. Cell differentiation trajectory analyses identify IL-1-NF-κB and WNT-ß-catenin as top signaling pathways that positively and negatively associate with adipogenesis, respectively. Upon wounding, activation of adipocyte progenitors and wound-induced adipogenesis are mediated in part by neutrophils through the IL-1R-NF-κB-CREB signaling axis. In contrast, WNT activation, by WNT ligand and/or ablation of Gsk3, inhibits the adipogenic potential of dFBs but promotes lipolysis and dedifferentiation of mature adipocytes, contributing to myofibroblast formation. Finally, sustained WNT activation and inhibition of adipogenesis is seen in human keloids. These data reveal molecular mechanisms underlying the plasticity of dermal adipocyte lineage cells, defining potential therapeutic targets for defective wound healing and scar formation.

Signalling by senescent melanocytes hyperactivates hair growth.

Niche signals maintain stem cells in a prolonged quiescence or transiently activate them for proper regeneration1. Altering balanced niche signalling can lead to regenerative disorders. Melanocytic skin nevi in human often display excessive hair growth, suggesting hair stem cell hyperactivity. Here, using genetic mouse models of nevi2,3, we show that dermal clusters of senescent melanocytes drive epithelial hair stem cells to exit quiescence and change their transcriptome and composition, potently enhancing hair renewal. Nevus melanocytes activate a distinct secretome, enriched for signalling factors. Osteopontin, the leading nevus signalling factor, is both necessary and sufficient to induce hair growth. Injection of osteopontin or its genetic overexpression is sufficient to induce robust hair growth in mice, whereas germline and conditional deletions of either osteopontin or CD44, its cognate receptor on epithelial hair cells, rescue enhanced hair growth induced by dermal nevus melanocytes. Osteopontin is overexpressed in human hairy nevi, and it stimulates new growth of human hair follicles. Although broad accumulation of senescent cells, such as upon ageing or genotoxic stress, is detrimental for the regenerative capacity of tissue4, we show that signalling by senescent cell clusters can potently enhance the activity of adjacent intact stem cells and stimulate tissue renewal. This finding identifies senescent cells and their secretome as an attractive therapeutic target in regenerative disorders.

IRX5 promotes DNA damage repair and activation of hair follicle stem cells.

The molecular mechanisms allowing hair follicles to periodically activate their stem cells (HFSCs) are incompletely characterized. Here, we identify the transcription factor IRX5 as a promoter of HFSC activation. Irx5-/- mice have delayed anagen onset, with increased DNA damage and diminished HFSC proliferation. Open chromatin regions form near cell cycle progression and DNA damage repair genes in Irx5-/- HFSCs. DNA damage repair factor BRCA1 is an IRX5 downstream target. Inhibition of FGF kinase signaling partially rescues the anagen delay in Irx5-/- mice, suggesting that the Irx5-/- HFSC quiescent phenotype is partly due to failure to suppress Fgf18 expression. Interfollicular epidermal stem cells also show decreased proliferation and increased DNA damage in Irx5-/-mice. Consistent with a role for IRX5 as a promoter of DNA damage repair, we find that IRX genes are upregulated in many cancer types and that there is a correlation between IRX5 and BRCA1 expression in breast cancer.

Convergent deployment of ancestral functions during the evolution of mammalian flight membranes.

Lateral flight membranes, or patagia, have evolved repeatedly in diverse mammalian lineages. While little is known about patagium development, its recurrent evolution may suggest a shared molecular basis. By combining transcriptomics, developmental experiments, and mouse transgenics, we demonstrate that lateral Wnt5a expression in the marsupial sugar glider (Petaurus breviceps) promotes the differentiation of its patagium primordium. We further show that this function of Wnt5a reprises ancestral roles in skin morphogenesis predating mammalian flight and has been convergently used during patagium evolution in eutherian bats. Moreover, we find that many genes involved in limb development have been redeployed during patagium outgrowth in both the sugar glider and bat. Together, our findings reveal that deeply conserved genetic toolkits contribute to the evolutionary transition to flight in mammals.