Widely used as anti-cancer and immunosuppressive agents, thiopurines have narrow therapeutic indices due to frequent toxicities, partly explained by TPMT genetic polymorphisms. Recent studies identified germline NUDT15 variation as another critical determinant of thiopurine intolerance, but the underlying molecular mechanisms and its clinical implications remain unknown. In 270 children enrolled in clinical trials for acute lymphoblastic leukemia in Guatemala, Singapore, and Japan, we identified 4 NUDT15 coding variants (p.Arg139Cys, p.Arg139His, p.Val18Ile, p.Val18_Val19insGlyVal) that resulted in 74.4%–100% loss of nucleotide diphosphatase activity. Loss-of-function NUDT15 diplotypes were consistently associated with thiopurine intolerance across three cohorts (P=0.021, 2.1×10−5, and 0.0054, respectively; meta-analysis P=4.45×10−8, allelic effect size=−11.5). Mechanistically, NUDT15 inactivated thiopurine metabolites and decreased its cytotoxicity in vitro, and patients with defective NUDT15 alleles showed excessive thiopurine active metabolites and toxicity. Taken together, our results indicate that a comprehensive pharmacogenetic model integrating NUDT15 variants may inform personalized thiopurine therapy.

Recent genomic profiling of childhood acute lymphoblastic leukemia (ALL) identified a novel high-risk subtype with a gene expression signature resembling Philadelphia chromosome-positive ALL and a poor prognosis (Ph-like ALL). However, the role of inherited genetic variation in Ph-like ALL pathogenesis remains unknown. In a genome-wide association study (GWAS) of 511 ALL cases and 6,661 non-ALL controls, we identified a single susceptibility locus for Ph-like ALL (GATA3, rs3824662, P=2.17×10−14, odds ratio [OR]=3.85, for Ph-like ALL vs. non-ALL; P=1.05×10−8, OR=3.25, for Ph-like ALL vs. non-Ph-like ALL) that was independently validated. The rs3824662 risk allele was associated with somatic lesions underlying Ph-like ALL (i.e., CRLF2 rearrangement, JAK mutation, and IKZF1 deletion) and directly influenced GATA3 transcription. Finally, GATA3 SNP genotype was also associated with early treatment response and the risk of ALL relapse. Our results provide insights into interactions between host and tumor genomes and their importance in ALL pathogenesis and prognosis.

As a prototype of genomics-guided precision medicine, individualized thiopurine dosing based on pharmacogenetics is a highly effective way to mitigate hematopoietic toxicity of this class of drugs. Recently, NUDT15 deficiency was identified as a genetic cause of thiopurine toxicity, and NUDT15-informed preemptive dose reduction was quickly adopted in clinical settings. To exhaustively identify pharmacogenetic variants in this gene, we developed massively parallel NUDT15 function assays to determine the variants’ effect on protein abundance and thiopurine cytotoxicity. Of the 3,097 possible missense variants, we characterized the abundance of 2,922 variants and found 54 hotspot residues at which variants resulted in complete loss of protein stability. Analyzing 2,935 variants in the thiopurine cytotoxicity-based assay, we identified 17 additional residues where variants altered NUDT15 activity without affecting protein stability. We identified structural elements key to NUDT15 stability and/or catalytical activity with single amino acid resolution. Functional effects for NUDT15 variants accurately predicted toxicity risk alleles in patients treated with thiopurines with far superior sensitivity and specificity compared to bioinformatic prediction algorithms. In conclusion, our massively parallel variant function assays identified 1,152 deleterious NUDT15 variants, providing a comprehensive reference of variant function and vastly improving the ability to implement pharmacogenetics-guided thiopurine treatment individualization.

Bridging the survival gap for children with cancer, between those (the great majority) in low and middle income countries (LMIC) and their economically advantaged counterparts, is a challenge that has been addressed by twinning institutions in high income countries with centers in LMIC. The long-established partnership between a Central American consortium--Asociación de Hemato-Oncología Pediátrica de Centro América (AHOPCA)--and institutions in Europe and North America provides a striking example of such a twinning program. The demonstrable success of this endeavor offers a model for improving the health outcomes of children with cancer worldwide. As this remarkable enterprise celebrates its 15th anniversary, it is appropriate to reflect on its origin, subsequent growth and development, and the lessons it provides for others embarking on or already engaged in similar journeys. Many challenges have been encountered and not all yet overcome. Commitment to the endeavor, collaboration in its achievements and determination to overcome obstacles collectively are the hallmarks that stamp AHOPCA as a particularly successful partnership in advancing pediatric oncology in the developing world.

RB1, TP53, and possibly other tumor suppressor genes located at 18q and other localizations are involved in pediatric osteosarcoma carcinogenesis, together with other genetic alterations not fully understood to date. Based on these results, the presence of an altered RB1 gene should be regarded as a poor prognostic factor for pediatric osteosarcoma.