Molecular identification and functional verification of SPL9 and SPL15 of Lilium.

The transformation of plants from juveniles to adults is a key process in plant growth and development, and the main regulatory factors are miR156 and SQUAMOSA promoter binding protein-like (SPL) transcription factors. Lilium is an ornamental bulb, but it has a long maturation time. In this experiment, Lilium bulbs were subjected to a temperature treatment of 15 °C for 4 weeks to initiate vegetative phase change. Transmission electron microscopy indicated the cell wall of bud core tissue undergoing vegetative phase change became thinner, the starch grains were reduced, and the growth of the juvenile stage was accelerated. The key transcription factors LbrSPL9 and LbrSPL15 were cloned, and the phylogenetic analysis showed they possessed high homology with other plant SPLs. Subcellular localization and transcription activation experiments confirmed LbrSPL9 and LbrSPL15 were mainly located in the nucleus and exhibited transcriptional activity. The results of in situ hybridization showed the expression levels of LbrSPL9 and LbrSPL15 were increased after temperature change treatment. The functional verification experiment of the transgenic plants confirmed that the overexpression of LbrSPL9 and LbrSPL15 could shorten maturation time. These findings help elucidate the regulatory mechanisms of phase transition in Lilium and provide a reference for breeding research in other bulbous flowers.

Developing an efficient DNA barcoding system to differentiate between Lilium species.

BACKGROUND: Lilium is an important ornamental bulb, possesses medicinal properties, and is also edible. Species within the Lilium genus share very similar morphology and macroscopic characteristics, thus they cannot be easily and clearly distinguished from one another. To date, no efficient species-specific markers have been developed for classifying wild lily species, which poses an issue with further characterizing its medicinal properties. RESULTS: To develop a simple and reliable identification system for Lilium, 45 representative species from 6 sections were used to develop a DNA barcoding system, which was based on DNA sequence polymorphisms. In this study, we assessed five commonly used DNA barcode candidates (ITS, rbcL, ycf1b, matK and psbA-trnH) and five novel barcode candidates obtained from highly variable chloroplast genomic regions (trnL-trnF, trnS-trnG, trnF-ndhJ, trnP-psaJ-rpI33 and psbB-psbH). We showed that a set of three novel DNA barcodes (ITS + trnP-psaJ-rpI33 + psbB-psbH) could be efficiently used as a genetic marker to distinguish between lily species, as assessed by methods including DNAsp, BI and ML tree, and Pair Wise Group (PWG). CONCLUSIONS: A rapid and reliable DNA barcoding method was developed for all 45 wild Lilium species by using ITS, trnP-psaJ-rpI33, and psbB-psbH as DNA barcoding markers. The method can be used in the classification of wild Lilium species, especially endangered species, and also provides an effective method for selective lily breeding.

Metabolome-Based Discrimination Analysis of Five Lilium Bulbs Associated with Differences in Secondary Metabolites.

The bulbs of several Lilium species are considered to be both functional foods and traditional medicine in northern and eastern Asia. Considering the limited information regarding the specific bioactive compounds contributing to the functional properties of these bulbs, we compared the secondary metabolites of ten Lilium bulb samples belonging to five different species, using an ultrahigh-performance liquid chromatography-electrospray ionization-tandem mass spectrometry (UPLC-ESI-MS/MS)-based secondary metabolomics approach. In total, 245 secondary metabolites were detected; further, more metabolites were detected from purple Lilium bulbs (217 compounds) than from white bulbs (123-171 compounds). Similar metabolite profiles were detected in samples within the same species irrespective of where they were collected. By combining herbal analysis and screening differential metabolites, steroid saponins were considered the key bioactive compounds in medicinal lilies. Of the 14 saponins detected, none were accumulated in the bulbs of L. davidii var. willmottiae, also called sweet lily. The purple bulbs of L. regale accumulated more secondary metabolites, and, notably, more phenolic acid compounds and flavonoids. Overall, this study elucidates the differential metabolites in lily bulbs with varying functions and colors and provides a reference for further research on functional foods and the medicinal efficacy of Lilium species.

Characterization of a Cis-Prenyltransferase from Lilium longiflorum Anther.

