Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Deep sequencing-based characterization of transcriptome of trifoliate orange (<italic>Poncirus trifoliata</italic>
(L.) Raf.) in response to cold stress</title>
<author><name sortKey="Wang, Min" sort="Wang, Min" uniqKey="Wang M" first="Min" last="Wang">Min Wang</name>
<affiliation><nlm:aff id="Aff1">Key Laboratory of Horticultural Plant Biology (MOE), College of Horticulture and Forestry Science, Huazhong Agricultural University, 430070 Wuhan, China</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Zhang, Xiaona" sort="Zhang, Xiaona" uniqKey="Zhang X" first="Xiaona" last="Zhang">Xiaona Zhang</name>
<affiliation><nlm:aff id="Aff1">Key Laboratory of Horticultural Plant Biology (MOE), College of Horticulture and Forestry Science, Huazhong Agricultural University, 430070 Wuhan, China</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Liu, Ji Hong" sort="Liu, Ji Hong" uniqKey="Liu J" first="Ji-Hong" last="Liu">Ji-Hong Liu</name>
<affiliation><nlm:aff id="Aff1">Key Laboratory of Horticultural Plant Biology (MOE), College of Horticulture and Forestry Science, Huazhong Agricultural University, 430070 Wuhan, China</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">26219960</idno>
<idno type="pmc">4518522</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4518522</idno>
<idno type="RBID">PMC:4518522</idno>
<idno type="doi">10.1186/s12864-015-1629-7</idno>
<date when="2015">2015</date>
<idno type="wicri:Area/Pmc/Corpus">000588</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Deep sequencing-based characterization of transcriptome of trifoliate orange (<italic>Poncirus trifoliata</italic>
(L.) Raf.) in response to cold stress</title>
<author><name sortKey="Wang, Min" sort="Wang, Min" uniqKey="Wang M" first="Min" last="Wang">Min Wang</name>
<affiliation><nlm:aff id="Aff1">Key Laboratory of Horticultural Plant Biology (MOE), College of Horticulture and Forestry Science, Huazhong Agricultural University, 430070 Wuhan, China</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Zhang, Xiaona" sort="Zhang, Xiaona" uniqKey="Zhang X" first="Xiaona" last="Zhang">Xiaona Zhang</name>
<affiliation><nlm:aff id="Aff1">Key Laboratory of Horticultural Plant Biology (MOE), College of Horticulture and Forestry Science, Huazhong Agricultural University, 430070 Wuhan, China</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Liu, Ji Hong" sort="Liu, Ji Hong" uniqKey="Liu J" first="Ji-Hong" last="Liu">Ji-Hong Liu</name>
<affiliation><nlm:aff id="Aff1">Key Laboratory of Horticultural Plant Biology (MOE), College of Horticulture and Forestry Science, Huazhong Agricultural University, 430070 Wuhan, China</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">BMC Genomics</title>
<idno type="eISSN">1471-2164</idno>
<imprint><date when="2015">2015</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><sec><title>Background</title>
<p>Trifoliate orange (<italic>Poncirus trifoliata</italic>
(L.) Raf.) is extremely cold hardy after a full acclimation; however the underlying molecular mechanisms underlying this economically valuable trait remain poorly understood. In this study, global transcriptome profiles of trifoliate orange under cold conditions (4 °C) over a time course were generated by high-throughput sequencing.</p>
</sec>
<sec><title>Results</title>
<p>More than 68 million high-quality reads were produced and assembled into a non-redundant data of 77,292 unigenes with an average length of 1112 bp (N50 = 1778 bp). Of these, 23,846 had significant sequence similarity to known genes and these were assigned to 61 gene ontology (GO) categories and 25 clusters of orthologous groups (COG) involved in 128 KEGG pathways. Sequences derived from cold-treated and control plants were mapped to the assembled transcriptome, resulting in the identification of 5549 differentially expressed genes (DEGs). These comprised 600 (462 up-regulated, 138 down-regulated), 2346 (1631 up-regulated, 715 down-regulated), and 5177 (2702 up-regulated, 2475 down-regulated) genes from the cold-treated samples at 6, 24 and 72 h, respectively. The accuracy of the RNA-seq derived transcript expression data was validated by analyzing the expression patterns of 17 DEGs by qPCR. Plant hormone signal transduction, plant-pathogen interaction, and secondary metabolism were the most significantly enriched GO categories amongst in the DEGs. A total of 60 transcription factors were shown to be cold responsive. In addition, a number of genes involved in the catabolism and signaling of hormones, such as abscisic acid, ethylene and gibberellin, were affected by the cold stress. Meanwhile, levels of putrescine progressively increased under cold, which was consistent with up-regulation of an arginine decarboxylase gene.</p>
</sec>
<sec><title>Conclusions</title>
<p>This dataset provides valuable information regarding the trifoliate orange transcriptome changes in response to cold stress and may help guide future identification and functional analysis of genes that are importnatn for enhancing cold hardiness.</p>
</sec>
<sec><title>Electronic supplementary material</title>
<p>The online version of this article (doi:10.1186/s12864-015-1629-7) contains supplementary material, which is available to authorized users.</p>
</sec>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Chinnusamy, V" uniqKey="Chinnusamy V">V Chinnusamy</name>
</author>
<author><name sortKey="Zhu, Jh" uniqKey="Zhu J">JH Zhu</name>
</author>
<author><name sortKey="Zhu, Jk" uniqKey="Zhu J">JK Zhu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Michael, Mf" uniqKey="Michael M">MF Michael</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zhu, Jh" uniqKey="Zhu J">JH Zhu</name>
</author>
<author><name sortKey="Dong, Ch" uniqKey="Dong C">CH Dong</name>
</author>
<author><name sortKey="Zhu, Jk" uniqKey="Zhu J">JK Zhu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Medina, J" uniqKey="Medina J">J Medina</name>
</author>
<author><name sortKey="Catala, R" uniqKey="Catala R">R Catalá</name>
</author>
<author><name sortKey="Salinas, J" uniqKey="Salinas J">J Salinas</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nakashima, K" uniqKey="Nakashima K">K Nakashima</name>
</author>
<author><name sortKey="Ito, Y" uniqKey="Ito Y">Y Ito</name>
</author>
<author><name sortKey="Yamaguchi Shinozaki, K" uniqKey="Yamaguchi Shinozaki K">K Yamaguchi-Shinozaki</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chinnusamy, V" uniqKey="Chinnusamy V">V Chinnusamy</name>
</author>
<author><name sortKey="Ohta, M" uniqKey="Ohta M">M Ohta</name>
</author>
<author><name sortKey="Kanrar, S" uniqKey="Kanrar S">S Kanrar</name>
</author>
<author><name sortKey="Lee, B" uniqKey="Lee B">B Lee</name>
</author>
<author><name sortKey="Hong, Xh" uniqKey="Hong X">XH Hong</name>
</author>
<author><name sortKey="Agarwal, M" uniqKey="Agarwal M">M Agarwal</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Dong, Ch" uniqKey="Dong C">CH Dong</name>
</author>
<author><name sortKey="Agarwal, M" uniqKey="Agarwal M">M Agarwal</name>
</author>
<author><name sortKey="Zhang, Yy" uniqKey="Zhang Y">YY Zhang</name>
</author>
<author><name sortKey="Xie, Q" uniqKey="Xie Q">Q Xie</name>
</author>
<author><name sortKey="Zhu, Jk" uniqKey="Zhu J">JK Zhu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Agarwal, M" uniqKey="Agarwal M">M Agarwal</name>
</author>
<author><name sortKey="Hao, Yj" uniqKey="Hao Y">YJ Hao</name>
</author>
<author><name sortKey="Kapoor, A" uniqKey="Kapoor A">A Kapoor</name>
</author>
<author><name sortKey="Dong, Ch" uniqKey="Dong C">CH Dong</name>
</author>
<author><name sortKey="Fujii, H" uniqKey="Fujii H">H Fujii</name>
</author>
<author><name sortKey="Zheng, Xy" uniqKey="Zheng X">XY Zheng</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Matsui, A" uniqKey="Matsui A">A Matsui</name>
</author>
<author><name sortKey="Ishida, J" uniqKey="Ishida J">J Ishida</name>
</author>
<author><name sortKey="Morosawa, T" uniqKey="Morosawa T">T Morosawa</name>
</author>
<author><name sortKey="Mochizuki, Y" uniqKey="Mochizuki Y">Y Mochizuki</name>
</author>
<author><name sortKey="Kaminuma, E" uniqKey="Kaminuma E">E Kaminuma</name>
</author>
<author><name sortKey="Endo, Ta" uniqKey="Endo T">TA Endo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Winfield, Mo" uniqKey="Winfield M">MO Winfield</name>
</author>
<author><name sortKey="Lu, C" uniqKey="Lu C">C Lu</name>
</author>
<author><name sortKey="Wilson, Id" uniqKey="Wilson I">ID Wilson</name>
</author>
<author><name sortKey="Coghill, Ja" uniqKey="Coghill J">JA Coghill</name>
</author>
<author><name sortKey="Edwards, Kj" uniqKey="Edwards K">KJ Edwards</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tian, Dq" uniqKey="Tian D">DQ Tian</name>
</author>
<author><name sortKey="Pan, Xy" uniqKey="Pan X">XY Pan</name>
</author>
<author><name sortKey="Yu, Ym" uniqKey="Yu Y">YM Yu</name>
</author>
<author><name sortKey="Wang, Wy" uniqKey="Wang W">WY Wang</name>
</author>
<author><name sortKey="Zhang, F" uniqKey="Zhang F">F Zhang</name>
</author>
<author><name sortKey="Ge, Yy" uniqKey="Ge Y">YY Ge</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wang, Xc" uniqKey="Wang X">XC Wang</name>
</author>
<author><name sortKey="Zhao, Qy" uniqKey="Zhao Q">QY Zhao</name>
</author>
<author><name sortKey="Ma, Cl" uniqKey="Ma C">CL Ma</name>
</author>
<author><name sortKey="Zhang, Zh" uniqKey="Zhang Z">ZH Zhang</name>
</author>
<author><name sortKey="Cao, Hl" uniqKey="Cao H">HL Cao</name>
</author>
<author><name sortKey="Kong, Ym" uniqKey="Kong Y">YM Kong</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wang, Jm" uniqKey="Wang J">JM Wang</name>
</author>
<author><name sortKey="Yang, Y" uniqKey="Yang Y">Y Yang</name>
</author>
<author><name sortKey="Liu, Xh" uniqKey="Liu X">XH Liu</name>
</author>
<author><name sortKey="Huang, J" uniqKey="Huang J">J Huang</name>
</author>
<author><name sortKey="Wang, Q" uniqKey="Wang Q">Q Wang</name>
</author>
<author><name sortKey="Gu, Jh" uniqKey="Gu J">JH Gu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Fowler, S" uniqKey="Fowler S">S Fowler</name>
</author>
<author><name sortKey="Thomashow, Mf" uniqKey="Thomashow M">MF Thomashow</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Martin, Ja" uniqKey="Martin J">JA Martin</name>
</author>
<author><name sortKey="Wang, Z" uniqKey="Wang Z">Z Wang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Huang, Jz" uniqKey="Huang J">JZ Huang</name>
</author>
<author><name sortKey="Lu, X" uniqKey="Lu X">X Lu</name>
</author>
<author><name sortKey="Yan, H" uniqKey="Yan H">H Yan</name>
</author>
<author><name sortKey="Chen, Sy" uniqKey="Chen S">SY Chen</name>
</author>
<author><name sortKey="Zhang, Wk" uniqKey="Zhang W">WK Zhang</name>
</author>
<author><name sortKey="Huang, Rf" uniqKey="Huang R">RF Huang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Shi, Y" uniqKey="Shi Y">Y Shi</name>
</author>
<author><name sortKey="Yan, X" uniqKey="Yan X">X Yan</name>
</author>
<author><name sortKey="Zhao, Ps" uniqKey="Zhao P">PS Zhao</name>
</author>
<author><name sortKey="Yin, Hx" uniqKey="Yin H">HX Yin</name>
</author>
<author><name sortKey="Zhao, X" uniqKey="Zhao X">X Zhao</name>
</author>
<author><name sortKey="Xiao, Hl" uniqKey="Xiao H">HL Xiao</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chen, J" uniqKey="Chen J">J Chen</name>
</author>
<author><name sortKey="Tian, Q" uniqKey="Tian Q">Q Tian</name>
</author>
<author><name sortKey="Pang, T" uniqKey="Pang T">T Pang</name>
</author>
<author><name sortKey="Jiang, Lb" uniqKey="Jiang L">LB Jiang</name>
</author>
<author><name sortKey="Wu, Rl" uniqKey="Wu R">RL Wu</name>
</author>
<author><name sortKey="Xia, Xl" uniqKey="Xia X">XL Xia</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kakumanu, A" uniqKey="Kakumanu A">A Kakumanu</name>
</author>
<author><name sortKey="Ambavaram, Mmr" uniqKey="Ambavaram M">MMR Ambavaram</name>
</author>
<author><name sortKey="Klumas, C" uniqKey="Klumas C">C Klumas</name>
</author>
<author><name sortKey="Krishnan, A" uniqKey="Krishnan A">A Krishnan</name>
</author>
<author><name sortKey="Batlang, U" uniqKey="Batlang U">U Batlang</name>
</author>
<author><name sortKey="Myers, E" uniqKey="Myers E">E Myers</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wu, Ga" uniqKey="Wu G">GA Wu</name>
</author>
<author><name sortKey="Prochnik, S" uniqKey="Prochnik S">S Prochnik</name>
</author>
<author><name sortKey="Jenkins, J" uniqKey="Jenkins J">J Jenkins</name>
</author>
<author><name sortKey="Salse, J" uniqKey="Salse J">J Salse</name>
</author>
<author><name sortKey="Hellsten, U" uniqKey="Hellsten U">U Hellsten</name>
</author>
<author><name sortKey="Murat, F" uniqKey="Murat F">F Murat</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zhang, Ck" uniqKey="Zhang C">CK Zhang</name>
</author>
<author><name sortKey="Lang, P" uniqKey="Lang P">P Lang</name>
</author>
<author><name sortKey="Dane, F" uniqKey="Dane F">F Dane</name>
</author>
<author><name sortKey="Ebel, Rc" uniqKey="Ebel R">RC Ebel</name>
</author>
<author><name sortKey="Singh, Nk" uniqKey="Singh N">NK Singh</name>
</author>
<author><name sortKey="Locy, Rd" uniqKey="Locy R">RD Locy</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Peng, T" uniqKey="Peng T">T Peng</name>
</author>
<author><name sortKey="Zhu, Xf" uniqKey="Zhu X">XF Zhu</name>
</author>
<author><name sortKey="Fan, Qj" uniqKey="Fan Q">QJ Fan</name>
</author>
<author><name sortKey="Sun, Pp" uniqKey="Sun P">PP Sun</name>
</author>
<author><name sortKey="Liu, Jh" uniqKey="Liu J">JH Liu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sahin Cevik, M" uniqKey="Sahin Cevik M">M Sahin-Çevik</name>
</author>
<author><name sortKey="Moore, Ga" uniqKey="Moore G">GA Moore</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wang, J" uniqKey="Wang J">J Wang</name>
</author>
<author><name sortKey="Sun, Pp" uniqKey="Sun P">PP Sun</name>
</author>
<author><name sortKey="Chen, Cl" uniqKey="Chen C">CL Chen</name>
</author>
<author><name sortKey="Wang, Y" uniqKey="Wang Y">Y Wang</name>
</author>
<author><name sortKey="Fu, Xz" uniqKey="Fu X">XZ Fu</name>
</author>
<author><name sortKey="Liu, Jh" uniqKey="Liu J">JH Liu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Huang, Xs" uniqKey="Huang X">XS Huang</name>
</author>
<author><name sortKey="Wang, W" uniqKey="Wang W">W Wang</name>
</author>
<author><name sortKey="Zhang, Q" uniqKey="Zhang Q">Q Zhang</name>
</author>
<author><name sortKey="Liu, Jh" uniqKey="Liu J">JH Liu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Grabherr, Mg" uniqKey="Grabherr M">MG Grabherr</name>
</author>
<author><name sortKey="Haas, Bj" uniqKey="Haas B">BJ Haas</name>
</author>
<author><name sortKey="Yassour, M" uniqKey="Yassour M">M Yassour</name>
</author>
<author><name sortKey="Levin, Jz" uniqKey="Levin J">JZ Levin</name>
</author>
<author><name sortKey="Thompson, Da" uniqKey="Thompson D">DA Thompson</name>
</author>
<author><name sortKey="Amit, I" uniqKey="Amit I">I Amit</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ye, J" uniqKey="Ye J">J Ye</name>
</author>
<author><name sortKey="Fang, L" uniqKey="Fang L">L Fang</name>
</author>
<author><name sortKey="Zheng, Hk" uniqKey="Zheng H">HK Zheng</name>
</author>
<author><name sortKey="Zhang, Y" uniqKey="Zhang Y">Y Zhang</name>
</author>
<author><name sortKey="Chen, J" uniqKey="Chen J">J Chen</name>
</author>
<author><name sortKey="Zhang, Zj" uniqKey="Zhang Z">ZJ Zhang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Conesa, A" uniqKey="Conesa A">A Conesa</name>
</author>
<author><name sortKey="Gotz, S" uniqKey="Gotz S">S Götz</name>
</author>
<author><name sortKey="Garcia G Mez, Jm" uniqKey="Garcia G Mez J">JM García-Gómez</name>
</author>
<author><name sortKey="Terol, J" uniqKey="Terol J">J Terol</name>
</author>
<author><name sortKey="Tal N, M" uniqKey="Tal N M">M Talón</name>
</author>
<author><name sortKey="Robles, M" uniqKey="Robles M">M Robles</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Iseli, C" uniqKey="Iseli C">C Iseli</name>
</author>
<author><name sortKey="Jongeneel, Cv" uniqKey="Jongeneel C">CV Jongeneel</name>
</author>
<author><name sortKey="Bucher, P" uniqKey="Bucher P">P Bucher</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Li, Rq" uniqKey="Li R">RQ Li</name>
</author>
<author><name sortKey="Yu, C" uniqKey="Yu C">C Yu</name>
</author>
<author><name sortKey="Li, Yr" uniqKey="Li Y">YR Li</name>
</author>
<author><name sortKey="Lam, Tw" uniqKey="Lam T">TW Lam</name>
</author>
<author><name sortKey="Yiu, Sm" uniqKey="Yiu S">SM Yiu</name>
</author>
<author><name sortKey="Kristianse, K" uniqKey="Kristianse K">K Kristianse</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kim, Yo" uniqKey="Kim Y">YO Kim</name>
</author>
<author><name sortKey="Kim, Js" uniqKey="Kim J">JS Kim</name>
</author>
<author><name sortKey="Kang, H" uniqKey="Kang H">H Kang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kim, Jy" uniqKey="Kim J">JY Kim</name>
</author>
<author><name sortKey="Park, Sj" uniqKey="Park S">SJ Park</name>
</author>
<author><name sortKey="Kwak, Kj" uniqKey="Kwak K">KJ Kwak</name>
</author>
<author><name sortKey="Jung, Ch" uniqKey="Jung C">CH Jung</name>
</author>
<author><name sortKey="Ahn, Sj" uniqKey="Ahn S">SJ Ahn</name>
</author>
<author><name sortKey="Goh, Ch" uniqKey="Goh C">CH Goh</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kovtun, Y" uniqKey="Kovtun Y">Y Kovtun</name>
</author>
<author><name sortKey="Chiu, Wl" uniqKey="Chiu W">WL Chiu</name>
</author>
<author><name sortKey="Tena, G" uniqKey="Tena G">G Tena</name>
</author>
<author><name sortKey="Sheen, J" uniqKey="Sheen J">J Sheen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jaspers, P" uniqKey="Jaspers P">P Jaspers</name>
</author>
<author><name sortKey="Kangasj Rvi, J" uniqKey="Kangasj Rvi J">J Kangasjärvi</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Noctor, G" uniqKey="Noctor G">G Noctor</name>
</author>
<author><name sortKey="Veljovic Jovanovic, S" uniqKey="Veljovic Jovanovic S">S Veljovic-Jovanovic</name>
</author>
<author><name sortKey="Driscoll, S" uniqKey="Driscoll S">S Driscoll</name>
</author>
<author><name sortKey="Novitskaya, L" uniqKey="Novitskaya L">L Novitskaya</name>
</author>
<author><name sortKey="Foyer, Ch" uniqKey="Foyer C">CH Foyer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Shinozaki, K" uniqKey="Shinozaki K">K Shinozaki</name>
</author>
<author><name sortKey="Yamaguchi Shinozaki, K" uniqKey="Yamaguchi Shinozaki K">K Yamaguchi-Shinozaki</name>
</author>
<author><name sortKey="Seki, M" uniqKey="Seki M">M Seki</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Luan, S" uniqKey="Luan S">S Luan</name>
</author>
<author><name sortKey="Kudla, J" uniqKey="Kudla J">J Kudla</name>
</author>
<author><name sortKey="Rodriguez Concepcion, M" uniqKey="Rodriguez Concepcion M">M Rodriguez-Concepcion</name>
</author>
<author><name sortKey="Yalovsky, S" uniqKey="Yalovsky S">S Yalovsky</name>
</author>
<author><name sortKey="Gruissem, W" uniqKey="Gruissem W">W Gruissem</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Valmonte, Gr" uniqKey="Valmonte G">GR Valmonte</name>
</author>
<author><name sortKey="Arthur, K" uniqKey="Arthur K">K Arthur</name>
</author>
<author><name sortKey="Higgins, Cm" uniqKey="Higgins C">CM Higgins</name>
</author>
<author><name sortKey="Macdiarmid, Rm" uniqKey="Macdiarmid R">RM MacDiarmid</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Franza, S" uniqKey="Franza S">S Franza</name>
</author>
<author><name sortKey="Ehlerta, B" uniqKey="Ehlerta B">B Ehlerta</name>
</author>
<author><name sortKey="Liesea, A" uniqKey="Liesea A">A Liesea</name>
</author>
<author><name sortKey="Kurtha, J" uniqKey="Kurtha J">J Kurtha</name>
</author>
<author><name sortKey="Cazale, Ac" uniqKey="Cazale A">AC Cazalé</name>
</author>
<author><name sortKey="Romeisa, T" uniqKey="Romeisa T">T Romeisa</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Xu, Yj" uniqKey="Xu Y">YJ Xu</name>
</author>
<author><name sortKey="Gao, S" uniqKey="Gao S">S Gao</name>
</author>
<author><name sortKey="Yang, Yj" uniqKey="Yang Y">YJ Yang</name>
</author>
<author><name sortKey="Huang, Ny" uniqKey="Huang N">NY Huang</name>
</author>
<author><name sortKey="Cheng, Ln" uniqKey="Cheng L">LN Cheng</name>
</author>
<author><name sortKey="Wei, Q" uniqKey="Wei Q">Q Wei</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Janska, A" uniqKey="Janska A">A Janská</name>
</author>
<author><name sortKey="Marsik, P" uniqKey="Marsik P">P Maršík</name>
</author>
<author><name sortKey="Zelenkova, S" uniqKey="Zelenkova S">S Zelenková</name>
</author>
<author><name sortKey="Ovesna, J" uniqKey="Ovesna J">J Ovesná</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Peleg, Z" uniqKey="Peleg Z">Z Peleg</name>
</author>
<author><name sortKey="Blumwald, E" uniqKey="Blumwald E">E Blumwald</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tuteja, N" uniqKey="Tuteja N">N Tuteja</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Heino, P" uniqKey="Heino P">P Heino</name>
</author>
<author><name sortKey="Sandman, G" uniqKey="Sandman G">G Sandman</name>
</author>
<author><name sortKey="L Ng, V" uniqKey="L Ng V">V Lång</name>
</author>
<author><name sortKey="Nordin, K" uniqKey="Nordin K">K Nordin</name>
</author>
<author><name sortKey="Palva, Et" uniqKey="Palva E">ET Palva</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lang, V" uniqKey="Lang V">V Lang</name>
</author>
<author><name sortKey="Mantyla, E" uniqKey="Mantyla E">E Mantyla</name>
</author>
<author><name sortKey="Welin, B" uniqKey="Welin B">B Welin</name>
</author>
<author><name sortKey="Palva, Et" uniqKey="Palva E">ET Palva</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Espasandin, Fd" uniqKey="Espasandin F">FD Espasandin</name>
</author>
<author><name sortKey="Maiale, Sj" uniqKey="Maiale S">SJ Maiale</name>
</author>
<author><name sortKey="Calzadilla, P" uniqKey="Calzadilla P">P Calzadilla</name>
</author>
<author><name sortKey="Ruiz, Oa" uniqKey="Ruiz O">OA Ruiz</name>
</author>
<author><name sortKey="Sansberro, Pa" uniqKey="Sansberro P">PA Sansberro</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="T Htiharju, S" uniqKey="T Htiharju S">S Tähtiharju</name>
</author>
<author><name sortKey="Palva, T" uniqKey="Palva T">T Palva</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cao, Wh" uniqKey="Cao W">WH Cao</name>
</author>
<author><name sortKey="Liu, J" uniqKey="Liu J">J Liu</name>
</author>
<author><name sortKey="He, Xj" uniqKey="He X">XJ He</name>
</author>
<author><name sortKey="Mu, Rl" uniqKey="Mu R">RL Mu</name>
</author>
<author><name sortKey="Zhou, Hl" uniqKey="Zhou H">HL Zhou</name>
</author>
<author><name sortKey="Chen, Sy" uniqKey="Chen S">SY Chen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jung, Jy" uniqKey="Jung J">JY Jung</name>
</author>
<author><name sortKey="Shin, R" uniqKey="Shin R">R Shin</name>
</author>
<author><name sortKey="Schachtman, Dp" uniqKey="Schachtman D">DP Schachtman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kendrick, Md" uniqKey="Kendrick M">MD Kendrick</name>
</author>
<author><name sortKey="Chang, C" uniqKey="Chang C">C Chang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chu, C" uniqKey="Chu C">C Chu</name>
</author>
<author><name sortKey="Lee, Tm" uniqKey="Lee T">TM Lee</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ciardi, Ja" uniqKey="Ciardi J">JA Ciardi</name>
</author>
<author><name sortKey="Deikman, J" uniqKey="Deikman J">J Deikman</name>
</author>
<author><name sortKey="Orzolek, Md" uniqKey="Orzolek M">MD Orzolek</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Yu, Xm" uniqKey="Yu X">XM Yu</name>
</author>
<author><name sortKey="Griffith, M" uniqKey="Griffith M">M Griffith</name>
</author>
<author><name sortKey="Wiseman, Sb" uniqKey="Wiseman S">SB Wiseman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zhang, Z" uniqKey="Zhang Z">Z Zhang</name>
</author>
<author><name sortKey="Huang, R" uniqKey="Huang R">R Huang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Catala, R" uniqKey="Catala R">R Catalá</name>
</author>
<author><name sortKey="L Pez Cobollo, R" uniqKey="L Pez Cobollo R">R López-Cobollo</name>
</author>
<author><name sortKey="Castellano, Mm" uniqKey="Castellano M">MM Castellano</name>
</author>
<author><name sortKey="Angostob, T" uniqKey="Angostob T">T Angostob</name>
</author>
<author><name sortKey="Alonsoc, Jm" uniqKey="Alonsoc J">JM Alonsoc</name>
</author>
<author><name sortKey="Eckerc, Jr" uniqKey="Eckerc J">JR Eckerc</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Shi, Y" uniqKey="Shi Y">Y Shi</name>
</author>
<author><name sortKey="Tian, S" uniqKey="Tian S">S Tian</name>
</author>
<author><name sortKey="Hou, L" uniqKey="Hou L">L Hou</name>
</author>
<author><name sortKey="Huang, Xz" uniqKey="Huang X">XZ Huang</name>
</author>
<author><name sortKey="Zhang, Xy" uniqKey="Zhang X">XY Zhang</name>
</author>
<author><name sortKey="Guo, Hw" uniqKey="Guo H">HW Guo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Colebrook, Eh" uniqKey="Colebrook E">EH Colebrook</name>
</author>
<author><name sortKey="Thomas, Sg" uniqKey="Thomas S">SG Thomas</name>
</author>
<author><name sortKey="Phillips, Al" uniqKey="Phillips A">AL Phillips</name>
</author>
<author><name sortKey="Hedden, P" uniqKey="Hedden P">P Hedden</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zentella, R" uniqKey="Zentella R">R Zentella</name>
</author>
<author><name sortKey="Zhang, Zl" uniqKey="Zhang Z">ZL Zhang</name>
</author>
<author><name sortKey="Park, M" uniqKey="Park M">M Park</name>
</author>
<author><name sortKey="Thomas, Sg" uniqKey="Thomas S">SG Thomas</name>
</author>
<author><name sortKey="Endo, A" uniqKey="Endo A">A Endo</name>
</author>
<author><name sortKey="Murase, K" uniqKey="Murase K">K Murase</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Achard, P" uniqKey="Achard P">P Achard</name>
</author>
<author><name sortKey="Gong, F" uniqKey="Gong F">F Gong</name>
</author>
<author><name sortKey="Cheminant, S" uniqKey="Cheminant S">S Cheminant</name>
</author>
<author><name sortKey="Alioua, M" uniqKey="Alioua M">M Alioua</name>
</author>
<author><name sortKey="Hedden, P" uniqKey="Hedden P">P Hedden</name>
</author>
<author><name sortKey="Genschik, P" uniqKey="Genschik P">P Genschik</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Miller, Gad" uniqKey="Miller G">GAD Miller</name>
</author>
<author><name sortKey="Mittler, R" uniqKey="Mittler R">R Mittler</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Suzuki, N" uniqKey="Suzuki N">N Suzuki</name>
</author>
<author><name sortKey="Bajad, S" uniqKey="Bajad S">S Bajad</name>
</author>
<author><name sortKey="Shuman, J" uniqKey="Shuman J">J Shuman</name>
</author>
<author><name sortKey="Shulaev, V" uniqKey="Shulaev V">V Shulaev</name>
</author>
<author><name sortKey="Mittler, R" uniqKey="Mittler R">R Mittler</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Agarwal, Pk" uniqKey="Agarwal P">PK Agarwal</name>
</author>
<author><name sortKey="Jha, B" uniqKey="Jha B">B Jha</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Liu, Jh" uniqKey="Liu J">JH Liu</name>
</author>
<author><name sortKey="Peng, T" uniqKey="Peng T">T Peng</name>
</author>
<author><name sortKey="Dai, W" uniqKey="Dai W">W Dai</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Zhang, Lc" uniqKey="Zhang L">LC Zhang</name>
</author>
<author><name sortKey="Liu, Gx" uniqKey="Liu G">GX Liu</name>
</author>
<author><name sortKey="Zhao, Gy" uniqKey="Zhao G">GY Zhao</name>
</author>
<author><name sortKey="Xia, C" uniqKey="Xia C">C Xia</name>
</author>
<author><name sortKey="Jia, Jz" uniqKey="Jia J">JZ Jia</name>
</author>
<author><name sortKey="Liu, X" uniqKey="Liu X">X Liu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tiburcio, Af" uniqKey="Tiburcio A">AF Tiburcio</name>
</author>
<author><name sortKey="Altabella, T" uniqKey="Altabella T">T Altabella</name>
</author>
<author><name sortKey="Bitrian, M" uniqKey="Bitrian M">M Bitrián</name>
</author>
<author><name sortKey="Alcazar, R" uniqKey="Alcazar R">R Alcázar</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pertea, G" uniqKey="Pertea G">G Pertea</name>
</author>
<author><name sortKey="Huang, X" uniqKey="Huang X">X Huang</name>
</author>
<author><name sortKey="Liang, F" uniqKey="Liang F">F Liang</name>
</author>
<author><name sortKey="Antonescu, V" uniqKey="Antonescu V">V Antonescu</name>
</author>
<author><name sortKey="Sultana, R" uniqKey="Sultana R">R Sultana</name>
</author>
<author><name sortKey="Karamycheva, S" uniqKey="Karamycheva S">S Karamycheva</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kanehisa, M" uniqKey="Kanehisa M">M Kanehisa</name>
</author>
<author><name sortKey="Goto, S" uniqKey="Goto S">S Goto</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mortazavi, A" uniqKey="Mortazavi A">A Mortazavi</name>
</author>
<author><name sortKey="Williams, Ba" uniqKey="Williams B">BA Williams</name>
</author>
<author><name sortKey="Mccue, K" uniqKey="Mccue K">K McCue</name>
</author>
<author><name sortKey="Schaeffer, L" uniqKey="Schaeffer L">L Schaeffer</name>
</author>
<author><name sortKey="Wold, B" uniqKey="Wold B">B Wold</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Fu, Xz" uniqKey="Fu X">XZ Fu</name>
</author>
<author><name sortKey="Chen, Cw" uniqKey="Chen C">CW Chen</name>
</author>
<author><name sortKey="Wang, Y" uniqKey="Wang Y">Y Wang</name>
</author>
<author><name sortKey="Liu, Jh" uniqKey="Liu J">JH Liu</name>
</author>
<author><name sortKey="Moriguchi, T" uniqKey="Moriguchi T">T Moriguchi</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">BMC Genomics</journal-id>
<journal-id journal-id-type="iso-abbrev">BMC Genomics</journal-id>
<journal-title-group><journal-title>BMC Genomics</journal-title>
</journal-title-group>
<issn pub-type="epub">1471-2164</issn>
<publisher><publisher-name>BioMed Central</publisher-name>
<publisher-loc>London</publisher-loc>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">26219960</article-id>
<article-id pub-id-type="pmc">4518522</article-id>
<article-id pub-id-type="publisher-id">1629</article-id>
<article-id pub-id-type="doi">10.1186/s12864-015-1629-7</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject>
</subj-group>
</article-categories>
<title-group><article-title>Deep sequencing-based characterization of transcriptome of trifoliate orange (<italic>Poncirus trifoliata</italic>
(L.) Raf.) in response to cold stress</article-title>
</title-group>
<contrib-group><contrib contrib-type="author"><name><surname>Wang</surname>
<given-names>Min</given-names>
</name>
<address><email>510178941@qq.com</email>
</address>
<xref ref-type="aff" rid="Aff1"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Zhang</surname>
<given-names>Xiaona</given-names>
</name>
<address><email>510178941@qq.com</email>
</address>
<xref ref-type="aff" rid="Aff1"></xref>
</contrib>
<contrib contrib-type="author" corresp="yes"><name><surname>Liu</surname>
<given-names>Ji-Hong</given-names>
</name>
<address><email>liujihong@mail.hzau.edu.cn</email>
</address>
<xref ref-type="aff" rid="Aff1"></xref>
</contrib>
<aff id="Aff1">Key Laboratory of Horticultural Plant Biology (MOE), College of Horticulture and Forestry Science, Huazhong Agricultural University, 430070 Wuhan, China</aff>
</contrib-group>
<pub-date pub-type="epub"><day>29</day>
<month>7</month>
<year>2015</year>
</pub-date>
<pub-date pub-type="pmc-release"><day>29</day>
<month>7</month>
<year>2015</year>
</pub-date>
<pub-date pub-type="collection"><year>2015</year>
</pub-date>
<volume>16</volume>
<issue>1</issue>
<elocation-id>555</elocation-id>
<history><date date-type="received"><day>12</day>
<month>2</month>
<year>2015</year>
</date>
<date date-type="accepted"><day>12</day>
<month>5</month>
<year>2015</year>
</date>
</history>
<permissions><copyright-statement>© Wang et al. 2015</copyright-statement>
<license license-type="open-access"><license-p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0">http://creativecommons.org/licenses/by/4.0</ext-link>
), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/publicdomain/zero/1.0/">http://creativecommons.org/publicdomain/zero/1.0/</ext-link>
) applies to the data made available in this article, unless otherwise stated.</license-p>
</license>
</permissions>
<abstract id="Abs1"><sec><title>Background</title>
<p>Trifoliate orange (<italic>Poncirus trifoliata</italic>
(L.) Raf.) is extremely cold hardy after a full acclimation; however the underlying molecular mechanisms underlying this economically valuable trait remain poorly understood. In this study, global transcriptome profiles of trifoliate orange under cold conditions (4 °C) over a time course were generated by high-throughput sequencing.</p>
</sec>
<sec><title>Results</title>
<p>More than 68 million high-quality reads were produced and assembled into a non-redundant data of 77,292 unigenes with an average length of 1112 bp (N50 = 1778 bp). Of these, 23,846 had significant sequence similarity to known genes and these were assigned to 61 gene ontology (GO) categories and 25 clusters of orthologous groups (COG) involved in 128 KEGG pathways. Sequences derived from cold-treated and control plants were mapped to the assembled transcriptome, resulting in the identification of 5549 differentially expressed genes (DEGs). These comprised 600 (462 up-regulated, 138 down-regulated), 2346 (1631 up-regulated, 715 down-regulated), and 5177 (2702 up-regulated, 2475 down-regulated) genes from the cold-treated samples at 6, 24 and 72 h, respectively. The accuracy of the RNA-seq derived transcript expression data was validated by analyzing the expression patterns of 17 DEGs by qPCR. Plant hormone signal transduction, plant-pathogen interaction, and secondary metabolism were the most significantly enriched GO categories amongst in the DEGs. A total of 60 transcription factors were shown to be cold responsive. In addition, a number of genes involved in the catabolism and signaling of hormones, such as abscisic acid, ethylene and gibberellin, were affected by the cold stress. Meanwhile, levels of putrescine progressively increased under cold, which was consistent with up-regulation of an arginine decarboxylase gene.</p>
</sec>
<sec><title>Conclusions</title>
<p>This dataset provides valuable information regarding the trifoliate orange transcriptome changes in response to cold stress and may help guide future identification and functional analysis of genes that are importnatn for enhancing cold hardiness.</p>
</sec>
<sec><title>Electronic supplementary material</title>
<p>The online version of this article (doi:10.1186/s12864-015-1629-7) contains supplementary material, which is available to authorized users.</p>
</sec>
</abstract>
<kwd-group xml:lang="en"><title>Keywords</title>
<kwd><italic>Poncirus trifoliata</italic>
</kwd>
<kwd>RNA-seq</kwd>
<kwd>Cold stress</kwd>
<kwd>Transcriptome profiling</kwd>
<kwd>Digital gene expression</kwd>
<kwd>Citrus</kwd>
</kwd-group>
<custom-meta-group><custom-meta><meta-name>issue-copyright-statement</meta-name>
<meta-value>© The Author(s) 2015</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body><sec id="Sec1" sec-type="introduction"><title>Background</title>
<p>Plants are frequently challenged by various environmental stresses, among which cold temperatures are one of the major factors that adversely influence plant growth and development, crop yield potential and geographic distribution [<xref ref-type="bibr" rid="CR1">1</xref>
]. Plants have evolved a diverse set of adaptive mechanisms in order to withstand cold stresses, one of these is acquired freezing tolerance following prior exposure to non-freezing low temperature, a phenomenon known as cold acclimation (CA) [<xref ref-type="bibr" rid="CR2">2</xref>
]. Studies of the mechanisms underlying CA-mediated improvement of freezing tolerance have suggested the importance of a wide range of physiological, biochemical, cellular and molecular processes, and these have been associated with modulation of gene transcription [<xref ref-type="bibr" rid="CR3">3</xref>
].</p>
<p>Significant progress has been made in identifying key components of the cold stress signaling network. Of particular note was the identification and characterization of C-repeat (CRT)-binding factors (CBFs), also known as the dehydration-responsive element-binding factors 1 (DREB1s) [<xref ref-type="bibr" rid="CR4">4</xref>
], which have been reported to be conserved among dicots and monocots. CBF genes play critical roles in mediating cold stress response by regulating a spectrum of cold-regulated (COR) genes, collectively called the CBF regulon, through binding to the <italic>cis</italic>
-acting element (CRT/DRE) within their promoters [<xref ref-type="bibr" rid="CR4">4</xref>
, <xref ref-type="bibr" rid="CR5">5</xref>
]. Regulation of COR genes by CBFs constitutes the predominant cold signaling pathway in plants. In addition, characterization of either positive or negative regulators of CBF genes, including <italic>ICE1</italic>
(Inducer of CBF Expression 1), <italic>HOS1</italic>
(High Expression of Osmotically Responsive Gene 1), and <italic>MYB15</italic>
has provided a more complete understanding of the complexity of CBF-mediated cold signaling [<xref ref-type="bibr" rid="CR6">6</xref>
–<xref ref-type="bibr" rid="CR8">8</xref>
].</p>
<p>The expression patterns of large numbers of cold-responsive genes, including both regulatory genes and structural genes, are altered during exposure to cold [<xref ref-type="bibr" rid="CR5">5</xref>
]. For example, in <italic>Arabidopsis thaliana</italic>
approximately 1000 genes were shown to exhibit different expression under cold stress, while more than 2 % of wheat genome were reported to show altered expression in response to cold [<xref ref-type="bibr" rid="CR9">9</xref>
, <xref ref-type="bibr" rid="CR10">10</xref>
]. Several other studies have reported similar changes in gene expression patterns in various plant species [<xref ref-type="bibr" rid="CR11">11</xref>
–<xref ref-type="bibr" rid="CR13">13</xref>
]. Although increasing numbers of genes involved in cold responses are being identified and characterized in an array of plant species, the proportions are still small relative to observed genome wide changes, and most genes involved in this process have yet to be functionally characterized. On the other hand, plants have also developed a CBF-independent pathway for adapting to cold stress, which is corroborated by the finding that only 12 % of cold-responsive genes are likely under the control of CBFs [<xref ref-type="bibr" rid="CR14">14</xref>
]. However, this CBF-independent pathway has not been well characterized.</p>
<p>Genome-wide transcriptome analysis represents a potentially valuable strategy for elucidating the breadth of molecular mechanisms underlying physiological processes and it can substantially increase the efficiency of identifying genes of interest. Recent developments in high-throughput DNA sequencing technologies and associated analytical approaches, such as RNA-Seq, allow the generation of large-scale transcriptome data for both model and non-model species [<xref ref-type="bibr" rid="CR15">15</xref>
]. Over the last few years, deep sequencing using Illumina-based RNA-Seq has been increasingly applied to capture an overview of RNA transcript profiles in a range of plant species following exposure to diverse and adverse environmental conditions, including high salinity, drought, cold temperatures, and pathogen attack [<xref ref-type="bibr" rid="CR11">11</xref>
–<xref ref-type="bibr" rid="CR13">13</xref>
, <xref ref-type="bibr" rid="CR16">16</xref>
–<xref ref-type="bibr" rid="CR18">18</xref>
]. For example, Illumina sequencing-based snapshots of the transcriptomes of non-model species, such as <italic>Anthrurim andraeanum</italic>
, <italic>Camellia sinensis</italic>
and <italic>Lilium lancifolium</italic>
under cold treatment have been reported, together with identification of many differentially expressed genes (DEGs) [<xref ref-type="bibr" rid="CR11">11</xref>
–<xref ref-type="bibr" rid="CR13">13</xref>
]. RNA-Seq has thus been proven to act as a powerful approach for creating transcriptomic data for non-model plants, particularly those without substantial existing sequence data [<xref ref-type="bibr" rid="CR19">19</xref>
].</p>
<p>Trifoliate orange (<italic>Poncirus trifoliata</italic>
(L.) Raf.), a diploid (2n = 2× = 18) with high degree of heterozygosity, is most closely related to <italic>Citrus</italic>
spp., the genome of which has been sequenced [<xref ref-type="bibr" rid="CR20">20</xref>
]. One of its notable attributes is the striking cold hardiness after a period of CA. To date, several cold-responsive <italic>P. trifoliata</italic>
genes have been isolated by homology-based cloning and low-throughput cDNA library screening. One such study described the characterization of cDNAs corresponding to eight up-regulated and six down-regulated <italic>P. trifoliata</italic>
genes following a gradual CA temperature regime, using mRNA differential display-reverse transcription (DDRT)-PCR [<xref ref-type="bibr" rid="CR21">21</xref>
]. Additionally, suppression subtractive hybridization (SSH) has been used to screen cDNA libraries, resulting in the identification of non-redundant differentially expressed sequence tags (ESTs) [<xref ref-type="bibr" rid="CR22">22</xref>
, <xref ref-type="bibr" rid="CR23">23</xref>
]. Based on these studies, several <italic>P. trifoliata</italic>
stress-responsive regulatory or structural genes, such as <italic>PtADC</italic>
, and <italic>PtrbHLH</italic>
, have been functionally characterized and shown to play key roles in cold tolerance [<xref ref-type="bibr" rid="CR24">24</xref>
, <xref ref-type="bibr" rid="CR25">25</xref>
]. However, such analytical approaches have not provided a global overview of the molecular mechanisms underlying the CA-mediated cold tolerance in <italic>P. trifoliata</italic>
. In this current study, a comprehensive transcriptome analysis of trifoliate orange under cold stress was performed, using the Illumina-based paired-end deep sequencing platform. More than 6.1 million pairs of high-quality sequence were assembled and annotated, and large sets associated with transcripts involved in diverse metabolic processes and signaling pathways were identified. In addition, a total of 5549 predicted genes with regulatory or protective roles were revealed to be either up- or down-regulated by cold at different time points. Exposure to cold temperatures resulted in an increase in the abundance of transcripts associated with several metabolic pathways and the expression data further suggested the involvement of both the CBF-dependent and independent pathways in the cold responses. The RNA-seq and digital expression profiling described here have yielded valuable insights into the understanding the molecular events related to cold responses in trifoliate orange and the dataset represents a significant resources for targeting genes in future breeding efforts.</p>
</sec>
<sec id="Sec2" sec-type="results"><title>Results</title>
<sec id="Sec3"><title>Transcriptome sequencing and <italic>de novo</italic>
assembly</title>
<p>A cDNA library was constructed using equal amounts of RNA extracted from <italic>P. trifoliata</italic>
seedlings that had been exposed to various stresses, including cold temperatures (4 °C), high salinity and drought. To characterize the trifoliate orange transcriptome, the cDNA library was subjected to paired-end (PE) read sequencing using the Illumina HiSeq2000 platform. After removing reads of low quality, adaptor sequences or reads with >5 % ambiguous nucleotides, we got a total of 68,041,582 clean PE reads consisting of 6,123,742,380 nucleotides (5.84 Gb), with an average GC content of 45.2 %. These high-quality reads were then <italic>de novo</italic>
assembled using the Trinity program [<xref ref-type="bibr" rid="CR26">26</xref>
], resulting in 115,114 contigs with an average length of 372 bp and an N50 length of 776 bp (Table <xref rid="Tab1" ref-type="table">1</xref>
). Based on the PE sequence information, the contigs were further assembled to give 77,292 unigenes, which accounted for 85,969,539 nucleotides (82.0 Mb), with an average length of 1112 bp. Of these unigenes, 48,256 (62.5 %) were longer than 500 bp, 31,674 (41.0 %) were longer than 1000 bp, 20,209 (26.2 %) were longer than 1500 bp and 12,309 (15.9 %) were longer than 2000 bp (Fig. <xref rid="Fig1" ref-type="fig">1</xref>
). Further analysis showed that sequencing depth of the unigenes ranged from 0.0175 to 93,708, with an average of 49.57, while the number of reads uniquely mapped to a unigene ranged from 1 to 1,675,941. The transcriptome data have been deposited at DDBJ/EMBL/GenBank under the accession number of GCVN00000000.<table-wrap id="Tab1"><label>Table 1</label>
<caption><p>Summary of <italic>P. trifoliata</italic>
RNA-Seq data</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th></th>
<th>High-quality sequences (n)</th>
<th>Total bases (base pairs, bp)</th>
<th>Average length (bp)</th>
<th>N50 (bp)</th>
</tr>
</thead>
<tbody><tr><td align="justify">Reads</td>
<td>68,041,582</td>
<td>6,123,742,380</td>
<td>2*90</td>
<td>-</td>
</tr>
<tr><td align="justify">Contigs</td>
<td>115,114</td>
<td>42,846,307</td>
<td>372</td>
<td>776</td>
</tr>
<tr><td align="justify">Unigenes</td>
<td>77,292</td>
<td>85,969,539</td>
<td>1112</td>
<td>1778</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="Fig1"><label>Fig. 1</label>
<caption><p>Length distributions of contigs (<bold>a</bold>
) and unigenes (<bold>b</bold>
) in trifoliate orange transcriptome</p>
</caption>
<graphic xlink:href="12864_2015_1629_Fig1_HTML" id="MO1"></graphic>
</fig>
</p>
<p>Gene expression levels were calculated based on FPKM (fragments per kb per million fragments) method [<xref ref-type="bibr" rid="CR26">26</xref>
], and among the unigenes, 28,393 (36.7 %) and 65,372 (84.5 %) had FPKM values smaller than 1.0 and 10.0, respectively, whereas only 1241 (1.6 %) had FPKM values > 100. This suggests that most of the unigenes were expressed at low levels.</p>
</sec>
<sec id="Sec4"><title>Functional annotation of unigenes</title>
<p>The assembled unigene sequences were used to search public databases, including the NCBI non-redundant protein database (nr), the Swissprot protein database, the Gene Ontology (GO) database, the Clusters of Orthologous Groups (COGs) database, and Kyoto Encyclopedia of Genes and Genomes (KEGG) database, using the BLASTX algorithm with an <italic>E</italic>
-value threshold of 1e<sup>−5</sup>
. In total, 59,777 (77.3 %) unigenes were matched to a sequence in at least one of the above-mentioned databases (Table <xref rid="Tab2" ref-type="table">2</xref>
) and 58,001 and 36,445 unigenes had significant hits (<italic>E</italic>
-value < 1e<sup>−5</sup>
) in the Nr and Swissprot non-redundant protein databases, respectively. In terms of <italic>E</italic>
-value distribution, 37.4 % of the homologs ranged between 1e<sup>−5</sup>
and 1e<sup>−45</sup>
, while a majority of the sequences (62.5 %) showed a threshold e-value < 1e<sup>−45</sup>
indicating strong homology (Fig. <xref rid="Fig2" ref-type="fig">2a</xref>
). Among the unigenes annotated using the Nr database, 74.6 % had more than 60 % similarity with the corresponding gene sequence (Fig. <xref rid="Fig2" ref-type="fig">2b</xref>
). Overall, the unigene sequences exhibited most similar BLASTx matches to gene sequences from <italic>Vitis vinifera</italic>
(31.8 %), followed by those from <italic>Ricinus communis</italic>
(26.5 %), <italic>Populus trichocarpa</italic>
(22.2 %), <italic>Glycine max</italic>
(6.3 %), <italic>Medicago truncatula</italic>
(2.0 %) and <italic>A. thaliana</italic>
(1.4 %) genes (Fig. <xref rid="Fig2" ref-type="fig">2c</xref>
).<table-wrap id="Tab2"><label>Table 2</label>
<caption><p>List of <italic>Poncirus trifoliata</italic>
transcriptome annotations</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th>Public database</th>
<th>No. of unigene hits</th>
<th>Percentage (%)</th>
</tr>
</thead>
<tbody><tr><td align="justify">Nt</td>
<td align="justify">53,702</td>
<td align="justify">69.5</td>
</tr>
<tr><td align="justify">Nr</td>
<td align="justify">58,001</td>
<td align="justify">75.0</td>
</tr>
<tr><td align="justify">SwissProt</td>
<td align="justify">36,445</td>
<td align="justify">47.2</td>
</tr>
<tr><td align="justify">GO</td>
<td align="justify">34,747</td>
<td align="justify">45.0</td>
</tr>
<tr><td align="justify">COG</td>
<td align="justify">23,846</td>
<td align="justify">30.9</td>
</tr>
<tr><td align="justify">KEGG</td>
<td align="justify">45,819</td>
<td align="justify">59.3</td>
</tr>
<tr><td align="justify">ALL</td>
<td align="justify">59,777</td>
<td align="justify">77.3</td>
</tr>
</tbody>
</table>
<table-wrap-foot><p>Nt: Nucleotide database; Nr: Non-redundant protein sequence database</p>
<p>SwissProt: Swiss-Prot protein sequence database; GO: Gene Onotology database</p>
<p>COG; Cluster of Orthologous Groups of proteins</p>
<p>KEGG: Kyoto Encyclopedia of Genes and Genomes</p>
</table-wrap-foot>
</table-wrap>
<fig id="Fig2"><label>Fig. 2</label>
<caption><p>Characteristics of homology search of unigenes against NCBI non-redundant (Nr) database. <bold>a</bold>
. <italic>E</italic>
-value distribution of BLAST hits for each unique sequence with a cut-off <italic>E</italic>
-value of 1e<sup>−5</sup>
. <bold>b</bold>
. Similarity distribution of top BLAST hits for each unigene. <bold>c</bold>
. Species distribution of the top BLAST hits for each unigene with a cut-off <italic>E</italic>
-value of 1e<sup>−5</sup>
</p>
</caption>
<graphic xlink:href="12864_2015_1629_Fig2_HTML" id="MO2"></graphic>
</fig>
</p>
<p>GO analysis was widely used to assign putative gene function to uncharacterized sequences. We used the information from the Nr annotation together with the WEGO program [<xref ref-type="bibr" rid="CR27">27</xref>
] to categorize 34,747 unigenes into 61 functional groups (Fig. <xref rid="Fig3" ref-type="fig">3</xref>
), and these were further divided among three GO terms: ‘including molecular function’, ‘cellular component’, and ‘biological process’ using BLAST2GO [<xref ref-type="bibr" rid="CR28">28</xref>
]. Most of the unigenes in the ‘molecular function’ category were sub-categorized into ‘catalytic activity’ (25,902) and ‘binding’ (24,250), followed by ‘transporter activity’ (3890) and ‘nucleic acid binding transcription factor activity’ (1400). In the ‘cellular component’ category, the GO terms ‘cell’ (31,049), ‘cell part’ (31,049), ‘organelle part’ (23,958), and ‘membrane’ (14,312) predominated, while ‘cellular process’ (27,062), ‘metabolic process’ (26,810), and ‘response to stimulus’ (12,031) were the three most represented GO terms in the ‘biological process’ category.<fig id="Fig3"><label>Fig. 3</label>
<caption><p>Histogram of Gene Ontology (GO) classifications. Percentages indicate the proportion of unigenes that have the relevant GO annotations, including biological process, cellular components and molecular function</p>
</caption>
<graphic xlink:href="12864_2015_1629_Fig3_HTML" id="MO3"></graphic>
</fig>
</p>
<p>We then performed a phylogenetic classification of the trifoliate orange sequences using the Cluster of Orthologous Groups (COG) database, which is composed of protein sequences derived from the genomes of bacteria, plants, and animals. Each protein in the COG database is assumed to have evolved from an ancestral protein. On the basis of sequence similarities, of the 36,445 SwissProt sequence matches to the trifoliate orange sequences, 23,846 unigenes could be grouped into 25 COG categories (Fig. <xref rid="Fig4" ref-type="fig">4</xref>
), Of these ‘general function prediction only’ (7797; 32.7 %) represented the largest class, followed by ‘transcription’ (3786; 15.9 %), ‘replication, recombination and repair’ (3695; 15.5 %), and ‘signal transduction mechanisms’ (3283; 13.8 %), whereas 1736 (5.8 %) were assigned to the category of ‘function unknown’ category. In addition, only a few unigenes grouped in the ‘extracellular structures’ (7, 0.03 %) and ‘nuclear structure’ (14, 0.06 %) categories.<fig id="Fig4"><label>Fig. 4</label>
<caption><p>Distribution of genes with COG functional classification in the transcriptome. A total of 23,846 sequences have a COG classification among 25 categories. A: RNA processing and modification; B: Chromatin structure and dynamics; C: Energy production and conversion; D: Cell cycle control, cell division, chromoso; E: Amino acid transport and metabolism; F: Nucleotide transport and metabolism; G: Carbohydrate transport and metabolism; H: Coenzyme transport and metabolism; I: Lipid transport and metabolism; J: Translation, ribosomal structure and biogenesis; K: Transcription; L: Replication, recombination and repair; M: Cell wall/membrane/envelope biogenesis; N: Cell motility; O: Posttranslational modification, protein turnover, chaperones; P: Inorganic ion transport and metabolism; Q: Secondary metabolites biosynthesis, transport and catabolism; R: General function prediction only; S: Function unknown; T: Signal transduction mechanisms; U: Intracellular trafficking, secretion, and vesicular transport; V: Defense mechanisms; W: Extracellular structures; Y: Nuclear structure; Z: Cytoskeleton</p>
</caption>
<graphic xlink:href="12864_2015_1629_Fig4_HTML" id="MO4"></graphic>
</fig>
</p>
<p>When the assembled unigenes were analyzed using the BLASTX function of the KEGG database using BLASTX. As a result, 45,819 (59.3 %) mapped to 128 KEGG pathways (Additional file <xref rid="MOESM1" ref-type="media">1</xref>
), of which ‘metabolic pathway’ (7987; 23.0 %), ‘biosynthesis of secondary metabolites’ (3923; 11.3 %), ‘plant-pathogen interaction’ (2111; 6.1 %), ‘plant hormone signal transduction’ (1693; 4.9 %) and ‘RNA transport’ (1067; 3.1 %) were the most abundant groups. Interestingly, we noted that a large number of unigenes involved in the biosynthesis of primary metabolites, such as starch, sucrose, arginine, and proline, all of which have been reported to be associated with stress tolerance, were also enriched in the dataset (Additional file <xref rid="MOESM2" ref-type="media">2</xref>
).</p>
<p><italic>P. trifoliata</italic>
protein-coding sequences were deduced using BLASTx and the unigene sequences to interrogate the Nr, SwissProt, KEGG, and COG, using BLASTx protein databases, resulting in a total of 57,934 predicted coding sequences (CDSs). Among the unigenes with CDSs, 49.6 % were 500–2000 bp, while 44.6 % were < 500 bp, and 811 were > 3000 bp (Fig. <xref rid="Fig5" ref-type="fig">5</xref>
). A total of 889 CDSs did not have a match in any of the protein databases using ESTScan program [<xref ref-type="bibr" rid="CR29">29</xref>
].<fig id="Fig5"><label>Fig. 5</label>
<caption><p>Coding sequence (CDS) of unigenes predicted by BLASTx and ESTScan. The size distribution of the CDSs and proteins based on BLASTx (<bold>a</bold>
) and ESTScan (<bold>b</bold>
)</p>
</caption>
<graphic xlink:href="12864_2015_1629_Fig5_HTML" id="MO5"></graphic>
</fig>
</p>
</sec>
<sec id="Sec5"><title>Identification of differentially expressed genes (DEGs)</title>
<p>In order to gain a global view of transcript expression in trifoliate orange in response to cold stress and to identify cold-responsive genes, we constructed four cDNA libraries from plants either grown under normal conditions or exposed to a cold treatment (4 °C for 6, 24, or 72 h). After data filtering, at least 7.1 million 49 nt single-end clean reads were obtained for each library, which represented ~99.3 % of the raw data (Table <xref rid="Tab3" ref-type="table">3</xref>
). An analysis of sequencing saturation indicated that the output of cleaned reads was sufficient to identify DEGs (Additional file <xref rid="MOESM3" ref-type="media">3</xref>
). The filtered cleaned reads were then mapped onto the assembled <italic>P. trifoliata</italic>
transcriptome using SOAPaligner/soap2 software [<xref ref-type="bibr" rid="CR30">30</xref>
]. The uniquely matched reads, which ranged from 45.5 % to 49.3 % of the total cleaned reads, were subsequently used to assess relative expression levels (Table <xref rid="Tab3" ref-type="table">3</xref>
). Each transcript had a range of read coverage, but more than 40 % of the detected unigenes had > 50 % coverage (data not shown). According to the criteria used for identifying DEGs, including low false discovery rate (FDR, less than 10<sup>−3</sup>
) and |log<sub>2</sub>
FPKM|> 1.0, a total of 5549 unigenes were considered to be differentially expressed between control and cold treated plants, with 600 at in the 6 h sample (462 up-regulated, 138 down-regulated), 2346 in the 24 h sample (1631 up-regulated, 715 down-regulated), and 5475 in the 72 h sample (2702 up-regulated, 2473 down-regulated) (Fig. <xref rid="Fig6" ref-type="fig">6a</xref>
, <xref rid="Fig6" ref-type="fig">b</xref>
, Additional file <xref rid="MOESM4" ref-type="media">4</xref>
and Additional file <xref rid="MOESM2" ref-type="media">2</xref>
). When the DEGs at different time points were evaluated (Fig. <xref rid="Fig6" ref-type="fig">6c</xref>
), we observed that their numbers increased with longer cold treatments, and that the number of DEGs between two different cold treatment samples (24 vs 6, 72 vs. 6 and 72 vs. 24) was relatively small.<table-wrap id="Tab3"><label>Table 3</label>
<caption><p>Sequencing and read mapping of the differentially expressed genes (DEGs)</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th>Summary</th>
<th>0 h</th>
<th>6 h</th>
<th>24 h</th>
<th>72 h</th>
<th>Average</th>
</tr>
</thead>
<tbody><tr><td align="justify">Raw reads</td>
<td align="justify">7,164,446</td>
<td align="justify">7,325,695</td>
<td align="justify">7,563,619</td>
<td align="justify">7,212,252</td>
<td align="justify">7,316,503</td>
</tr>
<tr><td align="justify">Clean reads</td>
<td align="justify">7,113,843</td>
<td align="justify">7,271,861</td>
<td align="justify">7,515,414</td>
<td align="justify">7,163,347</td>
<td align="justify">7,266,116</td>
</tr>
<tr><td align="justify">Total base pairs</td>
<td align="justify">348,578,307</td>
<td align="justify">356,321,189</td>
<td align="justify">368,255,286</td>
<td align="justify">351,004,003</td>
<td align="justify">356,039,696</td>
</tr>
<tr><td align="justify">Clean Reads (%)</td>
<td align="justify">99.29</td>
<td align="justify">99.27</td>
<td align="justify">99.36</td>
<td align="justify">99.32</td>
<td align="justify">99.31</td>
</tr>
<tr><td align="justify">Mapped reads</td>
<td align="justify">6,621,141</td>
<td align="justify">6,727,919</td>
<td align="justify">6,879,216</td>
<td align="justify">6,467,083</td>
<td align="justify">6,467,087</td>
</tr>
<tr><td align="justify">Mapped reads* (%)</td>
<td align="justify">93.07</td>
<td align="justify">92.52</td>
<td align="justify">91.53</td>
<td align="justify">90.28</td>
<td align="justify">91.85</td>
</tr>
<tr><td align="justify">Perfect match</td>
<td align="justify">5,583,419</td>
<td align="justify">5,707,629</td>
<td align="justify">5,813,295</td>
<td align="justify">5,505,742</td>
<td align="justify">5,652,521</td>
</tr>
<tr><td align="justify">Perfect match* (%)</td>
<td align="justify">78.49</td>
<td align="justify">78.49</td>
<td align="justify">77.35</td>
<td align="justify">76.83</td>
<td align="justify">77.79</td>
</tr>
<tr><td align="justify"><=2 Mismatch</td>
<td align="justify">1,037,722</td>
<td align="justify">1,020,290</td>
<td align="justify">1,065,921</td>
<td align="justify">961,341</td>
<td align="justify">1,021,319</td>
</tr>
<tr><td align="justify"><=2 Mismatch* (%)</td>
<td align="justify">14.59</td>
<td align="justify">14.03</td>
<td align="justify">14.18</td>
<td align="justify">13.42</td>
<td align="justify">14.06</td>
</tr>
<tr><td align="justify">Unique match</td>
<td align="justify">3,117,820</td>
<td align="justify">3,309,503</td>
<td align="justify">3,456,940</td>
<td align="justify">3,290,477</td>
<td align="justify">3,293,685</td>
</tr>
<tr><td align="justify">Unique match* (%)</td>
<td align="justify">43.83</td>
<td align="justify">45.51</td>
<td align="justify">46.00</td>
<td align="justify">45.93</td>
<td align="justify">45.33</td>
</tr>
<tr><td align="justify">Multi-Position match</td>
<td align="justify">3,503,321</td>
<td align="justify">3,418,416</td>
<td align="justify">3,422,276</td>
<td align="justify">3,176,606</td>
<td align="justify">3,380,155</td>
</tr>
<tr><td align="justify">Multi-Position match</td>
<td align="justify">49.25</td>
<td align="justify">47.01</td>
<td align="justify">45.54</td>
<td align="justify">44.35</td>
<td align="justify">46.54</td>
</tr>
<tr><td align="justify">Unmapped reads</td>
<td align="justify">492,702</td>
<td align="justify">543,942</td>
<td align="justify">636,198</td>
<td align="justify">696,264</td>
<td align="justify">592,277</td>
</tr>
<tr><td align="justify">Unmapped reads* (%)</td>
<td align="justify">6.93</td>
<td align="justify">7.48</td>
<td align="justify">8.47</td>
<td align="justify">9.72</td>
<td align="justify">8.15</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="Fig6"><label>Fig. 6</label>
<caption><p>Differentially expressed genes (DEGs) of <italic>Poncirus trifoliata</italic>
under cold stress. A-B. Venn maps of up-regulated (<bold>a</bold>
) and down-regulated genes (<bold>b</bold>
). C. The number of up-regulated and down-regulated DEGs compared between different time points under cold treatment at 0, 6, 24 and 72 h</p>
</caption>
<graphic xlink:href="12864_2015_1629_Fig6_HTML" id="MO6"></graphic>
</fig>
</p>
</sec>
<sec id="Sec6"><title>Gene Ontology clustering of DEGs</title>
<p>The DEGs were first mapped to the three main categories of the GO database: ‘molecular function’, ‘cellular component’ and ‘biological process’. Significantly enriched GO terms then determined using a hypergeometric test to identify functional clusters and biochemical pathways (corrected <italic>P</italic>
-value). In the ‘biological process’ category, 18 GO terms were significantly enriched in the 6 h cold treatment sample (Additional file <xref rid="MOESM5" ref-type="media">5</xref>
), of which the GO terms ‘response to stimulus’, ‘response to organic substance’, ‘response to stress’, ‘response to carbohydrate stimulus’, ‘response to chemical stimulus’ and ‘response to hormone stimulus’ were most highly enriched. We interpreted this to indicate an involvement of the genes mapping to these pathways in perception and transduction of the cold signal. The substantial up-regulation of certain kinases (CL6797.Contig3_All, Unigene16428_All, CL8004.Contig1_All) and TFs (Unigene20109_All, Unigene11424_All, Unigene14645_All, Unigene19023_All, Unigene21445_All) at 6 h was particularly notable. Also of interest was the fact that the GO terms ‘response to chitin’, ‘response to fungus’, ‘response to other organism’, ‘response to biotic stimulus’, ‘response to wounding’, ‘defense response’, ‘response to salicylic acid stimulus’ and ‘defense response to fungus’ were enriched in the 6 h sample, which suggests pathway crosstalk between biotic stress responses and the perception and transduction of cold induced signaling.</p>
<p>In the 24 h cold treatment sample the enriched GO terms in the ‘biological process’ category were similar to those at 6 h; however, more DEGs were identified in the various categories. Lastly, in the 72 h cold treatment sample, we observed that in addition to the stress-responsive categories, the two photosynthesis-related GO terms, ‘photosynthesis, light harvesting’ and ‘photosynthesis’ were over-represented. Interestingly, all of the DEGs related to photosynthesis were down-regulated upon exposure to the cold treatment, suggesting that the photosynthetic capacity was suppressed. In addition, the categories ‘RNA secondary structure unwinding’ and ‘translation’ were enriched in the 72 h sample, from which we infer that regulation of gene expression in response to cold temperatures may occur at the post-transcriptional and translational levels. Four unigenes enriched in ‘RNA secondary structure unwinding’ showed an increase in abundance, and three of these were annotated as glycine-rich RNA-binding protein s, which have been shown in <italic>A. thaliana</italic>
to be involved in biotic or abiotic stress responses [<xref ref-type="bibr" rid="CR31">31</xref>
, <xref ref-type="bibr" rid="CR32">32</xref>
]. The fourth gene was annotated as cellular nucleic acid binding protein.</p>
</sec>
<sec id="Sec7"><title>Pathways enrichment analysis of DEGs</title>
<p>Using the KEGG database, pathways displaying significant changes (Q value ≤ 0.05) in response to the cold treatment were identified (Table <xref rid="Tab4" ref-type="table">4</xref>
). In the 6 h cold treatment sample, 218 of the 600 DEGs were associated with KEGG pathways. This analysis revealed five significantly enriched pathways: ‘plant-pathogen interaction’, ‘plant hormone signal transduction’, ‘riboflavin metabolism’, ‘flavone and flavonol biosynthesis’ and ‘glycosyl phosphatidyl inositol-anchor biosynthesis’. In the 24 h cold treatment sample, enriched pathways included ‘plant-pathogen interaction’, ‘plant hormone signal transduction’, ‘ether lipid metabolism’, ‘GPI-anchor biosynthesis’, and ‘zeatin biosynthesis’ were significantly changed and more than 12 pathways were significantly changed after 72 h of cold treatment, among which ‘photosynthesis-antenna proteins’, ‘plant-pathogen interaction’, ‘photosynthesis’, ‘flavone and flavonol biosynthesis’ and ‘stilbenoid, diarylheptanoid and gingerol biosynthesis’ were the most significantly enriched. Interestingly, genes involved in the ‘plant-pathogen interaction’ pathway, such as <italic>CDPK</italic>
(CL4696.Contig2_All), calcium-binding protein (Unigene18341_All, Unigene4416_All), and <italic>WRKY</italic>
(CL3620.Contig1_All, Unigene11423_All) showed large changes in expression levels at all three points of the cold treatment (Table <xref rid="Tab4" ref-type="table">4</xref>
). Finally, the ‘plant hormone signal transduction’ pathway showed significant changes in all three samples, suggesting that hormone signaling is involved in the cold stress response of <italic>P. trifoliata</italic>
.<table-wrap id="Tab4"><label>Table 4</label>
<caption><p>Significantly enriched gene pathways involving differentially expressed genes (DEGs) following the cold stress treatment</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th></th>
<th>Pathway</th>
<th>DEGs with pathway annotation</th>
<th>All genes with pathway annotation (34747)</th>
<th>Q value</th>
<th>Pathway ID</th>
</tr>
</thead>
<tbody><tr><td>6 h vs 0 h</td>
<td>Plant-pathogen interaction</td>
<td>37</td>
<td>2111 (6.1 %)</td>
<td>8.66e-07</td>
<td>ko04626</td>
</tr>
<tr><td></td>
<td>Plant hormone signal transduction</td>
<td>30</td>
<td>1693 (4.9 %)</td>
<td>9.60e-06</td>
<td>ko04075</td>
</tr>
<tr><td></td>
<td>Riboflavin metabolism</td>
<td>5</td>
<td>101 (0.3 %)</td>
<td>9.47e-03</td>
<td>ko00740</td>
</tr>
<tr><td></td>
<td>Flavone and flavonol biosynthesis</td>
<td>6</td>
<td>164 (0.5 %)</td>
<td>9.70e-03</td>
<td>ko00944</td>
</tr>
<tr><td></td>
<td>Glycosylphosphatidylinositol(GPI)-anchor biosynthesis</td>
<td>6</td>
<td>229 (0.7 %)</td>
<td>4.23e-02</td>
<td>ko00563</td>
</tr>
<tr><td>24 h vs 0 h</td>
<td>Plant-pathogen interaction</td>
<td>125</td>
<td>2111 (6.1 %)</td>
<td>1.35e-16</td>
<td>ko04626</td>
</tr>
<tr><td></td>
<td>Plant hormone signal transduction</td>
<td>86</td>
<td>1693 (4.9 %)</td>
<td>6.46e-08</td>
<td>ko04075</td>
</tr>
<tr><td></td>
<td>Ether lipid metabolism</td>
<td>21</td>
<td>372 (1.1 %)</td>
<td>2.27e-02</td>
<td>ko00565</td>
</tr>
<tr><td></td>
<td>Zeatin biosynthesis</td>
<td>15</td>
<td>245 (0.7 %)</td>
<td>3.55e-02</td>
<td>ko00908</td>
</tr>
<tr><td></td>
<td>Glycosylphosphatidylinositol(GPI)-anchor biosynthesis</td>
<td>15</td>
<td>229 (0.7 %)</td>
<td>2.27e-02</td>
<td>ko00563</td>
</tr>
<tr><td>72 h vs 0 h</td>
<td>Photosynthesis - antenna proteins</td>
<td>13</td>
<td>26 (0.1 %)</td>
<td>1.12e-07</td>
<td>ko00196</td>
</tr>
<tr><td></td>
<td>Plant-pathogen interaction</td>
<td>190</td>
<td>2111 (6.1 %)</td>
<td>6.43e-06</td>
<td>ko04626</td>
</tr>
<tr><td></td>
<td>Photosynthesis</td>
<td>22</td>
<td>106 (0.3 %)</td>
<td>1.89e-05</td>
<td>ko00195</td>
</tr>
<tr><td></td>
<td>Flavone and flavonol biosynthesis</td>
<td>25</td>
<td>164 (0.5 %)</td>
<td>8.39e-04</td>
<td>ko00944</td>
</tr>
<tr><td></td>
<td>Diterpenoid biosynthesis</td>
<td>17</td>
<td>100 (0.3 %)</td>
<td>2.11e-03</td>
<td>ko00904</td>
</tr>
<tr><td></td>
<td>Ether lipid metabolism</td>
<td>42</td>
<td>372 (1.1 %)</td>
<td>2.11e-03</td>
<td>ko00565</td>
</tr>
<tr><td></td>
<td>Plant hormone signal transduction</td>
<td>142</td>
<td>1693 (4.9 %)</td>
<td>2.11e-03</td>
<td>ko04075</td>
</tr>
<tr><td></td>
<td>Ribosome</td>
<td>46</td>
<td>431 (1.2 %)</td>
<td>3.27e-03</td>
<td>ko03010</td>
</tr>
<tr><td></td>
<td>Cutin, suberine and wax biosynthesis</td>
<td>19</td>
<td>129 (0.4 %)</td>
<td>4.61e-03</td>
<td>ko00073</td>
</tr>
<tr><td></td>
<td>Limonene and pinene degradation</td>
<td>33</td>
<td>328 (0.9 %)</td>
<td>4.37e-02</td>
<td>ko00903</td>
</tr>
<tr><td></td>
<td>Endocytosis</td>
<td>61</td>
<td>664 (1.9 %)</td>
<td>1.58e-02</td>
<td>ko04144</td>
</tr>
<tr><td></td>
<td>Stilbenoid, diarylheptanoid and gingerol biosynthesis</td>
<td>39</td>
<td>330 (1.0 %)</td>
<td>2.10e-03</td>
<td>ko00945</td>
</tr>
</tbody>
</table>
</table-wrap>
</p>
</sec>
<sec id="Sec8"><title>Validation of RNA-seq based DEG results by quantitative real-time RT-PCR (qPCR)</title>
<p>To validate the RNA-seq data, qPCR analysis of transcript abundance was performed of 17 randomly selected DEGs, including 12 that were up-regulated and 5 that were shown to be down-regulated in the unigene dataset (Additional file <xref rid="MOESM2" ref-type="media">2</xref>
). As shown in Fig. <xref rid="Fig7" ref-type="fig">7</xref>
, the fold change values obtained by qPCR were highly consistent with those based on the RNA-Seq data for all of the tested genes (Fig. <xref rid="Fig7" ref-type="fig">7</xref>
), except for a difference in the expression patterns of two down-regulated genes at two time points. This result supports the reliability of the RNA-Seq analysis.<fig id="Fig7"><label>Fig. 7</label>
<caption><p>qPCR analysis of differentially expressed genes in <italic>Poncirus trifoliata</italic>
under low temperature. (<bold>a</bold>
, <bold>b</bold>
) Transcript levels of 17 randomly selected DEGs, including 12 up-regulated (<bold>a</bold>
) and 5 down-regulated (<bold>b</bold>
). The Y-axis on the left shows the relative gene expression levels analyzed by qPCR (red lines), while Y-axis on the right shows corresponding expression data of RNA-seq (gray histogram). The X-axis represents the time (hours) of 4 °C treatment. The bars represent SE (n = 3). (<bold>c</bold>
) Comparison between the log<sub>2</sub>
of gene expression ratios obtained from RNA-seq data and qRT-PCR</p>
</caption>
<graphic xlink:href="12864_2015_1629_Fig7_HTML" id="MO7"></graphic>
</fig>
</p>
</sec>
<sec id="Sec9"><title>Expression of genes involved in polyamine and ethylene metabolism</title>
<p>Polyamines (PAs), including diamine putrescine (Put), triamine spermidine (Spd) and tetramine spermine (Spm), are polycationic low molecular weight compounds that are ubiquitoius in most organisms. Several unigenes encoding polyamine biosynthetic enzymes, such as <italic>ADC</italic>
(arginine decarboxylase), <italic>ODC</italic>
(ornithine decarboxylase), <italic>SAMDC</italic>
(S-adenosylmethionine decarboxylase), <italic>SPDS</italic>
(spermidine synthase) and <italic>SPMS</italic>
(spermine synthase) were identified, and of these <italic>ADC</italic>
(Unigene25843_All), which encodes a key enzyme in Put biosynthesis, was highly expressed in the cold treated samples compared to the control samples. Meanwhile, the expression of <italic>ODC</italic>
(Unigene17207_All, Unigene11015_All, CL6768.Contig1_All, CL9042.Contig1_All), an alternative enzyme involved in Put biosynthesis, showed little change in response to cold exposure. This finding indicates that the <italic>ADC</italic>
pathway may be more important for the Put biosynthesis of <italic>P. trifoliata</italic>
than the <italic>ODC</italic>
one during cold stress. SAMDC is involved in the synthesis of spermidine and spermine by providing the aminopropyl moiety derived from a decarboxylation of <italic>S</italic>
-adenosyl methionine (SAM). We observed that the <italic>P. trifoliata SAMDC</italic>
unigene (Unigene20820_All) was expressed at higher levels after the low temperature treatment. Meanwhile, unigenes with homology to <italic>SPDS</italic>
, <italic>SPMS</italic>
, and <italic>ACL5</italic>
showed no increase in expression levels during the cold treatment. Furthermore, the expression of SAM synthases (CL5314.Contig1_All, Unigene24811_All, and Unigene21190_All) was substantially induced. Furthermore, aminocyclopropane carboxylic acid (ACC) synthase (ACS) converts SAM to ACC, which is then oxidized into ethylene by ACC oxidase (ACO). Notably, an <italic>ACS</italic>
gene (CL4785.Contig1_All) was down-regulated during cold stress, while expression of an <italic>ACO</italic>
(CL6929.Contig1_All) gene was barely detectable over the entire time course. We also measured the content of free polyamines in the cold-treated samples. Put levels were progressively elevated under cold treatment, while levels of Spd and Spm underwent only minor changes with the exception of a slight increase of Spd at 24 h and Spm at 6 h (Fig. <xref rid="Fig8" ref-type="fig">8</xref>
).<fig id="Fig8"><label>Fig. 8</label>
<caption><p>Quantification of free polyamines (PAs) of <italic>Poncirus trifoliata</italic>
under cold treatment. Three types of free polyamines, including putrescine, spermine, and spermidine, are assessed by HPLC. The bars represent SE (n = 3). Asterisks show that the values are significantly when compared with the levels at the onset of cold treatment (**<italic>P</italic>
< 0.01)</p>
</caption>
<graphic xlink:href="12864_2015_1629_Fig8_HTML" id="MO8"></graphic>
</fig>
</p>
</sec>
<sec id="Sec10"><title>Identification of transcription factors (TFs) in response to cold stress</title>
<p>We determined that the transcript abundance of 60 TFs exhibited highly dynamic changes in response to the cold treatment (Fig. <xref rid="Fig9" ref-type="fig">9</xref>
), of which those encoding AP2/ERF domain-containing proteins constituted the largest group (43.3 %), followed by WRKY proteins (1.7 %), NACs (1.5 %), zinc finger proteins (ZFPs, 0.7 %), MYBs (0.7 %), and bHLHs (0.7 %). In addition, the expression of two heat shock transcription factors (Hsfs) (Unigene12132, CL9882.Contig2) was found to be higher following the cold treatments (Additional file <xref rid="MOESM6" ref-type="media">6</xref>
). Transcript levels of Unigene12132_All were induced by 1.7, 3.8, and 4.5 folds at 6, 24 and 72 h in comparison with that of control, while CL9882.Contig2_All showed a progressive increase and was 2.8 fold higher than the control at 72 h.<fig id="Fig9"><label>Fig. 9</label>
<caption><p>A heat map indicating expression patterns of transcription factors under cold stress. Columns and rows in the heat map represent samples collected at differenet time points of cold stress and common differenttially expressed TFs, respectively</p>
</caption>
<graphic xlink:href="12864_2015_1629_Fig9_HTML" id="MO9"></graphic>
</fig>
</p>
<p>The differentially expressed TFs were classified into six subgroups based on their expression patterns (Fig. <xref rid="Fig9" ref-type="fig">9</xref>
). Subgroup I was characterized by genes showing low changes in transcript levels early during the cold treatment compared with the control plants, a peak of expression at 24 h and then a decrease at the 72 h time point. TFs in subgroup II showed a relatively high fold change at all three time points and five TFs belonged to this group, including three AP2/ERF domain-containing TFs, one NAC and one WRKY. The four TFs in subgroup III showed the highest expressional levels at 6 and 24 h, but lower levels at 72 h. The transcript abundance of TFs in subgroup IV, including four AP2/ERF domain-containing TFs, two zinc finger and a single WRKY, was greater at the early time point before declining, which suggests that they play roles in early transduction of the cold signal. TFs in subgroup V showed a peak in expression at 72 h, while TFs in subgroup VI, the largest subgroup, showed high levels of transcript abundance at 24 h and the levels remained high at 72 h.</p>
<p>Among the differentially expressed AP2/ERF TFs, four CBFs/DREBs-like unigenes (Unigene5199_All, Unigene36912_All, Unigene24431_All and Unigene14645_All) which were all clustered in subgroup I, and therefore had a similar expression pattern, were significantly induced by cold treatment. For example, Unigene24431_All, which was annotated as a <italic>DREB2A</italic>
-like gene, had a RPKM value of 35.79, 139.60, 497.48 and 1085.57 in the 0, 6, 24 and 72 h samples. In addition, the expression of an ICE1-like TF (Unigene24033_All) was also found to be up-regulated in the cold-treated plants.</p>
</sec>
<sec id="Sec11"><title>Identification of hormone and signal transduction-related unigenes</title>
<p>The GO and KEGG enrichment analyses of the DEGs both highlighted signal transduction as being particularly affected by the cold treatment. More than 60 unigenes predicted to encode protein kinases were differentially expressed in trifoliate orange grown under the cold conditions (Additional file <xref rid="MOESM7" ref-type="media">7</xref>
) and of these, genes encoding receptor-like protein kinases (RLKs) made up the largest group. More than 30 RLK genes were serine/threonine-protein kinase, three were wall-associated receptor kinases (Unigene20744_All, Unigene25402_All, Unigene22769_All), and one was a proline-rich receptor-like protein kinase (Unigene20499_All). Furthermore, the expression of <italic>BRASSINOSTEROIDS INSENSITIVE</italic>
1 (BRI1, CL3291.Contig2_All) was found to increase, while an annotated brassinosteroid LRR receptor kinase (Unigene19522_All) showed a decrease in transcript abundance under cold stress. Additionally, eight genes annotated as kinases were up-regulated. Mitogen-activated protein kinase (MAPK) cascades have been shown to play key roles in mediating stress-associated signaling pathway [<xref ref-type="bibr" rid="CR33">33</xref>
] and among the DEGs, a MAPK (CL8004.Contig1_All) and two MAPKKK (Unigene25596_All, CL3174.Contig_All) were up-regulated by the cold treatment. In addition, the expression of several kinases related to MAPK (mitogen-activated protein kinase) cascades was similarly induced.</p>
<p>We also observed that the expression of genes associated with biosynthesis or signal transduction of the phytohormones abscisic acid (ABA), gibberellin (GA) or ethylene showed substantial changes following cold temperatures. KEGG analysis of the DEGs suggested that carotenoid biosynthesis (ko00906), which is associated with ABA biosynthesis, was induced under cold condition. As an example, the expression of Unigene19385_All, which is predicted to encode a 9-<italic>cis</italic>
-epoxycarotenoid dioxygenase (NCED), a key rate-limiting enzyme in ABA synthesis, was progressively up-regulated by ~3-4 folds at 24 h and 72 h, respectively, compared to the 0 h time point. Interestingly, the expression of a unigene (CL851.Contig21_All) encoding an ABA 8′-hydroxylase 1, a key enzyme in ABA catabolism, was rapidly up-regulated after 6 h of cold treatment, but its transcript abundance declined to the same expressional level as 0 h after 72 h. This apparent regulation of genes involved in ABA biosynthesis and catabolism suggests dynamic and multifaceted control of ABA content in response to stress. Also we also noted genes involved in ABA signaling pathway showed similar changes in expression. For example, three genes encoding protein phosphatase 2C (PP2C), CL6797.Contig2_All, CL6797.Contig1_All and Unigene24667_All, were up-regulated at 6 h, and four additional <italic>PP2C</italic>
genes (CL3200.Contig4_All, CL9817.Contig1_All, Unigene24667_All, CL2118.Contig1_All and Unigene934_All) showed increase in transcript abundance after 24 h of cold stress. In contrast, three other <italic>PP2C</italic>
genes (Unigene22810_All, Unigene21083_All, and Unigene23511_All) and a unigene (Unigene23445_All) encoding PYR/PYL (Pyrabatin Resistance/Pyrabatin Resistance 1-Like), ABA receptors, were down-regulated at 72 h.</p>
<p>Although expression of the genes associated with ethylene biosynthesis was generally down-regulated by the cold treatment, we identified unigenes related to ethylene meditated signaling that showed the opposite trend. This was exemplified by 11 ethylene signal transduction-related genes, three ethylene receptors (ETRs), a MAPK and seven ERFs, all of which were up-regulated by cold exposure. However, we also determined that a unigene (Unigene13110_All) encoding an Ein3-binding F-box protein (EBF1/2) (Unigene13110_All) was down-regulated.</p>
<p>We also observed that the expression of six genes related to biosynthesis of the defense-related hormone jasmonic acid (JA) was upregulated: one encoding allene oxide cyclase (Unigene15722_All), two encoding 12-oxophytodienoate reductase 2 (CL5141.Contig4_All and CL5141.Contig2_All) and two encoding AMP dependent CoA ligase (CL2829.Contig1_All, Unigene24432_All). Furthermore, three genes encoding JAZ and five genes encoding MYC2, which have been shown to be involved in JA signaling pathway, were also significantly up-regulated by the cold conditions. The expression patterns of the genes involved in hormone metabolism and signal transduction were analyzed by qPCR (Fig. <xref rid="Fig10" ref-type="fig">10</xref>
). Although the samples were different from those used for RNA-Seq, expression patterns were quite similar between the two methods, which further confirmed reliability of screening the DEGs.<fig id="Fig10"><label>Fig. 10</label>
<caption><p>Expression patterns of genes involved in hormone metabolism and signal transduction. Expression of the selected genes was comfirmed by qPCR using samples treated at low temperature. The samples were different from those used for RNA-Seq. The bars represent SE (n = 3)</p>
</caption>
<graphic xlink:href="12864_2015_1629_Fig10_HTML" id="MO10"></graphic>
</fig>
</p>
<p>Reactive oxygen species (ROS) are known to provide a signal that triggers cellular changes in response to altered environmental conditions [<xref ref-type="bibr" rid="CR34">34</xref>
]. The chloroplast is considered to be a major producer of ROS during abiotic stress and inhibited photosynthesis is a key factor that results in the accumulation of ROS [<xref ref-type="bibr" rid="CR35">35</xref>
]. In this study, transcript levels of unigenes associated with the ‘Photosynthesis-antenna proteins’ and ‘Photosynthesis’ GO categories, such as chlorophyll A/B binding protein (Unigene12186_All, Unigene4387_All, Unigene12405_All, Unigene18912_All, CL216.Contig1_All), chloroplast pigment-binding protein CP24 (CL9248.Contig2_All), cytochrome b6-f complex iron-sulfur subunit (Unigene1165_All), PSI reaction center subunit II (Unigene28504_All), and photosystem II oxygen-evolving enhancer protein (Unigene14933_All), decreased following the exposure treatment. Conversely, the abundance of transcripts corresponding to ROS receptor proteins or redox-sensitive transcription factors, such as Hsfs (Unigene12132_All, CL9882.Contig2_All) was up-regulated (Additional file <xref rid="MOESM8" ref-type="media">8</xref>
). We noted that expression of the unigene CL8980.Contig1_All, which encodes MBF1c, a key regulator of thermotolerance in <italic>A. thaliana</italic>
and a potential regulatory hub linking ROS signaling with pathogen and abiotic stress responses, increased in response to the cold treatment. It is known that ROS promotes cellular signaling thought a ROS-induced MAPK signaling pathway [<xref ref-type="bibr" rid="CR33">33</xref>
]. We identified several MAPKs in this study that were significantly up-regulated by the cold treatment, although their associations with ROS mediated-signaling have yet to be established.</p>
<p>Calcium serves as a second messenger and plays vital roles in a number of signaling pathways [<xref ref-type="bibr" rid="CR36">36</xref>
, <xref ref-type="bibr" rid="CR37">37</xref>
]. In plants, oscillations in calcium levels are detected by calcium sensor proteins, such as calcium-related proteins, calmodulin and calmodulin-binding protein (CBP), calcium-dependent protein kinases (CDPKs) and calcium-regulated phosphatases, which relay signals and trigger downstream responses [<xref ref-type="bibr" rid="CR38">38</xref>
, <xref ref-type="bibr" rid="CR39">39</xref>
]. In this study, we identified four genes (CL9078.Contig2_All, CL4696.Contig2_All, Unigene5810_All, CL653.Contig20_All) encoding CDPKs, one gene (Unigene22833_All) encoding a CBP and four genes (Unigene25450_All, Unigene26697_All, Unigene14626_All, and CL3238.Contig1_All) encoding CBL-interacting protein kinases (CIPKs), all of which were up-regulated by the cold treatment (Additional file <xref rid="MOESM8" ref-type="media">8</xref>
).</p>
</sec>
</sec>
<sec id="Sec12" sec-type="discussion"><title>Discussion</title>
<sec id="Sec13"><title>Illumina paired-end sequencing, assembly, and functional annotation</title>
<p>In this study, we performed an RNA-seq analysis of the <italic>P. trifoliata</italic>
transcriptome using Illumina-based 2 × 90 bp paired-end reads sequencing. Illumina, a widely used NGS platform for transcriptome assembly, produces relatively shorter reads but generates higher transcriptome coverage at lower expenses compared with other platforms [<xref ref-type="bibr" rid="CR40">40</xref>
]. Approximately 6 Gb high-quality reads were obtained, which were <italic>de novo</italic>
assembled into 77,292 unigenes. The average length of the assembled unigenes was 1112 bp and 55.5 % of the unigenes were longer than 1 kb, which indicates a high quality assembly that is likely to include many full-length cDNAs. In addition, 77.3 % of the assembled unigenes had homologs in at least one of the public databases that we searched, and 70.5 % of the unigenes had a homolog that was determined with a high probability score in three species: <italic>V. vinifera</italic>
(31.8 %), <italic>R. communis</italic>
(26.5 %) and <italic>P. trichocarpa</italic>
(22.2 %). Significantly, only 274 unigenes aligned with previously reported <italic>P. trifoliata</italic>
sequences, underlining the sparse gene coverage of this species in the Nr database. Taken together, these results suggested that our <italic>P. trifoliata</italic>
EST dataset represents a valuable transcriptome resource for gene discovery and functional analysis.</p>
</sec>
<sec id="Sec14"><title>Analysis of DEGs that are responsive to cold temperatures</title>
<p>It has been well documented that plants can enhance their cold tolerance by modulating gene transcription, and examples of gene induction and repression in response to cold temperatures have been reported [<xref ref-type="bibr" rid="CR41">41</xref>
, <xref ref-type="bibr" rid="CR42">42</xref>
]. In this study, 5549 <italic>P. trifoliata</italic>
genes were differentially expressed as a consequence of the cold treatment. Of these, 462, 1631 and 2702 unigenes were up-regulated after 6, 24, and 72 h of cold treatment, respectively, while 138, 715 and 2475 unigenes were down-regulated, This predominance of gene induction during the cold stress has been found in other studies [<xref ref-type="bibr" rid="CR11">11</xref>
, <xref ref-type="bibr" rid="CR12">12</xref>
]. It is notable that the number of DEGs progressively increased over the course of the cold treatment, indicating ongoing adaptations over an extended time. More than 28 % of the <italic>P. trifoliata</italic>
DEGs had no annotated homologs in the Nr database and these may be specific to <italic>P. trifoliata</italic>
or represent cold-responsive genes with homologs that have not been identified in previous studies with other plant species.</p>
</sec>
<sec id="Sec15"><title>Phytohormone signals under low temperature</title>
<p>The phytohormones ABA, GA and ethylene are known to play key roles in a wide range of adaptive responses to abiotic stresses [<xref ref-type="bibr" rid="CR43">43</xref>
]. In this study, the transcriptome data suggested that the expression of genes involved in signal transduction pathways related to several plant hormones were modified within 6 h of cold treatment. Moreover, increasing numbers of genes related to these pathways showed changes in expression over the 72 h of the cold treatment time course, suggesting a critical role of the hormone pathways in conferring trifoliate orange response to the cold stress. ABA is perhaps the most well characterized phytohormone in terms of integrating various stress signals and modulating downstream stress responses, although its involvement in cold responses is still debated [<xref ref-type="bibr" rid="CR44">44</xref>
]. It has been suggested that plant response to the cold stress may be ABA-independent, but there is increasing evidence to the contrary [<xref ref-type="bibr" rid="CR45">45</xref>
, <xref ref-type="bibr" rid="CR46">46</xref>
]. NCED is a rate-limiting enzyme in ABA synthesis and earlier studies reported transcriptional regulation of <italic>NCED</italic>
following abiotic stress [<xref ref-type="bibr" rid="CR47">47</xref>
]. Accordingly, in our analysis, we noticed that an <italic>NCED</italic>
gene was up-regulated by approximately 3–4 folds after 24 and 72 h of cold treatment. Up-regulation of the <italic>NCED</italic>
gene implies that synthesis of ABA might be elevated in trifoliate orange under cold. In addition, we detected the increased expression of <italic>PP2C</italic>
, a key component in the ABA signaling pathway (Table <xref rid="Tab5" ref-type="table">5</xref>
). This result is consistent with studies of <italic>P. euphratica</italic>
and <italic>A. thaliana</italic>
[<xref ref-type="bibr" rid="CR18">18</xref>
, <xref ref-type="bibr" rid="CR48">48</xref>
], but differs from an equivalent analysis of on <italic>Lilium lancifolium</italic>
[<xref ref-type="bibr" rid="CR13">13</xref>
]<italic>.</italic>
Thus, whether and how ABA mediated signaling is involved in the cold responses of trifoliate orange remains to be further established.<table-wrap id="Tab5"><label>Table 5</label>
<caption><p>Hormone-related genes that were differentially expressed during the cold treatment</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th rowspan="2">Gene ID</th>
<th colspan="3">Fold change</th>
<th rowspan="2">Annotation</th>
</tr>
<tr><th>6 h vs.0 h</th>
<th>24 h vs. 0 h</th>
<th>72 h vs. 0 h</th>
</tr>
</thead>
<tbody><tr><td align="center">ABA</td>
<td></td>
<td></td>
<td></td>
<td align="center"></td>
</tr>
<tr><td align="center">Unigene19385_All</td>
<td>0.82</td>
<td>3</td>
<td>3.8</td>
<td>9-cis-epoxycarotenoid dioxygenase (<italic>NCED</italic>
) [<italic>Citrus sinensis</italic>
]</td>
</tr>
<tr><td align="center">CL927.Contig1_All</td>
<td>2.9</td>
<td>2.9</td>
<td>4</td>
<td>short chain alcohol dehydrogenase [<italic>Ricinus communis</italic>
]</td>
</tr>
<tr><td align="center">CL851.Contig21_All</td>
<td>5.7</td>
<td>4.1</td>
<td>0.51</td>
<td>ABA 8′-hydroxylase [<italic>Citrus sinensis</italic>
]</td>
</tr>
<tr><td align="center">Unigene33085_All</td>
<td>9.7</td>
<td>12.8</td>
<td>14.1</td>
<td>ABA 8′-hydroxylase [<italic>Glycine max</italic>
]</td>
</tr>
<tr><td align="center">Unigene15626_All</td>
<td>−1.1</td>
<td>2</td>
<td>4</td>
<td>ABA 8′-hydroxylase [<italic>Vitis vinifera</italic>
]</td>
</tr>
<tr><td align="center">CL6797.Contig1_All</td>
<td>11.2</td>
<td>11.6</td>
<td>9.9</td>
<td>protein phosphatase 2C 25 isoform 1 [<italic>Vitis vinifera</italic>
]</td>
</tr>
<tr><td align="center">CL3200.Contig4_All</td>
<td>7.9</td>
<td>12.2</td>
<td>10.5</td>
<td>protein phosphatase 2C 59 [<italic>Vitis vinifera</italic>
]</td>
</tr>
<tr><td align="center">CL9817.Contig1_All</td>
<td>3.4</td>
<td>4.1</td>
<td>2.7</td>
<td>protein phosphatase 2C (PP2C) [<italic>Fagus sylvatica</italic>
]</td>
</tr>
<tr><td align="center">Ethylene</td>
<td></td>
<td></td>
<td></td>
<td align="center"></td>
</tr>
<tr><td align="center">CL5314.Contig1_All</td>
<td>1.3</td>
<td>4.9</td>
<td>5</td>
<td>S-adenosylmethionine synthetase [<italic>Litchi chinensis</italic>
]</td>
</tr>
<tr><td align="center">Unigene20820_All</td>
<td>0.7</td>
<td>2.2</td>
<td>3</td>
<td>S-adenosylmethionine decarboxylase [x <italic>Citrofortunella microcarpa</italic>
]</td>
</tr>
<tr><td align="center">Unigene21338_All</td>
<td>−1</td>
<td>−1.7</td>
<td>−3.7</td>
<td>1-aminocyclopropane-1-carboxylate oxidase homolog 4-like [<italic>Vitis vinifera</italic>
]</td>
</tr>
<tr><td align="center">Unigene23036_All</td>
<td>0.4</td>
<td>1.5</td>
<td>2.6</td>
<td>probable 2-oxoglutarate/Fe(II)-dependent dioxygenase-like [<italic>Vitis vinifera</italic>
]</td>
</tr>
<tr><td align="center">Unigene28463_All</td>
<td>1.3</td>
<td>2.7</td>
<td>2.8</td>
<td>ethylene response 2 [<italic>Citrus sinensis</italic>
]</td>
</tr>
<tr><td align="center">CL9038.Contig1_All</td>
<td>11</td>
<td>10.4</td>
<td>9.6</td>
<td>ein3-binding f-box protein 4 [<italic>Populus trichocarpa</italic>
]</td>
</tr>
<tr><td align="center">Unigene5278_All</td>
<td>3.6</td>
<td>3.6</td>
<td>3.2</td>
<td>Ethylene-responsive transcription factor 1B, putative [<italic>Ricinus communis</italic>
]</td>
</tr>
<tr><td align="center">Unigene19854_All</td>
<td>2</td>
<td>3</td>
<td>3</td>
<td>Ethylene-responsive transcription factor, putative [<italic>Ricinus communis</italic>
]</td>
</tr>
<tr><td align="center">Unigene23980_All</td>
<td>1.1</td>
<td>1.3</td>
<td>2</td>
<td>AP2/ERF domain-containing transcription factor [<italic>Populus trichocarpa</italic>
]</td>
</tr>
<tr><td align="center">Unigene21687_All</td>
<td>2.8</td>
<td>3.2</td>
<td>3.6</td>
<td>ethylene responsive transcription factor 1a [<italic>Prunus salicina</italic>
]</td>
</tr>
<tr><td align="center">Unigene21465_All</td>
<td>1.1</td>
<td>2.1</td>
<td>1.9</td>
<td>AP2/ERF domain-containing transcription factor [<italic>Populus trichocarpa</italic>
]</td>
</tr>
<tr><td align="center">JA</td>
<td></td>
<td></td>
<td></td>
<td align="center"></td>
</tr>
<tr><td align="center">CL5141.Contig4_All</td>
<td>12.1</td>
<td>12.7</td>
<td>14.1</td>
<td>old-yellow-enzyme homolog [<italic>Catharanthus roseus</italic>
]</td>
</tr>
<tr><td align="center">Unigene15722_All</td>
<td>0.3</td>
<td>1.5</td>
<td>2</td>
<td>allene oxide cyclase 2 [<italic>Glycine max</italic>
]</td>
</tr>
<tr><td align="center">CL2829.Contig1_All</td>
<td>4</td>
<td>4.2</td>
<td>3.5</td>
<td>4-coumarate--CoA ligase-like 9-like [<italic>Vitis vinifera</italic>
]</td>
</tr>
<tr><td align="center">CL9506.Contig1_All</td>
<td>3.6</td>
<td>4.8</td>
<td>4.4</td>
<td>JAZ1 [<italic>Vitis rupestris</italic>
]</td>
</tr>
<tr><td align="center">Unigene17827_All</td>
<td>1.3</td>
<td>3.6</td>
<td>3.9</td>
<td>transcription factor bHLH041 [<italic>Vitis vinifera</italic>
]</td>
</tr>
<tr><td align="center">CL6722.Contig1_All</td>
<td>1.5</td>
<td>3.4</td>
<td>2.5</td>
<td>transcription factor bHLH30-like [<italic>Glycine max</italic>
]</td>
</tr>
<tr><td align="center">Unigene25537_All</td>
<td>0.9</td>
<td>1.8</td>
<td>1.9</td>
<td>transcription factor bHLH13 [<italic>Vitis vinifera</italic>
]</td>
</tr>
<tr><td align="center">Cytokinin</td>
<td></td>
<td></td>
<td></td>
<td align="center"></td>
</tr>
<tr><td align="center">Unigene25908_All</td>
<td>3.7</td>
<td>3.2</td>
<td>1.9</td>
<td>two-component response regulator ARR5-like [<italic>Vitis vinifera</italic>
]</td>
</tr>
<tr><td align="center">Unigene25909_All</td>
<td>1.5</td>
<td>2.3</td>
<td>2.2</td>
<td>two-component response regulator ARR5-like [<italic>Vitis vinifera</italic>
]</td>
</tr>
<tr><td align="center">Unigene20772_All</td>
<td>0.6</td>
<td>3.4</td>
<td>4.8</td>
<td>Glucan endo-1,3-beta-glucosidase precursor, putative [<italic>Ricinus communis</italic>
]</td>
</tr>
<tr><td align="center">CL8049.Contig1_All</td>
<td>0.6</td>
<td>1.5</td>
<td>0.8</td>
<td>glucan endo-1,3-beta-glucosidase-like protein 3-like isoform 2 [<italic>Vitis vinifera</italic>
]</td>
</tr>
<tr><td align="center">GA</td>
<td></td>
<td></td>
<td></td>
<td align="center"></td>
</tr>
<tr><td align="center">CL9460.Contig1_All</td>
<td>1.2</td>
<td>2.2</td>
<td>2.9</td>
<td>gibberellin 2-beta-dioxygenase 1-like [<italic>Vitis vinifera</italic>
]</td>
</tr>
<tr><td align="center">Unigene21338_All</td>
<td>−1</td>
<td>−1.7</td>
<td>−3.7</td>
<td>Gibberellin 20 oxidase [<italic>Medicago truncatula</italic>
]</td>
</tr>
<tr><td align="center">Unigene14999_All</td>
<td>3.3</td>
<td>4.7</td>
<td>5.7</td>
<td>GRAS family transcription factor [<italic>Populus trichocarpa</italic>
]</td>
</tr>
<tr><td align="center">Auxin</td>
<td></td>
<td></td>
<td></td>
<td align="center"></td>
</tr>
<tr><td align="center">CL7975.Contig1_All</td>
<td>4.2</td>
<td>6.1</td>
<td>6</td>
<td>cytochrome P450 83B1 [<italic>Vitis vinifera</italic>
]</td>
</tr>
<tr><td align="center">CL5737.Contig2_All</td>
<td>4.9</td>
<td>7</td>
<td>5.3</td>
<td>phenylalanine N-monooxygenase-like [<italic>Vitis vinifera</italic>
]</td>
</tr>
<tr><td align="center">Unigene22952_All</td>
<td>−0.4</td>
<td>−1</td>
<td>−1.9</td>
<td>protein AUXIN SIGNALING F-BOX 3 [<italic>Vitis vinifera</italic>
]</td>
</tr>
<tr><td align="center">CL9773.Contig2_All</td>
<td>0.9</td>
<td>4.9</td>
<td>6.2</td>
<td>IAA15 [<italic>Solanum lycopersicum</italic>
]</td>
</tr>
<tr><td align="center">Unigene20450_All</td>
<td>1.4</td>
<td>4.1</td>
<td>6</td>
<td>Auxin-responsive protein IAA1, putative [<italic>Ricinus communis</italic>
]</td>
</tr>
<tr><td align="center">Brassinosteroid</td>
<td></td>
<td></td>
<td></td>
<td align="center"></td>
</tr>
<tr><td align="center">Unigene19738_All</td>
<td>−0.1</td>
<td>4.6</td>
<td>4.6</td>
<td>cytochrome P450, putative [<italic>Ricinus communis</italic>
]</td>
</tr>
<tr><td align="center">Unigene33085_All</td>
<td>9.7</td>
<td>12.8</td>
<td>14.1</td>
<td>cytochrome P450, putative [<italic>Ricinus communis</italic>
]</td>
</tr>
<tr><td align="center">CL1502.Contig3_All</td>
<td>3.9</td>
<td>4.1</td>
<td>3.4</td>
<td>wall-associated receptor kinase-like 1-like [<italic>Vitis vinifera</italic>
]</td>
</tr>
<tr><td align="center">CL7113.Contig3_All</td>
<td>3.9</td>
<td>4.6</td>
<td>3.8</td>
<td>wall-associated receptor kinase-like 9-like [<italic>Vitis vinifera</italic>
]</td>
</tr>
<tr><td align="center">CL8664.Contig1_All</td>
<td>15.2</td>
<td>14.7</td>
<td>13.5</td>
<td>xyloglucan endotransglycosylase [<italic>Cucumis sativus</italic>
]</td>
</tr>
</tbody>
</table>
</table-wrap>
</p>
<p>Ethylene has been reported to participate in a wide range of cellular and developmental processes, and in both abiotic and biotic stress responses [<xref ref-type="bibr" rid="CR49">49</xref>
, <xref ref-type="bibr" rid="CR50">50</xref>
]. The molecular mechanisms underlying perception and transduction of the ethylene signal and activation of hundreds of relevant genes have been extensively studied [<xref ref-type="bibr" rid="CR51">51</xref>
]. Ethylene signaling has been implicated in the cold stress response [<xref ref-type="bibr" rid="CR52">52</xref>
–<xref ref-type="bibr" rid="CR55">55</xref>
], but different studies have attributed it with either positive or negative effects. For example, the positive role of ethylene in cold stress was demonstrated in rice [<xref ref-type="bibr" rid="CR52">52</xref>
], tomato [<xref ref-type="bibr" rid="CR53">53</xref>
], winter rye [<xref ref-type="bibr" rid="CR54">54</xref>
] and <italic>A. thaliana</italic>
[<xref ref-type="bibr" rid="CR56">56</xref>
], but Shi et al. [<xref ref-type="bibr" rid="CR57">57</xref>
] reported that ethylene signaling negatively influences freezing tolerance by repressing the expression of <italic>CBF</italic>
and type-A <italic>ARR</italic>
genes. The discrepancy may reflect differences in the species studied or to the growth conditions of the tissues used for the study [<xref ref-type="bibr" rid="CR56">56</xref>
]. In our study, we observed that the expression an ACC synthase gene (<italic>ACS</italic>
, CL4785.Contig1_All) was repressed following cold stress. ACS is involved in ethylene synthesis and so its down-regulation may lead to reduced ethylene production (Fig. <xref rid="Fig10" ref-type="fig">10</xref>
, Table <xref rid="Tab5" ref-type="table">5</xref>
). This result is consistent with an earlier study [<xref ref-type="bibr" rid="CR52">52</xref>
] and prompts us to propose that ethylene may be possibly a negative regulator of cold responses in trifoliate orange.</p>
<p>In addition to ABA and ethylene, other hormones may also play a role in cold stress response. The role of the growth-related hormone GA in the response to abiotic stress is becoming increasingly elucidated [<xref ref-type="bibr" rid="CR58">58</xref>
]. For example, DELLA proteins have been shown to suppress GA signaling and can promote GA degradation via the ubiquitin-proteasome [<xref ref-type="bibr" rid="CR59">59</xref>
]. Accordingly, exposure of <italic>A. thaliana</italic>
seedlings to cold stress promoted DELLA accumulation and resulted in DELLA-mediated growth restriction [<xref ref-type="bibr" rid="CR60">60</xref>
]. We observed that the expression of a unigene (Unigene17855_All) encoding a DELLA protein was up-regulated, suggesting that degradation of GA may be expedited by cold stress. In this regard, the cold-exposed plants may slow their growth in order to conserve energy for adapting to the adverse environments. Taken together, these results suggest that diverse hormonal signals are involved in responses and the expression modulation of large numbers of downstream genes.</p>
</sec>
<sec id="Sec16"><title>ROS-mediated signals associated with cold response</title>
<p>ROS produced during abiotic stresses are known to be toxic to cellular functions when they accumulate to high levels; however, there is increasing evidence that at low levels ROS can mediate an array of signals [<xref ref-type="bibr" rid="CR34">34</xref>
]. Inhibition of photosynthesis is a significant cause of ROS production [<xref ref-type="bibr" rid="CR35">35</xref>
]. In this work, ‘photosynthesis-antenna proteins’ and ‘photosynthesis’ are two major GO terms that were associated with substantial change in transcript levels in response to cold stress. Interestingly, transcript levels of genes involved in these pathways were decreased after exposure to the cold treatment indicating the photosynthetic capacity of <italic>P. trifoliata</italic>
was inhibited, in turn leading to ROS accumulation. The ROS signal is thought to be sensed by yet unidentified receptor proteins or redox-sensitive TFs, such as Hsfs [<xref ref-type="bibr" rid="CR61">61</xref>
]. In our study, two Hsf unigenes (Unigene12132_All, CL9882.Contig2_All) were differentially expressed in response to cold, but their potential roles in ROS signal sensing remain to be investigated. In addition, ROS signals are transduced to downstream targets by either a MAPK cascade or the action of TFs, such as zinc finger proteins (ZFPs), WRKYs and multiprotein bridging factor 1c (MBF1c) [<xref ref-type="bibr" rid="CR62">62</xref>
]. We observed in our study that the transcript levels of a MAPK (CL8004.Contig1_All) and two MAPKKKs (Unigene25596_All, CL3174.Contig_All) were up-regulated by the cold treatment, as were those corresponding to unigenes encoding several ZFP and WRKY proteins. We hypothesize that different components of the ROS signal transduction pathway are involved in mediating the cold stress response, and our dataset has revealed some potential candidates.</p>
</sec>
<sec id="Sec17"><title>Transcription factors involved in cold stress response</title>
<p>Numerous families of TFs are known to orchestrate the signals that are transduced when plants are challenged with various abiotic stresses, including the MYB, AP2/EREBP, WRKY, NAC, MYC and bHLH families, and individual members have been revealed to promote or suppress abiotic stress responses [<xref ref-type="bibr" rid="CR63">63</xref>
, <xref ref-type="bibr" rid="CR64">64</xref>
]. In the current study, 54 putative TFs were found to be induced either at a given time point or during the whole course of the cold treatment. The largest group of the cold-inducible TFs belonged to the AP2/ERF family and was composed of 26 members. Of these, three unigenes (Unigene5199, Unigene14645 and Unigene36912) were annotated as <italic>CBF/DREB</italic>
genes, which have been shown to play important roles in cold acclimation leading to freezing tolerance [<xref ref-type="bibr" rid="CR5">5</xref>
]. Our results therefore suggest that the CBF pathway is conserved in trifoliate orange responses to the cold stress. We also noticed that the cold stress induced the expression of TFs in the NAC (CL6682.Contig4, Unigene7546), bHLH (CL6722.Contig1, Unigene22237, Unigene8451), WRKY (Unigene5506, Unigene178), and ERF (Unigene20110) families. Earlier studies have shown that members of these families function in cold tolerance [<xref ref-type="bibr" rid="CR65">65</xref>
], but whether they work, independently or synergistically, to promote cold tolerance of trifoliate orange remains to be determined.</p>
</sec>
<sec id="Sec18"><title>Polyamine biosynthetic genes implicate in cold response</title>
<p>Upon exposure to various abiotic stresses, plants accumulate various metabolites that are thus considered as potential indices indicating the magnitude of stress tolerance [<xref ref-type="bibr" rid="CR66">66</xref>
]. As an example, polyamines (PAs) are viewed as important stress-related metabolites that can function to reduce damage associated with adverse environmental cues. The polyamines, being polycationic, stabilize cellular membranes through binding to anions or act as critical protective metabolites by mitigating osmotic stress [<xref ref-type="bibr" rid="CR67">67</xref>
]. A number of earlier studies reported that PA biosynthetic genes are up-regulated by cold conditions; and overexpression of these biosynthetic genes can confer cold tolerance [<xref ref-type="bibr" rid="CR24">24</xref>
]. In this study, we detected a change in the mRNA levels of annotated polyamine biosynthetic genes in response to cold stress, including <italic>ADC</italic>
(Unigene25843_All) and <italic>SAMDC</italic>
(Unigene20820_All), which encode two rate-limiting enzymes in the synthesis of putrescine, spermidine and spermine, respectively. The up-regulation of <italic>ADC</italic>
gene coincides with the increase in the putrescine level in the cold-treated samples. However, <italic>SAMDC</italic>
gene induction was not accompanied by a concurrent elevation of spermidine and spermine. These data suggest that synthesis of putrescine might be promoted, which may provide the plants to combat the cold stress.</p>
</sec>
</sec>
<sec id="Sec19" sec-type="conclusion"><title>Conclusions</title>
<p>In this study, a <italic>P. trifoliata</italic>
transcriptome dataset comprising 77,292 predicted transcripts was generated by high-throughput sequencing and dynamic changes in gene expression under cold treatments were observed. A large number of cold-responsive genes were revealed, including those encoding TFs, ROS and hormone signaling elements, and enzymes associated with the synthesis of protective metabolites. The transcriptome and digital expression profiling of <italic>P. trifoliata</italic>
provide a valuable resource for functional evaluation of the cold-responsive genes, which are assumed to hold great potential for genetic engineering of cold tolerance.</p>
</sec>
<sec id="Sec20" sec-type="materials|methods"><title>Methods</title>
<sec id="Sec21"><title>Plant materials and stress treatments</title>
<p>Three-month-old trifoliate orange (<italic>Poncirus trifoliata</italic>
(L.) Raf.) seedlings were collected from an experimental nursery at Huazhong Agricultural University. The plants were carefully and extensively washed and placed in glass beakers of tap water and kept for 5 d at 25 °C in a growth chamber. In order to promote the expression of stress-responsive genes, the seedlings were treated with different stresses, including cold, high salinity, and drought. For the cold treatment, the plants were placed at 4 °C and the leaves were harvested immediately (0 h time point) or after 6, 24 or 72 h. The salinity treatment was applied by placing the seedlings in 200 mM NaCl solution for same time periods as the cold treatment. For the drought treatment, the seedlings were placed on filter paper and dehydrated at ambient room temperature and the whole plants were collected at 0, 0.5, 3, and 6 h. For each time point, at least 20 seedlings were used, and the collected samples were immediately frozen in liquid nitrogen after sampling and stored at −80 °C until use for further use. In addition, a new set of seedlings were treated at 4 °C for 0, 6, 24 or 72 h; the leaves were sampled for analysis of expression patterns of several genes involved in hormone metabolism and signal transduction.</p>
</sec>
<sec id="Sec22"><title>RNA extraction</title>
<p>Total RNA was extracted using an RNAiso Plus kit (TaKaRa, Dalian, China) according to the manufacturer's instructions. The concentration and quality of total RNA was evaluated with a NanoDrop<sup>TM</sup>
2000 UV–vis Spectrophotometer (Thermo Scientific, Waltham, MA, USA). For transcriptome sequencing, equal amount of total RNA from the stress-treated samples were pooled.</p>
</sec>
<sec id="Sec23"><title>cDNA library construction and RNA-sequencing</title>
<p>The pooled cDNA library was constructed using an mRNA-Seq assay for paired-end transcriptome sequencing. In brief, poly(A) mRNA was enriched from 20 μg of total RNA using oligo(dT) magnetic beads, and sheared into short fragments (200–700 bp) with fragmentation buffer. The short mRNA fragments were synthesized into double-strand cDNAs using random hexamer primers, and purified with a QiaQuick PCR Purification Kit (Qiagen, CA, USA). Then cDNAs were washed with EB buffer (10 mM Tris–HCl, pH 8.5) for end-repairing and poly(A) addition. Illumina paired-end sequencing adapters were then ligated to the ends of the 3′-adenylated cDNA fragments. Agarose gel electrophoresis was carried out to select fragments of a suitable size that was enriched by PCR amplification to construct the cDNA library. The library was sequenced by BGI (Shenzhen, China) on an Illumina HiSeq<sup>™</sup>
2000 platform.</p>
</sec>
<sec id="Sec24"><title>Bioinformatics analysis of the transcriptome</title>
<p>Raw reads generated by the Hiseq<sup>TM</sup>
2000 were filtered to remove low quality reads (reads containing more than 50 % bases with Q-value ≤20), adaptor-containing reads, and reads containing more than 5 % ambiguous nucleotides. After the preprocessing, >5 Gb filtered short reads, termed ‘clean reads’, were obtained. These were then <italic>de novo</italic>
assembled using Trinity to generate a collection of non-redundant unigenes. To this end, the reads with overlapping sequences were merged to generate longer contiguous sequences (contigs), and the reads were then mapped back to the contigs. The relation and distance among these contigs were determined based on paired-end reads, which enabled the detection of contigs from the same transcripts as well as the distances between these contigs. The contigs were then assembled using Trinity to identify the sequences without end extension, which were defined as unigenes. TIGR Gene Indices clustering tools (TGICL) [<xref ref-type="bibr" rid="CR68">68</xref>
] was used to eliminate redundancy and further assemble the unigenes. The unigenes displaying >70 % sequence similarity were grouped into a cluster and the prefix CL was assigned to each cluster, while those from the other group was considered as singletons with the prefix unigene.</p>
<p>All non-redundant transcripts were subjected to BLASTx alignment (<italic>E</italic>
-value < 10<sup>−5</sup>
) against public protein databases, including Nr, Swiss-Prot, KEGG and COG. The assembled unigenes were annotated using Blast2GO to generate GO terms according to molecular function, biological process and cellular component, based on the best BLASTx hit from the NR database. Unigenes without an annotation were analyzed by ESTScan to predict coding regions and sequence direction. Metabolic pathways of the unigenes were determined according to Kyoto Encyclopedia of Genes and Genomes (KEGG) database [<xref ref-type="bibr" rid="CR69">69</xref>
].</p>
</sec>
<sec id="Sec25"><title>Analysis and mapping of DEG reads</title>
<p>Four cDNA libraries from trifoliate orange treated with cold for 0, 6, 24, and 72 h were sequenced and the raw data were processed as described above. The short cleaned reads were mapped to the assembled unigenes using SOAPaligner/soap2 software. The read counts, indicating unigene expression, in different libraries were normalized as RPKM (reads per kilobase per million mapped reads) as described by Mortazavi et al. [<xref ref-type="bibr" rid="CR70">70</xref>
]. Differentially expressed genes (DEGs) were identified using the R package. The DEGs were also annotated using the GO database, and the numbers of DEGs in each GO term were calculated. KEGG pathway analysis of the DEGs was also performed to identify the associated biochemical and signal transduction pathways. The data of this study have been deposited in NCBI’s Gene Expression Omnibus (GEO) under an accession number GSE67439 and are accessible through GEO Series (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE67439">http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE67439</ext-link>
).</p>
</sec>
<sec id="Sec26"><title>Quantitative real-time RT-PCR analysis</title>
<p>Quantitative real-time RT-PCR (qPCR) was performed with a Roche LightCycler 480 Real-Time System (Roche, Switzerland) to examine expression patterns of the selected unigenes. Each reaction contained 5 μl of 2x QuantiFast SYBR Green PCR Master Mix (Qiagen, Hilden, Germany), 1 μl of cDNA, and 1 μM of gene-specific primers in a final volume of 10 μl. The PCR reactions were performed under the following conditions: 95 °C for 5 min, followed by 45 cycles of 95 °C for 15 s, 58 °C for 20 s and 72 °C for 20 s. Relative expression levels of each gene were calculated using the 2<sup>-ΔΔCt</sup>
algorithm by normalizing to expression of expression of the <italic>P. trifoliata Actin</italic>
gene, which was used as an internal control. Four technical replicates were used for each sample and the data are shown as means ± standard errors (SE) (n = 3). The source of variation resulted from the technical errors, such as operational approach, equipment and reagent. The primer sequences used for qPCR are listed in Additional file <xref rid="MOESM9" ref-type="media">9</xref>
.</p>
</sec>
<sec id="Sec27"><title>Quantification of free PAs by high-performance liquid chromatography (HPLC)</title>
<p>Free polyamines were measured as described by Fu et al. [<xref ref-type="bibr" rid="CR71">71</xref>
] with some modifications. Briefly, leaf samples were ground into powder in liquid nitrogen with a pestle and mortar and ~0.1 g of powder was homogenized in 1 ml 5 % pre-cooled perchloric acid (PCA) containing 0.5 g/L dithiothreitol (DTT), and kept on ice for 30 min, followed by centrifugation at 13,400 × g for 15 min at 4 °C. This extraction was then repeated, and the supernatants were combined, and 1 ml was mixed with 50 μl of 1 M hexamethylene diamine, 1 ml of 2 M NaOH and 10 μl of benzoyl chloride. The mixture was vortexed for 30 s and placed in a water bath for 30 min. Subsequently, 2 ml of saturated NaCl and 2 ml of diethyl ether were added to the mixture and mixed well by vortexing before centrifugation at 6900 × g for 5 min. One ml of the supernatant was dried in a concentrator (Labconco, USA) and dissolved in 500 μl of methanol (Fisher, USA). The solution was filtered through 0.22 μm filter membrane and 20 μl was analyzed using an Agilent 1200 HPLC system (USA). The polyamines were determined in triplicate.</p>
</sec>
</sec>
</body>
<back><app-group><app id="App1"><sec id="Sec28"><title>Additional files</title>
<p><media position="anchor" xlink:href="12864_2015_1629_MOESM1_ESM.xlsx" id="MOESM1"><label>Additional file 1:</label>
<caption><p><bold>KEGG mapping of the</bold>
<bold><italic>P. trifoliata</italic>
</bold>
<bold> transcriptome (XLSX).</bold>
</p>
</caption>
</media>
<media position="anchor" xlink:href="12864_2015_1629_MOESM2_ESM.xlsx" id="MOESM2"><label>Additional file 2:</label>
<caption><p><bold>Differentially expressed genes (DEGs) between the cold treatment (6 h, 24 h and 72 h) and control (0 h, untreated plants).</bold>
RPKM: reads per kilobase per million mapped reads. FDR: false discovery rate (XLSX).</p>
</caption>
</media>
<media position="anchor" xlink:href="12864_2015_1629_MOESM3_ESM.pdf" id="MOESM3"><label>Additional file 3:</label>
<caption><p><bold>Saturation evaluation of different gene expressions (A) and output of cleaned reads (B) (PDF).</bold>
</p>
</caption>
</media>
<media position="anchor" xlink:href="12864_2015_1629_MOESM4_ESM.pdf" id="MOESM4"><label>Additional file 4:</label>
<caption><p><bold>Scatter diagram showing the distribution of up-regulated differentially expressed genes (DEGs) in red and down-regulated DEGs in green (PDF).</bold>
</p>
</caption>
</media>
<media position="anchor" xlink:href="12864_2015_1629_MOESM5_ESM.xlsx" id="MOESM5"><label>Additional file 5:</label>
<caption><p><bold>Significantly enriched GO terms amongst the differentially expressed genes (XLSX).</bold>
</p>
</caption>
</media>
<media position="anchor" xlink:href="12864_2015_1629_MOESM6_ESM.xlsx" id="MOESM6"><label>Additional file 6:</label>
<caption><p><bold>List of 60 differentially expressed genes (DEGs) encoding transcription factors (XLSX).</bold>
</p>
</caption>
</media>
<media position="anchor" xlink:href="12864_2015_1629_MOESM7_ESM.xlsx" id="MOESM7"><label>Additional file 7:</label>
<caption><p><bold>Differentially expressed genes (DEGs) encoding protein kinases and DEGs involved in phytohormone metabolism and signaling transduction pathway (XLSX).</bold>
</p>
</caption>
</media>
<media position="anchor" xlink:href="12864_2015_1629_MOESM8_ESM.pdf" id="MOESM8"><label>Additional file 8:</label>
<caption><p><bold>A list of some of the important differentially expressed genes (DEGs) involved in photosynthesis, ROS and Ca</bold>
<sup><bold>2+</bold>
</sup>
<bold>-mediated signal transduction (PDF).</bold>
</p>
</caption>
</media>
<media position="anchor" xlink:href="12864_2015_1629_MOESM9_ESM.pdf" id="MOESM9"><label>Additional file 9:</label>
<caption><p><bold>Primers used for qPCR analysis (PDF).</bold>
</p>
</caption>
</media>
</p>
</sec>
</app>
</app-group>
<glossary><title>Abbreviations</title>
<def-list><def-item><term>ACS</term>
<def><p>Aminocyclopropane carboxylic acid (ACC) synthase</p>
</def>
</def-item>
<def-item><term>ADC</term>
<def><p>Arginine decarboxylase</p>
</def>
</def-item>
<def-item><term>CA</term>
<def><p>Cold acclimation</p>
</def>
</def-item>
<def-item><term>CBF</term>
<def><p>C-repeat (CRT)-binding factor</p>
</def>
</def-item>
<def-item><term>CDSs</term>
<def><p>Coding sequences</p>
</def>
</def-item>
<def-item><term>COG</term>
<def><p>Clusters of orthologous groups</p>
</def>
</def-item>
<def-item><term>COR</term>
<def><p>Cold-regulated</p>
</def>
</def-item>
<def-item><term>DEGs</term>
<def><p>Differentially expressed genes</p>
</def>
</def-item>
<def-item><term>DREB1</term>
<def><p>Dehydration-responsive element-binding factors 1</p>
</def>
</def-item>
<def-item><term>FPKM</term>
<def><p>Fragments per kb per million fragments</p>
</def>
</def-item>
<def-item><term>GO</term>
<def><p>Gene ontology</p>
</def>
</def-item>
<def-item><term>HPLC</term>
<def><p>High-performance liquid chromatography</p>
</def>
</def-item>
<def-item><term>KEGG</term>
<def><p>Kyoto encyclopedia of genes and genomes</p>
</def>
</def-item>
<def-item><term>MAPK</term>
<def><p>Mitogen-activated protein kinase</p>
</def>
</def-item>
<def-item><term>qPCR</term>
<def><p>Quantitative real-time RT-PCR</p>
</def>
</def-item>
<def-item><term>ODC</term>
<def><p>Ornithine decarboxylase</p>
</def>
</def-item>
<def-item><term>ROS</term>
<def><p>Reactive oxygen species</p>
</def>
</def-item>
<def-item><term>SAMDC</term>
<def><p><italic>S</italic>
-adenosylmethionine decarboxylase</p>
</def>
</def-item>
<def-item><term>SPDS</term>
<def><p>Spermidine synthase</p>
</def>
</def-item>
<def-item><term>SPMS</term>
<def><p>Spermine synthase</p>
</def>
</def-item>
<def-item><term>TF</term>
<def><p>Transcription factor</p>
</def>
</def-item>
</def-list>
</glossary>
<fn-group><fn><p><bold>Competing interests</bold>
</p>
<p>The authors declare that they have no competing interests.</p>
</fn>
<fn><p><bold>Authors’ contributions</bold>
</p>
<p>Conceived and designed the experiments: JHL. Performed the experiments and analyzed the data: MW, XNZ. Wrote the manuscript: MW, JHL. Read and approved the final manuscript: JHL. All authors read and approved the final manuscript.</p>
</fn>
</fn-group>
<ack><title>Acknowledgements</title>
<p>This work was supported by the National High Technology Research and Development Program (863 Program) of China (2011AA100205), the Ministry of Agriculture, and the National Natural Science Foundation of China. We thank PlantScribe (<ext-link ext-link-type="uri" xlink:href="http://www.plantscribe.com">www.plantscribe.com</ext-link>
) for editing this manuscript.</p>
</ack>
<ref-list id="Bib1"><title>References</title>
<ref id="CR1"><label>1.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chinnusamy</surname>
<given-names>V</given-names>
</name>
<name><surname>Zhu</surname>
<given-names>JH</given-names>
</name>
<name><surname>Zhu</surname>
<given-names>JK</given-names>
</name>
</person-group>
<article-title>Cold stress regulation of gene expression in plants</article-title>
<source>Trends Plant Sci</source>
<year>2007</year>
<volume>12</volume>
<fpage>444</fpage>
<lpage>51</lpage>
<pub-id pub-id-type="doi">10.1016/j.tplants.2007.07.002</pub-id>
<pub-id pub-id-type="pmid">17855156</pub-id>
</element-citation>
</ref>
<ref id="CR2"><label>2.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Michael</surname>
<given-names>MF</given-names>
</name>
</person-group>
<article-title>Plant cold acclimation: freezing tolerance genes and regulatory mechanisms</article-title>
<source>Annu Rev Plant Biol</source>
<year>1999</year>
<volume>50</volume>
<fpage>571</fpage>
<lpage>99</lpage>
<pub-id pub-id-type="doi">10.1146/annurev.arplant.50.1.571</pub-id>
</element-citation>
</ref>
<ref id="CR3"><label>3.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhu</surname>
<given-names>JH</given-names>
</name>
<name><surname>Dong</surname>
<given-names>CH</given-names>
</name>
<name><surname>Zhu</surname>
<given-names>JK</given-names>
</name>
</person-group>
<article-title>Interplay between cold-responsive gene regulation, metabolism and RNA processing during plant cold acclimation</article-title>
<source>Curr Opin Plant Biol</source>
<year>2007</year>
<volume>10</volume>
<fpage>290</fpage>
<lpage>5</lpage>
<pub-id pub-id-type="doi">10.1016/j.pbi.2007.04.010</pub-id>
<pub-id pub-id-type="pmid">17468037</pub-id>
</element-citation>
</ref>
<ref id="CR4"><label>4.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Medina</surname>
<given-names>J</given-names>
</name>
<name><surname>Catalá</surname>
<given-names>R</given-names>
</name>
<name><surname>Salinas</surname>
<given-names>J</given-names>
</name>
</person-group>
<article-title>The CBFs: three Arabidopsis transcription factors to cold acclimate</article-title>
<source>Plant Sci</source>
<year>2011</year>
<volume>180</volume>
<fpage>3</fpage>
<lpage>11</lpage>
<pub-id pub-id-type="doi">10.1016/j.plantsci.2010.06.019</pub-id>
<pub-id pub-id-type="pmid">21421341</pub-id>
</element-citation>
</ref>
<ref id="CR5"><label>5.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nakashima</surname>
<given-names>K</given-names>
</name>
<name><surname>Ito</surname>
<given-names>Y</given-names>
</name>
<name><surname>Yamaguchi-Shinozaki</surname>
<given-names>K</given-names>
</name>
</person-group>
<article-title>Transcriptional regulatory networks in response to abiotic stresses in Arabidopsis and grasses</article-title>
<source>Plant Physiol</source>
<year>2009</year>
<volume>149</volume>
<fpage>88</fpage>
<lpage>95</lpage>
<pub-id pub-id-type="doi">10.1104/pp.108.129791</pub-id>
<pub-id pub-id-type="pmid">19126699</pub-id>
</element-citation>
</ref>
<ref id="CR6"><label>6.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chinnusamy</surname>
<given-names>V</given-names>
</name>
<name><surname>Ohta</surname>
<given-names>M</given-names>
</name>
<name><surname>Kanrar</surname>
<given-names>S</given-names>
</name>
<name><surname>Lee</surname>
<given-names>B</given-names>
</name>
<name><surname>Hong</surname>
<given-names>XH</given-names>
</name>
<name><surname>Agarwal</surname>
<given-names>M</given-names>
</name>
<etal></etal>
</person-group>
<article-title>ICE1: a regulator of cold-induced transcriptome and freezing tolerance in <italic>Arabidopsis</italic>
</article-title>
<source>Gene Dev</source>
<year>2003</year>
<volume>17</volume>
<fpage>1043</fpage>
<lpage>54</lpage>
<pub-id pub-id-type="doi">10.1101/gad.1077503</pub-id>
<pub-id pub-id-type="pmid">12672693</pub-id>
</element-citation>
</ref>
<ref id="CR7"><label>7.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dong</surname>
<given-names>CH</given-names>
</name>
<name><surname>Agarwal</surname>
<given-names>M</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>YY</given-names>
</name>
<name><surname>Xie</surname>
<given-names>Q</given-names>
</name>
<name><surname>Zhu</surname>
<given-names>JK</given-names>
</name>
</person-group>
<article-title>The negative regulator of plant cold responses, HOS1, is a RING E3 ligase that mediates the ubiquitination and degradation of ICE1</article-title>
<source>Proc Natl Acad Sci U S A</source>
<year>2006</year>
<volume>103</volume>
<fpage>8281</fpage>
<lpage>6</lpage>
<pub-id pub-id-type="doi">10.1073/pnas.0602874103</pub-id>
<pub-id pub-id-type="pmid">16702557</pub-id>
</element-citation>
</ref>
<ref id="CR8"><label>8.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Agarwal</surname>
<given-names>M</given-names>
</name>
<name><surname>Hao</surname>
<given-names>YJ</given-names>
</name>
<name><surname>Kapoor</surname>
<given-names>A</given-names>
</name>
<name><surname>Dong</surname>
<given-names>CH</given-names>
</name>
<name><surname>Fujii</surname>
<given-names>H</given-names>
</name>
<name><surname>Zheng</surname>
<given-names>XY</given-names>
</name>
<etal></etal>
</person-group>
<article-title>A R2R3 type MYB transcription factor is involved in the cold regulation of CBF genes and in acquired freezing tolerance</article-title>
<source>J Biol Chem</source>
<year>2006</year>
<volume>281</volume>
<fpage>37636</fpage>
<lpage>45</lpage>
<pub-id pub-id-type="doi">10.1074/jbc.M605895200</pub-id>
<pub-id pub-id-type="pmid">17015446</pub-id>
</element-citation>
</ref>
<ref id="CR9"><label>9.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Matsui</surname>
<given-names>A</given-names>
</name>
<name><surname>Ishida</surname>
<given-names>J</given-names>
</name>
<name><surname>Morosawa</surname>
<given-names>T</given-names>
</name>
<name><surname>Mochizuki</surname>
<given-names>Y</given-names>
</name>
<name><surname>Kaminuma</surname>
<given-names>E</given-names>
</name>
<name><surname>Endo</surname>
<given-names>TA</given-names>
</name>
<etal></etal>
</person-group>
<article-title><italic>Arabidopsis</italic>
transcriptome analysis under drought, cold, high-salinity and ABA treatment conditions using a tiling array</article-title>
<source>Plant Cell Physiol</source>
<year>2008</year>
<volume>49</volume>
<fpage>1135</fpage>
<lpage>49</lpage>
<pub-id pub-id-type="doi">10.1093/pcp/pcn101</pub-id>
<pub-id pub-id-type="pmid">18625610</pub-id>
</element-citation>
</ref>
<ref id="CR10"><label>10.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Winfield</surname>
<given-names>MO</given-names>
</name>
<name><surname>Lu</surname>
<given-names>C</given-names>
</name>
<name><surname>Wilson</surname>
<given-names>ID</given-names>
</name>
<name><surname>Coghill</surname>
<given-names>JA</given-names>
</name>
<name><surname>Edwards</surname>
<given-names>KJ</given-names>
</name>
</person-group>
<article-title>Plant responses to cold: transcriptome analysis of wheat</article-title>
<source>Plant Biotechnol J</source>
<year>2010</year>
<volume>8</volume>
<fpage>749</fpage>
<lpage>71</lpage>
<pub-id pub-id-type="doi">10.1111/j.1467-7652.2010.00536.x</pub-id>
<pub-id pub-id-type="pmid">20561247</pub-id>
</element-citation>
</ref>
<ref id="CR11"><label>11.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tian</surname>
<given-names>DQ</given-names>
</name>
<name><surname>Pan</surname>
<given-names>XY</given-names>
</name>
<name><surname>Yu</surname>
<given-names>YM</given-names>
</name>
<name><surname>Wang</surname>
<given-names>WY</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>F</given-names>
</name>
<name><surname>Ge</surname>
<given-names>YY</given-names>
</name>
<etal></etal>
</person-group>
<article-title>De novo characterization of the <italic>Anthurium</italic>
transcriptome and analysis of its digital gene expression under cold stress</article-title>
<source>BMC Genomics</source>
<year>2013</year>
<volume>14</volume>
<fpage>827</fpage>
<pub-id pub-id-type="doi">10.1186/1471-2164-14-827</pub-id>
<pub-id pub-id-type="pmid">24267953</pub-id>
</element-citation>
</ref>
<ref id="CR12"><label>12.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname>
<given-names>XC</given-names>
</name>
<name><surname>Zhao</surname>
<given-names>QY</given-names>
</name>
<name><surname>Ma</surname>
<given-names>CL</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>ZH</given-names>
</name>
<name><surname>Cao</surname>
<given-names>HL</given-names>
</name>
<name><surname>Kong</surname>
<given-names>YM</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Global transcriptome profiles of <italic>Camellia sinensis</italic>
during cold acclimation</article-title>
<source>BMC Genomics</source>
<year>2013</year>
<volume>14</volume>
<fpage>415</fpage>
<pub-id pub-id-type="doi">10.1186/1471-2164-14-415</pub-id>
<pub-id pub-id-type="pmid">23799877</pub-id>
</element-citation>
</ref>
<ref id="CR13"><label>13.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname>
<given-names>JM</given-names>
</name>
<name><surname>Yang</surname>
<given-names>Y</given-names>
</name>
<name><surname>Liu</surname>
<given-names>XH</given-names>
</name>
<name><surname>Huang</surname>
<given-names>J</given-names>
</name>
<name><surname>Wang</surname>
<given-names>Q</given-names>
</name>
<name><surname>Gu</surname>
<given-names>JH</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Transcriptome profiling of the cold response and signaling pathways in <italic>Lilium lancifolium</italic>
</article-title>
<source>BMC Genomics</source>
<year>2014</year>
<volume>15</volume>
<fpage>203</fpage>
<pub-id pub-id-type="doi">10.1186/1471-2164-15-203</pub-id>
<pub-id pub-id-type="pmid">24636716</pub-id>
</element-citation>
</ref>
<ref id="CR14"><label>14.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Fowler</surname>
<given-names>S</given-names>
</name>
<name><surname>Thomashow</surname>
<given-names>MF</given-names>
</name>
</person-group>
<article-title><italic>Arabidopsis</italic>
transcriptome profiling indicates that multiple regulatory pathways are activated during cold acclimation in addition to the CBF cold response pathway</article-title>
<source>Plant Cell</source>
<year>2002</year>
<volume>14</volume>
<fpage>1675</fpage>
<lpage>90</lpage>
<pub-id pub-id-type="doi">10.1105/tpc.003483</pub-id>
<pub-id pub-id-type="pmid">12172015</pub-id>
</element-citation>
</ref>
<ref id="CR15"><label>15.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Martin</surname>
<given-names>JA</given-names>
</name>
<name><surname>Wang</surname>
<given-names>Z</given-names>
</name>
</person-group>
<article-title>Next-generation transcriptome assembly</article-title>
<source>Nat Rev Genet</source>
<year>2011</year>
<volume>12</volume>
<fpage>671</fpage>
<lpage>82</lpage>
<pub-id pub-id-type="doi">10.1038/nrg3068</pub-id>
<pub-id pub-id-type="pmid">21897427</pub-id>
</element-citation>
</ref>
<ref id="CR16"><label>16.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Huang</surname>
<given-names>JZ</given-names>
</name>
<name><surname>Lu</surname>
<given-names>X</given-names>
</name>
<name><surname>Yan</surname>
<given-names>H</given-names>
</name>
<name><surname>Chen</surname>
<given-names>SY</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>WK</given-names>
</name>
<name><surname>Huang</surname>
<given-names>RF</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Transcriptome characterization and sequencing-based identification of salt-responsive genes in <italic>Millettia pinnata</italic>
, a semi-mangrove plant</article-title>
<source>DNA Res</source>
<year>2012</year>
<volume>19</volume>
<fpage>195</fpage>
<lpage>207</lpage>
<pub-id pub-id-type="doi">10.1093/dnares/dss004</pub-id>
<pub-id pub-id-type="pmid">22351699</pub-id>
</element-citation>
</ref>
<ref id="CR17"><label>17.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Shi</surname>
<given-names>Y</given-names>
</name>
<name><surname>Yan</surname>
<given-names>X</given-names>
</name>
<name><surname>Zhao</surname>
<given-names>PS</given-names>
</name>
<name><surname>Yin</surname>
<given-names>HX</given-names>
</name>
<name><surname>Zhao</surname>
<given-names>X</given-names>
</name>
<name><surname>Xiao</surname>
<given-names>HL</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Transcriptomic analysis of a tertiary relict plant, extreme xerophyte <italic>Reaumuria soongorica</italic>
to identify genes related to drought adaptation</article-title>
<source>PLoS One</source>
<year>2013</year>
<volume>8</volume>
<fpage>e63993</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pone.0063993</pub-id>
<pub-id pub-id-type="pmid">23717523</pub-id>
</element-citation>
</ref>
<ref id="CR18"><label>18.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname>
<given-names>J</given-names>
</name>
<name><surname>Tian</surname>
<given-names>Q</given-names>
</name>
<name><surname>Pang</surname>
<given-names>T</given-names>
</name>
<name><surname>Jiang</surname>
<given-names>LB</given-names>
</name>
<name><surname>Wu</surname>
<given-names>RL</given-names>
</name>
<name><surname>Xia</surname>
<given-names>XL</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Deep-sequencing transcriptome analysis of low temperature perception in a desert tree, <italic>Populus euphratica</italic>
</article-title>
<source>BMC Genomics</source>
<year>2014</year>
<volume>15</volume>
<fpage>326</fpage>
<pub-id pub-id-type="doi">10.1186/1471-2164-15-326</pub-id>
<pub-id pub-id-type="pmid">24884892</pub-id>
</element-citation>
</ref>
<ref id="CR19"><label>19.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kakumanu</surname>
<given-names>A</given-names>
</name>
<name><surname>Ambavaram</surname>
<given-names>MMR</given-names>
</name>
<name><surname>Klumas</surname>
<given-names>C</given-names>
</name>
<name><surname>Krishnan</surname>
<given-names>A</given-names>
</name>
<name><surname>Batlang</surname>
<given-names>U</given-names>
</name>
<name><surname>Myers</surname>
<given-names>E</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Effects of drought on gene expression in maize reproductive and leaf meristem tissue revealed by RNA-Seq</article-title>
<source>Plant Physiol</source>
<year>2012</year>
<volume>160</volume>
<fpage>846</fpage>
<lpage>67</lpage>
<pub-id pub-id-type="doi">10.1104/pp.112.200444</pub-id>
<pub-id pub-id-type="pmid">22837360</pub-id>
</element-citation>
</ref>
<ref id="CR20"><label>20.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wu</surname>
<given-names>GA</given-names>
</name>
<name><surname>Prochnik</surname>
<given-names>S</given-names>
</name>
<name><surname>Jenkins</surname>
<given-names>J</given-names>
</name>
<name><surname>Salse</surname>
<given-names>J</given-names>
</name>
<name><surname>Hellsten</surname>
<given-names>U</given-names>
</name>
<name><surname>Murat</surname>
<given-names>F</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Sequencing of diverse mandarin, pummelo and orange genomes reveals complex history of admixture during citrus domestication</article-title>
<source>Nat Biotechnol</source>
<year>2014</year>
<volume>32</volume>
<fpage>656</fpage>
<lpage>62</lpage>
<pub-id pub-id-type="doi">10.1038/nbt.2906</pub-id>
<pub-id pub-id-type="pmid">24908277</pub-id>
</element-citation>
</ref>
<ref id="CR21"><label>21.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname>
<given-names>CK</given-names>
</name>
<name><surname>Lang</surname>
<given-names>P</given-names>
</name>
<name><surname>Dane</surname>
<given-names>F</given-names>
</name>
<name><surname>Ebel</surname>
<given-names>RC</given-names>
</name>
<name><surname>Singh</surname>
<given-names>NK</given-names>
</name>
<name><surname>Locy</surname>
<given-names>RD</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Cold acclimation induced genes of trifoliate orange (<italic>Poncirus trifoliata</italic>
)</article-title>
<source>Plant Cell Rep</source>
<year>2005</year>
<volume>23</volume>
<fpage>764</fpage>
<lpage>9</lpage>
<pub-id pub-id-type="doi">10.1007/s00299-004-0883-y</pub-id>
<pub-id pub-id-type="pmid">15449021</pub-id>
</element-citation>
</ref>
<ref id="CR22"><label>22.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Peng</surname>
<given-names>T</given-names>
</name>
<name><surname>Zhu</surname>
<given-names>XF</given-names>
</name>
<name><surname>Fan</surname>
<given-names>QJ</given-names>
</name>
<name><surname>Sun</surname>
<given-names>PP</given-names>
</name>
<name><surname>Liu</surname>
<given-names>JH</given-names>
</name>
</person-group>
<article-title>Identification and characterization of low temperature stress responsive genes in <italic>Poncirus trifoliata</italic>
by suppression subtractive hybridization</article-title>
<source>Gene</source>
<year>2012</year>
<volume>492</volume>
<fpage>220</fpage>
<lpage>8</lpage>
<pub-id pub-id-type="doi">10.1016/j.gene.2011.10.025</pub-id>
<pub-id pub-id-type="pmid">22056698</pub-id>
</element-citation>
</ref>
<ref id="CR23"><label>23.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sahin-Çevik</surname>
<given-names>M</given-names>
</name>
<name><surname>Moore</surname>
<given-names>GA</given-names>
</name>
</person-group>
<article-title>Identification and expression analysis of cold-regulated genes from the cold-hardy Citrus relative <italic>Poncirus trifoliata</italic>
(L.) Raf</article-title>
<source>Plant Mol Biol</source>
<year>2006</year>
<volume>62</volume>
<fpage>83</fpage>
<lpage>97</lpage>
<pub-id pub-id-type="doi">10.1007/s11103-006-9005-2</pub-id>
<pub-id pub-id-type="pmid">16900323</pub-id>
</element-citation>
</ref>
<ref id="CR24"><label>24.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname>
<given-names>J</given-names>
</name>
<name><surname>Sun</surname>
<given-names>PP</given-names>
</name>
<name><surname>Chen</surname>
<given-names>CL</given-names>
</name>
<name><surname>Wang</surname>
<given-names>Y</given-names>
</name>
<name><surname>Fu</surname>
<given-names>XZ</given-names>
</name>
<name><surname>Liu</surname>
<given-names>JH</given-names>
</name>
</person-group>
<article-title>An arginine decarboxylase gene <italic>PtADC</italic>
from <italic>Poncirus trifoliata</italic>
confers abiotic stress tolerance and promotes primary root growth in <italic>Arabidopsis</italic>
</article-title>
<source>J Exp Bot</source>
<year>2011</year>
<volume>62</volume>
<fpage>2899</fpage>
<lpage>914</lpage>
<pub-id pub-id-type="doi">10.1093/jxb/erq463</pub-id>
<pub-id pub-id-type="pmid">21282323</pub-id>
</element-citation>
</ref>
<ref id="CR25"><label>25.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Huang</surname>
<given-names>XS</given-names>
</name>
<name><surname>Wang</surname>
<given-names>W</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>Q</given-names>
</name>
<name><surname>Liu</surname>
<given-names>JH</given-names>
</name>
</person-group>
<article-title>A basic helix-loop-helix transcription factor, <italic>PtrbHLH</italic>
, of <italic>Poncirus trifoliata</italic>
confers cold tolerance and modulates peroxidase-mediated scavenging of hydrogen peroxide</article-title>
<source>Plant Physiol</source>
<year>2013</year>
<volume>162</volume>
<fpage>1178</fpage>
<lpage>94</lpage>
<pub-id pub-id-type="doi">10.1104/pp.112.210740</pub-id>
<pub-id pub-id-type="pmid">23624854</pub-id>
</element-citation>
</ref>
<ref id="CR26"><label>26.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Grabherr</surname>
<given-names>MG</given-names>
</name>
<name><surname>Haas</surname>
<given-names>BJ</given-names>
</name>
<name><surname>Yassour</surname>
<given-names>M</given-names>
</name>
<name><surname>Levin</surname>
<given-names>JZ</given-names>
</name>
<name><surname>Thompson</surname>
<given-names>DA</given-names>
</name>
<name><surname>Amit</surname>
<given-names>I</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Full-length transcriptome assembly from RNA-Seq data without a reference genome</article-title>
<source>Nat Biotechnol</source>
<year>2011</year>
<volume>29</volume>
<fpage>644</fpage>
<lpage>52</lpage>
<pub-id pub-id-type="doi">10.1038/nbt.1883</pub-id>
<pub-id pub-id-type="pmid">21572440</pub-id>
</element-citation>
</ref>
<ref id="CR27"><label>27.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ye</surname>
<given-names>J</given-names>
</name>
<name><surname>Fang</surname>
<given-names>L</given-names>
</name>
<name><surname>Zheng</surname>
<given-names>HK</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>Y</given-names>
</name>
<name><surname>Chen</surname>
<given-names>J</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>ZJ</given-names>
</name>
<etal></etal>
</person-group>
<article-title>WEGO: a web tool for plotting GO annotations</article-title>
<source>Nucleic Acids Res</source>
<year>2006</year>
<volume>34</volume>
<fpage>293</fpage>
<lpage>7</lpage>
<pub-id pub-id-type="doi">10.1093/nar/gkl031</pub-id>
</element-citation>
</ref>
<ref id="CR28"><label>28.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Conesa</surname>
<given-names>A</given-names>
</name>
<name><surname>Götz</surname>
<given-names>S</given-names>
</name>
<name><surname>García-Gómez</surname>
<given-names>JM</given-names>
</name>
<name><surname>Terol</surname>
<given-names>J</given-names>
</name>
<name><surname>Talón</surname>
<given-names>M</given-names>
</name>
<name><surname>Robles</surname>
<given-names>M</given-names>
</name>
</person-group>
<article-title>Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research</article-title>
<source>Bioinformatics</source>
<year>2005</year>
<volume>21</volume>
<fpage>3674</fpage>
<lpage>6</lpage>
<pub-id pub-id-type="doi">10.1093/bioinformatics/bti610</pub-id>
<pub-id pub-id-type="pmid">16081474</pub-id>
</element-citation>
</ref>
<ref id="CR29"><label>29.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Iseli</surname>
<given-names>C</given-names>
</name>
<name><surname>Jongeneel</surname>
<given-names>CV</given-names>
</name>
<name><surname>Bucher</surname>
<given-names>P</given-names>
</name>
</person-group>
<article-title>ESTScan: a program for detecting, evaluating, and reconstructing potential coding regions in EST sequences</article-title>
<source>Proc Int Conf Intell Syst Mol Biol</source>
<year>1999</year>
<volume>7</volume>
<fpage>138</fpage>
<lpage>48</lpage>
<pub-id pub-id-type="pmid">10786296</pub-id>
</element-citation>
</ref>
<ref id="CR30"><label>30.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname>
<given-names>RQ</given-names>
</name>
<name><surname>Yu</surname>
<given-names>C</given-names>
</name>
<name><surname>Li</surname>
<given-names>YR</given-names>
</name>
<name><surname>Lam</surname>
<given-names>TW</given-names>
</name>
<name><surname>Yiu</surname>
<given-names>SM</given-names>
</name>
<name><surname>Kristianse</surname>
<given-names>K</given-names>
</name>
<etal></etal>
</person-group>
<article-title>SOAP2: an improved ultrafast tool for short read alignment</article-title>
<source>Bioinformatics</source>
<year>2009</year>
<volume>25</volume>
<fpage>1966</fpage>
<lpage>7</lpage>
<pub-id pub-id-type="doi">10.1093/bioinformatics/btp336</pub-id>
<pub-id pub-id-type="pmid">19497933</pub-id>
</element-citation>
</ref>
<ref id="CR31"><label>31.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kim</surname>
<given-names>YO</given-names>
</name>
<name><surname>Kim</surname>
<given-names>JS</given-names>
</name>
<name><surname>Kang</surname>
<given-names>H</given-names>
</name>
</person-group>
<article-title>Cold-inducible zinc finger-containing glycine-rich RNA-binding protein contributes to the enhancement of freezing tolerance in <italic>Arabidopsis thaliana</italic>
</article-title>
<source>Plant J</source>
<year>2005</year>
<volume>42</volume>
<fpage>890</fpage>
<lpage>900</lpage>
<pub-id pub-id-type="doi">10.1111/j.1365-313X.2005.02420.x</pub-id>
<pub-id pub-id-type="pmid">15941401</pub-id>
</element-citation>
</ref>
<ref id="CR32"><label>32.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kim</surname>
<given-names>JY</given-names>
</name>
<name><surname>Park</surname>
<given-names>SJ</given-names>
</name>
<name><surname>Kwak</surname>
<given-names>KJ</given-names>
</name>
<name><surname>Jung</surname>
<given-names>CH</given-names>
</name>
<name><surname>Ahn</surname>
<given-names>SJ</given-names>
</name>
<name><surname>Goh</surname>
<given-names>CH</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Cold shock domain proteins and glycine-rich RNA-binding proteins from <italic>Arabidopsis thaliana</italic>
can promote the cold adaptation process in <italic>Escherichia coli</italic>
</article-title>
<source>Nucleic Acids Res</source>
<year>2007</year>
<volume>35</volume>
<fpage>506</fpage>
<lpage>16</lpage>
<pub-id pub-id-type="doi">10.1093/nar/gkl1076</pub-id>
<pub-id pub-id-type="pmid">17169986</pub-id>
</element-citation>
</ref>
<ref id="CR33"><label>33.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kovtun</surname>
<given-names>Y</given-names>
</name>
<name><surname>Chiu</surname>
<given-names>WL</given-names>
</name>
<name><surname>Tena</surname>
<given-names>G</given-names>
</name>
<name><surname>Sheen</surname>
<given-names>J</given-names>
</name>
</person-group>
<article-title>Functional analysis of oxidative stress-activated mitogen-activated protein kinase cascade in plants</article-title>
<source>Proc Natl Acad Sci U S A</source>
<year>2000</year>
<volume>97</volume>
<fpage>2940</fpage>
<lpage>5</lpage>
<pub-id pub-id-type="doi">10.1073/pnas.97.6.2940</pub-id>
<pub-id pub-id-type="pmid">10717008</pub-id>
</element-citation>
</ref>
<ref id="CR34"><label>34.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jaspers</surname>
<given-names>P</given-names>
</name>
<name><surname>Kangasjärvi</surname>
<given-names>J</given-names>
</name>
</person-group>
<article-title>Reactive oxygen species in abiotic stress signaling</article-title>
<source>Physiol Plant</source>
<year>2010</year>
<volume>138</volume>
<fpage>405</fpage>
<lpage>13</lpage>
<pub-id pub-id-type="doi">10.1111/j.1399-3054.2009.01321.x</pub-id>
<pub-id pub-id-type="pmid">20028478</pub-id>
</element-citation>
</ref>
<ref id="CR35"><label>35.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Noctor</surname>
<given-names>G</given-names>
</name>
<name><surname>Veljovic-Jovanovic</surname>
<given-names>S</given-names>
</name>
<name><surname>Driscoll</surname>
<given-names>S</given-names>
</name>
<name><surname>Novitskaya</surname>
<given-names>L</given-names>
</name>
<name><surname>Foyer</surname>
<given-names>CH</given-names>
</name>
</person-group>
<article-title>Drought and oxidative load in the leaves of C3 plants: a predominant role for photorespiration?</article-title>
<source>Ann Bot</source>
<year>2002</year>
<volume>89</volume>
<fpage>841</fpage>
<lpage>50</lpage>
<pub-id pub-id-type="doi">10.1093/aob/mcf096</pub-id>
<pub-id pub-id-type="pmid">12102510</pub-id>
</element-citation>
</ref>
<ref id="CR36"><label>36.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Shinozaki</surname>
<given-names>K</given-names>
</name>
<name><surname>Yamaguchi-Shinozaki</surname>
<given-names>K</given-names>
</name>
<name><surname>Seki</surname>
<given-names>M</given-names>
</name>
</person-group>
<article-title>Regulatory network of gene expression in the drought and cold stress responses</article-title>
<source>Curr Opin Plant Biol</source>
<year>2003</year>
<volume>6</volume>
<fpage>410</fpage>
<lpage>7</lpage>
<pub-id pub-id-type="doi">10.1016/S1369-5266(03)00092-X</pub-id>
<pub-id pub-id-type="pmid">12972040</pub-id>
</element-citation>
</ref>
<ref id="CR37"><label>37.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Luan</surname>
<given-names>S</given-names>
</name>
<name><surname>Kudla</surname>
<given-names>J</given-names>
</name>
<name><surname>Rodriguez-Concepcion</surname>
<given-names>M</given-names>
</name>
<name><surname>Yalovsky</surname>
<given-names>S</given-names>
</name>
<name><surname>Gruissem</surname>
<given-names>W</given-names>
</name>
</person-group>
<article-title>Calmodulins and calcineurin B–like proteins calcium sensors for specific signal response coupling in plants</article-title>
<source>Plant Cell</source>
<year>2002</year>
<volume>14</volume>
<fpage>389</fpage>
<lpage>400</lpage>
</element-citation>
</ref>
<ref id="CR38"><label>38.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Valmonte</surname>
<given-names>GR</given-names>
</name>
<name><surname>Arthur</surname>
<given-names>K</given-names>
</name>
<name><surname>Higgins</surname>
<given-names>CM</given-names>
</name>
<name><surname>MacDiarmid</surname>
<given-names>RM</given-names>
</name>
</person-group>
<article-title>Calcium-dependent protein kinases in plants: evolution, expression and function</article-title>
<source>Plant Cell Physiol</source>
<year>2014</year>
<volume>55</volume>
<fpage>551</fpage>
<lpage>69</lpage>
<pub-id pub-id-type="doi">10.1093/pcp/pct200</pub-id>
<pub-id pub-id-type="pmid">24363288</pub-id>
</element-citation>
</ref>
<ref id="CR39"><label>39.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Franza</surname>
<given-names>S</given-names>
</name>
<name><surname>Ehlerta</surname>
<given-names>B</given-names>
</name>
<name><surname>Liesea</surname>
<given-names>A</given-names>
</name>
<name><surname>Kurtha</surname>
<given-names>J</given-names>
</name>
<name><surname>Cazalé</surname>
<given-names>AC</given-names>
</name>
<name><surname>Romeisa</surname>
<given-names>T</given-names>
</name>
</person-group>
<article-title>Calcium-dependent protein kinase CPK21 functions in abiotic stress response in <italic>Arabidopsis thaliana</italic>
</article-title>
<source>Mol Plant</source>
<year>2011</year>
<volume>4</volume>
<fpage>83</fpage>
<lpage>96</lpage>
<pub-id pub-id-type="doi">10.1093/mp/ssq064</pub-id>
<pub-id pub-id-type="pmid">20978086</pub-id>
</element-citation>
</ref>
<ref id="CR40"><label>40.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Xu</surname>
<given-names>YJ</given-names>
</name>
<name><surname>Gao</surname>
<given-names>S</given-names>
</name>
<name><surname>Yang</surname>
<given-names>YJ</given-names>
</name>
<name><surname>Huang</surname>
<given-names>NY</given-names>
</name>
<name><surname>Cheng</surname>
<given-names>LN</given-names>
</name>
<name><surname>Wei</surname>
<given-names>Q</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Transcriptome sequencing and whole genome expression profiling of chrysanthemum under dehydration stress</article-title>
<source>BMC Genomics</source>
<year>2013</year>
<volume>14</volume>
<fpage>662</fpage>
<pub-id pub-id-type="doi">10.1186/1471-2164-14-662</pub-id>
<pub-id pub-id-type="pmid">24074255</pub-id>
</element-citation>
</ref>
<ref id="CR41"><label>41.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Janská</surname>
<given-names>A</given-names>
</name>
<name><surname>Maršík</surname>
<given-names>P</given-names>
</name>
<name><surname>Zelenková</surname>
<given-names>S</given-names>
</name>
<name><surname>Ovesná</surname>
<given-names>J</given-names>
</name>
</person-group>
<article-title>Cold stress and acclimation–what is important for metabolic adjustment?</article-title>
<source>Plant Biol</source>
<year>2010</year>
<volume>12</volume>
<fpage>395</fpage>
<lpage>405</lpage>
<pub-id pub-id-type="doi">10.1111/j.1438-8677.2009.00299.x</pub-id>
<pub-id pub-id-type="pmid">20522175</pub-id>
</element-citation>
</ref>
<ref id="CR42"><label>42.</label>
<mixed-citation publication-type="other">Lee B, Henderson DA, Zhu JK. The <italic>Arabidopsis</italic>
cold-responsive transcriptome and its regulation by ICE1. Plant Cell. 2005;17:3155–75</mixed-citation>
</ref>
<ref id="CR43"><label>43.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Peleg</surname>
<given-names>Z</given-names>
</name>
<name><surname>Blumwald</surname>
<given-names>E</given-names>
</name>
</person-group>
<article-title>Hormone balance and abiotic stress tolerance in crop plants</article-title>
<source>Curr Opin Plant Biol</source>
<year>2011</year>
<volume>14</volume>
<fpage>290</fpage>
<lpage>5</lpage>
<pub-id pub-id-type="doi">10.1016/j.pbi.2011.02.001</pub-id>
<pub-id pub-id-type="pmid">21377404</pub-id>
</element-citation>
</ref>
<ref id="CR44"><label>44.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tuteja</surname>
<given-names>N</given-names>
</name>
</person-group>
<article-title>Abscisic acid and abiotic stress signaling</article-title>
<source>Plant Signal Behav</source>
<year>2007</year>
<volume>2</volume>
<fpage>135</fpage>
<lpage>8</lpage>
<pub-id pub-id-type="doi">10.4161/psb.2.3.4156</pub-id>
<pub-id pub-id-type="pmid">19516981</pub-id>
</element-citation>
</ref>
<ref id="CR45"><label>45.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Heino</surname>
<given-names>P</given-names>
</name>
<name><surname>Sandman</surname>
<given-names>G</given-names>
</name>
<name><surname>Lång</surname>
<given-names>V</given-names>
</name>
<name><surname>Nordin</surname>
<given-names>K</given-names>
</name>
<name><surname>Palva</surname>
<given-names>ET</given-names>
</name>
</person-group>
<article-title>Abscisic acid deficiency prevents development of freezing tolerance in <italic>Arabidopsis thaliana</italic>
(L.) Heynh</article-title>
<source>Theor Appl Genet</source>
<year>1990</year>
<volume>79</volume>
<fpage>801</fpage>
<lpage>6</lpage>
<pub-id pub-id-type="doi">10.1007/BF00224248</pub-id>
<pub-id pub-id-type="pmid">24226742</pub-id>
</element-citation>
</ref>
<ref id="CR46"><label>46.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lang</surname>
<given-names>V</given-names>
</name>
<name><surname>Mantyla</surname>
<given-names>E</given-names>
</name>
<name><surname>Welin</surname>
<given-names>B</given-names>
</name>
<name><surname>Palva</surname>
<given-names>ET</given-names>
</name>
</person-group>
<article-title>Alterations in water status, endogenous abscisic acid content, and expression of rab18 gene during the development of freezing tolerance in <italic>Arabidopsis thaliana</italic>
</article-title>
<source>Plant Physiol</source>
<year>1994</year>
<volume>104</volume>
<fpage>1341</fpage>
<lpage>9</lpage>
<pub-id pub-id-type="pmid">12232173</pub-id>
</element-citation>
</ref>
<ref id="CR47"><label>47.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Espasandin</surname>
<given-names>FD</given-names>
</name>
<name><surname>Maiale</surname>
<given-names>SJ</given-names>
</name>
<name><surname>Calzadilla</surname>
<given-names>P</given-names>
</name>
<name><surname>Ruiz</surname>
<given-names>OA</given-names>
</name>
<name><surname>Sansberro</surname>
<given-names>PA</given-names>
</name>
</person-group>
<article-title>Transcriptional regulation of 9-<italic>cis</italic>
-epoxycarotenoid dioxygenase (NCED) gene by putrescine accumulation positively modulates ABA synthesis and drought tolerance in <italic>Lotus tenuis</italic>
plants</article-title>
<source>Plant Physiol Biochem</source>
<year>2014</year>
<volume>76</volume>
<fpage>29</fpage>
<lpage>35</lpage>
<pub-id pub-id-type="doi">10.1016/j.plaphy.2013.12.018</pub-id>
<pub-id pub-id-type="pmid">24448322</pub-id>
</element-citation>
</ref>
<ref id="CR48"><label>48.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tähtiharju</surname>
<given-names>S</given-names>
</name>
<name><surname>Palva</surname>
<given-names>T</given-names>
</name>
</person-group>
<article-title>Antisense inhibition of protein phosphatase 2C accelerates cold acclimation in <italic>Arabidopsis thaliana</italic>
</article-title>
<source>Plant J</source>
<year>2001</year>
<volume>26</volume>
<fpage>461</fpage>
<lpage>70</lpage>
<pub-id pub-id-type="doi">10.1046/j.1365-313X.2001.01048.x</pub-id>
<pub-id pub-id-type="pmid">11439132</pub-id>
</element-citation>
</ref>
<ref id="CR49"><label>49.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cao</surname>
<given-names>WH</given-names>
</name>
<name><surname>Liu</surname>
<given-names>J</given-names>
</name>
<name><surname>He</surname>
<given-names>XJ</given-names>
</name>
<name><surname>Mu</surname>
<given-names>RL</given-names>
</name>
<name><surname>Zhou</surname>
<given-names>HL</given-names>
</name>
<name><surname>Chen</surname>
<given-names>SY</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Modulation of ethylene responses affects plant salt-stress responses</article-title>
<source>Plant Physiol</source>
<year>2007</year>
<volume>143</volume>
<fpage>707</fpage>
<lpage>19</lpage>
<pub-id pub-id-type="doi">10.1104/pp.106.094292</pub-id>
<pub-id pub-id-type="pmid">17189334</pub-id>
</element-citation>
</ref>
<ref id="CR50"><label>50.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jung</surname>
<given-names>JY</given-names>
</name>
<name><surname>Shin</surname>
<given-names>R</given-names>
</name>
<name><surname>Schachtman</surname>
<given-names>DP</given-names>
</name>
</person-group>
<article-title>Ethylene mediates response and tolerance to potassium deprivation in Arabidopsis</article-title>
<source>Plant Cell</source>
<year>2009</year>
<volume>21</volume>
<fpage>607</fpage>
<lpage>21</lpage>
<pub-id pub-id-type="doi">10.1105/tpc.108.063099</pub-id>
<pub-id pub-id-type="pmid">19190240</pub-id>
</element-citation>
</ref>
<ref id="CR51"><label>51.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kendrick</surname>
<given-names>MD</given-names>
</name>
<name><surname>Chang</surname>
<given-names>C</given-names>
</name>
</person-group>
<article-title>Ethylene signaling: new levels of complexity and regulation</article-title>
<source>Curr Opin Plant Boil</source>
<year>2008</year>
<volume>11</volume>
<fpage>479</fpage>
<lpage>85</lpage>
<pub-id pub-id-type="doi">10.1016/j.pbi.2008.06.011</pub-id>
</element-citation>
</ref>
<ref id="CR52"><label>52.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chu</surname>
<given-names>C</given-names>
</name>
<name><surname>Lee</surname>
<given-names>TM</given-names>
</name>
</person-group>
<article-title>The relationship between ethylene biosynthesis and chilling tolerance in seedlings of rice (<italic>Oryza sativa L.</italic>
)</article-title>
<source>Bot Bull Acad Sin</source>
<year>1989</year>
<volume>30</volume>
<fpage>263</fpage>
<lpage>73</lpage>
</element-citation>
</ref>
<ref id="CR53"><label>53.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ciardi</surname>
<given-names>JA</given-names>
</name>
<name><surname>Deikman</surname>
<given-names>J</given-names>
</name>
<name><surname>Orzolek</surname>
<given-names>MD</given-names>
</name>
</person-group>
<article-title>Increased ethylene synthesis enhances chilling tolerance in tomato</article-title>
<source>Physiol Plant</source>
<year>1997</year>
<volume>101</volume>
<fpage>333</fpage>
<lpage>40</lpage>
<pub-id pub-id-type="doi">10.1111/j.1399-3054.1997.tb01005.x</pub-id>
</element-citation>
</ref>
<ref id="CR54"><label>54.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yu</surname>
<given-names>XM</given-names>
</name>
<name><surname>Griffith</surname>
<given-names>M</given-names>
</name>
<name><surname>Wiseman</surname>
<given-names>SB</given-names>
</name>
</person-group>
<article-title>Ethylene induces antifreeze activity in winter rye leaves</article-title>
<source>Plant Physiol</source>
<year>2001</year>
<volume>126</volume>
<fpage>1232</fpage>
<lpage>40</lpage>
<pub-id pub-id-type="doi">10.1104/pp.126.3.1232</pub-id>
<pub-id pub-id-type="pmid">11457973</pub-id>
</element-citation>
</ref>
<ref id="CR55"><label>55.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname>
<given-names>Z</given-names>
</name>
<name><surname>Huang</surname>
<given-names>R</given-names>
</name>
</person-group>
<article-title>Enhanced tolerance to freezing in tobacco and tomato overexpressing transcription factor <italic>TERF2/LeERF2</italic>
is modulated by ethylene biosynthesis</article-title>
<source>Plant Mol Biol</source>
<year>2010</year>
<volume>73</volume>
<fpage>241</fpage>
<lpage>9</lpage>
<pub-id pub-id-type="doi">10.1007/s11103-010-9609-4</pub-id>
<pub-id pub-id-type="pmid">20135196</pub-id>
</element-citation>
</ref>
<ref id="CR56"><label>56.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Catalá</surname>
<given-names>R</given-names>
</name>
<name><surname>López-Cobollo</surname>
<given-names>R</given-names>
</name>
<name><surname>Castellano</surname>
<given-names>MM</given-names>
</name>
<name><surname>Angostob</surname>
<given-names>T</given-names>
</name>
<name><surname>Alonsoc</surname>
<given-names>JM</given-names>
</name>
<name><surname>Eckerc</surname>
<given-names>JR</given-names>
</name>
<etal></etal>
</person-group>
<article-title>The Arabidopsis 14-3-3 protein RARE COLD INDUCIBLE 1A links low-temperature response and ethylene biosynthesis to regulate freezing tolerance and cold acclimation</article-title>
<source>Plant Cell</source>
<year>2014</year>
<volume>26</volume>
<fpage>3326</fpage>
<lpage>42</lpage>
<pub-id pub-id-type="doi">10.1105/tpc.114.127605</pub-id>
<pub-id pub-id-type="pmid">25122152</pub-id>
</element-citation>
</ref>
<ref id="CR57"><label>57.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Shi</surname>
<given-names>Y</given-names>
</name>
<name><surname>Tian</surname>
<given-names>S</given-names>
</name>
<name><surname>Hou</surname>
<given-names>L</given-names>
</name>
<name><surname>Huang</surname>
<given-names>XZ</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>XY</given-names>
</name>
<name><surname>Guo</surname>
<given-names>HW</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Ethylene signaling negatively regulates freezing tolerance by repressing expression of CBF and type-A ARR genes in Arabidopsis</article-title>
<source>Plant Cell</source>
<year>2012</year>
<volume>24</volume>
<fpage>2578</fpage>
<lpage>95</lpage>
<pub-id pub-id-type="doi">10.1105/tpc.112.098640</pub-id>
<pub-id pub-id-type="pmid">22706288</pub-id>
</element-citation>
</ref>
<ref id="CR58"><label>58.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Colebrook</surname>
<given-names>EH</given-names>
</name>
<name><surname>Thomas</surname>
<given-names>SG</given-names>
</name>
<name><surname>Phillips</surname>
<given-names>AL</given-names>
</name>
<name><surname>Hedden</surname>
<given-names>P</given-names>
</name>
</person-group>
<article-title>The role of gibberellin signalling in plant responses to abiotic stress</article-title>
<source>J Exp Biol</source>
<year>2014</year>
<volume>217</volume>
<fpage>67</fpage>
<lpage>75</lpage>
<pub-id pub-id-type="doi">10.1242/jeb.089938</pub-id>
<pub-id pub-id-type="pmid">24353205</pub-id>
</element-citation>
</ref>
<ref id="CR59"><label>59.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zentella</surname>
<given-names>R</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>ZL</given-names>
</name>
<name><surname>Park</surname>
<given-names>M</given-names>
</name>
<name><surname>Thomas</surname>
<given-names>SG</given-names>
</name>
<name><surname>Endo</surname>
<given-names>A</given-names>
</name>
<name><surname>Murase</surname>
<given-names>K</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Global analysis of DELLA direct targets in early gibberellin signaling in <italic>Arabidopsis</italic>
</article-title>
<source>Plant Cell</source>
<year>2007</year>
<volume>19</volume>
<fpage>3037</fpage>
<lpage>57</lpage>
<pub-id pub-id-type="doi">10.1105/tpc.107.054999</pub-id>
<pub-id pub-id-type="pmid">17933900</pub-id>
</element-citation>
</ref>
<ref id="CR60"><label>60.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Achard</surname>
<given-names>P</given-names>
</name>
<name><surname>Gong</surname>
<given-names>F</given-names>
</name>
<name><surname>Cheminant</surname>
<given-names>S</given-names>
</name>
<name><surname>Alioua</surname>
<given-names>M</given-names>
</name>
<name><surname>Hedden</surname>
<given-names>P</given-names>
</name>
<name><surname>Genschik</surname>
<given-names>P</given-names>
</name>
</person-group>
<article-title>The cold-inducible CBF1 factor–dependent signaling pathway modulates the accumulation of the growth-repressing DELLA proteins via its effect on gibberellin metabolism</article-title>
<source>Plant Cell</source>
<year>2008</year>
<volume>20</volume>
<fpage>2117</fpage>
<lpage>29</lpage>
<pub-id pub-id-type="doi">10.1105/tpc.108.058941</pub-id>
<pub-id pub-id-type="pmid">18757556</pub-id>
</element-citation>
</ref>
<ref id="CR61"><label>61.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Miller</surname>
<given-names>GAD</given-names>
</name>
<name><surname>Mittler</surname>
<given-names>R</given-names>
</name>
</person-group>
<article-title>Could heat shock transcription factors function as hydrogen peroxide sensors in plants?</article-title>
<source>Ann Bot</source>
<year>2006</year>
<volume>98</volume>
<fpage>279</fpage>
<lpage>88</lpage>
<pub-id pub-id-type="doi">10.1093/aob/mcl107</pub-id>
<pub-id pub-id-type="pmid">16740587</pub-id>
</element-citation>
</ref>
<ref id="CR62"><label>62.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Suzuki</surname>
<given-names>N</given-names>
</name>
<name><surname>Bajad</surname>
<given-names>S</given-names>
</name>
<name><surname>Shuman</surname>
<given-names>J</given-names>
</name>
<name><surname>Shulaev</surname>
<given-names>V</given-names>
</name>
<name><surname>Mittler</surname>
<given-names>R</given-names>
</name>
</person-group>
<article-title>The transcriptional co-activator MBF1c is a key regulator of thermotolerance in <italic>Arabidopsis thaliana</italic>
</article-title>
<source>J Biol Chem</source>
<year>2008</year>
<volume>283</volume>
<fpage>9269</fpage>
<lpage>75</lpage>
<pub-id pub-id-type="doi">10.1074/jbc.M709187200</pub-id>
<pub-id pub-id-type="pmid">18201973</pub-id>
</element-citation>
</ref>
<ref id="CR63"><label>63.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Agarwal</surname>
<given-names>PK</given-names>
</name>
<name><surname>Jha</surname>
<given-names>B</given-names>
</name>
</person-group>
<article-title>Transcription factors in plants and ABA dependent and independent abiotic stress signalling</article-title>
<source>Biol Plant</source>
<year>2010</year>
<volume>54</volume>
<fpage>201</fpage>
<lpage>12</lpage>
<pub-id pub-id-type="doi">10.1007/s10535-010-0038-7</pub-id>
</element-citation>
</ref>
<ref id="CR64"><label>64.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname>
<given-names>JH</given-names>
</name>
<name><surname>Peng</surname>
<given-names>T</given-names>
</name>
<name><surname>Dai</surname>
<given-names>W</given-names>
</name>
</person-group>
<article-title>Critical <italic>cis</italic>
-acting elements and interacting transcription factors: key players associated with abiotic stress responses in plants</article-title>
<source>Plant Mol Biol Rep</source>
<year>2014</year>
<volume>32</volume>
<fpage>303</fpage>
<lpage>17</lpage>
<pub-id pub-id-type="doi">10.1007/s11105-013-0667-z</pub-id>
</element-citation>
</ref>
<ref id="CR65"><label>65.</label>
<mixed-citation publication-type="other">Chinnusamy V, Zhu JK, Sunkar R. Gene regulation during cold stress acclimation in plants. Plant Stress Tolerance. 2010;639:39–55.</mixed-citation>
</ref>
<ref id="CR66"><label>66.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname>
<given-names>LC</given-names>
</name>
<name><surname>Liu</surname>
<given-names>GX</given-names>
</name>
<name><surname>Zhao</surname>
<given-names>GY</given-names>
</name>
<name><surname>Xia</surname>
<given-names>C</given-names>
</name>
<name><surname>Jia</surname>
<given-names>JZ</given-names>
</name>
<name><surname>Liu</surname>
<given-names>X</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Characterization of a wheat R2R3-MYB transcription factor gene, <italic>TaMYB19</italic>
, involved in enhanced abiotic stresses in <italic>Arabidopsis</italic>
</article-title>
<source>Plant Cell Physiol</source>
<year>2014</year>
<volume>55</volume>
<fpage>1802</fpage>
<lpage>12</lpage>
<pub-id pub-id-type="doi">10.1093/pcp/pcu109</pub-id>
<pub-id pub-id-type="pmid">25146486</pub-id>
</element-citation>
</ref>
<ref id="CR67"><label>67.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tiburcio</surname>
<given-names>AF</given-names>
</name>
<name><surname>Altabella</surname>
<given-names>T</given-names>
</name>
<name><surname>Bitrián</surname>
<given-names>M</given-names>
</name>
<name><surname>Alcázar</surname>
<given-names>R</given-names>
</name>
</person-group>
<article-title>The roles of polyamines during the lifespan of plants: from development to stress</article-title>
<source>Planta</source>
<year>2014</year>
<volume>240</volume>
<fpage>1</fpage>
<lpage>18</lpage>
<pub-id pub-id-type="doi">10.1007/s00425-014-2055-9</pub-id>
<pub-id pub-id-type="pmid">24659098</pub-id>
</element-citation>
</ref>
<ref id="CR68"><label>68.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pertea</surname>
<given-names>G</given-names>
</name>
<name><surname>Huang</surname>
<given-names>X</given-names>
</name>
<name><surname>Liang</surname>
<given-names>F</given-names>
</name>
<name><surname>Antonescu</surname>
<given-names>V</given-names>
</name>
<name><surname>Sultana</surname>
<given-names>R</given-names>
</name>
<name><surname>Karamycheva</surname>
<given-names>S</given-names>
</name>
<etal></etal>
</person-group>
<article-title>TIGR Gene Indices clustering tools (TGICL): a software system for fast clustering of large EST datasets</article-title>
<source>Bioinformatics</source>
<year>2003</year>
<volume>19</volume>
<fpage>651</fpage>
<lpage>2</lpage>
<pub-id pub-id-type="doi">10.1093/bioinformatics/btg034</pub-id>
<pub-id pub-id-type="pmid">12651724</pub-id>
</element-citation>
</ref>
<ref id="CR69"><label>69.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kanehisa</surname>
<given-names>M</given-names>
</name>
<name><surname>Goto</surname>
<given-names>S</given-names>
</name>
</person-group>
<article-title>KEGG: kyoto encyclopedia of genes and genomes</article-title>
<source>Nucleic Acids Res</source>
<year>2000</year>
<volume>28</volume>
<fpage>27</fpage>
<lpage>30</lpage>
<pub-id pub-id-type="doi">10.1093/nar/28.1.27</pub-id>
<pub-id pub-id-type="pmid">10592173</pub-id>
</element-citation>
</ref>
<ref id="CR70"><label>70.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mortazavi</surname>
<given-names>A</given-names>
</name>
<name><surname>Williams</surname>
<given-names>BA</given-names>
</name>
<name><surname>McCue</surname>
<given-names>K</given-names>
</name>
<name><surname>Schaeffer</surname>
<given-names>L</given-names>
</name>
<name><surname>Wold</surname>
<given-names>B</given-names>
</name>
</person-group>
<article-title>Mapping and quantifying mammalian transcriptomes by RNA-Seq</article-title>
<source>Nat Methods</source>
<year>2008</year>
<volume>5</volume>
<fpage>621</fpage>
<lpage>8</lpage>
<pub-id pub-id-type="doi">10.1038/nmeth.1226</pub-id>
<pub-id pub-id-type="pmid">18516045</pub-id>
</element-citation>
</ref>
<ref id="CR71"><label>71.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Fu</surname>
<given-names>XZ</given-names>
</name>
<name><surname>Chen</surname>
<given-names>CW</given-names>
</name>
<name><surname>Wang</surname>
<given-names>Y</given-names>
</name>
<name><surname>Liu</surname>
<given-names>JH</given-names>
</name>
<name><surname>Moriguchi</surname>
<given-names>T</given-names>
</name>
</person-group>
<article-title>Ectopic expression of <italic>MdSPDS1</italic>
in sweet orange (<italic>Citrus sinensis</italic>
Osbeck) reduces canker susceptibility: involvement of H<sub>2</sub>
O<sub>2</sub>
production and transcriptional alteration</article-title>
<source>BMC Plant Biol</source>
<year>2011</year>
<volume>11</volume>
<fpage>55</fpage>
<pub-id pub-id-type="doi">10.1186/1471-2229-11-55</pub-id>
<pub-id pub-id-type="pmid">21439092</pub-id>
</element-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Bois/explor/OrangerV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 0005888 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 0005888 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Wicri/Bois |area= OrangerV1 |flux= Pmc |étape= Corpus |type= RBID |clé= |texte= }}
![]() | This area was generated with Dilib version V0.6.25. | ![]() |