Analysis of differentially expressed genes in relation to tiller development and plant height based on transcriptomic sequencing of two tall fescue cultivars

doi:10.11686/cyxb2020486

Abstract

Abstract:

Tiller number and plant height are important agronomic traits for grasses. Identification of genes controlling tiller development and plant height is of great significance for molecular assisted breeding of forage and turfgrass. In this study， two varieties of tall fescue （Festuca arundinacea） 'Kentucky-31' （K31） and 'Regenerate' were used to construct transcriptome maps and explore the differentially expressed genes （DEGs） regulating shoot growth and development. Transcriptomic sequencing was performed based on Illumina HiSeq 2500×Miseq 300， the data obtained were assembled de novo and a total of 77872 unigenes were obtained. The obtained unigenes were then blasted against the non-redundant protein （NR）， universal protein （Uniprot）， gene ontology （GO）， kyoto encyclopedia of genes and genomes （KEGG） and cluster of orthologous groups （COG） databases， and 59927， 40213， 44447， 15146 and 13767 unigenes， respectively， were successfully annotated. A total of 1573 up-regulated DEGs and 1441 down-regulated DEGs were detected in 'Regenerate' when compared to 'K31'. GO enrichment analysis found that DEGs are mainly enriched in biological processes such as cell and cell component and macromolecular complex assembly. A total of 42 differentially expressed transcription factors were noted among the identified DEGs， mainly belonging to teosinte branched 1， cycloidea， proliferating cell factors （TCP）， WRKY and auxin response factor （ARF） families. DEGs relating to eight plant hormone categories were also annotated， including auxin， cytokinin， abscisic acid， gibberellin， ethylene， brassinosteroids， salicylic acid and jasmonic acid. The expression pattern of DEGs verified by real-time PCR was consistent with RNA-Seq， indicating the RNA sequencing results were accurate. This study has advanced the transcriptome knowledge for tall fescue， obtained a preliminary identification of the candidate genes controlling plant height and tiller development， and provided theoretical support for further research on gene function and molecular breeding.

Key words: tall fescue, RNA-seq, differentially expressed genes, transcription factors, plant hormones

Zhi-min YANG, Rui XING, Yun-jia DING, Li-li ZHUANG. Analysis of differentially expressed genes in relation to tiller development and plant height based on transcriptomic sequencing of two tall fescue cultivars[J]. Acta Prataculturae Sinica, 2022, 31(1): 145-163.

Figures/Tables 14

References 58

1	Janssen B J， Drummond R S， Snowden K C. Regulation of axillary shoot development. Current Opinion in Plant Biology， 2014， 17： 28-35.
2	Teo Z W N， Song S， Wang Y Q， et al. New insights into the regulation of inflorescence architecture. Trends in Plant Science， 2014， 19（3）： 158-165.
3	Rameau C， Bertheloot J， Leduc N， et al. Multiple pathways regulate shoot branching. Frontiers in Plant Science， 2015， 5： 741.
4	Karin S， Thomas S， Mathias R， et al. The Lateral suppressor （Ls） gene of tomato encodes a new member of the VHIID protein family. Proceedings of the National Academy of Sciences of the United States of America， 1999， 96（1）： 290-295.
5	Li X Y， Qian Q， Fu Z M， et al. Control of tillering in rice. Nature， 2003， 422（6932）： 618.
6	Thomas G， Oliver C， Elisabeth S， et al. Molecular analysis of the LATERAL SUPPRESSOR gene in Arabidopsis reveals a conserved control mechanism for axillary meristem formation. Genes & Development， 2003， 17（9）： 1175-1187.
7	Raman S， Greb T， Peaucelle A， et al. Interplay of miR164， CUP-SHAPED COTYLEDON genes and LATERAL SUPPRESSOR controls axillary meristem formation in Arabidopsis thaliana. Plant Journal， 2008， 55（1）： 65-76.
8	Thomas K， Jessica A， Thomas M， et al. Arabidopsis regulator of axillary meristems1 controls a leaf axil stem cell niche and modulates vegetative development. The Plant Cell， 2006， 18（3）： 598-611.
9	Aguilar-Martínez J A， Poza-Carrión C， Cubas P. Arabidopsis BRANCHED1 acts as an integrator of branching signals within axillary buds. The Plant Cell， 2007， 19（2）： 458-472.
10	Minakuchi K， Kameoka H， Yasuno N， et al. FINE CULM1 （FC1） works downstream of strigolactones to inhibit the outgrowth of axillary buds in rice. Plant and Cell Physiology， 2010， 51（7）： 1127-1135.
11	Nils B， Alexandre de S G， Jean-Paul P， et al. The pea TCP transcription factor PsBRC1 acts downstream of strigolactones to control shoot branching. Plant Physiology， 2012， 158（1）： 225-238.
12	Shi P B， Guy K M， Wu W F， et al. Genome-wide identification and expression analysis of the ClTCP transcription factors in Citrullus lanatus. BMC Plant Biology， 2016， 16： 85-97.
13	Lee D J， Zeevaart J A. Molecular cloning of GA 2-oxidase3 from spinach and its ectopic expression in Nicotiana sylvestris. Plant Physiology， 2005， 138（1）： 243-254.
14	Schomburg F M. Overexpression of a novel class of gibberellin 2-oxidases decreases gibberellin levels and creates dwarf plants. Plant Cell， 2003， 15（1）： 151-163.
15	Guo Y L， Zou S H， Wang H P， et al. Cloning and preliminary functional analysis of three gibberellin 2-oxidase genes in petunia （Petunia hybrida）. Journal of Agricultural Biotechnology， 2013， 21（8）： 940-948.
	郭余龙，邹世慧，王会平，等. 矮牵牛3个赤霉素2-氧化酶基因的克隆及初步功能分析. 农业生物技术学报， 2013， 21（8）： 940-948.
