A new plant sex-linked gene with high sequence diversity and possible introgression of the X copy

We describe patterns of DNA sequence diversity in a newly identified sex-linked gene, SlX9/SlY9, in Silene latifolia (Caryophyllaceae). The copies on both sex chromosomes seem to be functional, and each maps close to the respective X- and Y-linked copy of another sex-linked gene pair, SlCypX/SlCypY. The Y-linked copy has low diversity, similar to what has been found for several other Y-linked genes in S. latifolia, and consistent with the theoretical expectations of hitch-hiking processes occurring on a non-recombining chromosome. However, SlX9 has higher diversity than other genes on the S. latifolia X chromosome. We evaluate the hypothesis of introgression from the closely related species S. dioica as an explanation for the high sequence diversity observed.     

Indirect evidence from DNA sequence diversity for genetic degeneration of the Y-chromosome in dioecious species of the plant Silene: the SlY4/SlX4 and DD44-X/DD44-Y gene pairs

The action of natural selection is expected to reduce the effective population size of a nonrecombining chromosome, and this is thought to be the chief factor leading to genetic degeneration of Y-chromosomes, which cease recombining during their evolution from ordinary chromosomes. Low effective population size of Y chromosomes can be tested by studying DNA sequence diversity of Y-linked genes. In the dioecious plant, Silene latifolia, which has sex chromosomes, one comparison (SlX1 vs. SlY1) indeed finds lower Y diversity compared with the homologous X-linked gene, and one Y-linked gene with no X-linked homologue has lower species-wide diversity than a homologous autosomal copy (SlAp3Y vs. SlAp3A). To test whether this is a general pattern for Y-linked genes, we studied two further recently described X and Y homologous gene pairs in samples from several populations of S. latifolia and S. dioica. Diversity is reduced for both Y-linked genes, compared with their X-linked homologues. Our new data are analysed to show that the low Y effective size cannot be explained by different levels of gene flow for the X vs. the Y chromosomes, either between populations or between these closely related species. Thus, all four Y-linked genes that have now been studied in these plants (the two studied here, and two previously studied genes, have low diversity). This supports other evidence for an ongoing degeneration process in these species.     

Extreme population structure and high interspecific divergence of the Silene Y chromosome

Previous studies have demonstrated that the diversity of Y-linked genes is substantially lower than that of their X-linked homologs in the plant Silene latifolia. This difference has been attributed to selective sweeps, Muller's ratchet, and background selection, processes that are predicted to severely affect the evolution of the nonrecombining Y chromosome. We studied the DNA diversity of a noncoding region of the homologous genes DD44Y and DD44X, sampling S. latifolia populations from a wide geographical area and also including the closely related species S. dioica, S. diclinis, and S. heuffelii. On the Y chromosome of S. latifolia, we found substantial DNA diversity. Geographical population structure was far higher than on the X chromosome and differentiation between the species was also higher for the Y than for the X chromosome. Our findings indicate that the loss of genetic diversity on the Y chromosome in Silene occurs within local populations rather than within entire species. These results are compatible with background selection, Muller's ratchet, and local selective sweeps, but not with species-wide selective sweeps. The higher interspecific divergence of DD44Y, compared to DD44X, supports the hypothesis that Y chromosome differentiation between incipient species precedes reproductive isolation of the entire genome, forming an early stage in the process of speciation.     

DNA polymorphism in recombining and non-recombining mating-type-specific loci of the smut fungus Microbotryum