A group of prenyltransferases catalyze chain elongation of farnesyl diphosphate (FPP) to designated lengths via consecutive condensation reactions with specific numbers of isopentenyl diphosphate (IPP). cis-Prenyltransferases, which catalyze cis-double bond formation during IPP condensation, usually synthesize long-chain products as lipid carriers to mediate peptidoglycan biosynthesis in prokaryotes and protein glycosylation in eukaryotes. Unlike only one or two cis-prenyltransferases in bacteria, yeast, and animals, plants have several cis-prenyltransferases and their functions are less understood. As reported here, a cis-prenyltransferase from Lilium longiflorum anther, named LLA66, was expressed in Saccharomyces cerevisiae and characterized to produce C40/C45 products without the capability to restore the growth defect from Rer2-deletion, although it was phylogenetically categorized as a long-chain enzyme. Our studies suggest that evolutional mutations may occur in the plant cis-prenyltransferase to convert it into a shorter-chain enzyme.

The Chloroplast Genome of Lilium henrici: Genome Structure and Comparative Analysis.

Lilium henrici Franchet, which belongs to the family Liliaceae, is an endangered plant native to China. The wild populations of L. henrici have been largely reduced by habitat degradation or loss. In our study, we determined the whole chloroplast genome sequence for L. henrici and compared its structure with other Lilium (including Nomocharis) species. The chloroplast genome of L. henrici is a circular structure and 152,784 bp in length. The large single copy and small single copy is 82,429 bp and 17,533 bp in size, respectively, and the inverted repeats are 26,411 bp in size. The L. henrici chloroplast genome contains 116 different genes, including 78 protein coding genes, 30 tRNA genes, 4 rRNA genes, and 4 pseudogenes. There were 51 SSRs detected in the L. henrici chloroplast genome sequence. Genic comparison among L. henrici with other Lilium (including Nomocharis) chloroplast genomes shows that the sequence lengths and gene contents show little variation, the only differences being in three pseudogenes. Phylogenetic analysis revealed that N. pardanthina was a sister species to L. henrici. Overall, this study, providing L. henrici genomic resources and the comparative analysis of Lilium chloroplast genomes, will be beneficial for the evolutionary study and phylogenetic reconstruction of the genus Lilium, molecular barcoding in population genetics.

Meiotic chromosome behavior of the male-fertile allotriploid lily cultivar 'Cocossa'.

KEY MESSAGE: Cytological observations of microsporogenesis in the allotriploid lily cultivar 'Cocossa' showed that viable pollen production could be attributed mainly to disoriented spindles, abnormal cytokinesis, and cytomixis during male meiosis. To identify the reasons why the allotriploid lily cultivar 'Cocossa' can produce aneuploid and euploid functional male gametes and can be used as the paternal parent in lily introgression breeding, we performed a detailed investigation of microsporogenesis using the conventional cytological methods. The allotriploid not only produced single pollen grains with variable sizes but also produced adherent pollen grains. Pollen viability was estimated at 50.1% based on staining and 30.8% based on germination. Based on the chromosomal analysis of BC2 plants derived from Oriental cultivars (♀) crossed with the OOT cultivar 'Cocossa' (♂), it was concluded that the objective allotriploid contributed haploid (x), diploid (2x), and aneuploid chromosome complements. Common meiotic abnormalities were observed, indicating the high genetic imbalance of this allotriploid. In addition to normally oriented metaphase II spindles (linear and perpendicular), abnormal spindles, such as parallel, tripolar, fused, and multiple spindles, accounted for 6.21, 6.41, 14.27, and 1.17%, respectively. Tripolar and fused spindles resulted in the production of triads and dyads, which contributed to unreduced pollen production. Some microsporocytes exhibited complete or partial absence of cytokinesis, which led to relatively high frequencies of monads, dyads, and triads. Furthermore, the phenomenon of cytomixis during microsporogenesis occurred mainly in the first meiotic prophase and early development of pollen grains, which we assume is a possible cause of unreduced gamete generation. Our study offers a new resource for lily introgression breeding.

Transcriptome comparison reveals key candidate genes in response to vernalization of Oriental lily.

