Phylogenomics of the genus Populus reveals extensive interspecific gene flow and balancing selection

The population chloroplast genomes of Populus reveal the phylogenetic relationship between three new taxa of sect. Leucoides and their parents

BACKGROUND

Poplars are important woody plants, which are widely distributed in the forests from the subtropics to the north of the Northern Hemisphere. Poplars have high ecological and economic value. However, there are frequent interspecific and intraspecific hybrids in Populus, resulting in a large number of intermediate taxa, which makes the morphological identification of Populus very challenging. Plastid genome is an important tool to study the evolutionary relationship of plants. Therefore, comparison and phylogenetic analysis were carried out based on the population chloroplast genomes of 34 individuals from 7 taxa.

RESULTS

In this study, seven newly assembled and annotated chloroplast genomes of Populus were reported. They all had typical quadripartite structures with the same GC content (37.6%), but there were differences within the population, and the genome size ranged from 155,736 bp to 156,812 bp. In all Populus species, 134 genes were identified, including 88 protein coding genes (PCGs), 37 tRNA and 8 rRNA genes. The gene sequences alignment of different taxa showed that the gene sequences and content were relatively conservative, there was no gene rearrangement, and only 3 highly variable regions (psbZ-trnG, ndhC-trnV and trnN-trnR) were identified, which can be used as molecular markers. Most PCGs had high codon usage bias and 3 positive selection genes (rps7, rps12 and rpl16) have been identified. The analysis of population genetic structure and phylogeny showed that the chloroplast genomes supported that Populus was a monophyletic taxon, which could be divided into four sections (Abaso, Turanga, Populus and ATL (Aigeiros, Tacamahaca and Leucoides)). Among them, P. dafengensis, P. butuoensis and P. szechuanica had the closest genetic relationship, P. gonggaensis and P. cathayana had the closest genetic relationship, it was speculated that the taxa of Sect. Tacamahaca may be the main female parent of the three new taxa from Sect. Leucoides.

CONCLUSION

In general, this study provides valuable insights for new species identification, phylogenetic relationships, breeding and resource development, and genetic diversity of Populus.

Species delimitation in the Populus laurifolia complex (Salicaceae) based on phylogenetic and morphometric evidence

Due to significant morphological differences and extensive interspecific hybridization, there are numerous species complexes with taxonomic challenges in the genus Populus. Integrative taxonomy, which combines evidence of morphology, molecular phylogeny, niche differentiation, and reproductive isolation, provides the most effective approaches for species delimitation. The Populus laurifolia complex, which belongs to Populus subg. Tacamahaca (Salicaceae), is distributed in the Altai Mountains and Tianshan Mountains. This complex exhibits morphological variability, making species delimitation challenging. Due to limited sampling and systematic studies, its taxonomy has remained unresolved. In this study, 337 specimens, along with online digital samples representing nearly all wild populations, were collected. Morphological analyses were performed to evaluate key traits and clarify species boundaries. Phylogenetic relationships were reconstructed using concatenation and coalescent methods based on 566,375 nuclear single-nucleotide polymorphisms (SNPs). Ecological niche differentiation was assessed, and ABBA-BABA analysis was used to examine interspecific hybridization. The results revealed that this complex, based on a series of significant character states, could be morphologically distinguished into three species-P. laurifolia (Populus pilosa considered a synonym of P. laurifolia), Populus talassica, and Populus pamirica-which also correspond to three well-supported clades in the phylogenetic trees. P. pamirica exhibits some degree of ecological niche differentiation from P. talassica and P. laurifolia, whereas the latter two show minimal differentiation. Gene flow within the complex remains limited. This research underscores the importance of integrating multiple lines of evidence in the classification of Populus, providing a framework for future taxonomic studies.

Trans-Specific Polymorphisms Between Cryptic Daphnia Species Affect Fitness and Behavior

Shared polymorphisms, loci with identical alleles across species, are of unique interest in evolutionary biology as they may represent cases of selection maintaining ancient genetic variation post-speciation, or contemporary selection promoting convergent evolution. In this study, we investigate the abundance of shared polymorphism between two members of the Daphnia pulex species complex. We test whether the presence of shared mutations is consistent with the action of balancing selection or alternative hypotheses such as hybridization, incomplete lineage sorting or convergent evolution. We analyzed over 2,000 genomes from six taxa in the D. pulex species group and examined the prevalence and distribution of shared alleles between the focal species pair, North American and European D. pulex. We show that North American and European D. pulex diverged over 10 million years ago, yet retained tens of thousands of shared polymorphisms. We suggest that the number of shared polymorphisms between North American and European D. pulex cannot be fully explained by hybridization or incomplete lineage sorting alone. We show that most shared polymorphisms could be the product of convergent evolution, that a limited number appear to be old trans-specific polymorphisms, and that balancing selection is affecting convergent and ancient mutations alike. Finally, we provide evidence that a blue wavelength opsin gene with trans-specific polymorphisms has functional effects on behavior and fitness in the wild.

Sequencing historical RNA: unrealized potential to increase understanding of the plant tree of life

Recent studies have demonstrated that it is a misconception that transcriptome sequencing requires tissue preserved at ultracold temperatures. Here, we outline the potential origins of this misconception and its possible role in biasing the geographic distribution of published plant transcriptomes. We highlight the importance of ensuring diverse sampling by providing an overview of the questions that transcriptomes can answer about the forces shaping the plant tree of life. We discuss how broadening transcriptome sequencing to include existing specimens will allow the field to grow and more fully utilize biological collections. We hope this article encourages the expansion of the current trend in 'herbariomics' research to include whole-transcriptome sequencing of historical RNA.

Integrating evolutionary genomics of forest trees to inform future tree breeding amid rapid climate change

Global climate change is leading to rapid and drastic shifts in environmental conditions, posing threats to biodiversity and nearly all life forms worldwide. Forest trees serve as foundational components of terrestrial ecosystems and play a crucial and leading role in combating and mitigating the adverse effects of extreme climate events, despite their own vulnerability to these threats. Therefore, understanding and monitoring how natural forests respond to rapid climate change is a key priority for biodiversity conservation. Recent progress in evolutionary genomics, driven primarily by cutting-edge multi-omics technologies, offers powerful new tools to address several key issues. These include precise delineation of species and evolutionary units, inference of past evolutionary histories and demographic fluctuations, identification of environmentally adaptive variants, and measurement of genetic load levels. As the urgency to deal with more extreme environmental stresses grows, understanding the genomics of evolutionary history, local adaptation, future responses to climate change, and conservation and restoration of natural forest trees will be critical for research at the nexus of global change, population genomics, and conservation biology. In this review, we explore the application of evolutionary genomics to assess the effects of global climate change using multi-omics approaches and discuss the outlook for breeding of climate-adapted trees.

Phylogenomics of mulberries (Morus, Moraceae) inferred from plastomes and single copy nuclear genes

Mulberries (genus Morus), belonging to the order Rosales, family Moraceae, are important woody plants due to their economic values in sericulture, as well as for nutritional benefits and medicinal values. However, the taxonomy and phylogeny of Morus, especially for the Asian species, remains challenging due to its wide geographical distribution, morphological plasticity, and interspecific hybridization. To better understand the evolutionary history of Morus, we combined plastomes and a large-scale nuclear gene analyses to investigate their phylogenetic relationships. We assembled the plastomes and screened 211 single-copy nuclear genes from 13 Morus species and related taxa. The plastomes of Morus species were relatively conserved in terms of genome size, gene content, synteny, IR boundary and codon usage. Using nuclear data, our results elucidated identical topologies based on coalescent and concatenation methods. The genus Morus was supported as monophyletic, with M. notabilis as the first diverging lineage and the two North American Morus species, M. microphylla and M. rubra, as sister to the other Asian species. In the Asian Morus species, interspecific relationships were completely resolved. However, cyto-nuclear discordances and gene tree-species tree conflicts were detected in the phylogenies of Morus, with multiple evidences supporting hybridization/introgression as the main cause of discordances between nuclear and plastid phylogenies, while gene tree-species tree conflicts were mainly caused by ILS.

