Research Article
Print
Research Article
Phylogenomic and morphological evidence reveal a new species of spider lily, Lycoris longifolia (Amaryllidaceae) from China
expand article infoYi-Lei Lou§, Dai-Kun Ma|, Ze-Tao Jin|, Hui Wang, Lu-Huan Lou, Shui-Hu Jin, Kun Liu#, Bin-Bin Liu|
‡ Zhejiang Agriculture and Forestry University, Hangzhou, China
§ Beijing Forestry University, Beijing, China
| Institute of Botany, Chinese Academy of Sciences, Beijing, China
¶ Nanjing Agricultural University, Nanjing, China
# Anhui Normal University, Wuhu, China
Open Access

Abstract

Lycoris longifolia, a new species from China, was described and illustrated here. Our phylogenomic evidence based on whole plastomes strongly supported the separate phylogenetic position of this new species, and morphologically it could also be distinguished by its long leaves with a distinct purplish-red midrib on the abaxial surface.

Keywords

Lycoris, morphological, phylogenomics, whole plastome

Introduction

The genus Lycoris Herb., including ca. 13–20 species of flowering plants in the family Amaryllidaceae, subfamily Amaryllidoideae, is native to eastern and southern Asia. Herbert described the first species in 1820, L. aurea (L’Hér.) Herb., which has important ornamental and medicinal values (Hsu et al. 1994). In the mid-20th century, an American horticulturist, Hayward, did much work on introducing and cultivating Lycoris species. Given the easily distinguished habit of the populations of Lycoris aurea distributed in northern Taiwan and southernmost Japan, i.e., the leaves appear in autumn, about a month later than L. aurea, and no remains of leaf bases (Hsu et al. 1994), Hayward described these populations as a new species, L. traubii W.Hayw. (Hayward 1957; Hsu et al. 1994; Kurita 1987). Having narrower perianth lobes and long-exserted stamens (Hsu et al. 1994; Ji and Meerow 2000), the populations from South Gansu (Kang Xian) and Northwest Hubei (Feng Xian) were described as a variety of Lycoris aurea, as L. aurea var. angustitepala P.S.Hsu, Kurita, Z.Z.Yu & J.Z.Lin (Hsu et al. 1994). In the last decades, numerous new species or hybrids of Lycoris have been published in its diversity center, i.e. mainland China, such as L. hunanensis M.H.Quan, L.J.Ou & C.W.She (Quan et al. 2013), L. × hubeiensis KunLiu (Meng et al. 2018), L. tsinlingensis P.C.Zhang, YiJunLu & TingWang (Lu et al. 2020), and L. wulingensis S.Y.Zhang (Zhang et al. 2021). Nowadays, more than 30 species and varieties have been recognized in the genus (Hsu et al. 1994; Ji and Meerow 2000; Kim 2004; Quan et al. 2013; Meng et al. 2018; Lu et al. 2020; Zhang et al. 2021), and nearly 20 of them are from China.

During our recent field explorations in Sichuan Province, China, we collected a wild flowering plant of Lycoris, which resembles L. aurea with yellow flowers. However, it could be easily distinguished from L. aurea by markedly long leaves with a distinct purplish-red midrib on the abaxial surface. Our morphological and molecular evidence strongly supported this population as a new Lycoris species.

Materials and methods

Total genomic DNAs were extracted from 15mg of silica gel dried leaves using a modified CTAB method (Li et al. 2013). The library was prepared at the Molecular Biology Experiment Center, Germplasm Bank of Wild Species in Southwest China using a NEBNext UltraTM II DNA Library Prep Kit (New England Biolabs, Ipswich, MA, USA). The paired-end (150 bp) reads have been generated on the HiSeq 2500 (Illumina, Inc., San Diego, CA, USA) platform in Beijing Genomics Institution (BGI) (Shenzhen, China), ca. 8 GB of raw data for this new species. The raw reads have been deposited in the NCBI Sequence Read Archive in the BioProject (PRJNA857321) with the Run number SRR20072320.

