A dominant-negative actin mutation alters corolla tube width and pollinator visitation in Mimulus lewisii

Evolutionary plant-pollinator responses to anthropogenic land-use change: impacts on ecosystem services

Agricultural intensification at field and landscape scales, including increased use of agrochemicals and loss of semi-natural habitats, is a major driver of insect declines and other community changes. Efforts to understand and mitigate these effects have traditionally focused on ecological responses. At the same time, adaptations to pesticide use and habitat fragmentation in both insects and flowering plants show the potential for rapid evolution. Yet we lack an understanding of how such evolutionary responses may propagate within and between trophic levels with ensuing consequences for conservation of species and ecological functions in agroecosystems. Here, we review the literature on the consequences of agricultural intensification on plant and animal evolutionary responses and interactions. We present a novel conceptualization of evolutionary change induced by agricultural intensification at field and landscape scales and emphasize direct and indirect effects of rapid evolution on ecosystem services. We exemplify by focusing on economically and ecologically important interactions between plants and pollinators. We showcase available eco-evolutionary theory and plant-pollinator modelling that can improve predictions of how agricultural intensification affects interaction networks, and highlight available genetic and trait-focused methodological approaches. Specifically, we focus on how spatial genetic structure affects the probability of propagated responses, and how the structure of interaction networks modulates effects of evolutionary change in individual species. Thereby, we highlight how combined trait-based eco-evolutionary modelling, functionally explicit quantitative genetics, and genomic analyses may shed light on conditions where evolutionary responses impact important ecosystem services.

A single amino acid substitution in MdLAZY1A dominantly impairs shoot gravitropism in Malus

Plant architecture is 1 of the most important factors that determines crop yield potential and productivity. In apple (Malus domestica), genetic improvement of tree architecture has been challenging due to a long juvenile phase and growth as complex trees composed of a distinct scion and a rootstock. To better understand the genetic control of apple tree architecture, the dominant weeping growth phenotype was investigated. We report the identification of MdLAZY1A (MD13G1122400) as the genetic determinant underpinning the Weeping (W) locus that largely controls weeping growth in Malus. MdLAZY1A is 1 of the 4 paralogs in apple that are most closely related to AtLAZY1 involved in gravitropism in Arabidopsis (Arabidopsis thaliana). The weeping allele (MdLAZY1A-W) contains a single nucleotide mutation c.584T>C that leads to a leucine to proline (L195P) substitution within a predicted transmembrane domain that colocalizes with Region III, 1 of the 5 conserved regions in LAZY1-like proteins. Subcellular localization revealed that MdLAZY1A localizes to the plasma membrane and nucleus in plant cells. Overexpressing the weeping allele in apple cultivar Royal Gala (RG) with standard growth habit impaired its gravitropic response and altered the growth to weeping-like. Suppressing the standard allele (MdLAZY1A-S) by RNA interference (RNAi) in RG similarly changed the branch growth direction to downward. Overall, the L195P mutation in MdLAZY1A is genetically causal for weeping growth, underscoring not only the crucial roles of residue L195 and Region III in MdLAZY1A-mediated gravitropic response but also a potential DNA base editing target for tree architecture improvement in Malus and other crops.

Integrating BSA-Seq with RNA-Seq Reveals a Novel Fasciated Ear5 Mutant in Maize

Increasing grain yield is required to meet the rapidly expanding demands for food, feed, and fuel. Inflorescence meristems are central to plant growth and development. However, the question concerning whether inflorescence development can be regulated to improve grain yield remains unclear. Here, we describe a naturally occurring single recessive mutation called fea5 that can increase grain yield in maize. Using bulk segregant analysis sequencing (BSA-seq), the candidate region was initially mapped to a large region on chromosome 4 (4.68 Mb-11.26 Mb). Transcriptome sequencing (RNA-seq) revealed a total of 1246 differentially expressed genes (DEGs), of which 835 were up-regulated and 411 were down-regulated. Further analysis revealed the enrichment of DEGs in phytohormone signal transduction. Consistently, phytohormone profiling indicated that auxin (IAA), jasmonic acid (JA), ethylene (ETH), and cytokinin (CK) levels increased significantly, whereas the gibberellin (GA) level decreased significantly in fea5. By integrating BSA-seq with RNA-seq, we identified Zm00001d048841 as the most likely candidate gene. Our results provide valuable insight into this new germplasm resource and the molecular mechanism underlying fasciated ears that produce a higher kernel row number in maize.