A new species of Populus and the extensive hybrid speciation arising from it on the Qinghai-Tibet Plateau

While the diversity of species formation is broadly acknowledged, significant debate exists regarding the universal nature of hybrid species formation. Through an 18-year comprehensive study of all Populus species on the Qinghai-Tibet Plateau, 23 previously recorded species and 8 new species were identified. Based on morphological characteristics, these can be classified into three groups: species in section Leucoides, species with large leaves, and species with small leaves in section Tacamahaca. By conducting whole-genome re-sequencing of 150 genotypes from these 31 species, 2.28 million single nucleotide polymorphisms (SNPs) were identified. Phylogenetic analysis utilizing these SNPs not only revealed a highly intricate evolutionary network within the large-leaf species of section Tacamahaca but also confirmed that a new species, P. curviserrata, naturally hybridized with P. cathayana, P. szechuanica, and P. ciliata, resulting in 11 hybrid species. These findings indicate the widespread occurrence of hybrid species formation within this genus, with hybridization serving as a key evolutionary mechanism for Populus on the plateau. A novel hypothesis, "Hybrid Species Exterminating Their Ancestral Species (HSEAS)," is introduced to explain the mechanisms of hybrid species formation at three different scales: the entire plateau, the southeastern mountain region, and individual river valleys.

Detection of Ghost Introgression Requires Exploiting Topological and Branch Length Information

In recent years, the study of hybridization and introgression has made significant progress, with ghost introgression-the transfer of genetic material from extinct or unsampled lineages to extant species-emerging as a key area for research. Accurately identifying ghost introgression, however, presents a challenge. To address this issue, we focused on simple cases involving 3 species with a known phylogenetic tree. Using mathematical analyses and simulations, we evaluated the performance of popular phylogenetic methods, including HyDe and PhyloNet/MPL, and the full-likelihood method, Bayesian Phylogenetics and Phylogeography (BPP), in detecting ghost introgression. Our findings suggest that heuristic approaches relying on site-pattern counts or gene-tree topologies struggle to differentiate ghost introgression from introgression between sampled non-sister species, frequently leading to incorrect identification of donor and recipient species. The full-likelihood method BPP uses multilocus sequence alignments directly-hence taking into account both gene-tree topologies and branch lengths, by contrast, is capable of detecting ghost introgression in phylogenomic datasets. We analyzed a real-world phylogenomic dataset of 14 species of Jaltomata (Solanaceae) to showcase the potential of full-likelihood methods for accurate inference of introgression.

The super-pangenome of Populus unveils genomic facets for its adaptation and diversification in widespread forest trees

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.

Haplotype-resolved genome assembly of poplar line NL895 provides a valuable tree genomic resource

Poplar line NL895 can potentially become a model plant for poplar study as it is a widely cultivated elite line. However, the lack of genome resources hindered the use of NL895 as the major plant material in poplar. In this study, we provided a high-quality genome assembly for poplar line NL895 with PacBio single molecule real-time (SMRT) sequencing and High-throughput chromosome conformation capture (Hi-C) technology. The raw assembly of NL895 for the diploid genome included 606 contigs with a total size of ~815 Mb, and the monoploid genome included 246 contigs with a total size of ~412 Mb. The haplotype-resolved chromosomes in the diploid genomes were also generated. All the monoploid, diploid, and haplotype-resolved genomes showed more than 97% completeness and they can largely improve the mapping efficiency in RNA-Seq analysis. By comprehensively comparing the two haplotype genomes we found the heterozygosity of NL895 is much higher than other poplar lines. We also found that NL895 harbors more genomic variants and more gene diversity. The haplotype-specific genes showed higher variable gene expression patterns. These characters would be attributed to the high heterosis of poplar line NL895. The allele-specific expression (ASE) was also investigated and lots of alleles showed biased expressions in different tissues or environmental conditions. Taken together, the genome sequence for NL895 is a valuable tree genomic resource and it would greatly facilitate studies in poplar.

Balancing selection and the functional effects of shared polymorphism in cryptic Daphnia species

The patterns of genetic variation within and between related taxa represent the genetic history of a species. Shared polymorphisms, loci with identical alleles across species, are of unique interest as they may represent cases of ancient selection maintaining functional variation post-speciation. In this study, we investigate the abundance of shared polymorphism in the Daphnia pulex species complex. We test whether shared mutations are consistent with the action of balancing selection or alternative hypotheses such as hybridization, incomplete lineage sorting, or convergent evolution. We analyzed over 2,000 genomes from North American and European D. pulex and several outgroup species to examine the prevalence and distribution of shared alleles between the focal species pair, North American and European D. pulex. We show that while North American and European D. pulex diverged over ten million years ago, they retained tens of thousands of shared alleles. We found that the number of shared polymorphisms between North American and European D. pulex cannot be explained by hybridization or incomplete lineage sorting alone. Instead, we show that most shared polymorphisms could be the product of convergent evolution, that a limited number appear to be old trans-specific polymorphisms, and that balancing selection is affecting young and ancient mutations alike. Finally, we provide evidence that a blue wavelength opsin gene with trans-specific polymorphisms has functional effects on behavior and fitness in the wild. Ultimately, our findings provide insights into the genetic basis of adaptation and the maintenance of genetic diversity between species.

High-quality haplotype-resolved genome assembly for ring-cup oak (Quercus glauca) provides insight into oaks demographic dynamics

Quercus section Cyclobalanopsis represents a dominant woody lineage in East Asian evergreen broadleaved forests. Regardless of its ecological and economic importance, little is known about the genomes of species in this unique oak lineage. Quercus glauca is one of the most widespread tree species in the section Cyclobalanopsis. In this study, a high-quality haplotype-resolved reference genome was assembled for Q. glauca from PacBio HiFi and Hi-C reads. The genome size, contig N50, and scaffold N50 measured 902.88, 7.60, and 69.28 Mb, respectively, for haplotype1, and 913.28, 7.20, and 71.53 Mb, respectively, for haplotype2. A total of 37,457 and 38,311 protein-coding genes were predicted in haplotype1 and haplotype2, respectively. Homologous chromosomes in the Q. glauca genome had excellent gene pair collinearity. The number of R-genes in Q. glauca was similar to most East Asian oaks but less than oak species from Europe and America. Abundant structural variation in the Q. glauca genome could contribute to environmental stress tolerance in Q. glauca. Sections Cyclobalanopsis and Cerris diverged in the Oligocene, in agreement with fossil records for section Cyclobalanopsis, which document its presence in East Asia since the early Miocene. The demographic dynamics of closely related oak species were largely similar. The high-quality reference genome provided here for the most widespread species in section Cyclobalanopsis will serve as an essential genomic resource for evolutionary studies of key oak lineages while also supporting studies of interspecific introgression, local adaptation, and speciation in oaks.

The role of introgressive hybridization in shaping the geographically isolated gene pools of wax palm populations (genus Ceroxylon)

