Genome-wide transcriptome analysis of the transition from primary to secondary stem development in Populus trichocarpa

Global transcriptomics identification and analysis of transcriptional factors in different tissues of the paper mulberry

BACKGROUND

The paper mulberry (Broussonetia papyifera) is one of the multifunctional tree species in agroforestry system and is also commonly utilized in traditional medicine in China and other Asian countries. To identify the transcription factors (TFs) and comprehensively understand their regulatory roles in the growth of the paper mulberry, a global transcriptomics TF prediction and the differential expression analysis among root, shoot and leaf were performed by using RNA-seq.

RESULTS

Results indicate that there is 1, 337 TFs encoded by the paper mulberry and they belong to the 55 well-characterized TF families. Based on the phylogenetic analysis, the TFs exist extensively in all organisms are more conservative than those exclusively exist in plant and the paper mulberry has the closest relationship with the mulberry. According to the results of differential expression analysis, there are 627 TFs which exhibit the differential expression profiles in root, shoot and leaf. ARR-Bs, ARFs, NACs and bHLHs together with other root-specific and highly expressed TFs might account for the developed lateral root and unconspicuous taproot in the paper mulberry. Meanwhile, five TCPs highly expressed in shoot of the paper mulberry might negatively regulate the expression of 12 LBDs in shoot. Besides, LBDs, which could directly or indirectly cooperate with ARFs, bHLHs and NACs, seem to be the center knot involving in the regulation of the shoot development in the paper mulberry.

CONCLUSIONS

Our study provides the comprehensive transcriptomics identification of TFs in the paper mulberry without genome reference. A large number of lateral organ growth regulation related TFs exhibiting the tissue differential expression may entitle the paper mulberry the developed lateral roots, more branches and rapid growth. It will increase our knowledge of the structure and composition of TFs in tree plant and it will substantially contribute to the improving of this tree.

Tissue regeneration after bark girdling: an ideal research tool to investigate plant vascular development and regeneration

Regeneration is a common strategy for plants to survive the intrinsic and extrinsic challenges they face through their life cycle, and it may occur upon wounding. Bark girdling is applied to improve fruit production or harvest bark as medicinal material. When tree bark is removed, the cambium and phloem will be peeled off. After a small strip of bark is removed from trees, newly formed periderm and wound cambium develop from the callus on the surface of the trunk, and new phloem is subsequently derived from the wound cambium. However, after large-scale girdling, the newly formed sieve elements (SEs) appear earlier than the regenerated cambium, and both of them derive from differentiating xylem cells rather than from callus. This secondary vascular tissue regeneration mainly involves three key stages: callus formation and xylem cell dedifferentiation; SEs appearance and wound cambium formation. The new bark is formed within 1 month in poplar, Eucommia; thus, it provides high temporal resolution of regenerated tissues at different stages. In this review, we will illustrate the morphology, gene expression and phytohormone regulation of vascular tissue regeneration after large-scale girdling in trees, and also discuss the potential utilization of the bark girdling system in studies of plant vascular development and tissue regeneration.

Physiological and genomic basis of mechanical-functional trade-off in plant vasculature

Some areas in plant abiotic stress research are not frequently addressed by genomic and molecular tools. One such area is the cross reaction of gravitational force with upward capillary pull of water and the mechanical-functional trade-off in plant vasculature. Although frost, drought and flooding stress greatly impact these physiological processes and consequently plant performance, the genomic and molecular basis of such trade-off is only sporadically addressed and so is its adaptive value. Embolism resistance is an important multiple stress- opposition trait and do offer scopes for critical insight to unravel and modify the input of living cells in the process and their biotechnological intervention may be of great importance. Vascular plants employ different physiological strategies to cope with embolism and variation is observed across the kingdom. The genomic resources in this area have started to emerge and open up possibilities of synthesis, validation and utilization of the new knowledge-base. This review article assesses the research till date on this issue and discusses new possibilities for bridging physiology and genomics of a plant, and foresees its implementation in crop science.

Secondary cell wall composition and candidate gene expression in developing willow (Salix purpurea) stems

The properties of the secondary cell wall (SCW) in willow largely determine the suitability of willow biomass feedstock for potential bioenergy and biofuel applications. SCW development has been little studied in willow and it is not known how willow compares with model species, particularly the closely related genus Populus. To address this and relate SCW synthesis to candidate genes in willow, a tractable bud culture-derived system was developed in Salix purpurea, and cell wall composition and RNA-Seq transcriptome were followed in stems during early development. A large increase in SCW deposition in the period 0-2 weeks after transfer to soil was characterised by a big increase in xylan content, but no change in the frequency of substitution of xylan with glucuronic acid, and increased abundance of putative transcripts for synthesis of SCW cellulose, xylan and lignin. Histochemical staining and immunolabeling revealed that increased deposition of lignin and xylan was associated with xylem, xylem fibre cells and phloem fibre cells. Transcripts orthologous to those encoding xylan synthase components IRX9 and IRX10 and xylan glucuronyl transferase GUX1 in Arabidopsis were co-expressed, and showed the same spatial pattern of expression revealed by in situ hybridisation at four developmental stages, with abundant expression in proto-xylem, xylem fibre and ray parenchyma cells and some expression in phloem fibre cells. The results show a close similarity with SCW development in Populus species, but also give novel information on the relationship between spatial and temporal variation in xylan-related transcripts and xylan composition.

