Recent Y chromosome divergence despite ancient origin of dioecy in poplars (Populus)

All species of the genus Populus (poplar, aspen) are dioecious, suggesting an ancient origin of this trait. Despite some empirical counter examples, theory suggests that nonrecombining sex‐linked regions should quickly spread, eventually becoming heteromorphic chromosomes. In contrast, we show using whole‐genome scans that the sex‐associated region in Populus trichocarpa is small and much younger than the age of the genus. This indicates that sex determination is highly labile in poplar, consistent with recent evidence of ‘turnover’ of sex‐determination regions in animals. We performed whole‐genome resequencing of 52 P. trichocarpa (black cottonwood) and 34 Populus balsamifera (balsam poplar) individuals of known sex. Genomewide association studies in these unstructured populations identified 650 SNPs significantly associated with sex. We estimate the size of the sex‐linked region to be ~100 kbp. All SNPs significantly associated with sex were in strong linkage disequilibrium despite the fact that they were mapped to six different chromosomes (plus 3 unmapped scaffolds) in version 2.2 of the reference genome. We show that this is likely due to genome misassembly. The segregation pattern of sex‐associated SNPs revealed this to be an XY sex‐determining system. Estimated divergence times of X and Y haplotype sequences (6–7 Ma) are much more recent than the divergence of P. trichocarpa (poplar) and Populus tremuloides (aspen). Consistent with this, in P. tremuloides, we found no XY haplotype divergence within the P. trichocarpa sex‐determining region. These two species therefore have a different genomic architecture of sex, suggestive of at least one turnover event in the recent past.     

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.     

History rather than hybridization determines population structure and adaptation in Populus balsamifera

Hybridization between species is known to greatly affect their genetic diversity and, therefore, their evolution. Also, within species, there may be genetic clusters between which gene flow is limited, which may impact natural selection. However, few studies have looked simultaneously at the influence of among‐species and within‐species gene flow. Here, we study the influence of hybridization between Populus balsamifera and Populus trichocarpa on population structure and adaptation in P. balsamifera. We did this by sampling a total of 1517 individuals from across the ranges of these two species, and by genotyping them using a combination of 93 nuclear and 17 cpDNA SNPs. We found that hybridization is mostly limited to the contact zone where the species’ distributions overlap. Within P. balsamifera, we found multiple levels of population structure. Interestingly, the border between the Eastern and Central clusters is very sharp, whereas the border between the Central and Western clusters is diffuse. Outlier analysis revealed that three loci associated with the sharp border were also associated with climate. We hypothesize that the observed clusters derive from three refugia during the Pleistocene ice ages. Between the Central and Western clusters, post‐glacial long‐distance gene flow has led to the diffusion of their border. In the Eastern cluster, we hypothesize that endogenous genomic barriers have developed, leading to the sharp border and a spurious climate association. We conclude that the large‐scale genetic structure of P. balsamifera is mostly shaped by historical factors and the influence of interspecific hybridization is limited.     

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.     

Riverscape-level patterns of riparian plant diversity along a successional gradient, Willamette river, Oregon

Riparian forest communities dominated by Populus balsamifera ssp. trichocarpa L. (Torr. and Gray ex Hook.) Brayshaw are important contributors to biodiversity in terrestrial and aquatic ecosystems of the Western United States. Species composition along a successional gradient from stand initiation to late-succession of P. balsamifera-dominated riparian forests was investigated along 145 km of the Willamette River, Oregon. There were 151 total species encountered across 28 stands and a mean species richness of 33.3 species per stand. Young stands were dominated by P. balsamifera and Salix tree spp. and opportunistic herbaceous species. Understory trees, shrubs, and herbaceous species as well as late-successional tree species established 12–15 years after stand initiation. Fraxinus latifolia Benth. was the dominant late-successional tree species. Vertical structural diversity, P. balsamifera mean diameter at breast height, large tree biomass, and stand age were strongly correlated with understory species presence and abundance based on non-metric multidimensional scaling (NMS) ordination. There were no young stands on mid and high terraces and this was reflected in geomorphic position being strongly correlated with the stand age gradient. Abundance of Phalaris arundinacea L. an invasive grass species, was also significantly correlated with plant species composition and abundance. This study indicates that Willamette River riparian forests are diverse and therefore important to the biodiversity of the Willamette River valley and that their presence as a mosaic of communities of different successional stages may be threatened by human interventions, including influences exerted by introduced plant species.     

