雌雄蕊原基分化期低温胁迫下2个小麦品种幼穗miRNA表达谱分析

doi:10.7668/hbnxb.20191253

摘要/Abstract

摘要： 为探讨春季低温（倒春寒）胁迫下抗倒春寒能力不同的小麦品种相关miRNA及其响应规律，提供miRNA的表达谱和差异表达miRNA的信息，揭示差异表达miRNA在小麦抗倒春寒中的作用，采用高通量测序（Small RNA-seq）技术对雌雄蕊原基分化期0℃低温胁迫72 h及未低温处理（对照）的矮抗58（AK58）和郑麦366（ZM366）幼穗进行测序，通过生物信息学分析，以padj<0.05且|log₂（Fold_change）|≥1为条件，筛选低温胁迫下差异表达显著的miRNA，利用psRobot（v1.2）对差异表达miRNA靶基因进行预测，然后对差异表达显著miRNA的靶基因进行GO富集分析和KEGG Pathway分析。结果表明，小麦幼穗测序共鉴定到112类miRNA，分属于38个不同的MIRNA家族。AK58、ZM366分别鉴定出85，88个差异表达miRNA，AK58有27个miRNA表达有显著差异，其中23个显著上调表达，4个显著下调表达，差异表达极显著的miRNA为2个，分别是tae-miR9672b、tae-miR9672a-3p；ZM366有48个差异表达显著的miRNA，13个显著上调表达，35个显著下调表达，差异表达极显著的miRNA为4个，分别为tae-miR9672b、tae-miR9672a-3p、tae-miR6201、tae-miR9674b-5p。对低温胁迫下AK58与ZM366差异表达的miRNA进行靶基因预测，AK58 85个差异表达miRNA对应848个靶基因；ZM366 88个差异表达miRNA对应6 537个靶基因。GO功能富集分析结果显示，AK58差异表达miRNA靶基因有20个生物学过程、6个细胞成分、6个分子功能类别极显著富集（FDR<0.01）；ZM366有10个生物学过程、14个细胞成分、19个分子功能类别极显著富集（FDR<0.01）。对2个小麦品种幼穗中差异表达miRNA靶基因进行KEGG代谢通路（Pathway）富集分析，AK58差异表达miRNA靶基因显著富集（P<0.05）到RNA聚合酶，嘌呤代谢，嘧啶代谢，光合生物固碳途径，自噬，托烷、哌啶和吡啶生物碱的生物合成，核苷酸切除修复，错配修复，DNA复制，核糖体，烟酸和烟酰胺代谢11个代谢通路；ZM366差异表达miRNA靶基因显著富集（P<0.05）在植物-病原互作，植物昼夜节律，缬氨酸、亮氨酸和异亮氨酸降解，细胞周期，苯丙氨酸代谢，苯丙烷生物合成，RNA聚合酶，植物激素信号转导，氰胺酸代谢，次生代谢产物的生物合成，光合生物中碳固定，丙酮酸代谢，类胡萝卜素生物合成13个代谢通路。miRNA调控mRNA参与小麦应答低温胁迫中，小麦可能通过miRNA171a、miRNA156、miRNA398等多个miRNA的介导调控来抵御倒春寒胁迫，植物昼夜节律代谢途径可能与低温胁迫密切相关。2个小麦品种抗倒春寒能力不同的原因可能与miRNA调控光形态建成转导途径和开花过程的靶基因表达模式不同有关。

关键词: 小麦, 幼穗, 低温胁迫, RNA-seq, miRNA

