Expansion and diversification of the Populus R2R3-MYB family of transcription factors

BAR expressolog identification: expression profile similarity ranking of homologous genes in plant species

Large numbers of sequences are now readily available for many plant species, allowing easy identification of homologous genes. However, orthologous gene identification across multiple species is made difficult by evolutionary events such as whole-genome or segmental duplications. Several developmental atlases of gene expression have been produced in the past couple of years, and it may be possible to use these transcript abundance data to refine ortholog predictions. In this study, clusters of homologous genes between seven plant species - Arabidopsis, soybean, Medicago truncatula, poplar, barley, maize and rice - were identified. Following this, a pipeline to rank homologs within gene clusters by both sequence and expression profile similarity was devised by determining equivalent tissues between species, with the best expression profile match being termed the 'expressolog'. Five electronic fluorescent pictograph (eFP) browsers were produced as part of this effort, to aid in visualization of gene expression data and to complement existing eFP browsers at the Bio-Array Resource (BAR). Within the eFP browser framework, these expression profile similarity rankings were incorporated into an Expressolog Tree Viewer to allow cross-species homolog browsing by both sequence and expression pattern similarity. Global analyses showed that orthologs with the highest sequence similarity do not necessarily exhibit the highest expression pattern similarity. Other orthologs may show different expression patterns, indicating that such genes may require re-annotation or more specific annotation. Ultimately, it is envisaged that this pipeline will aid in improvement of the functional annotation of genes and translational plant research.

Drought induces alterations in the stomatal development program in Populus

Much is known about the physiological control of stomatal aperture as a means by which plants adjust to water availability. By contrast, the role played by the modulation of stomatal development to limit water loss has received much less attention. The control of stomatal development in response to water deprivation in the genus Populus is explored here. Drought induced declines in stomatal conductance as well as an alteration in stomatal development in two genotypes of Populus balsamifera. Leaves that developed under water-deficit conditions had lower stomatal indices than leaves that developed under well-watered conditions. Transcript abundance of genes that could hypothetically underpin drought-responsive changes in stomatal development was examined, in two genotypes, across six time points, under two conditions, well-watered and with water deficit. Populus homologues of STOMAGEN, ERECTA (ER), STOMATA DENSITY AND DISTRIBUTION 1 (SDD1), and FAMA had variable transcript abundance patterns congruent with their role in the modulation of stomatal development in response to drought. Conversely, there was no significant variation in transcript abundance between genotypes or treatments for the Populus homologues of YODA (YDA) and TOO MANY MOUTHS (TMM). The findings highlight the role that could be played by stomatal development during leaf expansion as a longer term means by which to limit water loss from leaves. Moreover, the results point to the key roles played by the regulation of the homologues of STOMAGEN, ER, SDD1, and FAMA in the control of this response in poplar.

Genome-wide analysis of the MYB transcription factor superfamily in soybean

BACKGROUND

The MYB superfamily constitutes one of the most abundant groups of transcription factors described in plants. Nevertheless, their functions appear to be highly diverse and remain rather unclear. To date, no genome-wide characterization of this gene family has been conducted in a legume species. Here we report the first genome-wide analysis of the whole MYB superfamily in a legume species, soybean (Glycine max), including the gene structures, phylogeny, chromosome locations, conserved motifs, and expression patterns, as well as a comparative genomic analysis with Arabidopsis.

RESULTS

A total of 244 R2R3-MYB genes were identified and further classified into 48 subfamilies based on a phylogenetic comparative analysis with their putative orthologs, showed both gene loss and duplication events. The phylogenetic analysis showed that most characterized MYB genes with similar functions are clustered in the same subfamily, together with the identification of orthologs by synteny analysis, functional conservation among subgroups of MYB genes was strongly indicated. The phylogenetic relationships of each subgroup of MYB genes were well supported by the highly conserved intron/exon structures and motifs outside the MYB domain. Synonymous nucleotide substitution (dN/dS) analysis showed that the soybean MYB DNA-binding domain is under strong negative selection. The chromosome distribution pattern strongly indicated that genome-wide segmental and tandem duplication contribute to the expansion of soybean MYB genes. In addition, we found that ~ 4% of soybean R2R3-MYB genes had undergone alternative splicing events, producing a variety of transcripts from a single gene, which illustrated the extremely high complexity of transcriptome regulation. Comparative expression profile analysis of R2R3-MYB genes in soybean and Arabidopsis revealed that MYB genes play conserved and various roles in plants, which is indicative of a divergence in function.

CONCLUSIONS

In this study we identified the largest MYB gene family in plants known to date. Our findings indicate that members of this large gene family may be involved in different plant biological processes, some of which may be potentially involved in legume-specific nodulation. Our comparative genomics analysis provides a solid foundation for future functional dissection of this family gene.

Expression of the R2R3-MYB transcription factor TaMYB14 from Trifolium arvense activates proanthocyanidin biosynthesis in the legumes Trifolium repens and Medicago sativa

Proanthocyanidins (PAs) are oligomeric flavonoids and one group of end products of the phenylpropanoid pathway. PAs have been reported to be beneficial for human and animal health and are particularly important in pastoral agricultural systems for improved animal production and reduced greenhouse gas emissions. However, the main forage legumes grown in these systems, such as Trifolium repens and Medicago sativa, do not contain any substantial amounts of PAs in leaves. We have identified from the foliar PA-accumulating legume Trifolium arvense an R2R3-MYB transcription factor, TaMYB14, and provide evidence that this transcription factor is involved in the regulation of PA biosynthesis in legumes. TaMYB14 expression is necessary and sufficient to up-regulate late steps of the phenylpropanoid pathway and to induce PA biosynthesis. RNA interference silencing of TaMYB14 resulted in almost complete cessation of PA biosynthesis in T. arvense, whereas Nicotiana tabacum, M. sativa, and T. repens plants constitutively expressing TaMYB14 synthesized and accumulated PAs in leaves up to 1.8% dry matter. Targeted liquid chromatography-multistage tandem mass spectrometry analysis identified foliar PAs up to degree of polymerization 6 in leaf extracts. Hence, genetically modified M. sativa and T. repens plants expressing TaMYB14 provide a viable option for improving animal health and mitigating the negative environmental impacts of pastoral animal production systems.

Comprehensive analysis of CCCH zinc finger family in poplar (Populus trichocarpa)

BACKGROUND

CCCH zinc finger proteins contain a typical motif of three cysteines and one histidine residues and serve regulatory functions at all stages of mRNA metabolism. In plants, CCCH type zinc finger proteins comprise a large gene family represented by 68 members in Arabidopsis and 67 in rice. These CCCH proteins have been shown to play diverse roles in plant developmental processes and environmental responses. However, this family has not been studied in the model tree species Populus to date.

RESULTS

In the present study, a comprehensive analysis of the genes encoding CCCH zinc finger family in Populus was performed. Using a thorough annotation approach, a total of 91 full-length CCCH genes were identified in Populus, of which most contained more than one CCCH motif and a type of non-conventional C-X(11)-C-X(6)-C-X(3)-H motif was unique for Populus. All of the Populus CCCH genes were phylogeneticly clustered into 13 distinct subfamilies. In each subfamily, the gene structure and motif composition were relatively conserved. Chromosomal localization of these genes revealed that most of the CCCHs (81 of 90, 90 %) are physically distributed on the duplicated blocks. Thirty-four paralogous pairs were identified in Populus, of which 22 pairs (64.7 %) might be created by the whole genome segment duplication, whereas 4 pairs seem to be resulted from tandem duplications. In 91 CCCH proteins, we also identified 63 putative nucleon-cytoplasm shuttling proteins and 3 typical RNA-binding proteins. The expression profiles of all Populus CCCH genes have been digitally analyzed in six tissues across different developmental stages, and under various drought stress conditions. A variety of expression patterns of CCCH genes were observed during Populus development, of which 34 genes highly express in root and 22 genes show the highest level of transcript abundance in differentiating xylem. Quantitative real-time RT-PCR (RT-qPCR) was further performed to confirm the tissue-specific expression and responses to drought stress treatment of 12 selected Populus CCCH genes.

CONCLUSIONS