The raw data generated from the Illumina platform was trimmed by Trimmomatic v.0.40 (Bolger et al. 2014) with the default parameters. The clean data was checked by FastQC (Andrews 2010) for quality control. We used the successive assembly approach (Zhang et al. 2015), combining the reference-based and the de novo assembly methods to assemble the chloroplast genome; this method has been performed well in various angiosperm lineages (e.g., Liu et al. 2019, 2020a, 2020b, 2021, 2022; Wang et al. 2020). We annotated the assembled chloroplast genome with two reference genomes (MK353216 and MH118290) downloaded from GenBank, and checked the start and stop codons carefully by translating the coding sequences of plastome into proteins in Geneious Prime (Kearse et al. 2012). We also verified the boundary of two reverse complementary repeats in the plastome using Find Repeats embedded in Geneious Prime (Kearse et al. 2012). The assembled chloroplast genome has been submitted to GenBank with the accession number ON960856. The gene map of the new species Lycoris longifolia chloroplast genome was drawn by OrganellarGenomeDRAW (OGDRAW) version 1.3.1 (Greiner et al. 2019).

We downloaded 24 chloroplast genomes from GenBank as the ingroup and Narcissus poeticus L. as the outgroup for phylogenomic analysis. Given the potential effect of the missing data for the accurate phylogenetic inference, we used the whole plastome (WP) and 78 coding sequences (CDS) to estimate the phylogeny, respectively. Because of the nearly identical sequence of two inverted repeats (IR) in plastomes, we only included one repeat of IR region for downstream WP analyses. Each CDS sequence was extracted separately by Geneious Prime; the WP matrix was aligned with MAFFT v. 7.480 (Nakamura et al. 2018) with default parameters. The WP matrix was trimmed using trimAL v1.2 (Capella-Gutiérrez et al. 2009) with a heuristic method to decide on the best-automated method. All 78 CDS sequences of each plastome were concatenated by AMAS (Borowiec 2016). The best-fit partitioning schemes and/or nucleotide substitution models for the 78 CDS sequences were estimated using PartitionFinder2 (Stamatakis 2006; Lanfear et al. 2016), under the corrected Akaike information criterion (AICc) and linked branch lengths, as well as with rcluster (Lanfear et al. 2014) algorithm options. The resulting optimal partitioning schemes and evolutionary model for each CDS sequence were applied for the following tree inference. We used IQ-TREE2 v. 2.1.3 (Minh et al. 2020) with 1000 SH-aLRT and the ultrafast bootstrap replicates and RAxML 8.2.12 (Stamatakis 2014) with GTRGAMMA model for each partition and clade support assessed with 200 rapid BS replicates for the Maximum Likelihood (ML) analysis. The BI was performed with MrBayes 3.2.7 (Ronquist et al. 2012). The Markov Chain Monte Carlo (MCMC) analyses were run for 10,000,000 generations. The stationarity was regarded to be reached when the average standard deviation of split frequencies remained below 0.01. Trees were sampled every 1,000 generations, and the first 25% of samples were discarded as burn-in. The remaining trees were used to build a 50% majority-rule consensus tree. Considering the possible different evolutionary forces in the chloroplast genome, we also used ASTRAL-III (Zhang et al. 2018) for estimating a coalescent-based species tree based on the 78 CDS sequences.

Results

The chloroplast genome of Lycoris longifolia was 158,413 bp in length, with a typical quadripartite structure consisting of a large single copy region and a small single copy region separated by two long inverted repeats (Fig. 1B). And this structure has been nearly similar to other Lycoris chloroplast genomes released in GenBank. They contained the same number of coding sequences (78), tRNAs (30), and rRNAs (4).

Figure 1. 

