Widespread ancient whole-genome duplications in Malpighiales coincide with Eocene global climatic upheaval

Pan-genome analysis reveals the evolution and diversity of Malus

Malus Mill., a genus of temperate perennial trees with great agricultural and ecological value, has diversified through hybridization, polyploidy and environmental adaptation. Limited genomic resources for wild Malus species have hindered the understanding of their evolutionary history and genetic diversity. We sequenced and assembled 30 high-quality Malus genomes, representing 20 diploids and 10 polyploids across major evolutionary lineages and geographical regions. Phylogenomic analyses revealed ancient gene duplications and conversions, while six newly defined genome types, including an ancestral type shared by polyploid species, facilitated the detection of strong signals for extensive introgressions. The graph-based pan-genome captured shared and species-specific structural variations, facilitating the development of a molecular marker for apple scab resistance. Our pipeline for analyzing selective sweep identified a mutation in MdMYB5 having reduced cold and disease resistance during domestication. This study advances Malus genomics, uncovering genetic diversity and evolutionary insights while enhancing breeding for desirable traits.

Nanoparticles alter the nature and strength of intraploidy and interploidy interactions in plants

Engineered nanoparticles have profound impacts on organisms, yet there is limited understanding of how nanoparticle exposure shapes species interactions that are key to natural community dynamics. By growing plants of the same (intraploidy) and different ploidy levels (interploidy) of Fragaria in axenic microcosms, we examined the influence of nanoparticles on species interactions in polyploid and diploid plants. Under copper oxide (CuO) nanoparticle exposure, polyploids experienced reduced competition and a shift towards facilitation, when growing with both polyploids (the effect of polyploids on polyploids measured by the relative interaction index, RII8x,8x) and diploids (the effect of diploids on polyploids, RII8x,2x). This reduction in competitive interactions in polyploids, in line with the stress gradient hypothesis, was primarily driven by nanoscale effects. In contrast, the strength of competitive interactions (RII8x,8× and RII8x,2x) increased under CuO bulk particles compared to control conditions. Different from polyploids, diploids experienced neutral interactions (RII2x,2x and RII2x,8x) under both nanoparticles and bulk particles. These findings highlight ploidy-specific interaction dynamics and the need to consider species interactions when predicting organismal responses to nanoparticle pollution in ecological communities, providing critical insights for conservation strategies and sustainable nanotechnology applications.

Chromosome number variation and phylogenetic divergence of East Asian Cirsium sect. Onotrophe subsect. Nipponocirsium (Compositae), with a new species from Taiwan

BACKGROUND

This study explored chromosome number variation, phylogenetic divergence, and mechanisms underlying speciation in East Asian thistle Cirsium Mill. sect. Onotrophe (Cass.) DC. subsect. Nipponocirsium Kitam. (Compositae). The study focused on the newly identified species from Taiwan: Cirsium pengii Y.H. Tseng, P.C. Liao & Chih Y. Chang. Utilizing phylotranscriptomic data to reconstruct evolutionary relationships between the Taiwanese and Japanese taxa of Cirsium subsect. Nipponocirsium as well as their divergence times and chromosomal characteristics. Additionally, the chromosome number, morphology, and pollen morphology of the unknown Cirsium taxon are compared with other known subsect. Nipponocirsium taxa from Taiwan.

RESULTS

Phylotranscriptomic analysis reveals a division within subsect. Nipponocirsium into Japanese and Taiwanese clades. In the Taiwanese clade, C. pengii is basal, while C. tatakaense remains monophyletic with other Taiwanese species despite higher genetic diversity. The prevalent chromosome number in this subsection is tetraploid (2n = 4x = 68), common in Japanese taxa, while Taiwanese members have 2n = 4x = 64. Notably, C. pengii has a diploid number (2n = 32), indicating descending dysploidy followed by polyploidization in Taiwan. This polyploidization, driven by glaciations, likely shaped the evolution of Nipponocirsium. Divergence time estimates suggest the separation of Japanese and Taiwanese clades around 0.74 million years ago (Myr) during glacial periods. Cirsium pengii diverged around 0.47 Myr, while tetraploid species C. kawakamii and C. tatakaense diverged around 0.35 Myr. These species likely evolved in separate refugia, with distinct species boundaries confirmed through species delimitation analysis, karyotype, morphology, and pollen morphology comparisons.

CONCLUSIONS

These findings enhance our understanding of chromosomal evolution and speciation within subsect. Nipponocirsium and underscore the importance of integrating transcriptomic data in phylogenetic studies. This study provides a comprehensive framework for further investigations into the genetic diversity and adaptive mechanisms of this ecologically vital group.

Allopolyploidy enhances survival advantages for urban environments in the native plant genus Commelina

BACKGROUND AND AIMS

Urbanization-induced environmental changes affect the geographical distribution of natural plant species. This study focused on how polyploidization, a dynamic genome change, influences the survival and distribution of Commelina communis (Cc) and its subspecies C. communis f. ciliata (Ccfc), which have variable chromosome numbers (e.g. Cc, 2n = 88 for Cc; Ccfc, 2n = 46 for Ccfc). The aim was to investigate polyploidization effects on natural plant distribution in urban environments.

METHODS

The geographical distribution across urban-rural gradients was investigated at a total of 218 sites in Japan. Stomata size and density were measured and compared between Cc and Ccfc. Flow cytometry determined genome size and polyploidy. Chromosome karyotyping was performed using the genomic in situ hybridization (GISH) method.

KEY RESULTS

Urban areas were exclusively dominated by Cc, while Cc and Ccfc coexisted in rural areas. Cc had larger and fewer stomata and a genome size more than twice that of Ccfc. GISH results indicated that Cc possesses Ccfc and another unknown genome, suggesting allopolyploidy.

CONCLUSIONS

Our results show that the ploidy difference affects the geographical distribution, stomata traits and genome size between two distinct taxa in the genus Commelina, C. communis as a neo-tetraploid and C. communis f. ciliata, the diploid. Cc is an allopolyploid and is therefore not only polyploidy but also has an additional genome that provides new sets of genes and alleles, contributing to Cc having enhanced survival potentials in urban environments compared with Ccfc. This is the first investigation to clarify the distribution difference related to urban environments, the difference in stomata traits and genome size, and to study chromosome composition in Commelina species.

Phylogenetic insights into the Salicaceae: The evolution of willows and beyond

The Salicaceae includes approximately 54 genera and over 1,400 species with a cosmopolitan distribution. Members of the family are well-known for their diverse secondary plant metabolites, and they play crucial roles in tropical and temperate forest ecosystems. Phylogenetic reconstruction of the Salicaceae has been historically challenging due to the limitations of molecular markers and the extensive history of hybridization and polyploidy within the family. Our study employs whole-genome sequencing of 74 species to generate an extensive phylogeny of the Salicaceae. We generated two RAD-Seq enriched whole-genome sequence datasets and extracted two additional gene sets corresponding to the universal Angiosperms353 and Salicaceae-specific targeted-capture arrays. We reconstructed maximum likelihood-based molecular phylogenies using supermatrix and coalescent-based supertree approaches. Our fossil-calibrated phylogeny estimates that the Salicaceae originated around 128 million years ago and unravels the complex taxonomic relationships within the family. Our findings confirm the non-monophyly of the subgenus Salix s.l. and further support the merging of subgenera Chamaetia and Vetrix, both of which exhibit intricate patterns within and among different sections. Overall, our study not only enhances our understanding of the evolution of the Salicaceae, but also provides valuable insights into the complex relationships within the family.

Local adaptation, recombination, and the fate of neopolyploids

Polyploidy is widely recognized as an important speciation mechanism because it isolates tetraploids from their diploid progenitors. Polyploidy also provides new genetic material that may facilitate adaptive evolution. However, new mutations are more likely to arise after a neopolyploid has already successfully invaded a population. Thus, the role of adaptive forces in establishing a polyploid remains unclear. One solution to this apparent paradox may lie in the capacity of polyploids to suppress recombination among preexisting locally adapted alleles. The local adaptation mechanism requires that spatially heterogeneous selection acts on multiple loci and that gene flow introduces maladapted alleles to the population where the polyploid forms. The mechanism requires neither strong genetic drift nor any intrinsic benefit of genome doubling and can accommodate any mode of gene action. A unique prediction of the mechanism is that adaptive alleles should predate polyploidization, a pattern consistent with observations from a few well-studied polyploids. The mechanism is also consistent with the coexistence of both diploid and tetraploid cytotypes, fitness heterogeneity among independently derived polyploids, and the prevalence of outcrossing among older polyploids. The local adaptation mechanism also makes novel predictions about circumstances favoring polyploid invasions that can be tested using molecular genetic or comparative approaches.

