The role of phenotypic plasticity in shaping ecological networks.

Plasticity-mediated changes in interaction dynamics and structure may scale up and affect the ecological network in which the plastic species are embedded. Despite their potential relevance for understanding the effects of plasticity on ecological communities, these effects have seldom been analysed. We argue here that, by boosting the magnitude of intra-individual phenotypic variation, plasticity may have three possible direct effects on the interactions that the plastic species maintains with other species in the community: may expand the interaction niche, may cause a shift from one interaction niche to another or may even cause the colonization of a new niche. The combined action of these three factors can scale to the community level and eventually expresses itself as a modification in the topology and functionality of the entire ecological network. We propose that this causal pathway can be more widespread than previously thought and may explain how interaction niches evolve quickly in response to rapid changes in environmental conditions. The implication of this idea is not solely eco-evolutionary but may also help to understand how ecological interactions rewire and evolve in response to global change.

Hybridization and introgression are prevalent in Southern European Erysimum (Brassicaceae) species.

BACKGROUND AND AIMS: Hybridization is a common and important force in plant evolution. One of its outcomes is introgression - the transfer of small genomic regions from one taxon to another by hybridization and repeated backcrossing. This process is believed to be common in glacial refugia, where range expansions and contractions can lead to cycles of sympatry and isolation, creating conditions for extensive hybridization and introgression. Polyploidization is another genome-wide process with a major influence on plant evolution. Both hybridization and polyploidization can have complex effects on plant evolution. However, these effects are often difficult to understand in recently evolved species complexes. METHODS: We combined flow cytometry, analyses of transcriptomic sequences and pollen tube growth assays to investigate the consequences of polyploidization, hybridization and introgression on the recent evolution of several Erysimum (Brassicaceae) species from the South of the Iberian Peninsula, a well-known glacial refugium. This species complex differentiated in the last 2 million years, and its evolution has been hypothesized to be determined mainly by polyploidization, interspecific hybridization and introgression. KEY RESULTS: Our results support a scenario of widespread hybridization involving both extant and 'ghost' taxa. Several taxa studied here, most notably those with purple corollas, are polyploids, probably of allopolyploid origin. Moreover, hybridization in this group might be an ongoing phenomenon, as pre-zygotic barriers appeared weak in many cases. CONCLUSIONS: The evolution of Erysimum spp. has been determined by hybridization to a large extent. Species with purple (polyploids) and yellow flowers (mostly diploid) exhibit a strong signature of introgression in their genomes, indicating that hybridization occurred regardless of colour and across ploidy levels. Although the adaptive value of such genomic exchanges remains unclear, our results demonstrate the significance of hybridization for plant diversification, which should be taken into account when studying plant evolution.

Lack of ITS sequence homogenization in Erysimum species (Brassicaceae) with different ploidy levels.

The internal transcribed spacers (ITS) exhibit concerted evolution by the fast homogenization of these sequences at the intragenomic level. However, the rate and extension of this process are unclear and might be conditioned by the number and divergence of the different ITS copies. In some cases, such as hybrid species and polyploids, ITS sequence homogenization appears incomplete, resulting in multiple haplotypes within the same organism. Here, we studied the dynamics of concerted evolution in 85 individuals of seven plant species of the genus Erysimum (Brassicaceae) with multiple ploidy levels. We estimated the rate of concerted evolution and the degree of sequence homogenization separately for ITS1 and ITS2 and whether these varied with ploidy. Our results showed incomplete sequence homogenization, especially for polyploid samples, indicating a lack of concerted evolution in these taxa. Homogenization was usually higher in ITS2 than in ITS1, suggesting that concerted evolution operates more efficiently on the former. Furthermore, the hybrid origin of several species appears to contribute to the maintenance of high haplotype diversity, regardless of the level of ploidy. These findings indicate that sequence homogenization of ITS is a dynamic and complex process that might result in varying intra- and inter-genomic diversity levels.

Satellitome comparison of two oedipodine grasshoppers highlights the contingent nature of satellite DNA evolution.

