Characterization of NAC domain transcription factors implicated in control of vascular cell differentiation in Arabidopsis and Populus

Growth-regulating factor 15-mediated vascular cambium differentiation positively regulates wood formation in hybrid poplar (Populus alba × P. glandulosa)

Introduction

Hybrid poplars are industrial trees in China. An understanding of the molecular mechanism underlying wood formation in hybrid poplars is necessary for molecular breeding. Although the division and differentiation of vascular cambial cells is important for secondary growth and wood formation, the regulation of this process is largely unclear.

Methods

In this study, mPagGRF15 OE and PagGRF15-SRDX transgenic poplars were generated to investigate the function of PagGRF15. RNA-seq and qRT-PCR were conducted to analyze genome-wide gene expression, while ChIP‒seq and ChIP-PCR were used to identified the downstream genes regulated by PagGRF15.

Results and discussion

We report that PagGRF15 from hybrid poplar (Populus alba × P. glandulosa), a growth-regulating factor, plays a critical role in the regulation of vascular cambium activity. PagGRF15 was expressed predominantly in the cambial zone of vascular tissue. Overexpression of mPagGRF15 (the mutated version of GRF15 in the miR396 target sequence) in Populus led to decreased plant height and internode number. Further stem cross sections showed that the mPagGRF15 OE plants exhibited significant changes in vascular pattern with an increase in xylem and a reduction in phloem. In addition, cambium cell files were decreased in the mPagGRF15 OE plants. However, dominant suppression of the downstream genes of PagGRF15 using PagGRF15-SRDX showed an opposite phenotype. Based on the RNA-seq and ChIP-seq results, combining qRT-PCR and ChIP-PCR analysis, candidate genes, such as WOX4b, PXY and GID1.3, were obtained and found to be mainly involved in cambial activity and xylem differentiation. Accordingly, we speculated that PagGRF15 functions as a positive regulator mediating xylem differentiation by repressing the expression of the WOX4a and PXY genes to set the pace of cambial activity. In contrast, PagGRF15 mediated the GA signaling pathway by upregulating GID1.3 expression to stimulate xylem differentiation. This study provides valuable information for further studies on vascular cambium differentiation mechanisms and genetic improvement of the specific gravity of wood in hybrid poplars.

PbBPC4 involved in a xylem-deficient dwarf phenotype in pear by directly regulating the expression of PbXND1

Dwarfing is an important agronomic trait in fruit breeding. At present, dwarf cultivars or dwarfing rootstocks are used for high-density planting. Although some dwarf rootstocks have been used in the cultivation of pear (Pyrus bretschneideri Rehd), the breeding of dwarf pear rootstocks or cultivars is still sorely lacking. A previous study reported that PbXND1 results in a xylem-dwarf phenotype in pear trees. However, the regulatory mechanism upstream of PbXND1 is unclear. In this study, we identified PbBPC4 as an upstream regulatory factor of PbXND1 in yeast one-hybrid assays. In β-glucuronidase staining and dual-luciferase assays, PbBPC4 enhanced the activity of the PbXND1 promoter. Tobacco plants overexpressing PbBPC4 showed decreased plant height because of a reduced xylem size. Similar changes in the xylem was observed in transgenic pear roots; those overexpressing PbBPC4 showed reduced xylem size, and those with silencing PbBPC4 expression showed increased xylem size, greater density of xylem vessels, and a larger proportion of the xylem out of the total cross-section area. Expression analyses showed that PbBPC4 increases the transcription of PbXND1, leading to reduced transcript levels of genes involved in the positive regulation of xylem development, ultimately resulting in a xylem-deficient dwarf phenotype. Taken together, our results reveal the mechanism by which PbBPC4 participates in the regulation of xylem development via directly altering the expression of PbXND1, thus leading to the dwarf phenotype in pear. These findings have reference value for the breeding of dwarf pear trees.

The transcription factor ERF108 interacts with AUXIN RESPONSE FACTORs to mediate cotton fiber secondary cell wall biosynthesis

Phytohormones play indispensable roles in plant growth and development. However, the molecular mechanisms underlying phytohormone-mediated regulation of fiber secondary cell wall (SCW) formation in cotton (Gossypium hirsutum) remain largely underexplored. Here, we provide mechanistic evidence for functional interplay between the APETALA2/ethylene response factor (AP2/ERF) transcription factor GhERF108 and auxin response factors GhARF7-1 and GhARF7-2 in dictating the ethylene-auxin signaling crosstalk that regulates fiber SCW biosynthesis. Specifically, in vitro cotton ovule culture revealed that ethylene and auxin promote fiber SCW deposition. GhERF108 RNA interference (RNAi) cotton displayed remarkably reduced cell wall thickness compared with controls. GhERF108 interacted with GhARF7-1 and GhARF7-2 to enhance the activation of the MYB transcription factor gene GhMYBL1 (MYB domain-like protein 1) in fibers. GhARF7-1 and GhARF7-2 respond to auxin signals that promote fiber SCW thickening. GhMYBL1 RNAi and GhARF7-1 and GhARF7-2 virus-induced gene silencing (VIGS) cotton displayed similar defects in fiber SCW formation as GhERF108 RNAi cotton. Moreover, the ethylene and auxin responses were reduced in GhMYBL1 RNAi plants. GhMYBL1 directly binds to the promoters of GhCesA4-1, GhCesA4-2, and GhCesA8-1 and activates their expression to promote cellulose biosynthesis, thereby boosting fiber SCW formation. Collectively, our findings demonstrate that the collaboration between GhERF108 and GhARF7-1 or GhARF7-2 establishes ethylene-auxin signaling crosstalk to activate GhMYBL1, ultimately leading to the activation of fiber SCW biosynthesis.

Transcriptional reprogramming during recovery from drought stress in Eucalyptus grandis

The importance of drought as a constraint to agriculture and forestry is increasing with climate change. Genetic improvement of plants' resilience is one of the mitigation strategies to curb this threat. Although recovery from drought stress is important to long-term drought adaptation and has been considered as an indicator of dehydration tolerance in annual crops, this has not been well explored in forest trees. Thus, we aimed to investigate the physiological and transcriptional changes during drought stress and rewatering in Eucalyptus grandis W. Hill ex Maiden. We set up a greenhouse experiment where we imposed drought stress on 2-year-old seedlings and rewatered the recovery group after 17 days of drought. Our measurement of leaf stomatal conductance (gs) showed that, while gs was reduced by drought stress, it fully recovered after 5 days of rewatering. The RNA-seq analysis from stem samples revealed that genes related to known stress responses such as phytohormone and reactive oxygen species signaling were upregulated, while genes involved in metabolism and growth were downregulated due to drought stress. We observed reprogramming of signal transduction pathways and metabolic processes at 1 day of rewatering, indicating a quick response to rewatering. Our results suggest that recovery from drought stress may entail alterations in the jasmonic acid, salicylic acid, ethylene and brassinosteroid signaling pathways. Using co-expression network analysis, we identified hub genes, including the putative orthologs of ABI1, ABF2, ABF3, HAI2, BAM1, GolS2 and SIP1 during drought and CAT2, G6PD1, ADG1 and FD-1 during recovery. Taken together, by highlighting the molecular processes and identifying key genes, this study gives an overview of the mechanisms underlying the response of E. grandis to drought stress and recovery that trees may face repeatedly throughout their long life cycle. This provides a useful reference to the identification and further investigation of signaling pathways and target genes for future tree improvement.

ZmHOX32 is related to photosynthesis and likely functions in plant architecture of maize

HOX32, a member of the HD-ZIP III family, functions in the leaf morphogenesis and plant photosynthesis. However, the regulatory mechanism of HOX32 in maize has not been studied and the regulatory relationship in photosynthesis is unclear. We conducted a comprehensive study, including phylogenetic analysis, expression profiling at both transcriptome and translatome levels, subcellular localization, tsCUT&Tag, co-expression analysis, and association analysis with agronomic traits on HOX32 for the dissection of the functional roles of HOX32. ZmHOX32 shows conservation in plants. As expected, maize HOX32 protein is specifically expressed in the nucleus. ZmHOX32 showed constitutively expression at both transcriptome and translatome levels. We uncovered the downstream target genes of ZmHOX32 by tsCUT&Tag and constructed a cascaded regulatory network combining the co-expression networks. Both direct and indirect targets of ZmHOX32 showed significant gene ontology enrichment in terms of photosynthesis in maize. The association study suggested that ZmHOX32 plays an important role in regulation of plant architecture. Our results illustrate a complex regulatory network of HOX32 involving in photosynthesis and plant architecture, which deepens our understanding of the phenotypic variation in plants.

