Fibroblast inflammatory priming determines regenerative versus fibrotic skin repair in reindeer.

Adult mammalian skin wounds heal by forming fibrotic scars. We report that full-thickness injuries of reindeer antler skin (velvet) regenerate, whereas back skin forms fibrotic scar. Single-cell multi-omics reveal that uninjured velvet fibroblasts resemble human fetal fibroblasts, whereas back skin fibroblasts express inflammatory mediators mimicking pro-fibrotic adult human and rodent fibroblasts. Consequently, injury elicits site-specific immune responses: back skin fibroblasts amplify myeloid infiltration and maturation during repair, whereas velvet fibroblasts adopt an immunosuppressive phenotype that restricts leukocyte recruitment and hastens immune resolution. Ectopic transplantation of velvet to scar-forming back skin is initially regenerative, but progressively transitions to a fibrotic phenotype akin to the scarless fetal-to-scar-forming transition reported in humans. Skin regeneration is diminished by intensifying, or enhanced by neutralizing, these pathologic fibroblast-immune interactions. Reindeer represent a powerful comparative model for interrogating divergent wound healing outcomes, and our results nominate decoupling of fibroblast-immune interactions as a promising approach to mitigate scar.

Olig2 defines a subset of neural stem cells that produce specific olfactory bulb interneuron subtypes in the subventricular zone of adult mice.

Distinct neural stem cells (NSCs) reside in different regions of the subventricular zone (SVZ) and generate multiple olfactory bulb (OB) interneuron subtypes in the adult brain. However, the molecular mechanisms underlying such NSC heterogeneity remain largely unknown. Here, we show that the basic helix-loop-helix transcription factor Olig2 defines a subset of NSCs in the early postnatal and adult SVZ. Olig2-expressing NSCs exist broadly but are most enriched in the ventral SVZ along the dorsoventral axis complementary to dorsally enriched Gsx2-expressing NSCs. Comparisons of Olig2-expressing NSCs from early embryonic to adult stages using single cell transcriptomics reveal stepwise developmental changes in their cell cycle and metabolic properties. Genetic studies further show that cross-repression contributes to the mutually exclusive expression of Olig2 and Gsx2 in NSCs/progenitors during embryogenesis, but that their expression is regulated independently from each other in adult NSCs. Finally, lineage-tracing and conditional inactivation studies demonstrate that Olig2 plays an important role in the specification of OB interneuron subtypes. Altogether, our study demonstrates that Olig2 defines a unique subset of adult NSCs enriched in the ventral aspect of the adult SVZ.

IL-1/MyD88-Dependent G-CSF and IL-6 Secretion Mediates Postburn Anemia.

The anemia of critical illness (ACI) is a nearly universal pathophysiological consequence of burn injury and a primary reason burn patients require massive quantities of transfused blood. Inflammatory processes are expected to drive postburn ACI and prevent meaningful erythropoietic stimulation through iron or erythropoietin supplementation, but to this day no specific inflammatory pathways have been identified as a critical mechanism. In this study, we examined whether secretion of G-CSF and IL-6 mediates distinct features of postburn ACI and interrogated inflammatory mechanisms that could be responsible for their secretion. Our analysis of mouse and human skin samples identified the burn wound as a primary source of G-CSF and IL-6 secretion. We show that G-CSF and IL-6 are secreted independently through an IL-1/MyD88-dependent mechanism, and we ruled out TLR2 and TLR4 as critical receptors. Our results indicate that IL-1/MyD88-dependent G-CSF secretion plays a key role in impairing medullary erythropoiesis and IL-6 secretion plays a key role in limiting the access of erythroid cells to iron. Importantly, we found that IL-1α/ß neutralizing Abs broadly attenuated features of postburn ACI that could be attributed to G-CSF or IL-6 secretion and rescued deficits of circulating RBC counts, hemoglobin, and hematocrit caused by burn injury. We conclude that wound-based IL-1/MyD88 signaling mediates postburn ACI through induction of G-CSF and IL-6 secretion.

HoxD transcription factors define monosynaptic sensory-motor specificity in the developing spinal cord.

