Non-canonical bases differentially represented in the sex chromosomes of the dioecious plant Silene latifolia.

The oxidation of 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC), known as oxi-mCs, garners significant interest in plants as potential epigenetic marks. While research in mammals has established a role in cell reprogramming, carcinogenesis, and gene regulation, their functions in plants remain unclear. In rice, 5hmC has been associated with transposable elements (TEs) and heterochromatin. This study utilizes Silene latifolia, a dioecious plant with heteromorphic sex chromosomes and a genome with a large proportion of TEs, which provides a favourable environment for the study of oxi-mCs in individual sexes. Notably, we detected surprisingly high levels of oxi-mCs in S. latifolia comparable with mammals. Nuclei showed enrichment in heterochromatic regions, except for 5hmC whose signal was homogeneously distributed. Intriguingly, the same X chromosome in females displayed overall enrichment of 5hmC and 5fC compared with its counterpart. This fact is shared with 5mC, resembling dosage compensation. Co-localization showed higher correlation between 5mC and 5fC than with 5hmC, indicating no potential relationship between 5hmC and 5fC. Additionally, the promoter of several sex-linked genes and sex-biased TEs clustered in a clear sex-dependent way. Together, these findings unveil a hypothetical role for oxi-mCs in S. latifolia sex chromosome development, warranting further exploration.

Chemical genetics in Silene latifolia elucidate regulatory pathways involved in gynoecium development.

Dioecious plants possess diverse sex determination systems and unique mechanisms of reproductive organ development; however, little is known about how sex-linked genes shape the expression of regulatory cascades that lead to developmental differences between sexes. In Silene latifolia, a dioecious plant with stable dimorphism in floral traits, early experiments suggested that female-regulator genes act on the factors that determine the boundaries of the flower whorls. To identify these regulators, we sequenced the transcriptome of male flowers with fully developed gynoecia, induced by rapid demethylation in the parental generation. Eight candidates were found to have a positive role in gynoecium promotion, floral organ size, and whorl boundary, and affect the expression of class B MADS-box flower genes. To complement our transcriptome analysis, we closely examined the floral organs in their native state using field emission environmental scanning electron microscopy, and examined the differences between females and androhermaphrodites in their placenta and ovule organization. Our results reveal the regulatory pathways potentially involved in sex-specific flower development in the classical model of dioecy, S. latifolia. These pathways include previously hypothesized and unknown female-regulator genes that act on the factors that determine the flower boundaries, and a negative regulator of anther development, SUPERMAN-like (SlSUP).

Recent expansion of the non-recombining sex-linked region on Silene latifolia sex chromosomes.

Evolution of a non-recombining sex-specific region on the Y (or W) chromosome (NRY) is a key step in sex chromosome evolution, but how recombination suppression evolves is not well understood. Studies in many different organisms indicated that NRY evolution often involves several expansion steps. Why such NRY expansions occur remains unclear, although it is though that they are likely driven by sexually antagonistic selection. This paper describes a recent NRY expansion due to shift of the pseudoautosomal boundary on the sex chromosomes of a dioecious plant Silene latifolia. The shift resulted in inclusion of at least 16 pseudoautosomal genes into the NRY. This region is pseudoautosomal in closely related Silene dioica and Silene diclinis, indicating that the NRY expansion occurred in S. latifolia after it speciated from the other species ~120 thousand years ago. As S. latifolia and S. dioica actively hybridise across Europe, interspecific gene flow could blur the PAR boundary in these species. The pseudoautosomal genes have significantly elevated genetic diversity (π ~ 3% at synonymous sites), which is consistent with balancing selection maintaining diversity in this region. The recent shift of the PAR boundary in S. latifolia offers an opportunity to study the process of on-going NRY expansion.

Water availability drives population divergence and sex-specific responses in a dioecious plant.