16	Wang Y， Deng D. Molecular basis and evolutionary pattern of GA-GID1-DELLA regulatory module. Molecular Genetics & Genomics， 2014， 289（1）： 1-9.
17	Woodward A W， Bartel B. Auxin： Regulation， action， and interaction. Annals of Botany， 2005， 95（5）： 707-735.
18	Chung Y， Maharjan P M， Lee O， et al. Auxin stimulates DWARF4 expression and brassinosteroid biosynthesis in Arabidopsis. Plant Journal for Cell & Molecular Biology， 2011， 66（4）： 564-578.
19	Salchert K， Bhalerao R， Koncz-Kalman Z， et al. Control of cell elongation and stress responses by steroid hormones and carbon catabolic repression in plants. Philosophical Transactions of the Royal Society of London， 1998， 353（1374）： 1517-1520.
20	Hong Z， Ueguchi-Tanaka M， Shimizu-Sato S， et al. Loss-of-function of a rice brassinosteriod biosynthetic enzyme， C-6 oxidase， prevents the organized arrangement and polar elongation of cells in the leaves and stem. Plant Journal， 2002， 32（4）： 495-508.
21	Mian M A R， Zhang Y， Wang Z Y， et al. Analysis of tall fescue ESTs representing different abiotic stresses， tissue types and developmental stages. BMC Plant Biology， 2008， 8（1）： 27.
22	Hu T， Sun X， Zhang X， et al. An RNA sequencing transcriptome analysis of the high-temperature stressed tall fescue reveals novel insights into plant thermotolerance. BMC Genomics， 2014， 15（1）： 1147.
23	Su Y L， Liu M J， Li H F. Research progress on MADS-box gene in rice （Oryza sativa）. Journal of Henan Agricultural Sciences， 2016， 45（9）： 1-7.
	苏亚丽，刘梦佳，李海峰. 水稻MADS-box基因研究进展. 河南农业科学， 2016， 45（9）： 1-7.
