Positional cloning of the s haplotype determining the floral and incompatibility phenotype of the long-styled morph of distylous Turnera subulata

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.

Transcriptome and Network Analyses of Heterostyly in Turnera subulata Provide Mechanistic Insights: Are S-Loci a Red-Light for Pistil Elongation?

Heterostyly employs distinct hermaphroditic floral morphs to enforce outbreeding. Morphs differ structurally in stigma/anther positioning, promoting cross-pollination, and physiologically blocking self-fertilization. Heterostyly is controlled by a self-incompatibility (S)-locus of a small number of linked S-genes specific to short-styled morph genomes. Turnera possesses three S-genes, namely TsBAHD (controlling pistil characters), TsYUC6, and TsSPH1 (controlling stamen characters). Here, we compare pistil and stamen transcriptomes of floral morphs of T. subulata to investigate hypothesized S-gene function(s) and whether hormonal differences might contribute to physiological incompatibility. We then use network analyses to identify genetic networks underpinning heterostyly. We found a depletion of brassinosteroid-regulated genes in short styled (S)-morph pistils, consistent with hypothesized brassinosteroid-inactivating activity of TsBAHD. In S-morph anthers, auxin-regulated genes were enriched, consistent with hypothesized auxin biosynthesis activity of TsYUC6. Evidence was found for auxin elevation and brassinosteroid reduction in both pistils and stamens of S- relative to long styled (L)-morph flowers, consistent with reciprocal hormonal differences contributing to physiological incompatibility. Additional hormone pathways were also affected, however, suggesting S-gene activities intersect with a signaling hub. Interestingly, distinct S-genes controlling pistil length, from three species with independently evolved heterostyly, potentially intersect with phytochrome interacting factor (PIF) network hubs which mediate red/far-red light signaling. We propose that modification of the activities of PIF hubs by the S-locus could be a common theme in the evolution of heterostyly.

The long and short of the S-locus in Turnera (Passifloraceae)

Distyly is an intriguing floral adaptation that increases pollen transfer precision and restricts inbreeding. It has been a model system in evolutionary biology since Darwin. Although the S-locus determines the long- and short-styled morphs, the genes were unknown in Turnera. We have now identified these genes. We used deletion mapping to identify, and then sequence, BAC clones and genome scaffolds to construct S/s haplotypes. We investigated candidate gene expression, hemizygosity, and used mutants, to explore gene function. The s-haplotype possessed 21 genes collinear with a region of chromosome 7 of grape. The S-haplotype possessed three additional genes and two inversions. TsSPH1 was expressed in filaments and anthers, TsYUC6 in anthers and TsBAHD in pistils. Long-homostyle mutants did not possess TsBAHD and a short-homostyle mutant did not express TsSPH1. Three hemizygous genes appear to determine S-morph characteristics in T. subulata. Hemizygosity is common to all distylous species investigated, yet the genes differ. The pistil candidate gene, TsBAHD, differs from that of Primula, but both may inactivate brassinosteroids causing short styles. TsYUC6 is involved in auxin synthesis and likely determines pollen characteristics. TsSPH1 is likely involved in filament elongation. We propose an incompatibility mechanism involving TsYUC6 and TsBAHD.

Primula vulgaris (primrose) genome assembly, annotation and gene expression, with comparative genomics on the heterostyly supergene

Primula vulgaris (primrose) exhibits heterostyly: plants produce self-incompatible pin- or thrum-form flowers, with anthers and stigma at reciprocal heights. Darwin concluded that this arrangement promotes insect-mediated cross-pollination; later studies revealed control by a cluster of genes, or supergene, known as the S (Style length) locus. The P. vulgaris S locus is absent from pin plants and hemizygous in thrum plants (thrum-specific); mutation of S locus genes produces self-fertile homostyle flowers with anthers and stigma at equal heights. Here, we present a 411 Mb P. vulgaris genome assembly of a homozygous inbred long homostyle, representing ~87% of the genome. We annotate over 24,000 P. vulgaris genes, and reveal more genes up-regulated in thrum than pin flowers. We show reduced genomic read coverage across the S locus in other Primula species, including P. veris, where we define the conserved structure and expression of the S locus genes in thrum. Further analysis reveals the S locus has elevated repeat content (64%) compared to the wider genome (37%). Our studies suggest conservation of S locus genetic architecture in Primula, and provide a platform for identification and evolutionary analysis of the S locus and downstream targets that regulate heterostyly in diverse heterostylous species.

