A Highly Scalable Method for Joint Whole-Genome Sequencing and Gene-Expression Profiling of Single Cells.

To address how genetic variation alters gene expression in complex cell mixtures, we developed direct nuclear tagmentation and RNA sequencing (DNTR-seq), which enables whole-genome and mRNA sequencing jointly in single cells. DNTR-seq readily identified minor subclones within leukemia patients. In a large-scale DNA damage screen, DNTR-seq was used to detect regions under purifying selection and identified genes where mRNA abundance was resistant to copy-number alteration, suggesting strong genetic compensation. mRNA sequencing (mRNA-seq) quality equals RNA-only methods, and the low positional bias of genomic libraries allowed detection of sub-megabase aberrations at ultra-low coverage. Each cell library is individually addressable and can be re-sequenced at increased depth, allowing multi-tiered study designs. Additionally, the direct tagmentation protocol enables coverage-independent estimation of ploidy, which can be used to identify cell singlets. Thus, DNTR-seq directly links each cell's state to its corresponding genome at scale, enabling routine analysis of heterogeneous tumors and other complex tissues.

Joint single-cell genetic and transcriptomic analysis reveal pre-malignant SCP-like subclones in human neuroblastoma.

BACKGROUND: Neuroblastoma (NB) is a heterogeneous embryonal malignancy and the deadliest tumor of infancy. It is a complex disease that can result in diverse clinical outcomes. In some children, tumors regress spontaneously. Others respond well to existing treatments. But for the high-risk group, which constitutes approximately 40% of all patients, the prognosis remains dire despite collaborative efforts in basic and clinical research. While its exact cellular origin is still under debate, NB is assumed to arise from the neural crest cell lineage including multipotent Schwann cell precursors (SCPs), which differentiate into sympatho-adrenal cell states eventually producing chromaffin cells and sympathoblasts. METHODS: To investigate clonal development of neuroblastoma cell states, we performed haplotype-specific analysis of human tumor samples using single-cell multi-omics, including joint DNA/RNA sequencing of sorted single cells (DNTR-seq). Samples were also assessed using immunofluorescence stainings and fluorescence in-situ hybridization (FISH). RESULTS: Beyond adrenergic tumor cells, we identify subpopulations of aneuploid SCP-like cells, characterized by clonal expansion, whole-chromosome 17 gains, as well as expression programs of proliferation, apoptosis, and a non-immunomodulatory phenotype. CONCLUSION: Aneuploid pre-malignant SCP-like cells represent a novel feature of NB. Genetic evidence and tumor phylogeny suggest that these clones and malignant adrenergic populations originate from aneuploidy-prone cells of migrating neural crest or SCP origin, before lineage commitment to sympatho-adrenal cell states. Our findings expand the phenotypic spectrum of NB cell states. Considering the multipotency of SCPs in development, we suggest that the transformation of fetal SCPs may represent one possible mechanism of tumor initiation in NB with chromosome 17 aberrations as a characteristic element.

An overlooked subset of Cx3cr1wt/wt microglia in the Cx3cr1CreER-Eyfp/wt mouse has a repopulation advantage over Cx3cr1CreER-Eyfp/wt microglia following microglial depletion.