This study provides the first systematic analysis of the Populus CCCH proteins. Comprehensive genomic analyses suggested that segmental duplications contribute significantly to the expansion of Populus CCCH gene family. Transcriptome profiling provides first insights into the functional divergences among members of Populus CCCH gene family. Particularly, some CCCH genes may be involved in wood development while others in drought tolerance regulation. Our results presented here may provide a starting point for the functional dissection of this family of potential RNA-binding proteins.

Apple miRNAs and tasiRNAs with novel regulatory networks

BACKGROUND

MicroRNAs (miRNAs) and their regulatory functions have been extensively characterized in model species but whether apple has evolved similar or unique regulatory features remains unknown.

RESULTS

We performed deep small RNA-seq and identified 23 conserved, 10 less-conserved and 42 apple-specific miRNAs or families with distinct expression patterns. The identified miRNAs target 118 genes representing a wide range of enzymatic and regulatory activities. Apple also conserves two TAS gene families with similar but unique trans-acting small interfering RNA (tasiRNA) biogenesis profiles and target specificities. Importantly, we found that miR159, miR828 and miR858 can collectively target up to 81 MYB genes potentially involved in diverse aspects of plant growth and development. These miRNA target sites are differentially conserved among MYBs, which is largely influenced by the location and conservation of the encoded amino acid residues in MYB factors. Finally, we found that 10 of the 19 miR828-targeted MYBs undergo small interfering RNA (siRNA) biogenesis at the 3' cleaved, highly divergent transcript regions, generating over 100 sequence-distinct siRNAs that potentially target over 70 diverse genes as confirmed by degradome analysis.

CONCLUSIONS

Our work identified and characterized apple miRNAs, their expression patterns, targets and regulatory functions. We also discovered that three miRNAs and the ensuing siRNAs exploit both conserved and divergent sequence features of MYB genes to initiate distinct regulatory networks targeting a multitude of genes inside and outside the MYB family.

The R2R3-MYB transcription factor gene family in maize

MYB proteins comprise a large family of plant transcription factors, members of which perform a variety of functions in plant biological processes. To date, no genome-wide characterization of this gene family has been conducted in maize (Zea mays). In the present study, we performed a comprehensive computational analysis, to yield a complete overview of the R2R3-MYB gene family in maize, including the phylogeny, expression patterns, and also its structural and functional characteristics. The MYB gene structure in maize and Arabidopsis were highly conserved, indicating that they were originally compact in size. Subgroup-specific conserved motifs outside the MYB domain may reflect functional conservation. The genome distribution strongly supports the hypothesis that segmental and tandem duplication contribute to the expansion of maize MYB genes. We also performed an updated and comprehensive classification of the R2R3-MYB gene families in maize and other plant species. The result revealed that the functions were conserved between maize MYB genes and their putative orthologs, demonstrating the origin and evolutionary diversification of plant MYB genes. Species-specific groups/subgroups may evolve or be lost during evolution, resulting in functional divergence. Expression profile study indicated that maize R2R3-MYB genes exhibit a variety of expression patterns, suggesting diverse functions. Furthermore, computational prediction potential targets of maize microRNAs (miRNAs) revealed that miR159, miR319, and miR160 may be implicated in regulating maize R2R3-MYB genes, suggesting roles of these miRNAs in post-transcriptional regulation and transcription networks. Our comparative analysis of R2R3-MYB genes in maize confirm and extend the sequence and functional characteristics of this gene family, and will facilitate future functional analysis of the MYB gene family in maize.

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.

Transcription factors of Lotus: regulation of isoflavonoid biosynthesis requires coordinated changes in transcription factor activity

Isoflavonoids are a class of phenylpropanoids made by legumes, and consumption of dietary isoflavonoids confers benefits to human health. Our aim is to understand the regulation of isoflavonoid biosynthesis. Many studies have shown the importance of transcription factors in regulating the transcription of one or more genes encoding enzymes in phenylpropanoid metabolism. In this study, we coupled bioinformatics and coexpression analysis to identify candidate genes encoding transcription factors involved in regulating isoflavonoid biosynthesis in Lotus (Lotus japonicus). Genes encoding proteins belonging to 39 of the main transcription factor families were examined by microarray analysis of RNA from leaf tissue that had been elicited with glutathione. Phylogenetic analyses of each transcription factor family were used to identify subgroups of proteins that were specific to L. japonicus or closely related to known regulators of the phenylpropanoid pathway in other species. R2R3MYB subgroup 2 genes showed increased expression after treatment with glutathione. One member of this subgroup, LjMYB14, was constitutively overexpressed in L. japonicus and induced the expression of at least 12 genes that encoded enzymes in the general phenylpropanoid and isoflavonoid pathways. A distinct set of six R2R3MYB subgroup 2-like genes was identified. We suggest that these subgroup 2 sister group proteins and those belonging to the main subgroup 2 have roles in inducing isoflavonoid biosynthesis. The induction of isoflavonoid production in L. japonicus also involves the coordinated down-regulation of competing biosynthetic pathways by changing the expression of other transcription factors.

Expansion and diversification of the SET domain gene family following whole-genome duplications in Populus trichocarpa

Background

Histone lysine methylation modifies chromatin structure and regulates eukaryotic gene transcription and a variety of developmental and physiological processes. SET domain proteins are lysine methyltransferases containing the evolutionarily-conserved SET domain, which is known to be the catalytic domain.

Results

We identified 59 SET genes in the Populus genome. Phylogenetic analyses of 106 SET genes from Populus and Arabidopsis supported the clustering of SET genes into six distinct subfamilies and identified 19 duplicated gene pairs in Populus. The chromosome locations of these gene pairs and the distribution of synonymous substitution rates showed that the expansion of the SET gene family might be caused by large-scale duplications in Populus. Comparison of gene structures and domain architectures of each duplicate pair indicated that divergence took place at the 3'- and 5'-terminal transcribed regions and at the N- and C-termini of the predicted proteins, respectively. Expression profile analysis of Populus SET genes suggested that most Populus SET genes were expressed widely, many with the highest expression in young leaves. In particular, the expression profiles of 12 of the 19 duplicated gene pairs fell into two types of expression patterns.

Conclusions

The 19 duplicated SET genes could have originated from whole genome duplication events. The differences in SET gene structure, domain architecture, and expression profiles in various tissues of Populus suggest that members of the SET gene family have a variety of developmental and physiological functions. Our study provides clues about the evolution of epigenetic regulation of chromatin structure and gene expression.

The Class II KNOX gene KNAT7 negatively regulates secondary wall formation in Arabidopsis and is functionally conserved in Populus

• The formation of secondary cell walls in cell types such as tracheary elements and fibers is a defining characteristic of vascular plants. The Arabidopsis transcription factor KNAT7 is a component of a transcription network that regulates secondary cell wall biosynthesis, but its function has remained unclear. • We conducted anatomical, biochemical and molecular phenotypic analyses of Arabidopsis knat7 loss-of-function alleles, KNAT7 over-expression lines and knat7 lines expressing poplar KNAT7. • KNAT7 was strongly expressed in concert with secondary wall formation in Arabidopsis and poplar. Arabidopsis knat7 loss-of-function alleles exhibited irregular xylem phenotypes, but also showed increased secondary cell wall thickness in fibers. Increased commitment to secondary cell wall biosynthesis was accompanied by increased lignin content and elevated expression of secondary cell wall biosynthetic genes. KNAT7 over-expression resulted in thinner interfascicular fiber cell walls. • Taken together with data demonstrating that KNAT7 is a transcriptional repressor, we hypothesize that KNAT7 is a negative regulator of secondary wall biosynthesis, and functions in a negative feedback loop that represses metabolically inappropriate commitment to secondary wall formation, thereby maintaining metabolic homeostasis. The conservation of the KNAT7 regulatory module in poplar suggests new ways to manipulate secondary cell wall deposition for improvement of bioenergy traits in this tree.

Genome-wide identification, evolutionary expansion, and expression profile of homeodomain-leucine zipper gene family in poplar (Populus trichocarpa)

Background

Homeodomain-leucine zipper (HD-ZIP) proteins are plant-specific transcriptional factors known to play crucial roles in plant development. Although sequence phylogeny analysis of Populus HD-ZIPs was carried out in a previous study, no systematic analysis incorporating genome organization, gene structure, and expression compendium has been conducted in model tree species Populus thus far.

Principal Findings

