The flowering plant genus Aquilegia (columbine) is an important contributor to biodiversity and an example of both biotic and abiotic niche adaptation across much of the Northern Hemisphere, especially in California. Here we report a near-chromosome level draft genome assembly for A. eximia, a California endemic species. A. eximia is a serpentine-soil specialist and is very closely related to 2 columbine species also being studied for the California Conservation Genomics Project (CCGP), A. formosa (widespread) and A. pubescens (high alpine). Utilizing high throughput, long reads (PacBio) and chromatin capture (Omni-C), the A. eximia genome makes marked contiguity improvements compared to the existing reference genome for another North American columbine, A. coerulea "Goldsmith." The A. eximia genome will also be more useful for aligning whole genome resequencing data from California columbines than the genomes for more distantly related columbine species, the Asian A. oxysepala var. kansuensis and the European A. vulgaris. Notably, we found evidence that A. eximia, A. coerulea "Goldsmith," and A. vulgaris all share the same overall genome structure and differ from A. oxysepala var. kansuensis by the same reciprocal translocation. The A. eximia reference genome will be a valuable tool for identifying patterns of plant biodiversity across California for the CCGP, as well as for future population genomic and trait mapping studies.

Phyllospadix spp. (surfgrass) are flowering plants and keystone species in the rocky intertidal and subtidal environments of the North Pacific Ocean. Here we report a chromosome level assembly for P. torreyi, which occurs along the coast of California, sometimes in sympatry with P. scouleri. Both of these species and their putative hybrids are being studied as part of the California Conservation Genomics Project. Phyllospadix are dioecious, and males are exceptionally rare compared to females. Using high throughput, long reads (PacBio) and chromatin capture (Omni-C), we assembled a chromosome level genome for a male individual and a contig level assembly for a female individual. Comparison between the male and female assembly confirmed that the male is the heterogametic sex and has a massive Y chromosome at 124.8 megabases, which encompasses over 27% of the male genome. We also compared the male P. torreyi assembly to a genome from its sister genus, the monoecious Zostera marina, and found relatively high levels of synteny, that syntenic gene blocks on the P. torreyi sex chromosomes align to a single chromosome of Z. marina, and an estimated divergence time of ca. 25 million years ago. The Phyllospadix genome will be a powerful tool for studying marine dispersal, sex ratios, genetic diversity, sex chromosome evolution, and other dynamics in a keystone marine species.

The modification of fertile stamens into sterile staminodes has occurred independently many times in the flowering plant lineage. In the genus Aquilegia (columbine) and its closest relatives, the two stamen whorls closest to the carpels have been converted to staminodes. In Aquilegia, the only genetic analyses of staminode development have been reverse genetic approaches revealing that B-class floral identity genes are involved. A. jonesii, the only species of columbine where staminodes have reverted to fertile stamens, allows us to explore the genetic architecture of staminode development using a forward genetic approach. We performed QTL analysis using an outcrossed F2 population between A. jonesii and a horticultural variety that makes fully developed staminodes, A. coerulea 'Origami'. Our results reveal a polygenic basis for staminode loss where the two staminode whorls are under some level of independent control. We also discovered that staminode loss in A. jonesii is not complete, in which staminode-like traits sometimes occur in the inner fertile stamens, potentially representing a fading boundary of gene expression. The QTLs identified in this study provide a map to guide future reverse genetic and functional studies examining the genetic basis and evolutionary significance of this trait.

The California Conservation Genomics Project (CCGP) is a unique, critically important step forward in the use of comprehensive landscape genetic data to modernize natural resource management at a regional scale. We describe the CCGP, including all aspects of project administration, data collection, current progress, and future challenges. The CCGP will generate, analyze, and curate a single high-quality reference genome and 100-150 resequenced genomes for each of 153 species projects (representing 235 individual species) that span the ecological and phylogenetic breadth of California's marine, freshwater, and terrestrial ecosystems. The resulting portfolio of roughly 20 000 resequenced genomes will be analyzed with identical informatic and landscape genomic pipelines, providing a comprehensive overview of hotspots of within-species genomic diversity, potential and realized corridors connecting these hotspots, regions of reduced diversity requiring genetic rescue, and the distribution of variation critical for rapid climate adaptation. After 2 years of concerted effort, full funding ($12M USD) has been secured, species identified, and funds distributed to 68 laboratories and 114 investigators drawn from all 10 University of California campuses. The remaining phases of the CCGP include completion of data collection and analyses, and delivery of the resulting genomic data and inferences to state and federal regulatory agencies to help stabilize species declines. The aspirational goals of the CCGP are to identify geographic regions that are critical to long-term preservation of California biodiversity, prioritize those regions based on defensible genomic criteria, and provide foundational knowledge that informs management strategies at both the individual species and ecosystem levels.

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.

