The Eucalyptus linker histone variant EgH1.3 cooperates with the transcription factor EgMYB1 to control lignin biosynthesis during wood formation

Genomic and transcriptomic studies on flavonoid biosynthesis in Lagerstroemia indica

BACKGROUND

Lagerstroemia indica is a widely cultivated ornamental woody shrub/tree of the family Lythraceae that is used as a traditional medicinal plant in East Asia and Egypt. However, unlike other ornamental woody plants, its genome is not well-investigated, which hindered the discovery of the key genes that regulate important traits and the synthesis of bioactive compounds.

RESULTS

In this study, the genomic sequences of L. indica were determined using several next-generation sequencing technologies. Altogether, 324.01 Mb sequences were assembled and 98.21% (318.21 Mb) of them were placed in 24 pseudo-chromosomes. The heterozygosity, repeated sequences, and GC residues occupied 1.65%, 29.17%, and 38.64% of the genome, respectively. In addition, 28,811 protein-coding gene models, 327 miRNAs, 552 tRNAs, 214 rRNAs, and 607 snRNAs were identified. The intra- and interspecies synteny and Ks analysis revealed that L. indica exhibits a hexaploidy. The co-expression profiles of the genes involved in the phenylpropanoid (PA) and flavonoid/anthocyanin (ABGs) pathways with the R2R3 MYB genes (137 members) showed that ten R2R3 MYB genes positively regulate flavonoid/anthocyanin biosynthesis. The colors of flowers with white, purple (PB), and deep purplish pink (DPB) petals were found to be determined by the levels of delphinidin-based (Dp) derivatives. However, the substrate specificities of LiDFR and LiOMT probably resulted in the different compositions of flavonoid/anthocyanin. In L. indica, two LiTTG1s (LiTTG1-1 and LiTTG1-2) were found to be the homologs of AtTTG1 (WD40). LiTTG1-1 was found to repress anthocyanin biosynthesis using the tobacco transient transfection assay.

CONCLUSIONS

This study showed that the ancestor L. indica experienced genome triplication approximately 38.5 million years ago and that LiTTG1-1 represses anthocyanin biosynthesis. Furthermore, several genes such as LiDFR, LiOMTs, and R2R3 LiMYBs are related to anthocyanin biosynthesis. Further studies are required to clarify the mechanisms and alleles responsible for flower color development.

Wood of trees: Cellular structure, molecular formation, and genetic engineering

Wood is an invaluable asset to human society due to its renewable nature, making it suitable for both sustainable energy production and material manufacturing. Additionally, wood derived from forest trees plays a crucial role in sequestering a significant portion of the carbon dioxide fixed during photosynthesis by terrestrial plants. Nevertheless, with the expansion of the global population and ongoing industrialization, forest coverage has been substantially decreased, resulting in significant challenges for wood production and supply. Wood production practices have changed away from natural forests toward plantation forests. Thus, understanding the underlying genetic mechanisms of wood formation is the foundation for developing high-quality, fast-growing plantation trees. Breeding ideal forest trees for wood production using genetic technologies has attracted the interest of many. Tremendous studies have been carried out in recent years on the molecular, genetic, and cell-biological mechanisms of wood formation, and considerable progress and findings have been achieved. These studies and findings indicate enormous possibilities and prospects for tree improvement. This review will outline and assess the cellular and molecular mechanisms of wood formation, as well as studies on genetically improving forest trees, and address future development prospects.

Transcriptome analysis of the transition from primary to secondary growth of vertical stem in Eucalyptus grandis

Eucalyptus was one of the most cultivated hardwood species worldwide, with rapid growth, good wood properties and a wide range of adaptability. Eucalyptus stem undergoes primary growth (longitudinal growth) followed by secondary growth (radial growth), which produces biomass that is an important source of energy worldwide. In order to better understand the genetic regulation of secondary growth in Eucalyptus grandis, Transcriptome analyses in stem segments along a developmental gradient from the third internode to the eleventh internode of E. grandis that spanned primary to secondary growth were carried out. 5,149 genes that were differentially expressed during stem development were identified. Combining the trend analysis by the Mfuzz method and the module-trait correlation analysis by the Weighted Gene Co-expression Network Analysis method, a total of 70 differentially expressed genes (DEGs) selected from 868 DEGs with high connectivity were found to be closely correlated with secondary growth. Results revealed that the differential expression of these DEGs suggests that they may involve in the primary growth or secondary growth. AP1, YAB2 TFs and EXP genes are highly expressed in the IN3, whereas NAC, MYB TFs are likely to be important for secondary growth. These results will expand our understanding of the complex molecular and cellular events of secondary growth and provide a foundation for future studies on wood formation in Eucalyptus.

Cytokinin promotes anthocyanin biosynthesis via regulating sugar accumulation and MYB113 expression in Eucalyptus

Anthocyanins are flavonoid-like substances that play important roles in plants' adaptation to various environmental stresses. In this research, we discovered that cytokinin (CK) alone could effectively induce the anthocyanin biosynthesis in Eucalyptus and many other perennial woody plant species, but not in tobacco and Arabidopsis, suggesting a diverse role of CK in regulating anthocyanin biosynthesis in different species. Transcriptomic and metabolomic strategies were used to further clarify the specific role of CK in regulating anthocyanin biosynthesis in Eucalyptus. The results showed that 801 and 2241 genes were differentially regulated at 6 and 24 h, respectively, after CK treatment. Pathway analysis showed that most of the differentially expressed genes were categorized into pathways related to cellular metabolism or transport of metabolites, including amino acids and sugars. The metabolomic results well supported the transcriptome data, which showed that most of the differentially regulated metabolites were related to the metabolism of sugar, amino acids and flavonoids. Moreover, CK treatment significantly induced the accumulation of sucrose in the CK-treated leaves, while sugar starvation mimicked by either defoliation or shading treatment of the basal leaves significantly reduced the sugar increase of the CK-treated leaves and thus inhibited CK-induced anthocyanin biosynthesis. The results of in vitro experiment also suggested that CK-induced anthocyanin in Eucalyptus was sugar-dependent. Furthermore, we identified an early CK-responsive transcription factor MYB113 in Eucalyptus, the expression of which was significantly upregulated by CK treatment in Eucalyptus, but was inhibited in Arabidopsis. Importantly, the overexpression of EgrMYB113 in the Eucalyptus hairy roots was associated with significant anthocyanin accumulation and upregulation of most of the anthocyanin biosynthetic genes. In conclusion, our study demonstrates a key role of CK in the regulation of anthocyanin biosynthesis in Eucalyptus, providing a molecular basis for further understanding the regulatory mechanism and diversity of hormone-regulated anthocyanin biosynthesis in different plant species.

