Neuroblastoma Susceptibility Research Articles

An evaluation of genetic predisposition to congenital anomalies and pediatric cancer supports KAT6B as a novel neuroblastoma susceptibility gene.

PurposeWhile congenital anomalies are among the strongest risk factors for developing pediatric cancer, the genetic underpinnings remain unclear. Therefore, we evaluated germline susceptibility in children with congenital anomalies and cancer. MethodsThrough the Genetic Overlap Between Anomalies and Cancer in Kids (GOBACK) Study, we recruited 47 participants with anomalies and cancer, along with their biological families. Genome sequencing was performed, focusing on single nucleotide variants, indels, and structural variants. Pathogenic or likely pathogenic (P/LP) variants were identified by American College of Medical Genetics and Genomics (ACMG) classification. ResultsWe identified P/LP variants in 23.4% (11 of 47) of participants. These variants encompassed: 1) Four genes associated with both anomalies and cancer (WT1, USP9X, PTPN1, and LZTR1), 2) Two established cancer predisposition genes (TP53 in 2 participants and PAX5), and 3) Four genes that are associated with anomalies (MMUT, FBN1, COL3A1, and KAT6B). We further investigated the role of KAT6B on neuroblastoma in a gene-based analysis from 409 neuroblastoma cases and 952 controls. This analysis demonstrated a significant enrichment of rare, predicted deleterious variants (p = 0.017), with odds ratios ranging from 2 to 4 based on the conditions we applied. ConclusionThis study demonstrates a molecular diagnostic yield of 23.4% in participants with both anomalies and cancer. Additionally, the findings further implicate the role of KAT6B as a novel neuroblastoma predisposition gene.

Integrative genomic analyses identify neuroblastoma risk genes involved in neuronal differentiation

Genome-Wide Association Studies (GWAS) have been decisive in elucidating the genetic predisposition of neuroblastoma (NB). The majority of genetic variants identified in GWAS are found in non-coding regions, suggesting that they can be causative of pathogenic dysregulations of gene expression. Nonetheless, pinpointing the potential causal genes within implicated genetic loci remains a major challenge. In this study, we integrated NB GWAS and expression Quantitative Trait Loci (eQTL) data from adrenal gland to identify candidate genes impacting NB susceptibility. We found that ZMYM1, CBL, GSKIP and WDR81 expression was dysregulated by NB predisposing variants. We further investigated the functional role of the identified genes through computational analysis of RNA sequencing (RNA-seq) data from single-cell and whole-tissue samples of NB, neural crest, and adrenal gland tissues, as well as through in vitro differentiation assays in NB cell cultures. Our results indicate that dysregulation of ZMYM1, CBL, GSKIP, WDR81 may lead to malignant transformation by affecting early and late stages of normal program of neuronal differentiation. Our findings enhance the understanding of how specific genes contribute to NB pathogenesis by highlighting their influence on neuronal differentiation and emphasizing the impact of genetic risk variants on the regulation of genes involved in critical biological processes.

Super-enhancer-driven IRF2BP2 enhances ALK activity and promotes neuroblastoma cell proliferation.

Super-enhancers (SEs) typically govern the expression of critical oncogenes and play a fundamental role in the initiation and progression of cancer. Focusing on genes that are abnormally regulated by SE in cancer may be a new strategy for understanding pathogenesis. In the context of this investigation, we have identified a previously unreported SE-driven gene IRF2BP2 in neuroblastoma (NB). The expression and prognostic value of IRF2BP2 were detected in public databases and clinical samples. The effect of IRF2BP2 on NB cell growth and apoptosis was evaluated through in vivo and in vitro functional loss experiments. The molecular mechanism of IRF2BP2 was investigated by the study of chromatin regulatory regions and transcriptome sequencing. The sustained high expression of IRF2BP2 results from the activation of a novel SE established by NB master transcription factors MYCN, MEIS2, and HAND2, and they form a new complex that regulates the gene network associated with the proliferation of NB cell populations. We also observed a significant enrichment of the AP-1 family at the binding sites of IRF2BP2. Remarkably, within NB cells, AP-1 plays a pivotal role in shaping the chromatin accessibility landscape, thereby exposing the binding site for IRF2BP2. This orchestrated action enables AP-1 and IRF2BP2 to collaboratively stimulate the expression of the NB susceptibility gene ALK, thereby upholding the highly proliferative phenotype characteristic of NB. Our findings indicate that SE-driven IRF2BP2 can bind to AP-1 to maintain the survival of tumor cells via regulating chromatin accessibility of the NB susceptibility gene ALK.