Maximum likelihood phylogeny of Lycoris inferred from RAxML analysis of the whole plastome data A numbers above the branches indicate the SH-aLRT support and Ultrafast Bootstrap support (black) by IQ-TREE2, the bootstrap support (red) by RAxML, the posterior probabilities (green) by MrBayes, and the local posterior possibility (orange) by ASTRAL-III. The upper-left inset was a gene map of the new species Lycoris longifolia chloroplast genome B.

The WP matrix was 131,649 bp in length, with the poor sites trimmed by trimAL (Capella-Gutiérrez et al. 2009); the concatenated CDSs were 67,953 bp in length. These two matrices generated seven trees (Fig. 1A, Suppl. material 16). The four ML trees (Fig. 1A, Suppl. material 1, 3, 4), two Bayesian trees (Suppl. material 2, 5), and the species tree (Suppl. material 6) resulted in a consistent phylogenetic position, and this new species, Lycoris longifolia, formed a separate clade (Fig. 1A). This result showed that this new species has been distant from other species in Lycoris. The examined morphological characters, long leaves and purplish-red midrib abaxially, also supported its distinguished status.

Taxonomy

Lycoris longifolia L.H.Lou, sp. nov.

Chinese name: 长叶石蒜 Figs 2, 3

Diagnosis

Most similar to L. aurea but differs from it by markedly longer leaves, abaxially with a distinct purplish-red midrib.

Type

China. Sichuan: Ya’an, Yucheng, Bifengxia, Houyancun, Yanjiashan, under the shrub along the stream, elevation ca. 950 m, 10 May 2021, L.H. Lou & Y.L. Lou 8765 (holotype PE [barcode 02347459]!; isotypes KUN!, PE [barcode 02347457]!).

Additional Specimens examined

China. Sichuan: Ya’an, Yucheng, Bifengxia, Houyancun, Yanjiashan, under the shrub along the stream, elevation ca. 950 m, 30 July 2021, L.H. Lou & Y.L. Lou 8766 (paratype PE [barcode 02347458]!).

Description

Bulbous perennial. Bulbs subglobose, 3–6 cm diam., tunics membranous, dark brown. Leaves ligulate, acute at the apex, ca. 80–120 × 1.5–2 cm, absent at the flowering time and appearing in autumn, dark green, with a prominent midrib on the abaxial surface, abaxial midrib distinctly purplish-red. Inflorescence scapose, umbellate; scape solid, 70–75 cm long, ca. 2.0 cm diam. at base, light green with purplish-red base; involucral bracts 2, lanceolate, 5.0–9.0 cm long by 1.8 cm wide at base, membranous, light green; bracteoles membranous, lanceolate, 1.0–4.0 cm long. Flowers 5–7 per umbel; pedicels 2–2.5 cm long; perianth with 6 tepals; tube ca. 1.5 cm; lobes yellow, abaxially with white mid-vein, strongly recurved, narrowly oblanceolate, ca. 7 × 0.8–1.0 cm, margin strongly undulate. Stamen filaments 6, creamy-yellow, slightly longer than perianth; anther light purplish, dorsifixed, 8–10 mm long before anthesis. Style creamy-yellow but rose-red at apex, slightly exceeding filaments; stigma purplish-red; ovary green, ovoid, ca. 5 mm long.

Phenology

Scape produced from July to August, and vegetative growth from September to May next year. This new species grows along the forest edge near the riverside, and Quercus glauca Thunb. and Pinus massoniana Lamb. are the dominant associated species.

Etymology

The specific epithet alludes to length of leaf blades, a diagnostic character.

Distribution

This new species has been narrowly discovered in Ya’an, Sichuan, China. Some localities of Southwestern China have been poorly discovered, and a comprehensive floristic investigation will help elucidate the germplasm resources.

Figure 2. 

Field photos of Lycoris longifolia A-C flowers D bulb E the distinct purplish-red midrib abaxially F vegetative growth period, showing the long leaves.

Figure 3. 

Illustration of Lycoris longifolia, drawn by Ai-Li Li (PE).

Key to the species of Lycoris in China

