Spatiotemporal reconstruction of the Aquilegia rapid radiation through next-generation sequencing of rapidly evolving cpDNA regions

Characterizing complete mitochondrial genome of Aquilegia amurensis and its evolutionary implications

BACKGROUND

Aquilegia is a model system for studying the evolution of adaptive radiation. However, very few studies have been conducted on the Aquilegia mitochondrial genome. Since mitochondria play a key role in plant adaptation to abiotic stress, analyzing the mitochondrial genome may provide a new perspective for understanding adaptive evolution.

RESULTS

The Aquilegia amurensis mitochondrial genome was characterized by a circular chromosome and two linear chromosomes, with a total length of 538,736 bp; the genes included 33 protein-coding genes, 24 transfer RNA (tRNA) genes and 3 ribosomal RNA (rRNA) genes. We subsequently conducted a phylogenetic analysis based on single nucleotide polymorphisms (SNPs) in the mitochondrial genomes of 18 Aquilegia species, which were roughly divided into two clades: the European-Asian clade and the North American clade. Moreover, the genes mttB and rpl5 were shown to be positively selected in European-Asian species, and they may help European and Asian species adapt to environmental changes.

CONCLUSIONS

In this study, we assembled and annotated the first mitochondrial genome of the adaptive evolution model plant Aquilegia. The subsequent analysis provided us with a basis for further molecular studies on Aquilegia mitochondrial genomes and valuable information on adaptive evolution in Aquilegia.

Phylogeny, biogeography, and character evolution of Anaphalis (Gnaphalieae, Asteraceae)

The HAP clade, mainly including Helichrysum Mill, Anaphalis DC., and Pseudognaphalium Kirp., is a major component of tribe Gnaphalieae (Asteraceae). In this clade, Anaphalis represents the largest genus of Asian Gnaphalieae. The intergeneric relationships among Anaphalis and its related genera and the infrageneric taxonomy of this genus are complex and remain controversial. However, there are few studies that have focused on these issues. Herein, based on the current most comprehensive sampling of the HAP clade, especially Anaphalis, we conducted phylogenetic analyses using chloroplast (cp) genome and nuclear ribosomal DNA (nrDNA) to evaluate the relationships within HAP clade, test the monophyly of Anaphalis, and examine the infrageneric taxonomy of this genus. Meanwhile, the morphological characters were verified to determine the circumscription and infrageneric taxonomy system of Anaphalis. Additionally, the biogeographical history, diversification processes, and evolution of crucial morphological characters were estimated and inferred. Our phylogenetic analyses suggested that Anaphalis is polyphyletic because it nested with Helichrysum and Pseudognaphalium. Two and four main clades of Anaphalis were identified in cp genome and nrDNA trees, respectively. Compared with nrDNA trees, the cp genome trees were more effective for phylogenetic resolution. After comprehensively analyzing morphological and phylogenetic evidence, it was concluded that the achene surface ornamentation and leaf base showed less homoplasy and supported the two Anaphalis lineages that were inferred from cp genome. Our biogeographical analyses based on cp genome indicated that HAP clade underwent rapid diversification from late Miocene to Pliocene. The two Anaphalis lineages appeared to have originated in Africa, then spread to Western and Southern Asia, and subsequently moved into Southwestern China forming a diversity center. The dispersal patterns of the two Anaphalis lineages were different. One dispersed around the world, except in Africa and South America. The other one dispersed to Eastern and Southeastern Asia from the ancestral origin region.

Flower morphology as a predictor of pollination mode in a biotic to abiotic pollination continuum

BACKGROUND AND AIMS

Wind pollination has evolved repeatedly in flowering plants, yet the identification of a wind pollination syndrome as a set of integrated floral traits can be elusive. Thalictrum (Ranunculaceae) comprises temperate perennial herbs that have transitioned repeatedly from insect to wind pollination while also exhibiting mixed pollination, providing an ideal system to test for evolutionary correlation between floral morphology and pollination mode in a biotic to abiotic continuum. Moreover, the lack of floral organ fusion across this genus allows testing for specialization to pollination vectors in the absence of this feature.

METHODS

We expanded phylogenetic sampling in the genus from a previous study using six chloroplast loci, which allowed us to test whether species cluster into distinct pollination syndromes based on floral morphology. We then used multivariate analyses on floral traits followed by ancestral state reconstruction of the emerging flower morphotypes and determined whether these traits are evolutionarily correlated under a Bayesian framework with Brownian motion.

KEY RESULTS

Floral traits fell into five distinct clusters, which were reduced to three after considering phylogenetic relatedness and were largely consistent with flower morphotypes and associated pollination vectors. Multivariate evolutionary analyses found a positive correlation between the lengths of floral reproductive structures (styles, stigmas, filaments and anthers). Shorter reproductive structures tracked insect-pollinated species and clades in the phylogeny, whereas longer structures tracked wind-pollinated ones, consistent with selective pressures exerted by biotic vs. abiotic pollination vectors, respectively.

CONCLUSIONS

Although detectable suites of integrated floral traits across Thalictrum were correlated with wind or insect pollination at the extremes of the morphospace distribution, a presumed intermediate, mixed pollination mode morphospace was also detected. Thus, our data broadly support the existence of detectable flower morphotypes from convergent evolution underlying the evolution of pollination mode in Thalictrum, presumably via different paths from an ancestral mixed pollination state.

Phylogeny of Leontopodium (Asteraceae) in China-with a reference to plastid genome and nuclear ribosomal DNA

The infrageneric taxonomy system, species delimitation, and interspecies systematic relationships of Leontopodium remain controversial and complex. However, only a few studies have focused on the molecular phylogeny of this genus. In this study, the characteristics of 43 chloroplast genomes of Leontopodium and its closely related genera were analyzed. Phylogenetic relationships were inferred based on chloroplast genomes and nuclear ribosomal DNA (nrDNA). Finally, together with the morphological characteristics, the relationships within Leontopodium were identified and discussed. The results showed that the chloroplast genomes of Filago, Gamochaeta, and Leontopodium were well-conserved in terms of gene number, gene order, and GC content. The most remarkable differences among the three genera were the length of the complete chloroplast genome, large single-copy region, small single-copy region, and inverted repeat region. In addition, the chloroplast genome structure of Leontopodium exhibited high consistency and was obviously different from that of Filago and Gamochaeta in some regions, such as matk, trnK (UUU)-rps16, petN-psbM, and trnE (UUC)-rpoB. All the phylogenetic trees indicated that Leontopodium was monophyletic. Except for the subgeneric level, our molecular phylogenetic results were inconsistent with the previous taxonomic system, which was based on morphological characteristics. Nevertheless, we found that the characteristics of the leaf base, stem types, and carpopodium base were phylogenetically correlated and may have potential value in the taxonomic study of Leontopodium. In the phylogenetic trees inferred using complete chloroplast genomes, the subgen. Leontopodium was divided into two clades (Clades 1 and 2), with most species in Clade 1 having herbaceous stems, amplexicaul, or sheathed leaves, and constricted carpopodium; most species in Clade 2 had woody stems, not amplexicaul and sheathed leaves, and not constricted carpopodium.

Phylogeny of Trigonotis in China-with a special reference to its nutlet morphology and plastid genome