Functional investigation of five R2R3-MYB transcription factors associated with wood development in Eucalyptus using DAP-seq-ML

A multi-tiered transcriptional network regulates xylem differentiation and secondary cell wall (SCW) formation in plants, with evidence of both conserved and lineage-specific SCW network architecture. We aimed to elucidate the roles of selected R2R3-MYB transcription factors (TFs) linked to Eucalyptus wood formation by identifying genome-wide TF binding sites and direct target genes through an improved DAP-seq protocol combined with machine learning for target gene assignment (DAP-seq-ML). We applied this to five TFs including a well-studied SCW master regulator (EgrMYB2; homolog of AtMYB83), a repressor of lignification (EgrMYB1; homolog of AtMYB4), a TF affecting SCW thickness and vessel density (EgrMYB137; homolog of PtrMYB074) and two TFs with unclear roles in SCW regulation (EgrMYB135 and EgrMYB122). Each DAP-seq TF peak set (average 12,613 peaks) was enriched for canonical R2R3-MYB binding motifs. To improve the reliability of target gene assignment to peaks, a random forest classifier was developed from Arabidopsis DAP-seq, RNA-seq, chromatin, and conserved noncoding sequence data which demonstrated significantly higher precision and recall to the baseline method of assigning genes to proximal peaks. EgrMYB1, EgrMYB2 and EgrMYB137 predicted targets showed clear enrichment for SCW-related biological processes. As validation, EgrMYB137 overexpression in transgenic Eucalyptus hairy roots increased xylem lignification, while its dominant repression in transgenic Arabidopsis and Populus reduced xylem lignification, stunted growth, and caused downregulation of SCW genes. EgrMYB137 targets overlapped significantly with those of EgrMYB2, suggesting partial functional redundancy. Our results show that DAP-seq-ML identified biologically relevant R2R3-MYB targets supported by the finding that EgrMYB137 promotes SCW lignification in planta.

Histone variant HTB4 delays leaf senescence by epigenetic control of Ib bHLH transcription factor-mediated iron homeostasis

Leaf senescence is an orderly process regulated by multiple internal factors and diverse environmental stresses including nutrient deficiency. Histone variants are involved in regulating plant growth and development. However, their functions and underlying regulatory mechanisms in leaf senescence remain largely unclear. Here, we found that H2B histone variant HTB4 functions as a negative regulator of leaf senescence. Loss of function of HTB4 led to early leaf senescence phenotypes that were rescued by functional complementation. RNA-seq analysis revealed that several Ib subgroup basic helix-loop-helix (bHLH) transcription factors (TFs) involved in iron (Fe) homeostasis, including bHLH038, bHLH039, bHLH100, and bHLH101, were suppressed in the htb4 mutant, thereby compromising the expressions of FERRIC REDUCTION OXIDASE 2 (FRO2) and IRON-REGULATED TRANSPORTER (IRT1), two important components of the Fe uptake machinery. Chromatin immunoprecipitation-quantitative polymerase chain reaction analysis revealed that HTB4 could bind to the promoter regions of Ib bHLH TFs and enhance their expression by promoting the enrichment of the active mark H3K4me3 near their transcriptional start sites. Moreover, overexpression of Ib bHLH TFs or IRT1 suppressed the premature senescence phenotype of the htb4 mutant. Our work established a signaling pathway, HTB4-bHLH TFs-FRO2/IRT1-Fe homeostasis, which regulates the onset and progression of leaf senescence.

Transcriptome-Based Construction of the Gibberellin Metabolism and Signaling Pathways in Eucalyptus grandis × E. urophylla, and Functional Characterization of GA20ox and GA2ox in Regulating Plant Development and Abiotic Stress Adaptations

Gibberellins (GAs) are the key regulators controlling plant growth, wood production and the stress responses in perennial woody plants. The role of GA in regulating the above-mentioned processes in Eucalyptus remain largely unclear. There is still a lack of systematic identification and functional characterization of GA-related genes in Eucalyptus. In this study, a total of 59,948 expressed genes were identified from the major vegetative tissues of the E. grandis × E. urophylla using transcriptome sequencing. Then, the key gene families in each step of GA biosynthesis, degradation and signaling were investigated and compared with those of Arabidopsis, rice, and Populus. The expression profile generated using Real-time quantitative PCR showed that most of these genes exhibited diverse expression patterns in different vegetative organs and in response to abiotic stresses. Furthermore, we selectively overexpressed EguGA20ox1, EguGA20ox2 and EguGA2ox1 in both Arabidopsis and Eucalyptus via Agrobacterium tumefaciens or A. rhizogenes-mediated transformation. Though both Arabidopsis EguGA20ox1- and EguGA20ox2-overexpressing (OE) lines exhibited better vegetative growth performance, they were more sensitive to abiotic stress, unlike EguGA2ox1-OE plants, which exhibited enhanced stress resistance. Moreover, overexpression of EguGA20ox in Eucalyptus roots caused significantly accelerated hairy root initiation and elongation and improved root xylem differentiation. Our study provided a comprehensive and systematic study of the genes of the GA metabolism and signaling and identified the role of GA20ox and GA2ox in regulating plant growth, stress tolerance, and xylem development in Eucalyptus; this could benefit molecular breeding for obtaining high-yield and stress-resistant Eucalyptus cultivars.

