Ecological correlates and genetic consequences of evolutionary transitions from distyly to homostyly

Ecological correlates of population genetics in Linum suffruticosum, an heterostylous polyploid and taxonomic complex endemic to the Western Mediterranean Basin

Linum suffruticosum s.l. is a taxonomic complex widespread in the Western Mediterranean basin. The complex is characterized by a high phenotypic and cytogenetic diversity, and by a unique three-dimensional heterostyly system that makes it an obligate outcrosser. We studied the patterns of genetic diversity and structure of populations throughout the entire distribution of L. suffruticosum s.l. with microsatellite markers. We analysed their relationships with various biological and ecological variables, including the morph ratio and sex organ reciprocity of populations measured with a novel multi-dimensional method. Populations consistently showed an approximate 1:1 morph ratio with high sex organ reciprocity and high genetic diversity. We found high genetic differentiation of populations, showing a pattern of isolation by distance. The Rif mountains in NW Africa were the most important genetic barrier. The taxonomic treatment within the group was not related to the genetic differentiation of populations, but to their environmental differentiation. Genetic diversity was unrelated to latitude, elevation, population size, niche suitability or breeding system. However, there was a clear influence of ploidy level on the genetic diversity of populations, and a seeming centre-periphery pattern in its distribution. Our results suggest that polyploidization events, high outcrossing rates, isolation by distance and important geographical barriers to gene flow have played major roles in the microevolutionary history of this species complex.

Development of polymorphic microsatellite markers for distylous-homostylous Primula secundiflora (Primulaceae) using HiSeq sequencing

Primula secundiflora is an insect-pollinated, perennial herb belonging to the section Proliferae (Primulaceae) that exhibits considerable variation in its mating system, with predominantly outcrossing populations comprising long-styled and short-styled floral morphs and selfing populations comprising only homostyles. To facilitate future investigations of the population genetics and mating patterns of this species, we developed 25 microsatellite markers from P. secundiflora using next-generation sequencing and measured polymorphism and genetic diversity in a sample of 30 individuals from three natural populations. The markers displayed high polymorphism, with the number of observed alleles per locus ranging from three to 16 (mean = 8.36). The observed and expected heterozygosities ranged from 0.100 to 1.000 and 0.145 to 0.843, respectively. Twenty-one of the loci were also successfully amplified in P. denticulata. These microsatellite markers should provide powerful tools for investigating patterns of population genetic diversity and the evolutionary relationships between distyly and homostyly in P. secundiflora.

Genetic Causes and Genomic Consequences of Breakdown of Distyly in Linum trigynum

Distyly is an iconic floral polymorphism governed by a supergene, which promotes efficient pollen transfer and outcrossing through reciprocal differences in the position of sexual organs in flowers, often coupled with heteromorphic self-incompatibility. Distyly has evolved convergently in multiple flowering plant lineages, but has also broken down repeatedly, often resulting in homostylous, self-compatible populations with elevated rates of self-fertilization. Here, we aimed to study the genetic causes and genomic consequences of the shift to homostyly in Linum trigynum, which is closely related to distylous Linum tenue. Building on a high-quality genome assembly, we show that L. trigynum harbors a genomic region homologous to the dominant haplotype of the distyly supergene conferring long stamens and short styles in L. tenue, suggesting that loss of distyly first occurred in a short-styled individual. In contrast to homostylous Primula and Fagopyrum, L. trigynum harbors no fixed loss-of-function mutations in coding sequences of S-linked distyly candidate genes. Instead, floral gene expression analyses and controlled crosses suggest that mutations downregulating the S-linked LtWDR-44 candidate gene for male self-incompatibility and/or anther height could underlie homostyly and self-compatibility in L. trigynum. Population genomic analyses of 224 whole-genome sequences further demonstrate that L. trigynum is highly self-fertilizing, exhibits significantly lower genetic diversity genome-wide, and is experiencing relaxed purifying selection and less frequent positive selection on nonsynonymous mutations relative to L. tenue. Our analyses shed light on the loss of distyly in L. trigynum, and advance our understanding of a common evolutionary transition in flowering plants.

Asymmetrical disassortative pollination mediated by long-/short-tongued pollinators in a distylous Limonium myrianthum (Plumbaginaceae) with a short corolla tubular small flower

In heterostylous plants, short-tongued pollinators are often ineffective/inefficient owing to the limitations imposed by a long corolla tube. However, it is unclear how disassortative pollen transfer is achieved in small flowers. We investigated the pollination pattern and floral morph variation by analyzing heterostylous syndrome, pollinator groups, and pollen deposition after a single visitation in two Limonium myrianthum populations with short-corolla-tubular small flowers. The predominant pollinators in the Hutubi population were pollen-seeking short-tongued syrphids, which can only transfer pollen between high-level sexual organs. In the Xishan population, nectar-seeking short-tongued insects were efficient pollinators with symmetrical disassortative pollen transfer between high- and low-level sexual organs, whereas long-tongued pollinators had a low efficiency between high-level sexual organs due to the low contact probability with the stigma of long-styled flowers (L-morph), which no longer offered the same advantage observed in tubular flowers. Asymmetrical disassortative pollination may cause the female fitness of short-styled (S-morph) individuals in the Hutubi and L-morph individuals in the Xishan population to suffer greater selection pressure and exhibit a higher degree of floral morph variation. Limonium myrianthum exhibits an unusual pollination pattern in which the small flowers with short corolla tubes make it possible for short-tongued insects to become effective pollinators. However, factors such as the position of stigma-anther within the flower, pollinator species and their preference further caused asymmetrical disassortative pollen transfer. Therefore, more factors should be considered when evaluating the effectiveness of short- and long-tongued insects in pollination service.

