Comprehensive Multi-omics Analysis Reveals Mitochondrial Stress as a Central Biological Hub for Spaceflight Impact.

Spaceflight is known to impose changes on human physiology with unknown molecular etiologies. To reveal these causes, we used a multi-omics, systems biology analytical approach using biomedical profiles from fifty-nine astronauts and data from NASA's GeneLab derived from hundreds of samples flown in space to determine transcriptomic, proteomic, metabolomic, and epigenetic responses to spaceflight. Overall pathway analyses on the multi-omics datasets showed significant enrichment for mitochondrial processes, as well as innate immunity, chronic inflammation, cell cycle, circadian rhythm, and olfactory functions. Importantly, NASA's Twin Study provided a platform to confirm several of our principal findings. Evidence of altered mitochondrial function and DNA damage was also found in the urine and blood metabolic data compiled from the astronaut cohort and NASA Twin Study data, indicating mitochondrial stress as a consistent phenotype of spaceflight.

The ADP/ATP translocase drives mitophagy independent of nucleotide exchange.

Mitochondrial homeostasis depends on mitophagy, the programmed degradation of mitochondria. Only a few proteins are known to participate in mitophagy. Here we develop a multidimensional CRISPR-Cas9 genetic screen, using multiple mitophagy reporter systems and pro-mitophagy triggers, and identify numerous components of parkin-dependent mitophagy1. Unexpectedly, we find that the adenine nucleotide translocator (ANT) complex is required for mitophagy in several cell types. Whereas pharmacological inhibition of ANT-mediated ADP/ATP exchange promotes mitophagy, genetic ablation of ANT paradoxically suppresses mitophagy. Notably, ANT promotes mitophagy independently of its nucleotide translocase catalytic activity. Instead, the ANT complex is required for inhibition of the presequence translocase TIM23, which leads to stabilization of PINK1, in response to bioenergetic collapse. ANT modulates TIM23 indirectly via interaction with TIM44, which regulates peptide import through TIM232. Mice that lack ANT1 show blunted mitophagy and consequent profound accumulation of aberrant mitochondria. Disease-causing human mutations in ANT1 abrogate binding to TIM44 and TIM23 and inhibit mitophagy. Together, our findings show that ANT is an essential and fundamental mediator of mitophagy in health and disease.

Somatic structural variation targets neurodevelopmental genes and identifies SHANK2 as a tumor suppressor in neuroblastoma.

Neuroblastoma is a malignancy of the developing sympathetic nervous system that accounts for 12% of childhood cancer deaths. Like many childhood cancers, neuroblastoma shows a relative paucity of somatic single-nucleotide variants (SNVs) and small insertions and deletions (indels) compared to adult cancers. Here, we assessed the contribution of somatic structural variation (SV) in neuroblastoma using a combination of whole-genome sequencing (WGS) of tumor-normal pairs (n = 135) and single-nucleotide polymorphism (SNP) genotyping of primary tumors (n = 914). Our study design allowed for orthogonal validation and replication across platforms. SV frequency, type, and localization varied significantly among high-risk tumors. MYCN nonamplified high-risk tumors harbored an increased SV burden overall, including a significant excess of tandem duplication events across the genome. Genes disrupted by SV breakpoints were enriched in neuronal lineages and associated with phenotypes such as autism spectrum disorder (ASD). The postsynaptic adapter protein-coding gene, SHANK2, located on Chromosome 11q13, was disrupted by SVs in 14% of MYCN nonamplified high-risk tumors based on WGS and 10% in the SNP array cohort. Expression of SHANK2 was low across human-derived neuroblastoma cell lines and high-risk neuroblastoma tumors. Forced expression of SHANK2 in neuroblastoma cells resulted in significant growth inhibition (P = 2.6 × 10-2 to 3.4 × 10-5) and accelerated neuronal differentiation following treatment with all-trans retinoic acid (P = 3.1 × 10-13 to 2.4 × 10-30). These data further define the complex landscape of somatic structural variation in neuroblastoma and suggest that events leading to deregulation of neurodevelopmental processes, such as inactivation of SHANK2, are key mediators of tumorigenesis in this childhood cancer.

