Search

Despite the economic, ecological, and scientific importance of the generaPopulusL. (poplars, cottonwoods, and aspens) andSalixL. Salicaceae (willows), we know little about the sources of differences in species diversity between the genera and of the phylogenetic conflict that often confounds estimating phylogenetic trees.Salixsubgenera and sections, in particular, have been difficult to classify, with one recent attempt termed a ‘spectacular failure’ due to a speculated radiation of the subgeneraVetrixandChamaetia. Here we use targeted sequence capture to understand the evolutionary history of this portion of the Salicaceae plant family. Our phylogenetic hypothesis was based on 787 gene regions and identified extensive phylogenetic conflict among genes. Our analysis supported some previously described subgeneric relationships and confirmed polyphyly of others. Using an fbranchanalysis we identified several cases of hybridization in deep branches of the phylogeny, which likely contributed to discordance among gene trees. In addition, we identified a rapid increase in diversification rate near the origination of theVetrix-Chamaetiaclade inSalix. This region of the tree coincided with several nodes that lacked strong statistical support, indicating a possible increase in incomplete lineage sorting due to rapid diversification. The extraordinary level of both recent and ancient hybridization in bothPopulusandSalixhave played important roles in the diversification and diversity in these two genera.Supplementary data files will be provided by request tomatt.olson@ttu.edu

PREMISE The evolution of sex chromosomes is driven by sexual dimorphism, yet it can be challenging to document sexually dimorphic traits in dioecious plant species. At the genetic level, sexual dimorphism can be identified through sequence variation between females and males associated with sexually antagonistic traits and different fitness optima. This study aims to examine sexual dimorphism for 26 traits in three populations of Salix purpurea (a diversity panel and F1 and F2 populations) and determine the effect of the traits on biomass yield, a key trait in Salix bioenergy crops across multiple years, locations, and under manipulated growth conditions. METHODS Sexual dimorphism was evaluated for morphological, phenological, physiological, and wood composition traits in a diversity panel of unrelated S. purpurea accessions and in full-sib F1 and F2 families produced through controlled cross pollinations and grown in replicated field trials. RESULTS We observed sexual dimorphism in the timing of development for several traits that were highly predictive of biomass yield across three populations of S. purpurea. Across all populations and years surveyed, males had significantly shallower branching angle. Male plants highly predictive of biomass yield across three populations of S. purpurea also accumulated more nitrogen under fertilizer amendment as measured by SPAD in the diversity panel and had greater susceptibility to the rust fungus Melampsora americana in the F2 family. Allometric modelling of biomass yield showed an effect of sex and of location on the interaction between yield and stem height. CONCLUSIONS These results provide evidence of sexual dimorphism for certain traits in S. purpurea that may be involved in sex chromosome evolution.

16 Pág. Instituto de Ciencias Forestales (ICIFOR), Fine-scale meiotic recombination is fundamental to the outcome of natural and artificial selection. Here, dense genetic mapping and haplotype reconstruction were used to estimate recombination for a full factorial Populus trichocarpa cross of 7 males and 7 females. Genomes of the resulting 49 full-sib families (N = 829 offspring) were resequenced, and high-fidelity biallelic SNP/INDELs and pedigree information were used to ascertain allelic phase and impute progeny genotypes to recover gametic haplotypes. The 14 parental genetic maps contained 1,820 SNP/INDELs on average that covered 376.7 Mb of physical length across 19 chromosomes. Comparison of parental and progeny haplotypes allowed fine-scale demarcation of cross-over regions, where 38,846 cross-over events in 1,658 gametes were observed. Cross-over events were positively associated with gene density and negatively associated with GC content and long-terminal repeats. One of the most striking findings was higher rates of cross-overs in males in 8 out of 19 chromosomes. Regions with elevated male cross-over rates had lower gene density and GC content than windows showing no sex bias. High-resolution analysis identified 67 candidate cross-over hotspots spread throughout the genome. DNA sequence motifs enriched in these regions showed striking similarity to those of maize, Arabidopsis, and wheat. These findings, and recombination estimates, will be useful for ongoing efforts to accelerate domestication of this and other biomass feedstocks, as well as future studies investigating broader questions related to evolutionary history, perennial development, phenology, wood formation, vegetative propagation, and dioecy that cannot be studied using annual plant model systems., This research was supported by the U.S. Department of Energy (DOE), Office of Biological and Environmental Research through the Center for Bioenergy Innovation (CBI), and a DOE Bioenergy Research Center. The publisher, by accepting the article for publication, acknowledges that the U.S. Government retains a nonexclusive, paid-up, irrevocable, worldwide license to publish, or reproduce the published form of this work, or allow others to do so, for U.S. Government purposes. The views expressed in the article do not necessarily represent the views of the U.S. Department of Energy or the United States Government. Work on heterochiasmy was supported by the National Science Foundation (award 1542509 to SD). The breeding of the 7×7 population was initially funded by an award under the Western Sun Grant program to Washington State University (Subgrant T-0013A). Dr Macaya-Sanz was partially supported by “Atracción de Talento Investigador” of the Community of Madrid (ref. 2019-T2/BIO-12780). The project “Conservación y promoción de recursos genéticos forestales contempladas en el programa nacional de desarrollo rural” (ref. IMP2018-002), funded by the Spanish Ministry of Ecological Transition, also covered some expenses of this study.