The speciation continuum is the process by which genetic groups diverge until they reach reproductive isolation. It has become common in the literature to show that this process is gradual and flickering, with possibly many instances of secondary contact and introgression after divergence has started. The level of divergence might vary among genomic regions due to, among others, the different forces and roles of selection played by the shared regions. Through hybrid capture, we sequenced ca. 4,000 nuclear regions in populations of six species of wax palms, five of which form a monophyletic group (genus Ceroxylon, Arecaceae: Ceroxyloideae). We show that in this group, the different populations show varying degrees of introgressive hybridization, and two of them are backcrosses of the other three 'pure' species. This is particularly interesting because these three species are dioecious, have a shared main pollinator, and have slightly overlapping reproductive seasons but highly divergent morphologies. Our work supports shows wax palms diverge under positive and background selection in allopatry, and hybridize due to secondary contact and inefficient reproductive barriers, which sustain genetic diversity. Introgressed regions are generally not under positive selection. Peripheral populations are backcrosses of other species; thus, introgressive hybridization is likely modulated by demographic effects rather than selective pressures. In general, these species might function as an 'evolutionary syngameon' where expanding, peripheral, small, and isolated populations maintain diversity by crossing with available individuals of other wax palms. In the Andean context, species can benefit from gained variation from a second taxon or the enhancement of population sizes by recreating a common genetic pool.

Genome-Wide Identification and Expression Profiling of Velvet Complex Transcription Factors in Populus alba × Populus glandulosa

The discovery of new genes with novel functions is a major driver of adaptive evolutionary innovation in plants. Especially in woody plants, due to genome expansion, new genes evolve to regulate the processes of growth and development. In this study, we characterized the unique VeA transcription factor family in Populus alba × Populus glandulosa, which is associated with secondary metabolism. Twenty VeA genes were characterized systematically on their phylogeny, genomic distribution, gene structure and conserved motif, promoter binding site, and expression profiling. Furthermore, through ChIP-qPCR, Y1H, and effector-reporter assays, it was demonstrated that PagMYB128 directly regulated PagVeA3 to influence the biosynthesis of secondary metabolites. These results provide a basis for further elucidating the function of VeAs gene in poplar and its genetic regulation mechanism.

Nuclear phylogenomics of angiosperms and insights into their relationships and evolution

Angiosperms (flowering plants) are by far the most diverse land plant group with over 300,000 species. The sudden appearance of diverse angiosperms in the fossil record was referred to by Darwin as the "abominable mystery," hence contributing to the heightened interest in angiosperm evolution. Angiosperms display wide ranges of morphological, physiological, and ecological characters, some of which have probably influenced their species richness. The evolutionary analyses of these characteristics help to address questions of angiosperm diversification and require well resolved phylogeny. Following the great successes of phylogenetic analyses using plastid sequences, dozens to thousands of nuclear genes from next-generation sequencing have been used in angiosperm phylogenomic analyses, providing well resolved phylogenies and new insights into the evolution of angiosperms. In this review we focus on recent nuclear phylogenomic analyses of large angiosperm clades, orders, families, and subdivisions of some families and provide a summarized Nuclear Phylogenetic Tree of Angiosperm Families. The newly established nuclear phylogenetic relationships are highlighted and compared with previous phylogenetic results. The sequenced genomes of Amborella, Nymphaea, Chloranthus, Ceratophyllum, and species of monocots, Magnoliids, and basal eudicots, have facilitated the phylogenomics of relationships among five major angiosperms clades. All but one of the 64 angiosperm orders were included in nuclear phylogenomics with well resolved relationships except the placements of several orders. Most families have been included with robust and highly supported placements, especially for relationships within several large and important orders and families. Additionally, we examine the divergence time estimation and biogeographic analyses of angiosperm on the basis of the nuclear phylogenomic frameworks and discuss the differences compared with previous analyses. Furthermore, we discuss the implications of nuclear phylogenomic analyses on ancestral reconstruction of morphological, physiological, and ecological characters of angiosperm groups, limitations of current nuclear phylogenomic studies, and the taxa that require future attention.

Genomic Insights into High-Altitude Adaptation: A Comparative Analysis of Roscoea alpina and R. purpurea in the Himalayas

Environmental stress at high altitudes drives the development of distinct adaptive mechanisms in plants. However, studies exploring the genetic adaptive mechanisms of high-altitude plant species are scarce. In the present study, we explored the high-altitude adaptive mechanisms of plants in the Himalayas through whole-genome resequencing. We studied two widespread members of the Himalayan endemic alpine genus Roscoea (Zingiberaceae): R. alpina (a selfing species) and R. purpurea (an outcrossing species). These species are distributed widely in the Himalayas with distinct non-overlapping altitude distributions; R. alpina is distributed at higher elevations, and R. purpurea occurs at lower elevations. Compared to R. purpurea, R. alpina exhibited higher levels of linkage disequilibrium, Tajima's D, and inbreeding coefficient, as well as lower recombination rates and genetic diversity. Approximately 96.3% of the genes in the reference genome underwent significant genetic divergence (FST ≥ 0.25). We reported 58 completely divergent genes (FST = 1), of which only 17 genes were annotated with specific functions. The functions of these genes were primarily related to adapting to the specific characteristics of high-altitude environments. Our findings provide novel insights into how evolutionary innovations promote the adaptation of mountain alpine species to high altitudes and harsh habitats.

Phylogenomics reveal Populusgonggaensis as a hybrid between P.lasiocarpa and P.cathayana (Salicaceae)

High levels of intra-specific polymorphism and frequent hybridisation make it difficult to define species and correctly apply their scientific names. Populus L. is a challenging genus with plentiful natural and artificial hybrids. This study is a part of the project 'Flora of Pan-Himalaya' and aims to determine the taxonomic identity of P.gonggaensis N. Chao & J.R. He and to find out whether it is of hybrid origin. Whole-genome sequencing data were obtained from 57 samples. The SNP matrix was developed for phylogenetic reconstruction, ABBA-BABA statistics, PCA and ADMIXTURE analysis. The results indicate that P.gonggaensis is a spontaneous hybrid between P.lasiocarpa and P.cathayana. This study points out the importance of SNP data and comprehensive analyses for discovering the potential interspecific hybridisation and clarifies the usage of the name. In addition, the lectotype of P.gonggaensis was designated.

Editing Metabolism, Sex, and Microbiome: How Can We Help Poplar Resist Pathogens?

Poplar (Populus) is a genus of woody plants of great economic value. Due to the growing economic importance of poplar, there is a need to ensure its stable growth by increasing its resistance to pathogens. Genetic engineering can create organisms with improved traits faster than traditional methods, and with the development of CRISPR/Cas-based genome editing systems, scientists have a new highly effective tool for creating valuable genotypes. In this review, we summarize the latest research data on poplar diseases, the biology of their pathogens and how these plants resist pathogens. In the final section, we propose to plant male or mixed poplar populations; consider the genes of the MLO group, transcription factors of the WRKY and MYB families and defensive proteins BbChit1, LJAMP2, MsrA2 and PtDef as the most promising targets for genetic engineering; and also pay attention to the possibility of microbiome engineering.

A model of hybrid speciation process drawn from three new poplar species originating from distant hybridization between sections

Although numerous studies have been conducted on hybrid speciation, our understanding of this process remains limited. Through an 18-year systematic investigation of all taxa of Populus on the Qinghai-Tibet Plateau, we discovered three new taxa with clear characteristics of sect. Leucoides. Further evidence was gathered from morphology, whole-genome bioinformatics, biogeography, and breeding to demonstrate synthetically that they all originated from distant hybridization between sect. Leucoides and sect. Tacamahaca. P. gonggaensis originated from the hybridization of P. lasiocarpa with P. cathayana, P. butuoensis from the hybridization of P. wilsonii with P. szechuanica, and P. dafengensis from the hybridization of P. lasiocarpa with P. szechuanica. Due to heterosis, the three hybrid taxa possess greater ecological adaptability than their ancestral species. We propose a hybrid speciation process model that incorporates orthogonal, reverse, and backcrossing events. This model can adequately explain some crucial evolutionary concerns, such as the nuclear-cytoplasmic conflict on phylogeny and the extinction of ancestral species within the distribution range of hybrid species.

