Combination anti-PD-1 and anti-CTLA-4 therapy generates waves of clonal responses that include progenitor-exhausted CD8+ T cells.

Combination checkpoint blockade with anti-PD-1 and anti-CTLA-4 antibodies has shown promising efficacy in melanoma. However, the underlying mechanism in humans remains unclear. Here, we perform paired single-cell RNA and T cell receptor (TCR) sequencing across time in 36 patients with stage IV melanoma treated with anti-PD-1, anti-CTLA-4, or combination therapy. We develop the algorithm Cyclone to track temporal clonal dynamics and underlying cell states. Checkpoint blockade induces waves of clonal T cell responses that peak at distinct time points. Combination therapy results in greater magnitude of clonal responses at 6 and 9 weeks compared to single-agent therapies, including melanoma-specific CD8+ T cells and exhausted CD8+ T cell (TEX) clones. Focused analyses of TEX identify that anti-CTLA-4 induces robust expansion and proliferation of progenitor TEX, which synergizes with anti-PD-1 to reinvigorate TEX during combination therapy. These next generation immune profiling approaches can guide the selection of drugs, schedule, and dosing for novel combination strategies.

Mapping the cellular biogeography of human bone marrow niches using single-cell transcriptomics and proteomic imaging.

Non-hematopoietic cells are essential contributors to hematopoiesis. However, heterogeneity and spatial organization of these cells in human bone marrow remain largely uncharacterized. We used single-cell RNA sequencing (scRNA-seq) to profile 29,325 non-hematopoietic cells and discovered nine transcriptionally distinct subtypes. We simultaneously profiled 53,417 hematopoietic cells and predicted their interactions with non-hematopoietic subsets. We employed co-detection by indexing (CODEX) to spatially profile over 1.2 million cells. We integrated scRNA-seq and CODEX data to link predicted cellular signaling with spatial proximity. Our analysis revealed a hyperoxygenated arterio-endosteal neighborhood for early myelopoiesis, and an adipocytic localization for early hematopoietic stem and progenitor cells (HSPCs). We used our CODEX atlas to annotate new images and uncovered mesenchymal stromal cell (MSC) expansion and spatial neighborhoods co-enriched for leukemic blasts and MSCs in acute myeloid leukemia (AML) patient samples. This spatially resolved, multiomic atlas of human bone marrow provides a reference for investigation of cellular interactions that drive hematopoiesis.

Mapping the Cellular Biogeography of Human Bone Marrow Niches Using Single-Cell Transcriptomics and Proteomic Imaging.

The bone marrow is the organ responsible for blood production. Diverse non-hematopoietic cells contribute essentially to hematopoiesis. However, these cells and their spatial organization remain largely uncharacterized as they have been technically challenging to study in humans. Here, we used fresh femoral head samples and performed single-cell RNA sequencing (scRNA-Seq) to profile 29,325 enriched non-hematopoietic bone marrow cells and discover nine transcriptionally distinct subtypes. We next employed CO-detection by inDEXing (CODEX) multiplexed imaging of 18 individuals, including both healthy and acute myeloid leukemia (AML) samples, to spatially profile over one million single cells with a novel 53-antibody panel. We discovered a relatively hyperoxygenated arterio-endosteal niche for early myelopoiesis, and an adipocytic, but not endosteal or perivascular, niche for early hematopoietic stem and progenitor cells. We used our atlas to predict cell type labels in new bone marrow images and used these predictions to uncover mesenchymal stromal cell (MSC) expansion and leukemic blast/MSC-enriched spatial neighborhoods in AML patient samples. Our work represents the first comprehensive, spatially-resolved multiomic atlas of human bone marrow and will serve as a reference for future investigation of cellular interactions that drive hematopoiesis.

A longitudinal single-cell and spatial multiomic atlas of pediatric high-grade glioma.