Current Understanding of the Genetics and Molecular Mechanisms Regulating Wood Formation in Plants

Unlike herbaceous plants, woody plants undergo volumetric growth (a.k.a. secondary growth) through wood formation, during which the secondary xylem (i.e., wood) differentiates from the vascular cambium. Wood is the most abundant biomass on Earth and, by absorbing atmospheric carbon dioxide, functions as one of the largest carbon sinks. As a sustainable and eco-friendly energy source, lignocellulosic biomass can help address environmental pollution and the global climate crisis. Studies of Arabidopsis and poplar as model plants using various emerging research tools show that the formation and proliferation of the vascular cambium and the differentiation of xylem cells require the modulation of multiple signals, including plant hormones, transcription factors, and signaling peptides. In this review, we summarize the latest knowledge on the molecular mechanism of wood formation, one of the most important biological processes on Earth.

Functional Diversification of Populus FLOWERING LOCUS D-LIKE3 Transcription Factor and Two Paralogs in Shoot Ontogeny, Flowering, and Vegetative Phenology

Both the evolution of tree taxa and whole-genome duplication (WGD) have occurred many times during angiosperm evolution. Transcription factors are preferentially retained following WGD suggesting that functional divergence of duplicates could contribute to traits distinctive to the tree growth habit. We used gain- and loss-of-function transgenics, photoperiod treatments, and circannual expression studies in adult trees to study the diversification of three Populus FLOWERING LOCUS D-LIKE (FDL) genes encoding bZIP transcription factors. Expression patterns and transgenic studies indicate that FDL2.2 promotes flowering and that FDL1 and FDL3 function in different vegetative phenophases. Study of dominant repressor FDL versions indicates that the FDL proteins are partially equivalent in their ability to alter shoot growth. Like its paralogs, FDL3 overexpression delays short day-induced growth cessation, but also induces distinct heterochronic shifts in shoot development-more rapid phytomer initiation and coordinated delay in both leaf expansion and the transition to secondary growth in long days, but not in short days. Our results indicate that both regulatory and protein coding sequence variation contributed to diversification of FDL paralogs that has led to a degree of specialization in multiple developmental processes important for trees and their local adaptation.

PagGRF12a interacts with PagGIF1b to regulate secondary xylem development through modulating PagXND1a expression in Populus alba × P. glandulosa

Growth-regulating factors (GRFs) are important regulators of plant development and growth, but their possible roles in xylem development in woody plants remain unclear. Here, we report that Populus alba × Papulus glandulosa PagGRF12a negatively regulates xylem development in poplar. PagGRF12a is expressed in vascular tissues. Compared to non-transgenic control plants, transgenic poplar plants overexpressing PagGRF12a exhibited reduced xylem width and plants with repressed expression of PagGRF12a exhibited increased xylem width. Xylem NAC domain 1 (XND1) encodes a NAC domain transcription factor that regulates xylem development and transcriptional analyses revealed that PagXND1a is highly upregulated in PagGRF12a-overexpressing plants and downregulated in PagGRF12a-suppressed plants, indicating that PagGRF12a may regulate xylem development through PagXND1a. Transient transcriptional assays and chromatin immunoprecipitation-polymerase chain reaction assays confirmed that PagGRF12a directly upregulates PagXND1a. In addition, PagGRF12a interacts with the GRF-Interacting Factor (GIF) PagGIF1b, and this interaction enhances the effects of PagGRF12a on PagXND1a. Our results indicate that PagGRF12a inhibits xylem development by upregulating the expression of PagXND1a.

Chia (Salvia hispanica) Gene Expression Atlas Elucidates Dynamic Spatio-Temporal Changes Associated With Plant Growth and Development

Chia (Salvia hispanica L.), now a popular superfood and a pseudocereal, is one of the richest sources of dietary nutrients such as protein, fiber, and polyunsaturated fatty acids (PUFAs). At present, the genomic and genetic information available in the public domain for this crop are scanty, which hinders an understanding of its growth and development and genetic improvement. We report an RNA-sequencing (RNA-Seq)-based comprehensive transcriptome atlas of Chia sampled from 13 tissue types covering vegetative and reproductive growth stages. We used ~355 million high-quality reads of total ~394 million raw reads from transcriptome sequencing to generate de novo reference transcriptome assembly and the tissue-specific transcript assemblies. After the quality assessment of the merged assemblies and implementing redundancy reduction methods, 82,663 reference transcripts were identified. About 65,587 of 82,663 transcripts were translated into 99,307 peptides, and we were successful in assigning InterPro annotations to 45,209 peptides and gene ontology (GO) terms to 32,638 peptides. The assembled transcriptome is estimated to have the complete sequence information for ~86% of the genes found in the Chia genome. Furthermore, the analysis of 53,200 differentially expressed transcripts (DETs) revealed their distinct expression patterns in Chia's vegetative and reproductive tissues; tissue-specific networks and developmental stage-specific networks of transcription factors (TFs); and the regulation of the expression of enzyme-coding genes associated with important metabolic pathways. In addition, we identified 2,411 simple sequence repeats (SSRs) as potential genetic markers from the transcripts. Overall, this study provides a comprehensive transcriptome atlas, and SSRs, contributing to building essential genomic resources to support basic research, genome annotation, functional genomics, and molecular breeding of Chia.

Genome-wide identification and comprehensive analysis of NAC family genes involved in fruit development in kiwifruit (Actinidia)

BACKGROUND

NAC transcription factors (TFs) are plant-specific proteins encoded by a large gene family. They play important roles in diverse biological processes, such as plant growth and development, leaf senescence, and responses to biotic or abiotic stresses. Functions of a number of NAC TFs have been identified mainly in model plants. However, very few studies on NAC TFs have been conducted in the fruit tree of kiwifruit.

RESULTS

Genome-wide NAC genes were identified and their phylogeny, genomic structure, chromosomal location, synteny relationships, protein properties and conserved motifs were analyzed. In addition, the fruit developmental process was evaluated in a new kiwifruit cultivar of Actinidia eriantha 'Ganlu 1'. And expressions for all those NAC genes were analyzed by quantitative real-time PCR method in fruits of 'Ganlu 1' during its developmental process. Our research identified 142 NAC TFs which could be phylogenetically divided into 23 protein subfamilies. The genomic structures of those NAC genes indicated that their exons were between one and ten. Analysis of chromosomal locations suggested that 116 out of 142 NACs distributed on all the 29 kiwifruit chromosomes. In addition, genome-wide gene expression analysis showed that expressions of 125 out of 142 NAC genes could be detected in fruit samples.

CONCLUSION

Our comprehensive study provides novel information on NAC genes and expression patterns in kiwifruit fruit. This research would be helpful for future functional identification of NAC genes involved in kiwifruit fruit development.

Identification and Characterization of Secondary Wall-Associated NAC Genes and Their Involvement in Hormonal Responses in Tobacco (Nicotiana tabacum)

Secondary wall-associated NAC (SWN) genes are a subgroup of NAC (NAM, ATAF, and CUC) transcription factors (TF) that play a key role in regulating secondary cell wall biosynthesis in plants. However, this gene family has not been systematically characterized, and their potential roles in response to hormones are unknown in Nicotiana tabacum. In this study, a total of 40 SWN genes, of which 12 from Nicotiana tomentosiformis, 13 from Nicotiana sylvestris, and 15 from Nicotiana tabacum, were successfully identified. The 15 SWNs from Nicotiana tabacum were further classified into three groups, namely, vascular-related NAC domain genes (NtVNDs), NAC secondary wall thickening promoting factor genes (NtNSTs), and secondary wall-associated NAC domain genes (NtSNDs). The protein characteristic, gene structure, and chromosomal location of 15 NtSWNs (also named Nt1 to Nt15) were also analyzed. The NtVND and NtNST group genes had five conserved subdomains in their N-terminal regions and a motif (LP[Q/x] L[E/x] S[P/A]) in their diverged C- terminal regions. Some hormones, dark and low-temperature related cis-acting elements, were significantly enriched in the promoters of NtSWN genes. A comprehensive expression profile analysis revealed that Nt4 and Nt12 might play a role in vein development. Others might be important for stem development. Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) revealed that in the NtNST group, genes such as Nt7, Nt8, and Nt13 were more sensitive than the genes in NtVND and NtSND groups under abiotic stress conditions. A transactivation assay further suggested that Nt7, Nt8, and Nt13 showed a significant transactivation activity. Overall, SWN genes were finally identified and characterized in diploid and tetraploid tobacco, revealing new insights into their evolution, variation, and homology relationships. Transcriptome, cis-acting element, qRT-PCR, and transactivation assay analysis indicated the roles in hormonal and stress responses, which provided further resources in molecular mechanism and genetic improvement.

