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.     

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.     

Asymmetrical conspecific seed-siring advantage between Silene latifolia and S. dioica

Background and Aims

Silene dioica and S. latifolia experience only limited introgression despite overlapping flowering phenologies, geographical distributions, and some pollinator sharing. Conspecific pollen precedence and other reproductive barriers operating between pollination and seed germination may limit hybridization. This study investigates whether barriers at this stage contribute to reproductive isolation between these species and, if so, which mechanisms are responsible.

Methods

Pollen-tube lengths for pollen of both species in styles of both species were compared. Additionally, both species were pollinated with majority S. latifolia and majority S. dioica pollen mixes; then seed set, seed germination rates and hybridity of the resulting seedlings were determined using species-specific molecular markers.

Key Results

The longest pollen tubes were significantly longer for conspecific than heterospecific pollen in both species, indicating conspecific pollen precedence. Seed set but not seed germination was lower for flowers pollinated with pure heterospecific versus pure conspecific pollen. Mixed-species pollinations resulted in disproportionately high representation of nonhybrid offspring for pollinations of S. latifolia but not S. dioica flowers.

Conclusions

The finding of conspecific pollen precedence for pollen-tube growth but not seed siring in S. dioica flowers may be explained by variation in pollen-tube growth rates, either at different locations in the style or between leading and trailing pollen tubes. Additionally, this study finds a barrier to hybridization operating between pollination and seed germination against S. dioica but not S. latifolia pollen. The results are consistent with the underlying cause of this barrier being attrition of S. dioica pollen tubes or reduced success of heterospecifically fertilized ovules, rather than time-variant mechanisms. Post-pollination, pre-germination barriers to hybridization thus play a partial role in limiting introgression between these 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.     

Substitution rates in a new Silene latifolia sex-linked gene, SlssX/Y

Dioecious white campion Silene latifolia has sex chromosomal sex determination, with homogametic (XX) females and heterogametic (XY) males. This species has become popular in studies of sex chromosome evolution. However, the lack of genes isolated from the X and Y chromosomes of this species is a major obstacle for such studies. Here, I report the isolation of a new sex-linked gene, Slss, with strong homology to spermidine synthase genes of other species. The new gene has homologous intact copies on the X and Y chromosomes (SlssX and SlssY, respectively). Synonymous divergence between the SlssX and SlssY genes is 4.7%, and nonsynonymous divergence is 1.4%. Isolation of a homologous gene from nondioecious S. vulgaris provided a root to the gene tree and allowed the estimation of the silent and replacement substitution rates along the SlssX and SlssY lineages. Interestingly, the Y-linked gene has higher synonymous and nonsynonymous substitution rates. The elevated synonymous rate in the SlssY gene, compared with SlssX, confirms our previous suggestion that the S. latifolia Y chromosome has a higher mutation rate, compared with the X chromosome. When differences in silent substitution rate are taken into account, the Y-linked gene still demonstrates significantly faster accumulation of nonsynonymous substitutions, which is consistent with the theoretical prediction of relaxed purifying selection in Y-linked genes, leading to the accumulation of nonsynonymous substitutions and genetic degeneration of the Y-linked genes.     

Genome size, quantitative genetics and the genomic basis for flower size evolution in Silene latifolia

BACKGROUND AND AIMS: The overall goal of this paper is to construct an overview of the genetic basis for flower size evolution in Silene latifolia. It aims to examine the relationship between the molecular bases for flower size and the underlying assumption of quantitative genetics theory that quantitative variation is ultimately due to the impact of a number of structural genes. SCOPE: Previous work is reviewed on the quantitative genetics and potential for response to selection on flower size, and the relationship between flower size and nuclear DNA content in S. latifolia. These earlier findings provide a framework within which to consider more recent analyses of a joint quantitative trait loci (QTL) analysis of flower size and DNA content in this species. KEY RESULTS: Flower size is a character that fits the classical quantitative genetics model of inheritance very nicely. However, an earlier finding that flower size is correlated with nuclear DNA content suggested that quantitative aspects of genome composition rather than allelic substitution at structural loci might play a major role in the evolution of flower size. The present results reported here show that QTL for flower size are correlated with QTL for DNA content, further corroborating an earlier result and providing additional support for the conclusion that localized variations in DNA content underlie evolutionary changes in flower size. CONCLUSIONS: The search image for QTL should be broadened to include overall aspects of genome regulation. As we prepare to enter the much-heralded post-genomic era, we also need to revisit our overall models of the relationship between genotype and phenotype to encompass aspects of genome structure and composition beyond structural genes.     

