Inter-individual variation in ovarian reserve after gonadotoxic treatment in female childhood cancer survivors - a genome-wide association study: results from PanCareLIFE.

OBJECTIVE: We aimed to discover new variants associated with low ovarian reserve after gonadotoxic treatment among adult female childhood cancer survivors using a genome-wide association study approach. DESIGN: Genome-wide association study. SUBJECTS: A discovery cohort of adult female childhood cancer survivors, from the pan-European PanCareLIFE cohort (n=743; median age: 25.8 years), excluding those who received bilateral ovarian irradiation, bilateral oophorectomy, central nerve system or total body irradiation, or stem cell transplantation. Replication was attempted in the USA-based St. Jude Lifetime Cohort (n=391; median age: 31.3 years). EXPOSURE: Female childhood cancer survivors are at risk of therapy-related gonadal impairment. Alkylating agents are well-established risk factors, and the inter-individual variability in gonadotoxicity may be explained by genetic polymorphisms. Data were collected in real-life conditions and cyclophosphamide equivalent dose was used to quantify alkylation agent exposure. INTERVENTION: No intervention was performed. MAIN OUTCOME MEASURE: Anti-Müllerian hormone (AMH) levels served as a proxy for ovarian function and findings were combined in a meta-analysis. RESULTS: Three genome-wide significant (<5.0x10-8) and 16 genome-wide suggestive (<5.0x10-6) loci were associated with log-transformed AMH levels, adjusted for cyclophosphamide equivalent dose of alkylating agents, age at diagnosis, and age at study in the PanCareLIFE cohort. Based on effect allele frequency (EAF) (>0.01 if not genome-wide significant), p-value (<5.0×10-6), and biological relevance, 15 SNPs were selected for replication. None of the SNPs were statistically significantly associated with AMH levels. A meta-analysis indicated that rs78861946 was associated at borderline genome-wide statistical significance (Reference/effect allele: C/T; EAF: 0.04, Beta (SE): -0.484 (0.091), p-value= 9.39×10-8). CONCLUSION: This study found no genetic variants associated with a lower ovarian reserve after gonadotoxic treatment, as the findings of this GWAS were not statistically significant replicated in the replication cohort. Suggestive evidence for potential importance of one variant is briefly discussed, but the lack of statistical significance calls for larger cohort sizes. As the population of childhood cancer survivors is increasing, large-scale and systematic research is needed to identify genetic variants that could aid predictive risk models of gonadotoxicity and as well as fertility preservation options for childhood cancer survivors.

Trans-Ancestral Genetic Risk Factors for Treatment-Related Type 2 Diabetes Mellitus in Survivors of Childhood Cancer.

PURPOSE: Type 2 diabetes mellitus (T2D) is a prevalent long-term complication of treatment in survivors of childhood cancer, with marked racial/ethnic differences in burden. In this study, we investigated trans-ancestral genetic risks for treatment-related T2D. PATIENTS AND METHODS: Leveraging whole-genome sequencing data from the St Jude Lifetime Cohort (N = 3,676, 304 clinically ascertained cases), we conducted ancestry-specific genome-wide association studies among survivors of African and European genetic ancestry (AFR and EUR, respectively) followed by trans-ancestry meta-analysis. Trans-/within-ancestry replication including data from the Childhood Cancer Survivor Study (N = 5,965) was required for prioritization. Three external general population T2D polygenic risk scores (PRSs) were assessed, including multiancestry PRSs. Treatment risk effect modification was evaluated for prioritized loci. RESULTS: Four novel T2D risk loci showing trans-/within-ancestry replication evidence were identified, with three loci achieving genome-wide significance (P < 5 × 10-8). Among these, common variants at 5p15.2 (LINC02112), 2p25.3 (MYT1L), and 19p12 (ZNF492) showed evidence of modifying alkylating agent-related T2D risk in both ancestral groups, but showed disproportionately greater risk in AFR survivors (AFR odds ratios [ORs], 3.95-17.81; EUR ORs, 2.37-3.32). In survivor-specific RNA-sequencing data (N = 207), the 19p12 locus variant was associated with greater ZNF492 expression dysregulation after exposures to alkylators. Elevated T2D risks across ancestry groups were only observed with increasing values for multiancestry T2D PRSs and were especially increased among survivors treated with alkylators (top v bottom quintiles: ORAFR, 20.18; P = .023; OREUR, 13.44; P = 1.3 × 10-9). CONCLUSION: Our findings suggest therapy-related genetic risks contribute to the increased T2D burden among non-Hispanic Black childhood cancer survivors. Additional study of how therapy-related genetic susceptibility contributes to this disparity is needed.