Here, we describe a polymorphic population of Aquilegia coerulea with a naturally occurring floral homeotic mutant, A. coerulea var. daileyae, where the characteristic petals with nectar spurs are replaced with a second set of sepals. Although it would be expected that this loss of pollinator reward would be disadvantageous to the mutant, we find that it has reached relatively high frequency (-25%) and is under strong, positive selection across multiple seasons (s = 0.17-0.3) primarily due to reduced floral herbivory. We identify the underlying locus (APETALA3-3) and multiple causal loss-of-function mutations indicating an ongoing soft sweep. Elevated linkage disequilibrium around the two most common causal alleles indicates that positive selection has been occurring for many generations. Lastly, genotypic frequencies at AqAP3-3 indicate a degree of positive assortative mating by morphology. Together, these data provide both a compelling example that large-scale discontinuous morphological changes differentiating taxa can occur due to single mutations and a particularly clear example of linking genotype, phenotype, and fitness.

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.

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.

Based on morphological and molecular phylogenetic studies, Semiaquilegia danxiashanensis, a new species from Danxia Shan in northern Guangdong, southern China, is described and illustrated. This species is easily distinguishable from each of other three known species in the genus by characters of the flowers and fruits. In addition, molecular phylogenetic analyses of both the nuclear ITS and the plastid trnL-F region strongly supported S. danxiashanensis as a separate species from other species of Semiaquilegia. We provide a detailed morphological and habitat description, distribution, as well as colour photographs and illustrations of the new species.

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.

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.

*For correspondence: scott.hodges@lifesci.ucsb.edu (SH); magnus.nordborg@gmi.oeaw.ac.at (MN) 1Gregor Mendel Institute, Austrian Academy of Sciences, Vienna BioCenter (VBC), Vienna, 11 Austria; 2Department of Ecology, Evolution, and Marine Biology, University of California 12 Santa Barbara, Santa Barbara, California, USA; 3Central-European Institute of Technology 13 (CEITEC), Masaryk University, Brno, Czech Republic; 4Vienna Graduate School of 14 Population Genetics, Vienna, Austria; 5Department of Energy Joint Genome Institute, 15 Walnut Creek, California, USA; 6HudsonAlpha Institute of Biotechnology, Huntsville, 16 Alabama, USA; 7Institute of Experimental Botany, Centre of the Region Haná for 17 Biotechnological and Agricultural Research, Olomouc, Czech Republic; 8Department of 18 Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts, USA 19

The columbine genus Aquilegia has been grown as a garden ornamental for centuries. Because of its unusual floral shape with the sepals being petaloid and colorful and the petals forming slender nectar spurs, many collections were made and species have been available to scientists as well as horticulturalists. Given its widespread natural distribution throughout the northern hemisphere, this availability was essential for early studies establishing that species in the genus were often highly intercompatible and their hybrids could be established. This led to a number of studies showing the relatively simple genetic inheritance of many of the dramatic differences among flowers of various species. More recent work has established the phylogenetic relationships among most of the species and showing that the genus had a burst of diversification after evolving the distinctive nectar spurs. Other studies have shown that this burst of diversification was likely due to adaptation to different pollinators, especially in North America, with changes in flower color, orientation and spur length to match the preferences and tongue lengths of new pollinators. Further studies have dissected the genetic basis of some of these traits. Aquilegia has also been noteworthy for its contributions to our understanding of the evolution of floral organ identity and horticultural varieties have, and will continue, to aid in these efforts. Aquilegia has now been the subject of substantial effort to develop it as a new model genomic system with the establishment of a high-quality reference sequence, derived from a highly inbred horticultural line as well as other resources such as a functional assay to assess specific gene functions. These resources poise Aquilegia to continue to be at the forefront of plant evolutionary and ecological research.

Speciation can be described as a reduction, and the eventual cessation, in the ability to interbreed. Thus, determining how gene flow differs within and between nascent species can illuminate the relative stage the taxa have attained in the speciation process. Aquilegia formosa and A. pubescens are fully intercompatible, yet occur in different habitats and have flowers specialized for pollination by hummingbirds and hawkmoths, respectively. Using 79 SNP loci, we genotyped nearly 1000 individuals from populations of both species in close proximity to each other and from putative hybrid zones. The species shared all but one SNP polymorphism, and on average, allele frequencies differed by only 0.14. However, the species were clearly differentiated using Structure, and admixed individuals were primarily identified at putative hybrid zones. PopGraph identified a highly integrated network among all populations, but populations of each species and hybrid zones occupied distinct regions in the network. Using either conditional graph distance (cGD) or Fst/(1‐Fst), we found significant isolation by distance (IBD) among populations. Within species, IBD was strong, indicating high historic gene flow. IBD extended approximately 100 km in A. pubescens and 30 km in A. formosa. However, IBD between the species was very weak and extended only a few km beyond hybrid zones, suggesting little recent gene flow. The extensive sharing of SNP polymorphisms between these species suggests that they are very early in the speciation process while the low signal of IBD suggests that they have largely ceased gene exchange.

