Single-cell analysis reveals the Comma-1D cell line as a unique model for mammary gland development and breast cancer.

The mammary gland epithelial tree contains two distinct cell populations, luminal and basal. The investigation of how this heterogeneity is developed and how it influences tumorigenesis has been hampered by the need to perform studies on these populations using animal models. Comma-1D is an immortalized mouse mammary epithelial cell line that has unique morphogenetic properties. By performing single-cell RNA-seq studies, we found that Comma-1D cultures consist of two main populations with luminal and basal features, and a smaller population with mixed lineage and bipotent characteristics. We demonstrated that multiple transcription factors associated with the differentiation of the mammary epithelium in vivo also modulate this process in Comma-1D cultures. Additionally, we found that only cells with luminal features were able to acquire transformed characteristics after an oncogenic HER2 (also known as ERBB2) mutant was introduced in their genomes. Overall, our studies characterize, at a single-cell level, the heterogeneity of the Comma-1D cell line and illustrate how Comma-1D cells can be used as an experimental model to study both the differentiation and the transformation processes in vitro.

Impact of T-cell immunity on chemotherapy response in childhood acute lymphoblastic leukemia.

Although acute lymphoblastic leukemia (ALL) is highly responsive to chemotherapy, it is unknown how or which host immune factors influence the long-term remission of this cancer. To this end, we systematically evaluated the effects of T-cell immunity on Ph+ ALL therapy outcomes. Using a murine Arf-/-BCR-ABL1 B-cell ALL model, we showed that loss of T cells in the host drastically increased leukemia relapse after dasatinib or cytotoxic chemotherapy. Although ABL1 mutations emerged early during dasatinib treatment in both immunocompetent and immunocompromised hosts, T-cell immunity was essential for suppressing the outgrowth of drug-resistant leukemia. Bulk and single-cell transcriptome profiling of T cells during therapy pointed to the activation of type 1 immunity-related cytokine signaling being linked to long-term leukemia remission in mice. Consistent with these observations, interferon γ and interleukin 12 directly modulated dasatinib antileukemia efficacy in vivo. Finally, we evaluated peripheral blood immune cell composition in 102 children with ALL during chemotherapy and observed a significant association of T-cell abundance with treatment outcomes. Together, these results suggest that T-cell immunity plays pivotal roles in maintaining long-term remission of ALL, highlighting that the interplay between host immunity and drug resistance can be harnessed to improve ALL chemotherapy outcomes.

A developmentally prometastatic niche to hepatoblastoma in neonatal liver mediated by the Cxcl1/Cxcr2 axis.

BACKGROUND AND AIMS: Hepatoblastoma (HB) is the most common pediatric liver cancer. Its predominant occurrence in very young children led us to investigate whether the neonatal liver provides a protumorigenic niche to HB development. APPROACH AND RESULTS: HB development was compared between orthotopic transplantation models established in postnatal day 5 (P5) and 60 (P60) mice (P5Tx and P60Tx models). Single-cell RNA-sequencing (sc-RNAseq) was performed using tumor and liver tissues from both models and the top candidate cell types and genes identified are investigated for their roles in HB cell growth, migration, and survival. CONCLUSIONS: We found that various HB cell lines including HepG2 cells were consistently and considerably more tumorigenic and metastatic in the P5Tx model than in the P60Tx models. Sc-RNAseq of the P5Tx and P60Tx HepG2 models revealed that the P5Tx tumor was more hypoxic and had a larger number of activated hepatic stellate cells (aHSCs) in the tumor-surrounding liver that express significantly higher levels of Cxcl1 than those from the P60Tx model. We found these differences were developmentally present in normal P5 and P60 liver. We showed that the Cxcl1/Cxcr2 axis mediated HB cell migration and was critical to HB cell survival under hypoxia. Treating HepG2 P60Tx model with recombinant CXCL1 protein induced intrahepatic and pulmonary metastasis and CXCR2 knockout (KO) in HepG2 cells abolished their metastatic potential in the P5Tx model. Lastly, we showed that in tumors from patients with metastatic HB, there was a similar larger population of aHSCs in the tumor-surrounding liver than in localized tumors, and tumor hypoxia was uniquely associated with prognosis of patients with HB among pediatric cancers. We demonstrated that the neonatal liver provides a prometastatic niche to HB development through the Cxcl1/Cxcr2 axis.