BACKGROUND: The full catalog of satellite DNA (satDNA) within a same genome constitutes the satellitome. The Library Hypothesis predicts that satDNA in relative species reflects that in their common ancestor, but the evolutionary mechanisms and pathways of satDNA evolution have never been analyzed for full satellitomes. We compare here the satellitomes of two Oedipodine grasshoppers (Locusta migratoria and Oedaleus decorus) which shared their most recent common ancestor about 22.8 Ma ago. RESULTS: We found that about one third of their satDNA families (near 60 in every species) showed sequence homology and were grouped into 12 orthologous superfamilies. The turnover rate of consensus sequences was extremely variable among the 20 orthologous family pairs analyzed in both species. The satDNAs shared by both species showed poor association with sequence signatures and motives frequently argued as functional, except for short inverted repeats allowing short dyad symmetries and non-B DNA conformations. Orthologous satDNAs frequently showed different FISH patterns at both intra- and interspecific levels. We defined indices of homogenization and degeneration and quantified the level of incomplete library sorting between species. CONCLUSIONS: Our analyses revealed that satDNA degenerates through point mutation and homogenizes through partial turnovers caused by massive tandem duplications (the so-called satDNA amplification). Remarkably, satDNA amplification increases homogenization, at intragenomic level, and diversification between species, thus constituting the basis for concerted evolution. We suggest a model of satDNA evolution by means of recursive cycles of amplification and degeneration, leading to mostly contingent evolutionary pathways where concerted evolution emerges promptly after lineages split.

Phenotypic plasticity guides Moricandia arvensis divergence and convergence across the Brassicaceae floral morphospace.

Many flowers exhibit phenotypic plasticity. By inducing the production of several phenotypes, plasticity may favour the rapid exploration of different regions of the floral morphospace. We investigated how plasticity drives Moricandia arvensis, a species displaying within-individual floral polyphenism, across the floral morphospace of the entire Brassicaceae family. We compiled the multidimensional floral phenotype, the phylogenetic relationships, and the pollination niche of over 3000 species to construct a family-wide floral morphospace. We assessed the disparity between the two M. arvensis floral morphs (as the distance between the phenotypic spaces occupied by each morph) and compared it with the family-wide disparity. We measured floral divergence by comparing disparity with the most common ancestor, and estimated the convergence of each floral morph with other species belonging to the same pollination niches. Moricandia arvensis exhibits a plasticity-mediated floral disparity greater than that found between species, genera and tribes. The novel phenotype of M. arvensis moves outside the region occupied by its ancestors and relatives, crosses into a new region where it encounters a different pollination niche, and converges with distant Brassicaceae lineages. Our study suggests that phenotypic plasticity favours floral divergence and rapid appearance of convergent flowers, a process which facilitates the evolution of generalist pollination systems.

Asymmetric Reproductive Barriers and Gene Flow Promote the Rise of a Stable Hybrid Zone in the Mediterranean High Mountain.

Hybrid zones have the potential to shed light on evolutionary processes driving adaptation and speciation. Secondary contact hybrid zones are particularly powerful natural systems for studying the interaction between divergent genomes to understand the mode and rate at which reproductive isolation accumulates during speciation. We have studied a total of 720 plants belonging to five populations from two Erysimum (Brassicaceae) species presenting a contact zone in the Sierra Nevada mountains (SE Spain). The plants were phenotyped in 2007 and 2017, and most of them were genotyped the first year using 10 microsatellite markers. Plants coming from natural populations were grown in a common garden to evaluate the reproductive barriers between both species by means of controlled crosses. All the plants used for the field and greenhouse study were characterized by measuring traits related to plant size and flower size. We estimated the genetic molecular variances, the genetic differentiation, and the genetic structure by means of the F-statistic and Bayesian inference. We also estimated the amount of recent gene flow between populations. We found a narrow unimodal hybrid zone where the hybrid genotypes appear to have been maintained by significant levels of a unidirectional gene flow coming from parental populations and from weak reproductive isolation between them. Hybrid plants exhibited intermediate or vigorous phenotypes depending on the analyzed trait. The phenotypic differences between the hybrid and the parental plants were highly coherent between the field and controlled cross experiments and through time. The highly coherent results obtained by combining field, experimental, and genetic data demonstrate the existence of a stable and narrow unimodal hybrid zone between Erysimum mediohispanicum and Erysimum nevadense at the high elevation of the Sierra Nevada mountains.

Within-individual phenotypic plasticity in flowers fosters pollination niche shift.

