GENETIC ARCHITECTURE OF SEXUAL DIMORPHISM IN A SUBDIOECIOUS PLANT WITH A PROTO-SEX CHROMOSOME

Multi-model GWAS reveals key loci for horticultural traits in reconstructed garden strawberry

The cultivated garden strawberry (Fragaria × ananassa) has a rich history, originating from the hybridization of two wild octoploid strawberry species in the 18th century. Two-step reconstruction of Fragaria × ananassa through controlled crossings between pre-improved selections of its parental species is a promising approach for enriching the breeding germplasm of strawberry for wider adaptability. We created a population of reconstructed strawberry by hybridizing elite selections of F. virginiana and F. chiloensis. A replicated field experiment was conducted to evaluate the population's performance for eleven horticulturally important traits, over multiple years. Population structure analyses based on Fana-50 k SNP array data confirmed pedigree-based grouping of the progenies into four distinct groups. As complex traits are often influenced by environmental variables, and population structure can lead to spurious associations, we tested multiple genome-wide association study (GWAS) models. GWAS uncovered 39 quantitative trait loci (QTL) regions for eight traits distributed across twenty chromosomes, including 11 consistent and 28 putative QTLs. Candidate genes for traits including winter survival, flowering time, runnering vigor, and hermaphrodism were identified within the QTL regions. To our knowledge, this study marks the first comprehensive investigation of adaptive and horticultural traits in a large, multi-familial reconstructed strawberry population using SNP markers.

The role of mitochondria in sex- and age-specific gene expression in a species without sex chromosomes

Mitochondria perform an array of functions, many of which involve interactions with gene products encoded by the nucleus. These mitochondrial functions, particularly those involving energy production, can be expected to differ between sexes and across ages. Here, we measured mitochondrial effects on sex- and age-specific gene expression in parental and reciprocal F1 hybrids between allopatric populations of Tigriopus californicus with over 20% mitochondrial DNA divergence. Because the species lacks sex chromosomes, sex-biased mitochondrial effects are not confounded by the effects of sex chromosomes. Results revealed pervasive sex differences in mitochondrial effects, including effects on energetics and aging involving nuclear interactions throughout the genome. Using single-individual RNA sequencing, sex differences were found to explain more than 80% of the variance in gene expression. Males had higher expression of mitochondrial genes and mitochondrially targeted proteins (MTPs) involved in oxidative phosphorylation (OXPHOS), while females had elevated expression of non-OXPHOS MTPs, indicating strongly sex-dimorphic energy metabolism at the whole organism level. Comparison of reciprocal F1 hybrids allowed insights into the nature of mito-nuclear interactions, showing both mitochondrial effects on nuclear expression, and nuclear effects on mitochondrial expression. While based on a small set of crosses, sex-specific increases in mitochondrial expression with age were associated with longer life. Network analyses identified nuclear components of strong mito-nuclear interactions and found them to be sexually dimorphic. These results highlight the profound impact of mitochondria and mito-nuclear interactions on sex- and age-specific gene expression.

The role of mitochondria in sex- and age-specific gene expression in a species without sex chromosomes

Mitochondria perform an array of functions, many of which involve interactions with gene products encoded by the nucleus. These mitochondrial functions, particularly those involving energy production, can be expected to differ between sexes and across ages. Here we measured mitochondrial effects on sex- and age-specific gene expression in parental and reciprocal F1 hybrids between allopatric populations of Tigriopus californicus with over 20% mitochondrial DNA divergence. Because the species lacks sex chromosomes, sex-biased mitochondrial effects are not confounded by the effects of sex chromosomes. Using single-individual RNA sequencing, sex differences were found to explain more than 80% of the variance in gene expression. Males had higher expression of mitochondrial genes and mitochondrially targeted proteins (MTPs) involved in oxidative phosphorylation (OXPHOS), while females had elevated expression of non-OXPHOS MTPs, indicating strongly sex-dimorphic energy metabolism at the whole organism level. Comparison of reciprocal F1 hybrids allowed insights into the nature of mito-nuclear interactions, showing both mitochondrial effects on nuclear expression, as well as nuclear effects on mitochondrial expression. Across both sexes, increases in mitochondrial expression with age were associated with longer life. Network analyses identified nuclear components of strong mito-nuclear interactions, and found them to be sexually dimorphic. These results highlight the profound impact of mitochondria and mito-nuclear interactions on sex- and age-specific gene expression.

The maintenance of polygenic sex determination depends on the dominance of fitness effects which are predictive of the role of sexual antagonism

In species with polygenic sex determination (PSD), multiple male- and female-determining loci on different proto-sex chromosomes segregate as polymorphisms within populations. The extent to which these polymorphisms are at stable equilibria is not yet resolved. Previous work demonstrated that PSD is most likely to be maintained as a stable polymorphism when the proto-sex chromosomes have opposite (sexually antagonistic) fitness effects in males and females. However, these models usually consider PSD systems with only two proto-sex chromosomes, or they do not broadly consider the dominance of the alleles under selection. To address these shortcomings, I used forward population genetic simulations to identify selection pressures that can maintain PSD under different dominance scenarios in a system with more than two proto-sex chromosomes (modeled after the house fly). I found that overdominant fitness effects of male-determining proto-Y chromosomes are more likely to maintain PSD than dominant, recessive, or additive fitness effects. The overdominant fitness effects that maintain PSD tend to have proto-Y chromosomes with sexually antagonistic effects (male-beneficial and female-detrimental). In contrast, dominant fitness effects that maintain PSD tend to have sexually antagonistic multi-chromosomal genotypes, but the individual proto-sex chromosomes do not have sexually antagonistic effects. These results demonstrate that sexual antagonism can be an emergent property of the multi-chromosome genotype without individual sexually antagonistic chromosomes. My results further illustrate how the dominance of fitness effects has consequences for both the likelihood that PSD will be maintained as well as the role sexually antagonistic selection is expected to play in maintaining the polymorphism.

