Presence versus absence of CYP734A50 underlies the style-length dimorphism in primroses

Converging on long and short: The genetics, molecular biology and evolution of heterostyly

Heterostyly is a fascinating floral polymorphism that enhances outcrossing. In heterostylous species the flowers of the two or three morphs differ in multiple traits, including reciprocal reproductive-organ placement and self-incompatibility. These traits are controlled by individual genes within an S-locus supergene, whose suppressed recombination ensures the coordinated inheritance of the morph phenotypes. Recent breakthroughs about the genetic and molecular basis of heterostyly have resulted from studies on many independently evolved instances and include the following: The S-locus is a hemizygous region comprising several individual genes in multiple heterostylous taxa. In many systems, a single gene within the S-locus plays dual roles in regulating both female traits of style length and self-incompatibility type, often involving brassinosteroid signalling. The S-loci have evolved through stepwise or segmental duplication in different lineages. The frequent breakdown of heterostyly generally results from individual mutations at the S-locus and leads to a genomic selfing syndrome. These discoveries suggest convergent and genetically constrained evolution of heterostyly at the molecular level.

Spiral phyllotaxis predicts left-right asymmetric growth and style deflection in mirror-image flowers of Cyanella alba

Many animals and plants show left-right (LR) asymmetry. The LR asymmetry of mirror-image flowers has clear functional significance, with the reciprocal placement of male and female organs in left- versus right-handed flowers promoting cross-pollination. Here, we study how handedness of mirror-image flowers is determined and elaborated during development in the South African geophyte Cyanella alba. Inflorescences of C. alba produce flowers with a largely consistent handedness. However, this handedness has no simple genetic basis and individual plants can switch their predominant handedness between years. Rather, it is the direction of the phyllotactic spiral that predicts floral handedness. Style deflection is driven by increased cell expansion in the adaxial carpel facing the next oldest flower compared to the other adaxial carpel. The more expanding carpel shows transcriptional signatures of increased auxin signaling and auxin application can reverse the orientation of style deflection. We propose that a recently described inherent LR auxin asymmetry in the initiating organs of spiral phyllotaxis determines handedness in C. alba, creating a stable yet non-genetic floral polymorphism. This mechanism links chirality across different levels of plant development and exploits a developmental constraint in a core patterning process to produce morphological variation of ecological relevance.

Rapid detection of RNase-based self-incompatibility in Lysimachia monelli (Primulaceae)

PREMISE

Primroses famously employ a system that simultaneously expresses distyly and filters out self-pollen. Other species in the Primulaceae family, including Lysimachia monelli (blue pimpernel), also express self-incompatibility (SI), but involving a system with distinct features and an unknown molecular genetic basis.

METHODS

We utilize a candidate-based transcriptome sequencing (RNA-seq) approach, relying on candidate T2/S-RNase Class III and S-linked F-box-motif-containing genes and harnessing the unusual evolutionary and genetic features of SI, to examine whether an RNase-based mechanism underlies SI in L. monelli. We term this approach "SI detection with RNA-seq" (SIDR).

RESULTS

The results of sequencing, crossing, population genetics, and molecular evolutionary features each support a causal association linking the recovered genotypes with SI phenotypes. The finding of RNase-based SI in Primulaceae (Ericales) all but cements the long-held view that this mechanism was present in the ancestral pentapetal eudicot, whose descendants now comprise two-thirds of angiosperms. It also significantly narrows the plausible maximum age for the heterostyly evolution within the family.

CONCLUSIONS

SIDR is powerful, flexible, inexpensive, and most critically enables work in often-neglected species. It may be used with or without candidate genes to close enormous gaps in understanding the genetic basis of SI and the history of breeding system evolution.

Cell elongation and altered phytohormone levels play a role in establishing distyly in Averrhoa carambola

The flowers of distylous plants exhibit two distinct morphologies that facilitate precise pollen transfer. Averrhoa carambola, a woody plant characterized by distyly, has an unclear molecular regulatory mechanism underlying this trait. Its prolonged flowering period and substantial flower production render it an excellent model for investigating the distylous syndrome. This study aims to elucidate the mechanism of distyly in A. carambola and to identify the regulatory genes. The long-style cultivar 'Daguo Tianyangtao 1' and the short-style cultivar 'Daguo Tianyangtao 3' were selected as models for this investigation. We examined phenotypic characteristics, anatomical structures, and endogenous hormone content associated with distyly. Transcriptomic data were utilized to pinpoint candidate genes involved in the regulation of distyly, followed by a bioinformatics analysis these genes. The results indicate that variations in cell elongation contribute to the differential heights of stigmas and anthers in A. carambola, thereby resulting in the distylous syndrome. Auxins, Gibberellin A3 (GA3), Gibberellin A4 (GA4), and brassinolide (BL) were found to influence elongation of styles, whereas Gibberellin A1 (GA1) and GA4 affected filament elongation. Transcriptome sequencing analysis identified 34 hormone-related differentially expressed genes (DEGs) and 16 cell development-related DEGs in different morphs of pistils, and 29 hormone-related DEGs and 22 cell development-related DEGs were identified in different morphs of stamens. Four candidate genes-AcaBRU1, AcaPRE1, AcaXTH2, and AcaEXPA11-were found to possess conserved motifs characteristic of their respective families. Consequently, various plant hormones modulate the expression of response genes, leading to differences in elongation of style and filament cells between different flower types of A. carambola, thereby promoting the distylous syndrome. This study provides a theoretical basis for understanding the mechanisms of distyly formation in woody plants.