Genetic Factors and Long-term Treatment-Related Neurocognitive Deficits, Anxiety, and Depression in Childhood Leukemia Survivors: An Exome-Wide Association Study.

BACKGROUND: An increased risk of neurocognitive deficits, anxiety, and depression has been reported in childhood cancer survivors. METHODS: We analyzed associations of neurocognitive deficits, as well as anxiety and depression, with common and rare genetic variants derived from whole-exome sequencing data of acute lymphoblastic leukemia (ALL) survivors from the PETALE cohort. In addition, significant associations were assessed using stratified and multivariable analyses. Next, top-ranking common associations were analyzed in an independent SJLIFE replication cohort of ALL survivors. RESULTS: Significant associations were identified in the entire discovery cohort (N = 229) between the AK8 gene and changes in neurocognitive function, whereas PTPRZ1, MUC16, TNRC6C-AS1 were associated with anxiety. Following stratification according to sex, the ZNF382 gene was linked to a neurocognitive deficit in males, whereas APOL2 and C6orf165 were associated with anxiety and EXO5 with depression. Following stratification according to prognostic risk groups, the modulatory effect of rare variants on depression was additionally found in the CYP2W1 and PCMTD1 genes. In the replication SJLIFE cohort (N = 688), the male-specific association in the ZNF382 gene was not significant; however, a P value<0.05 was observed when the entire SJLIFE cohort was analyzed. ZNF382 was significant in males in the combined cohorts as shown by meta-analyses as well as the depression-associated gene EXO5. CONCLUSIONS: Further research is needed to confirm whether the current findings, along with other known risk factors, may be valuable in identifying patients at increased risk of these long-term complications. IMPACT: Our results suggest that specific genes may be related to increased neuropsychological consequences.

Development and validation of age-specific risk prediction models for primary ovarian insufficiency in long-term survivors of childhood cancer: a report from the Childhood Cancer Survivor Study and St Jude Lifetime Cohort.

BACKGROUND: Female survivors of childhood cancer are at risk for primary ovarian insufficiency (POI), defined as the cessation of gonadal function before the age of 40 years. We aimed to develop and validate models to predict age-specific POI risk among long-term survivors of childhood cancer. METHODS: To develop models to predict age-specific POI risk for the ages of 21-40 years, we used data from the Childhood Cancer Survivor Study (CCSS). Female survivors aged 18 years or older at their latest follow-up, with self-reported menstrual history information and free of subsequent malignant neoplasms within 5 years of diagnosis, were included. We evaluated models that used algorithms based on statistical or machine learning to consider all predictors, including cancer treatments. Cross-validated prediction performance metrics (eg, area under the receiver operating characteristic curve [AUROC]) were compared to select the best-performing models. For external validation of the models, we used data from 5-year survivors in the St Jude Lifetime Cohort (SJLIFE) with ovarian status clinically ascertained using hormone measurements (menopause defined by follicle stimulating hormone >30 mIU/mL and oestradiol <17 pg/mL) and medical chart or questionnaire review. We also evaluated an SJLIFE-based polygenic risk score for POI among 1985 CCSS survivors with genotype data available. FINDINGS: 7891 female CCSS survivors (922 with POI) were included in the development of the POI risk prediction model, and 1349 female SJLIFE survivors (101 with POI) were included in the validation study. Median follow-up from cancer diagnosis was 23·7 years (IQR 18·3-30·0) in CCSS and 15·1 years (10·4-22·9) in SJLIFE. Between the ages of 21 and 40 years, POI prevalence increased from 7·9% (95% CI 7·3-8·5) to 18·6% (17·3-20·0) in CCSS and 7·3% (5·8-8·9) to 14·9% (11·6-19·1) in SJLIFE. Age-specific logistic regression models considering ovarian radiation dosimetry or prescribed pelvic and abdominal radiation dose, along with individual chemotherapy predictors, performed well in CCSS. In the SJLIFE validation, the prescribed radiation dose model performed well (AUROC 0·88-0·95), as did a simpler model that considered any exposures to pelvic or abdominal radiotherapy or alkylators (0·82-0·90). Addition of the polygenic risk predictor significantly improved the average positive predictive value (from 0·76 [95% CI 0·63-0·89] to 0·87 [0·80-0·94]; p=0·029) among CCSS survivors treated with ovarian radiation and chemotherapy. INTERPRETATION: POI risk prediction models using treatment information showed robust prediction performance in adult survivors of childhood cancer. FUNDING: Canadian Institutes of Health Research, US National Cancer Institute.

