Influence of genetic factors on long-term treatment related neurocognitive complications, and on anxiety and depression in survivors of childhood acute lymphoblastic leukemia: The Petale study

Background A substantial number of survivors of childhood acute lymphoblastic leukemia suffer from treatment-related late adverse effects including neurocognitive impairment. While multiple studies have described neurocognitive outcomes in childhood acute lymphoblastic leukemia (ALL) survivors, relatively few have investigated their association with individual genetic constitution. Methods To further address this issue, genetic variants located in 99 genes relevant to the effects of anticancer drugs and in 360 genes implicated in nervous system function and predicted to affect protein function, were pooled from whole exome sequencing data of childhood ALL survivors (PETALE cohort) and analyzed for an association with neurocognitive complications, as well as with anxiety and depression. Variants that sustained correction for multiple testing were genotyped in entire cohort (n = 236) and analyzed with same outcomes. Results Common variants in MTR, PPARA, ABCC3, CALML5, CACNB2 and PCDHB10 genes were associated with deficits in neurocognitive tests performance, whereas a variant in SLCO1B1 and EPHA5 genes was associated with anxiety and depression. Majority of associations were modulated by intensity of treatment. Associated variants were further analyzed in an independent SJLIFE cohort of 545 ALL survivors. Two variants, rs1805087 in methionine synthase, MTR and rs58225473 in voltage-dependent calcium channel protein encoding gene, CACNB2 are of particular interest, since associations of borderline significance were found in replication cohort and remain significant in combined discovery and replication groups (OR = 1.5, 95% CI, 1–2.3; p = 0.04 and; OR = 3.7, 95% CI, 1.25–11; p = 0.01, respectively). Variant rs4149056 in SLCO1B1 gene also deserves further attention since previously shown to affect methotrexate clearance and short-term toxicity in ALL patients. Conclusions Current findings can help understanding of the influence of genetic component on long-term neurocognitive impairment. Further studies are needed to confirm whether identified variants may be useful in identifying survivors at increased risk of these complications.


Introduction
Acute lymphoblastic leukemia (ALL) is the most frequent childhood cancer [1] accounting for approximately 25% of all cases [2]. The five-year survival rate of childhood ALL is currently greater than 85% due to the optimization of multi-agent risk-adapted treatment strategies [2,3]. However, the exposure to specific chemotherapeutic agents and/or cranial radiation therapy during a susceptible period of child development results in late-adverse effects (LAEs) [4][5][6] including neurocognitive impairments [2]. Clinically significant deficits among ALL survivors are most commonly found in attention [7][8][9][10][11][12], working memory [13], processing speed [9,14,15] and executive functions, such as verbal fluency and cognitive flexibility [16]. Neurocognitive impairment in childhood ALL survivors persist for many years after treatment [17,18]. Large survey studies like the Childhood Cancer Survivors Study (CCSS) as well as other studies conducted in childhood ALL survivors [7,9,19] have demonstrated higher risk of depression, anxiety, behavioural difficulties, distress, as well as post-traumatic symptoms compared to siblings [20][21][22][23][24][25]. Longitudinal follow-up in long-term survivors have indicated that frequency of distress evolves over time, with more than 10% of survivors experiencing significant increase in depression as well as in anxiety [26].
The childhood ALL survivor population is increasing in size and lifespan, and this specific population needs an effective evaluation and targeted interventions [27]. Thus, better understanding of LAEs and factors contributing to their development is important to guide survivorship health surveillance and strategies to prevent or remediate treatment-related toxicities [2,3]. Here we assessed the role of genetic factors in neurocognitive impairments as well as in anxiety and depression by interrogating the relationship between above mentioned complications and genotypic profiling of 459 candidate genes obtained through whole exome sequencing (WES) of childhood ALL survivors.

