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

PdMYB118, isolated from a red leaf mutant of Populus deltoids, is a new transcription factor regulating anthocyanin biosynthesis in poplar

A new anthocyanin biosynthesis transcription factor PdMYB118, which could be used for the genetic engineering of colorful tree species, was indentified from a red leaf mutant of Populus deltoids. In higher plants, the biosynthesis of anthocyanins is regulated by several classes of transcription factors (TFs), including R2R3-MYB, bHLH and WD-repeat proteins. In this work, we isolated an MYB gene regulating anthocyanin biosynthesis from a red leaf mutant of Populus deltoids, which accumulated more anthocyanins in the leaves and showed higher expression levels of anthocyanin biosynthesis genes than did the wild type. Gene expression analyses of all TFs regulating anthocyanin biosynthesis demonstrated that only a MYB118 homologous gene, PdMYB118, was up-regulated in the mutant compared with the wide type. Subcellular localization analyses in poplar leaf mesophyll protoplasts showed that PdMYB118-YFP fusion protein was specifically located in nucleus. When transiently expressed in poplar leaf protoplasts, PdMYB118 specifically promoted the expression of anthocyanidin biosynthesis genes. Dual-luciferase assays revealed that PdMYB118 can directly activate the promoters of these genes. When overexpressed in Shanxin Yang (P. davidiana × P. bolleana), a hybrid clone commercially grown for landscaping in the northern part of China, transgenic plants overexpressing PdMYB118 produced more anthocyanins in the leaves and turned their color into redness when grown in both greenhouse and field. Consistently, transcripts of some important anthocyanidin biosynthesis genes were significantly increased in the leaves of transgenic plants. All these results indicate that PdMYB118 functions as an essential transcription factor regulating anthocyanin biosynthesis in poplar and could be used for the genetic engineering of colorful tree species.

MADS-box genes underground becoming mainstream: plant root developmental mechanisms

Plant growth is largely post-embryonic and depends on meristems that are active throughout the lifespan of an individual. Developmental patterns rely on the coordinated spatio-temporal expression of different genes, and the activity of transcription factors is particularly important during most morphogenetic processes. MADS-box genes constitute a transcription factor family in eukaryotes. In Arabidopsis, their proteins participate in all major aspects of shoot development, but their role in root development is still not well characterized. In this review we synthetize current knowledge pertaining to the function of MADS-box genes highly expressed in roots: XAL1, XAL2, ANR1 and AGL21, as well as available data for other MADS-box genes expressed in this organ. The role of Trithorax group and Polycomb group complexes on MADS-box genes' epigenetic regulation is also discussed. We argue that understanding the role of MADS-box genes in root development of species with contrasting architectures is still a challenge. Finally, we propose that MADS-box genes are key components of the gene regulatory networks that underlie various gene expression patterns, each one associated with the distinct developmental fates observed in the root. In the case of XAL1 and XAL2, their role within these networks could be mediated by regulatory feedbacks with auxin.

The R2R3 MYB transcription factor MYB189 negatively regulates secondary cell wall biosynthesis in Populus

Secondary cell wall (SCW) biosynthesis during wood formation in trees is controlled by a multilevel regulatory network that coordinates the expression of substantial genes. However, few transcription factors involved in the negative regulation of secondary wall biosynthesis have been characterized in tree species. In this study, we isolated an R2R3 MYB transcription factor MYB189 from Populus trichocarpa, which is expressed predominantly in secondary vascular tissues, especially in the xylem. A novel repression motif was identified in the C-terminal region of MYB189, which indicates this factor was a transcriptional repressor. Overexpression (OE) of MYB189 in Arabidopsis and poplar resulted in a significant reduction in the contents of lignin, cellulose and hemicelluloses. Vascular development in stems of MYB189 OE lines was markedly inhibited, leading to a dramatic decrease in SCW thickness of xylem fibers. Gene expression analyses showed that most of the structural genes involved in the biosynthesis of lignin, cellulose and xylans were significantly downregulated in MYB189-overexpressing poplars compared with the wild-type control. Chromatin immunoprecipitation-quantitative real-time polymerase chain reaction and transient expression assays revealed that MYB189 could directly bind to the promoters of secondary wall biosynthetic genes to repress their expression. Together, these data suggest that MYB189 acts as a repressor to regulate SCW biosynthesis in poplar.

Overexpression of CmMYB15 provides chrysanthemum resistance to aphids by regulating the biosynthesis of lignin

MYB transcription factors are widely involved in the development of and physiological processes in plants. Here, we isolated the chrysanthemum R2R3-MYB family transcription factor CmMYB15, a homologous gene of AtMYB15. It was demonstrated that CmMYB15 expression was induced by aphids and that CmMYB15 could bind to AC elements, which usually exist in the promoter of lignin biosynthesis genes. Overexpression of CmMYB15 in chrysanthemum enhanced the resistance of aphids. Additionally, the content of lignin and the expression of several lignin biosynthesis genes increased. In summary, the results indicate that CmMYB15 regulates lignin biosynthesis genes that enhance the resistance of chrysanthemum to aphids.

Genome-wide characterization of aspartic protease (AP) gene family in Populus trichocarpa and identification of the potential PtAPs involved in wood formation

BACKGROUND

Aspartic protease (AP) is one of four large proteolytic enzyme families that are involved in plant growth and development. Little is known about the AP gene family in tree species, although it has been characterized in Arabidopsis, rice and grape. The AP genes that are involved in tree wood formation remain to be determined.

RESULTS

A total of 67 AP genes were identified in Populus trichocarpa (PtAP) and classified into three categories (A, B and C). Chromosome mapping analysis revealed that two-thirds of the PtAP genes were located in genome duplication blocks, indicating the expansion of the AP family by segmental duplications in Populus. The microarray data from the Populus eFP browser demonstrated that PtAP genes had diversified tissue expression patterns. Semi-qRT-PCR analysis further determined that more than 10 PtAPs were highly or preferentially expressed in the developing xylem. When the involvement of the PtAPs in wood formation became the focus, many SCW-related cis-elements were found in the promoters of these PtAPs. Based on PtAP_promoter::GUS techniques, the activities of PtAP66 promoters were observed only in fiber cells, not in the vessels of stems as the xylem and leaf veins developed in the transgenic Populus tree, and strong GUS signals were detected in interfascicular fiber cells, roots, anthers and sepals of PtAP17_promoter::GUS transgenic plants. Intensive GUS activities in various secondary tissues implied that PtAP66 and PtAP17 could function in wood formation. In addition, most of the PtAP proteins were predicted to contain N- and (or) O-glycosylation sites, and the integration of PNGase F digestion and western blotting revealed that the PtAP17 and PtAP66 proteins were N-glycosylated in Populus.

CONCLUSIONS

Comprehensive characterization of the PtAP genes suggests their functional diversity during Populus growth and development. Our findings provide an overall understanding of the AP gene family in trees and establish a better foundation to further describe the roles of PtAPs in wood formation.

Genome-wide identification and expression analysis of R2R3, 3R- and 4R-MYB transcription factors during lignin biosynthesis in flax (Linum usitatissimum)

MYB transcription factors (TFs) have vital roles in regulating lignin or fiber development. Flax (Linum usitatissimum) is known as one of the plants with high fiber production capacity. However, no studies have been conducted to identify and characterize MYB TFs in the flax genome. Results showed that flax genome harbours 167 R2R3, seven 3R, and one 4R-type MYB TFs. 22 MYB genes (%13) were estimated to be tandem duplicated dated around 13.3-86.98 Mya. 130 flax MYB members have apparent orthologous with Arabidopsis, in which 17 R2R3 MYBs are associated with lignin biosynthesis. MYB062, MYB072, MYB096, MYB141, and MYB146 genes were up-regulated in tissues having higher lignin production capacity. In opposite, MYB012 and MYB113 genes were down-regulated which points out the involvement of those genes in the lignin biosynthesis mechanism. This comprehensive study can provide a basis for understanding the role of MYBs in fiber or lignin production in flax.

A novel MYB transcription factor regulates ascorbic acid synthesis and affects cold tolerance

Dehydroascorbate reductase (DHAR) plays an important role in stress responses, but the transcriptional regulation of DHAR in response to abiotic stress is still poorly understood. In this study, we isolated a novel R2R3-type MYB transcription factor from Pyrus betulaefolia by yeast one-hybrid screening, designated as PbrMYB5. PbrMYB5 was localized in the nucleus and could bind specifically to the promoter of PbrDHAR2. PbrMYB5 was greatly induced by cold and salt but slightly by dehydration. Overexpression of PbrMYB5 in tobacco conferred enhanced tolerance to chilling stresses, whereas down-regulation of PbrMYB5 in P. betulaefolia by virus-induced gene silencing resulted in elevated chilling sensitivity. Transgenic tobacco exhibited higher expression levels of NtDHAR2 and accumulated larger amount of ascorbic acid (AsA) than the wild-type plants. Virus-induced gene silencing of PbrMYB5 in P. betulaefolia down-regulated PbrDHAR2 abundance and decreased AsA level, accompanied by an increased sensitivity to the chilling stress. Taken together, these results demonstrated that PbrMYB5 was an activator of AsA biosynthesis and may play a positive role in chilling tolerance, at least in part, due to the modulation of AsA synthesis by regulating the PbrDHAR2 expression.