Comparative sequence analyses indicate that Coffea (Asterids) and Vitis (Rosids) derive from the same paleo-hexaploid ancestral genome

The complete sequence of Vitis vinifera revealed that the rosid clade derives from a hexaploid ancestor. At present, no analysis of complete genome sequence is available for an asterid, the other large eudicot clade, which includes the economically important species potato, tomato and coffee. To elucidate the genomic history of asterids, we compared the sequence of an 800 kb region of diploid Coffea genome to the orthologous regions of V. vinifera, Populus trichocarpa and Arabidopsis thaliana. We found a very high level of collinearity between around 80 genes of the three rosid species and Coffea. Collinearity comparisons between orthologous and paralogous regions indicates that (1) the Coffea (and consequently all asterids) and rosids share the same hexaploid ancestor; (2) the diploidization process (loss of duplicated and redundant copies from the whole genome duplication) was very advanced in the most recent common ancestor of rosids and asterids. Finally, no additional polyploidization events were detected in the Coffea lineage. Differences in gene loss rates were detected among the three rosid species and linked to the divergence in protein sequences.     

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.     

Species‐specific microsatellite markers to monitor gene flow between exotic poplars and their natural relatives in eastern North America

Species‐specific microsatellite markers were obtained for the unambiguous recognition of five poplar species of ecological and commercial importance to eastern North America: the native species Populus balsamifera and Populus deltoides, the exotic species Populus maximowiczii, Populus nigra, Populus trichocarpa and their interspecific hybrids. Forty‐four of 71 tested primer pairs amplified simple sequence repeat (SSR) loci for all five taxa. Six of these loci showed non‐overlapping allelic diversity between species, including fixed differences. Together, they were useful to identify unambiguously the five taxa and to validate parental contributions in a group of hybrid progeny. These markers will be invaluable to detect gene flow from plantations of exotic poplar into adjacent stands of native species and between the two potentially hybridizing native species P. balsamifera and P. deltoides.     

Fitness dynamics within a poplar hybrid zone: I. Prezygotic and postzygotic barriers impacting a native poplar hybrid stand

Hybridization and introgression are pervasive evolutionary phenomena that provide insight into the selective forces that maintain species boundaries, permit gene flow, and control the direction of evolutionary change. Poplar trees (Populus L.) are well known for their ability to form viable hybrids and maintain their distinct species boundaries despite this interspecific gene flow. We sought to quantify the hybridization dynamics and postzygotic fitness within a hybrid stand of balsam poplar (Populus balsamifera L.), eastern cottonwood (P. deltoides Marsh.), and their natural hybrids to gain insight into the barriers maintaining this stable hybrid zone. We observed asymmetrical hybrid formation with P. deltoides acting as the seed parent, but with subsequent introgression biased toward P. balsamifera. Native hybrids expressed fitness traits intermediate to the parental species and were not universally unfit. That said, native hybrid seedlings were absent from the seedling population, which may indicate additional selective pressures controlling their recruitment. It is imperative that we understand the selective forces maintaining this native hybrid zone in order to quantify the impact of exotic poplar hybrids on this native system.     

Fitness dynamics within a poplar hybrid zone: II. Impact of exotic sex on native poplars in an urban jungle