Parallel evolution of morphological and genomic selfing syndromes accompany the breakdown of heterostyly

Evolutionary transitions from outcrossing to selfing in flowering plants have convergent morphological and genomic signatures and can involve parallel evolution within related lineages. Adaptive evolution of morphological traits is often assumed to evolve faster than nonadaptive features of the genomic selfing syndrome. We investigated phenotypic and genomic changes associated with transitions from distyly to homostyly in the Primula oreodoxa complex. We determined whether the transition to selfing occurred more than once and investigated stages in the evolution of morphological and genomic selfing syndromes using 22 floral traits and both nuclear and plastid genomic data from 25 populations. Two independent transitions were detected representing an earlier and a more recently derived selfing lineage. The older lineage exhibited classic features of the morphological and genomic selfing syndrome. Although features of both selfing syndromes were less developed in the younger selfing lineage, they exhibited parallel development with the older selfing lineage. This finding contrasts with the prediction that some genomic changes should lag behind adaptive changes to morphological traits. Our findings highlight the value of comparative studies on the timing and extent of transitions from outcrossing to selfing between related lineages for investigating the tempo of morphological and molecular evolution.

Genomic analyses elucidate S-locus evolution in response to intra-specific losses of distyly in Primula vulgaris

Distyly, a floral dimorphism that promotes outcrossing, is controlled by a hemizygous genomic region known as the S-locus. Disruptions of genes within the S-locus are responsible for the loss of distyly and the emergence of homostyly, a floral monomorphism that favors selfing. Using whole-genome resequencing data of distylous and homostylous individuals from populations of Primula vulgaris and leveraging high-quality reference genomes of Primula we tested, for the first time, predictions about the evolutionary consequences of transitions to selfing on S-genes. Our results reveal a previously undetected structural rearrangement in CYPᵀ associated with the shift to homostyly and confirm previously reported, homostyle-specific, loss-of-function mutations in the exons of the S-gene CYPᵀ. We also discovered that the promoter and intronic regions of CYPᵀ in distylous and homostylous individuals are conserved, suggesting that down-regulation of CYPᵀ via mutations in its promoter and intronic regions is not a cause of the shift to homostyly. Furthermore, we found that hemizygosity is associated with reduced genetic diversity in S-genes compared with their paralogs outside the S-locus. Additionally, the shift to homostyly lowers genetic diversity in both the S-genes and their paralogs, as expected in primarily selfing plants. Finally, we tested, for the first time, long-standing theoretical models of changes in S-locus genotypes during early stages of the transition to homostyly, supporting the assumption that two copies of the S-locus might reduce homostyle fitness.

Floral morph variation mediated by clonal growth and pollinator functional groups of Limonium otolepis in a heterostylous fragmented population

Abstract. Heterostyly, a genetic style polymorphism, is linked to symmetric pollen transfer, vital for its maintenance. Clonal growth typically impacts sexual reproduction by influencing pollen transfer. However, the floral morph variation remains poorly understood under the combined effects of pollinators and clonal growth in heterostyly characterized by negative frequency-dependent selection and disassortative mating. We estimated morph ratios, ramets per genet and heterostylous syndrome and quantified legitimate pollen transfer via clonal growth, pollinators and reciprocal herkogamy between floral morphs in Limonium otolepis, a fragmented population composed of five subpopulations in the desert environment of northwestern China, with small flower and large floral morph variation. All subpopulations but one exhibited pollen-stigma morphology dimorphism. The compatibility between mating types with different pollen-stigma morphologies remained consistent regardless of reciprocal herkogamy. Biased ratios and ramets per genet of the two mating types with distinct pollen-stigma morphologies caused asymmetric pollen flow and varying fruit sets in all subpopulations. Short-tongued insects were the primary pollinators due to small flower sizes. However, pollen-feeding Syrphidae sp. triggered asymmetry in pollen flow between high and low sex organs, with short-styled morphs having lower stigma pollen depositions and greater variation. Clonal growth amplified this variation by reducing intermorph pollen transfer. All in all, pollinators and clonal growth jointly drive floral morph variation. H-morphs with the same stigma-anther position and self-incompatibility, which mitigate the disadvantages of sunken low sex organs with differing from the classical homostyly, might arise from long- and short-styled morphs through a 'relaxed selection'. This study is the first to uncover the occurrence of the H-morph and its associated influencing factors in a distylous plant featuring clonal growth, small flowers and a fragmented population.

Two floral forms in the same species-distyly

MAIN CONCLUSION

This paper reviews the progress of research on the morphology, physiology and molecular biology of distyly in plants. It will help to elucidate the mysteries of distyly in plants. Distyly is a unique representative type of heterostyly in plants, primarily characterized by the presence of long style and short style within the flowers of the same species. This interesting trait has always fascinated researchers. With the rapid development of molecular biology, the molecular mechanism for the production of dimorphic styles in plants is also gaining ground. Researchers have been studying plant dimorphic styles from various perspectives. The researchers are gradually unravelling the mechanisms by which plants produce distyly traits. This paper reviews advances in the study of plant dimorphic style characteristics, mainly in terms of the morphology, physiology and molecular biology of plants with dimorphic styles. The aim is to provide a theoretical basis for the study of the mechanism of distyly formation in plants.

Integrative taxonomy in a rapid speciation group associated with mating system transition: A case study in the Primula cicutariifolia complex