Overexpression of OsZHD1, a zinc finger homeodomain class homeobox transcription factor, induces abaxially curled and drooping leaf in rice

Leaf rolling is receiving considerable attention as an important agronomic trait in rice (Oryza sativa L.). However, little has been known on the molecular mechanism of rice leaf rolling, especially the abaxial rolling. We identified a novel abaxially curled and drooping leaf-dominant mutant from a T₁ transgenic rice line. The abaxially curled leaf phenotypes, co-segregating with the inserted transferred DNA, were caused by overexpression of a zinc finger homeodomain class homeobox transcription factor (OsZHD1). OsZHD1 exhibited a constitutive expression pattern in wild-type plants and accumulated in the developing leaves and panicles. Artificial overexpression of OsZHD1 or its closest homolog OsZHD2 induced the abaxial leaf curling. Histological analysis indicated that both the increased number and the abnormal arrangement of bulliform cells in leaf were responsible for the abaxially curled leaves. We herein reported OsZHD1 with key roles in rice morphogenesis, especially in the modulating of leaf rolling, which provided a novel insight into the molecular mechanism of leaf development in rice.

Transcriptome sequencing and profiling of expressed genes in cambial zone and differentiating xylem of Japanese cedar (Cryptomeria japonica)

BACKGROUND

Forest trees have ecological and economic importance, and Japanese cedar has highly valued wood attributes. Thus, studies of molecular aspects of wood formation offer practical information that may be used for screening and forward genetics approaches to improving wood quality.

RESULTS

After identifying expressed sequence tags in Japanese cedar tissue undergoing xylogenesis, we designed a custom cDNA microarray to compare expression of highly regulated genes throughout a growing season. This led to identification of candidate genes involved both in wood formation and later cessation of growth and dormancy. Based on homology to orthologous protein groups, the genes were assigned to functional classes. A high proportion of sequences fell into functional classes related to posttranscriptional modification and signal transduction, while transcription factors and genes involved in the metabolism of sugars, cell-wall synthesis and lignification, and cold hardiness were among other classes of genes identified as having a potential role in xylem formation and seasonal wood formation.

CONCLUSIONS

We obtained 55,051 unique sequences by next-generation sequencing of a cDNA library prepared from cambial meristem and derivative cells. Previous studies on conifers have identified unique sequences expressed in developing xylem, but this is the first comprehensive study utilizing a collection of expressed sequence tags for expression studies related to xylem formation in Japanese cedar, which belongs to a different lineage than the Pinaceae. Our characterization of these sequences should allow comparative studies of genome evolution and functional genetics of wood species.

Transcriptome sequencing and analysis of the fast growing shoots of moso bamboo (Phyllostachys edulis)

Background

The moso bamboo, a large woody bamboo with the highest ecological, economic, and cultural value of all bamboos, has one of the highest growth speeds in the world. Genetic research into moso bamboo has been scarce, partly because of the lack of previous genomic resources. In the present study, for the first time, we performed de novo transcriptome sequencing and mapped to the moso bamboo genomic resources (reference genome and genes) to produce a comprehensive dataset for the fast growing shoots of moso bamboo.

Results

The fast growing shoots mixed with six different heights and culms after leaf expansion of moso bamboo transcriptome were sequenced using the Illumina HiSeq™ 2000 sequencing platform, respectively. More than 80 million reads including 65,045,670 and 68,431,884 clean reads were produced in the two libraries. More than 81% of the reads were matched to the reference genome, and nearly 50% of the reads were matched to the reference genes. The genes with log 2 ratio > 2 or < −2 (P<0.001) were characterized as the most differentially expressed genes. 6,076 up-regulated and 4,613 down-regulated genes were classified into functional categories. Candidate genes which mainly involved transcript factors, plant hormones, cell cycle regulation, cell wall metabolism and cell morphogenesis genes were further analyzed and they may form a network that regulates the fast growth of moso bamboo shoots.

Conclusion

Firstly, our data provides the most comprehensive transcriptomic resource for moso bamboo to date. Candidate genes have been identified and they are potentially involved in the growth and development of moso bamboo. The results give a better insight into the mechanisms of moso bamboo shoots rapid growth and provide gene resources for improving plant growth.

Two poplar cellulose synthase-like D genes, PdCSLD5 and PdCSLD6, are functionally conserved with Arabidopsis CSLD3

