Differentiation-related epigenomic changes define clinically distinct keratinocyte cancer subclasses.

Keratinocyte cancers (KC) are the most prevalent malignancies in fair-skinned populations, posing a significant medical and economic burden to health systems. KC originate in the epidermis and mainly comprise basal cell carcinoma (BCC) and cutaneous squamous cell carcinoma (cSCC). Here, we combined single-cell multi-omics, transcriptomics, and methylomics to investigate the epigenomic dynamics during epidermal differentiation. We identified ~3,800 differentially accessible regions between undifferentiated and differentiated keratinocytes, corresponding to regulatory regions associated with key transcription factors. DNA methylation at these regions defined AK/cSCC subtypes with epidermal stem cell- or keratinocyte-like features. Using cell-type deconvolution tools and integration of bulk and single-cell methylomes, we demonstrate that these subclasses are consistent with distinct cells-of-origin. Further characterization of the phenotypic traits of the subclasses and the study of additional unstratified KC entities uncovered distinct clinical features for the subclasses, linking invasive and metastatic KC cases with undifferentiated cells-of-origin. Our study provides a thorough characterization of the epigenomic dynamics underlying human keratinocyte differentiation and uncovers novel links between KC cells-of-origin and their prognosis.

Activation of gp130 signaling in T cells drives TH17-mediated multi-organ autoimmunity.

The IL-6-gp130-STAT3 signaling axis is a major regulator of inflammation. Activating mutations in the gene encoding gp130 and germline gain-of-function mutations in STAT3 (STAT3GOF) are associated with multi-organ autoimmunity, severe morbidity, and adverse prognosis. To dissect crucial cellular subsets and disease biology involved in activated gp130 signaling, the gp130-JAK-STAT3 axis was constitutively activated using a transgene, L-gp130, specifically targeted to T cells. Activating gp130 signaling in T cells in vivo resulted in fatal, early onset, multi-organ autoimmunity in mice that resembled human STAT3GOF disease. Female mice had more rapid disease progression than male mice. On a cellular level, gp130 signaling induced the activation and effector cell differentiation of T cells, promoted the expansion of T helper type 17 (TH17) cells, and impaired the activity of regulatory T cells. Transcriptomic profiling of CD4+ and CD8+ T cells from these mice revealed commonly dysregulated genes and a gene signature that, when applied to human transcriptomic data, improved the segregation of patients with transcriptionally diverse STAT3GOF mutations from healthy controls. The findings demonstrate that increased gp130-STAT3 signaling leads to TH17-driven autoimmunity that phenotypically resembles human STAT3GOF disease.

Functionally distinct cancer-associated fibroblast subpopulations establish a tumor promoting environment in squamous cell carcinoma.

Cutaneous squamous cell carcinoma (cSCC) is a serious public health problem due to its high incidence and metastatic potential. It may progress from actinic keratosis (AK), a precancerous lesion, or the in situ carcinoma, Bowen's disease (BD). During this progression, malignant keratinocytes activate dermal fibroblasts into tumor promoting cancer-associated fibroblasts (CAFs), whose origin and emergence remain largely unknown. Here, we generate and analyze >115,000 single-cell transcriptomes from healthy skin, BD and cSCC of male donors. Our results reveal immunoregulatory and matrix-remodeling CAF subtypes that may derive from pro-inflammatory and mesenchymal fibroblasts, respectively. These CAF subtypes are largely absent in AK and interact with different cell types to establish a pro-tumorigenic microenvironment. These findings are cSCC-specific and could not be recapitulated in basal cell carcinomas. Our study provides important insights into the potential origin and functionalities of dermal CAFs that will be highly beneficial for the specific targeting of the cSCC microenvironment.

TET2 and TET3 loss disrupts small intestine differentiation and homeostasis.

TET2/3 play a well-known role in epigenetic regulation and mouse development. However, their function in cellular differentiation and tissue homeostasis remains poorly understood. Here we show that ablation of TET2/3 in intestinal epithelial cells results in a murine phenotype characterized by a severe homeostasis imbalance in the small intestine. Tet2/3-deleted mice show a pronounced loss of mature Paneth cells as well as fewer Tuft and more Enteroendocrine cells. Further results show major changes in DNA methylation at putative enhancers, which are associated with cell fate-determining transcription factors and functional effector genes. Notably, pharmacological inhibition of DNA methylation partially rescues the methylation and cellular defects. TET2/3 loss also alters the microbiome, predisposing the intestine to inflammation under homeostatic conditions and acute inflammation-induced death. Together, our results uncover previously unrecognized critical roles for DNA demethylation, possibly occurring subsequently to chromatin opening during intestinal development, culminating in the establishment of normal intestinal crypts.