A short story gets longer: recent insights into the molecular basis of heterostyly

Heterostyly is a fascinating adaptation to promote outbreeding and a classical paradigm of botany. In the most common type of heterostyly, plants either form flowers with long styles and short stamens, or short styles and long stamens. This reciprocal organ positioning reduces pollen wastage and promotes cross-pollination, thus increasing male fitness. In addition, in many heterostylous species selfing and the generation of unfit progeny due to inbreeding depression is limited by a self-incompatibility system, thus promoting female fitness. The two floral forms are genetically determined by the S locus as a complex supergene, namely a chromosomal region containing several individual genes that control the different traits, such as style or stamen length, and are held together by very tight linkage due to suppressed recombination. Recent molecular-genetic studies in several systems, including Turnera, Fagopyrum, Linum, and Primula have begun to identify and characterize the causal heterostyly genes residing at the S locus. An emerging theme from several families is that the dominant S haplotype represents a hemizygous region not present on the recessive s haplotype. This provides an explanation for the suppressed recombination and suggests a scenario for the chromosomal evolution of the S locus. In this review, we discuss the results from recent molecular-genetic analyses in light of the classical models on the genetics and evolution of heterostyly.

Sequencing the genomic regions flanking S-linked PvGLO sequences confirms the presence of two GLO loci, one of which lies adjacent to the style-length determinant gene CYP734A50

Primula vulgaris contains two GLOBOSA loci, one located adjacent to the style length determinant gene CYP734A50 which lies within the S -locus. Using a combination of BAC walking and PacBio sequencing, we have sequenced two substantial genomic contigs in and around the S-locus of Primula vulgaris. Using these data, we were able to demonstrate that two alleles of PvGlo ^P as well as PvGlo ^T can be present in the genome of a single plant, providing empirical evidence that these two forms of the MADS-box gene GLOBOSA are separate loci and not allelic as previously reported. We propose they should be renamed PvGLO1 and PvGLO2. BAC contigs extending from each GLOBOSA locus were identified and fully sequenced. No homologous genes were found between the contigs other than the GLOBOSA genes themselves, consistent with their identity as separate loci. Exons of the recently identified style-length determinant gene CYP734A50 were identified on one end of the contig containing PvGLO2 and these genes are adjacent in the genome, suggesting that PvGLO2 lies either within or at least very close to the S-locus. Current evidence suggests that both CYP734A50 and GLO2 are specific to the S-morph mating type and are hemizygous rather than heterozygous in the Primula genome. This finding contrasts classical models of the HSI locus, which propose that components of the S-locus are allelic, suggesting that these models may need to be reconsidered.

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

Heterostyly is a wide-spread floral adaptation to promote outbreeding, yet its genetic basis and evolutionary origin remain poorly understood. In Primula (primroses), heterostyly is controlled by the S-locus supergene that determines the reciprocal arrangement of reproductive organs and incompatibility between the two morphs. However, the identities of the component genes remain unknown. Here, we identify the Primula CYP734A50 gene, encoding a putative brassinosteroid-degrading enzyme, as the G locus that determines the style-length dimorphism. CYP734A50 is only present on the short-styled S-morph haplotype, it is specifically expressed in S-morph styles, and its loss or inactivation leads to long styles. The gene arose by a duplication specific to the Primulaceae lineage and shows an accelerated rate of molecular evolution. Thus, our results provide a mechanistic explanation for the Primula style-length dimorphism and begin to shed light on the evolution of the S-locus as a prime model for a complex plant supergene.

DOI:http://dx.doi.org/10.7554/eLife.17956.001

Flowers are highly specialized structures that many plants use to reproduce. Male organs called stamens on the flowers make pollen that can be transferred – usually by insect carriers or the wind – to a female structure called the stigma on another plant. However, since many flowers contain both male and female organs, it is also possible for the pollen to land on the stigma of the same flower, leading to a process called “self-fertilization”.