Root hairs are tip-growing long tubular outgrowths of specialized epidermal cells, and are important for nutrient and water uptake and interaction with the soil microflora. Here we characterized two poplar cellulose synthase-like D (CSLD) genes, PdCSLD5 and PdCSLD6, the most probable orthologs to the Arabidopsis AtCSLD3/KOJAK gene. Both PdCSLD5 and PdCSLD6 are strongly expressed in roots, including in the root hairs. Subcellular localization experiments showed that these two proteins are located not only in the polarized plasma membrane of root hair tips, but also in Golgi apparatus of the root hair and non-hair-forming cells. Overexpression of these two poplar genes in the atcsld3 mutant was able to rescue most of the defects caused by disruption of AtCSLD3, including root hair morphological changes, altered cell wall monosaccharide composition, increased non-crystalline β-1,4-glucan and decreased crystalline cellulose contents. Taken together, our results provide evidence indicating that PdCSLD5 and PdCSLD6 are functionally conserved with AtCSLD3 and support a role for PdCSLD5 and PdCSL6 specifically in crystalline cellulose production in poplar root hair tips. The results presented here also suggest that at least part of the mechanism of root hair formation is conserved between herbaceous and woody plants.

Next-generation sequencing-based mRNA and microRNA expression profiling analysis revealed pathways involved in the rapid growth of developing culms in Moso bamboo

BACKGROUND

As one of the fastest-growing lignocellulose-abundant plants on Earth, bamboos can reach their final height quickly due to the expansion of individual internodes already present in the buds; however, the molecular processes underlying this phenomenon remain unclear. Moso bamboo (Phyllostachys heterocycla cv. Pubescens) internodes from four different developmental stages and three different internodes within the same stage were used in our study to investigate the molecular processes at the transcriptome and post-transcriptome level.

RESULTS

Our anatomical observations indicated the development of culms was dominated by cell division in the initial stages and by cell elongation in the middle and late stages. The four major endogenous hormones appeared to actively promote culm development. Using next-generation sequencing-based RNA-Seq, mRNA and microRNA expression profiling technology, we produced a transcriptome and post-transcriptome in possession of a large fraction of annotated Moso bamboo genes, and provided a molecular basis underlying the phenomenon of sequentially elongated internodes from the base to the top. Several key pathways such as environmental adaptation, signal transduction, translation, transport and many metabolisms were identified as involved in the rapid elongation of bamboo culms.

CONCLUSIONS

This is the first report on the temporal and spatial transcriptome and gene expression and microRNA profiling in a developing bamboo culms. In addition to gaining more insight into the unique growth characteristics of bamboo, we provide a good case study to analyze gene, microRNA expression and profiling of non-model plant species using high-throughput short-read sequencing. Also, we demonstrate that the integrated analysis of our multi-omics data, including transcriptome, post-transcriptome, proteome, yield more complete representations and additional biological insights, especially the complex dynamic processes occurring in Moso bamboo culms.

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.

Unexplored potentials of epigenetic mechanisms of plants and animals-theoretical considerations

Morphological and functional changes of cells are important for adapting to environmental changes and associated with continuous regulation of gene expressions. Genes are regulated–in part–by epigenetic mechanisms resulting in alternating patterns of gene expressions throughout life. Epigenetic changes responding to the environmental and intercellular signals can turn on/off specific genes, but do not modify the DNA sequence. Most epigenetic mechanisms are evolutionary conserved in eukaryotic organisms, and several homologs of epigenetic factors are present in plants and animals. Moreover, in vitro studies suggest that the plant cytoplasm is able to induce a nuclear reassembly of the animal cell, whereas others suggest that the ooplasm is able to induce condensation of plant chromatin. Here, we provide an overview of the main epigenetic mechanisms regulating gene expression and discuss fundamental epigenetic mechanisms and factors functioning in both plants and animals. Finally, we hypothesize that animal genome can be reprogrammed by epigenetic factors from the plant protoplast.

Genome-wide analysis of the expansin gene superfamily reveals grapevine-specific structural and functional characteristics

BACKGROUND

Expansins are proteins that loosen plant cell walls in a pH-dependent manner, probably by increasing the relative movement among polymers thus causing irreversible expansion. The expansin superfamily (EXP) comprises four distinct families: expansin A (EXPA), expansin B (EXPB), expansin-like A (EXLA) and expansin-like B (EXLB). There is experimental evidence that EXPA and EXPB proteins are required for cell expansion and developmental processes involving cell wall modification, whereas the exact functions of EXLA and EXLB remain unclear. The complete grapevine (Vitis vinifera) genome sequence has allowed the characterization of many gene families, but an exhaustive genome-wide analysis of expansin gene expression has not been attempted thus far.

METHODOLOGY/PRINCIPAL FINDINGS

We identified 29 EXP superfamily genes in the grapevine genome, representing all four EXP families. Members of the same EXP family shared the same exon-intron structure, and phylogenetic analysis confirmed a closer relationship between EXP genes from woody species, i.e. grapevine and poplar (Populus trichocarpa), compared to those from Arabidopsis thaliana and rice (Oryza sativa). We also identified grapevine-specific duplication events involving the EXLB family. Global gene expression analysis confirmed a strong correlation among EXP genes expressed in mature and green/vegetative samples, respectively, as reported for other gene families in the recently-published grapevine gene expression atlas. We also observed the specific co-expression of EXLB genes in woody organs, and the involvement of certain grapevine EXP genes in berry development and post-harvest withering.