Long noncoding RNAs are spatially correlated with transcription factors and regulate lung development.

Long noncoding RNAs (lncRNAs) are thought to play important roles in regulating gene transcription, but few have well-defined expression patterns or known biological functions during mammalian development. Using a conservative pipeline to identify lncRNAs that have important biological functions, we identified 363 lncRNAs in the lung and foregut endoderm. Importantly, we show that these lncRNAs are spatially correlated with transcription factors across the genome. In-depth expression analyses of lncRNAs with genomic loci adjacent to the critical transcription factors Nkx2.1, Gata6, Foxa2 (forkhead box a2), and Foxf1 mimic the expression patterns of their protein-coding neighbor. Loss-of-function analysis demonstrates that two lncRNAs, LL18/NANCI (Nkx2.1-associated noncoding intergenic RNA) and LL34, play distinct roles in endoderm development by controlling expression of critical developmental transcription factors and pathways, including retinoic acid signaling. In particular, we show that LL18/NANCI acts upstream of Nkx2.1 and downstream from Wnt signaling to regulate lung endoderm gene expression. These studies reveal that lncRNAs play an important role in foregut and lung endoderm development by regulating multiple aspects of gene transcription, often through regulation of transcription factor expression.

A transcriptome-based classifier to determine molecular subtypes in medulloblastoma.

Medulloblastoma is a highly heterogeneous pediatric brain tumor with five molecular subtypes, Sonic Hedgehog TP53-mutant, Sonic Hedgehog TP53-wildtype, WNT, Group 3, and Group 4, defined by the World Health Organization. The current mechanism for classification into these molecular subtypes is through the use of immunostaining, methylation, and/or genetics. We surveyed the literature and identified a number of RNA-Seq and microarray datasets in order to develop, train, test, and validate a robust classifier to identify medulloblastoma molecular subtypes through the use of transcriptomic profiling data. We have developed a GPL-3 licensed R package and a Shiny Application to enable users to quickly and robustly classify medulloblastoma samples using transcriptomic data. The classifier utilizes a large composite microarray dataset (15 individual datasets), an individual microarray study, and an RNA-Seq dataset, using gene ratios instead of gene expression measures as features for the model. Discriminating features were identified using the limma R package and samples were classified using an unweighted mean of normalized scores. We utilized two training datasets and applied the classifier in 15 separate datasets. We observed a minimum accuracy of 85.71% in the smallest dataset and a maximum of 100% accuracy in four datasets with an overall median accuracy of 97.8% across the 15 datasets, with the majority of misclassification occurring between the heterogeneous Group 3 and Group 4 subtypes. We anticipate this medulloblastoma transcriptomic subtype classifier will be broadly applicable to the cancer research and clinical communities.

Hedgehog actively maintains adult lung quiescence and regulates repair and regeneration.

Postnatal tissue quiescence is thought to be a default state in the absence of a proliferative stimulus such as injury. Although previous studies have demonstrated that certain embryonic developmental programs are reactivated aberrantly in adult organs to drive repair and regeneration, it is not well understood how quiescence is maintained in organs such as the lung, which displays a remarkably low level of cellular turnover. Here we demonstrate that quiescence in the adult lung is an actively maintained state and is regulated by hedgehog signalling. Epithelial-specific deletion of sonic hedgehog (Shh) during postnatal homeostasis in the murine lung results in a proliferative expansion of the adjacent lung mesenchyme. Hedgehog signalling is initially downregulated during the acute phase of epithelial injury as the mesenchyme proliferates in response, but returns to baseline during injury resolution as quiescence is restored. Activation of hedgehog during acute epithelial injury attenuates the proliferative expansion of the lung mesenchyme, whereas inactivation of hedgehog signalling prevents the restoration of quiescence during injury resolution. Finally, we show that hedgehog also regulates epithelial quiescence and regeneration in response to injury via a mesenchymal feedback mechanism. These results demonstrate that epithelial-mesenchymal interactions coordinated by hedgehog actively maintain postnatal tissue homeostasis, and deregulation of hedgehog during injury leads to aberrant repair and regeneration in the lung.

