Evolution under tight linkage to mating type

The Evolution of Extreme Genetic Variability in a Parasite-Resistance Complex

Genomic regions that play a role in parasite defense are often found to be highly variable, with the major histocompatibility complex serving as an iconic example. Single nucleotide polymorphisms may represent only a small portion of this variability, with Indel polymorphisms and copy number variation further contributing. In extreme cases, haplotypes may no longer be recognized as orthologous. Understanding the evolution of such highly divergent regions is challenging because the most extreme variation is not visible using reference-assisted genomic approaches. Here we analyze the case of the Pasteuria Resistance Complex in the crustacean Daphnia magna, a defense complex in the host against the common and virulent bacterium Pasteuria ramosa. Two haplotypes of this region have been previously described, with parts of it being nonhomologous, and the region has been shown to be under balancing selection. Using pan-genome analysis and tree reconciliation methods to explore the evolution of the Pasteuria Resistance Complex and its characteristics within and between species of Daphnia and other Cladoceran species, our analysis revealed a remarkable diversity in this region even among host species, with many nonhomologous hyper-divergent haplotypes. The Pasteuria Resistance Complex is characterized by extensive duplication and losses of Fucosyltransferase (FuT) and Galactosyltransferase (GalT) genes that are believed to play a role in parasite defense. The Pasteuria Resistance Complex region can be traced back to common ancestors over 250 million years. The unique combination of an ancient resistance complex and a dynamic, hyper-divergent genomic environment presents a fascinating opportunity to investigate the role of such regions in the evolution and long-term maintenance of resistance polymorphisms. Our findings offer valuable insights into the evolutionary forces shaping disease resistance and adaptation, not only in the genus Daphnia, but potentially across the entire Cladocera class.

The genetic architecture of the load linked to dominant and recessive self-incompatibility alleles in Arabidopsis halleri and Arabidopsis lyrata

The long-term balancing selection acting on mating types or sex-determining genes is expected to lead to the accumulation of deleterious mutations in the tightly linked chromosomal segments that are locally 'sheltered' from purifying selection. However, the factors determining the extent of this accumulation are poorly understood. Here, we took advantage of variations in the intensity of balancing selection along a dominance hierarchy formed by alleles at the sporophytic self-incompatibility system of the Brassicaceae to compare the pace at which linked deleterious mutations accumulate among them. We first experimentally measured the phenotypic manifestation of the linked load at three different levels of the dominance hierarchy. We then sequenced and phased polymorphisms in the chromosomal regions linked to 126 distinct copies of S-alleles in two populations of Arabidopsis halleri and three populations of Arabidopsis lyrata. We find that linkage to the S-locus locally distorts phylogenies over about 10-30 kb along the chromosome. The more intense balancing selection on dominant S-alleles results in greater fixation of linked deleterious mutations, while recessive S-alleles accumulate more linked deleterious mutations that are segregating. Hence, the structure rather than the overall magnitude of the linked genetic load differs between dominant and recessive S-alleles. Our results have consequences for the long-term evolution of new S-alleles, the evolution of dominance modifiers between them, and raise the question of why the non-recombining regions of some sex and mating type chromosomes expand over evolutionary times while others, such as the S-locus of the Brassicaceae, remain restricted to small chromosomal regions.

Stepwise recombination suppression around the mating-type locus in an ascomycete fungus with self-fertile spores

Recombination is often suppressed at sex-determining loci in plants and animals, and at self-incompatibility or mating-type loci in plants and fungi. In fungal ascomycetes, recombination suppression around the mating-type locus is associated with pseudo-homothallism, i.e. the production of self-fertile dikaryotic sexual spores carrying the two opposite mating types. This has been well studied in two species complexes from different families of Sordariales: Podospora anserina and Neurospora tetrasperma. However, it is unclear whether this intriguing association holds in other species. We show here that Schizothecium tetrasporum, a fungus from a third family in the order Sordariales, also produces mostly self-fertile dikaryotic spores carrying the two opposite mating types. This was due to a high frequency of second meiotic division segregation at the mating-type locus, indicating the occurrence of a single and systematic crossing-over event between the mating-type locus and the centromere, as in P. anserina. The mating-type locus has the typical Sordariales organization, plus a MAT1-1-1 pseudogene in the MAT1-2 haplotype. High-quality genome assemblies of opposite mating types and segregation analyses revealed a suppression of recombination in a region of 1.47 Mb around the mating-type locus. We detected three evolutionary strata, indicating a stepwise extension of recombination suppression. The three strata displayed no rearrangement or transposable element accumulation but gene losses and gene disruptions were present, and precisely at the strata margins. Our findings indicate a convergent evolution of self-fertile dikaryotic sexual spores across multiple ascomycete fungi. The particular pattern of meiotic segregation at the mating-type locus was associated with recombination suppression around this locus, that had extended stepwise. This association between pseudo-homothallism and recombination suppression across lineages and the presence of gene disruption at the strata limits are consistent with a recently proposed mechanism of sheltering deleterious alleles to explain stepwise recombination suppression.

Epistatic selection on a selfish Segregation Distorter supergene - drive, recombination, and genetic load

Meiotic drive supergenes are complexes of alleles at linked loci that together subvert Mendelian segregation resulting in preferential transmission. In males, the most common mechanism of drive involves the disruption of sperm bearing one of a pair of alternative alleles. While at least two loci are important for male drive-the driver and the target-linked modifiers can enhance drive, creating selection pressure to suppress recombination. In this work, we investigate the evolution and genomic consequences of an autosomal, multilocus, male meiotic drive system, Segregation Distorter (SD) in the fruit fly, Drosophila melanogaster. In African populations, the predominant SD chromosome variant, SD-Mal, is characterized by two overlapping, paracentric inversions on chromosome arm 2R and nearly perfect (~100%) transmission. We study the SD-Mal system in detail, exploring its components, chromosomal structure, and evolutionary history. Our findings reveal a recent chromosome-scale selective sweep mediated by strong epistatic selection for haplotypes carrying Sd, the main driving allele, and one or more factors within the double inversion. While most SD-Mal chromosomes are homozygous lethal, SD-Mal haplotypes can recombine with other, complementing haplotypes via crossing over, and with wildtype chromosomes via gene conversion. SD-Mal chromosomes have nevertheless accumulated lethal mutations, excess non-synonymous mutations, and excess transposable element insertions. Therefore, SD-Mal haplotypes evolve as a small, semi-isolated subpopulation with a history of strong selection. These results may explain the evolutionary turnover of SD haplotypes in different populations around the world and have implications for supergene evolution broadly.

Ecological Load and Balancing Selection in Circumboreal Barnacles

