The Phenotypic and Genetic Underpinnings of Flower Size in Polemoniaceae

Towards a global perspective for Salvia L.: Phylogeny, diversification and floral evolution

Salvia is the most species-rich genus in Lamiaceae, encompassing approximately 1000 species distributed all over the world. We sought a new evolutionary perspective for Salvia by employing macroevolutionary analyses to address the tempo and mode of diversification. To study the association of floral traits with speciation and extinction, we modelled and explored the evolution of corolla length and the lever-mechanism pollination system across our Salvia phylogeny. We reconstructed a multigene phylogeny for 366 species of Salvia in the broad sense including all major recognized lineages and 50 species from Iran, a region previously overlooked in studies of the genus. Our comprehensive sampling of Iranian species of Salvia provides higher phylogenetic resolution for southwestern Asian species than obtained in previous studies. Our phylogenetic data in combination with divergence time estimates were used to examine the evolution of corolla length, woody versus herbaceous habit, and presence versus absence of a lever mechanism. We investigated the timing and dependence of Salvia diversification related to corolla length evolution through a disparity test and BAMM analysis. A HiSSE model was used to evaluate the dependency of diversification on the lever-mechanism pollination system in Salvia. A medium corolla length (15-18 mm) was reconstructed as the ancestral state for Salvia with multiple shifts to shorter and longer corollas. Macroevolutionary model analyses indicate that corolla length disparity is high throughout Salvia evolution, significantly different from expectations under a Brownian motion model during the last 28 million years of evolution. Our analyses show evidence of a higher diversification rate of corolla length for some Andean species of Salvia compared to other members of the genus. Based on our tests of diversification models, we reject the hypothesis of a direct effect of the lever mechanism on Salvia diversification. Therefore, we suggest caution in considering the lever-mechanism pollination system as one of the main drivers of speciation in Salvia.

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.

Parallelism in Flower Evolution and Development

Flower evolution is characterized by widespread repetition, with adaptations to pollinator environment evolving in parallel. Recent studies have expanded our understanding of the developmental basis of adaptive floral novelties—petal fusion, bilateral symmetry, heterostyly, and floral dimensions. In this article, we describe patterns of trait evolution and review developmental genetic mechanisms underlying floral novelties. We discuss the diversity of mechanisms for parallel adaptation, the evidence for constraints on these mechanisms, and how constraints help explain observed macroevolutionary patterns. We describe parallel evolution resulting from similarities at multiple hierarchical levels—genetic, developmental, morphological, functional—which indicate general principles in floral evolution, including the central role of hormone signaling. An emerging pattern is mutational bias that may contribute to rapid patterns of parallel evolution, especially if the derived trait can result from simple degenerative mutations. We argue that such mutational bias may be less likely to govern the evolution of novelties patterned by complex developmental pathways.

Differential Gene Expression with an Emphasis on Floral Organ Size Differences in Natural and Synthetic Polyploids of Nicotiana tabacum (Solanaceae)

Floral organ size, especially the size of the corolla, plays an important role in plant reproduction by facilitating pollination efficiency. Previous studies have outlined a hypothesized organ size pathway. However, the expression and function of many of the genes in the pathway have only been investigated in model diploid species; therefore, it is unknown how these genes interact in polyploid species. Although correlations between ploidy and cell size have been shown in many systems, it is unclear whether there is a difference in cell size between naturally occurring and synthetic polyploids. To address these questions comparing floral organ size and cell size across ploidy, we use natural and synthetic polyploids of Nicotiana tabacum (Solanaceae) as well as their known diploid progenitors. We employ a comparative transcriptomics approach to perform analyses of differential gene expression, focusing on candidate genes that may be involved in floral organ size, both across developmental stages and across accessions. We see differential expression of several known floral organ candidate genes including ARF2, BIG BROTHER, and GASA/GAST1. Results from linear models show that ploidy, cell width, and cell number positively influence corolla tube circumference; however, the effect of cell width varies by ploidy, and diploids have a significantly steeper slope than both natural and synthetic polyploids. These results demonstrate that polyploids have wider cells and that polyploidy significantly increases corolla tube circumference.

Out of Sight, Out of Mind: Widespread Nuclear and Plastid-Nuclear Discordance in the Flowering Plant Genus Polemonium (Polemoniaceae) Suggests Widespread Historical Gene Flow Despite Limited Nuclear Signal