Abstract: To explore the relevant miRNAs and their response laws of wheat varieties with different resistance to cold in spring under the stress of low temperature in spring(late spring cold), high-throughput sequencing(Small RNA-seq) technology was used to analyze the young panicles of Aikang 58(AK58) and Zhengmai 366(ZM366) in the differentiation stage of pistil and stamen primordium at 0℃ low temperature stress for 72 h and without low temperature treatment(control).The purpose of sequencing is provide miRNA expression profiles and information on differentially expressed miRNAs, and reveal the role of differentially expressed miRNAs in wheat resistance to late spring cold. Through bioinformatics analysis, under the condition of padj<0.05 and|log₂(Fold_change)| ≥ 1, the miRNAs with significant differential expression was screened under low temperature stress;using psRobot(v1.2) to predict the differentially expressed miRNA target genes, and then target genes with significantly differentially expressed miRNAs were subjected to GO enrichment analysis and KEGG Pathway analysis. The results showed that 112 types of miRNAs were identified by sequencing of wheat ears, belonging to 38 different MIRNA families.AK58 and ZM366 identified 85 and 88 differentially expressed miRNAs, respectively. AK58 had 27 miRNAs with significant differences in expression, 23 of which were significantly up-regulated, 4 were significantly down-regulated, and 2 miRNAs with extremely significant differential expression were tae-miR9672b, tae-miR9672a-3p;ZM366 had 48 miRNAs with significant differential expression, 13 with significantly up-regulated expression, 35 with significantly down-regulated expression, 4 miRNAs with extremely significant differential expression, tae-miR9672b, tae-miR9672a-3p, tae-miR6201, tae-miR9674b-5p.Target genes were predicted for the differentially expressed miRNAs of AK58 and ZM366 under low temperature stress. 85 differentially expressed miRNAs of AK58 correspond to 848 target genes;88 differentially expressed miRNAs of ZM366 correspond to 6 537 target genes.GO function enrichment analysis results showed that AK58 differentially expressed miRNA target genes had 20 biological processes, 6 cell components, and 6 molecular function categories extremely significantly enriched(FDR<0.01).ZM366 had 10 biological processes, 14 cell components, and 19 molecular functional categories extremely significantly enriched(FDR<0.01).The KEGG metabolic pathway(Pathway) enrichment analysis of the differentially expressed miRNA target genes in the young ears of two wheat varieties showed that AK58 differentially expressed miRNA target genes were significantly enriched(P<0.05) to RNA polymerase, purine metabolism, and pyrimidine metabolism, photosynthetic organism carbon fixation pathway, autophagy, tropane, piperidine and pyridine alkaloid biosynthesis, nucleotide excision repair, mismatch repair, DNA replication, ribosome, niacin and nicotinamide metabolism 11 metabolic pathways.ZM366 differentially expressed miRNA target genes were significantly enriched(P<0.05) in plant pathogen interactions, plant circadian rhythms, degradation of valine, leucine and isoleucine, cell cycle, phenylalanine metabolism, propionic acid biology synthesis, RNA polymerase, plant hormone signal transduction, cyanamide acid metabolism, biosynthesis of secondary metabolites, carbon fixation of photosynthetic organisms, pyruvate metabolism, carotenoid biosynthesis, 13 metabolic pathways.Wheat miRNA regulation of mRNA participates in wheat response to low temperature stress. It is possible that miRNAs such as miRNA171a, miRNA156, miRNA398 and other miRNAs can mediate regulation to resist cold stress. Plant circadian rhythms may be closely related to low temperature stress.The reason for the different resistance to spring cold of the two wheat varieties may be related to the different expression patterns of target genes in the miRNA-regulated photomorphogenesis pathway and flowering process.

Key words: Wheat, Young spike, Low-temperature stress, RNA-seq, miRNA

中图分类号:

S512.03
Q78

张自阳, 王斌, 王智煜, 王志伟, 朱启迪, 茹振钢, 刘明久. 雌雄蕊原基分化期低温胁迫下2个小麦品种幼穗miRNA表达谱分析[J]. 华北农学报, 2021, 36(1): 81-94. doi: 10.7668/hbnxb.20191253.

ZHANG Ziyang, WANG Bin, WANG Zhiyu, WANG Zhiwei, ZHU Qidi, RU Zhengang, LIU Mingjiu. Analysis of miRNA Expression Profile in Young Ear of Two Wheat Varieties Under Low Temperature Stress at the Female Stamen Primordium Differentiation Stage[J]. ACTA AGRICULTURAE BOREALI-SINICA, 2021, 36(1): 81-94. doi: 10.7668/hbnxb.20191253.

http://www.hbnxb.net/CN/Y2021/V36/I1/81

参考文献

