The Populus Class III HD ZIP transcription factor POPCORONA affects cell differentiation during secondary growth of woody stems

Class III HD-ZIPs govern vascular cell fate: an HD view on patterning and differentiation

Plant vasculature is required for the transport of water and solutes throughout the plant body. It is constituted of xylem, specialized for transport of water, and phloem, that transports photosynthates. These two differentiated tissues are specified early in development and arise from divisions in the procambium, which is the vascular meristem during primary growth. During secondary growth, the xylem and phloem are further expanded via differentiation of cells derived from divisions in the cambium. Almost all of the developmental fate decisions in this process, including vascular specification, patterning, and differentiation, are regulated by transcription factors belonging to the class III homeodomain-leucine zipper (HD-ZIP III) family. This review draws together the literature describing the roles that these genes play in vascular development, looking at how HD-ZIP IIIs are regulated, and how they in turn influence other regulators of vascular development. Themes covered vary, from interactions between HD-ZIP IIIs and auxin, cytokinin, and brassinosteroids, to the requirement for exquisite spatial and temporal regulation of HD-ZIP III expression through miRNA-mediated post-transcriptional regulation, and interactions with other transcription factors. The literature described places the HD-ZIP III family at the centre of a complex network required for initiating and maintaining plant vascular tissues.

Agronomic performance of Populus deltoides trees engineered for biofuel production

BACKGROUND

One of the major barriers to the development of lignocellulosic feedstocks is the recalcitrance of plant cell walls to deconstruction and saccharification. Recalcitrance can be reduced by targeting genes involved in cell wall biosynthesis, but this can have unintended consequences that compromise the agronomic performance of the trees under field conditions. Here we report the results of a field trial of fourteen distinct transgenic Populus deltoides lines that had previously demonstrated reduced recalcitrance without yield penalties under greenhouse conditions.

RESULTS

Survival and productivity of the trial were excellent in the first year, and there was little evidence for reduced performance of the transgenic lines with modified target gene expression. Surprisingly, the most striking phenotypic effects in this trial were for two empty-vector control lines that had modified bud set and bud flush. This is most likely due to somaclonal variation or insertional mutagenesis. Traits related to yield, crown architecture, herbivory, pathogen response, and frost damage showed few significant differences between target gene transgenics and empty vector controls. However, there were a few interesting exceptions. Lines overexpressing the DUF231 gene, a putative O-acetyltransferase, showed early bud flush and marginally increased height growth. Lines overexpressing the DUF266 gene, a putative glycosyltransferase, had significantly decreased stem internode length and slightly higher volume index. Finally, lines overexpressing the PFD2 gene, a putative member of the prefoldin complex, had a slightly reduced volume index.

CONCLUSIONS

This field trial demonstrates that these cell wall modifications, which decreased cell wall recalcitrance under laboratory conditions, did not seriously compromise first-year performance in the field, despite substantial challenges, including an outbreak of a stem boring insect (Gypsonoma haimbachiana), attack by a leaf rust pathogen (Melampsora spp.), and a late frost event. This bodes well for the potential utility of these lines as advanced biofuels feedstocks.

MicroRNA and Transcription Factor: Key Players in Plant Regulatory Network

Recent achievements in plant microRNA (miRNA), a large class of small and non-coding RNAs, are very exciting. A wide array of techniques involving forward genetic, molecular cloning, bioinformatic analysis, and the latest technology, deep sequencing have greatly advanced miRNA discovery. A tiny miRNA sequence has the ability to target single/multiple mRNA targets. Most of the miRNA targets are transcription factors (TFs) which have paramount importance in regulating the plant growth and development. Various families of TFs, which have regulated a range of regulatory networks, may assist plants to grow under normal and stress environmental conditions. This present review focuses on the regulatory relationships between miRNAs and different families of TFs like; NF-Y, MYB, AP2, TCP, WRKY, NAC, GRF, and SPL. For instance NF-Y play important role during drought tolerance and flower development, MYB are involved in signal transduction and biosynthesis of secondary metabolites, AP2 regulate the floral development and nodule formation, TCP direct leaf development and growth hormones signaling. WRKY have known roles in multiple stress tolerances, NAC regulate lateral root formation, GRF are involved in root growth, flower, and seed development, and SPL regulate plant transition from juvenile to adult. We also studied the relation between miRNAs and TFs by consolidating the research findings from different plant species which will help plant scientists in understanding the mechanism of action and interaction between these regulators in the plant growth and development under normal and stress environmental conditions.

Vascular Morphodynamics During Secondary Growth

Quantification of vascular morphodynamics during secondary growth has been hampered by the scale of the process. Even in the tiny model plant Arabidopsis thaliana, the xylem can include more than 2000 cells in a single cross section, rendering manual counting impractical. Moreover, due to its deep location, xylem is an inaccessible tissue, limiting live imaging. A novel method to visualize and measure secondary growth progression has been proposed: "the Quantitative Histology" approach. This method is based on a detailed anatomical atlas, and image segmentation coupled with machine learning to automatically extract cell shapes and identify cell type. Here we present a new version of this approach, with a user-friendly interface implemented in the open source software LithoGraphX.

