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.     

Insights into the genetic architecture of sexual dimorphism from an interspecific cross between two diverging Silene (Caryophyllaceae) species

The evolution of sexual dimorphism in species with separate sexes is influenced by the resolution of sexual conflicts creating sex differences through genetic linkage or sex‐biased expression. Plants with different degrees of sexual dimorphism are thus ideal to study the genetic basis of sexual dimorphism. In this study we explore the genetic architecture of sexual dimorphism between Silene latifolia and Silene dioica. These species have chromosomal sex determination and differ in the extent of sexual dimorphism. To test whether QTL for sexually dimorphic traits have accumulated on the sex chromosomes and to quantify their contribution to species differences, we create a linkage map and performed QTL analysis of life history, flower and vegetative traits using an unidirectional interspecific F2 hybrid cross. We found support for an accumulation of QTL on the sex chromosomes and that sex differences explained a large proportion of the variance between species, suggesting that both natural and sexual selection contributed to species divergence. Sexually dimorphic traits that also differed between species displayed transgressive segregation. We observed a reversal in sexual dimorphism in the F2 population, where males tended to be larger than females, indicating that sexual dimorphism is constrained within populations but not in recombinant hybrids. This study contributes to the understanding of the genetic basis of sexual dimorphism and its evolution in Silene.     

Assembly and annotation of a draft genome sequence for Glycine latifolia,a perennial wild relative of soybean

Glycine latifolia (Benth.) Newell & Hymowitz (2n = 40), one of the 27 wild perennial relatives of soybean, possesses genetic diversity and agronomically favorable traits that are lacking in soybean. Here, we report the 939‐Mb draft genome assembly of G. latifolia (PI 559298) using exclusively linked‐reads sequenced from a single Chromium library. We organized scaffolds into 20 chromosome‐scale pseudomolecules utilizing two genetic maps and the Glycine max (L.) Merr. genome sequence. High copy numbers of putative 91‐bp centromere‐specific tandem repeats were observed in consecutive blocks within predicted pericentromeric regions on several pseudomolecules. No 92‐bp putative centromeric repeats, which are abundant in G. max, were detected in G. latifolia or Glycine tomentella. Annotation of the assembled genome and subsequent filtering yielded a high confidence gene set of 54 475 protein‐coding loci. In comparative analysis with five legume species, genes related to defense responses were significantly overrepresented in Glycine‐specific orthologous gene families. A total of 304 putative nucleotide‐binding site (NBS)‐leucine‐rich‐repeat (LRR) genes were identified in this genome assembly. Different from other legume species, we observed a scarcity of TIR‐NBS‐LRR genes in G. latifolia. The G. latifolia genome was also predicted to contain genes encoding 367 LRR‐receptor‐like kinases, a family of proteins involved in basal defense responses and responses to abiotic stress. The genome sequence and annotation of G. latifolia provides a valuable source of alternative alleles and novel genes to facilitate soybean improvement. This study also highlights the efficacy and cost‐effectiveness of the application of Chromium linked‐reads in diploid plant genome de novo assembly.     

MK17, a specific marker closely linked to the gynoecium suppression region on the Y chromosome in Silene latifolia

The aim of this work was to isolate new DNA markers linked to the Silene latifolia Y chromosome. To do this we created a chromosome-specific plasmid library after DOP-PCR amplification of laser-microdissected Y-chromosomes. The library screening led to the isolation of several clones yielding mostly to exclusive male specific hybridization signals. Subsequent PCR confirmed the Y-unique linkage for one of the sequences. This DNA sequence called MK17 has no homology to any known DNA sequence and it is not expressed. Based on PCR and Southern analyses, MK17 is present only in dioecious species of the Elisanthe section of the genus Silene (S. latifolia, S. dioica, and S. diclinis) and it is absent in related gynodioecious and hermaphroditic species. The mapping analysis using a panel of deletion mutants showed that MK17 is closely linked to the region controlling suppression of gynoecium development. Hence MK17 represents a valuable marker to isolate genes controlling the gynoecium development suppression on the Y chromosome of S. latifolia.     

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.     

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.     