[1] 赵虹,王西成,胡卫国,曹廷杰,李博.黄淮南片麦区小麦倒春寒冻害成因及预防措施[J].河南农业科学,2014,43(8):34-38.doi:10.3969/j.issn.1004-3268.2014.08.008.
Zhao H,Wang X C,Hu W G,Cao T J,Li B.Genetic analysis and countermeasures of wheat late-spring-coldness injury in South Huang-huai wheat region[J]. Journal of Henan Agricultural Sciences,2014,43(8):34-38.
[2] Sanghera G S,Wani S H,Hussain W,Singh N B.Engineering cold stress tolerance in crop plants[J]. Current Genomics,2011,12(1):30-43.doi:10.2174/138920211794520178.
[3] 张自阳,王智煜,王斌,王志伟,朱启迪,霍云凤,茹振钢,刘明久.春季穗分化阶段低温处理对不同小麦品种幼穗结实性及生理特性的影响[J].华北农学报,2019,34(4):130-139.doi:10.7668/hbnxb.20190321.
Zhang Z Y,Wang Z Y,Wang B,Wang Z W,Zhu Q D,Huo Y F,Ru Z G,Liu M J.Effects of low temperature treatment at spring spike differentiation stage on young ear fruiting and physiological characteristics of different wheat varieties[J]. Acta Agriculturae Boreali-Sinica,2019,34(4):130-139.
[4] Fowler D B,Limin A E.Interactions among factors regulating phenological development and acclimation rate determine low-temperature tolerance in wheat[J]. Annals of Botany,2004,94(5):717-724.doi:10.1093/aob/mch196.
[5] Paulsen G M,Heyne E G.Grain production of winter wheat after spring freeze injury[J]. Agronomy Journal,1983,75(4):705-707.doi:10.2134/agronj 1983.00021962007500040031x.
[6] Maqbool A,Shafiq S,Lake L.Radiant frost tolerance in pulse crops-a review[J]. Euphytica,2010,172(1):1-12.doi:10.1007/s10681-009-0031-4.
[7] 高艳,唐建卫,殷贵鸿,韩玉林,黄峰,王丽娜,于海飞,李楠楠,张倩,杨光宇,李新平.倒春寒发生时期和次数对冬小麦产量性状的影响[J].麦类作物学报,2015,35(5):687-692.doi:10.7606/j.issn.1009-1041.2015.05.17 Gao Y,Tang J W,Yin G H,Han Y L,Hung F,Wang L N,Yu H F,Li N N,Zhang Q,Yang G Y,Li X P.Effect of different periods and frequency of late spring coldness on winter wheat yield related traits[J]. Journal of Triticeae Crops,2015,35(5):687-692.
[8] 刘璇,王瑞丽,周伟,方保停,郑宏远,张艳林,詹克慧.春季低温对冬小麦穗部发育和粒重的影响[J].河南农业大学学报,2013,47(4):373-380.doi:10.3969/j.issn.1000-2340.2013.04.001 Liu X,Wang R L,Zhou W,Fang B T,Zheng H Y,Zhang Y L,Zhan K H.Effect of spring low temperature on ear development and grain weight of winter wheat[J]. Journal of Henan Agricultural University,2013,47(4):373-380.
[9] Rinalducci S,Egidi M G,Karimzadeh G,Jazii F R,Zolla L.Proteomic analysis of a spring wheat cultivar in response to prolonged cold stress[J]. Electrophoresis,2011,32(14):1807-1818.doi:10.1002/elps.201000663.
[10] 吴青霞,杨林,邵慧,冉从福,杨子博,余静,李立群,李学军.药隔期低温胁迫对小麦生理及产量的影响[J].麦类作物学报,2013,33(4):752-757.doi:10.7606/j.issn.1009-1041.2013.04.023 Wu Q X,Yang L,Shao H,Ran C F,Yang Z B,Yu J,Li L Q,Li X J.Effects of low temperature stress at anther connective formation phase on physiological characteristics and yield of wheat[J]. Journal of Triticeae Crops,2013,33(4):752-757.
[11] Vítámvás P,Saalbach G,Prášil I T, apková V,Opatrná J,Ahmed J.WCS120 protein family and proteins soluble upon boiling in cold-acclimated winter wheat[J]. Journal of Plant Physiology,2007,164(9):1197-1207.doi:10.1016/j.jplph.2006.06.011.
[12] Zhang S J,Song G Q,Li Y L,Gao J,Wang J,Chen G J,Li H S,Li G Y, Zhao Z D.Comparative proteomic analysis of cold responsive proteins in two wheat cultivars with different tolerance to spring radiation frost[J]. Frontiers of Agricultural Science and Engineering,2014,1(1):37-45.doi:10.15302/J-FASE-2014008.
[13] Monroy A F,Dryanova A,Malette B,Oren D H,Farajalla M R,Liu W C,Danyluk J,Ubayasena L W C,Kane K,Scoles G J,Sarhan F, Gulick P J.Regulatory gene candidates and gene expression analysis of cold acclimation in winter and spring wheat[J]. Plant Molecular Biology,2007,64(4):409-423.doi:10.1007/s11103-007-9161-z.
[14] Winfield M O,Lu C G,Wilson I D,Coghill J A,Edwards K J.Plant responses to cold:Transcriptome analysis of wheat[J]. Plant Biotechnology Journal 2010,8(7):749-771.doi:10.1111/j.1467-7652.2010.00536.x.
[15] Herman E M,Rotter K,Premakumar R,Elwinger G,Bae R,Ehler-King L,Chen S X,Livingston D P.Additional freeze hardiness in wheat acquired by exposure to-3℃ is associated with extensive physiological,morphological,and molecular changes[J]. Journal of Experimental Botany,2006,57(14):3601-3618.doi:10.1093/jxb/erl111.
[16] Kang G Z,Li G Z,Yang W P,Han Q X,Ma H Z,Wang Y H,Ren J P, Zhu Y J, Guo T C.Transcriptional profile of the spring freeze response in the leaves of bread wheat (Triticum aestivum L.)[J]. Acta Physiologiae Plantarum,2013,35(2):575-587.doi:10.1007/s11738-012-1099-3.
[17] Fuller M P,Fuller A M,Kaniouras S,Christophers J,Fredericks T.The freezing characteristics of wheat at ear emergence[J]. European Journal of Agronomy,2007,26(4):435-441.doi:10.1016/j.eja.2007.01.001.
[18] 王瑞霞,闫长生,张秀英,孙果忠,钱兆国,亓晓蕾,牟秋焕,肖世和.春季低温对小麦产量和光合特性的影响[J].作物学报,2018,44(2):288-296.doi:10.3724/SP.J.1006.2018.00288.
Wang R X,Yan C S,Zhang X Y,Sun G Z,Qian Z G,Yuan X L,Mou Q H,Xiao S H.Effect of low temperature in spring on yield and photosynthetic characteristics of wheat[J]. Acta Agronomica Sinica,2018,44(2):288-296.
[19] Voinnet O.Origin,biogenesis,and activity of plant MicroRNAs[J]. Cell,2009,136(4):669-687.doi:10.1016/j.cell.2009.01.046.
[20] Budak H,Kantar M.Harnessing NGS and big data optimally:comparison of miRNA prediction from assembled versus non-assembled sequencing data:the case of the grass Aegilops tauschii complex genome[J]. Omics:A Journal of Integrative Biology,2015,19(7):407-415.doi:10.1089/omi.2015.0038.
[21] Budak H,Kantar M,Bulut R,Akpinar B A.Stress responsive miRNAs and isomiRs in cereals[J]. Plant Science,2015,235:1-13.doi:10.1016/j.plantsci.2015.02.008.
[22] Sunkar R,Zhu J K.Novel and stress-regulated microRNAs and other small RNAs from Arabidopsis[J]. Plant Cell,2004,16(8):2001-2019.doi:10.1105/tpc.104.022830.
[23] Zhang J Y,Xu Y Y,Huan Q,Chong K.Deep sequencing of Brachypodium small RNAs at the global genome level identifies microRNAs involved in cold stress response[J]. BMC Genomics,2009,10(1):449.doi:10.1186/1471-2164-10-449.
[24] Tang Z H,Zhang L P,Xu C G,Yuan S H,Zhang F T,Zheng Y L,Zhao C P.Uncovering small RNA-mediated responses to cold stress in a wheat thermosensitive genic male-sterile line by deep sequencing[J]. Plant Physiology,2012,159(2):721-738.doi:10.1104/pp.112.196048.
[25] Shu Y J,Liu Y,Li W,Song L L,Zhang J,Guo C H.Genome-wide investigation of microRNAs and their targets in response to freezing stress in Medicago sativa L.,based on high-throughput sequencing[J]. G3(Bethesda),2016,6(3):755-765.doi:10.1534/g3.115.025981.
[26] Yang C H,Li D Y,Mao D H,Liu X,Ji C J,Li X B,Zhao X F,Cheng Z K,Chen C Y,Zhu L H.Overexpression of microRNA319 impacts leaf morphogenesis and leads to enhanced cold tolerance in rice(Oryza sativa L.)[J]. Plant Cell & Environment,2013,36(12):2207-2218.doi:10.1111/pce.12130.
[27] Liu H H,Tian X,Li Y J,Wu C A,Zheng C C.Microarray-based analysis of stress-regulated microRNAs in Arabidopsis thaliana[J]. RNA,2008,14(5):836-843.doi:10.1261/rna.895308.
[28] Zhou X F,Wang G D,Sutoh K,Zhu J K,Zhang W X.Identification of cold-inducible microRNAs in plants by transcriptome analysis[J]. Biochimica et Biophysica Acta,2008,1779(11):780-788.doi:10.1016/j.bbagrm.2008.04.005.
[29] Hafner M,Landgraf P,Ludwig J,Rice A,Ojo T,Lin C,Holoch D,Lim C,Tuschl T.Identification of microRNAs and other small regulatory RNAs using cDNA library sequencing[J]. Methods, 2008,44(1):3-12.doi:10.1016/j.ymeth.2007.09.009.
[30] Quinlan A R,Hall I M.BEDTools:a flexible suite of utilities for comparing genomic features[J]. Bioinformatics,2010,26(6):841-842.doi:10.1093/bioinformatics/btq033.
[31] Wang L K,Feng Z X,Wang X,Wang X W,Zhang X G.DEGseq:an R package for identifying differentially expressed genes from RNA-seq data[J]. Bioinformatics,2010,26(1):136-138.doi:10.1093/bioinformatics/btp612.
[32] Wu H J,Ma Y K,Chen T,Wang M,Wang X J.PsRobot:a web-based plant small RNA meta-analysis toolbox[J]. Nucleic Acids Research,2012,40(W1):22-28.doi:10.1093/nar/gks554.
[33] Budak H,Khan Z,Kantar M.History and current status of wheat miRNAs using next-generation sequencing and their roles in development and stress[J]. Briefings in Functional Genomics,2015,14(3):189-198.doi:10.1093/bfgp/elu021.
[34] Cao X,Wu Z,Jiang F L,Zhou R,Yang Z.Identification of chilling stress-responsive tomato microRNAs and their target genes by high-throughput sequencing and degradome analysis[J]. BMC Genomics,2014,15(1):1130.doi:10.1186/1471-2164-15-1130.
[35] Xu S C,Liu N,Mao W H,Hu Q Z,Wang G F,Gong Y M.Identification of chilling-responsive microRNAs and their targets in vegetable soybean(Glycine max L.)[J]. Scientific Reports,2016,6:26619.doi:10.1038/srep26619.
[36] Sun X M,Fan G T,Su L Y,Wang W J,Liang Z C,Li S H,Xin H P.Identification of cold-inducible microRNAs in grapevine[J]. Frontiers in Plant Science,2015,46:595.doi:10.3389/fpls.2015.00595.
[37] Lü D K,Bai X,Li Y,Ding X D,Ge Y,Cai H,Ji W,Wu N,Zhu Y M.Profiling of cold-stress-responsive miRNAs in rice by microarrays[J]. Gene,2010,459(1-2):39-47.doi:10.1016/j.gene.2010.03.011.
[38] Song G Q,Zhang R Z, Zhang S J, Li Y L, Gao J, Han X D,Chen M L,Wang J,Li W, Li G Y.Response of microRNAs to cold treatment in the young spikes of common wheat[J]. BMC Genomics,2017,18(1):212.doi:10.1186/s12864-017-3556-2.
[39] Kruszka K,Pacak A,Swida-Barteczka A,Nuc P,Alaba S,Wroblewska Z,Karlowski W,Jarmolowski A,Szweykowska-Kulinska Z.Transcriptionally and post-transcriptionally regulated microRNAs in heat stress response in barley[J]. Journal of Experimental Botany,2014,65(20):6123-6135.doi:10.1093/jxb/eru353.
[40] Wang Q S,Liu N,Yang X Y,Tu L L,Zhang X L.Small RNA-mediated responses to low-and high-temperature stresses in cotton[J]. Scientific Reports,2016,6(1):35558.doi:10.1038/srep35558.
[41] Chen L,Luan Y S,Zhai J M.Sp-miR396a-5p acts as a stress-responsive genes regulator by conferring tolerance to abiotic stresses and susceptibility to Phytophthora nicotianae infection in transgenic tobacco[J]. Plant Cell Reports,2015,34(12):2013-2025.doi:10.1007/s00299-015-1847-0.
[42] Llave C,Xie Z X,Kasschau K D,Carrington J C. Cleavage of Scarecrow-like mRNA targets directed by a class of Arabidopsis miRNA[J]. Science,2002,297(5589):2053-2056.doi:10.1126/science.1076311.
[43] Huang W,Peng S Y,Xian Z Q,Lin D B,Hu G J,Yang L,Ren M Z,Li Z G.Overexpression of a tomato miR171 target gene SlGRAS24 impacts multiple agronomical traits via regulating gibberellin and auxin homeostasis[J]. Plant Biotechnology Journal,2017,15(4):472-488.doi:10.1111/pbi.12646.
[44] Ding X L,Li J J,Zhang H,He T T,Han S H,Li Y W,Yang S P,Gai J Y.Identification of miRNAs and their targets by high-throughput sequencing and degradome analysis in cytoplasmic male-sterile line NJCMS1A and its maintainer NJCMS1B of soybean[J]. BMC Genomics,2016,17(1):24.doi:10.1186/s12864-015-2352-0.
[45] Wang B,Sun Y F,Song N,Wei J P,Wang X J,Feng H,Yin Z Y,Kang Z S.MicroRNAs involving in cold,wounding and salt stresses in Triticum aestivum L.[J]. Plant Physiol Biochem,2014,80:90-96.doi:10.1016/j.plaphy.2014.03.020.
[46] Zhang Y,Zhu X J,Chen X,Song C N,Zou Z W,Wang Y H,Wang M L,Fang W P,Li X H.Identification and characterization of cold-responsive microRNAs in tea plant(Camellia sinensis) and their targets using high-throughput sequencing and degradome analysis[J]. BMC Plant Biology,2014,14(1):271.doi:10.1186/s12870-014-0271-x.
[47] 曹杰.小麦穗发育及过氧化氢胁迫相关miRNA的克隆与鉴定[D].北京:中国农业大学,2019.
Cao J.Characteriztion of miRNAs associated with spike development and H₂O₂response in wheat[D].Beijing:China Agricultural University,2019.
[48] Xie K B,Wu C Q,Xiong L Z.Genomic organization,differential expression,and interaction of SQUAMOSA promoter-Binding-like transcription factors and microRNA 156 in rice[J]. Plant Physiology,2006,142(1):280-293.doi:10.1104/pp.106.084475.
[49] Liu J,Cheng X L,Liu P,Sun J Q.miR156-targeted SBP-Box transcription factors interact with DWARF53 to regulate TEOSINTE BRANCHED1 and BARREN STALK1 expression in bread wheat[J]. Plant Physiology,2017,174(3):1931-1948.doi:10.1104/pp.17.00445.
[50] Sunkar R, Kapoor A,Zhu J K.Posttranscriptional induction of two Cu/Zn superoxide dismutase genes in Arabidopsis is mediated by downregulation of miR398 and important for oxidative stress tolerance[J]. The Plant Cell,2006,18(8):2051-2065.doi:10.1105/tpc.106.041673.
[51] Seo P J,Mas P.Stressing the role of the plant circadian clock[J]. Trends in Plant Science,2015,20(4):230-237.doi:10.1016/j.tplants.2015.01.001.
[52] Kant P, Gordon M, Kant S, Zolla G, Davydov O, Heimer Y M, Chalifa-Caspi V, Shaked R, Barak S.Functional-genomics-based identification of genes that regulate Arabidopsis responses to multiple abiotic stresses[J]. Plant Cell & Environment,2008,31(6):697-714.doi:10.1111/j.1365-3040.2008.01779.x.
[53] Gould P D, Locke J C W, Larue C, Southern M M, Davis S J, Hanano S, Moyle R, Milich R, Putterill J, Millar A J, Hall A.The molecular basis of temperature compensation in the Arabidopsis circadian clock[J]. The Plant Cell,2006,18(5):1177-1187.doi:10.1105/tpc.105.039990.
[54] Bieniawska Z,Espinoza C,Schlereth A,Sulpice R,Hincha D K,Hannah M A.Disruption of the Arabidopsis circadian clock is responsible for extensive variation in the cold-responsive transcriptome[J]. Plant Physiology,2008,147(1):263-279.doi:10.1104/pp.108.118059.

选择文件类型/文献管理软件名称

选择包含的内容