Structural, evolutionary and functional analysis of the NAC domain protein family in Eucalyptus

Genome-wide analysis of NAC transcription factors and exploration of candidate genes regulating selenium metabolism in Broussonetia papyrifera

MAIN CONCLUSION

Genome-wide identification revealed 79 BpNAC genes belonging to 16 subfamilies, and their gene structures and evolutionary relationships were characterized. Expression analysis highlighted their importance in plant selenium stress responses. Paper mulberry (Broussonetia papyrifera), a deciduous arboreal plant of the Moraceae family, is distinguished by its leaves, which are abundant in proteins, polysaccharides, and flavonoids, positioning it as a novel feedstock. NAC transcription factors, exclusive to plant species, are crucial in regulating growth, development, and response to biotic and abiotic stress. However, extensive characterization of the NAC family within paper mulberry is lacking. In this study, 79 BpNAC genes were identified from the paper mulberry genome, with an uneven distribution across 13 chromosomes. A comprehensive, genome-wide analysis of BpNACs was performed, including investigating gene structures, promoter regions, and chromosomal locations. Phylogenetic tree analysis, alongside comparisons with Arabidopsis thaliana NACs, allowed for categorizing these genes into 16 subfamilies in alignment with gene structure and motif conservation. Collinearity analysis suggested a significant homologous relationship between the NAC genes of paper mulberry and those in Morus notabilis, Ficus hispida, Antiaris toxicaria, and Cannabis sativa. Integrating transcriptome data and Se content revealed that 12 BpNAC genes were associated with selenium biosynthesis. Subsequent RT-qPCR analysis corroborated the correlation between BpNAC59, BpNAC62 with sodium selenate, and BpNAC55 with sodium selenite. Subcellular localization experiments revealed the nuclear functions of BpNAC59 and BpNAC62. This study highlights the potential BpNAC transcription factors involved in selenium metabolism, providing a foundation for strategically breeding selenium-fortified paper mulberry.

Gene co-expression network analysis in areca floral organ and the potential role of the AcMADS17 and AcMADS23 in transgenic Arabidopsis

Areca catechu L., a monocot belonging to the palm family, is monoecious, with female and male flowers separately distributed on the same inflorescence. To discover the molecular mechanism of flower development in Areca, we sequenced different floral samples to generate tissue-specific transcriptomic profiles. We conducted a comparative analysis of the transcriptomic profiles of apical sections of the inflorescence with male flowers and the basal section of the inflorescence with female flowers. Based on the RNA sequencing dataset, we applied weighted gene co-expression network analysis (WGCNA) to identify sepal, petal, stamen, stigma and other specific modules as well as hub genes involved in specific floral organ development. The syntenic and expression patterns of AcMADS-box genes were analyzed in detail. Furthermore, we analyzed the open chromatin regions and transcription factor PI binding sites in male and female flowers by assay for transposase-accessible chromatin sequencing (ATAC-seq) assay. Heterologous expression revealed the important role of AcMADS17 and AcMADS23 in floral organ development. Our results provide a valuable genomic resource for the functional analysis of floral organ development in Areca.

Phosphorylation of birch BpNAC90 improves the activation of gene expression to confer drought tolerance

The NAC transcription factors (TFs) play important roles in mediating abiotic stress tolerance; however, the mechanism is still not fully known. Here, an NAC gene (BpNAC90) from a gene regulatory network of Betula platyphylla (birch) that responded to drought was characterized. Overexpression and knockout of BpNAC90 displayed increased and reduced drought tolerance, respectively, relative to wild-type (WT) birch. BpNAC90 binds to different DNA motifs to regulate target genes in conferring drought tolerance, such as Eomes2, ABRE and Tgif2. BpNAC90 is phosphorylated by drought stress at Ser 205 by birch SNF1-related protein kinase 2 (BpSRK2A). Mutated BpNAC90 (termed S205A) with abolished phosphorylation, was transformed into birch for overexpression. The transgenic S205A plants displayed significantly reduced drought tolerance compared with plants overexpressing BpNAC90, but still showed increased drought tolerance relative to WT birch. At the same time, S205A showed a decreased capability to bind to motifs and reduced activation of target gene expression, which contributed to the reduced drought tolerance. Additionally, BpSRK2A and BpNAC90 can be induced by drought stress and form a complex to phosphorylate BpNAC90. The results together indicated that phosphorylation of BpNAC90 is necessary in conferring drought tolerance in birch.

Comprehensive Analysis of NAC Transcription Factors Reveals Their Evolution in Malvales and Functional Characterization of AsNAC019 and AsNAC098 in Aquilaria sinensis

NAC is a class of plant-specific transcription factors that are widely involved in the growth, development and (a)biotic stress response of plants. However, their molecular evolution has not been extensively studied in Malvales, especially in Aquilaria sinensis, a commercial and horticultural crop that produces an aromatic resin named agarwood. In this study, 1502 members of the NAC gene family were identified from the genomes of nine species from Malvales and three model plants. The macroevolutionary analysis revealed that whole genome duplication (WGD) and dispersed duplication (DSD) have shaped the current architectural structure of NAC gene families in Malvales plants. Then, 111 NAC genes were systemically characterized in A. sinensis. The phylogenetic analysis suggests that NAC genes in A. sinensis can be classified into 16 known clusters and four new subfamilies, with each subfamily presenting similar gene structures and conserved motifs. RNA-seq analysis showed that AsNACs presents a broad transcriptional response to the agarwood inducer. The expression patterns of 15 AsNACs in A. sinensis after injury treatment indicated that AsNAC019 and AsNAC098 were positively correlated with the expression patterns of four polyketide synthase (PKS) genes. Additionally, AsNAC019 and AsNAC098 were also found to bind with the AsPKS07 promoter and activate its transcription. This comprehensive analysis provides valuable insights into the molecular evolution of the NAC gene family in Malvales plants and highlights the potential mechanisms of AsNACs for regulating secondary metabolite biosynthesis in A. sinensis, especially for the biosynthesis of 2-(2-phenyl) chromones in agarwood.

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.

Genome-Wide Identification, Evolution and Expression Profile Analysis of NAC Transcription Factor in Simmondsia chinensis

