BackgroundAglaonema commutatum ‘Red Valentine’, as a foliage ornamental plant, is widely used for interior and exterior decoration because of its easy cultivation and management. However, reduced proportion of red foliage during large-scale production of A. commutatum seedlings is a frequent occurrence, which has considerable implications on the plant’s ornamental and market value. However, the molecular mechanisms underlying this phenomenon remain unclear.ResultsTo explore the molecular basis of the variation in leaf color of A. commutatum Red Valentine, we performed transcriptome sequencing with the Illumina platform using two different varieties of A. commutatum, namely Red Valentine and a green mutant, at three different stages of leaf development. We annotated 63,621 unigenes and 14,186 differentially expressed genes by pairwise comparison. Furthermore, we identified 26 anthocyanin biosynthesis structural genes. The transcript per million (TPM) values were significantly higher for Red Valentine than for the green mutant in all three developmental stages, consistent with the high anthocyanin content of Red Valentine leaves. We detected positive transcription factors that may be involved in the regulation of anthocyanin biosynthesis using BLAST and through correlation analysis. Downregulation of these transcription factors may downregulate the expression of anthocyanin genes. We obtained full-length cDNA of the anthocyanin biosynthesis and regulatory genes and constructed phylogenetic trees to ensure accuracy of the analysis.ConclusionsOur study provides insights into the molecular mechanisms underlying leaf variation in A. commutatum Red Valentine and may be used to facilitate the breeding of ornamental cultivars with high anthocyanin levels.