The LONP1 gene encodes Lon protease, which is responsible for degrading damaged or misfolded proteins and binding mitochondrial DNA. Previously, LONP1 variants have been identified in patients with cerebral, ocular, dental, auricular, and skeletal anomalies (CODAS syndrome) and mitochondrial diseases. Seizures were occasionally observed. However, the association between LONP1 and epilepsy remains elusive. In this study, we performed trio-based whole-exome sequencing in a cohort of 450 patients with unexplained epilepsy and identified four pairs of compound heterozygous LONP1 variants in four unrelated cases. All patients exhibited good responses to anti-seizure medications and demonstrated no developmental delay or intellectual disabilities. The variant allele frequencies observed in this study were absent or low in the general population and were significantly lower than those of benign variants. At least one variant in each biallelic pair affected hydrogen bonding and/or altered protein stability. The CODAS syndrome-associated variants were concentrated in the AAA+ module, especially the α domain. Four of the five mitochondrial disease-associated variants were located in the AAA + domain and the NTD5H and NTD3H subdomains. In contrast, each of the biallelic variants from the patients with pure epilepsy had one variant located in the linker domain, and the other variant located in the mitochondrial targeting sequence or P domain. This study suggested that LONP1 gene is potentially a novel candidate gene for pure epilepsy. The phenotypic variation is associated with the sub-regional effects of variants.