De novo assembly of the seed transcriptome and search for potential EST-SSR markers for an endangered, economically important tree species: Elaeagnus mollis Diels

Abstract

Seeds of Elaeagnus mollis (Elaeagnaceae) produce an edible oil and contain more vitamin E (Ve) than major oil-seed crops. Despite its economic value, there is no information on its genome sequence. Here, we used the Illumina platform to determine the seed transcriptome of E. mollis to identify the genes related to Ve biosynthesis and potential simple sequence repeat (SSR) markers. In total, 100,999 unigenes were obtained with an average length of 605 bp and N50 of 985 bp. Of these unigenes, 52,256 (51.7%) were annotated in at least one public database (NT, NR, PFAM, SwissProt, KOG, KO, and GO) in searches using blastn/x. The unigene annotation identified 15 unigenes encoding six enzymes (GGR, HPPD, HPT/VTE2, MPBQ-MT/VTE3, TC/VTE1, and γ-TMT/VTE4) putatively involved in Ve biosynthesis. In addition, 16,810 SSRs distributed in 14,057 unigenes were mined. Of these, 2820, 583, and 3423 SSRs were located in the 5′–UTR, coding sequence (CDS), and 3′–UTR regions, respectively, while the remaining 9984 SSRs had undetermined physical locations. The largest group of repeat motifs comprised mononucleotide repeats (70.76%), followed by dinucleotide (15.59%) and trinucleotide (12.10%) repeats. AG/CT (8.69%) and AAG/CTT (4.15%) were the main dinucleotide and trinucleotide repeats, respectively. Furthermore, 9597 SSR-specific primer pairs were designed. Among 100 primer pairs selected randomly to determine their usefulness, 53 proved to be efficient. To the our best of knowledge, this work is the first study of the E. mollis transcriptome and constitutes valuable genomics data for future genetic engineering studies to alter the amount of Ve. The identified potential EST-SSR markers can be used for population genetics studies and assisted-breeding of E. mollis.

Keywords

Elaeagnaceae / Elaeagnus mollis / Transcriptome / Vitamin E / Simple sequence repeats (SSR)

References

[1] Cao X, Zhang G. Application of the hyper-accurate mapping algorithm FANSe for next-generation sequencing in non-model organisms. Sci Sin Vitae, 2017, 47: 702-707.

[2] Chase MW, Christenhusz MJM, Fay MF, Byng JW, Judd WS, Soltis DE, Mabberley DJ, Sennikov AN, Soltis PS, Stevens PF. An update of the Angiosperm Phylogeny Group classification for the orders and families of flowering plants: APG IV. Bot J Linn Soc, 2016, 181: 1-20.

[3] Chen SY, Li HJ, Liu GS. Progress of vitamin E metabolic engineering in plants. Transgenic Res, 2006, 15: 655-665.

[4] Chen HL, Wang LX, Wang SH, Liu CJ, Blair MW, Cheng XZ. Transcriptome sequencing of mung bean (Vigna radiate L.) genes and the identification of EST-SSR markers. PLoS ONE, 2015, 10: e0120273.

[5] Chen LY, Cao YN, Yuan N, Nakamura K, Wang GM, Qiu YX. Characterization of transcriptome and development of novel EST-SSR makers based on next-generation sequencing technology in Neolitsea sericea (Lauraceae) endemic to East Asian land-bridge islands. Mol Breed, 2015, 35: 187.

[6] Dorp KV, Hölzl G, Plohmann C, Eisenhut M, Abraham M, Weber APM, Hanson AD, Dörmann P. Remobilization of phytol from chlorophyll degradation is essential for tocopherol synthesis and growth of Arabidopsis. Plant Cell, 2015, 27: 2846-2859.

[7] Du X, Hussain N, Li ZL, Chen XY, Hua SJ, Zhang DQ, Jiang LX. Effect of gibberellin on the biosynthesis of tocopherols in oilseed rape (Brassica napus L.) and Arabidopsis. J Agr Food Chem, 2015, 63: 360-369.

[8] Gailing O, Bodénès C, Finkeldey R, Kremer A, Plomion C. Genetic mapping of EST-derived simple sequence repeats (EST-SSRs) to identify QTL for leaf morphological characters in a Quercus robur full-sib family. Tree Genet Genomes, 2013, 9: 1361-1367.

[9] Georgiadou EC, Ntourou T, Goulas V, Manganaris GA, Kalaitzis P, Fotopoulos V. Temporal analysis reveals a key role for VTE5 in vitamin E biosynthesis in olive fruit during on-tree development. Front Plant Sci, 2015, 6: 871.

[10] Götz S, Garcíagómez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, Robles M, Talón M, Dopazo J, Conesa A. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res, 2008, 36: 3420.

[11] Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng QD, Chen ZH, Mauceli E, Hacohen N, Gnirke A, Rhind N, Palma F, Birren BW, Chad Nusbaum, Lindblad-Toh K, Friedman N, Regev A. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol, 2011, 29: 644-652.

[12] Harish MC, Dachinamoorthy P, Balamurugan S, Murugan SB, Sathishkumar R. Enhancement of α-tocopherol content through transgenic and cell suspension culture systems in tobacco. Acta Physiol Plant, 2013, 35: 1121-1130.

[13] Hu Z, Zhang T, Gao XX, Wang Y, Zhang Q, Zhou HJ, Zhao GF, Wang ML, Woeste KE, Zhao P. De novo assembly and characterization of the leaf, bud, and fruit transcriptome from the vulnerable tree Juglans mandshurica for the development of 20 new microsatellite markers using Illumina sequencing. Mol Genet Genom, 2016, 291: 849-862.