Acorn barnacle adults experience environmental heterogeneity at various spatial scales of their circumboreal habitat, raising the question of how adaptation to high environmental variability is maintained in the face of strong juvenile dispersal and mortality. Here, we show that 4% of genes in the barnacle genome experience balancing selection across the entire range of the species. Many of these genes harbor mutations maintained across 2 My of evolution between the Pacific and Atlantic oceans. These genes are involved in ion regulation, pain reception, and heat tolerance, functions which are essential in highly variable ecosystems. The data also reveal complex population structure within and between basins, driven by the trans-Arctic interchange and the last glaciation. Divergence between Atlantic and Pacific populations is high, foreshadowing the onset of allopatric speciation, and suggesting that balancing selection is strong enough to maintain functional variation for millions of years in the face of complex demography.

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.

Recombination suppression and evolutionary strata around mating-type loci in fungi: documenting patterns and understanding evolutionary and mechanistic causes

Genomic regions determining sexual compatibility often display recombination suppression, as occurs in sex chromosomes, plant self-incompatibility loci and fungal mating-type loci. Regions lacking recombination can extend beyond the genes determining sexes or mating types, by several successive steps of recombination suppression. Here we review the evidence for recombination suppression around mating-type loci in fungi, sometimes encompassing vast regions of the mating-type chromosomes. The suppression of recombination at mating-type loci in fungi has long been recognized and maintains the multiallelic combinations required for correct compatibility determination. We review more recent evidence for expansions of recombination suppression beyond mating-type genes in fungi ('evolutionary strata'), which have been little studied and may be more pervasive than commonly thought. We discuss testable hypotheses for the ultimate (evolutionary) and proximate (mechanistic) causes for such expansions of recombination suppression, including (1) antagonistic selection, (2) association of additional functions to mating-type, such as uniparental mitochondria inheritance, (3) accumulation in the margin of nonrecombining regions of various factors, including deleterious mutations or transposable elements resulting from relaxed selection, or neutral rearrangements resulting from genetic drift. The study of recombination suppression in fungi could thus contribute to our understanding of recombination suppression expansion across a broader range of organisms.

Evolution of the mating type gene pair and multiple sexes in Tetrahymena

The multiple mating type system of the Ciliate Tetrahymena thermophila is a self/non-self recognition system, whose specificity resides in a head-to-head, functionally distinct pair of genes, MTA and MTB. We have now sequenced and analyzed these mating type genes in nine additional Tetrahymena species. We conclude that MTA and MTB are derived from a common ancestral gene and have co-evolved for at least ∼150 Myr. We show that T. shanghaiensis, a perpetual selfer (unisexual) species, has a single mating type gene pair, whose MTA and MTB genes likely have different mating type specificity. We document the recent replacement of a complete different set of mating type specificities for another, illustrating how quickly this can happen. We discuss how varying conditions of reproductive stress could result in evolutionary co-adaptations of MTA and MTB genes and changes in mating type determination mechanisms.

•

The two Tetrahymena mating type proteins evolved from a common ancestor

•

Successive replacement waves generated the current diversity of mating type proteins

•

Well defined segments of both mating type proteins show differential lineage sorting

•

Perpetual selfer, T. shanghaiensis, has a heterospecific mating type gene pair

Footprints of natural selection at the mannose-6-phosphate isomerase locus in barnacles

The rocky intertidal is a natural laboratory to study how natural selection acts on the genes and proteins responsible for organismal survival and reproduction. Alternative forms of enzymes that differ across the intertidal have been known for decades and have provided examples of selection, but the genetic basis of such enzyme variation is known in only a few cases. In this paper, we present molecular evidence of natural selection at the Mpi gene, a key enzyme in energy metabolism that alters survival of barnacles living across the stress gradient imposed by the intertidal. Our study demonstrates how natural selection can facilitate survival in highly heterogeneous environments through the maintenance of multiple molecular solutions to ecological stresses.

The mannose-6-phosphate isomerase (Mpi) locus in Semibalanus balanoides has been studied as a candidate gene for balancing selection for more than two decades. Previous work has shown that Mpi allozyme genotypes (fast and slow) have different frequencies across Atlantic intertidal zones due to selection on postsettlement survival (i.e., allele zonation). We present the complete gene sequence of the Mpi locus and quantify nucleotide polymorphism in S. balanoides, as well as divergence to its sister taxon Semibalanus cariosus. We show that the slow allozyme contains a derived charge-altering amino acid polymorphism, and both allozyme classes correspond to two haplogroups with multiple internal haplotypes. The locus shows several footprints of balancing selection around the fast/slow site: an enrichment of positive Tajima’s D for nonsynonymous mutations, an excess of polymorphism, and a spike in the levels of silent polymorphism relative to silent divergence, as well as a site frequency spectrum enriched for midfrequency mutations. We observe other departures from neutrality across the locus in both coding and noncoding regions. These include a nonsynonymous trans-species polymorphism and a recent mutation under selection within the fast haplogroup. The latter suggests ongoing allelic replacement of functionally relevant amino acid variants. Moreover, predicted models of Mpi protein structure provide insight into the functional significance of the putatively selected amino acid polymorphisms. While footprints of selection are widespread across the range of S. balanoides, our data show that intertidal zonation patterns are variable across both spatial and temporal scales. These data provide further evidence for heterogeneous selection on Mpi.

Fitness differences suppress the number of mating types in evolving isogamous species

Sexual reproduction is not always synonymous with the existence of two morphologically different sexes; isogamous species produce sex cells of equal size, typically falling into multiple distinct self-incompatible classes, termed mating types. A long-standing open question in evolutionary biology is: what governs the number of these mating types across species? Simple theoretical arguments imply an advantage to rare types, suggesting the number of types should grow consistently; however, empirical observations are very different. While some isogamous species exhibit thousands of mating types, such species are exceedingly rare, and most have fewer than 10. In this paper, we present a mathematical analysis to quantify the role of fitness variation-characterized by different mortality rates-in determining the number of mating types emerging in simple evolutionary models. We predict that the number of mating types decreases as the variance of mortality increases.

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.

Little Evidence of Antagonistic Selection in the Evolutionary Strata of Fungal Mating-Type Chromosomes (Microbotryum lychnidis-dioicae)

