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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Functional and genetic analysis of choroid plexus development in zebrafish.

The choroid plexus, an epithelial-based structure localized in the brain ventricle, is the major component of the blood-cerebrospinal fluid barrier. The choroid plexus produces the cerebrospinal fluid and regulates the components of the cerebrospinal fluid. Abnormal choroid plexus function is associated with neurodegenerative diseases, tumor formation in the choroid plexus epithelium, and hydrocephaly. In this study, we used zebrafish (Danio rerio) as a model system to understand the genetic components of choroid plexus development. We generated an enhancer trap line, Et(cp:EGFP) (sj2), that expresses enhanced green fluorescent protein (EGFP) in the choroid plexus epithelium. Using immunohistochemistry and fluorescent tracers, we demonstrated that the zebrafish choroid plexus possesses brain barrier properties such as tight junctions and transporter activity. Thus, we have established zebrafish as a functionally relevant model to study choroid plexus development. Using an unbiased approach, we performed a forward genetic dissection of the choroid plexus to identify genes essential for its formation and function. Using Et(cp:EGFP) (sj2), we isolated 10 recessive mutant lines with choroid plexus abnormalities, which were grouped into five classes based on GFP intensity, epithelial localization, and overall choroid plexus morphology. We also mapped the mutation for two mutant lines to chromosomes 4 and 21, respectively. The mutants generated in this study can be used to elucidate specific genes and signaling pathways essential for choroid plexus development, function, and/or maintenance and will provide important insights into how these genetic mutations contribute to disease.

The genomic landscape of diffuse intrinsic pontine glioma and pediatric non-brainstem high-grade glioma.

Pediatric high-grade glioma (HGG) is a devastating disease with a less than 20% survival rate 2 years after diagnosis. We analyzed 127 pediatric HGGs, including diffuse intrinsic pontine gliomas (DIPGs) and non-brainstem HGGs (NBS-HGGs), by whole-genome, whole-exome and/or transcriptome sequencing. We identified recurrent somatic mutations in ACVR1 exclusively in DIPGs (32%), in addition to previously reported frequent somatic mutations in histone H3 genes, TP53 and ATRX, in both DIPGs and NBS-HGGs. Structural variants generating fusion genes were found in 47% of DIPGs and NBS-HGGs, with recurrent fusions involving the neurotrophin receptor genes NTRK1, NTRK2 and NTRK3 in 40% of NBS-HGGs in infants. Mutations targeting receptor tyrosine kinase-RAS-PI3K signaling, histone modification or chromatin remodeling, and cell cycle regulation were found in 68%, 73% and 59% of pediatric HGGs, respectively, including in DIPGs and NBS-HGGs. This comprehensive analysis provides insights into the unique and shared pathways driving pediatric HGG within and outside the brainstem.

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.

Generation of fertile and diploid fish, medaka (Oryzias latipes), from nuclear transplantation of blastula and four-somite-stage embryonic cells into nonenucleated unfertilized eggs.

In two experimental series of transplantation of embryonic cell nuclei into nonenucleated unfertilized eggs in medaka (Oryzias latipes), fertile and diploid nuclear transplants were successfully generated. In the first experiment, nuclei from blastula cells of a medaka stock with the wild-type body color were transplanted into 1722 eggs from the orange-red variety. Of 26 adult nuclear transplants with the wild-type body color, 22 were, as expected, triploid and sterile, but the other four were fertile. Three of the four were diploid, and the last one was tetraploid. They transmitted the wild-type body color to the F1 and F2 progenies in a Mendelian fashion. In the second experiment, cell nuclei from four-somite-stage embryos of the orangered variety carrying the green fluorescent protein (GFP) transgene were transplanted into 1688 recipients of the same strain. Three adult nuclear transplants expressing GFP were obtained. Two of them were triploid and sterile, but the remaining one was fertile and diploid. The transgene of the donor nuclei was transmitted to the F(1) and F(2) offspring in a Mendelian fashion. These observations that diploid and fertile nuclear transplants could be obtained without enucleation of the recipient eggs may have important implications for future nuclear transplantation in medaka.

Recapitulation of fast skeletal muscle development in zebrafish by transgenic expression of GFP under the mylz2 promoter.

A 1,934-bp muscle-specific promoter from the zebrafish mylz2 gene was isolated and characterized by transgenic analysis. By using a series of 5' promoter deletions linked to the green fluorescent protein (gfp) reporter gene, transient transgenic analysis indicated that the strength of promoter activity appeared to correlate to the number of muscle cis-elements in the promoter and that a minimal -77-bp region was sufficient for a relatively strong promoter activity in muscle cells. Stable transgenic lines were obtained from several mylz2-gfp constructs. GFP expression in the 1,934-bp promoter transgenic lines mimicked well the expression pattern of endogenous mylz2 mRNA in both somitic muscle and nonsomitic muscles, including fin, eye, jaw, and gill muscles. An identical pattern of GFP expression, although at a much lower level, was observed from a transgenic line with a shorter 871-bp promoter. Our observation indicates that there is no distinct cis-element for activation of mylz2 in different skeletal muscles. Furthermore, RNA encoding a dominant negative form of cAMP-dependent protein kinase A was injected into mylz2-gfp transgenic embryos and GFP expression was significantly reduced due to an expanded slow muscle development at the expense of GFP-expressing fast muscle. The mylz2-gfp transgene was also transferred into two zebrafish mutants, spadetail and chordino, and several novel phenotypes in muscle development in these mutants were discovered.

Generation of two-color transgenic zebrafish using the green and red fluorescent protein reporter genes gfp and rfp.

Two tissue-specific promoters were used to express both green fluorescent protein (GFP) and red fluorescent protein (RFP) in transgenic zebrafish embryos. One promoter (CK), derived from a cytokeratin gene, is active specifically in skin epithelia in embryos, and the other promoter (MLC) from a muscle-specific gene encodes a myosin light chain 2 polypeptide. When the 2 promoters drove the 2 reporter genes to express in the same embryos, both genes were faithfully expressed in the respective tissues, skin or muscle. When the 2 fluorescent proteins were expressed in the same skin or muscle cells under the same promoter, GFP fluorescence appeared earlier than RFP fluorescence in both skin and muscle tissues, probably owing to a higher detection sensitivity of GFP. However, RFP appeared to be more stable as its fluorescence steadily increased during development. Finally, F(1) transgenic offspring were obtained expressing GFP in skin cells under the CK promoter and RFP in muscle cells under the MLC promoter. Our study demonstrates the feasibility of monitoring expression of multiple genes in different tissues in the same transgenic organism.