Transcriptome-wide analysis reveals core transcriptional regulators associated with culm development and variation in Dendrocalamus sinicus, the strongest woody bamboo in the world

Transcription factors (TFs) play indispensable roles in plant development and stress responses. As the largest woody bamboo species in the world, Dendrocalamus sinicus is endemic to Yunnan Province, China, and possesses two natural variants characterized by culm shape, namely straight or bent culms. Understanding the transcriptional regulation network of D. sinicus provides a unique opportunity to clarify the growth and development characteristics of woody bamboos. In this study, 10,236 TF transcripts belonging to 57 families were identified from transcriptome data of two variants at different developmental stages, from which we constructed a transcriptional regulatory network and unigene-coding protein-TFs interactive network of culm development for this attractive species. Gene function enrichment analysis revealed that hormone signaling and MAPK signaling pathways were two most enriched pathways in TF-regulated network. Based on PPI analysis, 50 genes interacting with nine TFs were screened as the core regulation components related to culm development. Of them, 18 synergistic genes of seven TFs, including nuclear cap-binding protein subunit 1, transcription factor GTE9-like, and ATP-dependent DNA helicase DDX11 isoform X1, involved in culm-shape variation. Most of these genes would interact with MYB, C3H, and ARF transcription factors. Six members with two each from ARF, C3H, and MYB transcription factor families and six key interacting genes (IAA3, IAA19, leucine-tRNA ligase, nuclear cap-binding protein subunit 1, elongation factor 2, and coiled-coil domain-containing protein 94) cooperate with these transcription factors were differentially expressed at development stage of young culms, and were validated by quantitative PCR. Our results represent a crucial step towards understanding the regulatory mechanisms of TFs involved in culm development and variation of D. sinicus.

Transcriptome profiling of celery petiole tissues reveals peculiarities of the collenchyma cell wall formation

Transcriptome and biochemical analyses are applied to individual plant cell types to reveal potential players involved in the molecular machinery of cell wall formation in specialized cells such as collenchyma. Plant collenchyma is a mechanical tissue characterized by an irregular, thickened cell wall and the ability to support cell elongation. The composition of the collenchyma cell wall resembles that of the primary cell wall and includes cellulose, xyloglucan, and pectin; lignin is absent. Thus, the processes associated with the formation of the primary cell wall in the collenchyma can be more pronounced compared to other tissues due to its thickening. Primary cell walls intrinsic to different tissues may differ in structure and composition, which should be reflected at the transcriptomic level. For the first time, we conducted transcriptome profiling of collenchyma strands isolated from young celery petioles and compared them with other tissues, such as parenchyma and vascular bundles. Genes encoding proteins involved in the primary cell wall formation during cell elongation, such as xyloglucan endotransglucosylase/hydrolases, expansins, and leucine-rich repeat proteins, were significantly activated in the collenchyma. As the key players in the transcriptome orchestra of collenchyma, xyloglucan endotransglucosylase/hydrolase transcripts were characterized in more detail, including phylogeny and expression patterns. The comprehensive approach that included transcriptome and biochemical analyses allowed us to reveal peculiarities of collenchyma cell wall formation and modification, matching the abundance of upregulated transcripts and their potential substrates for revealed gene products. As a result, specific isoforms of multigene families were determined for further functional investigation.

Potato Periderm Development and Tuber Skin Quality

The periderm is a corky tissue that replaces the epidermis when the latter is damaged, and is critical for preventing pathogen invasion and water loss. The periderm is formed through the meristematic activity of phellogen cells (cork cambium). The potato skin (phellem cells) composes the outer layers of the tuber periderm and is a model for studying cork development. Early in tuber development and following tuber expansion, the phellogen becomes active and produces the skin. New skin layers are continuously added by division of the phellogen cells until tuber maturation. Some physiological disorders of the potato tuber are related to abnormal development of the skin, including skinning injuries and russeting of smooth-skinned potatoes. Thus, characterizing the potato periderm contributes to modeling cork development in plants and helps to resolve critical agricultural problems. Here, we summarize the data available on potato periderm formation, highlighting tissue characteristics rather than the suberization processes.

Overexpression of EgrIAA20 from Eucalyptus grandis, a Non-Canonical Aux/IAA Gene, Specifically Decouples Lignification of the Different Cell-Types in Arabidopsis Secondary Xylem

Wood (secondary xylem) formation is regulated by auxin, which plays a pivotal role as an integrator of developmental and environmental cues. However, our current knowledge of auxin-signaling during wood formation is incomplete. Our previous genome-wide analysis of Aux/IAAs in Eucalyptus grandis showed the presence of the non-canonical paralog member EgrIAA20 that is preferentially expressed in cambium. We analyzed its cellular localization using a GFP fusion protein and its transcriptional activity using transactivation assays, and demonstrated its nuclear localization and strong auxin response repressor activity. In addition, we functionally tested the role of EgrIAA20 by constitutive overexpression in Arabidopsis to investigate for phenotypic changes in secondary xylem formation. Transgenic Arabidopsis plants overexpressing EgrIAA20 were smaller and displayed impaired development of secondary fibers, but not of other wood cell types. The inhibition in fiber development specifically affected their cell wall lignification. We performed yeast-two-hybrid assays to identify EgrIAA20 protein partners during wood formation in Eucalyptus, and identified EgrIAA9A, whose ortholog PtoIAA9 in poplar is also known to be involved in wood formation. Altogether, we showed that EgrIAA20 is an important auxin signaling component specifically involved in controlling the lignification of wood fibers.

Full-length transcriptomic identification of R2R3-MYB family genes related to secondary cell wall development in Cunninghamia lanceolata (Chinese fir)

BACKGROUND