Sex differences in early transcriptomic responses to oxidative stress in the copepod Tigriopus californicus

Background

Patterns of gene expression can be dramatically different between males and females of the same species, in part due to genes on sex chromosomes. Here we test for sex differences in early transcriptomic response to oxidative stress in a species which lacks heteromorphic sex chromosomes, the copepod Tigriopus californicus.

Results

Male and female individuals were separately exposed to control conditions and pro-oxidant conditions (hydrogen peroxide and paraquat) for periods of 3 hours and 6 hours. Variance partitioning showed the greatest expression variance among individuals, highlighting the important information that can be obscured by the common practice of pooling individuals. Gene expression variance between sexes was greater than that among treatments, showing the profound effect of sex even when males and females share the same genome. Males exhibited a larger response to both pro-oxidants, differentially expressing more than four times as many genes, including up-regulation of more antioxidant genes, heat shock proteins and protease genes. While females differentially expressed fewer genes, the magnitudes of fold change were generally greater, indicating a more targeted response. Although females shared a smaller fraction of differentially expressed genes between stressors and time points, expression patterns of antioxidant and protease genes were more similar between stressors and more GO terms were shared between time points.

Conclusions

Early transcriptomic responses to the pro-oxidants H₂O₂ and paraquat in copepods revealed substantial variation among individuals and between sexes. The finding of such profound sex differences in oxidative stress response, even in the absence of sex chromosomes, highlights the importance of studying both sexes and the potential for developing sex-specific strategies to promote optimal health and aging in humans.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07179-5.

YY males of the dioecious plant Mercurialis annua are fully viable but produce largely infertile pollen

The suppression of recombination during sex-chromosome evolution is thought to be favoured by linkage between the sex-determining locus and sexually antagonistic loci, and leads to the degeneration of the chromosome restricted to the heterogametic sex. Despite substantial evidence for genetic degeneration at the sequence level, the phenotypic effects of the earliest stages of sex-chromosome evolution are poorly known. Here, we compare the morphology, viability and fertility between XY and YY individuals produced by crossing seed-producing males in the dioecious plant Mercurialis annua, which has young sex chromosomes with limited X-Y sequence divergence. We found no significant difference in viability or vegetative morphology between XY and YY males. However, electron microscopy revealed clear differences in pollen anatomy, and YY males were significantly poorer sires in competition with their XY counterparts. Our study suggests either that the X chromosome is required for full male fertility in M. annua, or that male fertility is sensitive to the dosage of relevant Y-linked genes. We discuss the possibility that the maintenance of male-fertility genes on the X chromosome might have been favoured in recent population expansions that selected for the ability of females to produce pollen in the absence of males.

Insights into the genetic architecture of sexual dimorphism from an interspecific cross between two diverging Silene (Caryophyllaceae) species

The evolution of sexual dimorphism in species with separate sexes is influenced by the resolution of sexual conflicts creating sex differences through genetic linkage or sex-biased expression. Plants with different degrees of sexual dimorphism are thus ideal to study the genetic basis of sexual dimorphism. In this study we explore the genetic architecture of sexual dimorphism between Silene latifolia and Silene dioica. These species have chromosomal sex determination and differ in the extent of sexual dimorphism. To test whether QTL for sexually dimorphic traits have accumulated on the sex chromosomes and to quantify their contribution to species differences, we create a linkage map and performed QTL analysis of life history, flower and vegetative traits using an unidirectional interspecific F2 hybrid cross. We found support for an accumulation of QTL on the sex chromosomes and that sex differences explained a large proportion of the variance between species, suggesting that both natural and sexual selection contributed to species divergence. Sexually dimorphic traits that also differed between species displayed transgressive segregation. We observed a reversal in sexual dimorphism in the F2 population, where males tended to be larger than females, indicating that sexual dimorphism is constrained within populations but not in recombinant hybrids. This study contributes to the understanding of the genetic basis of sexual dimorphism and its evolution in Silene.

Effects of oxidative stress on sex-specific gene expression in the copepod Tigriopus californicus revealed by single individual RNA-seq

Oxidative stress reflects the imbalance of pro-oxidants and antioxidants. Prolonged oxidative stress can induce cellular damage, diseases and aging, and the effects may be sex-specific. Tigriopus californicus has recently been proposed as an alternative model system for sex-specific studies due to the absence of sex chromosomes. In this study, we used comparative transcriptomic analyses to assess sex-specific transcriptional responses to oxidative stress. Male and female individuals were maintained separately in one of three treatments: 1) control conditions with an algae diet, 2) pro-oxidant (H₂O₂) conditions with an algae diet or 3) decreased antioxidant conditions (reduced carotenoids due to a yeast diet). Single individual RNA-seq was then conducted for twenty-four libraries using Ligation Mediated RNA sequencing (LM-Seq). Variance in gene expression was partitioned into 62.3% between sexes, 26.85% among individuals and 10.85% among treatments. Within each of the three treatments, expression was biased toward females. However, compared to the control treatment, males in both pro-oxidant and decreased antioxidant treatments differentially expressed more genes while females differentially expressed fewer genes but with a greater magnitude of fold change. As the first study of copepods to apply single individual RNA-seq, the findings will contribute to a better understanding of transcriptomic variation among individuals as well as sex-specific response mechanisms to oxidative stress in the absence of sex chromosomes.

