Lilium spp. pollen in China (Liliaceae): Taxonomic and Phylogenetic Implications and Pollen Evolution Related to Environmental Conditions

Abstract

Recent molecular and karyologic studies have significantly modified delimitation of Lilium. However, despite the importance of pollen evolution in the genus comprehensive studies with electron microscopy and evaluation of pollen evolution are lacking. Therefore, we studied pollen morphology in a sample of 65 individuals from 37 taxa covering all the sections distributed in the world, using scanning electron microscopy. Our collection of 49 individuals from 21 taxa covering all five sections in China was also included in the database. We found pollen tetrads in L. bakerianum. Based on present and previous studies, our results suggest that pollen from L. formosanum should be classified as a new type, Formosanum. Combined with morphological and molecular evidence, pollen sculpture patterns appear to reflect phylogenetic relationships and are useful for species or subsection delimitation. Based on a comprehensive survey and correlation with potential functional implications, we propose the following hypothesis: evolution of an exine sculpture shows pollen type trends from Martagon → Callose → Concolor → Formosanum. The evolutionary trend regarding pollen sculpture and size could be related to selective pressure to adapt to environmental conditions. Pollen size and shape showed a significantly positive correlation with annual precipitation, and smaller pollen grains appear to adapt better in habitats with extreme conditions. Evolution trends in exine sculpture do not appear to be definitively correlated with pollen size and shape.

About the Authors:

Yun-peng Du

Affiliation: College of Landscape Architecture, National Engineering Research Center for Flowers, Beijing Forestry University, Beijing, China

Chi Wei

Affiliation: College of Landscape Architecture, National Engineering Research Center for Flowers, Beijing Forestry University, Beijing, China

Zhong-xuan Wang

Affiliation: College of Landscape Architecture, National Engineering Research Center for Flowers, Beijing Forestry University, Beijing, China

Shuang Li

Affiliation: College of Landscape Architecture, National Engineering Research Center for Flowers, Beijing Forestry University, Beijing, China

Heng-bin He

Affiliation: College of Landscape Architecture, National Engineering Research Center for Flowers, Beijing Forestry University, Beijing, China

Gui-xia Jia

* E-mail: [email protected]

Affiliation: College of Landscape Architecture, National Engineering Research Center for Flowers, Beijing Forestry University, Beijing, China

Introduction

Approximately 110 to 115 Lilium species are distributed in the cold and temperate regions of the Northern Hemisphere [1], [2], particularly in East Asia, the Himalayas and Hengduan Mountains, North America and Europe. A total of 55 species occur in China [2]. De Jong [3] and Patterson and Givnish [4] consider southwest China and the Himalayas to be the center of origin of this genus.

Classification of this genus has been historically complicated. Several classifications for Lilium have been proposed based on morphological characters. Detailed studies have been performed for East Asia, European and North American Lilium species [5]–[7]. Based on 13 morphological characters and two germination types, Comber divided this genus into the following seven sections: Martagon Rchb., Pseudolirium Endl. which is limited to North America, Liriotypus Asch. and Graeb. which is distributed across Europe and the Caucasus, Archelirion Baker, Sinomartagon Comber, Leucolirion Wilson, and Daurolirion Comber, representing the most widely accepted taxonomical divisions [8]. Wang and Tang recognized sect. Lophophorum (Bur. et Franch.) Wang et Tang out of sect. Sinomartagon Comber and included campaniform-flowered species [9]. Liang [10] and Haw [11] modified sect. Lophophorum to accommodate the Nomocharis-like Lilium species in sect. Sinomartagon Comber. Chinese species were divided into five sections: Martagon, Archelirion, Sinomartagon, Leucolirion and Lophophorum.

Recently, molecular phylogenetic analyses and chromosome techniques have improved the understanding of several groups within the genus and modified the phylogenetic position of Comber’s classification, such as placement of sect. Daurolirion Comber in sect. Sinomartagon, L. henryi in subsect. Leucolirion 6b, modification of sect. Lophophorum and relationship confirmation in sect. Liriotypus [12]–[13], [16]–[21]. Preliminary research found that sect. Sinomartagon, which mainly occurred in China, was complicated and polyphyletic. As indicated by Patterson and Givnish [4], intercontinental dispersal details of the genus Lilium are not yet clear. The division of subsect. Sinomartagon 5c and sect. Lophophorum are controversial and will require further research. Nishikawa et al. [18], [19] suggested that L. henryi be classified into subsect. Leucolirion 6a and that it showed similar morphological features with L. rosthornii, thereby demonstrating that the phylogenetic position of L. rosthornii needs further study.

There are few relevant studies regarding pollen morphology which defines the taxonomic and reflect the evolution of the genus. According to the description within Lilium by Baranova [22] based on the number, shape and arrangement of columellae that form the muri, there are three morphological types of pollen: (1) Martagon (muri formed by rectangular columellae); (2) Callose (muri formed by rounded columellae); and (3) Concolor (muri formed by separated rounded and polygonal columellae). Previous studies found that most Lilium species have single pollen grains. However, pollen tetrads were found in L. sempervivoideum H. Lév. and L. amoenum E. H. Wilson ex Sealy, and the size and sculptural elements confirmed the taxa as two subspecies in L. sempervivoideum [23]. The pollen morphology of L. lophophorum (Bur. et Franch.) Franch., L. henrici Franch., L. souliei (Franch.) Sealy and L. nanum Klotz. et Garcke supported placement in Lilium, which differs from Nomocharis in apture and sculptural elements, and showed an evolutionary aperture trend from monocolpate to porate [24]. Muratović et al. [25] showed that two related European species, L. bosniacum and L. carniolicum, share similar pollen morphology. In addition, pollen morphology of some Chinese species under scanning electron microscopy (SEM), including 9 species described by Li and Qin [26], 10 species and 3 cultivars by Zhang et al. [27], and 12 species and 6 cultivars by Wu et al. [28], could provide taxonomic implications within Lilium:pollen has not only the commonness of genus, but also the specificity of single species. Interspecific pollen size and morphological characteristics has some difference, which has a certain reference value for the classification of Lilium. For example, there are differences in pollen size among L. cernnum, L. lanciflium and L. pumilum. Also there are significant differences in pollen ornamentation and morphology between L. leucanthum and its variety L. leucanthum var. centifolium from Qinling Mountains. Results obtained by Wang et al. [29] and Liu [30] indicate that pollen size parameters of L. pumilum and L. concolor from different provenances showed different degrees of variation. Determining whether this is a universal phenomenon in other species or if it is an important characteristic of species from different provenances requires a large amount of data. Therefore, a comprehensive study of Lilium pollen morphology in China, particularly a comparative study between populations or provenances using SEM, was needed. Environmental constraints act on patterns of differentiation, as all character states are the consequence of interactions between phylogenetic and environmental constraints [31]. Pollen could be considered a functional unit, with the exine structure as a compromise between four main functions [32]: protective, harmomegathic (the ability to absorb bending stresses, as occur during desiccation), reservoir (the role in producing an adhesive surface or as recognition substances) and clustering. Based on a comparative survey in the Rhaponticum group (Asteraceae, Centaureinae), Hidalgo et al. [33] hypothesized that occurrence in habitats with extreme conditions appears to be the main factor in the evolution from one form of pollen to another.

