Search

Understanding the underlying mechanisms and links between genome evolution and adaptive innovations stands as a key goal in evolutionary studies. Poplars, among the world's most widely distributed and cultivated trees, exhibit extensive phenotypic diversity and environmental adaptability. In this study, we present a genus-level super-pangenome comprising 19 Populus genomes, revealing the likely pivotal role of private genes in facilitating local environmental and climate adaptation. Through the integration of pangenomes with transcriptomes, methylomes, and chromatin accessibility mapping, we unveil that the evolutionary trajectories of pangenes and duplicated genes are closely linked to local genomic landscapes of regulatory and epigenetic architectures, notably CG methylation in gene-body regions. Further comparative genomic analyses have enabled the identification of 142 202 structural variants across species that intersect with a significant number of genes and contribute substantially to both phenotypic and adaptive divergence. We have experimentally validated a ∼180-bp presence/absence variant affecting the expression of the CUC2 gene, crucial for leaf serration formation. Finally, we developed a user-friendly web-based tool encompassing the multi-omics resources associated with the Populus super-pangenome (http://www.populus-superpangenome.com). Together, the present pioneering super-pangenome resource in forest trees not only aids in the advancement of breeding efforts of this globally important tree genus but also offers valuable insights into potential avenues for comprehending tree biology., (Copyright © 2024 The Author. Published by Elsevier Inc. All rights reserved.)