Sex-specific stress tolerance, proteolysis, and lifespan in the invertebrate Tigriopus californicus

Because stress tolerance and longevity are mechanistically and phenotypically linked, the sex with higher acute stress tolerance might be expected to also live longer. On the other hand, the association between stress tolerance and lifespan may be complicated by tradeoffs between acute tolerance and long-term survival. Here we use the copepod Tigriopus californicus to test for sex differences in stress resistance, proteolytic activity and longevity. Unlike many model organisms, this species does not have sex chromosomes. However, substantial sex differences were still observed. Females were found to have superior tolerance to a range of acute stressors (high temperature, high salinity, low salinity, copper and bisphenol A (BPA)) across a variety of treatments including different populations, pure vs. hybrid crosses, and different shading environments. Upregulation of proteolytic capacity - one molecular mechanism for responding to acute stress - was also found to be sexually dimorphic. In the combined stress treatment of chronic copper exposure followed by acute heat exposure, proteolytic capacity was suppressed for males. Females, however, maintained a robust proteolytic stress response. While females consistently showed greater tolerance to short-term stress, lifespan was largely equivalent between the two sexes under both benign conditions and mild thermal stress. Our findings indicate that short-term stress tolerance does not predict long-term survival under relatively mild conditions.

Repeated translocation of a gene cassette drives sex-chromosome turnover in strawberries

Turnovers of sex-determining systems represent important diversifying forces across eukaryotes. Shifts in sex chromosomes—but conservation of the master sex-determining genes—characterize distantly related animal lineages. Yet in plants, in which separate sexes have evolved repeatedly and sex chromosomes are typically homomorphic, we do not know whether such translocations drive sex-chromosome turnovers within closely related taxonomic groups. This phenomenon can only be demonstrated by identifying sex-associated nucleotide sequences, still largely unknown in plants. The wild North American octoploid strawberries (Fragaria) exhibit separate sexes (dioecy) with homomorphic, female heterogametic (ZW) inheritance, yet sex maps to three different chromosomes in different taxa. To characterize these turnovers, we identified sequences unique to females and assembled their reads into contigs. For most octoploid Fragaria taxa, a short (13 kb) sequence was observed in all females and never in males, implicating it as the sex-determining region (SDR). This female-specific “SDR cassette” contains both a gene with a known role in fruit and pollen production and a novel retrogene absent on Z and autosomal chromosomes. Phylogenetic comparison of SDR cassettes revealed three clades and a history of repeated translocation. Remarkably, the translocations can be ordered temporally due to the capture of adjacent sequence with each successive move. The accumulation of the “souvenir” sequence—and the resultant expansion of the hemizygous SDR over time—could have been adaptive by locking genes into linkage with sex. Terminal inverted repeats at the insertion borders suggest a means of movement. To our knowledge, this is the first plant SDR shown to be translocated, and it suggests a new mechanism (“move-lock-grow”) for expansion and diversification of incipient sex chromosomes.

Sex chromosomes frequently restructure themselves during organismal evolution, often becoming highly differentiated. This dynamic process is poorly understood for most taxa, especially during the early stages typical of many dioecious flowering plants. We show that in wild strawberries, a female-specific region of DNA is associated with sex and has repeatedly changed its genomic location, each time increasing the size of the hemizygous female-specific sequence on the W sex chromosome. This observation shows, for the first time to our knowledge, that plant sex regions can “jump” and suggests that this phenomenon may be adaptive by gathering and locking new genes into linkage with sex. This conserved and presumed causal sex-determining sequence, which varies in both genomic location and degree of differentiation, will facilitate future studies to understand how sex chromosomes first begin to differentiate.

Gene buddies: linked balanced polymorphisms reinforce each other even in the absence of epistasis

The fates of genetic polymorphisms maintained by balancing selection depend on evolutionary dynamics at linked sites. While coevolution across linked, epigenetically-interacting loci has been extensively explored, such supergenes may be relatively rare. However, genes harboring adaptive variation can occur in close physical proximity while generating independent effects on fitness. Here, I present a model in which two linked loci without epistasis are both under balancing selection for unrelated reasons. Using forward-time simulations, I show that recombination rate strongly influences the retention of adaptive polymorphism, especially for intermediate selection coefficients. A locus is more likely to retain adaptive variation if it is closely linked to another locus under balancing selection, even if the two loci have no interaction. Thus, two linked polymorphisms can both be retained indefinitely even when they would both be lost to drift if unlinked. While these results may be intuitive, they have important implications for genetic architecture: clusters of mutually reinforcing genes may underlie phenotypic variation in natural populations, and such genes cannot be assumed to be functionally associated. Future studies that measure selection coefficients and recombination rates among closely linked genes will be fruitful for characterizing the extent of this phenomenon.

