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.

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.

Decoding human fetal liver haematopoiesis.

Definitive haematopoiesis in the fetal liver supports self-renewal and differentiation of haematopoietic stem cells and multipotent progenitors (HSC/MPPs) but remains poorly defined in humans. Here, using single-cell transcriptome profiling of approximately 140,000 liver and 74,000 skin, kidney and yolk sac cells, we identify the repertoire of human blood and immune cells during development. We infer differentiation trajectories from HSC/MPPs and evaluate the influence of the tissue microenvironment on blood and immune cell development. We reveal physiological erythropoiesis in fetal skin and the presence of mast cells, natural killer and innate lymphoid cell precursors in the yolk sac. We demonstrate a shift in the haemopoietic composition of fetal liver during gestation away from being predominantly erythroid, accompanied by a parallel change in differentiation potential of HSC/MPPs, which we functionally validate. Our integrated map of fetal liver haematopoiesis provides a blueprint for the study of paediatric blood and immune disorders, and a reference for harnessing the therapeutic potential of HSC/MPPs.

Intra-tumour diversification in colorectal cancer at the single-cell level.

Every cancer originates from a single cell. During expansion of the neoplastic cell population, individual cells acquire genetic and phenotypic differences from each other. Here, to investigate the nature and extent of intra-tumour diversification, we characterized organoids derived from multiple single cells from three colorectal cancers as well as from adjacent normal intestinal crypts. Colorectal cancer cells showed extensive mutational diversification and carried several times more somatic mutations than normal colorectal cells. Most mutations were acquired during the final dominant clonal expansion of the cancer and resulted from mutational processes that are absent from normal colorectal cells. Intra-tumour diversification of DNA methylation and transcriptome states also occurred; these alterations were cell-autonomous, stable, and followed the phylogenetic tree of each cancer. There were marked differences in responses to anticancer drugs between even closely related cells of the same tumour. The results indicate that colorectal cancer cells experience substantial increases in somatic mutation rate compared to normal colorectal cells, and that genetic diversification of each cancer is accompanied by pervasive, stable and inherited differences in the biological states of individual cancer cells.

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.

BBKNN: fast batch alignment of single cell transcriptomes.

MOTIVATION: Increasing numbers of large scale single cell RNA-Seq projects are leading to a data explosion, which can only be fully exploited through data integration. A number of methods have been developed to combine diverse datasets by removing technical batch effects, but most are computationally intensive. To overcome the challenge of enormous datasets, we have developed BBKNN, an extremely fast graph-based data integration algorithm. We illustrate the power of BBKNN on large scale mouse atlasing data, and favourably benchmark its run time against a number of competing methods. AVAILABILITY AND IMPLEMENTATION: BBKNN is available at https://github.com/Teichlab/bbknn, along with documentation and multiple example notebooks, and can be installed from pip. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

ChIP-seq analysis reveals distinct H3K27me3 profiles that correlate with transcriptional activity.

Transcriptional control is dependent on a vast network of epigenetic modifications. One epigenetic mark of particular interest is tri-methylation of lysine 27 on histone H3 (H3K27me3), which is catalysed and maintained by Polycomb Repressive Complex 2 (PRC2). Although this histone mark is studied widely, the precise relationship between its local pattern of enrichment and regulation of gene expression is currently unclear. We have used ChIP-seq to generate genome-wide maps of H3K27me3 enrichment, and have identified three enrichment profiles with distinct regulatory consequences. First, a broad domain of H3K27me3 enrichment across the body of genes corresponds to the canonical view of H3K27me3 as inhibitory to transcription. Second, a peak of enrichment around the transcription start site (TSS) is commonly associated with 'bivalent' genes, where H3K4me3 also marks the TSS. Finally and most surprisingly, we identified an enrichment profile with a peak in the promoter of genes that is associated with active transcription. Genes with each of these three profiles were found in different proportions in each of the cell types studied. The data analysis techniques developed here will be useful for the identification of common enrichment profiles for other histone modifications that have important consequences for transcriptional regulation.

Circulatory HMGB1 is an early predictive and prognostic biomarker of ARDS and mortality in a swine model of polytrauma.

Acute respiratory distress syndrome (ARDS) is a leading cause of morbidity and mortality in polytrauma patients. Pharmacological treatments of ARDS are lacking, and ARDS patients rely on supportive care. Accurate diagnosis of ARDS is vital for early intervention and improved outcomes but is presently delayed up to days. The use of biomarkers for early identification of ARDS development is a potential solution. Inflammatory mediators high-mobility group box 1 (HMGB1), syndecan-1 (SDC-1), and C3a have been previously proposed as potential biomarkers. For this study, we analyzed these biomarkers in animals undergoing smoke inhalation and 40% total body surface area burns, followed by intensive care for 72 h post-injury (PI) to determine their association with ARDS and mortality. We found that the levels of inflammatory mediators in serum were affected, as well as the degree of HMGB1 and Toll-like receptor 4 (TLR4) signal activation in the lung. The results showed significantly increased HMGB1 expression levels in animals that developed ARDS compared with those that did not. Receiver operating characteristic (ROC) analysis showed that HMGB1 levels at 6 h PI were significantly associated with ARDS development (AUROC=0.77) and mortality (AUROC=0.82). Logistic regression analysis revealed that levels of HMGB1 ≥24.10 ng/ml are associated with a 13-fold higher incidence of ARDS [OR:13.57 (2.76-104.3)], whereas the levels of HMGB1 ≥31.39 ng/ml are associated with a 12-fold increase in mortality [OR: 12.00 (2.36-93.47)]. In addition, we found that mesenchymal stem cell (MSC) therapeutic treatment led to a significant decrease in systemic HMGB1 elevation but failed to block SDC-1 and C3a increases. Immunohistochemistry analyses showed that smoke inhalation and burn injury induced the expression of HMGB1 and TLR4 and stimulated co-localization of HMGB1 and TLR4 in the lung. Interestingly, MSC treatment reduced the presence of HMGB1, TLR4, and the HMGB1-TLR4 co-localization. These results show that serum HMGB1 is a prognostic biomarker for predicting the incidence of ARDS and mortality in swine with smoke inhalation and burn injury. Therapeutically blocking HMGB1 signal activation might be an effective approach for attenuating ARDS development in combat casualties or civilian patients.

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.

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.

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.

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.

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.

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.

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.