Skip to main content
  • Research Article
  • Open access
  • Published:

Comparative transcriptomics provide insight into the morphogenesis and evolution of fistular leaves in Allium

Abstract

Background

Fistular leaves frequently appear in Allium species, and previous developmental studies have proposed that the process of fistular leaf formation involves programmed cell death. However, molecular evidence for the role of programmed cell death in the formation of fistular leaf cavities has yet to be reported.

Results

In this study, we characterized the leaf transcriptomes of nine Allium species, including six fistular- and three solid-leaved species. In addition, we identified orthologous genes and estimated their Ka and Ks values, in order to ascertain their selective pattern. Phylogenetic analysis based on the transcriptomes revealed that A. tuberosum was the most ancestral among the nine species, and analysis of orthologous genes between A. tuberosum and the other eight species indicated that 149 genes were subject to positive selection; whereas >3000 had undergone purifying selection in each species.

Conclusions

We found that many genes that are potentially related to programmed cell death either exhibited rapid diversification in fistular-leaved species, or were conserved in solid-leaved species in evolutionary history. These genes potentially involved in programmed cell death might play important roles in the formation of fistular leaf cavities in Allium, and the differing selection patterns in fistular- and solid-leaved species may be responsible for the evolution of fistular leaves.

Background

As one of the largest genera of the petaloid monocotyledons, Allium (Amaryllidaceae) comprises more than 920 species [1] and includes several economically important crops that are cultivated for consumption or medicinal uses, such as garlic (A. sativum), welsh onion (A. fistulosum), leek (A. porrum), Chinese chives (A. tuberosum), onion (A. cepa), Chinese jiaotou (A. chinense), and shallot (A. ascalonicum). Varied leaf shapes can be observed among Allium species, including flat, columnar, solid, and fistular morphologies. Morphological and cellular studies have found that fistular leaves develop from solid precursors [2]. Developmental investigation of the leaves of A. fistulosum found that the process of fistular leaf formation involved programmed cell death (PCD) [2]. PCD is a spontaneous, programmed, self-destructive cellular process that plays a key role in tissue differentiation, homeostasis, and organ morphogenesis, including that of leaves [35]. However, molecular evidence for the involvement of PCD in the formation of fistular leaf cavities is still absent in Allium, owing to the limited availability of genetic resources for Allium spp.

The paucity of genetic resources in Allium is, in part, due to the fact that Allium spp. have the largest genomes among eukaryotes [6, 7], with genome sizes ranging from 6860 to 30,870 Mbp per C [8]. The enormous size of these genomes has been a major obstacle for their characterization and for gene mining in the family as a whole. In the past 10 years, the next generation sequencing technologies have undergone rapid development, and more than 80 plant species have had their complete genome [9]. However, none of genome of Allium species has been characterized, because of their enormous size.

Because transcriptome analysis by next generation sequencing is rapid, inexpensive, and unconstrained by genomic complexity, it has been widely used as a primary tool for gene discovery and expression profiling in hundreds of plant species [10, 11]. Moreover, transcriptome analysis can also be used as an important tool for investigating the domesticated patterns of crops [12, 13], as well as for investigating the mechanisms of development for specific traits [14]. In Allium, the transcriptome of four species, i.e. garlic [15, 16], onion [17], Chinese chive [18], and welsh onion [19], have been sequenced and de novo assembled. Among these, onion and Welsh onion are fistular-leaved species, whereas the other two produce solid leaves. However, despite the large numbers of generated expressed sequence tags, the molecular mechanism for the formation of fistular leave cavities is still uncharacterized.

Therefore, we analyzed the transcriptomes of nine economically important Allium species, including eight vegetable species and one herbal species (A. macrostemon). Thereafter, genes that had undergone significant selection were identified in both fistular- and solid-leaf species, respectively, and the genes with evolutionary divergence between fistular- and solid-leaf species were screened. As a result, the current study presents molecular evidence regarding a potential evolutionary mechanism for the development of fistular leaves.

Methods

Plant material and RNA extraction