Evolution of recombination rates between sex chromosomes

In species with genetic sex-determination, the chromosomes carrying the sex-determining genes have often evolved non-recombining regions and subsequently evolved the full set of characteristics denoted by the term 'sex chromosomes'. These include size differences, creating chromosomal heteromorphism, and loss of gene functions from one member of the chromosome pair. Such characteristics and changes have been widely reviewed, and underlie molecular genetic approaches that can detect sex chromosome regions. This review deals mainly with the evolution of new non-recombining regions, focusing on how certain evolutionary situations select for suppressed recombination (rather than the proximate mechanisms causing suppressed recombination between sex chromosomes). Particularly important is the likely involvement of sexually antagonistic polymorphisms in genome regions closely linked to sex-determining loci. These may be responsible for the evolutionary strata of sex chromosomes that have repeatedly formed by recombination suppression evolving across large genome regions. More studies of recently evolved non-recombining sex-determining regions should help to test this hypothesis empirically, and may provide evidence about whether other situations can sometimes lead to sex-linked regions evolving. Similarities with other non-recombining genome regions are discussed briefly, to illustrate common features of the different cases, though no general properties apply to all of them.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

Genetic Mapping and Phylogenetic Analysis Reveal Intraspecific Variation in Sex Chromosomes of the Virginian Strawberry

With their extraordinary diversity in sexual systems, flowering plants offer unparalleled opportunities to understand sex determination and to reveal generalities in the evolution of sex chromosomes. Comparative genetic mapping of related taxa with good phylogenetic resolution can delineate the extent of sex chromosome diversity within plant groups, and lead the way to understanding the evolutionary drivers of such diversity. The North American octoploid wild strawberries provide such an opportunity. We performed linkage mapping using targeted sequence capture for the subdioecious western Fragaria virginiana ssp. platypetala and compared the location of its sex-determining region (SDR) to those of 2 other (sub)dioecious species, the eastern subspecies, F. virginiana ssp. virginiana (whose SDR is at 0-5.5 Mb on chromosome VI of the B2 subgenome), and the sister species F. chiloensis (whose SDR is at 37 Mb on chromosome VI of the Av subgenome). Male sterility was dominant in F. virginiana ssp. platypetala and mapped to a chromosome also in homeologous group VI. Likewise, one major quantitative trait locus (QTL) for female fertility overlapped the male sterility region. However, the SDR mapped to yet another subgenome (B1), and to a different location (13 Mb), but similar to the location inferred in one population of the naturally occurring hybrid between F. chiloensis and F. virginiana (F. ×ananassa ssp. cuneifolia). Phylogenetic analysis of chromosomes across the octoploid taxa showed consistent subgenomic composition reflecting shared evolutionary history but also reinforced within-species variation in the SDR-carrying chromosome, suggesting either repeated evolution, or recent turnovers in SDR.

Present-day sympatry belies the evolutionary origin of a high-order polyploid

Disentangling the evolutionary histories of polyploids, especially those with high ploidies, can reveal fundamental processes in speciation. Despite occurring frequently during evolution, the origins of many extant polyploid plant species remain largely unknown. By integrating linkage mapping, polyploid phylogeny and sex-determining region (SDR) in a unified framework, we statistically evaluated evolutionary hypotheses concerning the origin of a recently recognized decaploid strawberry (Fragaria cascadensis). The maximum-likelihood phylogenies and topology tests across homeologous groups consistently rejected the seemingly parsimonious hypothesis of 'contemporary sympatric speciation' via hybridization between octoploid and diploid congeners. Instead, most chromosomes supported 'ancient hybrid speciation' between a maternal octoploid progenitor ancestral to extant octoploid strawberries and a paternal, extinct Fragaria iinumae-like diploid progenitor, probably in Beringia during the Pleistocene. The absence of a shared SDR between the decaploid and other Fragaria is also consistent with an older origin rather than a recent hybrid origin in situ. Our study reveals a long evolutionary history of the decaploid despite its recent discovery, and highlights the pitfalls of inferring polyploid origins from niche/range alone or combined with morphology. It can serve as an exemplary starting step towards building much-needed model systems of established polyploids that have been, and remain to be, recognized.

A Century of Sex Determination in Flowering Plants

Plants have evolved a diverse array of strategies for sexual reproduction, particularly through the modification of male and female organs at distinct points in development. The immense variation in sexual systems across the land plants provides a unique opportunity to study the genetic, epigenetic, phylogenetic, and ecological underpinnings of sex determination. Here, we reflect on more than a century of research into flowering plant sex determination, placing a particular focus on the foundational genetic and cytogenetic observations, experiments, and hypotheses. Building on the seminal work on the genetics of plant sex, modern comparative genomic analyses now allow us to address longstanding questions about sex determination and the origins of sex chromosomes.

Homomorphic ZW chromosomes in a wild strawberry show distinctive recombination heterogeneity but a small sex-determining region