Populus cathayana genome and population resequencing provide insights into its evolution and adaptation

Populus cathayana Rehder, an indigenous poplar species of ecological and economic importance, is widely distributed in a high-elevation range from southwest to northeast China. Further development of this species as a sustainable poplar resource has been hindered by a lack of genome information the at the population level. Here, we produced a chromosome-level genome assembly of P. cathayana, covering 406.55 Mb (scaffold N50 = 20.86 Mb) and consisting of 19 chromosomes, with 35 977 protein-coding genes. Subsequently, we made a genomic variation atlas of 438 wild individuals covering 36 representative geographic areas of P. cathayana, which were divided into four geographic groups. It was inferred that the Northwest China regions served as the genetic diversity centers and a population bottleneck happened during the history of P. cathayana. By genotype-environment association analysis, 947 environment-association loci were significantly associated with temperature, solar radiation, precipitation, and altitude variables. We identified local adaptation genes involved in DNA repair and UV radiation response, among which UVR8, HY5, and CUL4 had key roles in high-altitude adaptation of P. cathayana. Predictions of adaptive potential under future climate conditions showed that P. cathayana populations in areas with drastic climate change were anticipated to have greater maladaptation risk. These results provide comprehensive insights for understanding wild poplar evolution and optimizing adaptive potential in molecular breeding.

Phylogenomics reveals patterns of ancient hybridization and differential diversification that contribute to phylogenetic conflict in willows, poplars, and close relatives

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.

Identification and adaptive evolution analysis of glutaredoxin genes in Populus spp

Glutaredoxin (GRX) is a class of small redox proteins widely involved in cellular redox homeostasis and the regulation of various cellular processes. The role of GRX gene in the differentiation of Populus spp. is rarely reported. We compared the similarities and differences of GRX genes among four sections of poplar using bioinformatics, corrected the annotations of some GRX genes, and focused on analysing their transcript profiling and adaptive evolution in Populus spp. A total of 219 GRX genes were identified in four sections of poplar, among which annotations for 13 genes were corrected. Differences in GRX genes were found between sect. Turanga, represented by P. euphratica, and other poplar sections. Most notably, P. euphratica had the smallest number of duplication events for GRX genes (n = 9) and no tandem duplications, whereas there were >25 duplication events for all other poplars. Furthermore, we detected 18 pairs of GRX genes under positive selection pressure in various sections of poplar, and identified two groups of GRX genes in the Salicaceae that potentially underwent positive selection. Expression profiling results showed that the PtrGRX34 and its orthologous genes were upregulated under stress treatments. In summary, the GRX gene family underwent expansion during poplar differentiation, and some genes underwent rapid evolution during this process, which may be beneficial for Populus spp. to adapt to environmental changes. This study may provide more insights into the molecular mechanisms of Populus spp. adaptation to environmental changes and the adaptive evolution of GRX genes.

Haplotype-resolved genome assembly of Populus tremula × P. alba reveals aspen-specific megabase satellite DNA

Populus species play a foundational role in diverse ecosystems and are important renewable feedstocks for bioenergy and bioproducts. Hybrid aspen Populus tremula × P. alba INRA 717-1B4 is a widely used transformation model in tree functional genomics and biotechnology research. As an outcrossing interspecific hybrid, its genome is riddled with sequence polymorphisms which present a challenge for sequence-sensitive analyses. Here we report a telomere-to-telomere genome for this hybrid aspen with two chromosome-scale, haplotype-resolved assemblies. We performed a comprehensive analysis of the repetitive landscape and identified both tandem repeat array-based and array-less centromeres. Unexpectedly, the most abundant satellite repeats in both haplotypes lie outside of the centromeres, consist of a 147 bp monomer PtaM147, frequently span >1 megabases, and form heterochromatic knobs. PtaM147 repeats are detected exclusively in aspens (section Populus) but PtaM147-like sequences occur in LTR-retrotransposons of closely related species, suggesting their origin from the retrotransposons. The genomic resource generated for this transformation model genotype has greatly improved the design and analysis of genome editing experiments that are highly sensitive to sequence polymorphisms. The work should motivate future hypothesis-driven research to probe into the function of the abundant and aspen-specific PtaM147 satellite DNA.

Sex-determining region complements traditionally used in phylogenetic studies nuclear and chloroplast sequences in investigation of Aigeiros Duby and Tacamahaca Spach poplars (genus Populus L., Salicaceae)

Members of the genus Populus L. play an important role in the formation of forests in the northern hemisphere and are used in urban landscaping and timber production. Populus species of closely related sections show extensive hybridization. Therefore, the systematics of the genus is rather complicated, especially for poplars of hybrid origin. We aimed to assess the efficiency of application of the sex-determining region (SDR) in addition to the nuclear and chloroplast genome loci traditionally used in phylogenetic studies of poplars to investigate relationships in sections Aigeiros Duby and Tacamahaca Spach. Targeted deep sequencing of NTS 5S rDNA, ITS, DSH 2, DSH 5, DSH 8, DSH 12, DSH 29, 6, 15, 16, X18, trnG-psbK-psbI, rps2-rpoC2, rpoC2-rpoC1, as well as SDR and ARR17 gene was performed for 379 poplars. The SDR and ARR17 gene together with traditionally used multicopy and single-copy loci of nuclear and chloroplast DNA allowed us to obtain a clustering that is most consistent with poplar systematics based on morphological data and to shed light on several controversial hypotheses about the origin of the studied taxa (for example, the inexpediency of separating P. koreana, P. maximowiczii, and P. suaveolens into different species). We present a scheme of relationships between species and hybrids of sections Aigeiros and Tacamahaca based on molecular genetic, morphological, and geographical data. The geographical proximity of species and, therefore, the possibility of hybridization between them appear to be more important than the affiliation of species to the same section. We speculate that sections Aigeiros and Tacamahaca are distinguished primarily on an ecological principle (plain and mountain poplars) rather than on a genetic basis. Joint analysis of sequencing data for the SDR and chloroplast genome loci allowed us to determine the ancestors of P. × petrovskoe - P. laurifolia (female tree) × P. × canadensis (male tree), and P. × rasumovskoe - P. nigra (female tree) × P. suaveolens (male tree). Thus, the efficiency of using the SDR for the study of poplars of sections Aigeiros and Tacamahaca and the prospects of its use for the investigation of species of the genus Populus were shown.

Phylogenomics reveals extreme gene tree discordance in a lineage of dominant trees: hybridization, introgression, and incomplete lineage sorting blur deep evolutionary relationships despite clear species groupings in Eucalyptus subgenus Eudesmia

