Development of SSR Markers Associated with Biosynthesis Pathway of Steviol Glycosides in Stevia through De Novo Transcriptome Assembly

Khayam Nekoui, Mojtaba; Moazam Jazi, Maryam; Mardi, Mohsen; Kadkhodaei, Saeid

Development of SSR Markers Associated with Biosynthesis Pathway of Steviol Glycosides in Stevia through De Novo Transcriptome Assembly

Document Type : Original Research

Authors

Mojtaba Khayam Nekoui ¹

Maryam Moazam Jazi ²

Mohsen Mardi ³

Saeid Kadkhodaei ⁴

¹ Tarbiat Modares University

² Shahid Beheshti Medical University

³ Seed and Plant Certification and Registration Institute

⁴ Isfahan University of Technology

Abstract

In stevia (Stevia rebaudiana), breeding programs are mainly aimed at developing plants with high Rebaudioside-A (RA) content. To this end, in order to screen stevia plants and selection of varieties with the highest amount of desired sweeteners (RA) using molecular markers, the present study was conducted on RNA-seq data of varieties having different amounts of RA. We took advantage of CLC to make de novo transcriptome assembly for each variety with k-mer and contig length values of 20 and 200bp, respectively. The assembly was annotated using the latest Arabidopsis proteome release. To identify signatures of candidate polymorphic SSRs among the stevia varieties, the assembled sequences were used as an input for CandiSSR, followed by designing primer pairs for identified polymorphic SSRs. We identified 368 potential polymorphic SSRs based on the stevia transcriptome analysis, among which 360 were qualified for primer design. Almost 89% of the contig sequences possessing polymorphic SSRs had the best blast hit against Arabidopsis proteome. We found contigs similar to the UDP-Glycosyltransferase protein family and Deoxyxylulose-5-phosphate synthase which are involved in biosynthesis pathway of steviol glycosides. Also, gene set enrichment analysis using PlantGSE through Hypergeometric test (FDR<0.05) identified enriched metabolic pathways in the sequences contained polymorphic SSRs; It is therefore most likely that such connections exist between the SSRs and biosynthesis of steviol glycosides. Hence, it could conceivably be hypothesized that the SSR markers developed in this study would be reliable in molecular breeding of stevia toward selection of varieties with high RA content.

Keywords

De novo transcriptome assembly

Marker development

SSR

Stevia

Steviol glycoside biosynthesis pathway

UDP-Glycosyltransferase

20.1001.1.23222115.1399.11.2.12.3

Subjects

Bioinformatics

[1] Yadav, A. K., Singh, S., Dhyani, D., Ahuja, P. S. (2011). A review on the improvement of stevia [Stevia rebaudiana (Bertoni)]. Canadian Journal of Plant Science. 91(1), 1-27.
[2] Soufi, S., D’Urso, G., Pizza, C., Rezgui, S., Bettaieb, T. & Montoro, P. (2016). Steviol glycosides targeted analysis in leaves of Stevia rebaudiana (Bertoni) from plants cultivated under chilling stress conditions. Food chemistry, 190,572-580.
[3] Wölwer-Rieck, U. (2012). The leaves of Stevia rebaudiana (Bertoni), their constituents and the analyses thereof: a review. Journal of agricultural and food chemistry, 60(4), 886-895.
[4] Zalapa, J. E., Cuevas, H., Zhu, H., Steffan, S., Senalik, D., Zeldin, E., McCown, B., Harbut, R. & Simon, P. (2012). Using next-generation sequencing approaches to isolate simple sequence repeat (SSR) loci in the plant sciences. American journal of botany, 99(2), 193-208.
[5] Davey, J. W., Hohenlohe, P. A., Etter, P. D., Boone, J. Q., Catchen, J. M., & Blaxter, M. L. (2011). Genome-wide genetic marker discovery and genotyping using next-generation sequencing. Nature Reviews Genetics, 12(7),499.
[6] Zhang, H., Wei, L., Miao, H., Zhang, T. & Wang, C. (2012). Development and validation of genic-SSR markers in sesame by RNA-seq. BMC genomics, 13(1), 316.
[7] Huang, H., Xia, E.-H., Zhang, H.-B., Yao, Q.-Y., & Gao, L.-Z. (2017). De novo transcriptome sequencing of Camellia sasanqua and the analysis of major candidate genes related to floral traits. Plant Physiology and Biochemistry, 120, 103-111.
[8] Chen, J., Hou, K., Qin, P., Liu, H., Yi, B., Yang, W., & Wu, W. (2014). RNA-Seq for gene identification and transcript profiling of three Stevia rebaudiana genotypes. BMC genomics, 15(1), 571.
[9] Xia, E.-H., Yao, Q.-Y., Zhang, H.-B., Jiang, J.-J., Zhang, L.-P., & Gao, L.-Z. (2016). CandiSSR: an efficient pipeline used for identifying candidate polymorphic SSRs based on multiple assembled sequences. Frontiers in plant science, 6, 1171.
[10] Mathur, S., Bulchandani, N., Parihar, S., & Shekhawat, G. S. (2017). Critical Review on Steviol Glycosides: Pharmacological, Toxicological and Therapeutic Aspects of High Potency Zero Caloric Sweetener. International Journal of Pharmacology, 13(7), 916-928.
[11] Vieira, M. L. C., Santini, L., Diniz, A. L., & Munhoz, C. d. F. (2016). Microsatellite markers: what they mean and why they are so useful. Genetics and molecular biology, 39(3), 312-328.
[12] Sharma, N., Kaur, R., & Era, V. (2016). Potential of RAPD and ISSR markers for assessing genetic diversity among Stevia rebaudiana Bertoni accessions.
[13] Kaur, R., Sharma, N., & Raina, R. (2015). Identification and functional annotation of expressed sequencetags based SSR markers of Stevia rebaudiana. Turkish Journal of Agriculture and Forestry, 39(3), 439-450.
[14] Slater, G. S. C., & Birney, E. (2005). Automated generation of heuristics for biological sequence comparison. BMC bioinformatics, 6(1), 31.
[15] Shi, J., Huang, S., Zhan, J., Yu, J., Wang, X., Hua, W. & Wang, H. (2013). Genome-wide microsatellite characterization and marker development in the sequenced Brassica crop species. DNA research, 21(1), 53-68.
[16] Kim, M. J., Jin, J., Zheng, J., Wong, L., Chua, N.-H., & Jang, I.-C. (2015). Comparative transcriptomics unravel biochemical specialization of leaf tissues of Stevia (Stevia rebaudiana) for diterpenoid production. Plant physiology, pp. 01353.02015.
[17] Rodrıguez-Concepción, M., & Boronat, A. (2002). Elucidation of the methylerythritol phosphate pathway for isoprenoid biosynthesis in bacteria and plastids. A metabolic milestone achieved through genomics. Plant physiology, 130(3), 1079-1089.