Mining Myb transcription factors related to wood development in Larix olgensis

Abstract

Larix olgensis A.Henry is a fast-growing tree used for afforestation in northeastern China and has great ecological and economic value. For studying developmental genes in the xylem of this species, we investigated the Myb transcription factor family, one of the largest families of transcription factors in plants, which plays an important role in the regulation of lignification in plant secondary walls. By sequencing a L. olgensis cDNA library using the Illumina HiSeq2500 high-throughput sequencing platform, we obtained 58,683 unigene sequences, of which 16,554 unigenes were longer than 1000 bp, accounting for 28.2% of the total database. The alignment of these genes with the GO, COG, KEGG, Swiss-Prot and NR databases resulted in annotated 29,350 unigenes. We obtained a total of 1460 differentially expressed genes, of which 453 were upregulated and 1007 were downregulated at the two developmental stages analyzed. The gene annotations showed a wide range of biological functions and metabolic pathways. The 10 Myb transcription factors that were obtained from the differentially expressed genes were analyzed by real-time quantitative PCR (qRT-PCR). The results showed that four Myb transcription factors may be associated with xylem development in L. olgensis. Due to the large genome size of conifers, genomics research on these species has lagged behind that for other plant groups. Our data provide the basis for further studies on xylem development in L. olgensis.

Keywords

Larix olgensis / RNA-seq / Myb transcription factors / Differentially expressed genes

References

[1] Ahuja MR, Neale DB. Evolution of genome size in conifers. Silvae Genet, 2005, 54: 126-137.

[2] Allona I, Quinn M, Shoop E, Swope K, Cyr SS, Carlis J, Riedl J, Retzel E, Campbell MM, Sederoff R, Whetten RW. Analysis of xylem formation in pine by cDNA sequencing. Proc Natl Acad Sci USA, 1998, 5: 9693-9698.

[3] Altschul SF, Madden TL, Schäffer AA, Zhang JH, Zhang Z, Miller W, Lipman DJ. Gapped BLAST and PSI BLAST: a new generation of protein database search programs. Nucleic Acids Res, 1997, 25: 3389-3402.

[4] Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol, 2010, 11: 106.

[5] Apweiler R, Bairoch A, Wu CH, BarkerWC Boeckmann B, Ferro S, Gasteiger E, Huang HZ, Lopez R, Magrane M, Martin MJ, Natale DA, O’Donovan C, Redaschi N, Yeh LS. UniProt: the universal protein knowledgebase. Nucleic Acids Res, 2004, 32: 115-119.

[6] Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry M, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Tarver LI, Kasarskis A, Lewis S, Matese JC, Richiardson JE, Ringwald M, Rubin GM, Sherlock G. Gene ontology: tool for the unification of biology. Nat Genet, 2000, 25: 25-29.

[7] Bedon F, Grima-Pettenati J, Mackay J. Conifer R2R3-MYB transcription factors: sequence analyses and gene expression in wood-forming tissues of white spruce (Picea glauca). BMC Plant Biol, 2007, 7: 17.

[8] Bhargava A, Mansfield SD, Hall HC, Douglas CJ, Ellis BE. MYB75 functions in regulation of secondary cell wall formation in the Arabidopsis inflorescence stem. Plant Physiol, 2010, 154: 1428.

[9] Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL, Vandesompele J, Wittwer CT. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem, 2009, 55: 611-622.

[10] Chai GH, Wang ZG, Tang XF, Yu L, Qi G, Wang D, Yan XF, Kong YZ, Zhou GK. R2R3-MYB gene pairs in populus: evolution and contribution to secondary wall formation and flowering time. J Exp Bot, 2014, 65: 4255-4269.

[11] Chano V, Collada C, Soto A. Transcriptomic analysis of wound xylem formation in Pinus canariensis. BMC Plant Biol, 2017, 17: 234.

[12] Chen XR, Huang SX, Zhang Y, Sheng GL, Zhang BY, Li QY, Zhu F, Xu JY. Transcription profiling and identification of infection-related genes in Phytophthora cactorum. Mol Genet Genom, 2018, 293: 541-555.

[13] Chiaiese P, Ruotolo G, Matteo AD, De Santo AV, De Marco A, Filippone E. Cloning and expression analysis of kenaf (Hibiscus cannabinus L.) major lignin and cellulose biosynthesis gene sequences and polymer quantification during plant development. Ind Crops Prod, 2011, 34: 1072-1078.

