Evidence for widespread selection in shaping the genomic landscape during speciation of Populus

Increasing our understanding of how evolutionary processes drive the genomic landscape of variation is fundamental to a better understanding of the genomic consequences of speciation. However, genome‐wide patterns of within‐ and between‐ species variation have not been fully investigated in most forest tree species despite their global ecological and economic importance. Here, we use whole‐genome resequencing data from four Populus species spanning the speciation continuum to reconstruct their demographic histories and investigate patterns of diversity and divergence within and between species. Using Populus trichocarpa as an outgroup species, we further infer the genealogical relationships and estimate the extent of ancient introgression among the three aspen species (Populus tremula, Populus davidiana and Populus tremuloides) throughout the genome. Our results show substantial variation in these patterns along the genomes with this variation being strongly predicted by local recombination rates and the density of functional elements. This implies that the interaction between recurrent selection and intrinsic genomic features has dramatically sculpted the genomic landscape over long periods of time. In addition, our findings provide evidence that, apart from background selection, recent positive selection and long‐term balancing selection have also been crucial components in shaping patterns of genome‐wide variation during the speciation process.     

Scale and direction of adaptive introgression between black cottonwood (Populus trichocarpa) and balsam poplar (P. balsamifera)

Introgression can introduce novel genetic variation at a faster rate than mutation alone and result in adaptive introgression when adaptive alleles are maintained in the recipient genome over time by natural selection. A previous study from our group demonstrated adaptive introgression from Populus balsamifera into P. trichocarpa in a target genomic region. Here we expand our local ancestry analysis to the whole genome of both parents to provide a comprehensive view of introgression patterns and to identify additional candidate regions for adaptive introgression genomewide. Populus trichocarpa is a large, fast‐growing tree of mild coastal regions of the Pacific Northwest, whereas P. balsamifera is a smaller stature tree of continental and boreal regions with intense winter cold. The species hybridize where they are parapatric. We detected asymmetric patterns of introgression across the whole genome of these two poplar species adapted to contrasting environments, with stronger introgression from P. balsamifera to P. trichocarpa than vice versa. Admixed P. trichocarpa individuals contained more genomic regions with unusually high levels of introgression (19 regions) and also the largest introgressed genome fragment (1.02 Mb) compared with admixed P. balsamifera (nine regions). Our analysis also revealed numerous candidate regions for adaptive introgression with strong signals of selection, notably related to disease resistance, and enriched for genes that may play crucial roles in survival and adaptation. Furthermore, we detected a potential overrepresentation of subtelomeric regions in P. balsamifera introgressed into P. trichocarpa and possible protection of sex‐determining regions from interspecific gene flow.     

Genome sequence and genetic transformation of a widely distributed and cultivated poplar

Populus alba is widely distributed and cultivated in Europe and Asia. This species has been used for diverse studies. In this study, we assembled a de novo genome sequence of P. alba var. pyramidalis (= P. bolleana) and confirmed its high transformation efficiency and short transformation time by experiments. Through a process of hybrid genome assembly, a total of 464 M of the genome was assembled. Annotation analyses predicted 37 901 protein‐coding genes. This genome is highly collinear to that of P. trichocarpa, with most genes having orthologs in the two species. We found a marked expansion of gene families related to histone and the hormone auxin but loss of disease resistance genes in P. alba if compared with the closely related P. trichocarpa. The genome sequence presented here represents a valuable resource for further molecular functional analyses of this species as a new tree model, poplar breeding practices and comparative genomic analyses across different poplars.     

Multilocus analysis of nucleotide variation and speciation in three closely related Populus (Salicaceae) species