The specificity of monosynaptic connections between proprioceptive sensory neurons and their recipient spinal motor neurons depends on multiple factors, including motor neuron positioning and dendrite morphology, axon projection patterns of proprioceptive sensory neurons in the spinal cord, and the ligand-receptor molecules involved in cell-to-cell recognition. However, with few exceptions, the transcription factors engaged in this process are poorly characterized. Here, we show that members of the HoxD family of transcription factors play a crucial role in the specificity of monosynaptic sensory-motor connections. Mice lacking Hoxd9, Hoxd10 and Hoxd11 exhibit defects in locomotion but have no obvious defects in motor neuron positioning or dendrite morphology through the medio-lateral and rostro-caudal axes. However, we found that quadriceps motor neurons in these mice show aberrant axon development and receive inappropriate inputs from proprioceptive sensory axons innervating the obturator muscle. These genetic studies demonstrate that the HoxD transcription factors play an integral role in the synaptic specificity of monosynaptic sensory-motor connections in the developing spinal cord.

Physical interactions between Gsx2 and Ascl1 balance progenitor expansion versus neurogenesis in the mouse lateral ganglionic eminence.

The Gsx2 homeodomain transcription factor promotes neural progenitor identity in the lateral ganglionic eminence (LGE), despite upregulating the neurogenic factor Ascl1. How this balance in maturation is maintained is unclear. Here, we show that Gsx2 and Ascl1 are co-expressed in subapical progenitors that have unique transcriptional signatures in LGE ventricular zone (VZ) cells. Moreover, whereas Ascl1 misexpression promotes neurogenesis in dorsal telencephalic progenitors, the co-expression of Gsx2 with Ascl1 inhibits neurogenesis. Using luciferase assays, we found that Gsx2 reduces the ability of Ascl1 to activate gene expression in a dose-dependent and DNA binding-independent manner. Furthermore, Gsx2 physically interacts with the basic helix-loop-helix (bHLH) domain of Ascl1, and DNA-binding assays demonstrated that this interaction interferes with the ability of Ascl1 to bind DNA. Finally, we modified a proximity ligation assay for tissue sections and found that Ascl1-Gsx2 interactions are enriched within LGE VZ progenitors, whereas Ascl1-Tcf3 (E-protein) interactions predominate in the subventricular zone. Thus, Gsx2 contributes to the balance between progenitor maintenance and neurogenesis by physically interacting with Ascl1, interfering with its DNA binding and limiting neurogenesis within LGE progenitors.

Cross talk between LAM cells and fibroblasts may influence alveolar epithelial cell behavior in lymphangioleiomyomatosis.

Lymphangioleiomyomatosis (LAM) is a female-specific cystic lung disease in which tuberous sclerosis complex 2 (TSC2)-deficient LAM cells, LAM-associated fibroblasts (LAFs), and other cell types infiltrate the lungs. LAM lesions can be associated with type II alveolar epithelial (AT2) cells. We hypothesized that the behavior of AT2 cells in LAM is influenced locally by LAFs. We tested this hypothesis in the patient samples and in vitro. In human LAM lung, nodular AT2 cells show enhanced proliferation when compared with parenchymal AT2 cells, demonstrated by increased Ki67 expression. Furthermore, nodular AT2 cells express proteins associated with epithelial activation in other disease states including matrix metalloproteinase 7, and fibroblast growth factor 7 (FGF7). In vitro, LAF-conditioned medium is mitogenic and positively chemotactic for epithelial cells, increases the rate of epithelial repair, and protects against apoptosis. In vitro, LAM patient-derived TSC2 null cells cocultured with LAFs upregulate LAF expression of the epithelial chemokine and mitogen FGF7, a potential mediator of fibroblast-epithelial cross talk, in a mechanistic target of rapamycin (mTOR)-dependent manner. In a novel in vitro model of LAM, ex vivo cultured LAM lung-derived microtissues promote both epithelial migration and adhesion. Our findings suggest that AT2 cells in LAM display a proliferative, activated phenotype and fibroblast accumulation following LAM cell infiltration into the parenchyma contributes to this change in AT2 cell behavior. Fibroblast-derived FGF7 may contribute to the cross talk between LAFs and hyperplastic epithelium in vivo, but does not appear to be the main driver of the effects of LAFs on epithelial cells in vitro.