Aquaporins Respond to Chilling in the Phloem by Altering Protein and mRNA Expression

Previous experiments using heat exchangers (liquid cooled blocks) to chill a portion of plant stem have shown a transient stoppage in phloem translocation and an increase in measured phloem pressure. Although a chilled-induced stoppage of phloem transport has been known for over 100 years, the mechanism of this phenomenon is still poorly understood. Recently, work has highlighted that aquaporins occur within the plasma membrane of the sieve tubes along the entire source-to-sink pathway, and that isoforms of these water channel proteins may change dynamically. Aquaporins show regulatory roles in controlling tissue and cellular water status in response to environmental hardships. Thus, we tested if protein localization and mRNA transcript abundance changes occur in response to chilling in balsam poplar (Populus balsamifera) using immunohistochemistry and qrtPCR. The results of the immunolocalization experiments show that the labeling intensity of the sieve elements treated for only 2 min of chill time significantly increased for PIP2. After 10 min of chilling, this signal declined significantly to lower than that of the pre-chilled sieve elements. Overall, the abundance of mRNA transcript increased for the tested PIP2s following cold application. We discuss the implication that aquaporins are responsible for the alleviation of sieve tube pressure and the resumption of flow following a cold-induced blockage event.

Evolution of the modular, disordered stress proteins known as dehydrins

Dehydrins, plant proteins that are upregulated during dehydration stress conditions, have modular sequences that can contain three conserved motifs (the Y-, S-, and K-segments). The presence and order of these motifs are used to classify dehydrins into one of five architectures: Kn, SKn, KnS, YnKn, and YnSKn, where the subscript n describes the number of copies of that motif. In this study, an architectural and phylogenetic analysis was performed on 426 dehydrin sequences that were identified in 53 angiosperm and 3 gymnosperm genomes. It was found that angiosperms contained all five architectures, while gymnosperms only contained Kn and SKn dehydrins. This suggests that the ancestral dehydrin in spermatophytes was either Kn or SKn, and the Y-segment containing dehydrins first arose in angiosperms. A high-level split between the YnSKn dehydrins from either the Kn or SKn dehydrins could not be confidently identified, however, two lower level architectural divisions appear to have occurred after different duplication events. The first likely occurred after a whole genome duplication, resulting in the duplication of a Y3SK2 dehydrin; the duplicate subsequently lost an S- and K- segment to become a Y3K1 dehydrin. The second split occurred after a tandem duplication of a Y1SK2 dehydrin, where the duplicate lost both the Y- and S- segment and gained four K-segments, resulting in a K6 dehydrin. We suggest that the newly arisen Y3K1 dehydrin is possibly on its way to pseudogenization, while the newly arisen K6 dehydrin developed a novel function in cold protection.

JcMYB1, a Jatropha R2R3MYB Transcription Factor Gene, Modulates Lipid Biosynthesis in Transgenic Plants

The lipid biosynthesis pathway in plants has been studied in detail; however, the factors that regulate the pathway at the transcription level are largely unknown. LEAFY COTYLEDON1 (LEC1), WRINKLED1 (WRI1) and FUSCA3 (FUS3) are considered master regulators to control seed oil content in Arabidopsis. Beside these master regulators, several other transcription factors that may regulate the pathway in plants are poorly studied. In the present work, we have shown the involvement of an uncharacterized Jatropha curcas R2R3MYB gene (JcMYB1) in seed oil biosynthesis. Seed oil analysis and expression profiling of fatty acid (FA) and triacylglycerol (TAG) biosynthetic genes in transgenic Arabidopsis and tobacco plants revealed that JcMYB1 enhances seed oil accumulation and alters FA composition by regulating the expression of endogenous pathway genes in transgenics. Using virus-induced gene silencing (VIGS) in Jatropha, we demonstrated that the suppression of JcMYB1 reduced lipid content with altered FA composition. Agro-infiltration and yeast one-hybrid assay results showed that JcMYB1 protein directly binds to the diacylglycerol acyltransferase1 (DGAT1) promoter, a rate-limiting enzyme of TAG biosynthesis, and activates its expression. These results suggested that JcMYB1 may augment the lipid content by regulating lipid biosynthetic genes. Additionally, manipulation of JcMYB1 in oil crop plants may be used for the potential improvement of oil production and quality.

Integration of Transcriptome, Proteome, and Metabolome Provides Insights into How Calcium Enhances the Mechanical Strength of Herbaceous Peony Inflorescence Stems

Weak stem mechanical strength severely restrains cut flowers quality and stem weakness can be alleviated by calcium (Ca) treatment, but the mechanisms underlying Ca-mediated enhancement of stem mechanical strength remain largely unknown. In this study, we performed a comparative transcriptomic, proteomic, and metabolomic analysis of herbaceous peony (Paeonia lactiflora Pall.) inflorescence stems treated with nanometer Ca carbonate (Nano-CaCO₃). In total, 2643 differentially expressed genes (DEGs) and 892 differentially expressed proteins (DEPs) were detected between the Control and nano-CaCO₃ treatment. Among the 892 DEPs, 152 were coregulated at both the proteomic and transcriptomic levels, and 24 DEPs related to the secondary cell wall were involved in signal transduction, energy metabolism, carbohydrate metabolism and lignin biosynthesis, most of which were upregulated after nano-CaCO₃ treatment during the development of inflorescence stems. Among these four pathways, numerous differentially expressed metabolites (DEMs) related to lignin biosynthesis were identified. Furthermore, structural observations revealed the thickening of the sclerenchyma cell walls, and the main wall constitutive component lignin accumulated significantly in response to nano-CaCO₃ treatment, thereby indicating that Ca can enhance the mechanical strength of the inflorescence stems by increasing the lignin accumulation. These results provided insights into how Ca treatment enhances the mechanical strength of inflorescence stems in P. lactiflora.

Genome-wide identification and expression analysis of the MYB transcription factor in moso bamboo (Phyllostachys edulis)

The MYB family, one of the largest transcription factor (TF) families in the plant kingdom, plays vital roles in cell formation, morphogenesis and signal transduction, as well as responses to biotic and abiotic stresses. However, the underlying function of bamboo MYB TFs remains unclear. To gain insight into the status of these proteins, a total of 85 PeMYBs, which were further divided into 11 subgroups, were identified in moso bamboo (Phyllostachys edulis) by using a genome-wide search strategy. Gene structure analysis showed that PeMYBs were significantly different, with exon numbers varying from 4 to 13. Phylogenetic analysis indicated that PeMYBs clustered into 27 clades, of which the function of 18 clades has been predicted. In addition, almost all of the PeMYBs were differently expressed in leaves, panicles, rhizomes and shoots based on RNA-seq data. Furthermore, qRT-PCR analysis showed that 12 PeMYBs related to the biosynthesis and deposition of the secondary cell wall (SCW) were constitutively expressed, and their transcript abundance levels have changed significantly with increasing height of the bamboo shoots, for which the degree of lignification continuously increased. This result indicated that these PeMYBs might play fundamental roles in SCW thickening and bamboo shoot lignification. The present comprehensive and systematic study on the members of the MYB family provided a reference and solid foundation for further functional analysis of MYB TFs in moso bamboo.

Functional Characterization of Invertase Inhibitors PtC/VIF1 and 2 Revealed Their Involvements in the Defense Response to Fungal Pathogen in Populus trichocarpa

In higher plants, cell wall invertase (CWI) and vacuolar invertase (VI) were considered to be essential coordinators in carbohydrate partitioning, sink strength determination, and stress responses. An increasing body of evidence revealed that the tight regulation of CWI and VI substantially depends on the post-translational mechanisms, which were mediated by small proteinaceous inhibitors (C/VIFs, Inhibitor of β-Fructosidases). As yet, the extensive survey of the molecular basis and biochemical property of C/VIFs remains largely unknown in black cottonwood (Populus trichocarpa Torr. & A. Gray), a model species of woody plants. In the present work, we have initiated a systematic review of the genomic structures, phylogenies, cis-regulatory elements, and conserved motifs as well as the tissue-specific expression, resulting in the identification of 39 genes encoding C/VIF in poplar genome. We characterized two putative invertase inhibitors PtC/VIF1 and 2, showing predominant transcript levels in the roots and highly divergent responses to the selected stress cues including fusarium wilt, drought, ABA, wound, and senescence. In silico prediction of the signal peptide hinted us that they both likely had the apoplastic targets. Based on the experimental visualization via the transient and stable transformation assays, we confirmed that PtC/VIF1 and 2 indeed secreted to the extracellular compartments. Further validation of their recombinant enzymes revealed that they displayed the potent inhibitory affinities on the extracted CWI, supporting the patterns that act as the typical apoplastic invertase inhibitors. To our knowledge, it is the first report on molecular characterization of the functional C/VIF proteins in poplar. Our results indicate that PtC/VIF1 and 2 may exert essential roles in defense- and stress-related responses. Moreover, novel findings of the up- and downregulated C/VIF genes and functional enzyme activities enable us to further unravel the molecular mechanisms in the promotion of woody plant performance and adapted-biotic stress, underlying the homeostatic control of sugar in the apoplast.

