Gallbladder cancer (GBC) is a highly aggressive malignancy of the biliary tract. Most cases of GBC are diagnosed in low-income and middle-income countries, and research into this disease has long been limited. In this study we therefore investigate the epigenetic changes along the model of GBC carcinogenesis represented by the sequence gallstone disease → dysplasia → GBC in Chile, the country with the highest incidence of GBC worldwide. To perform epigenome-wide methylation profiling, genomic DNA extracted from sections of formalin-fixed, paraffin-embedded gallbladder tissue was analyzed using Illumina Infinium MethylationEPIC BeadChips. Preprocessed, quality-controlled data from 82 samples (gallstones n=32, low-grade dysplasia n=13, high-grade dysplasia n=9, GBC n=28) were available to identify differentially methylated markers, regions, and pathways as well as changes in copy number variations (CNVs). The number and magnitude of epigenetic changes increased with disease development and predominantly involved the hypermethylation of cytosine-guanine dinucleotide islands and gene promoter regions. The methylation of genes implicated in Wnt signaling, Hedgehog signaling, and tumor suppression increased with tumor grade. CNVs also increased with GBC development and affected cyclin-dependent kinase inhibitor 2A, MDM2 proto-oncogene, tumor protein P53, and cyclin D1 genes. Gains in the targetable Erb-B2 receptor tyrosine kinase 2 gene were detected in 14% of GBC samples. Our results indicate that GBC carcinogenesis comprises three main methylation stages: early (gallstone disease and low-grade dysplasia), intermediate (high-grade dysplasia), and late (GBC). The identified gradual changes in methylation and CNVs may help to enhance our understanding of the mechanisms underlying this aggressive disease and eventually lead to improved treatment and early diagnosis of GBC.