The population-genetic processes leading to the genetic degeneration of non-recombining regions have mainly been studied in animal and plant sex chromosomes. Here, we report population genetic analysis of the processes in the non-recombining mating-type-specific regions of the smut fungus Microbotryum violaceum. M. violaceum has A1 and A2 mating types, determined by mating-type-specific ‘sex chromosomes'' that contain 1–2 Mb long non-recombining regions. If genetic degeneration were occurring, then one would expect reduced DNA polymorphism in the non-recombining regions of this fungus. The analysis of DNA diversity among 19 M. violaceum strains, collected across Europe from Silene latifolia flowers, revealed that (i) DNA polymorphism is relatively low in all 20 studied loci (π∼0.15%), (ii) it is not significantly different between the two mating-type-specific chromosomes nor between the non-recombining and recombining regions, (iii) there is substantial population structure in M. violaceum populations, which resembles that of its host species, S. latifolia, and (iv) there is significant linkage disequilibrium, suggesting that widespread selfing in this species results in a reduction of the effective recombination rate across the genome. We hypothesise that selfing-related reduction of recombination across the M. violaceum genome negates the difference in the level of DNA polymorphism between the recombining and non-recombining regions, and may possibly lead to similar levels of genetic degeneration in the mating-type-specific regions of the non-recombining ‘sex chromosomes'' and elsewhere in the genome.     

DNA diversity in sex-linked and autosomal genes of the plant species Silene latifolia and Silene dioica

The relatively recent origin of sex chromosomes in the plant genus Silene provides an opportunity to study the early stages of sex chromosome evolution and, potentially, to test between the different population genetic processes likely to operate in nonrecombining chromosomes such as Y chromosomes. We previously reported much lower nucleotide polymorphism in a Y-linked gene (SlY1) of the plant Silene latifolia than in the homologous X-linked gene (SlX1). Here, we report a more extensive study of nucleotide diversity in these sex-linked genes, including a larger S. latifolia sample and a sample from the closely related species Silene dioica, and we also study the diversity of an autosomal gene, CCLS37.1. We demonstrate that nucleotide diversity in the Y-linked genes of both S. latifolia and S. dioica is very low compared with that of the X-linked gene. However, the autosomal gene also has low DNA polymorphism, which may be due to a selective sweep. We use a single individual of the related hermaphrodite species Silene conica, as an outgroup to show that the low SlY1 diversity is not due to a lower mutation rate than that for the X-linked gene. We also investigate several other possibilities for the low SlY1 diversity, including differential gene flow between the two species for Y-linked, X-linked, and autosomal genes. The frequency spectrum of nucleotide polymorphism on the Y chromosome deviates significantly from that expected under a selective-sweep model. However, we detect population subdivision in both S. latifolia and S. dioica, so it is not simple to test for selective sweeps. We also discuss the possibility that Y-linked diversity is reduced due to highly variable male reproductive success, and we conclude that this explanation is unlikely.     

DOES LOCAL ADAPTATION CAUSE HIGH POPULATION DIFFERENTIATION OF SILENE LATIFOLIA Y CHROMOSOMES?

Natural selection can reduce the effective population size of the nonrecombining Y chromosome, whereas local adaptation of Y-linked genes can increase the population divergence and overall intra-species polymorphism of Y-linked sequences. The plant Silene latifolia evolved a Y chromosome relatively recently, and most known X-linked genes have functional Y homologues, making the species interesting for comparisons of X- and Y-linked diversity and subdivision. Y-linked genes show higher population differentiation, compared to X-linked genes, and this might be maintained by local adaptation in Y-linked genes (or low sequence diversity). Here we attempt to test between these causes by investigating DNA polymorphism and population differentiation using a larger set of Y-linked and X-linked S. latifolia genes (than used previously), and show that net sequence divergence for Y-linked sequences (measured by D(a) , also known as δ) is low, and not consistently higher than X-linked genes. This does not support local adaptation, instead, the higher values of differentiation measures for the Y-linked genes probably result largely from reduced total variation on the Y chromosome, which in turn reflect deterministic processes lowering effective population sizes of evolving Y-chromosomes.     

Natural variation of the Y chromosome suppresses sex ratio distortion and modulates testis-specific gene expression in Drosophila simulans