The genus Trigonotis comprises nearly 60 species mainly distributed in East and Southeast Asia. China has the largest number of Trigonotis species in the world, with a total of 44 species, of which 38 are endemic. Nutlet morphology is useful for the taxonomic delimitation of Trigonotis. However, there are still controversial circumscriptions of nutlet shape in some species. In previous studies, interspecies phylogenetic relationships were inferred using few DNA markers and very few taxa, which possibly led to erroneous or incomplete conclusions. In this study, the nutlet morphology of 39 Trigonotis taxa and the characteristics of 34 complete chloroplast genomes (29 taxa) were investigated and analyzed. Then, the phylogenetic relationships were discussed within this genus based on complete chloroplast genomes. To the best of our knowledge, this study is the first comprehensive analysis of nutlet morphology and complete chloroplast genome of Trigonotis. Based on nutlet morphology, Trigonotis can be divided into two groups: Group 1, hemispherical or oblique tetrahedron with carpopodiums, and Group 2, inverted tetrahedron without carpopodiums. The chloroplast genome of Trigonotis exhibited a typical quadripartite structure, including 84-86 protein-coding, 37 transfer RNA, and 8 ribosomal RNA genes, with a total length of 147,247-148,986 bp. Genes in the junctions were well conserved in Trigonotis, similar to those in other Boraginaceae s.str. species. Furthermore, Trigonotis chloroplast genomes showed relatively high diversity, with more conserved genic regions than intergenic regions; in addition, we detected 14 hot spots (Pi > 0.005) in non-coding regions. Phylogenetic analyses based on chloroplast genome data identified highly resolved relationships between Trigonotis species. Specifically, Trigonotis was divided into two clades with strong support: one clade included species with hemispherical or oblique tetrahedron nutlets with carpopodiums and bracts, whereas the other clade included species with inverted tetrahedron nutlets without carpopodiums or bracts. Our results may inform future taxonomic, phylogenetic, and evolutionary studies on Boraginaceae.

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.

Geographic-genomic and geographic-phenotypic differentiation of the Aquilegia viridiflora complex

How species diverge into different lineages is a central issue in evolutionary biology. Despite the increasing evidence indicating that such divergences do not need geographic isolation, the correlation between lineage divergence and the adaptive ecological divergence of phenotype corresponding to distribution is still unknown. In addition, gene flow has been widely detected during and through such diverging processes. We used one widely distributed Aquilegia viridiflora complex as a model system to examine genomic differentiation and corresponding phenotypic variations along geographic gradients. Our phenotypic analyses of 20 populations from northwest to northeast China identified two phenotypic groups along the geographic cline. All examined traits are distinct from each other, although a few intermediate individuals occur in their contacting regions. We further sequenced the genomes of representative individuals of each population. However, four distinct genetic lineages were detected based on nuclear genomes. In particular, we recovered numerous genetic hybrids in the contact regions of four lineages. Gene flow is widespread and continuous between four lineages but much higher between contacting lineages than geographically isolated lineages. Gene flow and natural selection might result in inconsistency between heredity and phenotype. Moreover, many genes with fast lineage-specific mutations were identified to be involved in local adaptation. Our results suggest that both geographic isolation and local selection exerted by the environment and pollinators may together create geographic distributions of phenotypic variations as well as the underlying genomic divergences in numerous lineages.

The diversity of elaborate petals in Isopyreae (Ranunculaceae): a special focus on nectary structure

Elaborate petals are highly diverse in morphology, structure, and epidermal differentiation and play a key role in attracting pollinators. There have been few studies on the elaborate structure of petals in the tribe Isopyreae (Ranunculaceae). Seven genera in Isopyreae (Aquilegia, Semiaquilegia, Urophysa, Isopyrum, Paraquilegia, Dichocarpum, and Leptopyrum) have petals that vary in morphology, and two genera (Enemion and Thalictrum) have no petals. The petals of nine species belonged to 7 genera in the tribe were studied to reveal their nectary structure, epidermal micromorphology and ancestral traits. The petal nectaries of Isopyreae examined in this study were located at the tip of spurs (Aquilegia yabeana and A. rockii), or the bottom of shallow sacs (Semiaquilegia adoxoides, Urophysa henryi, Isopyrum manshuricum, and Paraquilegia microphylla), a cup-shaped structure (Dichocarpum fargesii) and a bilabiate structure (Leptopyrum fumarioides). The petal nectary of eight species in Isopyreae (except A. ecalcarata) was composed of secretory epidermis, nectary parenchyma, and vascular tissues, and some sieve tubes reached the secretory parenchyma cells. Among the eight species with nectaries examined in the present study, A. yabeana had the most developed nectaries, with 10-15 layers of secretory parenchyma cells. The epidermal cells of mature petals of the nine species were divided into 11 types. Among these 11 types, there were two types of secretory cells and two types of trichomes. Aquilegia yabeana and A. rockii had the highest number of cell types (eight types), and I. manshuricum and L. fumarioides had the lowest number of cell types (three types). Aquilegia ecalcarata had no secretory cells, and the papillose conical polygonal secretory cells of D. fargesii were different from those of the other seven species with nectaries. Trichomes were found only in Aquilegia, Semiaquilegia, Urophysa, and Paraquilegia. The ancestral mode of nectar presentation in Isopyreae was petals with hidden nectar (70.58%). The different modes of nectar presentation in petals may reflect adaptations to different pollinators in Isopyreae.

The Boechera model system for evolutionary ecology

Model systems in biology expand the research capacity of individuals and the community. Closely related to Arabidopsis, the genus Boechera has emerged as an important ecological model owing to the ability to integrate across molecular, functional, and eco-evolutionary approaches. Boechera species are broadly distributed in relatively undisturbed habitats predominantly in western North America and provide one of the few experimental systems for identification of ecologically important genes through genome-wide association studies and investigations of selection with plants in their native habitats. The ecologically, evolutionarily, and agriculturally important trait of apomixis (asexual reproduction via seeds) is common in the genus, and field experiments suggest that abiotic and biotic environments shape the evolution of sex. To date, population genetic studies have focused on the widespread species B. stricta, detailing population divergence and demographic history. Molecular and ecological studies show that balancing selection maintains genetic variation in ~10% of the genome, and ecological trade-offs contribute to complex trait variation for herbivore resistance, flowering phenology, and drought tolerance. Microbiome analyses have shown that host genotypes influence leaf and root microbiome composition, and the soil microbiome influences flowering phenology and natural selection. Furthermore, Boechera offers numerous opportunities for investigating biological responses to global change. In B. stricta, climate change has induced a shift of >2 weeks in the timing of first flowering since the 1970s, altered patterns of natural selection, generated maladaptation in previously locally-adapted populations, and disrupted life history trade-offs. Here we review resources and results for this eco-evolutionary model system and discuss future research directions.

Genetic architecture underlying variation in floral meristem termination in Aquilegia

Floral organs are produced by floral meristems (FMs), which harbor stem cells in their centers. Since each flower only has a finite number of organs, the stem cell activity of an FM will always terminate at a specific time point, a process termed floral meristem termination (FMT). Variation in the timing of FMT can give rise to floral morphological diversity, but how this process is fine-tuned at a developmental and evolutionary level is poorly understood. Flowers from the genus Aquilegia share identical floral organ arrangement except for stamen whorl number (SWN), making Aquilegia a well-suited system for investigation of this process: differences in SWN between species represent differences in the timing of FMT. By crossing A. canadensis and A. brevistyla, quantitative trait locus (QTL) mapping has revealed a complex genetic architecture with seven QTL. We explored potential candidate genes under each QTL and characterized novel expression patterns of select loci of interest using in situ hybridization. To our knowledge, this is the first attempt to dissect the genetic basis of how natural variation in the timing of FMT is regulated, and our results provide insight into how floral morphological diversity can be generated at the meristematic level.