Accurate species delimitation is the key to biodiversity conservation and is fundamental to most branches of biology. However, species delimitation remains challenging in those evolutionary radiations associated with mating system transition from outcrossing to self-fertilization, which have frequently occurred in angiosperms and are usually accompanied by rapid speciation. Here, using the Primula cicutariifolia complex as a case, we integrated molecular, morphological and reproductive isolation evidence to test and verify whether its outcrossing (distylous) and selfing (homostylous) populations have developed into independent evolutionary lineages. Phylogenetic trees based on whole plastomes and SNPs of the nuclear genome both indicated that the distylous and homostylous populations grouped into two different clades. Multispecies coalescent, gene flow and genetic structure analyses all supported such two clades as two different genetic entities. In morphology, as expected changes in selfing syndrome, homostylous populations have significantly fewer umbel layers and smaller flower and leaf sizes compared to distylous populations, and the variation range of some floral traits, such as corolla diameter and umbel layers, show obvious discontinuity. Furthermore, hand-pollinated hybridization between the two clades produced almost no seeds, indicating that well post-pollination reproductive isolation has been established between them. Therefore, the distylous and homostylous populations in this studied complex are two independent evolutionary lineages, and thus these distylous populations should be treated as a distinct species, here named Primula qiandaoensis W. Zhang & J.W. Shao sp. nov.. Our empirical study of the P. cicutariifolia complex highlights the importance of applying multiple lines of evidence, in particular genomic data, to delimit species in pervasive evolutionary plant radiations associated with mating system transition.

Inheritance of distyly and homostyly in self-incompatible Primula forbesii

The evolutionary transition from self-incompatible distyly to self-compatible homostyly frequently occurs in heterostylous taxa. Although the inheritance of distyly and homostyly has been deeply studied, our understanding on modifications of the classical simple Mendelian model is still lacking. Primula forbesii, a biennial herb native to southwest China, is a typical distylous species, but after about 20 years of cultivation with open pollination, self-compatible homostyly appeared, providing ideal material for the study of the inheritance of distyly and homostyly. In this study, exogenous homobrassinolide was used to break the heteromorphic incompatibility of P. forbesii. Furthermore, we performed artificial pollination and open-pollination experiments to observe the distribution of floral morphs in progeny produced by different crosses. The viability of seeds from self-pollination was always the lowest among all crosses, and the homozygous S-morph plants (S/S) occurred in artificial pollination experiments but may experience viability selection. The distyly of P. forbesii is governed by a single S-locus, with S-morph dominant hemizygotes (S/-) and L-morph recessive homozygotes (-/-). Homostylous plants have a genotype similar to L-morph plants, and homostyly may be caused by one or more unlinked modifier genes outside the S-locus. Open pollinations confirm that autonomous self-pollination occurs frequently in L-morphs and homostylous plants. This study deepens the understanding of the inheritance of distyly and details a case of homostyly that likely originated from one or more modifier genes.

Diverse mating consequences of the evolutionary breakdown of the sexual polymorphism heterostyly

Reproductive systems of flowering plants are evolutionarily fluid, with mating patterns changing in response to shifts in abiotic conditions, pollination systems, and population characteristics. Changes in mating should be particularly evident in species with sexual polymorphisms that become ecologically destabilized, promoting transitions to alternative reproductive systems. Here, we decompose female mating portfolios (incidence of selfing, outcross mate number, and intermorph mating) in eight populations of Primula oreodoxa, a self-compatible insect-pollinated herb. This species is ancestrally distylous, with populations subdivided into two floral morphs that usually mate with each other (disassortative mating). Stages in the breakdown of polymorphism also occur, including "mixed" populations of distylous and homostylous (self-pollinating) morphs and purely homostylous populations. Population morph ratios vary with elevation in association with differences in pollinator availability, providing an unusual opportunity to investigate changes in mating patterns accompanying transitions in reproductive systems. Unexpectedly, individuals mostly outcrossed randomly, with substantial disassortative mating in at most two distylous populations. As predicted, mixed populations had higher selfing rates than distylous populations, within mixed populations, homostyles selfed almost twice as much as the distylous morphs, and homostylous populations exhibited the highest selfing rates. Populations with homostyles outcrossed with fewer mates and mate number varied negatively with population selfing rates. These differences indicate maintenance of distyly at low elevation, transition to monomorphic selfing at high elevation, and uncertain, possibly variable fates at intermediate elevation. By quantifying the earliest changes in mating that initiate reproductive transitions, our study highlights the key role of mating in promoting evolutionary divergence.

Genomic evidence supports the genetic convergence of a supergene controlling the distylous floral syndrome

Heterostyly, a plant sexual polymorphism controlled by the S-locus supergene, has evolved numerous times among angiosperm lineages and represents a classic example of convergent evolution in form and function. Determining whether underlying molecular convergence occurs could provide insights on constraints to floral evolution. Here, we investigated S-locus genes in distylous Gelsemium (Gelsemiaceae) to determine whether there is evidence of molecular convergence with unrelated distylous species. We used several approaches, including anatomical measurements of sex-organ development and transcriptome and whole-genome sequencing, to identify components of the S-locus supergene. We also performed evolutionary analysis with candidate S-locus genes and compared them with those reported in Primula and Turnera. The candidate S-locus supergene of Gelsemium contained four genes, of which three appear to have originated from gene duplication events within Gelsemiaceae. The style-length genes GeCYP in Gelsemium and CYP734A50 in Primula likely arose from duplication of the same gene, CYP734A1. Three out of four S-locus genes in Gelsemium elegans were hemizygous, as previously reported in Primula and Turnera. We provide genomic evidence on the genetic convergence of the supergene underlying distyly among distantly related angiosperm lineages and help to illuminate the genetic architecture involved in the evolution of heterostyly.

