Starch is the main storage carbohydrate in plants and an important natural resource for food, feed and industrial raw materials. However, the details regarding the pathway for starch biosynthesis and the diversity of biosynthetic enzymes involved in this process are poorly understood. This study uses a comprehensive phylogenetic analysis of 74 sequenced plant genomes to revisit the evolutionary history of the genes encoding ADP-glucose pyrophosphorylase (AGPase), starch synthase (SS), starch branching enzyme (SBE) and starch de-branching enzyme (DBE). Additionally, the protein structures and expression patterns of these four core genes in starch biosynthesis were studied to determine their functional differences. The results showed that AGPase, SS, SBE and DBE have undergone complicated evolutionary processes in plants and that gene/genome duplications are responsible for the observed differences in isoform numbers. A structure analysis of these proteins suggested that the deletion/mutation of amino acids in some active sites resulted in not only structural variation but also sub-functionalization or neo-functionalization. Expression profiling indicated that AGPase-, SS-, SBE- and DBE-encoding genes exhibit spatio-temporally divergent expression patterns related to the composition of functional complexes in starch biosynthesis. This study provides a comprehensive atlas of the starch biosynthetic pathway, and these data should support future studies aimed at increasing understanding of starch biosynthesis and the functional evolutionary divergence of AGPase, SS, SBE, and DBE in plants.