[14] Iseli C, Jongeneel CV, Bucher P. ESTScan: a program for detecting, evaluating, and reconstructing potential coding regions in EST sequences. Proc Int Conf Intell Syst Mol Biol, 1999, 99: 138-148.

[15] Kan LN, Wang L, Ding QZ, Wu YW, Ouyang J. Flash extraction and physicochemical characterization of oil from Elaeagnus mollis Diels seeds. J Oleo Sci, 2017, 66: 345-352.

[16] Kong SL, Abdullah SNA, Chai LH, Amiruddin MD. Molecular cloning, gene expression profiling and in silico sequence analysis of vitamin E biosynthetic genes from the oil palm. Plant Gene, 2016, 5: 100-108.

[17] Li DJ, Zhi D, Bi Q, Liu XH, Men ZH. De novo assembly and characterization of bark transcriptome using Illumina sequencing and development of EST-SSR markers in rubber tree (Hevea brasiliensis Muell. Arg.). BMC Genom, 2012, 13: 192.

[18] Liang SH, Yang RN, Dong CW, Yang QP. Physicochemical properties and fatty acid profiles of Elaeagnus mollis Diels nut oils. J Oleo Sci, 2015, 64: 1267-1272.

[19] Liu FL, Hu ZM, Liu WH, Li JJ, Wang WJ, Liang ZR, Sun XT. Distribution, function and evolution characterization of microsatellite in Sargassum thunbergii (Fucales, Phaeophyta) transcriptome and their application in marker development. Sci Rep, 2016, 6: 18947.

[20] Liu CC, Ying D, Guan XL, Fu Q, Zhang Z, Hu ZH, Zheng J, Lu YZ, Li W. De novo transcriptomic analysis and development of EST-SSRs for Sorbus pohuashanensis (Hance) Hedl. PLoS ONE, 2017, 12: e0179219.

[21] Liu HC, Cao GL, Han YP, Jiang ZF, Zhao HW, Li WB. Identification of the QTL underlying the vitamin E content of soybean seeds. Plant Breed, 2017, 136: 147-154.

[22] Mao XZ, Tao C, Olyarchuk JG, Wei LP. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics, 2005, 21(19): 3787-3793.

[23] Mène-Saffrané L, Pellaud S. Current strategies for vitamin E biofortification of crops. Curr Opin Biotech, 2017, 44: 189-197.

[24] Wang YL, Qin YY, Du Z, Yan GQ. Genetic diversity and differentiation of the endangered tree Elaeagnus mollis Diels (Elaeagnus L.) as revealed by simple sequence repeat (SSR) markers. Biochem Syst Ecol, 2012, 40: 25-33.

[25] Wang L, Yang Y, Zhao Y, Yang S, Udikeri S, Liu T. De Novo Characterization of the root transcriptome and development of EST-SSR markers in Paris polyphylla Smith var. yunnanensis, an endangered medical plant. J Agric Sci Technol, 2016, 18: 437-452.

[26] Wang PL, Yang LP, Zhang EH, Qin ZH, Wang HY, Liao YY, Wang XY, Gao L. Characterization and development of EST-SSR markers from a cold-stressed transcriptome of Centipedegrass by Illumina paired-end sequencing. Plant Mol Biol Rep, 2016, 35: 1-9.

[27] Wang WC, Chen SY, Zhang XZ. Characterization of the complete chloroplast genome of Elaeagnus mollis, a rare and endangered oil plant. Conserv Gene Resour, 2017, 9: 439-442.

[28] Wang XD, Guo K, Guo JL, Juan LI, Lu JC. Chemical constituents from leaves of Elaeagnus mollis. Chin Tradit Herbal Drugs, 2017, 48: 236-240.

[29] Xing W, Liao JY, Cai MY, Xia QF, Liu Y, Zeng W, Jin XL. De novo assembly of transcriptome from Rhododendron latoucheae Franch using Illumina sequencing and development of new EST-SSR markers for genetic diversity analysis in Rhododendron. Tree Genet Genom, 2017, 13: 53.

[30] Yan GQ, Tian LH, Yang LY. Studies on the browning in callus induction of Elaeagnus mollis. Acta Bot Boreal-Occident Sin, 2004, 24: 1384-1389.

[31] Yan LP, Liu CL, Wu DJ, Li L, Shu J, Sun C, Xia Y, Zhao LJ. De novo transcriptome analysis of Fraxinus velutina using Illumina platform and development of EST-SSR markers. Biol Plant, 2016, 61: 210-218.

[32] Yang W, Cahoon RE, Hunter SC, Zhang C, Han J, Borgschulte T, Cahoon EB. Vitamin E biosynthesis: functional characterization of the monocot homogentisate geranylgeranyl transferase. Plant J Cell Mol Biol, 2011, 65: 206-217.

[33] Ye ZY, Wang ZS, Li YX, Rao GD, Zhang JG. Development and characterization of microsatellite markers in the Elaeagnus mollis Diels. Acta Bot Boreal-Occident Sin, 2016, 36: 274-279.

[34] Zhang JM, Zhang F. Population structure and genetic variation of the endangered species Elaeagnus mollis Diels (Elaeagnaceae). Genet Mol Res, 2015, 14: 5950-5957.

[35] Zhang LW, Li YR, Tao AF, Fang PP, Qi JM. Development and characterization of 1906 EST-SSR narkers from unigenes in jute (Corchorus spp.). PLoS ONE, 2015, 10: e0140861.

[36] Zhou Q, Luo D, Ma LC, Xie WG, Wang Y, Wang YR, Liu ZP. Development and cross-species transferability of EST-SSR markers in Siberian wildrye (Elymus sibiricus L.) using Illumina sequencing. Sci Rep, 2016, 6: 20549.