Supergene evolution via stepwise duplications and neofunctionalization of a floral-organ identity gene

The molecular basis of phenotypic evolution: beyond the usual suspects

It has been well documented that mutations in coding DNA or cis-regulatory elements underlie natural phenotypic variation in many organisms. However, the development of sophisticated functional tools in recent years in a wide range of traditionally non-model systems have revealed many 'unusual suspects' in the molecular bases of phenotypic evolution, including upstream open reading frames (uORFs), cryptic splice sites, and small RNAs. Furthermore, large-scale genome sequencing, especially long-read sequencing, has identified a cornucopia of structural variation underlying phenotypic divergence and elucidated the composition of supergenes that control complex multi-trait polymorphisms. In this review article we highlight recent studies that demonstrate this great diversity of molecular mechanisms producing adaptive genetic variation and the panoply of evolutionary paths leading to the 'grandeur of life'.

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.

A FLOWERING LOCUS T ortholog is associated with photoperiod-insensitive flowering in hemp (Cannabis sativa L.)

Hemp (Cannabis sativa L.) is an extraordinarily versatile crop, with applications ranging from medicinal compounds to seed oil and fibre products. Cannabis sativa is a short-day plant, and its flowering is highly controlled by photoperiod. However, substantial genetic variation exists for photoperiod sensitivity in C. sativa, and photoperiod-insensitive ("autoflower") cultivars are available. Using a bi-parental mapping population and bulked segregant analysis, we identified Autoflower2, a 0.5 Mbp locus significantly associated with photoperiod-insensitive flowering in hemp. Autoflower2 contains an ortholog of the central flowering time regulator FLOWERING LOCUS T (FT) from Arabidopsis thaliana which we termed CsFT1. We identified extensive sequence divergence between alleles of CsFT1 from photoperiod-sensitive and insensitive cultivars of C. sativa, including a duplication of CsFT1 and sequence differences, especially in introns. Furthermore, we observed higher expression of one of the CsFT1 copies found in the photoperiod-insensitive cultivar. Genotyping of several mapping populations and a diversity panel confirmed a correlation between CsFT1 alleles and photoperiod response, affirming that at least two independent loci involved in the photoperiodic control of flowering, Autoflower1 and Autoflower2, exist in the C. sativa gene pool. This study reveals the multiple independent origins of photoperiod insensitivity in C. sativa, supporting the likelihood of a complex domestication history in this species. By integrating the genetic relaxation of photoperiod sensitivity into novel C. sativa cultivars, expansion to higher latitudes will be permitted, thus allowing the full potential of this versatile crop to be reached.

PfPIN5 promotes style elongation by regulating cell length in Primula forbesii Franch

BACKGROUND AND AIMS

Style dimorphism is one of the polymorphic characteristics of flowers in heterostylous plants, which have two types of flowers: the pin morph, with long styles and shorter anthers, and the thrum morph, with short styles and longer anthers. The formation of dimorphic styles has received attention in the plant world. Previous studies showed that CYP734A50 in Primula determined style length and limited style elongation and that the brassinosteroid metabolic pathway was involved in regulation of style length. However, it is unknown whether there are other factors affecting the style length of Primula.

METHODS

Differentially expressed genes highly expressed in pin morph styles were screened based on Primula forbesii transcriptome data. Virus-induced gene silencing was used to silence these genes, and the style length and anatomical changes were observed 20 days after injection.

KEY RESULTS

PfPIN5 was highly expressed in pin morph styles. When PfPIN5 was silenced, the style length was shortened in pin and long-homostyle plants by shortening the length of style cells. Moreover, silencing CYP734A50 in thrum morph plants increased the expression level of PfPIN5 significantly, and the style length increased. The results indicated that PfPIN5, an auxin efflux transporter gene, contributed to regulation of style elongation in P. forbesii.

CONCLUSIONS

The results implied that the auxin pathway might also be involved in the formation of styles of P. forbesii, providing a new pathway for elucidating the molecular mechanism of style elongation in P. forbesii.

Genomic analyses elucidate S-locus evolution in response to intra-specific losses of distyly in Primula vulgaris

Distyly, a floral dimorphism that promotes outcrossing, is controlled by a hemizygous genomic region known as the S-locus. Disruptions of genes within the S-locus are responsible for the loss of distyly and the emergence of homostyly, a floral monomorphism that favors selfing. Using whole-genome resequencing data of distylous and homostylous individuals from populations of Primula vulgaris and leveraging high-quality reference genomes of Primula we tested, for the first time, predictions about the evolutionary consequences of transitions to selfing on S-genes. Our results reveal a previously undetected structural rearrangement in CYPᵀ associated with the shift to homostyly and confirm previously reported, homostyle-specific, loss-of-function mutations in the exons of the S-gene CYPᵀ. We also discovered that the promoter and intronic regions of CYPᵀ in distylous and homostylous individuals are conserved, suggesting that down-regulation of CYPᵀ via mutations in its promoter and intronic regions is not a cause of the shift to homostyly. Furthermore, we found that hemizygosity is associated with reduced genetic diversity in S-genes compared with their paralogs outside the S-locus. Additionally, the shift to homostyly lowers genetic diversity in both the S-genes and their paralogs, as expected in primarily selfing plants. Finally, we tested, for the first time, long-standing theoretical models of changes in S-locus genotypes during early stages of the transition to homostyly, supporting the assumption that two copies of the S-locus might reduce homostyle fitness.

