The Neuroblastoma Microenvironment, Heterogeneity and Immunotherapeutic Approaches.

Neuroblastoma is a peripheral nervous system tumor that almost exclusively occurs in young children. Although intensified treatment modalities have led to increased patient survival, the prognosis for patients with high-risk disease is still around 50%, signifying neuroblastoma as a leading cause of cancer-related deaths in children. Neuroblastoma is an embryonal tumor and is shaped by its origin from cells within the neural crest. Hence, neuroblastoma usually presents with a low mutational burden and is, in the majority of cases, driven by epigenetically deregulated transcription networks. The recent development of Omic techniques has given us detailed knowledge of neuroblastoma evolution, heterogeneity, and plasticity, as well as intra- and intercellular molecular communication networks within the neuroblastoma microenvironment. Here, we discuss the potential of these recent discoveries with emphasis on new treatment modalities, including immunotherapies which hold promise for better future treatment regimens.

Target Genes of c-MYC and MYCN with Prognostic Power in Neuroblastoma Exhibit Different Expressions during Sympathoadrenal Development.

Deregulation of the MYC family of transcription factors c-MYC (encoded by MYC), MYCN, and MYCL is prevalent in most human cancers, with an impact on tumor initiation and progression, as well as response to therapy. In neuroblastoma (NB), amplification of the MYCN oncogene and over-expression of MYC characterize approximately 40% and 10% of all high-risk NB cases, respectively. However, the mechanism and stage of neural crest development in which MYCN and c-MYC contribute to the onset and/or progression of NB are not yet fully understood. Here, we hypothesized that subtle differences in the expression of MYCN and/or c-MYC targets could more accurately stratify NB patients in different risk groups rather than using the expression of either MYC gene alone. We employed an integrative approach using the transcriptome of 498 NB patients from the SEQC cohort and previously defined c-MYC and MYCN target genes to model a multigene transcriptional risk score. Our findings demonstrate that defined sets of c-MYC and MYCN targets with significant prognostic value, effectively stratify NB patients into different groups with varying overall survival probabilities. In particular, patients exhibiting a high-risk signature score present unfavorable clinical parameters, including increased clinical risk, higher INSS stage, MYCN amplification, and disease progression. Notably, target genes with prognostic value differ between c-MYC and MYCN, exhibiting distinct expression patterns in the developing sympathoadrenal system. Genes associated with poor outcomes are mainly found in sympathoblasts rather than in chromaffin cells during the sympathoadrenal development.

Expression and activation of nuclear hormone receptors result in neuronal differentiation and favorable prognosis in neuroblastoma.

BACKGROUND: Neuroblastoma (NB), a childhood tumor derived from the sympathetic nervous system, presents with heterogeneous clinical behavior. While some tumors regress spontaneously without medical intervention, others are resistant to therapy, associated with an aggressive phenotype. MYCN-amplification, frequently occurring in high-risk NB, is correlated with an undifferentiated phenotype and poor prognosis. Differentiation induction has been proposed as a therapeutic approach for high-risk NB. We have previously shown that MYCN maintains an undifferentiated state via regulation of the miR-17 ~ 92 microRNA cluster, repressing the nuclear hormone receptors (NHRs) estrogen receptor alpha (ERα) and the glucocorticoid receptor (GR). METHODS: Cell viability was determined by WST-1. Expression of differentiation markers was analyzed by Western blot, RT-qPCR, and immunofluorescence analysis. Metabolic phenotypes were studied using Agilent Extracellular Flux Analyzer, and accumulation of lipid droplets by Nile Red staining. Expression of angiogenesis, proliferation, and neuronal differentiation markers, and tumor sections were assessed by immunohistochemistry. Gene expression from NB patient as well as adrenal gland cohorts were analyzed using GraphPad Prism software (v.8) and GSEA (v4.0.3), while pseudo-time progression on post-natal adrenal gland cells from single-nuclei transcriptome data was computed using scVelo. RESULTS: Here, we show that simultaneous activation of GR and ERα potentiated induction of neuronal differentiation, reduced NB cell viability in vitro, and decreased tumor burden in vivo. This was accompanied by a metabolic reprogramming manifested by changes in the glycolytic and mitochondrial functions and in lipid droplet accumulation. Activation of the retinoic acid receptor alpha (RARα) with all-trans retinoic acid (ATRA) further enhanced the differentiated phenotype as well as the metabolic switch. Single-cell nuclei transcriptome analysis of human adrenal glands indicated a sequential expression of ERα, GR, and RARα during development from progenitor to differentiated chromaffin cells. Further, in silico analysis revealed that patients with higher combined expression of GR, ERα, and RARα mRNA levels had elevated expression of neuronal differentiation markers and a favorable outcome. CONCLUSION: Together, our findings suggest that combination therapy involving activation of several NHRs could be a promising pharmacological approach for differentiation treatment of NB patients.