Overexpression of OsHox32 Results in Pleiotropic Effects on Plant Type Architecture and Leaf Development in Rice

BACKGROUND

The Class III homeodomain Leu zipper (HD-Zip III) gene family plays important roles in plant growth and development. Here, we analyze the function of OsHox32, an HD-Zip III family member, and show that it exhibits pleiotropic effects on regulating plant type architecture and leaf development in rice.

RESULTS

Transgenic lines overexpressing OsHox32 (OsHox32-OV) produce narrow leaves that roll towards the adaxial side. Histological analysis revealed a decreased number of bulliform cells in OsHox32-OV lines. In addition, the angle between the leaf and culm was reduced, resulting in an erect plant phenotype. The height of the plants was reduced, resulting in a semi-dwarf phenotype. In addition, the chlorophyll level was reduced, resulting in a decrease in the photosynthetic rate, but water use efficiency was significantly improved, presumably due to the rolled leaf phenotype. OsHox32 exhibited constitutive expression in different organs, with higher mRNA levels in the stem, leaf sheath, shoot apical meristems and young roots, suggesting a role in plant-type and leaf development. Moreover, OsHox32 mRNA levels were higher in light and lower in the dark under both long-day and short-day conditions, indicating that OsHox32 may be associated with light regulation. Photosynthesis-associated and chlorophyll biosynthesis-associated genes were down-regulated to result in the reduction of photosynthetic capacity in OsHox32-OV lines. mRNA level of six rice YABBY genes is up-regulated or down-regulated by OsHox32, suggesting that OsHox32 may regulate the architecture of plant type and leaf development by controlling the expression of YABBY genes in rice. In addition, OsHox32 mRNA level was induced by the phytohormones, indicating that OsHox32 may be involved in phytohormones regulatory pathways.

CONCLUSIONS

OsHox32, an HD-Zip III family member, plays pleiotropic effects on plant type architecture and leaf development in rice.

The SHORT-ROOT-like gene PtSHR2B is involved in Populus phellogen activity

SHORT-ROOT (SHR) is a GRAS transcription factor first characterized for its role in the specification of the stem cell niche and radial patterning in Arabidopsis thaliana (At) roots. Three SHR-like genes have been identified in Populus trichocarpa (Pt). PtSHR1 shares high similarity with AtSHR over the entire length of the coding sequence. The two other Populus SHR-like genes, PtSHR2A and PtSHR2B, are shorter in their 5' ends when compared with AtSHR. Unlike PtSHR1, that is expressed throughout the cambial zone of greenhouse-grown Populus trees, PtSHR2Bprom:uidA expression was detected in the phellogen. Additionally, PtSHR1 and PtSHR2B expression patterns markedly differ in the shoot apex and roots of in vitro plants. Transgenic hybrid aspen expressing PtSHR2B under the 35S constitutive promoter showed overall reduced tree growth while the proportion of bark increased relative to the wood. Reverse transcription-quantitative PCR (RT-qPCR) revealed increased transcript levels of cytokinin metabolism and response-related genes in the transgenic plants consistent with an increase of total cytokinin levels. This was confirmed by cytokinin quantification by LC-MS/MS. Our results indicate that PtSHR2B appears to function in the phellogen and therefore in the regulation of phellem and periderm formation, possibly acting through modulation of cytokinin homeostasis. Furthermore, this work points to a functional diversification of SHR after the divergence of the Populus and Arabidopsis lineages. This finding may contribute to selection and breeding strategies of cork oak in which, unlike Populus, the phellogen is active throughout the entire tree lifespan, being at the basis of a highly profitable cork industry.

Over-expression of AQUA1 in Populus alba Villafranca clone increases relative growth rate and water use efficiency, under Zn excess condition

Transgenic Populus alba over-expressing a TIP aquaporin ( aqua1) showed a higher growth rate under Zn excess, suggesting that aqua1 could be involved in water homeostasis, rather than in Zn homeostasis. Populus is the internationally accepted model for physiological and developmental studies of tree traits under stress. In plants, aquaporins facilitate and regulate the diffusion of water, however, few poplar aquaporins have been characterized to date. In this study, we reported for the first time an in vivo characterization of Populus alba clone Villafranca transgenic plants over-expressing a TIP aquaporin (aqua1) of P. x euramericana clone I-214. An AQUA1:GFP chimeric construct, over-expressed in P. alba Villafranca clones, shows a cytoplasmic localization in roots, and it localizes in guard cells in leaves. When over-expressed in transgenic plants, aqua1 confers a higher growth rate compared to wild-type (wt) plants, without affecting chlorophyll accumulation, relative water content (RWC), and fluorescence performances, but increasing the intrinsic Transpiration Efficiency. In response to Zn (1 mM), transgenic lines did not show a significant increase in Zn accumulation as compared to wt plants, even though the over-expression of this gene confers higher tolerance in root tissues. These results suggest that, in poplar plants, this gene could be principally involved in regulation of water homeostasis and biomass production, rather than in Zn homeostasis.

Identification and Comparative Analysis of microRNA in Wheat (Triticum aestivum L.) Callus Derived from Mature and Immature Embryos during In vitro Culture

