AutoGVP: a dockerized workflow integrating ClinVar and InterVar germline sequence variant classification.

SUMMARY: With the increasing rates of exome and whole genome sequencing, the ability to classify large sets of germline sequencing variants using up-to-date American College of Medical Genetics-Association for Molecular Pathology (ACMG-AMP) criteria is crucial. Here, we present Automated Germline Variant Pathogenicity (AutoGVP), a tool that integrates germline variant pathogenicity annotations from ClinVar and sequence variant classifications from a modified version of InterVar (PVS1 strength adjustments, removal of PP5/BP6). This tool facilitates large-scale, clinically focused classification of germline sequence variants in a research setting. AVAILABILITY AND IMPLEMENTATION: AutoGVP is an open source dockerized workflow implemented in R and freely available on GitHub at https://github.com/diskin-lab-chop/AutoGVP.

Pan-cancer genome and transcriptome analyses of 1,699 paediatric leukaemias and solid tumours.

Analysis of molecular aberrations across multiple cancer types, known as pan-cancer analysis, identifies commonalities and differences in key biological processes that are dysregulated in cancer cells from diverse lineages. Pan-cancer analyses have been performed for adult but not paediatric cancers, which commonly occur in developing mesodermic rather than adult epithelial tissues. Here we present a pan-cancer study of somatic alterations, including single nucleotide variants, small insertions or deletions, structural variations, copy number alterations, gene fusions and internal tandem duplications in 1,699 paediatric leukaemias and solid tumours across six histotypes, with whole-genome, whole-exome and transcriptome sequencing data processed under a uniform analytical framework. We report 142 driver genes in paediatric cancers, of which only 45% match those found in adult pan-cancer studies; copy number alterations and structural variants constituted the majority (62%) of events. Eleven genome-wide mutational signatures were identified, including one attributed to ultraviolet-light exposure in eight aneuploid leukaemias. Transcription of the mutant allele was detectable for 34% of protein-coding mutations, and 20% exhibited allele-specific expression. These data provide a comprehensive genomic architecture for paediatric cancers and emphasize the need for paediatric cancer-specific development of precision therapies.

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

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

Germline 16p11.2 Microdeletion Predisposes to Neuroblastoma.

Neuroblastoma is a cancer of the developing sympathetic nervous system. It is diagnosed in 600-700 children per year in the United States and accounts for 12% of pediatric cancer deaths. Despite recent advances in our understanding of this malignancy's complex genetic architecture, the contribution of rare germline variants remains undefined. Here, we conducted a genome-wide analysis of large (>500 kb), rare (<1%) germline copy number variants (CNVs) in two independent, multi-ethnic cohorts totaling 5,585 children with neuroblastoma and 23,505 cancer-free control children. We identified a 550-kb deletion on chromosome 16p11.2 significantly enriched in neuroblastoma cases (0.39% of cases and 0.03% of controls; p = 3.34 × 10-9). Notably, this CNV corresponds to a known microdeletion syndrome that affects approximately one in 3,000 children and confers risk for diverse developmental phenotypes including autism spectrum disorder and other neurodevelopmental disorders. The CNV had a substantial impact on neuroblastoma risk, with an odds ratio of 13.9 (95% confidence interval = 5.8-33.4). The association remained significant when we restricted our analysis to individuals of European ancestry in order to mitigate potential confounding by population stratification (0.42% of cases and 0.03% of controls; p = 4.10 × 10-8). We used whole-genome sequencing (WGS) to validate the deletion in paired germline and tumor DNA from 12 cases. Finally, WGS of four parent-child trios revealed that the deletion primarily arose de novo without maternal or paternal bias. This finding expands the clinical phenotypes associated with 16p11.2 microdeletion syndrome to include cancer, and it suggests that disruption of the 16p11.2 region may dysregulate neurodevelopmental pathways that influence both neurological phenotypes and neuroblastoma.

svpluscnv: analysis and visualization of complex structural variation data.