Influence of polyploidy on morphology and distribution of the Cypress Spurge (Euphorbia cyparissias, Euphorbiaceae)

Polyploidy can cause differences in phenotypic and physiological traits among different cytotypes of the same species. Polyploids may have larger organs or occupy different ecological niches than their diploid counterparts, therefore they are hypothesized to have larger distributions or prosper in stressful environments, such as higher elevations. The Cypress spurge (Euphorbia cyparissias L.; Euphorbiaceae) is a widespread European heteroploid species including di- (2x), tetra- (4x) and hexaploid (6x) cytotypes. We tested the hypotheses that polyploids are more widespread and more abundant at higher elevations and have larger organs than their diploid ancestors in the case of E. cyparissias. We also analysed whether genome downsizing had occurred after polyploidisation. We conducted a comprehensive geographic sampling of 617 populations of E. cyparissias throughout Europe. We estimated their relative genome size using flow cytometry and inferred ploidy level of each population. We scored 13 morphological traits of vegetative and seed characters and performed statistical analyses. The study indicates that polyploidisation facilitated colonisation of new areas in E. cyparissias, where the tetraploids are most widespread, whereas the diploids are limited to putative Pleistocene refugia, mostly in southern Europe. On the other hand, the three ploidies do not differ in their elevational distribution. Although some quantitative morphological traits exhibited an increasing trend with increasing ploidy, most traits did not differ significantly among the three ploidies, and there was no overall phenotypic differentiation among them. Given that individuals of different ploidies thrive in similar habitats across the same elevations, we suggest that ecological segregation following polyploidisation is a more important trigger for morphological differentiation than polyploidisation itself in autopolyploid plants. The study demonstrates that polyploidisation can be crucial for the colonisation of new areas and for range expansion, but it does not necessarily influence elevational distribution nor confer a different phenotype.

Rethinking pathways to the dioecy-polyploidy association: Genera with many dioecious species have fewer polyploids

PREMISE

Numerous studies have found a positive association between dioecy and polyploidy; however, this association presents a theoretical conflict: While polyploids are predicted to benefit from self-reproduction for successful establishment, dioecious species cannot self-reproduce. We propose a theoretical framework to resolve this apparent conflict. We hypothesize that the inability of dioecious species to self-reproduce hinders their establishment as polyploids. We therefore expect that genera with many dioecious species have fewer polyploids, leading to a negative association between polyploidy and dioecy across genera.

METHODS

We used three publicly available databases to determine ploidy and sexual systems for 131 genera and 546 species. We quantified (1) the relationship between the frequency of polyploid species and the frequency of dioecious species across genera, and (2) the proportion of polyploids with hermaphroditism and dioecy across species, adjusting for phylogenetic history.

RESULTS

Across genera, we found a negative relationship between the proportion of polyploids and the proportion of dioecious species, a consistent trend across clades. Across all species, we found that sexual system (dioecious or not) was not associated with polyploidy.

CONCLUSIONS

Polyploids are rare in genera in which the majority of species are dioecious, consistent with the theory that self-reproduction favors polyploid establishment. The low frequency of polyploidy among dioecious species indicates the association is not as widespread as previously suggested. Our findings are consistent with previous studies identifying a positive relationship between the two traits, but only if polyploidy promotes a transition to dioecy, and not the reverse.

Radiating diversification and niche conservatism jointly shape the inverse latitudinal diversity gradient of Potentilla L. (Rosaceae)

BACKGROUND

The latitudinal diversity gradient (LDG), characterized by an increase in species richness from the poles to the equator, is one of the most pervasive biological patterns. However, inverse LDGs, in which species richness peaks in extratropical regions, are also found in some lineages and their causes remain unclear. Here, we test the roles of evolutionary time, diversification rates, and niche conservatism in explaining the inverse LDG of Potentilla (ca. 500 species). We compiled the global distributions of ~ 90% of Potentilla species, and reconstructed a robust phylogenetic framework based on whole-plastome sequences. Next, we analyzed the divergence time, ancestral area, diversification rate, and ancestral niche to investigate the macroevolutionary history of Potentilla.

RESULTS

The genus originated in the Qinghai-Tibet Plateau during the late Eocene and gradually spread to other regions of the Northern Hemisphere posterior to the late Miocene. Rapid cooling after the late Pliocene promoted the radiating diversification of Potentilla. The polyploidization, as well as some cold-adaptive morphological innovations, enhanced the adaptation of Potentilla species to the cold environment. Ancestral niche reconstruction suggests that Potentilla likely originated in a relatively cool environment. The species richness peaks at approximately 45 °N, a region characterized by high diversification rates, and the environmental conditions are similar to the ancestral climate niche. Evolutionary time was not significantly correlated with species richness in the latitudinal gradient.

CONCLUSIONS

Our results suggest that the elevated diversification rates in middle latitude regions and the conservatism in thermal niches jointly determined the inverse LDG in Potentilla. This study highlights the importance of integrating evolutionary and ecological approaches to explain the diversity pattern of biological groups on a global scale.

Interspecific transfer of genetic information through polyploid bridges

Hybridization blurs species boundaries and leads to intertwined lineages resulting in reticulate evolution. Polyploidy, the outcome of whole genome duplication (WGD), has more recently been implicated in promoting and facilitating hybridization between polyploid species, potentially leading to adaptive introgression. However, because polyploid lineages are usually ephemeral states in the evolutionary history of life it is unclear whether WGD-potentiated hybridization has any appreciable effect on their diploid counterparts. Here, we develop a model of cytotype dynamics within mixed-ploidy populations to demonstrate that polyploidy can in fact serve as a bridge for gene flow between diploid lineages, where introgression is fully or partially hampered by the species barrier. Polyploid bridges emerge in the presence of triploid organisms, which despite critically low levels of fitness, can still allow the transfer of alleles between diploid states of independently evolving mixed-ploidy species. Notably, while marked genetic divergence prevents polyploid-mediated interspecific gene flow, we show that increased recombination rates can offset these evolutionary constraints, allowing a more efficient sorting of alleles at higher-ploidy levels before introgression into diploid gene pools. Additionally, we derive an analytical approximation for the rate of gene flow at the tetraploid level necessary to supersede introgression between diploids with nonzero introgression rates, which is especially relevant for plant species complexes, where interspecific gene flow is ubiquitous. Altogether, our results illustrate the potential impact of polyploid bridges on the (re)distribution of genetic material across ecological communities during evolution, representing a potential force behind reticulation.

Evolutionary divergence of subgenomes in common carp provides insights into speciation and allopolyploid success

Hybridization and polyploidization have made great contributions to speciation, heterosis, and agricultural production within plants, but there is still limited understanding and utilization in animals. Subgenome structure and expression reorganization and cooperation post hybridization and polyploidization are essential for speciation and allopolyploid success. However, the mechanisms have not yet been comprehensively assessed in animals. Here, we produced a high-fidelity reference genome sequence for common carp, a typical allotetraploid fish species cultured worldwide. This genome enabled in-depth analysis of the evolution of subgenome architecture and expression responses. Most genes were expressed with subgenome biases, with a trend of transition from the expression of subgenome A during the early stages to that of subgenome B during the late stages of embryonic development. While subgenome A evolved more rapidly, subgenome B contributed to a greater level of expression during development and under stressful conditions. Stable dominant patterns for homoeologous gene pairs both during development and under thermal stress suggest a potential fixed heterosis in the allotetraploid genome. Preferentially expressing either copy of a homoeologous gene at higher levels to confer development and response to stress indicates the dominant effect of heterosis. The plasticity of subgenomes and their shifting of dominant expression during early development, and in response to stressful conditions, provide novel insights into the molecular basis of the successful speciation, evolution, and heterosis of the allotetraploid common carp.

Ferns as facilitators of community recovery following biotic upheaval

The competitive success of ferns has been foundational to hypotheses about terrestrial recolonization following biotic upheaval, from wildfires to the Cretaceous-Paleogene asteroid impact (66 million years ago). Rapid fern recolonization in primary successional environments has been hypothesized to be driven by ferns' high spore production and wind dispersal, with an emphasis on their competitive advantages as so-called disaster taxa. We propose that a competition-based view of ferns is outdated and in need of reexamination in light of growing research documenting the importance of positive interactions (i.e., facilitation) between ferns and other species. Here, we integrate fossil and modern perspectives on fern ecology to propose that ferns act as facilitators of community assemblage following biotic upheaval by stabilizing substrates, enhancing soil properties, and mediating competition. Our reframing of ferns as facilitators has broad implications for both community ecology and ecosystem recovery dynamics, because of ferns' global distribution and habitat diversity.

The impact of multifactorial stress combination on plants, crops, and ecosystems: how should we prepare for what comes next?