Combinatorial Transcriptional Profiling of Mouse and Human Enteric Neurons Identifies Shared and Disparate Subtypes In Situ.

BACKGROUND & AIMS: The enteric nervous system (ENS) coordinates essential intestinal functions through the concerted action of diverse enteric neurons (ENs). However, integrated molecular knowledge of EN subtypes is lacking. To compare human and mouse ENs, we transcriptionally profiled healthy ENS from adult humans and mice. We aimed to identify transcripts marking discrete neuron subtypes and visualize conserved EN subtypes for humans and mice in multiple bowel regions. METHODS: Human myenteric ganglia and adjacent smooth muscle were isolated by laser-capture microdissection for RNA sequencing. Ganglia-specific transcriptional profiles were identified by computationally subtracting muscle gene signatures. Nuclei from mouse myenteric neurons were isolated and subjected to single-nucleus RNA sequencing, totaling more than 4 billion reads and 25,208 neurons. Neuronal subtypes were defined using mouse single-nucleus RNA sequencing data. Comparative informatics between human and mouse data sets identified shared EN subtype markers, which were visualized in situ using hybridization chain reaction. RESULTS: Several EN subtypes in the duodenum, ileum, and colon are conserved between humans and mice based on orthologous gene expression. However, some EN subtype-specific genes from mice are expressed in completely distinct morphologically defined subtypes in humans. In mice, we identified several neuronal subtypes that stably express gene modules across all intestinal segments, with graded, regional expression of 1 or more marker genes. CONCLUSIONS: Our combined transcriptional profiling of human myenteric ganglia and mouse EN provides a rich foundation for developing novel intestinal therapeutics. There is congruency among some EN subtypes, but we note multiple species differences that should be carefully considered when relating findings from mouse ENS research to human gastrointestinal studies.

Maturation of heart valve cell populations during postnatal remodeling.

Heart valve cells mediate extracellular matrix (ECM) remodeling during postnatal valve leaflet stratification, but phenotypic and transcriptional diversity of valve cells in development is largely unknown. Single cell analysis of mouse heart valve cells was used to evaluate cell heterogeneity during postnatal ECM remodeling and leaflet morphogenesis. The transcriptomic analysis of single cells from postnatal day (P)7 and P30 murine aortic (AoV) and mitral (MV) heart valves uncovered distinct subsets of melanocytes, immune and endothelial cells present at P7 and P30. By contrast, interstitial cell populations are different from P7 to P30. P7 valve leaflets exhibit two distinct collagen- and glycosaminoglycan-expressing interstitial cell clusters, and prevalent ECM gene expression. At P30, four interstitial cell clusters are apparent with leaflet specificity and differential expression of complement factors, ECM proteins and osteogenic genes. This initial transcriptomic analysis of postnatal heart valves at single cell resolution demonstrates that subpopulations of endothelial and immune cells are relatively constant throughout postnatal development, but interstitial cell subpopulations undergo changes in gene expression and cellular functions in primordial and mature valves.

Single-Cell Transcriptomic Analysis Identifies a Unique Pulmonary Lymphangioleiomyomatosis Cell.

Rationale: Lymphangioleiomyomatosis (LAM) is a metastatic neoplasm of reproductive-age women associated with mutations in tuberous sclerosis complex genes. LAM causes cystic remodeling of the lung and progressive respiratory failure. The sources and cellular characteristics of LAM cells underlying disease pathogenesis remain elusive.Objectives: Identification and characterization of LAM cells in human lung and uterus using a single-cell approach.Methods: Single-cell and single-nuclei RNA sequencing on LAM (n = 4) and control (n = 7) lungs, immunofluorescence confocal microscopy, ELISA, and aptamer proteomics were used to identify and validate LAMCORE cells and secreted biomarkers, predict cellular origins, and define molecular and cellular networks in LAM.Measurements and Main Results: A unique cell type termed LAMCORE was identified, which was distinct from, but closely related to, lung mesenchymal cells. LAMCORE cells expressing signature genes included known LAM markers such as PMEL, FIGF, CTSK, and MLANA and novel biomarkers validated by aptamer screening, ELISA, and immunofluorescence microscopy. LAM cells in lung and uterus are morphologically indistinguishable and share similar gene expression profiles and biallelic TSC2 mutations, supporting a potential uterine origin for the LAMCORE cell. Effects of LAM on resident pulmonary cell types indicated recruitment and activation of lymphatic endothelial cells.Conclusions: A unique population of LAMCORE cells was identified in lung and uterus of patients with LAM, sharing close transcriptomic identity. LAM cell selective markers, secreted biomarkers, and the predicted cellular molecular features provide new insights into the signaling and transcriptional programs that may serve as diagnostic markers and therapeutic targets to influence the pathogenesis of LAM.