Historical tectonism and climate oscillations can isolate and contract the geographical distributions of many plant species, and they are even known to trigger species divergence and ultimately speciation. Here, we estimated the nucleotide variation and speciation in three closely related Populus species, Populus tremuloides, P. tremula and P. davidiana, distributed in North America and Eurasia. We analysed the sequence variation in six single‐copy nuclear loci and three chloroplast (cpDNA) fragments in 497 individuals sampled from 33 populations of these three species across their geographic distributions. These three Populus species harboured relatively high levels of nucleotide diversity and showed high levels of nucleotide differentiation. Phylogenetic analysis revealed that P. tremuloides diverged earlier than the other two species. The cpDNA haplotype network result clearly illustrated the dispersal route from North America to eastern Asia and then into Europe. Molecular dating results confirmed that the divergence of these three species coincided with the sundering of the Bering land bridge in the late Miocene and a rapid uplift of the Qinghai‐Tibetan Plateau around the Miocene/Pliocene boundary. Vicariance‐driven successful allopatric speciation resulting from historical tectonism and climate oscillations most likely played roles in the formation of the disjunct distributions and divergence of these three Populus species.     

Herbivore‐induced poplar cytochrome P450 enzymes of the CYP71 family convert aldoximes to nitriles which repel a generalist caterpillar

Numerous plant species emit volatile nitriles upon herbivory, but the biosynthesis as well as the relevance of these nitrogenous compounds in plant–insect interactions remains unknown. Populus trichocarpa has been shown to produce a complex blend of nitrogenous volatiles, including aldoximes and nitriles, after herbivore attack. The aldoximes were previously reported to be derived from amino acids by the action of cytochrome P450 enzymes of the CYP79 family. Here we show that nitriles are derived from aldoximes by another type of P450 enzyme in P. trichocarpa. First, feeding of deuterium‐labeled phenylacetaldoxime to poplar leaves resulted in incorporation of the label into benzyl cyanide, demonstrating that poplar volatile nitriles are derived from aldoximes. Then two P450 enzymes, CYP71B40v3 and CYP71B41v2, were characterized that produce aliphatic and aromatic nitriles from their respective aldoxime precursors. Both possess typical P450 sequence motifs but do not require added NADPH or cytochrome P450 reductase for catalysis. Since both enzymes are expressed after feeding by gypsy moth caterpillars, they are likely to be involved in herbivore‐induced volatile nitrile emission in P. trichocarpa. Olfactometer experiments showed that these volatile nitriles have a strong repellent activity against gypsy moth caterpillars, suggesting they play a role in induced direct defense against poplar herbivores.     

Genotyping‐by‐sequencing and ecological niche modeling illuminate phylogeography,admixture, and Pleistocene range dynamics in quaking aspen (Populus tremuloides)

Populus tremuloides is the widest‐ranging tree species in North America and an ecologically important component of mesic forest ecosystems displaced by the Pleistocene glaciations. Using phylogeographic analyses of genome‐wide SNPs (34,796 SNPs, 183 individuals) and ecological niche modeling, we inferred population structure, ploidy levels, admixture, and Pleistocene range dynamics of P. tremuloides, and tested several historical biogeographical hypotheses. We found three genetic lineages located mainly in coastal–Cascades (cluster 1), east‐slope Cascades–Sierra Nevadas–Northern Rockies (cluster 2), and U.S. Rocky Mountains through southern Canadian (cluster 3) regions of the P. tremuloides range, with tree graph relationships of the form ((cluster 1, cluster 2), cluster 3). Populations consisted mainly of diploids (86%) but also small numbers of triploids (12%) and tetraploids (1%), and ploidy did not adversely affect our genetic inferences. The main vector of admixture was from cluster 3 into cluster 2, with the admixture zone trending northwest through the Rocky Mountains along a recognized phenotypic cline (Utah to Idaho). Clusters 1 and 2 provided strong support for the “stable‐edge hypothesis” that unglaciated southwestern populations persisted in situ since the last glaciation. By contrast, despite a lack of clinal genetic variation, cluster 3 exhibited “trailing‐edge” dynamics from niche suitability predictions signifying complete northward postglacial expansion. Results were also consistent with the “inland dispersal hypothesis” predicting postglacial assembly of Pacific Northwestern forest ecosystems, but rejected the hypothesis that Pacific‐coastal populations were colonized during outburst flooding from glacial Lake Missoula. Overall, congruent patterns between our phylogeographic and ecological niche modeling results and fossil pollen data demonstrate complex mixtures of stable‐edge, refugial locations, and postglacial expansion within P. tremuloides. These findings confirm and refine previous genetic studies, while strongly supporting a distinct Pacific‐coastal genetic lineage of quaking aspen.     