Feasible and efficient tissue culture plays an important role in plant genetic engineering. Wheat (Triticum aestivum L.) immature embryos (IMEs) are preferred for tissue culture to mature embryos (MEs) because IMEs easily generate embryogenic callus, producing large number of plants. The molecular mechanisms of regulation and the biological pathways involved in embryogenic callus formation in wheat remain unclear. Here, microRNAs (miRNAs) potentially involved in embryogenic callus formation and somatic embryogenesis were identified through deep sequencing of small RNAs (sRNAs) and analyzed with bioinformatics tools. Six sRNA libraries derived from calli of IMEs and MEs after 3, 6, or 15 d of culture (DC) were constructed and sequenced. A total of 85 known miRNAs were identified, of which 30, 33, and 18 were differentially expressed (P < 0.05) between the IME and ME libraries at 3, 6, and 15 DC, respectively. Additionally, 171 novel and 41 candidate miRNAs were also identified, of the novel miRNA, 69, 67, and 37 were differentially expressed (P < 0.05) between the two types of libraries at 3, 6, and 15 DC, respectively. The expression patterns of eight known and eight novel miRNAs were validated using quantitative real-time polymerase chain reaction. Gene ontology annotation of differentially expressed miRNA targets provided information regarding the underlying molecular functions, biological processes, and cellular components involved in embryogenic callus development. Functional miRNAs, such as miR156, miR164, miR1432, miR398, and miR397, differentially expressed in IMEs and MEs might be related to embryogenic callus formation and somatic embryogenesis. This study suggests that miRNA plays an important role in embryogenic callus formation and somatic embryogenesis in wheat, and our data provide a useful resource for further research.

Xylem development - from the cradle to the grave

The development and growth of plants, as well as their successful adaptation to a variety of environments, is highly dependent on the conduction of water, nutrients and other important molecules throughout the plant body. Xylem is a specialized vascular tissue that serves as a conduit of water and minerals and provides mechanical support for upright growth. Wood, also known as secondary xylem, constitutes the major part of mature woody stems and roots. In the past two decades, a number of key factors including hormones, signal transducers and (post)transcriptional regulators have been shown to control xylem formation. We outline the main mechanisms shown to be essential for xylem development in various plant species, with an emphasis on Arabidopsis thaliana, as well as several tree species where xylem has a long history of investigation. We also summarize the processes which have been shown to be instrumental during xylem maturation. This includes mechanisms of cell wall formation and cell death which collectively complete xylem cell fate.

Association mapping in Populus reveals the interaction between Pto-miR530a and its target Pto-KNAT1

We used transcript profiling and multi-SNP association to investigate the genetic regulatory relationship between miRNA Pto-miR530a and its target Pto-KNAT1, identifying additive, dominant, and epistatic effects. MicroRNAs (miRNAs) play crucial roles in the post-transcriptional regulation of plant growth and development; indeed, many studies have described the importance of miRNA-target interactions in herbaceous species. However, elucidation of the miRNA-target interactions in trees may require novel strategies. In the present study, we describe a strategy combining expression profiling by reverse transcription quantitative PCR (RT-qPCR) and association mapping with multiple single nucleotide polymorphisms (SNPs) to evaluate the interaction between Pto-miR530a and its target Pto-KNAT1 in Populus tomentosa. RT-qPCR analysis showed a negative correlation (r = -0.62, P < 0.05) between expression levels of Pto-miR530a and Pto-KNAT1 in eight tissues. We used a Bayesian hierarchical model to identify allelic variants of Pto-miR530a and Pto-KNAT1 that associated with eight traits related to growth and wood properties, in a population of 460 unrelated individuals of P. tomentosa. This analysis identified 27 associations, with the proportions of phenotypic variance (R (2)) contributed by each SNP ranging of 0.82-15.81 %, the additive effects of each SNP ranging of 0.16-18.09, and the dominant effects ranging from -14.09 to 19.00. Epistatic interaction models showed a strong interaction among SNPs in the miRNA target with R (2) of 0.1-3.56 %, and information gain of significant SNP pairs of -3.09 to 0.93 %, representing the regulatory interactions between the miRNA and the mRNA. Thus, we used a new strategy that combines association genetics and expression profiling based on SNPs to study the regulatory relationship between this miRNA and its target mRNA, thereby providing novel advances in our understanding of the genetic architecture of important traits.

Molecular control of wood formation in trees

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

A resource for characterizing genome-wide binding and putative target genes of transcription factors expressed during secondary growth and wood formation in Populus

Identifying transcription factor target genes is essential for modeling the transcriptional networks underlying developmental processes. Here we report a chromatin immunoprecipitation sequencing (ChIP-seq) resource consisting of genome-wide binding regions and associated putative target genes for four Populus homeodomain transcription factors expressed during secondary growth and wood formation. Software code (programs and scripts) for processing the Populus ChIP-seq data are provided within a publically available iPlant image, including tools for ChIP-seq data quality control and evaluation adapted from the human Encyclopedia of DNA Elements (ENCODE) project. Basic information for each transcription factor (including members of Class I KNOX, Class III HD ZIP, BEL1-like families) binding are summarized, including the number and location of binding regions, distribution of binding regions relative to gene features, associated putative target genes, and enriched functional categories of putative target genes. These ChIP-seq data have been integrated within the Populus Genome Integrative Explorer (PopGenIE) where they can be analyzed using a variety of web-based tools. We present an example analysis that shows preferential binding of transcription factor ARBORKNOX1 to the nearest neighbor genes in a pre-calculated co-expression network module, and enrichment for meristem-related genes within this module including multiple orthologs of Arabidopsis KNOTTED-like Arabidopsis 2/6.