X-linked sex-ratio distorters that disrupt spermatogenesis can cause a deficiency in functional Y-bearing sperm and a female-biased sex ratio. Y-linked modifiers that restore a normal sex ratio might be abundant and favored when a X-linked distorter is present. Here we investigated natural variation of Y-linked suppressors of sex-ratio in the Winters systems and the ability of these chromosomes to modulate gene expression in Drosophila simulans. Seventy-eight Y chromosomes of worldwide origin were assayed for their resistance to the X-linked sex-ratio distorter gene Dox. Y chromosome diversity caused males to sire ∼63% to ∼98% female progeny. Genome-wide gene expression analysis revealed hundreds of genes differentially expressed between isogenic males with sensitive (high sex ratio) and resistant (low sex ratio) Y chromosomes from the same population. Although the expression of about 75% of all testis-specific genes remained unchanged across Y chromosomes, a subset of post-meiotic genes was upregulated by resistant Y chromosomes. Conversely, a set of accessory gland-specific genes and mitochondrial genes were downregulated in males with resistant Y chromosomes. The D. simulans Y chromosome also modulated gene expression in XXY females in which the Y-linked protein-coding genes are not transcribed. The data suggest that the Y chromosome might exert its regulatory functions through epigenetic mechanisms that do not require the expression of protein-coding genes. The gene network that modulates sex ratio distortion by the Y chromosome is poorly understood, other than that it might include interactions with mitochondria and enriched for genes expressed in post-meiotic stages of spermatogenesis.     

Plant sex chromosomes: molecular structure and function

Recent molecular and genomic studies carried out in a number of model dioecious plant species, including Asparagus officinalis, Carica papaya, Silene latifolia, Rumex acetosa and Marchantia polymorpha, have shed light on the molecular structure of both homomorphic and heteromorphic sex chromosomes, and also on the gene functions they have maintained since their evolution from a pair of autosomes. The molecular structure of sex chromosomes in species from different plant families represents the evolutionary pathway followed by sex chromosomes during their evolution. The degree of Y chromosome degeneration that accompanies the suppression of recombination between the Xs and Ys differs among species. The primitive Ys of A. officinalis and C. papaya have only diverged from their homomorphic Xs in a short male-specific and non-recombining region (MSY), while the heteromorphic Ys of S. latifolia, R. acetosa and M. polymorpha have diverged from their respective Xs. As in the Y chromosomes of mammals and Drosophila, the accumulation of repetitive DNA, including both transposable elements and satellite DNA, has played an important role in the divergence and size enlargement of plant Ys, and consequently in reducing gene density. Nevertheless, the degeneration process in plants does not appear to have reached the Y-linked genes. Although a low gene density has been found in the sequenced Y chromosome of M. polymorpha, most of its genes are essential and are expressed in the vegetative and reproductive organs in both male and females. Similarly, most of the Y-linked genes that have been isolated and characterized up to now in S. latifolia are housekeeping genes that have X-linked homologues, and are therefore expressed in both males and females. Only one of them seems to be degenerate with respect to its homologous region in the X. Sequence analysis of larger regions in the homomorphic X and Y chromosomes of papaya and asparagus, and also in the heteromorphic sex chromosomes of S. latifolia and R. acetosa, will reveal the degenerative changes that the Y-linked gene functions have experienced during sex chromosome evolution.     

Accumulation of chloroplast DNA sequences on the Y chromosome of Silene latifolia

Silene latifolia is a model dioecious plant with heteromorphic sex chromosomes. The Y chromosome is the largest in this species. Theoretical models propose an accumulation of repetitive DNA sequences in non-recombining parts of the Y chromosome. In this study, we isolated a BAC7H5 clone preferentially hybridizing to the Y chromosome of S. latifolia. Sequence analysis revealed that this BAC7H5 contains part of the chloroplast genome, indicating that these chloroplast sequences have accumulated on the Y chromosome and also may contribute to its large size. We constructed Y chromosome- and X chromosome-specific libraries and screened them to find Y- and/or X-linked copies of chloroplast sequences. Sequence analysis revealed higher divergence of a non-genic region of the chloroplast sequences located on the Y chromosome while genic regions tested showed only very low (max 0.9%) divergence from their chloroplast homologues.     

Testing for the Footprint of Sexually Antagonistic Polymorphisms in the Pseudoautosomal Region of a Plant Sex Chromosome Pair

