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

Recent and Ancient Signature of Balancing Selection around the S-Locus in Arabidopsis halleri and A. lyrata

Balancing selection can maintain different alleles over long evolutionary times. Beyond this direct effect on the molecular targets of selection, balancing selection is also expected to increase neutral polymorphism in linked genome regions, in inverse proportion to their genetic map distances from the selected sites. The genes controlling plant self-incompatibility are subject to one of the strongest forms of balancing selection, and they show clear signatures of balancing selection. The genome region containing those genes (the S-locus) is generally described as nonrecombining, and the physical size of the region with low recombination has recently been established in a few species. However, the size of the region showing the indirect footprints of selection due to linkage to the S-locus is only roughly known. Here, we improved estimates of this region by surveying synonymous polymorphism and estimating recombination rates at 12 flanking region loci at known physical distances from the S-locus region boundary, in two closely related self-incompatible plants Arabidopsis halleri and A. lyrata. In addition to studying more loci than previous studies and using known physical distances, we simulated an explicit demographic scenario for the divergence between the two species, to evaluate the extent of the genomic region whose diversity departs significantly from neutral expectations. At the closest flanking loci, we detected signatures of both recent and ancient indirect effects of selection on the S-locus flanking genes, finding ancestral polymorphisms shared by both species, as well as an excess of derived mutations private to either species. However, these effects are detected only in a physically small region, suggesting that recombination in the flanking regions is sufficient to quickly break up linkage disequilibrium with the S-locus. Our approach may be useful for distinguishing cases of ancient versus recently evolved balancing selection in other systems

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

The effect of sheltered load on reproduction in Solanum carolinense, a species with variable self-incompatibility

In previous studies, we have investigated the strength of self-incompatibility (SI) in Solanum carolinense, a highly successful weed with a fully functional SI system that inhabits early successional and other disturbed habitats. We have found that the SI response in S. carolinense is a plastic trait—its strength being affected by the age of the flowers, and the presence of developing fruits and that there are genetic differences among families in their self-fertility. However, in species with a fully functional SI response, selfing would not be that common. As a result, deleterious recessives scattered though the genome of horsenettle are only occasionally exposed to selection. It has been suggested that deleterious recessives accumulate near S-alleles in strong SI species because the S-locus is located in a non-recombining region of the genome and because strong S-alleles are never in the homozygous state, thus sheltering some of the genetic load near the S-locus from selection. We performed a series of laboratory and greenhouse experiments to determine the extent to which sheltered load adds to the overall magnitude of inbreeding depression in horsenettle. Specifically, we amplified and sequenced the S-alleles from 16 genets collected from a large population in Pennsylvania and performed a series of controlled self-pollinations. We then grew the selfed progeny in the greenhouse; recorded various measures of growth and reproductive output; and amplified and sequenced their S-allele(s). We found that the heterozygous progeny of self-pollinations produce more flowers and have a greater ability to set both self and cross seed than S-homozygous progeny. We also found evidence of variation in the magnitude of load among S-alleles. These results suggest that sheltered load might slow the fixation of weak (partially compatible) S-alleles in this population, thus adding to the maintenance of a mixed mating system rather than leading to the fixation of the selfing allele

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

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

QTL analysis of heterostyly in Primula sieboldii and its application for morph identification in wild populations

Background and Aims Primula sieboldii is a perennial clonal herb that is distributed around the Sea of Japan and is endangered in Japan. Its breeding system is characterized by heteromorphic self-incompatibility, and the morph ratio within a population is very important for reproductive success. The aims of this study were to construct a linkage map, map the S locus as a qualitative trait and quantitative trait loci (QTLs) for floral morphological traits related to heterostyly, and predict the morph type in wild populations by using molecular markers for devising a conservation strategy.Methods A linkage map was constructed with 126 markers. The QTLs for four floral traits and the S locus were mapped. Using the genotypes of loci that were located near both the S locus and the QTLs with large effects, morphs of 59 wild genets were predicted.Key Results The linkage map consisted of 14 linkage groups (LGs). The S locus was mapped to LG 7. Major QTLs for stigma and anther heights were detected in the same region as the S locus. These QTLs exhibited high logarithm of the odds scores and explained a high percentage of the phenotypic variance (>85 %). By analysing these two traits within each morph, additional QTLs for each trait were detected. Using the four loci linked to the S locus, the morphs of 43 genets in three wild populations could be predicted.Conclusions This is the first report of a linkage map and QTL analysis for floral morphology related to heterostyly in P. sieboldii. Floral morphologies related to heterostyly are controlled by the S locus in LG 7 and by several QTLs in other LGs. Additionally, this study showed that molecular markers are effective tools for investigating morph ratios in a population containing the non-flowering individuals or during the non-flowering seasons