Different molecular changes underlie the same phenotypic transition: Origins and consequences of independent shifts to homostyly within species

The repeated transition from outcrossing to selfing is a key topic in evolutionary biology. However, the molecular basis of such shifts has been rarely examined due to lack of knowledge of the genes controlling these transitions. A classic example of mating system transition is the repeated shift from heterostyly to homostyly. Occurring in 28 angiosperm families, heterostyly is characterized by the reciprocal position of male and female sexual organs in two (or three) distinct, usually self-incompatible floral morphs. Conversely, homostyly is characterized by a single, self-compatible floral morph with reduced separation of male and female organs, facilitating selfing. Here, we investigate the origins of homostyly in Primula vulgaris and its microevolutionary consequences by integrating surveys of the frequency of homostyles in natural populations, DNA sequence analyses of the gene controlling the position of female sexual organs (CYPᵀ), and microsatellite genotyping of both progeny arrays and natural populations characterized by varying frequencies of homostyles. As expected, we found that homostyles displace short-styled individuals, but long-style morphs are maintained at low frequencies within populations. We also demonstrated that homostyles repeatedly evolved from short-styled individuals in association with different types of loss-of-function mutations in CYPᵀ. Additionally, homostyly triggers a shift to selfing, promoting increased inbreeding within and genetic differentiation among populations. Our results elucidate the causes and consequences of repeated transitions to homostyly within species, and the putative mechanisms precluding its fixation in P. vulgaris. This study represents a benchmark for future analyses of losses of heterostyly in other angiosperms.

The rapid appearance of homostyly in a cultivated distylous population of Primula forbesii

Evolutionary breakdown from rigorous outbreeding to self-fertilization frequently occurs in angiosperms. Since the pollinators are not necessary, self-compatible populations often reduce investment in floral display characteristics and pollination reward. Primula forbesii is a biennial herb with distribution restricted to southwest China; it was initially a self-incompatible distylous species, but after 20 years of artificial domestication, homostyly appeared. This change in style provides an ideal material to explore the time required for plant mating systems to adapt to new environmental changes and test whether flower attraction has reduced following transitions to selfing. We did a survey in wild populations of P. forbesii where its seeds were originally collected 20 years ago and recorded the floral morph frequencies and morphologies. The floral morphologies, self-incompatibility, floral scent, and pollinator visitation between distyly and homostyly were compared in greenhouse. Floral morph frequencies of wild populations did not change, while the cultivated population was inclined to L-morph and produced homostyly. Evidence from stigma papillae and pollen size supports the hypothesis that the homostyly possibly originated from mutations of large effect genes in distylous linkage region. Transitions to self-compatible homostyly are accompanied by smaller corolla size, lower amounts of terpenoids, especially linalool and higher amounts of fatty acid derivatives. The main pollinators in the greenhouse were short-tongued Apis cerana. However, homostyly had reduced visiting frequency. The mating system of P. forbesii changed rapidly in just about 20 years of domestication, and our findings confirm the hypothesis that the transition to selfing is accompanied by decreased flower attraction.

The S-Gene YUC6 Pleiotropically Determines Male Mating Type and Pollen Size in Heterostylous Turnera (Passifloraceae): A Novel Neofunctionalization of the YUCCA Gene Family

In heterostylous, self-incompatible Turnera species, a member of the YUCCA gene family, YUC6, resides at the S-locus and has been hypothesized to determine the male mating type. YUCCA gene family members synthesize the auxin, indole-3-acetic acid, via a two-step process involving the TAA gene family. Consequently, it has been speculated that differences in auxin concentration in developing anthers are the biochemical basis underlying the male mating type. Here, we provide empirical evidence that supports this hypothesis. Using a transgenic knockdown approach, we show that YUC6 acts pleiotropically to control both the male physiological mating type and pollen size, but not the filament length dimorphism associated with heterostyly in Turnera. Using qPCR to assess YUC6 expression in different transgenic lines, we demonstrate that the level of YUC6 knockdown correlates with the degree of change observed in the male mating type. Further assessment of YUC6 expression through anther development, in the knockdown lines, suggests that the male mating type is irreversibly determined during a specific developmental window prior to microsporogenesis, which is consistent with the genetically sporophytic nature of this self-incompatibility system. These results represent the first gene controlling male mating type to be characterized in any species with heterostyly.

Direct evidence supporting Darwin's hypothesis of cross-pollination promoted by sex organ reciprocity

The floral phenotype plays a main role in the attraction and fit of pollinators. Both perianth traits and the positioning of sex organs can be subjected to natural selection and determine nonrandom mating patterns in populations. In stylar-polymorphic species, the Darwinian hypothesis predicts increased mating success between individuals with sex organs at equivalent heights (i.e. with higher reciprocity). We used paternity analyses in experimental populations of a stylar-dimorphic species. By comparing the observed mating patterns with those expected under random mating, we tested the effects of sex organ reciprocity and perianth traits on mating success. We also analysed phenotypic selection on perianth traits through female and male functions. The (dis)similarity of parental perianth traits had no direct effects on the mating patterns. Sex organ reciprocity had a positive effect on mating success. Narrow floral tubes increased this effect in upper sex organs. Perianth traits showed little signs of phenotypic selection. Female and absolute fitness measures resulted in different patterns of phenotypic selection. We provide precise empirical evidence of the Darwinian hypothesis about the functioning of stylar polymorphisms, demonstrating that mating patterns are determined by sex organ reciprocity and only those perianth traits which are critical to pollinator fit.