In this study, a comprehensive analysis of Populus HD-ZIP gene family was performed. Sixty-three full-length HD-ZIP genes were found in Populus genome. These Populus HD-ZIP genes were phylogenetically clustered into four distinct subfamilies (HD-ZIP I–IV) and predominately distributed across 17 linkage groups (LG). Fifty genes from 25 Populus paralogous pairs were located in the duplicated blocks of Populus genome and then preferentially retained during the sequential evolutionary courses. Genomic organization analyses indicated that purifying selection has played a pivotal role in the retention and maintenance of Populus HD-ZIP gene family. Microarray analysis has shown that 21 Populus paralogous pairs have been differentially expressed across different tissues and under various stresses, with five paralogous pairs showing nearly identical expression patterns, 13 paralogous pairs being partially redundant and three paralogous pairs diversifying significantly. Quantitative real-time RT-PCR (qRT-PCR) analysis performed on 16 selected Populus HD-ZIP genes in different tissues and under both drought and salinity stresses confirms their tissue-specific and stress-inducible expression patterns.

Conclusions

Genomic organizations indicated that segmental duplications contributed significantly to the expansion of Populus HD-ZIP gene family. Exon/intron organization and conserved motif composition of Populus HD-ZIPs are highly conservative in the same subfamily, suggesting the members in the same subfamilies may also have conservative functionalities. Microarray and qRT-PCR analyses showed that 89% (56 out of 63) of Populus HD-ZIPs were duplicate genes that might have been retained by substantial subfunctionalization. Taken together, these observations may lay the foundation for future functional analysis of Populus HD-ZIP genes to unravel their biological roles.

Phylogenetic analyses unravel the evolutionary history of NAC proteins in plants

NAC (NAM/ATAF/CUC) proteins are one of the largest groups of transcription factors in plants. Although many NAC proteins based on Arabidopsis and rice genomes have been reported in a number of species, a complete survey and classification of all NAC genes in plant species from disparate evolutionary groups is lacking. In this study, we analyzed whole-genome sequences from nine major lineages of land plants to unveil the relationships between these proteins. Our results show that there are fewer than 30 NAC proteins present in both mosses and lycophytes, whereas more than 100 were found in most of the angiosperms. Phylogenetic analyses suggest that NAC proteins consist of 21 subfamilies, most of which have highly conserved non-NAC domain motifs. Six of these subfamilies existed in early-diverged land plants, whereas the remainder diverged only within the angiosperms. We hypothesize that NAC proteins probably originated sometime more than 400 million years ago and expanded together with the differentiation of plants into organisms of increasing complexity possibly after the divergence of lycophytes from the other vascular plants.

Acetylesterase-mediated deacetylation of pectin impairs cell elongation, pollen germination, and plant reproduction

Pectin is a major component of the primary cell wall of higher plants. Some galacturonyl residues in the backbone of pectinaceous polysaccharides are often O-acetylated at the C-2 or C-3 position, and the resulting acetylesters change dynamically during the growth and development of plants. The processes involve both enzymatic acetylation and deacetylation. Through genomic sequence analysis, we identified a pectin acetylesterase (PAE1) from black cottonwood (Populus trichocarpa). Recombinant Pt PAE1 exhibited preferential activity in releasing the acetate moiety from sugar beet (Beta vulgaris) and potato (Solanum tuberosum) pectin in vitro. Overexpressing Pt PAE1 in tobacco (Nicotiana tabacum) decreased the level of acetyl esters of pectin but not of xylan. Deacetylation engendered differential changes in the composition and/or structure of cell wall polysaccharides that subsequently impaired the cellular elongation of floral styles and filaments, the germination of pollen grains, and the growth of pollen tubes. Consequently, plants overexpressing PAE1 exhibited severe male sterility. Furthermore, in contrast to the conventional view, PAE1-mediated deacetylation substantially lowered the digestibility of pectin. Our data suggest that pectin acetylesterase functions as an important structural regulator in planta by modulating the precise status of pectin acetylation to affect the remodeling and physiochemical properties of the cell wall's polysaccharides, thereby affecting cell extensibility.

The R2R3-MYB transcription factor gene family in maize

MYB proteins comprise a large family of plant transcription factors, members of which perform a variety of functions in plant biological processes. To date, no genome-wide characterization of this gene family has been conducted in maize (Zea mays). In the present study, we performed a comprehensive computational analysis, to yield a complete overview of the R2R3-MYB gene family in maize, including the phylogeny, expression patterns, and also its structural and functional characteristics. The MYB gene structure in maize and Arabidopsis were highly conserved, indicating that they were originally compact in size. Subgroup-specific conserved motifs outside the MYB domain may reflect functional conservation. The genome distribution strongly supports the hypothesis that segmental and tandem duplication contribute to the expansion of maize MYB genes. We also performed an updated and comprehensive classification of the R2R3-MYB gene families in maize and other plant species. The result revealed that the functions were conserved between maize MYB genes and their putative orthologs, demonstrating the origin and evolutionary diversification of plant MYB genes. Species-specific groups/subgroups may evolve or be lost during evolution, resulting in functional divergence. Expression profile study indicated that maize R2R3-MYB genes exhibit a variety of expression patterns, suggesting diverse functions. Furthermore, computational prediction potential targets of maize microRNAs (miRNAs) revealed that miR159, miR319, and miR160 may be implicated in regulating maize R2R3-MYB genes, suggesting roles of these miRNAs in post-transcriptional regulation and transcription networks. Our comparative analysis of R2R3-MYB genes in maize confirm and extend the sequence and functional characteristics of this gene family, and will facilitate future functional analysis of the MYB gene family in maize.

Molecular characterization of 60 isolated wheat MYB genes and analysis of their expression during abiotic stress

The proteins of the MYB superfamily play central roles in developmental processes and defence responses in plants. Sixty unique wheat MYB genes that contain full-length cDNA sequences were isolated. These 60 genes were grouped into three categories, namely one R1R2R3-MYB, 22 R2R3-MYBs, and 37 MYB-related members. The sequence composition of the R2 and R3 repeats was conserved among the 22 wheat R2R3-MYB proteins. Phylogenetic comparison of the members of this superfamily among wheat, rice, and Arabidopsis revealed that the putative functions of some wheat MYB proteins were clustered into the Arabidopsis functional clades. Tissue-specific expression profiles showed that most of the wheat MYB genes were expressed in all of the tissues examined, suggesting that wheat MYB genes take part in multiple cellular processes. The expression analysis during abiotic stress identified a group of MYB genes that respond to one or more stress treatments. The overexpression of a salt-inducible gene, TaMYB32, enhanced the tolerance to salt stress in transgenic Arabidopsis. This study is the first comprehensive study of the MYB gene family in Triticeae.

Molecular features of secondary vascular tissue regeneration after bark girdling in Populus

Regeneration is a common strategy for plants to repair damage to their tissue after attacks from other organisms or physical assaults. However, how differentiating cells acquire regenerative competence and rebuild the pattern of new tissues remains largely unknown. Using anatomical observation and microarray analysis, we investigated the morphological process and molecular features of secondary vascular tissue regeneration after bark girdling in trees. After bark girdling, new phloem and cambium regenerate from differentiating xylem cells and rebuild secondary vascular tissue pattern within 1 month. Differentiating xylem cells acquire regenerative competence through epigenetic regulation and cell cycle re-entry. The xylem developmental program was blocked, whereas the phloem or cambium program was activated, resulting in the secondary vascular tissue pattern re-establishment. Phytohormones play important roles in vascular tissue regeneration. We propose a model describing the molecular features of secondary vascular tissue regeneration after bark girdling in trees. It provides information for understanding mechanisms of tissue regeneration and pattern formation of the secondary vascular tissues in plants.

Comparative analysis of GT14/GT14-like gene family in Arabidopsis, Oryza, Populus, Sorghum and Vitis

Glycosyltransferase family14 (GT14) belongs to the glycosyltransferase (GT) superfamily that plays important roles in the biosynthesis of cell walls, the most abundant source of cellulosic biomass for bioethanol production. It has been hypothesized that DUF266 proteins are a new class of GTs related to GT14. In this study, we identified 62 GT14 and 106 DUF266 genes (named GT14-like herein) in Arabidopsis, Oryza, Populus, Sorghum and Vitis. Our phylogenetic analysis separated GT14 and GT14-like genes into two distinct clades, which were further divided into eight and five groups, respectively. Similarities in protein domain, 3D structure and gene expression were uncovered between the two phylogenetic clades, supporting the hypothesis that GT14 and GT14-like genes belong to one family. Therefore, we proposed a new family name, GT14/GT14-like family that combines both subfamilies. Variation in gene expression and protein subcellular localization within the GT14-like subfamily were greater than those within the GT14 subfamily. One-half of the Arabidopsis and Populus GT14/GT14-like genes were found to be preferentially expressed in stem/xylem, indicating that they are likely involved in cell wall biosynthesis. This study provided new insights into the evolution and functional diversification of the GT14/GT14-like family genes.