Despite the economic, ecological, and scientific importance of the genera Salix L. (willows) and Populus L. (poplars, cottonwoods, and aspens) Salicaceae, 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. Salix subgenera and sections, in particular, have been difficult to classify, with one recent attempt termed a "spectacular failure" due to a speculated radiation of the subgenera Vetrix and Chamaetia. 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 the polyphyly of others. Using an fbranch analysis, 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 the Vetrix-Chamaetia clade in Salix. 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 both Salix and Populus have played important roles in the diversification and diversity in these two genera., (© The Author(s) 2023. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Speciation, the continuous process by which new species form, is often investigated by looking at the variation of nucleotide diversity and differentiation across the genome (hereafter genomic landscapes). A key challenge lies in how to determine the main evolutionary forces at play shaping these patterns. One promising strategy, albeit little used to date, is to comparatively investigate these genomic landscapes as progression through time by using a series of species pairs along a divergence gradient. Here, we resequenced 201 whole-genomes from eight closely related Populus species, with pairs of species at different stages along the divergence gradient to learn more about speciation processes. Using population structure and ancestry analyses, we document extensive introgression between some species pairs, especially those with parapatric distributions. We further investigate genomic landscapes, focusing on within-species (i.e. nucleotide diversity and recombination rate) and among-species (i.e. relative and absolute divergence) summary statistics of diversity and divergence. We observe relatively conserved patterns of genomic divergence across species pairs. Independent of the stage across the divergence gradient, we find support for signatures of linked selection (i.e. the interaction between natural selection and genetic linkage) in shaping these genomic landscapes, along with gene flow and standing genetic variation. We highlight the importance of investigating genomic patterns on multiple species across a divergence gradient and discuss prospects to better understand the evolutionary forces shaping the genomic landscapes of diversity and differentiation., (© 2023 John Wiley & Sons Ltd.)

Long noncoding ribonucleic acids (lncRNAs) play crucial roles in regulating key biological processes; however, our knowledge of lncRNAs' roles in plant adaptive evolution is still limited. Here, we determined the divergence of conserved lncRNAs in closely related poplar species that were either tolerant or sensitive to salt stress by comparative transcriptome analysis. Among the 34,363 identified lncRNAs, ~3% were shared among poplar species with conserved sequences but diversified in their function, copy number, originating genomic region and expression patterns. Further cluster analysis revealed that the conserved lncRNAs showed more similar expression patterns within salt-tolerant poplars (Populus euphratica and P. pruinosa) than between salt-tolerant and salt-sensitive poplars. Among these lncRNAs, the antisense lncRNA lncERF024 was induced by salt and the differentiated expression between salt-sensitive and salt-tolerant poplars. The overexpression of lncERF024 in P. alba var. pyramidalis enhanced poplar tolerance to salt stress. Furthermore, RNA pull-down and RNA-seq analysis showed that numerous candidate genes or proteins associated with stress response and photosynthesis might be involved in salt resistance in PeulncERF024-OE poplars. Altogether, our study provided a novel insight into how the diversification of lncRNA expression contributes to plant adaptation traits and showed that lncERF024 may be involved in the regulation of both gene expression and protein function conferring salt tolerance in Populus., (© The Author(s) 2023. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permission@oup.com.)

Salt and drought impair plant osmotic homeostasis and greatly limit plant growth and development. Plants decrease stomatal aperture to reduce water loss and maintain osmotic homeostasis, leading to improved stress tolerance. Herein, we identified the C 2 H 2 transcription factor gene OSMOTIC STRESS INDUCED C 2 H 2 1 (OSIC1) from Populus alba var. pyramidalis to be induced by salt, drought, polyethylene glycol 6000 (PEG6000) and abscisic acid (ABA). Overexpression of OSIC1 conferred transgenic poplar more tolerance to high salinity, drought and PEG6000 treatment by reducing stomatal aperture, while its mutant generated by the CRISPR/Cas9 system showed the opposite phenotype. Furthermore, OSIC1 directly up-regulates PalCuAOζ in vitro and in vivo, encoding a copper-containing polyamine oxidase, to enhance H 2 O 2 accumulation in guard cells and thus modulates stomatal closure when stresses occur. Additionally, ABA-, drought- and salt-induced PalMPK3 phosphorylates OSIC1 to increase its transcriptional activity to PalCuAOζ. This regulation of OSIC1 at the transcriptional and protein levels guarantees rapid stomatal closure when poplar responds to osmotic stress. Our results revealed a novel transcriptional regulatory mechanism of H 2 O 2 production in guard cells mediated by the OSIC1-PalCuAOζ module. These findings deepen our understanding of how perennial woody plants, like poplar, respond to osmotic stress caused by salt and drought and provide potential targets for breeding., (© 2023 The Authors. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd.)

Genomic divergence with gene flow is very common in both plants and animals. However, divergence and gene flow are two counteracting factors during speciation. Identifying the types of genes that are likely to be introgressed and what genetic factors restrict further effective reproduction of interspecific hybrids is of great interest to biologists. We aimed to address these issues using three related tree species, Populus alba (Pa), P. tremula (Pt), and P. tremuloides (Ps), and the interspecific hybrid of the former two species, P. × canescens (Pc). We collected 105 genomes for these four poplar lineages, including 28 Pa, 38Pt, 21 Ps, and 18 Pc individuals, to reconstruct their evolutionary histories. Our coalescence-based simulations indicated that Pa diverged earliest from Ps and Pt, and asymmetrical gene flow existed between any two lineages, with especially large ancient gene flow occurring between Pa and Pt. The genomic landscape of divergence between pairs of the three species are highly heterogeneous, which may have arisen through both divergent sorting of ancient polymorphisms and ongoing gene flow. We found that extant regions of the genome with introgressed ancestry reduced genetic divergence but elevated recombination rates and accounted for 5.76 % of the total genome. Introgressed genes were functionally associated with stress resistance, including innate immune response, anti-adversity response, and programmed cell death. However, candidate genes underlying postmating barriers of Pc were homozygous and resistant to introgression due to the incompatibility of alleles between loci after hybridization and were associated with endosperm and gamete formation and disease resistance. Our study revealed genomic dynamics during speciation with gene flow and identified regions of the genome that were likely introgressed and adaptive as well as candidate loci responsible for hybrid incompatibility that resulted in the formation of postmating barriers after hybridization., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Plants usually suffer from phosphorus starvation because of the low inorganic phosphate (Pi) status of most soils. To cope with this, plants have evolved an adaptive phosphate starvation response (PSR) which involves both developmental and metabolic changes regulated mainly by PHOSPHATE STARVATION RESPONSE1 (PHR1) and its homologs. Here, we elucidated how perennial woody plants, such as poplars (Populus spp.), respond to low-Pi stress. We first performed RNA-seq analysis of low-Pi-treated poplars and identified PtoWRKY40 is rapidly downregulated and protein degraded after stress. Overexpressing and knocking-down PtoWRKY40 downregulated and upregulated the expression of Pi starvation signaling genes, respectively, such as PHOSPHATE TRANSPORTER1 (PHT1)-type genes and PURPLE ACID PHOSPHATASE genes. PtoWRKY40 bound to the W box in the promoter of several PtoPHT1s and repressed their expression. Moreover, PtoWRKY40 interacted with PtoPHR1-LIKE3 (PtoPHL3), a PHR1 homolog in poplar, to inhibit the latter binding to the P1BS element and thus reduced PtoPHT1s' transcription under Pi-sufficient conditions. However, Pi deficiency decreased PtoWRKY40 abundance and therefore released its inhibition on PHT1s. In conclusion, we have uncovered a PSR mechanism mediated by PtoWRKY40 and PtoPHL3 which regulates Pi content in poplars, deepening our understanding of how poplars adapt to diverse Pi conditions and regulate appropriate responses to maintain Pi homeostasis., (© American Society of Plant Biologists 2022. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Background: Salicaceae species have diverse sex determination systems and frequent sex chromosome turnovers. However, compared with poplars, the diversity of sex determination in willows is poorly understood, and little is known about the evolutionary forces driving their turnover. Here, we characterized the sex determination in two Salix species, S. chaenomeloides and S. arbutifolia, which have an XY system on chromosome 7 and 15, respectively., Results: Based on the assemblies of their sex determination regions, we found that the sex determination mechanism of willows may have underlying similarities with poplars, both involving intact and/or partial homologs of a type A cytokinin response regulator (RR) gene. Comparative analyses suggested that at least two sex turnover events have occurred in Salix, one preserving the ancestral pattern of male heterogamety, and the other changing heterogametic sex from XY to ZW, which could be partly explained by the "deleterious mutation load" and "sexually antagonistic selection" theoretical models. We hypothesize that these repeated turnovers keep sex chromosomes of willow species in a perpetually young state, leading to limited degeneration., Conclusions: Our findings further improve the evolutionary trajectory of sex chromosomes in Salicaceae species, explore the evolutionary forces driving the repeated turnovers of their sex chromosomes, and provide a valuable reference for the study of sex chromosomes in other species., (© 2022. The Author(s).)

Background: Sensing and responding to stresses determine the tolerance of plants to adverse environments. The triploid Chinese white poplar is widely cultivated in North China because of its adaptation to a wide range of habitats including highly saline ones. However, its triploid genome complicates any detailed investigation of the molecular mechanisms underlying its adaptations., Results: We report a haplotype-resolved genome of this triploid poplar and characterize, using reverse genetics and biochemical approaches, a MYB gene, SALT RESPONSIVE MYB TRANSCRIPTION FACTOR (SRMT), which combines NUCLEAR FACTOR Y SUBUNIT C 9 (PtoNF-YC9) and RESPONSIVE TO DESICCATION 26 (PtoRD26), to regulate an ABA-dependent salt-stress response signaling. We reveal that the salt-inducible PtoRD26 is dependent on ABA signaling. We demonstrate that ABA or salt drives PtoNF-YC9 shuttling into the nucleus where it interacts with SRMT, resulting in the rapid expression of PtoRD26 which in turn directly regulates SRMT. This positive feedback loop of SRMT-PtoRD26 can rapidly amplify salt-stress signaling. Interference with either component of this regulatory module reduces the salt tolerance of this triploid poplar., Conclusion: Our findings reveal a novel ABA-dependent salt-responsive mechanism, which is mediated by the PtoNF-YC9-SRMT-PtoRD26 module that confers salt tolerance to this triploid poplar. These genes may therefore also serve as potential and important modification targets in breeding programs., (© 2022. The Author(s).)

Hybridization and resulting introgression are important processes shaping the tree of life and appear to be far more common than previously thought. However, how the genome evolution was shaped by various genetic and evolutionary forces after hybridization remains unresolved. Here we used whole-genome resequencing data of 227 individuals from multiple widespread Populus species to characterize their contemporary patterns of hybridization and to quantify genomic signatures of past introgression. We observe a high frequency of contemporary hybridization and confirm that multiple previously ambiguous species are in fact F1 hybrids. Seven species were identified, which experienced different demographic histories that resulted in strikingly varied efficacy of selection and burdens of deleterious mutations. Frequent past introgression has been found to be a pervasive feature throughout the speciation of these Populus species. The retained introgressed regions, more generally, tend to contain reduced genetic load and to be located in regions of high recombination. We also find that in pairs of species with substantial differences in effective population size, introgressed regions are inferred to have undergone selective sweeps at greater than expected frequencies in the species with lower effective population size, suggesting that introgression likely have higher potential to provide beneficial variation for species with small populations. Our results, therefore, illustrate that demography and recombination have interplayed with both positive and negative selection in determining the genomic evolution after hybridization., (© The Author(s) 2022. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

The abscisic acid (ABA) signalling pathway is involved in the plant response to osmotic stress caused by drought and/or salinity. Although the ABA signalling pathway has been elucidated in Arabidopsis, it remains elusive in woody poplars. In this study, genome-wide analyses of U-box genes in poplars revealed that a U-box E3 ubiquitin ligase gene, PalPUB79, is significantly induced following drought, salinity and ABA signalling. PalPUB79 overexpression enhanced drought tolerance in transgenic poplars, while PalPUB79 RNAi lines were more sensitive to drought. PalPUB79 positively regulated ABA signalling pathway. Furthermore, PalPUB79 interacted with PalWRKY77, a negative transcriptional regulator of ABA signalling, and mediated its ubiquitination for degradation, therefore counteracting its inhibitory effect on PalRD26 transcription. However, the finding that PalWRKY77 negatively regulates PalPUB79 expression was indicative of a negative feedback loop between PalWRKY77 and PalPUB79 during ABA signalling in poplar. These findings provide novel insight into the mechanism through which PalPUB79 enhances the ABA-mediated stress response in woody poplars., (© 2021 The Authors. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd.)

Background: Plant secondary growth depends on the activity of the vascular cambium, which produces xylem and phloem. Wood derived from xylem is the most abundant form of biomass globally and has played key socio-economic and subsistence roles throughout human history. However, despite intensive study of vascular development, the full diversity of cell types and the gene networks engaged are still poorly understood., Results: Here, we have applied an optimized protoplast isolation protocol and RNA sequencing to characterize the high-resolution single-cell transcriptional landscape of highly lignified poplar stems. We identify 20 putative cell clusters with a series of novel cluster-specific marker genes and find that these cells are highly heterogeneous based on the transcriptome. Analysis of these marker genes' expression dynamics enables reconstruction of the cell differentiation trajectories involved in phloem and xylem development. We find that different cell clusters exhibit distinct patterns of phytohormone responses and emphasize the use of our data to predict potential gene redundancy and identify candidate genes related to vascular development in trees., Conclusions: These findings establish the transcriptional landscape of major cell types of poplar stems at single-cell resolution and provide a valuable resource for investigating basic principles of vascular cell specification and differentiation in trees., (© 2021. The Author(s).)

High salinity, one of the most widespread abiotic stresses, inhibits photosynthesis, reduces vegetation growth, blocks respiration and disrupts metabolism in plants. In order to survive their long-term lifecycle, trees, such as Populus species, recruit the abscisic acid (ABA) signaling pathway to adapt to a saline environment. However, the molecular mechanism behind the ABA-mediated salt stress response in woody plants remains elusive. We have isolated a WRKY transcription factor gene, PalWRKY77, from Populus alba var. pyramidalis (poplar), the expression of which is repressed by salt stress. PalWRKY77 decreases salt tolerance in poplar. Furthermore, PalWRKY77 negatively regulated ABA-responsive genes and relieved ABA-mediated growth inhibition, indicating that PalWRKY77 is a repressor of the ABA response. In vivo and in vitro assays revealed that PalWRKY77 targets the ABA- and salt-induced PalNAC002 and PalRD26 genes by binding to the W-boxes in their promoters. In addition, overexpression of both PalNAC002 and PalRD26 could elevate salt tolerance in transgenic poplars. These findings reveal a novel negative regulation mechanism for the ABA signaling pathway mediated by PalWRKY77 that results in more sensitivity to salt stress in poplar. This deepens our understanding of the complex responses of woody species to salt stress., (© 2020 Society for Experimental Biology and John Wiley & Sons Ltd.)

Almost all plants in the genus Populus are dioecious (i.e. trees are either male or female), but it is unknown whether dioecy evolved in a common ancestor or independently in different subgenera. Here, we sequence the small peritelomeric X- and Y-linked regions of P. deltoides chromosome XIX. Two genes are present only in the Y-linked region. One is a duplication of a non-Y-linked, female-specifically expressed response regulator, which produces siRNAs that block this gene's expression, repressing femaleness. The other is an LTR/Gypsy transposable element family member, which generates long non-coding RNAs. Overexpression of this gene in A. thaliana promotes androecium development. We also find both genes in the sex-determining region of P. simonii, a different poplar subgenus, which suggests that they are both stable components of poplar sex-determining systems. By contrast, only the duplicated response regulator gene is present in the sex-linked regions of P. davidiana and P. tremula. Therefore, findings in our study suggest dioecy may have evolved independently in different poplar subgenera.

Many recent studies have addressed the mechanisms operating during the early stages of speciation, but surprisingly few studies have tested theoretical predictions on the evolution of strong reproductive isolation (RI). To help address this gap, we first undertook a quantitative review of the hybrid zone literature for flowering plants in relation to reproductive barriers. Then, using Populus as an exemplary model group, we analysed genome-wide variation for phylogenetic tree topologies in both early- and late-stage speciation taxa to determine how these patterns may be related to the genomic architecture of RI. Our plant literature survey revealed variation in barrier complexity and an association between barrier number and introgressive gene flow. Focusing on Populus, our genome-wide analysis of tree topologies in speciating poplar taxa points to unusually complex genomic architectures of RI, consistent with earlier genome-wide association studies. These architectures appear to facilitate the 'escape' of introgressed genome segments from polygenic barriers even with strong RI, thus affecting their relationships with recombination rates. Placed within the context of the broader literature, our data illustrate how phylogenomic approaches hold great promise for addressing the evolution and temporary breakdown of RI during late stages of speciation. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Species are becoming extinct at unprecedented rates as a consequence of human activity. Hence it is important to understand the evolutionary dynamics of species with already small population sizes. Populus ilicifolia is a vulnerable poplar species that is isolated from other poplar species and is uniquely adapted to the Tropics. It has a very limited size, reproduces partly clonally and is therefore an excellent case study for conservation genomics. We present here the first annotated draft genome of P. ilicifolia, characterize genome-wide patterns of polymorphisms and compare those to other poplar species with larger natural ranges. P. ilicifolia experienced a more prolonged and severe decline of effective population size (N e ) and signs of genetic erosion than any other poplar species with which it was compared. At present, the species has the lowest genome-wide genetic diversity, the highest abundance of long runs of homozygosity, high inbreeding levels as well as a high overall accumulation of deleterious variants. However, more effective purging of severely deleterious variants and adaptation to the Tropics may have contributed to its survival. Hence, in spite of its limited genetic variation, it is certainly worth pursuing the conservation efforts of this unique species., (© 2020 Society for Experimental Biology and John Wiley & Sons Ltd.)

Salinity restricts the growth of trees to varying extents, but the regulatory mechanisms involved in their varying salt tolerance are largely unknown. In an effort to elucidate these mechanisms, we identified a total of 99 genes in the Ethylene Responsive Factor (ERF) family of transcription factors and examined their expression patterns under salt stress in Populus alba var. pyramidalis. We found that a B4 group gene, PalERF109, was rapidly induced by salt treatment and preferentially expressed in stems and petioles, where it is probably involved in transport of ions and water in xylem. Overexpression of PalERF109 enhanced the salt tolerance of the poplar, and further analysis showed that it directly upregulated a high-affinity K+transporter (HKT) gene, PalHKT1;2. The results clearly indicate that PalERF109 enhances salt tolerance at least partially through direct activation of PalHKT1;2 and extends understanding of the roles of ERF genes in tree stress responses., (© The Author(s) 2020. Published by Oxford University Press. All rights reserved. For permissions, please email: journals.permission@oup.com.)

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., (© 2020 John Wiley & Sons Ltd.)

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., (© 2020 John Wiley & Sons Ltd.)

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., (© 2019 The Authors. New Phytologist © 2019 New Phytologist Trust.)

How genome divergence eventually leads to speciation is a topic of prime evolutionary interest. Genomic islands of elevated divergence are frequently reported between diverging lineages, and their size is expected to increase with time and gene flow under the speciation-with-gene-flow model. However, such islands can also result from divergent sorting of ancient polymorphisms, recent ecological selection regardless of gene flow, and/or recurrent background selection and selective sweeps in low-recombination regions. It is challenging to disentangle these nonexclusive alternatives, but here we attempt to do this in an analysis of what drove genomic divergence between four lineages comprising a species complex of desert poplar trees. Within this complex we found that two morphologically delimited species, Populus euphratica and Populus pruinosa , were paraphyletic while the four lineages exhibited contrasting levels of gene flow and divergence times, providing a good system for testing hypotheses on the origin of divergence islands. We show that the size and number of genomic islands that distinguish lineages are not associated with either rate of recent gene flow or time of divergence. Instead, they are most likely derived from divergent sorting of ancient polymorphisms and divergence hitchhiking. We found that highly diverged genes under lineage-specific selection and putatively involved in ecological and morphological divergence occur both within and outside these islands. Our results highlight the need to incorporate demography, absolute divergence measurement, and gene flow rate to explain the formation of genomic islands and to identify potential genomic regions involved in speciation., Competing Interests: The authors declare no conflict of interest., (Copyright © 2018 the Author(s). Published by PNAS.)

Populus pruinosa is a large tree that grows in deserts and shows distinct differences in both morphology and adaptation compared to its sister species, P. euphratica. Here we present a draft genome sequence for P. pruinosa and examine genomic variations between the 2 species. A total of 60 Gb of clean reads from whole-genome sequencing of a P. pruinosa individual were generated using the Illumina HiSeq2000 platform. The assembled genome is 479.3 Mb in length, with an N50 contig size of 14.0 kb and a scaffold size of 698.5 kb; 45.47% of the genome is composed of repetitive elements. We predicted 35 131 protein-coding genes, of which 88.06% were functionally annotated. Gene family clustering revealed 224 unique and 640 expanded gene families in the P. pruinosa genome. Further evolutionary analysis identified numerous genes with elevated values for pairwise genetic differentiation between P. pruinosa and P. euphratica. We provide the genome sequence and gene annotation for P. pruinosa. A large number of genetic variations were recovered by comparison of the genomes between P. pruinosa and P. euphratica. These variations will provide a valuable resource for studying the genetic bases for the phenotypic and adaptive divergence of the 2 sister species., (© The Authors 2017. Published by Oxford University Press.)

Background: The factors that contribute to and maintain hybrid zones between distinct species are highly variable, depending on hybrid origins, frequencies and fitness. In this study, we aimed to examine genetic origins, compositions and possible maintenance of Populus × jrtyschensis, an assumed natural hybrid between two distantly related species. This hybrid poplar occurs mainly on the floodplains along the river valleys between the overlapping distributions of the two putative parents., Results: We collected 566 individuals from 45 typical populations of P. × jrtyschensis, P. nigra and P. laurifolia. We genotyped them based on the sequence variations of one maternally inherited chloroplast DNA (cpDNA) fragment and genetic polymorphisms at 20 SSR loci. We further sequenced eight nuclear genes for 168 individuals from 31 populations. Two groups of cpDNA haplotypes characteristic of P. nigra and P. laurifolia respectively were both recovered for P. × jrtyschensis. Genetic structures and coalescent tests of two sets of nuclear population genetic data suggested that P. × jrtyschensis originated from hybridizations between the two assumed parental species. All examined populations of P. × jrtyschensis comprise mainly F1 hybrids from interspecific hybridizations between P. nigra and P. laurifolia. In the habitats of P. × jrtyschensis, there are lower concentrations of soil nitrogen than in the habitats occupied by the other two species., Conclusions: Our extensive examination of the genetic composition of P. × jrtyschensis suggested that it is typical of F1-dominated hybrid zones. This finding plus the low concentration of soil nitrogen in the floodplain soils support the F1-dominated bounded hybrid superiority hypothesis of hybrid zone maintenance for this particular hybrid poplar.

Background: Soil salinity is a significant factor that impairs plant growth and agricultural productivity, and numerous efforts are underway to enhance salt tolerance of economically important plants. Populus species are widely cultivated for diverse uses. Especially, they grow in different habitats, from salty soil to mesophytic environment, and are therefore used as a model genus for elucidating physiological and molecular mechanisms of stress tolerance in woody plants., Description: The Salinity Tolerant Poplar Database (STPD) is an integrative database for salt-tolerant poplar genome biology. Currently the STPD contains Populus euphratica genome and its related genetic resources. P. euphratica, with a preference of the salty habitats, has become a valuable genetic resource for the exploitation of tolerance characteristics in trees. This database contains curated data including genomic sequence, genes and gene functional information, non-coding RNA sequences, transposable elements, simple sequence repeats and single nucleotide polymorphisms information of P. euphratica, gene expression data between P. euphratica and Populus tomentosa, and whole-genome alignments between Populus trichocarpa, P. euphratica and Salix suchowensis. The STPD provides useful searching and data mining tools, including GBrowse genome browser, BLAST servers and genome alignments viewer, which can be used to browse genome regions, identify similar sequences and visualize genome alignments. Datasets within the STPD can also be downloaded to perform local searches., Conclusions: A new Salinity Tolerant Poplar Database has been developed to assist studies of salt tolerance in trees and poplar genomics. The database will be continuously updated to incorporate new genome-wide data of related poplar species. This database will serve as an infrastructure for researches on the molecular function of genes, comparative genomics, and evolution in closely related species as well as promote advances in molecular breeding within Populus. The STPD can be accessed at http://me.lzu.edu.cn/stpd/ .

Background: Populus euphratica Oliv and P. pruinosa Schrenk (Salicaceae) both grow in dry desert areas with high summer temperatures. However, P. euphratica is distributed in dry deserts with deep underground water whereas P. pruinosa occurs in deserts in which there is underground water close to the surface. We therefore hypothesized that these two sister species may have evolved divergent regulatory and metabolic pathways during their interaction with different salt habitats and other stresses. To test this hypothesis, we compared transcriptomes from callus exposed to 24 h of salt stress and control callus samples from both species and identified differentially expressed genes (DEGs) and alternative splicing (AS) events that had occurred under salt stress., Results: A total of 36,144 transcripts were identified and 1430 genes were found to be differentially expressed in at least one species in response to salt stress. Of these DEGs, 884 and 860 were identified in P. euphratica and P. pruinosa, respectively, while 314 DEGs were common to both species. On the basis of parametric analysis of gene set enrichment, GO enrichment in P. euphratica was found to be significantly different from that in P. pruinosa. Numerous genes involved in hormone biosynthesis, transporters and transcription factors showed clear differences between the two species in response to salt stress. We also identified 26,560 AS events which were mapped to 8380 poplar genomic loci from four libraries. GO enrichments for genes undergoing AS events in P. euphratica differed significantly from those in P. pruinosa., Conclusions: A number of salt-responsive genes in both P. euphratica and P. pruinosa were identified and candidate genes with potential roles in the salinity adaptation were proposed. Transcriptome comparisons of two sister desert poplar species under salt stress suggest that these two species may have developed different genetic pathways in order to adapt to different desert salt habitats. The DEGs that were found to be common to both species under salt stress may be especially important for future genetic improvement of cultivated poplars or other crops through transgenic approaches in order to increase tolerance of saline soil conditions.

Understanding which genes have evolved rapidly with the recent tree speciation in arid habitats can provide valuable insights into different adaptation mechanisms. We employed a comparative evolutionary analysis of expressed sequence tags (ESTs) from two desert poplars, Populus pruinosa and P. euphratica, which diverged in the recent past. Following an approach taken previously with P. euphratica, we conducted a deep transcriptomic analysis of P. pruinosa. To maximize representation of conditional transcripts, mRNA was obtained from living tissues of two types of callus and desert-grown trees. De novo assembly generated 114,866 high-quality unique sequences using Solexa sequence data. Following assembly we were able to identify, with high confidence, 2859 orthologous sequence pairs between the two species. Based on the ratio of nonsynonymous (Ka) to synonymous (Ks) substitutions, we identified a total of 84 (2.9%) ortholog pairs exhibiting rapid evolution with signs of strong selection (Ka/Ks>1). Genes homologous to these ortholog pairs in model species are mainly involved in 'responses to stress', 'ubiquitin-dependent protein catabolic processes', and 'biological regulation'. Finally, we examined the expression patterns of candidate genes with rapid evolution in response to salt stress. Only one pair of orthologs up-regulated their expression in both species while three and four genes were found to up-regulated in P. pruinosa and in P. euphratica respectively. Our findings together suggest that the genes at the same category or network but with differentiated expressions or functions may have evolved rapidly during adaptive divergence of the two species to differentiated salty desert habitats.

Despite the high economic and ecological importance of forests, our knowledge of the genomic evolution of trees under salt stress remains very limited. Here we report the genome sequence of the desert poplar, Populus euphratica, which exhibits high tolerance to salt stress. Its genome is very similar and collinear to that of the closely related mesophytic congener, P. trichocarpa. However, we find that several gene families likely to be involved in tolerance to salt stress contain significantly more gene copies within the P. euphratica lineage. Furthermore, genes showing evidence of positive selection are significantly enriched in functional categories related to salt stress. Some of these genes, and others within the same categories, are significantly upregulated under salt stress relative to their expression in another salt-sensitive poplar. Our results provide an important background for understanding tree adaptation to salt stress and facilitating the genetic improvement of cultivated poplars for saline soils.

Populus euphratica is well-adapted to extreme desert environments and is an important model species for studying the effects of abiotic stresses on trees. Here we present the first deep transcriptomic analysis of this species. To maximize representation of conditional transcripts, mRNA was obtained from living tissues of desert-grown trees and two types of callus (salt-stressed and unstressed). De novo assembly generated 86,777 Unigenes using Solexa sequence data. These sequences covered 92% of previously reported P. euphratica expressed sequence tags (ESTs) and 90% of the TIGR poplar ESTs, and a total of 58,499 high-quality unique sequences were annotated by BLAST similarity searches against public databases. We found that 27% of the total Unigenes were differentially expressed (up- or down-regulated) in response to salt stress in P. euphratica callus. These differentially expressed genes are mainly involved in transport, transcription, cellular communication and metabolism. In addition, we found that numerous putative genes involved in ABA regulation and biosynthesis were also differentially regulated. This study represents the deepest transcriptomic and gene-annotation analysis of P. euphratica to date. The genetic knowledge acquired should be very useful for future studies of the molecular adaptation of this tree species to abiotic stress and facilitate genetic manipulation of other poplar species.

Background: The fixed genetic differences between ecologically divergent species were found to change greatly depending on the markers examined. With such species it is difficult to differentiate between shared ancestral polymorphisms and past introgressions between the diverging species. In order to disentangle these possibilities and provide a further case for DNA barcoding of plants, we examine genetic differentiation between two ecologically divergent poplar species, Populus euphratica Oliver and P. pruinosa Schrenk using three different types of genetic marker., Methodology/principal Findings: We genotyped 290 individuals from 29 allopatric and sympatric populations, using chloroplast (cp) DNA, nuclear (nr) ITS sequences and eight simple sequence repeat (SSR) loci. Three major cpDNA haplotypes were widely shared between the two species and between-species cpDNA differentiation (F(CT)) was very low, even lower than among single species populations. The average SSR F(CT) values were higher. Bayesian clustering analysis of all loci allowed a clear delineation of the two species. Gene flow, determined by examining all SSR loci, was obvious but only slightly asymmetrical. However, the two species were almost fixed for two different nrITS genotypes that had the highest F(CT), although a few introgressed individuals were detected both in allopatric and sympatric populations., Conclusions: The two species shared numerous ancestral polymorphisms at cpDNA and a few SSR loci. Both ITS and a combination of nuclear SSR data could be used to differentiate between the two species. Introgressions and gene flow were obvious between the two species either during or after their divergence. Our findings underscore the complex genetic differentiations between ecologically diverged species and highlight the importance of nuclear DNA (especially ITS) differentiation for delimiting closely related plant species.

Populus euphratica is a salt-tolerant tree species growing in semi-arid saline areas. A Na(+)/H(+) antiporter gene was successfully isolated from this species through RACE cloning, and named PeSOS1. The isolated cDNA was 3665 bp long and contained a 3438 bp open reading frame that was predicted to encode a 127-kDa protein with 12 hypothetical transmembrane domains in the N-terminal part and a long hydrophilic cytoplasmic tail in the C-terminal part. The amino acid sequence of this PeSOS1 gene showed 64% identity with the previously isolated SOS1 gene from the glycophyte Arabidopsis thaliana. The level of protein expressed by PeSOS1 in the leaves of P. euphratica was significantly up-regulated in the presence of high (200 mM) concentrations of NaCl, while the mRNA level in the leaves remained relatively constant. Immunoanalysis suggested that the protein encoded by PeSOS1 is localized in the plasma membrane. Expression of PeSOS1 partially suppressed the salt sensitive phenotypes of the EP432 bacterial strain, which lacks the activity of the two Na(+)/H(+) antiporters EcNhaA and EcNhaB. These results suggest that PeSOS1 may play an essential role in the salt tolerance of P. euphratica and may be useful for improving salt tolerance in other tree species.

Nitric oxide (NO) and hydrogen peroxide (H2O2) function as signalling molecules in plants under abiotic and biotic stresses. Calluses from Populus euphratica, which show salt tolerance, were used to study the interaction of NO and H2O2 in plant adaptation to salt resistance. The nitric oxide synthase (NOS) activity was identified in the calluses, and this activity was induced under 150 mM NaCl treatment. Under 150 mM NaCl treatment, the sodium (Na) percentage decreased, but the potassium (K) percentage and the K/Na ratio increased in P. euphratica calluses. Application of glucose/glucose oxidase (G/GO, a H2O2 donor) and sodium nitroprusside (SNP, a NO donor) revealed that both H2O2 and NO resulted in increased K/Na ratio in a concentration-dependent manner. Diphenylene iodonium (DPI, an NADPH oxidase inhibitor) counteracted H2O2 and NO effect by increasing the Na percentage, decreasing the K percentage and K/Na ratio. NG-monomethyl-L-Arg monoacetate (NMMA, an NO synthase inhibitor) and 2-phenyl-4,4,5,5-tetramethyl-imidazoline-1-oxyl-3-oxyde (PTIO, a specific NO scavenger) only reversed NO effect, but did not block H2O2 effect. The increased activity of plasma membrane (PM) H+ -ATPase caused by salt stress was reversed by treatment with DPI and NMMA. Exogenous H2O2 increased the activity of PM H+ -ATPase, but the effect could not be diminished by NMMA and PTIO. The NO-induced increase of PM H+ -ATPase can be reversed by NMMA and PTIO, but not by DPI. Western blot analysis demonstrated that NO and H2O2 stimulated the expression of PM H+ -ATPase in P. euphratica calluses. These results indicate that NO and H2O2 served as intermediate molecules in inducing salt resistance in the calluses from P. euphratica under slat stress by increasing the K/Na ratio, which was dependent on the increased PM H+ -ATPase activity.

Dioecy, the presence of separate sexes on distinct individuals, has evolved repeatedly in multiple plant lineages. However, the specific mechanisms by 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 chromosomes in flowering plants by showing how transposition and rearrangement of a single gene can control sex in both XY and ZW systems. [ABSTRACT FROM AUTHOR]