High-resolution mapping of the S-locus in Turnera leads to the discovery of three genes tightly associated with the S-alleles

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.

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.

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.

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

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

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

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.

The multi-protein family of sulfotransferases in plants: composition, occurrence, substrate specificity, and functions

All members of the sulfotransferase (SOT, EC 2.8.2.-) protein family transfer a sulfuryl group from the donor 3'-phosphoadenosine 5'-phosphosulfate (PAPS) to an appropriate hydroxyl group of several classes of substrates. The primary structure of these enzymes is characterized by a histidine residue in the active site, defined PAPS binding sites and a longer SOT domain. Proteins with this SOT domain occur in all organisms from all three domains, usually as a multi-protein family. Arabidopsis thaliana SOTs, the best characterized SOT multi-protein family, contains 21 members. The substrates for several plant enzymes have already been identified, such as glucosinolates, brassinosteroids, jasmonates, flavonoids, and salicylic acid. Much information has been gathered on desulfo-glucosinolate (dsGl) SOTs in A. thaliana. The three cytosolic dsGl SOTs show slightly different expression patterns. The recombinant proteins reveal differences in their affinity to indolic and aliphatic dsGls. Also the respective recombinant dsGl SOTs from different A. thaliana ecotypes differ in their kinetic properties. However, determinants of substrate specificity and the exact reaction mechanism still need to be clarified. Probably, the three-dimensional structures of more plant proteins need to be solved to analyze the mode of action and the responsible amino acids for substrate binding. In addition to A. thaliana, more plant species from several families need to be investigated to fully elucidate the diversity of sulfated molecules and the way of biosynthesis catalyzed by SOT enzymes.

Changes in nucleosome position at transcriptional start sites of specific genes in Zea mays mediator of paramutation1 mutants

Nucleosomes facilitate compaction of DNA within the confines of the eukaryotic nucleus. This packaging of DNA and histone proteins must accommodate cellular processes, such as transcription and DNA replication. The repositioning of nucleosomes to facilitate cellular processes is likely regulated by several factors. In Zea mays, Mediator of paramutation1 (MOP1) has been demonstrated to be an epigenetic regulator of gene expression. Based on sequence orthology and mutant phenotypes, MOP1 is likely to function in an RNA-dependent pathway to mediate changes to chromatin. High-resolution microarrays were used to assay the distribution of nucleosomes across the transcription start sites (TSSs) of ~400 maize genes in wild type and mutant mop1-1 tissues. Analysis of nucleosome distribution in leaf, immature tassel and ear shoot tissues resulted in the identification of three genes showing consistent differences in nucleosome positioning and occupancy between wild type and mutant mop1-1. These specific changes in nucleosome distribution were located upstream as well as downstream of the TSS. No direct relationship between the specific changes in nucleosome distribution and transcription were observed through quantitative expression analysis in these tissues. In silico prediction suggests that nucleosome positioning is not dictated by intrinsic DNA sequence signals in the TSSs of two of the identified genes, suggesting a role for chromatin remodeling proteins in MOP1-mediated pathways. These results also indicate that MOP1 contributions to nucleosome position may be either separate from changes in gene expression, or cooperative with development and other levels of regulation in coordinating gene expression.

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.

The genetic architecture of natural variation in flower morphology

A pollination syndrome is defined as a suite of floral traits that are associated with the attraction of a specific group of animals as pollinators. Traits such as flower morphology, color, scent, and rewards contribute to the plant's reproductive success by attracting pollinators. Here we focus on the genetics of natural variation in flower morphology and how the adaptation between plants and their cognate pollinator class contributes to plant's reproductive success. We review recent work on the genetic basis of interspecific differences in reproductive organ morphology and discuss possible genetic mechanisms for coordinated changes in complex syndromes.

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

Heterostyly is a plant breeding system occurring in approximately 28 plant families and it has often been used as a model system in plant genetics and evolution. Although heterostyly has been studied for over a century beginning with Charles Darwin, the genes determining floral architecture and incompatibility are still unknown. To identify the genes residing at the S-locus of distylous Turnera subulata, we used a positional cloning strategy and assembled three BAC contigs across the S-locus region. In total, 31 overlapping BAC clones were assembled into contigs 1, 2 and SL. We developed and mapped numerous co-dominant markers from the ends of BAC clones across the S-locus region and assayed X-ray deletion mutants to delimit the region of the contig containing the S-locus. Deletion mapping revealed that a single BAC clone (L22s) within contig-SL contains the s haplotype, while two additional BAC clones (I1 and K15) may contain parts of the dominant S haplotype. Furthermore, we exploited the contigs assembled and investigated the rates of recombination at the S-locus as well as in two regions on either side of the S-locus. We found that recombination rates (estimated in kb/cM) are 2-5 times lower at the S-locus relative to flanking regions, although they are not statistically significant. The present study represents a landmark in the molecular characterization of the S-locus of a heterostylous species. We are now on the verge of identifying the genes that have remained elusive since Darwin's comprehensive study of heterostylous systems more than 125 years ago.

Darwin and genetics

Darwin's theory of natural selection lacked an adequate account of inheritance, making it logically incomplete. We review the interaction between evolution and genetics, showing how, unlike Mendel, Darwin's lack of a model of the mechanism of inheritance left him unable to interpret his own data that showed Mendelian ratios, even though he shared with Mendel a more mathematical and probabilistic outlook than most biologists of his time. Darwin's own "pangenesis" model provided a mechanism for generating ample variability on which selection could act. It involved, however, the inheritance of characters acquired during an organism's life, which Darwin himself knew could not explain some evolutionary situations. Once the particulate basis of genetics was understood, it was seen to allow variation to be passed intact to new generations, and evolution could then be understood as a process of changes in the frequencies of stable variants. Evolutionary genetics subsequently developed as a central part of biology. Darwinian principles now play a greater role in biology than ever before, which we illustrate with some examples of studies of natural selection that use DNA sequence data and with some recent advances in answering questions first asked by Darwin.