Evolution of putative barrier loci at an intermediate stage of speciation with gene flow in campions (Silene)

Understanding the origin of new species is a central goal in evolutionary biology. Diverging lineages often evolve highly heterogeneous patterns of genetic differentiation; however, the underlying mechanisms are not well understood. We investigated evolutionary processes governing genetic differentiation between the hybridizing campions Silene dioica (L.) Clairv. and S. latifolia Poiret. Demographic modelling indicated that the two species diverged with gene flow. The best‐supported scenario with heterogeneity in both migration rate and effective population size suggested that a small proportion of the loci evolved without gene flow. Differentiation (F _ST) and sequence divergence (d _XY) were correlated and both tended to peak in the middle of most linkage groups, consistent with reduced gene flow at highly differentiated loci. Highly differentiated loci further exhibited signatures of selection. In between‐species population pairs, isolation by distance was stronger for genomic regions with low between‐species differentiation than for highly differentiated regions that may contain barrier loci. Moreover, differentiation landscapes within and between species were only weakly correlated, suggesting that linked selection due to shared recombination and gene density landscapes is not the dominant determinant of genetic differentiation in these lineages. Instead, our results suggest that divergent selection shaped the genomic landscape of differentiation between the two Silene species, consistent with predictions for speciation in the face of gene flow.     

THE EVOLUTIONARY DYNAMICS OF SEXUALLY ANTAGONISTIC MUTATIONS IN PSEUDOAUTOSOMAL REGIONS OF SEX CHROMOSOMES

Sex chromosomes can evolve gene contents that differ from the rest of the genome, as well as larger sex differences in gene expression compared with autosomes. This probably occurs because fully sex‐linked beneficial mutations substitute at different rates from autosomal ones, especially when fitness effects are sexually antagonistic (SA). The evolutionary properties of genes located in the recombining pseudoautosomal region (PAR) of a sex chromosome have not previously been modeled in detail. Such PAR genes differ from classical sex‐linked genes by having two alleles at a locus in both sexes; in contrast to autosomal genes, however, variants can become associated with gender. The evolutionary fates of PAR genes may therefore differ from those of either autosomal or fully sex‐linked genes. Here, we model their evolutionary dynamics by deriving expressions for the selective advantages of PAR gene mutations under different conditions. We show that, unless selection is very strong, the probability of invasion of a population by an SA mutation is usually similar to that of an autosomal mutation, unless there is close linkage to the sex‐determining region. Most PAR genes should thus evolve similarly to autosomal rather than sex‐linked genes, unless recombination is very rare in the PAR.     

Consequences of genetic linkage for the maintenance of sexually antagonistic polymorphism in hermaphrodites

When selection differs between males and females, pleiotropic effects among genes expressed by both sexes can result in sexually antagonistic selection (SA), where beneficial alleles for one sex are deleterious for the other. For hermaphrodites, alleles with opposing fitness effects through each sex function represent analogous genetic constraints on fitness. Recent theory based on single‐locus models predicts that the maintenance of SA genetic variation should be greatly reduced in partially selfing populations. However, selfing also reduces the effective rate of recombination, which should facilitate selection on linked allelic combinations and expand opportunities for balancing selection in a multilocus context. Here, I develop a two‐locus model of SA selection for simultaneous hermaphrodites, and explore the joint influence of linkage, self‐fertilization, and dominance on the maintainance of SA polymorphism. I find that the effective reduction in recombination caused by selfing significantly expands the parameter space where SA polymorphism can be maintained relative to single‐locus models. In particular, linkage facilitates the invasion of male‐beneficial alleles, partially compensating for the “female‐bias” in the net direction of selection created by selfing. I discuss the implications of accounting for linkage among SA loci for the maintenance of SA genetic variation and mixed mating systems in hermaphrodites.     

EVOLUTION OF SEX DETERMINATION SYSTEMS WITH HETEROGAMETIC MALES AND FEMALES IN SILENE