Dissection of the transcriptional program regulating secondary wall biosynthesis during wood formation in poplar

Wood biomass is mainly made of secondary cell walls; hence, elucidation of the molecular mechanisms underlying the transcriptional regulation of secondary wall biosynthesis during wood formation will be instrumental to design strategies for genetic improvement of wood biomass. Here, we provide direct evidence demonstrating that the poplar (Populus trichocarpa) wood-associated NAC domain transcription factors (PtrWNDs) are master switches activating a suite of downstream transcription factors, and together, they are involved in the coordinated regulation of secondary wall biosynthesis during wood formation. We show that transgenic poplar plants with dominant repression of PtrWNDs functions exhibit a drastic reduction in secondary wall thickening in woody cells, and those with PtrWND overexpression result in ectopic deposition of secondary walls. Analysis of PtrWND2B overexpressors revealed up-regulation of the expression of a number of wood-associated transcription factors, the promoters of which were also activated by PtrWND6B and the Eucalyptus EgWND1. Transactivation analysis and electrophoretic mobility shift assay demonstrated that PtrWNDs and EgWND1 activated gene expression through direct binding to the secondary wall NAC-binding elements, which are present in the promoters of several wood-associated transcription factors and a number of genes involved in secondary wall biosynthesis and modification. The WND-regulated transcription factors PtrNAC150, PtrNAC156, PtrNAC157, PtrMYB18, PtrMYB74, PtrMYB75, PtrMYB121, PtrMYB128, PtrZF1, and PtrGATA8 were able to activate the promoter activities of the biosynthetic genes for all three major wood components. Our study has uncovered that the WND master switches together with a battery of their downstream transcription factors form a transcriptional network controlling secondary wall biosynthesis during wood formation.

The poplar-poplar rust interaction: insights from genomics and transcriptomics

Poplars are extensively cultivated worldwide, and their susceptibility to the leaf rust fungus Melampsora larici-populina leads to considerable damages in plantations. Despite a good knowledge of the poplar rust life cycle, and particularly the epidemics on poplar, the perennial status of the plant host and the obligate biotrophic lifestyle of the rust fungus are bottlenecks for molecular investigations. Following the completion of both M. larici-populina and Populus trichocarpa genome sequences, gene families involved in poplar resistance or in rust fungus virulence were investigated, allowing the identification of key genetic determinants likely controlling the outcome of the interaction. Specific expansions of resistance and defense-related genes in poplar indicate probable innovations in perennial species in relation with host-pathogen interactions. The genome of M. Larici-populina contains a strikingly high number of genes encoding small secreted proteins (SSPs) representing hundreds of candidate effectors. Transcriptome analyses of interacting partners in compatible and incompatible interactions revealed conserved set of genes involved in poplar defense reactions as well as timely regulated expression of SSP transcripts during host tissues colonisation. Ongoing functional studies of selected candidate effectors will be achieved mainly on the basis of recombinant protein purification and subsequent characterisation.

Shoot branching and leaf dissection in tomato are regulated by homologous gene modules

Aerial plant architecture is predominantly determined by shoot branching and leaf morphology, which are governed by apparently unrelated developmental processes, axillary meristem formation, and leaf dissection. Here, we show that in tomato (Solanum lycopersicum), these processes share essential functions in boundary establishment. Potato leaf (C), a key regulator of leaf dissection, was identified to be the closest paralog of the shoot branching regulator Blind (Bl). Comparative genomics revealed that these two R2R3 MYB genes are orthologs of the Arabidopsis thaliana branching regulator REGULATOR OF AXILLARY MERISTEMS1 (RAX1). Expression studies and complementation analyses indicate that these genes have undergone sub- or neofunctionalization due to promoter differentiation. C acts in a pathway independent of other identified leaf dissection regulators. Furthermore, the known leaf complexity regulator Goblet (Gob) is crucial for axillary meristem initiation and acts in parallel to C and Bl. Finally, RNA in situ hybridization revealed that the branching regulator Lateral suppressor (Ls) is also expressed in leaves. All four boundary genes, C, Bl, Gob, and Ls, may act by suppressing growth, as indicated by gain-of-function plants. Thus, leaf architecture and shoot architecture rely on a conserved mechanism of boundary formation preceding the initiation of leaflets and axillary meristems.

A NAC domain protein family contributing to the regulation of wood formation in poplar

Wood harvested from trees is one of the most widely utilized natural materials on our planet. Recent environmental issues have prompted an increase in the demand for wood, especially as a cost-effective and renewable resource for industry and energy, so it is important to understand the process of wood formation. In the present study, we focused on poplar (Populus trichocarpa) NAC domain protein genes which are homologous to well-known Arabidopsis transcription factors regulating the differentiation of xylem vessels and fiber cells. From phylogenetic analysis, we isolated 16 poplar NAC domain protein genes, and named them PtVNS (VND-, NST/SND- and SMB-related proteins) genes. Expression analysis revealed that 12 PtVNS (also called PtrWND) genes including both VND and NST groups were expressed in developing xylem tissue and phloem fiber, whereas in primary xylem vessels, only PtVNS/PtrWND genes of the VND group were expressed. By using the post-translational induction system of Arabidopsis VND7, a master regulator of xylem vessel element differentiation, many poplar genes functioning in xylem vessel differentiation downstream from NAC domain protein genes were identified. Transient expression assays showed the variation in PtVNS/PtrWND transactivation activity toward downstream genes, even between duplicate gene pairs. Furthermore, overexpression of PtVNS/PtrWND genes induced ectopic secondary wall thickening in poplar leaves as well as in Arabidopsis seedlings with different levels of induction efficiency according to the gene. These results suggest that wood formation in poplar is regulated by cooperative functions of the NAC domain proteins.

Gene coexpression network analysis as a source of functional annotation for rice genes

With the existence of large publicly available plant gene expression data sets, many groups have undertaken data analyses to construct gene coexpression networks and functionally annotate genes. Often, a large compendium of unrelated or condition-independent expression data is used to construct gene networks. Condition-dependent expression experiments consisting of well-defined conditions/treatments have also been used to create coexpression networks to help examine particular biological processes. Gene networks derived from either condition-dependent or condition-independent data can be difficult to interpret if a large number of genes and connections are present. However, algorithms exist to identify modules of highly connected and biologically relevant genes within coexpression networks. In this study, we have used publicly available rice (Oryza sativa) gene expression data to create gene coexpression networks using both condition-dependent and condition-independent data and have identified gene modules within these networks using the Weighted Gene Coexpression Network Analysis method. We compared the number of genes assigned to modules and the biological interpretability of gene coexpression modules to assess the utility of condition-dependent and condition-independent gene coexpression networks. For the purpose of providing functional annotation to rice genes, we found that gene modules identified by coexpression analysis of condition-dependent gene expression experiments to be more useful than gene modules identified by analysis of a condition-independent data set. We have incorporated our results into the MSU Rice Genome Annotation Project database as additional expression-based annotation for 13,537 genes, 2,980 of which lack a functional annotation description. These results provide two new types of functional annotation for our database. Genes in modules are now associated with groups of genes that constitute a collective functional annotation of those modules. Additionally, the expression patterns of genes across the treatments/conditions of an expression experiment comprise a second form of useful annotation.

Molecular dissection of xylan biosynthesis during wood formation in poplar