CONCLUSION

Our comprehensive analysis of the grapevine EXP superfamily confirmed and extended current knowledge about the structural and functional characteristics of this gene family, and also identified properties that are currently unique to grapevine expansin genes. Our data provide a model for the functional characterization of grapevine gene families by combining phylogenetic analysis with global gene expression profiling.

Gene expression patterns underlying changes in xylem structure and function in response to increased nitrogen availability in hybrid poplar

Nitrogen availability has a strong influence on plant growth and development. In this study, we examined the effect of nitrogen availability on xylogenesis in hybrid poplar (Populus trichocarpa x deltoides H11-11). Saplings of hybrid poplar were fertilized for 33 d with either high or adequate levels of ammonium nitrate. We observed enhanced radial growth, wider vessels and fibres and thinner fibre walls in the secondary xylem of high N relative to adequate N plants. These anatomical differences translated into altered hydraulic properties with xylem being more transport efficient but also more vulnerable to drought-induced cavitation in high N plants. The changes in xylem structure and function were associated with differences in gene expression as revealed by the transcriptome analysis of the developing xylem region. We found 388 genes differentially expressed (fold change ±1.5, P-value ≤ 0.05), including a number of genes putatively involved in nitrogen and carbohydrate metabolism and various aspects of xylem cell differentiation. Several genes encoding known transcriptional regulators of secondary cell wall deposition were down-regulated in high N plants, corresponding with thinner secondary cell walls in these plants. The results of this study provide us with gene candidates potentially affecting xylem hydraulic and structural traits.

Association genetics of chemical wood properties in black poplar (Populus nigra)

Black poplar (Populus nigra) is a potential feedstock for cellulosic ethanol production, although breeding for this specific end use is required. Our goal was to identify associations between single nucleotide polymorphism (SNP) markers within candidate genes encoding cellulose and lignin biosynthetic enzymes, with chemical wood property phenotypic traits, toward the aim of developing genomics-based breeding technologies for bioethanol production. Pyrolysis molecular beam mass spectrometry was used to determine contents of five- and six-carbon sugars, lignin, and syringyl : guaiacyl ratio. The association population included 599 clones from 17 half-sib families, which were successfully genotyped using 433 SNPs from 39 candidate genes. Statistical analyses were performed to estimate genetic parameters, linkage disequilibrium (LD), and single marker and haplotype-based associations. A moderate to high heritability was observed for all traits. The LD, across all candidate genes, showed a rapid decay with physical distance. Analysis of single marker-phenotype associations identified six significant marker-trait pairs, whereas nearly 280 haplotypes were associated with phenotypic traits, in both an individual and multiple trait-specific manner. The rapid decay of LD within candidate genes in this population and the genetic associations identified suggest a close relationship between the associated SNPs and the causative polymorphisms underlying the genetic variation of lignocellulosic traits in black poplar.

De novo characterization of the Chinese fir (Cunninghamia lanceolata) transcriptome and analysis of candidate genes involved in cellulose and lignin biosynthesis

BACKGROUND

Chinese fir (Cunninghamia lanceolata) is an important timber species that accounts for 20-30% of the total commercial timber production in China. However, the available genomic information of Chinese fir is limited, and this severely encumbers functional genomic analysis and molecular breeding in Chinese fir. Recently, major advances in transcriptome sequencing have provided fast and cost-effective approaches to generate large expression datasets that have proven to be powerful tools to profile the transcriptomes of non-model organisms with undetermined genomes.

RESULTS

In this study, the transcriptomes of nine tissues from Chinese fir were analyzed using the Illumina HiSeq™ 2000 sequencing platform. Approximately 40 million paired-end reads were obtained, generating 3.62 gigabase pairs of sequencing data. These reads were assembled into 83,248 unique sequences (i.e. Unigenes) with an average length of 449 bp, amounting to 37.40 Mb. A total of 73,779 Unigenes were supported by more than 5 reads, 42,663 (57.83%) had homologs in the NCBI non-redundant and Swiss-Prot protein databases, corresponding to 27,224 unique protein entries. Of these Unigenes, 16,750 were assigned to Gene Ontology classes, and 14,877 were clustered into orthologous groups. A total of 21,689 (29.40%) were mapped to 119 pathways by BLAST comparison against the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The majority of the genes encoding the enzymes in the biosynthetic pathways of cellulose and lignin were identified in the Unigene dataset by targeted searches of their annotations. And a number of candidate Chinese fir genes in the two metabolic pathways were discovered firstly. Eighteen genes related to cellulose and lignin biosynthesis were cloned for experimental validating of transcriptome data. Overall 49 Unigenes, covering different regions of these selected genes, were found by alignment. Their expression patterns in different tissues were analyzed by qRT-PCR to explore their putative functions.