Recombination suppression on sex chromosomes often extends in a stepwise manner, generating evolutionary strata of differentiation between sex chromosomes. Sexual antagonism is a widely accepted explanation for evolutionary strata, postulating that sets of genes beneficial in only one sex are successively linked to the sex-determining locus. The anther-smut fungus Microbotryum lychnidis-dioicae has mating-type chromosomes with evolutionary strata, only some of which link mating-type genes. Male and female roles are non-existent in this fungus, but mating-type antagonistic selection can also generate evolutionary strata, although the life cycle of the fungus suggests it should be restricted to few traits. Here, we tested the hypothesis that mating-type antagonism may have triggered recombination suppression beyond mating-type genes in M. lychnidis-dioicae by searching for footprints of antagonistic selection in evolutionary strata not linking mating-type loci. We found that these evolutionary strata (i) were not enriched in genes upregulated in the haploid phase, where cells are of alternative mating types, (ii) carried no gene differentially expressed between mating types, and (iii) carried no genes displaying footprints of specialization in terms of protein sequences (d_N/d_S) between mating types after recommended filtering. Without filtering, eleven genes showed signs of positive selection in the strata not linking mating-type genes, which constituted an enrichment compared to autosomes, but their functions were not obviously involved in antagonistic selection. Thus, we found no strong evidence that antagonistic selection has contributed to extending recombination suppression beyond mating-type genes. Alternative hypotheses should therefore be explored to improve our understanding of the sex-related chromosome evolution.

Rapid loss of self-incompatibility in experimental populations of the perennial outcrossing plant Linaria cavanillesii

Transitions from self-incompatibility to self-compatibility in angiosperms may be frequently driven by selection for reproductive assurance when mates or pollinators are rare, and are often succeeded by loss of inbreeding depression by purging. Here, we use experimental evolution to investigate the spread of self-compatibility from one such population of the perennial plant Linaria cavanillesii into self-incompatible (SI) populations that still have high inbreeding depression. We introduced self-compatible (SC) individuals at different frequencies into replicate experimental populations of L. cavanillesii that varied in access to pollinators. Our experiment revealed a rapid shift to self-compatibility in all replicates, driven by both greater seed set and greater outcross siring success of SC individuals. We discuss our results in the light of computer simulations that confirm the tendency of self-compatibility to spread into SI populations under the observed conditions. Our study illustrates the ease with which self-compatibility can spread among populations, a requisite for species-wide transitions from self-incompatibility to self-compatibility.

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.

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

Breaking linkage between mating compatibility factors: Tetrapolarity in Microbotryum

Linkage of genes determining separate self-incompatibility mechanisms is a general expectation of sexual eukaryotes that helps to resolve conflicts between reproductive assurance and recombination. However, in some organisms, multiple loci are required to be heterozygous in offspring while segregating independently in meiosis. This condition, termed "tetrapolarity" in basidiomycete fungi, originated in the ancestor to that phylum, and there have been multiple reports of subsequent transitions to "bipolarity" (i.e., linkage of separate mating factors). In the genus Microbotryum, we present the first report of the breaking of linkage between two haploid self-incompatibility factors and derivation of a tetrapolar breeding system. This breaking of linkage is associated with major alteration of genome structure, with the compatibility factors residing on separate mating-type chromosome pairs, reduced in size but retaining the structural dimorphism characteristic for regions of recombination suppression. The challenge to reproductive assurance from unlinked compatibility factors may be overcome by the automictic mating system in Microbotryum (i.e., mating among products of the same meiosis). As a curious outcome, this linkage transition and its effects upon outcrossing compatibility rates may reinforce automixis as a mating system. These observations contribute to understanding mating systems and linkage as fundamental principles of sexual life cycles, with potential impacts on conventional wisdom regarding mating-type evolution.

Convergent evolution at the gametophytic self-incompatibility system in Malus and Prunus

S-RNase-based gametophytic self-incompatibility (GSI) has evolved once before the split of the Asteridae and Rosidae. This conclusion is based on the phylogenetic history of the S-RNase that determines pistil specificity. In Rosaceae, molecular characterizations of Prunus species, and species from the tribe Pyreae (i.e., Malus, Pyrus, Sorbus) revealed different numbers of genes determining S-pollen specificity. In Prunus only one pistil and pollen gene determine GSI, while in Pyreae there is one pistil but multiple pollen genes, implying different specificity recognition mechanisms. It is thus conceivable that within Rosaceae the genes involved in GSI in the two lineages are not orthologous but possibly paralogous. To address this hypothesis we characterised the S-RNase lineage and S-pollen lineage genes present in the genomes of five Rosaceae species from three genera: M. × domestica (apple, self-incompatible (SI); tribe Pyreae), P. persica (peach, self-compatible (SC); Amygdaleae), P. mume (mei, SI; Amygdaleae), Fragaria vesca (strawberry, SC; Potentilleae), and F. nipponica (mori-ichigo, SI; Potentilleae). Phylogenetic analyses revealed that the Malus and Prunus S-RNase and S-pollen genes belong to distinct gene lineages, and that only Prunus S-RNase and SFB-lineage genes are present in Fragaria. Thus, S-RNase based GSI system of Malus evolved independently from the ancestral system of Rosaceae. Using expression patterns based on RNA-seq data, the ancestral S-RNase lineage gene is inferred to be expressed in pistils only, while the ancestral S-pollen lineage gene is inferred to be expressed in tissues other than pollen.

Exogenous selection rather than cytonuclear incompatibilities shapes asymmetrical fitness of reciprocal Arabidopsis hybrids

Reciprocal crosses between species often display an asymmetry in the fitness of F1 hybrids. This pattern, referred to as isolation asymmetry or Darwin's corollary to Haldane's rule, is a general feature of reproductive isolation in plants, yet factors determining its magnitude and direction remain unclear. We evaluated reciprocal species crosses between two naturally hybridizing diploid species of Arabidopsis to assess the degree of isolation asymmetry at different postmating life stages. We found that pollen from Arabidopsis arenosa will usually fertilize ovules from Arabidopsis lyrata; the reverse receptivity being less complete. Maternal A. lyrata parents set more F1 hybrid seed, but germinate at lower frequency, reversing the asymmetry. As predicted by theory, A. lyrata (the maternal parent with lower seed viability in crosses) exhibited accelerated chloroplast evolution, indicating that cytonuclear incompatibilities may play a role in reproductive isolation. However, this direction of asymmetrical reproductive isolation is not replicated in natural suture zones, where delayed hybrid breakdown of fertility at later developmental stages, or later-acting selection against A. arenosa maternal hybrids (unrelated to hybrid fertility, e.g., substrate adaptation) may be responsible for an excess of A. lyrata maternal hybrids. Exogenous selection rather than cytonuclear incompatibilities thus shapes the asymmetrical postmating isolation in nature.

The draft genome of Primula veris yields insights into the molecular basis of heterostyly

Background

