Transcriptional signals of transformation in human cancer.

BACKGROUND: As normal cells transform into cancers, their cell state changes, which may drive cancer cells into a stem-like or more primordial, foetal, or embryonic cell state. The transcriptomic profile of this final state may encode information about cancer's origin and how cancers relate to their normal cell counterparts. METHODS: Here, we used single-cell atlases to study cancer transformation in transcriptional terms. We utilised bulk transcriptomes across a wide spectrum of adult and childhood cancers, using a previously established method to interrogate their relationship to normal cell states. We extend and validate these findings using single-cell cancer transcriptomes and organ-specific atlases of colorectal and liver cancer. RESULTS: Our bulk transcriptomic data reveals that adult cancers rarely return to an embryonic state, but that a foetal state is a near-universal feature of childhood cancers. This finding was confirmed with single-cell cancer transcriptomes. CONCLUSIONS: Our findings provide a nuanced picture of transformation in human cancer, indicating cancer-specific rather than universal patterns of transformation pervade adult epithelial cancers.

Targetable NOTCH1 rearrangements in reninoma.

Reninomas are exceedingly rare renin-secreting kidney tumours that derive from juxtaglomerular cells, specialised smooth muscle cells that reside at the vascular inlet of glomeruli. They are the central component of the juxtaglomerular apparatus which controls systemic blood pressure through the secretion of renin. We assess somatic changes in reninoma and find structural variants that generate canonical activating rearrangements of, NOTCH1 whilst removing its negative regulator, NRARP. Accordingly, in single reninoma nuclei we observe excessive renin and NOTCH1 signalling mRNAs, with a concomitant non-excess of NRARP expression. Re-analysis of previously published reninoma bulk transcriptomes further corroborates our observation of dysregulated Notch pathway signalling in reninoma. Our findings reveal NOTCH1 rearrangements in reninoma, therapeutically targetable through existing NOTCH1 inhibitors, and indicate that unscheduled Notch signalling may be a disease-defining feature of reninoma.

An integrated single-cell analysis of human adrenal cortex development.

The adrenal glands synthesize and release essential steroid hormones such as cortisol and aldosterone, but many aspects of human adrenal gland development are not well understood. Here, we combined single-cell and bulk RNA sequencing, spatial transcriptomics, IHC, and micro-focus computed tomography to investigate key aspects of adrenal development in the first 20 weeks of gestation. We demonstrate rapid adrenal growth and vascularization, with more cell division in the outer definitive zone (DZ). Steroidogenic pathways favored androgen synthesis in the central fetal zone, but DZ capacity to synthesize cortisol and aldosterone developed with time. Core transcriptional regulators were identified, with localized expression of HOPX (also known as Hop homeobox/homeobox-only protein) in the DZ. Potential ligand-receptor interactions between mesenchyme and adrenal cortex were seen (e.g., RSPO3/LGR4). Growth-promoting imprinted genes were enriched in the developing cortex (e.g., IGF2, PEG3). These findings reveal aspects of human adrenal development and have clinical implications for understanding primary adrenal insufficiency and related postnatal adrenal disorders, such as adrenal tumor development, steroid disorders, and neonatal stress.

Transcriptional signatures associated with persisting CD19 CAR-T cells in children with leukemia.

