A systems genetics approach to deciphering the effect of dosage variation on leaf morphology in Populus.

Gene copy number variation is frequent in plant genomes of various species, but the impact of such gene dosage variation on morphological traits is poorly understood. We used a large population of Populus carrying genomically characterized insertions and deletions across the genome to systematically assay the effect of gene dosage variation on a suite of leaf morphology traits. A systems genetics approach was used to integrate insertion and deletion locations, leaf morphology phenotypes, gene expression, and transcriptional network data, to provide an overview of how gene dosage influences morphology. Dosage-sensitive genomic regions were identified that influenced individual or pleiotropic morphological traits. We also identified cis-expression quantitative trait loci (QTL) within these dosage QTL regions, a subset of which modulated trans-expression QTL as well. Integration of data types within a gene co-expression framework identified co-expressed gene modules that are dosage sensitive, enriched for dosage expression QTL, and associated with morphological traits. Functional description of these modules linked dosage-sensitive morphological variation to specific cellular processes, as well as candidate regulatory genes. Together, these results show that gene dosage variation can influence morphological variation through complex changes in gene expression, and suggest that frequently occurring gene dosage variation has the potential to likewise influence quantitative traits in nature.

Overcoming the field-of-view to diameter trade-off in microendoscopy via computational optrode-array microscopy.

High-resolution microscopy of deep tissue with large field-of-view (FOV) is critical for elucidating organization of cellular structures in plant biology. Microscopy with an implanted probe offers an effective solution. However, there exists a fundamental trade-off between the FOV and probe diameter arising from aberrations inherent in conventional imaging optics (typically, FOV < 30% of diameter). Here, we demonstrate the use of microfabricated non-imaging probes (optrodes) that when combined with a trained machine-learning algorithm is able to achieve FOV of 1x to 5x the probe diameter. Further increase in FOV is achieved by using multiple optrodes in parallel. With a 1 × 2 optrode array, we demonstrate imaging of fluorescent beads (including 30 FPS video), stained plant stem sections and stained living stems. Our demonstration lays the foundation for fast, high-resolution microscopy with large FOV in deep tissue via microfabricated non-imaging probes and advanced machine learning.

A comprehensive genomic scan reveals gene dosage balance impacts on quantitative traits in Populus trees.

Gene dosage variation and the associated changes in gene expression influence a wide variety of traits, ranging from cancer in humans to yield in plants. It is also expected to affect important traits of ecological and agronomic importance in forest trees, but this variation has not been systematically characterized or exploited. Here we performed a comprehensive scan of the Populus genome for dosage-sensitive loci affecting quantitative trait variation for spring and fall phenology and biomass production. The study population was a large collection of clonally propagated F1 hybrid lines of Populus that saturate the genome 10-fold with deletions and insertions (indels) of known sizes and positions. As a group, the phenotypic means of the indel lines consistently differed from control nonindel lines, with an overall negative effect of both insertions and deletions on all biomass-related traits but more diverse effects and an overall wider phenotypic distribution of the indel lines for the phenology-related traits. We also investigated the correlation between gene dosage at specific chromosomal locations and phenotype, to identify dosage quantitative trait loci (dQTL). Such dQTL were detected for most phenotypes examined, but stronger effect dQTL were identified for the phenology-related traits than for the biomass traits. Our genome-wide screen for dosage sensitivity in a higher eukaryote demonstrates the importance of global genomic balance and the impact of dosage on life history traits.

Evolutionary network genomics of wood formation in a phylogenetic survey of angiosperm forest trees.

Wood formation was present in early angiosperms, but has been highly modified through evolution to generate the anatomical diversity seen in extant angiosperm lineages. In this project, we modeled changes in gene coexpression relationships associated with the evolution of wood formation in a phylogenetic survey of 13 angiosperm tree species. Gravitropic stimulation was used as an experimental treatment to alter wood formation and also perturb gene expression. Gene transcript abundances were determined using RNA sequencing of developing wood tissues from upright trees, and from the top (tension wood) and bottom (opposite wood) tissues of gravistimulated trees. A network-based approach was employed to align gene coexpression networks across species based on orthologous relationships. A large-scale, multilayer network was modeled that identified both lineage-specific gene coexpression modules and modules conserved across multiple species. Functional annotation and analysis of modules identified specific regulatory processes associated with conserved modules, including regulation of hormones, protein phosphorylation, meristem development and epigenetic processes. Our results provide novel insights into the evolution and development of wood formation, and demonstrate the ability to identify biological processes and genes important for the evolution of a foundational trait in nonmodel, undomesticated forest trees.

Brassinosteroid regulation of wood formation in poplar.