Transcriptome analysis of salt-responsive and wood-associated NACs in Populus simonii × Populus nigra

BACKGROUND

NAC (NAM, ATAF1-2, and CUC2) family is one of the largest plant-specific transcription factor families known to play significant roles in plant development processes and stress responses.

RESULTS

In the study, a total of 112 NACs were identified to be differentially expressed in the comparisons of leaves and stems, leaves and roots, roots and stems of Populus simonii×P. nigra among 289 members by RNA-Seq. And 148, 144 and 134 NACs were detected to be salt-responsive in the roots, stems and leaves under 150 mM NaCl stress, respectively. Among them, a total of 53 salt-responsive NACs were shared across the three tissues. Under salt stress, 41/37 NACs were identified to be up/down-regulated in the leaves of Populus simonii × P.nigra among 170 non-redundant NACs by RT-qPCR, which was similar with RNA-Seq results. The expression pattern analysis of 6 NACs including four randomly up-regulated genes (NAC86, NAC105, NAC139 and NAC163) and two down-regulated genes (NAC15 and NAC149) indicated a few NACs showed specific temporal and spatial expression patterns in the three tissues of Populus simonii×P.nigra. Based on transcriptome screening and phylogenic analysis of differentially expressed NACs in different tissues under salt stress, 18 potential NACs associated with wood formation and 20 involved in stress responses were identified in Populus simonii×P.nigra.

CONCLUSIONS

The study further gains an understanding of the connection of tissue specificity and gene function in poplar, and lays the foundation of functional analysis of poplar NACs in stress responses.

Insights of Molecular Mechanism of Xylem Development in Five Black Poplar Cultivars

Black poplar (Populus deltoides, P. nigra, and their hybrids) is the main poplar cultivars in China. It offers interesting options of large-scale biomass production for bioenergy due to its rapid growth and high yield. Poplar wood properties were associated with chemical components and physical structures during wood formation. In this study, five poplar cultivars, P. euramericana 'Zhonglin46' (Pe1), P. euramericana 'Guariento' (Pe2), P. nigra 'N179' (Pn1), P. deltoides 'Danhong' (Pd1), and P. deltoides 'Nanyang' (Pd2), were used to explore the molecular mechanism of xylem development. We analyzed the structural differences of developing xylem in the five cultivars and profiled the transcriptome-wide gene expression patterns through RNA sequencing. The cross sections of the developing xylem showed that the cell wall thickness of developed fiber in Pd1 was thickest and the number of xylem vessels of Pn1 was the least. A total of 10,331 differentially expressed genes were identified among 10 pairwise comparisons of the five cultivars, most of them were related to programmed cell death and secondary cell wall thickening. K-means cluster analysis and Gene Ontology enrichment analysis showed that the genes highly expressed in Pd1 were related to nucleotide decomposition, metabolic process, transferase, and microtubule cytoskeleton; whereas the genes highly expressed in Pn1 were involved in cell wall macromolecule decomposition and polysaccharide binding processes. Based on a weighted gene co-expression network analysis, a large number of candidate regulators for xylem development were identified. And their potential regulatory roles to cell wall biosynthesis genes were validated by a transient overexpression system. This study provides a set of promising candidate regulators for genetic engineering to improve feedstock and enhance biofuel conversion in the bioenergy crop Populus.

KNAT2/6b, a class I KNOX gene, impedes xylem differentiation by regulating NAC domain transcription factors in poplar

Wood formation is the terminal differentiation of xylem mother cells derived from cambial initials, and negative regulators play important roles in xylem differentiation. The molecular mechanism of the negative regulator of xylem differentiation PagKNAT2/6b was investigated. PagKNAT2/6b is an ortholog of Arabidopsis KNAT2 and KNAT6 that is highly expressed in phloem and xylem. Compared to nontransgenic control plants, transgenic poplar plants overexpressing PagKNAT2/6b present with altered vascular patterns, characterized by decreased secondary xylem with thin cell walls containing less cellulose, xylose and lignin. RNA sequencing analyses revealed that differentially expressed genes are enriched in xylem differentiation and secondary wall synthesis functions. Expression of NAM/ATAF/CUC (NAC) domain genes including PagSND1-A1, PagSND1-A2, PagSND1-B2 and PagVND6-C1 is downregulated by PagKNAT2/6b, while PagXND1a is directly upregulated. Accordingly, the dominant repression form of PagKNAT2/6b leads to increased xylem width per stem diameter through downregulation of PagXND1a. PagKNAT2/6b can inhibit cell differentiation and secondary wall deposition during wood formation in poplar by modulating the expression of NAC domain transcription factors. Direct activation of PagXND1a by PagKNAT2/6b is a key node in the negative regulatory network of xylem differentiation by KNOXs.

Over-expression of poplar NAC15 gene enhances wood formation in transgenic tobacco

BACKGROUND

NAC (NAM/ATAF/CUC) is one of the largest plant-specific transcription factor (TF) families known to play significant roles in wood formation. Acting as master gene regulators, a few NAC genes can activate secondary wall biosynthesis during wood formation in woody plants.

RESULTS

In the present study, firstly, we screened 110 differentially expressed NAC genes in the leaves, stems, and roots of di-haploid Populus simonii×P. nigra by RNA-Seq. Then we identified a nucleus-targeted gene, NAC15 gene, which was one of the highly expressed genes in the stem among 110 NAC family members. Thirdly, we conducted expression pattern analysis of NAC15 gene, and observed NAC15 gene was most highly expressed in the xylem by RT-qPCR. Moreover, we transferred NAC15 gene into tobacco and obtained 12 transgenic lines overexpressing NAC15 gene (TLs). And the relative higher content of hemicellulose, cellulose and lignin was observed in the TLs compared to the control lines containing empty vector (CLs). It also showed darker staining in the culms of the TLs with phloroglucinol staining, compared to the CLs. Furthermore, the relative expression level of a few lignin- and cellulose-related genes was significantly higher in the TLs than that in the CLs.

CONCLUSIONS

The overall results indicated that NAC15 gene is highly expressed in the xylem of poplar and may be a potential candidate gene playing an important role in wood formation in transgenic tobacco.

Induction of PrMADS10 on the lower side of bent pine tree stems: potential role in modifying plant cell wall properties and wood anatomy

The molecular mechanisms underlying inclination responses in trees are unclear. In this study, we identified a MADS-box transcription factor differentially expressed early after inclination in the stems of Pinus radiata D. Don. PrMADS10 has a CDS of 582 bp and encodes a group II MADS-box transcription factor. We measured highest accumulation of this transcript on the lower side of inclined pine stems. In an effort to identify putative targets, we stably transformed Arabidopsis thaliana with a 35S::PrMADS10 construct. Transcriptome analysis revealed 1,219 genes differentially-expressed, with 690 and 529 genes up- and down-regulated respectively, when comparing the transgenic and wild-type. Differentially-expressed genes belong to different biological processes, but were enriched in cell wall remodeling and phenylpropanoid metabolic functions. Interestingly, lignin content was 30% higher in transgenic as compared to wild-type plants consistent with observed changes in gene expression. Differentially expressed transcription factors and phenylpropanoid genes were analyzed using STRING. Several MYB and NAC transcription factors showed interactions with genes of the phenylpropanoid pathway. Together, these results implicate PrMADS10 as a regulatory factor, triggering the expression of other transcription factors and genes involved in the synthesis of lignin.

Weighted gene co-expression network analysis unveils gene networks associated with the Fusarium head blight resistance in tetraploid wheat

BACKGROUND