Recombination in ancient, heteromorphic sex chromosomes is typically suppressed at the sex-determining region (SDR) and proportionally elevated in the pseudoautosomal region (PAR). However, little is known about recombination dynamics of young, homomorphic plant sex chromosomes. We examine male and female function in crosses and unrelated samples of the dioecious octoploid strawberry Fragaria chiloensis in order to map the small and recently evolved SDR controlling both traits and to examine recombination patterns on the incipient ZW chromosome. The SDR of this ZW system is located within a 280 kb window, in which the maternal recombination rate is lower than the paternal one. In contrast to the SDR, the maternal PAR recombination rate is much higher than the rates of the paternal PAR or autosomes, culminating in an elevated chromosome-wide rate. W-specific divergence is elevated within the SDR and a single polymorphism is observed in high species-wide linkage disequilibrium with sex. Selection for recombination suppression within the small SDR may be weak, but fluctuating sex ratios could favor elevated recombination in the PAR to remove deleterious mutations on the W. The recombination dynamics of this nascent sex chromosome with a modestly diverged SDR may be typical of other dioecious plants.

Plant Sex Chromosomes

Although individuals in most flowering plant species, and in many haploid plants, have both sex functions, dioecious species-in which individuals have either male or female functions only-are scattered across many taxonomic groups, and many species have genetic sex determination. Among these, some have visibly heteromorphic sex chromosomes, and molecular genetic studies are starting to uncover sex-linked markers in others, showing that they too have fully sex-linked regions that are either too small or are located in chromosomes that are too small to be cytologically detectable from lack of pairing, lack of visible crossovers, or accumulation of heterochromatin. Detailed study is revealing that, like animal sex chromosomes, plant sex-linked regions show evidence for accumulation of repetitive sequences and genetic degeneration. Estimating when recombination stopped confirms the view that many plants have young sex-linked regions, making plants of great interest for studying the timescale of these changes.

Multilocus Sex Determination Revealed in Two Populations of Gynodioecious Wild Strawberry, Fragaria vesca subsp. bracteata

Gynodioecy, the coexistence of females and hermaphrodites, occurs in 20% of angiosperm families and often enables transitions between hermaphroditism and dioecy. Clarifying mechanisms of sex determination in gynodioecious species can thus illuminate sexual system evolution. Genetic determination of gynodioecy, however, can be complex and is not fully characterized in any wild species. We used targeted sequence capture to genetically map a novel nuclear contributor to male sterility in a self-pollinated hermaphrodite of Fragaria vesca subsp. bracteata from the southern portion of its range. To understand its interaction with another identified locus and possibly additional loci, we performed crosses within and between two populations separated by 2000 km, phenotyped the progeny and sequenced candidate markers at both sex-determining loci. The newly mapped locus contains a high density of pentatricopeptide repeat genes, a class commonly involved in restoration of fertility caused by cytoplasmic male sterility. Examination of all crosses revealed three unlinked epistatically interacting loci that determine sexual phenotype and vary in frequency between populations. Fragaria vesca subsp. bracteata represents the first wild gynodioecious species with genomic evidence of both cytoplasmic and nuclear genes in sex determination. We propose a model for the interactions between these loci and new hypotheses for the evolution of sex determining chromosomes in the subdioecious and dioecious Fragaria.

Genetics of dioecy and causal sex chromosomes in plants

Dioecy (separate male and female individuals) ensures outcrossing and is more prevalent in animals than in plants. Although it is common in bryophytes and gymnosperms, only 5% of angiosperms are dioecious. In dioecious higher plants, flowers borne on male and female individuals are, respectively deficient in functional gynoecium and androecium. Dioecy is inherited via three sex chromosome systems: XX/XY, XX/X0 and WZ/ZZ, such that XX or WZ is female and XY, X0 or ZZ are males. The XX/XY system generates the rarer XX/X0 and WZ/ZZ systems. An autosome pair begets XY chromosomes. A recessive loss-of-androecium mutation (ana) creates X chromosome and a dominant gynoecium-suppressing (GYS) mutation creates Y chromosome. The ana/ANA and gys/GYS loci are in the sex-determining region (SDR) of the XY pair. Accumulation of inversions, deleterious mutations and repeat elements, especially transposons, in the SDR of Y suppresses recombination between X and Y in SDR, making Y labile and increasingly degenerate and heteromorphic from X. Continued recombination between X and Y in their pseudoautosomal region located at the ends of chromosomal arms allows survival of the degenerated Y and of the species. Dioecy is presumably a component of the evolutionary cycle for the origin of new species. Inbred hermaphrodite species assume dioecy. Later they suffer degenerate-Y-led population regression. Cross-hybridization between such extinguishing species and heterologous species, followed by genome duplication of segregants from hybrids, give rise to new species.

Evolutionary origins and dynamics of octoploid strawberry subgenomes revealed by dense targeted capture linkage maps

Whole-genome duplications are radical evolutionary events that have driven speciation and adaptation in many taxa. Higher-order polyploids have complex histories often including interspecific hybridization and dynamic genomic changes. This chromosomal reshuffling is poorly understood for most polyploid species, despite their evolutionary and agricultural importance, due to the challenge of distinguishing homologous sequences from each other. Here, we use dense linkage maps generated with targeted sequence capture to improve the diploid strawberry (Fragaria vesca) reference genome and to disentangle the subgenomes of the wild octoploid progenitors of cultivated strawberry, Fragaria virginiana and Fragaria chiloensis. Our novel approach, POLiMAPS (Phylogenetics Of Linkage-Map-Anchored Polyploid Subgenomes), leverages sequence reads to associate informative interhomeolog phylogenetic markers with linkage groups and reference genome positions. In contrast to a widely accepted model, we find that one of the four subgenomes originates with the diploid cytoplasm donor F. vesca, one with the diploid Fragaria iinumae, and two with an unknown ancestor close to F. iinumae. Extensive unidirectional introgression has converted F. iinumae-like subgenomes to be more F. vesca-like, but never the reverse, due either to homoploid hybridization in the F. iinumae-like diploid ancestors or else strong selection spreading F. vesca-like sequence among subgenomes through homeologous exchange. In addition, divergence between homeologous chromosomes has been substantially augmented by interchromosomal rearrangements. Our phylogenetic approach reveals novel aspects of the complicated web of genetic exchanges that occur during polyploid evolution and suggests a path forward for unraveling other agriculturally and ecologically important polyploid genomes.

