Transcriptome Sequencing of Two Kentucky Bluegrass (Poa pratensis L.) Genotypes in Response to Heat Stress


  • Qiong LI Shanghai University, School of Life Science, Shanghai 200444 (CN)
  • Yali HE Shanghai Jiaotong University, School of Agriculture and Biology, Shanghai 200240 (CN)
  • Mingyue TU Shanghai University, School of Life Science, Shanghai 200444 (CN)
  • Junhui YAN Shanghai University, School of Life Science, Shanghai 200444 (CN)
  • Liangliang YU Shanghai University, School of Life Science, Shanghai 200444 (CN)
  • Weiwei QI Shanghai University, School of Life Science, Shanghai 200444 (CN)
  • Xiaojun YUAN Shanghai University, School of Life Science, Shanghai 200444 (CN)



heat stress; transcriptome sequencing; turfgrass; unigenes


Kentucky bluegrass (Poa pratensis L.) (KBG) is a major cool-season turfgrass. However, as its complex genetic background and production modes, limited genomic and transcriptomic information of KBG was known so far. In this study, a transcriptome study between wild type material Ninglan (summer stress sensitive) and cultivar material KBG03 (summer stress tolerant) was conducted, under optimal (25 °C) and high (40 °C) temperatures. A total of 81.42 Gb clean reads were generated and de novo assembled into 110,784 unigenes with an average length of 1,105 bp. About 50% KBG unigenes were similar to the non-redundant (NR) database. Among the NR BLASTx top hits, 27.47% unigenes were matched to Brachypodium distachyon. After heat stress, a massive amount of unigenes showed significantly differential expression in both genotypes. After 2h heat stress, more specially up-regulated differentially expressed unigenes (DEGs) and less down-regulated DEGs were detected in Ninglan than in KBG03. At 24h, the expression of 1710 and 730 unigenes were increased and decreased uniquely in Ninglan, and 1361 up-regulated DEGs and 757 down-regulated DEGs were just found in KBG03. More heat shock proteins (HSPs) and heat shock transcription factors (HSFs) DEGs were also identified at 2h than 24h in both genotypes. In addition, by Gene Ontology (GO) enrichment analysis, three core terms (“protein folding”, “response to heat”, and “response to hydrogen peroxide”) of biological process (BP) ontology were found in both genotypes under different heat stress time. The DEGs shared in both genotypes might be related to the basic mechanism of thermal response in KBG.


Ahn JH, Kim JS, Kim S, Soh HY, Shin H, Jang H, … Huh JH (2015). De novo transcriptome analysis to identify anthocyanin biosynthesis genes responsible for tissue-specific pigmentation in zoysiagrass (Zoysia japonica Steud.). PLoS One 10(4):e0124497.

Albertini E, Marconi G, Barcaccia G, Raggi L, Falcinelli M (2004). Isolation of candidate genes for apomixis in Poa pratensis L. Plant Molecular Biology 56(6):879-894.

Amaradasa BS, Amundsen K (2016). Transcriptome profiling of buffalograss challenged with the leaf spot pathogen Curvularia inaequalis. Frontiers in Plant Science 7:715.

An J, Shen XF, Ma QB, Yang CY, Liu SM, Chen Y (2014). Transcriptome profiling to discover putative genes associated with paraquat resistance in goosegrass (Eleusine indica L.). PLoS One 9(6):e99940.

Bushman BS, Amundsen KL, Warnke SE, Robins JG, Johnson PG (2016). Transcriptome profiling of Kentucky bluegrass (Poa pratensis L.) accessions in response to salt stress. BMC Genomics 17:48.

Chen SY, Huang X, Yan XQ, Liang Y, Wang YZ, Li XF, … Liu GS (2013). Transcriptome analysis in sheepgrass (Leymus chinensis): a dominant perennial grass of the Eurasian Steppe. PLoS One 8(7):e67974.

Czaban A, Sharma S, Byrne SL, Spannagl M, Mayer KFX, Asp T (2015). Comparative transcriptome analysis within the Lolium/Festuca species complex reveals high sequence conservation. BMC Genomics 16:249.

Dinkins RD, Nagabhyru P, Graham MA, Boykin D, Schardl CL (2017). Transcriptome response of Lolium arundinaceum to its fungal endophyte Epichloe coenophiala. New Phytologist 213(1):324-337.

Driedonks N, Xu JM, Peters JL, Park S, Rieu I (2015). Multi-level interactions between heat shock factors, heat shock proteins, and the redox system regulate acclimation to heat. Frontiers in Plant Science 6:999.