Genetic architecture and evolution of the S locus supergene in Primula vulgaris

Darwin’s studies on heterostyly in Primula described two floral morphs, pin and thrum, with reciprocal anther and stigma heights that promote insect-mediated cross-pollination. This key innovation evolved independently in several angiosperm families. Subsequent studies on heterostyly in Primula contributed to the foundation of modern genetic theory and the neo-Darwinian synthesis. The established genetic model for Primula heterostyly involves a diallelic S locus comprising several genes, with rare recombination events that result in self-fertile homostyle flowers with anthers and stigma at the same height. Here we reveal the S locus supergene as a tightly-linked cluster of thrum-specific genes that are absent in pins. We show that thrums are hemizygous not heterozygous for the S locus, which suggests that homostyles do not arise by recombination between S locus haplotypes as previously proposed. Duplication of a floral homeotic gene 51.7 MYA, followed by its neofunctionalisation, created the current S locus assemblage which led to floral heteromorphy in Primula. Our findings provide new insights into the structure, function and evolution of this archetypal supergene

Self-(in)compatibility inheritance and allele-specific marker development in yellow mustard (Sinapis alba)

Yellow mustard (Sinapis alba) has a sporophytic self-incompatibility reproduction system. Genetically stable self-incompatible (SI) and self-compatible (SC) inbred lines have recently been developed in this crop. Understanding the S haplotype of different inbred lines and the inheritance of the self-(in)compatibility (SI/SC) trait is very important for breeding purposes. In this study, we used the S-locus gene-specific primers in Brassica rapa and Brassica oleracea to clone yellow mustard S-locus genes of SI lines Y514 and Y1130 and SC lines Y1499 and Y1501. The PCR amplification results and DNA sequences of the S-locus genes revealed that Y514 carried the class I S haplotype, while Y1130, Y1499, and Y1501 had the class II S haplotype. The results of our genetic studies indicated that self-incompatibility was dominant over self-compatibility and controlled by a one-gene locus in the two crosses of Y514 × Y1499 and Y1130 × Y1501. Of the five S-locus gene polymorphic primer pairs, Sal-SLGI and Sal-SRKI each generated one dominant marker for the SI phenotype of Y514; Sal-SLGII and Sal-SRKII produced dominant marker(s) for the SC phenotype of Y1501 and Y1499; Sal-SP11II generated one dominant marker for Y1130. These markers co-segregated with the SI/SC phenotype in the F(2) populations of the two crosses. In addition, co-dominant markers were developed by mixing the two polymorphic primer pairs specific for each parent in the multiplex PCR, which allowed zygosity to be determined in the F(2) populations. The SI/SC allele-specific markers have proven to be very useful for the selection of the desirable SC genotypes in our yellow mustard breeding program. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11032-013-9943-8) contains supplementary material, which is available to authorized users

S-locus diversity of sweet cherry varieties from Galicia, North Western Spain (abstract)

Publishe

Comparative gene mapping in Arabidopsis lyrata chromosomes 6 and 7 and A thaliana chromosome IV: evolutionary history, rearrangements and local recombination rates

We have increased the density of genetic markers on the Arabidopsis lyrata chromosomes AL6 and AL7 corresponding to the A. thaliana chromosome IV, in order to determine chromosome rearrangements between these two species, and to compare recombination fractions across the same intervals. We confirm the two rearrangements previously inferred (a reciprocal translocation and a large inversion, which we infer to be pericentric). By including markers around the centromere regions of A. thaliana chromosomes IV and V, we localize the AL6 centromere, and can localize the breakpoints of these chromosome rearrangements more precisely than previously. One translocation breakpoint was close to the centromere, and the other coincided with one end of the inversion, suggesting that a single event caused both rearrangements. At the resolution of our mapping, apart from these rearrangements, all other markers are in the same order in A. lyrata and A. thaliana. We could thus compare recombination rates in the two species. We found slightly higher values in A. thaliana, and a minimum estimate for regions not close to a centromere in A. lyrata is 4-5 centimorgans per megabase. The mapped region of AL7 includes the self-incompatibility loci (S-loci), and this region has been predicted to have lower recombination than elsewhere in the genome. We mapped 17 markers in a region of 1·23 Mb surrounding these loci, and compared the approximately 600 kb closest to the S-loci with the surrounding region of approximately the same size. There were significantly fewer recombination events in the closer than the more distant region, supporting the above prediction, but showing that the low recombination region is very limited in size. © 2006 Cambridge University Press.</p