Fusarium head blight (FHB) resistance in the durum wheat breeding gene pool is rarely reported. Triticum turgidum ssp. carthlicum line Blackbird is a tetraploid relative of durum wheat that offers partial FHB resistance. Resistance QTL were identified for the durum wheat cv. Strongfield × Blackbird population on chromosomes 1A, 2A, 2B, 3A, 6A, 6B and 7B in a previous study. The objective of this study was to identify the defense mechanisms underlying the resistance of Blackbird and report candidate regulator defense genes and single nucleotide polymorphism (SNP) markers within these genes for high-resolution mapping of resistance QTL reported for the durum wheat cv. Strongfield/Blackbird population.

RESULTS

Gene network analysis identified five networks significantly (P < 0.05) associated with the resistance to FHB spread (Type II FHB resistance) one of which showed significant correlation with both plant height and relative maturity traits. Two gene networks showed subtle differences between Fusarium graminearum-inoculated and mock-inoculated plants, supporting their involvement in constitutive defense. The candidate regulator genes have been implicated in various layers of plant defense including pathogen recognition (mainly Nucleotide-binding Leucine-rich Repeat proteins), signaling pathways including the abscisic acid and mitogen activated protein (MAP) kinase, and downstream defense genes activation including transcription factors (mostly with dual roles in defense and development), and cell death regulator and cell wall reinforcement genes. The expression of five candidate genes measured by quantitative real-time PCR was correlated with that of RNA-seq, corroborating the technical and analytical accuracy of RNA-sequencing.

CONCLUSIONS

Gene network analysis allowed identification of candidate regulator genes and genes associated with constitutive resistance, those that will not be detected using traditional differential expression analysis. This study also shed light on the association of developmental traits with FHB resistance and partially explained the co-localization of FHB resistance with plant height and maturity QTL reported in several previous studies. It also allowed the identification of candidate hub genes within the interval of three previously reported FHB resistance QTL for the Strongfield/Blackbird population and associated SNPs for future high resolution mapping studies.

Characterization of a NAC transcription factor involved in the regulation of pomegranate seed hardness (Punica granatum L.)

The pomegranate, Punica granatum L., which has been cultivated since antiquity, is known to be a superfruit, possessing an array of functional anti-oxidants and various other health benefits. The hardness of pomegranate seeds is an important indicator of fruit quality, which in turn affects economic value and market demand. However, the molecular mechanism underlying pomegranate seed hardness remains to be fully understood. In this study, we found a positive correlation between seed hardness and lignin content in two pomegranate varieties: "Tunisia" and "Sanbai". Specifically, genes associated with lignin biosynthesis were differentially expressed in soft-seed and hard-seed pomegranate varieties. Among these differential genes, we cloned and characterized the NAC transcription factor PgSND1-like. Sequence alignment found a single base replacement at the 166-bp position of CDS in the PgSND1-like gene from "Tunisia" and "Sanbai". Both PgSND1-like (Sanbai) and PgSND1-like (Tunisia) proteins are localized in the cell nucleus and have a transcription activation domain in the C-terminus. Yeast two-hybrid analysis indicated that PgSND1-like protein interacts with itself to form a homodimer. Overexpression of PgSND1-like (Sanbai) in Arabidopsis showed a higher lignin content in inflorescence stem and mature seed compared with wild-type Arabidopsis. Accordingly, the expression levels of several lignin biosynthesis-associated genes were upregulated in stem cells and mature seeds of transgenic plants. However, PgSND1-like (Tunisia) transgenic Arabidopsis showed no phenotypic differences with wild-type Arabidopsis. Taken together, we suggest that PgSND1-like may regulate at least two different functions in two pomegranate varieties, promoting lignin biosynthesis and seed hardness of pomegranate.

Branching out: new insights into the genetic regulation of shoot architecture in trees

Directional growth in all plants involves both phototropic and gravitropic responses. Accordingly, mechanisms controlling shoot architecture throughout the plant kingdom are likely similar. However, as forms vary between species due in part to gene copy number and functional divergence, some aspects of how plants predetermine and regulate architecture can differ. This is especially true when comparing annual herbaceous species (e.g. model plants) to woody perennials such as trees. In the past decade, inexpensive genomic sequencing and technological advances enabled gene discovery and functional analyses in trees. This led to the identification of genes associated with tree shoot architecture control. Here, we present recent discoveries on the regulation of shoot architectures for which causative genes have been identified, including dwarf, weeping, columnar, and pillar growth habits. We also discuss potential applications of these findings.

Natural variation at XND1 impacts root hydraulics and trade-off for stress responses in Arabidopsis

Soil water uptake by roots is a key component of plant performance and adaptation to adverse environments. Here, we use a genome-wide association analysis to identify the XYLEM NAC DOMAIN 1 (XND1) transcription factor as a negative regulator of Arabidopsis root hydraulic conductivity (Lp_r). The distinct functionalities of a series of natural XND1 variants and a single nucleotide polymorphism that determines XND1 translation efficiency demonstrate the significance of XND1 natural variation at species-wide level. Phenotyping of xnd1 mutants and natural XND1 variants show that XND1 modulates Lp_r through action on xylem formation and potential indirect effects on aquaporin function and that it diminishes drought stress tolerance. XND1 also mediates the inhibition of xylem formation by the bacterial elicitor flagellin and counteracts plant infection by the root pathogen Ralstonia solanacearum. Thus, genetic variation at XND1, and xylem differentiation contribute to resolving the major trade-off between abiotic and biotic stress resistance in Arabidopsis.

Soil water uptake is a major determinant of plant performance and stress tolerance. Here the authors show that, by affecting xylem formation in the root, natural variation at the Arabidopsis XND1 locus has contrasting effects on root hydraulics and drought tolerance versus pathogen resistance.

Overexpression of a Novel Apple NAC Transcription Factor Gene, MdNAC1, Confers the Dwarf Phenotype in Transgenic Apple (Malus domestica)

Plant height is an important trait for fruit trees. The dwarf characteristic is commonly associated with highly efficient fruit production, a major objective when breeding for apple (Malus domestica). We studied the function of MdNAC1, a novel NAC transcription factor (TF) gene in apple related to plant dwarfing. Localized primarily to the nucleus, MdNAC1 has transcriptional activity in yeast cells. Overexpression of the gene results in a dwarf phenotype in transgenic apple plants. Their reduction in size is manifested by shorter, thinner stems and roots, and a smaller leaf area. The transgenics also have shorter internodes and fewer cells in the stems. Levels of endogenous abscisic acid (ABA) and brassinosteroid (BR) are lower in the transgenic plants, and expression is decreased for genes involved in the biosynthesis of those phytohormones. All of these findings demonstrate that MdNAC1 has a role in plants dwarfism, probably by regulating ABA and BR production.

Identification of new protein-protein and protein-DNA interactions linked with wood formation in Populus trichocarpa

Cellular processes, such as signal transduction and cell wall deposition, are organized by macromolecule interactions. Experimentally determined protein-protein interactions (PPIs) and protein-DNA interactions (PDIs) relevant to woody plant development are sparse. To begin to develop a Populus trichocarpa Torr. & A. Gray wood interactome, we applied the yeast-two-hybrid (Y2H) assay in different ways to enable the discovery of novel PPIs and connected networks. We first cloned open reading frames (ORFs) for 361 genes markedly upregulated in secondary xylem compared with secondary phloem and performed a binary Y2H screen with these proteins. By screening a xylem cDNA library for interactors of a subset of these proteins and then recapitulating the process by using a subset of the interactors as baits, we ultimately identified 165 PPIs involving 162 different ORFs. Thirty-eight transcription factors (TFs) included in our collection of P. trichocarpa wood ORFs were used in a Y1H screen for binding to promoter regions of three genes involved in lignin biosynthesis resulting in 40 PDIs involving 20 different TFs. The network incorporating both the PPIs and PDIs included 14 connected subnetworks, with the largest having 132 members. Protein-protein interactions and PDIs validated previous reports and also identified new candidate wood formation proteins and modules through their interactions with proteins and promoters known to be involved in secondary cell wall synthesis. Selected examples are discussed including a PPI between Mps one binder (MOB1) and a mitogen-activated protein kinase kinase kinase kinase (M4K) that was further characterized by assays confirming the PPI as well as its effect on subcellular localization. Mapping of published transcriptomic data showing developmentally detailed expression patterns across a secondary stem onto the network supported that the PPIs and PDIs are relevant to wood formation, and also illustrated that wood-associated interactions involve gene products that are not upregulated in secondary xylem.