Eucalypts are a large and ecologically important group of plants on the Australian continent, and understanding their evolution is important in understanding evolution of the unique Australian flora. Previous phylogenies using plastome DNA, nuclear-ribosomal DNA, or random genome-wide SNPs, have been confounded by limited genetic sampling or by idiosyncratic biological features of the eucalypts, including widespread plastome introgression. Here we present phylogenetic analyses of Eucalyptus subgenus Eudesmia (22 species from western, northern, central and eastern Australia), in the first study to apply a target-capture sequencing approach using custom, eucalypt-specific baits (of 568 genes) to a lineage of Eucalyptus. Multiple accessions of all species were included, and target-capture data were supplemented by separate analyses of plastome genes (average of 63 genes per sample). Analyses revealed a complex evolutionary history likely shaped by incomplete lineage sorting and hybridization. Gene tree discordance generally increased with phylogenetic depth. Species, or groups of species, toward the tips of the tree are mostly supported, and three major clades are identified, but the branching order of these clades cannot be confirmed with confidence. Multiple approaches to filtering the nuclear dataset, by removing genes or samples, could not reduce gene tree conflict or resolve these relationships. Despite inherent complexities in eucalypt evolution, the custom bait kit devised for this research will be a powerful tool for investigating the evolutionary history of eucalypts more broadly.

Nutrient conditions mediate mycorrhizal effects on biomass production and cell wall chemistry in poplar

Large-scale biofuel production from lignocellulosic feedstock is limited by the financial and environmental costs associated with growing and processing lignocellulosic material and the resilience of these plants to environmental stress. Symbiotic associations with arbuscular (AM) and ectomycorrhizal (EM) fungi represent a potential strategy for expanding feedstock production while reducing nutrient inputs. Comparing AM and EM effects on wood production and chemical composition is a necessary step in developing biofuel feedstocks. Here, we assessed the productivity, biomass allocation and secondary cell wall (SCW) composition of greenhouse-grown Populus tremuloidesMichx. inoculated with either AM or EM fungi. Given the long-term goal of reducing nutrient inputs for biofuel production, we further tested the effects of nutrient availability and nitrogen:phosphorus stoichiometry on mycorrhizal responses. Associations with both AM and EM fungi increased plant biomass by 14-74% depending on the nutrient conditions but had minimal effects on SCW composition. Mycorrhizal plants, especially those inoculated with EM fungi, also allocated a greater portion of their biomass to roots, which could be beneficial in the field where plants are likely to experience both water and nutrient stress. Leaf nutrient content was weakly but positively correlated with wood production in mycorrhizal plants. Surprisingly, phosphorus played a larger role in EM plants compared with AM plants. Relative nitrogen and phosphorus availability were correlated with shifts in SCW composition. For AM associations, the benefit of increased wood biomass may be partially offset by increased lignin content, a trait that affects downstream processing of lignocellulosic tissue for biofuels. By comparing AM and EM effects on the productivity and chemical composition of lignocellulosic tissue, this work links broad functional diversity in mycorrhizal associations to key biofuel traits and highlights the importance of considering both biotic and abiotic factors when developing strategies for sustainable biofuel production.

Drivers of genomic landscapes of differentiation across a Populus divergence gradient

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.

Wind-dispersed seeds blur phylogeographic breaks: The complex evolutionary history of Populus lasiocarpa around the Sichuan Basin

The strength of phylogeographic breaks can vary among species in the same area despite being subject to the same geological and climate history due to differences in biological traits. Several important phylogeographic breaks exist around the Sichuan Basin in Southwest China but few studies have focused on wind-dispersed plants. Here, we investigated the phylogeographic patterns and the evolutionary history of Populus lasiocarpa, a wind-pollinated and wind-dispersed tree species with a circum-Sichuan Basin distribution in southwest China. We sequenced and analyzed three plastid DNA fragments (ptDNA) and eight nuclear microsatellites (nSSRs) of 265 individuals of P. lasiocarpa from 21 populations spanning the entire distribution range. Distribution patterns based on nSSR data revealed that there are three genetic groups in P. lasiocarpa. This is consistent with the three phylogeographic breaks (Sichuan Basin, the Kaiyong Line and the 105°E line), where the Sichuan basin acts as the main barrier to gene flow between western and eastern groups. However, the distribution pattern based on ptDNA haplotypes poorly matched the phylogeographic breaks, and wind-dispersed seeds may be one of the main contributing factors. Species distribution modelling suggested a larger potential distribution in the last glacial maximum with a severe bottleneck during the last interglacial. A DIYABC model also suggested a population contraction and expansion for both western and eastern lineages. These results indicate that biological traits are likely to affect the evolutionary history of plants, and that nuclear molecular markers, which experience higher levels of gene flow, might be better indicators of phylogeographic breaks.

Genomic divergence and introgression among three Populus species

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.

A chromosome-level Populus qiongdaoensis genome assembly provides insights into tropical adaptation and a cryptic turnover of sex determination

Populus species have long been used as model organisms to study the adaptability of trees and the evolution of sex chromosomes. As a species belonging to the section Populus and limited to tropical areas, the P. qiongdaoensis genome contains important information for tropical poplar studies and protection. Here, we report a chromosome-level genome assembly and annotation of a female P. qiongdaoensis. Gene family clustering, positive selection detection and historical reconstruction of population dynamics revealed the tropical adaptation of P. qiongdaoensis, and showed convergent evolution with another tropical poplar, P. ilicifolia, at the molecular level, especially on some functional genes (e.g., PIF3 and PIL1). In addition, we also identified a ZW sex determination system on chromosome 19 of P. qiongdaoensis, and inferred that it seems to have a similar sex determination mechanism to other poplars, controlled by a type-A cytokinin response regulator (RR) gene. However, comparison and phylogenetic analysis of the sex determination regions confirmed a cryptic sex turnover event in the section Populus, which may be caused by the translocation and duplication of the RR gene driven by Helitron-like transposable elements. Our study provides new insights into the environmental adaptation and sex chromosome evolution of poplars, and emphasizes the importance of using long read sequencing in ecological and evolutionary inferences of plants.

High-resolution mapping reveals hotspots and sex-biased recombination in Populus trichocarpa

Fine-scale meiotic recombination is fundamental to the outcome of natural and artificial selection. Here, dense genetic mapping and haplotype reconstruction were used to estimate recombination for a full factorial Populus trichocarpa cross of 7 males and 7 females. Genomes of the resulting 49 full-sib families (N = 829 offspring) were resequenced, and high-fidelity biallelic SNP/INDELs and pedigree information were used to ascertain allelic phase and impute progeny genotypes to recover gametic haplotypes. The 14 parental genetic maps contained 1,820 SNP/INDELs on average that covered 376.7 Mb of physical length across 19 chromosomes. Comparison of parental and progeny haplotypes allowed fine-scale demarcation of cross-over regions, where 38,846 cross-over events in 1,658 gametes were observed. Cross-over events were positively associated with gene density and negatively associated with GC content and long-terminal repeats. One of the most striking findings was higher rates of cross-overs in males in 8 out of 19 chromosomes. Regions with elevated male cross-over rates had lower gene density and GC content than windows showing no sex bias. High-resolution analysis identified 67 candidate cross-over hotspots spread throughout the genome. DNA sequence motifs enriched in these regions showed striking similarity to those of maize, Arabidopsis, and wheat. These findings, and recombination estimates, will be useful for ongoing efforts to accelerate domestication of this and other biomass feedstocks, as well as future studies investigating broader questions related to evolutionary history, perennial development, phenology, wood formation, vegetative propagation, and dioecy that cannot be studied using annual plant model systems.

Origin and evolutionary history of Populus (Salicaceae): Further insights based on time divergence and biogeographic analysis

INTRODUCTION

Populus (Salicaceae) species harbour rich biodiversity and are widely distributed throughout the Northern Hemisphere. However, the origin and biogeography of Populus remain poorly understood.

METHODS

We infer the divergence times and the historical biogeography of the genus Populus through phylogenetic analysis of 34 chloroplast fragments based on a large sample.

RESULTS AND DISCUSSION