Neuroblastoma susceptibility and association of N7-methylguanosine modification gene polymorphisms: multi-center case-control study.

Neuroblastoma (NB) is a common extracranial solid malignancy in children. The N7-methylguanosine (m7G) modification gene METTL1/WDR4 polymorphisms may serve as promising molecular markers for identifying populations susceptible to NB. TaqMan probes was usded to genotype METTL1/WDR4 single nucleotide polymorphisms (SNPs) in 898 NB patients and 1734 healthy controls. A logistic regression model was utilized to calculate the odds ratio (OR) and 95% confidence interval (CI), evaluating the association between genotype polymorphisms and NB susceptibility. The analysis was also stratified by age, sex, tumor origin site, and clinical stage. Individual polymorphism of the METTL1/WDR4 gene investigated in this study did not show significant associations with NB susceptibility. However, combined genotype analysis revealed that carrying all 5 WDR4 protective genotypes was associated with a significantly lower NB risk compared to having 0-4 protective genotypes (AOR = 0.82, 95% CI = 0.69-0.96, P = 0.014). Further stratified analyses revealed that carrying 1-3 METTL1 risk genotypes, the WDR4 rs2156316 CG/GG genotype, the WDR4 rs2248490 CG/GG genotype, and having all five WDR4 protective genotypes were all significantly correlated with NB susceptibility in distinct subpopulations. In conclusion, our findings suggest significant associations between m7G modification gene METTL1/WDR4 SNPs and NB susceptibility in specific populations. Genetic variation in m7G modification gene is associated with susceptibility to NB. Single nucleotide polymorphisms in METTL1/WDR4 are associated with susceptibility to NB. Single nucleotide polymorphisms of METTL1/WDR4 can be used as a biomarker for screening NB susceptible populations.

Integrated re-analysis of transcriptomic and proteomic datasets reveals potential mechanisms for Zika viral-based oncolytic therapy in neuroblastoma.

Paediatric neuroblastoma and brain tumours account for a third of all childhood cancer-related mortality. High-risk neuroblastoma is highly aggressive and survival is poor despite intensive multi-modal therapies with significant toxicity. Novel therapies are desperately needed. The Zika virus (ZIKV) can access the nervous system and there is growing interest in employing ZIKV as a potential therapy against paediatric nervous system tumours, including neuroblastoma. Here, we perform extensive data mining, integration and re-analysis of ZIKV infection datasets to highlight molecular mechanisms that may govern the oncolytic response in neuroblastoma cells. We collate infection data of multiple neuroblastoma cell lines by different ZIKV strains from a body of published literature to inform the susceptibility of neuroblastoma to the ZIKV oncolytic response. Integrating published transcriptomics, interaction proteomics, dependency factor and compound datasets we propose the involvement of multiple host systems during ZIKV infection. Through data mining of published literature, we observed most paediatric neuroblastoma cell lines to be highly susceptible to ZIKV infection and propose the PRVABC59 ZIKV strain to be the most promising candidate for neuroblastoma oncolytic virotherapy. ZIKV induces TNF signalling, lipid metabolism, the Unfolded Protein Response (UPR), and downregulates cell cycle and DNA replication processes. ZIKV infection is dependent on sterol regulatory element binding protein (SREBP)-regulated lipid metabolism and three protein complexes; V-ATPase, ER Membrane Protein Complex (EMC) and mammalian translocon. We propose ZIKV non-structural protein 4B (NS4B) as a likely mediator of ZIKVs interaction with IRE1-mediated UPR, lipid metabolism and mammalian translocon. Our work provides a significant understanding of ZIKV infection in neuroblastoma cells, which will facilitate the progression of ZIKV-based oncolytic virotherapy through pre-clinical research and clinical trials.

Exploring the role of HLA variants in neuroblastoma susceptibility through whole exome sequencing.