Impaired oxygen-sensitive regulation of mitochondrial biogenesis within the von Hippel-Lindau syndrome.

Mitochondria are the main consumers of oxygen within the cell. How mitochondria sense oxygen levels remains unknown. Here we show an oxygen-sensitive regulation of TFAM, an activator of mitochondrial transcription and replication, whose alteration is linked to tumours arising in the von Hippel-Lindau syndrome. TFAM is hydroxylated by EGLN3 and subsequently bound by the von Hippel-Lindau tumour-suppressor protein, which stabilizes TFAM by preventing mitochondrial proteolysis. Cells lacking wild-type VHL or in which EGLN3 is inactivated have reduced mitochondrial mass. Tumorigenic VHL variants leading to different clinical manifestations fail to bind hydroxylated TFAM. In contrast, cells harbouring the Chuvash polycythaemia VHLR200W mutation, involved in hypoxia-sensing disorders without tumour development, are capable of binding hydroxylated TFAM. Accordingly, VHL-related tumours, such as pheochromocytoma and renal cell carcinoma cells, display low mitochondrial content, suggesting that impaired mitochondrial biogenesis is linked to VHL tumorigenesis. Finally, inhibiting proteolysis by targeting LONP1 increases mitochondrial content in VHL-deficient cells and sensitizes therapy-resistant tumours to sorafenib treatment. Our results offer pharmacological avenues to sensitize therapy-resistant VHL tumours by focusing on the mitochondria.

Aberrant splicing in neuroblastoma generates RNA-fusion transcripts and provides vulnerability to spliceosome inhibitors.

The paucity of recurrent mutations has hampered efforts to understand and treat neuroblastoma. Alternative splicing and splicing-dependent RNA-fusions represent mechanisms able to increase the gene product repertoire but their role in neuroblastoma remains largely unexplored. Here we investigate the presence and possible roles of aberrant splicing and splicing-dependent RNA-fusion transcripts in neuroblastoma. In addition, we attend to establish whether the spliceosome can be targeted to treat neuroblastoma. Through analysis of RNA-sequenced neuroblastoma we show that elevated expression of splicing factors is a strong predictor of poor clinical outcome. Furthermore, we identified >900 primarily intrachromosomal fusions containing canonical splicing sites. Fusions included transcripts from well-known oncogenes, were enriched for proximal genes and in chromosomal regions commonly gained or lost in neuroblastoma. As a proof-of-principle that these fusions can generate altered gene products, we characterized a ZNF451-BAG2 fusion, producing a truncated BAG2-protein which inhibited retinoic acid induced differentiation. Spliceosome inhibition impeded neuroblastoma fusion expression, induced apoptosis and inhibited xenograft tumor growth. Our findings elucidate a splicing-dependent mechanism generating altered gene products in neuroblastoma and show that the spliceosome is a potential target for clinical intervention.

The long non-coding RNA Cerox1 is a post transcriptional regulator of mitochondrial complex I catalytic activity.