Many plants have developed mechanisms that prevent self-fertilization. For example, primroses produce two different types of flowers that arrange their stamens and stigmas differently. The stigma sits on the top of a stalk known as the style. Some primroses produce flowers with short stamens and a long style, resulting in the stigma being located high up in the flower (“pin” flowers), while others produce flowers with a short style and long stamens (“thrum” flowers). Primrose pollen is carried by insects and the different lengths of the styles and stamens make it more likely that pollen from a pin flower will land on the stigma of a thrum flower instead of a pin flower (and vice versa).

Although primrose flowers have fascinated botanists for centuries, the genes responsible for making the two types of flower had not been identified. Genetic studies indicated that different genes control the length of the stamens and style. However, these genes appear to be very close to each other on primrose DNA, which made it difficult to study them individually.

Huu et al. identified a gene called CYP734A50 that is responsible for the difference in style length in the flowers of a primrose called Primula veris. The gene is only present in the plants that have thrum flowers across a wide range of primrose species and genetic mutations that inactivate the gene in these plants result in flowers with longer styles. CYP734A50 encodes an enzyme that breaks down plant hormones called brassinosteroids, which normally promote growth. Treating thrum flowers with brassinosteroids increased the length of the styles. Future challenges are to identify the other genes that are responsible for producing pin and thrum flowers and to understand how this group of genes evolved.

DOI:http://dx.doi.org/10.7554/eLife.17956.002

A new physical mapping approach refines the sex-determining gene positions on the Silene latifolia Y-chromosome

Sex chromosomes are particularly interesting regions of the genome for both molecular genetics and evolutionary studies; yet, for most species, we lack basic information, such as the gene order along the chromosome. Because they lack recombination, Y-linked genes cannot be mapped genetically, leaving physical mapping as the only option for establishing the extent of synteny and homology with the X chromosome. Here, we developed a novel and general method for deletion mapping of non-recombining regions by solving “the travelling salesman problem”, and evaluate its accuracy using simulated datasets. Unlike the existing radiation hybrid approach, this method allows us to combine deletion mutants from different experiments and sources. We applied our method to a set of newly generated deletion mutants in the dioecious plant Silene latifolia and refined the locations of the sex-determining loci on its Y chromosome map.

Oakleaf: an S locus-linked mutation of Primula vulgaris that affects leaf and flower development

In Primula vulgaris outcrossing is promoted through reciprocal herkogamy with insect-mediated cross-pollination between pin and thrum form flowers. Development of heteromorphic flowers is coordinated by genes at the S locus. To underpin construction of a genetic map facilitating isolation of these S locus genes, we have characterised Oakleaf, a novel S locus-linked mutant phenotype. We combine phenotypic observation of flower and leaf development, with classical genetic analysis and next-generation sequencing to address the molecular basis of Oakleaf. Oakleaf is a dominant mutation that affects both leaf and flower development; plants produce distinctive lobed leaves, with occasional ectopic meristems on the veins. This phenotype is reminiscent of overexpression of Class I KNOX-homeodomain transcription factors. We describe the structure and expression of all eight P. vulgaris PvKNOX genes in both wild-type and Oakleaf plants, and present comparative transcriptome analysis of leaves and flowers from Oakleaf and wild-type plants. Oakleaf provides a new phenotypic marker for genetic analysis of the Primula S locus. We show that none of the Class I PvKNOX genes are strongly upregulated in Oakleaf leaves and flowers, and identify cohorts of 507 upregulated and 314 downregulated genes in the Oakleaf mutant.

Integration of genetic and physical maps of the Primula vulgaris S locus and localization by chromosome in situ hybridization

Heteromorphic flower development in Primula is controlled by the S locus. The S locus genes, which control anther position, pistil length and pollen size in pin and thrum flowers, have not yet been characterized. We have integrated S-linked genes, marker sequences and mutant phenotypes to create a map of the P. vulgaris S locus region that will facilitate the identification of key S locus genes. We have generated, sequenced and annotated BAC sequences spanning the S locus, and identified its chromosomal location. We have employed a combination of classical genetics and three-point crosses with molecular genetic analysis of recombinants to generate the map. We have characterized this region by Illumina sequencing and bioinformatic analysis, together with chromosome in situ hybridization. We present an integrated genetic and physical map across the P. vulgaris S locus flanked by phenotypic and DNA sequence markers. BAC contigs encompass a 1.5-Mb genomic region with 1 Mb of sequence containing 82 S-linked genes anchored to overlapping BACs. The S locus is located close to the centromere of the largest metacentric chromosome pair. These data will facilitate the identification of the genes that orchestrate heterostyly in Primula and enable evolutionary analyses of the S locus.