Complex developmental and transcriptional dynamics underlie pollinator-driven evolutionary transitions in nectar spur morphology in Aquilegia (columbine)

PREMISE

Determining the developmental programs underlying morphological variation is key to elucidating the evolutionary processes that generated the stunning biodiversity of the angiosperms. Here, we characterized the developmental and transcriptional dynamics of the elaborate petal nectar spur of Aquilegia (columbine) in species with contrasting pollination syndromes and spur morphologies.

METHODS

We collected petal epidermal cell number and length data across four Aquilegia species, two with short, curved nectar spurs of the bee-pollination syndrome and two with long, straight spurs of the hummingbird-pollination syndrome. We also performed RNA-seq on A. brevistyla (bee) and A. canadensis (hummingbird) distal and proximal spur compartments at multiple developmental stages. Finally, we intersected these data sets with a previous QTL mapping study on spur length and shape to identify new candidate loci.

RESULTS

The differential growth between the proximal and distal surfaces of curved spurs is primarily driven by differential cell division. However, independent transitions to straight spurs in the hummingbird syndrome have evolved by increasing differential cell elongation between spur surfaces. The RNA-seq data reveal these tissues to be transcriptionally distinct and point to auxin signaling as being involved with the differential cell elongation responsible for the evolution of straight spurs. We identify several promising candidate genes for future study.

CONCLUSIONS

Our study, taken together with previous work in Aquilegia, reveals the complexity of the developmental mechanisms underlying trait variation in this system. The framework we established here will lead to exciting future work examining candidate genes and processes involved in the rapid radiation of the genus.

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.

Evidence for two types of Aquilegia ecalcarata and its implications for adaptation to new environments

Spurs have played an important role in the radiation of the genus Aquilegia, but little is known about how the spurless state arose in A. ecalcarata. Here we aim to characterize the genetic divergence within A. ecalcarata and gain insights into the origin of this species. A total of 19 populations from A. ecalcarata and 23 populations from three of its closest relatives (Aquilegia kansuensis, Aquilegia rockii and Aquilegia yabeana) were sampled in this study. We sequenced fifteen nuclear gene fragments across the genome and three chloroplast loci to conduct phylogenetic, PCoA and STRUCTURE analyses. Our analyses indicate that A. ecalcarata may not be monophyletic and can be divided into two distinct lineages (A. ecalcarata I and A. ecalcarata II). A. ecalcarata I is genetically close to A. kansuensis, whereas A. ecalcarata II is close to A. rockii. Isolation-with-migration analysis suggested that historical gene flow was low between A. ecalcarata I and A. rockii, as well as between A. ecalcarata II and A. kansuensis. The two distinct lineages of A. ecalcarata show significant divergence in 13 floral traits and also have distinct distributions. In addition, both A. ecalcarata I and II are adapted to a stony environment that differs from that of their closest relatives, indicating a habitat shift may have driven new adaptations. Our findings enrich the understanding of how floral evolution contributes to species diversification.

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.

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.

Genetic architecture of floral traits in bee- and hummingbird-pollinated sister species of Aquilegia (columbine)

Interactions with animal pollinators have helped shape the stunning diversity of flower morphologies across the angiosperms. A common evolutionary consequence of these interactions is that some flowers have converged on suites of traits, or pollination syndromes, that attract and reward specific pollinator groups. Determining the genetic basis of these floral pollination syndromes can help us understand the processes that contributed to the diversification of the angiosperms. Here, we characterize the genetic architecture of a bee-to-hummingbird pollination shift in Aquilegia (columbine) using QTL mapping of 17 floral traits encompassing color, nectar composition, and organ morphology. In this system, we find that the genetic architectures underlying differences in floral color are quite complex, and we identify several likely candidate genes involved in anthocyanin and carotenoid floral pigmentation. Most morphological and nectar traits also have complex genetic underpinnings; however, one of the key floral morphological phenotypes, nectar spur curvature, is shaped by a single locus of large effect.

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.

Insights into the evo-devo of plant reproduction using next-generation sequencing approaches

The development of plant model organisms has traditionally been analyzed using resource-heavy, tailored applications that are not easily transferable to distantly related non-model taxa. Thus, our understanding of plant development has been limited to a subset of traits, and evolutionary studies conducted most effectively either across very wide [e.g. Arabidopsis thaliana and Oryza sativa (rice)] or narrow (i.e. population level) phylogenetic distances. As plant biologists seek to capitalize on natural diversity for crop improvement, enhance ecosystem functioning, and better understand plant responses to climate change, high-throughput and broadly applicable forms of existing molecular biology assays are becoming an invaluable resource. Next-generation sequencing (NGS) is increasingly becoming a powerful tool in evolutionary developmental biology (evo-devo) studies, particularly through its application to understanding trait evolution at different levels of gene regulation. Here, I review some of the most common and emerging NGS-based methods, using exemplar studies in reproductive plant evo-devo to illustrate their potential.

A pilot study of eDNA metabarcoding to estimate plant biodiversity by an alpine glacier core (Adamello glacier, North Italy)

Current biodiversity loss is a major concern and thus biodiversity assessment of modern ecosystems is compelling and needs to be contextualized on a longer timescale. High Throughput Sequencing (HTS) is progressively becoming a major source of data on biodiversity time series. In this multi proxy study, we tested, for the first time, the potential of HTS to estimate plant biodiversity archived in the surface layers of a temperate alpine glacier, amplifying the trnL barcode for vascular plants from eDNA of firn samples. A 573 cm long core was drilled by the Adamello glacier and cut into sections; produced samples were analyzed for physical properties, stable isotope ratio, and plant biodiversity by eDNA metabarcoding and conventional light microscopy analysis. Results highlighted the presence of pollen and plant remains within the distinct layers of snow, firn and ice. While stable isotope ratio showed a scarcely informative pattern, DNA metabarcoding described distinct plant species composition among the different samples, with a broad taxonomic representation of the biodiversity of the catchment area and a high-ranking resolution. New knowledge on climate and plant biodiversity changes of large catchment areas can be obtained by this novel approach, relevant for future estimates of climate change effects.

Taxonomy of Aquilegia glandulosa Fischer ex Link and related species (Ranunculaceae) in North Asia

The paper presents the annotated checklist and morphological differences of the Aquilegia glandulosa from the nearest related species and provides new data on the distribution of all species from this group. An identification key to the species of Aquilegia glandulosa group from the North Asia is provided and diagnostic characters are discussed.

POPOVICH, encoding a C2H2 zinc-finger transcription factor, plays a central role in the development of a key innovation, floral nectar spurs, in Aquilegia