BACKGROUND: Fluorescent reporter labeling and promoter-driven Cre-recombinant technologies have facilitated cellular investigations of physiological and pathological processes, including the widespread use of the Cx3cr1CreER-Eyfp/wt mouse strain for studies of microglia. METHODS: Immunohistochemistry, Flow Cytometry, RNA sequencing and whole-genome sequencing were used to identify the subpopulation of microglia in Cx3cr1CreER-Eyfp/wt mouse brains. Genetically mediated microglia depletion using Cx3cr1CreER-Eyfp/wtRosa26DTA/wt mice and CSF1 receptor inhibitor PLX3397 were used to deplete microglia. Primary microglia proliferation and migration assay were used for in vitro studies. RESULTS: We unexpectedly identified a subpopulation of microglia devoid of genetic modification, exhibiting higher Cx3cr1 and CX3CR1 expression than Cx3cr1CreER-Eyfp/wtCre+Eyfp+ microglia in Cx3cr1CreER-Eyfp/wt mouse brains, thus termed Cx3cr1highCre-Eyfp- microglia. This subpopulation constituted less than 1% of all microglia under homeostatic conditions, but after Cre-driven DTA-mediated microglial depletion, Cx3cr1highCre-Eyfp- microglia escaped depletion and proliferated extensively, eventually occupying one-third of the total microglial pool. We further demonstrated that the Cx3cr1highCre-Eyfp- microglia had lost their genetic heterozygosity and become homozygous for wild-type Cx3cr1. Therefore, Cx3cr1highCre-Eyfp- microglia are Cx3cr1wt/wtCre-Eyfp-. Finally, we demonstrated that CX3CL1-CX3CR1 signaling regulates microglial repopulation both in vivo and in vitro. CONCLUSIONS: Our results raise a cautionary note regarding the use of Cx3cr1CreER-Eyfp/wt mouse strains, particularly when interpreting the results of fate mapping, and microglial depletion and repopulation studies.

Elastic dosage compensation by X-chromosome upregulation.

X-chromosome inactivation and X-upregulation are the fundamental modes of chromosome-wide gene regulation that collectively achieve dosage compensation in mammals, but the regulatory link between the two remains elusive and the X-upregulation dynamics are unknown. Here, we use allele-resolved single-cell RNA-seq combined with chromatin accessibility profiling and finely dissect their separate effects on RNA levels during mouse development. Surprisingly, we uncover that X-upregulation elastically tunes expression dosage in a sex- and lineage-specific manner, and moreover along varying degrees of X-inactivation progression. Male blastomeres achieve X-upregulation upon zygotic genome activation while females experience two distinct waves of upregulation, upon imprinted and random X-inactivation; and ablation of Xist impedes female X-upregulation. Female cells carrying two active X chromosomes lack upregulation, yet their collective RNA output exceeds that of a single hyperactive allele. Importantly, this conflicts the conventional dosage compensation model in which naïve female cells are initially subject to biallelic X-upregulation followed by X-inactivation of one allele to correct the X dosage. Together, our study provides key insights to the chain of events of dosage compensation, explaining how transcript copy numbers can remain remarkably stable across developmental windows wherein severe dose imbalance would otherwise be experienced by the cell.

Identification of functionally distinct and interacting cancer cell subpopulations from glioblastoma with intratumoral genetic heterogeneity.

BACKGROUND: Glioblastomas display a high level of intratumoral heterogeneity with regard to both genetic and histological features. Within single tumors, subclones have been shown to communicate with each other to affect overall tumor growth. The aim of this study was to broaden the understanding of interclonal communication in glioblastoma. METHODS: We have used the U-343 model, consisting of U-343 MG, U-343 MGa, U-343 MGa 31L, and U-343 MGa Cl2:6, a set of distinct glioblastoma cell lines that have been derived from the same tumor. We characterized these with regard to temozolomide sensitivity, protein secretome, gene expression, DNA copy number, and cancer cell phenotypic traits. Furthermore, we performed coculture and conditioned media-based experiments to model cell-to-cell signaling in a setting of intratumoral heterogeneity. RESULTS: Temozolomide treatment of a coculture composed of all 4 U-343 cell lines presents a tumor relapse model where the least sensitive population, U-343 MGa 31L, outlives the others. Interestingly, the U-343 cell lines were shown to have distinct gene expression signatures and phenotypes although they were derived from a single tumor. The DNA copy number analysis revealed both common and unique alterations, indicating the evolutionary relationship between the cells. Moreover, these cells were found to communicate and affect each other's proliferation, both via contact-dependent and -independent interactions, where NOTCH1, TGFBI, and ADAMTS1 signaling effects were involved, respectively. CONCLUSIONS: These results provide insight into how complex the signaling events may prove to be in a setting of intratumoral heterogeneity in glioblastoma and provide a map for future studies.