The uncharacterized gene EVE contributes to vessel element dimensions in Populus.

The radiation of angiosperms led to the emergence of the vast majority of today's plant species and all our major food crops. Their extraordinary diversification occurred in conjunction with the evolution of a more efficient vascular system for the transport of water, composed of vessel elements. The physical dimensions of these water-conducting specialized cells have played a critical role in angiosperm evolution; they determine resistance to water flow, influence photosynthesis rate, and contribute to plant stature. However, the genetic factors that determine their dimensions are unclear. Here we show that a previously uncharacterized gene, ENLARGED VESSEL ELEMENT (EVE), contributes to the dimensions of vessel elements in Populus, impacting hydraulic conductivity. Our data suggest that EVE is localized in the plasma membrane and is involved in potassium uptake of differentiating xylem cells during vessel development. In plants, EVE first emerged in streptophyte algae, but expanded dramatically among vessel-containing angiosperms. The phylogeny, structure and composition of EVE indicates that it may have been involved in an ancient horizontal gene-transfer event.

Integration of genetic, genomic and transcriptomic information identifies putative regulators of adventitious root formation in Populus.

BACKGROUND: Adventitious roots (AR) develop from tissues other than the primary root, in a process physiologically regulated by phytohormones. Adventitious roots provide structural support and contribute to water and nutrient absorption, and are critical for commercial vegetative propagation of several crops. Here we quantified the number of AR, root architectural traits and root biomass in cuttings from a pseudo-backcross population of Populus deltoides and Populus trichocarpa. Quantitative trait loci (QTL) mapping and whole-transcriptome analysis of individuals with alternative QTL alleles for AR number were used to identify putative regulators of AR development. RESULTS: Parental individuals and progeny showed extensive segregation for AR developmental traits. Quantitative trait loci for number of AR mapped consistently in the same interval of linkage group (LG) II and LG XIV, explaining 7-10 % of the phenotypic variation. A time series transcriptome analysis identified 26,121 genes differentially expressed during AR development, particularly during the first 24 h after cuttings were harvested. Of those, 1929 genes were differentially regulated between individuals carrying alternative alleles for the two QTL for number of AR, in one or more time point. Eighty-one of these genes were physically located within the QTL intervals for number of AR, including putative homologs of the Arabidopsis genes SUPERROOT2 (SUR2) and TRYPTOPHAN SYNTHASE ALPHA CHAIN (TSA1), both of which are involved in the auxin indole-3-acetic acid (IAA) biosynthesis pathway. CONCLUSIONS: This study suggests the involvement of two genes of the tryptophan-dependent auxin biosynthesis pathway, SUR2 and TSA1, in the regulation of a critical trait for the clonal propagation of woody species. A possible model for this regulation is that poplar individuals that have poor AR formation synthesize auxin indole-3-acetic acid (IAA) primarily through the tryptophan (Trp) pathway. Much of the Trp pathway flux appears to be directed to the synthesis of indole glucosinolates (IG), as suggested by the over-expression of SUR2. Individuals that are efficient in AR formation may utilize alternative (non-Trp) pathways to synthesize IAA, based on the observation that they down-regulate the expression of TSA1, one of the critical steps in the synthesis of tryptophan.

Genetical genomics of Populus leaf shape variation.

