Intercellular extrachromosomal DNA copy-number heterogeneity drives neuroblastoma cell state diversity.

Neuroblastoma exhibits significant inter- and intra-tumor genetic heterogeneity and varying clinical outcomes. Extrachromosomal DNAs (ecDNAs) may drive this heterogeneity by independently segregating during cell division, leading to rapid oncogene amplification. While ecDNA-mediated oncogene amplification is linked to poor prognosis in various cancers, the effects of ecDNA copy-number heterogeneity on intermediate phenotypes are poorly understood. Here, we leverage DNA and RNA sequencing from the same single cells in cell lines and neuroblastoma patients to investigate these effects. By analyzing ecDNA amplicon structures, we reveal extensive intercellular ecDNA copy-number heterogeneity. We also provide direct evidence of how this heterogeneity influences the expression of cargo genes, including MYCN and its downstream targets, and the overall transcriptional state of neuroblastoma cells. Our findings highlight the role of ecDNA copy number in promoting rapid adaptability of cellular states within tumors, underscoring the need for ecDNA-specific treatment strategies to address tumor formation and adaptation.

Reconstructing extrachromosomal DNA structural heterogeneity from long-read sequencing data using Decoil.

Circular extrachromosomal DNA (ecDNA) is a form of oncogene amplification found across cancer types and associated with poor outcome in patients. ecDNA can be structurally complex and can contain rearranged DNA sequences derived from multiple chromosome locations. As the structure of ecDNA can impact oncogene regulation and may indicate mechanisms of its formation, disentangling it at high resolution from sequencing data is essential. Even though methods have been developed to identify and reconstruct ecDNA in cancer genome sequencing, it remains challenging to resolve complex ecDNA structures, in particular amplicons with shared genomic footprints. We here introduce Decoil, a computational method that combines a breakpoint-graph approach with LASSO regression to reconstruct complex ecDNA and deconvolve co-occurring ecDNA elements with overlapping genomic footprints from long-read nanopore sequencing. Decoil outperforms de novo assembly and alignment-based methods in simulated long-read sequencing data for both simple and complex ecDNAs. Applying Decoil on whole-genome sequencing data uncovered different ecDNA topologies and explored ecDNA structure heterogeneity in neuroblastoma tumors and cell lines, indicating that this method may improve ecDNA structural analyses in cancer.

ML-driven segmentation of microvascular features during histological examination of tissue-engineered vascular grafts.

Introduction: The development of next-generation tissue-engineered medical devices such as tissue-engineered vascular grafts (TEVGs) is a leading trend in translational medicine. Microscopic examination is an indispensable part of animal experimentation, and histopathological analysis of regenerated tissue is crucial for assessing the outcomes of implanted medical devices. However, the objective quantification of regenerated tissues can be challenging due to their unusual and complex architecture. To address these challenges, research and development of advanced ML-driven tools for performing adequate histological analysis appears to be an extremely promising direction. Methods: We compiled a dataset of 104 representative whole slide images (WSIs) of TEVGs which were collected after a 6-month implantation into the sheep carotid artery. The histological examination aimed to analyze the patterns of vascular tissue regeneration in TEVGs in situ. Having performed an automated slicing of these WSIs by the Entropy Masker algorithm, we filtered and then manually annotated 1,401 patches to identify 9 histological features: arteriole lumen, arteriole media, arteriole adventitia, venule lumen, venule wall, capillary lumen, capillary wall, immune cells, and nerve trunks. To segment and quantify these features, we rigorously tuned and evaluated the performance of six deep learning models (U-Net, LinkNet, FPN, PSPNet, DeepLabV3, and MA-Net). Results: After rigorous hyperparameter optimization, all six deep learning models achieved mean Dice Similarity Coefficients (DSC) exceeding 0.823. Notably, FPN and PSPNet exhibited the fastest convergence rates. MA-Net stood out with the highest mean DSC of 0.875, demonstrating superior performance in arteriole segmentation. DeepLabV3 performed well in segmenting venous and capillary structures, while FPN exhibited proficiency in identifying immune cells and nerve trunks. An ensemble of these three models attained an average DSC of 0.889, surpassing their individual performances. Conclusion: This study showcases the potential of ML-driven segmentation in the analysis of histological images of tissue-engineered vascular grafts. Through the creation of a unique dataset and the optimization of deep neural network hyperparameters, we developed and validated an ensemble model, establishing an effective tool for detecting key histological features essential for understanding vascular tissue regeneration. These advances herald a significant improvement in ML-assisted workflows for tissue engineering research and development.