R2R3-MYB is a class of transcription factor crucial in regulating secondary cell wall development during wood formation. The regulation of wood formation in gymnosperm has been understudied due to its large genome size. Using Single-Molecule Real-Time sequencing, we obtained full-length transcriptomic libraries from the developmental stem of Cunninghamia lanceolata, a perennial conifer known as Chinese fir. The R2R3-MYB of C. lanceolata (hereafter named as ClMYB) associated with secondary wall development were identified based on phylogenetic analysis, expression studies and functional study on transgenic line.

RESULTS

The evolutionary relationship of 52 ClMYBs with those from Arabidopsis thaliana, Eucalyptus grandis, Populus trichocarpa, Oryza sativa, two gymnosperm species, Pinus taeda, and Picea glauca were established by neighbour-joining phylogenetic analysis. A large number of ClMYBs resided in the woody-expanded subgroups that predominated with the members from woody dicots. In contrast, the woody-preferential subgroup strictly carrying the members of woody dicots contained only one candidate. The results suggest that the woody-expanded subgroup emerges before the gymnosperm/angiosperm split, while most of the woody-preferential subgroups are likely lineage-specific to woody dicots. Nine candidates shared the same subgroups with the A. thaliana orthologs, with known function in regulating secondary wall development. Gene expression analysis inferred that ClMYB1/2/3/4/5/26/27/49/51 might participate in secondary wall development, among which ClMYB1/2/5/26/27/49 were significantly upregulated in the highly lignified compression wood region, reinforcing their regulatory role associated with secondary wall development. ClMYB1 was experimentally proven a transcriptional activator that localised in the nucleus. The overexpression of ClMYB1 in Nicotiana benthamiana resulted in an increased lignin deposition in the stems. The members of subgroup S4, ClMYB3/4/5 shared the ERF-associated amphiphilic repression motif with AtMYB4, which is known to repress the metabolism of phenylpropanoid derived compounds. They also carried a core motif specific to gymnosperm lineage, suggesting divergence of the regulatory process compared to the angiosperms.

CONCLUSIONS

This work will enrich the collection of full-length gymnosperm-specific R2R3-MYBs related to stem development and contribute to understanding their evolutionary relationship with angiosperm species.

Histone methyltransferase ATX1 dynamically regulates fiber secondary cell wall biosynthesis in Arabidopsis inflorescence stem

Secondary wall thickening in the sclerenchyma cells is strictly controlled by a complex network of transcription factors in vascular plants. However, little is known about the epigenetic mechanism regulating secondary wall biosynthesis. In this study, we identified that ARABIDOPSIS HOMOLOG of TRITHORAX1 (ATX1), a H3K4-histone methyltransferase, mediates the regulation of fiber cell wall development in inflorescence stems of Arabidopsis thaliana. Genome-wide analysis revealed that the up-regulation of genes involved in secondary wall formation during stem development is largely coordinated by increasing level of H3K4 tri-methylation. Among all histone methyltransferases for H3K4me3 in Arabidopsis, ATX1 is markedly increased during the inflorescence stem development and loss-of-function mutant atx1 was impaired in secondary wall thickening in interfascicular fibers. Genetic analysis showed that ATX1 positively regulates secondary wall deposition through activating the expression of secondary wall NAC master switch genes, SECONDARY WALL-ASSOCIATED NAC DOMAIN PROTEIN1 (SND1) and NAC SECONDARY WALL THICKENING PROMOTING FACTOR1 (NST1). We further identified that ATX1 directly binds the loci of SND1 and NST1, and activates their expression by increasing H3K4me3 levels at these loci. Taken together, our results reveal that ATX1 plays a key role in the regulation of secondary wall biosynthesis in interfascicular fibers during inflorescence stem development of Arabidopsis.

A Novel NAC Transcription Factor From Eucalyptus, EgNAC141, Positively Regulates Lignin Biosynthesis and Increases Lignin Deposition

Wood formation is a complicated process under the control of a large set of transcription factors. NAC transcription factors are considered "master switches" in this process. However, few NAC members have been cloned and characterized in Eucalyptus, which is one of the most economically important woody plants. Here, we reported an NAC transcription factor from Eucalyptus grandis, EgNAC141, which has no Arabidopsis orthologs associated with xylogenesis-related processes. EgNAC141 was predominantly expressed in lignin-rich tissues, such as the stem and xylem. Overexpression of EgNAC141 in Arabidopsis resulted in stronger lignification, larger xylem, and higher lignin content. The expression of lignin biosynthetic genes in transgenic plants was significantly higher compared with wild-type plants. The transient expression of EgNAC141 activated the expression of Arabidopsis lignin biosynthetic genes in a dual-luciferase assay. Overall, these results showed that EgNAC141 is a positive regulator of lignin biosynthesis and may help us understand the regulatory mechanism of wood formation.

The chromatin-modifying protein HUB2 is involved in the regulation of lignin composition in xylem vessels

PIRIN2 (PRN2) was earlier reported to suppress syringyl (S)-type lignin accumulation of xylem vessels of Arabidopsis thaliana. In the present study, we report yeast two-hybrid results supporting the interaction of PRN2 with HISTONE MONOUBIQUITINATION2 (HUB2) in Arabidopsis. HUB2 has been previously implicated in several plant developmental processes, but not in lignification. Interaction between PRN2 and HUB2 was verified by β-galactosidase enzymatic and co-immunoprecipitation assays. HUB2 promoted the deposition of S-type lignin in the secondary cell walls of both stem and hypocotyl tissues, as analysed by pyrolysis-GC/MS. Chemical fingerprinting of individual xylem vessel cell walls by Raman and Fourier transform infrared microspectroscopy supported the function of HUB2 in lignin deposition. These results, together with a genetic analysis of the hub2 prn2 double mutant, support the antagonistic function of PRN2 and HUB2 in deposition of S-type lignin. Transcriptome analyses indicated the opposite regulation of the S-type lignin biosynthetic gene FERULATE-5-HYDROXYLASE1 by PRN2 and HUB2 as the underlying mechanism. PRN2 and HUB2 promoter activities co-localized in cells neighbouring the xylem vessel elements, suggesting that the S-type lignin-promoting function of HUB2 is antagonized by PRN2 for the benefit of the guaiacyl (G)-type lignin enrichment of the neighbouring xylem vessel elements.