Pediatric high-grade glioma (pHGG) is an incurable central nervous system malignancy that is a leading cause of pediatric cancer death. While pHGG shares many similarities to adult glioma, it is increasingly recognized as a molecularly distinct, yet highly heterogeneous disease. In this study, we longitudinally profiled a molecularly diverse cohort of 16 pHGG patients before and after standard therapy through single-nucleus RNA and ATAC sequencing, whole-genome sequencing, and CODEX spatial proteomics to capture the evolution of the tumor microenvironment during progression following treatment. We found that the canonical neoplastic cell phenotypes of adult glioblastoma are insufficient to capture the range of tumor cell states in a pediatric cohort and observed differential tumor-myeloid interactions between malignant cell states. We identified key transcriptional regulators of pHGG cell states and did not observe the marked proneural to mesenchymal shift characteristic of adult glioblastoma. We showed that essential neuromodulators and the interferon response are upregulated post-therapy along with an increase in non-neoplastic oligodendrocytes. Through in vitro pharmacological perturbation, we demonstrated novel malignant cell-intrinsic targets. This multiomic atlas of longitudinal pHGG captures the key features of therapy response that support distinction from its adult counterpart and suggests therapeutic strategies which are targeted to pediatric gliomas.

Transcriptional immunogenomic analysis reveals distinct immunological clusters in paediatric nervous system tumours.

BACKGROUND: Cancer immunotherapies including immune checkpoint inhibitors and Chimeric Antigen Receptor (CAR) T-cell therapy have shown variable response rates in paediatric patients highlighting the need to establish robust biomarkers for patient selection. While the tumour microenvironment in adults has been widely studied to delineate determinants of immune response, the immune composition of paediatric solid tumours remains relatively uncharacterized calling for investigations to identify potential immune biomarkers. METHODS: To inform immunotherapy approaches in paediatric cancers with embryonal origin, we performed an immunogenomic analysis of RNA-seq data from 925 treatment-naïve paediatric nervous system tumours (pedNST) spanning 12 cancer types from three publicly available data sets. RESULTS: Within pedNST, we uncovered four broad immune clusters: Paediatric Inflamed (10%), Myeloid Predominant (30%), Immune Neutral (43%) and Immune Desert (17%). We validated these clusters using immunohistochemistry, methylation immune inference and segmentation analysis of tissue images. We report shared biology of these immune clusters within and across cancer types, and characterization of specific immune cell frequencies as well as T- and B-cell repertoires. We found no associations between immune infiltration levels and tumour mutational burden, although molecular cancer entities were enriched within specific immune clusters. CONCLUSIONS: Given the heterogeneity of immune infiltration within pedNST, our findings suggest personalized immunogenomic profiling is needed to guide selection of immunotherapeutic strategies.

Genetic predisposition to neuroblastoma results from a regulatory polymorphism that promotes the adrenergic cell state.

Childhood neuroblastomas exhibit plasticity between an undifferentiated neural crest-like mesenchymal cell state and a more differentiated sympathetic adrenergic cell state. These cell states are governed by autoregulatory transcriptional loops called core regulatory circuitries (CRCs), which drive the early development of sympathetic neuronal progenitors from migratory neural crest cells during embryogenesis. The adrenergic cell identity of neuroblastoma requires LMO1 as a transcriptional cofactor. Both LMO1 expression levels and the risk of developing neuroblastoma in children are associated with a single nucleotide polymorphism, G/T, that affects a GATA motif in the first intron of LMO1. Here, we showed that WT zebrafish with the GATA genotype developed adrenergic neuroblastoma, while knock-in of the protective TATA allele at this locus reduced the penetrance of MYCN-driven tumors, which were restricted to the mesenchymal cell state. Whole genome sequencing of childhood neuroblastomas demonstrated that TATA/TATA tumors also exhibited a mesenchymal cell state and were low risk at diagnosis. Thus, conversion of the regulatory GATA to a TATA allele in the first intron of LMO1 reduced the neuroblastoma-initiation rate by preventing formation of the adrenergic cell state. This mechanism was conserved over 400 million years of evolution, separating zebrafish and humans.

Genetic Predisposition to Neuroblastoma Results from a Regulatory Polymorphism that Promotes the Adrenergic Cell State.

Childhood neuroblastomas exhibit plasticity between an undifferentiated neural crest-like "mesenchymal" cell state and a more differentiated sympathetic "adrenergic" cell state. These cell states are governed by autoregulatory transcriptional loops called core regulatory circuitries (CRCs), which drive the early development of sympathetic neuronal progenitors from migratory neural crest cells during embryogenesis. The adrenergic cell identity of neuroblastoma requires LMO1 as a transcriptional co-factor. Both LMO1 expression levels and the risk of developing neuroblastoma in children are associated with a single nucleotide polymorphism G/T that affects a G ATA motif in the first intron of LMO1. Here we show that wild-type zebrafish with the G ATA genotype develop adrenergic neuroblastoma, while knock-in of the protective T ATA allele at this locus reduces the penetrance of MYCN-driven tumors, which are restricted to the mesenchymal cell state. Whole genome sequencing of childhood neuroblastomas demonstrates that T ATA/ T ATA tumors also exhibit a mesenchymal cell state and are low risk at diagnosis. Thus, conversion of the regulatory G ATA to a T ATA allele in the first intron of LMO1 reduces the neuroblastoma initiation rate by preventing formation of the adrenergic cell state, a mechanism that is conserved over 400 million years of evolution separating zebrafish and humans.