Phosphoproteomic analysis of distylous Turnera subulata identifies pathways related to endoreduplication that correlate with reciprocal herkogamy

PREMISE

A multi-omic approach was used to explore proteins and networks hypothetically important for establishing filament dimorphisms in heterostylous Turnera subulata (Sm.) as an exploratory method to identify genes for future empirical research.

METHODS

Mass spectrometry (MS) was used to identify differentially expressed proteins and differentially phosphorylated peptides in the developing filaments between the L- and S-morphs. RNAseq was used to generate a co-expression network of the developing filaments, MS data were mapped to the co-expression network to identify hypothetical relationships between the S-gene responsible for filament dimorphisms and differentially expressed proteins.

RESULTS

Mapping all MS identified proteins to a co-expression network of the S-morph's developing filaments identified several clusters containing SPH1 and other differentially expressed or phosphorylated proteins. Co-expression analysis clustered CDKG2, a protein that induces endoreduplication, and SPH1-suggesting a shared biological function. MS analysis suggests that the protein is present and phosphorylated only in the S-morph, and thus active only in the S-morph. A series of CDKG2 regulators, including ATM1, and cell cycle regulators also correlated with the presence of reciprocal herkogamy, supporting our interest in the protein.

CONCLUSIONS

This work has built a foundation for future empirical work, specifically supporting the role of CDKG2 and ATM1 in promoting filament elongation in response to SPH1 perception.

Sex on Steroids: How Brassinosteroids Shape Reproductive Development in Flowering Plants

Since the discovery of brassinolide in the pollen of rapeseed, brassinosteroids (BRs) have consistently been associated with reproductive traits. However, compared to what is known for how BRs shape vegetative development, the understanding of how these hormones regulate reproductive traits is comparatively still lacking. Nevertheless, there is now considerable evidence that BRs regulate almost all aspects of reproduction, from ovule and pollen formation to seed and fruit development. Here, we review the current body of knowledge on how BRs regulate reproductive processes in plants and what is known about how these pathways are transduced at the molecular level. We also discuss how the manipulation of BR biosynthesis and signaling can be a promising avenue for improving crop traits that rely on efficient reproduction. We thus propose that BRs hold an untapped potential for plant breeding, which could contribute to attaining food security in the coming years.

Unveiling the Genome-Wide Consequences of Range Expansion and Mating System Transitions in Primula vulgaris

Genetic diversity is heterogeneously distributed among populations of the same species, due to the joint effects of multiple demographic processes, including range contractions and expansions, and mating systems shifts. Here, we ask how both processes shape genomic diversity in space and time in the classical Primula vulgaris model. This perennial herb originated in the Caucasus region and was hypothesized to have expanded westward following glacial retreat in the Quaternary. Moreover, this species is a long-standing model for mating system transitions, exemplified by shifts from heterostyly to homostyly. Leveraging a high-quality reference genome of the closely related Primula veris and whole-genome resequencing data from both heterostylous and homostylous individuals from populations encompassing a wide distribution of P. vulgaris, we reconstructed the demographic history of P. vulgaris. Results are compatible with the previously proposed hypothesis of range expansion from the Caucasus region approximately 79,000 years ago and suggest later shifts to homostyly following rather than preceding postglacial colonization of England. Furthermore, in accordance with population genetic theoretical predictions, both processes are associated with reduced genetic diversity, increased linkage disequilibrium, and reduced efficacy of purifying selection. A novel result concerns the contrasting effects of range expansion versus shift to homostyly on transposable elements, for the former, process is associated with changes in transposable element genomic content, while the latter is not. Jointly, our results elucidate how the interactions among range expansion, transitions to selfing, and Quaternary climatic oscillations shape plant evolution.

The molecular basis of phenotypic evolution: beyond the usual suspects

It has been well documented that mutations in coding DNA or cis-regulatory elements underlie natural phenotypic variation in many organisms. However, the development of sophisticated functional tools in recent years in a wide range of traditionally non-model systems have revealed many 'unusual suspects' in the molecular bases of phenotypic evolution, including upstream open reading frames (uORFs), cryptic splice sites, and small RNAs. Furthermore, large-scale genome sequencing, especially long-read sequencing, has identified a cornucopia of structural variation underlying phenotypic divergence and elucidated the composition of supergenes that control complex multi-trait polymorphisms. In this review article we highlight recent studies that demonstrate this great diversity of molecular mechanisms producing adaptive genetic variation and the panoply of evolutionary paths leading to the 'grandeur of life'.

The Primula edelbergii S-locus is an example of a jumping supergene

Research on supergenes, non-recombining genomic regions housing tightly linked genes that control complex phenotypes, has recently gained prominence in genomics. Heterostyly, a floral heteromorphism promoting outcrossing in several angiosperm families, is controlled by the S-locus supergene. The S-locus has been studied primarily in closely related Primula species and, more recently, in other groups that independently evolved heterostyly. However, it remains unknown whether genetic architecture and composition of the S-locus are maintained among species that share a common origin of heterostyly and subsequently diverged across larger time scales. To address this research gap, we present a chromosome-scale genome assembly of Primula edelbergii, a species that shares the same origin of heterostyly with Primula veris (whose S-locus has been characterized) but diverged from it 18 million years ago. Comparative genomic analyses between these two species allowed us to show, for the first time, that the S-locus can 'jump' (i.e. translocate) between chromosomes maintaining its function in controlling heterostyly. Additionally, we found that four S-locus genes were conserved but reshuffled within the supergene, seemingly without affecting their expression, thus we could not detect changes explaining the lack of self-incompatibility in P. edelbergii. Furthermore, we confirmed that the S-locus is not undergoing genetic degeneration. Finally, we investigated P. edelbergii evolutionary history within Ericales in terms of whole genome duplications and transposable element accumulation. In summary, our work provides a valuable resource for comparative analyses aimed at investigating the genetics of heterostyly and the pivotal role of supergenes in shaping the evolution of complex phenotypes.