Integrative analysis of wood biomass and developing xylem transcriptome provide insights into mechanisms of lignin biosynthesis in wood formation of Pinus massoniana

Lignin is an important renewable energy source as an excellent new battery fuel and ideal substitutes for the petrochemical industry. However, the molecular mechanism underlying lignin biosynthesis in wood formation of P. massoniana remains unexplored. Thus, an integrative analysis of wood biomass and the developing xylem transcriptome was performed to identify genes involved in lignin biosynthesis. A total of 1624 differentially expressed genes (DEGs) were identified, consisting of 797 upregulated and 827 downregulated genes (MaxG vs MinG). Additionally, 122 candidate genes and 17 DEGs were successfully annotated to the lignin biosynthesis pathway. All upregulated MYB and NAC genes were regulators of secondary cell wall formation. Moreover, the qRT-PCR analyses shown that 9 lignin biosynthesis-related genes and 7 transcription factor-encoding genes were upregulated (MaxG vs MinG), which indicated that the downregulation of lignin biosynthesis-related genes might be the possible causes of growth retardation and dwarf phenotype in some P. massoniana individuals. The identification of lignin biosynthesis-related genes can provide valuable genetic basis and resource for further researches on molecular mechanisms of lignin biosynthesis and contribute to the future investigations of bioengineering and synthetic biology to regulate lignin content in wood formation for the pulp and wood utilization industry.

Grapevine VlbZIP30 improves drought resistance by directly activating VvNAC17 and promoting lignin biosynthesis through the regulation of three peroxidase genes

Drought stress severely affects grapevine quality and yield, and recent reports have revealed that lignin plays an important role in protection from drought stress. Since little is known about lignin-mediated drought resistance in grapevine, we investigated its significance. Herein, we show that VlbZIP30 mediates drought resistance by activating the expression of lignin biosynthetic genes and increasing lignin deposition. Transgenic grapevine plants overexpressing VlbZIP30 exhibited lignin deposition (mainly G and S monomers) in the stem secondary xylem under control conditions, which resulted from the upregulated expression of VvPRX4 and VvPRX72. Overexpression of VlbZIP30 improves drought tolerance, characterized by a reduction in the water loss rate, maintenance of an effective photosynthesis rate, and increased lignin content (mainly G monomer) in leaves under drought conditions. Electrophoretic mobility shift assay, luciferase reporter assays, and chromatin immunoprecipitation-qPCR assays indicated that VlbZIP30 directly binds to the G-box cis-element in the promoters of lignin biosynthetic (VvPRX N1) and drought-responsive (VvNAC17) genes to regulate their expression. In summary, we report a novel VlbZIP30-mediated mechanism linking lignification and drought tolerance in grapevine. The results of this study may be of value for the development of molecular breeding strategies to produce drought-resistant fruit crops.

DNA methylation and its effects on gene expression during primary to secondary growth in poplar stems

Background

As an important epigenetic mark, 5-methylcytosine (5mC) methylation is involved in many DNA-dependent biological processes and plays a role during development and differentiation of multicellular organisms. However, there is still a lack of knowledge about the dynamic aspects and the roles of global 5mC methylation in wood formation in tree trunks. In this study, we not only scrutinized single-base resolution methylomes of primary stems (PS), transitional stems (TS), and secondary stems (SS) of Populus trichocarpa using a high-throughput bisulfite sequencing technique, but also analyzed the effects of 5mC methylation on the expression of genes involved in wood formation.

Results

The overall average percentages of CG, CHG, and CHH methylation in poplar stems were ~ 53.6%, ~ 37.7%, and ~ 8.5%, respectively, and the differences of 5mC in genome-wide CG/CHG/CHH contexts among PS, TS, and SS were statistically significant (p < 0.05). The evident differences in CG, CHG, and CHH methylation contexts among 2 kb proximal promoters, gene bodies, and 2 kb downstream regions were observed among PS, TS, and SS. Further analysis revealed a perceptible global correlation between 5mC methylation levels of gene bodies and transcript levels but failed to reveal a correlation between 5mC methylation levels of proximal promoter regions and transcript levels. We identified 653 and 858 DMGs and 4978 and 4780 DEGs in PS vs TS and TS vs SS comparisons, respectively. Only 113 genes of 653 DMGs and 4978 DEGs, and 114 genes of 858 DMGs and 4780 DEG were common. Counterparts of some of these common genes in other species, including Arabidopsis thaliana, are known to be involved in secondary cell wall biosynthesis and hormone signaling. This indicates that methylation may directly modulate wood formation genes and indirectly attune hormone signaling genes, which in turn impact wood formation.

Conclusions

DNA methylation only marginally affects pathway genes or regulators involved in wood formation, suggesting that further studies of wood formation should lean towards the indirect effects of methylation. The information and data we provide here will be instrumental for understanding the roles of methylation in wood formation in tree species.

Transcriptomic Profiling of Cryptomeria fortunei Hooibrenk Vascular Cambium Identifies Candidate Genes Involved in Phenylpropanoid Metabolism

Cryptomeria fortunei Hooibrenk (Chinese cedar) is a coniferous tree from southern China that has an important function in landscaping and timber production. Lignin is one of the key components of secondary cell walls, which have a crucial role in conducting water and providing mechanical support for the upward growth of plants. It is mainly biosynthesized via the phenylpropanoid metabolic pathway, of which the molecular mechanism remains so far unresolved in C. fortunei. In order to obtain further insight into this pathway, we performed transcriptome sequencing of the C. fortunei cambial zone at 5 successive growth stages. We generated 78,673 unigenes from transcriptome data, of which 45,214 (57.47%) were successfully annotated in the non-redundant protein database (NR). A total of 8975 unigenes were identified to be significantly differentially expressed between Sample_B and Sample_A after analyzing their expression profiles. Of the differentially expressed genes (DEGs), 6817 (75.96%) and 2158 (24.04%) were up- and down-regulated, respectively. 83 DEGs were involved in phenylpropanoid metabolism, 37 DEGs that encoded v-Myb avian myeloblastosis viral oncogene homolog (MYB) transcription factor (TF), and many candidates that encoded lignin synthesizing enzymes. These findings contribute to understanding the expression pattern of C. fortunei cambial zone transcriptome. Furthermore, our results provide additional insight towards understanding the molecular mechanisms of wood formation in C. fortunei.

