Lee E. Gunter

FLOWERING LOCUS T duplication coordinates reproductive and vegetative growth in perennial poplar

Annual plants grow vegetatively at early developmental stages and then transition to the reproductive stage, followed by senescence in the same year. In contrast, after successive years of vegetative growth at early ages, woody perennial shoot meristems begin repeated transitions between vegetative and reproductive growth at sexual maturity. However, it is unknown how these repeated transitions occur without a developmental conflict between vegetative and reproductive growth. We report that functionally diverged paralogs FLOWERING LOCUS T1 (FT1) and FLOWERING LOCUS T2 (FT2), products of whole-genome duplication and homologs of Arabidopsis thaliana gene FLOWERING LOCUS T (FT), coordinate the repeated cycles of vegetative and reproductive growth in woody perennial poplar (Populus spp.). Our manipulative physiological and genetic experiments coupled with field studies, expression profiling, and network analysis reveal that reproductive onset is determined by FT1 in response to winter temperatures, whereas vegetative growth and inhibition of bud set are promoted by FT2 in response to warm temperatures and long days in the growing season. The basis for functional differentiation between FT1 and FT2 appears to be expression pattern shifts, changes in proteins, and divergence in gene regulatory networks. Thus, temporal separation of reproductive onset and vegetative growth into different seasons via FT1 and FT2 provides seasonality and demonstrates the evolution of a complex perennial adaptive trait after genome duplication.

Genome structure and emerging evidence of an incipient sex chromosome in Populus

The genus Populus consists of dioecious woody species with largely unknown genetic mechanisms for gender determination. We have discovered genetic and genomic features in the peritelomeric region of chromosome XIX that suggest this region of the Populus genome is in the process of developing characteristics of a sex chromosome. We have identified a gender-associated locus that consistently maps to this region. Furthermore, comparison of genetic maps across multiple Populus families reveals consistently distorted segregation within this region. We have intensively characterized this region using an F(1) interspecific cross involving the female genotype that was used for genome sequencing. This region shows suppressed recombination and high divergence between the alternate haplotypes, as revealed by dense map-based genome assembly using microsatellite markers. The suppressed recombination, distorted segregation, and haplotype divergence were observed only for the maternal parent in this cross. Furthermore, the progeny of this cross showed a strongly male-biased sex ratio, in agreement with Haldanes rule that postulates that the heterogametic sex is more likely to be absent, rare, or sterile in interspecific crosses. Together, these results support the role of chromosome XIX in sex determination and suggest that sex determination in Populus occurs through a ZW system in which the female is the heterogametic gender.

Antisense Down-Regulation of 4CL Expression Alters Lignification, Tree Growth, and Saccharification Potential of Field-Grown Poplar

Transgenic down-regulation of the Pt4CL1 gene family encoding 4-coumarate:coenzyme A ligase (4CL) has been reported as a means for reducing lignin content in cell walls and increasing overall growth rates, thereby improving feedstock quality for paper and bioethanol production. Using hybrid poplar (Populus tremula × Populus alba), we applied this strategy and examined field-grown transformants for both effects on wood biochemistry and tree productivity. The reductions in lignin contents obtained correlated well with 4CL RNA expression, with a sharp decrease in lignin amount being observed for RNA expression below approximately 50% of the nontransgenic control. Relatively small lignin reductions of approximately 10% were associated with reduced productivity, decreased wood syringyl/guaiacyl lignin monomer ratios, and a small increase in the level of incorporation of H-monomers (p-hydroxyphenyl) into cell walls. Transgenic events with less than approximately 50% 4CL RNA expression were characterized by patches of reddish-brown discolored wood that had approximately twice the extractive content of controls (largely complex polyphenolics). There was no evidence that substantially reduced lignin contents increased growth rates or saccharification potential. Our results suggest that the capacity for lignin reduction is limited; below a threshold, large changes in wood chemistry and plant metabolism were observed that adversely affected productivity and potential ethanol yield. They also underline the importance of field studies to obtain physiologically meaningful results and to support technology development with transgenic trees.

Population genomics of Populus trichocarpa identifies signatures of selection and adaptive trait associations

Forest trees are dominant components of terrestrial ecosystems that have global ecological and economic importance. Despite distributions that span wide environmental gradients, many tree populations are locally adapted, and mechanisms underlying this adaptation are poorly understood. Here we use a combination of whole-genome selection scans and association analyses of 544 Populus trichocarpa trees to reveal genomic bases of adaptive variation across a wide latitudinal range. Three hundred ninety-seven genomic regions showed evidence of recent positive and/or divergent selection and enrichment for associations with adaptive traits that also displayed patterns consistent with natural selection. These regions also provide unexpected insights into the evolutionary dynamics of duplicated genes and their roles in adaptive trait variation.

Genome resequencing reveals multiscale geographic structure and extensive linkage disequilibrium in the forest tree Populus trichocarpa

