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

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.

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.

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.

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.

Genetic architecture of floral traits in Iris hexagona and Iris fulva.

The formation of hybrids among closely related species has been observed in numerous plant taxa. Selection by pollinators on floral traits can act as an early reproductive isolating barrier and may be especially important when there is overlap in distribution and flowering time. In this study, we use Quantitative Trait Locus (QTL) mapping based on 293 codominant SNP markers in an F2 population (n = 328) to assess the size, magnitude, and location of the genetic regions controlling floral traits known to be important for pollinator attraction in 2 species of Lousiana Irises, Iris fulva and Iris hexagona. We also evaluate correlations among F2 traits and identify transgression in the hybrid population. Overall, we observe that differences in most floral traits between I. fulva and I. hexagona are controlled by multiple QTLs and are distributed across several linkage groups. We also find evidence of transgression at several QTL, suggesting that hybridization can contribute to generating phenotypic variation, which may be adaptive in rapidly changing environments.

Transcriptome sequencing and phylogenetic analysis of floral and leaf MIKC(C) MADS-box and R2R3 MYB transcription factors from the monocot Iris fulva.

The Louisiana Irises serve as an important system for the study of the evolutionary processes of speciation, including reproductive isolation, hybridization, and adaptation. Sequencing methods today allow for the generation of resources key to elucidating the genetic basis of these phenomena. Here we describe the transcriptomes of floral and young leaf tissue from Iris fulva generated by massively parallel pyrosequencing. In order to identify potential candidates for the study of reproductive isolation and adaptation in the Louisiana Irises we phylogenetically analyzed the type II MIKC(C) MADS-box and R2R3 MYB transcription factors expressed in these tissues. A total of 25 Iris MIKC(C) MADS-box genes in 9 clades and 42 Iris R2R3 MYB genes in 19 clades were identified. Through the identification of eudicot and monocot specific clades, these analyses contribute to our understanding of the evolution of these prominent transcription factor families in the angiosperms.

QTL mapping reveals the genetic architecture of loci affecting pre- and post-zygotic isolating barriers in Louisiana Iris.

BACKGROUND: Hybridization among Louisiana Irises has been well established and the genetic architecture of reproductive isolation is known to affect the potential for and the directionality of introgression between taxa. Here we use co-dominant markers to identify regions where QTL are located both within and between backcross maps to compare the genetic architecture of reproductive isolation and fitness traits across treatments and years. RESULTS: QTL mapping was used to elucidate the genetic architecture of reproductive isolation between Iris fulva and Iris brevicaulis. Homologous co-dominant EST-SSR markers scored in two backcross populations between I. fulva and I. brevicaulis were used to generate genetic linkage maps. These were used as the framework for mapping QTL associated with variation in 11 phenotypic traits likely responsible for reproductive isolation and fitness. QTL were dispersed throughout the genome, with the exception of one region of a single linkage group (LG) where QTL for flowering time, sterility, and fruit production clustered. In most cases, homologous QTL were not identified in both backcross populations, however, homologous QTL for flowering time, number of growth points per rhizome, number of nodes per inflorescence, and number of flowers per node were identified on several linkage groups. CONCLUSIONS: Two different traits affecting reproductive isolation, flowering time and sterility, exhibit different genetic architectures, with numerous QTL across the Iris genome controlling flowering time and fewer, less distributed QTL affecting sterility. QTL for traits affecting fitness are largely distributed across the genome with occasional overlap, especially on LG 4, where several QTL increasing fitness and decreasing sterility cluster. Given the distribution and effect direction of QTL affecting reproductive isolation and fitness, we have predicted genomic regions where introgression may be more likely to occur (those regions associated with an increase in fitness and unlinked to loci controlling reproductive isolation) and those that are less likely to exhibit introgression (those regions linked to traits decreasing fitness and reproductive isolation).

Environmental and molecular analysis of the floral transition in the lower eudicot Aquilegia formosa.

BACKGROUND: Flowering is a critical transition in plant development, the timing of which can have considerable fitness consequences. Until recently, research into the genetic control of flowering time and its associated developmental changes was focused on core eudicots (for example, Arabidopsis) or monocots (for example, Oryza). Here we examine the flowering response of Aquilegia formosa, a member of the eudicot order Ranunculales that is emerging as an important model for the investigation of plant ecology and evolution. RESULTS: We have determined that A. formosa has a strong vernalization requirement but little or no photoperiod response, making it a day neutral (DN) plant. Consistent with this, the Aquilegia homolog of FLOWERING LOCUS T (AqFT) is expressed in both long and short days but surprisingly, the locus is expressed before the transition to flowering. In situ hybridizations with homologs of several Arabidopsis Floral Pathway Integrators (FPIs) do not suggest conserved functions relative to Arabidopsis, the potential exceptions being AqLFY and AqAGL24.2. CONCLUSIONS: In Aquilegia, vernalization is critical to flowering but this signal is not strictly required for the transcriptional activation of AqFT. The expression patterns of AqLFY and AqAGL24.2 suggest a hypothesis for the development of Aquilegia's determinate inflorescence whereby their differential expression controls the progression of each meristem from inflorescence to floral identity. Interestingly, none of the Aquilegia expression patterns are consistent with a function in floral repression which, combined with the lack of a FLC homolog, means that new candidate genes must be identified for the control of vernalization response in Aquilegia.

In the Light of Evolution: A Reevaluation of Conservation in the CO-FT Regulon and Its Role in Photoperiodic Regulation of Flowering Time.

In order to maximize reproductive success, plants have evolved different strategies to control the critical developmental shift marked by the transition to flowering. As plants have adapted to diverse environments across the globe, these strategies have evolved to recognize and respond to local seasonal cues through the induction of specific downstream genetic pathways, thereby ensuring that the floral transition occurs in favorable conditions. Determining the genetic factors involved in controlling the floral transition in many species is key to understanding how this trait has evolved. Striking genetic discoveries in Arabidopsis thaliana (Arabidopsis) and Oryza sativa (rice) revealed that similar genes in both species control flowering in response to photoperiod, suggesting that this genetic module could be conserved between distantly related angiosperms. However, as we have gained a better understanding of the complex evolution of these genes and their functions in other species, another possibility must be considered: that the genetic module controlling flowering in response to photoperiod is the result of convergence rather than conservation. In this review, we show that while data clearly support a central role of FLOWERING LOCUS T (FT) homologs in floral promotion across a diverse group of angiosperms, there is little evidence for a conserved role of CONSTANS (CO) homologs in the regulation of these loci. In addition, although there is an element of conserved function for FT homologs, even this component has surprising complexity in its regulation and evolution.