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

Profiling of the gene expression and alternative splicing landscapes of Eucalyptus grandis

Eucalyptus is a widely planted hardwood tree species due to its fast growth, superior wood properties and adaptability. However, the post-transcriptional regulatory mechanisms controlling tissue development and stress responses in Eucalyptus remain poorly understood. In this study, we performed a comprehensive analysis of the gene expression profile and the alternative splicing (AS) landscape of E. grandis using strand-specific RNA-Seq, which encompassed 201 libraries including different organs, developmental stages, and environmental stresses. We identified 10 416 genes (33.49%) that underwent AS, and numerous differentially expressed and/or differential AS genes involved in critical biological processes, such as primary-to-secondary growth transition of stems, adventitious root formation, aging and responses to phosphorus- or boron-deficiency. Co-expression analysis of AS events and gene expression patterns highlighted the potential upstream regulatory role of AS events in multiple processes. Additionally, we highlighted the lignin biosynthetic pathway to showcase the potential regulatory functions of AS events in the KNAT3 and IRL3 genes within this pathway. Our high-quality expression atlas and AS landscape serve as valuable resources for unravelling the genetic control of woody plant development, long-term adaptation, and understanding transcriptional diversity in Eucalyptus. Researchers can conveniently access these resources through the interactive ePlant browser (https://bar.utoronto.ca/eplant_eucalyptus).

Transcriptome analysis of the transition from primary to secondary growth of vertical stem in Eucalyptus grandis

Eucalyptus was one of the most cultivated hardwood species worldwide, with rapid growth, good wood properties and a wide range of adaptability. Eucalyptus stem undergoes primary growth (longitudinal growth) followed by secondary growth (radial growth), which produces biomass that is an important source of energy worldwide. In order to better understand the genetic regulation of secondary growth in Eucalyptus grandis, Transcriptome analyses in stem segments along a developmental gradient from the third internode to the eleventh internode of E. grandis that spanned primary to secondary growth were carried out. 5,149 genes that were differentially expressed during stem development were identified. Combining the trend analysis by the Mfuzz method and the module-trait correlation analysis by the Weighted Gene Co-expression Network Analysis method, a total of 70 differentially expressed genes (DEGs) selected from 868 DEGs with high connectivity were found to be closely correlated with secondary growth. Results revealed that the differential expression of these DEGs suggests that they may involve in the primary growth or secondary growth. AP1, YAB2 TFs and EXP genes are highly expressed in the IN3, whereas NAC, MYB TFs are likely to be important for secondary growth. These results will expand our understanding of the complex molecular and cellular events of secondary growth and provide a foundation for future studies on wood formation in Eucalyptus.

Genes expression profiles in vascular cambium of Eucalyptus urophylla × Eucalyptus grandis at different ages

BACKGROUND

Wood is a secondary xylem generated by vascular cambium. Vascular cambium activities mainly include cambium proliferation and vascular tissue formation through secondary growth, thereby producing new secondary phloem inward and secondary xylem outward and leading to continuous tree thickening and wood formation. Wood formation is a complex biological process, which is strictly regulated by multiple genes. Therefore, molecular level research on the vascular cambium of different tree ages can lead to the identification of both key and related genes involved in wood formation and further explain the molecular regulation mechanism of wood formation.

RESULTS

In the present study, RNA-Seq and Pac-Bio Iso-Seq were used for profiling gene expression changes in Eucalyptus urophylla × Eucalyptus grandis (E. urograndis) vascular cambium at four different ages. A total of 59,770 non-redundant transcripts and 1892 differentially expressed genes (DEGs) were identified. The expression trends of the DEGs related to cell division and differentiation, cell wall biosynthesis, phytohormone, and transcription factors were analyzed. The DEGs encoding expansin, kinesin, cycline, PAL, GRP9, KNOX, C2C2-dof, REV, etc., were highly expressed in E. urograndis at three years old, leading to positive effects on growth and development. Moreover, some gene family members, such as NAC, MYB, HD-ZIP III, RPK, and RAP, play different regulatory roles in wood formation because of their sophisticated transcriptional network and function redundantly.

CONCLUSIONS

These candidate genes are a potential resource to further study wood formation, especially in fast-growing and adaptable eucalyptus. The results may also serve as a basis for further research to unravel the molecular mechanism underlying wood formation.

The PtoTCP20-miR396d-PtoGRF15 module regulates secondary vascular development in Populus

Secondary vascular development is a key biological characteristic of woody plants and the basis of wood formation. Our understanding of gene expression regulation and dynamic changes in microRNAs (miRNAs) during secondary vascular development is still limited. Here we present an integrated analysis of the miRNA and mRNA transcriptome of six phase-specific tissues-the shoot apex, procambium, primary vascular tissue, cambium, secondary phloem, and secondary xylem-in Populus tomentosa. Several novel regulatory modules, including the PtoTCP20-miR396d-PtoGRF15 module, were identified during secondary vascular development in Populus. A series of biochemical and molecular experiments confirmed that PtoTCP20 activated transcription of the miR396d precursor gene and that miR396d targeted PtoGRF15 to downregulate its expression. Plants overexpressing miR396d (35S:miR396d) showed enhanced secondary growth and increased xylem production. Conversely, during the transition from primary to secondary vascular development, plants with downregulated PtoTCP20expression (PtoTCP20-SRDX), downregulated miR396 expression (35S:STTM396), and PtoGRF15 overexpression (35S:PtoGRF15) showed delayed secondary growth. Novel regulatory modules were identified by integrated analysis of the miRNA and mRNA transcriptome, and the regulatory role of the PtoTCP20-miR396d-PtoGRF15 signaling cascade in secondary vascular development was validated in Populus, providing information to support improvements in forest cultivation and wood properties.

Overexpression of miR390b promotes stem elongation and height growth in Populus

MicroRNA390 (miR390) is involved in plant growth and development by down-regulating the expression of the downstream genes trans-acting short interfering RNA3 (TAS3) and AUXIN RESPONSE FACTORs (ARFs). There is a scarcity of research on the involvement of the miR390-TAS3-ARFs pathway in the stem development of Populus. Here, differentially expressed miRNAs during poplar stem development were screened by small RNA sequencing analysis, and a novel function of miR390b in stem development was revealed. Overexpression of miR390b (OE-miR390b) resulted in a large increase in the number of xylem fiber cells and a slight decrease in the cell length at the longitudinal axis. Overall increases in stem elongation and plant height were observed in the OE-miR390b plants. According to transcriptome sequencing results and transient co-expression analysis, TAS3.1 and TAS3.2 were identified as the target genes of miR390 in poplar and were negatively regulated by miR390 in the apex. The transcription levels of ARF3.2 and ARF4 were significantly repressed in OE-miR390b plants and strongly negatively correlated with the number of xylem fiber cells along the longitudinal axis. These findings indicate that the conserved miR390-TAS3-ARFs pathway in poplar is involved in stem elongation and plant height growth.

Genome-Wide Identification and Characterization of Auxin Response Factor (ARF) Gene Family Involved in Wood Formation and Response to Exogenous Hormone Treatment in Populus trichocarpa

Auxin is a key regulator that virtually controls almost every aspect of plant growth and development throughout its life cycle. As the major components of auxin signaling, auxin response factors (ARFs) play crucial roles in various processes of plant growth and development. In this study, a total of 35 PtrARF genes were identified, and their phylogenetic relationships, chromosomal locations, synteny relationships, exon/intron structures, cis-elements, conserved motifs, and protein characteristics were systemically investigated. We also analyzed the expression patterns of these PtrARF genes and revealed that 16 of them, including PtrARF1, 3, 7, 11, 13-17, 21, 23, 26, 27, 29, 31, and 33, were preferentially expressed in primary stems, while 15 of them, including PtrARF2, 4, 6, 9, 10, 12, 18-20, 22, 24, 25, 28, 32, and 35, participated in different phases of wood formation. In addition, some PtrARF genes, with at least one cis-element related to indole-3-acetic acid (IAA) or abscisic acid (ABA) response, responded differently to exogenous IAA and ABA treatment, respectively. Three PtrARF proteins, namely PtrARF18, PtrARF23, and PtrARF29, selected from three classes, were characterized, and only PtrARF18 was a transcriptional self-activator localized in the nucleus. Moreover, Y2H and bimolecular fluorescence complementation (BiFC) assay demonstrated that PtrARF23 interacted with PtrIAA10 and PtrIAA28 in the nucleus, while PtrARF29 interacted with PtrIAA28 in the nucleus. Our results provided comprehensive information regarding the PtrARF gene family, which will lay some foundation for future research about PtrARF genes in tree development and growth, especially the wood formation, in response to cellular signaling and environmental cues.

Single-nuclei transcriptome analysis of the shoot apex vascular system differentiation in Populus

Differentiation of stem cells in the plant apex gives rise to aerial tissues and organs. Presently, we lack a lineage map of the shoot apex cells in woody perennials - a crucial gap considering their role in determining primary and secondary growth. Here, we used single-nuclei RNA-sequencing to determine cell type-specific transcriptomes of the Populus vegetative shoot apex. We identified highly heterogeneous cell populations clustered into seven broad groups represented by 18 transcriptionally distinct cell clusters. Next, we established the developmental trajectories of the epidermis, leaf mesophyll and vascular tissue. Motivated by the high similarities between Populus and Arabidopsis cell population in the vegetative apex, we applied a pipeline for interspecific single-cell gene expression data integration. We contrasted the developmental trajectories of primary phloem and xylem formation in both species, establishing the first comparison of vascular development between a model annual herbaceous and a woody perennial plant species. Our results offer a valuable resource for investigating the principles underlying cell division and differentiation conserved between herbaceous and perennial species while also allowing us to examine species-specific differences at single-cell resolution.

PagGRF11 Overexpression Promotes Stem Development and Dwarfing in Populus

Poplar is not only an important woody plant, but also a model species for molecular plant studies. We identified PagGRF11 (pAxG07Gg0005700), a homolog of the Arabidopsis AtGRF1 (AT4G37740) and AtGRF2 (AT2G22840) gene. We transformed the poplar clone "84K" with PagGRF11, and the transgenic overexpressed plants (PagGRF11-OE) showed plant height reduction (dwarfing), stem diameter increase, internode shortening, and larger leaf area. The Arabidopsis overexpression line grf-oe (Overexpression of PagGRF11 in Arabidopsis), mutant line atgrf (a loss-of-function mutant of the AtGRF1 gene of Arabidopsis thaliana), and mutant trans-complementary line atgrf+oe (overexpression of PagGRF11 in mutant plants (atgrf)) also showed different leaf size phenotypes. Further, tissue sections revealed that increased xylem production was the main cause of stem thickening. Transcriptome differential expression analysis of PagGRF11 overexpressed and control plants showed that PagGRF11 promoted CCCH39(C3H39) expression. The expression profile of CCCH39 in different tissues showed that it was highly expressed in xylem. Yeast single hybrid and instantaneous double luciferase assay results showed that PagGRF11 directly transcribed and activated CCCH39 expression through interaction with cis-acting element GARE (TCTGTTG), thus promoting xylem development. This is the first finding that GRF positively regulates xylem development through CCCH39 expression activation and further suggests that PagGRF11 is a potential target for increasing wood yield.

PdeHCA2 affects biomass in Populus by regulating plant architecture, the transition from primary to secondary growth, and photosynthesis

PdeHCA2 regulates the transition from primary to secondary growth, plant architecture, and affects photosynthesis by targeting PdeBRC1 and controlling the anatomy of the mesophyll, and intercellular space, respectively. Branching, secondary growth, and photosynthesis are vital developmental processes of woody plants that determine plant architecture and timber yield. However, the mechanisms underlying these processes are unknown. Here, we report that the Populus transcription factor High Cambium Activity 2 (PdeHCA2) plays a role in the transition from primary to secondary growth, vascular development, and branching. In Populus, PdeHCA2 is expressed in undifferentiated provascular cells during primary growth, in phloem cells during secondary growth, and in leaf veins, which is different from the expression pattern of its homolog in Arabidopsis. Overexpression of PdeHCA2 has pleiotropic effects on shoot and leaf development; overexpression lines showed delayed growth of shoots and leaves, reduced photosynthesis, and abnormal shoot branching. In addition, auxin-, cytokinin-, and photosynthesis-related genes were differentially regulated in these lines. Electrophoretic mobility shift assays and transcriptome analysis indicated that PdeHCA2 directly up-regulates the expression of BRANCHED1 and the MADS-box gene PdeAGL9, which regulate plant architecture, by binding to cis-elements in the promoters of these genes. Taken together, our findings suggest that HCA2 regulates several processes in woody plants including vascular development, photosynthesis, and branching by affecting the proliferation and differentiation of parenchyma cells.

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

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

A Multilayer Interactome Network Constructed in a Forest Poplar Population Mediates the Pleiotropic Control of Complex Traits

The effects of genes on physiological and biochemical processes are interrelated and interdependent; it is common for genes to express pleiotropic control of complex traits. However, the study of gene expression and participating pathways in vivo at the whole-genome level is challenging. Here, we develop a coupled regulatory interaction differential equation to assess overall and independent genetic effects on trait growth. Based on evolutionary game theory and developmental modularity theory, we constructed multilayer, omnigenic networks of bidirectional, weighted, and positive or negative epistatic interactions using a forest poplar tree mapping population, which were organized into metagalactic, intergalactic, and local interstellar networks that describe layers of structure between modules, submodules, and individual single nucleotide polymorphisms, respectively. These multilayer interactomes enable the exploration of complex interactions between genes, and the analysis of not only differential expression of quantitative trait loci but also previously uncharacterized determinant SNPs, which are negatively regulated by other SNPs, based on the deconstruction of genetic effects to their component parts. Our research framework provides a tool to comprehend the pleiotropic control of complex traits and explores the inherent directional connections between genes in the structure of omnigenic networks.

Transcriptome-wide identification and characterization of microRNAs in diverse phases of wood formation in Populus trichocarpa

We applied miRNA expression profiling method to Populus trichocarpa stems of the three developmental stages, primary stem (PS), transitional stem (TS), and secondary stem (SS), to investigate miRNA species and their regulation on lignocellulosic synthesis and related processes. We obtained 892, 872, and 882 known miRNAs and 1727, 1723, and 1597 novel miRNAs, from PS, TS, and SS, respectively. Comparisons of these miRNA species among different developmental stages led to the identification of 114, 306, and 152 differentially expressed miRNAs (DE-miRNAs), which had 921, 2639, and 2042 candidate target genes (CTGs) in the three respective stages of the same order. Correlation analysis revealed 47, 439, and 71 DE-miRNA-CTG pairs of high negative correlation in PS, TS, and SS, respectively. Through biological process analysis, we finally identified 34, 6, and 76 miRNA-CTG pairs from PS, TS, and SS, respectively, and the miRNA target genes in these pairs regulate or participate lignocellulosic biosynthesis-related biological processes: cell division and differentiation, cell wall modification, secondary cell wall biosynthesis, lignification, and programmed cell death processes. This is the first report on an integrated analysis of genome-wide mRNA and miRNA profilings during multiple phases of poplar stem development. Our analysis results imply that individual miRNAs modulate secondary growth and lignocellulosic biosynthesis through regulating transcription factors and lignocellulosic biosynthetic pathway genes, resulting in more dynamic promotion, suppression, or regulatory circuits. This study advanced our understanding of many individual miRNAs and their essential, diversified roles in the dynamic regulation of secondary growth in woody tree species.

Functional understanding of secondary cell wall cellulose synthases in Populus trichocarpa via the Cas9/gRNA-induced gene knockouts

Plant cellulose is synthesized by a large plasma membrane-localized cellulose synthase (CesA) complex. However, an overall functional determination of secondary cell wall (SCW) CesAs is still lacking in trees, especially one based on gene knockouts. Here, the Cas9/gRNA-induced knockouts of PtrCesA4, 7A, 7B, 8A and 8B genes were produced in Populus trichocarpa. Based on anatomical, immunohistochemical and wood composition evidence, we gained a comprehensive understanding of five SCW PtrCesAs at the genetic level. Complete loss of PtrCesA4, 7A/B or 8A/B led to similar morphological abnormalities, indicating similar and nonredundant genetic functions. The absence of the gelatinous (G) layer, one-layer-walled fibres and a 90% decrease in cellulose in these mutant woods revealed that the three classes of SCW PtrCesAs are essential for multilayered SCW structure and wood G-fibre. In addition, the mutant primary and secondary phloem fibres lost the n(G + L)- and G-layers and retained the thicker S-layers (L, lignified; S, secondary). Together with polysaccharide immunolocalization data, these findings suggest differences in the role of SCW PtrCesAs-synthesized cellulose in wood and phloem fibre wall structures. Overall, this functional understanding of the SCW PtrCesAs provides further insights into the impact of lacking cellulose biosynthesis on growth, SCW, wood G-fibre and phloem fibre wall structures in the tree.

Comparative Analyses of Full-Length Transcriptomes Reveal Gnetum luofuense Stem Developmental Dynamics

Genus Gnetum, of which the majority species are pantropical liana, have broad industrial uses including for string, nets, and paper production. Although numerous studies have investigated anatomical structures during stem development, the underlying molecular mechanisms that regulate this developmental trajectory in Gnetum species remain poorly understood. A total of 12 full-length transcriptomes were generated from four stem developmental stages of an arborescent representative of this genus, Gnetum luofuense, using Oxford Nanopore Technologies. The results of this analysis reveal a total of 24,151 alternative splicing (AS) and 134,391 alternative polyadenylation events. A remarkably dynamic pattern of AS events, especially in the case of intron retentions, was found across the four developmental stages while no dynamic pattern was found among transcript numbers with varied poly(A) sites. A total of 728 long non-coding RNAs were also detected; the number of cis-regulated target genes dramatically increased while no changes were found among trans-regulated target genes. In addition, a K-means clustering analysis of all full-length transcripts revealed that primary growth is associated with carbohydrate metabolism and fungi defense, while secondary growth is closely linked with photosynthesis, nitrogen transportation, and leaf ontogenesis. The use of weighted gene co-expression network analysis as well as differentially expressed transcripts reveals that bHLH, GRF, and MYB-related transcription factors are involved in primary growth, while AP2/ERF, MYB, NAC, PLAZ, and bZIP participate in G. luofuense stem secondary growth. The results of this study provide further evidence that Nanopore sequencing technology provides a cost-effective method for generating full-length transcriptome data as well as for investigating seed plant organ development.

Fasciclin-Like Arabinogalactan-Protein 16 (FLA16) Is Required for Stem Development in Arabidopsis

The predominant Fascilin 1 (FAS1)-containing proteins in plants belong to the Fasciclin-Like Arabinogalactan-protein (FLA) family of extracellular glycoproteins. In addition to FAS1 domains, these multi-domain FLA proteins contain glycomotif regions predicted to direct addition of large arabinogalactan (AG) glycans and many contain signal sequences for addition of a glycosylphosphatidylinositol (GPI)-anchor to tether them to the plasma membrane. FLAs are proposed to play both structural and signaling functions by forming a range of interactions in the plant extracellular matrix, similar to FAS1-containing proteins in animals. FLA group B members contain two FAS1 domains and are not predicted to be GPI-anchored. None of the group B members have been functionally characterized or their sub-cellular location resolved, limiting understanding of their function. We investigated the group B FLA16 in Arabidopsis that is predominantly expressed in inflorescence tissues. FLA16 is the most highly expressed FLA in the stem after Group A members FLA11 and FLA12 that are stem specific. A FLA16-YFP fusion protein driven by the endogenous putative FLA16 promoter in wild type background showed expression in cells with secondary cell walls, and FLA16 displayed characteristics of cell wall glycoproteins with moderate glycosylation. Investigation of a fla16 mutant showed loss of FLA16 leads to reduced stem length and altered biomechanical properties, likely as a result of reduced levels of cellulose. Immuno-labeling indicated support for FLA16 location to the plasma-membrane and (apoplastic) cell wall of interfascicular stem fiber cells. Together these results indicate FLA16, a two-FAS1 domain FLAs, plays a role in plant secondary cell wall synthesis and function.

Comparative transcriptome analyses define genes and gene modules differing between two Populus genotypes with contrasting stem growth rates

BACKGROUND

Wood provides an important biomass resource for biofuel production around the world. The radial growth of tree stems is central to biomass production for forestry and biofuels, but it is challenging to dissect genetically because it is a complex trait influenced by many genes. In this study, we adopted methods of physiology, transcriptomics and genetics to investigate the regulatory mechanisms of tree radial growth and wood development.

RESULTS

Physiological comparison showed that two Populus genotypes presented different rates of radial growth of stems and accumulation of woody biomass. A comparative transcriptional network approach was used to define and characterize functional differences between two Populus genotypes. Analyses of transcript profiles from wood-forming tissue of the two genotypes showed that 1542, 2295 and 2110 genes were differentially expressed in the pre-growth, fast-growth and post-growth stages, respectively. The co-expression analyses identified modules of co-expressed genes that displayed distinct expression profiles. Modules were further characterized by correlating transcript levels with genotypes and physiological traits. The results showed enrichment of genes that participated in cell cycle and division, whose expression change was consistent with the variation of radial growth rates. Genes related to secondary vascular development were up-regulated in the faster-growing genotype in the pre-growth stage. We characterized a BEL1-like (BELL) transcription factor, PeuBELL15, which was up-regulated in the faster-growing genotype. Analyses of transgenic Populus overexpressing as well as CRISPR/Cas9-induced mutants for BELL15 showed that PeuBELL15 improved accumulation of glucan and lignin, and it promoted secondary vascular growth by regulating the expression of genes relevant for cellulose synthases and lignin biosynthesis.

CONCLUSIONS

This study illustrated that active division and expansion of vascular cambium cells and secondary cell wall deposition of xylem cells contribute to stem radial increment and biomass accumulation, and it identified relevant genes for these complex growth traits, including a BELL transcription factor gene PeuBELL15. This provides genetic resources for improving and breeding elite genotypes with fast growth and high wood biomass.

TDIF overexpression in poplars retards internodal elongation and enhances leaf venation through interaction with other phytohormones

As a member of the CLAVATA3 (CLV3)/EMBRYO SURROUNDING REGION-related (CLE) peptide family, tracheary element differentiation inhibitory factor (TDIF) plays crucial roles in vascular meristem maintenance by promoting cell proliferation and inhibiting xylem cell differentiation. In Populus trichocarpa, six TDIF-encoding genes are all expressed in vascular tissues, and in Arabidopsis PtTDIFpro:GUS lines, the expression driven by PtTDIF promoters were predominantly detected in stem vascular bundles, initiating leaves and leaf veins. Although exogenous application of two poplar TDIF peptides did not evidently affect the shoot growth in vitro, overexpression of PtTDIF genes in hybrid poplar severely retarded the internodal elongation by upregulating the expression of GA2ox and GA20ox genes and thus decreasing the level of endogenous gibberellins (GAs), which phenotypic defect could be rescued by exogenously applied GA3. In addition, TDIF overexpression unexpectedly induced a more complex venation pattern in poplar leaves, which was underpinned by the elevated expression of WOX4 and WOX13 genes. Upon TDIF treatment, the DR5:GUS poplar leaves revealed a higher GUS activity and in TDIF-overexpressing leaves, the transcript abundances of several PIN-FORMED (PIN) genes, especially that of PIN1, were increased, which implied an integration of TDIF and auxin in mediating this process. Collectively, data of this work presented novel activities of TDIF involved in internode elongation and leaf vein formation, thus revealing the divergent functions of TDIF in perennial tree species from those in annual herbaceous Arabidopsis.

Allelic variation in candidate genes associated with wood properties of cultivated poplars (Populus)

INTRODUCTION

Since Populus has veritable value as timber, plywood, pulp, and paper, genomic research should create the sound basis for further breeding toward desirable wood quality attributes.

MATERIALS AND METHODS

In this study, we addressed the need for a research methodology that initially identifies and then characterize candidate genes encoding enzymes with wood property phenotypic traits, toward the aim of developing a genomics-based breeding technology.

RESULTS

On 23 different poplar species/hybrid samples, we successfully amplified 55 primers designed on Populus trichocarpa L. Considering the number of polymorphic sites, out of 73,206 bp, 51 SNPs and 31 indel events were found. Non-synonymous single base mutations could be detected in number of 30, 21 out of 164 sequences were the number of minimum recombination events and 41 significant pairwise comparisons between loci could be detected.

DISCUSSION AND CONCLUSION

Our results provide a roadmap for a future association genetic study between nucleotide diversity and precise evaluation of phenotype.

Expression profiles of cell-wall related genes vary broadly between two common maize inbreds during stem development

Background

The cellular machinery for cell wall synthesis and metabolism is encoded by members of large multi-gene families. Maize is both a genetic model for grass species and a potential source of lignocellulosic biomass from crop residues. Genetic improvement of maize for its utility as a bioenergy feedstock depends on identification of the specific gene family members expressed during secondary wall development in stems.

Results

High-throughput sequencing of transcripts expressed in developing rind tissues of stem internodes provided a comprehensive inventory of cell wall-related genes in maize (Zea mays, cultivar B73). Of 1239 of these genes, 854 were expressed among the internodes at ≥95 reads per 20 M, and 693 of them at ≥500 reads per 20 M. Grasses have cell wall compositions distinct from non-commelinid species; only one-quarter of maize cell wall-related genes expressed in stems were putatively orthologous with those of the eudicot Arabidopsis. Using a slope-metric algorithm, five distinct patterns for sub-sets of co-expressed genes were defined across a time course of stem development. For the subset of genes associated with secondary wall formation, fifteen sequence motifs were found in promoter regions. The same members of gene families were often expressed in two maize inbreds, B73 and Mo17, but levels of gene expression between them varied, with 30% of all genes exhibiting at least a 5-fold difference at any stage. Although presence-absence and copy-number variation might account for much of these differences, fold-changes of expression of a CADa and a FLA11 gene were attributed to polymorphisms in promoter response elements.

Conclusions

Large genetic variation in maize as a species precludes the extrapolation of cell wall-related gene expression networks even from one common inbred line to another. Elucidation of genotype-specific expression patterns and their regulatory controls will be needed for association panels of inbreds and landraces to fully exploit genetic variation in maize and other bioenergy grass species.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

The developmental dynamics of the Populus stem transcriptome

The Populus shoot undergoes primary growth (longitudinal growth) followed by secondary growth (radial growth), which produces biomass that is an important source of energy worldwide. We adopted joint PacBio Iso-Seq and RNA-seq analysis to identify differentially expressed transcripts along a developmental gradient from the shoot apex to the fifth internode of Populus Nanlin895. We obtained 87 150 full-length transcripts, including 2081 new isoforms and 62 058 new alternatively spliced isoforms, most of which were produced by intron retention, that were used to update the Populus annotation. Among these novel isoforms, there are 1187 long non-coding RNAs and 356 fusion genes. Using this annotation, we found 15 838 differentially expressed transcripts along the shoot developmental gradient, of which 1216 were transcription factors (TFs). Only a few of these genes were reported previously. The differential expression of these TFs suggests that they may play important roles in primary and secondary growth. AP2, ARF, YABBY and GRF TFs are highly expressed in the apex, whereas NAC, bZIP, PLATZ and HSF TFs are likely to be important for secondary growth. Overall, our findings provide evidence that long-read sequencing can complement short-read sequencing for cataloguing and quantifying eukaryotic transcripts and increase our understanding of the vital and dynamic process of shoot development.

Differential effects of nickel dosages on in vitro and in vivo seed germination and expression of a high affinity nickel-transport family protein (AT2G16800) in trembling aspen (Populus tremuloides)

It has been demonstrated that a number of metals including mercury (Hg), zinc (Zn), cadmium (Cd), cobalt (Co), lead (Pb), copper (Cu), and nickel (Ni) decrease seed germination rates and plant growth. The threshold levels of metal toxicity on seed germination, plant development, and gene regulation have not been studied in detail. The main objective of this study was to assess in vitro and in vivo the effects of different doses of nickel on Trembling aspen (Populus tremuloides) seed germination and regulation of the high affinity nickel transporter family protein (AT2G16800) gene. The in vitro assays showed that Nickel completely inhibited seed germination even at the lowest concentration of 0.401 mg Ni per mL (in media) tested. However, when the same concentration of nickel (150 mg Ni per 1 kg of dry soil) was added to soil samples, during the vivo assays, almost all of the seeds germinated. Significant inhibition of seed germination was observed when soil samples were treated with at least 400 mg/kg of Ni. No damages were observed on growing seedlings treated with 150, 400, and 800 mg/kg of Ni. Only the highest dose of 1, 600 mg/kg resulted in visible leaf and stem damages and reduced growth on 75% of seedlings. A significant repression of the AT2G16800 gene was observed for the 400, 800, and 1600 mg/kg of nickel treatments compared to the water control with the lowest level of expression observed in samples treated with 800 mg/kg of Ni. Results of this study suggest that P. tremuloides populations will likely be sustainable for long term in sites that are highly contaminated with Ni including mining regions since the bioavailable amount of this metal is usually below 400 mg/kg in Canada.

Identification and expression profiling of genes governing lignin biosynthesis in Casuarina equisetifolia L

Casuarina equisetifolia L. is an important multi-purpose, fast growing and widely planted tree species native to tropical and subtropical coastlines of Australia, Southeast Asia, Malaysia, Melanesia, Polynesia and New Caledonia. It is a nitrogen-fixing tree mainly used for charcoal making, construction poles, landscaping, timber, pulp, firewood, windbreaks, shelterbelts, soil erosion and sand dune stabilization. Casuarina wood is presently used for paper and pulp production. Raw material with reduced lignin is highly preferred to increase the pulp yield. Hence, understanding the molecular regulation of wood formation in this tree species is vital for selecting industrially suitable phenotypes for breeding programs. The lignin biosynthetic pathway has been extensively studied in tree species like Eucalypts, poplars, pines, Picea, Betula and Acacia sp. However, studies on wood formation at molecular level is presently lacking in casuarinas. Hence, in the present study, the transcriptome of the developing secondary tissues of 15 years old Casuarina equiseitfolia subsp. equisetifolia was sequenced, de novo assembled, annotated and mapped to functional pathways. Transcriptome sequencing generated a total of 26,985 transcripts mapped to 31 pathways. Mining of the annotated data identified nine genes involved in lignin biosynthesis pathway and relative expression of the transcripts in four tissues including scale-like leaves, needle-like brachlets, wood and root were documented. The expression of CeCCR1 and CeF5H were found to be significantly high in wood tissues, while maximum expression of CeHCT was documented in stem. Additionally, CeTUBA and CeH2A were identified as the most stable reference transcript for normalization of qRT-PCR data in C. equisetifolia. The present study is the first wood genomic resource in C. equisetifolia, which will be valuable for functional genomics research in this genus.

Plant lipid transfer proteins: are we finally closing in on the roles of these enigmatic proteins?

The nonspecific lipid transfer proteins (LTPs) are small compact proteins folded around a tunnel-like hydrophobic cavity, making them suitable for lipid binding and transport. LTPs are encoded by large gene families in all land plants, but they have not been identified in algae or any other organisms. Thus, LTPs are considered key proteins for plant survival on and colonization of land. LTPs are abundantly expressed in most plant tissues, both above and below ground. They are usually localized to extracellular spaces outside the plasma membrane. Although the in vivo functions of LTPs remain unclear, accumulating evidence suggests a role for LTPs in the transfer and deposition of monomers required for assembly of the waterproof lipid barriers, such as cutin and cuticular wax, suberin, and sporopollenin, formed on many plant surfaces. Some LTPs may be involved in other processes, such as signaling during pathogen attacks. Here, we present the current status of LTP research with a focus on the role of these proteins in lipid barrier deposition and cell expansion. We suggest that LTPs facilitate extracellular transfer of barrier materials and adhesion between barriers and extracellular materials. A growing body of research may uncover the true role of LTPs in plants.

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

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

Transcriptome analysis provides insights into xylogenesis formation in Moso bamboo (Phyllostachys edulis) shoot

Maturation-related changes in cell wall composition and the molecular mechanisms underlying cell wall changes were investigated from the apical, middle and basal segments in moso bamboo shoot (MBS). With maturation extent from apical to basal regions in MBS, lignin and cellulose content increased, whereas heteroxylan exhibited a decreasing trend. Activities of phenylalanine amonnialyase (PAL), cinnamyl alcohol dehydrogenase (CAD) and cinnamate-4-hydroxylase (C4H), which are involved in lignin biosynthesis, increased rapidly from the apex to the base sections. The comparative transcriptomic analysis was carried out to identify some key genes involved in secondary cell walls (SCW) formation underlying the cell wall compositions changes including 63, 8, 18, and 31 functional unigenes encoding biosynthesis of lignin, cellulose, xylan and NAC-MYB-based transcription factors, respectively. Genes related to secondary cell wall formation and lignin biosynthesis had higher expression levels in the middle and basal segments compared to those in the apical segments. Furthermore, the expression profile of PePAL gene showed positive relationships with cellulose-related gene PeCESA4, xylan-related genes PeIRX9 and PeIRX10. Our results indicated that lignification occurred in the more mature middle and basal segments in MBS at harvest while lignification of MBS were correlated with higher expression levels of PeCESA4, PeIRX9 and PeIRX10 genes.

Transcriptome analysis of genes involved in secondary cell wall biosynthesis in developing internodes of Miscanthus lutarioriparius

Miscanthus is a promising lignocellulosic bioenergy crop for bioethanol production. To identify candidate genes and regulation networks involved in secondary cell wall (SCW) development in Miscanthus, we performed de novo transcriptome analysis of a developing internode. According to the histological and in-situ histochemical analysis, an elongating internode of M. lutarioriparius can be divided into three distinct segments, the upper internode (UI), middle internode (MI) and basal internode (BI), each representing a different stage of SCW development. The transcriptome analysis generated approximately 300 million clean reads, which were de novo assembled into 79,705 unigenes. Nearly 65% of unigenes was annotated in seven public databases. Comparative profiling among the UI, MI and BI revealed four distinct clusters. Moreover, detailed expression profiling was analyzed for gene families and transcription factors (TFs) involved in SCW biosynthesis, assembly and modification. Based on the co-expression patterns, putative regulatory networks between TFs and SCW-associated genes were constructed. The work provided the first transcriptome analysis of SCW development in M. lutarioriparius. The results obtained provide novel insights into the biosynthesis and regulation of SCW in Miscanthus. In addition, the genes identified represent good candidates for further functional studies to unravel their roles in SCW biosynthesis and modification.

Characteristics of hepatic insulin-sensitive nonalcoholic fatty liver disease

Nonalcoholic fatty liver disease (NAFLD) plays a crucial role in type 2 diabetes and hepatocellular carcinoma. The major underlying pathogenesis is hepatic insulin resistance. The aim of the present study was to characterize patients with NAFLD with paradoxically normal hepatic insulin sensitivity relative to patients with NAFLD with hepatic insulin resistance. We recruited 26 patients with NAFLD and divided them into three groups ranked by the level of hepatic insulin sensitivity (HIS; high‐HIS, mid‐HIS, low‐HIS), as assessed by the hyperinsulinemic‐euglycemic clamp studies using stable isotope. Hepatic insulin sensitivity of the high‐HIS group was identical to that of the non‐NAFLD lean control (clamped percent suppression of endogenous glucose production, 91.1% ± 5.2% versus 91.0% ± 8.5%, respectively) and was significantly higher than that of the low‐HIS group (66.6% ± 7.5%; P < 0.01). Adiposity (subcutaneous, visceral, intrahepatic, and muscular lipid content), hepatic histopathology, and expression levels of various genes by using liver biopsies, muscle, and adipose tissue insulin sensitivity, plasma metabolites by metabolomics analysis, putative biomarkers, and lifestyles were assessed and compared between the high‐HIS and low‐HIS groups. Among these, adipose tissue insulin sensitivity assessed by clamped percent suppression of free fatty acid, serum high molecular weight adiponectin, and plasma tricarboxylic acid cycle metabolites, such as citric acid and cis‐aconitic acid, were significantly higher in the high‐HIS group compared to the low‐HIS group. In contrast, there were no differences in adiposity, including intrahepatic lipid content assessed by proton magnetic resonance spectroscopy (28.3% ± 16.1% versus 20.4% ± 9.9%, respectively), hepatic histopathology, other putative biomarkers, and lifestyles. Conclusion: High levels of adipose tissue insulin sensitivity, serum high molecular weight adiponectin, and plasma tricarboxylic acid cycle metabolites are unique characteristics that define patients with hepatic insulin‐sensitive NAFLD regardless of intrahepatic lipid content. (Hepatology Communications 2017;1:634–647)

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

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

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

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

Identifying gene coexpression networks underlying the dynamic regulation of wood-forming tissues in Populus under diverse environmental conditions

Trees modify wood formation through integration of environmental and developmental signals in complex but poorly defined transcriptional networks, allowing trees to produce woody tissues appropriate to diverse environmental conditions. In order to identify relationships among genes expressed during wood formation, we integrated data from new and publically available datasets in Populus. These datasets were generated from woody tissue and include transcriptome profiling, transcription factor binding, DNA accessibility and genome-wide association mapping experiments. Coexpression modules were calculated, each of which contains genes showing similar expression patterns across experimental conditions, genotypes and treatments. Conserved gene coexpression modules (four modules totaling 8398 genes) were identified that were highly preserved across diverse environmental conditions and genetic backgrounds. Functional annotations as well as correlations with specific experimental treatments associated individual conserved modules with distinct biological processes underlying wood formation, such as cell-wall biosynthesis, meristem development and epigenetic pathways. Module genes were also enriched for DNase I hypersensitivity footprints and binding from four transcription factors associated with wood formation. The conserved modules are excellent candidates for modeling core developmental pathways common to wood formation in diverse environments and genotypes, and serve as testbeds for hypothesis generation and testing for future studies.

Dynamic changes in the transcriptome of Populus hopeiensis in response to abscisic acid

Abscisic acid (ABA) plays a fundamental role in plant response and adaptation to abiotic stresses, such as drought, high salinity and low temperature. Populus hopeiensis exhibits exceptional tolerance to water-deficit environments and is therefore an excellent choice for studying drought tolerance in trees. This study provides a global view of transcriptome dynamics in P. hopeiensis in response to exogenous ABA using Illumina RNA-sequencing. Endogenous ABA content increased and reached a peak at 8 h after ABA treatment and then significantly decreased at latter time points. Differential expression analysis and Gene ontology enrichment revealed that the number of transcripts exhibited significant increase during the first 8 hours after ABA treatment, which then significantly decreased at 12 and 24 h. Transcription factors (TFs) analysis showed that six different patterns were observed based on the expression of the six TFs families (AP2/ERF, NAC, MYB, MYB-related, bZIP and WRKY) and the majority of differentially expressed TFs increased rapidly after ABA treatment. This study provides a robust resource for investigating the functions of genes induced by ABA and will help to develop a better understanding of the molecular regulatory mechanism in response to drought in poplar.

Field Guide to Plant Model Systems

For the past several decades, advances in plant development, physiology, cell biology, and genetics have relied heavily on the model (or reference) plant Arabidopsis thaliana. Arabidopsis resembles other plants, including crop plants, in many but by no means all respects. Study of Arabidopsis alone provides little information on the evolutionary history of plants, evolutionary differences between species, plants that survive in different environments, or plants that access nutrients and photosynthesize differently. Empowered by the availability of large-scale sequencing and new technologies for investigating gene function, many new plant models are being proposed and studied.

CsWRKY46, a WRKY transcription factor from cucumber, confers cold resistance in transgenic-plant by regulating a set of cold-stress responsive genes in an ABA-dependent manner

Plant WRKY transcription factors are trans-regulatory proteins that are involved in plant immune responses, development and senescence; however, their roles in abiotic stress are still not well understood, especially in the horticultural crop cucumber. In this study, a novel cucumber WRKY gene, CsWRKY46 was cloned and identified, which was up-regulated in response to cold stress and exogenous abscisic acid (ABA) treatment. CsWRKY46 is belonging to group II of the WRKY family, CsWRKY46 was found exclusively in the nucleus, as indicated by a transient expression assay. Yeast one-hybrid assay shown that CsWRKY46 interact with the W-box in the promoter of ABI5. Transgenic Arabidopsis lines over-expressing CsWRKY46, WRK46-OE1 and WRK46-OE5 had higher seedling survival rates upon freezing treatment compared with that of the wild-type. The above over-expression lines also showed much a higher proline accumulation, less electrolyte leakage and lower malondialdehyde (MDA) levels. Furthermore, the CsWRKY46 overexpression lines were hypersensitive to ABA during seed germination, but the seedlings were not. Quantitative RT-PCR analyses revealed that the expression levels of the ABA-responsive transcription factor ABI5 were higher in the WRKY46-OE lines than in wild-type and that the overexpression of CsWRKY46 increased the expression of stress-inducible genes, including RD29A and COR47. Taken together, our results demonstrated that CsWRKY46 from cucumber conferred cold tolerance to transgenic plants and positively regulated the cold signaling pathway in an ABA-dependent manner.

Transcriptional regulation of vascular cambium activity during the transition from juvenile to mature stages in Cunninghamia lanceolata

Cunninghamia lanceolata (Lamb.) Hook., an evergreen conifer distributed in southern China, has been recognized as the most commercially important timber species due to its rapid growth. However, the molecular mechanisms underlying growth alternation due to vascular cambium activity are poorly understood. Here, we used cryosectioning to isolate the vascular cambium tissue of C. lanceolata at three stages, namely, juvenile, transition and mature (3-, 13-, and 35-year-old trees respectively) for transcriptome-wide analysis. Through assembling and annotation of transcripts, 108,767 unigenes and some potential growth-regulated genes were identified. A total of 5213, 4873 and 2541 differentially expressed genes (DEGs) were identified in the three stages. DEGs related to cambial activity, cell division and cell wall modification were detected at various developmental stages of the vascular cambium. In addition, some putative genes involved in plant hormone biosynthesis were also differentially regulated. These results indicate that various cambium-related molecular activities result in alterations in the growth of C. lanceolata, particularly during the transition from juvenile to mature stages. The findings of the present study improve our understanding of cambium development and may aid in studies of the molecular mechanisms of wood production and provide fundamental insights into the establishment of the optimal rotation period for silvicultural trees.

Validation of Reference Genes for Gene Expression by Quantitative Real-Time RT-PCR in Stem Segments Spanning Primary to Secondary Growth in Populus tomentosa

The vertical segments of Populus stems are an ideal experimental system for analyzing the gene expression patterns involved in primary and secondary growth during wood formation. Suitable internal control genes are indispensable to quantitative real time PCR (qRT-PCR) assays of gene expression. In this study, the expression stability of eight candidate reference genes was evaluated in a series of vertical stem segments of Populus tomentosa. Analysis through software packages geNorm, NormFinder and BestKeeper showed that genes ribosomal protein (RP) and tubulin beta (TUBB) were the most unstable across the developmental stages of P. tomentosa stems, and the combination of the three reference genes, eukaryotic translation initiation factor 5A (eIF5A), Actin (ACT6) and elongation factor 1-beta (EF1-beta) can provide accurate and reliable normalization of qRT-PCR analysis for target gene expression in stem segments undergoing primary and secondary growth in P. tomentosa. These results provide crucial information for transcriptional analysis in the P. tomentosa stem, which may help to improve the quality of gene expression data in these vertical stem segments, which constitute an excellent plant system for the study of wood formation.

Large-scale transcriptional profiling of lignified tissues in Tectona grandis

BACKGROUND

Currently, Tectona grandis is one of the most valuable trees in the world and no transcript dataset related to secondary xylem is available. Considering how important the secondary xylem and sapwood transition from young to mature trees is, little is known about the expression differences between those successional processes and which transcription factors could regulate lignin biosynthesis in this tropical tree. Although MYB transcription factors are one of the largest superfamilies in plants related to secondary metabolism, it has not yet been characterized in teak. These results will open new perspectives for studies of diversity, ecology, breeding and genomic programs aiming to understand deeply the biology of this species.

RESULTS

We present a widely expressed gene catalog for T. grandis using Illumina technology and the de novo assembly. A total of 462,260 transcripts were obtained, with 1,502 and 931 genes differentially expressed for stem and branch secondary xylem, respectively, during age transition. Analysis of stem and branch secondary xylem indicates substantial similarity in gene ontologies including carbohydrate enzymes, response to stress, protein binding, and allowed us to find transcription factors and heat-shock proteins differentially expressed. TgMYB1 displays a MYB domain and a predicted coiled-coil (CC) domain, while TgMYB2, TgMYB3 and TgMYB4 showed R2R3-MYB domain and grouped with MYBs from several gymnosperms and flowering plants. TgMYB1, TgMYB4 and TgCES presented higher expression in mature secondary xylem, in contrast with TgMYB2, TgHsp1, TgHsp2, TgHsp3, and TgBi whose expression is higher in young lignified tissues. TgMYB3 is expressed at lower level in secondary xylem.

CONCLUSIONS

Expression patterns of MYB transcription factors and heat-shock proteins in lignified tissues are dissimilar when tree development was evaluated, obtaining more expression of TgMYB1 and TgMYB4 in lignified tissues of 60-year-old trees, and more expression in TgHsp1, TgHsp2, TgHsp3 and TgBi in stem secondary xylem of 12-year-old trees. We are opening a door for further functional characterization by reverse genetics and marker-assisted selection with those genes. Investigation of some of the key regulators of lignin biosynthesis in teak, however, could be a valuable step towards understanding how rigidity of teak wood and extractives content are different from most other woods. The obtained transcriptome data represents new sequences of T. grandis deposited in public databases, representing an unprecedented opportunity to discover several related-genes associated with secondary xylem such as transcription factors and stress-related genes in a tropical tree.

Quantitative Proteomic and Phosphoproteomic Approaches for Deciphering the Signaling Pathway for Tension Wood Formation in Poplar

Trees adjust their growth following forced changes in orientation to re-establish a vertical position. In angiosperms, this adjustment involves the differential regulation of vascular cambial activity between the lower (opposite wood) and upper (tension wood) sides of the leaning stem. We investigated the molecular mechanisms leading to the formation of differential wood types through a quantitative proteomic and phosphoproteomic analysis on poplar subjected to a gravitropic stimulus. We identified and quantified 675 phosphopeptides, corresponding to 468 phosphoproteins, and 3 763 nonphosphorylated peptides, corresponding to 1 155 proteins, in the differentiating xylem of straight-growing trees (control) and trees subjected to a gravitational stimulus during 8 weeks. About 1% of the peptides were specific to a wood type (straight, opposite, or tension wood). Proteins quantified in more than one type of wood were more numerous: a mixed linear model showed 389 phosphopeptides and 556 proteins to differ in abundance between tension wood and opposite wood. Twenty-one percent of the phosphoproteins identified here were described in their phosphorylated form for the first time. Our analyses revealed remarkable developmental molecular plasticity, with wood type-specific phosphorylation events, and highlighted the involvement of different proteins in the biosynthesis of cell wall components during the formation of the three types of wood.

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.

Aldehyde Dehydrogenase Gene Superfamily in Populus: Organization and Expression Divergence between Paralogous Gene Pairs

Aldehyde dehydrogenases (ALDHs) constitute a superfamily of NAD(P)+-dependent enzymes that catalyze the irreversible oxidation of a wide range of reactive aldehydes to their corresponding nontoxic carboxylic acids. ALDHs have been studied in many organisms from bacteria to mammals; however, no systematic analyses incorporating genome organization, gene structure, expression profiles, and cis-acting elements have been conducted in the model tree species Populus trichocarpa thus far. In this study, a comprehensive analysis of the Populus ALDH gene superfamily was performed. A total of 26 Populus ALDH genes were found to be distributed across 12 chromosomes. Genomic organization analysis indicated that purifying selection may have played a pivotal role in the retention and maintenance of PtALDH gene families. The exon-intron organizations of PtALDHs were highly conserved within the same family, suggesting that the members of the same family also may have conserved functionalities. Microarray data and qRT-PCR analysis indicated that most PtALDHs had distinct tissue-specific expression patterns. The specificity of cis-acting elements in the promoter regions of the PtALDHs and the divergence of expression patterns between nine paralogous PtALDH gene pairs suggested that gene duplications may have freed the duplicate genes from the functional constraints. The expression levels of some ALDHs were up- or down-regulated by various abiotic stresses, implying that the products of these genes may be involved in the adaptation of Populus to abiotic stresses. Overall, the data obtained from our investigation contribute to a better understanding of the complexity of the Populus ALDH gene superfamily and provide insights into the function and evolution of ALDH gene families in vascular plants.

Genome-wide identification, classification, and expression analysis of 14-3-3 gene family in Populus

Background

In plants, 14-3-3 proteins are encoded by a large multigene family and are involved in signaling pathways to regulate plant development and protection from stress. Although twelve Populus 14-3-3s were identified based on the Populus trichocarpa genome V1.1 in a previous study, no systematic analysis including genome organization, gene structure, duplication relationship, evolutionary analysis and expression compendium has been conducted in Populus based on the latest P. trichocarpa genome V3.0.

Principal Findings

Here, a comprehensive analysis of Populus 14-3-3 family is presented. Two new 14-3-3 genes were identified based on the latest P. trichocarpa genome. In P. trichocarpa, fourteen 14-3-3 genes were grouped into ε and non-ε group. Exon-intron organizations of Populus 14-3-3s are highly conserved within the same group. Genomic organization analysis indicated that purifying selection plays a pivotal role in the retention and maintenance of Populus 14-3-3 family. Protein conformational analysis indicated that Populus 14-3-3 consists of a bundle of nine α-helices (α1-α9); the first four are essential for formation of the dimer, while α3, α5, α7, and α9 form a conserved peptide-binding groove. In addition, α1, α3, α5, α7, and α9 were evolving at a lower rate, while α2, α4, and α6 were evolving at a relatively faster rate. Microarray analyses showed that most Populus 14-3-3s are differentially expressed across tissues and upon exposure to various stresses.

Conclusions

The gene structures and their coding protein structures of Populus 14-3-3s are highly conserved among group members, suggesting that members of the same group might also have conserved functions. Microarray and qRT-PCR analyses showed that most Populus 14-3-3s were differentially expressed in various tissues and were induced by various stresses. Our investigation provided a better understanding of the complexity of the 14-3-3 gene family in poplars.

Endo-β-1,4-glucanases impact plant cell wall development by influencing cellulose crystallization

Cell walls are vital to the normal growth and development of plants as they protect the protoplast and provide rigidity to the stem. Here, two poplar and Arabidopsis orthologous endoglucanases, which have been proposed to play a role in secondary cell wall development, were examined. The class B endoglucanases, PtGH9B5 and AtGH9B5, are secreted enzymes that have a predicted glycosylphosphatidylinositol anchor, while the class C endoglucanases, PtGH9C2 and AtGH9C2, are also predicted to be secreted but instead contain a carbohydrate-binding module. The poplar endoglucanases were expressed in Arabidopsis using both a 35S promoter and the Arabidopsis secondary cell wall-specific CesA8 promoter. Additionally, Arabidopsis t-DNA insertion lines and an RNAi construct was created to downregulate AtGH9C2 in Arabidopsis. All of the plant lines were examined for changes in cell morphology and patterning, growth and development, cell wall crystallinity, microfibril angle, and proportion of cell wall carbohydrates. Misregulation of PtGH9B5/AtGH9B5 resulted in changes in xylose content, while misregulation of PtGH9C2/AtGH9C2 resulted in changes in crystallinity, which was inversely correlated with changes in plant height and rosette diameter. Together, these results suggest that these endoglucanases affect secondary cell wall development by contributing to the cell wall crystallization process.

Hsf and Hsp gene families in Populus: genome-wide identification, organization and correlated expression during development and in stress responses

Background

Heat shock proteins (Hsps) are molecular chaperones that are involved in many normal cellular processes and stress responses, and heat shock factors (Hsfs) are the transcriptional activators of Hsps. Hsfs and Hsps are widely coordinated in various biological processes. Although the roles of Hsfs and Hsps in stress responses have been well characterized in Arabidopsis, their roles in perennial woody species undergoing various environmental stresses remain unclear.

Results

Here, a comprehensive identification and analysis of Hsf and Hsp families in poplars is presented. In Populus trichocarpa, we identified 42 paralogous pairs, 66.7% resulting from a whole genome duplication. The gene structure and motif composition are relatively conserved in each subfamily. Microarray and quantitative real-time RT-PCR analyses showed that most of the Populus Hsf and Hsp genes are differentially expressed upon exposure to various stresses. A coexpression network between Populus Hsf and Hsp genes was generated based on their expression. Coordinated relationships were validated by transient overexpression and subsequent qPCR analyses.

Conclusions

The comprehensive analysis indicates that different sets of PtHsps are downstream of particular PtHsfs and provides a basis for functional studies aimed at revealing the roles of these families in poplar development and stress responses.

Electronic supplementary material

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

Antisense expression of the fasciclin-like arabinogalactan protein FLA6 gene in Populus inhibits expression of its homologous genes and alters stem biomechanics and cell wall composition in transgenic trees

Fasciclin-like arabinogalactan proteins (FLAs) play important roles in the growth and development of roots, stems, and seeds in Arabidopsis. However, their biological functions in woody plants are largely unknown. In this work, we investigated the possible function of PtFLA6 in poplar. Quantitative real-time PCR, PtFLA6-yellow fluorescent protein (YFP) fusion protein subcellular localization, Western blotting, and immunohistochemical analyses demonstrated that the PtFLA6 gene was expressed specifically in the xylem of mature stem, and PtFLA6 protein was distributed ubiquitous in plant cells and accumulated predominantly in stem xylem fibres. Antisense expression of PtFLA6 in the aspen hybrid clone Poplar davidiana×Poplar bolleana reduced the transcripts of PtFLA6 and its homologous genes. Transgenic plants that showed a significant reduction in the transcripts of PtFLAs accumulated fewer PtFLA6 and arabinogalactan proteins than did the non-transgenic plants, leading to reduced stem flexural strength and stiffness. Further studies revealed that the altered stem biomechanics of transgenic plants could be attributed to the decreased cellulose and lignin composition in the xylem. In addition expression of some xylem-specific genes involved in cell wall biosynthesis was downregulated in these transgenic plants. All these results suggest that engineering the expression of PtFLA6 and its homologues could modulate stem mechanical properties by affecting cell wall composition in trees.

RNA sequencing of Populus x canadensis roots identifies key molecular mechanisms underlying physiological adaption to excess zinc

Populus x canadensis clone I-214 exhibits a general indicator phenotype in response to excess Zn, and a higher metal uptake in roots than in shoots with a reduced translocation to aerial parts under hydroponic conditions. This physiological adaptation seems mainly regulated by roots, although the molecular mechanisms that underlie these processes are still poorly understood. Here, differential expression analysis using RNA-sequencing technology was used to identify the molecular mechanisms involved in the response to excess Zn in root. In order to maximize specificity of detection of differentially expressed (DE) genes, we consider the intersection of genes identified by three distinct statistical approaches (61 up- and 19 down-regulated) and validate them by RT-qPCR, yielding an agreement of 93% between the two experimental techniques. Gene Ontology (GO) terms related to oxidation-reduction processes, transport and cellular iron ion homeostasis were enriched among DE genes, highlighting the importance of metal homeostasis in adaptation to excess Zn by P. x canadensis clone I-214. We identified the up-regulation of two Populus metal transporters (ZIP2 and NRAMP1) probably involved in metal uptake, and the down-regulation of a NAS4 gene involved in metal translocation. We identified also four Fe-homeostasis transcription factors (two bHLH38 genes, FIT and BTS) that were differentially expressed, probably for reducing Zn-induced Fe-deficiency. In particular, we suggest that the down-regulation of FIT transcription factor could be a mechanism to cope with Zn-induced Fe-deficiency in Populus. These results provide insight into the molecular mechanisms involved in adaption to excess Zn in Populus spp., but could also constitute a starting point for the identification and characterization of molecular markers or biotechnological targets for possible improvement of phytoremediation performances of poplar trees.

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.

Expression divergence of cellulose synthase (CesA) genes after a recent whole genome duplication event in Populus

Secondary cell wall-associated CesA genes in Populus have undergone a functional differentiation in expression pattern that may be attributable to evolutionary alteration of regulatory modules. Gene duplication is an important mechanism for functional divergence of genes. Secondary cell wall-associated cellulose synthase genes (CesA4, CesA7 and CesA8) are duplicated in Populus plants due to a recent whole genome duplication event. Here, we demonstrate that duplicate CesA genes show tissue-dependent expression divergence in Populus plants. Real-time PCR analysis of Populus CesA genes suggested that Pt × tCesA8-B was more highly expressed than Pt × tCesA8-A in phloem and secondary xylem tissue of mature stem. Histochemical and histological analyses of transformants expressing a GFP-GUS fusion gene driven by Populus CesA promoters revealed that the duplicate CesA genes showed different expression patterns in phloem fibers, secondary xylem, root cap and leaf trichomes. We predicted putative cis-regulatory motifs that regulate expression of secondary cell wall-associated CesA genes, and identified 19 motifs that are highly conserved in the CesA gene family of eudicotyledonous plants. Furthermore, a transient transactivation assay identified candidate transcription factors that affect levels and patterns of expression of Populus CesA genes. The present study reveals that secondary cell wall-associated CesA genes in Populus have undergone a functional differentiation in expression pattern that may be attributable to evolutionary alteration of regulatory modules.