CONCLUSIONS

A substantial fraction of transcript sequences was obtained from the deep sequencing of Chinese fir. The assembled Unigene dataset was used to discover candidate genes of cellulose and lignin biosynthesis. This transcriptome dataset will provide a comprehensive sequence resource for molecular genetics research of C. lanceolata.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Expression and functions of proteases in vascular tissues

With the emergence of new models for wood formation and the increasing emphasis on improving the efficiency of cellulosic biofuel production, research on vascular tissue biology has intensified in recent years. Some of the most active areas of research focus on manipulating activity of enzymes in the cellulose, hemicellulose, pectin and lignin pathways. In addition, great strides have been made in the characterization of transcriptional networks controlling genes that affect differentiation, secondary cell wall synthesis and programmed cell death in xylem. Less attention has been devoted to the characterization of proteases that may be important regulators of post-translational events that affect vascular cell differentiation and function and cell wall composition. Several genes for proteases and components of the ubiquitin/26S proteasome pathway are upregulated in xylem and phloem and in cell culture systems for studying the differentiation of xylem tracheary elements (TEs). Although small molecule protease inhibitors have been used to explore the roles of proteases during the differentiation of cultured TEs, only a small number of vascular tissue-associated protease genes have been directly tested to determine whether they play roles in vascular tissue biology. In this report, we review roles for proteases in vascular cell differentiation and function as determined through the use of protease inhibitors and genetic analyses and conclude by identifying opportunities for future research in this area.

Natural hypolignification is associated with extensive oligolignol accumulation in flax stems

Flax (Linum usitatissimum) stems contain cells showing contrasting cell wall structure: lignified in inner stem xylem tissue and hypolignified in outer stem bast fibers. We hypothesized that stem hypolignification should be associated with extensive phenolic accumulation and used metabolomics and transcriptomics to characterize these two tissues. (1)H nuclear magnetic resonance clearly distinguished inner and outer stem tissues and identified different primary and secondary metabolites, including coniferin and p-coumaryl alcohol glucoside. Ultrahigh-performance liquid chromatography-Fourier transform ion cyclotron resonance-mass spectrometry aromatic profiling (lignomics) identified 81 phenolic compounds, of which 65 were identified, to our knowledge, for the first time in flax and 11 for the first time in higher plants. Both aglycone forms and glycosides of monolignols, lignin oligomers, and (neo)lignans were identified in both inner and outer stem tissues, with a preponderance of glycosides in the hypolignified outer stem, indicating the existence of a complex monolignol metabolism. The presence of coniferin-containing secondary metabolites suggested that coniferyl alcohol, in addition to being used in lignin and (neo)lignan formation, was also utilized in a third, partially uncharacterized metabolic pathway. Hypolignification of bast fibers in outer stem tissues was correlated with the low transcript abundance of monolignol biosynthetic genes, laccase genes, and certain peroxidase genes, suggesting that flax hypolignification is transcriptionally regulated. Transcripts of the key lignan genes Pinoresinol-Lariciresinol Reductase and Phenylcoumaran Benzylic Ether Reductase were also highly abundant in flax inner stem tissues. Expression profiling allowed the identification of NAC (NAM, ATAF1/2, CUC2) and MYB transcription factors that are likely involved in regulating both monolignol production and polymerization as well as (neo)lignan production.

The multiple fuzzy origins of woodiness within Balsaminaceae using an integrated approach. Where do we draw the line?

BACKGROUND AND AIMS

The family Balsaminaceae is essentially herbaceous, except for some woodier species that can be described as 'woody' herbs or small shrubs. The family is nested within the so-called balsaminoid clade of Ericales, including the exclusively woody families Tetrameristaceae and Marcgraviaceae, which is sister to the remaining families of the predominantly woody order. A molecular phylogeny of Balsaminaceae is compared with wood anatomical observations to find out whether the woodier species are derived from herbaceous taxa (i.e. secondary woodiness), or whether woodiness in the family represents the ancestral state for the order (i.e. primary woodiness).

METHODS

Wood anatomical observations of 68 Impatiens species and Hydrocera triflora, of which 47 are included in a multigene phylogeny, are carried out using light and scanning electron microscopy and compared with the molecular phylogenetic insights.

KEY RESULTS

There is much continuous variation in wood development between the Impatiens species studied, making the distinction between herbaceousness and woodiness difficult. However, the most woody species, unambiguously considered as truly woody shrubs, all display paedomorphic wood features pointing to secondary woodiness. This hypothesis is further supported by the molecular phylogeny, demonstrating that these most woody species are derived from herbaceous (or less woody) species in at least five independent clades. Wood formation in H. triflora is mostly confined to the ribs of the stems and shows paedomorphic wood features as well, suggesting that the common ancestor of Balsaminaceae was probably herbaceous.

CONCLUSIONS