On the origins of observations of heterostyly in Primula

In 1862, Charles Darwin published his landmark study on the different forms of flower in Primula; he coined the term distyly and subsequently expanded his studies to other species, including those with tristyly. Darwin is widely recognized as the first to study pin and thrum flowers in Primula, and to provide an explanation for the functional significance of the two floral morphs. Our laboratory is pursuing the genes that underpin floral heteromorphy in Primula, work influenced by Darwin's observations. One day, while appreciating a print of Primula vulgaris from William Curtis' Flora Londinensis, I was struck by the fact that I was looking at images of dimorphic Primula flowers captured in a late-1700s copper-plate engraving that predated Darwin's observations by over 70 yr. This realization triggered a journey into archives of botanical texts, herbals and florilegia from the 16(th) to 19(th) Centuries, and correspondence archives, in search of earlier documents that could have influenced Darwin and the origins of an idea. Darwin was not the first to observe floral heteromorphy in Primula, but he was the first to realize the significance of the two floral morphs. Darwin's insight and exposition of purpose have underpinned all consequent work on the subject.

The status of supergenes in the 21st century: recombination suppression in Batesian mimicry and sex chromosomes and other complex adaptations

I review theoretical models for the evolution of supergenes in the cases of Batesian mimicry in butterflies, distylous plants and sex chromosomes. For each of these systems, I outline the genetic evidence that led to the proposal that they involve multiple genes that interact during ‘complex adaptations’, and at which the mutations involved are not unconditionally advantageous, but show advantages that trade‐off against some disadvantages. I describe recent molecular genetic studies of these systems and questions they raise about the evolution of suppressed recombination. Nonrecombining regions of sex chromosomes have long been known, but it is not yet fully understood why recombination suppression repeatedly evolved in systems in distantly related taxa, but does not always evolve. Recent studies of distylous plants are tending to support the existence of recombination‐suppressed genome regions, which may include modest numbers of genes and resemble recently evolved sex‐linked regions. For Batesian mimicry, however, molecular genetic work in two butterfly species suggests a new supergene scenario, with a single gene mutating to produce initial adaptive phenotypes, perhaps followed by modifiers specifically refining and perfecting the new phenotype.

Supergenes and their role in evolution

Adaptation is commonly a multidimensional problem, with changes in multiple traits required to match a complex environment. This is epitomized by balanced polymorphisms in which multiple phenotypes co-exist and are maintained in a population by a balance of selective forces. Consideration of such polymorphisms led to the concept of the supergene, where alternative phenotypes in a balanced polymorphism segregate as if controlled by a single genetic locus, resulting from tight genetic linkage between multiple functional loci. Recently, the molecular basis for several supergenes has been resolved. Thus, major chromosomal inversions have been shown to be associated with polymorphisms in butterflies, ants and birds, offering a mechanism for localised reduction in recombination. In several examples of plant self-incompatibility, the functional role of multiple elements within the supergene architecture has been demonstrated, conclusively showing that balanced polymorphism can be maintained at multiple coadapted and tightly linked elements. Despite recent criticism, we argue that the supergene concept remains relevant and is more testable than ever with modern molecular methods.

S-LOCUS EARLY FLOWERING 3 is exclusively present in the genomes of short-styled buckwheat plants that exhibit heteromorphic self-incompatibility