annoFuse: an R Package to annotate, prioritize, and interactively explore putative oncogenic RNA fusions.

BACKGROUND: Gene fusion events are significant sources of somatic variation across adult and pediatric cancers and are some of the most clinically-effective therapeutic targets, yet low consensus of RNA-Seq fusion prediction algorithms makes therapeutic prioritization difficult. In addition, events such as polymerase read-throughs, mis-mapping due to gene homology, and fusions occurring in healthy normal tissue require informed filtering, making it difficult for researchers and clinicians to rapidly discern gene fusions that might be true underlying oncogenic drivers of a tumor and in some cases, appropriate targets for therapy. RESULTS: We developed annoFuse, an R package, and shinyFuse, a companion web application, to annotate, prioritize, and explore biologically-relevant expressed gene fusions, downstream of fusion calling. We validated annoFuse using a random cohort of TCGA RNA-Seq samples (N = 160) and achieved a 96% sensitivity for retention of high-confidence fusions (N = 603). annoFuse uses FusionAnnotator annotations to filter non-oncogenic and/or artifactual fusions. Then, fusions are prioritized if previously reported in TCGA and/or fusions containing gene partners that are known oncogenes, tumor suppressor genes, COSMIC genes, and/or transcription factors. We applied annoFuse to fusion calls from pediatric brain tumor RNA-Seq samples (N = 1028) provided as part of the Open Pediatric Brain Tumor Atlas (OpenPBTA) Project to determine recurrent fusions and recurrently-fused genes within different brain tumor histologies. annoFuse annotates protein domains using the PFAM database, assesses reciprocality, and annotates gene partners for kinase domain retention. As a standard function, reportFuse enables generation of a reproducible R Markdown report to summarize filtered fusions, visualize breakpoints and protein domains by transcript, and plot recurrent fusions within cohorts. Finally, we created shinyFuse for algorithm-agnostic interactive exploration and plotting of gene fusions. CONCLUSIONS: annoFuse provides standardized filtering and annotation for gene fusion calls from STAR-Fusion and Arriba by merging, filtering, and prioritizing putative oncogenic fusions across large cancer datasets, as demonstrated here with data from the OpenPBTA project. We are expanding the package to be widely-applicable to other fusion algorithms and expect annoFuse to provide researchers a method for rapidly evaluating, prioritizing, and translating fusion findings in patient tumors.

Diagnosing Cornelia de Lange syndrome and related neurodevelopmental disorders using RNA sequencing.

PURPOSE: Neurodevelopmental disorders represent a frequent indication for clinical exome sequencing. Fifty percent of cases, however, remain undiagnosed even upon exome reanalysis. Here we show RNA sequencing (RNA-seq) on human B-lymphoblastoid cell lines (LCL) is highly suitable for neurodevelopmental Mendelian gene testing and demonstrate the utility of this approach in suspected cases of Cornelia de Lange syndrome (CdLS). METHODS: Genotype-Tissue Expression project transcriptome data for LCL, blood, and brain were assessed for neurodevelopmental Mendelian gene expression. Detection of abnormal splicing and pathogenic variants in these genes was performed with a novel RNA-seq diagnostic pipeline and using a validation CdLS-LCL cohort (n = 10) and test cohort of patients who carry a clinical diagnosis of CdLS but negative genetic testing (n = 5). RESULTS: LCLs share isoform diversity of brain tissue for a large subset of neurodevelopmental genes and express 1.8-fold more of these genes compared with blood (LCL, n = 1706; whole blood, n = 917). This enables testing of more than 1000 genetic syndromes. The RNA-seq pipeline had 90% sensitivity for detecting pathogenic events and revealed novel diagnoses such as abnormal splice products in NIPBL and pathogenic coding variants in BRD4 and ANKRD11. CONCLUSION: The LCL transcriptome enables robust frontline and/or reflexive diagnostic testing for neurodevelopmental disorders.

Correcting transcription factor gene sets for copy number and promoter methylation variations.