Single-Cell Profiling of AKI in a Murine Model Reveals Novel Transcriptional Signatures, Profibrotic Phenotype, and Epithelial-to-Stromal Crosstalk.

BACKGROUND: Current management of AKI, a potentially fatal disorder that can also initiate or exacerbate CKD, is merely supportive. Therefore, deeper understanding of the molecular pathways perturbed in AKI is needed to identify targets with potential to lead to improved treatment. METHODS: We performed single-cell RNA sequencing (scRNA-seq) with the clinically relevant unilateral ischemia-reperfusion murine model of AKI at days 1, 2, 4, 7, 11, and 14 after AKI onset. Using real-time quantitative PCR, immunofluorescence, Western blotting, and both chromogenic and single-molecule in situ hybridizations, we validated AKI signatures in multiple experiments. RESULTS: Our findings show the time course of changing gene expression patterns for multiple AKI stages and all renal cell types. We observed elevated expression of crucial injury response factors-including kidney injury molecule-1 (Kim1), lipocalin 2 (Lcn2), and keratin 8 (Krt8)-and of several novel genes (Ahnak, Sh3bgrl3, and Col18a1) not previously examined in kidney pathologies. AKI induced proximal tubule dedifferentiation, with a pronounced nephrogenic signature represented by Sox4 and Cd24a. Moreover, AKI caused the formation of "mixed-identity cells" (expressing markers of different renal cell types) that are normally seen only during early kidney development. The injured tubules acquired a proinflammatory and profibrotic phenotype; moreover, AKI dramatically modified ligand-receptor crosstalk, with potential pathologic epithelial-to-stromal interactions. Advancing age in AKI onset was associated with maladaptive response and kidney fibrosis. CONCLUSIONS: The scRNA-seq, comprehensive, cell-specific profiles provide a valuable resource for examining molecular pathways that are perturbed in AKI. The results fully define AKI-associated dedifferentiation programs, potential pathologic ligand-receptor crosstalk, novel genes, and the improved injury response in younger mice, and highlight potential targets of kidney injury.

Psychrophilic proteases dramatically reduce single-cell RNA-seq artifacts: a molecular atlas of kidney development.

Single-cell RNA-seq is a powerful technique. Nevertheless, there are important limitations, including the technical challenges of breaking down an organ or tissue into a single-cell suspension. Invariably, this has required enzymatic incubation at 37°C, which can be expected to result in artifactual changes in gene expression patterns. Here, we describe a dissociation method that uses a protease with high activity in the cold, purified from a psychrophilic microorganism. The entire procedure is carried out at 6°C or colder, at which temperature mammalian transcriptional machinery is largely inactive, thereby effectively 'freezing in' the in vivo gene expression patterns. To test this method, we carried out RNA-seq on 20,424 single cells from postnatal day 1 mouse kidneys, comparing the results of the psychrophilic protease method with procedures using 37°C incubation. We show that the cold protease method provides a great reduction in gene expression artifacts. In addition, the results produce a single-cell resolution gene expression atlas of the newborn mouse kidney, an interesting time in development when mature nephrons are present yet nephrogenesis remains extremely active.

Reduced Abd-B Hox function during kidney development results in lineage infidelity.