[14] Correll MJ, Pyle TP, Millar KD, Yao J, Edelmann RE, Kiss JZ. Transcriptome analyses of Arabidopsis thaliana seedlings grown in space: implications for gravity-responsive genes. Planta, 2013, 238: 519-533.

[15] Deng YY, Li JQ, Wu SF, Zhu YP, Chen YW, He FC. Integrated nr database in protein annotation system and its localization. Comput Eng, 2006, 32: 71-72.

[16] Desprez T, Juraniec M, Crowell EF, Jouy H, Pochylova Z, Parcy F, Höfte H, Gonneau M, Vernhettes S. Organization of cellulose synthase complexes involved in primary cell wall synthesis in Arabidopsis thaliana. Proc Natl Acad Sci USA, 2007, 104: 15572-15577.

[17] Dubos C, Stracke R, Grotewold E, Weisshaar B, Martin C, Lepiniec L. Myb transcription factors in Arabidopsis. Trends Plant Sci, 2010, 15: 573-581.

[18] Feng SQ, Wang YL, Yang S, Xu YT, Chen XS. Anthocyanin biosynthesis in pears is regulated by a R2R3-MYB transcription factor PyMYB10. Planta, 2010, 232: 245-255.

[19] Fernández-Pozo N, Canales J, Guerrero-Fernández D, Villalobos DP, Díaz-Moreno SM, Bautista R, Monterroso AF, Guevara MÁ, Perdiguero P, Collada C, Cervera MT, Soto Á, Ricardo Ordás, Cantón FR, Avila C, Cánovas FM, Claros MG. EuroPineDB: a high-coverage web database for maritime pine transcriptome. BMC Genom, 2011, 12: 366.

[20] Gibbs DJ, Voß U, Harding SA, Fannon J, Moody LA, Yamada E, Swarup K, Nibau C, Bassel GW, Choudhary A, Lavenus J, Bradshaw SJ, Stekel DJ, Bennett MJ, Coates J. AtMYB93 is a novel negative regulator of lateral root development in Arabidopsis. New Phytol, 2014, 203: 1194-1207.

[21] Goicoechea M, Lacombe E, Legay S, Mihaljevic S, Rech P, Jauneau A, Lapierre C, Pollet B, Verhaegen D, Chaubet-Gigot N, Grima-Pettenati J. EgMYB2, a new transcriptional activator from Eucalyptus xylem, regulates secondary cell wall formation and lignin biosynthesis. Plant J, 2005, 43: 553-567.

[22] Gomez R, Gonzalez J, Herrera R, Herrera R, Ramos P. MYB transcription factors and a putative flavonoid transporter ABCC-like are differentially expressed in radiata pine seedlings exposed to inclination. J Plant Growth Regul, 2018, 37: 64-75.

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

[24] Guo F, Guo YF, Wang P, Wang Y, Ni DJ. Transcriptional profling of catechins biosynthesis genes during tea plant leaf development. Planta, 2017, 246: 1139-1152.

[25] Guo HY, Wang YC, Wang LQ, Hu P, Wang YM, Jia YY, Zhang CR, Zhang Y, Zhang YM, Wang C, Yang CP. Expression of the MYB transcription factor gene BplMYB46 affects abiotic stress tolerance and secondary cell wall deposition in Betula platyphylla. Plant Biotechnol J, 2017, 15: 107-121.

[26] Guo XL, Yu C, Wan HH, Zhen N, Xu TL, Tan JR, Pan HT, Zhang QX. Transcriptome of the floral transition in Rosa chinensis ‘Old Blush’. BMC Genom, 2017, 18: 199.

[27] Hu M, Tan BY, Li XM, Tu JF, Yang FC, Zhu HY, Qin ZQ. Isolation and characterization of a cinnamoyl CoA reductase gene (CCR) in pear (Pyrus pyrifolia). Agric Sci Technol, 2014, 15: 926-932.

[28] Huang HH, Xu LL, Tong ZK, Lin EP, Liu QP, Cheng LJ, Zhu MY. De novo characterization of the Chinese fir (Cunninghamia lanceolata) transcriptome and analysis of candidate genes involved in cellulose and lignin biosynthesis. BMC Genom, 2012, 13: 648.

[29] Humphreys JM, Chapple C. Rewriting the lignin roadmap. Curr Opin Plant Biol, 2002, 5: 224-229.