Gene set analysis provides a method to generate statistical inferences across sets of linked genes, primarily using high-throughput expression data. Common gene sets include biological pathways, operons, and targets of transcriptional regulators. In higher eukaryotes, especially when dealing with diseases with strong genetic and epigenetic components such as cancer, copy number loss and gene silencing through promoter methylation can eliminate the possibility that a gene is transcribed. This, in turn, can adversely affect the estimation of transcription factor or pathway activity from a set of target genes, as some of the targets may not be responsive to transcriptional regulation. Here we introduce a simple filtering approach that removes genes from consideration if they show copy number loss or promoter methylation, and demonstrate the improvement in inference of transcription factor activity in a simulated dataset based on the background expression observed in normal head and neck tissue.

A proteogenomic surfaceome study identifies DLK1 as an immunotherapeutic target in neuroblastoma.

Cancer immunotherapies have produced remarkable results in B-cell malignancies; however, optimal cell surface targets for many solid cancers remain elusive. Here, we present an integrative proteomic, transcriptomic, and epigenomic analysis of tumor specimens along with normal tissues to identify biologically relevant cell surface proteins that can serve as immunotherapeutic targets for neuroblastoma, an often-fatal childhood cancer of the developing nervous system. We apply this approach to human-derived cell lines (N=9) and cell/patient-derived xenograft (N=12) models of neuroblastoma. Plasma membrane-enriched mass spectrometry identified 1,461 cell surface proteins in cell lines and 1,401 in xenograft models, respectively. Additional proteogenomic analyses revealed 60 high-confidence candidate immunotherapeutic targets and we prioritized Delta-like canonical notch ligand 1 (DLK1) for further study. High expression of DLK1 directly correlated with the presence of a super-enhancer spanning the DLK1 locus. Robust cell surface expression of DLK1 was validated by immunofluorescence, flow cytometry, and immunohistochemistry. Short hairpin RNA mediated silencing of DLK1 in neuroblastoma cells resulted in increased cellular differentiation. ADCT-701, a DLK1-targeting antibody-drug conjugate (ADC), showed potent and specific cytotoxicity in DLK1-expressing neuroblastoma xenograft models. Moreover, DLK1 is highly expressed in several adult cancer types, including adrenocortical carcinoma (ACC), pheochromocytoma/paraganglioma (PCPG), hepatoblastoma, and small cell lung cancer (SCLC), suggesting potential clinical benefit beyond neuroblastoma. Taken together, our study demonstrates the utility of comprehensive cancer surfaceome characterization and credentials DLK1 as an immunotherapeutic target. Highlights: Plasma membrane enriched proteomics defines surfaceome of neuroblastomaMulti-omic data integration prioritizes DLK1 as a candidate immunotherapeutic target in neuroblastoma and other cancersDLK1 expression is driven by a super-enhancer DLK1 silencing in neuroblastoma cells results in cellular differentiation ADCT-701, a DLK1-targeting antibody-drug conjugate, shows potent and specific cytotoxicity in DLK1-expressing neuroblastoma preclinical models.

OpenPBTA: The Open Pediatric Brain Tumor Atlas.

Pediatric brain and spinal cancers are collectively the leading disease-related cause of death in children; thus, we urgently need curative therapeutic strategies for these tumors. To accelerate such discoveries, the Children's Brain Tumor Network (CBTN) and Pacific Pediatric Neuro-Oncology Consortium (PNOC) created a systematic process for tumor biobanking, model generation, and sequencing with immediate access to harmonized data. We leverage these data to establish OpenPBTA, an open collaborative project with over 40 scalable analysis modules that genomically characterize 1,074 pediatric brain tumors. Transcriptomic classification reveals universal TP53 dysregulation in mismatch repair-deficient hypermutant high-grade gliomas and TP53 loss as a significant marker for poor overall survival in ependymomas and H3 K28-mutant diffuse midline gliomas. Already being actively applied to other pediatric cancers and PNOC molecular tumor board decision-making, OpenPBTA is an invaluable resource to the pediatric oncology community.

Evaluation of the DLL3-targeting antibody-drug conjugate rovalpituzumab tesirine in preclinical models of neuroblastoma.

