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.

Common variants upstream of MLF1 at 3q25 and within CPZ at 4p16 associated with neuroblastoma.

Neuroblastoma is a cancer of the developing sympathetic nervous system that most commonly presents in young children and accounts for approximately 12% of pediatric oncology deaths. Here, we report on a genome-wide association study (GWAS) in a discovery cohort or 2,101 cases and 4,202 controls of European ancestry. We identify two new association signals at 3q25 and 4p16 that replicated robustly in multiple independent cohorts comprising 1,163 cases and 4,396 controls (3q25: rs6441201 combined P = 1.2x10-11, Odds Ratio 1.23, 95% CI:1.16-1.31; 4p16: rs3796727 combined P = 1.26x10-12, Odds Ratio 1.30, 95% CI: 1.21-1.40). The 4p16 signal maps within the carboxypeptidase Z (CPZ) gene. The 3q25 signal resides within the arginine/serine-rich coiled-coil 1 (RSRC1) gene and upstream of the myeloid leukemia factor 1 (MLF1) gene. Increased expression of MLF1 was observed in neuroblastoma cells homozygous for the rs6441201 risk allele (P = 0.02), and significant growth inhibition was observed upon depletion of MLF1 (P < 0.0001) in neuroblastoma cells. Taken together, we show that common DNA variants within CPZ at 4p16 and upstream of MLF1 at 3q25 influence neuroblastoma susceptibility and MLF1 likely plays an important role in neuroblastoma tumorigenesis.

Correction: Targeting the mTOR Complex by Everolimus in NRAS Mutant Neuroblastoma.

[This corrects the article DOI: 10.1371/journal.pone.0147682.].

Genetic predisposition to neuroblastoma mediated by a LMO1 super-enhancer polymorphism.

Neuroblastoma is a paediatric malignancy that typically arises in early childhood, and is derived from the developing sympathetic nervous system. Clinical phenotypes range from localized tumours with excellent outcomes to widely metastatic disease in which long-term survival is approximately 40% despite intensive therapy. A previous genome-wide association study identified common polymorphisms at the LMO1 gene locus that are highly associated with neuroblastoma susceptibility and oncogenic addiction to LMO1 in the tumour cells. Here we investigate the causal DNA variant at this locus and the mechanism by which it leads to neuroblastoma tumorigenesis. We first imputed all possible genotypes across the LMO1 locus and then mapped highly associated single nucleotide polymorphism (SNPs) to areas of chromatin accessibility, evolutionary conservation and transcription factor binding sites. We show that SNP rs2168101 G>T is the most highly associated variant (combined P = 7.47 × 10(-29), odds ratio 0.65, 95% confidence interval 0.60-0.70), and resides in a super-enhancer defined by extensive acetylation of histone H3 lysine 27 within the first intron of LMO1. The ancestral G allele that is associated with tumour formation resides in a conserved GATA transcription factor binding motif. We show that the newly evolved protective TATA allele is associated with decreased total LMO1 expression (P = 0.028) in neuroblastoma primary tumours, and ablates GATA3 binding (P < 0.0001). We demonstrate allelic imbalance favouring the G-containing strand in tumours heterozygous for this SNP, as demonstrated both by RNA sequencing (P < 0.0001) and reporter assays (P = 0.002). These findings indicate that a recently evolved polymorphism within a super-enhancer element in the first intron of LMO1 influences neuroblastoma susceptibility through differential GATA transcription factor binding and direct modulation of LMO1 expression in cis, and this leads to an oncogenic dependency in tumour cells.

A LIN28B-RAN-AURKA Signaling Network Promotes Neuroblastoma Tumorigenesis.

A more complete understanding of aberrant oncogenic signaling in neuroblastoma, a malignancy of the developing sympathetic nervous system, is paramount to improving patient outcomes. Recently, we identified LIN28B as an oncogenic driver in high-risk neuroblastoma. Here, we identify the oncogene RAN as a LIN28B target and show regional gain of chromosome 12q24 as an additional somatic alteration resulting in increased RAN expression. We show that LIN28B influences RAN expression by promoting RAN Binding Protein 2 expression and by directly binding RAN mRNA. Further, we demonstrate a convergence of LIN28B and RAN signaling on Aurora kinase A activity. Collectively, these findings demonstrate that LIN28B-RAN-AURKA signaling drives neuroblastoma oncogenesis, suggesting that this pathway may be amenable to therapeutic targeting.

CASC15-S Is a Tumor Suppressor lncRNA at the 6p22 Neuroblastoma Susceptibility Locus.

Chromosome 6p22 was identified recently as a neuroblastoma susceptibility locus, but its mechanistic contributions to tumorigenesis are as yet undefined. Here we report that the most highly significant single-nucleotide polymorphism (SNP) associations reside within CASC15, a long noncoding RNA that we define as a tumor suppressor at 6p22. Low-level expression of a short CASC15 isoform (CASC15-S) associated highly with advanced neuroblastoma and poor patient survival. In human neuroblastoma cells, attenuating CASC15-S increased cellular growth and migratory capacity. Gene expression analysis revealed downregulation of neuroblastoma-specific markers in cells with attenuated CASC15-S, with concomitant increases in cell adhesion and extracellular matrix transcripts. Altogether, our results point to CASC15-S as a mediator of neural growth and differentiation, which impacts neuroblastoma initiation and progression.

Relapsed neuroblastomas show frequent RAS-MAPK pathway mutations.

The majority of patients with neuroblastoma have tumors that initially respond to chemotherapy, but a large proportion will experience therapy-resistant relapses. The molecular basis of this aggressive phenotype is unknown. Whole-genome sequencing of 23 paired diagnostic and relapse neuroblastomas showed clonal evolution from the diagnostic tumor, with a median of 29 somatic mutations unique to the relapse sample. Eighteen of the 23 relapse tumors (78%) showed mutations predicted to activate the RAS-MAPK pathway. Seven of these events were detected only in the relapse tumor, whereas the others showed clonal enrichment. In neuroblastoma cell lines, we also detected a high frequency of activating mutations in the RAS-MAPK pathway (11/18; 61%), and these lesions predicted sensitivity to MEK inhibition in vitro and in vivo. Our findings provide a rationale for genetic characterization of relapse neuroblastomas and show that RAS-MAPK pathway mutations may function as a biomarker for new therapeutic approaches to refractory disease.

CODEX: a normalization and copy number variation detection method for whole exome sequencing.

High-throughput sequencing of DNA coding regions has become a common way of assaying genomic variation in the study of human diseases. Copy number variation (CNV) is an important type of genomic variation, but detecting and characterizing CNV from exome sequencing is challenging due to the high level of biases and artifacts. We propose CODEX, a normalization and CNV calling procedure for whole exome sequencing data. The Poisson latent factor model in CODEX includes terms that specifically remove biases due to GC content, exon capture and amplification efficiency, and latent systemic artifacts. CODEX also includes a Poisson likelihood-based recursive segmentation procedure that explicitly models the count-based exome sequencing data. CODEX is compared to existing methods on a population analysis of HapMap samples from the 1000 Genomes Project, and shown to be more accurate on three microarray-based validation data sets. We further evaluate performance on 222 neuroblastoma samples with matched normals and focus on a well-studied rare somatic CNV within the ATRX gene. We show that the cross-sample normalization procedure of CODEX removes more noise than normalizing the tumor against the matched normal and that the segmentation procedure performs well in detecting CNVs with nested structures.