Ontogenetic stage-specific quantitative trait loci contribute to divergence in developmental trajectories of sexually dimorphic fins between medaka populations

Sexual dimorphism can evolve when males and females differ in phenotypic optima. Genetic constraints can, however, limit the evolution of sexual dimorphism. One possible constraint is derived from alleles expressed in both sexes. Because males and females share most of their genome, shared alleles with different fitness effects between sexes are faced with intralocus sexual conflict. Another potential constraint is derived from genetic correlations between developmental stages. Sexually dimorphic traits are often favoured at adult stages, but selected against as juvenile, so developmental decoupling of traits between ontogenetic stages may be necessary for the evolution of sexual dimorphism in adults. Resolving intralocus conflicts between sexes and ages is therefore a key to the evolution of age-specific expression of sexual dimorphism. We investigated the genetic architecture of divergence in the ontogeny of sexual dimorphism between two populations of the Japanese medaka (Oryzias latipes) that differ in the magnitude of dimorphism in anal and dorsal fin length. Quantitative trait loci (QTL) mapping revealed that few QTL had consistent effects throughout ontogenetic stages and the majority of QTL change the sizes and directions of effects on fin growth rates during ontogeny. We also found that most QTL were sex-specific, suggesting that intralocus sexual conflict is almost resolved. Our results indicate that sex- and age-specific QTL enable the populations to achieve optimal developmental trajectories of sexually dimorphic traits in response to complex natural and sexual selection.

Fragaria: a genus with deep historical roots and ripe for evolutionary and ecological insights

The cultivated strawberry, Fragaria ×ananassa, is one of the youngest domesticated plants. Its 18th century origin via hybridization in Europe between the North American F. virginiana and the South American F. chiloensis was documented by the botanist Antoine Nicolas Duchesne. His 1766 "Natural History of Strawberries" is an extraordinary work that integrates fundamental discoveries on the biology, ecology, and phylogeny of Fragaria with applied information on cultivation and ethnobotanical uses, serving as an inspiration for current research in the genus. Fragaria species exhibit the full range of sexual systems in the gynodioecy pathway from hermaphroditism to dioecy (and back again), as well as variation in self-compatibility, and evidence of sex chromosomes with female heterogamety. The genus is also characterized by interspecific hybridization and polyploidy, with a natural range of ploidy levels from diploids to decaploids. This biological diversity, combined with the availability of genomic resources and the ease of growing and experimenting with the plants, makes Fragaria a very attractive system for ecological and evolutionary genomics. The goal of this review is to introduce Fragaria as a model genus and to provide a roadmap for future integrative research. These research directions will deepen our understanding of the ecological and evolutionary context that shaped the ancestors of the cultivated strawberry, not only providing information that can be applied to efforts to shape the future of this important fruit crop but also our understanding of key transitions in plant evolution.

An angiosperm-wide analysis of the gynodioecy-dioecy pathway

BACKGROUND AND AIMS

About 6 % of an estimated total of 240 000 species of angiosperms are dioecious. The main precursors of this sexual system are thought to be monoecy and gynodioecy. A previous angiosperm-wide study revealed that many dioecious species have evolved through the monoecy pathway; some case studies and a large body of theoretical research also provide evidence in support of the gynodioecy pathway. If plants have evolved through the gynodioecy pathway, gynodioecious and dioecious species should co-occur in the same genera. However, to date, no large-scale analysis has been conducted to determine the prevalence of the gynodioecy pathway in angiosperms. In this study, this gap in knowledge was addressed by performing an angiosperm-wide survey in order to test for co-occurrence as evidence of the gynodioecy pathway.

METHODS

Data from different sources were compiled to obtain (to our knowledge) the largest dataset on gynodioecy available, with 275 genera that include at least one gynodioecious species. This dataset was combined with a dioecy dataset from the literature, and a study was made of how often dioecious and gynodioecious species could be found in the same genera using a contingency table framework.

KEY RESULTS

It was found that, overall, angiosperm genera with both gynodioecious and dioecious species occur more frequently than expected, in agreement with the gynodioecy pathway. Importantly, this trend holds when studying different classes separately (or sub-classes, orders and families), suggesting that the gynodioecy pathway is not restricted to a few taxa but may instead be widespread in angiosperms.

CONCLUSIONS

This work complements that previously carried out on the monoecy pathway and suggests that gynodioecy is also a common pathway in angiosperms. The results also identify angiosperm families where some (or all) dioecious species may have evolved from gynodioecious precursors. These families could be the targets of future small-scale studies on transitions to dioecy taking phylogeny explicitly into account.

The scope and strength of sex-specific selection in genome evolution