BACKGROUND: Oriental hybrid lily 'Sorbonne', a very important cut flower for lily, is enjoyed great popularity in the world, but it must experience a period of low winter temperature to initiate or accelerate the flowering process. To gain a better understanding of the temperature signaling pathway and the molecular metabolic reactions involved in the vernalization response, a genome-wide transcriptional analysis using RNA-Seq was performed. RESULTS: 188,447,956 sequencing reads was assembled into 66,327 unigenes and showed similarity to known proteins in the Swiss-Prot protein database, and 2,893, 30,406 and 60,737 unigenes aligned to existing sequences in the KEGG, COG, and GO databases. Based on qRT-PCR results, we studied the expression of three signal regulation pathways genes-the plant hormones signal transduction (LoAP2, LoIAA1, LoARF10), the DNA methylation (LoCMT, LoFLD), and vernalizatin pathway (LoFLC, LoVRN1, LoVRN2, LoFT, LoSOC1, LoLFY, LoSVP) in the immature flower buds of Oriental hybrid lily. In addition, we identified two vernalizaiton-related genes (LoSVP and LoVRN1) from the cDNA library, which appear to be promising candidates for playing key roles in the development and response of flowering in Oriental lily plants, and LoSVP had a function in delaying flowering but LoVRN1could promote flowering early. CONCLUSIONS: We collected a sample for transcriptome sequencing and comparison when the bulb's apical meristem was in the time of floral transition when the apical meristem had not converted into the morphological differentiation process, which helped to obtain more genes playing key roles in the floral induction pathways. The upstream and downstream relationship between different genes were forecasted by the analysis of genes' expression levels in a wide range of time. Future research that is targeted towards how genes interact on each other, which will promote establishing and perfecting the molecular mechanisms of floral induction pathway by vernalization.

Transcriptome Analysis Identifies Key Candidate Genes Mediating Purple Ovary Coloration in Asiatic Hybrid Lilies.

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.

Morphological and ecological divergence of Lilium and Nomocharis within the Hengduan Mountains and Qinghai-Tibetan Plateau may result from habitat specialization and hybridization.

BACKGROUND: Several previous studies have shown that some morphologically distinctive, small genera of vascular plants that are endemic to the Qinghai-Tibetan Plateau and adjacent Hengduan Mountains appear to have unexpected and complex phylogenetic relationships with their putative sisters, which are typically more widespread and more species rich. In particular, the endemic genera may form one or more poorly resolved paraphyletic clades within the sister group despite distinctive morphology. Plausible explanations for this evolutionary and biogeographic pattern include extreme habitat specialization and hybridization. One genus consistent with this pattern is Nomocharis Franchet. Nomocharis comprises 7-15 species bearing showy-flowers that are endemic to the H-D Mountains. Nomocharis has long been treated as sister to Lilium L., which is comprised of more than 120 species distributed throughout the temperate Northern Hemisphere. Although Nomocharis appears morphologically distinctive, recent molecular studies have shown that it is nested within Lilium, from which is exhibits very little sequence divergence. In this study, we have used a dated molecular phylogenetic framework to gain insight into the timing of morphological and ecological divergence in Lilium-Nomocharis and to preliminarily explore possible hybridization events. We accomplished our objectives using dated phylogenies reconstructed from nuclear internal transcribed spacers (ITS) and six chloroplast markers. RESULTS: Our phylogenetic reconstruction revealed several Lilium species nested within a clade of Nomocharis, which evolved ca. 12 million years ago and is itself nested within the rest of Lilium. Flat/open and horizon oriented flowers are ancestral in Nomocharis. Species of Lilium nested within Nomocharis diverged from Nomocharis ca. 6.5 million years ago. These Lilium evolved recurved and campanifolium flowers as well as the nodding habit by at least 3.5 million years ago. Nomocharis and the nested Lilium species had relatively low elevation ancestors (<1000 m) and underwent diversification into new, higher elevational habitats 3.5 and 5.5 million years ago, respectively. Our phylogeny reveals signatures of hybridization including incongruence between the plastid and nuclear gene trees, geographic clustering of the maternal (i.e., plastid) lineages, and divergence ages of the nuclear gene trees consistent with speciation and secondary contact, respectively. CONCLUSIONS: The timing of speciation and ecological and morphological evolutionary events in Nomocharis are temporally consistent with uplift in the Qinghai-Tibetan Plateau and of the Hengduan Mountains 7 and 3-4 million years ago, respectively. Thus, we speculate that the mountain building may have provided new habitats that led to specialization of morphological and ecological features in Nomocharis and the nested Lilium along ecological gradients. Additionally, we suspect that the mountain building may have led to secondary contact events that enabled hybridization in Lilium-Nomocharis. Both the habitat specialization and hybridization have probably played a role in generating the striking morphological differences between Lilium and Nomocharis.

Molecular phylogeny and genetic variation in the genus Lilium native to China based on the internal transcribed spacer sequences of nuclear ribosomal DNA.