The complexity of environmental conditions encountered by plants in the field, or in nature, is gradually increasing due to anthropogenic activities that promote global warming, climate change, and increased levels of pollutants. While in the past it seemed sufficient to study how plants acclimate to one or even two different stresses affecting them simultaneously, the complex conditions developing on our planet necessitate a new approach of studying stress in plants: Acclimation to multiple stress conditions occurring concurrently or consecutively (termed, multifactorial stress combination [MFSC]). In an initial study of the plant response to MFSC, conducted with Arabidopsis thaliana seedlings subjected to an MFSC of six different abiotic stresses, it was found that with the increase in the number and complexity of different stresses simultaneously impacting a plant, plant growth and survival declined, even if the effects of each stress involved in such MFSC on the plant was minimal or insignificant. In three recent studies, conducted with different crop plants, MFSC was found to have similar effects on a commercial rice cultivar, a maize hybrid, tomato, and soybean, causing significant reductions in growth, biomass, physiological parameters, and/or yield traits. As the environmental conditions on our planet are gradually worsening, as well as becoming more complex, addressing MFSC and its effects on agriculture and ecosystems worldwide becomes a high priority. In this review, we address the effects of MFSC on plants, crops, agriculture, and different ecosystems worldwide, and highlight potential avenues to enhance the resilience of crops to MFSC.

Genomic evidence for rediploidization and adaptive evolution following the whole-genome triplication

Whole-genome duplication (WGD), or polyploidy, events are widespread and significant in the evolutionary history of angiosperms. However, empirical evidence for rediploidization, the major process where polyploids give rise to diploid descendants, is still lacking at the genomic level. Here we present chromosome-scale genomes of the mangrove tree Sonneratia alba and the related inland plant Lagerstroemia speciosa. Their common ancestor has experienced a whole-genome triplication (WGT) approximately 64 million years ago coinciding with a period of dramatic global climate change. Sonneratia, adapting mangrove habitats, experienced extensive chromosome rearrangements post-WGT. We observe the WGT retentions display sequence and expression divergence, suggesting potential neo- and sub-functionalization. Strong selection acting on three-copy retentions indicates adaptive value in response to new environments. To elucidate the role of ploidy changes in genome evolution, we improve a model of the polyploidization-rediploidization process based on genomic evidence, contributing to the understanding of adaptive evolution during climate change.

Phylogenomic analyses and chromosome ploidy identification reveal multiple cryptic species in Allium sikkimense complex (Amaryllidaceae)

Polyploidization is a process that typically leads to instantaneous reproductive isolation and has, therefore, been considered as one of the major evolutionary forces in the species-rich Hengduan Mountains (HM), yet this topic remains poorly studied in the region. Allium sikkimense and its relatives (about eight species) compose a natural diploid-polyploid complex with the highest diversity in the HM and adjacent areas. A combination of nuclear ribosomal DNA (nrDNA), plastome, transcriptome, and ploidy identification through chromosome counting and flow cytometry is employed to reconstruct the phylogenetic relationships in this complex and to investigate the frequency and the evolutionary significance of polyploidy in the complex. The plastome failed to resolve the phylogenetic relationships of the different species in the A. sikkimense complex, and the phylogenetic tree based on nrDNA also has limited resolution. However, our study reveals a well-resolved phylogenetic framework for species in the A. sikkimense complex using more than 1,000 orthologous genes from the transcriptome data. Previously recognized morphospecies A. sikkimense are non-monophyletic and comprise at least two independently evolved lineages (i.e., cryptic species), each forming a clade with different diploid species in this complex. The embedded pattern of octoploid A. jichouense and tetraploid A. sp. nov. within different polyploid samples of A. sikkimense supports a possible scenario of budding speciation (via niche divergence). Furthermore, our results reveal that co-occurring species in the A. sikkimense complex usually have different ploidy levels, suggesting that polyploidy is an important process for reproductive isolation of sympatric Allium species. Phylogenetic network analyses suggested that the phylogenetic relationships of the A. sikkimense complex, allowing for reticulation events, always fit the dataset better than a simple bifurcating tree. In addition, the included or exserted filaments, which have long been used to delimit species, are highly unreliable taxonomically due to their extensive parallel and convergent evolution.

Past climate cooling and orogenesis of the Hengduan Mountains have influenced the evolution of Impatiens sect. Impatiens (Balsaminaceae) in the Northern Hemisphere

BACKGROUND

Impatiens sect. Impatiens is distributed across the Northern Hemisphere and has diversified considerably, particularly within the Hengduan Mountains (HDM) in southwest China. Yet, the infra-sectional phylogenetic relationships are not well resolved, largely due to limited taxon sampling and an insufficient number of molecular markers. The evolutionary history of its diversification is also poorly understood. In this study, plastome data and the most complete sampling to date were used to reconstruct a robust phylogenetic framework for this section. The phylogeny was then used to investigate its biogeographical history and diversification patterns, specifically with the aim of understanding the role played by the HDM and past climatic changes in its diversification.

RESULTS

A stable phylogeny was reconstructed that strongly supported both the monophyly of the section and its division into seven major clades (Clades I-VII). Molecular dating and ancestral area reconstruction suggest that sect. Impatiens originated in the HDM and Southeast China around 11.76 Ma, after which different lineages dispersed to Northwest China, temperate Eurasia, and North America, mainly during the Pliocene and Pleistocene. An intercontinental dispersal event from East Asia to western North America may have occurred via the Bering Land Bridge or Aleutian Islands. The diversification rate was high during its early history, especially with the HDM, but gradually decreased over time both within and outside the HDM. Multiple linear regression analysis showed that the distribution pattern of species richness was strongly associated with elevation range, elevation, and mean annual temperature. Finally, ancestral niche analysis indicated that sect. Impatiens originated in a relatively cool, middle-elevation area.

CONCLUSIONS

We inferred the evolutionary history of sect. Impatiens based on a solid phylogenetic framework. The HDM was the primary source or pump of its diversity in the Northern Hemisphere. Orogeny and climate change may have also shaped its diversification rates, as a steady decrease in the diversification rate coincided with the uplift of the HDM and climate cooling. These findings provide insights into the distribution pattern of sect. Impatiens and other plants in the Northern Hemisphere.

Genetic Diversity and Phylogenetic Analysis of Zygophyllum loczyi in Northwest China's Deserts Based on the Resequencing of the Genome

In order to study the genetics of local adaptation in all main deserts of northwest China, whole genomes of 169 individuals were resequenced, which covers 20 populations of Zygophyllum loczyi (Zygophyllales: Zygophylaceae). We describe more than 15 million single nucleotide polymorphisms and numerous InDels. The expected heterozygosity and PIC values associated with local adaptation varied significantly across biogeographic regions. Variation in environmental factors contributes largely to the population genetic structure of Z. loczyi. Bayesian analysis performed with STRUCTURE defined four genetic clusters, while the results of principle component analysis were similar. Our results shows that the Qaidam Desert group appears to be diverging into two branches characterized by significant geographic separation and gene flow with two neighboring deserts. Geological data assume that it is possible that the Taklamakan Desert was the original distribution site, and Z. loczyi could have migrated later on and expanded within other desert areas. The above findings provide insights into the processes involved in biogeography, phylogeny, and differentiation within the northwest deserts of China.

Angiosperm phylogenetic diversity is lower in Africa than South America

Although originating from a common Gondwanan flora, the diversity and composition of the floras of Africa and South America have greatly diverged since continental breakup of Africa from South America now having much higher plant species richness. However, the phylogenetic diversity of the floras and what this tells us about their evolution remained unexplored. We show that for a given species richness and considering land surface area, topography, and present-day climate, angiosperm phylogenetic diversity in South America is higher than in Africa. This relationship holds regardless of whether all climatically matched areas or only matched areas in tropical climates are considered. Phylogenetic diversity is high relative to species richness in refugial areas in Africa and in northwestern South America, once the gateway for immigration from the north. While species richness is strongly influenced by massive plant radiations in South America, we detect a pervasive influence of historical processes on the phylogenetic diversity of both the South American and African floras.

Phylogenomic Analyses Reveal an Allopolyploid Origin of Core Didymocarpinae (Gesneriaceae) Followed by Rapid Radiation