The current study utilized local varieties of nine economically important Allium species. Among the eight vegetable species, three [A. sativum (garlic; SAT), A. porrum (leek; POR), and A. tuberosum (Chinese chives; TUB)] possess flat, solid leaves, whereas four species, A. fistulosum (welsh onion; FIS), A. ascalonicum (shallot; ASC), A. cepa (onion; CEP), and A. cepa var. agrogarum (AGR), have cylindrical, fistular leaves, and A. chinense (Chinese jiaotou; CHI) has triangular, fistular leaves (Fig. 1). In addition, the leaves of the herbal species A. macrostemon (MAC) are also fistular, but internal cavities are very small (Fig. 1). The varieties of SAT, CEP, and MAC were collected from Chaling (Hunan, China), whereas those of ASC, POR, TUB, and CHI were collected from Ningxiang (Changsha, China), and the varieties of FIS and AGR were collected from Yuanjiang (Hunan, China) and Fuyu (Jinlin, China), respectively. The transverse section of leaves of nine Allium species was observed by Nikon AZ100 microscope (Nikon, Toyota, Japan).

Fig. 1
figure 1

The transverse section picture of leaves of nine Allium species

All the varieties were established in the experimental field of the Institute of Bast Fiber Crops, Chinese Academy of Agricultural Sciences, Changsha, China on Sept. 15, 2014, and on Mar. 10, 2015, leaf tissue was sampled from three individuals of each species, immediately frozen in liquid nitrogen, and stored at −80°C until used. The total RNA of each sample was extracted using an EZNA. Plant RNA Kit (OMEGA Bio-Tek, Norcross, GA, USA), according to the manufacturer’s protocol.

cDNA library construction, sequencing, and assembly

Illumina sequencing for the nine species was performed at Novogene Bioinformatics Technology Co., Ltd (Beijing, China; www.novogene.com). First, the total RNA of each sample was used to construct a cDNA library with fragment lengths of ~250 bp. Thereafter, paired-end sequencing was performed using the Illumina HiSeq 2500 sequencing platform (Illumina, San Diego, CA, USA), and the clean sequencing reads of each species were deposited in the NCBI Sequence Read Archive (Additional file 1: Table S1). After trimming the adapter sequences and filtering low-quality reads, the reads of each species were used to assemble their transcriptomes de novo using Trinity [20]. Sequences not extended on either end were categorized as unigenes. The sequences of unigenes assembled are deposited in the NCBI Transcriptome Shotgun Assembly (TSA) Database (Additional file 1: Table S1).

Gene annotation

To annotate the unigenes of each species, the corresponding sequences were searched against public databases, including the NCBI non-redundant protein sequences (NR) database, NCBI nucleotide sequences (NT) database, eukaryotic ortholog groups (KOG) database, Kyoto Encyclopedia of Genes and Genomes ortholog (KO) database, Swiss-Prot protein database, Gene Ontology (GO) database, and protein family (PFAM) database.

Ortholog search and construction of phylogenetic tree

The OrthoMCL method, which is based on the Markov cluster algorithm [21], was used to identify for orthologous genes among the nine species. According to the orthologous gene sequences, the phylogenetic relationships of nine species were reconstructed using Phyml 3.0, which is based on a Bayes algorithm [22], and the tree was visualized using MEGA 4.0 [23].

Ka/Ks analysis

Positive selection is a process that favors the retention of beneficial mutations, as well as the rapid diversification of the affected genes. If a gene has evolved under positive selection, the sequence is expected to contain more non-synonymous nucleotide substitutions (Ka) than synonymous nucleotide substitutions (Ks) and, accordingly, should also exhibit a high Ka/Ks ratio. In contrast, purifying selection is a process that removes deleterious alleles, and the sequences of genes subjected to purifying selection are relatively conserved. Accordingly, these sequences contain fewer, if any, non-synonymous substitutions, which results in a Ka that is much smaller than Ks, and a low Ka/Ks ratio. To identify genes subjected to significant selection, the Ka and Ks values were estimated using the Codeml program of Phylogenetic Analysis by Maximum Likelihood (PAML) based on the basic model [24], and genes with Ka/Ks ratios >1 were considered as under positive selection, whereas genes with Ka/Ks ratios <0.1 were considered as under purifying selection [25].