We present a comprehensive phylogeny derived from nuclear ribosomal DNA (nrDNA) for 214 samples representing 98 species and five varieties, including 44 species and five varieties native to China. Our collection of 25 species and five varieties (44 samples) covering all five sections of the genus (Comber) distributed in China also were included in the internal transcribed spacer (ITS) database. This study incorporates previous research with an emphasis on Chinese species, including the controversial subsection, Sinomartagon 5c Comber. In the phylogenetic tree obtained by maximum parsimony (PAUP) and maximum likelihood (RAxML) analyses, the samples were divided into four major groups. Our results suggest that the subsection (subsect.) 5c Comber should be classified into the true subsect. 5c and the section (sect.) Lophophorum. And the latter was divided into three subsections (subsect. Lophophorum I, subsect. Lophophorum II, and subsect. Lophophorum III). Based on molecular phylogenetic analysis and fluorescence in situ hybridization, we report that L. henryi and L. rosthornii are closely related, and we propose their classification into subsect. Leucolirion 6a. Our results support Comber's subdivision of sect. Leucolirion, which was primarily based on bulb color. Chinese species were divided into five sections: sect. Martagon, sect. Archelirion, sect. Leucolirion, sect. Sinomartagon, and sect. Lophophorum. These findings contribute to our understanding of the phylogeny, origin, and classification of Lilium.

Analysis of genetic relationships and identification of lily cultivars based on inter-simple sequence repeat markers.

Inter-simple sequence repeat (ISSR) markers were used to discriminate 62 lily cultivars of 5 hybrid series. Eight ISSR primers generated 104 bands in total, which all showed 100% polymorphism, and an average of 13 bands were amplified by each primer. Two software packages, POPGENE 1.32 and NTSYSpc 2.1, were used to analyze the data matrix. Our results showed that the observed number of alleles (NA), effective number of alleles (NE), Nei's genetic diversity (H), and Shannon's information index (I) were 1.9630, 1.4179, 0.2606, and 0.4080, respectively. The highest genetic similarity (0.9601) was observed between the Oriental x Trumpet and Oriental lilies, which indicated that the two hybrids had a close genetic relationship. An unweighted pair-group method with arithmetic means dendrogram showed that the 62 lily cultivars clustered into two discrete groups. The first group included the Oriental and OT cultivars, while the Asiatic, LA, and Longiflorum lilies were placed in the second cluster. The distribution of individuals in the principal component analysis was consistent with the clustering of the dendrogram. Fingerprints of all lily cultivars built from 8 primers could be separated completely. This study confirmed the effect and efficiency of ISSR identification in lily cultivars.

[Molecular identification in genus of Lilium based on DNA barcoding].

To establish a new method for identifying genus of Lilium by DNA barcoding technology, ITS, ITS2, psbA-trnH, matK and rbcL sequences were analyzed in term of variation of inter- and intra-species, barcoding gap, neighbor-joining tree to distinguish genus of Lilium based on 978 sequences from experimental and GenBank database, and identification efficiency was evaluated by Nearest distance and BLAST1 methods. The results showed that DNA barcoding could identify different species in genus of Lilium. ITS sequence performed higher identification efficiency, and had significant difference between intra- and inter-species. And NJ tree could also divide species into different clades. Results indicate that DNA barcoding can identify genus of Lilium accurately. ITS sequence can be the optimal barcode to identify species of Lilium.

[Genetic background of lily germplasm resources based on SSR markers].

To study the genetic background of lily (Lilium spp.) germplasm resources, and accurately evaluate and select excellent germplasm for genetic improvement of lily, we analyzed the genetic background of 62 lily germplasm accessions from 11 provinces of China by using simple sequence repeat (SSR) molecular markers. The results showed that 15 out of 83 pairs of lily SSR primers were polymorphic. A total of 157 allelic loci were amplified, with the number of alleles per locus ranging from 5 to 19 and the average number of effective alleles per locus being 4.162 8. The average observed heterozygosity and expected heterozygosity were 0.228 2 and 0.694 1, respectively. The average polymorphic information content was 0.678 8. The average Nei's diversity index and Shannon's information index were 0.694 1 and 1.594 9, respectively, indicating that the tested lily germplasm had high genetic diversity. The 62 germplasm accessions were classified into 5 groups by the unweighted pair group method with arithmetic mean (UPGMA) and into 3 groups by the principal component analysis. The two analyses revealed a geographic correlation among different groups. The majority of lily germplasm accessions from the same source tended to cluster together. The population structure analysis classified the lily accessions into 4 populations and 1 mixed population. The above results provide a theoretical basis and genetic resources for the precise identification and breeding of lily germplasm resources.

Cytogenetic studies on meiotic chromosome behaviors in sterile Oriental x Trumpet lily.

