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.

A spatially resolved atlas of the human lung characterizes a gland-associated immune niche.

Single-cell transcriptomics has allowed unprecedented resolution of cell types/states in the human lung, but their spatial context is less well defined. To (re)define tissue architecture of lung and airways, we profiled five proximal-to-distal locations of healthy human lungs in depth using multi-omic single cell/nuclei and spatial transcriptomics (queryable at lungcellatlas.org ). Using computational data integration and analysis, we extend beyond the suspension cell paradigm and discover macro and micro-anatomical tissue compartments including previously unannotated cell types in the epithelial, vascular, stromal and nerve bundle micro-environments. We identify and implicate peribronchial fibroblasts in lung disease. Importantly, we discover and validate a survival niche for IgA plasma cells in the airway submucosal glands (SMG). We show that gland epithelial cells recruit B cells and IgA plasma cells, and promote longevity and antibody secretion locally through expression of CCL28, APRIL and IL-6. This new 'gland-associated immune niche' has implications for respiratory health.

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 roadmap of human gonadal development.

Gonadal development is a complex process that involves sex determination followed by divergent maturation into either testes or ovaries1. Historically, limited tissue accessibility, a lack of reliable in vitro models and critical differences between humans and mice have hampered our knowledge of human gonadogenesis, despite its importance in gonadal conditions and infertility. Here, we generated a comprehensive map of first- and second-trimester human gonads using a combination of single-cell and spatial transcriptomics, chromatin accessibility assays and fluorescent microscopy. We extracted human-specific regulatory programmes that control the development of germline and somatic cell lineages by profiling equivalent developmental stages in mice. In both species, we define the somatic cell states present at the time of sex specification, including the bipotent early supporting population that, in males, upregulates the testis-determining factor SRY and sPAX8s, a gonadal lineage located at the gonadal-mesonephric interface. In females, we resolve the cellular and molecular events that give rise to the first and second waves of granulosa cells that compartmentalize the developing ovary to modulate germ cell differentiation. In males, we identify human SIGLEC15+ and TREM2+ fetal testicular macrophages, which signal to somatic cells outside and inside the developing testis cords, respectively. This study provides a comprehensive spatiotemporal map of human and mouse gonadal differentiation, which can guide in vitro gonadogenesis.

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.

Local and systemic responses to SARS-CoV-2 infection in children and adults.

It is not fully understood why COVID-19 is typically milder in children1-3. Here, to examine the differences between children and adults in their response to SARS-CoV-2 infection, we analysed paediatric and adult patients with COVID-19 as well as healthy control individuals (total n = 93) using single-cell multi-omic profiling of matched nasal, tracheal, bronchial and blood samples. In the airways of healthy paediatric individuals, we observed cells that were already in an interferon-activated state, which after SARS-CoV-2 infection was further induced especially in airway immune cells. We postulate that higher paediatric innate interferon responses restrict viral replication and disease progression. The systemic response in children was characterized by increases in naive lymphocytes and a depletion of natural killer cells, whereas, in adults, cytotoxic T cells and interferon-stimulated subpopulations were significantly increased. We provide evidence that dendritic cells initiate interferon signalling in early infection, and identify epithelial cell states associated with COVID-19 and age. Our matching nasal and blood data show a strong interferon response in the airways with the induction of systemic interferon-stimulated populations, which were substantially reduced in paediatric patients. Together, we provide several mechanisms that explain the milder clinical syndrome observed in children.

Mapping the temporal and spatial dynamics of the human endometrium in vivo and in vitro.

The endometrium, the mucosal lining of the uterus, undergoes dynamic changes throughout the menstrual cycle in response to ovarian hormones. We have generated dense single-cell and spatial reference maps of the human uterus and three-dimensional endometrial organoid cultures. We dissect the signaling pathways that determine cell fate of the epithelial lineages in the lumenal and glandular microenvironments. Our benchmark of the endometrial organoids reveals the pathways and cell states regulating differentiation of the secretory and ciliated lineages both in vivo and in vitro. In vitro downregulation of WNT or NOTCH pathways increases the differentiation efficiency along the secretory and ciliated lineages, respectively. We utilize our cellular maps to deconvolute bulk data from endometrial cancers and endometriotic lesions, illuminating the cell types dominating in each of these disorders. These mechanistic insights provide a platform for future development of treatments for common conditions including endometriosis and endometrial carcinoma.

Blood and immune development in human fetal bone marrow and Down syndrome.