Phenotypic plasticity, the ability of a genotype of producing different phenotypes when exposed to different environments, may impact ecological interactions. We study here how within-individual plasticity in Moricandia arvensis flowers modifies its pollination niche. During spring, this plant produces large, cross-shaped, UV-reflecting lilac flowers attracting mostly long-tongued large bees. However, unlike most co-occurring species, M. arvensis keeps flowering during the hot, dry summer due to its plasticity in key vegetative traits. Changes in temperature and photoperiod in summer trigger changes in gene expression and the production of small, rounded, UV-absorbing white flowers that attract a different assemblage of generalist pollinators. This shift in pollination niche potentially allows successful reproduction in harsh conditions, facilitating M. arvensis to face anthropogenic perturbations and climate change.

Promoting Individual and Collective Creativity in Science Students.

Creativity is a scientific skill necessary to develop a successful research career. We expose the importance of a growth mindset, divergent, lateral, and associative thinking, serendipity, and being part of a nonhierarchical and diverse research team to improve both individual and collective creativity.

Independent evolution of ancestral and novel defenses in a genus of toxic plants (Erysimum, Brassicaceae).

Phytochemical diversity is thought to result from coevolutionary cycles as specialization in herbivores imposes diversifying selection on plant chemical defenses. Plants in the speciose genus Erysimum (Brassicaceae) produce both ancestral glucosinolates and evolutionarily novel cardenolides as defenses. Here we test macroevolutionary hypotheses on co-expression, co-regulation, and diversification of these potentially redundant defenses across this genus. We sequenced and assembled the genome of E. cheiranthoides and foliar transcriptomes of 47 additional Erysimum species to construct a phylogeny from 9868 orthologous genes, revealing several geographic clades but also high levels of gene discordance. Concentrations, inducibility, and diversity of the two defenses varied independently among species, with no evidence for trade-offs. Closely related, geographically co-occurring species shared similar cardenolide traits, but not glucosinolate traits, likely as a result of specific selective pressures acting on each defense. Ancestral and novel chemical defenses in Erysimum thus appear to provide complementary rather than redundant functions.

Plants are often attacked by insects and other herbivores. As a result, they have evolved to defend themselves by producing many different chemicals that are toxic to these pests. As producing each chemical costs energy, individual plants often only produce one type of chemical that is targeted towards their main herbivore. Related species of plants often use the same type of chemical defense so, if a particular herbivore gains the ability to cope with this chemical, it may rapidly become an important pest for the whole plant family. To escape this threat, some plants have gained the ability to produce more than one type of chemical defense. Wallflowers, for example, are a group of plants in the mustard family that produce two types of toxic chemicals: mustard oils, which are common in most plants in this family; and cardenolides, which are an innovation of the wallflowers, and which are otherwise found only in distantly related plants such as foxglove and milkweed. The combination of these two chemical defenses within the same plant may have allowed the wallflowers to escape attacks from their main herbivores and may explain why the number of wallflower species rapidly increased within the last two million years. Züst et al. have now studied the diversity of mustard oils and cardenolides present in many different species of wallflower. This analysis revealed that almost all of the tested wallflower species produced high amounts of both chemical defenses, while only one species lacked the ability to produce cardenolides. The levels of mustard oils had no relation to the levels of cardenolides in the tested species, which suggests that the regulation of these two defenses is not linked. Furthermore, Züst et al. found that closely related wallflower species produced more similar cardenolides, but less similar mustard oils, to each other. This suggests that mustard oils and cardenolides have evolved independently in wallflowers and have distinct roles in the defense against different herbivores. The evolution of insect resistance to pesticides and other toxins is an important concern for agriculture. Applying multiple toxins to crops at the same time is an important strategy to slow the evolution of resistance in the pests. The findings of Züst et al. describe a system in which plants have naturally evolved an equivalent strategy to escape their main herbivores. Understanding how plants produce multiple chemical defenses, and the costs involved, may help efforts to breed crop species that are more resistant to herbivores and require fewer applications of pesticides.

Interpopulation spread of a parasitic B chromosome is unlikely through males in the grasshopper Eyprepocnemis plorans.