In order to determine the reasons for pollen sterility in lily hybrids, four diploid sterile Oriental x Trumpet (OT) lily cultivars ('Nymph', 'Gluhwein', 'Yelloween', and 'Shocking') were used to investigate the meiotic chromosome behaviors in pollen mother cells (PMCs), using genomic in situ hybridization and conventional cytological methods. At metaphase I, chromosome associations were quite variable, not only among different genotypes but also in different PMCs of the same genotype. In addition to bivalents, a certain amount of univalent, trivalents, and quadrivalents were observed in all of the investigated genotypes. In addition, ring octavalents and ring hexavalents were observed in 'Nymph'. Even dodecavalents were observed in 'Nymph'. These abnormal chromosome associations at metaphase I implied the occurrence of chromosome interchanges (translocation) in these intersectional hybrids. At anaphase-telophase, a large number of laggard chromosomes and different kinds of chromosome bridge configurations were observed. At the tetrad stage, micronuclei and polyads were also found in many PMCs. All of these abnormal chromosome behaviors in PMCs were responsible for the pollen sterility in lily hybrids.

Rapid identification of Lilium species and polysaccharide contents based on near infrared spectroscopy and weighted partial least square method.

Polysaccharide is the main active compound of Lilium, and showed many activities, such as hypoglycemic, antioxidant, immune-modulatory. There are three types' Lilium in China market, i.e. Lilium lancifolium Thunb (JD), Lilium davidiivar. Unicolor Salisb (L. davidii var)(LZBH), and Lilium brownii F.E. Brown var. viridulum Baker (BH). Near infrared spectroscopy (NIR) technique has become popular in the fields of quality control, due to its advantages, such as fast, non-destructive, and can detect several ingredients, simultaneously. In this study, a classification model was established based on NIR technique and random forest method to accurately distinguish three types' Lilium species, and the classification accuracy reached 94.37%. Furthermore, taking the effects of neighbor wavelength into account, a new weighted partial least square algorithm was proposed to establish an accurate and quantitative model for predicting the polysaccharide contents of these samples. In the model establishing process, some signal pre-treatment methods were optimized, and the validation results with highest determination coefficient (R2) and low root mean square errors of prediction (RMSEP) were, 0.9455 and 0.9098, respectively. The obtained results showed that combined NIR technique with chemometrics was an effective and green method for quality control.

Overexpression of LiTPS2 from a cultivar of lily (Lilium 'Siberia') enhances the monoterpenoids content in tobacco flowers.

Lily, a famous cut flower with highly fragrance, has high ornamental and economic values. Monoterpenes are the main components contributing to its fragrance, and terpene synthase (TPS) genes play critical roles in the biosynthesis of monoterpenoids. To understand the function of TPS and to explore the molecular mechanism of floral scent in cultivar Lilium 'Siberia', transcriptomes of petal at different flowering stages and leaf were obtained by RNA sequencing and three unigenes related to TPS genes were selected for further validation. Quantitative real-time PCR showed that the expression level of LiTPS2 was greater than that of the other two TPS genes. Phylogenetic analysis indicated that LiTPS2 belonged to the TPSb subfamily, which was responsible for monoterpenes synthesis. Subcellular localization demonstrated that LiTPS2 was located in the chloroplasts. Furthermore, functional characterization showed that LiTPS2 utilized both geranyl pyrophosphate (GPP) and farnesyl pyrophosphate (FPP) to produce monoterpenoids such as linalool and sesquiterpenes like trans-nerolidol, respectively. Ectopic expression in transgenic tobacco plants suggested that the amount of linalool from the flowers of transgenic plants was 2-3 fold higher than that of wild-type plants. And the emissions of myrcene and (E)-ß-ocimene were also accumulated from the flowers of LiTPS2 transgenic lines. Surprisingly, these three compounds were the main fragrance components of oriental lily hybrids. Our results indicated that LiTPS2 contributed to the production of monoterpenes and could effectively regulate the aroma of Lilium cultivars, laying the foundation for biotechnological modification of floral scent profiles.

Phylogenetic background of Orange lily (Lilium bulbiferum s.l.) cultivars from a genetically isolated environment.