Trans-ancestral genetic study of diabetes mellitus risk in survivors of childhood cancer: a report from the St. Jude Lifetime Cohort and the Childhood Cancer Survivor Study.

Type 2 diabetes mellitus (T2D) is an established late effect of treatment for childhood cancer. Leveraging detailed cancer treatment and whole-genome sequencing data among survivors of childhood cancer of European (EUR) and African (AFR) genetic ancestry in the St. Jude Lifetime Cohort (N=3,676; 304 cases), five novel diabetes mellitus (DM) risk loci were identified with independent trans-/within-ancestry replication, including in 5,965 survivors of the Childhood Cancer Survivor Study. Among these, common risk variants at 5p15.2 ( LINC02112 ), 2p25.3 ( MYT1L ), and 19p12 ( ZNF492 ) modified alkylating agent-related risks across ancestry groups, but AFR survivors with risk alleles experienced disproportionately greater risk of DM (AFR, variant ORs: 3.95-17.81; EUR, variant ORs: 2.37-3.32). Novel risk locus XNDC1N was identified in the first genome-wide DM rare variant burden association analysis in survivors (OR=8.65, 95% CI: 3.02-24.74, P=8.1×10 -6 ). Lastly, a general-population 338-variant multi-ancestry T2D polygenic risk score was informative for DM risk in AFR survivors, and showed elevated DM odds after alkylating agent exposures (quintiles: combined OR EUR =8.43, P=1.1×10 -8 ; OR AFR =13.85, P=0.033). This study supports future precision diabetes surveillance/survivorship care for all childhood cancer survivors, including those with AFR ancestry.

Inference on the Genetic Architecture of Breast Cancer Risk.

BACKGROUND: What are the major determinants of women's breast cancer risk? Rare mutations such as those in the BRCA1/2 genes, polygenic scores of common alleles identified by genome-wide association studies, or nongenetic factors? METHODS: The population-based Nordic Twin Study of Cancer, with 3,933 breast cancer cases among 21,054 monozygotic (MZ) and 30,939 dizygotic (DZ) female twin pairs, provides three key clues to this question: (i) the average lifetime risk, approximately 8%, does not differ by twin zygosity; (ii) the mean time interval between diagnoses when both twins develop disease (i.e., disease concordance) also does not differ by zygosity; but, (iii) conditioning on one twin having developed disease, the incidence rate in the co-twin is approximately 1% per year if the pair is MZ and 0.5% per year if DZ. RESULTS: Assuming that nongenetic risk factors are shared similarly between twins regardless of zygosity, we can draw two conclusions from (i) to (iii). CONCLUSIONS: First, (i) and (iii) imply that the chief determinant of risk is in the germline DNA, because the conditional incidence rate is several-fold higher than the average risk (8% lifetime) in MZ twins but only half as much in DZ twins. Second, the seeming inconsistency between the two-fold conditional incidence rate (iii) and the equality of the mean inter-twin disease intervals in disease concordance (ii) can be resolved if the risk factors in the germline DNA are rare variants, not common variants. IMPACT: This paper details simple deductive reasoning for these conclusions and draws a critical inference regarding breast cancer etiology. See related In the Spotlight, p. 1477.