NAC transcription factors (TFs) are one of the largest plant-specific gene families and play important roles in plant growth, development, and the biotic and abiotic stress response. Although the sequencing of Jojoba (Simmondsia chinensis) has been completed, the genome-wide identification and analysis of its NAC TFs has not been reported. In this study, a total of 57 genes were identified in Jojoba, which were divided into eight groups based on phylogenetic analysis. The genes clustered in the same groups have a similar gene structure and motif distribution. Based on the analysis of cis-elements in NAC TFs, nine cis-acting elements were identified in the promoter region that involved in light response, hormonal response, and stress response. Synteny analysis showed a greater collinearity between Jojoba and V. vinifera than Arabidopsis thaliana. The 24 genes in the Jojoba NAC TFs are derived from fragment replication, which may be the main source of NAC amplification. Gene expression analysis identified seven genes that were highly expressed in seeds. The differential expression analysis of NAC TFs in cotyledon and embryonic axis tissues showed that the expression of 10 genes was up-regulated and 1 gene was down-regulated. This study provides more information on the classification, gene structure, conserved motif, and evolution of NAC TFs in Jojoba, facilitating further exploration of their specific functional analysis in Jojoba seed development.

The NAC transcription factor family in Eucommia ulmoides: Genome-wide identification, characterization, and network analysis in relation to the rubber biosynthetic genes

The NAC transcription factor family is a large plant gene family, participating in plant growth and development, secondary metabolite synthesis, biotic and abiotic stresses responses, and hormone signaling. Eucommia ulmoides is a widely planted economic tree species in China that can produce trans-polyisoprene: Eucommia rubber (Eu-rubber). However, genome-wide identification of the NAC gene family has not been reported in E. ulmoides. In this study, 71 NAC proteins were identified based on genomic database of E. ulmoides. Phylogenetic analysis showed that the EuNAC proteins were distributed in 17 subgroups based on homology with NAC proteins in Arabidopsis, including the E. ulmoides-specific subgroup Eu_NAC. Gene structure analysis suggested that the number of exons varied from 1 to 7, and multitudinous EuNAC genes contained two or three exons. Chromosomal location analysis revealed that the EuNAC genes were unevenly distributed on 16 chromosomes. Three pairs of genes of tandem duplicates genes and 12 segmental duplications were detected, which indicated that segmental duplications may provide the primary driving force of expansion of EuNAC. Prediction of cis-regulatory elements indicated that the EuNAC genes were involved in development, light response, stress response and hormone response. For the gene expression analysis, the expression levels of EuNAC genes in various tissues were quite different. To explore the effect of EuNAC genes on Eu-rubber biosynthesis, a co-expression regulatory network between Eu-rubber biosynthesis genes and EuNAC genes was constructed, which indicated that six EuNAC genes may play an important role in the regulation of Eu-rubber biosynthesis. In addition, this six EuNAC genes expression profiles in E. ulmoides different tissues were consistent with the trend in Eu-rubber content. Quantitative real-time PCR analysis showed that EuNAC genes were responsive to different hormone treatment. These results will provide a useful reference for further studies addressing the functional characteristics of the NAC genes and its potential role in Eu-rubber biosynthesis.

Genome-Wide Characterization and Evolutionary Expansion of Poplar NAC Transcription Factors and Their Tissue-Specific Expression Profiles under Drought

The NAC (NAM, ATAF1/2 and CUC2) is a large gene family of plant-specific transcription factors that play a pivotal role in various physiological processes and abiotic stresses. Due to the lack of genome-wide characterization, intraspecific and interspecific synteny, and drought-responsive expression pattern of NAC genes in poplar, the functional characterization of drought-related NAC genes have been scarcely reported in Populus species. Here, we identified a total of 170 NAC domain-containing genes in the P. trichocarpa genome, 169 of which were unevenly distributed on its nineteen chromosomes. These NAC genes were phylogenetically divided into twenty subgroups, some of which exhibited a similar pattern of exon-intron architecture. The synteny and Ka/Ks analysis indicated that the expansion of NAC genes in poplar was mainly due to gene duplication events occurring before and after the divergence of Populus and Salix. Ten PdNAC (P. deltoids × P. euramericana cv.'Nanlin895') genes were randomly selected and cloned. Their drought-responsive expression profiles showed a tissue-specific pattern. The transcription factor PdNAC013 was verified to be localized in the nucleus. Our research results provide genomic information for the expansion of NAC genes in the poplar genome, and for further characterizing putative poplar NAC genes associated with water-deficit.

Comprehensive analysis of the NAC transcription factor gene family in Kandelia obovata reveals potential members related to chilling tolerance

BACKGROUND

Kandelia obovata is an important mangrove species extensively distributed in Eastern Asia that is susceptible to low-temperature stress. NAC (NAM, ATAF1/2 and CUC2) domain proteins are transcription factors (TFs) that play various roles in plant growth and development and in the plant response to environmental stresses. Nevertheless, genome-wide analyses of K. obovata NAC genes (KoNACs) and their responses to chilling stress have rarely been studied.

METHODS

The KoNAC gene family was identified and characterized using bioinformatic analysis, the subcellular location of some NAC proteins was confirmed using confocal microscopy analysis, and the KoNACs that responded to chilling stress were screened using RNA-seq and qRT-PCR analysis.

RESULTS

A total of 79 KoNACs were identified, and they were unequally distributed across all 18 chromosomes of K. obovata. The KoNAC proteins could be divided into 16 subgroups according to the phylogenetic tree based on NAC family members of Arabidopsis thaliana. The KoNACs exhibited greater synteny with A. thaliana sequences than with Oryza sativa sequences, indicating that KoNACs underwent extensive evolution after the divergence of dicotyledons and monocotyledons. Segmental duplication was the main driving force of the expansions of KoNAC genes. Confocal microscopy analysis verified that the four randomly selected KoNACs localized to the nucleus, indicating the accuracy of the bioinformatic predictions. Tissue expression pattern analysis demonstrated that some KoNAC genes showed tissue-specific expression, suggesting that these KoNACs might be important for plant development and growth. Additionally, the expression levels of 19 KoNACs were significantly (15 positively and 4 negatively) induced by cold treatment, demonstrating that these KoNACs might play important roles during cold stress responses and might be candidate genes for the genetic engineering of K. obovata with enhanced chilling stress tolerance. Coexpression network analysis revealed that 381 coexpressed pairs (between 13 KoNACs and 284 other genes) were significantly correlated.

CONCLUSIONS