Enrichment analysis

The enrichment of GO functional categories in the selected genes of each species were analyzed using GOseq, which is based on the Wallenius non-central hyper-geometric distribution [26], and enrichment analysis of the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways was performed using KOBAS [27]. Q values were used for determining the P-value threshold in multiple test and analysis [28], and pathways or GO categories with Q < 0.05 were considered significantly enriched.

Results

Transcriptome sequencing, assembly, and annotation

Illumina paired-end sequencing of leaf transcriptomes from nine Allium spp. yielded about 99.9–128.9 million clean sequencing reads with a mean length of ~125 bp for each species (Table 1). Except for AGR, in which only 83,186 unigenes were de novo assembled, all species yielded >117 thousand unigenes, with an average length of 529–641 bp. POR had the largest transcriptome, at 109.28 Mb, whereas the transcriptomes of AGR and CHI were relatively small, at 53.35 and 66.84 Mb, respectively, and the other six species possessed intermediate sized transcriptomes, ranging from 71.62 to 85.48 Mb (Table 1). In addition, except for the unigenes of POR and AGR, of which 25.62 and 41.04% were successfully annotated, respectively, ~30–35% of unigenes were typically annotated.

Table 1 Summary of transcriptome sequencing for nine Allium species

Phylogenetic analysis

According to our phylogenetic tree (Fig. 2), TUB was located at the root of the tree, and had far phylogenetic distance with the others, suggesting that TUB was more ancestral than the other species. This finding was consistent with the result of previous study [29]. The other eight species were split into four clades. Interestingly, the classification based on the transcriptomes was consistent with that based on leaf morphology, i.e., the species with flat, solid leaves (SAT and POR) were grouped together, the species with cylindrical, fistular leaves (ASC, FIS, CEP, and AGR) were placed in a monophyletic clade, and the species with small-gapped and triangular fistular leaves (MAC and CHI) were placed in a single branch, as well (Fig. 2).

Fig. 2
figure 2

Phylogenetic tree of nine Allium species, constructed on the basis of their transcriptomes

Identification of genes subjected to significant selection

A total of 29,990 orthologous groups were identified, with 9634 from the two flat, solid-leaved species (SAT and POR; Fig. 3a), 6995 from the four cylindrical, fistular-leaved species (ASC, FIS, CEP, and AGR; Fig. 3b), 3720 among all nine species (Fig. 3c).

Fig. 3
figure 3

Venn diagram of the orthologous genes in nine Allium species. a The orthologs of two solid-leaved species: SAT and POR; b the orthologs of four columnar fistular-leaved species; c the orthologs of all nine species, in which SAT_POR and ASC_CEP_AGR_FIS represent the orthologs of SAT and POR and the orthologs of ASC, CEP, AGR, and FIS, respectively; d the orthologs of TUB and the other eight species

Because the TUB transcriptome was ancestral to the transcriptomes of the other eight species, it was designated as the ancestral species in the present study. The orthologous genes between TUB and the other eight species were searched (Fig. 3d), and their Ka and Ks values were estimated and used to determine the selective pattern of them (Fig. 4). A total of 149 genes were found to undergo positive selection with 13 to 27 identified in each species (Fig. 5), suggesting that these 149 genes exhibited a rapid diversification in evolutionary history. Among these 149 positively selected genes, 39 genes (from 18 orthologous groups) were mitochondrial (Additional file 2: Table S2), 29 genes (from 13 orthologous groups) were chloroplast (Additional file 3: Table S3), and 81 genes (from 70 orthologous groups) were nuclear (Additional file 4: Table S4). In addition, >3000 genes were found to undergo purifying selection in each species (Fig. 5), which suggested that the sequences of these genes were relatively conserved, in evolutionary history.

Fig. 4
figure 4