The Complete Chloroplast Genomes of Primula obconica Provide Insight That Neither Species nor Natural Section Represent Monophyletic Taxa in Primula (Primulaceae)

The genus Primula (Primulaceae) comprises more than 500 species, with 300 species distributed in China. The contradictory results between systematic analyses and morphology-based taxonomy make taxonomy studies difficult. Furthermore, frequent introgression between closely related species of Primula can result in non-monophyletic species. In this study, the complete chloroplast genome of sixteen Primula obconica subsp. obconica individuals were assembled and compared with 84 accessions of 74 species from 21 sections of the 24 sections of the genus in China. The plastome sizes of P. obconica subsp. obconica range from 153,584 bp to 154,028 bp. Genome-wide variations were detected, and 1915 high-quality SNPs and 346 InDels were found. Most SNPs were detected in downstream and upstream gene regions (45.549% and 41.91%). Two cultivated accessions, ZP1 and ZP2, were abundant with SSRs. Moreover, 12 SSRs shared by 9 accessions showed variations that may be used as molecular markers for population genetic studies. The phylogenetic tree showed that P. obconica subsp. obconica cluster into two independent clades. Two subspecies have highly recognizable morphological characteristics, isolated geographical distribution areas, and distinct phylogenetic relationships compared with P. obconica subsp. obconica. We elevate the two subspecies of P. obconica to separate species. Our phylogenetic tree is largely inconsistent with morphology-based taxonomy. Twenty-one sections of Primula were mainly divided into three clades. The monophyly of Sect. Auganthus, Sect. Minutissimae, Sect. Sikkimensis, Sect. Petiolares, and Sect. Ranunculoides are well supported in the phylogenetic tree. The Sect. Obconicolisteri, Sect. Monocarpicae, Sect. Carolinella, Sect. Cortusoides, Sect. Aleuritia, Sect. Denticulata, Sect. Proliferae Pax, and Sect. Crystallophlomis are not a monophyletic group. The possible explanations for non-monophyly may be hybridization, polyploidization, recent introgression, incorrect taxonomy, or chloroplast capture. Multiple genomic data and population genetic studies are therefore needed to reveal the evolutionary history of Primula. Our results provided valuable information for intraspecific variation and phylogenetic relationships within Primula.

Evolution of Autonomous Selfing in Marginal Habitats: Spatiotemporal Variation in the Floral Traits of the Distylous Primula wannanensis

Outcrossing plant species are more likely to exhibit autonomous selfing in marginal habitats to ensure reproduction under conditions of limited pollinator and/or mate availability. Distyly is a classical paradigm that promotes outcrossing; however, little is known about the variation in floral traits associated with distylous syndrome in marginal populations. In this study, we compared the variation in floral traits including stigma and anther height, corolla tube length, herkogamy, and corolla diameter between the central and peripheral populations of the distylous Primula wannanensis, and assessed the variation of floral traits at early and late florescence stages for each population. To evaluate the potential consequences of the variation in floral traits on the mating system, we investigated seed set in each population under both open-pollinated and pollinator-excluded conditions. The flower size of both short- and long-styled morphs was significantly reduced in late-opening flowers compared with early opening flowers in both central and peripheral populations. Sex-organ reciprocity was perfect in early opening flowers; however, it was largely weakened in the late-opening flowers of peripheral populations compared with central populations. Of these flowers, disproportionate change in stigma height (elongated in S-morph and shortened in L-morph) was the main cause of reduced herkogamy, and seed set was fairly high under pollinator-excluded condition. Our results provide empirical support for the hypothesis on the evolution of delayed autonomous selfing in marginal populations of distylous species. Unsatisfactory pollinator service is likely to have promoted reproductive assurance of distylous plants with largely reduced herkogamy mimicking "homostyles."

Development of the first microsatellite markers using high-throughput sequencing for a hexaploid coastal species, Tournefortia argentea L. f. (Boraginaceae)

BACKGROUD

Tournefortia argentea L. f. is a hexaploid shrub or tree species with ecological and evolutionary significances, which forms the fringe of vegetation closest to the sea on tropical coral islands. Previous studies have never addressed on genetic information, and thus genomic resources remain scarce.

METHODS AND RESULTS

We used nine individuals from different islands to identify polymorphic microsatellites of T. argentea by Illumina high-throughput sequencing. Thirty-five polymorphic microsatellite markers were developed. Characteristics of each locus were tested using 48 individuals collected from three populations of T. argentea. A total of 320 alleles were found across the 35 microsatellite loci. The number of alleles per locus ranged from 5 to 15, with an average of 9.1. Observed and expected heterozygosities in each locus per population varied from 0.000 to 1.000 and from 0.000 to 0.893, respectively.

CONCLUSIONS

In this study, we report the development of 35 polymorphic microsatellite markers based on Illumina high-throughput sequencing. These markers will facilitate the investigations of genetic diversity, population structures and evolutionary history of T. argentea.

Ecological Interactions, Environmental Gradients, and Gene Flow in Local Adaptation

Despite long-standing interest in local adaptation of plants to their biotic and abiotic environment, existing theory, and many case studies, little work to date has addressed within-species evolution of concerted strategies and how these might contrast with patterns across species. Here we consider the interactions between pollinators, herbivores, and resource availability in shaping plant local adaptation, how these interactions impact plant phenotypes and gene flow, and the conditions where multiple traits align along major environmental gradients such as latitude and elevation. Continued work in emerging model systems will benefit from the melding of classic experimental approaches with novel population genetic analyses to reveal patterns and processes in plant local adaptation.