Plants employ a diverse set of defense mechanisms to mediate interactions with insects and fungi. These relationships can leave lasting impacts on host plant genome structure such as rapid expansion of gene families through tandem duplication. These genomic signatures provide important clues about the complexities of plant/biotic stress interactions and evolution. We used a pseudo‐backcross hybrid family to identify quantitative trait loci (QTL) controlling associations between Populus trees and several common Populus diseases and insects. Using whole‐genome sequences from each parent, we identified candidate genes that may mediate these interactions. Candidates were partially validated using mass spectrometry to identify corresponding QTL for defensive compounds. We detected significant QTL for two interacting fungal pathogens and three insects. The QTL intervals contained candidate genes potentially involved in physical and chemical mechanisms of host–plant resistance and susceptibility. In particular, we identified adjoining QTLs for a phenolic glycoside and Phyllocolpa sawfly abundance. There was also significant enrichment of recent tandem duplications in the genomic intervals of the native parent, but not the exotic parent. Tandem gene duplication may be an important mechanism for rapid response to biotic stressors, enabling trees with long juvenile periods to reach maturity despite many coevolving biotic stressors., This work focuses on elucidating the host genetic drivers of community composition in a hybrid Populus family. We mapped quantitative trait loci for pathogenic fungi and leaf modifying insects in plantations and connected these to genomic intervals of the two host species. We found that tandemly duplicated genes were significantly enriched in the intervals and that these genes have functions that are consistent with rapid host co‐evolution with the interacting organisms.

Populus euphratica is well adapted to extreme desert environments and is an important model species for elucidating the mechanisms of abiotic stress resistance in trees. The current assembly of P. euphratica genome is highly fragmented with many gaps and errors, thereby impeding downstream applications. Here, we report an improved chromosome-level reference genome of P. euphratica (v2.0) using single-molecule sequencing and chromosome conformation capture (Hi-C) technologies. Relative to the previous reference genome, our assembly represents a nearly 60-fold improvement in contiguity, with a scaffold N50 size of 28.59 Mb. Using this genome, we have found that extensive expansion of Gypsy elements in P. euphratica led to its rapid increase in genome size compared to any other Salicaceae species studied to date, and potentially contributed to adaptive divergence driven by insertions near genes involved in stress tolerance. We also detected a wide range of unique structural rearrangements in P. euphratica, including 2,549 translocations, 454 inversions, 121 tandem and 14 segmental duplications. Several key genes likely to be involved in tolerance to abiotic stress were identified within these regions. This high-quality genome represents a valuable resource for poplar breeding and genetic improvement in the future, as well as comparative genomic analysis with other Salicaceae species.

Gene-editing methods, particularly CRISPR, provide extraordinary opportunities for scientific insights and applications in the life sciences. However, the prospects for near-term applications to commercial forestry appear limited. Loss-of-function phenotypes that can be imparted by mutation of one or a few conserved genes offer the best opportunities in the near term. For traits with complex inheritance, there is insufficient science to guide gene-editing efforts, and Genome-Wide Association Studies (GWASs), without strong validation, typically cannot provide high-confidence gene identification. Other obstacles include the difficulty of transformation in many important genotypes, difficulties of transient editing or complete editor removal, and complexity of use in breeding programs. Gene edits that cause loss-of-function traits will generally be recessive, and thus not be expressed among outbred progeny, so vegetative propagules (clones) will be required in most cases. There are also important societal constraints, such as strict regulations for field trials in most countries, and market certification systems that do not allow any kinds of recombinant DNA-modified trees, including those produced by gene-editing, in certified production forests. We conclude that gene-editing applications will be extremely limited for the foreseeable future (i.e., at least 10 years). Nevertheless, gene-editing is a very powerful scientific tool that will be widely used by molecular forest scientists and can lead to important applications in the longer term, if research advances are made on key fronts and regulatory and market obstacles greatly attenuated.

The rapid analysis of biopolymers including lignin and sugars in lignocellulosic biomass cell walls is essential for the analysis of the large sample populations needed for identifying heritable genetic variation in biomass feedstocks for biofuels and bioproducts. In this study, we reported the analysis of cell wall lignin content, syringyl/guaiacyl (S/G) ratio, as well as glucose and xylose content by high-throughput pyrolysis-molecular beam mass spectrometry (py-MBMS) for3,600 samples derived from hundreds of accessions of

Sex dimorphism and gene expression were studied in developing catkins in 159 F2 individuals from the bioenergy crop Salix purpurea, and potential mechanisms and pathways for regulating sex development were explored. Differential expression, eQTL, bisulfite sequencing, and network analysis were used to characterize sex dimorphism, detect candidate master regulator genes, and identify pathways through which the sex determination region (SDR) may mediate sex dimorphism. Eleven genes are presented as candidates for master regulators of sex, supported by gene expression and network analyses. These include genes putatively involved in hormone signaling, epigenetic modification, and regulation of transcription. eQTL analysis revealed a suite of transcription factors and genes involved in secondary metabolism and floral development that were predicted to be under direct control of the sex determination region. Furthermore, data from bisulfite sequencing and small RNA sequencing revealed strong differences in expression between males and females that would implicate both of these processes in sex dimorphism pathways. These data indicate that the mechanism of sex determination in Salix purpurea is likely different from that observed in the related genus Populus. This further demonstrates the dynamic nature of SDRs in plants, which involves a multitude of mechanisms of sex determination and a high rate of turnover.