Quantitative trait locus mapping of Populus bark features and stem diameter

BACKGROUND

Bark plays important roles in photosynthate transport and storage, along with physical and chemical protection. Bark texture varies extensively among species, from smooth to fissured to deeply furrowed, but its genetic control is unknown. This study sought to determine the main genomic regions associated with natural variation in bark features and stem diameter. Quantitative trait loci (QTL) were mapped using an interspecific pseudo-backcross pedigree (Populus trichocarpa x P. deltoides and P. deltoides) for bark texture, bark thickness and diameter collected across three years, two sites and three biological replicates per site.

RESULTS

QTL specific to bark texture were highly reproducible in shared intervals across sites, years and replicates. Significant positive correlations and co-localization between trait QTL suggest pleiotropic regulators or closely linked genes. A list of candidate genes with related putative function, location close to QTL maxima and with the highest expression level in the phloem, xylem and cambium was identified.

CONCLUSION

Candidate genes for bark texture included an ortholog of Arabidopsis ANAC104 (PopNAC128), which plays a role in lignified fiber cell and ray development, as well as Pinin and Fasciclin (PopFLA) genes with a role in cell adhesion, cell shape and migration. The results presented in this study provide a basis for future genomic characterization of genes found within the QTL for bark texture, bark thickness and diameter in order to better understand stem and bark development in Populus and other woody perennial plants. The QTL mapping approach identified a list of prime candidate genes for further validation using functional genomics or forward genetics approaches.

XYLEM NAC DOMAIN1, an angiosperm NAC transcription factor, inhibits xylem differentiation through conserved motifs that interact with RETINOBLASTOMA-RELATED

The Arabidopsis thaliana gene XYLEM NAC DOMAIN1 (XND1) is upregulated in xylem tracheary elements. Yet overexpression of XND1 blocks differentiation of tracheary elements. The molecular mechanism of XND1 action was investigated. Phylogenetic and motif analyses indicated that XND1 and its homologs are present only in angiosperms and possess a highly conserved C-terminal region containing linear motifs (CKII-acidic, LXCXE, E2F^TD -like and LXCXE-mimic) predicted to interact with the cell cycle and differentiation regulator RETINOBLASTOMA-RELATED (RBR). Protein-protein interaction and functional analyses of XND1 deletion mutants were used to test the importance of RBR-interaction motifs. Deletion of either the LXCXE or the LXCXE-mimic motif reduced both the XND1-RBR interaction and XND1 efficacy as a repressor of differentiation, with loss of the LXCXE motif having the strongest negative impacts. The function of the XND1 C-terminal domain could be partially replaced by RBR fused to the N-terminal domain of XND1. XND1 also transactivated gene expression in yeast and plants. The properties of XND1, a transactivator that depends on multiple linear RBR-interaction motifs to inhibit differentiation, have not previously been described for a plant protein. XND1 harbors an apparently angiosperm-specific combination of interaction motifs potentially linking the general differentiation regulator RBR with a xylem-specific pathway for inhibition of differentiation.

De novo transcriptome assembly and quantification reveal differentially expressed genes between soft-seed and hard-seed pomegranate (Punica granatum L.)

Pomegranate (Punica granatum L.) belongs to Punicaceae, and is valued for its social, ecological, economic, and aesthetic values, as well as more recently for its health benefits. The 'Tunisia' variety has softer seeds and big arils that are easily swallowed. It is a widely popular fruit; however, the molecular mechanisms of the formation of hard and soft seeds is not yet clear. We conducted a de novo assembly of the seed transcriptome in P. granatum L. and revealed differential gene expression between the soft-seed and hard-seed pomegranate varieties. A total of 35.1 Gb of data were acquired in this study, including 280,881,106 raw reads. Additionally, de novo transcriptome assembly generated 132,287 transcripts and 105,743 representative unigenes; approximately 13,805 unigenes (37.7%) were longer than 1,000 bp. Using bioinformatics annotation libraries, a total of 76,806 unigenes were annotated and, among the high-quality reads, 72.63% had at least one significant match to an existing gene model. Gene expression and differentially expressed genes were analyzed. The seed formation of the two pomegranate cultivars involves lignin biosynthesis and metabolism, including some genes encoding laccase and peroxidase, WRKY, MYB, and NAC transcription factors. In the hard-seed pomegranate, lignin-related genes and cellulose synthesis-related genes were highly expressed; in soft-seed pomegranates, expression of genes related to flavonoids and programmed cell death was slightly higher. We validated selection of the identified genes using qRT-PCR. This is the first transcriptome analysis of P. granatum L. This transcription sequencing greatly enriched the pomegranate molecular database, and the high-quality SSRs generated in this study will aid the gene cloning from pomegranate in the future. It provides important insights into the molecular mechanisms underlying the formation of soft seeds in pomegranate.

SSGA and MSGA: two seed-growing algorithms for constructing collaborative subnetworks

The establishment of a collaborative network of transcription factors (TFs) followed by decomposition and then construction of subnetworks is an effective way to obtain sets of collaborative TFs; each set controls a biological process or a complex trait. We previously developed eight gene association methods for genome-wide coexpression analysis between each TF and all other genomic genes and then constructing collaborative networks of TFs but only one algorithm, called Triple-Link Algorithm, for building collaborative subnetworks. In this study, we developed two more algorithms, Single Seed-Growing Algorithm (SSGA) and Multi-Seed Growing Algorithm (MSGA), for building collaborative subnetworks of TFs by identifying the fully-linked triple-node seeds from a decomposed collaborative network and then growing them into subnetworks with two different strategies. The subnetworks built from the three algorithms described above were comparatively appraised in terms of both functional cohesion and intra-subnetwork association strengths versus inter-subnetwork association strengths. We concluded that SSGA and MSGA, which performed more systemic comparisons and analyses of edge weights and network connectivity during subnetwork construction processes, yielded more functional and cohesive subnetworks than Triple-Link Algorithm. Together, these three algorithms provide alternate approaches for acquiring subnetworks of collaborative TFs. We also presented a framework to outline how to use these three algorithms to obtain collaborative TF sets governing biological processes or complex traits.

Poplar stem transcriptome is massively remodelled in response to single or repeated mechanical stimuli

Background

Trees experience mechanical stimuli -like wind- that trigger thigmomorphogenetic syndrome, leading to modifications of plant growth and wood quality. This syndrome affects tree productivity but is also believed to improve tree acclimation to chronic wind. Wind is particularly challenging for trees, because of their stature and perenniality. Climate change forecasts are predicting that the occurrence of high wind will worsen, making it increasingly vital to understand the mechanisms regulating thigmomorphogenesis, especially in perennial plants. By extension, this also implies factoring in the recurring nature of wind episodes. However, data on the molecular processes underpinning mechanoperception and transduction of mechanical signals, and their dynamics, are still dramatically lacking in trees.

Results

Here we performed a genome-wide and time-series analysis of poplar transcriptional responsiveness to transitory and recurring controlled stem bending, mimicking wind. The study revealed that 6% of the poplar genome is differentially expressed after a single transient bending. The combination of clustering, Gene Ontology categorization and time-series expression approaches revealed the diversity of gene expression patterns and biological processes affected by stem bending. Short-term transcriptomic responses entailed a rapid stimulation of plant defence and abiotic stress signalling pathways, including ethylene and jasmonic acid signalling but also photosynthesis process regulation. Late transcriptomic responses affected genes involved in cell wall organization and/or wood development. An analysis of the molecular impact of recurring bending found that the vast majority (96%) of the genes differentially expressed after a first bending presented reduced or even net-zero amplitude regulation after the second exposure to bending.

Conclusion

This study constitutes the first dynamic characterization of the molecular processes affected by single or repeated stem bending in poplar. Moreover, the global attenuation of the transcriptional responses, observed from as early as after a second bending, indicates the existence of a mechanism governing a fine tuning of plant responsiveness. This points toward several mechanistic pathways that can now be targeted to elucidate the complex dynamics of wind acclimation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3670-1) contains supplementary material, which is available to authorized users.

Functional Characterization of Populus PsnSHN2 in Coordinated Regulation of Secondary Wall Components in Tobacco