Resolving the spatial architecture of myeloma and its microenvironment at the single-cell level.

In multiple myeloma spatial differences in the subclonal architecture, molecular signatures and composition of the microenvironment remain poorly characterized. To address this shortcoming, we perform multi-region sequencing on paired random bone marrow and focal lesion samples from 17 newly diagnosed patients. Using single-cell RNA- and ATAC-seq we find a median of 6 tumor subclones per patient and unique subclones in focal lesions. Genetically identical subclones display different levels of spatial transcriptional plasticity, including nearly identical profiles and pronounced heterogeneity at different sites, which can include differential expression of immunotherapy targets, such as CD20 and CD38. Macrophages are significantly depleted in the microenvironment of focal lesions. We observe proportional changes in the T-cell repertoire but no site-specific expansion of T-cell clones in intramedullary lesions. In conclusion, our results demonstrate the relevance of considering spatial heterogeneity in multiple myeloma with potential implications for models of cell-cell interactions and disease progression.

Single-Cell RNA Profiling of Human Skin Reveals Age-Related Loss of Dermal Sheath Cells and Their Contribution to a Juvenile Phenotype.

The dermal sheath (DS) is a population of mesenchyme-derived skin cells with emerging importance for skin homeostasis. The DS includes hair follicle dermal stem cells, which exhibit self-renewal and serve as bipotent progenitors of dermal papilla (DP) cells and DS cells. Upon aging, stem cells exhibit deficiencies in self-renewal and their number is reduced. While the DS of mice has been examined in considerable detail, our knowledge of the human DS, the pathways contributing to its self-renewal and differentiation capacity and potential paracrine effects important for tissue regeneration and aging is very limited. Using single-cell RNA sequencing of human skin biopsies from donors of different ages we have now analyzed the transcriptome of 72,048 cells, including 50,149 fibroblasts. Our results show that DS cells that exhibit stem cell characteristics were lost upon aging. We further show that HES1, COL11A1, MYL4 and CTNNB1 regulate DS stem cell characteristics. Finally, the DS secreted protein Activin A showed paracrine effects on keratinocytes and dermal fibroblasts, promoting proliferation, epidermal thickness and pro-collagen production. Our work provides a detailed description of human DS identity on the single-cell level, its loss upon aging, its stem cell characteristics and its contribution to a juvenile skin phenotype.

Single-cell transcriptomes of the human skin reveal age-related loss of fibroblast priming.

Fibroblasts are an essential cell population for human skin architecture and function. While fibroblast heterogeneity is well established, this phenomenon has not been analyzed systematically yet. We have used single-cell RNA sequencing to analyze the transcriptomes of more than 5,000 fibroblasts from a sun-protected area in healthy human donors. Our results define four main subpopulations that can be spatially localized and show differential secretory, mesenchymal and pro-inflammatory functional annotations. Importantly, we found that this fibroblast 'priming' becomes reduced with age. We also show that aging causes a substantial reduction in the predicted interactions between dermal fibroblasts and other skin cells, including undifferentiated keratinocytes at the dermal-epidermal junction. Our work thus provides evidence for a functional specialization of human dermal fibroblasts and identifies the partial loss of cellular identity as an important age-related change in the human dermis. These findings have important implications for understanding human skin aging and its associated phenotypes.

Loss of G9a preserves mutation patterns but increases chromatin accessibility, genomic instability and aggressiveness in skin tumours.

Mutations in, and the altered expression of, epigenetic modifiers are pervasive in human tumours, making epigenetic factors attractive antitumour targets. The open-versus-closed chromatin state within the cells-of-origin of cancer correlates with the uneven distribution of mutations. However, the long-term effect of targeting epigenetic modifiers on mutability in patients with cancer is unclear. Here, we increased chromatin accessibility by deleting the histone H3 lysine 9 (H3K9) methyltransferase G9a in murine epidermis and show that this does not alter the single nucleotide variant burden or global genomic distribution in chemical mutagen-induced squamous tumours. G9a-depleted tumours develop after a prolonged latency compared with their wild-type counterparts, but are more aggressive and have an expanded cancer progenitor pool, pronounced genomic instability and frequent loss-of-function p53 mutations. Thus, we call for caution when assessing long-term therapeutic benefits of chromatin modifier inhibitors, which may promote more aggressive disease.