Active medulloblastoma enhancers reveal subgroup-specific cellular origins.

Medulloblastoma is a highly malignant paediatric brain tumour, often inflicting devastating consequences on the developing child. Genomic studies have revealed four distinct molecular subgroups with divergent biology and clinical behaviour. An understanding of the regulatory circuitry governing the transcriptional landscapes of medulloblastoma subgroups, and how this relates to their respective developmental origins, is lacking. Here, using H3K27ac and BRD4 chromatin immunoprecipitation followed by sequencing (ChIP-seq) coupled with tissue-matched DNA methylation and transcriptome data, we describe the active cis-regulatory landscape across 28 primary medulloblastoma specimens. Analysis of differentially regulated enhancers and super-enhancers reinforced inter-subgroup heterogeneity and revealed novel, clinically relevant insights into medulloblastoma biology. Computational reconstruction of core regulatory circuitry identified a master set of transcription factors, validated by ChIP-seq, that is responsible for subgroup divergence, and implicates candidate cells of origin for Group 4. Our integrated analysis of enhancer elements in a large series of primary tumour samples reveals insights into cis-regulatory architecture, unrecognized dependencies, and cellular origins.

Latent cellular analysis robustly reveals subtle diversity in large-scale single-cell RNA-seq data.

Single-cell RNA sequencing (scRNA-seq) is a powerful tool for characterizing the cell-to-cell variation and cellular dynamics in populations which appear homogeneous otherwise in basic and translational biological research. However, significant challenges arise in the analysis of scRNA-seq data, including the low signal-to-noise ratio with high data sparsity, potential batch effects, scalability problems when hundreds of thousands of cells are to be analyzed among others. The inherent complexities of scRNA-seq data and dynamic nature of cellular processes lead to suboptimal performance of many currently available algorithms, even for basic tasks such as identifying biologically meaningful heterogeneous subpopulations. In this study, we developed the Latent Cellular Analysis (LCA), a machine learning-based analytical pipeline that combines cosine-similarity measurement by latent cellular states with a graph-based clustering algorithm. LCA provides heuristic solutions for population number inference, dimension reduction, feature selection, and control of technical variations without explicit gene filtering. We show that LCA is robust, accurate, and powerful by comparison with multiple state-of-the-art computational methods when applied to large-scale real and simulated scRNA-seq data. Importantly, the ability of LCA to learn from representative subsets of the data provides scalability, thereby addressing a significant challenge posed by growing sample sizes in scRNA-seq data analysis.

CNS angiogenesis and barriergenesis occur simultaneously.

The blood-brain barrier (BBB) plays a vital role in the central nervous system (CNS). A comprehensive understanding of BBB development has been hampered by difficulties in observing the differentiation of brain endothelial cells (BECs) in real-time. Here, we generated two transgenic zebrafish line, Tg(glut1b:mCherry) and Tg(plvap:EGFP), to serve as in vivo reporters of BBB development. We showed that barriergenesis (i.e. the induction of BEC differentiation) occurs immediately as endothelial tips cells migrate into the brain parenchyma. Using the Tg(glut1b:mCherry) transgenic line, we performed a genetic screen and identified a zebrafish mutant with a nonsense mutation in gpr124, a gene known to play a role in CNS angiogenesis and BBB development. We also showed that our transgenic plvap:EGFP line, a reporter of immature brain endothelium, is initially expressed in newly formed brain endothelial cells, but subsides during BBB maturation. Our results demonstrate the ability to visualize the in vivo differentiation of brain endothelial cells into the BBB phenotype and establish that CNS angiogenesis and barriergenesis occur simultaneously.

Zebrafish Cacna1fa is required for cone photoreceptor function and synaptic ribbon formation.