To generate energy efficiently, the cell is uniquely challenged to co-ordinate the abundance of electron transport chain protein subunits expressed from both nuclear and mitochondrial genomes. How an effective stoichiometry of this many constituent subunits is co-ordinated post-transcriptionally remains poorly understood. Here we show that Cerox1, an unusually abundant cytoplasmic long noncoding RNA (lncRNA), modulates the levels of mitochondrial complex I subunit transcripts in a manner that requires binding to microRNA-488-3p. Increased abundance of Cerox1 cooperatively elevates complex I subunit protein abundance and enzymatic activity, decreases reactive oxygen species production, and protects against the complex I inhibitor rotenone. Cerox1 function is conserved across placental mammals: human and mouse orthologues effectively modulate complex I enzymatic activity in mouse and human cells, respectively. Cerox1 is the first lncRNA demonstrated, to our knowledge, to regulate mitochondrial oxidative phosphorylation and, with miR-488-3p, represent novel targets for the modulation of complex I activity.

Genomic Variants Among Threatened Acropora Corals.

Genomic sequence data for non-model organisms are increasingly available requiring the development of efficient and reproducible workflows. Here, we develop the first genomic resources and reproducible workflows for two threatened members of the reef-building coral genus Acropora We generated genomic sequence data from multiple samples of the Caribbean A. cervicornis (staghorn coral) and A. palmata (elkhorn coral), and predicted millions of nucleotide variants among these two species and the Pacific A. digitifera A subset of predicted nucleotide variants were verified using restriction length polymorphism assays and proved useful in distinguishing the two Caribbean acroporids and the hybrid they form ("A. prolifera"). Nucleotide variants are freely available from the Galaxy server (usegalaxy.org), and can be analyzed there with computational tools and stored workflows that require only an internet browser. We describe these data and some of the analysis tools, concentrating on fixed differences between A. cervicornis and A. palmata In particular, we found that fixed amino acid differences between these two species were enriched in proteins associated with development, cellular stress response, and the host's interactions with associated microbes, for instance in the ABC transporters and superoxide dismutase. Identified candidate genes may underlie functional differences in how these threatened species respond to changing environments. Users can expand the presented analyses easily by adding genomic data from additional species, as they become available.

Unity in defence: honeybee workers exhibit conserved molecular responses to diverse pathogens.

BACKGROUND: Organisms typically face infection by diverse pathogens, and hosts are thought to have developed specific responses to each type of pathogen they encounter. The advent of transcriptomics now makes it possible to test this hypothesis and compare host gene expression responses to multiple pathogens at a genome-wide scale. Here, we performed a meta-analysis of multiple published and new transcriptomes using a newly developed bioinformatics approach that filters genes based on their expression profile across datasets. Thereby, we identified common and unique molecular responses of a model host species, the honey bee (Apis mellifera), to its major pathogens and parasites: the Microsporidia Nosema apis and Nosema ceranae, RNA viruses, and the ectoparasitic mite Varroa destructor, which transmits viruses. RESULTS: We identified a common suite of genes and conserved molecular pathways that respond to all investigated pathogens, a result that suggests a commonality in response mechanisms to diverse pathogens. We found that genes differentially expressed after infection exhibit a higher evolutionary rate than non-differentially expressed genes. Using our new bioinformatics approach, we unveiled additional pathogen-specific responses of honey bees; we found that apoptosis appeared to be an important response following microsporidian infection, while genes from the immune signalling pathways, Toll and Imd, were differentially expressed after Varroa/virus infection. Finally, we applied our bioinformatics approach and generated a gene co-expression network to identify highly connected (hub) genes that may represent important mediators and regulators of anti-pathogen responses. CONCLUSIONS: Our meta-analysis generated a comprehensive overview of the host metabolic and other biological processes that mediate interactions between insects and their pathogens. We identified key host genes and pathways that respond to phylogenetically diverse pathogens, representing an important source for future functional studies as well as offering new routes to identify or generate pathogen resilient honey bee stocks. The statistical and bioinformatics approaches that were developed for this study are broadly applicable to synthesize information across transcriptomic datasets. These approaches will likely have utility in addressing a variety of biological questions.

Giraffe genome sequence reveals clues to its unique morphology and physiology.