Xylan, being the second most abundant polysaccharide in dicot wood, is considered to be one of the factors contributing to wood biomass recalcitrance for biofuel production. To better utilize wood as biofuel feedstock, it is crucial to functionally characterize all the genes involved in xylan biosynthesis during wood formation. In this report, we investigated roles of poplar families GT43 and GT8 glycosyltransferases in xylan biosynthesis during wood formation. There exist seven GT43 genes in the genome of poplar (Populus trichocarpa), five of which, namely PtrGT43A, PtrGT43B, PtrGT43C, PtrGT43D, and PtrGT43E, were shown to be highly expressed in the developing wood and their encoded proteins were localized in the Golgi. Comprehensive genetic complementation coupled with chemical analyses demonstrated that overexpression of PtrGT43A/B/E but not PtrGT43C/D was able to rescue the xylan defects conferred by the Arabidopsis irx9 mutant, whereas overexpression of PtrGT43C/D but not PtrGT43A/B/E led to a complementation of the xylan defects in the Arabidopsis irx14 mutant. The essential roles of poplar GT43 members in xylan biosynthesis was further substantiated by RNAi down-regulation of GT43B in the hybrid poplar (Populus alba x tremula) leading to reductions in wall thickness and xylan content in wood, and an elevation in the abundance of the xylan reducing end sequence. Wood digestibility analysis revealed that cellulase digestion released more glucose from the wood of poplar GT43B RNAi lines than the control wood, indicating a decrease in wood biomass recalcitrance. Furthermore, RNAi down-regulation of another poplar wood-associated glycosyltransferase, PoGT8D, was shown to cause decreases in wall thickness and xylan content as well as in the abundance of the xylan reducing end sequence. Together, these findings demonstrate that the poplar GT43 members form two functionally non-redundant groups, namely PtrGT43A/B/E as functional orthologs of Arabidopsis IRX9 and PtrGT43C/D as functional orthologs of Arabidopsis IRX14, all of which are involved in the biosynthesis of xylan backbones, and that the poplar GT8D is essential for the biosynthesis of the xylan reducing end sequence.

Integration of molecular biology tools for identifying promoters and genes abundantly expressed in flowers of Oncidium Gower Ramsey

BACKGROUND

Orchids comprise one of the largest families of flowering plants and generate commercially important flowers. However, model plants, such as Arabidopsis thaliana do not contain all plant genes, and agronomic and horticulturally important genera and species must be individually studied.

RESULTS

Several molecular biology tools were used to isolate flower-specific gene promoters from Oncidium 'Gower Ramsey' (Onc. GR). A cDNA library of reproductive tissues was used to construct a microarray in order to compare gene expression in flowers and leaves. Five genes were highly expressed in flower tissues, and the subcellular locations of the corresponding proteins were identified using lip transient transformation with fluorescent protein-fusion constructs. BAC clones of the 5 genes, together with 7 previously published flower- and reproductive growth-specific genes in Onc. GR, were identified for cloning of their promoter regions. Interestingly, 3 of the 5 novel flower-abundant genes were putative trypsin inhibitor (TI) genes (OnTI1, OnTI2 and OnTI3), which were tandemly duplicated in the same BAC clone. Their promoters were identified using transient GUS reporter gene transformation and stable A. thaliana transformation analyses.

CONCLUSIONS

By combining cDNA microarray, BAC library, and bombardment assay techniques, we successfully identified flower-directed orchid genes and promoters.

Evolutionary and comparative analysis of MYB and bHLH plant transcription factors

The expansion of gene families encoding regulatory proteins is typically associated with the increase in complexity characteristic of multi-cellular organisms. The MYB and basic helix-loop-helix (bHLH) families provide excellent examples of how gene duplication and divergence within particular groups of transcription factors are associated with, if not driven by, the morphological and metabolic diversity that characterize the higher plants. These gene families expanded dramatically in higher plants; for example, there are approximately 339 and 162 MYB and bHLH genes, respectively, in Arabidopsis, and approximately 230 and 111, respectively, in rice. In contrast, the Chlamydomonas genome has only 38 MYB genes and eight bHLH genes. In this review, we compare the MYB and bHLH gene families from structural, evolutionary and functional perspectives. The knowledge acquired on the role of many of these factors in Arabidopsis provides an excellent reference to explore sequence-function relationships in crops and other plants. The physical interaction and regulatory synergy between particular sub-classes of MYB and bHLH factors is perhaps one of the best examples of combinatorial plant gene regulation. However, members of the MYB and bHLH families also interact with a number of other regulatory proteins, forming complexes that either activate or repress the expression of sets of target genes that are increasingly being identified through a diversity of high-throughput genomic approaches. The next few years are likely to witness an increasing understanding of the extent to which conserved transcription factors participate at similar positions in gene regulatory networks across plant species.

Sequence and function of basic helix-loop-helix proteins required for stomatal development in Arabidopsis are deeply conserved in land plants

Stomata are a broadly conserved feature of land plants with a crucial role regulating transpiration and gas exchange between the plant and atmosphere. Stereotyped cell divisions within a specialized cell lineage of the epidermis generate stomata and define the pattern of their distribution. The behavior of the stomatal lineage varies in its detail among different plant groups, but general features include asymmetric cell divisions and an immediate precursor (the guard mother cell [GMC]) that divides symmetrically to form the pair of cells that will differentiate into the guard cells. In Arabidopsis, the closely related basic helix-loop-helix (bHLH) subgroup Ia transcription factors SPEECHLESS, MUTE, and FAMA promote asymmetric divisions, the acquisition of GMC identity and guard cell differentiation, respectively. Genome sequence data indicate that these key positive regulators of stomatal development are broadly conserved among land plants. While orthologies can be established among individual family members within the angiosperms, more distantly related groups contain subgroup Ia bHLHs of unclear affinity. We demonstrate group Ia members from the moss Physcomitrella patens can partially complement MUTE and FAMA and recapitulate gain of function phenotypes of group Ia genes in multiple steps in the stomatal lineage in Arabidopsis. Our data are consistent with a mechanism whereby a multifunctional transcription factor underwent duplication followed by specialization to provide the three (now nonoverlapping) functions of the angiosperm stomatal bHLHs.

Arabidopsis as a model for wood formation

Wood (secondary xylem) is one of the most important sustainable energy sources for humans. Arabidopsis, despite its herbaceous nature, has become an excellent model to study wood formation. Recent progress has shown that conserved molecular mechanisms may exist in herbaceous plants and trees during vascular development and wood formation. Several transcription factor families and plant hormone species as well as other factors contribute to the regulation of xylem development in both Arabidopsis and woody plants. In this review, we highlight how information gained from the analysis of vascular development in Arabidopsis has improved our understanding of wood formation in trees.

Analysis of TETRAKETIDE α-PYRONE REDUCTASE function in Arabidopsis thaliana reveals a previously unknown, but conserved, biochemical pathway in sporopollenin monomer biosynthesis

The precise structure of the sporopollenin polymer that is the major constituent of exine, the outer pollen wall, remains poorly understood. Recently, characterization of Arabidopsis thaliana genes and corresponding enzymes involved in exine formation has demonstrated the role of fatty acid derivatives as precursors of sporopollenin building units. Fatty acyl-CoA esters synthesized by ACYL-COA SYNTHETASE5 (ACOS5) are condensed with malonyl-CoA by POLYKETIDE SYNTHASE A (PKSA) and PKSB to yield α-pyrone polyketides required for exine formation. Here, we show that two closely related genes encoding oxidoreductases are specifically and transiently expressed in tapetal cells during microspore development in Arabidopsis anthers. Mutants compromised in expression of the reductases displayed a range of pollen exine layer defects, depending on the mutant allele. Phylogenetic studies indicated that the two reductases belong to a large reductase/dehydrogenase gene family and cluster in two distinct clades with putative orthologs from several angiosperm lineages and the moss Physcomitrella patens. Recombinant proteins produced in bacteria reduced the carbonyl function of tetraketide α-pyrone compounds synthesized by PKSA/B, and the proteins were therefore named TETRAKETIDE α-PYRONE REDUCTASE1 (TKPR1) and TKPR2 (previously called DRL1 and CCRL6, respectively). TKPR activities, together with those of ACOS5 and PKSA/B, identify a conserved biosynthetic pathway leading to hydroxylated α-pyrone compounds that were previously unknown to be sporopollenin precursors.

Evolutionary conservation of the transcriptional network regulating secondary cell wall biosynthesis