The evolution of novel features, such as eyes or wings, that allow organisms to exploit their environment in new ways can lead to increased diversification rates. Therefore, understanding the genetic and developmental mechanisms involved in the origin of these key innovations has long been of interest to evolutionary biologists. In flowering plants, floral nectar spurs are a prime example of a key innovation, with the independent evolution of spurs associated with increased diversification rates in multiple angiosperm lineages due to their ability to promote reproductive isolation via pollinator specialization. As none of the traditional plant model taxa have nectar spurs, little is known about the genetic and developmental basis of this trait. Nectar spurs are a defining feature of the columbine genus Aquilegia (Ranunculaceae), a lineage that has experienced a relatively recent and rapid radiation. We use a combination of genetic mapping, gene expression analyses, and functional assays to identify a gene crucial for nectar spur development, POPOVICH (POP), which encodes a C2H2 zinc-finger transcription factor. POP plays a central role in regulating cell proliferation in the Aquilegia petal during the early phase (phase I) of spur development and also appears to be necessary for the subsequent development of nectaries. The identification of POP opens up numerous avenues for continued scientific exploration, including further elucidating of the genetic pathway of which it is a part, determining its role in the initial evolution of the Aquilegia nectar spur, and examining its potential role in the subsequent evolution of diverse spur morphologies across the genus.

A role for the Auxin Response Factors ARF6 and ARF8 homologs in petal spur elongation and nectary maturation in Aquilegia

The petal spur of the basal eudicot Aquilegia is a key innovation associated with the adaptive radiation of the genus. Previous studies have shown that diversification of Aquilegia spur length can be predominantly attributed to variation in cell elongation. However, the genetic pathways that control the development of petal spurs are still being investigated. Here, we focus on a pair of closely related homologs of the AUXIN RESPONSE FACTOR family, AqARF6 and AqARF8, to explore their roles in Aquileiga coerulea petal spur development. Expression analyses of the two genes show that they are broadly expressed in vegetative and floral organs, but have relatively higher expression in petal spurs, particularly at later stages. Knockdown of the two AqARF6 and AqARF8 transcripts using virus-induced gene silencing resulted in largely petal-specific defects, including a significant reduction in spur length due to a decrease in cell elongation. These spurs also exhibited an absence of nectar production, which was correlated with downregulation of STYLISH homologs that have previously been shown to control nectary development. This study provides the first evidence of ARF6/8 homolog-mediated petal development outside the core eudicots. The genes appear to be specifically required for cell elongation and nectary maturation in the Aquilegia petal spur.

A chromosome-scale reference genome of Aquilegia oxysepala var. kansuensis

The genus Aquilegia (Ranunculaceae) has been cultivated as ornamental and medicinal plants for centuries. With petal spurs of strikingly diverse size and shape, Aquilegia has also been recognized as an excellent system for evolutionary studies. Pollinator-mediated selection for longer spurs is believed to have shaped the evolution of this genus, especially the North American taxa. Recently, however, an opposite evolutionary trend was reported in an Asian lineage, where multiple origins of mini- or even nonspurred morphs have occurred. Interesting as it is, the lack of genomic resources has limited our ability to decipher the molecular and evolutionary mechanisms underlying spur reduction in this special lineage. Using long-read sequencing (PacBio Sequel), in combination with optical maps (BioNano DLS) and Hi-C, we assembled a high-quality reference genome of A. oxysepala var. kansuensis, a sister species to the nonspurred taxon. The final assembly is approximately 293.2 Mb, 94.6% (277.4 Mb) of which has been anchored to 7 pseudochromosomes. A total of 25,571 protein-coding genes were predicted, with 97.2% being functionally annotated. When comparing this genome with that of A. coerulea, we detected a large rearrangement between Chr1 and Chr4, which might have caused the Chr4 of A. oxysepala var. kansuensis to partly deviate from the "decaying" path that was taken before the split of Aquilegia and Semiaquilegia. This high-quality reference genome is fundamental to further investigations on the development and evolution of petal spurs and provides a strong foundation for the breeding of new horticultural Aquilegia cultivars.

Maximize Resolution or Minimize Error? Using Genotyping-By-Sequencing to Investigate the Recent Diversification of Helianthemum (Cistaceae)

A robust phylogenetic framework, in terms of extensive geographical and taxonomic sampling, well-resolved species relationships and high certainty of tree topologies and branch length estimations, is critical in the study of macroevolutionary patterns. Whereas Sanger sequencing-based methods usually recover insufficient phylogenetic signal, especially in recently diversified lineages, reduced-representation sequencing methods tend to provide well-supported phylogenetic relationships, but usually entail remarkable bioinformatic challenges due to the inherent trade-off between the number of SNPs and the magnitude of associated error rates. The genus Helianthemum (Cistaceae) is a species-rich and taxonomically complex Palearctic group of plants that diversified mainly since the Upper Miocene. It is a challenging case study since previous attempts using Sanger sequencing were unable to resolve the intrageneric phylogenetic relationships. Aiming to obtain a robust phylogenetic reconstruction based on genotyping-by-sequencing (GBS), we established a rigorous methodological workflow in which we i) explored how variable settings during dataset assembly have an impact on error rates and on the degree of resolution under concatenation and coalescent approaches, ii) assessed the effect of two extreme parameter configurations (minimizing error rates vs. maximizing phylogenetic resolution) on tree topology and branch lengths, and iii) evaluated the effects of these two configurations on estimates of divergence times and diversification rates. Our analyses produced highly supported topologically congruent phylogenetic trees for both configurations. However, minimizing error rates did produce more reliable branch lengths, critically affecting the accuracy of downstream analyses (i.e. divergence times and diversification rates). In addition to recommending a revision of intrageneric systematics, our results enabled us to identify three highly diversified lineages in Helianthemum in contrasting geographical areas and ecological conditions, which started radiating in the Upper Miocene.

Cell number explains the intraspecific spur-length variation in an Aquilegia species

Variations of nectar spur length allow pollinators to utilize resources in novel ways, leading to the different selective pressures on spurs and allowing taxa to diversify. However, the mechanisms underlying spur length variation remain unclear. Interspecific comparisons of spur length suggest that both cell division and anisotropic expansion could explain the changes of spur length, and that hormone-related genes contribute to the process of spur formation. In contrast, little is known about intraspecific spur length variation. In Aquilegia rockii, spur length varies strikingly, ranging from 1 mm to 18 mm. To examine the potential mechanisms underlying spur length variation in A. rockii, we observed cell morphology and analyzed RNA-seq of short- and long-spurred flowers. Scanning electron microscopy revealed that at two positions on spurs there were no differences in either cell density or cell anisotropy between short- and long-spurred flowers, suggesting that in A. rockii changes in cell number may explain variations in spur length. In addition, we screened transcriptomes of short- and long-spurred flowers for differentially expressed genes; this screen identified several genes linked to cell division (e.g., F-box, CDKB2-2, and LST8), a finding which is consistent with our analysis of the cellular morphology of spurs. However, we did not find any highly expressed genes involved in the hormone pathway in long-spurred flowers. In contrast to previous hypotheses that anisotropic cell expansion leads to interspecific spur variation in Aquilegia, our results suggest that cell number changes and related genes are mainly responsible for spur length variations of A. rockii. Furthermore, the underlying mechanisms of similar floral traits in morphology may be quite different, enriching our understanding of the mechanisms of flower diversity in angiosperms.

Comparative transcriptomics of early petal development across four diverse species of Aquilegia reveal few genes consistently associated with nectar spur development

BACKGROUND