The terms 'herbaceousness' and 'woodiness' are notoriously difficult to use in Balsaminaceae. However, anatomical observations and molecular sequence data show that the woodier species are derived from less woody or clearly herbaceous species, demonstrating that secondary woodiness has evolved in parallel.

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.

Dynamic DNA cytosine methylation in the Populus trichocarpa genome: tissue-level variation and relationship to gene expression

Background

DNA cytosine methylation is an epigenetic modification that has been implicated in many biological processes. However, large-scale epigenomic studies have been applied to very few plant species, and variability in methylation among specialized tissues and its relationship to gene expression is poorly understood.

Results

We surveyed DNA methylation from seven distinct tissue types (vegetative bud, male inflorescence [catkin], female catkin, leaf, root, xylem, phloem) in the reference tree species black cottonwood (Populus trichocarpa). Using 5-methyl-cytosine DNA immunoprecipitation followed by Illumina sequencing (MeDIP-seq), we mapped a total of 129,360,151 36- or 32-mer reads to the P. trichocarpa reference genome. We validated MeDIP-seq results by bisulfite sequencing, and compared methylation and gene expression using published microarray data. Qualitative DNA methylation differences among tissues were obvious on a chromosome scale. Methylated genes had lower expression than unmethylated genes, but genes with methylation in transcribed regions ("gene body methylation") had even lower expression than genes with promoter methylation. Promoter methylation was more frequent than gene body methylation in all tissues except male catkins. Male catkins differed in demethylation of particular transposable element categories, in level of gene body methylation, and in expression range of genes with methylated transcribed regions. Tissue-specific gene expression patterns were correlated with both gene body and promoter methylation.

Conclusions

We found striking differences among tissues in methylation, which were apparent at the chromosomal scale and when genes and transposable elements were examined. In contrast to other studies in plants, gene body methylation had a more repressive effect on transcription than promoter methylation.

Contrasting patterns of evolution following whole genome versus tandem duplication events in Populus

Comparative analysis of multiple angiosperm genomes has implicated gene duplication in the expansion and diversification of many gene families. However, empirical data and theory suggest that whole-genome and small-scale duplication events differ with respect to the types of genes preserved as duplicate pairs. We compared gene duplicates resulting from a recent whole genome duplication to a set of tandemly duplicated genes in the model forest tree Populus trichocarpa. We used a combination of microarray expression analyses of a diverse set of tissues and functional annotation to assess factors related to the preservation of duplicate genes of both types. Whole genome duplicates are 700 bp longer and are expressed in 20% more tissues than tandem duplicates. Furthermore, certain functional categories are over-represented in each class of duplicates. In particular, disease resistance genes and receptor-like kinases commonly occur in tandem but are significantly under-retained following whole genome duplication, while whole genome duplicate pairs are enriched for members of signal transduction cascades and transcription factors. The shape of the distribution of expression divergence for duplicated pairs suggests that nearly half of the whole genome duplicates have diverged in expression by a random degeneration process. The remaining pairs have more conserved gene expression than expected by chance, consistent with a role for selection under the constraints of gene balance. We hypothesize that duplicate gene preservation in Populus is driven by a combination of subfunctionalization of duplicate pairs and purifying selection favoring retention of genes encoding proteins with large numbers of interactions.

Analyses of the oligopeptide transporter gene family in poplar and grape

Background

Oligopeptide transporters (OPTs) are a group of membrane-localized proteins that have a broad range of substrate transport capabilities and that are thought to contribute to many biological processes. The OPT proteins belong to a small gene family in plants, which includes about 25 members in Arabidopsis and rice. However, no comprehensive study incorporating phylogeny, chromosomal location, gene structure, expression profiling, functional divergence and selective pressure analysis has been reported thus far for Populus and Vitis.

Results

In the present study, a comprehensive analysis of the OPT gene family in Populus (P. trichocarpa) and Vitis (V. vinifera) was performed. A total of 20 and 18 full-length OPT genes have been identified in Populus and Vitis, respectively. Phylogenetic analyses indicate that these OPT genes consist of two classes that can be further subdivided into 11 groups. Gene structures are considerably conserved among the groups. The distribution of OPT genes was found to be non-random across chromosomes. A high proportion of the genes are preferentially clustered, indicating that tandem duplications may have contributed significantly to the expansion of the OPT gene family. Expression patterns based on our analyses of microarray data suggest that many OPT genes may be important in stress response and functional development of plants. Further analyses of functional divergence and adaptive evolution show that, while purifying selection may have been the main force driving the evolution of the OPTs, some of critical sites responsible for the functional divergence may have been under positive selection.

Conclusions

Overall, the data obtained from our investigation contribute to a better understanding of the complexity of the Populus and Vitis OPT gene family and of the function and evolution of the OPT gene family in higher plants.

Genome-wide comparison of two poplar genotypes with different growth rates