Discovery cohort
The discovery cohort included 236 patients diagnosed and treated for childhood ALL according to Dana Farber Cancer Institute (DFCI) ALL 87-01 to 05-01 protocols at Sainte-Justine University Health Center (SJUHC), Montreal, (Quebec), Canada. The participants were recruited during 2013-2015 in the context of the PETALE study, a multidisciplinary research project with the goal to identify and comprehensively characterize associated predictive biomarkers of long-term treatment related complications in childhood ALL survivors [6]. Eligible participants were younger than 19 years old at diagnosis, at least 5 years after diagnosis of ALL and older than 12 years at evaluation, without history of relapse or refractory ALL or neurological condition or Down syndrome and had not received a hematopoietic stem cell transplant. The time from end of treatment to evaluation ranged from 5-26 year with a median of 13 years. The patients were classified to standard (SR) and high risk (HR) groups based on prognostic factors, including age, white blood cell count, immunophenotype, and central nervous system (CNS) status at diagnosis [28,29]. The frequency of patients assigned to SR and HR groups during the treatment was 46.6% and 53.4%, respectively. They were almost exclusively of reported French Canadian descent (>95%). The HapMap genotype reference data [30,31] was used for Principal component analysis (PCA) [30] to test and confirm predominant European ancestry (S1 Fig).

Replication cohort
The replication cohort consisted of 545 ALL survivors (274 males and 271 females) of European ancestry (based on genotype data) enrolled in the St. Jude Lifetime Cohort (SJLIFE) study and were evaluated using the same testing procedures as in PETALE cohort. Participants were younger than 19 years old at diagnosis and younger than 40 years old at SJLIFE evaluation, with no diagnosis of neurologic condition or Down syndrome, no history of relapse and had not received a hematopoietic stem cell transplant. The time from primary cancer diagnosis to the most recent date of neurocognitive evaluation ranged from 10.95-45.60 years with a median of 26.84 years. The risk group assignment during the treatment (SR and HR groups) was not available for this cohort.

Neuropsychological evaluation
A neurocognitive evaluation was performed using standardized testing procedures. It included three indices from two cognitive measures that reflect common impairments among childhood ALL survivors and are also good predictors of general neuropsychological outcomes [32]: Trail Making Test-Condition 4-Letter-Number Sequencing score and Verbal Fluency-Condition 1 -Letter fluency score from the Delis-Kaplan Executive Function System (D-KEFS) [33]; and Digit Span from the Wechsler Adult Intelligence Scale-Fourth Edition (WAIS-IV) [34] total score. Trail Making Test (D-KEFS) score is a measure that reflects processing speed, psychomotor speed, and cognitive flexibility [35]. Verbal Fluency (D-KEFS) score is a measure of phonological fluency in verbal modality [36]. Digit span (WAIS-IV) total score is a measure of verbal working memory [37]. All raw scores were converted to ageadjusted scaled scores based on population means [38]. The neurocognitive outcomes were transformed into dichotomous variables and studied accordingly. For each of these variables, scores lower than one and a half standard deviations below the mean of the normative dataset were indicative of impairment [39], all other scores were considered non-impaired. Inventory (BYI), a self-report instrument to document psychological status in children from 7 to 18 years old [41]. For older participants (�19 years) we used the Brief Symptom Inventory-18 (BSI-18 anxiety and depression score), an 18-item self-report questionnaire, assessing psychological distress in adults [40], previously also used in cohorts of young and older adult survivors of childhood cancer [42,43]. Internal consistency coefficients measured by Cronbach's alphas were all satisfactory, >.80 [44]. Age-adjusted scores one standard deviation above the population mean were considered as impaired.

Sequencing and quality control
Whole-exome sequencing (WES) was performed on germline DNA, extracted from peripheral blood samples from a subset of 191 participants of PETALE cohort, using standard protocols as described previously [6].

