Pervasive chromosomal instability and karyotype order in tumour evolution.

Chromosomal instability in cancer consists of dynamic changes to the number and structure of chromosomes1,2. The resulting diversity in somatic copy number alterations (SCNAs) may provide the variation necessary for tumour evolution1,3,4. Here we use multi-sample phasing and SCNA analysis of 1,421 samples from 394 tumours across 22 tumour types to show that continuous chromosomal instability results in pervasive SCNA heterogeneity. Parallel evolutionary events, which cause disruption in the same genes (such as BCL9, MCL1, ARNT (also known as HIF1B), TERT and MYC) within separate subclones, were present in 37% of tumours. Most recurrent losses probably occurred before whole-genome doubling, that was found as a clonal event in 49% of tumours. However, loss of heterozygosity at the human leukocyte antigen (HLA) locus and loss of chromosome 8p to a single haploid copy recurred at substantial subclonal frequencies, even in tumours with whole-genome doubling, indicating ongoing karyotype remodelling. Focal amplifications that affected chromosomes 1q21 (which encompasses BCL9, MCL1 and ARNT), 5p15.33 (TERT), 11q13.3 (CCND1), 19q12 (CCNE1) and 8q24.1 (MYC) were frequently subclonal yet appeared to be clonal within single samples. Analysis of an independent series of 1,024 metastatic samples revealed that 13 focal SCNAs were enriched in metastatic samples, including gains in chromosome 8q24.1 (encompassing MYC) in clear cell renal cell carcinoma and chromosome 11q13.3 (encompassing CCND1) in HER2+ breast cancer. Chromosomal instability may enable the continuous selection of SCNAs, which are established as ordered events that often occur in parallel, throughout tumour evolution.

Skin in the game: a review of single-cell and spatial transcriptomics in dermatological research.

Single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics (ST) are two emerging research technologies that uniquely characterize gene expression microenvironments on a cellular or subcellular level. The skin, a clinically accessible tissue composed of diverse, essential cell populations, serves as an ideal target for these high-resolution investigative approaches. Using these tools, researchers are assembling a compendium of data and discoveries in healthy skin as well as a range of dermatologic pathophysiologies, including atopic dermatitis, psoriasis, and cutaneous malignancies. The ongoing advancement of single-cell approaches, coupled with anticipated decreases in cost with increased adoption, will reshape dermatologic research, profoundly influencing disease characterization, prognosis, and ultimately clinical practice.

Spectral clustering of single-cell multi-omics data on multilayer graphs.

MOTIVATION: Single-cell sequencing technologies that simultaneously generate multimodal cellular profiles present opportunities for improved understanding of cell heterogeneity in tissues. How the multimodal information can be integrated to obtain a common cell type identification, however, poses a computational challenge. Multilayer graphs provide a natural representation of multi-omic single-cell sequencing datasets, and finding cell clusters may be understood as a multilayer graph partition problem. RESULTS: We introduce two spectral algorithms on multilayer graphs, spectral clustering on multilayer graphs and the weighted locally linear (WLL) method, to cluster cells in multi-omic single-cell sequencing datasets. We connect these algorithms through a unifying mathematical framework that represents each layer using a Hamiltonian operator and a mixture of its eigenstates to integrate the multiple graph layers, demonstrating in the process that the WLL method is a rigorous multilayer spectral graph theoretic reformulation of the popular Seurat weighted nearest neighbor (WNN) algorithm. Implementing our algorithms and applying them to a CITE-seq dataset of cord blood mononuclear cells yields results similar to the Seurat WNN analysis. Our work thus extends spectral methods to multimodal single-cell data analysis. AVAILABILITY AND IMPLEMENTATION: The code used in this study can be found at https://github.com/jssong-lab/sc-spectrum. All public data used in the article are accurately cited and described in Materials and Methods and in Supplementary Information. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Multiregion exome sequencing of ovarian immature teratomas reveals 2N near-diploid genomes, paucity of somatic mutations, and extensive allelic imbalances shared across mature, immature, and disseminated components.