Allopolyploid plants have long been regarded as possessing genetic advantages under certain circumstances due to the combined effects of their hybrid origins and duplicated genomes. However, the evolutionary consequences of allopolyploidy in lineage diversification remain to be fully understood. Here, we investigate the evolutionary consequences of allopolyploidy using 138 transcriptomic sequences of Gesneriaceae, including 124 newly sequenced, focusing particularly on the largest subtribe Didymocarpinae. We estimated the phylogeny of Gesneriaceae using concatenated and coalescent-based methods based on five different nuclear matrices and 27 plastid genes, focusing on relationships among major clades. To better understand the evolutionary affinities in this family, we applied a range of approaches to characterize the extent and cause of phylogenetic incongruence. We found that extensive conflicts between nuclear and chloroplast genomes and among nuclear genes were caused by both incomplete lineage sorting (ILS) and reticulation, and we found evidence of widespread ancient hybridization and introgression. Using the most highly supported phylogenomic framework, we revealed multiple bursts of gene duplication throughout the evolutionary history of Gesneriaceae. By incorporating molecular dating and analyses of diversification dynamics, our study shows that an ancient allopolyploidization event occurred around the Oligocene-Miocene boundary, which may have driven the rapid radiation of core Didymocarpinae.

The duplication of genomes and genetic networks and its potential for evolutionary adaptation and survival during environmental turmoil

The importance of whole-genome duplication (WGD) for evolution is controversial. Whereas some view WGD mainly as detrimental and an evolutionary dead end, there is growing evidence that polyploidization can help overcome environmental change, stressful conditions, or periods of extinction. However, despite much research, the mechanistic underpinnings of why and how polyploids might be able to outcompete or outlive nonpolyploids at times of environmental upheaval remain elusive, especially for autopolyploids, in which heterosis effects are limited. On the longer term, WGD might increase both mutational and environmental robustness due to redundancy and increased genetic variation, but on the short-or even immediate-term, selective advantages of WGDs are harder to explain. Here, by duplicating artificially generated Gene Regulatory Networks (GRNs), we show that duplicated GRNs-and thus duplicated genomes-show higher signal output variation than nonduplicated GRNs. This increased variation leads to niche expansion and can provide polyploid populations with substantial advantages to survive environmental turmoil. In contrast, under stable environments, GRNs might be maladaptive to changes, a phenomenon that is exacerbated in duplicated GRNs. We believe that these results provide insights into how genome duplication and (auto)polyploidy might help organisms to adapt quickly to novel conditions and to survive ecological uproar or even cataclysmic events.

Hybridization and polyploidization generate evolutionary diversification of the genus Buddleja in the Sino-Himalayan region. A commentary on 'Deciphering complex reticulate evolution of Asian Buddleja (Scrophulariaceae): insights into the taxonomy and speciation of polyploid taxa in the Sino-Himalayan region'

This article comments on:

Fengmao Yang, Jia Ge, Yongjie Guo, Richard Olmstead and Weibang Sun. Deciphering complex reticulate evolution of Asian Buddleja (Scrophulariaceae): insights into the taxonomy and speciation of polyploid taxa in the Sino-Himalayan region, Annals of Botany, Volume 132, Issue 1, 01 July 2023, Pages 15–28, https://doi.org/10.1093/aob/mcad022

Evolution of connective glands reveals a new synapomorphy for Malpighiaceae and the hidden potential of staminal glands for Malpighiales systematics

Connective glands are important morphological characters for the taxonomy of some genera of Malpighiaceae, with few recent studies having just elucidated these glands' anatomical and ecological functions. In order to test the systematic relevance of connective glands to the currently accepted phylogenetic informal clades of Malpighiaceae, we characterised the anatomy and/or histochemistry of two-thirds of Malpighiaceae genera and ten species from nine families of Malpighiales to test: 1. Do connective glands occur in the flowers of all informal clades of Malpighiaceae?; and 2. Are they taxonomically relevant to characterise those clades? We sampled 25 genera and 26 species of Malpighiaceae, processing their anthers using traditional anatomical methods and characterising their glands using light microscopy and SEM imaging. Selected species were subjected to histochemical tests, and an additional 21 genera and 33 species of Malpighiaceae and nine families (ten species) of Malpighiales were included in our sampling from the literature. Three anatomical characters were scored, coded and mapped using Maximum Likelihood methods onto the molecular phylogeny of Malpighiaceae. All sampled species of Malpighiaceae showed connective glands characterised as epidermal or trichomal elaiophores. Our character-mapping analyses recovered connective elaiophores as a new synapomorphy for Malpighiaceae. Different types of epidermal or trichomal elaiophores were recovered as homoplasies for the Christianella and Banisteriopsis clades and the genera Byrsonima, Camarea and Cottsia. Our analyses also recovered the glands' place of insertion in the stamen and the exudate type as potential new synapomorphies or homoplasies for the families of Malpighiales sampled. Our results propose the connective elaiophores as a new synapomorphy for Malpighiaceae and hypothesise the role that different staminal glands might play in the systematics of Malpighiales. Further comprehensive anatomical studies are still needed for the staminal glands of most families of this order to shed new light on the patterns recovered in our study.

Genomic analyses provide insights into the polyploidization-driven herbicide adaptation in Leptochloa weeds

Polyploidy confers a selective advantage under stress conditions; however, whether polyploidization mediates enhanced herbicide adaptation remains largely unknown. Tetraploid Leptochloa chinensis is a notorious weed in the rice ecosystem, causing severe yield loss in rice. In China, L. chinensis has only one sister species, the diploid L. panicea, whose damage is rarely reported. To gain insights into the effects of polyploidization on herbicide adaptation, we first assembled a high-quality genome of L. panicea and identified genome structure variations with L. chinensis. Moreover, we identified herbicide-resistance genes specifically expanded in L. chinensis, which may confer a greater herbicide adaptability in L. chinensis. Analysis of gene retention and loss showed that five herbicide target-site genes and several herbicide nontarget-site resistance gene families were retained during polyploidization. Notably, we identified three pairs of polyploidization-retained genes including LcABCC8, LcCYP76C1 and LcCYP76C4 that may enhance herbicide resistance. More importantly, we found that both copies of LcCYP76C4 were under herbicide selection during the spread of L. chinensis in China. Furthermore, we identified another gene potentially involved in herbicide resistance, LcCYP709B2, which is also retained during polyploidization and under selection. This study provides insights into the genomic basis of the enhanced herbicide adaptability of Leptochloa weeds during polyploidization and provides guidance for the precise and efficient control of polyploidy weeds.

Phylotranscriptomic analyses reveal multiple whole-genome duplication events, the history of diversification and adaptations in the Araceae

BACKGROUND AND AIMS

The Araceae are one of the most diverse monocot families with numerous morphological and ecological novelties. Plastid and mitochondrial genes have been used to investigate the phylogeny and to interpret shifts in the pollination biology and biogeography of the Araceae. In contrast, the role of whole-genome duplication (WGD) in the evolution of eight subfamilies remains unclear.

METHODS

New transcriptomes or low-depth whole-genome sequences of 65 species were generated through Illumina sequencing. We reconstructed the phylogenetic relationships of Araceae using concatenated and species tree methods, and then estimated the age of major clades using TreePL. We inferred the WGD events by Ks and gene tree methods. We investigated the diversification patterns applying time-dependent and trait-dependent models. The expansions of gene families and functional enrichments were analysed using CAFE and InterProScan.

KEY RESULTS

Gymnostachydoideae was the earliest diverging lineage followed successively by Orontioideae, Lemnoideae and Lasioideae. In turn, they were followed by the clade of 'bisexual climbers' comprised of Pothoideae and Monsteroideae, which was resolved as the sister to the unisexual flowers clade of Zamioculcadoideae and Aroideae. A special WGD event ψ (psi) shared by the True-Araceae clade occurred in the Early Cretaceous. Net diversification rates first declined and then increased through time in the Araceae. The best diversification rate shift along the stem lineage of the True-Araceae clade was detected, and net diversification rates were enhanced following the ψ-WGD. Functional enrichment analyses revealed that some genes, such as those encoding heat shock proteins, glycosyl hydrolase and cytochrome P450, expanded within the True-Araceae clade.

CONCLUSIONS

Our results improve our understanding of aroid phylogeny using the large number of single-/low-copy nuclear genes. In contrast to the Proto-Araceae group and the lemnoid clade adaption to aquatic environments, our analyses of WGD, diversification and functional enrichment indicated that WGD may play a more important role in the evolution of adaptations to tropical, terrestrial environments in the True-Araceae clade. These insights provide us with new resources to interpret the evolution of the Araceae.

Inference of Ancient Polyploidy Using Transcriptome Data