By combining previous and present results, our aim is to provide palynology information regarding the extent of pollen morphological diversity within and between species/sections, link this information with current phylogenetic data, and contribute to a comprehensive understanding of the systematics and evolution of Lilium by: (1) verifying pollen types in Chinese species; (2) discussing these findings within a phylogenetic framework; and (3) correlating pollen morphology and size with habitat.

Materials and Methods

Ethics Statement

For conservation reasons, we collected only small amounts of Lilium pollen. Our activities did not have any adverse effects of Lilium populations, and permission for collecting was obtained from National Engineering Research Center for Flowers of China (No.35 Tsinghua East Road Haidian District, Beijing, P.R. China).

Table 1 shows the origin of all species studied which were collected fresh. Our observations for 21 taxa represented by 49 samples were also included in the database. Division into the sections was based on studies by Comber [8], Gao et al. [13], Nishikawa et al. [18], [19] and Du et al. [16]) as follows: two species and one variety of sect. Martagon, nine species and three varieties of sect. Sinomartagon including L. dauricum Ker-Gawl (sect. Daurolirion Comber), eight species and one variety of sect. Leucolirion including L. brownii F. E. Brown ex Miellez (sect. Archelirion Comber), L. henryi Baker and L. rosthornii Diels (sect. Sinomartagon Comber), one species of sect. Archelirion, two species of sect. Liriotypus, three species of sect. Pseudolirium, and six species and one variety of sect. Lophophorum (Table 1).

[Figure omitted. See PDF.]

Table 1. Origin of the material and numerical results of the study.

https://doi.org/10.1371/journal.pone.0087841.t001

To compare pollen features at the genus level, a database with results of the present study and those of Halbritter [34], Kosenko [35], Liang et al. [23], [24], Muratović et al. [25], Wu et al. [28], Zhang et al. [27] and Zhang et al. [36] was established (Table 1). Because the results do not permit a calculation of all the parameters cited above, our observations from 49 samples representing 21 taxa was used only for statistical analyses.

Morphological observations were conducted using SEM. Pollen was directly attached to double-sided adhesive tape and examined under the microscope to locate the pollen. The sample was then taped to the object stage. Following gold spray coating, observation and image acquisition were conducted using a Hitachi S-3400 scanning electron microscope following the method described by Avetissian [37]. All microscopy procedures were performed at the Biotechnology Centre, Beijing Forestry University, China.

Biometric measurements were made using Image-Pro Plus 6.0 (Media Cybernetics, USA). For each sample, 10 (5 for L. brownii from Shennongjia, Hubei Province) fully developed pollen grains were measured. Parameters considered were the polar axis (P), equatorial diameter (E), P/E ratio, lumina diameter and muri width (Table 1). The description of pollen morphology was based on the shape and sculpturing classifications of Baranova [22]. Arithmetic mean and standard deviations are shown for each parameter.