In the context of relapsed and refractory childhood pre-B cell acute lymphoblastic leukemia (R/R B-ALL), CD19-targeting chimeric antigen receptor (CAR)-T cells often induce durable remissions, which requires the persistence of CAR-T cells. In this study, we systematically analyzed CD19 CAR-T cells of 10 children with R/R B-ALL enrolled in the CARPALL trial via high-throughput single-cell gene expression and T cell receptor sequencing of infusion products and serial blood and bone marrow samples up to 5 years after infusion. We show that long-lived CAR-T cells developed a CD4/CD8 double-negative phenotype with an exhausted-like memory state and distinct transcriptional signature. This persistence signature was dominant among circulating CAR-T cells in all children with a long-lived treatment response for which sequencing data were sufficient (4/4, 100%). The signature was also present across T cell subsets and clonotypes, indicating that persisting CAR-T cells converge transcriptionally. This persistence signature was also detected in two adult patients with chronic lymphocytic leukemia with decade-long remissions who received a different CD19 CAR-T cell product. Examination of single T cell transcriptomes from a wide range of healthy and diseased tissues across children and adults indicated that the persistence signature may be specific to long-lived CAR-T cells. These findings raise the possibility that a universal transcriptional signature of clinically effective, persistent CD19 CAR-T cells exists.

Spatial and molecular profiling of the mononuclear phagocyte network in classic Hodgkin lymphoma.

Classic Hodgkin lymphoma (cHL) has a rich immune infiltrate, which is an intrinsic component of the neoplastic process. Malignant Hodgkin Reed-Sternberg cells (HRSCs) create an immunosuppressive microenvironment by the expression of regulatory molecules, preventing T-cell activation. It has also been demonstrated that mononuclear phagocytes (MNPs) in the vicinity of HRSCs express similar regulatory mechanisms in parallel, and their presence in tissue is associated with inferior patient outcomes. MNPs in cHL have hitherto been identified by a small number of canonical markers and are usually described as tumor-associated macrophages. The organization of MNP networks and interactions with HRSCs remains unexplored at high resolution. Here, we defined the global immune-cell composition of cHL and nonlymphoma lymph nodes, integrating data across single-cell RNA sequencing, spatial transcriptomics, and multiplexed immunofluorescence. We observed that MNPs comprise multiple subsets of monocytes, macrophages, and dendritic cells (DCs). Classical monocytes, macrophages and conventional DC2s were enriched in the vicinity of HRSCs, but plasmacytoid DCs and activated DCs were excluded. Unexpectedly, cDCs and monocytes expressed immunoregulatory checkpoints PD-L1, TIM-3, and the tryptophan-catabolizing protein IDO, at the same level as macrophages. Expression of these molecules increased with age. We also found that classical monocytes are important signaling hubs, potentially controlling the retention of cDC2 and ThExh via CCR1-, CCR4-, CCR5-, and CXCR3-dependent signaling. Enrichment of the cDC2-monocyte-macrophage network in diagnostic biopsies is associated with early treatment failure. These results reveal unanticipated complexity and spatial polarization within the MNP compartment, further demonstrating their potential roles in immune evasion by cHL.

The pathophysiology of bilateral and multifocal Wilms tumors: What we can learn from the study of predisposition syndromes.

Approximately 5% of patients with Wilms tumor present with synchronous bilateral disease. The development of synchronous bilateral Wilms tumor (BWT) is highly suggestive of a genetic or epigenetic predisposition. Patients with known germline predisposition to Wilms tumor (WT1 variants, Beckwith Wiedemann spectrum, TRIM28 variants) have a higher incidence of BWT. This Children's Oncology Group (COG)-International Society for Pediatric Oncology (SIOP-) HARMONICA initiative review for pediatric renal tumors details germline genetic and epigenetic predisposition to BWT development, with an emphasis on alterations in 11p15.5 (ICR1 gain of methylation, paternal uniparental disomy, and postzygotic somatic mosaicism), WT1, TRIM28, and REST. Molecular mechanisms that result in BWT are often also present in multifocal Wilms tumor (multiple separate tumors in one or both kidneys). We identify priority areas for international collaborative research to better understand how predisposing genetic or epigenetic factors associate with response to neoadjuvant chemotherapy, oncologic outcomes, and long-term renal function outcomes.

Mapping single-cell transcriptomes in the intra-tumoral and associated territories of kidney cancer.