Human epigenetic and transcriptional T cell differentiation atlas for identifying functional T cell-specific enhancers.

The clinical benefit of T cell immunotherapies remains limited by incomplete understanding of T cell differentiation and dysfunction. We generated an epigenetic and transcriptional atlas of T cell differentiation from healthy humans that included exhausted CD8 T cells and applied this resource in three ways. First, we identified modules of gene expression and chromatin accessibility, revealing molecular coordination of differentiation after activation and between central memory and effector memory. Second, we applied this healthy molecular framework to three settings-a neoadjuvant anti-PD1 melanoma trial, a basal cell carcinoma scATAC-seq dataset, and autoimmune disease-associated SNPs-yielding insights into disease-specific biology. Third, we predicted genome-wide cis-regulatory elements and validated this approach for key effector genes using CRISPR interference, providing functional annotation and demonstrating the ability to identify targets for non-coding cellular engineering. These studies define epigenetic and transcriptional regulation of human T cells and illustrate the utility of interrogating disease in the context of a healthy T cell atlas.

Rapid induction of antigen-specific CD4+ T cells is associated with coordinated humoral and cellular immunity to SARS-CoV-2 mRNA vaccination.

SARS-CoV-2 mRNA vaccines have shown remarkable clinical efficacy, but questions remain about the nature and kinetics of T cell priming. We performed longitudinal antigen-specific T cell analyses on healthy SARS-CoV-2-naive and recovered individuals prior to and following mRNA prime and boost vaccination. Vaccination induced rapid antigen-specific CD4+ T cell responses in naive subjects after the first dose, whereas CD8+ T cell responses developed gradually and were variable in magnitude. Vaccine-induced Th1 and Tfh cell responses following the first dose correlated with post-boost CD8+ T cells and neutralizing antibodies, respectively. Integrated analysis revealed coordinated immune responses with distinct trajectories in SARS-CoV-2-naive and recovered individuals. Last, whereas booster vaccination improved T cell responses in SARS-CoV-2-naive subjects, the second dose had little effect in SARS-CoV-2-recovered individuals. These findings highlight the role of rapidly primed CD4+ T cells in coordinating responses to the second vaccine dose in SARS-CoV-2-naive individuals.

CD8+ T cells contribute to survival in patients with COVID-19 and hematologic cancer.

Patients with cancer have high mortality from coronavirus disease 2019 (COVID-19), and the immune parameters that dictate clinical outcomes remain unknown. In a cohort of 100 patients with cancer who were hospitalized for COVID-19, patients with hematologic cancer had higher mortality relative to patients with solid cancer. In two additional cohorts, flow cytometric and serologic analyses demonstrated that patients with solid cancer and patients without cancer had a similar immune phenotype during acute COVID-19, whereas patients with hematologic cancer had impairment of B cells and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-specific antibody responses. Despite the impaired humoral immunity and high mortality in patients with hematologic cancer who also have COVID-19, those with a greater number of CD8 T cells had improved survival, including those treated with anti-CD20 therapy. Furthermore, 77% of patients with hematologic cancer had detectable SARS-CoV-2-specific T cell responses. Thus, CD8 T cells might influence recovery from COVID-19 when humoral immunity is deficient. These observations suggest that CD8 T cell responses to vaccination might provide protection in patients with hematologic cancer even in the setting of limited humoral responses.

CD8 T cells compensate for impaired humoral immunity in COVID-19 patients with hematologic cancer.