Molecular insights into self-incompatibility systems: From evolution to breeding

Plants have evolved diverse self-incompatibility (SI) systems for outcrossing. Since Darwin's time, considerable progress has been made toward elucidating this unrivaled reproductive innovation. Recent advances in interdisciplinary studies and applications of biotechnology have given rise to major breakthroughs in understanding the molecular pathways that lead to SI, particularly the strikingly different SI mechanisms that operate in Solanaceae, Papaveraceae, Brassicaceae, and Primulaceae. These best-understood SI systems, together with discoveries in other "nonmodel" SI taxa such as Poaceae, suggest a complex evolutionary trajectory of SI, with multiple independent origins and frequent and irreversible losses. Extensive exploration of self-/nonself-discrimination signaling cascades has revealed a comprehensive catalog of male and female identity genes and modifier factors that control SI. These findings also enable the characterization, validation, and manipulation of SI-related factors for crop improvement, helping to address the challenges associated with development of inbred lines. Here, we review current knowledge about the evolution of SI systems, summarize key achievements in the molecular basis of pollen‒pistil interactions, discuss potential prospects for breeding of SI crops, and raise several unresolved questions that require further investigation.

A terpene synthase supergene locus determines chemotype in Melaleuca alternifolia (tea tree)

Leaf oil terpenes vary categorically in many plant populations, leading to discrete phenotypes of adaptive and economic significance, but for most species, a genetic explanation for the concerted fluctuation in terpene chemistry remains unresolved. To uncover the genetic architecture underlying multi-component terpene chemotypes in Melaleuca alternifolia (tea tree), a genome-wide association study was undertaken for 148 individuals representing all six recognised chemotypes. A number of single nucleotide polymorphisms in a genomic region of c. 400 kb explained large proportions of the variation in key monoterpenes of tea tree oil. The region contained a cluster of 10 monoterpene synthase genes, including four genes predicted to encode synthases for 1,8-cineole, terpinolene, and the terpinen-4-ol precursor, sabinene hydrate. Chemotype-dependent null alleles at some sites suggested structural variants within this gene cluster, providing a possible basis for linkage disequilibrium in this region. Genotyping in a separate domesticated population revealed that all alleles surrounding this gene cluster were fixed after artificial selection for a single chemotype. These observations indicate that a supergene accounts for chemotypes in M. alternifolia. A genetic model with three haplotypes, encompassing the four characterised monoterpene synthase genes, explained the six terpene chemotypes, and was consistent with available biparental cross-segregation data.

Haplotype-resolved genome assembly provides insights into the evolution of S-locus supergene in distylous Nymphoides indica

Distyly has evolved independently in numerous animal-pollinated angiosperm lineages. Understanding of its molecular basis has been restricted to a few species, primarily Primula. Here, we investigate the genetic architecture of the single diallelic locus (S-locus) supergene, a linkage group of functionally associated genes, and explore how it may have evolved in distylous Nymphoides indica, a lineage of flowering plants not previously investigated. We assembled haplotype-resolved genomes, used read-coverage-based genome-wide association study (rb-GWAS) to locate the S-locus supergene, co-expression network analysis to explore gene networks underpinning the development of distyly, and comparative genomic analyses to investigate the origins of the S-locus supergene. We identified three linked candidate S-locus genes - NinBAS1, NinKHZ2, and NinS1 - that were only evident in the short-styled morph and were hemizygous. Co-expression network analysis suggested that brassinosteroids contribute to dimorphic sex organs in the short-styled morph. Comparative genomic analyses indicated that the S-locus supergene likely evolved via stepwise duplications and has been affected by transposable element activities. Our study provides novel insight into the structure, regulation, and evolution of the supergene governing distyly in N. indica. It also provides high-quality genomic resources for future research on the molecular mechanisms underlying the striking evolutionary convergence in form and function across heterostylous taxa.

Colour polymorphism associated with a gene duplication in male wood tiger moths

Colour is often used as an aposematic warning signal, with predator learning expected to lead to a single colour pattern within a population. However, there are many puzzling cases where aposematic signals are also polymorphic. The wood tiger moth, Arctia plantaginis, displays bright hindwing colours associated with unpalatability, and males have discrete colour morphs which vary in frequency between localities. In Finland, both white and yellow morphs can be found, and these colour morphs also differ in behavioural and life-history traits. Here, we show that male colour is linked to an extra copy of a yellow family gene that is only present in the white morphs. This white-specific duplication, which we name valkea, is highly upregulated during wing development. CRISPR targeting valkea resulted in editing of both valkea and its paralog, yellow-e, and led to the production of yellow wings. We also characterise the pigments responsible for yellow, white, and black colouration, showing that yellow is partly produced by pheomelanins, while black is dopamine-derived eumelanin. Our results add to a growing number of studies on the genetic architecture of complex and seemingly paradoxical polymorphisms, and the role of gene duplications and structural variation in adaptive evolution.

A supergene in seaweed flies modulates male traits and female perception