Seventy-nine KoNACs were identified in K. obovata, nineteen of which displayed chilling-induced expression patterns. These genes may serve as candidates for functional analyses of KoNACs engaged in chilling stress. Our results lay the foundation for evolutionary analyses of KoNACs and their molecular mechanisms in response to environmental stress.

Cell wall integrity regulation across plant species

Plant cell walls are highly dynamic and chemically complex structures surrounding all plant cells. They provide structural support, protection from both abiotic and biotic stress as well as ensure containment of turgor. Recently evidence has accumulated that a dedicated mechanism exists in plants, which is monitoring the functional integrity of cell walls and initiates adaptive responses to maintain integrity in case it is impaired during growth, development or exposure to biotic and abiotic stress. The available evidence indicates that detection of impairment involves mechano-perception, while reactive oxygen species and phytohormone-based signaling processes play key roles in translating signals generated and regulating adaptive responses. More recently it has also become obvious that the mechanisms mediating cell wall integrity maintenance and pattern triggered immunity are interacting with each other to modulate the adaptive responses to biotic stress and cell wall integrity impairment. Here we will review initially our current knowledge regarding the mode of action of the maintenance mechanism, discuss mechanisms mediating responses to biotic stresses and highlight how both mechanisms may modulate adaptive responses. This first part will be focused on Arabidopsis thaliana since most of the relevant knowledge derives from this model organism. We will then proceed to provide perspective to what extent the relevant molecular mechanisms are conserved in other plant species and close by discussing current knowledge of the transcriptional machinery responsible for controlling the adaptive responses using selected examples.

FAR1/FHY3 Transcription Factors Positively Regulate the Salt and Temperature Stress Responses in Eucalyptus grandis

FAR-RED ELONGATED HYPOCOTYLS3 (FHY3) and its homolog FAR-RED IMPAIRED RESPONSE1 (FAR1), which play pivotal roles in plant growth and development, are essential for the photo-induced phyA nuclear accumulation and subsequent photoreaction. The FAR1/FHY3 family has been systematically characterized in some plants, but not in Eucalyptus grandis. In this study, genome-wide identification of FAR1/FHY3 genes in E. grandis was performed using bioinformatic methods. The gene structures, chromosomal locations, the encoded protein characteristics, 3D models, phylogenetic relationships, and promoter cis-elements were analyzed with this gene family. A total of 33 FAR1/FHY3 genes were identified in E. grandis, which were divided into three groups based on their phylogenetic relationships. A total of 21 pairs of duplicated repeats were identified by homology analysis. Gene expression analysis showed that most FAR1/FHY3 genes were differentially expressed in a spatial-specific manner. Gene expression analysis also showed that FAR1/FHY3 genes responded to salt and temperature stresses. These results and observation will enhance our understanding of the evolution and function of the FAR1/FHY3 genes in E. grandis and facilitate further studies on the molecular mechanism of the FAR1/FHY3 gene family in growth and development regulations, especially in response to salt and temperature.

Targeted re-sequencing and genome-wide association analysis for wood property traits in breeding population of Eucalyptus tereticornis × E. grandis

Globally, Eucalyptus plantations occupy 22 million ha area and is one of the preferred hardwood species due to their short rotation, rapid growth, adaptability and wood properties. In this study, we present results of GWAS in parents and 100 hybrids of Eucalyptus tereticornis × E. grandis using 762 genes presumably involved in wood formation. Comparative analysis between parents predicted 32,202 polymorphic SNPs with high average read depth of 269-562× per individual per nucleotide. Seventeen wood related traits were phenotyped across three diverse environments and GWAS was conducted using 13,610 SNPs. A total of 45 SNP-trait associations were predicted across two locations. Seven large effect markers were identified which explained more than 80% of phenotypic variation for fibre area. This study has provided an array of candidate genes which may govern fibre morphology in this genus and has predicted potential SNPs which can guide future breeding programs in tropical Eucalyptus.

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.

Systematic analysis of the Capsicum ERF transcription factor family: identification of regulatory factors involved in the regulation of species-specific metabolites

BACKGROUND

ERF transcription factors (TFs) belong to the Apetala2/Ethylene responsive Factor (AP2/ERF) TF family and play a vital role in plant growth and development processes. Capsorubin and capsaicinoids have relatively high economic and nutritional value, and they are specifically found in Capsicum. However, there is little understanding of how ERFs participate in the regulatory networks of capsorubin and capsaicinoids biosynthesis.

RESULTS

In this study, a total of 142 ERFs were identified in the Capsicum annuum genome. Subsequent phylogenetic analysis allowed us to divide ERFs into DREB (dehydration responsive element binding proteins) and ERF subfamilies, and further classify them into 11 groups with several subgroups. Expression analysis of biosynthetic pathway genes and CaERFs facilitated the identification of candidate genes related to the regulation of capsorubin and capsaicinoids biosynthesis; the candidates were focused in cluster C9 and cluster C10, as well as cluster L3 and cluster L4, respectively. The expression patterns of CaERF82, CaERF97, CaERF66, CaERF107 and CaERF101, which were found in cluster C9 and cluster C10, were consistent with those of accumulating of carotenoids (β-carotene, zeaxanthin and capsorubin) in the pericarp. In cluster L3 and cluster L4, the expression patterns of CaERF102, CaERF53, CaERF111 and CaERF92 were similar to those of the accumulating capsaicinoids. Furthermore, CaERF92, CaERF102 and CaERF111 were found to be potentially involved in temperature-mediated capsaicinoids biosynthesis.

CONCLUSION

This study will provide an extremely useful foundation for the study of candidate ERFs in the regulation of carotenoids and capsaicinoids biosynthesis in peppers.

Implementing the CRISPR/Cas9 Technology in Eucalyptus Hairy Roots Using Wood-Related Genes

Eucalypts are the most planted hardwoods worldwide. The availability of the Eucalyptus grandis genome highlighted many genes awaiting functional characterization, lagging behind because of the lack of efficient genetic transformation protocols. In order to efficiently generate knock-out mutants to study the function of eucalypts genes, we implemented the powerful CRISPR/Cas9 gene editing technology with the hairy roots transformation system. As proofs-of-concept, we targeted two wood-related genes: Cinnamoyl-CoA Reductase1 (CCR1), a key lignin biosynthetic gene and IAA9A an auxin dependent transcription factor of Aux/IAA family. Almost all transgenic hairy roots were edited but the allele-editing rates and spectra varied greatly depending on the gene targeted. Most edition events generated truncated proteins, the prevalent edition types were small deletions but large deletions were also quite frequent. By using a combination of FT-IR spectroscopy and multivariate analysis (partial least square analysis (PLS-DA)), we showed that the CCR1-edited lines, which were clearly separated from the controls. The most discriminant wave-numbers were attributed to lignin. Histochemical analyses further confirmed the decreased lignification and the presence of collapsed vessels in CCR1-edited lines, which are characteristics of CCR1 deficiency. Although the efficiency of editing could be improved, the method described here is already a powerful tool to functionally characterize eucalypts genes for both basic research and industry purposes.