Wood formation is a biological process during which the most abundant lignocellulosic biomass on earth is produced. Although a number of transcription factors have been linked to the regulation of wood formation process, none of them has been demonstrated to be a higher hierarchical regulator that coordinately regulates secondary wall biosynthesis genes. Here, we identified a Populus gene, PsnSHN2, a counterpart of the Arabidopsis AP2/ERF type transcription factor, SHINE2. PsnSHN2 is predominantly expressed in xylem tissues and acted evidently as a high hierarchical transcriptional activator. Overexpression of PsnSHN2 in tobacco significantly altered the expression of both transcription factors and biosynthesis genes involved in secondary wall formation, leading to the thickened secondary walls and the changed cell wall composition. The most significant changes occurred in the contents of cellulose and hemicellulose that increased 37% and 28%, respectively, whereas the content of lignin that decreased 34%. Furthermore, PsnSHN2 activated or repressed the promoter activities of transcription factors involved in secondary wall biosynthesis and bound to five cis-acting elements enriched in the promoter regions of these transcription factors. Taken together, our results suggest PsnSHN2 coordinately regulate secondary wall formation through selective up/down-regulation of its downstream transcription factors that control secondary wall formation.

Characterization, Expression, and Functional Analysis of a Novel NAC Gene Associated with Resistance to Verticillium Wilt and Abiotic Stress in Cotton

Elucidating the mechanism of resistance to biotic and abiotic stress is of great importance in cotton. In this study, a gene containing the NAC domain, designated GbNAC1, was identified from Gossypium barbadense L. Homologous sequence alignment indicated that GbNAC1 belongs to the TERN subgroup. GbNAC1 protein localized to the cell nucleus. GbNAC1 was expressed in roots, stems, and leaves, and was especially highly expressed in vascular bundles. Functional analysis showed that cotton resistance to Verticillium wilt was reduced when the GbNAC1 gene was silenced using the virus-induced gene silencing (VIGS) method. GbNAC1-overexpressing Arabidopsis showed enhanced resistance to Verticillium dahliae compared to wild-type. Thus, GbNAC1 is involved in the positive regulation of resistance to Verticillium wilt. In addition, analysis of GbNAC1-overexpressing Arabidopsis under different stress treatments indicated that it is involved in plant growth, development, and response to various abiotic stresses (ABA, mannitol, and NaCl). This suggests that GbNAC1 plays an important role in resistance to biotic and abiotic stresses in cotton. This study provides a foundation for further study of the function of NAC genes in cotton and other plants.

Breeding maize for silage and biofuel production, an illustration of a step forward with the genome sequence

The knowledge of the gene families mostly impacting cell wall digestibility variations would significantly increase the efficiency of marker-assisted selection when breeding maize and grass varieties with improved silage feeding value and/or with better straw fermentability into alcohol or methane. The maize genome sequence of the B73 inbred line was released at the end of 2009, opening up new avenues to identify the genetic determinants of quantitative traits. Colocalizations between a large set of candidate genes putatively involved in secondary cell wall assembly and QTLs for cell wall digestibility (IVNDFD) were then investigated, considering physical positions of both genes and QTLs. Based on available data from six RIL progenies, 59 QTLs corresponding to 38 non-overlapping positions were matched up with a list of 442 genes distributed all over the genome. Altogether, 176 genes colocalized with IVNDFD QTLs and most often, several candidate genes colocalized at each QTL position. Frequent QTL colocalizations were found firstly with genes encoding ZmMYB and ZmNAC transcription factors, and secondly with genes encoding zinc finger, bHLH, and xylogen regulation factors. In contrast, close colocalizations were less frequent with genes involved in monolignol biosynthesis, and found only with the C4H2, CCoAOMT5, and CCR1 genes. Close colocalizations were also infrequent with genes involved in cell wall feruloylation and cross-linkages. Altogether, investigated colocalizations between candidate genes and cell wall digestibility QTLs suggested a prevalent role of regulation factors over constitutive cell wall genes on digestibility variations.

Wheat NAC transcription factor TaNAC29 is involved in response to salt stress

Soil salinity is considered as one of the most severe abiotic stress factors, which limit plant growth and cause significant losses in crop yield. NAC transcription factors have been proven to play vital roles in abiotic stress signaling in plants. As a staple crop, wheat production is severely constrained by salt stress whereas only a few NAC genes have been characterized functionally. To promote the application of NAC genes in wheat improvement by genetic engineering, a NAC gene designated TaNAC29 was characterized in common wheat. Expression analysis showed that TaNAC29 gene was involved in response to salt, drought and ABA treatments. TaNAC29 protein displays transactivation activity. To determine its role, transgenic Arabidopsis overexpressing TaNAC29 controlled by the CaMV-35S promoter was generated and subjected to salt stress for morphological and physiological assays. Morphological analysis showed that transgenic plants had enhanced tolerance to salt stress, as indicated by improved physiological traits, including more green leaves, reduced H2O2 accumulation, strengthened cell membrane stability and higher SOD, POD, CAT and APX activities. Moreover, the transcript levels of stress-related genes were significantly higher in TaNAC29 overexpression line than those in WT under salt treatment. Taken together, our results demonstrate that TaNAC29 confers salt stress tolerance through reducing H2O2 accumulation and membrane damage by enhancing the antioxidant system, and participating in regulating the abiotic stress-responsive signaling pathway.

NAC Transcription Factors in Senescence: From Molecular Structure to Function in Crops

Within the last decade, NAC transcription factors have been shown to play essential roles in senescence, which is the focus of this review. Transcriptome analyses associate approximately one third of Arabidopsis NAC genes and many crop NAC genes with senescence, thereby implicating NAC genes as important regulators of the senescence process. The consensus DNA binding site of the NAC domain is used to predict NAC target genes, and protein interaction sites can be predicted for the intrinsically disordered transcription regulatory domains of NAC proteins. The molecular characteristics of these domains determine the interactions in gene regulatory networks. Emerging local NAC-centered gene regulatory networks reveal complex molecular mechanisms of stress- and hormone-regulated senescence and basic physiological steps of the senescence process. For example, through molecular interactions involving the hormone abscisic acid, Arabidopsis NAP promotes chlorophyll degradation, a hallmark of senescence. Furthermore, studies of the functional rice ortholog, OsNAP, suggest that NAC genes can be targeted to obtain specific changes in lifespan control and nutrient remobilization in crop plants. This is also exemplified by the wheat NAM1 genes which promote senescence and increase grain zinc, iron, and protein content. Thus, NAC genes are promising targets for fine-tuning senescence for increased yield and quality.

A metabolomic assessment of NAC154 transcription factor overexpression in field grown poplar stem wood

Several xylem-associated regulatory genes have been identified that control processes associated with wood formation in poplar. Prominent among these are the NAC domain transcription factors (NACs). Here, the putative involvement of Populus NAC154, a co-ortholog of the Arabidopsis gene SND2, was evaluated as a regulator of "secondary" biosynthetic processes in stem internode tissues by interrogating aqueous methanolic extracts from control and transgenic trees. Comprehensive untargeted metabolite profiling was accomplished with a liquid chromatography-mass spectrometry platform that utilized two different chromatographic supports (HILIC and reversed phase) and both positive and negative ionization modes. Evaluation of current and previous year tissues provided datasets for assessing the effects of NAC154 overexpression in wood maturation processes. Phenolic glycoside levels as well as those of oligolignols, sucrose and arginine were modulated with phenotypic and chemotypic traits exhibiting similar trends. Specifically, increased levels of arginine in the NAC154 overexpressing tissues supports a role for the transcription factor in senescence/dormancy-associated processes.

Molecular control of wood formation in trees

Wood (also termed secondary xylem) is the most abundant biomass produced by plants, and is one of the most important sinks for atmospheric carbon dioxide. The development of wood begins with the differentiation of the lateral meristem, vascular cambium, into secondary xylem mother cells followed by cell expansion, secondary wall deposition, programmed cell death, and finally heartwood formation. Significant progress has been made in the past decade in uncovering the molecular players involved in various developmental stages of wood formation in tree species. Hormonal signalling has been shown to play critical roles in vascular cambium cell proliferation and a peptide-receptor-transcription factor regulatory mechanism similar to that controlling the activity of apical meristems is proposed to be involved in the maintenance of vascular cambium activity. It has been demonstrated that the differentiation of vascular cambium into xylem mother cells is regulated by plant hormones and HD-ZIP III transcription factors, and the coordinated activation of secondary wall biosynthesis genes during wood formation is mediated by a transcription network encompassing secondary wall NAC and MYB master switches and their downstream transcription factors. Most genes encoding the biosynthesis enzymes for wood components (cellulose, xylan, glucomannan, and lignin) have been identified in poplar and a number of them have been functionally characterized. With the availability of genome sequences of tree species from both gymnosperms and angiosperms, and the identification of a suite of wood-associated genes, it is expected that our understanding of the molecular control of wood formation in trees will be greatly accelerated.

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

