Abstract Heterozygous variants in the glucocerebrosidase GBA1 gene are an increasingly recognized risk factor for Parkinson’s disease (PD). Due to the GBAP1 pseudogene, which shares 96% sequence homology with the GBA1 coding region, accurate variant calling by array-based or short-read sequencing methods remains a major challenge in understanding the genetic landscape of GBA1-associated PD. We analyzed 660 patients with PD, 100 patients with Parkinsonism and 808 healthy controls from the Luxembourg Parkinson’s study, sequenced using amplicon-based long-read DNA sequencing technology. We found that 12.1% (77/637) of PD patients carried GBA1 variants, with 10.5% (67/637) of them carrying known pathogenic variants (including severe, mild, risk variants). In comparison, 5% (34/675) of the healthy controls carried GBA1 variants, and among them, 4.3% (29/675) were identified as pathogenic variant carriers. We found four GBA1 variants in patients with atypical parkinsonism. Pathogenic GBA1 variants were 2.6-fold more frequently observed in PD patients compared to controls (OR = 2.6; CI = [1.6,4.1]). Three novel variants of unknown significance (VUS) were identified. Using a structure-based approach, we defined a potential risk prediction method for VUS. This study describes the full landscape of GBA1-related parkinsonism in Luxembourg, showing a high prevalence of GBA1 variants as the major genetic risk for PD. Although the long-read DNA sequencing technique used in our study may be limited in its effectiveness to detect potential structural variants, our approach provides an important advancement for highly accurate GBA1 variant calling, which is essential for providing access to emerging causative therapies for GBA1 carriers.