Petal nectar spurs, which facilitate pollination through animal attraction and pollen placement, represent a key innovation promoting diversification in the genus Aquilegia (Ranunculaceae). Identifying the genetic components that contribute to the development of these three-dimensional structures will inform our understanding of the number and types of genetic changes that are involved in the evolution of novel traits. In a prior study, gene expression between two regions of developing petals, the laminar blade and the spur cup, was compared at two developmental stages in the horticultural variety A. coerulea 'Origami'. Several hundred genes were differentially expressed (DE) between the blade and spur at both developmental stages. In order to narrow in on a set of genes crucial to early spur formation, the current study uses RNA sequencing (RNAseq) to conduct comparative expression analyses of petals from five developmental stages between four Aquilegia species, three with morphologically variable nectar spurs, A. sibirica, A. formosa, and A. chrysantha, and one that lacks nectar spurs, A. ecalcarata.

RESULTS

Petal morphology differed increasingly between taxa across the developmental stages assessed, with petals from all four taxa being indistinguishable pre-spur formation at developmental stage 1 (DS1) and highly differentiated by developmental stage 5 (DS5). In all four taxa, genes involved in mitosis were down-regulated over the course of the assessed developmental stages, however, many genes involved in mitotic processes remained expressed at higher levels later in development in the spurred taxa. A total of 690 genes were identified that were consistently DE between the spurred taxa and A. ecalcarata at all five developmental stages. By comparing these genes with those identified as DE between spur and blade tissue in A. coerulea 'Origami', a set of only 35 genes was identified that shows consistent DE between petal samples containing spur tissue versus those without spur tissue.

CONCLUSIONS

The results of this study suggest that expression differences in very few loci are associated with the presence and absence of spurs. In general, it appears that the spurless petals of A. ecalcarata cease cell divisions and enter the cell differentiation phase at an earlier developmental time point than those that produce spurs. This much more tractable list of 35 candidates genes will greatly facilitate targeted functional studies to assess the genetic control and evolution of petal spurs in Aquilegia.

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.

Macroevolutionary dynamics of nectar spurs, a key evolutionary innovation

Floral nectar spurs are widely considered a key innovation promoting diversification in angiosperms by means of pollinator shifts. We investigated the macroevolutionary dynamics of nectar spurs in the tribe Antirrhineae (Plantaginaceae), which contains 29 genera and 300-400 species (70-80% spurred). The effect of nectar spurs on diversification was tested, with special focus on Linaria, the genus with the highest number of species. We generated the most comprehensive phylogeny of Antirrhineae to date and reconstructed the evolution of nectar spurs. Diversification rate heterogeneity was investigated using trait-dependent and trait-independent methods, and accounting for taxonomic uncertainty. The association between changes in spur length and speciation was examined within Linaria using model testing and ancestral state reconstructions. We inferred four independent acquisitions of nectar spurs. Diversification analyses revealed that nectar spurs are loosely associated with increased diversification rates. Detected rate shifts were delayed by 5-15 Myr with respect to the acquisition of the trait. Active evolution of spur length, fitting a speciational model, was inferred in Linaria, which is consistent with a scenario of pollinator shifts driving diversification. Nectar spurs played a role in diversification of the Antirrhineae, but diversification dynamics can only be fully explained by the complex interaction of multiple biotic and abiotic factors.

On the origin of orphan hybrids between Aquilegia formosa and Aquilegia flavescens

We report the investigation of an Aquilegia flavescens × A. formosa population in British Columbia that is disjunct from its parents-the latter species is present locally but ecologically separated, while the former is entirely absent. To confirm hybridity, we used multivariate analysis of floral characters of field-sampled populations to ordinate phenotypes of putative hybrids in relation to those of the parental species. Microsatellite genotypes at 11 loci from 72 parental-type and putative hybrid individuals were analysed to assess evidence for admixture. Maternally inherited plastid sequences were analysed to infer the direction of hybridization and test hypotheses on the origin of the orphan hybrid population. Plants from the orphan hybrid population are on average intermediate between typical A. formosa and A. flavescens for most phenotypes examined and show evidence of genetic admixture. This population lies beyond the range of A. flavescens, but within the range of A. formosa. No pure A. flavescens individuals were observed in the vicinity, nor is this species known to occur within 200 km of the site. The hybrids share a plastid haplotype with local A. formosa populations. Alternative explanations for this pattern are evaluated. While we cannot rule out long-distance pollen dispersal followed by proliferation of hybrid genotypes, we consider the spread of an A. formosa plastid during genetic swamping of a historical A. flavescens population to be more parsimonious.

Morphological Innovations and Vast Extensions of Mountain Habitats Triggered Rapid Diversification Within the Species-Rich Irano-Turanian Genus Acantholimon (Plumbaginaceae)

The Irano-Turanian floristic region spans a topographically complex and climatically continental territory, which has served as a source of xerophytic taxa for neighboring regions and is represented by a high percent of endemics. Yet, a comprehensive picture of the abiotic and biotic factors that have driven diversification within this biota remains to be established due to the scarcity of phylogenetic studies. Acantholimon is an important component of the subalpine steppe flora of the Irano-Turanian region, containing c. 200 cushion-forming sub-shrubby pungent-leaved species. Our recent molecular phylogenetic study has led to enlarging the circumscription of this genus to include eight mono- or oligospecific genera lacking the characteristic life-form and leaves. Using the same molecular phylogeny, here we investigate the tempo and mode of diversification as well as the biogeographic patterns in this genus, to test the hypothesis that a combination of key morphological innovations and abiotic factors is behind Acantholimon high species diversity. Molecular dating analysis indicates that Acantholimon s.l. started to diversify between the Late Miocene and the Pliocene and the biogeographic analysis points to an Eastern Iran-Afghanistan origin. Macroevolutionary models support the hypothesis that the high diversity of the genus is explained by accelerated diversification rates in two clades associated with the appearance of morphological key innovations such as a cushion life-form and pungent leaves; this would have favored the colonization of water-stressed, substrate-poor mountainous habitats along the newly uplifted IT mountains during the Mio-Pliocene. Given the apparent similarity of mountain habitats for most species of Acantholimon, we hypothesize that its current high species diversity responds to a scenario of non-adaptive radiation fueled by allopatric speciation rather than evolutionary radiation driven by ecological opportunity. Similar scenarios might underlie the high diversity of other speciose genera in the topographically complex Irano-Turanian landscape, though this remains to be tested with fine-grained distribution and climatic data.

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.

Comparative Analysis of the Chloroplast Genomes of the Chinese Endemic Genus Urophysa and Their Contribution to Chloroplast Phylogeny and Adaptive Evolution

Urophysa is a Chinese endemic genus comprising two species, Urophysa rockii and Urophysa henryi. In this study, we sequenced the complete chloroplast (cp) genomes of these two species and of their relative Semiquilegia adoxoides. Illumina sequencing technology was used to compare sequences, elucidate the intra- and interspecies variations, and infer the phylogeny relationship with other Ranunculaceae family species. A typical quadripartite structure was detected, with a genome size from 158,473 to 158,512 bp, consisting of a pair of inverted repeats separated by a small single-copy region and a large single-copy region. We analyzed the nucleotide diversity and repeated sequences components and conducted a positive selection analysis by the codon-based substitution on single-copy coding sequence (CDS). Seven regions were found to possess relatively high nucleotide diversity, and numerous variable repeats and simple sequence repeats (SSR) markers were detected. Six single-copy genes (atpA, rpl20, psaA, atpB, ndhI, and rbcL) resulted to have high posterior probabilities of codon sites in the positive selection analysis, which means that the six genes may be under a great selection pressure. The visualization results of the six genes showed that the amino acid properties across each column of all species are variable in different genera. All these regions with high nucleotide diversity, abundant repeats, and under positive selection will provide potential plastid markers for further taxonomic, phylogenetic, and population genetics studies in Urophysa and its relatives. Phylogenetic analyses based on the 79 single-copy genes, the whole complete genome sequences, and all CDS sequences showed same topologies with high support, and U. rockii was closely clustered with U. henryi within the Urophysa genus, with S. adoxoides as their closest relative. Therefore, the complete cp genomes in Urophysa species provide interesting insights and valuable information that can be used to identify related species and reconstruct their phylogeny.