Association analyses
Two sets of candidate genes were selected-the first consisted of 99 genes implicated in the metabolism of methotrexate (MTX) and corticosteroids (CS) which are known to impact neurocognitive outcomes [54, 55] and the second consisted of 360 genes implicated in nervous system function, selected using the KEGG PATHWAY Database [56]. A total of 76 common variants (27 in MTX/CS pathway and 49 in nervous system function) and 1337 rare variants that satisfied all above filtering criteria were identified as functionally predicted and were used in association analyses. The analyses between common genetic variants and neurocognitive outcomes as well as with anxiety/depression were performed by the allelic chi-square or Fisher's exact test implemented in PLINK v.1.07 [57,58]. Analyses were performed in 191 sequenced patients and stratified by sex, risk groups with different treatment intensity, and treatment with chemotherapy alone or chemotherapy and cranial radiation because these factors have an established role in modulating neurocognitive outcomes [4,59]. The Benjamini-Hochberg procedure for false discovery rate (FDR) [60,61] was used to adjust for multiple testing with a cut-off value of < 5% considered statistically significant. Selective genotyping of top-ranking common SNPs (based additionally on Bonferroni p-value corrected for the number of variants tested, p<0.001 and p<0.0019 for the neural and MTX/CS pathways, respectively) was carried out on the Sequenom platform at the McGill University and Génome Québec Innovation Centre, Montreal, (Quebec), Canada, to confirm the results and extend the analysis to entire cohort (n = 236) including one hundred ninety-one patients analysed above (S1 Table). Associations of genotyped variants with the outcomes were assessed using chisquare or Fisher exact test in SPSS v.24.0.0.0 and appropriate genetic models, which were presented relative to the minor allele. Genotype-outcome association was represented as an odds ratio (OR) with a 95% confidence interval (CI). For the variants of MTX pathway for which the association showed similar trend in validation cohort, the modulation of the effect by cumulative drug dose was also analyzed. For that, cumulative drug dose was dichotomized to above and below the median and the association was analyzed in each subgroup. Additionally, logistic regression model was used in which main effect (genotype and drug dose) and interaction term were added. The detailed list of the studied polymorphisms (DOI: 10.6084/m9. figshare.8051573), as well as the summary statistics for all polymorphisms analyzed from the sequencing data beyond those already presented in the regular and supplemental tables (DOI: 10.6084/m9.figshare.8051825) are provided.
For rare variants associations, we used the SKAT-O test (Optimal Sequence Kernel Association Test) [62,63] implemented in SKAT package v.1.3.2.1 [64] with FDR < 5% considered statistically significant. Collapsing approach that combines several rare variants into a single variable [65,66], with iterative exclusion of each single variant, was additionally performed to allow weighting variant contributions to association signals. These analyses were performed as exploratory and associated variants were not further analyzed by genotyping.

Replication analysis
Genotype data for selected variants were obtained from a larger effort to sequence wholegenomes of over three thousand long-term survivors participating in the SJLIFE cohort. For this replication analysis, we restricted inclusion to 545 ALL survivors of European ancestry. Associations of selected variants with respective neurocognitive outcomes were examined using chi-square or Fisher's exact tests, as appropriate, implemented in PLINK 1.9 [57, 58].

Neurocognitive and emotional disturbances
The median age of ALL survivors at the time of evaluation was 21 years, with almost equal sex distribution, their demographics and clinical characteristics are presented in Table 1. The most prevalent deficit in neurocognitive test performance was noted for digit span (19.5%) followed by verbal fluency (18.6%) and trail making test (9.3%). Moderate-severe anxiety was noted in 10.1% survivors, whereas 11.5% of survivors were affected by moderate-severe depression, which was comparable to published normative groups on anxiety and depression [40-42].