Improved genome assembly provides new insights into genome evolution in a desert poplar (Populus euphratica)

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.     

Genetic linkage mapping in aspen (Populus tremula L. and Populus tremuloides Michx.)

A large number of simple sequence repeat (SSR) marker-containing genetic maps are available for several Populus species. For aspen however, no SSR-containing map has been published so far. In this study, genetic linkage mapping was carried out with an interspecific mapping pedigree of 61 full-sib hybrids of European × quaking aspen (Populus tremula L. × Populus tremuloides Michx.), using the two-way pseudo-testcross strategy. Amplified fragment-length polymorphism (AFLP) and SSR markers were used for mapping, resulting in the first SSR-containing genetic linkage maps for aspen. The maps allow comparisons with a Populus consensus map and other published genetic maps of the genus Populus. The maps showed good collinearity to each other and to the Populus consensus map and provide a direct link to the Populus trichocarpa genomic sequence. Sex as a morphological trait was assessed in the mapping population and mapped on a non-terminal position of linkage group XIX on the male P. tremuloides map.     

Secondary successional trajectories of structural and catabolic bacterial communities in oil‐polluted soil planted with hybrid poplar

Poplars have widely been used for rhizoremediation of a broad range of organic contaminants for the past two decades. Still, there is a knowledge gap regarding the rhizosphere‐associated bacterial communities of poplars and their dynamics during the remediation process. It is envisaged that a detailed understanding of rhizosphere‐associated microbial populations will greatly contribute to a better design and implementation of rhizoremediation. To investigate the long‐term succession of structural and catabolic bacterial communities in oil‐polluted soil planted with hybrid poplar, we carried out a 2‐year field study. Hybrid aspen (Populus tremula × Populus tremuloides) seedlings were planted in polluted soil excavated from an accidental oil‐spill site. Vegetated and un‐vegetated soil samples were collected for microbial community analyses at seven different time points during the course of 2 years and sampling time points were chosen to cover the seasonal variation in the boreal climate zone. Bacterial community structure was accessed by means of 16S rRNA gene amplicon pyrosequencing, whereas catabolic diversity was monitored by pyrosequencing of alkane hydroxylase and extradiol dioxygenase genes. We observed a clear succession of bacterial communities on both structural and functional levels from early to late‐phase communities. Sphingomonas type extradiol dioxygenases and alkane hydroxylase homologs of Rhodococcus clearly dominated the early‐phase communities. The high‐dominance/low‐diversity functional gene communities underwent a transition to low‐dominance/high‐diversity communities in the late phase. These results pointed towards increased catabolic capacities and a change from specialist to generalist strategy of bacterial communities during the course of secondary succession.     

Interspecific crossability studies in poplars

Summary Pollen-stigma interactions were studied with scanning electron microscopy in intraspecific and intersectional crosses of Populus tremuloides (Sect. LEUCE); P. deltoides, P. nigra (Sect. AIGEIROS); and P. trichocarpa (Sect. TACAMAHACA). Intraspecific variation in hydration, size of hydrated pollen grains, and tube morphology is described. Exine sculpture in P. tremuloides was densely spinulate or scabrate, in P. trichocarpa densely reticulate, in the other species somewhat intermediate. In six of seven intersectional cross combinations pollen-tube behavior on the stigmatic surface at 5–22 hours after pollination was the same as in conspecific matings; in P. tremuloides x P. trichocarpa crosses, the percentage of twisted and bulging tubes was significantly higher. Irradiated mentor pollen did not affect pollentube behavior of adjacent incompatible pollen, but increased hybrid seedling production in interspecific crosses. By contrast, hexane treatment of pollen or stigma, and the spraying of catkins with gibberellic acid or naphtalene acetic acid were ineffective. The results are discussed in relation to several mechanisms presumed to be involved in restricting gene exchange among poplar species; and in relation to Hogenboom's (1973) incongruity concept.     

