In the case of Newcastle disease virus, multiple factors such as host adaptation, immune response evasion, and selective pressures have been suggested to result in evolution of viruses and the emergence of genetic variants. Multiple studies on virus classification and global epidemiological links have yielded consistent data. Here, we have performed a molecular analysis study of circulating Newcastle disease viruses in Iran (1995-2016). According to evolutionary divergences, subgenotype VIg, VIj, VIIj, VIId, XIIIa and XIIId isolates have been circulating in the country during a 21-year period. Based on data analysis, VIg isolates shared highest sequence identity with Russian and Polish isolates of the VIg subgenotype, while VIj subgenotype isolates (2012) were most similar to a virus isolated in 2015 in India. Analysis of the evolutionary divergence of subgenotype VIIj suggests that Chinese and Ukrainian viruses may have played a crucial role in the emergence of VIIj isolates. Evolutionary difference studies also indicated that XIIIa isolates circulating in Iran may have caused the emergence of adapted variants of subgenotype XIIId. Therefore, we propose that the evolutionary and epidemiological study of virulent Newcastle disease viruses could help to provide accurate molecular data about variants circulating in the region, thus aiding in the design of more efficient recombinant vaccines.