The near-neutral model of B chromosome evolution predicts that population invasion is quite fast. To test this prediction, in 1994, we introduced males of the grasshopper Eyprepocnemis plorans from a B-carrying population into a B-lacking population and monitored the evolution of B-chromosome frequency up to 2013. We observed fluctuating very low B frequency across years but, remarkably, the B chromosome introduced (the B2 variant) was found up to 1996 only, whereas the B1 variant was present from 1996 onwards, presumably introduced by fishermen using E. plorans males as bait. Effective introgression of genetic material from the donor population was evidenced by the presence of a satellite DNA on autosome 9 (up to 1999) and the presence of one individual in 2006 showing an ISSR marker profile being highly similar to that found in the donor population. This indicated that the males introduced by us effectively mated with resident females, but donor genes rapidly decreased in frequency after this non-recurrent migration event. Taken together, our results indicated: (i) that the non-recurrent migration event had a slight, transient genetic effect on the recipient population, which was diluted in only a few generations; and (ii) that even with recurrent migration (forced by fishermen) the B chromosome failed to increase in frequency. Bearing in mind that B chromosomes in this species drive through females only, we hypothesize that B chromosomes most likely failed invasion in both migration events because the migrating sex shows no B-drive.

Pollination effectiveness in a generalist plant: adding the genetic component.

The pollination effectiveness of a flower visitor has traditionally been measured as the product of a quantity component that depends on the frequency of interaction and a quality component that measures the per-visit effects on plant reproduction. We propose that this could be complemented with a genetic component informing about each pollinator's contribution to the genetic diversity and composition of the plant progeny. We measured the quantity and quality components of effectiveness of most pollinator functional groups of the generalist herb Erysimum mediohispanicum. We used 10 microsatellite markers to calculate the genetic component as the diversity of sires among siblings and included it into the calculation of the pollination effectiveness. Functional groups varied in the quantity and quality components, which were shown to be decoupled. Functional groups also differed in the genetic component. This component changed the estimates of pollination effectiveness, increasing the differences between some functional groups and modifying the pollination effectiveness landscape. We demonstrate that including the genetic component in the calculation of the pollination effectiveness may allow a more complete quantification of the contribution of each pollinator to the reproductive success of a plant, providing information on its mating patterns and long-term fitness.

Anther Rubbing, a New Mechanism That Actively Promotes Selfing in Plants.

Self-fertilization has recurrently evolved in plants, involving different strategies and traits and often loss of attractive functions, collectively known as the selfing syndrome. However, few traits that actively promote self-fertilization have been described. Here we describe a novel mechanism promoting self-fertilization in the Brassicaceae species Erysimum incanum. This mechanism, which we called "anther rubbing," consists of autonomous, repeated, and coordinated movements of the stamens over the stigma during flower opening. We have documented anther rubbing by time-lapse videos and experimentally show that it causes self-pollen deposition on stigmas and is sufficient to achieve maximal reproductive output in E. incanum. We predict that these movements should occur in species with limited inbreeding depression, and indeed we find that inbreeding depression in seed production is negligible in this species. While many studies have documented complex floral traits that promote outcrossing, the occurrence of anther rubbing demonstrates that plants can evolve elaborate and underappreciated adaptations to promote self-fertilization.

Comparative assessment shows the reliability of chloroplast genome assembly using RNA-seq.

Chloroplast genomes (cp genomes) are widely used in comparative genomics, population genetics, and phylogenetic studies. Obtaining chloroplast genomes from RNA-Seq data seems feasible due to the almost full transcription of cpDNA. However, the reliability of chloroplast genomes assembled from RNA-Seq instead of genomic DNA libraries remains to be thoroughly verified. In this study, we assembled chloroplast genomes for three Erysimum (Brassicaceae) species from three RNA-Seq replicas and from one genomic library of each species, using a streamlined bioinformatics protocol. We compared these assembled genomes, confirming that assembled cp genomes from RNA-Seq data were highly similar to each other and to those from genomic libraries in terms of overall structure, size, and composition. Although post-transcriptional modifications, such as RNA-editing, may introduce variations in the RNA-seq data, the assembly of cp genomes from RNA-seq appeared to be reliable. Moreover, RNA-Seq assembly was less sensitive to sources of error such as the recovery of nuclear plastid DNAs (NUPTs). Although some precautions should be taken when producing reference genomes in non-model plants, we conclude that assembling cp genomes from RNA-Seq data is a fast, accurate, and reliable strategy.

Long-term monitoring of B-chromosome invasion and neutralization in a population of Prospero autumnale (Asparagaceae).