24	Guo S Y， Xu Y Y， Liu H H， et al. The interaction between OsMADS57 and OsTB1 modulates rice tillering via DWARF14. Nature Communications， 2013， 4（1）： 1566.
25	Fabio F， Veronica G， Nilla P， et al. The rice StMADS11-like genes OsMADS22 and OsMADS47 cause floral reversions in Arabidopsis without complementing the svp and agl24 mutants. Journal of Experimental Botany， 2008， 59（8）： 2181-2190.
26	Sakurai A， Yokota T， Clouse S D. Brassinosteroids： Steroidal plant hormones. Tokyo： Springer-Verlag， 1999： 21-45.
27	Kamuro Y， Takatsuto S. Capability for and problems of practical uses of brassinosteroids. Acs Symposium Series， 1991， 474： 292-297.
28	Cano-Delgado A， Yin Y H， Yu C， et al. BRL1 and BRL3 are novel brassinosteroid receptors that function in vascular differentiation in Arabidopsis. Development， 2004， 131（21）： 5341-5351.
29	Martín-Trillo M， Cubas P. TCP genes： A family snapshot ten years later. Trends in Plant Science， 2010， 15（1）： 31-39.
30	Li C X， Potuschak T， Colón-Carmona A， et al. Arabidopsis TCP20 links regulation of growth and cell division control pathways. Proceedings of the National Academy of Sciences of the United States of America， 2005， 102（36）： 12978-12983.
31	Steiner E， Efroni I， Gopalraj M， et al. The ArabidopsisO-linked N-acetylglucosamine transferase SPINDLY interacts with class I TCPs to facilitate cytokinin responses in leaves and flowers. Plant Cell， 2012， 24（1）： 96-108.
32	Daviere J M， Wild M， Regnault T， et al. Class I TCP-DELLA interactions in inflorescence shoot apex determine plant height. Current Biology， 2014， 24（16）： 1923-1928.
33	Luo X， Sun X L， Liu B H， et al. Ectopic expression of a WRKY homolog from Glycine soja alters flowering time in Arabidopsis. PLoS One， 2013， 8（8）： e73295.
34	Tian X J， Li X F， Zhou W J， et al. Transcription factor OsWRKY53 positively regulates brassinosteroid signaling and plant architecture. Plant Physiology， 2017， 175（3）： 1337-1349.
35	Zhang C Q， Xu Y， Lu Y， et al. The WRKY transcription factor OsWRKY78 regulates stem elongation and seed development in rice. Planta， 2011， 234（3）： 541-554.
36	Song Y， Chen L G， Zhang L P， et al. Overexpression of OsWRKY72 gene interferes in the abscisic acid signal and auxin transport pathway of Arabidopsis. Journal of Biosciences， 2010， 35（3）： 459-471.
37	Yu Y C. The molecular mechanism of WRKY71 transcription factor in regulating flower and branch development of Arabidopsis thialiana. Jinan： Shandong University， 2011.
	于延冲. 拟南芥转录因子WRKY71对花和分枝发育的调控机制研究. 济南：山东大学， 2011.
38	Vanderauwera S， Vandenbroucke K， Inze A， et al. AtWRKY15 perturbation abolishes the mitochondrial stress response that steers osmotic stress tolerance in Arabidopsis. Proceedings of the National Academy of Sciences， 2012， 109（49）： 20113-20118.
39	Yu F， Huaxia Y， Lu W， et al. GhWRKY15， a member of the WRKY transcription factor family identified from cotton （Gossypium hirsutum）， is involved in disease resistance and plant development. BMC Plant Biology， 2012， 12（1）： 144.
40	Cai Y H， Chen X J， Xie K， et al. Dlf1， a WRKY transcription factor， is involved in the control of flowering time and plant height in rice. PLoS One， 2014， 9（7）： e102529.
41	Qiu D Y， Xiao J， Xie W B， et al. Rice gene network inferred from expression profiling of plants overexpressing OsWRKY13， a positive regulator of disease resistance. Molecular Plant， 2008， 1（3）： 538-551.
42	Wang H H， Hao J J， Chen X J， et al. Overexpression of rice WRKY89 enhances ultraviolet B tolerance and disease resistance in rice plants. Plant Molecular Biology， 2007， 65（6）： 799-815.
43	Yao C， Finlayson S A. Abscisic acid is a general negative regulator of Arabidopsis axillary bud growth. Plant Physiology， 2015， 169（1）： 611-626.
44	Arney S E， Mitchell D L. The effect of abscisic acid on stem elongation and correlative inhibition. New Phytologist， 1969， 68（4）： 1001-1015.
45	Chatfield S P， Stirnberg P， Forde B G， et al. The hormonal regulation of axillary bud growth in Arabidopsis. Plant Journal， 2000， 24（2）： 159-169.
46	Cline M G， Choonseok O. Reappraisal of the role of abscisic acid and its interaction with auxin in apical dominance. Annals of Botany， 2006， 98（4）： 891-897.
47	Bris M L， Nicole M F， Jacob Y， et al. Regulation of bud dormancy by manipulation of ABA in isolated buds of Rosa hybrida cultured invitro. Functional Plant Biology， 1999， 26（3）： 273-281.
48	Wang R F， Zhang J W， Lv P， et al. Effects of endogenous hormones on tiller development process of different maize varieties. Scientia Agricultura Sinica， 2012， 45（5）： 840-847.
	王如芳，张吉旺，吕鹏，等. 不同类型玉米品种分蘖发生过程中内源激素的作用. 中国农业科学， 2012， 45（5）： 840-847.