BACKGROUND: Leaf morphology varies extensively among plant species and is under strong genetic control. Mutagenic screens in model systems have identified genes and established molecular mechanisms regulating leaf initiation, development, and shape. However, it is not known whether this diversity across plant species is related to naturally occurring variation at these genes. Quantitative trait locus (QTL) analysis has revealed a polygenic control for leaf shape variation in different species suggesting that loci discovered by mutagenesis may only explain part of the naturally occurring variation in leaf shape. Here we undertook a genetical genomics study in a poplar intersectional pseudo-backcross pedigree to identify genetic factors controlling leaf shape. The approach combined QTL discovery in a genetic linkage map anchored to the Populus trichocarpa reference genome sequence and transcriptome analysis. RESULTS: A major QTL for leaf lamina width and length:width ratio was identified in multiple experiments that confirmed its stability. A transcriptome analysis of expanding leaf tissue contrasted gene expression between individuals with alternative QTL alleles, and identified an ADP-ribosylation factor (ARF) GTPase (PtARF1) as a candidate gene for regulating leaf morphology in this pedigree. ARF GTPases are critical elements in the vesicular trafficking machinery. Disruption of the vesicular trafficking function of ARF by the pharmacological agent Brefeldin A (BFA) altered leaf lateral growth in the narrow-leaf P. trichocarpa suggesting a molecular mechanism of leaf shape determination. Inhibition of the vesicular trafficking processes by BFA interferes with cycling of PIN proteins and causes their accumulation in intercellular compartments abolishing polar localization and disrupting normal auxin flux with potential effects on leaf expansion. CONCLUSIONS: In other model systems, ARF proteins have been shown to control the localization of auxin efflux carriers, which function to establish auxin gradients and apical-basal cell polarity in developing plant organs. Our results support a model where PtARF1 transcript abundance changes the dynamics of endocytosis-mediated PIN localization in leaf cells, thus affecting lateral auxin flux and subsequently lamina leaf expansion. This suggests that evolution of differential cellular polarity plays a significant role in leaf morphological variation observed in subgenera of genus Populus.

Diversification in the genetic architecture of gene expression and transcriptional networks in organ differentiation of Populus.

A fundamental goal of systems biology is to identify genetic elements that contribute to complex phenotypes and to understand how they interact in networks predictive of system response to genetic variation. Few studies in plants have developed such networks, and none have examined their conservation among functionally specialized organs. Here we used genetical genomics in an interspecific hybrid population of the model hardwood plant Populus to uncover transcriptional networks in xylem, leaves, and roots. Pleiotropic eQTL hotspots were detected and used to construct coexpression networks a posteriori, for which regulators were predicted based on cis-acting expression regulation. Networks were shown to be enriched for groups of genes that function in biologically coherent processes and for cis-acting promoter motifs with known roles in regulating common groups of genes. When contrasted among xylem, leaves, and roots, transcriptional networks were frequently conserved in composition, but almost invariably regulated by different loci. Similarly, the genetic architecture of gene expression regulation is highly diversified among plant organs, with less than one-third of genes with eQTL detected in two organs being regulated by the same locus. However, colocalization in eQTL position increases to 50% when they are detected in all three organs, suggesting conservation in the genetic regulation is a function of ubiquitous expression. Genes conserved in their genetic regulation among all organs are primarily cis regulated (approximately 92%), whereas genes with eQTL in only one organ are largely trans regulated. Trans-acting regulation may therefore be the primary driver of differentiation in function between plant organs.

Quantitative genetic analysis of biomass and wood chemistry of Populus under different nitrogen levels.

The genetic control of carbon allocation and partitioning in woody perennial plants is poorly understood despite its importance for carbon sequestration, biofuels and other wood-based industries. It is also unclear how environmental cues, such as nitrogen availability, impact the genes that regulate growth, biomass allocation and wood composition in trees. We phenotyped 396 clonally replicated genotypes of an interspecific pseudo-backcross pedigree of Populus for wood composition and biomass traits in above- and below-ground organs. The loci that regulate growth, carbon allocation and partitioning under two nitrogen conditions were identified, defining the contribution of environmental cues to their genetic control. Sixty-three quantitative trait loci were identified for the 20 traits analyzed. The majority of quantitative trait loci are specific to one of the two nitrogen treatments, demonstrating significant nitrogen-dependent genetic control. A highly significant genetic correlation was observed between plant growth and lignin/cellulose composition, and quantitative trait loci co-localization identified the genomic position of potential pleiotropic regulators. Pleiotropic loci linking higher growth rates to wood with less lignin are excellent targets to engineer tree germplasm improved for pulp, paper and cellulosic ethanol production. The causative genes are being identified with a genetical genomics approach.