MOTIVATION: Despite widespread prevalence of somatic structural variations (SVs) across most tumor types, understanding of their molecular implications often remains poor. SVs are extremely heterogeneous in size and complexity, hindering the interpretation of their pathogenic role. Tools integrating large SV datasets across platforms are required to fully characterize the cancer's somatic landscape. RESULTS: svpluscnv R package is a swiss army knife for the integration and interpretation of orthogonal datasets including copy number variant segmentation profiles and sequencing-based structural variant calls. The package implements analysis and visualization tools to evaluate chromosomal instability and ploidy, identify genes harboring recurrent SVs and detects complex rearrangements such as chromothripsis and chromoplexia. Further, it allows systematic identification of hot-spot shattered genomic regions, showing reproducibility across alternative detection methods and datasets. AVAILABILITY AND IMPLEMENTATION: https://github.com/ccbiolab/svpluscnv. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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

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

Graphical Interpretation and Analysis of Proteins and their Ontologies (GiaPronto): A One-Click Graph Visualization Software for Proteomics Data Sets.

Here we present a free interactive web tool to process and visualize proteomics data sets with a single click. GiaPronto can process all proteomics quantification methods, i.e. label-free, SILAC and isobaric labeling, and analyze post-translational modifications (PTMs). The software performs normalization and statistics, assists determination of regulated proteins, biomarkers and Gene Ontology (GO) enrichment, and provides high resolution images and tables for further data analysis. We foresee that GiaPronto will become the most rapid and simple tool for assessing data quality and determining most relevant features of proteomic data sets. GiaPronto is available at giapronto.diskinlab.org.

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

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

Fine mapping of 2q35 high-risk neuroblastoma locus reveals independent functional risk variants and suggests full-length BARD1 as tumor-suppressor.

A previous genome-wide association study (GWAS) identified common variation at the BARD1 locus as being highly associated with susceptibility to high-risk neuroblastoma, but the mechanisms underlying this association have been not extensively investigated. Here, we performed a fine mapping analysis of BARD1 locus (2q35) using GWAS data from 556 high-risk neuroblastoma patients and 2,575 controls of European-American ancestry, and identified two independent genome-wide neuroblastoma-associated loci. Functional single-nucleotide polymorphism (SNP) prioritization identified two causative variants that independently contributed to neuroblastoma risk, and each replicated robustly in multiple independent cohorts comprising 445 high-risk cases and 3,170 controls (rs17489363: combined p = 1.07 × 10-31 , OR:1.79, 95% CI:1.62-1.98 and rs1048108: combined p = 7.27 × 10-14 , OR:0.65, 95% CI:0.58-0.73). Particularly, the T risk allele of rs17489363 in the canonical promoter region of full-length BARD1 altered binding site of the transcription factor HSF1 and correlated with low expression of full-length BARD1 mRNA and protein. Low-level expression of full-length BARD1 associated with advanced neuroblastoma. In human neuroblastoma cells, attenuating full-length BARD1 increased proliferation and invasion capacity. In conclusion, we have identified two potentially causative SNPs at the BARD1 locus associated with predisposition to high-risk neuroblastoma, and have shown that full-length BARD1 may act as tumor suppressor.

Genetic susceptibility to neuroblastoma: current knowledge and future directions.

Neuroblastoma, a malignancy of the developing peripheral nervous system that affects infants and young children, is a complex genetic disease. Over the past two decades, significant progress has been made toward understanding the genetic determinants that predispose to this often lethal childhood cancer. Approximately 1-2% of neuroblastomas are inherited in an autosomal dominant fashion and a combination of co-morbidity and linkage studies has led to the identification of germline mutations in PHOX2B and ALK as the major genetic contributors to this familial neuroblastoma subset. The genetic basis of "sporadic" neuroblastoma is being studied through a large genome-wide association study (GWAS). These efforts have led to the discovery of many common susceptibility alleles, each with modest effect size, associated with the development and progression of sporadic neuroblastoma. More recently, next-generation sequencing efforts have expanded the list of potential neuroblastoma-predisposing mutations to include rare germline variants with a predicted larger effect size. The evolving characterization of neuroblastoma's genetic basis has led to a deeper understanding of the molecular events driving tumorigenesis, more precise risk stratification and prognostics and novel therapeutic strategies. This review details the contemporary understanding of neuroblastoma's genetic predisposition, including recent advances and discusses ongoing efforts to address gaps in our knowledge regarding this malignancy's complex genetic underpinnings.