Transcriptional Regulation of Sorghum Stem Composition: Key Players Identified Through Co-expression Gene Network and Comparative Genomics Analyses

Most sorghum biomass accumulates in stem secondary cell walls (SCW). As sorghum stems are used as raw materials for various purposes such as feed, energy and fiber reinforced polymers, identifying the genes responsible for SCW establishment is highly important. Taking advantage of studies performed in model species, most of the structural genes contributing at the molecular level to the SCW biosynthesis in sorghum have been proposed while their regulatory factors have mostly not been determined. Validation of the role of several MYB and NAC transcription factors in SCW regulation in Arabidopsis and a few other species has been provided. In this study, we contributed to the recent efforts made in grasses to uncover the mechanisms underlying SCW establishment. We reported updated phylogenies of NAC and MYB in 9 different species and exploited findings from other species to highlight candidate regulators of SCW in sorghum. We acquired expression data during sorghum internode development and used co-expression analyses to determine groups of co-expressed genes that are likely to be involved in SCW establishment. We were able to identify two groups of co-expressed genes presenting multiple evidences of involvement in SCW building. Gene enrichment analysis of MYB and NAC genes provided evidence that while NAC SECONDARY WALL THICKENING PROMOTING FACTOR NST genes and SECONDARY WALL-ASSOCIATED NAC DOMAIN PROTEIN gene functions appear to be conserved in sorghum, NAC master regulators of SCW in sorghum may not be as tissue compartmentalized as in Arabidopsis. We showed that for every homolog of the key SCW MYB in Arabidopsis, a similar role is expected for sorghum. In addition, we unveiled sorghum MYB and NAC that have not been identified to date as being involved in cell wall regulation. Although specific validation of the MYB and NAC genes uncovered in this study is needed, we provide a network of sorghum genes involved in SCW both at the structural and regulatory levels.

Analysis of NAC Domain Transcription Factor Genes of Tectona grandis L.f. Involved in Secondary Cell Wall Deposition

NAC proteins are one of the largest families of plant-specific transcription factors (TFs). They regulate diverse complex biological processes, including secondary xylem differentiation and wood formation. Recent genomic and transcriptomic studies of Tectona grandis L.f. (teak), one of the most valuable hardwood trees in the world, have allowed identification and analysis of developmental genes. In the present work, T. grandis NAC genes were identified and analyzed regarding to their evolution and expression profile during wood formation. We analyzed the recently published T. grandis genome, and identified 130 NAC proteins that are coded by 107 gene loci. These proteins were classified into 23 clades of the NAC family, together with Populus, Eucalyptus, and Arabidopsis. Data on transcript expression revealed specific temporal and spatial expression patterns for the majority of teak NAC genes. RT-PCR indicated expression of VND genes (Tg11g04450-VND2 and Tg15g08390-VND4) related to secondary cell wall formation in xylem vessels of 16-year-old juvenile trees. Our findings open a way to further understanding of NAC transcription factor genes in T. grandis wood biosynthesis, while they are potentially useful for future studies aiming to improve biomass and wood quality using biotechnological approaches.

Analysis of Orthologous SECONDARY WALL-ASSOCIATED NAC DOMAIN1 (SND1) Promotor Activity in Herbaceous and Woody Angiosperms

SECONDARY WALL-ASSOCIATED NAC DOMAIN1 (SND1) is a master regulator of fibre secondary wall deposition in Arabidopsis thaliana (Arabidopsis), with homologs in other angiosperms and gymnosperms. However, it is poorly understood to what extent the fibre-specific regulation of the SND1 promoter, and that of its orthologs, is conserved between diverged herbaceous and woody lineages. We performed a reciprocal reporter gene analysis of orthologous SND1 promoters from Arabidopsis (AthSND1), Eucalyptus grandis (EgrNAC61) and Populus alba × P. grandidentata (PagWND1A) relative to secondary cell wall-specific Cellulose Synthase4 (CesA4) and CesA7 promoters, in both a non-woody (Arabidopsis) and a woody (poplar) system. β-glucuronidase (GUS) reporter analysis in Arabidopsis showed that the SND1 promoter was active in vascular tissues as previously reported and showed interfascicular and xylary fibre-specific expression in inflorescence stems, while reporter constructs of the woody plant-derived promoters were partial to the (pro)cambium-phloem and protoxylem. In transgenic P. tremula × P. alba plants, all three orthologous SND1 promoters expressed the GUS reporter similarly and preferentially in developing secondary xylem, ray parenchyma and cork cambium. Ours is the first study to reciprocally test orthologous SND1 promoter specificity in herbaceous and woody species, revealing diverged regulatory functions in the herbaceous system.

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.

Hardwood Tree Genomics: Unlocking Woody Plant Biology

Woody perennial angiosperms (i.e., hardwood trees) are polyphyletic in origin and occur in most angiosperm orders. Despite their independent origins, hardwoods have shared physiological, anatomical, and life history traits distinct from their herbaceous relatives. New high-throughput DNA sequencing platforms have provided access to numerous woody plant genomes beyond the early reference genomes of Populus and Eucalyptus, references that now include willow and oak, with pecan and chestnut soon to follow. Genomic studies within these diverse and undomesticated species have successfully linked genes to ecological, physiological, and developmental traits directly. Moreover, comparative genomic approaches are providing insights into speciation events while large-scale DNA resequencing of native collections is identifying population-level genetic diversity responsible for variation in key woody plant biology across and within species. Current research is focused on developing genomic prediction models for breeding, defining speciation and local adaptation, detecting and characterizing somatic mutations, revealing the mechanisms of gender determination and flowering, and application of systems biology approaches to model complex regulatory networks underlying quantitative traits. Emerging technologies such as single-molecule, long-read sequencing is being employed as additional woody plant species, and genotypes within species, are sequenced, thus enabling a comparative ("evo-devo") approach to understanding the unique biology of large woody plants. Resource availability, current genomic and genetic applications, new discoveries and predicted future developments are illustrated and discussed for poplar, eucalyptus, willow, oak, chestnut, and pecan.