Pollen transport: Illuminating a key mechanism of disassortative pollination

Floral sexual polymorphisms have evolved repeatedly in angiosperms and are thought to reduce self-pollination and increase pollen export. Using a powerful pollen-labeling technique, quantum dots, a new study shows that pollen placement on pollinator bodies plays a critical role in disassortative pollination.

The Genomic Selfing Syndrome Accompanies the Evolutionary Breakdown of Heterostyly

The evolutionary transition from outcrossing to selfing can have important genomic consequences. Decreased effective population size and the reduced efficacy of selection are predicted to play an important role in the molecular evolution of the genomes of selfing species. We investigated evidence for molecular signatures of the genomic selfing syndrome using 66 species of Primula including distylous (outcrossing) and derived homostylous (selfing) taxa. We complemented our comparative analysis with a microevolutionary study of P. chungensis, which is polymorphic for mating system and consists of both distylous and homostylous populations. We generated chloroplast and nuclear genomic data sets for distylous, homostylous, and distylous–homostylous species and identified patterns of nonsynonymous to synonymous divergence (d_N/d_S) and polymorphism (π_N/π_S) in species or lineages with contrasting mating systems. Our analysis of coding sequence divergence and polymorphism detected strongly reduced genetic diversity and heterozygosity, decreased efficacy of purifying selection, purging of large-effect deleterious mutations, and lower rates of adaptive evolution in samples from homostylous compared with distylous populations, consistent with theoretical expectations of the genomic selfing syndrome. Our results demonstrate that self-fertilization is a major driver of molecular evolutionary processes with genomic signatures of selfing evident in both old and relatively young homostylous populations.

Sexual conflict in protandrous flowers and the evolution of gynodioecy

Sexual interference between male and female function in hermaphrodite plants is reduced by protandry. In environments with insufficient pollinator service, prolongation of male function owing to limited pollen removal could restrict the duration of female function and lower seed production. We provide evidence that this form of sexual conflict has played a role in the spread of females in gynodioecious populations of Cyananthus delavayi in the pollen-limited environments in which this subalpine species occurs. Using field experiments involving artificial pollen removal from the strongly protandrous flowers of hermaphrodites, we demonstrated a trade-off between male- and female-phase duration with no influence on overall floral longevity. Pollen removal at the beginning of anthesis resulted in hermaphrodite seed production matching that of females. In contrast, restricted pollen removal increased the duration of male function at the expense of female function lowering maternal fertility compared to females. This pattern was evident in five populations with females experiencing a twofold average seed fertility advantage compared to hermaphrodites. Gynodioecy often appears to evolve from protandrous ancestors and pollen limitation is widespread in flowering plants suggesting that sexual conflict may play an unappreciated role in the evolution of this form of sexual dimorphism.

Dynamics behind disjunct distribution, hotspot-edge refugia, and discordant RADseq/mtDNA variability: insights from the Emei mustache toad

Background

The distribution of genetic diversity and the underlying processes are important for conservation planning but are unknown for most species and have not been well studied in many regions. In East Asia, the Sichuan Basin and surrounding mountains constitute an understudied region that exhibits a “ring” of high species richness overlapping the eastern edge of the global biodiversity hotspot Mountains of Southwest China. We examine the distributional history and genetic diversification of the Emei mustache toad Leptobrachium boringii, a typical “ring” element characterized by disjunct ranges in the mountains, by integrating time-calibrated gene tree, genetic variability, individual-level clustering, inference of population splitting and mixing from allele frequencies, and paleoclimatic suitability modeling.

Results

The results reveal extensive range dynamics, including secondary contact after long-term isolation via westward dispersal accompanied by variability loss. They allow the proposal of a model that combines recurrent contractions caused by Quaternary climatic changes and some failed expansions under suitable conditions for explaining the shared disjunct distribution pattern. Providing exceptional low-elevation habitats in the hotspot area, the eastern edge harbors both long-term refugial and young immigrant populations. This finding and a synthesis of evidence from other taxa demonstrate that a certain contributor to biodiversity, one that preserves and receives low-elevation elements of the east in this case, can be significant for only a particular part of a hotspot. By clarifying the low variability of these refugial populations, we show that discordant mitochondrial estimates of diversity can be obtained for populations that experienced admixture, which would have unlikely left proportional immigrant alleles for each locus.

Conclusions

Dispersal after long-term isolation can explain much of the spatial distribution of genetic diversity in this species, while secondary contact and long-term persistence do not guarantee a large variation. The model for the formation of disjunct ranges may apply to many other taxa isolated in the mountains surrounding the Sichuan Basin. Furthermore, this study provides insights into the heterogeneous nature of hotspots and discordant variability obtained from genome-wide and mitochondrial data.

Reproductive transitions in plants and animals: selfing syndrome, sexual selection and speciation

The evolution of predominant self-fertilisation frequently coincides with the evolution of a collection of phenotypes that comprise the 'selfing syndrome', in both plants and animals. Genomic features also display a selfing syndrome. Selfing syndrome traits often involve changes to male and female reproductive characters that were subject to sexual selection and sexual conflict in the obligatorily outcrossing ancestor, including the gametic phase for both plants and animals. Rapid evolution of reproductive traits, due to both relaxed selection and directional selection under the new status of predominant selfing, lays the genetic groundwork for reproductive isolation. Consequently, shifts in sexual selection pressures coupled to transitions to selfing provide a powerful paradigm for investigating the speciation process. Plant and animal studies, however, emphasise distinct selective forces influencing reproductive-mode transitions: genetic transmission advantage to selfing or reproductive assurance outweighing the costs of inbreeding depression vs the costs of males and meiosis. Here, I synthesise links between sexual selection, evolution of selfing and speciation, with particular focus on identifying commonalities and differences between plant and animal systems and pointing to areas warranting further synergy.