A microarray-based genotyping and genetic mapping approach for highly heterozygous outcrossing species enables localization of a large fraction of the unassembled Populus trichocarpa genome sequence.

Microarrays have demonstrated significant power for genome-wide analyses of gene expression, and recently have also revolutionized the genetic analysis of segregating populations by genotyping thousands of loci in a single assay. Although microarray-based genotyping approaches have been successfully applied in yeast and several inbred plant species, their power has not been proven in an outcrossing species with extensive genetic diversity. Here we have developed methods for high-throughput microarray-based genotyping in such species using a pseudo-backcross progeny of 154 individuals of Populus trichocarpa and P. deltoides analyzed with long-oligonucleotide in situ-synthesized microarray probes. Our analysis resulted in high-confidence genotypes for 719 single-feature polymorphism (SFP) and 1014 gene expression marker (GEM) candidates. Using these genotypes and an established microsatellite (SSR) framework map, we produced a high-density genetic map comprising over 600 SFPs, GEMs and SSRs. The abundance of gene-based markers allowed us to localize over 35 million base pairs of previously unplaced whole-genome shotgun (WGS) scaffold sequence to putative locations in the genome of P. trichocarpa. A high proportion of sampled scaffolds could be verified for their placement with independently mapped SSRs, demonstrating the previously un-utilized power that high-density genotyping can provide in the context of map-based WGS sequence reassembly. Our results provide a substantial contribution to the continued improvement of the Populus genome assembly, while demonstrating the feasibility of microarray-based genotyping in a highly heterozygous population. The strategies presented are applicable to genetic mapping efforts in all plant species with similarly high levels of genetic diversity.

High-throughput gene and SNP discovery in Eucalyptus grandis, an uncharacterized genome.

BACKGROUND: Benefits from high-throughput sequencing using 454 pyrosequencing technology may be most apparent for species with high societal or economic value but few genomic resources. Rapid means of gene sequence and SNP discovery using this novel sequencing technology provide a set of baseline tools for genome-level research. However, it is questionable how effective the sequencing of large numbers of short reads for species with essentially no prior gene sequence information will support contig assemblies and sequence annotation. RESULTS: With the purpose of generating the first broad survey of gene sequences in Eucalyptus grandis, the most widely planted hardwood tree species, we used 454 technology to sequence and assemble 148 Mbp of expressed sequences (EST). EST sequences were generated from a normalized cDNA pool comprised of multiple tissues and genotypes, promoting discovery of homologues to almost half of Arabidopsis genes, and a comprehensive survey of allelic variation in the transcriptome. By aligning the sequencing reads from multiple genotypes we detected 23,742 SNPs, 83% of which were validated in a sample. Genome-wide nucleotide diversity was estimated for 2,392 contigs using a modified theta (theta) parameter, adapted for measuring genetic diversity from polymorphisms detected by randomly sequencing a multi-genotype cDNA pool. Diversity estimates in non-synonymous nucleotides were on average 4x smaller than in synonymous, suggesting purifying selection. Non-synonymous to synonymous substitutions (Ka/Ks) among 2,001 contigs averaged 0.30 and was skewed to the right, further supporting that most genes are under purifying selection. Comparison of these estimates among contigs identified major functional classes of genes under purifying and diversifying selection in agreement with previous researches. CONCLUSION: In providing an abundance of foundational transcript sequences where limited prior genomic information existed, this work created part of the foundation for the annotation of the E. grandis genome that is being sequenced by the US Department of Energy. In addition we demonstrated that SNPs sampled in large-scale with 454 pyrosequencing can be used to detect evolutionary signatures among genes, providing one of the first genome-wide assessments of nucleotide diversity and Ka/Ks for a non-model plant species.