B chromosomes have been reported in about 15% of eukaryotes, but long-term dynamics of B chromosomes in a single natural population has rarely been analyzed. Prospero autumnale plants collected in 1981 and 1983 at Cuesta de La Palma population had shown the presence of B chromosomes. We analyze here seven additional samples collected between 1987 and 2015, and show that B frequency increased significantly during the 1980s and showed minor fluctuations between 2005 and 2015. A mother-offspring analysis of B chromosome transmission, at population level, showed significant drive on the male side (kB  = 0.65) and significant drag on the female side (kB  = 0.33), with average B transmission rate being very close to the Mendelian rate (0.5). No significant effects of B chromosomes were observed on a number of vigor and fertility-related traits. Within a parasite/host framework, these results suggest that B chromosomes' drive on the male side is the main pathway for B chromosome invasion, whereas B chromosome drag on the female side might be the main manifestation of host genome resistance in this species. Prospero autumnale thus illuminates a novel evolutionary pathway for B chromosome neutralization by means of a decrease in B transmission through the nondriving sex.

Molecular phylogeny and evolutionary history of Moricandia DC (Brassicaceae).

BACKGROUND: The phylogeny of tribe Brassiceae (Brassicaceae) has not yet been resolved because of its complex evolutionary history. This tribe comprises economically relevant species, including the genus Moricandia DC. This genus is currently distributed in North Africa, Middle East, Central Asia and Southern Europe, where it is associated with arid and semi-arid environments. Although some species of Moricandia have been used in several phylogenetic studies, the phylogeny of this genus is not well established. METHODS: Here we present a phylogenetic analysis of the genus Moricandia using a nuclear (the internal transcribed spacers of the ribosomal DNA) and two plastidial regions (parts of the NADH dehydrogenase subunit F gene and the trnT-trnF region). We also included in the analyses members of their sister genus Rytidocarpus and from the close genus Eruca. RESULTS: The phylogenetic analyses showed a clear and robust phylogeny of the genus Moricandia. The Bayesian inference tree was concordant with the maximum likelihood and timing trees, with the plastidial and nuclear trees showing only minor discrepancies. The genus Moricandia appears to be formed by two main lineages: the Iberian clade including three species, and the African clade including the four species inhabiting the Southern Mediterranean regions plus M. arvensis. DISCUSSION: We dated the main evolutionary events of this genus, showing that the origin of the Iberian clade probably occurred after a range expansion during the Messinian period, between 7.25 and 5.33 Ma. In that period, an extensive African-Iberian floral and faunal interchange occurred due to the existence of land bridges between Africa and Europa in what is, at present-days, the Strait of Gibraltar. We have demonstrated that a Spanish population previously ascribed to Rytidocarpus moricandioides is indeed a Moricandia species, and we propose to name it as M. rytidocarpoides sp. nov. In addition, in all the phylogenetic analyses, M. foleyi appeared outside the Moricandia lineage but within the genus Eruca. Therefore, M. foleyi should be excluded from the genus Moricandia and be ascribed, at least provisionally, to the genus Eruca.

Timing of the evolutionary history of Corallinaceae (Corallinales, Rhodophyta).

The temporal dimension of the most recent Corallinaceae (order Corallinales) phylogeny was presented here, based on first occurrence time estimates from the fossil record. Calibration of the molecular clock of the genetic marker SSU entailed a separation of Corallinales from Hapalidiales in the Albian (Early Cretaceous ~105 mya). Neither the calibration nor the fossil record resolved the succession of appearance of the first three emerging subfamilies: Mastophoroideae, Corallinoideae, and Neogoniolithoideae. The development of the tetra/bisporangial conceptacle roofs by filaments surrounding and interspersed among the sporangial initials was an evolutionary novelty emerging at the Cretaceous-Paleogene boundary (~66 mya). This novelty was shared by the subfamilies Hydrolithoideae, Metagoniolithoideae, and Lithophylloideae, which diverged in the early Paleogene. Subclades within the Metagoniolithoideae and Lithophylloideae diversified in the late Oligocene-middle Miocene (~28-12 mya). The most common reef corallinaceans (Hydrolithon, Porolithon, Harveylithon, "Pneophyllum" conicum, and subclades within Lithophylloideae) appeared in this interval in the Indo-Australian Archipelago.