Qualitative and quantitative analyses of flower scent in Silene latifolia

The quantitative and qualitative variability in floral scent of 98 specimens of the dioecious species Silene latifolia belonging to 15 European and 19 North American populations was determined. Floral scent was collected from single flowers using dynamic headspace methods, and analysed by Micro-SPE and GC-MS methods. The flowers showed a nocturnal rhythm, and scent was emitted only at night. The amount of emitted volatiles varied greatly during the season, from 400 ng/flower/2 min in June to 50 ng/flower/2 min in August and September. The qualitative variability in the floral scent was high and different chemotypes, characterised by specific scent compounds, were found. Female and male flowers emitted the same type and amount of volatiles. The differences in floral scent composition between European and North American populations were small. Typical compounds were isoprenoids like lilac aldehyde isomers, or trans-beta-ocimene, and benzenoids like benzaldehyde, phenyl acetaldehyde, or veratrole. Some of these compounds are known to attract nocturnal Lepidoptera species. The high qualitative variability is discussed in relation to the pollination biology of S. latifolia, and the results are compared with other studies investigating intraspecific variability of flower scent.     

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.     

Expression of ASK1-like genes in arrested stamens of female Silene latifolia plants

Silene latifolia is a dioecious plant in which sex is determined by heteromorphic sex chromosomes. In female plants, stamen development is arrested before microspore mother cells are formed. In this study, we isolated four cDNAs (SlSKP1-1 to 4) encoding ASK1-like protein as expression markers to reveal when expression levels are reduced in arrested stamens of female flowers. Expression patterns of the SlSKP1 genes were analyzed by in-situ hybridization. We use the flower development classification of Grant et al. (in Plant J 6:471–480, 1994). SlSKP1 genes were highly expressed in primary parietal cells and primary sporogenous cells that develop into microspore mother cells in male flowers. Expression levels started to be reduced in the external stamens of the female flowers when stamen development was arrested at stage 7. Although microspore mother cells could not be developed in female flowers and SlSKP1 expression may be unnecessary in arrested stamens, SlSKP1 genes were still expressed in sporogenous cells of degenerated stamens at stage 8. Parietal cells stopped differentiating earlier than sporogenous cells in arrested stamens. These results suggest that not all types of cell are arrested simultaneously at a particular stage of stamen development during stamen suppression of S. latifolia.     

The relationship between pollen grain size and progeny gender in dioecious Silene latifolia (Caryophyllaceae)

Summary Experiments and observations conducted during the past 90 years have provided conflicting evidence concerning the existence of a size difference between pollen grains containing an X chromosome (female-determining) and those containing a Y chromosome (maledetermining) in dioecious Silene latifolia. Were such a size difference to exist, this might explain, at least in part, the observation that X-bearing pollen tubes reach the ovary more quickly, on average, than Y-bearing pollen tubes. We tested for such a size difference by separating pollen collected from single anthers into three size classes: small, large, and random. Fruit set (number of pollinated flowers that set fruit) and seed set (number of seeds per capsule) did not differ for these three pollination treatments. Progeny sex ratios resulting from these three pollen size classes also did not differ significantly. Thus, pollen grain size is not affected by which of the two sex chromosomes is present. Our experiment is the first direct test of this relationship. Based on our results, size differences should not be invoked to explain competitive differences in male- and female-determining microgametophytes.     

Sex-specific physiology and source-sink relations in the dioecious plant Silene latifolia