Genome-Wide Analysis of Rare Haplotypes Associated with Breast Cancer Risk.

Numerous common genetic variants have been linked to breast cancer risk, but they only partially explain the total breast cancer heritability. Inference from Nordic population-based twin data indicates rare high-risk loci as the chief determinant of breast cancer risk. Here, we use haplotypes, rather than single variants, to identify rare high-risk loci for breast cancer. With computationally phased genotypes from 181,034 white British women in the UK Biobank, a genome-wide haplotype-breast cancer association analysis was conducted using sliding windows of 5 to 500 consecutive array-genotyped variants. In the discovery stage, haplotype-breast cancer associations were evaluated retrospectively in the prestudy-enrollment data including 5,487 breast cancer cases. Breast cancer hazard ratios (HR) for additive haplotypic effects were estimated using Cox regression. The replication analysis included a prospective cohort of women free of breast cancer at enrollment, of whom 3,524 later developed breast cancer. This two-stage analysis detected 13 rare loci (frequency <1%), each associated with an appreciable breast cancer-risk increase (discovery: HRs = 2.84-6.10, P < 5 × 10-8; replication: HRs = 2.08-5.61, P < 0.01). In contrast, the variants that formed these rare haplotypes individually exhibited much smaller effects. Functional annotation revealed extensive cis-regulatory DNA elements in breast cancer-related cells underlying the replicated rare haplotypes. Using phased, imputed genotypes from 30,064 cases and 25,282 controls in the DRIVE OncoArray case-control study, 6 of the 13 rare-loci associations were found generalizable (odds ratio estimates: 1.48-7.67, P < 0.05). This study demonstrates the complementary advantage of utilizing rare haplotypes to capture novel risk loci and suggests the potential for the discovery of more genetic elements contributing to cancer heritability as large data sets of germline whole-genome sequencing become available. SIGNIFICANCE: A genome-wide two-stage haplotype analysis identifies rare haplotypes associated with breast cancer risk and suggests that the rare risk haplotypes represent long-range interactions with regulatory consequences influencing cancer risk.

Germline Genetic and Treatment-Related Risk Factors for Diabetes Mellitus in Survivors of Childhood Cancer: A Report From the Childhood Cancer Survivor Study and St Jude Lifetime Cohorts.

PURPOSE: To characterize germline genetic risk factors of diabetes mellitus among long-term survivors of childhood cancer. METHODS: Adult survivors of childhood cancer from the Childhood Cancer Survivor Study (CCSS) Original Cohort (n = 5,083; 383 with diabetes) were used to conduct a discovery genome-wide association study. Replication was performed using the CCSS Expansion (n = 2,588; 40 with diabetes) and the St Jude Lifetime (SJLIFE; n = 3,351; 208 with diabetes) cohorts. Risk prediction models, stratified on exposure to abdominal radiation, were calculated using logistic regression including attained age, sex and body mass index, diagnosis, alkylating chemotherapy, age at cancer diagnosis, and a polygenic risk score (PRS) on the basis of 395 diabetes variants from the general population. Area under the receiver operating characteristic curve (AUC) was calculated for models on the basis of traditional risk factors, clinical risk factors, and PRS. RESULTS: There was a genome-wide significant association of rs55849673-A with diabetes among survivors (odds ratio, 2.9; 95% CI, 2.0 to 4.2; P = 3.7 × 10-8), which is related to expression of ERCC6L2 in the Genotype-Tissue Expression project. The association of rs55849673-A was observed largely among survivors not exposed to abdominal radiation (odds ratio = 3.5, P = 1.1 × 10-7) and the frequency of rs55849673-A was consistently higher among diabetic survivors in the CCSS Expansion and SJLIFE cohorts. Risk prediction models including traditional diabetes risk factors, clinical risk factors and PRS had an optimism-corrected AUC of 0.801, with an AUC of 0.751 in survivors treated with abdominal radiation versus 0.813 in survivors who did not receive abdominal radiation. CONCLUSION: There is evidence for a novel locus of diabetes among survivors not exposed to abdominal radiation. Further refinement and validation of clinic-based risk prediction models for diabetes among long-term survivors of childhood cancer is warranted.

