Phenotypic and genetic evidence for ecological speciation of Aquilegia japonica and A. oxysepala

Genetically- and environmentally-dependent processes drive interspecific and intraspecific divergence in the Chinese relict endemic genus Dipteronia

China is a hotspot of relict plant species that were once widespread throughout the Northern Hemisphere. Recent research has demonstrated that the occurrence of long-term stable refugia in the mountainous regions of central and south-western China allowed their persistence through the late Neogene climate fluctuations. One of these relict lineages is Dipteronia, an oligotypic tree genus with a fossil record extending to the Paleocene. Here, we investigated the genetic variability, demographic dynamics and diversification patterns of the two currently recognized Dipteronia species (D ipteronia sinensis and D . dyeriana). Molecular data were obtained from 45 populations of Dipteronia by genotyping three cpDNA regions, two single copy nuclear genes and 15 simple sequence repeat loci. The genetic study was combined with niche comparison analyses on the environmental space, ecological niche modeling, and landscape connectivity analysis. We found that the two Dipteronia species have highly diverged both in genetic and ecological terms. Despite the incipient speciation processes that can be observed in D. sinensis, the occurrence of long-term stable refugia and, particularly, a dispersal corridor along Daba Shan-west Qinling, likely ensured its genetic and ecological integrity to date. Our study will not only help us to understand how populations of Dipteronia species responded to the tectonic and climatic changes of the Cenozoic, but also provide insight into how Arcto-Tertiary relict plants in East Asia survived, evolved, and diversified.

A cornucopia of diversity-Ranunculales as a model lineage

The Ranunculales are a hyperdiverse lineage in many aspects of their phenotype, including growth habit, floral and leaf morphology, reproductive mode, and specialized metabolism. Many Ranunculales species, such as opium poppy and goldenseal, have a high medicinal value. In addition, the order includes a large number of commercially important ornamental plants, such as columbines and larkspurs. The phylogenetic position of the order with respect to monocots and core eudicots and the diversity within this lineage make the Ranunculales an excellent group for studying evolutionary processes by comparative studies. Lately, the phylogeny of Ranunculales was revised, and genetic and genomic resources were developed for many species, allowing comparative analyses at the molecular scale. Here, we review the literature on the resources for genetic manipulation and genome sequencing, the recent phylogeny reconstruction of this order, and its fossil record. Further, we explain their habitat range and delve into the diversity in their floral morphology, focusing on perianth organ identity, floral symmetry, occurrences of spurs and nectaries, sexual and pollination systems, and fruit and dehiscence types. The Ranunculales order offers a wealth of opportunities for scientific exploration across various disciplines and scales, to gain novel insights into plant biology for researchers and plant enthusiasts alike.

The past, present, and future of ecogeographic isolation between closely related Aquilegia plants

Quantifying the strength of the ecogeographic barrier is an important aspect of plant speciation research, and serves as a practical step to understanding the evolutionary trajectory of plants under climate change. Here, we quantified the extent of ecogeographic isolation in four closely related Aquilegia species that radiated in the Mountains of SW China and adjacent regions, often lacking intrinsic barriers. We used environmental niche models to predict past, present, and future species potential distributions and compared them to determine the degree of overlap and ecogeographic isolation. Our investigation found significant ecological differentiation in all studied species pairs except A. kansuensis and A. ecalacarata. The current strengths of ecogeographic isolation are above 0.5 in most cases. Compared with current climates, most species had an expanding range in the Last Glacial Maximum, the Mid Holocene, and under four future climate scenarios. Our results suggested that ecogeographic isolation contributes to the diversification and maintenance of Aquilegia species in the Mountains of northern and SW China and would act as an essential reproductive barrier in the future.

Genetic and Epigenetic Signatures Associated with the Divergence of Aquilegia Species

Widely grown in the Northern Hemisphere, the genus Aquilegia (columbine) is a model system in adaptive radiation research. While morphological variations between species have been associated with environmental factors, such as pollinators, how genetic and epigenetic factors are involved in the rapid divergence in this genus remains under investigated. In this study, we surveyed the genomes and DNA methylomes of ten Aquilegia species, representative of the Asian, European and North American lineages. Our analyses of the phylogeny and population structure revealed high genetic and DNA methylomic divergence across these three lineages. By multi-level genome-wide scanning, we identified candidate genes exhibiting lineage-specific genetic or epigenetic variation patterns that were signatures of inter-specific divergence. We demonstrated that these species-specific genetic variations and epigenetic variabilities are partially independent and are both functionally related to various biological processes vital to adaptation, including stress tolerance, cell reproduction and DNA repair. Our study provides an exploratory overview of how genetic and epigenetic signatures are associated with the diversification of the Aquilegia species.