Supergenes, tightly linked sets of alleles, offer some of the most spectacular examples of polymorphism persisting under long-term balancing selection. However, we still do not understand their evolution and persistence, especially in the face of accumulation of deleterious elements. Here, we show that an overdominant supergene in seaweed flies, Coelopa frigida, modulates male traits, potentially facilitating disassortative mating and promoting intraspecific polymorphism. Across two continents, the Cf-Inv(1) supergene strongly affected the composition of male cuticular hydrocarbons (CHCs) but only weakly affected CHC composition in females. Using gas chromatography-electroantennographic detection, we show that females can sense male CHCs and that there may be differential perception between genotypes. Combining our phenotypic results with RNA-seq data, we show that candidate genes for CHC biosynthesis primarily show differential expression for Cf-Inv(1) in males but not females. Conversely, candidate genes for odorant detection were differentially expressed in both sexes but showed high levels of divergence between supergene haplotypes. We suggest that the reduced recombination between supergene haplotypes may have led to rapid divergence in mate preferences as well as increasing linkage between male traits, and overdominant loci. Together this probably helped to maintain the polymorphism despite deleterious effects in homozygotes.

Two floral forms in the same species-distyly

MAIN CONCLUSION

This paper reviews the progress of research on the morphology, physiology and molecular biology of distyly in plants. It will help to elucidate the mysteries of distyly in plants. Distyly is a unique representative type of heterostyly in plants, primarily characterized by the presence of long style and short style within the flowers of the same species. This interesting trait has always fascinated researchers. With the rapid development of molecular biology, the molecular mechanism for the production of dimorphic styles in plants is also gaining ground. Researchers have been studying plant dimorphic styles from various perspectives. The researchers are gradually unravelling the mechanisms by which plants produce distyly traits. This paper reviews advances in the study of plant dimorphic style characteristics, mainly in terms of the morphology, physiology and molecular biology of plants with dimorphic styles. The aim is to provide a theoretical basis for the study of the mechanism of distyly formation in plants.

Genome sequencing reveals the genetic architecture of heterostyly and domestication history of common buckwheat

Common buckwheat, Fagopyrum esculentum, is an orphan crop domesticated in southwest China that exhibits heterostylous self-incompatibility. Here we present chromosome-scale assemblies of a self-compatible F. esculentum accession and a self-compatible wild relative, Fagopyrum homotropicum, together with the resequencing of 104 wild and cultivated F. esculentum accessions. Using these genomic data, we report the roles of transposable elements and whole-genome duplications in the evolution of Fagopyrum. In addition, we show that (1) the breakdown of heterostyly occurs through the disruption of a hemizygous gene jointly regulating the style length and female compatibility and (2) southeast Tibet was involved in common buckwheat domestication. Moreover, we obtained mutants conferring the waxy phenotype for the first time in buckwheat. These findings demonstrate the utility of our F. esculentum assembly as a reference genome and promise to accelerate buckwheat research and breeding.

Inheritance of distyly and homostyly in self-incompatible Primula forbesii

The evolutionary transition from self-incompatible distyly to self-compatible homostyly frequently occurs in heterostylous taxa. Although the inheritance of distyly and homostyly has been deeply studied, our understanding on modifications of the classical simple Mendelian model is still lacking. Primula forbesii, a biennial herb native to southwest China, is a typical distylous species, but after about 20 years of cultivation with open pollination, self-compatible homostyly appeared, providing ideal material for the study of the inheritance of distyly and homostyly. In this study, exogenous homobrassinolide was used to break the heteromorphic incompatibility of P. forbesii. Furthermore, we performed artificial pollination and open-pollination experiments to observe the distribution of floral morphs in progeny produced by different crosses. The viability of seeds from self-pollination was always the lowest among all crosses, and the homozygous S-morph plants (S/S) occurred in artificial pollination experiments but may experience viability selection. The distyly of P. forbesii is governed by a single S-locus, with S-morph dominant hemizygotes (S/-) and L-morph recessive homozygotes (-/-). Homostylous plants have a genotype similar to L-morph plants, and homostyly may be caused by one or more unlinked modifier genes outside the S-locus. Open pollinations confirm that autonomous self-pollination occurs frequently in L-morphs and homostylous plants. This study deepens the understanding of the inheritance of distyly and details a case of homostyly that likely originated from one or more modifier genes.

Annotation of the Turnera subulata (Passifloraceae) Draft Genome Reveals the S-Locus Evolved after the Divergence of Turneroideae from Passifloroideae in a Stepwise Manner

A majority of Turnera species (Passifloraceae) exhibit distyly, a reproductive system involving both self-incompatibility and reciprocal herkogamy. This system differs from self-incompatibility in Passiflora species. The genetic basis of distyly in Turnera is a supergene, restricted to the S-morph, and containing three S-genes. How supergenes and distyly evolved in Turnera, and the other Angiosperm families exhibiting distyly remain largely unknown. Unraveling the evolutionary origins in Turnera requires the generation of genomic resources and extensive phylogenetic analyses. Here, we present the annotated draft genome of the S-morph of distylous Turnera subulata. Our annotation allowed for phylogenetic analyses of the three S-genes' families across 56 plant species ranging from non-seed plants to eudicots. In addition to the phylogenetic analysis, we identified the three S-genes' closest paralogs in two species of Passiflora. Our analyses suggest that the S-locus evolved after the divergence of Passiflora and Turnera. Finally, to provide insights into the neofunctionalization of the S-genes, we compared expression patterns of the S-genes with close paralogs in Arabidopsis and Populus trichocarpa. The annotation of the T. subulata genome will provide a useful resource for future comparative work. Additionally, this work has provided insights into the convergent nature of distyly and the origin of supergenes.