The flowering plant Primula veris is a common spring blooming perennial that is widely cultivated throughout Europe. This species is an established model system in the study of the genetics, evolution, and ecology of heterostylous floral polymorphisms. Despite the long history of research focused on this and related species, the continued development of this system has been restricted due the absence of genomic and transcriptomic resources.

Results

We present here a de novo draft genome assembly of P. veris covering 301.8 Mb, or approximately 63% of the estimated 479.22 Mb genome, with an N50 contig size of 9.5 Kb, an N50 scaffold size of 164 Kb, and containing an estimated 19,507 genes. The results of a RADseq bulk segregant analysis allow for the confident identification of four genome scaffolds that are linked to the P. veris S-locus. RNAseq data from both P. veris and the closely related species P. vulgaris allow for the characterization of 113 candidate heterostyly genes that show significant floral morph-specific differential expression. One candidate gene of particular interest is a duplicated GLOBOSA homolog that may be unique to Primula (PveGLO2), and is completely silenced in L-morph flowers.

Conclusions

The P. veris genome represents the first genome assembled from a heterostylous species, and thus provides an immensely important resource for future studies focused on the evolution and genetic dissection of heterostyly. As the first genome assembled from the Primulaceae, the P. veris genome will also facilitate the expanded application of phylogenomic methods in this diverse family and the eudicots as a whole.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0567-z) contains supplementary material, which is available to authorized users.

Sustained heterozygosity across a self-incompatibility locus in an inbred ascidian

Because self-incompatibility loci are maintained heterozygous and recombination within self-incompatibility loci would be disadvantageous, self-incompatibility loci are thought to contribute to structural and functional differentiation of chromosomes. Although the hermaphrodite chordate, Ciona intestinalis, has two self-incompatibility genes, this incompatibility system is incomplete and self-fertilization occurs under laboratory conditions. Here, we established an inbred strain of C. intestinalis by repeated self-fertilization. Decoding genome sequences of sibling animals of this strain identified a 2.4-Mbheterozygous region on chromosome 7. A self-incompatibility gene, Themis-B, was encoded within this region. This observation implied that this self-incompatibility locus and the linkage disequilibrium of its flanking region contribute to the formation of the 2.4-Mb heterozygous region, probably through recombination suppression. We showed that different individuals in natural populations had different numbers and different combinations of Themis-B variants, and that the rate of self-fertilization varied among these animals. Our result explains why self-fertilization occurs under laboratory conditions. It also supports the concept that the Themis-B locus is preferentially retained heterozygous in the inbred line and contributes to the formation of the 2.4-Mb heterozygous region. High structural variations might suppress recombination, and this long heterozygous region might represent a preliminary stage of structural differentiation of chromosomes.

Species and population level molecular profiling reveals cryptic recombination and emergent asymmetry in the dimorphic mating locus of C. reinhardtii

Heteromorphic sex-determining regions or mating-type loci can contain large regions of non-recombining sequence where selection operates under different constraints than in freely recombining autosomal regions. Detailed studies of these non-recombining regions can provide insights into how genes are gained and lost, and how genetic isolation is maintained between mating haplotypes or sex chromosomes. The Chlamydomonas reinhardtii mating-type locus (MT) is a complex polygenic region characterized by sequence rearrangements and suppressed recombination between its two haplotypes, MT+ and MT−. We used new sequence information to redefine the genetic contents of MT and found repeated translocations from autosomes as well as sexually controlled expression patterns for several newly identified genes. We examined sequence diversity of MT genes from wild isolates of C. reinhardtii to investigate the impacts of recombination suppression. Our population data revealed two previously unreported types of genetic exchange in Chlamydomonas MT—gene conversion in the rearranged domains, and crossover exchanges in flanking domains—both of which contribute to maintenance of genetic homogeneity between haplotypes. To investigate the cause of blocked recombination in MT we assessed recombination rates in crosses where the parents were homozygous at MT. While normal recombination was restored in MT+×MT+ crosses, it was still suppressed in MT−×MT− crosses. These data revealed an underlying asymmetry in the two MT haplotypes and suggest that sequence rearrangements are insufficient to fully account for recombination suppression. Together our findings reveal new evolutionary dynamics for mating loci and have implications for the evolution of heteromorphic sex chromosomes and other non-recombining genomic regions.

Sex chromosomes and mating-type loci are often atypical in their structure and evolutionary dynamics. One distinguishing feature is the absence of recombination that results in genetic isolation and promotes rapid evolution and sometimes degeneration. We investigated gene content, sex-regulated expression, and recombination of mating locus (MT) genes in the unicellular alga Chlamydomonas reinhardtii. Despite the lack of observable recombination in and around Chlamydomonas MT, genes from its two mating types are far more similar to each other than expected for a non-recombining region. This discrepancy is explained by our finding evidence of genetic exchange between the two mating types within wild populations. In addition, we observed an unexpected asymmetry in the recombination behavior of the two mating types that may have contributed to the preferential expansion of one MT haplotype over the other through insertion of new genes. Our data suggest a mechanism to explain the emergence of heteromorphic sex chromosomes in haploid organisms by asymmetric expansion rather than by loss or degeneration as occurs in some Y or W chromosomes from diploid organisms. Our observations support a revised view of recombination in sex-determining regions as a quantitative phenomenon that can significantly affect rates of evolution and sex-linked genetic diversification.

High-density genetic maps for loci involved in nuclear male sterility (NMS1) and sporophytic self-incompatibility (S-locus) in chicory (Cichorium intybus L., Asteraceae)

High-density genetic maps were constructed for loci involved in nuclear male sterility (NMS1-locus) and sporophytic self-incompatibility (S-locus) in chicory (Cichorium intybus L.). The mapping population consisted of 389 F1' individuals derived from a cross between two plants, K28 (male-sterile) and K59 (pollen-fertile), both heterozygous at the S-locus. This F1' mapping population segregated for both male sterility (MS) and strong self-incompatibility (SI) phenotypes. Phenotyping F1' individuals for MS allowed us to map the NMS1-locus to linkage group (LG) 5, while controlled diallel and factorial crosses to identify compatible/incompatible phenotypes mapped the S-locus to LG2. To increase the density of markers around these loci, bulked segregant analysis was used. Bulks and parental plants K28 and K59 were screened using amplified fragment length polymorphism (AFLP) analysis, with a complete set of 256 primer combinations of EcoRI-ANN and MseI-CNN. A total of 31,000 fragments were generated, of which 2,350 showed polymorphism between K59 and K28. Thirteen AFLP markers were identified close to the NMS1-locus and six in the vicinity of the S-locus. From these AFLP markers, eight were transformed into sequence-characterized amplified region (SCAR) markers and of these five showed co-dominant polymorphism. The chromosomal regions containing the NMS1-locus and the S-locus were each confined to a region of 0.8 cM. In addition, we mapped genes encoding proteins similar to S-receptor kinase, the female determinant of sporophytic SI in the Brasicaceae, and also markers in the vicinity of the putative S-locus of sunflower, but none of these genes or markers mapped close to the chicory S-locus.