Haematopoiesis in the bone marrow (BM) maintains blood and immune cell production throughout postnatal life. Haematopoiesis first emerges in human BM at 11-12 weeks after conception1,2, yet almost nothing is known about how fetal BM (FBM) evolves to meet the highly specialized needs of the fetus and newborn. Here we detail the development of FBM, including stroma, using multi-omic assessment of mRNA and multiplexed protein epitope expression. We find that the full blood and immune cell repertoire is established in FBM in a short time window of 6-7 weeks early in the second trimester. FBM promotes rapid and extensive diversification of myeloid cells, with granulocytes, eosinophils and dendritic cell subsets emerging for the first time. The substantial expansion of B lymphocytes in FBM contrasts with fetal liver at the same gestational age. Haematopoietic progenitors from fetal liver, FBM and cord blood exhibit transcriptional and functional differences that contribute to tissue-specific identity and cellular diversification. Endothelial cell types form distinct vascular structures that we show are regionally compartmentalized within FBM. Finally, we reveal selective disruption of B lymphocyte, erythroid and myeloid development owing to a cell-intrinsic differentiation bias as well as extrinsic regulation through an altered microenvironment in Down syndrome (trisomy 21).

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.

Single-cell multi-omics analysis of the immune response in COVID-19.

Analysis of human blood immune cells provides insights into the coordinated response to viral infections such as severe acute respiratory syndrome coronavirus 2, which causes coronavirus disease 2019 (COVID-19). We performed single-cell transcriptome, surface proteome and T and B lymphocyte antigen receptor analyses of over 780,000 peripheral blood mononuclear cells from a cross-sectional cohort of 130 patients with varying severities of COVID-19. We identified expansion of nonclassical monocytes expressing complement transcripts (CD16+C1QA/B/C+) that sequester platelets and were predicted to replenish the alveolar macrophage pool in COVID-19. Early, uncommitted CD34+ hematopoietic stem/progenitor cells were primed toward megakaryopoiesis, accompanied by expanded megakaryocyte-committed progenitors and increased platelet activation. Clonally expanded CD8+ T cells and an increased ratio of CD8+ effector T cells to effector memory T cells characterized severe disease, while circulating follicular helper T cells accompanied mild disease. We observed a relative loss of IgA2 in symptomatic disease despite an overall expansion of plasmablasts and plasma cells. Our study highlights the coordinated immune response that contributes to COVID-19 pathogenesis and reveals discrete cellular components that can be targeted for therapy.

UBE2QL1 is disrupted by a constitutional translocation associated with renal tumor predisposition and is a novel candidate renal tumor suppressor gene.

Investigation of rare familial forms of renal cell carcinoma (RCC) has led to the identification of genes such as VHL and MET that are also implicated in the pathogenesis of sporadic RCC. In order to identify a novel candidate renal tumor suppressor gene, we characterized the breakpoints of a constitutional balanced translocation, t(5;19)(p15.3;q12), associated with familial RCC and found that a previously uncharacterized gene UBE2QL1 was disrupted by the chromosome 5 breakpoint. UBE2QL1 mRNA expression was downregulated in 78.6% of sporadic RCC and, although no intragenic mutations were detected, gene deletions and promoter region hypermethylation were detected in 17.3% and 20.3%, respectively, of sporadic RCC. Reexpression of UBE2QL1 in a deficient RCC cell line suppressed anchorage-independent growth. UBE2QL1 shows homology to the E2 class of ubiquitin conjugating enzymes and we found that (1) UBE2QL1 possesses an active-site cysteine (C88) that is monoubiquitinated in vivo, and (2) UBE2QL1 interacts with FBXW7 (an F box protein providing substrate recognition to the SCF E3 ubiquitin ligase) and facilitates the degradation of the known FBXW7 targets, CCNE1 and mTOR. These findings suggest UBE2QL1 as a novel candidate renal tumor suppressor gene.

A familial case with interstitial 2q36 deletion: variable phenotypic expression in full and mosaic state.

Submicroscopic chromosomal anomalies play an important role in the etiology of craniofacial malformations, including midline facial defects with hypertelorism (MFDH). MFDH is a common feature combination in several conditions, of which Frontonasal Dysplasia is the most frequently encountered manifestation; in most cases the etiology remains unknown. We identified a parent to child transmission of a 6.2 Mb interstitial deletion of chromosome region 2q36.1q36.3 by array-CGH and confirmed by FISH and microsatellite analysis. The patient and her mother both presented an MFDH phenotype although the phenotype in the mother was much milder than her daughter. Inspection of haplotype segregation within the family of 2q36.1 region suggests that the deletion arose on a chromosome derived from the maternal grandfather. Evidences based on FISH, microsatellite and array-CGH analysis point to a high frequency mosaicism for presence of a deleted region 2q36 occurring in blood of the mother. The frequency of mosaicism in other tissues could not be determined. We here suggest that the milder phenotype observed in the proband's mother can be explained by the mosaic state of the deletion. This most likely arose by an early embryonic deletion in the maternal embryo resulting in both gonadal and somatic mosaicism of two cell lines, with and without the deleted chromosome. The occurrence of gonadal mosaicism increases the recurrence risk significantly and is often either underestimated or not even taken into account in genetic counseling where new mutation is suspected.