Different molecular changes underlie the same phenotypic transition: Origins and consequences of independent shifts to homostyly within species

The repeated transition from outcrossing to selfing is a key topic in evolutionary biology. However, the molecular basis of such shifts has been rarely examined due to lack of knowledge of the genes controlling these transitions. A classic example of mating system transition is the repeated shift from heterostyly to homostyly. Occurring in 28 angiosperm families, heterostyly is characterized by the reciprocal position of male and female sexual organs in two (or three) distinct, usually self-incompatible floral morphs. Conversely, homostyly is characterized by a single, self-compatible floral morph with reduced separation of male and female organs, facilitating selfing. Here, we investigate the origins of homostyly in Primula vulgaris and its microevolutionary consequences by integrating surveys of the frequency of homostyles in natural populations, DNA sequence analyses of the gene controlling the position of female sexual organs (CYPᵀ), and microsatellite genotyping of both progeny arrays and natural populations characterized by varying frequencies of homostyles. As expected, we found that homostyles displace short-styled individuals, but long-style morphs are maintained at low frequencies within populations. We also demonstrated that homostyles repeatedly evolved from short-styled individuals in association with different types of loss-of-function mutations in CYPᵀ. Additionally, homostyly triggers a shift to selfing, promoting increased inbreeding within and genetic differentiation among populations. Our results elucidate the causes and consequences of repeated transitions to homostyly within species, and the putative mechanisms precluding its fixation in P. vulgaris. This study represents a benchmark for future analyses of losses of heterostyly in other angiosperms.

Genomic evidence supports the genetic convergence of a supergene controlling the distylous floral syndrome

Heterostyly, a plant sexual polymorphism controlled by the S-locus supergene, has evolved numerous times among angiosperm lineages and represents a classic example of convergent evolution in form and function. Determining whether underlying molecular convergence occurs could provide insights on constraints to floral evolution. Here, we investigated S-locus genes in distylous Gelsemium (Gelsemiaceae) to determine whether there is evidence of molecular convergence with unrelated distylous species. We used several approaches, including anatomical measurements of sex-organ development and transcriptome and whole-genome sequencing, to identify components of the S-locus supergene. We also performed evolutionary analysis with candidate S-locus genes and compared them with those reported in Primula and Turnera. The candidate S-locus supergene of Gelsemium contained four genes, of which three appear to have originated from gene duplication events within Gelsemiaceae. The style-length genes GeCYP in Gelsemium and CYP734A50 in Primula likely arose from duplication of the same gene, CYP734A1. Three out of four S-locus genes in Gelsemium elegans were hemizygous, as previously reported in Primula and Turnera. We provide genomic evidence on the genetic convergence of the supergene underlying distyly among distantly related angiosperm lineages and help to illuminate the genetic architecture involved in the evolution of heterostyly.

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 CYP^T 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.

The S-Gene YUC6 Pleiotropically Determines Male Mating Type and Pollen Size in Heterostylous Turnera (Passifloraceae): A Novel Neofunctionalization of the YUCCA Gene Family

In heterostylous, self-incompatible Turnera species, a member of the YUCCA gene family, YUC6, resides at the S-locus and has been hypothesized to determine the male mating type. YUCCA gene family members synthesize the auxin, indole-3-acetic acid, via a two-step process involving the TAA gene family. Consequently, it has been speculated that differences in auxin concentration in developing anthers are the biochemical basis underlying the male mating type. Here, we provide empirical evidence that supports this hypothesis. Using a transgenic knockdown approach, we show that YUC6 acts pleiotropically to control both the male physiological mating type and pollen size, but not the filament length dimorphism associated with heterostyly in Turnera. Using qPCR to assess YUC6 expression in different transgenic lines, we demonstrate that the level of YUC6 knockdown correlates with the degree of change observed in the male mating type. Further assessment of YUC6 expression through anther development, in the knockdown lines, suggests that the male mating type is irreversibly determined during a specific developmental window prior to microsporogenesis, which is consistent with the genetically sporophytic nature of this self-incompatibility system. These results represent the first gene controlling male mating type to be characterized in any species with heterostyly.

Genomic analyses of the Linum distyly supergene reveal convergent evolution at the molecular level

Supergenes govern multi-trait-balanced polymorphisms in a wide range of systems; however, our understanding of their origins and evolution remains incomplete. The reciprocal placement of stigmas and anthers in pin and thrum floral morphs of distylous species constitutes an iconic example of a balanced polymorphism governed by a supergene, the distyly S-locus. Recent studies have shown that the Primula and Turnera distyly supergenes are both hemizygous in thrums, but it remains unknown whether hemizygosity is pervasive among distyly S-loci. As hemizygosity has major consequences for supergene evolution and loss, clarifying whether this genetic architecture is shared among distylous species is critical. Here, we have characterized the genetic architecture and evolution of the distyly supergene in Linum by generating a chromosome-level genome assembly of Linum tenue, followed by the identification of the S-locus using population genomic data. We show that hemizygosity and thrum-specific expression of S-linked genes, including a pistil-expressed candidate gene for style length, are major features of the Linum S-locus. Structural variation is likely instrumental for recombination suppression, and although the non-recombining dominant haplotype has accumulated transposable elements, S-linked genes are not under relaxed purifying selection. Our findings reveal remarkable convergence in the genetic architecture and evolution of independently derived distyly supergenes, provide a counterexample to classic inversion-based supergenes, and shed new light on the origin and maintenance of an iconic floral polymorphism.