The ability to make secondary cell walls was a pivotal step for vascular plants in their conquest of dry land. Here, we review recent molecular and genetic studies that reveal that a group of Arabidopsis (Arabidopsis thaliana) secondary wall-associated NAC domain transcription factors are master switches regulating a cascade of downstream transcription factors, leading to activation of the secondary wall biosynthetic program. Close homologs of the Arabidopsis secondary wall NACs and their downstream transcription factors exist in diverse taxa of vascular plants and some are functional orthologs of their Arabidopsis counterparts. There is evidence to suggest that the secondary wall NAC-mediated transcriptional regulation of secondary wall biosynthesis is a conserved mechanism throughout vascular plants.

Genome-wide analysis of plant metal transporters, with an emphasis on poplar

The specific transport of metal ions, mediated by membrane-localized metal transporters, is of fundamental importance in all eukaryotes. Genome-wide analysis of metal transporters was undertaken, making use of whole genome sequences of the green alga Chlamydomonas reinhardtii, the moss Physcomitrella patens, the lycophyte Selaginella moellendorffii, the monocots rice and sorghum, and the dicots Arabidopsis thaliana, poplar, grapevine, as well as of the yeast Saccharomyces cerevisiae. A repertoire of 430 metal transporters was found in total across eight photosynthetic plants, as well as in S. cerevisiae. Seventy-two full-length metal transporter genes were identified in the Populus genome alone, which is the largest number of metal transporters genes identified in any single species to date. Diversification of some transporter family gene clusters appears to have occurred in a lineage-specific manner. Expression analysis of Populus metal transporters indicates that some family members show tissue-specific transcript abundance. Taken together, the data provide a picture into the diversification of these important gene families.

Endogenous overexpression of Populus MYB186 increases trichome density, improves insect pest resistance, and impacts plant growth

Trichomes are specialized epidermal cells that generally play a role in reducing transpiration and act as a deterrent to herbivory. In a screen of activation-tagged Populus tremula × Populus alba 717-1B4 trees, we identified a mutant line, fuzzy, with increased foliar trichome density. This mutant also had a 35% increase in growth rate and a 200% increase in the rate of photosynthesis as compared with wild-type poplar. The fuzzy mutant had significant resistance to feeding by larvae of the white-spotted tussock moth (Orgyia leucostigma), a generalist insect pest of poplar trees. The fuzzy trichome phenotype is attributable to activation tagging and increased expression of the gene encoding PtaMYB186, which is related to Arabidopsis thaliana MYB106, a known regulator of trichome initiation. The fuzzy phenotype can be recapitulated by overexpressing PtaMYB186 in poplar. PtaMYB186 overexpression results in reconfiguration of the poplar transcriptome, with changes in the transcript abundance of suites of genes that are related to trichome differentiation. It is notable that a plant with misexpression of a gene responsible for trichome development also had altered traits related to growth rate and pest resistance, suggesting that non-intuitive facets of plant development might be useful targets for plant improvement.

Members of the LATERAL ORGAN BOUNDARIES DOMAIN transcription factor family are involved in the regulation of secondary growth in Populus

Regulation of secondary (woody) growth is of substantial economic and environmental interest but is poorly understood. We identified and subsequently characterized an activation-tagged poplar (Populus tremula × Populus alba) mutant with enhanced woody growth and changes in bark texture caused primarily by increased secondary phloem production. Molecular characterization of the mutation through positioning of the tag and retransformation experiments shows that the phenotype is conditioned by activation of an uncharacterized gene that encodes a novel member of the LATERAL ORGAN BOUNDARIES DOMAIN (LBD) family of transcription factors. Homology analysis showed highest similarity to an uncharacterized LBD1 gene from Arabidopsis thaliana, and we consequently named it Populus tremula × Populus alba (Pta) LBD1. Dominant-negative suppression of Pta LBD1 via translational fusion with the repressor SRDX domain caused decreased diameter growth and suppressed and highly irregular phloem development. In wild-type plants, LBD1 was most highly expressed in the phloem and cambial zone. Two key Class I KNOTTED1-like homeobox genes that promote meristem identity in the cambium were downregulated, while an Altered Phloem Development gene that is known to promote phloem differentiation was upregulated in the mutant. A set of four LBD genes, including the LBD1 gene, was predominantly expressed in wood-forming tissues, suggesting a broader regulatory role of these transcription factors during secondary woody growth in poplar.

MYB transcription factors in Arabidopsis

The MYB family of proteins is large, functionally diverse and represented in all eukaryotes. Most MYB proteins function as transcription factors with varying numbers of MYB domain repeats conferring their ability to bind DNA. In plants, the MYB family has selectively expanded, particularly through the large family of R2R3-MYB. Members of this family function in a variety of plant-specific processes, as evidenced by their extensive functional characterization in Arabidopsis (Arabidopsis thaliana). MYB proteins are key factors in regulatory networks controlling development, metabolism and responses to biotic and abiotic stresses. The elucidation of MYB protein function and regulation that is possible in Arabidopsis will provide the foundation for predicting the contributions of MYB proteins to the biology of plants in general.

Subgroup 4 R2R3-MYBs in conifer trees: gene family expansion and contribution to the isoprenoid- and flavonoid-oriented responses

Transcription factors play a fundamental role in plants by orchestrating temporal and spatial gene expression in response to environmental stimuli. Several R2R3-MYB genes of the Arabidopsis subgroup 4 (Sg4) share a C-terminal EAR motif signature recently linked to stress response in angiosperm plants. It is reported here that nearly all Sg4 MYB genes in the conifer trees Picea glauca (white spruce) and Pinus taeda (loblolly pine) form a monophyletic clade (Sg4C) that expanded following the split of gymnosperm and angiosperm lineages. Deeper sequencing in P. glauca identified 10 distinct Sg4C sequences, indicating over-representation of Sg4 sequences compared with angiosperms such as Arabidopsis, Oryza, Vitis, and Populus. The Sg4C MYBs share the EAR motif core. Many of them had stress-responsive transcript profiles after wounding, jasmonic acid (JA) treatment, or exposure to cold in P. glauca and P. taeda, with MYB14 transcripts accumulating most strongly and rapidly. Functional characterization was initiated by expressing the P. taeda MYB14 (PtMYB14) gene in transgenic P. glauca plantlets with a tissue-preferential promoter (cinnamyl alcohol dehydrogenase) and a ubiquitous gene promoter (ubiquitin). Histological, metabolite, and transcript (microarray and targeted quantitative real-time PCR) analyses of PtMYB14 transgenics, coupled with mechanical wounding and JA application experiments on wild-type plantlets, allowed identification of PtMYB14 as a putative regulator of an isoprenoid-oriented response that leads to the accumulation of sesquiterpene in conifers. Data further suggested that PtMYB14 may contribute to a broad defence response implicating flavonoids. This study also addresses the potential involvement of closely related Sg4C sequences in stress responses and plant evolution.

Comprehensive analysis of NAC domain transcription factor gene family in Populus trichocarpa

BACKGROUND

NAC (NAM, ATAF1/2 and CUC2) domain proteins are plant-specific transcriptional factors known to play diverse roles in various plant developmental processes. NAC transcription factors comprise of a large gene family represented by more than 100 members in Arabidopsis, rice and soybean etc. Recently, a preliminary phylogenetic analysis was reported for NAC gene family from 11 plant species. However, no comprehensive study incorporating phylogeny, chromosomal location, gene structure, conserved motifs, and expression profiling analysis has been presented thus far for the model tree species Populus.

RESULTS

In the present study, a comprehensive analysis of NAC gene family in Populus was performed. A total of 163 full-length NAC genes were identified in Populus, and they were phylogenetically clustered into 18 distinct subfamilies. The gene structure and motif compositions were considerably conserved among the subfamilies. The distributions of 120 Populus NAC genes were non-random across the 19 linkage groups (LGs), and 87 genes (73%) were preferentially retained duplicates that located in both duplicated regions. The majority of NACs showed specific temporal and spatial expression patterns based on EST frequency and microarray data analyses. However, the expression patterns of a majority of duplicate genes were partially redundant, suggesting the occurrence of subfunctionalization during subsequent evolutionary process. Furthermore, quantitative real-time RT-PCR (RT-qPCR) was performed to confirm the tissue-specific expression patterns of 25 NAC genes.

CONCLUSION

Based on the genomic organizations, we can conclude that segmental duplications contribute significantly to the expansion of Populus NAC gene family. The comprehensive expression profiles analysis provides first insights into the functional divergence among members in NAC gene family. In addition, the high divergence rate of expression patterns after segmental duplications indicates that NAC genes in Populus are likewise to have been retained by substantial subfunctionalization. Taken together, our results presented here would be helpful in laying the foundation for functional characterization of NAC gene family and further gaining an understanding of the structure-function relationship between these family members.