Overexpression of ScMYBAS1 alternative splicing transcripts differentially impacts biomass accumulation and drought tolerance in rice transgenic plants

Drought is the most significant environmental stress for agricultural production worldwide, and tremendous efforts have been made to improve crop yield under the increasing water scarcity. Transcription factors are major players in the regulation of water stress-related genes in plants. Recently, different MYB transcription factors were characterized for their involvement in drought response. A sugarcane R2R3-MYB gene (ScMYBAS1) and its four alternative forms of transcript (ScMYAS1-2, ScMYBAS1-3, ScMYBAS1-4 and ScMYBAS1-5) were identified in this study. The subcellular localization, in Nicotiniana benthamiana, of the TFs fused in frame with GFP revealed that ScMYBAS1-2-GFP and ScMYBAS1-3-GFP were observed in the nucleus. The overexpression of ScMYBAS1-2 and ScMYBAS1-3 spliced transcripts in rice promoted change in plant growth under both well-watered and drought conditions. The ScMYBAS1-2 and ScMYBAS1-3 transgenic lines revealed a higher relative water content (RWC) compared to the wild type before maximum stress under drought conditions. The ScMYBAS1-2 transgenic lines showed a reduction in biomass (total dry weight). Conversely, ScMYBAS1-3 showed an increased biomass (total dry weight) relative to the wild-type. The overexpression of ScMYBAS1-3 in rice transgenic lines showed involvement with drought tolerance and biomass and, for this reason, was considered a good target for plant transformation, particularly for use in developing genotypes with drought tolerance and biomass accumulation.

Tomato MYB49 enhances resistance to Phytophthora infestans and tolerance to water deficit and salt stress

MYB49-overexpressing tomato plants showed significant resistance to Phytophthora infestans and tolerance to drought and salt stresses. This finding reveals the potential application of tomato MYB49 in future molecular breeding. Biotic and abiotic stresses severely reduce the productivity of tomato worldwide. Therefore, it is necessary to find key genes to simultaneously improve plant resistance to pathogens and tolerance to various abiotic stresses. In this study, based on homologous relationships with Arabidopsis R2R3-MYBs (AtMYBs) involved in responses to biotic and abiotic stresses, we identified a total of 24 R2R3-MYB transcription factors in the tomato genome. Among these tomato R2R3-MYBs, MYB49 (Solyc10g008700.1) was clustered into subgroup 11 by phylogenetic analysis, and its expression level was significantly induced after treatment with P. infestans, NaCl and PEG6000. Overexpression of MYB49 in tomato significantly enhanced the resistance of tomato to P. infestans, as evidenced by decreases in the number of necrotic cells, sizes of lesion, abundance of P. infestans, and disease index. Likewise, MYB49-overexpressing transgenic tomato plants also displayed increased tolerance to drought and salt stresses. Compared to WT plants, the accumulation of reactive oxygen species (ROS), malonaldehyde content, and relative electrolyte leakage was decreased, and peroxidase activity, superoxide dismutase activity, chlorophyll content, and photosynthetic rate were increased in MYB49-overexpressing tomato plants under P. infestans, salt or drought stress. These results suggested that tomato MYB49, as a positive regulator, could enhance the capacity to scavenge ROS, inhibit cell membrane damage and cell death, and protect chloroplasts, resulting in an improvement in resistance to P. infestans and tolerance to salt and drought stresses, and they provide a candidate gene for tomato breeding to enhance biotic stress resistance and abiotic stress tolerance.

Histone H3K9 demethylase JMJ25 epigenetically modulates anthocyanin biosynthesis in poplar

Anthocyanins are involved in several aspects of development and defence in poplar (Populus spp.). Although, over the past decades, significant progress has been made in uncovering these anthocyanin biosynthetic and regulatory mechanisms, the fundamental understanding of the epigenetic regulation in this pathway is still largely unclear. Here, we isolated a histone H3K9 demethylase gene JMJ25 from Populus and characterized its role in anthocyanin biosynthesis by genetic and biochemical approaches. JMJ25 was induced by continuous dark treatment. Overexpression of JMJ25 led to downregulated expression of anthocyanin biosynthetic genes in transgenic poplar, resulting in a significant reduction in anthocyanin content. ChIP-qPCR assays showed that JMJ25 could directly associate with MYB182 chromatin and dynamically demethylate at H3K9me2. Furthermore, JMJ25 also affected the DNA methylation levels of MYB182. By contrast, knockout of JMJ25 by CRISPR/Cas9 resulted in ectopic anthocyanin accumulation under dark condition and increased expression of anthocyanin biosynthetic genes. Our results support a model in which JMJ25 directly affects MYB182 expression by altering the histone methylation status of its chromatin and DNA methylation, resulting in repression of anthocyanin accumulation. This study uncovered an epigenetic mechanism that modulates anthocyanin biosynthesis in poplar.

SECONDARY WALL ASSOCIATED MYB1 is a positive regulator of secondary cell wall thickening in Brachypodium distachyon and is not found in the Brassicaceae

Grass biomass is comprised chiefly of secondary walls that surround fiber and xylem cells. A regulatory network of interacting transcription factors in part regulates cell wall thickening. We identified Brachypodium distachyon SECONDARY WALL ASSOCIATED MYB1 (SWAM1) as a potential regulator of secondary cell wall biosynthesis based on gene expression, phylogeny, and transgenic plant phenotypes. SWAM1 interacts with cellulose and lignin gene promoters with preferential binding to AC-rich sequence motifs commonly found in the promoters of cell wall-related genes. SWAM1 overexpression (SWAM-OE) lines had greater above-ground biomass with only a slight change in flowering time while SWAM1 dominant repressor (SWAM1-DR) plants were severely dwarfed with a striking reduction in lignin of sclerenchyma fibers and stem epidermal cell length. Cellulose, hemicellulose, and lignin genes were significantly down-regulated in SWAM1-DR plants and up-regulated in SWAM1-OE plants. There was no reduction in bioconversion yield in SWAM1-OE lines; however, it was significantly increased for SWAM1-DR samples. Phylogenetic and syntenic analyses strongly suggest that the SWAM1 clade was present in the last common ancestor between eudicots and grasses, but is not in the Brassicaceae. Collectively, these data suggest that SWAM1 is a transcriptional activator of secondary cell wall thickening and biomass accumulation in B. distachyon.

Genome-Wide Analysis of the NAC Transcription Factor Gene Family Reveals Differential Expression Patterns and Cold-Stress Responses in the Woody Plant Prunus mume

NAC transcription factors (TFs) participate in multiple biological processes, including biotic and abiotic stress responses, signal transduction and development. Cold stress can adversely impact plant growth and development, thereby limiting agricultural productivity. Prunus mume, an excellent horticultural crop, is widely cultivated in Asian countries. Its flower can tolerate freezing-stress in the early spring. To investigate the putative NAC genes responsible for cold-stress, we identified and analyzed 113 high-confidence PmNAC genes and characterized them by bioinformatics tools and expression profiles. These PmNACs were clustered into 14 sub-families and distributed on eight chromosomes and scaffolds, with the highest number located on chromosome 3. Duplicated events resulted in a large gene family; 15 and 8 pairs of PmNACs were the result of tandem and segmental duplicates, respectively. Moreover, three membrane-bound proteins (PmNAC59/66/73) and three miRNA-targeted genes (PmNAC40/41/83) were identified. Most PmNAC genes presented tissue-specific and time-specific expression patterns. Sixteen PmNACs (PmNAC11/19/20/23/41/48/58/74/75/76/78/79/85/86/103/111) exhibited down-regulation during flower bud opening and are, therefore, putative candidates for dormancy and cold-tolerance. Seventeen genes (PmNAC11/12/17/21/29/42/30/48/59/66/73/75/85/86/93/99/111) were highly expressed in stem during winter and are putative candidates for freezing resistance. The cold-stress response pattern of 15 putative PmNACs was observed under 4 °C at different treatment times. The expression of 10 genes (PmNAC11/20/23/40/42/48/57/60/66/86) was upregulated, while 5 genes (PmNAC59/61/82/85/107) were significantly inhibited. The putative candidates, thus identified, have the potential for breeding the cold-tolerant horticultural plants. This study increases our understanding of functions of the NAC gene family in cold tolerance, thereby potentially intensifying the molecular breeding programs of woody plants.

Regulation of Lignin Biosynthesis and Its Role in Growth-Defense Tradeoffs

Plant growth-defense tradeoffs are fundamental for optimizing plant performance and fitness in a changing biotic/abiotic environment. This process is thought to involve readjusting resource allocation to different pathways. It has been frequently observed that among secondary cell wall components, alteration in lignin biosynthesis results in changes in both growth and defense. How this process is regulated, leading to growth or defense, remains largely elusive. In this article, we review the canonical lignin biosynthesis pathway, the recently discovered tyrosine shortcut pathway, and the biosynthesis of unconventional C-lignin. We summarize the current model of the hierarchical transcriptional regulation of lignin biosynthesis. Moreover, the interface between recently identified transcription factors and the hierarchical model are also discussed. We propose the existence of a transcriptional co-regulation mechanism coordinating energy allowance among growth, defense and lignin biosynthesis.