Extensive trans-specific polymorphism at the mating type locus of the root decay fungus Heterobasidion

Incompatibility systems in which individuals bearing identical alleles reject each other favor the maintenance of a diversity of alleles. Mushroom mating type loci (MAT) encode for dozens or hundreds of incompatibility alleles whose loss from the population is greatly restricted through negative frequency selection, leading to a system of alleles with highly divergent sequences. Here, we use DNA sequences of homeodomain (HD) encoding genes at the MAT locus of five closely related species of the root rot basidiomycete Heterobasidion annosum sensu lato to show that the extended coalescence time of MAT alleles greatly predates speciation in the group, contrasting loci outside of MAT that show allele divergences largely consistent with the species phylogeny with those of MAT, which show rampant trans-species polymorphism. We observe a roughly 6-fold greater genealogical depth and polymorphism of MAT compared with non-MAT that argues for the maintenance of balanced polymorphism for a minimum duration of 24 My based on a molecular-clock calibrated species phylogeny. As with other basidiomycete HD genes, balancing selection appears to be concentrated at the specificity-determining region in the N-terminus of the protein based on identification of codons under selection and the absence of recombination within the region. However, the elevated polymorphism extends into the nonspecificity determining regions as well as a neighboring non-MAT gene, the mitochondrial intermediate peptidase (MIP). In doing so, increased divergence should decrease recombination among alleles and as a by-product create incompatibilities in the functional domains not involved in allele recognition but in regulating sexual development.

Mutualism and asexual reproduction influence recognition genes in a fungal symbiont

Mutualism between microbes and insects is common and alignment of the reproductive interests of microbial symbionts with this lifestyle typically involves clonal reproduction and vertical transmission by insect partners. Here the Amylostereum fungus-Sirex woodwasp mutualism was used to consider whether their prolonged association and predominance of asexuality have affected the mating system of the fungal partner. Nucleotide information for the pheromone receptor gene rab1, as well as the translation elongation factor 1α gene and ribosomal RNA internal transcribed spacer region were utilized. The identification of rab1 alleles in Amylostereum chailletii and Amylostereum areolatum populations revealed that this gene is more polymorphic than the other two regions, although the diversity of all three regions was lower than what has been observed in free-living Agaricomycetes. Our data suggest that suppressed recombination might be implicated in the diversification of rab1, while no evidence of balancing selection was detected. We also detected positive selection at only two codons, suggesting that purifying selection is important for the evolution of rab1. The symbiotic relationship with their insect partners has therefore influenced the diversity of this gene and influenced the manner in which selection drives and maintains this diversity in A. areolatum and A. chailletii.

Inbreeding depression in urban environments of the bird's nest fungus Cyathus stercoreus (Nidulariaceae: Basidiomycota)

Many organisms display codispersal of offspring, but fewer display codispersal of compatible gametes. This mechanism enhances the ability of a species to colonize after long distance dispersal as a mechanism of reproductive assurance, but it also fosters inbreeding and potential reduction in fitness. Here we investigated both long distance dispersal and inbreeding in the bird's nest fungus Cyathus stercoreus, a dung and mulch-associated fungus with a splash cup fruiting body appearing like a miniature bird's nest of 'eggs' or peridioles that contain thousands of mating compatible meiotic spores. To investigate the genetic structure in the species, six North American urban populations were hierarchically sampled and genotyped using 10 microsatellite markers. We detected significant levels of inbreeding through heterozygote deficiencies at four loci, with global FIS=0.061. Dispersal limitation was suggested by both spatial autocorrelation and the detection of population structure between Louisiana and Michigan using clustering and F-statistics. Although inbreeding may facilitate colonization by the fungus, it has a negative effect on the fitness of populations as estimated from a 15% reduction in growth rates of inbred strains relative to outcrossed. Mating tests revealed that C. stercoreus has a higher estimated number of mating-type alleles (MAT-A= 39, MAT-B= 24) than other species of bird's nest fungi, which would increase its outcrossing efficiency. We speculate that the increased number of mating-type alleles is the result of a recent range and population size expansion into urban environments.

Inbreeding depression in self-incompatible North-American Arabidopsis lyrata: disentangling genomic and S-locus-specific genetic load

Newly formed selfing lineages may express recessive genetic load and suffer inbreeding depression. This can have a genome-wide genetic basis, or be due to loci linked to genes under balancing selection. Understanding the genetic architecture of inbreeding depression is important in the context of the maintenance of self-incompatibility and understanding the evolutionary dynamics of S-alleles. We addressed this using North-American subspecies of Arabidopsis lyrata. This species is normally self-incompatible and outcrossing, but some populations have undergone a transition to selfing. The goals of this study were to: (1) quantify the strength of inbreeding depression in North-American populations of A. lyrata; and (2) disentangle the relative contribution of S-linked genetic load compared with overall inbreeding depression. We enforced selfing in self-incompatible plants with known S-locus genotype by treatment with CO(2), and compared the performance of selfed vs outcrossed progeny. We found significant inbreeding depression for germination rate (δ=0.33), survival rate to 4 weeks (δ=0.45) and early growth (δ=0.07), but not for flowering rate. For two out of four S-alleles in our design, we detected significant S-linked load reflected by an under-representation of S-locus homozygotes in selfed progeny. The presence or absence of S-linked load could not be explained by the dominance level of S-alleles. Instead, the random nature of the mutation process may explain differences in the recessive deleterious load among lineages.

Microevolution of S-allele frequencies in wild cherry populations: respective impacts of negative frequency dependent selection and genetic drift

Negative frequency dependent selection (NFDS) is supposed to be the main force controlling allele evolution at the gametophytic self-incompatibility locus (S-locus) in strictly outcrossing species. Genetic drift also influences S-allele evolution. In perennial sessile organisms, evolution of allelic frequencies over two generations is mainly shaped by individual fecundities and spatial processes. Using wild cherry populations between two successive generations, we tested whether S-alleles evolved following NFDS qualitative and quantitative predictions. We showed that allelic variation was negatively correlated with parental allelic frequency as expected under NFDS. However, NFDS predictions in finite population failed to predict more than half S-allele quantitative evolution. We developed a spatially explicit mating model that included the S-locus. We studied the effects of self-incompatibility and local drift within populations due to pollen dispersal in spatially distributed individuals, and variation in male fecundity on male mating success and allelic frequency evolution. Male mating success was negatively related to male allelic frequency as expected under NFDS. Spatial genetic structure combined with self-incompatibility resulted in higher effective pollen dispersal. Limited pollen dispersal in structured distributions of individuals and genotypes and unequal pollen production significantly contributed to S-allele frequency evolution by creating local drift effects strong enough to counteract the NFDS effect on some alleles.