Eurasia is the likely location of the early divergences of Salicaceae after the Cretaceous-Paleogene (K-Pg) mass extinction, followed by recurrent spread to the remainder of the Old World and the New World beginning in the Eocene; the extant Populus species began to diversity during the early Oligocene (approximately 27.24 Ma), climate changes during the Oligocene may have facilitated the diversification of modern poplar species; three separate lineages of Populus from Eurasia colonized North America in the Cenozoic via the Bering Land Bridges (BLB); We hypothesize that the present day disjunction in Populus can be explained by two scenarios: (i) Populus likely originated in Eurasia and subsequently colonized other regions, including North America; and (ii) the fact that the ancestor of the genus Populus that was once widely distributed in the Northern Hemisphere and eventually wiped out due to the higher extinction rates in North America, similar to the African Rand flora. We hypothesize that disparities in extinction across the evolutionary history of Populus in different regions shape the modern biogeography of Populus. Further studies with dense sampling and more evidence are required to test these hypotheses. Our research underscores the significance of combining phylogenetic analyses with biogeographic interpretations to enhance our knowledge of the origin, divergence, and distribution of biodiversity in temperate plant floras.

Pervasive hybridization during evolutionary radiation of Rhododendron subgenus Hymenanthes in mountains of southwest China

Radiations are especially important for generating species biodiversity in mountainous ecosystems. The contribution of hybridization to such radiations has rarely been examined. Here, we use extensive genomic data to test whether hybridization was involved in evolutionary radiation within Rhododendron subgenus Hymenanthes, whose members show strong geographic isolation in the mountains of southwest China. We sequenced genomes for 143 species of this subgenus and 93 species of four other subgenera, and found that Hymenanthes was monophyletic and radiated during the late Oligocene to middle Miocene. Widespread hybridization events were inferred within and between the identified clades and subclades. This suggests that hybridization occurred both early and late during diversification of subgenus Hymenanthes, although the extent to which hybridization, speciation through mixing-isolation-mixing or hybrid speciation, accelerated the diversification needs further exploration. Cycles of isolation and contact in such and other montane ecosystems may have together promoted species radiation through hybridization between diverging populations and species. Similar radiation processes may apply to other montane floras in this region and elsewhere.

Repeated turnovers keep sex chromosomes young in willows

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02769-w.

Duplication of NRAMP3 gene in poplars generated two homologous transporters with distinct functions

Transition metals are essential for a wealth of metabolic reactions, but their concentrations need to be tightly controlled across cells and cell compartments, as metal excess or imbalance has deleterious effects. Metal homeostasis is achieved by a combination of metal transport across membranes and metal binding to a variety of molecules. Gene duplication is a key process in evolution, as emergence of advantageous mutations on one of the copies can confer a new function. Here, we report that the poplar genome contains two paralogues encoding NRAMP3 metal transporters localized in tandem. All Populus species analyzed had two copies of NRAMP3, whereas only one could be identified in Salix species indicating that duplication occurred when the two genera separated. Both copies are under purifying selection and encode functional transporters, as shown by expression in the yeast heterologous expression system. However, genetic complementation revealed that only one of the paralogues has retained the original function in release of metals stored in the vacuole previously characterized in A. thaliana. Confocal imaging showed that the other copy has acquired a distinct localization to the Trans Golgi Network (TGN). Expression in poplar suggested that the copy of NRAMP3 localized on the TGN has a novel function in the control of cell-to-cell transport of manganese. This work provides a clear case of neo-functionalization through change in the subcellular localization of a metal transporter as well as evidence for the involvement of the secretory pathway in cell-to-cell transport of manganese.

Species divergence with gene flow and hybrid speciation on the Qinghai-Tibet Plateau

The Qinghai-Tibet Plateau (QTP) sensu lato (sl), comprising the platform, the Himalaya and the Hengduan Mountains, is characterized by a large number of endemic plant species. This evolutionary cradle may have arisen from explosive species diversification because of geographic isolation. However, gene flow has been widely detected during the speciation processes of all groups examined, suggesting that natural selection may have also played an important role during species divergence in this region. In addition, natural hybrids have been recovered in almost all species-rich genera. This suggests that numerous species in this region are still 'on the speciation pathway to complete reproductive isolation (RI)'. Such hybrids could directly develop into new species through hybrid polyploidization and homoploid hybrid speciation (HHS). HHS may take place more easily than previously thought through alternate inheritance of alleles of parents at multiple RI loci. Therefore, isolation, selection and hybridization could together have promoted species diversification of numerous plant genera on the QTP sl. We emphasize the need for identification and functional analysis of alleles of major genes for speciation, and especially encourage investigations of parallel adaptive divergence causing RI across different lineages within similar but specific habitats in this region.

Origin and diversification of the saguaro cactus (Carnegiea gigantea): a within-species phylogenomic analysis

Reconstructing accurate historical relationships within a species poses numerous challenges, not least in many plant groups in which gene flow is high enough to extend well beyond species boundaries. Nonetheless, the extent of tree-like history within a species is an empirical question on which it is now possible to bring large amounts of genome sequence to bear. We assess phylogenetic structure across the geographic range of the saguaro cactus, an emblematic member of Cactaceae, a clade known for extensive hybridization and porous species boundaries. Using 200 Gb of whole genome resequencing data from 20 individuals sampled from 10 localities, we assembled two data sets comprising 150,000 biallelic single nucleotide polymorphisms (SNPs) from protein coding sequences. From these we inferred within-species trees and evaluated their significance and robustness using five qualitatively different inference methods. Despite the low sequence diversity, large census population sizes, and presence of wide-ranging pollen and seed dispersal agents, phylogenetic trees were well resolved and highly consistent across both data sets and all methods. We inferred that the most likely root, based on marginal likelihood comparisons, is to the east and south of the region of highest genetic diversity, which lies along the coast of the Gulf of California in Sonora, Mexico. Together with striking decreases in marginal likelihood found to the north, this supports hypotheses that saguaro's current range reflects post-glacial expansion from the refugia in the south of its range. We conclude with observations about practical and theoretical issues raised by phylogenomic data sets within species, in which SNP-based methods must be used rather than gene tree methods that are widely used when sequence divergence is higher. These include computational scalability, inference of gene flow, and proper assessment of statistical support in the presence of linkage effects.

Phylogenomics and Biogeography of Populus Based on Comprehensive Sampling Reveal Deep-Level Relationships and Multiple Intercontinental Dispersals

Populus not only has significant economic and ecological values, but also serves as a model tree that is widely used in the basic research of tree growth, physiology, and genetics. However, high levels of morphological variation and extensive interspecific hybridization of Populus pose an obstacle for taxonomy, and also to the understanding of phylogenetic interspecific relationships and biogeographical history. In this study, a total of 103 accessions representing almost all the wild species of Populus were collected and whole-genome re-sequenced to examine the phylogenetic relationships and biogeography history. On the basis of 12,916,788 nuclear single nucleotide polymorphisms (SNPs), we reconstructed backbone phylogenies using concatenate and coalescent methods, we highly disentangled the species relationships of Populus, and several problematic taxa were treated as species complexes. Furthermore, the phylogeny of the chloroplast genome showed extensive discordance with the trees from the nuclear genome data, and due to extensive chloroplast capture and hybridization of Populus species, plastomes could not accurately evaluate interspecies relationships. Ancient gene flow between clades and some hybridization events were also identified by ABBA-BABA analysis. The reconstruction of chronogram and ancestral distributions suggested that North America was the original region of this genus, and subsequent long dispersal and migration across land bridges were contributed to the modern range of Populus. The diversification of Populus mainly occurred in East Asia in recent 15 Ma, possibly promoted by the uplift of the Tibetan Plateau. This study provided comprehensive evidence on the phylogeny of Populus and proposed a four-subgeneric classification and a new status, subgenus Abaso. Meanwhile, ancestral distribution reconstruction with nuclear data advanced the understanding of the biogeographic history of Populus.