PREMISE: Water availability is an important abiotic factor, resulting in differences between plant species growing in xeric and mesic habitats. Species with populations occurring in both habitat types allow examination of whether water availability has acted as a selective force at the intraspecific level. Investigating responses to water availability with a dioecious species allows determination of whether males and females, which often have different physiologies and life histories, respond differently. METHODS: An experiment varying water availability was performed under an outdoor rain-out shelter using plants from two mesic and two xeric populations of the dioecious plant Silene latifolia. Early growth rate, flowering propensity, flower size, and specific leaf area were measured. At the end of the season, the plants were harvested, aboveground and root biomass were measured, and the total number of flowers and fruit produced were counted. RESULTS: Compared to the two mesic populations, plants from the two xeric populations grew more slowly, were less likely to flower, took longer to flower, had thicker leaves, invested less in aboveground biomass and more in root biomass, produced fewer flowers and fruit, but were more likely to live. Many traits exhibited significant habitat type × treatment interactions. Compared to the xeric populations, males-but not females-from mesic populations had less root biomass and greatly reduced their flower production in response to low water availability. CONCLUSIONS: Mesic and xeric populations responded in ways congruent with water availability being a selective force for among-population divergence, especially for males.

The slowdown of Y chromosome expansion in dioecious Silene latifolia due to DNA loss and male-specific silencing of retrotransposons.

BACKGROUND: The rise and fall of the Y chromosome was demonstrated in animals but plants often possess the large evolutionarily young Y chromosome that is thought has expanded recently. Break-even points dividing expansion and shrinkage phase of plant Y chromosome evolution are still to be determined. To assess the size dynamics of the Y chromosome, we studied intraspecific genome size variation and genome composition of male and female individuals in a dioecious plant Silene latifolia, a well-established model for sex-chromosomes evolution. RESULTS: Our genome size data are the first to demonstrate that regardless of intraspecific genome size variation, Y chromosome has retained its size in S. latifolia. Bioinformatics study of genome composition showed that constancy of Y chromosome size was caused by Y chromosome DNA loss and the female-specific proliferation of recently active dominant retrotransposons. We show that several families of retrotransposons have contributed to genome size variation but not to Y chromosome size change. CONCLUSIONS: Our results suggest that the large Y chromosome of S. latifolia has slowed down or stopped its expansion. Female-specific proliferation of retrotransposons, enlarging the genome with exception of the Y chromosome, was probably caused by silencing of highly active retrotransposons in males and represents an adaptive mechanism to suppress degenerative processes in the haploid stage. Sex specific silencing of transposons might be widespread in plants but hidden in traditional hermaphroditic model plants.

DNA methylation and genetic degeneration of the Y chromosome in the dioecious plant Silene latifolia.

BACKGROUND: S. latifolia is a model organism for the study of sex chromosome evolution in plants. Its sex chromosomes include large regions in which recombination became gradually suppressed. The regions tend to expand over time resulting in the formation of evolutionary strata. Non-recombination and later accumulation of repetitive sequences is a putative cause of the size increase in the Y chromosome. Gene decay and accumulation of repetitive DNA are identified as key evolutionary events. Transposons in the X and Y chromosomes are distributed differently and there is a regulation of transposon insertion by DNA methylation of the target sequences, this points to an important role of DNA methylation during sex chromosome evolution in Silene latifolia. The aim of this study was to elucidate whether the reduced expression of the Y allele in S. latifolia is caused by genetic degeneration or if the cause is methylation triggered by transposons and repetitive sequences. RESULTS: Gene expression analysis in S. latifolia males has shown expression bias in both X and Y alleles. To determine whether these differences are caused by genetic degeneration or methylation spread by transposons and repetitive sequences, we selected several sex-linked genes with varying degrees of degeneration and from different evolutionary strata. Immunoprecipitation of methylated DNA (MeDIP) from promoter, exon and intron regions was used and validated through bisulfite sequencing. We found DNA methylation in males, and only in the promoter of genes of stratum I (older). The Y alleles in genes of stratum I were methylation enriched compared to X alleles. There was also abundant and high percentage methylation in the CHH context in most sequences, indicating de novo methylation through the RdDM pathway. CONCLUSIONS: We speculate that TE accumulation and not gene decay is the cause of DNA methylation in the S. latifolia Y sex chromosome with influence on the process of heterochromatinization.