The existence of sexually antagonistic (SA) polymorphism is widely considered the most likely explanation for the evolution of suppressed recombination of sex chromosome pairs. This explanation is largely untested empirically, and no such polymorphisms have been identified, other than in fish, where no evidence directly implicates these genes in events causing loss of recombination. We tested for the presence of loci with SA polymorphism in the plant Silene latifolia, which is dioecious (with separate male and female individuals) and has a pair of highly heteromorphic sex chromosomes, with XY males. Suppressed recombination between much of the Y and X sex chromosomes evolved in several steps, and the results in Bergero et al. (2013) show that it is still ongoing in the recombining or pseudoautosomal, regions (PARs) of these chromosomes. We used molecular evolutionary approaches to test for the footprints of SA polymorphisms, based on sequence diversity levels in S. latifolia PAR genes identified by genetic mapping. Nucleotide diversity is high for at least four of six PAR genes identified, and our data suggest the existence of polymorphisms maintained by balancing selection in this genome region, since molecular evolutionary (HKA) tests exclude an elevated mutation rate, and other tests also suggest balancing selection. The presence of sexually antagonistic alleles at a locus or loci in the PAR is suggested by the very different X and Y chromosome allele frequencies for at least one PAR gene.     

Identifying new sex-linked genes through BAC sequencing in the dioecious plant Silene latifolia

Background

Silene latifolia represents one of the best-studied plant sex chromosome systems. A new approach using RNA-seq data has recently identified hundreds of new sex-linked genes in this species. However, this approach is expected to miss genes that are either not expressed or are expressed at low levels in the tissue(s) used for RNA-seq. Therefore other independent approaches are needed to discover such sex-linked genes.

Results

Here we used 10 well-characterized S. latifolia sex-linked genes and their homologs in Silene vulgaris, a species without sex chromosomes, to screen BAC libraries of both species. We isolated and sequenced 4 Mb of BAC clones of S. latifolia X and Y and S. vulgaris genomic regions, which yielded 59 new sex-linked genes (with S. vulgaris homologs for some of them). We assembled sequences that we believe represent the tip of the Xq arm. These sequences are clearly not pseudoautosomal, so we infer that the S. latifolia X has a single pseudoautosomal region (PAR) on the Xp arm. The estimated mean gene density in X BACs is 2.2 times lower than that in S. vulgaris BACs, agreeing with the genome size difference between these species. Gene density was estimated to be extremely low in the Y BAC clones. We compared our BAC-located genes with the sex-linked genes identified in previous RNA-seq studies, and found that about half of them (those with low expression in flower buds) were not identified as sex-linked in previous RNA-seq studies. We compiled a set of ~70 validated X/Y genes and X-hemizygous genes (without Y copies) from the literature, and used these genes to show that X-hemizygous genes have a higher probability of being undetected by the RNA-seq approach, compared with X/Y genes; we used this to estimate that about 30 % of our BAC-located genes must be X-hemizygous. The estimate is similar when we use BAC-located genes that have S. vulgaris homologs, which excludes genes that were gained by the X chromosome.

Conclusions

Our BAC sequencing identified 59 new sex-linked genes, and our analysis of these BAC-located genes, in combination with RNA-seq data suggests that gene losses from the S. latifolia Y chromosome could be as high as 30 %, higher than previous estimates of 10-20 %.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1698-7) contains supplementary material, which is available to authorized users.     

Sex-Ratio Meiotic Drive and Y-Linked Resistance in Drosophila affinis

Genetic elements that cheat Mendelian segregation by biasing transmission in their favor gain a significant fitness benefit. Several examples of sex-ratio meiotic drive, where one sex chromosome biases its own transmission at the cost of the opposite sex chromosome, exist in animals and plants. While the distorting sex chromosome gains a significant advantage by biasing sex ratio, the autosomes, and especially the opposite sex chromosome, experience strong selection to resist this transmission bias. In most well-studied sex-ratio meiotic drive systems, autosomal and/or Y-linked resistance has been identified. We specifically surveyed for Y-linked resistance to sex-ratio meiotic drive in Drosophila affinis by scoring the sex ratio of offspring sired by males with a driving X and one of several Y chromosomes. Two distinct types of resistance were identified: a restoration to 50/50 sex ratios and a complete reversal of sex ratio to all sons. We confirmed that fathers siring all sons lacked a Y chromosome, consistent with previously published work. Considerable variation in Y-chromosome morphology exists in D. affinis, but we showed that morphology does not appear to be associated with resistance to sex-ratio meiotic drive. We then used two X chromosomes (driving and standard) and three Y chromosomes (susceptible, resistant, and lacking) to examine fertility effects of all possible combinations. We find that both the driving X and resistant and lacking Y have significant fertility defects manifested in microscopic examination of testes and a 48-hr sperm depletion assay. Maintenance of variation in this sex-ratio meiotic drive system, including both the X-linked distorter and the Y-resistant effects, appear to be mediated by a complex interaction between fertility fitness and transmission dynamics.     

