[1] Abe , H.,M.Nakano,A.Nakatsuka,M.Nakayama,M.Koshioka,M.Yamagishi, 2002.Genetic analysis of floral anthocyanin pigmentation traits in Asiatic hybrid lily using molecular linkage maps,Theor. Appl. Genet 105 1175 -1182
[2] Benjamini , Y.,D.Yekutieli, 2001.The control of the false discovery rate in multiple testing under dependency, 29 1165 -1188
[3] Conesa , A,S.Götz,J.M.García-Gómez,J.Terol,M.Talón,M.Robles, 2005.Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research, 21 3674 -3676
[4] Cooley , A.M.,J.H.Willis, 2009.Genetic divergence causes parallel evolution of flower color in Chilean Mimulus, 183 729 -739
[5] Grabherr , M.G.,B.J.Haas,M.Yassour,J.Z.Levin,D.A.Thompson,I.Amit,X.Adiconis,L.Fan,R.Raychowdhury,Q.Zeng,Z.Chen,E.Mauceli,N.Hacohen,A.Gnirke,N.Rhind,Palma , B.W.,Birren , C.,Nusbaum , K.,Lindblad-Toh , N.,Friedman and A.,Regev, 2011.Full-length transcriptome assembly from RNA-Seq data without a reference genome,Nat. Biotechnol 29 644 -652
[6] Huang , J.,X.H.Liu,J.M.Wang,Y.M.Lü, 2014.Transcriptomic analysis of Asiatic lily in the process of vernalization via RNA-seq,Mol. Biol. Rep 41 3839 -3852
[7] Hsu , C.C.,Y.Y.Chen,W.C.Tsai,W.H.Chen,H.H.Chen, 2015.Three R2R3-MYB transcription factors regulate distinct floral pigmentation patterning in Phalaenopsis spp,Plant 168 175 -191
[8] Iseli , C.,C.V.Jongeneel,P.Bucher, 1999.ESTScan: a program for detecting, evaluating, and reconstructing potential coding regions in EST sequences,Proc. Int. Conf. Intelligent Syst. Mol. Biol 99 138 -148
[9] Koes , R.,W.Verweij,F.Quattrocchio, 2005.Flavonoids: a colorful model for the regulation and evolution of biochemical pathways,Trends Plant Sci. 10 236 -242
[10] Lai , Y.S.,Y.Shimoyamada,M.Nakayama,M.Yamagishi, 2012.Pigment accumulation and transcription of LhMYB12 and anthocyanin biosynthesis genes during flower development in the Asiatic hybrid lily (Lilium spp,Plant Sci. 193 136 -147
[11] Li , X.Y.,C.X.Wang,Silva , X.Y.,Liu , X.,Duan , T.L.,Li and H.M.,Sun, 2014.Transcriptome analysis of carbohydrate metabolism during bulblet formation and development in Lilium davidii var. unicolor,BMC Plant Biol 14 358 -
[12] Lim , K.,J.M. van Tuyl,Anderson, 2006.Lilium hybrids,In: Flower Breeding and Genetics: Issues, Challenges and Opportunities for the 21st Century 513 -532
[13] Mortazavi , A.,B.A.Williams,K.McCue,L.Schaeffer,B.Wold, 2008.Mapping and quantifying mammalian transcriptomes by RNA-Seq,Nat. Meth 5 621 -628
[14] Shahin , A.,M.J.Smulders,J.M. van Tuyl , P.,Arens and F.T.,Bakker, 2014.Using multi-locus allelic sequence data to estimate genetic divergence among four Lilium (Liliaceae) cultivars, 5 567 -567
[15] Tamura , K.,D.Peterson,N.Peterson,M.Nei,S.Kumar, 2011.Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods,Mol. Biol 28 2731 -2739
[16] Wang , J.M.,Q.Wang,Y.Yang,X.H.Liu,J.H.Gu,W.Q.Li,S.L.Ma,Y.M.Lu, 2014.De novo assembly and characterization of stress transcriptome and regulatory networks under temperature, salt and hormone stresses in Lilium lancifolium,Mol. Biol. Rep 41 8231 -8245
[17] Winkel-Shirley , B., 2001.Flavonoid biosynthesis. A colorful model for genetics, biochemistry, cell biology, and biotechnology,Plant 126 485 -493
[18] Xu , W.J.,C.Dubos,L.Lepiniec, 2015.Transcriptional control of flavonoid biosynthesis by MYB-bHLH-WDR complexes,Trends Plant Sci. 20 176 -185
[19] Yamagishi , M., 2013.How genes paint lily flowers: Regulation of colouration and pigmentation patterning,Sci. Hortic 163 27 -36
[20] Yamagishi , M.,Y.Shimoyamada,T.Nakatsuka,K.Masuda, 2010.Two R2R3-MYB genes, homologs of petunia AN2, regulate anthocyanin biosynthesis in flower tepals, tepal spots and leaves of Asiatic hybrid lily,Plant Cell Physiol 51 463 -474
[21] Yamagishi , M., S.,Toda and K.,Tasaki, 2014.The novel allele of the LhMYB12 gene is involved in splatter‐ type spot formation on the flower tepals of Asiatic hybrid lilies (Lilium spp, 201 1009 -1020
[22] Ye , J.,L.Fang,H.K.Zheng,Y.Zhang,J.Chen,Z.J.Zhang,J.Wang,S.T.Li,R.Q.Li,L.Bolund,J.Wang, 2006.WEGO: a web tool for plotting GO annotations,Nucl. Acids Res 34 293 -297
[23] Yuan , Y.W.,J.M.Sagawa,L.Frost,J.P.Vela and H.D. Bradshaw, 2014.Transcriptional control of floral anthocyanin pigmentation in monkeyflowers (Mimulus), 204 1013 -1027
[24] Zhang , M.F.,L.M.Jiang,D.M.Zhang,G.X.Jia, .De novo transcriptome characterization of Lilium 'Sorbonne' and key enzymes related to the flavonoid biosynthesis,Mol. Genet 290 399 -412
[25] Zhang , Y., Y. Cheng, H. Ya, S.,Xu, .Transcriptome sequencing of purple petal spot region in tree peony reveals differentially expressed anthocyanin structural genes,Front. Plant Sci. 6 964 -
[26] Zhao , D.,J.Yao,C.Ning,J.S.Meng,S.S.Lin,W.Ding,J.Tao, 2014.Transcriptome sequencing of a chimaera reveals coordinated expression of anthocyanin biosynthetic genes mediating yellow formation in herbaceous peony (Paeonia lactiflora Pall,BMC Genomics 15 219 -224