Brassinosteroids have been implicated in the differentiation of vascular cell types in herbaceous plants, but their roles during secondary growth and wood formation are not well defined. Here we pharmacologically and genetically manipulated brassinosteroid levels in poplar trees and assayed the effects on secondary growth and wood formation, and on gene expression within stems. Elevated brassinosteroid levels resulted in increases in secondary growth and tension wood formation, while inhibition of brassinosteroid synthesis resulted in decreased growth and secondary vascular differentiation. Analysis of gene expression showed that brassinosteroid action is positively associated with genes involved in cell differentiation and cell-wall biosynthesis. The results presented here show that brassinosteroids play a foundational role in the regulation of secondary growth and wood formation, in part through the regulation of cell differentiation and secondary cell wall biosynthesis.

The auxin receptor TIR1 homolog (PagFBL 1) regulates adventitious rooting through interactions with Aux/IAA28 in Populus.

Adventitious roots occur naturally in many species and can also be induced from explants of some tree species including Populus, providing an important means of clonal propagation. Auxin has been identified as playing a crucial role in adventitious root formation, but the associated molecular regulatory mechanisms need to be elucidated. In this study, we examined the role of PagFBL1, the hybrid poplar (Populus alba × P. glandulosa clone 84K) homolog of Arabidopsis auxin receptor TIR1, in adventitious root formation in poplar. Similar to the distribution pattern of auxin during initiation of adventitious roots, PagFBL1 expression was concentrated in the cambium and secondary phloem in stems during adventitious root induction and initiation phases, but decreased in emerging adventitious root primordia. Overexpressing PagFBL1 stimulated adventitious root formation and increased root biomass, while knock-down of PagFBL1 transcript levels delayed adventitious root formation and decreased root biomass. Transcriptome analyses of PagFBL1 overexpressing lines indicated that an extensive remodelling of gene expression was stimulated by auxin signalling pathway during early adventitious root formation. In addition, PagIAA28 was identified as downstream targets of PagFBL1. We propose that the PagFBL1-PagIAA28 module promotes adventitious rooting and could be targeted to improve Populus propagation by cuttings.

Mosaic modularity: an updated perspective and research agenda for the evolution of vascular cambial growth.

Secondary growth from a vascular cambium, present today only in seed plants and isoetalean lycophytes, has a 400-million-yr evolutionary history that involves considerably broader taxonomic diversity, most of it hidden in the fossil record. Approaching vascular cambial growth as a complex developmental process, we review data from living plants and fossils that reveal diverse modes of secondary growth. These are consistent with a modular nature of secondary growth, when considered as a tracheophyte-wide structural feature. This modular perspective identifies putative constituent developmental modules of cambial growth, for which we review developmental anatomy and regulation. Based on these data, we propose a hypothesis that explains the sources of diversity of secondary growth, considered across the entire tracheophyte clade, and opens up new avenues for exploring the origin of secondary growth. In this hypothesis, various modes of secondary growth reflect a mosaic pattern of expression of different developmental-regulatory modules among different lineages. We outline an approach that queries three information systems (living seed plants, living seed-free plants, and fossils) and integrates data on developmental regulation, anatomy, gene evolution and phylogeny to test the mosaic modularity hypothesis and its implications, and to inform efforts aimed at understanding the evolution of secondary growth.

A System for Dosage-Based Functional Genomics in Poplar.

Altering gene dosage through variation in gene copy number is a powerful approach to addressing questions regarding gene regulation, quantitative trait loci, and heterosis, but one that is not easily applied to sexually transmitted species. Elite poplar (Populus spp) varieties are created through interspecific hybridization, followed by clonal propagation. Altered gene dosage relationships are believed to contribute to hybrid performance. Clonal propagation allows for replication and maintenance of meiotically unstable ploidy or structural variants and provides an alternative approach to investigating gene dosage effects not possible in sexually propagated species. Here, we built a genome-wide structural variation system for dosage-based functional genomics and breeding of poplar. We pollinated Populus deltoides with gamma-irradiated Populus nigra pollen to produce >500 F1 seedlings containing dosage lesions in the form of deletions and insertions of chromosomal segments (indel mutations). Using high-precision dosage analysis, we detected indel mutations in ∼55% of the progeny. These indels varied in length, position, and number per individual, cumulatively tiling >99% of the genome, with an average of 10 indels per gene. Combined with future phenotype and transcriptome data, this population will provide an excellent resource for creating and characterizing dosage-based variation in poplar, including the contribution of dosage to quantitative traits and heterosis.

Transcriptional and Hormonal Regulation of Gravitropism of Woody Stems in Populus.