Evidence for degeneration of the Y chromosome in the dioecious plant Silene latifolia

The human Y--probably because of its nonrecombining nature--has lost 97% of its genes since X and Y chromosomes started to diverge [1, 2]. There are clear signs of degeneration in the Drosophila miranda neoY chromosome (an autosome fused to the Y chromosome), with neoY genes showing faster protein evolution [3-6], accumulation of unpreferred codons [6], more insertions of transposable elements [5, 7], and lower levels of expression [8] than neoX genes. In the many other taxa with sex chromosomes, Y degeneration has hardly been studied. In plants, many genes are expressed in pollen [9], and strong pollen selection may oppose the degeneration of plant Y chromosomes [10]. Silene latifolia is a dioecious plant with young heteromorphic sex chromosomes [11, 12]. Here we test whether the S. latifolia Y chromosome is undergoing genetic degeneration by analyzing seven sex-linked genes. S. latifolia Y-linked genes tend to evolve faster at the protein level than their X-linked homologs, and they have lower expression levels. Several Y gene introns have increased in length, with evidence for transposable-element accumulation. We detect signs of degeneration in most of the Y-linked gene sequences analyzed, similar to those of animal Y-linked and neo-Y chromosome genes.     

The role of chromosomal rearrangements in the evolution of <Emphasis Type="Italic">Silene latifolia</Emphasis> sex chromosomes

Silene latifolia is a model plant for studies of the early steps of sex chromosome evolution. In comparison to mammalian sex chromosomes that evolved 300 mya, sex chromosomes of S. latifolia appeared approximately 20 mya. Here, we combine results from physical mapping of sex-linked genes using polymerase chain reaction on microdissected arms of the S. latifolia X chromosome, and fluorescence in situ hybridization analysis of a new cytogenetic marker, Silene tandem repeat accumulated on the Y chromosome. The data are interpreted in the light of current genetic linkage maps of the X chromosome and a physical map of the Y chromosome. Our results identify the position of the centromere relative to the mapped genes on the X chromosome. We suggest that the evolution of the S. latifolia Y chromosome has been accompanied by at least one paracentric and one pericentric inversion. These results indicate that large chromosomal rearrangements have played an important role in Y chromosome evolution in S. latifolia and that chromosomal rearrangements are an integral part of sex chromosome evolution.     

Dynamic gene order on the <Emphasis Type="Italic">Silene latifolia</Emphasis> Y chromosome

Dioecious Silene latifolia evolved heteromorphic sex chromosomes within the last ten million years, making it a species of choice for studies of the early stages of sex chromosome evolution in plants. About a dozen genes have been isolated from its sex chromosomes and basic genetic and deletion maps exist for the X and Y chromosomes. However, discrepancies between Y chromosome maps led to the proposal that individual Y chromosomes may differ in gene order. Here, we use an alternative approach, with fluorescence in situ hybridization (FISH), to locate individual genes on S. latifolia sex chromosomes. We demonstrate that gene order on the Y chromosome differs between plants from two populations. We suggest that dynamic gene order may be a general property of Y chromosomes in species with XY systems, in view of recent work demonstrating that the gene order on the Y chromosomes of humans and chimpanzees are dramatically different.     

Surprising Differences in the Variability of Y Chromosomes in African and Cosmopolitan Populations of Drosophila melanogaster