Integrated physiological analyses, transcriptome, and DNA methylation reveal superiority of pear stigma-style complex development regulation

Styles and stigmas are crucial components of the fertilization process that allows a pear tree to bear fruit. The information regarding the development mechanism of pear style and stigma is still unclear. Our results demonstrated that IAA, ABA, and BR are significantly increased at 1 DBF, while JA is decreased at 5 DBF. The fructose and starch contents significantly increased at 1 DBF when the style with stigma was ready for pollination. Transcriptome and DNA methylation analysis showed 8087 DEGs and 3771 DMRs were enriched in plant hormones biosynthesis, carbohydrate biosynthesis and metabolism, and TFs in 1 DBF as compared with 7 DBF. The CHH methylation type of DMRs accounts for 84.75%. Most DMRs of CHH upregulated in 1 DBF vs. 7 DBF. This study found for the first time that transcription factor ERFs and DNA methylation are involved in regulating the growth and development of fruit plant style and stigma.

A hemizygous supergene controls homomorphic and heteromorphic self-incompatibility systems in Oleaceae

Self-incompatibility (SI) has evolved independently multiple times and prevents self-fertilization in hermaphrodite angiosperms. Several groups of Oleaceae such as jasmines exhibit distylous flowers, with two compatibility groups each associated with a specific floral morph.1 Other Oleaceae species in the olive tribe have two compatibility groups without associated morphological variation.2,3,4,5 The genetic basis of both homomorphic and dimorphic SI systems in Oleaceae is unknown. By comparing genomic sequences of three olive subspecies (Olea europaea) belonging to the two compatibility groups, we first locate the genetic determinants of SI within a 700-kb hemizygous region present only in one compatibility group. We then demonstrate that the homologous hemizygous region also controls distyly in jasmine. Phylogenetic analyses support a common origin of both systems, following a segmental genomic duplication in a common ancestor. Examination of the gene content of the hemizygous region in different jasmine and olive species suggests that the mechanisms determining compatibility groups and floral phenotypes (whether homomorphic or dimorphic) in Oleaceae rely on the presence/absence of two genes involved in gibberellin and brassinosteroid regulation.

The homomorphic self-incompatibility system in Oleaceae is controlled by a hemizygous genomic region expressing a gibberellin pathway gene

In flowering plants, outcrossing is commonly ensured by self-incompatibility (SI) systems. These can be homomorphic (typically with many different allelic specificities) or can accompany flower heteromorphism (mostly with just two specificities and corresponding floral types). The SI system of the Oleaceae family is unusual, with the long-term maintenance of only two specificities but often without flower morphology differences. To elucidate the genomic architecture and molecular basis of this SI system, we obtained chromosome-scale genome assemblies of Phillyrea angustifolia individuals and related them to a genetic map. The S-locus region proved to have a segregating 543-kb indel unique to one specificity, suggesting a hemizygous region, as observed in all distylous systems so far studied at the genomic level. Only one of the predicted genes in this indel region is found in the olive tree, Olea europaea, genome, also within a segregating indel. We describe complete association between the presence/absence of this gene and the SI types determined for individuals of seven distantly related Oleaceae species. This gene is predicted to be involved in catabolism of the gibberellic acid (GA) hormone, and experimental manipulation of GA levels in developing buds modified the male and female SI responses of the two specificities in different ways. Our results provide a unique example of a homomorphic SI system, where a single conserved gibberellin-related gene in a hemizygous indel underlies the long-term maintenance of two groups of reproductive compatibility.

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.

PfPIN5 promotes style elongation by regulating cell length in Primula forbesii Franch

BACKGROUND AND AIMS

Style dimorphism is one of the polymorphic characteristics of flowers in heterostylous plants, which have two types of flowers: the pin morph, with long styles and shorter anthers, and the thrum morph, with short styles and longer anthers. The formation of dimorphic styles has received attention in the plant world. Previous studies showed that CYP734A50 in Primula determined style length and limited style elongation and that the brassinosteroid metabolic pathway was involved in regulation of style length. However, it is unknown whether there are other factors affecting the style length of Primula.

METHODS

Differentially expressed genes highly expressed in pin morph styles were screened based on Primula forbesii transcriptome data. Virus-induced gene silencing was used to silence these genes, and the style length and anatomical changes were observed 20 days after injection.

KEY RESULTS

PfPIN5 was highly expressed in pin morph styles. When PfPIN5 was silenced, the style length was shortened in pin and long-homostyle plants by shortening the length of style cells. Moreover, silencing CYP734A50 in thrum morph plants increased the expression level of PfPIN5 significantly, and the style length increased. The results indicated that PfPIN5, an auxin efflux transporter gene, contributed to regulation of style elongation in P. forbesii.

CONCLUSIONS

The results implied that the auxin pathway might also be involved in the formation of styles of P. forbesii, providing a new pathway for elucidating the molecular mechanism of style elongation in P. forbesii.

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.

The genetic control of herkogamy