Phylogenetic analysis is complicated by interspecific gene flow and the presence of shared ancestral polymorphisms, particularly those maintained by balancing selection. In this study, we aimed to examine the prevalence of these factors during the diversification of Populus, a model tree genus in the Northern Hemisphere. We constructed phylogenetic trees of 29 Populus taxa using 80 individuals based on re-sequenced genomes. Our species tree analyses recovered four main clades in the genus based on consensus nuclear phylogenies, but in conflict with the plastome phylogeny. A few interspecific relationships remained unresolved within the multiple-species clade because of inconsistent gene trees. Our results indicated that gene flow has been widespread within each clade and also occurred among the four clades during their early divergence. We identified 45 candidate genes with ancient polymorphisms maintained by balancing selection. These genes were mainly associated with mating compatibility, growth and stress resistance. Both gene flow and selection-mediated ancient polymorphisms are prevalent in the genus Populus. These are potentially important contributors to adaptive variation. Our results provide a framework for the diversification of model tree genus that will facilitate future comparative studies.

Summary • In an attempt to elucidate the molecular mechanisms of Melampsora rust resistance in Populus trichocarpa, we have mapped two resistance loci, MXC3 and MER, and intensively characterized the flanking genomic sequence for the MXC3 locus and the level of linkage disequilibrium (LD) in natural populations. • We used an interspecific backcross pedigree and a genetic map that was highly saturated with AFLP and SSR markers, and assembled shotgun-sequence data in the region containing markers linked to MXC3. • The two loci were mapped to different linkage groups. Linkage disequilibrium for MXC3 was confined to two closely linked regions spanning 34 and 16 kb, respectively. The MXC3 region also contained six disease-resistance candidate genes. • The MER and MXC3 loci are clearly distinct, and may have different mechanisms of resistance, as different classes of putative resistance genes were present near each locus. The suppressed recombination previously observed in the MXC3 region was possibly caused by extensive hemizygous rearrangements confined to the original parent tree. The relatively low observed LD may facilitate association studies using candidate genes for rust resistance, but will probably inhibit marker-aided selection.

Sex dimorphism and gene expression were studied in developing catkins in 159 F2individuals from the bioenergy cropSalix purpurea, and potential mechanisms and pathways for regulating sex development were explored. Differential expression, eQTL, bisulfite sequencing, and network analysis were used to characterize sex dimorphism, detect candidate master regulator genes, and identify pathways through which the sex determination region (SDR) may mediate sex dimorphism. Eleven genes are presented as candidates for master regulators of sex, supported by gene expression and network analyses. These include genes putatively involved in hormone signaling, epigenetic modification, and regulation of transcription. eQTL analysis revealed a suite of transcription factors and genes involved in secondary metabolism and floral development that were predicted to be under direct control of the sex determination region. Furthermore, data from bisulfite sequencing and small RNA sequencing revealed strong differences in expression between males and females that would implicate both of these processes in sex dimorphism pathways. These data indicate that the mechanism of sex determination inSalix purpureais likely different from that observed in the related genusPopulus. This further demonstrates the dynamic nature of SDRs in plants, which involves a multitude of mechanisms of sex determination and a high rate of turnover.

The development of non-recombining sex chromosomes has radical effects on the evolution of discrete sexes and sexual dimorphism. Although dioecy is rare in plants, sex chromosomes have evolved repeatedly throughout the diversification of angiosperms, and many of these sex chromosomes are relatively young compared to those found in vertebrates. In this study, we designed and used a sequence capture array to identify a novel sex-linked region (SLR) in Salix nigra, a basal species in the willow clade, and demonstrated that this species has XY heterogamety. We did not detect any genetic overlap with the previously characterized ZW SLRs in willows, which map to a different chromosome. The S. nigra SLR is characterized by strong recombination suppression across a 2 MB region and an excess of low-frequency alleles, resulting in a low Tajima's D compared to the remainder of the genome. We speculate that either a recent bottleneck in population size or factors related to positive or background selection generated this differential pattern of Tajima's D on the X and autosomes. This discovery provides insights into factors that may influence the evolution of sex chromosomes in plants and contributes to a large number of recent observations that underscore their dynamic nature.

The capability to generate densely sampled single nucleotide polymorphism (SNP) data is essential in diverse subdisciplines of biology, including crop breeding, pathology, forensics, forestry, ecology, evolution, and conservation. However, the wet-lab expertise and bioinformatics training required to conduct genome-scale variant discovery remain limiting factors for investigators with limited resources.Here we present ISSRseq, a PCR-based method for reduced representation of genomic variation using simple sequence repeats as priming sites to sequence inter simple sequence repeat (ISSR) regions. Briefly, ISSR regions are amplified with single primers, pooled, used to construct sequencing libraries with a commercially-available kit, and sequenced on the Illumina platform. We also present a flexible bioinformatic pipeline that assembles ISSR loci, calls and hard filters variants, outputs data matrices in common formats, and conducts population analyses using R.Using three angiosperm species as case studies, we demonstrate that ISSRseq is highly repeatable, necessitates only simple wet-lab skills and commonplace instrumentation, is flexible in terms of the number of single primers used, and can generate genomic-scale variant discovery on par with existing RRS methods which require more complex wet lab procedures.ISSRseq represents a straightforward approach to SNP genotyping in any organism, and we predict that this method will be particularly useful for those studying population genomics and phylogeography of non-model organisms. Furthermore, the ease of ISSRseq relative to other RRS methods should prove useful to those lacking advanced expertise in wet lab methods or bioinformatics.

Premise The family Salicaceae has proved taxonomically challenging, especially in the genus Salix, which is speciose and features frequent hybridization and polyploidy. Past efforts to reconstruct the phylogeny with molecular barcodes have failed to resolve the species relationships of many sections of the genus. Methods We used the wealth of sequence data in the family to design sequence capture probes to target regions of 300-1200 bp of exonic regions of 972 genes. Results We recovered sequence data for nearly all of the targeted genes in three species of Populus and three species of Salix. We present a species tree, discuss concordance among gene trees, and present population genomic summary statistics for these loci. Conclusions Our sequence capture array has extremely high capture efficiency within the genera Populus and Salix, resulting in abundant phylogenetic information. Additionally, these loci show promise for population genomic studies.

