Analysis of BAC clones containing homologous sequences on the end of the Xq arm and on chromosome 7 in the dioecious plant Silene latifolia

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.

Rapid degeneration of noncoding DNA regions surrounding SlAP3X/Y after recombination suppression in the dioecious plant Silene latifolia

Silene latifolia is a dioecious plant with heteromorphic XY sex chromosomes. Previous studies of sex chromosome–linked genes have suggested a gradual divergence between the X-linked and the Y-linked genes in proportion to the distance from the pseudoautosomal region. However, such a comparison has yet to be made for the noncoding regions. To better characterize the nonrecombining region of the X and Y chromosomes, we sequenced bacterial artificial chromosome clones containing the sex chromosome–linked paralogs SlAP3X and SlAP3Y, including 115 kb and 73 kb of sequences, respectively, flanking these genes. The synonymous nucleotide divergence between SlAP3X and SlAP3Y indicated that recombination stopped approximately 3.4 million years ago. Sequence homology analysis revealed the presence of six long terminal repeat retrotransposon-like elements. Using the nucleotide divergence calculated between left and right long terminal repeat sequences, insertion dates were estimated to be 0.083–1.6 million years ago, implying that all elements detected were inserted after recombination stopped. A reciprocal sequence homology search facilitated the identification of four homologous noncoding DNA regions between the X and Y chromosomes, spanning 6.7% and 10.6% of the X chromosome–derived and Y chromosome–derived sequences, respectively, investigated. Genomic Southern blotting and fluorescence in situ hybridization showed that the noncoding DNA flanking SlAP3X/Y has homology to many regions throughout the genome, regardless of whether they were homologous between the X and Y chromosomes. This finding suggests that most noncoding DNA regions rapidly lose their counterparts because of the introduction of transposable elements and indels (insertion–deletions) after recombination has stopped.

Dynamic gene order on the Silene latifolia 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.