Differences in reproductive demands between the sexes of dioecious plants could cause divergence in physiology between the sexes. We found that the reproductive effort of female Silene latifolia plants increased to more than twice that of male plants or female plants that were prevented from setting fruit by lack of pollination after 4 weeks of flowering. Whole-plant source/sink ratios of pollinated females were significantly lower than those of males or unpollinated females because of investment in fruit. We hypothesized that these differences in source/sink ratio between the sexes and within females, depending on pollination, would lead to differences in leaf photosynthetic rates. Within females, we found that photosynthetic capacity was consistent with measurement of whole-plant source/sink ratio. Females that were setting fruit had 30% higher light-saturated photosynthetic rates by 28 days after flowering than females that were not setting fruit. Males, however, had consistently higher photosynthetic rates than females from 10 days after flowering onwards. Males also had approximately twice the dark respiration rates of fruiting females. We found that female reproductive structures are longer-lived and contribute more carbon to their own support than male reproductive structures. Despite the higher rates of leaf dark respiration and lower calyx photosynthetic rates, males fix more carbon than do females. We conclude that females have a sink-regulated mechanism of photosynthesis that allows them to respond to variations in fruit set. This mechanism is not, however, sufficient to explain why male S. latifolia plants have higher rates of photosynthesis, higher source/sink ratios, and lower reproductive allocation, but fail to grow larger than female plants.     

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.     

Stereoisomeric pattern of lilac aldehyde in Silene latifolia, a plant involved in a nursery pollination system

The monoterpene lilac aldehyde is found in floral scent of several plants species, among them Silene latifolia. This plant is involved in a nursery pollination system, because a noctuid moth, Hadena bicruris, is not only pollinator but also seed predator. Lilac aldehyde is the key floral scent compound of S. latifolia for attracting Hadena. This monoterpene has three stereogenic centers, and eight different isomers are possible. Here, we analysed the ratio of lilac aldehyde isomers from plants originating from 18 different populations of S. latifolia using enantioselective multidimensional GC-MS (enantio-MDGC-MS), and compared resulting variability with variability found in total scent emitted by specimen under study. Though variability in total emitted scent was high, ratio of lilac aldehyde isomers was a more conservative trait. There was no correlation between the ratio of lilac aldehyde isomers and the total emitted floral scent pattern. Both, ratio of stereoisomers and total emitted scent were independent from the geographic origin of the plants. In conclusion, the ratio of lilac aldehyde stereoisomers in S. latifolia is a reliable trait, and may used by the nursery pollinator H. bicruris for host-plant detection.     

Tackling the diversity of sex determination

A workshop on ‘The evolution of sex determination systems’ was held at a remote place in the Swiss Alps from 17 to 20 June 2009. It brought together theoreticians and empiricists, the latter ranging from molecular geneticists to evolutionary ecologists, all trying to understand key aspects of sex determination. The topics discussed included the evolutionary origins of sex determination, the diversity of sex determination mechanisms in different taxa, and the transition from genotypic to environmental sex determination and vice versa.     

Disentangling the effects of mating systems and mutation rates on cytoplamic diversity in gynodioecious Silene nutans and dioecious Silene otites

Many flowering plant species exhibit a variety of distinct sexual morphs, the two most common cases being the co-occurrence of females and males (dioecy) or the co-occurrence of hermaphrodites and females (gynodioecy). In this study, we compared DNA sequence variability of the three genomes (nuclear, mitochondrial and chloroplastic) of a gynodioecious species, Silene nutans, with that of a closely related dioecious species, Silene otites. In the light of theoretical models, we expect cytoplasmic diversity to differ between the two species due to the selective dynamics that acts on cytoplasmic genomes in gynodioecious species: under an epidemic scenario, the gynodioecious species is expected to exhibit lower cytoplasmic diversity than the dioecious species, while the opposite is expected in the case of balancing selection maintaining sterility cytoplasms in the gynodioecious species. We found no difference between the species for nuclear gene diversity, but, for the cytoplasmic loci, the gynodioecious S. nutans had more haplotypes, and higher nucleotide diversity, than the dioecious relative, S. otites, even though the latter has a relatively high rate of mitochondrial synonymous substitutions, and therefore presumably a higher mutation rate. Therefore, as the mitochondrial mutation rate cannot account for the higher cytoplasmic diversity found in S. nutans, our findings support the hypothesis that gynodioecy in S. nutans has been maintained by balancing selection rather than by epidemic-like dynamics.     