Genomic and genic deletions of the FOX gene cluster on 16q24.1 and inactivating mutations of FOXF1 cause alveolar capillary dysplasia and other malformations.

Alveolar capillary dysplasia with misalignment of pulmonary veins (ACD/MPV) is a rare, neonatally lethal developmental disorder of the lung with defining histologic abnormalities typically associated with multiple congenital anomalies (MCA). Using array CGH analysis, we have identified six overlapping microdeletions encompassing the FOX transcription factor gene cluster in chromosome 16q24.1q24.2 in patients with ACD/MPV and MCA. Subsequently, we have identified four different heterozygous mutations (frameshift, nonsense, and no-stop) in the candidate FOXF1 gene in unrelated patients with sporadic ACD/MPV and MCA. Custom-designed, high-resolution microarray analysis of additional ACD/MPV samples revealed one microdeletion harboring FOXF1 and two distinct microdeletions upstream of FOXF1, implicating a position effect. DNA sequence analysis revealed that in six of nine deletions, both breakpoints occurred in the portions of Alu elements showing eight to 43 base pairs of perfect microhomology, suggesting replication error Microhomology-Mediated Break-Induced Replication (MMBIR)/Fork Stalling and Template Switching (FoSTeS) as a mechanism of their formation. In contrast to the association of point mutations in FOXF1 with bowel malrotation, microdeletions of FOXF1 were associated with hypoplastic left heart syndrome and gastrointestinal atresias, probably due to haploinsufficiency for the neighboring FOXC2 and FOXL1 genes. These differences reveal the phenotypic consequences of gene alterations in cis.

Characterization of a 3;6 translocation associated with renal cell carcinoma.

The most frequent cause of familial clear cell renal cell carcinoma (RCC) is von Hippel-Lindau disease and the VHL tumor suppressor gene (TSG) is inactivated in most sporadic clear cell RCC. Although there is relatively little information on the mechanisms of tumorigenesis of clear cell RCC without VHL inactivation, a subset of familial cases harbors a balanced constitutional chromosome 3 translocation. To date nine different chromosome 3 translocations have been associated with familial or multicentric clear cell RCC; and in three cases chromosome 6 was also involved. To identify candidate genes for renal tumorigenesis we characterized a constitutional translocation, t(3;6)(q22;q16.1) associated with multicentric RCC without evidence of VHL target gene dysregulation. Analysis of breakpoint sequences revealed a 1.3-kb deletion on chromosome 6 within the intron of a 2 exon predicted gene (NT_007299.434). However, RT-PCR analysis failed to detect the expression of this gene in lymphoblast, fibroblast, or kidney tumor cell lines. No known genes were disrupted by the translocation breakpoints but several candidate TSGs (e.g., EPHB1, EPHA7, PPP2R3A RNF184, and STAG1) map within close proximity to the breakpoints.

Requirement for PI 3-kinase gamma in macrophage migration to MCP-1 and CSF-1.

Phosphoinositide 3-kinases (PI3Ks) are important regulators of cell migration. The PI3K isoform gamma is primarily expressed in haematopoietic cells, and is activated by G protein-coupled receptors (GPCRs). Here, we investigate the contribution of PI3Kgamma to macrophage responses to chemoattractants, using bone marrow-derived macrophages from wild-type and PI3Kgamma-null mice. We observe that early membrane ruffling induced by MCP-1, which activates a GPCR, or by CSF-1, which activates a tyrosine kinase receptor, is unaltered in PI3Kgamma(-/-) mice, although by 30 min MCP-1-induced cell polarization was strongly reduced in PI3Kgamma(-/-) compared to wild-type macrophages. The migration behaviour of the macrophages was analysed by time-lapse microscopy in Dunn chemotaxis chambers. PI3Kgamma(-/-) macrophages showed reduced migration speed and translocation, and no chemotaxis to MCP-1. Interestingly, there was also a reduction in migration efficiency in PI3Kgamma(-/-) macrophages stimulated with CSF-1 although early CSF-1R signalling was normal. These results indicate that the initial actin reorganization induced by either a GPCR or tyrosine kinase receptor agonist is not dependent on PI3Kgamma, whereas PI3Kgamma is needed for optimal migration of macrophages to either agonist.