Cancer patients have increased morbidity and mortality from Coronavirus Disease 2019 (COVID-19), but the underlying immune mechanisms are unknown. In a cohort of 100 cancer patients hospitalized for COVID-19 at the University of Pennsylvania Health System, we found that patients with hematologic cancers had a significantly higher mortality relative to patients with solid cancers after accounting for confounders including ECOG performance status and active cancer status. We performed flow cytometric and serologic analyses of 106 cancer patients and 113 non-cancer controls from two additional cohorts at Penn and Memorial Sloan Kettering Cancer Center. Patients with solid cancers exhibited an immune phenotype similar to non-cancer patients during acute COVID-19 whereas patients with hematologic cancers had significant impairment of B cells and SARS-CoV-2-specific antibody responses. High dimensional analysis of flow cytometric data revealed 5 distinct immune phenotypes. An immune phenotype characterized by CD8 T cell depletion was associated with a high viral load and the highest mortality of 71%, among all cancer patients. In contrast, despite impaired B cell responses, patients with hematologic cancers and preserved CD8 T cells had a lower viral load and mortality. These data highlight the importance of CD8 T cells in acute COVID-19, particularly in the setting of impaired humoral immunity. Further, depletion of B cells with anti-CD20 therapy resulted in almost complete abrogation of SARS-CoV-2-specific IgG and IgM antibodies, but was not associated with increased mortality compared to other hematologic cancers, when adequate CD8 T cells were present. Finally, higher CD8 T cell counts were associated with improved overall survival in patients with hematologic cancers. Thus, CD8 T cells likely compensate for deficient humoral immunity and influence clinical recovery of COVID-19. These observations have important implications for cancer and COVID-19-directed treatments, immunosuppressive therapies, and for understanding the role of B and T cells in acute COVID-19.

Deep immune profiling of COVID-19 patients reveals distinct immunotypes with therapeutic implications.

Coronavirus disease 2019 (COVID-19) is currently a global pandemic, but human immune responses to the virus remain poorly understood. We used high-dimensional cytometry to analyze 125 COVID-19 patients and compare them with recovered and healthy individuals. Integrated analysis of ~200 immune and ~50 clinical features revealed activation of T cell and B cell subsets in a proportion of patients. A subgroup of patients had T cell activation characteristic of acute viral infection and plasmablast responses reaching >30% of circulating B cells. However, another subgroup had lymphocyte activation comparable with that in uninfected individuals. Stable versus dynamic immunological signatures were identified and linked to trajectories of disease severity change. Our analyses identified three immunotypes associated with poor clinical trajectories versus improving health. These immunotypes may have implications for the design of therapeutics and vaccines for COVID-19.

Alpha-Fetoprotein-Producing Lung Hepatoid Adenocarcinoma with Brain Metastasis Treated with S-1.

Lung hepatoid adenocarcinoma (HAC) is a rare primary lung carcinoma pathologically characterized by hepatocellular carcinoma-like tumor cells, the majority of which produce alpha-fetoprotein (AFP). The clinical prognosis of lung HAC is generally poor, and effective therapeutic regimens for inoperable or recurrent cases have not been established. Here, we report a case of AFP-producing lung HAC with brain metastasis with long-term disease control, treated with the 5-fluorouracil-derived regimen S-1. The patient was a 66-year-old male admitted to the hospital with alexia. Chest X-ray revealed a massive tumor in the left upper lobe, and a head CT scan revealed a metastasis in the left parietal lobe. The laboratory data showed a remarkably elevated AFP level (97,561 ng/mL). Pathological assessment of the resected brain tumor revealed HAC, which was compatible with the lung biopsies. Together with the absence of other metastatic lesions, a final diagnosis of primary lung HAC, stage IV T4N3M1b, was given. The patient first underwent non-small cell lung cancer chemotherapy regimens (carboplatin and paclitaxel as the first line, and pemetrexed as the second line), but had clinical progression. After third-line oral S-1 (tegafur/gimeracil/oteracil) administration, the serum AFP level significantly dropped and the patient achieved long-term disease control without relapse, surviving more than 19 months after disease presentation. The autopsy result was consistent with the diagnosis of primary lung HAC, and immunohistochemical staining was AFP+, glypican 3+, and spalt-like transcription factor 4+. Here, we report the case of a rare primary lung HAC with apparent disease control on S-1 therapy, together with a literature review.

CAR T-cell therapy is effective for CD19-dim B-lymphoblastic leukemia but is impacted by prior blinatumomab therapy.