Application of Optical Fluorescence Spectroscopy for Studying Bee Abundance in Tropaeolum majus L. (Tropaeolaceae)

Tropaeolum majus L. species produce flowers with all sorts of colors, from yellow to red. This work aimed to apply optical fluorescence spectroscopy to study bee abundance in T. majus, answering the following questions: (1) do corolla temperature and weather conditions affect the abundance of visiting bee species? (2) do flower color and corolla fluorescence affect the abundance of visiting bee species? (3) do red flowers attract more visiting bees? (4) is there a relationship between bee visits and flower compounds? The bee species Apis mellifera, Paratrigona lineata, and Trigona spinipes were the most observed in T. majus flowers. The latter was more active in the morning and preferred orange and yellow flowers. These colors also had higher temperatures and fluorescence emissions than did the red ones and those with yellow-red and orange-red nectar guides. Orange flowers emitted a broadband UV spectrum (between 475 and 800 nm). This range might be due to compounds such as hydroxycinnamic acid, flavonols, isoflavonoids, flavones, phenolic acid, and chlorophyll. Extracts from different T. majus corolla colors showed that flowers emit specific fluorescent signals, mainly related to bee color vision and learning, thus acting as a means of communication between bees and flowers. In this way, this information evidences the interaction between bees and T. majus flowers, allowing conservation actions for pollinators.

The leaf polarity factors SGS3 and YABBYs regulate style elongation through auxin signaling in Mimulus lewisii

Style length is a major determinant of breeding strategies in flowering plants and can vary dramatically between and within species. However, little is known about the genetic and developmental control of style elongation. We characterized the role of two classes of leaf adaxial-abaxial polarity factors, SUPPRESSOR OF GENE SILENCING3 (SGS3) and the YABBY family transcription factors, in the regulation of style elongation in Mimulus lewisii. We also examined the spatiotemporal patterns of auxin response during style development. Loss of SGS3 function led to reduced style length via limiting cell division, and downregulation of YABBY genes by RNA interference resulted in shorter styles by decreasing both cell division and cell elongation. We discovered an auxin response minimum between the stigma and ovary during the early stages of pistil development that marks style differentiation. Subsequent redistribution of auxin response to this region was correlated with style elongation. Auxin response was substantially altered when both SGS3 and YABBY functions were disrupted. We suggest that auxin signaling plays a central role in style elongation and that the way in which auxin signaling controls the different cell division and elongation patterns underpinning natural style length variation is a major question for future research.

Mapping of Genetic Locus for Leaf Trichome Formation in Chinese Cabbage Based on Bulked Segregant Analysis

Chinese cabbage is a leafy vegetable, and its leaves are the main edible organs. The formation of trichomes on the leaves can significantly affect its taste, so studying this phenomenon is of great significance for improving the quality of Chinese cabbage. In this study, two varieties of Chinese cabbage, W30 with trichome leaves and 082 with glabrous leaves, were crossed to generate F₁ and F₁ plants, which were self-fertilized to develop segregating populations with trichome or glabrous morphotypes. The two bulks of the different segregating populations were used to conduct bulked segregant analysis (BSA). A total of 293.4 M clean reads were generated from the samples, and plants from the trichome leaves (AL) bulk and glabrous leaves (GL) bulk were identified. Between the two DNA pools generated from the trichome and glabrous plants, 55,048 SNPs and 272 indels were generated. In this study, three regions (on chromosomes 6, 10 and scaffold000100) were identified, and the annotation revealed three candidate genes that may participate in the formation of leaf trichomes. These findings suggest that the three genes-Bra025087 encoding a cyclin family protein, Bra035000 encoding an ATP-binding protein/kinase/protein kinase/protein serine/threonine kinase and Bra033370 encoding a WD-40 repeat family protein-influence the formation of trichomes by participating in trichome morphogenesis (GO: 0010090). These results demonstrate that BSA can be used to map genes associated with traits and provide new insights into the molecular mechanism of leafy trichome formation in Chinese cabbage.

Comparative analysis of corolla tube development across three closely related Mimulus species with different pollination syndromes