• Plant population genomics informs evolutionary biology, breeding, conservation and bioenergy feedstock development. For example, the detection of reliable phenotype-genotype associations and molecular signatures of selection requires a detailed knowledge about genome-wide patterns of allele frequency variation, linkage disequilibrium and recombination. • We resequenced 16 genomes of the model tree Populus trichocarpa and genotyped 120 trees from 10 subpopulations using 29,213 single-nucleotide polymorphisms. • Significant geographic differentiation was present at multiple spatial scales, and range-wide latitudinal allele frequency gradients were strikingly common across the genome. The decay of linkage disequilibrium with physical distance was slower than expected from previous studies in Populus, with r(2) dropping below 0.2 within 3-6 kb. Consistent with this, estimates of recent effective population size from linkage disequilibrium (N(e) ≈ 4000-6000) were remarkably low relative to the large census sizes of P. trichocarpa stands. Fine-scale rates of recombination varied widely across the genome, but were largely predictable on the basis of DNA sequence and methylation features. • Our results suggest that genetic drift has played a significant role in the recent evolutionary history of P. trichocarpa. Most importantly, the extensive linkage disequilibrium detected suggests that genome-wide association studies and genomic selection in undomesticated populations may be more feasible in Populus than previously assumed.

A Multifactor Analysis of Fungal and Bacterial Community Structure in the Root Microbiome of Mature Populus deltoides Trees

Bacterial and fungal communities associated with plant roots are central to the host health, survival and growth. However, a robust understanding of the root-microbiome and the factors that drive host associated microbial community structure have remained elusive, especially in mature perennial plants from natural settings. Here, we investigated relationships of bacterial and fungal communities in the rhizosphere and root endosphere of the riparian tree species Populus deltoides, and the influence of soil parameters, environmental properties (host phenotype and aboveground environmental settings), host plant genotype (Simple Sequence Repeat (SSR) markers), season (Spring vs. Fall) and geographic setting (at scales from regional watersheds to local riparian zones) on microbial community structure. Each of the trees sampled displayed unique aspects to its associated community structure with high numbers of Operational Taxonomic Units (OTUs) specific to an individual trees (bacteria >90%, fungi >60%). Over the diverse conditions surveyed only a small number of OTUs were common to all samples within rhizosphere (35 bacterial and 4 fungal) and endosphere (1 bacterial and 1 fungal) microbiomes. As expected, Proteobacteria and Ascomycota were dominant in root communities (>50%) while other higher-level phylogenetic groups (Chytridiomycota, Acidobacteria) displayed greatly reduced abundance in endosphere compared to the rhizosphere. Variance partitioning partially explained differences in microbiome composition between all sampled roots on the basis of seasonal and soil properties (4% to 23%). While most variation remains unattributed, we observed significant differences in the microbiota between watersheds (Tennessee vs. North Carolina) and seasons (Spring vs. Fall). SSR markers clearly delineated two host populations associated with the samples taken in TN vs. NC, but overall host genotypic distances did not have a significant effect on corresponding communities that could be separated from other measured effects.

Connecting genes, coexpression modules, and molecular signatures to environmental stress phenotypes in plants.

BackgroundOne of the eminent opportunities afforded by modern genomic technologies is the potential to provide a mechanistic understanding of the processes by which genetic change translates to phenotypic variation and the resultant appearance of distinct physiological traits. Indeed much progress has been made in this area, particularly in biomedicine where functional genomic information can be used to determine the physiological state (e.g., diagnosis) and predict phenotypic outcome (e.g., patient survival). Ecology currently lacks an analogous approach where genomic information can be used to diagnose the presence of a given physiological state (e.g., stress response) and then predict likely phenotypic outcomes (e.g., stress duration and tolerance, fitness).ResultsHere, we demonstrate that a compendium of genomic signatures can be used to classify the plant abiotic stress phenotype in Arabidopsis according to the architecture of the transcriptome, and then be linked with gene coexpression network analysis to determine the underlying genes governing the phenotypic response. Using this approach, we confirm the existence of known stress responsive pathways and marker genes, report a common abiotic stress responsive transcriptome and relate phenotypic classification to stress duration.ConclusionLinking genomic signatures to gene coexpression analysis provides a unique method of relating an observed plant phenotype to changes in gene expression that underlie that phenotype. Such information is critical to current and future investigations in plant biology and, in particular, to evolutionary ecology, where a mechanistic understanding of adaptive physiological responses to abiotic stress can provide researchers with a tool of great predictive value in understanding species and population level adaptation to climate change.

The F-Box Gene Family Is Expanded in Herbaceous Annual Plants Relative to Woody Perennial Plants