Characterization of Terpene synthase variation in flowers of wild aquilegia species from Northeastern Asia

There are several causes for the great diversity in floral terpenes. The terpene products are determined by the catalytic fidelity, efficiency and plasticity of the active sites of terpene synthases (TPSs). However, the molecular mechanism of TPS in catalyzing terpene biosynthesis and its evolutionary fate in wild plant species remain largely unknown. In this study, the functionality of terpene synthases and their natural variants were assessed in two Northeastern Asia endemic columbine species and their natural hybrid. Synoptically, TPS7, TPS8, and TPS9 were highly expressed in these Aquilegia species from the Zuojia population. The in vitro and in vivo enzymatic assays revealed that TPS7 and TPS8 mainly produced (+)-limonene and β-sesquiphellandrene, respectively, whereas TPS9 produced pinene, similar to the major components released from Aquilegia flowers. Multiple sequence alignment of Aquilegia TPS7 and TPS8 in the Zuojia population revealed amino acid polymorphisms. Domain swapping and amino acid substitution assays demonstrated that 413A, 503I and 529D had impacts on TPS7 catalytic activity, whereas 420G, 538F and 545 L affected the ratio of β-sesquiphellandrene to β-bisabolene in TPS8. Moreover, these key polymorphic amino acid residues were found in Aquilegia species from the Changbai Mountain population. Interestingly, amino acid polymorphisms in TPSs were present in individuals with low expression levels, and nonsynonymous mutations could impact the catalytic activity or product specificity of these genes. The results of this study will shed new light on the function and evolution of TPS genes in wild plant species and are beneficial to the modification of plant fragrances.

Genetic divergence and local adaptation of Liriodendron driven by heterogeneous environments

Ecological adaptive differentiation alters both the species diversity and intraspecific genetic diversity in forests, thus affecting the stability of forest ecosystems. Therefore, knowledge of the genetic underpinnings of the ecological adaptive differentiation of forest species is critical for effective species conservation. In this study, single-nucleotide polymorphisms (SNPs) from population transcriptomes were used to investigate the spatial distribution of genetic variation in Liriodendron to assess whether environmental variables can explain genetic divergence. We examined the contributions of environmental variables to population divergence and explored the genetic underpinnings of local adaptation using a landscape genomic approach. Niche models and statistical analyses showed significant niche divergence between L. chinense and L. tulipifera, suggesting that ecological adaptation may play a crucial role in driving interspecific divergence. We detected a new fine-scale genetic structure in L. chinense, and divergence of the six groups occurred during the late Pliocene to early Pleistocene. Redundancy analysis (RDA) revealed significant associations between genetic variation and multiple environmental variables. Environmental association analyses identified 67 environmental association loci (EALs; nonsynonymous SNPs) that underwent interspecific or intraspecific differentiation, 28 of which were associated with adaptive genes. These 28 candidate adaptive loci provide substantial evidence for local adaptation in Liriodendron. Our findings reveal ecological adaptive divergence pattern between Liriodendron species and provide novel insight into the role of heterogeneous environments in shaping genetic structure and driving local adaptation among populations, informing future L. chinense conservation efforts.

Divergence in the Aquilegia ecalcarata complex is correlated with geography and climate oscillations: Evidence from plastid genome data

Quaternary climate oscillations and geographical heterogeneity play important roles in determining species and genetic diversity distribution patterns, but how these factors affect the migration and differentiation of East Asian plants species at the population level remains poorly understood. The Aquilegia ecalcarata complex, a group that originated in the Late Tertiary and is widely distributed throughout East Asia, displays high genetic variation that is suitable for studying elaborate phylogeographic patterns and demographic history related to the impact of Quaternary climate and geography. We used plastid genome data from 322 individuals in 60 populations of the A. ecalcarata complex to thoroughly explore the impact of Quaternary climate oscillations and geography on the phylogeographic patterns and demographic history of the A. ecalcarata complex through a series of phylogenetic, divergence time estimation, and demographic history analyses. The dry, cold climate and frequent climate oscillations that occurred during the early Pleistocene and the Mid-Pleistocene transition led to the differentiation of the A. ecalcarata complex, which was isolated in various areas. Geographically, the A. ecalcarata complex can be divided into Eastern and Western Clades and five subclades, which conform to the divergence of the East Asian flora. Our results clearly show the impact of Quaternary climate and geography on evolutionary history at the population level. These findings promote the understanding of the relationship between plant genetic differentiation and climate and geographical factors of East Asia at the population level.