Angiosperm trees reorient their woody stems by asymmetrically producing a specialized xylem tissue, tension wood, which exerts a strong contractile force resulting in negative gravitropism of the stem. Here, we show, in Populus trees, that initial gravity perception and response occurs in specialized cells through sedimentation of starch-filled amyloplasts and relocalization of the auxin transport protein, PIN3. Gibberellic acid treatment stimulates the rate of tension wood formation and gravibending and enhances tissue-specific expression of an auxin-responsive reporter. Gravibending, maturation of contractile fibers, and gibberellic acid (GA) stimulation of tension wood formation are all sensitive to transcript levels of the Class I KNOX homeodomain transcription factor-encoding gene ARBORKNOX2 (ARK2). We generated genome-wide transcriptomes for trees in which gene expression was perturbed by gravistimulation, GA treatment, and modulation of ARK2 expression. These data were employed in computational analyses to model the transcriptional networks underlying wood formation, including identification and dissection of gene coexpression modules associated with wood phenotypes, GA response, and ARK2 binding to genes within modules. We propose a model for gravitropism in the woody stem in which the peripheral location of PIN3-expressing cells relative to the cambium results in auxin transport toward the cambium in the top of the stem, triggering tension wood formation, while transport away from the cambium in the bottom of the stem triggers opposite wood formation.

Transcriptional and temporal response of Populus stems to gravi-stimulation.

Plants modify development in response to external stimuli, to produce new growth that is appropriate for environmental conditions. For example, gravi-stimulation of leaning branches in angiosperm trees results in modifications of wood development, to produce tension wood that pulls leaning stems upright. Here, we use gravi-stimulation and tension wood response to dissect the temporal changes in gene expression underlying wood formation in Populus stems. Using time-series analysis of seven time points over a 14-d experiment, we identified 8,919 genes that were differentially expressed between tension wood (upper) and opposite wood (lower) sides of leaning stems. Clustering of differentially expressed genes showed four major transcriptional responses, including gene clusters whose transcript levels were associated with two types of tissue-specific impulse responses that peaked at about 24-48 h, and gene clusters with sustained changes in transcript levels that persisted until the end of the 14-d experiment. Functional enrichment analysis of those clusters suggests they reflect temporal changes in pathways associated with hormone regulation, protein localization, cell wall biosynthesis and epigenetic processes. Time-series analysis of gene expression is an underutilized approach for dissecting complex developmental responses in plants, and can reveal gene clusters and mechanisms influencing development.

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.

Transcript profiling of a novel plant meristem, the monocot cambium.

While monocots lack the ability to produce a vascular cambium or woody growth, some monocot lineages evolved a novel lateral meristem, the monocot cambium, which supports secondary radial growth of stems. In contrast to the vascular cambium found in woody angiosperm and gymnosperm species, the monocot cambium produces secondary vascular bundles, which have an amphivasal organization of tracheids encircling a central strand of phloem. Currently there is no information concerning the molecular genetic basis of the development or evolution of the monocot cambium. Here we report high-quality transcriptomes for monocot cambium and early derivative tissues in two monocot genera, Yucca and Cordyline. Monocot cambium transcript profiles were compared to those of vascular cambia and secondary xylem tissues of two forest tree species, Populus trichocarpa and Eucalyptus grandis. Monocot cambium transcript levels showed that there are extensive overlaps between the regulation of monocot cambia and vascular cambia. Candidate regulatory genes that vary between the monocot and vascular cambia were also identified, and included members of the KANADI and CLE families involved in polarity and cell-cell signaling, respectively. We suggest that the monocot cambium may have evolved in part through reactivation of genetic mechanisms involved in vascular cambium regulation.

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

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

Gravitropisms and reaction woods of forest trees - evolution, functions and mechanisms.

Contents 790 I. 790 II. 792 III. 795 IV. 797 V. 798 VI. 800 VII. 800 800 References 800 SUMMARY: The woody stems of trees perceive gravity to determine their orientation, and can produce reaction woods to reinforce or change their position. Together, graviperception and reaction woods play fundamental roles in tree architecture, posture control, and reorientation of stems displaced by wind or other environmental forces. Angiosperms and gymnosperms have evolved strikingly different types of reaction wood. Tension wood of angiosperms creates strong tensile force to pull stems upward, while compression wood of gymnosperms creates compressive force to push stems upward. In this review, the general features and evolution of tension wood and compression wood are presented, along with descriptions of how gravitropisms and reaction woods contribute to the survival and morphology of trees. An overview is presented of the molecular and genetic mechanisms underlying graviperception, initial graviresponse and the regulation of tension wood development in the model angiosperm, Populus. Critical research questions and new approaches are discussed.