Tisagenlecleucel, a chimeric antigen receptor (CAR) T-cell product targeting CD19 is approved for relapsed/refractory B-cell acute lymphoblastic leukemia (B-ALL). However, the impact of pretreatment variables, such as CD19 expression level, on leukemic blasts, the presence of CD19- subpopulations, and especially prior CD19-targeted therapy, on the response to CAR T-cell therapy has not been determined. We analyzed 166 patients treated with CAR T-cell therapy at our institution. Eleven patients did not achieve a minimal residual disease (MRD)- deep remission, whereas 67 patients had a recurrence after achieving a MRD- deep remission: 28 patients with CD19+ leukemia and 39 patients with CD19- leukemia. Return of CD19+ leukemia was associated with loss of CAR T-cell function, whereas CD19- leukemia was associated with continued CAR T-cell function. There were no significant differences in efficacy of CAR T cells in CD19-dim B-ALL, compared with CD19-normal or -bright B-ALL. Consistent with this, CAR T cells recognized and lysed cells with very low levels of CD19 expression in vitro. The presence of dim CD19 or rare CD19- events by flow cytometry did not predict nonresponse or recurrence after CAR T-cell therapy. However, prior therapy with the CD19-directed, bispecific T-cell engager blinatumomab was associated with a significantly higher rate of failure to achieve MRD- remission or subsequent loss of remission with antigen escape. Finally, immunophenotypic heterogeneity and lineage plasticity were independent of underlying clonotype and cytogenetic abnormalities.

Molecular Neuropathology in Practice: Clinical Profiling and Integrative Analysis of Molecular Alterations in Glioblastoma.

Molecular profiling of glioblastoma has revealed complex cytogenetic, epigenetic, and molecular abnormalities that are necessary for diagnosis, prognosis, and treatment. Our neuro-oncology group has developed a data-driven, institutional consensus guideline for efficient and optimal workup of glioblastomas based on our routine performance of molecular testing. We describe our institution's testing algorithm, assay development, and genetic findings in glioblastoma, to illustrate current practices and challenges in neuropathology related to molecular and genetic testing. We have found that coordination of test requisition, tissue handling, and incorporation of results into the final pathologic diagnosis by the neuropathologist improve patient care. Here, we present analysis of O6-methylguanine-DNA-methyltransferase promoter methylation and next-generation sequencing results of 189 patients, obtained utilizing our internal processes led by the neuropathology team. Our institutional pathway for neuropathologist-driven molecular testing has streamlined the management of glioblastoma samples for efficient return of results for incorporation of genomic data into the pathological diagnosis and optimal patient care.

Differences in Genomic Profiles and Outcomes Between Thoracic and Adrenal Neuroblastoma.

BACKGROUND: Neuroblastoma is a biologically and clinically heterogeneous disease. Based on recent studies demonstrating an association between the primary tumor site, prognosis, and commonly measured tumor biological features, we hypothesized that neuroblastomas arising in different sites would show distinct genomic features reflective of the developmental biology of the sympathicoadrenal nervous system. METHODS: We first compared genomic and epigenomic data of primary diagnostic neuroblastomas originating in the adrenal gland (n = 646) compared to thoracic sympathetic ganglia (n = 118). We also evaluated association of common germline variation with these primary sites in 1027 European-American neuroblastoma patients. RESULTS: We observed higher rates of MYCN amplification, chromosome 1q gain, and chromosome 11q deletion among adrenal tumors, which were highly predictive of functional RNA signatures. Surprisingly, thoracic neuroblastomas were more likely to harbor ALK driver mutations than adrenal cases among all cases (odds ratio = 1.89, 95% confidence interval = 1.04 to 3.43), and among cases without MYCN amplification (odds ratio = 2.86, 95% confidence interval = 1.48 to 5.49). Common germline single nucleotide polymorphisms (SNPs) in BARD1 (previously associated with high-risk neuroblastoma) were found to be strongly associated with predisposition for origin at adrenal, rather than thoracic, sites. CONCLUSIONS: Neuroblastomas arising in the adrenal gland are more likely to harbor structural DNA aberrations including MYCN amplification, whereas thoracic tumors show defects in mitotic checkpoints resulting in hyperdiploidy. Despite the general association of ALK mutations with high-risk disease, thoracic tumors are more likely to harbor gain-of-function ALK aberrations. Site of origin is likely reflective of stage of sympathetic nervous system development when malignant transformation occurs and is a surrogate for underlying tumor biology.

Cross-Cohort Analysis Identifies a TEAD4-MYCN Positive Feedback Loop as the Core Regulatory Element of High-Risk Neuroblastoma.