'A most complex marriage arrangement': recent advances on heterostyly and unresolved questions

Heterostylous genetic polymorphisms provide paradigmatic systems for investigating adaptation and natural selection. Populations are usually comprised of two (distyly) or three (tristyly) mating types, maintained by negative frequency-dependent selection resulting from disassortative mating. Theory predicts this mating system should result in equal style-morph ratios (isoplethy) at equilibrium. Here, I review recent advances on heterostyly, focusing on examples challenging stereotypical depictions of the polymorphism and unresolved questions. Comparative analyses indicate multiple origins of heterostyly, often within lineages. Ecological studies demonstrate that structural components of heterostyly are adaptations improving the proficiency of animal-mediated cross-pollination and reducing pollen wastage. Both neutral and selective processes cause deviations from isoplethy in heterostylous populations, and, under some ecological and demographic conditions, cause breakdown of the polymorphism, resulting in either the evolution of autogamy and mixed mating, or transitions to alternative outcrossing systems, including dioecy. Earlier ideas on the genetic architecture of the S-locus supergene governing distyly have recently been overturned by discovery that the dominant S-haplotype is a hemizygous region absent from the s-haplotype. Ecological, phylogenetic and molecular genetic data have validated some features of theoretical models on the selection of the polymorphism. Although heterostyly is the best-understood floral polymorphism in angiosperms, many unanswered questions remain.

Phylogenomic analysis reveals multiple evolutionary origins of selfing from outcrossing in a lineage of heterostylous plants

Evolutionary transitions from outcrossing to selfing often occur in heterostylous plants. Selfing homostyles originate within distylous populations and frequently evolve to become reproductively isolated species. We investigated this process in 10 species of Primula section Obconicolisteri using phylogenomic approaches and inferred how often homostyly originated from distyly and its consequences for population genetic diversity and floral trait evolution. We estimated phylogenetic relationships and reconstructed character evolution using the whole plastome comprised of 76 protein-coding genes. To investigate mating patterns and genetic diversity we screened 15 microsatellite loci in 40 populations. We compared floral traits among distylous and homostylous populations to determine how phenotypically differentiated homostyles were from their distylous ancestors. Section Obconicolisteri was monophyletic and we estimated multiple independent transitions from distyly to homostyly. High selfing rates characterised homostylous populations and this was associated with reduced genetic diversity. Flower size and pollen production were reduced in homostylous populations, but pollen size was significantly larger in some homostyles than in distylous morphs. Repeated transitions to selfing in section Obconicolisteri are likely to have been fostered by the complex montane environments that species occupy. Unsatisfactory pollinator service is likely to have promoted reproductive assurance in homostyles leading to subsequent population divergence through isolation.

Global transcriptome and gene co-expression network analyses on the development of distyly in Primula oreodoxa

Distyly is a genetically controlled flower polymorphism that has intrigued both botanists and evolutionary biologists ever since Darwin's time. Despite extensive reports on the pollination and evolution of distylous systems, the genetic basis and mechanism of molecular regulation remain unclear. In the present study, comparative transcriptome profiling was conducted in primrose (Primula oreodoxa), the prime research model for heterostyly. Thirty-six transcriptomes were sequenced for styles at different stages and corolla tube in the three morphs of P. oreodoxa. Large numbers of differentially expressed genes (DEGs) were detected in the transcriptomes of styles across different morphs. Several transcription factors (TFs) and phytohormone metabolism-related genes were highlighted in S-morphs. A growing number of genes showed differential expression patterns along with the development of styles, suggesting that the genetic control of distyly may be more complicated than ever expected. Analysis of co-expression networks and module-trait relationships identified modules significantly associated with style development. CYP734A50, a key S-locus gene whose products degrade brassinosteroids, was co-expressed with many genes in the module and showed significant negative association with style length. In addition, crucial TFs involved in phytohormone signaling pathways were found to be connected with CYP734A50 in the co-expression module. Our global transcriptomic analysis has identified DEGs that are potentially involved in regulation of style length in P. oreodoxa, and may shed light on the evolution and broad biological processes of heterostyly.

Intraspecific variation of self-incompatibility in the distylous plant Primula merrilliana

Heteromorphic self-incompatibility can prevent self- and intramorph fertilization while favouring intermorph mating and the maintenance of morph-ratio stability in heterostylous populations. However, variation in the expression of self-incompatibility intraspecies has seldom been assessed. Through hand pollinations and microsatellite markers, the variation in the expression of self-incompatibility and genetic diversity were studied in distylous plant Primula merrilliana. We discovered that the strength of self-incompatibility varied extensively among individuals and populations, from pronounced to weaker self-incompatibility in distylous populations, all the way to strong self-compatibility in homostylous populations. Each distylous population included self-incompatible (SI), partly self-compatible (PSC) and self-compatible (SC) individuals, with the index of self-compatibility (ISC) ranging from 0.07 to 0.68 across populations. Self-compatible populations (ISC > 0.25) were not genetically clustered but were more closely related to populations with high SI and SC individuals were mixed with SI individuals within populations. The ISC and the proportions of SC and PSC individuals were higher in peripheral than in central populations, but no decrease of genetic diversity and no deviations of floral morph ratio from isoplethy were detected in peripheral populations. Additionally, the expression of self-incompatibility was stronger in long-styled flowers than in short-styled flowers. The variation in the strength of self-incompatibility documented in P. merrilliana cautions against the estimation of ISC from a few individuals or populations in distylous species and provides a more complex and nuanced understanding of the role of self-incompatibility in heterostyly.