Although a number of susceptibility loci for neuroblastoma (NB) have been identified by genome-wide association studies, it is still unclear whether variants in the HLA region contribute to NB susceptibility. In this study, we conducted a comprehensive genetic analysis of variants in the HLA region among 724 NB patients and 2863 matched controls from different cohorts. We exploited whole-exome sequencing data to accurately type HLA alleles with an ensemble approach on the results from three different typing tools, and carried out rigorous sample quality control to ensure a fine-scale ancestry matching. The frequencies of common HLA alleles were compared between cases and controls by logistic regression under additive and non-additive models. Population stratification was taken into account adjusting for ancestry-informative principal components. We detected significant HLA associations with NB. In particular, HLA-DQB1*05:02 (OR = 1.61; padj = 5.4 × 10-3) and HLA-DRB1*16:01 (OR = 1.60; padj = 2.3 × 10-2) alleles were associated to higher risk of developing NB. Conditional analysis highlighted the HLA-DQB1*05:02 allele and its residue Ser57 as key to this association. DQB1*05:02 allele was not associated to clinical features worse outcomes in the NB cohort. Nevertheless, a risk score derived from the allelic combinations of five HLA variants showed a substantial predictive value for patient survival (HR = 1.53; p = 0.032) that was independent from established NB prognostic factors. Our study leveraged powerful computational methods to explore WES data and HLA variants and to reveal complex genetic associations. Further studies are needed to validate the mechanisms of these interactions that contribute to the multifaceted pattern of factors underlying the disease initiation and progression.

Correlation of BARD1 gene polymorphisms with risk of neuroblastoma: a meta-analysis

BRCA1-associated RING domain protein 1 (BARD1) gene polymorphisms may be associated with neuroblastoma (NB) susceptibility. However, the results remain controversial. Relevant studies were identified by searching PubMed, Web of Science, Embase, China National Knowledge Infrastructure databases up to March 5, 2023. The strength of the association between BARD1 polymorphisms and susceptibility of NB was assessed by calculating odds ratios (ORs) and 95% confidence intervals (95% CIs) through the fixed- or random-effects model. Eight articles involving 12 studies were finally included. We found that rs6435862 T > G, rs3768716 A > G, rs17487792 C > T and rs7587476 C > T variant increase the risk of NB in allelic, dominant, recessive, homozygous and heterozygous genetic models, while rs7585356 G > A variant appeared protective against NB. When stratified by ethnicity, subgroup analysis indicated that the above association remained significant in Caucasian populations in all genetic models, except for rs7585356G > A polymorphism in Asians. In Asian populations, we found the similar results in the allelic and dominant model of rs6435862 T > G, rs3768716 A > G, rs17487792 C > T and rs7587476 C > T as in Caucasians, while there lacked a significant association in the other three model. In addition, rs7585356 G > A was not associated with an increased risk of NB in the Asian population. After Bonferroni correction, significant associations for rs7585356 G > A disappeared in both Asian and Caucasian populations, with no significant association found for rs7587476 in the allelic and dominant models among Asians. BARD1 polymorphisms might be significantly associated with NB susceptibility. It is crucial that these finding should be further confirmed through extensive and well-planned studies.

Genetic variants of m7G modification genes influence neuroblastoma susceptibility

ObjectiveNeuroblastoma is a life-threatening pediatric solid tumor whose etiology remains unclear. N7-methylguanosine (m7G) is one of the most important epigenetic modifications of RNA, which plays a crucial role in tumorigenesis. The m7G-mediated genes METTL1 and WDR4 also have been reported to be dysregulated in various cancers. However, the implications of METTL1 and WDR4 in neuroblastoma have not been clarified. MethodsGiven the oncogenic potential of m7G modification, we performed a case-control study to assess the association of METTL1 and WDR4 genes polymorphisms with neuroblastoma risk in a Chinese population consisting of 402 cases and 473 controls. Odds ratios (ORs) and 95 % confidence intervals (CIs) were applied to evaluate the associations between studied polymorphisms and neuroblastoma risk. The adjusted odds ratio (AOR) was adjusted for age and gender. ResultsOverall, four polymorphisms were significantly associated with neuroblastoma risk, including METTL1 rs2291617 (recessive model: adjusted OR = 1.59, 95 % CI = 1.08–2.34, P = 0.019), WDR4 rs2156316 (dominant model: adjusted OR = 0.74, 95 % CI = 0.57–0.97, P = 0.028), WDR4 rs6586250 (dominant model: adjusted OR = 0.59, 95 % CI = 0.42–0.84, P = 0.004) and WDR4 rs15736 (dominant model: adjusted OR = 0.60, 95 % CI = 0.42–0.85, P = 0.004). Stratified analysis showed stronger correlations between significant polymorphisms and neuroblastoma risk among subgroups divided by age, gender, tumor origin, and clinical stage. Furthermore, expression quantitative trait loci (eQTL) analysis revealed that significant polymorphisms were associated with the expression of the adjacent genes. ConclusionsOur study indicated that four polymorphisms in m7G-mediated genes contribute to neuroblastoma susceptibility in the eastern Chinese population. However, our findings should be verified further by large-scale and well-designed studies.

Functional TET2 gene polymorphisms increase the risk of neuroblastoma in Chinese children.