Fusion of petals to form a corolla tube is considered a key innovation contributing to the diversification of many flowering plant lineages. Corolla tube length often varies dramatically among species and is a major determinant of pollinator preference. However, our understanding of the developmental dynamics underlying corolla tube length variation is very limited. Here we examined corolla tube growth in the Mimulus lewisii species complex, an emerging model system for studying the developmental genetics and evo-devo of pollinator-associated floral traits. We compared developmental and cellular processes associated with corolla tube length variation among the bee-pollinated M. lewisii, the hummingbird-pollinated Mimulus verbenaceus, and the self-pollinated Mimulus parishii. We found that in all three species, cell size is non-uniformly distributed along the mature tube, with the longest cells just distal to the stamen insertion site. Differences in corolla tube length among the three species are not associated with processes of organogenesis or early development but are associated with variation in multiple processes occurring later in development, including the location and duration of cell division and cell elongation. The tube growth curves of the small-flowered M. parishii and large-flowered M. lewisii are essentially indistinguishable, except that M. parishii tubes stop growing earlier at a smaller size, suggesting a critical role of heterochrony in the shift from outcrossing to selfing. These results not only highlight the developmental process associated with corolla tube variation among species but also provide a baseline reference for future developmental genetic analyses of mutants or transgenic plants with altered corolla tube morphology in this emerging model system.

Developmental Genetics of Corolla Tube Formation: Role of the tasiRNA-ARF Pathway and a Conceptual Model

Over 80,000 angiosperm species produce flowers with petals fused into a corolla tube. The corolla tube contributes to the tremendous diversity of flower morphology and plays a critical role in plant reproduction, yet it remains one of the least understood plant structures from a developmental genetics perspective. Through mutant analyses and transgenic experiments, we show that the tasiRNA-ARF pathway is required for corolla tube formation in the monkeyflower species Mimulus lewisii Loss-of-function mutations in the M. lewisii orthologs of ARGONAUTE7 and SUPPRESSOR OF GENE SILENCING3 cause a dramatic decrease in abundance of TAS3-derived small RNAs and a moderate upregulation of AUXIN RESPONSE FACTOR3 (ARF3) and ARF4, which lead to inhibition of lateral expansion of the bases of petal primordia and complete arrest of the upward growth of the interprimordial regions, resulting in unfused corollas. Using the DR5 auxin-responsive promoter, we discovered that auxin signaling is continuous along the petal primordium base and the interprimordial region during the critical stage of corolla tube formation in the wild type, similar to the spatial pattern of MlARF4 expression. Auxin response is much weaker and more restricted in the mutant. Furthermore, exogenous application of a polar auxin transport inhibitor to wild-type floral apices disrupted petal fusion. Together, these results suggest a new conceptual model highlighting the central role of auxin-directed synchronized growth of the petal primordium base and the interprimordial region in corolla tube formation.

Characterization of four polymorphic genes controlling red leaf colour in lettuce that have undergone disruptive selection since domestication

Anthocyanins protect plants from biotic and abiotic stressors and provide great health benefits to consumers. In this study, we cloned four genes (Red Lettuce Leaves 1 to 4: RLL1 to RLL4) that contribute to colour variations in lettuce. The RLL1 gene encodes a bHLH transcription factor, and a 5-bp deletion in some cultivars abolishes its function to activate the anthocyanin biosynthesis pathway. The RLL2 gene encodes an R2R3-MYB transcription factor, which was derived from a duplication followed by mutations in its promoter region. The RLL3 gene encodes an R2-MYB transcription factor, which down-regulates anthocyanin biosynthesis through competing with RLL2 for interaction with RLL1; a mis-sense mutation compromises the capacity of RLL3 to bind RLL1. The RLL4 gene encodes a WD-40 transcription factor, homologous to the RUP genes suppressing the UV-B signal transduction pathway in Arabidopsis; a mis-sense mutation in rll4 attenuates its suppressing function, leading to a high concentration of anthocyanins. Sequence analysis of the RLL1-RLL4 genes from wild and cultivated lettuce showed that their function-changing mutations occurred after domestication. The mutations in rll1 disrupt anthocyanin biosynthesis, while the mutations in RLL2, rll3 and rll4 activate anthocyanin biosynthesis, showing disruptive selection for leaf colour during domestication of lettuce. The characterization of multiple polymorphic genes in this study provides the necessary molecular resources for the rational breeding of lettuce cultivars with distinct levels of red pigments and green cultivars with high levels of health-promoting flavonoids.

Nectary size is a pollination syndrome trait in Penstemon