49	Domagalska M A， Leyser O. Signal integration in the control of shoot branching. Nature Reviews Molecular Cell Biology， 2011， 12（4）： 211-221.
50	Reddy S K， Finlayson S A. Phytochrome B promotes branching in Arabidopsis by suppressing auxin signaling. Plant Physiology， 2014， 164（3）： 1542-1550.
51	Song Y L， You J， Xiong L Z. Characterization of OsIAA1 gene， a member of rice Aux/IAA family involved in auxin and brassinosteroid hormone responses and plant morphogenesis. Plant Molecular Biology， 2009， 70（3）： 297-309.
52	Fukaki H， Nakao Y， Okushima Y， et al. Tissue-specific expression of stabilized SOLITARY-ROOT/IAA14 alters lateral root development in Arabidopsis. Plant Journal， 2005， 44（3）： 382-395.
53	Sakamoto T， Morinaka Y， Ishiyama K， et al. Genetic manipulation of gibberellin metabolism in transgenic rice. Nature Biotechnology， 2003， 21（8）： 909-913.
54	Shan C， Mei Z L， Duan J L， et al. OsGA2ox5， a gibberellin metabolism enzyme， is involved in plant growth， the root gravity response and salt stress. PLoS One， 2014， 9（1）： e87110.
55	Huang J， Tang D， Shen Y， et al. Activation of gibberellin 2-oxidase 6 decreases active gibberellin levels and creates a dominant semi-dwarf phenotype in rice （Oryza sativa L.）. Journal of Genetics and Genomics， 2010， 37（1）： 23-36.
56	Hong X F， Jiang P Y， Zheng Z S， et al. Effects of gibberellin （GA₃） on tillering control and tillering rate improvement in rice （Oryza sativa L.） tillering stage. Journal of Zhejiang Agricultural Sciences， 1998（1）： 3-5.
	洪晓富，蒋鹏炎，郑寨生，等. 水稻分蘖期喷施赤霉素（GA₃）对控制分蘖和提高成穗率的效果. 浙江农业科学， 1998（1）： 3-5.
57	Barleben L， Ma X Y， Koepke J， et al. Expression， purification， crystallization and preliminary X-ray analysis of strictosidine glucosidase， an enzyme initiating biosynthetic pathways to a unique diversity of indole alkaloid skeletons. Biochimica et Biophysica Acta， 2005， 1747（1）： 89-92.
58	Kovi M R， Zhang Y， Yu S， et al. Candidacy of a chitin-inducible gibberellin-responsive gene for a major locus affecting plant height in rice that is closely linked to green revolution gene sd1. Theoretical & Applied Genetics， 2011， 123（5）： 705-714.

基因编号Unigene ID	正向引物（5′-3′）Forward primers （5′-3′）	反向引物（5′-3′）Reverse primers （5′-3′）
Unigene003394	GGTGGACGAGAAGCAGTACT	CGGGTAGTAGTTGCAGGTGA
Unigene054336	CTCAAGGGCAAGGAATCTGC	ATTCGATCGGGGAGAAAGCA
Unigene005703/4/5	AGAGTGTGATGGCGTTGAGA	GCAGATTCCTTGGCCTTGAG
Unigene005713	TGTCTCAACACCGTCACACT	ATGTTGAGGGAGATGGTGGC
Unigene074113	CAGAACTGGCAGAGCAACTC	GGGAGGAGGAAGACGACTAC
Unigene031989	CATCGTGCCAGTCAACTACAG	AGGTTGAAGTAGGAGTGGCC
Unigene055976	TCCCTCGTCAACTACCCCTA	TTGGGGAGCAGATGACACAT
Unigene063823	GCAAGCACCGATCCAGTATG	CGGACTTGTAGGATTTGGCG
Unigene036563	CTCAACTCCAACCGCTTCTG	GTTGTTGCTGAGGTCGAGC
Unigene067511	CTGCTCCTGTCCACCTGTAA	AGTCAGCGAAGAAGTCACCA
Unigene056675	GCGTCGTCACTGCTTTCTAT	CGTCGATGTACAGGTTGGTC
Unigene031117	CTCTCACTTGCCGTCTTTGG	GACGAGGCAGTGTCCATCT
Unigene059338	GAGCTGCCTTTCATCGCTAC	CTCTGGGATGGCAGAAAGGA
Unigene033181	GCTTACTCTCAGCCTCCACT	GGCCTGCAGAAATTGGGAAT
Unigene023003	TTGCTTTAGAAGGCAGGCAC	AGCAGCAAACACACCACTTT
Unigene033817	TGAAGCAGATGTGGAGCAGA	GCCTCCGCTTTATTGCTTGA
Unigene061235	AAGAACTGCATGCCGTTGTT	ATTCTTGGAGGCGTTCATGC
Unigene074427	CAACAAGGGTGCTGAGATCG	GGCGATGCTACCTTGTTCAG
Unigene037483	CTCGATAAGCTGTGGGACGA	CTAGTGCTGTGTTGATGGCC
Unigene003016	ATGGCATGTTCAACGGTGTC	GAGTATCTGTCCGTGGCTGA
Unigene066660	TCATCCCAGGTGAACACTCC	GCGCTGTACAAGTTCGGTAG
Unigene074432	ATCCGTCTCAGGTCCAAGTG	GCGAAATACTGTCCGCCAAT
Unigene022674	GACGGACGGAATCAAGAGGA	CATACCATTCGAAGTGCCGG
Unigene054725	GTAGCACCTCCGTCCCTATC	GGACTTGATCGAGCTGAGGA