1 Flowers actinomorphic 2
Flowers zygomorphic 6
2 Margin of perianth lobes not undulate 3
Margin of perianth lobes basally minutely undulate 4
3 Perianth pale purple but apically blue L. sprengeri
Perianth white or yellow L. longituba
4 Perianth purple L. squamigera
Perianth not purple 5
5 Perianth yellow L. anhuiensis
Perianth white, abaxially with purple midvein L. incarnata
6 Leaves appearing in autumn 7
Leaves appearing in spring 15
7 Perianth yellow or ocher-yellow 8
Perianth bright red, deep red, rose-red, or white 11
8 Perianth yellow; leaves 1.5–5 cm wide 9
Perianth ocher-yellow; leaves 1.0–1.5 cm wide 10
9 Leaves ensiform, ca. 60 × 2–5 cm L. aurea
Leaves ligulate, ca. 100 × 1.5–2 cm L. longifolia
10 Leaves ensiform, apex acuminate L. straminea
Leaves ligulate, apex obtuse L. hunanensis
11 Perianth bright red, deep red, or rose-red 12
Perianth white L. houdyshelii
12 Perianth bright red or deep red, lobes strongly recurved 13
Perianth rose-red, lobes slightly recurved 14
13 Perianth bright red L. radiata
Perianth deep red with white but faintly pale red filaments L. hubeiensis
14 Leaves ligulate, ca. 0.8 cm wide L. rosea
Leaves narrowly ligulate, ca. 0.5 cm wide L. wulingensis
15 Perianth white 16
Perianth yellow or orange-red 17
16 Perianth white without pink stripes L. caldwellii
Perianth white with pink stripes L. shaanxiensis
17 Perianth yellow in bud, becoming orange-red as buds develop L. tsinlingensis
Perianth yellow 18
18 Perianth lobes without red stripes L. chinensis
Perianth lobes abaxially with red stripes L. guangxiensis

Acknowledgements

All the phylogenomic analyses have been run on the Dell T7920 workstation (Institute of Botany, Chinese Academy of Sciences); this study is also supported by the Bioinformatics Center of Nanjing Agricultural University. National Natural Science Foundation of China (Grant No. 32000163, 32270216, and 31620103902) and Project of National Plant Specimen Resource Center (NPSRC) (E0117G1001) support this study.