NAC domain transcription factors regulate many developmental processes and stress responses in plants and vary widely in number and family structure. We analysed the characteristics and evolution of the NAC gene family of Eucalyptus grandis, a fast-growing forest tree in the rosid order Myrtales. NAC domain genes identified in the E. grandis genome were subjected to amino acid sequence, phylogenetic and motif analyses. Transcript abundance in developing tissues and abiotic stress conditions in E. grandis and E. globulus was quantified using RNA-seq and reverse transcription quantitative PCR (RT-qPCR). One hundred and eighty-nine E. grandis NAC (EgrNAC) proteins, arranged into 22 subfamilies, are extensively duplicated in subfamilies associated with stress response. Most EgrNAC genes form tandem duplicate arrays that frequently carry signatures of purifying selection. Sixteen amino acid motifs were identified in EgrNAC proteins, eight of which are enriched in, or unique to, Eucalyptus. New candidates for the regulation of normal and tension wood development and cold responses were identified. This first description of a Myrtales NAC domain family reveals an unique history of tandem duplication in stress-related subfamilies that has likely contributed to the adaptation of eucalypts to the challenging Australian environment. Several new candidates for the regulation of stress, wood formation and tree-specific development are reported.

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.

Genome-Wide Identification of VQ Motif-Containing Proteins and their Expression Profiles Under Abiotic Stresses in Maize

VQ motif-containing proteins play crucial roles in abiotic stress responses in plants. Recent studies have shown that some VQ proteins physically interact with WRKY transcription factors to activate downstream genes. In the present study, we identified and characterized genes encoding VQ motif-containing proteins using the most recent version of the maize genome sequence. In total, 61VQ genes were identified. In a cluster analysis, these genes clustered into nine groups together with their homologous genes in rice and Arabidopsis. Most of the VQ genes (57 out of 61 numbers) identified in maize were found to be single-copy genes. Analyses of RNA-seq data obtained using seedlings under long-term drought treatment showed that the expression levels of most ZmVQ genes (41 out of 61 members) changed during the drought stress response. Quantitative real-time PCR analyses showed that most of the ZmVQ genes were responsive to NaCl treatment. Also, approximately half of the ZmVQ genes were co-expressed with ZmWRKY genes. The identification of these VQ genes in the maize genome and knowledge of their expression profiles under drought and osmotic stresses will provide a solid foundation for exploring their specific functions in the abiotic stress responses of maize.

A role for OVATE FAMILY PROTEIN1 (OFP1) and OFP4 in a BLH6-KNAT7 multi-protein complex regulating secondary cell wall formation in Arabidopsis thaliana

Formation of secondary walls is a complex process that requires the coordinated and developmentally regulated expression of secondary wall biosynthetic genes. In Arabidopsis thaliana, a transcriptional network orchestrates the biosynthesis and deposition of the main SCW components in xylem and fiber cells. It was recently reported that interacting TALE homeodomain proteins BEL-LIKE HOMEODOMAIN6 (BLH6) and KNOTTED ARABIDOPSIS THALIANA7 (KNAT7) negatively regulate secondary cell wall formation in the interfascicular fibers of Arabidopsis inflorescence stems. Members of the Arabidopsis OVATE FAMILY PROTEIN (OFP) family of transcriptional regulators have been shown to physically interact in yeast with various KNAT and BLH proteins, forming a proposed TALE-OFP protein interaction network. This study presents molecular and genetic data indicating that OFP1 and OFP4, previously reported to interact with TALE homeodomain proteins, enhance the repression activity of BLH6, supporting a role for these OFPs as components of a putative multi-protein transcription regulatory complex containing BLH6 and KNAT7.

Functional repression of PtSND2 represses growth and development by disturbing auxin biosynthesis, transport and signaling in transgenic poplar

Using chimeric repressor silencing technology, we previously reported that functional repression of PtSND2 severely arrested wood formation in transgenic poplar (Populus). Here, we provide further evidence that auxin biosynthesis, transport and signaling were disturbed in these transgenic plants, leading to pleiotropic defects in their growth patterns, including inhibited leaf enlargement and vascular tissue development in the leaf central vein, suppressed cambial growth and fiber elongation in the stem, and arrested growth in the root system. Two transgenic lines, which displayed the most remarkable phenotypic deviation from the wild-type, were selected for detailed studies. In both transgenic lines, expression of genes for auxin biosynthesis, transport and signaling was down-regulated, and indole-3-acetic acid distribution was severely disturbed in the apical buds, leaves, stems and roots of field-grown transgenic plants. Transient transcription dual-luciferase assays of ProPtTYDC2::LUC, ProPttLAX2::LUC and ProPoptrIAA20.2::LUC in poplar protoplasts revealed that expression of auxin-related genes might be regulated by PtSND2 at the transcriptional level. All these results indicate that functional repression of PtSND2 altered auxin biosynthesis, transport and signaling, and thereby disturbed the normal growth and development of transgenic plants.

Genome-wide identification and expression pattern of drought-responsive members of the NAC family in maize

NAC proteins are plant-specific transcription factors (TFs). Although they play a pivotal role in regulating distinct biological processes, TFs in maize are yet to be investigated comprehensively. Within the maize genome, we identified 152 putative NAC domain-encoding genes (ZmNACs), including eight membrane-bound members, by systematic sequence analysis and physically mapped them onto ten chromosomes of maize. In silico analysis of the ZmNACs and comparison with similar genes in other plants such as Arabidopsis, rice, and soybean, revealed a similar NAC sequence architecture. Phylogenetically, the ZmNACs were arranged into six distinct subgroups (I-VI) possessing conserved motifs. Phylogenetic analysis using stress-related NAC TFs from Arabidopsis, rice, and soybean as seeding sequences identified 24 of the 152 ZmNACs (all from Group II) as putative stress-responsive genes, including one dehydration-responsive ZmSNAC1 gene reported earlier. One drought-tolerant genotype (HKI577) and one susceptible genotype (PC13T-3) were used for studying the expression pattern of the NAC genes during drought stress. qRT-PCR based expression profiles of 11 genes predicted to be related to stress confirmed strong differential gene expression during drought stress. Phylogenetic analyses revealed that ZmNAC18, ZmNAC51, ZmNAC145, and ZmNAC72, which were up-regulated in the tolerant genotype and down-regulated in the susceptible genotype, belonged to the same group to which also belong other drought-responsive genes, namely SNAC1, OsNAC6, ANAC019, and ANAC055, which act as a transcriptional activator and are strongly induced under stress from various abiotic sources. Differentially expressed ZmNAC genes, alone or in combination with each other or with other type(s) of TFs, may control the general cellular machinery and regulate stress-responsive downstream genes. Alternatively, they may serve as a platform to regulate a broad set of genes, which are subsequently fine-tuned by specific regulators. This genome-wide identification and expression profiling opens new avenues for systematic functional analysis of new members of the NAC gene family, which may be exploited in developing lines that are better adapted to drought.

Complexity of the transcriptional network controlling secondary wall biosynthesis