The ecologically dominant and economically important genus Populus, with its available full genome sequence, has become an ideal woody species for genomic study. Rapid growth is one of the primary advantageous features of Populus, and extensive physiological research has been carried out on the growth of Populus throughout the growing period among different clones. However, the molecular information related to the mechanisms of rapid growth is rather limited. In this study, an Affymetrix poplar genome array was employed to analyze the transcriptomic changes from the pre-growth to the fast-growth phase in two poplar clones (P.deltoides × P.nigra, DN2, and P.nigra × (P.deltoides × P. nigra), NE19) with different growth rates. A total of 1,695 differently expressed genes were identified between two time points in NE19 and DN2 (two-way ANOVA, P < 0.01 and fold change ≥2). Except for genes changing in common for both clones, many transcripts were regulated specifically in one genotype. After functional analysis of the differentially expressed genes, distinct biological strategies seemed to be utilized by the two genotypes to accommodate their fast-growth phase. The faster-growing clone NE19, which has a higher photosynthetic rate and larger total leaf area, emphasized growth-related primary metabolism. However, the slower-growing DN2 tended to have more up-regulated genes involved in defense-related secondary metabolism and stress response. Emphasis of such divergent biological processes in two clones may explain their significant growth differences during the fast-growth phase.

Transcriptomic analysis of grain amaranth (Amaranthus hypochondriacus) using 454 pyrosequencing: comparison with A. tuberculatus, expression profiling in stems and in response to biotic and abiotic stress

BACKGROUND

Amaranthus hypochondriacus, a grain amaranth, is a C4 plant noted by its ability to tolerate stressful conditions and produce highly nutritious seeds. These possess an optimal amino acid balance and constitute a rich source of health-promoting peptides. Although several recent studies, mostly involving subtractive hybridization strategies, have contributed to increase the relatively low number of grain amaranth expressed sequence tags (ESTs), transcriptomic information of this species remains limited, particularly regarding tissue-specific and biotic stress-related genes. Thus, a large scale transcriptome analysis was performed to generate stem- and (a)biotic stress-responsive gene expression profiles in grain amaranth.

RESULTS

A total of 2,700,168 raw reads were obtained from six 454 pyrosequencing runs, which were assembled into 21,207 high quality sequences (20,408 isotigs + 799 contigs). The average sequence length was 1,064 bp and 930 bp for isotigs and contigs, respectively. Only 5,113 singletons were recovered after quality control. Contigs/isotigs were further incorporated into 15,667 isogroups. All unique sequences were queried against the nr, TAIR, UniRef100, UniRef50 and Amaranthaceae EST databases for annotation. Functional GO annotation was performed with all contigs/isotigs that produced significant hits with the TAIR database. Only 8,260 sequences were found to be homologous when the transcriptomes of A. tuberculatus and A. hypochondriacus were compared, most of which were associated with basic house-keeping processes. Digital expression analysis identified 1,971 differentially expressed genes in response to at least one of four stress treatments tested. These included several multiple-stress-inducible genes that could represent potential candidates for use in the engineering of stress-resistant plants. The transcriptomic data generated from pigmented stems shared similarity with findings reported in developing stems of Arabidopsis and black cottonwood (Populus trichocarpa).

CONCLUSIONS

This study represents the first large-scale transcriptomic analysis of A. hypochondriacus, considered to be a highly nutritious and stress-tolerant crop. Numerous genes were found to be induced in response to (a)biotic stress, many of which could further the understanding of the mechanisms that contribute to multiple stress-resistance in plants, a trait that has potential biotechnological applications in agriculture.

Identification of genes involved in cell wall biogenesis in grasses by differential gene expression profiling of elongating and non-elongating maize internodes

Despite the economic importance of grasses as food, feed, and energy crops, little is known about the genes that control their cell wall synthesis, assembly, and remodelling. Here a detailed transcriptome analysis that allowed the identification of genes involved in grass cell wall biogenesis is provided. Differential gene expression profiling, using maize oligonucleotide arrays, was used to identify genes differentially expressed between an elongating internode, containing cells exhibiting primary cell wall synthesis, and an internode that had just ceased elongation and in which many cells were depositing secondary cell wall material. This is one of only a few studies specifically aimed at the identification of cell wall-related genes in grasses. Analysis identified new candidate genes for a role in primary and secondary cell wall biogenesis in grasses. The results suggest that many proteins involved in cell wall processes during normal development are also recruited during defence-related cell wall remodelling events. This work provides a platform for studies in which candidate genes will be functionally tested for involvement in cell wall-related processes, increasing our knowledge of cell wall biogenesis and its regulation in grasses. Since several grasses are currently being developed as lignocellulosic feedstocks for biofuel production, this improved understanding of grass cell wall biogenesis is timely, as it will facilitate the manipulation of traits favourable for sustainable food and biofuel production.

Characterization of the plasma membrane proteins and receptor-like kinases associated with secondary vascular differentiation in poplar