The nonrecombining Drosophila melanogaster Y chromosome is heterochromatic and has few genes. Despite these limitations, there remains ample opportunity for natural selection to act on the genes that are vital for male fertility and on Y factors that modulate gene expression elsewhere in the genome. Y chromosomes of many organisms have low levels of nucleotide variability, but a formal survey of D. melanogaster Y chromosome variation had yet to be performed. Here we surveyed Y-linked variation in six populations of D. melanogaster spread across the globe. We find surprisingly low levels of variability in African relative to Cosmopolitan (i.e., non-African) populations. While the low levels of Cosmopolitan Y chromosome polymorphism can be explained by the demographic histories of these populations, the staggeringly low polymorphism of African Y chromosomes cannot be explained by demographic history. An explanation that is entirely consistent with the data is that the Y chromosomes of Zimbabwe and Uganda populations have experienced recent selective sweeps. Interestingly, the Zimbabwe and Uganda Y chromosomes differ: in Zimbabwe, a European Y chromosome appears to have swept through the population.     

Rapid de novo evolution of X chromosome dosage compensation in Silene latifolia, a plant with young sex chromosomes

Silene latifolia is a dioecious plant with heteromorphic sex chromosomes that have originated only ～10 million years ago and is a promising model organism to study sex chromosome evolution in plants. Previous work suggests that S. latifolia XY chromosomes have gradually stopped recombining and the Y chromosome is undergoing degeneration as in animal sex chromosomes. However, this work has been limited by the paucity of sex-linked genes available. Here, we used 35 Gb of RNA-seq data from multiple males (XY) and females (XX) of an S. latifolia inbred line to detect sex-linked SNPs and identified more than 1,700 sex-linked contigs (with X-linked and Y-linked alleles). Analyses using known sex-linked and autosomal genes, together with simulations indicate that these newly identified sex-linked contigs are reliable. Using read numbers, we then estimated expression levels of X-linked and Y-linked alleles in males and found an overall trend of reduced expression of Y-linked alleles, consistent with a widespread ongoing degeneration of the S. latifolia Y chromosome. By comparing expression intensities of X-linked alleles in males and females, we found that X-linked allele expression increases as Y-linked allele expression decreases in males, which makes expression of sex-linked contigs similar in both sexes. This phenomenon is known as dosage compensation and has so far only been observed in evolutionary old animal sex chromosome systems. Our results suggest that dosage compensation has evolved in plants and that it can quickly evolve de novo after the origin of sex chromosomes.     

Molecular and evolutionary analysis of a plant Y chromosome.

Plants have evolved a great diversity of sex determination systems. Among these, the XY system, also found in mammals, is one of the most exciting since it gives the opportunity to compare the evolution of sex chromosomes in two different kingdoms. Whereas genetic and molecular mechanisms controlling sex determination in drosophila and mammals, have been well studied, very little is known about such processes in plants. White campion (Silene latifolia) is an example of plant with X and Y chromosomes. What is the origin of the X and Y chromosomes? How did they evolve from a pair of autosomes? In our laboratory, we have isolated the first active genes located on a plant Y chromosome. We are using them as markers to trace the origin and evolution of sex chromosomes in the Silene genus.     

X-linked genes evolve higher codon bias in Drosophila and Caenorhabditis

Comparing patterns of molecular evolution between autosomes and sex chromosomes (such as X and W chromosomes) can provide insight into the forces underlying genome evolution. Here we investigate patterns of codon bias evolution on the X chromosome and autosomes in Drosophila and Caenorhabditis. We demonstrate that X-linked genes have significantly higher codon bias compared to autosomal genes in both Drosophila and Caenorhabditis. Furthermore, genes that become X-linked evolve higher codon bias gradually, over tens of millions of years. We provide several lines of evidence that this elevation in codon bias is due exclusively to their chromosomal location and not to any other property of X-linked genes. We present two possible explanations for these observations. One possibility is that natural selection is more efficient on the X chromosome due to effective haploidy of the X chromosomes in males and persistently low effective numbers of reproducing males compared to that of females. Alternatively, X-linked genes might experience stronger natural selection for higher codon bias as a result of maladaptive reduction of their dosage engendered by the loss of the Y-linked homologs.