The plant genus Silene has become a model for evolutionary studies of sex chromosomes and sex‐determining mechanisms. A recent study performed in Silene colpophylla showed that dioecy and the sex chromosomes in this species evolved independently from those in Silene latifolia, the most widely studied dioecious Silene species. The results of this study show that the sex‐determining system in Silene otites, a species related to S. colpophylla, is based on female heterogamety, a sex determination system that is unique among the Silene species studied to date. Our phylogenetic data support the placing of S. otites and S. colpophylla in the subsection Otites and the analysis of ancestral states suggests that the most recent common ancestor of S. otites and S. colpophylla was most probably dioecious. These observations imply that a switch from XX/XY sex determination to a ZZ/ZW system (or vice versa) occurred in the subsection Otites. This is the first report of two different types of heterogamety within one plant genus of this mostly nondioecious plant family.     

The potential for sexually antagonistic polymorphism in different genome regions

Sex differences in the fitness effects of alleles at a single locus (intralocus sexual antagonism, or SA) have several evolutionary consequences. Among the consequences of SA, polymorphisms at genes partially linked to the sex-determining region of the sex chromosome pair potentially drive the evolution of suppressed recombination between the sex chromosomes. Understanding the conditions under which SA polymorphism can exist at such pseudo-autosomal (or PAR) loci should increase understanding of the evolution of recombination between sex chromosome pairs, and can help predict when we may expect potentially empirically detectable allele frequency differences between the sexes. Models so far published have concluded that PAR genes can maintain SA polymorphisms over a wider range of selection coefficients than autosomal ones, but have used restrictive assumptions. We expand the modeling of SA alleles at a single locus with the full range of degrees of linkage to the male-specific region, to include strong or weak selection and the possibility of different dominance coefficients in the two sexes. We confirm the previous major conclusion that SA polymorphisms are generally maintained in a larger region of parameter space if the locus is in the PAR than if it is autosomal.     

Evolution and selection of Rhg1, a copy‐number variant nematode‐resistance locus

The soybean cyst nematode (SCN) resistance locus Rhg1 is a tandem repeat of a 31.2 kb unit of the soybean genome. Each 31.2‐kb unit contains four genes. One allele of Rhg1, Rhg1‐b, is responsible for protecting most US soybean production from SCN. Whole‐genome sequencing was performed, and PCR assays were developed to investigate allelic variation in sequence and copy number of the Rhg1 locus across a population of soybean germplasm accessions. Four distinct sequences of the 31.2‐kb repeat unit were identified, and some Rhg1 alleles carry up to three different types of repeat unit. The total number of copies of the repeat varies from 1 to 10 per haploid genome. Both copy number and sequence of the repeat correlate with the resistance phenotype, and the Rhg1 locus shows strong signatures of selection. Significant linkage disequilibrium in the genome outside the boundaries of the repeat allowed the Rhg1 genotype to be inferred using high‐density single nucleotide polymorphism genotyping of 15 996 accessions. Over 860 germplasm accessions were found likely to possess Rhg1 alleles. The regions surrounding the repeat show indications of non‐neutral evolution and high genetic variability in populations from different geographic locations, but without evidence of fixation of the resistant genotype. A compelling explanation of these results is that balancing selection is in operation at Rhg1.     

Chloroplast Genome of Native Silene latifolia subsp. alba from Fennoscandia Shows High Level of Differences from Invasive White Campion

Silene latifolia is an herbaceous plant with great invasive potential. Spread along trade routes from Europe to almost all continents, white campion became particularly widespread in North America. We sequenced the chloroplast genome of S. latifolia subsp. alba from a native range in southeast Fennoscandia. The chloroplast genome of native S. latifolia subsp. alba forms a 151,747-bp circle, has two inverted repeat regions (25,993 bp each), large single copy (82,708 bp), and small single copy (17,106 bp) regions. It contains 77 protein-coding genes, 30 tRNA genes, and four rRNA genes. SSRs and long DNA repeats were identified. Comparison of a newly sequenced plastome of S. latifolia subsp. alba with plastomes of invasive specimens of species from North America and Japan revealed a high level of single nucleotide polymorphisms (SNPs) among them. A total of 214 SNPs were found, among which 110 were identified in intergenic spacers, 74 in exons, and 30 in introns. Intraspecific shifts in inverted repeat boundaries were identified. Our research suggests that high polymorphic regions may be potential molecular markers for population studies and that high intraspecific genetic polymorphism may contribute to a species’ invasive success.

     

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.     