Flowers of the lower eudicot Aquilegia (columbine) possess morphological innovations, namely elaborate petal spurs and a fifth distinct organ identity, the staminodium, that are well suited to the investigation of key questions in developmental,evolution. The recent evolution of these characteristics combined with a growing set of genetic and genomic resources has provided insight into how the traits arose and diversified. The petal spur appears to represent a key innovation that diversified largely via modification of specific aspects of cell expansion. In the case of the staminodium, gene duplication has played a role in allowing a novel organ identity to be carved out of the traditional ABC program.

Genetic marker studies can assist restoration practice through selection of seed sources that conserve historical levels of gene diversity and population genetic differentiation. We examined genetic variation and structure within and among mainland and island populations of Elymus glaucus, a perennial bunchgrass species native to western North American grasslands that is targeted for grassland restoration. Island populations of E. glaucus represent sensitive sites and potentially distinctive seed sources for reintroduction, and little is known of their genetic composition. Genetic diversity and structure were estimated using amplified fragment length polymorphism markers for 21 populations and 416 individuals distributed across two coastal California mainland locations and three California Channel Islands. Eight primer combinations resulted in 166 markers, of which 165 (99.4 %) were polymorphic. The number of polymorphic bands was significantly greater among mainland populations relative to island sites, and locally common alleles were present for each sampled island and mainland location. Population structure was high (62.9 %), with most variation (55.8 %) distributed among populations, 7.1 % between mainland and island locations, and the remainder (37.1 %) within populations. Isolation by distance was only apparent among islands. Using marker data to recommend appropriate seed sources for restoration, E. glaucus seeds are best derived within islands with collections representing a large number of individuals from matching environments. Given the limited gene flow and prior evidence of adaptive divergence among populations of this species, regional collections are recommended in all cases to maintain diversity and to avoid long-distance introductions of highly differentiated plant material.

Aquilegia is a well-known model system in the field of evolutionary biology, but obtaining a resolved and well-supported phylogenetic reconstruction for the genus has been hindered by its recent and rapid diversification. Here, we applied 454 next-generation sequencing to PCR amplicons of 21 of the most rapidly evolving regions of the plastome to generate c. 24kb of sequences from each of 84 individuals from throughout the genus. The resulting phylogeny has well-supported resolution of the main lineages of the genus, although recent diversification such as in the European taxa remains unresolved. By producing a chronogram of the whole Ranunculaceae family based on published data, we inferred calibration points for dating the Aquilegia radiation. The genus originated in the upper Miocene c.6.9 million yr ago (Ma) in Eastern Asia, and diversification occurred c. 4.8Ma with the split of two main clades, one colonizing North America, and the other Western Eurasia through the mountains of Central Asia. This was followed by a back-to-Asia migration, originating from the European stock using a North Asian route. These results provide the first backbone phylogeny and spatiotemporal reconstruction of the Aquilegia radiation, and constitute a robust framework to address the adaptative nature of speciation within the group.

Petals show enormous diversity in appearance and function, and play key roles in angiosperm evolution.Petals have been lost multiple times, giving rise to the occurrences of apetalous flowers in different lineages.

The role of petal spurs and specialized pollinator interactions has been studied since Darwin. Aquilegia petal spurs exhibit striking size and shape diversity, correlated with specialized pollinators ranging from bees to hawkmoths in a textbook example of adaptive radiation. Despite the evolutionary significance of spur length, remarkably little is known about Aquilegia spur morphogenesis and its evolution. Using experimental measurements, both at tissue and cellular levels, combined with numerical modelling, we have investigated the relative roles of cell divisions and cell shape in determining the morphology of the Aquilegia petal spur. Contrary to decades-old hypotheses implicating a discrete meristematic zone as the driver of spur growth, we find that Aquilegia petal spurs develop via anisotropic cell expansion. Furthermore, changes in cell anisotropy account for 99 per cent of the spur-length variation in the genus, suggesting that the true evolutionary innovation underlying the rapid radiation of Aquilegia was the mechanism of tuning cell shape.