Phylogenomic data from a rapidly increasing number of studies provide new evidence for resolving relationships in recently radiated clades, but they also pose new challenges for inferring evolutionary histories. Most existing methods for reconstructing phylogenetic hypotheses rely solely on algorithms that only consider incomplete lineage sorting (ILS) as a cause of intra- or intergenomic discordance. Here, we utilize a variety of methods, including those to infer phylogenetic networks, to account for both ILS and introgression as a cause for nuclear and cytoplasmic-nuclear discordance using phylogenomic data from the recently radiated flowering plant genus Polemonium (Polemoniaceae), an ecologically diverse genus in Western North America with known and suspected gene flow between species. We find evidence for widespread discordance among nuclear loci that can be explained by both ILS and reticulate evolution in the evolutionary history of Polemonium. Furthermore, the histories of organellar genomes show strong discordance with the inferred species tree from the nuclear genome. Discordance between the nuclear and plastid genome is not completely explained by ILS, and only one case of discordance is explained by detected introgression events. Our results suggest that multiple processes have been involved in the evolutionary history of Polemonium and that the plastid genome does not accurately reflect species relationships. We discuss several potential causes for this cytoplasmic-nuclear discordance, which emerging evidence suggests is more widespread across the Tree of Life than previously thought. [Cyto-nuclear discordance, genomic discordance, phylogenetic networks, plastid capture, Polemoniaceae, Polemonium, reticulations.].

A two-tier bioinformatic pipeline to develop probes for target capture of nuclear loci with applications in Melastomataceae

PREMISE

Putatively single-copy nuclear (SCN) loci, which are identified using genomic resources of closely related species, are ideal for phylogenomic inference. However, suitable genomic resources are not available for many clades, including Melastomataceae. We introduce a versatile approach to identify SCN loci for clades with few genomic resources and use it to develop probes for target enrichment in the distantly related Memecylon and Tibouchina (Melastomataceae).

METHODS

We present a two-tiered pipeline. First, we identified putatively SCN loci using MarkerMiner and transcriptomes from distantly related species in Melastomataceae. Published loci and genes of functional significance were then added (384 total loci). Second, using HybPiper, we retrieved 689 homologous template sequences for these loci using genome-skimming data from within the focal clades.

RESULTS

We sequenced 193 loci common to Memecylon and Tibouchina. Probes designed from 56 template sequences successfully targeted sequences in both clades. Probes designed from genome-skimming data within a focal clade were more successful than probes designed from other sources.

DISCUSSION

Our pipeline successfully identified and targeted SCN loci in Memecylon and Tibouchina, enabling phylogenomic studies in both clades and potentially across Melastomataceae. This pipeline could be easily applied to other clades with few genomic resources.

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

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

Developmental mechanisms involved in the diversification of flowers

We all appreciate the fantastic diversity of flowers. How flowers diversified, however, remains largely enigmatic. The mechanisms underlying the diversification of flowers are complex because the overall appearance of a flower is determined by many factors, such as the shape and size of its receptacle, and the arrangement, number, type, shape and colour of floral organs. Modifications of the developmental trajectories of a flower and its components, therefore, can lead to the generation of new floral types. In this Review, by summarizing the recent progress in studying the initiation, identity determination, morphogenesis and maturation of floral organs, we present our current understanding of the mechanisms underlying the diversification of flowers.

Orthologous nuclear markers and new transcriptomes that broadly cover the phylogenetic diversity of Acanthaceae

PREMISE

Information on orthologous groups of genes, their sequence variability, and annotation is required for project design in phylogenetic reconstruction. This resource is unavailable for the flowering plant family Acanthaceae (>4000 species).

METHODS

We compared transcriptome sequences spanning the extant diversity of Acanthaceae in order to provide a set of orthologous low-copy nuclear genes and assess their utility for reconstructing phylogenetic relationships within this group of plants.

RESULTS

We present new transcriptome assemblies for eight species representing all major clades of Acanthaceae. The assemblies of five of these species are entirely based on new sequence data. Of these five species, three are from subfamilies for which no genomic resources were previously available (Nelsonioideae and Thunbergioideae). These five new transcriptomes are more complete than all others from public databases. Furthermore, we provide alignments with sequence information, annotation, and statistics for potential phylogenetic utility of 1619 orthologous low-copy nuclear markers.

