Among 158 distinct variants, 80 (50.6%) were previously unreported, confirming broad allelic heterogeneity. Eleven patients showed more than one variant. Segregation analysis indicated biallelic disease in five patients, digenic in one, de novo variant with unknown phase in two. Furthermore, our NGS protocol allowed the identification of two patients with somatic mosaicism, which was undetectable with Sanger sequencing. Among patients without PKD1/PKD2 variants, we identified three with possible alternative diagnosis: a patient with biallelic mutations in PKHD1, confirming the overlap between recessive and dominant PKD, and two patients with variants in ALG8 and PRKCSH, respectively. Genotype-phenotype correlations showed that patients with PKD1 variants predicted to truncate (T) the protein experienced end-stage renal disease 9 years earlier than patients with PKD1 non-truncating (NT) mutations and >13 years earlier than patients with PKD2 mutations. ADPKD-PKD1 T cases showed a disease onset significantly earlier than ADPKD-PKD1 NT and ADPK-PKD2, as well as a significant earlier diagnosis. These data emphasize the need to combine clinical information with genetic data to achieve useful prognostic predictions.