Lymphogranuloma venereum (LGV) is already endemic in vulnerable populations in several European countries; however, molecular epidemiology data with improved accuracy are necessary to better understand LGV epidemic in these countries. Current strategies to study the molecular epidemiology of LGV cases involve schemes based on a few genetic fragments of Chlamydia trachomatis, which have demonstrated limited discriminatory power for LGV. Therefore, this study aimed to propose a new combination of molecular markers based on the most variable genes of L-genotype genomes to improve the characterization of the current LGV epidemic in Madrid, Spain. Four genes were selected according to their diversity index (CTLon_0054, CTLon_0087, CTLon_0243 and CTLon_0301) for use in combination with ompA. In silico and experimental studies were performed to compare the previously described multilocus sequence typing (MLST) schemes with our proposal. Moreover, the proposed scheme was applied (n = 68) to analyze the spatio-temporal spread of the LGV cases. Our proposal demonstrated higher diversity allowing the identification of three main groups compared to the previously published MLST based on hypervariable genes wherein only a single sequence type was identified. The temporal analysis showed that the major cluster was progressively diversifying, revealing a very active transmission chain. Furthermore, an L2b genome identical to that of the origin of the epidemic was detected, suggesting reintroductions or a low screening rate in vulnerable populations. The spatial distribution suggests that the selection and spread of new variants occurs from the central district to the peripheral regions. The scheme proposed in this study has proven to be useful for appropriate discrimination of LGV strains. This study, to our knowledge for the first time, demonstrates a spatio-temporal spread that increases our understanding and identifies areas with special susceptibility for maintenance of the endemic situation of LGV.