Herkogamy is the spatial separation of anthers and stigmas within complete flowers, and is a key floral trait that promotes outcrossing in many angiosperms. The degree of separation between pollen-producing anthers and receptive stigmas has been shown to influence rates of self-pollination amongst plants, with a reduction in herkogamy increasing rates of successful selfing in self-compatible species. Self-pollination is becoming a critical issue in horticultural crops grown in environments where biotic pollinators are limited, absent, or difficult to utilise. In these cases, poor pollination results in reduced yield and misshapen fruit. Whilst there is a growing body of work elucidating the genetic basis of floral organ development, the genetic and environmental control points regulating herkogamy are poorly understood. A better understanding of the developmental and regulatory pathways involved in establishing varying degrees of herkogamy is needed to provide insights into the production of flowers more adept at selfing to produce consistent, high-quality fruit. This review presents our current understanding of herkogamy from a genetics and hormonal perspective.

Heterostyly

Scharman and Lenhard introduce heterostyly, a phenomenon where individuals in a plant population produce flowers with more than one morphologically distinct form.

Molecular insights into self-incompatibility systems: From evolution to breeding

Plants have evolved diverse self-incompatibility (SI) systems for outcrossing. Since Darwin's time, considerable progress has been made toward elucidating this unrivaled reproductive innovation. Recent advances in interdisciplinary studies and applications of biotechnology have given rise to major breakthroughs in understanding the molecular pathways that lead to SI, particularly the strikingly different SI mechanisms that operate in Solanaceae, Papaveraceae, Brassicaceae, and Primulaceae. These best-understood SI systems, together with discoveries in other "nonmodel" SI taxa such as Poaceae, suggest a complex evolutionary trajectory of SI, with multiple independent origins and frequent and irreversible losses. Extensive exploration of self-/nonself-discrimination signaling cascades has revealed a comprehensive catalog of male and female identity genes and modifier factors that control SI. These findings also enable the characterization, validation, and manipulation of SI-related factors for crop improvement, helping to address the challenges associated with development of inbred lines. Here, we review current knowledge about the evolution of SI systems, summarize key achievements in the molecular basis of pollen‒pistil interactions, discuss potential prospects for breeding of SI crops, and raise several unresolved questions that require further investigation.

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.

The genomes of Darwin's primroses reveal chromosome-scale adaptive introgression and differential permeability of species boundaries

Introgression is an important source of genetic variation that can determine species adaptation to environmental conditions. Yet, definitive evidence of the genomic and adaptive implications of introgression in nature remains scarce. The widespread hybrid zones of Darwin's primroses (Primula elatior, Primula veris, and Primula vulgaris) provide a unique natural laboratory for studying introgression in flowering plants and the varying permeability of species boundaries. Through analysis of 650 genomes, we provide evidence of an introgressed genomic region likely to confer adaptive advantage in conditions of soil toxicity. We also document unequivocal evidence of chloroplast introgression, an important precursor to species-wide chloroplast capture. Finally, we provide the first evidence that the S-locus supergene, which controls heterostyly in primroses, does not introgress in this clade. Our results contribute novel insights into the adaptive role of introgression and demonstrate the importance of extensive genomic and geographical sampling for illuminating the complex nature of species boundaries.

Haplotype-resolved genome assembly provides insights into the evolution of S-locus supergene in distylous Nymphoides indica

Distyly has evolved independently in numerous animal-pollinated angiosperm lineages. Understanding of its molecular basis has been restricted to a few species, primarily Primula. Here, we investigate the genetic architecture of the single diallelic locus (S-locus) supergene, a linkage group of functionally associated genes, and explore how it may have evolved in distylous Nymphoides indica, a lineage of flowering plants not previously investigated. We assembled haplotype-resolved genomes, used read-coverage-based genome-wide association study (rb-GWAS) to locate the S-locus supergene, co-expression network analysis to explore gene networks underpinning the development of distyly, and comparative genomic analyses to investigate the origins of the S-locus supergene. We identified three linked candidate S-locus genes - NinBAS1, NinKHZ2, and NinS1 - that were only evident in the short-styled morph and were hemizygous. Co-expression network analysis suggested that brassinosteroids contribute to dimorphic sex organs in the short-styled morph. Comparative genomic analyses indicated that the S-locus supergene likely evolved via stepwise duplications and has been affected by transposable element activities. Our study provides novel insight into the structure, regulation, and evolution of the supergene governing distyly in N. indica. It also provides high-quality genomic resources for future research on the molecular mechanisms underlying the striking evolutionary convergence in form and function across heterostylous taxa.

Brassinosteroid catabolic enzyme CYP734A129 regulates the morphologies of leaves and floral organs in woodland strawberry

Brassinosteroids (BRs) play critical roles in plant growth and development and regulate many important agronomic traits. However, the functions of BRs in strawberry are unclear. This study identified two mutants, named P6 and R87, in woodland strawberry (Fragaria vesca) from EMS mutagenesis populations that exhibit narrow leaves, petals and sepals. Mapping by sequencing and genetic studies revealed that the F. vesca CYP734A129, encoding a putative BR catabolic enzyme, is the causative gene for both P6 and R87. Overexpression of CYP734A129 in both F. vesca and Arabidopsis causes a severe dwarf phenotype, and the BRI1-EMS-SUPPRESSOR 1 (BES1) protein is less abundant in the CYP734A129-overexpressing Arabidopsis seedlings. This suggests that CYP734A129 is functionally conserved with CYP734A1, as a BR-inactivating enzyme. Transcriptome analysis of young leaves revealed that four BR biosynthetic genes were significantly downregulated in P6 (cyp734a129), and photosynthesis-related genes were highly enriched among the up-regulated genes in P6 compared to the wild type. This further supports that CYP734A129 inactivates BRs in F. vesca. Furthermore, we showed that mutations in CYP734A129 do not affect fruit shape and color during ripening in strawberry. Overall, our results suggest that F. vesca CYP734A129 is a BR catabolic enzyme, and provide insights into the roles of CYP734A129 in strawberry.

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.