Common variants
Among common variants implicated in nervous system function obtained from WES data, significant associations were detected for four of them (CALML5, CACNB2, PCDHB10 and EPHA5) either in all survivors or following stratification according to sex, risk groups or CRT (S2 Table). These variants were further analyzed by genotyping in the entire PETALE cohort and the association was confirmed for all of them ( Table 2). The analyses were performed for the same subgroups for which association was noted for WES data, and additionally in all participants. The neurocognitive deficit related to digit span task was associated in an additive manner with the minor allele of rs58225473 in CACNB2 gene either in all patients (p = 0.02), or those who received chemotherapy only (p = 0.004). Homozygotes for the minor C allele of CALML5 rs10904516 were at higher risk of having deficit in verbal fluency score, whereas the neurocognitive deficit related to trail making test was associated with the minor allele of rs2907323 in PCDHB10 gene, both potentiated in HR participants (p = 0.03 and p = 0.01 respectively). The carriers of the minor C allele of EPHA5 rs33932471 were at higher risk of both moderate-severe anxiety and depression, with the strongest effects seen in females (p = 0.02 and p = 0.003, respectively).
Among common variants implicated in MTX/CS pathway obtained from WES data, the significant associations were detected for 6 of them (MTR, PPARA, ABCC3, SHMT1 and High risk  Table). The variants in MTR, PPARA, ABCC3 and SLCO1B1 genes were further analyzed by genotyping in entire PETALE cohort and the association was confirmed for all of them ( Table 3). The association between deficit in verbal fluency score and GG genotype of MTR rs1805087 was seen for all survivors (p = 0.01) and male participants (p = 0.002). Deficits in verbal fluency performance were also associated with GG genotype of ABCC3 rs12604031 among HR patients (p = 0.001) as well as with rs1800206 in PPARA gene in low risk groups (p = 0.008). The risk of moderate-severe depression was highest among carriers of the minor G allele of SLCO1B1 rs4149056 who received chemotherapy only (p = 0.002). All variants found significantly associated with tested outcome (except those initially confined to risk subgroup such as those in ABCC3 and PCDHB10 genes) were further analyzed for an association with respective outcomes in an independent cohort of ALL survivors (SJLIFE cohort) (Tables 4 and 5). Two associations were noticeable. The association of borderline significance between deficit in verbal fluency score and the minor allele of MTR rs1805087 was seen in all survivors (OR = 1.7; 95% CI, 1.0-2.8; p = 0.05). The association between deficit in digit span score and GG genotype of CACNB2 rs58225473 showed similar trend as in PETALE cohort for all participants (OR = 3.7, 95% CI, 1.0-13.9), as well as for patients who received chemotherapy only (OR = 3.8, 95% CI, 0.9-16.5), however they did not reach significance (p = 0.08 and 0.09, respectively). The associations for rs58225473 and rs1805087 variants in CACNB2 and MTR genes were significant for combined PETALE and SJLIFE cohorts (S4 Table). The association between deficit in digit span score and GG genotype of CACNB2 rs58225473 was significant for all participants (OR = 3.7; 95% CI, 1.25-11; p = 0.01) and for patients who received chemotherapy only (OR = 7.2, 95% CI, 2.1-25; p = 0.0004). The association between deficit in verbal fluency score and the minor allele of MTR rs1805087 was seen in all survivors (OR = 1.5; 95% CI, 1-2.3; p = 0.04) and in male participants (OR = 1.8; 95% CI, 1-3.1; p = 0.04). Given that the MTR belongs to the MTX pathway, we further explored whether the effect of rs1805087 was modulated by cumulative MTX doses, for which such data were available in the discovery group ( Table 1). The relationship with the deficit in verbal fluency score was particularly obvious in patients who received higher overall cumulative doses (Fig 1, p = 0.01 for patients with cumulative doses above median vs. p = 0.3, for patients with cumulative doses below median). Genetics of neurocognitive complications in survivors of childhood ALL