Mutations in the human CACNA1F gene cause incomplete congenital stationary night blindness type 2 (CSNB2), a non-progressive, clinically heterogeneous retinal disorder. However, the molecular mechanisms underlying CSNB2 have not been fully explored. Here, we describe the positional cloning of a blind zebrafish mutant, wait until dark (wud), which encodes a zebrafish homolog of human CACNA1F. We identified two zebrafish cacna1f paralogs and showed that the cacna1fa transcript (the gene mutated in wud) is expressed exclusively in the photoreceptor layer. We demonstrated that Cacna1fa localizes at the photoreceptor synapse and is absent from wud mutants. Electroretinograms revealed abnormal cone photoreceptor responses from wud mutants, indicating a defect in synaptic transmission. Although there are no obvious morphological differences, we found that wud mutants lacked synaptic ribbons and that wud is essential for the development of synaptic ribbons. We found that Ribeye, the most prominent synaptic ribbon protein, was less abundant and mislocalized in adult wud mutants. In addition to cloning wud, we identified synaptojanin 1 (synj1) as the defective gene in slacker (slak), a blind mutant with floating synaptic ribbons. We determined that Cacna1fa was expressed in slak photoreceptors and that Synj1 was initially expressed wud photoreceptors, but was absent by 5 days postfertilization. Collectively, our data demonstrate that Cacna1fa is essential for cone photoreceptor function and synaptic ribbon formation and reveal a previously unknown yet critical role of L-type voltage-dependent calcium channels in the expression and/or distribution of synaptic ribbon proteins, providing a new model to study the clinical variability in human CSNB2 patients.

Oncogenic KRAS promotes malignant brain tumors in zebrafish.

BACKGROUND: Zebrafish have been used as a vertebrate model to study human cancers such as melanoma, rhabdomyosarcoma, liver cancer, and leukemia as well as for high-throughput screening of small molecules of therapeutic value. However, they are just emerging as a model for human brain tumors, which are among the most devastating and difficult to treat. In this study, we evaluated zebrafish as a brain tumor model by overexpressing a human version of oncogenic KRAS (KRAS(G12V)). METHODS: Using zebrafish cytokeratin 5 (krt5) and glial fibrillary acidic protein (gfap) gene promoters, we activated Ras signaling in the zebrafish central nervous system (CNS) through transient and stable transgenic overexpression. Immunohistochemical analyses were performed to identify activated pathways in the resulting brain tumors. The effects of the MEK inhibitor U0126 on oncogenic KRAS were evaluated. RESULTS: We demonstrated that transient transgenic expression of KRAS(G12V) in putative neural stem and/or progenitor cells induced brain tumorigenesis. When expressed under the control of the krt5 gene promoter, KRAS(G12V) induced brain tumors in ventricular zones (VZ) at low frequency. The majority of other tumors were composed mostly of spindle and epithelioid cells, reminiscent of malignant peripheral nerve sheath tumors (MPNSTs). In contrast, when expressed under the control of the gfap gene promoter, KRAS(G12V) induced brain tumors in both VZs and brain parenchyma at higher frequency. Immunohistochemical analyses indicated prominent activation of the canonical RAS-RAF-ERK pathway, variable activation of the mTOR pathway, but no activation of the PI3K-AKT pathway. In a krt5-derived stable and inducible transgenic line, expression of oncogenic KRAS resulted in skin hyperplasia, and the MEK inhibitor U0126 effectively suppressed this pro-proliferative effects. In a gfap-derived stable and inducible line, expression of oncogenic KRAS led to significantly increased mitotic index in the spinal cord. CONCLUSIONS: Our studies demonstrate that zebrafish could be explored to study cellular origins and molecular mechanisms of brain tumorigenesis and could also be used as a platform for studying human oncogene function and for discovering oncogenic RAS inhibitors.

Conditional control of gene function by an invertible gene trap in zebrafish.

Conditional mutations are essential for determining the stage- and tissue-specific functions of genes. Here we achieve conditional mutagenesis in zebrafish using FT1, a gene-trap cassette that can be stably inverted by both Cre and Flp recombinases. We demonstrate that intronic insertions in the gene-trapping orientation severely disrupt the expression of the host gene, whereas intronic insertions in the neutral orientation do not significantly affect host gene expression. Cre- and Flp-mediated recombination switches the orientation of the gene-trap cassette, permitting conditional rescue in one orientation and conditional knockout in the other. To illustrate the utility of this system we analyzed the functional consequence of intronic FT1 insertion in supv3l1, a gene encoding a mitochondrial RNA helicase. Global supv311 mutants have impaired mitochondrial function, embryonic lethality, and agenesis of the liver. Conditional rescue of supv311 expression in hepatocytes specifically corrected the liver defects. To test whether the liver function of supv311 is required for viability we used Flp-mediated recombination in the germline to generate a neutral allele at the locus. Subsequently, tissue-specific expression of Cre conditionally inactivated the targeted locus. Hepatocyte-specific inactivation of supv311 caused liver degeneration, growth retardation, and juvenile lethality, a phenotype that was less severe than the global disruption of supv311. Thus, supv311 is required in multiple tissues for organismal viability. Our mutagenesis approach is very efficient and could be used to generate conditional alleles throughout the zebrafish genome. Furthermore, because FT1 is based on the promiscuous Tol2 transposon, it should be applicable to many organisms.