To understand the genetic mechanisms underlying wood anatomical and morphological traits in Populus trichocarpa, we used 869 unrelated genotypes from a common garden in Clatskanie, Oregon that were previously collected from across the distribution range in western North America. Using GEMMA mixed model analysis, we tested for the association of 25 phenotypic traits and nine multitrait combinations with 6.741 million SNPs covering the entire genome. Broad-sense trait heritabilities ranged from 0.117 to 0.477. Most traits were significantly correlated with geoclimatic variables suggesting a role of climate and geography in shaping the variation of this species. Fifty-seven SNPs from single trait GWAS and 11 SNPs from multitrait GWAS passed an FDR threshold of 0.05, leading to the identification of eight and seven nearby candidate genes, respectively. The percentage of phenotypic variance explained (PVE) by the significant SNPs for both single and multitrait GWAS ranged from 0.01% to 6.18%. To further evaluate the potential roles of candidate genes, we used a multi-omic network containing five additional data sets, including leaf and wood metabolite GWAS layers and coexpression and comethylation networks. We also performed a functional enrichment analysis on coexpression nearest neighbors for each gene model identified by the wood anatomical and morphological trait GWAS analyses. Genes affecting cell wall composition and transport related genes were enriched in wood anatomy and stomatal density trait networks. Signaling and metabolism related genes were also common in networks for stomatal density. For leaf morphology traits (leaf dry and wet weight) the networks were significantly enriched for GO terms related to photosynthetic processes as well as cellular homeostasis. The identified genes provide further insights into the genetic control of these traits, which are important determinants of the suitability and sustainability of improved genotypes for lignocellulosic biofuel production.

It is increasingly realized that homoploid hybrid speciation (HHS), which involves no change in chromosome number, is an important mechanism of speciation. HHS will likely increase in frequency as ecological and geographical barriers between species are continuing to be disrupted by human activities. HHS requires the establishment of reproductive isolation between a hybrid and its parents, but the underlying genes and genetic mechanisms remain largely unknown. In this study, we reveal by integrated approaches that reproductive isolation originates in one homoploid hybrid plant species through the inheritance of alternate alleles at genes that determine parental premating isolation. The parent species of this hybrid species are reproductively isolated by differences in flowering time and survivorship on soils containing high concentrations of iron. We found that the hybrid species inherits alleles of parental isolating major genes related to flowering time from one parent and alleles of major genes related to iron tolerance from the other parent. In this way, it became reproductively isolated from one parent by the difference in flowering time and from the other by habitat adaptation (iron tolerance). These findings and further modeling results suggest that HHS may occur relatively easily via the inheritance of alternate parental premating isolating genes and barriers.

The ages and sizes of a sex-determination region (SDR) are difficult to determine in non-model species. Due to the lack of recombination and enrichment of repetitive elements in SDRs, the quality of assembly with short sequencing reads is universally low. Unique features present in the SDRs help provide clues about how SDRs are established and how they evolve in the absence of recombination. Several Populus species have been reported with a male heterogametic configuration of sex (XX/XY system) mapped on chromosome 19, but the exact location of the SDR has been inconsistent among species, and thus far, none of these SDRs has been fully assembled in a genomic context. Here we identify the Y-SDR from a Y-linked contig directly from a long-read PacBio assembly of a Populus trichocarpa male individual. We also identified homologous gene sequences in the SDR of P. trichocarpa and the SDR of the W chromosome in Salix purpurea. We show that inverted repeats (IRs) found in the Y-SDR and the W-SDR are lineage-specific. We hypothesize that, although the two IRs are derived from the same orthologous gene within each species, they likely have independent evolutionary histories. Furthermore, the truncated inverted repeats in P. trichocarpa may code for small RNAs that target the homologous gene for RNA-directed DNA methylation. These findings support the hypothesis that diverse sex-determining systems may be achieved through similar evolutionary pathways, thereby providing a possible mechanism to explain the lability of sex-determination systems in plants in general.

Background Pyrolysis-molecular beam mass spectrometry (py-MBMS) analysis of a pedigree of Populus trichocarpa was performed to study the phenotypic plasticity and heritability of lignin content and lignin monomer composition. Instrumental and microspatial environmental variability were observed in the spectral features and corrected to reveal underlying genetic variance of biomass composition. Results Lignin-derived ions were highly impacted by microspatial environmental variation which demonstrates phenotypic plasticity of lignin composition in Populus trichocarpa biomass. Broad-sense heritability of lignin composition after correcting for microspatial and instrumental variation was determined to be H2 = 0.56 based on py-MBMS based on ions known to derive from lignin. Heritability of lignin monomeric syringyl/guaiacyl ratio (S/G) was H2 = 0.81. Broad-sense heritability was also high (up to H2 = 0.79) for ions derived from other components of the biomass including phenolics (e.g., salicylates) and C5 sugars (e.g., xylose). Lignin and phenolic ion abundances were primarily driven by maternal effects, and paternal effects were either similar or stronger for the most heritable carbohydrate-derived ions. Conclusions We have shown that many biopolymer-derived ions from py-MBMS show substantial phenotypic plasticity in response to microenvironmental variation in plantations. Nevertheless, broad-sense heritability for biomass composition can be quite high after correcting for spatial environmental variation. This work outlines the importance in accounting for instrumental and microspatial environmental variation in biomass composition data for applications in heritability measurements and genomic selection for breeding poplar for renewable fuels and materials.