Polyploidizations, or whole-genome duplications (WGDs), in plants have increased biological complexity, facilitated evolutionary innovation, and likely enabled adaptation under harsh conditions. Besides genomic data, transcriptome data have been widely employed to detect WGDs, due to their efficient accessibility to the gene space of a species. Age distributions based on synonymous substitutions (so-called K_S age distributions) for paralogs assembled from transcriptome data have identified numerous WGDs in plants, paving the way for further studies on the importance of WGDs for the evolution of seed and flowering plants. However, it is still unclear how transcriptome-based age distributions compare to those based on genomic data. In this chapter, we implemented three different de novo transcriptome assembly pipelines with two popular assemblers, i.e., Trinity and SOAPdenovo-Trans. We selected six plant species with published genomes and transcriptomes to evaluate how assembled transcripts from different pipelines perform when using K_S distributions to detect previously documented WGDs in the six species. Further, using genes predicted in each genome as references, we evaluated the effects of missing genes, gene family clustering, and de novo assembled transcripts on the transcriptome-based K_S distributions. Our results show that, although the transcriptome-based K_S distributions differ from the genome-based ones with respect to their shapes and scales, they are still reasonably reliable for unveiling WGDs, except in species where most duplicates originated from a recent WGD. We also discuss how to overcome some possible pitfalls when using transcriptome data to identify WGDs.

Detecting signals of adaptive evolution in grape plastomes with a focus on the Cretaceous-Palaeogene (K/Pg) transition

BACKGROUND AND AIMS

Although plastid genes are widely used in phylogenetic studies, signals of positive selection have been scarcely investigated in the grape family. The plastomes from 91 accessions of Vitaceae were examined to understand the extent to which positive selection is present and to identify which genes are involved. Moreover, the changes through time of genes under episodic positive selection were investigated and the hypothesis of an adaptive process following the Cretaceous-Palaeogene (K/Pg) transition about 66 million years ago was tested.

METHODS

Different codon-substitution models were used to assess pervasive and episodic positive selection events on 70 candidate plastid genes. Divergence times between lineages were estimated and stochastic character mapping analysis was used to simulate variation over time of the genes found to be under episodic positive selection.

KEY RESULTS

A total of 20 plastid genes (29 %) showed positive selection. Among them, 14 genes showed pervasive signatures of positive selection and nine genes showed episodic signatures of positive selection. In particular, four of the nine genes (psbK, rpl20, rpoB, rps11) exhibited a similar pattern showing an increase in the rate of variation close to the K/Pg transition.

CONCLUSION

Multiple analyses have shown that the grape family has experienced ancient and recent positive selection events and that the targeted genes are involved in essential functions such as photosynthesis, self-replication and metabolism. Our results are consistent with the idea that the K/Pg transition has favoured an increased rate of change in some genes. Intense environmental perturbations have influenced the rapid diversification of certain lineages, and new mutations arising on some plastid genes may have been fixed by natural selection over the course of many generations.

Tree2GD: a phylogenomic method to detect large-scale gene duplication events

MOTIVATION

Whole-genome duplication events have long been discovered throughout the evolution of eukaryotes, contributing to genome complexity and biodiversity and leaving traces in the descending organisms. Therefore, an accurate and rapid phylogenomic method is needed to identify the retained duplicated genes on various lineages across the target taxonomy.

RESULTS

Here, we present Tree2GD, an integrated method to identify large-scale gene duplication events by automatically perform multiple procedures, including sequence alignment, recognition of homolog, gene tree/species tree reconciliation, Ks distribution of gene duplicates and synteny analyses. Application of Tree2GD on 2 datasets, 12 metazoan genomes and 68 angiosperms, successfully identifies all reported whole-genome duplication events exhibited by these species, showing effectiveness and efficiency of Tree2GD on phylogenomic analyses of large-scale gene duplications.

AVAILABILITY AND IMPLEMENTATION

Tree2GD is written in Python and C++ and is available at https://github.com/Dee-chen/Tree2gd.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Convergent evolution of AP2/ERF III and IX subfamilies through recurrent polyploidization and tandem duplication during eudicot adaptation to paleoenvironmental changes

Whole-genome duplication (WGD or polyploidization) has been suggested as a genetic contributor to angiosperm adaptation to environmental changes. However, many eudicot lineages did not undergo recent WGD (R-WGD) around and/or after the Cretaceous-Paleogene (K-Pg) boundary, times of severe environmental changes; how those plants survived has been largely ignored. Here, we collected 22 plants from major branches of the eudicot phylogeny and classified them into two groups according to the occurrence or absence of R-WGD: 12 R-WGD-containing plants (R-WGD-Y) and 10 R-WGD-lacking plants (R-WGD-N). Subsequently, we identified 496 gene-rich families in R-WGD-Y and revealed that members of the AP2/ERF transcription factor family were convergently over-retained after R-WGDs and showed exceptional cold stimulation. The evolutionary trajectories of the AP2/ERF family were then compared between R-WGD-Y and R-WGD-N to reveal convergent expansions of the AP2/ERF III and IX subfamilies through recurrent independent WGDs and tandem duplications (TDs) after the radiation of the plants. The expansions showed coincident enrichments in- times around and/or after the K-Pg boundary, when global cooling was a major environmental stressor. Consequently, convergent expansions and co-retentions of AP2/ERF III C-repeat binding factor (CBF) duplicates and their regulons in different eudicot lineages contributed to the rewiring of cold-specific regulatory networks. Moreover, promoter analysis of cold-responsive AP2/ERF genes revealed an underlying cis-regulatory code (G-box: CACGTG). We propose a seesaw model of WGDs and TDs in the convergent expansion of AP2/ERF III and IX genes that has contributed to eudicot adaptation during paleoenvironmental changes, and we suggest that TD may be a reciprocal/alternative mechanism for genetic innovation in plants that lack WGD.

Phylogenomics as an effective approach to untangle cross-species hybridization event: A case study in the family Nymphaeaceae

Hybridization is common and considered as an important evolutionary force to increase intraspecific genetic diversity. Detecting hybridization events is crucial for understanding the evolutionary history of species and further improving molecular breeding. The studies on identifying hybridization events through the phylogenomic approach are still limited. We proposed the conception and method of identifying allopolyploidy events by phylogenomics. The reconciliation and summary of nuclear multi-labeled gene family trees were adopted to untangle hybridization events from next-generation data in our novel phylogenomic approach. Given horticulturalists’ relatively clear cultivated crossbreeding history, the water lily family is a suitable case for examining recent allopolyploidy events. Here, we reconstructed and confirmed the well-resolved nuclear phylogeny for the Nymphaeales family in the context of geological time as a framework for identifying hybridization signals. We successfully identified two possible allopolyploidy events with the parental lineages for the hybrids in the family Nymphaeaceae based on summarization from multi-labeled gene family trees of Nymphaeales. The lineages where species Nymphaea colorata and Nymphaea caerulea are located may be the progenitors of horticultural cultivated species Nymphaea ‘midnight’ and Nymphaea ‘Woods blue goddess’. The proposed hybridization hypothesis is also supported by horticultural breeding records. Our methodology can be widely applied to identify hybridization events and theoretically facilitate the genome breeding design of hybrid plants.

Benefits of alignment quality-control processing steps and an Angiosperms353 phylogenomics pipeline applied to the Celastrales

We examined the impact of successive alignment quality-control steps on downstream phylogenomic analyses. We applied a recently published phylogenomics pipeline that was developed for the Angiosperms353 target-sequence-capture probe set to the flowering plant order Celastrales. Our final dataset consists of 158 species, including at least one exemplar from all 109 currently recognized Celastrales genera. We performed nine quality-control steps and compared the inferred resolution, branch support, and topological congruence of the inferred gene and species trees with those generated after each of the first six steps. We describe and justify each of our quality-control steps, including manual masking, in detail so that they may be readily applied to other lineages. We found that highly supported clades could generally be relied upon even if stringent orthology and alignment quality-control measures had not been applied. But separate instances were identified, for both concatenation and coalescence, wherein a clade was highly supported before manual masking but then subsequently contradicted. These results are generally reassuring for broad-scale analyses that use phylogenomics pipelines, but also indicate that we cannot rely exclusively on these analyses to conclude how challenging phylogenetic problems are best resolved.

Diversity, phylogeny, and adaptation of bryophytes: insights from genomic and transcriptomic data

Bryophytes including mosses, liverworts, and hornworts are among the earliest land plants, and occupy a crucial phylogenetic position to aid in the understanding of plant terrestrialization. Despite their small size and simple structure, bryophytes are the second largest group of extant land plants. They live ubiquitously in various habitats and are highly diversified, with adaptive strategies to modern ecosystems on Earth. More and more genomes and transcriptomes have been assembled to address fundamental questions in plant biology. Here, we review recent advances in bryophytes associated with diversity, phylogeny, and ecological adaptation. Phylogenomic studies have provided increasing supports for the monophyly of bryophytes, with hornworts sister to the Setaphyta clade including liverworts and mosses. Further comparative genomic analyses revealed that multiple whole-genome duplications might have contributed to the species richness and morphological diversity in mosses. We highlight that the biological changes through gene gain or neofunctionalization that primarily evolved in bryophytes have facilitated the adaptation to early land environments; among the strategies to adapt to modern ecosystems in bryophytes, desiccation tolerance is the most remarkable. More genomic information for bryophytes would shed light on key mechanisms for the ecological success of these 'dwarfs' in the plant kingdom.