Immature teratoma is a subtype of malignant germ cell tumor of the ovary that occurs most commonly in the first three decades of life, frequently with bilateral ovarian disease. Despite being the second most common malignant germ cell tumor of the ovary, little is known about its genetic underpinnings. Here we performed multiregion whole-exome sequencing to interrogate the genetic zygosity, clonal relationship, DNA copy number, and mutational status of 52 pathologically distinct tumor components from ten females with ovarian immature teratomas, with bilateral tumors present in five cases and peritoneal dissemination in seven cases. We found that ovarian immature teratomas are genetically characterized by 2N near-diploid genomes with extensive loss of heterozygosity and an absence of genes harboring recurrent somatic mutations or known oncogenic variants. All components within a single ovarian tumor (immature teratoma, mature teratoma with different histologic patterns of differentiation, and yolk sac tumor) were found to harbor an identical pattern of loss of heterozygosity across the genome, indicating a shared clonal origin. In contrast, the four analyzed bilateral teratomas showed distinct patterns of zygosity changes in the right versus left sided tumors, indicating independent clonal origins. All disseminated teratoma components within the peritoneum (including gliomatosis peritonei) shared a clonal pattern of loss of heterozygosity with either the right or left primary ovarian tumor. The observed genomic loss of heterozygosity patterns indicate that diverse meiotic errors contribute to the formation of ovarian immature teratomas, with 11 out of the 15 genetically distinct clones determined to result from nondisjunction errors during meiosis I or II. Overall, these findings suggest that copy-neutral loss of heterozygosity resulting from meiotic abnormalities may be sufficient to generate ovarian immature teratomas from germ cells.

Loss of ZNF750 in ocular and cutaneous sebaceous carcinoma.

BACKGROUND: Sebaceous carcinoma (SeC) is an uncommon malignancy arising from sebaceous glands of the conjunctiva and skin. Recurrent mutations in the ZNF750 were recently identified in ocular SeC. We assessed whether ZNF750 loss is a specific feature of ocular SeC or a general feature of sebaceous tumors. METHODS: Immunostaining for ZNF750 expression was performed in 54 benign and malignant sebocytic proliferations. Staining for ZNF750 was scored on a three-tier scale: positive (>75%), partially positive (5%-74%), and negative (<5%). RESULTS: ZNF750 expression was negative in 4/11 ocular SeC, and partially positive in 4/11 ocular SeC and 6/13 cutaneous SeC. No extraocular tumors were negative. No loss was found in sebaceous adenoma or sebaceous hyperplasia. In nine previously sequenced ocular SeCs, two lacked detectable somatic mutations in ZNF750, but showed complete loss of staining, indicating non-mutational inactivation of ZNF750. CONCLUSION: We show complete loss of the ZNF750 epidermal differentiation regulator in about half of ocular SeC, highlighting the most common genetic defect in this cancer type. Loss of ZNF750 expression is seen even in tumors without truncating mutations and reduced in many of the remaining ocular and cutaneous SeC. In contrast, no ZNF750 loss was detected in benign sebaceous proliferations.

Phylogenetic analyses of melanoma reveal complex patterns of metastatic dissemination.

Melanoma is difficult to treat once it becomes metastatic. However, the precise ancestral relationship between primary tumors and their metastases is not well understood. We performed whole-exome sequencing of primary melanomas and multiple matched metastases from eight patients to elucidate their phylogenetic relationships. In six of eight patients, we found that genetically distinct cell populations in the primary tumor metastasized in parallel to different anatomic sites, rather than sequentially from one site to the next. In five of these six patients, the metastasizing cells had themselves arisen from a common parental subpopulation in the primary, indicating that the ability to establish metastases is a late-evolving trait. Interestingly, we discovered that individual metastases were sometimes founded by multiple cell populations of the primary that were genetically distinct. Such establishment of metastases by multiple tumor subpopulations could help explain why identical resistance variants are identified in different sites after initial response to systemic therapy. One primary tumor harbored two subclones with different oncogenic mutations in CTNNB1, which were both propagated to the same metastasis, raising the possibility that activation of wingless-type mouse mammary tumor virus integration site (WNT) signaling may be involved, as has been suggested by experimental models.