Cell Death and Cell Cycle Arrest of Silene latifolia Stamens and Pistils After Microbotryum lychnidis-dioicae Infection.

Mechanisms of suppression of pistil primordia in male flowers and of stamen primordia in female flowers differ in diclinous plants. In this study, we investigated how cell death and cell cycle arrest are related to flower organ formation in Silene latifolia. Using in situ hybridization and a TUNEL (terminal deoxynucleotidyl transferase dUTP nick end labeling) assay, we detected both cell cycle arrest and cell death in suppressed stamens of female flowers and suppressed pistils of male flowers in S. latifolia. In female flowers infected with Microbotryum lychnidis-dioicae, developmental suppression of stamens is released, and cell cycle arrest and cell death do not occur. Smut spores are formed in S. latifolia anthers infected with M. lychnidis-dioicae, followed by cell death in the endothelium, middle layer, tapetal cells and pollen mother cells. Cell death is difficult to detect using a fluorescein isothiocyanate-labeled TUNEL assay due to strong autofluorescence in the anther. We therefore combined a TUNEL assay in an infrared region with transmission electron microscopy to detect cell death in anthers. We show that following infection by M. lychnidis-dioicae, a TUNEL signal was not detected in the endothelium, middle layer or pollen mother cells, and cell death with outflow of cell contents, including the nucleoplast, was observed in tapetal cells.

Differential adaptation drives ecological speciation in campions (Silene): evidence from a multi-site transplant experiment.

In order to investigate the role of differential adaptation for the evolution of reproductive barriers, we conducted a multi-site transplant experiment with the dioecious sister species Silene dioica and S. latifolia and their hybrids. Crosses within species as well as reciprocal first-generation (F1 ) and second-generation (F2 ) interspecific hybrids were transplanted into six sites, three within each species' habitat. Survival and flowering were recorded over 4 yr. At all transplant sites, the local species outperformed the foreign species, reciprocal F1 hybrids performed intermediately and F2 hybrids underperformed in comparison to F1 hybrids (hybrid breakdown). Females generally had slightly higher cumulative fitness than males in both within- and between-species crosses and we thus found little evidence for Haldane's rule acting on field performance. The strength of selection against F1 and F2 hybrids as well as hybrid breakdown increased with increasing strength of habitat adaptation (i.e. the relative fitness difference between the local and the foreign species) across sites. Our results suggest that differential habitat adaptation led to ecologically dependent post-zygotic reproductive barriers and drives divergence and speciation in this Silene system.

Sequence diversity patterns suggesting balancing selection in partially sex-linked genes of the plant Silene latifolia are not generated by demographic history or gene flow.

DNA sequence diversity in genes in the partially sex-linked pseudoautosomal region (PAR) of the sex chromosomes of the plant Silene latifolia is higher than expected from within-species diversity of other genes. This could be the footprint of sexually antagonistic (SA) alleles that are maintained by balancing selection in a PAR gene (or genes) and affect polymorphism in linked genome regions. SA selection is predicted to occur during sex chromosome evolution, but it is important to test whether the unexpectedly high sequence polymorphism could be explained without it, purely by the combined effects of partial linkage with the sex-determining region and the population's demographic history, including possible introgression from Silene dioica. To test this, we applied approximate Bayesian computation-based model choice to autosomal sequence diversity data, to find the most plausible scenario for the recent history of S. latifolia and then to estimate the posterior density of the most relevant parameters. We then used these densities to simulate variation to be expected at PAR genes. We conclude that an excess of variants at high frequencies at PAR genes should arise in S. latifolia populations only for genes with strong associations with fully sex-linked genes, which requires closer linkage with the fully sex-linked region than that estimated for the PAR genes where apparent deviations from neutrality were observed. These results support the need to invoke selection to explain the S. latifolia PAR gene diversity, and encourage further work to test the possibility of balancing selection due to sexual antagonism.