Genome-wide mapping of histone H3 lysine 4 trimethylation in Eucalyptus grandis developing xylem.

BACKGROUND: Histone modifications play an integral role in plant development, but have been poorly studied in woody plants. Investigating chromatin organization in wood-forming tissue and its role in regulating gene expression allows us to understand the mechanisms underlying cellular differentiation during xylogenesis (wood formation) and identify novel functional regions in plant genomes. However, woody tissue poses unique challenges for using high-throughput chromatin immunoprecipitation (ChIP) techniques for studying genome-wide histone modifications in vivo. We investigated the role of the modified histone H3K4me3 (trimethylated lysine 4 of histone H3) in gene expression during the early stages of wood formation using ChIP-seq in Eucalyptus grandis, a woody biomass model. RESULTS: Plant chromatin fixation and isolation protocols were optimized for developing xylem tissue collected from field-grown E. grandis trees. A "nano-ChIP-seq" procedure was employed for ChIP DNA amplification. Over 9 million H3K4me3 ChIP-seq and 18 million control paired-end reads were mapped to the E. grandis reference genome for peak-calling using Model-based Analysis of ChIP-Seq. The 12,177 significant H3K4me3 peaks identified covered ~1.5% of the genome and overlapped some 9,623 protein-coding genes and 38 noncoding RNAs. H3K4me3 library coverage, peaking ~600 - 700 bp downstream of the transcription start site, was highly correlated with gene expression levels measured with RNA-seq. Overall, H3K4me3-enriched genes tended to be less tissue-specific than unenriched genes and were overrepresented for general cellular metabolism and development gene ontology terms. Relative expression of H3K4me3-enriched genes in developing secondary xylem was higher than unenriched genes, however, and highly expressed secondary cell wall-related genes were enriched for H3K4me3 as validated using ChIP-qPCR. CONCLUSIONS: In this first genome-wide analysis of a modified histone in a woody tissue, we optimized a ChIP-seq procedure suitable for field-collected samples. In developing E. grandis xylem, H3K4me3 enrichment is an indicator of active transcription, consistent with its known role in sustaining pre-initiation complex formation in yeast. The H3K4me3 ChIP-seq data from this study paves the way to understanding the chromatin landscape and epigenomic architecture of xylogenesis in plants, and complements RNA-seq evidence of gene expression for the future improvement of the E. grandis genome annotation.

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

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

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

BACKGROUND: One of the great advantages of next generation sequencing is the ability to generate large genomic datasets for virtually all species, including non-model organisms. It should be possible, in turn, to apply advanced computational approaches to these datasets to develop models of biological processes. In a practical sense, working with non-model organisms presents unique challenges. In this paper we discuss some of these challenges for ChIP-seq and RNA-seq experiments using the undomesticated tree species of the genus Populus. RESULTS: We describe specific challenges associated with experimental design in Populus, including selection of optimal genotypes for different technical approaches and development of antibodies against Populus transcription factors. Execution of the experimental design included the generation and analysis of Chromatin immunoprecipitation-sequencing (ChIP-seq) data for RNA polymerase II and transcription factors involved in wood formation. We discuss criteria for analyzing the resulting datasets, determination of appropriate control sequencing libraries, evaluation of sequencing coverage needs, and optimization of parameters. We also describe the evaluation of ChIP-seq data from Populus, and discuss the comparison between ChIP-seq and RNA-seq data and biological interpretations of these comparisons. CONCLUSIONS: These and other "lessons learned" highlight the challenges but also the potential insights to be gained from extending next generation sequencing-supported network analyses to undomesticated non-model species.

Modeling transcriptional networks regulating secondary growth and wood formation in forest trees.

The complex interactions among the genes that underlie a biological process can be modeled and presented as a transcriptional network, in which genes (nodes) and their interactions (edges) are shown in a graphical form similar to a wiring diagram. A large number of genes have been identified that are expressed during the radial woody growth of tree stems (secondary growth), but a comprehensive understanding of how these genes interact to influence woody growth is currently lacking. Modeling transcriptional networks has recently been made tractable by next-generation sequencing-based technologies that can comprehensively catalog gene expression and transcription factor-binding genome-wide, but has not yet been extensively applied to undomesticated tree species or woody growth. Here we discuss basic features of transcriptional networks, approaches for modeling biological networks, and examples of biological network models developed for forest trees to date. We discuss how transcriptional network research is being developed in the model forest tree genus, Populus, and how this research area can be further developed and applied. Transcriptional network models for forest tree secondary growth and wood formation could ultimately provide new predictive models to accelerate hypothesis-driven research and develop new breeding applications.