DISCUSSION

Our method of inferring assemblies from multiple pooled tissue samples delivers more complete transcriptomes than any available ones from Acanthaceae. We make available to the community new resources (e.g., sequence information, variability, and annotation of orthologous low-copy nuclear genes) that will help phylogenetic reconstruction in Acanthaceae.

A Universal Probe Set for Targeted Sequencing of 353 Nuclear Genes from Any Flowering Plant Designed Using k-Medoids Clustering

Sequencing of target-enriched libraries is an efficient and cost-effective method for obtaining DNA sequence data from hundreds of nuclear loci for phylogeny reconstruction. Much of the cost of developing targeted sequencing approaches is associated with the generation of preliminary data needed for the identification of orthologous loci for probe design. In plants, identifying orthologous loci has proven difficult due to a large number of whole-genome duplication events, especially in the angiosperms (flowering plants). We used multiple sequence alignments from over 600 angiosperms for 353 putatively single-copy protein-coding genes identified by the One Thousand Plant Transcriptomes Initiative to design a set of targeted sequencing probes for phylogenetic studies of any angiosperm group. To maximize the phylogenetic potential of the probes, while minimizing the cost of production, we introduce a k-medoids clustering approach to identify the minimum number of sequences necessary to represent each coding sequence in the final probe set. Using this method, 5-15 representative sequences were selected per orthologous locus, representing the sequence diversity of angiosperms more efficiently than if probes were designed using available sequenced genomes alone. To test our approximately 80,000 probes, we hybridized libraries from 42 species spanning all higher-order groups of angiosperms, with a focus on taxa not present in the sequence alignments used to design the probes. Out of a possible 353 coding sequences, we recovered an average of 283 per species and at least 100 in all species. Differences among taxa in sequence recovery could not be explained by relatedness to the representative taxa selected for probe design, suggesting that there is no phylogenetic bias in the probe set. Our probe set, which targeted 260 kbp of coding sequence, achieved a median recovery of 137 kbp per taxon in coding regions, a maximum recovery of 250 kbp, and an additional median of 212 kbp per taxon in flanking non-coding regions across all species. These results suggest that the Angiosperms353 probe set described here is effective for any group of flowering plants and would be useful for phylogenetic studies from the species level to higher-order groups, including the entire angiosperm clade itself.

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.

Contrasting lengths of Pelargonium floral nectar tubes result from late differences in rate and duration of growth

BACKGROUND AND AIMS

Much of morphological evolution in flowers has arisen from pollinator-mediated selection, often manifest as a match between the length of the pollinator's proboscis and the depth of tubular corollas or spurs. We investigate development, growth and homology of the unique nectar tube of Pelargonium, frequently described as 'a spur adnate to the pedicel'.

METHODS

We focused on two species. The nectar tube of P. ionidiflorum is three times longer than that of P. odoratissimum. Light and scanning electron microscopy were carried out, and daily growth measurements were used to compare nectar tube development and vascular patterns.

KEY RESULTS

Nectar tubes in both species are initiated centripetally to the dorsal sepal in a space created by lateral displacement of two antepetalous stamens. The cavity deepens through subsequent intercalary growth of the receptacle that proceeds at the same rate in both species until tubes reach approx. 10 mm in length. Differences in final nectar tube lengths arise via an increase in the rate and duration of growth of the receptacle that begins just before anthesis (floral opening) and continues for several days past anthesis in P. ionidiflorum but does not occur in P. odoratissimum. Epidermal cells of the dorsal surface of the nectar tube in P. ionidiflorum are approx. 1.6 times longer than those in P. odoratissimum. Histological sections show no evidence that the nectar tube is a spur that became evolutionarily fused to the pedicel.

CONCLUSIONS

Nectar tubes in Pelargonium are localized cavities that form in the receptacle via intercalary growth. Differences in the rate and duration of growth just prior to and following anthesis underlie differences in final tube lengths. Because differences in cell lengths do not fully account for differences in nectar tube lengths, evolutionary diversification must involve changes in both cell cycle and cell expansion.

Comparative transcriptomic analysis of the evolution and development of flower size in Saltugilia (Polemoniaceae)

BACKGROUND