Characterization of 30 microsatellite markers in distylous Primula sinolisteri (Primulaceae) using HiSeq sequencing

PREMISE OF THE STUDY

Microsatellite markers were developed for Primula sinolisteri, a perennial distylous herb belonging to section Obconicolisteri (Primulaceae), to facilitate future investigations of the population genetics and mating patterns of populations in this species.

METHODS AND RESULTS

We developed 30 microsatellite markers for P. sinolisteri using HiSeq X-Ten sequencing and measured polymorphism and genetic diversity in a sample of 36 individuals from three natural populations. The markers displayed relatively high polymorphism, with the number of observed alleles per locus ranging from one to 19 (mean = 4.42). The observed and expected heterozygosity ranged from 0-1.000 and 0.083-0.882, respectively. Twenty-nine of the loci were also successfully amplified in homostylous P. sinolisteri var. aspera.

CONCLUSIONS

The microsatellite markers we have identified in P. sinolisteri provide powerful tools for investigating patterns of population genetic diversity and the evolutionary relationships between heterostyly and homostyly in this species.

Pollinator-mediated selection on floral traits varies in space and between morphs in Primula secundiflora

Elucidating how variation in selection shapes the evolution of flowers is key to understanding adaptive differentiation processes. We estimated pollinator-mediated selection through female function in L-morph (long-style and short-anther phenotype) and S-morph (short-style and long-anther phenotype) flowers among four Primula secundiflora populations with different pollinator assemblages. Variation in pollinator assemblage strongly contributed to differences in reproductive success among populations and between morphs of the primrose species. A wider corolla tube width was selected in the bumblebee-dominated populations, whereas shorter corolla tube length and wider corolla tube width were selected in the syrphid fly-dominated populations. Morph-specific variation in pollinator-mediated selection on corolla tube length was detected in the syrphid fly-dominated populations. A shorter corolla tube was selected in the L-morph flowers. However, similar selective pressure on this trait was not observed in the S-morph flowers. These results show that variation in pollinator assemblage leads to variation in selection in space and between morphs. The findings highlight the potential forces of different pollinator agents in driving floral evolution in this primrose species.

Genetics of distyly and homostyly in a self-compatible Primula

The transition from outcrossing to selfing through the breakdown of distyly to homostyly has occurred repeatedly among families of flowering plants. Homostyles can originate by major gene changes at the S-locus linkage group, or by unlinked polygenic modifiers. Here, we investigate the inheritance of distyly and homostyly in Primula oreodoxa, a subalpine herb endemic to Sichuan, China. Controlled self- and cross-pollinations confirmed that P. oreodoxa unlike most heterostylous species is fully self-compatible. Segregation patterns indicated that the inheritance of distyly is governed by a single Mendelian locus with the short-styled morph carrying at least one dominant S-allele (S-) and long-styled plants homozygous recessive (ss). Crossing data were consistent with a model in which homostyly results from genetic changes at the distylous linkage group, with the homostylous allele (S^h) dominant to the long-styled allele (s), but recessive to the short-styled allele (S). Progeny tests of open-pollinated seed families revealed high rates of intermorph mating in the L-morph but considerable selfing and possibly intramorph mating in the S-morph and in homostyles. S-morph plants homozygous at the S-locus (SS) occurred in several populations but may experience viability selection. The crossing data from distylous and homostylous plants are consistent with either recombination at the S-locus governing distyly, or mutation at gene(s) controlling sex-organ height; both models predict the same patterns of segregation. Recent studies on the molecular genetics of distyly in Primula demonstrating the hemizygous nature of genes at the S-locus make it more likely that homostyles have resulted from mutation rather than recombination.

Development of microsatellite markers for Primula oreodoxa (Primulaceae), a distylous-homostylous species

PREMISE OF THE STUDY

Microsatellite markers were developed for a distylous-homostylous species, Primula oreodoxa (Primulaceae), to investigate the mating patterns and gene flow in the species.

METHODS AND RESULTS

Using RAD sequencing, 42,777 contigs were generated. A total of 1566 putative simple sequence repeat (SSR) loci were identified by MISA, and 1433 primer sets were designed. After initial screening of 107 SSR loci, 24 loci displayed polymorphism. The number of alleles per locus detected ranged from one to 11 in the sampled populations of P. oreodoxa. Levels of observed and expected heterozygosity for each locus ranged from 0 to 0.917 and 0 to 0.816, respectively. Fifteen and 13 of these loci could be successfully amplified in the two congeneric species P. obconica and P. heucherifolia, respectively.

CONCLUSIONS

The SSR markers developed here will be valuable for genetic analysis and elucidation of mating patterns in P. oreodoxa.

Morphology, adaptation and speciation

The study of plant evolution and development in a phylogenetic context has accelerated research advances in both areas over the last decade. The addition of a robust phylogeny for plant taxa based on DNA as well as morphology has given a strong context for this research. Genetics and genomics, including sequencing of many genes, and a better understanding of non-genetic, responsive changes, by plants have increased knowledge of how the different body forms of plants have arisen. Here, I overview the papers in this Special Issue of Annals of Botany on Morphological Adaptation, bringing together a range of papers that link phylogeny and morphology. These lead to models of development and functional adaptation across a range of plant systems, with implications for ecology and ecosystems, as well as development and evolution.