A total of 53 plant species accessions from different geographic regions, including four melatonin precursor-coding genes obtained from Arachis hypogaea (ASMT1, 2, 3 and T5H) underwent extensive molecular evolutionary analyses. Evolutionary relationships were inferred and showed that dichotomous bifurcating trees did not reflect the true phylogeny since reticulate events took place due likely to recombination. Thus, a phylogenetic network was reconstructed for each type of enzyme and highlighted the presence of such incompatibilities. GARD algorithm pointed out that ASMT1, 2, and 3-coding gene sequences contained recombination sites with significant topological incongruence on both sides of the breakpoints (for ASMT1, and 2), while only on one side of the breakpoints for ASMT3. In contrast, no statistically recombination signal was recorded in T5H-coding gene. Furthermore, gene duplication was localized in the ancestor of a monophyletic group of Populus accessions. Selection pressure was assessed using several statistical models incorporated in HyPhy package through the datamonkey web server. It was demonstrated that numerous individual sites and tree branches experienced predominantly purifying selection. In contrast, the BUSTED model evidenced a gene-wide episodic diversifying selection in the phylogeny of only three enzyme-coding genes (ASMT, and 2, and T5H). Likewise, it was shown that Mixed Effects Model of Episodic Selection (MEME) model detected only episodic positively selected sites in all four melatonin enzymes-coding genes; whereas, REL model failed to detect neither positive nor negative selection in tested individual sites of ASMT3-coding gene.