Dioecy, the presence of separate sexes on distinct individuals, has evolved repeatedly in multiple plant lineages. However, the specific mechanisms through which sex systems evolve and their commonalities among plant species remain poorly understood. With both XY and ZW sex systems, the family Salicaceae provides a system to uncover the evolutionary forces driving sex chromosome turnovers. In this study, we performed a genome-wide association study to characterize sex determination in two Populus species, P. euphratica and P. alba. Our results reveal an XY system of sex determination on chromosome 14 of P. euphratica, and a ZW system on chromosome 19 of P. alba. We further assembled the corresponding sex determination regions, and found that their sex chromosome turnovers may be driven by the repeated translocations of a Helitron-like transposon. During the translocation, this factor may have captured partial or intact sequences that are orthologous to a type-A cytokinin response regulator gene. Based on results from this and other recently published studies, we hypothesize that this gene may act as a master regulator of sex determination for the entire family. We propose a general model to explain how the XY and ZW sex systems in this family can be determined by the same RR gene. Our study provides new insights into the diversification of incipient sex chromosome in flowering plants by showing how transposition and rearrangement of a single gene can control sex in both XY and ZW systems.

PremiseSexual dimorphism in dioecious plant species is often not obvious or is absent. Dioecious species populations also often exhibit deviations from expected sex ratios. Previous studies on members of the Salicaceae family have shown strong, partial, and no sexual dimorphism. Some studies have shown sex-biased ratios in severalSalixspp., however,S. purpureahas never been examined for evidence of sexual dimorphism or for the presence of sex-ratio bias, and therefore a comprehensive phenotypic study is needed to fill this knowledge gap.MethodsThis study examined a suite of morphological, phenological, physiological and wood composition traits from multi-environment and multi-year replicated field trials in a diversity panel of unrelatedS. purpureaaccessions and in full-sib F1and F2families produced through controlled cross pollinations to test for sexual dimorphism and sex ratio bias.Key ResultsSignificant evidence of sexual dimorphism was found in vegetative traits with greater means for many traits in male genotypes compared to females across three populations ofS. purpurea, measured across multiple years that were highly predictive of biomass yield. Male plants exhibited greater nitrogen accumulation under fertilizer amendment as measured by SPAD in the diversity panel, and males showed greater susceptibility to fungal infection byMelampsoraspp in the F2family. There were also consistent female-biased sex ratios in both the F1and F2families.ConclusionsThese results provide the first evidence of sexual dimorphism inS. purpureaand also confirm the prevalence of female-biased sex ratios previously found in otherSalixspecies.

The development of non-recombining sex chromosomes has radical effects on the evolution of discrete sexes and sexual dimorphism. Although dioecy is rare in plants, sex chromosomes have evolved repeatedly throughout the diversification of angiosperms, and many of these sex chromosomes are relatively young compared to those found in vertebrates. In this study, we designed and used a sequence capture array to identify a novel sex-linked region (SLR) inSalix nigra, a basal species in the willow clade, and demonstrated that this species has XY heterogamety. We did not detect any genetic overlap with the previously characterized ZW SLRs in willows, which map to a different chromosome. TheS. nigraSLR is characterized by strong recombination suppression across a 2 MB region and an excess of low frequency alleles, resulting in a low Tajima’s D compared to the remainder of the genome. We speculate that either a recent bottleneck in population size or factors related to positive or background selection generated this differential pattern of Tajima’s D on the X and autosomes. This discovery provides insights into factors that may influence the evolution of sex chromosomes in plants and contributes to a large number of recent observations that underscore their dynamic nature.

BackgroundSex chromosomes in a wide variety of species share common characteristics, including the presence of suppressed recombination surrounding sex determination loci. They have arisen independently in numerous lineages, providing a conclusive example of convergent evolution. Mammalian sex chromosomes contain multiple palindromic repeats across the non-recombining region that facilitate sequence conservation through gene conversion, and contain genes that are crucial for sexual reproduction. Plant sex chromosomes are less well understood, and in particular it is not clear how coding sequence conservation is maintained in the absence of homologous recombination.ResultsHere we present the first evidence of large palindromic structures in a plant sex chromosome, based on a highly contiguous assembly of the W chromosome of the dioecious shrubSalix purpurea. Two consecutive palindromes span over a region of 200 kb, with conspicuous 20 kb stretches of highly conserved sequences among the four arms. The closely-related speciesS. suchowensisalso has two copies of a portion of the palindrome arm and provides strong evidence for gene conversion. Four genes in the palindrome are homologous to genes in the SDR of the closely-related genusPopulus, which is located on a different chromosome. These genes show distinct, floral-biased expression patterns compared to paralogous copies on autosomes.ConclusionThe presence of palindromic structures in sex chromosomes of mammals and plants highlights the intrinsic importance of these features in adaptive evolution in the absence of recombination. Convergent evolution is driving both the independent establishment of sex chromosomes as well as their fine-scale sequence structure.

In hybrid zones occurring in marginal environments, adaptive introgression from one species into the genomic background of another may constitute a mechanism facilitating adaptation at range limits. Although recent studies have improved our understanding of adaptive introgression in widely distributed tree species, little is known about the dynamics of this process in populations at the margins of species ranges. We investigated the extent of introgression between three species of the genus Populus sect. Tacamahaca (P. balsamifera, P. angustifolia and P. trichocarpa) at the margins of their distributions in the Rocky Mountain region of the United States and Canada. Using genotyping by sequencing (GBS), we analysed ~ 83,000 single nucleotide polymorphisms genotyped in 296 individuals from 29 allopatric and sympatric populations of the three species. We found a trispecies hybrid complex present throughout the zone of range overlap, including early as well as advanced generation backcross hybrids, indicating recurrent gene flow in this hybrid complex. Using genomic cline analysis, we found evidence of non-neutral patterns of introgression at 23% of loci in hybrids, of which 47% and 8% represented excess ancestry from P. angustifolia and P. balsamifera, respectively. Gene ontology analysis suggested these genomic regions were enriched for genes associated with photoperiodic regulation, metal ion transport, maintenance of redox homeostasis and cell wall metabolites involved in regulation of seasonal dormancy. Our study demonstrates the role of adaptive introgression in a multispecies hybrid complex in range-edge populations and has implications for understanding the evolutionary dynamics of adaptation in hybrid zones, especially at the margins of species distributions.