Trees bearing novel or exotic gene components are poised to contribute to the bioeconomy for a variety of purposes such as bioenergy production, phytoremediation, and carbon sequestration within the forestry sector, but sustainable release of trees with novel traits in large‐scale plantations requires the quantification of risks posed to native tree populations. Over the last century, exotic hybrid poplars produced through artificial crosses were planted throughout eastern Canada as ornamentals or windbreaks and these exotics provide a proxy by which to examine the fitness of exotic poplar traits within the natural environment to assess risk of exotic gene escape, establishment, and spread into native gene pools. We assessed postzygotic fitness traits of native and exotic poplars within a naturally regenerated stand in eastern Canada (Quebec City, QC). Pure natives (P. balsamifera and P. deltoides spp. deltoides), native hybrids (P. deltoides × P. balsamifera), and exotic hybrids (trees bearing Populus nigra and P. maximowiczii genetic components) were screened for reproductive biomass, yield, seed germination, and fungal disease susceptibility. Exotic hybrids expressed fitness traits intermediate to pure species and were not significantly different from native hybrids. They formed fully viable seed and backcrossed predominantly with P. balsamifera. These data show that exotic hybrids were not unfit and were capable of establishing and competing within the native stand. Future research will seek to examine the impact of exotic gene regions on associated biotic communities to fully quantify the risk exotic poplars pose to native poplar forests.     

Whole-genome sequencing uncovers the structural and transcriptomic landscape of hexaploid wheat/Ambylopyrum muticum introgression lines

Wheat is a globally vital crop, but its limited genetic variation creates a challenge for breeders aiming to maintain or accelerate agricultural improvements over time. Introducing novel genes and alleles from wheat's wild relatives into the wheat breeding pool via introgression lines is an important component of overcoming this low variation but is constrained by poor genomic resolution and limited understanding of the genomic impact of introgression breeding programmes. By sequencing 17 hexaploid wheat/Ambylopyrum muticum introgression lines and the parent lines, we have precisely pinpointed the borders of introgressed segments, most of which occur within genes. We report a genome assembly and annotation of Am. muticum that has facilitated the identification of Am. muticum resistance genes commonly introgressed in lines resistant to stripe rust. Our analysis has identified an abundance of structural disruption and homoeologous pairing across the introgression lines, likely caused by the suppressed Ph1 locus. mRNAseq analysis of six of these introgression lines revealed that novel introgressed genes are rarely expressed and those that directly replace a wheat orthologue have a tendency towards downregulation, with no discernible compensation in the expression of homoeologous copies. This study explores the genomic impact of introgression breeding and provides a schematic that can be followed to characterize introgression lines and identify segments and candidate genes underlying the phenotype. This will facilitate more effective utilization of introgression pre-breeding material in wheat breeding programmes.     

Parallel introgression and selection on introduced alleles in a native species

As humans cause the redistribution of species ranges, hybridization between previously allopatric species is on the rise. Such hybridization can have complex effects on overall fitness of native species as new allelic combinations are tested. Widespread species introductions provide a unique opportunity to study selection on introgressed alleles in independent, replicated populations. We examined selection on alleles that repeatedly introgressed from introduced rainbow trout (Oncorhynchus mykiss) into native westslope cutthroat trout (Oncorhynchus clarkii lewisi) populations in western Canada. We found that the degree of introgression of individual single nucleotide polymorphisms from the invasive species into the native is correlated between independent watersheds. A number of rainbow trout alleles have repeatedly swept to high frequency in native populations, suggesting parallel adaptive advantages. Using simulations, we estimated large selection coefficients up to 0.05 favoring several rainbow trout alleles in the native background. Although previous studies have found reduced hybrid fitness and genome‐wide resistance to introgression in westslope cutthroat trout, our results suggest that some introduced genomic regions are strongly favored by selection. Our study demonstrates the utility of replicated introductions as case studies for understanding parallel adaptation and the interactions between selection and introgression across the genome. We suggest that understanding this variation, including consideration of beneficial alleles, can inform management strategies for hybridizing species.     

The effects of spacing on growth, morphology and biomass production and allocation in two hybrid poplar clones growing in the boreal region of Canada