The functional variant rs34330 of CDKN1B is associated with risk of neuroblastoma.

The genetic aetiology of sporadic neuroblastoma is still largely unknown. We have identified diverse neuroblastoma susceptibility loci by genomewide association studies (GWASs); however, additional SNPs that likely contribute to neuroblastoma susceptibility prompted this investigation for identification of additional variants that are likely hidden among signals discarded by the multiple testing corrections used in the analysis of genomewide data. There is evidence suggesting the CDKN1B, coding for the cycle inhibitor p27Kip1, is involved in neuroblastoma. We thus assess whether genetic variants of CDKN1B are associated with neuroblastoma. We imputed all possible genotypes across CDKN1B locus on a discovery case series of 2101 neuroblastoma patients and 4202 genetically matched controls of European ancestry. The most significantly associated rs34330 was analysed in an independent Italian cohort of 311 cases and 709 controls. In vitro functional analysis was carried out in HEK293T and in neuroblastoma cell line SHEP-2, both transfected with pGL3-CDKN1B-CC or pGL3-CDKN1B-TT constructs. We identified an association of the rs34330 T allele (-79C/T) with the neuroblastoma risk (Pcombined = 0.002; OR = 1.17). The risk allele (T) of this single nucleotide polymorphism led to a lower transcription rate in cells transfected with a luciferase reporter driven by the polymorphic p27Kip1 promoter (P < 0.05). Three independent sets of neuroblastoma tumours carrying -79TT genotype showed a tendency towards lower CDKN1B mRNA levels. Our study shows that a functional variant, associated with a reduced CDKN1B gene transcription, influences neuroblastoma susceptibility.

MIBG avidity correlates with clinical features, tumor biology, and outcomes in neuroblastoma: A report from the Children's Oncology Group.

BACKGROUND: Prior studies suggest that neuroblastomas that do not accumulate metaiodobenzylguanidine (MIBG) on diagnostic imaging (MIBG non-avid) may have more favorable features compared with MIBG avid tumors. We compared clinical features, biologic features, and clinical outcomes between patients with MIBG nonavid and MIBG avid neuroblastoma. PROCEDURE: Patients had metastatic high- or intermediate-risk neuroblastoma and were treated on Children's Oncology Group protocols A3973 or A3961. Comparisons of clinical and biologic features according to MIBG avidity were made with chi-squared or Fisher exact tests. Event-free (EFS) and overall (OS) survival compared using log-rank tests and modeled using Cox models. RESULTS: Thirty of 343 patients (8.7%) had MIBG nonavid disease. Patients with nonavid tumors were less likely to have adrenal primary tumors (34.5 vs. 57.2%; P = 0.019), bone metastases (36.7 vs. 61.7%; P = 0.008), or positive urine catecholamines (66.7 vs. 91.0%; P < 0.001) compared with patients with MIBG avid tumors. Nonavid tumors were more likely to be MYCN amplified (53.8 vs. 32.6%; P = 0.030) and had lower norepinephrine transporter expression. Patients with MIBG nonavid disease had a 5-year EFS of 50.0% compared with 38.7% for patients with MIBG avid disease (P = 0.028). On multivariate testing in high-risk patients, MIBG avidity was the sole adverse prognostic factor for EFS identified (hazard ratio 1.77; 95% confidence interval 1.04-2.99; P = 0.034). CONCLUSIONS: Patients with MIBG nonavid neuroblastoma have lower rates of adrenal primary tumors, bone metastasis, and catecholamine secretion. Despite being more likely to have MYCN-amplified tumors, these patients have superior outcomes compared with patients with MIBG avid disease.