The 5-methylcytosine (m5C) is the key chemical modification in RNAs. As one of the demethylases in m5C, TET2 has been shown as a tumor suppressor. However, the impact of TET2 gene polymorphisms on neuroblastoma has not been elucidated. 402 neuroblastoma patients and 473 controls were genotyped for TET2 gene polymorphisms using the TaqMan method. The impact of TET2 gene polymorphisms on neuroblastoma susceptibility was determined using multivariate logistic regression analysis. We also adopted genotype-tissue expression database to explore the impact of TET2 gene polymorphisms on the expression of host and nearby genes. We used the R2 platform and Sangerbox tool to analyze the relationship between gene expression and neuroblastoma risk and prognosis through non-parametric testing and Kaplan-Meier analysis, respectively. We found the TET2 gene polymorphisms (rs10007915 G > C and rs7670522 A > C) and the combined 2-5 risk genotypes can significantly increase neuroblastoma risk. Stratification analysis showed that these significant associations were more prominent in certain subgroups. TET2 rs10007915 G > C and rs7670522 A > C are significantly associated with reduced expression of TET2 mRNA. Moreover, lower expression of TET2 gene is associated with high risk, MYCN amplification, and poor prognosis of neuroblastoma. The rs10007915 G > C and rs7670522 A > C are significantly related to the increased expression of inorganic pyrophosphatase 2 mRNA, and higher expression of PPA2 gene is associated with high risk, MYCN amplification, and poor prognosis of neuroblastomas. In summary, TET2 rs10007915 G > C and rs7670522 A > C significantly confer neuroblastoma susceptibility, and further research is needed to investigate the underlying mechanisms.

Integrated re-analysis of transcriptomic and proteomic datasets reveals potential mechanisms for Zika viral-based oncolytic therapy in neuroblastoma

Background Paediatric neuroblastoma and brain tumours account for a third of all childhood cancer-related mortality. High-risk neuroblastoma is highly aggressive and survival is poor despite intensive multi-modal therapies with significant toxicity. Novel therapies are desperately needed. The Zika virus (ZIKV) can access the nervous system and there is growing interest in employing ZIKV as a potential therapy against paediatric nervous system tumours, including neuroblastoma. Methods Here, we perform extensive data mining, integration and re-analysis of ZIKV infection datasets to highlight molecular mechanisms that may govern the oncolytic response in neuroblastoma cells. We collate infection data of multiple neuroblastoma cell lines by different ZIKV strains from a body of published literature to inform the susceptibility of neuroblastoma to the ZIKV oncolytic response. Integrating published transcriptomics, interaction proteomics, dependency factor and compound datasets we propose the involvement of multiple host systems during ZIKV infection. Results Through data mining of published literature, we observed most paediatric neuroblastoma cell lines to be highly susceptible to ZIKV infection and propose the PRVABC59 ZIKV strain to be the most promising candidate for neuroblastoma oncolytic virotherapy. ZIKV induces TNF signalling, lipid metabolism, the Unfolded Protein Response (UPR), and downregulates cell cycle and DNA replication processes. ZIKV infection is dependent on sterol regulatory element binding protein (SREBP)-regulated lipid metabolism and three protein complexes; V-ATPase, ER Membrane Protein Complex (EMC) and mammalian translocon. We propose ZIKV non-structural protein 4B (NS4B) as a likely mediator of ZIKVs interaction with IRE1-mediated UPR, lipid metabolism and mammalian translocon. Conclusions Our work provides a significant understanding of ZIKV infection in neuroblastoma cells, which will facilitate the progression of ZIKV-based oncolytic virotherapy through pre-clinical research and clinical trials.

MiR-100 rs1834306 a > G polymorphism decreases neuroblastoma risk in Chinese children.

Neuroblastoma is a common malignant tumor stemming from the sympathetic nervous system in children, which is often life-threatening. The genetics of neuroblastoma remains unclear. Studies have shown that miRNAs participate in the regulation of a broad spectrum of biological pathways. The abnormity in the miRNA is associated with the risk of various cancers, including neuroblastoma. However, research on the relationship of miRNA polymorphisms with neuroblastoma susceptibility is still in the initial stage. In this research, a retrospective case-control study was conducted to explore whether miR-100 rs1834306 A > G polymorphism is associated with neuroblastoma susceptibility. We enrolled 402 cases and 473 controls for the study. The logistic regression analysis was adopted to calculate odds ratios (ORs) and 95% confidence intervals (CIs) for the association between miR-100 rs1834306 A > G and neuroblastoma risk. Our results elucidated that the miR-100 rs1834306 A > G polymorphism was associated with the decreased risk of neuroblastoma (AG versus AA: adjusted OR = 0.72, 95% CI = 0.53-0.98, and P = 0.038). The subsequent stratified analysis further found that rs1834306 AG/GG genotype reduced the risk of neuroblastoma in the subgroup with tumors of the mediastinum origin (adjusted OR = 0.63, 95% CI = 0.41-0.95, and P = 0.029). In summary, miR-100 rs1834306 A > G polymorphism was shown to associate with decreased neuroblastoma risk in Chinese children, especially for neuroblastoma of mediastinum origin. This conclusion needs to be verified in additional large-size case-control studies.

