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.     

Nucleotide diversity in Silene latifolia autosomal and sex-linked genes

The plant Silene latifolia has separate sexes and sex chromosomes, and is of interest for studying the early stages of sex chromosome evolution, especially the evolution of non-recombining regions on the Y chromosome. Hitch-hiking processes associated with ongoing genetic degeneration of the non-recombining Y chromosome are predicted to reduce Y-linked genes'' effective population sizes, and S. latifolia Y-linked genes indeed have lower diversity than X-linked ones. We tested whether this represents a true diversity reduction on the Y, versus the alternative possibility, elevated diversity at X-linked genes, by collecting new data on nucleotide diversity for autosomal genes, which had previously been little studied. We find clear evidence that Y-linked genes have reduced diversity. However, another alternative explanation for a low Y effective size is a high variance in male reproductive success. Autosomal genes should then also have lower diversity than expected, relative to the X, but this is not found in our loci. Taking into account the higher mutation rate of Y-linked genes, their low sequence diversity indicates a strong effect of within-population hitch-hiking on the Y chromosome.     

Comparison of the X and Y chromosome organization in Silene latifolia

Here we compare gene orders on the Silene latifolia sex chromosomes. On the basis of the deletion mapping results (11 markers and 23 independent Y chromosome deletion lines used), we conclude that a part of the Y chromosome (covering a region corresponding to at least 23.9 cM on the X chromosome) has been inverted. The gradient in silent-site divergence suggests that this inversion took place after the recombination arrest in this region. Because recombination arrest events followed by Y chromosome rearrangements also have been found in the human Y chromosome, this process seems to be a general evolutionary pathway.     

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.     

Duplicative transfer of a MADS box gene to a plant Y chromosome

Y chromosomes carry genes with functions in male reproduction and often have few other loci. Their evolution and the causes of genetic degeneration are of great interest. In addition to genetic degeneration, the acquisition of autosomal genes may be important in Y chromosome evolution. We here report that the dioecious plant Silene latifolia harbors a complete MADS box gene, SlAP3Y, duplicated onto the Y chromosome. This gene has no X-linked homologs but only an autosomal paralog, SlAP3A, and sequence divergence suggests that the duplication is a quite old event that occurred soon after the evolution of the sex chromosomes. Evolutionary sequence analyses using homologs of closely related species, including hermaphroditic Silene conica and dioecious Silene dioica and Silene diclinis, suggest that both SlAP3A and SlAP3Y genes encode functional proteins. Indeed, quantitative RT-PCR and in situ hybridization analyses showed that SlAP3A is expressed specifically in developing petals, but SlAP3Y is much more strongly expressed in developing stamens. The S. conica homolog, ScAP3A, is expressed in developing petals, suggesting subfunctionalization with evolution of male-specific functions, possibly due to evolutionary change in regulatory elements. Our results suggest that the acquisition of autosomal genes is an important event in the evolution of plant Y chromosomes.     

A sex-chromosome mutation in Silene latifolia

Silene latifolia is dioecious, yet rare hermaphrodites have been found, and such natural mutants can provide valuable insight into genetic mechanisms. Here, we describe a hermaphrodite-inducing mutation that is almost certainly localized to the gynoecium-suppression region of the Y chromosome in S. latifolia. The mutant Y chromosome was passed through the megaspore, and the presence of two X chromosomes was not necessary for seed development in the parent. This result supports a lack of degeneration of the Y chromosome in S. latifolia, consistent with the relatively recent formation of the sex chromosomes in this species. When crossed to wild-type plants, hermaphrodites performed poorly as females, producing low seed numbers. When hermaphrodites were pollen donors, the sex ratio of offspring they produced through crosses was biased towards females. This suggests that hermaphroditic S. latifolia would fail to thrive and potentially explains the rarity of hermaphrodites in natural populations of S. latifolia. These results indicate that the Y chromosome in Silene latifolia remains very similar to the X, perhaps mostly differing in the primary sex determination regions.     

An accumulation of tandem DNA repeats on the Y chromosome in Silene latifolia during early stages of sex chromosome evolution

Sex chromosomes in mammals are about 300 million years old and typically have a highly degenerated Y chromosome. The sex chromosomes in the dioecious plant Silene latifolia in contrast, represent an early stage of evolution in which functional X–Y gene pairs are still frequent. In this study, we characterize a novel tandem repeat called TRAYC, which has accumulated on the Y chromosome in S. latifolia. Its presence demonstrates that processes of satellite accumulation are at work even in this early stage of sex chromosome evolution. The presence of TRAYC in other species of the Elisanthe section suggests that this repeat had spread after the sex chromosomes evolved but before speciation within this section. TRAYC possesses a palindromic character and a strong potential to form secondary structures, which could play a role in satellite evolution. TRAYC accumulation is most prominent near the centromere of the Y chromosome. We propose a role for the centromere as a starting point for the cessation of recombination between the X and Y 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.     

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.     

Structure and evolution of <Emphasis Type="Italic">Apetala3</Emphasis>, a sex-linked gene in <Emphasis Type="Italic">Silene latifolia</Emphasis>

Background  

The evolution of sex chromosomes is often accompanied by gene or chromosome rearrangements. Recently, the gene AP3 was characterized in the dioecious plant species Silene latifolia. It was suggested that this gene had been transferred from an autosome to the Y chromosome.     