RNA-seq analysis of lignocellulose-related genes in hybrid Eucalyptus with contrasting wood basic density

BACKGROUND

Wood basic density (WBD), the biomass of plant cell walls per unit volume, is an important trait for elite tree selection in kraft pulp production. Here, we investigated the correlation between WBD and wood volumes or wood properties using 98 open-pollinated, 2.4 to 2.8 year-old hybrid Eucalyptus (Eucalyptus urophylla x E. grandis). Transcript levels of lignocellulose biosynthesis-related genes were studied.

RESULTS

The progeny plants had average WBD of 516 kg/m³ with normal distribution and did not show any correlations between WBD and wood volume or components of α-cellulose, hemicellulose and Klason lignin content. Transcriptomic analysis of two groups of five plants each with high (570-609 kg/m³) or low (378-409 kg/m³) WBD was carried out by RNA-Seq analysis with total RNAs extracted from developing xylem tissues at a breast height. Lignocellulose biosynthesis-related genes, such as cellulose synthase, invertase, cinnamate-4-hydroxylase and cinnamoyl-CoA reductase showed higher transcript levels in the high WBD group. Among plant cell wall modifying genes, increased transcript levels of several expansin and xyloglucan endo-transglycosylase/hydrolase genes were also found in high WBD plants. Interestingly, strong transcript levels of several cytoskeleton genes encoding tubulin, actin and myosin were observed in high WBD plants. Furthermore, we also found elevated transcript levels of genes encoding NAC, MYB, basic helix-loop-helix, homeodomain, WRKY and LIM transcription factors in the high WBD plants. All these results indicate that the high WBD in plants has been associated with the increased transcription of many genes related to lignocellulose formation.

CONCLUSIONS

Most lignocellulose biosynthesis related genes exhibited a tendency to transcribe at relatively higher level in high WBD plants. These results suggest that lignocellulose biosynthesis-related genes may be associated with WBD.

Genome-wide analysis of Eucalyptus grandis WRKY genes family and their expression profiling in response to hormone and abiotic stress treatment

The WRKY transcription factors, a large family of proteins in plants, are involved in multiple developmental and biological processes including response to phytohormones and abiotic stress. However, little information is available regarding the WRKY family in Eucalyptus, which has been the most widely planted hardwood trees in tropical and subtropical areas. In this study, a total of 79 WRKY genes (named as EgrWRKY1-79) were identified from the Eucalyptus grandis genome and classified into three main groups according to the phylogenetic analysis, which was further supported by their gene structure and conserved motifs. Of which, 28 EgrWRKYs were involved in tandem duplication but none for segmental duplication, indicating that tandem duplication was the main cause for the expansion of WRKY gene family in E. grandis. Subsequently, expression profiles of EgrWRKY genes in eight different tissues and in response to treatments of three hormones (SA, JA, and BR) and two abiotic stresses (salt and cold) were analyzed. The results revealed that the EgrWRKY genes had differential expression in their transcript abundance and they were differentially expressed in response to plant hormones and salt and cold stresses, suggesting their contributions to plant developmental processes as well as abiotic stresses with the involvement of hormone signaling transduction. Taken together, these findings will increase our understanding of EgrWRKY gene family involved in abiotic stresses and hormone signaling transduction, and also will provide some stress-responsive candidate EgrWRKY genes for further characterization of their functions in Eucalyptus.

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.

PpNAC1, a main regulator of phenylalanine biosynthesis and utilization in maritime pine

The transcriptional regulation of phenylalanine metabolism is particularly important in conifers, long-lived species that use large amounts of carbon in wood. Here, we show that the Pinus pinaster transcription factor, PpNAC1, is a main regulator of phenylalanine biosynthesis and utilization. A phylogenetic analysis classified PpNAC1 in the NST proteins group and was selected for functional characterization. PpNAC1 is predominantly expressed in the secondary xylem and compression wood of adult trees. Silencing of PpNAC1 in P. pinaster results in the alteration of stem vascular radial patterning and the down-regulation of several genes associated with cell wall biogenesis and secondary metabolism. Furthermore, transactivation and EMSA analyses showed that PpNAC1 is able to activate its own expression and PpMyb4 promoter, while PpMyb4 is able to activate PpMyb8, a transcriptional regulator of phenylalanine and lignin biosynthesis in maritime pine. Together, these results suggest that PpNAC1 is a functional ortholog of the ArabidopsisSND1 and NST1 genes and support the idea that key regulators governing secondary cell wall formation could be conserved between gymnosperms and angiosperms. Understanding the molecular switches controlling wood formation is of paramount importance for fundamental tree biology and paves the way for applications in conifer biotechnology.

Long cold exposure induces transcriptional and biochemical remodelling of xylem secondary cell wall in Eucalyptus

Although eucalypts are the most planted hardwood trees worldwide, the majority of them are frost sensitive. The recent creation of frost-tolerant hybrids such as Eucalyptus gundal plants (E. gunnii × E. dalrympleana hybrids), now enables the development of industrial plantations in northern countries. Our objective was to evaluate the impact of cold on the wood structure and composition of these hybrids, and on the biosynthetic and regulatory processes controlling their secondary cell-wall (SCW) formation. We used an integrated approach combining histology, biochemical characterization and transcriptomic profiling as well as gene co-expression analyses to investigate xylem tissues from Eucalyptus hybrids exposed to cold conditions. Chilling temperatures triggered the deposition of thicker and more lignified xylem cell walls as well as regulation at the transcriptional level of SCW genes. Most genes involved in lignin biosynthesis, except those specifically dedicated to syringyl unit biosynthesis, were up-regulated. The construction of a co-expression network enabled the identification of both known and potential new SCW transcription factors, induced by cold stress. These regulators at the crossroads between cold signalling and SCW formation are promising candidates for functional studies since they may contribute to the tolerance of E. gunnii × E. dalrympleana hybrids to cold.

Phosphorus application reduces aluminum toxicity in two Eucalyptus clones by increasing its accumulation in roots and decreasing its content in leaves