Inter-annual maintenance of the fine-scale genetic structure in a biennial plant.

Within plant populations, space-restricted gene movement, together with environmental heterogeneity, can result in a spatial variation in gene frequencies. In biennial plants, inter-annual flowering migrants can homogenize gene frequencies between consecutive cohorts. However, the actual impact of these migrants on spatial genetic variation remains unexplored. Here, we used 10 nuclear microsatellite and one plastid genetic marker to characterize the spatial genetic structure within two consecutive cohorts in a population of the biennial plant Erysimum mediohispanicum (Brassicaceae). We explored the maintenance of this structure between consecutive flowering cohorts at different levels of complexity, and investigated landscape effects on gene flow. We found that cohorts were not genetically differentiated and showed a spatial genetic structure defined by a negative genetic-spatial correlation at fine scale that varied in intensity with compass directions. This spatial genetic structure was maintained when comparing plants from different cohorts. Additionally, genotypes were consistently associated with environmental factors such as light availability and soil composition, but to a lesser extent compared with the spatial autocorrelation. We conclude that inter-annual migrants, in combination with limited seed dispersal and environmental heterogeneity, play a major role in shaping and maintaining the spatial genetic structure among cohorts in this biennial plant.

Phylogenetic relationships of corallinaceae (Corallinales, Rhodophyta): taxonomic implications for reef-building corallines.

A new, more complete, five-marker (SSU, LSU, psbA, COI, 23S) molecular phylogeny of the family Corallinaceae, order Corallinales, shows a paraphyletic grouping of seven well-supported monophyletic clades. The taxonomic implications included the amendment of two subfamilies, Neogoniolithoideae and Metagoniolithoideae, and the rejection of Porolithoideae as an independent subfamily. Metagoniolithoideae contained Harveylithon gen. nov., with H. rupestre comb. nov. as the generitype, and H. canariense stat. nov., H. munitum comb. nov., and H. samoënse comb. nov. Spongites and Pneophyllum belonged to separate clades. The subfamily Neogoniolithoideae included the generitype of Spongites, S. fruticulosus, for which an epitype was designated. Pneophyllum requires reassesment. The generitype of Hydrolithon, H. reinboldii, was a younger heterotypic synonym of H. boergesenii. The evolutionary novelty of the subfamilies Hydrolithoideae, Metagoniolithoideae, and Lithophylloideae was the development of tetra/bisporangial conceptacle roofs by filaments surrounding and interspersed among the sporangial initials.

The role of pollinators in the evolution of corolla shape variation, disparity and integration in a highly diversified plant family with a conserved floral bauplan.

BACKGROUND AND AIMS: Brassicaceae is one of the most diversified families in the angiosperms. However, most species from this family exhibit a very similar floral bauplan. In this study, we explore the Brassicaceae floral morphospace, examining how corolla shape variation (an estimation of developmental robustness), integration and disparity vary among phylogenetically related species. Our aim is to check whether these floral attributes have evolved in this family despite its apparent morphological conservation, and to test the role of pollinators in driving this evolution. METHODS: Using geometric morphometric tools, we calculated the phenotypic variation, disparity and integration of the corolla shape of 111 Brassicaceae taxa. We subsequently inferred the phylogenetic relationships of these taxa and explored the evolutionary lability of corolla shape. Finally, we sampled the pollinator assemblages of every taxon included in this study, and determined their pollination niches using a modularity algorithm. We explore the relationship between pollination niche and the attributes of corolla shape. KEY RESULTS: Phylogenetic signal was weak for all corolla shape attributes. All taxa had generalized pollination systems. Nevertheless, they belong to different pollination niches. There were significant differences in corolla shape among pollination niches even after controlling for the phylogenetic relationship of the plant taxa. Corolla shape variation and disparity was significantly higher in those taxa visited mostly by nocturnal moths, indicating that this pollination niche is associated with a lack of developmental robustness. Corolla integration was higher in those taxa visited mostly by hovering long-tongued flies and long-tongued large bees. CONCLUSIONS: Corolla variation, integration and disparity were evolutionarily labile and evolved very recently in the evolutionary history of the Brassicaceae. These floral attributes were strongly related to the pollination niche. Even in a plant clade having a very generalized pollination system and exhibiting a conserved floral bauplan, pollinators can drive the evolution of important developmental attributes of corolla shape.