DNA probes for the Y-chromosome ofSilene latifolia,a dioecious angiosperm

Summary In order to obtain markers for the Y chromosome ofSilene latifolia, we pooled equal weights of leaf tissue from 18 female siblings into one sample and repeated the process with 18 male siblings. Pooling was intended to provide a common genetic background for each sample, leaving the absence or presence of the Y chromosome as the primary difference between the two samples. DNA was extracted from each sample and subjected to polymerase chain reaction (PCR) amplification with arbitrary 10 bp primers. Four of 60 primers used gave an amplification with the male DNA not found among those from the female DNA. Each of these was subsequently shown to provide a reliable marker for the Y chromosome.     

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.     

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.     

Semi-automatic laser beam microdissection of the Y chromosome and analysis of Y chromosome DNA in a dioecious plant, Silene latifolia

Silene latifolia has heteromorphic sex chromosomes, the X and Y chromosomes. The Y chromosome, which is thought to carry the male determining gene, was isolated by UV laser microdissection and amplified by degenerate oligonucleotide-primed PCR. In situ chromosome suppression of the amplified Y chromosome DNA in the presence of female genomic DNA as a competitor showed that the microdissected Y chromosome DNA did not specifically hybridize to the Y chromosome, but hybridized to all chromosomes. This result suggests that the Y chromosome does not contain Y chromosome-enriched repetitive sequences. A repetitive sequence in the microdissected Y chromosome, RMY1, was isolated while screening repetitive sequences in the amplified Y chromosome. Part of the nucleotide sequence shared a similarity to that of X-43.1, which was isolated from microdissected X chromosomes. Since fluorescence in situ hybridization analysis with RMY1 demonstrated that RMY1 was localized at the ends of the chromosome, RMY1 may be a subtelomeric repetitive sequence. Regarding the sex chromosomes, RMY1 was detected at both ends of the X chromosome and at one end near the pseudoautosomal region of the Y chromosome. The different localization of RMY1 on the sex chromosomes provides a clue to the problem of how the sex chromosomes arose from autosomes.     

The sex chromosomes of Silene latifolia revisited and revised

Classical studies have established that, during meiosis, the X and Y chromosomes of the model dioecious plant Silene latifolia pair over a region at the ends of their q arms. We used fluorescence in situ hybridization of two molecular markers to demonstrate that this widely accepted model is incorrect. From these data we conclude that the homologous arm of the X chromosome is the p arm and that of the Y chromosome is the q arm. The establishment of the proper orientation of the pseudoautosomal region is essential for mapping and evolutionary studies.     

Genetic and functional analysis of DD44, a sex-linked gene from the dioecious plant Silene latifolia,provides clues to early events in sex chromosome evolution

Silene latifolia is a dioecious plant with heteromorphic sex chromosomes. The sex chromosomes of S. latifolia provide an opportunity to study the early events in sex chromosome evolution because of their relatively recent emergence. In this article, we present the genetic and physical mapping, expression analysis, and molecular evolutionary analysis of a sex-linked gene from S. latifolia, DD44 (Differential Display 44). DD44 is homologous to the oligomycin sensitivity-conferring protein, an essential component of the mitochondrial ATP synthase, and is ubiquitously expressed in both sexes. We have been able to genetically map DD44 to a region of the Y chromosome that is genetically linked to the carpel-suppressing locus. Although we have physically mapped DD44 to the distal end of the long arm of the X chromosome using fluorescence in situ hybridization (FISH), DD44 maps to the opposite arm of the Y chromosome as determined by our genetic map. These data suggest that chromosomal rearrangements have occurred on the Y chromosome, which may have contributed to the genetic isolation of the Y chromosome. We discuss the implications of these results with respect to the structural and functional evolution of the S. latifolia Y chromosome.     

C-banded karyotypes of two Silene species with heteromorphic sex chromosomes.

Mitotic metaphase chromosomes of Silene latifolia (white campion) and Silene dioica (red campion) were studied and no substantial differences between the conventional karyotypes of these two species were detected. The classification of chromosomes into three distinct groups proposed for S. latifolia by Ciupercescu and colleagues was considered and discussed. Additionally, a new small satellite on the shorter arm of homobrachial chromosome 5 was found. Giemsa C-banded chromosomes of the two analysed species show many fixed and polymorphic heterochromatic bands, mainly distally and centromerically located. Our C-banding studies provided an opportunity to better characterize the sex chromosomes and some autosome types, and to detect differences between the two Silene karyotypes. It was shown that S. latifolia possesses a larger amount of polymorphic heterochromatin, especially of the centromeric type. The two Silene sex chromosomes are easily distinguishable not only by length or DNA amount differences but also by their Giemsa C-banding patterns. All Y chromosomes invariably show only one distally located band, and no other fixed or polymorphic bands on this chromosome were observed in either species. The X chromosomes possess two terminally located fixed bands, and some S. latifolia X chromosomes also have an extra-centric segment of variable length. The heterochromatin amount and distribution revealed by our Giemsa C-banding studies provide a clue to the problem of sex chromosome and karyotype evolution in these two closely related dioecious Silene species.     

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.