Integrative genomics identifies LMO1 as a neuroblastoma oncogene.

Neuroblastoma is a childhood cancer of the sympathetic nervous system that accounts for approximately 10% of all paediatric oncology deaths. To identify genetic risk factors for neuroblastoma, we performed a genome-wide association study (GWAS) on 2,251 patients and 6,097 control subjects of European ancestry from four case series. Here we report a significant association within LIM domain only 1 (LMO1) at 11p15.4 (rs110419, combined P = 5.2 × 10(-16), odds ratio of risk allele = 1.34 (95% confidence interval 1.25-1.44)). The signal was enriched in the subset of patients with the most aggressive form of the disease. LMO1 encodes a cysteine-rich transcriptional regulator, and its paralogues (LMO2, LMO3 and LMO4) have each been previously implicated in cancer. In parallel, we analysed genome-wide DNA copy number alterations in 701 primary tumours. We found that the LMO1 locus was aberrant in 12.4% through a duplication event, and that this event was associated with more advanced disease (P < 0.0001) and survival (P = 0.041). The germline single nucleotide polymorphism (SNP) risk alleles and somatic copy number gains were associated with increased LMO1 expression in neuroblastoma cell lines and primary tumours, consistent with a gain-of-function role in tumorigenesis. Short hairpin RNA (shRNA)-mediated depletion of LMO1 inhibited growth of neuroblastoma cells with high LMO1 expression, whereas forced expression of LMO1 in neuroblastoma cells with low LMO1 expression enhanced proliferation. These data show that common polymorphisms at the LMO1 locus are strongly associated with susceptibility to developing neuroblastoma, but also may influence the likelihood of further somatic alterations at this locus, leading to malignant progression.

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

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

A family-based study of gene variants and maternal folate and choline in neuroblastoma: a report from the Children's Oncology Group.

PURPOSE: Neuroblastoma is a childhood cancer of the sympathetic nervous system with embryonic origins. Previous epidemiologic studies suggest maternal vitamin supplementation during pregnancy reduces the risk of neuroblastoma. We hypothesized offspring and maternal genetic variants in folate-related and choline-related genes are associated with neuroblastoma and modify the effects of maternal intake of folate, choline, and folic acid. METHODS: The Neuroblastoma Epidemiology in North America (NENA) study recruited 563 affected children and their parents through the Children's Oncology Group's Childhood Cancer Research Network. We used questionnaires to ascertain pre-pregnancy supplementation and estimate usual maternal dietary intake of folate, choline, and folic acid. We genotyped 955 genetic variants related to folate or choline using DNA extracted from saliva samples and used a log-linear model to estimate both child and maternal risk ratios and stratum-specific risk ratios for gene-environment interactions. RESULTS: Overall, no maternal or offspring genotypic results met criteria for a false discovery rate (FDR) Q-value <0.2. Associations were also null for gene-environment interaction with pre-pregnancy vitamin supplementation, dietary folic acid, and folate. FDR-significant gene-choline interactions were found for offspring SNPs rs10489810 and rs9966612 located in MTHFD1L and TYMS, respectively, with maternal choline dietary intake dichotomized at the first quartile. CONCLUSION: These results suggest that variants related to one-carbon metabolism are not strongly associated with neuroblastoma. Choline-related variants may play a role; however, the functional consequences of the interacting variants are unknown and require independent replication.

Copy number variation at 1q21.1 associated with neuroblastoma.