Genome-wide classification and evolutionary analysis of the bHLH family of transcription factors in Arabidopsis, poplar, rice, moss, and algae

Basic helix-loop-helix proteins (bHLHs) are found throughout the three eukaryotic kingdoms and constitute one of the largest families of transcription factors. A growing number of bHLH proteins have been functionally characterized in plants. However, some of these have not been previously classified. We present here an updated and comprehensive classification of the bHLHs encoded by the whole sequenced genomes of Arabidopsis (Arabidopsis thaliana), Populus trichocarpa, Oryza sativa, Physcomitrella patens, and five algae species. We define a plant bHLH consensus motif, which allowed the identification of novel highly diverged atypical bHLHs. Using yeast two-hybrid assays, we confirm that (1) a highly diverged bHLH has retained protein interaction activity and (2) the two most conserved positions in the consensus play an essential role in dimerization. Phylogenetic analysis permitted classification of the 638 bHLH genes identified into 32 subfamilies. Evolutionary and functional relationships within subfamilies are supported by intron patterns, predicted DNA-binding motifs, and the architecture of conserved protein motifs. Our analyses reveal the origin and evolutionary diversification of plant bHLHs through differential expansions, domain shuffling, and extensive sequence divergence. At the functional level, this would translate into different subfamilies evolving specific DNA-binding and protein interaction activities as well as differential transcriptional regulatory roles. Our results suggest a role for bHLH proteins in generating plant phenotypic diversity and provide a solid framework for further investigations into the role carried out in the transcriptional regulation of key growth and developmental processes.

Coordinated transcriptional regulation of two key genes in the lignin branch pathway--CAD and CCR--is mediated through MYB- binding sites

BACKGROUND

Cinnamoyl CoA reductase (CCR) and cinnamyl alcohol dehydrogenase (CAD) catalyze the final steps in the biosynthesis of monolignols, the monomeric units of the phenolic lignin polymers which confer rigidity, imperviousness and resistance to biodegradation to cell walls. We have previously shown that the Eucalyptus gunnii CCR and CAD2 promoters direct similar expression patterns in vascular tissues suggesting that monolignol production is controlled, at least in part, by the coordinated transcriptional regulation of these two genes. Although consensus motifs for MYB transcription factors occur in most gene promoters of the whole phenylpropanoid pathway, functional evidence for their contribution to promoter activity has only been demonstrated for a few of them. Here, in the lignin-specific branch, we studied the functional role of MYB elements as well as other cis-elements identified in the regulatory regions of EgCAD2 and EgCCR promoters, in the transcriptional activity of these gene promoters.

RESULTS

By using promoter deletion analysis and in vivo footprinting, we identified an 80 bp regulatory region in the Eucalyptus gunnii EgCAD2 promoter that contains two MYB elements, each arranged in a distinct module with newly identified cis-elements. A directed mutagenesis approach was used to introduce block mutations in all putative cis-elements of the EgCAD2 promoter and in those of the 50 bp regulatory region previously delineated in the EgCCR promoter. We showed that the conserved MYB elements in EgCAD2 and EgCCR promoters are crucial both for the formation of DNA-protein complexes in EMSA experiments and for the transcriptional activation of EgCAD2 and EgCCR promoters in vascular tissues in planta. In addition, a new regulatory cis-element that modulates the balance between two DNA-protein complexes in vitro was found to be important for EgCAD2 expression in the cambial zone.

CONCLUSIONS

Our assignment of functional roles to the identified cis-elements clearly demonstrates the importance of MYB cis-elements in the transcriptional regulation of two genes of the lignin-specific pathway and support the hypothesis that MYB elements serve as a common means for the coordinated regulation of genes in the entire lignin biosynthetic pathway.

The poplar MYB transcription factors, PtrMYB3 and PtrMYB20, are involved in the regulation of secondary wall biosynthesis

Dicot wood is mainly composed of cellulose, xylan and lignin, and its formation requires the coordinated regulation of their biosynthesis. In this report, we demonstrate that the poplar wood-associated MYB transcriptional activators, PtrMYB3 and PtrMYB20, activate the biosynthetic pathways of cellulose, xylan and lignin when overexpressed in Arabidopsis and they are also able to activate the promoter activities of poplar wood biosynthetic genes. We also show that PtrMYB3 and PtrMYB20 are functional orthologs of Arabidopsis MYB46 and MYB83, and their expression is directly activated by poplar PtrWND2, suggesting their involvement in the regulation of wood formation in poplar.

The poplar PtrWNDs are transcriptional activators of secondary cell wall biosynthesis

Secondary cell walls, consisting of cellulose, hemicelluloses and lignin, make up the bulk of wood biomass. It is therefore expected that dissection of the molecular mechanisms underlying secondary wall biosynthesis and its regulation will be instrumental to unravel the process of wood formation in tree species. Wood formation requires the coordinated activation of genes in the secondary wall biosynthetic program that is essential for the biosynthesis and assembly of wood components. It has recently been discovered that a group of poplar (Populus trichocarpa) wood-associated NAC domain transcription factors, PtrWNDs, which are functional orthologs of the Arabidopsis SND1, are capable of turning on the entire secondary wall biosynthetic program when expressed in Arabidopsis. In addition, two of the PtrWNDs were found to be able to activate the promoters of poplar wood biosynthetic genes and a number of other poplar wood-associated transcription factors. Further testing reveals that the promoters of these poplar wood-associated transcription factors are also activated by other PtrWNDs. It is therefore proposed that PtrWNDs are master transcriptional switches regulating a cascade of downstream transcription factors and thereby mediate the coordinated activation of wood biosynthetic genes during wood formation.

R2R3-NaMYB8 regulates the accumulation of phenylpropanoid-polyamine conjugates, which are essential for local and systemic defense against insect herbivores in Nicotiana attenuata

Although phenylpropanoid-polyamine conjugates (PPCs) occur ubiquitously in plants, their biological roles remain largely unexplored. The two major PPCs of Nicotiana attenuata plants, caffeoylputrescine (CP) and dicaffeoylspermidine, increase dramatically in local and systemic tissues after herbivore attack and simulations thereof. We identified NaMYB8, a homolog of NtMYBJS1, which in BY-2 cells regulates PPC biosynthesis, and silenced its expression by RNA interference in N. attenuata (ir-MYB8), to understand the ecological role(s) of PPCs. The regulatory role of NaMYB8 in PPC biosynthesis was validated by a microarray analysis, which revealed that transcripts of several key biosynthetic genes in shikimate and polyamine metabolism accumulated in a NaMYB8-dependent manner. Wild-type N. attenuata plants typically contain high levels of PPCs in their reproductive tissues; however, NaMYB8-silenced plants that completely lacked CP and dicaffeoylspermidine showed no changes in reproductive parameters of the plants. In contrast, a defensive role for PPCs was clear; both specialist (Manduca sexta) and generalist (Spodoptera littoralis) caterpillars feeding on systemically preinduced young stem leaves performed significantly better on ir-MYB8 plants lacking PPCs compared with wild-type plants expressing high levels of PPCs. Moreover, the growth of M. sexta caterpillars was significantly reduced when neonates were fed ir-MYB8 leaves sprayed with synthetic CP, corroborating the role of PPCs as direct plant defense. The spatiotemporal accumulation and function of PPCs in N. attenuata are consistent with the predictions of the optimal defense theory: plants preferentially protect their most fitness-enhancing and vulnerable parts, young tissues and reproductive organs, to maximize their fitness.

Functional characterization of poplar wood-associated NAC domain transcription factors