The origins of giraffe's imposing stature and associated cardiovascular adaptations are unknown. Okapi, which lacks these unique features, is giraffe's closest relative and provides a useful comparison, to identify genetic variation underlying giraffe's long neck and cardiovascular system. The genomes of giraffe and okapi were sequenced, and through comparative analyses genes and pathways were identified that exhibit unique genetic changes and likely contribute to giraffe's unique features. Some of these genes are in the HOX, NOTCH and FGF signalling pathways, which regulate both skeletal and cardiovascular development, suggesting that giraffe's stature and cardiovascular adaptations evolved in parallel through changes in a small number of genes. Mitochondrial metabolism and volatile fatty acids transport genes are also evolutionarily diverged in giraffe and may be related to its unusual diet that includes toxic plants. Unexpectedly, substantial evolutionary changes have occurred in giraffe and okapi in double-strand break repair and centrosome functions.

Elephantid Genomes Reveal the Molecular Bases of Woolly Mammoth Adaptations to the Arctic.

Woolly mammoths and living elephants are characterized by major phenotypic differences that have allowed them to live in very different environments. To identify the genetic changes that underlie the suite of woolly mammoth adaptations to extreme cold, we sequenced the nuclear genome from three Asian elephants and two woolly mammoths, and we identified and functionally annotated genetic changes unique to woolly mammoths. We found that genes with mammoth-specific amino acid changes are enriched in functions related to circadian biology, skin and hair development and physiology, lipid metabolism, adipose development and physiology, and temperature sensation. Finally, we resurrected and functionally tested the mammoth and ancestral elephant TRPV3 gene, which encodes a temperature-sensitive transient receptor potential (thermoTRP) channel involved in thermal sensation and hair growth, and we show that a single mammoth-specific amino acid substitution in an otherwise highly conserved region of the TRPV3 channel strongly affects its temperature sensitivity.

Genome-wide analysis of signatures of selection in populations of African honey bees (Apis mellifera) using new web-based tools.