Farrell JD, Byrne S, Paina C, Asp T (2014). De novo assembly of the perennial ryegrass transcriptome using an RNA-Seq strategy. PLoS One 9(8):e103567.

Fox SE, Preece J, Kimbrel JA, Marchini GL, Sage A, Youens-Clark K, . . . Jaiswal P (2013). Sequencing and de novo transcriptome assembly of Brachypodium sylvaticum (Poaceae). Applications in Plant Science 1(3):1200011.

Gan L, Di R, Chao YH, Han LB, Chen XW, Wu C, Yin SX (2016). De novo transcriptome analysis for Kentucky bluegrass bwarf mutants induced by space mutation. PLoS One 11(5):e0151768.

Giorno F, Wolters-Arts M, Grillo S, Scharf KD, Vriezen WH, Mariani C (2010). Developmental and heat stress-regulated expression of HsfA2 and small heat shock proteins in tomato anthers. Journal of Experimental Botany 61(2):453-462.

Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, … Regev A (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nature Biotechnology 29(7):644-652.

Hahn A, Bublak D, Schleiff E, Scharf KD (2011). Crosstalk between Hsp90 and Hsp70 chaperones and heat stress transcription factors in tomato. Plant Cell 23(2):741-755.

Howarth C (1991). Molecular responses of plants to an increased incidence of heat shock. Plant, Cell & Environment 14(8):831-841.

Howarth C (2005). Genetic improvements of tolerance to high temperature. In: Ashraf M, Harris PJC (Eds). Abiotic stresses – plant resistance through breeding and molecular approaches. The Haworth Press, New York, pp 277-300.

Hu T, Sun XY, Zhang XZ, Nevo E, Fu JM (2014). An RNA sequencing transcriptome analysis of the high-temperature stressed tall fescue reveals novel insights into plant thermotolerance. BMC Genomics 15:1147.

Huang LK, Yan HD, Zhao XX, Zhang XQ, Wang J, Frazier T, … Liu W (2015). Identifying differentially expressed genes under heat stress and developing molecular markers in orchardgrass (Dactylis glomerata L.) through transcriptome analysis. Molecular Ecology Resources 15(6):1497-1509.

Huang X, Yan HD, Zhang XQ, Zhang J, Frazier TP, Huang DJ, . . . Yan YH (2016). De novo transcriptome analysis and molecular marker development of two Hemarthria species. Frontiers in Plant Science 7:496.

Huff DR (2003). Kentucky bluegrass. Turfgrass biology, genetics, and breeding. John Wiley & Sons Hoboken, NJ, pp 27-38.

Kovi MR, Abdelhalim M, Kunapareddy A, Ergon A, Tronsmo AM, Brurberg MB, … Rognli OA (2016). Global transcriptome changes in perennial ryegrass during early infection by pink snow mould. Scientific Reports 6:28702.

Li SF, Li F, Wang JW, Zhang W, Meng QW, Chen THH, … Yang XH (2011). Glycinebetaine enhances the tolerance of tomato plants to high temperature during germination of seeds and growth of seedlings. Plant Cell and Environment 34(11):1931-1943.

Li YF, Wang YX, Tang YH, Kakani VG, Mahalingam R (2013). Transcriptome analysis of heat stress response in switchgrass (Panicum virgatum L.). BMC Plant Biology 13:153.

Meyer WA (1982). Breeding disease-resistant cool-season turfgrass cultivars for the United States. Plant Disease 66(4):341-344.

Mishra SK, Tripp J, Winkelhaus S, Tschiersch B, Theres K, Nover L, Scharf KD (2002). In the complex family of heat stress transcription factors, HSfA1 has a unique role as master regulator of thermotolerance in tomato. Genes & Development 16(12):1555-1567.

Ni Y, Guo N, Zhao QL, Guo YJ (2016). Identification of candidate genes involved in wax deposition in Poa pratensis by RNA-seq. BMC Genomics 17:314.

Nover L, Bharti K, Doring P, Mishra SK, Ganguli A, Scharf KD (2001). Arabidopsis and the heat stress transcription factor world: how many heat stress transcription factors do we need? Cell Stress Chaperon 6(3):177-189.

Orshinsky AM, Hu JN, Opiyo SO, Reddyvari-Channarayappa V, Mitchell TK, Boehm MJ (2012). RNA-Seq analysis of the Sclerotinia homoeocarpa - creeping bentgrass pathosystem. PLoS One 7(8):e41150.

Priest HD, Fox SE, Rowley ER, Murray JR, Michael TP, Mockler TC (2014). Analysis of global gene expression in Brachypodium distachyon reveals extensive network plasticity in response to abiotic stress. PLoS One 9(1):e87499.

Puyang XH, An MY, Han LB, Zhang XZ (2015). Protective effect of spermidine on salt stress induced oxidative damage in two Kentucky bluegrass (Poa pratensis L.) cultivars. Ecotoxicology and Environmental Safety 117:96-106.

Sun J, Ren LP, Cheng Y, Gao JJ, Dong B, Chen SM, Chen FD, Jiang JF (2014). Identification of differentially expressed genes in Chrysanthemum nankingense (Asteraceae) under heat stress by RNA Seq. Gene 552(1):59-66.

Sung DY, Kaplan F, Lee KJ, Guy CL (2003). Acquired tolerance to temperature extremes. Trends in Plant Science 8(4):179-187.

Tian J, Belanger FC, Huang BR (2009). Identification of heat stress-responsive genes in heat-adapted thermal Agrostis scabra by suppression subtractive hybridization. Journal of Plant Physiology 166(6):588-601.

Toledo-Silva G, Cardoso-Silva CB, Jank L, Souza AP (2013). De novo transcriptome assembly for the tropical grass Panicum maximum Jacq. PLoS One 8(7):e70781.

Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, … Pachter L (2010). Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nature Biotechnology 28(5):511-515.

Vierling E (1991). The roles of heat shock proteins in plants. Annual review of plant biology 42(1):579-620.

Wachholtz M, Heng-Moss T, Twigg P, Baird L, Lu GQ, Amundsen K (2013). Transcriptome analysis of two buffalograss cultivars. BMC Genomics 14:613.

Wahid A, Gelani S, Ashraf M, Foolad MR (2007). Heat tolerance in plants: An overview. Environmental and Experimental Botany 61(3):199-223.

Wan XL, Zhou Q, Wang YY, Wang WE, Bao MZ, Zhang JW (2015). Identification of heat-responsive genes in carnation (Dianthus caryophyllus L.) by RNA-seq. Frontiers in Plant Science 6:519.

Wang DF, Luthe DS (2003). Heat sensitivity in a bentgrass variant. Failure to accumulate a chloroplast heat shock protein isoform implicated in heat tolerance. Plant Physiology 133(1):319-327.

Wei SJ, Du ZL, Gao F, Ke X, Li J, Liu JX, Zhou YJ (2015). Global transcriptome profiles of 'Meyer' Zoysiagrass in response to cold stress. PLoS One 10(6):e0131153.

Weis E, Berry JA (1988). Plants and high temperature stress. Symposia of the Society for Experimental Biology 42:329-346.

Xu KJ, Sun FL, Chai GQ, Wang YF, Shi LL, Liu SD, Xi YJ (2015). De novo assembly and transcriptome analysis of two contrary tillering mutants to learn the mechanisms of tillers outgrowth in switchgrass (Panicum virgatum L.). Frontiers in Plant Science 6:749.

Xu Y, Zhan C, Huang B (2011). Heat shock proteins in association with heat tolerance in grasses. International Journal of Proteomics 2011:529648.

Yang X, Yu XY, Li YF (2013). De novo assembly and characterization of the barnyardgrass (Echinochloa crus-galli) transcriptome using next-generation pyrosequencing. PLoS One 8(7):e69168.

Yokotani N, Ichikawa T, Kondou Y, Matsui M, Hirochika H, Iwabuchi M, Oda K (2008). Expression of rice heat stress transcription factor OsHsfA2e enhances tolerance to environmental stresses in transgenic Arabidopsis. Planta 227(5):957-967.

Zhang L, Ren J, Li T, Wang A, Tan D (2016). De novo transcriptome sequencing of cold-treated Kentucky bluegrass (Poa pratensis) and analysis of the genes involved in cold tolerance. Journal of Horticulture 23:182.

Zhang Y, Mian MAR, Chekhovskiy K, So S, Kupfer D, Lai H, Roe BA (2005). Differential gene expression in Festuca under heat stress conditions. Journal of Experimental Botany 56(413):897-907.




How to Cite

LI, Q., HE, Y., TU, M., YAN, J., YU, L., QI, W., & YUAN, X. (2018). Transcriptome Sequencing of Two Kentucky Bluegrass (Poa pratensis L.) Genotypes in Response to Heat Stress. Notulae Botanicae Horti Agrobotanici Cluj-Napoca, 47(2), 328–338.



Research Articles
DOI: 10.15835/nbha47111365