Hox genes can function as key drivers of segment identity, with Hox mutations in Drosophila often resulting in dramatic homeotic transformations. In addition, however, they can serve other essential functions. In mammals, the study of Hox gene roles in development is complicated by the presence of four Hox clusters with a total of 39 genes showing extensive functional overlap. In this study, in order to better understand shared core Hox functions, we examined kidney development in mice with frameshift mutations of multiple Abd-B type Hox genes. The resulting phenotypes included dramatically reduced branching morphogenesis of the ureteric bud, premature depletion of nephron progenitors and abnormal development of the stromal compartment. Most unexpected, however, we also observed a cellular level lineage infidelity in nephron segments. Scattered cells within the proximal tubules, for example, expressed genes normally expressed only in collecting ducts. Multiple combinations of inappropriate nephron segment specific marker expression were found. In some cases, cells within a tubule showed incorrect identity, while in other cases cells showed ambiguous character, with simultaneous expression of genes associated with more than one nephron segment. These results give evidence that Hox genes have an overlapping core function at the cellular level in driving and/or maintaining correct differentiation decisions.

Cross-platform single cell analysis of kidney development shows stromal cells express Gdnf.

The developing kidney provides a useful model for study of the principles of organogenesis. In this report we use three independent platforms, Drop-Seq, Chromium 10x Genomics and Fluidigm C1, to carry out single cell RNA-Seq (scRNA-Seq) analysis of the E14.5 mouse kidney. Using the software AltAnalyze, in conjunction with the unsupervised approach ICGS, we were unable to identify and confirm the presence of 16 distinct cell populations during this stage of active nephrogenesis. Using a novel integrative supervised computational strategy, we were able to successfully harmonize and compare the cell profiles across all three technological platforms. Analysis of possible cross compartment receptor/ligand interactions identified the nephrogenic zone stroma as a source of GDNF. This was unexpected because the cap mesenchyme nephron progenitors had been thought to be the sole source of GDNF, which is a key driver of branching morphogenesis of the collecting duct system. The expression of Gdnf by stromal cells was validated in several ways, including Gdnf in situ hybridization combined with immunohistochemistry for SIX2, and marker of nephron progenitors, and MEIS1, a marker of stromal cells. Finally, the single cell gene expression profiles generated in this study confirmed and extended previous work showing the presence of multilineage priming during kidney development. Nephron progenitors showed stochastic expression of genes associated with multiple potential differentiation lineages.

Recombineering-based dissection of flanking and paralogous Hox gene functions in mouse reproductive tracts.

Hox genes are key regulators of development. In mammals, the study of these genes is greatly confounded by their large number, overlapping functions and interspersed shared enhancers. Here, we describe the use of a novel recombineering strategy to introduce simultaneous frameshift mutations into the flanking Hoxa9, Hoxa10 and Hoxa11 genes, as well as their paralogs on the HoxD cluster. The resulting Hoxa9,10,11 mutant mice displayed dramatic synergistic homeotic transformations of the reproductive tracts, with the uterus anteriorized towards oviduct and the vas deferens anteriorized towards epididymis. The Hoxa9,10,11 mutant mice also provided a genetic setting that allowed the discovery of Hoxd9,10,11 redundant reproductive tract patterning function. Both shared and distinct Hox functions were defined. Hoxd9,10,11 play a crucial role in the regulation of uterine immune function. Non-coding non-polyadenylated RNAs were among the key Hox targets, with dramatic downregulation in mutants. We observed Hox cross-regulation of transcription and splicing. In addition, we observed a surprising anti-dogmatic apparent posteriorization of the uterine epithelium. In caudal regions of the uterus, the normal simple columnar epithelium flanking the lumen was replaced by a pseudostratified transitional epithelium, normally found near the more posterior cervix. These results identify novel molecular functions of Hox genes in the development of the male and female reproductive tracts.

Key pathways regulated by HoxA9,10,11/HoxD9,10,11 during limb development.