Genome-Wide Association Study of Pregnancy in Childhood Cancer Survivors: A Report from the Childhood Cancer Survivor Study.

BACKGROUND: Gonadotoxic treatment-related infertility has a significant impact on quality of life in childhood cancer survivors. Genome-wide association analyses to delineate the risk of infertility in childhood cancer survivors have not been previously reported. METHODS: Leveraging genotype data from a large survivor cohort, the Childhood Cancer Survivor Study (CCSS), we investigated the role of SNPs on future pregnancy or siring a pregnancy in survivors without pelvic, testicular, or brain radiation who had ever been married. We calculated sex-stratified hazard ratios, using Cox proportional hazards modeling, adjusting for birth cohort (before 1965 vs. 1965 or later) and doses of relevant chemotherapies; replication was attempted in the independent St. Jude Lifetime Cohort study (SJLIFE). RESULTS: In the CCSS cohort, nine SNPs were found to be suggestive (P < 10-7) or statistically significantly (P < 5 × 10-8) associated with pregnancy, however, none of the SNPs were replicated in SJLIFE. Cohorts differed based on the overall pregnancy rate, frequency of sterilizing procedures, and birth cohort. CONCLUSIONS: We were not able to replicate our findings of SNPs associated with pregnancy in childhood cancer survivors. IMPACT: Future attempts at replication should be considered in cohorts treated in a comparable era. In addition, understanding the role of genetics in fertility in childhood cancer survivors may be better approached using more advanced sequencing techniques.

Forgoing needed medical care among long-term survivors of childhood cancer: racial/ethnic-insurance disparities.

PURPOSE: To investigate racial/ethnic-related disparities by insurance status in "forgoing needed medical care in the last year due to finances" in childhood cancer survivors. METHODS: Our study included 3310 non-Hispanic/Latinx White, 562 non-Hispanic/Latinx Black, and 92 Hispanic/Latinx survivors from the St. Jude Lifetime Cohort Study. Logistic regression analyses, guided by Andersen's Healthcare Utilization Model, were adjusted for "predisposing" (survey age, sex, childhood cancer diagnosis and treatment, and treatment era) and "need" (perceived health status) factors. Additional adjustment for household income/education and chronic health conditions was considered. RESULTS: Risk of forgoing care was highest among non-Hispanic/Latinx Blacks and lowest among Hispanics/Latinxs for each insurance status. Among privately insured survivors, relative to non-Hispanic/Latinx Whites, non-Hispanic/Latinx Blacks were more likely to forgo care (adjusted OR: 1.82, 95% CI: 1.30-2.54): this disparity remained despite additional adjustment for household income/education (adjusted OR: 1.43, 95% CI: 1.01-2.01). In contrast, publicly insured survivors, regardless of race/ethnicity, had similar risk of forgoing care as privately insured non-Hispanic/Latinx Whites. All uninsured survivors had high risk of forgoing care. Additional chronic health condition adjustment did not alter these results. CONCLUSIONS: Provision of public insurance to all childhood cancer survivors may diminish racial/ethnic disparities in forgoing care that exist among the privately insured and reduce the risk of forgoing care among uninsured survivors to that of privately insured non-Hispanic/Latinx Whites. IMPLICATIONS FOR CANCER SURVIVORS: Under public insurance, childhood cancer survivors had low risk of forgoing care, at the similar level to privately insured non-Hispanic/Latinx Whites, regardless of race/ethnicity.

Effect of Genetic Variation in CYP450 on Gonadal Impairment in a European Cohort of Female Childhood Cancer Survivors, Based on a Candidate Gene Approach: Results from the PanCareLIFE Study.