Demographic history and natural selection shape patterns of deleterious mutation load and barriers to introgression across Populus genome

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 F₁ 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.

Trans-specific polymorphism and the convergent evolution of supertypes in major histocompatibility complex class II genes in darters (Etheostoma)

Major Histocompatibility Complex (MHC) genes are one of the most polymorphic gene groups known in vertebrates. MHC genes also exhibit allelic variants that are shared among taxa, referred to as trans-specific polymorphism (TSP). The role that selection plays in maintaining such high diversity within species, as well as TSP, is an ongoing discussion in biology. In this study, we used deep-sequencing techniques to characterize MHC class IIb gene diversity in three sympatric species of darters. We found at least 5 copies of the MHC gene in darters, with 126 genetic variants encoding 122 unique amino acid sequences. We identified four supertypes based on the binding properties of proteins encoded by the sequences. Although each species had a unique pool of variants, many variants were shared between species pairs and across all three species. Phylogenetic analysis showed that the variants did not group together monophyletically based on species identity or on supertype. An expanded phylogenetic analysis showed that some darter alleles grouped together with alleles from other percid fishes. Our findings show that TSP occurs in darters, which suggests that balancing selection is acting at the genotype level. Supertypes, however, are most likely evolving convergently, as evidenced by the fact that alleles do not form monophyletic groups based on supertype. Our research demonstrates that selection may be acting differently on MHC genes at the genotype and supertype levels, selecting for the maintenance of high genotypic diversity while driving the convergent evolution of similar MHC phenotypes across different species.

When tropical and subtropical congeners met: Multiple ancient hybridization events within Eriobotrya in the Yunnan-Guizhou Plateau, a tropical-subtropical transition area in China

Global climate changes during the Miocene may have created ample opportunities for hybridization between members of tropical and subtropical biomes at the boundary between these zones. Yet, very few studies have explored this possibility. The Yunnan-Guizhou Plateau (YGP) in Southwest China is a biodiversity hotspot for vascular plants, located in a transitional area between the floristic regions of tropical Southeast Asia and subtropical East Asia. The genus Eriobotrya (Rosaceae) comprises both tropical and subtropical taxa, with 12 species recorded in the YGP, making it a suitable basis for testing the hypothesis of between-biome hybridization. Therefore, we surveyed the evolutionary history of Eriobotrya by examining three chloroplast regions and five nuclear genes for 817 individuals (47 populations) of 23 Eriobotrya species (including 19 populations of 12 species in the YGP), plus genome re-sequencing of 33 representative samples. We concluded that: (1) phylogenetic positions for 16 species exhibited strong cytonuclear conflicts, most probably due to ancient hybridization; (2) the YGP is a hotspot for hybridization, with 11 species showing clear evidence of chloroplast capture; and (3) Eriobotrya probably originated in tropical Asia during the Eocene. From the Miocene onwards, the intensification of the Eastern Asia monsoon and global cooling may have shifted the tropical-subtropical boundary and caused secondary contact between species, thus providing ample opportunity for hybridization and diversification of Eriobotrya, especially in the YGP. Our study highlights the significant role that paleoclimate changes probably played in driving hybridization and generating rich species diversity in climate transition zones.

Candidates for Balancing Selection in Leishmania donovani Complex Parasites

The Leishmania donovani species complex is the causative agent of visceral leishmaniasis, which cause 20–40,000 fatalities a year. Here, we conduct a screen for balancing selection in this species complex. We used 384 publicly available L. donovani and L. infantum genomes, and sequence 93 isolates of L. infantum from Brazil to describe the global diversity of this species complex. We identify five genetically distinct populations that are sufficiently represented by genomic data to search for signatures of selection. We find that signals of balancing selection are generally not shared between populations, consistent with transient adaptive events, rather than long-term balancing selection. We then apply multiple diversity metrics to identify candidate genes with robust signatures of balancing selection, identifying a curated set of 24 genes with robust signatures. These include zeta toxin, nodulin-like, and flagellum attachment proteins. This study highlights the extent of genetic divergence between L. donovani complex parasites and provides genes for further study.

Resequencing of 388 cassava accessions identifies valuable loci and selection for variation in heterozygosity

Background

Heterozygous genomes are widespread in outcrossing and clonally propagated crops. However, the variation in heterozygosity underlying key agronomic traits and crop domestication remains largely unknown. Cassava is a staple crop in Africa and other tropical regions and has a highly heterozygous genome.

Results

We describe a genomic variation map from 388 resequenced genomes of cassava cultivars and wild accessions. We identify 52 loci for 23 agronomic traits through a genome-wide association study. Eighteen allelic variations in heterozygosity for nine candidate genes are significantly associated with seven key agronomic traits. We detect 81 selective sweeps with decreasing heterozygosity and nucleotide diversity, harboring 548 genes, which are enriched in multiple biological processes including growth, development, hormone metabolisms and responses, and immune-related processes. Artificial selection for decreased heterozygosity has contributed to the domestication of the large starchy storage root of cassava. Selection for homozygous GG allele in MeTIR1 during domestication contributes to increased starch content. Selection of homozygous AA allele in MeAHL17 is associated with increased storage root weight and cassava bacterial blight (CBB) susceptibility. We have verified the positive roles of MeTIR1 in increasing starch content and MeAHL17 in resistance to CBB by transient overexpression and silencing analysis. The allelic combinations in MeTIR1 and MeAHL17 may result in high starch content and resistance to CBB.

Conclusions

This study provides insights into allelic variation in heterozygosity associated with key agronomic traits and cassava domestication. It also offers valuable resources for the improvement of cassava and other highly heterozygous crops.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02524-7.

Evolutionary origin of a tetraploid Allium species on the Qinghai-Tibet Plateau

Extinct taxa may be detectable if they were ancestors to extant hybrid species, which retain their genetic signature. In this study, we combined phylogenomics, population genetics and fluorescence in situ hybridization (GISH and FISH) analyses to trace the origin of the alpine tetraploid Allium tetraploideum (2n = 4x = 32), one of the five known members in the subgenus Cyathophora. We found that A. tetraploideum was an obvious allotetrapoploid derived from ancestors including at least two closely related diploid species, A. farreri and A. cyathophorum, from which it differs by multiple ecological and genomic attributes. However, these two species cannot account for the full genome of A. tetraploideum, indicating that at least one extinct diploid is also involved in its ancestry. Furthermore, A. tetraploideum appears to have arisen via homoploid hybrid speciation (HHS) from two extinct allotetraploid parents, which derived in turn from the aforementioned diploids. Other modes of origin were possible, but all were even more complex and involved additional extinct ancestors. Our study together highlights how some polyploid species might have very complex origins, involving both HHS and polyploid speciation and also extinct ancestors.

Genome-wide comparative analyses of GATA transcription factors among seven Populus genomes

GATA transcription factors (TFs) are widespread eukaryotic regulators whose DNA-binding domain is a class IV zinc finger motif (CX2CX17-20CX2C) followed by a basic region. We identified 262 GATA genes (389 GATA TFs) from seven Populus genomes using the pipeline of GATA-TFDB. Alternative splicing forms of Populus GATA genes exhibit dynamics of GATA gene structures including partial or full loss of GATA domain and additional domains. Subfamily III of Populus GATA genes display lack CCT and/or TIFY domains. 21 Populus GATA gene clusters (PCs) were defined in the phylogenetic tree of GATA domains, suggesting the possibility of subfunctionalization and neofunctionalization. Expression analysis of Populus GATA genes identified the five PCs displaying tissue-specific expression, providing the clues of their biological functions. Amino acid patterns of Populus GATA motifs display well conserved manner of Populus GATA genes. The five Populus GATA genes were predicted as membrane-bound GATA TFs. Biased chromosomal distributions of GATA genes of three Populus species. Our comparative analysis approaches of the Populus GATA genes will be a cornerstone to understand various plant TF characteristics including evolutionary insights.