Long-lived perennial plants, with distinctive habits of inter-annual growth, defense, and physiology, are of great economic and ecological importance. However, some biological mechanisms resulting from genome duplication and functional divergence of genes in these systems remain poorly studied. Here, we discovered an association between a poplar (Populus trichocarpa) 5-enolpyruvylshikimate 3-phosphate synthase gene (PtrEPSP) and lignin biosynthesis. Functional characterization of PtrEPSP revealed that this isoform possesses a helix-turn-helix motif in the N terminus and can function as a transcriptional repressor that regulates expression of genes in the phenylpropanoid pathway in addition to performing its canonical biosynthesis function in the shikimate pathway. We demonstrated that this isoform can localize in the nucleus and specifically binds to the promoter and represses the expression of a SLEEPER-like transcriptional regulator, which itself specifically binds to the promoter and represses the expression of PtrMYB021 (known as MYB46 in Arabidopsis thaliana), a master regulator of the phenylpropanoid pathway and lignin biosynthesis. Analyses of overexpression and RNAi lines targeting PtrEPSP confirmed the predicted changes in PtrMYB021 expression patterns. These results demonstrate that PtrEPSP in its regulatory form and PtrhAT form a transcriptional hierarchy regulating phenylpropanoid pathway and lignin biosynthesis in Populus.

Genomic studies have been used to identify genes underlying many important plant secondary metabolic pathways. However, genes for salicinoid phenolic glycosides (SPGs)—ecologically important compounds with significant commercial, cultural, and medicinal applications—remain largely undescribed. We used a linkage map derived from a full‐sib population of hybrid cottonwoods (Populus spp.) to search for quantitative trait loci (QTL) for the SPGs salicortin and HCH‐salicortin. SSR markers and primer sequences were used to anchor the map to the V3.0 P. trichocarpa genome. We discovered 21 QTL for the two traits, including a major QTL for HCH‐salicortin (R 2 = .52) that colocated with a QTL for salicortin on chr12. Using the V3.0 Populus genome sequence, we identified 2,983 annotated genes and 1,480 genes of unknown function within our QTL intervals. We note ten candidate genes of interest, including a BAHD‐type acyltransferase that has been potentially linked to Populus SPGs. Our results complement other recent studies in Populus with implications for gene discovery and the evolution of defensive chemistry in a model genus. To our knowledge, this is the first study to use a full‐sib mapping population to identify QTL intervals and gene lists associated with SPGs.

Populus has been investigated as a promising biomass feedstock for alternative fuels and chemicals. Physicochemical characteristics and genomic information of biomass feedstocks are among the essential information that can help not only advance our understanding of biomass recalcitrance but also in its efficient utilization. Herein, the recalcitrance of natural Populus variants was elucidated in three aspects: (1) sugar release, (2) physicochemical properties, and (3) relative variation of gene expression within natural poplar variants. The sugar release performance of natural Populus variants was evaluated with their correlation with biomass recalcitrance-related characteristics. Among the physicochemical properties of poplar, the lignin content, lignin molecular weight, lignin S/G ratio, and cellulose accessibility were found to correlate with sugar release. The results demonstrated that the lignin content was negatively correlated with sugar release, whereas the lignin molecular weight, lignin S/G ratio, and cellulose accessibility were positively associated with poplar sugar release. The trend of differential gene expression of each variant also supports the characterization results and their effects on biomass conversion.

The presence of thousands of independent origins of dioecy in angiosperms provides a unique opportunity to address the parallel evolution of the molecular pathways underlying unisexual flowers. Recent progress towards identifying sex determination genes has identified hormone response pathways, mainly associated with cytokinin and ethylene response pathways, as having been recruited multiple times independently to control unisexuality. Moreover, transcriptomics has begun to identify commonalities among intermediate sections of signal transduction pathways. These recent advances set the stage for development of a comparative evolutionary development research program to identify the shared and unique aspects of the genetic pathways of unisexual flower development in angiosperms.

Salix purpurea is a woody perennial that is bred as a high-yielding bioenergy crop in North America. To gain an understanding of the genotypic variation associated with phenotypic diversity, this study characterized the population structure and genetic diversity of S. purpurea from its native range of Europe and naturalized range of the Northeastern United States (US). A total of 273 genotypes of S. purpurea were analyzed, which included 95 naturalized accessions and 19 horticultural cultivars from the US and 159 accessions collected from the native range of four European countries. All individuals were evaluated using a filtered set of 2287 genotyping-by-sequencing (GBS) single nucleotide polymorphism (SNP) markers. Using five clustering techniques (PCA, neighbor joining, STRUCTURE, DAPC, and affinity propagation), population structure was resolved into three broadly classified groups. Further analysis revealed seven to eight subpopulations which corresponded to geographical collection sites, where performance of the DAPC and affinity propagation methods was superior to STRUCTURE analysis for the purposes of characterizing population structure and performing population assignment. The native European accessions exhibited greater diversity and subpopulation structure than the US naturalized accessions, where there was a clear geographical delineation between the alpine/subalpine collections and the lowland collections at the Baltic Sea and Oder River. We also show that a subset of the horticultural cultivars had a higher likelihood of similarity to US naturalized populations which display hybrid ancestry, where both naturalized and cultivated genotypes appear to share a common ancestor. Additionally, several accessions collected from different sites were found to be clonal. Ongoing and future conservation and association studies will benefit from these known substructures and diversity assessments.