Recent selection for self‐compatibility in a population of Leavenworthia alabamica

The evolution of self‐compatibility (SC) is the first step in the evolutionary transition in plants from outcrossing enforced by self‐incompatibility (SI) to self‐fertilization. In the Brassicaceae, SI is controlled by alleles of two tightly linked genes at the S‐locus. Despite permitting inbreeding, mutations at the S‐locus leading to SC may be selected if they provide reproductive assurance and/or gain a transmission advantage in a population when SC plants self‐ and outcross. Positive selection can leave a genomic signature in the regions physically linked to the focus of selection when selection has occurred recently. From an SC population of Leavenworthia alabamica with a known nonfunctional mutation at the S‐locus, we collected sequence data from a ～690 Kb region surrounding the S‐locus, as well as from regions not linked to the S‐locus. To test for recent positive selection acting at the S‐locus, we examined polymorphism and the site‐frequency spectra. Using forward simulations, we demonstrate that recent selection of the strength expected for SC at a locus formerly under balancing selection can generate patterns similar to those seen in our empirical data.     

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.     

Differences in style length confer prezygotic isolation between two dioecious species of Silene in sympatry

One fundamental signature of reinforcement is elevated prezygotic reproductive isolation between related species in sympatry relative to allopatry. However, this alone is inadequate evidence for reinforcement, as traits conferring reproductive isolation can occur as a by‐product of other forces. We conducted crosses between Silene latifolia and S. diclinis, two closely related dioecious flowering plant species. Crosses with S. latifolia mothers from sympatry exhibited lower seed set than mothers from five allopatric populations when S. diclinis was the father. However, two other allopatric populations also exhibited low seed set. A significant interaction between style length and sire species revealed that seed set declined as style length increased when interspecific, but not intraspecific, fathers where used. Moreover, by varying the distance pollen tubes had to traverse, we found interspecific pollen placement close to the ovary resulted in seed set in both long‐ and short‐styled S. latifolia mothers. Our results reveal that the long styles of S. latifolia in sympatry with S. diclinis contribute to the prevention of hybrid formation. We argue that forces other than reinforcing selection are likely to be responsible for the differences in style length seen in sympatry.     

Clonal genetic structure and diversity in populations of an aquatic plant with combined vs. separate sexes

Clonality is often implicated in models of the evolution of dioecy, but few studies have explicitly compared clonal structure between plant sexual systems, or between the sexes in dioecious populations. Here, we exploit the occurrence of monoecy and dioecy in clonal Sagittaria latifola (Alismataceae) to evaluate two main hypotheses: (i) clone sizes are smaller in monoecious than dioecious populations, because of constraints imposed on clone size by costs associated with geitonogamy; (ii) in dioecious populations, male clones are larger and flower more often than female clones because of sex‐differential reproductive costs. Differences in clone size and flowering could result in discordance between ramet‐ and genet‐based sex ratios. We used spatially explicit sampling to address these hypotheses in 10 monoecious and 11 dioecious populations of S. latifolia at the northern range limit in Eastern North America. In contrast to our predictions, monoecious clones were significantly larger than dioecious clones, probably due to their higher rates of vegetative growth and corm production, and in dioecious populations, there was no difference in clone size between females and males; ramet‐ and genet‐based sex ratios were therefore highly correlated. Genotypic diversity declined with latitude for both sexual systems, but monoecious populations exhibited lower genotypic richness. Differences in life history between the sexual systems of S. latifolia appear to be the most important determinants of clonal structure and diversity.     

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.