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

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

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

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

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

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

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

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

Comparative Genomics Elucidates the Origin of a Supergene Controlling Floral Heteromorphism.

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

Divergent evolution and niche differentiation within the common peatmoss Sphagnum magellanicum.

PREMISE OF THE STUDY: Populations with phenotypic polymorphism in discrete characters may be good models for investigating genome evolution and speciation. Sphagnum magellanicum Brid. is found throughout the northern hemisphere, and despite considerable variation in morphological characters, it is considered one of the least taxonomically controversial peatmoss species. We have observed two main morphs of the species associated with different microhabitats. Here we investigated the genomic and environmental basis of this intraspecific morphological variation. METHODS: We conducted transplant and common garden experiments to test whether the two morphs are genetically differentiated. We then used RAD-sequencing to quantify the genomic divergence between the morphs and approximate Bayesian computation (ABC) to infer the most likely demographic scenario explaining the genome-wide differentiation of the two morphs. KEY RESULTS: We found that genomic differentiation between the two morphs is unexpectedly high and that several of the differentiated morphological characters have a genetic basis. Using simulation approaches, we found support for a scenario of ancient divergence followed by recent secondary contact. CONCLUSIONS: We show that the two morphs represent the two main genetic clusters previously found worldwide. Our results demonstrate that relatively minor morphological differentiation in a presumed phenotypically plastic peatmoss may be associated with massive divergence across the genome.

Spatial Genetic Structure of the Abundant and Widespread Peatmoss Sphagnum magellanicum Brid.

Spore-producing organisms have small dispersal units enabling them to become widespread across continents. However, barriers to gene flow and cryptic speciation may exist. The common, haploid peatmoss Sphagnum magellanicum occurs in both the Northern and Southern hemisphere, and is commonly used as a model in studies of peatland ecology and peatmoss physiology. Even though it will likely act as a rich source in functional genomics studies in years to come, surprisingly little is known about levels of genetic variability and structuring in this species. Here, we assess for the first time how genetic variation in S. magellanicum is spatially structured across its full distribution range (Northern Hemisphere and South America). The morphologically similar species S. alaskense was included for comparison. In total, 195 plants were genotyped at 15 microsatellite loci. Sequences from two plastid loci (trnG and trnL) were obtained from 30 samples. Our results show that S. alaskense and almost all plants of S. magellanicum in the northern Pacific area are diploids and share the same gene pool. Haploid plants occur in South America, Europe, eastern North America, western North America, and southern Asia, and five genetically differentiated groups with different distribution ranges were found. Our results indicate that S. magellanicum consists of several distinct genetic groups, seemingly with little or no gene flow among them. Noteworthy, the geographical separation of diploids and haploids is strikingly similar to patterns found within other haploid Sphagnum species spanning the Northern Hemisphere. Our results confirm a genetic division between the Beringian and the Atlantic that seems to be a general pattern in Sphagnum taxa. The pattern of strong genetic population structuring throughout the distribution range of morphologically similar plants need to be considered in future functional genomic studies of S. magellanicum.