Common copy number variations (CNVs) represent a significant source of genetic diversity, yet their influence on phenotypic variability, including disease susceptibility, remains poorly understood. To address this problem in human cancer, we performed a genome-wide association study of CNVs in the childhood cancer neuroblastoma, a disease in which single nucleotide polymorphism variations are known to influence susceptibility. We first genotyped 846 Caucasian neuroblastoma patients and 803 healthy Caucasian controls at approximately 550,000 single nucleotide polymorphisms, and performed a CNV-based test for association. We then replicated significant observations in two independent sample sets comprised of a total of 595 cases and 3,357 controls. Here we describe the identification of a common CNV at chromosome 1q21.1 associated with neuroblastoma in the discovery set, which was confirmed in both replication sets. This CNV was validated by quantitative polymerase chain reaction, fluorescent in situ hybridization and analysis of matched tumour specimens, and was shown to be heritable in an independent set of 713 cancer-free parent-offspring trios. We identified a previously unknown transcript within the CNV that showed high sequence similarity to several neuroblastoma breakpoint family (NBPF) genes and represents a new member of this gene family (NBPF23). This transcript was preferentially expressed in fetal brain and fetal sympathetic nervous tissues, and the expression level was strictly correlated with CNV state in neuroblastoma cells. These data demonstrate that inherited copy number variation at 1q21.1 is associated with neuroblastoma and implicate a previously unknown neuroblastoma breakpoint family gene in early tumorigenesis of this childhood cancer.

Phenotype restricted genome-wide association study using a gene-centric approach identifies three low-risk neuroblastoma susceptibility Loci.

Neuroblastoma is a malignant neoplasm of the developing sympathetic nervous system that is notable for its phenotypic diversity. High-risk patients typically have widely disseminated disease at diagnosis and a poor survival probability, but low-risk patients frequently have localized tumors that are almost always cured with little or no chemotherapy. Our genome-wide association study (GWAS) has identified common variants within FLJ22536, BARD1, and LMO1 as significantly associated with neuroblastoma and more robustly associated with high-risk disease. Here we show that a GWAS focused on low-risk cases identified SNPs within DUSP12 at 1q23.3 (P = 2.07 × 10⁻6), DDX4 and IL31RA both at 5q11.2 (P = 2.94 × 10⁻6 and 6.54 × 10⁻7 respectively), and HSD17B12 at 11p11.2 (P = 4.20 × 10⁻7) as being associated with the less aggressive form of the disease. These data demonstrate the importance of robust phenotypic data in GWAS analyses and identify additional susceptibility variants for neuroblastoma.

RNAi screen of the protein kinome identifies checkpoint kinase 1 (CHK1) as a therapeutic target in neuroblastoma.

Neuroblastoma is a childhood cancer that is often fatal despite intense multimodality therapy. In an effort to identify therapeutic targets for this disease, we performed a comprehensive loss-of-function screen of the protein kinome. Thirty kinases showed significant cellular cytotoxicity when depleted, with loss of the cell cycle checkpoint kinase 1 (CHK1/CHEK1) being the most potent. CHK1 mRNA expression was higher in MYC-Neuroblastoma-related (MYCN)-amplified (P < 0.0001) and high-risk (P = 0.03) tumors. Western blotting revealed that CHK1 was constitutively phosphorylated at the ataxia telangiectasia response kinase target site Ser345 and the autophosphorylation site Ser296 in neuroblastoma cell lines. This pattern was also seen in six of eight high-risk primary tumors but not in control nonneuroblastoma cell lines or in seven of eight low-risk primary tumors. Neuroblastoma cells were sensitive to the two CHK1 inhibitors SB21807 and TCS2312, with median IC(50) values of 564 nM and 548 nM, respectively. In contrast, the control lines had high micromolar IC(50) values, indicating a strong correlation between CHK1 phosphorylation and CHK1 inhibitor sensitivity (P = 0.0004). Furthermore, cell cycle analysis revealed that CHK1 inhibition in neuroblastoma cells caused apoptosis during S-phase, consistent with its role in replication fork progression. CHK1 inhibitor sensitivity correlated with total MYC(N) protein levels, and inducing MYCN in retinal pigmented epithelial cells resulted in CHK1 phosphorylation, which caused growth inhibition when inhibited. These data show the power of a functional RNAi screen to identify tractable therapeutical targets in neuroblastoma and support CHK1 inhibition strategies in this disease.