Myanmar aims to eliminate malaria by 2030. However, recent increase of malaria incidence is a great challenge to archive that goal. Increasing prevalence of Plasmodium vivax also hinders this endeavor. Monitoring genetic structure of the parasite is necessary to understand genetic nature and evolutionary aspect of P. vivax population in Myanmar. Partial fragment flanking blocks I and II of merozoite surface protein-3 alpha of P. vivax (pvmsp-3α) was amplified from P. vivax isolates collected in Pyin Oo Lwin, Mandalay Region, Myanmar in 2013–2015. Sequence analysis of pvmsp-3α was performed to determine genetic diversity and natural selection of this gene. Spatio-temporal genetic changes of pvmsp-3α in Myanmar P. vivax population were also investigated via comparative analysis of gene sequences obtained in this study and previously reported Myanmar pvmsp-3α sequences. Genetic diversity of Myanmar pvmsp-3α was detected in P. vivax isolates analyzed. Size polymorphisms in block I and amino acid changes and recombination events in block II were main factors contributing to the genetic diversity of pvmsp-3α. Comparative spatio-temporal analysis with previously reported Myanmar pvmsp-3α populations revealed the presence of genetic differences by population with moderate genetic differentiation between populations. Similar pattern of natural selection was also detected in Myanmar pvmsp-3α populations. These suggested that enough size of the P. vivax population sufficient to generate or maintain the genetic diversity remains in the population. Thus, continuous molecular surveillance of genetic structure of Myanmar P. vivax is necessary.