Evolution of the S-locus region in Arabidopsis relatives

The S locus, a single polymorphic locus, is responsible for self-incompatibility (SI) in the Brassicaceae family and many related plant families. Despite its importance, our knowledge of S-locus evolution is largely restricted to the causal genes encoding the S-locus receptor kinase (SRK) receptor and S-locus cysteine-rich protein (SCR) ligand of the SI system. Here, we present high-quality sequences of the genomic region of six S-locus haplotypes: Arabidopsis (Arabidopsis thaliana; one haplotype), Arabidopsis lyrata (four haplotypes), and Capsella rubella (one haplotype). We compared these with reference S-locus haplotypes of the self-compatible Arabidopsis and its SI congener A. lyrata. We subsequently reconstructed the likely genomic organization of the S locus in the most recent common ancestor of Arabidopsis and Capsella. As previously reported, the two SI-determining genes, SCR and SRK, showed a pattern of coevolution. In addition, consistent with previous studies, we found that duplication, gene conversion, and positive selection have been important factors in the evolution of these two genes and appear to contribute to the generation of new recognition specificities. Intriguingly, the inactive pseudo-S-locus haplotype in the self-compatible species C. rubella is likely to be an old S-locus haplotype that only very recently became fixed when C. rubella split off from its SI ancestor, Capsella grandiflora.

DNA polymorphism in recombining and non-recombining mating-type-specific loci of the smut fungus Microbotryum

The population-genetic processes leading to the genetic degeneration of non-recombining regions have mainly been studied in animal and plant sex chromosomes. Here, we report population genetic analysis of the processes in the non-recombining mating-type-specific regions of the smut fungus Microbotryum violaceum. M. violaceum has A1 and A2 mating types, determined by mating-type-specific 'sex chromosomes' that contain 1-2 Mb long non-recombining regions. If genetic degeneration were occurring, then one would expect reduced DNA polymorphism in the non-recombining regions of this fungus. The analysis of DNA diversity among 19 M. violaceum strains, collected across Europe from Silene latifolia flowers, revealed that (i) DNA polymorphism is relatively low in all 20 studied loci (π∼0.15%), (ii) it is not significantly different between the two mating-type-specific chromosomes nor between the non-recombining and recombining regions, (iii) there is substantial population structure in M. violaceum populations, which resembles that of its host species, S. latifolia, and (iv) there is significant linkage disequilibrium, suggesting that widespread selfing in this species results in a reduction of the effective recombination rate across the genome. We hypothesise that selfing-related reduction of recombination across the M. violaceum genome negates the difference in the level of DNA polymorphism between the recombining and non-recombining regions, and may possibly lead to similar levels of genetic degeneration in the mating-type-specific regions of the non-recombining 'sex chromosomes' and elsewhere in the genome.

High variability in a mating type linked region in the dry rot fungus Serpula lacrymans caused by frequency-dependent selection?

BACKGROUND

The mating type loci that govern the mating process in fungi are thought to be influenced by negative frequency-dependent selection due to rare allele advantage. In this study we used a mating type linked DNA marker as a proxy to indirectly study the allelic richness and geographic distribution of mating types of one mating type locus (MAT A) in worldwide populations of the dry rot fungus Serpula lacrymans. This fungus, which causes serious destruction to wooden constructions in temperate regions worldwide, has recently expanded its geographic range with a concomitant genetic bottleneck.

RESULTS

High allelic richness and molecular variation was detected in the mating type linked marker as compared to other presumably neutral markers. Comparable amounts of genetic variation appeared in the mating type linked marker in populations from nature and buildings, which contrast the pattern observed with neutral genetic markers where natural populations were far more variable. Some geographic structuring of the allelic variation in the mating type linked marker appeared, but far less than that observed with neutral markers. In founder populations of S. lacrymans, alleles co-occurring in heterokaryotic individuals were more divergent than expected by chance, which agrees with the expectation for populations where few mating alleles exists. The analyzed DNA marker displays trans-species polymorphism wherein some alleles from the closely related species S. himantoides are more similar to those of S. lacrymans than other alleles from S. himantoides.

CONCLUSIONS

Our results support the idea that strong negative frequency-dependent selection maintains high levels of genetic variation in MAT-linked genomic regions, even in recently bottlenecked populations of S. lacrymans.

The sheltered genetic load linked to the s locus in plants: new insights from theoretical and empirical approaches in sporophytic self-incompatibility

Inbreeding depression and mating systems evolution are closely linked, because the purging of deleterious mutations and the fitness of individuals may depend on outcrossing vs. selfing rates. Further, the accumulation of deleterious mutations may vary among genomic regions, especially for genes closely linked to loci under balancing selection. Sporophytic self-incompatibility (SSI) is a common genetic mechanism in angiosperm that enables hermaphrodite plants to avoid selfing and promote outcrossing. The SSI phenotype is determined by the S locus and may depend on dominance relationships among alleles. Since most individuals are heterozygous at the S locus and recombination is suppressed in the S-locus region, it has been suggested that deleterious mutations could accumulate at genes linked to the S locus, generating a "sheltered load." In this article, we first theoretically investigate the conditions generating sheltered load in SSI. We show that deleterious mutations can accumulate in linkage with specific S alleles, and particularly if those S alleles are dominant. Second, we looked for the presence of sheltered load in Arabidopsis halleri using CO(2) gas treatment to overcome self-incompatibility. By examining the segregation of S alleles and measuring the relative fitness of progeny, we found significant sheltered load associated with the most dominant S allele (S15) of three S alleles tested. This sheltered load seems to be expressed at several stages of the life cycle and to have a larger effect than genomic inbreeding depression.

A new theory of MHC evolution: beyond selection on the immune genes

