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.     

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.     

Heterogeneous Histories of Recombination Suppression on Stickleback Sex Chromosomes

How consistent are the evolutionary trajectories of sex chromosomes shortly after they form? Insights into the evolution of recombination, differentiation, and degeneration can be provided by comparing closely related species with homologous sex chromosomes. The sex chromosomes of the threespine stickleback (Gasterosteus aculeatus) and its sister species, the Japan Sea stickleback (G. nipponicus), have been well characterized. Little is known, however, about the sex chromosomes of their congener, the blackspotted stickleback (G. wheatlandi). We used pedigrees to obtain experimentally phased whole genome sequences from blackspotted stickleback X and Y chromosomes. Using multispecies gene trees and analysis of shared duplications, we demonstrate that Chromosome 19 is the ancestral sex chromosome and that its oldest stratum evolved in the common ancestor of the genus. After the blackspotted lineage diverged, its sex chromosomes experienced independent and more extensive recombination suppression, greater X–Y differentiation, and a much higher rate of Y degeneration than the other two species. These patterns may result from a smaller effective population size in the blackspotted stickleback. A recent fusion between the ancestral blackspotted stickleback Y chromosome and Chromosome 12, which produced a neo-X and neo-Y, may have been favored by the very small size of the recombining region on the ancestral sex chromosome. We identify six strata on the ancestral and neo-sex chromosomes where recombination between the X and Y ceased at different times. These results confirm that sex chromosomes can evolve large differences within and between species over short evolutionary timescales.     

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.     

Extensive Divergence Between Mating-Type Chromosomes of the Anther-Smut Fungus

Genomic regions that determine mating compatibility are subject to distinct evolutionary forces that can lead to a cessation of meiotic recombination and the accumulation of structural changes between members of the homologous chromosome pair. The relatively recent discovery of dimorphic mating-type chromosomes in fungi can aid the understanding of sex chromosome evolution that is common to dioecious plants and animals. For the anther-smut fungus, Microbotryum lychnidis-dioicae (= M. violaceum isolated from Silene latifolia), the extent of recombination cessation on the dimorphic mating-type chromosomes has been conflictingly reported. Comparison of restriction digest optical maps for the two mating-type chromosomes shows that divergence extends over 90% of the chromosome lengths, flanked at either end by two pseudoautosomal regions. Evidence to support the expansion of recombination cessation in stages from the mating-type locus toward the pseudoautosomal regions was not found, but evidence of such expansion could be obscured by ongoing processes that affect genome structure. This study encourages the comparison of forces that may drive large-scale recombination suppression in fungi and other eukaryotes characterized by dimorphic chromosome pairs associated with sexual life cycles.     

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.     

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

1. Download : Download high-res image (180KB)
2. Download : Download full-size image

     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Avian Neo-Sex Chromosomes Reveal Dynamics of Recombination Suppression and W Degeneration

How the avian sex chromosomes first evolved from autosomes remains elusive as 100 million years (My) of divergence and degeneration obscure their evolutionary history. The Sylvioidea group of songbirds is interesting for understanding avian sex chromosome evolution because a chromosome fusion event ∼24 Ma formed “neo-sex chromosomes” consisting of an added (new) and an ancestral (old) part. Here, we report the complete female genome (ZW) of one Sylvioidea species, the great reed warbler (Acrocephalus arundinaceus). Our long-read assembly shows that the added region has been translocated to both Z and W, and whereas the added-Z has retained its gene order the added-W part has been heavily rearranged. Phylogenetic analyses show that recombination between the homologous added-Z and -W regions continued after the fusion event, and that recombination suppression across this region took several million years to be completed. Moreover, recombination suppression was initiated across multiple positions over the added-Z, which is not consistent with a simple linear progression starting from the fusion point. As expected following recombination suppression, the added-W show signs of degeneration including repeat accumulation and gene loss. Finally, we present evidence for nonrandom maintenance of slowly evolving and dosage-sensitive genes on both ancestral- and added-W, a process causing correlated evolution among orthologous genes across broad taxonomic groups, regardless of sex linkage.