Integrative transcriptomic and proteomic study of Zika viral infection reveals potential mechanisms for oncolytic therapy in neuroblastoma

Background: Paediatric neuroblastoma and brain tumours account for a third of all childhood cancer-related mortality. High-risk neuroblastoma is highly aggressive and survival is poor despite intensive multi-modal therapies with significant toxicity. Novel therapies are desperately needed. The Zika virus (ZIKV) is neurotropic and there is growing interest in employing ZIKV as a potential therapy against paediatric nervous system tumours, including neuroblastoma. Methods: Here, we perform extensive analysis of ZIKV infection studies to identify molecular mechanisms that may govern the oncolytic response in neuroblastoma cells. We summarise the neuroblastoma cell lines and ZIKV strains utilised and re-evaluate the infection data to deduce the susceptibility of neuroblastoma to the ZIKV oncolytic response. Integrating transcriptomics, interaction proteomics, dependency factor and compound datasets we show the involvement of multiple host systems during ZIKV infection. Results: We identified that most paediatric neuroblastoma cell lines are highly susceptible to ZIKV infection and that the PRVABC59 ZIKV strain is the most promising candidate for neuroblastoma oncolytic virotherapy. ZIKV induces TNF signalling, lipid metabolism, the Unfolded Protein Response (UPR), and downregulates cell cycle and DNA replication processes. ZIKV is dependent on sterol regulatory element binding protein (SREBP)-regulated lipid metabolism and three protein complexes; V-ATPase, ER Membrane Protein Complex (EMC) and mammalian translocon. We propose ZIKV non-structural protein 4B (NS4B) as a likely mediator of ZIKVs interaction with IRE1-mediated UPR, lipid metabolism and mammalian translocon. Conclusions: Our work provides a significant understanding of ZIKV infection in neuroblastoma cells, which will facilitate the progression of ZIKV-based oncolytic virotherapy through pre-clinical research and clinical trials.

Abstract 911: Butyrophilin 3A2 expression plays a critical role in phosphoantigen-mediated γδ T cell cytotoxicity of neuroblastoma cells