Rare variants
The analysis of functionally predicted rare variants in PETALE cohort led to the detection of an association between the deficit in trail making test score and rare variants enrichment in SLCO2B1, HSPA4 and GSTT1 genes (p = 0.0002, p = 0.004 and p = 0.003, respectively, Table 6). Using the collapsing approach, we explored variant combinations that contributed to the observed association signal, identifying two variants in GSTT1, three in HSPA4 and four in SLCO2B1 gene. Replication analyses were not performed for these findings because information regarding these variants was not available in the SJLIFE cohort.

Discussion
Functionally predicted germline common variants in MTR, PPARA, SLCO1B1, ABCC3, CALML5, CACNB2 and PCDHB10 genes were found to be significantly associated with deficits in neurocognitive tests performance, whereas a variant in EPHA5 gene was significantly associated with both anxiety and depression.

Neurocognitive performance
Among common variants associated with an impairment in neurocognitive function, rs10904516 in the MTR gene, which was associated with a deficit in verbal fluency, seems particularly interesting given similar observation in the SJLIFE cohort. The MTR gene encodes a B12 dependent methionine synthase involved in remethylation of homocysteine (Hcy), which is the crucial step in methionine production in all types of cells [67]. Mutations in the MTR gene as well as severe deficiency of vitamin B12 could result in elevated concentration of Hcy in plasma and cerebrospinal fluid. Studies have shown that Hcy exerts a neurotoxic action and may participate in the mechanisms of neurodegeneration, such as excitotoxicity, oxidative stress, calcium accumulation, and apoptosis [68][69][70]. MTR gene is involved in the metabolic pathway of MTX. Administration of MTX was associated with acute and subacute neurotoxic effects; these detrimental effects may accumulate over time [69]. The detected common functional polymorphism (rs1805087) leading to Asp919Gly amino acid replacement in the MTR gene could affect enzymatic activity, thus increasing the level of Hcy [68,69,71]. Indeed, we have shown interaction between MTR rs1805087 and cumulative MTX dose in survivors with the deficit in verbal fluency score. This confirms previous finding: this variant together with polymorphisms of other genes that are implicated in the Hcy pathway were already studied in the context of MTX long-term neurotoxicity and has been found to affect neurocognitive function in childhood ALL survivors [72,73]. CACNB2 rs58225473 variant was associated with the neurocognitive deficit as defined by digit span test, which measures working memory. Similar risk values, although not significant, were noted in the SJLIFE cohort in all participants and in the group of survivors who received chemotherapy only. The CACNB2 gene encodes an auxiliary voltage-dependent subunit of Ltype calcium-channel that is mainly expressed in brain and heart tissue. Voltage-dependent calcium channels are crucial for neuronal differentiation and maturation; they induce large number of intracellular events such as neurotransmitter release, neuronal excitability, synaptic plasticity and gene regulation [74]. Calcium influx mediated by those channels has both spatial and temporal components and encodes important signaling information [75]. Moreover, in a recent GWAS study CACNB2 was identified among four significant risk loci underlying genetic effects shared between five major psychiatric disorders that included schizophrenia, autism spectrum disorder, attention deficit-hyperactivity disorder, bipolar and major depressive disorders [76]. Additionally, the rare variants in this gene were found in affected members of families with autism spectrum disease [77]. Given important role that CACNB2 can play, it is not surprising that it was studied as a possible pharmacological target in treatment of mental disorders [78]. The rs58225473 is c.1803T>G substitution (NM_201590.2) leading to Asp601Glu replacement, which is predicted to affect protein function and possibly calcium channel function. Common variants in several other genes influenced neurocognitive decline in PETALE cohort. PPARA gene belongs to PPARs receptor family of ligand-activated transcription factors involved in the regulation of inflammation [79]. Effects of glucocorticoids can be reinforced by PPAR ligands [80]. The enhanced heterodimer formations of PPARA could be associated with increased expression of brain and glial cell-derived neurotrophic factors [81]. The rs1800206 variant in PPARA gene was associated with lower verbal fluency score in females and survivors assigned to SR group or chemotherapy only. Similar association was noted for ABCC3 rs12604031. ABCC3 is a member of the superfamily of ATP-binding cassette (ABC) transporters, and the bioavailability of MTX may be affected by this transporter [82].
CALML5 gene is related to the calmodulin family of calcium binding proteins highly implicated in CNS function [83]. Its protective role and implication in inhibition of neuronal death was described in Alzheimer disease [84,85]. We observed the significant association of variant rs10904516 and deficit in verbal fluency score in all survivors with the stronger effect seen in the HR group. Similar association mostly confined to HR group was noted between lower score on trial making test and rs2907323 in the PCDHB10 gene. PCDH (protocadherin) genes, are expressed in the central and peripheral nervous systems and are required for their normal development. They mediate a variety of processes, including neuronal survival, morphogenesis and connectivity, synaptic maintenance, and spatial patterning of axons and dendrites [86]. The variants in PCDH genes have been reported to be associated with dyslexia and bipolar disorder [87,88].
The neurotoxic effects of treatments in childhood ALL have been the subject of multiple investigations. These effects consist of central neurotoxicity clearly noticeable by encephalopathy and/or neurodevelopmental cognitive deficits [89][90][91], particularly in survivors exposed to a highly intensified treatment protocols with CNS-directed chemotherapy, even in the absence of CRT [4,19,92,93]. Cognitive impairment and information processing have been associated with intensity and duration of CS treatment [94,95]. Female survivors were reported to have more severe short-term memory impairment [59] and lower scores on attentional indices, cognitive flexibility [96,97] and visuomotor control [4]. Female childhood ALL survivors are more likely to present cerebral white matter damage [98] that may affect cognitive functioning. Congruent with these previous observations, several associations detected in the present study were modulated by sex and treatment intensity (reflected by presence or not of CRT or risk groups).