Flower size varies dramatically across angiosperms, representing innovations over the course of >130 million years of evolution and contributing substantially to relationships with pollinators. However, the genetic underpinning of flower size is not well understood. Saltugilia (Polemoniaceae) provides an excellent non-model system for extending the genetic study of flower size to interspecific differences that coincide with variation in pollinators.

RESULTS

Using targeted gene capture methods, we infer phylogenetic relationships among all members of Saltugilia to provide a framework for investigating the genetic control of flower size differences via RNA-Seq de novo assembly. Nuclear concatenation and species tree inference methods provide congruent topologies. The inferred evolutionary trajectory of flower size is from small flowers to larger flowers. We identified 4 to 10,368 transcripts that are differentially expressed during flower development, with many unigenes associated with cell wall modification and components of the auxin and gibberellin pathways.

CONCLUSIONS

Saltugilia is an excellent model for investigating covarying floral and pollinator evolution. Four candidate genes from model systems (BIG BROTHER, BIG PETAL, GASA, and LONGIFOLIA) show differential expression during development of flowers in Saltugilia, and four other genes (FLOWERING-PROMOTING FACTOR 1, PECTINESTERASE, POLYGALACTURONASE, and SUCROSE SYNTHASE) fit into hypothesized organ size pathways. Together, these gene sets provide a strong foundation for future functional studies to determine their roles in specifying interspecific differences in flower size.

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

A third of all angiosperm species produce flowers with petals fused into a corolla tube. The various elaborations of corolla tube attributes, such as length, width and curvature, have enabled plants to exploit many specialized pollinator groups. These elaborations often differ dramatically among closely related species, contributing to pollinator shift and pollinator-mediated reproductive isolation and speciation. However, very little is known about the genetic and developmental control of these corolla tube attributes. Here we report the characterization of a semi-dominant mutant in the monkeyflower species Mimulus lewisii, with a substantial decrease in corolla tube width but no change in tube length. This morphological alteration leads to a ˜ 70% decrease in bumblebee visitation rate for the homozygous mutant compared to the wild-type. Through bulk segregant analysis and transgenic experiment, we show that the mutant phenotype is caused by a dominant-negative mutation in an actin gene. This mutation decreases epidermal cell width but not length, and probably also reduces the number of lateral cell divisions. These results suggest a surprising potential role for a 'housekeeping' gene in fine-tuning the development of an ecologically important floral trait.

De novo Sequencing and Comparative Transcriptomics of Floral Development of the Distylous Species Lithospermum multiflorum

Genes controlling the morphological, micromorphological, and physiological components of the breeding system distyly have been hypothesized, but many of the genes have not been investigated throughout development of the two floral morphs. To this end, the present study is an examination of comparative transcriptomes from three stages of development for the floral organs of the morphs of Lithospermum multiflorum. Transcriptomes of flowers of the two morphs, from various stages of development, were sequenced using an Illumina HiSeq 2000. The floral transcriptome of L. multiflorum was assembled, and differential gene expression (DE) was identified between morphs, throughout development. Additionally, Gene Ontology (GO) terms for DE genes were determined. Fewer genes were DE early in development compared to later in development, with more genes highly expressed in the gynoecium of the SS morph and the corolla and androecium of the LS morph. A reciprocal pattern was observed later in development, and many more genes were DE during this latter stage. During early development, DE genes appear to be involved in growth and floral development, and during later development, DE genes seem to affect physiological functions. Interestingly, many genes involved in response to stress were identified as DE between morphs.

The report of my death was an exaggeration: A review for researchers using microsatellites in the 21st century

Microsatellites, or simple sequence repeats (SSRs), have long played a major role in genetic studies due to their typically high polymorphism. They have diverse applications, including genome mapping, forensics, ascertaining parentage, population and conservation genetics, identification of the parentage of polyploids, and phylogeography. We compare SSRs and newer methods, such as genotyping by sequencing (GBS) and restriction site associated DNA sequencing (RAD-Seq), and offer recommendations for researchers considering which genetic markers to use. We also review the variety of techniques currently used for identifying microsatellite loci and developing primers, with a particular focus on those that make use of next-generation sequencing (NGS). Additionally, we review software for microsatellite development and report on an experiment to assess the utility of currently available software for SSR development. Finally, we discuss the future of microsatellites and make recommendations for researchers preparing to use microsatellites. We argue that microsatellites still have an important place in the genomic age as they remain effective and cost-efficient markers.