Questioning the Clinical Utility of Exome Sequencing in Developing Countries.

The availability of whole-exome sequencing has revolutionized the study of genetic disease in recent years, particularly in dermatology, where clinical phenotypes are readily recognized. As this technology becomes increasingly affordable and accessible, questions are emerging regarding the clinical and ethical responsibilities of physicians who determine variants underlying disease, especially with regard to children, for whom treatment may be warranted and clinical course improved based on a known genotype. These responsibilities are accentuated in the developing countries, which harbor most consanguineous populations and thus bear the brunt of monogenic genodermatoses. Although many genetic disorders are identified in these populations, limited educational and clinical infrastructure rarely offers opportunities to improve the course of disease. Here we report a genetic study that illustrates these challenges.

The double-edged sword of Notch signaling in cancer.

Recent deep sequencing of cancer genomes has produced an explosion of new data implicating Notch signaling in several human cancers. Unlike most other pathways, these data indicate that Notch signaling can be either oncogenic or tumor suppressive, depending on the cellular context. In some instances, these relationships were predicted from mouse models or presaged by developmental roles for Notch, but in other cases were unanticipated. This review discusses the pathogenic and translational significance of these new findings.

Methods and challenges in timing chromosomal abnormalities within cancer samples.

MOTIVATION: Tumors acquire many chromosomal amplifications, and those acquired early in the lifespan of the tumor may be not only important for tumor growth but also can be used for diagnostic purposes. Many methods infer the order of the accumulation of abnormalities based on their occurrence in a large cohort of patients. Recently, Durinck et al. (2011) and Greenman et al. (2012) developed methods to order a single tumor's chromosomal amplifications based on the patterns of mutations accumulated within those regions. This method offers an unprecedented opportunity to assess the etiology of a single tumor sample, but has not been widely evaluated. RESULTS: We show that the model for timing chromosomal amplifications is limited in scope, particularly for regions with high levels of amplification. We also show that the estimation of the order of events can be sensitive for events that occur early in the progression of the tumor and that the partial maximum likelihood method of Greenman et al. (2012) can give biased estimates, particularly for moderate read coverage or normal contamination. We propose a maximum-likelihood estimation procedure that fully accounts for sequencing variability and show that it outperforms the partial maximum-likelihood estimation method. We also propose a Bayesian estimation procedure that stabilizes the estimates in certain settings. We implement these methods on a small number of ovarian tumors, and the results suggest possible differences in how the tumors acquired amplifications. AVAILABILITY AND IMPLEMENTATION: We provide implementation of these methods in an R package cancerTiming, which is available from the Comprehensive R Archive Network (CRAN) at http://CRAN.R-project.org/.

Loss-of-function mutations in Notch receptors in cutaneous and lung squamous cell carcinoma.

Squamous cell carcinomas (SCCs) are one of the most frequent forms of human malignancy, but, other than TP53 mutations, few causative somatic aberrations have been identified. We identified NOTCH1 or NOTCH2 mutations in ~75% of cutaneous SCCs and in a lesser fraction of lung SCCs, defining a spectrum for the most prevalent tumor suppressor specific to these epithelial malignancies. Notch receptors normally transduce signals in response to ligands on neighboring cells, regulating metazoan lineage selection and developmental patterning. Our findings therefore illustrate a central role for disruption of microenvironmental communication in cancer progression. NOTCH aberrations include frameshift and nonsense mutations leading to receptor truncations as well as point substitutions in key functional domains that abrogate signaling in cell-based assays. Oncogenic gain-of-function mutations in NOTCH1 commonly occur in human T-cell lymphoblastic leukemia/lymphoma and B-cell chronic lymphocytic leukemia. The bifunctional role of Notch in human cancer thus emphasizes the context dependency of signaling outcomes and suggests that targeted inhibition of the Notch pathway may induce squamous epithelial malignancies.