Wood is the most abundant biomass produced by land plants. Dissection of the molecular mechanisms underlying the transcriptional regulation of wood formation is a fundamental issue in plant biology and has important implications in tree biotechnology. Although a number of transcription factors in tree species have been shown to be associated with wood formation and some of them are implicated in lignin biosynthesis, none of them have been demonstrated to be key regulators of the biosynthesis of all three major components of wood. In this report, we have identified a group of NAC domain transcription factors, PtrWNDs, that are preferentially expressed in developing wood of poplar (Populus trichocarpa). Expression of PtrWNDs in the Arabidopsis (Arabidopsis thaliana) snd1 nst1 double mutant effectively complemented the secondary wall defects in fibers, indicating that PtrWNDs are capable of activating the entire secondary wall biosynthetic program. Overexpression of PtrWND2B and PtrWND6B in Arabidopsis induced the expression of secondary wall-associated transcription factors and secondary wall biosynthetic genes and, concomitantly, the ectopic deposition of cellulose, xylan, and lignin. Furthermore, PtrWND2B and PtrWND6B were able to activate the promoter activities of a number of poplar wood-associated transcription factors and wood biosynthetic genes. Together, these results demonstrate that PtrWNDs are functional orthologs of SND1 and suggest that PtrWNDs together with their downstream transcription factors form a transcriptional network involved in the regulation of wood formation in poplar.

Extensive structural renovation of retrogenes in the evolution of the Populus genome

Retroposition, as an important copy mechanism for generating new genes, was believed to play a negligible role in plants. As a representative dicot, the genomic sequences of Populus (poplar; Populus trichocarpa) provide an opportunity to investigate this issue. We identified 106 retrogenes and found the majority (89%) of them are associated with functional signatures in sequence evolution, transcription, and (or) translation. Remarkably, examination of gene structures revealed extensive structural renovation of these retrogenes: we identified 18 (17%) of them undergoing either chimerization to form new chimerical genes and (or) intronization (transformation into intron sequences of previously exonic sequences) to generate new intron-containing genes. Such a change might occur at a high speed, considering eight out of 18 such cases occurred recently after divergence between Arabidopsis (Arabidopsis thaliana) and Populus. This pattern also exists in Arabidopsis, with 15 intronized retrogenes occurring after the divergence between Arabidopsis and papaya (Carica papaya). Thus, the frequency of intronization in dicots revealed its importance as a mechanism in the evolution of exon-intron structure. In addition, we also examined the potential impact of the Populus nascent sex determination system on the chromosomal distribution of retrogenes and did not observe any significant effects of the extremely young sex chromosomes.

Genome-wide analysis of major intrinsic proteins in the tree plant Populus trichocarpa: characterization of XIP subfamily of aquaporins from evolutionary perspective

BACKGROUND

Members of major intrinsic proteins (MIPs) include water-conducting aquaporins and glycerol-transporting aquaglyceroporins. MIPs play important role in plant-water relations. The model plants Arabidopsis thaliana, rice and maize contain more than 30 MIPs and based on phylogenetic analysis they can be divided into at least four subfamilies. Populus trichocarpa is a model tree species and provides an opportunity to investigate several tree-specific traits. In this study, we have investigated Populus MIPs (PtMIPs) and compared them with their counterparts in Arabidopsis, rice and maize.

RESULTS

Fifty five full-length MIPs have been identified in Populus genome. Phylogenetic analysis reveals that Populus has a fifth uncharacterized subfamily (XIPs). Three-dimensional models of all 55 PtMIPs were constructed using homology modeling technique. Aromatic/arginine (ar/R) selectivity filters, characteristics of loops responsible for solute selectivity (loop C) and gating (loop D) and group conservation of small and weakly polar interfacial residues have been analyzed. Majority of the non-XIP PtMIPs are similar to those in Arabidopsis, rice and maize. Additional XIPs were identified from database search and 35 XIP sequences from dicots, fungi, moss and protozoa were analyzed. Ar/R selectivity filters of dicots XIPs are more hydrophobic compared to fungi and moss XIPs and hence they are likely to transport hydrophobic solutes. Loop C is longer in one of the subgroups of dicot XIPs and most probably has a significant role in solute selectivity. Loop D in dicot XIPs has higher number of basic residues. Intron loss is observed on two occasions: once between two subfamilies of eudicots and monocot and in the second instance, when dicot and moss XIPs diverged from fungi. Expression analysis of Populus MIPs indicates that Populus XIPs don't show any tissue-specific transcript abundance.

CONCLUSION

Due to whole genome duplication, Populus has the largest number of MIPs identified in any single species. Non-XIP MIPs are similar in all four plant species considered in this study. Small and weakly polar residues at the helix-helix interface are group conserved presumably to maintain the hourglass fold of MIP channels. Substitutions in ar/R selectivity filter, insertion/deletion in loop C, increasing basic nature of loop D and loss of introns are some of the events occurred during the evolution of dicot XIPs.

Genotype and time of day shape the Populus drought response

As exposure to episodic drought can impinge significantly on forest health and the establishment of productive tree plantations, there is great interest in understanding the mechanisms of drought response in trees. The ecologically dominant and economically important genus Populus, with its sequenced genome, provides an ideal opportunity to examine transcriptome level changes in trees in response to a drought stimulus. The transcriptome level drought response of two commercially important Populus clones (P. deltoides x P. nigra, DN34, and P. nigra x P. maximowiczii, NM6) was characterized over a diurnal period using a 4 x 2 x 2 complete randomized factorial anova experimental design (four time points, two genotypes and two treatment conditions), using Affymetrix Poplar GeneChip microarrays. Notably, the specific genes that exhibited changes in transcript abundance in response to drought differed between the genotypes and/or the time of day that they exhibited their greatest differences. This study emphasizes the fact that it is not possible to draw simple, generalized conclusions about the drought response of the genus Populus on the basis of one species, nor on the basis of results collected at a single time point. The data derived from our studies provide insights into the variety of genetic mechanisms underpinning the Populus drought response, and provide candidates for future experiments aimed at understanding this response across this economically and ecologically important genus.

Post-translational modification of an R2R3-MYB transcription factor by a MAP Kinase during xylem development

Despite the pivotal role played by R2R3-MYB family members in the regulation of plant gene expression, little is known about post-translational regulation of these proteins. In animals, the MYB family member, c-MYB, is post-translationally modified by a mitogen-activated protein kinase (MAPK), p42(mapk). In order to test the hypothesis that R2R3-MYB proteins may be regulated by MAPK activity, interplay between a R2R3-MYB family member expressed in differentiating pine xylem (Pinus taeda MYB4, PtMYB4) and MAPK proteins expressed in the same tissue was examined. One of the MAPK proteins expressed in pine xylem, PtMAPK6, phosphorylated PtMYB4. Recombinant PtMAPK6 phosphorylated PtMYB4 on serine-236, located in the C-terminal activation domain of this transcription factor in a context that is found in other plant MYB proteins. Modification of the PtMAPK6 target serine in PtMYB4 did not appear to alter DNA binding in vitro but did alter the ability of PtMYB4 to promote transcriptional activation in yeast. PtMAPK6 activity was detected in developing xylem cells that had ceased cell division and formed secondary walls. Together, the data support a role for PtMAPK6 during early xylem development and suggest a function for this kinase in regulating gene expression through phosphorylation of PtMYB4.

The wound-, pathogen-, and ultraviolet B-responsive MYB134 gene encodes an R2R3 MYB transcription factor that regulates proanthocyanidin synthesis in poplar

In poplar (Populus spp.), the major defense phenolics produced in leaves are the flavonoid-derived proanthocyanidins (PAs) and the salicin-based phenolic glycosides. Transcriptional activation of PA biosynthetic genes leading to PA accumulation in leaves occurs following herbivore damage and mechanical wounding as well as infection by the fungal biotroph Melampsora medusae. In this study, we have identified a poplar R2R3 MYB transcription factor gene, MYB134, that exhibits close sequence similarity to the Arabidopsis (Arabidopsis thaliana) PA regulator TRANSPARENT TESTA2 and that is coinduced with PA biosynthetic genes following mechanical wounding, M. medusae infection, and exposure to elevated ultraviolet B light. Overexpression of MYB134 in poplar resulted in transcriptional activation of the full PA biosynthetic pathway and a significant plant-wide increase in PA levels, and electrophoretic mobility shift assays showed that recombinant MYB134 protein is able to bind to promoter regions of PA pathway genes. MYB134-overexpressing plants exhibited a concomitant reduction in phenolic glycoside concentrations and other minor alterations to levels of small phenylpropanoid metabolites. Our data provide insight into the regulatory mechanisms controlling stress-induced PA metabolism in poplar, and the identification of a regulator of stress-responsive PA biosynthesis constitutes a valuable tool for manipulating PA metabolism in poplar and investigating the biological functions of PAs in resistance to biotic and abiotic stresses.