Isolation and characterization of microsatellites in trembling aspen (Populus tremuloides)

 We have identified, isolated, and characterized microsatellite/simple sequence repeat (SSR) loci in trembling aspen (Populus tremuloides) by screening partial genomic libraries. We have also examined the compatibility and use of the P. tremuloides SSR primers to resolve microsatellites in other Populus species. Fourteen microsatellites were identified from 1600 clones screened. The TC/AG microsatellites were the most abundant. A total of 29 alleles were detected in 36 P. tremuloides individuals at the four SSR loci (two each of di- and tri-nucleotide repeats) characterized. The number of alleles at the SSR loci ranged from 5 to 11, with an average of 7.25 alleles per locus, and the observed heterozygosity ranged from 0.19 to 0.82, with a mean of 0.46 per locus. Although the highest polymorphism was observed for a dinucleotide SSR locus, the trinucleotide SSR loci showed substantial polymorphism. There were 34 unique multilocus genotypes among the 36 P. tremuloides individuals examined, and 89% of the individuals had unique multilocus genotypes. Two pairs of SSR primers were successful in PCR, amplifying genomic DNA and resolving microsatellites of comparable size from Populus deltoides, P. nigra, P.×canadensis, and P. maximowiczii. The microsatellite DNA markers developed could be used for clonal fingerprinting, certification of controlled crosses, genome mapping, marker-assisted early selection, genetic diversity assessments, and conservation and sustainable management of poplar genetic resources. Received: 14 November 1997 / Accepted: 17 November 1997     

New resources for genetic studies in Populus nigra: genome‐wide SNP discovery and development of a 12k Infinium array

Whole genome resequencing of 51 Populus nigra (L.) individuals from across Western Europe was performed using Illumina platforms. A total number of 1 878 727 SNPs distributed along the P. nigra reference sequence were identified. The SNP calling accuracy was validated with Sanger sequencing. SNPs were selected within 14 previously identified QTL regions, 2916 expressional candidate genes related to rust resistance, wood properties, water‐use efficiency and bud phenology and 1732 genes randomly spread across the genome. Over 10 000 SNPs were selected for the construction of a 12k Infinium Bead‐Chip array dedicated to association mapping. The SNP genotyping assay was performed with 888 P. nigra individuals. The genotyping success rate was 91%. Our high success rate was due to the discovery panel design and the stringent parameters applied for SNP calling and selection. In the same set of P. nigra genotypes, linkage disequilibrium throughout the genome decayed on average within 5–7 kb to half of its maximum value. As an application test, ADMIXTURE analysis was performed with a selection of 600 SNPs spread throughout the genome and 706 individuals collected along 12 river basins. The admixture pattern was consistent with genetic diversity revealed by neutral markers and the geographical distribution of the populations. These newly developed SNP resources and genotyping array provide a valuable tool for population genetic studies and identification of QTLs through natural‐population based genetic association studies in P. nigra.     

Thaumatin-like proteins are differentially expressed and localized in phloem tissues of hybrid poplar

Background  

Two thaumatin-like proteins (TLPs) were previously identified in phloem exudate of hybrid poplar (Populus trichocarpa × P. deltoides) using proteomics methods, and their sieve element localization confirmed by immunofluorescence. In the current study, we analyzed different tissues to further understand TLP expression and localization in poplar, and used immunogold labelling to determine intracellular localization.     

Population genomics of the eastern cottonwood (Populus deltoides)