Short cell cycle duration is a phenotype of human epidermal stem cells.

BACKGROUND: A traditional view is that stem cells (SCs) divide slowly. Meanwhile, both embryonic and pluripotent SCs display a shorter cell cycle duration (CCD) in comparison to more committed progenitors (CPs). METHODS: We examined the in vitro proliferation and cycling behavior of somatic adult human cells using live cell imaging of passage zero keratinocytes and single-cell RNA sequencing. RESULTS: We found two populations of keratinocytes: those with short CCD and protracted near exponential growth, and those with long CCD and terminal differentiation. Applying the ergodic principle, the comparative numbers of cycling cells in S phase in an enriched population of SCs confirmed a shorter CCD than CPs. Further, analysis of single-cell RNA sequencing of cycling adult human keratinocyte SCs and CPs indicated a shortening of both G1 and G2M phases in the SC. CONCLUSIONS: Contrary to the pervasive paradigm, SCs progress through cell cycle more quickly than more differentiated dividing CPs. Thus, somatic human adult keratinocyte SCs may divide infrequently, but divide rapidly when they divide. Additionally, it was found that SC-like proliferation persisted in vitro.

A single-cell atlas of IL-23 inhibition in cutaneous psoriasis distinguishes clinical response.

Psoriasis vulgaris and other chronic inflammatory diseases improve markedly with therapeutic blockade of interleukin-23 (IL-23) signaling, but the genetic mechanisms underlying clinical responses remain poorly understood. Using single-cell transcriptomics, we profiled immune cells isolated from lesional psoriatic skin before and during IL-23 blockade. In clinically responsive patients, a psoriatic transcriptional signature in skin-resident memory T cells was strongly attenuated. In contrast, poorly responsive patients were distinguished by persistent activation of IL-17-producing T (T17) cells, a mechanism distinct from alternative cytokine signaling or resistance isolated to epidermal keratinocytes. Even in IL-23 blockade-responsive patients, we detected a recurring set of recalcitrant, disease-specific transcriptional abnormalities. This irreversible immunological state may necessitate ongoing IL-23 inhibition. Spatial transcriptomic analyses also suggested that successful IL-23 blockade requires dampening of >90% of IL-17-induced response in lymphocyte-adjacent keratinocytes, an unexpectedly high threshold. Collectively, our data establish a patient-level paradigm for dissecting responses to immunomodulatory treatments.

Isolation of human cutaneous immune cells for single-cell RNA sequencing.

Single-cell RNA sequencing (scRNA-seq) allows for high-resolution analysis of transcriptionally dysregulated cell subpopulations in inflammatory diseases. However, it can be challenging to properly isolate viable immune cells from human skin for scRNA-seq due to its barrier properties. Here, we present a protocol to isolate high-viability human cutaneous immune cells. We describe steps for obtaining and enzymatically dissociating a skin biopsy specimen and isolating immune cells using flow cytometry. We then provide an overview of downstream computational techniques to analyze sequencing data. For complete details on the use and execution of this protocol, please refer to Cook et al. (2022)1 and Liu et al. (2022).2.

STmut: a framework for visualizing somatic alterations in spatial transcriptomics data of cancer.

Spatial transcriptomic technologies, such as the Visium platform, measure gene expression in different regions of tissues. Here, we describe new software, STmut, to visualize somatic point mutations, allelic imbalance, and copy number alterations in Visium data. STmut is tested on fresh-frozen Visium data, formalin-fixed paraffin-embedded (FFPE) Visium data, and tumors with and without matching DNA sequencing data. Copy number is inferred on all conditions, but the chemistry of the FFPE platform does not permit analyses of single nucleotide variants. Taken together, we propose solutions to add the genetic dimension to spatial transcriptomic data and describe the limitations of different datatypes.