Evolutionary history of Hemerocallis in Japan inferred from chloroplast and nuclear phylogenies and levels of interspecific gene flow

The perennial herb genus Hemerocallis (Asphodelaceae) shows four flowering types: diurnal half-day, diurnal one-day, nocturnal half-day, and nocturnal one-day flowering. These flowering types are corresponding to their main pollinators, and probably act as a primary mechanism of reproductive isolation. To examine how the four flowering types diverged, we reconstructed the phylogeny of the Japanese species of Hemerocallis using 1615 loci of nuclear genome-wide SNPs and 2078 bp sequences of four cpDNA regions. We also examined interspecific gene flows among taxa by an Isolation-with-Migration model and a population structure analysis. Our study revealed an inconsistency between chloroplast and nuclear genome phylogenies, which may have resulted from chloroplast capture. Each of the following five clusters is monophyletic and clearly separated on the nuclear genome-wide phylogenetic tree: (I) two nocturnal flowering species with lemon-yellow flowers, H. citrina (half-day flowering) and H. lilioasphodelus (one-day flowering); (II) a diurnal one-day flowering species with yellow-orange flowers, H. middendorffii; (III) a variety of a diurnal half-day flowering species with reddish orange flowers, H. fulva var. disticha; (IV) another variety of a diurnal half-day flowering species with reddish orange flowers, H. fulva var. aurantiaca, and a diurnal one-day flowering species with yellow-orange flowers, H. major; (V) a diurnal half-day flowering species with yellow-orange flowers, H. hakuunensis. The five clusters are consistent with traditional phenotype-based taxonomy (cluster I, cluster II, and clusters III-V correspond to Hemerocallis sect. Hemerocallis, Capitatae, and Fulvae, respectively). These findings could indicate that three flowering types (nocturnal flowering, diurnal one-day flowering, and diurnal half-day flowering) diverged in early evolutionary stages of Hemerocallis and subsequently a change from diurnal half-day flowering to diurnal one-day flowering occurred in a lineage of H. major. While genetic differentiation among the five clusters was well maintained, significant gene flow was detected between most pairs of taxa, suggesting that repeated hybridization played a role in the evolution of those taxa.

Comparative chloroplast genomes and phylogenetic analysis of Aquilegia

PREMISE

Aquilegia is an ideal taxon for studying the evolution of adaptive radiation. Current phylogenies of Aquilegia based on different molecular markers are inconsistent, and therefore a clear and accurate phylogeny remains uncertain. Analyzing the chloroplast genome, with its simple structure and low recombination rate, may help solve this problem.

METHODS

Next-generation sequencing data were generated or downloaded for Aquilegia species, enabling their chloroplast genomes to be assembled. The assemblies were used to estimate the genome characteristics and infer the phylogeny of Aquilegia.

RESULTS

In this study, chloroplast genome sequences were assembled for Aquilegia species distributed across Asia, North America, and Europe. Three of the genes analyzed (petG, rpl36, and atpB) were shown to be under positive selection and may be related to adaptation. The phylogenetic tree of Aquilegia showed that its member species formed two clades with high support, North American and European species, with the Asian species being paraphyletic; A. parviflora and A. amurensis clustered with the North American species, while the remaining Asian species were found in the European clade. In addition, A. oxysepala var. kansuensis should be considered as a separate species rather than a variety.

DISCUSSION

The complete chloroplast genomes of these Aquilegia species provide new insights into the reconstruction of the phylogeny of related species and contribute to the further study of this genus.

Leaf size variations in a dominant desert shrub, Reaumuria soongarica, adapted to heterogeneous environments