Radiation and hybridization underpin the spread of the fire ant social supergene

Some of the most striking polymorphisms in nature are regulated by “supergenes,” which are clusters of tightly linked genes that coordinately control complex phenotypes. Here, we study the evolutionary history of a supergene regulating colony social organization in fire ants. We show that the three inversions constituting the social supergene emerged sequentially during the separation of the ancestral lineages of Solenopsis invicta and Solenopsis richteri. Once completely assembled in S. richteri, the supergene introgressed into multiple closely related species despite recent hybridization being uncommon between several of the species. These findings provide a rare and striking example of how introgression can lead to the rapid spread of a novel variant controlling complex traits.

Supergenes are clusters of tightly linked genes that jointly produce complex phenotypes. Although widespread in nature, how such genomic elements are formed and how they spread are in most cases unclear. In the fire ant Solenopsis invicta and closely related species, a “social supergene controls whether a colony maintains one or multiple queens. Here, we show that the three inversions constituting the Social b (Sb) supergene emerged sequentially during the separation of the ancestral lineages of S. invicta and Solenopsis richteri. The two first inversions arose in the ancestral population of both species, while the third one arose in the S. richteri lineage. Once completely assembled in the S. richteri lineage, the supergene first introgressed into S. invicta, and from there into the other species of the socially polymorphic group of South American fire ant species. Surprisingly, the introgression of this large and important genomic element occurred despite recent hybridization being uncommon between several of the species. These results highlight how supergenes can readily move across species boundaries, possibly because of fitness benefits they provide and/or expression of selfish properties favoring their transmission.

Supergene potential of a selfish centromere

Selfishly evolving centromeres bias their transmission by exploiting the asymmetry of female meiosis and preferentially segregating to the egg. Such female meiotic drive systems have the potential to be supergenes, with multiple linked loci contributing to drive costs or enhancement. Here, we explore the supergene potential of a selfish centromere (D) in Mimulus guttatus, which was discovered in the Iron Mountain (IM) Oregon population. In the nearby Cone Peak population, D is still a large, non-recombining and costly haplotype that recently swept, but shorter haplotypes and mutational variation suggest a distinct population history. We detected D in five additional populations spanning more than 200 km; together, these findings suggest that selfish centromere dynamics are widespread in M. guttatus. Transcriptome comparisons reveal elevated differences in expression between driving and non-driving haplotypes within, but not outside, the drive region, suggesting large-scale cis effects of D's spread on gene expression. We use the expression data to refine linked candidates that may interact with drive, including Nuclear Autoantigenic Sperm Protein (NASP^SIM3), which chaperones the centromere-defining histone CenH3 known to modify Mimulus drive. Together, our results show that selfishly evolving centromeres may exhibit supergene behaviour and lay the foundation for future genetic dissection of drive and its costs. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

Association mapping of colour variation in a butterfly provides evidence that a supergene locks together a cluster of adaptive loci

Supergenes are genetic architectures associated with discrete and concerted variation in multiple traits. It has long been suggested that supergenes control these complex polymorphisms by suppressing recombination between sets of coadapted genes. However, because recombination suppression hinders the dissociation of the individual effects of genes within supergenes, there is still little evidence that supergenes evolve by tightening linkage between coadapted genes. Here, combining a landmark-free phenotyping algorithm with multivariate genome-wide association studies, we dissected the genetic basis of wing pattern variation in the butterfly Heliconius numata. We show that the supergene controlling the striking wing pattern polymorphism displayed by this species contains several independent loci associated with different features of wing patterns. The three chromosomal inversions of this supergene suppress recombination between these loci, supporting the hypothesis that they may have evolved because they captured beneficial combinations of alleles. Some of these loci are, however, associated with colour variations only in a subset of morphs where the phenotype is controlled by derived inversion forms, indicating that they were recruited after the formation of the inversions. Our study shows that supergenes and clusters of adaptive loci in general may form via the evolution of chromosomal rearrangements suppressing recombination between co-adapted loci but also via the subsequent recruitment of linked adaptive mutations. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

Characterization of PISTILLATA-like Genes and Their Promoters from the Distyly Fagopyrum esculentum