GLUT3 promotes macrophage signaling and function via RAS-mediated endocytosis in atopic dermatitis and wound healing.

The facilitative GLUT1 and GLUT3 hexose transporters are expressed abundantly in macrophages, but whether they have distinct functions remains unclear. We confirmed that GLUT1 expression increased after M1 polarization stimuli and found that GLUT3 expression increased after M2 stimulation in macrophages. Conditional deletion of Glut3 (LysM-Cre Glut3fl/fl) impaired M2 polarization of bone marrow-derived macrophages. Alternatively activated macrophages from the skin of patients with atopic dermatitis showed increased GLUT3 expression, and a calcipotriol-induced model of atopic dermatitis was rescued in LysM-Cre Glut3fl/fl mice. M2-like macrophages expressed GLUT3 in human wound tissues as assessed by transcriptomics and costaining, and GLUT3 expression was significantly decreased in nonhealing, compared with healing, diabetic foot ulcers. In an excisional wound healing model, LysM-Cre Glut3fl/fl mice showed significantly impaired M2 macrophage polarization and delayed wound healing. GLUT3 promoted IL-4/STAT6 signaling, independently of its glucose transport activity. Unlike plasma membrane-localized GLUT1, GLUT3 was localized primarily to endosomes and was required for the efficient endocytosis of IL-4Rα subunits. GLUT3 interacted directly with GTP-bound RAS in vitro and in vivo through its intracytoplasmic loop domain, and this interaction was required for efficient STAT6 activation and M2 polarization. PAK activation and macropinocytosis were also impaired without GLUT3, suggesting broader roles for GLUT3 in the regulation of endocytosis. Thus, GLUT3 is required for efficient alternative macrophage polarization and function, through a glucose transport-independent, RAS-mediated role in the regulation of endocytosis and IL-4/STAT6 activation.

Use of an expanded immunohistochemical panel to distinguish cutaneous Hodgkin lymphoma from histopathologic imitators.

In lymph nodes, classical Hodgkin lymphoma can typically be distinguished from non-Hodgkin lymphoma (NHL) by the presence of Hodgkin and Reed-Sternberg cells that co-express CD30 and CD15. However, anaplastic large cell lymphoma (ALCL) and diffuse large B-cell lymphoma (DLBCL) can show identical features, and some cases of classical Hodgkin lymphoma lack CD15 expression, rendering them difficult to differentiate from CD30-positive NHL. The differential diagnosis of cutaneous Hodgkin lymphoma similarly includes ALCL and DLBCL, and, additionally, tumors of mycosis fungoides. Recent studies have shown that classical Hodgkin lymphoma is of B-cell origin in virtually all cases, and shows at least focal weak expression of the B-cell marker PAX5 and often focal weak expression and no expression of the B-cell markers Oct-2 and BOB.1, respectively. All three of these markers are almost invariably absent in T-cell lymphomas and are strongly expressed in B-cell lymphomas. We report a 40-year-old man with classical Hodgkin lymphoma who developed cutaneous nodules. A biopsy from one revealed Hodgkin/Reed-Sternberg cells with a similar immunophenotype to the diagnostic lymph node biopsy, namely CD30+/CD15+, diffusely but weakly PAX5+, focally weakly Oct-2+ and lacking BOB.1 expression, thereby confirming a diagnosis of cutaneous Hodgkin lymphoma. To our knowledge, this is the first report of the expression pattern of the combination of PAX5, Oct-2 and BOB.1 in the context of cutaneous involvement by Hodgkin lymphoma.

Defining Patient-Level Molecular Heterogeneity in Psoriasis Vulgaris Based on Single-Cell Transcriptomics.