Males and females share the vast majority of their genomes and yet are often subject to different, even conflicting, selection. Genomic and transcriptomic developments have made it possible to assess sex-specific selection at the molecular level, and it is clear that sex-specific selection shapes the evolutionary properties of several genomic characteristics, including transcription, post-transcriptional regulation, imprinting, genome structure and gene sequence. Sex-specific selection is strongly influenced by mating system, which also causes neutral evolutionary changes that affect different regions of the genome in different ways. Here, we synthesize theoretical and molecular work in order to provide a cohesive view of the role of sex-specific selection and mating system in genome evolution. We also highlight the need for a combined approach, incorporating both genomic data and experimental phenotypic studies, in order to understand precisely how sex-specific selection drives evolutionary change across the genome.

Targeted sequence capture provides insight into genome structure and genetics of male sterility in a gynodioecious diploid strawberry, Fragaria vesca ssp. bracteata (Rosaceae)

Gynodioecy is a sexual system wherein females coexist with hermaphrodites. It is of interest not only because male-sterile plants are advantageous in plant breeding but also because it can be a crucial step in the evolutionary transition to entirely separate sexes (dioecy) from a hermaphroditic ancestor. The gynodioecious diploid wild strawberry, Fragaria vesca ssp. bracteata (Rosaceae), is a member of a clade with both dioecious and cultivated species, making it an ideal model in which to study the genetics of male sterility. To create a genetic map of F. v. ssp. bracteata, we identified informative polymorphisms from genomic sequencing (3−5x coverage) of two outbred plants from the same population. Using targeted enrichment, we sequenced 200 bp surrounding each of 6575 polymorphisms in 48 F1 offspring, yielding genotypes at 98% of targeted sites with mean coverage >100x, plus more than 600-kb high-coverage nontargeted sequence. With the resulting linkage map of 7802 stringently filtered markers (5417 targeted), we assessed recombination rates and genomic incongruities. Consistent with past work in strawberries, male sterility is dominant, segregates 1:1, and maps to a single location in the female. Further mapping an additional 55 offspring places male sterility in a gene-dense, 338-kb region of chromosome 4. The region is not syntenic with the sex-determining regions in the closely related octoploids, F. chiloensis and F. virginiana, suggesting either independent origins or translocation. The 57 genes in this region do not include protein families known to control male sterility and thus suggest alternate mechanisms for the suppression of male function.

Sex-determining chromosomes and sexual dimorphism: insights from genetic mapping of sex expression in a natural hybrid Fragaria × ananassa subsp. cuneifolia

We studied the natural hybrid (Fragaria × ananassa subsp. cuneifolia) between two sexually dimorphic octoploid strawberry species (Fragaria virginiana and Fragaria chiloensis) to gain insight into the dynamics of sex chromosomes and the genesis of sexual dimorphism. Male sterility is dominant in both the parental species and thus will be inherited maternally, but the chromosome that houses the sex-determining region differs. Thus, we asked whether (1) the cytotypic composition of hybrid populations represents one or both maternal species, (2) the sex-determining chromosome of the hybrid reflects the location of male sterility within the maternal donor species and (3) crosses from the hybrid species show less sexual dimorphism than the parental species. We found that F. × ananassa subsp. cuneifolia populations consisted of both parental cytotypes but one predominated within each population. Genetic linkage mapping of two crosses showed dominance of male sterility similar to the parental species, however, the map location of male sterility reflected the maternal donor in one cross, but not the other. Moreover, female function mapped to a single region in the first cross, but to two regions in the second cross. Aside from components of female function (fruit set and seed set), other traits that have been found to be significantly sexually dimorphic in the pure species were either not dimorphic or were dimorphic in the opposite direction to the parental species. These results suggest that hybrids experience some disruption of dimorphism in secondary sexual traits, as well as novel location and number of quantitative trait locus (QTL) affecting sex function.

Control of flower size

Flowers exhibit amazing morphological diversity in many traits, including their size. In addition to interspecific flower size differences, many species maintain significant variation in flower size within and among populations. Flower size variation can contribute to reproductive isolation of species and thus has clear evolutionary consequences. In this review we integrate information on flower size variation from both evolutionary and developmental biology perspectives. We examine the role of flower size in the context of mating system evolution. In addition, we describe what is currently known about the genetic basis of flower size based on quantitative trait locus (QTL) mapping in several different plant species and molecular genetic studies in model plants, primarily Arabidopsis thaliana. Work in Arabidopsis suggests that many independent pathways regulate floral organ growth via effects on cell proliferation and/or cell expansion.

Divergence in Defence against Herbivores between Males and Females of Dioecious Plant Species

Defensive traits may evolve differently between sexes in dioecious plant species. Our current understanding of this process hinges on a partial view of the evolution of resistance traits that may result in male-biased herbivory in dioecious populations. Here, we present a critical summary of the current state of the knowledge of herbivory in dioecious species and propose alternative evolutionary scenarios that have been neglected. These scenarios consider the potential evolutionary and functional determinants of sexual dimorphism in patterns of resource allocation to reproduction, growth, and defence. We review the evidence upon which two previous reviews of sex-biased herbivory have concluded that male-biased herbivory is a rule for dioecious species, and we caution readers about a series of shortcomings of many of these studies. Lastly, we propose a minimal standard protocol that should be followed in any studies that intend to elucidate the (co)evolution of interactions between dioecious plants and their herbivores.