The different forms of flowers in a species have attracted the attention of many evolutionary biologists, including Charles Darwin. In Fagopyrum esculentum (common buckwheat), the occurrence of dimorphic flowers, namely short-styled and long-styled flowers, is associated with a type of self-incompatibility (SI) called heteromorphic SI. The floral morphology and intra-morph incompatibility are both determined by a single genetic locus named the S-locus. Plants with short-styled flowers are heterozygous (S/s) and plants with long-styled flowers are homozygous recessive (s/s) at the S-locus. Despite recent progress in our understanding of the molecular basis of flower development and plant SI systems, the molecular mechanisms underlying heteromorphic SI remain unresolved. By examining differentially expressed genes from the styles of the two floral morphs, we identified a gene that is expressed only in short-styled plants. The novel gene identified was completely linked to the S-locus in a linkage analysis of 1,373 plants and had homology to EARLY FLOWERING 3. We named this gene S-LOCUS EARLY FLOWERING 3 (S-ELF3). In an ion-beam-induced mutant that harbored a deletion in the genomic region spanning S-ELF3, a phenotype shift from short-styled flowers to long-styled flowers was observed. Furthermore, S-ELF3 was present in the genome of short-styled plants and absent from that of long-styled plants both in world-wide landraces of buckwheat and in two distantly related Fagopyrum species that exhibit heteromorphic SI. Moreover, independent disruptions of S-ELF3 were detected in a recently emerged self-compatible Fagopyrum species and a self-compatible line of buckwheat. The nonessential role of S-ELF3 in the survival of individuals and the prolonged evolutionary presence only in the genomes of short-styled plants exhibiting heteromorphic SI suggests that S-ELF3 is a suitable candidate gene for the control of the short-styled phenotype of buckwheat plants.

Isolation of the floral morph-related genes in heterostylous flax (Linum grandiflorum): the genetic polymorphism and the transcriptional and post-transcriptional regulations of the S locus

Heterostylous species have two types of flowers, thrum and pin morphs, and these are controlled by a single diallelic locus designated the 'S locus'; fertilization between these two types of flowers is successful. The S gene and the molecular mechanism by which it operates remain to be uncovered, although heterostyly has been studied since the time of Darwin. We compared transcripts and proteins of the thrum and pin flowers of heterostylous flax (Linum grandiflorum) to characterize the molecular differences between them and to elucidate the molecular machinery of heterostyly. Twelve floral morph-related genes were eventually isolated by an integrated study of subtraction and 2D-PAGE analyses, and four genes, TSS1, LgAP1, LgMYB21 and LgSKS1, were predicted to be related to heterostyly. TSS1, a thrum style-specific gene, showed some features suitable for the S gene. Although its biological function is unclear, TSS1 was expressed only in the thrum style and is probably linked to the S locus. LgMYB21, another thrum style gene, would be involved in floral morphogenesis. LgMYB21 was highly expressed in the thrum style, which is shorter than the pin style, and its overexpression in Arabidopsis reduced pistil length. Furthermore, a comparison of transcript and protein accumulations showed no differences in the mRNA accumulation of some thrum-specific proteins, including LgSKS1, suggesting that these are regulated by floral morph-specific post-transcriptional regulation. The Linum S locus regulates not only S specificity but also many floral phenotypes. Dynamic regulation of transcripts and proteins would be necessary for the pleiotropic function of the Linum S locus.

Floral heteromorphy in Primula vulgaris: progress towards isolation and characterization of the S locus

BACKGROUND

The common primrose, Primula vulgaris, along with many other species of the Primulaceae, exhibits floral heteromorphy in which different individuals develop one of two possible forms of flower, known as pin and thrum. Both flower types are hermaphrodite and exhibit reciprocal positions of male and female reproductive structures, which together with a sporophytic incompatibility system, prevent self-pollination and promote out-crossing. The development of the two different forms of flower is controlled by a co-adapted linkage group of genes known as the S locus.

SCOPE

Here progress towards identification and characterization of these genes is described to provide a molecular genetic explanation of the different floral characteristics that define heterostyly in Primula as observed and described by Charles Darwin. Previous work to identify and characterize developmental mutations linked to the P. vulgaris S locus, together with the isolation of S locus-linked genes and polymorphic DNA sequences markers, is summarized. The development of tools are described which will facilitate isolation and characterization of the S locus and its environs, including the creation of two expressed sequence tag libraries from pin and thrum flowers, as well as the construction and screening of two bacterial artificial chromosome (BAC) libraries containing thrum genomic DNA. Screening of these libraries with four S locus-linked sequences has enabled us to assemble four BAC contigs representing over 40 individual overlapping BAC clones which represent over 2·2 Mb of S locus-linked genomic sequence. PCR-based approaches for identification of the allelic origin of these BACs are described as well as identification of an additional 14 S locus-linked genes within BAC-end sequences.

CONCLUSIONS

On-going work to assemble the four S locus-linked contigs into one contiguous sequence spanning the S locus is outlined in preparation for sequence analysis and characterization of the genes located within this region.