Transcriptional analysis of arogenate dehydratase genes identifies a link between phenylalanine biosynthesis and lignin biosynthesis

Biogenesis of the secondary cell wall in trees involves the massive biosynthesis of the phenylalanine-derived polymer lignin. Arogenate dehydratase (ADT) catalyzes the last, and rate-limiting, step of the main pathway for phenylalanine biosynthesis. In this study, we found that transcript levels for several members of the large ADT gene family, including ADT-A and ADT-D, were enhanced in compression wood of maritime pine, a xylem tissue enriched in lignin. Transcriptomic analysis of maritime pine silenced for PpMYB8 revealed that this gene plays a critical role in coordinating the deposition of lignin with the biosynthesis of phenylalanine. Specifically, it was found that ADT-A and ADT-D were strongly down-regulated in PpMYB8-silenced plants and that they were transcriptionally regulated through direct interaction of this transcription factor with regulatory elements present in their promoters. Another transcription factor, PpHY5, exhibited an expression profile opposite to that of PpMYB8 and also interacted with specific regulatory elements of ADT-A and ADT-D genes, suggesting that it is involved in transcriptional regulation of phenylalanine biosynthesis. Taken together, our results reveal that PpMYB8 and PpHY5 are involved in the control of phenylalanine formation and its metabolic channeling for lignin biosynthesis and deposition during wood formation in maritime pine.

A Molecular Blueprint of Lignin Repression

Although lignin is essential to ensure the correct growth and development of land plants, it may be an obstacle to the production of lignocellulosics-based biofuels, and reduces the nutritional quality of crops used for human consumption or livestock feed. The need to tailor the lignocellulosic biomass for more efficient biofuel production or for improved plant digestibility has fostered considerable advances in our understanding of the lignin biosynthetic pathway and its regulation. Most of the described regulators are transcriptional activators of lignin biosynthesis, but considerably less attention has been devoted to the repressors of this pathway. We provide a comprehensive overview of the molecular factors that negatively impact on the lignification process at both the transcriptional and post-transcriptional levels.

Regulation of Cell Wall Thickening by a Medley of Mechanisms

To provide physical support for developing structures and to withstand the pressures associated with water and nutrient transport, some cells deposit a secondary cell wall, a rigid matrix of polysaccharide and phenolic biopolymers. The biosynthesis and deposition of these materials and the patterning of secondary wall-forming cells is controlled by a network of transcription factors. However, recent work suggests that this network forms the core of a more complex, multilevel regulatory system. This expanded system includes epigenetic, post-transcriptional, and post-translational regulation, and is coordinated with other pathways controlling primary growth and responses to environmental stimuli. New findings expand the set of transcription factors identified as secondary cell wall regulators and reveal novel regulatory processes that further govern secondary wall biogenesis.

GhbHLH18 negatively regulates fiber strength and length by enhancing lignin biosynthesis in cotton fibers

Cotton fibers are developed epidermal cells of the seed coat and contain large amounts of cellulose and minor lignin-like components. Lignin in the cell walls of cotton fibers effectively provides mechanical strength and is also presumed to restrict fiber elongation and secondary cell wall synthesis. To analyze the effect of lignin and lignin-like phenolics on fiber quality and the transcriptional regulation of lignin synthesis in cotton fibers, we characterized the function of a bHLH transcription factor, GhbHLH18, during fiber elongation stage. GhbHLH18 knock-down plants have longer and stronger fibers, and accumulate less lignin-like phenolics in mature cotton fibers than control plants. By mining public transcriptomic data for developing fibers, we discovered that GhbHLH18 is coexpressed with most lignin synthesis pathway genes. Furthermore, we showed that GhbHLH18 strongly binds to the E-box in the promoter region of GhPER8 and activates its expression. Transient over expression of GhPER8 protein in tobacco leaves significantly decreased the content of coniferyl alcohol and sinapic alcohol-the substrate respectively for G-lignin and S-lignin biosynthesis. These results suggest that GhbHLH18 is negatively associated with fiber quality by activating peroxidase-mediated lignin metabolism, thus the paper represents an alternative strategy to improve fiber quality.

A Standardized Synthetic Eucalyptus Transcription Factor and Promoter Panel for Re-engineering Secondary Cell Wall Regulation in Biomass and Bioenergy Crops

Re-engineering of transcriptional networks regulating secondary cell wall formation may allow the improvement of plant biomass in widely grown plantation crops such as Eucalyptus. However, there is currently a scarcity of freely available standardized biological parts (e.g., Phytobricks) compatible with Type IIS assembly approaches from forest trees, and there is a need to accelerate transcriptional network inference in nonmodel biomass crops. Here we describe the design and synthesis of a versatile three-panel biological parts collection of 221 secondary cell wall-related Eucalyptus grandis transcription factor coding sequences and 65 promoters that are compatible with GATEWAY, Golden Gate, MoClo, and GoldenBraid DNA assembly methods and generally conform to accepted Phytobrick syntaxes. This freely available resource is intended to accelerate synthetic biology applications in multiple plant biomass crops and enable reconstruction of secondary cell wall transcriptional networks using high-throughput assays such as DNA affinity purification sequencing (DAP-seq) and enhanced yeast one-hybrid (eY1H) screening.

RNA-Seq analysis of differential gene expression in Betula luminifera xylem during the early stages of tension wood formation

Background

