Abstract
Lilium lankongense Franchet is a lily species found on the Qinghai-Tibet Plateau. It is pink with deep red spots, has a high ornamental value, and is used in hybrid breeding of horticultural lily varieties. We have insufficient knowledge of the genetic resources of L. lankongense and its phylogenetic relationships with related species. Recent molecular phylogenetic studies have shown a very close phylogenetic relationship between L. lankongense and the five species L. duchartrei, L. stewartianum, L. matangense, L. lophophorum, and L. nanum. However, molecular markers still lack sufficient signals for population-level research of the genus Lilium. We sequenced and compared the complete plastid sequences of L. lankongense and its five related species. The genomes ranged from 152,307 bp to 152,611 bp. There was a slight inconsistency detected in inverted repeat and single copy boundaries and there were 53 to 63 simple sequence repeats in the six species. Two of the 12 highly variable regions (trnC-petN and rpl32-trnL) were verified in 11 individuals and are promising for population-level studies. We used the complete sequence of 33 plastid genomes, the protein-coding region sequence, and the nuclear ITS sequence to reconstruct the phylogenetic tree of Lilium species. Our results showed that the plastid gene tree and nuclear gene tree were not completely congruent, which may be caused by hybridization, insufficient information contained in the nuclear ITS, or the small number of samples. The results of phylogenetic analysis based on plastid genomes indicated that the six Lilium species were closely related. Our study provides a preliminarily rebuilt backbone phylogeny that is significant for future molecular and morphological studies of Lilium.