Arabidopsis PISTILLATA (PI) encodes B-class MADS-box transcription factor (TF), and works together with APETALA3 (AP3) to specify petal and stamen identity. However, a small-scale gene duplication event of PI ortholog was observed in common buckwheat and resulted in FaesPI_1 and FaesPI_2. FaesPI_1 and FaesPI_2 were expressed only in the stamen of dimorphic flower (thrum and pin) of Fagopyrum esculentum. Moreover, intense beta-glucuronidase (GUS) staining was found in the entire stamen (filament and anther) in pFaesPI_1::GUS transgenic Arabidopsis, while GUS was expressed only in the filament of pFaesPI_2::GUS transgenic Arabidopsis. In addition, phenotype complementation analysis suggested that pFaesPI_1::FaesPI_1/pFaesPI_2::FaesPI_2 transgenic pi-1 Arabidopsis showed similar a flower structure with stamen-like organs or filament-like organs in the third whorl. This suggested that FaesPI_2 only specified filament development, but FaesPI_1 specified stamen development. Meanwhile, FaesPI_1 and FaesPI_2 were shown to function redundantly in regulating filament development, and both genes work together to require a proper stamen identity. The data also provide a clue to understanding the roles of PI-like genes involved in floral organ development during the early evolution of core eudicots and also suggested that FaesPI_1 and FaesPI_2 hold the potential application in bioengineering to develop a common buckwheat male sterile line.

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?

Origin, loss, and regain of self-incompatibility in angiosperms

The self-incompatibility (SI) system with the broadest taxonomic distribution in angiosperms is based on multiple S-locus F-box genes (SLFs) tightly linked to an S-RNase termed type-1. Multiple SLFs collaborate to detoxify nonself S-RNases while being unable to detoxify self S-RNases. However, it is unclear how such a system evolved, because in an ancestral system with a single SLF, many nonself S-RNases would not be detoxified, giving low cross-fertilization rates. In addition, how the system has been maintained in the face of whole-genome duplications (WGDs) or lost in other lineages remains unclear. Here we show that SLFs from a broad range of species can detoxify S-RNases from Petunia with a high detoxification probability, suggestive of an ancestral feature enabling cross-fertilization and subsequently modified as additional SLFs evolved. We further show, based on its genomic signatures, that type-1 was likely maintained in many lineages, despite WGD, through deletion of duplicate S-loci. In other lineages, SI was lost either through S-locus deletions or by retaining duplications. Two deletion lineages regained SI through type-2 (Brassicaceae) or type-4 (Primulaceae), and one duplication lineage through type-3 (Papaveraceae) mechanisms. Thus, our results reveal a highly dynamic process behind the origin, maintenance, loss, and regain of SI.

Female self-incompatibility type in heterostylous Primula is determined by the brassinosteroid-inactivating cytochrome P450 CYP734A50

Most flowering plants are hermaphrodites, with flowers having both male and female reproductive organs. One widespread adaptation to limit self-fertilization is self-incompatibility (SI), where self-pollen fails to fertilize ovules.^1,2 In homomorphic SI, many morphologically indistinguishable mating types are found, although in heteromorphic SI, the two or three mating types are associated with different floral morphologies.^3-6 In heterostylous Primula, a hemizygous supergene determines a short-styled S-morph and a long-styled L-morph, corresponding to two different mating types, and full seed set only results from intermorph crosses.^7-9 Style length is controlled by the brassinosteroid (BR)-inactivating cytochrome P450 CYP734A50,¹⁰ yet it remains unclear what defines the male and female incompatibility types. Here, we show that CYP734A50 also determines the female incompatibility type. Inactivating CYP734A50 converts short S-morph styles into long styles with the same incompatibility behavior as L-morph styles, and this effect can be mimicked by exogenous BR treatment. In vitro responses of S- and L-morph pollen grains and pollen tubes to increasing BR levels could only partly explain their different in vivo behavior, suggesting both direct and indirect effects of the different BR levels in S- versus L-morph stigmas and styles in controlling pollen performance. This BR-mediated SI provides a novel mechanism for preventing self-fertilization. The joint control of morphology and SI by CYP734A50 has important implications for the evolutionary buildup of the heterostylous syndrome and provides a straightforward explanation for why essentially all of the derived self-compatible homostylous Primula species are long homostyles.¹¹.

Physalis floridana CRABS CLAW mediates neofunctionalization of GLOBOSA genes in carpel development

Floral B-function MADS-box genes, such as GLOBOSA (GLO), function in corolla and stamen organ identity specification. The functions of these genes outside these floral whorls are rarely reported. DOLL1 is a GLO gene controlling corolla and androecium organ identity. In this study we found that, in Physalis floridana double-layered-lantern 1 (doll1) mutant pollinated with wild-type pollen, fruit set was extremely low, indicating that doll1 females are dysfunctional. Stigma and style structure, stigma receptivity, pollen tube guidance, and embryo sac development were also impaired in doll1. P. floridana CRABS CLAW (PFCRC), predominantly expressed in carpels, was repressed in doll1 native carpels. Loss-of-function of PFCRC altered carpel meristem determinacy, carpel closure, and ovule number, and the resultant 'pistil' consisted of multiple spirally-arranged dorsiventral carpels occasionally with 1-2 naked ovules on the margin and trichomes at each mutated carpel tip, implying an alteration of carpel organ identity. Regulatory and genetic interactions between B-class MADS-box genes and PFCRC were revealed in a context-dependent manner in floral development. Our work reveals a new role for the B-function genes in carpel and ovule development via regulating PFCRC, providing a new understanding of genetic regulatory networks between MADS-domain and CRC transcription factors in mediating carpel organ specification, functionality, and origin.