Betula luminifera H. Winkler, which is widely distributed in southern China, is an economically important broadleaf tree species. However, little genomic information of B. luminifera is available, and little is known about the molecular mechanisms of wood formation in this species. Meanwhile, few efforts have focused on investigating the early transcriptional changes during tension wood formation in woody plants.

Results

A reference transcriptome dataset was first generated containing 45,700 Unigenes, and 35,135 (76.9%) Unigenes were annotated by a BLAST similarity search against four public databases. Then, based on an anatomical investigation, the global gene expression changes during the early stages of tension wood formation were analyzed. Gene expression profiling showed that a total of 13,273 Unigenes were differentially regulated during the early stages of tension wood formation. Most genes involved in cellulose and lignin biosynthesis were highlighted to reveal their biological importance in tension wood formation. In addition, the transcription levels of many genes involved in the auxin response pathway were significantly changed during the early stages of tension wood formation. Furthermore, 18 TFs co-expressed with key enzymes of cellulose synthesis were identified.

Conclusions

Our results revealed the transcriptional changes associated with TW formation and identified potential key genes in the regulation of this process. These results will help to dissect the molecular mechanism of wood formation and provide key candidate genes for marker-assisted selection in B. luminifera.

Temporal analysis of Arabidopsis genes activated by Eucalyptus grandis NAC transcription factors associated with xylem fibre and vessel development

Secondary cell wall (SCW) deposition in Arabidopsis is regulated among others by NAC transcription factors, where SND1 chiefly initiates xylem fibre differentiation while VND6 controls metaxylem vessel SCW development, especially programmed cell death and wall patterning. The translational relevance of Arabidopsis SCW regulation theory and the utility of characterized transcription factors as modular synthetic biology tools for improving commercial fibre crops is unclear. We investigated inter-lineage gene activation dynamics for potential fibre and vessel differentiation regulators from the widely grown hardwood Eucalyptus grandis (Myrtales). EgrNAC26, a VND6 homolog, and EgrNAC61, an SND1 homolog, were transiently expressed in Arabidopsis mesophyll protoplasts in parallel to determine early and late (i.e. 7 and 14 hours post-transfection) gene targets. Surprisingly, across the time series EgrNAC26 activated only a subset of SCW-related transcription factors and biosynthetic genes activated by EgrNAC61, specializing instead in targeting vessel-specific wall pit and programmed cell death markers. Promoters of EgrNAC26 and EgrNAC61 both induced reporter gene expression in vessels of young Arabidopsis plants, with EgrNAC61 also conferring xylem- and cork cambium-preferential expression in Populus. Our results demonstrate partial conservation, with notable exceptions, of SND1 and VND6 homologs in Eucalyptus and a first report of cork cambium expression for EgrNAC61.

The R2R3MYB Gene Family in Phyllostachys edulis: Genome-Wide Analysis and Identification of Stress or Development-Related R2R3MYBs

The MYB transcription factor (TF) is one of the largest gene families in plants and involved to multiple biological processes. However, little is known about the MYB family and its functional role in the genome of moso bamboo. In the present study, a total of 114 R2R3MYB genes were first identified from moso bamboo genome and full-length non-chimeric (FLNC) reads. Phylogenetic analysis coupled with gene structure analysis and motif determination resulted in the division of these PheR2R3MYBs into 17 subgroups. The position of eight proteins along an external branch in the phylogenetic tree suggested their relatively ancient origin. The genes in this group were all substituted by (Met, M)/(Arg, R) at conservative W residues in both R2 and R3 repeats, and half were found to possess no transcriptional activation activity. The analysis of evolutionary patterns and divergence suggests that the expansion of PheMYBs was mainly attributable to whole genome duplication (WGD) under different selection pressures. Expressional analysis based on microarray and qRT-PCR data performed diverse expression patterns of R2R3MYBs in response to both various abiotic stimuli and flower development. Furthermore, the co-expression analysis of R2R3MYBs suggested an intricate interplay of growth- and stress-related responses. Finally, we found a hub gene, PheMYB4, was involved in a complex proteins interaction network. Further functional analysis indicated that ectopic overexpression of its homologous gene, PheMYB4-1, could increase tolerance to cold treatment and sensitivity to drought and salt treatment of transgenic Arabidopsis seedlings. These findings provide comprehensive insights into the MYB family members in moso bamboo and offer candidate MYB genes for further studies on their roles in stress resistance.

Co-expression network of transcription factors reveal ethylene-responsive element-binding factor as key regulator of wood phenotype in Eucalyptus tereticornis

Suitability of wood biomass for pulp production is dependent on the cellular architecture and composition of secondary cell wall. Presently, systems genetics approach is being employed to understand the molecular basis of trait variation and co-expression network analysis has enabled holistic understanding of complex trait such as secondary development. Transcription factors (TFs) are reported as key regulators of meristematic growth and wood formation. The hierarchical TF network is a multi-layered system which interacts with downstream structural genes involved in biosynthesis of cellulose, hemicelluloses and lignin. Several TFs have been associated with wood formation in tree species such as Populus, Eucalyptus, Picea and Pinus. However, TF-specific co-expression networks to understand the interaction between these regulators are not reported. In the present study, co-expression network was developed for TFs expressed during wood formation in Eucalyptus tereticornis and ethylene-responsive element-binding factor, EtERF2, was identified as the major hub transcript which co-expressed with other secondary cell wall biogenesis-specific TFs such as EtSND2, EtVND1, EtVND4, EtVND6, EtMYB70, EtGRAS and EtSCL8. This study reveals a probable role of ethylene in determining natural variation in wood properties in Eucalyptus species. Understanding this transcriptional regulation underpinning the complex bio-processing trait of wood biomass will complement the Eucalyptus breeding program through selection of industrially suitable phenotypes by marker-assisted selection.

A comparative transcriptomic approach to understanding the formation of cork