BACKGROUND: The 39 mammalian Hox genes show problematic patterns of functional overlap. In order to more fully define the developmental roles of Hox genes it is necessary to remove multiple combinations of paralogous and flanking genes. In addition, the downstream molecular pathways regulated by Hox genes during limb development remain incompletely delineated. RESULTS: In this report we examine limb development in mice with frameshift mutations in six Hox genes, Hoxa9,10,11 and Hoxd9,10,11. The mice were made with a novel recombineering method that allows the simultaneous targeting of frameshift mutations into multiple flanking genes. The Hoxa9,10,11 (-/-) /Hoxd9,10,11 (-/-) mutant mice show a reduced ulna and radius that is more severe than seen in Hoxa11 (-/-)/Hoxd11 (-/-) mice, indicating a minor role for the flanking Hox9,10 genes in zeugopod development, as well as their primary function in stylopod development. The mutant mice also show severe reduction of Shh expression in the zone of polarizing activity, and decreased Fgf8 expression in the apical ectodermal ridge, thereby better defining the roles of these specific Hox genes in the regulation of critical signaling centers during limb development. Importantly, we also used laser capture microdissection coupled with RNA-Seq to characterize the gene expression programs in wild type and mutant limbs. Resting, proliferative and hypertrophic compartments of E15.5 forelimb zeugopods were examined. The results provide an RNA-Seq characterization of the progression of gene expression patterns during normal endochondral bone formation. In addition of the Hox mutants showed strongly altered expression of Pknox2, Zfp467, Gdf5, Bmpr1b, Dkk3, Igf1, Hand2, Shox2, Runx3, Bmp7 and Lef1, all of which have been previously shown to play important roles in bone formation. CONCLUSIONS: The recombineering based frameshift mutation of the six flanking and paralogous Hoxa9,10,11 and Hoxd9,10,11 genes provides a resource for the analysis of their overlapping functions. Analysis of the Hoxa9,10,11 (-/-) /Hoxd9,10,11 (-/-) mutant limbs confirms and extends the results of previous studies using mice with Hox mutations in single paralogous groups or with entire Hox cluster deletions. The RNA-Seq analysis of specific compartments of the normal and mutant limbs defines the multiple key perturbed pathways downstream of these Hox genes.

Tissue-resident macrophages specifically express Lactotransferrin and Vegfc during ear pinna regeneration in spiny mice.

The details of how macrophages control different healing trajectories (regeneration vs. scar formation) remain poorly defined. Spiny mice (Acomys spp.) can regenerate external ear pinnae tissue, whereas lab mice (Mus musculus) form scar tissue in response to an identical injury. Here, we used this dual species system to dissect macrophage phenotypes between healing modes. We identified secreted factors from activated Acomys macrophages that induce a pro-regenerative phenotype in fibroblasts from both species. Transcriptional profiling of Acomys macrophages and subsequent in vitro tests identified VEGFC, PDGFA, and Lactotransferrin (LTF) as potential pro-regenerative modulators. Examining macrophages in vivo, we found that Acomys-resident macrophages secreted VEGFC and LTF, whereas Mus macrophages do not. Lastly, we demonstrate the requirement for VEGFC during regeneration and find that interrupting lymphangiogenesis delays blastema and new tissue formation. Together, our results demonstrate that cell-autonomous mechanisms govern how macrophages react to the same stimuli to differentially produce factors that facilitate regeneration.

Single-cell sequencing dissects the transcriptional identity of activated fibroblasts and identifies novel persistent distal tubular injury patterns in kidney fibrosis.

Examining kidney fibrosis is crucial for mechanistic understanding and developing targeted strategies against chronic kidney disease (CKD). Persistent fibroblast activation and tubular epithelial cell (TEC) injury are key CKD contributors. However, cellular and transcriptional landscapes of CKD and specific activated kidney fibroblast clusters remain elusive. Here, we analyzed single cell transcriptomic profiles of two clinically relevant kidney fibrosis models which induced robust kidney parenchymal remodeling. We dissected the molecular and cellular landscapes of kidney stroma and newly identified three distinctive fibroblast clusters with "secretory", "contractile" and "vascular" transcriptional enrichments. Also, both injuries generated failed repair TECs (frTECs) characterized by decline of mature epithelial markers and elevation of stromal and injury markers. Notably, frTECs shared transcriptional identity with distal nephron segments of the embryonic kidney. Moreover, we identified that both models exhibited robust and previously unrecognized distal spatial pattern of TEC injury, outlined by persistent elevation of renal TEC injury markers including Krt8 and Vcam1, while the surviving proximal tubules (PTs) showed restored transcriptional signature. We also found that long-term kidney injuries activated a prominent nephrogenic signature, including Sox4 and Hox gene elevation, which prevailed in the distal tubular segments. Our findings might advance understanding of and targeted intervention in fibrotic kidney disease.