Under acidic conditions, aluminum (Al) toxicity is an important factor limiting plant productivity; however, the application of phosphorus (P) might alleviate the toxic effects of Al. In this study, seedlings of two vegetatively propagated Eucalyptus clones, E. grandis × E. urophylla 'G9' and E. grandis × E. urophylla 'DH32-29'were subjected to six treatments (two levels of Al stress and three levels of P). Under excessive Al stress, root Al content was higher, whereas shoot and leaf Al contents were lower with P application than those without P application. Further, Al accumulation was higher in the roots, but lower in the shoots and leaves of G9 than in those of DH32-29. The secretion of organic acids was higher under Al stress than under no Al stress. Further, under Al stress, the roots of G9 secreted more organic acids than those of DH32-29. With an increase in P supply, Al-induced secretion of organic acids from roots decreased. Under Al stress, some enzymes, including PEPC, CS, and IDH, played important roles in organic acid biosynthesis and degradation. Thus, our results indicate that P can reduce Al toxicity via the fixation of elemental Al in roots and restriction of its transport to stems and leaves, although P application cannot promote the secretion of organic acid anions. Further, the higher Al-resistance of G9 might be attributed to the higher Al accumulation in and organic acid anion secretion from roots and the lower levels of Al in leaves.

Asymmetric Evolution and Expansion of the NAC Transcription Factor in Polyploidized Cotton

Polyploidy in Gossypium hirsutum conferred different properties from its diploid ancestors under the regulation of transcription factors. The NAC transcription factor is a plant-specific family that can be related to plant growth and development. So far, little is known about the NAC family in cotton. This study identified 495 NAC genes in three cotton species and investigated the evolution and expansion of different genome-derived NAC genes in cotton. We revealed 15 distinct NAC subfamilies in cotton. Different subfamilies had different gene proportions, expansion rate, gene loss rate, and orthologous exchange rate. Paleohexaploidization (35%) and cotton-specific decaploidy (32%) might have primarily led to the expansion of the NAC family in cotton. Half of duplication events in G. hirsutum were inherited from its diploid ancestor, and others might have occurred after interspecific hybridization. In addition, NAC genes in the At and Dt subgenomes displayed asymmetric molecular evolution, as evidenced by their different gene loss rates, orthologous exchange, evolutionary rates, and expression levels. The dominant duplication event was different during the cotton evolutionary history. Different genome-derived NACs might have interacted with each other, which ultimately resulted in morphogenetic evolution. This study delineated the expansion and evolutionary history of the NAC family in cotton and illustrated the different fates of NAC genes during polyploidization.

Genome-wide analysis and expression patterns of the NAC transcription factor family in Medicago truncatula

NAC transcription factor (TF) family proteins are expressed in various developmental stages and following various stresses. NAC TFs are involved in mediating various physiological functions of plants and participate in various signaling pathways under biotic or abiotic stress. The present study provided a comprehensive functional analysis of members of the MtNAC TF family. Via screening of Medicago truncatula genome information, we identified 97 MtNAC TFs in M. truncatula and compared the phylogenetic analysis of 14 conserved groups with their Arabidopsis and rice counterparts. The NAC TFs were categorized into 14 groups based on their conserved motifs and gene structure. The predicted M. truncatula NAC genes were distributed among eight chromosomes, and in addition, we found that these genes showed mass gene duplication. Through expression profiling of RNA-seq data analysis, we determined that NAC family members were expressed significantly under different abiotic stresses. This indicates that the NAC TF shows different functions in M. truncatula. Together, this genome-wide analysis of the NAC gene family in M. truncatula, could be applied to improving stress tolerance in plants.

A New Insight into the Evolution and Functional Divergence of SWEET Transporters in Chinese White Pear (Pyrus bretschneideri)

SWEET genes are a recently identified plant gene family that play an indispensable role in sugar efflux. However, no systematic study has been performed in pear. In this research, 18 SWEET transporters identified in pear, almost twice the number found in woodland strawberry and Japanese apricot, were divided into four clades. Conserved motifs and six exons of the SWEET transporters were found in six species. SWEET transporters contained seven transmembrane segments (TMSs) that evolved from an internal duplication of an ancestral three-TMSs unit, connected by TMS4. This is the first direct evidence identifying internal repeats through bioinformatics analysis. Whole-genome duplication (WGD) or segmental duplication and dispersed duplication represent the main driving forces for SWEET family evolution in six species, with former duplications more important in pear. Gene expression results suggested that PbSWEET15 and PbSWEET17 have no expression in any tissues because of critical lost residues and that 62.5% of PbSWEET duplicate gene pairs have functional divergence. Additionally, PbSWEET6, PbSWEET7 and PbSWEET14 were found to play important roles in sucrose efflux from leaves, and the high expression of PbSWEET1 and PbSWEET2 might contribute to unloading sucrose from the phloem in the stem. Finally, PbSWEET5, PbSWEET9 and PbSWEET10 might contribute to pollen development. Overall, our study provides important insights into the evolution of the SWEET gene family in pear and four other Rosaceae, and the important candidate PbSWEET genes involved in the development of different tissues were identified in pear.

Special trends in CBF and DREB2 groups in Eucalyptus gunnii vs Eucalyptus grandis suggest that CBF are master players in the trade-off between growth and stress resistance

Annotation of the Eucalyptus grandis genome showed a large amplification of the dehydration-responsive element binding 1/C-repeat binding factor (DREB1/CBF) group without recent DREB2 gene duplication compared with other plant species. The present annotation of the CBF and DREB2 genes from a draft of the Eucalyptus gunnii genome sequence reveals at least one additional CBF copy in the E. gunnii genome compared with E. grandis, suggesting that this group is still evolving, unlike the DREB2 group. This study aims to investigate the redundancy/neo- or sub-functionalization of the duplicates and the relative involvement of the two groups in abiotic stress responses in both E. grandis and E. gunnii (lower growth but higher cold resistance). A comprehensive transcriptional analysis using high-throughput quantitative real-time polymerase chain reaction (qRT-PCR) was performed on leaves, stems and roots from the two Eucalyptus species after cold, heat or drought treatment. A large CBF cluster accounted for most of the cold response in all the organs, whereas heat and drought responses mainly involved a small CBF cluster and the DREB2 genes. In addition, CBF putative target genes, known to be involved in plant tolerance and development, were found to be cold-regulated. The higher transcript amounts of both the CBF and target genes in the cold tolerant E. gunnii contrasted with the higher CBF induction rates in the fast growing E. grandis. Altogether, the present results, in agreement with previous data about Eucalyptus transgenic lines over-expressing CBF, suggest that these factors, which promote both stress protection and growth limitation, participate in the trade-off between growth and resistance in this woody species.