Genomic background sequences systematically outperform synthetic ones in de novo motif discovery for ChIP-seq data.

Efficient de novo motif discovery from the results of wide-genome mapping of transcription factor binding sites (ChIP-seq) is dependent on the choice of background nucleotide sequences. The foreground sequences (ChIP-seq peaks) represent not only specific motifs of target transcription factors, but also the motifs overrepresented throughout the genome, such as simple sequence repeats. We performed a massive comparison of the 'synthetic' and 'genomic' approaches to generate background sequences for de novo motif discovery. The 'synthetic' approach shuffled nucleotides in peaks, while in the 'genomic' approach selected sequences from the reference genome randomly or only from gene promoters according to the fraction of A/T nucleotides in each sequence. We compiled the benchmark collections of ChIP-seq datasets for mouse, human and Arabidopsis, and performed de novo motif discovery. We showed that the genomic approach has both more robust detection of the known motifs of target transcription factors and more stringent exclusion of the simple sequence repeats as possible non-specific motifs. The advantage of the genomic approach over the synthetic approach was greater in plants compared to mammals. We developed the AntiNoise web service (https://denovosea.icgbio.ru/antinoise/) that implements a genomic approach to extract genomic background sequences for twelve eukaryotic genomes.

Computational Reconstruction of the Transcription Factor Regulatory Network Induced by Auxin in Arabidopsis thaliana L.

In plant hormone signaling, transcription factor regulatory networks (TFRNs), which link the master transcription factors to the biological processes under their control, remain insufficiently characterized despite their crucial function. Here, we identify a TFRN involved in the response to the key plant hormone auxin and define its impact on auxin-driven biological processes. To reconstruct the TFRN, we developed a three-step procedure, which is based on the integrated analysis of differentially expressed gene lists and a representative collection of transcription factor binding profiles. Its implementation is available as a part of the CisCross web server. With the new method, we distinguished two transcription factor subnetworks. The first operates before auxin treatment and is switched off upon hormone application, the second is switched on by the hormone. Moreover, we characterized the functioning of the auxin-regulated TFRN in control of chlorophyll and lignin biosynthesis, abscisic acid signaling, and ribosome biogenesis.

Impact of primary care system interventions on healthcare worker satisfaction and intention to stay in the Philippines: a follow-up study.

BACKGROUND: The Philippines has a shortage and uneven distribution of healthcare workers (HCWs). Job satisfaction is an important element to HCW retention and attracting new HCWs into the health system. OBJECTIVE: This study measured HCWs' intent to stay and HCWs' satisfaction after implementation of multiple interventions intended to strengthen the primary care system, and determine factors significantly associated with HCWs' intent to stay. METHODOLOGY: This is a serial cross-sectional study in urban, rural and remote primary care sites in the Philippines. All physicians, nurses, midwives, dentists, community health workers and support staff were invited to participate. Baseline HCWs' intent to stay and satisfaction were obtained using a self-administered questionnaire prior to implementation of interventions. The same survey was again conducted in the years 2021 and 2022, corresponding to 5 and 6 years after initial implementation for the urban site, and 2 and 3 years for the rural and remote sites. We used multiple logistic regression to determine factors associated with intent to stay. RESULTS: There were 430 survey respondents (89.4% response rate) for year 2021, and 417 survey respondents (97.4% response rate) for year 2022. The urban and rural sites had significant increase in several HCW satisfaction domains, while the remote site had significant decrease in several HCW satisfaction domains. There was no significant difference in the intent to stay in the three sites. Factors that decreased intent to stay included length of employment, job involvement and employment as a nurse, while factors that increased intent to stay included job satisfaction, enjoyment and working in the urban site. CONCLUSION: HCW satisfaction improved in the urban site and rural site, while HCW satisfaction declined in the remote site. Intention to stay of primary care HCWs did not significantly change.

Observation of Seven Astrophysical Tau Neutrino Candidates with IceCube.

We report on a measurement of astrophysical tau neutrinos with 9.7 yr of IceCube data. Using convolutional neural networks trained on images derived from simulated events, seven candidate ν_{τ} events were found with visible energies ranging from roughly 20 TeV to 1 PeV and a median expected parent ν_{τ} energy of about 200 TeV. Considering backgrounds from astrophysical and atmospheric neutrinos, and muons from π^{±}/K^{±} decays in atmospheric air showers, we obtain a total estimated background of about 0.5 events, dominated by non-ν_{τ} astrophysical neutrinos. Thus, we rule out the absence of astrophysical ν_{τ} at the 5σ level. The measured astrophysical ν_{τ} flux is consistent with expectations based on previously published IceCube astrophysical neutrino flux measurements and neutrino oscillations.

Safari with an Electron Gun: Visualization of Protein and Membrane Interactions in Mitochondria in Natural Environment.

This paper presents new structural data about mitochondria using correlative light and electron microscopy (CLEM) and cryo-electron tomography. These state-of-the-art structural biology methods allow studying biological objects at nanometer scales under natural conditions. Non-invasiveness of these methods makes them comparable to observing animals in their natural environment on a safari. The paper highlights two areas of research that can only be accomplished using these methods. The study visualized location of the Aß42 amyloid aggregates in relation to mitochondria to test a hypothesis of development of mitochondrial dysfunction in Alzheimer's disease. The results showed that the Aß42 aggregates do not interact with mitochondria, although some of them are closely located. Therefore, the study demonstrated that mitochondrial dysfunction is not directly associated with the effects of aggregates on mitochondrial structure. Other processes should be considered as sources of mitochondrial dysfunction. Second unique area presented in this work is high-resolution visualization of the mitochondrial membranes and proteins in them. Analysis of the cryo-ET data reveals toroidal holes in the lamellar structures of cardiac mitochondrial cristae, where ATP synthases are located. The study proposes a new mechanism for sorting and clustering protein complexes in the membrane based on topology. According to this suggestion, position of the OXPHOS system proteins in the membrane is determined by its curvature. High-resolution tomography expands and complements existing ideas about the structural and functional organization of mitochondria. This makes it possible to study the previously inaccessible structural interactions of proteins with each other and with membranes in vivo.

An enteric ultrastructural surface atlas of the model insect Manducasexta.

The tobacco hornworm is a laboratory model that is particularly suitable for analyzing gut inflammation, but a physiological reference standard is currently unavailable. Here, we present a surface atlas of the healthy hornworm gut generated by scanning electron microscopy and nano-computed tomography. This comprehensive overview of the gut surface reveals morphological differences between the anterior, middle, and posterior midgut, allowing the screening of aberrant gut phenotypes while accommodating normal physiological variations. We estimated a total resorptive midgut surface of 0.42 m2 for L5d6 larvae, revealing its remarkable size. Our data will support allometric scaling and dose conversion from Manduca sexta to mammals in preclinical research, embracing the 3R principles. We also observed non-uniform gut colonization by enterococci, characterized by dense biofilms in the pyloric cone and downstream of the pylorus associated with pore and spine structures in the hindgut intima, indicating a putative immunosurveillance function in the lepidopteran hindgut.

Enhancing the Spin Hall Effect of Cylindrically Polarized Beams.

Two linked gear wheels in a micromachine can be simultaneously rotated in opposite directions by using a laser beam that has in its section areas the spin angular momentum (SAM) of the opposite sign. However, for instance, a cylindrical vector beam has zero SAM in the focus. We alter a cylindrical vector beam so as to generate areas in its focus where the SAM is of opposite signs. The first alteration is adding to the cylindrical vector beam a linearly polarized beam. Thus, we study superposition of two rotationally symmetric beams: those with cylindrical and linear polarization. We obtain an expression for the SAM and prove two of its properties. The first property is that changing superposition coefficients does not change the shape of the SAM density distribution, whereas the intensity changes. The second property is that maximal SAM density is achieved when both beams in the superposition have the same energy. The second perturbation is adding a spatial carrier frequency. We study the SAM density of a cylindrical vector beam with a spatial carrier frequency. Due to periodic modulation, upon propagation in space, such a beam is split into two beams, having left and right elliptic polarization. Thus, in the beam transverse section, areas with the spin of different signs are separated in space, which is a manifestation of the spin Hall effect. We demonstrate that such light beams can be generated by metasurfaces, with the transmittance depending periodically on one coordinate.

SkQ3 Exhibits the Most Pronounced Antioxidant Effect on Isolated Rat Liver Mitochondria and Yeast Cells.

Oxidative stress is involved in a wide range of age-related diseases. A critical role has been proposed for mitochondrial oxidative stress in initiating or promoting these pathologies and the potential for mitochondria-targeted antioxidants to fight them, making their search and testing a very urgent task. In this study, the mitochondria-targeted antioxidants SkQ1, SkQ3 and MitoQ were examined as they affected isolated rat liver mitochondria and yeast cells, comparing SkQ3 with clinically tested SkQ1 and MitoQ. At low concentrations, all three substances stimulated the oxidation of respiratory substrates in state 4 respiration (no ADP addition); at higher concentrations, they inhibited the ADP-triggered state 3 respiration and the uncoupled state, depolarized the inner mitochondrial membrane, contributed to the opening of the mPTP (mitochondrial permeability transition pore), did not specifically affect ATP synthase, and had a pronounced antioxidant effect. SkQ3 was the most active antioxidant, not possessing, unlike SkQ1 or MitoQ, prooxidant activity with increasing concentrations. In yeast cells, all three substances reduced prooxidant-induced intracellular oxidative stress and cell death and prevented and reversed mitochondrial fragmentation, with SkQ3 being the most efficient. These data allow us to consider SkQ3 as a promising potential therapeutic agent to mitigate pathologies associated with oxidative stress.

Passenger Gene Coamplifications Create Collateral Therapeutic Vulnerabilities in Cancer.

DNA amplifications in cancer do not only harbor oncogenes. We sought to determine whether passenger coamplifications could create collateral therapeutic vulnerabilities. Through an analysis of >3,000 cancer genomes followed by the interrogation of CRISPR-Cas9 loss-of-function screens across >700 cancer cell lines, we determined that passenger coamplifications are accompanied by distinct dependency profiles. In a proof-of-principle study, we demonstrate that the coamplification of the bona fide passenger gene DEAD-Box Helicase 1 (DDX1) creates an increased dependency on the mTOR pathway. Interaction proteomics identified tricarboxylic acid (TCA) cycle components as previously unrecognized DDX1 interaction partners. Live-cell metabolomics highlighted that this interaction could impair TCA activity, which in turn resulted in enhanced mTORC1 activity. Consequently, genetic and pharmacologic disruption of mTORC1 resulted in pronounced cell death in vitro and in vivo. Thus, structurally linked coamplification of a passenger gene and an oncogene can result in collateral vulnerabilities. SIGNIFICANCE: We demonstrate that coamplification of passenger genes, which were largely neglected in cancer biology in the past, can create distinct cancer dependencies. Because passenger coamplifications are frequent in cancer, this principle has the potential to expand target discovery in oncology. This article is featured in Selected Articles from This Issue, p. 384.

Association of 2-oxoacid dehydrogenase complexes with respirasomes in mitochondria.

In the present study, cryo-electron tomography was used to investigate the localization of 2-oxoacid dehydrogenase complexes (OADCs) in cardiac mitochondria and mitochondrial inner membrane samples. Two classes of ordered OADC inner cores with different symmetries were distinguished and their quaternary structures modeled. One class corresponds to pyruvate dehydrogenase complexes and the other to dehydrogenase complexes of α-ketoglutarate and branched-chain α-ketoacids. OADCs were shown to be localized in close proximity to membrane-embedded respirasomes, as observed both in densely packed lamellar cristae of cardiac mitochondria and in ruptured mitochondrial samples where the dense packing is absent. This suggests the specificity of the OADC-respirasome interaction, which allows localized NADH/NAD+ exchange between OADCs and complex I of the respiratory chain. The importance of this local coupling is based on OADCs being the link between respiration, glycolysis and amino acid metabolism. The coupling of these basic metabolic processes can vary in different tissues and conditions and may be involved in the development of various pathologies. The present study shows that this important and previously missing parameter of mitochondrial complex coupling can be successfully assessed using cryo-electron tomography.

Elimusertib has Antitumor Activity in Preclinical Patient-Derived Pediatric Solid Tumor Models.

The small-molecule inhibitor of ataxia telangiectasia and Rad3-related protein (ATR), elimusertib, is currently being tested clinically in various cancer entities in adults and children. Its preclinical antitumor activity in pediatric malignancies, however, is largely unknown. We here assessed the preclinical activity of elimusertib in 38 cell lines and 32 patient-derived xenograft (PDX) models derived from common pediatric solid tumor entities. Detailed in vitro and in vivo molecular characterization of the treated models enabled the evaluation of response biomarkers. Pronounced objective response rates were observed for elimusertib monotherapy in PDX, when treated with a regimen currently used in clinical trials. Strikingly, elimusertib showed stronger antitumor effects than some standard-of-care chemotherapies, particularly in alveolar rhabdomysarcoma PDX. Thus, elimusertib has strong preclinical antitumor activity in pediatric solid tumor models, which may translate to clinically meaningful responses in patients.

Evaluation of Mechanical Versus Manual Root Canal Preparation in Primary Molars-A Comparative In Vitro Study.

The endodontic treatment of primary teeth is to maintain the function of the tooth free of symptoms until its physiological exfoliation. A critical factor for success is how quickly and effectively the root canal preparation can be performed. Therefore, the aim of this comparative in vitro study was to analyze the efficiency of two mechanical root canal preparation systems FM (FlexMaster) and HF (HyFlex EDM) to manual KF (K-file) on extracted primary molars. A total of 45 teeth were divided into three groups (n = 15): KF (#15-35), FM (04#30) and HF (25/~ OneFile). Root canal preparation was performed, and the preparation time was measured. All root canals were non-destructively analyzed by micro-computed tomography in the cervical, middle and apical thirds before and after preparation with regard to the parameters of canal transport (in µm) and centering ratio (0-1). Statistical analysis was performed at a 5% significance level using non-parametric tests. HF caused the lowest canal transport in the apical third (p = 0.008). The centering ratio value of HF was significantly higher in the middle third of the root canals than in the other two groups (p < 0.01). The mean instrumentation time was significantly higher for KF (6.67 min) than for FM (4.69 min) and HF (4.03 min, p < 0.01). HF can be recommended for primary molar root canal treatment.

How Poor Is Aphyllophoroid Fungi Diversity in the Boreal Urban Greenhouses of Eastern Europe?

It is generally accepted that mycobiota diversity in urban greenhouses is poorer than in natural ecosystems, but our knowledge on this field of research is fragmentary. Here, we present the results of a long-term study of aphyllophoroid macrofungi (Basidiomycota) forming fruitbodies on non-native sub/tropical woody and herbaceous plants in the greenhouses of Saint Petersburg, Moscow, and Ekaterinburg botanical gardens located in the hemiboreal vegetation subzone of Eastern Europe. Over 20 years of research, fruitbodies of 58 species of aphyllophoroid fungi have been identified. Fungal species that developed on the wooden structures of greenhouses and building materials made of local wood are discussed separately. The list of fungi on non-native substrates is dominated by saprobes (93.1% of total list) as well as mycorrhizal with basidiomata on plants (8.6%). Phytopathogens have the lowest number (7.0%), and ¾ of species are widespread locally. Non-native plants are dominated by native fungal species (78.9%), while the percentage of non-native species is low (21.1%). In the three surveyed cities, the area of the studied greenhouses is 2.8 ha, and not a single species of fungi has been found twice on the same substrate. Half of the identified species are characterized by a single specimen (29 species/50.9%). Hymenochaete rheicolor was discovered in Russia for the first time and its known distribution is discussed. Only six (Antrodia gossypium, Hyphodontia arguta, Lyomyces sambuci, Peniophora cinerea, Ramariopsis kunzei, and Trechispora farinacea) local species (10.5%) were collected in all the three cities. The α-diversity of mycobiota (mean number of species per site, Shannon Index, and Menhinick Index) in the Ekaterinburg's greenhouses is 1.2-3.0 times lower compared to suburban forest parks and old-growth natural forests, while ß-diversity (Whittaker Index, Jaccard Index, and Morisita-Horn Index), on the contrary, is 2.1-7.7 times higher. With the plants' age, the probability of detecting fungi on them increases significantly. In greenhouses, phytopathogenic aphyllophoroid macrofungi are collected on woody plants only, but the probability of their development is not related to the plants' age.

Mutational topography reflects clinical neuroblastoma heterogeneity.

Neuroblastoma is a pediatric solid tumor characterized by strong clinical heterogeneity. Although clinical risk-defining genomic alterations exist in neuroblastomas, the mutational processes involved in their generation remain largely unclear. By examining the topography and mutational signatures derived from all variant classes, we identified co-occurring mutational footprints, which we termed mutational scenarios. We demonstrate that clinical neuroblastoma heterogeneity is associated with differences in the mutational processes driving these scenarios, linking risk-defining pathognomonic variants to distinct molecular processes. Whereas high-risk MYCN-amplified neuroblastomas were characterized by signs of replication slippage and stress, homologous recombination-associated signatures defined high-risk non-MYCN-amplified patients. Non-high-risk neuroblastomas were marked by footprints of chromosome mis-segregation and TOP1 mutational activity. Furthermore, analysis of subclonal mutations uncovered differential activity of these processes through neuroblastoma evolution. Thus, clinical heterogeneity of neuroblastoma patients can be linked to differences in the mutational processes that are active in their tumors.

Doxorubicin Changes the Spatial Organization of the Genome around Active Promoters.

In this study, we delve into the impact of genotoxic anticancer drug treatment on the chromatin structure of human cells, with a particular focus on the effects of doxorubicin. Using Hi-C, ChIP-seq, and RNA-seq, we explore the changes in chromatin architecture brought about by doxorubicin and ICRF193. Our results indicate that physiologically relevant doses of doxorubicin lead to a local reduction in Hi-C interactions in certain genomic regions that contain active promoters, with changes in chromatin architecture occurring independently of Top2 inhibition, cell cycle arrest, and differential gene expression. Inside the regions with decreased interactions, we detected redistribution of RAD21 around the peaks of H3K27 acetylation. Our study also revealed a common structural pattern in the regions with altered architecture, characterized by two large domains separated from each other. Additionally, doxorubicin was found to increase CTCF binding in H3K27 acetylated regions. Furthermore, we discovered that Top2-dependent chemotherapy causes changes in the distance decay of Hi-C contacts, which are driven by direct and indirect inhibitors. Our proposed model suggests that doxorubicin-induced DSBs cause cohesin redistribution, which leads to increased insulation on actively transcribed TAD boundaries. Our findings underscore the significant impact of genotoxic anticancer treatment on the chromatin structure of the human genome.

Spin Hall Effect in Paraxial Vectorial Light Beams with an Infinite Number of Polarization Singularities.

Elements of micromachines can be driven by light, including structured light with phase and/or polarization singularities. We investigate here a paraxial vector Gaussian beam with an infinite number of polarization singularities residing evenly on a straight line. The intensity distribution is derived analytically and the polarization singularities are shown to exist only in the initial plane and in the far field. The azimuthal angle of the polarization singularities is shown to increase in the far field by π/2. We obtain the longitudinal component of the spin angular momentum (SAM) density and show that it is independent of the azimuthal angle of the polarization singularities. Upon propagation in free space, an infinite number of C-points is generated, where polarization is circular. We show that the SAM density distribution has a shape of four spots, two with left and two with right elliptic polarization. The distance to the transverse plane with the maximal SAM density decreases with decreasing distance between the polarization singularities in the initial plane. Generating such alternating areas with positive and negative SAM density, despite linear polarization in the initial plane, manifests the optical spin Hall effect. Application areas of the obtained results include designing micromachines with optically driven elements.

Defining the landscape of circular RNAs in neuroblastoma unveils a global suppressive function of MYCN.

Circular RNAs (circRNAs) are a regulatory RNA class. While cancer-driving functions have been identified for single circRNAs, how they modulate gene expression in cancer is not well understood. We investigate circRNA expression in the pediatric malignancy, neuroblastoma, through deep whole-transcriptome sequencing in 104 primary neuroblastomas covering all risk groups. We demonstrate that MYCN amplification, which defines a subset of high-risk cases, causes globally suppressed circRNA biogenesis directly dependent on the DHX9 RNA helicase. We detect similar mechanisms in shaping circRNA expression in the pediatric cancer medulloblastoma implying a general MYCN effect. Comparisons to other cancers identify 25 circRNAs that are specifically upregulated in neuroblastoma, including circARID1A. Transcribed from the ARID1A tumor suppressor gene, circARID1A promotes cell growth and survival, mediated by direct interaction with the KHSRP RNA-binding protein. Our study highlights the importance of MYCN regulating circRNAs in cancer and identifies molecular mechanisms, which explain their contribution to neuroblastoma pathogenesis.