The major histocompatibility complex (MHC) is a dense region of immune genes with high levels of polymorphism, which are arranged in haplotype blocks. Traditional models of balancing selection (i.e. overdominance and negative frequency dependence) were developed to study the population genetics of single genes. However, the MHC is a multigene family surrounded by linked (non-neutral) polymorphisms, and not all of its features are well explained by these models. For example, (i) the high levels of polymorphism in small populations, (ii) the unexpectedly large genetic differentiation between populations, (iii) the shape of the allelic genealogy associated with trans-species evolution, and (iv) the close associations between particular MHC (human leucocyte antigen, HLA) haplotypes and the approximately 100 pathologies in humans. Here, I propose a new model of MHC evolution named Associative Balancing Complex evolution that can explain these phenomena. The model proposes that recessive deleterious mutations accumulate as a 'sheltered load' nearby MHC genes. These mutations can accumulate because (i) they are rarely expressed as homozygotes given the high MHC gene diversity and (ii) purifying selection is inefficient with low recombination rates (cf. Muller's ratchet). Once fixed, these mutations add to balancing selection and further reinforce linkage through epistatic selection against recombinants.

Genetic linkage map for Amylostereum areolatum reveals an association between vegetative growth and sexual and self-recognition

Amylostereum areolatum is a filamentous fungus that grows through tip extension, branching and hyphal fusion. In the homokaryotic phase, the hyphae of different individuals are capable of fusing followed by heterokaryon formation, only if they have dissimilar allelic specificities at their mating-type (mat) loci. In turn, hyphal fusion between heterokaryons persists only when they share the same alleles at all of their heterokaryon incompatibility (het) loci. In this study we present the first genetic linkage map for A. areolatum, onto which the mat and het loci, as well as quantitative trait loci (QTLs) for mycelial growth rate are mapped. The recognition loci (mat-A and het-A) are positioned near QTLs associated with mycelial growth, suggesting that the genetic determinants influencing recognition and growth rate in A. areolatum are closely associated. This was confirmed when isolates associated with specific mat and het loci displayed significantly different mycelial growth rates. Although the link between growth and sexual recognition has previously been observed in other fungi, this is the first time that an association between growth and self-recognition has been shown.

Ancient trans-specific polymorphism at pheromone receptor genes in basidiomycetes

In the majority of sexual organisms, reproduction occurs almost exclusively through the combination of distinct and alternate forms, called sexes or mating types. In some fungi, there can be dozens to hundreds of alternate alleles that determine compatible mating types. Such extensive polymorphism is expected to be maintained by balancing selection, and in extreme cases may give rise to trans-specific polymorphism. Here, we analyzed sequences of two pheromone receptors in the Microbotryum fungal species complex (Basidiomycota), which has only two alternate mating types. Several lines of evidence strongly suggest that the pheromone receptors are two allelic sequences acting to determine the alternate A1 and A2 mating types required for mating in Microbotryum. Phylogenetic trees of pheromone receptors in the Microbotryum species complex indicated a trans-specific polymorphism: the Microbotryum sequences from a given mating type were all more similar to the pheromone receptors of distantly related classes of fungi than to the alternate pheromone receptor in the Microbotryum species. A phylogenetic tree built using other known pheromone receptors from basidiomycetes showed that trans-specific polymorphism is widespread. The pheromone receptor alleles from Microbotryum appeared as the oldest, being at least 370 million years old. This represents the oldest known trans-specific polymorphism known in any organism so far, which may be due to the existence of sex chromosomes, obligate sexuality, mitochondrial inheritance linked to the mating type, and a highly selfing mating system in Microbotryum.

Recombination at Prunus S-locus region SLFL1 gene

In Prunus, the self-incompatibility (S-) locus region is <70 kb. Two genes--the S-RNase, which encodes the functional female recognition component, and the SFB gene, which encodes the pollen recognition component--must co-evolve as a genetic unit to maintain functional incompatibility. Therefore, recombination must be severely repressed at the S-locus. Levels of recombination at genes in the vicinity of the S-locus have not yet been rigorously tested; thus it is unknown whether recombination is also severely repressed at these loci. In this work, we looked at variability levels and patterns at the Prunus spinosa SLFL1 gene, which is physically close to the S-RNase gene. Our results suggest that the recombination level increases near the SLFL1 coding region. These findings are discussed in the context of theoretical models predicting an effect of linked weakly deleterious mutations on the relatedness of S-locus specificities. Moreover, we show that SLFL1 belongs to a gene family of at least five functional genes and that SLFL1 pseudogenes are frequently found in the S-locus region.

Preservation of a pseudogene by gene conversion and diversifying selection

Interlocus gene conversion is considered a crucial mechanism for generating novel combinations of polymorphisms in duplicated genes. The importance of gene conversion between duplicated genes has been recognized in the major histocompatibility complex and self-incompatibility genes, which are likely subject to diversifying selection. To theoretically understand the potential role of gene conversion in such situations, forward simulations are performed in various two-locus models. The results show that gene conversion could significantly increase the number of haplotypes when diversifying selection works on both loci. We find that the tract length of gene conversion is an important factor to determine the efficacy of gene conversion: shorter tract lengths can more effectively generate novel haplotypes given the gene conversion rate per site is the same. Similar results are also obtained when one of the duplicated genes is assumed to be a pseudogene. It is suggested that a duplicated gene, even after being silenced, will contribute to increasing the variability in the other locus through gene conversion. Consequently, the fixation probability and longevity of duplicated genes increase under the presence of gene conversion. On the basis of these findings, we propose a new scenario for the preservation of a duplicated gene: when the original donor gene is under diversifying selection, a duplicated copy can be preserved by gene conversion even after it is pseudogenized.

Unequal segregation of SRK alleles at the S locus in Brassica cretica

In the Brassicaceae plant family, which includes the Arabidopsis and Brassica genera, self-incompatibility (SI) is controlled by genes at the S locus. Using experimental crosses, we studied the pattern of inheritance of S-locus alleles in the wild species Brassica cretica. Four full-sib families were established and unequal segregation of alleles at the SRK SI gene was found in one family. The segregation distortion acted in favour of a recessive (class II) allele and was best explained by some form of gametic-level selection. Our findings are discussed in the light of theoretical predictions of differential accumulation of deleterious mutations among S-locus alleles.

Genomic consequences of selection on self-incompatibility genes

Frequency-dependent selection at plant self-incompatibility systems is inherent and well understood theoretically. A self-incompatibility locus leads to a strong peak of diversity in the genome, to a unique distribution of diversity across the species and possibly to increased introgression between closely related species. We review recent empirical studies demonstrating these features and relate the empirical findings to theoretical predictions. We show how these features are being exploited in searches for other genes under multi-allelic balancing selection and for inference on recent breakdown of self-incompatibility.

Mechanism of self-sterility in a hermaphroditic chordate

Hermaphroditic organisms avoid inbreeding by a system of self-incompatibility (SI). A primitive chordate (ascidian) Ciona intestinalis is an example of such an organism, but the molecular mechanism underlying its SI system is not known. Here, we show that the SI system is governed by two gene loci that act cooperatively. Each locus contains a tightly linked pair of polycystin 1-related receptor (s-Themis) and fibrinogen-like ligand (v-Themis) genes, the latter of which is located in the first intron of s-Themis but transcribed in the opposite direction. These genes may encode male- and female-side self-recognition molecules. The SI system of C. intestinalis has a similar framework to that of flowering plants but utilizing different molecules.

Expression of 10 S-class SLF-like genes in Nicotiana alata pollen and its implications for understanding the pollen factor of the S locus

The S locus of Nicotiana alata encodes a polymorphic series of ribonucleases (S-RNases) that determine the self-incompatibility (SI) phenotype of the style. The pollen product of the S locus (pollen S) in N. alata is unknown, but in species from the related genus Petunia and in self-incompatible members of the Plantaginaceae and Rosaceae, this function has been assigned to an F-box protein known as SLF or SFB. Here we describe the identification of 10 genes (designated DD1-10) encoding SLF-related proteins that are expressed in N. alata pollen. Because our approach to cloning the DD genes was based on sequences of SLFs from other species, we presume that one of the DD genes encodes the N. alata SLF ortholog. Seven of the DD genes were exclusively expressed in pollen and a low level of sequence variation was found in alleles of each DD gene. Mapping studies confirmed that all 10 DD genes were linked to the S locus and that at least three were located in the same chromosomal segment as pollen S. Finally, the different topologies of the phylogenetic trees produced using available SLF-related sequences and those produced using S-RNase sequences suggests that pollen S and the S-RNase have different evolutionary histories.

Effects of recombination on hitchhiking diversity in the Brassica self-incompatibility locus complex

In self-incompatibility, a number of S haplotypes are maintained by frequency-dependent selection, which results in trans-specific S haplotypes. The region of several kilobases (approximately 40-60 kb) from SP6 to SP2, including self-incompatibility-related genes and some adjacent genes in Brassica rapa, has high nucleotide diversity due to the hitchhiking effect, and therefore we call this region the "S-locus complex." Recombination in the S-locus complex is considered to be suppressed. We sequenced regions of >50 kb of the S-locus complex of three S haplotypes in B. rapa and found higher nucleotide diversity in intergenic regions than in coding regions. Two highly similar regions of >10 kb were found between BrS-8 and BrS-46. Phylogenetic analysis using trans-specific S haplotypes (called interspecific pairs) of B. rapa and B. oleracea suggested that recombination reduced the nucleotide diversity in these two regions and that the genes not involved in self-incompatibility in the S-locus complex and the kinase domain, but not the S domain, of SRK have also experienced recombination. Recombination may reduce hitchhiking diversity in the S-locus complex, whereas the region from the S domain to SP11 would disfavor recombination.

Sex determination in Chlamydomonas

The sex-determination system of the unicellular green alga, Chlamydomonas reinhardtii, is governed by genes in the mating-type (MT) locus and entails additional genes located in autosomes. Gene expression is initiated by nitrogen starvation, and cells differentiate into plus or minus gametes within 6h. Reviewed is our current understanding of gametic differentiation and fertilization, initiation of zygote development, and the uniparental inheritance of organelle genomes.

S locus genes and the evolution of self-fertility in Arabidopsis thaliana

Loss of self-incompatibility (SI) in Arabidopsis thaliana was accompanied by inactivation of genes required for SI, including S-LOCUS RECEPTOR KINASE (SRK) and S-LOCUS CYSTEINE-RICH PROTEIN (SCR), coadapted genes that constitute the SI specificity-determining S haplotype. Arabidopsis accessions are polymorphic for PsiSRK and PsiSCR, but it is unknown if the species harbors structurally different S haplotypes, either representing relics of ancestral functional and structurally heteromorphic S haplotypes or resulting from decay concomitant with or subsequent to the switch to self-fertility. We cloned and sequenced the S haplotype from C24, in which self-fertility is due solely to S locus inactivation, and show that this haplotype was produced by interhaplotypic recombination. The highly divergent organization and sequence of the C24 and Columbia-0 (Col-0) S haplotypes demonstrate that the A. thaliana S locus underwent extensive structural remodeling in conjunction with a relaxation of selective pressures that once preserved the integrity and linkage of coadapted SRK and SCR alleles. Additional evidence for this process was obtained by assaying 70 accessions for the presence of C24- or Col-0-specific sequences. Furthermore, analysis of SRK and SCR polymorphisms in these accessions argues against the occurrence of a selective sweep of a particular allele of SCR, as previously proposed.

The transition to self-compatibility in Arabidopsis thaliana and evolution within S-haplotypes over 10 Myr

A recent investigation found evidence that the transition of Arabidopsis thaliana from ancestral self-incompatibility (SI) to full self-compatibility occurred very recently and suggested that this occurred through a selective fixation of a nonfunctional allele (PsiSCR1) at the SCR gene, which determines pollen specificity in the incompatibility response. The main evidence is the lack of polymorphism at the SCR locus in A. thaliana. However, the nearby SRK gene, which determines stigma specificity in self-incompatible Brassicaceae species, has extremely high sequence diversity, with 3 very divergent SRK haplotypes, 2 of them present in multiple strains. Such high diversity is extremely unusual in this species, and it suggests the possibility that multiple, different SRK haplotypes may have been preserved from A. thaliana's self-incompatible ancestor. To study the evolution of S-haplotypes in the A. thaliana lineage, we searched the 2 most closely related Arabidopsis species Arabidopsis lyrata and Arabidopsis halleri, in which most populations have retained SI, and found SRK sequences corresponding to all 3 A. thaliana haplogroup sequences. Our molecular evolutionary analyses of these 3 S-haplotypes provide an independent estimate of the timing of the breakdown of SI and again exclude an ancient transition to selfing in A. thaliana. Comparing sequences of each of the 3 haplogroups between species, we find that 2 of the 3 SRK sequences (haplogroups A and B) are similar throughout their length, suggesting that little or no recombination with other SRK alleles has occurred since these species diverged. The diversity difference between the SCR and SRK loci in A. thaliana, however, suggests crossing-over, either within SRK or between the SCR and SRK loci. If the loss of SI involved fixation of the PsiSCR1 sequence, the exchange must have occurred during its fixation. Divergence between the species is much lower at the S-locus, compared with reference loci, and we discuss two contributory possibilities. Introgression may have occurred between A. lyrata and A. halleri and between their ancestral lineage and A. thaliana, at least for some period after their split. In addition, the coalescence times of sequences of individual S-haplogroups are expected to be less than those of alleles at non-S-loci.