Pistil Mating Type and Morphology Are Mediated by the Brassinosteroid Inactivating Activity of the S-Locus Gene BAHD in Heterostylous Turnera Species

Heterostyly is a breeding system that promotes outbreeding through a combination of morphological and physiological floral traits. In Turnera these traits are governed by a single, hemizygous S-locus containing just three genes. We report that the S-locus gene, BAHD, is mutated and encodes a severely truncated protein in a self-compatible long homostyle species. Further, a self-compatible long homostyle mutant possesses a T. krapovickasii BAHD allele with a point mutation in a highly conserved domain of BAHD acyl transferases. Wild type and mutant TkBAHD alleles were expressed in Arabidopsis to assay for brassinosteroid (BR) inactivating activity. The wild type but not mutant allele caused dwarfism, consistent with the wild type possessing, but the mutant allele having lost, BR inactivating activity. To investigate whether BRs act directly in self-incompatibility, BRs were added to in vitro pollen cultures of the two mating types. A small morph specific stimulatory effect on pollen tube growth was found with 5 µM brassinolide, but no genotype specific inhibition was observed. These results suggest that BAHD acts pleiotropically to mediate pistil length and physiological mating type through BR inactivation, and that in regard to self-incompatibility, BR acts by differentially regulating gene expression in pistils, rather than directly on pollen.

The maintenance of polymorphism in an ancient social supergene

Supergenes, regions of the genome with suppressed recombination between sets of functional mutations, contribute to the evolution of complex phenotypes in diverse systems. Excluding sex chromosomes, most supergenes discovered so far appear to be young, being found in one species or a few closely related species. Here, we investigate how a chromosome harbouring an ancient supergene has evolved over about 30 million years (Ma). The Formica supergene underlies variation in colony queen number in at least five species. We expand previous analyses of sequence divergence on this chromosome to encompass about 90 species spanning the Formica phylogeny. Within the nonrecombining region, the gene knockout contains 22 single nucleotide polymorphisms (SNPs) that are consistently differentiated between two alternative supergene haplotypes in divergent European Formica species, and we show that these same SNPs are present in most Formica clades. In these clades, including an early diverging Nearctic Formica clade, individuals with alternative genotypes at knockout also have higher differentiation in other portions of this chromosome. We identify hotspots of SNPs along this chromosome that are present in multiple Formica clades to detect genes that may have contributed to the emergence and maintenance of the genetic polymorphism. Finally, we infer three gene duplications on one haplotype, based on apparent heterozygosity within these genes in the genomes of haploid males. This study strengthens the evidence that this supergene originated early in the evolution of Formica and that just a few loci in this large region of suppressed recombination retain strongly differentiated alleles across contemporary Formica lineages.

A Draft Genome of the Ginger Species Alpinia nigra and New Insights into the Genetic Basis of Flexistyly

Angiosperms possess various strategies to ensure reproductive success, such as stylar polymorphisms that encourage outcrossing. Here, we investigate the genetic basis of one such dimorphism that combines both temporal and spatial separation of sexual function, termed flexistyly. It is a floral strategy characterised by the presence of two morphs that differ in the timing of stylar movement. We performed a de novo assembly of the genome of Alpinia nigra using high-depth genomic sequencing. We then used Pool-seq to identify candidate regions for flexistyly based on allele frequency or coverage differences between pools of anaflexistylous and cataflexistylous morphs. The final genome assembly size was 2 Gb, and showed no evidence of recent polyploidy. The Pool-seq did not reveal large regions with high FST values, suggesting large structural chromosomal polymorphisms are unlikely to underlie differences between morphs. Similarly, no region had a 1:2 mapping depth ratio which would be indicative of hemizygosity. We propose that flexistyly is governed by a small genomic region that might be difficult to detect with Pool-seq, or a complex genomic region that proved difficult to assemble. Our genome will be a valuable resource for future studies of gingers, and provides the first steps towards characterising this complex floral phenotype.

Whole-genome sequencing and genome regions of special interest: Lessons from major histocompatibility complex, sex determination, and plant self-incompatibility

Whole‐genome sequencing of non‐model organisms is now widely accessible and has allowed a range of questions in the field of molecular ecology to be investigated with greater power. However, some genomic regions that are of high biological interest remain problematic for assembly and data‐handling. Three such regions are the major histocompatibility complex (MHC), sex‐determining regions (SDRs) and the plant self‐incompatibility locus (S‐locus). Using these as examples, we illustrate the challenges of both assembling and resequencing these highly polymorphic regions and how bioinformatic and technological developments are enabling new approaches to their study. Mapping short‐read sequences against multiple alternative references improves genotyping comprehensiveness at the S‐locus thereby contributing to more accurate assessments of allelic frequencies. Long‐read sequencing, producing reads of several tens to hundreds of kilobase pairs in length, facilitates the assembly of such regions as single sequences can span the multiple duplicated gene copies of the MHC region, and sequence through repetitive stretches and translocations in SDRs and S‐locus haplotypes. These advances are adding value to short‐read genome resequencing approaches by allowing, for example, more accurate haplotype phasing across longer regions. Finally, we assessed further technical improvements, such as nanopore adaptive sequencing and bioinformatic tools using pangenomes, which have the potential to further expand our knowledge of a number of genomic regions that remain challenging to study with classical resequencing approaches.