Integrome signatures of lentiviral gene therapy for SCID-X1 patients.

Lentiviral vector (LV)-based gene therapy holds promise for a broad range of diseases. Analyzing more than 280,000 vector integration sites (VISs) in 273 samples from 10 patients with X-linked severe combined immunodeficiency (SCID-X1), we discovered shared LV integrome signatures in 9 of 10 patients in relation to the genomics, epigenomics, and 3D structure of the human genome. VISs were enriched in the nuclear subcompartment A1 and integrated into super-enhancers close to nuclear pore complexes. These signatures were validated in T cells transduced with an LV encoding a CD19-specific chimeric antigen receptor. Intriguingly, the one patient whose VISs deviated from the identified integrome signatures had a distinct clinical course. Comparison of LV and gamma retrovirus integromes regarding their 3D genome signatures identified differences that might explain the lower risk of insertional mutagenesis in LV-based gene therapy. Our findings suggest that LV integrome signatures, shaped by common features such as genome organization, may affect the efficacy of LV-based cellular therapies.

Convergent evolution and multi-wave clonal invasion in H3 K27-altered diffuse midline gliomas treated with a PDGFR inhibitor.

The majority of diffuse midline gliomas, H3 K27-altered (DMG-H3 K27-a), are infiltrating pediatric brain tumors that arise in the pons with no effective treatment. To understand how clonal evolution contributes to the tumor's invasive spread, we performed exome sequencing and SNP array profiling on 49 multi-region autopsy samples from 11 patients with pontine DMG-H3 K27-a enrolled in a phase I clinical trial of PDGFR inhibitor crenolanib. For each patient, a phylogenetic tree was constructed by testing multiple possible clonal evolution models to select the one consistent with somatic mutations and copy number variations across all tumor regions. The tree was then used to deconvolute subclonal composition and prevalence at each tumor region to study convergent evolution and invasion patterns. Somatic variants in the PI3K pathway, a late event, are enriched in our cohort, affecting 70% of patients. Convergent evolution of PI3K at distinct phylogenetic branches was detected in 40% of the patients. 24 (~ 50%) of tumor regions were occupied by subclones of mixed lineages with varying molecular ages, indicating multiple waves of invasion across the pons and extrapontine. Subclones harboring a PDGFRA amplicon, including one that amplified a PDGRFAY849C mutant allele, were detected in four patients; their presence in extrapontine tumor and normal brain samples imply their involvement in extrapontine invasion. Our study expands the current knowledge on tumor invasion patterns in DMG-H3 K27-a, which may inform the design of future clinical trials.

The miR-424(322)/503 gene cluster regulates pro- versus anti-inflammatory skin DC subset differentiation by modulating TGF-ß signaling.

Transforming growth factor ß (TGF-ß) family ligands are key regulators of dendritic cell (DC) differentiation and activation. Epidermal Langerhans cells (LCs) require TGF-ß family signaling for their differentiation, and canonical TGF-ß1 signaling secures a non-activated LC state. LCs reportedly control skin inflammation and are replenished from peripheral blood monocytes, which also give rise to pro-inflammatory monocyte-derived DCs (moDCs). By studying mechanisms in inflammation, we previously screened LCs versus moDCs for differentially expressed microRNAs (miRNAs). This revealed that miR-424/503 is the most strongly inversely regulated (moDCs > LCs). We here demonstrate that miR-424/503 is induced during moDC differentiation and promotes moDC differentiation in human and mouse. Inversely, forced repression of miR-424 during moDC differentiation facilitates TGF-ß1-dependent LC differentiation. Mechanistically, miR-424/503 deficiency in monocyte/DC precursors leads to the induction of TGF-ß1 response genes critical for LC differentiation. Therefore, the miR-424/503 gene cluster plays a decisive role in anti-inflammatory LC versus pro-inflammatory moDC differentiation from monocytes.