Secondary walls in the form of wood and fibers are the most abundant biomass produced by vascular plants, and are important raw materials for many industrial uses. Understanding how secondary walls are constructed is of significance in basic plant biology and also has far-reaching implications in genetic engineering of plant biomass better suited for various end uses, such as biofuel production. Secondary walls are composed of three major biopolymers, i.e., cellulose, hemicelluloses and lignin, the biosynthesis of which requires the coordinated transcriptional regulation of all their biosynthesis genes. Genomic and molecular studies have identified a number of transcription factors, whose expression is associated with secondary wall biosynthesis. We comprehensively review how these secondary wall-associated transcription factors function together to turn on the secondary wall biosynthetic program, which leads to secondary wall deposition in vascular plants. The transcriptional network regulating secondary wall biosynthesis employs a multi-leveled feed-forward loop regulatory structure, in which the top-level secondary wall NAC (NAM, ATAF1/2 and CUC2) master switches activate the second-level MYB master switches and they together induce the expression of downstream transcription factors and secondary wall biosynthesis genes. Secondary wall NAC master switches and secondary wall MYB master switches bind to and activate the SNBE (secondary wall NAC binding element) and SMRE (secondary wall MYB-responsive element) sites, respectively, in their target gene promoters. Further investigation of what and how developmental signals trigger the transcriptional network to regulate secondary wall biosynthesis and how different secondary wall-associated transcription factors function cooperatively in activating secondary wall biosynthetic pathways will lead to a better understanding of the molecular mechanisms underlying the transcriptional control of secondary wall biosynthesis.

A cotton (Gossypium hirsutum) gene encoding a NAC transcription factor is involved in negative regulation of plant xylem development

NAC proteins that compose of one large family of plant specific transcription factors (TF) play the important roles in many biological processes (such as morphogenesis, development, senescence and stress signal transduction). In this study, a gene (designated as GhXND1) encoding a NAC transcription factor was identified in cotton. Sequence analysis indicated that GhXND1 gene contains two introns inserted in its open reading frame (ORF). GhXND1 protein is localized in the cell nucleus, and displays the transactivation activity. GhXND1 transcripts were mainly detected in cotyledons, petals, roots, hypocotyls and stems, but little or no signals of GhXND1 expression were found in the other tissues. Ectopic expression of GhXND1 in Arabidopsis resulted in a reduction in number of xylem vessel cells and cell wall thickness of interfascicular fibers in the transgenic plants, compared with those of wild type. And expression of some cell wall biosynthesis-related genes was down-regulated in the GhXND1 transgenic plants. Collectively, the data presented in this study suggested that GhXND1 may be involved in regulation of plant xylem development.

Identification of predominant genes involved in regulation and execution of senescence-associated nitrogen remobilization in flag leaves of field grown barley

The transcriptomes of senescing flag leaves collected from barley field plots with standard or high nitrogen supply were compared to identify genes specifically associated with nitrogen remobilization during leaf senescence under agronomically relevant conditions. In flag leaves collected in field plots with high nitrogen supply, the decline in chlorophyll content was delayed. By comparing changes in gene expression for the two nitrogen levels, it was possible to discriminate genes related to nitrogen remobilization during senescence and genes involved in other processes associated with the late development of leaves under field conditions. Predominant genes that were more strongly upregulated during senescence of flag leaves from plants with standard nitrogen supply included genes encoding the transcription factor HvNAC026, serine type protease SCPL51, and the autophagy factors APG7 and ATG18F. Elevated expression of these genes in senescing leaves from plants with standard nitrogen supply indicates important roles of the corresponding proteins in nitrogen remobilization. In comparison, the genes upregulated in both flag leaf samples might have roles in general senescence processes associated with late leaf development. Among these genes were the transcription factor genes HvNAC001, HvNAC005, HvNAC013, HvWRKY12 and MYB, genes encoding the papain-like cysteine peptidases HvPAP14 and HvPAP20, as well as a subtilase gene.

Navigating the transcriptional roadmap regulating plant secondary cell wall deposition

The current status of lignocellulosic biomass as an invaluable resource in industry, agriculture, and health has spurred increased interest in understanding the transcriptional regulation of secondary cell wall (SCW) biosynthesis. The last decade of research has revealed an extensive network of NAC, MYB and other families of transcription factors regulating Arabidopsis SCW biosynthesis, and numerous studies have explored SCW-related transcription factors in other dicots and monocots. Whilst the general structure of the Arabidopsis network has been a topic of several reviews, they have not comprehensively represented the detailed protein-DNA and protein-protein interactions described in the literature, and an understanding of network dynamics and functionality has not yet been achieved for SCW formation. Furthermore the methodologies employed in studies of SCW transcriptional regulation have not received much attention, especially in the case of non-model organisms. In this review, we have reconstructed the most exhaustive literature-based network representations to date of SCW transcriptional regulation in Arabidopsis. We include a manipulable Cytoscape representation of the Arabidopsis SCW transcriptional network to aid in future studies, along with a list of supporting literature for each documented interaction. Amongst other topics, we discuss the various components of the network, its evolutionary conservation in plants, putative modules and dynamic mechanisms that may influence network function, and the approaches that have been employed in network inference. Future research should aim to better understand network function and its response to dynamic perturbations, whilst the development and application of genome-wide approaches such as ChIP-seq and systems genetics are in progress for the study of SCW transcriptional regulation in non-model organisms.

Chimeric repressor of PtSND2 severely affects wood formation in transgenic Populus

NAC domain transcription factors are important regulators that activate the secondary wall biosynthesis in wood formation. In this work, we investigated the possible functions of an NAC family member SECONDARY WALL-ASSOCIATED NAC DOMAIN PROTEIN2 (PtSND2) using chimeric repressor silencing technology. Reverse transcription-polymerase chain reaction, subcellular localization and transcriptional activation analyses indicated that PtSND2 is a wood-associated transcriptional factor with the predicted transcriptional activation activity, which could be inhibited by the repression domain SUPERMAN REPRESSION DOMAIN X (SRDX) in yeast. Wood formation was severely repressed in transgenic poplar plants overexpressing PtSND2-SRDX. Meanwhile, the secondary cell wall thickness of xylem fibers was restrained, and the contents of cellulose and lignin were obviously decreased in the stems of transgenic plants. Further studies indicated that expressions of a number of wood-associated genes were down-regulated in the stems of transgenic plants. Our results suggest that PtSND2 may play important roles during the secondary growth of stems in poplar.

Integrated -omics: a powerful approach to understanding the heterogeneous lignification of fibre crops

Lignin and cellulose represent the two main components of plant secondary walls and the most abundant polymers on Earth. Quantitatively one of the principal products of the phenylpropanoid pathway, lignin confers high mechanical strength and hydrophobicity to plant walls, thus enabling erect growth and high-pressure water transport in the vessels. Lignin is characterized by a high natural heterogeneity in its composition and abundance in plant secondary cell walls, even in the different tissues of the same plant. A typical example is the stem of fibre crops, which shows a lignified core enveloped by a cellulosic, lignin-poor cortex. Despite the great value of fibre crops for humanity, however, still little is known on the mechanisms controlling their cell wall biogenesis, and particularly, what regulates their spatially-defined lignification pattern. Given the chemical complexity and the heterogeneous composition of fibre crops' secondary walls, only the use of multidisciplinary approaches can convey an integrated picture and provide exhaustive information covering different levels of biological complexity. The present review highlights the importance of combining high throughput -omics approaches to get a complete understanding of the factors regulating the lignification heterogeneity typical of fibre crops.

Tissue-type-specific transcriptome analysis identifies developing xylem-specific promoters in poplar

Plant biotechnology offers a means to create novel phenotypes. However, commercial application of biotechnology in crop improvement programmes is severely hindered by the lack of utility promoters (or freedom to operate the existing ones) that can drive gene expression in a tissue-specific or temporally controlled manner. Woody biomass is gaining popularity as a source of fermentable sugars for liquid fuel production. To improve the quantity and quality of woody biomass, developing xylem (DX)-specific modification of the feedstock is highly desirable. To develop utility promoters that can drive transgene expression in a DX-specific manner, we used the Affymetrix Poplar Genome Arrays to obtain tissue-type-specific transcriptomes from poplar stems. Subsequent bioinformatics analysis identified 37 transcripts that are specifically or strongly expressed in DX cells of poplar. After further confirmation of their DX-specific expression using semi-quantitative PCR, we selected four genes (DX5, DX8, DX11 and DX15) for in vivo confirmation of their tissue-specific expression in transgenic poplars. The promoter regions of the selected DX genes were isolated and fused to a β-glucuronidase (GUS)-reported gene in a binary vector. This construct was used to produce transgenic poplars via Agrobacterium-mediated transformation. The GUS expression patterns of the resulting transgenic plants showed that these promoters were active in the xylem cells at early seedling growth and had strongest expression in the developing xylem cells at later growth stages of poplar. We conclude that these DX promoters can be used as a utility promoter for DX-specific biomass engineering.