Genome sequencing reveals the genetic architecture of heterostyly and domestication history of common buckwheat

Common buckwheat, Fagopyrum esculentum, is an orphan crop domesticated in southwest China that exhibits heterostylous self-incompatibility. Here we present chromosome-scale assemblies of a self-compatible F. esculentum accession and a self-compatible wild relative, Fagopyrum homotropicum, together with the resequencing of 104 wild and cultivated F. esculentum accessions. Using these genomic data, we report the roles of transposable elements and whole-genome duplications in the evolution of Fagopyrum. In addition, we show that (1) the breakdown of heterostyly occurs through the disruption of a hemizygous gene jointly regulating the style length and female compatibility and (2) southeast Tibet was involved in common buckwheat domestication. Moreover, we obtained mutants conferring the waxy phenotype for the first time in buckwheat. These findings demonstrate the utility of our F. esculentum assembly as a reference genome and promise to accelerate buckwheat research and breeding.

PbrBZR1 interacts with PbrARI2.3 to mediate brassinosteroid-regulated pollen tube growth during self-incompatibility signaling in pear

S-RNase-mediated self-incompatibility (SI) prevents self-fertilization and promotes outbreeding to ensure genetic diversity in many flowering plants, including pear (Pyrus sp.). Brassinosteroids (BRs) have well-documented functions in cell elongation, but their molecular mechanisms in pollen tube growth, especially in the SI response, remain elusive. Here, exogenously applied brassinolide (BL), an active BR, countered incompatible pollen tube growth inhibition during the SI response in pear. Antisense repression of BRASSINAZOLE-RESISTANT1 (PbrBZR1), a critical component of BR signaling, blocked the positive effect of BL on pollen tube elongation. Further analyses revealed that PbrBZR1 binds to the promoter of EXPANSIN-LIKE A3 (PbrEXLA3) to activate its expression. PbrEXLA3 encodes an expansin that promotes pollen tube elongation in pear. The stability of dephosphorylated PbrBZR1 was substantially reduced in incompatible pollen tubes, where it is targeted by ARIADNE2.3 (PbrARI2.3), an E3 ubiquitin ligase that is strongly expressed in pollen. Our results show that during the SI response, PbrARI2.3 accumulates and negatively regulates pollen tube growth by accelerating the degradation of PbrBZR1 via the 26S proteasome pathway. Together, our results show that an ubiquitin-mediated modification participates in BR signaling in pollen and reveal the molecular mechanism by which BRs regulate S-RNase-based SI.

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.

Brassinosteroid signals cooperate with katanin-mediated microtubule severing to control stamen filament elongation

Proper stamen filament elongation is essential for pollination and plant reproduction. Plant hormones are extensively involved in every stage of stamen development; however, the cellular mechanisms by which phytohormone signals couple with microtubule dynamics to control filament elongation remain unclear. Here, we screened a series of Arabidopsis thaliana mutants showing different microtubule defects and revealed that only those unable to sever microtubules, lue1 and ktn80.1234, displayed differential floral organ elongation with less elongated stamen filaments. Prompted by short stamen filaments and severe decrease in KTN1 and KTN80s expression in qui-2 lacking five BZR1-family transcription factors (BFTFs), we investigated the crosstalk between microtubule severing and brassinosteroid (BR) signaling. The BFTFs transcriptionally activate katanin-encoding genes, and the microtubule-severing frequency was severely reduced in qui-2. Taken together, our findings reveal how BRs can regulate cytoskeletal dynamics to coordinate the proper development of reproductive organs.

Primula luquanensis sp. nov. (Primulaceae), a New Species from Southwestern China, Reveals a Novel Floral Form in the Heterostyly-Prevailing Genus

A new species, Primula luquanensis Z.K.Wu and Wei Zhou sp. nov. (Primulaceae) is described and illustrated from Yunnan Province, China. It is morphologically assigned to P. sect Aleuritia based on its dwarf and hairless habit and coverage by farina on both sides of the leaf blade and scape. This new species is similar to P. nutantiflora and P. yunnanensis, but it is easily distinguished by its stolons, solitary bract, bell-shaped corolla and monomorphic floral form. The new species also has a substantially reduced corolla tube, presenting a unique floral form in a genus where heterostyly typically prevails.

Annotation of the Turnera subulata (Passifloraceae) Draft Genome Reveals the S-Locus Evolved after the Divergence of Turneroideae from Passifloroideae in a Stepwise Manner

A majority of Turnera species (Passifloraceae) exhibit distyly, a reproductive system involving both self-incompatibility and reciprocal herkogamy. This system differs from self-incompatibility in Passiflora species. The genetic basis of distyly in Turnera is a supergene, restricted to the S-morph, and containing three S-genes. How supergenes and distyly evolved in Turnera, and the other Angiosperm families exhibiting distyly remain largely unknown. Unraveling the evolutionary origins in Turnera requires the generation of genomic resources and extensive phylogenetic analyses. Here, we present the annotated draft genome of the S-morph of distylous Turnera subulata. Our annotation allowed for phylogenetic analyses of the three S-genes' families across 56 plant species ranging from non-seed plants to eudicots. In addition to the phylogenetic analysis, we identified the three S-genes' closest paralogs in two species of Passiflora. Our analyses suggest that the S-locus evolved after the divergence of Passiflora and Turnera. Finally, to provide insights into the neofunctionalization of the S-genes, we compared expression patterns of the S-genes with close paralogs in Arabidopsis and Populus trichocarpa. The annotation of the T. subulata genome will provide a useful resource for future comparative work. Additionally, this work has provided insights into the convergent nature of distyly and the origin of supergenes.

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.

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.