Intra-clonal competition was studied in young hybrid poplar plantations to assess the effects of spacing on growth, biomass production and allocation, and morphological characteristics of above- and below-ground tree parts. Three spacings were used as whole-plots (1 × 1 m, 3 × 3 m and 5 × 5 m), with two hybrid poplar clones as subplots (BT747, Populus balsamifera L. × P. trichocarpa Torr. & Gray; MB915, P. maximowiczii A. Henry × P. balsamifera L.) in a split-plot design. After six growing seasons, diameter at breast height (dbh) increased by about 120% from the 1 × 1 m to the 5 × 5 m spacing for clone MB915, while there was no significant change in dbh for the other clone. The effect of spacing on height growth was opposite for the clones; it increased by about 175% from the narrowest to the widest spacing for clone MB915, while it decreased by about 27% for clone BT747. Estimates of above-ground biomass production after six growing seasons were significantly reduced with increasing spacing, with 29.6, 4.9 and 3.2 MgDM ha⁻¹ on average from the narrowest to the widest spacing. Branch traits and the vertical distribution of leaf area were the most affected by spacing for both clones, while live crown ratio and percentage of syllepsis did not change. Spacing also affected proportions of biomass allocated to stem, leaves, and branches, but allocation to roots did not change.     

Parallel plumage colour evolution and introgressive hybridization in wheatears

Genetic and phenotypic mosaics, in which various phenotypes and different genomic regions show discordant patterns of species or population divergence, offer unique opportunities to study the role of ancestral and introgressed genetic variation in phenotypic evolution. Here, we investigated the evolution of discordant phenotypic and genetic divergence in a monophyletic clade of four songbird taxa—pied wheatear (O. pleschanka), Cyprus wheatear (Oenanthe cypriaca), and western and eastern subspecies of black‐eared wheatear (O. h. hispanica and O. h. melanoleuca). Phenotypically, black back and neck sides distinguish pied and Cyprus wheatears from the white‐backed/necked black‐eared wheatears. Meanwhile, mitochondrial variation only distinguishes western black‐eared wheatear. In the absence of nuclear genetic data, and given frequent hybridization among eastern black‐eared and pied wheatear, it remains unclear whether introgression is responsible for discordance between mitochondrial divergence patterns and phenotypic similarities, or whether plumage coloration evolved in parallel. Multispecies coalescent analyses of about 20,000 SNPs obtained from RAD data mapped to a draft genome assembly resolve the species tree, provide evidence for the parallel evolution of colour phenotypes and establish western and eastern black‐eared wheatears as independent taxa that should be recognized as full species. The presence of the entire admixture spectrum in the Iranian hybrid zone and the detection of footprints of introgression from pied into eastern black‐eared wheatear beyond the hybrid zone despite strong geographic structure of ancestry proportions furthermore suggest a potential role for introgression in parallel plumage colour evolution. Our results support the importance of standing heterospecific and/or ancestral variation in phenotypic evolution.     

When a clonal genome finds its way back to a sexual species: evidence from ongoing but rare introgression in the hybridogenetic water frog complex

Besides several exceptions, asexual metazoans are usually viewed as ephemeral sinks for genomes, which become ‘frozen’ in clonal lineages after their emergence from ancestral sexual species. Here, we investigated whether and at what rate the asexuals are able to introgress their genomes back into the parental sexual population, thus more or less importantly affecting the gene pools of sexual species. We focused on hybridogenetic hybrids of western Palaearctic water frogs (Pelophylax esculentus), which originate through hybridization between P. ridibundus and P. lessonae, but transmit only clonal ridibundus genome into their gametes. Although usually mating with P. lessonae, P. esculentus may upon mating with P. ridibundus or another hybrid produce sexually reproducing P. ridibundus offspring with the introgressed ex‐clonal genome. We compared the rate of nuclear amplified fragment length polymorphism (AFLP) and mitochondrial introgression in two types of populations, that is, those where P. ridibundus occurs in isolation and those where it lives with the hybridogens. Although significant differentiation (Φpt) between sexual and clonal ridibundus genomes suggested limited gene flow between sexuals and hybridogens, a non‐negligible (~5%) proportion of P. ridibundus bore introgressed mtDNA and AFLP markers. Whereas transfer of mtDNA was exclusively unidirectional, introgression of nuclear markers was bidirectional. The proportion of introgressed P. ridibundus was highest in syntopic populations with P. esculentus, proving an ongoing and site‐specific interspecific genetic transfer mediated by hybridogenetic hybrids. It turns out that asexual hybrids are not just a sink for genes of sexual species, but may significantly influence the genetic architecture of their sexual counterparts.     