Genetic dissection of fruit quality traits in the octoploid cultivated strawberry highlights the role of homoeo-QTL in their control

Fruit quality traits are major breeding targets in the Rosaceae. Several of the major Rosaceae species are current or ancient polyploids. To dissect the inheritance of fruit quality traits in polyploid fleshy fruit species, we used a cultivated strawberry segregating population comprising a 213 full-sibling F1 progeny from a cross between the variety 'Capitola' and the genotype 'CF1116'. We previously developed the most comprehensive strawberry linkage map, which displays seven homoeology groups (HG), including each four homoeology linkage groups (Genetics 179:2045-2060, 2008). The map was used to identify quantitative trait loci (QTL) for 19 fruit traits related to fruit development, texture, colour, anthocyanin, sugar and organic acid contents. Analyses were carried out over two or three successive years on field-grown plants. QTL were detected for all the analysed traits. Because strawberry is an octopolyploid species, QTL controlling a given trait and located at orthologous positions on different homoeologous linkage groups within one HG are considered as homoeo-QTL. We found that, for various traits, about one-fourth of QTL were putative homoeo-QTL and were localised on two linkage groups. Several homoeo-QTL could be detected the same year, suggesting that several copies of the gene underlying the QTL are functional. The detection of some other homoeo-QTL was year-dependent. Therefore, changes in allelic expression could take place in response to environmental changes. We believe that, in strawberry as in other polyploid fruit species, the mechanisms unravelled in the present study may play a crucial role in the variations of fruit quality.

About PAR: the distinct evolutionary dynamics of the pseudoautosomal region

Sex chromosomes differ from other chromosomes in the striking divergence they often show in size, structure, and gene content. Not only do they possess genes controlling sex determination that are restricted to either the X or Y (or Z or W) chromosomes, but in many taxa they also include recombining regions. In these 'pseudoautosomal regions' (PARs), sequence homology is maintained by meiotic pairing and exchange in the heterogametic sex. PARs are unique genomic regions, exhibiting some features of autosomes, but they are also influenced by their partial sex linkage. Here we review the distribution and structure of PARs among animals and plants, the theoretical predictions concerning their evolutionary dynamics, the reasons for their persistence, and the diversity and content of genes that reside within them. It is now clear that the evolution of the PAR differs in important ways from that of genes in either the non-recombining regions of sex chromosomes or the autosomes.

Preservation of the Y transcriptome in a 10-million-year-old plant sex chromosome system

Classical genetic studies discovered loss of genes from the ancient sex chromosome systems of several animals (genetic degeneration), and complete genome sequencing confirms that the heterogametic sex is hemizygous for most sex-linked genes. Genetic degeneration is thought to result from the absence of recombination between the sex chromosome pair (reviewed by [1]) and is very rapid after sex chromosome-autosome fusions in Drosophila [2-4]. Plant sex chromosome systems allow study of the time course of degeneration, because they evolved from a state wholly without sex chromosomes (rather than after a large genome region fused to a preexisting sex chromosome), and, in several taxa, recombination stopped very recently. However, despite increasing genetic and physical mapping of plant nonrecombining sex-determining regions [5-8], it remains very difficult to discover sex-linked genes, and it is unclear whether Y-linked genes are losing full function. We therefore developed a high-throughput method using RNA-Seq to identify sex linkage in Silene latifolia. Recombination suppression between this plant's XY sex chromosome pair started only about 10 million years ago [9]. Our approach identifies several hundred new sex-linked genes, and we show that this young Y chromosome retains many genes, yet these already have slightly reduced gene expression and are accumulating changes likely to reduce protein functions.

Gynodioecy to dioecy: are we there yet?

BACKGROUND

The 'gynodioecy-dioecy pathway' is considered to be one of the most important evolutionary routes from hermaphroditism to separate sexes (dioecy). Despite a large accumulation of evidence for female seed fertility advantages in gynodioecious species (females and hermaphrodites coexist) in support of the first step in the gynodioecy-dioecy pathway, we still have very little evidence for the second step, i.e. the transition from gynodioecy to dioecy.

SCOPE

We review the literature to evaluate whether basic predictions by theory are supported. To establish whether females' seed fertility advantage and frequencies are sufficient to favour the invasion of males, we review these for species along the gynodioecy-dioecy pathway published in the last 5 years. We then review the empirical evidence for predictions deriving from the second step, i.e. hermaphrodites' male fertility increases with female frequency, selection favours greater male fertility in hermaphrodites in gynodioecious species, and, where males and hermaphrodites coexist with females (subdioecy), males have greater male fertility than hermaphrodites. We review how genetic control and certain ecological features (pollen limitation, selfing, plasticity in sex expression and antagonists) influence the trajectory of a population along the gynodioecy-dioecy pathway.

CONCLUSIONS

Females tend to have greater seed fertility advantages over hermaphrodites where the two coexist, and this advantage is positively correlated with female frequency across species, as predicted by theory. A limited number of studies in subdioecious species have demonstrated that males have an advantage over hermaphrodites, as also predicted by theory. However, less evidence exists for phenotypic selection to increase male traits of hermaphrodites or for increasing male function of hermaphrodites in populations with high female frequency. A few key case studies underline the importance of examining multiple components of male fertility and the roles of pollen limitation, selfing and plasticity, when evaluating advantages. We conclude that we do not yet have a full understanding of the transition from gynodioecy to dioecy.