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.

BackgroundSpinocerebellar ataxias (SCAs) are a group of clinically and genetically diverse and autosomal-dominant disorders characterised by neurological deficits in the cerebellum. At present, there is no cure for SCAs. Of the different distinct subtypes of autosomal-dominant SCAs identified to date, causative genes for only a fraction of them are currently known. In this study, we investigated the cause of an autosomal-dominant SCA phenotype in a family that exhibits cerebellar ataxia and pontocerebellar atrophy along with a global reduction in brain volume.Methods and resultsWhole-exome analysis revealed a missense mutation c.G1391A (p.R464H) in the coding region of the coiled-coil domain containing 88C (CCDC88C) gene in all affected individuals. Functional studies showed that the mutant form of CCDC88C activates the c-Jun N-terminal kinase (JNK) pathway, induces caspase 3 cleavage and triggers apoptosis.ConclusionsThis study expands our understanding of the cause of autosomal-dominant SCAs, a group of heterogeneous congenital neurological conditions in humans, and unveils a link between the JNK stress pathway and cerebellar atrophy.

Supplemental Figure.1. 6-MP Dose Intensity in Patients with both ABCC4 c.912G>T and ITPA c.94C>A variants. 6-MP dose intensity were compared in patients with (n=11) and without both variants (n=134) using Mann Whitney U test. "+" indicates mean values.

No abstract

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 novel genetic cause of thiopurine toxicity, and NUDT15-informed preemptive dose reduction is quickly adopted in clinical settings. To exhaustively identify pharmacogenetic variants in this gene, we developed massively parallel NUDT15 function assays to determine 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 aminoacid resolution. Functional effects for NUDT15 variants accurately predicted toxicity risk alleles in 2,398 patients treated with thiopurines, with 100% sensitivity and specificity, in contrast with poor performance of bioinformatic prediction algorithms. In conclusion, our massively parallel variant function assays identified 1,103 deleterious NUDT15 variants, providing a comprehensive reference of variant function and vastly improving the ability to implement pharmacogeneticsguided thiopurine treatment individualization.