Scatter diagrams of Ka and Ks values. The Ka and Ks values of the orthologs of TUB and the other eight species were estimated. The x and y axes represent the Ka and Ks values, respectively. Orthologous genes with Ka > 1 or Ks > 4 are not shown

Fig. 5
figure 5

The gene number subjected to positive and purifying selection in each species

Enrichment analysis for selected genes

None of the GO categories were enriched by positively selected genes in any species, although 51–72 terms were significantly enriched by genes subjected to purifying selection in each species (Q < 0.05; Additional file 5: Table S5). Among the enriched terms, 23 were enriched by genes of purifying selection in all eight species, suggesting that the functional terms are essential for growth in the species. In addition, one KEGG pathway, oxidative phosphorylation, exhibited significant enrichment by positively selected genes in all species, except FIS (Q < 0.05); however, no pathways were markedly enriched by genes of purifying selection (Q < 0.05).

Identification of selected genes potentially related to PCD

Among 149 identified positively selected genes, 25 from 11 orthologous groups were annotated as proteins involved in plant PCD (Table 2), including 15 genes related to the metabolism of cytochrome c, which is an important factor for inducing PCD [30], and two genes encoding PCD-related kinase proteins. In addition, there were 4970 orthologous groups in which the orthologous gene of at least two species had undergone purifying selection. Of the 4,970 orthologous groups, there were 116 groups in which only the orthologous genes from the flat, solid-leaved species(SAT and POR) were subjected to purifying selection (Additional file 6: Table S6), suggesting that the genes are possibly related to leaf morphology. Among the 116 orthologous groups, five were PCD-related (Table 3), especially OG01095 that encodes a metacaspase-like protein, which is a key regulator of plant PCD [31].

Table 2 Positively selected genes potentially related to programmed cell death
Table 3 Programmed cell death-related orthologous groups in which only the orthologous genes from the flat, solid-leaved species A. sativum and A. porrum were subjected to purifying selection

Discussion

Although Allium comprises a large number of species and even includes a few species of significant economic importance, the transcriptomes of only four species (i.e. garlic, onion, welsh onion, and Chinese chive) have been analyzed [1519]. Therefore, in the present study, the transcriptomes of nine Allium species were sequenced and assembled de novo, and large numbers of genes were identified in each species. These results will help us investigate the major genes responsible for important agronomic traits in Allium species, and will also further our understanding of the mechanisms underlying the formation and development of such traits, including leaf morphology.

Leaves are the primary site of plant photosynthesis and respiration. Diverse leaf forms can be observed in angiosperms, and they can be categorized into two major categories: bifacial or unifacial. Development studies have proposed that the diversity of leaf morphology is caused by differential and asymmetrical growth along the adaxial–abaxial, central–marginal, and proximal–distal polarities [3234]. Allium spp. exhibit a wide variety of leaf shapes, including cylindrical and hollow, flat and solid, and hollow and triangular types, the differences of which can be observed in transverse sections [35]. As in other angiosperms, there is a great possibility that the diversity of leaf appearance in Allium is caused by these three polarities. For example, only the abaxial identity existed in leaf blades results in a unifacial fistular leaf in Welsh onion [33, 34] However, the formation of fistular leaf cavities is still poorly understood. Although that a cellular investigation of the leaves of A. fistulosum found that the process of fistular leaf formation involved PCD [2], molecular evidences for the involvement of PCD in the formation of fistular leaf cavities is still unknown in Allium.

Previous studies have shown that mitochondria play a central role in controlling PCD in plants [30], by releasing molecules, including cytochrome c that drive the destruction of the cell [36]. In animals, cytochrome c triggers PCD by forming a complex that induces the formation of apoptosomes, and those apoptosomes then activate downstream effector caspases, which are key regulators of PCD [30]. Although no caspase homologs have been found in plants, caspase-like proteins have been described in cells undergoing PCD, such as the cysteine endopeptidases, which can be subdivided into two groups: vacuolar processing enzymes and metacaspases [37, 38]. In the present study, we identified 15 positively selected genes from five orthologous groups that were potentially related to the metabolism of cytochrome c and were found to undergo positive selection. Among these 15 genes, 12 were subjected positive selection in fistular-leaved species, except 3 genes from OG13835 that underwent positive selection in the species with solid or small-gapped fistular leaves. In addition, a KDEL-tailed cysteine endopeptidase-encoding gene from the fistular-leaved ASC exhibited rapid evolution, and a metacaspase-4-like protein-encoding gene was found to undergo positive selection in solid-leaved species, but evolved neutrally in fistular-leaved species.