Gucy1α1 specifically marks kidney, heart, lung and liver fibroblasts.

Fibrosis is a common outcome of numerous pathologies, including chronic kidney disease (CKD), a progressive renal function deterioration. Current approaches to target activated fibroblasts, key effector contributors to fibrotic tissue remodeling, lack specificity. Here, we report Gucy1α1 as a specific kidney fibroblast marker. Gucy1α1 levels significantly increased over the course of two clinically relevant murine CKD models and directly correlated with established fibrosis markers. Immunofluorescent (IF) imaging showed that Gucy1α1 comprehensively labelled cortical and medullary quiescent and activated fibroblasts in the control kidney and throughout injury progression, respectively. Unlike traditionally used markers platelet derived growth factor receptor beta (Pdgfrß) and vimentin (Vim), Gucy1α1 did not overlap with off-target populations such as podocytes. Notably, Gucy1α1 labelled kidney fibroblasts in both male and female mice. Furthermore, we observed elevated GUCY1α1 expression in the human fibrotic kidney and lung. Studies in the murine models of cardiac and liver fibrosis revealed Gucy1α1 elevation in activated Pdgfrß-, Vim- and alpha smooth muscle actin (αSma)-expressing fibroblasts paralleling injury progression and resolution. Overall, we demonstrate Gucy1α1 as an exclusive fibroblast marker in both sexes. Due to its multiorgan translational potential, GUCY1α1 might provide a novel promising strategy to specifically target and mechanistically examine fibroblasts.

Single-cell sequencing dissects the transcriptional identity of activated fibroblasts and identifies novel persistent distal tubular injury patterns in kidney fibrosis.

Examining kidney fibrosis is crucial for mechanistic understanding and developing targeted strategies against chronic kidney disease (CKD). Persistent fibroblast activation and tubular epithelial cell (TEC) injury are key CKD contributors. However, cellular and transcriptional landscapes of CKD and specific activated kidney fibroblast clusters remain elusive. Here, we analyzed single cell transcriptomic profiles of two clinically relevant kidney fibrosis models which induced robust kidney parenchymal remodeling. We dissected the molecular and cellular landscapes of kidney stroma and newly identified three distinctive fibroblast clusters with "secretory", "contractile" and "vascular" transcriptional enrichments. Also, both injuries generated failed repair TECs (frTECs) characterized by decline of mature epithelial markers and elevation of stromal and injury markers. Notably, frTECs shared transcriptional identity with distal nephron segments of the embryonic kidney. Moreover, we identified that both models exhibited robust and previously unrecognized distal spatial pattern of TEC injury, outlined by persistent elevation of renal TEC injury markers including Krt8, while the surviving proximal tubules (PTs) showed restored transcriptional signature. Furthermore, we found that long-term kidney injuries activated a prominent nephrogenic signature, including Sox4 and Hox gene elevation, which prevailed in the distal tubular segments. Our findings might advance understanding of and targeted intervention in fibrotic kidney disease.

Single cell transcriptomic analysis of HPV16-infected epithelium identifies a keratinocyte subpopulation implicated in cancer.

Persistent HPV16 infection is a major cause of the global cancer burden. The viral life cycle is dependent on the differentiation program of stratified squamous epithelium, but the landscape of keratinocyte subpopulations which support distinct phases of the viral life cycle has yet to be elucidated. Here, single cell RNA sequencing of HPV16 infected compared to uninfected organoids identifies twelve distinct keratinocyte populations, with a subset mapped to reconstruct their respective 3D geography in stratified squamous epithelium. Instead of conventional terminally differentiated cells, an HPV-reprogrammed keratinocyte subpopulation (HIDDEN cells) forms the surface compartment and requires overexpression of the ELF3/ESE-1 transcription factor. HIDDEN cells are detected throughout stages of human carcinogenesis including primary human cervical intraepithelial neoplasias and HPV positive head and neck cancers, and a possible role in promoting viral carcinogenesis is supported by TCGA analyses. Single cell transcriptome information on HPV-infected versus uninfected epithelium will enable broader studies of the role of individual keratinocyte subpopulations in tumor virus infection and cancer evolution.