To investigate pollen evolution, pollen features from the database were superimposed onto a phylogenetic framework constructed by internal transcribed spacer (ITS) sequences. We selected 23 taxa, including Cardiocrinum gigantum as an outgroup based on preliminary study [16]. The other 16 taxa, including an outgroup of Northolirion bulbiliferum, were cited [13], [18]–[20]. All of the sequence data from the ITS regions were double-checked visually, edited in BioEdit 5.0.6 [38], and aligned in Clustal X 1.83 [39] using default settings. Phylogenetic analyses were performed using maximum parsimony (MP) and maximum likelihood (ML). The MP tree was constructed with PAUP* 4.0b10 [40]. ML phylogenetic analysis was performed using RAxML v. 7.0.4 with unique model parameters for ITS sequences [41]. A general time reversible model (GTR) was applied with discrete gamma distribution. Bootstrap pseudo replicates were performed 1000 times using the fast bootstrapping option and the best scoring ML tree. Phylogenetic trees were visualized using Treeview [42]. The best scoring tree was visualized with FigTree v1.3.1 (http://tree.bio.ed.ac.uk/).

SPSS 18.0 was used for statistical analyses. One-way analysis of variance (ANOVA) was used to evaluate whether or not differences in pollen characteristics from different provenances were significant. In those cases in which ANOVA revealed significant differences, a least significant difference (LSD) test was performed (Table 2). Correlation analysis was used to evaluate the relationship between pollen parameters and habitat (Our collection of 21 taxa, including 49 samples) (Table 3 and 4).

[Figure omitted. See PDF.]

Table 2. Variance analysis of pollen morphology traits among species from different provenances.

https://doi.org/10.1371/journal.pone.0087841.t002

[Figure omitted. See PDF.]

Table 3. Geographical and climate factors of provenances studied.

https://doi.org/10.1371/journal.pone.0087841.t003

[Figure omitted. See PDF.]

Table 4. Correlation coefficient between the main characteristic of pollen morphology and the geographical and climate factors of provenances.

https://doi.org/10.1371/journal.pone.0087841.t004

Results

Numerical results of the palynological study are summarized in Table 1. According to the Lilium classification by Comber [8], Nishikawa et al. [18], [19] and Du et al. (2013) [16], palynological results for each section are as follows:

Martagon

Representatives of Martagon (Figure 1A–F) have pollen size ranging from L. distichum Nakai et Kamibayashi (91.56×34.57 µm, P/E of 2.66) to L. tsingtauense Gilg (104.70×44.52 µm, P/E of 2.36). The shape can be long-ellipsoidal. L. tsingtauense and L. distichum pollen show similar exine sculptures with a slight difference in that the intersection of the muri presents a protuberant and enlarged columellae in L. tsingtauense pollen (Figure 1E–F). L. tsingtauense and L. distichum correspond to Martagon type pollen, which is in agreement with Baranova [22].

[Figure omitted. See PDF.]

Figure 1. SEM photographs of pollen grains of Lilium.

Figs A, B, L. distichum (Jia 07001). Figs C, D, L. distichum (Du 11005). Figs E, F, L. tsingtauense. Figs G, H, L. brownii (Wang 11006). Figs I, J, L. brownii (Du 09004). Figs K, L, L. brownii (Du 11014). Figs M, N, L. lancifolium (Du 08001). Figs O, P, L. lancifolium (Wang 11001). Scale bars: 10 µm.

https://doi.org/10.1371/journal.pone.0087841.g001

Sinomartagon

Triploid L. lancifolium pollen is relatively larger than the other diploid species in section Sinomartagon. The shape can be ellipsoid to long-ellipsoidal. Species belonging to Sinomartagon (Figure 1M–P; 2A–P; 3A–P; 4A–P) show a decreasing tendency in size with regard to pollen features from subsect. 5a to subsect. 5c as follows: L. lancifolium Thunb (111.20×41.16 µm, P/E of 2.72); L. davidii Duchartre ex Elwes (93.09×33.55 µm, P/E of 2.79); L. davidii var. willmottiae (E. H. Wilson) Raffill (91.89×33.59, P/E of 2.75); L. duchartrei Franch (85.07×33.31, P/E of 2.56); L. leichtlinii Hook. f. var. maximowiczii (Regel) Baker (82.81×33.04, P/E of 2.52); L. pumilum DC (87.39×32.66, P/E of 2.67); L. concolor Salisb. var. pulchellum (Fisch.) Regel (70.77×30.58, P/E of 2.33); L. wardii Stapf ex Stearn (87.77×35.26, P/E of 2.49); and L. taliense Franch (74.70×29.70, P/E of 2.50) (Table 1).

[Figure omitted. See PDF.]

Figure 2. SEM photographs of pollen grains of Lilium.

Figs A, B, L. lancifolium (Du 09001). Figs C, D, L. lancifolium(Du 11001). Figs E, F, L. lancifolium (Du 09002). Figs G, H, L. lancifolium (Du 11002). Figs I, J, L. lancifolium (Du 11003). Figs K, L, L. lancifolium (Du 09003). Figs M, N, L. lancifolium (Wang 11002). Figs O, P, L. lancifolium (Du 11004). Scale bars: 10 µm.

https://doi.org/10.1371/journal.pone.0087841.g002

[Figure omitted. See PDF.]

Figure 3. SEM photographs of pollen grains of Lilium.

Figs A, B, L. dauricum. Figs C, D, L. davidii (Wang 11002). Figs E, F. L. davidii (Du 11008). Figs G, H, L. davidii (Wang 11003). Figs I, J, L. davidii (Du 11009). Figs K, L, L. davidii var. willmottiae (Jia 07002). Figs M, N, L. davidii var. willmottiae (Wang 11004). Figs O, P, L. leichtlinii var. maximowiczii (Du 12001). Scale bars: 10 µm.

https://doi.org/10.1371/journal.pone.0087841.g003

[Figure omitted. See PDF.]

Figure 4. SEM photographs of pollen grains of Lilium.

Figs A, B, L. leichtlinii var. maximowiczii (Du 11010). Figs C, D, L.duchartrei. Figs E, F, L.pumilum. Figs G, H, L.concolor var. pulchellum (Du 11011). Figs I, J, L.concolor var. pulchellum (Du 11012). Figs K, L, L. taliense (Jia 09002). Figs M, N, L.taliense (Du 12002). Figs O, P, L.wardii. Scale bars: 10 µm.

https://doi.org/10.1371/journal.pone.0087841.g004

SEM of pollen from this section showed differences in exine sculpture. L. concolor var. pulchellum clearly shows Concolor type pollen (Figure 4E–H), whereas pollen from other Sinomartagon species show rectangular columellae that correspond to Martagon type.

Leucolirion

Representatives of Leucolirion (Figure 1G–L; 5C–P; 6A–P; 7A–B) have pollen sizes ranging from 87.56×38.69, P/E of 2.28 for L. sulphureum Baker apud Hook. f. to 120.13×53.16, P/E of 2.29 for L. brownii from Shennongjia in the Hubei Province. The exine is ornamented with reticulation, and the shape is ellipsoidal to long-ellipsoidal. Pollen from L. brownii, L. henryi and L. rosthornii correspond to the Martagon type (Figure 1G–L; 5C–J). The pollen exine of L. regale Wilson, L. leucanthum (Baker) Baker, L. sargentiae Wilson and L. sulphureum show rounded columellae, corresponding to the Callose type (Figure 5M–P; 6A–P). Pollen of L. formosanum Wallace does not correspond to any of the pollen types proposed by Baranova (1985), and the exine surface is reticulate with solid muri or irregular-rugulate (Figure 7A–B). We suggest creating a new category of Formosanum type pollen.

[Figure omitted. See PDF.]

Figure 5. SEM photographs of pollen grains of Lilium.

Figs A, B, L.bakerianum. Figs C, D, L.henryi (Du 11015). Figs E, F, L.henryi (Du 12004). Figs G, H, L.rosthornii (Jia 10001). Figs I, J, L.rosthornii (Du 12005). Figs K, L, L. regale (Du 11015). Figs M, N, L. regale (Wang 11007). Figs O, P, L. regale (Du 11016). Scale bars: 10 µm.

https://doi.org/10.1371/journal.pone.0087841.g005

[Figure omitted. See PDF.]

Figure 6. SEM photographs of pollen grains of Lilium.

Figs A, B, L. regale (Du 12006). Figs C, D, L.sargentiae (Jia 11001). Figs E, F, L.sargentiae (Du 11017). Figs G, H, L.sargentiae (Du 11018). Figs I, J, L.sargentiae (Wang 11008). Figs K, L, L. sargentiae (Wang 11009). Figs M, N, L. leucanthum. Figs O, P. L. sulphureum. Scale bars: 10 µm.

https://doi.org/10.1371/journal.pone.0087841.g006

[Figure omitted. See PDF.]

Figure 7. SEM photographs of pollen grains of Lilium.

Figs A, B, L. formosanum. Scale bars: 10 µm.

https://doi.org/10.1371/journal.pone.0087841.g007

Lophophorum

The division of Lophophorum remains controversial. Due to resource constraints, we selected L. bakerianum as the only representative. In a preliminary study, L. bakerianum was placed in the subsect. Lophophorum III [16]. Pollen tetrads were found in L. bakerianum (Figure 5A–B), observed as blunt quadrangles in a polar view with a pollen size of 78.40×79.71 µm, P/E of 1.03. The exine ornamentation was reticulate with muri formed by separated rounded columellae, corresponding to the Callose type. This is the first study to report that pollen grains of L. bakerianum are tetrads.

Discussion

Pollen Morphology and Phylogenetic Implications Based on ITS Sequences

Within the genus Lilium, several diversified pollen exine sculptures were observed. Pollen sculpture patterns appear to reflect phylogenetic relationships and are useful for species or subsection delimitations within a section.

The phylogenetic tree was resolved into six groups based on ITS sequences (Figure 8). Group I is comprised of three clades. Present and previous molecular analyses [18], [19] and cytology studies [21], [43]–[44] indicated that sect. Sinomartagon is complicated and polyphyletic, as supported by the pollen data (Figure 8). Species of sect. Sinomartagon (5a&5b) and L. formosanum and L. brownii (sect. Leucolirion 6b) formed a clade with strong support ([ML] BS = 81, [MP] BS = 77). Sect. Sinomartagon species show two pollen types, Martagon and Concolor. Differences in the shape and arrangement of the types of columellae could have implications on species delimitation. L. concolor var. pulchellum clearly shows Concolor type pollen (Figure 4G–J), and the remainder show Martagon type pollen. Meanwhile, three subtypes were showed within Martagon type pollen (Figure 8A, B, C). Exine sculptures of L. davidii and L. davidii var. willmottiae showed spherical but dorsally and ventrally compressed, and loosely arranged columellae with protuberant and enlarged columellae at the intersection of the muri (Figure 3B–N; 99B). Muri of L. leichtlinii var. maximowiczii have regular long-rectangular columellae (Figure 3O–P; 8C).

[Figure omitted. See PDF.]

Figure 8. Mapping pollen exine sculpture characters on the ML tree of Lilium based on nuclear internal transcribed spacer (ITS) sequence data (adapted from Du et al. unpublished data).

Values along branches represent bootstrap (BS) of ML and MP, respectively. A, B, C, D, Representatives of subsection 5a&5b. (A) L.pumilum (Jia 09001). (B) L.davidii (Du 11008). (C) L. leichtlinii var. maximowiczii (Du 11010). (D) L. concolor var. pulchellum (Du 11012). E, F, Representatives of subsection Leucolirion 6b. (E) L.brownii (Wang 11006). (F) L.formosanum (Du 12008). G, Representatives of subsection 5c (G) L.taliense (Jia 09002). H, I, Representatives of subsection Lophophorum III. (H) L.sempervivoideum (Liang and Zhang 1985). (I) L.bakerianum (Du 12003). J, K, Section Martagon. (J) L.distichum (Jia 07001). (K) L.tsingtauense (Du 11006). L, Representatives of subsection Lophophorum I. (L) L. Lophophorum (Liang and Zhang 1985). M, Representatives of subsection Lophophorum II. (M) L. henricii (Liang and Zhang 1985). N, Representatives of section Pseudolirium. (N) L.canadense (Kosenko, 1999). O, Representatives of section Liriotypus. (O) L. bosniacum (Muratović et al. 2010). P, Representatives of section Sinomartagon I. (P) L.duchartrei (Wang 11005). Q, Representatives of section Archelirion. (Q) L. speciosum var. gloriosoides (Wu et al. 2007). R, Representatives of subsection Leucolirion 6b. (R) L. regale (Du 11015).

https://doi.org/10.1371/journal.pone.0087841.g008

In addition, L. formosanum and L. brownii (Leucolirion 6b) formed a strongly supported clade ([ML] BS = 87, [MP] BS = 86) (Figure 8). In the cytological study, both L. formosanum and L. brownii lacked intercalary satellites on the first two pairs of chromosomes [43], [46] However, the pollen of L. formosanum corresponds to the new Formosanum type, while L. brownii clearly shows Martagon type pollen.

In our analysis, L. taliense and L. wardii ([ML] BS = 80, [MP] BS = 85) formed a clade with L. amoenum, L. sempervivoideum and L. bakerianum ([ML] BS = 79, [MP] BS = 84) that was moderately supported ([ML] BS = 56, [MP] BS<50) (Figure 8). In cytological studies, the L. taliense karyotype resembles that of L. wardii [21], [43]. Morphologically, L. taliense and L. wardii showed a close relationship with each other, with a black line in the central groove of the tepals [2]. Pollen evidence from our study supports this result, as the pollen morphology of the two species is generally similar Martagon type (Figure 4K–P). Morphologically, L. amoenum, L. sempervivoideum and L. bakerianum show nodding or horizontal campanulate flowers, nectaries of inner tepals neither papillose nor with fimbriate projections, tepals without blotch at the base adaxially, and a papillose stem [2]. The L. bakerianum var. delavayi karyotype resembles that of L. sempervivoideum [43]. Pollen tetrads are present in L. amoenum, L. sempervivoideum [23] and L. bakerianum (Figure 5A–B). The similarity in pollen morphology also suggests a close relationship with this clade.

The Martagon clade was monophyletic with strong support ([ML] BS = 98,[MP] BS = 96), which is consistent with previous studies by Gao et al. [13] and Nishikawa et al. [18], [19]. Morphologically, they all show whorled leaves. Although the Martagon clade shows Martagon type pollen, there are minor differences among species in the Martagon clade in the protuberant and enlarged columellae at the intersection of the muri in L. tsingtauense (Figure 1E–F, 8K).

Group II was weakly supported ([ML] BS<50, [MP] BS<50) and consisted of two Lophophorum clades. Lophophorum clade I consisted of L. lophophorum and L. nanum ([ML] BS = 77 and [MP] BS = 71). In palynological research, L. lophophorum and L. nanum are monocolpate or 2–3 porate and show similar pollen morphologies (Figure 8L) [24]. Representatives of L. henricii in subsect. Lophophorum II show Concolor type pollen with much smaller rectangular columellae between two separate irregular round columellae compared with L. concolor var. pulchellum (Figure 8M) [24]. Morphologically, representatives of the Lophophorum clade I show nectaries on the inner tepals with fimbriate/cristate projections on both surfaces, while representatives of subsect. Lophophorum II do not have this feature. The distinct pollen morphology of the species suggests a distinct relationship to support the molecular data [2].

In Group III, the Pseudolirium clade was resolved with strong support ([ML] BS = 77, [MP] BS = 86). Pollen morphologies from representatives of sect. Pseudolirium were similar, with loosely arranged columellae and protuberant and enlarged columellae at the intersection of the muri (Figure 8N) [35].

Group IV, consisting of L. kesselringianum, L. bosniacum and L. carniolicum, was robustly supported ([ML] BS = 97, [MP] BS = 97). The three representatives of European species (sect. Liriotypus) showed Martagon type pollen with densely arranged columellae forming the muri (Figure 8O) [25].

Group V was only comprised of L. duchartrei (Sinomartagon clade I) (Figure 8). Based on morphology, Comber classified it in subsection Sinomartagon 5a. However, in our ITS analyses, L. duchartrei showed a distant relationship with subsect. 5a species as well as other sect. Sinomartagon species (Figure 8), which is in accordance with previous studies [13], [18]–[19]. L. duchartrei (Sinomartagon clade I) karyotypes have been shown to be dissimilar to subsect. 5a Sinomartagon clade I species [21], [43]. L. duchartrei pollen resembles that of L. taliense and L. wardii in exine sculpture (Figure 4C–D; K–P). Given that the phylogenetic position of L. duchartrei is uncertain, classification as a subsection of sect. Sinomartagon may be reasonable.

In Group VI, L. speciosum var. gloriosoides in sect. Archelirion formed a weakly supported group consisting of representatives of Leucolirion 6a ([ML] BS<50 and [MP] BS<50). Leucolirion 6a comprised two subclades ([ML] BS = 74 and [MP] BS = 75). L. speciosum var. gloriosoides showed a distinct relationship with L. henryi and L. rosthornii in the ITS tree, although they all have Martagon-type pollen (Figure 8). In addition, L. henryi and L. rosthornii, which have orange reflexed flowers with prominent papillae and pubescent nectaries, were resolved with strong support ([ML] BS = 98, [MP] BS = 100). Species of L. leucanthum, L. sargentiae, L. sulphureum and L. regale formed a clade with strong support ([ML] BS = 95, [MP] BS = 95). Pollen morphology from our study supports this result, as the pollen morphology of L. henryi resembles L. rosthornii as Martagon type (Figure 5C–H), while L. leucanthum, L. sargentiae, L. sulphureum and L. regale have all generally similar Callose-type pollen (Figure 5K–L; 6A–P).

Within Lilium, we noted that one section may show two or more different pollen types (Figure 8). For instance, sect. Sinomartagon shows two exine sculpture pollen types, Martagon and Concolor, and the former contains three different subtypes. No unique and uniform characters can be used to clearly distinguish sections or some subsections due to shared or overlapping characters.

Trends in Exine Sculpture Evolution and Factors Implicated in Pollen Evolutionary Trends

Changes in pollen morphology occur in response to selective pressure. Pollen grains should be considered a functional unit, with the exine structure as a compromise between the following four main functions [32]: protective, harmomegathic, reservoir and clustering. These functions are in response to physical components of the environment, such as water, nutrients, temperatures and growing season.

The genus Lilium is distributed in East Asia, the Himalayas and Hengduan Mountains, Europe and North America [1], [2]. Patterson and Givnish [4] concluded that Lilium evolved in the Himalayas and then dispersed into Eurasia and North America, although intercontinental dispersal details are not clear [4]. In the karyotype study, the relatively primitive group (sect. Lophophorum (Bur. et Franch.) Wang et Tang) have no intercalary satellites on the first two pairs of chromosomes. As expected, the genus Notholirion which is commonly considered to be closely related to and more ancient than Lilium have also been found to lack intercalary satellites [4], [21], [43], [45], [46], these species are distributed mainly in the Himalayas, which makes this region more likely as the origin of the genus Lilium. Besides, Stewart [45] pointed out that secondary constrictions correlated with chromatin distribution as well as with geographic distribution. In his study, most East Asian and Himalayan species were found to have intercalary satellites on the first two pairs of submedian centromeres, while all North American species and European species lack such feature [21], [25], [43], [45,]. Hence, the situation among the Himalayan and Henduan Mountains species appears more complicated since they seem to have both types, even within the same section indicating that the Hengduan Mountains are more likely as the differentiation centre of the genus [43], [45], [47]. Pollen evidence from our study supports this hypothesis, as the pollen morphology of the species from the Himalayas and Hengduan Mountains contain all of the pollen types in the genus Lilium. It should be emphasized that the Formosanum type of pollen can be found in the genus Nomocharis, which is restricted to the Himalayas and Hengduan Mountains and it has been suggested that it be included in Lilium [13]–[16]. In contrast, the other regions possess a relatively singular pollen type, which led us to investigate Lilium pollen-type evolution.

Based on a pollen study of 30 genera in Liliaceae, including 69 species, using light microscopy (LM), Nair and Sharma [46] proposed the scheme of pollen exine sculpture evolution (Liliaceae) as follows: ornate → retipilate → reticulate → ring-shaped reticulate → scattered reticulate → exineless ornamentation. Regarding pollen morphology of the genus Nomocharis and Lilium under LM and SEM, Liang and Zhang [24] suggested the following evolution trend of pollen exine sculpture in the genus Lilium: ornate → retipilate → reticulate, which is consistent with the hypothesis of Nair and Sharma [46]. In addition, retipilate was divided into three exine sculptures to form three Lilium pollen types, including Martagon, Callose and Concolor. These results led us to propose the following hypothesis of pollen evolution in exine sculptures within Lilium: Martagon type → Callose type → Concolor type → Formosanum type. This also shows a reduction and simplification tendency.

In monocots, reduction and simplification of the exine structure is even more extreme than in dicots and culminate in certain evolutionary orders [48]. Zavada proposed major evolutionary trends of wall structure types in monocots as follows: primitive tectate-columellate (perforate or imperforate) wall structures give rise to monocotyledonous atectate or granular walls and eventual extreme reduction of the exine, which may be completely absent [48]. All morphological characteristics are the consequence of interactions between phylogenetic and environmental constraints [32]. There is no doubt that environmental constraints act on patterns of differentiation in Lilium [49].

In addition to environmental factors, are there other factors such as pollinator affecting pollen differentiation? Numerous recent findings on the intricate flower pollinator networks suggest that differentiation of flora structures and functions in plants are tightly connected with those of pollinators [50]–[53]. Flowers of Lilium are insect-pollinated. Brantjes and Bos’s [54] report that both diurnal and nocturnal hawkmoths visit L. martogon. Skinner [55] observed pollinators of 13 American species of Lilium and recognized five types of pollination: butterfly pollination (L. humboldtii, L. kelleyanum and L. kelloggii), butterfly and hummingbird pollination (L. wigginsii, L. pardalinum and L. pardalinum ssp. vollmeri), hummingbird pollination (L. columbianum, L. occidentale and L. bolanderi), hummingbird and bumblebee pollination (L. parvum and L. maritimum) and hawkmoth pollination (L. washingtonianum and L. parryi). Additional studies reported various flower visitors, including fritillaries to L. concolor var. pulchellum [56], swallowtails to L. dauricum [56], nocturnal hawkmoths to L. formosanum [57] and L. japonicum [58], and hawkmoths and swallowtails to L. auratum [56], [59]. Nevertheless, wild Lilium species are self-compatible. Autogamy might have ensured reproductive success for the species in the environments where pollinators are comparatively rare [60] and the growing season short as well as where other conditions are selective. Pollen of L. martagon and L. dauricum with different pollinators are Martagon type (Table 1) [54], [56]. Thus, for Lilium, the pollinators have little effect on changes of pollen morphology and size. Moreover, in Sonneratiaceae, the genus Duabanga has much smaller pollen grains than Sonneratia while both genera are bat-pollinated, hence absence of correlation with pollinator [32]. However, how floral traits as well as pollen grains of lilies are adapted to various pollinators remains to be studied in more detail.

As mentioned, environmental factors constitute a selective pressure that may produce changes in the pollen structure. The number and arrangement of the types of columellae could influence exine flexibility. Exines with loose arrangement of round columellae (Callose, Concolor) are relatively more flexible compared to exines with tight, rectangular columellae (Martagon). This exine structure corresponds to harmomegathic function and the ability to absorb bending stresses that may occur during desiccation [32].

Relationship between Pollen Size and Conditions in the Genus Lilium

Species at high elevations with relatively extreme conditions in the Himalayas and Hengduan Mountains have a relatively small pollen size (Table 1). Plants growing at high elevations with extreme conditions often have reduced morphological features [61], [62], and pollen size also appears to show this tendency [33], [63]. Furthermore, species native to Northeast China, such as L. leichtlinii var. maximowiczii, L. concolor var. pulchellum and L. dauricum, also have a relatively small pollen size (Table 1). A long winter and short growing season forces rapid completion of plant growth cycles. Hence, small pollen size also appears to be more efficient in tolerating relatively extreme conditions. In contrast, L. brownii and L. formosanum have relatively bigger pollen grains, likely because they occur in an environment suitable for slow development that allow pollen grains to grow to a larger size, as discussed by Muller [32]. Therefore, pollen size could be related to the selective pressure to adapt to environmental conditions. In addition, smaller pollen grains appear to survive better in habitats with extreme conditions.

Within the same species, pollen parameters from different provenances showed some variation (Table 1). Based on one-way ANOVA and LSD tests in 11 species (Table 2), results showed significant or extremely significant differences in pollen parameters in the same species from different habitats within the provenances, except for several single parameters in some species (Table 2). This may be explained as a consequence of interactions between phylogenetic and environmental constraints, as emphasized by Kawano and Kato [31].

Due to the wide variety of habitats in the distribution of Lilium, correlation analysis was performed to test which environmental factor is critical in influencing pollen parameters. The correlation analysis showed a significantly positive correlation between the polar axis and annual precipitation (P<0.05) and an extremely significantly positive correlation between the P/E ratio and annual precipitation (P<0.01) (Table 4). Hence, pollen size and shape showed a significantly positive correlation with annual precipitation. The absolute size of pollen grains influences harmomegathic function [64], which can be reduced to such a degree that exine elasticity alone can accommodate the necessary changes in volume to adapt to environmental conditions [30].

Pollen exine sculptures of different species in similar environmental conditions may be different types. For instance, L. concolor var. pulchellum possesses Concolor type pollen, whereas L. dauricum shows Martagon type pollen, and both occur in Northeast China under similar environmental conditions. Species native to the Himalayas and Hengduan Mountains show all of the pollen types within Lilium. The evolutionary trend of exine sculptures is likely not definitively correlated with pollen size.

Conclusion

In conclusion, we suggest recognising a new pollen type, Formosanum, to accommodate pollen from L. formosanum. Pollen sculpture patterns appear to reflect phylogenetic relationships and are useful for species or subsection delimitation within sections. The reduction and simplification evolutionary trend of pollen sculpture and size could be related to the selective pressure to better adapt to environment conditions, especially extreme environmental conditions. In addition, the evolutionary trends of exine sculpture and pollen size and shape are not definitively correlated. Pollen size and shape show a significantly positive correlation with annual precipitation. However, additional studies are required to confirm these evolution hypotheses in a broad sample of the genus using techniques such as transmission electron microscopy (TEM) to observe strata of the exine.

Author Contributions

Conceived and designed the experiments: GXJ YPD. Performed the experiments: YPD ZXW. Analyzed the data: YPD SL HBH. Contributed reagents/materials/analysis tools: YPD ZXW CW. Wrote the paper: YPD.

Citation: Du Y-p, Wei C, Wang Z-x, Li S, He H-b, Jia G-x (2014) Lilium spp. pollen in China (Liliaceae): Taxonomic and Phylogenetic Implications and Pollen Evolution Related to Environmental Conditions. PLoS ONE 9(1): e87841. https://doi.org/10.1371/journal.pone.0087841

References

1. Mcrae EA (1998) Lily species. In: Lilies. Portland: Timber Press. 105–204.

2. Liang SY, Tamura M (2000) In: Wu ZY, Raven PH (eds.) Flora of China. Beijing: Science Press/Missouri Botanical Garden Press. Volume 24: 135–159.

3. De Jong PC (1974) Some notes on the evolution of lilies. Lily Year Book of the North American Lily Society. 27: 23–28.

4. Patterson TB, Givnish TJ (2002) Phylogeny, concerted convergence, and phylogenetic niche conservatism in the core Liliales: Insights from rbcL and ndhF sequence data. Evolution 56: 233–252.

5. Endlicher M (1836–1840) Genera plantarum secundum ordines naturales disposita. Vienna: Beck.

6. Baker JG (1871) A new synopsis of all the known lilies. The gardeners’ chronicle. 28,104.

7. Wilson E H (1925) The lilies of Eastern Asia. Dulau, London.

8. Comber HF (1949) A new classification of the genus Lilium. Lily Year Book RHS. 13, 86–105.

9. Wang FZ, Tang J (1980) Lilium L. In: Flora Reipublicae Popularis Sinicae. Beijing: Science Press. Volume 14: 116–157.

10. Liang SY (1980) Flora Reipublicae Popularis Sinicae. Anagiospermae, Monocotyledoneae Liliaceae (I). Beijing: Science Press. 14: 116–157.

11. Haw SG (1986) The lilies of China. Portland: Timber Press.

12. Dubouzet JG, Shinoda K (1999) Phylogenetic analysis of the internal transcribed spacer region of Japanese Lilium species. Theoretical and Applied Genetics 98: 954–960.

13. Gao YD, Hohenegger Markus, Harris AJ, Zhou SD, He XJ, et al. (2012) A new species in the genus Nomocharis Franchet (Liliaceae): evidence that brings the genus Nomocharis into Lilium. Plant Systematics and Evolution 298: 69–85.

14. Hayashi K, Kawano S (2000) Molecular systematics of Lilium and allied genera (Liliaceae): phylogenetic relationships among Lilium and related genera based on the rbcL and matK gene sequence data. Plant Species Biology 15: 73–93.

15. Peruzzi L, Leitch IJ, Caparelli KF (2009) Chromosome diversity and evolution in Liliaceae. Annals of Botany103: 459–475.

16. Du YP, He HB, Wang ZX, Li S, Wei C et al.. (2013) Molecular phylogeny and genetic variation in the genus Lilium native to China based on nrDNA ITS sequences. DOI: 10.1007/s10265-013-0600-4.

17. İkinci N, Oberprieler C, Güner A (2006) On the origin of European lilies: phylogenetic analysis of Lilium section Liriotypus (Liliaceae) using sequences of the nuclear ribosomal transcribed spacers. Willdenowia 36: 647–656.

18. Nishikawa T, Okazaki K, Uchino T, Arakawa K, Nagamine T (1999) A molecular phylogeny of Lilium in the internal transcribed spacer region of nuclear ribosomal DNA. Journal of Molecular Evolution 49: 238–249.

19. Nishikawa T, Okazaki K, Uchino T, Arakawa K, Nagamine T (2001) Phylogenetic analysis of section Sinomartagon in genus Lilium using sequences of the internal transcribed spacer region in nuclear ribosomal DNA. Breeding Science 51: 39–46.

20. Rešetnik I, Liber Z, Satovic Z, Cigić P, Nikolić T (2007) Molecular phylogeny and systematics of the Lilium carniolicum group (Liliaceae) based on nuclear ITS sequences. Plant Systematics and Evolution 265: 45–58.

21. Smyth DR, Kongsuwan K, Wisudharomn S (1989) A survey of C-band patterns in chromosomes of Lilium (Liliaceae). Plant Systematics and Evolution 163: 53–69.

22. Baranova M V (1985) Palynoderm ultrastructure and morphological types of pollen grains in the genus Lilium (Liliaceae). – Botanicheskiy Zhurnal (Leningrad) 70: 297–304.

23. Liang SY, Zhang WX (1984) Pollen tetrads in Lilium with a discussion on the delimitation between L. semperviroideum and its Ally. Acta Phytotaxonomica Sinica 22: 297–300.

24. Liang SY, Zhang WX (1985) Pollen morphology of the genus Nomocharis and its delimitation with Lilium. Acta Phytotaxonomica Sinica 23: 405–417.

25. Muratović E, Bogunić F, Šoljan D, Martin J, Valles J, et al. (2010) Stomata and pollen grain characteristics of two endemic lilies: Lilium bosniacum and L. carniolicum (Liliaceae). Phytologia Balcanica 16 (2): 285–292.

26. Li JQ, Qin XP (1993) Studies on pollen morphological of some Lilies. Journal of Wuhan Botanical Research 11(2): 120–126.

27. Zhang LN, Lei JJ, Han L, Hu XY (2006) Study on Pollen Morphology and Ultrastructure of Lilium spp. Acta Agriculturae Boreali-occidentalis Sinica 15(6): 174–178.

28. Wu ZH, Shi JS, Xi ML, Liu GX (2007) Pollen characteristics and relationships of 12 species and 6 cultivars of Lilium. Journal of Zhejiang Forestry College 24(4): 406–412.

29. Wang X, Lei JJ, Guan JZ (2007) Study on pollen diversity within Lilium concolor Salisb. Heilongjiang Agricultural Science 1: 47–49.

30. Liu DY (2011) Study on the genetic diversity and utilization of the germplasm resource of Lilium pumilum DC. Phylogeny. Ph.D. Thesis, Beijing Forestry University.

31. Kawano S, Kato H (1995) Evolutionary biology of Trillium and related genera (Trilliaceae) II. Cladistic analyses on gross morphological characters and phylogeny of the genus Trillium. Plant Species Biology 10: 169–183.

32. Muller J (1979) Form and function in angiosperm pollen. Annals of the Missouri Botanical Garden 66: 593–632.

33. Hidalgo O, Susanna A, Garcia-Jacas N, Martin J (2008) From acaveate to caveate: evolution of pollen types in the Rhaponticum group (Asteraceae, Centaureinae) related to extreme conditions. Botanical Journal of the Linnean Society 158: 499–510.

34. Halbritter H, Hesse M (1993) Sulcus morphology in some monocot families. Grana 32(2): 87–99.

35. Kosenko V N (1999) Contributions to the pollen morphology and taxonomy of the Liliaceae. Grana 38: 20–30.

36. Zhang YL, Zhang QX, Xie SL (2010) Pollen Morphology of 8 Species in Lilium from Qinba Mountain Areas. Acta Agriculturae Boreal–occidentalis Sinica 19(1): 144–146.

37. Avetissian EM (1950) Simplified method of pollen preparation by acetolysis. – Zhurn. Russk. Bot. Obshch. Akad. Nauk 35(385): 6–7.

38. Hall TA (1999) BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acid Symposium Series 41: 95–98.

39. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG (1997) CLUSTAL-X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Research 25: 4876–4882.

40. Swofford DL (2000) PAUP* 4.0: Phylogenetic analysis using parsimony (* and other methods). Version 4.0 b10. Sunderland, MA: Sinauer Associates.

41. Stamatakis A (2006) RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22: 2688–2690.

42. Page RMD (1996) Treeview: An application to display phylogenetic trees on personal computers. Computer Applications in the Biosciences 12: 357–358.

43. Gao YD, Zhou SD, He XJ (2011) Karyotype studies in thirty-two species of Lilium (Liliaceae) from China. Nordic Journal of Botany 29: 746–761.

44. Gao YD, Zhou SD, Jin X, Wan HJ (2012) Chromosome diversity and evolution in tribe Lilieae (Liliaceae) with emphasis on Chinese species. Journal of Plant Research 125: 55–69.

45. Stewart RN (1947) The morphology of somatic chromosomes in Lilium. American Journal of Botany 34: 9–26.

46. Nair P K, Sharma S (1965) Pollen morphology of Liliaceae. Journal of Palynology 1: 38–61.

47. Zhou SD(2008)The phylogenetic classification and evolution of trib Lilieae (Liliaceae s.str.) in China. – PhD thesis, Sichuan University.

48. Zavada MS (1983) Comparative morphology of monocot pollen and evolutionary trends of apertures and wall structures. The Botanical Review 49: 333–378.

49. Hayashi K, Kawano S (2000) Molecular systematics of Lilium and allied genera (Liliaceae): phylogenetic relationships among Lilium and related genera based on the rbcL and matK gene sequence data. Plant Species Biology 15: 73–93.

50. Thien LB., Kawano S., Azuma H., Latimer S., Devall MS. et al.. (1998) The floral biology of the Magnoliaceae. In: Magnolias and Their Allies (ed. D. Hunt) 37–58. International Dendrology Society and the Magnolia Society, Milborne Port, UK.

51. Gottsberger G (1999) Pollination and evolution in neotropical Annonaceae. Plant Species Biology 14: 143–152.

52. Knudsen JT (1999) Floral scent differentiation among coflowering, sympatric species of Geonoma(Arecaceae). Plant Species Biology14: 137–142.

53. Williams NH, Whitten WM (1999) Molecular phylogeny and floral fragrances of male euglossine bee-pollinated orchids: a study of Stanhopea (Orchidaceae). Plant Species Biology 14: 129–136.

54. Brantjes NBM, Bos JJ (1980) Hawkmoth behaviour and flower adaptation reducing self pollination in two Liliflorae. New Phytolorogist 84: 139–143.

55. Skinner MW (1988) Comparative pollination ecology and floral evolution in Pacific coast Lilium (PhD dissertation). Harvard University, Cambridge.

56. Hayashi K, Kawano S (2007)Lilium dauricumKer-Gawler. In: Kawano S. (ed.).Life-History Monographs of Japanese Plants III. Hokkaido University Press, Sapporo, 1–23. (In Japanese with an English summary.).

57. Rodger JG, van Kleunen M, Johnson SD (2010) Does specialized pollination impede plant invasions? International Journal of Plant Sciences171: 382–391.

58. Chiba S, Shimizu T (2004) Life history of the near threatened Lilium japonicumThunb. in Nagano and pollinators. Bulletin of the Botanical Society of Nagano37: 1–8 (In Japanese)..

59. Morinaga SI, Kumano Y, Ota A, Yamaoka R, Sakai S (2009) Day–night fluctuations in floral scent and their effects on reproductive success in Lilium auratum. Population Ecology 51: 187–195.

60. Arroyo MTK, Armesto J, Primack R (1985) Community studies in pollination ecology in the high temperate Andes of central Chile. II. Effect of temperature on visitation rates and pollination possibilities. Plant Systematics and Evolution 149: 187–203.

61. Hedberg O (1964) Features of Afroalpine plant ecology. Acta Phytogeographica Suecica 49: 5–144.

62. Cabrera AL (1968) Ecología vegetal de la puna. In: Troll C, ed. Geoecology of the mountainous regions of the tropical Americas. Bonn: Dummer in Kommission. 91–116.

63. Tate JA, Simpson BB (2004) Breeding system evolution in Tarasa (Malvaceae) and selection for reduced pollen grain size in the polyploid species. American Journal of Botany 91: 207–213.

64. Wodehouse RP (1935) Pollen Grains. New York: Mc Graw-Hill Book Company Incorporated. 323–340.