High-risk neuroblastomas show a paucity of recurrent somatic mutations at diagnosis. As a result, the molecular basis for this aggressive phenotype remains elusive. Recent progress in regulatory network analysis helped us elucidate disease-driving mechanisms downstream of genomic alterations, including recurrent chromosomal alterations. Our analysis identified three molecular subtypes of high-risk neuroblastomas, consistent with chromosomal alterations, and identified subtype-specific master regulator proteins that were conserved across independent cohorts. A 10-protein transcriptional module-centered around a TEAD4-MYCN positive feedback loop-emerged as the regulatory driver of the high-risk subtype associated with MYCN amplification. Silencing of either gene collapsed MYCN-amplified (MYCNAmp) neuroblastoma transcriptional hallmarks and abrogated viability in vitro and in vivo Consistently, TEAD4 emerged as a robust prognostic marker of poor survival, with activity independent of the canonical Hippo pathway transcriptional coactivators YAP and TAZ. These results suggest novel therapeutic strategies for the large subset of MYCN-deregulated neuroblastomas.Significance: Despite progress in understanding of neuroblastoma genetics, little progress has been made toward personalized treatment. Here, we present a framework to determine the downstream effectors of the genetic alterations sustaining neuroblastoma subtypes, which can be easily extended to other tumor types. We show the critical effect of disrupting a 10-protein module centered around a YAP/TAZ-independent TEAD4-MYCN positive feedback loop in MYCNAmp neuroblastomas, nominating TEAD4 as a novel candidate for therapeutic intervention. Cancer Discov; 8(5); 582-99. ©2018 AACR.This article is highlighted in the In This Issue feature, p. 517.

QuantumClone: clonal assessment of functional mutations in cancer based on a genotype-aware method for clonal reconstruction.

Motivation: In cancer, clonal evolution is assessed based on information coming from single nucleotide variants and copy number alterations. Nonetheless, existing methods often fail to accurately combine information from both sources to truthfully reconstruct clonal populations in a given tumor sample or in a set of tumor samples coming from the same patient. Moreover, previously published methods detect clones from a single set of variants. As a result, compromises have to be done between stringent variant filtering [reducing dispersion in variant allele frequency estimates (VAFs)] and using all biologically relevant variants. Results: We present a framework for defining cancer clones using most reliable variants of high depth of coverage and assigning functional mutations to the detected clones. The key element of our framework is QuantumClone, a method for variant clustering into clones based on VAFs, genotypes of corresponding regions and information about tumor purity. We validated QuantumClone and our framework on simulated data. We then applied our framework to whole genome sequencing data for 19 neuroblastoma trios each including constitutional, diagnosis and relapse samples. We confirmed an enrichment of damaging variants within such pathways as MAPK (mitogen-activated protein kinases), neuritogenesis, epithelial-mesenchymal transition, cell survival and DNA repair. Most pathways had more damaging variants in the expanding clones compared to shrinking ones, which can be explained by the increased total number of variants between these two populations. Functional mutational rate varied for ancestral clones and clones shrinking or expanding upon treatment, suggesting changes in clone selection mechanisms at different time points of tumor evolution. Availability and implementation: Source code and binaries of the QuantumClone R package are freely available for download at https://CRAN.R-project.org/package=QuantumClone. Contact: gudrun.schleiermacher@curie.fr or valentina.boeva@inserm.fr. Supplementary information: Supplementary data are available at Bioinformatics online.

LMO1 Synergizes with MYCN to Promote Neuroblastoma Initiation and Metastasis.

A genome-wide association study identified LMO1, which encodes an LIM-domain-only transcriptional cofactor, as a neuroblastoma susceptibility gene that functions as an oncogene in high-risk neuroblastoma. Here we show that dßh promoter-mediated expression of LMO1 in zebrafish synergizes with MYCN to increase the proliferation of hyperplastic sympathoadrenal precursor cells, leading to a reduced latency and increased penetrance of neuroblastomagenesis. The transgenic expression of LMO1 also promoted hematogenous dissemination and distant metastasis, which was linked to neuroblastoma cell invasion and migration, and elevated expression levels of genes affecting tumor cell-extracellular matrix interaction, including loxl3, itga2b, itga3, and itga5. Our results provide in vivo validation of LMO1 as an important oncogene that promotes neuroblastoma initiation, progression, and widespread metastatic dissemination.