Degeneration of the nonrecombining regions in the mating-type chromosomes of the anther-smut fungi.

Dimorphic mating-type chromosomes in fungi are excellent models for understanding the genomic consequences of recombination suppression. Their suppressed recombination and reduced effective population size are expected to limit the efficacy of natural selection, leading to genomic degeneration. Our aim was to identify the sequences of the mating-type chromosomes (a1 and a2) of the anther-smut fungi and to investigate degeneration in their nonrecombining regions. We used the haploid a1 Microbotryum lychnidis-dioicae reference genome sequence. The a1 and a2 mating-type chromosomes were both isolated electrophoretically and sequenced. Integration with restriction-digest optical maps identified regions of recombination and nonrecombination in the mating-type chromosomes. Genome sequence data were also obtained for 12 other Microbotryum species. We found strong evidence of degeneration across the genus in the nonrecombining regions of the mating-type chromosomes, with significantly higher rates of nonsynonymous substitution (dN/dS) than in nonmating-type chromosomes or in recombining regions of the mating-type chromosomes. The nonrecombining regions of the mating-type chromosomes also showed high transposable element content, weak gene expression, and gene losses. The levels of degeneration did not differ between the a1 and a2 mating-type chromosomes, consistent with the lack of homogametic/heterogametic asymmetry between them, and contrasting with X/Y or Z/W sex chromosomes.

RAD mapping reveals an evolving, polymorphic and fuzzy boundary of a plant pseudoautosomal region.

How loss of genetic exchanges (recombination) evolves between sex chromosomes is a long-standing question. Suppressed recombination may evolve when a sexually antagonistic (SA) polymorphism occurs in a partially sex-linked 'pseudoautosomal' region (or 'PAR'), maintaining allele frequency differences between the two sexes, and creating selection for closer linkage with the fully sex-linked region of the Y chromosome in XY systems, or the W in ZW sex chromosome systems. Most evidence consistent with the SA polymorphism hypothesis is currently indirect, and more studies of the genetics and population genetics of PAR genes are clearly needed. The sex chromosomes of the plant Silene latifolia are suitable for such studies, as they evolved recently and the loss of recombination could still be ongoing. Here, we used RAD sequencing to genetically map sequences in this plant, which has a large genome (c. 3 gigabases) and no available whole-genome sequence. We mapped 83 genes on the sex chromosomes, and comparative mapping in the related species S. vulgaris supports previous evidence for additions to an ancestral PAR and identified at least 12 PAR genes. We describe evidence that recombination rates have been reduced in meiosis of both sexes, and differences in recombination between S. latifolia families suggest ongoing recombination suppression. Large allele frequency differences between the sexes were found at several loci closely linked to the PAR boundary, and genes in different regions of the PAR showed striking sequence diversity patterns that help illuminate the evolution of the PAR.

The population biology of fungal invasions.

Fungal invasions are increasingly recognized as a significant component of global changes, threatening ecosystem health and damaging food production. Invasive fungi also provide excellent models to evaluate the generality of results based on other eukaryotes. We first consider here the reasons why fungal invasions have long been overlooked: they tend to be inconspicuous, and inappropriate methods have been used for species recognition. We then review the information available on the patterns and mechanisms of fungal invasions. We examine the biological features underlying invasion success of certain fungal species. We review population structure analyses, revealing native source populations and strengths of bottlenecks. We highlight the documented ecological and evolutionary changes in invaded regions, including adaptation to temperature, increased virulence, hybridization, shifts to clonality and association with novel hosts. We discuss how the huge census size of most fungi allows adaptation even in bottlenecked, clonal invaders. We also present new analyses of the invasion of the anther-smut pathogen on white campion in North America, as a case study illustrating how an accurate knowledge of species limits and phylogeography of fungal populations can be used to decipher the origin of invasions. This case study shows that successful invasions can occur even when life history traits are particularly unfavourable to long-distance dispersal and even with a strong bottleneck. We conclude that fungal invasions are valuable models to contribute to our view of biological invasions, in particular by providing insights into the traits as well as ecological and evolutionary processes allowing successful introductions.