Integrative network analysis reveals USP7 haploinsufficiency inhibits E-protein activity in pediatric T-lineage acute lymphoblastic leukemia (T-ALL).

USP7, which encodes a deubiquitylating enzyme, is among the most frequently mutated genes in pediatric T-ALL, with somatic heterozygous loss-of-function mutations (haploinsufficiency) predominantly affecting the subgroup that has aberrant TAL1 oncogene activation. Network analysis of > 200 T-ALL transcriptomes linked USP7 haploinsufficiency with decreased activities of E-proteins. E-proteins are also negatively regulated by TAL1, leading to concerted down-regulation of E-protein target genes involved in T-cell development. In T-ALL cell lines, we showed the physical interaction of USP7 with E-proteins and TAL1 by mass spectrometry and ChIP-seq. Haploinsufficient but not complete CRISPR knock-out of USP7 showed accelerated cell growth and validated transcriptional down-regulation of E-protein targets. Our study unveiled the synergistic effect of USP7 haploinsufficiency with aberrant TAL1 activation on T-ALL, implicating USP7 as a haploinsufficient tumor suppressor in T-ALL. Our findings caution against a universal oncogene designation for USP7 while emphasizing the dosage-dependent consequences of USP7 inhibitors currently under development as potential cancer therapeutics.

Network-based systems pharmacology reveals heterogeneity in LCK and BCL2 signaling and therapeutic sensitivity of T-cell acute lymphoblastic leukemia.

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive hematological malignancy, and novel therapeutics are much needed. Profiling patient leukemia' drug sensitivities ex vivo, we discovered that 44.4% of childhood and 16.7% of adult T-ALL cases exquisitely respond to dasatinib. Applying network-based systems pharmacology analyses to examine signal circuitry, we identified preTCR-LCK activation as the driver of dasatinib sensitivity, and T-ALL-specific LCK dependency was confirmed in genome-wide CRISPR-Cas9 screens. Dasatinib-sensitive T-ALLs exhibited high BCL-XL and low BCL2 activity and venetoclax resistance. Discordant sensitivity of T-ALL to dasatinib and venetoclax is strongly correlated with T-cell differentiation, particularly with the dynamic shift in LCK vs. BCL2 activation. Finally, single-cell analysis identified leukemia heterogeneity in LCK and BCL2 signaling and T-cell maturation stage, consistent with dasatinib response. In conclusion, our results indicate that developmental arrest in T-ALL drives differential activation of preTCR-LCK and BCL2 signaling in this leukemia, providing unique opportunities for targeted therapy.

A comparison of methods accounting for batch effects in differential expression analysis of UMI count based single cell RNA sequencing.

Accounting for batch effects, especially latent batch effects, in differential expression (DE) analysis is critical for identifying true biological effects. Single-cell RNA sequencing (scRNA-seq) is a powerful tool for quantifying cell-to-cell variation in transcript abundance and characterizing cellular dynamics. Although many scRNA-seq DE analysis methods accommodate known batch variables, their performance has not been systematically evaluated. Moreover, the challenge of accounting for latent batch variables in scRNA-seq DE analysis is largely unmet. In contrast, many methods have been developed to account for batch variables (either known or latent) in other high-dimensional data, especially bulk RNA-seq. We extensively evaluate 11 methods for batch variables in different scRNA-seq DE analysis scenarios, with a primary focus on latent batch variables. We demonstrate that for known batch variables, incorporating them as covariates into a regression model outperformed approaches using a batch-corrected matrix. For latent batches, fixed effects models have inflated FDRs, whereas aggregation-based methods and mixed effects models have significant power loss. Surrogate variable based methods generally control the FDR well while achieving good power with small group effects. However, their performance (except that of SVA) deteriorated substantially in scenarios involving large group effects and/or group label impurity. In these settings, SVA achieves relatively good performance despite an occasionally inflated FDR (up to 0.2). Finally we make the following recommendations for scRNA-seq DE analysis: 1) incorporate known batch variables instead of using batch-corrected data; and 2) employ SVA for latent batch correction. However, better methods are still needed to fully unleash the power of scRNA-seq.