Unveiling the Patterns of Reticulated Evolutionary Processes with Phylogenomics: Hybridization and Polyploidy in the genus Rosa

Reticulation, caused by hybridization and allopolyploidization, is considered an important and frequent phenomenon in the evolution of numerous plant lineages. Although both processes represent important driving forces of evolution, they are mostly ignored in phylogenetic studies involving a large number of species. Indeed only a scattering of methods exists to recover a comprehensive reticulated evolutionary history for a broad taxon sampling. Among these methods, comparisons of topologies obtained from plastid markers with those from a few nuclear sequences are favored, even though they restrict in-depth studies of hybridization and polyploidization. The genus Rosa encompasses c. 150 species widely distributed throughout the northern hemisphere and represents a challenging taxonomic group in which hybridization and polyploidization are prominent. Our main objective was to develop a general framework that would take patterns of reticulation into account in the study of the phylogenetic relationships among Rosa species. Using amplicon sequencing we targeted allele variation in the nuclear genome as well as haploid sequences in the chloroplast genome. We successfully recovered robust plastid and nuclear phylogenies and performed in-depth tests for several scenarios of hybridization using a maximum pseudo-likelihood approach on taxon subsets. Our diploid-first approach followed by hybrid and polyploid grafting resolved most of the evolutionary relationships among Rosa subgenera, sections, and selected species. Based on these results, we provide new directions for a future revision of the infrageneric classification in Rosa. The stepwise strategy proposed here can be used to reconstruct the phylogenetic relationships of other challenging taxonomic groups with large numbers of hybrid and polyploid taxa.

Species delimitation and hybridization history of a hazel species complex

BACKGROUND AND AIMS

Hybridization increases species adaptation and biodiversity but also obscures species boundaries. In this study, species delimitation and hybridization history were examined within one Chinese hazel species complex (Corylus chinensis-Corylus fargesii). Two species including four varieties have already been described for this complex, with overlapping distributions.

METHODS

A total of 322 trees from 44 populations of these four varieties across their ranges were sampled for morphological and molecular analyses. Climatic datasets based on 108 geographical locations were used to evaluate their niche differentiations. Flowering phenology was also observed for two co-occurring species or varieties.

KEY RESULTS

Four statistically different phenotypic clusters were revealed, but these clusters were highly inconsistent with the traditional taxonomic groups. All the clusters showed statistically distinct niches, with complete or partial geographical isolation. Only two clusters displayed a distributional overlap, but they had distinct flowering phenologies at the site where they co-occurred. Population-level evidence based on the genotypes of ten simple sequence repeat loci supported four phenotypic clusters. In addition, one cluster was shown to have an admixed genetic composition derived from the other three clusters through repeated historical hybridizations.

CONCLUSIONS

Based on our new evidence, it is better to treat the four clusters identified here as four independent species. One of them was shown to have an admixed genetic composition derived from the other three through repeated historical hybridizations. This study highlights the importance of applying integrative and statistical methods to infer species delimitations and hybridization history. Such a protocol should be adopted widely for future taxonomic studies.

A Cautionary Note on the Use of Genotype Callers in Phylogenomics

Next-generation-sequencing genotype callers are commonly used in studies to call variants from newly sequenced species. However, due to the current availability of genomic resources, it is still common practice to use only one reference genome for a given genus, or even one reference for an entire clade of a higher taxon. The problem with traditional genotype callers, such as the one from GATK, is that they are optimized for variant calling at the population level. However, when these callers are used at the phylogenetic level, the consequences for downstream analyses can be substantial. Here, we performed simulations to compare the performance between the genotype callers of GATK and ATLAS, and present their differences at various phylogenetic scales. We show that the genotype caller of GATK substantially underestimates the number of variants at the phylogenetic level, but not at the population level. We also found that the accuracy of heterozygote calls declines with increasing distance to the reference genome. We quantified this decline and found that it is very sharp in GATK, while ATLAS maintains high accuracy even at moderately divergent species from the reference. We further suggest that efforts should be taken towards acquiring more reference genomes per species, before pursuing high-scale phylogenomic studies. [ATLAS; efficiency of SNP calling; GATK; heterozygote calling; next-generation sequencing; reference genome; variant calling.].

Chloroplast genomes in Populus (Salicaceae): comparisons from an intensively sampled genus reveal dynamic patterns of evolution

The chloroplast is one of two organelles containing a separate genome that codes for essential and distinct cellular functions such as photosynthesis. Given the importance of chloroplasts in plant metabolism, the genomic architecture and gene content have been strongly conserved through long periods of time and as such are useful molecular tools for evolutionary inferences. At present, complete chloroplast genomes from over 4000 species have been deposited into publicly accessible databases. Despite the large number of complete chloroplast genomes, comprehensive analyses regarding genome architecture and gene content have not been conducted for many lineages with complete species sampling. In this study, we employed the genus Populus to assess how more comprehensively sampled chloroplast genome analyses can be used in understanding chloroplast evolution in a broadly studied lineage of angiosperms. We conducted comparative analyses across Populus in order to elucidate variation in key genome features such as genome size, gene number, gene content, repeat type and number, SSR (Simple Sequence Repeat) abundance, and boundary positioning between the four main units of the genome. We found that some genome annotations were variable across the genus owing in part from errors in assembly or data checking and from this provided corrected annotations. We also employed complete chloroplast genomes for phylogenetic analyses including the dating of divergence times throughout the genus. Lastly, we utilized re-sequencing data to describe the variations of pan-chloroplast genomes at the population level for P. euphratica. The analyses used in this paper provide a blueprint for the types of analyses that can be conducted with publicly available chloroplast genomes as well as methods for building upon existing datasets to improve evolutionary inference.

Plant sex chromosomes defy evolutionary models of expanding recombination suppression and genetic degeneration

Hundreds of land plant lineages have independently evolved separate sexes in either gametophytes (dioicy) or sporophytes (dioecy), but 43% of all dioecious angiosperms are found in just 34 entirely dioecious clades, suggesting that their mode of sex determination evolved a long time ago. Here, we review recent insights on the molecular mechanisms that underlie the evolutionary change from individuals that each produce male and female gametes to individuals specializing in the production of just one type of gamete. The canonical model of sex chromosome evolution in plants predicts that two sex-determining genes will become linked in a sex-determining region (SDR), followed by expanding recombination suppression, chromosome differentiation and, ultimately, degeneration. Experimental work, however, is showing that single genes function as master regulators in model systems, such as the liverwort Marchantia and the angiosperms Diospyros and Populus. In Populus, this type of regulatory function has been demonstrated by genome editing. In other systems, including Actinidia, Asparagus and Vitis, two coinherited factors appear to independently regulate female and male function, yet sex chromosome differentiation has remained low. We discuss the best-understood systems and evolutionary pathways to dioecy, and present a meta-analysis of the sizes and ages of SDRs. We propose that limited sexual conflict explains why most SDRs are small and sex chromosomes remain homomorphic. It appears that models of increasing recombination suppression with age do not apply because selection favours mechanisms in which sex determination depends on minimal differences, keeping it surgically precise.