Lily tepals have a short lifespan. Once the tepals senesce, the ornamental value of the flower is lost. Some cultivars have attractive purple ovaries and fruits which greatly enhance the ornamental value of Asiatic hybrid lilies. However, little is known about the molecular mechanisms of anthocyanin biosynthesis in Asiatic hybrid lily ovaries. To investigate the transcriptional network that governs purple ovary coloration in Asiatic hybrid lilies, we obtained transcriptome data from green ovaries (S1) and purple ovaries (S2) of Asiatic “Tiny Padhye”. Comparative transcriptome analysis revealed 4228 differentially expressed genes. Differential expression analysis revealed that ten unigenes including four CHS genes, one CHI gene, one F3H gene, one F3′H gene, one DFR gene, one UFGT gene, and one 3RT gene were significantly up-regulated in purple ovaries. One MYB gene, LhMYB12-Lat, was identified as a key transcription factor determining the distribution of anthocyanins in Asiatic hybrid lily ovaries. Further qPCR results showed unigenes related to anthocyanin biosynthesis were highly expressed in purple ovaries of three purple-ovaried Asiatic hybrid lilies at stages 2 and 3, while they showed an extremely low level of expression in ovaries of three green-ovaried Asiatic hybrid lilies during all developmental stages. In addition, shading treatment significantly decreased pigment accumulation by suppressing the expression of several unigenes related to anthocyanin biosynthesis in ovaries of Asiatic “Tiny Padhye”. Lastly, a total of 15,048 Simple Sequence Repeats (SSRs) were identified in 13,710 sequences, and primer pairs for SSRs were designed. The results could further our understanding of the molecular mechanisms of anthocyanin biosynthesis in Asiatic hybrid lily ovaries.

Previous studies have suggested that multidrug and toxic compound extrusion (MATE) proteins might be involved in flavonoid transportation. However, whether MATE proteins are involved in anthocyanin accumulation in Lilium is unclear. Here, a flavonoid transport-related MATE candidate gene, LhDTX35, was cloned from the Asiatic hybrid lily cultivar ‘Tiny Padhye’ by rapid amplification of 5’ and 3’ cDNA ends (RACE) and found to encode 507 amino acids. BLASTx results indicated that LhDTX35 showed high homology to the DTX35 genes of other species. Bioinformatics analysis predicted that the protein encoded by LhDTX35 possessed 12 typical transmembrane segments and had functional domains typical of the MATE-like superfamily. Phylogenetic analysis grouped LhDTX35 in the same clade as the DTX35 of other species. Notably, the expression pattern of LhDTX35 was positively correlated with floral anthocyanin accumulation in ‘Tiny Padhye’. A subcellular localization assay showed that the protein encoded by LhDTX35 was plasmalemma localized but not nuclear, indicating that the LhDTX35 gene may function as a carrier protein to transport anthocyanins in Lilium. Functional complementation of the Arabidopsis DTX35 gene demonstrated that LhDTX35 could restore silique-infertility and the anthocyaninless phenotype of an Arabidopsis DTX35 mutant. These results indicated that LhDTX35 might be involved in anthocyanin accumulation in Lilium.