Neuroblastomas have neuroendocrine features and often show similar gene expression patterns to small cell lung cancer including high expression of delta-like ligand 3 (DLL3). Here we determine the efficacy of rovalpituzumab tesirine (Rova-T), an antibody drug conjugated (ADC) with a pyrrolobenzodiazepine (PBD) dimer toxin targeting DLL3, in preclinical models of human neuroblastoma. We evaluated DLL3 expression in RNA sequencing data sets and performed immunohistochemistry (IHC) on neuroblastoma patient derived xenograft (PDX), human neuroblastoma primary tumor and normal childhood tissue microarrays (TMAs). We then evaluated the activity of Rova-T against 11 neuroblastoma PDX models using varying doses and schedules and compared anti-tumor activity to expression levels. DLL3 mRNA was differentially overexpressed in neuroblastoma at comparable levels to small cell lung cancer, as well as Wilms and rhabdoid tumors. DLL3 protein was robustly expressed across the neuroblastoma PDX array, but membranous staining was variable. The human neuroblastoma array, however, showed staining in only 44% of cases, whereas no significant staining was observed in the normal childhood tissue array. Rova-T showed a clear dose response effect across the 11 models tested, with a single dose inducing a complete or partial response in 3/11 and stable disease in another 3/11 models. No overt signs of toxicity were observed, and there was no treatment-related mortality. Strong membranous staining was necessary, but not sufficient, for anti-tumor activity. Rova-T has activity in a subset of neuroblastoma preclinical models, but heterogeneous expression in these models and the near absence of expression seen in human tumors suggests that any DLL3-targeting clinical trial should be only performed with a robust companion diagnostic to evaluate DLL3 expression for patient selection.

Targeted gene expression profiling of inverted papilloma and squamous cell carcinoma.

BACKGROUND: Inverted papilloma (IP) is a sinonasal tumor with a well-known potential for malignant transformation. The purpose of this study was to identify the genes and pathways associated with IP, with progression to carcinoma-in-situ and invasive carcinoma. METHODS: To determine genes and molecular pathways that may indicate progression and correlate with histologic changes, we analyzed six IP without dysplasia, five IP with carcinoma-in-situ, and 13 squamous cell carcinoma ex-IP by targeted sequencing. The HTG EdgeSeq Oncology Biomarker Panel coupled with next-generation sequencing was used to evaluate 2560 transcripts associated with solid tumors. RESULTS: Progressive upregulation of 11 genes were observed (CALD1, COL1A1, COL3A1, COL4A2, COL5A2, FN1, ITGA5, LGALS1, MMP11, SERPINH1, SPARC) in the order of invasive carcinoma > carcinoma-in-situ > IP without dysplasia. When compared with IP without dysplasia, more genes are differentially expressed in invasive carcinoma than carcinoma-in-situ samples (341 downregulated/333 upregulated vs. 195 downregulated/156 upregulated). Gene set enrichment analysis determined three gene sets in common between the cohorts (epithelial mesenchymal transition, extracellular matrix organization, and coagulation). CONCLUSIONS: Progressive upregulation of genes specific to IP malignant degeneration has significant clinical implications. This panel of 11 genes will improve concordance of histologic classification, which can directly impact treatment and patient outcomes. Additionally, future studies on larger tumor sets may observe upregulation in the gene panel that preceded histologic changes, which may be useful for further risk stratification.

MYC Levels Regulate Metastatic Heterogeneity in Pancreatic Adenocarcinoma.

The degree of metastatic disease varies widely among patients with cancer and affects clinical outcomes. However, the biological and functional differences that drive the extent of metastasis are poorly understood. We analyzed primary tumors and paired metastases using a multifluorescent lineage-labeled mouse model of pancreatic ductal adenocarcinoma (PDAC)-a tumor type in which most patients present with metastases. Genomic and transcriptomic analysis revealed an association between metastatic burden and gene amplification or transcriptional upregulation of MYC and its downstream targets. Functional experiments showed that MYC promotes metastasis by recruiting tumor-associated macrophages, leading to greater bloodstream intravasation. Consistent with these findings, metastatic progression in human PDAC was associated with activation of MYC signaling pathways and enrichment for MYC amplifications specifically in metastatic patients. Collectively, these results implicate MYC activity as a major determinant of metastatic burden in advanced PDAC. SIGNIFICANCE: Here, we investigate metastatic variation seen clinically in patients with PDAC and murine PDAC tumors and identify MYC as a major driver of this heterogeneity.This article is highlighted in the In This Issue feature, p. 275.