Heterogeneity and concordance in locus‐specific differentiation and introgression between species of towhees

The maintenance or breakdown of reproductive isolation is an observable outcome of secondary contact between species. In cases where hybrids beyond the F1 are formed, the representation of each species' ancestry can vary dramatically among genomic regions. This genomic heterogeneity in ancestry and introgression can offer insight into evolutionary processes, particularly if introgression is compared in multiple hybrid zones. Similarly, considerable heterogeneity exists across the genome in the extent to which populations and species have diverged, reflecting the combined effects of different evolutionary processes on genetic variation. We studied hybridization across two hybrid zones of two phenotypically well‐differentiated bird species in Mexico (Pipilo maculatus and P. ocai), to investigate genomic heterogeneity in differentiation and introgression. Using genotyping‐by‐sequencing (GBS) and hierarchical Bayesian models, we genotyped 460 birds at over 41 000 single nucleotide polymorphism (SNP) loci. We identified loci exhibiting extreme introgression relative to the genome‐wide expectation using a Bayesian genomic cline model. We also estimated locus‐specific F_ST and identified loci with exceptionally high genetic divergence between the parental species. We found some concordance of locus‐specific introgression in the two independent hybrid zones (6–20% of extreme loci shared across zones), reflecting areas of the genome that experience similar gene flow when the species interact. Additionally, heterogeneity in introgression and divergence across the genome revealed another subset of loci under the influence of locally specific factors. These results are consistent with a history in which reproductive isolation has been influenced by a common set of loci in both hybrid zones, but where local environmental and stochastic factors also lead to genomic differentiation.     

The effect of industrial pollution on copper,lead, zinc and cadmium concentration in xylem rings of resistant (Populus marilandica) and sensitive (P. balsamifera) species of poplar

Summary The patterns of radial distribution of copper, lead, zinc and cadmium in two species of poplar, growing under polluted conditions near a copper smelter, were investigated. Populus marilandica is very resistant to such pollution, but P. balsamifera is sensitive. Soils were examined for concentrations of these metals, and only a surface layer of soil was polluted with copper and lead. P. balsamifera absorbed much more cadmium than P. marilandica under the same conditions. P. marilandica shows a pattern of copper radial distribution with significantly increased concentration of the meral in the layers in the middle of the trunk compared to P. balsamifera. Zinc and cadmium have a pattern with increased metal concentration from the bark towards the pith of the tree. In the case of copper and lead the pattern indicates increased concentrations of these metals in the rings adjacent to the pith while the pattern in other parts of the tree remains unconfirmed due to poor precision. Poplar xylem seems to accumulate zinc and cadmium while the content of lead in xylem is much lower than in the adjacent soil.This paper is dedicated to the memory of Professor Stefan Bialobok, Institute of Dendrology, Polish Academy of Sciences in Kornik     

LANDSCAPE GENOMICS OF POPULUS TRICHOCARPA: THE ROLE OF HYBRIDIZATION,LIMITED GENE FLOW,AND NATURAL SELECTION IN SHAPING PATTERNS OF POPULATION STRUCTURE

Populus trichocarpa is an ecologically important tree across western North America. We used a large population sample of 498 accessions over a wide geographical area genotyped with a 34K Populus SNP array to quantify geographical patterns of genetic variation in this species (landscape genomics). We present evidence that three processes contribute to the observed patterns: (1) introgression from the sister species P. balsamifera, (2) isolation by distance (IBD), and (3) natural selection. Introgression was detected only at the margins of the species’ distribution. IBD was significant across the sampled area as a whole, but no evidence of restricted gene flow was detected in a core of drainages from southern British Columbia (BC). We identified a large number of F_ST outliers. Gene Ontology analyses revealed that F_ST outliers are overrepresented in genes involved in circadian rhythm and response to red/far‐red light when the entire dataset is considered, whereas in southern BC heat response genes are overrepresented. We also identified strong correlations between geoclimate variables and allele frequencies at F_ST outlier loci that provide clues regarding the selective pressures acting at these loci.