Functional plasticity of miR165/166 in plant development revealed by small tandem target mimic

MicroRNA 165 and 166 (miR165/166) is composed of nine members and targets five members (PHB, PHV, REV, ATHB8 and ATHB15) of the HD-ZIP III transcription factor family. Mutants generated by traditional methods could hardly reveal the overall functions of miR165/166 in plant development. In this study, the expressions of all miR165/166 members were simultaneously blocked by over-expressing STTM165/166-31 in Arabidopsis and tomato for functional dissection of miR165/166 family. Following a down-regulation of over 90% endogenous miR165/166, the target HD-ZIP III genes were correspondingly up-regulated in the STTM transgenic Arabidopsis and tomato plants. Notably, the STTM165/166-31 over-expressed Arabidopsis and tomato displayed pleiotropic effects on development which were not frequently observed in previously identified genetic mutants of either individual miR165/166 gene or any of the five target genes. Furthermore, the transgenic Arabidopsis showed increased IAA content and decreased IAA sensitivity accompanied by enhanced expressions of genes responsible for auxin biosynthesis and signaling, suggesting possible roles of auxin in mediation of miR165/166-regulated processes. Importantly, the transgenic Arabidopsis exhibited the improved behavior under salt stress. Overall, such diverse variations in plant development and physiological process revealed by STTM165/166 demonstrate a key role of miR165/166-mediated network in regulating plant development and responses to abiotic stresses.

Identification of differentially expressed transcripts associated with bast fibre development in Corchorus capsularis by suppression subtractive hybridization

The present study documented the predominant role of WRKY transcription factor in controlling genes of different pathways related to fibre formation in jute and could be a candidate gene for the improvement of jute fiber. Understanding of molecular mechanism associated with bast fibre development is of immense significance to achieve desired improvement in jute (Corchorus sp.). Therefore, suppression subtractive hybridization was successfully applied to identify genes involved in fibre developmental process in jute. The subtracted library of normal Corchorus capsularis as tester with respect to its fibre-deficient mutant as driver resulted in 2,685 expressed sequence tags which were assumed to represent the differentially expressed genes between two genotypes. The identified expressed sequence tags were assembled and clustered into 225 contigs and 231 singletons. Among these 456 unigenes, 377 were classified into 15 different functional categories while others were of unknown functional category. Reverse Northern analysis of the unigenes showed distinct variation in hybridization intensity of 11 transcripts between two genotypes tested. The findings were also documented by Northern and real-time PCR analysis. Varied expression level of these transcripts suggested their crucial involvement in fibre development in this species. Among these transcripts, WRKY transcription factor was documented to be a most important transcript which was in agreement with its known role in other plant species in possible regulation related to cell wall biosynthesis, expansion and lignification. This report constitutes first systematic analysis of genes involved in fibre development process in jute. It may be suggested that the information generated in this study would be useful for genetic improvement of fibre traits in this plant species.

The Populus ARBORKNOX1 homeodomain transcription factor regulates woody growth through binding to evolutionarily conserved target genes of diverse function

The class I KNOX homeodomain transcription factor ARBORKNOX1 (ARK1) is a key regulator of vascular cambium maintenance and cell differentiation in Populus. Currently, basic information is lacking concerning the distribution, functional characteristics, and evolution of ARK1 binding in the Populus genome. Here, we used chromatin immunoprecipitation sequencing (ChIP-seq) technology to identify ARK1 binding loci genome-wide in Populus. Computational analyses evaluated the distribution of ARK1 binding loci, the function of genes associated with bound loci, the effect of ARK1 binding on transcript levels, and evolutionary conservation of ARK1 binding loci. ARK1 binds to thousands of loci which are highly enriched proximal to the transcriptional start sites of genes of diverse functions. ARK1 target genes are significantly enriched in paralogs derived from the whole-genome salicoid duplication event. Both ARK1 and a maize (Zea mays) homolog, KNOTTED1, preferentially target evolutionarily conserved genes. However, only a small portion of ARK1 target genes are significantly differentially expressed in an ARK1 over-expression mutant. This study describes the functional characteristics and evolution of DNA binding by a transcription factor in an undomesticated tree, revealing complexities similar to those shown for transcription factors in model animal species.

The role of HD-ZIP III transcription factors and miR165/166 in vascular development and secondary cell wall formation