Identification and expression analyses of new potential regulators of xylem development and cambium activity in cassava (Manihot esculenta)

We have identified new potential regulators of xylem cell-type determination and cellular proliferation in cassava and studied their expression in roots. Results are highly relevant for cassava biotechnology. Cassava's root system is composed of two types of root that coexist in every individual: the fibrous and the storage roots. Whether a root becomes fibrous or storage depends on the xylem cell types that it develops: fibrous roots develop xylem fibres and vessels while storage roots develop parenchyma xylem, the starch-storing tissue. A crucial question in cassava root development is how the specific xylem cell types differentiate and proliferate in the fibrous and storage roots. Using phylogenetic, protein sequence and synteny analyses we identified (1) MeVND6, MeVND7.1, MeVND7.2, MeNST3.1 and MeNST3.2 as the potential cassava orthologues of the Arabidopsis regulators of xylem cell type determination AtVND6, AtVND7 and AtNST3; and (2) MeWOX4.1 and MeWOX4.2 as the potential cassava orthologues of the Arabidopsis cambium regulator AtWOX4. Fibrous and storage roots were anatomically characterised and tested for the expression of the identified genes. Results revealed that (1) MeVND7.1 and MeVND7.2 are expressed in the fibrous but not in the storage roots; (2) MeVND6 shows low expression in both root types; (3) MeNST3.1 is not expressed in the fibrous or storage roots, while MeNST3.2 is highly expressed in both root-types and (4) MeWOX4.1 and, to a higher level, MeWOX4.2 are expressed in both the fibrous and storage roots. Results open new avenues for research in cassava root development and for food security-oriented biotechnology programmes.

Evolution of plant conducting cells: perspectives from key regulators of vascular cell differentiation

One crucial problem that plants faced during their evolution, particularly during the transition to growth on land, was how to transport water, nutrients, metabolites, and small signaling molecules within a large, multicellular body. As a solution to this problem, land plants developed specific tissues for conducting molecules, called water-conducting cells (WCCs) and food-conducting cells (FCCs). The well-developed WCCs and FCCs in extant plants are the tracheary elements and sieve elements, respectively, which are found in vascular plants. Recent molecular genetic studies revealed that transcriptional networks regulate the differentiation of tracheary and sieve elements, and that the networks governing WCC differentiation are largely conserved among land plant species. In this review, we discuss the molecular evolution of plant conducting cells. By focusing on the evolution of the key transcription factors that regulate vascular cell differentiation, the NAC transcription factor VASCULAR-RELATED NAC-DOMAIN for WCCs and the MYB-coiled-coil (CC)-type transcription factor ALTERED PHLOEM DEVELOPMENT for sieve elements, we describe how land plants evolved molecular systems to produce the specialized cells that function as WCCs and FCCs.

Biosynthesis and Metabolic Fate of Phenylalanine in Conifers

The amino acid phenylalanine (Phe) is a critical metabolic node that plays an essential role in the interconnection between primary and secondary metabolism in plants. Phe is used as a protein building block but it is also as a precursor for numerous plant compounds that are crucial for plant reproduction, growth, development, and defense against different types of stresses. The metabolism of Phe plays a central role in the channeling of carbon from photosynthesis to the biosynthesis of phenylpropanoids. The study of this metabolic pathway is particularly relevant in trees, which divert large amounts of carbon into the biosynthesis of Phe-derived compounds, particularly lignin, an important constituent of wood. The trunks of trees are metabolic sinks that consume a considerable percentage of carbon and energy from photosynthesis, and carbon is finally immobilized in wood. This paper reviews recent advances in the biosynthesis and metabolic utilization of Phe in conifer trees. Two alternative routes have been identified: the ancient phenylpyruvate pathway that is present in microorganisms, and the arogenate pathway that possibly evolved later during plant evolution. Additionally, an efficient nitrogen recycling mechanism is required to maintain sustained growth during xylem formation. The relevance of phenylalanine metabolic pathways in wood formation, the biotic interactions, and ultraviolet protection is discussed. The genetic manipulation and transcriptional regulation of the pathways are also outlined.

Eucalyptus hairy roots, a fast, efficient and versatile tool to explore function and expression of genes involved in wood formation

Eucalyptus are of tremendous economic importance being the most planted hardwoods worldwide for pulp and paper, timber and bioenergy. The recent release of the Eucalyptus grandis genome sequence pointed out many new candidate genes potentially involved in secondary growth, wood formation or lineage-specific biosynthetic pathways. Their functional characterization is, however, hindered by the tedious, time-consuming and inefficient transformation systems available hitherto for eucalypts. To overcome this limitation, we developed a fast, reliable and efficient protocol to obtain and easily detect co-transformed E. grandis hairy roots using fluorescent markers, with an average efficiency of 62%. We set up conditions both to cultivate excised roots in vitro and to harden composite plants and verified that hairy root morphology and vascular system anatomy were similar to wild-type ones. We further demonstrated that co-transformed hairy roots are suitable for medium-throughput functional studies enabling, for instance, protein subcellular localization, gene expression patterns through RT-qPCR and promoter expression, as well as the modulation of endogenous gene expression. Down-regulation of the Eucalyptus cinnamoyl-CoA reductase1 (EgCCR1) gene, encoding a key enzyme in lignin biosynthesis, led to transgenic roots with reduced lignin levels and thinner cell walls. This gene was used as a proof of concept to demonstrate that the function of genes involved in secondary cell wall biosynthesis and wood formation can be elucidated in transgenic hairy roots using histochemical, transcriptomic and biochemical approaches. The method described here is timely because it will accelerate gene mining of the genome for both basic research and industry purposes.

Comprehensive analysis of NAC transcription factors in diploid Gossypium: sequence conservation and expression analysis uncover their roles during fiber development