Plastome sequences and exploration of tree-space help to resolve the phylogeny of riceflowers (Thymelaeaceae: Pimelea)

Data sets comprising small numbers of genetic markers are not always able to resolve phylogenetic relationships. This has frequently been the case in molecular systematic studies of plants, with many analyses being based on sequence data from only two or three chloroplast genes. An example of this comes from the riceflowers Pimelea Banks & Sol. ex Gaertn. (Thymelaeaceae), a large genus of flowering plants predominantly distributed in Australia. Despite the considerable morphological variation in the genus, low sequence divergence in chloroplast markers has led to the phylogeny of Pimelea remaining largely uncertain. In this study, we resolve the backbone of the phylogeny of Pimelea in comprehensive Bayesian and maximum-likelihood analyses of plastome sequences from 41 taxa. However, some relationships received only moderate to poor support, and the Pimelea clade contained extremely short internal branches. By using topology-clustering analyses, we demonstrate that conflicting phylogenetic signals can be found across the trees estimated from individual chloroplast protein-coding genes. A relaxed-clock dating analysis reveals that Pimelea arose in the mid-Miocene, with most divergences within the genus occurring during a subsequent rapid diversification. Our new phylogenetic estimate offers better resolution and is more strongly supported than previous estimates, providing a platform for future taxonomic revisions of both Pimelea and the broader subfamily. Our study has demonstrated the substantial improvements in phylogenetic resolution that can be achieved using plastome-scale data sets in plant molecular systematics.

Evolutionary dynamism in bryophytes: Phylogenomic inferences confirm rapid radiation in the moss family Funariaceae

Rapid diversifications of plants are primarily documented and studied in angiosperms, which are perceived as evolutionarily dynamic. Recent studies have, however, revealed that bryophytes have also undergone periods of rapid radiation. The speciose family Funariaceae, including the model taxon Physcomitrella patens, is one such lineage. Here, we infer relationships among major lineages within the Entosthodon-Physcomitrium complex from virtually complete organellar exomes (i.e., 123 genes) obtained through high throughput sequencing of genomic libraries enriched in these loci via targeted locus capture. Based on these extensive exonic data we (1) reconstructed a robust backbone topology of the Funariaceae, (2) confirmed the monophyly of Funaria and the polyphyly of Entosthodon, Physcomitrella, and Physcomitrium, and (3) argue for the occurrence of a rapid radiation within the Entosthodon-Physcomitrium complex that began 28 mya and gave rise more than half of the species diversity of the family. This diversification may have been triggered by a whole genome duplication and coincides with global Eocene cooling that continued through the Oligocene and Miocene. The Funariaceae join a growing list of bryophyte lineages whose history is marked by at least one burst of diversification, and our study thereby strengthens the view that bryophytes are evolutionarily dynamic lineages and that patterns and processes characterizing the evolution of angiosperms may be universal among land plants.

Next-generation polyploid phylogenetics: rapid resolution of hybrid polyploid complexes using PacBio single-molecule sequencing

Difficulties in generating nuclear data for polyploids have impeded phylogenetic study of these groups. We describe a high-throughput protocol and an associated bioinformatics pipeline (Pipeline for Untangling Reticulate Complexes (Purc)) that is able to generate these data quickly and conveniently, and demonstrate its efficacy on accessions from the fern family Cystopteridaceae. We conclude with a demonstration of the downstream utility of these data by inferring a multi-labeled species tree for a subset of our accessions. We amplified four c. 1-kb-long nuclear loci and sequenced them in a parallel-tagged amplicon sequencing approach using the PacBio platform. Purc infers the final sequences from the raw reads via an iterative approach that corrects PCR and sequencing errors and removes PCR-mediated recombinant sequences (chimeras). We generated data for all gene copies (homeologs, paralogs, and segregating alleles) present in each of three sets of 50 mostly polyploid accessions, for four loci, in three PacBio runs (one run per set). From the raw sequencing reads, Purc was able to accurately infer the underlying sequences. This approach makes it easy and economical to study the phylogenetics of polyploids, and, in conjunction with recent analytical advances, facilitates investigation of broad patterns of polyploid evolution.

High-throughput multiplex cpDNA resequencing clarifies the genetic diversity and genetic relationships among Brassica napus, Brassica rapa and Brassica oleracea

Brassica napus (rapeseed) is a recent allotetraploid plant and the second most important oilseed crop worldwide. The origin of B. napus and the genetic relationships with its diploid ancestor species remain largely unresolved. Here, chloroplast DNA (cpDNA) from 488 B. napus accessions of global origin, 139 B. rapa accessions and 49 B. oleracea accessions were populationally resequenced using Illumina Solexa sequencing technologies. The intraspecific cpDNA variants and their allelic frequencies were called genomewide and further validated via EcoTILLING analyses of the rpo region. The cpDNA of the current global B. napus population comprises more than 400 variants (SNPs and short InDels) and maintains one predominant haplotype (Bncp1). Whole-genome resequencing of the cpDNA of Bncp1 haplotype eliminated its direct inheritance from any accession of the B. rapa or B. oleracea species. The distribution of the polymorphism information content (PIC) values for each variant demonstrated that B. napus has much lower cpDNA diversity than B. rapa; however, a vast majority of the wild and cultivated B. oleracea specimens appeared to share one same distinct cpDNA haplotype, in contrast to its wild C-genome relatives. This finding suggests that the cpDNA of the three Brassica species is well differentiated. The predominant B. napus cpDNA haplotype may have originated from uninvestigated relatives or from interactions between cpDNA mutations and natural/artificial selection during speciation and evolution. These exhaustive data on variation in cpDNA would provide fundamental data for research on cpDNA and chloroplasts.

Using targeted enrichment of nuclear genes to increase phylogenetic resolution in the neotropical rain forest genus Inga (Leguminosae: Mimosoideae)

Evolutionary radiations are prominent and pervasive across many plant lineages in diverse geographical and ecological settings; in neotropical rainforests there is growing evidence suggesting that a significant fraction of species richness is the result of recent radiations. Understanding the evolutionary trajectories and mechanisms underlying these radiations demands much greater phylogenetic resolution than is currently available for these groups. The neotropical tree genus Inga (Leguminosae) is a good example, with ~300 extant species and a crown age of 2-10 MY, yet over 6 kb of plastid and nuclear DNA sequence data gives only poor phylogenetic resolution among species. Here we explore the use of larger-scale nuclear gene data obtained though targeted enrichment to increase phylogenetic resolution within Inga. Transcriptome data from three Inga species were used to select 264 nuclear loci for targeted enrichment and sequencing. Following quality control to remove probable paralogs from these sequence data, the final dataset comprised 259,313 bases from 194 loci for 24 accessions representing 22 Inga species and an outgroup (Zygia). Bayesian phylogenies reconstructed using either all loci concatenated or a gene-tree/species-tree approach yielded highly resolved phylogenies. We used coalescent approaches to show that the same targeted enrichment data also have significant power to discriminate among alternative within-species population histories within the widespread species I. umbellifera. In either application, targeted enrichment simplifies the informatics challenge of identifying orthologous loci associated with de novo genome sequencing. We conclude that targeted enrichment provides the large volumes of phylogenetically-informative sequence data required to resolve relationships within recent plant species radiations, both at the species level and for within-species phylogeographic studies.