Evolution of complex phenotypes depends on the adaptive importance of individual traits, and the developmental changes required to modify traits. Floral syndromes are complex adaptations to pollinators that include color, nectar, and shape variation. Hummingbird-adapted flowers have evolved a remarkable number of times from bee-adapted ancestors in Penstemon, and previous work demonstrates that color over shape better distinguishes bee from hummingbird syndromes. Here, we examined the relative importance of nectar volume and nectary development in defining Penstemon pollination syndromes. We tested the evolutionary association of nectar volume and nectary area with pollination syndrome across 19 Penstemon species. In selected species, we assessed cellular-level processes shaping nectary size. Within a segregating population from an intersyndrome cross, we assessed trait correlations between nectar volume, nectary area, and the size of stamens on which nectaries develop. Nectar volume and nectary area displayed an evolutionary association with pollination syndrome. These traits were correlated within a genetic cross, suggesting a mechanistic link. Nectary area evolution involves parallel processes of cell expansion and proliferation. Our results demonstrate that changes to nectary patterning are an important contributor to pollination syndrome diversity and provide further evidence that repeated origins of hummingbird adaptation involve parallel developmental processes in Penstemon.

Fine-mapping and transcriptome analysis of BoGL-3, a wax-less gene in cabbage (Brassica oleracea L. var. capitata)

The great majority of terrestrial plants produce epicuticular wax that is used to protect plants from a variety of biotic and abiotic stresses. Cabbage epicuticular wax is a white crystalline compound of various lipids. Wax-less cabbage has the characteristics of lustrous green leaves and beautiful exterior, which facilitates the brilliant green cabbage breeding. CGL-3 is a spontaneous wax-less mutant identified from cabbage. Genetic analysis indicated that the waxy deficiency of the mutant was controlled by a single dominant gene. To clarify the mechanism of the waxy deficiency, fine-mapping and transcriptome analysis of the wax-less gene, BoGL-3, were carried out in this study. The result of fine mapping showed that the wax-less gene, BoGL-3, was delimited in a 33.5-kb interval which is between the flanking marker C08-98 and the end of chromosome 8. Two cDNA libraries, constructed with wax-less cabbage CGL-3 and the wild-type cabbage WT, were sequenced for screening of the target gene BoGL-3. A total of 8340 genes were identified with significant differential expression between CGL-3 and WT. Among these genes, 3187 were up-regulated and 5153 were down-regulated in CGL-3. With homologous analysis, four differential expressed genes related to wax metabolism were obtained. Among these four genes, only Bol018504 is located within the region of fine-mapping. Bol08504 is homologous to CER1, which encodes fatty acid hydroxylase and plays an important role in wax synthesis in Arabidopsis. However, there was no difference of Bol08504 sequence between CGL-3 and WT. We suggested that Bol018504 was regulated by BoGL-3. The suppression of Bol018504 leads to the reduction of wax. These findings will be helpful to reveal the mechanism of the wax metabolism in cabbage and develop lustrous green cabbage germplasm material.

Monkeyflowers (Mimulus): new model for plant developmental genetics and evo-devo

Contents Summary 694 I. Introduction 694 II. The system 695 III. Regulation of carotenoid pigmentation 695 IV. Formation of periodic pigmentation patterns 696 V. Developmental genetics of corolla tube formation and elaboration 697 VI. Molecular basis of floral trait variation underlying pollinator shift 698 VII. Outlook 699 Acknowledgements 699 References 699 SUMMARY: Monkeyflowers (Mimulus) have long been recognized as a classic ecological and evolutionary model system. However, only recently has it been realized that this system also holds great promise for studying the developmental genetics and evo-devo of important plant traits that are not found in well-established model systems such as Arabidopsis. Here, I review recent progress in four different areas of plant research enabled by this new model, including transcriptional regulation of carotenoid biosynthesis, formation of periodic pigmentation patterns, developmental genetics of corolla tube formation and elaboration, and the molecular basis of floral trait divergence underlying pollinator shift. These examples suggest that Mimulus offers ample opportunities to make exciting discoveries in plant development and evolution.

Of puzzles and pavements: a quantitative exploration of leaf epidermal cell shape