The Arabidopsis vascular system is composed of xylem and phloem, which form a well-defined collateral pattern in vascular bundles. Xylary element and fibers develop secondary cell walls (SCWs) that provide mechanical strength to support plant growth and to transport water and minerals to all above ground organs. SCWs also constitute the majority of terrestrial biomass for biofuel production. The biosynthesis of secondary cell walls are known to be under transcriptional regulation. Transcription factors, such as NAC (NAM, ATAF1/2 and CUC2) and MYB domain proteins, serve as master regulators in SCW development. Recent studies indicated that Class III homeodomain leucine zipper transcription factors (HD-ZIP III TFs) and microRNA 165/166 (miR165/166) may play important roles in SCW formation. Here we discuss the diverse functions of miR165/166 and HD-ZIPIII in vascular development and their interaction with the regulatory pathways of SCW biosynthesis.

Involvement of ethylene-responsive microRNAs and their targets in increased latex yield in the rubber tree in response to ethylene treatment

The rubber tree is an economically important plant that produces natural rubber for various industrial uses. The application of ethylene contributes to increased latex production in rubber trees; however, the molecular biology behind the effects of ethylene on latex yield remains to be elucidated. Recently, the intersection between microRNA (miRNA) regulation and phytohormone responses has been revealed. Insight into the regulation of miRNAs and their target genes should help to determine the functional importance of miRNAs as well as the role of miRNAs in signaling under ethylene stimulation in the rubber tree. In this study, hbr-miR159 and hbr-miR166 were down-regulated in bark under ethylene treatment. The ethylene also down-regulated ATHB15-like (Class III Homeodomain Leucine Zipper, HD-ZIP III) which have been extensively implicated in the regulation of primary and secondary vascular tissue pattern formation. The strong negative-regulation of ARF6/ARF8 caused by hbr-miR167 involved in an attenuation of vascular development and may gradually lead to bark dryness syndrome in the long term ethylene treatment. The negative correlation of hbr-miR172 and its target REF3 in the inner soft bark under ethylene treatment results in dramatic increases in latex yield in the ethylene-sensitive clone of the rubber tree. The overall results suggested that the differential expression of HD-ZIP III, miR167/ARF6, ARF8, and miR172/REF3 and related genes may play possible roles in the response to ethylene treatment, resulting in longer latex flow and increased latex yield.

Deep sequencing on a genome-wide scale reveals diverse stage-specific microRNAs in cambium during dormancy-release induced by chilling in poplar

BACKGROUND

Trees in temperate zones show periodicity by alternating active and dormant states to adapt to environmental conditions. Although phytohormones and transcriptional regulation were found to be involved in growth cessation and dormancy transition, little is known about the mechanisms of the dormancy-active growth transition, especially dormancy maintenance and release. Small RNAs are a group of short non-coding RNAs regulating gene expressions at the post-transcriptional level during plant development and the responses to environmental stress. No report on the expression profiling of small RNAs in the cambial meristem during the dormancy-active growth transition has been reported to date.

RESULTS

Three small RNA libraries from the cambium of poplar, representing endodormancy induced by short day conditions, ecodormancy induced by chilling and active growth induced by long day conditions, respectively, were generated and sequenced by Illumina high-throughput sequencing technology. This yielded 123 known microRNAs (miRNAs) with significant expression changes, which included developmental-, phytohormone- and stress-related miRNAs. Interestingly, miR156 and miR172 showed opposite expression patterns in the cambial dormancy-active growth transition. Additionally, miR160, which is involved in the auxin signaling pathway, was expressed specifically during endodormancy release by chilling. Furthermore, 275 novel miRNAs expressed in the cambial zone were identified, and 34 of them had high detection frequencies and unique expression patterns. Finally, the target genes of these novel miRNAs were predicted and some were validated experimentally by 5'RACE.

CONCLUSIONS

Our results provided a comprehensive analysis of small RNAs in the cambial meristem during dormancy-release at the genome-wide level and novel evidence of miRNAs involved in the regulation of this biological process.

Differential expression of microRNAs and their targets reveals a possible dual role in physiological bark disorder in rubber tree

Trunk phloem necrosis (TPN), a physiological bark disorder of the rubber tree (Hevea brasiliensis), is a serious problem that affects the yield of natural rubber. The resultant bark dryness occurs in up to half of a plantation's trees in almost every rubber tree plantation region, causing a great annual loss of dry rubber for natural rubber production. Different types of injury and physical damage caused by mechanical activation as well as environmental stresses cause physiological bark disorder in tree. Due to the essential role of miR166, miR393 and miR167 in vascular development and abiotic stress response in diverse plant species, it was interesting to investigate the role of these miRNAs in rubber trees, particularly during development of a physiological bark disorder. In this study, the expression pattern of miR166, miR393 and miR167; and their target genes, HD-ZIP III; TIR1 and ARF8, respectively; was demonstrated in healthy tree and different TPN trees. Their existence and function in vivo was validated using RNA ligase-mediated 5' rapid amplification of cDNA ends. Taken together, the results suggest a possible dual role of these three miRNAs in maintaining normal bark regeneration in healthy trees, coping with overtapping by affecting the wound healing system leading to abnormal bark regeneration in overtapped-TPN trees, and act as additional forces that enhance the attenuation of vascular development resulting in bark necrosis and cell death in the natural-TPN tree. This is the first study to address the molecular events of miRNAs involved in the physiological bark disorder TPN in rubber tree. Further study will open the possibility to better understanding of physiological and molecular perspectives during TPN development, and lead to improvement of monitoring the exploitation of rubber tree plantations.