[30] Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M. The KEGG resource for deciphering the genome. Nucleic Acids Res, 2004, 32: 277-280.

[31] Kasirajan L, Charles S, Aruchamy K. Cloning and expression analysis of 4-coumarate coa ligase (4CL) gene in sugarcane clones varying in lignin content. Proc Natl Acad Sci India, 2017, 4: 1-6.

[32] Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol, 2009, 10: 25.

[33] Lee MM, Schiefelbein J. Cell pattern in the arabidopsis root epidermis determined by lateral inhibition with feedback. Plant Cell, 2002, 14: 611-618.

[34] Legay S, Sivadon P, Blervacq AS, Pavy N, Baghdady A, Tremblay L, Levasseur C, Ladouce N, Lapierre C, Séguin A, Hawkins S, Mackay J, Grima-Pettenati J. EgMYB1, an R2R3 MYB transcription factor from eucalyptus negatively regulates secondary cell wall formation in Arabidopsis and poplar. New Phytol, 2010, 188: 774-786.

[35] Li XG, Harry XW, Shannon KD, Simon GS. Generation and analysis of expressed sequence tags from six developing xylem libraries in Pinus radiata D. Don. BMC Genom, 2009, 10: 41.

[36] Li XG, Harry XW, Simon GS. Transcriptome profiling of Pinus radiata juvenile wood with contrasting stiffness identifies putative candidate genes involved in microfibril orientation and cell wall mechanics. BMC Genom, 2011, 12: 480.

[37] Li CF, Wang XQ, Lu WX, Liu R, Tian Q, Sun YM, Luo KM. A poplar R2R3-MYB transcription factor, PtrMYB152, is involved in regulation of lignin biosynthesis during secondary cell wall formation. Plant Cell Tissue Organ Cult, 2014, 119: 553-563.

[38] Li CF, Wang XQ, Ran LY, Tian QY, Fan D, Luo K. PtoMYB92 is a transcriptional activator of the lignin biosynthetic pathway during secondary cell wall formation in Populus tomentosa. Plant Cell Physiol, 2015, 56: 2436-2446.

[39] Li WF, Yang WH, Zhang SG, Han SY, Qi LW. Transcriptome analysis provides insights into wood formation during larch tree aging. Tree Genet Genomes, 2017, 13: 19.

[40] Liu YX, Zou DM, Wu BS, Lin DH, Zhang ZH, Wu JC. Cloning and expression analysis of a CCoAOMT, homolog in loquat fruit in response to low-temperature storage. Postharvest Biol Technol, 2015, 105: 45-50.

[41] Machado A, Wu YR, Yang YM, Llewellyn DJ, Dennis ES. The MYB transcription factor GhMYB25 regulates early fibre and trichome development. Plant J, 2009, 59: 52-62.

[42] McCarthy RL, Zhong RQ, Ye ZH. MYB83 is a direct target of SND1 and acts redundantly with MYB46 in the regulation of secondary cell wall biosynthesis in Arabidopsis. Plant Cell Physiol, 2009, 50: 1950-1964.

[43] McCarthy RL, Zhong RQ, Fowler S, Lyskowski D, Piyasena H, Carleton K, Spicer C, Ye ZH. The poplar MYB transcription factors, PtrMYB3 and PtrMYB20, are involved in the regulation of secondary wall biosynthesis. Plant Cell Physiol, 2010, 51: 1084-1090.

[44] Meng XX, Xu F, Song QL, Ye JB, Liao YL, Zhang WW. Isolation, characterization and functional analysis of a novel 3-hydroxy-3-methylglutaryl-coenzyme A synthase gene (GbHMGS2) from Ginkgo biloba. Acta Physiol Plant, 2018, 40: 72.

[45] Mishima K, Fujiwara T, Iki T, Kuroda K, Yamashita K, Tamura M, Fujisawa Y, Watanabe A. Transcriptome sequencing and profiling of expressed genes in cambial zone and differentiating xylem of Japanese cedar (Cryptomeria japonica). BMC Genom, 2014, 15: 219.

[46] Patzlaff A, McInnis S, Courtenay A, Surman C, Newman LJ, Smith C, Bevan MW, Mansfield S, Whetten RW, Sederoff RR, Campbell MM. Characterisation of a pine MYB that regulates lignification. Plant J, 2003, 36: 743-754.