Confluence, synnovation, and depauperons in plant diversification

We review the empirical phylogenetic literature on plant diversification, highlighting challenges in separating the effects of speciation and extinction, in specifying diversification mechanisms, and in making convincing arguments. In recent discussions of context dependence, key opportunities and landscapes, and indirect effects and lag times, we see a distinct shift away from single-point/single-cause 'key innovation' hypotheses toward more nuanced explanations involving multiple interacting causal agents assembled step-wise through a tree. To help crystalize this emerging perspective we introduce the term 'synnovation' (a hybrid of 'synergy' and 'innovation') for an interacting combination of traits with a particular consequence ('key synnovation' in the case of increased diversification rate), and the term 'confluence' for the sequential coming together of a set of traits (innovations and synnovations), environmental changes, and geographic movements along the branches of a phylogenetic tree. We illustrate these concepts using the radiation of Bromeliaceae. We also highlight the generality of these ideas by considering how rate heterogeneity associated with a confluence relates to the existence of particularly species-poor lineages, or 'depauperons.' Many challenges are posed by this re-purposed research framework, including difficulties associated with partial taxon sampling, uncertainty in divergence time estimation, and extinction.

Demographic History, Population Structure, and Local Adaptation in Alpine Populations of Cardamine impatiens and Cardamine resedifolia

Species evolution depends on numerous and distinct forces, including demography and natural selection. For example, local adaptation and population structure affect the evolutionary history of species living along environmental clines. This is particularly relevant in plants, which are often characterized by limited dispersal ability and the need to respond to abiotic and biotic stress factors specific to the local environment. Here we study the demographic history and the possible existence of local adaptation in two related species of Brassicaceae, Cardamine impatiens and Cardamine resedifolia, which occupy separate habitats along the elevation gradient. Previous genome-wide analyses revealed the occurrence of distinct selective pressures in the two species, with genes involved in cold response evolving particularly fast in C. resedifolia. In this study we surveyed patterns of molecular evolution and genetic variability in a set of 19 genes, including neutral and candidate genes involved in cold response, across 10 populations each of C. resedifolia and C. impatiens from the Italian Alps (Trentino). We inferred the population structure and demographic history of the two species, and tested the occurrence of signatures of local adaptation in these genes. The results indicate that, despite a slightly higher population differentiation in C. resedifolia than in C. impatiens, both species are only weakly structured and that populations sampled at high altitude experience less gene flow than low-altitude ones. None of the genes showed signatures of positive selection, suggesting that they do not seem to play relevant roles in the current evolutionary processes of adaptation to alpine environments of these species.

Molecular basis for three-dimensional elaboration of the Aquilegia petal spur

By enforcing specific pollinator interactions, Aquilegia petal nectar spurs maintain reproductive isolation between species. Spur development is the result of three-dimensional elaboration from a comparatively two-dimensional primordium. Initiated by localized, oriented cell divisions surrounding the incipient nectary, this process creates a pouch that is extended by anisotropic cell elongation. We hypothesized that the development of this evolutionary novelty could be promoted by non-mutually exclusive factors, including (i) prolonged, KNOX-dependent cell fate indeterminacy, (ii) localized organ sculpting and/or (iii) redeployment of hormone-signalling modules. Using cell division markers to guide transcriptome analysis of microdissected spur tissue, we present candidate mechanisms underlying spur outgrowth. We see dynamic expression of factors controlling cell proliferation and hormone signalling, but no evidence of contribution from indeterminacy factors. Transcriptome dynamics point to a novel recruitment event in which auxin-related factors that normally function at the organ margin were co-opted to this central structure. Functional perturbation of the transition between cell division and expansion reveals an unexpected asymmetric component of spur development. These findings indicate that the production of this three-dimensional form is an example of organ sculpting via localized cell division with novel contributions from hormone signalling, rather than a product of prolonged indeterminacy.

Complete plastome sequence of Thalictrum coreanum (Ranunculaceae) and transfer of the rpl32 gene to the nucleus in the ancestor of the subfamily Thalictroideae

BACKGROUND

Plastids originated from cyanobacteria and the majority of the ancestral genes were lost or functionally transferred to the nucleus after endosymbiosis. Comparative genomic investigations have shown that gene transfer from plastids to the nucleus is an ongoing evolutionary process but molecular evidence for recent functional gene transfers among seed plants have only been documented for the four genes accD, infA, rpl22, and rpl32.

RESULTS

The complete plastid genome of Thalictrum coreanum, the first from the subfamily Thalictroideae (Ranunculaceae), was sequenced and revealed the losses of two genes, infA and rpl32. The functional transfer of these two genes to the nucleus in Thalictrum was verified by examination of nuclear transcriptomes. A survey of the phylogenetic distribution of the rpl32 loss was performed using 17 species of Thalictrum and representatives of related genera in the subfamily Thalictroideae. The plastid-encoded rpl32 gene is likely nonfunctional in members of the subfamily Thalictroideae (Aquilegia, Enemion, Isopyrum, Leptopyrum, Paraquilegia, and Semiaquilegia) including 17 Thalictrum species due to the presence of indels that disrupt the reading frame. A nuclear-encoded rpl32 with high sequence identity was identified in both Thalictrum and Aquilegia. The phylogenetic distribution of this gene loss/transfer and the high level of sequence similarity in transit peptides suggest a single transfer of the plastid-encoded rpl32 to the nucleus in the ancestor of the subfamily Thalictroideae approximately 20-32 Mya.

CONCLUSIONS

The genome sequence of Thalictrum coreanum provides valuable information for improving the understanding of the evolution of plastid genomes within Ranunculaceae and across angiosperms. Thalictrum is unusual among the three sequenced Ranunculaceae plastid genomes in the loss of two genes infA and rpl32, which have been functionally transferred to the nucleus. In the case of rpl32 this represents the third documented independent transfer from the plastid to the nucleus with the other two transfers occurring in the unrelated angiosperm families Rhizophoraceae and Salicaceae. Furthermore, the transfer of rpl32 provides additional molecular evidence for the monophyly of the subfamily Thalictroideae.

Using targeted enrichment of nuclear genes to increase phylogenetic resolution in the neotropical rain forest genus Inga (Leguminosae: Mimosoideae)