Size Does Matter: Staging of Silene latifolia Floral Buds for Transcriptome Studies.

Dioecious plants in the Caryophyllaceae family are susceptible to infection by members of the anthericolous smut fungi. In our studies of the Silene latifolia/Microbotryum lychnidis-dioicae pathosystem, we were interested in characterizing the plant-pathogen interaction at the molecular level before and during teliosporogenesis. This takes place during floral bud development, and we hoped to capture the interaction by Illumina Next-Gen RNA-Sequencing. Using previous literature that documented the stages of the floral buds for S. latifolia, we examined the floral buds from plants grown and infected under growth chamber conditions, using the disserting microscope to determine the stage of floral buds based on the morphology. We compiled the information and determined the size of floral buds that correspond to the desired stages of development for tissue collection, for the purpose of RNA-sequencing. This offers a practical approach for researchers who require a large number of floral buds/tissue categorized by stages of development, ascertaining whether infected/uninfected buds are at comparable stages of development and whether this also holds true for male vs. female buds. We also document our experience in infecting the plants and some of the unusual morphologies we observed after infection.

Possible mechanisms responsible for absence of a retrotransposon family on a plant Y chromosome.

Some transposable elements (TEs) show extraordinary variance in abundance along sex chromosomes but the mechanisms responsible for this variance are unknown. Here, we studied Ogre long terminal repeat (LTR) retrotransposons in Silene latifolia, a dioecious plant with evolutionarily young heteromorphic sex chromosomes. Ogre elements are ubiquitous in the S. latifolia genome but surprisingly absent on the Y chromosome. Bacterial artificial chromosome (BAC) library analysis and fluorescence in situ hybridization (FISH) were used to determine Ogre structure and chromosomal localization. Next generation sequencing (NGS) data were analysed to assess the transcription level and abundance of small RNAs. Methylation of Ogres was determined by bisulphite sequencing. Phylogenetic analysis was used to determine mobilization time and selection forces acting on Ogre elements. We characterized three Ogre families ubiquitous in the S. latifolia genome. One family is nearly absent on the Y chromosome despite all the families having similar structures and spreading mechanisms. We showed that Ogre retrotransposons evolved before sex chromosomes appeared but were mobilized after formation of the Y chromosome. Our data suggest that the absence of one Ogre family on the Y chromosome may be caused by 24-nucleotide (24-nt) small RNA-mediated silencing leading to female-specific spreading. Our findings highlight epigenetic silencing mechanisms as potentially crucial factors in sex-specific spreading of some TEs, but other possible mechanisms are also discussed.

Neo-sex chromosome inheritance across species in Silene hybrids.

Neo-sex chromosomes, which form through the major restructuring of ancestral sex chromosome systems, have evolved in various taxa. Such restructuring often consists of the fusion of an autosome to an existing sex chromosome, resulting in novel sex chromosome formations (e.g. X1X2Y or XY1Y2.). Comparative studies are often made between restructured sex chromosome systems of closely related species, and here we evaluate the consequences of variable sex chromosome systems to hybrids. If neo-sex chromosomes are improperly inherited across species, this could lead to aberrant development and reproductive isolation. In this study, we examine the fate of neo-sex chromosomes in hybrids of the flowering plants Silene diclinis and Silene latifolia. Silene diclinis has a neo-sex chromosome system (XY1Y2) that is thought to have evolved from an ancestral XY system that is still present in S. latifolia. These species do not hybridize naturally, and improper sex chromosome inheritance could contribute to reproductive isolation. We investigated whether this major restructuring of sex chromosomes prevents their proper inheritance in a variety of hybrid crosses, including some F2 - and later-generation hybrids, with sex chromosome-linked, species-specific, polymorphic markers and chromosome squashes. We discovered that despite the differences in sex chromosomes that exist between these two species, proper segregation had occurred in hybrids that made it to flowering, including later-generation hybrids, indicating that neo-sex chromosome formation alone does not result in complete reproductive isolation between these two species. Additionally, hybrids with aberrant sex expression (e.g. neuter, hermaphrodite) also inherited the restructured sex chromosomes properly, highlighting that issues with sexual development in hybrids can be caused by intrinsic genetic incompatibility rather than improper sex chromosome inheritance.