Like the mitochondria, the chloroplasts are also important in triggering PCD, by mediating reactive oxygen species signaling. In plants, photosynthesis is performed in photosystem II of chloroplasts. When photosynthesis is less efficient, or chloroplasts receive excessive levels of light, chloroplasts produce more reactive oxygen species, and then the reactive oxygen species drives cell death finally [30, 39, 40]. In the present study, we identified five chloroplast genes that are potentially involved in photosynthesis in photosystem II and exhibited rapid diversification. Among them, four of the genes were from fistular-leaved species, and one was from a flat, solid-leaved species (i.e., SAT). In addition, two photosystem II protein-encoding genes were also found to have undergone purfying selection in solid-leaved species, whereas they were evolving neutrally in fistular-leaved species.

After the “death signal” is generated by mitochondria and/or chloroplasts, it is transduced by kinase proteins. Previous studies have shown that PCD signaling involves the activation of mitogen-activated protein kinases that induce hypersensitive responses and cell apoptosis [41]. Recently, the phosphatidyl-3 kinases ATM and ATR were reported as the main sensors of DNA damage-induced PCD [42, 43]. In the present study, the ATM and ATR kinase-encoding genes in the fistular-leaved species FIS, CEP, and CHI exhibited evidence of rapid diversification, but not in solid-leaved species, and mitogen-activated protein kinase-encoding gene underwent positive selection in solid-leaved species, whereas it was evolving neutrally in fistular-leaved species.

In addition, two genes, LSD-one-like protein 2 (LOL2) and constitutively activated cell death 1 (CAD1), have also been reported to play important roles in regulating PCD in Arabidopsis [44, 45]. In the present study, the orthologs of these two genes were also found to exhibit different selective patterns among the fistular- and solid-leaved species. These potential PCD-related genes with evolutionary divergence might have roles in regulating the formation of cavity of fistular leaf in Allium species.

Conclusion

In this study, we compared the leaf transcriptomes of six fistular- and three solid-leaf species of Allium. The selective patterns of orthologous genes in nine species were analyzed, and the result showed that many mitochondrial, chloroplast and nuclear genes that were potentially involved in plant PCD and were either subject to positive selection in fistular-leaved species but not in solid-leaved species, or evolved neutrally in fistular-leaved species but purifying selection in solid-leaved species. These potential PCD-related genes might play roles in the formation of the inner cavities of fistular leaves in Allium, and the differing selection patterns in fistular- and solid-leaved species may be responsible for the evolution of fistular leaves.

Abbreviations

AGR:

A. cepa var. agrogarum

ASC:

A. ascalonicum

CEP:

A. cepa

CHI:

A. chinense

FIS:

A. fistulosum

GO:

Gene Ontology database

Ka:

Non-synonymous nucleotide substitutions

KO:

KEGG ortholog database

KOG:

Eukaryotic ortholog groups database

Ks:

Synonymous nucleotide substitutions

MAC:

A. macrostemon

NR:

NCBI non-redundant protein sequencesdatabase

NT:

NCBI nucleotide sequences database

PCD:

Programmed cell death

PFAM:

Protein family database

POR:

A. porrum

SAT:

A. sativum

TUB:

A. tuberosum