Abstract The clinical use of antibodies targeting the diaslganglioside GD2 improved survival and clearly identified a role for immunotherapy in the aggressive pediatric solid tumor neuroblastoma (NB). Even so, the 5-year survival rate for high risk patients remains 50%, warranting novel immunotherapeutic approaches. Given their MHC-independent cytotoxic activity towards transformed cells, γδ T cells are attractive tools for adoptive cellular immunotherapy. γδ T cell cytotoxicity is mediated through the expression of various activating receptors that contribute to multiple mechanisms of anti-tumor activity. Specifically, their expression of the γδ-TCR works in concert with tumor expression of members of the butyrophilin-3 (BTN3A) subfamily for cytotoxicity in response to phosphoantigens (pAg). In other cancers, stimulation of pAg production through zoledronate (ZOL) supplementation promotes γδ T cell recognition/killing. We have identified that γδ T cells are differentially cytotoxic towards NB cell lines, suggesting variations in expression of key players responsible for γδ T cell susceptibility. This work serves to understand the significance of BTN3A members in NB, and their impact on the susceptibility of NB to γδ-TCR mediated killing. It is hypothesized that the pAg-mediated response can be augmented by increasing intracellular pAg via use of ZOL or by manipulating BTN3A members. To test this, NB cells were pre-treated with either vehicle or ZOL for 24 hours and then co-incubated with expanded γδ T cells for 4 hours at various effector:target (E:T) ratios. Apoptotic death was then determined for NB cells using Annexin-V/7-AAD staining. Across 8 tested NB cell lines, Kelly, NB-1643, NGP, SK-N-AS, SH-SY5Y, NLF, and SMS-SAN cells were resistant to γδ T cells at a 10:1 E:T ratio when untreated with ZOL. IMR-5 demonstrated 12-18% lysis at a 5:1 E:T. When NB cells were pre-treated with ZOL, γδ T cell cytotoxicity was enhanced towards all models except SK-N-AS. We identified three NB cell lines with differential susceptibility to γδ T cell killing, SK-N-AS<NLF<IMR5. Using a high-risk NB gene expression dataset, expression of important BTN members was queried, revealing IMR5 to have BT3.2low expression, NLF to have BT3.2mod expression, and SK-N-AS to have BT3.2high expression. BT3.2 expression levels was confirmed via western blotting. It has been postulated that BT3.2 may be a decoy receptor that interacts with γδ T cells without intracellular signaling, given it lacks the intracellular pAg sensor. Due to its high expression in SK-N-AS, a lentiviral system was used to knockdown BT3.2. SK-N-A-S shBT3.2 models had 48-55% lysis at a 5:1 E:T when untreated and 65-72% lysis when pre-treated with ZOL. In summary, the anti-NB γδ T cell response can be augmented with ZOL pre-treatment. Additionally, the expression of BT3.2 may be responsible for differential sensitivity of NB cell lines to γδ T cells. Citation Format: Andrew Ho, Hunter C. Jonus, Kelly C. Goldsmith. Butyrophilin 3A2 expression plays a critical role in phosphoantigen-mediated γδ T cell cytotoxicity of neuroblastoma cells [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 911.

Variant rs8400 enhances ALKBH5 expression through disrupting miR-186 binding and promotes neuroblastoma progression.

AlkB homolog 5 (ALKBH5) has been proven to be closely related to tumors. However, the role and molecular mechanism of ALKBH5 in neuroblastomas have rarely been reported. The potential functional single-nucleotide polymorphisms (SNPs) in ALKBH5 were identified by National Center for Biotechnology Information (NCBI) dbSNP screening and SNPinfo software. TaqMan probes were used for genotyping. A multiple logistic regression model was used to evaluate the effects of different SNP loci on the risk of neuroblastoma. The expression of ALKBH5 in neuroblastoma was evaluated by Western blotting and immunohistochemistry (IHC). Cell counting kit-8 (CCK-8), plate colony formation and 5-ethynyl-2'-deoxyuridine (EdU) incorporation assays were used to evaluate cell proliferation. Wound healing and Transwell assays were used to compare cell migration and invasion. Thermodynamic modelling was performed to predict the ability of miRNAs to bind to ALKBH5 with the rs8400 G/A polymorphism. RNA sequencing, N6-methyladenosine (m6A) sequencing, m6A methylated RNA immunoprecipitation (MeRIP) and a luciferase assay were used to identify the targeting effect of ALKBH5 on SPP1. ALKBH5 was highly expressed in neuroblastoma. Knocking down ALKBH5 inhibited the proliferation, migration and invasion of cancer cells. miR-186-3p negatively regulates the expression of ALKBH5, and this ability is affected by the rs8400 polymorphism. When the G nucleotide was mutated to A, the ability of miR-186-3p to bind to the 3'-UTR of ALKBH5 decreased, resulting in upregulation of ALKBH5. SPP1 is the downstream target gene of the ALKBH5 oncogene. Knocking down SPP1 partially restored the inhibitory effect of ALKBH5 downregulation on neuroblastoma. Downregulation of ALKBH5 can improve the therapeutic efficacy of carboplatin and etoposide in neuroblastoma. We first found that the rs8400 G>A polymorphism in the m6A demethylase-encoding gene ALKBH5 increases neuroblastoma susceptibility and determines the related mechanisms. The aberrant regulation of ALKBH5 by miR-186-3p caused by this genetic variation in ALKBH5 promotes the occurrence and development of neuroblastoma through the ALKBH5-SPP1 axis.

CFDP1 is a neuroblastoma susceptibility gene that regulates transcription factors of the noradrenergic cell identity.

Pleiotropic genetic factors (e.g., DNA polymorphisms) may be involved in the initiation of neuroblastoma (NB) and coronary artery disease (CAD) given their common origin from defects in neural crest development. To discover novel NB susceptibility genes, we conducted a three-stage survey including a meta-analysis of NB and CAD genome-wide association data, prioritization of NB causal variants, and validation in an independent cohort of affected individuals-control subjects. The lead SNP, rs13337397 at the 16q23.1 locus, associated with both diseases in the meta-analysis and with NB in the validation study. All the SNPs in linkage disequilibrium with rs13337397 were annotated using the H3K27ac epigenetic marker of neural crest cells (NCC) and NB cell lines. Indeed, we identified the functional SNP rs13337017, mapping within an enhancer of NCCs and NB cell lines and showing long-range interactions with CFDP1 by Hi-C analysis. Luciferase assays indicated that the risk allele of rs13337017 increased CFDP1 expression in NB cell lines. Of note, CFDP1 high expression associated with unfavorable prognostic markers in an analysis including 498 NB transcriptomes. Moreover, depletion of CFDP1 markedly decreased viability and migration and increased apoptotic rates in NB cell lines. Finally, transcriptome and qPCR analyses revealed that the depletion of CFDP1 may affect noradrenergic neuron differentiation by downregulating master regulators of sympathetic noradrenergic identity, including PHOX2B, HAND2, and GATA3. Our data strongly suggest that CFDP1 acts as oncogene in NB. In addition, we provide evidence that genetic predisposition to NB can be mediated by the alteration of noradrenergic lineage-specific gene expression.

Inherited rare variants in homologous recombination and neurodevelopmental genes are associated with increased risk of neuroblastoma.

Neuroblastoma (NB) is the most common solid extracranial paediatric tumour. Genome-wide association studies have driven the discovery of common risk variants, but no large study has investigated the contribution of rare variants to NB susceptibility. Here, we conducted a whole-exome sequencing (WES) of 664 NB cases and 822 controls and used independent validation datasets to identify genes with rare risk variants and involved pathways. WES was performed at 50×depth and variants were jointly called in cases and controls. We developed two models to identify mutations with high clinical impact (P/LP model) and to discover less penetrant risk mutations affecting non-canonical cancer pathways (RPV model). We performed a gene-level collapsing test using Firth's logistic regression in 242 selected cancer predisposition genes (CPGs) and a gene-sets burden analysis of biologically-informed pathways. Twelve percent of patients carried P/LP variants in CPGs and showed a significant enrichment (P=2.3×10-4) compared to controls (6%). We identified P/LP variants in 45 CPGs enriched in homologous recombination (HR) pathway. The most P/LP enriched genes in NB were BRCA1, ALK and RAD51C. Additionally, we found higher RPV burden in gene-sets of neuron differentiation, neural tube development and synapse assembly, and in gene-sets associated with neurodevelopmental disorders (NDD). The high fraction of NB patients with P/LP variants indicates the need of genetic counselling. Furthermore, inherited rare variants predispose to NB development by affecting mechanisms related to HR and neurodevelopmental processes, and demonstrate that NDD genes are altered in NB at the germline level. Associazione Italiana per la Ricerca sul Cancro, Fondazione Italiana per la Lotta al Neuroblastoma, Associazione Oncologia Pediatrica e Neuroblastoma, Regione Campania, Associazione Giulio Adelfio onlus, and Italian Health Ministry.

NSUN2 gene rs13181449 C>T polymorphism reduces neuroblastoma risk

Neuroblastoma is the most common tumor in infants. RNA m5C modification regulates the survival, differentiation, and migration of cells affecting RNA function. However, the effects of the m5C modification methyltransferase gene NSUN2 polymorphism on neuroblastoma susceptibility have not been reported. TaqMan method was used to determine genotypes of four NSUN2 polymorphisms (rs4702373 C>T, rs13181449 C>T, rs166049 T>G, and rs8192120 A>C) in 402 patients with neuroblastoma and 473 cancer-free controls from Jiangsu province, China. Odds ratio (OR) and 95% confidence interval (CI) were used to evaluate the association of NSUN2 polymorphisms with neuroblastoma susceptibility. The association was also further assessed in subgroups stratified by age, sex, tumor origin, and stage. GTEx was used to analyze the effect of these polymorphisms on NSUN2 expression. We found the rs13181449 C>T was significantly associated with reduced neuroblastoma risk (CT vs. CC: adjusted OR = 0.68, 95% CI = 0.51–0.92, P = 0.012; CT/TT vs. CC: adjusted OR = 0.70, 95% CI = 0.53–0.92, P = 0.010). Compared with 0–2 protective genotypes, those with 3–4 protective genotypes could significantly reduce the neuroblastoma risk (adjusted OR = 0.68, 95% CI = 0.52 to 0.90, P = 0.006). Stratification analysis showed that the protective effect of rs13181449 polymorphism remained significant in children with age >18 months, boys, and those with early INSS stages. Moreover, children with more protective genotypes in the same subgroups also exhibited significantly reduced neuroblastoma risk. GTEx analysis showed that the rs13181449 T genotype was related with decreased NSUN2 gene expression. In conclusions, NSUN2 rs13181449 polymorphism is associated with decreased neuroblastoma risk, and the underlying mechanism in neuroblastoma needs further study.

Genetic Analysis in African American Children Supports Ancestry-Specific Neuroblastoma Susceptibility.

Neuroblastoma is rarer in African American (AA) children compared with American children of European descent. AA children affected with neuroblastoma, however, more frequently develop the high-risk form of the disease. We have genotyped an AA cohort of 629 neuroblastoma cases (254 high-risk) and 2,990 controls to investigate genetic susceptibility to neuroblastoma in AAs. We confirmed the known neuroblastoma susceptibility gene BARD1 at genome-wide significance in the subset of high-risk cases. We also estimated local admixture across the autosomal genome in the AA cases and controls and detected a signal at 4q31.22 where cases show an increase in European ancestry. A region at 17p13.1 showed increased African ancestry in the subgroup of high-risk cases with respect to intermediate- and low-risk cases. Using results from our published European American (EA) genome-wide association study (GWAS), we found that a polygenic score that included all independent SNPs showed a highly significant association (P value = 1.8 × 10-73) and explained 19% of disease risk variance in an independent EA cohort. In contrast, the best fit polygenic score (P value = 3.2 × 10-11) in AAs included only 22 independent SNPs with association P value < 2.75 × 10-6 in the EA GWAS, and explained 2% of neuroblastoma risk variance. The significance of the polygenic score dropped rapidly with inclusion of additional SNPs. These findings suggest that several common variants contribute to risk of neuroblastoma in an ancestry-specific fashion. This work supports the need for GWAS to be performed in populations of all races and ethnicities.

YTHDF2 Gene rs3738067 A>G Polymorphism Decreases Neuroblastoma Risk in Chinese Children: Evidence From an Eight-Center Case-Control Study.

Neuroblastoma is a primary malignancy mainly occurring in children. We have reported that polymorphisms of several N6-methyladenosine (m6A) RNA modification-related genes contributed to neuroblastoma risk in previous studies. YTHDF2, a “reader” of RNA m6A modification, is involved in cancer progression. Here, we estimated the association between a YTHDF2 gene rs3738067 A>G polymorphism and neuroblastoma susceptibility in 898 neuroblastoma patients and 1,734 healthy individuals from China. We found that the rs3738067 A>G could decrease neuroblastoma risk [AG vs. AA: adjusted odds ratio (OR) = 0.76, 95% confidence interval (CI) = 0.64–0.90, P = 0.002; AG/GG vs. AA: adjusted OR = 0.81, 95% CI = 0.69–0.95, P = 0.011). Besides, the rs3738067 AG/GG genotype was related to reduced neuroblastoma risk in the following subgroups: children aged 18 months and under, boys, patients with tumors originating from retroperitoneal, patients at clinical stage IV, and cases at clinical stages III plus IV. Importantly, false-positive report probability analysis proved our significant results worthy of close attention of. The expression quantitative trait locus analysis results revealed that the rs3738067 was associated with the expression of YTHDF2.

YTHDC1 gene polymorphisms and neuroblastoma susceptibility in Chinese children.

Neuroblastoma (NB) is the most common extracranial tumor in children. YTHDC1, a member of RNA methylation modification binding proteins, plays critical roles in tumor occurrence and metastasis. However, it is unclear whether YTHDC1 gene polymorphisms are related to NB susceptibility. Herein, we aimed to evaluate the association between YTHDC1 gene polymorphisms (rs2293596 T>C, rs2293595 T>C, rs3813832 T>C) and susceptibility of NB by logistic regression models. In this eight-center case-control study, 898 patients with NB and 1734 healthy controls were genotyped by TaqMan assay. The results showed that rs3813832 TC genotype could significantly reduce the susceptibility of NB compared with the TT genotype [adjusted odds ratio (AOR) = 0.81, 95% confidence interval (CI) = 0.68–0.96, P = 0.018]. Combined genotype analysis revealed that individuals with 3 protective genotypes had a prominently lower NB risk than those with 0-2 protective genotypes (AOR = 0.80, 95% CI = 0.68–0.94, P = 0.006). The stratified analysis also demonstrated the protective effect of rs3813832 TC/CC and 3 protective genotypes in certain subgroups. Further functional experiments revealed that YTHDC1 siRNA-554, targeting the area near the rs3813832 T>C polymorphism site, could observably inhibit the proliferation and migration of NB cells. In conclusion, our findings highlight the involvement of YTHDC1 gene and its genetic variants in the etiology of NB.