F-box proteins are generally responsible for substrate recognition in the Skp1-Cullin-F-box complexes that are involved in protein degradation via the ubiquitin-26S proteasome pathway. In plants, F-box genes influence a variety of biological processes, such as leaf senescence, branching, self-incompatibility, and responses to biotic and abiotic stresses. The number of F-box genes in Populus (Populus trichocarpa; approximately 320) is less than half that found in Arabidopsis (Arabidopsis thaliana; approximately 660) or Oryza (Oryza sativa; approximately 680), even though the total number of genes in Populus is equivalent to that in Oryza and 1.5 times that in Arabidopsis. We performed comparative genomics analysis between the woody perennial plant Populus and the herbaceous annual plants Arabidopsis and Oryza in order to explicate the functional implications of this large gene family. Our analyses reveal interspecific differences in genomic distribution, orthologous relationship, intron evolution, protein domain structure, and gene expression. The set of F-box genes shared by these species appear to be involved in core biological processes essential for plant growth and development; lineage-specific differences primarily occurred because of an expansion of the F-box genes via tandem duplications in Arabidopsis and Oryza. The number of F-box genes in the newly sequenced woody species Vitis (Vitis vinifera; 156) and Carica (Carica papaya; 139) is similar to that in Populus, supporting the hypothesis that the F-box gene family is expanded in herbaceous annual plants relative to woody perennial plants. This study provides insights into the relationship between the structure and composition of the F-box gene family in herbaceous and woody species and their associated developmental and physiological features.

Discovery and annotation of small proteins using genomics, proteomics, and computational approaches

Small proteins (10-200 amino acids [aa] in length) encoded by short open reading frames (sORF) play important regulatory roles in various biological processes, including tumor progression, stress response, flowering, and hormone signaling. However, ab initio discovery of small proteins has been relatively overlooked. Recent advances in deep transcriptome sequencing make it possible to efficiently identify sORFs at the genome level. In this study, we obtained ~2.6 million expressed sequence tag (EST) reads from Populus deltoides leaf transcriptome and reconstructed full-length transcripts from the EST sequences. We identified an initial set of 12,852 sORFs encoding proteins of 10-200 aa in length. Three computational approaches were then used to enrich for bona fide protein-coding sORFs from the initial sORF set: (1) coding-potential prediction, (2) evolutionary conservation between P. deltoides and other plant species, and (3) gene family clustering within P. deltoides. As a result, a high-confidence sORF candidate set containing 1469 genes was obtained. Analysis of the protein domains, non-protein-coding RNA motifs, sequence length distribution, and protein mass spectrometry data supported this high-confidence sORF set. In the high-confidence sORF candidate set, known protein domains were identified in 1282 genes (higher-confidence sORF candidate set), out of which 611 genes, designated as highest-confidence candidate sORF set, were supported by proteomics data. Of the 611 highest-confidence candidate sORF genes, 56 were new to the current Populus genome annotation. This study not only demonstrates that there are potential sORF candidates to be annotated in sequenced genomes, but also presents an efficient strategy for discovery of sORFs in species with no genome annotation yet available.

Analysis of 4,664 high-quality sequence-finished poplar full-length cDNA clones and their utility for the discovery of genes responding to insect feeding

BackgroundThe genus Populus includes poplars, aspens and cottonwoods, which will be collectively referred to as poplars hereafter unless otherwise specified. Poplars are the dominant tree species in many forest ecosystems in the Northern Hemisphere and are of substantial economic value in plantation forestry. Poplar has been established as a model system for genomics studies of growth, development, and adaptation of woody perennial plants including secondary xylem formation, dormancy, adaptation to local environments, and biotic interactions.ResultsAs part of the poplar genome sequencing project and the development of genomic resources for poplar, we have generated a full-length (FL)-cDNA collection using the biotinylated CAP trapper method. We constructed four FLcDNA libraries using RNA from xylem, phloem and cambium, and green shoot tips and leaves from the P. trichocarpa Nisqually-1 genotype, as well as insect-attacked leaves of the P. trichocarpa × P. deltoides hybrid. Following careful selection of candidate cDNA clones, we used a combined strategy of paired end reads and primer walking to generate a set of 4,664 high-accuracy, sequence-verified FLcDNAs, which clustered into 3,990 putative unique genes. Mapping FLcDNAs to the poplar genome sequence combined with BLAST comparisons to previously predicted protein coding sequences in the poplar genome identified 39 FLcDNAs that likely localize to gaps in the current genome sequence assembly. Another 173 FLcDNAs mapped to the genome sequence but were not included among the previously predicted genes in the poplar genome. Comparative sequence analysis against Arabidopsis thaliana and other species in the non-redundant database of GenBank revealed that 11.5% of the poplar FLcDNAs display no significant sequence similarity to other plant proteins. By mapping the poplar FLcDNAs against transcriptome data previously obtained with a 15.5 K cDNA microarray, we identified 153 FLcDNA clones for genes that were differentially expressed in poplar leaves attacked by forest tent caterpillars.ConclusionThis study has generated a high-quality FLcDNA resource for poplar and the third largest FLcDNA collection published to date for any plant species. We successfully used the FLcDNA sequences to reassess gene prediction in the poplar genome sequence, perform comparative sequence annotation, and identify differentially expressed transcripts associated with defense against insects. The FLcDNA sequences will be essential to the ongoing curation and annotation of the poplar genome, in particular for targeting gaps in the current genome assembly and further improvement of gene predictions. The physical FLcDNA clones will serve as useful reagents for functional genomics research in areas such as analysis of gene functions in defense against insects and perennial growth. Sequences from this study have been deposited in NCBI GenBank under the accession numbers EF144175 to EF148838.