References

  1. Herden T, Hanelt P, Friesen N. Phylogeny of Allium L. subgenus Anguinum (G. Don. ex W.D.J. Koch) N. Friesen (Amaryllidaceae). Mol Phylogenet Evol. 2016;95:79–93.

    Article  PubMed  Google Scholar 

  2. Ni X, Su H, Zhou Y, Wang F, Liu W. Leaf-shape remodeling: programmed cell death in fistular leaves of Allium fistulosum. Physiol Plant. 2015;153:419–31.

    Article  CAS  PubMed  Google Scholar 

  3. Meier P, Finch A, Evan G. Apoptosis in development. Nature. 2000;407:796–801.

    Article  CAS  PubMed  Google Scholar 

  4. Gunawardena A, Greenwood J, Dengler N. Programmed cell death remodels lace plant leaf shape during development. Plant Cell. 2004;16:60–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Gunawardena A, Sault K, Donnelly P, Greenwood J, Dengler N. Programmed cell death and leaf morphogenesis in Monstera obliqua (Araceae). Planta. 2005;221:607–18.

    Article  CAS  PubMed  Google Scholar 

  6. Labani R, Elkington T. Nuclear DNA variation in the genus Allium L. (Liliaceae). Heredity. 1987;59:119–28.

    Article  Google Scholar 

  7. Baranyi M, Greilhuber J. Genome size in Allium: in quest of reproducible data. Ann Bot. 1999;83:687–95.

    Article  Google Scholar 

  8. Ricroch A, Yockteng R, Brown SC, Nadot S. Evolution of genome size across some cultivated Allium species. Genome. 2005;48:511–20.

    Article  CAS  PubMed  Google Scholar 

  9. Kitts P, Church D, Thibaud-Nissen F, Choi J, Hem V, Sapojnikov V, Smith R, Tatusova T, Xiang C, Zherikov A, et al. Assembly: a resource for assembled genomes at NCBI. Nucleic Acids Res. 2016;44:D73–80.

    Article  PubMed  Google Scholar 

  10. Liu T, Zhu S, Tang Q, Chen P, Yu Y, Tang S. De novo assembly and characterization of transcriptome using Illumina paired-end sequencing and identification of CesA gene in ramie (Boehmeria nivea L. Gaud). BMC Genomics. 2013;14:125.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Mutz K, Heilkenbrinker A, Lonne M, Walter J, Stahl F. Transcriptome analysis using next-generation sequencing. Curr Opin Biotechnol. 2013;24:22–30.

    Article  CAS  PubMed  Google Scholar 

  12. Liu T, Tang S, Zhu S, Tang Q, Zheng X. Transcriptome comparison reveals the patterns of selection in domesticated and wild ramie (Boehmeria nivea L. Gaud). Plant Mol Biol. 2014;86:85–92.

    Article  CAS  PubMed  Google Scholar 

  13. Koenig D, Jiménez-Gómez J, Kimura S, Fulopa D, Chitwooda D, Headlanda L, Kumara R, Covingtona M, Devisettya U, Tata A, et al. Comparative transcriptomics reveals patterns of selection in domesticated and wild tomato. Proc Natl Acad Sci U S A. 2013;110:E2655–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Lou Q, Liu Y, Qi Y, Jiao S, Tian F, Jiang L, Wang Y. Transcriptome sequencing and metabolite analysis reveals the role of delphinidin metabolism in flower colour in grape hyacinth. J Exp Bot. 2014;65:3157–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Kamenetsky R, Faigenboim A, Mayer E, Michael T, Gershberg C, Kimhi S, Esquira I, Shalom S, Eshel D, Rabinowitch H, et al. Integrated transcriptome catalogue and organspecific profiling of gene expression in fertile garlic (Allium sativum L.). BMC Genomics. 2015;16:12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Liu T, Zeng L, Zhu S, Chen X, Tang Q, Mei S, Tang S. Large-scale development of expressed sequence tag-derived simple sequence repeat markers by deep transcriptome sequencing in garlic (Allium sativum L.). Mol Breed. 2015;35:204.

    Article  Google Scholar 

  17. Rajkumar H, Ramagoni R, Anchoju V, Vankudavath R, Syed A. De novo transcriptome analysis of Allium cepa L. (Onion) bulb to identify allergens and epitopes. PLoS One. 2015;10:e0135387.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Zhou S, Chen L, Liu S, Wang X, Sun X. De novo assembly and annotation of the Chinese chive (Allium tuberosum Rottler ex Spr.) transcriptome using the Illumina platform. PLoS One. 2015;10:e0133312.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Sun X, Yu X, Zhou S, Liu S. De novo assembly and characterization of the Welsh onion (Allium fistulosum L.) transcriptome using Illumina technology. Mol Genet Genomics. 2016;291:647–59.

    Article  CAS  PubMed  Google Scholar 

  20. Grabherr M, Haas B, Yassour M, Levin J, Thompson D, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, et al. Full length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Li L, Stoeckert Jr C, Roos D. OrthoMCL: Identification of ortholog groups for eukaryotic genomes. Genome Res. 2003;13:2178–89.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Guindon S, Dufayard J, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59:307–21.

    Article  CAS  PubMed  Google Scholar 

  23. Tamura K, Dudley J, Nei M, Kumar S. MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007;24:1596–9.

    Article  CAS  PubMed  Google Scholar 

  24. Yang Z. PAML 4: Phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007;24:1586–91.

    Article  CAS  PubMed  Google Scholar 

  25. Shah C, VanGompel MJW, Naeem V, Chen Y, Lee T, Angeloni N, Wang Y, Xu E. Widespread presence of human BOULE homologs among animals and conservation of their ancient reproductive function. PLoS Genet. 2010;6:e1001022.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Young M, Wakefield M, Smyth G, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11:R14.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, Kong L, Gao G, Li C, Wei L. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011;39:W316–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc B. 1995;57:289–300.

    Google Scholar 

  29. Friesen N, Fritsch R, Blattner F. Phylogeny and new intrageneric classification of Allium (Alliaceae) based on nuclear ribosomal DNA ITS sequences. Aliso. 2006;22:372–95.

    Google Scholar 

  30. Reape T, Brogan N, McCabe P. Mitochondrion and chloroplast regulation of plant programmed cell death. In: Gunawardena A, McCabe P, editors. Plant Programmed Cell Death. New York: Springer; 2015. p. 33–53.

    Chapter  Google Scholar 

  31. Lam E, Zhang Y. Regulating the reapers: activating metacaspases for programmed cell death. Trends Plant Sci. 2012;17:487–94.

    Article  CAS  PubMed  Google Scholar 

  32. Yamaguchi T, Tsukaya H. Evolutionary and developmenta studies of unifacial leaves in monocots: Juncus as a model system. J Plant Res. 2010;123:35–41.

    Article  PubMed  Google Scholar 

  33. Yamaguchi T, Nukazuka A, Tsukaya H. Leaf adaxial–abaxial polarity specification and lamina outgrowth: evolution and development. Plant Cell Physiol. 2012;53:1180–94.

    Article  CAS  PubMed  Google Scholar 

  34. Yamaguchi T, Yano S, Tsukaya H. Genetic framework for flattened leaf blade formation in unifacial leaves of Juncus prismatocarpus. Plant Cell. 2010;22:2141–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Lee SH. A comparative anatomical study on the leaves of some bulb crops in China. J Integr Plant Biol. 1952;1:1–7.

    Google Scholar 

  36. Diamond M, McCabe P. Mitochondrial regulation of plant programmed cell death. In: Kempken F, editor. Plant mitochondria, vol 1, Advances in plant biology. New York: Springer; 2011. p. 439–65.

    Google Scholar 

  37. Helm M, Schmid M, Hierl G, Terneus K, Tan L, Lottspeich F, Kieliszewski M, Gietl C. KDEL-tailed cysteine endopeptidases involved in programmed cell death, intercalation of new cells, and dismantling of extension scaffolds. Am J Bot. 2008;95:1049–62.

    Article  CAS  PubMed  Google Scholar 

  38. Lord C, Dauphinee A, Watts R, Gunawardena. A unveiling interactions among mitochondria, caspase-Like proteases, and the actin cytoskeleton during plant programmed cell death (PCD). PLoS ONE. 2013;8:e57110.

  39. Kim C, Meskauskiene R, Zhang S, Lee K, Ashok M, Blajecka K, Herrfurth C, Feussner I, Apela K. Chloroplasts of Arabidopsis are the source and a primary target of a plant-specific programmed cell death signaling pathway. Plant Cell. 2012;24:3026–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Aken O, Breusegem F. Licensed to kill: mitochondria, chloroplasts, and cell death. Trends Plant Sci. 2015;20:754–66.

    Article  PubMed  Google Scholar 

  41. Kuriyama H, Fukuda H. Developmental programmed cell death in plants. Curr Opin Plant Biol. 2002;5:568–73.

    Article  CAS  PubMed  Google Scholar 

  42. Amiard S, Depeiges A, Allain E, White C, Gallego M. Arabidopsis ATM and ATR kinases prevent propagation of genome damage caused by telomere dysfunction. Plant Cell. 2011;23:254–4265.

    Article  Google Scholar 

  43. Smetana O, Siroky J, Houlne G, Opatrny Z, Chaboute M. Non-apoptotic programmed cell death with paraptotic-like features in bleomycin-treated plant cells is suppressed by inhibition of ATM/ATR pathways or NtE2F overexpression. J Exp Bot. 2012;63:2631–44.

    Article  CAS  PubMed  Google Scholar 

  44. Morita-Yamamuro C, Tsutsui T, Sato M, Yoshioka H, Tamaoki M, Ogawa D, Matsuura H, Yoshihara T, Ikeda A, Uyeda I, et al. The Arabidopsis gene CAD1 controls programmed cell death in the plant immune system and encodes a protein containing a MACPF domain. Plant Cell Physiol. 2005;46:902–12.

    Article  CAS  PubMed  Google Scholar 

  45. Coll N, Epple P, Dangl J. Programmed cell death in the plant immune system. Cell Death Differ. 2011;18:1247–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We acknowledge the Novogene Bioinformatics Institute for its assistance in original data processing and related bioinformatics analysis.