Genome-wide association studies (GWAS) have great promise for identifying the loci that contribute to adaptive variation, but the complex genetic architecture of many quantitative traits presents a substantial challenge. We measured 14 morphological and physiological traits and identified single nucleotide polymorphism (SNP)-phenotype associations in a Populus trichocarpa population distributed from California, USA to British Columbia, Canada. We used whole-genome resequencing data of 882 trees with more than 6.78 million SNPs, coupled with multitrait association to detect polymorphisms with potentially pleiotropic effects. Candidate genes were validated with functional data. Broad-sense heritability (H2 ) ranged from 0.30 to 0.56 for morphological traits and 0.08 to 0.36 for physiological traits. In total, 4 and 20 gene models were detected using the single-trait and multitrait association methods, respectively. Several of these associations were corroborated by additional lines of evidence, including co-expression networks, metabolite analyses, and direct confirmation of gene function through RNAi. Multitrait association identified many more significant associations than single-trait association, potentially revealing pleiotropic effects of individual genes. This approach can be particularly useful for challenging physiological traits such as water-use efficiency or complex traits such as leaf morphology, for which we were able to identify credible candidate genes by combining multitrait association with gene co-expression and co-methylation data.

Dioecy has evolved numerous times in plants, but heteromorphic sex chromosomes are apparently rare. Sex determination has been studied in multiple Salix and Populus (Salicaceae) species, and P. trichocarpa has an XY sex determination system on chromosome 19, while S. suchowensis and S. viminalis have a ZW system on chromosome 15. Here we use whole genome sequencing coupled with quantitative trait locus mapping and a genome-wide association study to characterize the genomic composition of the non-recombining portion of the sex determination region. We demonstrate that Salix purpurea also has a ZW system on chromosome 15. The sex determination region has reduced recombination, high structural polymorphism, an abundance of transposable elements, and contains genes that are involved in sex expression in other plants. We also show that chromosome 19 contains sex-associated markers in this S. purpurea assembly, along with other autosomes. This raises the intriguing possibility of a translocation of the sex determination region within the Salicaceae lineage, suggesting a common evolutionary origin of the Populus and Salix sex determination loci.

Biological organisms are complex systems that are composed of functional networks of interacting molecules and macromolecules. Complex phenotypes are the result of orchestrated, hierarchical, heterogeneous collections of expressed genomic variants. However, the effects of these variants are the result of historic selective pressure and current environmental and epigenetic signals, and, as such, their co-occurrence can be seen as genome-wide correlations in a number of different manners. Biomass recalcitrance (i.e., the resistance of plants to degradation or deconstruction, which ultimately enables access to a plant’s sugars) is a complex polygenic phenotype of high importance to biofuels initiatives. This study makes use of data derived from the re-sequenced genomes from over 800 different Populus trichocarpa genotypes in combination with metabolomic and pyMBMS data across this population, as well as co-expression and co-methylation networks in order to better understand the molecular interactions involved in recalcitrance, and identify target genes involved in lignin biosynthesis/degradation. A Lines Of Evidence (LOE) scoring system is developed to integrate the information in the different layers and quantify the number of lines of evidence linking genes to lignin-related lignin-phenotypes across the network layers. The resulting Genome Wide Association Study networks, integrated with Single Nucleotide Polymorphism (SNP) correlation, co-methylation and co-expression networks through the LOE scores are proving to be a powerful approach to determine the pleiotropic and epistatic relationships underlying cellular functions and, as such, the molecular basis for complex phenotypes, such as recalcitrance.

Background One of the major barriers to the development of lignocellulosic feedstocks is the recalcitrance of plant cell walls to deconstruction and saccharification. Recalcitrance can be reduced by targeting genes involved in cell wall biosynthesis, but this can have unintended consequences that compromise the agronomic performance of the trees under field conditions. Here we report the results of a field trial of fourteen distinct transgenic Populus deltoides lines that had previously demonstrated reduced recalcitrance without yield penalties under greenhouse conditions. Results Survival and productivity of the trial were excellent in the first year, and there was little evidence for reduced performance of the transgenic lines with modified target gene expression. Surprisingly, the most striking phenotypic effects in this trial were for two empty-vector control lines that had modified bud set and bud flush. This is most likely due to somaclonal variation or insertional mutagenesis. Traits related to yield, crown architecture, herbivory, pathogen response, and frost damage showed few significant differences between target gene transgenics and empty vector controls. However, there were a few interesting exceptions. Lines overexpressing the DUF231 gene, a putative O-acetyltransferase, showed early bud flush and marginally increased height growth. Lines overexpressing the DUF266 gene, a putative glycosyltransferase, had significantly decreased stem internode length and slightly higher volume index. Finally, lines overexpressing the PFD2 gene, a putative member of the prefoldin complex, had a slightly reduced volume index. Conclusions This field trial demonstrates that these cell wall modifications, which decreased cell wall recalcitrance under laboratory conditions, did not seriously compromise first-year performance in the field, despite substantial challenges, including an outbreak of a stem boring insect (Gypsonoma haimbachiana), attack by a leaf rust pathogen (Melampsora spp.), and a late frost event. This bodes well for the potential utility of these lines as advanced biofuels feedstocks. Electronic supplementary material The online version of this article (doi:10.1186/s13068-017-0934-6) contains supplementary material, which is available to authorized users.