Whole-genome analyses disentangle reticulate evolution of primroses in a biodiversity hotspot

Biodiversity hotspots, such as the Caucasus mountains, provide unprecedented opportunities for understanding the evolutionary processes that shape species diversity and richness. Therefore, we investigated the evolution of Primula sect. Primula, a clade with a high degree of endemism in the Caucasus. We performed phylogenetic and network analyses of whole-genome resequencing data from the entire nuclear genome, the entire chloroplast genome, and the entire heterostyly supergene. The different characteristics of the genomic partitions and the resulting phylogenetic incongruences enabled us to disentangle evolutionary histories resulting from tokogenetic vs cladogenetic processes. We provide the first phylogeny inferred from the heterostyly supergene that includes all species of Primula sect. Primula. Our results identified recurrent admixture at deep nodes between lineages in the Caucasus as the cause of non-monophyly in Primula. Biogeographic analyses support the 'out-of-the-Caucasus' hypothesis, emphasizing the importance of this hotspot as a cradle for biodiversity. Our findings provide novel insights into causal processes of phylogenetic discordance, demonstrating that genome-wide analyses from partitions with contrasting genetic characteristics and broad geographic sampling are crucial for disentangling the diversification of species-rich clades in biodiversity hotspots.

Comparative transcriptomics reveals commonalities and differences in the genetic underpinnings of a floral dimorphism

Distyly, a floral dimorphism associated with heteromorphic self-incompatibility and controlled by the S-locus supergene, evolved independently multiple times. Comparative analyses of the first transcriptome atlas for the main distyly model, Primula veris, with other distylous species produced the following findings. A set of 53 constitutively expressed genes in P. veris did not include any of the housekeeping genes commonly used to normalize gene expression in qPCR experiments. The S-locus gene CYP^T acquired its role in controlling style elongation via a change in expression profile. Comparison of genes differentially expressed between floral morphs revealed that brassinosteroids and auxin are the main hormones controlling style elongation in P. veris and Fagopyrum esculentum, respectively. Furthermore, shared biochemical pathways might underlie the expression of distyly in the distantly related P. veris, F. esculentum and Turnera subulata, suggesting a degree of correspondence between evolutionary convergence at phenotypic and molecular levels. Finally, we provide the first evidence supporting the previously proposed hypothesis that distyly supergenes of distantly related species evolved via the recruitment of genes related to the phytochrome-interacting factor (PIF) signaling network. To conclude, this is the first study that discovered homologous genes involved in the control of distyly in distantly related taxa.

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.

Genomic analyses of the Linum distyly supergene reveal convergent evolution at the molecular level

Supergenes govern multi-trait-balanced polymorphisms in a wide range of systems; however, our understanding of their origins and evolution remains incomplete. The reciprocal placement of stigmas and anthers in pin and thrum floral morphs of distylous species constitutes an iconic example of a balanced polymorphism governed by a supergene, the distyly S-locus. Recent studies have shown that the Primula and Turnera distyly supergenes are both hemizygous in thrums, but it remains unknown whether hemizygosity is pervasive among distyly S-loci. As hemizygosity has major consequences for supergene evolution and loss, clarifying whether this genetic architecture is shared among distylous species is critical. Here, we have characterized the genetic architecture and evolution of the distyly supergene in Linum by generating a chromosome-level genome assembly of Linum tenue, followed by the identification of the S-locus using population genomic data. We show that hemizygosity and thrum-specific expression of S-linked genes, including a pistil-expressed candidate gene for style length, are major features of the Linum S-locus. Structural variation is likely instrumental for recombination suppression, and although the non-recombining dominant haplotype has accumulated transposable elements, S-linked genes are not under relaxed purifying selection. Our findings reveal remarkable convergence in the genetic architecture and evolution of independently derived distyly supergenes, provide a counterexample to classic inversion-based supergenes, and shed new light on the origin and maintenance of an iconic floral polymorphism.

Iterative evolution of supergene-based social polymorphism in ants

Species commonly exhibit alternative morphs, with individual fate being determined during development by either genetic factors, environmental cues or a combination thereof. Ants offer an interesting case study because many species are polymorphic in their social structure. Some colonies contain one queen while others contain many queens. This variation in queen number is generally associated with a suite of phenotypic and life-history traits, including mode of colony founding, queen lifespan, queen-worker dimorphism and colony size. The basis of this social polymorphism has been studied in five ant lineages, and remarkably social morph seems to be determined by a supergene in all cases. These 'social supergenes' tend to be large, having formed through serial inversions, and to comprise hundreds of linked genes. They have persisted over long evolutionary timescales, in multiple lineages following speciation events, and have spread between closely related species via introgression. Their evolutionary dynamics are unusually complex, combining recessive lethality, spatially variable selection, selfish genetic elements and non-random mating. Here, we synthesize the five cases of supergene-based social polymorphism in ants, highlighting interesting commonalities, idiosyncrasies and implications for the evolution of polymorphisms in general. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

The origin and evolution of RNase T2 family and gametophytic self-incompatibility system in plants