Tumor behavior is intricately dependent on the oncogenic properties of cancer cells and their multi-cellular interactions. To understand these dependencies within the wider microenvironment, we studied over 270,000 single-cell transcriptomes and 100 microdissected whole exomes from 12 patients with kidney tumors, prior to validation using spatial transcriptomics. Tissues were sampled from multiple regions of the tumor core, the tumor-normal interface, normal surrounding tissues, and peripheral blood. We find that the tissue-type location of CD8+ T cell clonotypes largely defines their exhaustion state with intra-tumoral spatial heterogeneity that is not well explained by somatic heterogeneity. De novo mutation calling from single-cell RNA-sequencing data allows us to broadly infer the clonality of stromal cells and lineage-trace myeloid cell development. We report six conserved meta-programs that distinguish tumor cell function, and find an epithelial-mesenchymal transition meta-program highly enriched at the tumor-normal interface that co-localizes with IL1B-expressing macrophages, offering a potential therapeutic target.

Precise identification of cancer cells from allelic imbalances in single cell transcriptomes.

A fundamental step of tumour single cell mRNA analysis is separating cancer and non-cancer cells. We show that the common approach to separation, using shifts in average expression, can lead to erroneous biological conclusions. By contrast, allelic imbalances representing copy number changes directly detect the cancer genotype and accurately separate cancer from non-cancer cells. Our findings provide a definitive approach to identifying cancer cells from single cell mRNA sequencing data.

Clonal diversification and histogenesis of malignant germ cell tumours.

Germ cell tumours (GCTs) are a collection of benign and malignant neoplasms derived from primordial germ cells. They are uniquely able to recapitulate embryonic and extraembryonic tissues, which carries prognostic and therapeutic significance. The developmental pathways underpinning GCT initiation and histogenesis are incompletely understood. Here, we study the relationship of histogenesis and clonal diversification in GCTs by analysing the genomes and transcriptomes of 547 microdissected histological units. We find no correlation between genomic and histological heterogeneity. However, we identify unifying features including the retention of fetal developmental transcripts across tissues, expression changes on chromosome 12p, and a conserved somatic evolutionary sequence of whole genome duplication followed by clonal diversification. While this pattern is preserved across all GCTs, the developmental timing of the duplication varies between prepubertal and postpubertal cases. In addition, tumours of younger children exhibit distinct substitution signatures which may lend themselves as potential biomarkers for risk stratification. Our findings portray the extensive diversification of GCT tissues and genetic subclones as randomly distributed, while identifying overarching transcriptional and genomic features.

Single-cell transcriptomics reveals a distinct developmental state of KMT2A-rearranged infant B-cell acute lymphoblastic leukemia.

KMT2A-rearranged infant ALL is an aggressive childhood leukemia with poor prognosis. Here, we investigated the developmental state of KMT2A-rearranged infant B-cell acute lymphoblastic leukemia (B-ALL) using bulk messenger RNA (mRNA) meta-analysis and examination of single lymphoblast transcriptomes against a developing bone marrow reference. KMT2A-rearranged infant B-ALL was uniquely dominated by an early lymphocyte precursor (ELP) state, whereas less adverse NUTM1-rearranged infant ALL demonstrated signals of later developing B cells, in line with most other childhood B-ALLs. We compared infant lymphoblasts with ELP cells and revealed that the cancer harbored hybrid myeloid-lymphoid features, including nonphysiological antigen combinations potentially targetable to achieve cancer specificity. We validated surface coexpression of exemplar combinations by flow cytometry. Through analysis of shared mutations in separate leukemias from a child with infant KMT2A-rearranged B-ALL relapsing as AML, we established that KMT2A rearrangement occurred in very early development, before hematopoietic specification, emphasizing that cell of origin cannot be inferred from the transcriptional state.

Single cell derived mRNA signals across human kidney tumors.