A rapid means of sex identification inSilene latifolia by use of flow cytometry

Sex identification in dioecious plants using nonflowering material would have broad applications in both basic and applied research. We present a method using flow cytometry for diagnosing the sex of the dioecious speciesSilene latifolia Poiret (Caryophyllaceae) by means of sexual differences in nuclear DNA content and base-pair composition. Males have a significantly larger genome, attributable to the known sex-chromosome heteromorphism. Males and females also differ in the AT/GC composition, attributable to differences in non-recombining portions of the sex chromosomes. The two measures enable assignment of individuals to sex with a combined error rate of 9%. These results forS. latifolia indicate useful directions for future research into sex diagnostics for other dioecious species.     

Ultrastructural analysis of the behavior of the dimorphic fungus Microbotryum violaceum in fungus-induced anthers of female Silene latifolia flowers

Summary. The development of male organs is induced in female flowers of the dioecious plant Silene latifolia by infection with the fungus Microbotryum violaceum. Stamens in a healthy female flower grow only to stage 6, whereas those in an infected female flower develop to the mature stage (stage 12), at which the stamens are filled with fungal teliospores instead of pollen grains. To investigate these host–parasite interactions, young floral buds and fungus-induced anthers of infected female flowers were examined by electron microscopy following fixation by a high-pressure freezing method. Using this approach, we found that parasitic hyphae of this fungus contain several extracellular vesicles and have a consistent appearance up to stage 8. At that stage, parasitic hyphae are observed adjacent to dying sporogenous cells in the infected female anther. At stage 9, an increased number of dead and dying sporogenous cells is observed, among which the sporogenous hyphae of the fungus develop and form initial teliospores. Several types of electron-dense material are present in proximity to some fungi at this stage. The initial teliospores contain two types of vacuoles, and the fungus cell wall contains abundant carbohydrate, as revealed by silver protein staining. The sporogenous cell is probably sensitive to infection by the fungus, resulting in disruption. In addition, the fungus accelerates cell death in the anther and utilizes constituents of the dead host cell to form the mature teliospore. Correspondence and reprints (present address): Molecular Membrane Biology Laboratory, RIKEN, 2-1, Hirosawa, Wako, Saitama 351-0198, Japan.     

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.     

Can the spread of agriculture in Europe be followed by tracing the spread of the weed Silene latifolia. A RAPD study

On the basis of gene frequency data of three flavone glycosylating genes, populations of the agricultural weed Silene latifolia (Caryophyllaceae) in Europe can be divided into two chemical races: an eastern and a western race. Morphological data also show a clear east-west division. When the two datasets are combined at least nine different geographical races can be distinguished using cluster analysis. Because these observations are hard to explain by selection, it has been proposed that these different races probably originated as a consequence of migration during the spread of agriculture over Europe in the past. To discriminate between selection and genetic drift many more selectively neutral easy-to-score characters are needed. In order to test whether random amplified polymorphic DNAs (RAPDs) might be suitable for this purpose, we performed a small-scale RAPD analysis on 16 geographical different populations. Using Jaccard's coefficient of similarity, we calculated genetic distances by pair-wise comparisons of both unique and shared amplification products, and a dendrogram was subsequently constructed using an unweighted pair-group method with arithmetical averages (UPGMA). On the basis of the dendrogram two clusters were discerned that clearly coincide with the aforementioned east-west division in populations. As there has been little or no artificial selection on this weed, its migration routes may be a good reflection of the different geographical routes agriculture has taken. We propose that a phylogenetic analysis of RAPD data of many more populations may provide additional information on the spread of agriculture over Europe.