BACKGROUND: With the development of inexpensive, high-throughput sequencing technologies, it has become feasible to examine questions related to population genetics and molecular evolution of non-model species in their ecological contexts on a genome-wide scale. Here, we employed a newly developed suite of integrated, web-based programs to examine population dynamics and signatures of selection across the genome using several well-established tests, including F ST, pN/pS, and McDonald-Kreitman. We applied these techniques to study populations of honey bees (Apis mellifera) in East Africa. In Kenya, there are several described A. mellifera subspecies, which are thought to be localized to distinct ecological regions. RESULTS: We performed whole genome sequencing of 11 worker honey bees from apiaries distributed throughout Kenya and identified 3.6 million putative single-nucleotide polymorphisms. The dense coverage allowed us to apply several computational procedures to study population structure and the evolutionary relationships among the populations, and to detect signs of adaptive evolution across the genome. While there is considerable gene flow among the sampled populations, there are clear distinctions between populations from the northern desert region and those from the temperate, savannah region. We identified several genes showing population genetic patterns consistent with positive selection within African bee populations, and between these populations and European A. mellifera or Asian Apis florea. CONCLUSIONS: These results lay the groundwork for future studies of adaptive ecological evolution in honey bees, and demonstrate the use of new, freely available web-based tools and workflows ( http://usegalaxy.org/r/kenyanbee ) that can be applied to any model system with genomic information.

Polar bears exhibit genome-wide signatures of bioenergetic adaptation to life in the arctic environment.

Polar bears (Ursus maritimus) face extremely cold temperatures and periods of fasting, which might result in more severe energetic challenges than those experienced by their sister species, the brown bear (U. arctos). We have examined the mitochondrial and nuclear genomes of polar and brown bears to investigate whether polar bears demonstrate lineage-specific signals of molecular adaptation in genes associated with cellular respiration/energy production. We observed increased evolutionary rates in the mitochondrial cytochrome c oxidase I gene in polar but not brown bears. An amino acid substitution occurred near the interaction site with a nuclear-encoded subunit of the cytochrome c oxidase complex and was predicted to lead to a functional change, although the significance of this remains unclear. The nuclear genomes of brown and polar bears demonstrate different adaptations related to cellular respiration. Analyses of the genomes of brown bears exhibited substitutions that may alter the function of proteins that regulate glucose uptake, which could be beneficial when feeding on carbohydrate-dominated diets during hyperphagia, followed by fasting during hibernation. In polar bears, genes demonstrating signatures of functional divergence and those potentially under positive selection were enriched in functions related to production of nitric oxide (NO), which can regulate energy production in several different ways. This suggests that polar bears may be able to fine-tune intracellular levels of NO as an adaptive response to control trade-offs between energy production in the form of adenosine triphosphate versus generation of heat (thermogenesis).

Aye-aye population genomic analyses highlight an important center of endemism in northern Madagascar.

We performed a population genomics study of the aye-aye, a highly specialized nocturnal lemur from Madagascar. Aye-ayes have low population densities and extensive range requirements that could make this flagship species particularly susceptible to extinction. Therefore, knowledge of genetic diversity and differentiation among aye-aye populations is critical for conservation planning. Such information may also advance our general understanding of Malagasy biogeography, as aye-ayes have the largest species distribution of any lemur. We generated and analyzed whole-genome sequence data for 12 aye-ayes from three regions of Madagascar (North, West, and East). We found that the North population is genetically distinct, with strong differentiation from other aye-ayes over relatively short geographic distances. For comparison, the average FST value between the North and East aye-aye populations--separated by only 248 km--is over 2.1-times greater than that observed between human Africans and Europeans. This finding is consistent with prior watershed- and climate-based hypotheses of a center of endemism in northern Madagascar. Taken together, these results suggest a strong and long-term biogeographical barrier to gene flow. Thus, the specific attention that should be directed toward preserving large, contiguous aye-aye habitats in northern Madagascar may also benefit the conservation of other distinct taxonomic units. To help facilitate future ecological- and conservation-motivated population genomic analyses by noncomputational biologists, the analytical toolkit used in this study is available on the Galaxy Web site.

Galaxy tools to study genome diversity.

BACKGROUND: Intra-species genetic variation can be used to investigate population structure, selection, and gene flow in non-model vertebrates; and due to the plummeting costs for genome sequencing, it is now possible for small labs to obtain full-genome variation data from their species of interest. However, those labs may not have easy access to, and familiarity with, computational tools to analyze those data. RESULTS: We have created a suite of tools for the Galaxy web server aimed at handling nucleotide and amino-acid polymorphisms discovered by full-genome sequencing of several individuals of the same species, or using a SNP genotyping microarray. In addition to providing user-friendly tools, a main goal is to make published analyses reproducible. While most of the examples discussed in this paper deal with nuclear-genome diversity in non-human vertebrates, we also illustrate the application of the tools to fungal genomes, human biomedical data, and mitochondrial sequences. CONCLUSIONS: This project illustrates that a small group can design, implement, test, document, and distribute a Galaxy tool collection to meet the needs of a particular community of biologists.

Polar and brown bear genomes reveal ancient admixture and demographic footprints of past climate change.

Polar bears (PBs) are superbly adapted to the extreme Arctic environment and have become emblematic of the threat to biodiversity from global climate change. Their divergence from the lower-latitude brown bear provides a textbook example of rapid evolution of distinct phenotypes. However, limited mitochondrial and nuclear DNA evidence conflicts in the timing of PB origin as well as placement of the species within versus sister to the brown bear lineage. We gathered extensive genomic sequence data from contemporary polar, brown, and American black bear samples, in addition to a 130,000- to 110,000-y old PB, to examine this problem from a genome-wide perspective. Nuclear DNA markers reflect a species tree consistent with expectation, showing polar and brown bears to be sister species. However, for the enigmatic brown bears native to Alaska's Alexander Archipelago, we estimate that not only their mitochondrial genome, but also 5-10% of their nuclear genome, is most closely related to PBs, indicating ancient admixture between the two species. Explicit admixture analyses are consistent with ancient splits among PBs, brown bears and black bears that were later followed by occasional admixture. We also provide paleodemographic estimates that suggest bear evolution has tracked key climate events, and that PB in particular experienced a prolonged and dramatic decline in its effective population size during the last ca. 500,000 years. We demonstrate that brown bears and PBs have had sufficiently independent evolutionary histories over the last 4-5 million years to leave imprints in the PB nuclear genome that likely are associated with ecological adaptation to the Arctic environment.