Correlational selection in the age of genomics

Ecologists and evolutionary biologists are well aware that natural and sexual selection do not operate on traits in isolation, but instead act on combinations of traits. This long-recognized and pervasive phenomenon is known as multivariate selection, or-in the particular case where it favours correlations between interacting traits-correlational selection. Despite broad acknowledgement of correlational selection, the relevant theory has often been overlooked in genomic research. Here, we discuss theory and empirical findings from ecological, quantitative genetic and genomic research, linking key insights from different fields. Correlational selection can operate on both discrete trait combinations and quantitative characters, with profound implications for genomic architecture, linkage, pleiotropy, evolvability, modularity, phenotypic integration and phenotypic plasticity. We synthesize current knowledge and discuss promising research approaches that will enable us to understand how correlational selection shapes genomic architecture, thereby linking quantitative genetic approaches with emerging genomic methods. We suggest that research on correlational selection has great potential to integrate multiple fields in evolutionary biology, including developmental and functional biology, ecology, quantitative genetics, phenotypic polymorphisms, hybrid zones and speciation processes.

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.

Evolutionary Divergence and Biased Expression of NAC Transcription Factors in Hexaploid Bread Wheat (Triticum aestivum L.)

The NAC genes, a large plant-specific family of transcription factors, regulate a wide range of pathways involved in development and response to biotic and abiotic stress. In this study, the NAC transcription factors were identified in 27 green plants, and the results showed that NAC transcription factors in plants undergo an appearance stage from water to land and a number expansion stage from gymnosperm to angiosperm. Investigating the evolutionary process of the NAC transcription factors from diploid species to hexaploid wheat revealed that tandem replications during the polyploidization process is an important event for increasing the number of NAC transcription factors in wheat. Then, the molecular characteristics, phylogenetic relationships, and expression patterns of 462 NAC transcription factors of hexaploid wheat (TaNACs) were analyzed. The protein structure results showed that TaNAC was relatively conservative at the N-terminal that contains five subdomains. All these TaNACs were divided into Group I and Group II by phylogenetic analysis, and the TaNACs in Group I should undergo strong artificial selection based on single nucleotide polymorphism (SNP) analysis. Through genome synteny and phylogenetic analysis, these TaNACs were classified into 88 groups and 9 clusters. The biased expression results of these TaNACs showed that there are 24 groups and 67 groups of neofunctionalization genes under biotic and abiotic stress, respectively, and 16 groups and 59 groups of subfunctionalization genes. This shows that neofunctionalization plays an important role in coping with different stresses. Our study provides new insights into the evolution of NAC transcription factors in hexaploid wheat.

Characterization of 30 microsatellite markers for distylous Primula denticulata (Primulaceae) using HiSeq sequencing

Primula denticulata exhibits considerable variation in floral morphology and flowering phenology along elevational gradients in SW China. We isolated 30 microsatellite markers from P. denticulata to facilitate further investigation of population genetics and floral evolution in this species. We used the HiSeq X-Ten sequencing system to develop a set of markers, and measured polymorphism and genetic diversity in a sample of 72 individuals from three natural populations of P. denticulata subsp. denticulata. The markers displayed relatively high polymorphism, with the number of alleles ranging from two to seven (mean = 3.567). The observed and expected heterozygosity ranged from 0 to 1.000 and 0.041 to 0.702, respectively. Twenty-eight of the loci were also successfully amplified in P. denticulata subsp. sinodenticulata. The microsatellite markers we have identified will provide valuable tools for investigations of the population genetic structure, mating systems and phylogeography of the P. denticulata complex, and will help to address questions concerning the ecological and genetic mechanisms responsible for the evolution of reproductive traits in the species.

The Turnera Style S-Locus Gene TsBAHD Possesses Brassinosteroid-Inactivating Activity When Expressed in Arabidopsis thaliana

Heterostyly distinct hermaphroditic floral morphs enforce outbreeding. Morphs differ structurally, promote cross-pollination, and physiologically block self-fertilization. In Turnera the self-incompatibility (S)-locus controlling heterostyly possesses three genes specific to short-styled morph genomes. Only one gene, TsBAHD, is expressed in pistils and this has been hypothesized to possess brassinosteroid (BR)-inactivating activity. We tested this hypothesis using heterologous expression in Arabidopsis thaliana as a bioassay, thereby assessing growth phenotype, and the impacts on the expression of endogenous genes involved in BR homeostasis and seedling photomorphogenesis. Transgenic A. thaliana expressing TsBAHD displayed phenotypes typical of BR-deficient mutants, with phenotype severity dependent on TsBAHD expression level. BAS1, which encodes an enzyme involved in BR inactivation, was downregulated in TsBAHD-expressing lines. CPD and DWF, which encode enzymes involved in BR biosynthesis, were upregulated. Hypocotyl growth of TsBAHD dwarfs responded to application of brassinolide in light and dark in a manner typical of plants over-expressing genes encoding BR-inactivating activity. These results provide empirical support for the hypothesis that TsBAHD possesses BR-inactivating activity. Further this suggests that style length in Turnera is controlled by the same mechanism (BR inactivation) as that reported for Primula, but using a different class of enzyme. This reveals interesting convergent evolution in a biochemical mechanism to regulate floral form in heterostyly.