The climate in arid Central Asia (ACA) has changed rapidly in recent decades, but the ecological consequences of this are far from clear. To predict the impacts of climate change on ecosystem functioning, greater attention should be given to the relationships between leaf functional traits and environmental heterogeneity. As a dominant constructive shrub widely distributed in ACA, Reaumuria soongarica provided us with an ideal model to understand how leaf functional traits of desert ecosystems responded to the heterogeneous environments of ACA. Here, to determine the influences of genetic and ecological factors, we characterized species-wide variations in leaf traits among 30 wild populations of R. soongarica and 16 populations grown in a common garden. We found that the leaf length, width, and leaf length to width ratio (L/W) of the northern lineage were significantly larger than those of other genetic lineages, and principal component analysis based on the in situ environmental factors distinguished the northern lineage from the other lineages studied. With increasing latitude, leaf length, width, and L/W in the wild populations increased significantly. Leaf length and L/W were negatively correlated with altitude, and first increased and then decreased with increasing mean annual temperature (MAT) and mean annual precipitation (MAP). Stepwise regression analyses further indicated that leaf length variation was mainly affected by latitude. However, leaf width was uncorrelated with altitude, MAT, or MAP. The common garden trial showed that leaf width variation among the eastern populations was caused by both local adaptation and phenotypic plasticity. Our findings suggest that R. soongarica preferentially changes leaf length to adjust leaf size to cope with environmental change. We also reveal phenotypic evidence for ecological speciation of R. soongarica. These results will help us better understand and predict the consequences of climate change for desert ecosystem functioning.

Rapid Divergence Followed by Adaptation to Contrasting Ecological Niches of Two Closely Related Columbine Species Aquilegia japonica and A. oxysepala

Elucidating the mechanisms underlying the genetic divergence between closely related species is crucial to understanding the origin and evolution of biodiversity. The genus Aquilegia L. has undergone rapid adaptive radiation, generating about 70 well-recognized species that are specialized to distinct habitats and pollinators. In this study, to address the underlying evolutionary mechanisms that drive the genetic divergence, we analyzed the whole genomes of two ecologically isolated Aquilegia species, A. oxysepala and A. japonica as well as their putative hybrid. Our comparative genomic analyses reveal that while the two species diverged only recently and experienced recurrent gene flow, a high level of genetic divergence is observed in their nuclear genomes. In particular, candidate genomic regions that show signature of selection differ dramatically between the two species. Given that the splitting time of the two species is broadly matched with the decrease in effective population sizes, we propose that allopatric isolation together with natural selection have preceded the interspecific gene flow in the process of speciation. The observed high genetic divergence is likely an outcome of combined effects of natural selection, genetic drift and divergent sorting of ancestral polymorphisms. Our study provides a genome-wide view of how genetic divergence has evolved between closely related species.

The Aquilegia genome provides insight into adaptive radiation and reveals an extraordinarily polymorphic chromosome with a unique history

The columbine genus Aquilegia is a classic example of an adaptive radiation, involving a wide variety of pollinators and habitats. Here we present the genome assembly of A. coerulea 'Goldsmith', complemented by high-coverage sequencing data from 10 wild species covering the world-wide distribution. Our analyses reveal extensive allele sharing among species and demonstrate that introgression and selection played a role in the Aquilegia radiation. We also present the remarkable discovery that the evolutionary history of an entire chromosome differs from that of the rest of the genome - a phenomenon that we do not fully understand, but which highlights the need to consider chromosomes in an evolutionary context.

Ecological and Genetic Divergences with Gene Flow of Two Sister Species (Leucomeris decora and Nouelia insignis) Driving by Climatic Transition in Southwest China

Understanding of the processes of divergence and speciation is a major task for biodiversity researches and may offer clearer insight into mechanisms generating biological diversity. Here, we employ an integrative approach to explore genetic and ecological differentiation of Leucomeris decora and Nouelia insignis distributed allopatrically along the two sides of the biogeographic boundary 'Tanaka Line' in Southwest China. We addressed these questions using ten low-copy nuclear genes and nine plastid DNA regions sequenced among individuals sampled from 28 populations across their geographic ranges in China. Phylogenetic, coalescent-based population genetic analyses, approximate Bayesian computation (ABC) framework and ecological niche models (ENMs) were conducted. We identified a closer phylogenetic relationship in maternal lineage of L. decora with N. insignis than that between L. decora and congeneric Leucomeris spectabilis. A deep divergence between the two species was observed and occurred at the boundary between later Pliocene and early Pleistocene. However, the evidence of significant chloroplast DNA gene flow was also detected between the marginal populations of L. decora and N. insignis. Niche models and statistical analyses showed significant ecological differentiation, and two nuclear loci among the ten nuclear genes may be under divergent selection. These integrative results imply that the role of climatic shift from Pliocene to Pleistocene may be the prominent factor for the divergence of L. decora and N. insignis, and population expansion after divergence may have given rise to chloroplast DNA introgression. The divergence was maintained by differential selection despite in the face of gene flow.