Ribonuclease (RNase) T2 genes are found widely in both eukaryotes and prokaryotes, and genes from this family have been revealed to have various functions in plants. In particular, S-RNase is known to be the female determinant in the S-RNase-based gametophytic self-incompatibility (GSI) system. However, the origin and evolution of the RNase T2 gene family and GSI system are not well understood. In this study, 785 RNase T2 genes were identified in 81 sequenced plant genomes representing broad-scale diversity and divided into three subgroups (Class I, II, and III) based on phylogenetic and synteny network analysis. Class I was found to be of ancient origin and to emerge in green algae, Class II was shown to originate with the appearance of angiosperms, while Class III was discovered to be eudicot-specific. Each of the three major classes could be further classified into several subclasses of which some subclasses were found to be lineage-specific. Furthermore, duplication, deletion, or inactivation of the S/S-like-locus was revealed to be linked to repeated loss and gain of self-incompatibility in different species from distantly related plant families with GSI. Finally, the origin and evolutionary history of S-locus in Rosaceae species was unraveled with independent loss and gain of S-RNase occurred in different subfamilies of Rosaceae. Our findings provide insights into the origin and evolution of the RNase T2 family and the GSI system in plants.

Comparative Genomics Elucidates the Origin of a Supergene Controlling Floral Heteromorphism

Supergenes are nonrecombining genomic regions ensuring the coinheritance of multiple, coadapted genes. Despite the importance of supergenes in adaptation, little is known on how they originate. A classic example of supergene is the S locus controlling heterostyly, a floral heteromorphism occurring in 28 angiosperm families. In Primula, heterostyly is characterized by the cooccurrence of two complementary, self-incompatible floral morphs and is controlled by five genes clustered in the hemizygous, ca. 300-kb S locus. Here, we present the first chromosome-scale genome assembly of any heterostylous species, that of Primula veris (cowslip). By leveraging the high contiguity of the P. veris assembly and comparative genomic analyses, we demonstrated that the S-locus evolved via multiple, asynchronous gene duplications and independent gene translocations. Furthermore, we discovered a new whole-genome duplication in Ericales that is specific to the Primula lineage. We also propose a mechanism for the origin of S-locus hemizygosity via nonhomologous recombination involving the newly discovered two pairs of CFB genes flanking the S locus. Finally, we detected only weak signatures of degeneration in the S locus, as predicted for hemizygous supergenes. The present study provides a useful resource for future research addressing key questions on the evolution of supergenes in general and the S locus in particular: How do supergenes arise? What is the role of genome architecture in the evolution of complex adaptations? Is the molecular architecture of heterostyly supergenes across angiosperms similar to that of Primula?

Origin, loss, and regain of self-incompatibility in angiosperms

The self-incompatibility (SI) system with the broadest taxonomic distribution in angiosperms is based on multiple S-locus F-box genes (SLFs) tightly linked to an S-RNase termed type-1. Multiple SLFs collaborate to detoxify nonself S-RNases while being unable to detoxify self S-RNases. However, it is unclear how such a system evolved, because in an ancestral system with a single SLF, many nonself S-RNases would not be detoxified, giving low cross-fertilization rates. In addition, how the system has been maintained in the face of whole-genome duplications (WGDs) or lost in other lineages remains unclear. Here we show that SLFs from a broad range of species can detoxify S-RNases from Petunia with a high detoxification probability, suggestive of an ancestral feature enabling cross-fertilization and subsequently modified as additional SLFs evolved. We further show, based on its genomic signatures, that type-1 was likely maintained in many lineages, despite WGD, through deletion of duplicate S-loci. In other lineages, SI was lost either through S-locus deletions or by retaining duplications. Two deletion lineages regained SI through type-2 (Brassicaceae) or type-4 (Primulaceae), and one duplication lineage through type-3 (Papaveraceae) mechanisms. Thus, our results reveal a highly dynamic process behind the origin, maintenance, loss, and regain of SI.

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

Female self-incompatibility type in heterostylous Primula is determined by the brassinosteroid-inactivating cytochrome P450 CYP734A50

Most flowering plants are hermaphrodites, with flowers having both male and female reproductive organs. One widespread adaptation to limit self-fertilization is self-incompatibility (SI), where self-pollen fails to fertilize ovules.^1,2 In homomorphic SI, many morphologically indistinguishable mating types are found, although in heteromorphic SI, the two or three mating types are associated with different floral morphologies.^3-6 In heterostylous Primula, a hemizygous supergene determines a short-styled S-morph and a long-styled L-morph, corresponding to two different mating types, and full seed set only results from intermorph crosses.^7-9 Style length is controlled by the brassinosteroid (BR)-inactivating cytochrome P450 CYP734A50,¹⁰ yet it remains unclear what defines the male and female incompatibility types. Here, we show that CYP734A50 also determines the female incompatibility type. Inactivating CYP734A50 converts short S-morph styles into long styles with the same incompatibility behavior as L-morph styles, and this effect can be mimicked by exogenous BR treatment. In vitro responses of S- and L-morph pollen grains and pollen tubes to increasing BR levels could only partly explain their different in vivo behavior, suggesting both direct and indirect effects of the different BR levels in S- versus L-morph stigmas and styles in controlling pollen performance. This BR-mediated SI provides a novel mechanism for preventing self-fertilization. The joint control of morphology and SI by CYP734A50 has important implications for the evolutionary buildup of the heterostylous syndrome and provides a straightforward explanation for why essentially all of the derived self-compatible homostylous Primula species are long homostyles.¹¹.

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

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

Pistil Mating Type and Morphology Are Mediated by the Brassinosteroid Inactivating Activity of the S-Locus Gene BAHD in Heterostylous Turnera Species