Discovery of regulatory noncoding variants in individual cancer genomes by using cis-X.

We developed cis-X, a computational method for discovering regulatory noncoding variants in cancer by integrating whole-genome and transcriptome sequencing data from a single cancer sample. cis-X first finds aberrantly cis-activated genes that exhibit allele-specific expression accompanied by an elevated outlier expression. It then searches for causal noncoding variants that may introduce aberrant transcription factor binding motifs or enhancer hijacking by structural variations. Analysis of 13 T-lineage acute lymphoblastic leukemias identified a recurrent intronic variant predicted to cis-activate the TAL1 oncogene, a finding validated in vivo by chromatin immunoprecipitation sequencing of a patient-derived xenograft. Candidate oncogenes include the prolactin receptor PRLR activated by a focal deletion that removes a CTCF-insulated neighborhood boundary. cis-X may be applied to pediatric and adult solid tumors that are aneuploid and heterogeneous. In contrast to existing approaches, which require large sample cohorts, cis-X enables the discovery of regulatory noncoding variants in individual cancer genomes.

Co-activation of hedgehog and AKT pathways promote tumorigenesis in zebrafish.

The zebrafish has become an important model for cancer research. Several cancer models have been established by transgenic expression of human or mouse oncogenes in zebrafish. Since it is amenable to efficient transgenesis, zebrafish have immense potential to be used for studying interaction of oncogenes and pathways at the organismal level. Using the Gal4VP16-UAS binary transgenic expression approach, we established stable transgenic lines expressing an EGFP fusion protein of an activated zebrafish Smoothened (Smoa1-EGFP). Expression of the zebrafish Smoa1-EGFP itself did not lead to tumor formation either in founder fish or subsequent generations, however, co-expressing a constitutively active human AKT1 resulted in several tumor types, including spindle cell sarcoma, rhabdomyoma, ocular melanoma, astrocytoma, and myxoma. All tumor types showed GFP expression and increased Patched 1 levels, suggesting involvement of zebrafish Smoa1 in tumorigenesis. Immunofluorescence studies showed that tumors also expressed elevated levels of phosphorylated AKT, indicating activation of the PI3K-AKT pathway. These results suggest that co-activation of the hedgehog and AKT pathways promote tumorigenesis, and that the binary transgenic approach is a useful tool for studying interaction of oncogenes and oncogenic pathways in zebrafish.

Long-read sequencing unveils IGH-DUX4 translocation into the silenced IGH allele in B-cell acute lymphoblastic leukemia.

IGH@ proto-oncogene translocation is a common oncogenic event in lymphoid lineage cancers such as B-ALL, lymphoma and multiple myeloma. Here, to investigate the interplay between IGH@ proto-oncogene translocation and IGH allelic exclusion, we perform long-read whole-genome and transcriptome sequencing along with epigenetic and 3D genome profiling of Nalm6, an IGH-DUX4 positive B-ALL cell line. We detect significant allelic imbalance on the wild-type over the IGH-DUX4 haplotype in expression and epigenetic data, showing IGH-DUX4 translocation occurs on the silenced IGH allele. In vitro, this reduces the oncogenic stress of DUX4 high-level expression. Moreover, patient samples of IGH-DUX4 B-ALL have similar expression profile and IGH breakpoints as Nalm6, suggesting a common mechanism to allow optimal dosage of non-toxic DUX4 expression.

Cloned zebrafish by nuclear transfer from long-term-cultured cells.

Although mammals have been cloned from genetically manipulated cultured cells, a comparable achievement has not been realized in lower vertebrates. Here we report that fertile transgenic zebrafish can be obtained by nuclear transfer using embryonic fibroblast cells from long-term cultures. The donor nuclei, modified by retroviral insertions expressing green fluorescent protein (GFP), were transplanted into manually enucleated eggs. Overall, a 2% success rate was achieved, resulting in 11 adult transgenic zebrafish expressing GFP. These nuclear transplants produced fertile, diploid offspring, and their F1/F2 progeny continued to express GFP in a pattern identical to that of the founder fish. This finding demonstrates that slowly dividing nuclei from cultured cells can be reprogrammed to support rapid embryonic development and sets up a foundation for targeted genetic manipulation in zebrafish.