BACKGROUND: Female childhood cancer survivors (CCSs) carry a risk of therapy-related gonadal dysfunction. Alkylating agents (AA) are well-established risk factors, yet inter-individual variability in ovarian function is observed. Polymorphisms in CYP450 enzymes may explain this variability in AA-induced ovarian damage. We aimed to evaluate associations between previously identified genetic polymorphisms in CYP450 enzymes and AA-related ovarian function among adult CCSs. METHODS: Anti-Müllerian hormone (AMH) levels served as a proxy for ovarian function in a discovery cohort of adult female CCSs, from the pan-European PanCareLIFE cohort (n = 743; age (years): median 25.8, interquartile range (IQR) 22.1-30.6). Using two additive genetic models in linear and logistic regression, nine genetic variants in three CYP450 enzymes were analyzed in relation to cyclophosphamide equivalent dose (CED) score and their impact on AMH levels. The main model evaluated the effect of the variant on AMH and the interaction model evaluated the modifying effect of the variant on the impact of CED score on log-transformed AMH levels. Results were validated, and meta-analysis performed, using the USA-based St. Jude Lifetime Cohort (n = 391; age (years): median 31.3, IQR 26.6-37.4). RESULTS: CYP3A4*3 was significantly associated with AMH levels in the discovery and replication cohort. Meta-analysis revealed a significant main deleterious effect (Beta (95% CI): -0.706 (-1.11--0.298), p-value = 7 × 10-4) of CYP3A4*3 (rs4986910) on log-transformed AMH levels. CYP2B6*2 (rs8192709) showed a significant protective interaction effect (Beta (95% CI): 0.527 (0.126-0.928), p-value = 0.01) on log-transformed AMH levels in CCSs receiving more than 8000 mg/m2 CED. CONCLUSIONS: Female CCSs CYP3A4*3 carriers had significantly lower AMH levels, and CYP2B6*2 may have a protective effect on AMH levels. Identification of risk-contributing variants may improve individualized counselling regarding the treatment-related risk of infertility and fertility preservation options.

Genetic factors in treatment-related cardiovascular complications in survivors of childhood acute lymphoblastic leukemia.

Aim: Cardiovascular disease represents one of the main causes of secondary morbidity and mortality in patients with childhood cancer. Patients & methods: To further address this issue, we analyzed cardiovascular complications in relation to common and rare genetic variants derived through whole-exome sequencing from childhood acute lymphoblastic leukemia survivors (PETALE cohort). Results: Significant associations were detected among common variants in the TTN gene, left ventricular ejection fraction (p ≤ 0.0005), and fractional shortening (p ≤ 0.001). Rare variants enrichment in the NOS1, ABCG2 and NOD2 was observed in relation to left ventricular ejection fraction, and in NOD2 and ZNF267 genes in relation to fractional shortening. Following stratification according to risk groups, the modulatory effect of rare variants was additionally found in the CBR1, ABCC5 and AKR1C3 genes. None of the associations was replicated in St-Jude Lifetime Cohort Study. Conclusion: Further studies are needed to confirm whether the described genetic markers may be useful in identifying patients at increased risk of these complications.

Clinical and genetic risk factors for radiation-associated ototoxicity: A report from the Childhood Cancer Survivor Study and the St. Jude Lifetime Cohort.