Establishing links between phenotypes and molecular variants is of central importance to accelerate genetic improvement of economically important plant species. Our work represents the first genome-wide association study to the inherently complex and currently poorly understood genetic architecture of industrially relevant wood traits. Here, we employed an Illumina Infinium 34K single nucleotide polymorphism (SNP) genotyping array that generated 29,233 high-quality SNPs in c. 3500 broad-based candidate genes within a population of 334 unrelated Populus trichocarpa individuals to establish genome-wide associations. The analysis revealed 141 significant SNPs (α ≤ 0.05) associated with 16 wood chemistry/ultrastructure traits, individually explaining 3-7% of the phenotypic variance. A large set of associations (41% of all hits) occurred in candidate genes preselected for their suggested a priori involvement with secondary growth. For example, an allelic variant in the FRA8 ortholog explained 21% of the total genetic variance in fiber length, when the trait's heritability estimate was considered. The remaining associations identified SNPs in genes not previously implicated in wood or secondary wall formation. Our findings provide unique insights into wood trait architecture and support efforts for population improvement based on desirable allelic variants.

Knowledge of genetic variation in response of Populus to heavy metals like cadmium (Cd) is an important step in understanding the underlying mechanisms of tolerance. In this study, a pseudo-backcross pedigree of Populus trichocarpa and Populus deltoides was characterized for Cd exposure. The pedigree showed significant variation for Cd tolerance thus enabling the identification of relatively tolerant and susceptible genotypes for intensive characterization. A total of 16 QTLs at logarithm of odds (LOD) ratio > 2.5, were found to be associated with total dry weight, its components, and root volume. Four major QTLs for total dry weight were mapped to different linkage groups in control (LG III) and Cd conditions (LG XVI) and had opposite allelic effects on Cd tolerance, suggesting that these genomic regions were differentially controlled. The phenotypic variation explained by Cd QTL for all traits under study varied from 5.9% to 11.6% and averaged 8.2% across all QTL. Leaf Cd contents also showed significant variation suggesting the phytoextraction potential of Populus genotypes, though heritability of this trait was low (0.22). A whole-genome microarray study was conducted by using two genotypes with extreme responses for Cd tolerance in the above study and differentially expressed genes were identified. Candidatemore » genes including CAD2 (CADMIUM SENSITIVE 2), HMA5 (HEAVY METAL ATPase5), ATGTST1 (Arabidopsis thaliana Glutathione S-Transferase1), ATGPX6 (Glutathione peroxidase 6), and ATMRP 14 (Arabidopsis thaliana Multidrug Resistance associated Protein 14) were identified from QTL intervals and microarray study. Functional characterization of these candidate genes could enhance phytoremediation capabilities of Populus.« less

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

Recent advances in Populus genomics have greatly expanded its popularity as a reference for fundamental as well as applied research in woody plants. In this review, we provide an overview of the state-of-the-science in Populus genomics research, including experimental and computational genomics. We have surveyed and summarized the following: 1) pioneering as well as more recent reports of genetics- and genomics-based investigations in Populus, 2) the positive impact of technological improvements, 3) findings from phylogenetic analyses of gene families, and 4) genomic databases. In the area of Populus experimental genomics, genetic approaches have been advanced to the genome scale with resolution to the gene and/or single nucleotide level. On the other hand, the modern “omics” approaches have been successfully applied to analysis of gene function, such as transcriptome profiling using microarrays as well as the next-generation DNA sequencing technology, and characterization of proteome and metabolome using...

The 18S-28S rDNA and 5S rDNA loci in Populus trichocarpa were localized using fluorescent in situ hybridization (FISH). Two 18S-28S rDNA sites and one 5S rDNA site were identified and located at the ends of 3 different chromosomes. FISH signals from the Arabidopsis-type telomere repeat sequence were observed at the distal ends of each chromosome. Six BAC clones selected from 2 linkage groups based on genome sequence assembly (LG-I and LG-VI) were localized on 2 chromosomes, as expected. BACs from LG-I hybridized to the longest chromosome in the complement. All BAC positions were found to be concordant with sequence assembly positions. BAC-FISH will be useful for delineating each of the Populus trichocarpa chromosomes and improving the sequence assembly of this model angiosperm tree species.

We investigated the efficiency of RNA interference (RNAi) in Arabidopsis using transitive and homologous inverted repeat (hIR) vectors. hIR constructs carry self-complementary intron-spliced fragments of the target gene whereas transitive vectors have the target sequence fragment adjacent to an intron-spliced, inverted repeat of heterologous origin. Both transitive and hIR constructs facilitated specific and heritable silencing in the three genes studied (AP1, ETTIN and TTG1). Both types of vectors produced a phenotypic series that phenocopied reduction of function mutants for the respective target gene. The hIR yielded up to fourfold higher proportions of events with strongly manifested reduction of function phenotypes compared to transitive RNAi. We further investigated the efficiency and potential off-target effects of AP1 silencing by both types of vectors using genome-scale microarrays and quantitative RT-PCR. The depletion of AP1 transcripts coincided with reduction of function phenotypic changes among both hIR and transitive lines and also showed similar expression patterns among differentially regulated genes. We did not detect significant silencing directed against homologous potential off-target genes when constructs were designed with minimal sequence similarity. Both hIR and transitive methods are useful tools in plant biotechnology and genomics. The choice of vector will depend on specific objectives such as cloningmore » throughput, number of events and degree of suppression required.« less