Funding

This work was supported by grants from the Agricultural Science and Technology Innovation Program of China (CAAS-ASTIP-IBFC). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Availability of data and materials

The sequencing reads is available in the NCBI SRA database under accession SRX1560526, SRX1560538, SRX1560563, SRX1560658, SRX1560673, SRX1560692, SRX1560712, SRX1560726, and SRX1560727; the transcripts assembled are deposited in the NCBI Transcriptome Shotgun Assembly (TSA) Database under accession number GEOY00000000, GFAP00000000, GFAR00000000, GFAN00000000, GFAO00000000, GFAL00000000, GFAJ00000000, GFAK00000000, GFAM00000000.

Authors’ contributions

SZ and QD performed the experiment. SZ revised the manuscript, and submitted the data to database. ST planned and designed the research. YY conducted the fieldwork, and ZT prepared the reagents of experiments. TL analyzed the data and wrote the manuscript. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Touming Liu.

Additional files

Additional file 1: Table S1.

Total length of clean sequencing reads from each species and their accession numbers in database. (DOCX 16 kb)

Additional file 2: Table S2.

Positively selected mitochondrial genes in evolution. (DOCX 13 kb)

Additional file 3: Table S3.

Positively selected chloroplast genes in evolution. (DOCX 13 kb)

Additional file 4: Table S4.

Positively selected nuclear genes in evolution. (XLSX 14 kb)

Additional file 5: Table S5.

Table S5 Gene ontology (GO) terms enriched by genes of purifying selection in each species. (XLSX 16 kb)

Additional file 6: Table S6.

Orthologous groups in which only the orthologous genes from the flat, solid-leaved species A. sativum and A. porrum were subjected to purifying selection. (XLSX 17 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhu, S., Tang, S., Tan, Z. et al. Comparative transcriptomics provide insight into the morphogenesis and evolution of fistular leaves in Allium . BMC Genomics 18, 60 (2017). https://doi.org/10.1186/s12864-016-3474-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-016-3474-8

Keywords