BACKGROUND: Cranial radiation therapy (CRT) is associated with ototoxicity, which manifests as hearing loss and tinnitus. The authors sought to identify clinical determinants and genetic risk factors for ototoxicity among adult survivors of pediatric cancer treated with CRT. METHODS: Logistic regression evaluated associations of tinnitus (n = 1991) and hearing loss (n = 2198) with nongenetic risk factors and comorbidities among CRT-treated survivors in the Childhood Cancer Survivor Study. Genome-wide association studies (GWASs) of CRT-related tinnitus and hearing loss were also performed. RESULTS: Males were more likely to report CRT-related tinnitus (9.4% vs 5.4%; P = 5.1 × 10-4 ) and hearing loss (14.0% vs 10.7%; P = .02) than females. Survivors with tinnitus or hearing loss were more likely to experience persistent dizziness or vertigo (tinnitus: P < 2 × 10-16 ; hearing loss: P = 6.4 × 10-9 ), take antidepressants (tinnitus: P = .02; hearing loss: P = .01), and report poorer overall health (tinnitus: P = 1.5 × 10-6 ; hearing loss: P = 1.7 × 10-6 ) in comparison with controls. GWAS of CRT-related tinnitus revealed a genome-wide significant signal in chromosome 1 led by rs203248 (P = 1.5 × 10-9 ), whereas GWAS of CRT-related hearing loss identified rs332013 (P = 5.8 × 10-7 ) in chromosome 8 and rs67522722 (P = 7.8 × 10-7 ) in chromosome 6 as nearly genome-wide significant. A replication analysis identified rs67522722, intronic to ATXN1, as being significantly associated with CRT-related hearing loss (P = .03) and de novo hearing loss (P = 3.6 × 10-4 ). CONCLUSIONS: CRT-associated ototoxicity was associated with sex, several neuro-otological symptoms, increased antidepressant use, and poorer self-reported health. GWAS of CRT-related hearing loss identified rs67522722, which was supported in an independent cohort of survivors. LAY SUMMARY: Hearing loss and subjective tinnitus (the perception of noise or ringing in the ear) are long-term side effects of cancer treatment and are common in children treated with radiation to the brain. These toxicities can affect childhood development and potentially contribute to serious learning and behavioral difficulties. This study's data indicate that males are at greater risk for hearing loss and tinnitus than females after radiation therapy to the brain. Those who develop these toxicities are more likely to use antidepressants and report poorer overall health. Health care providers can improve the management of survivors by informing patients and/or their parents of these risks.

A Novel Locus Predicts Spermatogenic Recovery among Childhood Cancer Survivors Exposed to Alkylating Agents.

Exposure to high doses of alkylating agents is associated with increased risk of impaired spermatogenesis among nonirradiated male survivors of childhood cancer, but there is substantial variation in this risk. Here we conducted a genetic study for impaired spermatogenesis utilizing whole-genome sequencing data from 167 nonirradiated male childhood cancer survivors of European ancestry from the St. Jude Lifetime Cohort treated with cyclophosphamide equivalent dose (CED) ≥4,000 mg/m2. Sperm concentration from semen analysis was assessed as the primary outcome. Common variants (MAF > 0.05) were adjusted for age at cancer diagnosis, CED, and top principal components. Rare/low-frequency variants (MAF ≤ 0.05) were evaluated jointly by various functional annotations and 4-kb sliding windows. A novel locus at 7q21.3 containing TAC1/ASNS was associated with decreased sperm concentration (rs7784118: P = 3.5 × 10-8). This association was replicated in two independent samples of SJLIFE survivors of European ancestry, including 34 nonirradiated male survivors treated with 0 < CED < 4,000 mg/m2 (P = 3.1 × 10-4) and 24 male survivors treated with CED ≥4,000 mg/m2 and radiotherapy <40 Gray (P = 0.012). No association was observed among survivors not exposed to alkylating agents included in the CED (P > 0.29). rs7784118 conferred 3.48- and 9.73-fold increases in risk for clinically defined oligospermia and azoospermia and improved prediction of normospermic, oligospermic, and azoospermic states by 13.7%, 5.3%, and 21.7%. rs7784118 was associated with decreased testosterone level, increased levels of follicle stimulating and luteinizing hormones, and 8.52-fold increased risk of Leydig cell failure. Additional research is warranted to determine how this SNP influences spermatogenesis and to assess its clinical utility in characterizing high-risk survivors and guiding intervention strategies. SIGNIFICANCE: The identified genetic markers harbor potential clinical utility in characterizing high-risk survivors and guiding intervention strategies including pretreatment patient counseling and use of fertility preservation services.