Epidermal cells of leaves are diverse: tabular pavement cells, trichomes, and stomatal complexes. Pavement cells from the monocot Zea mays (maize) and the eudicot Arabidopsis thaliana (Arabidopsis) have highly undulate anticlinal walls. The molecular basis for generating these undulating margins has been extensively investigated in these species. This has led to two assumptions: first, that particular plant lineages are characterized by particular pavement cell shapes; and second, that undulatory cell shapes are common enough to be model shapes. To test these assumptions, we quantified pavement cell shape in epidermides from the leaves of 278 vascular plant taxa. We found that monocot pavement cells tended to have weakly undulating margins, fern cells had strongly undulating margins, and eudicot cells showed no particular undulation degree. Cells with highly undulating margins, like those of Arabidopsis and maize, were in the minority. We also found a trend towards more undulating cell margins on abaxial leaf surfaces; and that highly elongated leaves in ferns, monocots and gymnosperms tended to have highly elongated cells. Our results reveal the diversity of pavement cell shapes, and lays the quantitative groundwork for testing hypotheses about pavement cell form and function within a phylogenetic context.

Of puzzles and pavements: a quantitative exploration of leaf epidermal cell shape

Epidermal cells of leaves are diverse: tabular pavement cells, trichomes, and stomatal complexes. Pavement cells from the monocot Zea mays (maize) and the eudicot Arabidopsis thaliana (Arabidopsis) have highly undulate anticlinal walls. The molecular basis for generating these undulating margins has been extensively investigated in these species. This has led to two assumptions: first, that particular plant lineages are characterized by particular pavement cell shapes; and second, that undulatory cell shapes are common enough to be model shapes. To test these assumptions, we quantified pavement cell shape in epidermides from the leaves of 278 vascular plant taxa. We found that monocot pavement cells tended to have weakly undulating margins, fern cells had strongly undulating margins, and eudicot cells showed no particular undulation degree. Cells with highly undulating margins, like those of Arabidopsis and maize, were in the minority. We also found a trend towards more undulating cell margins on abaxial leaf surfaces; and that highly elongated leaves in ferns, monocots and gymnosperms tended to have highly elongated cells. Our results reveal the diversity of pavement cell shapes, and lays the quantitative groundwork for testing hypotheses about pavement cell form and function within a phylogenetic context.

Cloning and characterization of a Mimulus lewisii NPR1 gene involved in regulating plant resistance to Rhizoctonia solani

The monkey flower Mimulus lewisii is a new emerging model plant for the study in corolla tube formation, pigmentation patterns and pollinator selection, etc. However, the cultivation and management of this plant are difficult due to its susceptibility to a wide range of pathogens and the lack of rigid varieties with high levels of resistance to pathogens. In this regard, genetic engineering is a promising tool that may possibly allow us to enhance the M. lewisii disease resistance against pathogens. Here, we reported the isolation and characterization of non-expressor of pathogenesis related gene 1 (NPR1) gene from M. lewisii. The phylogenetic tree constructed based on the deduced sequence of MlNPR1 with homologs from other species revealed that MlNPR1 grouped together with other known NPR1 proteins of Scrophulariaceae family, and was nearest to Mimulus guttatus. Furthermore, expression analysis showed that MlNPR1 was upregulated after SA treatment and fungal infection. To understand the defensive role of this gene, we overexpressed MlNPR1 in M. lewisii. The transgenic lines showed slight phenotypic abnormalities, but constitutive expression of MlNPR1 activates defense signaling pathways by priming the expression of antifungal PR genes. Moreover, MlNPR1 transgenic lines showed enhanced resistance to Rhizoctonia solani there was delay in symptoms and reduced disease severity than non-transgenic plants. Altogether, the present study suggests that increasing the expression level of MlNPR1 may be a promising approach for development of monkey flower cultivars with enhanced resistance to diseases.

Bulked-Segregant Analysis Coupled to Whole Genome Sequencing (BSA-Seq) for Rapid Gene Cloning in Maize

Forward genetics remains a powerful method for revealing the genes underpinning organismal form and function, and for revealing how these genes are tied together in gene networks. In maize, forward genetics has been tremendously successful, but the size and complexity of the maize genome made identifying mutant genes an often arduous process with traditional methods. The next generation sequencing revolution has allowed for the gene cloning process to be significantly accelerated in many organisms, even when genomes are large and complex. Here, we describe a bulked-segregant analysis sequencing (BSA-Seq) protocol for cloning mutant genes in maize. Our simple strategy can be used to quickly identify a mapping interval and candidate single nucleotide polymorphisms (SNPs) from whole genome sequencing of pooled F2 individuals. We employed this strategy to identify narrow odd dwarf as an enhancer of teosinte branched1, and to identify a new allele of defective kernel1 Our method provides a quick, simple way to clone genes in maize.