Heterostyly is a breeding system that promotes outbreeding through a combination of morphological and physiological floral traits. In Turnera these traits are governed by a single, hemizygous S-locus containing just three genes. We report that the S-locus gene, BAHD, is mutated and encodes a severely truncated protein in a self-compatible long homostyle species. Further, a self-compatible long homostyle mutant possesses a T. krapovickasii BAHD allele with a point mutation in a highly conserved domain of BAHD acyl transferases. Wild type and mutant TkBAHD alleles were expressed in Arabidopsis to assay for brassinosteroid (BR) inactivating activity. The wild type but not mutant allele caused dwarfism, consistent with the wild type possessing, but the mutant allele having lost, BR inactivating activity. To investigate whether BRs act directly in self-incompatibility, BRs were added to in vitro pollen cultures of the two mating types. A small morph specific stimulatory effect on pollen tube growth was found with 5 µM brassinolide, but no genotype specific inhibition was observed. These results suggest that BAHD acts pleiotropically to mediate pistil length and physiological mating type through BR inactivation, and that in regard to self-incompatibility, BR acts by differentially regulating gene expression in pistils, rather than directly on pollen.

A Draft Genome of the Ginger Species Alpinia nigra and New Insights into the Genetic Basis of Flexistyly

Angiosperms possess various strategies to ensure reproductive success, such as stylar polymorphisms that encourage outcrossing. Here, we investigate the genetic basis of one such dimorphism that combines both temporal and spatial separation of sexual function, termed flexistyly. It is a floral strategy characterised by the presence of two morphs that differ in the timing of stylar movement. We performed a de novo assembly of the genome of Alpinia nigra using high-depth genomic sequencing. We then used Pool-seq to identify candidate regions for flexistyly based on allele frequency or coverage differences between pools of anaflexistylous and cataflexistylous morphs. The final genome assembly size was 2 Gb, and showed no evidence of recent polyploidy. The Pool-seq did not reveal large regions with high FST values, suggesting large structural chromosomal polymorphisms are unlikely to underlie differences between morphs. Similarly, no region had a 1:2 mapping depth ratio which would be indicative of hemizygosity. We propose that flexistyly is governed by a small genomic region that might be difficult to detect with Pool-seq, or a complex genomic region that proved difficult to assemble. Our genome will be a valuable resource for future studies of gingers, and provides the first steps towards characterising this complex floral phenotype.

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.

The Genomic Architecture and Evolutionary Fates of Supergenes

Supergenes are genomic regions containing sets of tightly linked loci that control multi-trait phenotypic polymorphisms under balancing selection. Recent advances in genomics have uncovered significant variation in both the genomic architecture as well as the mode of origin of supergenes across diverse organismal systems. Although the role of genomic architecture for the origin of supergenes has been much discussed, differences in the genomic architecture also subsequently affect the evolutionary trajectory of supergenes and the rate of degeneration of supergene haplotypes. In this review, we synthesize recent genomic work and historical models of supergene evolution, highlighting how the genomic architecture of supergenes affects their evolutionary fate. We discuss how recent findings on classic supergenes involved in governing ant colony social form, mimicry in butterflies, and heterostyly in flowering plants relate to theoretical expectations. Furthermore, we use forward simulations to demonstrate that differences in genomic architecture affect the degeneration of supergenes. Finally, we discuss implications of the evolution of supergene haplotypes for the long-term fate of balanced polymorphisms governed by supergenes.

Asymmetries of reproductive isolation are reflected in directionalities of hybridization: integrative evidence on the complexity of species boundaries

The complex nature of species boundaries has been a central topic in evolutionary biology ever since Darwin. Despite numerous separate studies on reproductive isolation and hybridization, their relationship remains underinvestigated. Are the strengths and asymmetries of reproductive barriers reflected in the extent and directionalities of interspecific genetic exchange? We combined field, experimental, and molecular data to quantify strengths and asymmetries of sympatric reproductive barriers and hybridization between florally heteromorphic primroses. We also assessed whether generalist pollinators discriminate between different floral cues and contribute to reproductive isolation, a long-debated topic. Sympatric reproductive isolation is high but incomplete, and most phenotypic intermediates are genetic F₁ hybrids, whereas backcrosses are rare, revealing low interspecific gene flow. Species integrity rests on multiple barriers, but ethological isolation is among the strongest, demonstrating that even generalist pollinators crucially contribute to the maintenance of species boundaries. Furthermore, reproductive barriers are weaker for Primula veris and short-styled plants, results corroborated by molecular data. Thus, in florally heteromorphic systems, both species- and morph-dependent asymmetries affect permeability of species boundaries. Our study illustrates how the interactions between complex floral syndromes and pollinators shape species boundaries in unique, previously undescribed ways.

Characterization of 30 microsatellite markers for distylous Primula denticulata (Primulaceae) using HiSeq sequencing

Primula denticulata exhibits considerable variation in floral morphology and flowering phenology along elevational gradients in SW China. We isolated 30 microsatellite markers from P. denticulata to facilitate further investigation of population genetics and floral evolution in this species. We used the HiSeq X-Ten sequencing system to develop a set of markers, and measured polymorphism and genetic diversity in a sample of 72 individuals from three natural populations of P. denticulata subsp. denticulata. The markers displayed relatively high polymorphism, with the number of alleles ranging from two to seven (mean = 3.567). The observed and expected heterozygosity ranged from 0 to 1.000 and 0.041 to 0.702, respectively. Twenty-eight of the loci were also successfully amplified in P. denticulata subsp. sinodenticulata. The microsatellite markers we have identified will provide valuable tools for investigations of the population genetic structure, mating systems and phylogeography of the P. denticulata complex, and will help to address questions concerning the ecological and genetic mechanisms responsible for the evolution of reproductive traits in the species.