Evaluation of experimental design and computational parameter choices affecting analyses of ChIP-seq and RNA-seq data in undomesticated poplar trees

Background

One of the great advantages of next generation sequencing is the ability to generate large genomic datasets for virtually all species, including non-model organisms. It should be possible, in turn, to apply advanced computational approaches to these datasets to develop models of biological processes. In a practical sense, working with non-model organisms presents unique challenges. In this paper we discuss some of these challenges for ChIP-seq and RNA-seq experiments using the undomesticated tree species of the genus Populus.

Results

We describe specific challenges associated with experimental design in Populus, including selection of optimal genotypes for different technical approaches and development of antibodies against Populus transcription factors. Execution of the experimental design included the generation and analysis of Chromatin immunoprecipitation-sequencing (ChIP-seq) data for RNA polymerase II and transcription factors involved in wood formation. We discuss criteria for analyzing the resulting datasets, determination of appropriate control sequencing libraries, evaluation of sequencing coverage needs, and optimization of parameters. We also describe the evaluation of ChIP-seq data from Populus, and discuss the comparison between ChIP-seq and RNA-seq data and biological interpretations of these comparisons.

Conclusions

These and other "lessons learned" highlight the challenges but also the potential insights to be gained from extending next generation sequencing-supported network analyses to undomesticated non-model species.

Small but thick enough--the Arabidopsis hypocotyl as a model to study secondary growth

The continuous production of vascular tissues through secondary growth results in radial thickening of plant organs and is pivotal for various aspects of plant growth and physiology, such as water transport capacity or resistance to mechanical stress. It is driven by the vascular cambium, which produces inward secondary xylem and outward secondary phloem. In the herbaceous plant Arabidopsis thaliana (Arabidopsis), secondary growth occurs in stems, in roots and in the hypocotyl. In the latter, radial growth is most prominent and not obscured by parallel ongoing elongation growth. Moreover, its progression is reminiscent of the secondary growth mode of tree trunks. Thus, the Arabidopsis hypocotyl is a very good model to study basic molecular mechanisms of secondary growth. Genetic approaches have succeeded in the identification of various factors, including peptides, receptors, transcription factors and hormones, which appear to participate in a complex network that controls radial growth. Many of these players are conserved between herbaceous and woody plants. In this review, we will focus on what is known about molecular mechanisms and regulators of vascular secondary growth in the Arabidopsis hypocotyl.

The formation of wood and its control

Wood continues to increase in importance as a sustainable source of energy and shelter. Wood formation is a dynamic process derived from plant secondary (radial) growth. Several experimental systems have been employed to study wood formation and its regulation. The use of genetic manipulation approaches and genome-wide analyses in model plants have significantly advanced our understanding of wood formation. In this review, we provide an update of our knowledge of the genetic and hormonal regulation of wood formation based on research in different plants systems, as well as considering the subject from an evo-devo perspective.

Transcriptome characteristics and six alternative expressed genes positively correlated with the phase transition of annual cambial activities in Chinese Fir (Cunninghamia lanceolata (Lamb.) Hook)

Background

The molecular mechanisms that govern cambial activity in angiosperms are well established, but little is known about these molecular mechanisms in gymnosperms. Chinese fir (Cunninghamia lanceolata (Lamb.) Hook), a diploid (2n  = 2x  = 22) gymnosperm, is one of the most important industrial and commercial timber species in China. Here, we performed transcriptome sequencing to identify the repertoire of genes expressed in cambium tissue of Chinese fir.

Methodology/Principal Findings

Based on previous studies, the four stage-specific cambial tissues of Chinese fir were defined using transmission electron microscopy (TEM). In total, 20 million sequencing reads (3.6 Gb) were obtained using Illumina sequencing from Chinese fir cambium tissue collected at active growth stage, with a mean length of 131 bp and a N50 of 90 bp. SOAPdenovo software was used to assemble 62,895 unigenes. These unigenes were further functionally annotated by comparing their sequences to public protein databases. Expression analysis revealed that the altered expression of six homologous genes (ClWOX1, ClWOX4, ClCLV1-like, ClCLV-like, ClCLE12, and ClPIN1-like) correlated positively with changes in cambial activities; moreover, these six genes might be directly involved in cambial function in Chinese fir. Further, the full-length cDNAs and DNAs for ClWOX1 and ClWOX4 were cloned and analyzed.

Conclusions

In this study, a large number of tissue/stage-specific unigene sequences were generated from the active growth stage of Chinese fir cambium. Transcriptome sequencing of Chinese fir not only provides extensive genetic resources for understanding the molecular mechanisms underlying cambial activities in Chinese fir, but also is expected to be an important foundation for future genetic studies of Chinese fir. This study indicates that ClWOX1 and ClWOX4 could be possible reverse genetic target genes for revealing the molecular mechanisms of cambial activities in Chinese fir.

The poplar MYB master switches bind to the SMRE site and activate the secondary wall biosynthetic program during wood formation