Molecular Basis of Overdominance at a Flower Color Locus

Single-gene overdominance is one of the major mechanisms proposed to explain heterosis (i.e., hybrid vigor), the phenomenon that hybrid offspring between two inbred lines or varieties show superior phenotypes to both parents. Although sporadic examples of single-gene overdominance have been reported over the decades, the molecular nature of this phenomenon remains poorly understood and it is unclear whether any generalizable principle underlies the various cases. Through bulk segregant analysis, chemical profiling, and transgenic experiments, we show that loss-of-function alleles of the FLAVONE SYNTHASE (FNS) gene cause overdominance in anthocyanin-based flower color intensity in the monkeyflower species Mimulus lewisii FNS negatively affects flower color intensity by competing with the anthocyanin biosynthetic enzymes for the same substrates, yet positively affects flower color intensity by producing flavones, the colorless copigments required for anthocyanin stabilization, leading to enhanced pigmentation in the heterozyote (FNS/fns) relative to both homozygotes (FNS/FNS and fns/fns). We suggest that this type of antagonistic pleiotropy (i.e., alleles with opposing effects on different components of the phenotypic output) might be a general principle underlying single-gene overdominance.

Association between Petal Form Variation and CYC2-like Genotype in a Hybrid Line of Sinningia speciosa

This study used three-dimensional (3D) micro-computed tomography (μCT) imaging to examine petal form variation in a hybrid cross of Sinningia speciosa between a cultivar with actinomorphic flowers and a variety with zygomorphic flowers. The major objectives were to determine the genotype-phenotype associations between the petal form variation and CYCLOIDEA2-like alleles in S. speciosa (SsCYC) and to morphologically investigate the differences in petal types between actinomorphic and zygomorphic flowers. In this study, μCT was used to accurately acquire 3D floral images. Landmark-based geometric morphometrics (GM) was applied to evaluate the major form variations of the petals. Nine morphological traits of the petals were defined according to the form variations quantified through the GM analysis. The results indicated that the outward curvature of dorsal petals, the midrib asymmetry of lateral petals, and the dilation of ventral region of the tube were closely associated with the SsCYC genotype. Multiple analyses of form similarity between the petals suggested that the dorsal and ventral petals of actinomorphic plants resembled the ventral petals of zygomorphic plants. This observation indicated that the transition from zygomorphic to actinomorphic flowers in S. speciosa might be caused by the ventralization of the dorsal petals. We demonstrated that the 3D-GM approach can be used to determine genotype-phenotype associations and to provide morphological evidence for the transition of petal types between actinomorphic and zygomorphic flowers in S. speciosa.

Comparative transcriptome analyses of flower development in four species of Achimenes (Gesneriaceae)

BACKGROUND

Flowers have an amazingly diverse display of colors and shapes, and these characteristics often vary significantly among closely related species. The evolution of diverse floral form can be thought of as an adaptive response to pollination and reproduction, but it can also be seen through the lens of morphological and developmental constraints. To explore these interactions, we use RNA-seq across species and development to investigate gene expression and sequence evolution as they relate to the evolution of the diverse flowers in a group of Neotropical plants native to Mexico-magic flowers (Achimenes, Gesneriaceae).

RESULTS

The assembled transcriptomes contain between 29,000 and 42,000 genes expressed during development. We combine sequence orthology and coexpression clustering with analyses of protein evolution to identify candidate genes for roles in floral form evolution. Over 25% of transcripts captured were distinctive to Achimenes and overrepresented by genes involved in transcription factor activity. Using a model-based clustering approach we find dynamic, temporal patterns of gene expression among species. Selection tests provide evidence of positive selection in several genes with roles in pigment production, flowering time, and morphology. Combining these approaches to explore genes related to flower color and flower shape, we find distinct patterns that correspond to transitions of floral form among Achimenes species.

CONCLUSIONS

The floral transcriptomes developed from four species of Achimenes provide insight into the mechanisms involved in the evolution of diverse floral form among closely related species with different pollinators. We identified several candidate genes that will serve as an important and useful resource for future research. High conservation of sequence structure, patterns of gene coexpression, and detection of positive selection acting on few genes suggests that large phenotypic differences in floral form may be caused by genetic differences in a small set of genes. Our characterized floral transcriptomes provided here should facilitate further analyses into the genomics of flower development and the mechanisms underlying the evolution of diverse flowers in Achimenes and other Neotropical Gesneriaceae.