Anxiety and depression
We studied dimensions of internalized symptoms which are frequent in normative populations, namely anxiety and depression [99], highlighting impairment in mental quality of life of childhood ALL survivors [21,22,100]. Although anxiety and depression measures are not equivalent to clinical diagnosis derived from the gold standard systematic interview [101], moderate-severe levels are generally interpreted as a risk for having clinically relevant anxiety or depression.
The rs11556218 in the EPHA5 gene was associated with higher risk of both anxiety and depression that was further potentiated in female patients. The EPHA5 gene codes for a brainspecific kinase that is selectively expressed in a subset of serotonin neurons during embryonic and postnatal development [102]. Receptors in the EPH subfamily modify the strength of existing synapses in the adult brain [103]. Divergent vulnerabilities between females and males could be explained by gender differences in brain maturation [104,105], which might make females more vulnerable to the neurotoxic effects of chemotherapy. Other assumptions, such as endocrine factors, have also been hypothesized to explain sex differences in the susceptibility [59].
Moreover, we identified the association between depression and the presence of variant rs4149056 in SLCO1B1 gene. This association was detected in the group of survivors that received chemotherapy without CRT. SLCO1B1 gene encodes a liver-specific member of the organic anion transporter family involved in hepatic uptake of MTX. This association deserved further attention given that the same variant was detected through genome wide association studies to contribute to inter-individual variability in the clearance of high-dose MTX [106]. It was subsequently replicated in independent cohorts and shown also as a predictor of shortterm toxicity following MTX treatment [106][107][108][109][110][111]. MTX treatment has been associated with adverse emotional or behavioral outcomes [20], thus these results could justify further studies of SLCO1B1 gene in related contexts.