Tumor cells may share some patterns of gene expression with their cell of origin, providing clues into the differentiation state and origin of cancer. Here, we study the differentiation state and cellular origin of 1300 childhood and adult kidney tumors. Using single cell mRNA reference maps of normal tissues, we quantify reference "cellular signals" in each tumor. Quantifying global differentiation, we find that childhood tumors exhibit fetal cellular signals, replacing the presumption of "fetalness" with a quantitative measure of immaturity. By contrast, in adult cancers our assessment refutes the suggestion of dedifferentiation towards a fetal state in most cases. We find an intimate connection between developmental mesenchymal populations and childhood renal tumors. We demonstrate the diagnostic potential of our approach with a case study of a cryptic renal tumor. Our findings provide a cellular definition of human renal tumors through an approach that is broadly applicable to human cancer.

Inherent mosaicism and extensive mutation of human placentas.

Placentas can exhibit chromosomal aberrations that are absent from the fetus1. The basis of this genetic segregation, which is known as confined placental mosaicism, remains unknown. Here we investigated the phylogeny of human placental cells as reconstructed from somatic mutations, using whole-genome sequencing of 86 bulk placental samples (with a median weight of 28 mg) and of 106 microdissections of placental tissue. We found that every bulk placental sample represents a clonal expansion that is genetically distinct, and exhibits a genomic landscape akin to that of childhood cancer in terms of mutation burden and mutational imprints. To our knowledge, unlike any other healthy human tissue studied so far, the placental genomes often contained changes in copy number. We reconstructed phylogenetic relationships between tissues from the same pregnancy, which revealed that developmental bottlenecks genetically isolate placental tissues by separating trophectodermal lineages from lineages derived from the inner cell mass. Notably, there were some cases with full segregation-within a few cell divisions of the zygote-of placental lineages and lineages derived from the inner cell mass. Such early embryonic bottlenecks may enable the normalization of zygotic aneuploidy. We observed direct evidence for this in a case of mosaic trisomic rescue. Our findings reveal extensive mutagenesis in placental tissues and suggest that mosaicism is a typical feature of placental development.

Somatic mutations and single-cell transcriptomes reveal the root of malignant rhabdoid tumours.

Malignant rhabdoid tumour (MRT) is an often lethal childhood cancer that, like many paediatric tumours, is thought to arise from aberrant fetal development. The embryonic root and differentiation pathways underpinning MRT are not firmly established. Here, we study the origin of MRT by combining phylogenetic analyses and single-cell mRNA studies in patient-derived organoids. Comparison of somatic mutations shared between cancer and surrounding normal tissues places MRT in a lineage with neural crest-derived Schwann cells. Single-cell mRNA readouts of MRT differentiation, which we examine by reverting the genetic driver mutation underpinning MRT, SMARCB1 loss, suggest that cells are blocked en route to differentiating into mesenchyme. Quantitative transcriptional predictions indicate that combined HDAC and mTOR inhibition mimic MRT differentiation, which we confirm experimentally. Our study defines the developmental block of MRT and reveals potential differentiation therapies.

Tumor to normal single-cell mRNA comparisons reveal a pan-neuroblastoma cancer cell.

Neuroblastoma is a childhood cancer that resembles developmental stages of the neural crest. It is not established what developmental processes neuroblastoma cancer cells represent. Here, we sought to reveal the phenotype of neuroblastoma cancer cells by comparing cancer (n = 19,723) with normal fetal adrenal single-cell transcriptomes (n = 57,972). Our principal finding was that the neuroblastoma cancer cell resembled fetal sympathoblasts, but no other fetal adrenal cell type. The sympathoblastic state was a universal feature of neuroblastoma cells, transcending cell cluster diversity, individual patients, and clinical phenotypes. We substantiated our findings in 650 neuroblastoma bulk transcriptomes and by integrating canonical features of the neuroblastoma genome with transcriptional signals. Overall, our observations indicate that a pan-neuroblastoma cancer cell state exists, which may be attractive for novel immunotherapeutic and targeted avenues.