Seed shape and size of Silene latifolia, differences between sexes, and influence of the parental genome in hybrids with Silene dioica.

Introduction: Plants undergo various natural changes that dramatically modify their genomes. One is polyploidization and the second is hybridization. Both are regarded as key factors in plant evolution and result in phenotypic differences in different plant organs. In Silene, we can find both examples in nature, and this genus has a seed shape diversity that has long been recognized as a valuable source of information for infrageneric classification. Methods: Morphometric analysis is a statistical study of shape and size and their covariations with other variables. Traditionally, seed shape description was limited to an approximate comparison with geometric figures (rounded, globular, reniform, or heart-shaped). Seed shape quantification has been based on direct measurements, such as area, perimeter, length, and width, narrowing statistical analysis. We used seed images and processed them to obtain silhouettes. We performed geometric morphometric analyses, such as similarity to geometric models and elliptic Fourier analysis, to study the hybrid offspring of S. latifolia and S. dioica. Results: We generated synthetic tetraploids of Silene latifolia and performed controlled crosses between diploid S. latifolia and Silene dioica to analyze seed morphology. After imaging capture and post-processing, statistical analysis revealed differences in seed size, but not in shape, between S. latifolia diploids and tetraploids, as well as some differences in shape among the parentals and hybrids. A detailed inspection using fluorescence microscopy allowed for the identification of shape differences in the cells of the seed coat. In the case of hybrids, differences were found in circularity and solidity. Overal seed shape is maternally regulated for both species, whereas cell shape cannot be associated with any of the sexes. Discussion: Our results provide additional tools useful for the combination of morphology with genetics, ecology or taxonomy. Seed shape is a robust indicator that can be used as a complementary tool for the genetic and phylogenetic analyses of Silene hybrid populations.

Early research on anther-smut disease: A fuller view of science?

The anther-smut host-pathogen system has provided extensive insights into the evolutionary ecology of disease resistance, transmission modes, host shifts, pathogen specialization, and disease evolution in metapopulations. It also has led to unexpected insights into sex ratio distorters, sex chromosome evolution, and transposable elements in fungi. In addition, anther-smut disease played a major role in Linnaeus' germ theory and the correspondence on parasitic castration between Darwin and Becker, one of the first female botanists. Here, we explicitly highlight some of the realities in the process of science, using an unusual autobiographical approach to describe how we came to collaborate on this system in the 1980s. Using perspectives from our different career stages, we present a surprising narrative that could not be deduced from merely reading the published papers. While our work was grounded in previous ecological and evolutionary theory, it was the product as much of empirical failures and intellectual roadblocks, as the result of a progressive scientific method. Our experiences illustrate not only the "human dimension of science" but more importantly show that linear sequences of hypothesis testing do not necessarily lead to new study systems and new ideas. We suggest there is a need to re-evaluate the scientific method in ecology and evolution, especially where the challenge is to engage in a productive dialog between natural history and theory.

High temperatures reduce growth, infection, and transmission of a naturally occurring fungal plant pathogen.