Wood is mainly composed of secondary walls, which constitute the most abundant stored carbon produced by vascular plants. Understanding the molecular mechanisms controlling secondary wall deposition during wood formation is not only an important issue in plant biology but also critical for providing molecular tools to custom-design wood composition suited for diverse end uses. Past molecular and genetic studies have revealed a transcriptional network encompassing a group of wood-associated NAC and MYB transcription factors that are involved in the regulation of the secondary wall biosynthetic program during wood formation in poplar trees. Here, we report the functional characterization of poplar orthologs of MYB46 and MYB83 that are known to be master switches of secondary wall biosynthesis in Arabidopsis. In addition to the two previously-described PtrMYB3 and PtrMYB20, two other MYBs, PtrMYB2 and PtrMYB21, were shown to be MYB46/MYB83 orthologs by complementation and overexpression studies in Arabidopsis. The functional roles of these PtrMYBs in regulating secondary wall biosynthesis were further demonstrated in transgenic poplar plants showing an ectopic deposition of secondary walls in PtrMYB overexpressors and a reduction of secondary wall thickening in their dominant repressors. Furthermore, PtrMYB2/3/20/21 together with two other tree MYBs, the Eucalyptus EgMYB2 and the pine PtMYB4, were shown to differentially bind to and activate the eight variants of the 7-bp SMRE consensus sequence, composed of ACC(A/T)A(A/C)(T/C). Together, our results indicate that the tree MYBs, PtrMYB2/3/20/21, EgMYB2 and PtMYB4, are master transcriptional switches that activate the SMRE sites in the promoters of target genes and thereby regulate secondary wall biosynthesis during wood formation.

PtrHB7, a class III HD-Zip gene, plays a critical role in regulation of vascular cambium differentiation in Populus

A key question in the secondary growth of trees is how differentiation of the vascular cambium cells is directed to concurrently form two different tissues: xylem or phloem. class III homeodomain-leucine zipper (HD-Zip III) genes are known to play critical roles in the initiation, patterning, and differentiation of the vascular system in the process of primary and secondary growth. However, the mechanism of how these genes control secondary vascular differentiation is unknown. Here, we show that a Populus class III HD-Zip gene, PtrHB7, was preferentially expressed in cambial zone. PtrHB7-suppressed plants displayed significant changes in vascular tissues with a reduction in xylem but increase in phloem. Transcriptional analysis revealed that genes regulating xylem differentiation were down-regulated, whereas genes regulating phloem differentiation were up-regulated. Correspondingly, PtrHB7 overexpression enhanced differentiation of cambial cells toward xylem cells but inhibited phloem differentiation. PtrHB7 regulation on cambial cell differentiation was associated with its transcript abundance. Together, the results demonstrated that PtrHB7 plays a critical role in controlling a balanced differentiation between secondary xylem and phloem tissues in the process of Populus secondary growth in a dosage-dependent manner.

Hormone interactions in xylem development: a matter of signals

Xylem provides long-distance transport of water and nutrients as well as structural support in plants. The development of the xylem tissues is modulated by several internal signals. In the last decades, the bloom of genetic and genomic tools has led to increased understanding of the molecular mechanisms underlying the function of the traditional plant hormones in xylem specification and differentiation. Critical functions have been assigned to novel signaling molecules, such as thermospermine. These signals do not function independently, but interact in a manner we are only now beginning to understand. We review the current knowledge of hormone signaling pathways and their crosstalk in cambial cell initiation and maintenance, and in xylem specification and differentiation.

DAMBE5: a comprehensive software package for data analysis in molecular biology and evolution

Since its first release in 2001 as mainly a software package for phylogenetic analysis, data analysis for molecular biology and evolution (DAMBE) has gained many new functions that may be classified into six categories: 1) sequence retrieval, editing, manipulation, and conversion among more than 20 standard sequence formats including MEGA, NEXUS, PHYLIP, GenBank, and the new NeXML format for interoperability, 2) motif characterization and discovery functions such as position weight matrix and Gibbs sampler, 3) descriptive genomic analysis tools with improved versions of codon adaptation index, effective number of codons, protein isoelectric point profiling, RNA and protein secondary structure prediction and calculation of minimum folding energy, and genomic skew plots with optimized window size, 4) molecular phylogenetics including sequence alignment, testing substitution saturation, distance-based, maximum parsimony, and maximum-likelihood methods for tree reconstructions, testing the molecular clock hypothesis with either a phylogeny or with relative-rate tests, dating gene duplication and speciation events, choosing the best-fit substitution models, and estimating rate heterogeneity over sites, 5) phylogeny-based comparative methods for continuous and discrete variables, and 6) graphic functions including secondary structure display, optimized skew plot, hydrophobicity plot, and many other plots of amino acid properties along a protein sequence, tree display and drawing by dragging nodes to each other, and visual searching of the maximum parsimony tree. DAMBE features a graphic, user-friendly, and intuitive interface and is freely available from http://dambe.bio.uottawa.ca (last accessed April 16, 2013).

Identification of molecular processes needed for vascular formation through transcriptome analysis of different vascular systems

Background