Evolutionary radiations are prominent and pervasive across many plant lineages in diverse geographical and ecological settings; in neotropical rainforests there is growing evidence suggesting that a significant fraction of species richness is the result of recent radiations. Understanding the evolutionary trajectories and mechanisms underlying these radiations demands much greater phylogenetic resolution than is currently available for these groups. The neotropical tree genus Inga (Leguminosae) is a good example, with ~300 extant species and a crown age of 2-10 MY, yet over 6 kb of plastid and nuclear DNA sequence data gives only poor phylogenetic resolution among species. Here we explore the use of larger-scale nuclear gene data obtained though targeted enrichment to increase phylogenetic resolution within Inga. Transcriptome data from three Inga species were used to select 264 nuclear loci for targeted enrichment and sequencing. Following quality control to remove probable paralogs from these sequence data, the final dataset comprised 259,313 bases from 194 loci for 24 accessions representing 22 Inga species and an outgroup (Zygia). Bayesian phylogenies reconstructed using either all loci concatenated or a gene-tree/species-tree approach yielded highly resolved phylogenies. We used coalescent approaches to show that the same targeted enrichment data also have significant power to discriminate among alternative within-species population histories within the widespread species I. umbellifera. In either application, targeted enrichment simplifies the informatics challenge of identifying orthologous loci associated with de novo genome sequencing. We conclude that targeted enrichment provides the large volumes of phylogenetically-informative sequence data required to resolve relationships within recent plant species radiations, both at the species level and for within-species phylogeographic studies.

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

Natural selection is thought to be a driving force that can cause the evolution of reproductive isolation. The genus Aquilegia is a model system to address how natural selection promotes the process of speciation. Morphological differences between A. oxysepala, A. japonica and their hybrids were quantified for two vegetative (plant height and leaf area) and three floral morphological (sepal area, corolla length and diameter) traits. We also evaluated the genetic variability of the two species and their hybrids based on two chloroplast (1225 bp), four nuclear (5811 bp) genes and 15 microsatellites. Our results revealed that differentiation of A. japonica and A. oxysepala at the ecological and morphological levels also involved divergence at the genetic level. In addition, the analysis of nucleotide variation patterns showed that the two species possessed numerous fixation sites at nuclear genes gAA4, gA7 and gAA12. Furthermore, we found that all of the phenotypic hybrids also showed a genetically admixed ancestry. These findings suggest that natural selection has indeed facilitated the formation of distinct genetic variation patterns in the two Aquilegia species and habitat adaptation has been driving the ecologically based evolution of reproductive isolation.

Variability among the most rapidly evolving plastid genomic regions is lineage-specific: implications of pairwise genome comparisons in Pyrus (Rosaceae) and other angiosperms for marker choice

Plastid genomes exhibit different levels of variability in their sequences, depending on the respective kinds of genomic regions. Genes are usually more conserved while noncoding introns and spacers evolve at a faster pace. While a set of about thirty maximum variable noncoding genomic regions has been suggested to provide universally promising phylogenetic markers throughout angiosperms, applications often require several regions to be sequenced for many individuals. Our project aims to illuminate evolutionary relationships and species-limits in the genus Pyrus (Rosaceae)—a typical case with very low genetic distances between taxa. In this study, we have sequenced the plastid genome of Pyrus spinosa and aligned it to the already available P. pyrifolia sequence. The overall p-distance of the two Pyrus genomes was 0.00145. The intergenic spacers between ndhC–trnV, trnR–atpA, ndhF–rpl32, psbM–trnD, and trnQ–rps16 were the most variable regions, also comprising the highest total numbers of substitutions, indels and inversions (potentially informative characters). Our comparative analysis of further plastid genome pairs with similar low p-distances from Oenothera (representing another rosid), Olea (asterids) and Cymbidium (monocots) showed in each case a different ranking of genomic regions in terms of variability and potentially informative characters. Only two intergenic spacers (ndhF–rpl32 and trnK–rps16) were consistently found among the 30 top-ranked regions. We have mapped the occurrence of substitutions and microstructural mutations in the four genome pairs. High AT content in specific sequence elements seems to foster frequent mutations. We conclude that the variability among the fastest evolving plastid genomic regions is lineage-specific and thus cannot be precisely predicted across angiosperms. The often lineage-specific occurrence of stem-loop elements in the sequences of introns and spacers also governs lineage-specific mutations. Sequencing whole plastid genomes to find markers for evolutionary analyses is therefore particularly useful when overall genetic distances are low.

Chloroplast DNA sequence utility for the lowest phylogenetic and phylogeographic inferences in angiosperms: the tortoise and the hare IV

PREMISE OF THE STUDY

Noncoding chloroplast DNA (NC-cpDNA) sequences are the staple data source of low-level phylogeographic and phylogenetic studies of angiosperms. We followed up on previous papers (tortoise and hare II and III) that sought to identify the most consistently variable regions of NC-cpDNA. We used an exhaustive literature review and newly available whole plastome data to assess applicability of previous conclusions at low taxonomic levels.

METHODS

We aligned complete plastomes of 25 species pairs from across angiosperms, comparing the number of genetic differences found in 107 NC-cpDNA regions and matK. We surveyed Web of Science for the plant phylogeographic literature between 2007 and 2013 to assess how NC-cpDNA has been used at the intraspecific level.

KEY RESULTS

Several regions are consistently the most variable across angiosperm lineages: ndhF-rpl32, rpl32-trnL((UAG)), ndhC-trnV((UAC)), 5'rps16-trnQ((UUG)), psbE-petL, trnT((GGU))-psbD, petA-psbJ, and rpl16 intron. However, there is no universally best region. The average number of regions applied to low-level studies is ∼2.5, which may be too little to access the full discriminating power of this genome.

CONCLUSIONS

Plastome sequences have been used successfully at lower and lower taxonomic levels. Our findings corroborate earlier works, suggesting that there are regions that are most likely to be the most variable. However, while NC-cpDNA sequences are commonly used in plant phylogeographic studies, few of the most variable regions are applied in that context. Furthermore, it appears that in most studies too few NC-cpDNAs are used to access the discriminating power of the cpDNA genome.

Phylogenetic signal detection from an ancient rapid radiation: Effects of noise reduction, long-branch attraction, and model selection in crown clade Apocynaceae

Crown clade Apocynaceae comprise seven primary lineages of lianas, shrubs, and herbs with a diversity of pollen aggregation morphologies including monads, tetrads, and pollinia, making them an ideal group for investigating the evolution and function of pollen packaging. Traditional molecular systematic approaches utilizing small amounts of sequence data have failed to resolve relationships along the spine of the crown clade, a likely ancient rapid radiation. The previous best estimate of the phylogeny was a five-way polytomy, leaving ambiguous the homology of aggregated pollen in two major lineages, the Periplocoideae, which possess pollen tetrads, and the milkweeds (Secamonoideae plus Asclepiadoideae), which possess pollinia. To assess whether greatly increased character sampling would resolve these relationships, a plastome sequence data matrix was assembled for 13 taxa of Apocynaceae, including nine newly generated complete plastomes, one partial new plastome, and three previously reported plastomes, collectively representing all primary crown clade lineages and outgroups. The effects of phylogenetic noise, long-branch attraction, and model selection (linked versus unlinked branch lengths among data partitions) were evaluated in a hypothesis-testing framework based on Shimodaira-Hasegawa tests. Discrimination among alternative crown clade resolutions was affected by all three factors. Exclusion of the noisiest alignment positions and topologies influenced by long-branch attraction resulted in a trichotomy along the spine of the crown clade consisting of Rhabdadenia+the Asian clade, Baisseeae+milkweeds, and Periplocoideae+the New World clade. Parsimony reconstruction on all optimal topologies after noise exclusion unambiguously supports parallel evolution of aggregated pollen in Periplocoideae (tetrads) and milkweeds (pollinia). Our phylogenomic approach has greatly advanced the resolution of one of the most perplexing radiations in Apocynaceae, providing the basis for study of convergent floral morphologies and their adaptive value.