The constituents of plasma membrane proteins, particularly the integral membrane proteins, are closely associated with the differentiation of plant cells. Secondary vascular differentiation, which gives rise to the increase in plant stem diameter, is the key process by which the volume of the plant body grows. However, little is known about the plasma membrane proteins that specifically function in the vascular differentiation process. Proteomic analysis of the membrane proteins in poplar differentiating secondary vascular tissues led to the identification 226 integral proteins in differentiating xylem and phloem tissues. A majority of the integral proteins identified were receptors (55 proteins), transporters (34 proteins), cell wall formation related (27 proteins) or intracellular trafficking (17 proteins) proteins. Gene expression analysis in developing vascular cells further demonstrated that cambium differentiation involves the expression of a group of receptor kinases which mediate an array of signaling pathways during secondary vascular differentiation. This paper provides an outline of the protein composition of the plasma membrane in differentiating secondary vascular tissues and sheds light on the role of receptor kinases during secondary vascular development.

Sequence information on simple sequence repeats and single nucleotide polymorphisms through transcriptome analysis of mungbean

Mungbean (Vigna radiata (L.) Wilczek) is a unique species in its ability to fix atmospheric nitrogen, with early maturity, and relatively good drought resistance. We used 454 sequencing technology for transcriptome sequencing. A total of 150 159 and 142 993 reads produced 5 254 and 6 374 large contigs (≥ 500 bp) with an average length of 833 and 853 for Sunhwa and Jangan, respectively. Functional annotation to known sequences yielded 41.34% and 41.74% unigenes for Jangan and Sunhwa. A higher number of simple sequence repeat (SSR) motifs was identified in Jangan (1 630) compared with that of Sunhwa (1 334). A similar SSR distribution pattern was observed in both varieties. A total of 8 249 single nucleotide polymorphisms (SNPs) and indels with 2 098 high-confidence candidates were identified in the two mungbean varieties. The average distance between individual SNPs was approximately 860 bp. Our report demonstrates the utility of transcriptomic data for implementing a functional annotation and development of genetic markers. We also provide large resource sequence data for mungbean improvement programs.

De novo assembled expressed gene catalog of a fast-growing Eucalyptus tree produced by Illumina mRNA-Seq

Background

De novo assembly of transcript sequences produced by short-read DNA sequencing technologies offers a rapid approach to obtain expressed gene catalogs for non-model organisms. A draft genome sequence will be produced in 2010 for a Eucalyptus tree species (E. grandis) representing the most important hardwood fibre crop in the world. Genome annotation of this valuable woody plant and genetic dissection of its superior growth and productivity will be greatly facilitated by the availability of a comprehensive collection of expressed gene sequences from multiple tissues and organs.

Results

We present an extensive expressed gene catalog for a commercially grown E. grandis × E. urophylla hybrid clone constructed using only Illumina mRNA-Seq technology and de novo assembly. A total of 18,894 transcript-derived contigs, a large proportion of which represent full-length protein coding genes were assembled and annotated. Analysis of assembly quality, length and diversity show that this dataset represent the most comprehensive expressed gene catalog for any Eucalyptus tree. mRNA-Seq analysis furthermore allowed digital expression profiling of all of the assembled transcripts across diverse xylogenic and non-xylogenic tissues, which is invaluable for ascribing putative gene functions.

Conclusions

De novo assembly of Illumina mRNA-Seq reads is an efficient approach for transcriptome sequencing and profiling in Eucalyptus and other non-model organisms. The transcriptome resource (Eucspresso, http://eucspresso.bi.up.ac.za/) generated by this study will be of value for genomic analysis of woody biomass production in Eucalyptus and for comparative genomic analysis of growth and development in woody and herbaceous plants.

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

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

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

Wood has a wide variety of uses and is arguably the most important renewable raw material. The composition of xylem cell types in wood determines the utility of different types of wood for distinct commercial applications. Using expression profiling and phylogenetic analysis, we identified many xylem-associated regulatory genes that may control the differentiation of cells involved in wood formation in Arabidopsis and poplar. Prominent among these are NAC domain transcription factors (NACs). We studied NACs with putative involvement as negative (XND1 from Arabidopsis and its poplar orthologs PopNAC118, PopNAC122, PopNAC128, PopNAC129), or positive (SND2 and SND3 from Arabidopsis and their poplar orthologs PopNAC105, PopNAC154, PopNAC156, PopNAC157) regulators of secondary cell wall synthesis. Using quantitative PCR and in situ hybridization, we evaluated expression of these Populus NACs in a developmental gradient and in association with reaction wood and found that representatives from both groups were associated with wood-forming tissue and phloem fibers. Additionally, XND1 orthologs were expressed in mesophyll cells of developing leaves. We prepared transgenic Arabidopsis and poplar plants for overexpression of selected NACs. XND1 overexpression in poplar resulted in severe stunting. Additionally, poplar XND1 overexpressors lacked phloem fibers and showed reductions in cell size and number, vessel number, and frequency of rays in the xylem. Overexpression of PopNAC122, an XND1 ortholog, yielded an analogous phenotype in Arabidopsis. Overexpression of PopNAC154 in poplar reduced height growth and increased the relative proportion of bark versus xylem.