Preemptive targeted pharmacogenetic testing of candidate variations in DPYD is currently being used to limit toxicity associated with fluoropyrimidines. The use of innovative next generation sequencing (NGS) approaches could unveil additional rare (minor allele frequency <1%) genetic risk variants. However, their predictive value and management in clinical practice are still controversial, at least partly due to the challenges associated with functional analyses of rare variants. The aim of this study was to define the predictive power of rare DPYD variants burden on the risk of severe fluoropyrimidine-related toxicity.The DPYD coding sequence and untranslated regions were analyzed by NGS in 120 patients developing grade 3–5 (NCI-CTC vs3.0) fluoropyrimidine-related toxicity and 104 matched controls (no-toxicity). The functional impact of rare variants was assessed using two different in silico predictive tools (i.e., Predict2SNP and ADME Prediction Framework) and structural modeling. Plasma concentrations of uracil (U) and dihydrouracil (UH2) were quantified in carriers of the novel variants.Here, we demonstrate that the burden of rare variants was significantly higher in patients with toxicity compared to controls (p = 0.007, Mann-Whitney test). Carriers of at least one rare missense DPYD variant had a 16-fold increased risk in the first cycle and an 11-fold increased risk during the entire course of chemotherapy of developing a severe adverse event compared to controls (p = 0.013 and p = 0.0250, respectively by multinomial regression model). Quantification of plasmatic U/UH2 metabolites and in silico visualization of the encoded protein were consistent with the predicted functional effect for the novel variations.Analysis and consideration of rare variants by DPYD-sequencing could improve prevention of severe toxicity of fluoropyrimidines and improve patients’ quality of life.