Climate change is rapidly altering the distribution of suitable habitats for many species as well as their pathogenic microbes. For many pathogens, including vector-borne diseases of humans and agricultural pathogens, climate change is expected to increase transmission and lead to pathogen range expansions. However, if pathogens have a lower heat tolerance than their host, increased warming could generate so-called thermal refugia for hosts. Predicting the outcomes of warming on disease transmission requires detailed knowledge of the thermal tolerances of both the host and the pathogen. Such thermal tolerance studies are generally lacking for fungal pathogens of wild plant populations, despite the fact that plants form the base of all terrestrial communities. Here, we quantified three aspects of the thermal tolerance (growth, infection, and propagule production) of the naturally occurring fungal pathogen Microbotryum lychnidis-dioicae, which causes a sterilizing anther-smut disease on the herbaceous plant Silene latifolia. We also quantified two aspects of host thermal tolerance: seedling survival and flowering rate. We found that temperatures >30°C reduced the ability of anther-smut spores to germinate, grow, and conjugate in vitro. In addition, we found that high temperatures (30°C) during or shortly after the time of inoculation strongly reduced the likelihood of infection in seedlings. Finally, we found that high summer temperatures in the field temporarily cured infected plants, likely reducing transmission. Notably, high temperatures did not reduce survival or flowering of the host plants. Taken together, our results show that the fungus is considerably more sensitive to high temperatures than its host plant. A warming climate could therefore result in reduced disease spread or even local pathogen extirpation, leading to thermal refugia for the host.

Heterochiasmy and Sex Chromosome Evolution in Silene.

The evolution of a non-recombining sex-specific region is a key step in sex chromosome evolution. Suppression of recombination between the (proto-) X- and Y-chromosomes in male meiosis creates a non-recombining Y-linked region (NRY), while the X-chromosome continues to recombine in females. Lack of recombination in the NRY defines its main properties-genetic degeneration and accumulation of repetitive DNA, making X and Y chromosomes very different from each other. How and why recombination suppression on sex chromosomes evolves remains controversial. A strong difference in recombination rates between the sexes (heterochiasmy) can facilitate or even cause recombination suppression. In the extreme case-complete lack of recombination in the heterogametic sex (achiasmy)-the entire sex-specific chromosome is automatically non-recombining. In this study, I analyse sex-specific recombination rates in a dioecious plant Silene latifolia (Caryophyllaceae), which evolved separate sexes and sex chromosomes ~11 million years ago. I reconstruct high-density RNAseq-based genetic maps including over five thousand genic markers for the two sexes separately. The comparison of the male and female maps reveals only modest heterochiasmy across the genome, with the exception of the sex chromosomes, where recombination is suppressed in males. This indicates that heterochiasmy likely played only a minor, if any, role in NRY evolution in S. latifolia, as recombination suppression is specific to NRY rather than to the entire genome in males. Other mechanisms such as structural rearrangements and/or epigenetic modifications were likely involved, and comparative genome analysis and genetic mapping in multiple Silene species will help to shed light on the mechanism(s) of recombination suppression that led to the evolution of sex chromosomes.

The origin and evolution of sex chromosomes, revealed by sequencing of the Silene latifolia female genome.

White campion (Silene latifolia, Caryophyllaceae) was the first vascular plant where sex chromosomes were discovered. This species is a classic model for studies on plant sex chromosomes due to presence of large, clearly distinguishable X and Y chromosomes that originated de novo about 11 million years ago (mya), but lack of genomic resources for this relatively large genome (∼2.8 Gb) remains a significant hurdle. Here we report S. latifolia female genome assembly integrated with sex-specific genetic maps of this species, focusing on sex chromosomes and their evolution. The analysis reveals a highly heterogeneous recombination landscape with strong reduction in recombination rate in the central parts of all chromosomes. Recombination on the X chromosome in female meiosis primarily occurs at the very ends, and over 85% of the X chromosome length is located in a massive (∼330 Mb) gene-poor, rarely recombining pericentromeric region (Xpr). The results indicate that the non-recombining region on the Y chromosome (NRY) initially evolved in a relatively small (∼15 Mb), actively recombining region at the end of the q-arm, possibly as a result of inversion on the nascent X chromosome. The NRY expanded about 6 mya via linkage between the Xpr and the sex-determining region, which may have been caused by expanding pericentromeric recombination suppression on the X chromosome. These findings shed light on the origin of sex chromosomes in S. latifolia and yield genomic resources to assist ongoing and future investigations into sex chromosome evolution.