Vascular system formation has been studied through molecular and genetic approaches in Arabidopsis, a herbaceous dicot that is used as a model system. Different vascular systems have developed in other plants such as crops and trees. Uncovering shared mechanisms underlying vascular development by transcriptome analysis of different vascular systems may help to transfer knowledge acquired from Arabidopsis to other economically important species.

Results

Conserved vascular genes and biological processes fundamental to vascular development were explored across various plants. Through comparative transcriptome analysis, 226 genes from Arabidopsis, 217 genes from poplar and 281 genes from rice were identified as constituting 107 conserved vascular gene groups. These gene groups are expressed mainly in vascular tissues and form a complex coexpression network with multiple functional connections. To date, only half of the groups have been experimentally investigated. The conserved vascular gene groups were classified into 9 essential processes for vascular development. 18 groups (17%) lack of annotations were classified as having unknown functions.

Conclusion

The study provides a map of fundamental biological processes conserved across different vascular systems. It identifies gaps in the experimental investigation of pathways active in vascular formation, which if explored, could lead to a more complete understanding of vascular development.

KNOX1 genes regulate lignin deposition and composition in monocots and dicots

Plant secondary cell walls are deposited mostly in vascular tissues such as xylem vessels, tracheids, and fibers. These cell walls are composed of a complex matrix of compounds including cellulose, hemicellulose, and lignin. Lignin functions primarily to maintain the structural and mechanical integrity of both the transport vessel and the entire plant itself. Since lignin has been identified as a major source of biomass for biofuels, regulation of secondary cell wall biosynthesis has been a topic of much recent investigation. Biosynthesis and patterning of lignin involves many developmental and environmental cues including evolutionarily conserved transcriptional regulatory modules and hormonal signals. Here, we investigate the role of the class I Knotted1-like-homeobox (KNOX) genes and gibberellic acid in the lignin biosynthetic pathway in a representative monocot and a representative eudicot. Knotted1 overexpressing mutant plants showed a reduction in lignin content in both maize and tobacco. Expression of four key lignin biosynthesis genes was analyzed and revealed that KNOX1 genes regulate at least two steps in the lignin biosynthesis pathway. The negative regulation of lignin both in a monocot and a eudicot by the maize Kn1 gene suggests that lignin biosynthesis may be preserved across large phylogenetic distances. The evolutionary implications of regulation of lignification across divergent species are discussed.

Characterization of microRNAs expressed during secondary wall biosynthesis in Acacia mangium

MicroRNAs (miRNAs) play critical regulatory roles by acting as sequence specific guide during secondary wall formation in woody and non-woody species. Although thousands of plant miRNAs have been sequenced, there is no comprehensive view of miRNA mediated gene regulatory network to provide profound biological insights into the regulation of xylem development. Herein, we report the involvement of six highly conserved amg-miRNA families (amg-miR166, amg-miR172, amg-miR168, amg-miR159, amg-miR394, and amg-miR156) as the potential regulatory sequences of secondary cell wall biosynthesis. Within this highly conserved amg-miRNA family, only amg-miR166 exhibited strong differences in expression between phloem and xylem tissue. The functional characterization of amg-miR166 targets in various tissues revealed three groups of HD-ZIP III: ATHB8, ATHB15, and REVOLUTA which play pivotal roles in xylem development. Although these three groups vary in their functions, -psRNA target analysis indicated that miRNA target sequences of the nine different members of HD-ZIP III are always conserved. We found that precursor structures of amg-miR166 undergo exhaustive sequence variation even within members of the same family. Gene expression analysis showed three key lignin pathway genes: C4H, CAD, and CCoAOMT were upregulated in compression wood where a cascade of miRNAs was downregulated. This study offers a comprehensive analysis on the involvement of highly conserved miRNAs implicated in the secondary wall formation of woody plants.

MicroRNAs in trees

MicroRNAs (miRNAs) are 20-24 nucleotide long molecules processed from a specific class of RNA polymerase II transcripts that mainly regulate the stability of mRNAs containing a complementary sequence by targeted degradation in plants. Many features of tree biology are regulated by miRNAs affecting development, metabolism, adaptation and evolution. MiRNAs may be modified and harnessed for controlled suppression of specific genes to learn about gene function, or for practical applications through genetic engineering. Modified (artificial) miRNAs act as dominant suppressors and are particularly useful in tree genetics because they bypass the generations of inbreeding needed for fixation of recessive mutations. The purpose of this review is to summarize the current status of information on miRNAs in trees and to guide future studies on the role of miRNAs in the biology of woody perennials and to illustrate their utility in directed genetic modification of trees.

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.

From thin to thick: major transitions during stem development

The variability of shoot architecture in plants is striking and one of the most extreme examples of adaptive growth in higher organisms. Mediated by the differential activity of apical and lateral meristems, flexibility in stem growth essentially contributes to this variability. In spite of this importance, the regulation of major events in stem development is largely unexplored. Recently, however, novel approaches exploiting knowledge from root and leaf development are starting to shed light on molecular mechanisms that regulate this essential plant organ. In this review, we summarize our understanding of initial patterning events in stems, discuss prerequisites for the initiation of lateral stem growth and highlight the burning questions in this context.