AgriSeqDB: an online RNA-Seq database for functional studies of agriculturally relevant plant species

BACKGROUND

The genome-wide expression profile of genes in different tissues/cell types and developmental stages is a vital component of many functional genomic studies. Transcriptome data obtained by RNA-sequencing (RNA-Seq) is often deposited in public databases that are made available via data portals. Data visualization is one of the first steps in assessment and hypothesis generation. However, these databases do not typically include visualization tools and establishing one is not trivial for users who are not computational experts. This, as well as the various formats in which data is commonly deposited, makes the processes of data access, sharing and utility more difficult. Our goal was to provide a simple and user-friendly repository that meets these needs for data-sets from major agricultural crops.

DESCRIPTION

AgriSeqDB ( https://expression.latrobe.edu.au/agriseqdb ) is a database for viewing, analysing and interpreting developmental and tissue/cell-specific transcriptome data from several species, including major agricultural crops such as wheat, rice, maize, barley and tomato. The disparate manner in which public transcriptome data is often warehoused and the challenge of visualizing raw data are both major hurdles to data reuse. The popular eFP browser does an excellent job of presenting transcriptome data in an easily interpretable view, but previous implementation has been mostly on a case-by-case basis. Here we present an integrated visualisation database of transcriptome data-sets from six species that did not previously have public-facing visualisations. We combine the eFP browser, for gene-by-gene investigation, with the Degust browser, which enables visualisation of all transcripts across multiple samples. The two visualisation interfaces launch from the same point, enabling users to easily switch between analysis modes. The tools allow users, even those without bioinformatics expertise, to mine into data-sets and understand the behaviour of transcripts of interest across samples and time. We have also incorporated an additional graphic download option to simplify incorporation into presentations or publications.

CONCLUSION

Powered by eFP and Degust browsers, AgriSeqDB is a quick and easy-to-use platform for data analysis and visualization in five crops and Arabidopsis. Furthermore, it provides a tool that makes it easy for researchers to share their data-sets, promoting research collaborations and data-set reuse.

Comprehensive cross-genome survey and phylogeny of glycoside hydrolase family 16 members reveals the evolutionary origin of EG16 and XTH proteins in plant lineages

Carbohydrate-active enzymes (CAZymes) are central to the biosynthesis and modification of the plant cell wall. An ancient clade of bifunctional plant endo-glucanases (EG16 members) was recently revealed and proposed to represent a transitional group uniting plant xyloglucan endo-transglycosylase/hydrolase (XTH) gene products and bacterial mixed-linkage endo-glucanases in the phylogeny of glycoside hydrolase family 16 (GH16). To gain broader insights into the distribution and frequency of EG16 and other GH16 members in plants, the PHYTOZOME, PLAZA, NCBI and 1000 PLANTS databases were mined to build a comprehensive census among 1289 species, spanning the broad phylogenetic diversity of multiple algae through recent plant lineages. EG16, newly identified EG16-2 and XTH members appeared first in the green algae. Extant EG16 members represent the early adoption of the β-jellyroll protein scaffold from a bacterial or early-lineage eukaryotic GH16 gene, which is characterized by loop deletion and extension of the N terminus (in EG16-2 members) or C terminus (in XTH members). Maximum-likelihood phylogenetic analysis of EG16 and EG16-2 sequences are directly concordant with contemporary estimates of plant evolution. The lack of expansion of EG16 members into multi-gene families across green plants may point to a core metabolic role under tight control, in contrast to XTH genes that have undergone the extensive duplications typical of cell-wall CAZymes. The present census will underpin future studies to elucidate the physiological role of EG16 members across plant species, and serve as roadmap for delineating the closely related EG16 and XTH gene products in bioinformatic analyses of emerging genomes and transcriptomes.

Transcriptional regulation of osmotic stress tolerance in wheat (Triticum aestivum L.)

The current review provides an updated, new insights into the regulation of transcription mediated underlying mechanisms of wheat plants to osmotic stress perturbations. Osmotic stress tolerance mechanisms being complex are governed by multiple factors at physiological, biochemical and at the molecular level, hence approaches like "OMICS" that can underpin mechanisms behind osmotic tolerance in wheat is of paramount importance. The transcription factors (TFs) are a class of molecular proteins, which are involved in regulation, modulation and orchestrating the responses of plants to a variety of environmental stresses. Recent reports have provided novel insights on the role of TFs in osmotic stress tolerance via direct molecular links. However, our knowledge on the regulatory role TFs during osmotic stress tolerance in wheat remains limited. The present review in its first part sheds light on the importance of studying the role of osmotic stress tolerance in wheat plants and second aims to decipher molecular mechanisms of TFs belonging to several classes, including DREB, NAC, MYB, WRKY and bHLH, which have been reported to engage in osmotic stress mediated gene expression in wheat and third part covers the systems biology approaches to understand the transcriptional regulation of osmotic stress and the role of long non-coding RNAs in response to osmotic stress with special emphasis on wheat. The current concept may lead to an understanding in molecular regulation and signalling interaction of TFs under osmotic stress to clarify challenges and problems for devising potential strategies to improve complex regulatory events involved in plant tolerance to osmotic stress adaptive pathways in wheat.

A 5-Enolpyruvylshikimate 3-Phosphate Synthase Functions as a Transcriptional Repressor in Populus

Long-lived perennial plants, with distinctive habits of inter-annual growth, defense, and physiology, are of great economic and ecological importance. However, some biological mechanisms resulting from genome duplication and functional divergence of genes in these systems remain poorly studied. Here, we discovered an association between a poplar (Populus trichocarpa) 5-enolpyruvylshikimate 3-phosphate synthase gene (PtrEPSP) and lignin biosynthesis. Functional characterization of PtrEPSP revealed that this isoform possesses a helix-turn-helix motif in the N terminus and can function as a transcriptional repressor that regulates expression of genes in the phenylpropanoid pathway in addition to performing its canonical biosynthesis function in the shikimate pathway. We demonstrated that this isoform can localize in the nucleus and specifically binds to the promoter and represses the expression of a SLEEPER-like transcriptional regulator, which itself specifically binds to the promoter and represses the expression of PtrMYB021 (known as MYB46 in Arabidopsis thaliana), a master regulator of the phenylpropanoid pathway and lignin biosynthesis. Analyses of overexpression and RNAi lines targeting PtrEPSP confirmed the predicted changes in PtrMYB021 expression patterns. These results demonstrate that PtrEPSP in its regulatory form and PtrhAT form a transcriptional hierarchy regulating phenylpropanoid pathway and lignin biosynthesis in Populus.

Genome-wide identification and classification of MYB superfamily genes in peach

The MYB transcription factor superfamily is one of the largest superfamilies modulating various biological processes in plants. Over the past few decades, many MYB superfamily genes have been identified and characterized in some plant species. However, genes belonging to the MYB superfamily in peach (Prunus persica) have not been comprehensively identified and characterized although the genome sequences of peach were released several years ago. In total, this study yielded a set of 256 MYB superfamily genes that was divided into five subfamilies: the R2R3-MYB (2R-MYB), R1R2R3-MYB (3R-MYB), MYB-related (1R-MYB), 4R-MYB, and Atypical-MYB subfamilies. These subfamilies contained 128, 4, 109, 1, and 14 members, respectively. The 128 R2R3-MYB subfamily genes in peach were further clustered into 35 groups, and the 109 MYB-related subfamily genes were further clustered into 6 groups: the CCA1-like, CPC-like, TBP-like, I-box-binding-like, R-R-type, and Peach-specific groups. The motif compositions and exon/intron structures within each group within the R2R3-MYB or MYB-related subfamily in peach were highly conserved. The logo sequences of the R2 and R3 repeats of R2R3-MYB subfamily members were highly conserved with those in these repeats of several other plant species. Except for 48 novel peach-specific MYB genes, the remaining 208 out of 256 MYB genes in peach were conserved with the corresponding 198 MYB genes in A. thaliana. Additionally, the 256 MYB genes unevenly distributed on chromosomes 1 to 8 of the peach genome. Eighty-one orthologous pairs of peach/A. thaliana MYB genes were identified among 256 MYB genes in peach and 198 MYB genes in A. thaliana in this study. In addition, 146 pairs of paralogous MYB genes were identified on the eight chromosomes of peach. The expression levels of some of the 51 MYB genes selected for qRT-PCR analysis decreased or increased with red-fleshed fruit development, while the expression patterns of some genes followed no clear rules over the five developmental stages of fruits. This study laid the foundation for further functional analysis of MYB superfamily genes in peach and enriched the knowledge of MYB superfamily genes in plant species.

Genome-wide expression profiling of leaves and roots of watermelon in response to low nitrogen

BACKGROUND

Nitrogen (N) is a key macronutrient required for plant growth and development. In this study, watermelon plants were grown under hydroponic conditions at 0.2 mM N, 4.5 mM N, and 9 mM N for 14 days.

RESULTS

Dry weight and photosynthetic assimilation at low N (0.2 mM) was reduced by 29 and 74% compared with high N (9 mM). The photochemical activity (Fv/Fm) was also reduced from 0.78 at high N to 0.71 at low N. The N concentration in the leaf, stem, and root of watermelon under low N conditions was reduced by 68, 104, and 108%, respectively compared with 9 mM N treatment after 14 days of N treatment. In the leaf tissues of watermelon grown under low N conditions, 9598 genes were differentially expressed, out of which 4533 genes (47.22%) were up-regulated whereas, 5065 genes (52.78%) were down-regulated compared with high N. Similarly in the root tissues, 3956 genes were differentially expressed, out of which 1605 genes were up-regulated (40.57%) and 2351 genes were down-regulated (59.43%), compared with high N. Our results suggest that leaf tissues are more sensitive to N deficiency compared with root tissues. The gene ontology (GO) analysis showed that the availability of N significantly affected 19 biological processes, 8 cell component metabolic pathways, and 3 molecular functions in the leaves; and 13 biological processes, 12 molecular functions, and 5 cell component metabolic pathways in the roots of watermelon. The low affinity nitrate transporters, high affinity nitrate transporters, ammonium transporters, genes related with nitrogen assimilation, and chlorophyll and photosynthesis were expressed differentially in response to low N. Three nitrate transporters (Cla010066, Cla009721, Cla012765) substantially responded to low nitrate supply in the root and leaf tissues. Additionally, a large number of transcription factors (1365) were involved in adaptation to low N availability. The major transcription factor families identified in this study includes MYB, AP2-EREBP, bHLH, C2H2 and NAC.

CONCLUSION

Candidate genes identified in this study for nitrate uptake and transport can be targeted and utilized for further studies in watermelon breeding and improvement programs to improve N uptake and utilization efficiency.

A WD40 Repeat Protein from Camellia sinensis Regulates Anthocyanin and Proanthocyanidin Accumulation through the Formation of MYB⁻bHLH⁻WD40 Ternary Complexes

Flavan-3-ols and oligomeric proanthocyanidins (PAs) are the main nutritional polyphenols in green tea (Camellia sinensis), which provide numerous benefits to human health. To date, the regulatory mechanism of flavan-3-ol biosynthesis in green tea remains open to study. Herein, we report the characterization of a C. sinensis tryptophan-aspartic acid repeat protein (CsWD40) that interacts with myeloblastosis (MYB) and basic helix-loop-helix (bHLH) transcription factors (TFs) to regulate the biosynthesis of flavan-3-ols. Full length CsWD40 cDNA was cloned from leaves and was deduced to encode 342 amino acids. An in vitro yeast two-hybrid assay demonstrated that CsWD40 interacted with two bHLH TFs (CsGL3 and CsTT8) and two MYB TFs (CsAN2 and CsMYB5e). The overexpression of CsWD40 in Arabidopsis thaliana transparent testa glabra 1 (ttg1) restored normal trichome and seed coat development. Ectopic expression of CsWD40 alone in tobacco resulted in a significant increase in the anthocyanins of transgenic petals. CsWD40 was then coexpressed with CsMYB5e in tobacco plants to increase levels of both anthocyanins and PAs. Furthermore, gene expression analysis revealed that CsWD40 expression in tea plants could be induced by several abiotic stresses. Taken together, these data provide solid evidence that CsWD40 partners with bHLH and MYB TFs to form ternary WBM complexes to regulate anthocyanin, PA biosynthesis, and trichome development.

Poplar MYB transcription factor PtrMYB012 and its Arabidopsis AtGAMYB orthologs are differentially repressed by the Arabidopsis miR159 family

A phenotype-based screening of the T1 transgenic Arabidopsis population transformed by overexpression constructs of the entire poplar MYB transcription factor family found that overexpression of a poplar MYB transcription factor, PtrMYB012, in Arabidopsis resulted in upwardly curled rosette leaves, dwarfism and male sterility. Sequence analysis identified that PtrMYB012 is homologous to the Arabidopsis GAMYB genes (e.g., AtMYB65 and AtMYB33). Gene expression analysis revealed that PtrMYB012 is specifically expressed in floral tissues, especially in male catkins, similar to AtMYB65. It was well known that Arabidopsis GAMYBs are negatively regulated by microRNA159 (miR159) during vegetative growth; thus, the typical phenotypes of upwardly curled leaves, dwarfism and male sterility were only shown in overexpression of GAMYBs with mutations in the miR159 target sequence. To confirm our phenotypic consequences, we independently re-produced transgenic Arabidopsis plants overexpressing PtrMYB012 without mutations in the miR159 target sequence. The resulting 35 S::PtrMYB012 Arabidopsis plants phenocopied the previous transgenic Arabidopsis plants, suggesting that PtrMYB012 is probably not a target of Arabidopsis miR159 despite containing the conserved miR159 target sequence. To gain further insight, we produced transgenic poplars overexpressing the intact PtrMYB012. As a result, no conspicuous phenotype was found in 35 S::PtrMYB012 poplar plants. These results suggest that PtrMYB012 transcripts are down-regulated by miR159 in poplar but not in Arabidopsis. Indeed, subsequent 5'-RACE analysis confirmed that PtrMYB012 transcripts are completely degraded in poplar, probably by miR159, but not in Arabidopsis. These results suggest that species-specific family members of miR159 are important for the regulation of normal growth and development in plants.

A CsMYB6-CsTRY module regulates fruit trichome initiation in cucumber

Fruit epidermal features such as the number and size of trichomes or spines are important fruit quality traits in cucumber production. Little is known about the molecular mechanisms underlying fruit spine formation in cucumber. Here, we report functional characterization of the cucumber CsMYB6 gene, which encodes a MIXTA-like MYB transcription factor that plays an important role in regulating fruit trichome development. Spatial-temporal expression analyses revealed high-level expression of CsMYB6 in the epidermis of cucumber ovaries during fruit spine initiation, which was similar to the expression of CsTRY, a homolog of the Arabidopsis TRY gene that also plays a key role in trichome development. Overexpression of CsMYB6 and CsTRY in cucumber and Arabidopsis revealed that CsMYB6 and CsTRY act as negative regulators of trichome initiation in both species, and that CsMYB6 acted upstream of CsTRY in this process. CsMYB6 was found to bind to the three MYB binding sites inside the promoter region of CsTRY, and protein-protein interaction assays suggested that CsTRY also directly interacted with CsMYB6 protein. The results also revealed conserved and divergent roles of CsMYB6 and its Arabidopsis homolog AtMYB106 in trichome development. Collectively, our results reveal a novel mechanism in which the CsMYB6-CsTRY complex negatively regulates fruit trichome formation in cucumber.

A group of Populus trichocarpa DUF231 proteins exhibit differential O-acetyltransferase activities toward xylan

Wood represents the most abundant biomass produced by plants and one of its major components is acetyl xylan. Acetylation in xylan can occur at O-2 or O-3 of a xylosyl residue, at both O-2 and O-3 of a xylosyl residue, and at O-3 of a xylosyl residue substituted at O-2 with glucuronic acid. Acetyltransferases responsible for the regiospecific acetylation of xylan in tree species have not yet been characterized. Here we report the biochemical characterization of twelve Populus trichocarpa DUF231-containing proteins, named PtrXOATs, for their roles in the regiospecific acetylation of xylan. The PtrXOAT genes were found to be differentially expressed in Populus organs and among them, PtrXOAT1, PtrXOAT2, PtrXOAT9 and PtrXOAT10 exhibited the highest level of expression in stems undergoing wood formation. Activity assays of recombinant proteins demonstrated that all twelve PtrXOAT proteins were able to transfer acetyl groups from acetyl CoA onto a xylohexaose acceptor with PtrXOAT1, PtrXOAT2, PtrXOAT3, PtrXOAT11 and PtrXOAT12 having the highest activity. Structural analysis of the PtrXOAT-catalyzed reaction products using ¹H NMR spectroscopy revealed that PtrXOAT1, PtrXAOT2 and PtrXOAT3 mediated 2-O- and 3-O-monoacetylation and 2,3-di-O-acetylation of xylosyl residues and PtrXOAT11 and PtrXOAT12 only catalyzed 2-O- and 3-O-monoacetylation of xylosyl residues. Of the twelve PtrXOATs, only PtrXOAT9 and PtrXOAT10 were capable of transferring acetyl groups onto the O-3 position of 2-O-glucuronic acid-substituted xylosyl residues. Furthermore, when expressed in the Arabidopsis eskimo1 mutant, PtrXOAT1, PtrXAOT2 and PtrXOAT3 were able to rescue the defects in xylan acetylation. Together, these results demonstrate that the twelve PtrXOATs are acetyltransferases with different roles in xylan acetylation in P. trichocarpa.

Identification of MdDof genes in apple and analysis of their response to biotic or abiotic stress

As a classic plant-specific transcription factor family - the Dof domain proteins - are involved in a variety of biological processes in organisms ranging from unicellular Chlamydomonas to higher plants. However, there are limited reports of MdDof (Malus domestica Borkh. DNA-binding One Zinc Finger) domain proteins in fruit trees, especially in apple. In this study we identified 54 putative Dof transcription factors in the apple genome. We analysed the gene structures, protein motifs, and chromosome locations of each of the MdDof genes. Next, we characterised all 54 MdDofs their expression patterns under different abiotic and biotic stress conditions. It was found that MdDof6,26 not only played an important role in the biotic/abiotic stress but may also be involved in many molecular functions. Further, both in flower development and pollen tube growth it was found that the relative expression of MdDof24 increased rapidly, also with gene ontology analysis it was indicated that MdDof24 was involved in the chemical reaction and flower development pathways. Taken together, our results provide useful clues as to the function of MdDof genes in apple and serve as a reference for studies of Dof zinc finger genes in other plants.

Annotation of the Corymbia terpene synthase gene family shows broad conservation but dynamic evolution of physical clusters relative to Eucalyptus

Terpenes are economically and ecologically important phytochemicals. Their synthesis is controlled by the terpene synthase (TPS) gene family, which is highly diversified throughout the plant kingdom. The plant family Myrtaceae are characterised by especially high terpene concentrations, and considerable variation in terpene profiles. Many Myrtaceae are grown commercially for terpene products including the eucalypts Corymbia and Eucalyptus. Eucalyptus grandis has the largest TPS gene family of plants currently sequenced, which is largely conserved in the closely related E. globulus. However, the TPS gene family has been well studied only in these two eucalypt species. The recent assembly of two Corymbia citriodora subsp. variegata genomes presents an opportunity to examine the conservation of this important gene family across more divergent eucalypt lineages. Manual annotation of the TPS gene family in C. citriodora subsp. variegata revealed a similar overall number, and relative subfamily representation, to that previously reported in E. grandis and E. globulus. Many of the TPS genes were in physical clusters that varied considerably between Eucalyptus and Corymbia, with several instances of translocation, expansion/contraction and loss. Notably, there was greater conservation in the subfamilies involved in primary metabolism than those involved in secondary metabolism, likely reflecting different selective constraints. The variation in cluster size within subfamilies and the broad conservation between the eucalypts in the face of this variation are discussed, highlighting the potential contribution of selection, concerted evolution and stochastic processes. These findings provide the foundation to better understand terpene evolution within the ecologically and economically important Myrtaceae.

Reconstructing the suberin pathway in poplar by chemical and transcriptomic analysis of bark tissues

The tree bark periderm confers the first line of protection against pathogen invasion and abiotic stresses. The phellogen (cork cambium) externally produces cork (phellem) cells that are dead at maturity; while metabolically active, these tissues synthesize cell walls, as well as cell wall modifications, namely suberin and waxes. Suberin is a heteropolymer with aliphatic and aromatic domains, composed of acylglycerols, cross-linked polyphenolics and solvent-extractable waxes. Although suberin is essentially ubiquitous in vascular plants, the biochemical functions of many enzymes and the genetic regulation of its synthesis are poorly understood. We have studied suberin and wax composition in four developmental stages of hybrid poplar (Populus tremula x Populus alba) stem periderm. The amounts of extracellular ester-linked acyl lipids per unit area increased with tissue age, a trend not observed with waxes. We used RNA-Seq deep-sequencing technology to investigate the cork transcriptome at two developmental stages. The transcript analysis yielded 455 candidates for the biosynthesis and regulation of poplar suberin, including genes with proven functions in suberin metabolism, genes highlighted as candidates in other plant species and novel candidates. Among these, a gene encoding a putative lipase/acyltransferase of the GDSL-motif family emerged as a suberin polyester synthase candidate, and specific isoforms of peroxidase and laccase genes were preferentially expressed in cork, suggesting that their corresponding proteins may be involved in cross-linking aromatics to form lignin-like polyphenolics. Many transcriptional regulators with possible roles in meristem identity, cork differentiation and acyl-lipid metabolism were also identified. Our work provides the first large-scale transcriptomic dataset on the suberin-synthesizing tissue of poplar bark, contributing to our understanding of tree bark development at the molecular level. Based on these data, we have proposed a number of hypotheses that can be used in future research leading to novel biological insights into suberin biosynthesis and its physiological function.

The function of two type II metacaspases in woody tissues of Populus trees

Metacaspases (MCs) are cysteine proteases that are implicated in programmed cell death of plants. AtMC9 (Arabidopsis thaliana Metacaspase9) is a member of the Arabidopsis MC family that controls the rapid autolysis of the xylem vessel elements, but its downstream targets in xylem remain uncharacterized. PttMC13 and PttMC14 were identified as AtMC9 homologs in hybrid aspen (Populus tremula × tremuloides). A proteomic analysis was conducted in xylem tissues of transgenic hybrid aspen trees which carried either an overexpression or an RNA interference construct for PttMC13 and PttMC14. The proteomic analysis revealed modulation of levels of both previously known targets of metacaspases, such as Tudor staphylococcal nuclease, heat shock proteins and 14-3-3 proteins, as well as novel proteins, such as homologs of the PUTATIVE ASPARTIC PROTEASE3 (PASPA3) and the cysteine protease RD21 by PttMC13 and PttMC14. We identified here the pathways and processes that are modulated by PttMC13 and PttMC14 in xylem tissues. In particular, the results indicate involvement of PttMC13 and/or PttMC14 in downstream proteolytic processes and cell death of xylem elements. This work provides a valuable reference dataset on xylem-specific metacaspase functions for future functional and biochemical analyses.

Divergence and evolution of cotton bHLH proteins from diploid to allotetraploid

BACKGROUND

Polyploidy is considered a major driving force in genome expansion, yielding duplicated genes whose expression may be conserved or divergence as a consequence of polyploidization.

RESULTS

We compared the genome sequences of tetraploid cotton (Gossypium hirsutum) and its two diploid progenitors, G. arboreum and G. raimondii, and found that the bHLH genes were conserved over the polyploidization. Oppositely, the expression of the homeolgous gene pairs was diversified. The biased homeologous proportion for bHLH family is significantly higher (64.6%) than the genome wide homeologous expression bias (40%). Compared with cacao (T. cacao), orthologous genes only accounted for a small proportion (41.7%) of whole cotton bHLHs family. The further Ks analysis indicated that bHLH genes underwent at least two distinct episodes of whole genome duplication: a recent duplication (1.0-60.0 million years ago, MYA, 0.005 < Ks < 0.312) and an old duplication (> 60.0 MYA, 0.312 < Ks < 3.0). The old duplication event might have played a key role in the expansion of the bHLH family. Both recent and old duplicated pairs (68.8%) showed a divergent expression profile, indicating specialized functions. The expression diversification of the duplicated genes suggested it might be a universal feature of the long-term evolution of cotton.

CONCLUSIONS

Overview of cotton bHLH proteins indicated a conserved and divergent evolution from diploids to allotetraploid. Our results provided an excellent example for studying the long-term evolution of polyploidy.

The Bio-Analytic Resource for Plant Biology

Bioinformatic tools have become part of the way plant researchers undertake investigations. Large data sets encompassing genomes, transcriptomes, proteomes, epigenomes, and other "-omes" that have been generated in the past decade may be easily accessed with such tools, such that hypotheses may be generated at the click of a mouse. In this chapter, we'll cover the use of bioinformatic tools available at the Bio-Analytic Resource for Plant Biology at http://bar.utoronto.ca for exploring gene expression and coexpression patterns, undertaking promoter analyses, performing functional classification enrichment analyses for sets of genes, and examining protein-protein interactions. We also touch on some newer bioinformatic tools that allow integration of data from several sources for improved hypothesis generation, both for Arabidopsis and translationally. Most of the data sets come from Arabidopsis, but useful BAR tools for other species will be mentioned where appropriate.

GFDP: the gene family database in poplar

A gene family is formed by duplication of a single original gene. Poplar trees (genus Populus) are important, principally because of their ecological and economic benefits, and are one of the most widely distributed and adaptable trees in the world. Systematic identification and annotation of gene family members are primary steps in studying the function and evolution of poplar genomes. Here, we describe the construction of the Gene Family Database in Poplar (GFDP), which contains information that systematically describes 6551 genes distributed in 145 gene families. GFDP is designed to present important biological information, such as gene structure, protein length, isoelectric point and functional and evolutionary information, using highly visual displays. Data and graphs are visualized by a web-based interface. Users can browse and download data through all the major browsers. GFDP provides a comprehensive platform with a solid foundation for further study of poplar gene families. GFDP is free available.

Comprehensive Analysis of MYB Gene Family and Their Expressions Under Abiotic Stresses and Hormone Treatments in Tamarix hispida

The MYB transcription factors (TFs) is a plant TF families, which involves in hormone signal transduction, and abiotic stress tolerance, etc. However, there are few studies on the MYB TFs family and its regulatory mechanism in Tamarix hispida. In this study, 14 MYB genes (named ThMYB1 - ThMYB14) were cloned and characterized from T. hispida. The transcription profiles of ThMYBs in T. hispida under different abiotic stress conditions were monitored using qRT-PCR. Most of studied ThMYBs were significantly downregulated and/or upregulated by salt and osmotic stress, ABA, GA3 and JA treatments in at least one organ. Especially, ThMYB13 was induced in the leaves and roots of T. hispida when exposed to NaCl treatment at all study periods, indicating that it may involve in salt stress. To further study ThMYB13 function, ThMYB13 overexpression and knock-down plants and control plants transformed with an empty pROKII were obtained using a transient transformation system. Overexpression of ThMYB13 in T. hispida displayed the lowest O^2-, H₂O₂ and MDA accumulation, minimal cell death, the most stable K⁺/Na⁺ ratio and the lowest electrolyte leakage rate among the three kinds of transient expression in T. hispida. Conversely, the RNAi-silencing, transiently transformed plants displayed the opposite physiological changes. Therefore, ThMYB13 might play a role in salt stress tolerance in transgenic T. hispida plants.

PtoMYB170 positively regulates lignin deposition during wood formation in poplar and confers drought tolerance in transgenic Arabidopsis

Wood formation is a complex developmental process under multi-level transcriptional control executed by a large set of transcription factors. However, only limited members have been characterized to be key regulators of lignin biosynthesis in poplar. Here we report the conserved and unique functions of PtoMYB170, a transcription factor identified from Populus tomentosa (Chinese white poplar), in lignin deposition and drought tolerance in comparison with its duplicate paralog PtoMYB216. PtoMYB170 is preferentially expressed in young leaves and xylem tissues. Overexpression of PtoMYB170 in transgenic poplar plants resulted in stronger lignification and more thickened secondary wall in xylem compared with wild-type plants, whereas the CRISPR/Cas9-generated mutation of PtoMYB170 weakened lignin deposition, thereby leading to a more flexible and collapsed xylem phenotype. Transient expression experiments demonstrated that PtoMYB170 specifically activated the expression of lignin biosynthetic genes, consistent with the function of PtoMYB216. However, GUS staining assays revealed that PtoMYB170 was specifically expressed in guard cells of transgenic Arabidopsis while PtoMYB216 was not. Heterologous expression of PtoMYB170 in Arabidopsis enhanced stomatal closure in the dark and resulted in drought tolerance of the transgenic plants through reduced water loss, indicating a diversified role from PtoMYB216. These results revealed the PtoMYB170-dependent positive transcriptional regulation on lignin deposition in poplar and its coordinated function in enhancing drought tolerance by promoting dark-induced stomatal closure.

GmMYB58 and GmMYB205 are seed-specific activators for isoflavonoid biosynthesis in Glycine max

GmMYB58 and GmMYB205 are key positive regulators that are involved in isoflavonoid biosynthesis in seeds of Glycine max, and they activate the expression of several structural genes in the isoflavonoid pathway. MYB transcription factors (TFs) are major regulators involved in flavonoid/isoflavonoid biosynthesis in many plant species. However, functions of most MYB TFs remain unknown in flavonoid/isoflavonoid pathway in Glycine max. In this study, we identified 321 MYB TFs by genome-wide searching, and further isolated and functionally characterized two MYB TFs, GmMYB58 and GmMYB205. The deduced GmMYB58 and GmMYB205 proteins contain highly conserved R2R3 repeat domain at the N-terminal region that is the signature motif of R2R3-type MYB TFs. GmMYB58 and GmMYB205 were highly expressed in early seed development stages than in the other tested organs. GmMYB58 and GmMYB205 GFP fusion proteins were found to be localized in the nucleus when they were transiently expressed in Arabidopsis thaliana mesophyll protoplast. Both GmMYB58 and GmMYB205 can activate the promoter activities of GmCHS, GmIFS2, and GmHID in the transient trans-activation assays, and the activation of GmHID by both GmMYB58 and GmMYB205 was further confirmed by yeast one-hybrid assay. In addition, over-expression of GmMYB58 and GmMYB205 resulted in significant increases in expression levels of several pathway genes in soybean hairy roots, in particular, IFS2 by more than fivefolds in GmMYB205-over-expressing lines. Moreover, isoflavonoid contents were remarkably enhanced in the GmMYB58 and GmMYB205 over-expressing hairy roots than in the control. Our results suggest that GmMYB58 and GmMYB205 are seed-specific TFs, and they can enhance isoflavonoid biosynthesis mainly through the regulation of GmIFS2 and GmHID in G. max.

Multiple Copies of a Simple MYB-Binding Site Confers Trans-regulation by Specific Flavonoid-Related R2R3 MYBs in Diverse Species

In apple, the MYB transcription factor MYB10 controls the accumulation of anthocyanins. MYB10 is able to auto-activate its expression by binding its own promoter at a specific motif, the R1 motif. In some apple accessions a natural mutation, termed R6, has more copies of this motif within the MYB10 promoter resulting in stronger auto-activation and elevated anthocyanins. Here we show that other anthocyanin-related MYBs selected from apple, pear, strawberry, petunia, kiwifruit and Arabidopsis are able to activate promoters containing the R6 motif. To examine the specificity of this motif, members of the R2R3 MYB family were screened against a promoter harboring the R6 mutation. Only MYBs from subgroups 5 and 6 activate expression by binding the R6 motif, with these MYBs sharing conserved residues in their R2R3 DNA binding domains. Insertion of the apple R6 motif into orthologous promoters of MYB10 in pear (PcMYB10) and Arabidopsis (AtMY75) elevated anthocyanin levels. Introduction of the R6 motif into the promoter region of an anthocyanin biosynthetic enzyme F3'5'H of kiwifruit imparts regulation by MYB10. This results in elevated levels of delphinidin in both tobacco and kiwifruit. Finally, an R6 motif inserted into the promoter the vitamin C biosynthesis gene GDP-L-Gal phosphorylase increases vitamin C content in a MYB10-dependent manner. This motif therefore provides a tool to re-engineer novel MYB-regulated responses in plants.

Drought Tolerance Conferred in Soybean (Glycine max. L) by GmMYB84, a Novel R2R3-MYB Transcription Factor

MYB-type transcription factors (MYB TFs) play diverse roles in plant development and stress responses. However, the mechanisms underlying the actions of MYB TFs during stress response remain unclear. In this study we identified a R2R3-MYB TF in soybean (Glycine max), denoted GmMYB84, which contributes to drought resistance. Expression of GmMYB84 was induced by drought, salt stress, H2O2 and ABA. Compared with the wild type (WT), GmMYB84-overexpressing soybean mutants (OE lines) exhibited enhanced drought resistance with a higher survival rate, longer primary root length, greater proline and reactive oxygen species (ROS) contents, higher antioxidant enzyme activities [peroxidase (POD), catalase (CAT) and superoxide dismutase (SOD)], a lower dehydration rate and reduced malondialdehyde (MDA) content. We also found that ROS could induce SOD/POD/CAT activity in OE lines. In particular, we found that the optimal level of ROS is required for GmMYB84 to modulate primary root elongation. Some ROS-related genes were up-regulated under abiotic stress in GmMYB84 transgenic plants compared with the WT. Furthermore, electrophoretic mobility shift assay and luciferase reporter analysis demonstrated that GmMYB84 binds directly to the promoter of GmRBOHB-1 and GmRBOHB-2 genes. Based on this evidence, we propose a model for how GmMYB84, H2O2 and antioxidant enzymes work together to control root growth under both optimal and drought stress conditions.

Diurnal Cycling Transcription Factors of Pineapple Revealed by Genome-Wide Annotation and Global Transcriptomic Analysis

Circadian clock provides fitness advantage by coordinating internal metabolic and physiological processes to external cyclic environments. Core clock components exhibit daily rhythmic changes in gene expression, and the majority of them are transcription factors (TFs) and transcription coregulators (TCs). We annotated 1,398 TFs from 67 TF families and 80 TCs from 20 TC families in pineapple, and analyzed their tissue-specific and diurnal expression patterns. Approximately 42% of TFs and 45% of TCs displayed diel rhythmic expression, including 177 TF/TCs cycling only in the nonphotosynthetic leaf tissue, 247 cycling only in the photosynthetic leaf tissue, and 201 cycling in both. We identified 68 TF/TCs whose cycling expression was tightly coupled between the photosynthetic and nonphotosynthetic leaf tissues. These TF/TCs likely coordinate key biological processes in pineapple as we demonstrated that this group is enriched in homologous genes that form the core circadian clock in Arabidopsis and includes a STOP1 homolog. Two lines of evidence support the important role of the STOP1 homolog in regulating CAM photosynthesis in pineapple. First, STOP1 responds to acidic pH and regulates a malate channel in multiple plant species. Second, the cycling expression pattern of the pineapple STOP1 and the diurnal pattern of malate accumulation in pineapple leaf are correlated. We further examined duplicate-gene retention and loss in major known circadian genes and refined their evolutionary relationships between pineapple and other plants. Significant variations in duplicate-gene retention and loss were observed for most clock genes in both monocots and dicots.

Identification and expression analysis under abiotic stress of the R2R3-MYB genes in Ginkgo biloba L

The R2R3-MYB gene family is the largest MYB subfamily in plants and is involved in the regulation of plant secondary metabolism and specific morphogenesis, as well as the response to biotic and abiotic stress. However, a systematic identification and characterization of this gene family has not been carried out in Ginkgo biloba. In this study, we performed a transcriptome-wide survey from four tissues of G. biloba to determine the genetic variation and expression pattern of the R2R3-MYB genes. We analyzed 45 GbMYBs and identified 42 with a complete coding sequence via conserved motif searches. The MYB domain and other motifs in GbMYBs are highly conserved with Arabidopsis thaliana AtMYBs. Phylogenetic analysis of the GbMYBs and AtMYBs categorized the R2R3-MYBs into 26 subgroups, of which 11 subgroups included proteins from both G. biloba and Arabidopsis, and 1 subgroup was specific to G. biloba. Moreover, the GbMYBs expression patterns were analyzed in different tissues and abiotic stress conditions. The results revealed that GbMYBs were differentially expressed in various tissues and following abiotic stresses and phytohormone treatments, indicating their possible roles in biological processes and abiotic stress tolerance and adaptation. Our study demonstrated the functional diversity of the GbMYBs and will provide a foundation for future research into their biological and molecular functions.

The transcription factor MYB115 contributes to the regulation of proanthocyanidin biosynthesis and enhances fungal resistance in poplar

Proanthocyanidins (PAs) are major defense phenolic compounds in the leaves of poplar (Populus spp.) in response to abiotic and biotic stresses. Transcriptional regulation of PA biosynthetic genes by the MYB-basic helix-loop-helix (bHLH)-WD40 complexes in poplar is not still fully understood. Here, an Arabidopsis TT2-like gene MYB115 was isolated from Populus tomentosa and characterized by various molecular, genetic and biochemical approaches. MYB115 restored PA productions in the seed coat of the Arabidopsis tt2 mutant. Overexpression of MYB115 in poplar activated expression of PA biosynthetic genes, resulting in a significant increase in PA concentrations. By contrast, the CRISPR/Cas9-generated myb115 mutant exhibited reduced PA content and decreased expression of PA biosynthetic genes. MYB115 directly activated the promoters of PA-specific structural genes. MYB115 interacted with poplar TT8. Coexpression of MYB115, TT8 and poplar TTG1 significantly enhanced the expression of ANR1 and LAR3. Additionally, transgenic plants overexpressing MYB115 had increased resistance to the fungal pathogen Dothiorella gregaria, whereas myb115 mutant exhibited greater sensitivity compared with wild-type plants. Our data provide insight into the regulatory mechanisms controlling PA biosynthesis by MYB115 in poplar, which could be effectively employed for metabolic engineering of PAs to improve resistance to fungal pathogens.

Characterization and Functional Analysis of the Poplar Pectate Lyase-Like Gene PtPL1-18 Reveal Its Role in the Development of Vascular Tissues

Pectin is a major component of plant cell walls, and the structure of pectin impacts on the properties of wood. Although we know that pectate lyase (PL, EC 4.2.2.2) has a major influence on the structure of pectin, our knowledge of Pectate lyase-like genes (PLL) in tree species remains limited. To better understand the characteristics of PLL genes in trees and to identify novel PLL genes that are potentially involved in the development of wood, we performed comprehensive analyses of gene structures, phylogenetic relationships, chromosomal locations, gene duplication events, conserved protein motifs, and gene expression patterns of 30 PLLs in Populus trichocarpa (PtPL1s). We performed an in silico gene expression profiling and quantitative real-time PCR analysis and found that most of the PtPL1 genes from subgroups Ia and Ib were highly expressed in xylem. PtPL1-18 from subgroup Ia was preferentially expressed in developing primary xylem and in xylem cells that were developing secondary walls. Overexpression of PtPL1-18 in poplar reduced plant growth and xylem development. Reduced secondary cell wall thickening and irregular xylem cells were observed in the transgenic trees, probably due to their lower pectin content. Although pectin is not a major component of plant secondary cell walls, our results are consistent with the PtPL1 genes performing important functions during wood formation.

Proteomics analysis reveals the molecular mechanism underlying the transition from primary to secondary growth of poplar

Wood is the most important natural source of energy and also provides fuel and fiber. Considering the significant role of wood, it is critical to understand how wood is formed. Integration of knowledge about wood development at the cellular and molecular levels will allow more comprehensive understanding of this complex process. In the present study, we used a comparative proteomic approach to investigate the differences in protein profiles between primary and secondary growth in young poplar stems using tandem mass tag (TMT)-labeling. More than 10,816 proteins were identified, and, among these, 3106 proteins were differentially expressed during primary to secondary growth. Proteomic data were validated using a combination of histochemical staining, enzyme activity assays, and quantitative real-time PCR. Bioinformatics analysis revealed that these differentially expressed proteins are related to various metabolic pathways, mainly including signaling, phytohormones, cell cycle, cell wall, secondary metabolism, carbohydrate and energy metabolism, and protein metabolism as well as redox and stress pathways. This large proteomics dataset will be valuable for uncovering the molecular changes occurring during the transition from primary to secondary growth. Further, it provides new and accurate information for tree breeding to modify wood properties.

Poplar MYB115 and MYB134 Transcription Factors Regulate Proanthocyanidin Synthesis and Structure

The accumulation of proanthocyanidins is regulated by a complex of transcription factors composed of R2R3 MYB, basic helix-loop-helix, and WD40 proteins that activate the promoters of biosynthetic genes. In poplar (genus Populus), MYB134 is known to regulate proanthocyanidin biosynthesis by activating key flavonoid genes. Here, we characterize a second MYB regulator of proanthocyanidins, MYB115. Transgenic poplar overexpressing MYB115 showed a high-proanthocyanidin phenotype and reduced salicinoid accumulation, similar to the effects of MYB134 overexpression. Transcriptomic analysis of MYB115- and MYB134-overexpressing poplar plants identified a set of common up-regulated genes encoding proanthocyanidin biosynthetic enzymes and several novel uncharacterized MYB transcriptional repressors. Transient expression experiments demonstrated the capacity of both MYB134 and MYB115 to activate flavonoid promoters, but only in the presence of a basic helix-loop-helix cofactor. Yeast two-hybrid experiments confirmed the direct interaction of these transcription factors. The unexpected identification of dihydromyricetin in leaf extracts of both MYB115- and MYB134-overexpressing poplar led to the discovery of enhanced flavonoid B-ring hydroxylation and an increased proportion of prodelphinidins in proanthocyanidin of the transgenics. The dramatic hydroxylation phenotype of MYB115 overexpressors is likely due to the up-regulation of both flavonoid 3',5'-hydroxylases and cytochrome b5 Overall, this work provides new insight into the complexity of the gene regulatory network for proanthocyanidin synthesis in poplar.

A salt-stress-regulator from the Poplar R2R3 MYB family integrates the regulation of lateral root emergence and ABA signaling to mediate salt stress tolerance in Arabidopsis

The roles of most MYB transcription factors (TFs) in the poplar remain unclear. Here, we demonstrate that PtrSSR1, a salt-stress-regulator in the Populus trichocarpa R2R3 MYB gene family, mediates the tolerance of transgenic Arabidopsis plants to salt stress. The transcripts of PtrSSR1 could be induced by salt stress rapidly in poplar. Subcellular localization and yeast assays indicated that PtrSSR1 encoded a nuclear protein with transactivation activity. The Arabidopsis transformants overexpressing PtrSSR1 clearly displayed lateral root emergence (LRE) inhibition compared with wild-type (Wt) under normal conditions; while upon NaCl treatment, the transformants showed improved tolerance, and the LRs emerged faster from salt-induced inhibition. A strong correlation could exist between the LRE mediated by PtrSSR1 and abscisic acid (ABA), mainly because the transformants displayed more sensitivity to exogenous ABA during both seed germination and LRE, and had a distinctly increased level of endogenous ABA. Furthermore, several ABA- and salt-related genes, such as NCED3, ABI1 and CBL1, were up-regulated. Thus, our results suggest that elevation in the endogenous ABA content bring alteration of plant LR development, and that the poplar R2R3 MYB TF PtrSSR1 vitally improve salt stress tolerance by integrating the regulation of LRE and ABA signaling in Arabidopsis.

PtoMYB156 is involved in negative regulation of phenylpropanoid metabolism and secondary cell wall biosynthesis during wood formation in poplar

Some R2R3 MYB transcription factors have been shown to be major regulators of phenylpropanoid biosynthetic pathway and impact secondary wall formation in plants. In this study, we describe the functional characterization of PtoMYB156, encoding a R2R3-MYB transcription factor, from Populus tomentosa. Expression pattern analysis showed that PtoMYB156 is widely expressed in all tissues examined, but predominantly in leaves and developing wood cells. PtoMYB156 localized to the nucleus and acted as a transcriptional repressor. Overexpression of PtoMYB156 in poplar repressed phenylpropanoid biosynthetic genes, leading to a reduction in the amounts of total phenolic and flavonoid compounds. Transgenic plants overexpressing PtoMYB156 also displayed a dramatic decrease in secondary wall thicknesses of xylem fibers and the content of cellulose, lignin and xylose compared with wild-type plants. Transcript accumulation of secondary wall biosynthetic genes was down-regulated by PtoMYB156 overexpression. Transcriptional activation assays revealed that PtoMYB156 was able to repress the promoter activities of poplar CESA17, C4H2 and GT43B. By contrast, knockout of PtoMYB156 by CRISPR/Cas9 in poplar resulted in ectopic deposition of lignin, xylan and cellulose during secondary cell wall formation. Taken together, these results show that PtoMYB156 may repress phenylpropanoid biosynthesis and negatively regulate secondary cell wall formation in poplar.