Rare variants' analysis
The association between deficits in the trail making test score was identified in relation to rare variants enrichment in HSPA4, SLCO2B1 and GSTT1 genes, with a very strong individual contribution of rs61745470 in HSPA4 gene. This variant was recently associated with familial genetic risk for suicide (as well as with risk for psychiatric or substance abuse conditions) [112]. The SLCO2B1 and GSTT1 genes are highly implicated in physiological and pharmacological distribution of drugs and endogenous molecules. The SLCO2B1, a member of the organic anion transporting polypeptide (OATP) family, is involved in steroid hormone uptake and transport of steroid conjugates [113,114]. GSTT1 was recently associated with higher risk for early onset of severe mental and bipolar disorders [115]. We also evaluated the association of deletion polymorphisms of GSTT1 [116] (found with the frequency of 23.7% in discovery cohort) with the deficits in the trail making test score. There was no association of GSTT1 null genotype with the deficits in the trail making test score.
The impact of here identified rare functional variants requires further investigation.

Concluding remarks
Our study has certain limitations. Its limited sample size may affect the accuracy of the results, particularly in the context of the stratified analysis. Other unmeasured factors in this study, for example, inflammation and oxidative stress, could modulate or potentiate association with genetic factors. The candidate gene approach may have missed genetic markers potentially involved in neurocognitive decline and mood disturbances that could have been detected through unbiased approaches. Among associations detected in the PETALE cohort only two showed a similar trend in SJLIFE cohort. Despite matching both outcomes and patients' characteristics between the two cohorts, it is possible that small sample size, differences in treatment protocols or time of ALL diagnosis [117][118][119][120] contributed to the observed discrepancies. Likewise, stratification by risk group designation was not available for the SJLIFE cohort, precluding replication of the risk-based stratified analyses. Although the analyses in PETALE cohort were corrected for multiple testing, and confounding was reduced due to homogeneous population and uniform treatment, we cannot exclude that some of the associations have been obtained by chance.
In conclusion, using a comprehensive candidate gene approach and whole exome sequencing data we identified a panel of functionally predicted genetic variants significantly associated with neurocognitive deficits, anxiety and depression in childhood ALL survivors. Additional exome wide analysis might lead to the discovery of novel genes and genetic variants associated with neurocognitive LAEs as well as with the mood disorders.
While we acknowledge that the identified germline variants still need to be evaluated and validated through replication and functional studies, the current findings can help further understanding of the influence of genetic component on long-term complications related to cancer therapy.
Supporting information S1 Fig. Principal component analysis (PCA). PCA analysis comparing sequencing data of 400 leukemia patients (including PETALE cohort) from Sainte-Justine University Health Center (SJUHC) to the HapMap genotype reference data (release 23) for Europeans (EUR), East Asians (EAS) and Africans (AFR). PC1, Principal Component 1; PC2, Principal Component 2. (TIF) S1 Table. Genotyping: Identity of polymorphisms, details of PCR and ASO hybridization. R, reverse, F, forward. The base substitution that distinguishes the two variants of each polymorphism is given in bold for ASO probes. dbSNP number is provided. Ancestral allele is given in bold and minor allele is underlined. The polymorphisms are presented as a change from ancestral to derived allele, unless ancestral allele is not known, when the change is given from major to minor allele. SNPs in coding region leading or not to amino-acid substitutions are indicated. (DOCX) S2 Table. Significant results of association study of common variants from the candidate genes of relevance for nervous system function, PETALE cohort, WES data (n = 191). a Association test based on comparing allele frequencies between cases and controls. All associations have FDR-BH (Benjamini-Hochberg false discovery rate) lower than 5%. All also have p value lower than 0.001, which is Bonferroni cut-off value for the number of variants tested in nervous system function pathway.  Table. Results of association study of common variants from methotrexate and corticosteroids pathways, PETALE cohort, WES data (n = 191). a Association test based on comparing allele frequencies between cases and controls. All associations have FDR-BH (Benjamini-Hochberg false discovery rate) lower than 5%. b Stratified analyses according to sex and treatment intensity (standard vs high risk); chemotherapy only vs chemotherapy and cranial radiation therapy (CRT). c SNPs or associations that did not qualify for genotyping with p value higher than 0.0019