Single-cell sequencing reveals clonal expansions of pro-inflammatory synovial CD8 T cells expressing tissue-homing receptors in psoriatic arthritis.

Psoriatic arthritis (PsA) is a debilitating immune-mediated inflammatory arthritis of unknown pathogenesis commonly affecting patients with skin psoriasis. Here we use complementary single-cell approaches to study leukocytes from PsA joints. Mass cytometry demonstrates a 3-fold expansion of memory CD8 T cells in the joints of PsA patients compared to peripheral blood. Meanwhile, droplet-based and plate-based single-cell RNA sequencing of paired T cell receptor alpha and beta chain sequences show pronounced CD8 T cell clonal expansions within the joints. Transcriptome analyses find these expanded synovial CD8 T cells to express cycling, activation, tissue-homing and tissue residency markers. T cell receptor sequence comparison between patients identifies clonal convergence. Finally, chemokine receptor CXCR3 is upregulated in the expanded synovial CD8 T cells, while two CXCR3 ligands, CXCL9 and CXCL10, are elevated in PsA synovial fluid. Our data thus provide a quantitative molecular insight into the cellular immune landscape of psoriatic arthritis.

Embryonal precursors of Wilms tumor.

Adult cancers often arise from premalignant clonal expansions. Whether the same is true of childhood tumors has been unclear. To investigate whether Wilms tumor (nephroblastoma; a childhood kidney cancer) develops from a premalignant background, we examined the phylogenetic relationship between tumors and corresponding normal tissues. In 14 of 23 cases studied (61%), we found premalignant clonal expansions in morphologically normal kidney tissues that preceded tumor development. These clonal expansions were defined by somatic mutations shared between tumor and normal tissues but absent from blood cells. We also found hypermethylation of the H19 locus, a known driver of Wilms tumor development, in 58% of the expansions. Phylogenetic analyses of bilateral tumors indicated that clonal expansions can evolve before the divergence of left and right kidney primordia. These findings reveal embryonal precursors from which unilateral and multifocal cancers develop.

Spatiotemporal immune zonation of the human kidney.

Tissue-resident immune cells are important for organ homeostasis and defense. The epithelium may contribute to these functions directly or by cross-talk with immune cells. We used single-cell RNA sequencing to resolve the spatiotemporal immune topology of the human kidney. We reveal anatomically defined expression patterns of immune genes within the epithelial compartment, with antimicrobial peptide transcripts evident in pelvic epithelium in the mature, but not fetal, kidney. A network of tissue-resident myeloid and lymphoid immune cells was evident in both fetal and mature kidney, with postnatal acquisition of transcriptional programs that promote infection-defense capabilities. Epithelial-immune cross-talk orchestrated localization of antibacterial macrophages and neutrophils to the regions of the kidney most susceptible to infection. Overall, our study provides a global overview of how the immune landscape of the human kidney is zonated to counter the dominant immunological challenge.

Rearrangement bursts generate canonical gene fusions in bone and soft tissue tumors.

Sarcomas are cancers of the bone and soft tissue often defined by gene fusions. Ewing sarcoma involves fusions between EWSR1, a gene encoding an RNA binding protein, and E26 transformation-specific (ETS) transcription factors. We explored how and when EWSR1-ETS fusions arise by studying the whole genomes of Ewing sarcomas. In 52 of 124 (42%) of tumors, the fusion gene arises by a sudden burst of complex, loop-like rearrangements, a process called chromoplexy, rather than by simple reciprocal translocations. These loops always contained the disease-defining fusion at the center, but they disrupted multiple additional genes. The loops occurred preferentially in early replicating and transcriptionally active genomic regions. Similar loops forming canonical fusions were found in three other sarcoma types. Chromoplexy-generated fusions appear to be associated with an aggressive form of Ewing sarcoma. These loops arise early, giving rise to both primary and relapse Ewing sarcoma tumors, which can continue to evolve in parallel.