Epigenetic regulator BMI1 promotes alveolar rhabdomyosarcoma proliferation and constitutes a novel therapeutic target.

Rhabdomyosarcoma (RMS) is an aggressive pediatric soft tissue sarcoma. There are two main subtypes of RMS, alveolar rhabdomyosarcoma (ARMS) and embryonal rhabdomyosarcoma. ARMS typically encompasses fusion-positive rhabdomyosarcoma, which expresses either PAX3-FOXO1 or PAX7-FOXO1 fusion proteins. There are no targeted therapies for ARMS; however, recent studies have begun to illustrate the cooperation between epigenetic proteins and the PAX3-FOXO1 fusion, indicating that epigenetic proteins may serve as targets in ARMS. Here, we investigate the contribution of BMI1, given the established role of this epigenetic regulator in sustaining aggression in cancer. We determined that BMI1 is expressed across ARMS tumors, patient-derived xenografts, and cell lines. We depleted BMI1 using RNAi and inhibitors (PTC-209 and PTC-028) and found that this leads to a decrease in cell growth/increase in apoptosis in vitro, and delays tumor growth in vivo. Our data suggest that BMI1 inhibition activates the Hippo pathway via phosphorylation of LATS1/2 and subsequent reduction in YAP levels and YAP/TAZ target genes. These results identify BMI1 as a potential therapeutic vulnerability in ARMS and warrant further investigation of BMI1 in ARMS and other sarcomas.

Macrophages in SHH subgroup medulloblastoma display dynamic heterogeneity that varies with treatment modality.

Tumor-associated macrophages (TAMs) play an important role in tumor immunity and comprise of subsets that have distinct phenotype, function, and ontology. Transcriptomic analyses of human medulloblastoma, the most common malignant pediatric brain cancer, showed that medulloblastomas (MBs) with activated sonic hedgehog signaling (SHH-MB) have significantly more TAMs than other MB subtypes. Therefore, we examined MB-associated TAMs by single-cell RNA sequencing of autochthonous murine SHH-MB at steady state and under two distinct treatment modalities: molecular-targeted inhibitor and radiation. Our analyses reveal significant TAM heterogeneity, identify markers of ontologically distinct TAM subsets, and show the impact of brain microenvironment on the differentiation of tumor-infiltrating monocytes. TAM composition undergoes dramatic changes with treatment and differs significantly between molecular-targeted and radiation therapy. We identify an immunosuppressive monocyte-derived TAM subset that emerges with radiation therapy and demonstrate its role in regulating T cell and neutrophil infiltration in MB.

A GPC2 antibody-drug conjugate is efficacious against neuroblastoma and small-cell lung cancer via binding a conformational epitope.

Glypican 2 (GPC2) is a MYCN-regulated, differentially expressed cell-surface oncoprotein and target for immune-based therapies in neuroblastoma. Here, we build on GPC2's immunotherapeutic attributes by finding that it is also a highly expressed, MYCN-driven oncoprotein on small-cell lung cancers (SCLCs), with significantly enriched expression in both the SCLC and neuroblastoma stem cell compartment.By solving the crystal structure of the D3-GPC2-Fab/GPC2 complex at 3.3 Å resolution, we further illustrate that the GPC2-directed antibody-drug conjugate (ADC; D3-GPC2-PBD), that links a human GPC2 antibody (D3) to DNA-damaging pyrrolobenzodiazepine (PBD) dimers, binds a tumor-specific, conformation-dependent epitope of the core GPC2 extracellular domain. We then show that this ADC induces durable neuroblastoma and SCLC tumor regression via induction of DNA damage, apoptosis, and bystander cell killing, notably with no signs of ADC-induced in vivo toxicity. These studies provide preclinical data to support the clinical translation of ADCs targeting GPC2.

