Введение 4
1 Обзор литературы 7
1.1 Актуальность исследования 7
1.2 Гомеобокс-содержащие гены 8
1.3 Гомеобокс-содержащие гены растений 9
1.4 Классификация гомеобокс-содержащих генов растений ... 10
1.5 Изучение гомеобокс-содержащих генов растений 10
1.6 Биоинформатический анализ 12
1.6.1 Анализ данных РНК секвенирования 12
1.6.2 Данные РНК секвенирования 13
1.6.3 Сборка транскриптома de novo 13
1.6.4 Дифференциальная экспрессия генов 14
2 Материалы и методы 17
2.1 Предобработка данных и de novoсборка транскриптома P.sylvestris 17
2.2 Оценка полноты сборки транскриптома 19
2.3 Отбор гомеобокс-содержащих транскриптов 19
2.4 Анализ дифференциальной экспрессии генов и построение
тепловой карты 20
2.5 Аннотация дифференциально экспрессирующихся транскриптов 22
3 Результаты и обсуждение 23
3.1 Сборка транскриптома и оценка качества 23
3.2 Классификация гомеобокс-содержащих транскриптов 23
3.3 Дифференциальная экспрессия гомеобокс-содержащих генов
в различных тканях Pinus sylvestris 25
Список использованных источников 35
Сосна обыкновенная Pinus sylvestrisшироко распространена по всей территории Российской Федерации и имеет важное экономическое и экологическое значение. P. sylvestrisотносится к группе голосеменных растений, семейство сосновые, рода сосна. На уровне геномики и транскриптомики данная группа изучена значительно меньше, чем группа покрытосеменных растений, поэтому изучение голосеменных представляет особый интерес.
Развитие P. sylvestris y других растений, принадлежащих семейству сосновых, изучается рядом ученых, однако на данный момент имеется сравнительно мало данных о генетических механизмах развития этого вида [1¬3]. Комплексный анализ профиля дифференциальной экспрессии гомеобокс- содержащих генов является важным инструментом для функционального анализа и определения роли данных генов в процессе развития и роста растения.
Актуальность работы связана с тем, что подобный анализ данных РНК-секвенирования P. sylvestrisранее не проводился.
Объектом настоящей работы являются взаимосвязи между уровнями экспрессии гомеобокс-содержащих транскриптов и особенностями морфологического развития у разных видов хвойных.
Предметом данного исследования являются нуклеотидные последовательности, полученные в результате ассемблирования транскриптома P. sylvestrisи обнаруженные в них гомеобокс-содержащие гены, их дифференциальная экспрессия в различных тканях.
Целью данной работы является определение функций гомеобокс- содержащих генов в различных тканях Pinus sylvestris.Иными словами, определить в каких процессах развития и дифференцировки тканей могут принимать участие гомеобокс-содержащие гены у сосны обыкновенной. Многие функции данных генов всё еще неизвестны (особенно у хвойных видов). Полученные данные помогут расширить и углубить понимание процессов развития и морфогенеза данного вида и семейства сосновых в целом. Для достижения поставленной цели был выполнен ряд задач, а именно:
1) Произведен отбор данных РНК-секвенирования различных тканей сосны обыкновенной;
2) Выполнены предобработка данных и оценка качества ридов;
3) Осуществлена de novoсборка транскриптома P. sylvestris;
4) Произведен отбор транскриптов, содержащих гомеобокс-содержащие гены;
5) Выполнен анализ дифференциальной экспрессии генов и их аннотация;
6) Определены функции гомеобокс-содержащих транскриптов.
Работа и связанные с ней исследования докладывалась на следующих конференциях:
1) 8thInternational Work-Conference on Bioinformatics and Biomedical En-gineering (IWBBIO), 2020 г, Гранада, устный доклад;
2) Международная конференция студентов, аспирантов и молодых ученых «Проспект Свободный — 2020», 2020 г, Красноярск, устный доклад,
II место;
3) XI международная конференция «Dynamical Systems Applied to Biology and Natural Sciences» (DSABNS), 2020 г, Тренто, устный доклад;
4) Международная конференция летней школы по биоинформатике, 2020 г, Москва, Институт биоинформатики, устный доклад;
5) 7thInternational Work-Conference on Bioinformatics and Biomedical Engineering (IWBBIO), 2019 г, Гранада, устный доклад;
6) Международная конференция студентов, аспирантов и молодых ученых «Проспект Свободный - 2019», 2019 г, Красноярск, устный доклад, III место;
7) X международная конференция «Dynamical Systems Applied to Biology and Natural Sciences» (DSABNS), 2019 г, Неаполь, устный доклад;
8) 56-я Международная научная студенческая конференция, 2018 г, Новосибирск, устный доклад.
Результаты работы и связанных с ней исследований опубликованы в следующих научных журналах и сборниках научных мероприятий:
• Гусева Т.А., Бирюков В.В. Гомеобокс-содержащие гены в изучении процесса развития сосны обыкновенной Pinus sylvestris// Проспект Свободный - 2020: материалы XVI Международной конференции студентов, аспирантов и молодых ученых. — 2020;
• Guseva T., Biriukov V., Sadovsky M. Role of Homeobox Genes in the Development of Pinus sylvestris// Lecture Notes in Computer Science, Springer Verlag. — 2020. — Vol. 12108 — Pp. 429-437;
• Guseva T., Sadovsky M., Biriukov V. Homeobox Genes: Investigating the Development of Pinus sylvestris (Scots pine) // 11thInternational Conference Dynaamical Systems Applied to Biology and Natural Sciences: Book of Abstracts. — 2020. — Vol. 11. — Pp. 113-115.
• Садовский М.Г., Гусева Т.А., Бирюков В.В. Определение тканеспецифичности в тотальном транскриптоме лиственницы сибирской с помощью условной энтропии // Проспект Свободный - 2019: материалы XV Международной конференции студентов, аспирантов и молодых ученых, посвященной Международному году Периодической таблицы химических элементов Д. И. Менделеева. — 2019. — Стр. 682-684;
• Sadovsky M., Guseva T., Biriukov V. Entropy Approach to Identify Tissue Specificity in Total Transcriptome // 10thInternational Conference Dy- naamical Systems Applied to Biology and Natural Sciences: Book of Ab¬stracts. — 2019. — Vol. 10. — Pp. 98-99;
• Sadovsky M., Guseva T., Biriukov V. Triplet Frequencies Implementation in Total Transcriptome Analysis// Lecture Notes in Computer Science, Sprin¬ger Verlag. — 2019. — Vol. 11465 — Pp. 370-378;
• Гусева Т.А. Исследование особенностей структуры тотального транскриптома лиственницы сибирской на основе статистических свойств контигов // Материалы 56-й Международной научной студенческой конференции (МНСК). — 2018. — Т. 56. — Стр. 10.
В ходе данной работы были выполнены все задачи, и поставленная цель была полностью достигнута:
1. Произведен отбор данных РНК-секвенирования различных тканей сосны обыкновенной из базы NCBI;
2. Выполнена предобработка данных и оценка качества видов;
3. Собран de novoтранскриптом Pinus sylvestris,содержащий 775 502 транскрипта общая длина которого составила 553 398 248 со средним содержанием GC-контента 40.19 %. сборка TrinityСФУ по результатам анализа является более полной, чем сборка TrinityФинляндия
4. Отобрано и аннотировано 243 гомеобокс-содержащих транскрипта;
5. Определены дифференциально экспрессирующиеся гомеобокс-содержащие гены пяти различных тканей сосны обыкновенной и построена тепловая карта (рисунок 5);
6. Предположительно определены функции некоторых гомеобокс-содержащих генов:
• Гены HDG2, PDF2, ROC1, ROC2, ROC8класса HD-ZIP IV участвуют в развитии кутикулы растений, аккумуляции антоцианов и дифференцировки эпидермальных клеток тканей почки и эмбриона; ген ROC8также участвует в дифференцировке эпидермальных клеток хвои;
• Ген WOX3семейства WOX экспрессируется в почках и играет роль в развитии боковых органов сосны обыкновенной;
• Ген WOX4характерен для развития сосудистых тканей Pinus sylvest- ris, в том числе и флоэмы;
• Гены WOX8и WOX9участвуют в развитии эмбриона Pinus sylvest¬ris;
• Ген HAT2экспрессируется в почке и эмбрионе и участвует в процессах, связанных с реакцией на свет и сигнализацией ауксина;
• Ген BLH4участвует в формировании клеточных стенок хвои, почки и флоэмы;
• Во флоэме и почке имеет место взаимодействие классов генов KNOX и BEL и вероятно ген KNAT2участвует в развитии проводящей ткани, гены KN-1, KNAT1, KNAT3имеют аналогичные функции в тканях эмбриона, почки и флоэмы;
• Гены HOX5, HOX20, HOX21, HAT7принадлежат к классу генов HD-ZIP I учувствуют в реакциях организма на стресс и регуляции развития почек сосны обыкновенной;
• Гены ATHB-13, PDF2, ROC1, ROC2, ROC8класса HD-ZIP ТУ участвуют в развитии кутикулы растений, аккумуляции антоцианов и дифференцировки эпидермальных клеток тканей почки и эмбриона.
Подводя итог, можно сказать, что была получена более качественная сборка транскриптома сосны обыкновенной (775 502 транскрипта, N50 1 273 н.о., BUSCO C:85.9 %) по сравнению со сборкой, полученной финскими коллегами (TrinityФинляндия) и впервые описаны и предположительно определены функции гомеобокс-содержащих генов в пяти тканях P. sylvestris.
[1] Transcript profiling for early stages during embryo development in Scots pine / Irene Merino, Malin Abrahamsson, Lieven Sterck et al. // BMC plant biology. — 2016. — Vol. 16, no. 1. — P. 255.
[2] Comparative analysis of transcript abundance in Pinus sylvestris after chal¬lenge with a saprotrophic, pathogenic or mutualistic fungus / Aleksan¬dra Adomas, Gregory Heller, Ake Olson et al. // Tree physiology. — 2008. — Vol. 28, no. 6. — Pp. 885-897.
[3] De Heredia, Unai Lopez. RNA-seq analysis in forest tree species: bioinfor- matic problems and solutions / Unai Lopez De Heredia, Jose Luis Vazquez- Poletti // Tree Genetics & Genomes. — 2016. — Vol. 12, no. 2. — P. 30.
[4] Viola, Ivana L. Structure and evolution of plant homeobox genes / Ivana L Viola, Daniel H Gonzalez // Plant Transcription Factors. — Else¬vier, 2016. — Pp. 101-112.
[5] What Do We Know about Homeodomain-Leucine Zipper I Transcription Factors? Functional and Biotechnological Considerations / Pamela A Ri- bone, Matias Capella, Agustin L Arce, Raquel L Chan // Plant Transcrip¬tion Factors. — Elsevier, 2016. — Pp. 343-356.
[6] The rice WUSCHEL-related homeobox genes are involved in reproduc¬tive organ development, hormone signaling and abiotic stress response / Saifeng Cheng, Yulan Huang, Ning Zhu, Yu Zhao // Gene. — 2014. — Vol. 549, no. 2. — Pp. 266-274.
[7] Analysis of the WUSCHEL-RELATED HOMEOBOX gene family in the conifer Picea abiesreveals extensive conservation as well as dynamic pat¬terns / Harald Hedman, Tianqing Zhu, Sara von Arnold, Joel J Sohlberg // BMC Plant Biology. — 2013. — Vol. 13, no. 1. — P. 89.
[8] Comparative expression pattern analysis of WUSCHEL-related homeobox 2 (WOX2) and WOX8/9 in developing seeds and somatic embryos of the gymnosperm Picea abies/ Joakim Palovaara, Henrik Hallberg, Clau¬dio Stasolla, Inger Hakman // New Phytologist. — 2010. — Vol. 188, no. 1. — Pp. 122-135.
[9] Palovaara, Joakim. Conifer WOX-related homeodomain transcription fac¬tors, developmental consideration and expression dynamic of WOX2 during Picea abies somatic embryogenesis / Joakim Palovaara, Inger Hakman // Plant Molecular Biology. — 2008. — Vol. 66, no. 5. — Pp. 533-549.
[10] Guo, Mei. Down-regulation of a homeodomain-leucine zipper i-class home- obox gene for improved plant performance. — 2018. — US Patent App. 15/599,830.
[11] Genome-wide analysis of homeobox gene family in legumes: identification, gene duplication and expression profiling / Annapurna Bhattacharjee, Ra- jesh Ghangal, Rohini Garg, Mukesh Jain // PLoS One. — 2015. — Vol. 10, no. 3.
[12] Homeodomain-leucine zipper transcription factors: structural features of these proteins, unique to plants / Matias Capella, Pamela A Ribone, Agustin L Arce, Raquel L Chan // Plant transcription factors. — Elsevier, 2016. — Pp. 113-126.
[13] Gehring, Walter J. Homeo boxes in the study of development / Wal¬ter J Gehring // Science. — 1987. — Vol. 236, no. 4806. — Pp. 1245-1252.
[14] Ades, Sarah E. Specificity of minor-groove and major-groove interactions in a homeodomain-DNA complex / Sarah E Ades, Robert T Sauer // Biochemistry. — 1995. — Vol. 34, no. 44. — Pp. 14601-14608.
[15] Homeodomain-DNA recognition / Walter J Gehring, Yan Qiu Qian, Mar¬tin Billeter et al. // Cell. — 1994. — Vol. 78, no. 2. — Pp. 211-223.
[16] Wolberger, Cynthia. Homeodomain interactions / Cynthia Wolberger // Current opinion in structural biology. — 1996. — Vol. 6, no. 1. — Pp. 62¬68.
[17] Moens, Cecilia B. Hox cofactors in vertebrate development / Cecil¬ia B Moens, Licia Selleri // Developmental biology. — 2006. — Vol. 291, no. 2. — Pp. 193-206.
[18] The developmental gene Knotted-1 is a member of a maize homeobox gene family / Erik Vollbrecht, Bruce Veit, Neelima Sinha, Sarah Hake // Nature. — 1991. — Vol. 350, no. 6315. — Pp. 241-243.
[19] The true story of the HD-Zip family / Federico D Ariel, Pablo A Manavella, Carlos A Dezar, Raquel L Chan // Trends in plant science. — 2007. — Vol. 12, no. 9. — Pp. 419-426.
[20] Mukherjee, Krishanu. A comprehensive classification and evolutionary analysis of plant homeobox genes / Krishanu Mukherjee, Luciano Broc- chieri, Thomas R Biirglin // Molecular biology and evolution. — 2009. — Vol. 26, no. 12. — Pp. 2775-2794.
[21] Jain, Mukesh. Genome-wide identification, classification, evolutionary ex-pansion and expression analyses of homeobox genes in rice / Mukesh Jain, Akhilesh K Tyagi, Jitendra P Khurana // The FEBS journal. — 2008. — Vol. 275, no. 11. — Pp. 2845-2861.
[22] Genome-wide identification, expansion, and evolution analysis of home- obox genes and their expression profiles during root development in carrot / Feng Que, Guang-Long Wang, Tong Li et al. // Functional & integrative genomics. — 2018. — Vol. 18, no. 6. — Pp. 685-700.
[23] Genome-wide identification, classification, and expression pattern of home- obox gene family in Brassica rapa under various stresses / Nadeem Khan, Chun-mei Hu, Waleed Amjad Khan et al. // Scientific reports. — 2018. — Vol. 8, no. 1. — Pp. 1-17.
[24] Gene family structure, expression and functional analysis of HD-Zip III genes in angiosperm and gymnosperm forest trees / Caroline L Cote, Fran¬cis Boileau, Vicky Roy et al. // BMC Plant Biology. — 2010. — Vol. 10, no. 1. — P. 273.
[25] KNOTTEDl-like homeobox genes of a gymnosperm, Norway spruce, expressed during somatic embryogenesis / Helena I Hjortswang, Anni¬ka Sundas Larsson, Geeta Bharathan et al. // Plant Physiology and Bio¬chemistry. — 2002. — Vol. 40, no. 10. — Pp. 837-843.
[26] Ahuja, M Raj. Evolution of genome size in conifers / M Raj Ahuja, David B Neale // Silvae genetica. — 2005. — Vol. 54, no. 1-6. — Pp. 126¬137.
[27] Unique features of the loblolly pine (Pinus taeda L.) megagenome re¬vealed through sequence annotation / Jill L Wegrzyn, John D Liechty, Kristian A Stevens et al. // Genetics. — 2014. — Vol. 196, no. 3. — Pp. 891-909.
[28] Decoding the massive genome of loblolly pine using haploid DNA and novel assembly strategies / David B Neale, Jill L Wegrzyn, Kristian A Stevens et al. // Genome biology. — 2014. — Vol. 15, no. 3. — P. R59.
[29] Haas, BJ. De novo transcript sequence reconstruction from RNA-Seq: ref-erence generation and analysis with trinity. Nat Protocol 8 (8): 1494-1512. — 2013.
[30] Identification of differentially expressed genes in rna-seq data of arabidopsis thaliana: A compound distribution approach / Arfa Anjum, Seema Jaggi, Eldho Varghese et al. // Journal of Computational Biology. — 2016. — Vol. 23, no. 4. — Pp. 239-247.
[31] Robinson, Mark D. A scaling normalization method for differential ex¬pression analysis of RNA-seq data / Mark D Robinson, Alicia Oshlack // Genome biology. — 2010. — Vol. 11, no. 3. — P. R25.
[32] Wilkinson, Leland. The history of the cluster heat map / Leland Wilkinson, Michael Friendly // The American Statistician. — 2009. — Vol. 63, no. 2. — Pp. 179-184.
[33] Utilization of tissue ploidy level variation in de novo transcriptome assem¬bly of Pinus sylvestris / Dario I Ojeda, Tiina M Mattila, Tom Ruttink et al. // G3: Genes, Genomes, Genetics. — 2019. — Vol. 9, no. 10. — Pp. 3409-3421.
[34] Andrews, Simon. A quality control tool for high throughput sequence data. 2010 / Simon Andrews, A FastQC // Google Scholar. — 2015.
[35] Bolger, Anthony M. Trimmomatic: a flexible trimmer for Illumina sequence data / Anthony M Bolger, Marc Lohse, Bjoern Usadel // Bioinformatics. — 2014. — Vol. 30, no. 15. — Pp. 2114-2120.
[36] BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs / Felipe A Simao, Robert M Waterhouse, Panagi¬otis loannidis et al. // Bioinformatics. — 2015. — Vol. 31, no. 19. — Pp. 3210-3212.
[37] Finn, Robert D. HMMER web server: interactive sequence similarity searching / Robert D Finn, Jody Clements, Sean R Eddy // Nucleic acids research. — 2011. — Vol. 39, no. suppl_2. — Pp. W29-W37.
[38] Eddy, Sean. HMMER user’s guide. biological sequence analysis using pro¬file hidden Markov models / Sean Eddy. — 2003.
[39] The Pfam protein families database / Alex Bateman, Ewan Birney, Loren¬zo Cerruti et al. // Nucleic acids research. — 2002. — Vol. 30, no. 1. — Pp. 276-280.
[40] Robinson, Mark D. edgeR: a Bioconductor package for differential ex¬pression analysis of digital gene expression data / Mark D Robinson, Davis J McCarthy, Gordon K Smyth // Bioinformatics. — 2010. — Vol. 26, no. 1. — Pp. 139-140.
[41] Li, Bo. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome / Bo Li, Colin N Dewey // BMC bioinfor¬matics. — 2011. — Vol. 12, no. 1. — P. 323.
[42] A new transcriptome and transcriptome profiling of adult and larval tissue in the box jellyfish Alatina alata: an emerging model for studying venom, vision and sex / Cheryl Lewis Ames, Joseph F Ryan, Alexandra E Bely et al. // BMC genomics. — 2016. — Vol. 17, no. 1. — P. 650.
[43] Basic local alignment search tool / Stephen F Altschul, Warren Gish, Webb Miller et al. // Journal of molecular biology. — 1990. — Vol. 215, no. 3. — Pp. 403-410.
[44] Consortium, UniProt. Activities at the universal protein resource (UniProt) / UniProt Consortium // Nucleic acids research. — 2014. — Vol. 42, no. D1. — Pp. D191-D198.
[45] Du Plessis, Louis. The what, where, how and why of gene ontolo- gy—a primer for bioinformaticians / Louis Du Plessis, Nives Skunca, Christophe Dessimoz // Briefings in bioinformatics. — 2011. — Vol. 12, no. 6. — Pp. 723-735.
[46] Elhiti, Mohamed. Structure and function of homodomain-leucine zipper (HD-Zip) proteins / Mohamed Elhiti, Claudio Stasolla // Plant signaling & behavior. — 2009. — Vol. 4, no. 2. — Pp. 86-88.
[47] Arabidopsis HD-Zip II transcription factors control apical embryo devel-opment and meristem function / Luana Turchi, Monica Carabelli, Valenti¬no Ruzza et al. // Development. — 2013. — Vol. 140, no. 10. — Pp. 2118¬2129.
[48] Comprehensive analysis of WOX genes uncovers that WOX13 is in¬volved in phytohormone-mediated fiber development in cotton / Peng He, Yuzhou Zhang, Hao Liu et al. // BMC plant biology. — 2019. — Vol. 19, no. 1. — P. 312.
[49] Chew, William. Role of homeodomain leucine zipper (HD-Zip) IV tran-scription factors in plant development and plant protection from deleterious environmental factors / William Chew, Maria Hrmova, Sergiy Lopato // International journal of molecular sciences. — 2013. — Vol. 14, no. 4. — Pp. 8122-8147.
[50] Expression dynamics of WOX genes mark cell fate decisions during early embryonic patterning in Arabidopsis thaliana / Achim Haecker, Rita Grob- Hardt, Bernd Geiges et al. // Development. — 2004. — Vol. 131, no. 3. — Pp. 657-668.
[51] Analysis of the WUSCHEL-RELATED HOMEOBOX gene family in Pinus pinaster: New insights into the gene family evolution / Jose M Alvarez, Na¬talia Bueno, Rafael A Canas et al. // Plant Physiology and Biochemistry. — 2018. — Vol. 123. — Pp. 304-318.
[52] The WUSCHEL-RELATED HOMEOBOX 3 gene Pa WOX 3 regulates lateral organ formation in Norway spruce / Jose M Alvarez, Joel Sohlberg, Peter Engstrom et al. // New Phytologist. — 2015. — Vol. 208, no. 4. — Pp. 1078-1088.
[53] Conserved factors regulate signalling in Arabidopsis thaliana shoot and root stem cell organizers / Ananda K Sarkar, Marijn Luijten, Shun- suke Miyashima et al. // Nature. — 2007. — Vol. 446, no. 7137. — Pp. 811¬814.
[54] Molecular evolution and gene expression differences within the HD-Zip transcription factor family of Zea mays L. / Hude Mao, Lijuan Yu, Zhan- jie Li et al. // Genetica. — 2016. — Vol. 144, no. 2. — Pp. 243-257.
[55] The Arabidopsis Athb-2 and-4 genes are strongly induced by far-red-rich light / Monica Carabelli, Giovanna Sessa, Simona Baima et al. // The Plant Journal. — 1993. — Vol. 4, no. 3. — Pp. 469-479.
[56] Shade avoidance responses are mediated by the ATHB-2 HD-zip protein, a negative regulator of gene expression / Corinna Steindler, Antonella Mat- teucci, Giovanna Sessa et al. // Development. — 1999. — Vol. 126, no. 19. — Pp. 4235-4245.
[57] The HAT2 gene, a member of the HD-Zip gene family, isolated as an auxin inducible gene by DNA microarray screening, affects auxin response in Arabidopsis / Shinichiro Sawa, Maki Ohgishi, Hideki Goda et al. // The Plant Journal. — 2002. — Vol. 32, no. 6. — Pp. 1011-1022.
[58] A dynamic balance between gene activation and repression regulates the shade avoidance response in Arabidopsis / Giovanna Sessa, Monica Cara¬belli, Massimiliano Sassi et al. // Genes & development. — 2005. — Vol. 19, no. 23. — Pp. 2811-2815.
[59] Characterization of the class IV homeodomain-leucine zipper gene family in Arabidopsis / Miyuki Nakamura, Hiroshi Katsumata, Mitsutomo Abe et al. // Plant physiology. — 2006. — Vol. 141, no. 4. — Pp. 1363-1375.
[60] Transcription Factors BLH2 and BLH4 Regulate Demethylesterification of Homogalacturonan in Seed Mucilage / Yan Xu, Yiping Wang, Xi- aoyu Wang et al. // Plant Physiology. — 2020. — Vol. 183, no. 1. — Pp. 96-111.
[61] WUSCHEL-RELATED HOMEOBOX 4 (WOX 4)-like genes regulate cam- bial cell division activity and secondary growth in Populus trees / Melis Ku- cukoglu, Jeanette Nilsson, Bo Zheng et al. // New Phytologist. — 2017. — Vol. 215, no. 2. — Pp. 642-657.
[62] WOX4 promotes procambial development / Jiabing Ji, Josh Strable, Re¬na Shimizu et al. // Plant physiology. — 2010. — Vol. 152, no. 3. — Pp. 1346-1356.
[63] KNAT2/6b, a class I KNOX gene, impedes xylem differentiation by regu¬lating NAC domain transcription factors in poplar. / Yanqiu Zhao, Xue- qin Song, Houjun Zhou et al. // The New phytologist. — 2019.
[64] Hay, Angela. KNOX genes: versatile regulators of plant development and diversity / Angela Hay, Miltos Tsiantis // Development. — 2010. — Vol. 137, no. 19. — Pp. 3153-3165.
[65] Phyllotactic pattern and stem cell fate are determined by the Arabidop- sis homeobox gene BELLRINGER / Mary E Byrne, Andrew T Groover,
Joseph R Fontana, Robert A Martienssen // Development. — 2003. — Vol.
130, no. 17. — Pp. 3941-3950.