The transcriptome comparison of two oak species reveals possible candidates accounting for the exceptionally thick and pure cork oak phellem, such as those involved in secondary metabolism and phellogen activity. Cork oak, Quercus suber, differs from other Mediterranean oaks such as holm oak (Quercus ilex) by the thickness and organization of the external bark. While holm oak outer bark contains sequential periderms interspersed with dead secondary phloem (rhytidome), the cork oak outer bark only contains thick layers of phellem (cork rings) that accumulate until reaching a thickness that allows industrial uses. Here we compare the cork oak outer bark transcriptome with that of holm oak. Both transcriptomes present similitudes in their complexity, but whereas cork oak external bark is enriched with upregulated genes related to suberin, which is the main polymer responsible for the protective function of periderm, the upregulated categories of holm oak are enriched in abiotic stress and chromatin assembly. Concomitantly with the upregulation of suberin-related genes, there is also induction of regulatory and meristematic genes, whose predicted activities agree with the increased number of phellem layers found in the cork oak sample. Further transcript profiling among different cork oak tissues and conditions suggests that cork and wood share many regulatory mechanisms, probably reflecting similar ontogeny. Moreover, the analysis of transcripts accumulation during the cork growth season showed that most regulatory genes are upregulated early in the season when the cork cambium becomes active. Altogether our work provides the first transcriptome comparison between cork oak and holm oak outer bark, which unveils new regulatory candidate genes of phellem development.

Integrated analysis and transcript abundance modelling of H3K4me3 and H3K27me3 in developing secondary xylem

Despite the considerable contribution of xylem development (xylogenesis) to plant biomass accumulation, its epigenetic regulation is poorly understood. Furthermore, the relative contributions of histone modifications to transcriptional regulation is not well studied in plants. We investigated the biological relevance of H3K4me3 and H3K27me3 in secondary xylem development using ChIP-seq and their association with transcript levels among other histone modifications in woody and herbaceous models. In developing secondary xylem of the woody model Eucalyptus grandis, H3K4me3 and H3K27me3 genomic spans were distinctly associated with xylogenesis-related processes, with (late) lignification pathways enriched for putative bivalent domains, but not early secondary cell wall polysaccharide deposition. H3K27me3-occupied genes, of which 753 (~31%) are novel targets, were enriched for transcriptional regulation and flower development and had significant preferential expression in roots. Linear regression models of the ChIP-seq profiles predicted ~50% of transcript abundance measured with strand-specific RNA-seq, confirmed in a parallel analysis in Arabidopsis where integration of seven additional histone modifications each contributed smaller proportions of unique information to the predictive models. This study uncovers the biological importance of histone modification antagonism and genomic span in xylogenesis and quantifies for the first time the relative correlations of histone modifications with transcript abundance in plants.

Mutant Transcriptome Sequencing Provides Insights into Pod Development in Peanut (Arachis hypogaea L.)

Pod size is the major yield component and a key target trait that is selected for in peanut breeding. However, although numerous quantitative trait loci (QTLs) for peanut pod size have been described, the molecular mechanisms underlying the development of this characteristic remain elusive. A peanut mutant with a narrower pod was developed in this study using ethyl methanesulfonate (EMS) mutagenesis and designated as the "pod width" mutant line (pw). The fresh pod weight of pw was only about 40% of that seen in the wild-type (WT) Zhonghua16, while the hull and seed filling of the mutant both also developed at earlier stages. Pods from both pw and WT lines were sampled 20, 40, and 60 days after flowering (DAF) and used for RNA-Seq analysis; the results revealed highly differentially expressed lignin metabolic pathway genes at all three stages, but especially at DAF 20 and DAF 40. At the same time, expression of genes related to auxin signal transduction was found to be significantly repressed during the pw early pod developmental stage. A genome-wide comparative analysis of expression profiles revealed 260 differentially expressed genes (DEGs) across all three stages, and two candidate genes, c26901_g1 (CAD) and c37339_g1 (ACS), responsible for pod width were identified by integrating expression patterns and function annotation of the common DEGs within the three stages. Taken together, the information provided in this study illuminates the processes underlying peanut pod development, and will facilitate further identification of causal genes and the development of improved peanut varieties with higher yields.

The Woody-Preferential Gene EgMYB88 Regulates the Biosynthesis of Phenylpropanoid-Derived Compounds in Wood

Comparative phylogenetic analyses of the R2R3-MYB transcription factor family revealed that five subgroups were preferentially found in woody species and were totally absent from Brassicaceae and monocots (Soler et al., 2015). Here, we analyzed one of these subgroups (WPS-I) for which no gene had been yet characterized. Most Eucalyptus members of WPS-I are preferentially expressed in the vascular cambium, the secondary meristem responsible for tree radial growth. We focused on EgMYB88, which is the most specifically and highly expressed in vascular tissues, and showed that it behaves as a transcriptional activator in yeast. Then, we functionally characterized EgMYB88 in both transgenic Arabidopsis and poplar plants overexpressing either the native or the dominant repression form (fused to the Ethylene-responsive element binding factor-associated Amphiphilic Repression motif, EAR). The transgenic Arabidopsis lines had no phenotype whereas the poplar lines overexpressing EgMYB88 exhibited a substantial increase in the levels of the flavonoid catechin and of some salicinoid phenolic glycosides (salicortin, salireposide, and tremulacin), in agreement with the increase of the transcript levels of landmark biosynthetic genes. A change in the lignin structure (increase in the syringyl vs. guaiacyl, S/G ratio) was also observed. Poplar lines overexpressing the EgMYB88 dominant repression form did not show a strict opposite phenotype. The level of catechin was reduced, but the levels of the salicinoid phenolic glycosides and the S/G ratio remained unchanged. In addition, they showed a reduction in soluble oligolignols containing sinapyl p-hydroxybenzoate accompanied by a mild reduction of the insoluble lignin content. Altogether, these results suggest that EgMYB88, and more largely members of the WPS-I group, could control in cambium and in the first layers of differentiating xylem the biosynthesis of some phenylpropanoid-derived secondary metabolites including lignin.