基因编号Unigene ID	正向引物（5′-3′）Forward primers （5′-3′）	反向引物（5′-3′）Reverse primers （5′-3′）
Unigene003394	GGTGGACGAGAAGCAGTACT	CGGGTAGTAGTTGCAGGTGA
Unigene054336	CTCAAGGGCAAGGAATCTGC	ATTCGATCGGGGAGAAAGCA
Unigene005703/4/5	AGAGTGTGATGGCGTTGAGA	GCAGATTCCTTGGCCTTGAG
Unigene005713	TGTCTCAACACCGTCACACT	ATGTTGAGGGAGATGGTGGC
Unigene074113	CAGAACTGGCAGAGCAACTC	GGGAGGAGGAAGACGACTAC
Unigene031989	CATCGTGCCAGTCAACTACAG	AGGTTGAAGTAGGAGTGGCC
Unigene055976	TCCCTCGTCAACTACCCCTA	TTGGGGAGCAGATGACACAT
Unigene063823	GCAAGCACCGATCCAGTATG	CGGACTTGTAGGATTTGGCG
Unigene036563	CTCAACTCCAACCGCTTCTG	GTTGTTGCTGAGGTCGAGC
Unigene067511	CTGCTCCTGTCCACCTGTAA	AGTCAGCGAAGAAGTCACCA
Unigene056675	GCGTCGTCACTGCTTTCTAT	CGTCGATGTACAGGTTGGTC
Unigene031117	CTCTCACTTGCCGTCTTTGG	GACGAGGCAGTGTCCATCT
Unigene059338	GAGCTGCCTTTCATCGCTAC	CTCTGGGATGGCAGAAAGGA
Unigene033181	GCTTACTCTCAGCCTCCACT	GGCCTGCAGAAATTGGGAAT
Unigene023003	TTGCTTTAGAAGGCAGGCAC	AGCAGCAAACACACCACTTT
Unigene033817	TGAAGCAGATGTGGAGCAGA	GCCTCCGCTTTATTGCTTGA
Unigene061235	AAGAACTGCATGCCGTTGTT	ATTCTTGGAGGCGTTCATGC
Unigene074427	CAACAAGGGTGCTGAGATCG	GGCGATGCTACCTTGTTCAG
Unigene037483	CTCGATAAGCTGTGGGACGA	CTAGTGCTGTGTTGATGGCC
Unigene003016	ATGGCATGTTCAACGGTGTC	GAGTATCTGTCCGTGGCTGA
Unigene066660	TCATCCCAGGTGAACACTCC	GCGCTGTACAAGTTCGGTAG
Unigene074432	ATCCGTCTCAGGTCCAAGTG	GCGAAATACTGTCCGCCAAT
Unigene022674	GACGGACGGAATCAAGAGGA	CATACCATTCGAAGTGCCGG
Unigene054725	GTAGCACCTCCGTCCCTATC	GGACTTGATCGAGCTGAGGA

样品Sample	原始数据Raw reads	测序数据Clean reads	原始下机数据Raw bases	测序下机数据Clean bases	测序数据比率 Clean data rate （%）
Regenerate	88475512	88345910	11834878498	11576251799	97.81
K31	76673214	76533900	10284187805	10037407880	97.60

样品Sample	原始数据Raw reads	测序数据Clean reads	原始下机数据Raw bases	测序下机数据Clean bases	测序数据比率 Clean data rate （%）
Regenerate	88475512	88345910	11834878498	11576251799	97.81
K31	76673214	76533900	10284187805	10037407880	97.60

RNA-Seq信息 RNA-Seq information	数据 Data
Unigene 数目 Unigene number	77872
最大unigene bases Max unigene bases （bp）	13402
最小unigene bases Min unigene bases （bp）	200
全部数据长度 Whole dataset lengths （bp）	47586469
平均unigene长度 Average unigene lengths （bp）	611.086
（G+C）/（A+T+G+C）（%）	53.2
N₅₀ （bp）	920
N₉₀ （bp）	262