GAS7 Deficiency Promotes Metastasis in MYCN-Driven Neuroblastoma.

One of the greatest barriers to curative treatment of neuroblastoma is its frequent metastatic outgrowth prior to diagnosis, especially in cases driven by amplification of the MYCN oncogene. However, only a limited number of regulatory proteins that contribute to this complex MYCN-mediated process have been elucidated. Here we show that the growth arrest-specific 7 (GAS7) gene, located at chromosome band 17p13.1, is preferentially deleted in high-risk MYCN-driven neuroblastoma. GAS7 expression was also suppressed in MYCN-amplified neuroblastoma lacking 17p deletion. GAS7 deficiency led to accelerated metastasis in both zebrafish and mammalian models of neuroblastoma with overexpression or amplification of MYCN. Analysis of expression profiles and the ultrastructure of zebrafish neuroblastoma tumors with MYCN overexpression identified that GAS7 deficiency led to (i) downregulation of genes involved in cell-cell interaction, (ii) loss of contact among tumor cells as critical determinants of accelerated metastasis, and (iii) increased levels of MYCN protein. These results provide the first genetic evidence that GAS7 depletion is a critical early step in the cascade of events culminating in neuroblastoma metastasis in the context of MYCN overexpression. SIGNIFICANCE: Heterozygous deletion or MYCN-mediated repression of GAS7 in neuroblastoma releases an important brake on tumor cell dispersion and migration to distant sites, providing a novel mechanism underlying tumor metastasis in MYCN-driven neuroblastoma.See related commentary by Menard, p. 2815.

NASA GeneLab RNA-seq consensus pipeline: standardized processing of short-read RNA-seq data.

With the development of transcriptomic technologies, we are able to quantify precise changes in gene expression profiles from astronauts and other organisms exposed to spaceflight. Members of NASA GeneLab and GeneLab-associated analysis working groups (AWGs) have developed a consensus pipeline for analyzing short-read RNA-sequencing data from spaceflight-associated experiments. The pipeline includes quality control, read trimming, mapping, and gene quantification steps, culminating in the detection of differentially expressed genes. This data analysis pipeline and the results of its execution using data submitted to GeneLab are now all publicly available through the GeneLab database. We present here the full details and rationale for the construction of this pipeline in order to promote transparency, reproducibility, and reusability of pipeline data; to provide a template for data processing of future spaceflight-relevant datasets; and to encourage cross-analysis of data from other databases with the data available in GeneLab.

CAMKV Is a Candidate Immunotherapeutic Target in MYCN Amplified Neuroblastoma.

We developed a computational pipeline designed to use RNA sequencing (n = 136) and gene expression profiling (n = 250) data from neuroblastoma tumors to identify cell surface proteins predicted to be highly expressed in MYCN amplified neuroblastomas and with little or no expression in normal human tissues. We then performed ChIP-seq in the MYCN amplified cell lines KELLY, NB-1643, and NGP to identify gene promoters that are occupied by MYCN protein to define the intersection with the differentially-expressed gene list. We initially identified 116 putative immunotherapy targets with predicted transmembrane domains, with the most significant differentially-expressed of these being the calmodulin kinase-like vesicle-associated gene (CAMKV, p = 2 × 10-6). CAMKV encodes a protein that binds calmodulin in the presence of calcium, but lacks the kinase activity of other calmodulin kinase family members. We confirmed that CAMKV is selectively expressed in 7/7 MYCN amplified neuroblastoma cell lines and showed that the transcription of CAMKV is directly controlled by MYCN. From membrane fractionation and immunohistochemistry, we verified that CAMKV is membranous in MYCN amplified neuroblastoma cell lines and patient-derived xenografts. Finally, immunohistochemistry showed that CAMKV is not expressed on normal tissues outside of the central nervous system. Together, these data demonstrate that CAMKV is a differentially-expressed cell surface protein that is transcriptionally regulated by MYCN, making it a candidate for targeting with antibodies or antibody-drug conjugates that do not cross the blood brain barrier.