Despite its economic importance as a bioenergy crop and key role in riparian ecosystems, little is known about genetic diversity and adaptation of the eastern cottonwood (Populus deltoides). Here, we report the first population genomics study for this species, conducted on a sample of 425 unrelated individuals collected in 13 states of the southeastern United States. The trees were genotyped by targeted resequencing of 18,153 genes and 23,835 intergenic regions, followed by the identification of single nucleotide polymorphisms (SNPs). This natural P. deltoides population showed low levels of subpopulation differentiation (F_ST = 0.022–0.106), high genetic diversity (θ_W = 0.00100, π = 0.00170), a large effective population size (N_e ≈ 32,900), and low to moderate levels of linkage disequilibrium. Additionally, genomewide scans for selection (Tajima's D), subpopulation differentiation (X^TX), and environmental association analyses with eleven climate variables carried out with two different methods (LFMM and BAYENV2) identified genes putatively involved in local adaptation. Interestingly, many of these genes were also identified as adaptation candidates in another poplar species, Populus trichocarpa, indicating possible convergent evolution. This study constitutes the first assessment of genetic diversity and local adaptation in P. deltoides throughout the southern part of its range, information we expect to be of use to guide management and breeding strategies for this species in future, especially in the face of climate change.     

Factors affecting fall down rates of dead aspen (Populus tremuloides) biomass following severe drought in west‐central Canada

Increases in mortality of trembling aspen (Populus tremuloides Michx.) have been recorded across large areas of western North America following recent periods of exceptionally severe drought. The resultant increase in standing, dead tree biomass represents a significant potential source of carbon emissions to the atmosphere, but the timing of emissions is partially driven by dead‐wood dynamics which include the fall down and breakage of dead aspen stems. The rate at which dead trees fall to the ground also strongly influences the period over which forest dieback episodes can be detected by aerial surveys or satellite remote sensing observations. Over a 12‐year period (2000–2012), we monitored the annual status of 1010 aspen trees that died during and following a severe regional drought within 25 study areas across west‐central Canada. Observations of stem fall down and breakage (snapping) were used to estimate woody biomass transfer from standing to downed dead wood as a function of years since tree death. For the region as a whole, we estimated that >80% of standing dead aspen biomass had fallen after 10 years. Overall, the rate of fall down was minimal during the year following stem death, but thereafter fall rates followed a negative exponential equation with k = 0.20 per year. However, there was high between‐site variation in the rate of fall down (k = 0.08–0.37 per year). The analysis showed that fall down rates were positively correlated with stand age, site windiness, and the incidence of decay fungi (Phellinus tremulae (Bond.) Bond. and Boris.) and wood‐boring insects. These factors are thus likely to influence the rate of carbon emissions from dead trees following periods of climate‐related forest die‐off episodes.     

Occupancy of red‐naped sapsuckers in a coniferous forest: using LiDAR to understand effects of vegetation structure and disturbance

Red‐naped sapsuckers (Sphyrapicus nuchalis) are functionally important because they create sapwells and cavities that other species use for food and nesting. Red‐naped sapsucker ecology within aspen (Populus tremuloides) has been well studied, but relatively little is known about red‐naped sapsuckers in conifer forests. We used light detection and ranging (LiDAR) data to examine occupancy patterns of red‐naped sapsuckers in a conifer‐dominated system. We surveyed for sapsuckers at 162 sites in northern Idaho, USA, during 2009 and 2010. We used occupancy models and an information‐theoretic approach to model sapsucker occupancy as a function of four LiDAR‐based metrics that characterized vegetation structure and tree harvest, and one non‐LiDAR metric that characterized distance to major roads. We evaluated model support across a range of territory sizes using Akaike's information criterion. Top model support was highest at the 4‐ha extent, which suggested that 4 ha was the most relevant scale describing sapsucker occupancy. Sapsuckers were positively associated with variation of canopy height and harvested area, and negatively associated with shrub and large tree density. These results suggest that harvest regimes and structural diversity of vegetation at moderate extents (e.g., 4 ha) largely influence occurrence of red‐naped sapsuckers in conifer forests. Given the current and projected declines of aspen populations, it will be increasingly important to assess habitat relationships, as well as demographic characteristics, of aspen‐associated species such as red‐naped sapsuckers within conifer‐dominated systems to meet future management and conservation goals.