Species Tree Estimation and the Impact of Gene Loss Following Whole-Genome Duplication

Whole-genome duplication (WGD) occurs broadly and repeatedly across the history of eukaryotes and is recognized as a prominent evolutionary force, especially in plants. Immediately following WGD, most genes are present in two copies as paralogs. Due to this redundancy, one copy of a paralog pair commonly undergoes pseudogenization and is eventually lost. When speciation occurs shortly after WGD; however, differential loss of paralogs may lead to spurious phylogenetic inference resulting from the inclusion of pseudoorthologs–paralogous genes mistakenly identified as orthologs because they are present in single copies within each sampled species. The influence and impact of including pseudoorthologs versus true orthologs as a result of gene extinction (or incomplete laboratory sampling) are only recently gaining empirical attention in the phylogenomics community. Moreover, few studies have yet to investigate this phenomenon in an explicit coalescent framework. Here, using mathematical models, numerous simulated data sets, and two newly assembled empirical data sets, we assess the effect of pseudoorthologs on species tree estimation under varying degrees of incomplete lineage sorting (ILS) and differential gene loss scenarios following WGD. When gene loss occurs along the terminal branches of the species tree, alignment-based (BPP) and gene-tree-based (ASTRAL, MP-EST, and STAR) coalescent methods are adversely affected as the degree of ILS increases. This can be greatly improved by sampling a sufficiently large number of genes. Under the same circumstances, however, concatenation methods consistently estimate incorrect species trees as the number of genes increases. Additionally, pseudoorthologs can greatly mislead species tree inference when gene loss occurs along the internal branches of the species tree. Here, both coalescent and concatenation methods yield inconsistent results. These results underscore the importance of understanding the influence of pseudoorthologs in the phylogenomics era. [Coalescent method; concatenation method; incomplete lineage sorting; pseudoorthologs; single-copy gene; whole-genome duplication.]

A phylotranscriptome study using silica gel-dried leaf tissues produces an updated robust phylogeny of Ranunculaceae

The utility of transcriptome data in plant phylogenetics has gained popularity in recent years. However, because RNA degrades much more easily than DNA, the logistics of obtaining fresh tissues has become a major limiting factor for widely applying this method. Here, we used Ranunculaceae to test whether silica-dried plant tissues could be used for RNA extraction and subsequent phylogenomic studies. We sequenced 27 transcriptomes, 21 from silica gel-dried (SD-samples) and six from liquid nitrogen-preserved (LN-samples) leaf tissues, and downloaded 27 additional transcriptomes from GenBank. Our results showed that although the LN-samples produced slightly better reads than the SD-samples, there were no significant differences in RNA quality and quantity, assembled contig lengths and numbers, and BUSCO comparisons between two treatments. Using these data, we conducted phylogenomic analyses, including concatenated- and coalescent-based phylogenetic reconstruction, molecular dating, coalescent simulation, phylogenetic network estimation, and whole genome duplication (WGD) inference. The resulting phylogeny was consistent with previous studies with higher resolution and statistical support. The 11 core Ranunculaceae tribes grouped into two chromosome type clades (T- and R-types), with high support. Discordance among gene trees is likely due to hybridization and introgression, ancient genetic polymorphism and incomplete lineage sorting. Our results strongly support one ancient hybridization event within the R-type clade and three WGD events in Ranunculales. Evolution of the three Ranunculaceae chromosome types is likely not directly related to WGD events. By clearly resolving the Ranunculaceae phylogeny, we demonstrated that SD-samples can be used for RNA-seq and phylotranscriptomic studies of angiosperms.

Phylogenomic Analysis Reconstructed the Order Matoniales from Paleopolyploidy Veil

Phylogenetic conflicts limit our understanding of the evolution of terrestrial life under multiple whole genome duplication events, and the phylogeny of early terrestrial plants remains full of controversy. Although much incongruence has been solved with so-called robust topology based on single or lower copy genes, the evolutionary mechanisms behind phylogenetic conflicts such as polyploidization remain poorly understood. Here, through decreasing the effects of polyploidization and increasing the samples of species, which represent all four orders and eight families that comprise early leptosporangiate ferns, we have reconstructed a robust phylogenetic tree and network with 1125 1-to-1 orthologs based on both coalescent and concatenation methods. Our data consistently suggest that Matoniales, as a monophyletic lineage including Matoniaceae and Dipteridaceae, should be redefined as an ordinal rank. Furthermore, we have identified and located at least 11 whole-genome duplication events within the evolutionary history of four leptosporangiates lineages, and associated polyploidization with higher speciation rates and mass extinction events. We hypothesize that paleopolyploidization may have enabled leptosporangiate ferns to survive during mass extinction events at the end Permian period and then flourish throughout the Mesozoic era, which is supported by extensive fossil records. Our results highlight how ancient polyploidy can result in rapid species radiation, thus causing phylogenetic conflicts yet allowing plants to survive and thrive during mass extinction events.

Polyploidy in urban environments

Polyploidy is a major driver of evolutionary change in plants, but many aspects of polyploidy in natural systems remain enigmatic. We argue that urban landscapes present an unprecedented opportunity to observe polyploidy in action. Integrating polyploid biology and urban evolutionary ecology, we discuss multiple factors expected to promote polyploid formation, establishment, and persistence in urban systems. We develop a predictive framework for the contemporary ecology and evolution of polyploid plants in cities, and through this novel perspective propose that studying polyploidy in an urban context could lead to breakthroughs in understanding fundamental processes in polyploid evolution. We conclude by highlighting the potential consequences of polyploidy in urban environments, and outline a roadmap for research into this currently unexplored field.

Lycophyte transcriptomes reveal two whole-genome duplications in Lycopodiaceae: Insights into the polyploidization of Phlegmariurus

Lycophytes are an ancient clade of the non-flowering vascular plants with chromosome numbers that vary from tens to hundreds. They are an excellent study system for examining whole-genome duplications (WGDs), or polyploidization, in spore-dispersed vascular plants. However, a lack of genome sequence data limits the reliable detection of very ancient WGDs, small-scale duplications (SSDs), and recent WGDs. Here, we integrated phylogenomic analysis and the distribution of synonymous substitutions per synonymous sites (Ks) of the transcriptomes of 13 species of lycophytes to identify, locate, and date multiple WGDs in the lycophyte family Lycopodiaceae. Additionally, we examined the genus Phlegmariurus for signs of genetic discordance, which can provide valuable insight into the underlying causes of such conflict (e.g., hybridization, incomplete lineage sorting, or horizontal gene transfer).We found strong evidence that two WGD events occurred along the phylogenetic backbone of Lycopodiaceae, with one occurring in the common ancestor of extant Phlegmariurus (Lycopodiaceae) approximately 22-23 million years ago (Mya) and the other occurring in the common ancestor of Lycopodiaceae around 206-214 Mya. Interestingly, we found significant genetic discordance in the genus Phlegmariurus, indicating that the genus has a complex evolutionary history. This study provides molecular evidence for multiple WGDs in Lycopodiaceae and offers phylogenetic clues to the evolutionary history of Lycopodiaceae.

A Reappraisal of Polyploidy Events in Grasses (Poaceae) in a Rapidly Changing World

Around 80% of megaflora species became extinct at the Cretaceous–Paleogene (K–Pg) boundary. Subsequent polyploidy events drove the survival of thousands of plant species and played a significant historical role in the development of the most successful modern cereal crops. However, current and rapid global temperature change poses an urgent threat to food crops worldwide, including the world’s big three cereals: rice, wheat, and maize, which are members of the grass family, Poaceae. Some minor cereals from the same family (such as teff) have grown in popularity in recent years, but there are important knowledge gaps regarding the similarities and differences between major and minor crops, including how polyploidy affects their biological processes under natural and (a)biotic stress conditions and thus the potential to harness polyploidization attributes for improving crop climate resilience. This review focuses on the impact of polyploidy events on the Poaceae family, which includes the world’s most important food sources, and discusses the past, present, and future of polyploidy research for major and minor crops. The increasing accessibility to genomes of grasses and their wild progenitors together with new tools and interdisciplinary research on polyploidy can support crop improvement for global food security in the face of climate change.

Untargeted Metabolomics Sheds Light on the Diversity of Major Classes of Secondary Metabolites in the Malpighiaceae Botanical Family