References

  • Capella-Gutiérrez S, Silla-Martínez JM, Gabaldón T (2009) trimAl: A tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 25(15): 1972–1973. https://doi.org/10.1093/bioinformatics/btp348
  • Greiner S, Lehwark P, Bock R (2019) OrganellarGenomeDRAW (OGDRAW) version 1.3.1: Expanded toolkit for the graphical visualization of organellar genomes. Nucleic Acids Research 47(W1): W59–W64. https://doi.org/10.1093/nar/gkz238
  • Hayward W (1957) Lycoris traubii, sp. nov. Plant Life 13: 40–41.
  • Hsu BS, Kurita S, Yu ZZ, Lin JZ (1994) Synopsis of the genus Lycoris (Amaryllidaceae). SIDA, Contributions to Botany 16: 301–331. https://doi.org/10.2307/41967120
  • Ji ZH, Meerow AW (2000) Amaryllidaceae. In: Wu ZY, Raven PH, Hong DY (Eds) Flora of China. Science Press & St. Louis: Missouri Botanical Garden Press, Beijing, 264–273.
  • Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, Buxton S, Cooper A, Markowitz S, Duran C, Thierer T, Ashton B, Meintjes P, Drummond A (2012) Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28(12): 1647–1649. https://doi.org/10.1093/bioinformatics/bts199
  • Kurita S (1987) Variation and evolution in the karyotype of Lycoris, Amaryllidaceae III. Intraspecific variation in the karyotype of L. traubii Hayward. Cytologia 52(1): 117–128. https://doi.org/10.1508/cytologia.52.117
  • Lanfear R, Calcott B, Kainer D, Mayer C, Stamatakis A (2014) Selecting optimal partitioning schemes for phylogenomic datasets. BMC Evolutionary Biology 14(1): 82. https://doi.org/10.1186/1471-2148-14-82
  • Lanfear R, Frandsen PB, Wright AM, Senfeld T, Calcott B (2016) PartitionFinder 2: New methods for selecting partitioned models of evolution for molecular and morphological phylogenetic analyses. Molecular Biology and Evolution 34: 772–773. https://doi.org/10.1093/molbev/msw260
  • Liu BB, Hong DY, Zhou SL, Xu C, Dong WP, Johnson G, Wen J (2019) Phylogenomic analyses of the Photinia complex support the recognition of a new genus Phippsiomeles and the resurrection of a redefined Stranvaesia in Maleae (Rosaceae). Journal of Systematics and Evolution 57(6): 678–694. https://doi.org/10.1111/jse.12542
  • Liu BB, Campbell CS, Hong DY, Wen J (2020a) Phylogenetic relationships and chloroplast capture in the Amelanchier-Malacomeles-Peraphyllum clade (Maleae, Rosaceae): Evidence from chloroplast genome and nuclear ribosomal DNA data using genome skimming. Molecular Phylogenetics and Evolution 147: e106784. https://doi.org/10.1016/j.ympev.2020.106784
  • Liu BB, Liu GN, Hong DY, Wen J (2020b) Eriobotrya belongs to Rhaphiolepis (Maleae, Rosaceae): Evidence from chloroplast genome and nuclear ribosomal DNA data. Frontiers in Plant Science 10: e1731. https://doi.org/10.3389/fpls.2019.01731
  • Liu BB, Ma ZY, Ren C, Hodel RGJ, Sun M, Liu XQ, Liu GN, Hong DY, Zimmer EA, Wen J (2021) Capturing single‐copy nuclear genes, organellar genomes, and nuclear ribosomal DNA from deep genome skimming data for plant phylogenetics: A case study in Vitaceae. Journal of Systematics and Evolution 59(5): 1124–1138. https://doi.org/10.1111/jse.12806
  • Liu BB, Ren C, Kwak M, Hodel RGJ, Xu C, He J, Zhou WB, Huang CH, Ma H, Qian GZ, Hong DY, Wen J (2022) Phylogenomic conflict analyses in the apple genus Malus s.l. reveal widespread hybridization and allopolyploidy driving diversification, with insights into the complex biogeographic history in the Northern Hemisphere. Journal of Integrative Plant Biology 64(5): 1020–1043. https://doi.org/10.1111/jipb.13246
  • Lu YJ, Wang T, Wang YC, Zhang PC (2020) Lycoris tsinlingensis (Amaryllidaceae), a new species from Shaanxi, China. Annales Botanici Fennici 57(4–6): 193–196. https://doi.org/10.5735/085.057.0424
  • Meng WQ, Zheng L, Shao JW, Zhou SB, Liu K (2018) A new natural allotriploid, Lycoris × hubeiensis hybr. nov. (Amaryllidaceae), identified by morphological, karyological and molecular data. Nordic Journal of Botany 36(6): njb-01780. https://doi.org/10.1111/njb.01780
  • Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, Von Haeseler A, Lanfear R (2020) IQ-TREE 2: New models and efficient methods for phylogenetic inference in the genomic era. Molecular Biology and Evolution 37(5): 1530–1534. https://doi.org/10.1093/molbev/msaa015
  • Quan MH, Ou LJ, She CW (2013) A new species of Lycoris (Amaryllidaceae) from Hunan, China. Novon: A Journal for Botanical Nomenclature 22(3): 307–310. https://doi.org/10.3417/2011013
  • Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, Larget B, Liu L, Suchard MA, Huelsenbeck JP (2012) MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Systematic Biology 61(3): 539–542. https://doi.org/10.1093/sysbio/sys029
  • Wang YB, Liu BB, Nie ZL, Chen HF, Chen FJ, Figlar RB, Wen J (2020) Major clades and a revised classification of Magnolia and Magnoliaceae based on whole plastid genome sequences via genome skimming. Journal of Systematics and Evolution 58(5): 673–695. https://doi.org/10.1111/jse.12588
  • Zhang N, Wen J, Zimmer EA (2015) Congruent deep relationships in the grape family (Vitaceae) based on sequences of chloroplast genomes and mitochondrial genes via genome skimming. PLoS ONE 10(12): e0144701. https://doi.org/10.1371/journal.pone.0144701
  • Zhang C, Rabiee M, Sayyari E, Mirarab S (2018) ASTRAL-III: Polynomial time species tree reconstruction from partially resolved gene trees. BMC Bioinformatics 19(S6): 153. https://doi.org/10.1186/s12859-018-2129-y