Determining how function evolves following gene duplication is necessary for understanding gene expansion. Transcription factors (TFs) are a class of proteins that regulate gene expression by binding to specific cis-acting elements in the promoters of target genes, subsequently activating or repressing their transcription. In the present study, we systematically examined the functional diversification of the NAC transcription factor (NAC-TFs) family by analyzing their chromosomal location, structure, phylogeny, and expression pattern in Gossypium raimondii (Gr) and G. arboreum (Ga). The 145 and 141 NAC genes identified in the Gr and Ga genomes, respectively, were annotated and divided into 18 subfamilies, which showed distinct divergence in gene structure and expression patterns during fiber development. In addition, when the functional parameters were examined, clear divergence was observed within tandem clusters, which suggested that subfunctionalization had occurred among duplicate genes. The expression patterns of homologous gene pairs also changed, suggestive of the diversification of gene function during the evolution of diploid cotton. These findings provide insights into the mechanisms underlying the functional differentiation of duplicated NAC-TFs genes in two diploid cotton species.

Structural, Evolutionary, and Functional Analysis of the Class III Peroxidase Gene Family in Chinese Pear (Pyrus bretschneideri)

Peroxidases (PRXs) are widely existed in various organisms and could be divided into different types according to their structures and functions. Specifically, the Class III Peroxidase, a plant-specific multi-gene family, involves in many physiological processes, such as the metabolism of auxin, the extension and thickening of cell wall, as well as the formation of lignin. By searching the pear genome database, 94 non-redundant PRXs from Pyrus bretschneideri (PbPRXs) were identified. Subsequently, analysis of phylogenetic relationships, gene structures, conserved motifs, and microsynteny was performed. These PbPRXs were unevenly distributed among 17 chromosomes of pear. In addition, 26 segmental duplication events but only one tandem duplication were occurred in these PbPRXs, implying segmental duplication was the main contributor to the expansion of the PbPRX family. By the Ka/Ks analysis, 26 out of 27 duplicated PbPRXs has experienced purifying selection. Twenty motifs were identified in PbPRXs based on the MEME analysis, 11 of which were enriched in pear. A total of 41 expressed genes were identified from ESTs of pear fruit. According to qRT-PCR, the expression trends of five PbPRXs in subgroup C were consistent with the change of lignin content during pear fruit development. So we inferred that the five PbPRXs were candidate genes involved in the lignin synthesis pathway. These results provided useful information for further researches of PRX genes in pear.

The NAC transcription factor family in maritime pine (Pinus Pinaster): molecular regulation of two genes involved in stress responses

BACKGROUND

NAC transcription factors comprise a large plant-specific gene family involved in the regulation of diverse biological processes. Despite the growing number of studies on NAC transcription factors in various species, little information is available about this family in conifers. The goal of this study was to identify the NAC transcription family in maritime pine (Pinus pinaster), to characterize ATAF-like genes in response to various stresses and to study their molecular regulation.

METHODS

We have isolated two maritime pine NAC genes and using a transient expression assay in N. benthamiana leaves estudied the promoter jasmonate response.

RESULTS

In this study, we identified 37 NAC genes from maritime pine and classified them into six main subfamilies. The largest group includes 12 sequences corresponding to stress-related genes. Two of these NAC genes, PpNAC2 and PpNAC3, were isolated and their expression profiles were examined at various developmental stages and in response to various types of stress. The expression of both genes was strongly induced by methyl jasmonate (MeJA), mechanical wounding, and high salinity. The promoter regions of these genes were shown to contain cis-elements involved in the stress response and plant hormonal regulation, including E-boxes, which are commonly found in the promoters of genes that respond to jasmonate, and binding sites for bHLH proteins. Using a transient expression assay in N. benthamiana leaves, we found that the promoter of PpNAC3 was rapidly induced upon MeJA treatment, while this response disappeared in plants in which the transcription factor NbbHLH2 was silenced.

CONCLUSION

Our results suggest that PpNAC2 and PpNAC3 encode stress-responsive NAC transcription factors involved in the jasmonate response in pine. Furthermore, these data also suggest that the jasmonate signaling pathway is conserved between angiosperms and gymnosperms. These findings may be useful for engineering stress tolerance in pine via biotechnological approaches.

Genome-wide analysis of the lignin toolbox of Eucalyptus grandis

Lignin, a major component of secondary cell walls, hinders the optimal processing of wood for industrial uses. The recent availability of the Eucalyptus grandis genome sequence allows comprehensive analysis of the genes encoding the 11 protein families specific to the lignin branch of the phenylpropanoid pathway and identification of those mainly involved in xylem developmental lignification. We performed genome-wide identification of putative members of the lignin gene families, followed by comparative phylogenetic studies focusing on bona fide clades inferred from genes functionally characterized in other species. RNA-seq and microfluid real-time quantitative PCR (RT-qPCR) expression data were used to investigate the developmental and environmental responsive expression patterns of the genes. The phylogenetic analysis revealed that 38 E. grandis genes are located in bona fide lignification clades. Four multigene families (shikimate O-hydroxycinnamoyltransferase (HCT), p-coumarate 3-hydroxylase (C3H), caffeate/5-hydroxyferulate O-methyltransferase (COMT) and phenylalanine ammonia-lyase (PAL)) are expanded by tandem gene duplication compared with other plant species. Seventeen of the 38 genes exhibited strong, preferential expression in highly lignified tissues, probably representing the E. grandis core lignification toolbox. The identification of major genes involved in lignin biosynthesis in E. grandis, the most widely planted hardwood crop world-wide, provides the foundation for the development of biotechnology approaches to develop tree varieties with enhanced processing qualities.

NAC-MYB-based transcriptional regulation of secondary cell wall biosynthesis in land plants

Plant cells biosynthesize primary cell walls (PCW) in all cells and produce secondary cell walls (SCWs) in specific cell types that conduct water and/or provide mechanical support, such as xylem vessels and fibers. The characteristic mechanical stiffness, chemical recalcitrance, and hydrophobic nature of SCWs result from the organization of SCW-specific biopolymers, i.e., highly ordered cellulose, hemicellulose, and lignin. Synthesis of these SCW-specific biopolymers requires SCW-specific enzymes that are regulated by SCW-specific transcription factors. In this review, we summarize our current knowledge of the transcriptional regulation of SCW formation in plant cells. Advances in research on SCW biosynthesis during the past decade have expanded our understanding of the transcriptional regulation of SCW formation, particularly the functions of the NAC and MYB transcription factors. Focusing on the NAC-MYB-based transcriptional network, we discuss the regulatory systems that evolved in land plants to modify the cell wall to serve as a key component of structures that conduct water and provide mechanical support.