Natural products produced by plants are one of the most investigated natural sources, which substantially contributed to the development of the natural products field. Even though these compounds are widely explored, the literature still lacks comprehensive investigations aiming to explore the evolution of secondary metabolites produced by plants, especially if classical methodologies are employed. The development of sensitive hyphenated techniques and computational tools for data processing has enabled the study of large datasets, being valuable assets for chemosystematic studies. Here, we describe a strategy for chemotaxonomic investigations using the Malpighiaceae botanical family as a model. Our workflow was based on MS/MS untargeted metabolomics, spectral searches, and recently described in silico classification tools, which were mapped into the latest molecular phylogeny accepted for this family. The metabolomic analysis revealed that different ionization modes and extraction protocols significantly impacted the chemical profiles, influencing the chemotaxonomic results. Spectral searches within public databases revealed several clades or genera-specific molecular families, being potential chemical markers for these taxa, while the in silico classification tools were able to expand the Malpighiaceae chemical space. The classes putatively annotated were used for ancestral character reconstructions, which recovered several classes of metabolites as homoplasies (i.e., non-exclusive) or synapomorphies (i.e., exclusive) for all sampled clades and genera. Our workflow combines several approaches to perform a comprehensive evolutionary chemical study. We expect it to be used on further chemotaxonomic investigations to expand chemical knowledge and reveal biological insights for compounds classes in different biological groups.

Phylotranscriptomics Resolves the Phylogeny of Pooideae and Uncovers Factors for Their Adaptive Evolution

Adaptation to cool climates has occurred several times in different angiosperm groups. Among them, Pooideae, the largest grass subfamily with ∼3,900 species including wheat and barley, have successfully occupied many temperate regions and play a prominent role in temperate ecosystems. To investigate possible factors contributing to Pooideae adaptive evolution to cooling climates, we performed phylogenetic reconstruction using five gene sets (with 1,234 nuclear genes and their subsets) from 157 transcriptomes/genomes representing all 15 tribes and 24 of 26 subtribes. Our phylogeny supports the monophyly of all tribes (except Diarrheneae) and all subtribes with at least two species, with strongly supported resolution of their relationships. Molecular dating suggests that Pooideae originated in the late Cretaceous, with subsequent divergences under cooling conditions first among many tribes from the early middle to late Eocene and again among genera in the middle Miocene and later periods. We identified a cluster of gene duplications (CGD5) shared by the core Pooideae (with 80% Pooideae species) near the Eocene–Oligocene transition, coinciding with the transition from closed to open habitat and an upshift of diversification rate. Molecular evolutionary analyses homologs of CBF for cold resistance uncovered tandem duplications during the core Pooideae history, dramatically increasing their copy number and possibly promoting adaptation to cold habitats. Moreover, duplication of AP1/FUL-like genes before the Pooideae origin might have facilitated the regulation of the vernalization pathway under cold environments. These and other results provide new insights into factors that likely have contributed to the successful adaptation of Pooideae members to temperate regions.

Pre- and postzygotic mechanisms preventing hybridization in co-occurring species of the Impatiens purpureoviolacea complex

In the species-rich genus Impatiens, few natural hybrids are known, even though closely related species often occur sympatrically. In this study, we aim to bridge the gap between micro- and macro-evolution to disentangle pre- and postzygotic mechanisms that may prevent hybridization in the Impatiens purpureoviolacea complex from Central Africa. We analyzed habitat types, species distribution, pollination syndromes, pollinator dependency, genome sizes, and chromosome numbers of seven out of the ten species of the complex as well as of one natural hybrid and reconstructed the ancestral chromosome numbers of the complex. Several species of the complex occur in sympatry or geographically very close to each other. All of them are characterized by pre- and/or postzygotic mechanisms potentially preventing hybridization. We found four independent polyploidization events within the complex. The only known natural hybrid always appears as single individual and is self-fertile. But the plants resulting from self-pollinated seeds often die shortly after first flowering. These results indicate that the investigated mechanisms in combination may effectively but not absolutely prevent hybridization in Impatiens and probably occur in other genera with sympatric species as well.

Frequent ploidy changes in Salicaceae indicates widespread sharing of the salicoid whole genome duplication by the relatives of Populus L. and Salix L

BACKGROUNDS

Populus and Salix belong to Salicaceae and are used as models to investigate woody plant physiology. The variation of karyotype and nuclear DNA content can partly reflect the evolutionary history of the whole genome, and can provide critical information for understanding, predicting, and potentially ameliorating the woody plant traits. Therefore, it is essential to study the chromosome number (CN) and genome size in detail to provide information for revealing the evolutionary process of Salicaceae.

RESULTS

In this study, we report the somatic CNs of seventeen species from eight genera in Salicaceae. Of these, CNs for twelve species and for five genera are reported for the first time. Among the three subfamilies of Salicaceae, the available data indicate CN in Samydoideae is n = 21, 22, 42. The only two genera, Dianyuea and Scyphostegia, in Scyphostegioideae respectively have n = 9 and 18. In Salicoideae, Populus, Salix and five genera closely related to them (Bennettiodendron, Idesia, Carrierea, Poliothyrsis, Itoa) are based on relatively high CNs from n = 19, 20, 21, 22 to n = 95 in Salix. However, the other genera of Salicoideae are mainly based on relatively low CNs of n = 9, 10, 11. The genome sizes of 35 taxa belonging to 14 genera of Salicaceae were estimated. Of these, the genome sizes of 12 genera and all taxa except Populus euphratica are first reported. Except for Dianyuea, Idesia and Bennettiodendron, all examined species have relatively small genome sizes of less than 1 pg, although polyploidization exists.

CONCLUSIONS

The variation of CN and genome size across Salicaceae indicates frequent ploidy changes and a widespread sharing of the salicoid whole genome duplication (WGD) by the relatives of Populus and Salix. The shrinkage of genome size after WGD indicates massive loss of genomic components. The phylogenetic asymmetry in clade of Populus, Salix, and their close relatives suggests that there is a lag-time for the subsequent radiations after the salicoid WGD event. Our results provide useful data for studying the evolutionary events of Salicaceae.

A genome sequence resource for the genus Passiflora, the genome of the wild diploid species Passiflora organensis

The genus Passiflora comprises a large group of plants popularly known as passionfruit, much appreciated for their exotic flowers and edible fruits. The species (∼500) are morphologically variable (e.g., growth habit, size, and color of flowers) and are adapted to distinct tropical ecosystems. In this study, we generated the genome of the wild diploid species Passiflora organensis Gardner by adopting a hybrid assembly approach. Passiflora organensis has a small genome of 259 Mbp and a heterozygosity rate of 81%, consistent with its reproductive system. Most of the genome sequences could be integrated into its chromosomes with cytogenomic markers (satellite DNA) as references. The repeated sequences accounted for 58.55% of the total DNA analyzed, and the Tekay lineage was the prevalent retrotransposon. In total, 25,327 coding genes were predicted. Passiflora organensis retains 5,609 singletons and 15,671 gene families. We focused on the genes potentially involved in the locus determining self-incompatibility and the MADS-box gene family, allowing us to infer expansions and contractions within specific subfamilies. Finally, we recovered the organellar DNA. Structural rearrangements and two mitoviruses, besides relics of other mobile elements, were found in the chloroplast and mt-DNA molecules, respectively. This study presents the first draft genome assembly of a wild Passiflora species, providing a valuable sequence resource for genomic and evolutionary studies on the genus, and support for breeding cropped passionfruit species.

Genotyping of polyploid plants using quantitative PCR: application in the breeding of white-fleshed triploid loquats (Eriobotrya japonica)

BACKGROUND

Ploidy manipulation is effective in seedless loquat breeding, in which flesh color is a key agronomic and economic trait. Few techniques are currently available for detecting the genotypes of polyploids in plants, but this ability is essential for most genetic research and molecular breeding.

RESULTS

We developed a system for genotyping by quantitative PCR (qPCR) that allowed flesh color genotyping in multiple tetraploid and triploid loquat varieties (lines). The analysis of 13 different ratios of DNA mixtures between two homozygous diploids (AA and aa) showed that the proportion of allele A has a high correlation (R² = 0.9992) with parameter b [b = a₁/(a₁ + a₂)], which is derived from the two normalized allele signals (a₁ and a₂) provided by qPCR. Cluster analysis and variance analysis from simulating triploid and tetraploid hybrids provided completely correct allelic configurations. Four genotypes (AAA, AAa, Aaa, aaa) were found in triploid loquats, and four (AAAA, AAAa, AAaa, Aaaa; absence of aaaa homozygotes) were found in tetraploid loquats. DNA markers analysis showed that the segregation of flesh color in all F₁ hybrids conformed to Mendel's law. When tetraploid B431 was the female parent, more white-fleshed triploids occurred among the progeny.

CONCLUSIONS