Supplementary materials

Supplementary material 1 

Figure S1

Yi-Lei Lou, Dai-Kun Ma, Ze-Tao Jin, Hui Wang, Lu-Huan Lou, Shui-Hu Jin, Kun Liu, Bin-Bin Liu

Data type: Images.

Explanation note: Maximum likelihood phylogeny of Lycoris inferred from IQ-TREE2 analysis of the whole plastome data. Numbers above the branches indicate the SH-aLRT support and Ultrafast Bootstrap support.

This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.
Download file (2.14 MB)
Supplementary material 2 

Figure S2

Yi-Lei Lou, Dai-Kun Ma, Ze-Tao Jin, Hui Wang, Lu-Huan Lou, Shui-Hu Jin, Kun Liu, Bin-Bin Liu

Data type: Images.

Explanation note: Bayesian inference phylogeny of Lycoris inferred from MrBayes analysis of the complete chloroplast genome data. Numbers above the branches indicate the posterior probabilities (PP).

This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.
Download file (2.14 MB)
Supplementary material 3 

Figure S3

Yi-Lei Lou, Dai-Kun Ma, Ze-Tao Jin, Hui Wang, Lu-Huan Lou, Shui-Hu Jin, Kun Liu, Bin-Bin Liu

Data type: Images.

Explanation note: Maximum likelihood phylogeny of Lycoris inferred from IQ-TREE2 analysis of the concatenated 78 plastid coding genes. Numbers above the branches indicate the SH-aLRT support and Ultrafast Bootstrap support.

This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.
Download file (2.15 MB)
Supplementary material 4 

Figure S4

Yi-Lei Lou, Dai-Kun Ma, Ze-Tao Jin, Hui Wang, Lu-Huan Lou, Shui-Hu Jin, Kun Liu, Bin-Bin Liu

Data type: Images.

Explanation note: Maximum likelihood phylogeny of Lycoris inferred from RAxML analysis of the concatenated 78 plastid coding genes. Numbers above the branches indicate the bootstrap support.

This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.
Download file (2.14 MB)
Supplementary material 5 

Figure S5

Yi-Lei Lou, Dai-Kun Ma, Ze-Tao Jin, Hui Wang, Lu-Huan Lou, Shui-Hu Jin, Kun Liu, Bin-Bin Liu

Data type: Images.

Explanation note: Bayesian inference phylogeny of Lycoris inferred from MrBayes analysis of the complete chloroplast genome data. Numbers above the branches indicate the posterior probabilities (PP).

This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.
Download file (2.14 MB)
Supplementary material 6 

Figure S6

Yi-Lei Lou, Dai-Kun Ma, Ze-Tao Jin, Hui Wang, Lu-Huan Lou, Shui-Hu Jin, Kun Liu, Bin-Bin Liu

Data type: Images.

Explanation note: Species tree of Lycoris inferred from ASTRAL-III of the 78 plastid coding genes. Numbers above the branches indicate the branch support values measuring the support for a local posterior possibility.

This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.
Download file (2.05 MB)
login to comment