[47] Patzlaff A, Newman LJ, Dubos C, Whetten RW, Smith C, McInnis S, Bevan MW, Sederoff RR, Campbell MM. Characterisation of Pt MYB1, an R2R3-MYB from pine xylem. Plant Mol Biol, 2003, 53: 597-608.

[48] Qiu LL, Jiang B, Fang J, Shen Y, Fang ZX, Saravana Kumar RM, Yi KK, Shen CJ, Yan DL, Zheng BS. Analysis of transcriptome in hickory (Carya cathayensis), and uncover the dynamics in the hormonal signaling pathway during graft process. BMC Genom, 2016, 17: 935.

[49] Sekeroglu A, Jacobsen JM, Jansenolesen I, Gupta S, Sheykhzade M, Olesen J, Bhatt DK. Effect of PGD2 on middle meningeal artery and mRNA expression profile of L-PGD2 synthase and DP receptors in trigeminovascular system and other pain processing structures in rat brain. Pharmacol Rep, 2016, 69: 50-56.

[50] Seo PJ, Xiang FN, Qiao M, Park JY, Lee YN, Kim SG, Lee YH, Park WJ, Park CM. The MYB96 transcription factor mediates abscisic acid signaling during drought stress response in Arabidopsis. Plant Physiol, 2009, 151: 275-289.

[51] Tatusov RL, Galperin MY, Natale DA, Koonin EV. The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res, 2000, 28: 33-36.

[52] Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol, 2010, 28: 511-515.

[53] Vanholme R, Morreel K, Ralph J, Boerjan W. Lignin engineering. Curr Opin Plant Biol, 2008, 11: 278-285.

[54] Vanholme R, Demedts B, Morreel K, Ralph J, Boerjan W. Lignin biosynthesis and structure. Plant Physiol, 2010, 153: 895-905.

[55] Visser EA, Wegrzyn JL, Steenkmap ET, Myburg AA, Naidoo S. Combined de novo and genome guided assembly and annotation of the Pinus patula juvenile shoot transcriptome. BMC Genom, 2015, 16: 1057.

[56] Wong MM, Cannon CH, Wickneswari R. Identification of lignin genes and regulatory sequences involved in secondary cell wall formation in Acacia auriculiformis and Acacia mangium via de novo transcriptome sequencing. BMC Genom, 2011, 12: 342.

[57] Yang L, Zhao X, Ran LY, Li CF, Di F, Keming L. PtoMYB156 is involved in negative regulation of phenylpropanoid metabolism and secondary cell wall biosynthesis during wood formation in poplar. Sci Rep, 2017, 7: 41209.

[58] Yu SW, Huang AN, Li J, Gao L, Feng YN, Pemberton E, Chen CL. OsNAC45, plays complex roles by mediating POD activity and the expression of development-related genes under various abiotic stresses in rice root. Plant Growth Regul, 2018, 84: 519-531.

[59] Zhang B, Hu ZZ, Zhang YJ, Li YL, Zhou S, Chen GF. A putative functional myb transcription factor induced by low temperature regulates anthocyanin biosynthesis in purple kale (Brassica oleracea, var. acephala, f. tricolor). Plant Cell Rep, 2012, 31: 281-289.

[60] Zhang Y, Zhang SG, Han SY, Li XM, Qi LW. Transcriptome profiling and in silico analysis of somatic embryos in Japanese larch (Larix leptolepis). Plant Cell Rep, 2012, 31: 1637-1657.

[61] Zhong R, Richardson EA, Ye ZH. Two NAC domain transcription factors, SND1 and NST1, function redundantly in regulation of secondary wall synthesis in fibers of Arabidopsis. Planta, 2007, 225: 1603-1611.

[62] Zhong RQ, Richardson EA, Ye ZH. The MYB46 transcription factor is a direct target of SND1 and regulates secondary wall biosynthesis in Arabidopsis. Plant Cell, 2007, 19: 2776-2792.

[63] Zhong RQ, Lee CH, Zhou JL, McCarthy Ye ZH. A battery of transcription factors involved in the regulation of secondary cell wall biosynthesis in Arabidopsis. Plant Cell, 2008, 20: 2763-2782.

[64] Zhou JL, Lee CH, Zhong RQ, Ye ZH. MYB58 and MYB63 are transcriptional activators of the lignin biosynthetic pathway during secondary cell wall formation in Arabidopsis. Plant Cell, 2009, 21: 248-266.