qPCR can detect the flesh color genotypes of loquat polyploids and provides an alternative method for analyzing polyploid genotype and breeding, dose effects and allele-specific expression.

Polyploidy: an evolutionary and ecological force in stressful times

Polyploidy has been hypothesized to be both an evolutionary dead-end and a source for evolutionary innovation and species diversification. Although polyploid organisms, especially plants, abound, the apparent nonrandom long-term establishment of genome duplications suggests a link with environmental conditions. Whole-genome duplications seem to correlate with periods of extinction or global change, while polyploids often thrive in harsh or disturbed environments. Evidence is also accumulating that biotic interactions, for instance, with pathogens or mutualists, affect polyploids differently than nonpolyploids. Here, we review recent findings and insights on the effect of both abiotic and biotic stress on polyploids versus nonpolyploids and propose that stress response in general is an important and even determining factor in the establishment and success of polyploidy.

Which factors contribute most to genome size variation within angiosperms?

Genome size varies greatly across the flowering plants and has played an important role in shaping their evolution. It has been reported that many factors correlate with the variation in genome size, but few studies have systematically explored this at the genomic level. Here, we scan genomic information for 74 species from 74 families in 38 orders covering the major groups of angiosperms (the taxonomic information was acquired from the latest Angiosperm Phylogeny Group (APG IV) system) to evaluate the correlation between genome size variation and different genome characteristics: polyploidization, different types of repeat sequence content, and the dynamics of long terminal repeat retrotransposons (LTRs). Surprisingly, we found that polyploidization shows no significant correlation with genome size, while LTR content demonstrates a significantly positive correlation. This may be due to genome instability after polyploidization, and since LTRs occupy most of the genome content, it may directly result in most of the genome variation. We found that the LTR insertion time is significantly negatively correlated with genome size, which may reflect the competition between insertion and deletion of LTRs in each genome, and that the old insertions are usually easy to recognize and eliminate. We also noticed that most of the LTR burst occurred within the last 3 million years, a timeframe consistent with the violent climate fluctuations in the Pleistocene. Our findings enhance our understanding of genome size evolution within angiosperms, and our methods offer immediate implications for corresponding research in other datasets.

The chromosome-based lavender genome provides new insights into Lamiaceae evolution and terpenoid biosynthesis

The aromatic shrub Lavandula angustifolia produces various volatile terpenoids that serve as resources for essential oils and function in plant-insect communication. To better understand the genetic basis of the terpenoid diversity in lavender, we present a high-quality reference genome for the Chinese lavender cultivar "Jingxun 2" using PacBio and Hi-C technologies to anchor the 894.50 Mb genome assembly into 27 pseudochromosomes. In addition to the γ triplication event, lavender underwent two rounds of whole-genome duplication (WGD) during the Eocene-Oligocene (29.6 MYA) and Miocene-Pliocene (6.9 MYA) transitions. As a result of tandem duplications and lineage-specific WGDs, gene families related to terpenoid biosynthesis in lavender are substantially expanded compared to those of five other species in Lamiaceae. Many terpenoid biosynthesis transcripts are abundant in glandular trichomes. We further integrated the contents of ecologically functional terpenoids and coexpressed terpenoid biosynthetic genes to construct terpenoid-gene networks. Typical gene clusters, including TPS-TPS, TPS-CYP450, and TPS-BAHD, linked with compounds that primarily function as attractants or repellents, were identified by their similar patterns of change during flower development or in response to methyl jasmonate. Comprehensive analysis of the genetic basis of the production of volatiles in lavender could serve as a foundation for future research into lavender evolution, phytochemistry, and ecology.

The Perfect Storm: Gene Tree Estimation Error, Incomplete Lineage Sorting, and Ancient Gene Flow Explain the Most Recalcitrant Ancient Angiosperm Clade, Malpighiales

The genomic revolution offers renewed hope of resolving rapid radiations in the Tree of Life. The development of the multispecies coalescent model and improved gene tree estimation methods can better accommodate gene tree heterogeneity caused by incomplete lineage sorting (ILS) and gene tree estimation error stemming from the short internal branches. However, the relative influence of these factors in species tree inference is not well understood. Using anchored hybrid enrichment, we generated a data set including 423 single-copy loci from 64 taxa representing 39 families to infer the species tree of the flowering plant order Malpighiales. This order includes 9 of the top 10 most unstable nodes in angiosperms, which have been hypothesized to arise from the rapid radiation during the Cretaceous. Here, we show that coalescent-based methods do not resolve the backbone of Malpighiales and concatenation methods yield inconsistent estimations, providing evidence that gene tree heterogeneity is high in this clade. Despite high levels of ILS and gene tree estimation error, our simulations demonstrate that these two factors alone are insufficient to explain the lack of resolution in this order. To explore this further, we examined triplet frequencies among empirical gene trees and discovered some of them deviated significantly from those attributed to ILS and estimation error, suggesting gene flow as an additional and previously unappreciated phenomenon promoting gene tree variation in Malpighiales. Finally, we applied a novel method to quantify the relative contribution of these three primary sources of gene tree heterogeneity and demonstrated that ILS, gene tree estimation error, and gene flow contributed to 10.0$\%$, 34.8$\%$, and 21.4$\%$ of the variation, respectively. Together, our results suggest that a perfect storm of factors likely influence this lack of resolution, and further indicate that recalcitrant phylogenetic relationships like the backbone of Malpighiales may be better represented as phylogenetic networks. Thus, reducing such groups solely to existing models that adhere strictly to bifurcating trees greatly oversimplifies reality, and obscures our ability to more clearly discern the process of evolution. [Coalescent; concatenation; flanking region; hybrid enrichment, introgression; phylogenomics; rapid radiation, triplet frequency.].

Genome-Wide Duplication of Allotetraploid Brassica napus Produces Novel Characteristics and Extensive Ploidy Variation in Self-Pollinated Progeny

Whole genome duplications (WGDs) have played a major role in angiosperm species evolution. Polyploid plants have undergone multiple cycles of ancient WGD events during their evolutionary history. However, little attention has been paid to the additional WGD of the existing allopolyploids. In this study, we explored the influences of additional WGD on the allopolyploid Brassica napus Compared to tetraploid B. napus, octoploid B. napus (AAAACCCC, 2n = 8x =76) showed significant differences in phenotype, reproductive ability and the ploidy of self-pollinated progeny. Genome duplication also altered a key reproductive organ feature in B. napus, that is, increased the number of pollen apertures. Unlike autopolyploids produced from the diploid Brassica species, the octoploid B. napus produced from allotetraploid B. napus had a relatively stable meiotic process, high pollen viability and moderate fertility under self-pollination conditions, indicating that sub-genomic interactions may be important for the successful establishment of higher-order polyploids. Doubling the genome of B. napus provided us with an opportunity to gain insight into the flexibility of the Brassica genomes. The genome size of self-pollinated progeny of octoploid B. napus varied greatly, and was accompanied by extensive genomic instability, such as aneuploidy, mixed-ploidy and mitotic abnormality. The octoploid B. napus could go through any of genome reduction, equilibrium or expansion in the short-term, thus providing a novel karyotype library for the Brassica genus. Our results reveal the short-term evolutionary consequences of recurrent polyploidization events, and help to deepen our understanding of polyploid plant evolution.

Spatial, climate and ploidy factors drive genomic diversity and resilience in the widespread grass Themeda triandra

Global climate change poses a significant threat to natural communities around the world, with many plant species showing signs of climate stress. Grassland ecosystems are not an exception, with climate change compounding contemporary pressures such as habitat loss and fragmentation. In this study, we assess the climate resilience of Themeda triandra, a foundational species and the most widespread plant in Australia, by assessing the relative contributions of spatial, environmental and ploidy factors to contemporary genomic variation. Reduced-representation genome sequencing on 472 samples from 52 locations was used to test how the distribution of genomic variation, including ploidy polymorphism, supports adaptation to hotter and drier climates. We explicitly quantified isolation by distance (IBD) and isolation by environment (IBE) and predicted genomic vulnerability of populations to future climates based on expected deviation from current genomic composition. We found that a majority (54%) of genomic variation could be attributed to IBD, while an additional 22% (27% when including ploidy information) could be explained by two temperature and two precipitation climate variables demonstrating IBE. Ploidy polymorphisms were common within populations (31/52 populations), indicating that ploidy mixing is characteristic of T. triandra populations. Genomic vulnerabilities were found to be heterogeneously distributed throughout the landscape, and our analysis suggested that ploidy polymorphism, along with other factors linked to polyploidy, reduced vulnerability to future climates by 60% (0.25-0.10). Our data suggests that polyploidy may facilitate adaptation to hotter climates and highlight the importance of incorporating ploidy in adaptive management strategies to promote the resilience of this and other foundation species.