Domesticated Orange lily ( LILIUM BULBIFERUM s.l.) cultivars do not typically produce seeds, and Orange lily is not native to Finland. Therefore, back crossing of the cultivars with wild species has not been possible. Genetic variability and genuineness of eight Finnish traditionally-grown Orange lily cultivars was studied. RAPD patterns were compared between the cultivars and genuine Orange lily ( LILIUM BULBIFERUM L.), and a related Dauricum group. The results showed partition of tested genotypes into four groups, L. CANADENSE as the outgroup. The cultivars were divided into two subgroups where the trait to form bulbils was characteristic to subgroup I. The cultivated strains differed from each other as much as from the seedling strains, but were genetically closer to genuine Orange lily than the Dauricum group. This indicates that the cultivars are genuine forms of Orange lily species. The special morphological features of the cultivars have likely been formed during centuries-long genetic isolation from natural populations.

Frequent gene flow blurred taxonomic boundaries of sections in Lilium L. (Liliaceae).

Gene flow between species may last a long time in plants. Reticulation inevitably causes difficulties in phylogenetic reconstruction. In this study, we looked into the genetic divergence and phylogeny of 20 Lilium species based on multilocus analyses of 8 genes of chloroplast DNA (cpDNA), the internally transcribed nuclear ribosomal DNA (nrITS) spacer and 20 loci extracted from the expressed sequence tag (EST) libraries of L. longiflorum Thunb. and L. formosanum Wallace. The phylogeny based on the combined data of the maternally inherited cpDNA and nrITS was largely consistent with the taxonomy of Lilium sections. This phylogeny was deemed the hypothetical species tree and uncovered three groups, i.e., Cluster A consisting of 4 taxa from the sections Pseudolirium and Liriotypus, Cluster B consisting of the 4 taxa from the sections Leucolirion, Archelirion and Daurolirion, and Cluster C comprising 10 taxa mostly from the sections Martagon and Sinomartagon. In contrast, systematic inconsistency occurred across the EST loci, with up to 19 genes (95%) displaying tree topologies deviating from the hypothetical species tree. The phylogenetic incongruence was likely attributable to the frequent genetic exchanges between species/sections, as indicated by the high levels of genetic recombination and the IMa analyses with the EST loci. Nevertheless, multilocus analysis could provide complementary information among the loci on the species split and the extent of gene flow between the species. In conclusion, this study not only detected frequent gene flow among Lilium sections that resulted in phylogenetic incongruence but also reconstructed a hypothetical species tree that gave insights into the nature of the complex relationships among Lilium species.

Chloroplast genomes of Lilium lancifolium, L. amabile, L. callosum, and L. philadelphicum: Molecular characterization and their use in phylogenetic analysis in the genus Lilium and other allied genera in the order Liliales.

Chloroplast (cp) genomes of Lilium amabile, L. callosum, L. lancifolium, and L. philadelphicum were fully sequenced. Using these four novel cp genome sequences and five other previously sequenced cp genomes, features of the cp genomes were characterized in detail among species in the genus Lilium and other related genera in the order Liliales. The lengths and nucleotide composition showed little variation. No structural variation was found among the cp genomes in Liliales. Gene contents were conserved among four newly sequenced cp genome in Lilium species, the only differences being in two pseudogenes. We identified 112 genes in 13 functional categories, 18 of which carried introns that were conserved among the species in Liliales. There were 16-21 SSR loci (>12 bp, >3 repeats) in the cp genomes in Lilium and the genomic locations of these loci were highly variable among the species. Average mutations were 15 SNPs per 1kb and 5 indels per 1kb, respectively, in the cp genomes of the newly sequenced four Lilium species. Phylogenetic classifications revealed some discrepancies between trees based on the cp genomes and previous classifications based on the morphology and geographic distributions.

Complete chloroplast genome sequences of Lilium: insights into evolutionary dynamics and phylogenetic analyses.

Lilium is a large genus that includes approximately 110 species distributed throughout cold and temperate regions of the Northern Hemisphere. The species-level phylogeny of Lilium remains unclear; previous studies have found universal markers but insufficient phylogenetic signals. In this study, we present the use of complete chloroplast genomes to explore the phylogeny of this genus. We sequenced nine Lilium chloroplast genomes and retrieved seven published chloroplast genomes for comparative and phylogenetic analyses. The genomes ranged from 151,655 bp to 153,235 bp in length and had a typical quadripartite structure with a conserved genome arrangement and moderate divergence. A comparison of sixteen Lilium chloroplast genomes revealed ten mutation hotspots. Single nucleotide polymorphisms (SNPs) for any two Lilium chloroplast genomes ranged from 8 to 1,178 and provided robust data for phylogeny. Except for some of the shortest internodes, phylogenetic relationships of the Lilium species inferred from the chloroplast genome obtained high support, indicating that chloroplast genome data will be useful to help resolve the deeper branches of phylogeny.