Identifying genetic variation underlying human diseases establishes targets for therapeutic development and helps tailor treatments to individual patients. Large-scale transcriptomic profiling has extended the study of such molecular heterogeneity between patients to somatic tissues. However, the lower resolution of bulk RNA profiling, especially in a complex, composite tissue such as the skin, has limited its success. Here we demonstrate approaches to interrogate patient-level molecular variance in a chronic skin inflammatory disease, psoriasis vulgaris, leveraging single-cell RNA-sequencing of CD45+ cells isolated from active lesions. Highly psoriasis-specific transcriptional abnormalities display greater than average inter-individual variance, nominating them as potential sources of clinical heterogeneity. We find that one of these chemokines, CXCL13, demonstrates significant correlation with severity of lesions within our patient series. Our analyses also establish that genes elevated in psoriatic skin-resident memory T cells are enriched for programs orchestrating chromatin and CDC42-dependent cytoskeleton remodeling, specific components of which are distinctly correlated with and against Th17 identity on a single-cell level. Collectively, these analyses describe systematic means to dissect cell type- and patient-level differences in cutaneous psoriasis using high-resolution transcriptional profiles of human inflammatory disease.

Classification of human chronic inflammatory skin disease based on single-cell immune profiling.

Inflammatory conditions represent the largest class of chronic skin disease, but the molecular dysregulation underlying many individual cases remains unclear. Single-cell RNA sequencing (scRNA-seq) has increased precision in dissecting the complex mixture of immune and stromal cell perturbations in inflammatory skin disease states. We single-cell-profiled CD45+ immune cell transcriptomes from skin samples of 31 patients (7 atopic dermatitis, 8 psoriasis vulgaris, 2 lichen planus (LP), 1 bullous pemphigoid (BP), 6 clinical/histopathologically indeterminate rashes, and 7 healthy controls). Our data revealed active proliferative expansion of the Treg and Trm components and universal T cell exhaustion in human rashes, with a relative attenuation of antigen-presenting cells. Skin-resident memory T cells showed the greatest transcriptional dysregulation in both atopic dermatitis and psoriasis, whereas atopic dermatitis also demonstrated recurrent abnormalities in ILC and CD8+ cytotoxic lymphocytes. Transcript signatures differentiating these rash types included genes previously implicated in T helper cell (TH2)/TH17 diatheses, segregated in unbiased functional networks, and accurately identified disease class in untrained validation data sets. These gene signatures were able to classify clinicopathologically ambiguous rashes with diagnoses consistent with therapeutic response. Thus, we have defined major classes of human inflammatory skin disease at the molecular level and described a quantitative method to classify indeterminate instances of pathologic inflammation. To make this approach accessible to the scientific community, we created a proof-of-principle web interface (RashX), where scientists and clinicians can visualize their patient-level rash scRNA-seq-derived data in the context of our TH2/TH17 transcriptional framework.

A single-cell transcriptional gradient in human cutaneous memory T cells restricts Th17/Tc17 identity.

The homeostatic mechanisms that fail to restrain chronic tissue inflammation in diseases, such as psoriasis vulgaris, remain incompletely understood. We profiled transcriptomes and epitopes of single psoriatic and normal skin-resident T cells, revealing a gradated transcriptional program of coordinately regulated inflammation-suppressive genes. This program, which is sharply suppressed in lesional skin, strikingly restricts Th17/Tc17 cytokine and other inflammatory mediators on the single-cell level. CRISPR-based deactivation of two core components of this inflammation-suppressive program, ZFP36L2 and ZFP36, replicates the interleukin-17A (IL-17A), granulocyte macrophage-colony-stimulating factor (GM-CSF), and interferon gamma (IFNγ) elevation in psoriatic memory T cells deficient in these transcripts, functionally validating their influence. Combinatoric expression analysis indicates the suppression of specific inflammatory mediators by individual program members. Finally, we find that therapeutic IL-23 blockade reduces Th17/Tc17 cell frequency in lesional skin but fails to normalize this inflammatory-suppressive program, suggesting how treated lesions may be primed for recurrence after withdrawal of treatment.